Global‐scale evaluation of coastal ocean alkalinity enhancement in a fully coupled Earth system model

The Paris Agreement plans for “net-zero” carbon dioxide (CO2) emissions during the second half of the 21st century. However, reducing emissions from some sectors is challenging, and “net-zero” permits carbon dioxide removal (CDR) activities. One CDR scheme is ocean alkalinity enhancement (OAE), which proposes dissolving basic minerals into seawater to increase its buffering capacity for CO2. While modeling studies have often investigated OAE at basin or global scale, some proposals focus on readily accessible coastal shelves, with TA added through the dissolution of seafloor olivine sands. Critically, by settling and dissolving sands on shallow seafloors, this retains the added TA in near-surface waters in direct contact with atmospheric CO2. To investigate this, we add dissolved TA at a rate of ∼29 Teq y−1 to the global shelves (<100m) of an Earth system model (UKESM1) running a high emissions scenario. As UKESM1 is fully coupled, wider effects of OAE-mediated increase in ocean CO2 uptake –e.g. atmospheric xCO2, air temperature and marine pH– are fully quantified. Applying OAE from 2020 to 2100 decreases atmospheric xCO2 ∼10 ppm, and increases air-to-sea CO2 uptake ∼8%. In-line with other studies, CO2 uptake per unit of TA added occurs at a rate of ∼0.8 mol C (mol TA)−1. Significantly for monitoring, advection of added TA results in ∼50% of CO2 uptake occurring remotely from OAE operations, and the model also exhibits noticeable land carbon reservoir changes. While practical uncertainties and model representation caveats remain, this analysis estimates the effectiveness of this specific OAE scheme to assist with net-zero planning

BGC-val: a model- and grid-independent Python toolkit to evaluate marine biogeochemical models

The biogeochemical evaluation toolkit, BGC-val, is a model- and grid-independent Python toolkit that has been built to evaluate marine biogeochemical models using a simple interface. Here, we present the ideas that motivated the development of the BGC-val software framework, introduce the code structure, and show some applications of the toolkit using model results from the Fifth Climate Model Intercomparison Project (CMIP5). A brief outline of how to access and install the repository is presented in Appendix A, but the specific details on how to use the toolkit are kept in the code repository. The key ideas that directed the toolkit design were model and grid independence, front-loading analysis functions and regional masking, interruptibility, and ease of use. We present each of these goals, why they were important, and what we did to address them. We also present an outline of the code structure of the toolkit illustrated with example plots produced by the toolkit. After describing BGC-val, we use the toolkit to investigate the performance of the marine physical and biogeochemical quantities of the CMIP5 models and highlight some predictions about the future state of the marine ecosystem under a business-as-usual CO2 concentration scenario (RCP8.5)

Modelling of the anthropogenic tritium transient and its decay product helium-3 in the Mediterranean Sea using a high-resolution regional model

International audienceThis numerical study provides the first simulation of the anthropogenic tritium invasion and its decay product helium-3 (3 He) in the Mediterranean Sea. The simulation covers the entire tritium (3 H) transient generated by the atmospheric nuclear weapons tests performed in the 1950s and early 1960s and is run till 2011. Tritium, helium-3 and their derived age estimates are particularly suitable for studying intermediate and deep-water ventilation and spreading of water masses at intermediate/deep levels. The simulation is made using a high-resolution regional model NEMO (Nucleus for European Modelling of the Ocean), in a regional configuration for the Mediterranean Sea called MED12, forced at the surface with prescribed tritium evolution derived from observations. The simulation is compared to measurements of tritium and helium-3 performed along large-scale transects in the Mediterranean Sea during the last few decades on cruises of R/V Meteor: M5/6, M31/1, M44/4, M51/2, M84/3, and R/V Poseidon: 234. The results show that the input function used for the tritium generates a realistic distribution of the main hydrographic features of the Mediterranean Sea circulation. In the eastern basin, the results highlight the weak formation of Adriatic Deep Water in the model, which explains its weak contribution to the Eastern Mediterranean Deep Water (EMDW) in the Ionian sub-basin. It produces a realistic representation of the Eastern Mediterranean Transient (EMT) signal, simulating a deep-water formation in the Aegean sub-basin at the beginning of 1993, with a realistic timing of deep-water renewal in the eastern basin

Earth system music: music generated from the United Kingdom Earth System Model (UKESM1)

