The impact of aragonite saturation variability on shelled pteropods: An attribution study in the California Current System.

Observations from the California Current System (CalCS) indicate that the long-term trend in ocean acidification (OA) and the naturally occurring corrosive conditions for the CaCO3 mineral aragonite (saturation state Ω < 1) have a damaging effect on shelled pteropods, a keystone group of calcifying organisms in the CalCS. Concern is heightened by recent findings suggesting that shell formation and developmental progress are already impacted when Ω falls below 1.5. Here, we quantify the impact of low Ω conditions on pteropods using an individual-based model (IBM) with life-stage-specific mortality, growth, and behavior in a high-resolution regional hindcast simulation of the CalCS between 1984 and 2019. Special attention is paid to attributing this impact to different processes that lead to such low Ω conditions, namely natural variability, long-term trend, and extreme events. We find that much of the observed damage in the CalCS, and specifically >70% of the shell CaCO3 loss, is due to the pteropods' exposure to naturally occurring low Ω conditions as a result of their diel vertical migration (DVM). Over the hindcast period, their exposure to damaging waters (Ω < 1.5) increases from 9% to 49%, doubling their shell CaCO3 loss, and increasing their mortality by ~40%. Most of this increased exposure is due to the shoaling of low Ω waters driven by the long-term trend in OA. Extreme OA events amplify this increase by ~40%. Our approach can quantify the health of pteropod populations under shifting environmental conditions, and attribute changes in fitness or population structure to changes in the stressor landscape across hierarchical time scales.

Coastal eutrophication drives acidification, oxygen loss, and ecosystem change in a major oceanic upwelling system.

Global change is leading to warming, acidification, and oxygen loss in the ocean. In the Southern California Bight, an eastern boundary upwelling system, these stressors are exacerbated by the localized discharge of anthropogenically enhanced nutrients from a coastal population of 23 million people. Here, we use simulations with a high-resolution, physical-biogeochemical model to quantify the link between terrestrial and atmospheric nutrients, organic matter, and carbon inputs and biogeochemical change in the coastal waters of the Southern California Bight. The model is forced by large-scale climatic drivers and a reconstruction of local inputs via rivers, wastewater outfalls, and atmospheric deposition; it captures the fine scales of ocean circulation along the shelf; and it is validated against a large collection of physical and biogeochemical observations. Local land-based and atmospheric inputs, enhanced by anthropogenic sources, drive a 79% increase in phytoplankton biomass, a 23% increase in primary production, and a nearly 44% increase in subsurface respiration rates along the coast in summer, reshaping the biogeochemistry of the Southern California Bight. Seasonal reductions in subsurface oxygen, pH, and aragonite saturation state, by up to 50 mmol m-3, 0.09, and 0.47, respectively, rival or exceed the global open-ocean oxygen loss and acidification since the preindustrial period. The biological effects of these changes on local fisheries, proliferation of harmful algal blooms, water clarity, and submerged aquatic vegetation have yet to be fully explored.

Pelagic calcifiers face increased mortality and habitat loss with warming and ocean acidification.

Global change is impacting the oceans in an unprecedented way, and multiple lines of evidence suggest that species distributions are changing in space and time. There is increasing evidence that multiple environmental stressors act together to constrain species habitat more than expected from warming alone. Here, we conducted a comprehensive study of how temperature and aragonite saturation state act together to limit Limacina helicina, globally distributed pteropods that are ecologically important pelagic calcifiers and an indicator species for ocean change. We co-validated three different approaches to evaluate the impact of ocean warming and acidification (OWA) on the survival and distribution of this species in the California Current Ecosystem. First, we used colocated physical, chemical, and biological data from three large-scale west coast cruises and regional time series; second, we conducted multifactorial experimental incubations to evaluate how OWA impacts pteropod survival; and third, we validated the relationships we found against global distributions of pteropods and carbonate chemistry. OWA experimental work revealed mortality increases under OWA, while regional habitat suitability indices and global distributions of L. helicina suggest that a multi-stressor framework is essential for understanding pteropod distributions. In California Current Ecosystem habitats, where pteropods are living close to their thermal maximum already, additional warming and acidification through unabated fossil fuel emissions (RCP 8.5) are expected to dramatically reduce habitat suitability.

Natural Analogues in pH Variability and Predictability across the Coastal Pacific Estuaries: Extrapolation of the Increased Oyster Dissolution under Increased pH Amplitude and Low Predictability Related to Ocean Acidification.

