Marine N₂O cycling from high spatial resolution concentration, stable isotopic and isotopomer measurements along a meridional transect in the eastern Pacific Ocean

Abstract

Nitrous oxide (N₂O) is a potent greenhouse gas and ozone depleting substance, with the ocean accounting for about one third of global emissions. In marine environments, a significant amount of N₂O is produced by biological processes in Oxygen Deficient Zones (ODZs). While recent technological advances are making surface N₂O concentration more available, high temporal and spatial resolution water-column N₂O concentration data are relatively scarce, limiting global N₂O ocean models’ predictive capability. We present a N₂O concentration, stable isotopic composition and isotopomer dataset of unprecedently large spatial coverage and depth resolution in the broader Pacific, crossing both the eastern tropical South and North Pacific Ocean ODZs collected as part of the GO-SHIP P18 repeat hydrography program in 2016/2017. We complement these data with dissolved gases (nitrogen, oxygen, argon) and nitrate isotope data to investigate the pathways controlling N₂O production in relation to apparent oxygen utilization and fixed nitrogen loss. N₂O yield significantly increased under low oxygen conditions near the ODZs. Keeling plot analysis revealed different N₂O sources above the ODZs under different oxygen regimes. Our stable isotopic data and relationships between the N₂O added by microbial processes (ΔN₂O) and dissolved inorganic nitrogen (DIN) deficit confirm increased N₂O production by denitrification under low oxygen conditions near the oxycline where the largest N₂O accumulations were observed. The slope for δ¹⁸O-N₂O versus site preference (SP, the difference between the central (α) and outer (β) N atoms in the linear N₂O molecule) in the eastern tropical North Pacific ODZ was lower than expected for pure N₂O reduction, likely because of the observed decrease in δ¹⁵N^β. This trend is consistent with prior ODZ studies and attributed to concurrent production of N₂O from nitrite with a low δ¹⁵N or denitrification with a SP >0‰. We estimated apparent isotope effects for N₂O consumption in the ETNP ODZ of 3.6‰ for ¹⁵N^bulk, 9.4‰ for ¹⁵N^α, -2.3‰ for ¹⁵N^β, 12.0‰ for ¹⁸O, and 11.7‰ for SP. These values were generally within ranges previously reported for previous laboratory and field experiments.

1 Introduction

Nitrous oxide (N₂O) is mainly produced by microbial processes and contributes to climate change as a tropospheric greenhouse gas approximately 275 times more potent than CO₂ on a per molecule basis (Forster et al., 2021). N₂O produced at the Earth’s surface has a long atmospheric residence time of more than 100 years and is the main source of ozone-depleting nitric-oxide radicals in the stratosphere (Nevison and Holland, 1997; Ravishankara et al., 2009). Biogeochemical models estimate that the ocean accounts for about one third of global N₂O emissions (Freing et al., 2012; Ciais et al., 2013; Tian et al., 2020; Forster et al., 2021) with a global oceanic flux estimate of 4.2 ± 1.0 Tg N y^-1 (Yang et al., 2020). Coastal upwelling systems associated with Oxygen Deficient Zones (ODZs) account for up to 35% of the oceanic N₂O total emissions (Bakker et al., 2014; Arévalo-Martínez et al., 2015; Yang et al., 2020), yet temporal variability from these productive coastal waters is still not well quantified.

N₂O is produced in oxic ocean waters as a by-product of nitrification through the oxidation of hydroxylamine (NH₂OH), an intermediate during ammonia (NH₃) oxidation to nitrite (NO₂^-) by both archaeal and bacterial nitrifiers (Caranto and Lancaster, 2017) (Figure 1). Under low oxygen (O₂)-conditions, ammonia oxidizing bacteria (AOB) produce N₂O by nitrifier-denitrification, the sequential NO₂^- reduction to N₂O during respiration (Goreau et al., 1980; Wrage et al., 2001; Frame and Casciotti, 2010). N₂O is also produced by ammonia oxidizing archaea (AOA) following a hybrid pathway where one atom in the N₂O molecule is derived from NH₃ (e.g., NH₂OH) and the other from NO₂^- (Trimmer et al., 2016). N₂O production by nitrification leads to strong positive correlations between Apparent Oxygen Utilization (AOU) and ΔN₂O (i.e., the difference between N₂O measured and at atmospheric equilibrium) and nitrate (NO₃^-) concentrations (Yoshinari and Knowles, 1976; Cohen and Gordon, 1979; Nevison et al., 2003). N₂O production yield per NO₂^- generated during nitrification by either AOA or AOB is generally low, varying from 0 to 2% of NO₃^- production (Yoshida et al., 1989; Frame and Casciotti, 2010; Santoro et al., 2010; Santoro et al., 2011; Löscher et al., 2012; Ryabenko et al., 2012) but is generally higher for AOB (Hink et al., 2017a; Hink et al., 2017b) and enhanced under low-O₂ conditions according to both culture and field observations (up to 10% at low O₂; Goreau et al., 1980; Ji et al., 2015).

Figure 1

Bourbonnais et al. (2021).

Under anoxic conditions, denitrification, the sequential reduction of NO₃^-, NO₂^-, and NO to N₂O, both produces and consumes N₂O. Enhanced N₂O production occurs under low-O₂ conditions (generally below 10 µmol kg^-1; Frey et al., 2020) due to different O₂ sensitivities for the different steps involved in denitrification, with nitrous oxide reductase being inhibited by even nanomolar O₂ concentrations (Dalsgaard et al., 2014). N₂O produced by incomplete denitrification significantly accumulates in low-O₂ conditions near the upper oxycline in the eastern tropical North Pacific (ETNP) and eastern tropical South Pacific (ETSP) ODZs (e.g., Ji et al., 2015; Kock et al., 2016; Bourbonnais et al., 2017; Casciotti et al., 2018; Frey et al., 2020; Kelly et al., 2021; Monreal et al., 2022). High N₂O accumulations of up to 1 μmol L^-1 were observed in surface waters off the Peru coast, where anoxic waters from the ETSP ODZ are continuously upwelled (Arévalo-Martínez et al., 2015). N₂O accumulation is caused by a decoupling between N₂O production and its reduction to N₂ by denitrification, the latter process being less O₂ tolerant (Dalsgaard et al., 2014). N₂O is generally completely reduced to non-bioavailable N₂ in anoxic waters. ODZs are thus generally net N₂O sinks, with concentrations near or below atmospheric equilibrium concentrations (Bange et al., 2001; Yamagishi et al., 2007; Kock et al., 2016).

Natural stable isotopes are broadly used as tracers of N-cycle processes in the ocean which integrate over space and time (e.g., Sigman et al., 2005; Altabet, 2006; Bourbonnais et al., 2009; Bourbonnais et al., 2015). The N₂O molecule contains both bulk (δ¹⁵N and δ¹⁸O) and, given its asymmetric arrangement of atoms, ¹⁵N site specific signatures that are valuable for identifying production and consumption processes. Bulk δ¹⁵N and δ¹⁸O are expressed as:

Units are in parts per thousand or per mil (‰) and R is the ratio of ¹⁵N/¹⁴N or ¹⁸O/¹⁶O. Reference materials are atmospheric N₂ for N (scale AIR-N₂) and mean ocean water for O (scale Vienna Standard Mean Ocean Water, V-SMOW). The bulk isotopic composition (δ¹⁵N and δ¹⁸O) of N₂O depends in part on the isotopic composition of its substrates. For instance, for hydroxylamine oxidation, bulk N₂O δ¹⁵N and δ¹⁸O depend on the δ¹⁵N of the source NH₃ and δ¹⁸O of dissolved O₂. Conversely for nitrifier-denitrification and denitrification, N₂O δ¹⁵N and δ¹⁸O is dependent on the δ¹⁵N and δ¹⁸O of source NO₃^- and/or NO₂^- (Ostrom and Ostrom, 2012). Additionally, significant O exchange usually occurs between NO₂^- and H₂O during N₂O production by nitrifier-denitrification or denitrification (Buchwald and Casciotti, 2010; Kool et al., 2011), which acts to decouple the δ¹⁸O values of source and product.

Isotopic fractionation during nitrification and denitrification is the other major influence on the δ¹⁵N and δ¹⁸O of N₂O. Kinetic isotope fractionation occurs as the molecules containing the lighter isotopes (e.g., ¹⁴N, ¹⁶O) react more quickly leaving the residual substrate enriched in heavier isotopes (e.g., ¹⁵N and ¹⁸O). The isotope effect (ε) is defined by:

where k₁ and k₂ are the specific reaction rates for the lighter and heavier isotope, respectively. N and O isotope effects (¹⁵ε, ¹⁸ε) during N₂O production and consumption vary substantially in laboratory culture as well as in the environment (Bourbonnais et al., 2017 and references therein).

In contrast to bulk isotope values, N₂O Site Preference (SP) is independent of initial isotopic composition of the substrate (Toyoda et al., 2002; Schmidt et al., 2004; Sutka et al., 2006). Thus, SP is generally process-dependent and can be used as a tracer to identify the source of N₂O. SP is calculated from the difference in δ¹⁵N between the central (α) and outer (β) N atoms in the linear, asymmetrical N₂O molecule (NNO):

Non-zero SP arises from the differential biochemical bond making and breaking experienced by each of the two N atoms as a consequence of their different molecular positions. Low SP isotopic signatures (-11 to 0‰) are associated with N₂O production via NO₂^- reduction by nitrifier-denitrification or denitrification. Much higher SP values are indicative of abiotic N₂O formation (Heil et al., 2014) or production by hydroxylamine oxidation (30-36‰) (Sutka et al., 2006; Frame and Casciotti, 2010). N₂O isotopic signatures have consequently been classified into distinct compositional fields used to evaluate the relative contribution from different processes (e.g., Wankel et al., 2017 and reference therein). However, SP does increase as a result of isotope fractionation during consumption by denitrification as discussed above (Yamagishi et al., 2005; Ostrom et al., 2007; Yamagishi et al., 2007).

The main objectives of this study are to better understand nitrous oxide (N₂O) distribution and production mechanisms along the meridional P18 transect in the eastern Pacific sampled as part of the U.S. GO-SHIP program (Figure 2). We sampled at an unprecedently high spatial resolution to investigate how N₂O dynamics respond to changes in O₂ concentrations and dissolved inorganic nitrogen (DIN) loss. We then evaluated pathways responsible for N₂O production in sub-oxic waters (O₂ > 5 μmol kg^-1) overlying the ETNP and ETSP ODZs as well as deep waters using stable isotopes and isotopomers and Keeling plot analysis. We also investigated N₂O cycling within the ETNP ODZ and estimated apparent isotope effects for N₂O consumption.

Figure 2

2 Methods

The U.S. GO-SHIP Program conducted a hydrographic survey along the P18 section in 2016/2017 (Figure 2). Core physical and biogeochemical data from the cruise (e.g., temperature and salinity, O₂ and nutrient concentrations, transient tracers, radiocarbon) are publicly available at https://cchdo.ucsd.edu/cruise/33RO20161119. Samples for N₂/Ar were collected every ~2 degrees latitude during leg 1 at 25 stations in total. N₂O stable isotope and isotopomer samples were collected at 44 stations, every ~2 degrees latitude during leg 1 and every 3 to 4 degrees during leg 2, except between 3°N and 3°S where the resolution was increased to ~1 degree. Samples for δ¹⁵N of NO₃^- analysis were collected during leg 1 at every degree of latitude, but only analyzed at station 32 in the ETNP and stations 53, 55, 62, 68, 74, 78, 90, 92, 96, 100, 104, 108 and 112 near the equator and the ETNP. We usually collected water column profiles from the surface to 2000 m depth, with deeper profiles at 3 stations.

2.1 N₂/Ar

