# Revisiting the distribution of oceanic N2 fixation and estimating diazotrophic contribution to marine production

Weiyi Tang, Seaver Wang, Debany Fonseca-Batista, Frank Dehairs, Scott

Gifford, Aridane G. Gonzalez, Morgane Gallinari, Hélène Planquette,

Géraldine Sarthou, Nicolas Cassar

► To cite this version:

Weiyi Tang, Seaver Wang, Debany Fonseca-Batista, Frank Dehairs, Scott Gifford, et al.. Revisiting the distribution of oceanic N2 fixation and estimating diazotrophic contribution to marine production. Nature Communications, 2019, 10, pp.831. 10.1038/s41467-019-08640-0 . hal-02322624

HAL Id: hal-02322624

https://hal.science/hal-02322624v1

Submitted on 21 Oct 2019

**HAL** is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
ARTICLE

https://doi.org/10.1038/s41467-019-08640-0

OPEN

# Revisiting the distribution of oceanic N₂ fixation and estimating diazotrophic contribution to marine production

Weiyi Tang<fnref n="1" />, Seaver Wang ID <fnref n="1" />, Debany Fonseca-Batista²,<fnref n="6" />, Frank Dehairs², Scott Gifford³, Aridane G. Gonzalez ID ⁴,⁵, Morgane Gallinari⁴, Hélène Planquette⁴, Géraldine Sarthou⁴ & Nicolas Cassar ID <fnref n="1" />,⁴

Marine N₂ fixation supports a significant portion of oceanic primary production by making N₂ bioavailable to planktonic communities, in the process influencing atmosphere-ocean carbon fluxes and our global climate. However, the geographical distribution and controlling factors of marine N₂ fixation remain elusive largely due to sparse observations. Here we present unprecedented high-resolution underway N₂ fixation estimates across over 6000 kilometers of the western North Atlantic. Unexpectedly, we find increasing N₂ fixation rates from the oligotrophic Sargasso Sea to North America coastal waters, driven primarily by cyano-bacterial diazotrophs. N₂ fixation is best correlated to phosphorus availability and chlorophyll-a concentration. Globally, intense N₂ fixation activity in the coastal oceans is validated by a meta-analysis of published observations and we estimate the annual coastal N₂ fixation flux to be 16.7 Tg N. This study broadens the biogeography of N₂ fixation, highlights the interplay of regulating factors, and reveals thriving diazotrophic communities in coastal waters with potential significance to the global nitrogen and carbon cycles.

<footnote n="1">¹Division of Earth and Ocean Sciences, Nicholas School of the Environment, Duke University, Durham, NC 27708, USA. ²Analytical, Environmental and Geo-Chemistry, Vrije Universiteit Brussel, 1050 Brussels, Belgium. ³Department of Marine Sciences, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA. ⁴Laboratoire des Sciences de l'Environnement Marin (LEMAR), UMR 6539 UBO CNRS IRD IFREMER, Institut Universitaire Européen de la Mer (IUEM), 29280 Brest, France. ⁵Instituto de Oceanografía y Cambio Global (IOCAG), Universidad de Las Palmas de Gran Canaria, 35214 Las Palmas, Spain.</footnote>

<footnote n="6">⁶Present address: Department of Biology, Dalhousie University, Halifax, NS B3H 4R2, Canada. Correspondence and requests for materials should be addressed to N.C. (email: Nicolas.Cassar@duke.edu)</footnote>
A pproximately half of global primary production occurs in the oceans¹. In the vast expanse of the oligotrophic oceans, marine primary production is limited by nitrogen². In these regions, nitrogen (N₂) fixation by diazotrophs has been hypothesized to be an important source of new nitrogen ultimately influencing the uptake and sequestration of CO₂³⁻⁵. For instance, in the North Atlantic subtropical gyre, a significant seasonal carbon drawdown in the absence of measurable nutrients has been attributed to episodic and patchy N₂ fixation events⁶,⁷. These events are, however, difficult to capture using current methods, which rely on discrete sampling. Furthermore, most observations to date have been collected in tropical and subtropical open oceans, overlooking the potential role of coastal regions⁸,⁹. The limited number of observations impede our ability to close regional marine nitrogen budgets, scale estimates globally, and identify factors controlling N₂ fixation¹⁰.

Marine N₂ fixation is generally believed to be regulated by various factors including light, temperature, nutrients, and trace-metal availability¹¹. While low light, low temperature, and high bioavailable nitrogen have traditionally been assumed to limit N₂ fixation, recent studies have reported significant diazotrophic activities in darker, colder, and more nitrogen-rich environments, thereby broadening the putative biogeography of marine N₂ fixation¹²⁻¹⁶. The discovery of these new niches has been accompanied by increasing appreciation for the large diversity of species fixing N₂, ranging from the well-known *Trichodesmium* and diatom-diazotroph associations (DDA) to more recently-recognized unicellular cyanobacteria¹⁷ and non-cyanobacterial diazotrophs¹⁸.

Taken together, the large uncertainty in regional and global budgets of marine N₂ fixation, associated with its patchy and recently broadened biogeography, demands new tools to adequately map this important biogeochemical process. We recently deployed a new method across more than 6000 km of the North Atlantic to revisit the geographical distribution and assess the controlling factors of N₂ fixation. Our near-real-time, continuous high-resolution measurements allowed us to locate hotspots of N₂ fixation and adapt our sampling strategy to characterize plankton communities and environmental properties¹⁹ (Methods). The contribution of N₂ fixation to net community production (NCP) was simultaneously evaluated at high resolution. Our new observations reveal hotspots of marine N₂ fixation along the Eastern Seaboard, and highlight the overlooked significance of N₂ fixation to coastal and global nitrogen and carbon cycling.

## Results and Discussion

**N₂ fixation distribution and controlling factors.** Our survey revealed substantial variability and diel cycling behavior in surface N₂ fixation rates, which ranged from less than 0.01 to nearly 15 nmol N L⁻¹ h⁻¹ (Fig. 1a and Supplementary Fig. 1). When integrated over 24-h N₂ fixation diel cycles, continuous estimates of daily surface N₂ fixation rates (≤0.19–97.6 nmol N L⁻¹ d⁻¹) were in line with discrete N₂ fixation rates concurrently determined by ¹⁵N₂ incubations (n = 7, r = 0.97, p < 0.01, Fig. 1b and Methods). To extrapolate surface measurements to the entire euphotic zone, we derived an empirical relationship between surface and depth-integrated N₂ fixation rates (Supplementary Fig. 2). Our lower-end measurements of less than 10 µmol N m⁻² d⁻¹ are within the range of published rates near Bermuda²⁰ (Fig. 1c). In contrast, high N₂ fixation rates reaching 3000 µmol N m⁻² d⁻¹ near the New Jersey coast are among the top 2% of rates ever reported in the global ocean²¹, further underscoring the high N₂ fixation along the continental shelf of the eastern seaboard²², and the value of high-frequency observations for identifying hotspots. The coastal region (bathymetry ≥ −200 m)

stands in sharp contrast to open ocean areas, with depth-integrated N₂ fixation rates in the coastal sectors (geometric mean of 577 µmol N m⁻² d⁻¹) being on average an order of magnitude larger than open ocean rates (geometric mean of 85 µmol N m⁻² d⁻¹).

Sea surface temperature did not appear to drive the spatial variability of N₂ fixation in our study area (Supplementary Fig. 3 and Supplementary Fig. 4). In open ocean regions, N₂ fixation rates varied substantially even when the temperature range was narrow. Moreover, high N₂ fixation rates were observed at temperatures ranging from 23 °C off New Jersey to 30 °C near the Florida coast. Fixed nitrogen is known to suppress N₂ fixation, but the threshold of inhibition differs among diazotrophs and can be fairly high²³. Observed dissolved inorganic nitrogen did not effectively regulate N₂ fixation rates (Supplementary Fig. 4). We also note that a recent study shows that these coastal waters may be nitrogen limited in summer²⁴. In addition, the excess of nitrogen in subsurface waters (N* = [NO₃⁻] - [PO₄³⁻] × 16), commonly used as a geochemical proxy for the distribution of N₂ fixation²⁵, was not a strong predictor of overlying N₂ fixation rates (Fig. 1c and Supplementary Fig. 4). This comparison should be interpreted with caution, as N* integrates over longer spatial and temporal scales than our observations. N* is also not well resolved in coastal waters and may be affected by other processes such as atmospheric nitrogen deposition²⁶. In contrast, some of the N₂ fixation hotspots coincided with high phosphorus concentration and excess phosphorus (P* = [PO₄³⁻] - [NO₃⁻]/16) at the ocean surface (Fig. 1b and Supplementary Fig. 4). Regions of N₂ fixation have been hypothesized to be coupled to areas of denitrification via the upwelling of waters deficient in nitrogen relative to phosphorus²⁷. The phytoplankton bloom near the New Jersey coast, where the highest N₂ fixation rates were observed, recurs almost every summer and may be associated with local upwelling²⁸. High N₂ fixation rates have also been reported in other upwelling systems worldwide, including the equatorial Atlantic Ocean²⁹, the northwest African coastal upwelling³⁰, and the Benguela Upwelling System³¹. The excess phosphorus may also result from terrestrial and/or riverine runoff. For example, N₂ fixation and carbon sequestration in the tropical North Atlantic were shown to be enhanced by the Amazon River plume³².