Coastal-estuarine habitats are rapidly changing due to global climate change, with impacts influenced by the variability of carbonate chemistry conditions. However, our understanding of the responses of ecologically and economically important calcifiers to pH variability and temporal variation is limited, particularly with respect to shell-building processes. We investigated the mechanisms driving biomineralogical and physiological responses in juveniles of introduced (Pacific; Crassostrea gigas) and native (Olympia; Ostrea lurida) oysters under flow-through experimental conditions over a six-week period that simulate current and future conditions: static control and low pH (8.0 and 7.7); low pH with fluctuating (24-h) amplitude (7.7 ± 0.2 and 7.7 ± 0.5); and high-frequency (12-h) fluctuating (8.0 ± 0.2) treatment. The oysters showed physiological tolerance in vital processes, including calcification, respiration, clearance, and survival. However, shell dissolution significantly increased with larger amplitudes of pH variability compared to static pH conditions, attributable to the longer cumulative exposure to lower pH conditions, with the dissolution threshold of pH 7.7 with 0.2 amplitude. Moreover, the high-frequency treatment triggered significantly greater dissolution, likely because of the oyster's inability to respond to the unpredictable frequency of variations. The experimental findings were extrapolated to provide context for conditions existing in several Pacific coastal estuaries, with time series analyses demonstrating unique signatures of pH predictability and variability in these habitats, indicating potentially benefiting effects on fitness in these habitats. These implications are crucial for evaluating the suitability of coastal habitats for aquaculture, adaptation, and carbon dioxide removal strategies.

Biogeography and genetic diversity of the atlantid heteropods.

The atlantid heteropods are regularly encountered, but rarely studied marine planktonic gastropods. Relying on a small (<14 mm), delicate aragonite shell and living in the upper ocean means that, in common with pteropods, atlantids are likely to be affected by imminent ocean changes. Variable shell morphology and widespread distributions indicate that the family is more diverse than the 23 currently known species. Uncovering this diversity is fundamental to determining the distribution of atlantids and to understanding their environmental tolerances. Here we present phylogenetic analyses of all described species of the family Atlantidae using 437 new and 52 previously published cytochrome c oxidase subunit 1 mitochondrial DNA (mtCO1) sequences. Specimens and published sequences were gathered from 32 Atlantic Ocean stations, 14 Indian Ocean stations and 21 Pacific Ocean stations between 35°N and 43°S. DNA barcoding and Automatic Barcode Gap Discovery (ABGD) proved to be valuable tools for the identification of described atlantid species, and also revealed ten additional distinct clades, suggesting that the diversity within this family has been underestimated. Only two of these clades displayed obvious morphological characteristics, demonstrating that much of the newly discovered diversity is hidden from morphology-based identification techniques. Investigation of six large atlantid collections demonstrated that 61% of previously described (morpho) species have a circumglobal distribution. Of the remaining 39%, two species were restricted to the Atlantic Ocean, five occurred in the Indian and Pacific oceans, one species was only found in the northeast Pacific Ocean, and one occurred only in the Southern Subtropical Convergence Zone. Molecular analysis showed that seven of the species with wide distributions were comprised of two or more clades that occupied distinct oceanographic regions. These distributions may suggest narrower environmental tolerances than the described morphospecies. Results provide an updated biogeography and mtCO1 reference dataset of the Atlantidae that may be used to identify atlantid species and provide a first step in understanding their evolutionary history and accurate distribution, encouraging the inclusion of this family in future plankton research.

Predictable patterns within the kelp forest can indirectly create temporary refugia from ocean acidification.

Kelps are recognized for providing many ecosystem services in coastal areas and considered in ocean acidification (OA) mitigation. However, assessing OA modification requires an understanding of the multiple parameters involved in carbonate chemistry, especially in highly dynamic systems. We studied the effects of sugar kelp (Saccharina latissima) on an experimental farm at the north end of Hood Canal, Washington-a low retentive coastal system. In this field mesocosm study, two oyster species (Magallana gigas, Ostrea lurida) were exposed at locations in the mid, edge, and outside the kelp array. The Hood Head Sugar Kelp Farm Model outputs were used to identify dominating factors in spatial and temporal kelp dynamics, while wavelet spectrum analyses helped in understanding predictability patterns. This was linked to the measured biological responses (dissolution, growth, isotopes) of the exposed organisms. Positioned in an area of high (sub)-diel tidal fluxes with low retention potential, there were no measurable alterations of the seawater pH at the study site, demonstrating that the kelp array could not induce a direct mitigating effect against OA. However, beneficial responses in calcifiers were still observed, which are linked to two causes: increased pH predictability and improved provisioning through kelp-derived particulate organic resource utilization and as such, kelp improved habitat suitability and indirectly created refugia against OA. This study can serve as an analogue for many coastal bay habitats where prevailing physical forcing drives chemical changes. Future macrophyte studies that investigate OA mitigating effects should focus also on the importance of predictability patterns, which can additionally improve the conditions for marine calcifiers and ecosystem services vulnerable to or compromised by OA, including aquaculture sustainability.