N₂/Ar and δ¹⁵N₂ samples were collected and preserved as in Charoenpong et al. (2014) and analyzed at the University of Massachusetts Dartmouth (UMass Dartmouth). In total, approximately 600 samples were collected in 60 mL serum bottles. Water samples were pumped, at 5 to 10 ml min^-1, through a continuous sparger which transfers dissolved gases quantitatively to a continuous flow of helium (He) carrier gas. Analysis time was about 10 minutes. Carrier gas was passed through water, carbon dioxide (CO₂), and software selectable hot-copper O₂ traps before admittance via an open split to an isotope-ratio mass spectrometer (IRMS). The GV IsoPrime IRMS was fitted with collectors for simultaneous measurement of N₂ (masses 28 and 29), O₂ (masses 32, 33, and 34), and Ar (mass 40). Gas and isotopic ratios were measured against artificial compressed gas mixtures of N₂, O₂, and Ar close to expected dissolved gas ratios. Reproducibility of N₂/Ar measurements were better than 0.5‰. Daily calibration against an atmosphere reference (water equilibrated with air) was carried out at precisely controlled temperatures of 10.0 and 20.0°C. Excess (biogenic) N₂ was calculated against equilibrium values expected from in situ temperature and salinity as in Charoenpong et al. (2014). Subtraction of the background N₂ excess was done as in Bourbonnais et al. (2015) using N₂/Ar values in source waters near the equator.

2.2 N₂O concentrations, stable isotopes and isotopomers

Samples for dissolved N₂O were collected in a similar fashion as for dissolved O₂/N₂/Ar samples. Tygon tubing was attached to the Niskin bottle and a 165 mL serum glass bottle was filled and overflowed with seawater at least 2 times before capping with a butyl stopper and crimp sealed with aluminum. This procedure was executed underwater in a plastic container to avoid air bubbles. After collection, 0.2 mL of a saturated HgCl₂ solution was injected to prevent biological activity. In total, approximately 800 samples were collected.

Samples were analyzed at UMass Dartmouth using a GV IsoPrime Continuous Flow, MultiCollector, IRMS (CF-MC-IRMS) coupled to an automated gas extraction as described in Bourbonnais et al. (2017). The collector configuration included masses 30, 31 for the NO⁺ fragment of N₂O (used for the determination of δ¹⁵N^α) and 44, 45, and 46 (used for the determination of δ¹⁵N^bulk and δ¹⁸O). Briefly, dissolved N₂O was pumped through a gas extractor, completely extracted through bubbling with a continuous flow of He and concentrated and purified in a purge-trap system. CO₂ and H₂O vapor were removed using chemical and cryogenic traps. N₂O was cryo-focused with two liquid N₂ traps and passed through a capillary gas chromatography (GC) column prior to IRMS analysis. Results were reproducible even at low N₂O concentration (down to ~5 nmol L^-1). N₂O concentrations in our samples were calculated from relative peak heights between the samples and a seawater standard of known N₂O concentration equilibrated with seawater at 5°C (12.5 nmol L^-1 at salinity 34 as calculated using the Weiss and Price (1980) equation). Equilibrium N₂O concentrations at depth were calculated using historical reconstructions of atmospheric N₂O and water mass age distributions estimated from chlorofluorocarbon and sulfur hexafluoride tracers as described in section 2.4.

The N₂O concentrations measured with our IRMS agreed well with those measured independently using gas chromatography and an electron capture detector (ECD) at sea during the same research expedition (data available at https://cchdo.ucsd.edu/cruise/33RO20161119).

The reproducibility of δ¹⁵N^bulk, the average of δ¹⁵N^α and δ¹⁵N^β, δ¹⁸O and SP as well as any instrumental drift were determined from measurements of the 5°C seawater standard distributed throughout an analytical run. We calibrated measurements and corrected for scrambling between the α and β positions (Westley et al., 2007) using four N₂O standards covering a large range of SP (as well as δ¹⁵N^bulk and δ¹⁸O) composition (see Mohn et al., 2014). These standards were analyzed in duplicate for each run to quantify the scrambling effect and potential offsets, and we iteratively solved for the different calibration parameters as described in Frame and Casciotti (2010) and Mohn et al. (2014). Correction for isobaric interference from ¹⁷O was included in these procedures. Standard deviations for triplicate measurements of our N₂O standards were typically below 0.1‰ for δ¹⁵N^bulk N₂O, 0.1‰ for δ¹⁸O-N₂O and 1.0‰ for SP, which were comparable to values reported by Mohn et al. (2014).

2.3 Nitrate isotopes

The δ¹⁵N of NO₃^- was measured using the Ti (III) reduction method for nitrate conversion to N₂O (Altabet et al., 2019). The product N₂O was purified and analyzed at UMass Dartmouth using a GV Isoprime continuous flow IRMS equipped with a custom purge trap extraction system and a PAL autosampler. Target sample size was 10 nmoles. Samples were standardized using a three-point correction with the international standards IAEA-N3 (δ¹⁵N = 4.7‰ vs AIR), USGS-34 (δ¹⁵N = -1.8‰ vs AIR) and USGS-35 (2.7‰ vs AIR). The isotope ratios are expressed in δ (‰) relative to AIR for N. The average standard deviation for duplicate samples was better than ± 0.3‰ for δ¹⁵N.

2.4 Calculation of N₂O at equilibrium using transit time distributions

In the same water parcel as N₂O, three transient tracers [chlorofluorocarbon (CFC)-11, CFC-12, and sulfur hexafluoride (SF₆)] were measured, each of which are affected by mixing differently due to their differing atmospheric histories. This concurrence provides a means of estimating, from multiple tracer ages, the impact of mixing on the ages. To estimate mean ventilation timescales, age distributions, otherwise known as transit-time distributions (TTDs), were approximated using the inverse Gaussian form:

where Γ is the mean age and Δ is the width, or standard deviation, of the distribution (Waugh et al., 2003).

Each water sample’s measured CFC-11, CFC-12, and SF₆ were matched to a TTD lookup table (Stanley et al., 2012; Sonnerup et al., 2015) that allowed for Δ/Γ to range, in 0.1 increments, from 0.3 up to 1.8, the range that can be constrained by those three tracers (Stöven et al., 2015), and included each TTD’s atmospheric N₂O, below. The matching criterion was within ± 3 σ of the combined tracer measurement and tracer saturation level uncertainties. The saturation levels dominate the uncertainty. Saturation levels of CFC-11, CFC-12, and SF₆ were estimated from the outcropping region of the respective isopycnals.

2.4.1 N₂O atmospheric history

The N₂O atmospheric history post-1800 was estimated using the polynomial fit presented in Freing et al. (2009), updated to the present using annual global mean N₂O atmospheric mixing ratios from the NOAA Earth System Research Laboratory (ESRL) Global Monitoring Division (https://www.esrl.noaa.gov/gmd/hats/combined/N2O.html). Pre-1800, a constant N₂O mixing ratio of 274.81 ppb was used.

2.4.2 Using TTDs to constrain biogenic N₂O

The total N₂O in a water parcel (N₂O_obs) is the sum of N₂O from two different sources: 1) equilibrated from the atmosphere at the time it was last at the surface (N₂O_bkg) and 2) added from microbial processes (ΔN₂O or N₂O_prod):

where X_N2O is the atmospheric mixing ratio of N₂O, P is the atmospheric pressure, and F_T,S is the temperature and salinity dependent solubility of N₂O in seawater (Weiss and Price, 1980). The N₂O mixing ratio of the TTDs for each water parcel were calculated using the atmospheric history of N₂O. The TTDs’ N₂O mixing ratio were used to calculate the N₂O concentration at atmospheric equilibrium (N₂O_bkg) for each water sample. For values where the CFC-11 was below detection (< 0.001 pmol kg^-1), an atmospheric N₂O mixing ratio of 274.81 ppb was used.

2.5 Keeling plot analysis

In order to determine the sources of high ΔN₂O near the oxycline in the ETNP and ETSP ODZs, the isotopic and isotopomer compositions of the N₂O produced within given water masses and at different O₂ concentrations were estimated using Keeling plot analysis (Keeling, 1961; Pataki et al., 2003; Yamagishi et al., 2007; Fujii et al., 2013; Casciotti et al., 2018). The isotopic composition of source N₂O was determined using equation 5 and:

where the obs, prod and bkg subscripts refer to N₂O measured, newly produced and at equilibrium, respectively.

Rearranging equations 5 and 7 gives:

The intercepts of the linear regression between the inverse of measured N₂O concentration and the δ¹⁵N^bulk, δ¹⁵N^α, δ¹⁵N^β, δ¹⁸O and SP thus represents the isotopic composition of produced N₂O.

2.6 Isotope effects during N₂O consumption

We estimated isotope effects associated with N₂O consumption in the ETNP ODZ using a closed-system Rayleigh model (Mariotti et al., 1981). Only data in the ETNP ODZ at [O₂] < 5 µmol kg^-1 where N₂O consumption occurs were selected for this analysis.

where δ¹⁵N₂O is the δ¹⁵N^bulk, δ¹⁵N^α, δ¹⁵N^β, δ¹⁸O or SP for samples at O₂ < 5 µmol kg^-1 (mostly in the ETNP ODZ), and the subscript refers to N₂O measured and initial concentrations and isotopic values before the onset of N₂O consumption. The isotope effects were estimated as the slopes of the linear regressions between δN₂O versus -ln[N₂O], -ln[N₂O] being an approximation of the [N₂O]_obs/[N₂O]_initial term.

3 Results

3.1 Water mass characterization

Water masses were characterized according to Emery (2001). The main surface water (0 – 500 m depth) masses along the P18 transects, listed from north to south with potential temperature and absolute salinity ranges in brackets, were Eastern North Pacific Central Water (ENPCW; 12.0-20°C; 34.2-35.0), Eastern North Pacific Transition Water (ENPTW; 11.0-20.0°C, 33.3-34.3), Pacific Equatorial Water (PEW; 7.0-23.0°C, 34.5-36.0), Eastern South Pacific Central Water (ESPCW; 8.0-24°C; 34.4-36.4), Eastern South Pacific Transition Water (ESPTW; 14.0-20.0°C; 34.6-35.2), and Antarctic Surface Water (AASW; -1.0-1.0°C, 34.0-34.6). Intermediate water masses (500 – 1500 m depth) include Eastern South Pacific Intermediate Water (ESPIW; 10.0–12.0°C, 34.0-34.3) and Antarctic Intermediate Water (2-10°C, 33.8-34.5). Circumpolar Deep Water (CDW; 0.1-2.0°C, 34.62-34.73) was observed below 1500 m depth (Figure 3).

Figure 3

3.2 Distribution of O₂, nutrient concentrations, NO₃^- isotopes, DIN deficit and biogenic N

The P18 transects crossed the ETNP ODZ and the fringe of the ETSP ODZ. O₂ concentrations decreased to < 5 µmol kg^-1 between 120 m and 845 m depth in the northernmost part of the transect in the ETNP ODZ. In the ETSP, O₂ concentrations generally remained > 5 µmol kg^-1, except at few stations/depths between 4.5°S and 9.7°S and 195 and 495 m depth. The oxycline depth varied between 70 m at 10.5°N to 285 m at 16°S (Figure 4A).

Figure 4

Nitrite concentrations were generally low (<2.7 µmol kg^-1), with highest concentrations between 13°N and 16°N and 110 m to 285 m depth in the ETNP ODZ. Nitrite also accumulated at the fringe of the ETSP ODZ, with maximum concentrations between 7.5°S and 16.5°S near 100 m depth. Modest nitrite accumulations (<0.5 µmol kg^-1) were also observed in the upper 200 m from 34°S to 60°S (Figure 4B). Nitrate concentrations were generally depleted in surface waters, with the deepest maximum nitracline found at 350 m depth in the oligotrophic subtropical South Pacific from 21°S to 31.5°S (Figure 4C). A high δ¹⁵N of NO₃^- of up to ~15‰ was observed in surface waters (Figure 4D), as expected following fractionation during NO₃^- assimilation by phytoplankton (Altabet, 2001; Granger et al., 2004). Nitrate concentrations were moderate in the upper 500 m of the ETNP ODZ (minimum [NO₃^-] of 21 µmol kg^-1), corresponding to ENPCW. The δ¹⁵N of NO₃^- was also elevated in the ODZ (up to 16.7‰ at 8.5°N and 425 m depth), consistent with fractionation during denitrification (Cline and Kaplan, 1975; Granger et al., 2008). Nitrate accumulation of up to 45 µmol kg^-1were observed between 625 and 1750 m depth and was more pronounced in the ETNP. Nitrate concentrations remained relatively constant below 2000 m depth, with higher concentrations (~38 µmol kg^-1) in the northern portion of the transect compared to the South (~33 µmol kg^-1) (Figures 4C, D).

The DIN deficit (Howell et al., 1997) was calculated using the following equations:

