The riddle of eastern tropical Pacific Ocean oxygen levels: the role of the supply by intermediate-depth waters

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

Simulated effects of southern hemispheric wind changes on the Pacific oxygen minimum zone

A coupled ocean biogeochemistry-circulation model is used to investigate the impact of observed past and anticipated future wind changes in the southern hemisphere on the oxygen minimum zone in the tropical Pacific. We consider the industrial period until the end of the 21st century and distinguish effects due to a strengthening of the westerlies from effects of a southward shift of the westerlies that is accompanied by a poleward expansion of the tropical trade winds. Our model results show that a strengthening of the westerlies counteracts part of the warming-induced decline in the global marine oxygen inventory. A poleward shift of the trade-westerlies boundary, however, triggers a significant decrease of oxygen in the tropical oxygen minimum zone. In a business-as-usual CO2 emission scenario, the poleward shift of the trade-westerlies boundary and warming-induced increase in stratification contribute equally to the expansion of suboxic waters in the tropical Pacific

MOMSO 1.0 – an eddying Southern Ocean model configuration with fairly equilibrated natural carbon

We present a new near-global coupled biogeochemical ocean-circulation model configuration. The configuration features a horizontal discretization with a grid spacing of less than 11 km in the Southern Ocean and gradually coarsens in meridional direction to more than 200 km at 64∘ N, where the model is bounded by a solid wall. The underlying code framework is the Geophysical Fluid Dynamics Laboratory (GFDL)'s Modular Ocean Model coupled to the Biogeochemistry with Light, Iron, Nutrients and Gases (BLING) ecosystem model of Galbraith et al. (2010). The configuration is unique in that it features both a relatively equilibrated oceanic carbon inventory and an eddying ocean circulation based on a realistic model geometry/bathymetry – a combination that has been precluded by prohibitive computational cost in the past. Results from a simulation with climatological forcing and a sensitivity experiment with increasing winds suggest that the configuration is sufficiently equilibrated to explore Southern Ocean carbon uptake dynamics on decadal timescales. The configuration is dubbed MOMSO, a Modular Ocean Model Southern Ocean configuration

Effects of increased isopycnal diffusivity mimicking the unresolved equatorial intermediate current system in an earth system climate model

Earth system climate models generally underestimate dissolved oxygen concentrations in the deep eastern equatorial Pacific. This problem is associated with the "nutrient trapping" problem, described by Najjar et al. [1992], and is, at least partially, caused by a deficient representation of the Equatorial Intermediate Current System (EICS). Here we emulate the unresolved EICS in the UVic earth system climate model by locally increasing the zonal isopycnal diffusivity. An anisotropic diffusivity of ∼50,000 m 2 s-1 yields an improved global representation of temperature, salinity and oxygen. In addition, it (1) resolves most of the local "nutrient trapping" and associated oxygen deficit in the eastern equatorial Pacific and (2) reduces spurious zonal temperature gradients on isopycnals without affecting other physical metrics such as meridional overturning or air-sea heat fluxes. Finally, climate projections of low-oxygenated waters and associated denitrification change sign and apparently become more plausibl

Meeting climate targets by direct CO2 injections: what price would the ocean have to pay?

We investigate the climate mitigation potential and collateral effects of direct injections of captured CO2 into the deep ocean as a possible means to close the gap between an intermediate CO2 emissions scenario and a specific temperature target, such as the 1.5 ∘C target aimed for by the Paris Agreement. For that purpose, a suite of approaches for controlling the amount of direct CO2 injections at 3000 m water depth are implemented in an Earth system model of intermediate complexity. Following the representative concentration pathway RCP4.5, which is a medium mitigation CO2 emissions scenario, cumulative CO2 injections required to meet the 1.5 ∘C climate goal are found to be 390 Gt C by the year 2100 and 1562 Gt C at the end of simulations, by the year 3020. The latter includes a cumulative leakage of 602 Gt C that needs to be reinjected in order to sustain the targeted global mean temperature. CaCO3 sediment and weathering feedbacks reduce the required CO2 injections that comply with the 1.5 ∘C target by about 13 % in 2100 and by about 11 % at the end of the simulation. With respect to the injection-related impacts we find that average pH values in the surface ocean are increased by about 0.13 to 0.18 units, when compared to the control run. In the model, this results in significant increases in potential coral reef habitats, i.e., the volume of the global upper ocean (0 to 130 m depth) with omega aragonite > 3.4 and ocean temperatures between 21 and 28 ∘C, compared to the control run. The potential benefits in the upper ocean come at the expense of strongly acidified water masses at depth, with maximum pH reductions of about −2.37 units, relative to preindustrial levels, in the vicinity of the injection sites. Overall, this study demonstrates that massive amounts of CO2 would need to be injected into the deep ocean in order to reach and maintain the 1.5 ∘C climate target in a medium mitigation scenario on a millennium timescale, and that there is a trade-off between injection-related reductions in atmospheric CO2 levels accompanied by reduced upper-ocean acidification and adverse effects on deep-ocean chemistry, particularly near the injection sites

