ALBERT

Description: There is currently no consensus on how humans are affecting the marine nitrogen (N) cycle, which limits marine biological production and CO2 uptake. Anthropogenic changes in ocean warming, deoxygenation, and atmospheric N deposition can all individually affect the marine N cycle and the oceanic production of the greenhouse gas nitrous oxide (N2O). However, the combined effect of these perturbations on marine N cycling, ocean productivity, and marine N2O production is poorly understood. Here we use an Earth system model of intermediate complexity to investigate the combined effects of estimated 21st century CO2 atmospheric forcing and atmospheric N deposition. Our simulations suggest that anthropogenic perturbations cause only a small imbalance to the N cycle relative to preindustrial conditions (∼+5 Tg N y−1 in 2100). More N loss from water column denitrification in expanded oxygen minimum zones (OMZs) is counteracted by less benthic denitrification, due to the stratification-induced reduction in organic matter export. The larger atmospheric N load is offset by reduced N inputs by marine N2 fixation. Our model predicts a decline in oceanic N2O emissions by 2100. This is induced by the decrease in organic matter export and associated N2O production and by the anthropogenically driven changes in ocean circulation and atmospheric N2O concentrations. After comprehensively accounting for a series of complex physical-biogeochemical interactions, this study suggests that N flux imbalances are limited by biogeochemical feedbacks that help stabilize the marine N inventory against anthropogenic changes. These findings support the hypothesis that strong negative feedbacks regulate the marine N inventory on centennial time scales.

Description: Nitrogen (N) fixation by specialized microorganisms (diazotrophs) influences global plankton productivity because it provides the ocean with most of its bio-available N. However, its global rate and large-scale spatial distribution is still regarded with considerable uncertainty. Here we use a global ocean nitrogen isotope model, in comparison with δ15NO3− observations, to constrain the pattern of N2 fixation across the Pacific Ocean. N2 fixation introduces isotopically light atmospheric N2 from to the ocean (δ15N = 0‰) relative to the oceanic average near 5‰, which makes nitrogen isotopes suitable to infer patterns of N2 fixation. Including atmospheric iron limitation of diazotrophy in the model shifts the pattern of simulated N2 fixation from the South Pacific to the North Pacific and from the East Pacific westward. These changes considerably improve the agreement with meridional transects of available δ15NO3− observations, as well as excess P (PO43− − NO3−/16), suggesting that atmospheric iron deposition is indeed important for N fixation in the Pacific Ocean. This study highlights the potential for using δ15N observations and model simulations to constrain patterns and rates of N fixation in the ocean.

Description: Dissolved organic nitrogen (DON) and phosphorus (DOP) represent the most abundant form of their respective nutrient pool in the surface layer of the oligotrophic oceans and play an important role in nutrient cycling and productivity. Since DOP is generally more labile than DON, it provides additional P that may stimulate growth of nitrogen-fixing diazotrophs that supply fixed nitrogen to balance denitrification in the ocean. In this study, we introduce semirecalcitrant components of DON and DOP as state variables in an existing global ocean-atmosphere-sea ice-biogeochemistry model of intermediate complexity to assess their impact on the spatial distribution of nitrogen fixation and the size of the marine fixed nitrogen inventory. Large-scale surface data sets of global DON and Atlantic Ocean DOP are used to constrain the model. Our simulations suggest that both preferential DOP remineralization and phytoplankton DOP uptake are important "non-Redfield" processes (i.e., deviate from molar N:P=16) that need to be accounted for to explain the observed patterns of DOP. Additional non-Redfield DOP sensitivity experiments testing dissolved organic matter (DOM) production rate uncertainties that best reproduce the observed spatial patterns of DON and DOP stimulate additional nitrogen fixation that increases the size of the global marine fixed nitrogen inventory by 4.7±1.7% compared to the simulation assuming Redfield DOM stoichiometry that underestimates the observed nitrogen inventory. The extra 8Tgyr-1 of nitrogen fixation stimulated in the Atlantic Ocean is mainly responsible for this increase due to its large spatial separation from water column denitrification, which buffers any potential nitrogen surplus in the Pacific Ocean. Our study suggests that the marine fixed nitrogen budget is sensitive to non-Redfield DOP dynamics because access to the relatively labile DOP pool expands the ecological niche for nitrogen-fixing diazotrophs.

