Samples for the quantification of nitrate (NO3-) and dissolved inorganic phosphorus (DIP) concentrations were collected at 12 depths between 0 and 200 m in acid-washed polyethylene bottles, poisoned with mercuric chloride (HgCl2, final concentration 20 mg L-1) and stored at 4 ∘C until analysis. Concentrations were determined using standard colorimetric techniques on a Bran Luebbe AA3 autoanalyzer. Detection limits for the procedures were 0.05 µmolL-1 for NO3- and DIP.

Samples for determining dissolved Fe concentrations were collected and analyzed as described in . Briefly, samples were collected using a Titane rosette mounted with 24 Teflon-coated 12 L GoFlos deployed with a Kevlar cable. Dissolved Fe concentrations were measured by flow injection with online preconcentration and chemiluminescence detection according to . The reliability of the method was monitored by analyzing the D1 SAFe seawater standard , and an internal acidified seawater standard was measured daily to monitor the stability of the analysis.

The sampling and analytical methods used to analyze the other parameters reported in the correlation table (Table 2) are described in detail in the methods sections for related papers in this issue .

Whole water (bulk) N2 fixation rates were measured in triplicate at all 17 stations using the 15N2 isotopic tracer technique (adapted from ). The 15N2 bubble technique was intentionally chosen due to the time limitation on making enriched 15N2 seawater inoculates (e.g., 6–9 depths = 6–9 inoculates) and the larger sample bottles required for making proper estimates of activity in oligotrophic environments. In addition, we aimed to avoid any potential overestimation due to trace metal and dissolved organic matter (DOM) contaminations often associated with the preparation of the 15N2-enriched seawater in our incubation bottles as Fe and DOM have been found to control N2 fixation or nifH gene expression in this region . However, the 15N / 14N ratio of the N2 pool available for N2 fixation (the term AN2 used in ) was measured in all incubation bottles by membrane inlet mass spectrometry (MIMS) to ensure accurate rate calculations (see below).

Seawater samples were collected from Niskin bottles into 10 % HCl-washed, sample-rinsed (three times) light-transparent polycarbonate (2.3 L) bottles from six depths (75, 50, 20, 10, 1, and 0.1 % surface irradiance levels) at all short-duration stations SD1 to SD15 and nine depths (75, 50, 35, 20, 10, 3, 1, 0.3, and 0.1 % surface irradiance levels) at LD A, LDB and LD C, corresponding to the sub-surface (5 m) down to 80 to 180 m, depending on the station. Bottles were sealed with caps fitted with silicon septa and amended with 2 mL of 98.9 at. % 15N2 (Cambridge isotopes). The purity of the 15N2 Cambridge isotope stocks was previously checked by and more recently by and . They were found to be lower than 2×10-8 mol : mol of 15N2, leading to a potential N2 fixation rate overestimation of < 1 %. Each bottle was shaken 20 times to break the 15N2 bubble and facilitate its dissolution, and was incubated for 24 h. At SD stations, bottles were incubated in on-deck incubators connected to surface circulating seawater at the specified irradiances using blue screening as the duration of the station (8 h) was too short to deploy in situ mooring lines. At LD stations (7 days), one profile was incubated following the same methodology in on-deck incubators and another replicate profile was incubated in situ for comparison on a drifting mooring line located at the same depth from which the samples were collected. Incubations were stopped by filtering the entire incubation bottle onto pre-combusted (450 ∘C, 4 h) 25 mm diameter glass fiber filters (GF/F, Whatman, 0.7 µm nominal pore size). Filters were subsequently dried at 60 ∘C for 24 h before analysis of 15N / 14N ratios and particulate N (PN) determinations using an elemental analyzer coupled to a mass spectrometer (EA-IRMS, Integra CN, SerCon Ltd) as described in .

