ALBERT

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: FS Maria S Merian – MSM129/2 07.06.2024 – 06.07.2024 St. John’s (Kanada) – Reykjavik (Island) 3. Wochenbericht (17.06. – 23.06.2024)

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.

Description: In September 2021, volcanic aerosol (mainly freshly formed sulfate plumes) originating from the eruption of Cumbre Vieja on La Palma, Canary Islands, Spain, crossed Cabo Verde at altitudes below 2 km. On 24 September 2021, an extraordinary large aerosol optical depth (AOD) close to 1 (daily mean at 500 nm) was observed at Mindelo, Cabo Verde. This event provided favorable conditions to obtain lidar-derived profiles of extinction and backscatter coefficients, lidar ratio, and depolarization ratio at 355, 532 and 1064 nm in the sulfate aerosol plume. A novel feature of the lidar system operated at Mindelo is the availability of extinction, lidar ratio and depolarization measurements at 1064 nm in addition to the standard wavelengths of 355 and 532 nm. Having measurements of these parameters at all three wavelengths is a major advantage for the aerosol characterization and in aerosol typing efforts as the lidar ratio and the particle linear depolarization ratio are key parameters for this purpose. In this article, we present the key results of the lidar observations obtained on one specific day, namely on 24 September 2021 at 04:38-05:57 UTC, including the first ever measurements of the particle extinction coefficient, the lidar ratio and the depolarization ratio at 1064 nm for volcanic sulfate, and discuss the findings in terms of aerosol optical properties and mass concentrations by comparison with a reference observation (16 September 2021) representing the typical background conditions before the start of the eruptions. We found an unusual high particle extinction coefficient of 721 +/- 51, 549 +/- 38 and 178 +/- 13 Mm - 1 , as well as an enhanced lidar ratio of 66.9 +/- 10.1, 60.2 +/- 9.2 and 30.8 +/- 8.7 sr at 355, 532 and 1064 nm, respectively, in the sulfate-dominated planetary boundary layer (PBL). The particle linear depolarization ratio was 〈= 0.9 % at all respective wavelengths. It is the first time that lidar-derived intensive aerosol optical properties could be derived for volcanic sulfate at all three wavelengths, and thus it is a highly valuable data set for global aerosol characterization. The lidar analysis also revealed a sulfate-related AOD of about 0.35 +/- 0.03 at 532 nm of the total PBL-related AOD of 0.43. The rest of the AOD contribution was caused by a lofted Saharan dust layer extending from 1.4 to 5 km and leading to a total AOD of 0.79 at 532 nm. Volcanic ash contribution to the observed aerosol plumes could be mostly excluded based on trajectory analysis and the observed optical properties. Peak mass concentration was 178.5 +/- 44.6 mu g m - 3 in the volcanic-influenced and sulfate-dominated polluted PBL, showing the hazardous potential of such sulfate plumes to significantly worsen local air quality even at remote locations.

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: Fjord-like estuaries are hotspots of biogeochemical cycling due to steep physicochemical gradients. The spatiotemporal distribution of nitrous oxide (N2O) within many of these systems is poorly described, especially in the southern hemisphere. The goal of this study is to describe the spatiotemporal distribution of N2O within a southern hemisphere fjord-like estuary, describe the main environmental drivers of this distribution, the air/sea flux of N2O, and the main drivers of N2O production. Cruises were undertaken in Macquarie Harbour, Tasmania to capture N2O concentrations and water column physicochemical profiles in winter (July 2022), spring (October 2022), summer (February 2023), and autumn (April 2023). N2O samples were collected at one depth at system end members, and at 5 depths at 4 stations within the harbour. Results indicate that N2O is consistently supersaturated (reaching 170 % saturation) below the system’s freshwater lens where oxygen concentrations are often hypoxic, but infrequently anoxic. In the surface lens, levels of N2O saturation vary with estimated river flow and with proximity to the system’s main freshwater endmember. The linear relationship between AOU and ΔN2O saturation indicates that nitrification is the process generating N2O in the system. When river flow was high (July and October 2022), surface water N2O was undersaturated (as low as 70 %) throughout most of the harbour. When river flow was low (February and April 2023) N2O was observed to be supersaturated at most stations. Calculated air/sea fluxes of N2O indicated that the system is generally a source of N2O to the atmosphere under weak river flow conditions and a sink during strong river flow conditions. The diapycnal flux was a minor contributor to surface water N2O concentrations, and subhalocline N2O is intercepted by the riverine surface lens and transported out of the system to the ocean during strong river flow conditions. In a changing climate, Western Tasmania is expected to receive higher winter rainfall and lower summer rainfall which may augment the source and sink dynamics of this system by enhancing the summer / autumn efflux of N2O to the atmosphere. This study is the first to report observations of N2O distribution, generation processes, and estimated diapycnal / surface N2O fluxes from this system.

Description: Numerical models are a powerful tool for investigating the dynamic processes in the interior of the Earth and other planets, but the reliability and predictive power of these discretized models depends on the numerical method as well as an accurate representation of material properties in space and time. In the specific context of geodynamic models, particle methods have been applied extensively because of their suitability for advection-dominated processes and have been used in applications such as tracking the composition of solid rock and melt in the Earth's mantle, fluids in lithospheric- and crustal-scale models, light elements in the liquid core, and deformation properties like accumulated finite strain or mineral grain size, along with many applications outside the Earth sciences.There have been significant benchmarking efforts to measure the accuracy and convergence behavior of particle methods, but these efforts have largely been limited to instantaneous solutions, or time-dependent models without analytical solutions. As a consequence, there is little understanding about the interplay of particle advection errors and errors introduced in the solution of the underlying transient, nonlinear flow equations. To address these limitations, we present two new dynamic benchmarks for transient Stokes flow with analytical solutions that allow us to quantify the accuracy of various advection methods in nonlinear flow. We use these benchmarks to measure the accuracy of our particle algorithm as implemented in the ASPECT geodynamic modeling software against commonly employed field methods and analytical solutions. In particular, we quantify if an algorithm that is higher-order accurate in time will allow for better overall model accuracy and verify that our algorithm reaches its intended optimal convergence rate. We then document that the observed increased accuracy of higher-order algorithms matters for geodynamic applications with an example of modeling small-scale convection underneath an oceanic plate and show that the predicted place and time of onset of small-scale convection depends significantly on the chosen particle advection method.Descriptions and implementations of our benchmarks are openly available and can be used to verify other advection algorithms. The availability of accurate, scalable, and efficient particle methods as part of the widely used open-source code ASPECT will allow geodynamicists to investigate complex time-dependent geodynamic processes such as elastic deformation, anisotropic fabric development, melt generation and migration, and grain damage.