Description: We present a new nitrogen isotope model incorporated into the three-dimensional ocean component of a global Earth system climate model designed for millennial timescale simulations. The model includes prognostic tracers for the two stable nitrogen isotopes, 14N and 15N, in the nitrate (NO3−), phytoplankton, zooplankton, and detritus variables of the marine ecosystem model. The isotope effects of algal NO3− uptake, nitrogen fixation, water column denitrification, and zooplankton excretion are considered as well as the removal of NO3− by sedimentary denitrification. A global database of δ15NO3− observations is compiled from previous studies and compared to the model results on a regional basis where sufficient observations exist. The model is able to qualitatively and quantitatively reproduce many of the observed patterns such as high subsurface values in water column denitrification zones and the meridional and vertical gradients in the Southern Ocean. The observed pronounced subsurface minimum in the Atlantic is underestimated by the model presumably owing to too little simulated nitrogen fixation there. Sensitivity experiments reveal that algal NO3− uptake, nitrogen fixation, and water column denitrification have the strongest effects on the simulated distribution of nitrogen isotopes, whereas the effect from zooplankton excretion is weaker. Both water column and sedimentary denitrification also have important indirect effects on the nitrogen isotope distribution by reducing the fixed nitrogen inventory, which creates an ecological niche for nitrogen fixers and, thus, stimulates additional N2 fixation in the model. Important model deficiencies are identified, and strategies for future improvement and possibilities for model application are outlined.

Description: The impact of increasing anthropogenic atmospheric nitrogen deposition on marine biogeochemistry is uncertain. We performed simulations to quantify its effect on nitrogen cycling and marine productivity in a global 3-D ocean biogeochemistry model. Nitrogen fixation provides an efficient feedback by decreasing immediately to deposition, whereas water column denitrification increases more gradually in the slowly expanding oxygen deficient zones. Counterintuitively, nitrogen deposition near oxygen deficient zones causes a net loss of marine nitrogen due to the stoichiometry of denitrification. In our idealized atmospheric deposition simulations that only account for nitrogen cycle perturbations, these combined stabilizing feedbacks largely compensate deposition and suppress the increase in global marine productivity to 〈2%, in contrast to a simulation that neglects nitrogen cycle feedbacks that predicts an increase of 〉15%. Our study emphasizes including the dynamic response of nitrogen fixation and denitrification to atmospheric nitrogen deposition to predict future changes of the marine nitrogen cycle and productivity.

Description: Key Points: Use of sedimentary nitrogen isotopes is examined; On average, sediment 15N/14N increases approx. 2 per mil during early burial; Isotopic alteration scales with water depth Abstract: Nitrogen isotopes are an important tool for evaluating past biogeochemical cycling from the paleoceanographic record. However, bulk sedimentary nitrogen isotope ratios, which can be determined routinely and at minimal cost, may be altered during burial and early sedimentary diagenesis, particularly outside of continental margin settings. The causes and detailed mechanisms of isotopic alteration are still under investigation. Case studies of the Mediterranean and South China Seas underscore the complexities of investigating isotopic alteration. In an effort to evaluate the evidence for alteration of the sedimentary N isotopic signal and try to quantify the net effect, we have compiled and compared data demonstrating alteration from the published literature. A 〉100 point comparison of sediment trap and surface sedimentary nitrogen isotope values demonstrates that, at sites located off of the continental margins, an increase in sediment 15N/14N occurs during early burial, likely at the seafloor. The extent of isotopic alteration appears to be a function of water depth. Depth-related differences in oxygen exposure time at the seafloor are likely the dominant control on the extent of N isotopic alteration. Moreover, the compiled data suggest that the degree of alteration is likely to be uniform through time at most sites so that bulk sedimentary isotope records likely provide a good means for evaluating relative changes in the global N cycle.

Description: A three-dimensional, process-based model of the ocean's carbon and nitrogen cycles, including 13C and 15N isotopes, is used to explore effects of idealized changes in the soft-tissue biological pump. Results are presented from one preindustrial control run (piCtrl) and six simulations of the Last Glacial Maximum (LGM) with increasing values of the spatially constant maximum phytoplankton growth rate μmax, which accelerates biological nutrient utilization mimicking iron fertilization. The default LGM simulation, without increasing μmax and with a shallower and weaker Atlantic Meridional Overturning Circulation and increased sea ice cover, leads to 280 Pg more respired organic carbon (Corg) storage in the deep ocean with respect to piCtrl. Dissolved oxygen concentrations in the colder glacial thermocline increase, which reduces water column denitrification and, with delay, nitrogen fixation, thus increasing the ocean's fixed nitrogen inventory and decreasing δ15NNO3 almost everywhere. This simulation already fits sediment reconstructions of carbon and nitrogen isotopes relatively well, but it overestimates deep ocean δ13CDIC and underestimates δ15NNO3 at high latitudes. Increasing μmax enhances Corg and lowers deep ocean δ13CDIC, improving the agreement with sediment data. In the model's Antarctic and North Pacific Oceans modest increases in μmax result in higher δ15NNO3 due to enhanced local nutrient utilization, improving the agreement with reconstructions there. Models with moderately increased μmax fit both isotope data best, whereas large increases in nutrient utilization are inconsistent with nitrogen isotopes although they still fit the carbon isotopes reasonably well. The best fitting models reproduce major features of the glacial δ13CDIC, δ15N, and oxygen reconstructions while simulating increased Corg by 510–670 Pg compared with the preindustrial ocean. These results are consistent with the idea that the soft-tissue pump was more efficient during the LGM. Both circulation and biological nutrient utilization could contribute. However, these conclusions are preliminary given our idealized experiments, which do not consider changes in benthic denitrification and spatially inhomogenous changes in aeolian iron fluxes. The analysis illustrates interactions between the carbon and nitrogen cycles as well as the complementary constraints provided by their isotopes. Whereas carbon isotopes are sensitive to circulation changes and indicate well the three-dimensional Corg distribution, nitrogen isotopes are more sensitive to biological nutrient utilization.

