ALBERT

Description: Ocean alkalinity enhancement (OAE) is a negative emissions technology (NET) that shows significant potential for climate change mitigation. By increasing the bicarbonate ion concentration in ocean water, OAE could enhance long-term carbon storage and mitigate ocean acidification. However, the side effects and/or potential co-benefits of OAE on natural planktonic communities remain poorly understood. To address this knowledge gap, a mesocosm experiment was conducted in the oligotrophic waters of Gran Canaria. A CO2-equilibrated total alkalinity (TA) gradient was employed in increments of 300 µmol L−1, ranging from ∼ 2400 to ∼ 4800 µmol L−1. This study represents the first attempt to evaluate the potential impacts of OAE on planktonic communities under natural conditions. The results show that net community production (NCP), gross production (GP), community respiration (CR) rates, and the metabolic balance (GP:CR) did not exhibit a linear response to the whole alkalinity gradient. Instead, significant polynomial and linear regression models were observed for all rates up to ΔTA 1800 µmol L−1, in relation to the dissolved inorganic carbon (DIC) concentrations. Notably, the ΔTA 1500 and 1800 µmol L−1 treatments showed peaks in NCP shifting from a heterotrophic to an autotrophic state, with NCP values of 4 and 8 µmol O2 kg−1 d−1, respectively. These peaks and the optimum curve were also reflected in the nanoplankton abundance, size-fractionated chlorophyll a, and 14C uptake data. Furthermore, abiotic precipitation occurred in the highest treatment after day 21, but no impact on the measured parameters was detected. Overall, a damaging effect of CO2-equilibrated OAE in the range applied here on phytoplankton primary production, community metabolism, and composition could not be inferred. In fact, a potential co-benefit to OAE was observed in the form of the positive curvilinear response to the DIC gradient up to the ΔTA 1800 treatment. Further experimental research at this scale is key to gain a better understanding of the short- and long-term effects of OAE on planktonic communities.

