Articles | Volume 18, issue 8
Research article
22 Apr 2021
Research article |  | 22 Apr 2021

Decoupling salinity and carbonate chemistry: low calcium ion concentration rather than salinity limits calcification in Baltic Sea mussels

Trystan Sanders, Jörn Thomsen, Jens Daniel Müller, Gregor Rehder, and Frank Melzner

The Baltic Sea has a salinity gradient decreasing from fully marine (> 25) in the west to below 7 in the central Baltic Proper. Habitat-forming and ecologically dominant mytilid mussels exhibit decreasing growth when salinity < 11; however, the mechanisms underlying reduced calcification rates in dilute seawater are not fully understood. Both [HCO3-] and [Ca2+] also decrease with salinity, challenging calcifying organisms through CaCO3 undersaturation (Ω≤1) and unfavourable ratios of calcification substrates ([Ca2+] and [HCO3-]) to the inhibitor (H+), expressed as the extended substrate–inhibitor ratio (ESIR). This study combined in situ monitoring of three southwest Baltic mussel reefs with two laboratory experiments to assess how various environmental conditions and isolated abiotic factors (salinity, [Ca2+], [HCO3-] and pH) impact calcification in mytilid mussels along the Baltic salinity gradient. Laboratory experiments rearing juvenile Baltic Mytilus at a range of salinities (6, 11 and 16), HCO3- concentrations (300–2100 µmol kg−1) and Ca2+ concentrations (0.5–4 mmol kg−1) reveal that as individual factors, low [HCO3-], pH and salinity cannot explain low calcification rates in the Baltic Sea. Calcification rates are impeded when Ωaragonite 1 or ESIR  0.7 primarily due to [Ca2+] limitation which becomes relevant at a salinity of ca. 11 in the Baltic Sea. Field monitoring of carbonate chemistry and calcification rates suggest increased food availability may be able to mask the negative impacts of periodic sub-optimal carbonate chemistry, but not when seawater conditions are permanently adverse, as observed in two Baltic reefs at salinities < 11. Regional climate models predict a rapid desalination of the southwest and central Baltic over the next century and potentially a reduction in [Ca2+] which may shift the distribution of marine calcifiers westward. It is therefore vital to understand the mechanisms by which the ionic composition of seawater impacts bivalve calcification for better predicting the future of benthic Baltic ecosystems.

1 Introduction

1.1 Baltic Sea hydrochemistry

The Baltic Sea is a semi-enclosed brackish water body with a decreasing salinity (S) gradient from > 25 in the Kattegat to salinities < 3 in the Gulf of Bothnia and the Gulf of Finland (Meier, 2006; Neumann, 2010; Fig. 1). Salinity and seawater carbonate chemistry are highly variable on both spatial and temporal scales with extreme fluctuations resulting from occasional inflow events of highly saline North Sea water, wind-driven upwelling of hypoxic and hypercapnic deep water, lateral transport of water masses, and pronounced regional salinity gradients (Thomas and Schneider, 1999; Melzner et al., 2013; Saderne et al., 2013; Mohrholz et al., 2015). Along the Baltic Sea salinity gradient, seawater total alkalinity (AT) decreases linearly from  2300 µmol kg−1 in the North Sea as the marine endmember to  1600 µmol kg−1 at a salinity of 6 in the Baltic Proper (Müller et al., 2016). In seawater that is near equilibrium with current atmospheric pCO2, the dissolved bases HCO3- and CO32- comprise ca. 96 % of AT. These two carbon species constitute ca. 99 % of all total dissolved inorganic carbon (CT) and subsequently both AT and CT are tightly coupled (Zeebe and Wolf-Gladrow, 2001). Along with CT, seawater [Ca2+] also decreases linearly with salinity from the North Sea to the Baltic Proper (Kremling and Wilhelm, 1997). This decreasing availability of seawater Ca2+ and CT puts pressure on calcifying organisms that extract these 2 substrates (Ca2+ and CT) from seawater for calcification.

Figure 1Map of the southwest and central Baltic Sea depicting the decreasing salinity gradient from the Kattegat (west) to the Baltic Proper (east) with the three field sites shown (Kiel, Ahrenshoop and Usedom). The inset map shows the position of the southwest Baltic Sea (box in inset) in relation to other regions: (1) Kattegat, (2) southwest Baltic, (3) Gulf of Riga, (4) central Baltic Proper, (5) Gulf of Finland and (6) Bothnian Sea. Salinity data were taken from open-source monitoring data (EU, Copernicus Marine Services, 2018; details in the Supplement).

1.2 Ecophysiology of calcifying mussels

Despite reduced [Ca2+] and CT along the Baltic Sea salinity gradient, calcifying mussels in the genus Mytilus dominate the Baltic benthos particularly below salinities of 11, being the primary biogenic reef-forming species in the central and eastern Baltic basins (Westerbom et al., 2002). In the Baltic Sea, these biogenic mussel reefs contribute more to ecosystem services than both seagrass meadows and macroalgae beds, being responsible for biodiversity hotspots and significant carbon and nutrient recycling in the benthos (Norling and Kautsky, 2008; Koivisto and Westerbom, 2010; Heckwolf et al., 2021; Attard et al., 2020). Models suggest that the contribution of filter feeding mussels to the region's biogeochemistry is so pronounced that increased bivalve aquaculture may significantly reduce levels of eutrophication in the central Baltic Sea (Kotta et al., 2020). Despite the ecological importance of Baltic Mytilus, calcification and growth rates of these mussels are drastically reduced below salinities of 16–11, resulting in extremely slow-growing mussels and reduced maximum body size (Vuorinen et al., 2002; Riisgård et al., 2014; Sanders et al., 2018; Westerbom et al., 2019). However, high abundances of mussels at salinities of 11–6 facilitate levels of ecosystem function comparable to that at higher salinities of > 20 (Elmgren and Hill, 1997). Below a threshold salinity of < 5, Baltic Mytilus can no longer survive and ecosystem function is severely reduced (Elmgren and Hill, 1997; Vuorinen et al., 2015).

Steep declines in Baltic Mytilus calcification below salinities of 11 have been attributed to unfavourable protein metabolism from osmotic stress and high energetic costs of CaCO3 biomineralization related to sub-optimal seawater chemistry (Tedengren and Kautsky, 1986; Maar et al., 2015; Sanders et al., 2018; Thomsen et al., 2018). Despite the challenging environment, high levels of eutrophication resulting in increased phytoplankton food availability has been shown to mask the negative impacts of sub-optimal carbonate chemistry on calcification in Baltic Mytilus in Kiel Fjord (Melzner et al., 2011; Thomsen et al., 2013). Similarly, macrophyte habitats can mitigate the negative effects of seawater acidification on calcifiers by localized pH increases in adjacent seawater (Wahl et al., 2017). The complexity of environmental control on calcification in Baltic mussels makes it difficult to predict how environmental change will impact calcifying mussels and associated ecosystem services.

1.3 Carbonate chemistry and calcification

Since salinity, [Ca+] and CT are tightly coupled in the Baltic Sea it is difficult to decipher which of these factors primarily controls calcification. The stability of CaCO3 structures in seawater is dependent on the saturation state (Ω) of the CaCO3 polymorph in question (predominantly aragonite or calcite in bivalves), which is defined as follows:

(1) Ω = Ca 2 + CO 3 2 - K sp ,

with Ksp being the temperature-, pressure- and salinity-dependent solubility product of the CaCO3 polymorph under consideration. CaCO3 dissolution becomes more thermodynamically favourable when Ω< 1 (undersaturation), which has major implications for marine organisms which biomineralize external CaCO3 structures (Waldbusser et al., 2014; Ries et al., 2016). Extended periods of undersaturation of the CaCO3 polymorph aragonite (Ωaragonite≤1) are commonplace in the Baltic Sea, theoretically imposing constraints on calcifying organisms that precipitate aragonite CaCO3 crystals, as mytilid mussels do for production of inner shell nacre (Tyrrell et al., 2008; Melzner et al., 2011). While Ωaragonite is a good predictor of biological calcification especially at fully marine conditions, there is little evidence supporting its mechanistic role in calcification physiology. Even with a rudimentary understanding of calcification mechanisms in molluscs, we know that increased [H+] inhibits extracellular calcification by diminishing the proton gradient between the slightly alkalized calcifying fluid in larvae and ambient seawater (Ramesh et al., 2017). Additionally, HCO3- rather than CO32- is the primary source of inorganic carbon for calcification in marine organisms (Bach, 2015). This is because in seawater at equilibrium with current atmospheric pCO2, HCO3- is ca. 10-fold more abundant in seawater than CO32-, and HCO3- transporters have been localized in the calcifying tissues of marine organisms (Roleda et al., 2012; Zoccola et al., 2015; Fassbender et al., 2016; Hu et al., 2018). Since the availability of HCO3- stimulates calcification, a more accurate and mechanistically relevant predictor for calcification is the substrate (HCO3-) inhibitor (H+) ratio or SIR (Jokiel, 2013; Bach, 2015; Thomsen et al., 2015; Cyronak et al., 2016). The effects of [Ca2+] are generally excluded in calculating SIR (Eq. 2) as [Ca2+] is high and generally constant at oceanic salinity (∼10 mmol kg−1).

(2) SIR = [ HCO 3 - ] [ H + ]

At salinities < 8, however, limitation of seawater Ca2+ has been shown to reduce calcification in larval Baltic Mytilus (Thomsen et al., 2018). Therefore, adopting an extended SIR (ESIR) to include [Ca2+] may provide an even more accurate predictor of calcification under low-salinity conditions.

(3) ESIR = [ Ca 2 + ] [ HCO 3 - ] [ H + ]

