Technical Note: A generic law-of-the-minimum flux limiter for simulating substrate limitation in biogeochemical models

DOE PAGES

Tang, J. Y.; Riley, W. J.

2016-02-05

We present a generic flux limiter to account for mass limitations from an arbitrary number of substrates in a biogeochemical reaction network. The flux limiter is based on the observation that substrate (e.g., nitrogen, phosphorus) limitation in biogeochemical models can be represented as to ensure mass conservative and non-negative numerical solutions to the governing ordinary differential equations. Application of the flux limiter includes two steps: (1) formulation of the biogeochemical processes with a matrix of stoichiometric coefficients and (2) application of Liebig's law of the minimum using the dynamic stoichiometric relationship of the reactants. This approach contrasts with the ad hoc down-regulationmore » approaches that are implemented in many existing models (such as CLM4.5 and the ACME (Accelerated Climate Modeling for Energy) Land Model (ALM)) of carbon and nutrient interactions, which are error prone when adding new processes, even for experienced modelers. Through an example implementation with a CENTURY-like decomposition model that includes carbon, nitrogen, and phosphorus, we show that our approach (1) produced almost identical results to that from the ad hoc down-regulation approaches under non-limiting nutrient conditions, (2) properly resolved the negative solutions under substrate-limited conditions where the simple clipping approach failed, (3) successfully avoided the potential conceptual ambiguities that are implied by those ad hoc down-regulation approaches. We expect our approach will make future biogeochemical models easier to improve and more robust.« less

Development of a 3D coupled physical-biogeochemical model for the Marseille coastal area (NW Mediterranean Sea): what complexity is required in the coastal zone?

PubMed

Fraysse, Marion; Pinazo, Christel; Faure, Vincent Martin; Fuchs, Rosalie; Lazzari, Paolo; Raimbault, Patrick; Pairaud, Ivane

2013-01-01

Terrestrial inputs (natural and anthropogenic) from rivers, the atmosphere and physical processes strongly impact the functioning of coastal pelagic ecosystems. The objective of this study was to develop a tool for the examination of these impacts on the Marseille coastal area, which experiences inputs from the Rhone River and high rates of atmospheric deposition. Therefore, a new 3D coupled physical/biogeochemical model was developed. Two versions of the biogeochemical model were tested, one model considering only the carbon (C) and nitrogen (N) cycles and a second model that also considers the phosphorus (P) cycle. Realistic simulations were performed for a period of 5 years (2007-2011). The model accuracy assessment showed that both versions of the model were able of capturing the seasonal changes and spatial characteristics of the ecosystem. The model also reproduced upwelling events and the intrusion of Rhone River water into the Bay of Marseille well. Those processes appeared to greatly impact this coastal oligotrophic area because they induced strong increases in chlorophyll-a concentrations in the surface layer. The model with the C, N and P cycles better reproduced the chlorophyll-a concentrations at the surface than did the model without the P cycle, especially for the Rhone River water. Nevertheless, the chlorophyll-a concentrations at depth were better represented by the model without the P cycle. Therefore, the complexity of the biogeochemical model introduced errors into the model results, but it also improved model results during specific events. Finally, this study suggested that in coastal oligotrophic areas, improvements in the description and quantification of the hydrodynamics and the terrestrial inputs should be preferred over increasing the complexity of the biogeochemical model.

Development of a 3D Coupled Physical-Biogeochemical Model for the Marseille Coastal Area (NW Mediterranean Sea): What Complexity Is Required in the Coastal Zone?

PubMed Central

Fraysse, Marion; Pinazo, Christel; Faure, Vincent Martin; Fuchs, Rosalie; Lazzari, Paolo; Raimbault, Patrick; Pairaud, Ivane

2013-01-01

Terrestrial inputs (natural and anthropogenic) from rivers, the atmosphere and physical processes strongly impact the functioning of coastal pelagic ecosystems. The objective of this study was to develop a tool for the examination of these impacts on the Marseille coastal area, which experiences inputs from the Rhone River and high rates of atmospheric deposition. Therefore, a new 3D coupled physical/biogeochemical model was developed. Two versions of the biogeochemical model were tested, one model considering only the carbon (C) and nitrogen (N) cycles and a second model that also considers the phosphorus (P) cycle. Realistic simulations were performed for a period of 5 years (2007–2011). The model accuracy assessment showed that both versions of the model were able of capturing the seasonal changes and spatial characteristics of the ecosystem. The model also reproduced upwelling events and the intrusion of Rhone River water into the Bay of Marseille well. Those processes appeared to greatly impact this coastal oligotrophic area because they induced strong increases in chlorophyll-a concentrations in the surface layer. The model with the C, N and P cycles better reproduced the chlorophyll-a concentrations at the surface than did the model without the P cycle, especially for the Rhone River water. Nevertheless, the chlorophyll-a concentrations at depth were better represented by the model without the P cycle. Therefore, the complexity of the biogeochemical model introduced errors into the model results, but it also improved model results during specific events. Finally, this study suggested that in coastal oligotrophic areas, improvements in the description and quantification of the hydrodynamics and the terrestrial inputs should be preferred over increasing the complexity of the biogeochemical model. PMID:24324589

A dynamic organic soil biogeochemical model for simulating the effects of wildfire on soil environmental conditions and carbon dynamics of black spruce forests

Treesearch

Shuhua Yi; A. David McGuire; Eric Kasischke; Jennifer Harden; Kristen Manies; Michelle Mack; Merritt Turetsky

2010-01-01

Ecosystem models have not comprehensively considered how interactions among fire disturbance, soil environmental conditions, and biogeochemical processes affect ecosystem dynamics in boreal forest ecosystems. In this study, we implemented a dynamic organic soil structure in the Terrestrial Ecosystem Model (DOS-TEM) to investigate the effects of fire on soil temperature...

Dimensionless Numbers For Morphological, Thermal And Biogeochemical Controls Of Hyporheic Processes

NASA Astrophysics Data System (ADS)

Bellin, Alberto; Marzadri, Alessandra; Tonina, Daniele

2013-04-01

Transport of solutes and heat within the hyporheic zone are interface processes that gained growing attention in the last decade, when several modelling strategies have been proposed, mainly at the local or reach scale. We propose to upscale local hyporheic biogeochemical processes to reach and network scales by means of a Lagrangian modelling framework, which allows to consider the impact of the flow structure on the processes modelled. This analysis shows that geochemical processes can be parametrized through two new Damköhler numbers, DaO, and DaT. DaO = ?up,50-?lim is defined as the ratio between the median hyporheic residence time, ?up,50 and the time of consuming dissolved oxygen to a prescribed threshold concentration, ?lim, below which reductive reactions are activated. It quantifies the biogeochemical status of the hyporheic zone and could be a metric for upscaling local hyporheic biogeochemical processes to reach and river-network scale processes. In addition, ?up,50 is the time scale of hyporheic advection; while ?lim is the representative time scale of biogeochemical reactions and indicates the distance along the streamline, measured as the time needed to travel that distance, that a particle of water travels before the dissolved oxygen concentration declines to [DO]lim, the value at which denitrification is activated. We show that DaO is representative of the redox status and indicates whether the hyporheic zone is a source or a sink of nitrate. Values of DaO larger than 1 indicate prevailing anaerobic conditions, whereas values smaller than 1 prevailing aerobic conditions. Similarly, DaT quantifies the importance of the temperature daily oscillations of the stream water on the hyporheic environment. It is defined as the ratio between ?up,50, and the time limit at which the ratio between the amplitude of the temperature oscillation within the hyporheic zone (evaluated along the streamline) and in the stream water is smaller than e-1. We show that values of DaT > 1 indicate a thermally stable hyporheic zone, where organism metabolism is not influenced by surface water thermal oscillations and biogeochemical reaction rates depend on the mean daily stream water temperature. Values smaller than 1 suggest that organisms need to adapt to the daily thermal variations and biogeochemical reaction rates will depend on the daily fluctuations induced by stream water.

A Unified Multi-scale Model for Cross-Scale Evaluation and Integration of Hydrological and Biogeochemical Processes

NASA Astrophysics Data System (ADS)

Liu, C.; Yang, X.; Bailey, V. L.; Bond-Lamberty, B. P.; Hinkle, C.

2013-12-01

Mathematical representations of hydrological and biogeochemical processes in soil, plant, aquatic, and atmospheric systems vary with scale. Process-rich models are typically used to describe hydrological and biogeochemical processes at the pore and small scales, while empirical, correlation approaches are often used at the watershed and regional scales. A major challenge for multi-scale modeling is that water flow, biogeochemical processes, and reactive transport are described using different physical laws and/or expressions at the different scales. For example, the flow is governed by the Navier-Stokes equations at the pore-scale in soils, by the Darcy law in soil columns and aquifer, and by the Navier-Stokes equations again in open water bodies (ponds, lake, river) and atmosphere surface layer. This research explores whether the physical laws at the different scales and in different physical domains can be unified to form a unified multi-scale model (UMSM) to systematically investigate the cross-scale, cross-domain behavior of fundamental processes at different scales. This presentation will discuss our research on the concept, mathematical equations, and numerical execution of the UMSM. Three-dimensional, multi-scale hydrological processes at the Disney Wilderness Preservation (DWP) site, Florida will be used as an example for demonstrating the application of the UMSM. In this research, the UMSM was used to simulate hydrological processes in rooting zones at the pore and small scales including water migration in soils under saturated and unsaturated conditions, root-induced hydrological redistribution, and role of rooting zone biogeochemical properties (e.g., root exudates and microbial mucilage) on water storage and wetting/draining. The small scale simulation results were used to estimate effective water retention properties in soil columns that were superimposed on the bulk soil water retention properties at the DWP site. The UMSM parameterized from smaller scale simulations were then used to simulate coupled flow and moisture migration in soils in saturated and unsaturated zones, surface and groundwater exchange, and surface water flow in streams and lakes at the DWP site under dynamic precipitation conditions. Laboratory measurements of soil hydrological and biogeochemical properties are used to parameterize the UMSM at the small scales, and field measurements are used to evaluate the UMSM.

Dynamic Biological Functioning Important for Simulating and Stabilizing Ocean Biogeochemistry

NASA Astrophysics Data System (ADS)

Buchanan, P. J.; Matear, R. J.; Chase, Z.; Phipps, S. J.; Bindoff, N. L.

2018-04-01

The biogeochemistry of the ocean exerts a strong influence on the climate by modulating atmospheric greenhouse gases. In turn, ocean biogeochemistry depends on numerous physical and biological processes that change over space and time. Accurately simulating these processes is fundamental for accurately simulating the ocean's role within the climate. However, our simulation of these processes is often simplistic, despite a growing understanding of underlying biological dynamics. Here we explore how new parameterizations of biological processes affect simulated biogeochemical properties in a global ocean model. We combine 6 different physical realizations with 6 different biogeochemical parameterizations (36 unique ocean states). The biogeochemical parameterizations, all previously published, aim to more accurately represent the response of ocean biology to changing physical conditions. We make three major findings. First, oxygen, carbon, alkalinity, and phosphate fields are more sensitive to changes in the ocean's physical state. Only nitrate is more sensitive to changes in biological processes, and we suggest that assessment protocols for ocean biogeochemical models formally include the marine nitrogen cycle to assess their performance. Second, we show that dynamic variations in the production, remineralization, and stoichiometry of organic matter in response to changing environmental conditions benefit the simulation of ocean biogeochemistry. Third, dynamic biological functioning reduces the sensitivity of biogeochemical properties to physical change. Carbon and nitrogen inventories were 50% and 20% less sensitive to physical changes, respectively, in simulations that incorporated dynamic biological functioning. These results highlight the importance of a dynamic biology for ocean properties and climate.

NIGHTHAWK - A Program for Modeling Saturated Batch and Column Experiments Incorporating Equilibrium and Kinetic Biogeochemistry

EPA Science Inventory

NIGHTHAWK simulates the fate and transport of biogeochemically reactive contaminants in the saturated subsurface. Version 1.2 supports batch and one- dimensional advective-dispersive-reactive transport involving a number of biogeochemical processes, including: microbially-mediate...

DOE Office of Scientific and Technical Information (OSTI.GOV)

Hammond, Glenn Edward; Yang, Xiaofan; Song, Xuehang

The groundwater-surface water interaction zone (GSIZ) plays an important role in riverine and watershed ecosystems as the exchange of waters of variable composition and temperature (hydrologic exchange flows) stimulate microbial activity and associated biogeochemical reactions. Variable temporal and spatial scales of hydrologic exchange flows, heterogeneity of the subsurface environment, and complexity of biogeochemical reaction networks in the GSIZ present challenges to incorporation of fundamental process representations and model parameterization across a range of spatial scales (e.g. from pore-scale to field scale). This paper presents a novel hybrid multiscale simulation approach that couples hydrologic-biogeochemical (HBGC) processes between two distinct length scalesmore » of interest.« less

The Water, Energy, and Biogeochemical Model (WEBMOD): A TOPMODEL application developed within the Modular Modeling System

NASA Astrophysics Data System (ADS)