Scientific data are almost always represented graphically in figures or in videos. With the ever-growing interest from the general public in understanding climate sciences, it is becoming increasingly important that scientists present this information in ways that are both accessible and engaging to non-experts. In this pilot study, we use time series data from the first United Kingdom Earth System Model (UKESM1) to create six procedurally generated musical pieces. Each of these pieces presents a unique aspect of the ocean component of the UKESM1, either in terms of a scientific principle or a practical aspect of modelling. In addition, each piece is arranged using a different musical progression, style and tempo. These pieces were created in the Musical Instrument Digital Interface (MIDI) format and then performed by a digital piano synthesiser. An associated video showing the time development of the data in time with the music was also created. The music and video were published on the lead author's YouTube channel. A brief description of the methodology was also posted alongside the video. We also discuss the limitations of this pilot study and describe several approaches to extend and expand upon this work

Evaluating the physical and biogeochemical state of the global ocean component of UKESM1 in CMIP6 historical simulations

The ocean plays a key role in modulating the climate of the Earth system (ES). At the present time it is also a major sink both for the carbon dioxide (CO2) released by human activities and for the excess heat driven by the resulting atmospheric greenhouse effect. Understanding the ocean's role in these processes is critical for model projections of future change and its potential impacts on human societies. A necessary first step in assessing the credibility of such future projections is an evaluation of their performance against the present state of the ocean. Here we use a range of observational fields to validate the physical and biogeochemical performance of the ocean component of UKESM1, a new Earth system model (ESM) for CMIP6 built upon the HadGEM3-GC3.1 physical climate model. Analysis focuses on the realism of the ocean's physical state and circulation, its key elemental cycles, and its marine productivity. UKESM1 generally performs well across a broad spectrum of properties, but it exhibits a number of notable biases. Physically, these include a global warm bias inherited from model spin-up, excess northern sea ice but insufficient southern sea ice and sluggish interior circulation. Biogeochemical biases found include shallow remineralization of sinking organic matter, excessive iron stress in regions such as the equatorial Pacific, and generally lower surface alkalinity that results in decreased surface and interior dissolved inorganic carbon (DIC) concentrations. The mechanisms driving these biases are explored to identify consequences for the behaviour of UKESM1 under future climate change scenarios and avenues for model improvement. Finally, across key biogeochemical properties, UKESM1 improves in performance relative to its CMIP5 precursor and performs well alongside its fellow members of the CMIP6 ensemble

Reproducible and relocatable regional ocean modelling: Fundamentals and practices

In response to an increasing demand for bespoke or tailored regional ocean modelling configurations, we outline fundamental principles and practices that can expedite the process to generate new configurations. The paper develops the principle of Reproducibility and advocates adherence by presenting benefits to the community and user. The elements to this principle are reproducible workflows and standardised assessment, with additional effort over existing working practices being balanced against the added value generated. The paper then decomposes the complex build process, for a new regional ocean configuration, into stages and presents guidance, advice and insight on each component. This advice is compiled from across the user community, is presented in the context of NEMOv4, though aims to transcend NEMO version. Detail and region specific worked examples are linked in companion repositories and DOIs. The aim is to broaden the user community skill base, and to accelerate development of new configurations in order to increase available time exploiting the configurations

Global Surface Ocean Acidification Indicators From 1750 to 2100

Accurately predicting future ocean acidification (OA) conditions is crucial for advancing OA research at regional and global scales, and guiding society's mitigation and adaptation efforts. This study presents a new model-data fusion product covering 10 global surface OA indicators based on 14 Earth System Models (ESMs) from the Coupled Model Intercomparison Project Phase 6 (CMIP6), along with three recent observational ocean carbon data products. The indicators include fugacity of carbon dioxide, pH on total scale, total hydrogen ion content, free hydrogen ion content, carbonate ion content, aragonite saturation state, calcite saturation state, Revelle Factor, total dissolved inorganic carbon content, and total alkalinity content. The evolution of these OA indicators is presented on a global surface ocean 1° × 1° grid as decadal averages every 10 years from preindustrial conditions (1750), through historical conditions (1850–2010), and to five future Shared Socioeconomic Pathways (2020–2100): SSP1-1.9, SSP1-2.6, SSP2-4.5, SSP3-7.0, and SSP5-8.5. These OA trajectories represent an improvement over previous OA data products with respect to data quantity, spatial and temporal coverage, diversity of the underlying data and model simulations, and the provided SSPs. The generated data product offers a state-of-the-art research and management tool for the 21st century under the combined stressors of global climate change and ocean acidification. The gridded data product is available in NetCDF at the National Oceanic and Atmospheric Administration (NOAA) National Centers for Environmental Information: https://www.ncei.noaa.gov/data/oceans/ncei/ocads/metadata/0259391.html, and global maps of these indicators are available in jpeg at: https://www.ncei.noaa.gov/access/ocean-carbon-acidification-data-system/synthesis/surface-oa-indicators.html

