Engineering Papers⌕ Search

Engineering topics

Murray, Christopher J.

Publications and source records attributed to Murray, Christopher J..

11 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

A novel construct for scaling groundwater–river interactions based on machine-guided hydromorphic classification

Hydrologic exchange between river channels and adjacent subsurface environments is a key process that influences water quality and ecosystem function in river corridors. Predictive numerical models are needed to understand responses of river corridors to environmental change and to support sustainable watershed management. We posit that systematic hydromorphic classification provides a scaling construct that facilitates extrapolation of outputs from local-scale mechanistic models to reduced-order models applicable at reach and watershed scales. This in turn offers the potential to improve large-scale predictions of river corridor hydrobiogeochemical processes. Here we present a new machine-guided hydromorphic classification methodology that addresses the key requirements of this objective, and we demonstrate its application to a segment of the Columbia River in the northwestern United States. The resulting hydromorphic classes form spatially coherent and physically interpretable hydromorphic units that exhibit distinct behaviors in terms of distributions of subsurface residence times (a primary control on critical biogeochemical reactions). This approach forms the basis of ongoing research that is evaluating the formulation of reduced-order models and transferability of results to other river reaches and larger scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Kilometer-Scale Hydrologic Exchange Flows in a Gravel Bed River Corridor and Their Implications to Solute Migration

In this work, a well-characterized field site along a major river corridor was used to investigate the dynamic pathways and impacts of subsurface hydrogeologic structures on kilometer-scale hydrologic exchange flows between river water and groundwater. A aqueous uranium (U aq ) plume exists at the site that discharges to the Columbia River. We performed temporally intensive monitoring of specific conductance (SpC) and U aq concentrations within the plume at increasing distances from the river shoreline and at locations both within and outside a presumed subsurface pathway of lateral hydrologic exchange. Specific conductance and U aq were utilized as in-situ tracers of hydrologic exchange and associated groundwater-surface water mixing. Monitoring data was collected for an approximate two year period during which seasonal river elevations varied by more than 2 m causing distinct events of river water intrusion and retreat from the near-shore aquifer, and large temporal variations in well-water analyte concentrations. Statistical clustering analysis was applied to identify well groupings exhibiting common time series behaviors. PFLOTRAN was implemented within a modeling domain of 3.3 x 1.9 km with dynamic flow boundary conditions and accurate subsurface topography for the basal aquitard underlying the aquifer. In-silico calculations of river water tracer behavior were performed to assess predominant flow directions and velocities influenced by river stage variations and variable aquitard topography, and their collective influences on SpC and U aq concentrations. We demonstrate complicated multidirectional flow behaviors at the kilometer scale resulting from hydrologic exchange and subsurface structure that influence plume dynamics and determine the overall system behavior.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Identification and mapping of riverbed sediment facies in the Columbia River through integration of field observations and numerical simulations

In the Hanford Reach of the Columbia River, a thin layer of recent alluvium overlies the sedimentary formations that comprise the unconfined groundwater aquifer. Experimental and modeling studies have demonstrated that this alluvial layer exerts significant control on the exchange of groundwater and surface water (hydrologic exchange flux), and is associated with elevated levels of biogeochemical activity. This layer is also observed to be strongly heterogeneous, and quantifying the spatial distribution of properties over the range of scales of interest is challenging. Facies are elements of a sediment classification scheme that groups complex geologic materials into a set of discrete classes according to distinguishing features. Facies classifications have been used as a framework for assigning heterogeneous material properties to grid cells of numerical models of flow and reactive transport in subsurface media. The usefulness of such an approach hinges on being able to relate facies to quantitative properties needed for flow and reactive transport modeling, and on being able to map facies over the domain of interest using readily available information. While aquifer facies have been used in various modeling contexts, application of this concept to riverbed sediments is relatively new. Here we describe an approach for categorizing and mapping recent alluvial (riverbed) sediments based on the integration of diverse observations with numerical simulations of river hydrodynamics. The facies have distinct distributions of sediment texture which correspond to variations in hydraulic properties, and therefore provide a useful framework for assigning heterogeneous properties in numerical simulations of hydrologic exchange flows and biogeochemical processes.