Pteropods make thinner shells in the upwelling region of the California Current Ecosystem.

Shelled pteropods are widely regarded as bioindicators for ocean acidification, because their fragile aragonite shells are susceptible to increasing ocean acidity. While short-term incubations have demonstrated that pteropod calcification is negatively impacted by ocean acidification, we know little about net calcification in response to varying ocean conditions in natural populations. Here, we examine in situ calcification of Limacina helicina pteropods collected from the California Current Ecosystem, a coastal upwelling system with strong spatial gradients in ocean carbonate chemistry, dissolved oxygen and temperature. Depth-averaged pH ranged from 8.03 in warmer offshore waters to 7.77 in cold CO2-rich waters nearshore. Based on high-resolution micro-CT technology, we showed that shell thickness declined by ~ 37% along the upwelling gradient from offshore to nearshore water. Dissolution marks covered only ~ 2% of the shell surface area and were not associated with the observed variation in shell thickness. We thus infer that pteropods make thinner shells where upwelling brings more acidified and colder waters to the surface. Probably the thinner shells do not result from enhanced dissolution, but are due to a decline in calcification. Reduced calcification of pteropods is likely to have major ecological and biogeochemical implications for the cycling of calcium carbonate in the oceans.

Severe biological effects under present-day estuarine acidification in the seasonally variable Salish Sea.

Estuaries are recognized as one of the habitats most vulnerable to coastal ocean acidification due to seasonal extremes and prolonged duration of acidified conditions. This is combined with co-occurring environmental stressors such as increased temperature and low dissolved oxygen. Despite this, evidence of biological impacts of ocean acidification in estuarine habitats is largely lacking. By combining physical, biogeochemical, and biological time-series observations over relevant seasonal-to-interannual time scales, this study is the first to describe both the spatial and temporal variation of biological response in the pteropod Limacina helicina to estuarine acidification in association with other stressors. Using clustering and principal component analyses, sampling sites were grouped according to their distribution of physical and biogeochemical variables over space and time. This identified the most exposed habitats and time intervals corresponding to the most severe negative biological impacts across three seasons and three years. We developed a cumulative stress index as a means of integrating spatial-temporal OA variation over the organismal life history. Our findings show that over the 2014-2016 study period, the severity of low aragonite saturation state combined with the duration of exposure contributed to overall cumulative stress and resulted in severe shell dissolution. Seasonally-variable estuaries such as the Salish Sea (Washington, U.S.A.) predispose sensitive organisms to more severe acidified conditions than those of coastal and open-ocean habitats, yet the sensitive organisms persist. We suggest potential environmental factors and compensatory mechanisms that allow pelagic calcifiers to inhabit less favorable habitats and partially offset associated stressors, for instance through food supply, increased temperature, and adaptation of their life history. The novel metric of cumulative stress developed here can be applied to other estuarine environments with similar physical and chemical dynamics, providing a new tool for monitoring biological response in estuaries under pressure from accelerating global change.

Natural and Anthropogenic Drivers of Acidification in Large Estuaries.

Oceanic uptake of anthropogenic carbon dioxide (CO2) from the atmosphere has changed ocean biogeochemistry and threatened the health of organisms through a process known as ocean acidification (OA). Such large-scale changes affect ecosystem functions and can have impacts on societal uses, fisheries resources, and economies. In many large estuaries, anthropogenic CO2-induced acidification is enhanced by strong stratification, long water residence times, eutrophication, and a weak acid-base buffer capacity. In this article, we review how a variety of processes influence aquatic acid-base properties in estuarine waters, including coastal upwelling, river-ocean mixing, air-water gas exchange, biological production and subsequent aerobic and anaerobic respiration, calcium carbonate (CaCO3) dissolution, and benthic inputs. We emphasize the spatial and temporal dynamics of partial pressure of CO2 (pCO2), pH, and calcium carbonate mineral saturation states. Examples from three large estuaries-Chesapeake Bay, the Salish Sea, and Prince William Sound-are used to illustrate how natural and anthropogenic processes and climate change may manifest differently across estuaries, as well as the biological implications of OA on coastal calcifiers.

Exoskeleton dissolution with mechanoreceptor damage in larval Dungeness crab related to severity of present-day ocean acidification vertical gradients.