A recent study in the Eastern South Pacific suggested Fe, rather than phosphorus, may limit N₂ fixation³³, questioning the spatial coupling between N₂ fixation and denitrification³⁴. However, we did not find a strong relationship between dissolved Fe and N₂ fixation rates across our study (Supplementary Fig. 4). Dissolved Fe ranged from 0.5 nmol L⁻¹ near Bermuda to around 1.8 nmol L⁻¹ along the Florida coast, which is higher than Fe measured in the Eastern South Pacific. Fe concentration is admittedly a poor predictor of Fe availability, as concentrations merely reflect snapshots of the complex interactions between sources and sinks³⁵. Consequently, we calculated Fe* (Fe* = Fe - R<sub>Fe</sub>PO₄; where R<sub>Fe</sub> = 0.47 mmol Fe : 1 mol PO₄) to evaluate whether Fe potentially limits phytoplankton growth³⁶. Positive values of Fe* across the study area indicate Fe was not limiting. The coupling between N₂ fixation and N loss may be a dominant factor in regions where Fe is abundant, notably in coastal oceans. We hypothesize that, in contrast to the Eastern South Pacific, North American coastal waters support substantial N₂ fixation due to high sedimentary nitrogen loss³⁷ and high Fe input (e.g., from sediment and atmospheric dust deposition)³⁸ (Supplementary Fig. 5). Interestingly, N₂ fixation rates correlated well with dissolved manganese (Mn) concentrations. While the correlation of N₂ fixation to Mn could be coincidental or symptomatic of other factors, it deserves further investigation as the physiological requirement for Mn in marine diazotrophs is poorly characterized.
<figure><img src="image_2.png" /><figcaption>**Fig. 1** N₂ fixation rates (NF) measured and calculated in August 2015 and 2016 over the western North Atlantic Ocean. **a** In situ hourly surface volumetric N₂ fixation rates determined by the continuous underway incubation method, overlaid on chlorophyll-a concentrations measured by the MODIS satellite during the respective cruise periods (Aug 3-12 in both 2015 and 2016). Note that chlorophyll-a concentrations are shown on a logarithmic scale.</figcaption></figure>

**b** Calculated daily surface volumetric N₂ fixation rates overlaid on the August climatology of surface excess phosphate P* (P* = [PO₄³⁻] - [NO₃⁻]/16). Nutrient data were obtained from the World Ocean Atlas 2013 version 2. Red vertical bars represent N₂ fixation rates determined by the discrete dissolved ¹⁵N₂ incubation method during the 2015 cruise. **c** Estimated depth-integrated N₂ fixation rates overlaid on the August climatology of subsurface excess nitrogen N* (N* = [NO₃⁻] - [PO₄³⁻] × 16).

Our N₂ fixation measurements strongly correlated to satellite estimates of chlorophyll-a concentrations ([Chl]) (Fig. 1a and Supplementary Fig. 4). This is unexpected as N₂ fixation is generally believed to be most significant where nitrogen is limited, such as the low biomass regions of the subtropical gyres. The low-[Chl] waters of the Sargasso Sea, typically viewed as N₂ fixation hotspots³⁹, exhibited lower N₂ fixation rates than those measured in the Mid-Atlantic Bight. This pattern was further supported by a meta-analysis, which showed that N₂ fixation is correlated to [Chl] in the global ocean (Supplementary Fig. 6). This pattern may be related to the stimulation of non-autotrophic N₂ fixation by organic matter in high [Chl] waters⁴⁰. However, while our field observations identify phosphorus, P*, and [Chl] as predictors of spatial variations in N₂ fixation in our study area, a meta-analysis of published results shows that none of the putative regulating factors of N₂ fixation can satisfactorily explain variations in volumetric rates globally (Supplementary Fig. 6). We posit that

N₂ fixation is likely driven by a complex interplay of spatially-variable environmental factors, also reflecting the heterogeneity and the large diversity of marine diazotrophs and their niches.

**Distribution of diazotrophic phylotypes.** Diazotrophs and their potential hosts were identified via high-throughput quantitative 16S rRNA and 18S rRNA gene sequencing from our 2015 cruise (Methods). Although the 16S rRNA gene approach differs from the *nifH* method for characterizing diazotrophs in terms of specificity and coverage⁴¹, it provides some insights into the broad distribution of diazotrophs. To address the cases where organisms may not be capable of N₂ fixation despite sharing a similar 16S rRNA gene with diazotrophs, we only searched for diazotrophs known to fix N₂ among our 16S rRNA gene sequences. Distinct diazotrophic communities were found to dominate in different ecological domains (Fig. 2a). Heterotrophic groups, which include
<figure><img src="image_3.png" /><figcaption>**Fig. 2** Absolute rRNA gene abundances of various diazotrophs and their potential hosts in surface seawater collected during the 2015 cruise. The inset map shows the locations where molecular samples were collected. The distributions of diazotrophs (a) and hosts (b) are compared with surface daily N₂ fixation rates (solid red line). Diazotrophs and hosts below detection limits at different sampling locations are not shown in the figure. "Others" in panel a includes diazotrophic proteobacteria belonging to *Bradyrhizobium*, *Mesorhizobium*, *Novosphingbium*, and *Paenibacillus*</figcaption></figure>

<figure><img src="image_4.png" /><figcaption>**Fig. 3** Distribution of Net Community Production (NCP) and contribution of N₂ fixation to NCP. **a** Color-coded NCP with negative values shown in grey color scale. **b** Distribution of the ratio of N₂ fixation-supported carbon production to NCP. Circle size in **b** represents the magnitude of NCP, contextualizing the large uncertainty in calculated ratios when the denominator (i.e., NCP) is small</figcaption></figure>

members known to be diazotrophs⁴², were more abundant than diazotrophic cyanobacteria in the open ocean, where N₂ fixation was relatively low. These observations are consistent with recent recognition of the widespread distribution of non-cyanobacterial diazotrophs, whose activities remain poorly constrained¹⁸. *Trichodesmium* peaked off South Carolina. *Richelia intracellularis* showed relatively high abundances in the oligotrophic open ocean and peaked in coastal waters, where its hosts—*Rhizosolenia* and *Hemiaulus*—were also found at relatively high abundances (Fig. 2b). The most striking N₂ fixation hotspot off the New Jersey coast was likely driven by a UCYN-A bloom that reached 2 × 10⁷ 16S rRNA gene copies L⁻¹, which is of comparable magnitude to UCYN-A abundances (2.5–3.5 × 10⁷ L⁻¹) observed previously in the same region using the nifH method²². One of UCYN-A's hypothesized hosts, *Braarudosphaera bigelowii* also flourished in this region. Across all samples, the ratio of UCYN-A (16S rRNA gene):*B. bigelowii* (18S rRNA gene) varied from 0.1 to 35 with a median of 0.24. The two organisms co-occurred in most

