Ocean acidification and warming are challenging marine organisms and ecosystems around the world. The synergetic effects of these two climate change stressors on jellyfish remain still understudied. Here, we examine the independent and combined effects of these two environmental variables on polyp population dynamics of the Mediterranean jellyfish Cotylorhiza tuberculata. An experiment was conducted to examine asexual reproduction by budding and strobilation considering current and ca. 2100 winter (Trial 1, 36 days) and summer (Trial 2, 36 days) conditions under the RCP8.5 (IPCC 2013). In Trial 1, a temperature of 18°C and two pH levels (current: 7.9 and, reduced: 7.7) were tested. Trial 2 considered two temperature levels 24°C and 30°C, under current and reduced acidification conditions (8.0 and 7.7, respectively). Ephyrae size and statolith formation of released ephyrae from polyps exposed to summer temperatures under both acidification treatment was also analyzed. Zooxanthellae density inside the polyps throughout the experiment was measured. C. tuberculata polyps could cope with the conditions mimicked in all experimental treatments and no significant effect of pH, temperature, or the combination of both variables on the abundance of polyps was observed. At 18°C, strobilation was reduced under high PCO2 conditions. Under summer treatments (24°C and 30°C), percentage strobilation was very low and several released ephyrae suffered malformations and reduced size, as a consequence of reduced pH and elevated temperatures, separately. The number of statoliths was not affected by pH or temperature, however, bigger statoliths were formed at elevated temperatures (30°C). Finally, zooxanthellae density was not affected by experimental conditions, even if, the duration of the experiment significantly affected symbiont concentration. Our results show that even though polyps of C. tuberculata would thrive the future worst scenario predicted for the Mediterranean Sea, their capacity to undergo a proper strobilation and to produce healthy ephyrae will be more vulnerable to climate induced environmental conditions, thereby affecting medusae recruitment and, therefore, population dynamics of the species.

In order to allow full comparability with other ocean acidification data sets, the R package seacarb (Gattuso et al, 2021) was used to compute a complete and consistent set of carbonate system variables, as described by Nisumaa et al. (2010). In this dataset the original values were archived in addition with the recalculated parameters (see related PI). The date of carbonate chemistry calculation by seacarb is 2022-3-14.