Ocean acidification and changes in biological production in the western subarctic region of the North Pacific over the quarter century, 1999–2023

Following the Last Glacial Maximum (LGM; ca. 23-19 calibrated [cal.] kyr before present [BP]), atmospheric and oceanic warming, together with global sea-level rise, drove widespread deglaciation of the Antarctic Ice Sheet, increasing the flux of freshwater to the ocean and leading to substantial changes in marine biological productivity. On the Antarctic continental shelf, periods of elevated biological productivity, often preserved in the sediment record as laminated (and sometimes varved) diatomaceous oozes (LDO), have been reported from several locations and are typically associated with the formation of calving bay re-entrants during ice sheet retreat. Understanding what drives the formation and deposition of LDOs, and the impact of deglacial processes on biogenic productivity more generally, can help inform how Antarctic coastal environments will respond to current and future ice sheet melting. In this study we utilise a suite of sediment cores recovered from Anvers-Hugo Trough (AHT), western Antarctic Peninsula shelf, which documents the transition from subglacial to glacimarine conditions following retreat of an expanded ice stream after the LGM. We present quantitative absolute diatom abundance (ADA) and species assemblage data, to investigate changes in biological productivity during the Last Glacial Transition (19-11 cal kyr BP). In combination with radiocarbon dating, we show that seasonally open marine conditions were established on the mid-shelf by 13.6 cal kyr BP, but LDOs did not start to accumulate until ∼11.5 cal kyr BP. The ∼1.4 kyr delay between the onset of seasonally open marine conditions and LDO deposition indicates that physiographic changes, and specifically the establishment of a calving bay in AHT, is insufficient to explain LDO deposition alone. LDO deposition in AHT coincides with the early Holocene climatic optimum (∼11.5 – 9.0 kyr) and is therefore explained in terms of increased atmospheric/ocean temperatures, high rates of sea and glacial ice melt and the formation of a well-stratified water column in the austral spring. An implication of our study is that extensive bathymetric mapping in conjunction with detailed core analyses is required to reliably infer environmental controls on LDO deposition.

Changes in biological productivity and ocean-climatic fluctuations during the last ~ 1.5 kyr in the Humboldt ecosystem off northern Chile (27°S): A multiproxy approach

Abstract. Southern Ocean organic carbon export plays an important role in the global carbon cycle, yet its basin-scale climatology and variability are uncertain due to limited coverage of in situ observations. In this study, a neural network approach based on the self-organizing map (SOM) is adopted to construct weekly gridded (1° × 1°) maps of organic carbon export for the Southern Ocean from 1998 to 2009. The SOM is trained with in situ measurements of O2 / Ar-derived net community production (NCP) that are tightly linked to the carbon export in the mixed layer on timescales of one to two weeks and with six potential NCP predictors: photosynthetically available radiation (PAR), particulate organic carbon (POC), chlorophyll (Chl), sea surface temperature (SST), sea surface height (SSH), and mixed layer depth (MLD). This nonparametric approach is based entirely on the observed statistical relationships between NCP and the predictors and, therefore, is strongly constrained by observations. A thorough cross-validation yields three retained NCP predictors, Chl, PAR, and MLD. Our constructed NCP is further validated by good agreement with previously published, independent in situ derived NCP of weekly or longer temporal resolution through real-time and climatological comparisons at various sampling sites. The resulting November–March NCP climatology reveals a pronounced zonal band of high NCP roughly following the Subtropical Front in the Atlantic, Indian, and western Pacific sectors, and turns southeastward shortly after the dateline. Other regions of elevated NCP include the upwelling zones off Chile and Namibia, the Patagonian Shelf, the Antarctic coast, and areas surrounding the Islands of Kerguelen, South Georgia, and Crozet. This basin-scale NCP climatology closely resembles that of the satellite POC field and observed air–sea CO2 flux. The long-term mean area-integrated NCP south of 50° S from our dataset, 17.9 mmol C m−2 d−1, falls within the range of 8.3 to 24 mmol C m−2 d−1 from other model estimates. A broad agreement is found in the basin-wide NCP climatology among various models but with significant spatial variations, particularly in the Patagonian Shelf. Our approach provides a comprehensive view of the Southern Ocean NCP climatology and a potential opportunity to further investigate interannual and intraseasonal variability.

The contribution of the marine biota to air‐sea fluxes of CO2 and O2 is often described in terms of biological production concepts, such as new production, export production, and net community production. We evaluate these three quantities using a basin‐scale ecosystem‐circulation model of the North Atlantic Ocean based on Redfield stoichiometry into which we introduce an artificial tracer which records the biotic contribution to air‐sea exchange of gases like O2 and CO2. It is found that on average the biological production rates overestimate the biotically effected air‐sea flux by some 20% and, in some regions, even predict the wrong direction. With primary production restricted to the euphotic zone, but respiration extending to farther below, the discrepancy can largely be attributed to the different integration depths used in the different concepts (euphotic zone, surface mixed layer), and on annual and longer timescales, all rates converge when using the base of the winter mixed layer rather than that of the euphotic zone as the reference depth. For the surface carbon budget, which ultimately controls air‐sea exchange of CO2, it is irrelevant whether carbon atoms cross this boundary in organic or inorganic speciation. Hence the transports of biotically generated surpluses or deficits of dissolved inorganic matter must also be accounted for. While their contribution amounts to only a few percent on the basin scale, the subduction of newly remineralized inorganic matter can locally account for about half of the biotically effected air‐sea flux, for example, in regions of mode‐water formation.