To ensure accurate rate calculations, the 15N / 14N ratio of the N2 pool in the incubation bottles was measured on each profile from triplicate surface incubation bottles from SD1 to SD14 and at all depths at SD15 and LD stations. Briefly, 12 mL was subsampled after incubation into Exetainers fixed with HgCl2 (final concentration 20 mg L-1) that were preserved upside down in the dark at 4 ∘C until analyzed using a MIMS according to . Lastly, we collected time zero samples at each station to determine the natural N isotopic signature of ambient particulate N (PN). The minimum quantifiable rates (quantification limit, QL) calculated using standard propagation of errors via the observed variability between replicate samples measured according to were 0.035 nmol N L-1 d-1.

Discrete N2 fixation rate measurements were depth integrated over the photic layer using trapezoidal integration procedures. Briefly, the N2 fixation at each pair of depths is averaged, and then multiplied by the difference between the two depths to get a total N2 fixation in that depth interval. These depth interval values are then summed over the entire depth range to get the integrated N2 fixation rate. The rate nearest the surface is assumed to be constant up to 0 m .

Spearman's rank correlation was used to examine the potential relationships between N2 fixation rates and hydrological, biogeochemical, and biological parameters across the longitudinal transect (n=102, α=0.05). A non-parametric Mann–Whitney test (α=0.05) was used to compare the MIMS data obtained following on-deck versus in situ incubations, and to compare nutrient and Chl a distributions between the western part and the eastern part of the transect.

Seawater temperature ranged from 21.4 to 30.0 ∘C in the sampled photic layer (0 to ∼ 80–180 m) over the cruise transect (Fig. 2a). The mixed layer depth (MLD) calculated according to the method was located around 20–40 m throughout the zonal transect: maximum temperatures were measured in the surface mixed layer (∼ 0–20/40 m) and remained almost constant along the longitudinal transect with 29.1 ± 0.3 ∘C in MA waters and 29.5 ± 0.4 ∘C in GY waters.

Based on other hydrographic measurements (dissolved nutrients, dissolved Fe and Chl a concentrations), the longitudinal transect was divided into two sub-regions: (1) the MA region from station SD1 (160∘ E) to LDB (165∘ W), and (2) the GY sub-region from station LDB (165 ∘ W) to SD15 (160∘ W). Chl a concentration in the upper 50 m was significantly (p< 0.05) higher in MA waters (0.17 µgL-1 on average) than in GY waters (0.06 µgL-1 on average) (Figs. 1 and 2b). The DCM was located around 80–100 m in MA waters and deepened to ∼ 150 m in GY waters, indicating higher oligotrophy in the GY region. Surface NO3- concentrations (Fig. 2c) were consistently close to or below the detection limit (0.05 µmolL-1) in the upper water column (0–50 m) throughout the transect and the depth of the nitracline gradually deepened from ∼ 75–100 m in MA waters to ∼ 115 m in GY waters. DIP concentrations were slightly higher than or close to detection limits (0.05 µmolL-1) in MA surface waters (0–50 m), and the phosphacline was shallower (20–45 m) than the nitracline and DIP concentrations increased significantly (p < 0.05) in GY waters and ranged from 0.13 to 0.17 µmolL-1 (Fig. 2d).

The 15N enrichment of the N2 pool after 24 h of incubation with the 15N2 tracer was on average 6.145 ± 0.798 at. % (n=54) in bottles incubated in on-deck incubators and significantly higher (p < 0.05) in bottles incubated on the mooring line (7.548 ± 0.557 at. % (n=44), Fig. 3a). However, the depth of incubation on the mooring line (between 5 and 180 m) did not have any significant effect (p > 0.05) on the isotopic signature of the N2 pool at LDB and LDC, which remained constant over the water column (Fig. 3b).

The natural N isotopic signature of suspended particles measured over the photic layer was on average -0.41 ‰ in MA waters and 8.06 ‰ in GY waters (Table 1). Those isotopic values were used as time zero samples to calculate N2 fixation rates.