Webb, R. M.; Wolock, D. M.; Linard, J. I.; Wieczorek, M. E.

2004-12-01

Process-based flow and transport simulation models can help increase understanding of how hydrologic flow paths affect biogeochemical mixing and reactions in watersheds. This presentation describes the Water, Energy, and Biogeochemical Model (WEBMOD), a new model designed to simulate water and chemical transport in both pristine and agricultural watersheds. WEBMOD simulates streamflow using TOPMODEL algorithms and also simulates irrigation, canopy interception, snowpack, and tile-drain flow; these are important processes for successful multi-year simulations of agricultural watersheds. In addition, the hydrologic components of the model are linked to the U.S. Geological Survey's (USGS) geochemical model PHREEQC such that solute chemistry for the hillslopes and streams also are computed. Model development, execution, and calibration take place within the USGS Modular Modeling System. WEBMOD is being validated at ten research watersheds. Five of these watersheds are nearly pristine and comprise the USGS Water, Energy, and Biogeochemical Budget (WEBB) Program field sites: Loch Vale, Colorado; Trout Lake, Wisconsin; Sleepers River, Vermont; Panola Mountain, Georgia; and the Luquillo Experimental Forest, Puerto Rico. The remaining five watersheds contain intensely cultivated fields being studied by USGS National Water Quality Assessment Program: Merced River, California; Granger Drain, Washington; Maple Creek, Nebraska; Sugar Creek, Indiana; and Morgan Creek, Delaware. Model calibration improved understanding of observed variations in soil moisture, solute concentrations, and stream discharge at the five WEBB watersheds and is now being set up to simulate the processes at the five agricultural watersheds that are now ending their first year of data collection.

Catchment hydro-biogeochemical response to forest harvest intensity and spatial pattern

EPA Science Inventory

We apply a new model, Visualizing Ecosystems for Land Management Assessment (VELMA), to Watershed 10 (WS10) in the H.J. Andrews Experimental Forest to simulate the effects of harvest intensity and spatial pattern on catchment hydrological and biogeochemical processes. Specificall...

Validation of an intermediate-complexity model for simulating marine biogeochemistry under anoxic conditions in the modern Black Sea

NASA Astrophysics Data System (ADS)

Romaniello, Stephen J.; Derry, Louis A.

2010-08-01

We test the ability of a new 1-D intermediate-complexity box model (ICBM) that includes process-based C, N, P, O, and S biogeochemistry to simulate profiles and fluxes of biogeochemically reactive species across a wide range of ocean redox states. The ICBM was developed to simulate whole ocean processes for paleoceanographic applications and has been tested with data from the modern global ocean. Here we adapt the circulation submodel of the ICBM to simulate water mass exchange and eddy diffusion processes in the Black Sea but make only very minor changes to the biogeochemical submodel. We force the model with estimated natural and anthropogenic inputs of tracers and nutrients to the Black Sea and compare the results of the simulations to modern observations. Ventilation of the Black Sea is modeled by depth-dependent entrainment of Cold Intermediate Layer water into Bosphorus plume water and subsequent intrusion into deep layers. The simulated profiles of circulation tracers θ, salinity, CFC-12, and radiocarbon agree well with available data, suggesting that the model does a reasonable job of representing physical exchange. Vertical profiles of biogeochemically active components are in good overall agreement with observations. The lack of trace metal (Mn and Fe) cycling in the model results in some discrepancies between the simulated profiles and observation across the suboxic zone; however, the overall redox balance is not sensitive to this difference. We compare modeled basin-wide biogeochemical fluxes to available estimates, but in a number of cases uncertainties in modern budgets limit our ability to test the model rigorously. In agreement with earlier work we find that fixed N losses via thiodenitrification are likely a major pathway in the Black Sea N cycle. Overall, the same biogeochemical submodel used to simulate the modern global ocean appears to perform well in simulating Black Sea processes without requiring significant modification. The ability of a single model to perform across a wide range of redox states is an important prerequisite for applying the ICBM to deep time paleoceanographic problems. The model source code is available as MATLAB™ 7 m-files provided as auxiliary material.

High resolution modelling of the biogeochemical processes in the eutrophic Loire River (France)

NASA Astrophysics Data System (ADS)

Minaudo, Camille; Moatar, Florentina; Curie, Florence; Gassama, Nathalie; Billen, Gilles

2016-04-01

A biogeochemical model was developed, coupling a physically based water temperature model (T-NET) with a semi-mechanistic biogeochemical model (RIVE, used in ProSe and Riverstrahler models) in order to assess at a fine temporal and spatial resolution the biogeochemical processes in the eutrophic Middle Loire hydrosystem (≈10 000 km², 3361 river segments). The code itself allows parallelized computing, which decreased greatly the calculation time (5 hours for simulating 3 years hourly). We conducted a daily survey during the period 2012-2014 at 2 sampling stations located in the Middle Loire of nutrients, chlorophyll pigments, phytoplankton and physic-chemical variables. This database was used as both input data (upstream Loire boundary) and validation data of the model (basin outlet). Diffuse and non-point sources were assessed based on a land cover analysis and WWTP datasets. The results appeared very sensible to the coefficients governing the dynamic of suspended solids and of phosphorus (sorption/desorption processes) within the model and some parameters needed to be estimated numerically. Both the Lagrangian point of view and fluxes budgets at the seasonal and event-based scale evidenced the biogeochemical functioning of the Loire River. Low discharge levels set up favorable physical conditions for phytoplankton growth (long water travel time, limited water depth, suspended particles sedimentation). Conversely, higher discharge levels highly limited the phytoplankton biomass (dilution of the colony, washing-out, limited travel time, remobilization of suspended sediments increasing turbidity), and most biogeochemical species were basically transferred downstream. When hydrological conditions remained favorable for phytoplankton development, P-availability was the critical factor. However, the model evidenced that most of the P in summer was recycled within the water body: on one hand it was assimilated by the algae biomass, and on the other hand it was released by mineralization of the dead cells. The high resolution of the model allowed understanding some fine temporal scale events, especially during some minor flood events occurring in summer. Paradoxically such events played two opposite roles: first it was disturbing the phytoplankton by diluting the biomass and remobilizing suspended sediments; then, it indirectly re-supplied the system with more available phosphorus, mainly because the washed-out phytoplankton could not assimilate the P available upstream. The model also pointed out the significant role played by Corbicula invasive clams in the river biogeochemical functioning, substantially reducing the phytoplankton biomass, and thus impacting the nutrients, oxygen and carbon cycles. However, the temporal and spatial distribution of Corbicula was questioned, and highlighted the need for data collection on this topic.

Scientific Discovery through Advanced Computing (SciDAC-3) Partnership Project Annual Report

DOE Office of Scientific and Technical Information (OSTI.GOV)

Hoffman, Forest M.; Bochev, Pavel B.; Cameron-Smith, Philip J..

The Applying Computationally Efficient Schemes for BioGeochemical Cycles ACES4BGC Project is advancing the predictive capabilities of Earth System Models (ESMs) by reducing two of the largest sources of uncertainty, aerosols and biospheric feedbacks, with a highly efficient computational approach. In particular, this project is implementing and optimizing new computationally efficient tracer advection algorithms for large numbers of tracer species; adding important biogeochemical interactions between the atmosphere, land, and ocean models; and applying uncertainty quanti cation (UQ) techniques to constrain process parameters and evaluate uncertainties in feedbacks between biogeochemical cycles and the climate system.

Simulation Based Exploration of Critical Zone Dynamics in Intensively Managed Landscapes

NASA Astrophysics Data System (ADS)

Kumar, P.

2017-12-01

The advent of high-resolution measurements of topographic and (vertical) vegetation features using areal LiDAR are enabling us to resolve micro-scale ( 1m) landscape structural characteristics over large areas. Availability of hyperspectral measurements is further augmenting these LiDAR data by enabling the biogeochemical characterization of vegetation and soils at unprecedented spatial resolutions ( 1-10m). Such data have opened up novel opportunities for modeling Critical Zone processes and exploring questions that were not possible before. We show how an integrated 3-D model at 1m grid resolution can enable us to resolve micro-topographic and ecological dynamics and their control on hydrologic and biogeochemical processes over large areas. We address the computational challenge of such detailed modeling by exploiting hybrid CPU and GPU computing technologies. We show results of moisture, biogeochemical, and vegetation dynamics from studies in the Critical Zone Observatory for Intensively managed Landscapes (IMLCZO) in the Midwestern United States.

Using NEON Data to Test and Refine Conceptual and Numerical Models of Soil Biogeochemical and Microbial Dynamics

NASA Astrophysics Data System (ADS)

Weintraub, S. R.; Stanish, L.; Ayers, E.

2017-12-01

Recent conceptual and numerical models have proposed new mechanisms that underpin key biogeochemical phenomena, including soil organic matter storage and ecosystem response to nitrogen deposition. These models seek to explicitly capture the ecological links among biota, especially microbes, and their physical and chemical environment to represent belowground pools and fluxes and how they respond to perturbation. While these models put forth exciting new concepts, their broad predictive abilities are unclear as some have been developed and tested against only small or regional datasets. The National Ecological Observatory Network (NEON) presents new opportunities to test and validate these models with multi-site data that span wide climatic, edaphic, and ecological gradients. NEON is measuring surface soil biogeochemical pools and fluxes along with diversity, abundance, and functional potential of soil microbiota at 47 sites distributed across the United States. This includes co-located measurements of soil carbon and nitrogen concentrations and stable isotopes, net nitrogen mineralization and nitrification rates, soil moisture, pH, microbial biomass, and community composition via 16S and ITS rRNA sequencing and shotgun metagenomic analyses. Early NEON data demonstrates that these wide edaphic and climatic gradients are related to changes in microbial community structure and functional potential, as well as element pools and process rates. Going forward, NEON's suite of standardized soil data has the potential to advance our understanding of soil communities and processes by allowing us to test the predictions of new soil biogeochemical frameworks and models. Here, we highlight several recently developed models that are ripe for this kind of data validation, and discuss key insights that may result. Further, we explore synergies with other networks, such as (i)LTER and (i)CZO, which may increase our ability to advance the frontiers of soil biogeochemical modeling.

Water, Carbon, and Nutrient Cycling Following Insect-induced Tree Mortality: How Well Do Plot-scale Observations Predict Ecosystem-Scale Response?

NASA Astrophysics Data System (ADS)

Brooks, P. D.; Barnard, H. R.; Biederman, J. A.; Borkhuu, B.; Edburg, S. L.; Ewers, B. E.; Gochis, D. J.; Gutmann, E. D.; Harpold, A. A.; Hicke, J. A.; Pendall, E.; Reed, D. E.; Somor, A. J.; Troch, P. A.

2011-12-01

Widespread tree mortality caused by insect infestations and drought has impacted millions of hectares across western North America in recent years. Although previous work on post-disturbance responses (e.g. experimental manipulations, fire, and logging) provides insight into how water and biogeochemical cycles may respond to insect infestations and drought, we find that the unique nature of these drivers of tree mortality complicates extrapolation to larger scales. Building from previous work on forest disturbance, we present a conceptual model of how temporal changes in forest structure impact the individual components of energy balance, hydrologic partitioning, and biogeochemical cycling and the interactions among them. We evaluate and refine this model using integrated observations and process modeling on multiple scales including plot, stand, flux tower footprint, hillslope, and catchment to identify scaling relationships and emergent patterns in hydrological and biogeochemical responses. Our initial results suggest that changes in forest structure at point or plot scales largely have predictable effects on energy, water, and biogeochemical cycles that are well captured by land surface, hydrological, and biogeochemical models. However, observations from flux towers and nested catchments suggest that both the hydrological and biogeochemical effects observed at tree and plot scales may be attenuated or exacerbated at larger scales. Compensatory processes are associated with attenuation (e.g. as transpiration decreases, evaporation and sublimation increase), whereas both attenuation and exacerbation may result from nonlinear scaling behavior across transitions in topography and ecosystem structure that affect the redistribution of energy, water, and solutes. Consequently, the effects of widespread tree mortality on ecosystem services of water supply and carbon sequestration will likely depend on how spatial patterns in mortality severity across the landscape affect large-scale hydrological partitioning.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Tang, J. Y.; Riley, W. J.

We present a generic flux limiter to account for mass limitations from an arbitrary number of substrates in a biogeochemical reaction network. The flux limiter is based on the observation that substrate (e.g., nitrogen, phosphorus) limitation in biogeochemical models can be represented as to ensure mass conservative and non-negative numerical solutions to the governing ordinary differential equations. Application of the flux limiter includes two steps: (1) formulation of the biogeochemical processes with a matrix of stoichiometric coefficients and (2) application of Liebig's law of the minimum using the dynamic stoichiometric relationship of the reactants. This approach contrasts with the ad hoc down-regulationmore » approaches that are implemented in many existing models (such as CLM4.5 and the ACME (Accelerated Climate Modeling for Energy) Land Model (ALM)) of carbon and nutrient interactions, which are error prone when adding new processes, even for experienced modelers. Through an example implementation with a CENTURY-like decomposition model that includes carbon, nitrogen, and phosphorus, we show that our approach (1) produced almost identical results to that from the ad hoc down-regulation approaches under non-limiting nutrient conditions, (2) properly resolved the negative solutions under substrate-limited conditions where the simple clipping approach failed, (3) successfully avoided the potential conceptual ambiguities that are implied by those ad hoc down-regulation approaches. We expect our approach will make future biogeochemical models easier to improve and more robust.« less