The role of sea ice formation in cycling of aluminium in northern Marguerite Bay, Antarctica

We investigate the sensitivity of the oxygen content and true oxygen utilization of key low‐oxygen regions Ω to pointwise changes in biological production. To understand how the combined water and biogenic particle transport controls the sensitivity patterns and the fate of oxygen in the ocean, we develop new relationships that link the steady‐state oxygen content and deficit of Ω to the downstream and upstream oxygen utilization rate (OUR), respectively. We find that the amount of oxygen from Ω that will be lost per unit volume at point r is linked to OUR(r) through the mean oxygen age accumulated in Ω. The geographic sensitivity pattern of the Ω‐integrated oxygen deficit is shaped by where the utilization occurs that causes this deficit. The contribution to the oxygen deficit of Ω from utilization at r is controlled by the mean time that water at r spends in Ω before next ventilation at the surface. We illustrate these relationships and the new transport timescales using a simple steady‐state data‐constrained carbon and oxygen model. We focus on Ω being the global ocean, the Pacific Hypoxic Zone (PHZ, [O2] < 62.5 µM), and the North Pacific oxygen minimum zone. The oxygen deficit of the PHZ is most sensitive where mode and intermediate waters form and where increased organic‐matter production directly increases the PHZ's oxygen demand. The fraction of the local oxygen concentration that will be utilized in respiration is as high as 90% in the PHZ and up to 70% in the water column beneath it.

Rates of summertime biological productivity in the Beaufort Gyre: A comparison between the low and record-low ice conditions of August 2011 and 2012

Abstract. The dissolved oxygen-to-argon ratio (O2∕Ar) in the oceanic mixed layer has been widely used to estimate net community production (NCP), which is the difference between gross primary production and community respiration; it is a measure of the strength of the biological pump. In order to obtain the high-resolution distribution of NCP and improve our understanding of its regulating factors in the slope region of the northern South China Sea (SCS), we conducted continuous measurements of dissolved O2, Ar, and CO2 with membrane inlet mass spectrometry (MIMS) during two cruises in October 2014 and June 2015. An overall autotrophic condition was observed in the study region in both cruises with an average Δ(O2∕Ar) of 1.1 % ± 0.9 % in October 2014 and 2.7 % ± 2.8 % in June 2015. NCP was on average 11.5 ± 8.7 mmol C m−2 d−1 in October 2014 and 11.6 ± 12.7 mmol C m−2 d−1 in June 2015. Correlations between dissolved inorganic nitrogen (DIN), Δ(O2∕Ar), and NCP were observed in both cruises, indicating that NCP is subject to the nitrogen limitation in the study region. In June 2015, we observed a rapid response of the ecosystem to the episodic nutrient supply induced by eddies. Eddy-entrained shelf water intrusion, which supplied large amounts of terrigenous nitrogen to the study region, promoted NCP in the study region by potentially more than threefold. In addition, upwelling brought large uncertainties to the estimation of NCP in the core region of the cold eddy (cyclone) in June 2015. The deep euphotic depth in the SCS and the absence of correlation between NCP and the average photosynthetically available radiation (PAR) in the mixed layer in the autumn indicate that light availability may not be a significant limitation on NCP in the SCS. This study helps us to understand the carbon cycle in the highly dynamic shelf system.

Mixed‐layer dynamics exert a first order control on nutrient and light availability for phytoplankton. In this study, we examine the influence of mixed‐layer dynamics on net community production (NCP) in the Southern Ocean on intra‐seasonal, seasonal, interannual, and decadal timescales, using biogeochemical Argo floats and satellite‐derived NCP estimates during the period from 1997 to 2020. On intraseasonal timescales, the shoaling of the mixed layer is more likely to enhance NCP in austral spring and winter, suggesting an alleviation of light limitation. As expected, NCP generally increases with light availability on seasonal timescales. On interannual timescales, NCP is correlated with mixed layer depth (MLD) and mixed‐layer‐averaged photosynthetically active radiation (PAR) in austral spring and winter, especially in regions with deeper mixed layers. Though recent studies have argued that winter MLD controls the subsequent growing season's iron and light availability, the limited number of Argo float observations contemporaneous with our satellite observations do not show a significant correlation between NCP and the previous‐winter's MLD on interannual timescales. Over the 1997–2020 period, we observe regional trends in NCP (e.g., increasing around S. America), but no trend for the entire Southern Ocean. Overall, our results show that the dependence of NCP on MLD is a complex function of timescales.

Carbon-13 constraints on the seasonal inorganic carbon budget at the BATS site in the northwestern Sargasso Sea

Depth-dependent fate of biologically-consumed dimethylsulfide in the Sargasso Sea