N2 fixation was detected at all 17 sampled stations and the average measured N2 fixation rates in the two previously defined sub-regions were (1) 8.9 ± 10 nmol N L-1 d-1 (range QL-48 nmol N L-1 d-1) over the photic layer in MA waters, and 2) 0.5 ± 0.6 nmol N L-1 d-1 (range QL-4.0 nmol N L-1 d-1) in GY waters (Fig. 2e). In MA waters, N2 fixation was largely restricted to the mixed layer, where average rates were 15 nmol N L-1 d-1, with local maxima (> 20 nmol N L-1 d-1) at stations SD1, SD6 and LDB and local minima (< 5 nmol N L-1 d-1) at SD8 and SD10. In GY waters, maximum rates reached 1–2 nmol N L-1 d-1 and were located deeper in the water column (∼ 50 m). When integrated over the photic layer, N2 fixation represented an average net N addition of 631 ± 286 µmolNm-2d-1 (range 196–1153 µmolNm-2d-1) in MA waters and of 85 ± 79 µmolNm-2d-1 (range 18–172 µmolNm-2d-1) in GY waters.

N2 fixation rates were significantly positively correlated with seawater temperature and photosynthetically active radiation (PAR) (p < 0.05), significantly negatively correlated with depth and salinity (p < 0.05) and not significantly correlated with dissolved oxygen concentrations (p > 0.05) (Table 2). N2 fixation rates were significantly positively correlated with dissolved Fe, dissolved organic N (DON), phosphorus (DOP), carbon (DOC), particulate organic N (PON), particulate organic carbon (POC), biogenic silica (BSi), Chl a concentrations, and primary and bacterial production (p < 0.05), and significantly negatively correlated with concentration of dissolved nutrients: NO3-, NH4+, DIP and silicate (p < 0.05). N2 fixation rates were significantly positively correlated with nifH abundances for Trichodesmium spp., UCYN-B and the symbionts of DDAs (het-1, het-2) (p < 0.05) and not significantly correlated with UCYN-A1 and UCYN-A2 abundances (p > 0.05) (Table 2). Regarding non-diazotrophic plankton, N2 fixation rates were significantly positively correlated with Prochlorococcus spp., Synechococcus spp., heterotrophic bacteria and protist abundances (p < 0.05) and significantly negatively correlated with picoeukaryotes (p < 0.05).

At the three stations where cell-specific N2 fixation rates were estimated by nanoSIMS (SD2, SD6 and LDB), the most abundant diazotroph phylotype was Trichodesmium with 1.3 × 105, 3.3 × 105 and 1.2 × 105 cells L-1, respectively, followed by UCYN-B, whose abundances were 2.0 × 104, 1.5 × 105 and 3.8 × 102 nifH copies L-1, respectively. Het-1 and het-2 combined were 1 to 2 orders of magnitude lower, ranging from 1.0 to 9.9 × 103 nifH copies L-1, and UCYN-A1 were below detection at the three stations. In summary, Trichodesmium and UCYN-B accounted for 98.2, 99.8 and 92.1 % of the total diazotroph community (based on the phylotypes targeted here) at SD2, SD6 and LDB, respectively (Table 3).

The 15N / 14N ratios of individual cells/trichomes of UCYN-B and Trichodesmium were measured via nanoSIMS analyses and used to estimate single-cell N2 fixation rates. A summary of the enrichment values and cell-specific N2 fixation is provided in Table 3. Individual trichomes exhibited significant 15N enrichments (0.610 ±  0.269, 0.637 ± 0.355 and 0.981 ± 0.466 at. % at stations SD2, SD6 and LDB, respectively) compared with time zero samples (0.369 ± 0.002 at. %). UCYN-B were also significantly 15N-enriched, with 1.163 ± 0.531 and 0.517 ± 0.237 at. % at SD2 and SD6, respectively (note that no UCYN-B could be sorted and analyzed by nanoSIMS at LDB as they accounted for only 0.3 % of the diazotroph community). Cell-specific N2 fixation rates of Trichodesmium were 38.9 ± 8.1, 29.3 ± 5.4 and 123.8 ± 24.8 fmol N cell-1 d-1 at SD2, SD6 and LDB. Cell-specific N2 fixation of UCYN-B was 30.0 ± 6.4 and 6.1 ± 1.2 fmol N cell-1 d-1 at SD1 and SD6. The contribution of Trichodesmium to bulk N2 fixation was 83.8, 47.1 and 52.9 % at stations SD2, SD6 and LDB, respectively. The contribution of UCYN-B was 10.1 and 6.1 % at SD2 and SD6, respectively (Table 3).

