Efficient dissolved organic carbon production and export in the oligotrophic ocean.

Biologically fixed carbon is transferred from the surface to deep ocean as sinking particles or dissolved organic carbon (DOC). DOC is estimated to account for ~20% of global export production, but the degree to which this varies regionally has not been assessed at a global scale. Here we present the first observationally based global-scale assessment of DOC production and export, obtained by combining an artificial neural network estimate of the global DOC distribution, and a data-constrained ocean circulation model. Our results demonstrate that the efficiency of DOC production and export varies more than threefold across oceanographic regions. DOC production and export display a pronounced peak in the oligotrophic subtropical oceans, where DOC accounts for roughly half of the total organic carbon export. These stratified nutrient-depleted regions are expected to expand with future warming, amplifying the role of DOC in the biological pump, and magnifying the need to improve DOC cycling in climate models.

Introduction

Organic carbon produced by phytoplankton in the surface ocean is partitioned into a spectrum of sizes, each of which is exported to the subsurface ocean via different mechanisms. In contrast to gravity-controlled export of particulate n class="Chemical">organic carbon (POC), export of dissolved organic carbon (DOC) from the surface to the subsurface is controlled by mixing. The contribution of DOC to the export of biogenic carbon from the surface ocean has been estimated at ~20% globally[1,2], and some studies suggest that DOC can account for more than 50% of the total carbon export on local to regional scales[3,4]. Mechanisms for regional differences in DOC export efficiency (the relative contribution of DOC to the total organic carbon export) are well documented and include variations in ambient nutrient concentrations and microbial community structure[5,6]. Nonetheless, the only basin-scale study of DOC production dynamics found that although the ratio of net DOC production to new production (as measured by the net uptake of nitrate, NO3 −, in the surface ocean) can vary significantly with latitude, a constant ratio was sufficient to explain the surface DOC distribution of the Atlantic Ocean[7,8]. Since the ratio of C:N in sinking particles is also relatively constant[9], it is not clear how the observed variability in DOC export efficiency arises. A global picture of DOC production and export would help to resolve this puzzle, but is currently lacking.

Here we develop a global picture of DOC production and export efficiency by applying an artificial neural network (ANN) to a global ocean n class="Chemical">DOC data set[10]. The ANN takes advantage of the fact that the DOC distribution in the ocean is controlled by the same physical and biological processes which also control well-sampled parameters such as salinity, temperature, macronutrients, chlorophyll, light penetration, and dissolved oxygen. These latter parameters are sufficiently highly sampled that global maps at 1° horizontal resolution have already been produced. By training an ANN using these global maps and the available DOC data, we are able to fill in gaps in the DOC observations and produce global maps of the DOC concentration at 1° horizontal resolution (see Methods). By coupling the ANN-mapped DOC concentrations to a global ocean circulation model, we are able to diagnose substantial regional variability in the efficiency of DOC production and export.

Results

Distribution of DOC

The global surface DOC fields estimated by our ANN model (Fig. 1a) reveal the highest n class="Chemical">DOC concentrations in the tropics and subtropics of ~80 µmol kg−1. Lowest surface DOC concentrations of ~40 μmol kg−1 are found in the Southern Ocean, due to upwelling of DOC-depleted deep-ocean waters. At intermediate depths of 300–600 m, DOC concentrations are the highest in the North Atlantic, where deep-water formation brings DOC to a depth, and in the subtropical gyres where mode waters form(Fig. 1b, c).

Fig. 1

Distribution of ANN-derived and observational DOC. Color map is the artificial neural network (ANN)-derived dissolved organic carbon (DOC) concentration, and colored dots are the observed DOC concentration at 20 m (a), 300 m (b), and 600 m (c). ANN-derived DOC correlates with ~30,000 observed DOC data points with slope = 1.00, R 2 = 0.95, and RMSE = 2.37 μmol kg−1. A correlation study of data points in the upper 74 m indicates slope = 1.00, R 2 = 0.85, and RMSE = 4.36 μmol kg−1