Hou, Zhangshuan↗
DOE OSTI2018

Riverbed Hydrologic Exchange Dynamics in a Large Regulated River Reach

Hydrologic exchange flux (HEF) is an important hydrologic component in river corridors that includes both bidirectional (hyporheic) and unidirectional (gaining/losing) surface water – groundwater exchanges. Quantifying HEF rates in a large regulated river is difficult due to the large spatial domains, complexity of geomorphologic features and subsurface properties, and the great stage variations created by dam operations at multiple time scales. In this study, we developed a method that combined numerical modeling and field measurements for estimating HEF rates across the river bed in a 7-km long reach of the highly regulated Columbia River. A high-resolution computational fluid dynamics (CFD) modeling framework was developed and validated by field measurements and other modeling results to characterize the HEF dynamics across the river bed. We found that about 85% of the time from 2008-2014 the river was losing water with an annual average net HEF rates across the river bed (Qz) of -2.3 m 3 s -1 (negative indicating downwelling). June was the only month that the river gained water, with monthly averaged Qz of 0.8 m 3 s -1 . We also found that the daily dam operations increased the hourly gross gaining and losing rate over an average year of 8% and 2%, respectively. By investigating the HEF feedbacks at various time scales, we suggest that the dam operations could reduce the HEF at seasonal time scale by decreasing the seasonal flow variations, while also enhance the HEF at sub-daily time scale by generating high frequency discharge variations. These changes could generate significant impacts on biogeochemical processes in the hyporheic zone.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Geochemical and Microbial Community Attributes in Relation to Hyporheic Zone Geological Facies

The hyporheic zone (HZ) is the active ecotone between the surface stream and groundwater, where exchanges of water, nutrients, and organic matter occur in response to variations in discharge and riverbed properties. Within this region, a confluence of surface-derived organic carbon and subsurface nitrogen (in the form of nitrate) has been shown to stimulate microbial activity and transformations of carbon and nitrogen species. For example, production of gases such as CO 2 , N 2 and N 2 O indicate hyporheic zones might have a significant effect on energy and nutrient flows between the atmosphere and the subsurface. Managed and seasonal river stage changes and geomorphology-controlled sediment texture drive water flow within the HZ of the Columbia River. To examine the relationship between sediment texture, biogeochemistry, and biological activity in the HZ, the grain size distributions for sediment samples taken across 320 m of shoreline were characterized to define geological facies, and the relationships among physical properties of the facies, physicochemical attributes of the local environment, and the structure and activity of associated microbial communities were examined. Mud and sand content and the presence of carbon and nitrogen oxidizers were found to explain the variability in many biogeochemical attributes. Microbial community analysis revealed a high relative abundance of putative ammonia-oxidizing Thaumarchaeota and nitrite-oxidizing Nitrospirae, together comprising ~20% of the total community across all samples, but scant ammonia-oxidizing Bacteria. Network analysis of operational taxonomic units and the measured geophysical, chemical, and functional parameters showed negative relationships between abundance-based modules of organisms and sand and mud contents, and positive relationships with total organic carbon. The relationships identified in this work indicate grain size distribution is a good predictor of biogeochemical properties, and that subsets of the overall microbial community respond to different sediment texture. Some member populations of these sub-communities appear to respond directly to environmental conditions, while others may be dependent on the function of the first group. For example, nitrification is a strong primary response to the observed conditions, and this activity appears to support a larger heterotrophic community. Relationships between facies and hydrobiogeochemical properties enables facies-based conditional simulation/mapping of these properties to inform multiscale modeling of hyporheic exchange and biogeochemical processes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Coupling among Microbial Communities, Biogeochemistry, and Mineralogy across Biogeochemical Facies