The in situ nifH expression for all diazotroph groups targeted by qPCR was estimated using a TaqMAN quantitative reverse transcription PCR (RT-qPCR) (Table 4). The sampling and filtering time (17:00–21:00 h) was not optimal for quantifying the nifH gene expression for all diazotrophs; however, it provides useful information about which diazotrophs were potentially active during the experiment and complements the nanoSIMS analysis which measures the in situ activity. Both Trichodesmium and UCYN-B dominated the biomass , as did their nifH gene expression at all three stations, especially SD2 and SD6. Of the two DDAs, het-1 had a higher nifH gene expression, which was consistent with its higher nifH abundance by DNA qPCR . UCYN-A1 was consistently below detection for the nifH gene expression and was also the least detected diazotroph by nifH qPCR .

Our understanding of the marine N cycle relies on accurate estimates of N fluxes to and from the ocean. Here we decided to use the “bubble addition method” to minimize potential trace metal and organic matter contaminations, which may have resulted in overestimating rates . Moreover, a recent extensive meta-analysis (13 studies, 368 observations) between bubble and enriched amendment experiments to measure 15N2 rates reported that underestimation of N2 fixation is negligible in experiments that last 12–24 h (e.g., error is -0.2 %); hence our 24 h based experiments should be within a small amount of error . However, we paid careful attention to accurately measure the term AN2 to avoid any potential underestimation and reveal that the way bottles are incubated (on-deck versus in situ) has a great influence on the AN2 value, and thus on N2 fixation rates.

Our MIMS results measured a significantly (p < 0.05) lower 15N enrichment of the N2 pool (6.145 ± 0.798 at. %) when bottles were incubated in on-deck incubators compared to when bottles were incubated on the mooring line (7.548 ± 0.557 at. %). This suggests that the 15N2 dissolution is much more efficient when bottles are incubated in situ, probably due to the higher pressure in seawater at the depth of incubation (1.5 to 19 bars between 5 and 180 m) compared to the pressure in the on-deck incubators (1 bar). The seawater temperature checked regularly in the on-deck incubators was equivalent to ambient surface temperature and likely did not explain the differences observed. This result highlights the need to perform routine MIMS measurements to use the most accurate AN2 value for rate calculations, independently of the 15N2 approach used (gas or dissolved). In our study, the theoretical AN2 value based on gas constant calculations was ∼ 8.2 at. %, so the deviation from this value is more important when bottles are incubated in on-deck incubators as compared to when they are incubated in situ. This suggests that the use of the bubble addition method without MIMS measurement potentially leads to higher underestimations when bottles are incubated in on-deck incubators, which is the case in the great majority of marine N2 fixation studies published so far . We are aware that the dissolution kinetics of 15N2 in the incubation bottles may have been progressive along the 24 h of incubation ; therefore, the N2 fixation rates provided here represent conservative values.

Despite the AN2 value being different according to the incubation mode, it did not change with the depth of incubation on the mooring line, indicating that a slightly higher pressure than atmospheric pressure (1.5 bar at 5 m depth) is enough to promote the 15N2 dissolution. It also indicates that the slightly lower seawater temperature (22–24 ∘C) recorded at ∼ 100–180 m where the deepest samples were incubated likely did not affect the solubilization of the 15N2 gas. In our study, the vertical profiles performed at LD stations and incubated either on-deck in triplicate or in situ in triplicate reveal identical (p > 0.05) N2 fixation rates regardless of the incubation method used . This indicates that in situ incubations and on-deck incubations that simulate appropriate light levels are a valid methodology for 15N2 fixation rate measurements on cruises during which mooring lines cannot be deployed, as long as routine measurements of the isotopic ratio of the N2 pool are performed in incubation bottles.