The sea‐air biological O2 flux assessed from measurements of surface O2 supersaturation in excess of Ar supersaturation (“O2 bioflux”) is increasingly being used to constrain net community production (NCP) in the upper ocean mixed layer. In making these calculations, one generally assumes that NCP is at steady state, mixed layer depth is constant, and there is no O2 exchange across the base of the mixed layer. The object of this paper is to evaluate the magnitude of errors introduced by violations of these assumptions. Therefore, we examine the differences between the sea‐air biological O2 flux and NCP in the Southern Ocean mixed layer as calculated using two ocean biogeochemistry general circulation models. In this approach, NCP is considered a known entity in the prognostic model, whereas O2 bioflux is estimated using the model‐predicted O2/Ar ratio to compute the mixed layer biological O2 saturation and the gas transfer velocity to calculate flux. We find that the simulated biological O2 flux gives an accurate picture of the regional‐scale patterns and trends in model NCP. However, on local scales, violations of the assumptions behind the O2/Ar method lead to significant, non‐uniform differences between model NCP and biological O2 flux. These errors arise from two main sources. First, venting of biological O2 to the atmosphere can be misaligned from NCP in both time and space. Second, vertical fluxes of oxygen across the base of the mixed layer complicate the relationship between NCP and the biological O2 flux. Our calculations show that low values of O2 bioflux correctly register that NCP is also low (<10 mmol m−2 day−1), but fractional errors are large when rates are this low. Values between 10 and 40 mmol m−2 day−1 in areas with intermediate mixed layer depths of 30 to 50 m have the smallest absolute and relative errors. Areas with O2 bioflux higher than 30 mmol m−2 day−1 and mixed layers deeper than 40 m tend to underestimate NCP by up to 20 mmol m−2 day−1. Excluding time periods when mixed layer biological O2 is undersaturated, O2 bioflux underestimates time‐averaged NCP by 5%–15%. If these time periods are included, O2 bioflux underestimates mixed layer NCP by 20%–35% in the Southern Ocean. The higher error estimate is relevant if one wants to estimate seasonal NCP since a significant amount of biological production takes place when mixed layer biological O2 is undersaturated.

The role of spring biological production for the air—sea CO2 flux was quantified in the Western Subarctic Gyre (48°N, 165°E), where the vertical profile of temperature revealed the existence of a temperature minimum (Tmin) layer in the North Pacific. The vertical profiles of temperature, salinity, dissolved oxygen, nutrients and dissolved inorganic carbon, DIC, in the upper water column were significantly variable year by year in spring, 1996–2000. Correspondingly, surface seawater at this site in spring was supersaturated with CO2 in 1997, 1999 and 2000, but was undersaturated in 1996 and 1998. The concentrations of DIC and nutrients in the winter mixed layer were estimated from those in the Tmin layer in spring with a correction for particle decomposition based on the apparent oxygen utilization. The net community production (NCP) and air—sea CO2 flux from winter to spring were calculated from the vertically integrated deficits of DIC and nutrients in the upper water column between the two seasons. The calculation of the carbon budget indicated large interannual variations of NCP (0–13 mmol m−2 d−1) and CO2 efflux (4–16 mmol m−2 d−1) for this period. The CO2 efflux was generally low in the year when NCP was high. The close coupling between biological production and CO2 efflux suggested the important role of the changes in the mixed-layer depth, as a key process controlling both processes, especially of the timing, so that a decrease in the mixed-layer depth could result in the activation of biological production. The early biological consumption of the surface DIC concentration could shorten the period for acting as a source for atmospheric CO2 and depress the CO2 efflux in the Western Subarctic Gyre from winter to spring in 1996 and 1998. On the contrary, in 1997, persistently deep vertical mixing until late spring could suppress the biological activity and give rise to long-lasting CO2 efflux.

Measurements of net community production (NCP) provide an upper constraint on the strength of the oceanic biological pump, the dominant mechanism for removing CO2 from the ocean surface and sequestering it at depth. In this investigation, our objectives were to describe the spatial and temporal variability of NCP associated with the spring ice-edge bloom in Baffin Bay and to identify the key environmental drivers controlling its variability. Using data collected between June 9 and July 10, 2016, we estimated NCP based on (1) underway measurements of surface water oxygen to argon ratios (O2:Ar), (2) underway measurements of the partial pressure of CO2, and (3) seasonal nitrate drawdown from discrete samples. These multiple approaches displayed high NCP (up to 5.7 mol C m–2) in eastern Baffin Bay, associated with modified Atlantic waters, and low NCP (<1 mol C m–2) in the presence of Arctic outflow waters in western Baffin Bay. Arctic outflow waters were characterized by low surface salinities and nitrate concentrations, suggesting that high freshwater content may have limited the nutrient availability of these waters. Different integration depths and timescales associated with each NCP approach were exploited to understand the temporal progression and succession of the bloom, revealing that the bloom was initiated under ice up to 15 days prior to ice retreat and that a large portion of NCP in eastern Baffin Bay (potentially up to 70%) was driven by primary production occurring below the surface-mixed layer.

Changes in biological production in the mixed water region (MWR) of the northwestern North Pacific during the last 27 kyr