Natural and drought scenarios in an east central Amazon forest: Fidelity of the Community Land Model 3.5 with three biogeochemical models

NASA Astrophysics Data System (ADS)

Sakaguchi, Koichi; Zeng, Xubin; Christoffersen, Bradley J.; Restrepo-Coupe, Natalia; Saleska, Scott R.; Brando, Paulo M.

2011-03-01

Recent development of general circulation models involves biogeochemical cycles: flows of carbon and other chemical species that circulate through the Earth system. Such models are valuable tools for future projections of climate, but still bear large uncertainties in the model simulations. One of the regions with especially high uncertainty is the Amazon forest where large-scale dieback associated with the changing climate is predicted by several models. In order to better understand the capability and weakness of global-scale land-biogeochemical models in simulating a tropical ecosystem under the present day as well as significantly drier climates, we analyzed the off-line simulations for an east central Amazon forest by the Community Land Model version 3.5 of the National Center for Atmospheric Research and its three independent biogeochemical submodels (CASA', CN, and DGVM). Intense field measurements carried out under Large Scale Biosphere-Atmosphere Experiment in Amazonia, including forest response to drought from a throughfall exclusion experiment, are utilized to evaluate the whole spectrum of biogeophysical and biogeochemical aspects of the models. Our analysis shows reasonable correspondence in momentum and energy turbulent fluxes, but it highlights three processes that are not in agreement with observations: (1) inconsistent seasonality in carbon fluxes, (2) biased biomass size and allocation, and (3) overestimation of vegetation stress to short-term drought but underestimation of biomass loss from long-term drought. Without resolving these issues the modeled feedbacks from the biosphere in future climate projections would be questionable. We suggest possible directions for model improvements and also emphasize the necessity of more studies using a variety of in situ data for both driving and evaluating land-biogeochemical models.

Comparative Biogeochemical Cycles of Bioenergy Crops Reveal Nitrogen-Fixation and Low GHG Emissions in a Miscanthus x giganteus Agro-ecosystem

USDA-ARS?s Scientific Manuscript database

We evaluated the relative greenhouse gas mitigation potential of plant species considered as biofuel feedstock crops by simulating the biogeochemical processes associated with Miscanthus x giganteus, Panicum virgatum, Zea mays, and a mixed prairie community. DayCent model simulations for Miscanthus ...

Linking Chaotic Advection with Subsurface Biogeochemical Processes

NASA Astrophysics Data System (ADS)

Mays, D. C.; Freedman, V. L.; White, S. K.; Fang, Y.; Neupauer, R.

2017-12-01

This work investigates the extent to which groundwater flow kinematics drive subsurface biogeochemical processes. In terms of groundwater flow kinematics, we consider chaotic advection, whose essential ingredient is stretching and folding of plumes. Chaotic advection is appealing within the context of groundwater remediation because it has been shown to optimize plume spreading in the laminar flows characteristic of aquifers. In terms of subsurface biogeochemical processes, we consider an existing model for microbially-mediated reduction of relatively mobile uranium(VI) to relatively immobile uranium(IV) following injection of acetate into a floodplain aquifer beneath a former uranium mill in Rifle, Colorado. This model has been implemented in the reactive transport code eSTOMP, the massively parallel version of STOMP (Subsurface Transport Over Multiple Phases). This presentation will report preliminary numerical simulations in which the hydraulic boundary conditions in the eSTOMP model are manipulated to simulate chaotic advection resulting from engineered injection and extraction of water through a manifold of wells surrounding the plume of injected acetate. This approach provides an avenue to simulate the impact of chaotic advection within the existing framework of the eSTOMP code.

Demonstrating the Value of Fine-resolution Optical Data for Minimising Aliasing Impacts on Biogeochemical Models of Surface Waters

NASA Astrophysics Data System (ADS)

Chappell, N. A.; Jones, T.; Young, P.; Krishnaswamy, J.

2015-12-01

There is increasing awareness that under-sampling may have resulted in the omission of important physicochemical information present in water quality signatures of surface waters - thereby affecting interpretation of biogeochemical processes. For dissolved organic carbon (DOC) and nitrogen this under-sampling can now be avoided using UV-visible spectroscopy measured in-situ and continuously at a fine-resolution e.g. 15 minutes ("real time"). Few methods are available to extract biogeochemical process information directly from such high-frequency data. Jones, Chappell & Tych (2014 Environ Sci Technol: 13289-97) developed one such method using optically-derived DOC data based upon a sophisticated time-series modelling tool. Within this presentation we extend the methodology to quantify the minimum sampling interval required to avoid distortion of model structures and parameters that describe fundamental biogeochemical processes. This shifting of parameters which results from under-sampling is called "aliasing". We demonstrate that storm dynamics at a variety of sites dominate over diurnal and seasonal changes and that these must be characterised by sampling that may be sub-hourly to avoid aliasing. This is considerably shorter than that used by other water quality studies examining aliasing (e.g. Kirchner 2005 Phys Rev: 069902). The modelling approach presented is being developed into a generic tool to calculate the minimum sampling for water quality monitoring in systems driven primarily by hydrology. This is illustrated with fine-resolution, optical data from watersheds in temperate Europe through to the humid tropics.

Implementation ambiguity: The fifth element long lost in uncertainty budgets for land biogeochemical modeling

NASA Astrophysics Data System (ADS)

Tang, J.; Riley, W. J.

2015-12-01

Previous studies have identified four major sources of predictive uncertainty in modeling land biogeochemical (BGC) processes: (1) imperfect initial conditions (e.g., assumption of preindustrial equilibrium); (2) imperfect boundary conditions (e.g., climate forcing data); (3) parameterization (type I equifinality); and (4) model structure (type II equifinality). As if that were not enough to cause substantial sleep loss in modelers, we propose here a fifth element of uncertainty that results from implementation ambiguity that occurs when the model's mathematical description is translated into computational code. We demonstrate the implementation ambiguity using the example of nitrogen down regulation, a necessary process in modeling carbon-climate feedbacks. We show that, depending on common land BGC model interpretations of the governing equations for mineral nitrogen, there are three different implementations of nitrogen down regulation. We coded these three implementations in the ACME land model (ALM), and explored how they lead to different preindustrial and contemporary land biogeochemical states and fluxes. We also show how this implementation ambiguity can lead to different carbon-climate feedback estimates across the RCP scenarios. We conclude by suggesting how to avoid such implementation ambiguity in ESM BGC models.

Potential for real-time understanding of coupled hydrologic and biogeochemical processes in stream ecosystems: Future integration of telemetered data with process models for glacial meltwater streams

NASA Astrophysics Data System (ADS)

McKnight, Diane M.; Cozzetto, Karen; Cullis, James D. S.; Gooseff, Michael N.; Jaros, Christopher; Koch, Joshua C.; Lyons, W. Berry; Neupauer, Roseanna; Wlostowski, Adam

2015-08-01

While continuous monitoring of streamflow and temperature has been common for some time, there is great potential to expand continuous monitoring to include water quality parameters such as nutrients, turbidity, oxygen, and dissolved organic material. In many systems, distinguishing between watershed and stream ecosystem controls can be challenging. The usefulness of such monitoring can be enhanced by the application of quantitative models to interpret observed patterns in real time. Examples are discussed primarily from the glacial meltwater streams of the McMurdo Dry Valleys, Antarctica. Although the Dry Valley landscape is barren of plants, many streams harbor thriving cyanobacterial mats. Whereas a daily cycle of streamflow is controlled by the surface energy balance on the glaciers and the temporal pattern of solar exposure, the daily signal for biogeochemical processes controlling water quality is generated along the stream. These features result in an excellent outdoor laboratory for investigating fundamental ecosystem process and the development and validation of process-based models. As part of the McMurdo Dry Valleys Long-Term Ecological Research project, we have conducted field experiments and developed coupled biogeochemical transport models for the role of hyporheic exchange in controlling weathering reactions, microbial nitrogen cycling, and stream temperature regulation. We have adapted modeling approaches from sediment transport to understand mobilization of stream biomass with increasing flows. These models help to elucidate the role of in-stream processes in systems where watershed processes also contribute to observed patterns, and may serve as a test case for applying real-time stream ecosystem models.

A Hybrid Multiscale Framework for Subsurface Flow and Transport Simulations

DOE PAGES

Scheibe, Timothy D.; Yang, Xiaofan; Chen, Xingyuan; ...

2015-06-01

Extensive research efforts have been invested in reducing model errors to improve the predictive ability of biogeochemical earth and environmental system simulators, with applications ranging from contaminant transport and remediation to impacts of biogeochemical elemental cycling (e.g., carbon and nitrogen) on local ecosystems and regional to global climate. While the bulk of this research has focused on improving model parameterizations in the face of observational limitations, the more challenging type of model error/uncertainty to identify and quantify is model structural error which arises from incorrect mathematical representations of (or failure to consider) important physical, chemical, or biological processes, properties, ormore » system states in model formulations. While improved process understanding can be achieved through scientific study, such understanding is usually developed at small scales. Process-based numerical models are typically designed for a particular characteristic length and time scale. For application-relevant scales, it is generally necessary to introduce approximations and empirical parameterizations to describe complex systems or processes. This single-scale approach has been the best available to date because of limited understanding of process coupling combined with practical limitations on system characterization and computation. While computational power is increasing significantly and our understanding of biological and environmental processes at fundamental scales is accelerating, using this information to advance our knowledge of the larger system behavior requires the development of multiscale simulators. Accordingly there has been much recent interest in novel multiscale methods in which microscale and macroscale models are explicitly coupled in a single hybrid multiscale simulation. A limited number of hybrid multiscale simulations have been developed for biogeochemical earth systems, but they mostly utilize application-specific and sometimes ad-hoc approaches for model coupling. We are developing a generalized approach to hierarchical model coupling designed for high-performance computational systems, based on the Swift computing workflow framework. In this presentation we will describe the generalized approach and provide two use cases: 1) simulation of a mixing-controlled biogeochemical reaction coupling pore- and continuum-scale models, and 2) simulation of biogeochemical impacts of groundwater – river water interactions coupling fine- and coarse-grid model representations. This generalized framework can be customized for use with any pair of linked models (microscale and macroscale) with minimal intrusiveness to the at-scale simulators. It combines a set of python scripts with the Swift workflow environment to execute a complex multiscale simulation utilizing an approach similar to the well-known Heterogeneous Multiscale Method. User customization is facilitated through user-provided input and output file templates and processing function scripts, and execution within a high-performance computing environment is handled by Swift, such that minimal to no user modification of at-scale codes is required.« less

One multi-media environmental system with linkage between meteorology/ hydrology/ air quality models and water quality model

NASA Astrophysics Data System (ADS)

Tang, C.; Lynch, J. A.; Dennis, R. L.

2016-12-01

The biogeochemical processing of nitrogen and associated pollutants is driven by meteorological and hydrological processes in conjunction with pollutant loading. There are feedbacks between meteorology and hydrology that will be affected by land-use change and climate change. Changes in meteorology will affect pollutant deposition. It is important to account for those feedbacks and produce internally consistent simulations of meteorology, hydrology, and pollutant loading to drive the (watershed/water quality) biogeochemical models. In this study, the ecological response to emission reductions in streams in the Potomac watershed was evaluated. Firstly, we simulated the deposition by using the fully coupled Weather Research & Forecasting (WRF) model and the Community Multiscale Air Quality (CAMQ) model; secondly, we created the hydrological data by the offline linked Variable Infiltration Capacity (VIC) model and the WRF model. Lastly, we investigated the water quality by one comprehensive/environment model, namely the linkage of CMAQ, WRF, VIC and the Model of Acidification of Groundwater In Catchment (MAGIC) model from 2002 to 2010.The simulated results (such as NO3, SO4, and SBC) fit well to the observed values. The linkage provides a generally accurate, well-tested tool for evaluating sensitivities to varying meteorology and environmental changes on acidification and other biogeochemical processes, with capability to comprehensively explore strategic policy and management design.

Modelling the fate of six common pharmaceuticals in a small stream: quantification of attenuation and retention in different stream-specific environments

NASA Astrophysics Data System (ADS)

Riml, Joakim; Wörman, Anders; Kunkel, Uwe; Radke, Michael

2013-04-01