Distribution of ANN-derived and observational DOC. Color map is the artificial neural network (ANN)-derived dissolved n class="Chemical">organic carbon (DOC) concentration, and colored dots are the observed DOC concentration at 20 m (a), 300 m (b), and 600 m (c). ANN-derived DOC correlates with ~30,000 observed DOC data points with slope = 1.00, R 2 = 0.95, and RMSE = 2.37 μmol kg−1. A correlation study of data points in the upper 74 m indicates slope = 1.00, R 2 = 0.85, and RMSE = 4.36 μmol kg−1

DOC production and export

To calculate the net DOC production in the surface ocean (hereafter abbreviated as n class="Chemical">NDP, also referred to as accumulation = [gross production − decomposition]), we used a diagnostic approach similar to the one previously applied to global observational fields of macronutrients[11,12]. In this approach, the ANN-produced DOC field is coupled to a data-constrained global ocean circulation model[13,14], and the production of DOC in the surface ocean is simulated using a restoring model (see Methods section; ref. [12]). We used the same method to calculate net NO3 − uptake in the surface ocean, which at steady state is equivalent to nitrate-fueled new production (see Methods). This method of calculating NDP and NO3 − uptake relies directly on the observed NO3 − and DOC fields, and circumvents the need to model DOC production and decomposition using parameterizations of food-web processes, as done in prognostic ocean biogeochemistry models (e.g., ref. [15]). Our diagnostic approach can thus be used to test some of the assumptions regarding the relationship between new production and NDP that are built into prognostic climate models[16]. Since the export of DOC might take place at a different location from where it is produced, we also performed a separate simulation to diagnose net DOC export from the surface ocean (hereafter abbreviated as NDX; see Methods). One important limitation of this approach is that the majority of the DOC observations are made in summer months, when DOC tends to accumulate at the surface[3], while our steady-state circulation model imposes winter-mixed layer depths in order to achieve proper ventilation of the deep ocean. Our approach thus assumes that all the DOC measured on the surface is preserved to be exported to the subsurface during winter, and thus has a lifetime of at least 3−6 months.

Our results show that net DOC production differs from net n class="Chemical">DOC export at 74 m (the bottom of the second layer of our ocean model), but both show the largest values in the subtropical oceans, roughly between 20° and 40° (Fig. 2). Our diagnostic model finds that the global NDX to below 74 m (which is equal to global NDP within the upper 74 m at steady state) is 2.31 ± 0.60 PgC yr−1, similar to previous estimates[1,2,15]. However, the spatial variability of DOC production shows large deviations from that expected due to nutrient availability. The subtropical gyres are the regions with the lowest rates of net primary productivity (NPP) and POC export[17-19]. However, DOC production is as large as ~15 gC m−2 yr−1 in the subtropics (Fig. 2a), suggesting a decoupling of net DOC production from nutrient availability in these oligotrophic regions. Indeed, our diagnostic calculations of NDP and NO3 − uptake show large spatial variability across the ocean, with zonally averaged NDP:NO3 − uptake ratios (always in mol C:mol N if omitted hereafter) of >15:1 for production above 36 m in the subtropics (~30° N and ~30° S), and ~10:1 for production above 74 m (Fig. 3a). These ratios are even larger than the C:N ratios of exported particulate organic matter (POM), which is roughly constant throughout the ocean at about 6.3−7.7 mol C:mol N[9]. This suggests that DOC production significantly amplifies the total C production (the net community production) in the oligotrophic oceans.

Fig. 2

Maps of net DOC production and export fluxes. Net DOC production (NDP) in the upper 74 m (a) and net DOC export (NDX) below 74 m (b). At steady state, the global summation of NDX is equal to that of NDP, and is 2.31 ± 0.60 PgC yr−1

Fig. 3

Latitudinal variability of DOC production and export to NO3 − uptake and C export. a Zonally averaged ratios of net DOC production (NDP) to NO3 − uptake within the upper 36 m (solid red line with dashed red lines indicating ±1 standard deviation), 74 m (solid blue line with dashed blue lines indicating ±1 standard deviation), and the deepest euphotic zone over the course of the year (Z eu,max; solid green line with dashed light green lines indicating ±1 standard deviation). Also shown is the global average ratio of C:N in POM (gray shading; ref. [9]) and the mean NDP:NO3 − uptake ratio estimated in a previous study for the Atlantic Ocean (solid brown line; refs. [7,8]). b Ratios of zonally averaged NDP (solid red line with dashed light red lines indicating ±1 standard deviation) and net DOC export (NDX; solid blue line with dashed light blue lines indicating ±1 standard deviation) to satellite-derived C export flux