Previous work reveals that calcification in Baltic larval bivalves drops significantly when ESIR  0.7 (Thomsen et al., 2018). Application of the ESIR is particularly important in coastal and estuarine zones such as the Baltic, where freshwater input, photosynthetic activity and wind-driven upwelling can delineate carbonate system parameters (pH and pCO2) from salinity and weaken the predictive power of SIR for calcification (Melzner et al., 2013; Fassbender et al., 2016). Severe desalination and warming are predicted in the Baltic Sea over the next 50–100 years, due to changes in weather patterns and precipitation, yet our ability to predict future carbonate chemistry changes remains somewhat limited due to potential antagonistic effects of desalination and continental weathering (Gräwe et al., 2013; Müller et al., 2016). The combination of our poor understanding of (1) how Baltic Sea carbonate chemistry will change in the future and (2) how environmental factors individually control calcification makes it difficult to predict the impacts of environmental change on Baltic Sea calcifiers and, by extension, Baltic Sea ecosystems.

1.4 Scope of this study

This study addresses the impacts of environmental change on Baltic Sea calcifying mussels, by combining a multi-year environmental monitoring programme and field growth study and at three Baltic Mytilus reefs and two complementary laboratory experiments. Field salinity, food availability, carbonate chemistry and calcification rates were monitored in three mussel reefs along the southwest Baltic Sea salinity gradient. This was complemented by laboratory experiments in which salinity was decoupled from calcification substrate availability (HCO3- and Ca2+) and pH to investigate how these individual factors control calcification rate. The results from this study shed light on the primary environmental factors driving growth rates in ecologically dominant calcifying Baltic mussels. This will support biogeochemical model projections to better predict changes in distribution and function of Baltic calcifying bivalves in the future.

2 Material and methods

2.1 Environmental monitoring

Three monitoring sites were selected along the southwest Baltic Sea coast based on previous population comparison studies spanning the geographic range of the steepest section of the Baltic Sea salinity gradient (Meier, 2006; Sanders et al., 2018; Fig. 1). These sites were Kiel (541949′′ N, 10860′′ E), Ahrenshoop (54237′′ N, 122524′′ E) and Usedom (54321′′ N, 14040′′ E). Mini CTD (conductivity, temperature, depth) loggers (Star-Oddi, Iceland), logging data every 3 h, were deployed on wooden pillars at each site between July and September 2015 (0.5–1.5 m depth) and replaced every 2–3 months until January 2018. Field [Ca2+] was calculated from salinity values using [Ca2+]–chlorinity relationships from Kremling and Wilhelm (1997) and a chlorinity–salinity relationship from Millero (1984). Reference measurements from field water samples concurred with these calculations from the literature (Fig. S1 in the Supplement). Additionally, every 2–3 months, two 500 mL water samples were taken from Ahrenshoop and Usedom (1 m depth) and immediately poisoned with saturated HgCl2 solution (0.5 ‰ final concentration) for pHtotal, AT and CT measurements at the Leibniz Institute for Baltic Sea Research, Warnemünde. Kiel carbonate chemistry data were taken from published monitoring data for January–December 2015 (Hiebenthal et al., 2017). Water pHtotal was measured spectrophotometrically at 25 C with non-purified m-cresol purple as indicator dye (pH calculated according to Müller and Rehder (2018), dye pH perturbation corrected after Hammer et al., 2014). CT was determined with a SOMMA system (Single Operator Multi-Parameter Metabolic Analyzer) operated at 15 C and AT was determined by open-cell titration at 20 C (Dickson et al., 2007). Field carbonate chemistry parameters were calculated using CT and pHtotal as inputs using the CO2Sys_v2.1 program (KHSO4 constants from Dickson, 1990; K1 and K2 dissociation constants from Millero, 2010).

2.2 Field calcification

In mussel populations originating from the three monitoring sites chosen in this study, calcification rates have been shown to drop significantly below salinities of 11 (Sanders et al., 2018). To investigate the underlying cause of this in the field, mussel settlement structures (n= three per site) were placed at each site in March 2016 to monitor the growth of the spring cohort of settled juveniles over the course of 15 months. The structures were made from a 35 cm long cylinder of grey PVC pipe with a diameter of 10 cm. In total, 20 × 2 cm diameter holes were drilled into the sides and a 0.2 mm gridded nylon net (standard seed sock for mussel aquaculture) was placed diagonally across the inside to increase surface area allowing for maximum settlement (Fig. S2). These settlement structures were cable-tied to wooden pillars at each site at ca. 0.5–1 m depth, and samples of 50 juvenile mussels were taken for shell length measurements from each site every 2–3 months. Individual shell mass (SM, µg CaCO3) was calculated from individual shell lengths using specific SL–CaCO3 relationships developed for each of the three populations (Sanders et al., 2018; Fig. S3). To allow direct comparison between sites, individual calcification rates (µg CaCO3 d−1) are expressed as the linear slopes of CaCO3 mass over time (days) over the first three sampling time points before growth rates deviate from linear.

2.3 Food availability

Chlorophyll-a (chl-a) concentrations (µg L−1) were used as a proxy for phytoplankton food availability in Baltic Sea surface waters (Wasmund et al., 2011). Values for chl-a concentrations at the three monitoring sites were obtained from field monitoring data supplied by the Landesamt für Landwirtschaft, Umwelt und ländliche Räume (LLUR) sampled from Mönkeberg (54217′′ N, 101036′′ E), 2.5 km from the field monitoring site in Kiel Fjord. Chl-a monitoring data for Ahrenshoop and Usedom were supplied by the Landesamt für Umwelt, Naturschutz und Geologie (LUNG) sampled from near Fischland (542043′′ N, 122655′′ E), 1 km from the sampling site in Ahrenshoop and Zinnowitz (54457′′ N, 135457′′ E), 4 km from the monitoring site in Usedom.

2.4 Laboratory experiments: animal collection and maintenance

A genetic hybrid zone and gradient in allele frequencies exists in the southwest Baltic Sea from Baltic Mytilus edulis-dominated hybrids at higher salinities (> 16) to Baltic Mytilus trossulus hybrids at low salinities (< 8), with a genetic transition zone occurring at salinities of 10–12 (Stuckas et al., 2017). Two cohorts of freshly settled Baltic Mytilus from the genetic transition zone (mean shell length SL = 0.6 mm for cohort 1 and 3.0 mm for cohort 2) were collected from wooden pillars in Ahrenshoop, Germany (542313′′ N, 122537′′ E, mean salinity = 10.8, temperature = 19.3 C, depth < 1 m) in July 2014 and 2016 (Fig. 1). Animals were transported within 6 h in chilled, aerated cool boxes to GEOMAR, Kiel, stored in 20 L plastic acclimation aquaria (16 C) for 3 weeks prior to experiments and fed twice daily with live Rhodomonas salina to achieve phytoplankton concentrations of > 10 000 cells mL−1, which is considered saturated food availability for Baltic Mytilus (Riisgård et al., 2013). Water was exchanged in acclimation aquaria every 3 d and the correct salinity was obtained by mixing 0.22 µm filtered seawater (FSW) from Kiel Fjord (salinity  16) with distilled water and adjusting AT to site-specific values via addition of 1 M NaHCO3 solution.

2.5 [HCO3-] and [Ca2+] manipulation experiments

Two laboratory experiments were conducted to study the impacts of calcification substrate availability (HCO3- and Ca2+) on calcification. The first experiment exposed Baltic Mytilus juveniles (cohort 1, mean SL = 0.6 mm) to three salinities (16, 11 and 6, reflecting mean salinities at the three monitoring sites) and six bicarbonate ion concentrations (300, 600, 900, 1500, and 2100 µmol kg−1), reflecting AT of full-strength seawater down to salinities < 3 in the northern Baltic Sea (Müller et al., 2016). The experiment was conducted in 2 L plastic aquaria, each containing  1600 juvenile mussels with four replicate aquaria per treatment (15 treatments, 60 aquaria in total). Experimental aquaria were continuously aerated with pressurized ambient air throughout the 70 d experimental duration. Salinity and carbonate chemistry manipulation was carried out by diluting FSW from Kiel Fjord to experimental salinities (16, 11 and 6) with distilled water in 800 L containers with CT and pHNBS determined weekly in this stock FSW (Sect. 2.6). To achieve experimental [HCO3-] values, 1 M NaHCO3 or 1 M HCl was added to increase or decrease [HCO3-], respectively (mean water parameter values in experimental aquaria for the duration of the experiment are given in Table 2). Water was prepared individually for each treatment aquarium 24 h before water changes and equilibrated with pressurized air at atmospheric pCO2. Water changes were conducted twice weekly for the first 30 d and 3 times weekly for the latter 40 d to prevent significant deviations in target water chemistry conditions as animals grew. Dead animals were removed and counted during water changes with declines in animal abundances accounted for when calculating volume of food addition. CT and pHNBS were determined in experimental aquaria immediately before and after water changes to ensure that [HCO3-] and pH did not deviate by > 200 µmol kg−1 or > 0.1 units, respectively, due to biological activity or the addition of phytoplankton food culture. Phytoplankton food was added twice daily during the first 30 d and 3 times daily during the latter 40 d to compensate for increasing biomass, ensuring saturated availability of food during experiments.