Description: Nitrous oxide (N2O) is a long-lived potent greenhouse gas and stratospheric ozone-depleting substance, which has been accumulating in the atmosphere since the pre-industrial period. The mole fraction of atmospheric N2O has increased by nearly 25 % from 270 parts per billion (ppb) in 1750 to 336 ppb in 2022, with the fastest annual growth rate since 1980 of more than 1.3 ppb yr-1 in both 2020 and 2021. As a core component of our global greenhouse gas assessments coordinated by the Global Carbon Project (GCP), we present a global N2O budget that incorporates both natural and anthropogenic sources and sinks, and accounts for the interactions between nitrogen additions and the biochemical processes that control N2O emissions. We use Bottom-Up (BU: inventory, statistical extrapolation of flux measurements, process-based land and ocean modelling) and Top-Down (TD: atmospheric measurement-based inversion) approaches. We provide a comprehensive quantification of global N2O sources and sinks in 21 natural and anthropogenic categories in 18 regions between 1980 and 2020. We estimate that total annual anthropogenic N2O emissions increased 40 % (or 1.9 Tg N yr-1) in the past four decades (1980–2020). Direct agricultural emissions in 2020, 3.9 Tg N yr−1 (best estimate) represent the large majority of anthropogenic emissions, followed by other direct anthropogenic sources (including ‘Fossil fuel and industry’, ‘Waste and wastewater’, and ‘Biomass burning’ (2.1 Tg N yr−1), and indirect anthropogenic sources (1.3 Tg N yr−1). For the year 2020, our best estimate of total BU emissions for natural and anthropogenic sources was 18.3 (lower-upper bounds: 10.5–27.0) Tg N yr-1, close to our TD estimate of 17.0 (16.6–17.4) Tg N yr-1. For the period 2010–2019, the annual BU decadal-average emissions for natural plus anthropogenic sources were 18.1 (10.4–25.9) Tg N yr-1 and TD emissions were 17.4 (15.8–19.20 Tg N yr-1. The once top emitter Europe has reduced its emissions since the 1980s by 31 % while those of emerging economies have grown, making China the top emitter since the 2010s. The observed atmospheric N2O concentrations in recent years have exceeded projected levels under all scenarios in the Coupled Model Intercomparison Project Phase 6 (CMIP6), underscoring the urgency to reduce anthropogenic N2O emissions. To evaluate mitigation efforts and contribute to the Global Stocktake of the United Nations Framework Convention on Climate Change, we propose establishing a global network for monitoring and modeling N2O from the surface through the stratosphere. The data presented in this work can be downloaded from https://doi.org/10.18160/RQ8P-2Z4R (Tian et al. 2023).

Description: In marine ecosystems, most physiological, ecological, or physical processes are size dependent. These include metabolic rates, the uptake of carbon and other nutrients, swimming and sinking velocities, and trophic interactions, which eventually determine the stocks of commercial species, as well as biogeochemical cycles and carbon sequestration. As such, broad-scale observations of plankton size distribution are important indicators of the general functioning and state of pelagic ecosystems under anthropogenic pressures. Here, we present the first global datasets of the Pelagic Size Structure database (PSSdb), generated from plankton imaging devices. This release includes the bulk particle normalized biovolume size spectrum (NBSS) and the bulk particle size distribution (PSD), along with their related parameters (slope, intercept, and R2) measured within the epipelagic layer (0–200 m) by three imaging sensors: the Imaging FlowCytobot (IFCB), the Underwater Vision Profiler (UVP), and benchtop scanners. Collectively, these instruments effectively image organisms and detrital material in the 7–10 000 µm size range. A total of 92 472 IFCB samples, 3068 UVP profiles, and 2411 scans passed our quality control and were standardized to produce consistent instrument-specific size spectra averaged to 1° × 1° latitude and longitude and by year and month. Our instrument-specific datasets span most major ocean basins, except for the IFCB datasets we have ingested, which were exclusively collected in northern latitudes, and cover decadal time periods (2013–2022 for IFCB, 2008–2021 for UVP, and 1996–2022 for scanners), allowing for a further assessment of the pelagic size spectrum in space and time. The datasets that constitute PSSdb's first release are available at https://doi.org/10.5281/zenodo.11050013 (Dugenne et al., 2024b). In addition, future updates to these data products can be accessed at https://doi.org/10.5281/zenodo.7998799.

Description: Upwelling of nutrient-rich deep waters make eastern boundary upwelling systems (EBUSs), such as the Humboldt Current system, hot spots of marine productivity. Associated settling of organic matter to depth and consecutive aerobic decomposition results in large subsurface water volumes being oxygen depleted. Under these circumstances, organic matter remineralisation can continue via denitrification, which represents a major loss pathway for bioavailable nitrogen. Additionally, anaerobic ammonium oxidation can remove significant amounts of nitrogen in these areas. Here we assess the interplay of suboxic water upwelling and nitrogen cycling in a manipulative offshore mesocosm experiment. Measured denitrification rates in incubations with water from the oxygen-depleted bottom layer of the mesocosms (via 15N label incubations) mostly ranged between 5.5 and 20 nmol N2 L−1 h−1 (interquartile range), reaching up to 80 nmol N2 L−1 h−1. However, actual in situ rates in the mesocosms, estimated via Michaelis–Menten kinetic scaling, did most likely not exceed 0.2–4.2 nmol N2 L−1 h−1 (interquartile range) due to substrate limitation. In the surrounding Pacific, measured denitrification rates were similar, although indications of substrate limitation were detected only once. In contrast, anammox (anaerobic ammonium oxidation) made only a minor contribution to the overall nitrogen loss when encountered in both the mesocosms and the Pacific Ocean. This was potentially related to organic matter C / N stoichiometry and/or process-specific oxygen and hydrogen sulfide sensitivities. Over the first 38 d of the experiment, total nitrogen loss calculated from in situ rates of denitrification and anammox was comparable to estimates from a full nitrogen budget in the mesocosms and ranged between ∼ 1 and 5.5 µmol N L−1. This represents up to ∼  20 % of the initially bioavailable inorganic and organic nitrogen standing stocks. Interestingly, this loss is comparable to the total amount of particulate organic nitrogen that was exported into the sediment traps at the bottom of the mesocosms at about 20 m depth. Altogether, this suggests that a significant portion, if not the majority of nitrogen that could be exported to depth, is already lost, i.e. converted to N2 in a relatively shallow layer of the surface ocean, provided that there are oxygen-deficient conditions like those during coastal upwelling in our study. Published data for primary productivity and nitrogen loss in all EBUSs reinforce such conclusion.

Description: This study contributes to the inaugural exploration of non-equilibrated Ocean Alkalinity Enhancement (OAE). The manipulation of Total Alkalinity (TA), involving silicate and calcium-based ∆TA gradients ranging from 0 to 600 µmol · L-1, was conducted without prior CO2 sequestration, under natural conditions and at a mesocosm scale (~60 m3). The resulting impact included an increase in pH and a decrease in pCO2, sustained across the experiment, as full natural equilibration via sea-gas exchange did not occur. Implemented in a neritic system under post-bloom conditions, a midway mixing event was simulated. Following an inorganic nutrient addition, discernible delays in bloom formation, as indicated by the Gross Production (GP) and Net Community Production (NCP) rates, as well as by the chlorophyll-a (Chla) concentrations, in relation to the ∆TA gradient, were observed. Notably, the delay was more pronounced for the calcium treatment set compared to the silicate one, where low TA treatments exhibited earlier responses than high TA ones. This delay is likely attributed to the previously documented, species-specific negative relationships between high pH/lowCO2 levels and phytoplankton growth rates. This study underscores the need for further investigation into the implications of this response pattern in terms of trophic transfer and seasonal suitability. Further, it is anticipated that a wider delay in bloom formation would be evident with a larger non equilibrated TA gradient. Thus, highlighting the importance of exploring variations in TA limits for a comprehensive understanding of the OAE’s impacts.

Description: Observed oxygen minimum zones (OMZs) in the tropical Pacific Ocean are located above intermediate-depth waters (IDWs), defined here as the 500–1500 m water layer. Typical climate models do not represent IDW properties well and are characterized by OMZs that are too deep-reaching. We analyze the role of the IDW in the misrepresentation of oxygen levels in a heterogeneous subset of ocean models characterized by a horizontal resolution ranging from 0.1 to 2.8∘. First, we show that forcing the extratropical boundaries (30∘ S and N) to observed oxygen values results in a significant increase in oxygen levels in the intermediate eastern tropical region. Second, we highlight the fact that the Equatorial Intermediate Current System (EICS) is a key feature connecting the western and eastern part of the basin. Typical climate models lack in representing crucial aspects of this supply at intermediate depth, as the EICS is basically absent in models characterized by a resolution lower than 0.25∘. These two aspects add up to a “cascade of biases” that hampers the correct representation of oxygen levels at intermediate depth in the eastern tropical Pacific Ocean and potentially future OMZ projections.

Description: In this study we look beyond the previously studied effects of oceanic CO2 injections on atmospheric and oceanic reservoirs and also account for carbon cycle and climate feedbacks between the atmosphere and the terrestrial biosphere. Considering these additional feedbacks is important since backfluxes from the terrestrial biosphere to the atmosphere in response to reducing atmospheric CO2 can further offset the targeted reduction. To quantify these dynamics we use an Earth system model of intermediate complexity to simulate direct injection of CO2 into the deep ocean as a means of emissions mitigation during a high CO2 emission scenario. In three sets of experiments with different injection depths, we simulate a 100-year injection period of a total of 70 GtC and follow global carbon cycle dynamics over another 900 years. In additional parameter perturbation runs, we varied the default terrestrial photosynthesis CO2 fertilization parameterization by ±50 % in order to test the sensitivity of this uncertain carbon cycle feedback to the targeted atmospheric carbon reduction through direct CO2 injections. Simulated seawater chemistry changes and marine carbon storage effectiveness are similar to previous studies. As expected, by the end of the injection period avoided emissions fall short of the targeted 70 GtC by 16–30 % as a result of carbon cycle feedbacks and backfluxes in both land and ocean reservoirs. The target emissions reduction in the parameter perturbation simulations is about 0.2 and 2 % more at the end of the injection period and about 9 % less to 1 % more at the end of the simulations when compared to the unperturbed injection runs. An unexpected feature is the effect of the model's internal variability of deep-water formation in the Southern Ocean, which, in some model runs, causes additional oceanic carbon uptake after injection termination relative to a control run without injection and therefore with slightly different atmospheric CO2 and climate. These results of a model that has very low internal climate variability illustrate that the attribution of carbon fluxes and accounting for injected CO2 may be very challenging in the real climate system with its much larger internal variability.

Description: The addition of carbonate minerals to seawater through an artificial ocean alkalinity enhancement (OAE) process increases the concentrations of hydroxide, bicarbonate, and carbonate ions. This leads to changes in the pH and the buffering capacity of the seawater. Consequently, OAE could have relevant effects on marine organisms and in the speciation and concentration of trace metals that are essential for their physiology. During September and October 2021, a mesocosm experiment was carried out in the coastal waters of Gran Canaria (Spain), consisting on the controlled variation of total alkalinity (TA). Different concentrations of carbonate salts (NaHCO3 and Na2CO3) previously homogenized were added to each mesocosm to achieve an alkalinity gradient between Δ0 to Δ2400 µmol L−1. The lowest point of the gradient was 2400 µmol kg−1, being the natural alkalinity of the medium, and the highest point was 4800 µmol kg−1. Iron (Fe) speciation was monitored during this experiment to analyse total dissolved iron (TdFe, unfiltered samples), dissolved iron (dFe, filtered through a 0.2 µm pore size filter), soluble iron (sFe, filtered through a 0.02 µm pore size filter), dissolved labile iron (dFe′), iron-binding ligands (LFe), and their conditional stability constants () because of change due to OAE and the experimental conditions in each mesocosm. Observed iron concentrations were within the expected range for coastal waters, with no significant increases due to OAE. However, there were variations in Fe size fractionation during the experiment. This could potentially be due to chemical changes caused by OAE, but such an effect is masked by the stronger biological interactions. In terms of size fractionation, sFe was below 1.0 nmol L−1, dFe concentrations were within 0.5–4.0 nmol L−1, and TdFe was within 1.5–7.5 nmol L−1. Our results show that over 99 % of Fe was complexed, mainly by L1 and L2 ligands with ranging between 10.92 ± 0.11 and 12.68 ± 0.32, with LFe ranging from 1.51 ± 0.18 to 12.3 ± 1.8 nmol L−1. Our data on iron size fractionation, concentration, and iron-binding ligands substantiate that the introduction of sodium salts in this mesocosm experiment did not modify iron dynamics. As a consequence, phytoplankton remained unaffected by alterations in this crucial element.

Description: The Weddell Sea Polynya (WSP) is a large opening within the sea ice cover of the Weddell Sea sector. It has been a rare event in the satellite period, appearing between 1973 and 1976 and again in 2016/2017. Coupled modelling studies have suggested that there may be a large-scale atmospheric response to the WSP. Here, the direct atmospheric response to the WSP is estimated from atmosphere-only numerical experiments. Three different models, the HadGEM3 UK Met Office model, the ECHAM5 Max Planck Institute model, and the OpenIFS ECMWF model, each at two different resolutions, are used to test the robustness of our results. The use of large ensembles reduces the weather variability and isolates the atmospheric response. Results show a large (∼100-200 Wm-2) turbulent air-sea flux anomaly above the polynya. The response to the WSP is local and of short duration (barely outlasting the WSP) with a similar magnitude and spatial pattern of lower-tropospheric warming and increase in precipitation in all six configurations. All models show a weak decrease in surface pressure over the WSP, but this response is small (∼2 hPa) in comparison to internal variability. The dynamic response is inconsistent between models and resolutions above the boundary layer, suggesting a weak or null response that is covered by internal variability aloft. The higher resolution does not alter the pattern of the response but increases its magnitude by up ∼50% in two of the three models. The response is influenced by natural variability in the westerly jet. The models perform well against ERA5 reanalysis data for the 1974 WSP in spatial response and magnitude, showing a turbulent heat flux of approximately 150 W m-2.