Previous studies have shown that the ratio of carbon to nutrients in DOM is higher than the standard Redfield C:N ratio of ~6.6:1[20]. Prognostic models have used a C:N ratio ~11.8:1 for DOM[15]. Field studies also indicate a high C:N ratio of ~10 in DOM[21]. Our study offers the first global data-based assessment of the ratio between n class="Chemical">NDP and NO3 − uptake, and shows that the C:N ratio in DOM can vary widely between regions and vertically within the euphotic zone. The NDP:NO3 − uptake ratio is highest near the sea surface and decreases with depth in the euphotic zone. NDP:NO3 − uptake within the upper 36 m is ~15–20 in the subtropical gyres, and reaches a maximum of >30 in the subtropical South Pacific (Supplementary Fig. 1). The NDP:NO3 − uptake ratio is lower when averaged over the top 74 m, but still reaches values of ~7–10 in the subtropical gyres (Fig. 3a). Inasmuch as a portion of new NO3 − input is exported as POM, the C:N ratio in the DOM pool will be even higher than these NDP:NO3 − uptake ratios in the subtropics. Indeed, a recent observational study demonstrated that the C:N ratio of total organic matter was at least twice that of POM in the North Atlantic subtropical gyre[22], suggesting that the C:N ratio of DOM is >13:1 in this region. In all, our results show that the nutrient efficiency of DOC accumulation in the upper ocean varies substantially with depth, region, and ocean basin (Supplementary Fig. 1). These results demonstrate that conditions such as temperature[23,24], ambient nutrient concentration[6], light availability[25,26], or bacterial community composition[5] might regulate the effect of new nutrients[7] in fueling DOC accumulation.

DOC contribution to total carbon export

To quantify the relative importance of DOC export to the biological pump, both globally and regionally, we compared our diagnosed net n class="Chemical">DOC production and export fluxes with C export fluxes derived from satellite measurements and empirical models (see Methods; refs. [17-19]). The ratios of NDP and NDX to satellite-derived C export show a pattern similar to the diagnostically derived NDP:NO3 − uptake, with the highest values at subtropical latitudes (~15–30° N and S; Fig. 3b). In the subtropics, the ratio of NDP:C export is ~0.2–0.4, while in the equatorial region, this ratio is ~0.15, and in the high latitudes, it is ~0.1 or less (Fig. 3b). The latitudinal variation in NDP:C export is not quite as extreme as that of NDP:NO3 − uptake, which could reflect an additional source of new nutrients to the subtropical gyres from nitrogen fixation[27]. The pattern of NDX:C export is similar to that of NDP:C export, but export ratios are lower in the tropics and higher in the subtropics (Fig. 3b). This reflects the large-scale wind-driven circulation in the low latitudes: upwelling prevents DOC export near the equator, while poleward surface flows transport DOC accumulated in the tropics to the subtropics, where it is subducted. This sustains the balance between global DOC production and export at steady state (Supplementary Fig. 2). In the high latitudes, the production and export of DOC are nearly identical, reflecting the export of locally produced DOC (Fig. 3b). The globally-averaged ratio of DOC export to C export is 0.20 ± 0.09, which is in agreement with previous studies which suggest that DOC accounts for about 20% of organic matter export[1]. Our results are also in broad agreement with tracer-based estimates of total carbon export at two time-series study locations in the subtropical North Atlantic and subtropical North Pacific, both of which show a substantial contribution of DOC export total C export[4].

Discussion

What are the causes of the efficient DOC export in the subtropical gyres? One hypothesis that correlates well with the observations presented here is the inhibition of n class="Chemical">DOC consumption by bacteria under low-nutrient conditions[5,6]. To examine the correlation of NDP efficiency with nutrient availability, we partitioned the euphotic zone into 10 different regions based on their nutrient concentration and ocean basin (Fig. 4a; see Methods section). Averaging the NDP:C export ratio over the resulting regions reveals that the most oligotrophic ocean regions, the subtropical gyres which have NO3 − concentrations less than 0.5 μM, have NDP:C export ratios of ~0.35–0.5 (Fig. 4b). The ratio of DOC export:C export is even higher (~0.45–0.7) than NDP:C export in subtropical regions, reflecting the additional subduction of DOC produced in the tropics. Throughout most of the rest of the ocean, the NDP:C export ratio ranges from only 0.1 to 0.2 (Fig. 4b). Finally, the ratio of NDP:NPP ranges from 0.02 to 0.08, also attaining the highest values in the subtropical gyres (Fig. 4b).

