Final Technical Report: Viral Infection of Subsurface Microorganisms and Metal/Radionuclide Transport

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

Weber, Karrie A.; Bender, Kelly S.; Li, Yusong

Microbially mediated metabolisms have been identified as a significant factor either directly or indirectly impacting the fate and transport of heavy metal/radionuclide contaminants. To date microorganisms have been isolated from contaminated environments. Examination of annotated finished genome sequences of many of these subsurface isolates from DOE sites, revealed evidence of prior viral infection. To date the role that viruses play influencing microbial mortality and the resulting community structure which directly influences biogeochemical cycling in soils and sedimentary environments remains poorly understood. The objective of this exploratory study was to investigate the role of viral infection of subsurface bacteria and themore » formation of contaminant-bearing viral particles. This objective was approached by examining the following working hypotheses: (i) subsurface microorganisms are susceptible to viral infections by the indigenous subsurface viral community, and (ii) viral surfaces will adsorb heavy metals and radionuclides. Our results have addressed basic research needed to accomplish the BER Long Term Measure to provide sufficient scientific understanding such that DOE sites would be able to incorporate coupled physical, chemical and biological processes into decision making for environmental remediation or natural attenuation and long-term stewardship by establishing viral-microbial relationships on the subsequent fate and transport of heavy metals and radionuclides. Here we demonstrated that viruses play a significant role in microbial mortality and community structure in terrestrial subsurface sedimentary systems. The production of viral-like particles within subsurface sediments in response to biostimulation with dissolved organic carbon and a terminal electron acceptor resulted in the production of viral-like particles. Organic carbon alone did not result in significant viral production and required the addition of a terminal electron acceptor (nitrate), indicating that nutrients are not limiting viral production, but rather substrates that can be converted into energy for host metabolism. Our results also revealed that cell abundance was not correlated to the mineralization of organic carbon, but rather viruses were positively correlated with carbon mineralization. This is a result of viral-mediated cell lysis and demonstrates that viruses are sensitive indicators of microbial activity. Viruses as an indicator of microbial activity was not unique to batch culture studies as results obtained from an in situ field experiment conducted at the DOE Old Rifle Field site. This study revealed that viral abundance increased in response to the injection of oxygenated groundwater and influx of dissolved organic carbon whereas cell abundance changes were minimal. However, the extent to which viral-mediated cell lysis alters organic matter pools subsequently influencing microbial community structure and biogeochemical function remains a critical question in subsurface biogeochemical cycling. The production of significant numbers of viruses in groundwater has implications for nanoparticulate metal as well as carbon transport in groundwater. We have demonstrated that the virus surface is reactive and will adsorb heavy metals. Thus viruses can promote colloidal contaminant mobility. Interestingly, the presence of heavy metals has a positive effect on infectivity of the phage, increasing phage infection which could lead to further production of viruses. Together, the results indicate that the sorption of metals to the surface of viruses could not only contribute to nanoparticulate metal as well as carbon transport but could also enhance infectivity further contributing to cell lysis which could subsequently influence biogeochemical cycling. As more viruses infect host microbial populations the high concentration of metals would enhance infection, resulting in cell lysis, and decreasing the metabolically active host population while yielding greater numbers of viruses capable of transporting contaminats. Additional studies will be necessary to further establish the potential relationship(s) between viruses, cells, carbon, and metals/radionuclides to provide sufficient scientific understanding to incorporate coupled physical, chemical, and biological processes into agent based and reactive transport models.« less

Quantifying the surface subsurface biogeochemical coupling during the VERTIGO ALOHA and K2 studies

NASA Astrophysics Data System (ADS)

Boyd, Philip W.; Gall, Mark P.; Silver, Mary W.; Coale, Susan L.; Bidigare, Robert R.; Bishop, James L. K. B.

2008-07-01

A central question addressed by the VERtical Transport In the Global Ocean (VERTIGO) study was 'What controls the efficiency of particle export between the surface and subsurface ocean'? Here, we present data from sites at ALOHA (N Central Pacific Gyre) and K2 (NW subarctic Pacific) on phytoplankton processes, and relate them via a simple planktonic foodweb model, to subsurface particle export (150-500 m). Three key factors enable quantification of the surface-subsurface coupling: a sampling design to overcome the temporal lag and spatial displacement between surface and subsurface processes; data on the size partitioning of net primary production (NPP) and subsequent transformations prior to export; estimates of the ratio of algal- to faecal-mediated vertical export flux. At ALOHA, phytoplankton were characterized by low stocks, NPP, Fv/ Fm (N-limited), and were dominated by picoplankton. The HNLC waters at K2 were characterized by both two-fold changes in NPP and floristic shifts (high to low proportion of diatoms) between deployment 1 and 2. Prediction of export exiting the euphotic zone was based on size partitioning of NPP, a copepod-dominated foodweb and a ratio of 0.2 (ALOHA) and 0.1 (K2) for algal:faecal particle flux. Predicted export was 20-22 mg POC m -2 d -1 at ALOHA (i.e. 10-11% NPP (0-125 m); 1.1-1.2×export flux at 150 m ( E150). At K2, export was 111 mg C m -2 d -1 (21% NPP (0-50 m); 1.8× E150) and 33 mg POC m -2 d -1 (11% NPP, 0-55 m); 1.4× E150) for deployments 1 and 2, respectively. This decrease in predicted export at K2 matches the observed trend for E150. Also, the low attenuation of export flux from 60 to 150 m is consistent with that between 150 and 500 m. This strong surface-subsurface coupling suggests that phytoplankton productivity and floristics play a key role at K2 in setting export flux, and moreover that pelagic particle transformations by grazers strongly influence to what extent sinking particles are further broken down in the underlying waters of the Twilight Zone.

Control of topography gradients on residence time distributions, mixing dynamics and reactive hotspot development

NASA Astrophysics Data System (ADS)

Bandopadhyay, Aditya; Le Borgne, Tanguy; Davy, Philippe

2017-04-01

Topography-driven subsurface flows are thought to play a central role in determining solute turnover and biogeochemical processes at different scales in the critical zone, including river-hyporheic zone exchanges, hillslope solute transport and reactions, and catchment biogeochemical cycles. Hydraulic head gradients, induced by topography gradients at different scales, generate a distribution of streamlines at depth, dictating the spatial distribution of redox sensitive species, the magnitude of surface water - ground water exchanges and ultimately the source/sink function of the subsurface. Flow velocities generally decrease with depth, leading to broad residence time distributions, which have been shown to affect river chemistry and geochemical reactions in catchments. In this presentation, we discuss the impact of topography-driven flows on mixing processes and the formation of localized reactive hotspots. For this, we solve analytically the coupled flow, mixing and reaction equations in two-dimensional vertical cross-sections of subsurface domains with different topography gradients. For a given topography gradient, we derive the spatial distribution of subsurface velocities, the rates of solute mixing accross streamlines and the induced kinetics of redox, precipitation and dissolution reactions using a Lagrangian approach (Le Borgne et al. 2014). We demonstrate that vertical velocity profiles driven by topography variations, act effectively as shear flows, hence stretching continuously the mixing fronts between recently infiltrated and resident water (Bandopadhyay et al. 2017). We thus derive analytical expressions for residence time distributions, mixing rates and kinetics of chemical reactions as a function of the topography gradients. We show that the rates dissolution and precipitation reactions are significantly enhanced by the existence of vertical velocity gradients and that reaction rates reach a maximum in a localized subsurface reactive layer, whose location and intensity depends on topography gradients. As a consequence of these findings, we discuss the links between topography variations, subsurface velocity gradients and biogeochemical processes in the critical zone. References: Bandopadhyay A., T. Le Borgne, Y. Méheust and M. Dentz (2017) Enhanced reaction kinetics and reactive mixing scale dynamics in mixing fronts under shear flow for arbitrary Damkohler numbers, Adv. in Water Resour. Vol. 100, p. 78-95 Le Borgne T., T. Ginn and M. Dentz (2014) Impact of Fluid Deformation on Mixing-Induced Chemical Reactions in Heterogeneous Flows, Geophys. Res. Lett., Vol. 41, 22, p. 7898-790

Microfluidic Experiments Studying Pore Scale Interactions of Microbes and Geochemistry

NASA Astrophysics Data System (ADS)

Chen, M.; Kocar, B. D.

2016-12-01

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

A spatially explicit hydro-ecological modeling framework (BEPS-TerrainLab V2.0): Model description and test in a boreal ecosystem in Eastern North America

NASA Astrophysics Data System (ADS)

Govind, Ajit; Chen, Jing Ming; Margolis, Hank; Ju, Weimin; Sonnentag, Oliver; Giasson, Marc-André

2009-04-01

SummaryA spatially explicit, process-based hydro-ecological model, BEPS-TerrainLab V2.0, was developed to improve the representation of ecophysiological, hydro-ecological and biogeochemical processes of boreal ecosystems in a tightly coupled manner. Several processes unique to boreal ecosystems were implemented including the sub-surface lateral water fluxes, stratification of vegetation into distinct layers for explicit ecophysiological representation, inclusion of novel spatial upscaling strategies and biogeochemical processes. To account for preferential water fluxes common in humid boreal ecosystems, a novel scheme was introduced based on laboratory analyses. Leaf-scale ecophysiological processes were upscaled to canopy-scale by explicitly considering leaf physiological conditions as affected by light and water stress. The modified model was tested with 2 years of continuous measurements taken at the Eastern Old Black Spruce Site of the Fluxnet-Canada Research Network located in a humid boreal watershed in eastern Canada. Comparison of the simulated and measured ET, water-table depth (WTD), volumetric soil water content (VSWC) and gross primary productivity (GPP) revealed that BEPS-TerrainLab V2.0 simulates hydro-ecological processes with reasonable accuracy. The model was able to explain 83% of the ET, 92% of the GPP variability and 72% of the WTD dynamics. The model suggests that in humid ecosystems such as eastern North American boreal watersheds, topographically driven sub-surface baseflow is the main mechanism of soil water partitioning which significantly affects the local-scale hydrological conditions.

Lateral, vertical, and longitudinal connectivity of runoff source areas drive stream hydro-biogeochemical signals across a low relief drainage network

NASA Astrophysics Data System (ADS)

Zimmer, M. A.; McGlynn, B. L.

2017-12-01

Our understanding of the balance between longitudinal, lateral, and vertical expansion and contraction of reactive flowpaths and source areas in headwater catchments is limited. To address this, we utilized an ephemeral-to-perennial stream network in the Piedmont region of North Carolina, USA to gain new understanding about critical zone mechanisms that drive runoff generation and biogeochemical signals in both groundwater and stream water. Here, we used chemical and hydrometric data collected from zero through second order catchments to characterize spatial and temporal runoff and overland, shallow soil, and deep subsurface flow across characteristic landscape positions. Our results showed that the active stream network was driven by two superimposed runoff generation regimes that produced distinct hydro-biogeochemical signals at the catchment outlet. The baseflow runoff generation regime expanded and contracted the stream network seasonally through the rise and fall of the seasonal water table. Superimposed on this, event-activated source area contributions were driven by surficial and shallow subsurface flowpaths. The subsurface critical zone stratigraphy in this landscape coupled with the precipitation regime activated these shallow flowpaths frequently. This drove an increase in dissolved organic carbon (DOC) concentrations with increases in runoff across catchment scales. DOC-runoff relationship variability and spread was driven by the balance between runoff regimes as well as a seasonal depletion of DOC from shallow subsurface flowpath activation and annual replenishment from litterfall. From this, we suggest that the hydro-biogeochemical signals at larger catchment outlets can be driven by a balance of longitudinal, lateral, and vertical source area contributions, critical zone structure, and complex hydrological processes.

Establishment of quantitative hydrological indexes for studies of hydro-biogeochemical interactions at the subsurface.

NASA Astrophysics Data System (ADS)

Alves Meira Neto, A.; Sengupta, A.; Wang, Y.; Volkmann, T.; Chorover, J.; Troch, P. A. A.

2017-12-01

