Search NASA⌕ Search

SEARCH · Search NASA

Results for “subsurface flow”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 721 records · Page 40

Understanding coastal wetland hydrology with a new regional‐scale, process‐based hydrological model

Abstract Coastal wetlands represent an ecotone between ocean and terrestrial ecosystems, providing important services, including flood mitigation, fresh water supply, erosion control, carbon sequestration, and wildlife habitat. The environmental setting of a wetland and the hydrological connectivity between a wetland and adjacent terrestrial and aquatic systems together determine wetland hydrology. Yet little is known about regional‐scale hydrological interactions among uplands, coastal wetlands, and coastal processes, such as tides, sea level rise, and saltwater intrusion, which together control the dynamics of wetland hydrology. This study presents a new regional‐scale, physically based, distributed wetland hydrological model, PIHM‐Wetland, which integrates the surface and subsurface hydrology with coastal processes and accounts for the influence of wetland inundation on energy budgets and evapotranspiration (ET). The model was validated using in situ hydro‐meteorological measurements and Moderate Resolution Imaging Spectroradiometer (MODIS) ET data for a forested and herbaceous wetland in North Carolina, USA, which confirmed that the model accurately represents the major wetland hydrological behaviours. Modelling results indicate that topographic gradient is a primary control of groundwater flow direction in adjacent uplands. However, seasonal climate patterns become the dominant control of groundwater flow at lower coastal plain and land–ocean interface. We found that coastal processes largely influence groundwater table (GWT) dynamics in the coastal zone, 300 to 800 m from the coastline in our study area. Among all the coastal processes, tides are the dominant control on GWT variation. Because of inundation, forested and herbaceous wetlands absorb an additional 6% and 10%, respectively, of shortwave radiation annually, resulting in a significant increase in ET. Inundation alters ET partitioning through canopy evaporation, transpiration, and soil evaporation, the effect of which is stronger in cool seasons than in warm seasons. The PIHM‐Wetland model provides a new tool that improves the understanding of wetland hydrological processes on a regional scale. Insights from this modelling study provide benchmarks for future research on the effects of sea level rise and climate change on coastal wetland functions and services.

Zhang, Yu↗

Hyporheic Exchange in Sand Dunes Under a Freely Deforming River Water Surface

Bedform-driven hyporheic exchange is conditioned by the head gradients at the sediment-water interface. Local exchange phenomena between the surface and the subsurface are often driven by the dynamic forces related to the velocity distributions around sediment waves. In open channels, the static forces are represented by the water depth, yet most computational fluid dynamics models use a rigid-lid approximation. Here, we investigated whether and when the deformation of the river’s free-surface influences bedform-driven hyporheic exchange. This was done through simulations of coupled open channel and hyporheic flows with the air-water interface modeled either as a free-surface or a rigid-lid across increasing subcritical Froude numbers from ~0.05 to ~0.95. The normalized hyporheic flux was higher when considering free-surface deformation across most of the range of Froude numbers considered. When the Froude number was larger than 0.6 and smaller than 0.85, both hydrostatic and nonhydrostatic-driven fluxes increased significantly compared to the rigid-lid approach and the total flux was about five times that predicted with the rigid-lid. These results indicate that bedform-driven hyporheic fluxes are underestimated in most studies that use a rigid-lid assumption in subcritical flows.

58 GEOSCIENCES↗

Uncertainties in Local Intense Precipitation Flood Modeling