A second experiment was conducted to investigate the impact of [Ca2+] on calcification rates. This experiment utilized the same experimental design as the first experiment, exposing juvenile Baltic Mytilus (cohort 2, mean SL = 3.0 mm) to the same three salinities (16, 11 and 6) and five different [Ca2+] values (0.5, 1, 2, 3, 4 mmol kg−1) covering the natural range of [Ca2+] in the southwest and central Baltic (Thomsen et al., 2018). Experimental seawater was manipulated through the preparation of Ca2+ free artificial seawater (CFASW) in three stock solutions (salinity 6, 11 and 16) by adding NaCl, NaSO4, KCl, NaHCO3, KBr, H3BO3, MgCl2 and SrCl2 to deionized water (Table S1 in the Supplement; Kester et al., 1967). Aquarium volumes in the Ca2+ experiment could not be as high as in the HCO3- experiment due to an inability to produce large volumes of CFASW. Consequently, the Ca2+ experiment was conducted in smaller 50 mL plastic aquaria with two animals per tank and four replicates (15 treatments, 60 aquaria in total) and lasted 37 d. Fewer animals were used per aquarium to minimize the impact of biomineralization on seawater carbonate chemistry due to smaller water volumes and larger experimental animals (Table S2). No mortality was observed in any aquaria during this experiment. Stock CFASW was prepared 24 h before water changes and continuously aerated with pressurized ambient air to equilibrate water with atmospheric pO2 and pCO2. For preparation of experimental water, [Ca2+] was adjusted through the addition of a 1 M CaCl2 stock solution to achieve target Ca2+ concentrations. CT and [Ca2+] were measured twice weekly in stock CFASW (Sect. 2.6). [Ca2+] was measured using a flame photometer (EFOX 5053, Eppendorf) calibrated with urine standards (Biorapid; see Table 3 for mean water parameters in experimental treatments for the duration of experiment). Aquaria were not aerated during experimentation, but twice-weekly monitoring of pHNBS and CT before and after water changes revealed mean pH did not change by more than 0.1 pH units and [Ca2+] did not deviate by more than 0.2 mmol kg−1. Twice-weekly water changes were sufficient to prevent respiratory build-up of CO2 between water changes (Tables 2 and 3), and animals were fed once daily with R. salina to ensure comparable levels of food per unit body mass between both experiments.

While it may be argued that the disparity in aquarium volumes and animal numbers makes comparisons between both experiments more difficult, measures were taken to ensure comparable conditions between both experiments (Table S2). Both laboratory experiments used juvenile mussels from the same sample population at the same developmental stage. Regular monitoring of carbonate chemistry in experimental aquaria ensured the frequency of water changes was sufficient to minimize the impacts of both biological activity and phytoplankton food culture addition on target carbonate chemistry values. Live microalgal food (R. salina) was added at pre-determined frequencies in both experiments to maintain saturated feeding conditions and a comparable level of daily food ration (no. of phytoplankton cells mg−1 body mass) between both experiments (Table S2). As biomass changed throughout the course of both experiments either through growth or mortality, feeding frequency was adjusted to ensure comparable concentrations of phytoplankton cells per milligram of body mass (Table S2) across both experiments and between treatments. Cell concentrations in phytoplankton cultures were monitored daily using a Multisizer 3 Coulter Counter (Beckman, Germany). In both experiments, aquaria were held in 16 C water baths and mean body mass per litre in experimental aquaria (average over duration of experiments) was comparable between both experiments (Table S2), minimizing the differential impacts of biological activity on experimental seawater conditions. Finally, unlike many other studies, after dilution of FSW with distilled water, experimental AT was adjusted to site-specific levels across all salinities to accurately imitate natural seawater conditions present in low-salinity (< 11) Baltic Sea habitats.

2.6 Carbonate chemistry

For carbonate chemistry measurements in experiments at GEOMAR, pHNBS was measured 3 times weekly using a WTW pH 3110 probe at 25 C calibrated with two standard precision NBS buffers (Radiometer analytical) at pH 7.0 and 10.0. Experimental temperature and salinity were monitored daily using a WTW Cond 315i probe. CT was analysed using an AIRICA CO2 analyser (Marianda, Kiel, Germany) calibrated by measuring certified reference material (Dickson et al., 2003). Experimental seawater carbonate system parameters (AT, [H+], [HCO3-], [CO32-] and Ωaragonite) were calculated with CT and pHNBS as inputs using the same methods described in Sect. 2.1. Values for Ωaragonite in the calcium experiment were calculated according to Eq. (1) from measured [Ca2+] in each technical replicate, as well as calculated [CO32-] and Ksp. The SIR and ESIR were calculated according to Eqs. (2) and (3).

2.7 Experimental calcification rates

Individual calcification rates in laboratory experiments were determined by calculating the change in CaCO3 mass (a proxy for shell mass; SM) over time between two time points (t0 and t1). Individual SM (mg) was calculated separately for each experiment at the beginning of experiments (t0) from mean SL per replicate tank and using a population-specific SL–CaCO3 mass relationship for the experimental population (Sanders et al., 2018; Fig. S3). Individual SM at the termination of each experiment (t1) was calculated by removing body tissue from individuals using forceps under a stereomicroscope and placing empty shells individually into pre-weighed and pre-dried aluminium foil cups which were then placed in a muffle furnace and ashed 450 C for 4 h to remove all organic content. The remaining inorganic mass was used as a proxy for CaCO3 mass. Individual calcification rates were calculated with the following equation:

(4) Calcification rate ( µ g CaCO 3 d - 1 ) = ( SM ) t 1 - ( SM ) t 0 no . of days / 1000 ,

where SM = shell mass (mg) at time points t0 and t1 represents SM at the beginning and end of each experiment, respectively, and number of days represents the duration of the experiments – 70 and 37 d for the HCO3- and Ca2+ experiments, respectively.

2.8 Data analysis