Detection of pharmaceutical residues in streaming waters is common in urbanized areas. Although the occurrence and source of these micropollutants is known, their behavior in these aquatic ecosystems is still only partly understood. Specifically, quantitative information of biogeochemical processes in stream-specific environments where predominant reactions occur is often missing. In an attempt to address this knowledge gap, we performed simultaneous tracer tests in Säva Brook, Sweden, with bezafibrate, clofibric acid, diclofenac, ibuprofen, metoprolol and naproxen, as well as with the more inert solutes uranine and Rhodamine WT. The breakthrough curves at five successive sampling stations along a 16 km long stream reach were evaluated using a coupled physical-biogeochemical model framework containing surface water transport together with a representation of transient storage in slow/immobile zones of the stream. The multi-tracer experiment opens for decoupling of hydrological and biogeochemical contribution to the fate, and by linking impact and sensitivity analyses to relative significance of model parameters the most important processes for each contaminant were elucidated. Specifically for Säva Brook, the proposed methodology revealed that the pharmaceutical-contaminated stream water remained in the storage zones for times corresponding to 5-25% of the flow time of the stream. Furthermore, the results indicate a great variability in terms of predominant biogeochemical processes between the different contaminants. Rapid reactions occurring in the transient storage zone attenuated both ibuprofen and clofibric acid, and we conclude that a major degradation pathway for these contaminants was biodegradation in the hyporheic zone. In contrast, bezafibrate, metoprolol, and naproxen were mainly affected by sorption both in the storage zone and the main channel, while diclofenac displayed negligible effects of biogeochemical reactions.

Biogeochemical Coupling between Ocean and Sea Ice

NASA Astrophysics Data System (ADS)

Wang, S.; Jeffery, N.; Maltrud, M. E.; Elliott, S.; Wolfe, J.

2016-12-01

Biogeochemical processes in ocean and sea ice are tightly coupled at high latitudes. Ongoing changes in Arctic and Antarctic sea ice domain likely influence the coupled system, not only through physical fields but also biogeochemical properties. Investigating the system and its changes requires representation of ocean and sea ice biogeochemical cycles, as well as their coupling in Earth System Models. Our work is based on ACME-HiLAT, a new offshoot of the Community Earth System Model (CESM), including a comprehensive representation of marine ecosystems in the form of the Biogeochemical Elemental Cycling Module (BEC). A full vertical column sea ice biogeochemical module has recently been incorporated into the sea ice component. We have further introduced code modifications to couple key growth-limiting nutrients (N, Si, Fe), dissolved and particulate organic matter, and phytoplankton classes that are important in polar regions between ocean and sea ice. The coupling of ocean and sea ice biology-chemistry will enable representation of key processes such as the release of important climate active constituents or seeding algae from melting sea ice into surface waters. Sensitivity tests suggest sea ice and ocean biogeochemical coupling influences phytoplankton competition, biological production, and the CO2 flux. Sea ice algal seeding plays an important role in determining phytoplankton composition of Arctic early spring blooms, since different groups show various responses to the seeding biomass. Iron coupling leads to increased phytoplankton biomass in the Southern Ocean, which also affects carbon uptake via the biological pump. The coupling of macronutrients and organic matter may have weaker influences on the marine ecosystem. Our developments will allow climate scientists to investigate the fully coupled responses of the sea ice-ocean BGC system to physical changes in polar climate.

Modeling seasonal variability of carbonate system parameters at the sediment -water interface in the Baltic Sea (Gdansk Deep)

NASA Astrophysics Data System (ADS)

Protsenko, Elizaveta; Yakubov, Shamil; Lessin, Gennady; Yakushev, Evgeniy; Sokołowski, Adam

2017-04-01

A one-dimensional fully-coupled benthic pelagic biogeochemical model BROM (Bottom RedOx Model) was used for simulations of seasonal variability of biogeochemical parameters in the upper sediment, Bottom Boundary Layer and the water column in the Gdansk Deep of the Baltic Sea. This model represents key biogeochemical processes of transformation of C, N, P, Si, O, S, Mn, Fe and the processes of vertical transport in the water column and the sediments. The hydrophysical block of BROM was forced by the output calculated with model GETM (General Estuarine Transport Model). In this study we focused on parameters of carbonate system at Baltic Sea, and mainly on their distributions near the sea-water interface. For validating of BROM we used field data (concentrations of main nutrients at water column and porewater of upper sediment) from the Gulf of Gdansk. The model allowed us to simulate the baseline ranges of seasonal variability of pH, Alkalinity, TIC and calcite/aragonite saturation as well as vertical fluxes of carbon in a region potentially selected for the CCS storage. This work was supported by project EEA CO2MARINE and STEMM-CCS.

Biogeochemistry of natural ponds in agricultural landscape: Lessons learned from modeling a kettle hole in Northeast Germany.

PubMed

Onandia, Gabriela; Lischeid, Gunnar; Kalettka, Thomas; Kleeberg, Andreas; Omari, Mohamed; Premke, Katrin; Arhonditsis, George B

2018-09-01

Kettle holes, small shallow ponds of glacial origin, represent hotspots for biodiversity and biogeochemical cycling. They abound in the young moraine landscape of Northeast Germany, potentially modulating element fluxes in a region where intensive agriculture prevails. The Rittgarten kettle hole, with semi-permanent hydroperiod and a surrounding reed belt, can be considered as a representative case study for such systems. Aiming to provide insights into the biogeochemical processes driving nutrient and primary producer dynamics in the Rittgarten kettle hole, we developed a mechanistic model that simulates the carbon, nitrogen, phosphorus and oxygen, phytoplankton, and free-floating macrophyte biomass dynamics. After model calibration and sensitivity analysis, our modeling exercise quantified the simulated nutrient fluxes associated with all the major biogeochemical processes considered by the model. Seasonality of nutrient concentrations, magnitude of primary productivity rates, and biogeochemical process characterization in the pond were reasonably reproduced by the model from July 2013 to July 2014. Our results suggest that the establishment of a phytoplankton community well-adapted to low light availability, together with the differential use of N and P from free-floating macrophytes and phytoplankton can explain their coexistence in kettle holes. Sediment nutrient release along with decomposition of decaying submerged macrophyte are essential drivers of internal nutrient cycling in kettle holes. Our results also suggest that the Rittgarten kettle hole act as a net source of CO 2 to the atmosphere on an annual scale, which offers a testable hypothesis for kettle holes with structural and functional similarities. We conclude by discussing the need to shed light on the effects of water level fluctuations on nutrient dynamics and biological succession patterns, as well as the relative importance of external sources and internal nutrient recycling mechanisms. Copyright © 2018 Elsevier B.V. All rights reserved.

Modeling and Assimilating Ocean Color Radiances

NASA Technical Reports Server (NTRS)

Gregg, Watson

2012-01-01

Radiances are the source of information from ocean color sensors to produce estimates of biological and geochemical constituents. They potentially provide information on various other aspects of global biological and chemical systems, and there is considerable work involved in deriving new information from these signals. Each derived product, however, contains errors that are derived from the application of the radiances, above and beyond the radiance errors. A global biogeochemical model with an explicit spectral radiative transfer model is used to investigate the potential of assimilating radiances. The results indicate gaps in our understanding of radiative processes in the oceans and their relationships with biogeochemical variables. Most important, detritus optical properties are not well characterized and produce important effects of the simulated radiances. Specifically, there does not appear to be a relationship between detrital biomass and its optical properties, as there is for chlorophyll. Approximations are necessary to get beyond this problem. In this reprt we will discuss the challenges in modeling and assimilation water-leaving radiances and the prospects for improving our understanding of biogeochemical process by utilizing these signals.

Microfluidic Experiments Studying Pore Scale Interactions of Microbes and Geochemistry

NASA Astrophysics Data System (ADS)

Chen, M.; Kocar, B. D.

2016-12-01

Understanding how physical phenomena, chemical reactions, and microbial behavior interact at the pore-scale is crucial to understanding larger scale trends in groundwater chemistry. Recent studies illustrate the utility of microfluidic devices for illuminating pore-scale physical-biogeochemical processes and their control(s) on the cycling of iron, uranium, and other important elements 1-3. These experimental systems are ideal for examining geochemical reactions mediated by microbes, which include processes governed by complex biological phenomenon (e.g. biofilm formation, etc.)4. We present results of microfluidic experiments using a model metal reducing bacteria and varying pore geometries, exploring the limitations of the microorganisms' ability to access tight pore spaces, and examining coupled biogeochemical-physical controls on the cycling of redox sensitive metals. Experimental results will provide an enhanced understanding of coupled physical-biogeochemical processes transpiring at the pore-scale, and will constrain and compliment continuum models used to predict and describe the subsurface cycling of redox-sensitive elements5. 1. Vrionis, H. A. et al. Microbiological and geochemical heterogeneity in an in situ uranium bioremediation field site. Appl. Environ. Microbiol. 71, 6308-6318 (2005). 2. Pearce, C. I. et al. Pore-scale characterization of biogeochemical controls on iron and uranium speciation under flow conditions. Environ. Sci. Technol. 46, 7992-8000 (2012). 3. Zhang, C., Liu, C. & Shi, Z. Micromodel investigation of transport effect on the kinetics of reductive dissolution of hematite. Environ. Sci. Technol. 47, 4131-4139 (2013). 4. Ginn, T. R. et al. Processes in microbial transport in the natural subsurface. Adv. Water Resour. 25, 1017-1042 (2002). 5. Scheibe, T. D. et al. Coupling a genome-scale metabolic model with a reactive transport model to describe in situ uranium bioremediation. Microb. Biotechnol. 2, 274-286 (2009).

Wetland eco-engineering: measuring and modeling feedbacks of oxidation processes between plants and clay-rich material

NASA Astrophysics Data System (ADS)

Saaltink, Rémon; Dekker, Stefan C.; Griffioen, Jasper; Wassen, Martin J.

2016-09-01

Interest is growing in using soft sediment as a foundation in eco-engineering projects. Wetland construction in the Dutch lake Markermeer is an example: here, dredging some of the clay-rich lake-bed sediment and using it to construct wetland will soon begin. Natural processes will be utilized during and after construction to accelerate ecosystem development. Knowing that plants can eco-engineer their environment via positive or negative biogeochemical plant-soil feedbacks, we conducted a 6-month greenhouse experiment to identify the key biogeochemical processes in the mud when Phragmites australis is used as an eco-engineering species. We applied inverse biogeochemical modeling to link observed changes in pore water composition to biogeochemical processes. Two months after transplantation we observed reduced plant growth and shriveling and yellowing of foliage. The N : P ratios of the plant tissue were low, and these were affected not by hampered uptake of N but by enhanced uptake of P. Subsequent analyses revealed high Fe concentrations in the leaves and roots. Sulfate concentrations rose drastically in our experiment due to pyrite oxidation; as reduction of sulfate will decouple Fe-P in reducing conditions, we argue that plant-induced iron toxicity hampered plant growth, forming a negative feedback loop, while simultaneously there was a positive feedback loop, as iron toxicity promotes P mobilization as a result of reduced conditions through root death, thereby stimulating plant growth and regeneration. Given these two feedback mechanisms, we propose the use of Fe-tolerant species rather than species that thrive in N-limited conditions. The results presented in this study demonstrate the importance of studying the biogeochemical properties of the situated sediment and the feedback mechanisms between plant and soil prior to finalizing the design of the eco-engineering project.

The roles of resuspension, diffusion and biogeochemical processes on oxygen dynamics offshore of the Rhône River, France: a numerical modeling study

NASA Astrophysics Data System (ADS)

Moriarty, Julia M.; Harris, Courtney K.; Fennel, Katja; Friedrichs, Marjorie A. M.; Xu, Kehui; Rabouille, Christophe

2017-04-01

Observations indicate that resuspension and associated fluxes of organic material and porewater between the seabed and overlying water can alter biogeochemical dynamics in some environments, but measuring the role of sediment processes on oxygen and nutrient dynamics is challenging. A modeling approach offers a means of quantifying these fluxes for a range of conditions, but models have typically relied on simplifying assumptions regarding seabed-water-column interactions. Thus, to evaluate the role of resuspension on biogeochemical dynamics, we developed a coupled hydrodynamic, sediment transport, and biogeochemical model (HydroBioSed) within the Regional Ocean Modeling System (ROMS). This coupled model accounts for processes including the storage of particulate organic matter (POM) and dissolved nutrients within the seabed; fluxes of this material between the seabed and the water column via erosion, deposition, and diffusion at the sediment-water interface; and biogeochemical reactions within the seabed. A one-dimensional version of HydroBioSed was then implemented for the Rhône subaqueous delta in France. To isolate the role of resuspension on biogeochemical dynamics, this model implementation was run for a 2-month period that included three resuspension events; also, the supply of organic matter, oxygen, and nutrients to the model was held constant in time. Consistent with time series observations from the Rhône Delta, model results showed that erosion increased the diffusive flux of oxygen into the seabed by increasing the vertical gradient of oxygen at the seabed-water interface. This enhanced supply of oxygen to the seabed, as well as resuspension-induced increases in ammonium availability in surficial sediments, allowed seabed oxygen consumption to increase via nitrification. This increase in nitrification compensated for the decrease in seabed oxygen consumption due to aerobic remineralization that occurred as organic matter was entrained into the water column. Additionally, entrainment of POM into the water column during resuspension events, and the associated increase in remineralization there, also increased oxygen consumption in the region of the water column below the pycnocline. During these resuspension events, modeled rates of oxygen consumption increased by factors of up to ˜ 2 and ˜ 8 in the seabed and below the pycnocline, respectively. When averaged over 2 months, the intermittent cycles of erosion and deposition led to a ˜ 16 % increase of oxygen consumption in the seabed, as well as a larger increase of ˜ 140 % below the pycnocline. These results imply that observations collected during quiescent periods, and biogeochemical models that neglect resuspension or use typical parameterizations for resuspension, may underestimate net oxygen consumption at sites like the Rhône Delta. Local resuspension likely has the most pronounced effect on oxygen dynamics at study sites with a high oxygen concentration in bottom waters, only a thin seabed oxic layer, and abundant labile organic matter.