In the Task 1 report for the U.S. Nuclear Regulatory Commission (NRC) funded Local Intense Precipitation (LIP) Probabilistic Flood Hazard Assessment (PFHA) Pilot Study, the Pacific Northwest National Laboratory (PNNL) reviewed hydrologic and hydraulic modeling approaches for estimation of LIP flooding and their implementation in readily available simulation software packages (Prasad and Yuan 2020). NRC’s interest is in developing guidance for performing LIP PFHA at nuclear power plant (NPP) sites. This report describes the work done under Task 2 of the study that focused on describing uncertainties associated with LIP flood modeling. Two main types of uncertainties exist in LIP flood modeling: (1) aleatory uncertainties that arise from the inherent natural variability of the hydrometeorological system and (2) epistemic uncertainties that arise from the analysts’ incomplete knowledge of the hydrometeorological and site configuration. Principal sources for aleatory variability in LIP flood simulations at NPP sites include precipitation (i.e., magnitude, duration, and temporal distribution of LIP events), initial conditions (e.g., soil moisture content, stormwater drainage discharge, surface storage), boundary conditions (e.g., upstream discharge, downstream water surface elevations), and long-term temporal trends (e.g., climate change). This report describes approaches used for estimating the aleatory variability in precipitation including precipitation-frequency analysis, numerical weather prediction, and stochastic weather generation. Data sources that can be used to estimate aleatory variability in initial and boundary conditions are also described. Approaches to include effects of long-term trends like those from climate change into estimation of aleatory variability are summarized. Sources of epistemic uncertainty in LIP flood simulations at NPP sites include process representation (e.g., multiple approaches to represent runoff generation, stormwater drainage, and hydraulic routing), site configurations (e.g., site layout, flow features, status of temporary flood protection, blockage of drains), model resolution, and long-term temporal trends (e.g., known/planned site alterations, land-use changes at and in the vicinity of the site). The report describes alternative process representations (methods and models), particularly those implemented in the LIP flood simulation software packages reviewed in Task 1 report, and lists the associated model parameters. Approaches for estimating model parameters when surface and subsurface water exchanges occur are also described.

42 ENGINEERING↗

Seasonal and Interannual Patterns and Controls of Hydrological Fluxes in an Amazon Floodplain Lake With a Surface-Subsurface Process Model

Floodplain lakes represent important aquatic ecosystems, and field-based estimates of their water budgets are difficult to obtain, especially over multiple years. We examine the hydrological fluxes for an Amazon floodplain lake connected to the Solimões River using a process-based hydrologic model. Water exchanges between the river and lake agree well with field estimates, including the timing of different hydrological phases. However, beyond available field data, modeling results show that the seven simulated years all differed from each other. These interannual differences were caused by the interplay between phases when water levels were rising with river-water flowing into the lake (RWRI), versus rising with lake-water flowing out to the river (RWLO). This exchange determines the river-water content in the lake ( C L ). Maximum C L occurred before river levels peaked because local catchment contributions can be sufficient to push lake-water out to the river, even as river levels rise. Numerical experiments show that the seasonal distribution of local rainfall, local catchment size, and interannual variability in both climate and river stage can contribute to differing dynamics of C L in a floodplain lake. Their impacts vary among phases: river-rise dominates the RWRI, whereas local hydrological processes dominate the RWLO and receding-water phases. Intermediate-to-long-term rainfall accumulation controls C L during the RWLO phase, whereas annual precipitation accumulation is important for C L during low water. Finally, our model generalizes beyond limited available field studies and offers potential to better understand floodplain lakes in other areas and how regional versus local changes in climate may affect their hydrological dynamics.

54 ENVIRONMENTAL SCIENCES↗

Ecosystems and Networks Integrated with Genes and Molecular Assemblies (Final Scientific/Technical Report)

The ENIGMA science focus area (SFA) is a multi-disciplinary, multi-institutional research effort focused on addressing foundational knowledge gaps in environmental microbial communities by studying groundwater and sediment microbiomes in the shallow subsurface at the contaminated Oak Ridge Field Research Site (FRC). We seek to discover and characterize the reciprocal interactions between the microbial communities and the geochemical and geophysical parameters of the shallow subsurface within the contamination plume. The goal of the Fields lab is to help manage and conduct laboratory experiments informed by field observations using ecological and physiological approaches to study microbial populations in situ as well as ex situ (in the laboratory). In particular, we aim to study and characterize the impact of hydrological constraints on free-living and biofilm biomass and activity with increasing spatial and temporal resolution under static and flow conditions. We also aim to help characterize novel microbial groups that are present and active under relevant field conditions, including the development of molecular techniques for in situ detection as well as metabolic interactions that underlie pertinent physiology and ecology.