N2 fixation rates measured in MA waters (average 631 ± 286 µmolNm-2d-1) are 3 to 4 times higher than model predictions for this area (150–200 µmolNm-2d-1, ). They are in the upper range of the higher category (100–1000 µmolNm-2d-1) of rates defined by in the N2 fixation MAREDAT database for the global ocean and thus identify the WTSP as an area for high N2 fixation in the global ocean. Recent studies performed in the western part of the WTSP, i.e., in the Solomon, Bismarck and Arafura seas, also reveal extremely high rates (> 600 µmolNm-2d-1), indicating that this high N2 fixation activity area extends geographically west–east from Australia to Tonga and north–south from the Equator to 25–30∘ S, covering a vast ocean area of ∼ 13×106 km2 (i.e., ∼ 20 % of the South Pacific Ocean area). However, the driver(s) for diazotrophy in this region is(are) still poorly resolved and raise(s) the question of which factors influence the distribution and activity of N2 fixation in the ocean. In a global-scale study conducted by , which investigated the correlations between N2 fixation and a variety of environmental parameters commonly accepted to control this process, they concluded that SST (or surface solar radiation) was the best predictor to explain the spatial distribution of N2 fixation in the surface ocean. Below we highlight the most plausible factors explaining such high N2 fixation rates in our study area.

Seawater temperature. Seawater temperature was unlikely to be the factor explaining the differences in N2 fixation rates observed between MA and GY waters, as it was consistently high (> 28 ∘C in the surface mixed layer) and optimal for the growth and nitrogenase activity of most diazotrophs all along the cruise transect. This indicates that other factors such as nutrient availability may explain the distribution of N2 fixation.

DIP availability. The ∼ 4000 km transect was divided into two main sub-regions: (1) the MA waters, harboring typical oligotrophic conditions with surface NO3- and DIP concentrations close to detection limits (0.05 µmolL-1), a nitracline located at 75–100 m, moderate surface Chl a concentrations (∼ 0.17 µgL-1), a DCM located at ∼ 80–100 m and very high N2 fixation rates (631 ± 286 µmolNm-2d-1 on average), and (2) the GY waters harboring ultra-oligotrophic conditions with undetectable NO3-, a deeper nitracline (115 m), and comparatively high DIP concentrations (0.15 µmolL-1), very low Chl a concentrations (0.06 µgL-1, DCM ∼ 150 m) and low N2 fixation rates (85 ± 79 µmolNm-2d-1).

In the NO3--depleted MA waters, low DIP concentrations are indicative of the consumption of DIP by the planktonic community, including diazotrophs. This is consistent with the negative correlation found between N2 fixation and DIP turnover time (the ratio between DIP concentrations and DIP uptake rates) (Table 2) and suggests a higher DIP limitation when N2 fixation is high and consumes DIP. The high DIP concentrations (> 0.1 µmolL-1) in GY surface waters compared to MA waters are consistent with former studies that consider the South Pacific Gyre as a high phosphate, low chlorophyll ecosystem , in which DIP accumulates in the absence of NO3- and low N2 fixation activity. In the high phosphate, low chlorophyll scenario, the community is limited by temperature and/or Fe availability . During the OUTPACE cruise, the DIP turnover time was variable but close to or below 2 days in MA waters , indicating a potential limitation by DIP at some stations. Trichodesmium, the most abundant and major contributor to N2 fixation during the cruise, is known to synthesize hydrolytic enzymes in order to acquire P from the dissolved organic phosphorus pool (DOP) . Moreover, Trichodesmium spp. differs from the other major diazotroph UCYN-B enumerated on the cruise in the forms of organic P it can synthesize. It is thus likely that DOP species that favored Trichodesmium over UCYN-B played a role in maintaining high Trichodesmium biomass in MA waters. It has to be noted that average DIP turnover times in MA waters were always much higher than those typically measured in severely DIP-limited environments such as the Mediterranean and Sargasso seas (e.g., ), suggesting that DIP concentrations are generally favorable for the development of certain diazotrophs in the WTSP, and do not alone explain why N2 fixation is high in MA waters and low in GY waters. However, it is likely that the depletion of DIP stocks at the end of the austral summer season forces the decline of diazotrophic blooms in the WTSP , concomitantly with the decline of SST.

