The accumulation of anthropogenic CO₂ in the ocean has major ecological, socioeconomic, and biogeochemical impacts, with repercussions for the ocean as a critical carbon sink. Ocean acidification (OA) disproportionately affects marine calcifiers, among which pelagic zooplanktonic pteropods play a significant role in carbonate export. The pteropod, due to the susceptibility of its aragonite shell to rapid dissolution, is one of most vulnerable groups and a key indicator for OA regional monitoring, but its regional sensitivities have not yet been extrapolated over global scales. To delineate spatial and temporal changes in pteropod shell dissolution, global OA status and the OA rate of change were evaluated, based on gridded climatologies of observations and using a Regional Ocean Modeling System (ROMS) biogeochemical/​ecosystem model. Pteropods, which dominate in the polar and subpolar regions, are characterized by low aragonite saturation state and low buffering capacity, where extended pteropod subsurface dissolution is projected. We show that pteropods are most susceptible to OA in the polar regions, subpolar North Pacific, and eastern boundary upwelling system regions, particularly the California and Humboldt Current Systems. Rates of acidification and corresponding increases in pteropod shell dissolution are projected to be the fastest in the North and South Equatorial Currents. These are the regions with the greatest impacts on pteropods and biogeochemical implications.

 

Example caption text > High res figure

Introduction

Accumulations of anthropogenic CO₂ in Earth’s atmosphere and ocean pose a major threat to global climate and society (Friedlingstein et al., 2020, 2022). The crucial role of the ocean in attenuating the increase in atmospheric CO₂, and thus global warming, is related to the excess atmospheric CO₂ that is absorbed by the ocean as dissolved CO₂, which alters the seawater carbonate system (Zeebe and Wolf-Gladrow, 2001). Increased CO₂ uptake causes a decline in pH, shoaling of the CaCO₃ saturation state horizon (Ω), and reduced carbonate ion concentration in a process known as ocean acidification (OA; Feely et al., 2004, and 2023, in this issue). Among global ocean change processes, OA represents one of the greatest threats to ocean ecosystems, marine-related socioeconomic activities (IPCC, 2022), and biogeochemical processes, with significant repercussions for the ocean carbon sink, because it alters transfer of carbon from the surface to the ocean depths (Feely et al., 2002; Lee and Feely, 2021).

The major mechanisms of carbon transport in the ocean’s interior include the inorganic solubility pump and the biological pump, within which the carbonate pump is an important transport process that leads to enhanced carbon sequestration on timescales of several centuries to millennia. Pelagic calcifiers, including coccolithophores, foraminifera, and shelled thecosomata (referred to as pteropods), are significant climate regulators (Milliman, 1993; Sarmiento and Gruber, 2006). Through their contributions to CaCO₃ export, calcifiers affect the rate, magnitude, and strength of the ocean as a carbon sink and its potential for affecting climate change (Feely et al., 2004; Gehlen et al., 2007). OA impacts the carbonate pump through several major pathways (Figure 1). The OA-induced decrease in Ω increases CaCO₃ dissolution, reduces calcification, and reduces the sinking speed of biogenic aggregates. The dissolution of biogenic carbonates at depth leads to decreased availability of CaCO₃ ballast material to accelerate sinking velocities of particulate organic carbon (Honjo et al., 2008). To understand the scale of the changes, regional biogeochemical processes should be extrapolated to global scale.

 

FIGURE 1. Provisional conceptualization of the major processes affecting pteropods under (a) current conditions, and (b) intensifying ocean acidification (OA). Arrows indicate direction, and arrow boldness indicates potential change in magnitude. Only processes relevant to the carbonate pump are shown. The thin horizontal line marks the aragonite saturation horizon, which increases over time due to ocean acidification. TA = Total alkalinity. > High res figure

 

Pteropods are zooplanktonic pelagic calcifiers with ubiquitous biogeographic distribution across the global ocean; their highest biomass occurs in the upper 200 m (Bednaršek et al., 2012a) during their diel vertical migration. They build shells of aragonite, a metastable form of calcium carbonate that is 50% more soluble than calcite (Mucci, 1983), which makes them more sensitive to OA than calcite-shelled organisms (Fabry et al., 2008; Bednaršek et al., 2014b; Manno et al., 2017). Observations and modeling suggest that pteropods contribute significantly more than previously thought to global carbonate export. Bednaršek et al. (2012a) estimated that pteropod contributions range between 20% and 42%, while a modeling study by Buitenhuis et al. (2019) found that pteropods dominate CaCO₃ export, ranging between 33% and 89% in the upper 200 m.