All data analysis was conducted using R Software version 3.6.3 (R Core Team, 2020) with all packages used listed in Table S3. Data distributions were tested for normality (Shapiro–Wilk test) and homogeneity of variances (Levene's test) before running parametric tests. In the bicarbonate experiment, a negative exponential decay model was fit to the relationship between [HCO3-] and calcification rate using a nonlinear least squares method for parameter optimization. The Akaike information criterion (AIC) was used to designate the most parsimonious model (Table S5). Linear regression models were fit to calcification rates with [Ca2+] as a co-variate and salinity as a factor in the calcium experiment. Negative exponential decay models were fit to explain the relationship between ESIR and Ωaragonite on laboratory calcification rates with model parameters (Cmax and K) compared (showing 95 % confidence intervals) between ESIR and Ωaragonite as predictors of calcification. The residual sum of squares (RSS) was used as a measure of model fit. The impacts of salinity were analysed statistically by comparing the Cmax parameter (±95 % confidence interval) of the non-linear model in the bicarbonate experiment and ANCOVA in the calcium experiment. Field calcification rates (monthly means) and environmental parameters (salinity, pH, chl-a, [Ca2+], [HCO3-], Ωaragonite, and ESIR) were analysed using parametric tests (ANCOVA) with time as a co-variate for field calcification rates and site as a fixed factor. Pairwise comparisons were conducted on significant factors by way of a Tukey test. Field temperatures for the given monitoring period (2015–2017) were analysed non-parametrically using a Kruskal–Wallis test. Finally, field calcification rates were calculated as the values for the linear slopes of cumulative calcification over time for the first three sampling periods. Log-transformed mean calcification rates for each site were plotted against mean values for individual field environmental parameters.

3 Results

3.1 Environmental monitoring

Salinity was different between all three monitoring sites (ANCOVA, F(2, 17 253)= 38 518, P<0.001, Table S6), being highest and more variable in Kiel and the lowest and most stable in Usedom (Fig. 2). As salinity and [Ca2+] exhibit a linear relationship (Fig. S1), [Ca2+] also differed to the same degree as salinity between the three sites (Fig. 2). Calculated mean [HCO3-] at the three field sites (Table 1) were significantly higher in Kiel compared to Ahrenshoop and Usedom (ANOVA, F(2, 55)= 38.80, P< 0.001, Table S6) with values never dropping below 1600 µmol kg−1 at any of the three sites during the monitoring period (Fig. 3c). While pHtotal was highly variable, particularly in Kiel (Fig. 3a), no difference was observed between all 3 sites (ANOVA, F(2,55)=1.217, P=0.304). Salinity and AT exhibited a linear relationship similar in slope to the SAT relationship reported by Müller et al. (2016) for the Kattegat in 2014 (Fig. 4). However, AT values measured in this study were ca. 100 µmol kg−1 higher at a given salinity than those reported by Müller et al. (2016) for the Kattegat and central Baltic Proper (model parameters given in Table S7). Calculated Ωaragonite was higher in Kiel compared to the other two sites (ANOVA, F(2,55)=7.22, P=0.002, Fig. 3e) with periods above and below Ωaragonite=1 in Kiel (mean = 1.17) and mean values below saturation (Ωaragonite < 1) in Ahrenshoop or Usedom at 0.83 and 0.76, respectively (Fig. 3e, Table 1). Similar patterns where observed in the calculated values for ESIR at the three sites (Fig. 3g) with Kiel exhibiting significantly higher mean ESIR values than both Ahrenshoop and Usedom (ANOVA F(2,55)=10.88, P< 0.001), with a mean value of 1.05±0.6 (Table 1). Mean ESIR values were below the threshold of 0.7 proposed by Thomsen et al. (2018) in Ahrenshoop and Usedom (means of 0.65 ± 0.19 and 0.55 ± 0.26, respectively), although periods above and below this threshold were observed at all sites (Fig. 3g). Mean chl-a concentrations at each site (Fig. 3i, Table 1) derived from monitoring data were significantly lower at Ahrenshoop compared to Kiel and Usedom (ANOVA, F(2,77)=13.8, P< 0.001, Table S6) and not significantly different between Kiel and Usedom (Tukey post hoc, P=0.35).

Table 1Summary of environmental conditions and field calcification rates for the experimental monitoring period (2015–2018) at the three Baltic Sea localities (Kiel, Ahrenshoop, Usedom). Salinity and temperature were recorded using mini-CTDs, and pHtotal, CT and AT were measured in laboratory conditions with pCO2, [HCO3-] and Ωaragonite calculated on the CO2Sys_v2.1 programme using pHtotal and CT as inputs. [Ca2+] and the ESIR were calculated from salinity and carbonate chemistry parameters, respectively. The chl-a data presented originate from physical monitoring data (Sect. 2.3; Kiel – 28 samples, Ahrenshoop – 25 samples, Usedom – 25 samples). Calcification rates represent the linear slope coefficients of calcification over time (days) at each site with linear model parameters given in Table S10. Mean values are presented ± standard deviation (below).

Download Print Version | Download XLSX

Figure 2Field salinity and [Ca2+] at the three Baltic Sea sites from August 2015–December 2017 (a). Salinity data are derived from deployed CTD loggers and [Ca2+] was calculated from these values (Sect. 2.1). Boxplots depict median salinity and interquartile ranges excluding outliers (b).


Figure 3Environmental monitoring data from Kiel, Ahrenshoop and Usedom. Time series for pHtotal (-log10[H+]) at in situ temperatures (a) and associated boxplots (b); [HCO3-(c), and boxplots (d); Ωaragonite with the horizontal dashed line depicting the theoretical saturation threshold of Ωaragonite=1 (e), and boxplots (f); extended substrate–inhibitor ratio ([Ca2+][HCO3-] / [H+]) with the horizontal dashed line depicting the theoretical saturation threshold of ESIR = 0.7 (g) and boxplots (h); chl-a as a proxy for food availability (i) and boxplots (j). Details describing carbonate chemistry calculations are described in Sect. 2.1. Boxplots display median values and interquartile ranges excluding outliers.


Figure 4ATS relationship from field monitoring data over the monitoring period for each site showing a linear fit with model parameters given in Table S7. Linear ATS relationships (grey dashed lines) for the Gulf of Bothnia, Gulf of Finland and the Kattegat are also depicted, using 2D linear model parameters for 2014 from Müller et al. (2016). The linear ATS relationship for the Gulf of Riga is shown using model parameters for 2008 from Beldowski et al. (2010).


3.2 Calcification rates

Calcification rates in the bicarbonate experiment exhibited a significant negative exponential decay relationship with decreasing [HCO3-] (Table S8), with calcification rates decreasing most abruptly below ca. 1000 µmol kg−1 (Fig. 5a) at salinities of 11 and 16. Maximum calcification rates (Cmax) at salinity 6 were only  13 % of those at 11 and 16 and not different between salinities of 11 and 16 (Fig. 5c). The subsequent experiment, which decoupled [Ca2+] from salinity, revealed a significant linear decrease (parameters given in Table S8) in calcification rate with [Ca2+] (ANCOVA, F(2,54)=106.9, P< 0.001) across the range of experimental [Ca2+] (Fig. 5b). Salinity as a single factor had no significant impact on calcification rates (ANCOVA, F(2,54)=0.83, P=0.442); however, a significant interaction effect between salinity and [Ca2+] suggests that the effects of [Ca2+] on calcification rate vary at different salinities. Combining both experiments, calcification rates did not exhibit any significant correlation with [Ca2+] or [HCO3-] alone under the given conditions (Fig. S4, Table S9). Calcification rates plotted separately against Ωaragonite and ESIR across both experiments revealed statistically significant estimations for both parameters (Cmax and K) in the negative exponential decay model (Fig. 6, Table S6). Estimated model parameters were not significantly different between Ωaragonite and ESIR as calcification predictors, (95 % CI; Fig. S6), although, a lower residual sum of squares (residual SS = 2866) in the ESIR model compared to the Ωaragonite model (residual SS = 3159) indicates slightly better performance of the ESIR model as a predictor for calcification in the southwest and central Baltic Sea. In the bicarbonate manipulation experiment, mean experimental [HCO3-] values were within 200 µmol kg−1 of target values except for 1 out of 15 treatments (Table 2). In the calcium manipulation experiment, [Ca2+] was within 0.2 mmol kg−1 of target values with the exception of 2 out of 15 treatments (Table 3). Despite different aquarium volumes and animal numbers per aquarium, mean biomass per litre in experimental aquaria in both experiments were comparable at 24.1 and 13.2 mg L−1 for the bicarbonate and calcium ion manipulation experiments, respectively (Table S2). Mean numbers of phytoplankton cells per unit biomass were also comparable between both experiments (Table S2). Despite experimental aquaria not being aerated during the calcium manipulation experiment, mean pH values fluctuated by 0.04–0.12 units between water changes and were comparable to field pH values at all three monitoring sites (Tables 1 and 3) suggesting minimal impacts of animal respiration on carbonate chemistry during the experiment. Calcification rates were comparable between both laboratory experiments (Fig. 5) and with field calcification rates, at least for Ahrenshoop (mean salinity: 10.9) and Usedom (mean salinity: 7.0) populations (Table 1). Mussels from Kiel (mean salinity: 15.2) exhibited significantly higher calcification rates (ANCOVA, F(2,12)=570, P< 0.001) compared to the two other field sites, being 1–2 orders of magnitude higher (Fig. 7). This contrasts with measured laboratory calcification rates where salinities 16 and 11 treatments exhibited comparable calcification rates (µg CaCO3 d−1).

Table 2Water chemistry in the bicarbonate ion manipulation experiment: mean parameters for experimental treatments over the course of the experiment. Columns from left to right present target and experimental salinity, pHtotal, partial pressure of carbon dioxide, total alkalinity, total dissolved inorganic carbon, target and actual bicarbonate ion concentration, carbonate ion concentration, carbon dioxide concentration, calcium ion concentration, aragonite saturation state, substrate–inhibitor ratio, and extended substrate–inhibitor ratio. All carbonate chemistry parameters were calculated using pHNBS and CT as inputs, and [Ca2+] was calculated as described in Sect. 2.1. Mean values are presented for each experimental treatment (n=4) ± standard deviation with measurements taken immediately before and after water changes.

Download Print Version | Download XLSX

Table 3Water chemistry in the calcium ion manipulation experiment: mean parameters for experimental treatments over the course of the experiment. Columns from left to right present target and experimental salinity, pHtotal, partial pressure of carbon dioxide, total alkalinity, total dissolved inorganic carbon, bicarbonate ion concentration, carbonate ion concentration, carbon dioxide concentration, target and experimental calcium ion concentrations, aragonite saturation state, substrate–inhibitor ratio, and extended substrate–inhibitor ratio. All carbonate chemistry parameters were calculated using pHNBS and CT as inputs. Mean values are presented for each experimental treatment (n=4) ± standard deviation with measurements taken immediately before and after water changes.

Download Print Version | Download XLSX

Figure 5Individual calcification rates (µg CaCO3 d−1) in both laboratory experiments. Calcification rates across [HCO3-] in the bicarbonate experiment are shown with negative exponential decay models for each salinity treatment (a). Calcification rates across [Ca2+] in the calcium experiment are fitted with linear trend lines for each salinity treatment (b). The model parameter Cmax (µg CaCO3 d−1) for each salinity in the bicarbonate experiment is depicted ± 95 % confidence interval (c).


Figure 6Mean calcification rates in each salinity treatment from both laboratory experiments ([HCO3-] – circles and [Ca2+] – triangles) across calculated Ωaragonite (a) and ESIR ([Ca2+][HCO3-] / [H+]) (b). Non-linear negative exponential decay models are graphically presented across both laboratory experiments, and vertical dotted lines represent the theoretical saturation thresholds, below which calcification rates are severely impeded (Ωaragonite=1, ESIR = 0.7).


Figure 7Field calcification over 15 months at the three monitoring sites in the Baltic Sea (Kiel, Ahrenshoop and Usedom). Shell mass was calculated from shell SL–CaCO3 mass relationships for each population (Fig. S3).


4 Discussion

4.1 Calcification in the field

The aim of this study was to identify the primary abiotic drivers responsible for decreasing calcification rates in mussel reefs along the Baltic Sea salinity gradient. Our results revealed both ESIR and Ωaragonite values were almost permanently below critical thresholds and saturation levels, respectively, in Ahrenshoop and Usedom and extended periods of undersaturation were observed in Kiel Fjord. Yet despite these conditions, net positive calcification rates were recorded at all sites with mussels in Kiel Fjord exhibiting calcification rates 1–2 orders of magnitude higher than the two low-salinity sites. This is interesting, as extremely low and variable pH (Fig. 3b) and other carbonate chemistry parameters have been well documented in Kiel Fjord with extended periods of Ωaragonite< 1 (undersaturation), particularly in early autumn during seasonal upwelling of hypercapnic water (Thomsen et al., 2010; Melzner et al., 2013; Saderne et al., 2013). These drastically higher rates of calcification in Kiel Fjord are unlikely to result solely from high ESIR and Ωaragonite values, as the laboratory experiments in this study clearly reveal minor differences in calcification rates above and below ESIR saturation (Fig. 6), compared to field calcification rates. It is more probable that the high levels of eutrophication in Kiel Fjord and consequently high phytoplankton food concentrations are responsible for the rapid calcification observed in this study (Melzner et al., 2011; Thomsen et al., 2013). Chl-a monitoring data in this study partially support this with mean values in Kiel significantly higher than those at Ahrenshoop. However, chl-a values in Usedom (low salinity) were comparable with those measured in Kiel despite drastically lower calcification rates in Usedom. Two-fold differences in particulate organic carbon (POC) have been documented between the inner and outer Kiel Fjord which cause higher food availability in the inner Fjord and thereby facilitate higher calcification rates (Thomsen et al., 2013). The growth monitoring site of this study was also located in the inner Fjord; however, chl-a measurements were taken in the outer Fjord and subsequently may not be fully representative of the Kiel mussel monitoring site. Values for chl-a and food availability to inner Kiel Fjord mussels may in fact be significantly higher than reported here, which likely explains the extremely high growth rates in Kiel compared to the other two monitoring sites. Conversely, the fact that calcification rates in Usedom were considerably lower than Kiel is likely due to limited [Ca2+] and sub-optimal ESIR values. The laboratory experiments presented here indicate that calcification rates decrease linearly when [Ca2+] < 4 mmol kg−1, and monitoring data revealed mean [Ca2+] at Ahrenshoop and Usedom, but not in Kiel, to be well below 4 mmol kg−1 (Fig. 2). Although pH is intrinsically linked to both Ω and ESIR, mean pHtotal values were not different between the three sites (Fig. 3a) despite significant differences between Ω and ESIR. [HCO3-] was  10 % lower in Ahrenshoop and Usedom compared to Kiel, whereas [Ca2+] was 22 % and 48 % lower in Ahrenshoop and Usedom than in Kiel, respectively. This suggests that low calcification rates at salinities  11 in the southwest Baltic result from low ESIR or Ω values primarily due to low seawater Ca2+ availability, rather than high [H+] and low [HCO3-] (Fig. S9).

4.2 Carbonate chemistry

This study aimed to investigate the mechanisms by which low salinity negatively impacts calcification in Baltic mussels. Calcification relies on the uptake of both calcium and inorganic carbon from seawater to precipitate CaCO3 crystals, and in the low saline southwest Baltic Sea, the availability of both substrates is significantly lower than in the North Sea (Kremling and Wilhelm, 1997; Müller et al., 2016). The laboratory experiments presented in this study demonstrate that HCO3- availability begins to limit calcification below  1000 µmol kg−1 (Fig. 5a), in line with similar experiments on Mytilus edulis juveniles (Thomsen et al., 2015). This corresponds to an AT of  1050 µmol kg−1 (at the experimental pH of 7.71). Such low [HCO3-] values are unlikely under fully marine conditions, due to high and stable AT in ocean waters of  2300 µmol kg−1 (Millero et al., 1998). However, estuarine and brackish water environments are often characterized by low AT due the linear relationship between salinity and AT and characteristically have a weaker pH buffering capacity and lower availability of inorganic carbon (Miller et al., 2009). Despite the low salinity, the Baltic Sea regionally exhibits relatively high AT due to high riverine AT from the southern drainage basin (Beldowski et al., 2010; Müller et al., 2016). Monitoring results from three coastal sites in the southwest Baltic revealed higher values of AT for any given salinity, when compared to the 2014 SAT relationship presented in Müller et al. (2016) (Fig. 4). As samples were taken in direct proximity to the coastline, this is likely due to the influence of local freshwater endmembers (either rivers or submarine groundwater discharge). Local freshwater endmembers in the southern drainage basin are known to have a higher AT than the mean volume-weighted Baltic Sea freshwater endmember, which is reduced by the influence of low AT rivers in Scandinavia (Müller et al., 2016). The proximity to the southern coastline and higher fraction of freshwater from the southern drainage basin likely explains the offset of the ATS relationships from our monitoring sites with respect to the relationships reported by Müller et al. (2016) for example. Dissolved organic carbon (DOC) species also constitute 1.5 %–2.3 % of AT in the southwest and central Baltic, and DOC from the Oder river discharge may contribute to the high coastal AT observed in the southwest Baltic in this study (Kuliński et al., 2014). The high AT observed across all sites in this study (1856–2071 µmol kg−1) suggests that HCO3- limitation at AT values < 1050 µmol kg−1 is likely not the primary reason for reduced calcification in Baltic mussels.

It is not possible to disentangle the effects of inorganic carbon availability from pH in the bicarbonate limitation experiment. This is because the methodologies applied in this study to manipulate [HCO3-] involved reducing AT in parallel and consequently lowering seawater pH as pCO2 remained unchanged (Table 2). Previous experiments utilizing methods which manipulate individual parameters of the carbonate system have suggested pH as a sole factor does not correlate strongly with molluscan calcification when compared to SIR, [CO32+] and Ωaragonite (Waldbusser et al., 2014; Thomsen et al., 2015). Reduced seawater pH through experimentally increasing pCO2 has been demonstrated to reduce calcification rates in bivalve larvae and coral due to a reduction in the proton gradient between ambient seawater and the calcifying space leading to reduced Ωaragonite at the site of calcification (Allison et al., 2014; Waldbusser et al., 2014; Ramesh et al., 2017). However, it is not possible to differentiate between the impacts of ambient seawater Ω, or alterations in carbonate chemistry at the site of calcification on net CaCO3 precipitation/dissolution in marine bivalves, since many marine organisms actively modify the carbonate chemistry at the site of calcification (Allison et al., 2014; Cyronak et al., 2016; Ramesh et al., 2017). Calcification is stimulated by the availability of calcification substrate (Ca2+ and HCO3-) and inhibited by H+; therefore the application of the SIR provides a more physiologically oriented predictor of calcification than Ω alone. However, the nature of the carbonate system results in the SIR co-correlating strongly with both [CO32+] and Ω at a constant salinity, temperature and pressure (Bach, 2015; Thomsen et al., 2015; Cyronak et al., 2016; Fassbender et al., 2016).

In our experiments SIR ([HCO3-] / [H+]) proved to be a poor predictor of calcification (Fig. S5) in comparison to Ωaragonite (Fig. 5a) when [Ca2+] was experimentally manipulated. Building on the findings from the HCO3- limitation experiment, this second experiment individually isolated the impacts of salinity and [Ca2+] revealing that salinity as a sole factor had no significant impact on calcification rates, whilst [Ca2+] correlated linearly with calcification across all [Ca2+]. This suggests Ca2+ to be the limiting factor for calcification at concentrations < 4 mmol kg−1 or corresponding to a salinity of  11 on the natural salinity gradient, in line with previous studies on Baltic mussel larvae and adults (Kossak, 2006; Sanders et al., 2018; Thomsen et al., 2018). Negative exponential decay models did not return significant parameter estimations for Cmax (maximum calcification rates) in the calcium experiment as in the bicarbonate experiment (Table S8). This indicates that threshold saturation values of [Ca2+] are above the maximum experimental levels tested here (ca. 3.8 mmol kg−1). A threshold of  4 mmol kg−1 for calcification has been observed for larvae and seems to be comparable for juveniles as well (Thomsen et al., 2018). Further experiments spanning a wider range of [Ca2+], e.g. 0.5–8 mmol kg−1, may reveal significant effects of salinity on calcification when [Ca2+] > 4 mmol kg−1. As a single factor, [Ca2+] was a poor predictor of calcification across both laboratory experiments (Fig. S4b). However, in this study, extending the SIR by including [Ca2+] as a substrate (ESIR) yielded a significant relationship with calcification, in line with Ωaragonite (Fig. 5). There was no statistically significant difference between both ESIR and Ωaragonite models as calcification predictors in this study, which may result from the low values of both in the Baltic Sea. Ωaragonite is particularly low in winter when low temperatures (Fig. S7) and high pCO2 from upwelling and water mixing increase Ksp and decrease [CO3+], respectively (Eq. 1). Salinity and temperature also impact aragonite saturation state (Ωaragonite) due to the salinity- and temperature-dependent nature of the solubility product of aragonite, Ksp (Eq. 1) and the CO2 dissociation constants. Conversely, temperature has almost no impact on ESIR when salinity and pressure are constant (see Bach, 2015). ESIR is also unimpacted by salinity per se, but rather [Ca2+] which correlates linearly with salinity in the Baltic Sea. Although future trends in the S–Ca2+ relationship in the southwest and central Baltic may deviate from current relationships (Sect. 4.5). Both ESIR and Ω generally correlate well in marine systems; however, this linear relationship is skewed in coastal zones where variability in carbonate chemistry parameters is high (Fassbender et al., 2016). This is highlighted in Fig. S8, where a given Ωaragonite of 1.0 can correspond to ESIR values of between 0.5–1.3 (well above and below the limiting ESIR threshold of 0.7). ESIR and Ωaragonite are certainly stronger predictors of calcification than calcification substrate availability alone in the southwest Baltic. However, the strong interdependency between both parameters makes it impossible to differentiate one from the other as a more powerful predictor of calcification in this study.

4.3 Salinity and calcification

Salinity is known to have a strong impact on bivalve growth and survival in the Baltic Sea, drastically reducing growth rates at salinities  11 (Kautsky et al., 1990; Kossak, 2006: Riisgård et al., 2014; Sanders et al., 2018). It has long been assumed that the underlying physiological mechanism of this stems from intracellular osmotic stress at low salinities and inefficient protein metabolism (Tedengren and Kautsky, 1986; Maar et al., 2015). However, recent work suggests that a reduced ability to biomineralize CaCO3 may be limiting growth rates at low salinities, rather than inefficient growth associated with osmotic stress (Riisgård et al., 2014; Sanders et al., 2018; Thomsen et al., 2018; Sillanpää et al., 2020). This is the first study to decouple salinity from both [HCO3-] and [Ca2+] and empirically test the impacts of each individual factor on bivalve calcification. In line with the findings of previous studies, both field monitoring and the bicarbonate laboratory experiment show that at salinities < 11, calcification rate is significantly reduced. However, when decoupled from [Ca2+] in the calcium experiment, salinity as a single factor did not significantly impact calcification. Interestingly, calcification rates were significantly reduced at a salinity of 6 in the bicarbonate experiment despite salinity having no impact on calcification in the calcium experiment. This could be a result of maximum [Ca2+] in the calcium experiment ( 3.8 mmol kg−1) being slightly below saturation levels ( 4 mmol kg−1). Thus, mild Ca2+ limitation may have masked the impacts of low salinity on calcification. The observed interaction between [Ca2+] and salinity (Table S6) on calcification supports this and, as discussed in Sect. 4.2, further experiments at saturated calcium concentrations (> 4 mmol kg−1) may reveal significant impacts of salinity on calcification in line with the bicarbonate experiment (Fig. 5a). The sampled experimental population at Ahrenshoop is also known to have high genetic diversity (Stuckas et al., 2017), and differences in the genetic composition between juvenile cohorts in both laboratory experiments in this study may explain the better performance at salinity 6 in the calcium experiment compared to the bicarbonate experiment. This coupled with slightly smaller animals in the bicarbonate experiment may be responsible for lower calcification rates at a salinity of 6 and higher mortality rates in the bicarbonate ion experiment. Future work should investigate how population-specific changes in genetics and physiological adaptation impact sensitivity of calcifying mussels to future desalination in the Baltic Sea.

The metabolic costs associated with osmotic stress at low salinities result from cellular mechanisms of volume regulation during salinity changes (Neufeld and Wright, 1996). Intracellular organic osmolytes (primarily free amino acids and quaternary ammonium compounds) are excreted from cells or mobilized from endogenous protein reserves during hypo- and hyper-salinity exposure, respectively, and associated with an energetic cost (Hawkins and Hilbish, 1992). Marine bivalves become iso-osmotic with surrounding seawater once long-term (> 2 weeks) acclimation to low salinity is complete and subsequently the energetic costs of protein excretion/mobilization should not persist (Willmer, 1978; Neufeld and Wright, 1996). Both laboratory experiments in this study lasted > 30 d, long enough to overcome the initial energetic costs of salinity acclimation. Despite both experiments having different durations (70 and 37 d for the bicarbonate and calcium experiments, respectively), it is unlikely that dissimilar exposure times would result in differential impacts of osmotic stress accumulation over time, as experimental durations were long enough to allow physiological acclimation to experimental salinities. Accordingly, the results here suggest that at salinities of 7–16, osmotic stress is not the primary cause of reduced growth in Baltic mussels.

4.4 Calcification physiology

Reduced calcification rates under sub-optimal carbonate chemistry conditions in this study may be explained in terms of calcification physiology. Adult bivalve molluscs are poor at regulating the ion composition of their extracellular fluid; thus low seawater [Ca2+] directly translates to low haemolymph [Ca2+] with the same concept applying to seawater pCO2, pH and [HCO3-]. However, extracellular carbonate chemistry is even less favourable for calcification than ambient seawater as production of metabolic CO2 maintains extracellular pCO2 higher (and pH lower) than ambient seawater to aid diffusive excretion of metabolic CO2 (Melzner et al., 2009; Heinemann et al., 2012). Sub-optimal carbonate chemistry in the extracellular fluid will therefore be in direct contact with the mantle epithelia (calcifying tissue) and likely impact energy-demanding ion transport processes involved in calcification. Increased costs of ion transport and maintenance of transmembrane ion gradients may be the likely driver of increasing mortality rates at lower pH values in the bicarbonate ion experiment (Fig. S10), compared to no mortality in the calcium ion experiment where pH values were similar to field values across all treatments (Tables 1 and 3). At a cellular level, low seawater Ca2+ availability may impact calcification by reducing rates of calcium transport across the mantle epithelia by plasma membrane Ca2+ ATPases (PMCA's) (Niggli et al., 1982; McConnaughey and Whelan, 1997). Although PMCA activity and gene expression has not been shown to be impacted by salinity in Mercenaria mercenaria and Crassoestrea gigas, respectively, actual rates of Ca2+ transport across the mantle epithelia have been shown to decrease at salinities of 14 compared to 28 in C. gigas (Ivanina et al., 2020; Sillanpää et al., 2020). This suggests that the reduced calcification rates observed in this study at low [Ca2+] may stem from a reduced supply of Ca2+ to the calcification site. Conversely, protons (H+) inhibit calcification in the calcifying space by lowering Ω and imposing constraints on mineral formation (Waldbusser et al., 2014; Cyronak et al., 2016). Organisms may increase the rate of H+ extrusion from the calcifying space (via the V-type H+-ATPase), but this requires more energy in the form of ATP, increasing calcification costs and potentially resulting in less energy available for other processes, e.g. protein synthesis (Waldbusser et al., 2013; Tresguerres, 2016). Higher energetic costs of calcification have been documented at low salinities (Sanders et al., 2018), and this may be the underlying mechanism explaining the reduced calcification rates at low salinities and [Ca2+]. Increased costs of calcification may also arise from changes in the organic content of the shell. The energetic cost of shell organic matrix protein synthesis is  20 times higher than the costs of CaCO3 precipitation (Palmer, 1992), and subsequently changes in the proportion of shell organic would have major impacts on the metabolic costs of shell production. Shell organic content was not quantified in this study; however, previous studies on Baltic Mytilus adults have observed higher organic content of shells in mussels living at low salinity (Telesca et al., 2019). Indeed, if this pattern is also true for juvenile mussels then this may be responsible for increased costs of calcification at low salinities. It remains elusive how the organic content of shells may be affected by long-term changes in seawater carbonate chemistry, and future work on calcification in the oceans should focus on how changes in shell organic content may influence the mechanisms and energetic costs of calcification.

4.5 Future environmental trends

Desalination of the southwest and central Baltic Proper is projected over the next century due to increased freshwater supply from river discharge resulting from changes in precipitation patterns (Meier, 2006; Gräwe et al., 2013). Depending on the composition of the additional freshwater input, this desalination might also suggest a reduction in AT and [Ca2+]. However, AT has in fact increased in the central Baltic Sea by  3.4 µmol kg−1 yr−1 over the past two decades, likely due to changes in precipitation and continental weathering (Müller et al., 2016), which has also increased [Ca2+] and Ca–AT ratios over a similar period (Kremling and Wilhelm, 1997). Although future changes in AT and [Ca2+] may be unclear at present, the magnitude of predicted desalination, coupled with a projected reduction in eutrophication and increasing atmospheric pCO2, points towards a trend of decreasing pH in the central Baltic Sea over the next century (Gustafsson et al., 2019; Gustafsson and Gustafsson, 2020). Future acidification combined with a potential mismatch in future S–Ca2+ relationships highlight the importance of utilizing ESIR as well as Ω in predicting vulnerability of marine calcifiers in the southwest and central Baltic Sea. In the absence of adaptation, projected environmental change may result in a westward distribution shift of marine mussels towards the Kattegat, and/or replacement of calcifying mussels in the southwest and central Baltic by brackish or freshwater calcifying bivalves such as the invasive euryhaline freshwater bivalve, Dreissena spp. These changes may have potential knock-on effects for benthic biodiversity and ecosystem function in the southwest and central Baltic Sea (Vuorinen et al., 2015).

5 Conclusions

The environmental gradients in the Baltic Sea provide an excellent system to investigate the impacts of seawater carbonate chemistry on marine calcifying organisms. The results presented here show that extremely slow calcification rates in the Baltic Sea arise from extended periods of Ωaragonite< 1 and low ESIR undersaturation. These sub-optimal conditions for calcification arise primarily from low seawater [Ca2+] rather than low [HCO3-] and low pH. Although high food availability and high AT near the coast may alleviate these negative impacts at salinities > 11, this is likely not possible when [Ca2+] is ultimately limiting calcification < 4 mm kg−1. Predicted acidification and desalination in the southwest Baltic is likely to impact Baltic Mytilus calcification and distributions; however, this will depend on future [Ca2+] trends. Upcoming work investigating the impacts of climate change on Baltic Sea calcifying organisms should (a) focus on better understanding future changes in [Ca2+] in the southwest and central Baltic Proper and (b) further investigate the mechanisms and costs associated with cellular Ca2+ transport in bivalves and assess their potential to adapt their calcification machinery to a stressful future Baltic Sea environment.

Data availability

All raw data are available via the Pangaea database: (Sanders et al., 2020).


The supplement related to this article is available online at:

Author contributions

TS, FM and JT designed the laboratory experiments conducted by TS and JT. TS, FM, GR and JDM implemented the environmental monitoring, and TS collected and analysed environmental data. JDM and GR analysed carbonate chemistry samples. TS wrote the paper with all authors providing significant contributions throughout.

Competing interests

The authors declare that they have no conflict of interest.


The authors thank Stefan Otto at the IOW for carbonate chemistry analysis of samples and Thomas Stegmann at the Christian Albrechts Universität zu Kiel for Ca2+ measurements. The authors also thank the Landesamt für Landwirtschaft, Umwelt und Ländliche Räume (LLUR, Kiel) and the Landesamt für Umwelt, Naturschutz und Geologie (LUNG, Güstrow) for chl-a monitoring data. Additionally, the authors would like to thank Luca Telesca (Cambridge University) for assisting with the analysis of environmental data.

Financial support

This research was supported by the Marie Curie ITN network “CACHE” (Calcium in a Changing Environment) European Union Seventh Framework Programme under grant agreement no. 605051. The position of Jens Daniel Müller was funded by BONUS, the joint Baltic Sea research and development programme (Art 185), and jointly by the European Union's Seventh Programme for research, technological development and demonstration and the German Federal Ministry of Education and Research through grant no. 03F0689A (BONUS PINBAL) and grant no. 03F0773A (BONUS INTEGRAL).

The article processing charges for this open-access publication were covered by a Research Centre of the Helmholtz Association.

Review statement

This paper was edited by Jean-Pierre Gattuso and reviewed by Laura Ramajo and one anonymous referee.


Allison, N., Cohen, I., Finch, A. A., Erez, J., Tudhope, A. W., and Edinburgh Ion Microprobe Facility: Corals concentrate dissolved inorganic carbon to facilitate calcification, Nat. Commun., 5, 5741,, 2014. 

Attard, K. M., Rodil, I. F., Berg, P., Mogg, A. O. M., Westerbom, M., Norkko, A., and Gludd, R. N.: Metabolism of a subtidal rocky mussel reef in a high-temperate setting: pathways of organic C flow, Mar. Ecol.-Prog. Ser., 645, 41–54,, 2020. 

Bach, L. T.: Reconsidering the role of carbonate ion concentration in calcification by marine organisms, Biogeosciences, 12, 4939–4951,, 2015. 

Beldowski, J., Löffler, A., Schneider, B., and Joensuu, L.: Distribution and biogeochemical control of the total CO2 and total alkalinity in the Baltic Sea, J. Marine Syst., 81, 252–259,, 2010. 

Cyronak, T., Schulz, K. G., and Jokiel, P. L.: The Omega myth: what really drives lower calcification rates in an acidifying ocean, ICES J. Mar. Sci., 73, 558–562,, 2016. 

Dickson, A. G.: Standard potential of the reaction – AgCl(s) + 12H2= Ag(s) + HCl(aq) and the standard acidity constant of the ion HSO4- in synthetic sea-water from 273.15-K to 318.15-K, J. Chem. Thermodyn., 22, 113–127,, 1990. 

Dickson, A. G., Afgan, J. D., and Anderson, G. C.: Reference materials for oceanic CO2 analysis: a method for the certification of total alkalinity, Mar. Chem., 80, 185–197,, 2003. 

Dickson, A. G., Sabine, C.vL., and Christian, J. R. (Eds.): Guide to Best Practices for Ocean CO2 Measurements, PICES Special Publication 3, IOCCP Report No. 8, 191 pp., 2007. 

Elmgren, R. and Hill, C. (Eds.): Ecosystem function at low biodiversity – the Baltic example, in: Marine Biodiversity: Patterns and Processes, edited by: Ormond, R. F. G., Gage, J. D., and Angel, M. V., Cambridge University Press, Cambridge, UK, 319–336,, 1997. 

EU, Copernicus Marine Service: Copernicus Marine Environment Monitoring Service – CMEMS, available at:, last access: 7 February 2018. 

Fassbender, A. J., Sabine, C. L., and Feifel, K. M.: Consideration of coastal carbonate chemistry in understanding biological calcification, Geophys. Res. Lett., 43, 4467–4476,, 2016. 

Gräwe, U., Friedland, R., and Burchard, H.: The future of the western Baltic Sea: two possible scenarios, Ocean Dynam., 63, 901–921,, 2013. 

Gustafsson, E. and Gustafsson, B. G.: Future acidification of the Baltic Sea – A sensitivity study, J. Marine Syst., 211, 103397,, 2020. 

Gustafsson, E., Hagens, M., Sun, X., Reed, D. C., Humborg, C., Slomp, C. P., and Gustafsson, B. G.: Sedimentary alkalinity generation and long-term alkalinity development in the Baltic Sea, Biogeosciences, 16, 437–456,, 2019. 

Hammer, K., Schneider, B., Kuliński, K., and Schulz-Bull, D. E.: Precision and accuracy of spectrophotometric pH measurements at environmental conditions in the Baltic Sea, Estuar. Coast. Shelf S., 146, 24–32,, 2014. 

Hawkins, A. J. S. and Hilbish, T. J.: The cost of cell volume regulation: protein metabolism during hyperosmotic adjustment, J. Mar. Biol. Assoc. UK, 72, 569–578,, 1992. 

Heckwolf, M. J., Peterson, A., Jänes, H., Horne, P., Künne, J., Liversage, K., Sajeva, M., Reusch, T. B. H., and Kotta, J.: From ecosystems to socio-economic benefits: A systematic review of coastal ecosystem services in the Baltic Sea, Sci. Total Environ., 755, 142565,, 2021. 

Heinemann, A., Fietzke, J., Melzner, F., Böhm, F., Thomsen, J., Garbe-Schönberg, D., and Eisenhauer, A.: Conditions of Mytilus edulis extracellular body fluids and shell composition in a pH-treatment experiment: Acid-base status, trace elements and δ11B, Geochem. Geophy. Geosy., 13, Q01005,, 2012. 

Hiebenthal, C., Fietzek, P., Müller, J. D., Otto, S., Rehder, G., Paulsen, M., Stuhr, A., Clemmesen, C., and Melzner, F.: Kiel Fjord carbonate chemistry data between 2015 (January) and 2016 (January), GEOMAR – Helmholtz Centre for Ocean Research Kiel, PANGEA,, 2017. 

Hu, H. Y., Yan, J.-J., Petersen, I., Himmerkus, N., Bleich, M., and Stumpp, M.: A SLC4 family bicarbonate transporter is critical for intracellular pH regulation and biomineralization in sea urchin embryos, eLife, 7, e36600,, 2018. 

Ivanina, A. V., Jarrett, A., Bell, T., Rimkevicius, T., Beniash, E., and Sokolova, I. M.: Effects of seawater salinity and pH on cellular metabolism and enzyme activities in biomineralizing tissues of marine bivalves, Comp. Biochem. Phys. A, 248, 110748,, 2020. 

Jokiel, P. L.: Coral reef calcification: carbonate, bicarbonate and proton flux under conditions of increasing ocean acidification, P. Roy. Soc. B, 280, 20130031,, 2013. 

Kautsky, N., Johannesson, K., and Tedengren, M.: Genotypic and phenotypic differences between Baltic and North Sea populations of Mytilus edulis evaluated through reciprocal transplantations. I. Growth and morphology, Mar. Ecol.-Prog. Ser., 59, 203–210,, 1990. 

Kester, D. R., Duedall, I. W., Connors, D. N., and Pytkowicz, R. M.: Preparation of artificial seawater, Limnol. Oceanogr., 12, 176–179,, 1967. 

Koivisto, M. E. and Westerbom, M.: Habitat structure and complexity as determinants of biodiversity in blue mussel beds on sublittoral rocky shores, Mar. Biol., 157, 1463–1474,, 2010. 

Kossak, U.: How climate change translates into ecological change: Impacts of warming and desalination on prey properties and predator-prey interactions in the Baltic Sea, PhD Thesis, Mathematics and Natural Sciences faculty of Christian Albrechts University, Kiel, 2006. 

Kotta, J., Futter, M., Kaasik, A., Liversage, K., Rätsep, M., Barboza, F. R., Bergström, L., Bergström, P., Bobsien, I., Díaz, E., Herkül, K., Jonsson, P. R., Korpinen, S., Kraufvelin, P., Krost, P., Lindahl, O., Lindegarth, M., Lyngsgaard, M. M., Mühl, M., Sandman, A. N., Orav-Kotta, H., Orlova, M., Skov, H., Rissanen, J., Šiaulys, A., Vidakovic, A., and Virtanen, E.: Cleaning up seas using blue growth initiatives: Mussel farming for eutrophication control in the Baltic Sea, Sci. Total Environ., 709, 136144,, 2020. 

Kremling, K. and Wilhelm, G.: Recent increase of the calcium concentrations in Baltic Sea waters, Mar. Pollut. Bull., 34, 763–767,, 1997. 

Kuliński, K., Schneider, B., Hammer, K., Machulik, U., and Schulz-Bull, D.: The influence of dissolved organic matter on the acid-base system of the Baltic Sea, J. Marine Syst., 132, 106–115,, 2014. 

Maar, M., Saurel, C., Landes, A., Dolmer, P., and Petersen, J. K.: Growth potential of blue mussels (M. edulis) exposed to different salinities evaluated by a Dynamic Energy Budget model, J. Marine Syst., 148, 48–55,, 2015. 

McConnaughey, T. A. and Whelan, J. F.: Calcification generates protons for nutrient and bicarbonate uptake, Earth-Sci. Rev., 42, 95–117,, 1997. 

Meier, H. E. M.: Baltic Sea climate in the late twenty-first century: a dynamical downscaling approach using two global models and two emission scenarios, Clim. Dynam., 27, 39–68,, 2006. 

Melzner, F., Gutowska, M. A., Langenbuch, M., Dupont, S., Lucassen, M., Thorndyke, M. C., Bleich, M., and Pörtner, H.-O.: Physiological basis for high CO2 tolerance in marine ectothermic animals: pre-adaptation through lifestyle and ontogeny?, Biogeosciences, 6, 2313–2331,, 2009. 

Melzner, F., Stange, P., Trübenbach, K., Thomsen, J., Casties, I., Panknin, U., Gorb, S. N., and Gutowska, M.: Food supply and seawater pCO2 impact calcification and internal shell dissolution in the blue mussel Mytilus edulis, PLoS ONE, 6, e24223,, 2011. 

Melzner, F., Thomsen, J., Koeve, W., Oschlies, W., Gutowska, M. A., Bange, H. W., Hansen, H. P., and Körtzinger, A.: Future ocean acidification will be amplified by hypoxia in coastal habitats, Mar. Biol., 160, 1875–1888,, 2013. 

Miller, A. W., Reynolds, A. C., Sobrino, C., and Riedel, G. F.: Shellfish face uncertain future in high CO2 world: Influence of acidification on oyster larvae calcification and growth in estuaries, PLoS ONE, 4, E5661,, 2009. 

Millero, F.: The conductivity-density-salinity-chlorinity relationships for estuarine water, Limnol. Oceanogr., 29, 1317–1321,, 1984. 

Millero, F. J.: Carbonate constants for estuarine waters, Mar. Freshwater Res., 61, 139–142,, 2010. 

Millero, F. J., Lee, K., and Roche, M.: Distribution of alkalinity in the surface waters of the major oceans, Mar. Chem., 60, 111–130,, 1998. 

Mohrholz, V., Naumann, M., Nausch, G., Krüger, S., and Gräwe, U.: Fresh oxygen for the Baltic Sea – an exceptional saline inflow after a decade of stagnation, J. Marine Syst. 148, 152–166,, 2015. 

Müller, J. D. and Rehder, G.: Metrology of pH Measurements in Brackish Waters – Part 2: Experimental Characterization of Purified meta-Cresol Purple for Spectrophotometric pHT Measurements, Front. Mar. Sci., 5, 177,, 2018. 

Müller, J. D., Schneider, B., and Rehder, G.: Long-term alkalinity trends in the Baltic Sea and their implications for CO2-induced acidification, Limnol. Oceanogr., 61, 1984–2002,, 2016. 

Neufeld, D. S. and Wright, S. H.: Response of cell volume in Mytilus gill to acute salinity change, J. Exp. Biol., 199, 473–484, 1996. 

Neumann, T.: Climate-change effects on the Baltic Sea ecosystem: A model study, J. Marine Syst., 81, 213–224,, 2010. 

Niggli, V., Sigel, E., and Carafoli, E.: The purified Ca2+ pump of human erythrocyte membranes catalyzes an electroneutral Ca2+-H+ exchange in reconstituted liposomal systems, Cell Calcium, 257, 2350–2356,, 1982. 

Norling, P. and Kautsky, N.: Patches of the mussel Mytilus sp. are islands of high biodiversity in subtidal sediment habitats in the Baltic Sea, Aquat. Biol., 4, 75–87,, 2008. 

Palmer, A. R.: Calcification in marine molluscs: How costly is it?, P. Natl. Acad. Sci. USA, 89, 1379–1382,, 1992. 

R Core Team: R: A language and environment for statistical computing, R Foundation for Statistical Computing, Vienna, Austria, (last access: 12 March 2021), 2020. 

Ramesh, K., Hu, M. Y., Thomsen, J., Bleich, M., and Melzner, F.: Mussel larvae modify calcifying fluid carbonate chemistry to promote calcification, Nat. Commun., 8, 1709,, 2017. 

Ries, J. B., Ghazaleh, M. N., Connolly, B., Westfiel, I., and Castillo, K. D.: Impacts of seawater saturation state (ΩA = 0.4–4.6) and temperature (10, 25 C) on the dissolution kinetics of whole-shell biogenic carbonates, Geochim. Cosmochim. Ac., 192, 318–337,, 2016. 

Riisgård, H. U., Pleissner, D., Lundgreen, K., and Larse, P. S.: Growth of mussels Mytlus edulis, at algal (Rhodomonas salina) concentrations below and above saturation levels for reduced filtration rate, Mar. Biol. Res., 9, 1005–1017,, 2013. 

Riisgård, H. U., Larsen, P. S., Turja, R., and Lundgreen, K.: Dwarfism of blue mussels in the low saline Baltic Sea – growth to the lower salinity limit, Mar. Ecol.-Prog. Ser., 517, 181–192,, 2014. 

Roleda, M. Y., Boyd, P. W., and Hurd, C. L.: Before ocean acidification: calcifier chemistry lessons (1), J. Phycol., 48, 840–843,, 2012. 

Saderne, V., Fietzek, P., and Herman, P. M. J.: Extreme variations of pCO2 and pH in a macrophyte meadow of the Baltic Sea in summer: evidence of the effect of photosynthesis and local upwelling, PLoS ONE, 8, e62689,, 2013. 

Sanders, T., Schmittmann, L., Nascimento-Schulze, J., and Melzner, F.: High calcification costs limit mussel growth at low salinity, Front. Mar. Sci., 5, 352,, 2018. 

Sanders, T., Thomsen, J., Müller, J. D., Rehder, G., and Melzner, F.: Mussel calcification in the Baltic Sea under natural and manipulated carbonate chemistry parameters, PANGAEA,, 2020. 

Sillanpää, J. K., Cardoso, J. C. R., Félix, R. C., Anjos, L., Power, D. M., and Sundell, K.: Dilution of seawater affects the Ca2+ transport in the outer mantle epithelium of Crassostrea gigas, Front. Physiol., 11, 1,, 2020. 

Stuckas, H., Knöbel, L., Schade, H., Breusing, C., Hinrichsen, H. H., Bartel, M., Langguth, C., and Melzner, F.: Combining hydrodynamic modelling with genetics: can passive larval drift shape the genetic structure of Baltic Mytilus populations?, Mol. Ecol., 26, 2765–2782,, 2017. 

Tedengren, M. and Kautsky, N.: Comparative study of the physiology and its probable effect on size in blue mussels (Mytilus edulis L.) from the North Sea and the Northern Baltic Proper, Ophelia, 25, 147–155,, 1986. 

Telesca, L., Peck, L. S., Sanders, T., Thyrring, J., Sejr, M. K., and Harper, E. M.: Biomineralization plasticity and environmental heterogeneity predict geographical resilience patterns of foundation species to future change, Glob. Change Biol., 25, 4179–4193,, 2019. 

Thomas, H. and Schneider, B.: The seasonal cycle of carbon dioxide in Baltic Sea surface waters, J. Marine Syst., 22, 53–67,, 1999. 

Thomsen, J., Gutowska, M. A., Saphörster, J., Heinemann, A., Trübenbach, K., Fietzke, J., Hiebenthal, C., Eisenhauer, A., Körtzinger, A., Wahl, M., and Melzner, F.: Calcifying invertebrates succeed in a naturally CO2-rich coastal habitat but are threatened by high levels of future acidification, Biogeosciences, 7, 3879–3891,, 2010. 

Thomsen, J., Casties, I., Pansch, C., Körtzinger, A., and Melzner, F.: Food availability outweighs ocean acidification effects in juvenile Mytilus edulis: laboratory and field experiments, Glob. Change Biol., 19, 1017–1027,, 2013. 

Thomsen, J., Haynert, K., Wegner, K. M., and Melzner, F.: Impact of seawater carbonate chemistry on the calcification of marine bivalves, Biogeosciences, 12, 4209–4220,, 2015. 

Thomsen, J., Ramesh, K., Sanders, T., Bleich, M., and Melzner, F.: Calcification in a marginal sea – influence of seawater [Ca2+] and carbonate chemistry on bivalve shell formation, Biogeosciences, 15, 1469–1482,, 2018. 

Tresguerres, M.: Novel and potential physiological roles of vacuolar-type H+-ATPase in marine organisms, J. Exp. Biol., 219, 2088–2097,, 2016. 

Tyrrell, T., Schneider, B., Charalampopoulou, A., and Riebesell, U.: Coccolithophores and calcite saturation state in the Baltic and Black Seas, Biogeosciences, 5, 485–494,, 2008. 

Vuorinen, I., Antsulevich, A. E., and Maximovich, N. V.: Spatial distribution and growth of the common mussel Mytilus edulis L. in the archipelago of SW-Finland, northern Baltic Sea, Boreal Environ. Res., 7, 41–52, 2002. 

Vuorinen, I., Hänninen, J., Rajasilta, M., Laine, P., Eklund, J., Montesine-Pouzols, F., Corona, F., Junker, K., Meier, H. E. M., and Dippner, J. W.: Scenario simulations of future salinity and ecological consequences in the Baltic Sea and adjacent North Sea areas – implications for environmental monitoring, Ecol. Indic., 50, 196–205,, 2015. 

Wahl, M., Schneider Covachã, S., Saderne, V., Hiebenthal, C., Müller, J. D., Pansch, C., and Sawall, Y.: Macroalgae may mitigate ocean acidification effects on mussel calcification by increasing pH and its fluctuations, Limnol. Oceanogr., 63, 3–21,, 2017. 

Waldbusser, G. G., Brunner, E. L., Haley, B. A., Hales, B., Langdon, C. J., and Prahl, F. G.: A developmental and energetic basis linking larval oyster shell formation to acidification sensitivity, Geophys. Res. Lett., 40, 2171–2176,, 2013. 

Waldbusser, G. G., Hales, B., Langdon, C. J., Haley, B. A., Schrader, P., Brunner, E. L., Gray, M. W., Miller, C. A., and Gimenez, I.: Saturation-state sensitivity of marine bivalve larvae to ocean acidification, Nat. Clim. Change, 5, 273–280,, 2014. 

Wasmund, N., Tuimala, J., Suikkanen, S., Vandepitte, L., and Kraberg, A.: Long-term trends in phytoplankton composition in the western and central Baltic Sea, J. Marine Syst., 87, 145–159,, 2011. 

Westerbom, M., Kilpi, M., and Mustonen, O.: Blue mussels, Mytilus edulis, at the edge of the range: population structure, growth and biomass along a salinity gradient in the north-eastern Baltic Sea, Mar. Biol., 140, 991–999,, 2002. 

Westerbom, M., Mustonen, O., Jaatinen, K., Kilpi, M., and Norkko, A.: Population dynamics at the range margin: Implications of climate change on sublittoral blue mussels (Mytilus trossulus), Front. Mar. Sci., 6, 292,, 2019. 

Willmer, P. G.: Volume regulation and solute balance in the nervous tissue of an osmoconforming bivalve (Mytilus edulis), J. Exp. Biol., 77, 157–179, 1978. 

Zeebe, R. E. and Wolf-Gladrow, D. A.: CO2 in seawater: Equilibrium, kinetics, isotopes, Amsterdam, Elsevier Oceanography Series, 2001. 

Zoccola, D., Ganot, P., Bertucci, A., Caminiti-Segonds, N., Techer, N., Voolstra, C. R., Aranda, M., Tambutté, É., Allemand, D., Casey, J. R., and Tambutté, S.: Bicarbonate transporters in corals point towards a key step in the evolution of cnidarian calcification, Sci. Rep.-UK, 5, 9983,, 2015. 

Short summary
The Baltic Sea is expected to experience a rapid drop in salinity and increases in acidity and warming in the next century. Calcifying mussels dominate Baltic Sea seafloor ecosystems yet are sensitive to changes in seawater chemistry. We combine laboratory experiments and a field study and show that a lack of calcium causes extremely slow growth rates in mussels at low salinities. Subsequently, climate change in the Baltic may have drastic ramifications for Baltic seafloor ecosystems.
Final-revised paper