59 BASIC BIOLOGICAL SCIENCES↗

A multirate mass transfer model to represent the interaction of multicomponent biogeochemical processes between surface water and hyporheic zones (SWAT-MRMT-R 1.0)

Abstract. Surface water quality along river corridors can be modulated by hyporheic zones (HZs) that are ubiquitous and biogeochemically active. Watershed management practices often ignore the potentially important role of HZs as a natural reactor. To investigate the effect of hydrological exchange and biogeochemical processes on the fate of nutrients in surface water and HZs, a novel model, SWAT-MRMT-R, was developed coupling the Soil and Water Assessment Tool (SWAT) watershed model and the reaction module from a flow and reactive transport code (PFLOTRAN). SWAT-MRMT-R simulates concurrent nonlinear multicomponent biogeochemical reactions in both the channel water and its surrounding HZs, connecting the channel water and HZs through hyporheic exchanges using multirate mass transfer (MRMT) representation. Within the model, HZs are conceptualized as transient storage zones with distinguished exchange rates and residence times. The biogeochemical processes within HZs are different from those in the channel water. Hyporheic exchanges are modeled as multiple first-order mass transfers between the channel water and HZs. As a numerical example, SWAT-MRMT-R is applied to the Hanford Reach of the Columbia River, a large river in the United States, focusing on nitrate dynamics in the channel water. Major nitrate contaminants entering the Hanford Reach include those from the legacy waste, irrigation return flows (irrigation water that is not consumed by crops and runs off as point sources to the stream), and groundwater seepage resulting from irrigated agriculture. A two-step reaction sequence for denitrification and an aerobic respiration reaction is assumed to represent the biogeochemical transformations taking place within the HZs. The spatially variable hyporheic exchange rates and residence times in this example are estimated with the basin-scale Networks with EXchange and Subsurface Storage (NEXSS) model. Our simulation results show that (1), given a residence time distribution, how the exchange fluxes to HZs are approximated when using MRMT can significantly change the amount of nitrate consumption in HZs through denitrification and (2) source locations of nitrate have a different impact on surface water quality due to the spatially variable hyporheic exchanges.

58 GEOSCIENCES↗

Unclassified Source Term and Radionuclide Data for the Groundwater Flow and Contaminant Transport Model for Corrective Action Units 101 and 102: Central and Western Pahute Mesa, Nye County, Nevada, Revision 1

This report documents the evaluation of information and data available on the unclassified source term and radionuclide (RN) contamination for Central and Western Pahute Mesa (PM): Corrective Action Units (CAUs) 101 and 102. The Central and Western PM CAUs constitute two areas of the Nevada National Security Site (NNSS) used for underground nuclear testing (Figure 1-1). Eighty-five underground nuclear tests were conducted at the PM CAUs between 1965 and 1992 (NNSA/NFO, 2015b). Underground nuclear tests deposit radioactive elements in the subsurface, resulting in groundwater contamination in the vicinity of the underground test locations. At the PM CAUs, the testing resulted in more than 2.68 × 10 7 curies (Ci) of radioactivity (decay-corrected to September 30, 2012) remaining underground (Finnegan et al., 2016). Of the 85 underground nuclear tests conducted at the PM CAUs, three were shallow tests detonated within 200 meters (m) of the surface and referred to as the “near-surface tests”; the remaining 82 were detonated deeper than 200 m below ground surface (bgs) and are referred to as the “deep underground tests” in this document.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

The Influence of Aqueous Dispersivity on Transport of Radiological Material at LANL MDA G - 20450

