Quantifying the Impact of Climate Change on Marine Diazotrophy: Insights From Earth System Models
Lewis Wrightson* 
Alessandro Tagliabue*- Department of Earth, Ocean and Ecological Sciences, University of Liverpool, Liverpool, United Kingdom
Nitrogen fixation is a major source of new nitrogen to the ocean, supporting biological productivity in the large nitrogen-limited tropical oceans. In Earth System Models, the response of nitrogen fixation to climate change acts in concert with projected changes to physical nitrogen supply to regulate the response of primary productivity in nitrogen-limited regions. We examine the response of diazotrophy from nine Earth System Models and find large variability in the magnitude and spatial pattern of nitrogen fixation in both contemporary periods and future projections. Although Earth System Models tend to agree that nitrogen fixation will decrease over the next century, strong regional variations exist, especially in the tropical Pacific which may counteract the response of the Atlantic and Indian oceans. As the climate driven trend of nitrogen fixation emerges by mid-century in the RCP8.5 scenario, on regional scales it may modulate the broad climate trends in productivity that emerge later in the century. The generally poor skill and lack of agreement amongst Earth System Models indicates that the climate response of nitrogen fixation is a key uncertainty in projections of future ocean primary production in the tropical oceans. Overall, we find that the future evolution of nitrogen fixation plays an important role in shaping future trends in net primary production in the tropics, but the poor skill of models highlights significant uncertainty, especially considering the role of multiple concurrent drivers.
Introduction
Future warming is projected to enhance ocean stratification, reducing the vertical supply of nutrient-rich subsurface waters to the nutrient-depleted surface waters of the euphotic zone in the tropical ocean in particular (Steinacher et al., 2010; Laufkötter et al., 2015). Two resources of particular importance in the tropics are phosphorus (P) and nitrogen (N), which are required for the production of DNA, ATP, proteins, pigments, and enzymes. A reduction in the vertical supply of N and P to the surface ocean due to ocean warming will cause phytoplankton to become nutrient limited, thereby restricting their growth and net primary production (NPP), with impacts on the biological carbon pump (Bopp et al., 2013; Moore C.M. et al., 2013; Moore et al., 2018). A common way to illustrate the impact of climate change on ocean properties is by using the concept of “time of emergence,” which assesses the time when the climate change signal in a given property emerges from the background variability. Such analyses illustrate that changes linked to NPP can emerge much later than those associated with sea-surface temperature (SST) and pH (Keller et al., 2014; Frölicher et al., 2016; Henson et al., 2017).
As dinitrogen (N2) fixation is a globally significant N supply mechanism (68–164 Tg N yr–1, Gruber and Sarmiento, 1997; Jickells et al., 2017; Luo et al., 2014; Luo et al., 2012; Tang et al., 2019a; Wang et al., 2019), future changes in N2 fixation by diazotrophs may modulate any N limitation arising from enhanced stratification. This means that future changes to diazotrophy, and the supply of new fixed N, will be an important component of the response of tropical ocean productivity to climate change in Earth System Model projections. The process of marine N2 fixation involves the reduction of N2 gas to ammonia (NH3) by marine organisms known as diazotrophs (Sohm et al., 2011). Fixing N2 gas requires the iron-rich enzyme nitrogenase, which requires 5 to 100-fold more iron (Fe) than non-diazotrophs (Raven, 1988; Kustka et al., 2003). N2 fixation also has an associated high energetic cost (Großkopf and Laroche, 2012), and like all phytoplankton, diazotrophs are susceptible to P limitation (Sañudo-Wilhelmy et al., 2001; Sohm et al., 2008). Overall, the need for light, Fe, and P, tends to confine diazotrophy to regions relatively rich in Fe and high light availability (such as the subtropical North Atlantic), and regions with ample P (such as the Western Pacific and subtropical South Atlantic, Falkowski, 1997; Mills et al., 2004; Goebel et al., 2008; Moore et al., 2009; Sohm et al., 2011). In N limited regions of the ocean, diazotrophs are only able to outcompete and proliferate provided that sufficient Fe and P are available (Moore et al., 2009; Dutkiewicz et al., 2012).
The most well studied diazotroph is the colony forming cyanobacterium, Trichodesmium, which occurs predominantly in the subtropical and tropical ocean (Sohm et al., 2011). Trichodesmium is one of a few diazotrophs that is able to be cultured in the laboratory allowing for detailed analysis of its physiology (Hutchins et al., 2013) and, it is one of the major contributors to global marine nitrogen fixation supplying ∼60–80 Tg N yr–1 (Mahaffey et al., 2005). For these reasons, Trichodesmium has been used as a model organism to parameterize diazotrophy in several global ocean models. These include the Earth System Models used for climate projections, which tend to restrict rates of diazotrophy based on light, Fe, P, and N availability. Recent laboratory work has demonstrated that diazotrophs are likely to be responding to a whole suite of interacting drivers such as temperature, CO2 and nutrient co-limitation in the future (Hutchins et al., 2007, 2013; Walworth et al., 2016). Moreover, the contributions of unicellular diazotrophs such as Crocosphaera and the unicellular cyanobacterial group A (UCYN-A) which have high abundances and broad ecological ranges (Moisander et al., 2010), have not been incorporated into Earth System Models. UCYN-A for example has been shown to occupy a larger ecological range than Trichodesmium and may be a more important player in marine N2 fixation away from the tropics (Martinez-Perez et al., 2016). It has recently been identified that heterotrophic diazotrophs are widespread and highly diverse (Bombar et al., 2016; Moisander et al., 2017), with the ability to fix N2 in environments that are traditionally believed to exclude diazotrophy (Rahav et al., 2013; Benavides et al., 2016). Heterotrophic diazotrophs are also not currently incorporated into Earth System Models, which may indicate that these models are underestimating the areal extent of N2 fixation by neglecting the niche occupied by heterotrophic diazotrophs.
As part of the IPCC 5th assessment report (IPCC-AR5), 9 of the Earth System Models used to project future changes in ocean biogeochemistry included marine N2 fixation (Ciais et al., 2014). These models represent diazotrophy with a range of complexity (Supplementary Table 1), at the simplest only including the effect of temperature and N inhibition, while others incorporate singular limitation by Fe and P. More recently, climate projections have also been made using models with variable non-diazotroph cellular stoichiometry (IPSL-Quota, Kwiatkowski et al., 2018). In this study, we compare 8 IPCC-AR5 models and the more recent variable cellular stoichiometry model IPSL-Quota to examine how ocean diazotrophy is projected to respond to climate. While previous studies have conducted inter-model comparisons for N2 fixation, they used arbitrary scaling factors for some models (Riche and Christian, 2018; Tang et al., 2019a) and did not assess the timescales and regionality of change. In this study, we update the reported global trends using the correct units for all models, address the regional distribution of N2 fixation and explore how the time of emergence of N2 fixation varies on global and regional scales both within and across models.
Methods
Model outputs of depth integrated N2 fixation from the coupled model intercomparison project 5 (CMIP5) that compiled model outputs were obtained from 8 models, alongside a more recent quota-based formulation of the IPSL-CM5A-LR model from Kwiatkowski et al. (2018). For each model, data from the historical (1850–2005) and future (RCP 8.5, 2006–2100) scenarios was re-gridded to fit a standard 1°× 1° global grid. To calculate differences in N2 fixation rates, both spatially and globally integrated, we used 20-year averaged time slices from the historical (1986–2005) and RCP 8.5 (2081–2100) simulations for all 9 models. We further investigated regional patterns of depth integrated N2 fixation rates using an established biogeochemical province approach (Longhurst, 2007) and derived the average trend of depth integrated N2 fixation for each province. We then compared the percentage of models that agreed with the average trend to interpret model agreement. We note here that prior work used an arbitrary scaling of 0.17 for the two IPSL model in the CMIP5 archive (Riche and Christian, 2018; Tang et al., 2019a). This is due to IPSL output being stored in the CMIP5 archive in carbon, rather than nitrogen, units. We here report results from those models using the correct adjustment of 1/7.625 (or 0.131), which reflects the model N/C ratio (Aumont and Bopp, 2006). Models vary markedly in the way in which they represent N2 fixation. Some represent an explicit diazotroph functional group, which obtains fixed N from N2 fixation. Others include N2 fixation implicitly as a source of NH3 under certain conditions. While others simply apply a N flux that restores the N:P ratio back to Redfield values of 16:1 (Supplementary Table 1).
Time of emergence (ToE) analysis was performed to detect when the climate driven N2 fixation trend emerges from the natural variability within the models (e.g., Keller et al., 2014). The ToE rests on the concept of comparing the temporal climate change signal to the background noise, which is determined using different thresholds. Similar to Frölicher et al. (2016), we assume that the N2 fixation signal had emerged when it was greater or less than the mean of the preindustrial control simulation ± 1 standard deviation consistently for 4 years in the direction of the overall trend (to avoid the false diagnosis of shorter oscillations). The sensitivity of the time of emergence to our noise threshold was assessed by also using 2 standard deviations from the preindustrial control simulation. The ToE at each model grid cell was calculated by comparing the annual mean global N2 fixation time series from the RCP8.5 scenario (2006–2100) to the parallel preindustrial control simulations which had been run for at least 300 years. We averaged the ToE across each biogeochemical province to facilitate comparisons.
Results and Discussion
Contemporary Nitrogen Fixation: Spatial Patterns and Global Trends
As in prior work (Riche and Christian, 2018; Tang et al., 2019a), there are large variations between models regarding the magnitude and spatial pattern of diazotrophy in the contemporary ocean (1986–2005). Across the range of models, the updated magnitude of total global N2 fixation using the correct units for all models, ranged 2.5-fold from 76.5 to 185.6 Tg N yr–1. Despite this range, the models broadly fell into two groups based upon their spatial patterns. The first group of models placed emphasis on higher rates of N2 fixation in the upwelling regions of the ocean, while the second group of models placed the majority of N2 fixation in the low nutrient subtropical gyres (Figure 1). Some models (CanESM2, CESM1-BEC, IPSL-CM5A-LR, IPSL-CM5A-MR, and IPSL-Quota) imposed a minimum temperature restriction, which restricts diazotrophs to low latitudes, but for those models without thermal limits, diazotrophy can occur over the entire global ocean (albeit at very low rates at high latitudes, Figure 1) this may be an imperfect solution to account for diazotrophy in the high latitudes. These models without thermal limits on diazotrophy tended to fall into the upwelling group which on the whole had larger global N2 fixation values than gyre models. One exception to this general rule is the N only model CanESM2, which does not include the role of secondary limiting factors, such as Fe and P which exert more control on diazotrophy than N. The large disagreement between the models suggests that our understanding of spatial distribution of N2 fixation is limited. Interestingly, the recent re-evaluation of global ocean N2 fixation tends to support the concept of the greatest rates of N2 fixation occurring in low nutrient subtropical gyres rather than in upwelling regions (Wang et al., 2019).
Figure 1. Average nitrogen fixation rates for the final 20-year time slice (1986–2005) of the historical scenario for 9 global nitrogen fixation models: (A) IPSL-CM5A-LR, (B) CESM1-BEC, (C) CanESM2, (D) IPSL-CM5A-MR, (E) GFDL-ESM2G, (F) GFDL-ESM2M, (G) IPSL-Quota, (H) MPI-ESM-LR, and (I) MPI-ESM-MR. Values above maps indicate the total globally integrated nitrogen fixation rates for the same time period. Red boxes indicate “upwelling” models and blue boxes indicate “gyre” models.
The difference in the spatial pattern between gyre models and upwelling models may be reflected in differences in the rate of total net primary production. For instance, when the ratio of N2 fixation to net primary productivity is plotted across the models we find stronger spatial agreement than for N2 fixation alone (Supplementary Figure 1). N2 fixation is a small fraction of net primary production in the upwelling regions for all models, increasing to a larger fraction in the gyres. This suggests that the higher rates of N2 fixation associated with upwelling regions in the GFDL and MPI models is a result of a much larger rate of upwelling zone NPP in these models, which is allied to a stronger contrast in net primary production between the upwelling and gyre regions.
Future Changes to Marine Diazotrophy: Spatial Patterns and Global Trends
Across the models evaluated, 8 predict a decrease or very little change in global N2 fixation between 2006 and 2100 (between +1 and −31%), with only one model (CanESM2) projecting an increase in diazotrophy. This is because CanESM2 has no limitation by Fe and P, and in this model N2 fixation responds positively to future ocean warming (Figure 2A). The N2 fixation trend is not consistent between models when assessed over the Atlantic, Pacific and Indian ocean basins, such that the global trend can mask important inter-basin variability. For example, stronger trends were predicted for the Atlantic and Indian ocean basins and were relatively consistent with the direction of the global trend (Figures 2B,D). However, the Pacific Ocean exhibited the most variability, with the sign of the trend being opposite to the global trend for the majority of models (Figure 2C). This suggests that diazotrophy in the Atlantic and Indian oceans will be negatively affected by climate change, but the Pacific Ocean may counteract the decrease occurring in the other basins in some cases. It is, however, difficult to assess how realistic these trends are as currently there are very few observations of N2 fixation in the Indian Ocean basin (Luo et al., 2012; Tang et al., 2019b).
Figure 2. Percentage anomaly of annual nitrogen fixation over the time period 2006–2100 for 4 regions: (A) the global ocean, (B) Atlantic Ocean, (C) Pacific Ocean, and (D) Indian Ocean. 9 nitrogen fixation models: CanESM2 (red), CESM1-BEC (blue), GFDL-ESM2G (light green), GFDL-ESM2M (gray), IPSL-CM5A-LR (cyan), IPSL-CM5A-MR (magenta), IPSL-Quota (orange), MPI-ESM-LR (yellow) and MPI-ESM-MR (dark green). Black dashed line indicates the zero line.
The changes in the total global N2 fixation rate also exhibited large variability once updated from prior work to use the correct units for all models, ranging from −50.1 Tg yr–1 to +58.0 Tg yr–1 (Table 1). When models are grouped based on their contemporary pattern of diazotrophy, upwelling models had very little coherency for the spatial pattern of diazotrophy (comparing the 2 MPI and 2 GFDL models, Supplementary Figure 2). For gyre models, however, moderate coherence existed between the 3 IPSL models and the CESM1-BEC model (CanESM2 was overlooked hereafter as an outlier). This shows that although upwelling models display a similar spatial pattern for contemporary diazotrophy, there is a little agreement on how the spatial pattern of diazotrophy will change, suggesting that mechanisms controlling diazotrophy in these models may actually be very different. For example, the 2 GFDL models impose Fe limitation on diazotrophy but the MPI models do not. On the other hand, the relatively good agreement between the gyre based models matches their common underlying drivers (Supplementary Table 1).
Table 1. Summary of model analysis: Total N2 fixation anomaly [RCP8.5 (2081:2100) – historical (1986:2005)], time of emergence of global N2 fixation when the natural variability is calculated using 1 standard deviation (σ) and 2 standard deviations (2σ) of the control simulation, model skill represents the correlation between model depth integrated N2 fixation (1986:2005) and observed depth integrated N2 fixation (Tang et al., 2019b).
Taking a province-based view to synthesize the temporal and spatial model comparison revealed regions where the average N2 fixation trend for the 9 models was projected to be strongly modified. The strongest decreases in diazotrophy occurred in the equatorial Atlantic, western equatorial Pacific and northern Indian. Provinces with the strongest decreases in diazotrophy coincided with provinces with the highest inter-model agreement (Figures 3A,B). The strongest overall increase in diazotrophy was found in the Pacific upwelling region, however, this province coincided with relatively low inter-model agreement (Figures 3A,B). In general, the provinces with lowest agreement corresponded to regions of low average trend, suggesting that these lower average trends were a consequence of divergence between models.
Figure 3. Regional comparison of the 9 global nitrogen fixation models using Longhurst provinces: (A) average trend of nitrogen fixation across all 9 models in each province, (B) percentage of models that agree with the average trend for each province, (C) average time of emergence across all 9 models for each province using 1 standard deviation of the control and (D) the contribution of the nitrogen fixation anomaly to the NO3 anomaly.
Future Changes to Marine Diazotrophy: Time of Emergence
The N2 fixation trend due to climate change emerges from the natural variability across the models between 2018 and 2040 when the natural variability is calculated using one standard deviation of the control simulation (Table 1). If we allow greater natural variability (by using 2 standard deviations), then the emergence is delayed by 5–20 years (Table 1 and Supplementary Figure 3). There is no clear link between the magnitude of the projected change in future N2 fixation and the time of emergence, however, when considering both sensitivities (1 and 2 standard deviations) the gyre based models have better inter model agreement of emergence (2018–2042) than the upwelling based models, which have a wider range of emergence times (2020–2054). Regionally, the N2 fixation trend emerges earliest (∼2030–2050) in the lower latitudes and much later (∼2080–>2100) in the high latitude Southern Hemisphere (Figure 3C), suggesting that N2 fixation in the low latitudes is most susceptible to change. The high latitude emergence of the N2 fixation signal may be difficult to assess as not all models allow high latitude N2 fixation in their mean state. On the whole N2 fixation decreases (e.g., in the tropical Atlantic) or shows very little change (e.g., in the subtropical gyres) in the tropical ocean, with a marked exception being the equatorial Pacific, which shows a significant increase in N2 fixation. This is likely due to climate driven changes in upwelling, which reduces N supply therein. Any change to the magnitude of “new” N being supplied by N2 fixation in the tropics will serve to amplify or dampen any stratification driven reduction in N supply. N2 fixation trends emerge later than trends in both sea surface temperature and pH (2016–2035), but earlier than trends in both oxygen and NPP (>2095, Frölicher et al., 2016; Henson et al., 2017; Keller et al., 2014; Rodgers et al., 2015).
Although climate change will directly impact NPP, it is important to also investigate the processes that indirectly affect NPP, such as N2 fixation and how it may modulate the NPP signal. For instance, in regions where N2 fixation declines (western and subtropical N. Pacific, Indian Ocean and the subtropical N. Atlantic) there is a strong contribution (between 4 and 73%) from the decrease in N2 fixation to the overall decrease of the upper ocean nitrate (NO3) inventory (Figures 3A,D). On the other hand, areas of increased N2 fixation act against the loss of NO3 from the upper ocean due to enhanced stratification (Figures 3A,D). In Figure 3D, a positive ratio indicates that the N2 fixation trend is in the same direction as the overall upper NO3 inventory trend, whereas a negative ratio indicates opposing trends (Figure 3A shows the overall directional trend in N2 fixation). As NPP in the topics is largely responding to changes to upper ocean N, the impact of climate driven changes in stratification may be amplified or dampened by the more rapidly emerging trends in rates of N2 fixation. Future work should also focus on evaluating the processes underpinning this potentially important feedback (e.g., rates of nitrification, zooplankton recycling and the sinking and remineralisation of diazotroph organic material) within the nitrogen cycle in Earth System Models.
Mechanistic Limitations of Current Models
Our analysis has shown that decadal trends in N2 fixation in response to climate change are likely to contribute to the changes in NPP simulated in Earth System Models. But we have found that a range of complexity exists for how N2 fixation is parameterized within the current models and large disagreement occurs in the spatial pattern and trends of N2 fixation projected by the different models. However, despite these differences, all models show poor skill compared to compiled N2 fixation observations (Table 1). We took the observations from the nitrogen fixation database of Tang et al. (2019b), regridded them onto the model horizontal grid and compared them to the model fields averaged over 1985:2005. Overall, model skill, taken as the correlation coefficient, ranges from as low as 0.06 for GFDL-ESM2M or up to 0.35 for MPI-ESM-LR (Table 1). This suggests that the change in N2 fixation in response to climate change and the associated impact on nitrogen supply to primary producers in the oligotrophic ocean remains a very poorly constrained component of Earth System Models. The overall trend in N2 fixation across the majority of the ESMs is at odds with other perturbation experiments predicting the opposite (Boyd and Doney, 2002; Dutkiewicz et al., 2014; Jiang et al., 2018; Tang et al., 2019a), highlighting the need for more consideration of the competition between diazotrophs and non-diazotrophs within ESMs in the future.
There are a range of potentially significant drivers not yet included in the parameterization of ocean N2 fixation in Earth System Models, which could affect the response to climate change and contribute to the general lack of model skill. Firstly, increasing CO2 has been shown to produce a “nutrient-like” response for growth and N2 fixation rates for both Trichodesmium and Crocosphaera spp. (Hutchins et al., 2013; Hutchins and Fu, 2017), with CO2 potentially defining an upper limit for maximum N2 fixation rates (Hutchins et al., 2013). However, increasing CO2 will increase ocean acidification, which has been shown to negatively impact diazotrophy by reducing the efficiency of the nitrogenase enzyme (Hong et al., 2017). Secondly, warming may push tropical diazotrophs past their thermal optima and/or facilitate the expansion of recently observed cold water diazotrophs occurring below 20°C (Needoba et al., 2007) and even in polar waters (Harding et al., 2018; Shiozaki et al., 2018). Breitbarth et al. (2007) found that the maximum rate of growth and N2 fixation in Trichodesmium occurs between 24 and 30°C and both Trichodesmium and Crocosphaera growth has been shown to rapidly decline at temperatures above 28 and 30°C, respectively (Boyd et al., 2013). Most models do not include an upper thermal threshold and so depending on how growth is parameterized, rates of diazotrophy could continue to increase past observed thermal thresholds. If the ecological niche for diazotrophy moves poleward due to excessive warming, light limitation may become a more important factor than previously (Boyd and Doney, 2002; Boatman et al., 2017). Thirdly, laboratory studies show that when Fe and P are both present at low concentrations, Trichodesmium cell sizes decreased, allowing uptake rates to be maintained and promoting both growth and diazotrophy, suggesting that Fe-P co-limitation may be as important in shaping N2 fixation rates (Walworth et al., 2016, 2018). However, dust deposition is predicted to increase due to enhanced desertification (Mahowald et al., 2017). This may promote diazotrophy by supplying more Fe, but would enhance P limitation (Mills et al., 2004) and negatively impact the effect of Fe:P co-limitation.
Previously it was believed that the contribution of coastal regions to global marine N2 fixation was very small, but over the last decade the importance of coastal diazotrophy has been highlighted, with it being estimated to contribute 16.7 ± 14.3 Tg N yr–1 (Mulholland et al., 2012; Tang et al., 2019b). Due to the coarse resolution of climate models, coastal environments are poorly represented and will need to be improved in the future to gain a more complete understanding of global marine N2 fixation. The majority of models have an explicit parameterization of diazotrophy use only one diazotroph for the entire ocean, but actual diazotroph biodiversity is far more complex (Tang et al., 2019b). Similarly, the reliance of models on the colonial cyanobacterium Trichodesmium as a model organism neglects other important groups such as the UCYN-A, which may be predominant globally (Martinez-Perez et al., 2016), and heterotrophic diazotrophs which occupy niches previously believed to be unfavorable to diazotrophs (Rahav et al., 2013; Benavides et al., 2016). Model studies tend to focus upon the impact of bottom-up control on diazotrophs and generally overlook top-down controls. For instance, Wang et al. (2019) showed that when the zooplankton maximum grazing rate on diazotrophs was decreased to that of diatoms, global N2 fixation increased by 62%. Diazotroph growth rates are also not well constrained, and models are hence biased toward relying on observed growth rates of Trichodesmium, which will restrict the modeled diazotroph to their niche (Dutkiewicz et al., 2009). In the future, a new generation of N2 fixation models that include these multiple interacting effects and account for diazotroph biodiversity would help increase confidence in our projections of how climate change will affect the ocean N cycle (Tang and Cassar, 2019) and net primary production in the tropics.
Conclusion
Nitrogen fixation will be crucial to supporting productivity in the future tropical ocean that is dustier and more stratified. Current models project different spatial patterns and magnitudes of diazotrophy in the contemporary ocean, models, however, broadly fall into two groups based upon their spatial distribution of N2 fixation. Generally models project that globally diazotrophy will decrease in the future, at the regional scale the models reveal more complexity in the response to climate change with the Pacific counteracting the decline projected for the Atlantic and Indian Oceans. Gyre models exhibited more coherency in the response of N2 fixation at the regional scale than upwelling models. For the models as a whole, the strongest inter-model agreement coincided with regions of strong decreases in diazotrophy in the Atlantic and Indian Oceans, whilst low model agreement occurred with strong increases in diazotrophy in the Pacific. The change in new nitrogen supply from diazotrophy may alleviate any changes to N limitation in the oligotrophic ocean due to increased stratification. The rapid emergence of the climate driven trend of diazotrophy from background variability by around 2040, implies that changes in new nitrogen supply from diazotrophy is likely to contribute to the modulation of the later emerging trends in net primary productivity on a regional scale. Overall, there is a need for a new generation of N2 fixation Earth System Models that account for the multiple interacting factors which drive diazotrophy including both top-down and bottom-up control. Understanding how climate will impact diazotrophy is essential to gain insight into how patterns of nutrient limitation and net primary productivity are likely to respond to future climate perturbations.
Data Availability Statement
The model data that was analyzed in this study was generated as part of CMIP5 and is freely available from the Earth System Grid Federation (https://esgf-node.llnl.gov/projects/cmip5/).
Author Contributions
LW and AT designed the study and wrote the manuscript. LW performed the data analysis. Both authors contributed to the article and approved the submitted version.
Funding
This work contained in this publication was conducted during a Ph.D. study supported by the Natural Environment Research Council (NERC) EAO Doctoral Training Partnership and is fully-funded by NERC (NE/L002469/1).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
The authors would like to thank Thomas Frölicher for discussion on the time of emergence calculation, Laurent Bopp, and Lester Kwiatkowski for providing model output, as well as all CMIP5 participants for making model output freely available. The authors also thank two reviewers for their comments that improved the manuscript.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmars.2020.00635/full#supplementary-material
References
Aumont, O., and Bopp, L. (2006). Globalizing results from ocean in situ iron fertilization studies. Global Biogeochem. Cycles 20:GB2017. doi: 10.1029/2005gb002591
Benavides, M., Bonnet, S., Hernández, N., Martínez-Pérez, A. M., Nieto-Cid, M., Álvarez-Salgado, X. A., et al. (2016). Basin-wide N2 fixation in the deep waters of the Mediterranean Sea. Global Biogeochem. Cycles 30, 952–961. doi: 10.1002/2015gb005326
Boatman, T. G., Lawson, T., and Geider, R. J. (2017). A key marine diazotroph in a changing ocean: the interacting effects of temperature, CO2 and light on the growth of Trichodesmium erythraeum IMS101. PLoS One 12:e0168796. doi: 10.1371/journal.pone.0168796
Bombar, D., Paerl, R. W., and Riemann, L. (2016). Marine non-cyanobacterial diazotrophs: moving beyond molecular detection. Trends Microbiol. 24, 916–927. doi: 10.1016/j.tim.2016.07.002
Bopp, L., Resplandy, L., Orr, J. C., Doney, S. C., Dunne, J. P., Gehlen, M., et al. (2013). Multiple stressors of ocean ecosystems in the 21st century: projections with CMIP5 models. Biogeosciences 10, 6225–6245. doi: 10.5194/bg-10-6225-2013
Boyd, P. W., and Doney, S. C. (2002). Modelling regional responses by marine pelagic ecosystems to global climate change. Geophys. Res. Lett. 29, 53–51. doi: 10.1029/2001gl014130
Boyd, P. W., Rynearson, T. A., Armstrong, E. A., Fu, F., Hayashi, K., Hu, Z., et al. (2013). Marine phytoplankton temperature versus growth responses from polar to tropical waters–outcome of a scientific community-wide study. PLoS One 8:e63091. doi: 10.1371/journal.pone.0063091
Breitbarth, E., Oschlies, A., and LaRoche, J. (2007). Physiological constraints on the global distribution of Trichodesmium? effect of temperature on diazotrophy. Biogeosciences 4, 53–61. doi: 10.5194/bg-4-53-2007
Ciais, P., Sabine, C., Bala, G., Bopp, L., Brovkin, V., Canadell, J., et al. (2014). “Carbon and other biogeochemical cycles,” in Climate change 2013: the Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, (Cambridge, MA: Cambridge University Press), 465–570.
Dunne, J. P., John, J. G., Shevliakova, E., Stouffer, R. J., Krasting, J. P., Malyshev, S. L., et al. (2013). GFDL’s ESM2 global coupled climate–carbon earth system models. Part II: carbon system formulation and baseline simulation characteristics. J. Clim. 26, 2247–2267. doi: 10.1175/JCLI-D-12-00150.1
Dutkiewicz, S., Follows, M. J., and Bragg, J. G. (2009). Modeling the coupling of ocean ecology and biogeochemistry. Global Biogeochem. Cycles 23, 1–15. doi: 10.1029/2008gb003405
Dutkiewicz, S., Ward, B. A., Monteiro, F., and Follows, M. J. (2012). Interconnection of nitrogen fixers and iron in the Pacific Ocean: theory and numerical simulations. Global Biogeochem. Cycles 26, 1–16. doi: 10.1029/2011GB004039
Dutkiewicz, S., Ward, B. A., Scott, J. R., and Follows, M. J. (2014). Understanding predicted shifts in diazotroph biogeography using resource competition theory. Biogeosciences 11, 5445–5461. doi: 10.5194/bg-11-5445-2014
Falkowski, P. G. (1997). Evolution of the nitrogen cycle and its influence on the biological sequestration of CO2 in the ocean. Nature 387, 272–275. doi: 10.1038/387272a0
Frölicher, T. L., Rodgers, K. B., Stock, C. A., and Cheung, W. W. L. (2016). Sources of uncertainties in 21st century projections of potential ocean ecosystem stressors. Global Biogeochem. Cycles 30, 1224–1243. doi: 10.1002/2015GB005338
Goebel, N. L., Edwards, C. A., Carter, B. J., Achilles, K. M., and Zehr, J. P. (2008). Growth and carbon content of three different-sized diazotrophic cyanobacteria observed in the subtropical North Pacific. J. Phycol. 44, 1212–1220. doi: 10.1111/j.1529-8817.2008.00581.x
Großkopf, T., and Laroche, J. (2012). Direct and indirect costs of dinitrogen fixation in Crocosphaera watsonii WH8501 and possible implications for the nitrogen cycle. Front. Microbiol. 3:236. doi: 10.3389/fmicb.2012.00236
Gruber, N., and Sarmiento, J. L. (1997). Global patterns of marine nitrogen fixation and denitrification. Global Biogeochem. Cycles 11, 235–266. doi: 10.1029/97GB00077
Harding, K., Turk-Kubo, K. A., Sipler, R. E., Mills, M. M., Bronk, D. A., and Zehr, J. P. (2018). Symbiotic unicellular cyanobacteria fix nitrogen in the Arctic Ocean. Proc. Natl. Acad. Sci. U.S.A. 115, 13371–13375. doi: 10.1073/pnas.1813658115
Henson, S. A., Beaulieu, C., Ilyina, T., John, J. G., Long, M., Séférian, R., et al. (2017). Rapid emergence of climate change in environmental drivers of marine ecosystems. Nat. Commun. 8:14682. doi: 10.1038/ncomms14682
Hong, H., Shen, R., Zhang, F., Wen, Z., Chang, S., Lin, W., et al. (2017). The complex effects of ocean acidification on the prominent N2-fixing cyanobacterium Trichodesmium. Science 356, 527–531. doi: 10.1126/science.aal2981
Hutchins, D. A., and Fu, F.-X. (2017). Microorganisms and ocean global change. Nat. Microbiol. 2:17058. doi: 10.1038/nmicrobiol.2017.58
Hutchins, D. A., Fu, F.-X., Webb, E. A., Walworth, N. G., and Tagliabue, A. (2013). Taxon-specific response of marine nitrogen fixers to elevated carbon dioxide concentrations. Nat. Geosci. 6, 790–795. doi: 10.1038/ngeo1858
Hutchins, D. A., Fu, F.-X., Zhang, Y., Warner, M., Feng, Y., Portune, K., et al. (2007). CO2 control of Trichodesmium N2 fixation, photosynthesis, growth rates, and elemental ratios: implications for past, present, and future ocean biogeochemistry. Limnol. Oceanogr. 52, 1293–1304. doi: 10.4319/lo.2007.52.4.1293
Ilyina, T., Six, K. D., Segschneider, J., Maier-Reimer, E., Li, H., and Núñez-Riboni, I. (2013). Global ocean biogeochemistry model HAMOCC: model architecture and performance as component of the MPI-Earth system model in different CMIP5 experimental realizations. J. Adv. Model. Earth Syst. 5, 287–315. doi: 10.1029/2012MS000178
Jiang, H.-B., Fu, F.-X., Rivero-Calle, S., Levine, N. M., Sañudo-Wilhelmy, S. A., Qu, P.-P., et al. (2018). Ocean warming alleviates iron limitation of marine nitrogen fixation. Nat. Clim. Change 8, 709–712. doi: 10.1038/s41558-018-0216-8
Jickells, T. D., Buitenhuis, E., Altieri, K., Baker, A. R., Capone, D., Duce, R. A., et al. (2017). A reevaluation of the magnitude and impacts of anthropogenic atmospheric nitrogen inputs on the ocean. Global Biogeochem. Cycles 31, 289–305. doi: 10.1002/2016gb005586
Keller, K. M., Joos, F., and Raible, C. C. (2014). Time of emergence of trends in ocean biogeochemistry. Biogeosciences 11, 3647–3659. doi: 10.5194/bg-11-3647-2014
Kustka, A. B., Sañudo-Wilhelmy, S. A., Carpenter, E. J., Capone, D., Burns, J., and Sunda, W. G. (2003). Iron requirements for dinitrogen- and ammonium-supported growth in cultures of Trichodesmium (IMS 101): comparison with nitrogen fixation rates and iron: carbon ratios of field populations. Limnol. Oceanogr. 48, 1869–1884. doi: 10.4319/lo.2003.48.5.1869
Kwiatkowski, L., Aumont, O., Bopp, L., and Ciais, P. (2018). The impact of variable phytoplankton stoichiometry on projections of primary production, food quality, and Carbon Uptake in the Global Ocean. Global Biogeochem. Cycles 32, 516–528. doi: 10.1002/2017GB005799
Laufkötter, C., Vogt, M., Gruber, N., Aita-Noguchi, M., Aumont, O., Bopp, L., et al. (2015). Drivers and uncertainties of future global marine primary production in marine ecosystem models. Biogeosciences 12, 6955–6984. doi: 10.5194/bg-12-6955-2015
Luo, Y., Doney, S., Anderson, L., Benavides, M., Berman-Frank, I., Bode, A., et al. (2012). Database of diazotrophs in global ocean: abundance, biomass, and nitrogen fixation rates. Earth Syst. Sci. Data 4, 47–73. doi: 10.5194/essd-4-47-2012
Luo, Y. W., Lima, I. D., Karl, D. M., Deutsch, C. A., and Doney, S. C. (2014). Data-based assessment of environmental controls on global marine nitrogen fixation. Biogeosciences 11, 691–708. doi: 10.5194/bg-11-691-2014
Mahaffey, C., Michaels, A. F., and Capone, D. G. (2005). The conundrum of marine N2 fixation. Am. J. Sci. 305, 546–595. doi: 10.2475/ajs.305.6-8.546
Mahowald, N. M., Scanza, R., Brahney, J., Goodale, C. L., Hess, P. G., Moore, J. K., et al. (2017). Aerosol deposition Impacts on Land and Ocean Carbon Cycles. Curr. Clim. Change Rep. 3, 16–31. doi: 10.1007/s40641-017-0056-z
Martinez-Perez, C., Mohr, W., Löscher, C. R., Dekaezemacker, J., Littmann, S., Yilmaz, P., et al. (2016). The small unicellular diazotrophic symbiont, UCYN-A, is a key player in the marine nitrogen cycle. Nat. Microbiol. 1:16163. doi: 10.1038/nmicrobiol.2016.163
Mills, M. M., Ridame, C., Davey, M., La Roche, J., and Geider, R. J. (2004). Iron and phosphorus co-limit nitrogen fixation in the eastern tropical North Atlantic. Nature 429, 292–294. doi: 10.1038/nature02550
Moisander, P. H., Beinart, R. A., Hewson, I., White, A. E., Johnson, K. S., Carlson, C. A., et al. (2010). Unicellular cyanobacterial distributions broaden the oceanic N2 fixation domain. Science 327, 1512–1514. doi: 10.1126/science.1185468
Moisander, P. H., Benavides, M., Bonnet, S., Berman-Frank, I., White, A. E., and Riemann, L. (2017). Chasing after non-cyanobacterial nitrogen fixation in marine pelagic environments. Front. Microbiol. 8:1736. doi: 10.3389/fmicb.2017.01736
Moore, C. M., Mills, M. M., Achterberg, E. P., Geider, R. J., LaRoche, J., Lucas, M. I., et al. (2009). Large-scale distribution of Atlantic nitrogen fixation controlled by iron availability. Nat. Geosci. 2, 867–871. doi: 10.1038/ngeo667
Moore, C. M., Mills, M. M., Arrigo, K. R., Berman-Frank, I., Bopp, L., Boyd, P. W., et al. (2013). Processes and patterns of oceanic nutrient limitation. Nat. Geosci. 6, 701–710. doi: 10.1038/ngeo1765
Moore, J. K., Fu, W., Primeau, F., Britten, G. L., Lindsay, K., Long, M., et al. (2018). Sustained climate warming drives declining marine biological productivity. Science 359, 1139–1143. doi: 10.1126/science.aao6379
Moore, J. K., Lindsay, K., Doney, S. C., Long, M. C., and Misumi, K. (2013). Marine ecosystem dynamics and biogeochemical cycling in the community earth system model [CESM1(BGC)]: comparison of the 1990s with the 2090s under the RCP4.5 and RCP8.5 Scenarios. J. Clim. 26, 9291–9312. doi: 10.1175/jcli-d-12-00566.1
Mulholland, M. R., Bernhardt, P. W., Blanco-Garcia, J. L., Mannino, A., Hyde, K., Mondragon, E., et al. (2012). 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. doi: 10.4319/lo.2012.57.4.1067
Needoba, J. A., Foster, R. A., Sakamoto, C., Zehr, J. P., and Johnson, K. S. (2007). Nitrogen fixation by unicellular diazotrophic cyanobacteria in the temperate oligotrophic North Pacific Ocean. Limnol. Oceanogr. 52, 1317–1327. doi: 10.4319/lo.2007.52.4.1317
Rahav, E., Bar-Zeev, E., Ohayion, S., Elifantz, H., Belkin, N., Herut, B., et al. (2013). Dinitrogen fixation in aphotic oxygenated marine environments. Front. Microbiol. 4:227. doi: 10.3389/fmicb.2013.00227
Raven, J. A. (1988). The iron and molybdenum use efficiencies of plant growth with different energy, carbon and nitrogen sources. New Phytol. 109, 279–287. doi: 10.1111/j.1469-8137.1988.tb04196.x
Riche, O., and Christian, J. (2018). Ocean dinitrogen fixation and its potential effects on ocean primary production in Earth system model simulations of anthropogenic warming. Elem Sci. Anth. 6:16. doi: 10.1525/elementa.277
Rodgers, K. B., Lin, J., and Frölicher, T. L. (2015). Emergence of multiple ocean ecosystem drivers in a large ensemble suite with an Earth system model. Biogeosciences 12, 3301–3320. doi: 10.5194/bg-12-3301-2015
Sañudo-Wilhelmy, S. A., Kustka, A. B., Gobler, C. J., Hutchins, D. A., Yang, M., Lwiza, K., et al. (2001). Phosphorus limitation of nitrogen fixation by Trichodesmium in the central Atlantic Ocean. Nature 411, 66–69. doi: 10.1038/35075041
Shiozaki, T., Fujiwara, A., Ijichi, M., Harada, N., Nishino, S., Nishi, S., et al. (2018). Diazotroph community structure and the role of nitrogen fixation in the nitrogen cycle in the Chukchi Sea (western Arctic Ocean). Limnol. Oceanogr. 63, 2191–2205. doi: 10.1002/lno.10933
Sohm, J. A., Mahaffey, C., and Capone, D. G. (2008). Assessment of relative phosphorus limitation of Trichodesmium spp. in the North Pacific, North Atlantic, and the north coast of Australia. Limnol. Oceanogr. 53, 2495–2502. doi: 10.4319/lo.2008.53.6.2495
Sohm, J. A., Webb, E. A., and Capone, D. G. (2011). Emerging patterns of marine nitrogen fixation. Nat. Rev. Microbiol. 9, 499–508. doi: 10.1038/nrmicro2594
Steinacher, M., Joos, F., Frölicher, T., Bopp, L., Cadule, P., Cocco, V., et al. (2010). Projected 21st century decrease in marine productivity: a multi-model analysis. Biogeosciences 7, 979–1005. doi: 10.5194/bg-7-979-2010
Tang, W., and Cassar, N. (2019). Data-driven modeling of the distribution of diazotrophs in the global ocean. Geophys. Res. Lett. 46, 12258–12269. doi: 10.1029/2019gl084376
Tang, W., Li, Z., and Cassar, N. (2019a). Machine learning estimates of global marine nitrogen fixation. J. Geophys. Res. Biogeosci. 124, 717–730. doi: 10.1029/2018jg004828
Tang, W., Wang, S., Fonseca-Batista, D., Dehairs, F., Gifford, S., Gonzalez, A. G., et al. (2019b). Revisiting the distribution of oceanic N2 fixation and estimating diazotrophic contribution to marine production. Nat. Commun. 10:831. doi: 10.1038/s41467-019-08640-0
Walworth, N. G., Fu, F.-X., Webb, E. A., Saito, M. A., Moran, D., Mcllvin, M. R., et al. (2016). Mechanisms of increased Trichodesmium fitness under iron and phosphorus co-limitation in the present and future ocean. Nat. Commun. 7:12081. doi: 10.1038/ncomms12081
Walworth, N. G., Fu, F.-X., Lee, M. D., Cai, X., Saito, M. A., Webb, E. A., et al. (2018). Nutrient-colimited Trichodesmium as a nitrogen source or sink in a future ocean. Appl. Environ. Microbiol. 84:e02137-17. doi: 10.1128/aem.02137-17
Wang, W.-L., Moore, J. K., Martiny, A. C., and Primeau, F. W. (2019). Convergent estimates of marine nitrogen fixation. Nature 566, 205–211. doi: 10.1038/s41586-019-0911-2
Keywords: diazotrophy, climate change, biogeochemistry, earth system models, nitrogen cycle
Citation: Wrightson L and Tagliabue A (2020) Quantifying the Impact of Climate Change on Marine Diazotrophy: Insights From Earth System Models. Front. Mar. Sci. 7:635. doi: 10.3389/fmars.2020.00635
Received: 27 January 2020; Accepted: 10 July 2020;
Published: 29 July 2020.
Edited by:
Wolfgang Koeve, GEOMAR Helmholtz Centre for Ocean Research Kiel, GermanyReviewed by:
Weiyi Tang, Princeton University, United StatesMar Benavides, Institut de Recherche Pour le Développement (IRD), France
Copyright © 2020 Wrightson and Tagliabue. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Lewis Wrightson, l.wrightson@liverpool.ac.uk; Alessandro Tagliabue, A.Tagliabue@liverpool.ac.uk