Nitrogen Cycle Model Research Articles

Ammonium “nutrient capacitor” model for δ15N signatures associated with marine anoxic events

Abstract Geochemical records of ancient periods of warm climate can be useful to help understand the looming effects of modern anthropogenic warming, including changes to biogeochemical nutrient cycles. Stable nitrogen isotope compositions of marine sediments archive the balance of processes in the global nitrogen cycle. However, the unusual isotopic signals of Mesozoic oceanic anoxic events (OAEs) remain enigmatic, thus hindering our understanding of nitrogen cycle processes and dynamics under conditions of ocean deoxygenation. Here, we present an ammonium “nutrient capacitor” model of the water-column nitrogen cycle to explain the anomalously negative isotopic compositions seen in Mesozoic OAE sediments. Our model applies isotopic inferences derived from high-resolution records of Lake Kivu sediments to show how periodic chemocline overturning of redox-stratified water columns during Mesozoic OAEs may have delivered ammonium to the photic zone in excess of primary producer requirements. Smoothed, stochastic sampling of the changing fluxes within the nitrogen cycle across these events can simulate OAE nitrogen isotope records.

A century of nitrogen dynamics in agricultural watersheds of Denmark

Intensive agriculture has been linked to increased nitrogen loads and adverse effects on downstream aquatic ecosystems. Sustained large net nitrogen surpluses have been shown in several contexts to form legacies in soil or waters, which delay the effects of reduction measures. In this study, detailed land use and agricultural statistics were used to reconstruct the annual nitrogen surpluses in three agriculture-dominated watersheds of Denmark (600–2700 km2) with well-drained loamy soils. These surpluses and long-term hydrological records were used as inputs to the process model ELEMeNT to quantify the nitrogen stores and fluxes for 1920–2020. A multi-objective calibration using timeseries of river nitrate loads, as well as other non-conventional data sources, allowed to explore the potential of these different data to constrain the nitrogen cycling model. We found the flux-weighted nitrate concentrations in the root zone percolate below croplands, a dataset not commonly used in calibrating watershed models, to be critical in reducing parameter uncertainty. Groundwater nitrate legacies built up in all three studied watersheds during 1950–1990 corresponding to ∼2% of the surplus (or ∼1 kg N ha yr−1) before they went down at a similar rate during 1990–2015. Over the same periods active soil nitrogen legacies first accumulated by approximately 10% of the surplus (∼5 kg N ha yr−1), before undergoing a commensurate reduction. Both legacies appear to have been the drivers of hysteresis in the diffuse load at the catchments’ outlet and hindrances to reaching water quality goals. Results indicate that the low cropland surpluses enforced during 2008–2015 had a larger impact on the diffuse river loads than the European Union’s untargeted grass set-aside policy of 1993–2008. Collectively, the measures of 1990–2015 are estimated to have reset the diffuse load regimes of the watersheds back to the situation prevailing in the 1960s.

Quantifying co‐benefits of water quality policies: An integrated assessment model of land and nitrogen management

AbstractDue to the nature of nitrogen cycling, policies designed to address water quality concerns have the potential to provide benefits beyond the targeted water quality improvements. For example, actions to protect water quality by reducing nitrate leaching from agriculture also reduce emissions of nitrous oxide, a potent greenhouse gas. These positive effects, which are incidental to the regulation's intended target, are termed “co‐benefits.” To quantify the co‐benefits associated with reduced nitrate leaching, we integrate an economic model of farmer decision making with a model of terrestrial nitrogen cycling for the watershed surrounding Lake Mendota, Wisconsin, USA. Our modeling approach provides a framework that links air and water pollutants in an agri‐environmental system and offers a direction for future studies. Our model results highlight the finding that the co‐benefits from nitrous oxide abatement are substantial, and their inclusion increases the benefit–cost ratio of water quality policies. Consideration of these co‐benefits has the potential to reverse the conclusions of benefit–cost analysis in the assessment of current water quality policies.

Pathways of Nitrous Oxide Production in the Eastern Tropical South Pacific Oxygen Minimum Zone

