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.

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).

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.

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.

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.

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.

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).

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).

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).

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).

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.

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.

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).

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).

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).

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.

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.

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).

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.