Apart from their biogeochemical role, pteropods are also an important component of the food web because they channel energy through the trophic levels and provide essential food resources for economically important fish. To date, experimental studies have demonstrated significant rate changes in biomineralization responses to reduced aragonite saturation state (Ω_ar), showing increased shell dissolution and reduced calcification (Lischka et al., 2011; Comeau et al., 2009, 2012; Bednaršek et al., 2014a, 2017a). This is further supported by expert-​consensus selected thresholds that delineate the magnitude and duration of OA exposure involved in negative responses (Bednaršek et al., 2019). Research on NOAA Pacific Marine Envi­ron­mental Laboratory (PMEL) cruises demonstrated that severe shell dissolution was already observable in the sensitive subpolar, polar, and upwelling regions under current conditions (Bednaršek et al., 2012b, 2014b, 2021; Niemi et al., 2021). This makes pteropods one of the most vulnerable groups to OA and a key indicator for OA vulnerability assessment and regional monitoring (Bednaršek et al., 2017b). Such rapid changes in biomineralization processes can have significant biogeochemical implications for carbon export, especially in regions where their biomass is high—the Arctic and Southern Oceans, the subpolar North Pacific, the equatorial Pacific, and the coastal Eastern Boundary Upwelling Systems, particularly the California Current System and the Humboldt Current System (Knecht et al., 2023; Bednaršek et al., 2012a).

Although pteropods are known to be important players in the global carbon budget, there are still a number of uncertainties regarding the drivers, the strength, and the significance of pteropods in the carbonate pump. Our understanding of pteropods’ aragonite response remains limited, as reflected in the lack of their representation in CaCO₃ cycle assessments and modeling. Refining our understanding of these processes means developing significantly better parameterizations of critical pteropod responses to OA and their changes across global scales. This paper, which is based on research conducted over the last 12 years as part of PMEL climate and acidification research cruises, summarizes the current status and 35-year trends for the best understood and most observed process of pteropod shell dissolution, which has direct consequences for carbon export and sequestration. By using the quantified sensitivity of biological responses (Bednaršek et al., 2014a,b; Feely et al., 2016), integrated chemical observations, and model outputs, we align the changes in Ω_ar in the upper 200 m with corresponding changes in shell dissolution and then scale up from regional to global scales. We hypothesize that the greatest changes in shell dissolution, and thus carbon export, will be observed in the areas where the highest pteropod biomass overlaps with the lowest Ω_ar. Here, we quantify the changes in pteropod shell dissolution, aiming to improve our understanding of biogeochemical processes across both temporal and spatial scales and to provide insights on how they are affected by human-caused climate change. The biological responses to ocean acidification described here address one of PMEL’s overarching missions: to better understand the impacts of climate change and ocean acidification on marine ecosystems.

 

Methodology

We calculated climatologies of observed Ω_ar from published observations in the GLODAPv2.2019 observational data set of dissolved inorganic carbon (DIC), total alkalinity (TA), silicate, phosphate, temperature, and salinity. The horizontal resolution is 1°×1°, with 57 depth levels (0–1,500 m) and monthly resolution (12 months). CO2SYS was used in MATLAB to calculate Ω_ar. Input data for CO2SYS were obtained from the following published climatologies: DIC from Broullón et al. (2020); TA, silicate, and phosphate from Broullón et al. (2019); and WOA-2018 climatologies for temperature (Locarnini et al., 2018) and salinity (Zweng et al., 2018). The average date for observations of DIC and TA in GLODAPv2.2019 used by Broullón et al. (2019, 2020) is the year 2005, with a standard deviation of ± 9 years.

Therefore, our calculated climatologies of observed Ω_ar represent a mean year of 2005, and the range of years of about 1996–2014 based on the interval of the mean ± standard deviation of the dates when the observed DIC and TA data were collected. The annual average climatology of observed Ω_ar is estimated as the mean of the monthly average climatologies. Multivariate ENSO Index Version 2 (MEI.v2) was calculated based on the sources described in https://psl.noaa.gov/enso/mei/.

We use monthly output from a hindcast simulation (1984–2019) with a high-​resolution Regional Ocean Modeling System (ROMS) model coupled to a biogeochemical/​ecosystem model BEC (Marchesiello et al., 2003; Shchepetkin and McWilliams, 2005; Moore et al., 2013; Desmet et al., 2022). The model runs on a telescopic grid covering the Pacific basin and centered on one pole at the US West Coast, which allows mesoscale dynamics at the pole to be fully resolved and basin-wide oceanic and atmospheric teleconnection to be captured (Frischknecht et al., 2018). Ω_ar is computed offline using the MOCSY 2.0 routine (Orr and Epitalon, 2015) and the constants recommended by Dickson et al. (2007). The model output was validated against Ω_ar observations from the NOAA Ocean Acidification Program West Coast Cruises and the climatological GLODAPv2 1° × 1° gridded product over the latitudinal gradient of 25°–62°N (Desmet et al., 2022). We further evaluated the skill of hindcast Ω_ar from 0–200 m compared with observations in GLODAPv2.2022 from 20°S–65°N and found excellent agreement with GLODAPv2.2022 for both the Ω_ar using the climatology of gridded observations in the global ocean (mean bias = 0.099, RMSE = 0.26, R2 = 0.94, Nash-Sutcliffe Efficiency = 0.93), and the predicted Ω_ar using the ROMS-BEC model of the northeast Pacific Ocean by Desmet et al. (2022) (mean bias = 0.015, RMSE = 0.33, R² = 0.91, Nash-Sutcliffe Efficiency = 0.90). In addition, the mean rate of decrease in surface Ω_ar predicted by the model from 25°S–65°N and 70°–185°W from 1984 to 2019 equals the published mean rate of decrease in global surface Ω_ar from 1982 to 2020 (i.e., –0.07 per decade; Gregor and Gruber, 2021).