Advances in the understanding of processes in the critical zone (CZ) are dependent on studies coupling the fields of hydrology, microbiology, geochemistry and soil development. At the same time, better insights are needed to integrate hydrologic information into biogeochemical analysis of subsurface environments. This study investigated potential hydrological indexes that help explaining spatiotemporal biogeochemical patterns. The miniLEO is a 2 m3, 10 degree sloping lysimeter located at Biosphere 2 - University of Arizona. The lysimeter was initially filled with pristine basaltic soil and subject to intermittent rainfall applications throughout the period of 18 months followed by its excavation, resulting in a grid-based sample collection at 324 locations. As a result, spatially distributed microbiological and geochemical patterns as well as soil physical properties were obtained. A hydrologic model was then developed in order to simulate the history of the system until the excavation. After being calibrated against sensor data to match its observed input-state-output behavior, the resulting distributed fields of flow velocities and moisture states were retrieved. These results were translated into several hydrological indexes to be used in with distributed microbiological and geochemical signatures. Our study attempts at conciliating sound hydrological modelling with an investigation of the subsurface biological signatures, thus providing a unique opportunity for understanding of fine-scale hydro-biological interactions.

Thinking outside the channel: Modeling nitrogen cycling in networked river ecosystems

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

Helton, Ashley; Poole, Geoffrey C.; Meyer, Judy

2011-01-01

Agricultural and urban development alters nitrogen and other biogeochemical cycles in rivers worldwide. Because such biogeochemical processes cannot be measured empirically across whole river networks, simulation models are critical tools for understanding river-network biogeochemistry. However, limitations inherent in current models restrict our ability to simulate biogeochemical dynamics among diverse river networks. We illustrate these limitations using a river-network model to scale up in situ measures of nitrogen cycling in eight catchments spanning various geophysical and land-use conditions. Our model results provide evidence that catchment characteristics typically excluded from models may control river-network biogeochemistry. Based on our findings, we identify importantmore » components of a revised strategy for simulating biogeochemical dynamics in river networks, including approaches to modeling terrestrial-aquatic linkages, hydrologic exchanges between the channel, floodplain/riparian complex, and subsurface waters, and interactions between coupled biogeochemical cycles.« less

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

NASA Astrophysics Data System (ADS)

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

2018-06-01

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

Numerical modeling of watershed-scale radiocesium transport coupled with biogeochemical cycling in forests

NASA Astrophysics Data System (ADS)

Mori, K.; Tada, K.; Tawara, Y.; Tosaka, H.; Ohno, K.; Asami, M.; Kosaka, K.

2015-12-01

Since the Fukushima Dai-ichi Nuclear Power Plant (FDNPP) accident, intensive monitoring and modeling works on radionuclide transfer in environment have been carried out. Although Cesium (Cs) concentration has been attenuating due to both physical and environmental half-life (i.e., wash-off by water and sediment), the attenuation rate depends clearly on the type of land use and land cover. In the Fukushima case, studying the migration in forest land use is important for predicting the long-term behavior of Cs because most of the contaminated region is covered by forests. Atmospheric fallout is characterized by complicated behavior in biogeochemical cycle in forests which can be described by biotic/abiotic interactions between many components. In developing conceptual and mathematical model on Cs transfer in forest ecosystem, defining the dominant components and their interactions are crucial issues (BIOMASS, 1997-2001). However, the modeling of fate and transport in geosphere after Cs exports from the forest ecosystem is often ignored. An integrated watershed modeling for simulating spatiotemporal redistribution of Cs that includes the entire region from source to mouth and surface to subsurface, has been recently developed. Since the deposited Cs can migrate due to water and sediment movement, the different species (i.e., dissolved and suspended) and their interactions are key issues in the modeling. However, the initial inventory as source-term was simplified to be homogeneous and time-independent, and biogeochemical cycle in forests was not explicitly considered. Consequently, it was difficult to evaluate the regionally-inherent characteristics which differ according to land uses, even if the model was well calibrated. In this study, we combine the different advantages in modeling of forest ecosystem and watershed. This enable to include more realistic Cs deposition and time series of inventory can be forced over the land surface. These processes are integrated into the watershed simulator GETFLOWS coupled with biogeochemical cycling in forests. We present brief a overview of the simulator and an application for reservoir basin.

Preliminary forecasts of Pacific bigeye tuna population trends under the A2 IPCC scenario

NASA Astrophysics Data System (ADS)

Lehodey, P.; Senina, I.; Sibert, J.; Bopp, L.; Calmettes, B.; Hampton, J.; Murtugudde, R.

2010-07-01

An improved version of the spatial ecosystem and population dynamics model SEAPODYM was used to investigate the potential impacts of global warming on tuna populations. The model included an enhanced definition of habitat indices, movements, and accessibility of tuna predators to different vertically migrant and non-migrant micronekton functional groups. The simulations covered the Pacific basin (model domain) at a 2° × 2° geographic resolution. The structure of the model allows an evaluation from multiple data sources, and parameterization can be optimized by adjoint techniques and maximum likelihood using fishing data. A first such optimized parameterization was obtained for bigeye tuna ( Thunnus obesus) in the Pacific Ocean using historical catch data for the last 50 years and a hindcast from a coupled physical-biogeochemical model driven by the NCEP atmospheric reanalysis. The parameterization provided very plausible biological parameter values and a good fit to fishing data from the different fisheries, both within and outside the time period used for optimization. We then employed this model to forecast the future of bigeye tuna populations in the Pacific Ocean. The simulation was driven by the physical-biogeochemical fields predicted from a global marine biogeochemistry - climate simulation. This global simulation was performed with the IPSL climate model version 4 (IPSL-CM4) coupled to the oceanic biogeochemical model PISCES and forced by atmospheric CO 2, from historical records over 1860-2000, and under the SRES A2 IPCC scenario for the 21st century (i.e. atmospheric CO 2 concentration reaching 850 ppm in the year 2100). Potential future changes in distribution and abundance under the IPCC scenario are presented but without taking into account any fishing effort. The simulation showed an improvement in bigeye tuna spawning habitat both in subtropical latitudes and in the eastern tropical Pacific (ETP) where the surface temperature becomes optimal for bigeye tuna spawning. The adult feeding habitat also improved in the ETP due to the increase of dissolved oxygen concentration in the sub-surface allowing adults to access deeper forage. Conversely, in the Western Central Pacific the temperature becomes too warm for bigeye tuna spawning. The decrease in spawning is compensated by an increase of larvae biomass in subtropical regions. However, natural mortality of older stages increased due to lower habitat values (too warm surface temperatures, decreasing oxygen concentration in the sub-surface and less food). This increased mortality and the displacement of surviving fish to the eastern region led to stable then declining adult biomass at the end of the century.

Hydrogeochemical zonation in intertidal salt marsh sediments: evidence of positive plant-soil feedback?

NASA Astrophysics Data System (ADS)

Moffett, K. B.; Dittmar, J.; Seyfferth, A.; Fendorf, S.; Gorelick, S.

2012-12-01

Surface and subsurface environments are linked by the biogeochemical activity in near-surface sediment and by the hydrological fluxes that mobilize its reagents and products. A particularly dynamic and interesting setting to study near-surface hydrogeochemistry is the intertidal zone. Here, the very strong tidal hydraulic forcing is often thought to dominate water and solute transport. However, we demonstrated the importance of two additional subsurface drivers: groundwater flow and plant root water uptake. A high-resolution, coupled surface water-groundwater model of an intertidal salt marsh in San Francisco Bay, CA showed that these three drivers vary over different spatial scales: tidal flooding varies over 10's of meters; groundwater flow varies over meters, particularly within channel banks; and plant root water uptake varies in 3D at the sub-meter scale. Expanding on this third driver, we investigated whether the spatial variations in soil-water-plant hydraulic interactions that occur due to vegetation zonation also cause distinct geochemical zonation in salt marsh sediment pore waters. The existence of such geochemical zonation was verified and mapped by detailed field observations of the chemical composition of sediments, pore waters, surface waters, and vegetation. The field data and the coupled hydrologic model were then further analyzed to evaluate potential causal mechanisms for the geochemical zonation, including testing the hypothesis that the vegetation affects pore water geochemistry via a positive feedback beneficial to itself. If further supported by future studies, this geochemical feedback may complement known physical ecosystem engineering mechanisms to help stabilize and organize intertidal wetlands.

Hyporheic Interfaces Serve as Ecological Control Points for Mountainous Landscape Biological Productivity

NASA Astrophysics Data System (ADS)

Newcomer, M. E.; Dwivedi, D.; Raberg, J.; Fox, P. M.; Nico, P. S.; Wainwright, H. M.; Conrad, M. E.; Bill, M.; Bouskill, N.; Williams, K. H.; Hubbard, S.; Steefel, C. I.

2017-12-01

Riverine systems in snow-dominated mountainous regions often express complex biogeochemistry and river nutrient indicators as a function of hydrologic variability. In early spring, meltwater infiltration from a ripened snowpack creates a hydrological gradient through hillslopes, floodplains, and hyporheic zones. During this time, these systems are more-or-less a passive filter that allows the rising limb of the hydrograph to display chemo-dynamic relationships (inversely proportional) with solutes and nutrients. During the growing season, temperatures, plants, microbes, and hydrologic gradients shift dramatically and activate hyporheic-zone biogeochemistry as a major control on water nutrient degradation. Hyporheic biogeochemical reliance on the timing of meltwater infiltration and the possibility of a longer vernal window under future climate change indicates the importance of hyporheic cycling as the dominant ecological control point on carbon and nitrogen fluxes and transformations. The objective of our study is to develop a predictive understanding of the subsurface and surface controls on hyporheic biogeochemical behavior through data-model integration. Data from our 2017 field campaign in the East River, Colorado, a pristine, mountainous watershed, were taken at key times during the rising, peak, falling, and dry limb of the hydrograph. Throughout multiple locations across this spatial and temporal gradient, we measured surface and subsurface gases, geochemistry, isotopes, and hydrological flow conditions and used this data to constrain a numerical flow and reactive transport model of the hyporheic zone that included microbial and flow feedback dynamics. Our data coupled with the predictive power of our numerical model reveal that the hyporheic zone serves dual roles throughout the year—as a net source of nutrients and solutes during the early vernal phase, shifting to a net sink of nutrients during the summer dry season. The possibility of a future lengthened vernal window motivates a better understanding of the role each ecological control point plays in processing landscape biological productivity and for understanding biogeochemical cycling in riverine systems.

Quantifying Linkages between Biogeochemical Processes in a Contaminated Aquifer-Wetland System Using Multivariate Statistics and HP1

NASA Astrophysics Data System (ADS)

Arora, B.; Mohanty, B. P.; McGuire, J. T.

2009-12-01

Fate and transport of contaminants in saturated and unsaturated zones in the subsurface is controlled by complex biogeochemical processes such as precipitation, sorption-desorption, ion-exchange, redox, etc. In dynamic systems such as wetlands and anaerobic aquifers, these processes are coupled and can interact non-linearly with each other. Variability in measured hydrological, geochemical and microbiological parameters thus corresponds to multiple processes simultaneously. To infer the contributing processes, it is important to eliminate correlations and to identify inter-linkages between factors. The objective of this study is to develop quantitative relationships between hydrological (initial and boundary conditions, hydraulic conductivity ratio, and soil layering), geochemical (mineralogy, surface area, redox potential, and organic matter) and microbiological factors (MPN) that alter the biogeochemical processes at the column scale. Data used in this study were collected from controlled flow experiments in: i) two homogeneous soil columns, ii) a layered soil column, iii) a soil column with embedded clay lenses, and iv) a soil column with embedded clay lenses and one central macropore. The soil columns represent increasing level of soil structural heterogeneity to better mimic the Norman Landfill research site. The Norman Landfill is a closed municipal facility with prevalent organic contamination. The sources of variation in the dataset were explored using multivariate statistical techniques and dominant biogeochemical processes were obtained using principal component analysis (PCA). Furthermore, artificial neural networks (ANN) coupled with HP1 was used to develop mathematical rules identifying different combinations of factors that trigger, sustain, accelerate/decelerate, or discontinue the biogeochemical processes. Experimental observations show that infiltrating water triggers biogeochemical processes in all soil columns. Similarly, slow release of water from low permeability clay lenses sustain biogeochemical cycling for a longer period of time than in homogeneous soil columns. Preliminary results indicate: i) certain variables (anion, cation concentrations, etc.) do not follow normal or lognormal distributions even at the column scale, ii) strong correlations exist between parameters related to redox geochemistry (pH with S2- concentrations), and iii) PCA can identify dominant processes (e.g. iron and sulfate reduction) occurring in the system by grouping together causative variables (e.g. dominant TEAPs).

Improving Permafrost Hydrology Prediction Through Data-Model Integration

NASA Astrophysics Data System (ADS)

Wilson, C. J.; Andresen, C. G.; Atchley, A. L.; Bolton, W. R.; Busey, R.; Coon, E.; Charsley-Groffman, L.

2017-12-01

