Iron is a key micronutrient controlling phytoplankton growth in vast regions of the global ocean. Despite its importance, uncertainties remain high regarding external iron source fluxes and internal cycling on a global scale. In this study, we used a global dissolved iron dataset, including GEOTRACES measurements, to constrain source and scavenging fluxes in the marine iron component of a global ocean biogeochemical model. Our model simulations tested three key uncertainties: source inputs of atmospheric soluble iron deposition (varying from 1.4–3.4 Gmol/yr), reductive sedimentary iron release (14–117 Gmol/yr), and compared a variable ligand parameterization to a constant distribution. In each simulation, scavenging rates were tuned to reproduce the observed global mean iron inventory for consistency. The variable ligand parameterization improved the global model-data misfit the most, suggesting that heterotrophic bacteria are an important source of ligands to the ocean. Model simulations containing high source fluxes of atmospheric soluble iron deposition (3.4 Gmol/yr) and reductive sedimentary iron release (114 Gmol/yr) further improved the model most notably in the surface ocean. High scavenging rates were then required to maintain the iron inventory resulting in relatively short surface and global ocean residence times of 0.83 and 7.5 years, respectively. The model simulates a tight spatial coupling between source inputs and scavenging rates, which may be too strong due to underrepresented ligands near source inputs, contributing to large uncertainties when constraining individual fluxes with dissolved iron concentrations. Model biases remain high and are discussed to help improve global marine iron cycle models.

Biological productivity in the ocean directly influences the partitioning of carbon between the atmosphere and ocean interior. Through this carbon cycle feedback, changing ocean productivity has long been hypothesized as a key pathway for modulating past atmospheric carbon dioxide levels and hence global climate. Because phytoplankton preferentially assimilate the light isotopes of carbon and the major nutrients nitrate and silicic acid, stable isotopes of carbon (C), nitrogen (N), and silicon (Si) in seawater and marine sediments can inform on ocean carbon and nutrient cycling, and by extension the relationship with biological productivity and global climate. Here we compile water column C, N, and Si stable isotopes from GEOTRACES-era data in four key ocean regions to review geochemical proxies of oceanic carbon and nutrient cycling based on the C, N, and Si isotopic composition of marine sediments. External sources and sinks as well as internal cycling (including assimilation, particulate matter export, and regeneration) are discussed as likely drivers of observed C, N, and Si isotope distributions in the ocean. The potential for C, N, and Si isotope measurements in sedimentary archives to record aspects of past ocean C and nutrient cycling is evaluated, along with key uncertainties and limitations associated with each proxy. Constraints on ocean C and nutrient cycling during late Quaternary glacial-interglacial cycles and over the Cenozoic are examined. This review highlights opportunities for future research using multielement stable isotope proxy applications and emphasizes the importance of such applications to reconstructing past changes in the oceans and climate system.

Phytoplankton productivity and export sequester climatically significant quantities of atmospheric carbon dioxide as particulate organic carbon through a suite of processes termed the biological pump. How the biological pump operated in the past is therefore important for understanding past atmospheric carbon dioxide concentrations and Earth’s climate history. However, reconstructing the history of the biological pump requires proxies. Due to their intimate association with biological processes, several bioactive trace metals and their isotopes are potential proxies for past phytoplankton productivity, including: iron, zinc, copper, cadmium, molybdenum, barium, nickel, chromium, and silver. Here we review the oceanic distributions, driving processes, and depositional archives for these nine metals and their isotopes based on GEOTRACES-era datasets. We offer an assessment of the overall maturity of each isotope system to serve as a proxy for diagnosing aspects of past ocean productivity and identify priorities for future research. This assessment reveals that cadmium, barium, nickel, and chromium isotopes offer the most promise as tracers of paleoproductivity, whereas iron, zinc, copper, and molybdenum do not. Too little is known about silver to make a confident determination. Intriguingly, the elements that are least sensitive to productivity may be used to trace other aspects of ocean chemistry, such as nutrient sources, particle scavenging, organic complexation, and ocean redox state. These complementary sensitivities suggest new opportunities for combining perspectives from multiple proxies that will ultimately enable painting a more complete picture of marine paleoproductivity, biogeochemical cycles, and Earth’s climate history.

Polar marine ecosystems are particularly vulnerable to the effects of climate change. Warming temperatures, freshening seawater and disruption to sea ice formation potentially all have detrimental cascading effects on food webs. New approaches are needed to better understand spatio-temporal interactions among biogeochemical processes at the base of Southern Ocean food webs, and how these interactions vary seasonally. In marine systems, isoscapes (models of the spatial variation in the stable isotopic composition) of carbon and nitrogen identify the spatial expression of varying biogeochemical processes on nutrient utilization by phytoplankton. Isoscapes also provide a baseline for interpreting stable isotope compositions of higher trophic level animals in movement, migration and diet research. Here we produce carbon and nitrogen isoscapes across the entire Southern Ocean (>40°S) using surface particulate organic matter (POM) isotope data, collected from multiple sources over the past 50 years and throughout the annual cycle. We use Integrated Nested Laplace Approximation (INLA)-based approaches to predict mean annual isoscapes and four seasonal isoscapes using a suite of environmental data as predictor variables. Clear spatial gradients in δ13C and δ15N values were predicted across the Southern Ocean, consistent with previous statistical and mechanistic isoscape views of isotopic variability in this region. We identify strong seasonal variability in both carbon and nitrogen isoscapes, with key implications for the use of static or annual average isoscape baselines in animal studies attempting to document seasonal migratory or foraging behaviours.