where DIN_exp is the concentrations of dissolved inorganic nitrogen expected assuming Redfield stoichiometry (typically 16N:1P), DIN_obs is the concentration of nitrate plus nitrite measured and m and b are the slope and intercept, respectively, of the relationship between DIN and phosphate (PO₄^3-) concentrations for source waters outside of the ETNP and ETSP ODZs. DIN_exp was calculated according to Chang et al. (2010; Chang et al., 2012). DIN deficit was particularly elevated (up to 13.75 µmol kg^-1) between 50 and 580 m depth in the ETNP ODZ. DIN deficit was also elevated at the fringes of the ETSP ODZ, with values of up to 12.4 µmol kg^-1 at 9.7°S and 205 m depth. Biogenic N derived from N₂/Ar data (as µmol N kg^-1) were generally comparable to DIN deficit, although generally a bit lower. Discrepancy between these two measurements could be caused by collection (N₂/Ar was not a prioritized gas sample during the P18 cruise) and/or storage effects for the N₂/Ar samples (Figures 4E, F).

3.3 ΔN₂O, N₂O stable isotopes and isotopomers

The largest ΔN₂O accumulation of 75.7 nmol kg^-1 at [O₂] < 5 µmol kg^-1 was observed at a depth of 95 m (σ_θ = 25.65), near the oxycline in the ENPCW water mass in the ETNP ODZ (14.5°N) (Figure 5A). High ΔN₂O of up to 73 nmol kg^-1 were also observed under low-O₂ conditions (15 µmol kg^-1) at the fringe of the ETSP ODZ (7-9°S) at 180-195 m depth (σ_θ = 26.4), which corresponds to the PEW water mass (Figures 3, 5). Overall, elevated ΔN₂O concentrations (>30 nmol kg^-1 were observed from the northern portion of the P18 transect to about 20°S, clearly associated with low-O₂ waters from the ETNP and ETSP ODZs. In contrast, slightly negative ΔN₂O (undersaturation) (~-1 nmol kg^-1) values were observed in surface waters at the southern part of the transect associated with the subduction of AASW forming AAIW (Emery, 2001). Rapid cooling and sinking of surface water, with insufficient time for re-equilibration with the atmosphere, is likely causing the observed N₂O undersaturation in near surface waters in this region. A slight increase in ΔN₂O (up to ~20 nmol kg^-1) was observed from 60°S to 20°S in AAIW waters. Negative ΔN₂O (minimum of ~-5 nmol kg^-1) were observed within the anoxic zone (O₂ < 5 µmol kg^-1) of the ETNP (Figure 5A).

Figure 5

The δ¹⁵N^bulk-N₂O was lowest in the first 500 m depth near the equator, with values less than 5‰, indicating production (Figure 5B). The highest δ¹⁵N^bulk-N₂O values (up to 18‰) were observed within the ETNP ODZ, consistent with the observed low ΔN₂O and indicating net consumption (Ostrom et al., 2007). δ¹⁵N^bulk-N₂O remained approximately 8 to 10‰ in waters deeper than 1000 m depth. The δ¹⁵N^α-N₂O nearly followed the same distribution pattern as δ¹⁵N^bulk-N₂O, with overall higher values (up to 42.5‰) in the ETNP ODZ (Figure 5C). The δ¹⁵N^β-N₂O contrasted with the δ¹⁵N^bulk-N₂O and δ¹⁵N^α-N₂O distributions, with highest values (up to 9.4‰) within the first 500 m flanking the equator and lowest values (minimum of -11.3‰) in the ETNP ODZ (Figure 5D). The δ¹⁸O-N₂O decreased along the north-south gradient of the P18 transect, with highest values (up 100‰) in the ETNP ODZ (Figure 5E). δ¹⁸O-N₂O remained fairly constant below 2000 m depth and increased from ~48 to 56‰ in CDW from south to north along the P18 transect. SP distribution was similar to that of δ¹⁵N^α-N₂O, with lowest values (minimum of -6.1‰) in low O₂ waters (>5 µmol kg^-1) in the upper 500 m near the equator and highest values (up to 49‰) in the ETNP ODZ (Figure 5F). In AASW and AAIW at the southern portion of the transect near Antarctica, δ¹⁵N^bulk, δ¹⁵N^α, δ¹⁵N^β, δ¹⁸O and SP were close to atmospheric N₂O values (6.2‰ for δ¹⁵N^bulk, 15.8‰ for δ¹⁵N^α, -3.4‰ for δ¹⁵N^β, 44.3‰ for δ¹⁸O and 19.2‰ for SP; Kelly et al., 2021). The isotopic signatures of N₂O, together with the observed slight N₂O undersaturation, indicate no N₂O production in this region.

3.4 Relationships between ΔN₂O, AOU and N deficit

The ratio between ΔN₂O/AOU, indicative of cumulative N₂O production yield, increased at low O₂ concentration to 0.3 (Figure 6A). Other high ΔN₂O/AOU ratios observed at higher O₂ concentrations were clearly influenced by the contrasting effect of photosynthesis on AOU, as indicated by higher pH at these shallower isopycnal ranges (Supplementary materials, Figure S1). The lowest SPs were observed at highest ΔN₂O and ΔN₂O/AOU ratios (Figure 6B). Positive relationships between ΔN₂O and AOU were observed in surface waters σ_θ < 26 kg m^-3 in the ETNP and ETSP (Supplementary materials, Figure S2). No significant relationship between ΔN₂O and AOU was observed at 26 < σ_θ < 27 kg m^-3 where highest ΔN₂O values were generally observed.

Figure 6

Significant positive relationships between ΔN₂O and DIN deficit were observed in the ETNP (18.5°N – Equator), with an increased slope toward deeper isopycnals where highest ΔN₂O were observed. A significant positive relationship between ΔN₂O and DIN deficit was also observed in the ETSP (Equator – 13°S) for the isopycnal range 26 < σ_θ < 27 kg m^-3. No such relationship was observed in the ETSP surface waters (σ_θ < 26 kg m^-3) (Figure 7).

Figure 7

4 Discussion

4.1 Sources of the highest ΔN₂O accumulations near the oxycline in the ETNP and ETSP

We observed relatively high N₂O concentrations of up to 84 nmol kg^-1 (940% supersaturation) near the upper oxycline in offshore waters along the P18 section crossing the ETNP ODZ and the fringe of the ETSP ODZ. High N₂O concentrations of up to ~100 nmol kg^-1 were also observed offshore near the oxycline in the ETNP ODZ (Trimmer et al., 2016; Kelly et al., 2021). At offshore stations in the ETSP ODZ, Casciotti et al. (2018) observed N₂O concentrations were up to ~70 nmol kg^-1, comparable to observations along the P18 line. In contrast, Bourbonnais et al. (2017) observed much higher N₂O concentrations of up to ~190 nmol kg^-1 in coastal surface waters off Peru in the ETSP ODZ. Differences in N₂O accumulation observed between these different studies appear to mostly relate to productivity (e.g., coastal versus offshore), but there is also evidence for dynamic N₂O cycling at offshore stations. For instance, the El Niño-Southern Oscillation and mesoscale processes such as eddies have been shown to influence N₂O distribution (e.g., Arévalo-Martínez et al., 2016; Ji et al., 2019; Babbin et al., 2020; Monreal et al., 2022)

Yield of N₂O from nitrification increased toward lower [O₂], as also observed in previous studies (e.g., Ryabenko et al., 2012; Ji et al., 2015; Ji et al., 2018; Frey et al., 2020) (Figure 6). If we assume that 138 moles of O₂ are reduced for 16 moles of NH₄⁺ nitrified (Ward, 2008), a yield of up to 0.5% N-N₂O mol produced per mol NH₄⁺ oxidized is obtained considering the highest ΔN₂O/AOU slope of 0.3. A 0.5% yield falls in the range observed in the Pacific Ocean (Yoshida et al., 1989), but is lower than the maximum N₂O production yield for nitrification reported for ¹⁵N-labeled incubation studies (up to ~3% in the ETSP ODZ) (Ji et al., 2018; Frey et al., 2020). However, while nitrification is clearly occurring in surface waters along the P18 transect, as shown by correlations between ΔN₂O and AOU, the data deviate from this relationship on isopycnals where the highest ΔN₂O were observed near the ODZ oxycline (Supplementary materials, Figure S2).

Two different sources of N₂O were observed from the Keeling plot analysis in both the ETNP and ETSP for [O₂] > 5 µmol kg^-1 and σ_θ < 27 kg m^-3 (Figures 8, 9). A break-point analysis was performed in R using the package “segmented” (Muggeo, 2003; Muggeo and Muggeo, 2017) as described in Kelly et al. (2021). The breakpoint was around 33 nmol kg^-1 (1/[N₂O] = 0.03) in both regions. Coincidentally, this analysis divided the data into two main O₂ regimes: 5 µmol kg^-1 < [O₂] < 100 µmol kg^-1 and [O₂] > 100 µmol kg^-1 (Table 1). It should be noted that highest ΔN₂O in the ETNP were observed at [O₂] below 5 µmol kg^-1 and were thus not included in the Keeling plot analysis. Kelly et al. (2021) found a clear relationship between distinct N₂O sources and [NO₂^-], with highest N₂O accumulation at elevated [NO₂^-] (up to ~1 µmol kg^-1) in the ETNP. In contrast to Kelly et al. (2021), no such relationship was observed when restricting the dataset to the same isopycnal range (σ_θ ≤ 25 kg m^-3).

Figure 8

Kelly et al., 2021). See Table 1 for intercept, R² and p-value for the linear regressions.

Figure 9

Kelly et al., 2021). See Table 1 for intercept, R² and p-value for the linear regressions.

The δ¹⁵N^bulk of source N₂O varied from 4.0 to 7.2‰ (Table 1). δ¹⁵N^bulk is affected by the δ¹⁵N of the substrate and fractionation effects during N₂O production (e.g., Sutka et al., 2006; Frame et al., 2014; Bourbonnais et al., 2017). Similar δ¹⁵N^bulk values for source N₂O were reported in the ETSP (Casciotti et al., 2018). These values fall within the expected range for the δ¹⁵N of remineralized NH₄⁺, the substrate for nitrification, assuming a δ¹⁵N of sinking particulate matter of 3-15‰ (Altabet et al., 1999; Altabet, 2001; Berelson et al., 2015) in the eastern Pacific and a modest isotope effect during ammonification (Altabet, 1988). N₂O could also be produced by denitrification as the δ¹⁵N of NO₃^- observed in this study (4.5-15.4‰) was also similar to observed δ¹⁵N^bulk-N₂O values. However, the upper range in δ¹⁵N values for these substrates is needed to reproduce the observed δ¹⁵N^bulk-N₂O during nitrification and denitrification as both processes are expected to add a lower δ¹⁵N due to isotopic fractionation, unless the substrate is completely consumed, or the reaction is diffusion limited (see Table 1 in Bourbonnais et al., 2017 for a compilation of ε associated with these processes and Ostrom and Ostrom, 2012; Frame et al., 2014).

A high δ¹⁸O (52.3‰) and relatively low δ¹⁵N^α (2.7‰) and SP (-9.0‰) were estimated for source N₂O in the ETNP just above the ODZ corresponding to highest ΔN₂O accumulations. In contrast, δ¹⁵N^bulk, δ¹⁵N^α, δ¹⁵N^β, δ¹⁸O and SP of source N₂O in surface waters ([O₂] higher than 100 µmol kg^-1) were more comparable to atmospheric signatures (Table 1). In the ETSP, the δ¹⁵N^α (11.0‰), δ¹⁸O (55.7‰) and SP (9.5‰) of source N₂O were higher than in the ETNP near the oxycline and the SP in oxic surface waters (11.3‰) was significantly lower. The isotopic values of source N₂O observed in surface waters overlying the ODZ in the ETSP contrasted with the values estimated by Casciotti et al. (2018). For instance, the δ¹⁸O was significantly lower and the SP was higher than the values reported by Casciotti et al. (2018) (Table 1). These differences were likely caused by spatial and/or temporal heterogeneity. Moreover, P18 only crossed the fringe of the ETSP ODZ, with [O₂] generally higher than 5 µmol kg^-1, except at a few stations/depths.