Aqueous dispersivity is a phenomenon by which a solute moving in groundwater mixes with water that does not contain the solute and becomes diluted. Solutes take paths of varying lengths through pore networks or experience slower velocities due to interactions with pore walls. Aqueous dispersivity leads to breakthrough curves with earlier arrival than expected in purely advective systems and persistent 'tails,' or late-arriving, dilute concentrations of the solute. Dispersivity can be longitudinal (in the direction of flow) or transverse (orthogonal to flow). LANL has disposed of radioactive material at Material Disposal Area (MDA) G since the late 1950's. Most of the disposed materials at MDA G are low-level waste (LLW), but other types of radioactive wastes are disposed or stored at the site. For example, there is a collection of remote-handled transuranic waste (RH-TRU) waste in 33 shafts, which were intended to be stored retrievably. Aqueous-phase radionuclides included in models of contaminant transport at MDA G are subject to dispersive processes in porous media of the vadose zone and saturated zone. To model dispersivity at MDA G probabilistically, we created site-specific distributions of aqueous dispersivity. First, we performed a comprehensive review of peer-reviewed publications, reports from other performance assessments of MDA G, and documentation from sites that were considered similar based on expert opinion, such as Yucca Mountain. Data were included from a wide range of sources, so that distributions developed based on the data reflect the current state of knowledge in properties of the subsurface environment. Data were grouped based on value type (e.g., determined via experiment or used in a different model), rock type, zone type (saturated or unsaturated), and site, and were weighted in the distributions according to their relevance to the site and zone type. We combined Monte Carlo simulation and bootstrapping techniques to produce approximate distributions of mean aqueous dispersivity. Distributions of dispersivity for saturated and unsaturated zone transport are differentiated by data weight and the length scales over which transport occurs. We find that underestimating or disregarding transverse dispersivity may lead to overstating simulated well capture in the saturated zone at MDA G. Neglecting transverse dispersivity is a common assumption due in part to use of one-dimensional transport models that do not consider solute spreading orthogonal to flow. Also, we find that variability in longitudinal dispersivity is associated with substantial travel time alterations in the MDA G unsaturated zone and in simulated contaminant breakthrough times at a downstream observation plane in the MDA G saturated zone. Finally, at large scales (>1000 m), a linear trend between the dispersivity coefficient and length over which transport occurs may overestimate the dispersivity coefficient. Future modeling work with application to solute spreading could include interacting advective and diffusive zones in the highly fractured basalt below MDA G. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Martian Plains Volcanism in Syria Planum and Tempe Mareotis as Analogs to the Eastern Snake River Plains, Idaho: Similarities and Possible Petrologic Contributions to Topography

Prior to the Mars Global Surveyor (MGS) and Mars Odyssey (MO) missions, The Syria Planum region of Mars was noted for several clusters of small (5-100 km) shield volcanoes and collapse craters, long tube and fissure-fed lava flows, and possible volcanic vents that were thought to be nearly contemporaneous with the volcanism in the Tempe- Mareotis province, which has long been known for volcanic shields and vents analogous to those of the Eastern Snake River Plains (ESRP) in Idaho. Recent MGS-based work on regional and global populations of martian small shields has revealed significant global trends in edifice attributes that are well-explained by eruption models with latitudinal variations in subsurface water/ice abundance, consistent with recent MO evidence for significant amounts of subsurface water that varies in latitude abundance s, and topographic and morphologic evidence for more geologically recent lava-ice relationships. However, while the global trends in small volcano data can be at least partially explained by volatile interactions with volcanism, some global and regional characteristics appear to be perhaps better explained by possible compositional, crystallinity or eruption style variations. This study expands the sampling of shields done in martian initial global studies for the Syria Planum and Tempe-Mareotis regions, which display a newly visible breadth and number of features in image and topography data. We compare these features to a similar range of features visible in the ESRP where both compositional and eruption style variations can quantitatively be shown to contribute to morphologic and topographic differences.

Sakimoto, S. E. H.↗

Approaches for the Simulation of Coupled Processes in Evolving Fractured Porous Media Enabled by Exascale Computing