Thresholds of understanding: Exploring assumptions of scale invariance vs. scale dependence in global biogeochemical models

NASA Astrophysics Data System (ADS)

Wieder, W. R.; Bradford, M.; Koven, C.; Talbot, J. M.; Wood, S.; Chadwick, O.

2016-12-01

High uncertainty and low confidence in terrestrial carbon (C) cycle projections reflect the incomplete understanding of how best to represent biologically-driven C cycle processes at global scales. Ecosystem theories, and consequently biogeochemical models, are based on the assumption that different belowground communities function similarly and interact with the abiotic environment in consistent ways. This assumption of "Scale Invariance" posits that environmental conditions will change the rate of ecosystem processes, but the biotic response will be consistent across sites. Indeed, cross-site comparisons and global-scale analyses suggest that climate strongly controls rates of litter mass loss and soil organic matter turnover. Alternatively, activities of belowground communities are shaped by particular local environmental conditions, such as climate and edaphic conditions. Under this assumption of "Scale Dependence", relationships generated by evolutionary trade-offs in acquiring resources and withstanding environmental stress dictate the activities of belowground communities and their functional response to environmental change. Similarly, local edaphic conditions (e.g. permafrost soils or reactive minerals that physicochemically stabilize soil organic matter on mineral surfaces) may strongly constrain the availability of substrates that biota decompose—altering the trajectory of soil biogeochemical response to perturbations. Identifying when scale invariant assumptions hold vs. where local variation in biotic communities or edaphic conditions must be considered is critical to advancing our understanding and representation of belowground processes in the face of environmental change. Here we introduce data sets that support assumptions of scale invariance and scale dependent processes and discuss their application in global-scale biogeochemical models. We identify particular domains over which assumptions of scale invariance may be appropriate and potential thresholds where shifts in ecosystem function may be expected. Finally, we discuss the mechanistic insight that can be applied in process-based models and datasets that can evaluate models across spatial and temporal scales.

The GRASP project - a multidisciplinary study of hydrology and biogeochemistry in a periglacial catchment area

NASA Astrophysics Data System (ADS)

Johansson, Emma; Lindborg, Tobias

2017-04-01

The Arctic region is sensitive to global warming, and permafrost thaw and release of old carbon are examples of processes that may have a positive feedback effect to the global climate system. Quantification and assumptions on future change are often based on model predictions. Such models require cross-disciplinary data of high quality that often is lacking. Biogeochemical processes in the landscape are highly influenced by the hydrology, which in turn is intimately related to permafrost processes. Thus, a multidisciplinary approach is needed when collecting data and setting up field experiments aiming at increase the understanding of these processes. Here we summarize and present data collected in the GRASP, Greenland Analogue Surface Project. GRASP is a catchment-scale field study of the periglacial area in the Kangerlussuaq region, West Greenland, focusing on hydrological and biogeochemical processes in the landscape. The site investigations were initiated in 2010 and have since then resulted in three separate data sets published in ESSD (Earth system and Science Data) each one focusing on i) meteorological data and hydrology, ii) biogeochemistry and iii) geometries of sediments and the active layer. The three data-sets, which are freely available via the PANGAEA data base, enable conceptual and coupled numerical modeling of hydrological and biogeochemical processes. An important strength with the GRASP data is that all data is collected within the same, relatively small, catchment area. This implies that measurements are more easily linked to the right source area or process. Despite the small catchment area it includes the major units of the periglacial hydrological system; a lake, a talik, a supra- and subpermafrost aquifer and, consequently, biogeochemical processes in each of these units may be studied. The new data from GRASP is both used with the aim to increase the knowledge of present day periglacial hydrology and biogeochemistry but also in order to predict consequences within these subjects of future climate change.

A review on vegetation models and applicability to climate simulations at regional scale

NASA Astrophysics Data System (ADS)

Myoung, Boksoon; Choi, Yong-Sang; Park, Seon Ki

2011-11-01