The stable isotopic and isotopomer signatures observed at the highest ΔN₂O and lowest [O₂] near the oxycline in both the ETNP and ETSP are consistent with production from denitrification or nitrifier-denitrification. SP is particularly useful for differentiating N₂O production processes as it is mainly pathway dependent and independent of the isotopic composition of the substrate (Schmidt et al., 2004; Sutka et al., 2004). For instance, formation of N₂O during bacterial denitrification occurs by combining two NO molecules following a trans mechanism (Toyoda et al., 2005; Magyar, 2017). This mechanism involves the asymmetrical intermediate trans-hyponitrite (-ONNO-) bridging two iron centers; the α N atom and the O comes from NO bonded to Fe_heme and the β N atom comes from NO bonded to Fe_B. If the precursor NO molecules are derived from the same substrate, this mechanism causes little difference between the δ¹⁵N of the α and β N atoms, resulting in a low SP (see Magyar, 2017 for more detail). In contrast, during N₂O formation according to a cis mechanism, a first NO molecule binds to one or the other iron center with the second NO molecule binding directly to the first NO, forming a symmetrical intermediate. Cleavage of ¹⁴N-O bond is preferred over ¹⁵N-O bond, leading to enrichment of the α position and a higher SP. Thus, N₂O produced by nitrification, either by archaea and bacteria, is generally associated with a high SP of 30-38‰ consistent with a cis-formation mechanism (Sutka et al., 2003; Sutka et al., 2004; Frame and Casciotti, 2010; Santoro et al., 2011; Löscher et al., 2012) whereas denitrification and nitrifier-nitrification are associated with much lower SPs (~ - 10 to 0‰) (Sutka et al., 2003; Sutka et al., 2004; Toyoda et al., 2005; Frame and Casciotti, 2010).

N₂O produced by denitrification is associated with a low SP of ~ -5 to 0‰ and adds a relatively low δ¹⁵N (especially at the α position) with an isotope effect (¹⁵ε) ranging from 13-37‰ (Barford et al., 1999; Sutka et al., 2003; Sutka et al., 2004; Toyoda et al., 2005; Sutka et al., 2006; Frame and Casciotti, 2010). On the other hand, branching fractionation, i.e., the preferential loss of ¹⁶O relative to ¹⁸O during NO₃^- reduction to N₂O, is expected to lead to high δ¹⁸O-N₂O values (Casciotti et al., 2002; Frame et al., 2014). The δ¹⁸O-N₂O is affected by both the branching isotope effects (¹⁸ϵ) of 25-30‰ for NO₃^- reduction to NO₂^- and 10-12‰ during NO₂^- reduction to N₂O during denitrification as well as the equilibration of NO₂^- O isotope with water (Casciotti et al., 2007; Casciotti and Buchwald, 2012). Nitrifier-denitrification is associated with a SP of -10-0‰ and adds a lower δ¹⁵N compared to denitrification due to a larger ¹⁵ε ranging from 31-58‰. Furthermore, an effective O isotope effect (¹⁸ϵ) of 8-12‰ was reported for NO₂^- reduction to N₂O during nitrifier denitrification (Sutka et al., 2003; Sutka et al., 2004; Sutka et al., 2006; Frame and Casciotti, 2010). Newly produced δ¹⁸O-N₂O values during nitrification and nitrifier denitrification are thus generally lower (13-35‰; Snider et al., 2012). Accordingly, low source δ¹⁸O-N₂O values were observed at [O₂] > 100 µmol kg^-1 in the ETNP and ETSP whereas higher values were observed at deeper isopycnal ranges overlying the ODZs (Table 1), consistent with the dominance of N₂O production by nitrification in surface waters and denitrification at lower O₂ concentrations near the oxycline.

Kelly et al. (2021) similarly observed two distinctive sources of N₂O in the ETNP above the σ_θ < 25 kg m^-3 isopycnal. SPs of 6-8‰ were observed. Based on an isotopic mass balance and assuming that the SPs for N₂O produced during nitrification and denitrification are 30-38‰ and 0‰, respectively (Sutka et al., 2003; Sutka et al., 2004; Toyoda et al., 2005; Frame and Casciotti, 2010; Santoro et al., 2011; Löscher et al., 2012), about 80% of the N₂O production was attributed to denitrification or nitrifier-denitrification. Bourbonnais et al. (2017) and Casciotti et al. (2018) also observed low SPs (-3 to 10‰) corresponding to large ΔN₂O accumulation near the oxycline above the ETSP ODZ. The much lower SP observed in the ETNP in this study (-9‰) preclude any contribution from nitrification. The higher SP (9.5‰) for source N₂O observed at low [O₂] near the fringe of the ETSP ODZ suggests a relatively minor contribution from nitrification (~1/3) and is consistent with the value reported by Casciotti et al. (2018). However, hybrid N₂O formation by archaea, where one atom is derived from nitric oxide (NO; from NO₂^-) and the other from hydroxylamine (NH₂OH; from NH₄⁺) was shown to be an important production pathway in marine environments and could possibly contribute to endmember signatures estimated from the Keeling plot analysis (Stieglmeier et al., 2014; Kozlowski et al., 2016; Trimmer et al., 2016; Frame et al., 2017). Frey et al. (2020) showed that hybrid N₂O production accounts for 70-85% of the total N₂O production from NH₄⁺ oxidation using ¹⁵N-labeled incubation experiments in the ETSP. SP may, in part, reflect the relative δ¹⁵N of the substrates for the α and β positions during hybrid archaeal N₂O production rather than being indicative of a particular pathway (Casciotti et al., 2018). Therefore, if lower SP could be produced by archaeal hybrid N₂O production, the fraction of denitrification needed to explain low SP signatures associated with high N₂O supersaturations observed above or in the ODZs in this and previous studies would further decrease (Bourbonnais et al., 2017; Casciotti et al., 2018; Kelly et al., 2021).

The observed relationships between ΔN₂O and DIN deficit at isopycnals σ_θ < 26 kg m^-3 and the absence of clear relationships between ΔN₂O and [NO₂^-] further support a role for denitrification (rather than nitrifier-denitrification) as the dominant N₂O production pathway in both the ETSP and the ETNP. The decreased slope for the ΔN₂O versus DIN deficit relationship toward surface isopycnal ranges observed in the ETNP suggests decreased N₂O yield at higher [O₂] concentrations or mixing/dilution of N₂O produced by denitrification advected from below, especially for upwelling waters near the equator. The δ¹⁵N of NO₃^- provides more insights into N₂O sources yet these measurements were mostly available at ETSP stations. A plot of SP versus Δ¹⁵N (δ¹⁵N-NO₃^- - δ¹⁵N-N₂O) showed that some of the highest ΔN₂O datapoint indeed fell within the expected compositional fields for N₂O production by bacterial denitrification (Δδ¹⁵N= 0-35‰, SP: -5-0‰) (Wankel et al. (2017) and references therein, Supplementary Figure S3). Some high ΔN₂O were associated with relatively high SPs (> 10‰), showing the overprinting effect of N₂O consumption during denitrification (Ostrom et al., 2007). Compositional field analysis (Figure S3) also suggests production by either archaeal and bacterial ammonia oxidation at higher [O₂] and lower ΔN₂O. Bourbonnais et al. (2017) also invoked incomplete denitrification as a major pathway for extreme N₂O accumulation in newly upwelled surface waters off Peru based on the absence of a relationship between ΔN₂O and apparent O₂ utilization and significant relationships between NO₃^- and N₂O isotopes. The dominance of denitrification for N₂O production near the oxycline is also supported by ¹⁵N-tracer incubation studies in the ETSP ODZ. N₂O production rates were indeed one order of magnitude higher for denitrification compared to ammonia oxidation under low-O₂ conditions at the ΔN₂O maximum just above the ETSP ODZ (Ji et al., 2015; Frey et al., 2020).

4.2 Sources of N₂O below the ETNP and ETSP ODZs and in AAIW

N₂O sources in deep ETNP and ETSP were also investigated using Keeling plot analysis (Supplementary materials, Figure S4 and Table 1). A deeper isopycnal range (σ_θ > 27.3 kg m^-3was selected for the ETNP due to the deeper ODZ at this location as in a previous study (Kelly et al., 2021). An isopycnal range with σ_θ > 27 kg m^-3 was selected for the ETSP comparable to the σ_θ range in Casciotti et al. (2018). N₂O source values of 9.1‰ (δ¹⁵N^bulk), 16.8‰ (δ¹⁵N^α), 1.4‰ (δ¹⁵N^β), 52.7‰ (δ¹⁸O) and 15.4‰ (SP) were estimated below the ETNP ODZ (Table 1). In the ETSP, the Keeling plot analysis revealed values of 7.0‰ (δ¹⁵N^bulk), 17.0‰ (δ¹⁵N^α), -3.0‰ (δ¹⁵N^β), 45.5‰ (δ¹⁸O) and 20.0‰ (SP). The generally higher δ¹⁵N^β as well as lower δ¹⁸O and SP compared to prior studies (Casciotti et al., 2018; Kelly et al., 2021) might be the result of a lesser influence from N₂O production and consumption in the ODZ since our analysis also included equatorial waters. Overall, these isotopic signatures suggest the dominance of N₂O production by nitrification in deeper waters below the ETNP and ESTP ODZs. In fact, our results are more in line with a N₂O source from nitrification as a relatively low δ¹⁸O is expected for this process (i.e., 13-35‰; Snider et al., 2012).

The sources of N₂O were investigated by restricting the Keeling plot analysis for the absolute salinity/potential temperature ranges characteristic of AAIW as described in section 3.1. The analysis was restricted to 60°S to 20°S, even if the northernmost extent of AAIW is found below the ETNP, to eliminate possible effects of the ETNP and ETSP ODZs on N₂O production in this water mass.

ΔN₂O clearly increased from -0.89 nmol kg^-1 (-7.4% supersaturation) at 52°S to up to 20.6 nmol kg^-1 (180% supersaturation) in the intermediate water mass AAIW (Figure 5), which is comparable to values observed by Carrasco et al. (2017). Keeling plot analysis revealed a N₂O source with a δ¹⁵N^bulk = 9.4‰, δ¹⁵N^α = 18.1‰, δ¹⁵N^β = 0.8‰, δ¹⁸O = 52.5‰, and SP = 17.2‰ (Table 1). The relatively high SP suggests that N₂O is mainly derived from nitrification, consistent with previous studies (Casciotti et al., 2018; Toyoda et al., 2019; Kelly et al., 2021). The positive significant relationship observed between ΔN₂O and AOU in the AAIW water mass observed in this study and Carrasco et al. (2017) further supports this interpretation (Supplementary materials, Figure S2). The lower δ¹⁵N^β and higher δ¹⁸O and SP might reflect the influence of N₂O consumption in the ODZs in these prior studies.

4.3 Isotopic signatures of N₂O consumption in the ETNP

High values of δ¹⁵N^α, δ¹⁸O and SP and low or negative ΔN₂O were observed at [O₂] < 5 µmol kg^-1 in the ETNP ODZ, which are clear signatures of N₂O consumption (Figure 5). During N₂O reduction to N₂ gas, the N-O bond is broken, leaving the remaining substrate (N₂O) enriched in ¹⁵N and ¹⁸O. The α position in N₂O is preferentially enriched in ¹⁵N compared to the β position since it is directly attached to the O atom being cleaved (e.g., Popp et al., 2002; Toyoda et al., 2002). Thus, ε for the α N atom (6.6-9.1‰) and δ¹⁸O (10.9-15‰) are relatively large with only a small or negligible ε for β N atom (Ostrom et al., 2007). Notably, while δ¹⁵N^α and δ¹⁸O increased in the ODZ, δ¹⁵N^β generally decreased, as also observed in previous ETNP and ETSP ODZ studies (Bourbonnais et al., 2017; Casciotti et al., 2018; Kelly et al., 2021) (Figure 5). This trend is not expected during pure N₂O consumption.

The slope for the relationship between δ¹⁸O versus and δ¹⁵N^α (1.6 ± 0.1) was indistinguishable from the expected slope (i.e., 1.7) during N₂O consumption in soils and pure denitrifier cultures (Ostrom et al., 2007). However, the observed slope for δ¹⁸O-N₂O versus SP (0.8 ± 0.1) clearly deviated from the expected value of 2.2 for pure N₂O reduction (Ostrom et al., 2007) (Figure 10). A wider range of 1.3-3.3 for δ¹⁸O versus SP during consumption by denitrification was also reported for soil mesocosms, which is still higher than the slope observed in this study (Lewicka-Szczebak et al., 2017). This observation is consistent with previous studies in marine ODZs, reporting slopes for δ¹⁸O versus SP ranging between (0.9-1.8) (Bourbonnais et al., 2017; Casciotti et al., 2018; Kelly et al., 2021). For instance, in the ODZ off Peru, the slope for increase in δ¹⁸O-N₂O versus SP deviates from what is expected during pure denitrification, mostly due to a decrease in δ¹⁵N^β within the ODZ (Bourbonnais et al., 2017; Casciotti et al., 2018). This trend showing a decreasing δ¹⁵N^β associated with high δ¹⁸O-N₂O has also been observed in other ODZs and marine anoxic environments (Yamagishi et al., 2005; Westley et al., 2006; Yamagishi et al., 2007; Farías et al., 2009; Kelly et al., 2021). The cause for the decreasing δ¹⁵N^β in the ODZ will be further discussed below.

Figure 10

Ostrom et al., 2007). ΔN₂O concentrations are shown as the color of the symbols (see color bar in top panel). Linear regressions for O₂ < 5 µmol kg^-1 are shown with black lines.

