The Ocean Model Intercomparison Project (OMIP) focuses on the physics and biogeochemistry of the ocean component of Earth system models participating in the sixth phase of the Coupled Model Intercomparison Project (CMIP6). OMIP aims to provide standard protocols and diagnostics for ocean models, while offering a forum to promote their common assessment and improvement. It also offers to compare solutions of the same ocean models when forced with reanalysis data (OMIP simulations) vs. when integrated within fully coupled Earth system models (CMIP6). Here we detail simulation protocols and diagnostics for OMIP's biogeochemical and inert chemical tracers. These passive-tracer simulations will be coupled to ocean circulation models, initialized with observational data or output from a model spin-up, and forced by repeating the 1948–2009 surface fluxes of heat, fresh water, and momentum. These so-called OMIP-BGC simulations include three inert chemical tracers (CFC-11, CFC-12, SF₆) and biogeochemical tracers (e.g., dissolved inorganic carbon, carbon isotopes, alkalinity, nutrients, and oxygen). Modelers will use their preferred prognostic BGC model but should follow common guidelines for gas exchange and carbonate chemistry. Simulations include both natural and total carbon tracers. The required forced simulation (omip1) will be initialized with gridded observational climatologies. An optional forced simulation (omip1-spunup) will be initialized instead with BGC fields from a long model spin-up, preferably for 2000 years or more, and forced by repeating the same 62-year meteorological forcing. That optional run will also include abiotic tracers of total dissolved inorganic carbon and radiocarbon, C_T^abio and ¹⁴C_T^abio, to assess deep-ocean ventilation and distinguish the role of physics vs. biology. These simulations will be forced by observed atmospheric histories of the three inert gases and CO₂ as well as carbon isotope ratios of CO₂. OMIP-BGC simulation protocols are founded on those from previous phases of the Ocean Carbon-Cycle Model Intercomparison Project. They have been merged and updated to reflect improvements concerning gas exchange, carbonate chemistry, and new data for initial conditions and atmospheric gas histories. Code is provided to facilitate their implementation.

Constraints on the Mediterranean Sea's storage of anthropogenic CO2 are limited, coming only from data-based approaches that disagree by more than a factor of two. Here we simulate this marginal sea's anthropogenic carbon storage by applying a perturbation approach in a high-resolution regional model. Our model simulates that, between 1800 and 2001, basin-wide CO2 storage by the Mediterranean Sea has increased by 1.0 Pg C, a lower limit based on the model's weak deep-water ventilation, as revealed by evaluation with CFC-12. Furthermore, by testing a data-based approach (transit time distribution) in our model, comparing simulated anthropogenic CO2 to values computed from simulated CFC-12 and physical variables, we conclude that the associated basin-wide storage of 1.7 Pg, published previously, must be an upper bound. Out of the total simulated storage of 1.0 Pg C, 75% comes from the air–sea flux into the Mediterranean Sea and 25% comes from net transport from the Atlantic across the Strait of Gibraltar. Sensitivity tests indicate that the Mediterranean Sea's higher total alkalinity, relative to the global-ocean mean, enhances the Mediterranean's total inventory of anthropogenic carbon by 10%. Yet the corresponding average anthropogenic change in surface pH does not differ significantly from the global-ocean average, despite higher total alkalinity. In Mediterranean deep waters, the pH change is estimated to be between −0.005 and −0.06 pH units.

²³⁰Th normalization is a valuable paleoceanographic tool for reconstructing high‐resolution sediment fluxes during the late Pleistocene (last ~500,000 years). As its application has expanded to ever more diverse marine environments, the nuances of²³⁰Th systematics, with regard to particle type, particle size, lateral advective/diffusive redistribution, and other processes, have emerged. We synthesized over 1000 sedimentary records of²³⁰Th from across the global ocean at two time slices, the late Holocene (0–5,000 years ago, or 0–5 ka) and the Last Glacial Maximum (18.5–23.5 ka), and investigated the spatial structure of²³⁰Th‐normalized mass fluxes. On a global scale, sedimentary mass fluxes were significantly higher during the Last Glacial Maximum (1.79–2.17 g/cm²kyr, 95% confidence) relative to the Holocene (1.48–1.68 g/cm²kyr, 95% confidence). We then examined the potential confounding influences of boundary scavenging, nepheloid layers, hydrothermal scavenging, size‐dependent sediment fractionation, and carbonate dissolution on the efficacy of²³⁰Th as a constant flux proxy. Anomalous²³⁰Th behavior is sometimes observed proximal to hydrothermal ridges and in continental margins where high particle fluxes and steep continental slopes can lead to the combined effects of boundary scavenging and nepheloid interference. Notwithstanding these limitations, we found that²³⁰Th normalization is a robust tool for determining sediment mass accumulation rates in the majority of pelagic marine settings (>1,000 m water depth).

