Buffering and Amplifying Interactions among OAW (Ocean Acidification & Warming) and Nutrient Enrichment on Early Life-Stage Fucus vesiculosus L. (Phaeophyceae) and Their Carry Over Effects to Hypoxia Impact

Ocean acidification and warming (OAW) are occurring globally. Additionally, at a more local scale the spreading of hypoxic conditions is promoted by eutrophication and warming. In the semi-enclosed brackish Baltic Sea, occasional upwelling in late summer and autumn may expose even shallow-water communities including the macroalga Fucus vesiculosus to particularly acidified, nutrient-rich and oxygen-poor water bodies. During summer 2014 (July–September) sibling groups of early life-stage F. vesiculosus were exposed to OAW in the presence and absence of enhanced nutrient levels and, subsequently to a single upwelling event in a near-natural scenario which included all environmental fluctuations in the Kiel Fjord, southwestern Baltic Sea, Germany (54°27 ´N, 10°11 ´W). We strove to elucidate the single and combined impacts of these potential stressors, and how stress sensitivity varies among genetically different sibling groups. Enhanced by a circumstantial natural heat wave, warming and acidification increased mortalities and reduced growth in F. vesiculosus germlings. This impact, however, was mitigated by enhanced nutrient conditions. Survival under OAW conditions strongly varied among sibling groups hinting at a substantial adaptive potential of the natural Fucus populations in the Western Baltic. A three-day experimental upwelling caused severe mortality of Fucus germlings, which was substantially more severe in those sibling groups which previously had been exposed to OAW. Our results show that global (OAW), regional (nutrient enrichment) and local pressures (upwelling), both alone and co-occurring may have synergistic and antagonistic effects on survival and/or growth of Fucus germlings. This result emphasizes the need to consider combined stress effects.