In this study, apparent isotope effects (ε_app) were derived for all isopycnal ranges within the ETNP ODZ at [O₂] < 5 µmol kg^-1, where N₂O consumption is known to occur (Dalsgaard et al., 2014). This approach was used in Casciotti et al. (2018) to further investigate if a decrease in ε¹⁸O or increase in εSP could explain the deviation from the expected slope for δ¹⁸O-N₂O versus SP during pure N₂O consumption in the ETSP ODZ. Apparent isotope effects calculated for all isopycnal ranges were 3.6‰ for ¹⁵N^bulk, 9.4‰ for ¹⁵N^α, -2.3‰ for ¹⁵N^β, 12.0‰ for ¹⁸O and 11.7‰ for SP (Table 2 and Figure 11). The ε_app for ¹⁵N^bulk was slightly lower than the range reported for pure culture by Ostrom et al. (2007). Yet, the ε_app for ¹⁵N^α and ¹⁵N^β were comparable to values (11.8‰ and -2‰, respectively) estimated by Kelly et al. (2021) in the ETNP ODZ. Notably, the negative isotope effect for ¹⁵N^β indicates a decrease in δ¹⁵N^β as N₂O is reduced to N₂ within the ODZ. The ε for ¹⁸O calculated in this study was significantly lower than the values (14.5 – 20.2‰) estimated in other ODZ marine field studies (Casciotti et al., 2018; Kelly et al., 2021), but comprised within the range observed for pure culture (10.9 – 15‰) (Ostrom et al., 2007). The ε for SP was higher than for pure culture (5.0 – 6.8‰) but comparable to values estimated in the ETSP ODZ (11.6 – 17.5‰) (Casciotti et al., 2018). Casciotti et al. (2018) observed different apparent isotope effects at different isopycnal ranges within the ETSP ODZ. They observed a decrease in the ε¹⁸O within the ODZ (from 25.6‰ in the upper, to 19.6‰ in the middle and 14.5‰ in the lower ODZ) while the εSP increased slightly and then decreased (17.5‰ in the upper to 21.7‰ in the middle to 11.6‰ in the lower ODZ). A more detailed analysis by isopycnal ranges (upper, middle, and lower ODZ) could however not reproduce observed trends in Casciotti et al. (2018) for the P18 transect. For instance, ε_app for ¹⁸O (13.9‰ in the upper ODZ to 11.4‰ in the middle, and 16.2‰ in the lower ODZ) and SP (14.2‰ in the upper ODZ to 7.5‰ in the middle, and 17.9‰ in the lower ODZ) both decreased toward the middle ODZ and then increased in the deeper ODZ (Supplementary materials, Figure S5). The slope for ε¹⁸O versus εSP also gradually decreased from 1.6 to 0.65 instead of being lowest in the middle ODZ. Irrespective of the observed trends, using a modeling approach and isotope values of substrates (NO₃^- and N¹⁸O₂^-) during denitrification, previous studies attributed the lower ε¹⁸O:εSP to an increase in SP due to a decrease in δ¹⁵N^β rather than a decrease in ε¹⁸O or increase in εSP in both the ETNP and ESTP ODZs (Casciotti et al., 2018; Kelly et al., 2021).

Figure 11

Several hypotheses have been proposed to explain the decreasing δ¹⁵N^β in ODZ waters, including concurrent N₂O production via denitrification of NO₃^- with a site preference >0‰ (Schmidt et al., 2004; Casciotti et al., 2018; Kelly et al., 2021). Rate experiments as well as prior stable abundance isotopic studies suggest concurrent N₂O production and consumption as well as a rapid N₂O turnover in the ODZ (Farías et al., 2009; Babbin et al., 2015; Ji et al., 2015; Bourbonnais et al., 2017; Casciotti et al., 2018; Frey et al., 2020; Kelly et al., 2021). Casciotti et al. (2018) modeled N₂O cycling at steady-state in the ODZ using both δ¹⁵N and δ¹⁸O of substrate molecules (NO₂^- and NO₃^-) during denitrification. N₂O isotopic signatures were only reproduced when NO₃^- was used as the substrate and implying a site preference >0‰ (i.e., greater fractionation at the ¹⁵N^β position relative to the ¹⁵N^α) during N₂O production by denitrification. This suggests that an internal (rather than ambient) pool of NO₂^- is used during denitrification as further discussed in Casciotti et al. (2018). Accordingly, a relatively long residence time was estimated for ambient NO₂^- (in the order of months) in the primary NO₂^- maximum and low-O₂ waters based on rates of abiotic O isotope exchange between nitrite and water (Buchwald and Casciotti, 2013; Bourbonnais et al., 2015). In another modeling study, a non-steady-state N₂O cycling as well as an ambient NO₂^- substrate source (with a δ¹⁵N as low as -30‰) or a SP of ~25‰ during concurrent N₂O production by denitrification were required to explain the observed low δ¹⁵N^β values in the ETNP ODZ (Kelly et al., 2021). The low δ¹⁵N^β of N₂O could also result from an alternate mechanism in the ODZ, such as N₂O production from AO archaea and bacteria, which is associated with a high SP (30-38‰; Sutka et al., 2003; Sutka et al., 2004; Frame and Casciotti, 2010; Santoro et al., 2011; Löscher et al., 2012). For instance, an archaeal AO isolated from ODZs (Nitrosopumilus maritimus) produced both N₂ and O₂ under dark anaerobic conditions, following a pathway that involves N₂O as an intermediate (Kraft et al., 2022). Yet, N₂O production rates from AO based on ¹⁵N-labeled experiments were relatively low (up to 0.1 nmol L^-1 d^-1) and generally at least one order of magnitude lower than N₂O production rates from NO₃^- (denitrification) under anoxic conditions in the ETSP ODZ (Frey et al., 2020).

The idea of denitrification with a site preference >0‰ is not new (i.e., see Schmidt et al., 2004). Several studies invoked N₂O production during denitrification with a site preference >0‰ in the ETNP and ETSP ODZs and Canadian Arctic bottom waters influenced by sedimentary processes (Casciotti et al., 2018; Lehmann et al., 2019; Kelly et al., 2021). Some variability in the SP of N₂O produced by denitrification was reported for different bacterial strains. For instance, Toyoda et al. (2005) observed a high SP of 22 to 24‰ for N₂O production during denitrification by Pseudomonas fluorescens, suggesting a symmetrical intermediate (i.e., cis formation mechanism). These observations for denitrifying bacteria are analogous to the wide range of SPs (15.8 – 37.1‰) that have been reported for different fungal species and strains (Maeda et al., 2015; Lazo-Murphy et al., 2022). These data challenge the conventional view that fungal denitrification fits into narrow compositional fields with relatively high SP values. Similarly, it is possible that denitrifying microbial consortia thriving under different O₂ regimes are associated with variable SPs, but this hypothesis remains to be verified.

5 Concluding remarks

This study presents oceanic N₂O concentration, stable isotope and isotopomer data of unprecedently high spatial resolution along the P18 line in the eastern Pacific sampled in 2016/2017. Highest ΔN₂O accumulations (up to 940% supersaturation) were observed close to the oxycline in both the ETNP and the fringe of the ETSP ODZ. N₂O yield from AOU increased at lower [O₂], consistent with previous studies. Keeling plot analysis identified two distinct sources of N₂O at different [O₂] regimes. At lower [O₂] concentrations close to the oxycline in both ODZs, where highest ΔN₂O were observed, SP of source N₂O was relatively low, suggesting production from denitrification (or nitrifier-denitrification). Relationships between ΔN₂O and DIN deficit and the isotopic composition of δ¹⁵N of the substrate (NO₃^-) further suggest that denitrification is the dominant process at low [O₂] concentrations. SP generally increased in more oxygenated surface waters, suggesting a greater contribution from nitrification. The isotopic composition of source N₂O in deeper waters of the ETNP and ETSP as well as in AAIW also suggested that nitrification was the main pathway for N₂O formation.

The δ¹⁸O versus δ¹⁵N^α relationship showed a slope characteristic of N₂O consumption during denitrification in low O₂ waters. Isotope effects calculated for N₂O consumption were consistent with previous field and laboratory studies (Ostrom et al., 2007; Casciotti et al., 2018; Kelly et al., 2021). Yet, δ¹⁵N^β decreased (rather than the predicted no change or slight increase during denitrification) in the ODZ, which was also observed by previous studies (e.g., Bourbonnais et al., 2017; Casciotti et al., 2018; and Kelly et al., 2021). These signatures are best explained by concurrent N₂O production (from NO₃^- or NO₂^-) with a site preference >0‰, an ambient NO₂^- source with a low δ¹⁵N and non-steady-state conditions (e.g., Kelly et al., 2021). Clearly, more research is needed to elucidate the SP and its variability for diverse denitrifying bacterial strains living in marine ODZs as current values are mostly from terrestrial environments (e.g., Toyoda et al., 2005). Some contribution from an alternative N₂O production pathway is also possible, for example anaerobic AO (see Kraft et al., 2022).

This study establishes a benchmark against which to evaluate changes in N₂O cycling for future decadal occupations of the P18 line. ODZ are currently expanding (Stramma et al., 2013), with unknown impacts on N₂O cycling. The vertical expansion of ODZs not only increases the volume of low-O₂ waters where N₂O is potentially produced but also increases N₂O exchange with the atmosphere. However, co-occurring warming causes stronger stratification, which could reduce mixing and N₂O outgassing to the atmosphere. At this point, it is unclear which mechanism will dominate. Thus, more observational data is needed to evaluate the impacts of ODZs expansion on marine N₂O cycling and atmospheric emissions.

References

• 1

  AltabetM. A. (1988). Variations in nitrogen isotopic composition between sinking and suspended particles: Implications for nitrogen cycling and particle transformation in the open ocean. Deep Sea Res. Part A. Oceanographic Res. Papers35 (4), 535–554. doi: 10.1016/0198-0149(88)90130-6

  • CrossRef
  • Google Scholar
  • View reference in article

• 2

  AltabetM. A. (2001). Nitrogen isotopic evidence for micronutrient control of fractional NO₃^– utilization in the equatorial pacific. Limnol. Oceanogr.46 (2), 368–380. doi: 10.4319/lo.2001.46.2.0368

  • CrossRef
  • Google Scholar
  • View reference in article

• 3

  AltabetM. A. (2006). “Isotopic tracers of the marine nitrogen cycle: present and past,” in Marine organic matter: Chemical and biological markers. the handbook of environmental chemistry. Ed. VolkmanJ. (Berlin: Springer-Verlag), pp251–pp293.

  • Google Scholar
  • View reference in article

• 4

  AltabetM. A.PilskalnC.ThunellR.PrideC.SigmanD.ChavezF.et al. (1999). The nitrogen isotope biogeochemistry of sinking particles from the margin of the Eastern north pacific. Deep Sea Res. Part I: Oceanographic Res. Papers46 (4), 655–679. doi: 10.1016/S0967-0637(98)00084-3

  • CrossRef
  • Google Scholar
  • View reference in article