Dissolved manganese (Mn) is a biologically essential element. Moreover, its oxidised form is involved in removing itself and several other trace elements from ocean waters. Here we report the longest thus far (17 500 km length) full-depth ocean section of dissolved Mn in the west Atlantic Ocean, comprising 1320 data values of high accuracy. This is the GA02 transect that is part of the GEOTRACES programme, which aims to understand trace element distributions. The goal of this study is to combine these new observations with new, state-of-the-art, modelling to give a first assessment of the main sources and redistribution of Mn throughout the ocean. To this end, we simulate the distribution of dissolved Mn using a global-scale circulation model. This first model includes simple parameterisations to account for the sources, processes and sinks of Mn in the ocean. Oxidation and (photo)reduction, aggregation and settling, as well as biological uptake and remineralisation by plankton are included in the model. Our model provides, together with the observations, the following insights: – The high surface concentrations of manganese are caused by the combination of photoreduction and sources contributing to the upper ocean. The most important sources are sediments, dust, and, more locally, rivers. – Observations and model simulations suggest that surface Mn in the Atlantic Ocean moves downwards into the southward-flowing North Atlantic Deep Water (NADW), but because of strong removal rates there is no elevated concentration of Mn visible any more in the NADW south of 40° N. – The model predicts lower dissolved Mn in surface waters of the Pacific Ocean than the observed concentrations. The intense oxygen minimum zone (OMZ) in subsurface waters is deemed to be a major source of dissolved Mn also mixing upwards into surface waters, but the OMZ is not well represented by the model. Improved high-resolution simulation of the OMZ may solve this problem. – There is a mainly homogeneous background concentration of dissolved Mn of about 0.10–0.15 nM throughout most of the deep ocean. The model reproduces this by means of a threshold on particulate manganese oxides of 25 pM, suggesting that a minimal concentration of particulate Mn is needed before aggregation and removal become efficient. – The observed distinct hydrothermal signals are produced by assuming both a strong source and a strong removal of Mn near hydrothermal vents.

The marine biological carbon pump is dominated by the vertical transfer of particulate organic carbon (POC) from the surface ocean to its interior. The efficiency of this transfer plays an important role in controlling the amount of atmospheric carbon that is sequestered in the ocean. Furthermore, the abundance and composition of POC is critical for the removal of numerous trace elements by scavenging, a number of which, such as iron, are essential for the growth of marine organisms, including phytoplankton. Observations and laboratory experiments have shown that POC is composed of numerous organic compounds that can have very different reactivities. However, this variable reactivity of POC has never been extensively considered, especially in modelling studies. Here, we introduced in the global ocean biogeochemical model NEMO-PISCES a description of the variable composition of POC based on the theoretical reactivity continuum model proposed by Boudreau and Ruddick (1991). Our model experiments show that accounting for a variable lability of POC increases POC concentrations in the ocean's interior by 1 to 2 orders of magnitude. This increase is mainly the consequence of a better preservation of small particles that sink slowly from the surface. Comparison with observations is significantly improved both in abundance and in size distribution. Furthermore, the amount of carbon that reaches the sediments is increased by more than a factor of 2, which is in better agreement with global estimates of the sediment oxygen demand. The impact on the major macronutrients (nitrate and phosphate) remains modest. However, iron (Fe) distribution is strongly altered, especially in the upper mesopelagic zone as a result of more intense scavenging: vertical gradients in Fe are milder in the upper ocean, which appears to be closer to observations. Thus, our study shows that the variable lability of POC can play a critical role in the marine biogeochemical cycles which advocates for more dedicated in situ and laboratory experiments.

Daily modeled fields of phosphate deposition to the Mediterranean from natural dust, anthropogenic combustion and wildfires were used to assess the effect of this external nutrient on marine biogeochemistry. The ocean model used is a high-resolution (1∕12°) regional coupled dynamical–biogeochemical model of the Mediterranean Sea (NEMO-MED12/PISCES). The input fields of phosphorus are for 2005, which are the only available daily resolved deposition fields from the global atmospheric chemical transport model LMDz-INCA. Traditionally, dust has been suggested to be the main atmospheric source of phosphorus, but the LMDz-INCA model suggests that combustion is dominant over natural dust as an atmospheric source of phosphate (PO₄, the bioavailable form of phosphorus in seawater) for the Mediterranean Sea. According to the atmospheric transport model, phosphate deposition from combustion (Pcomb) brings on average 40.5×10⁻⁶ mol PO₄ m⁻² yr⁻¹ over the entire Mediterranean Sea for the year 2005 and is the primary source over the northern part (e.g., 101×10⁻⁶ mol PO₄ m⁻² yr⁻¹ from combustion deposited in 2005 over the north Adriatic against 12.4×10⁻⁶ from dust). Lithogenic dust brings 17.2×10⁻⁶ mol PO₄ m⁻² yr⁻¹ on average over the Mediterranean Sea in 2005 and is the primary source of atmospheric phosphate to the southern Mediterranean Basin in our simulations (e.g., 31.8×10⁻⁶ mol PO₄ m⁻² yr⁻¹ from dust deposited in 2005 on average over the south Ionian basin against 12.4×10⁻⁶ from combustion). The evaluation of monthly averaged deposition flux variability of Pdust and Pcomb for the 1997–2012 period indicates that these conclusions may hold true for different years. We examine separately the two atmospheric phosphate sources and their respective flux variability and evaluate their impacts on marine surface biogeochemistry (phosphate concentration, chlorophyll a, primary production). The impacts of the different phosphate deposition sources on the biogeochemistry of the Mediterranean are found localized, seasonally varying and small, but yet statistically significant. Differences in the geographical deposition patterns between phosphate from dust and from combustion will cause contrasted and significant changes in the biogeochemistry of the basin. We contrast the effects of combustion in the northern basin (Pcomb deposition effects are found to be 10 times more important in the northern Adriatic, close to the main source region) to the effects of dust in the southern basin. These different phosphorus sources should therefore be accounted for in modeling studies.