Models have historically represented fractured porous media with continuum descriptions that characterize the media using bulk parameters. The impact of small-scale features is not captured in these models, although they may be controlling the performance of subsurface applications. Pore-scale models can simulate processes in small-scale features by representing the pore space geometry explicitly but are computationally expensive for large domains. The alternative multiscale approach entails the combination of pore-scale and continuum-scale descriptions in a single framework. We use Chombo-Crunch, a computational capability that discretizes complex geometries with an adaptive, embedded boundary method to contrast these two approaches. Chombo-Crunch takes advantage of recent computational performance and memory bandwidth improvements resulting from the emergence of exascale computing resources. These combined improvements enable the efficient simulation of reactive transport in fractured media with a high degree of fidelity and the ability to capture the control small-scale processes exert on the overall medium evolution.

42 ENGINEERING↗

Pump-and-Treat Impacts on Aquifer Conditions and the Development of System-Wide Performance Strategies - 20377

The pump-and-treat (P and T) system on the Hanford Site's Central Plateau currently treats multiple contaminants within the 200-ZP-1 operable unit, including nitrate through its biological treatment processes. Although it was initially designed to support the 200-ZP-1 Remedial Action Objectives, the P and T system has evolved in use to address contaminants from other regions of the aquifer as well. The P and T influents are biogeochemically distinct, treating groundwaters with elevated concentrations of cyanide, uranium, technetium-99, and nitrate. Given the varying geochemical signatures from current and potential future sources, operational conditions, and amendments considered for treatment processes, there is a need to identify the impacts that the P and T effluent has on the aquifer, and how the aquifer affects P and T performance. To this end, laboratory efforts have focused on characterizing P and T effluent and its biogeochemical effects on the sediments through a series repacked column experiments. The P and T effluent is rich in particulate matter. A systematic inspection of this particulate matter in the P and T effluent showed significant levels of manganese/iron oxides, calcite crystals, relatively large organic matter molecules, and biomass. Therefore, the column experiments were designed and conducted for four different treatment systems: artificial groundwater, filtered P and T effluent, non-filtered P and T effluent, and chlorinated P and T effluent. Different levels of microbiologically- and geochemically-induced changes in the columns occurred and modified the flow systems. For example, physical filtration occurred in the sediments where the particulate material had accumulated in the top portion of the column. Reactive transport modeling of the laboratory experiments is being carried out so that the reaction network can be extrapolated to the field scale. These simulations will provide predictive evaluations of system performance based on P and T operational changes and aquifer mass loading and capacity estimates. This information can be used to address combined subsurface and operational challenges, such as the need for frequent well redevelopment to maintain injection capacity. This information also provides the technical basis for decisions related to the P and T system operations that support remedy optimization efforts and short- and long-term remedy decisions related to multiple operable units on the Hanford Site Central Plateau. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Corrosion of Al 2 O 3 /Cr and Ti 2 O 3 /Cr composites in flowing air and CO 2 at 750°C

Hot-pressed Al 2 O 3 /Cr and Ti 2 O 3 /Cr composites were exposed to synthetic air or CO 2 at 750 °C. The Al 2 O 3 /Cr composites formed external Cr 2 O 3 scales, and exhibited slow parabolic oxidation kinetics, in both gases. The Cr 2 O 3 scales formed in CO 2 possessed finer grains, and thickened at a faster rate, than Cr 2 O 3 scales formed in air. The Ti 2 O 3 /Cr composites formed thick external TiO 2 scales with internal oxidation zones that also contained Cr 2 O 3 scales on subsurface Cr grains. After 20 h at 750 °C, the Ti 2 O 3 /Cr composites exhibited parabolic oxidation kinetics in both gases, although the oxidation rates were significantly higher than for the Al 2 O 3 /Cr composites.

36 MATERIALS SCIENCE↗

Paragenetic controls on CO 2 -fluid-rock interaction and weakening in a macroporous-dominated sandstone