Fig. 4

Regional variability of DOC production and export to NPP and C export. a Distribution of depth-integrated mean of nitrate above 74 m, which is used to regionalize the global ocean into distinct regions (solid borderlines): Antarctic Zone (AAZ), Sub-Antarctic Zone (SAZ), Sub-Tropical Pacific (STP), Atlantic (STA), and Indian (STI), and Tropical Pacific (TP), Atlantic (TA), and Indian (TI), and North Pacific (NP) and Atlantic (NA). b Ratio of regionally averaged NDP and NDX fluxes to satellite-derived net primary productivity (NPP) and C export fluxes for the regions shown in panel a. Error bars in b indicate ±1 standard deviation

Regional variability of DOC production and export to n class="Chemical">NPP and C export. a Distribution of depth-integrated mean of nitrate above 74 m, which is used to regionalize the global ocean into distinct regions (solid borderlines): Antarctic Zone (AAZ), Sub-Antarctic Zone (SAZ), Sub-Tropical Pacific (STP), Atlantic (STA), and Indian (STI), and Tropical Pacific (TP), Atlantic (TA), and Indian (TI), and North Pacific (NP) and Atlantic (NA). b Ratio of regionally averaged NDP and NDX fluxes to satellite-derived net primary productivity (NPP) and C export fluxes for the regions shown in panel a. Error bars in b indicate ±1 standard deviation

Besides nutrient concentrations, other factors which can influence DOC production and export efficiency include temperature[23,24] and ecosystem structure[5]. To investigate the effects of all of these factors on n class="Chemical">DOC export efficiency, we regressed the DOC production efficiency against the export-weighted mean temperature, NO3 −, and fractional picoplankton abundance (F pico) averaged over the top 74 m within each of our 10 regions (Fig. 5). Here, picoplankton abundance stands as a proxy for ecosystem structure. Picoplankton are small plankton (<2 μm in diameter) which are abundant in nutrient-poor stably stratified ecosystems. Some picoplankton like Prochlorococcus have even lost the ability to utilize nitrate[28], but remain efficient producers of many types of dissolved organic matter compounds[29].

Fig. 5

Regional relationship between DOC production:C export and several environmental parameters. Our diagnosed NDP:C export in the upper 74 m is plotted against regionally averaged temperature (a), F pico (fraction of picoplankton) (b) and log10(nitrate) (c). The diagnosed NDP:C export ratios can be explained using a multilinear regression (MLR) model with two predictor variables (F pico, and log10(nitrate)) (d). Error bars in a–c indicate ±1 standard deviation. Dashed envelope in d is the 95% confidence interval of the MLR prediction

Regional relationship between DOC production:C export and several environmental parameters. Our diagnosed n class="Chemical">NDP:C export in the upper 74 m is plotted against regionally averaged temperature (a), F pico (fraction of picoplankton) (b) and log10(nitrate) (c). The diagnosed NDP:C export ratios can be explained using a multilinear regression (MLR) model with two predictor variables (F pico, and log10(nitrate)) (d). Error bars in a–c indicate ±1 standard deviation. Dashed envelope in d is the 95% confidence interval of the MLR prediction

While previous studies have found increased partitioning of primary production to the DOC pool with increased temperature[23,24], our regional analysis shows that net n class="Chemical">DOC production:C export does not increase monotonically as temperature increases, due to the fact that NDP:C export is low in tropical regions which have the highest sea-surface temperatures (Fig. 5a). On the other hand, our analysis shows a strong positive correlation of the NDP:C export ratio with F pico (Fig. 5b), and a negative correlation between NDP:C export and NO3 − concentration (Fig. 5c). To combine the effect of all these parameters, we performed a multilinear regression (MLR) analysis using different combinations of these predictors. The MLR tests indicate that F pico and log10(NO3 −) are the most important predictors of NDP:C export ratio. Adding temperature as a predictor does not substantially improve the regression model (Supplementary Table 1). The derived regression equation NDP:C export = 0.710(±0.042) F pico – 0.101(±0.021) log10(NO3 −) (Supplementary Table 1) reproduces the regionally averaged NDP efficiency with an R 2 of 0.89 (Fig. 5d). The ratio of standardized coefficients of this MLR equation (i.e.,  = −1.17) indicates that the two predictors are almost equally important, but have opposite correlations with NDP:C export.