Trends in surface Ω_ar in the global ocean from 1982 to 2020 were calculated from OceanSODA-ETHZ (Gregor and Gruber, 2021), a gridded data set that was created by extrapolating surface ocean observations of pCO₂ from SOCAT (Bakker et al., 2016) and TA from GLODAP (Olsen et al., 2019) using a Geospatial Random Cluster Ensemble Regression (GRaCER) method (Gregor and Gruber, 2021). Gregor and Gruber (2021) reported near-zero negligible biases of calculated carbonate system variables compared with the observed data in GLODAP.

As pteropod biomass is mostly distributed in the upper 200 m (Bednaršek al., 2012a), we have used observational data and model outputs for the upper 200 m to delineate the changes in shell dissolution to Ω_ar. We relied on previous studies that determined pteropod spatial dominance (Bednaršek et al., 2012a; Knecht et al., 2023). We used the equation on the percent of pteropod individuals with severe shell dissolution = –66.29 × ln (Ω_ar) + 61.2 from Feely et al. (2016), which relies on two large-scale NOAA-supported West Coast cruise observations (2011 and 2013) of pteropod shell dissolution, previously evaluated for one cruise (in 2011) by Bednaršek et al. (2014b). Shell dissolution was described for the mid-stages (late juvenile and mid adult), referring to the juvenile and adult stages of the pteropod population (Bednaršek et al., 2016).

 

Results

Status and Trends of Aragonite Saturation State in the Upper 200 m

Climatologies based on observational data over the upper 200 m show a very large range of Ω_ar values across global scales (Figure 2). Some of the lowest Ω_ar values, which range from 0.95 to 1.5, are observed in the polar regions of the Arctic, including the Kara, Laptev, and Beaufort Seas, and the Arctic Ocean, as well as in the subpolar regions of the North Pacific between 45°N and 60°N, particularly the Okhotsk Sea, the Bering Sea, and the Gulf of Alaska, and in the Southern Ocean (Figures 2a and 3a,c). In the upwelling systems characterized by lower Ω_ar, particularly the California and Humboldt Current Systems, the annual Ω_ar values in closest coastal regions drop down to 1.2 and 1.5, respectively (Figure 1). Other locations of lower Ω_ar are the North and South Equatorial Current (NEC and SEC) regions between 5°N and 20°N and 0° and 15°S, where 0–200 m annual Ω_ar values are as low as 1.5–2 Ω_ar, thus exhibiting strong correlation with the temperature and the ENSO index, particularly in the equatorial Pacific, Monterey Bay, and Humboldt Current regions (Table 1).

 

FIGURE 2. The panels show (a) aragonite saturation state (Ωar) averaged over the upper 200 m based on observed data, and (b) projected percent of pteropods likely to show severe dissolution averaged over the 200 m water column. > High res figure

 

FIGURE 3. (left) Observation-based aragonite saturation state is shown (Ωar) averaged over the upper 200 m in the (a) Arctic Ocean from 50°N to 90°N, and (c) the Southern Ocean from 50°S to 90°S. (right) Percent of pteropods exhibiting severe dissolution averaged over the 200 m water column are charted for (b) the Arctic Ocean and (d) the Southern Ocean. > High res figure

 

TABLE 1. Annual averages for various locations in the North and South Pacific Oceans for the period 1984–2019 from Desmet et al. (2022), with Ωar ± stdev included over the upper 200 m, rate change of Ωar per decade, % of pteropods affected with dissolution ± stdev, and rate change of dissolution per decade. R2 values indicate goodness of fit, and p values indicate significance. Light gray values are p >0.05. > High res table

 

On the global level, the predicted trends for surface Ω_ar obtained from OceanSODA-ETHZ (Gregor and Gruber, 2021) for the relevant period show strong regional decreases (Figure 4a), which are regionally similar but spatially smaller in extent over the upper 200 m based on the ROMS-BEC model (Desmet et al., 2022; Figure 4b), resulting in additional percent pteropod shell dissolution increase (Figure 4c). The trends in surface Ω_ar for the California and Humboldt Current Systems are decreasing at rates that are generally more than 0.06 per decade (Figure 4a). Comparatively, in pteropod-relevant habitat over the upper 200 m, the slopes of the decline are significant for several regions (Table 1, Figure 4c). The tropical and subtropical regions from about 30°N to 30°S show the greatest Ω_ar decline of just over 0.1 Ω_ar unit per decade. In contrast, in the subpolar and polar regions the Ω_ar decline decreases to values <0.06 due, in part, to decreasing buffering capacity at the higher northern and southern latitudes. In subsurface waters, the narrow band of the equatorial Pacific from about 5°N to 5°S exhibits a decrease in Ω_ar of <0.05 over time as a result of upwelling of older water at deeper depths.

 