Fe availability. Before OUTPACE, our knowledge on Fe sources and concentrations in the WTSP was limited, especially in MA waters. During OUTPACE, reported high dissolved Fe concentrations in MA waters (range 0.2–66.2 nmol L-1, 1.7 nmol L-1 on average over the photic layer), i.e., significantly (p < 0.05) higher than those reported in GY waters (range 0.2–0.6 nmol L-1, 0.3 nM on average over the photic layer). The low dissolved Fe concentrations measured in the GY waters agree well with previous reports for the same region . However, the high dissolved Fe concentrations measured in MA waters were previously undocumented and reveal several maxima (> 50 nmol L-1) between stations SD7 and SD11, indicative of intense fertilization processes taking place in this region. found that atmospheric deposition measured during the cruise in this region was too low to explain the observed dissolved Fe concentrations in the surface water column. The seafloor of the WTSP hosts the Tonga–Kermadec subduction zone which stretches 2500 km from New Zealand to the Tonga archipelago (Fig. 1). It has among the highest densities of submarine volcanoes associated with hydrothermal vents recorded in the ocean (2.6 vents / 100 km, ), which discharge large quantities of material into the water column, including biogeochemically relevant elements such as Fe and Mn. used hydrological data recorded by Argo float in situ measurements, atlas data and simulations from a general ocean circulation model to argue that the high dissolved Fe concentrations may be sustained by a submarine source. They show that such Fe inputs could spread throughout the WTSP through mesoscale activity mainly westward through the South Equatorial Current, SEC (Fig. 1). hypothesize that the high dissolved Fe concentrations in MA waters compared to the GY ones is due to shallow inputs of hydrothermal origin together with potential Fe input from islands themselves . In our study, dissolved Fe concentrations were significantly positively correlated with N2 fixation and help to explain the distribution of N2 fixation rates measured across the OUTPACE transect.

In summary, our hypothesis to explain the spatial distribution of N2 fixation in MA waters is the following: when high DIP waters flow westward from the ETSP through the SEC and cross the South Pacific Gyre, N2-fixing organisms do not develop despite optimal SST (> 25 ∘C), likely because GY waters are Fe-depleted . When the high DIP, low DIN (dissolved inorganic N) waters from the gyre are advected west of the Tonga trench in Fe-rich and warm (> 25 ∘C) waters, all environmental conditions are fulfilled for diazotrophs to bloom extensively. According to , the strong depth difference between the nitracline and the phosphacline in MA waters associated with winter mixing allows a seasonal replenishment of DIP, which creates an excess of P relative to N and thus also favors N2 fixation in this region. Further investigations are required to better quantify Fe input from both islands and shallow volcanoes and associated hydrothermal activity along the Tonga volcanic arc for the upper mixed layer, study the fate of hydrothermal plumes in the water column at the local and regional scales, and investigate the potential impact of such hydrothermal inputs on diazotrophic communities at the scale of the whole WTSP.