AbstractOceanic emissions of nitrous oxide (N2O) account for roughly one‐third of all natural sources to the atmosphere. Hot‐spots of N2O outgassing occur over oxygen minimum zones (OMZs), where the presence of steep oxygen gradients surrounding anoxic waters leads to enhanced N2O production from both nitrification and denitrification. However, the relative contributions from these pathways to N2O production and outgassing in these regions remains poorly constrained, in part due to shared intermediary nitrogen tracers, and the tight coupling of denitrification sources and sinks. To shed light on this problem, we embed a new, mechanistic model of the OMZ nitrogen cycle within a three‐dimensional eddy‐resolving physical‐biogeochemical model of the Eastern Tropical South Pacific (ETSP), tracking contributions from remote advection, atmospheric exchange, and local nitrification and denitrification. The model indicates that net N2O production from denitrification is approximately one order of magnitude greater than nitrification within the ETSP OMZ. However, only ∼32% of denitrification‐derived N2O production ultimately outgasses to the atmosphere in this region (contributing ∼36% of the air‐sea N2O flux on an annual basis), while the remaining is exported out of the domain. Instead, remotely produced N2O advected into the OMZ region accounts for roughly half (∼57%) of the total N2O outgassing, with smaller contributions from nitrification (∼7%). Our results suggests that, together with enhanced production by denitrification, upwelling of remotely derived N2O contributes the most to N2O outgassing over the ETSP OMZ.

Towards an integrated phosphorus, carbon and nitrogen cycling model for topographically diverse grasslands