Although our diagnostic model does not reveal the mechanisms for the remarkably high DOC production efficiency in the oligotrophic subtropical ocean, high n class="Chemical">NDP efficiency appears to correlate significantly with ecosystem structure, as measured by the fractional picoplankton abundance. This suggests that ecological factors drive high DOC export from these regions. The stratified subtropical gyres are characterized by an efficient microbial loop, which processes NPP through the marine food web in the euphotic zone, building up a pool of more refractory DOC as fresh organic matter is degraded[30]. This refractory DOC only becomes available to marine microbes upon export out of the euphotic zone[3]. In contrast, the tropical and subpolar regions where nutrients are abundant are characterized by larger plankton communities, which promote the formation of fast-sinking particle aggregates and fecal pellets[31] and a less-efficient food web. As well, high availability of nutrients in these regions facilitates the microbial consumption of DOC[6]. Finally, the substantial increase in DOC production efficiency in the upper water column (Fig. 3a) indicates that high light levels may also play an important role in stimulating DOC production[25,26], and/or in inhibiting NO3 − uptake[28].

The role of DOC export with respect to the surface bioassimilation of nutrients needs to be better considered in prognostic climate models[16]. Some of the current climate models are able to produce a pattern of surface n class="Chemical">DOC accumulation similar to that of observations[1]. However, a closer look at the ratio of simulated surface DOC to nutrient concentrations from these models shows a several-fold mismatch among the models, which are all significantly at odds with the observations (Supplementary Fig. 3). Our results highlight several mechanisms that could improve model estimates of carbon export, including tying DOC production and remineralization to light availability, nutrient concentration, and phytoplankton community structure.

The remarkably efficient production and export of DOC in the nutrient-depleted waters of the subtropical ocean suggest a vital role for n class="Chemical">DOC in the future ocean. As nutrient-poor regions expand with global warming and the associated stratification of surface waters[32-34], efficient low-nutrient-demanding export of DOC from these regions will be an ever more important contributor to the oceanic biological sink of atmospheric CO2. Indeed, observations in the North Pacific subtropical gyre are already showing a long-term increase in surface DOC accumulation coincident with enhanced stratification[35]. Our analysis suggests that these trends will hold globally, amplifying the role of DOC in the ocean’s biological carbon sink in the future.

Methods

Observational data sets