N2 fixation rates were significantly negatively correlated with NO3- concentrations, consistent with the high energetic cost of N2 fixation compared to NO3- assimilation (Falkowski, 1983). They were also negatively correlated with depth and logically positively correlated with PAR and seawater temperature, two parameters which are depth dependent. Most N2 fixation took place in the surface mixed layer and rates were ∼ 15 nmol N L-1 d-1 in MA waters with local maxima at stations SD1, SD6 and LDB and local minima at SD8 and SD10. Trichodesmium, the most abundant diazotroph enumerated at those stations , is buoyant, and horizontal advection is well known to result in patchy distributions of Trichodesmium in the surface ocean , with huge surface accumulations named slicks, as observed during OUTPACE , surrounded by areas with lower accumulations. These physical processes may explain the differences between stations rather than local enrichments of nutrients due to islands as those three stations where the highest rates were measured are not located close to islands. However, the huge surface bloom observed at LDB (Fig. 1) and extensively studied by was mainly sustained by N2 fixation (secondary fueling picoplankton and diatoms, ), rather than deep nutrient inputs . This bloom had been drifting eastwards for several months and initially originated from Fiji and Tonga archipelagoes (https://outpace.mio.univ-amu.fr/spip.php?article160, last access: February 2016), which may have provided sufficient Fe to alleviate limitation and trigger this exceptional diazotroph bloom.

In MA waters, the dominant diazotroph phylotypes quantified using nifH quantitative PCR assays were Trichodesmium spp. and UCYN-B , which commonly peaked at > 106 nifH copies L-1 in surface (0–50 m) waters. DDAs (mainly het-1, but het-2 and het-3 were also detected) were the next most abundant diazotrophs . This result is consistent with the fact that abundances of those phylotypes co-varied and were significantly positively correlated with N2 fixation rates (Table 2). The two UCYN-A lineages (UCYN-A1 and UCYN-A2) were less abundant (< 1.0–1.5 % of total nifH copies, ) and not significantly correlated with N2 fixation rates (Table 2).

The relative contribution of different diazotroph phylotypes to bulk N2 fixation has been largely investigated through bulk and size fractionation measurements (usually comparing > and < 10 µm size fraction N2 fixation rates), which may be misleading since some small-size diazotrophs are attached to large-size particles or form colonies or symbioses with diatoms (e.g., UCYN-B, ) and some diazotrophic-derived N released by diazotrophs is assimilated by small and large non-diazotrophic plankton (e.g., ). Here we directly measured the in situ cell-specific N2 fixation activity of the two dominating diazotroph groups in MA waters: Trichodesmium and UCYN-B.

At all three studied stations, Trichodesmium dominated, accounting for 68.0–91.8 % of the diazotroph community, followed by UCYN-B, accounting for 0.3–31.7 %. In addition, Trichodesmium and UCYN-B had the highest measured gene expression (102–105 cDNA nifH copies L-1). It was not surprising that UCYN-B had a high gene expression given that the sampling time occurred later in the day (17:00–21:00); however, both Trichodesmium and het-1 (which typically reduce N2 and express nifH highest during the day, ) had a detectable and often equally as high expression as UCYN-B. Cell-specific N2 fixation rates reported here are on the same order of magnitude as those reported for field populations of Trichodesmium and UCYN-B . Trichodesmium was always the major contributor to N2 fixation, accounting for 47.1–83.8 % of bulk N2 fixation, while UCYN-B never exceeded 6.1–10.1 %, despite accounting for > 30 % of the diazotroph community at SD6. This may be linked with the lower 15N enrichment at SD6 (0.517 ± 0.237 at. %), which is due to a high proportion of inactive cells (at. % close to natural abundance) compared to SD2, where the majority of cells were active and highly 15N-enriched (1.163 ± 0.531 at. %). Such heterogeneity in N2 fixation rates among UCYN-B-like cells has already been reported by . Overall, these results show that the most abundant phylotype (Trichodesmium) accounts for the majority of N2 fixation, but not in the same proportion, further indicating that the abundance of micro-organisms in seawater cannot be equated to activity, which has already been reported for other functional groups such as bacteria (Boutrif et al., 2011). In the North Pacific Gyre (station ALOHA), Foster et al. (2013) report a higher contribution of UCYN-B to daily bulk N2 fixation (24–63 %) during the summer season, indicating that this group likely contributes more to the N budget at station ALOHA than in the WTSP, where Trichodesmium seems to be the major player.