The injection and storage of anthropogenic CO 2 in the subsurface is being deployed as a climate change mitigation tool; however, diagenetic-paragenetic heterogeneity in sandstone reservoirs often contributes to interval specific chemomechanical changes that affect injection and can increase leakage risk.In this paper we address reservoir heterogeneities’ impact on chemomechanical changes in a macroporous-dominated lithofacies of Morrow B sandstone, a formation containing several diagenetically-distinct hydraulic facies while undergoing enhanced oil recovery (EOR) and carbon dioxide (CO 2 ) sequestration. We performed three flow-through experiments using a CO 2 -charged or uncharged formation water combined with four indirect tensile strength tests per post-test sample. We then used the microstructure and paragenetic sequence to understand chemomechanical weakening with key observations as follows: dissolution of carbonates and feldspars changed porosity; increased permeability led to reclassifying each sample in a different hydraulic flow unit; decreased ultrasonic velocity; and did not lead to a loss of tensile strength. Tensile strength maintenance occurred due to the low abundance and minor dissolution of siderite, the stability of quartz, and the relative position of diagenetic ankerite within feldspar. This macroporous-dominated lithofacies is the primary reservoir for the Morrow B Sandstone, and is analogous to other porous sandstone reservoirs. It represents an end-member of a chemomechanically low-risk siliceous CO 2 sequestration and CO 2 -EOR reservoir.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

CZCS view of an oceanic acid waste dump

Plumes from the acid waste dump in the New York Bight were visible in all nine cloud-free Coastal Zone Color Scanner images from April and May 1981. The CZCS subsurface radiance channels displayed consistent spectral characteristics, which consisted of a strong increase in the 550 nm channel and a moderate increase in the 520 nm channel relative to the surrounding coastal waters. The 443 nm channel showed no change or a slight decrease in radiance within the plumes. These anomalous radiances preclude the calculation of pigment in the dump plumes using existing algorithms. However, the high radiances of the 550 nm channel can be used to examine the flow patterns of surface water in the vicinity of the dump.

Elrod, Jane A.↗

Evaluating simplifications of subsurface process representations for field-scale permafrost hydrology models

Abstract. Permafrost degradation within a warming climate poses a significant environmental threat through both the permafrost carbon feedback and damage to human communities and infrastructure. Understanding this threat relies on better understanding and numerical representation of thermo-hydrological permafrost processes and the subsequent accurate prediction of permafrost dynamics. All models include simplified assumptions, implying a tradeoff between model complexity and prediction accuracy. The main purpose of this work is to investigate this tradeoff when applying the following commonly made assumptions: (1) assuming equal density of ice and liquid water in frozen soil, (2) neglecting the effect of cryosuction in unsaturated freezing soil, and (3) neglecting advective heat transport during soil freezing and thaw. This study designed a set of 62 numerical experiments using the Advanced Terrestrial Simulator (ATS v1.2) to evaluate the effects of these choices on permafrost hydrological outputs, including both integrated and pointwise quantities. Simulations were conducted under different climate conditions and soil properties from three different sites in both column- and hillslope-scale configurations. Results showed that amongst the three physical assumptions, soil cryosuction is the most crucial yet commonly ignored process. Neglecting cryosuction, on average, can cause 10 %–20 % error in predicting evaporation, 50 %–60 % error in discharge, 10 %–30 % error in thaw depth, and 10 %–30 % error in soil temperature at 1 m beneath the surface. The prediction error for subsurface temperature and water saturation is more obvious at hillslope scales due to the presence of lateral flux. By comparison, using equal ice–liquid density has a minor impact on most hydrological metrics of interest but significantly affects soil water saturation with an averaged 5 %–15 % error. Neglecting advective heat transport presents the least error, 5 % or even much lower, in most metrics of interest for a large-scale Arctic tundra system without apparent influence caused by localized groundwater flow, and it can decrease the simulation time at hillslope scales by 40 %–80 %. By challenging these commonly made assumptions, this work provides permafrost hydrology scientists an important context for understanding the underlying physical processes, including allowing modelers to better choose the appropriate process representation for a given modeling experiment.

54 ENVIRONMENTAL SCIENCES↗

Numerical investigations to identify environmental factors for field-scale reactive transport of pathogens at riverbank filtration sites