The CMIP5 Earth System Models were unable to adequately predict the fate of the 16GT of permafrost carbon in a warming climate due to poor representation of Arctic ecosystem processes. The DOE Office of Science Next Generation Ecosystem Experiment, NGEE-Arctic project aims to reduce uncertainty in the Arctic carbon cycle and its impact on the Earth's climate system by improved representation of the coupled physical, chemical and biological processes that drive how much buried carbon will be converted to CO2 and CH4, how fast this will happen, which form will dominate, and the degree to which increased plant productivity will offset increased soil carbon emissions. These processes fundamentally depend on permafrost thaw rate and its influence on surface and subsurface hydrology through thermal erosion, land subsidence and changes to groundwater flow pathways as soil, bedrock and alluvial pore ice and massive ground ice melts. LANL and its NGEE colleagues are co-developing data and models to better understand controls on permafrost degradation and improve prediction of the evolution of permafrost and its impact on Arctic hydrology. The LANL Advanced Terrestrial Simulator was built using a state of the art HPC software framework to enable the first fully coupled 3-dimensional surface-subsurface thermal-hydrology and land surface deformation simulations to simulate the evolution of the physical Arctic environment. Here we show how field data including hydrology, snow, vegetation, geochemistry and soil properties, are informing the development and application of the ATS to improve understanding of controls on permafrost stability and permafrost hydrology. The ATS is being used to inform parameterizations of complex coupled physical, ecological and biogeochemical processes for implementation in the DOE ACME land model, to better predict the role of changing Arctic hydrology on the global climate system. LA-UR-17-26566.

Linking Metabolism, Elemental Cycles, and Environmental Conditions in the Deep Biosphere: Growth of a Model Extremophile, Archaeoglobus fulgidus, Under High-Pressure Conditions

NASA Astrophysics Data System (ADS)

Oliver, G. C. M.; Cario, A.; Rogers, K. L.

2015-12-01

A majority of Earth's biosphere is hosted in subsurface environments where global-scale biogeochemical and energy cycles are driven by diverse microbial communities that operate on and are influenced by micro-scale environmental variables. While the subsurface hosts a variety of geochemical and geothermal conditions, elevated pressures are common to all subsurface ecosystems. Understanding how microbes adapt to and thrive in high-pressure environments is essential to linking microbial subsurface processes with global-scale cycles. Here we are using a model extremophile, Archaeoglobus fulgidus, to determine how elevated pressures affect the growth, metabolism, and physiology of subsurface microorganisms. A. fulgidus cycles carbon and sulfur via heterotrophic and autotrophic sulfate reduction in various high temperature and high-pressure niches including shallow marine vents, deep-sea hydrothermal vents, and deep oil reservoirs. Here we report the results of A. fulgidus growth experiments at optimum temperature, 83°C, and pressures up to 600 bars. Exponential growth was observed over the entire pressure range, though growth rates were diminished at 500 and 600 bars compared to ambient pressure experimental controls. At pressures up to 400 bars, cell density yields and growth rates were at least as high as ambient pressure controls. Elevated pressures and extended incubation times stimulated cell flocculation, a common stress response in this strain, and cellular morphology was affected at pressures exceeding 400 bars. These results suggest that A. fulgidus continues carbon, sulfur and energy cycling unaffected by elevated pressures up to 400 bars, representing a variety of subsurface environments. The ability of subsurface organisms to drive biogeochemical cycles at elevated pressures is a critical link between the surface and subsurface biospheres and understanding how species-scale processes operate under these conditions is a vital part of global-scale biogeochemical models.

Advances in interpretation of subsurface processes with time-lapse electrical imaging

USGS Publications Warehouse

Singha, Kaminit; Day-Lewis, Frederick D.; Johnson, Tim B.; Slater, Lee D.

2015-01-01

Electrical geophysical methods, including electrical resistivity, time-domain induced polarization, and complex resistivity, have become commonly used to image the near subsurface. Here, we outline their utility for time-lapse imaging of hydrological, geochemical, and biogeochemical processes, focusing on new instrumentation, processing, and analysis techniques specific to monitoring. We review data collection procedures, parameters measured, and petrophysical relationships and then outline the state of the science with respect to inversion methodologies, including coupled inversion. We conclude by highlighting recent research focused on innovative applications of time-lapse imaging in hydrology, biology, ecology, and geochemistry, among other areas of interest.

Advances in interpretation of subsurface processes with time-lapse electrical imaging

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

Singha, Kamini; Day-Lewis, Frederick D.; Johnson, Timothy C.

2015-03-15

Electrical geophysical methods, including electrical resistivity, time-domain induced polarization, and complex resistivity, have become commonly used to image the near subsurface. Here, we outline their utility for time-lapse imaging of hydrological, geochemical, and biogeochemical processes, focusing on new instrumentation, processing, and analysis techniques specific to monitoring. We review data collection procedures, parameters measured, and petrophysical relationships and then outline the state of the science with respect to inversion methodologies, including coupled inversion. We conclude by highlighting recent research focused on innovative applications of time-lapse imaging in hydrology, biology, ecology, and geochemistry, among other areas of interest.

Numerical Modeling of Artificial Recharge: Determining Spatial/Temporal Sampling Resolution to Quantify Infiltration Rates and Effective Hydraulic Conductivity

NASA Astrophysics Data System (ADS)

Glose, T. J.; Hausner, M. B.; Lowry, C.

2016-12-01

The accurate, fine scale quantification of groundwater-surface water (GW-SW) interactions over large expanses in hydrologic systems is a fundamental need in order to accurately characterize critical zones of biogeochemical transformation and fluxes, as well as to provide insight into near-surface geologic heterogeneity. Paired fiber-optic distributed temperature sensing (FO-DTS) is a tool that is capable of synoptically sampling hydrologic systems, allowing GW-SW interactions to be examined at a fine scale over large distances. Within managed aquifer recharge (MAR) sites, differential recharge dynamics controlled by bed clogging and subsurface heterogeneity dictate the effectiveness of these sites at infiltrating water. Numerical modeling indicates that the use of paired FO-DTS in an MAR site can provide accurate quantification of flux at the GW-SW interface, as well as provide insight to the areal extent of geologic heterogeneity in the subsurface. However, the lateral and vertical separation of the fiber-optic cables is of vital importance. Here we present a 2-D, fully coupled groundwater flow and heat transport model with prescribed heterogeneity. Following a forward modeling approach, realizations simulating varying fiber-optic cable positioning, differential bed clogging, and hydraulic conductivity variability were analyzed over a suite of scenarios. The results from the model were then used as observations to calculate groundwater recharge rates and calibration targets for an inverse model to estimate subsurface heterogeneity.

Using lagged dependence to identify (de)coupled surface and subsurface soil moisture values

NASA Astrophysics Data System (ADS)

Carranza, Coleen D. U.; van der Ploeg, Martine J.; Torfs, Paul J. J. F.

2018-04-01

Recent advances in radar remote sensing popularized the mapping of surface soil moisture at different spatial scales. Surface soil moisture measurements are used in combination with hydrological models to determine subsurface soil moisture values. However, variability of soil moisture across the soil column is important for estimating depth-integrated values, as decoupling between surface and subsurface can occur. In this study, we employ new methods to investigate the occurrence of (de)coupling between surface and subsurface soil moisture. Using time series datasets, lagged dependence was incorporated in assessing (de)coupling with the idea that surface soil moisture conditions will be reflected at the subsurface after a certain delay. The main approach involves the application of a distributed-lag nonlinear model (DLNM) to simultaneously represent both the functional relation and the lag structure in the time series. The results of an exploratory analysis using residuals from a fitted loess function serve as a posteriori information to determine (de)coupled values. Both methods allow for a range of (de)coupled soil moisture values to be quantified. Results provide new insights into the decoupled range as its occurrence among the sites investigated is not limited to dry conditions.

Statistics of chemical gradients in heterogeneous porous media

NASA Astrophysics Data System (ADS)

Le Borgne, T.; Huck, P. D.; Dentz, M.; Villermaux, E.

2017-12-01

As they create chemical disequilibrium and drive mixing fluxes, spatial gradients in solute concentrations exert a strong control on mixing and biogeochemical reactions in the subsurface. Large concentration gradients may develop in particular at interfaces between surface water and groundwater bodies, such as hyporheic zones, sea water - surface water interfaces or recharge areas. They also develop around contaminant plumes and fluids injected in subsurface operations. While macrodispersion theories predict smooth gradients, decaying in time due to dispersive dissipation, we show that concentration gradients are sustained by flow heterogeneity and have broadly distributed values. We present a general theory predicting the statistics of concentration gradients from the flow heterogeneity (Le Borgne et al., 2017). Analytical predictions are validated from high resolution simulations of transport in heterogeneous Darcy fields ranging from low to high permeability variances and low to high Peclet numbers. This modelling framework hence opens new perspectives for quantifying the dynamics of chemical gradients and the kinetics of associated biogeochemical reactions in heterogeneous subsurface environments.Reference:Le Borgne T., P.D. Huck, M. Dentz and E. Villermaux (2017) Scalar gradients in stirred mixtures and the deconstruction of random fields, J. of Fluid Mech. vol. 812, pp. 578-610 doi:10.1017/jfm.2016.799

Concentration-Discharge Relationships, Nested Reaction Fronts, and the Balance of Oxidative and Acid-Base Weathering Fluxes in an Alpine Catchment, East River, Colorado

NASA Astrophysics Data System (ADS)

Winnick, M.; Carroll, R. W. H.; Williams, K. H.; Maxwell, R. M.; Maher, K.

2016-12-01

Although important for solute production and transport, the varied interactions between biogeochemical processes and subsurface hydrology remain poorly characterized. We investigate these couplings in the headwaters of the East River, CO, a high-elevation shale-dominated catchment system in the Rocky Mountains, using concentration-discharge (C-Q) relationships for major cations, anions, and organic carbon. Dissolved organic carbon (DOC) displays a positive C-Q relationship with well-defined clockwise hysteresis, indicating the mobilization and depletion of DOC in the upper soil horizons and highlighting the importance of shallow flowpaths through the snowmelt period. Cation and anion concentrations demonstrate that carbonate weathering, which dominates solute fluxes, is promoted by both carbonic acid and sulfuric acid derived from oxidation of pyrite in the shale bedrock. Sulfuric acid weathering in the deep subsurface dominates during base flow conditions when waters have infiltrated below the hypothesized pyrite oxidation front, whereas carbonic acid weathering plays a dominant role during the snowmelt period as a result of shallow flowpaths. Differential C-Q relationships between solutes suggest that infiltrating waters approach calcite saturation before reaching the pyrite oxidation front, after which sulfuric acid reduces carbonate alkalinity. This increase in CO2(aq) at the expense of HCO3- results in outgassing of CO2 when waters equilibrate to surface conditions, and reduces the export of carbon and alkalinity from the East River by roughly 33% annually. Future changes in snowmelt dynamics that control the balance of carbonic and sulfuric acid weathering therefore have the capacity to substantially alter the cycling of carbon in the East River catchment. Ultimately, we demonstrate that differential C-Q relationships between major solutes can provide unique insights into the complex subsurface flow and biogeochemical dynamics that operate at catchment scales.

Snowmelt controls on concentration-discharge relationships and the balance of oxidative and acid-base weathering fluxes in an alpine catchment, East River, Colorado: ACID-BASE VERSUS OXIDATIVE WEATHERING FLUXES

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

Winnick, Matthew J.; Carroll, Rosemary W. H.; Williams, Kenneth H.

Although important for riverine solute and nutrient fluxes, the connections between biogeochemical processes and subsurface hydrology remain poorly characterized. We investigate these couplings in the East River, CO, a high-elevation shale-dominated catchment in the Rocky Mountains, using concentration-discharge (C-Q) relationships for major cations, anions, and organic carbon. Dissolved organic carbon (DOC) displays a positive C-Q relationship with clockwise hysteresis, indicating mobilization and depletion of DOC in the upper soil horizons and emphasizing the importance of shallow flowpaths during snowmelt. Cation and anion concentrations demonstrate that carbonate weathering, which dominates solute fluxes, is promoted by both sulfuric acid derived from pyritemore » oxidation in the shale bedrock and carbonic acid derived from subsurface respiration. Sulfuric acid weathering dominates during baseflow conditions when waters infiltrate below the inferred pyrite oxidation front, whereas carbonic acid weathering plays a dominant role during snowmelt as a result of shallow flowpaths. Differential C-Q relationships between solutes suggest that infiltrating waters approach calcite saturation before reaching the pyrite oxidation front, after which sulfuric acid reduces carbonate alkalinity. This reduction in alkalinity results in CO 2 outgassing when waters equilibrate to surface conditions, and reduces the riverine export of carbon and alkalinity by roughly 33% annually. In conclusion, future changes in snowmelt dynamics that control the balance of carbonic and sulfuric acid weathering may substantially alter carbon cycling in the East River. Ultimately, we demonstrate that differential C-Q relationships between major solutes can provide unique insights into the complex subsurface flow and biogeochemical dynamics that operate at catchment scales.« less