FIGURE 4. Shading shows changes in (a) surface Ωar from 1982 to 2020 based on gridded observations in Ocean-SODA-ETHZ (Gregor and Gruber, 2021), (b) Ωar per decade from 1984 to 2019 in the eastern Pacific averaged over the upper 200 m based on the ROMS-BEC model hindcast of Desmet et al. (2022), and (c) pteropod rate of severe dissolution averaged over the 200 m water column over the eastern Pacific Ocean based on the ROMS-BEC model hindcast of Desmet et al. (2022) ranging from 25°S to 65°N. The yellow dots in panels b and c indicate the stations for which more extensive analyses have been made (see Table 1). > High res figure

 

Pteropod Shell Dissolution: Status and Trends

Projections for the current status and rate of change in the percent of pteropods affected by severe shell dissolution are based on observed climatologies and ROMS-BEC model outputs, respectively. The extent of pteropod shell dissolution strongly overlaps with regions where current Ω_ar reaches its lowest values and where buffering capacities are low. In the tropical and subtropical Pacific, the dissolution rate also coincides with the fastest rate of Ω_ar decline in the NEC (Figure 4). As a result of exposure to Ω_ar in the upper 200 m, the extent of individuals currently undergoing severe shell dissolution is projected to be as much as 40% in the Southern Ocean south of 60°S, 50% in the North Pacific at 50°N, and 60% in the Arctic. Corresponding to Ω_ar conditions in the upwelling regimes of the California and Humboldt Current Systems, the projections show an offshore-onshore gradient of biological impacts, with up to 40% of individuals affected by dissolution in the highly productive coastal regions of both systems. It is worth noting that previous regional dissolution observations in the polar regions and the California Current System (Bednaršek et al., 2012b, 2014b, 2021; Niemi et al., 2021) closely align with the projections delineated in this study. Projections in the regions of the NEC and SEC, Western North Atlantic, and Benguela Current show that 10%–30% of pteropods will be impacted by dissolution. The projected increase in pteropods impacted by dissolution ranges between 1% and 4%, indicating that an additional 4% to 14% of the total number of pteropods exhibited severe dissolution over the 1984–2019 period. Increases in dissolution are closely correlated with the decline in Ω_ar recorded over the 36-year period investigated; the regression over time is significant for most of the regions (Table 1) and explains from 24% to 53% of the variability in Ω_ar.

 

Discussion

Global projections show strong agreement between Ω_ar and the extent of pteropods impacted by severe shell dissolution under current conditions. This suggests that some inferences can be made about how increasing OA has impacted dissolution, noting the importance of current Ω_ar conditions, the rate of change of these conditions, and the corresponding biological responses. The results show that the regions with currently low Ω_ar are not the same as the regions with the greatest change; the latter actually are on the margins of the low Ω_ar regions, for example, in the offshore-onshore transition waters of the eastern North Pacific (Figure 4b,c), where large gradients in the buffer capacity occur (Jiang et al 2019, 2023). This indicates an expansion of the regions with low Ω_ar over time, and thus greater reduction of suitable pteropod habitat.

Polar/subpolar regions and eastern boundary upwelling systems are characterized by high pteropod biomass (Knecht et al., 2023). Rapid change in Ω_ar is evident near the frontal boundaries of both, so we may expect to see significant increases in dissolution in these areas over space and time. It is important to note that our projections concern the globally dominant pteropod species Limacina helicina. The few studies conducted so far that investigate the decline of CaCO₃ calcification at lower Ω_ar support the findings either across different regions (Lischka et al., 2011; Comeau et al., 2012; Bednaršek et al., 2014a) or in different species (Moya et al., 2016; Maas et al., 2018; Mekkes et al., 2021), indicating that dissolution estimates could be extrapolated among different pteropod species and regions. The rate of shell increase per decade might not seem substantial, but it is important to understand that this represents an additional percentage of the total number of pteropods that will exhibit severe shell dissolution. For example, a 4% increase per decade in polar regions indicates that an additional 12% of the total number of pteropods will exhibit severe shell dissolution by 2050, a surge from 50% to 62% of individuals affected, which could be substantial in population terms, especially in combination with reduced calcification. Shell dissolution projections are made for the early and mid-stages of the population, the demographics of which are critical to population dynamics (Bednaršek et al., 2016). However, current estimates of the relationship between shell dissolution and mortality are fairly uncertain given only one study on this (Bednaršek et al., 2017b). To date, no quantitative evaluation exists for a shell repair process (Peck et al., 2018; Niemi et al., 2021) that would offset extensive shell dissolution.

There is a much greater certainty that the population changes would accelerate in time in regions with overlapping OA and warming and marine heatwaves. Such stressor co-occurrence conditions may prove to be the most detrimental to pteropod populations (Lischka et al., 2011; Manno et al., 2017; Bednaršek et al., 2022; Elizondo and Vogt, 2022; Knecht at al., 2023). Co-linear deoxygenation and OA impacts (Feely et al., 2023, in this issue) could also have a relevant impact, but further study on this is needed.