samples, consistent with previous studies that suggest obligate symbiosis⁴³,⁴⁴. There is a growing interest in UCYN-A's unusual physiological and ecological traits⁴⁵. Its genetic diversity, evolution, and the niches inhabited by its different lineages deserve further investigation. Meanwhile, its presence in the coastal oceans provides new opportunities to study this unique organism. The divergent geographic distribution of different diazotrophs in our research area likely reflects their respective niches. For example, warm seawater is more favorable to *Trichodesmium*⁴⁶, while UCYN-A prefers temperate environments¹⁷. This difference may partly explain why *Trichodesmium* abundances peaked off the South Carolina coast (~30 °C) while UCYN-A dominated in regions with lower temperatures (~23 °C). Overall, the quantitative 16S rRNA and 18S rRNA gene sequencing methods revealed spatial patterns of diazotrophs and their hosts despite the assumptions of our quantitative sequencing methods, e.g., equivalent recovery efficiency for both the standard and the natural sequences in the sample (Methods).
<figure><img src="image_5.png" /><figcaption>**Fig. 4** Histograms comparing N₂ fixation rates in open and coastal oceans. **a** The magnitude of depth-integrated coastal N₂ fixation rates (red) is slightly larger than N₂ fixation rates in the open oceans (blue). Volumetric N₂ fixation rates at all depths (**b**) and at the surface (**c**) are generally higher in coastal regions than in the open ocean. **d** Integration depths used to calculate the depth-integrated N₂ fixation rates from the volumetric N₂ fixation rates. Geometric means and geometric standard deviations are shown in the legends</figcaption></figure>

**Variable contribution of N₂ fixation to new production.** N₂ fixation has been estimated to be an important source of new nitrogen in oligotrophic waters, supporting as much as 50% of new production⁴,⁴⁷, yet the contribution of N₂ fixation to productivity in coastal oceans remains relatively understudied⁴⁸. To assess the proportion of production fueled by N₂ fixation, we combined N₂ fixation observations with underway estimates of NCP based on high-frequency measurements of the dissolved O₂/Ar ratio (Methods). NCP was mostly positive, with higher rates along the North American coast where high [Chl] and N₂ fixation rates were also observed (Fig. 3a). NCP was relatively low in the open ocean likely due to nutrient limitation. As a rough estimate, we converted N₂ fixation rates and NCP to their carbon equivalents using a theoretical C:N:O₂ stoichiometry of 106:16:138. Regional differences in stoichiometry would modify but not erase the large gradients of N₂ fixation and NCP observed over our cruise transects. We found that the contribution of N₂ fixation to NCP varied substantially over the western North Atlantic (Fig. 3b). Across large portions of the oligotrophic subtropical ocean, no more than 20% NCP was generally fueled by N₂ fixation (with some high excursions). Other mechanisms of nutrients supply, such as revised estimates of vertical nitrate flux⁴⁹, must therefore be invoked to support the NCP rates we observed and those that have been reported in the Sargasso Sea. In contrast, the ratio of N₂ fixation to NCP exceeded 50% in some regions off the Cape Hatteras and New Jersey coasts (Fig. 3b). The high contribution of N₂ fixation to primary production off the New Jersey coast is supported by dual-tracer ¹⁵N₂ and ¹³C incubations (Supplementary Fig. 7). Despite methodological differences

between these techniques, both methods independently capture similar spatial patterns of contribution of N₂ fixation to biological carbon fixation (see NCP estimates in Methods). Our results highlight that N₂ fixation is not only high in coastal regions, but may also contribute significantly to marine production.

**Updated N₂ fixation budget and global implications.** To contextualize our findings, we performed a meta-analysis combining ¹⁵N₂ incubations collected during our cruises with discrete N₂ fixation measurements compiled from the literature, not including our underway continuous measurements (Supplementary Fig. 8). Our updated database contains over 80% more depth-integrated observations (1172 points in total) than the most up-to-date database currently available in the literature (630 points)²¹. Less than 15% of observations reported in the literature were collected in coastal waters. However, these observations support our findings of high N₂ fixation rates in the neritic environment (Fig. 4), notably on the eastern American coast⁵,²², eastern Arabian Sea⁵⁰, and estuaries of the Baltic Sea¹⁴. The similar magnitude of depth-integrated N₂ fixation rates in coastal and open oceans leads to significantly higher volumetric N₂ fixation rates in the coastal oceans due to shallower integration depths (Fig. 4 and Methods). However, recent reports of N₂ fixation in the deep ocean¹²,¹³ may reverse the pattern of marginally higher depth-integrated N₂ fixation rates in coastal waters, if N₂ fixation is integrated to the aphotic zone of the open ocean.

In light of these new observations, we reassessed the budgets of marine N₂ fixation globally as well as separately for the neritic
<table><caption>Table 1 N₂ fixation rates in the open and coastal ocean compared with previous studies.</caption><thead><tr><th>Estimates</th><th>Measurements</th><th>Area<br>(10<sup>6</sup> km<sup>2</sup>)</th><th>Range of rates<sup>a</sup><br>(µmol N m<sup>-2</sup> d<sup>-1</sup>)</th><th>Flux<br>(Tg N yr<sup>-1</sup>)</th></tr></thead><tbody><tr><td>Coastal ocean</td><td>143</td><td>24.5</td><td>1.2-397.3<sup>b</sup></td><td>16.7 (14.3)<sup>b</sup></td></tr><tr><td>Open ocean</td><td>857</td><td>290.2</td><td>1.2-160.2<sup>c</sup></td><td>6.6 (0.3)<sup>c</sup></td></tr><tr><td>Global sum</td><td>1000</td><td>314.7</td><td>7.5-288.9</td><td>179.4 (143.4)</td></tr><tr><td>Luo et al., 2012<sup>21</sup></td><td>630</td><td>288.5</td><td>7.4-90.9</td><td>64.2 (2.2)</td></tr><tr><td>Großkopf et al., 2012<sup>51</sup></td><td>335</td><td>215.1</td><td>1.2-397.3<br>1.2-160.2<br>13-590<br>7.9-120<br>19-115<sup>d</sup></td><td>196.1 (144.1)<br>70.8 (2.2)<br>137 (9.2)<br>62 (52-73)<br>177<sup>d</sup></td></tr></tbody></table>

and oceanic regions (Supplementary Fig. 9). Our updated geometric and arithmetic mean estimates of marine N₂ fixation for the global ocean at 70.8 Tg N yr⁻¹ and 196.1 Tg N yr⁻¹, respectively, are slightly higher than other estimates²¹,⁵¹ (Table 1 and Supplementary Table 1). Earlier studies report coastal N₂ fixation rates of 15 Tg N yr⁻¹, with most activity associated with benthic diazotrophs⁵². Our updated analysis shows that a significant portion of N₂ fixation is also occurring in the water column of coastal regions, contributing an additional 6.6 (geometric) or 16.7 (arithmetic) Tg N yr⁻¹ to the global budget. This nitrogen input could support the equivalent of 95 Tg C yr⁻¹ of primary production. These updated N₂ fixation budgets have large uncertainties since they are sensitive to the extrapolation method to scale the sparse data to the global ocean. Recent studies have also identified methodological issues in historical observations, which could lead to under⁵³ or over⁵⁴ estimations. However, the number of observations resulting from the revised dissolved ¹⁵N₂ incubation method is currently too limited to derive robust global estimates. Nevertheless, after accounting for both aquatic and sedimentary N₂ fixation, the coastal oceans may play a larger role in the nitrogen cycle than previously considered.

Using an underway method recently developed in our lab for continuous high-resolution N₂ fixation measurements, we mapped N₂ fixation at unprecedented scales across the western North Atlantic, identifying hotspots of N₂ fixation in the Mid-Atlantic Bight. Our study challenges the classic paradigm of N₂ fixation distribution and further underscores the central role coastal regions play in the global cycling of nutrients and carbon⁵⁵. With coastal regions being exposed to ever-increasing anthropogenic disturbances⁵⁶, expansion of coastal monitoring efforts using high-resolution methods will be critical to evaluate ongoing biogeochemical changes and their global repercussions.

## Methods