Physical properties of sediments are commonly used to define subsurface lithofacies and these same physical properties influence subsurface microbial communities. This suggests an (unexploited) opportunity to use the spatial distribution of facies to predict spatial variation in biogeochemically relevant microbial attributes. Here, we characterize three biogeochemical facies—oxidized, reduced, and transition—within one lithofacies and elucidate relationships among facies features and microbial community biomass, diversity, and community composition. Consistent with previous observations of biogeochemical hotspots at environmental transition zones, we find elevated biomass within a biogeochemical facies that occurred at the transition between oxidized and reduced biogeochemical facies. Microbial diversity—the number of microbial taxa—was lower within the reduced facies and was well-explained by a combination of pH and mineralogy. Null modeling revealed that microbial community composition was influenced by ecological selection imposed by redox state and mineralogy, possibly due to effects on nutrient availability or transport. As an illustrative case, we predict microbial biomass concentration across a three-dimensional spatial domain by coupling the spatial distribution of subsurface biogeochemical facies with biomass-facies relationships revealed here. We expect that merging such an approach with hydro-biogeochemical models will provide important constraints on simulated dynamics, thereby reducing uncertainty in model predictions.

Stegen, James C.↗
DOE OSTI2015

Bicarbonate impact on U(VI) bioreduction in a shallow alluvial aquifer

Field-scale biostimulation and desorption tracer experiments conducted in a uranium (U) contaminated, shallow alluvial aquifer have provided insight into the coupling of microbiology, biogeochemistry, and hydrogeology that control U mobility in the subsurface. Initial experiments successfully tested the concept that Fe-reducing bacteria such as Geobacter sp. could enzymatically reduce soluble U(VI) to insoluble U(IV) during in situ electron donor amendment (Anderson et al., 2003; Williams et al., 2011). In parallel, in situ desorption tracer tests using bicarbonate amendment demonstrated rate-limited U(VI) desorption (Fox et al., 2012). These results and prior laboratory studies underscored the importance of enzymatic U(VI)-reduction and suggested the ability to combine desorption and bioreduction of U(VI). Here in this paper we report the results of a new field experiment in which bicarbonate-promoted uranium desorption and acetate amendment were combined and compared to an acetate amendment-only experiment in the same experimental plot. Results confirm that bicarbonate amendment to alluvial aquifer sediments desorbs U(VI) and increases the abundance of Ca-uranyl-carbonato complexes. At the same time, the rate of acetate-promoted enzymatic U(VI) reduction was greater in the presence of added bicarbonate in spite of the increased dominance of Ca-uranyl-carbonato aqueous complexes. A model-simulated peak rate of U(VI) reduction was ~3.8 times higher during acetate-bicarbonate treatment than under acetate-only conditions. Lack of consistent differences in microbial community structure between acetate-bicarbonate and acetate-only treatments suggest that a significantly higher rate of U(VI) reduction in the bicarbonate-impacted sediment may be due to a higher intrinsic rate of microbial reduction induced by elevated concentrations of the bicarbonate oxyanion. The findings indicate that bicarbonate amendment may be useful in improving the engineered bioremediation of uranium in aquifers.

58 GEOSCIENCES↗
DOE OSTI2014

A high-performance workflow system for subsurface simulation

Subsurface modeling applications typically neglect uncertainty in the conceptual models, past or future scenarios, and attribute most or all uncertainty to errors in model parameters. In this contribution, uncertainty in technetium-99 transport in a heterogeneous, deep vadose zone is explored with respect to the conceptual model using a next generation user environment called Akuna. Akuna provides a range of tools to manage environmental modeling projects, from managing simulation data to visualizing results from high-performance computational simulators. Core toolsets accessible through the user interface include model setup, grid generation, parameter estimation, and uncertainty quantification. The BC Cribs site at Hanford in southeastern Washington State is used to demonstrate Akuna capabilities. At the BC Cribs site, conceptualization of the system is highly uncertain because only sparse information is available for the geologic conceptual model, the physical and chemical properties of the sediments, and the history of waste disposal operations. Using the Akuna toolset to perform an analysis of conservative solute transport, significant prediction uncertainty in simulated concentrations is demonstrated by conceptual model variation. This demonstrates that conceptual model uncertainty is an important consideration in sparse data environments such as BC Cribs. It is also demonstrated that Akuna and the underlying toolset provides an integrated modeling environment that streamlines model setup, parameter optimization, and uncertainty analyses for high-performance computing applications.

Akuna↗