Abiotic/biotic coupling in the rhizosphere: a reactive transport modeling analysis

Science.gov (United States)

Lawrence, Corey R.; Steefel, Carl; Maher, Kate

2014-01-01

A new generation of models is needed to adequately simulate patterns of soil biogeochemical cycling in response changing global environmental drivers. For example, predicting the influence of climate change on soil organic matter storage and stability requires models capable of addressing complex biotic/abiotic interactions of rhizosphere and weathering processes. Reactive transport modeling provides a powerful framework simulating these interactions and the resulting influence on soil physical and chemical characteristics. Incorporation of organic reactions in an existing reactive transport model framework has yielded novel insights into soil weathering and development but much more work is required to adequately capture root and microbial dynamics in the rhizosphere. This endeavor provides many advantages over traditional soil biogeochemical models but also many challenges.

Incorporating nitrogen fixing cyanobacteria in the global biogeochemical model HAMOCC

Science.gov (United States)

Paulsen, Hanna; Ilyina, Tatiana; Six, Katharina

2015-04-01

Nitrogen fixation by marine diazotrophs plays a fundamental role in the oceanic nitrogen and carbon cycle as it provides a major source of 'new' nitrogen to the euphotic zone that supports biological carbon export and sequestration. Since most global biogeochemical models include nitrogen fixation only diagnostically, they are not able to capture its spatial pattern sufficiently. Here we present the incorporation of an explicit, dynamic representation of diazotrophic cyanobacteria and the corresponding nitrogen fixation in the global ocean biogeochemical model HAMOCC (Hamburg Ocean Carbon Cycle model), which is part of the Max Planck Institute for Meteorology Earth system model (MPI-ESM). The parameterization of the diazotrophic growth is thereby based on available knowledge about the cyanobacterium Trichodesmium spp., which is considered as the most significant pelagic nitrogen fixer. Evaluation against observations shows that the model successfully reproduces the main spatial distribution of cyanobacteria and nitrogen fixation, covering large parts of the tropical and subtropical oceans. Besides the role of cyanobacteria in marine biogeochemical cycles, their capacity to form extensive surface blooms induces a number of bio-physical feedback mechanisms in the Earth system. The processes driving these interactions, which are related to the alteration of heat absorption, surface albedo and momentum input by wind, are incorporated in the biogeochemical and physical model of the MPI-ESM in order to investigate their impacts on a global scale. First preliminary results will be shown.

Biogeochemical consequences of vertical and lateral transport of particulate organic matter in the southern North Sea: A multiproxy approach

NARCIS (Netherlands)

Le Guitton, M.; Soetaert, K.; Sinninghe Damsté, J.S.; Middelburg, J.J.

2015-01-01

Vertical and lateral transports are of importance in continental shelf systems such as the North Sea andplay a major role in the processing of organic matter. We investigated the biogeochemical consequencesof these transports on particulate organic matter at the molecular level in the southern North

Linking Chaotic Advection with Subsurface Biogeochemical Processes

Science.gov (United States)

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.

Biogeochemical consequences of vertical and lateral transport of particulate organic matter in the southern North Sea: A multiproxy approach

NARCIS (Netherlands)

le Guitton, M.; Soetaert, K.; Sinninghe Damsté, J.S.; Middelburg, J.J.

2015-01-01

Vertical and lateral transports are of importance in continental shelf systems such as the North Sea and play a major role in the processing of organic matter. We investigated the biogeochemical consequences of these transports on particulate organic matter at the molecular level in the southern

Impact of urban effluents on summer hypoxia in the highly turbid Gironde Estuary, applying a 3D model coupling hydrodynamics, sediment transport and biogeochemical processes

Science.gov (United States)

Lajaunie-Salla, Katixa; Wild-Allen, Karen; Sottolichio, Aldo; Thouvenin, Bénédicte; Litrico, Xavier; Abril, Gwenaël

2017-10-01

Estuaries are increasingly degraded due to coastal urban development and are prone to hypoxia problems. The macro-tidal Gironde Estuary is characterized by a highly concentrated turbidity maximum zone (TMZ). Field observations show that hypoxia occurs in summer in the TMZ at low river flow and a few days after the spring tide peak. In situ data highlight lower dissolved oxygen (DO) concentrations around the city of Bordeaux, located in the upper estuary. Interactions between multiple factors limit the understanding of the processes controlling the dynamics of hypoxia. A 3D biogeochemical model was developed, coupled with hydrodynamics and a sediment transport model, to assess the contribution of the TMZ and the impact of urban effluents through wastewater treatment plants (WWTPs) and sewage overflows (SOs) on hypoxia. Our model describes the transport of solutes and suspended material and the biogeochemical mechanisms impacting oxygen: primary production, degradation of all organic matter (i.e. including phytoplankton respiration, degradation of river and urban watershed matter), nitrification and gas exchange. The composition and the degradation rates of each variable were characterized by in situ measurements and experimental data from the study area. The DO model was validated against observations in Bordeaux City. The simulated DO concentrations show good agreement with field observations and satisfactorily reproduce the seasonal and neap-spring time scale variations around the city of Bordeaux. Simulations show a spatial and temporal correlation between the formation of summer hypoxia and the location of the TMZ, with minimum DO centered in the vicinity of Bordeaux. To understand the contribution of the urban watershed forcing, different simulations with the presence or absence of urban effluents were compared. Our results show that in summer, a reduction of POC from SO would increase the DO minimum in the vicinity of Bordeaux by 3% of saturation. Omitting

Biogeochemical Impact of Long-Range Transported Dust over Northern South China Sea

Science.gov (United States)

Tsay, Si-Chee; Wang, S. H.; Hsu, N. C.

2011-01-01

Transpacific transport and impact of Asian dust aerosols have been well documented (e.g., results from ACE-Asia and regional follow-on campaigns), but little is known about dust invasion to the South China Sea (SCS). On 19-21 March 2010, a fierce Asian dust storm affected large areas from the Gobi deserts to the West Pacific, including Taiwan and Hong Kong. As a pilot study of the 7-SEAS (Seven South East Asian Studies) in the northern SCS, detailed characteristics of long-range transported dust aerosols were first observed by a comprehensive set of ground-based instruments deployed at the Dongsha islands (20deg42'52" N, 116deg43'51" E). Aerosol measurements such as particle mass concentrations, size distribution, optical properties, hygroscopicity, and vertical profiles help illustrate the evolution of this dust outbreak. Our results indicate that these dust particles were mixed with anthropogenic and marine aerosols, and transported near the surface. Satellite assessment of biogeochemical impact of dust deposition into open oceans is hindered by our current inability in retrieving areal dust properties and ocean colors over an extensive period of time, particularly under the influence of cloudy conditions. In this paper, we analyze the changes of retrieved Chlorophyll-a (Chl-a) concentration over the northern SCS, considered as oligotophic waters in the spring, from long-term SeaWiFS measurements since 1997. Over the past decade, six long-range transported dust events are identified based on spatiotemporal evolutions of PM10 measurements from regional monitoring stations, with the aid of trajectory analysis. Multi-year composites of Chl-a imagery for dust event and non-dust background during March-April are applied to overcome insufficient retrievals of Chl-a due to cloudy environment. Due to anthropogenic modification within a shallow boundary layer off the densely populated and industrial southeast coast of China, the iron ion activation of deliquescent dust

A biogeochemical transport model to simulate the attenuation of chlorinated hydrocarbon contaminant fluxes across the groundwater-surface water interface

DEFF Research Database (Denmark)

Malaguerra, Flavio; Binning, Philip John; Albrechtsen, Hans-Jørgen

2009-01-01

and biogeochemical transformation model of the discharge of a TCE plume into a stream, and to determine which parameters most strongly affect pollutant discharge concentrations. Here biological kinetics and the interaction with the soil matrix are implemented in PHREEQC. The ability of PHREEQC to deal with a large...

Physical/biogeochemical coupled model : impact of an offline vs online strategy

Science.gov (United States)

Hameau, Angélique; Perruche, Coralie; Bricaud, Clément; Gutknecht, Elodie; Reffray, Guillaume

2014-05-01

Mercator-Ocean, the French ocean forecasting center, has been developing several operational forecasting systems and reanalysis of the physical and biogeochemical 3D-Ocean. Here we study the impact of an offline vs online strategy to couple the physical (OPA) and biogeochemical (PISCES) modules included in the NEMO platform. For this purpose, we perform global one-year long simulations at 1° resolution. The model was initialized with global climatologies. The spin-up involved 10 years of biogeochemical off-line simulation forced by a climatology of ocean physics. The online mode consists in running physical and biogeochemical models simultaneously whereas in the offline mode, the biogeochemical model is launched alone, forced by averaged physical forcing (1 day, 7 days,… ). The Mercator operational biogeochemical system is currently using the offline mode with a weekly physical forcing. A special treatment is applied to the vertical diffusivity coefficient (Kz): as it varies of several orders of magnitude, we compute the mean of the LOG10 of Kz. Moreover, a threshold value is applied to remove the highest values corresponding to enhanced convection. To improve this system, 2 directions are explored. First, 3 physical forcing frequencies are compared to quantify errors due to the offline mode: 1 hour (online mode), 1 day and 1 week (offline modes). Secondly, sensitivity tests to the threshold value applied to Kz are performed. The simulations are evaluated by systematically comparing model fields to observations (Globcolour product and World Ocean Atlas 2005) at global and regional scales. We show first that offline simulations are in good agreement with online simulation. As expected, the lower the physical forcing frequency is, the closer to the online solution is the offline simulation. The threshold value on the vertical diffusivity coefficient manages the mixing strength within the mixed layer. A value of 1 m2.s-1 appears to be a good compromise to approach

Dynamic modeling of nitrogen losses in river networks unravels the coupled effects of hydrological and biogeochemical processes

Science.gov (United States)

Alexander, Richard B.; Böhlke, John Karl; Boyer, Elizabeth W.; David, Mark B.; Harvey, Judson W.; Mulholland, Patrick J.; Seitzinger, Sybil P.; Tobias, Craig R.; Tonitto, Christina; Wollheim, Wilfred M.

2009-01-01

The importance of lotic systems as sinks for nitrogen inputs is well recognized. A fraction of nitrogen in streamflow is removed to the atmosphere via denitrification with the remainder exported in streamflow as nitrogen loads. At the watershed scale, there is a keen interest in understanding the factors that control the fate of nitrogen throughout the stream channel network, with particular attention to the processes that deliver large nitrogen loads to sensitive coastal ecosystems. We use a dynamic stream transport model to assess biogeochemical (nitrate loadings, concentration, temperature) and hydrological (discharge, depth, velocity) effects on reach-scale denitrification and nitrate removal in the river networks of two watersheds having widely differing levels of nitrate enrichment but nearly identical discharges. Stream denitrification is estimated by regression as a nonlinear function of nitrate concentration, streamflow, and temperature, using more than 300 published measurements from a variety of US streams. These relations are used in the stream transport model to characterize nitrate dynamics related to denitrification at a monthly time scale in the stream reaches of the two watersheds. Results indicate that the nitrate removal efficiency of streams, as measured by the percentage of the stream nitrate flux removed via denitrification per unit length of channel, is appreciably reduced during months with high discharge and nitrate flux and increases during months of low-discharge and flux. Biogeochemical factors, including land use, nitrate inputs, and stream concentrations, are a major control on reach-scale denitrification, evidenced by the disproportionately lower nitrate removal efficiency in streams of the highly nitrate-enriched watershed as compared with that in similarly sized streams in the less nitrate-enriched watershed. Sensitivity analyses reveal that these important biogeochemical factors and physical hydrological factors contribute nearly

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

Science.gov (United States)

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

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

Directory of Open Access Journals (Sweden)

V. Sauerland

2018-03-01

Full Text Available 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.

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

Science.gov (United States)

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.

Development of an advanced eco-hydrologic and biogeochemical coupling model aimed at clarifying the missing role of inland water in the global biogeochemical cycle

Science.gov (United States)

Nakayama, Tadanobu

2017-04-01

Recent research showed that inland water including rivers, lakes, and groundwater may play some role in carbon cycling, although its contribution has remained uncertain due to limited amount of reliable data available. In this study, the author developed an advanced model coupling eco-hydrology and biogeochemical cycle (National Integrated Catchment-based Eco-hydrology (NICE)-BGC). This new model incorporates complex coupling of hydrologic-carbon cycle in terrestrial-aquatic linkages and interplay between inorganic and organic carbon during the whole process of carbon cycling. The model could simulate both horizontal transports (export from land to inland water 2.01 ± 1.98 Pg C/yr and transported to ocean 1.13 ± 0.50 Pg C/yr) and vertical fluxes (degassing 0.79 ± 0.38 Pg C/yr, and sediment storage 0.20 ± 0.09 Pg C/yr) in major rivers in good agreement with previous researches, which was an improved estimate of carbon flux from previous studies. The model results also showed global net land flux simulated by NICE-BGC (-1.05 ± 0.62 Pg C/yr) decreased carbon sink a little in comparison with revised Lund-Potsdam-Jena Wetland Hydrology and Methane (-1.79 ± 0.64 Pg C/yr) and previous materials (-2.8 to -1.4 Pg C/yr). This is attributable to CO2 evasion and lateral carbon transport explicitly included in the model, and the result suggests that most previous researches have generally overestimated the accumulation of terrestrial carbon and underestimated the potential for lateral transport. The results further implied difference between inverse techniques and budget estimates suggested can be explained to some extent by a net source from inland water. NICE-BGC would play an important role in reevaluation of greenhouse gas budget of the biosphere, quantification of hot spots, and bridging the gap between top-down and bottom-up approaches to global carbon budget.

Multiscale Investigation on Biofilm Distribution and Its Impact on Macroscopic Biogeochemical Reaction Rates

Science.gov (United States)

Yan, Zhifeng; Liu, Chongxuan; Liu, Yuanyuan; Bailey, Vanessa L.

2017-11-01

Biofilms are critical locations for biogeochemical reactions in the subsurface environment. The occurrence and distribution of biofilms at microscale as well as their impacts on macroscopic biogeochemical reaction rates are still poorly understood. This paper investigated the formation and distributions of biofilms in heterogeneous sediments using multiscale models and evaluated the effects of biofilm heterogeneity on local and macroscopic biogeochemical reaction rates. Sediment pore structures derived from X-ray computed tomography were used to simulate the microscale flow dynamics and biofilm distribution in the sediment column. The response of biofilm formation and distribution to the variations in hydraulic and chemical properties was first examined. One representative biofilm distribution was then utilized to evaluate its effects on macroscopic reaction rates using nitrate reduction as an example. The results revealed that microorganisms primarily grew on the surfaces of grains and aggregates near preferential flow paths where both electron donor and acceptor were readily accessible, leading to the heterogeneous distribution of biofilms in the sediments. The heterogeneous biofilm distribution decreased the macroscopic rate of biogeochemical reactions as compared with those in homogeneous cases. Operationally considering the heterogeneous biofilm distribution in macroscopic reactive transport models such as using dual porosity domain concept can significantly improve the prediction of biogeochemical reaction rates. Overall, this study provided important insights into the biofilm formation and distribution in soils and sediments as well as their impacts on the macroscopic manifestation of reaction rates.

Nitrous Oxide Emissions from Biofuel Crops and Parameterization in the EPIC Biogeochemical Model

Science.gov (United States)

This presentation describes year 1 field measurements of N2O fluxes and crop yields which are used to parameterize the EPIC biogeochemical model for the corresponding field site. Initial model simulations are also presented.

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

Science.gov (United States)

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.

A mercury transport and fate model (LM2-mercury) for mass budget assessment of mercury cycling in Lake Michigan

Science.gov (United States)

LM2-Mercury, a mercury mass balance model, was developed to simulate and evaluate the transport, fate, and biogeochemical transformations of mercury in Lake Michigan. The model simulates total suspended solids (TSS), disolved organic carbon (DOC), and total, elemental, divalent, ...

Extracellular Electron Transport Coupling Biogeochemical Processes Centimeters

DEFF Research Database (Denmark)

Risgaard-Petersen, Nils; Fossing, Henrik; Christensen, Peter Bondo

2010-01-01

of the oxygen uptake in laboratory incubations of initially homogenized and stabilized sediment. Using microsensors and process rate measurements we further investigated the effect of the electric currents on sediment biogeochemistry. Dissolved sulfide readily donated electrons to the networks and could...... confirmed the depth range of the electric communication and indicated donation of electrons directly from organotrophic bacteria. The separation of oxidation and reduction processes created steep pH gradients eventually causing carbonate precipitation at the surface. The results indicate that electron...... exchanging organisms have major biogeochemical importance as they allow widely separated electron donors and acceptors to react with one another....

Possible impacts of global warming on tundra and boreal forest ecosystems - comparison of some biogeochemical models

Energy Technology Data Exchange (ETDEWEB)

Ploechl, M.; Cramer, W.

1995-06-01

Global warming affects the magnitude of carbon, water and nitrogen fluxes between biosphere and atmosphere as well as the distribution of vegetation types. Biogeochemical models, global as well as patch models, can be used to estimate the differences between the mean values of annual net primary production (NPP) for the present and for future climate scenarios. Both approaches rely on the prescribed pattern of vegetation types. Structural, rule based models can predict such patterns, provided that vegetation and climate are in equilibrium. The coupling of biogeochemical and structural models gives the opportunity to test the sensitivity of biogeochemical processes not only to climatic change but also to biome shifts. Whether the annual mean NPP of a vegetation type increses or decreases depends strongly on the assumptions about a CO{sub 2} fertilization effect and nitrogen cycling. Results from our coupled model show that, given that direct CO{sub 2} effects are uncertain, (i) average NPP of these northern biomes might decrease under global warming, but (ii) total NPP of the region would increase, due to the northward shift of the taiga biome. (orig.)

Modeling multi-component transport and enhanced anaerobic dechlorination processes in a single fracture-clay matrix system

DEFF Research Database (Denmark)

Chambon, Julie Claire Claudia; Broholm, Mette Martina; Binning, Philip John

2010-01-01

Clayey tills contaminated with chlorinated solvents are a threat to groundwater and are difficult to remediate. A numerical model is developed for assessing leaching processes and for simulating the remediation via enhanced anaerobic dechlorination. The model simulates the transport...... to the physical processes, mainly diffusion in the matrix, than to the biogeochemical processes, when dechlorination is assumed to take place in a limited reaction zone only. The inclusion of sequential dechlorination in clay fracture transport models is crucial, as the contaminant flux to the aquifer...

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

Science.gov (United States)

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

Simulation of glacial ocean biogeochemical tracer and isotope distributions based on the PMIP3 suite of climate models

Science.gov (United States)

Khatiwala, Samar; Muglia, Juan; Kvale, Karin; Schmittner, Andreas

2016-04-01

In the present climate system, buoyancy forced convection at high-latitudes together with internal mixing results in a vigorous overturning circulation whose major component is North Atlantic Deep Water. One of the key questions of climate science is whether this "mode" of circulation persisted during glacial periods, and in particular at the Last Glacial Maximum (LGM; 21000 years before present). Resolving this question is both important for advancing our understanding of the climate system, as well as a critical test of numerical models' ability to reliably simulate different climates. The observational evidence, based on interpreting geochemical tracers archived in sediments, is conflicting, as are simulations carried out with state-of-the-art climate models (e.g., as part of the PMIP3 suite), which, due to the computational cost involved, do not by and large include biogeochemical and isotope tracers that can be directly compared with proxy data. Here, we apply geochemical observations to evaluate the ability of several realisations of an ocean model driven by atmospheric forcing from the PMIP3 suite of climate models to simulate global ocean circulation during the LGM. This results in a wide range of circulation states that are then used to simulate biogeochemical tracer and isotope (13C, 14C and Pa/Th) distributions using an efficient, "offline" computational scheme known as the transport matrix method (TMM). One of the key advantages of this approach is the use of a uniform set of biogeochemical and isotope parameterizations across all the different circulations based on the PMIP3 models. We compare these simulated distributions to both modern observations and data from LGM ocean sediments to identify similarities and discrepancies between model and data. We find, for example, that when the ocean model is forced with wind stress from the PMIP3 models the radiocarbon age of the deep ocean is systematically younger compared with reconstructions. Changes in

Using a Hydrodynamic and Biogeochemical Model to Investigate the Effects of Nutrient Loading from a Wastewater Treatment Plant into Lake Michigan

Science.gov (United States)

Khazaei, B.; Bravo, H.; Bootsma, H.

2017-12-01

There is clear evidence that excessive nutrient, in particular phosphorus (P), loading into Lake Michigan has produced significant problems, such as algal blooms, hypoxia, and reduced water quality. Addressing those problems requires understanding the transport and fate of P in the lake. The dominance of mixing and dispersion processes on the P transport has been demonstrated, yet recent research has shown the remarkable influence of dreissenid mussels and Cladophora on water clarity and the P budget. Since mussels and Cladophora tend to concentrate near the coastlines, nearshore-offshore P exchange is of a big importance. In this research, a computer model was developed to simulate the P cycle by incorporating the biogeochemical processes relevant to the transport of P into a 3D high-resolution hydrodynamic model. The near-bottom biogeochemical model consists of three linked modules: Cladophora, mussel, and sediment storage modules. The model was applied to the Milwaukee Metropolitan Sewerage District South Shore Wastewater Treatment Plant, between June and October of 2013 and 2015, as a case study. The plant outfall introduces a point source of P into the study domain—the nearshore zone of Lake Michigan adjacent to Milwaukee County. The model was validated against field observations of water temperature, dissolved phosphorus (DP), particulate phosphorus (PP), Cladophora biomass, and P content. The model simulations showed reasonably good agreement with field measurements. Model results showed a) different temporal patterns in 2013 and 2015, b) a larger range of fluctuations in DP than that in PP, and c) that the effects of mussels and Cladophora could explain the differences in patterns and ranges. PP concentrations showed more frequent spikes of concentration in 2013 due to resuspension events during that year because of stronger winds. The model is being applied as a management tool to test scenarios of nutrient loading to determine effluent P limits for the

Using a spatially-distributed hydrologic biogeochemistry model with nitrogen transport to study the spatial variation of carbon stocks and fluxes in a Critical Zone Observatory

Science.gov (United States)

Shi, Y.; Eissenstat, D. M.; He, Y.; Davis, K. J.

2017-12-01

Most current biogeochemical models are 1-D and represent one point in space. Therefore, they cannot resolve topographically driven land surface heterogeneity (e.g., lateral water flow, soil moisture, soil temperature, solar radiation) or the spatial pattern of nutrient availability. A spatially distributed forest biogeochemical model with nitrogen transport, Flux-PIHM-BGC, has been developed by coupling a 1-D mechanistic biogeochemical model Biome-BGC (BBGC) with a spatially distributed land surface hydrologic model, Flux-PIHM, and adding an advection dominated nitrogen transport module. Flux-PIHM is a coupled physically based 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, and is augmented by adding a topographic solar radiation module. Flux-PIHM is able to represent the link between groundwater and the surface energy balance, as well as land surface heterogeneities caused by topography. In the coupled Flux-PIHM-BGC model, each Flux-PIHM model grid couples a 1-D BBGC model, while nitrogen is transported among model grids via surface and subsurface water flow. In each grid, Flux-PIHM provides BBGC with soil moisture, soil temperature, and solar radiation, while BBGC provides Flux-PIHM with spatially-distributed leaf area index. The coupled Flux-PIHM-BGC model has been implemented at the Susquehanna/Shale Hills Critical Zone Observatory. The model-predicted aboveground vegetation carbon and soil carbon distributions generally agree with the macro patterns observed within the watershed. The importance of abiotic variables (including soil moisture, soil temperature, solar radiation, and soil mineral nitrogen) in predicting aboveground carbon distribution is calculated using a random forest. The result suggests that the spatial pattern of aboveground carbon is controlled by the distribution of soil mineral nitrogen. A Flux-PIHM-BGC simulation

A novel modeling tool with multi-stressor functionality for organic contaminant transport and fate in the Baltic Sea

Energy Technology Data Exchange (ETDEWEB)