**Underway N₂ fixation rate measurements.** N₂ fixation rates were estimated at high resolution using a continuous underway method of Flow-through incubations for Acetylene-Reduction Assays by Cavity ring-down laser Absorption Spectroscopy (FARACAS). A description of the method is presented in Cassar et al. (2018)¹⁹ with a brief outline below. Nitrogenase activity in seawater is estimated based on the conventional technique of acetylene (C₂H₂) reduction to ethylene (C₂H₄)⁵⁷⁻⁵⁹. Instead of injecting C₂H₂ gas directly into the incubation bottle, C₂H₂ gas produced from high-purity calcium carbide (Acros Organics) is first dissolved in 0.2-µm filtered seawater that is collected from a trace-metal clean towfish (Geofish) to make a C₂H₂-H₂O tracer. The dissolved C₂H₂ approach was previously applied for measuring nitrogenase activity in estuarine sediments⁶⁰. The C₂H₂-H₂O tracer is then mixed at a constant ratio using a two-channel peristaltic pump (Masterflex) with a continuous stream of seawater supplied by the Geofish. The mixture of C₂H₂-H₂O tracer and seawater is continuously pumped into a 9-L flow-through incubation reactor (Chemglass). The flow-through incubation reactor (Chemglass) at a flow rate of ~100 mL min⁻¹. The short flow-through incubation, with an e-folding residence time of 90 min (i.e.

~63% of the seawater in the incubation reactor replaced in 90 min), minimizes the effects of C₂H₂ on metabolic processes and on microbial community structure⁶¹,⁶². The incubation reactor is lit by a strand of cool-light LEDs fitted with blue filters to simulate the light quality at the ocean surface. The light intensity is instantly calculated and adjusted based on the ship's location and local time. A water jacket on the incubation reactor is flushed with a high-flow rate of continuous surface seawater to mimic the in situ sea surface temperature. Downstream of the flow-through incubation reactor, the seawater flows into a gas extraction chamber. This gas extraction chamber consists of a glass frit with medium-size pores (Chemglass) and a custom-built gas-water separation system. A flow of 35 mL min⁻¹ of C₂H₄-free air controlled by a mass flow controller (OMEGA) continuously purges the incubated seawater, extracting ethylene out of the seawater, and carrying it to a Cavity Ring-Down laser absorption Spectrometer (CRDS, Picarro) for analyses. This CRDS ethylene analyzer measures ethylene concentrations in real time at ppb levels with high frequency and accuracy⁶³. Approximately every 3 h, the incubation reactor is bypassed to determine the background ethylene concentration in the mixture of C₂H₂-H₂O tracer and seawater. The difference between the incubation ethylene and background ethylene concentrations represents ethylene production rates during the incubation period. Finally, ethylene production rates are converted to N₂ fixation rates using a conversion factor of 4:1⁵⁸,⁶⁴⁻⁶⁶. We acknowledge that 4:1 is a theoretical ratio with uncertainties. However, our comparison of FARACAS to the ¹⁵N₂ addition method shows good agreement when applying this conversion factor¹⁹. In the current configuration, the detection limit of FARACAS is 0.19 nmol N L⁻¹ d⁻¹, which is also comparable to the ¹⁵N₂ addition method.

**Discrete N₂ fixation and primary production incubations.** For comparison to our underway survey of N₂ fixation, discrete ¹⁵N₂ incubation experiments in parallel with ¹³C additions were also conducted at eight stations during the 2015 cruise using methods detailed in previous studies³⁰,⁵¹. Seawater samples were collected from each station at three levels in the euphotic zone, including the surface (5 m), an intermediate depth above the Deep Chlorophyll Maximum (DCM), and the DCM. Four liters of seawater were immediately filtered onto glass microfiber filters (MGF, 0.7 µm, Sartorius) to determine natural concentrations and isotopic signatures of particulate organic carbon (POC) and particulate nitrogen (PN). For incubation experiments, 4.5-L Nalgene polycarbonate bottles were first partly filled with natural seawater. Then, ¹⁵N₂-enriched filtered seawater (98% ¹⁵N atom%, Euriotop, batch number 23/051301) and NaH¹³CO₃ solution (99%, Euriotop) were added into the incubation bottles, reaching approximate final enrichments of 2 ¹⁵N atom% and 7 ¹³C atom%. Finally, 4.5-L Nalgene bottles were topped off with natural seawater from sampled depths and capped with septum-fitted screw caps. Incubations were subsequently performed for 24 h in on-deck incubators covered by blue light filters simulating light intensity at the sampled depths. Incubators were flushed with surface seawater to avoid heating due to sunlight. Finally, incubated seawater was filtered onto MGF filters, which were stored at -20 °C until further analysis on land. Filters were treated and analyzed for POC, PN, δ¹³Cₚₚₐ, and δ¹⁵Nₚₙ using an Elemental Analyzer-Isotope Ratio Mass Spectrometer (EA-IRMS; EuroVector Euro EA 3000 coupled to a Delta V Plus, Thermo Scientific) to calculate corresponding carbon uptake and N₂ fixation rates. The N₂ fixation rates measured by our underway method closely match the results obtained from discrete incubation experiments (n = 7, r = 0.97, p < 0.01). A more comprehensive inter-method comparison can be found in Figure 5 of Cassar et al. (2018)¹⁹, showing good agreement between the two methods.

**Nutrients and trace-metal analyses.** Nutrient samples were collected from a CTD rosette equipped with 24 12-L Niskin bottles. Seawater was subsampled in acid-washed 15-mL polypropylene vials and immediately preserved at -20 °C.

<footnote><sup>a</sup>Range of mean N₂ fixation rates for different latitudinal bands. Mean N₂ fixation rates in latitudinal bands are shown in Supplementary Table 1</footnote>

<footnote><sup>b</sup>Arithmetic estimates</footnote>

<footnote><sup>c</sup>Geometric estimates in italics</footnote>

<footnote><sup>d</sup>Station weighted estimates. 177 Tg N yr⁻¹ is obtained after correcting for the underestimation due to incomplete dissolution of ¹⁵N₂ in incubations</footnote>

<footnote>Uncertainties are shown in parentheses</footnote>
Nitrate + Nitrite and phosphate were analyzed on land using an Automatic Nutrients Analyzer with detection limits of 0.03 µM and 0.014 µM, respectively.

For trace-metal analyses, all reagents, standards, and blanks were prepared in acid-cleaned low-density polyethylene (LDPE) or Teflon-fluorinated ethylene propylene (FEP) bottles. Bottles were cleaned following GEOTRACES protocols. Trace metal samples were collected in surface seawater (~5 m) using a towed fish (UCSC) deployed along side of the ship while underway⁶⁷. During stops, the towhfish was recovered from seawater to avoid contamination. Surface seawater was pumped through Teflon tubing to a sink located in a home-made clean plastic bubble installed within the chemistry lab on the ship. There, seawater was filtered in-line from the Teflon tubing outlet using 0.22-µm pore-size Acropak filter cartridges and collected in acid-washed 60-mL LDPE bottles that were triple-rinsed with ~20 mL of filtered seawater before final sample collection. Samples for dissolved trace metal were then acidified to pH = 2 with concentrated HCl (Ultrapur grade, Merck) under a laminar flow hood equipped with HEPA filter. Samples were then double-bagged and stored in the dark, at room temperature, until analysis. On land, all analyses were performed in cleanroom environments at the Pôle Spectrométrie Océans (Brest, France).

Seawater samples were introduced to a PFA-ST nebulizer and a cyclonic spray chamber via a SeaFASTpico introduction system (Elemental Scientific Incorporated, Omaha, NE), following the protocol of Lagerström et al. (2013)⁶⁸. High-purity grade solutions and water (Milli-Q, 18.2 MΩ cm) were used to prepare the following reagents on a daily basis. Buffer was made from 0.5 M acetic acid (Ultrapur grade, Merck) and 0.6 M ammonium hydroxide (Ultrapur grade, Merck) and was adjusted to pH = 8.3. Elution acid was made of 1.6 M HNO₃ (Ultrapur grade, Merck) in Milli-Q water and spiked with 1 µg mL⁻¹ In (PlasmaCAL calibration standards) to allow for drift correction. Autosampler and column rinsing solutions were made from 0.012 M HCl (Ultrapur grade, Merck) in Milli-Q water.

Mixed element standard solution was prepared gravimetrically using high-purity standards (Fe, Mn, Cd, Co, Zn, Cu, Pb; PlasmaCAL calibration standards) in 0.8 M HNO₃ (Ultrapur grade, Merck). A six-point calibration curve was prepared by standard additions of the mixed element standard to our in-house standard (North Atlantic filtered seawater, collected at 55.87445° N/48.09345° W, 40 m depth, 0.15 nM) and run at the beginning, the middle and the end of each run. Final concentrations of samples and procedural blanks were calculated from In-normalized data. Precision was assessed through replicate samples (every tenth sample was a replicate) and accuracy was determined from analysis of consensus seawater (SAFe S1 and D2, and GSP, GSC).