Description: The nitrogen cycle is fundamental to Earth's biogeochemistry. Yet major uncertainties of quantification remain, particularly regarding the global oceanic nitrogen fixation rate. Hydrogen is produced during nitrogen fixation and will become supersaturated in surface waters if there is net release from diazotrophs. Ocean surveys of hydrogen supersaturation thus have the potential to illustrate the spatial and temporal distribution of nitrogen fixation, and to guide the far more onerous but quantitative methods for measuring it. Here we present the first transect of high resolution measurements of hydrogen supersaturations in surface waters along a meridional 10,000 km cruise track through the Atlantic. We compare measured saturations with published measurements of nitrogen fixation rates and also with model-derived values. If the primary source of excess hydrogen is nitrogen fixation and has a hydrogen release ratio similar to Trichodesmium, our hydrogen measurements would point to similar rates of fixation in the North and South Atlantic, roughly consistent with modelled fixation rates but not with measured rates, which are lower in the south. Possible explanations would include any substantial nitrogen fixation by newly discovered diazotrophs, particularly any having a hydrogen release ratio similar to or exceeding that of Trichodesmium; under-sampling of nitrogen fixation south of the equator related to excessive focus on Trichodesmium; and methodological shortcomings of nitrogen fixation techniques that cause a bias towards colonial diazotrophs relative to unicellular forms. Alternatively our data are affected by an unknown hydrogen source that is greater in the southern half of the cruise track than the northern.

Description: Changes in the ocean iron cycle could help explain the low atmospheric CO2 during the Last Glacial Maximum (LGM). Previous modeling studies have mostly considered changes in aeolian iron fluxes, although it is known that sedimentary and hydrothermal fluxes are important iron sources for today's ocean. Here we explore effects of preindustrial‐to‐LGM changes in atmospheric dust, sedimentary, and hydrothermal fluxes on the ocean's iron and carbon cycles in a global coupled biogeochemical‐circulation model. Considering variable atmospheric iron solubility decreases LGM surface soluble iron fluxes compared with assuming constant solubility. This limits potential increases in productivity and export production due to surface iron fertilization, lowering atmospheric CO2 by only 4 ppm. The effect is countered by a decrease in sedimentary flux due to lower sea level, which increases CO2 by 15 ppm. Assuming a 10 times higher iron dust solubility in the Southern Ocean, combined with changes in sedimentary flux, we obtain an atmospheric CO2 reduction of 13 ppm. The high uncertainty in the iron fluxes does not allow us to determine the net direction and magnitude of variations in atmospheric CO2 due to changes in the iron cycle. Our model does not account for changes to iron‐binding ligand concentrations that could modify the results. We conclude that when evaluating glacial‐interglacial changes in the ocean iron cycle, not only surface but also seafloor fluxes must be taken into account.

Description: Changes in the ocean iron cycle could help explain the low atmospheric CO2 during the Last Glacial Maximum (LGM). Previous modeling studies have mostly considered changes in aeolian iron fluxes, although it is known that sedimentary and hydrothermal fluxes are important iron sources for today's ocean. Here we explore effects of preindustrial-to-LGM changes in atmospheric dust, sedimentary, and hydrothermal fluxes on the ocean's iron and carbon cycles in a global coupled biogeochemical-circulation model. Considering variable atmospheric iron solubility decreases LGM surface soluble iron fluxes compared with assuming constant solubility. This limits potential increases in productivity and export production due to surface iron fertilization, lowering atmospheric CO2 by only 4 ppm. The effect is countered by a decrease in sedimentary flux due to lower sea level, which increases CO2 by 15 ppm. Assuming a 10 times higher iron dust solubility in the Southern Ocean, combined with changes in sedimentary flux, we obtain an atmospheric CO2 reduction of 13 ppm. The high uncertainty in the iron fluxes does not allow us to determine the net direction and magnitude of variations in atmospheric CO2 due to changes in the iron cycle. Our model does not account for changes to iron-binding ligand concentrations that could modify the results. We conclude that when evaluating glacial-interglacial changes in the ocean iron cycle, not only surface but also seafloor fluxes must be taken into account.