While induced bank filtration is a proven method for facilitating sustainable drinking water production, it is at risk from surface water contaminations (e.g., pathogens). Induced bank filtration and pathogen transport in groundwater have been studied extensively. However, long-term studies that consider real-world conditions are missing. These conditions include seasonal changes to environmental conditions and waterworks operations. Therefore, to analyze the effect of seasonal changes on the transport of human pathogenic viruses and their indicators in induced bank filtration, concentrations of adenoviruses and pathogen indicators were monitored over 16 months at an active bank filtration plant at the Rhine River, in Düsseldorf (Germany). Based on this data, a 2D groundwater model was created in PFLOTRAN that simulated flow, heat transport, conservative transport of chloride and the resulting electrical conductivity, reactive transport of oxygen and nitrate, and colloid-based transport of coliforms, somatic coliphages, and adenoviruses. The results show that reduced travel time was the key factor determining periods with a low removal of coliforms and somatic coliphages in the aquifer. Travel time was controlled by river level variations during rainy seasons, and the waterworks extraction rates during dry seasons. Further, for adenovirus transport, travel distance in the subsurface appeared to be the key factor, while travel time had no significant impact. Coliform removal increased when the colmation layer permeability decreased, while coliphage and adenovirus removal was unaffected by the colmation layer permeability. Seasonal changes in temperature and oxygen content did not significantly impact the removal of coliphages and adenoviruses in groundwater. Denitrifying conditions correlated with a lowered coliform removal, but the modelling could not establish a connection between denitrifying conditions and coliform removal. Our study showed that removal of pathogens and pathogen indicators at induced bank filtration plants varies greatly in time and space (e.g., for coliforms from 1 to 4 log-levels at 20m travel distance), and that adenovirus transport differs considerably from transport of coliforms and somatic coliphages.

59 BASIC BIOLOGICAL SCIENCES↗

Advanced Detection of Wellbore Failure for Safe and Secure Utilization of Subsurface Infrastructure

The main goal of this project was to create a state-of-the-art predictive capability that screens and identifies wellbores that are at the highest risk of catastrophic failure. This capability is critical to a host of subsurface applications, including gas storage, hydrocarbon extraction and storage, geothermal energy development, and waste disposal, which depend on seal integrity to meet U.S. energy demands in a safe and secure manner. In addition to the screening tool, this project also developed several other supporting capabilities to help understand fundamental processes involved in wellbore failure. This included novel experimental methods to characterize permeability and porosity evolution during compressive failure of cement, as well as methods and capabilities for understanding two-phase flow in damaged wellbore systems, and novel fracture-resistant cements made from recycled fibers.

36 MATERIALS SCIENCE↗

Effects of surface runoff and infiltration partition methods on hydrological modeling: A comparison of four schemes in two watersheds in the Northeastern US

The ratio according to which the precipitation/snowmelt is divided into surface and subsurface water is of primary importance to watershed scale hydrological modeling. The present study compared four different surface runoff and infiltration partition methods in the Richards-equation-based SWAT model (RSWAT), i.e., daily curve number method (DC-RSWAT), hourly curve number method (HC-RSWAT), Green-Ampt method (GA-RSWAT), and effective infiltration capacity method (EIC-RSWAT). These four versions of RSWAT and a daily-curve-number-based SWAT (DC-SWAT) were used to simulate daily flow rate (from 2001 to 2015) at the outlets of the Tuckahoe Creek Watershed (TCW) and Greensboro Watershed (GW) of the Chesapeake Bay in Maryland. Parameter sensitivity analysis, calibration, and prediction uncertainty analysis were conducted. The results indicate that surface runoff and infiltration partition methods have pronounced impacts on model performance. Comparative analysis showed that DC-RSWAT outperformed other versions of RSWAT for both calibration (from 2001 to2007) and validation (from 2008 to 2015) periods in the two test watersheds. In general, DC-RSWAT achieved more accurate deterministic flow simulation and reliable uncertainty estimation than other versions. As such, daily-curve-number method is recommended for surface runoff and infiltration partition in RSWAT for simulating stream flow.

hydrology, modeling, streamflow↗