Contemporary science on how livestock influence nutrient cycling in grazing systems is limited, particularly in topographically complex (i.e., slopes and aspects) hill country landscapes. Prominent slope and aspect variation affects primary production, animal behaviour and nutrient return. Here, we embed recent scientific advancements in nutrient dynamics across complex landscapes to (1) set up a soil organic carbon (SOC) saturation function to an existing SOC and total soil phosphorus (TSP) model (Bilotto et al. J N Z Grassl 81:171–178, 2019. https://doi.org/10.33584/jnzg.2019.81.397), (2) include total soil nitrogen (TSN) dynamics, and (3) establish if the model (herein the Grass-NEXT model) can simulate the spatial and temporal changes of TSP, SOC and TSN in hill country. A long-term P fertiliser experiment with contrasting different P fertilisation levels and associated sheep stocking regimes (herein, ‘farmlets’) was used for model testing. The Grass-NEXT model predicted TSP and SOC stocks with strong accuracy and precision (model performance), and TSN with a moderate performance across farmlets [Concordance Correlation Coefficient (CCC), 0.75, 0.72 and 0.49, respectively]. Grass-NEXT model simulated TSP, SOC and TSN distribution with moderate/strong performance across slopes (CCC, 0.94, 0.80 and 0.70) and aspects (CCC, 0.83, 0.67 and 0.51). Consistent with observed data, modelled changes in TSP and TSN were greater on low slopes and eastern aspects, but no clear pattern was observed for SOC stocks. The Grass-NEXT model provides an intuitive research tool for exploring management options for increasing SOC and TSN, as well as an instrument for monitoring and reporting on nutrient dynamics in complex landscapes.

Interpretation of abnormally negative carbon isotope excursions in Eifelian–Givetian (Middle Devonian) sediments of South China and implications for paleo-environmental variation

Many negative carbon isotope excursions (NCIEs) have been employed to interpret global biogeochemical events in Earth history, but they are also influenced by local depositional settings and diagenetic alteration, thus making it difficult to link them to global events. Like the NCIEs in Middle Devonian sediments of North America, similar NCIEs were also observed in the age-equivalent sediments of South China, and therefore it is necessary to further address the logical significance of these intercontinentally comparable NCIEs. Here, we present major and trace elements as well as δ15N, δ13Ccarb, and δ13Corg records across the Middle Devonian Eifelian–Givetian transition in the intra-platform basin of South China. The changes in stable isotopes (δ13Ccarb, δ13Corg, δ18Ocarb, and δ15N) and their relationships may jointly indicate limited influence of diagenetic processes. Our abnormally negative δ13Ccarb values (−7 to −5‰), even lighter than the mantle input end-member value (∼ −5‰), cannot be easily explained by dramatic changes in the δ13C values of global dissolved inorganic carbon (DIC) reservoir. In a restricted platform basin, enhanced nutrient input may stimulate high paleo-productivity, resulting in the depletion of DIC concentrations in surface waters. CO2 is then replenished via invasion from the atmosphere, and intense evaporative conditions may promote the large isotopic fractionation of CO2 invasion owing to hydration reactions of aqueous CO2. Our NCIEs are most likely related to the local 13C-depleted values of platform DIC owing to CO2 invasion associated with large isotopic fractionation due to enhanced evaporation and primary productivity. The “productivity-anoxia feedback” mechanism can explain the coupled changes in our terrestrial clastics input, productivity, and redox proxies. By combining redox proxies and C and N isotopic records, we further propose a model for local carbon and nitrogen cycling in a restricted and evaporative intra-platform basin. Our study highlights local depositional controls on NCIEs, and these local factors should be considered before NCIEs are applied to understand the global carbon cycle perturbations during the Middle Devonian.

Dam-induced flow velocity decrease leads to the transition from heterotrophic to autotrophic system through modifying microbial food web dynamics

The impoundment of reservoirs changes the river from a riverine heterotrophic system to a lacustrine autotrophic system, which could be attributed to the shift of pelagic microbial food webs in response to the dam-induced disturbances. However, little is known about what is the key factor controlling this variation and how different underlying interactions affect the food web dynamics. This study investigated the effects of flow velocity and nutrient supply on microbial plankton using a microcosm experiment. The results showed that flow velocity decrease was the main factor inducing the detritus-based food web transformed to the autotroph-based food web, with heterotrophic bacteria and protozoan dominated at high velocity, whereas phytoplankton and metazoan were prevalent in the lentic environment. The lentic-acclimated genera, such as Chlorella sp., Mallomonas sp. and Microcystis sp., showed hysteresis after the velocity recovery, suggesting the potential of algae bloom in reservoirs and even downstream of dams. We further conducted a flow-velocity manipulating experiment and constructed a multi-trophic nitrogen cycling model to provide a mechanistic explanation for the microbial food web dynamics and the nitrogen transformation performances. As indicated in model prediction and sensitivity analysis, the abiotic and biotic variations were directly or indirectly controlled by nutrient utilization and predator-prey interactions. Quantification of these bottom-up and top-down forces revealed the buffer role of predators in mitigating the positive effects of nutrient availability on autotrophs at low velocity and on heterotrophs at high velocity, respectively. This study highlights the importance of mastering the whole information of different trophic levels, in order to better capture the complex microbial food web interactions and the consequent biogeochemical processes in river-reservoir systems.

Substantial carbon drawdown potential from enhanced rock weathering in the United Kingdom

Achieving national targets for net-zero carbon emissions will require atmospheric carbon dioxide removal strategies compatible with rising agricultural production. One possible method for delivering on these goals is enhanced rock weathering, which involves modifying soils with crushed silicate rocks, such as basalt. Here we use dynamic carbon budget modelling to assess the carbon dioxide removal potential and agricultural benefits of implementing enhanced rock weathering strategies across UK arable croplands. We find that enhanced rock weathering could deliver net carbon dioxide removal of 6–30 MtCO2 yr−1 for the United Kingdom by 2050, representing up to 45% of the atmospheric carbon removal required nationally to meet net-zero emissions. This suggests that enhanced rock weathering could play a crucial role in national climate mitigation strategies if it were to gain acceptance across national political, local community and farm scales. We show that it is feasible to eliminate the energy-demanding requirement for milling rocks to fine particle sizes. Co-benefits of enhanced rock weathering include substantial mitigation of nitrous oxide, the third most important greenhouse gas, widespread reversal of soil acidification and considerable cost savings from reduced fertilizer usage. Our analyses provide a guide for other nations to pursue their carbon dioxide removal ambitions and decarbonize agriculture—a key source of greenhouse gases. Enhancing rock weathering across UK croplands could deliver substantial atmospheric carbon dioxide removal alongside agricultural co-benefits, according to coupled climate–carbon–nitrogen cycle model simulations.

An intense precipitation event causes a temperate forested drainage network to shift from N2O source to sink

AbstractNitrous oxide (N2O) evasion from streams and rivers is a significant, yet highly uncertain, flux in nitrogen cycle models. Most global estimates of lotic N2O emission assume that evasion rates are proportional to inorganic nitrogen inputs to a stream or river. However, many field studies do not detect relationships between lotic N2O evasion and dissolved nitrogen concentration, highlighting the need for better understanding of process‐based controls on this flux. This study reports 4‐yr time series of pN2O and N2O evasion from eight nested streams and rivers and detects an abrupt change in N2O dynamics associated with an intense rainstorm. This rainstorm, and the associated hydrologic flood event, pushed forested reaches across the watershed from consistent N2O sources to prolonged N2O sinks. We attribute this shift to disturbance of incomplete denitrification in the stream network and surrounding watershed, although alternate hypotheses are also discussed. There was continued availability of nitrate () for in‐stream processing, eliminating the possibility that ‐availability limited N2O production, and post‐storm N2O‐to‐nitrate ratios were lower than pre‐storm ratios suggesting that the large storm affected in‐situ nitrogen processing rates. The sustained period of post‐storm N2O undersaturation resulted in net negative evasion for five of the eight study sites in 2018, which mitigated emissions over the 4‐yr study. This nonlinear response in N2O dynamics illustrates the potential importance of storm events to control lotic N2O production and emissions.

Biogeochemical model of nitrogen cycling in Ahe (French Polynesia), a South Pacific coral atoll with pearl farming

A biogeochemical model (ECO3M-Atoll) was configured to simulate the lower food web in Ahe Atoll lagoon where phytoplankton is mostly nitrogen limited. Understanding the dynamics of phytoplankton – the main food source for oysters – is crucial for the management and the allocation of new pearl farming sites. After parametrizing the model with in situ observations, we tested different hypotheses about nitrogen cycling (benthic remineralization, atmospheric N fixation, etc.) and compared the results to a large observational dataset. Model results show that simulated (pico- and nano-) phytoplankton biomass and nitrogen concentrations are close to in situ data. The simulated biogeochemical processes (uptake and primary production) are also very similar to the observed values. In the model, primary production ranged from 1.00 to 2.00 mg C m−3 h−1 for pico- and 0.40 to 1.00 mg C m−3 h−1 for nanophytoplankton; mean N uptake was 2.02 μmol N m−3 h−1 for pico- and 1.25 μmol N m−3 h−1 for nanophytoplankton.

Implementation of nitrogen cycle in the CLASSIC land model

Abstract. A terrestrial nitrogen (N) cycle model is coupled to the carbon (C) cycle in the framework of the Canadian Land Surface Scheme Including Biogeochemical Cycles (CLASSIC). CLASSIC currently models physical and biogeochemical processes and simulates fluxes of water, energy, and CO2 at the land–atmosphere boundary. CLASSIC is similar to most models and its gross primary productivity increases in response to increasing atmospheric CO2 concentration. In the current model version, a downregulation parameterization emulates the effect of nutrient constraints and scales down potential photosynthesis rates, using a globally constant scalar, as a function of increasing CO2. In the new model when nitrogen (N) and carbon (C) cycles are coupled, cycling of N through the coupled soil–vegetation system facilitates the simulation of leaf N amount and maximum carboxylation capacity (Vcmax) prognostically. An increase in atmospheric CO2 decreases leaf N amount and therefore Vcmax, allowing the simulation of photosynthesis downregulation as a function of N supply. All primary N cycle processes that represent the coupled soil–vegetation system are modelled explicitly. These include biological N fixation; treatment of externally specified N deposition and fertilization application; uptake of N by plants; transfer of N to litter via litterfall; mineralization; immobilization; nitrification; denitrification; ammonia volatilization; leaching; and the gaseous fluxes of NO, N2O, and N2. The interactions between terrestrial C and N cycles are evaluated by perturbing the coupled soil–vegetation system in CLASSIC with one forcing at a time over the 1850–2017 historical period. These forcings include the increase in atmospheric CO2, change in climate, increase in N deposition, and increasing crop area and fertilizer input, over the historical period. An increase in atmospheric CO2 increases the C:N ratio of vegetation; climate warming over the historical period increases N mineralization and leads to a decrease in the vegetation C:N ratio; N deposition also decreases the vegetation C:N ratio. Finally, fertilizer input increases leaching, NH3 volatilization, and gaseous losses of N2, N2O, and NO. These model responses are consistent with conceptual understanding of the coupled C and N cycles. The simulated terrestrial carbon sink over the 1959–2017 period, from the simulation with all forcings, is 2.0 Pg C yr−1 and compares reasonably well with the quasi observation-based estimate from the 2019 Global Carbon Project (2.1 Pg C yr−1). The contribution of increasing CO2, climate change, and N deposition to carbon uptake by land over the historical period (1850–2017) is calculated to be 84 %, 2 %, and 14 %, respectively.

Global Dam‐Driven Changes to Riverine N:P:Si Ratios Delivered to the Coastal Ocean

AbstractRiver damming alters nutrient fluxes along the land‐ocean aquatic continuum as a result of biogeochemical processes in reservoirs. Both the changes in riverine nutrient fluxes and nutrient ratios impact ecosystem functioning of receiving water bodies. We utilize spatially distributed mechanistic models of nitrogen (N), phosphorus (P), and silicon (Si) cycling in reservoirs to quantify changes in nutrient stoichiometry of river discharge to coastal waters. The results demonstrate that the growing number of dams decouples the riverine fluxes of N, P, and Si. Worldwide, preferential removal of P over N in reservoirs increases N:P ratios delivered to the ocean, raising the potential for P limitation of coastal productivity. By midcentury, more than half of the rivers discharging to the coastal zone will experience a higher removal of reactive Si relative to reactive P and total N, in response to the rapid pace at which new hydroelectric dams are being built.

Insight into the microbiology of nitrogen cycle in the dairy manure composting process revealed by combining high-throughput sequencing and quantitative PCR

Nitrogen cycling during composting process is not yet fully understood. This study explored the key genes involved in nitrogen cycling during dairy manure composting process using high-throughput sequencing and quantitative PCR technologies. Results showed that nitrogen fixation occurred mainly during the thermophilic and cooling phases, and significantly enhanced the nitrogen content of compost. Thermoclostridium stercorarium was the main diazotroph. Ammonia oxidation occurred during the maturation phase and Nitrosomonas sp. was the most abundant ammonia oxidizing bacteria. Denitrification contributed to the greatest nitrogen loss during the composting process. The nirK community was dominated by Luteimonas sp. and Achromobacter sp., while the nirS community was dominated by Alcaligenes faecalis and Pseudomonas stutzeri. The nosZ community varied in a succession of Halomonas ilicicola, Pseudomonas flexibili and Labrenzia alba dominated communities according to different composting phases. Based on these results, nitrogen cycling models for different phases of the dairy manure composting process were established.

Elevated CO2 improved soil nitrogen mineralization capacity of rice paddy

Elevated CO2 would increase rice yields and may lead to nitrogen limitation and potentially influence the sustainability of agricultural production. Blindly increasing the amount of chemical fertilizer will damage the environment and is very unwise. Therefore, clarifying the response of soil nitrogen mineralization capacity to elevated CO2 is critical for both sustainable agriculture production and environmental protection. Here, we relied on Free-Air CO2 Enrichment (FACE) platform and used a waterlogged incubation method to investigate the effects of elevated CO2 on soil nitrogen mineralization capacity under different fertilization levels when planted different rice cultivars (strong and weak-CO2 response rice). According to the first-order kinetic equation fitting, compared with Ambient, elevated CO2 increased soil potential mineralized nitrogen (Np) by 16.18%. Path analysis indicated that fertilization status, rice cultivar, soil organic carbon and soil C: N ratio might affect Np. There was a significant positive correlation between soil nitrogen mineralization rate and Np. Under different fertilization conditions and rice cultivars, the improvement degree of soil nitrogen mineralization capacity (Np and soil nitrogen mineralization rate) by elevated CO2 was different. These findings suggest that more parameters and influencing factors should be taken into account when studying soil nitrogen cycle models under the condition of global change.

N<sub>2</sub>O changes from the Last Glacial Maximum to the preindustrial – Part 1: Quantitative reconstruction of terrestrial and marine emissions using N<sub>2</sub>O stable isotopes in ice cores

Abstract. Using high-precision and centennial-resolution ice core information on atmospheric nitrous oxide concentrations and its stable nitrogen and oxygen isotopic composition, we quantitatively reconstruct changes in the terrestrial and marine N2O emissions over the last 21 000 years. Our reconstruction indicates that N2O emissions from land and ocean increased over the deglaciation largely in parallel by 1.7±0.3 and 0.7±0.3 TgN yr−1, respectively, relative to the Last Glacial Maximum level. However, during the abrupt Northern Hemisphere warmings at the onset of the Bølling–Allerød warming and the end of the Younger Dryas, terrestrial emissions respond more rapidly to the northward shift in the Intertropical Convergence Zone connected to the resumption of the Atlantic Meridional Overturning Circulation. About 90 % of these large step increases were realized within 2 centuries at maximum. In contrast, marine emissions start to slowly increase already many centuries before the rapid warmings, possibly connected to a re-equilibration of subsurface oxygen in response to previous changes. Marine emissions decreased, concomitantly with changes in atmospheric CO2 and δ13C(CO2), at the onset of the termination and remained minimal during the early phase of Heinrich Stadial 1. During the early Holocene a slow decline in marine N2O emission of 0.4 TgN yr−1 is reconstructed, which suggests an improvement of subsurface water ventilation in line with slowly increasing Atlantic overturning circulation. In the second half of the Holocene total emissions remain on a relatively constant level, but with significant millennial variability. The latter is still difficult to attribute to marine or terrestrial sources. Our N2O emission records provide important quantitative benchmarks for ocean and terrestrial nitrogen cycle models to study the influence of climate on nitrogen turnover on timescales from several decades to glacial–interglacial changes.

Study on ecological dynamic model for phytoremediation of farmland drainage water

Eutrophication of surface water bodies resulted from pollution of the farmland drainage water is a serious water pollution problem recently, which has been paid more and more attention. Phytoremediation is to use specific aquatic or terrestrial plants for water pollutants such as nitrogen in situ absorption, degradation and fixing, with low cost and good effect of removing nitrogen pollutants, which is gradually used in sewage disposal. However, impact factor of phytoremediation is not known comprehensively enough since people have done research on it recently, and study on related dynamics is not much either. So many numerical models are not accurate enough, and the selection of reaction parameters is not unified. In this paper, two kinds of phytoremediation constructed wetlands such as surface wetland named Cell A and vertical flow constructed wetland named Cell B are used. Through experiments design, operation of the field and sampling analysis, combined with the principle of nitrogen transformation model, phytoremediation of farmland in the water cycle nitrogen ecology dynamic conceptual numerical model is developed. A kinetic model of nitrogen physical-chemical biodegradation reaction is established in this study. And genetic algorithm with the least-squares method of nonlinear optimization is used via MATLAB, which is an independent program to determine the unknown parameters in the developed model. Simulated results showed that the concentration of ammonia nitrogen, nitrate nitrogen and organic nitrogen in the effluent of Cell A is better than other pools. The ecological dynamic model of nitrogen cycle and the first-order kinetic model of removing total nitrogen (TN) were compared and analyzed in the simulation of the Cell A. It is indicated that the developed ecological dynamic model in this study could accurately predict the change of the effluent concentration in the constructed wetland. And the model in this study can provide theoretical support for the design and operation of the constructed wetland and the prediction of the effluent concentration to facilitate the mastery for the actual purification effect of the constructed wetland.

Assessing Marine Nitrogen Cycle Rates and Process Sensitivities With a Global 3‐D Inverse Model

AbstractWe present results from a global inverse marine nitrogen (N) cycle model that include nitrate (NO3−) and nitrite (NO2−) concentrations and their N isotopic compositions as constraints on N cycle process rates in marine oxygen deficient zones (ODZs). NO2− is an important intermediate in the N cycle, particularly in ODZs where it is a substrate in the N loss processes, denitrification, and anammox. Similar to earlier work, our model yields a total water column N loss rate of 61 ± 10 Tg N/year. However, by including NO2− and its N isotopic composition, we are able to assess the relative contributions of denitrification and anammox to N loss and examine some of the potential drivers of that balance. We find that anammox contributes 60% of global water column N loss, dominating N loss along the edges of ODZs, while denitrification is more important in the anoxic ODZ cores. The decoupling of anammox and denitrification is supported by NO2− oxidation, which co‐occurs with NO3− reduction and anammox in ODZs. High rates of NO2− oxidation (up to 400 nM/day), which are tightly coupled to heterotrophic NO3− reduction, are required to match NO3− and NO2− concentration and isotope observations in marine ODZs. Lowering the rate of NO2− oxidation in ODZs by adjusting O2‐sensitive parameters results in higher rates of water column N loss, highlighting the role of NO2− oxidation in maintaining the marine fixed N inventory.

Fraction of nitrous oxide production in nitrification and its effect on total soil emission: A meta-analysis and global-scale sensitivity analysis using a process-based model.

Nitrification in terrestrial soils is one of the major processes of emission of nitrous oxide (N2O), a potent greenhouse gas and stratospheric-ozone-depleting substance. We assessed the fraction of N2O emission associated with nitrification in soil through a meta-analysis and sensitivity analysis using a process-based model. We corrected observational values of gross nitrification and associated N2O emission rates from 71 records for various soils in the world spanning from 0.006% to 29.5%. We obtained a median value of 0.14%, and then assessed how the nitrification-associated N2O emission fraction has been considered in terrestrial nitrogen cycle models. Using a process-based biogeochemical model, we conducted a series of sensitivity analyses for the effects of different values of nitrification-associated N2O emission fraction on soil N2O emission. Using an empirical relationship between soil pH and nitrification-associated N2O emission fraction, the model well simulated global emission patterns (global total in the 2000s, 16.8 Tg N2O yr–1). Differences in the nitrification-associated N2O emission fraction caused differences in total N2O emission of as much as 2.5 Tg N2O yr–1. Therefore, to obtain reliable estimation of soil N2O emission for nitrogen and climate management, it is important to constrain the parameterization in models by ensuring extensive and accurate observations.

An incubation study of temperature sensitivity of greenhouse gas fluxes in three land-cover types near Sydney, Australia

Greenhouse gas (GHG) fluxes play crucial roles in regulating the Earth surface temperature. However, our understanding of the effect of land-cover and soil depth on the potential GHG fluxes and their temperature sensitivities (Q10) is limited, which consequently increases the uncertainty to predict GHG exchange between soils and the atmosphere. In the present study, we sampled soils with contrasting characteristics from three land-cover types (wetland, grassland, and forest) and soil depths (0–10, 10–20, and 20–30 cm) from the Cumberland Plain near Sydney, Australia, and incubated at optimal (60%) water holding capacity at three temperatures (15, 25, and 35 °C). Overall, GHG fluxes and Q10 values differed significantly among land-cover types and soil depths. CO2 and N2O emissions were highest in wetland followed by grassland and forest soils, and they decreased with soil depth. In contrast, CH4 uptake was highest in grassland followed by forest and wetland soils, and it increased with soil depth. Combining the three major GHGs, the global warming potential in soil from wetland was higher than that from grassland and forest. Moreover, Q10 values of CO2 and N2O emissions were: wetland > grassland > forest, while Q10 value of CH4 uptake showed the opposite pattern. Q10 values of CO2 and N2O emissions and CH4 uptake all increased with soil depth, demonstrating that subsoil has a higher potential for CO2 and N2O emissions and CH4 uptake in a warming climate. While these experiments were conducted under ideally controlled laboratory conditions, results suggest that the large carbon stocks in wetland soils are vulnerable to loss and thus may amplify climate warming; upland soils are crucial CH4 sinks and thus potentially mitigate climate change. In addition, the greater temperature sensitivities of CO2 and N2O emissions and CH4 uptake in subsoil should be accounted for in carbon and nitrogen cycling models.

A high-resolution biogeochemical model (ROMS 3.4 + bio_Fennel) of the East Australian Current system

Abstract. Understanding phytoplankton dynamics is critical across a range of topics, spanning from fishery management to climate change mitigation. It is particularly interesting in the East Australian Current (EAC) system, as the region's eddy field strongly conditions nutrient availability and therefore phytoplankton growth. Numerical models provide unparalleled insight into these biogeochemical dynamics. Yet, to date, modelling efforts off southeastern Australia have either targeted case studies (small spatial and temporal scales) or encompassed the whole EAC system but focused on climate change effects at the mesoscale (with a spatial resolution of 1/10∘). Here we couple a model of the pelagic nitrogen cycle (bio_Fennel) to a 10-year high-resolution (2.5–5 km horizontal) three-dimensional ocean model (ROMS) to resolve both regional and finer-scale biogeochemical processes occurring in the EAC system. We use several statistical metrics to compare the simulated surface chlorophyll to an ocean colour dataset (Copernicus-GlobColour) for the 2003–2011 period and show that the model can reproduce the observed phytoplankton surface patterns with a domain-wide RMSE of approximately 0.2 mg Chl a m−3 and a correlation coefficient of 0.76. This coupled configuration will provide a much-needed framework to examine phytoplankton variability in the EAC system providing insight into important ecosystem dynamics such as regional nutrient supply mechanisms and biogeochemical cycling occurring in EAC eddies.