• 5

  AltabetM. A.WassenaarL. I.DouenceC.RoyR. (2019). A Ti (III) reduction method for one-step conversion of seawater and freshwater nitrate into N2O for stable isotopic analysis of ¹⁵N/¹⁴N, ¹⁸O/¹⁶O and ¹⁷O/¹⁶O. Rapid Commun. Mass Spectrometry33 (15), 1227–1239. doi: 10.1002/rcm.8454

  • CrossRef
  • Google Scholar
  • View reference in article

• 6

  Arévalo-MartínezD. L.KockA.LöscherC. R.SchmitzR. A.BangeH. W. (2015). Massive nitrous oxide emissions from the tropical south pacific ocean. Nat. Geosci.8 (7), 530–533. doi: 10.1038/ngeo2469

  • CrossRef
  • Google Scholar
  • View reference in article

• 7

  Arévalo-MartínezD. L.KockA.LöscherC. R.SchmitzR. A.StrammaL.BangeH. W. (2016). Influence of mesoscale eddies on the distribution of nitrous oxide in the eastern tropical south pacific. Biogeosciences13 (4), 1105–1118. doi: 10.5194/bg-13-1105-2016

  • CrossRef
  • Google Scholar
  • View reference in article

• 8

  BabbinA. R.BianchiD.JayakumarA.WardB. B. (2015). Rapid nitrous oxide cycling in the suboxic ocean. Science348 (6239), 1127–1129. doi: 10.1126/science.aaa8380

  • CrossRef
  • Google Scholar
  • View reference in article

• 9

  BabbinA. R.BolesE. L.MühleJ.WeissR. F. (2020). On the natural spatio-temporal heterogeneity of south pacific nitrous oxide. Nat. Commun.11 (1), 1–9. doi: 10.1038/s41467-020-17509-6

  • CrossRef
  • Google Scholar
  • View reference in article

• 10

  BakkerD. C.BangeH. W.GruberN.JohannessenT.Upstill-GoddardR. C.BorgesA. V.et al. (2014). Air-sea interactions of natural long-lived greenhouse gases (CO₂, N₂O, CH₄) in a changing climate. Ocean-Atmosphere Interactions of Gases and Particles. (Switzerland AG: Springer Nature), pp113–169.

  • Google Scholar
  • View reference in article

• 11

  BangeH. W.AndreaeM. O.LalS.LawC. S.NaqviS. W. A.PatraP. K.et al. (2001). Nitrous oxide emissions from the Arabian Sea: A synthesis. Atmospheric Chem. Phys.1 (1), 61–71. doi: 10.5194/acp-1-61-2001

  • CrossRef
  • Google Scholar
  • View reference in article

• 12

  BarfordC. C.MontoyaJ. P.AltabetM. A.MitchellR. (1999). Steady-state nitrogen isotope effects of N₂ and N₂O production in paracoccus denitrificans. Appl. Environ. Microbiol.65 (3), 989–994. doi: 10.1128/AEM.65.3.989-994.1999

  • CrossRef
  • Google Scholar
  • View reference in article

• 13

  BerelsonW. M.HaskellW. Z.IIProkopenkoM.KnappA. N.HammondD. E.RollinsN.et al. (2015). Biogenic particle flux and benthic remineralization in the Eastern tropical south pacific. Deep Sea Res. Part I: Oceanographic Res. Papers99, 23–34. doi: 10.1016/j.dsr.2014.12.006

  • CrossRef
  • Google Scholar
  • View reference in article

• 14

  BourbonnaisA.AltabetM. A.CharoenpongC. N.LarkumJ.HuH.BangeH. W.et al. (2015). N-loss isotope effects in the Peru oxygen minimum zone studied using a mesoscale eddy as a natural tracer experiment. Global Biogeochemical Cycles29 (6), 793–811. doi: 10.1002/2014GB005001

  • CrossRef
  • Google Scholar
  • View reference in article

• 15

  BourbonnaisA.LehmannM. F.WaniekJ. J.Schulz-BullD. E. (2009). Nitrate isotope anomalies reflect N₂ fixation in the Azores front region (subtropical NE Atlantic). J. Geophys. Res.: Oceans114 (C3). doi: 10.1029/2007JC004617

  • CrossRef
  • Google Scholar
  • View reference in article

• 16

  BourbonnaisA.LetscherR. T.BangeH. W.EchevinV.LarkumJ.MohnJ.et al. (2017). N₂O production and consumption from stable isotopic and concentration data in the Peruvian coastal upwelling system. Global Biogeochemical Cycles31 (4), 678–698. doi: 10.1002/2016GB005567

  • CrossRef
  • Google Scholar
  • View reference in article

• 17

  BourbonnaisA.FreyC.SunX.BristowL. A.JayakumarA.OstromN. E.et al. (2021). Protocols for assessing transformation rates of nitrous oxide in the water column. Front. Mar. Sci.8, 611937. doi: 10.3389/fmars.2021.611937

  • CrossRef
  • Google Scholar
  • View reference in article

• 18

  BuchwaldC.CasciottiK. L. (2010). Oxygen isotopic fractionation and exchange during bacterial nitrite oxidation. Limnol. Oceanogr.55 (3), 1064–1074. doi: 10.4319/lo.2010.55.3.1064

  • CrossRef
  • Google Scholar
  • View reference in article

• 19

  BuchwaldC.CasciottiK. L. (2013). Isotopic ratios of nitrite as tracers of the sources and age of oceanic nitrite. Nat. Geosci.6(4), 308–313.

  • Google Scholar
  • View reference in article

• 20

  CarantoJ. D.LancasterK. M. (2017). Nitric oxide is an obligate bacterial nitrification intermediate produced by hydroxylamine oxidoreductase. Proc. Natl. Acad. Sci.114 (31), 8217–8222. doi: 10.1073/pnas.1704504114

  • CrossRef
  • Google Scholar
  • View reference in article

• 21

  CarrascoC.KarstensenJ.FariasL. (2017). On the nitrous oxide accumulation in intermediate waters of the eastern south pacific ocean. Front. Mar. Sci.4, 24. doi: 10.3389/fmars.2017.00024

  • CrossRef
  • Google Scholar
  • View reference in article

• 22

  CasciottiK. L.BöhlkeJ. K.McIlvinM. R.MroczkowskiS. J.HannonJ. E. (2007). Oxygen isotopes in nitrite: Analysis, calibration, and equilibration. Analytical Chem.79 (6), 2427–2436. doi: 10.1021/ac061598h

  • CrossRef
  • Google Scholar
  • View reference in article

• 23

  CasciottiK. L.BuchwaldC. (2012). Insights on the marine microbial nitrogen cycle from isotopic approaches to nitrification. Front. Microbiol.3, 356. doi: 10.3389/fmicb.2012.00356

  • CrossRef
  • Google Scholar
  • View reference in article

• 24

  CasciottiK. L.ForbesM.VedamatiJ.PetersB. D.MartinT. S.MordyC. W. (2018). Nitrous oxide cycling in the Eastern tropical south pacific as inferred from isotopic and isotopomeric data. Deep Sea Res. Part II: Topical Stud. Oceanogr.156, 155–167. doi: 10.1016/j.dsr2.2018.07.014

  • CrossRef
  • Google Scholar
  • View reference in article

• 25

  CasciottiK. L.SigmanD. M.HastingsM. G.BöhlkeJ. K.HilkertA. (2002). Measurement of the oxygen isotopic composition of nitrate in seawater and freshwater using the denitrifier method. Analytical Chem.74 (19), 4905–4912. doi: 10.1021/ac020113w

  • CrossRef
  • Google Scholar
  • View reference in article

• 26

  ChangB. X.DevolA. H.EmersonS. R. (2010). Denitrification and the nitrogen gas excess in the eastern tropical south pacific oxygen deficient zone. Deep Sea Res. Part I: Oceanographic Res. Papers57 (9), 1092–1101. doi: 10.1016/j.dsr.2010.05.009

  • CrossRef
  • Google Scholar
  • View reference in article

• 27

  ChangB. X.DevolA. H.EmersonS. R. (2012). Fixed nitrogen loss from the eastern tropical north pacific and Arabian Sea oxygen deficient zones determined from measurements of N₂: Ar. Global Biogeochem. Cycles26 (3). doi: 10.1029/2011GB004207

  • CrossRef
  • Google Scholar
  • View reference in article

• 28

  CharoenpongC. N.BristowL. A.AltabetM. A. (2014). A continuous flow isotope ratio mass spectrometry method for high precision determination of dissolved gas ratios and isotopic composition. Limnol. Oceanogr.: Methods12 (5), 323–337. doi: 10.4319/lom.2014.12.323

  • CrossRef
  • Google Scholar
  • View reference in article

• 29

  CiaisP.SabineC.BalaG.BoppL.BrovkinV.CanadellJ.et al. (2013). Carbon and Other Biogeochemical Cycles. In: Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change [StockerT. F.QinD.PlattnerG.-K.TignorM.AllenS. K.BoschungJ.et al (eds.)]. Cambridge University Press, Cambridge, United Kingdom and New York, NY, USA, pp465, pp570.

  • Google Scholar
  • View reference in article

• 30

  ClineJ. D.KaplanI. R. (1975). Isotopic fractionation of dissolved nitrate during denitrification in the eastern tropical north pacific ocean. Mar. Chem.3 (4), 271–299. doi: 10.1016/0304-4203(75)90009-2

  • CrossRef
  • Google Scholar
  • View reference in article

• 31

  CohenY.GordonL. I. (1979). Nitrous oxide production in the ocean. J. Geophys. Res.: Oceans84 (C1), 347–353. doi: 10.1029/JC084iC01p00347

  • CrossRef
  • Google Scholar
  • View reference in article

• 32

  DalsgaardT.StewartF. J.ThamdrupB.De BrabandereL.RevsbechN. P.UlloaO.et al. (2014). Oxygen at nanomolar levels reversibly suppresses process rates and gene expression in anammox and denitrification in the oxygen minimum zone off northern Chile. MBio5 (6), e01966–e01914. doi: 10.1128/mBio.01966-14

  • CrossRef
  • Google Scholar
  • View reference in article

• 33

  EmeryW. J. (2001). Water types and water masses. Encyclopedia ocean Sci.6, 3179–3187. doi: 10.1006/rwos.2001.0108

  • CrossRef
  • Google Scholar
  • View reference in article

• 34

  FaríasL.Castro-GonzálezM.CornejoM.CharpentierJ.FaúndezJ.BoontanonN.et al. (2009). Denitrification and nitrous oxide cycling within the upper oxycline of the eastern tropical south pacific oxygen minimum zone. Limnol. Oceanogr.54 (1), 132–144. doi: 10.4319/lo.2009.54.1.0132

  • CrossRef
  • Google Scholar
  • View reference in article