Description and evaluation of aerosol in UKESM1 and HadGEM3-GC3.1 CMIP6 historical simulations

We document and evaluate the aerosol schemes as implemented in the physical and Earth system models, HadGEM3-GC3.1 (GC3.1) and UKESM1, which are contributing to the 6th Coupled Model Intercomparison Project (CMIP6). The simulation of aerosols in the present-day period of the historical ensemble of these models is evaluated against a range of observations. Updates to the aerosol microphysics scheme are documented as well as differences in the aerosol representation between the physical and Earth system configurations. The additional Earth-system interactions included in UKESM1 leads to differences in the emissions of natural aerosol sources such as dimethyl sulfide, mineral dust and organic aerosol and subsequent evolution of these species in the model. UKESM1 also includes a stratospheric-tropospheric chemistry scheme which is fully coupled to the aerosol scheme, while GC3.1 employs a simplified aerosol chemistry mechanism driven by prescribed monthly climatologies of the relevant oxidants. Overall, the simulated speciated aerosol mass concentrations compare reasonably well with observations. Both models capture the negative trend in sulfate aerosol concentrations over Europe and the eastern United States of America (US) although the models tend to underestimate the sulfate concentrations in both regions. Interactive emissions of biogenic volatile organic compounds in UKESM1 lead to an improved agreement of organic aerosol over the US. Simulated dust burdens are similar in both models despite a 2-fold difference in dust emissions. Aerosol optical depth is biased low in dust source and outflow regions but performs well in other regions compared to a number of satellite and ground-based retrievals of aerosol optical depth. Simulated aerosol number concentrations are generally within a factor of 2 of the observations with both models tending to overestimate number concentrations over remote ocean regions, apart from at high latitudes, and underestimate over Northern Hemisphere continents. Finally, a new primary marine organic aerosol source is implemented in UKESM1 for the first time. The impact of this new aerosol source is evaluated. Over the pristine Southern Ocean, it is found to improve the seasonal cycle of organic aerosol mass and cloud droplet number concentrations relative to GC3.1 although underestimations in cloud droplet number concentrations remain. This paper provides a useful characterization of the aerosol climatology in both models facilitating the understanding of the numerous aerosol-climate interaction studies that will be conducted as part of CMIP6 and beyond

DMTs and Covid-19 severity in MS: a pooled analysis from Italy and France

We evaluated the effect of DMTs on Covid-19 severity in patients with MS, with a pooled-analysis of two large cohorts from Italy and France. The association of baseline characteristics and DMTs with Covid-19 severity was assessed by multivariate ordinal-logistic models and pooled by a fixed-effect meta-analysis. 1066 patients with MS from Italy and 721 from France were included. In the multivariate model, anti-CD20 therapies were significantly associated (OR&nbsp;=&nbsp;2.05, 95%CI&nbsp;=&nbsp;1.39–3.02, p&nbsp;&lt;&nbsp;0.001) with Covid-19 severity, whereas interferon indicated a decreased risk (OR&nbsp;=&nbsp;0.42, 95%CI&nbsp;=&nbsp;0.18–0.99, p&nbsp;=&nbsp;0.047). This pooled-analysis confirms an increased risk of severe Covid-19 in patients on anti-CD20 therapies and supports the protective role of interferon

SARS-CoV-2 susceptibility and COVID-19 disease severity are associated with genetic variants affecting gene expression in a variety of tissues

Variability in SARS-CoV-2 susceptibility and COVID-19 disease severity between individuals is partly due to genetic factors. Here, we identify 4 genomic loci with suggestive associations for SARS-CoV-2 susceptibility and 19 for COVID-19 disease severity. Four of these 23 loci likely have an ethnicity-specific component. Genome-wide association study (GWAS) signals in 11 loci colocalize with expression quantitative trait loci (eQTLs) associated with the expression of 20 genes in 62 tissues/cell types (range: 1:43 tissues/gene), including lung, brain, heart, muscle, and skin as well as the digestive system and immune system. We perform genetic fine mapping to compute 99% credible SNP sets, which identify 10 GWAS loci that have eight or fewer SNPs in the credible set, including three loci with one single likely causal SNP. Our study suggests that the diverse symptoms and disease severity of COVID-19 observed between individuals is associated with variants across the genome, affecting gene expression levels in a wide variety of tissue types