Introduction
Global climate change will expose marine populations to increased eutrophication and upwelling events at a regional scale and to ocean acidification and warming (OAW) at a more global scale [1]. For the Baltic Sea, current models predict an increase in pCO 2 from currently almost 400 μatm to 1000 μatm, an increase in sea surface temperature by up to 5°C and considerable enhancements of eutrophication and hypoxia during the next 100 years [2]. The combined effects of these global change factors will likely affect Baltic ecosystems [3]. Anthropogenic activities, such as agriculture, lead to eutrophic conditions in the largest part of the Baltic Sea, including the southwestern part with Kiel Bight [4]. Eutrophication may further intensify due to increased precipitation and river runoff [5,6]. Hypoxic conditions (< 2 mg O 2 L -1 ) are predicted to increase in the Baltic Sea during this century as modelled by Meier and Andersson [7]. In Kiel Bay, seasonal oxygen depletion has been known for decades [8], but is likely to intensify further.
Macrophytes play a key role in ecosystem services by the retention of excessive nutrients [9,10], acting as 'nutrient filters' [11] and providing the benthic ecosystem with oxygen [12]. A combination of global, regional and local stressors endangers macroalgae worldwide [13,14]. Among these, the bladder wrack Fucus vesiculosus is a dominant perennial macroalga in the intertidal and shallow subtidal of the Western, Central and Eastern Baltic Sea [15]. As a foundation species, F. vesiculosus provides habitat and food for a large variety of invertebrate assemblages [16,17]. A drastic decline in F. vesiculosus during the last five decades was observed, e.g. by almost 95% in Kiel Bay, Western Baltic Sea [18,19]. The shoaling of F. vesiculosus was mainly assigned to the severe indirect effects of eutrophication [18,[20][21][22]. These are i.a. the increased turbidity [20], competition with ephemeral algae [23] and increased palatability [24]. Also, enhanced sedimentation reduces the attachment and survival of F. vesiculosus zygotes [25,26]. Regarding direct effects, elevated nutrient concentrations positively affect photosynthesis and growth in adult F. vesiculosus [27].
As a consequence of eutrophication in the Baltic Sea, sedimentation and bacterial re-mineralization are enhanced, leading to seasonal hypoxia and hypercapnia in sub-surface waters [28]. Further spread of (periodic) hypoxic areas in the Baltic Sea is predicted until the end of the 21 st century due to the interplay between eutrophication and lower oxygen solubility under warming [7]. Hypoxia induced mortality of benthic organisms was reviewed by Gray et al. [29] and recognized at a Baltic [30] and a global scale [31]. Also, higher frequencies of local upwelling in the last decades have been documented in the Baltic Sea [32].
In addition to hypoxia, warming was observed to impact growth, survival and photosynthetic efficiency in F. vesiculosus, while acidification showed weaker effects in the early [33] and adult life-stage [34]. As a consequence of warming, poleward range shifts of seaweed populations have been observed worldwide [14,35,36]. Also in the Baltic Sea, northward range shifts of F. vesiculosus populations have been observed and are predicted to continue during the 21 st and 22 nd century [37]. In contrast, ocean acidification was reported to increase growth in non-calcifying macroalgae [27,38]. Assumed reasons are enhanced availability of CO 2 to saturate the carbon demand during photosynthesis or the saved energy when carbon concentrating mechanisms are downregulated [39]. However, physiological responses of the macroalga Macrocystis pyrifera to acidified conditions showed that increased CO 2 conditions did not affect growth or photosynthesis [40]. Early life-stage macroalgae may be particularly threatened by global change [41]: Warming lowers the germination success of F. vesiculosus at 25°C [42] and reduces survival in F. serratus germlings more severely than in adults [43]. Enhanced sedimentation caused by eutrophication is more detrimental for young [23,44] than for adult Fucus [15,45]. Despite recognized differences in stress sensitivity among life-stages of a species, past research has mainly focussed on adult forms [46,47].
Adaptation of marine populations to global change is favoured by genetic diversity [48,49]. Conversely, low intraspecific genetic diversity and high phenotypic plasticity and gene flow restrict adaptation to environmental stress [50][51][52][53]. The increased tolerance of the eelgrass Zostera marina with higher genetic diversity to a summer heat wave [54] showed that genetic diversity may buffer warming stress at the population level [55].
Baltic F. vesiculosus populations show reduced genetic variation compared to Atlantic populations probably due to isolation and bottlenecks as well as the eroded genetic variation due to selection [56]. It has been argued that environmental stress, e.g. osmotic stress in the brackish Baltic, lower genetic diversity and limited dispersal capacity of F. vesiculosus gametes [57] may favour local extinctions [58]. However, this study is the first one to test genetic variation in Baltic F. vesiculosus with regard to the sensitivity towards OAW, nutrient enrichment and subsequent hypoxia.
The aim of the present study was to investigate (1) how OAW, (2) nutrient enrichment and (3) upwelling events affect the survival and growth of F. vesiculosus germlings, (4) how OAW interacts with simultaneous exposure to high nutrient concentrations, (5) how these treatments modify hypoxia sensitivity and (6) whether sibling groups vary in their tolerance towards these environmental parameters. Our experimental concept, hence, consisted in the exposure of genetically different sibling groups of F. vesiculosus germlings to increased temperature, pCO 2 and nutrient conditions and to a final upwelling event while maintaining the natural variations.

Material and Methods
Collection, gamete acquisition, experimental design A total of 64 fertile F. vesiculosus were collected in a wave exposed area with mixed hard substrate and sand bottom in the southwestern Baltic Sea (Bülk, Germany, 54°27.327´N, 10°11.977 W) in mid-June 2014. To avoid the collection of siblings and ensuring for genetic variability, individuals sampled were distanced by at least 2 meters, which is the estimated maximum dispersal distance of most F. vesiculosus gametes [57]. After collection, algae were transported to the lab in cooler boxes. Fertile receptacles were cut from these dioecious algae and gender was determined (46 females, 18 males) under the microscope at 100 x magnification (Olympus BH-2). Receptacles were rinsed with tap water, blotted dry and stored in the dark for 6 days at 14°C. Before gamete release, all receptacles from one female and one male individual (i.e. one parental pair) were put in a small dish. Gamete release followed by egg fertilisation was induced by immerging receptacles into sand-filtered seawater (15-16 psu) and exposing them to light (110 μmol photons m -2 s -1 ) for 3 hours. In this way, gametes were obtained from16 parental pairs. No specific permits were required for this study, the location is not privately-owned or protected and the study did not involve an endangered or protected species.
One mL of homogeneously suspended fertilised eggs was pipetted onto the upper surface of each of 2 x 2 cm sandstone cubes. Each cube with its settled germlings represented one experimental population. 16 different populations, each composed of germlings stemming from one parental pair, were thus produced. Culture and monitoring of germlings took place in a room with windows approximating natural light conditions during 3 weeks with weekly water exchange (15-16 psu) at 15°C until introducing them to the Kiel Outdoor Benthocosms (KOB). This facility maintains the natural in situ fluctuations and simultaneously allows manipulating environmental variables (e.g. temperature, pCO 2 ) on top of these fluctuations, i.e. "delta treatments". Target temperatures were obtained and maintained by computer controlled heaters and coolers, while acidification was achieved by increasing the atmospheric pCO 2 within the enclosed head space above the tanks by injecting pre-mixed gas to maintain an atmospheric pCO 2 of 1100 μatm. Details of the experimental set-up of the KOB and the logged conditions in the tanks are given in Wahl et al. [59].

OAW x nutrient and upwelling treatment
The OAW x nutrient experiment took place during 2 summer months from mid-July until mid-September 2014. One PVC box (70 cm x 40 cm x 12 cm) was placed within each of the 12 Benthocosm main tanks, each of them containing all 16 experimental populations of germlings. Since the upper rim of the lid-less boxes was a few centimetres above the water surface, the water body within the PVC boxes was separated from the water body of the main tank but open to the atmosphere. Thus, the boxes experienced the same treatments as the main tank regarding OAW but were insulated from gene flow from the adult Fucus population in the main tank. Twice a week the water of the PVC box was exchanged by water from the main tank which had been filtered through a 50 μm mesh to prevent the introduction of F. vesiculosus eggs (100 μm diameter). As the single-factor and combined effects of temperature and pCO 2 on F. vesiculosus germlings have been investigated previously [33], we combined warming and acidification into a single factor (OAW) in the present run. The two fully crossed factors OAW and nutrients were applied at two levels each (ambient and future). Ambient and predicted future levels of OAW and nutrients were simulated by adding the expected shift to the natural fluctuations of the ambient fjord conditions as delta treatments [59].
The "ambient" condition represents the natural fjord conditions transported into the main tanks of the KOB by a continuous flow-through (1 tank-volume per day, i.e. 1500 L/ 24 h) of Kiel Fjord water pumped from 1 m depth. "Future" conditions were simulated by adding 5°C to the actual temperature of the Kiel Fjord and by increasing the pCO 2 concentration in the hooded headspace of the tanks to 1100 μatm according to the predictions for the year 2110 for the Baltic Sea [2]. The bi-weekly nutrient enrichment (2x ambient) was achieved by doubling the "ambient" concentration which was taken as the seven years (2006-2013) mean for each specific date of nutrient addition (Table 1, S1 Fig). NaNO 2 (Merck, Germany), NaNO 3 (Carl Roth, Germany) and H 2 NaO 4 P.H 2 O (ACROS organics, Germany) were dissolved in fjord water 10 minutes before addition to the nutrient treatments. The ratio P: N of the Kiel Fjord is approximately 1: 1.5 and does not match the Redfield Ratio, probably due to the nutrient input of the nearby located river Schwentine. Additionally, high organism activity during summer months decreases the overall nutrient availability in shallow water. This ambient P: N ratio was not altered when NO 3 , NO 2 and PO 4 concentrations (μmol L -1 ) were doubled under "future" conditions. Our analysis on CN ratio in Fucus tissue (see result section below) showed no differences under "future" conditions, indicating that carbon did not become a limiting factor under nutrient enrichment.
The upwelling experiment was performed immediately after the end of the OAW x nutrient treatment phase, i.e. when all the treatments in the tanks were set back to "ambient". During Table 1. Nutrient concentration under ambient and future conditions. "Ambient" nutrient concentrations of PO 4 , NO 2 , NO 3 (μmol L -1 ) for the respective summer months and "future" nutrient concentrations (μmol L -1 ).

July
August September

Abiotic conditions
During the OAW x nutrient experiment, the bi-weekly measured oxygen concentrations (mean = 8.91 ± 0.38 mg L -1 ) varied in the low OAW treatment with a minimum of 8.01 ± 0.23 mg L -1 and a maximum of 9.44 ± 0.16 mg L -1 (mean ± SD) and within the high OAW treatment (mean = 7.76 ± 0.56 mg L -1 ) with a minimum of 6.75 ± 0.56 mg L -1 and a maximum of 8.72 ± 0.21 mg L -1 . Day-night fluctuations measured in two hours intervals showed that O 2 concentrations varied between a minimum of 7.64 mg L -1 at 4:30 am and a maximum of 9.97 mg L -1 at 16:30 pm under ambient conditions. Ambient temperatures (factor level 'ambient') (mean = 19.26 ± 2.38°C) varied between biweekly measurements with a minimum of 15.37 ± 0.19°C and a maximum of 22.7 ± 0.08°C. Elevated temperatures (factor level 'future') (mean = 23.64 ± 2.68°C) varied between a minimum of 19.53 ± 0.52°C and a maximum of 27.47 ± 0.26°C, which occurred during a natural summer heat wave (mean ± SD) (Fig 1). Temperatures during day and night fluctuated between a minimum of 16.1°C at 10:30 pm and a maximum of 17.7°C at 8:30 am under ambient conditions.

Response variables
Growth. For growth measurements, digital images were taken of 10-15 randomly chosen individual germlings per population at 40 x magnification (SteREO Discovery. V8 -Carl Zeiss Jena GmbH) similar to Steen and Scrosati [61]. Measured individuals were chosen randomly, since germlings were too small for labelling. The projected side-view of the single germlings was measured with the image analysis software Image J 1.45s (National Institutes of Health, USA) and the mean of germlings' area of the perpendicular projection was calculated for each population. Germlings' area was measured at the beginning (area t 0 ) and after 8 weeks (area t) at the end of the OAW x nutrient experiment. Relative growth rate (RGR) in % d -1 was calculated as exponential growth: Where Δt is the time period between t 0 and t in days. Survival. Germling number was counted under a binocular at 25x magnification between the start (number t 0 ) and after 8 weeks (number t) of the OAW x nutrient experiment. Survival of germlings was expressed as the percent of surviving germlings and calculated as: For determining survival during the final upwelling experiment, the germling number at the end of the preceding OAW x nutrient experiment was set as t 0 and number t was the germling number after the upwelling treatment.
Log-effect ratio. Log effect ratios were performed to show the direction and the strength of the sibling groups' phenotypical responses to the different OAW and nutrients enhancements relative to the respective ambient conditions. These differences among sibling groups' responses were determined in order to assess whether higher genetic diversity also increases the variance in responses. Sibling groups' sensitivity to high OAW and N were calculated at ambient and high OAW. The same procedure was used for calculating the log effect ratios for survival. Negative growth rates observed in 3 populations were attributable to the mortality of the bigger sized individuals and a decrease in mean area t compared to t 0 . To avoid negative values for log effect ratio analysis, growth was measured as ratios of area t and area t 0 as: CN ratios. Germlings of five sibling groups (1, 3, 7, 12 and 14) of the two treatment combinations OAW-N-and OAW-N+ were pooled for CN analysis. For the analysis of carbon and nitrogen content, freeze-dried algal material was ground to powder and three subsamples of 2 mg from each treatment was packed and loaded into tin cartridges (6×6×12 mm). Then, the packages were combusted at 950°C and the absolute C and N contents in % dry weight (% DW) were automatically quantified in an elemental analyser (Elementar Vario EL III, Germany) using acetanilide as standard according to Verardo et al. [62].

Statistical Analysis
Growth (RGR, % d -1 ) and survival (%) were analysed using a split-plot ANOVA with the fixed factors 'OAW' (with two levels: OAW-and OAW+) and 'Nutrients' (with two levels: N-and N +) as well as the random factors 'Mesocosm' (with 12 levels) and 'Sibling Group' (with 16 levels). 'Mesocosm' was nested in the 'OAW' x 'Nutrients' interaction. This model allowed us to analyse the influence of the fixed factors. Furthermore, we were able to identify possible random-by-fixed factor interactions between 'Sibling Group', 'OAW' and 'Nutrients'. Survival data were arcsine transformed prior to the analysis to overcome data truncation. Differences in carbon and nitrogen contents (% DW) as well as in CN ratios between five single sibling groups were identified in a further split-plot ANOVA with the fixed factor 'Nutrients' (two levels: see above), which was combined with OAW only, and the random factors 'Mesocosm' (12 levels) and 'Sibling Group' (16 levels), while 'Mesocosm' was nested in 'Nutrients'. Split-plot ANO-VAs were performed by using Satterthwaite's method for denominator synthesis [63], which calculates appropriate error terms for the F-ratios of the respective effects. Normality of errors and homogeneity of variances were verified by using residual plots (Q-Q Plot and Standardized Residuals Plot, respectively) in STATISTICA and R. The significance level of all analyses was α = 0.05. Post-hoc tests were performed using Tukey's HSD. All statistical analyses were conducted using the software STATISTICA v. 12 [64] and R.
CN ratio. Germling carbon and nitrogen content (% DW) as well as the CN ratio did not differ significantly among the nutrient treatments (Split-plot ANOVA, p-value > 0.05, Table 3). The CN ratio was 10.56 ± 1.13 under ambient nutrient conditions and 7.18 ± 3.06 (mean ± SD) under high nutrient conditions.   Survival. In all sibling groups with different treatment histories, the three days of hypoxia reduced germling survival significantly. Survival was highest for germlings pre-treated at ambient conditions (OAW-N-) and lowest for germlings previously treated at high OAW and ambient nutrients (OAW+N-) (Fig 5). Thus, high OAW significantly doubled the sensitivity to subsequent hypoxia as compared to a non-warmed, non-acidified regime (Split-plot ANOVA, OAW, p-value < 0.05, Fig 5, Table 2C). Nutrient addition during the preceding experiment enhanced sensitivity to hypoxia in germlings stemming from a regime of ambient temperature and CO 2 conditions but not for those from a OAW+ regime as reflected in the significant interaction between 'OAW' and 'Nutrients' (Split-plot ANOVA, OAW x Nutrients, p-value < 0.05, Fig 5, Table 2C).

Discussion
The simulated OAW as expected for 2110 strongly reduced survival and growth of F. vesiculosus germlings. The analysis of the single factors warming and acidification in a previous study at the KOB showed that warming was the main driver of mortality of F. vesiculosus germlings, while acidification played a minor role [33]. At temperatures > 27°C, reduction in growth and photosynthetic efficiency were observed in adult F. vesiculosus, finally resulting in necrosis [34]. Reduced growth of germlings under heat stress (> 25°C) was also observed in this work. Table 3. Nutrient effect on carbon and nitrogen content (% DW) and on the CN ratio. Split-plot ANOVA with the fixed factor 'Nutrient' and the random factors 'Mesocosm' and 'Sibling group'. Effects of the nutrient treatment were analysed for (A) the carbon content (% DW), (B) the nitrogen content (% DW) and (C) the CN ratio. Df: degrees of freedom, SS: sums of squares and MS: mean squares. 'Den. Syn. Error df' and 'Den. Syn. Error MS' describe the denominator synthesis of degrees of freedom and mean squares, respectively. This shows that the upper limits of thermotolerance of F. vesiculosus performance are similar in early and adult life-stages. However, in our multi-factorial design, the negative effects of warming on germling survival and growth rates were strongly mitigated by high nutrient concentrations. Earlier studies on Baltic adult F. vesiculosus have shown that nutrient enrichment increases nutrient uptake [65] and enhances photosynthetic efficiency [27]. Similar responses were also observed in other algal species. Ulva rigida cultured under nutrient enrichment reacted with higher nitrogen uptake, higher nitrate reductase activity and higher growth rates. Moreover, the nitrogen reductase activity was enhanced under future (1000 μatm) compared to ambient pCO 2 (400 μatm) conditions [66]. Similarly, high pCO 2 enhanced nitrogen assimilation in the brown alga Hizikia fusiforme [67] and may have decreased the relative investment in the nitrogenintensive protein biosynthesis [67]. Consequently, nitrogen may have been freed for other processes, such as growth and nitrogen storage [68]. Possibly, the F. vesiculosus germlings in our experiments also took up more nitrogen under nutrient enrichment when additional CO 2 was provided under acidified (and warmed) conditions, resulting in increased growth and survival. Our experimental design did not allow for disentangling the different possible mechanisms of mitigating effects of nutrients on either warming or acidification. Photosynthesis is regarded as one of the most heat sensitive metabolic activities in the plant cell [69,70], with at least three major heat-stress sensitive sites in the photosynthetic machinery: the photosystems (mainly photosystem II with its oxygen-evolving complex), the ATP generating and the carbon assimilation processes [71]. Moreover, respiration rates are increased under warming [72]. As our response variables growth and survival represent responses integrating over many metabolic processes, several compensatory effects caused by high nutrient levels appear possible. In conclusion, nutrient enrichment compensated to some degree the severe negative effects of future heat stress on F. vesiculosus germlings, which may be further mitigated by higher carbon availability under acidified conditions. However, such direct beneficial effects of nutrient enrichment may be overridden by indirect detrimental effects of eutrophication (such as increases in water turbidity, sedimentation, grazing and abundance of epibiotic filamentous algae) at field conditions, as reviewed by Berger et al. [45]. During the OAW x nutrient experiment, the epibiota under ambient and enriched nutrient conditions have not been determined. Regular filtration of the water content of our experimental boxes kept the fouling load relatively low under both, ambient and high nutrient conditions, assuming that epibiota had no strong effects on the F. vesiculosus germlings.

Source of variation
Increased nitrogen uptake is accompanied by higher CO 2 uptake (even at ambient pCO 2 ), hence a constant CN ratio is maintained [66]. Our findings show that the CN ratio in F. vesiculosus germlings was lower in the high nutrient treatment. Although this difference was not significant, it suggests higher uptake rates of nitrogen under nutrient enrichment. Since the nitrogen concentration in Baltic F. vesiculosus thalli is lowest in summer [73], nutrient enrichment effects may be most conspicuous in this season. This may have contributed to the observed stress-mitigating effect of nutrient enrichment in F. vesiculosus germlings. The buffering of stress impact by additional resources was also observed in juvenile blue mussel Mytilus edulis, when high food conditions enhanced the tolerance to ocean acidification [74].
The different sibling groups showed high variations in survival under warming and acidification, indicating the enhanced potential for adaptation in genetically diverse populations [48]. The crucial role in genetic variation for recovery from disturbances has also been reported in estuarine macrophytes [75] allowing for adaptation under global change stress [76].
The three day hypoxia experiment in the KOB simulating an upwelling event induced substantial germling mortality. During a local upwelling event, deep water with low oxygen concentration, low temperatures an increased pCO 2 and high salinity is shoaling [77,78]. In our upwelling treatment, temperature did not decrease considerably (16.41 ± 0.33°C) compared to previous values of 19.26 ± 2.38°C. Likewise, salinity (22.8 psu) did increase only slightly relative to the previous condition ( Fig 1C) and the elevated CO 2 (as associated with hypoxic upwelling) has minor effects on germling survival [33] (Fig 1C). This leaves the low oxygen concentration (2.75 ± 0.41 mg L -1 ) during the upwelling event as the most likely driver of germling mortality. Mortality during the three days of upwelling was considerably higher than during the two months of the preceding experiment, illustrating the high susceptibility of F. vesiculosus germlings to hypoxia. This susceptibility is probably due to reduced respiration rates under dark conditions, as it was also observed under hypoxic conditions in Cladophora vagabunda and Gracilaria tikvahiae [79]. Reduced respiration rates are accompanied by decreased provision of ATP and biosynthetic precursors leading to higher stress sensitivity [80] and to a reduced metabolism [79]. Susceptibility to hypoxia impacts was highest on germlings which previously experienced warmed and acidified conditions. Thus, the impact of hypoxic upwelling events in the future may be amplified by synchronous OAW. The assumed increased respiration under warming [72] may have further increased the O 2 debt, which could not be balanced under hypoxic conditions. Consequently, F. vesiculosus germlings grown under high compared to ambient temperatures were less tolerant to hypoxia. Thus, germlings grown under high nutrient levels experienced higher mortality under hypoxia compared to those grown at low nutrient levels. Zou et al. [68] demonstrated that under high-nitrogen conditions, respiration was enhanced by high CO 2 compared to ambient CO 2 conditions in the macroalga H. fusiforme. Consequently, in algae growing under nutrient enrichment increased respiration might be necessary to support higher maintenance demands (e.g. due to increased RUBISCO contents) and greater uptake of extra nitrogen [68]. This nutrient-driven higher metabolism may have rendered these germlings more susceptible to hypoxia.
In summary, the responses to hypoxia depended on the preceding OAW x nutrient treatments we applied. This suggests that there are different protective mechanisms in F. vesiculosus germlings that vary with the type of stressor. Future expansions of hypoxic areas in the Baltic Sea [7] will have severe effects on F. vesiculosus recruitments, as observed in this experiment, as well as on the benthic community in general [81]. We demonstrated that the net impact of global change including warming, acidification, eutrophication and hypoxia may depend on the interaction among these global and regional factors. This finding underscores the importance for analysing the combined effects of multiple stressors and their interconnectivity for accurate predictions of future scenarios [3]. Moreover, the indirect effects of global change may be more significant than the direct effects [13]. Scaling up multiple stressors is crucial for predicting the fate of F. vesiculosus populations [82].
Supporting Information S1 Fig. Nutrient concentrations in the experimental boxes. Nutrient concentrations within the experimental germling boxes before (A, B) and after (C, D) the bi-weekly addition of the nutrients (PO 4 , NO 3 , NO 2 ) as well as initial NH 4 conditions in μmol L -1 in July and August. Initial nutrient concentrations were measured six times per month in the main KOB tank before the water addition to the boxes. Nutrient concentrations after additions were determined by adding the sum of the initial and additional nutrient concentration. (TIFF)