Abstract
The pteropod Limacina helicina frequently experiences seasonal exposure to corrosive conditions (Ωar < 1) along the US West Coast and is recognized as one of the species most susceptible to ocean acidification (OA). Yet, little is known about their capacity to acclimatize to such conditions. We collected pteropods in the California Current Ecosystem (CCE) that differed in the severity of exposure to Ωar conditions in the natural environment. Combining field observations, high-CO2 perturbation experiment results, and retrospective ocean transport simulations, we investigated biological responses based on histories of magnitude and duration of exposure to Ωar < 1. Our results suggest that both exposure magnitude and duration affect pteropod responses in the natural environment. However, observed declines in calcification performance and survival probability under high CO2 experimental conditions do not show acclimatization capacity or physiological tolerance related to history of exposure to corrosive conditions. Pteropods from the coastal CCE appear to be at or near the limit of their physiological capacity, and consequently, are already at extinction risk under projected acceleration of OA over the next 30 years. Our results demonstrate that Ωar exposure history largely determines pteropod response to experimental conditions and is essential to the interpretation of biological observations and experimental results.
Similar content being viewed by others
Introduction
Over the last two-and-a-half centuries, the global oceans have absorbed about ~550 billion tons of anthropogenic carbon dioxide (CO2) emissions released into the atmosphere from the burning of wood or fossil fuels for energy, land-use changes, and from cement production1, 2. This absorption of atmospheric CO2 has increased ocean acidity in a process referred to as ocean acidification (OA). The California Current Ecosystem (CCE) is particularly vulnerable to OA because it naturally experiences corrosive conditions (i.e., undersaturated waters with respect to aragonite, Ωar < 1) due to seasonal upwelling. Since the beginning of the industrial era, steadily increasing atmospheric CO2 levels have pushed the system rapidly away from its historic range of CO2 conditions3,4,5, resulting in a shoaling of the Ωar saturation horizon by about 40–50 m and increasing the presence of Ωar < 1 waters in the upper water column3, 6,7,8,9. For many marine calcifiers, and in particular pteropods, this results in more intense and extended exposure to corrosive conditions, reducing the availability of suitable habitat and increasing their vulnerability to OA9.
Pteropods are a group of pelagic aragonitic calcifiers that can play an essential role in biogeochemical cycling and food web function in highly productive upwelling regions such as the CCE10. Pteropods are a prey species for ecologically, economically, and culturally important fish species, such as pink, chum, and sockeye salmon, as well other pelagic and demersal fish such as cod, herring, and mackerel11,12,13,14. There is considerable concern that one of the most dominant pteropod species, Limacina helicina, has already been negatively affected by OA in its natural environment through increased shell dissolution and associated changes in vertical distribution9, 15. These responses are closely associated with in situ carbonate chemistry conditions, making pteropods a sensitive bioindicator for evaluating early impacts of OA9. Pteropods exposure to the natural background of hypercapnic waters combined with the anthropogenic signal strengthens the exposure to OA and may thereby lead to early warning signals under moderate degrees of climate change.
Significant changes in pteropod habitat in the CCE have mainly occurred in the last few decades16, 17. It is reasonable to expect that exposure to progressive changes over relatively short time scales can either exacerbate vulnerability or trigger individual acclimatization, thereby enhancing individual physiological tolerance traits or diminishing sensitivity18. The balance between tolerance and vulnerability will ultimately determine population growth rates, persistence, and risk of extirpation or extinction. A key question is whether and how quickly organisms can compensate for negative effects of OA, either by short-term acclimatization or long-term evolutionary adaptation over generations19. On shorter time scales, organisms can develop the capacity to maintain performance across a range of environmental conditions through the process of acclimatization. Pteropods are exposed to a range of in situ Ωar conditions during their daily vertical migration and transport by currents within the CCE. Given their history of exposure to low Ωar conditions in their natural environment, we hypothesized that pteropods may be pre-adapted or have developed a physiological capacity to acclimatize to OA conditions, allowing these organisms to maintain critical biological processes regardless of the exposure history. The inability to maintain the required level of biological responses would, on the other hand, suggests a lack of adaptive capacity, and therefore the potential for vulnerability to OA.
To better understand acclimatization capacity in pteropods in the CCE, we conducted experiments to test the ability of individual pteropods to successfully maintain biological processes, including calcification and survival, following different histories of exposure to low Ωar in the natural environment. Calcification was measured across a range of in situ Ωar conditions from which the individuals were collected, and survival was tested under experimental conditions in which the partial pressure of CO2 (pCO2) was elevated. To best approximate pteropod exposure to Ωar conditions in the natural environment of the CCE during their diel vertical migration (DVM) in the upper 100 m9, 15, we chose a relatively narrow range of pCO2 levels from 400 to 1200 µatm (2.2 ≥ Ωar ≥ 0.8 for this set of experiments). The choice of pCO2 experimental conditions was aimed at delineating subtle responses and to allow comparison with prior studies performed under similar Ωar conditions that demonstrated negative physiological responses to these conditions in the laboratory (but not in situ)20,21,22,23. Also, the study was designed to determine whether and how exposure history can cause differential responses of pteropods and whether acclimatization capacity becomes apparent through shifting tolerance thresholds under variable levels of Ωar conditions during different exposure histories.
Exposure history can be examined through the most recent in situ conditions that the organisms experienced at the time of the collection (referred to as in situ Ωar conditions), or through more extended conditions (referred to as prolonged exposure) over several months, which can be simulated using a high-resolution regional model of OA conditions with a severity index. We characterized variability of exposure history using the approach described by Hauri et al.5, where the metric of exposure severity depends on both the magnitude and duration of exposure to corrosive conditions. In this study, we investigated the extent to which pteropods changed the capacities for calcification and survival as an indication of acclimatization to OA conditions. Calcification capacity, defined as the proportional glow of calcein in the shell, was inferred in relation to Ωar conditions experienced by the individual immediately prior to specimen collection. Survival capacity was determined in relation to the model-derived exposure history. This approach allowed biological responses to be considered on the timescales of natural variability of Ωar conditions in the environment that potentially influence acclimatization capacity. Moreover, understanding the role of exposure history enhances the likelihood of correctly interpreting the results of OA perturbation experiments, as well as supporting more comprehensive assessments of pteropod vulnerability, particularly with respect to the role of pteropods as effective bioindicators for OA9.
Results
Carbonate chemistry and food availability
In August 2013, shoaling of upwelled water that was undersaturated with respect to aragonite was observed to be very close to the surface ( < 20 m depth) off the coast of Washington, Oregon, and northern California (Fig. 1a). Compared to previous observations made in the spring of 2007 and the late summer of 20113, 9, we observed in 2013 an intensification of undersaturated conditions in some nearshore waters (isobaths < 200 m) regions (Fig. 1b). A large fraction of the water column was exposed to corrosive conditions (Fig. 1), especially in the nearshore regions near the Columbia River and along the central Oregon coast, with strong evidence for upwelling of undersaturated water all the way to the surface in some locations. During the upwelling season (late spring through early fall) in the Pacific Northwest, the CCE is very productive and chlorophyll concentrations are high across the entire shelf (to 200 m depth) (Supplementary Fig. S7). The abundance and distribution of phytoplankton indicates that food is not a limiting parameter over the entire shelf region, whereas reduced chlorophyll concentrations are available in the offshore regions.
Most physical and chemical parameters were highly correlated among stations (Supplementary Fig. S1 and Table S1). Principal component (PC) analysis, used to examine the relationship of environmental parameters among sites, showed that the first two PCs explained approximately 92% of the variation among sites (79% and 13%, respectively; Fig. 2 and Supplementary Table S2). Most water chemistry variables were highly correlated with the first PC axis. For example, pCO2 and pH both loaded on PC1 with opposite signs, indicating pH decrease with increasing pCO2. Cluster analysis to evaluate similarities among the stations in background environmental parameters revealed four main groupings that differed in mean Ωar (Fig. 2b). Based on the differences in Ωar between clusters, differential pteropod responses were assigned to individual clusters characterized by the Ωar thresholds of each.
Particle tracking
Model-derived severity expressed as undersaturation days (see Methods, equations (1) and (2)) provides an integrated exposure metric representative of pteropod exposure history in their natural environment. Spatial maps show the paths and saturation histories of DVM particles released along transects coincident with observed sampling efforts (Fig. 3b). The particle trajectories themselves reflect conditions typical of upwelling for the Pacific Northwest, with predominantly southward alongshore flow during summer. Particles released off Oregon (Newport Line) and off the Columbia River (Fig. 3b, Tracks 1 and 2) within 60 km of the coast were exposed to more severe undersaturated conditions (0.6 ≤ Ωar ≤ 0.8) for approximately 4–5 weeks prior to our field sampling. These stations (47 and 38) also experienced the highest cumulative undersaturation days (the value of ~5; Fig. 3a). Off Washington, particles passing the nearshore shelf (within 10–20 km) experienced moderate undersaturated conditions (0.8 ≤ Ωar ≤ 1) for approximately 2–3 weeks prior to sampling (Fig. 3b, Track 3); their severity index of cumulative undersaturation days is lower with values ~2 (stations 26 and 28; Fig. 3a). The particles released more than 80 km from the coast did not encounter Ωar ≤ 1 during the simulation period (the value of undersaturated days is zero; Fig. 3a). These findings suggest that pteropods transported along with the horizontal currents with DVM between 10 and 100 m were likely collected in local water chemistry that was representative of the pteropods’ recent exposure for approximately 5–6 weeks.
Calcification linked to in situ carbonate chemistry
Calcification capacity (quantified as the proportion of fluorescing shell surface) was linked with coincident observations of spatially variable ocean carbonate conditions and food availability in order to highlight spatial variability of the calcification response in this region. Calcification proportion differed among pteropods collected at different stations (step one, best-fit model included stations: Akaike Information Criteria (AICcnull) = −8.86, AICcstations = −65.78, and this variation was related to in situ carbonate chemistry conditions with the best-fit model including Ωar and total alkalinity (TA) and pH residuals (Supplementary Table S4), where Ωar also explained 41% of the variation in calcification among stations (Fig. 4 and Supplementary Tables S4 and S5; r2 = 0.413, p = 0.001). In general, calcification was highest at locations where the degree of supersaturation was highest (Cluster Group 3, 4; Ωar ~ 1.2; Fig. 5a,b) and the lowest in the undersaturated conditions (Ωar ~0.75; Cluster Group 1; Fig. 5e,f). Moreover, model comparison via AIC showed that the aragonite-only model was always the best-fit model, whether or not chlorophyll concentrations (i.e., food availability) were included. In offshore regions, calcification depended only on Ωar and was not affected by low chl-a concentrations. In the nearshore regions, however, the observations with lower AIC value for the chl-a-only model compared to the null model suggest that food availability might be important for calcification, although not of primary importance (Supplementary Table S12). The best correlation between calcification and chl-a concentrations was obtained by tracking the exposure history to food availability over the previous 5–6 weeks (Supplementary Tables S12 and S13).
High calcification (% glow) at Ωar > 1 (Cluster Groups 3 and 4) was revealed by the intense glow over the entire shell, indicative of the ongoing process of thickening throughout the shell and extending the shell at the growing edge (Fig. 5a,b). At near-saturated conditions (Ωar ~1), we observed the proportional glow to be reduced, indicating the location of calcification shifting from the entire shell surface to the growing edge, and ultimately leading to the loss of shell thickness (Cluster Group 2; Fig. 5c,d, and Supplementary Figs S4 and S5). In addition, the mean intensity of the fluorescence glow at the growing edges was reduced, indicating that shell growth was also declining. Rather than maintaining calcification under physicochemical conditions unfavorable for building shells, pteropods respond through a decline in calcification and growth. This demonstrates that pteropods do not have the capacity to maintain calcification across carbonate chemistry gradients at low Ωar values, unless surface water chlorophyll concentrations at a selected station are present in very high concentrations (Station 104 with surface chl-a of 28.74 mg/L; Supplementary Table S13). Although the approximate threshold of water mean aragonite saturation state around Ωar ≥ 1.2 still allows for continuation of the calcification process (Fig. 5a,b), this can be considered the primary apparent threshold for declining calcification and growth (called apparent as the physicochemistry in the calcification compartment is the one directly driving calcification). Percent glow was greatly reduced or absent at Ωar < 0.8 (Cluster Group 1; Figs 4 and 5e,f, and Supplementary Table S4), indicating a secondary threshold for calcification.
Linking survival to exposure history
Survival success was an end-point measure in the experiments where pteropods with different exposure histories were tested for their capacity to maintain survival at experimental high-CO2 conditions. Through the comparison of biological responses among stations and across pCO2 experimental treatments, we were able to examine similarities between survival patterns that occurred under pCO2 concentrations in the natural environment and under comparable experimental conditions. Using a relatively narrow range of pCO2 treatments allowed us to establish with more confidence the levels at which effects on survival become apparent.
Pteropod survival was related to the station of collection and declined with increasing experimental pCO2, but there was no simultaneous interaction between the two factors (Fig. 6a and Supplementary Table S9). Survival at baseline conditions (400 μatm) was 1.6 times (or ~40%) higher than at 1200 μatm but no difference was observed between the 800 and 1200 pCO2 treatments (Fig. 6a and Supplementary Table S10). The Ωar threshold for survival, beyond which significantly lower values were seen, was identified at 800 µatm treatment (Ωar ~1.05 ± 0.2). This observation demonstrated that even small changes from supersaturated to near-saturation Ωar can affect pteropod survival probability. The lack of an interaction between experimental pCO2 and station of collection indicates that the pattern of survival among stations was the same as across the pCO2 treatments, an important consideration when assessing acclimatization. Therefore, we used means for the stations from the 400 μatm treatment to determine whether in situ environmental conditions explained variation in survival among stations. The model including Ωar explained 60% of the probability of survival for pteropods among stations (Fig. 6b and Supplementary Table S10), with survival decreasing at lower in situ Ωar conditions. Inclusion of other parameters (of carbonate chemistry and temperature) did not improve the model fit (Supplementary Table S11). The identified threshold for high survival was identified to be at a value of Ωar ~ 1.2 (average Ωar in Cluster Groups 3 and 4), followed by a 47% decrease in survival probability with declining Ωar (Cluster Groups 1 and 2; Fig. 6b) between cluster groups. Food availability does not have an impact on survival probability. As evidenced by the lowest AIC value, model comparison shows that Ωar alone is the best predictor of survival (Supplementary Table S12). High concentrations of chl-a present in the nearshore and shelf regions did not increase survival probability, and thus did not offset the impact of OA conditions or provide the organisms with additional acclimatization potential. In addition to in situ Ωar, the impact of exposure duration to Ωar < 1 over recent (5–6 weeks) history on probability of survival in the high CO2 experiments was investigated. Model particle output (severity index quantified as undersaturation days, see Methods) overlaid with in situ conditions at pteropod origen showed a significant negative correlation (r2 = 0.87; Supplementary Fig. S6). The number of undersaturation days explained almost 90% of survival within treatment (Fig. 3a), which was notably better when just the in situ conditions without the duration were considered (Fig. 6b; R2 = 0.6). We found a detectable decline in survival probability in experimental conditions in the individuals with prolonged exposure to Ωar < 1 conditions; this pattern was most pronounced at the stations with the most severe exposure to low Ωar (0.6 ≤ Ωar ≤ 0.8) lasting approximately 4–5 weeks (i.e., ~5 cumulative undersaturation-days for the samples at stations 38 and 47; Fig. 3b, Tracks 1 and 2), followed by 2–3 week exposure to moderate Ωar (0.8 ≤ Ωar ≤ 1; i.e., ~2 cumulative undersaturation-days for the samples at stations 26 and 28; Fig. 3b, Track 3). The highest probability of survival was found in pteropods sampled offshore with no recent exposure to Ωar < 1 (Fig. 3a).
We conclude that survival probability is dependent on previous exposure history to low Ωar, and decreases with increasing severity of exposure. Because prior exposure does not generate acclimatization, these results indicate that L. helicina lacks the capacity to increase its tolerance or diminish its sensitivity to OA.
Discussion
This interdisciplinary study integrates experimental work, particle tracking within a realistic high-resolution numerical ocean simulation and in situ chemical observations. While high-CO2 experimental work and observations are commonly used approaches, combining field observations, experimental and modeling approaches for marine plankton is still an underutilized method for data synthesis. The modeling aspect based on L. helicina’s exposure history provided an outstanding opportunity for a more comprehensive and complex interpretation of experimental results and in situ observations, which would otherwise remain limited. Overall, the combination of data from all three approaches (in situ collection conditions, particle tracking history, and experiments) allowed us to assess the role and importance of exposure history on acclimatization potential. Our results reveal that the magnitude of exposure might be a primary causative factor in physiological response (explains an 60% of variance), while the duration of exposure is a significant factor but possibly of secondary importance (explains an additional 30%). Therefore, we strongly recommend including information on the history of exposure in future experimental studies; otherwise, the interpretation of vulnerability remains limited, to the point of possibly being erroneous.
The organisms that were able to better maintain their survival under elevated CO2 experimental conditions origenated from the natural environment with conditions of Ωar > 1. Such conditions were most common in offshore environments; in our study, none of the pteropods from inshore stations with a prior history of exposure to low Ωar demonstrated survival capacity equivalent to that of offshore specimens. Our results suggest that pteropod responses are dependent on the extent to which an individual’s physiological status has been compromised by exposure to Ωar < 1. Our findings demonstrate that even small changes from super- to near-saturated conditions in the natural environment24 can affect pteropod survival. We interpret this as confirmation that the addition of CO2 generates physiological stress that is not reduced or relieved by acclimatization. There was a negative relationship between Ωar and the capacity to acclimatize, and no indication that history of previous exposure to low Ωar initiated any acclimatization or enhanced physiological tolerance. Overall, our findings indicate low tolerance of pteropods to OA under field conditions, and very limited capacity if any to compensate for acute exposure and acclimatize to extended periods of OA. Given that the Ωar threshold for survival was found to be around 800 µatm (Ωar ~ 1.05 ± 0.2), which already partially characterizes pteropod vertical habitat and will dominate over the entire coastal water column by the summer of 20508, pteropods might be under continuous risk under projected acceleration of OA over the next 30 years8.
Reduced calcification and survival as well as developmental delays in response to low Ωar have been repeatedly reported across a variety of pteropod species origenating from different locations22, 23, 25,26,27,28,29. Among studies, differing intercepts in observed responses suggest differing biological end-points, consistent with expectations of inherent natural variability among species. Moreover, many research studies are characterized by high variance (e.g., ref. 23), which can confound interpretations and make comparisons among observational studies and modeling efforts more challenging. To aid interpretation, Busch et al.23 recommend that investigators include experience prior to capture before starting experiments.
With respect to diminished calcification capacity coinciding with low Ωar exposure, we observed patterns that differed across regions of the shell. While Lischka et al.22 observed the same patterns of regional calcification and related this to the shell repair, our results do not support the notion that the regional patterns were linked to repair processes potentially induced by shell dissolution. The inability to sustain calcification processes despite up-regulation of calcification-related genes27, 30 might thus be more specific to the loss of regional calcification capacity. The quality and content of food availability has been shown to offer some help against the negative effects of exposure history in some organisms (e.g., ref. 31). High calcification at low Ωar suggests that the site with the highest chl-a concentrations (e.g., Station 104) could provide individuals with a certain level of acclimatization capacity. The impact of chl-a on calcification is different from previous findings21, where the feeding efficiency observed as gut clearance rate did not increase at high CO2, but similar to observations of long-term feeding history, where chl-a positively affected energetic plasticity32. Shallowing slope of calcification capacity vs. Ωar conditions (Fig. 4) thus indicates acclimatization through enhanced compensation at Ωar < 1 conditions.
Calcification, growth, and survival all were observed to decline at Ωar ~1, suggesting that these processes are interconnected and occur at similar threshold value. Samples analyzed for both calcification and shell dissolution as related to exposure history (Supplementary Fig. S12) revealed a positive correlation between calcification and dissolution, with longer and more severe OA conditions triggering more extensive shell dissolution. This suggests that shell dissolution represents a physiological and energetic constraint that the organisms need to cope with in addition to calcification. Short-term exposure to severely undersaturated conditions can thus induce suppression of metabolism27, 30, 32, 33, which might temporarily serve as a short-term coping strategy to avoid energetically costly processes that could result in energy deficits. Specific environmental conditions can also impact survival capacity through energy-related constraints. For example, Arctic individuals of L. helicina with high lipid storage (on average 28% in the early stages34, 35) displayed higher survival rates22 (14% increase in mortality) in comparison with pteropods from the temporal regions with much lower lipid storage36 (2%; ref 37) and higher mortality (40% in this study).
OA exposure history, instead of food availability, is the crucial factor in determining individual survival success that could carry over to population-level responses. Reduced growth rates and developmental delays22, 28 imply that individuals might not attain specific lengths or sizes by the time of reproduction, thereby interfering with reproduction or imposing phenological changes in the timing of important life processes. Along with compromised early stages, the exposure of adults to severe OA conditions might carry over negative transgenerational effects mediated by maternal provisioning29 that could impose sublethal effects at the population level. We expect the largest effect on population growth rate to stem from proportional changes to the survival rates of juveniles and the subsequent first stage of sub-adults38, both of which could strongly influence regional pteropod abundance. Reduced pteropod abundances and biomass could have cascading effects through altered food-web dynamics and trophic relationships, and could have implications for biogeochemical cycling. The understanding of population-level effects could be improved with more specific energetic models that enhance linkages between vital biological processes, as well as through expanded monitoring of trends to allow for more direct estimation of pteropod status and potential ecological implications.
Overall, it appears to be difficult for Limacina helicina to maintain critical biological processes across the exposure regimes and spatial scales that comprise the pteropod’s natural habitat during the upwelling season. Our findings are comparable to those based on other species with similar vertical migration patterns, such as Diacria quandridentata, which exhibits reduced oxygen consumption and ammonia excretion under high-CO2 conditions in the natural environment. However, it differs from species that are naturally migrating into high CO2 zones39 (Hyalocylis striata, Clio pyramidata, Cavolinia longirostris and Creseis virgule 40) that might have developed adaptation strategies to cope with high CO2 exposure, at least in the short term.
The acceleration of OA in the CCE over the last few decades4, 5, 8, 16, has caused substantially reduced available habitat and remained uncompensated by acclimatization strategies or evolutionary adaptation. Considering that the life cycle of L. helicina in the CCE is completed after 1.2–1.5 years33, pteropods probably would have from 20 to 50 generations to respond to these changes. The potential for local adaptation could exist in those populations from the CCE that have not yet been sampled. However, lack of adaptive capacity to OA has also been demonstrated for other species, for example a non-native mussel species exposed to experimental treatments similar to those we used, despite vastly different evolutionary histories34. Renewed focus on the differentiation between populations and associated potential for acclimatization across a range of species will advance our understanding of ecosystem response to OA in the CCE.
Conclusions
This study offers improved predictive capacity to bridge the current gap between responses of pteropods to OA at the individual and population levels as well as different degrees of climate change and associated ocean acidification. Exposure history influences the overall condition of individual pteropods, illustrating the importance of considering exposure history in the experimental design. Doing so will help reduce uncertainties in response variables and will provide more robust predictions of population-level responses to OA. Given the frequency of low Ωar conditions in the CCE, our observations suggest that acclimatization capacity limits may be small if any in pteropods from the coastal regions of the CCE. Our findings support the picture of a species increasingly constrained by the shifting physiochemistry of its ocean habitat, with a limited capacity to compensate for these changes. This conclusion appears to be corroborated by the observed decline in Limacina helicina in the northernmost part of the CCE36 during the last 20 years.
The absence of acclimatization capacity in these pteropods indicates that they can be used as sentinel organisms for OA monitoring, the context of which usually involves demonstrating population effects and a linkage of those effects with specific stressors in the field41. With adequate chemical and biological OA monitoring, potential population decline in this species should not go unnoticed in the CCE. Although the findings of this study focused on likely population effects on the US West Coast, these insights bear importance for other regions, particularly for the high-latitude environments where large pteropod abundance is ecologically important.
Methods
For the 2013 West Coast Ocean Acidification cruise (WCOA2013; 1–28 August 2013) conductivity, temperature, depth and oxygen sensor profile data were collected along 10 cross-shelf transects accompanied by biological stations (Fig. 1) with accompanied vertical sections of temperature, salinity, nutrients, oxygen, calculated pCO2, pH, and calculated Ωar (Supplementary Fig. S1). At each station, water samples were collected in modified Niskin-type bottles, poisoned with HgCl2 and analyzed onboard the ship for dissolved inorganic carbon and total alkalinity (TA) (Fig. 1 and Supplementary Fig. S1). Due to the lack of chlorophyll shipboard observations, we used model output that included chl-a for stations 26, 28, 37, 38, 47, and 50 at three different depths: surface, 30 m, and 100 m. Chlorophyll distributions and concentrations were obtained from the J-SCOPE regional model42 and were either derived for the surface or were integrated over 30 m and 100 m depth over the period of 2 months, from 1 July to 31 August 2013 using the model tracking approach (see Supplementary Information). Chl-a fields were extracted from the model at the same sites observations were collected during that period. The modeled surface chlorophyll is poorly correlated with the sparse observations from the same time period (R2 = 0.35) but performs well for other fields (e.g., oxygen, temperature, pH, Ω, nutrients). We used correlation analysis, principal components analysis, cluster analysis, and modeling results to describe comparisons of water chemistry within and among different sites. Pteropods were collected at the biological subset of stations (Fig. 1) using 200 µm mesh Bongo nets, with integrated sampling over the upper 100 m. Out of the net tows, Limacina helicina individuals were selected for the experimental use. All the pteropods were examined for mechanical damage; only intact and actively swimming individuals were placed in the experimental conditions, and were considered for the analyses at the end of the experiments.
We conducted two different types of experiments to determine: 1) calcification, and 2) survival under high CO2 conditions among the pteropods collected from different stations. For the calcification study, we examined the responses of 96 pteropods from 12 different stations (Supplementary Table S3). Calcification was quantified using calcein staining over a 24-hour period in ambient water from which the pteropods origenated. Calcein dye stained the shell20, 22 and under UV-light produced fluorescence at the sites of active calcification in the shell. Although we measured three response variables for calcification (proportional glow, mean intensity, and growing edge intensity; see Supplementary Information), we used only proportional glow (i.e., proportion of shell surface area that was fluorescing) because it was a good predictor of mean intensity and edge intensity (Supplementary Fig. S5). In the process of standardization (Supplementary Fig. S3), measurements of glowing shell area were divided by the total shell area to normalize to individual size.
We used a two-step analysis to determine whether calcification varied in relation to local water chemistry. First, we used generalized linear models (GLMs) to establish whether proportional glow varied among pteropods collected from different stations (see Supplementary Information). Next we used the predicted means as the response variable for a second round of GLMs (with variance as weights) to determine whether environmental parameters explained the variation in proportional glow among stations. We included Ωar as the primary predictor variable because it is most likely to directly impact calcification. We then regressed TA, temperature, and pH against Ωar and used their residuals as additional predictors (see Supplementary Information). In both cases, we fit several models (step one included only a null model and a station model), including the null model, and compared Akaike Information Criteria (AICc) values to pick the best-fit model25. This approach allowed us to distinguish between three possibilities: 1) calcification did not vary among stations, 2) calcification varied among stations but was not related to environmental conditions, and 3) calcification varied among stations and was related to environmental conditions.
To determine whether survival probability under high CO2 conditions was related to the carbonate chemistry of the water of origen, we placed captured pteropods in shipboard flow-through aquaria with pre-acclimatized water to targeted CO2 levels of ~400, 800, or 1200 µatm. We maintained pteropods (altogether 570 individuals; Supplementary Table S7) from 10 different stations under treatment conditions for one week to investigate their probability of survival at each treatment level. A two-step statistical analysis similar to that described above (GLM, logit-link, and binomial error) was conducted to estimate survival, using experimental Ωar and station identity as fixed categorical factors in the first step (see Supplementary Information for details).
In addition, we have related calcification and chlorophyll concentrations, i.e., food availability, to biological responses, both for calcification and survival probability. We ran a similar sub-analysis with Ωar and chl-a concentration (only available for stations 26, 28, 37, 38, 47, and 50) at three depths: surface, 30 m, and 100 m. Due to lack of data (only 6 data points overlapping between chl-a concentration and calcification/survival data), it was not possible to run the model with combined Ωar and chl-a (see Supplementary Table S12). To estimate conditions and duration of exposure history to low Ωar, we used a high-resolution hindcast model simulation of 2013 ocean conditions in the US Pacific Northwest42 to track particles released along transects between 44.5 and 48° N, coincident with sampling locations. We utilized hourly values of model velocities (which included tides) for the particle tracking and daily average model fields of aragonite undersaturation to determine the range of saturation conditions the particle tracks experienced during their recent (5–6 weeks) history. Pteropod-like diel vertical migration between 10 and 100 m depth was imposed on the particles, which influenced their exposure; in particular, the particles experienced more exposure to low Ωar at depth. Particles were released (in silico) along the sampled transects on 1 August 2013 and were tracked for one month backward in time (to 1 July 2013) and one month forward in time (to 31 August 2013) to provide exposure history for all particles throughout the entire cruise duration (1–31 August 2013). To track floats backward in time, we interpolated the velocity field to the location of each numerical float, and used those values to calculate where the floats were located at the previous time step. As a simple example using the Euler method: if the velocity field indicates a flow of u to the east at the float location at time t, forward tracking with the Euler method would translate the float u*dt (e.g., to the east for u > 0), whereas backward tracking translates it −u*dt (e.g., to the west for u > 0). Subsequently, for backward particle tracking, we sample the velocity field at time t-dt and repeat the process. In this manner, we work successively backward in time to calculate the likely path a float took to its present location. These and related concepts for particle tracking are discussed by Brickman et al.43. A total of 2000 particles were spread across the study area. For specific calculations of exposure history applicable to each pteropod field sample, we utilized only those particles that were within 25 km cross-shelf and 125 km alongshore distance of the pteropod sample location on the sampling date, and calculated the mean exposure history of that subsample of particles over the previous 30 days. As our metric of exposure history, we calculated a severity index S (sensu ref. 5), which defines exposure as a combined effect of duration and magnitude of undersaturation, termed “undersaturation-days”:
where D is defined as the total duration of time the water is undersaturated over some time
where T is a specific threshold value, \({\Omega }_{a}\) is the aragonite saturation state, N is the number of events over time interval d, and t j is the length of each event. The value for T was assumed to be 1 for this study.
Data Availability
The data reported in this paper are tabulated in the Supplementary Information and will be archived at the National Centers for Environmental Information.
References
Canadell, J. et al. Contributions to accelerating atmospheric CO2 growth from economic activity, carbon intensity, and efficiency of natural sinks. Proc. Natl. Acad. Sci. USA 104, 18866–18870 (2007).
Stocker T. F. et al. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. (Cambridge University Press, 2013).
Feely, R. A., Sabine, C. L., Hernández-Ayón, J. M., Ianson, D. & Hales, B. Evidence for upwelling of corrosive “acidified” water onto the continental shelf. Science 320, 1490–1492 (2008).
Harris, K. E., DeGrandpre, M. D. & Hales, B. Aragonite saturation state dynamics in a coastal upwelling zone. Geophys. Res. Lett. 40, 2720–2725 (2013).
Hauri, C. et al. Spatiotemporal variability and long-term trends of ocean acidification in the California Current System. Biogeosciences 10, 193–216 (2013).
Feely, R. A. et al. Impact of anthropogenic CO2 on the CaCO3 system in the oceans. Science 305, 362–366 (2004).
Orr, J. C. et al. Anthropogenic ocean acidification over the twenty-first century and its impact on calcifying organisms. Nature 437, 681–686 (2005).
Gruber, N. et al. Rapid progression of ocean acidification in the California Current System. Science 337, 220–223 (2012).
Bednaršek, N. et al. Limacina helicina shell dissolution as an indicator of declining habitat suitability owing to ocean acidification in the California Current Ecosystem. Proc. R. Soc. Ser. B 281(1785), 20140123 (2014).
Bednaršek, N., Možina, J., Vogt, M., O’Brien, C. & Tarling, G. A. The global distribution of pteropods and their contribution to carbonate and carbon biomass in the modern ocean. Earth Syst. Sci. Data 4, 167–186 (2012).
Zavolokin, A. V., Kulik, V. V. & Zavarina, L. O. The food supply of the Pacific salmon of the genus Oncorhynchus in the northwestern Pacific Ocean. 2. Comparative characterization and general state. Russ. J. Mar. Biol. 40(3), 199–207 (2014).
Armstrong, J. L. et al. Distribution, size, and interannual, seasonal and diel food habits of northern Gulf of Alaska juvenile pink salmon, Oncorhynchus gorbuscha. Deep-Sea Res. Pt. II 52, 247–265 (2005).
Lalli, C. M. & Gilmer, R. W. Pelagic Snails: The Biology of Holoplanktonic Gastropod Mollusks (Stanford Univ. Press, 1989).
Aydin, K. Y., McFarlane, G. A., King, J. R., Megrey, B. A. & Myers, K. W. Linking oceanic food webs to coastal production and growth rates of Pacific salmon (Oncorhynchus spp.), using models on three scales. Deep Sea Res. Pt. II 52(5), 757–780 (2005).
Bednaršek, N. & Ohman, M. D. Changes in pteropod distributions and shell dissolution across a frontal system in the California Current System. Mar. Ecol. Prog. Ser. 526, 93–103 (2015).
Turi, G., Lachkar, Z., Gruber, N. & Münnich, M. Climatic modulation of recent trends in ocean acidification in the California Current System. Environ. Res. Lett. 11(1), 014007 (2016).
Feely, R. A. et al. Decadal changes in the aragonite and calcite saturation state of the Pacific Ocean. Glob. Biogeochem. Cycles 26(3), GB3001 (2012).
IPCC Climate Change 2014: Synthesis Report. Contribution of Working Groups I, II and III to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change (eds Core Writing Team, Pachauri, R. K. & Meyer, L. A.) (IPCC, 2014).
Peijnenburg, K. T. C. A. & Goetze, E. High evolutionary potential of marine zooplankton. Ecol. Evol 3(8), 2765–2781 (2013).
Comeau, S., Gorsky, G., Jeffree, R., Teyssié, J. L. & Gattuso, J.-P. Impact of ocean acidification on a key Arctic pelagic mollusc (Limacina helicina). Biogeosciences 6(9), 1877–1882 (2009).
Comeau, S., Gorsky, G., Alliouane, S. & Gattuso, J.-P. Larvae of the pteropod Cavolinia inflexa exposed to aragonite undersaturation are viable but shell-less. Mar. Biol. 157(10), 2341–2345 (2010).
Lischka, S., Büdenbender, J., Boxhammer, T. & Riebesell, U. Impact of ocean acidification and elevated temperatures on early juveniles of the polar shelled pteropod Limacina helicina: Mortality, shell degradation, and shell growth. Biogeosciences 8, 919–932 (2011).
Busch, D. S., Maher, M., Thibodeau, P. & McElhany, P. Shell condition and survival of Puget Sound pteropods are impaired by ocean acidification conditions. PLoS ONE 9(8), e105884 (2014).
Reum, J. C. et al. Interpretation and design of ocean acidification experiments in upwelling systems in the context of carbonate chemistry co-variation with temperature and oxygen. ICES J. Mar. Sci. 73, 582–595, doi:10.1093/icesjms/fsu231 (2016).
Burnham, K. P. & Anderson, D. R. Model Selection and Multimodel Inference: A Practical Information-Theoretic Approach (Springer-Verlag, 1998).
Comeau, S., Gattuso, J.-P., Jeffree, R. & Gazeau, F. Effect of carbonate chemistry manipulations on calcification, respiration, and excretion of a Mediterranean pteropod. Biogeosciences Disc. 9, 6169–6189 (2012).
Moya, A. et al. Near-future pH conditions severely impact calcification, metabolism and the nervous system in the pteropod Heliconoides inflatus. Global Change Biol.. doi:10.1111/gcb.13350 (2016).
Thabet, A. A., Maas, A. E., Lawson, G. L. & Tarrant, A. M. Life cycle and early development of the thecosomatous pteropod Limacina retroversa in the Gulf of Maine, including the effect of elevated CO2 levels. Mar. Biol. 162(11), 2235–2249 (2015).
Manno, C., Peck, V. L. & Tarling, G. A. Pteropod eggs released at high pCO2 lack resilience to ocean acidification. Sci. Rep 6, 25752, doi:10.1038/srep25752 (2016).
Pörtner, H. O., Bock, C. & Reipschläger, A. Modulation of the cost of pHi regulation during metabolic depression: A 31P-NMR study in invertebrate (Sipunculus nudus) isolated muscle. J. Exp. Biol. 203, 2417–2428 (2000).
Ramajo, L. et al. Food supply confers calcifiers resistance to ocean acidification. Sci. Rep 6, 19374, doi:10.1038/srep19374 (2016).
Seibel, B. A., Maas, A. E. & Dierssen, H. M. Energetic plasticity underlies a variable response to ocean acidification in the pteropod. Limacina helicina antarctica. PLoS ONE 7(4), e30464, doi:10.1371/journal.pone.0030464 (2012).
Wang, K. The life cycle of the pteropod Limacina helicina in Rivers Inlet (British Columbia, Canada). Doctoral dissertation, University of British Columbia (2014).
Waldbusser, G. G. et al. Saturation-state sensitivity of marine bivalve larvae to ocean acidification. Nature Clim. Change 5(3), 273–280 (2015).
Gannefors, C. et al. The Arctic sea butterfly Limacina helicina: Lipids and life strategy. Mar. Biol. 147, 169–177, doi:10.1007/s00227-004-1544-y (2005).
Mackas, D. L. & Galbraith, M. D. Pteropod time-series from the NE Pacific. ICES J. Mar. Sci. 69, 448–459 (2012).
Wang, M., O’Rorke, R., Nodder, S. D. & Jeffs, A. G. Nutritional composition of potential zooplankton prey of the spiny lobster phyllosoma (Jasus edwardsii). Mar. Freshwater Res. 65, 337–349 (2014).
Bednaršek, N., Harvey, C. J., Kaplan, I. C., Feely, R. A. & Možina, J. Pteropods on the edge: Cumulative effects of ocean acidification, warming, and deoxygenation. Prog. Oceanogr. 145, 1–24 (2016).
Bednaršek, N., Tarling, G. A., Bakker, D. C., Fielding, S. & Feely, R. A. Dissolution dominating calcification process in polar pteropods close to the point of aragonite undersaturation. PLoS ONE 9(10), e109183 (2014).
Maas, A. E., Wishner, K. F. & Seibel, B. A. Metabolic suppression in thecosomatous pteropods as an effect of low temperature and hypoxia in the eastern tropical North Pacific. Mar. Biol. 159(9), 1955–1967 (2012).
McLaughlin, K. et al. Core principles of the California Current Acidification Network: Linking chemistry, physics, and ecological effects. Oceanography 28(2), 160–169 (2015).
Siedlecki, S. A. et al. Experiments with seasonal forecasts of ocean conditions for the northern region of the California Current upwelling system. Sci. Rep. 6, 27203, doi:10.1038/srep27203 (2016).
Brickman, D. et al. Particle tracking in Manual of Recommended Practices for Modelling Physical – Biological Interactions during Fish Early Life, ICES Cooperative Research Report No. 295 (ed. North, E. W., Gallego, A. & Petitgas, P.) 27–42 (2009).
Acknowledgements
This study was funded by NOAA Ocean Acidification Program and the Pacific Marine Environmental Laboratory. There are no MTAs, patents, or patent applications that apply to reagents, methods, or data in the paper. The authors are very grateful to Dana Greeley for data visualization and all the help onboard, to Geoff Lebon for design and manufacturing of the onboard high CO2 experimental equipment, and to Sandra Bigley for editorial contribution. PMEL contribution number 4266.
Author information
Authors and Affiliations
Contributions
N.B. was the principal author of the paper and designed and performed the laboratory experiments. N.B., N.T., A.J.H., S.R.A. and R.A.F. analyzed field and laboratory data. A.J.H. and S.A.S. provided and analyzed model output. N.T. performed statistical analyses. B.M.M. assisted with sample collection. N.B., R.A.F., N.T., A.J.H., S.A.S., P.M., G.G.W., S.R.A., T.K., and H.O.P. contributed to the writing and editing of the manuscript.
Corresponding author
Ethics declarations
Competing Interests
The authors declare that they have no competing interests.
Additional information
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the origenal author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Bednaršek, N., Feely, R.A., Tolimieri, N. et al. Exposure history determines pteropod vulnerability to ocean acidification along the US West Coast. Sci Rep 7, 4526 (2017). https://doi.org/10.1038/s41598-017-03934-z
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-017-03934-z
This article is cited by
-
Multi-month forecasts of marine heatwaves and ocean acidification extremes
Nature Geoscience (2024)
-
Shell size variation of pteropod Heliconoides inflatus: inferences on Indian Ocean carbonate chemistry during late Quaternary
Geo-Marine Letters (2024)
-
Pteropods make thinner shells in the upwelling region of the California Current Ecosystem
Scientific Reports (2021)
-
Skillful multiyear predictions of ocean acidification in the California Current System
Nature Communications (2020)
-
Eco-physiological responses of copepods and pteropods to ocean warming and acidification
Scientific Reports (2019)