Patterns of deoxygenation: sensitivity to natural and anthropogenic drivers

Observational estimates and numerical models both indicate a significant overall decline in marine oxygen levels over the past few decades. Spatial patterns of oxygen change, however, differ considerably between observed and modelled estimates. Particularly in the tropical thermocline that hosts open-ocean oxygen minimum zones, observations indicate a general oxygen decline, whereas most of the state-of-the-art models simulate increasing oxygen levels. Possible reasons for the apparent model-data discrepancies are examined. In order to attribute observed historical variations in oxygen levels, we here study mechanisms of changes in oxygen supply and consumption with sensitivity model simulations. Specifically, the role of equatorial jets, of lateral and diapycnal mixing processes, of changes in the wind-driven circulation and atmospheric nutrient supply, and of some poorly constrained biogeochemical processes are investigated. Predominantly wind-driven changes in the low-latitude oceanic ventilation are identified as a possible factor contributing to observed oxygen changes in the low-latitude thermocline during the past decades, while the potential role of biogeochemical processes remains difficult to constrain. We discuss implications for the attribution of observed oxygen changes to anthropogenic impacts and research priorities that may help to improve our mechanistic understanding of oxygen changes and the quality of projections into a changing future

Pilot Study on Potential Impacts of Fisheries-Induced Changes in Zooplankton Mortality on Marine Biogeochemistry

In this pilot study we link the yield of industrial fisheries to changes in the zooplankton mortality in an idealized way accounting for different target species (planktivorous fish—decreased zooplankton mortality; large predators—increased zooplankton mortality). This indirect approach is used in a global coupled biogeochemistry circulation model to estimate the range of the potential impact of industrial fisheries on marine biogeochemistry. The simulated globally integrated response on phytoplankton and primary production is in line with expectations—a high (low) zooplankton mortality results in a decrease (increase) of zooplankton and an increase (decrease) of phytoplankton. In contrast, the local response of zooplankton and phytoplankton depends on the region under consideration: In nutrient-limited regions, an increase (decrease) in zooplankton mortality leads to a decrease (increase) in both zooplankton and phytoplankton biomass. In contrast, in nutrient-replete regions, such as upwelling regions, we find an opposing response: an increase (decrease) of the zooplankton mortality leads to an increase (decrease) in both zooplankton and phytoplankton biomass. The results are further evaluated by relating the potential fisheries-induced changes in zooplankton mortality to those driven by CO2 emissions in a business-as-usual 21st century emission scenario. In our idealized case, the potential fisheries-induced impact can be of similar size as warming-induced changes in marine biogeochemistry

Diagnostics of diapycnal diffusion in z-level ocean models

In general ocean circulation models (OGCMs) diapycnal diffusion arises notonly from the discretisation of the explicit diffusion, but also by numericallyinduced diffusion, caused, e.g., by common discretisations of advectivetransport.In the present study, three different diagnostics to analyse the meandiapycnal diffusivities of individual tracers (vertically and horizontally) areintroduced: (i) The divergence method based on the work of Ledwellet al. (1998) infers the mean diapycnal diffusivity from theadvection-diffusion equation. (ii) The tracer flux method based on the workof Griffies et al. (2000), that determines the diapycnal flux crossing anisopycnal layer, is modified for the analysis of mean diapycnal diffusivitiesof a passive tracer. (iii) The variance method based on the work ofMorales Maqueda and Holloway (2006) is a more general approachas the diapycnal diffusion is analysed by the variance decay of the totaltracer concentration.These methods can be used for the analysis of the diffusivity of passivetracer independent of the model set-up, e.g. the advection scheme used, butsupport only information about mean diapycnal diffusivity of that tracerfield rather than for each individual layer. The applicability of thesemethods is tested in a set of 1- and 2-dimensional case studies. The effectof vertical advection and of diverging and converging isopycnals is shownseparately. In all three methods used, the transformation of the tracer ontoisopycnals leads to errors in the diagnosed diffusivities. It turns out thatthe tracer flux method is the most robust method and therefore the methodof choice. In order to keep the errors as small as possible, longer time meanvalues should be analysed