Description: 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, SF6) 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, CTabio and 14CTabio, 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 CO2 as well as carbon isotope ratios of CO2. 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.

Description: Iron limits phytoplankton growth and hence the biological carbon pump in the Southern Ocean1. Models assessing the impacts of iron on the global carbon cycle generally rely on dust input and sediment resuspension as the predominant sources2, 3. Although it was previously thought that most iron from deep-ocean hydrothermal activity was inaccessible to phytoplankton because of the formation of particulates4, it has been suggested that iron from hydrothermal activity5, 6, 7 may be an important source of oceanic dissolved iron8, 9, 10, 11, 12, 13. Here we use a global ocean model to assess the impacts of an annual dissolved iron flux of approximately 9×108 mol, as estimated from regional observations of hydrothermal activity11, 12, on the dissolved iron inventory of the world’s oceans. We find the response to the input of hydrothermal dissolved iron is greatest in the Southern Hemisphere oceans. In particular, observations of the distribution of dissolved iron in the Southern Ocean3 (Chever et al., manuscript in preparation; Bowie et al., manuscript in preparation) can be replicated in our simulations only when our estimated iron flux from hydrothermal sources is included. As the hydrothermal flux of iron is relatively constant over millennial timescales14, we propose that hydrothermal activity can buffer the oceanic dissolved iron inventory against shorter-term fluctuations in dust deposition.

Description: A high-resolution dynamical model (Nucleus for European Modelling of the Ocean, Mediterranean configuration – NEMO-MED12) was used to give the first simulation of the distribution of radiocarbon (14C) across the whole Mediterranean Sea. The simulation provides a descriptive overview of both the natural pre-bomb 14C and the entire anthropogenic radiocarbon transient generated by the atmospheric bomb tests performed in the 1950s and early 1960s. The simulation was run until 2011 to give the post-bomb distribution. The results are compared to available in situ measurements and proxy-based reconstructions. The radiocarbon simulation allows an additional and independent test of the dynamical model, NEMO-MED12, and its performance to produce the thermohaline circulation and deep-water ventilation. The model produces a generally realistic distribution of radiocarbon when compared with available in situ data. The results demonstrate the major influence of the flux of Atlantic water through the Strait of Gibraltar on the inter-basin natural radiocarbon distribution and characterize the ventilation of intermediate and deep water especially through the propagation of the anthropogenic radiocarbon signal. We explored the impact of the interannual variability on the radiocarbon distribution during the Eastern Mediterranean Transient (EMT) event. It reveals a significant increase in 14C concentration (by more than 60 ‰) in the Aegean deep water and at an intermediate level (value up to 10 ‰) in the western basin. The model shows that the EMT makes a major contribution to the accumulation of radiocarbon in the eastern Mediterranean deep waters.

Description: The distribution of dissolved aluminium in the West Atlantic Ocean shows a mirror image with that of dissolved silicic acid, hinting at intricate interactions between the ocean cycling of Al and Si. The marine biogeochemistry of Al is of interest because of its potential impact on diatom opal remineralisation, hence Si availability. Furthermore, the dissolved Al concentration at the surface ocean has been used as a tracer for dust input, dust being the most important source of the bio-essential trace element iron to the ocean. Previously, the dissolved concentration of Al was simulated reasonably well with only a dust source, and scavenging by adsorption on settling biogenic debris as the only removal process. Here we explore the impacts of (i) a sediment source of Al in the Northern Hemisphere (especially north of ~ 40° N), (ii) the imposed velocity field, and (iii) biological incorporation of Al on the modelled Al distribution in the ocean. The sediment source clearly improves the model results, and using a different velocity field shows the importance of advection on the simulated Al distribution. Biological incorporation appears to be a potentially important removal process. However, conclusive independent data to constrain the Al / Si incorporation ratio by growing diatoms are missing. Therefore, this study does not provide a definitive answer to the question of the relative importance of Al removal by incorporation compared to removal by adsorptive scavenging.

Description: This archive contains all output from the PISCES-v2 model, containing 24 carbon, nitrate, phosphate and iron tracers. The model is extended with models of lithogenic particles and thorium-230 and protactinium-231, adding eight additional tracers, namely small and big lithogenic particles, and dissolved, small and big Th-230 and Pa-231 (referred to as ProThorP). Thus the model contains 32 different tracers (ptrc) for two different simulations, as well as diagnostic data (diad). The model is executed as part of the global ocean general cirulation model NEMO in the ORCA2 configuration (2° x 2° cos(phi) x 31 layers).

Description: 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 (17500 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.

Description: In this dataset we present a global compilation of over 1000 sedimentary records of 230Th from across the global ocean at two time slices, the Late Holocene (0-5000 years ago, or 0-5 ka) and the Last Glacial Maximum (18.5-23.5 ka). Data have been screened for age control, errors, and lithogenic corrections. Overall quality levels were computed by summing each record's scores on the individual criteria. A record is optimal if it is based on a chronology that is constrained by δ18O or 14C and it provides both the raw nuclide concentrations and the associated errors. About one quarter of the records in the database achieved this highest quality level. The large majority of the records in the database are good, passing two of the three criteria, while the remaining quarter are of fair or poor quality.

Description: 230Th‐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 230Th systematics, with regards to particle type, particle size, lateral advective/diffusive redistribution, and other processes, have emerged. We synthesized over 1000 sedimentary records of 230Th from across the global ocean at two time slices, the Late Holocene (0‐5000 years ago, or 0‐5 ka) and the Last Glacial Maximum (18.5‐23.5 ka), and investigated the spatial structure of 230Th‐normalized mass fluxes. On a global scale, sedimentary mass fluxes were significantly higher during the Last Glacial Maximum (1.79‐2.17 g/cm2kyr, 95% confidence) relative to the Holocene (1.48‐1.68 g/cm2kyr, 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 230Th as a constant flux proxy. Anomalous 230Th 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 230Th‐normalization is a robust tool for determining sediment mass accumulation rates in the majority of pelagic marine settings (〉 1000 m water depth).

Description: 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, SF6) 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, Cabio T and 14Cabio T , 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 CO2 as well as carbon isotope ratios of CO2. 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.