**Diazotrophic community structure analysis.** Diazotrophic phylotypes were identified and quantified using data obtained from 16S rRNA amplicon sequencing of environmental DNA, targeting the V4 region⁶⁹. Eukaryotic hosts of some diazotrophic taxa were similarly detected using amplicon sequencing of the V4 region of the 18S rRNA gene. Detailed experimental protocols are described in Wang et al. (2018) including sample collection, addition of internal controls for quantitative sequencing, DNA extraction, primer sequences, PCR amplification steps, and procedures for the analysis of sequencing data⁷⁰. This quantitative analysis has previously been described and applied in other environments⁷¹,⁷². Here, the processes are described briefly. From 0.2 to 1 L of seawater (average of 0.8 L) pumped from a towed fish were filtered onto a 0.22-µm polycarbonate filter using a peristaltic pump. The low-volume samples were typically collected in coastal waters, where high biomass led to clogging of filters. The volume filtered was recorded for each sample. The filter was flash-frozen immediately in liquid nitrogen and stored at -80 °C. Following the cruise, DNA extraction was performed using the Qiagen DNeasy Plant Mini Kit according to the manufacturer's instructions, with several modifications adapted from Moisander et al. (2008)⁷³. Prior to bead beating, 3.04 ng of *Thermus thermophilus* (ATCC #27634D-5) genomic DNA and 0.679 ng of *Schizosaccharomyces pombe* (ATCC #24843D-5) genomic DNA were added to each sample as internal DNA standards, each in 50 µL volumes. These additions introduced ~5,780,000 and 2,800,000 copies of *S. pombe* and *T. thermophilus* rRNA sequences sample⁻¹, amounts expected to constitute ≤1% of total reads sample⁻¹ following sequencing⁷⁰. PCR cycle parameters are detailed in Wang et al. (2018). Following PCR purification using the Qiagen QIAquick PCR Purification Kit, samples were pooled at equimolar concentrations. Illumina MiSeq sequencing (300 bp PE reads, V3 chemistry) was performed at the Sequencing and Genomic Technologies Shared Resources core facility at the Duke Center for Genomic and Computational Biology (Durham, USA). Raw rRNA sequences and metadata are available from the NCBI Sequence Read Archive under accession number SRP126177.

We used QIIME to process and analyze our Illumina sequencing data following the pipeline described in Fadrosh et al. (2014)⁷⁴,⁷⁵. Taxonomy tables reporting raw counts of 16S rRNA gene and 18S rRNA gene were produced by open-reference operational taxonomic unit (OTU) picking at the 97% threshold using the Usearch 6.1 algorithm and the SILVA ribosomal RNA database⁷⁶⁻⁷⁸. The SILVA ribosomal RNA database was supplemented with the addition of full length 16S rRNA gene sequences of UCYN-A1 and UCYN-A2 (accession: NC_013771, CP001842, JPSP01000003, and JPSP01000022). Absolute abundances of the 16S rRNA gene or 18S rRNA gene for each OTU were subsequently calculated based on the number of internal standard sequences recovered⁷¹. Finally, the concentrations of 16S and 18S rRNA genes in the environment were calibrated for the volume of seawater sample filtered. Common diazotrophs observed from clone library studies across

the global ocean and their eukaryotic hosts were picked out from our 16S and 18S taxonomy tables, respectively⁴²,⁷³,⁷⁹⁻⁸¹.

The internal standard method is subject to a number of limitations and caveats⁷¹. A key assumption of the approach is that the recovery efficiency of the standard is equivalent to the recovery efficiency of the natural sequences in the sample. Variation of rRNA gene copy number is also an important consideration. However, while recovery of the standard may differ from recovery of natural taxa in the same sample, variation in standard recovery efficiency from sample to sample will reflect differences in starting material, losses during elution, and other processes as long as the same PCR and library preparation protocol is followed. In that case, any biases in the quantitative measurement due to amplification biases or DNA extraction recoveries should be consistent across the samples. Therefore, the 16S rRNA approach is informative when providing the spatial distribution and abundance patterns of diazotrophic taxa.

**NCP estimates.** NCP reflects the balance between plankton community photosynthesis and respiration. An excess of photosynthesis leads to a net production of particulate and dissolved organic carbon, which can either accumulate at the ocean surface or be exported to depth. To estimate the proportion of NCP fueled by N₂ fixation, we measured NCP underway using the O₂/Ar method⁸². Oxygen concentrations in the surface ocean are influenced by biological processes, such as photosynthesis and respiration, as well as physical processes including bubble injection, temperature, and pressure changes. Due to the similar solubility properties of O₂ and Ar, the biological O₂ supersaturation ([O₂]sat) can be calculated by removing the effects of physical processes determined from Ar supersaturation ([Ar]sat)⁸³. Biological O₂ supersaturation and undersaturation reflect the metabolic state of the surface ocean, suggesting autotrophic or heterotrophic conditions, respectively⁸⁴,⁸⁵. Under steady-state conditions within the mixed layer and vertical mixing is negligible over the residence time of O₂ at the ocean surface, NCP can be estimated based on the exchange of biological O₂ with the atmosphere using the equations below.

$$ \text{NCP} \approx k_{\text{O}_2} * [\text{O}_2]_{\text{sat}} * \Delta(\text{O}_2/\text{Ar}) \quad (1) $$

$$ \Delta(\text{O}_2/\text{Ar}) = \left[ \frac{([\text{O}_2]/[\text{Ar}])}{([\text{O}_2]/[\text{Ar}])_{\text{sat}}} - 1 \right] \quad (2) $$

$k_{\text{O}_2}$ is the gas exchange velocity for oxygen⁸⁶,⁸⁷. The uncertainties in the NCP estimate are mainly from errors associated with $k_{\text{O}_2}$, and vertical mixing of O₂. Dissolved O₂/Ar ratios in surface seawater were continuously measured by Equilibrator Inlet Mass Spectrometry (EIMS) during the 2015 and 2016 cruises. O₂/Ar-NCP estimates were converted to carbon-NCP assuming a constant O₂/C stoichiometry⁸⁸,⁸⁹.

We note that our observations of NCP fueled by N₂ fixation should be interpreted with caution mainly due to differences in timescales of integration. Our O₂/Ar-NCP observations integrate productivity over 3–4 days in this region, while the N₂ fixation measurements reflect hourly or daily rates. We cannot rule out the possibility that high N₂ fixation rates occurred during late-stages of a phytoplankton bloom when nitrogen was exhausted²⁹ or, conversely that release of N by diazotrophs relieved N starvation and initiated rapid growth of non-N₂ fixers⁹⁰. This artefact in integration timescales is circumvented with dual-tracer ¹⁵N₂ and ¹³C incubations, which also show a high contribution of N₂ fixation to primary production off the coast of New Jersey (Supplementary Fig. 7). In addition, negative NCP values accompanied by detectable N₂ fixation and heterotrophic diazotrophs were observed over a large portion of the transition zone between the neritic and open ocean regions. These observations may be attributed to transient net heterotrophy, advective transport of organic matter, or vertical mixing of O₂-depleted waters. Further studies within a Lagrangian framework will be required to explore the coupling between N₂ fixation and the net metabolic status of marine systems.

The contribution of N₂ fixation to NCP measured by FARACAS and O₂/Ar method shows a similar spatial pattern as the contribution of N₂ fixation to NPP measured by ¹⁵N₂/¹³C incubation. ¹³C-based primary production measures yield rates closer to NPP than NCP⁹¹. Therefore, the ¹⁵N₂/¹³C-based approach to assessing N₂ fixation's contribution to biological production should relate to our FARACAS-O₂/Ar according to the following equation:

$$ \frac{\text{N}_2 \text{ fixation}}{\text{NCP}} * \text{export ratio} = \frac{\text{N}_2 \text{ fixation}}{\text{NPP}} \quad (3) $$