Snowmelt controls on concentration-discharge relationships and the balance of oxidative and acid-base weathering fluxes in an alpine catchment, East River, Colorado: ACID-BASE VERSUS OXIDATIVE WEATHERING FLUXES

DOE PAGES

Winnick, Matthew J.; Carroll, Rosemary W. H.; Williams, Kenneth H.; ...

2017-03-01

Although important for riverine solute and nutrient fluxes, the connections between biogeochemical processes and subsurface hydrology remain poorly characterized. We investigate these couplings in the East River, CO, a high-elevation shale-dominated catchment in the Rocky Mountains, using concentration-discharge (C-Q) relationships for major cations, anions, and organic carbon. Dissolved organic carbon (DOC) displays a positive C-Q relationship with clockwise hysteresis, indicating mobilization and depletion of DOC in the upper soil horizons and emphasizing the importance of shallow flowpaths during snowmelt. Cation and anion concentrations demonstrate that carbonate weathering, which dominates solute fluxes, is promoted by both sulfuric acid derived from pyritemore » oxidation in the shale bedrock and carbonic acid derived from subsurface respiration. Sulfuric acid weathering dominates during baseflow conditions when waters infiltrate below the inferred pyrite oxidation front, whereas carbonic acid weathering plays a dominant role during snowmelt as a result of shallow flowpaths. Differential C-Q relationships between solutes suggest that infiltrating waters approach calcite saturation before reaching the pyrite oxidation front, after which sulfuric acid reduces carbonate alkalinity. This reduction in alkalinity results in CO 2 outgassing when waters equilibrate to surface conditions, and reduces the riverine export of carbon and alkalinity by roughly 33% annually. In conclusion, future changes in snowmelt dynamics that control the balance of carbonic and sulfuric acid weathering may substantially alter carbon cycling in the East River. Ultimately, we demonstrate that differential C-Q relationships between major solutes can provide unique insights into the complex subsurface flow and biogeochemical dynamics that operate at catchment scales.« less

CALIBRATION OF SUBSURFACE BATCH AND REACTIVE-TRANSPORT MODELS INVOLVING COMPLEX BIOGEOCHEMICAL PROCESSES

EPA Science Inventory

In this study, the calibration of subsurface batch and reactive-transport models involving complex biogeochemical processes was systematically evaluated. Two hypothetical nitrate biodegradation scenarios were developed and simulated in numerical experiments to evaluate the perfor...