The temperature, salinity, oxygen, and nutrient data were obtained from World Ocean Atlas 2013_v2 1° × 1° annual climatology (https://www.nodc.noaa.gov/OC5/woa13; refs. [36-39]). Monthly SeaWiFS n class="Chemical">chlorophyll maps were downloaded from http://www.science.oregonstate.edu/ocean.productivity/index.php for the period of 1997–2009, and averaged to obtain a monthly climatology. Monthly SeaWIFS chlorophyll data were also used to calculate the depth of the euphotic zone according to the relationship of ref. [40]. DOC data were obtained from a recent compilation[10] plus CLIVAR repeat section A10 (completed in 2011) which was downloaded from http://cdiac.ornl.gov. Ocean bathymetric data were downloaded from http://www.gebco.net.

ANN implementation

The global data set of DOC observations compiled by Letscher and Moore[10] was binned into the World Ocean Atlas 2013_v2 1° × 1° grid (https://www.nodc.noaa.gov/OC5/woa13/). Approximately 30,000 grid points were represented by one or more observations. These n class="Chemical">DOC data points were taken as the target to train an ANN. Input data used for the ANN included nitrate, phosphate, apparent oxygen utilization, silicate, dissolved oxygen, salinity, temperature, potential density anomaly, depth, bottom depth, depth of the euphotic zone, and chlorophyll concentration, each provided on the same 1° × 1° grid as the DOC observations. In each training cycle, 70% of data were sampled and used for training and 30% were used for test and validation. We used a single hidden layer with 10–20 neurons fully connected, in a feed-forward architecture, to a single-node output layer. A sigmoid response function was used for the hidden layer, and a linear response function was used for the output layer. Levenberg–Marquardt and Bayesian regularization methods were used for back propagation. This architecture was determined to be the simplest ANN that could be trained to reproduce validation data sets with acceptable correlation metrics, and with no sign of overfitting.

An ensemble of 20 individual ANNs were trained by randomly sampling 70% of the data set, each of which were able to produce the validation subset (30% of withheld observations) with R 2 ≥ 0.95. Figure 1 shows the average of these 20 ANN-produced DOC fields. In addition to excluding 30% of the data during each training process for validation, we also excluded n class="Chemical">DOC data along CLIVAR A10 section during all the ANN training sessions. The correlation between the ANN-derived DOC and observational DOC along CLIVAR A10 section yields a slope, R2, and RMSE of 1.02, 0.87, and 3.47 μmol kg–1, respectively (Supplementary Fig. 4). Globally averaged profiles of mean difference and root-mean-square error (RMSE) between observed and ANN-derived DOC reveal almost no bias, with RMSE ranging from ~1 μmol kg–1 in the deep ocean to ~5 μmol kg–1 on the surface (Supplementary Fig. 5).

In a critical examination of the ANN approach, we excluded data from whole-ocean basins from our training subset, and examined the fit of the resulting model to data from the excluded basin. In these strong tests of the ANN approach, the ANN-predicted and observed DOC compare favorably in each basin (Supplementary Fig. 6). R 2 values for the ANN-predicted vs. observed n class="Chemical">DOC concentrations in the excluded basins vary from 0.8 (North Atlantic and Indian Ocean) to >0.9 (South Pacific and South Atlantic), indicating the robustness of the ANN approach. We should note that this test is especially strict for the Indian Ocean, since it contributes more than 30% of the total DOC observations.

Diagnostic modeling

The ocean circulation model provides an off-line transport matrix that has been constrained by observational fields of salinity, potential temperature, natural ∆14C, CFC-11, sea-surface n class="Disease">height anomaly, sea-surface heat, and freshwater fluxes[13,14]. After interpolation of the ANN-derived DOC and climatological nitrate fields from WOA13 grids (1° × 1° and 102 depth levels) to our ocean circulation model (2° × 2° and 24 depth levels), the following matrix equation that is a steady-state mass balance was solved:

where T ss, T si, T is, and T ii (units of yr−1) are the partitions of the transport matrix that circulate the tracer C (DOC or n class="Chemical">NO3 −, units of µmol kg−1) from surface to surface, surface to interior, interior to surface, and interior to interior grid points, respectively. C s and C i are the tracer concentration at the surface and interior grid points, respectively. J and L are the net biological source or sink (i.e., production–decomposition for DOC, or uptake– regeneration for NO3 −) of C on the surface and interior, respectively (units of µmol kg−1 yr−1). Since we seek only the J term, but not L, we can substitute the observed tracer concentration for C i and write the first row of Eq. (1) separately, as

Here, J is the net uptake rate (for NO3 −; J up, nitrate) or net production rate (for DOC; J pr, DOC; NDP rate) on the surface which is modeled by restoring to observations with a timescale τ (units of yr):

Cs, obs and C i, obs were set equal to the World Ocean Atlas climatology for NO3 − and to our ANN-derived climatology for DOC, and upon choosing a value of (see below), we solved Eqs. 2–4 for C s and for J up, nitrate and J pr, DOC for each grid point in the surface ocean.

To investigate the depth variability of net DOC production within the euphotic zone (results shown in Fig. 3a and Supplementary Fig. 1), we ran three separate models using different depths to separate the surface and interior grid points. In the first model, we set the depth of the surface at the base of the first model layer (36 m). In the second model, we set the depth of the surface at the base of the second model layer (74 m). In the third model, we set the depth of the surface at the base of the model layer that is closest to (but not deeper than) the maximum depth of the euphotic zone. In this third model, the depth of the surface varies spatially due to variability in the depth of the euphotic zone. We emphasize that this model is entirely diagnostic, and does not rely on any assumptions regarding ecological or chemical processes in the upper ocean.

Wherever applicable, we used the density of seawater (from the equation of state of seawater using temperature, salinity, and pressure) to convert seawater mass-based units to seawater volume- or area-based units.

Net DOC export calculation

To assess the net export of DOC (n class="Chemical">NDX) from the euphotic zone, we solved the equationwhere J pr, DOC is the rate of net production of DOC (NDP) derived from our diagnostic model (Eqs. 2 and 4), and k is set to zero everywhere except in the layer directly below the export depth (74 m in Figs. 2 and 3), where it is set to a very large value (1012 yr–1) in order to effectively remove all tracer (C) within the grid box directly below the export depth. Solving this equation for C at steady state yields the rate of export (; in µmol kg−1 yr−1) needed to balance net DOC production above the export depth. The molar rate was converted to mass flux (gC m−2 yr−1) by applying the thickness of the third layer in the model, seawater density, and the atomic weight of carbon.

Monte Carlo uncertainty assessment

The above equations were solved 2,000 times with random inputs of ANN-derived DOC field (i.e., one n class="Chemical">DOC field out of 20 possible DOC fields was chosen for each run), circulation matrix (10 different circulation matrices from ref. [13]), and τ between 3 and 9 months. We note that using the 20 different DOC fields produced by the ANN in our model calculations allows us to propagate the spatially correlated errors of the ANN-predicted DOC fields. The results from these 2000 different models were taken to calculate averages and standard deviations of DOC production and NO3 − uptake. Uncertainty assessment of the export calculation was also performed in the same manner where the above production solutions and 10 different circulation matrices were used to solve Eq. 5.

Satellite-based estimation of C export and Fpico

Satellite-based C export was calculated using three empirical equations which require only temperature and NPP[17-19]. SeaWiFS n class="Chemical">NPP was downloaded from http://www.science.oregonstate.edu/ocean.productivity/index.php for all three NPP algorithms (VGPM[41], Eppley-VGPM[42], and CbPM[43]) and averaged over all months and years to produce an annual climatology. Overall, nine different satellite-derived C export maps were calculated using different combinations of NPP and export algorithm. These were averaged for Figs. 3b and 4b in the main text. F pico was calculated from our derived annual NPP from 1997 to 2009 using a previously established relationship[31].

Regionalization

The regions in Fig. 4 were divided based on the depth-integrated mean of NO3 − concentration between 0 and 74 m. Subtropical regions in the Atlantic, Pacific, and Indian Oceans (STA, STP, and STI, respectively) were divided from tropical (TA, TP, and TI), northern (NA and NP), and southern (SAZ) regions at a mean n class="Chemical">NO3 − concentration of 0.5 μmol kg–1. We note that we used 0.5 μmol kg−1 nitrate concentration only as a convenient marker for the edge of the oligotrophic gyres in the annual nitrate climatology, and not to imply that this particular value of NO3 − results in nitrate-limited ecosystems. In the Southern Ocean, the Sub-Antarctic zone (SAZ) was divided from the Antarctic Zone (AAZ) at a mean NO3 − concentration of 20 μmol kg−1.

CMIP5 result analysis

The results of CMIP5 models GFDL-ESM2M[44], CNRM-CM5[45], IPSL-CM5A-LR[46], and CESM1-BGC[47] were downloaded from http://pcmdi9.llnl.gov. Annual data sets of n class="Chemical">DOC and NO3 − from historical simulations were averaged over the time period 1956–2005. The resulting DOC and NO3 − fields were then averaged vertically for the upper 100 m, and then zonally. The ratios of the final average for DOC and NO3 − were then calculated and plotted in Supplementary Fig. 3. Before the ratio calculation, 38 μmol kg−1 DOC (the concentration of refractory DOC observed in the deep North Pacific Ocean) was added to the zonal average of DOC for the CNRM-CM5, IPSL-CM5A-LR, and CESM1-BGC models, which do not include a refractory DOC pool[1].

Map visualization

Color maps in Figs. 1, 2, and 4 were produced using Ocean Data View 4.7.7 software[48]. The gridding was performed through weighted averaging applying 15‰ for both X and Y axes scale-lengths.

Data availability

The data produced in this study are available from the corresponding author upon request.