Where the export ratio = $\frac{\text{NCP}}{\text{NPP}}$. In some cases, we estimate a contribution of N₂ fixation to NCP of 80–100%. Should we assume an export ratio of 8.4% for the oligotrophic Sargasso Sea based on BATS estimates⁹², the contribution of N₂ fixation to NPP implied by our approach (6–8%) is approximately in line with our discrete incubation-based estimates (Supplementary Fig. 7). Coastal environments likely exhibit higher export ratios of 0.2–0.3⁹³, which would yield a range of contribution of N₂ fixation to NPP of up to 16–30%. Thus, while discrepancies exist between ¹⁵N₂/¹³C and FARACAS-O₂/Ar-based approaches, the broad relationship between these quantities is as expected.
**Nitrogen budget via N₂ fixation in the global ocean.** An updated database of depth-integrated and volumetric N₂ fixation rates over the global ocean is presented in Supplementary Fig. 8. The complete dataset of global N₂ fixation is shown in the Supplementary Data 1, which includes 1172 depth-integrated and 4299 volumetric N₂ fixation measurements. We conducted a Welch's *t*-test to evaluate whether the N₂ fixation rates in the coastal oceans (bathymetry ≥ −200 m) are significantly higher than in the open ocean. This one-tailed hypothesis was examined at the 0.01 significance level. N₂ fixation rates were first log-transformed since they are approximately log-normally distributed (Fig. 4a–c). Based on these analyses, the volumetric N₂ fixation rates at different depths and at surface are significantly larger in coastal regions than in the open ocean (*p* < 0.01), while depth-integrated N₂ fixation rates appear to be similar in both systems.

Nitrogen inputs through N₂ fixation were further evaluated for coastal and open oceans separately by scaling to the surface areas of the respective regions. The areal extents of the coastal and open oceans were calculated using ArcGIS. Land (orange), coastal (cyan), and open ocean (blue) regions were delineated using bathymetric contour lines (GEBCO One Minute Grid), with depth criteria of 0 m and −200 m as shown in Supplementary Fig. 9. Surface areas were calculated under a Cylindrical Equal Area Projection (World) with 10° latitudinal bands except at latitudes higher than 50° N/50° S. We calculated global budgets after removing outliers identified using a Tukey's test. These outliers include some extremely high values in the Indian Ocean⁵⁰ and some extremely low rates in the eastern North Atlantic Ocean⁹⁴ (Supplementary Fig. 8). Global budgets including outliers were also computed but may be substantially biased and are not shown here. Geometric means were used because depth-integrated and volumetric N₂ fixation rates are approximately log-normally distributed (Fig. 4a–c). Flux budgets were determined by multiplying the geometric mean of N₂ fixation rates by the area of each latitudinal band. Uncertainties were estimated based on the propagation of errors. Regional and global N₂ fixation rates are presented in Table 1 and Supplementary Table 1. For comparison, we also present budgets based on the arithmetic mean of N₂ fixation rates. Large uncertainties are expected in high latitude regions because of the limited number of observations.

**Data availability**

The updated database of depth-integrated and volumetric N₂ fixation rates over the global ocean is provided in Supplementary Data 1. In addition, data used in this study are available from the corresponding author on reasonable request.

Received: 25 July 2018 Accepted: 22 January 2019

Published online: 19 February 2019

**References**

1. Field, C. B., Behrenfeld, M. J., Randerson, J. T. & Falkowski, P. Primary production of the biosphere: integrating terrestrial and oceanic components. *Science* **281**, 237–240 (1998).

2. Moore, C. et al. Processes and patterns of oceanic nutrient limitation. *Nat. Geosci.* **6**, 701–710 (2013).

3. Falkowski, P. G. Evolution of the nitrogen cycle and its influence on the biological sequestration of CO₂ in the ocean. *Nature* **387**, 272–275 (1997).

4. Karl, D. et al. The role of nitrogen fixation in biogeochemical cycling in the subtropical North Pacific Ocean. *Nature* **388**, 533–538 (1997).

5. Capone, D. G. et al. Nitrogen fixation by *Trichodesmium* spp.: an important source of new nitrogen to the tropical and subtropical North Atlantic Ocean. *Glob. Biogeochem. Cycles* **19**, GB2024 (2005).

6. Michaels, A. F., Bates, N. R., Buesseler, K. O., Carlson, C. A. & Knap, A. H. Carbon-cycle imbalances in the Sargasso Sea. *Nature* **372**, 537–539 (1994).

7. Lipschultz, F., Bates, N. R., Carlson, C. A. & Hansell, D. A. New production in the Sargasso Sea: history and current status. *Glob. Biogeochem. Cycles* **16**, 1001 (2002).

8. Holl, C. M., Waite, A. M., Pesant, S., Thompson, P. A. & Montoya, J. P. Unicellular diazotrophy as a source of nitrogen to Leeuwin Current coastal eddies. *Deep Sea Res. Part II Top. Stud. Oceanogr.* **54**, 1045–1054 (2007).

9. Rees, A. P., Gilbert, J. A. & Kelly-Gerreyn, B. A. Nitrogen fixation in the western English Channel (NE Atlantic ocean). *Mar. Ecol. Prog. Ser.* **374**, 7–12 (2009).

10. Gruber, N. Elusive marine nitrogen fixation. *Proc. Natl Acad. Sci. USA* **113**, 4246–4248 (2016).

11. Sohm, J. A., Webb, E. A. & Capone, D. G. Emerging patterns of marine nitrogen fixation. *Nat. Rev. Microbiol.* **9**, 499–508 (2011).

12. Fernandez, C., Farias, L. & Ulloa, O. Nitrogen fixation in denitrified marine waters. *PloS One* **6**, e20539 (2011).

13. Bonnet, S. et al. Aphotic N₂ fixation in the eastern tropical South Pacific Ocean. *PloS One* **8**, e81265 (2013).

14. Bentzon-Tilia, M. et al. Significant N₂ fixation by heterotrophs, photoheterotrophs and heterocystous cyanobacteria in two temperate estuaries. *ISME J.* **9**, 273–285 (2015).

15. Needoba, J. A., Foster, R. A., Sakamoto, C., Zehr, J. P. & Johnson, K. S. Nitrogen fixation by unicellular diazotrophic cyanobacteria in the temperate oligotrophic North Pacific Ocean. *Limnol. Oceanogr.* **52**, 1317–1327 (2007).

16. Voss, M., Croot, P., Lochte, K., Mills, M., & Peeken, I. Patterns of nitrogen fixation along 10° N in the tropical Atlantic. *Geophys. Res. Lett.* **31**, L23S09 (2004).

17. Moisander, P. H. et al. Unicellular cyanobacterial distributions broaden the Oceanic N₂ fixation domain. *Science* **327**, 1512–1514 (2010).

18. Bombar, D., Paerl, R. W. & Riemann, L. Marine non-cyanobacterial diazotrophs: moving beyond molecular detection. *Trends Microbiol.* **24**, 916–927 (2016).

19. Cassar, N., Tang, W., Gabathuler, H. & Huang, K. Method for high frequency underway N₂ fixation measurements: flow-through incubation acetylene reduction assays by cavity ring down laser absorption spectroscopy (FARACAS). *Anal. Chem.* **90**, 2839–2851 (2018).

20. Orcutt, K. M. et al. A seasonal study of the significance of N₂ fixation by *Trichodesmium* spp. at the Bermuda Atlantic Time-series Study (BATS) site. *Deep Sea Res. Part II Top. Stud. Oceanogr.* **48**, 1583–1608 (2001).

21. Luo, Y. W. et al. Database of diazotrophs in global ocean: abundance, biomass and nitrogen fixation rates. *Earth Syst. Sci. Data* **4**, 47–73 (2012).

22. Mulholland, M. R. et al. Rates of dinitrogen fixation and the abundance of diazotrophs in North American coastal waters between Cape Hatteras and Georges Bank. *Limnol. Oceanogr.* **57**, 1067–1083 (2012).

23. Knapp, A. N. The sensitivity of marine N₂ fixation to dissolved inorganic nitrogen. *Front. Microbiol.* **3**, 374 (2012).

24. Sedwick, P. et al. Assessing phytoplankton nutritional status and potential impact of wet deposition in seasonally oligotrophic waters of the Mid-Atlantic Bright. *Geophys. Res. Lett.* **45**, 3203–3211 (2018).

25. Gruber, N. & Sarmiento, J. L. Global patterns of marine nitrogen fixation and denitrification. *Glob. Biogeochem. Cycles* **11**, 235–266 (1997).

26. Benavides, M. & Voss, M. Five decades of N₂ fixation research in the North Atlantic Ocean. *Front. Mar. Sci.* **2**, 40 (2015).

27. Deutsch, C., Sarmiento, J. L., Sigman, D. M., Gruber, N. & Dunne, J. P. Spatial coupling of nitrogen inputs and losses in the ocean. *Nature* **445**, 163–167 (2007).