Multi `omics reveals role of phenotypic plasticity in governing biogeochemical hotspots within the groundwater-surface water (hyporheic) mixing zone

NASA Astrophysics Data System (ADS)

Graham, E.; Tfaily, M. M.; Crump, A.; Arntzen, E.; Romero, E. B.; Goldman, A. E.; Resch, T.; Kennedy, D.; Nelson, W. C.; Stegen, J.

2017-12-01

Subsurface groundwater-surface water mixing zones (hyporheic zones) contain spatially heterogeneous hotspots of enhanced biogeochemical activity that contribute disproportionately to river corridor function. We have a poor understanding of the processes governing hotspots, but recent advances have enabled greater mechanistic understanding. We employ a suite of ultra-high resolution measurements to investigate the mechanisms underlying biogeochemical cycles in hyporheic zone hotspots. We use Fourier transform ion cyclotron resonance mass spectrometry (FTICR-MS), metagenomic shotgun sequencing, and mass spectrometry of metaproteomes to characterize metabolite structure and metabolic transformations, microbiome structure and functional potential, and expressed microbiome functions in hyporheic sediments from the Columbia River in central Washington State. Surprisingly, microbiome structure and function in biogeochemical hotspots were indistinguishable from low-activity sediments. Metabolites were uncorrelated to protein expression but strongly related to aerobic respiration. Hotspot metabolites were distinguished by high molecular weight compounds and protein-, lignin-, and lipid-like molecules. Although the most common metabolic transformations were similar between hotspots and low-activity samples, hotspots contained a greater proportion of rare pathways, which in turn were correlated to metabolism. Our results contradicted our expectations that hotspots would be characterized by a unique microbiome with distinct physiology. Instead, our results indicate that microbial phenotypic plasticity underlies elevated hyporheic zone function, whereby the activity of rare pathways is stimulated by substrate availability. We therefore hypothesize that microbiome plasticity couples meso- (e.g., local root distribution) and macro-scale (e.g., landscape vegetation) resource heterogeneity to ecosystem-scale function. This indicates a need to mechanistically understand and represent microbiome physiological plasticity in predictive hydrobiogeochemical models that include the hyporheic zone.

Groundwater mixing at fracture intersections triggers massive iron-rich microbial mats

NASA Astrophysics Data System (ADS)

Bochet, O.; Le Borgne, T.; Bethencourt, L.; Aquilina, L.; Dufresne, A.; Pédrot, M.; Farasin, J.; Abbott, B. W.; Labasque, T.; Chatton, E.; Lavenant, N.; Petton, C.

2017-12-01

While most freshwater on Earth resides and flows in groundwater systems, these deep subsurface environments are often assumed to have little biogeochemical activity compared to surface environments. Here we report a massive microbial mat of iron-oxidizing bacteria, flourishing 60 meters below the surface, far below the mixing zone where most microbial activity is believed to occur. The abundance of microtubular structures in the mat hinted at the prevalence of of Leptothrix ochracea, but metagenomic analysis revealed a diverse consortium of iron-oxidizing bacteria dominated by unknown members of the Gallionellaceae family. This deep biogeochemical hot spot formed at the intersection of bedrock fractures, which maintain redox gradients by mixing water with different residence times and chemical compositions. Using measured fracture properties and hydrological conditions we developed a quantitative model to simulate the reactive zone where such deep hot spots could occur. While seasonal fluctuations are generally thought to decrease with depth, we found that meter-scale changes in water table level moved the depth of the reactive zone hundreds of meters because the microaerophilic threshold for ironoxidizers is highly sensitive to changes in mixing rates at fracture intersections. These results demonstrate that dynamic microbial communities can be sustained deep below the surface in bedrock fractures. Given the ubiquity of fractures at multiple scales in Earth's subsurface, such deep hot spots may strongly influence global biogeochemical cycles.

Geophysical characterisation of the groundwater-surface water interface

NASA Astrophysics Data System (ADS)

McLachlan, P. J.; Chambers, J. E.; Uhlemann, S. S.; Binley, A.

2017-11-01

Interactions between groundwater (GW) and surface water (SW) have important implications for water quantity, water quality, and ecological health. The subsurface region proximal to SW bodies, the GW-SW interface, is crucial as it actively regulates the transfer of nutrients, contaminants, and water between GW systems and SW environments. However, geological, hydrological, and biogeochemical heterogeneity in the GW-SW interface makes it difficult to characterise with direct observations. Over the past two decades geophysics has been increasingly used to characterise spatial and temporal variability throughout the GW-SW interface. Geophysics is a powerful tool in evaluating structural heterogeneity, revealing zones of GW discharge, and monitoring hydrological processes. Geophysics should be used alongside traditional hydrological and biogeochemical methods to provide additional information about the subsurface. Further integration of commonly used geophysical techniques, and adoption of emerging techniques, has the potential to improve understanding of the properties and processes of the GW-SW interface, and ultimately the implications for water quality and environmental health.

Microbial Metagenomics Reveals Climate-Relevant Subsurface Biogeochemical Processes.

PubMed

Long, Philip E; Williams, Kenneth H; Hubbard, Susan S; Banfield, Jillian F

2016-08-01

Microorganisms play key roles in terrestrial system processes, including the turnover of natural organic carbon, such as leaf litter and woody debris that accumulate in soils and subsurface sediments. What has emerged from a series of recent DNA sequencing-based studies is recognition of the enormous variety of little known and previously unknown microorganisms that mediate recycling of these vast stores of buried carbon in subsoil compartments of the terrestrial system. More importantly, the genome resolution achieved in these studies has enabled association of specific members of these microbial communities with carbon compound transformations and other linked biogeochemical processes-such as the nitrogen cycle-that can impact the quality of groundwater, surface water, and atmospheric trace gas concentrations. The emerging view also emphasizes the importance of organism interactions through exchange of metabolic byproducts (e.g., within the carbon, nitrogen, and sulfur cycles) and via symbioses since many novel organisms exhibit restricted metabolic capabilities and an associated extremely small cell size. New, genome-resolved information reshapes our view of subsurface microbial communities and provides critical new inputs for advanced reactive transport models. These inputs are needed for accurate prediction of feedbacks in watershed biogeochemical functioning and their influence on the climate via the fluxes of greenhouse gases, CO2, CH4, and N2O. Copyright © 2016 Elsevier Ltd. All rights reserved.

The impact of multi-decadal sub-surface circulation changes on sea surface chlorophyll patterns in the tropical Pacific

NASA Astrophysics Data System (ADS)

Schollaert Uz, S.; Busalacchi, A. J.; Smith, T. M.; Evans, M. N.; Brown, C.; Hackert, E. C.; Wang, X.

2016-12-01

The tropical Pacific is a region of strong forcing where physical oceanography primarily controls biological variability over the seasonal to interannual time scales observed since dedicated ocean color satellite remote sensing began in 1997. To quantify how multi-decadal, climate-scale changes impact marine biological dynamics, we used the correlation with sea-surface temperature and height to reconstruct a 50-year time series of surface chlorophyll concentrations. The reconstruction demonstrates greatest skill away from the coast and within 10o of the equator where chlorophyll variance is greatest and primarily associated with El Niño Southern Oscillation (ENSO) dynamics and secondarily associated with decadal variability. We observe significant basin-wide differences between east and central Pacific events when the El Niño events are strong: chlorophyll increases with La Niña and decreases with El Niño, with larger declines east of 180o for remotely-forced east Pacific events and west of 180o for locally-forced central Pacific events. Chlorophyll variations also reflect the physical dynamics of Pacific decadal variability with small but significant differences between cool and warm eras: consistent with advection variability west of 180o and likely driven by subsurface changes in the nutricline depth between 110-140oW. Comparisons with output from a fully-coupled biogeochemical model support the hypothesis that this anomalous region is controlled by lower frequency changes in subsurface circulation patterns that transport nutrients to the surface. Basin-wide chlorophyll distributions exhibiting spatial heterogeneity in response to multi-decadal climate forcing imply similar long-term changes in phytoplankton productivity, with implications for the marine food web and the ocean's role as a carbon sink.

Robust Representation of Integrated Surface-subsurface Hydrology at Watershed Scales

NASA Astrophysics Data System (ADS)

Painter, S. L.; Tang, G.; Collier, N.; Jan, A.; Karra, S.

2015-12-01

A representation of integrated surface-subsurface hydrology is the central component to process-rich watershed models that are emerging as alternatives to traditional reduced complexity models. These physically based systems are important for assessing potential impacts of climate change and human activities on groundwater-dependent ecosystems and water supply and quality. Integrated surface-subsurface models typically couple three-dimensional solutions for variably saturated flow in the subsurface with the kinematic- or diffusion-wave equation for surface flows. The computational scheme for coupling the surface and subsurface systems is key to the robustness, computational performance, and ease-of-implementation of the integrated system. A new, robust approach for coupling the subsurface and surface systems is developed from the assumption that the vertical gradient in head is negligible at the surface. This tight-coupling assumption allows the surface flow system to be incorporated directly into the subsurface system; effects of surface flow and surface water accumulation are represented as modifications to the subsurface flow and accumulation terms but are not triggered until the subsurface pressure reaches a threshold value corresponding to the appearance of water on the surface. The new approach has been implemented in the highly parallel PFLOTRAN (www.pflotran.org) code. Several synthetic examples and three-dimensional examples from the Walker Branch Watershed in Oak Ridge TN demonstrate the utility and robustness of the new approach using unstructured computational meshes. Representation of solute transport in the new approach is also discussed. Notice: This manuscript has been authored by UT-Battelle, LLC, under Contract No. DE-AC0500OR22725 with the U.S. Department of Energy. The United States Government retains and the publisher, by accepting the article for publication, acknowledges that the United States Government retains a non-exclusive, paid-up, irrevocable, world-wide license to publish or reproduce the published form of this manuscript, or allow others to do so, for the United States Government purposes.

Spectral induced polarization as a tool to map subsurface biogeochemical hot spots: a first laboratory evaluation in the Fe-S system

NASA Astrophysics Data System (ADS)

Nordsiek, Sven; Gilfedder, Ben; Frei, Sven

2017-04-01

Zones of intense biogeochemical reactivity (hot spots) arise in the saturated subsurface at the interface between regions with oxidizing and reducing conditions. Hot spots are both sinks and sources of different chemical compounds, thus they are of particular importance for element cycling in the subsurface. However, the investigation of hot spot structures is difficult, because they are not directly identifiable from the surface and can only be investigated by invasive methods in the subsurface. Additionally, they often form in sensitive wetland ecosystems where only non-destructive measurements are applicable to avoid significant degradation of these sensitive environments. Under these circumstances, geophysical methods may provide useful tools to identify biogeochemically active regions. One of the most important biogeochemical reactions in wetlands is the reduction of sulphate and formation and accumulation of FexSy minerals (where x and y delineate mineral stoichiometry). These reactions only occur in specific hot spots where specific chemical and microbial conditions are met. Within a research project concerning biogeochemical transformations and turnover in wetlands, we investigate the applicability of the geoelectrical method of spectral induced polarization (SIP) to locate and monitor regions containing polarizing FexSy particles as indicator for biogeochemical hot spots. After developing and testing a sample holder and a set of non-polarizing electrodes for laboratory SIP measurements, we performed experiments on natural soil samples taken from the hyporheic zone of a local river channel. The collected material originates from a location known for biogeochemical activity. The sample contains a high percentage of dark grayish/black sediment interpreted as FexSy, and possibly pyrite (FeS2). The material was homogenized and split into four samples. The FexSy concentration was adjusted to three different levels by oxidation using H2O2. For all samples we recorded the SIP spectra in the frequency range between 0.26 Hz and 25 Hz repeatedly. The sensitivity of SIP measurements to the FexSy content is promising. For a sample of the sediment saturated with its natural fluid (conductivity > 600 μS/cm), we observed phase shifts of more than 7 mrad at frequencies between 1 and 10 Hz. We expect a lower fluid conductivity and thus a larger maximum phase shift for measurements in wetlands. The application of the SIP method in field surveys at places with biogeochemical activity will be the next step in our study.

Acidification of subsurface coastal waters enhanced by eutrophication

EPA Science Inventory

Uptake of fossil-fuel carbon dioxide (CO2) from the atmosphere has acidified the surface ocean by ~0.1 pH units and driven down the carbonate saturation state. Ocean acidification is a threat to marine ecosystems and may alter key biogeochemical cycles. Coastal oceans have also b...

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

NASA Astrophysics Data System (ADS)

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 agriculture and public interests.

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

PubMed

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% (497kg) and 95% (205kg) 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 304kg (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% (240kg) 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 agriculture and public interests. Copyright © 2017. Published by Elsevier B.V.

Biogeochemical Coupling between Ocean and Sea Ice

NASA Astrophysics Data System (ADS)

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

2016-12-01

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

Coupled Spatiotemporal Dynamics of Microbial Community Ecology, Biogeochemistry, and Hydrologic Mixing

NASA Astrophysics Data System (ADS)

Stegen, J.; Johnson, T. C.; Fredrickson, J.; Wilkins, M.; Konopka, A.; Nelson, W.; Arntzen, E.; Chrisler, W.; Chu, R. K.; Fansler, S.; Kennedy, D.; Resch, T.; Tfaily, M. M.

2015-12-01

The hyporheic zone (HZ) is a critical ecosystem component that links terrestrial, surface water, and groundwater ecosystems. A dominant feature of the HZ is groundwater-surface water mixing and the input of terrestrially—as well as aquatically—derived organic carbon. In many systems the HZ has a relatively small spatial extent, but in larger riverine systems groundwater-surface water mixing can occur 100s of meters beyond the surface water shoreline; we consider these more distal locations to be within the 'subsurface interaction zone' (SIZ) as they are beyond the traditional HZ. Microbial communities in the HZ and SIZ drive biogeochemical processes in these system components, yet relatively little is known about the ecological processes that drive HZ and SIZ microbial communities. Here, we applied ecological theory, aqueous biogeochemistry, DNA sequencing, and ultra-high resolution organic carbon profiling to field samples collected through space (400m spatial extent) and time (7 month temporal extent) within the Hanford Site 300 Area. These data streams were integrated to evaluate how the influence of groundwater-surface water mixing on microbial communities changes when moving from the HZ to the broader SIZ. Our results indicate that groundwater-surface water mixing (i) consistently stimulated heterotrophic respiration, but only above a threshold of surface water intrusion, (ii) did not stimulate denitrification, (iii) caused deterministic shifts in HZ microbial communities due to changes in organic carbon composition, and (iv) did not cause shifts in SIZ microbial communities. These results suggest that microbial communities and the biogeochemical processes they drive are impacted by groundwater-surface water mixing primarily in the HZ and to a lesser extent in the SIZ.

Integrated Coupling of Surface and Subsurface Flow with HYDRUS-2D

NASA Astrophysics Data System (ADS)

Hartmann, Anne; Šimůnek, Jirka; Wöhling, Thomas; Schütze, Niels

2016-04-01

Describing interactions between surface and subsurface flow processes is important to adequately define water flow in natural systems. Since overland flow generation is highly influenced by rainfall and infiltration, both highly spatially heterogeneous processes, overland flow is unsteady and varies spatially. The prediction of overland flow needs to include an appropriate description of the interactions between the surface and subsurface flow. Coupling surface and subsurface water flow is a challenging task. Different approaches have been developed during the last few years, each having its own advantages and disadvantages. A new approach by Weill et al. (2009) to couple overland flow and subsurface flow based on a generalized Richards equation was implemented into the well-known subsurface flow model HYDRUS-2D (Šimůnek et al., 2011). This approach utilizes the one-dimensional diffusion wave equation to model overland flow. The diffusion wave model is integrated in HYDRUS-2D by replacing the terms of the Richards equation in a pre-defined runoff layer by terms defining the diffusion wave equation. Using this approach, pressure and flux continuity along the interface between both flow domains is provided. This direct coupling approach provides a strong coupling of both systems based on the definition of a single global system matrix to numerically solve the coupled flow problem. The advantage of the direct coupling approach, compared to the loosely coupled approach, is supposed to be a higher robustness, when many convergence problems can be avoided (Takizawa et al., 2014). The HYDRUS-2D implementation was verified using a) different test cases, including a direct comparison with the results of Weill et al. (2009), b) an analytical solution of the kinematic wave equation, and c) the results of a benchmark test of Maxwell et al. (2014), that included several known coupled surface subsurface flow models. Additionally, a sensitivity analysis evaluating the effects of various model parameters on simulated overland flow (while considering or neglecting the effects of subsurface flow) was carried out to verify the applicability of the model to different problems. The model produced reasonable results in describing the diffusion wave approximation and its interactions with subsurface flow processes. The model could handle coupled surface-subsurface processes for conditions involving runoff generated by infiltration excess, saturation excess, or run-on, as well as a combination of these runoff generating processes. Several standard features of the HYDRUS 2D model, such as root water uptake and evaporation from the soil surface, as well as evaporation from runoff layer, can still be considered by the new model. The code required relatively small time steps when overland flow was active, resulting in long simulation times, and sometimes produced poor mass balance. The model nevertheless showed potential to be a useful tool for addressing various issues related to irrigation research and to natural generation of overland flow at the hillslope scale. Maxwell, R., Putti, M., Meyerhoff, S., Delf, J., Ferguson, I., Ivanov, V., Kim, J., Kolditz, O., Kollet, S., Kumar, M., Lopez, S., Niu, J., Paniconi, C., Park, Y.-J., Phanikumar, M., Shen, C., Sudicky, E., and Sulis, M. (2014). Surface-subsurface model intercomparison: A first set of benchmark results to diagnose integrated hydrology and feedbacks. Water Resourc. Res., 50:1531-1549. Šimůnek, J., van Genuchten, M. T., and Šejna, M. (2011). The HYDRUS Software Package for Simulating Two- and Three-Dimensional Movement of Water, Heat, and Multiple Solutes in Variably-Saturated Media. Technical Manual, Version 2.0, PC Progress, Prague, Czech Republic. Takizawa, K., Bazilevs Y., Tezduyar, T. E., Long, C.C., Marsden, A. L. and Schjodt.K., Patient-Specific Cardiovascular Fluid Mechanics Analysis with the ST and ALE-VMS Method in Idelsohn, S. R. (2014). Numerical Simulations of Coupled Problems in Engineering. Springer. Weill, S., Mouche, E., and Patin, J. (2009). A generalized Richards equation for surface/subsurface flow modelling. Journal of Hydrology, 366:9-20.

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

NASA Astrophysics Data System (ADS)

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

2014-12-01

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

Linking sediment structure, hydrological functioning and biogeochemical cycling in disturbed coastal saltmarshes and implications for vegetation development

NASA Astrophysics Data System (ADS)

Spencer, Kate; Harvey, Gemma; James, Tempest; Simon, Carr; Michelle, Morris

2014-05-01

Saltmarsh restoration undoubtedly provides environmental enhancement, with vegetation quickly re-establishing following the breach of sea walls and subsequent tidal inundation of previously defended areas. Yet evidence increasingly suggests that the restored saltmarshes do not have the same biological characteristics as their natural counterparts (Mossman et al. 2012) and this may be in part be due to physicochemical parameters at the site including anoxia and poor drainage. Hence, restored saltmarshes may not offer the range and quality of ecosystem services anticipated. These environments will have been 'disturbed' by previous land use and there is little understanding of the impacts of this disturbance on the wider hydrogeomorphic and biogeochemical functioning in restored saltmarshes and the implications for saltmarsh vegetation development. This study examines linkages between physical sediment characteristics, sediment structure (using X-ray microtomography), sub-surface hydrology (using pressure transducers and time series analysis), and sediment and porewater geochemistry (major and trace elements, major anions) in sediment cores collected from undisturbed saltmarshes and those restored by de-embankment. Sub-surface sediments in restored saltmarshes have lower organic matter content, lower moisture content and higher bulk density than undisturbed sites. Using X-ray tomography a clear horizon can be observed which separates relict agricultural soils at depth with less dense and structureless sediments deposited since de-embankment. Ratios of open to closed pore space suggest that while undisturbed saltmarshes have the highest porosity, restored saltmarshes have larger void spaces, but limited pore connectivity. Sub-surface hydrological response to tidal flooding was subdued in the restored compared to the undisturbed site, suggesting that porewater flow may be impeded. Time series analysis indicated that flow pathways differ in restored saltmarsh sediments with preferential horizontal flows. The undisturbed saltmarsh displayed typical vertical geochemical sediment profiles. However, in the restored sites total Fe and Mn are elevated at depth indicating an absence of diagenetic cycling, whilst porewater sulphate and nitrate increased at depth suggesting that vertical solute transport is impeded in restored sites. In surface sediments, though total Hg concentrations are similar, Hg methylation rates are significantly higher than in the undisturbed saltmarsh suggesting that surface anoxia and poor drainage may result in increased mobilization and bioavailability of Hg. These findings have implications for the wider biogeochemical ecosystem services offered by saltmarsh restoration and the water-logged, anoxic conditions produced are unsuitable for seedling germination and plant growth. This highlights the need for integrated understanding of physical and biogeochemical processes.

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

NASA Astrophysics Data System (ADS)

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.

Role of Subsurface Physics in the Assimilation of Surface Soil Moisture Observations

NASA Technical Reports Server (NTRS)

Reichle, R. H.

2010-01-01

Root zone soil moisture controls the land-atmosphere exchange of water and energy and exhibits memory that may be useful for climate prediction at monthly scales. Assimilation of satellite-based surface soil moisture observations into a land surface model is an effective way to estimate large-scale root zone soil moisture. The propagation of surface information into deeper soil layers depends on the model-specific representation of subsurface physics that is used in the assimilation system. In a suite of experiments we assimilate synthetic surface soil moisture observations into four different models (Catchment, Mosaic, Noah and CLM) using the Ensemble Kalman Filter. We demonstrate that identical twin experiments significantly overestimate the information that can be obtained from the assimilation of surface soil moisture observations. The second key result indicates that the potential of surface soil moisture assimilation to improve root zone information is higher when the surface to root zone coupling is stronger. Our experiments also suggest that (faced with unknown true subsurface physics) overestimating surface to root zone coupling in the assimilation system provides more robust skill improvements in the root zone compared with underestimating the coupling. When CLM is excluded from the analysis, the skill improvements from using models with different vertical coupling strengths are comparable for different subsurface truths. Finally, the skill improvements through assimilation were found to be sensitive to the regional climate and soil types.

Opportunities and challenges in studies of deep life (Invited)

NASA Astrophysics Data System (ADS)

Edwards, K. J.

2010-12-01

Over the past two decades, there has been an increasing awareness within the geological, microbiological, and oceanographic communities of the potentially vast microbial biosphere that is harbored beneath the surface of the Earth. With this awareness has come a mounting effort to study this potential biome - to better quantify biomass abundance, activity, and biogeochemical activity. In the Earth system, the largest deep subsurface biome is also the least accessible - the deep ocean subsurface biosphere. The oceanic deep biosphere also has greatest potential for influencing global scale biogeochemical processes -the carbon and energy cycles for example, and other elemental cycles. To address these topics and mount interdisciplinary efforts to study the deep subsurface marine biosphere, we have recently formed a center in support integrative, collaborative investigations. The national science foundation Center for Dark Biosphere Investigations (C-DEBI), has been initiated for the explicit purpose of resolving the extent, function, dynamics and implications of the subseafloor biosphere. This talk will discuss C-DEBI science, with focus on some of the opportunities and challenges in the study of deep life in the ocean, and the role that C-DEBI will play in meeting them

Reactive Transport Modeling of Induced Calcite Precipitation Reaction Fronts in Porous Media Using A Parallel, Fully Coupled, Fully Implicit Approach

NASA Astrophysics Data System (ADS)

Guo, L.; Huang, H.; Gaston, D.; Redden, G. D.; Fox, D. T.; Fujita, Y.

2010-12-01

Inducing mineral precipitation in the subsurface is one potential strategy for immobilizing trace metal and radionuclide contaminants. Generating mineral precipitates in situ can be achieved by manipulating chemical conditions, typically through injection or in situ generation of reactants. How these reactants transport, mix and react within the medium controls the spatial distribution and composition of the resulting mineral phases. Multiple processes, including fluid flow, dispersive/diffusive transport of reactants, biogeochemical reactions and changes in porosity-permeability, are tightly coupled over a number of scales. Numerical modeling can be used to investigate the nonlinear coupling effects of these processes which are quite challenging to explore experimentally. Many subsurface reactive transport simulators employ a de-coupled or operator-splitting approach where transport equations and batch chemistry reactions are solved sequentially. However, such an approach has limited applicability for biogeochemical systems with fast kinetics and strong coupling between chemical reactions and medium properties. A massively parallel, fully coupled, fully implicit Reactive Transport simulator (referred to as “RAT”) based on a parallel multi-physics object-oriented simulation framework (MOOSE) has been developed at the Idaho National Laboratory. Within this simulator, systems of transport and reaction equations can be solved simultaneously in a fully coupled, fully implicit manner using the Jacobian Free Newton-Krylov (JFNK) method with additional advanced computing capabilities such as (1) physics-based preconditioning for solution convergence acceleration, (2) massively parallel computing and scalability, and (3) adaptive mesh refinements for 2D and 3D structured and unstructured mesh. The simulator was first tested against analytical solutions, then applied to simulating induced calcium carbonate mineral precipitation in 1D columns and 2D flow cells as analogs to homogeneous and heterogeneous porous media, respectively. In 1D columns, calcium carbonate mineral precipitation was driven by urea hydrolysis catalyzed by urease enzyme, and in 2D flow cells, calcium carbonate mineral forming reactants were injected sequentially, forming migrating reaction fronts that are typically highly nonuniform. The RAT simulation results for the spatial and temporal distributions of precipitates, reaction rates and major species in the system, and also for changes in porosity and permeability, were compared to both laboratory experimental data and computational results obtained using other reactive transport simulators. The comparisons demonstrate the ability of RAT to simulate complex nonlinear systems and the advantages of fully coupled approaches, over de-coupled methods, for accurate simulation of complex, dynamic processes such as engineered mineral precipitation in subsurface environments.

Hyporheic flow and transport processes: mechanisms, models, and biogeochemical implications

USGS Publications Warehouse

Boano, Fulvio; Harvey, Judson W.; Marion, Andrea; Packman, Aaron I.; Revelli, Roberto; Ridolfi, Luca; Anders, Wörman

2014-01-01

Fifty years of hyporheic zone research have shown the important role played by the hyporheic zone as an interface between groundwater and surface waters. However, it is only in the last two decades that what began as an empirical science has become a mechanistic science devoted to modeling studies of the complex fluid dynamical and biogeochemical mechanisms occurring in the hyporheic zone. These efforts have led to the picture of surface-subsurface water interactions as regulators of the form and function of fluvial ecosystems. Rather than being isolated systems, surface water bodies continuously interact with the subsurface. Exploration of hyporheic zone processes has led to a new appreciation of their wide reaching consequences for water quality and stream ecology. Modern research aims toward a unified approach, in which processes occurring in the hyporheic zone are key elements for the appreciation, management, and restoration of the whole river environment. In this unifying context, this review summarizes results from modeling studies and field observations about flow and transport processes in the hyporheic zone and describes the theories proposed in hydrology and fluid dynamics developed to quantitatively model and predict the hyporheic transport of water, heat, and dissolved and suspended compounds from sediment grain scale up to the watershed scale. The implications of these processes for stream biogeochemistry and ecology are also discussed."

Hyporheic flow and transport processes: Mechanisms, models, and biogeochemical implications

NASA Astrophysics Data System (ADS)

Boano, F.; Harvey, J. W.; Marion, A.; Packman, A. I.; Revelli, R.; Ridolfi, L.; Wörman, A.

2014-12-01

Fifty years of hyporheic zone research have shown the important role played by the hyporheic zone as an interface between groundwater and surface waters. However, it is only in the last two decades that what began as an empirical science has become a mechanistic science devoted to modeling studies of the complex fluid dynamical and biogeochemical mechanisms occurring in the hyporheic zone. These efforts have led to the picture of surface-subsurface water interactions as regulators of the form and function of fluvial ecosystems. Rather than being isolated systems, surface water bodies continuously interact with the subsurface. Exploration of hyporheic zone processes has led to a new appreciation of their wide reaching consequences for water quality and stream ecology. Modern research aims toward a unified approach, in which processes occurring in the hyporheic zone are key elements for the appreciation, management, and restoration of the whole river environment. In this unifying context, this review summarizes results from modeling studies and field observations about flow and transport processes in the hyporheic zone and describes the theories proposed in hydrology and fluid dynamics developed to quantitatively model and predict the hyporheic transport of water, heat, and dissolved and suspended compounds from sediment grain scale up to the watershed scale. The implications of these processes for stream biogeochemistry and ecology are also discussed.

Genome-to-Watershed Predictive Understanding of Terrestrial Environments

NASA Astrophysics Data System (ADS)

Hubbard, S. S.; Agarwal, D.; Banfield, J. F.; Beller, H. R.; Brodie, E.; Long, P.; Nico, P. S.; Steefel, C. I.; Tokunaga, T. K.; Williams, K. H.

2014-12-01

Although terrestrial environments play a critical role in cycling water, greenhouse gasses, and other life-critical elements, the complexity of interactions among component microbes, plants, minerals, migrating fluids and dissolved constituents hinders predictive understanding of system behavior. The 'Sustainable Systems 2.0' project is developing genome-to-watershed scale predictive capabilities to quantify how the microbiome affects biogeochemical watershed functioning, how watershed-scale hydro-biogeochemical processes affect microbial functioning, and how these interactions co-evolve with climate and land-use changes. Development of such predictive capabilities is critical for guiding the optimal management of water resources, contaminant remediation, carbon stabilization, and agricultural sustainability - now and with global change. Initial investigations are focused on floodplains in the Colorado River Basin, and include iterative model development, experiments and observations with an early emphasis on subsurface aspects. Field experiments include local-scale experiments at Rifle CO to quantify spatiotemporal metabolic and geochemical responses to O2and nitrate amendments as well as floodplain-scale monitoring to quantify genomic and biogeochemical response to natural hydrological perturbations. Information obtained from such experiments are represented within GEWaSC, a Genome-Enabled Watershed Simulation Capability, which is being developed to allow mechanistic interrogation of how genomic information stored in a subsurface microbiome affects biogeochemical cycling. This presentation will describe the genome-to-watershed scale approach as well as early highlights associated with the project. Highlights include: first insights into the diversity of the subsurface microbiome and metabolic roles of organisms involved in subsurface nitrogen, sulfur and hydrogen and carbon cycling; the extreme variability of subsurface DOC and hydrological controls on carbon and nitrogen cycling; geophysical identification of floodplain hotspots that are useful for model parameterization; and GEWaSC demonstration of how incorporation of identified microbial metabolic processes improves prediction of the larger system biogeochemical behavior.

Active Serpentinization and the Potential for a Diverse Subsurface Biosphere

NASA Astrophysics Data System (ADS)

Canovas, P. A.; Shock, E.

2013-12-01

The ubiquitous nature of serpentinization and the unique fluids it generates have major consequences for habitat generation, abiotic organic synthesis, and biosynthesis. The production of hydrogen from the anaerobic hydrolysis of ultramafic minerals sets the redox state of serpentinizing fluids to be thermodynamically favorable for these processes. Consequently, a host of specialized microbial populations and metabolisms can be sustained. Active low-temperature serpentinizing systems, such as the Samail ophiolite in Oman, offer an ideal opportunity to investigate biogeochemical processes during the alteration of ultramafic minerals. At the Samail ophiolite in particular, serpentinization may provide the potential for an active subsurface microbial community shielded from potentially unfavorable surface conditions. Support for this assertion comes from geochemical data including Mg, Ca, CH4 (aq), and H2 (aq) abundances indicating that methane is a product of serpentinization. To further investigate viable metabolic strategies, affinity calculations were performed on both the surface waters and the hyperalkaline springs, which may be considered as messengers of processes occurring in the subsurface. Almost all sites yield positive affinities (i.e., are thermodynamically favorable) for a diverse suite of serpentinization metabolisms including methanogenesis, anammox, and carbon monoxide, nitrate, and sulfate reduction with hydrogen, as well as anaerobic methanotrophy coupled to nitrate, nitrite, and sulfate reduction. Reaction path modeling was performed to ascertain the extent to which serpentinization and mixing of surface waters with hyperalkaline spring waters in the subsurface can generate suitable habitats. The serpentinization model simulates the reaction of pristine Oman harzburgite with surface water to quantify the redox state and generation of hyperalkaline spring water. Preliminary results show that water-rock ratios as high as 100 could effectively reduce the system and create a thermodynamic drive sufficient to convert all of the dissolved inorganic carbon into methane. This indicates that the system is poised to create the reducing conditions necessary to support a subsurface biosphere very early in the serpentinizing process, and that the subsurface biosphere could extend upwards to very near the surface. The mixing model simulates the percolation of surface water into the active serpentinization zone. During the mixing process, methane is calculated to be more stable than carbonate species until approximately 100g of surface water have been added to 1 kg of the serpentinizing fluid. These results suggest that unreacted surface water flowing directly into the serpentinizing zone can create the disequilibria necessary for methanogenesis, and possibly other metabolisms, to proceed while still maintaining the low redox state of the system. As long as the recharge to the hyperalkaline reservoir does not exceed ten percent of the reservoir, methanogenesis and other serpentinization metabolisms can thrive off the disequilibria generated through mixing.

First-order exchange coefficient coupling for simulating surface water-groundwater interactions: Parameter sensitivity and consistency with a physics-based approach

USGS Publications Warehouse

Ebel, B.A.; Mirus, B.B.; Heppner, C.S.; VanderKwaak, J.E.; Loague, K.

2009-01-01

Distributed hydrologic models capable of simulating fully-coupled surface water and groundwater flow are increasingly used to examine problems in the hydrologic sciences. Several techniques are currently available to couple the surface and subsurface; the two most frequently employed approaches are first-order exchange coefficients (a.k.a., the surface conductance method) and enforced continuity of pressure and flux at the surface-subsurface boundary condition. The effort reported here examines the parameter sensitivity of simulated hydrologic response for the first-order exchange coefficients at a well-characterized field site using the fully coupled Integrated Hydrology Model (InHM). This investigation demonstrates that the first-order exchange coefficients can be selected such that the simulated hydrologic response is insensitive to the parameter choice, while simulation time is considerably reduced. Alternatively, the ability to choose a first-order exchange coefficient that intentionally decouples the surface and subsurface facilitates concept-development simulations to examine real-world situations where the surface-subsurface exchange is impaired. While the parameters comprising the first-order exchange coefficient cannot be directly estimated or measured, the insensitivity of the simulated flow system to these parameters (when chosen appropriately) combined with the ability to mimic actual physical processes suggests that the first-order exchange coefficient approach can be consistent with a physics-based framework. Copyright ?? 2009 John Wiley & Sons, Ltd.

Performance Analysis and Scaling Behavior of the Terrestrial Systems Modeling Platform TerrSysMP in Large-Scale Supercomputing Environments

NASA Astrophysics Data System (ADS)

Kollet, S. J.; Goergen, K.; Gasper, F.; Shresta, P.; Sulis, M.; Rihani, J.; Simmer, C.; Vereecken, H.

2013-12-01

In studies of the terrestrial hydrologic, energy and biogeochemical cycles, integrated multi-physics simulation platforms take a central role in characterizing non-linear interactions, variances and uncertainties of system states and fluxes in reciprocity with observations. Recently developed integrated simulation platforms attempt to honor the complexity of the terrestrial system across multiple time and space scales from the deeper subsurface including groundwater dynamics into the atmosphere. Technically, this requires the coupling of atmospheric, land surface, and subsurface-surface flow models in supercomputing environments, while ensuring a high-degree of efficiency in the utilization of e.g., standard Linux clusters and massively parallel resources. A systematic performance analysis including profiling and tracing in such an application is crucial in the understanding of the runtime behavior, to identify optimum model settings, and is an efficient way to distinguish potential parallel deficiencies. On sophisticated leadership-class supercomputers, such as the 28-rack 5.9 petaFLOP IBM Blue Gene/Q 'JUQUEEN' of the Jülich Supercomputing Centre (JSC), this is a challenging task, but even more so important, when complex coupled component models are to be analysed. Here we want to present our experience from coupling, application tuning (e.g. 5-times speedup through compiler optimizations), parallel scaling and performance monitoring of the parallel Terrestrial Systems Modeling Platform TerrSysMP. The modeling platform consists of the weather prediction system COSMO of the German Weather Service; the Community Land Model, CLM of NCAR; and the variably saturated surface-subsurface flow code ParFlow. The model system relies on the Multiple Program Multiple Data (MPMD) execution model where the external Ocean-Atmosphere-Sea-Ice-Soil coupler (OASIS3) links the component models. TerrSysMP has been instrumented with the performance analysis tool Scalasca and analyzed on JUQUEEN with processor counts on the order of 10,000. The instrumentation is used in weak and strong scaling studies with real data cases and hypothetical idealized numerical experiments for detailed profiling and tracing analysis. The profiling is not only useful in identifying wait states that are due to the MPMD execution model, but also in fine-tuning resource allocation to the component models in search of the most suitable load balancing. This is especially necessary, as with numerical experiments that cover multiple (high resolution) spatial scales, the time stepping, coupling frequencies, and communication overheads are constantly shifting, which makes it necessary to re-determine the model setup with each new experimental design.

Identifying biogeochemical processes beneath stormwater infiltration ponds in support of a new best management practice for groundwater protection

USGS Publications Warehouse

O'Reilly, Andrew M.; Chang, Ni-Bin; Wanielista, Martin P.; Xuan, Zhemin; Schirmer, Mario; Hoehn, Eduard; Vogt, Tobias

2011-01-01

 When applying a stormwater infiltration pond best management practice (BMP) for protecting the quality of underlying groundwater, a common constituent of concern is nitrate. Two stormwater infiltration ponds, the SO and HT ponds, in central Florida, USA, were monitored. A temporal succession of biogeochemical processes was identified beneath the SO pond, including oxygen reduction, denitrification, manganese and iron reduction, and methanogenesis. In contrast, aerobic conditions persisted beneath the HT pond, resulting in nitrate leaching into groundwater. Biogeochemical differences likely are related to soil textural and hydraulic properties that control surface/subsurface oxygen exchange. A new infiltration BMP was developed and a full-scale application was implemented for the HT pond. Preliminary results indicate reductions in nitrate concentration exceeding 50% in soil water and shallow groundwater beneath the HT pond.

Modeling Global Ocean Biogeochemistry With Physical Data Assimilation: A Pragmatic Solution to the Equatorial Instability

NASA Astrophysics Data System (ADS)

Park, Jong-Yeon; Stock, Charles A.; Yang, Xiaosong; Dunne, John P.; Rosati, Anthony; John, Jasmin; Zhang, Shaoqing

2018-03-01

Reliable estimates of historical and current biogeochemistry are essential for understanding past ecosystem variability and predicting future changes. Efforts to translate improved physical ocean state estimates into improved biogeochemical estimates, however, are hindered by high biogeochemical sensitivity to transient momentum imbalances that arise during physical data assimilation. Most notably, the breakdown of geostrophic constraints on data assimilation in equatorial regions can lead to spurious upwelling, resulting in excessive equatorial productivity and biogeochemical fluxes. This hampers efforts to understand and predict the biogeochemical consequences of El Niño and La Niña. We develop a strategy to robustly integrate an ocean biogeochemical model with an ensemble coupled-climate data assimilation system used for seasonal to decadal global climate prediction. Addressing spurious vertical velocities requires two steps. First, we find that tightening constraints on atmospheric data assimilation maintains a better equatorial wind stress and pressure gradient balance. This reduces spurious vertical velocities, but those remaining still produce substantial biogeochemical biases. The remainder is addressed by imposing stricter fidelity to model dynamics over data constraints near the equator. We determine an optimal choice of model-data weights that removed spurious biogeochemical signals while benefitting from off-equatorial constraints that still substantially improve equatorial physical ocean simulations. Compared to the unconstrained control run, the optimally constrained model reduces equatorial biogeochemical biases and markedly improves the equatorial subsurface nitrate concentrations and hypoxic area. The pragmatic approach described herein offers a means of advancing earth system prediction in parallel with continued data assimilation advances aimed at fully considering equatorial data constraints.

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

NASA Astrophysics Data System (ADS)

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 new characterization and inversion approaches to quantify co-variability between above and below ground dynamics and the susceptibility of vegetation to drought stress. More information is provided at: Watershed.lbl.gov

A New Approach to Predict Microbial Community Assembly and Function Using a Stochastic, Genome-Enabled Modeling Framework

NASA Astrophysics Data System (ADS)

King, E.; Brodie, E.; Anantharaman, K.; Karaoz, U.; Bouskill, N.; Banfield, J. F.; Steefel, C. I.; Molins, S.

2016-12-01

Characterizing and predicting the microbial and chemical compositions of subsurface aquatic systems necessitates an understanding of the metabolism and physiology of organisms that are often uncultured or studied under conditions not relevant for one's environment of interest. Cultivation-independent approaches are therefore important and have greatly enhanced our ability to characterize functional microbial diversity. The capability to reconstruct genomes representing thousands of populations from microbial communities using metagenomic techniques provides a foundation for development of predictive models for community structure and function. Here, we discuss a genome-informed stochastic trait-based model incorporated into a reactive transport framework to represent the activities of coupled guilds of hypothetical microorganisms. Metabolic pathways for each microbe within a functional guild are parameterized from metagenomic data with a unique combination of traits governing organism fitness under dynamic environmental conditions. We simulate the thermodynamics of coupled electron donor and acceptor reactions to predict the energy available for cellular maintenance, respiration, biomass development, and enzyme production. While `omics analyses can now characterize the metabolic potential of microbial communities, it is functionally redundant as well as computationally prohibitive to explicitly include the thousands of recovered organisms into biogeochemical models. However, one can derive potential metabolic pathways from genomes along with trait-linkages to build probability distributions of traits. These distributions are used to assemble groups of microbes that couple one or more of these pathways. From the initial ensemble of microbes, only a subset will persist based on the interaction of their physiological and metabolic traits with environmental conditions, competing organisms, etc. Here, we analyze the predicted niches of these hypothetical microbes and assess the ability of a stochastically assembled community of organisms to predict subsurface biogeochemical dynamics.