The lack of accurate representations of biospheric components and their biophysical and biogeochemical processes is a great source of uncertainty in current climate models. The interactions between terrestrial ecosystems and the climate include exchanges not only of energy, water and momentum, but also of carbon and nitrogen. Reliable simulations of these interactions are crucial for predicting the potential impacts of future climate change and anthropogenic intervention on terrestrial ecosystems. In this paper, two biogeographical (Neilson's rule-based model and BIOME), two biogeochemical (BIOME-BGC and PnET-BGC), and three dynamic global vegetation models (Hybrid, LPJ, and MC1) were reviewed and compared in terms of their biophysical and physiological processes. The advantages and limitations of the models were also addressed. Lastly, the applications of the dynamic global vegetation models to regional climate simulations have been discussed.

Preface [to special section on recent Loch Vale Watershed research

USGS Publications Warehouse

Baron, Jill S.; Williams, Mark W.

2000-01-01

Catchment-scale intensive and extensive research conducted over the last decade shows that our understanding of the biogeochemical and hydrologic processes in subalpine and alpine basins is not yet sufficiently mature to model and predict how biogeochemical transformations and surface water quality will change in response to climatic or human-driven changes in energy, water, and chemicals. A better understanding of these processes is needed for input to decision-making regulatory agencies and federal land managers. In recognition of this problem the National Research Council [1998] has identified as a critical research need an improved understanding of how global change will affect biogeochemical interactions with the hydrologic cycle and biogeochemical controls over the transport of water, nutrients, and materials from land to freshwater ecosystems. Improved knowledge of alpine and subalpine ecosystems is particularly important since high-elevation catchments are very sensitive to small changes in the flux of energy, chemicals, and water. Furthermore, alpine ecosystems may act as early warning indicators for ecosystem changes at lower elevations.

Insights on multivariate updates of physical and biogeochemical ocean variables using an Ensemble Kalman Filter and an idealized model of upwelling

NASA Astrophysics Data System (ADS)

Yu, Liuqian; Fennel, Katja; Bertino, Laurent; Gharamti, Mohamad El; Thompson, Keith R.

2018-06-01

Effective data assimilation methods for incorporating observations into marine biogeochemical models are required to improve hindcasts, nowcasts and forecasts of the ocean's biogeochemical state. Recent assimilation efforts have shown that updating model physics alone can degrade biogeochemical fields while only updating biogeochemical variables may not improve a model's predictive skill when the physical fields are inaccurate. Here we systematically investigate whether multivariate updates of physical and biogeochemical model states are superior to only updating either physical or biogeochemical variables. We conducted a series of twin experiments in an idealized ocean channel that experiences wind-driven upwelling. The forecast model was forced with biased wind stress and perturbed biogeochemical model parameters compared to the model run representing the "truth". Taking advantage of the multivariate nature of the deterministic Ensemble Kalman Filter (DEnKF), we assimilated different combinations of synthetic physical (sea surface height, sea surface temperature and temperature profiles) and biogeochemical (surface chlorophyll and nitrate profiles) observations. We show that when biogeochemical and physical properties are highly correlated (e.g., thermocline and nutricline), multivariate updates of both are essential for improving model skill and can be accomplished by assimilating either physical (e.g., temperature profiles) or biogeochemical (e.g., nutrient profiles) observations. In our idealized domain, the improvement is largely due to a better representation of nutrient upwelling, which results in a more accurate nutrient input into the euphotic zone. In contrast, assimilating surface chlorophyll improves the model state only slightly, because surface chlorophyll contains little information about the vertical density structure. We also show that a degradation of the correlation between observed subsurface temperature and nutrient fields, which has been an issue in several previous assimilation studies, can be reduced by multivariate updates of physical and biogeochemical fields.

Exploring the Influence of Topography on Belowground C Processes Using a Coupled Hydrologic-Biogeochemical Model

NASA Astrophysics Data System (ADS)

Shi, Y.; Davis, K. J.; Eissenstat, D. M.; Kaye, J. P.; Duffy, C.; Yu, X.; He, Y.

2014-12-01

Belowground carbon processes are affected by soil moisture and soil temperature, but current biogeochemical models are 1-D and cannot resolve topographically driven hill-slope soil moisture patterns, and cannot simulate the nonlinear effects of soil moisture on carbon processes. Coupling spatially-distributed physically-based hydrologic models with biogeochemical models may yield significant improvements in the representation of topographic influence on belowground C processes. We will couple the Flux-PIHM model to the Biome-BGC (BBGC) model. Flux-PIHM is a coupled physically-based land surface hydrologic model, which incorporates a land-surface scheme into the Penn State Integrated Hydrologic Model (PIHM). The land surface scheme is adapted from the Noah land surface model. Because PIHM is capable of simulating lateral water flow and deep groundwater, Flux-PIHM is able to represent the link between groundwater and the surface energy balance, as well as the land surface heterogeneities caused by topography. The coupled Flux-PIHM-BBGC model will be tested at the Susquehanna/Shale Hills critical zone observatory (SSHCZO). The abundant observations, including eddy covariance fluxes, soil moisture, groundwater level, sap flux, stream discharge, litterfall, leaf area index, above ground carbon stock, and soil carbon efflux, make SSHCZO an ideal test bed for the coupled model. In the coupled model, each Flux-PIHM model grid will couple a BBGC cell. Flux-PIHM will provide BBGC with soil moisture and soil temperature information, while BBGC provides Flux-PIHM with leaf area index. Preliminary results show that when Biome- BGC is driven by PIHM simulated soil moisture pattern, the simulated soil carbon is clearly impacted by topography.

Microbial community dynamics in soil aggregates shape biogeochemical gas fluxes from soil profiles - upscaling an aggregate biophysical model.

PubMed

Ebrahimi, Ali; Or, Dani

2016-09-01

Microbial communities inhabiting soil aggregates dynamically adjust their activity and composition in response to variations in hydration and other external conditions. These rapid dynamics shape signatures of biogeochemical activity and gas fluxes emitted from soil profiles. Recent mechanistic models of microbial processes in unsaturated aggregate-like pore networks revealed a highly dynamic interplay between oxic and anoxic microsites jointly shaped by hydration conditions and by aerobic and anaerobic microbial community abundance and self-organization. The spatial extent of anoxic niches (hotspots) flicker in time (hot moments) and support substantial anaerobic microbial activity even in aerated soil profiles. We employed an individual-based model for microbial community life in soil aggregate assemblies represented by 3D angular pore networks. Model aggregates of different sizes were subjected to variable water, carbon and oxygen contents that varied with soil depth as boundary conditions. The study integrates microbial activity within aggregates of different sizes and soil depth to obtain estimates of biogeochemical fluxes from the soil profile. The results quantify impacts of dynamic shifts in microbial community composition on CO2 and N2 O production rates in soil profiles in good agreement with experimental data. Aggregate size distribution and the shape of resource profiles in a soil determine how hydration dynamics shape denitrification and carbon utilization rates. Results from the mechanistic model for microbial activity in aggregates of different sizes were used to derive parameters for analytical representation of soil biogeochemical processes across large scales of practical interest for hydrological and climate models. © 2016 John Wiley & Sons Ltd.

Biogeochemical cycling in terrestrial ecosystems - Modeling, measurement, and remote sensing

NASA Technical Reports Server (NTRS)

Peterson, D. L.; Matson, P. A.; Lawless, J. G.; Aber, J. D.; Vitousek, P. M.

1985-01-01

The use of modeling, remote sensing, and measurements to characterize the pathways and to measure the rate of biogeochemical cycling in forest ecosystems is described. The application of the process-level model to predict processes in intact forests and ecosystems response to disturbance is examined. The selection of research areas from contrasting climate regimes and sites having a fertility gradient in that regime is discussed, and the sites studied are listed. The use of remote sensing in determining leaf area index and canopy biochemistry is analyzed. Nitrous oxide emission is investigated by using a gas measurement instrument. Future research projects, which include studying the influence of changes on nutrient cycling in ecosystems and the effect of pollutants on the ecosystems, are discussed.

A process-oriented hydro-biogeochemical model enabling simulation of gaseous carbon and nitrogen emissions and hydrologic nitrogen losses from a subtropical catchment.

PubMed

Zhang, Wei; Li, Yong; Zhu, Bo; Zheng, Xunhua; Liu, Chunyan; Tang, Jialiang; Su, Fang; Zhang, Chong; Ju, Xiaotang; Deng, Jia

2018-03-01

Quantification of nitrogen losses and net greenhouse gas (GHG) emissions from catchments is essential for evaluating the sustainability of ecosystems. However, the hydrologic processes without lateral flows hinder the application of biogeochemical models to this challenging task. To solve this issue, we developed a coupled hydrological and biogeochemical model, Catchment Nutrients Management Model - DeNitrification-DeComposition Model (CNMM-DNDC), to include both vertical and lateral mass flows. By incorporating the core biogeochemical processes (including decomposition, nitrification, denitrification and fermentation) of the DNDC into the spatially distributed hydrologic framework of the CNMM, the simulation of lateral water flows and their influences on nitrogen transportation can be realized. The CNMM-DNDC was then calibrated and validated in a small subtropical catchment belonged to Yanting station with comprehensive field observations. Except for the calibration of water flows (surface runoff and leaching water) in 2005, stream discharges of water and nitrate in 2007, the model validations of soil temperature, soil moisture, crop yield, water flows in 2006 and associated nitrate loss, fluxes of methane, ammonia, nitric oxide and nitrous oxide, and stream discharges of water and nitrate in 2008 were statistically in good agreement with the observations. Meanwhile, our initial simulation of the catchment showed scientific predictions. For instance, it revealed the following: (i) dominant ammonia volatilization among the losses of nitrogenous gases (accounting for 11-21% of the applied annual fertilizer nitrogen in croplands); (ii) hotspots of nitrate leaching near the main stream; and (iii) a net GHG sink function of the catchment. These results implicate the model's promising capability of predicting ecosystem productivity, hydrologic nitrogen loads, losses of gaseous nitrogen and emissions of GHGs, which could be used to provide strategies for establishing sustainable catchments. In addition, the model's capability would be further proved by applying in other catchments with different backgrounds. Copyright © 2017. Published by Elsevier B.V.

Diel biogeochemical processes and their effect on the aqueous chemistry of streams: A review

USGS Publications Warehouse

Nimick, David A.; Gammons, Christopher H.; Parker, Stephen R.

2011-01-01

This review summarizes biogeochemical processes that operate on diel, or 24-h, time scales in streams and the changes in aqueous chemistry that are associated with these processes. Some biogeochemical processes, such as those producing diel cycles of dissolved O2 and pH, were the first to be studied, whereas processes producing diel concentration cycles of a broader spectrum of chemical species including dissolved gases, dissolved inorganic and organic carbon, trace elements, nutrients, stable isotopes, and suspended particles have received attention only more recently. Diel biogeochemical cycles are interrelated because the cyclical variations produced by one biogeochemical process commonly affect another. Thus, understanding biogeochemical cycling is essential not only for guiding collection and interpretation of water-quality data but also for geochemical and ecological studies of streams. Expanded knowledge of diel biogeochemical cycling will improve understanding of how natural aquatic environments function and thus lead to better predictions of how stream ecosystems might react to changing conditions of contaminant loading, eutrophication, climate change, drought, industrialization, development, and other factors.

Benthic-Pelagic Coupling in Biogeochemical and Climate Models: Existing Approaches, Recent developments and Roadblocks

NASA Astrophysics Data System (ADS)

Arndt, Sandra

2016-04-01

Marine sediments are key components in the Earth System. They host the largest carbon reservoir on Earth, provide the only long term sink for atmospheric CO2, recycle nutrients and represent the most important climate archive. Biogeochemical processes in marine sediments are thus essential for our understanding of the global biogeochemical cycles and climate. They are first and foremost, donor controlled and, thus, driven by the rain of particulate material from the euphotic zone and influenced by the overlying bottom water. Geochemical species may undergo several recycling loops (e.g. authigenic mineral precipitation/dissolution) before they are either buried or diffuse back to the water column. The tightly coupled and complex pelagic and benthic process interplay thus delays recycling flux, significantly modifies the depositional signal and controls the long-term removal of carbon from the ocean-atmosphere system. Despite the importance of this mutual interaction, coupled regional/global biogeochemical models and (paleo)climate models, which are designed to assess and quantify the transformations and fluxes of carbon and nutrients and evaluate their response to past and future perturbations of the climate system either completely neglect marine sediments or incorporate a highly simplified representation of benthic processes. On the other end of the spectrum, coupled, multi-component state-of-the-art early diagenetic models have been successfully developed and applied over the past decades to reproduce observations and quantify sediment-water exchange fluxes, but cannot easily be coupled to pelagic models. The primary constraint here is the high computation cost of simulating all of the essential redox and equilibrium reactions within marine sediments that control carbon burial and benthic recycling fluxes: a barrier that is easily exacerbated if a variety of benthic environments are to be spatially resolved. This presentation provides an integrative overview of the benthic-pelagic coupling that accounts for the complex process interplay from the euphotic ocean to the deep sediment. It explores the intensity of the benthic-pelagic coupling across different environments and from the seasonal to the geological timescale. Different modelling approaches of coupling sediment and water column dynamics in regional/global biogeochemical models and (paleo)climate models are critically evaluated and their most important limitations, as well as the implications for our ability to predict the response of the global carbon cycle to past or future perturbations is discussed. Finally, the presentation identifies major roadblocks to the development of new model approaches and highlights how new techniques, new observational and laboratory data, as well as a close interdisciplinary collaboration can overcome these roadblocks.

A review of ion and metal pollutants in urban green water infrastructures.

PubMed

Kabir, Md Imran; Daly, Edoardo; Maggi, Federico

2014-02-01

In urban environments, the breakdown of chemicals and pollutants, especially ions and metal compounds, can be favoured by green water infrastructures (GWIs). The overall aim of this review is to set the basis to model GWIs using deterministic approaches in contrast to empirical ones. If a better picture of chemicals and pollutant input and an improved understanding of hydrological and biogeochemical processes affecting these pollutants were known, GWIs could be designed to efficiently retain these pollutants for site-specific meteorological patterns and pollutant load. To this end, we surveyed the existing literature to retrieve a comprehensive dataset of anions and cations, and alkaline and transition metal pollutants incoming to urban environments. Based on this survey, we assessed the pollution load and ecological risk indexes for metals. The existing literature was then surveyed to review the metal retention efficiency of GWIs, and possible biogeochemical processes related to inorganic metal compounds were proposed that could be integrated in biogeochemical models of GWIs. © 2013.

The interplay between estuarine transport and biogeochemical processes in determining the nutrient conditions in bottom layers of non-tidal Gulf of Finland

NASA Astrophysics Data System (ADS)

Kõuts, Mariliis; Raudsepp, Urmas; Maljutenko, Ilja

2017-04-01

In coastal areas, especially estuaries, spatial distribution and seasonal cycling of chemical and biological variables is largely determined by local biogeochemical processes and water transport of different properties. In tidal estuaries, however, biogeochemical processes are affected by tides as frequent water exchange alters nutrient and oxygen concentrations. In wide and deep non-tidal estuary-type marginal seas spatial distribution and seasonal cycling are determined by the mixture of water transport and local biogeochemistry. The Baltic Sea is a stratified water basin where halocline divides the water column into two parts: upper layer, which is horizontally uniform in terms of distribution of chemical and biological parameters, and has clear seasonal cycle; and bottom part, where nutrient and oxygen dynamics is more complex. There water transport and sediment-water interface fluxes play a major role. Our prime focus is the Gulf of Finland in the Baltic Sea. It is a wide, non-tidal and stratified sub-basin known for its high nutrient concentrations and severe oxygen deficiency in summer. We modelled the Baltic Sea (including Gulf of Finland) using ERGOM, a biogeochemical model coupled with circulation model GETM. Seasonal cycling and water circulation were observed with a 40-year simulation from 1966 to 2006. Our results show that in shallow areas above halocline the seasonal cycle of phytoplankton, nutrients and oxygen concentrations is uniform in space. Water circulation does not create inhomogeneous distribution pattern of biogeochemical parameters and their seasonal cycle. The circulation in the Gulf of Finland is strongly modulated by the seasonality of estuarine transport. Below the halocline saline low-oxygen and nutrient-rich water is transported from the open Baltic Proper to the Gulf of Finland in spring and early summer. This results in the highest nutrient concentrations and the poorest oxygen conditions by the end of August. In the shallow area nutrients have high concentrations in March-April before the spring bloom of diatoms starts. Low oxygen and nutrient concentrations are observed at the end of August. There is a qualitative difference of nutrient dynamics between shallow and deep layers but quantification of the role of transport and local biogeochemical processes is still challenging.

Emergent Archetype Hydrological-Biogeochemical Response Patterns in Heterogeneous Catchments

NASA Astrophysics Data System (ADS)

Jawitz, J. W.; Gall, H. E.; Rao, P.

2013-12-01

What can spatiotemporally integrated patterns observed in stream hydrologic and biogeochemical signals generated in response to transient hydro-climatic and anthropogenic forcing tell us about the interactions between spatially heterogeneous soil-mediated hydrological and biogeochemical processes? We seek to understand how the spatial structure of solute sources coupled with hydrologic responses affect observed concentration-discharge (C-Q) patterns. These patterns are expressions of the spatiotemporal structure of solute loads exported from managed catchments, and their likely ecological consequences manifested in receiving water bodies (e.g., wetlands, rivers, lakes, and coastal waters). We investigated the following broad questions: (1) How does the correlation between flow-generating areas and biogeochemical source areas across a catchment evolve under stochastic hydro-climatic forcing? (2) What are the feasible hydrologic and biogeochemical responses that lead to the emergence of the observed archetype C-Q patterns? and; (3) What implications do these coupled dynamics have for catchment monitoring and implementation of management practices? We categorize the observed temporal signals into three archetypical C-Q patterns: dilution; accretion, and constant concentration. We introduce a parsimonious stochastic model of heterogeneous catchments, which act as hydrologic and biogeochemical filters, to examine the relationship between spatial heterogeneity and temporal history of solute export signals. The core concept of the modeling framework is considering the types and degree of spatial correlation between solute source zones and flow generating zones, and activation of different portions of the catchments during rainfall events. Our overarching hypothesis is that each of the archetype C-Q patterns can be generated by explicitly linking landscape-scale hydrologic responses and spatial distributions of solute source properties within a catchment. The model simulations reproduce the three major C-Q patterns observed in published data, offering valuable insight into coupled catchment processes. The findings have important implications for effective catchment management for water quality improvement, and stream monitoring strategies.

Global changes in biogeochemical cycles in response to human activities

NASA Technical Reports Server (NTRS)

Moore, Berrien, III; Melillo, Jerry

1994-01-01

The main objective of our research was to characterize biogeochemical cycles at continental and global scales in both terrestrial and aquatic ecosystems. This characterization applied to both natural ecosystems and those disturbed by human activity. The primary elements of interest were carbon and nitrogen and the analysis sought to quantify standing stocks and dynamic cycling processes. The translocation of major nutrients from the terrestrial landscape to the atmosphere (via trace gases) and to fluvial systems (via leaching, erosional losses, and point source pollution) were of particular importance to this study. Our aim was to develop the first generation of Earth System Models. Our research was organized around the construction and testing of component biogeochemical models which treated terrestrial ecosystem processes, aquatic nutrient transport through drainage basins, and trace gas exchanges at the continental and global scale. A suite of three complementary models were defined within this construct. The models were organized to operate at a 1/2 degree latitude by longitude level of spatial resolution and to execute at a monthly time step. This discretization afforded us the opportunity to understand the dynamics of the biosphere down to subregional scales, while simultaneously placing these dynamics into a global context.

Use of the Water, Energy, and Biogeochemical Model (WEBMOD) to Simulate Water Quality at Five U.S. Geological Survey Research Watersheds

NASA Astrophysics Data System (ADS)

Webb, R. M.; Leavesley, G. H.; Shanley, J. B.; Peters, N. E.; Aulenbach, B. T.; Blum, A. E.; Campbell, D. H.; Clow, D. W.; Mast, M. A.; Stallard, R. F.; Larsen, M. C.; Troester, J. W.; Walker, J. F.; White, A. F.

2003-12-01

The Water, Energy, and Biogeochemical Model (WEBMOD) was developed as an aid to compare and contrast basic hydrologic and biogeochemical processes active in the diverse hydroclimatic regions represented by the five U.S. Geological Survey (USGS) Water, Energy, and Biogeochemical Budget (WEBB) sites: Loch Vale, Colorado; Trout Lake, Wisconsin; Sleepers River, Vermont; Panola Mountain, Georgia; and Luquillo Experimental Forest, Puerto Rico. WEBMOD simulates solute concentrations for vegetation canopy, snow pack, impermeable ground, leaf litter, unsaturated and saturated soil zones, riparian zones and streams using selected process modules coupled within the USGS Modular Modeling System (MMS). Source codes for the MMS hydrologic modules include the USGS Precipitation Runoff Modeling System, the National Weather Service Hydro-17 snow model, and TOPMODEL. The hydrologic modules distribute precipitation and temperature to predict evapotranspiration, snow accumulation, snow melt, and streamflow. Streamflow generation mechanisms include infiltration excess, saturated overland flow, preferential lateral flow, and base flow. Input precipitation chemistry, and fluxes and residence times predicted by the hydrologic modules are input into the geochemical module where solute concentrations are computed for a series of discrete well-mixed reservoirs using calls to the geochemical engine PHREEQC. WEBMOD was used to better understand variations in water quality observed at the WEBB sites from October 1991 through September 1997. Initial calibrations were completed by fitting the simulated hydrographs with those measured at the watershed outlets. Model performance was then refined by comparing the predicted export of conservative chemical tracers such as chloride, with those measured at the watershed outlets. The model succeeded in duplicating the temporal variability of net exports of major ions from the watersheds.

SHIMMER (1.0): a novel mathematical model for microbial and biogeochemical dynamics in glacier forefield ecosystems

NASA Astrophysics Data System (ADS)

Bradley, J. A.; Anesio, A. M.; Singarayer, J. S.; Heath, M. R.; Arndt, S.

2015-10-01

SHIMMER (Soil biogeocHemIcal Model for Microbial Ecosystem Response) is a new numerical modelling framework designed to simulate microbial dynamics and biogeochemical cycling during initial ecosystem development in glacier forefield soils. However, it is also transferable to other extreme ecosystem types (such as desert soils or the surface of glaciers). The rationale for model development arises from decades of empirical observations in glacier forefields, and enables a quantitative and process focussed approach. Here, we provide a detailed description of SHIMMER, test its performance in two case study forefields: the Damma Glacier (Switzerland) and the Athabasca Glacier (Canada) and analyse sensitivity to identify the most sensitive and unconstrained model parameters. Results show that the accumulation of microbial biomass is highly dependent on variation in microbial growth and death rate constants, Q10 values, the active fraction of microbial biomass and the reactivity of organic matter. The model correctly predicts the rapid accumulation of microbial biomass observed during the initial stages of succession in the forefields of both the case study systems. Primary production is responsible for the initial build-up of labile substrate that subsequently supports heterotrophic growth. However, allochthonous contributions of organic matter, and nitrogen fixation, are important in sustaining this productivity. The development and application of SHIMMER also highlights aspects of these systems that require further empirical research: quantifying nutrient budgets and biogeochemical rates, exploring seasonality and microbial growth and cell death. This will lead to increased understanding of how glacier forefields contribute to global biogeochemical cycling and climate under future ice retreat.

Concentration-discharge relationships to understand the interplay between hydrological and biogeochemical processes: insights from data analysis and numerical experiments in headwater catchments.

NASA Astrophysics Data System (ADS)

De Dreuzy, J. R.; Marçais, J.; Moatar, F.; Minaudo, C.; Courtois, Q.; Thomas, Z.; Longuevergne, L.; Pinay, G.

2017-12-01

Integration of hydrological and biogeochemical processes led to emerging patterns at the catchment scale. Monitoring in rivers reflects the aggregation of these effects. While discharge time series have been measured for decades, high frequency water quality monitoring in rivers now provides prominent measurements to characterize the interplay between hydrological and biogeochemical processes, especially to infer the processes that happen in the heterogeneous subsurface. However, we still lack frameworks to relate observed patterns to specific processes, because of the "organized complexity" of hydrological systems. Indeed, it is unclear what controls, for example, patterns in concentration-discharge (C/Q) relationships due to non-linear processes and hysteresis effects. Here we develop a non-intensive process-based model to test how the integration of different landforms (i.e. geological heterogeneities and structures, topographical features) with different biogeochemical reactivity assumptions (e.g. reactive zone locations) can shape the overall water quality time series. With numerical experiments, we investigate typical patterns in high frequency C/Q relationships. In headwater basins, we found that typical hysteretic patterns in C/Q relationships observed in data time series can be attributed to differences in water and solute locations stored across the hillslope. At the catchment scale though, these effects tend to average out by integrating contrasted hillslopes' landforms. Together these results suggest that information contained in headwater water quality monitoring can be used to understand how hydrochemical processes determine downstream conditions.

Adapting to life: ocean biogeochemical modelling and adaptive remeshing

NASA Astrophysics Data System (ADS)

Hill, J.; Popova, E. E.; Ham, D. A.; Piggott, M. D.; Srokosz, M.

2013-11-01

An outstanding problem in biogeochemical modelling of the ocean is that many of the key processes occur intermittently at small scales, such as the sub-mesoscale, that are not well represented in global ocean models. As an example, state-of-the-art models give values of primary production approximately two orders of magnitude lower than those observed in the ocean's oligotrophic gyres, which cover a third of the Earth's surface. This is partly due to their failure to resolve sub-mesoscale phenomena, which play a significant role in nutrient supply. Simply increasing the resolution of the models may be an inefficient computational solution to this problem. An approach based on recent advances in adaptive mesh computational techniques may offer an alternative. Here the first steps in such an approach are described, using the example of a~simple vertical column (quasi 1-D) ocean biogeochemical model. We present a novel method of simulating ocean biogeochemical behaviour on a vertically adaptive computational mesh, where the mesh changes in response to the biogeochemical and physical state of the system throughout the simulation. We show that the model reproduces the general physical and biological behaviour at three ocean stations (India, Papa and Bermuda) as compared to a high-resolution fixed mesh simulation and to observations. The simulations capture both the seasonal and inter-annual variations. The use of an adaptive mesh does not increase the computational error, but reduces the number of mesh elements by a factor of 2-3, so reducing computational overhead. We then show the potential of this method in two case studies where we change the metric used to determine the varying mesh sizes in order to capture the dynamics of chlorophyll at Bermuda and sinking detritus at Papa. We therefore demonstrate adaptive meshes may provide a~suitable numerical technique for simulating seasonal or transient biogeochemical behaviour at high spatial resolution whilst minimising computational cost.

Adapting to life: ocean biogeochemical modelling and adaptive remeshing

NASA Astrophysics Data System (ADS)

Hill, J.; Popova, E. E.; Ham, D. A.; Piggott, M. D.; Srokosz, M.

2014-05-01

An outstanding problem in biogeochemical modelling of the ocean is that many of the key processes occur intermittently at small scales, such as the sub-mesoscale, that are not well represented in global ocean models. This is partly due to their failure to resolve sub-mesoscale phenomena, which play a significant role in vertical nutrient supply. Simply increasing the resolution of the models may be an inefficient computational solution to this problem. An approach based on recent advances in adaptive mesh computational techniques may offer an alternative. Here the first steps in such an approach are described, using the example of a simple vertical column (quasi-1-D) ocean biogeochemical model. We present a novel method of simulating ocean biogeochemical behaviour on a vertically adaptive computational mesh, where the mesh changes in response to the biogeochemical and physical state of the system throughout the simulation. We show that the model reproduces the general physical and biological behaviour at three ocean stations (India, Papa and Bermuda) as compared to a high-resolution fixed mesh simulation and to observations. The use of an adaptive mesh does not increase the computational error, but reduces the number of mesh elements by a factor of 2-3. Unlike previous work the adaptivity metric used is flexible and we show that capturing the physical behaviour of the model is paramount to achieving a reasonable solution. Adding biological quantities to the adaptivity metric further refines the solution. We then show the potential of this method in two case studies where we change the adaptivity metric used to determine the varying mesh sizes in order to capture the dynamics of chlorophyll at Bermuda and sinking detritus at Papa. We therefore demonstrate that adaptive meshes may provide a suitable numerical technique for simulating seasonal or transient biogeochemical behaviour at high vertical resolution whilst minimising the number of elements in the mesh. More work is required to move this to fully 3-D simulations.

Belowground Carbon Cycling Processes at the Molecular Scale: An EMSL Science Theme Advisory Panel Workshop

DOE Office of Scientific and Technical Information (OSTI.GOV)

Hess, Nancy J.; Brown, Gordon E.; Plata, Charity

2014-02-21

As part of the Belowground Carbon Cycling Processes at the Molecular Scale workshop, an EMSL Science Theme Advisory Panel meeting held in February 2013, attendees discussed critical biogeochemical processes that regulate carbon cycling in soil. The meeting attendees determined that as a national scientific user facility, EMSL can provide the tools and expertise needed to elucidate the molecular foundation that underlies mechanistic descriptions of biogeochemical processes that control carbon allocation and fluxes at the terrestrial/atmospheric interface in landscape and regional climate models. Consequently, the workshop's goal was to identify the science gaps that hinder either development of mechanistic description ofmore » critical processes or their accurate representation in climate models. In part, this report offers recommendations for future EMSL activities in this research area. The workshop was co-chaired by Dr. Nancy Hess (EMSL) and Dr. Gordon Brown (Stanford University).« less

Characterizing multiple timescales of stream and storage zone interaction that affect solute fate and transport in streams

USGS Publications Warehouse

Choi, Jungyill; Harvey, Judson W.; Conklin, Martha H.

2000-01-01

The fate of contaminants in streams and rivers is affected by exchange and biogeochemical transformation in slowly moving or stagnant flow zones that interact with rapid flow in the main channel. In a typical stream, there are multiple types of slowly moving flow zones in which exchange and transformation occur, such as stagnant or recirculating surface water as well as subsurface hyporheic zones. However, most investigators use transport models with just a single storage zone in their modeling studies, which assumes that the effects of multiple storage zones can be lumped together. Our study addressed the following question: Can a single‐storage zone model reliably characterize the effects of physical retention and biogeochemical reactions in multiple storage zones? We extended an existing stream transport model with a single storage zone to include a second storage zone. With the extended model we generated 500 data sets representing transport of nonreactive and reactive solutes in stream systems that have two different types of storage zones with variable hydrologic conditions. The one storage zone model was tested by optimizing the lumped storage parameters to achieve a best fit for each of the generated data sets. Multiple storage processes were categorized as possessing I, additive; II, competitive; or III, dominant storage zone characteristics. The classification was based on the goodness of fit of generated data sets, the degree of similarity in mean retention time of the two storage zones, and the relative distributions of exchange flux and storage capacity between the two storage zones. For most cases (>90%) the one storage zone model described either the effect of the sum of multiple storage processes (category I) or the dominant storage process (category III). Failure of the one storage zone model occurred mainly for category II, that is, when one of the storage zones had a much longer mean retention time (ts ratio > 5.0) and when the dominance of storage capacity and exchange flux occurred in different storage zones. We also used the one storage zone model to estimate a “single” lumped rate constant representing the net removal of a solute by biogeochemical reactions in multiple storage zones. For most cases the lumped rate constant that was optimized by one storage zone modeling estimated the flux‐weighted rate constant for multiple storage zones. Our results explain how the relative hydrologic properties of multiple storage zones (retention time, storage capacity, exchange flux, and biogeochemical reaction rate constant) affect the reliability of lumped parameters determined by a one storage zone transport model. We conclude that stream transport models with a single storage compartment will in most cases reliably characterize the dominant physical processes of solute retention and biogeochemical reactions in streams with multiple storage zones.

Uncertainty in Earth System Models: Benchmarks for Ocean Model Performance and Validation

NASA Astrophysics Data System (ADS)

Ogunro, O. O.; Elliott, S.; Collier, N.; Wingenter, O. W.; Deal, C.; Fu, W.; Hoffman, F. M.

2017-12-01

The mean ocean CO2 sink is a major component of the global carbon budget, with marine reservoirs holding about fifty times more carbon than the atmosphere. Phytoplankton play a significant role in the net carbon sink through photosynthesis and drawdown, such that about a quarter of anthropogenic CO2 emissions end up in the ocean. Biology greatly increases the efficiency of marine environments in CO2 uptake and ultimately reduces the impact of the persistent rise in atmospheric concentrations. However, a number of challenges remain in appropriate representation of marine biogeochemical processes in Earth System Models (ESM). These threaten to undermine the community effort to quantify seasonal to multidecadal variability in ocean uptake of atmospheric CO2. In a bid to improve analyses of marine contributions to climate-carbon cycle feedbacks, we have developed new analysis methods and biogeochemistry metrics as part of the International Ocean Model Benchmarking (IOMB) effort. Our intent is to meet the growing diagnostic and benchmarking needs of ocean biogeochemistry models. The resulting software package has been employed to validate DOE ocean biogeochemistry results by comparison with observational datasets. Several other international ocean models contributing results to the fifth phase of the Coupled Model Intercomparison Project (CMIP5) were analyzed simultaneously. Our comparisons suggest that the biogeochemical processes determining CO2 entry into the global ocean are not well represented in most ESMs. Polar regions continue to show notable biases in many critical biogeochemical and physical oceanographic variables. Some of these disparities could have first order impacts on the conversion of atmospheric CO2 to organic carbon. In addition, single forcing simulations show that the current ocean state can be partly explained by the uptake of anthropogenic emissions. Combined effects of two or more of these forcings on ocean biogeochemical cycles and ecosystems are challenging to predict since additive or antagonistic effects may occur. A benchmarking tool for accurate assessment and validation of marine biogeochemical outputs will be indispensable as the model community continues to improve ESM developments. It will provide a first order tool in understanding climate-carbon cycle feedbacks.

Marine Biogeochemistry Under The Influence of Fish And Fisheries: An Ecosystem Modeling Study

NASA Astrophysics Data System (ADS)

Disa, Deniz; Akoglu, Ekin; Salihoglu, Baris

2017-04-01

The ocean and the marine ecosystems are important controllers of the global carbon cycle. They play a pivotal role in capturing atmospheric carbon into the ocean body, transforming it into organic carbon through photosynthesis and transporting it to the depths of the ocean. Fish, which has a significant role in the marine food webs, is thought to have a considerable impact on carbon export. More specifically, fish has a control on plankton dynamics as a predator, it provides nutrient to the ecosystem by its metabolic activities and it has the ability of moving actively and transporting materials. Fishing is also expected to impact carbon cycle because it directly changes the fish biomasses. However, how fish impacts the biogeochemistry of marine ecosystems is not studied extensively. The aim of this study is to analyze the impact of fish and fisheries on marine biogeochemical processes by setting up an end-to-end model, which simulates lower and higher tropic levels of marine ecosystems simultaneously. For this purpose, a one dimensional biogeochemical model simulating lower tropic level dynamics (e.g. carbon export, nutrient cycles) and an food web model simulating fisheries exploitation and higher tropic level dynamics were online and two-way coupled. Representing the marine ecosystem from one end to the other, the coupled model served as a tool for the analysis of fishing impacts on marine biogeochemical dynamics. Results obtained after incorporation of higher trophic level model changed the plankton compositions and enhanced detritus pools and increased carbon export. Additionally, our model showed that active movement of fish contributed to transport of carbon from surface to the deeper parts of the ocean. Moreover, results after applying different fishing intensities indicated that changes in fisheries exploitation levels directly influence the marine nutrient cycles and hence, the carbon export. Depending on the target and the intensity of fisheries, considerable changes in the biogeochemical responses observed. In conclusion, unlike the models that do not represent the fish explicitly, we demonstrate how marine biogeochemical processes are impacted by the activity of fish assemblages and fisheries exploitation.

Predicting Biological Information Flow in a Model Oxygen Minimum Zone

NASA Astrophysics Data System (ADS)

Louca, S.; Hawley, A. K.; Katsev, S.; Beltran, M. T.; Bhatia, M. P.; Michiels, C.; Capelle, D.; Lavik, G.; Doebeli, M.; Crowe, S.; Hallam, S. J.

2016-02-01

Microbial activity drives marine biochemical fluxes and nutrient cycling at global scales. Geochemical measurements as well as molecular techniques such as metagenomics, metatranscriptomics and metaproteomics provide great insight into microbial activity. However, an integration of molecular and geochemical data into mechanistic biogeochemical models is still lacking. Recent work suggests that microbial metabolic pathways are, at the ecosystem level, strongly shaped by stoichiometric and energetic constraints. Hence, models rooted in fluxes of matter and energy may yield a holistic understanding of biogeochemistry. Furthermore, such pathway-centric models would allow a direct consolidation with meta'omic data. Here we present a pathway-centric biogeochemical model for the seasonal oxygen minimum zone in Saanich Inlet, a fjord off the coast of Vancouver Island. The model considers key dissimilatory nitrogen and sulfur fluxes, as well as the population dynamics of the genes that mediate them. By assuming a direct translation of biocatalyzed energy fluxes to biosynthesis rates, we make predictions about the distribution and activity of the corresponding genes. A comparison of the model to molecular measurements indicates that the model explains observed DNA, RNA, protein and cell depth profiles. This suggests that microbial activity in marine ecosystems such as oxygen minimum zones is well described by DNA abundance, which, in conjunction with geochemical constraints, determines pathway expression and process rates. Our work further demonstrates how meta'omic data can be mechanistically linked to environmental redox conditions and biogeochemical processes.

Residence-time framework for modeling multicomponent reactive transport in stream hyporheic zones

NASA Astrophysics Data System (ADS)

Painter, S. L.; Coon, E. T.; Brooks, S. C.

2017-12-01

Process-based models for transport and transformation of nutrients and contaminants in streams require tractable representations of solute exchange between the stream channel and biogeochemically active hyporheic zones. Residence-time based formulations provide an alternative to detailed three-dimensional simulations and have had good success in representing hyporheic exchange of non-reacting solutes. We extend the residence-time formulation for hyporheic transport to accommodate general multicomponent reactive transport. To that end, the integro-differential form of previous residence time models is replaced by an equivalent formulation based on a one-dimensional advection dispersion equation along the channel coupled at each channel location to a one-dimensional transport model in Lagrangian travel-time form. With the channel discretized for numerical solution, the associated Lagrangian model becomes a subgrid model representing an ensemble of streamlines that are diverted into the hyporheic zone before returning to the channel. In contrast to the previous integro-differential forms of the residence-time based models, the hyporheic flowpaths have semi-explicit spatial representation (parameterized by travel time), thus allowing coupling to general biogeochemical models. The approach has been implemented as a stream-corridor subgrid model in the open-source integrated surface/subsurface modeling software ATS. We use bedform-driven flow coupled to a biogeochemical model with explicit microbial biomass dynamics as an example to show that the subgrid representation is able to represent redox zonation in sediments and resulting effects on metal biogeochemical dynamics in a tractable manner that can be scaled to reach scales.

Global Biology Research Program: Biogeochemical Processes in Wetlands

NASA Technical Reports Server (NTRS)

Bartlett, D. S. (Editor)

1984-01-01

The results of a workshop examining potential NASA contributions to research on wetland processes as they relate to global biogeochemical cycles are summarized. A wetlands data base utilizing remotely sensed inventories, studies of wetland/atmosphere exchange processes, and the extrapolation of local measurements to global biogeochemical cycling processes were identified as possible areas for NASA support.

Constructing wetlands: measuring and modeling feedbacks of oxidation processes between plants and clay-rich material

NASA Astrophysics Data System (ADS)

Saaltink, Rémon; Dekker, Stefan C.; Griffioen, Jasper; Wassen, Martin J.

2016-04-01

Interest is growing in using soft sediment as a building material in eco-engineering projects. Wetland construction in the Dutch lake Markermeer is an example: here the option of dredging some of the clay-rich lake-bed sediment and using it to construct 10.000 ha of wetland will soon go under construction. Natural processes will be utilized during and after construction to accelerate ecosystem development. Knowing that plants can eco-engineer their environment via positive or negative biogeochemical plant-soil feedbacks, we conducted a six-month greenhouse experiment to identify the key biogeochemical processes in the mud when Phragmites australis is used as an eco-engineering species. We applied inverse biogeochemical modeling to link observed changes in pore water composition to biogeochemical processes. Two months after transplantation we observed reduced plant growth and shriveling as well as yellowing of foliage. The N:P ratios of plant tissue were low and were affected not by hampered uptake of N but by enhanced uptake of P. Plant analyses revealed high Fe concentrations in the leaves and roots. Sulfate concentrations rose drastically in our experiment due to pyrite oxidation; as reduction of sulfate will decouple Fe-P in reducing conditions, we argue that plant-induced iron toxicity hampered plant growth, forming a negative feedback loop, while simultaneously there was a positive feedback loop, as iron toxicity promotes P mobilization as a result of reduced conditions through root death, thereby stimulating plant growth and regeneration. Given these two feedback mechanisms, we propose that when building wetlands from these mud deposits Fe-tolerant species are used rather than species that thrive in N-limited conditions. The results presented in this study demonstrate the importance of studying the biogeochemical properties of the building material and the feedback mechanisms between plant and soil prior to finalizing the design of the eco-engineering project.

Error assessment of biogeochemical models by lower bound methods (NOMMA-1.0)

NASA Astrophysics Data System (ADS)

Sauerland, Volkmar; Löptien, Ulrike; Leonhard, Claudine; Oschlies, Andreas; Srivastav, Anand

2018-03-01

Biogeochemical models, capturing the major feedbacks of the pelagic ecosystem of the world ocean, are today often embedded into Earth system models which are increasingly used for decision making regarding climate policies. These models contain poorly constrained parameters (e.g., maximum phytoplankton growth rate), which are typically adjusted until the model shows reasonable behavior. Systematic approaches determine these parameters by minimizing the misfit between the model and observational data. In most common model approaches, however, the underlying functions mimicking the biogeochemical processes are nonlinear and non-convex. Thus, systematic optimization algorithms are likely to get trapped in local minima and might lead to non-optimal results. To judge the quality of an obtained parameter estimate, we propose determining a preferably large lower bound for the global optimum that is relatively easy to obtain and that will help to assess the quality of an optimum, generated by an optimization algorithm. Due to the unavoidable noise component in all observations, such a lower bound is typically larger than zero. We suggest deriving such lower bounds based on typical properties of biogeochemical models (e.g., a limited number of extremes and a bounded time derivative). We illustrate the applicability of the method with two real-world examples. The first example uses real-world observations of the Baltic Sea in a box model setup. The second example considers a three-dimensional coupled ocean circulation model in combination with satellite chlorophyll a.

Modeling the Oxygen Cycle in the Equatorial Pacific: Regulation of Physical and Biogeochemical Processes

NASA Astrophysics Data System (ADS)

Wang, X.; Murtugudde, R. G.; Zhang, D.

2016-12-01

Photosynthesis and respiration are important processes in all ecosystems on the Earth, in which carbon and oxygen are the two main elements. However, the oxygen cycle has received much less attention (relative to the carbon cycle) despite its big role in the earth system. Oxygen is a sensitive indicator of physical and biogeochemical processes in the ocean thus a key parameter for understanding the ocean's ecosystem and biogeochemistry. The Oxygen-Minimum-Zone (OMZ), often seen below 200 m, is a profound feature in the world oceans. There has been evidence of OMZ expansion over the past few decades in the tropical oceans. Climate models project that there would be a continued decline in dissolved oxygen (DO) and an expansion of the tropical OMZs under future warming conditions, which is of great concern because of the implications for marine organisms. We employ a validated three-dimensional model that simulates physical transport (circulation and vertical mixing), biological processes (O2 production and consumption) and ocean-atmosphere O2 exchange to quantify various sources and sinks of DO over 1980-2015. We show how we use observational data to improve our model simulation. Then we assess the spatial and temporal variability in simulated DO in the tropical Pacific Ocean, and explore the impacts of physical and biogeochemical processes on the DO dynamics, with a focus on the MOZ. Our analyses indicate that DO in the OMZ has a positive relationship with the 13ºC isotherm depth and a negative relationship with the concentration of dissolved organic material.

PEATBOG: a biogeochemical model for analyzing coupled carbon and nitrogen dynamics in northern peatlands

NASA Astrophysics Data System (ADS)

Wu, Y.; Blodau, C.

2013-08-01

Elevated nitrogen deposition and climate change alter the vegetation communities and carbon (C) and nitrogen (N) cycling in peatlands. To address this issue we developed a new process-oriented biogeochemical model (PEATBOG) for analyzing coupled carbon and nitrogen dynamics in northern peatlands. The model consists of four submodels, which simulate: (1) daily water table depth and depth profiles of soil moisture, temperature and oxygen levels; (2) competition among three plants functional types (PFTs), production and litter production of plants; (3) decomposition of peat; and (4) production, consumption, diffusion and export of dissolved C and N species in soil water. The model is novel in the integration of the C and N cycles, the explicit spatial resolution belowground, the consistent conceptualization of movement of water and solutes, the incorporation of stoichiometric controls on elemental fluxes and a consistent conceptualization of C and N reactivity in vegetation and soil organic matter. The model was evaluated for the Mer Bleue Bog, near Ottawa, Ontario, with regards to simulation of soil moisture and temperature and the most important processes in the C and N cycles. Model sensitivity was tested for nitrogen input, precipitation, and temperature, and the choices of the most uncertain parameters were justified. A simulation of nitrogen deposition over 40 yr demonstrates the advantages of the PEATBOG model in tracking biogeochemical effects and vegetation change in the ecosystem.

PEATBOG: a biogeochemical model for analyzing coupled carbon and nitrogen dynamics in northern peatlands

NASA Astrophysics Data System (ADS)

Wu, Y.; Blodau, C.

2013-03-01

Elevated nitrogen deposition and climate change alter the vegetation communities and carbon (C) and nitrogen (N) cycling in peatlands. To address this issue we developed a new process-oriented biogeochemical model (PEATBOG) for analyzing coupled carbon and nitrogen dynamics in northern peatlands. The model consists of four submodels, which simulate: (1) daily water table depth and depth profiles of soil moisture, temperature and oxygen levels; (2) competition among three plants functional types (PFTs), production and litter production of plants; (3) decomposition of peat; and (4) production, consumption, diffusion and export of dissolved C and N species in soil water. The model is novel in the integration of the C and N cycles, the explicit spatial resolution belowground, the consistent conceptualization of movement of water and solutes, the incorporation of stoichiometric controls on elemental fluxes and a consistent conceptualization of C and N reactivity in vegetation and soil organic matter. The model was evaluated for the Mer Bleue Bog, near Ottawa, Ontario, with regards to simulation of soil moisture and temperature and the most important processes in the C and N cycles. Model sensitivity was tested for nitrogen input, precipitation, and temperature, and the choices of the most uncertain parameters were justified. A simulation of nitrogen deposition over 40 yr demonstrates the advantages of the PEATBOG model in tracking biogeochemical effects and vegetation change in the ecosystem.

Biogeochemical carbon coupling influences global precipitation in geoengineering experiments

NASA Astrophysics Data System (ADS)

Fyfe, J. C.; Cole, J. N. S.; Arora, V. K.; Scinocca, J. F.

2013-02-01