Ocean acidification (OA) along the US West Coast is intensifying faster than observed in the global ocean. This is particularly true in nearshore regions (<200 m) that experience a lower buffering capacity while at the same time providing important habitats for ecologically and economically significant species. While the literature on the effects of OA from laboratory experiments is voluminous, there is little understanding of present-day OA in-situ effects on marine life. Dungeness crab (Metacarcinus magister) is perennially one of the most valuable commercial and recreational fisheries. We focused on establishing OA-related vulnerability of larval crustacean based on mineralogical and elemental carapace to external and internal carapace dissolution by using a combination of different methods ranging from scanning electron microscopy, energy dispersive X-ray spectroscopy, elemental mapping and X-ray diffraction. By integrating carapace features with the chemical observations and biogeochemical model hindcast, we identify the occurrence of external carapace dissolution related to the steepest Ω calcite gradients (∆Ωcal,60) in the water column. Dissolution features are observed across the carapace, pereopods (legs), and around the calcified areas surrounding neuritic canals of mechanoreceptors. The carapace dissolution is the most extensive in the coastal habitats under prolonged (1-month) long exposure, as demonstrated by the use of the model hindcast. Such dissolution has a potential to destabilize mechanoreceptors with important sensory and behavioral functions, a pathway of sensitivity to OA. Carapace dissolution is negatively related to crab larval width, demonstrating a basis for energetic trade-offs. Using a retrospective prediction from a regression models, we estimate an 8.3% increase in external carapace dissolution over the last two decades and identified a set of affected OA-related sublethal pathways to inform future risk assessment studies of Dungeness crabs.

Towards Integrating Evolution, Metabolism, and Climate Change Studies of Marine Ecosystems.

Global environmental changes are challenging the structure and functioning of ecosystems. However, a mechanistic understanding of how global environmental changes will affect ecosystems is still lacking. The complex and interacting biological and physical processes spanning vast temporal and spatial scales that constitute an ecosystem make this a formidable problem. A unifying framework based on ecological theory, that considers fundamental and realized niches, combined with metabolic, evolutionary, and climate change studies, is needed to provide the mechanistic understanding required to evaluate and forecast the future of marine communities, ecosystems, and their services.

Dissolution dominating calcification process in polar pteropods close to the point of aragonite undersaturation.

Thecosome pteropods are abundant upper-ocean zooplankton that build aragonite shells. Ocean acidification results in the lowering of aragonite saturation levels in the surface layers, and several incubation studies have shown that rates of calcification in these organisms decrease as a result. This study provides a weight-specific net calcification rate function for thecosome pteropods that includes both rates of dissolution and calcification over a range of plausible future aragonite saturation states (Ω(ar)). We measured gross dissolution in the pteropod Limacina helicina antarctica in the Scotia Sea (Southern Ocean) by incubating living specimens across a range of aragonite saturation states for a maximum of 14 days. Specimens started dissolving almost immediately upon exposure to undersaturated conditions (Ω(ar) ∼ 0.8), losing 1.4% of shell mass per day. The observed rate of gross dissolution was different from that predicted by rate law kinetics of aragonite dissolution, in being higher at Ω(ar) levels slightly above 1 and lower at Ω(ar) levels of between 1 and 0.8. This indicates that shell mass is affected by even transitional levels of saturation, but there is, nevertheless, some partial means of protection for shells when in undersaturated conditions. A function for gross dissolution against Ω(ar) derived from the present observations was compared to a function for gross calcification derived by a different study, and showed that dissolution became the dominating process even at Ω(ar) levels close to 1, with net shell growth ceasing at an Ω(ar) of 1.03. Gross dissolution increasingly dominated net change in shell mass as saturation levels decreased below 1. As well as influencing their viability, such dissolution of pteropod shells in the surface layers will result in slower sinking velocities and decreased carbon and carbonate fluxes to the deep ocean.

Poles apart: the "bipolar" pteropod species Limacina helicina is genetically distinct between the Arctic and Antarctic oceans.

The shelled pteropod (sea butterfly) Limacina helicina is currently recognised as a species complex comprising two sub-species and at least five "forma". However, at the species level it is considered to be bipolar, occurring in both the Arctic and Antarctic oceans. Due to its aragonite shell and polar distribution L. helicina is particularly vulnerable to ocean acidification. As a key indicator of the acidification process, and a major component of polar ecosystems, L. helicina has become a focus for acidification research. New observations that taxonomic groups may respond quite differently to acidification prompted us to reassess the taxonomic status of this important species. We found a 33.56% (+/-0.09) difference in cytochrome c oxidase subunit I (COI) gene sequences between L. helicina collected from the Arctic and Antarctic oceans. This degree of separation is sufficient for ordinal level taxonomic separation in other organisms and provides strong evidence for the Arctic and Antarctic populations of L. helicina differing at least at the species level. Recent research has highlighted substantial physiological differences between the poles for another supposedly bipolar pteropod species, Clione limacina. Given the large genetic divergence between Arctic and Antarctic L. helicina populations shown here, similarly large physiological differences may exist between the poles for the L. helicina species group. Therefore, in addition to indicating that L. helicina is in fact not bipolar, our study demonstrates the need for acidification research to take into account the possibility that the L. helicina species group may not respond in the same way to ocean acidification in Arctic and Antarctic ecosystems.