Linking Chaotic Advection with Subsurface Biogeochemical Processes

NASA Astrophysics Data System (ADS)

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

2017-12-01

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

Drogue performance evaluation. Part 1: Data acquisition. [for accurate description of subsurface water circulation

NASA Technical Reports Server (NTRS)

Johnson, R. E.

1975-01-01

Surface buoy/subsurface drogue drag coupling was investigated. Data acquisition methods and techniques derived from several experimental cruises on the Chesapeake Bay are presented. Four buoys were utilized: three coupled to steel plates rigidly attached to each other at right angles and at various depths; and one spar type that did not require drag plates. Data from these surface floats and the drogue depth combinations were processed. Errors in tracking the surface buoys are discussed.

Hierarchical Bayesian method for mapping biogeochemical hot spots using induced polarization imaging

DOE PAGES

Wainwright, Haruko M.; Flores Orozco, Adrian; Bucker, Matthias; ...

2016-01-29

In floodplain environments, a naturally reduced zone (NRZ) is considered to be a common biogeochemical hot spot, having distinct microbial and geochemical characteristics. Although important for understanding their role in mediating floodplain biogeochemical processes, mapping the subsurface distribution of NRZs over the dimensions of a floodplain is challenging, as conventional wellbore data are typically spatially limited and the distribution of NRZs is heterogeneous. In this work, we present an innovative methodology for the probabilistic mapping of NRZs within a three-dimensional (3-D) subsurface domain using induced polarization imaging, which is a noninvasive geophysical technique. Measurements consist of surface geophysical surveys andmore » drilling-recovered sediments at the U.S. Department of Energy field site near Rifle, CO (USA). Inversion of surface time domain-induced polarization (TDIP) data yielded 3-D images of the complex electrical resistivity, in terms of magnitude and phase, which are associated with mineral precipitation and other lithological properties. By extracting the TDIP data values colocated with wellbore lithological logs, we found that the NRZs have a different distribution of resistivity and polarization from the other aquifer sediments. To estimate the spatial distribution of NRZs, we developed a Bayesian hierarchical model to integrate the geophysical and wellbore data. In addition, the resistivity images were used to estimate hydrostratigraphic interfaces under the floodplain. Validation results showed that the integration of electrical imaging and wellbore data using a Bayesian hierarchical model was capable of mapping spatially heterogeneous interfaces and NRZ distributions thereby providing a minimally invasive means to parameterize a hydrobiogeochemical model of the floodplain.« less

Effect of Freeze-Thaw Cycles on Soil Nitrogen Reactive Transport in a Polygonal Arctic Tundra Ecosystem at Barrow AK Using 3-D Coupled ALM-PFLOTRAN

NASA Astrophysics Data System (ADS)

Yuan, F.; Wang, G.; Painter, S. L.; Tang, G.; Xu, X.; Kumar, J.; Bisht, G.; Hammond, G. E.; Mills, R. T.; Thornton, P. E.; Wullschleger, S. D.

2017-12-01

In Arctic tundra ecosystem soil freezing-thawing is one of dominant physical processes through which biogeochemical (e.g., carbon and nitrogen) cycles are tightly coupled. Besides hydraulic transport, freezing-thawing can cause pore water movement and aqueous species gradients, which are additional mechanisms for soil nitrogen (N) reactive-transport in Tundra ecosystem. In this study, we have fully coupled an in-development ESM(i.e., Advanced Climate Model for Energy, ACME)'s Land Model (ALM) aboveground processes with a state-of-the-art massively parallel 3-D subsurface thermal-hydrology and reactive transport code, PFLOTRAN. The resulting coupled ALM-PFLOTRAN model is a Land Surface Model (LSM) capable of resolving 3-D soil thermal-hydrological-biogeochemical cycles. This specific version of PFLOTRAN has incorporated CLM-CN Converging Trophic Cascade (CTC) model and a full and simple but robust soil N cycle. It includes absorption-desorption for soil NH4+ and gas dissolving-degasing process as well. It also implements thermal-hydrology mode codes with three newly-modified freezing-thawing algorithms which can greatly improve computing performance in regarding to numerical stiffness at freezing-point. Here we tested the model in fully 3-D coupled mode at the Next Generation Ecosystem Experiment-Arctic (NGEE-Arctic) field intensive study site at the Barrow Environmental Observatory (BEO), AK. The simulations show that: (1) synchronous coupling of soil thermal-hydrology and biogeochemistry in 3-D can greatly impact ecosystem dynamics across polygonal tundra landscape; and (2) freezing-thawing cycles can add more complexity to the system, resulting in greater mobility of soil N vertically and laterally, depending upon local micro-topography. As a preliminary experiment, the model is also implemented for Pan-Arctic region in 1-D column mode (i.e. no lateral connection), showing significant differences compared to stand-alone ALM. The developed ALM-PFLOTRAN coupling codes embeded within ESM will be used for Pan-Arctic regional evaluation of climate change-caused ecosystem responses and their feedbacks to climate system at various scales.

Long-Term Cultivation and Metagenomics Reveal Ecophysiology of Previously Uncultivated Thermophiles Involved in Biogeochemical Nitrogen Cycle

PubMed Central

Kato, Shingo; Sakai, Sanae; Hirai, Miho; Tasumi, Eiji; Nishizawa, Manabu; Suzuki, Katsuhiko; Takai, Ken

2018-01-01

Many thermophiles thriving in a natural high-temperature environment remain uncultivated, and their ecophysiological functions in the biogeochemical cycle remain unclear. In the present study, we performed long-term continuous cultivation at 65°C and 70°C using a microbial mat sample, collected from a subsurface geothermal stream, as the inoculum, and reconstructed the whole genome of the maintained populations using metagenomics. Some metagenome-assembled genomes (MAGs), affiliated into phylum-level bacterial and archaeal clades without cultivated representatives, contained genes involved in nitrogen metabolism including nitrification and denitrification. Our results show genetic components and their potential interactions for the biogeochemical nitrogen cycle in a subsurface geothermal environment. PMID:29459499

Optical method and apparatus for detection of surface and near-subsurface defects in dense ceramics

DOEpatents

Ellingson, William A.; Brada, Mark P.

1995-01-01

A laser is used in a non-destructive manner to detect surface and near-subsurface defects in dense ceramics and particularly in ceramic bodies with complex shapes such as ceramic bearings, turbine blades, races, and the like. The laser's wavelength is selected based upon the composition of the ceramic sample and the laser can be directed on the sample while the sample is static or in dynamic rotate or translate motion. Light is scattered off surface and subsurface defects using a preselected polarization. The change in polarization angle is used to select the depth and characteristics of surface/subsurface defects. The scattered light is detected by an optical train consisting of a charge coupled device (CCD), or vidicon, television camera which, in turn, is coupled to a video monitor and a computer for digitizing the image. An analyzing polarizer in the optical train allows scattered light at a given polarization angle to be observed for enhancing sensitivity to either surface or near-subsurface defects. Application of digital image processing allows subtraction of digitized images in near real-time providing enhanced sensitivity to subsurface defects. Storing known "feature masks" of identified defects in the computer and comparing the detected scatter pattern (Fourier images) with the stored feature masks allows for automatic classification of detected defects.

Wellhead with non-ferromagnetic materials

DOEpatents

Hinson, Richard A [Houston, TX; Vinegar, Harold J [Bellaire, TX

2009-05-19

Wellheads for coupling to a heater located in a wellbore in a subsurface formation are described herein. At least one wellhead may include a heater located in a wellbore in a subsurface formation; and a wellhead coupled to the heater. The wellhead may be configured to electrically couple the heater to one or more surface electrical components. The wellhead may include at least one non-ferromagnetic material such that ferromagnetic effects are inhibited in the wellhead. Systems and methods for using such wellheads for treating a subsurface formation are described herein.

Coupled charge migration and fluid mixing in reactive fronts

NASA Astrophysics Data System (ADS)

Ghosh, Uddipta; Bandopadhyay, Aditya; Jougnot, Damien; Le Borgne, Tanguy; Meheust, Yves

2017-04-01

Quantifying fluid mixing in subsurface environments and its consequence on biogeochemical reactions is of paramount importance owing to its role in processes such as contaminant migration, aquifer remediation, CO2 sequestration or clogging processes, to name a few (Dentz et al. 2011). The presence of strong velocity gradients in porous media is expected to lead to enhanced diffusive mixing and augmented reaction rates (Le Borgne et al. 2014). Accurate in situ imaging of subsurface reactive solute transport and mixing remains to date a challenging proposition: the opacity of the medium prevents optical imaging and field methods based on tracer tests do not provide spatial information. Recently developed geophysical methods based on the temporal monitoring of electrical conductivity and polarization have shown promises for mapping and monitoring biogeochemical reactions in the subsurface although it remains challenging to decipher the multiple sources of electrical signals (e.g. Knight et al. 2010). In this work, we explore the coupling between fluid mixing, reaction and charge migration in porous media to evaluate the potential of mapping reaction rates from electrical measurements. To this end, we develop a new theoretical framework based on a lamellar mixing model (Le Borgne et al. 2013) to quantify changes in electrical mobility induced by chemical reactions across mixing fronts. Electrical conductivity and induced polarization are strongly dependent on the concentration of ionic species, which in turn depend on the local reaction rates. Hence, our results suggest that variation in real and complex electrical conductivity may be quantitatively related to the mixing and reaction dynamics. Thus, the presented theory provides a novel upscaling framework for quantifying the coupling between mixing, reaction and charge migration in heterogeneous porous media flows. References: Dentz. et al., Mixing, spreading and reaction in heterogeneous media: A brief review J. Contam. Hydrol. 120-121, 1 (2011). Le Borgne et al. Impact of Fluid Deformation on Mixing-Induced Chemical Reactions in heterogeneous Flows, Geophys. Res. Lett. 41, 7898 (2014). Knight, et al., Geophysics at the interface: Response of geophysical properties to solid-fluid, fluid-fluid, and solid-solid interfaces. Rev. Geophys. 48, (2010). Le Borgne et al. (2013) Stretching, coalescence and mixing in porous media, Phys. Rev. Lett., 110, 204501

Subsurface flow pathway dynamics in the active layer of coupled permafrost-hydrogeological systems under seasonal and annual temperature variability.

NASA Astrophysics Data System (ADS)

Frampton, Andrew

2017-04-01

There is a need for improved understanding of the mechanisms controlling subsurface solute transport in the active layer in order to better understand permafrost-hydrological-carbon feedbacks, in particular with regards to how dissolved carbon is transported in coupled surface and subsurface terrestrial arctic water systems under climate change. Studying solute transport in arctic systems is also relevant in the context of anthropogenic pollution which may increase due to increased activity in cold region environments. In this contribution subsurface solute transport subject to ground surface warming causing permafrost thaw and active layer change is studied using a physically based model of coupled cryotic and hydrogeological flow processes combined with a particle tracking method. Changes in subsurface water flows and solute transport travel times are analysed for different modelled geological configurations during a 100-year warming period. Results show that for all simulated cases, the minimum and mean travel times increase non-linearly with warming irrespective of geological configuration and heterogeneity structure. The timing of the start of increase in travel time depends on heterogeneity structure, combined with the rate of permafrost degradation that also depends on material thermal and hydrogeological properties. These travel time changes are shown to depend on combined warming effects of increase in pathway length due to deepening of the active layer, reduced transport velocities due to a shift from horizontal saturated groundwater flow near the surface to vertical water percolation deeper into the subsurface, and pathway length increase and temporary immobilization caused by cryosuction-induced seasonal freeze cycles. The impact these change mechanisms have on solute and dissolved substance transport is further analysed by integrating pathway analysis with a Lagrangian approach, incorporating considerations for both dissolved organic and inorganic carbon releases. Further model development challenges are also highlighted and discussed, including coupling between subsurface and surface runoff, soil deformations, as well as site applications and larger system scales.

AquaDiva: Understanding the Link between the Surface and Subsurface Biogeosphere

NASA Astrophysics Data System (ADS)

Trumbore, S.; Küsel, K.; Totsche, K. U.; Schwab, V.; Herrmann, M.; Nowak, M. E.; Gleixner, G.

2017-12-01

In the collaborative research project AquaDiva, we combine hydrogeochemical, metagemonic and biogeochemical tools to understand how the complex interactions between geologic setting and surface land use influence the function and biodiversity of the subsurface, especially ground water ecosystems. At the Hainich Critical Zone Exploratory in central Germany, we investigate soil and seepage waters in recharge areas and aquifers in a fractured limestone setting characterized by a dynamic water infiltration regime. Within the Exploratory, we have so far identified three distinct biogeochemical zones in which land use and lithologic differences combine to give rise to surprisingly different biotic communities and hydrogeochemical properties with different degrees of connection to the surface. Here we will focus on how we have combined carbon isotopic, organic biomarkers such as phospholipid fatty acids, and `omics' approaches to determine (i) how deep soil-borne microorganisms can be traced into the subsurface, and (ii) which energy sources sustain microbial life in oligotrophic limestone aquifers. With increasing travel distance to the surface, there is a decline in the abundance of microbes, with less than 5% of the taxa identified overlapping with those identified in the soils. Dissolved organic matter also is altered as it passes through soils, demonstrating an overall increase in molecular weight and a change in molecular makeup as well as radiocarbon content. Using the radiocarbon signature (corrected for the influence of carbonate dissolution), as a way to identify if organic C is being supplied recently by plants, atmosphere sources of energy for communities within the aquifer differ for the identified biogeochemical zones and include (i) dominance of inputs of fresh organic carbon from the surface feeding heterotrophy in oxygenic environments; (ii) CO2 fixation linked to nitrogen and sulfur cycling in anoxic environments and (iii) rock-derived organic matter that enters the food web and supplies up to a quarter of the carbon recycled by heterotrophic groundwater microbial communities. The degree of importance of these metabolisms, and therefore their connection to the surface is different in the three zones, with the importance of internal cycling and rock-derived organic matter as a C source surprisingly higher in anoxic aquifers. The second major goal of AquaDiva is to determine the impact of individual precipitation events on the groundwater. Initial results show dramatic effects of singular precipitation events, particularly of events during snowmelt season when the major recharge takes place. Observed large fluctuations in the water table have been linked to changes in the vertical distribution of oxygen and events of carbonate dissolution/precipitation can be identified based on variations of CO2 and O2. Ongoing research in AquaDiva seeks to better quantify the rates of groundwater movement and mixing as a way of further understanding how the emergent subsurface properties can be better linked to processes of transport and transformation within the Critical Zone.

From terrestrial to aquatic fluxes: Integrating stream dynamics within a dynamic global vegetation modeling framework

NASA Astrophysics Data System (ADS)

Hoy, Jerad; Poulter, Benjamin; Emmett, Kristen; Cross, Molly; Al-Chokhachy, Robert; Maneta, Marco

2016-04-01

Integrated terrestrial ecosystem models simulate the dynamics and feedbacks between climate, vegetation, disturbance, and hydrology and are used to better understand biogeography and biogeochemical cycles. Extending dynamic vegetation models to the aquatic interface requires coupling surface and sub-surface runoff to catchment routing schemes and has the potential to enhance how researchers and managers investigate how changes in the environment might impact the availability of water resources for human and natural systems. In an effort towards creating such a coupled model, we developed catchment-based hydrologic routing and stream temperature model to pair with LPJ-GUESS, a dynamic global vegetation model. LPJ-GUESS simulates detailed stand-level vegetation dynamics such as growth, carbon allocation, and mortality, as well as various physical and hydrologic processes such as canopy interception and through-fall, and can be applied at small spatial scales, i.e., 1 km. We demonstrate how the coupled model can be used to investigate the effects of transient vegetation dynamics and CO2 on seasonal and annual stream discharge and temperature regimes. As a direct management application, we extend the modeling framework to predict habitat suitability for fish habitat within the Greater Yellowstone Ecosystem, a 200,000 km2 region that provides critical habitat for a range of aquatic species. The model is used to evaluate, quantitatively, the effects of management practices aimed to enhance hydrologic resilience to climate change, and benefits for water storage and fish habitat in the coming century.

GEOCHEMISTRY OF SUBSURFACE REACTIVE BARRIERS FOR REMEDIATION OF CONTAMINATED GROUND WATER

EPA Science Inventory

Reactive barriers that couple subsurface fluid flow with a passive chemical treatment zone are emerging, cost effective approaches for in-situ remediation of contaminated groundwater. Factors such as the build-up of surface precipitates, bio-fouling, and changes in subsurface tr...

Biogeochemical carbon coupling influences global precipitation in geoengineering experiments

NASA Astrophysics Data System (ADS)

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

2013-02-01