• 35

  ForsterP.StorelvmoT.ArmourK.CollinsW.DufresneJ.-L.FrameD.et al. (2021). “The earth’s energy budget, climate feedbacks, and climate sensitivity,” in Climate change 2021: The physical science basis. contribution of working group I to the sixth assessment report of the intergovernmental panel on climate change. Eds. [Masson-DelmotteV.ZhaiP.PiraniA.ConnorsS. L.PéanC.BergerS.CaudN.ChenY.GoldfarbL.GomisM. I.HuangM.LeitzellK.LonnoyE.MatthewsJ. B. R.MaycockT. K.WaterfieldT.YelekçiO.YuR.ZhouB. (Cambridge, United Kingdom and New York, NY, USA: Cambridge University Press), pp923–p1054. doi: 10.1017/9781009157896.009

  • CrossRef
  • Google Scholar
  • View reference in article

• 36

  FrameC. H.CasciottiK. L. (2010). Biogeochemical controls and isotopic signatures of nitrous oxide production by a marine ammonia-oxidizing bacterium. Biogeosciences7 (9), 2695–2709. doi: 10.5194/bg-7-2695-2010

  • CrossRef
  • Google Scholar
  • View reference in article

• 37

  FrameC. H.DealE.NevisonC. D.CasciottiK. L. (2014). N₂O production in the eastern south Atlantic: Analysis of N₂O stable isotopic and concentration data. Global Biogeochemical Cycles28 (11), 1262–1278. doi: 10.1002/2013GB004790

  • CrossRef
  • Google Scholar
  • View reference in article

• 38

  FrameC. H.LauE.NolanE. J.IVGoepfertT. J.LehmannM. F. (2017). Acidification enhances hybrid N₂O production associated with aquatic ammonia-oxidizing microorganisms. Front. Microbiol.7, 2104. doi: 10.3389/fmicb.2016.02104

  • CrossRef
  • Google Scholar
  • View reference in article

• 39

  FreingA.WallaceD. W.BangeH. W. (2012). Global oceanic production of nitrous oxide. Philos. Trans. R. Soc. B: Biol. Sci.367 (1593), 1245–1255. doi: 10.1098/rstb.2011.0360

  • CrossRef
  • Google Scholar
  • View reference in article

• 40

  FreingA.WallaceD. W. R.TanhuaT.WalterS.BangeH. W. (2009). North Atlantic production of nitrous oxide in the context of changing atmospheric levels. Global Biogeochemical Cycles23, 13. doi: 10.1029/2009gb003472

  • CrossRef
  • Google Scholar
  • View reference in article

• 41

  FreyC.BangeH. W.AchterbergE. P.JayakumarA.LöscherC. R.Arévalo-MartínezD. L.et al. (2020). Regulation of nitrous oxide production in low-oxygen waters off the coast of Peru. Biogeosciences17 (8), 2263–2287. doi: 10.5194/bg-17-2263-2020

  • CrossRef
  • Google Scholar
  • View reference in article

• 42

  FujiiA.ToyodaS.YoshidaO.WatanabeS.SasakiK. I.YoshidaN. (2013). Distribution of nitrous oxide dissolved in water masses in the eastern subtropical north pacific and its origin inferred from isotopomer analysis. J. Oceanogr.69 (2), 147–157. doi: 10.1007/s10872-012-0162-4

  • CrossRef
  • Google Scholar
  • View reference in article

• 43

  GoreauT. J.KaplanW. A.WofsyS. C.McElroyM. B.ValoisF. W.WatsonS. W. (1980). Production of NO₂^- and N₂O by nitrifying bacteria at reduced concentrations of oxygen. Appl. Environ. Microbiol.40 (3), 526–532. doi: 10.1128/aem.40.3.526-532.1980

  • CrossRef
  • Google Scholar
  • View reference in article

• 44

  GrangerJ.SigmanD. M.LehmannM. F.TortellP. D. (2008). Nitrogen and oxygen isotope fractionation during dissimilatory nitrate reduction by denitrifying bacteria. Limnol. Oceanogr.53 (6), 2533–2545. doi: 10.4319/lo.2008.53.6.2533

  • CrossRef
  • Google Scholar
  • View reference in article

• 45

  GrangerJ.SigmanD. M.NeedobaJ. A.HarrisonP. J. (2004). Coupled nitrogen and oxygen isotope fractionation of nitrate during assimilation by cultures of marine phytoplankton. Limnol. Oceanogr.49 (5), 1763–1773. doi: 10.4319/lo.2004.49.5.1763

  • CrossRef
  • Google Scholar
  • View reference in article

• 46

  HeilJ.WolfB.BrüggemannN.EmmeneggerL.TuzsonB.VereeckenH.et al. (2014). Site-specific ¹⁵N isotopic signatures of abiotically produced N₂O. Geochimica Cosmochimica Acta139, 72–82. doi: 10.1016/j.gca.2014.04.037

  • CrossRef
  • Google Scholar
  • View reference in article

• 47

  HinkL.LycusP.Gubry-RanginC.FrostegårdÅ.NicolG. W.ProsserJ. I.et al. (2017a). Kinetics of NH₃-oxidation, NO-turnover, N₂O-production and electron flow during oxygen depletion in model bacterial and archaeal ammonia oxidisers. Environ. Microbiol.19 (12), 4882–4896. doi: 10.1111/1462-2920.13914

  • CrossRef
  • Google Scholar
  • View reference in article

• 48

  HinkL.NicolG. W.ProsserJ. I. (2017b). Archaea produce lower yields of N₂O than bacteria during aerobic ammonia oxidation in soil. Environ. Microbiol.19 (12), 4829–4837.

  • Google Scholar
  • View reference in article

• 49

  HowellE. A.DoneyS. C.FineR. A.OlsonD. B. (1997). Geochemical estimates of denitrification rates for the Arabian Sea and bay of Bengal during WOCE. Geophys. Res. Lett.24, 2549–2552. doi: 10.1029/97GL01538

  • CrossRef
  • Google Scholar
  • View reference in article

• 50

  JiQ.AltabetM. A.BangeH. W.GracoM. I.MaX.Arévalo-MartínezD. L.et al. (2019). Investigating the effect of El niño on nitrous oxide distribution in the eastern tropical south pacific. Biogeosciences16 (9), 2079–2093. doi: 10.5194/bg-16-2079-2019

  • CrossRef
  • Google Scholar
  • View reference in article

• 51

  JiQ.BabbinA. R.JayakumarA.OleynikS.WardB. B. (2015). Nitrous oxide production by nitrification and denitrification in the Eastern tropical south pacific oxygen minimum zone. Geophys. Res. Lett.42 (24), 10–755. doi: 10.1002/2015GL066853

  • CrossRef
  • Google Scholar
  • View reference in article

• 52

  JiQ.BuitenhuisE.SuntharalingamP.SarmientoJ. L.WardB. B. (2018). Global nitrous oxide production determined by oxygen sensitivity of nitrification and denitrification. Global Biogeochemical Cycles32 (12), 1790–1802. doi: 10.1029/2018GB005887

  • CrossRef
  • Google Scholar
  • View reference in article

• 53

  KeelingC. D. (1961). The concentration and isotopic abundances of carbon dioxide in rural and marine air. Geochimica Cosmochimica Acta24 (3-4), 277–298. doi: 10.1016/0016-7037(61)90023-0

  • CrossRef
  • Google Scholar
  • View reference in article

• 54

  KellyC. L.TravisN. M.BayaP. A.CasciottiK. L. (2021). Quantifying nitrous oxide cycling regimes in the Eastern tropical north pacific ocean with isotopomer analysis. Global Biogeochemical Cycles35 (2), e2020GB006637. doi: 10.1029/2020GB006637

  • CrossRef
  • Google Scholar
  • View reference in article

• 55

  KockA.Arévalo-MartínezD. L.LöscherC. R.BangeH. W. (2016). Extreme N₂O accumulation in the coastal oxygen minimum zone off Peru. Biogeosciences13 (3), 827–840. doi: 10.5194/bg-13-827-2016

  • CrossRef
  • Google Scholar
  • View reference in article

• 56

  KoolD. M.DolfingJ.WrageN.Van GroenigenJ. W. (2011). Nitrifier denitrification as a distinct and significant source of nitrous oxide from soil. Soil Biol. Biochem.43 (1), 174–178. doi: 10.1016/j.soilbio.2010.09.030

  • CrossRef
  • Google Scholar
  • View reference in article

• 57

  KozlowskiJ. A.StieglmeierM.SchleperC.KlotzM. G.SteinL. Y. (2016). Pathways and key intermediates required for obligate aerobic ammonia-dependent chemolithotrophy in bacteria and thaumarchaeota. ISME J.10 (8), 1836–1845. doi: 10.1038/ismej.2016.2

  • CrossRef
  • Google Scholar
  • View reference in article

• 58

  KraftB.JehmlichN.LarsenM.BristowL. A.KönnekeM.ThamdrupB.et al. (2022). Oxygen and nitrogen production by an ammonia-oxidizing archaeon. Science375 (6576), 97–100. doi: 10.1126/science.abe6733

  • CrossRef
  • Google Scholar
  • View reference in article

• 59

  Lazo-MurphyB. M.LarsonS.StainesS.BruckH.McHenryJ.BourbonnaisA.et al. (2022). Nitrous oxide production and isotopomer composition by fungi isolated from salt marsh sediments. Front. Mar. Sci.9, 2645. doi: 10.3389/fmars.2022.1098508

  • CrossRef
  • Google Scholar
  • View reference in article

• 60

  LehmannN.KienastM.GrangerJ.BourbonnaisA.AltabetM. A.TremblayJ.É. (2019). Remote western Arctic nutrients fuel remineralization in deep Baffin bay. Global Biogeochemical Cycles33 (6), 649–667. doi: 10.1029/2018GB006134

  • CrossRef
  • Google Scholar
  • View reference in article

• 61

  Lewicka-SzczebakD.AugustinJ.GiesemannA.WellR. (2017). Quantifying N₂O reduction to N₂ based on N₂O isotopocules–validation with independent methods (helium incubation and ¹⁵N gas flux method). Biogeosciences14 (3), 711–732. doi: 10.5194/bg-14-711-2017

  • CrossRef
  • Google Scholar
  • View reference in article

• 62

  LöscherC. R.KockA.KönnekeM.LaRocheJ.BangeH. W.SchmitzR. A. (2012). Production of oceanic nitrous oxide by ammonia-oxidizing archaea. Biogeosciences9 (7), 2419–2429. doi: 10.5194/bg-9-2419-2012

  • CrossRef
  • Google Scholar
  • View reference in article

• 63

  MaedaK.SporA.Edel-HermannV.HeraudC.BreuilM. C.BizouardF.et al. (2015). N₂O production, a widespread trait in fungi. Sci. Rep.5 (1), 9697. doi: 10.1038/srep09697

  • CrossRef
  • Google Scholar
  • View reference in article

• 64

  MagyarP. M. (2017). Insights into pathways of nitrous oxide generation from novel isotopologue measurements. Dissertation (Ph.D.), California Institute of Technology, 119 pp. doi: 10.7907/Z93776RJ

  • CrossRef
  • Google Scholar
  • View reference in article

• 65

  MariottiA.GermonJ. C.HubertP.KaiserP.LetolleR.TardieuxA.et al. (1981). Experimental determination of nitrogen kinetic isotope fractionation: some principles; illustration for the denitrification and nitrification processes. Plant Soil62 (3), 413–430. doi: 10.1007/BF02374138

  • CrossRef
  • Google Scholar
  • View reference in article

• 66

  MohnJ.WolfB.ToyodaS.LinC. T.LiangM. C.BrüggemannN.et al. (2014). Interlaboratory assessment of nitrous oxide isotopomer analysis by isotope ratio mass spectrometry and laser spectroscopy: current status and perspectives. Rapid Commun. Mass Spectrometry28 (18), 1995–2007. doi: 10.1002/rcm.6982

  • CrossRef
  • Google Scholar
  • View reference in article

• 67

  MonrealP. J.KellyC. L.TravisN. M.CasciottiK. L. (2022). Identifying the sources and drivers of nitrous oxide accumulation in the eddy-influenced Eastern tropical north pacific oxygen-deficient zone. Glob. Biogeochem. Cycles36 (6), e2022GB007310. doi: 10.1007/BF02374138

  • CrossRef
  • Google Scholar
  • View reference in article

• 68

  MuggeoV. M. (2003). Estimating regression models with unknown break‐points. Stat. Med.22(19), 3055–3071. doi: 10.1002/sim.1545

  • CrossRef
  • Google Scholar
  • View reference in article

• 69

  MuggeoV. M.MuggeoM. V. M. (2017). Package ‘segmented’. Biometrika58 (525-534), 516.

  • Google Scholar
  • View reference in article

• 70

  NevisonC.ButlerJ. H.ElkinsJ. W. (2003). Global distribution of N₂O and the ΔN₂O-AOU yield in the subsurface ocean. Global Biogeochemical Cycles17 (4). doi: 10.1029/2003GB002068

  • CrossRef
  • Google Scholar
  • View reference in article

• 71

  NevisonC.HollandE. (1997). A reexamination of the impact of anthropogenically fixed nitrogen on atmospheric N₂O and the stratospheric O₃ layer. J. Geophys. Res.: Atmospheres102 (D21), 25519–25536. doi: 10.1029/97JD02391

  • CrossRef
  • Google Scholar
  • View reference in article

• 72

  OstromN. E.OstromP. H. (2012). “The isotopomers of nitrous oxide: analytical considerations and application to resolution of microbial production pathways,” in Handbook of environmental isotope geochemistry (Berlin, Heidelberg: Springer), pp453–pp476. doi: 10.1029/2006JG000287

  • CrossRef
  • Google Scholar
  • View reference in article

• 73

  OstromN. E.PittA.SutkaR.OstromP. H.GrandyA. S.HuizingaK. M.et al. (2007). Isotopologue effects during N₂O reduction in soils and in pure cultures of denitrifiers. J. Geophys. Res.: Biogeosciences112 (G2). doi: 10.1029/2006JG000287

  • CrossRef
  • Google Scholar
  • View reference in article

• 74

  PatakiD. E.EhleringerJ. R.FlanaganL. B.YakirD.BowlingD. R.StillC. J.et al. (2003). The application and interpretation of keeling plots in terrestrial carbon cycle research. Global Biogeochem. Cycles17 (1). doi: 10.1029/2001GB001850

  • CrossRef
  • Google Scholar
  • View reference in article

• 75

  PoppB. N.WestleyM. B.ToyodaS.MiwaT.DoreJ. E.YoshidaN.et al. (2002). Nitrogen and oxygen isotopomeric constraints on the origins and sea-to-air flux of N₂O in the oligotrophic subtropical north pacific gyre. Global Biogeochemical Cycles16 (4), 12–11. doi: 10.1029/2001GB001806

  • CrossRef
  • Google Scholar
  • View reference in article

• 76

  RavishankaraA. R.DanielJ. S.PortmannR. W. (2009). Nitrous oxide (N₂O): the dominant ozone-depleting substance emitted in the 21^st century. Science326 (5949), 123–125. doi: 10.1126/science.1176985

  • CrossRef
  • Google Scholar
  • View reference in article

• 77

  RyabenkoE.KockA.BangeH. W.AltabetM. A.WallaceD. W. R. (2012). Contrasting biogeochemistry of nitrogen in the Atlantic and pacific oxygen minimum zones. Biogeosciences9, 203–215. doi: 10.5194/bg-9-203-2012

  • CrossRef
  • Google Scholar
  • View reference in article

• 78

  SantoroA. E.BuchwaldC.McIlvinM. R.CasciottiK. L. (2011). Isotopic signature of N₂O produced by marine ammonia-oxidizing archaea. Science333 (6047), 1282–1285. doi: 10.1126/science.1208239

  • CrossRef
  • Google Scholar
  • View reference in article

• 79

  SantoroA. E.CasciottiK. L.FrancisC. A. (2010). Activity, abundance and diversity of nitrifying archaea and bacteria in the central California current. Environ. Microbiol.12 (7), 1989–2006. doi: 10.1111/j.1462-2920.2010.02205.x

  • CrossRef
  • Google Scholar
  • View reference in article

• 80

  SchmidtH. L.WernerR. A.YoshidaN.WellR. (2004). Is the isotopic composition of nitrous oxide an indicator for its origin from nitrification or denitrification? a theoretical approach from referred data and microbiological and enzyme kinetic aspects. Rapid Commun. Mass Spectrometry18 (18), 2036–2040. doi: 10.1002/rcm.1586

  • CrossRef
  • Google Scholar
  • View reference in article

• 81

  SigmanD. M.GrangerJ.DiFioreP. J.LehmannM. M.HoR.CaneG.et al. (2005). Coupled nitrogen and oxygen isotope measurements of nitrate along the eastern north pacific margin. Global Biogeochemical Cycles19 (4). doi: 10.1029/2005GB002458

  • CrossRef
  • Google Scholar
  • View reference in article

• 82

  SniderD. M.VenkiteswaranJ. J.SchiffS. L.SpoelstraJ. (2012). Deciphering the oxygen isotope composition of nitrous oxide produced by nitrification. Global Change Biol.18 (1), 356–370. doi: 10.1111/j.1365-2486.2011.02547.x

  • CrossRef
  • Google Scholar
  • View reference in article

• 83

  SonnerupR. E.MeckingS.BullisterJ. L.WarnerM. J. (2015). Transit time distributions and oxygen utilization rates from chlorofluorocarbons and sulfur hexafluoride in the southeast pacific ocean. J. Geophys. Res. Ocean.120, doi:10.1002/2015JC010781. doi: 10.1002/2015JC010781

  • CrossRef
  • Google Scholar
  • View reference in article

• 84

  StanleyR. H. R.DoneyS. C.JenkinsW. J.LottD. E. (2012). Apparent oxygen utilization rates calculated from tritium and helium-3 profiles at the Bermuda Atlantic time-series study site. Biogeosciences9 (6), 1969–1983. doi: 10.5194/bg-9-1969-2012

  • CrossRef
  • Google Scholar
  • View reference in article

• 85

  StieglmeierM.MooshammerM.KitzlerB.WanekW.Zechmeister-Boltenste. rnS.RichterA.et al. (2014). Aerobic nitrous oxide production through n-nitrosating hybrid formation in ammonia-oxidizing archaea. ISME J.8 (5), 1135–1146. doi: 10.1038/ismej.2013.220

  • CrossRef
  • Google Scholar
  • View reference in article

• 86

  StövenT.TanhuaT.HoppemaM.BullisterJ. L. (2015). Perspectives of transient tracer applications and limiting cases. Ocean Sci.11 (5), 699–718. doi: 10.5194/os-11-699-2015

  • CrossRef
  • Google Scholar
  • View reference in article

• 87

  StrammaL.BangeH. W.CzeschelR.LorenzoA.FrankM. (2013). On the role of mesoscale eddies for the biological productivity and biogeochemistry in the eastern tropical Pacific Ocean off Peru. Biogeosciences10 (11), 7293–7306. doi: 10.5194/bg-10-7293-2013

  • CrossRef
  • Google Scholar
  • View reference in article

• 88

  SutkaR. L.OstromN. E.OstromP. H.BreznakJ. A.GandhiH.PittA. J.et al. (2006). Distinguishing nitrous oxide production from nitrification and denitrification on the basis of isotopomer abundances. Appl. Environ. Microbiol.72 (1), 638–644. doi: 10.1128/AEM.72.1.638-644.2006

  • CrossRef
  • Google Scholar
  • View reference in article

• 89

  SutkaR. L.OstromN. E.OstromP. H.GandhiH.BreznakJ. A. (2003). Nitrogen isotopomer site preference of N₂O produced by nitrosomonas europaea and methylococcus capsulatus bath. Rapid Commun. Mass Spectrometry17 (7), 738–745. doi: 10.1002/rcm.968

  • CrossRef
  • Google Scholar
  • View reference in article

• 90

  SutkaR. L.OstromN. E.OstromP. H.PhanikumarM. S. (2004). Stable nitrogen isotope dynamics of dissolved nitrate in a transect from the north pacific subtropical gyre to the Eastern tropical north pacific. Geochimica Cosmochimica Acta68 (3), 517–527. doi: 10.1016/S0016-7037(03)00483-6

  • CrossRef
  • Google Scholar
  • View reference in article

• 91

  TianH.XuR.CanadellJ. G.ThompsonR. L.WiniwarterW.SuntharalingamP.et al. (2020). A comprehensive quantification of global nitrous oxide sources and sinks. Nature586, 248–256. doi: 10.1038/s41586-020-2780-0

  • CrossRef
  • Google Scholar
  • View reference in article

• 92

  ToyodaS.MutobeH.YamagishiH.YoshidaN.TanjiY. (2005). Fractionation of N₂O isotopomers during production by denitrifier. Soil Biol. Biochem.37 (8), 1535–1545. doi: 10.1016/j.soilbio.2005.01.009

  • CrossRef
  • Google Scholar
  • View reference in article

• 93

  ToyodaS.YoshidaN.MiwaT.MatsuiY.YamagishiH.TsunogaiU.et al. (2002). Production mechanism and global budget of N₂O inferred from its isotopomers in the western north pacific. Geophys. Res. Lett.29 (3), 7–1. doi: 10.1029/2001GL014311

  • CrossRef
  • Google Scholar
  • View reference in article

• 94

  ToyodaS.YoshidaO.YamagishiH.FujiiA.YoshidaN.WatanabeS. (2019). Identifying the origin of nitrous oxide dissolved in deep ocean by concentration and isotopocule analyses. Sci. Rep.9 (1), 1–9. doi: 10.1038/s41598-019-44224-0

  • CrossRef
  • Google Scholar
  • View reference in article

• 95

  TrimmerM.ChronopoulouP. M.MaanojaS. T.Upstill-GoddardR. C.KitidisV.PurdyK. J. (2016). Nitrous oxide as a function of oxygen and archaeal gene abundance in the north pacific. Nat. Commun.7 (1), 1–10. doi: 10.1038/ncomms13451

  • CrossRef
  • Google Scholar
  • View reference in article

• 96

  WankelS. D.ZiebisW.BuchwaldC.CharoenpongC.de BeerD.DentingerJ.et al. (2017). Evidence for fungal and chemodenitrification based N₂O flux from nitrogen impacted coastal sediments. Nat. Commun.8 (1), 1–11. doi: 10.1038/ncomms15595

  • CrossRef
  • Google Scholar
  • View reference in article

• 97

  WardB. B. (2008). “Nitrification in marine systems,” in Nitrogen in the marine environment. Eds. CaponeD. G.BronkD. A.MulhollandM. R.CarpenterE. J. (Elsevier, San Diego, California: Academic Press), pp199–pp261.

  • Google Scholar
  • View reference in article

• 98

  WaughD. W.HallT. M.HaineT. W. N. (2003). Relationships among tracer ages. J. Geophys. Research-Oceans108 (C5), 16. doi: 10.1029/2002jc001325

  • CrossRef
  • Google Scholar
  • View reference in article

• 99

  WeissR. F.PriceB. A. (1980). Nitrous oxide solubility in water and seawater. Mar. Chem.8 (4), 347–359. doi: 10.1016/0304-4203(80)90024-9

  • CrossRef
  • Google Scholar
  • View reference in article

• 100

  WestleyM. B.PoppB. N.RustT. M. (2007). The calibration of the intramolecular nitrogen isotope distribution in nitrous oxide measured by isotope ratio mass spectrometry. Rapid Commun. Mass Spectrometry21 (3), 391–405. doi: 10.1002/rcm.2828

  • CrossRef
  • Google Scholar
  • View reference in article

• 101

  WestleyM. B.YamagishiH.PoppB. N.YoshidaN. (2006). Nitrous oxide cycling in the black Sea inferred from stable isotope and isotopomer distributions. Deep Sea Res. Part II: Topical Stud. Oceanogr.53, 1802–1816. doi: 10.1016/j.dsr2.2006.03.012

  • CrossRef
  • Google Scholar
  • View reference in article

• 102

  WrageN.VelthofG. L.Van BeusichemM. L.OenemaO. (2001). Role of nitrifier denitrification in the production of nitrous oxide. Soil Biol. Biochem.33 (12-13), 1723–1732. doi: 10.1016/S0038-0717(01)00096-7

  • CrossRef
  • Google Scholar
  • View reference in article

• 103

  YamagishiH.WestleyM. B.PoppB. N.ToyodaS.YoshidaN.WatanabeS.et al. (2007). Role of nitrification and denitrification on the nitrous oxide cycle in the eastern tropical north pacific and gulf of California. J. Geophys. Res.: Biogeosciences112 (G2). doi: 10.1029/2006JG000227

  • CrossRef
  • Google Scholar
  • View reference in article

• 104

  YamagishiH.YoshidaN.ToyodaS.PoppB. N.WestleyM. B.WatanabeS. (2005). Contributions of denitrification and mixing on the distribution of nitrous oxide in the north pacific. Geophys. Res. Lett.32 (4). doi: 10.1029/2004GL021458

  • CrossRef
  • Google Scholar
  • View reference in article

• 105

  YangS.ChangB. X.WarnerM. J.WeberT. S.BourbonnaisA. M.SantoroA. E.et al. (2020). Global reconstruction reduces the uncertainty of oceanic nitrous oxide emissions and reveals a vigorous seasonal cycle. Proc. Natl. Acad. Sci.117 (22), 11954–11960. doi: 10.1073/pnas.1921914117

  • CrossRef
  • Google Scholar
  • View reference in article

• 106

  YoshidaN.MorimotoH.HiranoM.KoikeI.MatsuoS.WadaE.et al. (1989). Nitrification rates and ¹⁵N abundances of N₂O and NO₃^- in the western north pacific. Nature342 (6252), 895–897.

  • Google Scholar
  • View reference in article

• 107

  YoshinariT.KnowlesR. (1976). Acetylene inhibition of nitrous oxide reduction by denitrifying bacteria. Biochem. Biophys. Res. Commun.69 (3), 705–710. doi: 10.1016/0006-291X(76)90932-3

  • CrossRef
  • Google Scholar
  • View reference in article