28. Glenn, S. et al. Biogeochemical impact of summertime coastal upwelling on the New Jersey Shelf. *J. Geophys. Res.* **109**, C12S02 (2004).

29. Subramaniam, A., Mahaffey, C., Johns, W. & Mahowald, N. Equatorial upwelling enhances nitrogen fixation in the Atlantic Ocean. *Geophys. Res. Lett.* **40**, 1766–1771 (2013).

30. Fonseca-Batista, D. et al. Nitrogen fixation in the eastern Atlantic reaches similar levels in the southern and northern hemisphere. *J. Geophys. Res. Oceans* **122**, 587–601 (2017).

31. Sohm, J. A. et al. Nitrogen fixation in the South Atlantic Gyre and the Benguela upwelling system. *Geophys. Res. Lett.* **38**, L16608 (2011).

32. Subramaniam, A. et al. Amazon River enhances diazotrophy and carbon sequestration in the tropical North Atlantic Ocean. *Proc. Natl. Acad. Sci.* **105**, 10460–10465 (2008).

33. Knapp, A. N., Casciotti, K. L., Berelson, W. M., Prokopenko, M. G. & Capone, D. G. Low rates of nitrogen fixation in eastern tropical South Pacific surface waters. *Proc. Natl. Acad. Sci. USA* **113**, 4398–4403 (2016).

34. Bonnet, S., Caffin, M., Berthelot, H. & Moutin, T. Hot spot of N₂ fixation in the western tropical South Pacific pleads for a spatial decoupling between N₂ fixation and denitrification. *Proc. Natl. Acad. Sci. USA* **114**, 2800–2801 (2017).

35. Tagliabue, A. et al. The integral role of iron in ocean biogeochemistry. *Nature* **543**, 51–59 (2017).

36. Parekh, P., Follows, M. J. & Boyle, E. A. Decoupling of iron and phosphate in the global ocean. *Glob. Biogeochem. Cycles* **19**, GB2020 (2005).

37. Fennel, K. The role of continental shelves in nitrogen and carbon cycling: Northwestern North Atlantic case study. *Ocean Sci.* **6**, 539–548 (2010).

38. Conway, T. M. & John, S. G. Quantification of dissolved iron sources to the North Atlantic Ocean. *Nature* **511**, 212–215 (2014).

39. Lee, K., Karl, D. M., Wanninkhof, R. & Zhang, J. Z. Global estimates of net carbon production in the nitrate-depleted tropical and subtropical oceans. *Geophys. Res. Lett.* **29**, 1907 (2002).

40. Rahav, E., Giannetto, M. & Bar-Zeev, E. Contribution of mono and polysaccharides to heterotrophic N₂ fixation at the eastern Mediterranean coastline. *Sci. Rep.* **6**, 27858 (2016).

41. Gaby, J. C. & Buckley, D. H. A comprehensive aligned aligned *nifH* gene database: a multipurpose tool for studies of nitrogen-fixing bacteria. *Database* **2014**, bau001 (2014).

42. Turk-Kubo, K. A., Karamchandani, M., Capone, D. G. & Zehr, J. P. The paradox of marine heterotrophic nitrogen fixation: abundances of heterotrophic diazotrophs do not account for nitrogen fixation rates in the Eastern Tropical South Pacific. *Environ. Microbiol.* **16**, 3095–3114 (2014).

43. Thompson, A. W. et al. Unicellular cyanobacterium symbiotic with a single-celled eukaryotic alga. *Science* **337**, 1546–1550 (2012).

44. Cabello, A. M. et al. Global distribution and vertical patterns of a prymnesiophyte-cyanobacteria obligate symbiosis. *ISME J.* **10**, 693–706 (2015).
45. Zehr, J. P., Shilova, I. N., Farnelid, H. M., Muñoz-Marín, M. D. C. & Turk-Kubo, K. A. Unusual marine unicellular symbiosis with the nitrogen-fixing cyanobacterium UCYN-A. *Nat. Microbiol.* **2**, 16214 (2016).

46. Breitbarth, E., Oschlies, A. & LaRoche, J. Physiological constraints on the global distribution of *Trichodesmium*—effects of temperature on diazotrophy. *Biogeosciences* **4**, 53–61 (2008).

47. Mourino-Carballido, B. et al. Importance of N₂ fixation vs. nitrate eddy diffusion along a latitudinal transect in the Atlantic Ocean. *Limnol. Oceanogr.* **56**, 999–1007 (2011).

48. Benavides, M. et al. Low contribution of N₂ fixation to new production and excess nitrogen in the subtropical northeast Atlantic margin. *Deep Sea Res. Part I Oceanogr. Res. Pap.* **81**, 36–48 (2013).

49. Stanley, R. H., Jenkins, W. J., Doney, S. C., & Lott, D. III. The ³He flux gauge in the Sargasso Sea: a determination of physical nutrient fluxes to the euphotic zone at the Bermuda Atlantic Time-series Site. *Biogeosciences* **12**, 5199–5210 (2015).

50. Gandhi, N. et al. First direct measurements of N₂ fixation during a *Trichodesmium* bloom in the eastern Arabian Sea. *Glob. Biogeochem. Cycles* **25**, GB4014 (2011).

51. Großkopf, T. et al. Doubling of marine dinitrogen-fixation rates based on direct measurements. *Nature* **488**, 361–364 (2012).

52. Capone, D. G. in *Nitrogen Cycling in Coastal Marine Environments* (eds Blackburn, T. H. & Sørensen, J.) 85–123 (Wiley, Chichester, 1988).

53. Mohr, W., Großkopf, T., Wallace, D. W. & LaRoche, J. Methodological underestimation of oceanic nitrogen fixation rates. *PloS One* **5**, e12583 (2010).

54. Dabundo, R. et al. The contamination of commercial ¹⁵N₂ gas stocks with ¹⁵N-labeled nitrate and ammonium and consequences for nitrogen fixation measurements. *PloS One* **9**, e110335 (2014).

55. Bauer, J. E. et al. The changing carbon cycle of the coastal ocean. *Nature* **504**, 61–70 (2013).

56. Doney, S. C. The growing human footprint on coastal and open-ocean biogeochemistry. *Science* **328**, 1512–1516 (2010).

57. Hardy, R. W., Holsten, R., Jackson, E. & Burns, R. The acetylene-ethylene assay for N₂ fixation: laboratory and field evaluation. *Plant Physiol.* **43**, 1185–1207 (1968).

58. Capone, D. G. & Montoya, J. P. Nitrogen fixation and denitrification. *Methods Microbiol.* **30**, 501–515 (2001).

59. Wilson, S. T., Bottjer, D., Church, M. J. & Karl, D. M. Comparative assessment of nitrogen fixation methodologies, conducted in the oligotrophic North Pacific Ocean. *Appl. Environ. Microbiol.* **78**, 6516–6523 (2012).

60. Capone, D. G. & Carpenter, E. J. Perfusion method for assaying microbial activities in sediments: applicability to studies of N₂ fixation by C₂H₂ reduction. *Appl. Environ. Microbiol.* **43**, 1400–1405 (1982).

61. Fulweiler, R. W. et al. Examining the impact of acetylene on N-fixation and the active sediment microbial community. *Front. Microbiol.* **6**, 418 (2015).

62. Payne, W. Influence of acetylene on microbial and enzymatic assays. *J. Microbiol. Methods* **2**, 117–133 (1984).

63. Cassar, N., Bellenger, J. P., Jackson, R. B., Karr, J. & Barnett, B. A. N₂ fixation estimates in real-time by cavity ring-down laser absorption spectroscopy. *Oecologia* **168**, 335–342 (2012).

64. Montoya, J. P., Voss, M., Kahler, P. & Capone, D. G. A simple, high-precision, high-sensitivity tracer assay for N₂ fixation. *Appl. Environ. Microbiol.* **62**, 986–993 (1996).

65. Staal, M. et al. Nitrogen fixation along a north-south transect in the eastern Atlantic Ocean. *Limnol. Oceanogr.* **52**, 1305–1316 (2007).

66. Mulholland, M. R., Bronk, D. A. & Capone, D. G. Dinitrogen fixation and release of ammonium and dissolved organic nitrogen by *Trichodesmium* IMS101. *Aquat. Microb. Ecol.* **37**, 85–94 (2004).