Undeman, E., E-mail: emma.undeman@itm.su.se [Baltic Nest Institute, Baltic Sea Centre, Stockholm University, 10691 Stockholm (Sweden); Department of Applied Environmental Science, Stockholm University, 11418 Stockholm (Sweden); Gustafsson, E., E-mail: erik.gustafsson@su.se [Baltic Nest Institute, Baltic Sea Centre, Stockholm University, 10691 Stockholm (Sweden); Gustafsson, B.G., E-mail: bo.gustafsson@su.se [Baltic Nest Institute, Baltic Sea Centre, Stockholm University, 10691 Stockholm (Sweden)

2014-11-01

The coupled physical–biogeochemical model BALTSEM, previously used to assess nutrient/carbon cycles and eutrophication in the Baltic Sea, has been expanded to include algorithms for calculations of organic contaminant environmental transport and fate. This novel model version (BALTSEM-POP) is evaluated for polychlorinated biphenyls (PCBs), polychlorinated dibenzo-p-dioxins/dibenzofurans (PCDD/Fs) and hexachlorobenzene (HCB) in Baltic Sea surface water and sediment. Modeled dissolved concentrations are usually within a factor of 2–4 of observed concentrations, however with larger deviations for furans. Calculated concentrations in particulate organic matter are less accurate (within factors of 1–700), likely due to errors in estimated pelagic biomass, particulate matter–water partitioning, and large natural variability in field data. Concentrations in sediments are usually predicted within a factor of 6. The good performance of the model illustrates its usefulness for exploration of contaminant fate in response to variations in nutrient input and climatic conditions in the Baltic Sea marine environment. - Highlights: • A new model for organic chemical transport and fate in the Baltic Sea is presented. • Physical and biogeochemical processes are linked to organic contaminant transport. • The model is evaluated for PCBs, HCB and PCDD/Fs. • The model can predict dissolved concentrations within a factor of ca 2–4. • Predictions for concentrations in particulate matter and sediment are less accurate.

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

Science.gov (United States)

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.

A Thermodynamically-consistent FBA-based Approach to Biogeochemical Reaction Modeling

Science.gov (United States)

Shapiro, B.; Jin, Q.

2015-12-01

Microbial rates are critical to understanding biogeochemical processes in natural environments. Recently, flux balance analysis (FBA) has been applied to predict microbial rates in aquifers and other settings. FBA is a genome-scale constraint-based modeling approach that computes metabolic rates and other phenotypes of microorganisms. This approach requires a prior knowledge of substrate uptake rates, which is not available for most natural microbes. Here we propose to constrain substrate uptake rates on the basis of microbial kinetics. Specifically, we calculate rates of respiration (and fermentation) using a revised Monod equation; this equation accounts for both the kinetics and thermodynamics of microbial catabolism. Substrate uptake rates are then computed from the rates of respiration, and applied to FBA to predict rates of microbial growth. We implemented this method by linking two software tools, PHREEQC and COBRA Toolbox. We applied this method to acetotrophic methanogenesis by Methanosarcina barkeri, and compared the simulation results to previous laboratory observations. The new method constrains acetate uptake by accounting for the kinetics and thermodynamics of methanogenesis, and predicted well the observations of previous experiments. In comparison, traditional methods of dynamic-FBA constrain acetate uptake on the basis of enzyme kinetics, and failed to reproduce the experimental results. These results show that microbial rate laws may provide a better constraint than enzyme kinetics for applying FBA to biogeochemical reaction modeling.

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

Science.gov (United States)

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.

Sensitivity Analysis and Parameter Estimation for a Reactive Transport Model of Uranium Bioremediation

Science.gov (United States)

Meyer, P. D.; Yabusaki, S.; Curtis, G. P.; Ye, M.; Fang, Y.

2011-12-01

A three-dimensional, variably-saturated flow and multicomponent biogeochemical reactive transport model of uranium bioremediation was used to generate synthetic data . The 3-D model was based on a field experiment at the U.S. Dept. of Energy Rifle Integrated Field Research Challenge site that used acetate biostimulation of indigenous metal reducing bacteria to catalyze the conversion of aqueous uranium in the +6 oxidation state to immobile solid-associated uranium in the +4 oxidation state. A key assumption in past modeling studies at this site was that a comprehensive reaction network could be developed largely through one-dimensional modeling. Sensitivity analyses and parameter estimation were completed for a 1-D reactive transport model abstracted from the 3-D model to test this assumption, to identify parameters with the greatest potential to contribute to model predictive uncertainty, and to evaluate model structure and data limitations. Results showed that sensitivities of key biogeochemical concentrations varied in space and time, that model nonlinearities and/or parameter interactions have a significant impact on calculated sensitivities, and that the complexity of the model's representation of processes affecting Fe(II) in the system may make it difficult to correctly attribute observed Fe(II) behavior to modeled processes. Non-uniformity of the 3-D simulated groundwater flux and averaging of the 3-D synthetic data for use as calibration targets in the 1-D modeling resulted in systematic errors in the 1-D model parameter estimates and outputs. This occurred despite using the same reaction network for 1-D modeling as used in the data-generating 3-D model. Predictive uncertainty of the 1-D model appeared to be significantly underestimated by linear parameter uncertainty estimates.

IIASA's climate-vegetation-biogeochemical cycle module as a part of an integrated model for climate change

International Nuclear Information System (INIS)

Ganopolski, A.V.; Jonas, M.; Krabec, J.; Olendrzynski, K.; Petoukhov, V.K.; Venevsky, S.V.

1994-01-01

The main objective of this study is the development of a hierarchy of coupled climate biosphere models with a full description of the global biogeochemical cycles. These models are planned for use as the core of a set of integrated models of climate change and they will incorporate the main elements of the Earth system (atmosphere, hydrosphere, pedosphere and biosphere) linked with each other (and eventually with the antroposphere) through the fluxes of heat, momentum, water and through the global biogeochemical cycles of carbon and nitrogen. This set of integrated models can be considered to fill the gap between highly simplified integrated models of climate change and very sophisticated and computationally expensive coupled models, developed on the basis of general circulation models (GCMs). It is anticipated that this range of integrated models will be an effective tool for investigating the broad spectrum of problems connected with the coexistence of human society and biosphere

Multiscale Investigation on Biofilm Distribution and Its Impact on Macroscopic Biogeochemical Reaction Rates: BIOFILM DISTRIBUTION AND RATE SCALING

Energy Technology Data Exchange (ETDEWEB)

Yan, Zhifeng [Institute of Surface-Earth System Science, Tianjin University, Tianjin China; Pacific Northwest National Laboratory, Richland WA USA; Liu, Chongxuan [Pacific Northwest National Laboratory, Richland WA USA; School of Environmental Science and Engineering, Southern University of Science and Technology, Shenzhen China; Liu, Yuanyuan [Pacific Northwest National Laboratory, Richland WA USA; School of Earth Science and Engineering, Nanjing University, Nanjing China; Bailey, Vanessa L. [Pacific Northwest National Laboratory, Richland WA USA

2017-11-01

Biofilms are critical locations for biogeochemical reactions in the subsurface environment. The occurrence and distribution of biofilms at microscale as well as their impacts on macroscopic biogeochemical reaction rates are still poorly understood. This paper investigated the formation and distributions of biofilms in heterogeneous sediments using multiscale models, and evaluated the effects of biofilm heterogeneity on local and macroscopic biogeochemical reaction rates. Sediment pore structures derived from X-ray computed tomography were used to simulate the microscale flow dynamics and biofilm distribution in the sediment column. The response of biofilm formation and distribution to the variations in hydraulic and chemical properties was first examined. One representative biofilm distribution was then utilized to evaluate its effects on macroscopic reaction rates using nitrate reduction as an example. The results revealed that microorganisms primarily grew on the surfaces of grains and aggregates near preferential flow paths where both electron donor and acceptor were readily accessible, leading to the heterogeneous distribution of biofilms in the sediments. The heterogeneous biofilm distribution decreased the macroscopic rate of biogeochemical reactions as compared with those in homogeneous cases. Operationally considering the heterogeneous biofilm distribution in macroscopic reactive transport models such as using dual porosity domain concept can significantly improve the prediction of biogeochemical reaction rates. Overall, this study provided important insights into the biofilm formation and distribution in soils and sediments as well as their impacts on the macroscopic manifestation of reaction rates.

Biogeochemical processes in a clay formation in situ experiment: Part F - Reactive transport modelling

Energy Technology Data Exchange (ETDEWEB)

Tournassat, Christophe, E-mail: c.tournassat@brgm.fr [BRGM, French Geological Survey, Orleans (France); Alt-Epping, Peter [Rock-Water Interaction Group, Institute of Geological Sciences, University of Bern (Switzerland); Gaucher, Eric C. [BRGM, French Geological Survey, Orleans (France); Gimmi, Thomas [Rock-Water Interaction Group, Institute of Geological Sciences, University of Bern (Switzerland)] [Laboratory for Waste Management, Paul Scherrer Institut, Villigen (Switzerland); Leupin, Olivier X. [NAGRA, CH-5430 Wettingen (Switzerland); Wersin, Paul [Gruner Ltd., CH-4020 Basel (Switzerland)

2011-06-15

Highlights: > Reactive transport modelling was used to simulate simultaneously solute transport, thermodynamic reactions, ion exchange and biodegradation during an in-situ experiment in a clay-rock formation. > Opalinus clay formation has a high buffering capacity in terms of chemical perturbations caused by bacterial activity. > Buffering capacity is mainly attributed to the carbonate system and to the reactivity of clay surfaces (cation exchange, pH buffering). - Abstract: Reactive transport modelling was used to simulate solute transport, thermodynamic reactions, ion exchange and biodegradation in the Porewater Chemistry (PC) experiment at the Mont Terri Rock Laboratory. Simulations show that the most important chemical processes controlling the fluid composition within the borehole and the surrounding formation during the experiment are ion exchange, biodegradation and dissolution/precipitation reactions involving pyrite and carbonate minerals. In contrast, thermodynamic mineral dissolution/precipitation reactions involving alumo-silicate minerals have little impact on the fluid composition on the time-scale of the experiment. With the accurate description of the initial chemical condition in the formation in combination with kinetic formulations describing the different stages of bacterial activities, it has been possible to reproduce the evolution of important system parameters, such as the pH, redox potential, total organic C, dissolved inorganic C and SO{sub 4} concentration. Leaching of glycerol from the pH-electrode may be the primary source of organic material that initiated bacterial growth, which caused the chemical perturbation in the borehole. Results from these simulations are consistent with data from the over-coring and demonstrate that the Opalinus Clay has a high buffering capacity in terms of chemical perturbations caused by bacterial activity. This buffering capacity can be attributed to the carbonate system as well as to the reactivity of

Biogeochemical Controls on Technetium Mobility in Biogeochemical Controls on Technetium Mobility in FRC Sediments

International Nuclear Information System (INIS)

Lloyd, J.R.; McBeth, J.M.; Livens, F.R.; Bryan, N.D.; Ellis, B.; Sharma, H.; Burke, I.T.; Morris, K.

2004-01-01

Technetium-99 is a priority pollutant at numerous DOE sites, due to its long half-life (2.1 x 10 5 years), high mobility as Tc(VII) in oxic waters, and bioavailability as a sulfate analog. 99 Tc is far less mobile under anaerobic conditions, forming insoluble Tc(IV) precipitates. As anaerobic microorganisms can reduce soluble Tc(VII) to insoluble Tc(IV), microbial metabolism may have the potential to treat sediments and waters contaminated with Tc. Baseline studies of fundamental mechanisms of Tc(VII) bioreduction and precipitation (reviewed by Lloyd et al, 2002) have generally used pure cultures of metal-reducing bacteria, in order to develop conceptual models for the biogeochemical cycling of Tc. There is, however, comparatively little known about interactions of metal-reducing bacteria with environmentally relevant trace concentrations of Tc, against a more complex biogeochemical background provided by mixed microbial communities in the subsurface. The objective of this new NABIR project is to probe the site specific biogeochemical conditions that control the mobility of Tc at the FRC (Oak Ridge, TN). This information is required for the rational design of in situ bioremediation strategies for technetium-contaminated subsurface environments. We will use a combination of geochemical, mineralogical, microbiological and spectroscopic techniques to determine the solubility and phase associations of Tc in FRC sediments, and characterize the underpinning biogeochemical controls. A key strength of this project is that many of the techniques we are using have already been optimized by our research team, who are also studying the biogeochemical controls on Tc mobility in marine and freshwater sediments in the UK in a NERC funded companion study.

Traceable components of terrestrial carbon storage capacity in biogeochemical models.

Science.gov (United States)

Xia, Jianyang; Luo, Yiqi; Wang, Ying-Ping; Hararuk, Oleksandra

2013-07-01

Biogeochemical models have been developed to account for more and more processes, making their complex structures difficult to be understood and evaluated. Here, we introduce a framework to decompose a complex land model into traceable components based on mutually independent properties of modeled biogeochemical processes. The framework traces modeled ecosystem carbon storage capacity (Xss ) to (i) a product of net primary productivity (NPP) and ecosystem residence time (τE ). The latter τE can be further traced to (ii) baseline carbon residence times (τ'E ), which are usually preset in a model according to vegetation characteristics and soil types, (iii) environmental scalars (ξ), including temperature and water scalars, and (iv) environmental forcings. We applied the framework to the Australian Community Atmosphere Biosphere Land Exchange (CABLE) model to help understand differences in modeled carbon processes among biomes and as influenced by nitrogen processes. With the climate forcings of 1990, modeled evergreen broadleaf forest had the highest NPP among the nine biomes and moderate residence times, leading to a relatively high carbon storage capacity (31.5 kg cm(-2) ). Deciduous needle leaf forest had the longest residence time (163.3 years) and low NPP, leading to moderate carbon storage (18.3 kg cm(-2) ). The longest τE in deciduous needle leaf forest was ascribed to its longest τ'E (43.6 years) and small ξ (0.14 on litter/soil carbon decay rates). Incorporation of nitrogen processes into the CABLE model decreased Xss in all biomes via reduced NPP (e.g., -12.1% in shrub land) or decreased τE or both. The decreases in τE resulted from nitrogen-induced changes in τ'E (e.g., -26.7% in C3 grassland) through carbon allocation among plant pools and transfers from plant to litter and soil pools. Our framework can be used to facilitate data model comparisons and model intercomparisons via tracking a few traceable components for all terrestrial carbon

Modeling the nitrogen fluxes in the Black Sea using a 3D coupledhydrodynamical-biogeochemical model: transport versus biogeochemicalprocesses, exchanges across the shelf break and comparison of the shelf anddeep sea ecodynamics

Directory of Open Access Journals (Sweden)

M. Grégoire

2004-01-01

Full Text Available A 6-compartment biogeochemical model of nitrogen cycling and plankton productivity has been coupled with a 3D general circulation model in an enclosed environment (the Black Sea so as to quantify and compare, on a seasonal and annual scale, the typical internal biogeochemical functioning of the shelf and of the deep sea as well as to estimate the nitrogen and water exchanges at the shelf break. Model results indicate that the annual nitrogen net export to the deep sea roughly corresponds to the annual load of nitrogen discharged by the rivers on the shelf. The model estimated vertically integrated gross annual primary production is 130gCm-2yr-1 for the whole basin, 220gCm-2yr-1 for the shelf and 40gCm-2yr-1 for the central basin. In agreement with sediment trap observations, model results indicate a rapid and efficient recycling of particulate organic matter in the sub-oxic portion of the water column (60-80m of the open sea. More than 95% of the PON produced in the euphotic layer is recycled in the upper 100m of the water column, 87% in the upper 80 m and 67% in the euphotic layer. The model estimates the annual export of POC towards the anoxic layer to 4 1010molyr-1. This POC is definitely lost for the system and represents 2% of the annual primary production of the open sea.

Assessing the utility of frequency dependent nudging for reducing biases in biogeochemical models

Science.gov (United States)

Lagman, Karl B.; Fennel, Katja; Thompson, Keith R.; Bianucci, Laura

2014-09-01

Bias errors, resulting from inaccurate boundary and forcing conditions, incorrect model parameterization, etc. are a common problem in environmental models including biogeochemical ocean models. While it is important to correct bias errors wherever possible, it is unlikely that any environmental model will ever be entirely free of such errors. Hence, methods for bias reduction are necessary. A widely used technique for online bias reduction is nudging, where simulated fields are continuously forced toward observations or a climatology. Nudging is robust and easy to implement, but suppresses high-frequency variability and introduces artificial phase shifts. As a solution to this problem Thompson et al. (2006) introduced frequency dependent nudging where nudging occurs only in prescribed frequency bands, typically centered on the mean and the annual cycle. They showed this method to be effective for eddy resolving ocean circulation models. Here we add a stability term to the previous form of frequency dependent nudging which makes the method more robust for non-linear biological models. Then we assess the utility of frequency dependent nudging for biological models by first applying the method to a simple predator-prey model and then to a 1D ocean biogeochemical model. In both cases we only nudge in two frequency bands centered on the mean and the annual cycle, and then assess how well the variability in higher frequency bands is recovered. We evaluate the effectiveness of frequency dependent nudging in comparison to conventional nudging and find significant improvements with the former.

Functional Enzyme-Based Approach for Linking Microbial Community Functions with Biogeochemical Process Kinetics

Energy Technology Data Exchange (ETDEWEB)

Li, Minjing [School; Qian, Wei-jun [Pacific Northwest National Laboratory, Richland, Washington 99354, United States; Gao, Yuqian [Pacific Northwest National Laboratory, Richland, Washington 99354, United States; Shi, Liang [School; Liu, Chongxuan [Pacific Northwest National Laboratory, Richland, Washington 99354, United States; School

2017-09-28

The kinetics of biogeochemical processes in natural and engineered environmental systems are typically described using Monod-type or modified Monod-type models. These models rely on biomass as surrogates for functional enzymes in microbial community that catalyze biogeochemical reactions. A major challenge to apply such models is the difficulty to quantitatively measure functional biomass for constraining and validating the models. On the other hand, omics-based approaches have been increasingly used to characterize microbial community structure, functions, and metabolites. Here we proposed an enzyme-based model that can incorporate omics-data to link microbial community functions with biogeochemical process kinetics. The model treats enzymes as time-variable catalysts for biogeochemical reactions and applies biogeochemical reaction network to incorporate intermediate metabolites. The sequences of genes and proteins from metagenomes, as well as those from the UniProt database, were used for targeted enzyme quantification and to provide insights into the dynamic linkage among functional genes, enzymes, and metabolites that are necessary to be incorporated in the model. The application of the model was demonstrated using denitrification as an example by comparing model-simulated with measured functional enzymes, genes, denitrification substrates and intermediates

How to `Elk-test' biogeochemical models in a data rich world? (Invited)

Science.gov (United States)

Reichstein, M.; Ciais, P.; Seneviratne, S. I.; Carvalhais, N.; Dalmonech, D.; Jung, M.; Luo, Y.; Mahecha, M. D.; Moffat, A. M.; Tomelleri, E.; Zaehle, S.

2010-12-01

Process-oriented biogeochemical models are a primary tool that has been used to project future states of climate and ecosystems in the earth system in response to anthropogenic and other forcing, and receive tremendous attention also in the context us the planned assessment report AR5 by the IPCC. However, model intercomparison and data-model comparison studies indicate large uncertainties regarding predictions of global interactions between atmosphere and biosphere. Rigorous scientific testing of these models is essential but very challenging, largely because neither it is technically and ethically possible to perform global earth-scale experiments, nor do we have replicate Earths for hypothesis testing. Hence, model evaluations have to rely on monitoring data such as ecological observation networks, global remote sensing or short-term and small-scale experiments. Here, we critically examine strategies of how model evaluations have been performed with a particular emphasis on terrestrial ecosystems. Often weak ‘validations’ are being presented which do not take advantage of all the relevant information in the observed data, but also apparent falsifications are made, that are hampered by a confusion of system processes with system behavior. We propose that a stronger integration of recent advances in pattern-oriented and system-oriented methodologies will lead to more satisfying earth system model evaluation and development, and show a few enlightening examples from terrestrial biogeochemical modeling and other disciplines. Moreover it is crucial to take advantage of the multidimensional nature of arising earth observation data sets which should be matched by models simultaneously, instead of relying on univariate simple comparisons. A new critical model evaluation is needed to improve future IPCC assessments in order to reduce uncertainties by distinguishing plausible simulation trajectories from fairy tales.

Deriving forest fire ignition risk with biogeochemical process modelling.

Science.gov (United States)

Eastaugh, C S; Hasenauer, H

2014-05-01

Climate impacts the growth of trees and also affects disturbance regimes such as wildfire frequency. The European Alps have warmed considerably over the past half-century, but incomplete records make it difficult to definitively link alpine wildfire to climate change. Complicating this is the influence of forest composition and fuel loading on fire ignition risk, which is not considered by purely meteorological risk indices. Biogeochemical forest growth models track several variables that may be used as proxies for fire ignition risk. This study assesses the usefulness of the ecophysiological model BIOME-BGC's 'soil water' and 'labile litter carbon' variables in predicting fire ignition. A brief application case examines historic fire occurrence trends over pre-defined regions of Austria from 1960 to 2008. Results show that summer fire ignition risk is largely a function of low soil moisture, while winter fire ignitions are linked to the mass of volatile litter and atmospheric dryness.

Kinetic modeling of microbially-driven redox chemistry of radionuclides in subsurface environments: Coupling transport, microbial metabolism and geochemistry

International Nuclear Information System (INIS)

Wang, Yifeng; Papenguth, Hans W.

2000-01-01

Microbial degradation of organic matter is a driving force in many subsurface geochemical systems, and therefore may have significant impacts on the fate of radionuclides released into subsurface environments. In this paper, the authors present a general reaction-transport model for microbial metabolism, redox chemistry, and radionuclide migration in subsurface systems. The model explicitly accounts for biomass accumulation and the coupling of radionuclide redox reactions with major biogeochemical processes. Based on the consideration that the biomass accumulation in subsurface environments is likely to achieve a quasi-steady state, they have accordingly modified the traditional microbial growth kinetic equation. They justified the use of the biogeochemical models without the explicit representation of biomass accumulation, if the interest of modeling is in the net impact of microbial reactions on geochemical processes. They then applied their model to a scenario in which an oxic water flow containing both uranium and completing organic ligands is recharged into an oxic aquifer in a carbonate formation. The model simulation shows that uranium can be reduced and therefore immobilized in the anoxic zone created by microbial degradation

Kinetic modeling of microbially-driven redox chemistry of radionuclides in subsurface environments: Coupling transport, microbial metabolism and geochemistry

Energy Technology Data Exchange (ETDEWEB)

WANG,YIFENG; PAPENGUTH,HANS W.

2000-05-04

Microbial degradation of organic matter is a driving force in many subsurface geochemical systems, and therefore may have significant impacts on the fate of radionuclides released into subsurface environments. In this paper, the authors present a general reaction-transport model for microbial metabolism, redox chemistry, and radionuclide migration in subsurface systems. The model explicitly accounts for biomass accumulation and the coupling of radionuclide redox reactions with major biogeochemical processes. Based on the consideration that the biomass accumulation in subsurface environments is likely to achieve a quasi-steady state, they have accordingly modified the traditional microbial growth kinetic equation. They justified the use of the biogeochemical models without the explicit representation of biomass accumulation, if the interest of modeling is in the net impact of microbial reactions on geochemical processes. They then applied their model to a scenario in which an oxic water flow containing both uranium and completing organic ligands is recharged into an oxic aquifer in a carbonate formation. The model simulation shows that uranium can be reduced and therefore immobilized in the anoxic zone created by microbial degradation.

Stream biogeochemical resilience in the age of Anthropocene

Science.gov (United States)

Dong, H.; Creed, I. F.

2017-12-01

Recent evidence indicates that biogeochemical cycles are being pushed beyond the tolerance limits of the earth system in the age of the Anthropocene placing terrestrial and aquatic ecosystems at risk. Here, we explored the question: Is there empirical evidence of global atmospheric changes driving losses in stream biogeochemical resilience towards a new normal? Stream biogeochemical resilience is the process of returning to equilibrium conditions after a disturbance and can be measured using three metrics: reactivity (the highest initial response after a disturbance), return rate (the rate of return to equilibrium condition after reactive changes), and variance of the stationary distribution (the signal to noise ratio). Multivariate autoregressive models were used to derive the three metrics for streams along a disturbance gradient - from natural systems where global drivers would dominate, to relatively managed or modified systems where global and local drivers would interact. We observed a loss of biogeochemical resilience in all streams. The key biogeochemical constituent(s) that may be driving loss of biogeochemical resilience were identified from the time series of the stream biogeochemical constituents. Non-stationary trends (detected by Mann-Kendall analysis) and stationary cycles (revealed through Morlet wavelet analysis) were removed, and the standard deviation (SD) of the remaining residuals were analyzed to determine if there was an increase in SD over time that would indicate a pending shift towards a new normal. We observed that nitrate-N and total phosphorus showed behaviours indicative of a pending shift in natural and managed forest systems, but not in agricultural systems. This study provides empirical support that stream ecosystems are showing signs of exceeding planetary boundary tolerance levels and shifting towards a "new normal" in response to global changes, which can be exacerbated by local management activities. Future work will consider

Biogeochemical speciation of Fe in ocean water

NARCIS (Netherlands)

Hiemstra, T.; Riemsdijk, van W.H.

2006-01-01

The biogeochemical speciation of Fe in seawater has been evaluated using the consistent Non-Ideal Competitive Adsorption model (NICA¿Donnan model). Two types of data sets were used, i.e. Fe-hydroxide solubility data and competitive ligand equilibration/cathodic stripping voltammetry (CLE/CSV) Fe

Predictive Understanding of Mountainous Watershed Hydro-Biogeochemical Function and Response to Perturbations

Science.gov (United States)

Hubbard, S. S.; Williams, K. H.; Agarwal, D.; Banfield, J. F.; Beller, H. R.; Bouskill, N.; Brodie, E.; Maxwell, R. M.; Nico, P. S.; Steefel, C. I.; Steltzer, H.; Tokunaga, T. K.; Wainwright, H. M.; Dwivedi, D.; Newcomer, M. E.

2017-12-01

Recognizing the societal importance, vulnerability and complexity of mountainous watersheds, the `Watershed Function' project is developing a predictive understanding of how mountainous watersheds retain and release downgradient water, nutrients, carbon, and metals. In particular, the project is exploring how early snowmelt, drought, floods and other disturbances will influence mountainous watershed dynamics at seasonal to decadal timescales. Located in the 300km2 East River headwater catchment of the Upper Colorado River Basin, the project is guided by several constructs. First, the project considers the integrated role of surface and subsurface flow and biogeochemical reactions - from bedrock to the top of the vegetative canopy, from terrestrial through aquatic compartments, and from summit to receiving waters. The project takes a system-of-systems perspective, focused on developing new methods to quantify the cumulative watershed hydrobiogeochemical response to perturbations based on information from select subsystems within the watershed, each having distinct vegetation-subsurface biogeochemical-hydrological characteristics. A `scale-adaptive' modeling capability, in development using adaptive mesh refinement methods, serves as the organizing framework for the SFA. The scale-adaptive approach is intended to permit simulation of system-within-systems behavior - and aggregation of that behavior - from genome through watershed scales. This presentation will describe several early project discoveries and advances made using experimental, observational and numerical approaches. Among others, examples may include:quantiying how seasonal hydrological perturbations drive biogeochemical responses across critical zone compartments, with a focus on N and C transformations; metagenomic documentation of the spatial variability in floodplain meander microbial ecology; 3D reactive transport simulations of couped hydrological-biogeochemical behavior in the hyporheic zone; and

Carbon sequestration by patch fertilization: A comprehensive assessment using coupled physical-ecological-biogeochemical models

Energy Technology Data Exchange (ETDEWEB)

Sarmiento, Jorge L. [Princeton Univ., NJ (United States); Gnanadesikan, Anand [Princeton Univ., NJ (United States); Gruber, Nicolas [Univ. of California, Los Angeles, CA (United States); Jin, Xin [Univ. of California, Los Angeles, CA (United States); Armstrong, Robert [State Univ. of New York (SUNY), Plattsburgh, NY (United States)

2007-06-21

This final report summarizes research undertaken collaboratively between Princeton University, the NOAA Geophysical Fluid Dynamics Laboratory on the Princeton University campus, the State University of New York at Stony Brook, and the University of California, Los Angeles between September 1, 2000, and November 30, 2006, to do fundamental research on ocean iron fertilization as a means to enhance the net oceanic uptake of CO2 from the atmosphere. The approach we proposed was to develop and apply a suite of coupled physical-ecological-biogeochemical models in order to (i) determine to what extent enhanced carbon fixation from iron fertilization will lead to an increase in the oceanic uptake of atmospheric CO2 and how long this carbon will remain sequestered (efficiency), and (ii) examine the changes in ocean ecology and natural biogeochemical cycles resulting from iron fertilization (consequences). The award was funded in two separate three-year installments: September 1, 2000 to November 30, 2003, for a project entitled “Ocean carbon sequestration by fertilization: An integrated biogeochemical assessment.” A final report was submitted for this at the end of 2003 and is included here as Appendix 1; and, December 1, 2003 to November 30, 2006, for a follow-on project under the same grant number entitled “Carbon sequestration by patch fertilization: A comprehensive assessment using coupled physical-ecological-biogeochemical models.” This report focuses primarily on the progress we made during the second period of funding subsequent to the work reported on in Appendix 1. When we began this project, we were thinking almost exclusively in terms of long-term fertilization over large regions of the ocean such as the Southern Ocean, with much of our focus being on how ocean circulation and biogeochemical cycling would interact to control the response to a given fertilization scenario. Our research on these types of scenarios, which was carried out largely during the

Enhanced biogeochemical cycling and subsequent reduction of hydraulic conductivity associated with soil-layer interfaces in the vadose zone

Science.gov (United States)

Hansen, David J.; McGuire, Jennifer T.; Mohanty, Binayak P.

2013-01-01

Biogeochemical dynamics in the vadose zone are poorly understood due to the transient nature of chemical and hydrologic conditions, but are nonetheless critical to understanding chemical fate and transport. This study explored the effects of a soil layer on linked geochemical, hydrological, and microbiological processes. Three laboratory soil columns were constructed: a homogenized medium-grained sand, a homogenized organic-rich loam, and a sand-over-loam layered column. Upward and downward infiltration of water was evaluated during experiments to simulate rising water table and rainfall events respectively. In-situ collocated probes measured soil water content, matric potential, and Eh while water samples collected from the same locations were analyzed for Br−, Cl−, NO3−, SO42−, NH4+, Fe2+, and total sulfide. Compared to homogenous columns, the presence of a soil layer altered the biogeochemistry and water flow of the system considerably. Enhanced biogeochemical cycling was observed in the layered column over the texturally homogeneous soil columns. Enumerations of iron and sulfate reducing bacteria showed 1-2 orders of magnitude greater community numbers in the layered column. Mineral and soil aggregate composites were most abundant near the soil-layer interface; the presence of which, likely contributed to an observed order-of-magnitude decrease in hydraulic conductivity. These findings show that quantifying coupled hydrologic-biogeochemical processes occurring at small-scale soil interfaces is critical to accurately describing and predicting chemical changes at the larger system scale. Findings also provide justification for considering soil layering in contaminant fate and transport models because of its potential to increase biodegradation and/or slow the rate of transport of contaminants. PMID:22031578

Connections between physical, optical and biogeochemical processes in the Pacific Ocean

Science.gov (United States)

Xiu, Peng; Chai, Fei

2014-03-01

A new biogeochemical model has been developed and coupled to a three-dimensional physical model in the Pacific Ocean. With the explicitly represented dissolved organic pools, this new model is able to link key biogeochemical processes with optical processes. Model validation against satellite and in situ data indicates the model is robust in reproducing general biogeochemical and optical features. Colored dissolved organic matter (CDOM) has been suggested to play an important role in regulating underwater light field. With the coupled model, physical and biological regulations of CDOM in the euphotic zone are analyzed. Model results indicate seasonal variability of CDOM is mostly determined by biological processes, while the importance of physical regulation manifests in the annual mean terms. Without CDOM attenuating light, modeled depth-integrated primary production is about 10% higher than the control run when averaged over the entire basin, while this discrepancy is highly variable in space with magnitudes reaching higher than 100% in some locations. With CDOM dynamics integrated in physical-biological interactions, a new mechanism by which physical processes affect biological processes is suggested, namely, physical transport of CDOM changes water optical properties, which can further modify underwater light field and subsequently affect the distribution of phytoplankton chlorophyll. This mechanism tends to occur in the entire Pacific basin but with strong spatial variability, implying the importance of including optical processes in the coupled physical-biogeochemical model. If ammonium uptake is sufficient to permit utilization of DOM, that is, UB∗⩾-U{U}/{U}-{(1-r_b)}/{RB}, then bacteria uptake of DOM has the form of FB=(1-r_b){U}/{RB}, bacteria respiration, SB=r_b×U, remineralization by bacteria, EB=UC{UN}/{UC}-{(1-r_b)}/{RB}. If EB > 0, then UB = 0; otherwise, UB = -EB. If there is insufficient ammonium, that is, UB∗CO2 is calculated using the

Assessment of the GHG Reduction Potential from Energy Crops Using a Combined LCA and Biogeochemical Process Models: A Review

Directory of Open Access Journals (Sweden)

Dong Jiang

2014-01-01

Full Text Available The main purpose for developing biofuel is to reduce GHG (greenhouse gas emissions, but the comprehensive environmental impact of such fuels is not clear. Life cycle analysis (LCA, as a complete comprehensive analysis method, has been widely used in bioenergy assessment studies. Great efforts have been directed toward establishing an efficient method for comprehensively estimating the greenhouse gas (GHG emission reduction potential from the large-scale cultivation of energy plants by combining LCA with ecosystem/biogeochemical process models. LCA presents a general framework for evaluating the energy consumption and GHG emission from energy crop planting, yield acquisition, production, product use, and postprocessing. Meanwhile, ecosystem/biogeochemical process models are adopted to simulate the fluxes and storage of energy, water, carbon, and nitrogen in the soil-plant (energy crops soil continuum. Although clear progress has been made in recent years, some problems still exist in current studies and should be addressed. This paper reviews the state-of-the-art method for estimating GHG emission reduction through developing energy crops and introduces in detail a new approach for assessing GHG emission reduction by combining LCA with biogeochemical process models. The main achievements of this study along with the problems in current studies are described and discussed.

Study of the seasonal cycle of the biogeochemical processes in the Ligurian Sea using a 1D interdisciplinary model

NARCIS (Netherlands)

Raick, C.; Delhez, E.J.M.; Soetaert, K.E.R.; Grégoire, M.

2005-01-01

A one-dimensional coupled physical–biogeochemical model has been built to study the pelagic food web of the Ligurian Sea (NW Mediterranean Sea). The physical model is the turbulent closure model (version 1D) developed at the GeoHydrodynamics and Environmental Laboratory (GHER) of the University of

Quantifying the Variability of CH4 Emissions from Pan-Arctic Lakes with Lake Biogeochemical and Landscape Evolution Models

Science.gov (United States)

Tan, Z.; Zhuang, Q.

2014-12-01

Recent studies in the arctic and subarctic show that CH4 emissions from pan-arctic lakes are playing much more significant roles in the regional carbon cycling than previously estimated. Permafrost thawing due to pronounced warming at northern high latitudes affects lake morphology, changing its CH4 emissions. Thermokarst can enlarge the extent of artic lakes, exposing stable ancient carbon buried in the permafrost zone for degradation and changing a previously known carbon sink to a large carbon source. In some areas, the thawing of subarctic discontinuous and isolated permafrost can diminish thermokarst lakes. To date, few models have considered these important hydrological and biogeochemical processes to provide adequate estimation of CH4 emissions from these lakes. To fill this gap, we have developed a process-based climate-sensitive lake biogeochemical model and a landscape evolution model, which have been applied to quantify the state and variability of CH4 emissions from this freshwater system. Site-level experiments show the models are capable to capture the spatial and temporal variability of CH4 emissions from lakes across Siberia and Alaska. With the lake biogeochemical model solely, we estimate that the magnitude of CH4 emissions from lakes is 13.2 Tg yr-1 in the north of 60 ºN at present, which is on the same order of CH4 emissions from northern high-latitude wetlands. The maximum increment is 11.8 Tg CH4 yr-1 by the end of the 21st century when the worst warming scenario is assumed. We expect the landscape evolution model will improve the existing estimates.

Biogeochemical Modeling of the Second Rise of Oxygen

Science.gov (United States)

Smith, M. L.; Catling, D.; Claire, M.; Zahnle, K.

2014-03-01

The rise of atmospheric oxygen set the tempo for the evolution of complex life on Earth. Oxygen levels are thought to have increased in two broad steps: one step occurred in the Archean ~ 2.45 Ga (the Great Oxidation Event or GOE), and another step occured in the Neoproterozoic ~750-580 Ma (the Neoprotoerozoic Oxygenation Event or NOE). During the NOE, oxygen levels increased from ~1-10% of the present atmospheric level (PAL) (Holland, 2006), to ~15% PAL in the late Neoproterozoic, to ~100% PAL later in the Phanerozoic. Complex life requires O2, so this transition allowed complex life to evolve. We seek to understand what caused the NOE. To explore causes for the NOE, we build upon the biogeochemical model of Claire et al. (2006), which calculates the redox evolution of the atmosphere, ocean, biosphere, and crust in the Archean through to the early Proterozoic. In this model, the balance between oxygenconsuming and oyxgen-producing fluxes evolves over time such that at ~2.4 Ga, the rapidly acting sources of oxygen outweigh the rapidly-acting sinks. Or, in other words, at ~2.4 Ga, the flux of oxygen from organic carbon burial exceeds the sinks of oxygen from reaction with reduced volcanic and metamoprphic gases. The model is able to drive oxygen levels to 1-10% PAL in the Proterozoic; however, the evolving redox fluxes in the model cannot explain how oxygen levels pushed above 1-10% in the late Proterozoic. The authors suggest that perhaps another buffer, such as sulfur, is needed to describe Proterozoic and Phanerozoic redox evolution. Geologic proxies show that in the Proterozoic, up to 10% of the deep ocean may have been sulfidic. With this ocean chemistry, the global sulfur cycle would have worked differently than it does today. Because the sulfur and oxygen cycles interact, the oxygen concentration could have permanently changed due to an evolving sulfur cycle (in combination with evolving redox fluxes associated with other parts of the oxygen cycle and carbon

Biogeochemical Reactive Transport Model of the Redox Zone Experiment of the sp Hard Rock Laboratory in Sweden

International Nuclear Information System (INIS)

Molinero-Huguet, Jorge; Samper-Calvete, F. Javier; Zhang, Guoxiang; Yang, Changbing

2004-01-01

Underground facilities are being operated by several countries around the world for performing research and demonstration of the safety of deep radioactive waste repositories. The ''sp'' Hard Rock Laboratory is one such facility launched and operated by the Swedish Nuclear Fuel and Waste Management Company where various in situ experiments have been performed in fractured granites. One such experiment is the redox zone experiment, which aimed at evaluating the effects of the construction of an access tunnel on the hydrochemical conditions of a fracture zone. Dilution of the initially saline groundwater by fresh recharge water is the dominant process controlling the hydrochemical evolution of most chemical species, except for bicarbonate and sulfate, which unexpectedly increase with time. We present a numerical model of water flow, reactive transport, and microbial processes for the redox zone experiment. This model provides a plausible quantitatively based explanation for the unexpected evolution of bicarbonate and sulfate, reproduces the breakthrough curves of other reactive species, and is consistent with previous hydrogeological and solute transport models

Ecotoxicological, ecophysiological, and biogeochemical fundamentals of risk assessment

International Nuclear Information System (INIS)

Bashkin, V.N.; Kozlov, M.Ya.; Evstafjeva, E.V.

1993-01-01

Risk assessment (RA) influenced by different factors in radionuclide polluted regions is carried out by determining the biogeochemical structure of a region. Consequently, ecological-biogeochemical regionalization, ecotoxicological and ecophysiological monitoring of human population health are the important approach to RA. These criteria should conjugate with LCA of various industrial and agricultural products. Given fundamentals and approaches are needed for areas where traditional pollutants (heavy metals, pesticides, fertilizers, POPs etc) are enforced sharply by radioactive pollution. For RA of these complex pollutants, the methods of human adaptability to a polluted environment have been carried out. These techniques include biogeochemical, ecotoxicological, and ecophysiological analyses of risk factors as well as quantitative analysis of uncertainties using expert-modeling systems. Furthermore, the modern statistical methods are used for quantitative assessment of human adaptability to radioactive and nonradioactive pollutants. The results obtained in Chernobyl regions show the acceptability of these methods for risk assessment

Coupled Modeling of Rhizosphere and Reactive Transport Processes

Science.gov (United States)

Roque-Malo, S.; Kumar, P.

2017-12-01

The rhizosphere, as a bio-diverse plant root-soil interface, hosts many hydrologic and biochemical processes, including nutrient cycling, hydraulic redistribution, and soil carbon dynamics among others. The biogeochemical function of root networks, including the facilitation of nutrient cycling through absorption and rhizodeposition, interaction with micro-organisms and fungi, contribution to biomass, etc., plays an important role in myriad Critical Zone processes. Despite this knowledge, the role of the rhizosphere on watershed-scale ecohydrologic functions in the Critical Zone has not been fully characterized, and specifically, the extensive capabilities of reactive transport models (RTMs) have not been applied to these hydrobiogeochemical dynamics. This study uniquely links rhizospheric processes with reactive transport modeling to couple soil biogeochemistry, biological processes, hydrologic flow, hydraulic redistribution, and vegetation dynamics. Key factors in the novel modeling approach are: (i) bi-directional effects of root-soil interaction, such as simultaneous root exudation and nutrient absorption; (ii) multi-state biomass fractions in soil (i.e. living, dormant, and dead biological and root materials); (iii) expression of three-dimensional fluxes to represent both vertical and lateral interconnected flows and processes; and (iv) the potential to include the influence of non-stationary external forcing and climatic factors. We anticipate that the resulting model will demonstrate the extensive effects of plant root dynamics on ecohydrologic functions at the watershed scale and will ultimately contribute to a better characterization of efflux from both agricultural and natural systems.

Modeling the fate of nitrogen on the catchment scale using a spatially explicit hydro-biogeochemical simulation system

Science.gov (United States)

Klatt, S.; Butterbach-Bahl, K.; Kiese, R.; Haas, E.; Kraus, D.; Molina-Herrera, S. W.; Kraft, P.

2015-12-01

The continuous growth of the human population demands an equally growing supply for fresh water and food. As a result, available land for efficient agriculture is constantly diminishing which forces farmers to cultivate inferior croplands and intensify agricultural practices, e.g., increase the use of synthetic fertilizers. This intensification of marginal areas in particular will cause a dangerous rise in nitrate discharge into open waters or even drinking water resources. In order to reduce the amount of nitrate lost by surface runoff or lateral subsurface transport, bufferstrips have proved to be a valuable means. Current laws, however, promote rather static designs (i.e., width and usage) even though a multitude of factors, e.g., soil type, slope, vegetation and the nearby agricultural management, determines its effectiveness. We propose a spatially explicit modeling approach enabling to assess the effects of those factors on nitrate discharge from arable lands using the fully distributed hydrology model CMF coupled to the complex biogeochemical model LandscapeDNDC. Such a modeling scheme allows to observe the displacement of dissolved nutrients in both vertical and horizontal directions and serves to estimate both their uptake by the vegetated bufferstrip and loss to the environment. First results indicate a significant reduction of nitrate loss in the presence of a bufferstrip (2.5 m). We show effects induced by various buffer strip widths and plant cover on the nitrate retention.

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

Science.gov (United States)

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.

Wetland biogeochemical processes and simulation modeling

Science.gov (United States)

Bai, Junhong; Huang, Laibin; Gao, Haifeng; Jia, Jia; Wang, Xin

2018-02-01

As the important landscape with rich biodiversity and high productivity, wetlands can provide numerous ecological services including playing an important role in regulating global biogeochemical cycles, filteringpollutants from terrestrial runoff and atmospheric deposition, protecting and improving water quality, providing living habitats for plants and animals, controlling floodwaters, and retaining surface water flow during dry periods (Reddy and DeLaune, 2008; Qin and Mitsch, 2009; Zhao et al., 2016). However, more than 50% of the world's wetlands had been altered, degraded or lost through a wide range of human activities in the past 150 years, and only a small percentage of the original wetlands remained around the world after over two centuries of intensive development and urbanization (O'connell, 2003; Zhao et al., 2016).

Final Report Coupling in silico microbial models with reactive transport models to predict the fate of contaminants in the subsurface.

Energy Technology Data Exchange (ETDEWEB)

Lovley, Derek R.

2012-10-31

This project successfully accomplished its goal of coupling genome-scale metabolic models with hydrological and geochemical models to predict the activity of subsurface microorganisms during uranium bioremediation. Furthermore, it was demonstrated how this modeling approach can be used to develop new strategies to optimize bioremediation. The approach of coupling genome-scale metabolic models with reactive transport modeling is now well enough established that it has been adopted by other DOE investigators studying uranium bioremediation. Furthermore, the basic principles developed during our studies will be applicable to much broader investigations of microbial activities, not only for other types of bioremediation, but microbial metabolism in diversity of environments. This approach has the potential to make an important contribution to predicting the impact of environmental perturbations on the cycling of carbon and other biogeochemical cycles.

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

Science.gov (United States)

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.

Observationally-based Metrics of Ocean Carbon and Biogeochemical Variables are Essential for Evaluating Earth System Model Projections

Science.gov (United States)

Russell, J. L.; Sarmiento, J. L.

2017-12-01

The Southern Ocean is central to the climate's response to increasing levels of atmospheric greenhouse gases as it ventilates a large fraction of the global ocean volume. Global coupled climate models and earth system models, however, vary widely in their simulations of the Southern Ocean and its role in, and response to, the ongoing anthropogenic forcing. Due to its complex water-mass structure and dynamics, Southern Ocean carbon and heat uptake depend on a combination of winds, eddies, mixing, buoyancy fluxes and topography. Understanding how the ocean carries heat and carbon into its interior and how the observed wind changes are affecting this uptake is essential to accurately projecting transient climate sensitivity. Observationally-based metrics are critical for discerning processes and mechanisms, and for validating and comparing climate models. As the community shifts toward Earth system models with explicit carbon simulations, more direct observations of important biogeochemical parameters, like those obtained from the biogeochemically-sensored floats that are part of the Southern Ocean Carbon and Climate Observations and Modeling project, are essential. One goal of future observing systems should be to create observationally-based benchmarks that will lead to reducing uncertainties in climate projections, and especially uncertainties related to oceanic heat and carbon uptake.

Biogeochemical protocols and diagnostics for the CMIP6 Ocean Model Intercomparison Project (OMIP

Directory of Open Access Journals (Sweden)

J. C. Orr

2017-06-01

Full Text Available The Ocean Model Intercomparison Project (OMIP focuses on the physics and biogeochemistry of the ocean component of Earth system models participating in the sixth phase of the Coupled Model Intercomparison Project (CMIP6. OMIP aims to provide standard protocols and diagnostics for ocean models, while offering a forum to promote their common assessment and improvement. It also offers to compare solutions of the same ocean models when forced with reanalysis data (OMIP simulations vs. when integrated within fully coupled Earth system models (CMIP6. Here we detail simulation protocols and diagnostics for OMIP's biogeochemical and inert chemical tracers. These passive-tracer simulations will be coupled to ocean circulation models, initialized with observational data or output from a model spin-up, and forced by repeating the 1948–2009 surface fluxes of heat, fresh water, and momentum. These so-called OMIP-BGC simulations include three inert chemical tracers (CFC-11, CFC-12, SF6 and biogeochemical tracers (e.g., dissolved inorganic carbon, carbon isotopes, alkalinity, nutrients, and oxygen. Modelers will use their preferred prognostic BGC model but should follow common guidelines for gas exchange and carbonate chemistry. Simulations include both natural and total carbon tracers. The required forced simulation (omip1 will be initialized with gridded observational climatologies. An optional forced simulation (omip1-spunup will be initialized instead with BGC fields from a long model spin-up, preferably for 2000 years or more, and forced by repeating the same 62-year meteorological forcing. That optional run will also include abiotic tracers of total dissolved inorganic carbon and radiocarbon, CTabio and 14CTabio, to assess deep-ocean ventilation and distinguish the role of physics vs. biology. These simulations will be forced by observed atmospheric histories of the three inert gases and CO2 as well as carbon isotope ratios of CO2. OMIP-BGC simulation

Biogeochemical Protocols and Diagnostics for the CMIP6 Ocean Model Intercomparison Project (OMIP)

Science.gov (United States)

Orr, James C.; Najjar, Raymond G.; Aumont, Olivier; Bopp, Laurent; Bullister, John L.; Danabasoglu, Gokhan; Doney, Scott C.; Dunne, John P.; Dutay, Jean-Claude; Graven, Heather;

Biogeochemical protocols and diagnostics for the CMIP6 Ocean Model Intercomparison Project (OMIP)

Science.gov (United States)

Orr, James C.; Najjar, Raymond G.; Aumont, Olivier; Bopp, Laurent; Bullister, John L.; Danabasoglu, Gokhan; Doney, Scott C.; Dunne, John P.; Dutay, Jean-Claude; Graven, Heather; Griffies, Stephen M.; John, Jasmin G.; Joos, Fortunat; Levin, Ingeborg; Lindsay, Keith; Matear, Richard J.; McKinley, Galen A.; Mouchet, Anne; Oschlies, Andreas; Romanou, Anastasia; Schlitzer, Reiner; Tagliabue, Alessandro; Tanhua, Toste; Yool, Andrew

2017-06-01

The Ocean Model Intercomparison Project (OMIP) focuses on the physics and biogeochemistry of the ocean component of Earth system models participating in the sixth phase of the Coupled Model Intercomparison Project (CMIP6). OMIP aims to provide standard protocols and diagnostics for ocean models, while offering a forum to promote their common assessment and improvement. It also offers to compare solutions of the same ocean models when forced with reanalysis data (OMIP simulations) vs. when integrated within fully coupled Earth system models (CMIP6). Here we detail simulation protocols and diagnostics for OMIP's biogeochemical and inert chemical tracers. These passive-tracer simulations will be coupled to ocean circulation models, initialized with observational data or output from a model spin-up, and forced by repeating the 1948-2009 surface fluxes of heat, fresh water, and momentum. These so-called OMIP-BGC simulations include three inert chemical tracers (CFC-11, CFC-12, SF6) and biogeochemical tracers (e.g., dissolved inorganic carbon, carbon isotopes, alkalinity, nutrients, and oxygen). Modelers will use their preferred prognostic BGC model but should follow common guidelines for gas exchange and carbonate chemistry. Simulations include both natural and total carbon tracers. The required forced simulation (omip1) will be initialized with gridded observational climatologies. An optional forced simulation (omip1-spunup) will be initialized instead with BGC fields from a long model spin-up, preferably for 2000 years or more, and forced by repeating the same 62-year meteorological forcing. That optional run will also include abiotic tracers of total dissolved inorganic carbon and radiocarbon, CTabio and 14CTabio, to assess deep-ocean ventilation and distinguish the role of physics vs. biology. These simulations will be forced by observed atmospheric histories of the three inert gases and CO2 as well as carbon isotope ratios of CO2. OMIP-BGC simulation protocols are

Application of a SEEK filter to a 1D biogeochemical model of the Ligurian Sea: Twin experiments and real data assimilation

NARCIS (Netherlands)

Raick, C.; Alvera-Azcarate, A.; Barth, A.; Brankart, J.-M.; Soetaert, K.E.R.; Grégoire, M.

2007-01-01

The Singular Evolutive Extended Kalman (SEEK) filter has been implemented to assimilate in-situ data in a 1D coupled physical-ecosystem model of the Ligurian Sea. The biogeochemical model describes the partly decoupled nitrogen and carbon cycles of the pelagic food web. The GHER hydrodynamic 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?

Directory of Open Access Journals (Sweden)

Marion Fraysse

Full Text Available 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.

Biogeochemical provinces in the global ocean based on phytoplankton growth limitation

Science.gov (United States)

Hashioka, T.; Hirata, T.; Aita, M. N.; Chiba, S.

2016-02-01

The biogeochemical province is one of the useful concepts for the comprehensive understanding of regional differences of the marine ecosystem. Various biogeochemical provinces for lower-trophic level ecosystem have been proposed using a similarity-based classification of seasonal variations of chl-a concentration typified by Longhurst 1995 and 2006. Such categorizations well capture the regional differences of seasonality as "total phytoplankton". However, background biogeochemical mechanism to characterize the province boundary is not clear. Namely, the dominant phytoplankton group is different among regions and seasons, and their physiological characteristics are significantly different among groups. Recently some pieces of new biogeochemical information are available. One is an estimation of phytoplankton community structure from satellite observation, and it makes clear the key phytoplankton type in each region. Another is an estimation of limitation factors for phytoplankton growth (e.g., nutrients, temperature, light) in each region from modeling studies. In this study, we propose new biogeochemical provinces as a combination between the dominance of phytoplankton (i.e., diatoms, nano-, pico-phytoplankton or coexistence of two/three types) and their growth limitation factors (particularly we focused on nutrient limitation; N, P, Si or Fe). In this combination, we classified the global ocean into 23 biogeochemical provinces. The result suggests that even if the same type of phytoplankton dominates, the background mechanism could be different among regions. On the contrary, even if the regions geographically separate, the background mechanism could be similar among regions. This is important to understand that region/boundary does respond to environmental change. This biogeochemical province is useful for identification of key areas for future observation.

Upwelling events, coastal offshore exchange, links to biogeochemical processes - Highlights from the Baltic Sea Science Congress

Directory of Open Access Journals (Sweden)

Bogdan Ołdakowski

2008-03-01

Full Text Available The Baltic Sea Science Congress was held at Rostock University, Germany, from 19 to 22 March 2007. In the session entitled"Upwelling events, coastal offshore exchange, links to biogeochemical processes" 20 presentations were given,including 7 talks and 13 posters related to the theme of the session.This paper summarises new findings of the upwelling-related studies reported in the session. It deals with investigationsbased on the use of in situ and remote sensing measurements as well as numerical modelling tools. The biogeochemicalimplications of upwelling are also discussed.Our knowledge of the fine structure and dynamic considerations of upwelling has increased in recent decades with the advent ofhigh-resolution modern measurement techniques and modelling studies. The forcing and the overall structure, duration and intensity ofupwelling events are understood quite well. However, the quantification of related transports and the contribution to the overall mixingof upwelling requires further research. Furthermore, our knowledge of the links between upwelling and biogeochemical processes is stillincomplete. Numerical modelling has advanced to the extent that horizontal resolutions of c. 0.5 nautical miles can now be applied,which allows the complete spectrum of meso-scale features to be described. Even the development of filaments can be describedrealistically in comparison with high-resolution satellite data.But the effect of upwelling at a basin scale and possible changes under changing climatic conditions remain open questions.

CMS: Simulated Physical-Biogeochemical Data, SABGOM Model, Gulf of Mexico, 2005-2010

Data.gov (United States)

National Aeronautics and Space Administration — This dataset contains monthly mean ocean surface physical and biogeochemical data for the Gulf of Mexico simulated by the South Atlantic Bight and Gulf of Mexico...

Characterization of mixing errors in a coupled physical biogeochemical model of the North Atlantic: implications for nonlinear estimation using Gaussian anamorphosis

Directory of Open Access Journals (Sweden)

D. Béal

2010-02-01

Full Text Available In biogeochemical models coupled to ocean circulation models, vertical mixing is an important physical process which governs the nutrient supply and the plankton residence in the euphotic layer. However, vertical mixing is often poorly represented in numerical simulations because of approximate parameterizations of sub-grid scale turbulence, wind forcing errors and other mis-represented processes such as restratification by mesoscale eddies. Getting a sufficient knowledge of the nature and structure of these errors is necessary to implement appropriate data assimilation methods and to evaluate if they can be controlled by a given observation system.

In this paper, Monte Carlo simulations are conducted to study mixing errors induced by approximate wind forcings in a three-dimensional coupled physical-biogeochemical model of the North Atlantic with a 1/4° horizontal resolution. An ensemble forecast involving 200 members is performed during the 1998 spring bloom, by prescribing perturbations of the wind forcing to generate mixing errors. The biogeochemical response is shown to be rather complex because of nonlinearities and threshold effects in the coupled model. The response of the surface phytoplankton depends on the region of interest and is particularly sensitive to the local stratification. In addition, the statistical relationships computed between the various physical and biogeochemical variables reflect the signature of the non-Gaussian behaviour of the system. It is shown that significant information on the ecosystem can be retrieved from observations of chlorophyll concentration or sea surface temperature if a simple nonlinear change of variables (anamorphosis is performed by mapping separately and locally the ensemble percentiles of the distributions of each state variable on the Gaussian percentiles. The results of idealized observational updates (performed with perfect observations and neglecting horizontal correlations

Modeling reactive geochemical transport of concentrated aqueous solutions in variably saturated media

Energy Technology Data Exchange (ETDEWEB)

Zhang, Guoxiang; Zheng, Zuoping; Wan, Jiamin

2004-01-28

Concentrated aqueous solutions (CAS) have unique thermodynamic and physical properties. Chemical components in CAS are incompletely dissociated, especially those containing divalent or polyvalent ions. The problem is further complicated by the interaction between CAS flow processes and the naturally heterogeneous sediments. As the CAS migrates through the porous media, the composition may be altered subject to fluid-rock interactions. To effectively model reactive transport of CAS, we must take into account ion-interaction. A combination of the Pitzer ion-interaction and the ion-association model would be an appropriate way to deal with multiple-component systems if the Pitzer' parameters and thermodynamic data of dissolved components and the related minerals are available. To quantify the complicated coupling of CAS flow and transport, as well as the involved chemical reactions in natural and engineered systems, we have substantially extended an existing reactive biogeochemical transport code, BIO-CORE{sup 2D}{copyright}, by incorporating a comprehensive Pitzer ion-interaction model. In the present paper, the model, and two test cases against measured data were briefly introduced. Finally we present an application to simulate a laboratory column experiment studying the leakage of the high alkaline waste fluid stored in Hanford (a site of the U.S. Department of Energy, located in Washington State, USA). With the Pitzer ion-interaction ionic activity model, our simulation captures measured pH evolution. The simulation indicates that all the reactions controlling the pH evolution, including cation exchanges, mineral precipitation and dissolution, are coupled.

Reactive transport modeling of nitrogen in Seine River sediments

Science.gov (United States)

Akbarzadeh, Z.; Laverman, A.; Raimonet, M.; Rezanezhad, F.; Van Cappellen, P.

2016-02-01

Biogeochemical processes in sediments have a major impact on the fate and transport of nitrogen (N) in river systems. Organic matter decomposition in bottom sediments releases inorganic N species back to the stream water, while denitrification, anammox and burial of organic matter remove bioavailable N from the aquatic environment. To simulate N cycling in river sediments, a multi-component reactive transport model has been developed in MATLAB®. The model includes 3 pools of particulate organic N, plus pore water nitrate, nitrite, nitrous oxide and ammonium. Special attention is given to the production and consumption of nitrite, a N species often neglected in early diagenetic models. Although nitrite is usually considered to be short-lived, elevated nitrite concentrations have been observed in freshwater streams, raising concerns about possible toxic effects. We applied the model to sediment data sets collected at two locations in the Seine River, one upstream, the other downstream, of the largest wastewater treatment plant (WWTP) of the Paris conurbation. The model is able to reproduce the key features of the observed pore water depth profiles of the different nitrogen species. The modeling results show that the presence of oxygen in the overlying water plays a major role in controlling the exchanges of nitrite between the sediments and the stream water. In August 2012, sediments upstream of the WWTP switch from being a sink to a source of nitrite as the overlying water becomes anoxic. Downstream sediments remain a nitrite sink in oxic and anoxic conditions. Anoxic bottom waters at the upstream location promote denitrification, which produces nitrite, while at the downstream site, anammox and DNRA are important removal processes of nitrite.

Modelling freight transport

NARCIS (Netherlands)

Tavasszy, L.A.; Jong, G. de

2014-01-01

Freight Transport Modelling is a unique new reference book that provides insight into the state-of-the-art of freight modelling. Focusing on models used to support public transport policy analysis, Freight Transport Modelling systematically introduces the latest freight transport modelling

Projecting the long-term biogeochemical impacts of a diverse agroforestry system in the Midwest

Science.gov (United States)

Wolz, K. J.; DeLucia, E. H.; Paul, R. F.

2014-12-01

Annual, monoculture cropping systems have become the standard agricultural model in the Midwestern US. Unintended consequences of these systems include surface and groundwater pollution, greenhouse gas emissions, loss of biodiversity, and soil erosion. Diverse agroforestry (DA) systems dominated by fruit and nut trees/shrubs have been proposed as an agricultural model for the Midwestern US that can restore ecosystem services while simultaneously providing economically viable and industrially relevant staple food crops. A DA system including six species of fruit and nut crops was established on long-time conventional agricultural land at the University of Illinois at Urbana-Champaign in 2012, with the conventional corn-soybean rotation (CSR) as a control. Initial field measurements of the nitrogen and water cycles during the first two years of transition have indicated a significant decrease in N losses and modification of the seasonal evapotranspiration (ET) pattern. While these early results suggest that the land use transition from CSR to DA can have positive biogeochemical consequences, models must be utilized to make long-term biogeochemical projections in agroforestry systems. Initial field measurements of plant phenology, net N2O flux, nitrate leaching, soil respiration, and soil moisture were used to parameterize the DA system within the DayCENT biogeochemical model as the "savanna" ecosystem type. The model was validated with an independent subset of field measurements and then run to project biogeochemical cycling in the DA system for 25 years past establishment. Model results show that N losses via N2O emission or nitrate leaching reach a minimum within the first 5 years and then maintain this tight cycle into the future. While early ET field measurements revealed similar magnitudes between the DA and CSR systems, modeled ET continued to increase for the DA system throughout the projected time since the trees would continue to grow larger. These modeling

Effects of Model Resolution and Ocean Mixing on Forced Ice-Ocean Physical and Biogeochemical Simulations Using Global and Regional System Models

Science.gov (United States)

Jin, Meibing; Deal, Clara; Maslowski, Wieslaw; Matrai, Patricia; Roberts, Andrew; Osinski, Robert; Lee, Younjoo J.; Frants, Marina; Elliott, Scott; Jeffery, Nicole; Hunke, Elizabeth; Wang, Shanlin

2018-01-01

The current coarse-resolution global Community Earth System Model (CESM) can reproduce major and large-scale patterns but is still missing some key biogeochemical features in the Arctic Ocean, e.g., low surface nutrients in the Canada Basin. We incorporated the CESM Version 1 ocean biogeochemical code into the Regional Arctic System Model (RASM) and coupled it with a sea-ice algal module to investigate model limitations. Four ice-ocean hindcast cases are compared with various observations: two in a global 1° (40˜60 km in the Arctic) grid: G1deg and G1deg-OLD with/without new sea-ice processes incorporated; two on RASM's 1/12° (˜9 km) grid R9km and R9km-NB with/without a subgrid scale brine rejection parameterization which improves ocean vertical mixing under sea ice. Higher-resolution and new sea-ice processes contributed to lower model errors in sea-ice extent, ice thickness, and ice algae. In the Bering Sea shelf, only higher resolution contributed to lower model errors in salinity, nitrate (NO3), and chlorophyll-a (Chl-a). In the Arctic Basin, model errors in mixed layer depth (MLD) were reduced 36% by brine rejection parameterization, 20% by new sea-ice processes, and 6% by higher resolution. The NO3 concentration biases were caused by both MLD bias and coarse resolution, because of excessive horizontal mixing of high NO3 from the Chukchi Sea into the Canada Basin in coarse resolution models. R9km showed improvements over G1deg on NO3, but not on Chl-a, likely due to light limitation under snow and ice cover in the Arctic Basin.

MOPS-1.0: towards a model for the regulation of the global oceanic nitrogen budget by marine biogeochemical processes

Directory of Open Access Journals (Sweden)

I. Kriest

2015-09-01

Analysis of the model misfit with respect to observed biogeochemical tracer distributions and fluxes suggests a particle flux profile close to the one suggested by Martin et al. (1987. Simulated pelagic denitrification best agrees with the lower values between 59 and 84 Tg N yr−1 recently estimated by other authors.

Microbial and Organic Fine Particle Transport Dynamics in Streams - a Combined Experimental and Stochastic Modeling Approach

Science.gov (United States)

Drummond, Jen; Davies-Colley, Rob; Stott, Rebecca; Sukias, James; Nagels, John; Sharp, Alice; Packman, Aaron

2014-05-01

Transport dynamics of microbial cells and organic fine particles are important to stream ecology and biogeochemistry. Cells and particles continuously deposit and resuspend during downstream transport owing to a variety of processes including gravitational settling, interactions with in-stream structures or biofilms at the sediment-water interface, and hyporheic exchange and filtration within underlying sediments. Deposited cells and particles are also resuspended following increases in streamflow. Fine particle retention influences biogeochemical processing of substrates and nutrients (C, N, P), while remobilization of pathogenic microbes during flood events presents a hazard to downstream uses such as water supplies and recreation. We are conducting studies to gain insights into the dynamics of fine particles and microbes in streams, with a campaign of experiments and modeling. The results improve understanding of fine sediment transport, carbon cycling, nutrient spiraling, and microbial hazards in streams. We developed a stochastic model to describe the transport and retention of fine particles and microbes in rivers that accounts for hyporheic exchange and transport through porewaters, reversible filtration within the streambed, and microbial inactivation in the water column and subsurface. This model framework is an advance over previous work in that it incorporates detailed transport and retention processes that are amenable to measurement. Solute, particle, and microbial transport were observed both locally within sediment and at the whole-stream scale. A multi-tracer whole-stream injection experiment compared the transport and retention of a conservative solute, fluorescent fine particles, and the fecal indicator bacterium Escherichia coli. Retention occurred within both the underlying sediment bed and stands of submerged macrophytes. The results demonstrate that the combination of local measurements, whole-stream tracer experiments, and advanced modeling

A Study on the Evaluation Techniques of Pollutant Transport in Aquatic Environment

International Nuclear Information System (INIS)

Suh, Kyung Suk; Kim, Jin Hyuk; Kim, Ki Chel; Park, Geon Hyeong

2010-06-01

A study on the numerical modeling has been performed to evaluate the characteristics of hydrodynamics and pollutant transport in aquatic environment such as river, lake, estuary and coastal areas. A EFDC among the lots of codes was selected to analyze the characteristics of code and the applications. A EFDC is a general-purpose modeling package for simulating three-dimensional flow, transport, and biogeochemical processes in surface water systems including rivers, lakes, estuaries, reservoirs, wetlands and coastal regions. The structure of the EFDC model includes four major modules: a hydrodynamic model, a water quality model, a sediment transport model, and a toxic model. This report was described the characteristics of the EFDC model and examples of application in lake, river, estuary and ocean

Modeled Watershed Runoff Associated with Variations in Precipitation Data with Implications for Contaminant Fluxes

Science.gov (United States)

Watershed-scale fate and transport models are important tools for estimating the sources, transformation, and transport of contaminants to surface water systems. Precipitation is one of the primary inputs to watershed biogeochemical models, influencing changes in the water budge...

Targeted quantification of functional enzyme dynamics in environmental samples for microbially mediated biogeochemical processes: Targeted quantification of functional enzyme dynamics

Energy Technology Data Exchange (ETDEWEB)

Li, Minjing [School of Environmental Studies, China University of Geosciences, Wuhan 430074 People' s Republic of China; Gao, Yuqian [Pacific Northwest National Laboratory, Richland, WA 99354 USA; Qian, Wei-Jun [Pacific Northwest National Laboratory, Richland, WA 99354 USA; Shi, Liang [Pacific Northwest National Laboratory, Richland, WA 99354 USA; Liu, Yuanyuan [Pacific Northwest National Laboratory, Richland, WA 99354 USA; Nelson, William C. [Pacific Northwest National Laboratory, Richland, WA 99354 USA; Nicora, Carrie D. [Pacific Northwest National Laboratory, Richland, WA 99354 USA; Resch, Charles T. [Pacific Northwest National Laboratory, Richland, WA 99354 USA; Thompson, Christopher [Pacific Northwest National Laboratory, Richland, WA 99354 USA; Yan, Sen [School of Environmental Studies, China University of Geosciences, Wuhan 430074 People' s Republic of China; Fredrickson, James K. [Pacific Northwest National Laboratory, Richland, WA 99354 USA; Zachara, John M. [Pacific Northwest National Laboratory, Richland, WA 99354 USA; Liu, Chongxuan [Pacific Northwest National Laboratory, Richland, WA 99354 USA; School of Environmental Science and Engineering, Southern University of Science and Technology, Shenzhen 518055 People' s Republic of China

2017-07-13

Microbially mediated biogeochemical processes are catalyzed by enzymes that control the transformation of carbon, nitrogen, and other elements in environment. The dynamic linkage between enzymes and biogeochemical species transformation has, however, rarely been investigated because of the lack of analytical approaches to efficiently and reliably quantify enzymes and their dynamics in soils and sediments. Herein, we developed a signature peptide-based technique for sensitively quantifying dissimilatory and assimilatory enzymes using nitrate-reducing enzymes in a hyporheic zone sediment as an example. Moreover, the measured changes in enzyme concentration were found to correlate with the nitrate reduction rate in a way different from that inferred from biogeochemical models based on biomass or functional genes as surrogates for functional enzymes. This phenomenon has important implications for understanding and modeling the dynamics of microbial community functions and biogeochemical processes in environments. Our results also demonstrate the importance of enzyme quantification for the identification and interrogation of those biogeochemical processes with low metabolite concentrations as a result of faster enzyme-catalyzed consumption of metabolites than their production. The dynamic enzyme behaviors provide a basis for the development of enzyme-based models to describe the relationship between the microbial community and biogeochemical processes.

The Effects of Chlorophyll Assimilation on Carbon Fluxes in a Global Biogeochemical Model. [Technical Report Series on Global Modeling and Data Assimilation

Science.gov (United States)

Koster, Randal D. (Editor); Rousseaux, Cecile Severine; Gregg, Watson W.

2014-01-01

In this paper, we investigated whether the assimilation of remotely-sensed chlorophyll data can improve the estimates of air-sea carbon dioxide fluxes (FCO2). Using a global, established biogeochemical model (NASA Ocean Biogeochemical Model, NOBM) for the period 2003-2010, we found that the global FCO2 values produced in the free-run and after assimilation were within -0.6 mol C m(sup -2) y(sup -1) of the observations. The effect of satellite chlorophyll assimilation was assessed in 12 major oceanographic regions. The region with the highest bias was the North Atlantic. Here the model underestimated the fluxes by 1.4 mol C m(sup -2) y(sup -1) whereas all the other regions were within 1 mol C m(sup -2) y(sup -1) of the data. The FCO2 values were not strongly impacted by the assimilation, and the uncertainty in FCO2 was not decreased, despite the decrease in the uncertainty in chlorophyll concentration. Chlorophyll concentrations were within approximately 25% of the database in 7 out of the 12 regions, and the assimilation improved the chlorophyll concentration in the regions with the highest bias by 10-20%. These results suggest that the assimilation of chlorophyll data does not considerably improve FCO2 estimates and that other components of the carbon cycle play a role that could further improve our FCO2 estimates.

Multimillennium changes in dissolved oxygen under global warming: results from an AOGCM and offline ocean biogeochemical model

Science.gov (United States)

Yamamoto, A.; Abe-Ouchi, A.; Shigemitsu, M.; Oka, A.; Takahashi, K.; Ohgaito, R.; Yamanaka, Y.

2016-12-01

Long-term oceanic oxygen change due to global warming is still unclear; most future projections (such as CMIP5) are only performed until 2100. Indeed, few previous studies using conceptual models project oxygen change in the next thousands of years, showing persistent global oxygen reduction by about 30% in the next 2000 years, even after atmospheric carbon dioxide stops rising. Yet, these models cannot sufficiently represent the ocean circulation change: the key driver of oxygen change. Moreover, considering serious effect oxygen reduction has on marine life and biogeochemical cycling, long-term oxygen change should be projected for higher validity. Therefore, we used a coupled atmosphere-ocean general circulation model (AOGCM) and an offline ocean biogeochemical model, investigating realistic long-term changes in oceanic oxygen concentration and ocean circulation. We integrated these models for 2000 years under atmospheric CO2 doubling and quadrupling. After global oxygen reduction in the first 500 years, oxygen concentration in deep ocean globally recovers and overshoots, despite surface oxygen decrease and weaker Atlantic Meridional Overturning Circulation. Deep ocean convection in the Weddell Sea recovers and overshoots, after initial cessation. Thus, enhanced deep convection and associated Antarctic Bottom Water supply oxygen-rich surface waters to deep ocean, resulting global deep ocean oxygenation. We conclude that the change in ocean circulation in the Southern Ocean potentially drives millennial-scale oxygenation in the deep ocean; contrary to past reported long-term oxygen reduction and general expectation. In presentation, we will discuss the mechanism of response of deep ocean convection in the Weddell Sea and show the volume changes of hypoxic waters.

Biogeochemical modelling vs. tree-ring data - comparison of forest ecosystem productivity estimates

Science.gov (United States)

Zorana Ostrogović Sever, Maša; Barcza, Zoltán; Hidy, Dóra; Paladinić, Elvis; Kern, Anikó; Marjanović, Hrvoje

2017-04-01

Forest ecosystems are sensitive to environmental changes as well as human-induce disturbances, therefore process-based models with integrated management modules represent valuable tool for estimating and forecasting forest ecosystem productivity under changing conditions. Biogeochemical model Biome-BGC simulates carbon, nitrogen and water fluxes, and it is widely used for different terrestrial ecosystems. It was modified and parameterised by many researchers in the past to meet the specific local conditions. In this research, we used recently published improved version of the model Biome-BGCMuSo (BBGCMuSo), with multilayer soil module and integrated management module. The aim of our research is to validate modelling results of forest ecosystem productivity (NPP) from BBGCMuSo model with observed productivity estimated from an extensive dataset of tree-rings. The research was conducted in two distinct forest complexes of managed Pedunculate oak in SE Europe (Croatia), namely Pokupsko basin and Spačva basin. First, we parameterized BBGCMuSo model at a local level using eddy-covariance (EC) data from Jastrebarsko EC site. Parameterized model was used for the assessment of productivity on a larger scale. Results of NPP assessment with BBGCMuSo are compared with NPP estimated from tree ring data taken from trees on over 100 plots in both forest complexes. Keywords: Biome-BGCMuSo, forest productivity, model parameterization, NPP, Pedunculate oak

Biogeochemical Attributes That Affect the Fate and Transport of Military Relevant Contaminants Under Freeze-thaw Conditions

Science.gov (United States)

LeMonte, J.; Price, C. L.; Seiter, J.; Crocker, F. H.; Douglas, T.; Chappell, M. A.

2017-12-01

The roles and missions that the U.S. Department of Defense (DoD) undertakes in the Arctic are being reshaped by significant changes in the operational environment as a result of rising global temperatures and increased development of the vast training ranges available in Alaska. The Arctic is warming faster than any other region on Earth resulting in changing seasonality and precipitation patterns that, in turn, are leading to alterations in above ground vegetation, permafrost stability and summer sea ice extent. Collectively, these poorly defined ecosystem changes play critical roles in affecting the transport and eventual fate of persistent military relevant contaminants through unique Arctic and Subarctic terrestrial environments. As a result, management of military contaminants in a changing Arctic represents a unique and potentially significant liability to the Army and the DoD. The United States footprint in the Arctic region falls within the state of Alaska and U.S. Army Alaska manages 10% of all active Army training lands worldwide, which cover nearly 2,500 square miles in total land area. Primary recalcitrant contaminants of concern at active training ranges and at legacy sites include energetics (i.e. RDX and 2,4-dinitrotoluene) and heavy metals (i.e. antimony and lead). Through a series of field sampling and laboratory experiments, the objectives of this work are to: 1) quantify soil biogeochemical attributes that effect the physical fate and transport of military relevant contaminants in Arctic and subarctic soils under freeze-thaw conditions with a focus on near surface processes, and 2) quantify microbial diversity in Arctic and subarctic soils and the environmental constraints on community activity while exploring the effects of amendments on community function as they relate to contaminant transformation.

Skill assessment of the coupled physical-biogeochemical operational Mediterranean Forecasting System

Science.gov (United States)

Cossarini, Gianpiero; Clementi, Emanuela; Salon, Stefano; Grandi, Alessandro; Bolzon, Giorgio; Solidoro, Cosimo

2016-04-01

The Mediterranean Monitoring and Forecasting Centre (Med-MFC) is one of the regional production centres of the European Marine Environment Monitoring Service (CMEMS-Copernicus). Med-MFC operatively manages a suite of numerical model systems (3DVAR-NEMO-WW3 and 3DVAR-OGSTM-BFM) that provides gridded datasets of physical and biogeochemical variables for the Mediterranean marine environment with a horizontal resolution of about 6.5 km. At the present stage, the operational Med-MFC produces ten-day forecast: daily for physical parameters and bi-weekly for biogeochemical variables. The validation of the coupled model system and the estimate of the accuracy of model products are key issues to ensure reliable information to the users and the downstream services. Product quality activities at Med-MFC consist of two levels of validation and skill analysis procedures. Pre-operational qualification activities focus on testing the improvement of the quality of a new release of the model system and relays on past simulation and historical data. Then, near real time (NRT) validation activities aim at the routinely and on-line skill assessment of the model forecast and relays on the NRT available observations. Med-MFC validation framework uses both independent (i.e. Bio-Argo float data, in-situ mooring and vessel data of oxygen, nutrients and chlorophyll, moored buoys, tide-gauges and ADCP of temperature, salinity, sea level and velocity) and semi-independent data (i.e. data already used for assimilation, such as satellite chlorophyll, Satellite SLA and SST and in situ vertical profiles of temperature and salinity from XBT, Argo and Gliders) We give evidence that different variables (e.g. CMEMS-products) can be validated at different levels (i.e. at the forecast level or at the level of model consistency) and at different spatial and temporal scales. The fundamental physical parameters temperature, salinity and sea level are routinely validated on daily, weekly and quarterly base

Patterns of Transcript Abundance of Eukaryotic Biogeochemically-Relevant Genes in the Amazon River Plume.

Directory of Open Access Journals (Sweden)

Brian L Zielinski

Full Text Available The Amazon River has the largest discharge of all rivers on Earth, and its complex plume system fuels a wide array of biogeochemical processes, across a large area of the western tropical North Atlantic. The plume thus stimulates microbial processes affecting carbon sequestration and nutrient cycles at a global scale. Chromosomal gene expression patterns of the 2.0 to 156 μm size-fraction eukaryotic microbial community were investigated in the Amazon River Plume, generating a robust dataset (more than 100 million mRNA sequences that depicts the metabolic capabilities and interactions among the eukaryotic microbes. Combining classical oceanographic field measurements with metatranscriptomics yielded characterization of the hydrographic conditions simultaneous with a quantification of transcriptional activity and identity of the community. We highlight the patterns of eukaryotic gene expression for 31 biogeochemically significant gene targets hypothesized to be valuable within forecasting models. An advantage to this targeted approach is that the database of reference sequences used to identify the target genes was selectively constructed and highly curated optimizing taxonomic coverage, throughput, and the accuracy of annotations. A coastal diatom bloom highly expressed nitrate transporters and carbonic anhydrase presumably to support high growth rates and enhance uptake of low levels of dissolved nitrate and CO2. Diatom-diazotroph association (DDA: diatoms with nitrogen fixing symbionts blooms were common when surface salinity was mesohaline and dissolved nitrate concentrations were below detection, and hence did not show evidence of nitrate utilization, suggesting they relied on ammonium transporters to aquire recently fixed nitrogen. These DDA blooms in the outer plume had rapid turnover of the photosystem D1 protein presumably caused by photodegradation under increased light penetration in clearer waters, and increased expression of silicon

Diffusion Dominant Solute Transport Modelling in Fractured Media Under Deep Geological Environment - 12211

Energy Technology Data Exchange (ETDEWEB)

Kwong, S. [National Nuclear Laboratory (United Kingdom); Jivkov, A.P. [Research Centre for Radwaste and Decommissioning and Modelling and Simulation Centre, University of Manchester (United Kingdom)

2012-07-01

Deep geologic disposal of high activity and long-lived radioactive waste is gaining increasing support in many countries, where suitable low permeability geological formation in combination with engineered barriers are used to provide long term waste contaminant and minimise the impacts to the environment and risk to the biosphere. This modelling study examines the solute transport in fractured media under low flow velocities that are relevant to a deep geological environment. In particular, reactive solute transport through fractured media is studied using a 2-D model, that considers advection and diffusion, to explore the coupled effects of kinetic and equilibrium chemical processes. The effects of water velocity in the fracture, matrix porosity and diffusion on solute transport are investigated and discussed. Some illustrative modelled results are presented to demonstrate the use of the model to examine the effects of media degradation on solute transport, under the influences of hydrogeological (diffusion dominant) and microbially mediated chemical processes. The challenges facing the prediction of long term degradation such as cracks evolution, interaction and coalescence are highlighted. The potential of a novel microstructure informed modelling approach to account for these effects is discussed, particularly with respect to investigating multiple phenomena impact on material performance. The GRM code is used to examine the effects of media degradation for a geological waste disposal package, under the combined hydrogeological (diffusion dominant) and chemical effects in low groundwater flow conditions that are typical of deep geological disposal systems. An illustrative reactive transport modelling application demonstrates the use of the code to examine the interplay of kinetic controlled biogeochemical reactive processes with advective and diffusive transport, under the influence of media degradation. The initial model results are encouraging which show the

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

Science.gov (United States)

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...

Capturing optically important constituents and properties in a marine biogeochemical and ecosystem model

Science.gov (United States)

Dutkiewicz, S.; Hickman, A. E.; Jahn, O.; Gregg, W. W.; Mouw, C. B.; Follows, M. J.

2015-07-01

We present a numerical model of the ocean that couples a three-stream radiative transfer component with a marine biogeochemical-ecosystem component in a dynamic three-dimensional physical framework. The radiative transfer component resolves the penetration of spectral irradiance as it is absorbed and scattered within the water column. We explicitly include the effect of several optically important water constituents (different phytoplankton functional types; detrital particles; and coloured dissolved organic matter, CDOM). The model is evaluated against in situ-observed and satellite-derived products. In particular we compare to concurrently measured biogeochemical, ecosystem, and optical data along a meridional transect of the Atlantic Ocean. The simulation captures the patterns and magnitudes of these data, and estimates surface upwelling irradiance analogous to that observed by ocean colour satellite instruments. We find that incorporating the different optically important constituents explicitly and including spectral irradiance was crucial to capture the variability in the depth of the subsurface chlorophyll a (Chl a) maximum. We conduct a series of sensitivity experiments to demonstrate, globally, the relative importance of each of the water constituents, as well as the crucial feedbacks between the light field, the relative fitness of phytoplankton types, and the biogeochemistry of the ocean. CDOM has proportionally more importance at attenuating light at short wavelengths and in more productive waters, phytoplankton absorption is relatively more important at the subsurface Chl a maximum, and water molecules have the greatest contribution when concentrations of other constituents are low, such as in the oligotrophic gyres. Scattering had less effect on attenuation, but since it is important for the amount and type of upwelling irradiance, it is crucial for setting sea surface reflectance. Strikingly, sensitivity experiments in which absorption by any of the

Using geochemical indicators to distinguish high biogeochemical activity in floodplain soils and sediments

Energy Technology Data Exchange (ETDEWEB)

Kenwell, Amy [Hydrologic Sciences and Engineering Program, Colorado School of Mines, 1500 Illinois St., Golden, CO 80401 (United States); Navarre-Sitchler, Alexis, E-mail: asitchle@mines.edu [Hydrologic Sciences and Engineering Program, Colorado School of Mines, 1500 Illinois St., Golden, CO 80401 (United States); Prugue, Rodrigo [Hydrologic Sciences and Engineering Program, Colorado School of Mines, 1500 Illinois St., Golden, CO 80401 (United States); Spear, John R. [Department of Civil and Environmental Engineering, Colorado School of Mines, 1500 Illinois St., Golden, CO 80401 (United States); Hering, Amanda S. [Department of Applied Mathematics and Statistics, Colorado School of Mines, 1500 Illinois St., Golden, CO 80401 (United States); Maxwell, Reed M. [Hydrologic Sciences and Engineering Program, Colorado School of Mines, 1500 Illinois St., Golden, CO 80401 (United States); Carroll, Rosemary W.H. [Desert Research Institute, Division of Hydrologic Sciences, 2215 Raggio Parkway, Reno, NV 89512 (United States); Williams, Kenneth H. [Lawrence Berkeley National Laboratory, Berkeley, CA 94720 (United States)

2016-09-01

A better understanding of how microbial communities interact with their surroundings in physically and chemically heterogeneous subsurface environments will lead to improved quantification of biogeochemical reactions and associated nutrient cycling. This study develops a methodology to predict potential elevated rates of biogeochemical activity (microbial “hotspots”) in subsurface environments by correlating microbial DNA and aspects of the community structure with the spatial distribution of geochemical indicators in subsurface sediments. Multiple linear regression models of simulated precipitation leachate, HCl and hydroxylamine extractable iron and manganese, total organic carbon (TOC), and microbial community structure were used to identify sample characteristics indicative of biogeochemical hotspots within fluvially-derived aquifer sediments and overlying soils. The method has been applied to (a) alluvial materials collected at a former uranium mill site near Rifle, Colorado and (b) relatively undisturbed floodplain deposits (soils and sediments) collected along the East River near Crested Butte, Colorado. At Rifle, 16 alluvial samples were taken from 8 sediment cores, and at the East River, 46 soil/sediment samples were collected across and perpendicular to 3 active meanders and an oxbow meander. Regression models using TOC and TOC combined with extractable iron and manganese results were determined to be the best fitting statistical models of microbial DNA (via 16S rRNA gene analysis). Fitting these models to observations in both contaminated and natural floodplain deposits, and their associated alluvial aquifers, demonstrates the broad applicability of the geochemical indicator based approach. - Highlights: • Biogeochemical characterization of alluvial floodplain soils and sediments was performed to investigate parameters that may indicate microbial hot spot formation. • A correlation between geochemical parameters (total organic carbon and

Using geochemical indicators to distinguish high biogeochemical activity in floodplain soils and sediments

International Nuclear Information System (INIS)

Kenwell, Amy; Navarre-Sitchler, Alexis; Prugue, Rodrigo; Spear, John R.; Hering, Amanda S.; Maxwell, Reed M.; Carroll, Rosemary W.H.; Williams, Kenneth H.

2016-01-01

A better understanding of how microbial communities interact with their surroundings in physically and chemically heterogeneous subsurface environments will lead to improved quantification of biogeochemical reactions and associated nutrient cycling. This study develops a methodology to predict potential elevated rates of biogeochemical activity (microbial “hotspots”) in subsurface environments by correlating microbial DNA and aspects of the community structure with the spatial distribution of geochemical indicators in subsurface sediments. Multiple linear regression models of simulated precipitation leachate, HCl and hydroxylamine extractable iron and manganese, total organic carbon (TOC), and microbial community structure were used to identify sample characteristics indicative of biogeochemical hotspots within fluvially-derived aquifer sediments and overlying soils. The method has been applied to (a) alluvial materials collected at a former uranium mill site near Rifle, Colorado and (b) relatively undisturbed floodplain deposits (soils and sediments) collected along the East River near Crested Butte, Colorado. At Rifle, 16 alluvial samples were taken from 8 sediment cores, and at the East River, 46 soil/sediment samples were collected across and perpendicular to 3 active meanders and an oxbow meander. Regression models using TOC and TOC combined with extractable iron and manganese results were determined to be the best fitting statistical models of microbial DNA (via 16S rRNA gene analysis). Fitting these models to observations in both contaminated and natural floodplain deposits, and their associated alluvial aquifers, demonstrates the broad applicability of the geochemical indicator based approach. - Highlights: • Biogeochemical characterization of alluvial floodplain soils and sediments was performed to investigate parameters that may indicate microbial hot spot formation. • A correlation between geochemical parameters (total organic carbon and

Biogeochemical cycles of Chernobyl-born radionuclides in the contaminated forest ecosystems: long-term dynamics of the migration processes

Science.gov (United States)

Shcheglov, Alexey; Tsvetnova, Ol'ga; Klyashtorin, Alexey

2013-04-01

Biogeochemical migration is a dominant factor of the radionuclide transport through the biosphere. In the early XX century, V.I. Vernadskii, a Russian scientist known, noted about a special role living things play in transport and accumulation of natural radionuclide in various environments. The role of biogeochemical processes in migration and redistribution of technogenic radionuclides is not less important. In Russia, V. M. Klechkovskii and N.V. Timofeev-Ressovskii showed some important biogeochemical aspects of radionuclide migration by the example of global fallout and Kyshtym accident. Their followers, R.M. Alexakhin, M.A. Naryshkin, N.V. Kulikov, F.A. Tikhomirov, E.B. Tyuryukanova, and others also contributed a lot to biogeochemistry of radionuclides. In the post-Chernobyl period, this area of knowledge received a lot of data that allowed building the radioactive element balance and flux estimation in various biogeochemical cycles [Shcheglov et al., 1999]. Regrettably, many of recent radioecological studies are only focused on specific radionuclide fluxes or pursue some applied tasks, missing the holistic approach. Most of the studies consider biogeochemical fluxes of radioactive isotopes in terms of either dose estimation or radionuclide migration rates in various food chains. However, to get a comprehensive picture and develop a reliable forecast of environmental, ecological, and social consequences of radioactive pollution in a vast contaminated area, it is necessary to investigate all the radionuclide fluxes associated with the biogeochemical cycles in affected ecosystems. We believe such an integrated approach would be useful to study long-term environmental consequences of the Fukushima accident as well. In our long-term research, we tried to characterize the flux dynamics of the Chernobyl-born radionuclides in the contaminated forest ecosystems and landscapes as a part of the integrated biogeochemical process. Our field studies were started in June of

Introducing mixotrophy into a biogeochemical model describing an eutrophied coastal ecosystem: The Southern North Sea

Science.gov (United States)

Ghyoot, Caroline; Lancelot, Christiane; Flynn, Kevin J.; Mitra, Aditee; Gypens, Nathalie

2017-09-01

Most biogeochemical/ecological models divide planktonic protists between phototrophs (phytoplankton) and heterotrophs (zooplankton). However, a large number of planktonic protists are able to combine several mechanisms of carbon and nutrient acquisition. Not representing these multiple mechanisms in biogeochemical/ecological models describing eutrophied coastal ecosystems can potentially lead to different conclusions regarding ecosystem functioning, especially regarding the success of harmful algae, which are often reported as mixotrophic. This modelling study investigates the implications for trophic dynamics of including 3 contrasting forms of mixotrophy, namely osmotrophy (using alkaline phosphatase activity, APA), non-constitutive mixotrophy (acquired phototrophy by microzooplankton) and also constitutive mixotrophy. The application is in the Southern North Sea, an ecosystem that faced, between 1985 and 2005, a significant increase in the nutrient supply N:P ratio (from 31 to 81 mol N:P). The comparison with a traditional model shows that, when the winter N:P ratio in the Southern North Sea is above 22 molN molP-1 (as occurred from mid-1990s), APA allows a 3-32% increase of annual gross primary production (GPP). In result of the higher GPP, the annual sedimentation increases as well as the bacterial production. By contrast, APA does not affect the export of matter to higher trophic levels because the increased GPP is mainly due to Phaeocystis colonies, which are not grazed by copepods. Under high irradiance, non-constitutive mixotrophy appreciably increases annual GPP, transfer to higher trophic levels, sedimentation, and nutrient remineralisation. In this ecosystem, non-constitutive mixotrophy is also observed to have an indirect stimulating effect on diatoms. Constitutive mixotrophy in nanoflagellates appears to have little influence on this ecosystem functioning. An important conclusion from this work is that contrasting forms of mixotrophy have different

Ecohydrological Interfaces as Dynamic Hotspots of Biogeochemical Cycling

Science.gov (United States)

Krause, Stefan; Lewandowski, Joerg; Hannah, David; McDonald, Karlie; Folegot, Silvia; Baranov, Victor

2016-04-01

Ecohydrological interfaces, represent the boundaries between water-dependent ecosystems that can alter substantially the fluxes of energy and matter. There is still a critical gap of understanding the organisational principles of the drivers and controls of spatially and temporally variable ecohydrological interface functions. This knowledge gap limits our capacity to efficiently quantify, predict and manage the services provided by complex ecosystems. Many ecohydrological interfaces are characterized by step changes in microbial metabolic activity, steep redox gradients and often even thermodynamic phase shifts, for instance at the interfaces between atmosphere and water or soil matrix and macro-pores interfaces. This paper integrates investigations from point scale laboratory microcosm experiments with reach and subcatchment scale tracer experiments and numerical modeling studies to elaborate similarities in the drivers and controls that constitute the enhanced biogeochemical activity of different types of ecohydrologica interfaces across a range of spatial and temporal scales. We therefore combine smart metabolic activity tracers to quantify the impact of bioturbating benthic fauna onto ecosystem respiration and oxygen consumption and investigate at larger scale, how microbial metabolic activity and carbon turnover at the water-sediment interface are controlled by sediment physical and chemical properties as well as water temperatures. Numerical modeling confirmed that experimentally identified hotspots of streambed biogeochemical cycling were controlled by patterns of physical properties such as hydraulic conductivities or bioavailability of organic matter, impacting on residence time distributions and hence reaction times. In contrast to previous research, our investigations thus confirmed that small-scale variability of physical and chemical interface properties had a major impact on biogeochemical processing at the investigated ecohydrological interfaces

Comparison of alternative spatial resolutions in the application of a spatially distributed biogeochemical model over complex terrain

Science.gov (United States)

Turner, D.P.; Dodson, R.; Marks, D.

1996-01-01

Spatially distributed biogeochemical models may be applied over grids at a range of spatial resolutions, however, evaluation of potential errors and loss of information at relatively coarse resolutions is rare. In this study, a georeferenced database at the 1-km spatial resolution was developed to initialize and drive a process-based model (Forest-BGC) of water and carbon balance over a gridded 54976 km2 area covering two river basins in mountainous western Oregon. Corresponding data sets were also prepared at 10-km and 50-km spatial resolutions using commonly employed aggregation schemes. Estimates were made at each grid cell for climate variables including daily solar radiation, air temperature, humidity, and precipitation. The topographic structure, water holding capacity, vegetation type and leaf area index were likewise estimated for initial conditions. The daily time series for the climatic drivers was developed from interpolations of meteorological station data for the water year 1990 (1 October 1989-30 September 1990). Model outputs at the 1-km resolution showed good agreement with observed patterns in runoff and productivity. The ranges for model inputs at the 10-km and 50-km resolutions tended to contract because of the smoothed topography. Estimates for mean evapotranspiration and runoff were relatively insensitive to changing the spatial resolution of the grid whereas estimates of mean annual net primary production varied by 11%. The designation of a vegetation type and leaf area at the 50-km resolution often subsumed significant heterogeneity in vegetation, and this factor accounted for much of the difference in the mean values for the carbon flux variables. Although area wide means for model outputs were generally similar across resolutions, difference maps often revealed large areas of disagreement. Relatively high spatial resolution analyses of biogeochemical cycling are desirable from several perspectives and may be particularly important in the

Characteristics and Source Apportionment of Marine Aerosols over East China Sea Using a Source-oriented Chemical Transport Model

Science.gov (United States)

Kang, M.; Zhang, H.; Fu, P.

2017-12-01

Marine aerosols exert a strong influence on global climate change and biogeochemical cycling, as oceans cover beyond 70% of the Earth's surface. However, investigations on marine aerosols are relatively limited at present due to the difficulty and inconvenience in sampling marine aerosols as well as their diverse sources. East China Sea (ECS), lying over the broad shelf of the western North Pacific, is adjacent to the Asian mainland, where continental-scale air pollution could impose a heavy load on the marine atmosphere through long-range atmospheric transport. Thus, contributions of major sources to marine aerosols need to be identified for policy makers to develop cost effective control strategies. In this work, a source-oriented version of the Community Multiscale Air Quality (CMAQ) model, which can directly track the contributions from multiple emission sources to marine aerosols, is used to investigate the contributions from power, industry, transportation, residential, biogenic and biomass burning to marine aerosols over the ECS in May and June 2014. The model simulations indicate significant spatial and temporal variations of concentrations as well as the source contributions. This study demonstrates that the Asian continent can greatly affect the marine atmosphere through long-range transport.

Hyporheic zone as a bioreactor: sediment heterogeneity influencing biogeochemical processes

Science.gov (United States)

Perujo, Nuria; Romani, Anna M.; Sanchez-Vila, Xavier

2017-04-01

Mediterranean fluvial systems are characterized by frequent periods of low flow or even drought. During low flow periods, water from wastewater treatment plants (WWTPs) is proportionally large in fluvial systems. River water might be vertically transported through the hyporheic zone, and then porous medium acts as a complementary treatment system since, as water infiltrates, a suite of biogeochemical processes occurs. Subsurface sediment heterogeneity plays an important role since it influences the interstitial fluxes of the medium and drives biomass growing, determining biogeochemical reactions. In this study, WWTP water was continuously infiltrated for 3 months through two porous medium tanks: one consisting of 40 cm of fine sediment (homogeneous); and another comprised of two layers of different grain size sediments (heterogeneous), 20 cm of coarse sediment in the upper part and 20 cm of fine one in the bottom. Several hydrological, physicochemical and biological parameters were measured periodically (weekly at the start of the experiment and biweekly at the end). Analysed parameters include dissolved nitrogen, phosphorus, organic carbon, and oxygen all measured at the surface, and at 5, 20 and 40 cm depth. Variations in hydraulic conductivity with time were evaluated. Sediment samples were also analysed at three depths (surface, 20 and 40 cm) to determine bacterial density, chlorophyll content, extracellular polymeric substances, and biofilm function (extracellular enzyme activities and carbon substrate utilization profiles). Preliminary results suggest hydraulic conductivity to be the main driver of the differences in the biogeochemical processes occurring in the subsurface. At the heterogeneous tank, a low nutrient reduction throughout the whole medium is measured. In this medium, high hydraulic conductivity allows for a large amount of infiltrating water, but with a small residence time. Since some biological processes are largely time-dependent, small water

Diffusion Dominant Solute Transport Modelling In Deep Repository Under The Effect of Emplacement Media Degradation - 13285

International Nuclear Information System (INIS)

Kwong, S.; Jivkov, A.P.

2013-01-01

Deep geologic disposal of high activity and long-lived radioactive waste is being actively considered and pursued in many countries, where low permeability geological formations are used to provide long term waste contaminant with minimum impact to the environment and risk to the biosphere. A multi-barrier approach that makes use of both engineered and natural barriers (i.e. geological formations) is often used to further enhance the containment performance of the repository. As the deep repository system subjects to a variety of thermo-hydro-chemo-mechanical (THCM) effects over its long 'operational' lifespan (e.g. 0.1 to 1.0 million years, the integrity of the barrier system will decrease over time (e.g. fracturing in rock or clay)). This is broadly referred as media degradation in the present study. This modelling study examines the effects of media degradation on diffusion dominant solute transport in fractured media that are typical of deep geological environment. In particular, reactive solute transport through fractured media is studied using a 2-D model, that considers advection and diffusion, to explore the coupled effects of kinetic and equilibrium chemical processes, while the effects of degradation is studied using a pore network model that considers the media diffusivity and network changes. Model results are presented to demonstrate the use of a 3D pore-network model, using a novel architecture, to calculate macroscopic properties of the medium such as diffusivity, subject to pore space changes as the media degrade. Results from a reactive transport model of a representative geological waste disposal package are also presented to demonstrate the effect of media property change on the solute migration behaviour, illustrating the complex interplay between kinetic biogeochemical processes and diffusion dominant transport. The initial modelling results demonstrate the feasibility of a coupled modelling approach (using pore-network model and reactive

Ecotoxicological, ecophysiological and biogeochemical fundamentals of risk assessment

International Nuclear Information System (INIS)

Bashkin, V.; Evstafjeva, E.

1995-01-01

A quantitative risk assessment (RA) for complex influence of different factors in heavy polluted regions is possible to carry out only on a basis of determination of various links of biogeochemical trophical chains and analysis of the whole biogeochemical structure of the region under study. As an integrative assessment, the human adaptability should be chosen because the majority of trophical chains are closed by man. The given integrative criteria includes biogeochemical, ecophysiological and ecotoxicological assessment of risk factors. Consequently, ecological-biogeochemical regionalization, ecophysiological and ecotoxicological monitoring of human population health are the important approaches to RA. These criteria should be conjugated with LCA of various industrial and agricultural products. At the ultimate degree, the given approaches are needed for areas where traditional pollutants (heavy metals, POPS, pesticides, fertilizers) are enforced sharply by radioactive pollution. Due to the complex influence of pollutants, it is impossible to use individual guidelines. For RA of these complex pollutants, the methods of human adaptability assessment to a polluted environment have to be carried out. These methods include biogeochemical, ecotoxicological and ecophysiological analysis of risk factors as well as quantitative uncertainty analysis. Furthermore, the modern statistical methods such as correlative graphs etc., have to be used for quantitative assessment of human adaptability to complex influence of pollutants. The results obtained in the Chernobyl region have shown the acceptability of suggested methods

Disturbance decouples biogeochemical cycles across forests of the southeastern US

Science.gov (United States)

Ashley D. Keiser; Jennifer D. Knoepp; Mark A. Bradford

2016-01-01

Biogeochemical cycles are inherently linked through the stoichiometric demands of the organisms that cycle the elements. Landscape disturbance can alter element availability and thus the rates of biogeochemical cycling. Nitrification is a fundamental biogeochemical process positively related to plant productivity and nitrogen loss from soils to aquatic systems, and the...

Improving National Capability in Biogeochemical Flux Modelling: the UK Environmental Virtual Observatory (EVOp)

Science.gov (United States)

Johnes, P.; Greene, S.; Freer, J. E.; Bloomfield, J.; Macleod, K.; Reaney, S. M.; Odoni, N. A.

2012-12-01

The best outcomes from watershed management arise where policy and mitigation efforts are underpinned by strong science evidence, but there are major resourcing problems associated with the scale of monitoring needed to effectively characterise the sources rates and impacts of nutrient enrichment nationally. The challenge is to increase national capability in predictive modelling of nutrient flux to waters, securing an effective mechanism for transferring knowledge and management tools from data-rich to data-poor regions. The inadequacy of existing tools and approaches to address these challenges provided the motivation for the Environmental Virtual Observatory programme (EVOp), an innovation from the UK Natural Environment Research Council (NERC). EVOp is exploring the use of a cloud-based infrastructure in catchment science, developing an exemplar to explore N and P fluxes to inland and coastal waters in the UK from grid to catchment and national scale. EVOp is bringing together for the first time national data sets, models and uncertainty analysis into cloud computing environments to explore and benchmark current predictive capability for national scale biogeochemical modelling. The objective is to develop national biogeochemical modelling capability, capitalising on extensive national investment in the development of science understanding and modelling tools to support integrated catchment management, and supporting knowledge transfer from data rich to data poor regions, The AERC export coefficient model (Johnes et al., 2007) has been adapted to function within the EVOp cloud environment, and on a geoclimatic basis, using a range of high resolution, geo-referenced digital datasets as an initial demonstration of the enhanced national capacity for N and P flux modelling using cloud computing infrastructure. Geoclimatic regions are landscape units displaying homogenous or quasi-homogenous functional behaviour in terms of process controls on N and P cycling

Investigating the Role of Biogeochemical Processes in the Northern High Latitudes on Global Climate Feedbacks Using an Efficient Scalable Earth System Model

Energy Technology Data Exchange (ETDEWEB)

Jain, Atul K. [Univ. of Illinois, Urbana-Champaign, IL (United States)

2016-09-14

The overall objectives of this DOE funded project is to combine scientific and computational challenges in climate modeling by expanding our understanding of the biogeophysical-biogeochemical processes and their interactions in the northern high latitudes (NHLs) using an earth system modeling (ESM) approach, and by adopting an adaptive parallel runtime system in an ESM to achieve efficient and scalable climate simulations through improved load balancing algorithms.

The significance of GW-SW interactions for biogeochemical processes in sandy streambeds

Science.gov (United States)

Arnon, Shai; De Falco, Natalie; Fox, Aryeh; Laube, Gerrit; Schmidt, Christian; Fleckenstein, Jan; Boano, Fulvio

2015-04-01

Stream-groundwater interactions have a major impact on hyporheic exchange fluxes in sandy streambeds. However, the physical complexity of natural streams has limited our ability to study these types of interactions systematically, and to evaluate their importance to biogeochemical processes and nutrient cycling. In this work we were able to quantify the effect of losing and gaining fluxes on hyporheic exchange and nutrient cycling in homogeneous and heterogeneous streambeds by combining experiments in laboratory flumes and modeling. Tracer experiments for measuring hyporheic exchange were done using dyes and NaCl under various combinations of overlying water velocity and losing or gaining fluxes. Nutrient cycling experiments were conducted after growing a benthic biofilm by spiking with Sodium Benzoate (as a source of labile dissolved organic carbon, DOC) and measuring DOC and oxygen dynamics. The combination of experimental observations and modeling revealed that interfacial transport increases with the streambed hydraulic conductivity and proportional to the square of the overlying water velocity. Hyporheic exchange fluxes under losing and gaining flow conditions were similar, and became smaller when the losing or gaining flux increases. Increasing in streambed hydraulic conductivity led to higher hyporheic fluxes and reduction in the effects of losing and gaining flow conditions to constrain exchange. Despite the evident effect of flow conditions on hyporheic exchange, labile DOC uptake was positively linked to increasing overlying water velocity but was not affected by losing and gaining fluxes. This is because microbial aerobic activity was taking place at the upper few millimeters of the streambed as shown by local oxygen consumption rates, which was measured using microelectrodes. Based on modeling work, it is expected that GW-SW interaction will be more significant for less labile DOC and anaerobic processes. Our results enable us to study systematically

Advanced transport modeling of toroidal plasmas with transport barriers

International Nuclear Information System (INIS)

Fukuyama, A.; Murakami, S.; Honda, M.; Izumi, Y.; Yagi, M.; Nakajima, N.; Nakamura, Y.; Ozeki, T.

2005-01-01

Transport modeling of toroidal plasmas is one of the most important issue to predict time evolution of burning plasmas and to develop control schemes in reactor plasmas. In order to describe the plasma rotation and rapid transition self-consistently, we have developed an advanced scheme of transport modeling based on dynamical transport equation and applied it to the analysis of transport barrier formation. First we propose a new transport model and examine its behavior by the use of conventional diffusive transport equation. This model includes the electrostatic toroidal ITG mode and the electromagnetic ballooning mode and successfully describes the formation of internal transport barriers. Then the dynamical transport equation is introduced to describe the plasma rotation and the radial electric field self-consistently. The formation of edge transport barriers is systematically studied and compared with experimental observations. The possibility of kinetic transport modeling in velocity space is also examined. Finally the modular structure of integrated modeling code for tokamaks and helical systems is discussed. (author)

Geophysical monitoring and reactive transport modeling of ureolytically-driven calcium carbonate precipitation

Energy Technology Data Exchange (ETDEWEB)

Wu, Y.; Ajo-Franklin, J.B.; Spycher, N.; Hubbard, S.S.; Zhang, G.; Williams, K.H.; Taylor, J.; Fujita, Y.; Smith, R.

2011-07-15

Ureolytically-driven calcium carbonate precipitation is the basis for a promising in-situ remediation method for sequestration of divalent radionuclide and trace metal ions. It has also been proposed for use in geotechnical engineering for soil strengthening applications. Monitoring the occurrence, spatial distribution, and temporal evolution of calcium carbonate precipitation in the subsurface is critical for evaluating the performance of this technology and for developing the predictive models needed for engineering application. In this study, we conducted laboratory column experiments using natural sediment and groundwater to evaluate the utility of geophysical (complex resistivity and seismic) sensing methods, dynamic synchrotron x-ray computed tomography (micro-CT), and reactive transport modeling for tracking ureolytically-driven calcium carbonate precipitation processes under site relevant conditions. Reactive transport modeling with TOUGHREACT successfully simulated the changes of the major chemical components during urea hydrolysis. Even at the relatively low level of urea hydrolysis observed in the experiments, the simulations predicted an enhanced calcium carbonate precipitation rate that was 3-4 times greater than the baseline level. Reactive transport modeling results, geophysical monitoring data and micro-CT imaging correlated well with reaction processes validated by geochemical data. In particular, increases in ionic strength of the pore fluid during urea hydrolysis predicted by geochemical modeling were successfully captured by electrical conductivity measurements and confirmed by geochemical data. The low level of urea hydrolysis and calcium carbonate precipitation suggested by the model and geochemical data was corroborated by minor changes in seismic P-wave velocity measurements and micro-CT imaging; the latter provided direct evidence of sparsely distributed calcium carbonate precipitation. Ion exchange processes promoted through NH{sub 4}{sup

User manual of the multicompenent variably - saturated flow and transport model HP1

International Nuclear Information System (INIS)

Jacques, D.; Simunek, J.

2005-06-01

This report describes a new comprehensive simulation tool HP1 (HYDRUS1D-PHREEQC) that was obtained by coupling the HYDRUS-1D one-dimensional variably-saturated water flow and solute transport model with the PHREEQC geochemical code. The HP1 code incorporates modules simulating (1) transient water flow in variably-saturated media, (2) transport of multiple components, and (3) mixed equilibrium/kinetic geochemical reactions. The program numerically solves the Richards equation for variably-saturated water flow and advection-dispersion type equations for heat and solute transport. The flow equation incorporates a sink term to account for water uptake by plant roots. The heat transport equation considers transport due to conduction and convection with flowing water. The solute transport equations consider advective-dispersive transport in the liquid phase. The program can simulate a broad range of low-temperature biogeochemical reactions in water, soil and ground water systems including interactions with minerals, gases, exchangers, and sorption surfaces, based on thermodynamic equilibrium, kinetics, or mixed equilibrium-kinetic reactions. The program may be used to analyze water and solute movement in unsaturated, partially saturated, or fully saturated porous media. The flow region may be composed of nonuniform soils or sediments. Flow and transport can occur in the vertical, horizontal, or a generally inclined direction. The water flow part of the model can deal with prescribed head and flux boundaries, boundaries controlled by atmospheric conditions, as well as free drainage boundary conditions. The governing flow and transport equations were solved numerically using Galerkin-type linear finite element schemes. To test the accuracy of the coupling procedures implemented in HP1, simulation results were compared with (i) HYDRUS-1D for transport problems of multiple components subject to sequential first-order decay, (ii) PHREEQC for steady-state flow conditions, and

The interaction of the flux errors and transport errors in modeled atmospheric carbon dioxide concentrations

Science.gov (United States)

Feng, S.; Lauvaux, T.; Butler, M. P.; Keller, K.; Davis, K. J.; Jacobson, A. R.; Schuh, A. E.; Basu, S.; Liu, J.; Baker, D.; Crowell, S.; Zhou, Y.; Williams, C. A.

2017-12-01

Regional estimates of biogenic carbon fluxes over North America from top-down atmospheric inversions and terrestrial biogeochemical (or bottom-up) models remain inconsistent at annual and sub-annual time scales. While top-down estimates are impacted by limited atmospheric data, uncertain prior flux estimates and errors in the atmospheric transport models, bottom-up fluxes are affected by uncertain driver data, uncertain model parameters and missing mechanisms across ecosystems. This study quantifies both flux errors and transport errors, and their interaction in the CO2 atmospheric simulation. These errors are assessed by an ensemble approach. The WRF-Chem model is set up with 17 biospheric fluxes from the Multiscale Synthesis and Terrestrial Model Intercomparison Project, CarbonTracker-Near Real Time, and the Simple Biosphere model. The spread of the flux ensemble members represents the flux uncertainty in the modeled CO2 concentrations. For the transport errors, WRF-Chem is run using three physical model configurations with three stochastic perturbations to sample the errors from both the physical parameterizations of the model and the initial conditions. Additionally, the uncertainties from boundary conditions are assessed using four CO2 global inversion models which have assimilated tower and satellite CO2 observations. The error structures are assessed in time and space. The flux ensemble members overall overestimate CO2 concentrations. They also show larger temporal variability than the observations. These results suggest that the flux ensemble is overdispersive. In contrast, the transport ensemble is underdispersive. The averaged spatial distribution of modeled CO2 shows strong positive biogenic signal in the southern US and strong negative signals along the eastern coast of Canada. We hypothesize that the former is caused by the 3-hourly downscaling algorithm from which the nighttime respiration dominates the daytime modeled CO2 signals and that the latter

Reanalysis of biogeochemical properties in the Mediterranean Sea

Science.gov (United States)

Cossarini, Gianpiero; Teruzzi, Anna; Salon, Stefano; Solidoro, Cosimo

2014-05-01

In the 3D variational (3DVAR) assimilation approach the error covariance matrix can be decomposed in a series of operators. The decomposition makes the 3DVAR particularly suitable for marine biogeochemistry data assimilation, because of the reduced computational costs of the method and its modularity, which allows to define the covariance among the biogeochemical variables in a specific operator. In the present work, the results of 3DVAR assimilation of surface chlorophyll concentration in a multi-annual simulation of the Mediterranean Sea biogeochemistry are presented. The assimilated chlorophyll concentrations are obtained from satellite observations (Volpe et al. 2012). The multi-annual simulation is carried out using the OPATM-BFM model (Lazzari et al. 2012), which describes the low trophic web dynamics and is offline coupled with the MFS physical model (Oddo et al. 2009). In the OPATM-BFM four types of phytoplankton are simulated in terms of their content in carbon, nitrogen, phosphorous, silicon and chlorophyll. In the 3DVAR the error covariance matrix has been decomposed in three different operators, which account for the vertical, the horizontal and the biogeochemical covariance (Teruzzi et al. 2014). The biogeochemical operator propagates the result of the assimilation to the OPATM-BFM variables, providing innovation for the components of the four phytoplankton types. The biogeochemical covariance has been designed supposing that the assimilation preserves the physiological status and the relative abundances of phytoplankton types. Practically, the assimilation preserves the internal quotas of the components for each phytoplankton as long as the optimal growth rate condition are maintained. The quotas preservation is not applied when the phytoplankton is in severe declining growth phase, and the correction provided by the assimilation is set equal to zero. Moreover, the relative abundances among the phytoplankton functional types are preserved. The 3DVAR

Reactive solute transport in acidic streams

Science.gov (United States)

Broshears, R.E.

1996-01-01

Spatial and temporal profiles of Ph and concentrations of toxic metals in streams affected by acid mine drainage are the result of the interplay of physical and biogeochemical processes. This paper describes a reactive solute transport model that provides a physically and thermodynamically quantitative interpretation of these profiles. The model combines a transport module that includes advection-dispersion and transient storage with a geochemical speciation module based on MINTEQA2. Input to the model includes stream hydrologic properties derived from tracer-dilution experiments, headwater and lateral inflow concentrations analyzed in field samples, and a thermodynamic database. Simulations reproduced the general features of steady-state patterns of observed pH and concentrations of aluminum and sulfate in St. Kevin Gulch, an acid mine drainage stream near Leadville, Colorado. These patterns were altered temporarily by injection of sodium carbonate into the stream. A transient simulation reproduced the observed effects of the base injection.

Morphological, hydrological, biogeochemical and ecological changes and challenges in river restoration - the Thur River case study

Science.gov (United States)

Schirmer, M.; Luster, J.; Linde, N.; Perona, P.; Mitchell, E. A. D.; Barry, D. A.; Hollender, J.; Cirpka, O. A.; Schneider, P.; Vogt, T.; Radny, D.; Durisch-Kaiser, E.

2014-06-01

River restoration can enhance river dynamics, environmental heterogeneity and biodiversity, but the underlying processes governing the dynamic changes need to be understood to ensure that restoration projects meet their goals, and adverse effects are prevented. In particular, we need to comprehend how hydromorphological variability quantitatively relates to ecosystem functioning and services, biodiversity as well as ground- and surface water quality in restored river corridors. This involves (i) physical processes and structural properties, determining erosion and sedimentation, as well as solute and heat transport behavior in surface water and within the subsurface; (ii) biogeochemical processes and characteristics, including the turnover of nutrients and natural water constituents; and (iii) ecological processes and indicators related to biodiversity and ecological functioning. All these aspects are interlinked, requiring an interdisciplinary investigation approach. Here, we present an overview of the recently completed RECORD (REstored CORridor Dynamics) project in which we combined physical, chemical, and biological observations with modeling at a restored river corridor of the perialpine Thur River in Switzerland. Our results show that river restoration, beyond inducing morphologic changes that reshape the river bed and banks, triggered complex spatial patterns of bank infiltration, and affected habitat type, biotic communities and biogeochemical processes. We adopted an interdisciplinary approach of monitoring the continuing changes due to restoration measures to address the following questions: How stable is the morphological variability established by restoration? Does morphological variability guarantee an improvement in biodiversity? How does morphological variability affect biogeochemical transformations in the river corridor? What are some potential adverse effects of river restoration? How is river restoration influenced by catchment-scale hydraulics

Implications of a More Comprehensive Nitrogen Cycle in a Global Biogeochemical Ocean Model

Science.gov (United States)

Six, K. D.; Ilyina, T.

2016-02-01

Nitrogen plays a crucial role for nearly all living organisms in the Earth system. Changes in the marine nitrogen cycle not only alter the marine biota, but will also have an impact on the marine carbon cycle and, in turn, on climate due to the close coupling of the carbon-nitrogen cycle. The understanding of processes and controls of the marine nitrogen cycle is therefore a prerequisite to reduce uncertainties in the prediction of future climate. Nevertheless, most ocean biogeochemical components of modern Earth system models have a rather simplistic representation of marine N-cycle mainly focusing on nitrate. Here we present results of the HAMburg Ocean Carbon Cycle model (HAMOCC) as part of the MPI-ESM which was extended by a prognostic representation of ammonium and nitrite to resolve important processes of the marine N-cycle such as nitrification and anaerobic ammonium oxidation (anammox). Additionally, we updated the production of nitrous oxide, an important greenhouse gas, allowing for two sources from oxidation of ammonium (nitrification) and from reduction of nitrite (nitrifier-denitrification) at low oxygen concentrations. Besides an extended model data comparison we discuss the following aspects of the N-cycle by model means: (1) contribution of anammox to the loss of fixed nitrogen, and (2) production and emission of marine nitrous oxide.

A Hydrological Concept including Lateral Water Flow Compatible with the Biogeochemical Model ForSAFE

Directory of Open Access Journals (Sweden)

Giuliana Zanchi

2016-03-01

Full Text Available The study presents a hydrology concept developed to include lateral water flow in the biogeochemical model ForSAFE. The hydrology concept was evaluated against data collected at Svartberget in the Vindeln Research Forest in Northern Sweden. The results show that the new concept allows simulation of a saturated and an unsaturated zone in the soil as well as water flow that reaches the stream comparable to measurements. The most relevant differences compared to streamflow measurements are that the model simulates a higher base flow in winter and lower flow peaks after snowmelt. These differences are mainly caused by the assumptions made to regulate the percolation at the bottom of the simulated soil columns. The capability for simulating lateral flows and a saturated zone in ForSAFE can greatly improve the simulation of chemical exchange in the soil and export of elements from the soil to watercourses. Such a model can help improve the understanding of how environmental changes in the forest landscape will influence chemical loads to surface waters.

UZ Colloid Transport Model

International Nuclear Information System (INIS)

McGraw, M.

2000-01-01

The UZ Colloid Transport model development plan states that the objective of this Analysis/Model Report (AMR) is to document the development of a model for simulating unsaturated colloid transport. This objective includes the following: (1) use of a process level model to evaluate the potential mechanisms for colloid transport at Yucca Mountain; (2) Provide ranges of parameters for significant colloid transport processes to Performance Assessment (PA) for the unsaturated zone (UZ); (3) Provide a basis for development of an abstracted model for use in PA calculations

A state-space Bayesian framework for estimating biogeochemical transformations using time-lapse geophysical data

Energy Technology Data Exchange (ETDEWEB)

Chen, J.; Hubbard, S.; Williams, K.; Pride, S.; Li, L.; Steefel, C.; Slater, L.

2009-04-15

We develop a state-space Bayesian framework to combine time-lapse geophysical data with other types of information for quantitative estimation of biogeochemical parameters during bioremediation. We consider characteristics of end-products of biogeochemical transformations as state vectors, which evolve under constraints of local environments through evolution equations, and consider time-lapse geophysical data as available observations, which could be linked to the state vectors through petrophysical models. We estimate the state vectors and their associated unknown parameters over time using Markov chain Monte Carlo sampling methods. To demonstrate the use of the state-space approach, we apply it to complex resistivity data collected during laboratory column biostimulation experiments that were poised to precipitate iron and zinc sulfides during sulfate reduction. We develop a petrophysical model based on sphere-shaped cells to link the sulfide precipitate properties to the time-lapse geophysical attributes and estimate volume fraction of the sulfide precipitates, fraction of the dispersed, sulfide-encrusted cells, mean radius of the aggregated clusters, and permeability over the course of the experiments. Results of the case study suggest that the developed state-space approach permits the use of geophysical datasets for providing quantitative estimates of end-product characteristics and hydrological feedbacks associated with biogeochemical transformations. Although tested here on laboratory column experiment datasets, the developed framework provides the foundation needed for quantitative field-scale estimation of biogeochemical parameters over space and time using direct, but often sparse wellbore data with indirect, but more spatially extensive geophysical datasets.

MEDUSA-2.0: an intermediate complexity biogeochemical model of the marine carbon cycle for climate change and ocean acidification studies

Directory of Open Access Journals (Sweden)

A. Yool

2013-10-01

Full Text Available MEDUSA-1.0 (Model of Ecosystem Dynamics, nutrient Utilisation, Sequestration and Acidification was developed as an "intermediate complexity" plankton ecosystem model to study the biogeochemical response, and especially that of the so-called "biological pump", to anthropogenically driven change in the World Ocean (Yool et al., 2011. The base currency in this model was nitrogen from which fluxes of organic carbon, including export to the deep ocean, were calculated by invoking fixed C:N ratios in phytoplankton, zooplankton and detritus. However, due to anthropogenic activity, the atmospheric concentration of carbon dioxide (CO2 has significantly increased above its natural, inter-glacial background. As such, simulating and predicting the carbon cycle in the ocean in its entirety, including ventilation of CO2 with the atmosphere and the resulting impact of ocean acidification on marine ecosystems, requires that both organic and inorganic carbon be afforded a more complete representation in the model specification. Here, we introduce MEDUSA-2.0, an expanded successor model which includes additional state variables for dissolved inorganic carbon, alkalinity, dissolved oxygen and detritus carbon (permitting variable C:N in exported organic matter, as well as a simple benthic formulation and extended parameterizations of phytoplankton growth, calcification and detritus remineralisation. A full description of MEDUSA-2.0, including its additional functionality, is provided and a multi-decadal spin-up simulation (1860–2005 is performed. The biogeochemical performance of the model is evaluated using a diverse range of observational data, and MEDUSA-2.0 is assessed relative to comparable models using output from the Coupled Model Intercomparison Project (CMIP5.

A Process-Based Transport-Distance Model of Aeolian Transport

Science.gov (United States)

Naylor, A. K.; Okin, G.; Wainwright, J.; Parsons, A. J.

2017-12-01

We present a new approach to modeling aeolian transport based on transport distance. Particle fluxes are based on statistical probabilities of particle detachment and distributions of transport lengths, which are functions of particle size classes. A computational saltation model is used to simulate transport distances over a variety of sizes. These are fit to an exponential distribution, which has the advantages of computational economy, concordance with current field measurements, and a meaningful relationship to theoretical assumptions about mean and median particle transport distance. This novel approach includes particle-particle interactions, which are important for sustaining aeolian transport and dust emission. Results from this model are compared with results from both bulk- and particle-sized-specific transport equations as well as empirical wind tunnel studies. The transport-distance approach has been successfully used for hydraulic processes, and extending this methodology from hydraulic to aeolian transport opens up the possibility of modeling joint transport by wind and water using consistent physics. Particularly in nutrient-limited environments, modeling the joint action of aeolian and hydraulic transport is essential for understanding the spatial distribution of biomass across landscapes and how it responds to climatic variability and change.

Biogeochemical controls of uranium bioavailability from the dissolved phase in natural freshwaters

Science.gov (United States)

Croteau, Marie-Noele; Fuller, Christopher C.; Cain, Daniel J.; Campbell, Kate M.; Aiken, George R.

2016-01-01

To gain insights into the risks associated with uranium (U) mining and processing, we investigated the biogeochemical controls of U bioavailability in the model freshwater speciesLymnaea stagnalis (Gastropoda). Bioavailability of dissolved U(VI) was characterized in controlled laboratory experiments over a range of water hardness, pH, and in the presence of complexing ligands in the form of dissolved natural organic matter (DOM). Results show that dissolved U is bioavailable under all the geochemical conditions tested. Uranium uptake rates follow first order kinetics over a range encompassing most environmental concentrations. Uranium uptake rates in L. stagnalis ultimately demonstrate saturation uptake kinetics when exposure concentrations exceed 100 nM, suggesting uptake via a finite number of carriers or ion channels. The lack of a relationship between U uptake rate constants and Ca uptake rates suggest that U does not exclusively use Ca membrane transporters. In general, U bioavailability decreases with increasing pH, increasing Ca and Mg concentrations, and when DOM is present. Competing ions did not affect U uptake rates. Speciation modeling that includes formation constants for U ternary complexes reveals that the aqueous concentration of dicarbonato U species (UO2(CO3)2–2) best predicts U bioavailability to L. stagnalis, challenging the free-ion activity model postulate.

Modeling the biogeochemical impact of atmospheric phosphate deposition from desert dust and combustion sources to the Mediterranean Sea

Science.gov (United States)

Richon, Camille; Dutay, Jean-Claude; Dulac, François; Wang, Rong; Balkanski, Yves

2018-04-01

Daily modeled fields of phosphate deposition to the Mediterranean from natural dust, anthropogenic combustion and wildfires were used to assess the effect of this external nutrient on marine biogeochemistry. The ocean model used is a high-resolution (1/12°) regional coupled dynamical-biogeochemical model of the Mediterranean Sea (NEMO-MED12/PISCES). The input fields of phosphorus are for 2005, which are the only available daily resolved deposition fields from the global atmospheric chemical transport model LMDz-INCA. Traditionally, dust has been suggested to be the main atmospheric source of phosphorus, but the LMDz-INCA model suggests that combustion is dominant over natural dust as an atmospheric source of phosphate (PO4, the bioavailable form of phosphorus in seawater) for the Mediterranean Sea. According to the atmospheric transport model, phosphate deposition from combustion (Pcomb) brings on average 40.5×10-6 mol PO4 m-2 yr-1 over the entire Mediterranean Sea for the year 2005 and is the primary source over the northern part (e.g., 101×10-6 mol PO4 m-2 yr-1 from combustion deposited in 2005 over the north Adriatic against 12.4×10-6 from dust). Lithogenic dust brings 17.2×10-6 mol PO4 m-2 yr-1 on average over the Mediterranean Sea in 2005 and is the primary source of atmospheric phosphate to the southern Mediterranean Basin in our simulations (e.g., 31.8×10-6 mol PO4 m-2 yr-1 from dust deposited in 2005 on average over the south Ionian basin against 12.4×10-6 from combustion). The evaluation of monthly averaged deposition flux variability of Pdust and Pcomb for the 1997-2012 period indicates that these conclusions may hold true for different years. We examine separately the two atmospheric phosphate sources and their respective flux variability and evaluate their impacts on marine surface biogeochemistry (phosphate concentration, chlorophyll a, primary production). The impacts of the different phosphate deposition sources on the biogeochemistry of the

Biogeochemical Cycling and Environmental Stability of Pu Relevant to Long-Term Stewardship of DOE Sites

Energy Technology Data Exchange (ETDEWEB)

Francis, Arokiasamy J.; Santschi, Peter H.; Honeyman, Bruce D.

2005-06-01

The overall objective of this proposed research is to understand the biogeochemical cycling of Pu in environments of interest to long-term DOE stewardship issues. Central to Pu cycling (transport initiation to immobilization) is the role of microorganisms. The hypothesis underlying this proposal is that microbial activity is the causative agent in initiating the mobilization of Pu in near-surface environments: through the transformation of Pu associated with solid phases, production of extracellular polymeric substances (EPS) carrier phases, and the creation of microenvironments. Also, microbial processes are central to the immobilization of Pu species, through the metabolism of organically complexed Pu species and Pu associated with extracellular carrier phases and the creation of environments favorable for Pu transport retardation.

Biogeochemical Cycling and Environmental Stability of Pu Relevant to Long-Term Stewardship of DOE Sites

Energy Technology Data Exchange (ETDEWEB)

Honeyman, Bruce D.

2006-06-01

The overall objective of this proposed research is to understand the biogeochemical cycling of Pu in environments of interest to long-term DOE stewardship issues. Central to Pu cycling (transport initiation to immobilization) is the role of microorganisms. The hypothesis underlying this proposal is that microbial activity is the causative agent in initiating the mobilization of Pu in near-surface environments: through the transformation of Pu associated with solid phases, production of extracellular polymeric substances (EPS) carrier phases, and the creation of microenvironments. Also, microbial processes are central to the immobilization of Pu species, through the metabolism of organically complexed Pu species and Pu associated with extracellular carrier phases and the creation of environments favorable for Pu transport retardation.

Groundwater Impacts of Radioactive Wastes and Associated Environmental Modeling Assessment

Energy Technology Data Exchange (ETDEWEB)

Ma, Rui; Zheng, Chunmiao; Liu, Chongxuan

2012-11-01

This article provides a review of the major sources of radioactive wastes and their impacts on groundwater contamination. The review discusses the major biogeochemical processes that control the transport and fate of radionuclide contaminants in groundwater, and describe the evolution of mathematical models designed to simulate and assess the transport and transformation of radionuclides in groundwater.

Transport of nutrients from land to sea: Global modeling approaches and uncertainty analyses

NARCIS (Netherlands)

Beusen, A.H.W.

2014-01-01

This thesis presents four examples of global models developed as part of the Integrated Model to Assess the Global Environment (IMAGE). They describe different components of global biogeochemical cycles of the nutrients nitrogen (N), phosphorus (P) and silicon (Si), with a focus on approaches to

Application of a hybrid multiscale approach to simulate hydrologic and biogeochemical processes in the river-groundwater interaction zone.

Energy Technology Data Exchange (ETDEWEB)

Hammond, Glenn Edward; Yang, Xiaofan; Song, Xuehang; Song, Hyun-Seob; Hou, Zhangshuan; Chen, Xingyuan; Liu, Yuanyuan; Scheibe, Tim

2017-03-01

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 scales of interest.

Volume reduction outweighs biogeochemical processes in controlling phosphorus treatment in aged detention systems

Science.gov (United States)

Shukla, Asmita; Shukla, Sanjay; Annable, Michael D.; Hodges, Alan W.

2017-08-01

Stormwater detention areas (SDAs) play an important role in treating end-of-the-farm runoff in phosphorous (P) limited agroecosystems. Phosphorus transport from the SDAs, including those through subsurface pathways, are not well understood. The prevailing understanding of these systems assumes that biogeochemical processes play the primary treatment role and that subsurface losses can be neglected. Water and P fluxes from a SDA located in a row-crop farm were measured for two years (2009-2011) to assess the SDA's role in reducing downstream P loads. The SDA treated 55% (497 kg) and 95% (205 kg) of the incoming load during Year 1 (Y1, 09-10) and Year 2 (Y2, 10-11), respectively. These treatment efficiencies were similar to surface water volumetric retention (49% in Y1 and 84% in Y2) and varied primarily with rainfall. Similar water volume and P retentions indicate that volume retention is the main process controlling P loads. A limited role of biogeochemical processes was supported by low to no remaining soil P adsorption capacity due to long-term drainage P input. The fact that outflow P concentrations (Y1 = 368.3 μg L- 1, Y2 = 230.4 μg L- 1) could be approximated by using a simple mixing of rainfall and drainage P input further confirmed the near inert biogeochemical processes. Subsurface P losses through groundwater were 304 kg (27% of inflow P) indicating that they are an important source for downstream P. Including subsurface P losses reduces the treatment efficiency to 35% (from 61%). The aboveground biomass in the SDA contained 42% (240 kg) of the average incoming P load suggesting that biomass harvesting could be a cost-effective alternative for reviving the role of biogeochemical processes to enhance P treatment in aged, P-saturated SDAs. The 20-year present economic value of P removal through harvesting was estimated to be 341,000, which if covered through a cost share or a payment for P treatment services program could be a positive outcome for both

Mapping rice ecosystem dynamics and greenhouse gas emissions using multiscale imagery and biogeochemical models

Science.gov (United States)

Salas, W.; Torbick, N.

2017-12-01

Rice greenhouse gas (GHG) emissions in production hot spots have been mapped using multiscale satellite imagery and a processed-based biogeochemical model. The multiscale Synthetic Aperture Radar (SAR) and optical imagery were co-processed and fed into a machine leanring framework to map paddy attributes that are tuned using field observations and surveys. Geospatial maps of rice extent, crop calendar, hydroperiod, and cropping intensity were then used to parameterize the DeNitrification-DeComposition (DNDC) model to estimate emissions. Results, in the Red River Detla for example, show total methane emissions at 345.4 million kgCH4-C equivalent to 11.5 million tonnes CO2e (carbon dioxide equivalent). We further assessed the role of Alternative Wetting and Drying and the impact on GHG and yield across production hot spots with uncertainty estimates. The approach described in this research provides a framework for using SAR to derive maps of rice and landscape characteristics to drive process models like DNDC. These types of tools and approaches will support the next generation of Monitoring, Reporting, and Verification (MRV) to combat climate change and support ecosystem service markets.

Reply to 'Comment on kinetic modeling of microbially-driven redox chemistry of subsurface environments: coupling transport, microbial metabolism and geochemistry' by J. Griffioen

Science.gov (United States)

Hunter, K. S.; Van Cappellen, P.

2000-01-01

Our paper, 'Kinetic modeling of microbially-driven redox chemistry of subsurface environments: coupling transport, microbial metabolism and geochemistry' (Hunter et al., 1998), presents a theoretical exploration of biogeochemical reaction networks and their importance to the biogeochemistry of groundwater systems. As with any other model, the kinetic reaction-transport model developed in our paper includes only a subset of all physically, biologically and chemically relevant processes in subsurface environments. It considers aquifer systems where the primary energy source driving microbial activity is the degradation of organic matter. In addition to the primary biodegradation pathways of organic matter (i.e. respiration and fermentation), the redox chemistry of groundwaters is also affected by reactions not directly involving organic matter oxidation. We refer to the latter as secondary reactions. By including secondary redox reactions which consume reduced reaction products (e.g., Mn2+, FeS, H2S), and in the process compete with microbial heterotrophic populations for available oxidants (i.e. O2, NO3-, Mn(IV), Fe(III), SO42-), we predict spatio-temporal distributions of microbial activity which differ significantly from those of models which consider only the biodegradation reactions. That is, the secondary reactions have a significant impact on the distributions of the rates of heterotrophic and chemolithotrophic metabolic pathways. We further show that secondary redox reactions, as well as non-redox reactions, significantly influence the acid-base chemistry of groundwaters. The distributions of dissolved inorganic redox species along flowpaths, however, are similar in simulations with and without secondary reactions (see Figs. 3(b) and 7(b) in Hunter et al., 1998), indicating that very different biogeochemical reaction dynamics may lead to essentially the same chemical redox zonation of a groundwater system.

Towards an assessment of riverine dissolved organic carbon in surface waters of the western Arctic Ocean based on remote sensing and biogeochemical modeling

Science.gov (United States)

Le Fouest, Vincent; Matsuoka, Atsushi; Manizza, Manfredi; Shernetsky, Mona; Tremblay, Bruno; Babin, Marcel

2018-03-01

Future climate warming of the Arctic could potentially enhance the load of terrigenous dissolved organic carbon (tDOC) of Arctic rivers due to increased carbon mobilization within watersheds. A greater flux of tDOC might impact the biogeochemical processes of the coastal Arctic Ocean (AO) and ultimately its capacity to absorb atmospheric CO2. In this study, we show that sea-surface tDOC concentrations simulated by a physical-biogeochemical coupled model in the Canadian Beaufort Sea for 2003-2011 compare favorably with estimates retrieved by satellite imagery. Our results suggest that, over spring-summer, tDOC of riverine origin contributes to 35 % of primary production and that an equivalent of ˜ 10 % of tDOC is exported westwards with the potential of fueling the biological production of the eastern Alaskan nearshore waters. The combination of model and satellite data provides promising results to extend this work to the entire AO so as to quantify, in conjunction with in situ data, the expected changes in tDOC fluxes and their potential impact on the AO biogeochemistry at basin scale.

Nitrous oxide emissions from cropland: a procedure for calibrating the DayCent biogeochemical model using inverse modelling

Science.gov (United States)

Rafique, Rashad; Fienen, Michael N.; Parkin, Timothy B.; Anex, Robert P.

2013-01-01

DayCent is a biogeochemical model of intermediate complexity widely used to simulate greenhouse gases (GHG), soil organic carbon and nutrients in crop, grassland, forest and savannah ecosystems. Although this model has been applied to a wide range of ecosystems, it is still typically parameterized through a traditional “trial and error” approach and has not been calibrated using statistical inverse modelling (i.e. algorithmic parameter estimation). The aim of this study is to establish and demonstrate a procedure for calibration of DayCent to improve estimation of GHG emissions. We coupled DayCent with the parameter estimation (PEST) software for inverse modelling. The PEST software can be used for calibration through regularized inversion as well as model sensitivity and uncertainty analysis. The DayCent model was analysed and calibrated using N2O flux data collected over 2 years at the Iowa State University Agronomy and Agricultural Engineering Research Farms, Boone, IA. Crop year 2003 data were used for model calibration and 2004 data were used for validation. The optimization of DayCent model parameters using PEST significantly reduced model residuals relative to the default DayCent parameter values. Parameter estimation improved the model performance by reducing the sum of weighted squared residual difference between measured and modelled outputs by up to 67 %. For the calibration period, simulation with the default model parameter values underestimated mean daily N2O flux by 98 %. After parameter estimation, the model underestimated the mean daily fluxes by 35 %. During the validation period, the calibrated model reduced sum of weighted squared residuals by 20 % relative to the default simulation. Sensitivity analysis performed provides important insights into the model structure providing guidance for model improvement.

Role of zooplankton dynamics for Southern Ocean phytoplankton biomass and global biogeochemical cycles

DEFF Research Database (Denmark)

Le Quéré, Corinne; Buitenhuis, Erik T.; Moriarty, Róisín

2016-01-01

zooplankton community, despite iron limitation of phytoplankton community growth rates. This result has implications for the representation of global biogeochemical cycles in models as zooplankton faecal pellets sink rapidly and partly control the carbon export to the intermediate and deep ocean....

A positive and multi-element conserving time stepping scheme for biogeochemical processes in marine ecosystem models

Science.gov (United States)

Radtke, H.; Burchard, H.

2015-01-01

In this paper, an unconditionally positive and multi-element conserving time stepping scheme for systems of non-linearly coupled ODE's is presented. These systems of ODE's are used to describe biogeochemical transformation processes in marine ecosystem models. The numerical scheme is a positive-definite modification of the Runge-Kutta method, it can have arbitrarily high order of accuracy and does not require time step adaption. If the scheme is combined with a modified Patankar-Runge-Kutta method from Burchard et al. (2003), it also gets the ability to solve a certain class of stiff numerical problems, but the accuracy is restricted to second-order then. The performance of the new scheme on two test case problems is shown.

How do persistent organic pollutants be coupled with biogeochemical cycles of carbon and nutrients in terrestrial ecosystems under global climate change?

Energy Technology Data Exchange (ETDEWEB)

Teng, Ying [Chinese Academy of Sciences, Nanjing (China). Key Lab. of Soil Environment and Pollution Remediation; Griffith Univ., Nathan, QLD (Australia). Environmetnal Futures Centre and School of Biomolecular and Physical Sciences; Xu, Zhihong; Reverchon, Frederique [Griffith Univ., Nathan, QLD (Australia). Environmetnal Futures Centre and School of Biomolecular and Physical Sciences; Luo, Yongming [Chinese Academy of Sciences, Nanjing (China). Key Lab. of Soil Environment and Pollution Remediation

2012-03-15

Global climate change (GCC), especially global warming, has affected the material cycling (e.g., carbon, nutrients, and organic chemicals) and the energy flows of terrestrial ecosystems. Persistent organic pollutants (POPs) were regarded as anthropogenic organic carbon (OC) source, and be coupled with the natural carbon (C) and nutrient biogeochemical cycling in ecosystems. The objective of this work was to review the current literature and explore potential coupling processes and mechanisms between POPs and biogeochemical cycles of C and nutrients in terrestrial ecosystems induced by global warming. Global warming has caused many physical, chemical, and biological changes in terrestrial ecosystems. POPs environmental fate in these ecosystems is controlled mainly by temperature and biogeochemical processes. Global warming may accelerate the re-emissions and redistribution of POPs among environmental compartments via soil-air exchange. Soil-air exchange is a key process controlling the fate and transportation of POPs and terrestrial ecosystem C at regional and global scales. Soil respiration is one of the largest terrestrial C flux induced by microbe and plant metabolism, which can affect POPs biotransformation in terrestrial ecosystems. Carbon flow through food web structure also may have important consequences for the biomagnification of POPs in the ecosystems and further lead to biodiversity loss induced by climate change and POPs pollution stress. Moreover, the integrated techniques and biological adaptation strategy help to fully explore the coupling mechanisms, functioning and trends of POPs and C and nutrient biogeochemical cycling processes in terrestrial ecosystems. There is increasing evidence that the environmental fate of POPs has been linked with biogeochemical cycles of C and nutrients in terrestrial ecosystems under GCC. However, the relationships between POPs and the biogeochemical cycles of C and nutrients are still not well understood. Further

Mangrove forests: a potent nexus of coastal biogeochemical cycling

Science.gov (United States)

Barr, J. G.; Fuentes, J. D.; Shoemaker, B.; O'Halloran, T. L.; Lin, G., Sr.; Engel, V. C.

2014-12-01

Mangrove forests cover just 0.1% of the Earth's terrestrial surface, yet they provide a disproportionate source (~10 % globally) of terrestrially derived, refractory dissolved organic carbon to the oceans. Mangrove forests are biogeochemical reactors that convert biomass into dissolved organic and inorganic carbon at unusually high rates, and many studies recognize the value of mangrove ecosystems for the substantial amounts of soil carbon storage they produce. However, questions remain as to how mangrove forest ecosystem services should be valuated and quantified. Therefore, this study addresses several objectives. First, we demonstrate that seasonal and annual net ecosystem carbon exchange in three selected mangrove forests, derived from long-term eddy covariance measurements, represent key quantities in defining the magnitude of biogeochemical cycling and together with other information on carbon cycle parameters serves as a proxy to estimate ecosystem services. Second, we model ecosystem productivity across the mangrove forests of Everglades National Park and southern China by relating net ecosystem exchange values to remote sensing data. Finally, we develop a carbon budget for the mangrove forests in the Everglades National Park for the purposes of demonstrating that these forests and adjacent estuaries are sites of intense biogeochemical cycling. One conclusion from this study is that much of the carbon entering from the atmosphere as net ecosystem exchange (~1000 g C m-2 yr-1) is not retained in the net ecosystem carbon balance. Instead, a substantial fraction of the carbon entering the system as net ecosystem exchange is ultimately exported to the oceans or outgassed as reaction products within the adjacent estuary.

Water, Energy, and Biogeochemical Model (WEBMOD), user’s manual, version 1

Science.gov (United States)

Webb, Richard M.T.; Parkhurst, David L.

2017-02-08

The Water, Energy, and Biogeochemical Model (WEBMOD) uses the framework of the U.S. Geological Survey (USGS) Modular Modeling System to simulate fluxes of water and solutes through watersheds. WEBMOD divides watersheds into model response units (MRU) where fluxes and reactions are simulated for the following eight hillslope reservoir types: canopy; snowpack; ponding on impervious surfaces; O-horizon; two reservoirs in the unsaturated zone, which represent preferential flow and matrix flow; and two reservoirs in the saturated zone, which also represent preferential flow and matrix flow. The reservoir representing ponding on impervious surfaces, currently not functional (2016), will be implemented once the model is applied to urban areas. MRUs discharge to one or more stream reservoirs that flow to the outlet of the watershed. Hydrologic fluxes in the watershed are simulated by modules derived from the USGS Precipitation Runoff Modeling System; the National Weather Service Hydro-17 snow model; and a topography-driven hydrologic model (TOPMODEL). Modifications to the standard TOPMODEL include the addition of heterogeneous vertical infiltration rates; irrigation; lateral and vertical preferential flows through the unsaturated zone; pipe flow draining the saturated zone; gains and losses to regional aquifer systems; and the option to simulate baseflow discharge by using an exponential, parabolic, or linear decrease in transmissivity. PHREEQC, an aqueous geochemical model, is incorporated to simulate chemical reactions as waters evaporate, mix, and react within the various reservoirs of the model. The reactions that can be specified for a reservoir include equilibrium reactions among water; minerals; surfaces; exchangers; and kinetic reactions such as kinetic mineral dissolution or precipitation, biologically mediated reactions, and radioactive decay. WEBMOD also simulates variations in the concentrations of the stable isotopes deuterium and oxygen-18 as a result of

Inter-annual variability of the atmospheric carbon dioxide concentrations as simulated with global terrestrial biosphere models and an atmospheric transport model

Energy Technology Data Exchange (ETDEWEB)

Fujita, Daisuke; Saeki, Tazu; Nakazawa, Takakiyo [Tohoku Univ., Sendai (Japan). Center for Atmospheric and Oceanic Studies; Ishizawa, Misa; Maksyutov, Shamil [Inst. for Global Change Research, Yokohama (Japan). Frontier Research System for Global Change; Thornton, Peter E. [National Center for Atmospheric Research, Boulder, CO (United States). Climate and Global Dynamics Div.

2003-04-01

Seasonal and inter-annual variations of atmospheric CO{sub 2} for the period from 1961 to 1997 have been simulated using a global tracer transport model driven by a new version of the Biome BioGeochemical Cycle model (Biome-BGC). Biome-BGC was forced by daily temperature and precipitation from the NCEP reanalysis dataset, and the calculated monthly-averaged CO{sub 2} fluxes were used as input to the global transport model. Results from an inter-comparison with the Carnegie-Ames-Stanford Approach model (CASA) and the Simulation model of Carbon CYCLE in Land Ecosystems (Sim-CYCLE) model are also reported. The phase of the seasonal cycle in the Northern Hemisphere was reproduced generally well by Biome-BGC, although the amplitude was smaller compared to the observations and to the other biosphere models. The CO{sub 2} time series simulated by Biome-BGC were compared to the global CO{sub 2} concentration anomalies from the observations at Mauna Loa and the South Pole. The modeled concentration anomalies matched the phase of the inter-annual variations in the atmospheric CO{sub 2} observations; however, the modeled amplitude was lower than the observed value in several cases. The result suggests that a significant part of the inter-annual variability in the global carbon cycle can be accounted for by the terrestrial biosphere models. Simulations performed with another climate-based model, Sim-CYCLE, produced a larger amplitude of inter-annual variability in atmospheric CO{sub 2}, making the amplitude closer to the observed range, but with a more visible phase mismatch in a number of time periods. This may indicate the need to increase the Biome-BGC model sensitivity to seasonal and inter-annual changes in temperature and precipitation.

Inter-annual variability of the atmospheric carbon dioxide concentrations as simulated with global terrestrial biosphere models and an atmospheric transport model

International Nuclear Information System (INIS)

Fujita, Daisuke; Saeki, Tazu; Nakazawa, Takakiyo; Ishizawa, Misa; Maksyutov, Shamil; Thornton, Peter E.

2003-01-01

Seasonal and inter-annual variations of atmospheric CO 2 for the period from 1961 to 1997 have been simulated using a global tracer transport model driven by a new version of the Biome BioGeochemical Cycle model (Biome-BGC). Biome-BGC was forced by daily temperature and precipitation from the NCEP reanalysis dataset, and the calculated monthly-averaged CO 2 fluxes were used as input to the global transport model. Results from an inter-comparison with the Carnegie-Ames-Stanford Approach model (CASA) and the Simulation model of Carbon CYCLE in Land Ecosystems (Sim-CYCLE) model are also reported. The phase of the seasonal cycle in the Northern Hemisphere was reproduced generally well by Biome-BGC, although the amplitude was smaller compared to the observations and to the other biosphere models. The CO 2 time series simulated by Biome-BGC were compared to the global CO 2 concentration anomalies from the observations at Mauna Loa and the South Pole. The modeled concentration anomalies matched the phase of the inter-annual variations in the atmospheric CO 2 observations; however, the modeled amplitude was lower than the observed value in several cases. The result suggests that a significant part of the inter-annual variability in the global carbon cycle can be accounted for by the terrestrial biosphere models. Simulations performed with another climate-based model, Sim-CYCLE, produced a larger amplitude of inter-annual variability in atmospheric CO 2 , making the amplitude closer to the observed range, but with a more visible phase mismatch in a number of time periods. This may indicate the need to increase the Biome-BGC model sensitivity to seasonal and inter-annual changes in temperature and precipitation

Greenland's glacial fjords and their role in regional biogeochemical dynamics.

Science.gov (United States)

Crosby, J.; Arndt, S.

2017-12-01

Greenland's coastal fjords serve as important pathways that connect the Greenland Ice Sheet (GrIS) and the surrounding oceans. They export seasonal glacial meltwater whilst being significant sites of primary production. These fjords are home to some of the most productive ecosystems in the world and possess high socio-economic value via fisheries. A growing number of studies have proposed the GrIS as an underappreciated yet significant source of nutrients to surrounding oceans. Acting as both transfer routes and sinks for glacial nutrient export, fjords have the potential to act as significant biogeochemical processors, yet remain underexplored. Critically, an understanding of the quantitative contribution of fjords to carbon and nutrient budgets is lacking, with large uncertainties associated with limited availability of field data and the lack of robust upscaling approaches. To close this knowledge gap we developed a coupled 2D physical-biogeochemical model of the Godthåbsfjord system, a sub-Arctic sill fjord in southwest Greenland, to quantitatively assess the impact of nutrients exported from the GrIS on fjord primary productivity and biogeochemical dynamics. Glacial meltwater is found to be a key driver of fjord-scale circulation patterns, whilst tracer simulations reveal the relative nutrient contributions from meltwater-driven upwelling and meltwater export from the GrIS. Hydrodynamic circulation patterns and freshwater transit times are explored to provide a first understanding of the glacier-fjord-ocean continuum, demonstrating the complex pattern of carbon and nutrient cycling at this critical land-ocean interface.

Biogeochemical impact of a model western iron source in the Pacific Equatorial Undercurrent

OpenAIRE

Slemons, L.; Gorgues, T.; Aumont, Olivier; Menkès, Christophe; Murray, J. W.

2009-01-01

Trace element distributions in the source waters of the Pacific Equatorial Undercurrent (EUC) show the existence of elevated total acid-soluble iron concentrations. This region has been suggested to contribute enough bioavailable iron to regulate interannual and interglacial variability in biological productivity downstream in the high-nitrate low-chlorophyll upwelling zone of the eastern equatorial Pacific. We investigated the advection and first-order biogeochemical impact of an imposed, da...

Probabilistic transport models for fusion

International Nuclear Information System (INIS)

Milligen, B.Ph. van; Carreras, B.A.; Lynch, V.E.; Sanchez, R.

2005-01-01

A generalization of diffusive (Fickian) transport is considered, in which particle motion is described by probability distributions. We design a simple model that includes a critical mechanism to switch between two transport channels, and show that it exhibits various interesting characteristics, suggesting that the ideas of probabilistic transport might provide a framework for the description of a range of unusual transport phenomena observed in fusion plasmas. The model produces power degradation and profile consistency, as well as a scaling of the confinement time with system size reminiscent of the gyro-Bohm/Bohm scalings observed in fusion plasmas, and rapid propagation of disturbances. In the present work we show how this model may also produce on-axis peaking of the profiles with off-axis fuelling. It is important to note that the fluid limit of a simple model like this, characterized by two transport channels, does not correspond to the usual (Fickian) transport models commonly used for modelling transport in fusion plasmas, and behaves in a fundamentally different way. (author)

Regional impacts of iron-light colimitation in a global biogeochemical model

Science.gov (United States)

Galbraith, E. D.; Gnanadesikan, A.; Dunne, J. P.; Hiscock, M. R.

2010-03-01

Laboratory and field studies have revealed that iron has multiple roles in phytoplankton physiology, with particular importance for light-harvesting cellular machinery. However, although iron-limitation is explicitly included in numerous biogeochemical/ecosystem models, its implementation varies, and its effect on the efficiency of light harvesting is often ignored. Given the complexity of the ocean environment, it is difficult to predict the consequences of applying different iron limitation schemes. Here we explore the interaction of iron and nutrient cycles in an ocean general circulation model using a new, streamlined model of ocean biogeochemistry. Building on previously published parameterizations of photoadaptation and export production, the Biogeochemistry with Light Iron Nutrients and Gasses (BLING) model is constructed with only four explicit tracers but including macronutrient and micronutrient limitation, light limitation, and an implicit treatment of community structure. The structural simplicity of this computationally-inexpensive model allows us to clearly isolate the global effect that iron availability has on maximum light-saturated photosynthesis rates vs. the effect iron has on photosynthetic efficiency. We find that the effect on light-saturated photosynthesis rates is dominant, negating the importance of photosynthetic efficiency in most regions, especially the cold waters of the Southern Ocean. The primary exceptions to this occur in iron-rich regions of the Northern Hemisphere, where high light-saturated photosynthesis rates allow photosynthetic efficiency to play a more important role. In other words, the ability to efficiently harvest photons has little effect in regions where light-saturated growth rates are low. Additionally, we speculate that the phytoplankton cells dominating iron-limited regions tend to have relatively high photosynthetic efficiency, due to reduced packaging effects. If this speculation is correct, it would imply that

Regional impacts of iron-light colimitation in a global biogeochemical model

Directory of Open Access Journals (Sweden)

E. D. Galbraith

2010-03-01

Full Text Available Laboratory and field studies have revealed that iron has multiple roles in phytoplankton physiology, with particular importance for light-harvesting cellular machinery. However, although iron-limitation is explicitly included in numerous biogeochemical/ecosystem models, its implementation varies, and its effect on the efficiency of light harvesting is often ignored. Given the complexity of the ocean environment, it is difficult to predict the consequences of applying different iron limitation schemes. Here we explore the interaction of iron and nutrient cycles in an ocean general circulation model using a new, streamlined model of ocean biogeochemistry. Building on previously published parameterizations of photoadaptation and export production, the Biogeochemistry with Light Iron Nutrients and Gasses (BLING model is constructed with only four explicit tracers but including macronutrient and micronutrient limitation, light limitation, and an implicit treatment of community structure. The structural simplicity of this computationally-inexpensive model allows us to clearly isolate the global effect that iron availability has on maximum light-saturated photosynthesis rates vs. the effect iron has on photosynthetic efficiency. We find that the effect on light-saturated photosynthesis rates is dominant, negating the importance of photosynthetic efficiency in most regions, especially the cold waters of the Southern Ocean. The primary exceptions to this occur in iron-rich regions of the Northern Hemisphere, where high light-saturated photosynthesis rates allow photosynthetic efficiency to play a more important role. In other words, the ability to efficiently harvest photons has little effect in regions where light-saturated growth rates are low. Additionally, we speculate that the phytoplankton cells dominating iron-limited regions tend to have relatively high photosynthetic efficiency, due to reduced packaging effects. If this speculation is correct

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

Science.gov (United States)

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...

Tropical Pacific Climate, Carbon, and Ocean Biogeochemical Response to the Central American Seaway in a GFDL Earth System Model

Science.gov (United States)

Sentman, L. T.; Dunne, J. P.; Stouffer, R. J.; Krasting, J. P.; Wittenberg, A. T.; Toggweiler, J. R.; Broccoli, A. J.

2017-12-01

To explore the tropical Pacific climate, carbon, and ocean biogeochemical response to the shoaling and closure of the Central American Seaway during the Pliocene (5.3-2.6 Ma), we performed a suite of sensitivity experiments using the Geophysical Fluid Dynamics Laboratory Earth System Model, GFDL-ESM2G, varying only the seaway widths and sill depths. These novel ESM simulations include near-final closure of the seaway with a very narrow, 1º grid cell wide opening. Net mass transport through the seaway into the Caribbean is 20.5-23.1 Sv with a deep seaway, but only 14.1 Sv for the wide/shallow seaway because of the inter-basin bi-directional horizontal mass transport. Seaway transport originates from the Antarctic Circumpolar Current in the Pacific and rejoins it in the South Atlantic, reducing the Indonesian Throughflow and transporting heat and salt southward into the South Atlantic, in contrast to present-day and previous seaway simulations. Tropical Pacific mean climate and interannual variability is sensitive to the seaway shoaling, with the largest response to the wider/deeper seaway. In the tropical Pacific, the top 300-m warms 0.4-0.8°C, the equatorial east-west sea surface temperature gradient increases, the north-south sea surface temperature asymmetry at 110°W decreases, thermocline deepens 5-11 m, and the east-west thermocline gradient increases. In the Niño-3 region, ENSO amplitude increases, skewed toward more cold (La Niña) events, El Niño and La Niña develops earlier ( 3 months), the annual cycle weakens and the semi-annual and interannual cycles strengthen from increased symmetry of the north-south sea surface temperature gradient, and atmospheric global teleconnections strengthen with the seaway. The increase in global ocean overturning with the seaway results in a younger average ocean ideal age, reduced dissolved inorganic carbon inventory and marine net primary productivity, and altered inter-basin patterns of surface sediment carbonate

Development and Testing of Embedded Gridding Within the Regional Ocean Modeling System: Interaction Between Near-Shore Off-Shore Currents and Material Properties

National Research Council Canada - National Science Library

McWilliams, James

2003-01-01

...) continuing the development of the Regional Ocean Modeling System (ROMS) with respect to its hydrodynamic algorithms, physical transport parameterizations, and range of biogeochemical processes; (2...

PISCES-v2: an ocean biogeochemical model for carbon and ecosystem studies

Directory of Open Access Journals (Sweden)

O. Aumont

2015-08-01

of marine ecosystems (phytoplankton, microzooplankton and mesozooplankton and the biogeochemical cycles of carbon and of the main nutrients (P, N, Fe, and Si. The model is intended to be used for both regional and global configurations at high or low spatial resolutions as well as for short-term (seasonal, interannual and long-term (climate change, paleoceanography analyses. There are 24 prognostic variables (tracers including two phytoplankton compartments (diatoms and nanophytoplankton, two zooplankton size classes (microzooplankton and mesozooplankton and a description of the carbonate chemistry. Formulations in PISCES-v2 are based on a mixed Monod–quota formalism. On the one hand, stoichiometry of C / N / P is fixed and growth rate of phytoplankton is limited by the external availability in N, P and Si. On the other hand, the iron and silicon quotas are variable and the growth rate of phytoplankton is limited by the internal availability in Fe. Various parameterizations can be activated in PISCES-v2, setting, for instance, the complexity of iron chemistry or the description of particulate organic materials. So far, PISCES-v2 has been coupled to the Nucleus for European Modelling of the Ocean (NEMO and Regional Ocean Modeling System (ROMS systems. A full description of PISCES-v2 and of its optional functionalities is provided here. The results of a quasi-steady-state simulation are presented and evaluated against diverse observational and satellite-derived data. Finally, some of the new functionalities of PISCES-v2 are tested in a series of sensitivity experiments.

Biogeochemical cycling in the Taiwan Strait

Digital Repository Service at National Institute of Oceanography (India)

Naik, H.; Chen, C-T.A.

Based on repeat observations made during 2001-2003 along two transects in the Taiwan Strait this study aims at understanding factors controlling primary productivity with an emphasis on biogeochemical cycling of nitrogen, the major bio...

Reconstructing the Nd oceanic cycle using a coupled dynamical – biogeochemical model

Directory of Open Access Journals (Sweden)

T. Arsouze

2009-12-01

Full Text Available The decoupled behaviour observed between Nd isotopic composition (Nd IC, also referred as ε_Nd and Nd concentration cycles has led to the notion of a "Nd paradox". While ε_Nd behaves in a quasi-conservative way in the open ocean, leading to its broad use as a water-mass tracer, Nd concentration displays vertical profiles that increase with depth, together with a deep-water enrichment along the global thermohaline circulation. This non-conservative behaviour is typical of nutrients affected by scavenging in surface waters and remineralisation at depth. In addition, recent studies suggest the only way to reconcile both concentration and Nd IC oceanic budgets, is to invoke a "Boundary Exchange" process (BE, defined as the co-occurrence of transfer of elements from the margin to the sea with removal of elements from the sea by Boundary Scavenging as a source-sink term. However, these studies do not simulate the input/output fluxes of Nd to the ocean, and therefore prevents from crucial information that limits our understanding of Nd decoupling. To investigate this paradox on a global scale, this study uses for the first time a fully prognostic coupled dynamical/biogeochemical model with an explicit representation of Nd sources and sinks to simulate the Nd oceanic cycle. Sources considered include dissolved river fluxes, atmospheric dusts and margin sediment re-dissolution. Sinks are scavenging by settling particles. This model simulates the global features of the Nd oceanic cycle well, and produces a realistic distribution of Nd concentration (correct order of magnitude, increase with depth and along the conveyor belt, 65% of the simulated values fit in the ±10 pmol/kg envelop when compared to the data and isotopic composition (inter-basin gradient, characterization of the main water-masses, more than 70% of the simulated values fit in the ±3 ε_Nd envelop when compared to the data, though a slight overestimation of

Modelling of transport phenomena

International Nuclear Information System (INIS)

Itoh, Kimitaka; Itoh, Sanae; Fukuyama, Atsushi.

1993-09-01

In this review article, we discuss key features of the transport phenomena and theoretical modelling to understand them. Experimental observations have revealed the nature of anomalous transport, i.e., the enhancement of the transport coefficients by the gradients of the plasma profiles, the pinch phenomena, the radial profile of the anomalous transport coefficients, the variation of the transport among the Bohm diffusion, Pseudo-classical confinement, L-mode and variety of improved confinement modes, and the sudden jumps such as L-H transition. Starting from the formalism of the transport matrix, the modelling based on the low frequency instabilities are reviewed. Theoretical results in the range of drift wave frequency are examined. Problems in theories based on the quasilinear and mixing-length estimates lead to the renewal of the turbulence theory, and the physics picture of the self-sustained turbulence is discussed. The theory of transport using the fluid equation of plasma is developed, showing that the new approach is very promising in explaining abovementioned characteristics of anomalous transport in both L-mode and improved confinement plasmas. The interference of the fluxes is the key to construct the physics basis of the bifurcation theory for the L-H transition. The present status of theories on the mechanisms of improved confinement is discussed. Modelling on the nonlocal nature of transport is briefly discussed. Finally, the impact of the anomalous transport on disruptive phenomena is also described. (author) 95 refs

Calibration of a biome-biogeochemical cycles model for modeling the net primary production of teak forests through inverse modeling of remotely sensed data

Science.gov (United States)

Imvitthaya, Chomchid; Honda, Kiyoshi; Lertlum, Surat; Tangtham, Nipon

2011-01-01

In this paper, we present the results of a net primary production (NPP) modeling of teak (Tectona grandis Lin F.), an important species in tropical deciduous forests. The biome-biogeochemical cycles or Biome-BGC model was calibrated to estimate net NPP through the inverse modeling approach. A genetic algorithm (GA) was linked with Biome-BGC to determine the optimal ecophysiological model parameters. The Biome-BGC was calibrated by adjusting the ecophysiological model parameters to fit the simulated LAI to the satellite LAI (SPOT-Vegetation), and the best fitness confirmed the high accuracy of generated ecophysioligical parameter from GA. The modeled NPP, using optimized parameters from GA as input data, was evaluated using daily NPP derived by the MODIS satellite and the annual field data in northern Thailand. The results showed that NPP obtained using the optimized ecophysiological parameters were more accurate than those obtained using default literature parameterization. This improvement occurred mainly because the model's optimized parameters reduced the bias by reducing systematic underestimation in the model. These Biome-BGC results can be effectively applied in teak forests in tropical areas. The study proposes a more effective method of using GA to determine ecophysiological parameters at the site level and represents a first step toward the analysis of the carbon budget of teak plantations at the regional scale.

Factors Influencing Divergent Patterns of Phosphorus Availability in NY and PA Biogeochemical `Hotspots'

Science.gov (United States)

Saia, S. M.; Hofmeister, K.; Regan, J. M.; Buda, A. R.; Carrick, H. J.; Walter, M. T.

2016-12-01

Anthropogenic alteration of the soil phosphorus (P) cycle leads to subsequent water quality issues in agricultural dominated watersheds. In the humid Northeastern United States (NE US), variably saturated areas can generate surface runoff that transports P and stimulates biogeochemical processes; these hydrologically dynamic locations are often called biogeochemical `hotspots'. Many studies have evaluated nitrogen and carbon cycling in biogeochemical hot spots but few have focused on P. We hypothesized seasonally wet parts of the landscape (i.e., hotspots) have smaller biologically available P pools because runoff events frequently carry away nutrients like P. To test this hypothesis, we generated soil wetness index (SWI) maps from soil (SURRGO) and elevation (LiDAR rescaled to 3 m) data and used these maps to direct seasonal soil sampling near Klingerstown, Pennsylvania (PA) and Ithaca, New York (NY). We collected 5cm deep soil samples in PA (bimonthly) and NY (monthly) along soil moisture gradients for a range of land cover types (forest, fallow, and cropped) from May through October. We measured soil moisture in the field and percent organic matter (OM), pH, and three increasingly strong soil P extractions (dilute-salt-extractable P, oxalate-extractable P, and total-extractable P) in the laboratory. Our results indicated a negative relationship between dilute-salt-extractable P concentrations and SWI in PA and no relationship between these same variables in NY. We also found positive relationships between each of the three P extractions in PA but only a positive relationship between oxalate-extractable P and total-extractable P in NY. Our findings in PA support our hypothesis; namely, less biologically available P (i.e. dilute-salt-extractable P) is found in wetter areas of the landscape. However, divergent P availability patterns in NY point to further complexities and confounding variables in our understanding in soil P processes. Further studies will look

Biogeochemical Processes Responsible for the Enhanced Transport of Plutonium Under transient Unsaturated Ground Water Conditions

Energy Technology Data Exchange (ETDEWEB)

Fred J. Molz, III

2010-05-28

To better understand longer-term vadose zone transport in southeastern soils, field lysimeter experiments were conducted at the Savannah River Site (SRS) near Aiken, SC, in the 1980s. Each of the three lysimeters analyzed herein contained a filter paper spiked with different Pu solutions, and they were left exposed to natural environmental conditions (including the growth of annual weed grasses) for 11 years. The resulting Pu activity measurements from each lysimeter core showed anomalous activity distributions below the source, with significant migration of Pu above the source. Such results are not explainable by adsorption phenomena alone. A transient variably saturated flow model with root water uptake was developed and coupled to a soil reactive transport model. Somewhat surprisingly, the fully transient analysis showed results nearly identical to those of a much simpler steady flow analysis performed previously. However, all phenomena studied were unable to produce the upward Pu transport observed in the data. This result suggests another transport mechanism such as Pu uptake by roots and upward transport due to transpiration. Thus, the variably saturated flow and reactive transport model was extended to include uptake and transport of Pu within the root xylem, along with computational methodology and results. In the extended model, flow velocity in the soil was driven by precipitation input along with transpiration and drainage. Water uptake by the roots determined the flow velocity in the root xylem, and this along with uptake of Pu in the transpiration stream drove advection and dispersion of the two Pu species in the xylem. During wet periods with high potential evapotranspiration, maximum flow velocities through the xylem would approached 600 cm/hr, orders of magnitude larger that flow velocities in the soil. Values for parameters and the correct conceptual viewpoint for Pu transport in plant xylem was uncertain. This motivated further experiments devoted

Biogeochemical Processes Responsible for the Enhanced Transport of Plutonium Under transient Unsaturated Ground Water Conditions

International Nuclear Information System (INIS)

Molz, Fred J. III

2010-01-01

To better understand longer-term vadose zone transport in southeastern soils, field lysimeter experiments were conducted at the Savannah River Site (SRS) near Aiken, SC, in the 1980s. Each of the three lysimeters analyzed herein contained a filter paper spiked with different Pu solutions, and they were left exposed to natural environmental conditions (including the growth of annual weed grasses) for 11 years. The resulting Pu activity measurements from each lysimeter core showed anomalous activity distributions below the source, with significant migration of Pu above the source. Such results are not explainable by adsorption phenomena alone. A transient variably saturated flow model with root water uptake was developed and coupled to a soil reactive transport model. Somewhat surprisingly, the fully transient analysis showed results nearly identical to those of a much simpler steady flow analysis performed previously. However, all phenomena studied were unable to produce the upward Pu transport observed in the data. This result suggests another transport mechanism such as Pu uptake by roots and upward transport due to transpiration. Thus, the variably saturated flow and reactive transport model was extended to include uptake and transport of Pu within the root xylem, along with computational methodology and results. In the extended model, flow velocity in the soil was driven by precipitation input along with transpiration and drainage. Water uptake by the roots determined the flow velocity in the root xylem, and this along with uptake of Pu in the transpiration stream drove advection and dispersion of the two Pu species in the xylem. During wet periods with high potential evapotranspiration, maximum flow velocities through the xylem would approached 600 cm/hr, orders of magnitude larger that flow velocities in the soil. Values for parameters and the correct conceptual viewpoint for Pu transport in plant xylem was uncertain. This motivated further experiments devoted

Effects of Solar UV Radiation and Climate Change on Biogeochemical Cycling: Interactions and Feedbacks

Science.gov (United States)

Solar UV radiation, climate and other drivers of global change are undergoing significant changes and models forecast that these changes will continue for the remainder of this century. Here we assess the effects of solar UV radiation on biogeochemical cycles and the interactions...

Numerical simulation of in-situ chemical oxidation (ISCO) and biodegradation of petroleum hydrocarbons using a coupled model for bio-geochemical reactive transport

Science.gov (United States)

Marin, I. S.; Molson, J. W.

2013-05-01

Petroleum hydrocarbons (PHCs) are a major source of groundwater contamination, being a worldwide and well-known problem. Formed by a complex mixture of hundreds of organic compounds (including BTEX - benzene, toluene, ethylbenzene and xylenes), many of which are toxic and persistent in the subsurface and are capable of creating a serious risk to human health. Several remediation technologies can be used to clean-up PHC contamination. In-situ chemical oxidation (ISCO) and intrinsic bioremediation (IBR) are two promising techniques that can be applied in this case. However, the interaction of these processes with the background aquifer geochemistry and the design of an efficient treatment presents a challenge. Here we show the development and application of BIONAPL/Phreeqc, a modeling tool capable of simulating groundwater flow, contaminant transport with coupled biological and geochemical processes in porous or fractured porous media. BIONAPL/Phreeqc is based on the well-tested BIONAPL/3D model, using a powerful finite element simulation engine, capable of simulating non-aqueous phase liquid (NAPL) dissolution, density-dependent advective-dispersive transport, and solving the geochemical and kinetic processes with the library Phreeqc. To validate the model, we compared BIONAPL/Phreeqc with results from the literature for different biodegradation processes and different geometries, with good agreement. We then used the model to simulate the behavior of sodium persulfate (NaS2O8) as an oxidant for BTEX degradation, coupled with sequential biodegradation in a 2D case and to evaluate the effect of inorganic geochemistry reactions. The results show the advantages of a treatment train remediation scheme based on ISCO and IBR. The numerical performance and stability of the integrated BIONAPL/Phreeqc model was also verified.

Preface to: Indian Ocean biogeochemical processes and ecological variability

Digital Repository Service at National Institute of Oceanography (India)

Hood, R.R.; Naqvi, S.W.A.; Wiggert, J.D.

monsoonal in fluence. The biogeochemical and ecological impacts of this complex physical forcing are not yet fully understood. The Indian Ocean is truly one of the last great frontiers of ocea- nographic research. In addition, it appears... to be particularly vulnerable to climate change and anthropogenic impacts, yet it has been more than a decade since the last coordinated international study of biogeochemical and ecological proc esses was undertaken in this region. To obtain a better un...

Stochastic models of intracellular transport

KAUST Repository

Bressloff, Paul C.

2013-01-09

The interior of a living cell is a crowded, heterogenuous, fluctuating environment. Hence, a major challenge in modeling intracellular transport is to analyze stochastic processes within complex environments. Broadly speaking, there are two basic mechanisms for intracellular transport: passive diffusion and motor-driven active transport. Diffusive transport can be formulated in terms of the motion of an overdamped Brownian particle. On the other hand, active transport requires chemical energy, usually in the form of adenosine triphosphate hydrolysis, and can be direction specific, allowing biomolecules to be transported long distances; this is particularly important in neurons due to their complex geometry. In this review a wide range of analytical methods and models of intracellular transport is presented. In the case of diffusive transport, narrow escape problems, diffusion to a small target, confined and single-file diffusion, homogenization theory, and fractional diffusion are considered. In the case of active transport, Brownian ratchets, random walk models, exclusion processes, random intermittent search processes, quasi-steady-state reduction methods, and mean-field approximations are considered. Applications include receptor trafficking, axonal transport, membrane diffusion, nuclear transport, protein-DNA interactions, virus trafficking, and the self-organization of subcellular structures. © 2013 American Physical Society.

Mechanistic site-based emulation of a global ocean biogeochemical model (MEDUSA 1.0 for parametric analysis and calibration: an application of the Marine Model Optimization Testbed (MarMOT 1.1

Directory of Open Access Journals (Sweden)

J. C. P. Hemmings

2015-03-01

Full Text Available Biogeochemical ocean circulation models used to investigate the role of plankton ecosystems in global change rely on adjustable parameters to capture the dominant biogeochemical dynamics of a complex biological system. In principle, optimal parameter values can be estimated by fitting models to observational data, including satellite ocean colour products such as chlorophyll that achieve good spatial and temporal coverage of the surface ocean. However, comprehensive parametric analyses require large ensemble experiments that are computationally infeasible with global 3-D simulations. Site-based simulations provide an efficient alternative but can only be used to make reliable inferences about global model performance if robust quantitative descriptions of their relationships with the corresponding 3-D simulations can be established. The feasibility of establishing such a relationship is investigated for an intermediate complexity biogeochemistry model (MEDUSA coupled with a widely used global ocean model (NEMO. A site-based mechanistic emulator is constructed for surface chlorophyll output from this target model as a function of model parameters. The emulator comprises an array of 1-D simulators and a statistical quantification of the uncertainty in their predictions. The unknown parameter-dependent biogeochemical environment, in terms of initial tracer concentrations and lateral flux information required by the simulators, is a significant source of uncertainty. It is approximated by a mean environment derived from a small ensemble of 3-D simulations representing variability of the target model behaviour over the parameter space of interest. The performance of two alternative uncertainty quantification schemes is examined: a direct method based on comparisons between simulator output and a sample of known target model "truths" and an indirect method that is only partially reliant on knowledge of the target model output. In general, chlorophyll

The Role of External Inputs and Internal Cycling in Shaping the Global Ocean Cobalt Distribution: Insights From the First Cobalt Biogeochemical Model

Science.gov (United States)

Tagliabue, Alessandro; Hawco, Nicholas J.; Bundy, Randelle M.; Landing, William M.; Milne, Angela; Morton, Peter L.; Saito, Mak A.

2018-04-01

Cobalt is an important micronutrient for ocean microbes as it is present in vitamin B12 and is a co-factor in various metalloenzymes that catalyze cellular processes. Moreover, when seawater availability of cobalt is compared to biological demands, cobalt emerges as being depleted in seawater, pointing to a potentially important limiting role. To properly account for the potential biological role for cobalt, there is therefore a need to understand the processes driving the biogeochemical cycling of cobalt and, in particular, the balance between external inputs and internal cycling. To do so, we developed the first cobalt model within a state-of-the-art three-dimensional global ocean biogeochemical model. Overall, our model does a good job in reproducing measurements with a correlation coefficient of >0.7 in the surface and >0.5 at depth. We find that continental margins are the dominant source of cobalt, with a crucial role played by supply under low bottom-water oxygen conditions. The basin-scale distribution of cobalt supplied from margins is facilitated by the activity of manganese-oxidizing bacteria being suppressed under low oxygen and low temperatures, which extends the residence time of cobalt. Overall, we find a residence time of 7 and 250 years in the upper 250 m and global ocean, respectively. Importantly, we find that the dominant internal resupply process switches from regeneration and recycling of particulate cobalt to dissolution of scavenged cobalt between the upper ocean and the ocean interior. Our model highlights key regions of the ocean where biological activity may be most sensitive to cobalt availability.

Coupling ANIMO and MT3DMS for 3D regional-scale modeling of nutrient transport in soil and groundwater

Science.gov (United States)

Janssen, G.; Del Val Alonso, L.; Groenendijk, P.; Griffioen, J.

2012-12-01

We developed an on-line coupling between the 1D/quasi-2D nutrient transport model ANIMO and the 3D groundwater transport model code MT3DMS. ANIMO is a detailed, process-oriented model code for the simulation of nitrate leaching to groundwater, N- and P-loads on surface waters and emissions of greenhouse gasses. It is the leading nutrient fate and transport code in the Netherlands where it is used primarily for the evaluation of fertilization related legislation. In addition, the code is applied frequently in international research projects. MT3DMS is probably the most commonly used groundwater solute transport package worldwide. The on-line model coupling ANIMO-MT3DMS combines the state-of-the-art descriptions of the biogeochemical cycles in ANIMO with the advantages of using a 3D approach for the transport through the saturated domain. These advantages include accounting for regional lateral transport, considering groundwater-surface water interactions more explicitly, and the possibility of using MODFLOW to obtain the flow fields. An additional merit of the on-line coupling concept is that it preserves feedbacks between the saturated and unsaturated zone. We tested ANIMO-MT3DMS by simulating nutrient transport for the period 1970-2007 in a Dutch agricultural polder catchment covering an area of 118 km2. The transient groundwater flow field had a temporal resolution of one day and was calculated with MODFLOW-MetaSWAP. The horizontal resolution of the model grid was 100x100m and consisted of 25 layers of varying thickness. To keep computation times manageable, we prepared MT3DMS for parallel computing, which in itself is a relevant development for a large community of groundwater transport modelers. For the parameterization of the soil, we applied a standard classification approach, representing the area by 60 units with unique combinations of soil type, land use and geohydrological setting. For the geochemical parameterization of the deeper subsurface, however, we

Parameter Sensitivity and Laboratory Benchmarking of a Biogeochemical Process Model for Enhanced Anaerobic Dechlorination

Science.gov (United States)

Kouznetsova, I.; Gerhard, J. I.; Mao, X.; Barry, D. A.; Robinson, C.; Brovelli, A.; Harkness, M.; Fisher, A.; Mack, E. E.; Payne, J. A.; Dworatzek, S.; Roberts, J.

2008-12-01

A detailed model to simulate trichloroethene (TCE) dechlorination in anaerobic groundwater systems has been developed and implemented through PHAST, a robust and flexible geochemical modeling platform. The approach is comprehensive but retains flexibility such that models of varying complexity can be used to simulate TCE biodegradation in the vicinity of nonaqueous phase liquid (NAPL) source zones. The complete model considers a full suite of biological (e.g., dechlorination, fermentation, sulfate and iron reduction, electron donor competition, toxic inhibition, pH inhibition), physical (e.g., flow and mass transfer) and geochemical processes (e.g., pH modulation, gas formation, mineral interactions). Example simulations with the model demonstrated that the feedback between biological, physical, and geochemical processes is critical. Successful simulation of a thirty-two-month column experiment with site soil, complex groundwater chemistry, and exhibiting both anaerobic dechlorination and endogenous respiration, provided confidence in the modeling approach. A comprehensive suite of batch simulations was then conducted to estimate the sensitivity of predicted TCE degradation to the 36 model input parameters. A local sensitivity analysis was first employed to rank the importance of parameters, revealing that 5 parameters consistently dominated model predictions across a range of performance metrics. A global sensitivity analysis was then performed to evaluate the influence of a variety of full parameter data sets available in the literature. The modeling study was performed as part of the SABRE (Source Area BioREmediation) project, a public/private consortium whose charter is to determine if enhanced anaerobic bioremediation can result in effective and quantifiable treatment of chlorinated solvent DNAPL source areas. The modelling conducted has provided valuable insight into the complex interactions between processes in the evolving biogeochemical systems

SATURATED ZONE FLOW AND TRANSPORT MODEL ABSTRACTION

International Nuclear Information System (INIS)

B.W. ARNOLD

2004-01-01

The purpose of the saturated zone (SZ) flow and transport model abstraction task is to provide radionuclide-transport simulation results for use in the total system performance assessment (TSPA) for license application (LA) calculations. This task includes assessment of uncertainty in parameters that pertain to both groundwater flow and radionuclide transport in the models used for this purpose. This model report documents the following: (1) The SZ transport abstraction model, which consists of a set of radionuclide breakthrough curves at the accessible environment for use in the TSPA-LA simulations of radionuclide releases into the biosphere. These radionuclide breakthrough curves contain information on radionuclide-transport times through the SZ. (2) The SZ one-dimensional (I-D) transport model, which is incorporated in the TSPA-LA model to simulate the transport, decay, and ingrowth of radionuclide decay chains in the SZ. (3) The analysis of uncertainty in groundwater-flow and radionuclide-transport input parameters for the SZ transport abstraction model and the SZ 1-D transport model. (4) The analysis of the background concentration of alpha-emitting species in the groundwater of the SZ

Reconciling surface ocean productivity, export fluxes and sediment composition in a global biogeochemical ocean model

Directory of Open Access Journals (Sweden)

M. Gehlen

2006-01-01

Full Text Available This study focuses on an improved representation of the biological soft tissue pump in the global three-dimensional biogeochemical ocean model PISCES. We compare three parameterizations of particle dynamics: (1 the model standard version including two particle size classes, aggregation-disaggregation and prescribed sinking speed; (2 an aggregation-disaggregation model with a particle size spectrum and prognostic sinking speed; (3 a mineral ballast parameterization with no size classes, but prognostic sinking speed. In addition, the model includes a description of surface sediments and organic carbon early diagenesis. Model output is compared to data or data based estimates of ocean productivity, pe-ratios, particle fluxes, surface sediment bulk composition and benthic O2 fluxes. Model results suggest that different processes control POC fluxes at different depths. In the wind mixed layer turbulent particle coagulation appears as key process in controlling pe-ratios. Parameterization (2 yields simulated pe-ratios that compare well to observations. Below the wind mixed layer, POC fluxes are most sensitive to the intensity of zooplankton flux feeding, indicating the importance of zooplankton community composition. All model parameters being kept constant, the capability of the model to reproduce yearly mean POC fluxes below 2000 m and benthic oxygen demand does at first order not dependent on the resolution of the particle size spectrum. Aggregate formation appears essential to initiate an intense biological pump. At great depth the reported close to constant particle fluxes are most likely the result of the combined effect of aggregate formation and mineral ballasting.

Hydrological and associated biogeochemical consequences of rapid global warming during the Paleocene-Eocene Thermal Maximum

Science.gov (United States)

Carmichael, Matthew J.; Inglis, Gordon N.; Badger, Marcus P. S.; Naafs, B. David A.; Behrooz, Leila; Remmelzwaal, Serginio; Monteiro, Fanny M.; Rohrssen, Megan; Farnsworth, Alexander; Buss, Heather L.; Dickson, Alexander J.; Valdes, Paul J.; Lunt, Daniel J.; Pancost, Richard D.

2017-10-01

The Paleocene-Eocene Thermal Maximum (PETM) hyperthermal, 56 million years ago (Ma), is the most dramatic example of abrupt Cenozoic global warming. During the PETM surface temperatures increased between 5 and 9 °C and the onset likely took hydrological and associated biogeochemical feedbacks, and proxy data from the PETM can provide constraints on changes in warm climate hydrology simulated by general circulation models (GCMs). In this paper, we provide a critical review of biological and geochemical signatures interpreted as direct or indirect indicators of hydrological change at the PETM, explore the importance of adopting multi-proxy approaches, and present a preliminary model-data comparison. Hydrological records complement those of temperature and indicate that the climatic response at the PETM was complex, with significant regional and temporal variability. This is further illustrated by the biogeochemical consequences of inferred changes in hydrology and, in fact, changes in precipitation and the biogeochemical consequences are often conflated in geochemical signatures. There is also strong evidence in many regions for changes in the episodic and/or intra-annual distribution of precipitation that has not widely been considered when comparing proxy data to GCM output. Crucially, GCM simulations indicate that the response of the hydrological cycle to the PETM was heterogeneous - some regions are associated with increased precipitation - evaporation (P - E), whilst others are characterised by a decrease. Interestingly, the majority of proxy data come from the regions where GCMs predict an increase in PETM precipitation. We propose that comparison of hydrological proxies to GCM output can be an important test of model skill, but this will be enhanced by further data from regions of model-simulated aridity and simulation of extreme precipitation events.

Glacial-interglacial variability in ocean oxygen and phosphorus in a global biogeochemical model

Directory of Open Access Journals (Sweden)

V Palastanga

2013-02-01

Full Text Available Increased transfer of particulate matter from continental shelves to the open ocean during glacials may have had a major impact on the biogeochemistry of the ocean. Here, we assess the response of the coupled oceanic cycles of oxygen, carbon, phosphorus, and iron to the input of particulate organic carbon and reactive phosphorus from shelves. We use a biogeochemical ocean model and specifically focus on the Last Glacial Maximum (LGM. When compared to an interglacial reference run, our glacial scenario with shelf input shows major increases in ocean productivity and phosphorus burial, while mean deep-water oxygen concentrations decline. There is a downward expansion of the oxygen minimum zones (OMZs in the Atlantic and Indian Ocean, while the extension of the OMZ in the Pacific is slightly reduced. Oxygen concentrations below 2000 m also decline but bottom waters do not become anoxic. The model simulations show when shelf input of particulate organic matter and particulate reactive P is considered, low oxygen areas in the glacial ocean expand, but concentrations are not low enough to generate wide scale changes in sediment biogeochemistry and sedimentary phosphorus recycling. Increased reactive phosphorus burial in the open ocean during the LGM in the model is related to dust input, notably over the southwest Atlantic and northwest Pacific, whereas input of material from shelves explains higher burial fluxes in continental slope and rise regions. Our model results are in qualitative agreement with available data and reproduce the strong spatial differences in the response of phosphorus burial to glacial-interglacial change. Our model results also highlight the need for additional sediment core records from all ocean basins to allow further insight into changes in phosphorus, carbon and oxygen dynamics in the ocean on glacial-interglacial timescales.

A Eulerian nutrient to fish model of the Baltic Sea — A feasibility-study

Science.gov (United States)

Radtke, Hagen; Neumann, Thomas; Fennel, Wolfgang

2013-09-01

A nutrient-to-fish-model with an explicit two-way interaction between a biogeochemical model of the lower food web and a fish model component is presented for the example of the Baltic Sea, demonstrating the feasibility of a consistent coupling of the upper and lower parts of the food web in a Eulerian model system. In the Baltic Sea, the fish stock is dominated by two prey species (sprat and herring) and one predator (cod). The dynamics of the fish model is driven by size (mass-class) dependent predator-prey interactions while the interaction between the biogeochemical and Fish model component is established through feeding of prey fish on zooplankton and recycling of fish biomass to nutrients and detritus. The fish model component is coupled to an advanced three dimensional biogeochemical model (ERGOM, Neumann et al., 2002). A horizontally explicit representation of fish requires the implementation of fish behavior. As a first step, we propose an algorithm to stimulate fish migration by letting the fish follow the food. Moreover, fish species are guided to their respective spawning areas. Results of first three-dimensional simulations are presented with emphasis on the transport of matter by moving fish. The spawning areas of cod and sprat are in the deep basins, which are not well reached by advective transport. Hence the deposition of matter in these areas by spawning fish could play some role in the distribution of matter. The approach is not limited to applications for the Baltic and the model can be transferred also to other systems.

Radionuclide release from simulated waste material after biogeochemical leaching of uraniferous mineral samples

International Nuclear Information System (INIS)

Williamson, Aimee Lynn; Caron, François; Spiers, Graeme

2014-01-01

Biogeochemical mineral dissolution is a promising method for the released of metals in low-grade host mineralization that contain sulphidic minerals. The application of biogeochemical mineral dissolution to engineered leach heap piles in the Elliot Lake region may be considered as a promising passive technology for the economic recovery of low grade Uranium-bearing ores. In the current investigation, the decrease of radiological activity of uraniferous mineral material after biogeochemical mineral dissolution is quantified by gamma spectroscopy and compared to the results from digestion/ICP-MS analysis of the ore materials to determine if gamma spectroscopy is a simple, viable alternative quantification method for heavy nuclides. The potential release of Uranium (U) and Radium-226 ( 226 Ra) to the aqueous environment from samples that have been treated to represent various stages of leaching and passive closure processes are assessed. Dissolution of U from the solid phase has occurred during biogeochemical mineral dissolution in the presence of Acidithiobacillus ferrooxidans, with gamma spectroscopy indicating an 84% decrease in Uranium-235 ( 235 U) content, a value in accordance with the data obtained by dissolution chemistry. Gamma spectroscopy data indicate that only 30% of the 226 Ra was removed during the biogeochemical mineral dissolution. Chemical inhibition and passivation treatments of waste materials following the biogeochemical mineral dissolution offer greater protection against residual U and 226 Ra leaching. Pacified samples resist the release of 226 Ra contained in the mineral phase and may offer more protection to the aqueous environment for the long term, compared to untreated or inhibited residues, and should be taken into account for future decommissioning. - Highlights: • Gamma counting showed an 84% decrease in 235 U after biogeochemical mineral leaching. • Chemical digestion/ICP-MS analysis also showed an 84% decrease in total U. • Over

Transport properties site descriptive model. Guidelines for evaluation and modelling

International Nuclear Information System (INIS)

Berglund, Sten; Selroos, Jan-Olof

2004-04-01

This report describes a strategy for the development of Transport Properties Site Descriptive Models within the SKB Site Investigation programme. Similar reports have been produced for the other disciplines in the site descriptive modelling (Geology, Hydrogeology, Hydrogeochemistry, Rock mechanics, Thermal properties, and Surface ecosystems). These reports are intended to guide the site descriptive modelling, but also to provide the authorities with an overview of modelling work that will be performed. The site descriptive modelling of transport properties is presented in this report and in the associated 'Strategy for the use of laboratory methods in the site investigations programme for the transport properties of the rock', which describes laboratory measurements and data evaluations. Specifically, the objectives of the present report are to: Present a description that gives an overview of the strategy for developing Site Descriptive Models, and which sets the transport modelling into this general context. Provide a structure for developing Transport Properties Site Descriptive Models that facilitates efficient modelling and comparisons between different sites. Provide guidelines on specific modelling issues where methodological consistency is judged to be of special importance, or where there is no general consensus on the modelling approach. The objectives of the site descriptive modelling process and the resulting Transport Properties Site Descriptive Models are to: Provide transport parameters for Safety Assessment. Describe the geoscientific basis for the transport model, including the qualitative and quantitative data that are of importance for the assessment of uncertainties and confidence in the transport description, and for the understanding of the processes at the sites. Provide transport parameters for use within other discipline-specific programmes. Contribute to the integrated evaluation of the investigated sites. The site descriptive modelling of

Progress in transport modelling of internal transport barrier plasmas in JET

International Nuclear Information System (INIS)

Tala, T.; Bourdelle, C.; Imbeaux, F.; Moreau, D.; Garbet, X.; Joffrin, E.; Laborde, L.; Litaudon, X.; Mazon, D.; Parail, V.; Corrigan, G.; Heading, D.; Crisanti, F.; Mantica, P.; Salmi, A.; Strand, P.; Weiland, J.

2005-01-01

This paper will report on the recent progress in transport modelling of Internal Transport Barrier (ITB) plasmas. Two separate issues will be covered, fully predictive transport modelling of ITBs in the multi-tokamak database, including micro-stability analyses of ITBs, and predictive closed-loop (i.e. real-time control) transport simulations of the q-profile and ITBs. For the first time, the predictive capabilities of the mixed Bohm/GyroBohm and Weiland transport models are investigated with discharges from the ITPA ITB database by fully predictive transport simulations. The predictive transport simulations with the Bohm/GyroBohm model agree very well with experimental results from JET and JT-60U. In order to achieve a good agreement in DIII-D, the stabilisation had to be included into the model, showing the significant role played by the stabilisation in governing the physics of the ITBs. The significant role of the stabilisation is also emphasised by the gyrokinetic analysis. The Weiland transport model shows only limited agreement between the model predictions and experimental results with respect to the formation and location of the ITB. The fully predictive closed-loop simulations with real-time control of the q-profile and ITB show that it is possible to reach various set-point profiles for q and ITB and control them for longer than a current diffusion time in JET using the same real-time control technique as in the experiments. (author)

Transport modelling for ergodic configurations

International Nuclear Information System (INIS)

Runov, A.; Kasilov, S.V.; McTaggart, N.; Schneider, R.; Bonnin, X.; Zagorski, R.; Reiter, D.

2004-01-01

The effect of ergodization, either by additional coils like in TEXTOR-dynamic ergodic divertor (DED) or by intrinsic plasma effects like in W7-X, defines the need for transport models that are able to describe the ergodic configuration properly. A prerequisite for this is the concept of local magnetic coordinates allowing a correct discretization with minimized numerical errors. For these coordinates the appropriate full metric tensor has to be known. To study the transport in complex edge geometries (in particular for W7-X) two possible methods are used. First, a finite-difference discretization of the transport equations on a custom-tailored grid in local magnetic coordinates is used. This grid is generated by field-line tracing to guarantee an exact discretization of the dominant parallel transport (thus also minimizing the numerical diffusion problem). The perpendicular fluxes are then interpolated in a plane (a toroidal cut), where the interpolation problem for a quasi-isotropic system has to be solved by a constrained Delaunay triangulation (keeping the structural information for magnetic surfaces if they exist) and discretization. All toroidal terms are discretized by finite differences. Second, a Monte Carlo transport model originally developed for the modelling of the DED configuration of TEXTOR is used. A generalization and extension of this model was necessary to be able to handle W7-X. The model solves the transport equations with Monte Carlo techniques making use of mappings of local magnetic coordinates. The application of this technique to W7-X in a limiter-like configuration is presented. The decreasing dominance of parallel transport with respect to radial transport for electron heat, ion heat and particle transport results in increasingly steep profiles for the respective quantities within the islands. (author)

Innovative mathematical modeling in environmental remediation

Energy Technology Data Exchange (ETDEWEB)

Yeh, Gour T. [Taiwan Typhoon and Flood Research Institute (Taiwan); National Central Univ. (Taiwan); Univ. of Central Florida (United States); Gwo, Jin Ping [Nuclear Regulatory Commission (NRC), Rockville, MD (United States); Siegel, Malcolm D. [Sandia National Laboratories, Albuquerque, NM (United States); Li, Ming-Hsu [National Central Univ. (Taiwan); ; Fang, Yilin [Pacific Northwest National Laboratory (PNNL), Richland, WA (United States); Zhang, Fan [Inst. of Tibetan Plateau Research, Chinese Academy of Sciences (China); Luo, Wensui [Inst. of Tibetan Plateau Research, Chinese Academy of Sciences (China); Yabusaki, Steven B. [Pacific Northwest National Laboratory (PNNL), Richland, WA (United States)

2013-05-01

There are two different ways to model reactive transport: ad hoc and innovative reaction-based approaches. The former, such as the Kd simplification of adsorption, has been widely employed by practitioners, while the latter has been mainly used in scientific communities for elucidating mechanisms of biogeochemical transport processes. It is believed that innovative mechanistic-based models could serve as protocols for environmental remediation as well. This paper reviews the development of a mechanistically coupled fluid flow, thermal transport, hydrologic transport, and reactive biogeochemical model and example-applications to environmental remediation problems. Theoretical bases are sufficiently described. Four example problems previously carried out are used to demonstrate how numerical experimentation can be used to evaluate the feasibility of different remediation approaches. The first one involved the application of a 56-species uranium tailing problem to the Melton Branch Subwatershed at Oak Ridge National Laboratory (ORNL) using the parallel version of the model. Simulations were made to demonstrate the potential mobilization of uranium and other chelating agents in the proposed waste disposal site. The second problem simulated laboratory-scale system to investigate the role of natural attenuation in potential off-site migration of uranium from uranium mill tailings after restoration. It showed inadequacy of using a single Kd even for a homogeneous medium. The third example simulated laboratory experiments involving extremely high concentrations of uranium, technetium, aluminum, nitrate, and toxic metals (e.g.,Ni, Cr, Co).The fourth example modeled microbially-mediated immobilization of uranium in an unconfined aquifer using acetate amendment in a field-scale experiment. The purposes of these modeling studies were to simulate various mechanisms of mobilization and immobilization of radioactive wastes and to illustrate how to apply reactive transport models

Innovative mathematical modeling in environmental remediation

International Nuclear Information System (INIS)

Yeh, Gour T.; Gwo, Jin Ping; Siegel, Malcolm D.; Li, Ming-Hsu; Fang, Yilin; Zhang, Fan; Luo, Wensui; Yabusaki, Steven B.

2013-01-01

There are two different ways to model reactive transport: ad hoc and innovative reaction-based approaches. The former, such as the Kd simplification of adsorption, has been widely employed by practitioners, while the latter has been mainly used in scientific communities for elucidating mechanisms of biogeochemical transport processes. It is believed that innovative mechanistic-based models could serve as protocols for environmental remediation as well. This paper reviews the development of a mechanistically coupled fluid flow, thermal transport, hydrologic transport, and reactive biogeochemical model and example-applications to environmental remediation problems. Theoretical bases are sufficiently described. Four example problems previously carried out are used to demonstrate how numerical experimentation can be used to evaluate the feasibility of different remediation approaches. The first one involved the application of a 56-species uranium tailing problem to the Melton Branch Subwatershed at Oak Ridge National Laboratory (ORNL) using the parallel version of the model. Simulations were made to demonstrate the potential mobilization of uranium and other chelating agents in the proposed waste disposal site. The second problem simulated laboratory-scale system to investigate the role of natural attenuation in potential off-site migration of uranium from uranium mill tailings after restoration. It showed inadequacy of using a single Kd even for a homogeneous medium. The third example simulated laboratory experiments involving extremely high concentrations of uranium, technetium, aluminum, nitrate, and toxic metals (e.g.,Ni, Cr, Co). The fourth example modeled microbially-mediated immobilization of uranium in an unconfined aquifer using acetate amendment in a field-scale experiment. The purposes of these modeling studies were to simulate various mechanisms of mobilization and immobilization of radioactive wastes and to illustrate how to apply reactive transport models

Structure of peat soils and implications for biogeochemical processes and hydrological flow

Science.gov (United States)

Rezanezhad, F.; McCarter, C. P. R.; Gharedaghloo, B.; Kleimeier, C.; Milojevic, T.; Liu, H.; Weber, T. K. D.; Price, J. S.; Quinton, W. L.; Lenartz, B.; Van Cappellen, P.

2017-12-01

Permafrost peatlands contain globally important amounts of soil organic carbon and play major roles in global water, nutrient and biogeochemical cycles. The structure of peatland soils (i.e., peat) are highly complex with unique physical and hydraulic properties; where significant, and only partially reversible, shrinkage occurs during dewatering (including water table fluctuations), compression and/or decomposition. These distinct physical and hydraulic properties controls water flow, which in turn affect reactive and non-reactive solute transport (such as, sorption or degradation) and biogeochemical functions. Additionally, peat further attenuates solute migration through molecular diffusion into the inactive pores of Sphagnum dominated peat. These slow, diffusion-limited solute exchanges between the pore regions may give rise to pore-scale chemical gradients and heterogeneous distributions of microbial habitats and activity in peat soils. Permafrost peat plateaus have the same essential subsurface characteristics as other widely organic soil-covered peatlands, where the hydraulic conductivity is related to the degree of decomposition and soil compression. Increasing levels of decomposition correspond with a reduction of effective pore diameter and consequently restrict water and solute flow (by several orders of magnitude in hydraulic conductivity between the ground surface and a depth of 50 cm). In this presentation, we present the current knowledge of key physical and hydraulic properties related to the structure of globally available peat soils and discuss their implications for water storage, flow and the migration of solutes.

Diel biogeochemical processes in terrestrial waters

Science.gov (United States)

Nimick, David A.; Gammons, Christopher H.

2011-01-01

scientists typically do not encompass the wide diversity and range of processes that produce diel cycles, and (3) the logistics of making field measurements for 24-h periods has limited recognition and understanding of these important cycles. Thus, the topical session brought together hydrologists, biologists, geochemists, and ecologists to discuss field studies, laboratory experiments, theoretical modeling, and measurement techniques related to diel cycling. Hopefully with the cross-disciplinary synergy developed at the session as well as by this special issue, a more comprehensive understanding of the interrelationships between the diel processes will be developed. Needless to say, understanding diel processes is critical for regulatory agencies and the greater scientific community. And perhaps more importantly, expanded knowledge of biogeochemical cycling may lead to better predictions of how aquatic ecosystems might react to changing conditions of contaminant loading, eutrophication, climate change, drought, industrialization, development, and other variables.

A high-resolution hydrodynamic-biogeochemical coupled model of the Gulf of Cadiz – Alboran Sea region.

Directory of Open Access Journals (Sweden)

D. M. MACIAS

2014-12-01

Full Text Available The southern Iberia regional seas comprise the Gulf of Cadiz and the Alboran Sea sub-basins connected by the narrow Strait of Gibraltar. Both basins are very different in their hydrological and biological characteristics but are, also, tightly connected to each other. Integrative studies of the whole regional oceanic system are scarce and difficult to perform due to the relative large area to cover and the different relevant time-scales of the main forcings in each sub-basin. Here we propose, for the first time, a fully coupled, 3D, hydrodynamic-biogeochemical model that covers, in a single domain (~2km resolution both marine basins for a 20 years simulation (1989-2008. Model performance is assessed against available data in terms of spatial and temporal distributions of biological variables. In general, the proposed model is able to represent the climatological distributions of primary and secondary producers and also the main seasonality of primary production in the different sub-regions of the analyzed basins. Potential causes of the observed mismatches between model and data are identified and some solutions are proposed for future model development. We conclude that most of these mismatches could be attributed to the missing tidal forcing in the actual model configuration. This model is a first step to obtain a meaningful tool to study past and future oceanographic conditions in this important marine region constituting the unique connection of the Mediterranean Sea with the open world’s ocean.

Notes on HP1 a software package for simulating variably-saturated water flow, heat transport, solute transport, and biogeochemistry in porous media. HP1 Version 2.2

Energy Technology Data Exchange (ETDEWEB)

Jacques, D.; Simunek, J.

2010-01-15

HP1 is a comprehensive modeling tool in terms of processes and reactions for simulating reactive transport and biogeochemical processes in variably-saturated porous media. HP1 results from coupling the water and solute transport model HYDRUS-1D (Simunek et al., 2009a) and PHREEQC-2 (Parkhurst and Appelo, 1999). This note provides an overview of how to set up and execute a HP1 project using version 2.2.002 of HP1 and version 4.13 of the graphical user interface (GUI) of HYDRUS-1D. A large part of this note are step-by-step instructions for selected examples involving mineral dissolution and precipitation, cation exchange, surface complexation and kinetic degradation networks. The implementation of variably-saturated flow conditions, changing boundary conditions, a layered soil profile or immobile water is also illustrated.

Notes on HP1 a software package for simulating variably-saturated water flow, heat transport, solute transport, and biogeochemistry in porous media. HP1 Version 2.2

International Nuclear Information System (INIS)

Jacques, D.; Simunek, J.

2010-01-01

HP1 is a comprehensive modeling tool in terms of processes and reactions for simulating reactive transport and biogeochemical processes in variably-saturated porous media. HP1 results from coupling the water and solute transport model HYDRUS-1D (Simunek et al., 2009a) and PHREEQC-2 (Parkhurst and Appelo, 1999). This note provides an overview of how to set up and execute a HP1 project using version 2.2.002 of HP1 and version 4.13 of the graphical user interface (GUI) of HYDRUS-1D. A large part of this note are step-by-step instructions for selected examples involving mineral dissolution and precipitation, cation exchange, surface complexation and kinetic degradation networks. The implementation of variably-saturated flow conditions, changing boundary conditions, a layered soil profile or immobile water is also illustrated.

Biogeochemical cycles and continental ecosystems - Report on Science and Technology no. 27

International Nuclear Information System (INIS)

Pedro, Georges; Blanzat, Bernard; Albrecht, Pierre; Berthelin, Jacques; Boudot, Jean-Pierre; Munier-Lamy, Colette; Cossa, Daniel; Feix, Isabelle; Guillaumont, Robert; HUC, Alain Yves; Lavelle, Patrick; Lebrun, Michel; Lucas, Yves; Metivier, Henri; Ourisson, Guy; Raimbault, Patrick; Ranger, Jacques; Gerard, Frederic; Schmidt-Laine, Claudine; Dercourt, Jean; Gaillardet, Jerome; Bourrie, Guilhem; Trolard, Fabienne; Gerard, Frederic; Dambrine, Etienne; Meunier, Jean Dominique; Benoit, Marc; Breda, Nathalie; Dupouey, Jean-Luc; Granier, Andre; Franc, Alain; GARBAYE, Jean; Martin, Francis; Landmann, Guy; Loustau, Denis; Martinez, Jose; Crochon, Philippe; Gay, Jean-Didier; Peres, Jean-Marc; Tamponnet, Christian; Andreux, Francis; Tusseauvuillemin, Marie-Helene; Barker, Evelyne; Bouisset, Patrick; Germain, Pierre; Masson, Olivier; Boust, Dominique; Bailly du Bois, Pascal; Abdelouas, Abdesselam; Grambow, Bernd; Ansoborlo, Erich; Chiappini, Remo; Lobinski, Ryzsard; Montavon, Gilles; Moulin, Christophe; Moulin, Valerie; Ollivier, Bernard; Haeseler, Franck; Prieur, Daniel; Magot, Michel; Charmasson, Sabine; Poss, Roland; Grimaldi, Catherine; Grimaldi, Michel; Malet, Caroline

2007-11-01

environments is described very briefly to show that in this field everything at the surface of the planet is linked (chapter 5.2). In environments that are highly influenced by Man where the excess and overload of certain elements can lead to dysfunctions, we have focused on three types of environments: - the agricultural sector, where we looked at the practice of manuring using exogenous organic matter (EOM) which is found in areas of intensive breeding farms (chapter 6.1); - the industrial sector; we present a number of environments where certain elements are in excess, such as can happen in mining and metallurgy areas (chapter 6.2); - Finally, we looked at estuary environments because they are the collection point for numerous contaminants coming from the continents (chapter 6.3). The second unit proceeds from all we have learnt in the first, briefly that biogeochemical cycles are the consequence of complex systems that necessitate many scientific fields and numerous data. This consideration implies: a) further studies in certain scientific areas; b) a particular interest in data collection; c) the reliance on model building. Point a) constitutes the third part which stresses the development of all venues concerning the role of the soil as a living environment. Knowledge is still lacking in this area. Concerning this issue, two points are considered in particular: - research on micro-organisms that are crucial for transforming waste and other residues from living beings and their return to the mineral state (chapter 7); - research on soil organic matter, which plays a fundamental role because it is an intermediate between the mineral world and the living world (chapter 8). Points b) and c) are treated in the fourth part of the report which discusses model operations (chapter 9) and problems linked to a generalized collection of data (for instance, observatories, networks) (chapter 10)

Electronic transport in VO2—Experimentally calibrated Boltzmann transport modeling

International Nuclear Information System (INIS)

Kinaci, Alper; Rosenmann, Daniel; Chan, Maria K. Y.; Kado, Motohisa; Ling, Chen; Zhu, Gaohua; Banerjee, Debasish

2015-01-01

Materials that undergo metal-insulator transitions (MITs) are under intense study, because the transition is scientifically fascinating and technologically promising for various applications. Among these materials, VO 2 has served as a prototype due to its favorable transition temperature. While the physical underpinnings of the transition have been heavily investigated experimentally and computationally, quantitative modeling of electronic transport in the two phases has yet to be undertaken. In this work, we establish a density-functional-theory (DFT)-based approach with Hubbard U correction (DFT + U) to model electronic transport properties in VO 2 in the semiconducting and metallic regimes, focusing on band transport using the Boltzmann transport equations. We synthesized high quality VO 2 films and measured the transport quantities across the transition, in order to calibrate the free parameters in the model. We find that the experimental calibration of the Hubbard correction term can efficiently and adequately model the metallic and semiconducting phases, allowing for further computational design of MIT materials for desirable transport properties

Biological transportation networks: Modeling and simulation

KAUST Repository

Albi, Giacomo

2015-09-15

We present a model for biological network formation originally introduced by Cai and Hu [Adaptation and optimization of biological transport networks, Phys. Rev. Lett. 111 (2013) 138701]. The modeling of fluid transportation (e.g., leaf venation and angiogenesis) and ion transportation networks (e.g., neural networks) is explained in detail and basic analytical features like the gradient flow structure of the fluid transportation network model and the impact of the model parameters on the geometry and topology of network formation are analyzed. We also present a numerical finite-element based discretization scheme and discuss sample cases of network formation simulations.

Biogeochemical reactive-transport modelling of the interactions of medium activity long-lived nuclear waste in fractured argillite and the effect on redox conditions

International Nuclear Information System (INIS)

Small, J.S.; Steele, H.; Kwong, S.; Albrecht, A.

2010-01-01

Document available in extended abstract form only. The role of anaerobic microbial processes in mediating gas generation and redox reactions in organic (cellulose) containing low level activity nuclear wastes (LLW) is well established through monitoring of operational near-surface LLW disposal sites and municipal waste disposal sites. Modelling approaches based on Monod kinetic growth models to represent the complex suite of anaerobic processes have been developed and these models are able to reproduce the evolving biogeochemistry and gas generation of large scale and long term (10 year) experiments on cellulose waste degradation. In the case of geological disposal of medium activity long-lived nuclear waste (MAVL) microbial processes have the potential to exploit metabolic energy sources present in the waste, engineered barriers and host geological formation and as a consequence influence redox potential. Several electron donors and electron acceptors may be present in MAVL. Electron donors include; hydrogen (resulting from radiolysis and anaerobic corrosion of metals), and hydrolysis products of organic waste materials. Sulphate, nitrate and Fe(III) containing minerals and corrosion products are examples of electron acceptors present in intermediate level wastes. Significant amounts of organic matter, sulphate and iron minerals may also be present in host geological formations and have the potential to act as microbial energy sources once the system is perturbed by electron donors/acceptors from the waste. The construction of a geological disposal facility will physically disturb the host formation, potentially causing fracturing of the excavation damage zone (EDZ). The EDZ may thus provide environmental conditions, such as space and free water that together with nutrient and energy sources to promote microbial activity. In this study the Generalised Repository Model (GRM) developed to simulate the coupled microbiological, chemical and transport processes in near

Model Comparison for Electron Thermal Transport

Science.gov (United States)

Moses, Gregory; Chenhall, Jeffrey; Cao, Duc; Delettrez, Jacques

2015-11-01

Four electron thermal transport models are compared for their ability to accurately and efficiently model non-local behavior in ICF simulations. Goncharov's transport model has accurately predicted shock timing in implosion simulations but is computationally slow and limited to 1D. The iSNB (implicit Schurtz Nicolai Busquet electron thermal transport method of Cao et al. uses multigroup diffusion to speed up the calculation. Chenhall has expanded upon the iSNB diffusion model to a higher order simplified P3 approximation and a Monte Carlo transport model, to bridge the gap between the iSNB and Goncharov models while maintaining computational efficiency. Comparisons of the above models for several test problems will be presented. This work was supported by Sandia National Laboratory - Albuquerque and the University of Rochester Laboratory for Laser Energetics.

Biogeochemical cycle of boron in a forest ecosystem: the case study of Montiers beech-stand

International Nuclear Information System (INIS)

Roux, Philippe

2016-01-01

This thesis aims at establishing and understanding the biogeochemical cycle of boron and its isotopes within a forest ecosystem. In that context, many questions remain concerning the dynamics of boron within terrestrial ecosystems: - What are the major sources of boron? - What type of transfer occurs between the compartments of the environment? - What mechanisms are controlling those transfers? In order to establish this biogeochemical cycle, we quantified the different stocks (vegetation, humus and soil) and fluxes (atmospheric dust and dissolved deposition, throughfall, stem-flows, litterfall and drainage) of boron in the study site of Montiers. The use of boron isotopes will give us insight concerning the mechanisms controlling the dynamics of boron. This thesis is divided in 4 main parts: 1. The first part aims at establishing a new method of extraction, purification and measurement of boron and its isotopes within vegetation samples. 2. The second part focuses on the sources and mechanisms controlling boron within atmospheric dust and dissolved deposition on the study site of Montiers. 3. The third part aims at establishing the stocks and fluxes of boron on two distinct soils: a rendisoil (basic pH) and an alocrisoil (acid pH). The goal is to determine the influence of different soil properties on boron dynamics within its biogeochemical cycle. 4. The last part aims at establishing a model of boron and boron isotopes dynamics in the soil plant system. This model is mainly based of the measurement made in 2012. (author) [fr

Scaling Hydrologic Exchange Flows and Biogeochemical Reactions from Bedforms to Basins

Science.gov (United States)

Harvey, J. W.; Gomez-Velez, J. D.

2015-12-01

River water moves in and out of the main channel along pathways that are perpendicular to the channel's main axis that flow across or beneath the ground surface. These hydrologic exchange flows (HEFs) are difficult to measure, yet no less important than a river's downstream flow, or exchanges with the atmosphere and deeper groundwater (Harvey and Gooseff, 2015, WRR). There are very few comprehensive investigations of exchange fluxes to understand patterns with river size and relative importance of specific types of exchanges. We used the physically based model NEXSS to simulate multiple scales of hyporheic flow and their cumulative effects on solute reaction in large basins (on the order of Chesapeake Bay basin or larger). Our goal was to explain where and when particular types of hyporheic flow are important in enhancing key biogeochemical reactions, such as organic carbon respiration and denitrification. Results demonstrate that hyporheic flux (expressed per unit area of streambed) varies surprisingly little across the continuum of first-order streams to eighth-order rivers, and vertical exchange beneath small bedforms dominates in comparison with lateral flow beneath gravel bars and meanders. Also, the river's entire volume is exchanged many times with hyporheic flow within a basin, and the turnover length (after one entire river volume is exchanged) is strongly influenced by hydrogeomorphic differences between physiographic regions as well as by river size. The cumulative effects on biogeochemical reactions were assessed using a the reaction significance factor, RSF, which computes the cumulative potential for hyporheic reactions using a dimensionless index that balances reaction progress in a single hyporheic flow path against overall processing efficiency of river turnover through hyporheic flow paths of that type. Reaction significance appears to be strongly dominated by hydrologic factors rather than biogeochemical factors, and seems to be dominated by

Biogeochemical aspects of aquifer thermal energy storage

NARCIS (Netherlands)

Brons, H.J.

1992-01-01

During the process of aquifer thermal energy storage the in situ temperature of the groundwater- sediment system may fluctuate significantly. As a result the groundwater characteristics can be considerably affected by a variety of chemical, biogeochemical and microbiological