Both acute subsurface shell dissolution and long-term population decline, which could have regionally important implications for the carbon system, are detectable through modulation in the subsurface ocean carbonate system, TA perturbations, and carbonate export (Feely et al., 2004; Sarmiento and Gruber, 2006; Bednaršek et al., 2012a; 2014a). Synthesizing observational data, biomass distribution, and dissolution, we predict substantial declines in carbon export and sequestration in the polar and subpolar regions of the North Pacific between 50° and 60°N, in the California and Humboldt Current Systems, and in the NEC and SEC regions, and strongly suggest targeted monitoring of carbon and sinking fluxes in these regions. Insights from this study are directly relevant for carbon sequestration and marine CO₂ removal strategies.

 

Acknowledgments

The authors would like to thank Nicolas Gruber and Matthias Münnich for their support with respect to the model simulation that was performed at the HPC cluster of ETH Zürich. We also want to express our most sincere appreciation to Dana Greeley and Julian Herndon for their help with pteropod sampling on the NOAA West Coast cruises. GP acknowledges Luke Gregor (ETH Zürich) for help in conceptualizing the method for creating a three-dimensional monthly climatology for aragonite saturation. Funding for RAF and NB was provided by the NOAA Ocean Acidification and Global Ocean Observing and Monitoring Programs and the Pacific Marine Environmental Laboratory (OAP Contract Number NRDD 20848 and GOMO Fund Reference Number 100018302). NB also acknowledges support from the Slovene Research Agency (ARRS “Biomarkers of subcellular stress in the Northern Adriatic under global environmental change,” project # J12468), as well as NOAA’s Multistressor project (project number NA22NOS4780171). This is PMEL contribution number 5495.

Bednaršek, N., R.A. Feely, G. Pelletier, and F. Desmet. 2023. Global synthesis of the status and trends of ocean acidification impacts on shelled pteropods. Oceanography 36(2–3):130–137, https://doi.org/10.5670/oceanog.2023.210.

Bakker, D.C., B. Pfeil, C.S. Landa, N. Metzl, K.M. O’Brien, A. Olsen, K. Smith, C. Cosca, S. Harasawa, S.D. Jones, and others. 2016. A multi-decade record of high-​quality fCO₂ data in version 3 of the Surface Ocean CO₂ Atlas (SOCAT). Earth System Science Data 8:383–413, https://doi.org/​10.5194/​essd-8-383-2016.
1. Bednaršek, N., J. Možina, M. Vogt, C. O’Brien, and G.A. Tarling. 2012a. The global distribution of pteropods and their contribution to carbonate and carbon biomass in the modern ocean. Earth System Science Data 4(1):167–186, https://doi.org/10.5194/essd-4-167-2012.
2. Bednaršek, N., G.A. Tarling, D.C.E. Bakker, S. Fielding, E.M. Jones, H.J. Venables, P. Ward, A. Kuzirian, B. Lézé, R.A. Feely, and E.J. Murphy. 2012b. Extensive dissolution of live pteropods in the Southern Ocean. Nature Geoscience 5(12):881–885, https://doi.org/​10.1038/​ngeo1635.
3. Bednaršek, N., G.A. Tarling, D.C.E. Bakker, S. Fielding, and R.A. Feely. 2014a. Dissolution dominating calcification process in polar pteropods close to the point of aragonite undersaturation. PLoS One 9(10):e109183, https://doi.org/10.1371/journal.pone.0109183.
4. Bednaršek, N., R.A. Feely, J.C.P. Reum, B. Peterson, J. Menkel, S.R. Alin, and B. Hales. 2014b. Limacina helicina shell dissolution as an indicator of declining habitat suitability owing to ocean acidification in the California Current Ecosystem. Proceedings of the Royal Society B: Biological Sciences 281(1785):20140123, https://doi.org/​10.1098/​rspb.2014.0123.
5. Bednaršek, N., C.J. Harvey, I.C. Kaplan, R.A. Feely, and J. Možina. 2016. Pteropods on the edge: Cumulative effects of ocean acidification, warming, and deoxygenation. Progress in Oceanography 145:1–24, https://doi.org/10.1016/​j.pocean.2016.04.002.
6. Bednaršek, N., R.A. Feely, N. Tolimieri, A.J. Hermann, S.A. Siedlecki, G.G. Waldbusser, P. McElhany, S.R. Alin, T. Klinger, B. Moore-Maley, and H.O. Pörtner. 2017a. Exposure history determines pteropod vulnerability to ocean acidification along the US West Coast. Scientific Reports 7(1):4526, https://doi.org/​10.1038/s41598-017-03934-z.
7. Bednaršek, N., T. Klinger, C.J. Harvey, S. Weisberg, R.M. McCabe, R.A. Feely, J. Newton, and N. Tolimieri. 2017b. New ocean, new needs: Application of pteropod shell dissolution as a biological indicator for marine resource management. Ecological Indicators 76:240–244, https://doi.org/​10.1016/​j.ecolind.2017.01.025.
8. Bednaršek, N., R.A. Feely, E.L. Howes, B.P.V. Hunt, F. Kessouri, P. León, S. Lischka, A.E. Maas, K. McLaughlin, N.P. Nezlin, and others. 2019. Systematic review and meta-analysis toward synthesis of thresholds of ocean acidification impacts on calcifying pteropods and interactions with warming. Frontiers in Marine Science 6:227, https://doi.org/10.3389/fmars.2019.00227.
9. Bednaršek, N., K.A. Naish, R.A. Feely, C. Hauri, K. Kimoto, A.J. Hermann, C. Michel, A. Niemi, and D. Pilcher. 2021. Integrated assessment of ocean acidification risks to pteropods in the northern high latitudes: Regional comparison of exposure, sensitivity and adaptive capacity. Frontiers in Marine Science 8:671497, https://doi.org/10.3389/fmars.2021.671497.
10. Bednaršek, N., B.R. Carter, R.M. McCabe, R.A. Feely, E. Howard, F.P. Chavez, M. Elliott, J.L. Fisher, J. Jahncke, and Z. Siegrist. 2022. Pelagic calcifiers face increased mortality and habitat loss with warming and ocean acidification. Ecological Applications 32(7):e2674, https://doi.org/10.1002/eap.2674.
11. Broullón, D., F.F. Pérez, A. Velo, M. Hoppema, A. Olsen, T. Takahashi, R.M. Key, T. Tanhua, M. González-Dávila, E. Jeansson, and others. 2019. A global monthly climatology of total alkalinity: A neural network approach. Earth System Science Data 11:1,109–1,127, https://doi.org/10.5194/essd-11-1109-2019.
12. Broullón, D., F.F. Pérez, A. Velo, M. Hoppema, A. Olsen, T. Takahashi, R.M. Key, T. Tanhua, J.M. Santana-Casiano, and A. Kozyr. 2020. A global monthly climatology of oceanic total dissolved inorganic carbon: A neural network approach. Earth System Science Data 12:1,725–1,743, https://doi.org/​10.5194/essd-12-1725-2020.
13. Buitenhuis, E.T., C. Le Quere, N. Bednaršek, and R. Schiebel. 2019. Large contribution of pteropods to shallow CaCO₃ export. Global Biogeochemical Cycles 33(3):458–468, https://doi.org/​10.1029/​2018GB006110.
14. Comeau, S., G. Gorsky, R. Jeffree, J.L. Teyssié, and J.P. Gattuso. 2009. Impact of ocean acidification on a key Arctic pelagic mollusc (Limacina helicina). Biogeosciences 6(9):1,877–1,882, https://doi.org/​10.5194/bg-6-1877-2009.
15. Comeau, S., S. Alliouane, and J.P. Gattuso. 2012. Effects of ocean acidification on overwintering juvenile Arctic pteropods Limacina helicina. Marine Ecology Progress Series 456:279–284, https://doi.org/​10.3354/meps09696.
16. Desmet, F., N. Gruber, E.E. Köhn, M. Münnich, and M. Vogt. 2022. Tracking the space-time evolution of ocean acidification extremes in the California Current System and Northeast Pacific. Journal of Geophysical Research: Oceans 127:e2021JC018159, https://doi.org/​10.1029/​2021JC018159.
17. Dickson, A.G., C.L. Sabine, and J.R. Christian. 2007. Guide to best practices for ocean CO₂ measurements. PICES Special Publication 3(8):191.
18. Elizondo, U.H., and M. Vogt. 2022. Individual-based modeling of shelled pteropods. Ecological Modelling 468:109944, https://doi.org/10.1016/​j.ecolmodel.2022.109944.
19. Fabry, V.J., B.A. Seibel, R.A. Feely, and J.C. Orr. 2008. Impacts of ocean acidification on marine fauna and ecosystem processes. ICES Journal of Marine Science 65(3):414–432, https://doi.org/10.1093/icesjms/fsn048.
20. Feely, R.A., C.L. Sabine, K. Lee, F.J. Millero, M.F. Lamb, D. Greeley, J.L. Bullister, R.M. Key, T.-H. Peng, A. Kozyr, T. Ono, and C.S. Wong, 2002. In situ calcium carbonate dissolution in the Pacific Ocean. Global Biogeochemical Cycles 16(4):91-1–91-12, https://doi.org/​10.1029/2002GB001866.
21. Feely, R.A., C.L. Sabine, K. Lee, W. Berelson, J. Kleypas, V.J. Fabry, and F.J. Millero. 2004. Impact of anthropogenic CO₂ on the CaCO₃ system in the oceans. Science 305(5682):362–366, https://doi.org/​10.1126/science.1097329.
22. Feely, R.A., S.R. Alin, B. Carter, N. Bednaršek, B. Hales, F. Chan, T.M. Hill, B. Gaylord, E. Sanford, R.H. Byrne, and C.L. Sabine. 2016. Chemical and biological impacts of ocean acidification along the west coast of North America. Estuarine, Coastal and Shelf Science 183:260–270, https://doi.org/​10.1016/​j.ecss.2016.08.043.
23. Feely, R.A., L.-Q. Jiang, R. Wanninkhof, B.R. Carter, S.R. Alin, N. Bednaršek, and C.E. Cosca. 2023. Acidification of the global surface ocean: What we have learned from observations. Oceanography 36(2–3):120–129, https://doi.org/​10.5670/oceanog.2023.222.
24. Friedlingstein, P., M. O’Sullivan, M.W. Jones, R.M. Andrew, J. Hauck, A. Olsen, G.P. Peters, W. Peters, J. Pongratz, S. Sitch, and others. 2020. Global carbon budget 2020. Earth System Science Data 12(4):3,269–3,340, https://doi.org/10.5194/essd-12-3269-2020.
25. Friedlingstein, P., M.W. Jones, M. O’Sullivan, R.M. Andrew, D.C.E. Bakker, J. Hauck, C. Le Quéré, G.P. Peters, W. Peters, J. Pongratz, and others. 2022. Global carbon budget 2021. Earth System Science Data 14(4):1,917–2,005, https://doi.org/​10.5194/​essd-14-1917-2022.
26. Frischknecht, M., M. Münnich, and N. Gruber. 2018. Origin, transformation, and fate: The three-​dimensional biological pump in the California Current System. Journal of Geophysical Research: Oceans 123(11):7,939–7,962, https://doi.org/​10.1029/​2018JC013934.
27. Gehlen, M., R. Gangstø, B. Schneider, L. Bopp, O. Aumont, and C. Éthé. 2007. The fate of pelagic CaCO₃ production in a high CO₂ ocean: A model study. Biogeosciences 4(4):505–519, https://doi.org/​10.5194/bg-4-505-2007.
28. Gregor, L., and N. Gruber. 2021. OceanSODA-ETHZ: A global gridded data set of the surface ocean carbonate system for seasonal to decadal studies of ocean acidification. Earth System Science Data 13:777–808, https://doi.org/10.5194/essd-13-777-2021.
29. Honjo, S., S.J. Manganini, R.A. Krishfield, and R. Francois. 2008. Particulate organic carbon fluxes to the ocean interior and factors controlling the biological pump: A synthesis of global sediment trap programs since 1983. Progress in Oceanography 76(3):217–285, https://doi.org/​10.1016/j.pocean.2007.11.003.
30. IPCC (Intergovernmental Panel on Climate Change). 2022. Climate Change 2022: Impacts, Adaptation, and Vulnerability. Contribution of Working Group II to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change. H.-O. Pörtner, D.C. Roberts, M. Tignor, E.S. Poloczanska, K. Mintenbeck, A. Alegría, M. Craig, S. Langsdorf, S. Löschke, V. Möller, A. Okem, B. Rama, eds, Cambridge University Press, Cambridge University Press, Cambridge, UK, and New York, NY, USA, 3,056 pp.
31. Jiang, L.-Q., B.R. Carter, R.A. Feely, S.K. Lauvset, and A. Olsen. 2019. Surface ocean pH and buffer capacity: Past, present and future. Scientific Reports 9:18624, https://doi.org/10.1038/s41598-019-55039-4.
32. Jiang, L.-Q., J. Dunne, B.R. Carter, J.F. Tjiputra, J. Terhaar, J.D. Sharp, A. Olsen, S. Alin, D.C.E. Bakker, R.A. Feely, P. Hogan, and others. 2023. Global surface ocean acidification indicators from 1750 to 2100. Journal of Advances in Modeling Earth Systems 15(3):e2022MS003563, https://doi.org/​10.1029/2022MS003563.
33. Knecht, N.S., F. Benedetti, U.H. Elizondo, N. Bednaršek, S. Chaabane, C. de Weerd, K.T.C.A. Peijnenburg, R. Schiebel, and M. Vogt. 2023. The impact of zooplankton calcifiers on the marine carbon cycle. Global Biogeochemical Cycles 37:e2022GB007685, https://doi.org/​10.1029/2022GB007685.
34. Lee, K., and R.A. Feely. 2021. Dissolution resolution. Nature Geosciences 14:356–358, https://doi.org/​10.1038/s41561-021-00765-6.
35. Lischka, S., J. Büdenbender, T. Boxhammer, and U. Riebesell. 2011. 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(4):919–932, https://doi.org/​10.5194/bg-8-919-2011.
36. Locarnini, R.A., A.V. Mishonov, O.K. Baranova, T.P. Boyer, M.M. Zweng, H.E. Garcia, J.R. Reagan, D. Seidov, K. Weathers, C.R. Paver, and I. Smolyar. 2018. World Ocean Atlas 2018, Volume 1: Temperature. A. Mishonov, ed., NOAA Atlas NESDIS 81, 52 pp.
37. Manno, C., N. Bednaršek, G.A. Tarling, V.L. Peck, S. Comeau, D. Adhikari, D.C. Bakker, E. Bauerfeind, A.J. Bergan, M.I. Berning, and E. Buitenhuis. 2017. Shelled pteropods in peril: Assessing vulnerability in a high CO₂ ocean. Earth-Science Reviews 169:132–145, https://doi.org/10.1016/​j.earscirev.2017.04.005.
38. Marchesiello, P., J.C. McWilliams, and A. Shchepetkin. 2003. Equilibrium structure and dynamics of the California Current System. Journal of Physical Oceanography 33(4):753–783, https://doi.org/​10.1175/1520-0485(2003)33<753:esadot>2.0.co;2.
39. Mekkes, L., G. Sepúlveda-Rodríguez, G. Bielkinaitė, D. Wall-Palmer, G.J.A. Brummer, L.K. Dämmer, J. Huisman, E. van Loon, W. Renema, and K.T.C.A. Peijnenburg. 2021. Effects of ocean acidification on calcification of the sub-​Antarctic pteropod Limacina retroversa. Frontiers in Marine Science 8:581432, https://doi.org/10.3389/fmars.​2021.581432.
40. Milliman, J.D. 1993. Production and accumulation of calcium carbonate in the ocean: Budget of a nonsteady state. Global Biogeochemical Cycles 7(4):927–957, https://doi.org/​10.1029/​93GB02524.
41. Moore, J.K., K. Lindsay, S.C. Doney, M.C. Long, and K. Misumi. 2013. Marine ecosystem dynamics and biogeochemical cycling in the Community Earth system model [CESM1(BGC)]: Comparison of the 1990s with the 2090s under the RCP4.5 and RCP8.5 scenarios. Journal of Climate 26:9,291–9,312, https://doi.org/10.1175/JCLI-D-12-00566.1.
42. Mucci, A. 1983. The solubility of calcite and aragonite in seawater at various salinities, temperatures, and one atmosphere total pressure. American Journal of Science 283(7):780–799, https://doi.org/10.2475/ajs.283.7.780.
43. Niemi, A., N. Bednaršek, C. Michel, R.A. Feely, W. Williams, K. Azetsu-Scott, W. Walkusz, and J.D. Reist. 2021. Biological impact of ocean acidification in the Canadian Arctic: Widespread severe pteropod shell dissolution in Amundsen Gulf. Frontiers in Marine Science 8:600184, https://doi.org/​10.3389/​fmars.2021.600184.
44. Olsen, A., N. Lange, R.M. Key, T. Tanhua, M. Álvarez, S. Becker, H.C. Bittig, B.R. Carter, L. Cotrim Da Cunha, R.A. Feely, and others. 2019. GLODAPv2.2019 – An update of GLODAPv2. Earth System Science Data 11:1,437–1,461, https://doi.org/​10.5194/​essd-11-1437-2019.
45. Orr, J.C., and J.M. Epitalon. 2015. Improved routines to model the ocean carbonate system: Mocsy 2.0. Geoscientific Model Development 8(3):485–499, https://doi.org/10.5194/gmd-8-485-2015.
46. Peck, V.L., R.L. Oakes, E.M. Harper, C. Manno, and G.A. Tarling, 2018. Pteropods counter mechanical damage and dissolution through extensive shell repair. Nature Communications 9(1):264, https://doi.org/​10.1038/s41467-017-02692-w.
47. Sarmiento, J.L., and N. Gruber. 2006. Ocean Biogeochemical Dynamics. Princeton University Press, 528 pp.
48. Shchepetkin, A.F., and J.C. McWilliams. 2005. The regional oceanic modeling system (ROMS): A split-​explicit, free-surface, topography-​following-​​coordinate oceanic model. Ocean Modelling 9(4):347–404, https://doi.org/10.1016/​j.ocemod.2004.08.002.
49. Zeebe, R.E., and D. Wolf-Gladrow. 2001. CO₂ in Seawater: Equilibrium, Kinetics, Isotopes. Elsevier Oceanography Series, vol. 65, 360 pp.
50. Zweng, M.M., J.R. Reagan, D. Seidov, T.P. Boyer, R.A. Locarnini, H.E. Garcia, A.V. Mishonov, O.K. Baranova, K. Weathers, C.R. Paver, and I. Smolyar. 2018. World Ocean Atlas 2018, Volume 2: Salinity. A. Mishonov, ed., NOAA Atlas NESDIS 82, 50 pp.

This is an open access article made available under the terms of the Creative Commons Attribution 4.0 International License (https://creativecommons.org/licenses/by/4.0/), which permits use, sharing, adaptation, distribution, and reproduction in any medium or format as long as users cite the materials appropriately, provide a link to the Creative Commons license, and indicate the changes that were made to the original content. Images, animations, videos, or other third-party material used in articles are included in the Creative Commons license unless indicated otherwise in a credit line to the material. If the material is not included in the article’s Creative Commons license, users will need to obtain permission directly from the license holder to reproduce the material.