67. Bruland, K. W., Rue, E. L., Smith, G. J. & DiTullio, G. R. Iron, macronutrients and diatom blooms in the Peru upwelling regime: brown and blue waters of Peru. *Mar. Chem.* **93**, 81–103 (2005).

68. Lagerström, M. et al. Automated on-line flow-injection ICP-MS determination of trace metals (Mn, Fe, Co, Ni, Cu and Zn) in open ocean seawater: application to the GEOTRACES program. *Mar. Chem.* **155**, 71–80 (2013).

69. Martínez-Pérez, C. et al. The small unicellular diazotrophic symbiont, UCYN-A, is a key player in the marine nitrogen cycle. *Nat. Microbiol.* **1**, 16163 (2016).

70. Wang, S., Lin, Y., Gifford, S., Eveleth, R. & Cassar, N. Linking patterns of net community production and marine microbial community structure in the western North Atlantic. *ISME J.* **12**, 2582–2595 (2018).

71. Satinsky, B. M., Gifford, S. M., Crump, B. C. & Moran, M. A. Use of internal standards for quantitative metatranscriptome and metagenome analysis. *Methods Enzymol.* **531**, 237–250 (2013).

72. Smets, W. et al. A method for simultaneous measurement of soil bacterial abundances and community composition via 16S rRNA gene sequencing. *Soil Biol. Biochem.* **96**, 145–151 (2016).

73. Moisander, P. H., Beinart, R. A., Voss, M. & Zehr, J. P. Diversity and abundance of diazotrophic microorganisms in the South China Sea during intermonsoon. *ISME J.* **2**, 954–967 (2008).

74. Caporaso, J. G. et al. QIIME allows analysis of high-throughput community sequencing data. *Nat. Methods* **7**, 335–336 (2010).

75. Fadrosh, D. W. et al. An improved dual-indexing approach for multiplexed 16S rRNA gene sequencing on the Illumina MiSeq platform. *Microbiome* **2**, 6 (2014).

76. Pruesse, E. et al. SILVA: a comprehensive online resource for quality checked and aligned ribosomal RNA sequence data compatible with ARB. *Nucleic Acids Res.* **35**, 7188–7196 (2007).

77. Edgar, R. C. Search and clustering orders of magnitude faster than BLAST. *Bioinformatics* **26**, 2460–2461 (2010).

78. Edgar, R. C., Haas, B. J., Clemente, J. C., Quince, C. & Knight, R. UCHIME improves sensitivity and speed of chimera detection. *Bioinformatics* **27**, 2194–2200 (2011).

79. Shiozaki, T., Ijichi, M., Kodama, T., Takeda, S. & Furuya, K. Heterotrophic bacteria as major nitrogen fixers in the euphotic zone of the Indian Ocean. *Glob. Biogeochem. Cycles* **28**, 1096–1110 (2014).

80. Halm, H. et al. Heterotrophic organisms dominate nitrogen fixation in the South Pacific Gyre. *ISME J.* **6**, 1238–1249 (2012).

81. Zehr, J. P. Nitrogen fixation by marine cyanobacteria. *Trends Microbiol.* **19**, 162–173 (2011).

82. Cassar, N. et al. Continuous high-frequency dissolved O₂/Ar measurements by equilibrator inlet mass spectrometry. *Anal. Chem.* **81**, 1855–1864 (2009).

83. Craig, H. & Hayward, T. Oxygen supersaturation in the ocean: biological versus physical contributions. *Science* **235**, 199–202 (1987).

84. Emerson, S., Quay, P. & Wheeler, P. A. Biological productivity determined from oxygen mass balance and incubation experiments. *Deep Sea Res. Part I Oceanogr. Res. Pap.* **40**, 2351–2358 (1993).

85. Spitzer, W. S. & Jenkins, W. J. Rates of vertical mixing, gas exchange and new production: Estimates from seasonal gas cycles in the upper ocean near Bermuda. *J. Mar. Res.* **47**, 169–196 (1989).

86. Wanninkhof, R. Relationship between wind speed and gas exchange over the ocean. *J. Geophys. Res.* **97**, 7373–7382 (1992).

87. Reuer, M. K., Barnett, B. A., Bender, M. L., Falkowski, P. G. & Hendricks, M. B. New estimates of Southern Ocean biological production rates from O₂/Ar ratios and the triple isotope composition of O₂. *Deep Sea Res. Part I Oceanogr. Res. Pap.* **54**, 951–974 (2007).

88. Shulenberger, E. & Reid, J. L. The Pacific shallow oxygen maximum, deep chlorophyll maximum, and primary productivity, reconsidered. *Deep Sea Res. Part A Oceanogr. Res. Pap.* **28**, 901–919 (1981).

89. Laws, E. A. Photosynthetic quotients, new production and net community production in the open ocean. *Deep Sea Res. Part A Oceanogr. Res. Pap.* **38**, 143–167 (1991).

90. Berthelot, H., Bonnet, S., Grosso, O., Cornet, V. & Barani, A. Transfer of diazotroph-derived nitrogen towards non-diazotrophic planktonic communities: a comparative study between *Trichodesmium erythraeum*, *Crocosphaera watsonii* and *Cyanothoe* sp. *Biogeosciences* **13**, 4005–4021 (2016).

91. Regaudie-de-Gioux, A., Lasternas, S., Agusti, S. & Duarte, C. M. Comparing marine primary production estimates through different methods and development of conversion equations. *Front. Mar. Sci.* **1**, 19 (2014).

92. Siegel, D. A. et al. Global assessment of ocean carbon export by combining satellite observations and food-web models. *Glob. Biogeochem. Cycles* **28**, 181–196 (2014).

93. Laws, E. A., D'Sa, E. & Naik, P. Simple equations to estimate ratios of new or export production to total production from satellite-derived estimates of sea surface temperature and primary production. *Limnol. Oceanogr. Methods* **9**, 593–601 (2011).

94. Benavides, M., Agawin, N. S., Arístegui, J., Ferriol, P. & Stal, L. J. Nitrogen fixation by *Trichodesmium* and small diazotrophs in the subtropical northeast Atlantic. *Aquat. Microb. Ecol.* **65**, 43–53 (2011).

## Acknowledgements

We would like to thank the captains, crew, and marine technicians of the R/V *Atlantic Explorer*. We further thank Geoffrey Smith (UCSC) for his help with deployment of the towhish system. We also thank Zuchuan Li (Duke University, now at WHOI) for providing satellite images for adaptive sampling during the cruises. We are grateful to Yawei Luo (Xiamen University), Run Zhang (Xiamen University), Lasse Riemann (University of Copenhagen), Eyal Rahav (Israel Oceanographic and Limnological Research), Ilana Berman-Frank (Bar Ilan University), Carolin Loscher (Helmholtz Center for Ocean Research Kiel), Camila Fernandez (CNRS), Takuhei Shiozaki (Japan Agency for Marine-Earth Science and Technology), and Ajit Subramaniam (Lamont-Doherty Earth Observatory) for providing access to their datasets. We would like to acknowledge NASA for processing and distributing satellite data on OceanColor Web. This work was funded by an NSF-CAREER grant (#1350710) to N.C. and a Link Foundation Ocean Engineering & Instrumentation Fellowship to W. T. N.C. was also in part supported by the "Laboratoire
d'Excellence" LabexMER (ANR-10-LABX-19) and co-funded by a grant from the French government under the program "Investissements d'Avenir". H.P. was funded by RPDOC BITMAP (ANR-12-PDOC-0025-01). D.F.-B. and F.D. were supported by the VUB R&D, Strategic Research Plan "Tracers of Past & Present Global Changes".

**Author contributions**

N.C. designed the study. W.T. and N.C. performed the measurements and analysis of N₂ fixation and wrote the manuscript with contributions from all co-authors. S.W. performed and analyzed the NCP measurements. S.W., W.T., and S.G. performed molecular sequencing data collection and analyses. D.F.-B. and F.D. performed the ¹⁵N₂ and ¹³C incubation experiments and analyses. A.G.G, H.P., G.S., and M.G. collected and conducted the trace-metal analyses.

**Additional information**

**Supplementary Information** accompanies this paper at https://doi.org/10.1038/s41467-019-08640-0.

**Competing interests:** The authors declare no competing interests.

**Reprints and permission** information is available online at http://npg.nature.com/reprintsandpermissions/

**Journal peer review information:** *Nature Communications* thanks Douglas Capone and the other anonymous reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_6.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2019