Hydrology of two slopes in subarctic Yukon, Canada

NASA Astrophysics Data System (ADS)

Carey, Sean K.; Woo, Ming-Ko

1999-11-01

Two subarctic forested slopes in central Wolf Creek basin, Yukon, were studied in 1996-1997 to determine the seasonal pattern of the hydrologic processes. A south-facing slope has a dense aspen forest on silty soils with seasonal frost only and a north-facing slope has open stands of black spruce and an organic layer on top of clay sediments with permafrost. Snowmelt is advanced by approximately one month on the south-facing slope due to greater radiation receipt. Meltwater infiltrates its seasonally frozen soil with low ice content, recharging the soil moisture reservoir but yielding no lateral surface or subsurface flow. Summer evaporation depletes this recharged moisture and any additional rainfall input, at the expense of surface or subsurface flow. The north-facing slope with an ice rich substrate hinders deep percolation. Snow meltwater is impounded within the organic layer to produce surface runoff in rills and gullies, and subsurface flow along pipes and within the matrix of the organic soil. During the summer, most subsurface flows are confined to the organic layer which has hydraulic conductivities orders of magnitudes larger than the underlying boulder-clay. Evaporation on the north-facing slope declines as both the frost table and the water table descend in the summer. A water balance of the two slopes demonstrates that vertical processes of infiltration and evaporation dominate moisture exchanges on the south-facing slope, whereas the retardation of deep drainage by frost and by clayey soil on the permafrost slope promotes a strong lateral flow component, principally within the organic layer. These results have the important implication that permafrost slopes and organic horizons are the principal controls on streamflow generation in subarctic catchments.

Interpretation of Ground Temperature Anomalies in Hydrothermal Discharge Areas

NASA Astrophysics Data System (ADS)

Price, A. N.; Lindsey, C.; Fairley, J. P., Jr.

2017-12-01

Researchers have long noted the potential for shallow hydrothermal fluids to perturb near-surface temperatures. Several investigators have made qualitative or semi-quantitative use of elevated surface temperatures; for example, in snowfall calorimetry, or for tracing subsurface flow paths. However, little effort has been expended to develop a quantitative framework connecting surface temperature observations with conditions in the subsurface. Here, we examine an area of shallow subsurface flow at Burgdorf Hot Springs, in the Payette National Forest, north of McCall, Idaho USA. We present a simple analytical model that uses easily-measured surface data to infer the temperatures of laterally-migrating shallow hydrothermal fluids. The model is calibrated using shallow ground temperature measurements and overburden thickness estimates from seismic refraction studies. The model predicts conditions in the shallow subsurface, and suggests that the Biot number may place a more important control on the expression of near-surface thermal perturbations than previously thought. In addition, our model may have application in inferring difficult-to-measure parameters, such as shallow subsurface discharge from hydrothermal springs.

Understanding the rapidity of subsurface storm flow response from a fracture-oriented shallow vadose through a new perspective

NASA Astrophysics Data System (ADS)

Zhao, Peng; Zhao, Pei; Liang, Chuan; Li, Tianyang; Zhou, Baojia

2017-01-01

Velocity and celerity in hydrologic systems are controlled by different mechanisms. Efforts were made through joint sample collection and the use of hydrographs and tracers to understand the rapidity of the subsurface flow response to rainstorms on hourly time scales. Three deep subsurface flows during four natural rainstorm events were monitored. The results show that (1) deeper discharge was observed early in responding rainfall events and yielded a high hydrograph amplitude; (2) a ratio index, k, reflecting the dynamic change of the rainfall perturbation intensity in subsurface flow, might reveal inner causal relationships between the flow index and the tracer signal index. Most values of k were larger than 1 at the perturbation stage but approximated 1 at the no-perturbation stage; and (3) for statistical analysis of tracer signals in subsurface flows, the total standard deviation was 17.2, 11.9, 7.4 and 3.5 at perturbation stages and 4.4, 2.5, 1.1, and 0.95 at the non-perturbation stage for observed events. These events were 3-7 times higher in the former rather than the later, reflecting that the variation of tracer signals primarily occurred under rainfall perturbation. Thus, we affirmed that the dynamic features of rainfall have a key effect on rapid processes because, besides the gravity, mechanical waves originating from dynamic rainfall features are another driving factor for conversion between different types of rainfall mechanical energy. A conceptual model for pressure wave propagation was proposed, in which virtual subsurface flow processes in a heterogeneous vadose zone under rainfall are analogous to the water hammer phenomenon in complex conduit systems. Such an analogy can allow pressure in a shallow vadose to increase and decrease and directly influence the velocity and celerity of the flow reflecting a mechanism for rapid subsurface hydrologic response processes in the shallow vadose zone.

Connectivity in modelling subsurface stormflow and tracer transport in a forested hillslope

NASA Astrophysics Data System (ADS)

Laine-Kaulio, Hanne; Koivusalo, Harri

2016-04-01

Preferential flow characterizes subsurface water movement in forested hillslopes. Connectivity of preferential flowpaths in soil, as well as connections between the preferential flowpaths and the soil matrix, control the flow patterns in hillslopes. We simulated subsurface stormflow and tracer transport within a forested hillslope section in Eastern Finland with a modifiable and detailed, physically-based 3-D dual-permeability model. The objective was to investigate the role of flowpath connectivity in capturing the observed migration velocity and concentration status of a tracer plume during the different stages of an irrigation experiment. Hypotheses to explain the observed transport behaviour were formulated based on observed spatial variations in the soil properties. The hypotheses were then tested by running corresponding model versions against spatial high-frequency data from the irrigation experiment. The 3-D distribution of preferential flowpaths, and the connections and disconnections between them, had the most crucial role in capturing the observed transport event with the model. Lateral by-pass flow in the preferential flow network and the transmissivity feedback phenomenon dominated the tracer transport. However, disconnections in lateral preferential flow restrained the transport, making the transport route devious instead of being directly lateral. These disconnections were caused by spots of cemented soil material. Stones had an opposite effect on the flow route formation compared to the cemented spots because preferential flowpaths were present on stone surfaces. Preferential flowpaths related to stones were linked with other types of preferential flowpaths in soil, facilitating the connectivity of preferential flowpaths within the entire hillslope. To take into account the discontinuity in lateral preferential flow and the tortuosity of the dominant flow routes in the model, an explicit representation of the 3-D distribution of saturated hydraulic conductivity in the preferential flow network was required. In addition, with-depth changing porosity values of both pore domains of soil, i.e., the preferential flowpaths and the soil matrix, as well as explicit descriptions of local water losses into the underlying bedrock, were found important for successful simulations. Other investigated model features, including complexity in the values of the rest of the model parameters and in the calculation routines used in describing the flow connections between the preferential flowpaths and the soil matrix, were not found beneficial for the model outcome.

A reference data set of hillslope rainfall-runoff response, Panola Mountain Research Watershed, United States

USGS Publications Warehouse

Tromp-van, Meerveld; James, A.L.; McDonnell, Jeffery J.; Peters, N.E.

2008-01-01

Although many hillslope hydrologic investigations have been conducted in different climate, topographic, and geologic settings, subsurface stormflow remains a poorly characterized runoff process. Few, if any, of the existing data sets from these hillslope investigations are available for use by the scientific community for model development and validation or conceptualization of subsurface stormflow. We present a high-resolution spatial and temporal rainfall-runoff data set generated from the Panola Mountain Research Watershed trenched experimental hillslope. The data set includes surface and subsurface (bedrock surface) topographic information and time series of lateral subsurface flow at the trench, rainfall, and subsurface moisture content (distributed soil moisture content and groundwater levels) from January to June 2002. Copyright 2008 by the American Geophysical Union.

Flow pathways in the Slapton Wood catchment using temperature as a tracer

NASA Astrophysics Data System (ADS)

Birkinshaw, Stephen J.; Webb, Bruce

2010-03-01

SummaryThis study investigates the potential of temperature as a tracer to provide insights into flow pathways. The approach couples fieldwork and modelling experiments for the Eastergrounds Hollow within the Slapton Wood catchment, South Devon, UK. Measurements in the Eastergrounds Hollow were carried out for soil temperature, spring temperature, and the stream temperature and use was made of an existing 1989-1991 data set for the entire Slapton Wood catchment. The predominant flow in this hollow is a result of subsurface stormflow, and previous work has suggested that the water flows vertically down through the soil and then subsurface stormflow occurs at the soil/bedrock interface where the water is deflected laterally. The depth of the subsurface stormflow was previously thought to be around 2.2 m. However, analysis of the new spring, stream and soil temperature data suggests a deeper pathway for the subsurface stormflow. Modelling of water flow and heat transport was carried out using SHETRAN and this was calibrated to reproduce the water flow in the entire Slapton Wood catchment and soil temperatures in the Eastergrounds Hollow. The model was tested for the entire Eastergrounds Hollow with two different soil depths. A depth of 2.2 m, based on previous knowledge, was unable to reproduce the Eastergrounds spring temperature. A depth of 3.7 m produced an excellent comparison between measured and simulated stream and spring temperatures in the Eastergrounds Hollow. This work suggests that the depth of the flow pathways that produce the subsurface stormflow are deeper than previously thought. It also provides a demonstration on the use of temperature as a tracer to understand flow pathways.

Hillslope threshold response to rainfall: (1) a field based forensic approach

Treesearch

Chris B. Graham; Ross A. Woods; Jeffrey J. McDonnell

2010-01-01

Hillslope threshold response to storm rainfall is poorly understood. Basic questions regarding the type, location, and flow dynamics of lateral, subsurface flow remain unanswered, even at our most intensively studied field sites. Here we apply a forensic approach where we combined irrigation and excavation experiments at the well studied Maimai hillslope to determine...

Spatially explicit simulation of peatland hydrology and carbon dioxide exchange: Influence of mesoscale topography

NASA Astrophysics Data System (ADS)

Sonnentag, O.; Chen, J. M.; Roulet, N. T.; Ju, W.; Govind, A.

2008-06-01

Carbon dynamics in peatlands are controlled, in large part, by their wetness as defined by water table depth and volumetric liquid soil moisture content. A common type of peatland is raised bogs that typically have a multiple-layer canopy of vascular plants over a Sphagnum moss ground cover. Their convex form restricts water supply to precipitation and water is shed toward the margins, usually by lateral subsurface flow. The hydraulic gradient for lateral subsurface flow is governed by the peat surface topography at the mesoscale (˜200 m to 5 km). To investigate the influence of mesoscale topography on wetness, evapotranspiration (ET), and gross primary productivity (GPP) in a bog during the snow-free period, we compare the outputs of a further developed version of the daily Boreal Ecosystem Productivity Simulator (BEPS) with observations made at the Mer Bleue peatland, located near Ottawa, Canada. Explicitly considering mesoscale topography, simulated total ET and GPP correlate well with measured ET (r = 0.91) and derived gross ecosystem productivity (GEP; r = 0.92). Both measured ET and derived GEP are simulated similarly well when mesoscale topography is neglected, but daily simulated values are systematically underestimated by about 10% and 12% on average, respectively, due to greater wetness resulting from the lack of lateral subsurface flow. Owing to the differences in moss surface conductances of water vapor and carbon dioxide with increasing moss water content, the differences in the spatial patterns of simulated total ET and GPP are controlled by the mesotopographic position of the moss ground cover.

Spatially Explicit Simulation of Mesotopographic Controls on Peatland Hydrology and Carbon Fluxes

NASA Astrophysics Data System (ADS)

Sonnentag, O.; Chen, J. M.; Roulet, N. T.

2006-12-01

A number of field carbon flux measurements, paleoecological records, and model simulations have acknowledged the importance of northern peatlands in terrestrial carbon cycling and methane emissions. An important parameter in peatlands that influences both net primary productivity, the net gain of carbon through photosynthesis, and decomposition under aerobic and anaerobic conditions, is the position of the water table. Biological and physical processes involved in peatland carbon dynamics and their hydrological controls operate at different spatial scales. The highly variable hydraulic characteristics of the peat profile and the overall shape of the peat body as defined by its surface topography at the mesoscale (104 m2) are of major importance for peatland water table dynamics. Common types of peatlands include bogs with a slightly domed centre. As a result of the convex profile, their water supply is restricted to atmospheric inputs, and water is mainly shed by shallow subsurface flow. From a modelling perspective the influence of mesotopographic controls on peatland hydrology and thus carbon balance requires that process-oriented models that examine the links between peatland hydrology, ecosystem functioning, and climate must incorporate some form of lateral subsurface flow consideration. Most hydrological and ecological modelling studies in complex terrain explicitly account for the topographic controls on lateral subsurface flow through digital elevation models. However, modelling studies in peatlands often employ simple empirical parameterizations of lateral subsurface flow, neglecting the influence of peatlands low relief mesoscale topography. Our objective is to explicitly simulate the mesotopographic controls on peatland hydrology and carbon fluxes using the Boreal Ecosystem Productivity Simulator (BEPS) adapted to northern peatlands. BEPS is a process-oriented ecosystem model in a remote sensing framework that takes into account peatlands multi-layer canopy through vertically stratified mapped leaf area index. Model outputs are validated against multi-year measurements taken at an eddy-covariance flux tower located within Mer Bleue bog, a typical raised bog near Ottawa, Ontario, Canada. Model results for seasonal water table dynamics and evapotranspiration at daily time steps in 2003 are in good agreement with measurements with R2=0.74 and R2=0.79, respectively, and indicate the suitability of our pursued approach.

Interpretation of Ground Temperature Anomalies in Hydrothermal Discharge Areas

NASA Astrophysics Data System (ADS)

Price, Adam N.; Lindsey, Cary R.; Fairley, Jerry P.

2017-12-01

Researchers have long noted the potential for shallow hydrothermal fluids to perturb near-surface temperatures. Several investigators have made qualitative or semiquantitative use of elevated surface temperatures; for example, in snowfall calorimetry, or for tracing subsurface flow paths. However, a quantitative framework connecting surface temperature observations with conditions in the subsurface is currently lacking. Here, we model an area of shallow subsurface flow at Burgdorf Hot Springs, a rustic commercial resort in the Payette National Forest, north of McCall, ID, USA. We calibrate the model using shallow (0.2 m depth) ground temperature measurements and overburden thickness estimates from seismic refraction studies. The calibrated model predicts negligible loss of heat energy from the laterally migrating fluids at the Burgdorf site, in spite of the fact that thermal anomalies are observed in the unconsolidated near-surface alluvium. Although elevated near-surface ground temperatures are commonly assumed to result from locally high heat flux, this conflicts with the small apparent heat loss during lateral flow inferred at the Burgdorf site. We hypothesize an alternative explanation for near-surface temperature anomalies that is only weakly dependent on heat flux, and more strongly controlled by the Biot number, a dimensionless parameter that compares the rate at which convection carries heat away from the land surface to the rate at which it is supplied by conduction to the interface.

The use of distributed hydrological models for the Gard 2002 flash flood event: Analysis of associated hydrological processes

NASA Astrophysics Data System (ADS)

Braud, Isabelle; Roux, Hélène; Anquetin, Sandrine; Maubourguet, Marie-Madeleine; Manus, Claire; Viallet, Pierre; Dartus, Denis

2010-11-01

SummaryThis paper presents a detailed analysis of the September 8-9, 2002 flash flood event in the Gard region (southern France) using two distributed hydrological models: CVN built within the LIQUID® hydrological platform and MARINE. The models differ in terms of spatial discretization, infiltration and water redistribution representation, and river flow transfer. MARINE can also account for subsurface lateral flow. Both models are set up using the same available information, namely a DEM and a pedology map. They are forced with high resolution radar rainfall data over a set of 18 sub-catchments ranging from 2.5 to 99 km2 and are run without calibration. To begin with, models simulations are assessed against post field estimates of the time of peak and the maximum peak discharge showing a fair agreement for both models. The results are then discussed in terms of flow dynamics, runoff coefficients and soil saturation dynamics. The contribution of the subsurface lateral flow is also quantified using the MARINE model. This analysis highlights that rainfall remains the first controlling factor of flash flood dynamics. High rainfall peak intensities are very influential of the maximum peak discharge for both models, but especially for the CVN model which has a simplified overland flow transfer. The river bed roughness also influences the peak intensity and time. Soil spatial representation is shown to have a significant role on runoff coefficients and on the spatial variability of saturation dynamics. Simulated soil saturation is found to be strongly related with soil depth and initial storage deficit maps, due to a full saturation of most of the area at the end of the event. When activated, the signature of subsurface lateral flow is also visible in the spatial patterns of soil saturation with higher values concentrating along the river network. However, the data currently available do not allow the assessment of both patterns. The paper concludes with a set of recommendations for enhancing field observations in order to progress in process understanding and gather a larger set of data to improve the realism of distributed models.

Field study and simulation of diurnal temperature effects on infiltration and variably saturated flow beneath an ephemeral stream

USGS Publications Warehouse

Dudek Ronan, Anne; Prudic, David E.; Thodal, Carl E.; Constantz, Jim

1998-01-01

Two experiments were performed to investigate flow beneath an ephemeral stream and to estimate streambed infiltration rates. Discharge and stream-area measurements were used to determine infiltration rates. Stream and subsurface temperatures were used to interpret subsurface flow through variably saturated sediments beneath the stream. Spatial variations in subsurface temperatures suggest that flow beneath the streambed is dependent on the orientation of the stream in the canyon and the layering of the sediments. Streamflow and infiltration rates vary diurnally: Streamflow is lowest in late afternoon when stream temperature is greatest and highest in early morning when stream temperature is least. The lower afternoon Streamflow is attributed to increased infiltration rates; evapotranspiration is insufficient to account for the decreased Streamflow. The increased infiltration rates are attributed to viscosity effects on hydraulic conductivity from increased stream temperatures. The first set of field data was used to calibrate a two-dimensional variably saturated flow model that includes heat transport. The model was calibrated to (1) temperature fluctuations in the subsurface and (2) infiltration rates determined from measured Streamflow losses. The second set of field data was to evaluate the ability to predict infiltration rates on the basis of temperature measurements alone. Results indicate that the variably saturated subsurface flow depends on downcanyon layering of the sediments. They also support the field observations in indicating that diurnal changes in infiltration can be explained by temperature dependence of hydraulic conductivity. Over the range of temperatures and flows monitored, diurnal stream temperature changes can be used to estimate streambed infiltration rates. It is often impractical to maintain equipment for determining infiltration rates by traditional means; however, once a model is calibrated using both infiltration and temperature data, only relatively inexpensive temperature monitoring can later yield infiltration rates that are within the correct order of magnitude.

The soil hydrologic response to forest regrowth: a case study from southwestern Amazonia

NASA Astrophysics Data System (ADS)

Godsey, Sarah; Elsenbeer, Helmut

2002-05-01

As a large and dynamic land-use category, tropical secondary forests may affect climate, soils, and hydrology in a manner different from primary forests or agricultural areas. We investigated the saturated hydraulic conductivity Ksat of a Kandiudult under different land uses in Rondonia, Brazil. We measured Ksat at four depths (12·5, 20, 30 and 50 cm) under (a) primary forest, (b) a former banana-cacao plantation (SF1), and (c) an abandoned pasture (SF2). At 12·5 cm, all three land uses differ significantly ( = 0·1), but not at the 20 and 30 cm depths. At 50 cm, Ksat was significantly greater in the former pasture than in other land uses. Lateral subsurface flow is expected during intense rainfall (about 30 times per year) at 30 cm depth in SF1 and at 50 cm depth in the forest, whereas the relatively low permeability at shallow 12·5 cm in the SF2 may result not only in lateral subsurface flow, but also saturation overland flow. For modelling purposes, recovering systems seem to have Ksat values distinct from primary forest at shallow depths, whereas at deeper layers (>20 cm) they may be considered similar to forests.

Demonstration of the Tilting of the Gas-Water Interface under Hydrodynamic Conditions.

ERIC Educational Resources Information Center

Gretener, P. E.

1979-01-01

Describes the construction of an apparatus to demonstrate the tilting of an oil-water, gas-water, or gas-oil interface when the subsurface reservoir is under hydrodynamic conditions (i.e., when conditions of lateral flow exist). The model can be constructed of readily-available materials. (RE)

Evaluation of nutrient removal efficiency and microbial enzyme activity in a baffled subsurface-flow constructed wetland system

Treesearch

Lihua Cui; Ying Ouyang; Wenjie Gu; Weozhi Yang; Qiaoling Xu

2013-01-01

In this study, the enzyme activities and their relationships to domestic wastewater purification are investigated in four different types of subsurface-flow constructed wetlands (CWs), namely the traditional horizontal subsurface-flow, horizontal baffled subsurface-flow, vertical baffled subsurface-flow, and composite baffled subsurface-flow CWs. Results showed that...

Evaluating post-wildfire hydrologic recovery using ParFlow in southern California

NASA Astrophysics Data System (ADS)

Lopez, S. R.; Kinoshita, A. M.; Atchley, A. L.

2016-12-01

Wildfires are naturally occurring hazards that can have catastrophic impacts. They can alter the natural processes within a watershed, such as surface runoff and subsurface water storage. Generally, post-fire hydrologic models are either one-dimensional, empirically-based models, or two-dimensional, conceptually-based models with lumped parameter distributions. These models are useful in providing runoff measurements at the watershed outlet; however, do not provide distributed hydrologic simulation at each point within the watershed. This research demonstrates how ParFlow, a three-dimensional, distributed hydrologic model can simulate post-fire hydrologic processes by representing soil burn severity (via hydrophobicity) and vegetation recovery as they vary both spatially and temporally. Using this approach, we are able to evaluate the change in post-fire water components (surface flow, lateral flow, baseflow, and evapotranspiration). This model is initially developed for a hillslope in Devil Canyon, burned in 2003 by the Old Fire in southern California (USA). The domain uses a 2m-cell size resolution over a 25 m by 25 m lateral extent. The subsurface reaches 2 m and is assigned a variable cell thickness, allowing an explicit consideration of the soil burn severity throughout the stages of recovery and vegetation regrowth. Vegetation regrowth is incorporated represented by satellite-based Enhanced Vegetation Index (EVI) products. The pre- and post-fire surface runoff, subsurface storage, and surface storage interactions are evaluated and will be used as a basis for developing a watershed-scale model. Long-term continuous simulations will advance our understanding of post-fire hydrological partitioning between water balance components and the spatial variability of watershed processes, providing improved guidance for post-fire watershed management.

Subsurface lateral preferential flow network revealed by time-lapse ground-penetrating radar in a hillslope

NASA Astrophysics Data System (ADS)

Guo, Li; Chen, Jin; Lin, Henry

2014-12-01

Subsurface lateral preferential flow (LPF) has been observed to contribute substantially to hillslope and catchment runoff. However, the complex nature of LPF and the lack of an appropriate investigation method have hindered direct LPF observation in the field. Thus, the initiation, persistence, and dynamics of LPF networks remain poorly understood. This study explored the application of time-lapse ground-penetrating radar (GPR) together with an artificial infiltration to shed light on the nature of LPF and its dynamics in a hillslope. Based on our enhanced field experimental setup and carefully refined GPR data postprocessing algorithms, we developed a new protocol to reconstruct LPF networks with centimeter resolution. This is the first time that a detailed LPF network and its dynamics have been revealed noninvasively along a hillslope. Real-time soil water monitoring and field soil investigation confirmed the locations of LPF mapped by time-lapse GPR surveys. Our results indicated the following: (1) Increased spatial variations of radar signals after infiltration suggested heterogeneous soil water changes within the studied soil, which reflected the generation and dynamics of LPF; (2) Two types of LPF networks were identified, the network at the location of soil permeability contrasts and that formed via a series of connected preferential flow paths; and (3) The formation and distribution of LPF networks were influenced by antecedent soil water condition. Overall, this study demonstrates clearly that carefully designed time-lapse GPR surveys with enhanced data postprocessing offer a practical and nondestructive way of mapping LPF networks in the field, thereby providing a potentially significant enhancement in our ability to study complex subsurface flow processes across the landscape.

Lateral, Vertical, and Longitudinal Source Area Connectivity Drive Runoff and Carbon Export Across Watershed Scales

NASA Astrophysics Data System (ADS)

Zimmer, Margaret A.; McGlynn, Brian L.

2018-03-01

Watersheds are three-dimensional hydrologic systems where the longitudinal expansion/contraction of stream networks, vertical connection/disconnection between shallow and deep groundwater systems, and lateral connectivity of these water sources to streams mediate runoff production and nutrient export. The connectivity of runoff source areas during both baseflow and stormflow conditions and their combined influence on biogeochemical fluxes remain poorly understood. Here we focused on a set of 3.3 and 48.4 ha nested watersheds (North Carolina, USA). These watersheds comprise ephemeral and intermittent runoff-producing headwaters and perennial runoff-producing lowlands. Within these landscape elements, we characterized the timing and magnitude of precipitation, runoff, and runoff-generating flow paths. The active surface drainage network (ASDN) reflected connectivity to, and contributions from, source areas that differed under baseflow and stormflow conditions. The baseflow-associated ASDN expanded and contracted seasonally, driven by the rise and fall of the seasonal water table. Superimposed on this were event-activated source area contributions driven by connectivity to surficial and shallow subsurface flow paths. Frequently activated shallow flow paths also caused increased in-stream dissolved organic carbon (DOC) concentrations with increases in runoff across both watershed scales. The spread and variability within this DOC-runoff relationship was driven by a seasonal depletion of DOC from continual shallow subsurface flow path activation and subsequent replenishment from autumn litterfall. Our findings suggest that hydrobiogeochemical signals at larger watershed outlets can be driven by the expansion, contraction, and connection of lateral, longitudinal, and vertical source areas with distinct runoff generation processes.

Dynamics of nonreactive solute transport in the permafrost environment

NASA Astrophysics Data System (ADS)

Svyatskiy, D.; Coon, E. T.; Moulton, J. D.

2017-12-01

As part of the DOE Office of Science Next Generation Ecosystem Experiment, NGEE-Arctic, researchers are developing process-rich models to understand and predict the evolution of water sources and hydrologic flow pathways resulting from degrading permafrost. The sources and interaction of surface and subsurface water and flow paths are complex in space and time due to strong interplay between heterogeneous subsurface parameters, the seasonal to decadal evolution of the flow domain, climate driven melting and release of permafrost ice as a liquid water source, evolving surface topography and highly variable meteorological data. In this study, we seek to characterize the magnitude of vertical and lateral subsurface flows in a cold, wet tundra, polygonal landscape characteristic of the Barrow Peninsula, AK. To better understand the factors controlling water flux partitioning in these low gradient landscapes, NGEE researchers developed and are applying the Advanced Terrestrial Simulator (ATS), which fully couples surface and subsurface flow and energy processes, snow distribution and atmospheric forcing. Here we demonstrate the integration of a new solute transport model within the ATS, which enables the interpretation of applied and natural tracer experiments and observations aimed at quantifying water sources and flux partitioning. We examine the role of ice wedge polygon structure, freeze-thaw processes and soil properties on the seasonal transport of water within and through polygons features, and compare results to tracer experiments on 2D low-centered and high-centered transects corresponding to artificial as well as realistic topographical data from sites in polygonal tundra. These simulations demonstrate significant difference between flow patterns between permafrost and non-permafrost environments due to active layer freeze-thaw processes.

Geoelectrical monitoring of simulated subsurface leakage to support high-hazard nuclear decommissioning at the Sellafield Site, UK.

PubMed

Kuras, Oliver; Wilkinson, Paul B; Meldrum, Philip I; Oxby, Lucy S; Uhlemann, Sebastian; Chambers, Jonathan E; Binley, Andrew; Graham, James; Smith, Nicholas T; Atherton, Nick

2016-10-01

A full-scale field experiment applying 4D (3D time-lapse) cross-borehole Electrical Resistivity Tomography (ERT) to the monitoring of simulated subsurface leakage was undertaken at a legacy nuclear waste silo at the Sellafield Site, UK. The experiment constituted the first application of geoelectrical monitoring in support of decommissioning work at a UK nuclear licensed site. Images of resistivity changes occurring since a baseline date prior to the simulated leaks revealed likely preferential pathways of silo liquor simulant flow in the vadose zone and upper groundwater system. Geophysical evidence was found to be compatible with historic contamination detected in permeable facies in sediment cores retrieved from the ERT boreholes. Results indicate that laterally discontinuous till units forming localized hydraulic barriers substantially affect flow patterns and contaminant transport in the shallow subsurface at Sellafield. We conclude that only geophysical imaging of the kind presented here has the potential to provide the detailed spatial and temporal information at the (sub-)meter scale needed to reduce the uncertainty in models of subsurface processes at nuclear sites. Copyright © 2016 British Geological Survey, NERC. Published by Elsevier B.V. All rights reserved.

Factors Influencing the Sahelian Paradox at the Local Watershed Scale: Causal Inference Insights

NASA Astrophysics Data System (ADS)

Van Gordon, M.; Groenke, A.; Larsen, L.

2017-12-01

While the existence of paradoxical rainfall-runoff and rainfall-groundwater correlations are well established in the West African Sahel, the hydrologic mechanisms involved are poorly understood. In pursuit of mechanistic explanations, we perform a causal inference analysis on hydrologic variables in three watersheds in Benin and Niger. Using an ensemble of techniques, we compute the strength of relationships between observational soil moisture, runoff, precipitation, and temperature data at seasonal and event timescales. Performing analysis over a range of time lags allows dominant time scales to emerge from the relationships between variables. By determining the time scales of hydrologic connectivity over vertical and lateral space, we show differences in the importance of overland and subsurface flow over the course of the rainy season and between watersheds. While previous work on the paradoxical hydrologic behavior in the Sahel focuses on surface processes and infiltration, our results point toward the importance of subsurface flow to rainfall-runoff relationships in these watersheds. The hypotheses generated from our ensemble approach suggest that subsequent explorations of mechanistic hydrologic processes in the region include subsurface flow. Further, this work highlights how an ensemble approach to causal analysis can reveal nuanced relationships between variables even in poorly understood hydrologic systems.

Effects of rainfall patterns and land cover on the subsurface flow generation of sloping Ferralsols in southern China

PubMed Central

Yang, Jie; Tang, Chongjun; Chen, Lihua; Liu, Yaojun; Wang, Lingyun

2017-01-01

Rainfall patterns and land cover are two important factors that affect the runoff generation process. To determine the surface and subsurface flows associated with different rainfall patterns on sloping Ferralsols under different land cover types, observational data related to surface and subsurface flows from 5 m × 15 m plots were collected from 2010 to 2012. The experiment was conducted to assess three land cover types (grass, litter cover and bare land) in the Jiangxi Provincial Soil and Water Conservation Ecological Park. During the study period, 114 natural rainfall events produced subsurface flow and were divided into four groups using k-means clustering according to rainfall duration, rainfall depth and maximum 30-min rainfall intensity. The results showed that the total runoff and surface flow values were highest for bare land under all four rainfall patterns and lowest for the covered plots. However, covered plots generated higher subsurface flow values than bare land. Moreover, the surface and subsurface flows associated with the three land cover types differed significantly under different rainfall patterns. Rainfall patterns with low intensities and long durations created more subsurface flow in the grass and litter cover types, whereas rainfall patterns with high intensities and short durations resulted in greater surface flow over bare land. Rainfall pattern I had the highest surface and subsurface flow values for the grass cover and litter cover types. The highest surface flow value and lowest subsurface flow value for bare land occurred under rainfall pattern IV. Rainfall pattern II generated the highest subsurface flow value for bare land. Therefore, grass or litter cover are able to convert more surface flow into subsurface flow under different rainfall patterns. The rainfall patterns studied had greater effects on subsurface flow than on total runoff and surface flow for covered surfaces, as well as a greater effect on surface flows associated with bare land. PMID:28792507

On the characterization of subsurface flow and hydraulic conductivity from surface SP measurements: correcting for electrical heterogeneities.

NASA Astrophysics Data System (ADS)

Sailhac, P.; Marquis, G.; Darnet, M.; Szalai, S.

2003-04-01

Surface self potential measurements (SP) are useful to characterize underground fluid flow or chemical reactions (as redox) and can be used in addition to NMR and electrical prospecting in hydrological investigations. Assuming that the SP anomalies have an electrokinetic origin, the source of SP data is the divergence of underground fluid flow; one important problem with surface SP data is then its interpretation in terms of fluid flow geometry. Some integral transform techniques have been shown to be powerful for SP interpretation (e.g. Fournier 1989, Patella, 1997; Sailhac &Marquis 2001). All these techniques are based upon Green’{ }s functions to characterize underground water flow, but they assume a constant electrical conductivity in the subsurface. This unrealistic approximation results in the appearance of non-electrokinetic sources at strong lateral electrical conductivity contrasts. We present here new Green’{ }s functions suitable for media of heterogeneous electrical conductivity. This new approach allows the joint interpretation of electrical resistivity tomography and SP measurements to detect electrokinetic sources caused by fluid flow. Tests on synthetic examples show that it gives more realistic results that when a constant electrical conductivity is assumed.

Improving National Water Modeling: An Intercomparison of two High-Resolution, Continental Scale Models, CONUS-ParFlow and the National Water Model

NASA Astrophysics Data System (ADS)

Tijerina, D.; Gochis, D.; Condon, L. E.; Maxwell, R. M.

2017-12-01

Development of integrated hydrology modeling systems that couple atmospheric, land surface, and subsurface flow is growing trend in hydrologic modeling. Using an integrated modeling framework, subsurface hydrologic processes, such as lateral flow and soil moisture redistribution, are represented in a single cohesive framework with surface processes like overland flow and evapotranspiration. There is a need for these more intricate models in comprehensive hydrologic forecasting and water management over large spatial areas, specifically the Continental US (CONUS). Currently, two high-resolution, coupled hydrologic modeling applications have been developed for this domain: CONUS-ParFlow built using the integrated hydrologic model ParFlow and the National Water Model that uses the NCAR Weather Research and Forecasting hydrological extension package (WRF-Hydro). Both ParFlow and WRF-Hydro include land surface models, overland flow, and take advantage of parallelization and high-performance computing (HPC) capabilities; however, they have different approaches to overland subsurface flow and groundwater-surface water interactions. Accurately representing large domains remains a challenge considering the difficult task of representing complex hydrologic processes, computational expense, and extensive data needs; both models have accomplished this, but have differences in approach and continue to be difficult to validate. A further exploration of effective methodology to accurately represent large-scale hydrology with integrated models is needed to advance this growing field. Here we compare the outputs of CONUS-ParFlow and the National Water Model to each other and with observations to study the performance of hyper-resolution models over large domains. Models were compared over a range of scales for major watersheds within the CONUS with a specific focus on the Mississippi, Ohio, and Colorado River basins. We use a novel set of approaches and analysis for this comparison to better understand differences in process and bias. This intercomparison is a step toward better understanding how much water we have and interactions between surface and subsurface. Our goal is to advance our understanding and simulation of the hydrologic system and ultimately improve hydrologic forecasts.

Modeling streamflow in a snow-dominated forest watershed using the Water Erosion Prediction Project (WEPP) model

Treesearch

A. Srivastava; J. Q. Wu; W. J. Elliot; E. S. Brooks; D. C. Flanagan

2017-01-01

The Water Erosion Prediction Project (WEPP) model was originally developed for hillslope and small watershed applications. Recent improvements to WEPP have led to enhanced computations for deep percolation, subsurface lateral flow, and frozen soil. In addition, the incorporation of channel routing has made the WEPP model well suited for large watersheds with perennial...

Paleomagnetic correlation of basalt flows in selected coreholes near the Advanced Test Reactor Complex, the Idaho Nuclear Technology and Engineering Center, and along the southern boundary, Idaho National Laboratory, Idaho

USGS Publications Warehouse

Hodges, Mary K.V.; Champion, Duane E.

2016-10-03

The U.S. Geological Survey, in cooperation with the U.S. Department of Energy, used paleomagnetic data from 18 coreholes to construct three cross sections of subsurface basalt flows in the southern part of the Idaho National Laboratory (INL). These cross sections, containing descriptions of the subsurface horizontal and vertical distribution of basalt flows and sediment layers, will be used in geological studies, and to construct numerical models of groundwater flow and contaminant transport.Subsurface cross sections were used to correlate surface vents to their subsurface flows intersected by coreholes, to correlate subsurface flows between coreholes, and to identify possible subsurface vent locations of subsurface flows. Correlations were identified by average paleomagnetic inclinations of flows, and depth from land surface in coreholes, normalized to the North American Datum of 1927. Paleomagnetic data were combined, in some cases, with other data, such as radiometric ages of flows. Possible vent locations of buried basalt flows were identified by determining the location of the maximum thickness of flows penetrated by more than one corehole.Flows from the surface volcanic vents Quaking Aspen Butte, Vent 5206, Mid Butte, Lavatoo Butte, Crater Butte, Pond Butte, Vent 5350, Vent 5252, Tin Cup Butte, Vent 4959, Vent 5119, and AEC Butte are found in coreholes, and were correlated to the surface vents by matching their paleomagnetic inclinations, and in some cases, their stratigraphic positions.Some subsurface basalt flows that do not correlate to surface vents, do correlate over several coreholes, and may correlate to buried vents. Subsurface flows which correlate across several coreholes, but not to a surface vent include the D3 flow, the Big Lost flow, the CFA buried vent flow, the Early, Middle, and Late Basal Brunhes flows, the South Late Matuyama flow, the Matuyama flow, and the Jaramillo flow. The location of vents buried in the subsurface by younger basalt flows can be inferred if their flows are penetrated by several coreholes, by tracing the flows in the subsurface, and determining where the greatest thickness occurs.

Performance of the subsurface flow constructed wetlands for pretreatment of slightly polluted source water.

PubMed

Yang, Xu; Zhang, Xueping; Wang, Jifu; Zhao, Guangying; Wang, Baojian

2014-05-01

The slightly polluted source water of Yellow River was pretreated in a horizontal subsurface flow constructed wetland (HSFCW) and a lateral subsurface flow constructed wetland (LSFCW) in the Ji'nan city Reservoir, Shandong, China. During almost one years run, the results showed that at the hydraulic loading rate of 1 m/day, the removal efficiencies of chemical oxygen demand (COD), total nitrogen (TN), ammonium nitrogen (NH4 (+)-N) and total phosphorus (TP) in the HSFCW were 48.9, 51.4, 48.7 and 48.9 %, respectively, and the corresponding removal efficiencies in the LSFCW were 50.51, 53.12, 50.44 and 50.83 %, respectively. The HSFCW and LSFCW had a similar high potential for nutrients removal and LSFCW was slightly better. According to the China standard for surface water resources (GB3838-2002), mean effluent COD can reach the Class I (≤ 15 mg/L), and NH4 (+)-N and TP and TN can reach nearly the Class I (≤ 0.015 mg/L), the Class III (≤ 0.05 mg/L) and the Class IV (≤ 1.5 mg/L), respectively. It can be concluded that the slightly polluted source water from Reservoir was pretreated well by the constructed wetland.

Modeling post-wildfire hydrological processes with ParFlow

NASA Astrophysics Data System (ADS)

Escobar, I. S.; Lopez, S. R.; Kinoshita, A. M.

2017-12-01

Wildfires alter the natural processes within a watershed, such as surface runoff, evapotranspiration rates, and subsurface water storage. Post-fire hydrologic models are typically one-dimensional, empirically-based models or two-dimensional, conceptually-based models with lumped parameter distributions. These models are useful for modeling and predictions at the watershed outlet; however, do not provide detailed, distributed hydrologic processes at the point scale within the watershed. This research uses ParFlow, a three-dimensional, distributed hydrologic model to simulate post-fire hydrologic processes by representing the spatial and temporal variability of soil burn severity (via hydrophobicity) and vegetation recovery. Using this approach, we are able to evaluate the change in post-fire water components (surface flow, lateral flow, baseflow, and evapotranspiration). This work builds upon previous field and remote sensing analysis conducted for the 2003 Old Fire Burn in Devil Canyon, located in southern California (USA). This model is initially developed for a hillslope defined by a 500 m by 1000 m lateral extent. The subsurface reaches 12.4 m and is assigned a variable cell thickness to explicitly consider soil burn severity throughout the stages of recovery and vegetation regrowth. We consider four slope and eight hydrophobic layer configurations. Evapotranspiration is used as a proxy for vegetation regrowth and is represented by the satellite-based Simplified Surface Energy Balance (SSEBOP) product. The pre- and post-fire surface runoff, subsurface storage, and surface storage interactions are evaluated at the point scale. Results will be used as a basis for developing and fine-tuning a watershed-scale model. Long-term simulations will advance our understanding of post-fire hydrological partitioning between water balance components and the spatial variability of watershed processes, providing improved guidance for post-fire watershed management. In reference to the presenter, Isabel Escobar: Research is funded by the NASA-DIRECT STEM Program. Travel expenses for this presentation is funded by CSU-LSAMP. CSU-LSAMP is supported by the National Science Foundation under Grant # HRD-1302873 and the CSU Office of Chancellor.

Hillslope-scale experiment demonstrates role of convergence during two-step saturation

USGS Publications Warehouse

Gevaert, A. I.; Teuling, A. J.; Uijlenhoet, R.; DeLong, Stephen B.; Huxman, T. E.; Pangle, L. A.; Breshears, David D.; Chorover, J.; Pelletier, John D.; Saleska, S. R.; Zeng, X.; Troch, Peter A.

2014-01-01

Subsurface flow and storage dynamics at hillslope scale are difficult to ascertain, often in part due to a lack of sufficient high-resolution measurements and an incomplete understanding of boundary conditions, soil properties, and other environmental aspects. A continuous and extreme rainfall experiment on an artificial hillslope at Biosphere 2's Landscape Evolution Observatory (LEO) resulted in saturation excess overland flow and gully erosion in the convergent hillslope area. An array of 496 soil moisture sensors revealed a two-step saturation process. First, the downward movement of the wetting front brought soils to a relatively constant but still unsaturated moisture content. Second, soils were brought to saturated conditions from below in response to rising water tables. Convergent areas responded faster than upslope areas, due to contributions from lateral subsurface flow driven by the topography of the bottom boundary, which is comparable to impermeable bedrock in natural environments. This led to the formation of a groundwater ridge in the convergent area, triggering saturation excess runoff generation. This unique experiment demonstrates, at very high spatial and temporal resolution, the role of convergence on subsurface storage and flow dynamics. The results bring into question the representation of saturation excess overland flow in conceptual rainfall-runoff models and land-surface models, since flow is gravity-driven in many of these models and upper layers cannot become saturated from below. The results also provide a baseline to study the role of the co-evolution of ecological and hydrological processes in determining landscape water dynamics during future experiments in LEO.

Modeling solute transport in a heterogeneous unsaturated porous medium under dynamic boundary conditions on different spatial scales

NASA Astrophysics Data System (ADS)

Cremer, Clemens; Neuweiler, Insa; Bechtold, Michel

2013-04-01

Understanding transport of solutes/contaminants through unsaturated soil in the shallow subsurface is vital to assess groundwater quality, nutrient cycling or to plan remediation projects. Alternating precipitation and evaporation conditions causing upward and downward flux with differing flow paths, changes in saturation and related structural heterogeneity make the description of transport in the unsaturated zone near the soil-surface a complex problem. Preferential flow paths strongly depend, among other things, on the saturation of a medium. Recent studies (e.g. Bechtold et al., 2011) showed lateral flow and solute transport during evaporation conditions (upward flux) in vertically layered sand columns. Results revealed that during evaporation water and solute are redistributed laterally from coarse to fine media deeper in the soil, and towards zones of lowest hydraulic head near to the soil surface. These zones at the surface can be coarse or fine grained depending on saturation status and evaporation flux. However, if boundary conditions are reversed and precipitation is applied, the flow field is not reversed in the same manner, resulting in entirely different transport patterns for downward and upward flow. Therefore, considering net-flow rates alone is misleading when describing transport in the shallow unsaturated zone. In this contribution, we analyze transport of a solute in the shallow subsurface to assess effects resulting from the superposition of heterogeneous soil structures and dynamic flow conditions on various spatial scales. Two-dimensional numerical simulations of unsaturated flow and transport in heterogeneous porous media under changing boundary conditions are carried out using a finite-volume code coupled to a particle tracking algorithm to quantify solute transport and leaching rates. In order to validate numerical simulations, results are qualitatively compared to those of a physical experiment (Bechtold et al., 2011). Numerical simulations differ in lateral scale reaching from 0.2 m to 1.5 m, while the height of the domain is kept constant to 1.5m. Strong material heterogeneity is realized through vertical layers of coarse and fine sand. Both materials remain permanently under liquid-flow-dominated ('stage1') evaporation conditions. Spatial moments as well as the dilution index (Kitanidis, 1994) are used for quantification of transport behaviour. Results show that, while all simulations led to anomalous transport, infiltration-evaporation cycles lead to faster solute leaching rates than solely infiltration at the same net-infiltration rate in both homogeneous and heterogeneous media. Flow and transport-paths significantly differed between infiltration and evaporation, resulting in lateral water fluxes and hence lateral solute transport. Variation of the width of the model domain shows faster leaching rates for domains with small horizontal extent.

Redistribution of soil metals and organic carbon via lateral flowpaths at the catchment scale in a glaciated upland setting

Treesearch

Rebecca R. Bourgault; Donald S. Ross; Scott W. Bailey; Thomas D. Bullen; Kevin J. McGuire; John P. Gannon

2017-01-01

Emerging evidence shows that interactions between soils and subsurface flow paths contribute to spatial variations in stream water chemistry in headwater catchments. However, few have yet attempted to quantify chemical variations in soils at catchment and hillslope scales. Watershed 3 (WS3) at Hubbard Brook Experimental Forest, New Hampshire, USA, was studied in order...

Contrasting styles of large-scale displacement of unconsolidated sand: examples from the early Jurassic Navajo Sandstone on the Colorado Plateau, USA

NASA Astrophysics Data System (ADS)

Bryant, Gerald

2015-04-01

Large-scale soft-sediment deformation features in the Navajo Sandstone have been a topic of interest for nearly 40 years, ever since they were first explored as a criterion for discriminating between marine and continental processes in the depositional environment. For much of this time, evidence for large-scale sediment displacements was commonly attributed to processes of mass wasting. That is, gravity-driven movements of surficial sand. These slope failures were attributed to the inherent susceptibility of dune sand responding to environmental triggers such as earthquakes, floods, impacts, and the differential loading associated with dune topography. During the last decade, a new wave of research is focusing on the event significance of deformation features in more detail, revealing a broad diversity of large-scale deformation morphologies. This research has led to a better appreciation of subsurface dynamics in the early Jurassic deformation events recorded in the Navajo Sandstone, including the important role of intrastratal sediment flow. This report documents two illustrative examples of large-scale sediment displacements represented in extensive outcrops of the Navajo Sandstone along the Utah/Arizona border. Architectural relationships in these outcrops provide definitive constraints that enable the recognition of a large-scale sediment outflow, at one location, and an equally large-scale subsurface flow at the other. At both sites, evidence for associated processes of liquefaction appear at depths of at least 40 m below the original depositional surface, which is nearly an order of magnitude greater than has commonly been reported from modern settings. The surficial, mass flow feature displays attributes that are consistent with much smaller-scale sediment eruptions (sand volcanoes) that are often documented from modern earthquake zones, including the development of hydraulic pressure from localized, subsurface liquefaction and the subsequent escape of fluidized sand toward the unconfined conditions of the surface. The origin of the forces that produced the lateral, subsurface movement of a large body of sand at the other site is not readily apparent. The various constraints on modeling the generation of the lateral force required to produce the observed displacement are considered here, along with photodocumentation of key outcrop relationships.

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.

Retention and transport of nutrients in a third-order stream in northwestern California; hyporheic processes

USGS Publications Warehouse

Triska, F.J.; Kennedy, V.C.; Avanzino, R.J.; Zellweger, G.W.; Bencala, K.E.

1989-01-01

Chloride and nitrate were coinjected into the surface waters of a third-order stream for 20 d to exmaine solute retention, and the fate of nitrate during subsurface transport. A series of wells (shallow pits) 0.5-10 m from the adjacent channel were sampled to estimate the lateral interflow of water. Two subsurface return flows beneath the wetted channel were also examined. Results indicated that the capacity of the hyporheic zone for transient solute storage and as potential biological habitat varies with channel morphology, bed roughness, and permeability. A conceptual model that considers the groundwater-stream water interface as the fluvial boundary is proposed. -from Authors

Enhanced recharge rates and altered recharge sensitivity to climate variability through subsurface heterogeneity

NASA Astrophysics Data System (ADS)

Hartmann, Andreas; Gleeson, Tom; Wada, Yoshihide; Wagener, Thorsten

2017-04-01

Karst aquifers in Europe are an important source of fresh water contributing up to half of the total drinking water supply in some countries. Karstic groundwater recharge is one of the most important components of the water balance of karst systems as it feeds the karst aquifers. Presently available large-scale hydrological models do not consider karst heterogeneity adequately. Projections of current and potential future groundwater recharge of Europe's karst aquifers are therefore unclear. In this study we compare simulations of present (1991-2010) and future (2080-2099) recharge using two different models to simulate groundwater recharge processes. One model includes karst processes (subsurface heterogeneity, lateral flow and concentrated recharge), while the other is based on the conceptual understanding of common hydrological systems (homogeneous subsurface, saturation excess overland flow). Both models are driven by the bias-corrected 5 GCMs of the ISI-MIP project (RCP8.5). To further assess sensitivity of groundwater recharge to climate variability, we calculate the elasticity of recharge rates to annual precipitation, temperature and average intensity of rainfall events, which is the median change of recharge that corresponds to the median change of these climate variables within the present and future time period, respectively. Our model comparison shows that karst regions over Europe have enhanced recharge rates with greater inter-annual variability compared to those with more homogenous subsurface properties. Furthermore, the heterogeneous representation shows stronger elasticity concerning climate variability than the homogeneous subsurface representation. This difference tends to increase towards the future. Our results suggest that water management in regions with heterogeneous subsurface can expect a higher water availability than estimated by most of the current large-scale simulations, while measures should be taken to prepare for increasingly variable groundwater recharge rates.

Subsurface flow recently triggered the development of taliks under a sub-Arctic road embankment: a prelude to the collapse of the Alaska Highway?

NASA Astrophysics Data System (ADS)

Chen, L.; Fortier, D.; Sliger, M.; McKenzie, J. M.; Murchison, P.

2017-12-01

The Alaska Highway extends over 2200 km between central Alaska, U.S.A. and northern British-Columbia, Canada. This transportation corridor is crucial for the economy of Alaska as it is the only terrestrial link between mainland Alaska and the contiguous United States. Northern British Columbia and southwestern Yukon also greatly benefit from this highway for the transportation of goods and people across this remote corner of Canada. About a quarter of the Alaska Highway is built on permafrost, which is typically ice-rich and at a temperature near the point of thawing. Degradation of the permafrost under the embankment has led to severe structural damages to the highway such as deep longitudinal cracks, extended depressions, potholes and sinkholes. Here we present thermal data from the Beaver Creek experimental road test section in southwestern Yukon. Our study investigates convective heat transfers linked to subsurface water flow under the road embankment based on seven years (2009 to 2016) of thermal monitoring. Observation results demonstrate that snowmelt water infiltration in the spring causes rapid temperature increase of the upper portion of the embankment. Later in the summer, subsurface flow under the highway embankment can lead to step temperature-increase rates, which can be 200 times larger than those via conductive heat transfers. In the fall water trapped under the road significantly delays freeze back of the active layer and contributes to higher permafrost temperature. During the monitoring period, we observed the initiation and growth of taliks along sub-surface flow paths. Positive feedback mechanisms related to water flow through the taliks significantly increased permafrost degradation. Such taliks represent an un-precedent and presumably irreversible thermal state of the highway. Similar terrain conditions which severely threaten the structural integrity of the infrastructure on the short term are numerous along the Alaska Highway corridor.

Interactions among hydraulic conductivity distributions, subsurface topography, and transport thresholds revealed by a multitracer hillslope irrigation experiment

DOE PAGES

Jackson, C. Rhett; Du, Enhao; Klaus, Julian; ...

2016-08-12

Interactions among hydraulic conductivity distributions, subsurface topography, and lateral flow are poorly understood. We applied 407 mm of water and a suite of tracers over 51 h to a 12 by 16.5 m forested hillslope segment to determine interflow thresholds, preferential pathway pore velocities, large-scale conductivities, the time series of event water fractions, and the fate of dissolved nutrients. The 12% hillslope featured loamy sand A and E horizons overlying a sandy clay loam Bt at 1.25 m average depth. Interflow measured from two drains within an interception trench commenced after 131 and 208 mm of irrigation. Cumulative interflow equaledmore » 49% of applied water. Conservative tracer differences between the collection drains indicated differences in flow paths and storages within the plot. Event water fractions rose steadily throughout irrigation, peaking at 50% sixteen h after irrigation ceased. Data implied that tightly held water exchanged with event water throughout the experiment and a substantial portion of preevent water was released from the argillic layer. Surface-applied dye tracers bypassed the matrix, with peak concentrations measured shortly after flow commencement, indicating preferential network conductivities of 864–2240 mm/h, yet no macropore flow was observed. Near steady-state flow conditions indicated average conductivities of 460 mm/h and 2.5 mm/h for topsoils and the Bt horizon, respectively. Low ammonium and phosphorus concentrations in the interflow suggested rapid uptake or sorption, while higher nitrate concentrations suggested more conservative transport. Lastly, these results reveal how hydraulic conductivity variation and subsurface topographic complexity explain otherwise paradoxical solute and flow behaviors.« less

Role of surface and subsurface lateral water flows on summer precipitation in a complex terrain region: A WRF-Hydro case-study for Southern Germany

NASA Astrophysics Data System (ADS)

Rummler, Thomas; Arnault, Joel; Gochis, David; Kunstmann, Harald

2017-04-01

Recent developments in hydrometeorological modeling aim towards more sophisticated treatment of terrestrial hydrologic processes. The standard version of the Weather Research and Forecasting (WRF) model describes terrestrial water transport as a purely vertical process. The hydrologically enhanced version of WRF, namely WRF-Hydro, does account for lateral terrestrial water flows, which allows for a more comprehensive process description of the interdependencies between water- and energy fluxes at the land-atmosphere interface. In this study, WRF and WRF-Hydro are applied to the Bavarian Alpine region in southern Germany, a complex terrain landscape in a relatively humid, mid-latitude climate. Simulation results are validated with gridded and station observation of precipitation, temperature and river discharge. Differences between WRF and WRF-Hydro results are investigated with a joint atmospheric-terrestrial water budget analysis. Changes in the partitioning in (near-) surface runoff and percolation are prominent. However, values for evapotranspiration ET feature only marginal variations, suggesting that soil moisture content is not a limiting factor of ET in this specific region. Simulated precipitation fields during isolated summertime events still show appreciable differences, while differences in large-scale, multi-day rainy periods are less substantial. These differences are mainly related to differences in the moisture in- and outflow terms of the atmospheric water budget induced by the surface and sub-surface lateral redistribution of soil moisture in WRF-Hydro.

Improving Long-term Post-wildfire hydrologic simulations using ParFlow

NASA Astrophysics Data System (ADS)

Lopez, S. R.; Kinoshita, A. M.

2015-12-01

Wildfires alter the natural hydrologic processes within a watershed. After vegetation is burned, the combustion of organic material and debris settles into the soil creating a hydrophobic layer beneath the soil surface with varying degree of thickness and depth. Vegetation regrowth rates vary as a function of radiative exposure, burn severity, and precipitation patterns. Hydrologic models used by the Burned Area Emergency Response (BAER) teams use input data and model calibration constraints that are generally either one-dimensional, empirically-based models, or two-dimensional, conceptually-based models with lumped parameter distributions. These models estimate runoff measurements at the watershed outlet; however, do not provide a distributed hydrologic simulation at each point within the watershed. This work uses ParFlow, a three-dimensional, distributed hydrologic model to (1) correlate burn severity with hydrophobicity, (2) evaluate vegetation recovery rate on water components, and (3) improve flood prediction for managers to help with resource allocation and management operations in burned watersheds. ParFlow is applied to Devil Canyon (43 km2) in San Bernardino, California, which was 97% burned in the 2003 Old Fire. The model set-up uses a 30m-cell size resolution over a 6.7 km by 6.4 km lateral extent. The subsurface reaches 30 m and is assigned a variable cell thickness. Variable subsurface thickness allows users to explicitly consider the degree of recovery throughout the stages of regrowth. Burn severity maps from remotely sensed imagery are used to assign initial hydrophobic layer parameters and thickness. Vegetation regrowth is represented with satellite an Enhanced Vegetation Index. Pre and post-fire hydrologic response is evaluated using runoff measurements at the watershed outlet, and using water component (overland flow, lateral flow, baseflow) measurements.

Comparison of Two Conceptually Different Physically-based Hydrological Models - Looking Beyond Streamflows

NASA Astrophysics Data System (ADS)

Rousseau, A. N.; Álvarez; Yu, X.; Savary, S.; Duffy, C.

2015-12-01

Most physically-based hydrological models simulate to various extents the relevant watershed processes occurring at different spatiotemporal scales. These models use different physical domain representations (e.g., hydrological response units, discretized control volumes) and numerical solution techniques (e.g., finite difference method, finite element method) as well as a variety of approximations for representing the physical processes. Despite the fact that several models have been developed so far, very few inter-comparison studies have been conducted to check beyond streamflows whether different modeling approaches could simulate in a similar fashion the other processes at the watershed scale. In this study, PIHM (Qu and Duffy, 2007), a fully coupled, distributed model, and HYDROTEL (Fortin et al., 2001; Turcotte et al., 2003, 2007), a pseudo-coupled, semi-distributed model, were compared to check whether the models could corroborate observed streamflows while equally representing other processes as well such as evapotranspiration, snow accumulation/melt or infiltration, etc. For this study, the Young Womans Creek watershed, PA, was used to compare: streamflows (channel routing), actual evapotranspiration, snow water equivalent (snow accumulation and melt), infiltration, recharge, shallow water depth above the soil surface (surface flow), lateral flow into the river (surface and subsurface flow) and height of the saturated soil column (subsurface flow). Despite a lack of observed data for contrasting most of the simulated processes, it can be said that the two models can be used as simulation tools for streamflows, actual evapotranspiration, infiltration, lateral flows into the river, and height of the saturated soil column. However, each process presents particular differences as a result of the physical parameters and the modeling approaches used by each model. Potentially, these differences should be object of further analyses to definitively confirm or reject modeling hypotheses.

Topographically Driven Lateral Water Fluxes and Their Influence on Carbon Assimilation of a Black Spruce Ecosystem.

NASA Astrophysics Data System (ADS)

Govind, A.; Chen, J. M.; Margolis, H.; Bernier, P. Y.

2006-12-01

Current estimates of ecophysiological indicators overlook the effects of topographically-driven lateral flow of soil water. We hypothesize that topographically driven lateral water flows over the landscape have significant influence on the terrestrial carbon cycle. To this end, we simulated the hydrological controls on carbon cycle processes in a black spruce forest in central Quebec, Canada, using the Boreal Ecosystem Productivity Simulator (BEPS) at a daily time step. We accounted for lateral surface and subsurface flows in BEPS by incorporating a distributed, process-oriented hydrological procedure. The results show that modeled dynamics of ecophysiological processes such as evapotranspiration (ET) and photosynthesis (GPP) are consistent with the spatial variation of land cover, topography, soil texture, and leaf area index. Simulated ET and GPP averaged within the footprint of an eddy covariance tower in the watershed agree well with flux measurements with R2=0.77 and 0.83 for ET and GPP, respectively. For ET simulation, much of the remaining discrepancies are found in the winter when the model underestimates snow sublimation. For GPP, there is an underestimation in the fall coinciding with a mid growing season drought, showing the high sensitivity of the model to the soil water status. The key processes controlling primary production were hydraulic limitations for water transfer from soil, roots, stems and leaves through stomatal conductance. Therefore, a further understanding of soil water dynamics is warranted. Comparison with the soil water content of the footprint- averaged unsaturated zone showed that the model captured the annual trend. We also simulated the variations in the water table as well as the mid growing season drought, with a reasonable accuracy(R2=0.68). The foot print average water budget reveals that the annual precipitation of 835mm is partitioned into 282mm of ET, 541 mm of subsurface runoff, and 6 mm of storage change. To test the influence of topographically driven lateral water flow on the carbon cycle, we made three hydrological modeling scenarios viz. 1) explicit hydrological simulation including lateral water routing, 2) bucket model with implicit runoff calculations and 3) a control run, where the lateral water flow was turned off in the model. Bucket model overestimated GPP as much as 25% as opposed to explicit simulations because there was no topographical constrain on runoff. Flat areas dominated with mineral soils shows the highest overestimation because of an increase in stomatal conductance. Control simulation, on the other hand, underestimated GPP as much as 15% as opposed to explicit routing because of rapid soil saturation, which decreases stomatal conductance. These results suggest that lateral water flow does play a significant role in the terrestrial carbon cycle and should be accounted for in ecological models. For details please see http://ajit.govind.googlepages.com/agu2006

Model reduction of the numerical analysis of Low Impact Developments techniques

NASA Astrophysics Data System (ADS)

Brunetti, Giuseppe; Šimůnek, Jirka; Wöhling, Thomas; Piro, Patrizia

2017-04-01

Mechanistic models have proven to be accurate and reliable tools for the numerical analysis of the hydrological behavior of Low Impact Development (LIDs) techniques. However, their widespread adoption is limited by their complexity and computational cost. Recent studies have tried to address this issue by investigating the application of new techniques, such as surrogate-based modeling. However, current results are still limited and fragmented. One of such approaches, the Model Order Reduction (MOR) technique, can represent a valuable tool for reducing the computational complexity of a numerical problems by computing an approximation of the original model. While this technique has been extensively used in water-related problems, no studies have evaluated its use in LIDs modeling. Thus, the main aim of this study is to apply the MOR technique for the development of a reduced order model (ROM) for the numerical analysis of the hydrologic behavior of LIDs, in particular green roofs. The model should be able to correctly reproduce all the hydrological processes of a green roof while reducing the computational cost. The proposed model decouples the subsurface water dynamic of a green roof in a) one-dimensional (1D) vertical flow through a green roof itself and b) one-dimensional saturated lateral flow along the impervious rooftop. The green roof is horizontally discretized in N elements. Each element represents a vertical domain, which can have different properties or boundary conditions. The 1D Richards equation is used to simulate flow in the substrate and drainage layers. Simulated outflow from the vertical domain is used as a recharge term for saturated lateral flow, which is described using the kinematic wave approximation of the Boussinesq equation. The proposed model has been compared with the mechanistic model HYDRUS-2D, which numerically solves the Richards equation for the whole domain. The HYDRUS-1D code has been used for the description of vertical flow, while a Finite Volume Scheme has been adopted for lateral flow. Two scenarios involving flat and steep green roofs were analyzed. Results confirmed the accuracy of the reduced order model, which was able to reproduce both subsurface outflow and the moisture distribution in the green roof, significantly reducing the computational cost.

Relative Impacts of Low Permeability Subsurface Deposits on Recharge Basin Infiltration Rates

NASA Astrophysics Data System (ADS)

Oconnell, P.; Becker, M.; Pham, C.; Rodriguez, G.; Hutchinson, A.; Plumlee, M.

2017-12-01

Artificial recharge of aquifers through spreading basins has become an important component of water management in semi-arid climates. The rate at which water can be recharged in these basins is limited by the natural vertical permeability of the underlying deposits which may be highly variable both laterally and vertically. To help understand hydrostratigraphic controls on recharge, a newly constructed basin was surveyed and instrumented. Prior to flooding the basin, lithology was characterized by shallow hand coring, direct push coring, ground penetrating radar, and electrical resistivity. After flooding, recharge was monitored through piezometers, electrical resistivity, and a network of fiber optic distributed temperature sensing (DTS). The DTS network used temperature as a tracer to measure infiltration rate on 25 cm intervals both laterally and vertically. Several hundred paired DTS time series datasets (from fiber optic cables located at 0 and 0.5 meters below ground surface) were processed with the cross-wavelet transform (XWT) to calculate spatially and temporally continuous infiltration rates, which can be interpolated and animated to visualize heterogeneity. Time series data from 8-meter deep, vertically oriented DTS cables reveal depth intervals where infiltration rates vary. Inverted resistivity sections from repeated dipole-dipole surveys along the sidewall of a spreading basin exhibit a positive correlation with the distribution of relatively high and low infiltration rates, indicating zones of preferential downward (efficient) and lateral (inefficient) flow, respectively. In contrast to other monitored basins, no perching was observed in the vertically oriented DTS cables. The variation in recharge across the basin and the appearance of subsurface lateral flow can be explained in context of the alluvial depositional environment.

Subsurface drainage processes and management impacts

Treesearch

Elizabeth T. Keppeler; David Brown

1998-01-01

Storm-induced streamflow in forested upland watersheds is linked to rainfall by transient, variably saturated flow through several different flow paths. In the absence of exposed bedrock, shallow flow-restrictive layers, or compacted soil surfaces, virtually all of the infiltrated rainfall reaches the stream as subsurface flow. Subsurface runoff can occur within...

An optimization model to design and manage subsurface drip irrigation system for alfalfa

NASA Astrophysics Data System (ADS)

Kandelous, M.; Kamai, T.; Vrugt, J. A.; Simunek, J.; Hanson, B.; Hopmans, J. W.

2010-12-01

Subsurface drip irrigation (SDI) is one of the most efficient and cost-effective methods for watering alfalfa plants. Lateral installation depth and distance, emitter discharge, and irrigation time and frequency of SDI, in addition to soil and climatic conditions affect alfalfa’s root water uptake and yield. Here we use a multi-objective optimization approach to find optimal SDI strategies. Our approach uses the AMALGAM evolutionary search method, in combination with the HYDRUS-2D unsaturated flow model to maximize water uptake by alfalfa’s plant roots, and minimize loss of irrigation and drainage water to the atmosphere or groundwater. We use a variety of different objective functions to analyze SDI. These criteria include the lateral installation depth and distance, the lateral discharge, irrigation duration, and irrigation frequency. Our framework includes explicit recognition of the soil moisture status during the simulation period to make sure that the top soil is dry for harvesting during the growing season. Initial results show a wide spectrum of optimized SDI strategies for different root distributions, soil textures and climate conditions. The developed tool should be useful in helping farmers optimize their irrigation strategy and design.

Imaging and locating paleo-channels using geophysical data from meandering system of the Mun River, Khorat Plateau, Northeastern Thailand

NASA Astrophysics Data System (ADS)

Nimnate, P.; Thitimakorn, T.; Choowong, M.; Hisada, K.

2017-12-01

The Khorat Plateau from northeast Thailand, the upstream part of the Mun River flows through clastic sedimentary rocks. A massive amount of sand was transported. We aimed to understand the evolution of fluvial system and to discuss the advantages of two shallow geophysical methods for describing subsurface morphology of modern and paleo-channels. We applied Electrical Resistivity Tomography (ERT) and Ground Penetrating Radar (GPR) to characterize the lateral, vertical morphological and sedimentary structures of paleo-channels, floodplain and recent point bars. Both methods were interpreted together with on-sites boreholes to describe the physical properties of subsurface sediments. As a result, we concluded that four radar reflection patterns including reflection free, shingled, inclined and hummocky reflections were appropriated to apply as criteria to characterize lateral accretion, the meandering rivers with channel-filled sequence and floodplain were detected from ERT profiles. The changes in resistivity correspond well with differences in particle size and show relationship with ERT lithological classes. Clay, silt, sand, loam and bedrock were classified by the resistivity data. Geometry of paleo-channel embayment and lithological differences can be detected by ERT, whereas GPR provides detail subsurface facies for describing point bar sand deposit better than ERT.

Characteristics of Nitrogen Loss through Surface-Subsurface Flow on Red Soil Slopes of Southeast China

NASA Astrophysics Data System (ADS)

Zheng, Haijin; Liu, Zhao; Zuo, Jichao; Wang, Lingyun; Nie, Xiaofei

2017-12-01

Soil nitrogen (N) loss related to surface flow and subsurface flow (including interflow and groundwater flow) from slope lands is a global issue. A lysimetric experiment with three types of land cover (grass cover, GC; litter cover, LC; and bare land, BL) were carried out on a red soil slope land in southeast China. Total Nitrogen (TN) loss through surface flow, interflow and groundwater flow was observed under 28 natural precipitation events from 2015 to 2016. TN concentrations from subsurface flow on BL and LC plots were, on average, 2.7-8.2 and 1.5-4.4 times greater than TN concentrations from surface flow, respectively; the average concentration of TN from subsurface flow on GC was about 36-56% of that recorded from surface flow. Surface flow, interflow and groundwater flow contributed 0-15, 2-9 and 76-96%, respectively, of loss load of TN. Compared with BL, GC and LC intercepted 83-86% of TN loss through surface runoff; GC intercepted 95% of TN loss through subsurface flow while TN loss through subsurface flow on LC is 2.3 times larger than that on BL. In conclusion, subsurface flow especially groundwater flow is the dominant hydrological rout for N loss that is usually underestimated. Grass cover has the high retention of N runoff loss while litter mulch will increase N leaching loss. These findings provide scientific support to control N runoff loss from the red soil slope lands by using suitable vegetation cover and mulching techniques.

A Generalized Subsurface Flow Parameterization Considering Subgrid Spatial Variability of Recharge and Topography

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

Huang, Maoyi; Liang, Xu; Leung, Lai R.

2008-12-05

Subsurface flow is an important hydrologic process and a key component of the water budget, especially in humid regions. In this study, a new subsurface flow formulation is developed that incorporates spatial variability of both topography and recharge. It is shown through theoretical derivation and case studies that the power law and exponential subsurface flow parameterizations and the parameterization proposed by Woods et al.[1997] are all special cases of the new formulation. The subsurface flows calculated using the new formulation compare well with values derived from observations at the Tulpehocken Creek and Walnut Creek watersheds. Sensitivity studies show that whenmore » the spatial variability of topography or recharge, or both is increased, the subsurface flows increase at the two aforementioned sites and the Maimai hillslope. This is likely due to enhancement of interactions between the groundwater table and the land surface that reduce the flow path. An important conclusion of this study is that the spatial variability of recharge alone, and/or in combination with the spatial variability of topography can substantially alter the behaviors of subsurface flows. This suggests that in macroscale hydrologic models or land surface models, subgrid variations of recharge and topography can make significant contributions to the grid mean subsurface flow and must be accounted for in regions with large surface heterogeneity. This is particularly true for regions with humid climate and relatively shallow groundwater table where the combined impacts of spatial variability of recharge and topography are shown to be more important. For regions with arid climate and relatively deep groundwater table, simpler formulations, especially the power law, for subsurface flow can work well, and the impacts of subgrid variations of recharge and topography may be ignored.« less

Stratigraphic and structural characterization of the OU-1 area at the former George Air Force Base, Adelanto, Southern California

USGS Publications Warehouse

Catchings, R.D.; Gandhok, G.; Goldman, M.R.

2001-01-01

The former George Air Force Base (GAFB), now known as the Southern California Logistics Airport (SCLA), is located in the town of Adelanto, approximately 100 km northeast of Los Angeles, California (Fig. 1). In this report, we present acquisition parameters, data, and interpretations of seismic images that were acquired in the OU-1 area of GAFB during July 1999 (Fig. 2). GAFB is scheduled for conversion to civilian use, however, during its years as an Air Force base, trichlorethylene (TCE) was apparently introduced into the subsurface as a result of spills during normal aircraft maintenance operations. To comply with congressional directives, TCE contaminant removal has been ongoing since the early-tomid 1990s. However, only a small percentage of the TCE believed to have been introduced into the subsurface has been recovered, due largely to difficulty in locating the TCE within the subsurface. Because TCE migrates within the subsurface by ground water movement, attempts to locate the TCE contaminants in the subsurface have employed an array of ground-water monitoring and extraction wells. These wells primarily sample within a shallow-depth (~40 m) aquifer system. Cores obtained from the monitoring and extraction wells indicate that the aquifer, which is composed of sand and gravel channels, is bounded by aquitards composed largely of clay and other fine-grained sediments. Based on well logs, the aquifer is about 3 to 5 m thick along the seismic profiles. A more thorough understanding of the lateral variations in the depth and thickness of the aquifer system may be a key to finding and removing the remaining TCE. However, due to its complex depositional and tectonic history, the structural and stratigraphic sequences are not easily characterized. An indication of the complex nature of the structure and stratigraphy is the appreciable variation in stratigraphic sequences observed in some monitoring wells that are only a few tens of meters apart. To better characterize the shallow (upper 100 m) stratigraphy beneath GAFB, the US Environmental Protection Agency (USEPA) contracted the US Geological Survey (USGS) to acquire three seismic reflection/refraction profiles within an area known as Operational Unit #1 (OU-1). The principal objective of the seismic survey was to laterally characterize the subsurface with respect to structure and stratigraphy. In particular, we desired to (1) laterally “map” stratigraphic units (particularly aquifer layers) that were previously identified in monitoring wells within the OU-1 area and (2) identify structures, such as faults and folds, that affect the movement of ground water. Knowledge of lateral variations in stratigraphic units and structures that may affect those units is useful in constructing ground-water flow models, which aid in identifying possible TCE migration paths within the subsurface. Stratigraphic and structural characterization may also be useful in identifying surface locations and target depths for future wells (Catchings et al., 1996). Proper siting of wells is important because a welldefined aquifer is apparently not present in all locations at GAFB, as indicated by lithologic logs from existing wells (Montgomery Watson, 1995). Proper depth placement of monitoring and extraction wells is important because wells that are too shallow will not sample within the aquifer, and wells that are too deep risk puncturing the aquitard and allowing contaminants to flow to deeper levels.

Impact of transient stream flow on water exchange and reactions in the hyporheic zone of an in-stream gravel bar

NASA Astrophysics Data System (ADS)

Trauth, Nico; Schmidt, Christian; Fleckenstein, Jan H.

2015-04-01

Groundwater-surface water exchange is an important process that can facilitate the degradation of critical substances like nitrogen-species and contaminants, supporting a healthy status of the aquatic ecosystem. In our study, we simulate water exchange, solute transport and reactions within a natural in-stream gravel bar using a coupled surface and subsurface numerical model. Stream water flow is simulated by computational fluid dynamics software that provides hydraulic head distributions at the streambed, which are used as an upper boundary condition for a groundwater model. In the groundwater model water exchange, solute transport, aerobic respiration and denitrification in the subsurface are simulated. Ambient groundwater flow is introduced by lateral upstream and downstream hydraulic head boundaries that generate neutral, losing or gaining stream conditions. Stream water transports dissolved oxygen, organic carbon (as the dominant electron donor) and nitrate into the subsurface, whereas an additional nitrate source exists in the ambient groundwater. Scenarios of stream flow events varying in duration and stream stage are simulated and compared with steady state scenarios with respect to water fluxes, residence times and the solute turn-over rates. Results show, that water exchange and solute turn-over rates highly depend on the interplay between event characteristics and ambient groundwater levels. For scenarios, where the stream flow event shifts the hydraulic system to a net-neutral hydraulic gradient between the average stream stage and the ambient groundwater level (minimal exchange between ground- and surface water), solute consumption is higher, compared to the steady losing or gaining case. In contrast, events that induce strong losing conditions lead to a lower potential of solute consumption.

Modeling GPR data to interpret porosity and DNAPL saturations for calibration of a 3-D multiphase flow simulation

USGS Publications Warehouse

Sneddon, Kristen W.; Powers, Michael H.; Johnson, Raymond H.; Poeter, Eileen P.

2002-01-01

Dense nonaqueous phase liquids (DNAPLs) are a pervasive and persistent category of groundwater contamination. In an effort to better understand their unique subsurface behavior, a controlled and carefully monitored injection of PCE (perchloroethylene), a typical DNAPL, was performed in conjunction with the University of Waterloo at Canadian Forces Base Borden in 1991. Of the various geophysical methods used to monitor the migration of injected PCE, the U.S. Geological Survey collected 500-MHz ground penetrating radar (GPR) data. These data are used in determining calibration parameters for a multiphase flow simulation. GPR data were acquired over time on a fixed two-dimensional surficial grid as the DNAPL was injected into the subsurface. Emphasis is on the method of determining DNAPL saturation values from this time-lapse GPR data set. Interactive full-waveform GPR modeling of regularized field traces resolves relative dielectric permittivity versus depth profiles for pre-injection and later-time data. Modeled values are end members in recursive calculations of the Bruggeman-Hanai-Sen (BHS) mixing formula, yielding interpreted pre-injection porosity and post-injection DNAPL saturation values. The resulting interpreted physical properties of porosity and DNAPL saturation of the Borden test cell, defined on a grid spacing of 50 cm with 1-cm depth resolution, are used as observations for calibration of a 3-D multiphase flow simulation. Calculated values of DNAPL saturation in the subsurface at 14 and 22 hours after the start of injection, from both the GPR and the multiphase flow modeling, are interpolated volumetrically and presented for visual comparison.

Solute transport with time-variable flow paths during upward and downward flux in a heterogeneous unsaturated porous medium

NASA Astrophysics Data System (ADS)

Cremer, Clemens; Neuweiler, Insa; Bechtold, Michel; Vanderborght, Jan

2014-05-01

To acquire knowledge of solute transport through the unsaturated zone in the shallow subsurface is decisive to assess groundwater quality, nutrient cycling or to plan remediation strategies. The shallow subsurface is characterized by structural heterogeneity and strongly influenced by atmospheric conditions. This leads to changing flow directions, strong temporal changes in saturation and heterogeneous water fluxes during infiltration and evaporation events. Recent studies (e.g. Lehmann and Or, 2009; Bechtold et al.,2011) demonstrated the importance of lateral flow and solute transport during evaporation conditions (upward flux). The heterogeneous structure in these studies was constructed using two types of sand with strong material contrasts and arranged in parallel with a vertical orientation. Lateral transport and redistribution of solute from coarse to fine media was observed deeper in the soil column and from fine to coarse close to the soil surface. However, if boundary conditions are reversed due to precipitation, the flow field is not necessarily reversed in the same manner, resulting in entirely different transport patterns for downward and upward flow. Therefore, considering net-flow rates alone is misleading when describing transport under those conditions. In this contribution we analyze transport of a solute in the shallow subsurface to assess effects resulting from the temporal change of heterogeneous soil structures due to dynamic flow conditions. Two-dimensional numerical simulations of unsaturated flow and transport are conducted using a coupled finite volume and random walk particle tracking algorithm to quantify solute transport and leaching rates. Following previous studies (Lehmann and Or, 2009; Bechtold et al., 2011), the chosen domain is composed of two materials, coarse and fine sand, arranged in parallel with a vertical orientation. Hence, one sharp interface of strong material heterogeneity is induced. During evaporation both sands are assumed to stay under liquid-flow dominated evaporation conditions ("stage 1"). Simulations considering dynamic (infiltration-evaporation) and steady (solely infiltration) boundary conditions are carried out. The influence of dynamic boundary conditions (intensity and duration of precipitation and evaporation events) is examined in a multitude of simulations. If flow rates smaller than the saturated hydraulic conductivity of both materials are chosen to be applied as boundary condition, simulation results indicate that the flow field within the domain is exactly reversed. However, if applied flow rates exceed the saturated hydraulic conductivity of one material, the flow field is not just reversed, but different flow paths during downward and upward flow are observed. Results show the tendency of faster solute leaching under dynamic boundary conditions compared to steady infiltration conditions with the same net-infiltration rate. We use a double domain transport method as an upscaled model to reproduce vertically averaged concentration profiles with net flux only and compare the model parameters for information about flow dynamics and soil heterogeneity.

Simulated ground-water flow and sources of water in the Killbuck Creek Valley near Wooster, Wayne County, Ohio

USGS Publications Warehouse

Breen, K.J.; Kontis, A.L.; Rowe, G.L.; Haefner, R.J.

1995-01-01

The stratified-drift aquifer in the 3,000-ft (feet)-wide and 100-ft-deep buried valley of Killbuck Creek near Wooster in northeastern Ohio was studied. The stratified drift with adjacent sandstone and shale bedrock produce a system of ground-water flow representative of the western part of the glaciated north-eastern United States. The stratified-drift aquifer is an excellent source of water for municipal and industrial wells. The aquifer is recharged locally by water from precipitation on the valley floor and uplands, by infiltration from streams, and by lateral flow to the valley from the uplands. As a result, the aquifer is vulnerable to surface or subsurface spills of contaminants in the valley or the adjacent uplands. Quality of water in the stratified drift is affected by influx of water from bedrock lateral to or beneath the valley. This influx is controlled, in part, by the pumping stress placed on the stratified-drift aquifer. Hydrogeologic and aqueous-geochemical data were analyzed to establish the framework necessary for stead-state and transient simulations of ground-water flow in stratified drift and bedrock with a three-layer ground-water-flow model. A new model routine, the Variable-Recharge procedure, was developed to simulate areal recharge and the contribution of the uplands to the drift system. This procedure allows for water applied to land surface to infiltrate or to be rejected. Rejected recharge and ground water discharged when the water table is at land surface form surface runoff-this excess upland water can be redirected as runoff to other parts of the model. Infiltration of streamwater, areal recharge to uplands and valley, and lateral subsurface flow from the uplands to the valley are sources of water to the stratufued0druft aquifer. Water is removed from the stratified-drift aquifer at Wooster primarily by production wells pumping at a rate of approximately 8.5 ft3/s (cubic feet per second). The ground-water budget resulting from two types of simulations of ground-water flow in this study indicates the primary sources of water to the wells are recharge at or near land surface and lateral subsurface flow from the shale and sandstone bedrock. Components of recharge at land surface include induced infiltration from streams, precipitation on the valley floor, and infiltration of unchanneled upland runoff that reaches the valley floor. The steady-state simulation was designed to represent conditions during the fall of 1984. The transient simulation was designed to represent an 11-day snowmelt event, 23 February to 5 March 1985, that caused water levels to rise significantly throughout the valley. Areal recharge to the valley and flow from the uplands to the valley were determined through the Variable-Recharge procedure. The total steady-state recharge to the valley was 12.5 ft3/s. Upland sources, areal valley recharge, and induced infiltration from Killnuck Creek accounted for 63, 23, and 8 percent, respectively, of the valley recharge. An analysis of the simulated vertical flow to the buried stratified drift through surficial slit, clay, and fine sand indicates that about 75 percent of the total recharge to the buried deposits is the sum of areally extensive, relatively small flows less than about 0.01 ft? /s per model node), whereas about 25 percent of the recharge results from a really restricted, relatively large flows (greater than about 0.01 ft? /s per model node). The large-magnitude flows are located primarily beneath Clear and Little Killbuck Creeks where seepage provides abundant recharge and the surficial sediments grade into coarser alluvial-fan deposits. Chemical and isotopic studies of ground water and streamwater combined with measurements of stream infiltration provide independent support for the conclusions derived from computer simulation of ground-water flow. In addition, the chemical and isotopic studies helped quantity the rate and pathways of infiltrating water from

Simulating land-atmosphere feedbacks and response to widespread forest disturbance: The role of lower boundary configuration and dynamic water table in meteorological modeling

NASA Astrophysics Data System (ADS)

Forrester, M.; Maxwell, R. M.; Bearup, L. A.; Gochis, D.

2017-12-01

Numerical meteorological models are frequently used to diagnose land-atmosphere interactions and predict large-scale response to extreme or hazardous events, including widespread land disturbance or perturbations to near-surface moisture. However, few atmospheric modeling platforms consider the impact that dynamic groundwater storage, specifically 3D subsurface flow, has on land-atmosphere interactions. In this study, we use the Weather Research and Forecasting (WRF) mesoscale meteorological model to identify ecohydrologic and land-atmosphere feedbacks to disturbance by the mountain pine beetle (MPB) over the Colorado Headwaters region. Disturbance simulations are applied to WRF with various lower boundary configurations: Including default Noah land surface model soil moisture representation; a version of WRF coupled to ParFlow (PF), an integrated groundwater-surface water model that resolves variably saturated flow in the subsurface; and WRF coupled to PF in a static water table version, simulating only vertical and no lateral subsurface flow. Our results agree with previous literature showing MPB-induced reductions in canopy transpiration in all lower boundary scenarios, as well as energy repartitioning, higher water tables, and higher planetary boundary layer over infested regions. Simulations show that expanding from local to watershed scale results in significant damping of MPB signal as unforested and unimpacted regions are added; and, while deforestation appears to have secondary feedbacks to planetary boundary layer and convection, these slight perturbations to cumulative summer precipitation are insignificant in the context of ensemble methodologies. Notably, the results suggest that groundwater representation in atmospheric modeling affects the response intensity of a land disturbance event. In the WRF-PF case, energy and atmospheric processes are more sensitive to disturbance in regions with higher water tables. Also, when dynamic subsurface hydrology is removed, WRF simulates a greater response to MPB at the land-atmosphere interface, including greater changes to daytime skin temperature, Bowen ratio and near-surface humidity. These findings highlight lower boundary representations in computational meteorology and numerical land-atmosphere modeling.

Effect of subalpine canopy removal on snowpack, soil solution, and nutrient export, Fraser Experimental Forest, CO

USGS Publications Warehouse

Stottlemyer, R.; Troendle, C.A.

1999-01-01

Research on the effects of vegetation manipulation on snowpack, soil water, and streamwater chemistry and flux has been underway at the Fraser Experimental Forest (FEF), CO, since 1982. Greater than 95% of FEF snowmelt passes through watersheds as subsurface flow where soil processes significantly alter meltwater chemistry. To better understand the mechanisms accounting for annual variation in watershed streamwater ion concentration and flux with snowmelt, we studied subsurface water flow, its ion concentration, and flux in conterminous forested and clear cut plots. Repetitive patterns in subsurface flow and chemistry were apparent. Control plot subsurface flow chemistry had the highest ion concentrations in late winter and fall. When shallow subsurface flow occurred, its Ca2+, SO42-, and HCO3- concentrations were lower and K+ higher than deep flow. The percentage of Ca2+, NO3-, SO42-, and HCO3- flux in shallow depths was less and K+ slightly greater than the percentage of total flow. Canopy removal increased precipitation reaching the forest floor by about 40%, increased peak snowpack water equivalent (SWE) > 35%, increased the average snowpack Ca2+, NO3-, and NH4+ content, reduced the snowpack K+ content, and increased the runoff four-fold. Clear cutting doubled the percentage of subsurface flow at shallow depths, and increased K+ concentration in shallow subsurface flow and NO3- concentrations in both shallow and deep flow. The percentage change in total Ca2+, SO42-, and HCO3- flux in shallow depths was less than the change in water flux, while that of K+ and NO3- flux was greater. Relative to the control, in the clear cut the percentage of total Ca2+ flux at shallow depths increased from 5 to 12%, SO42- 5.4 to 12%, HCO3- from 5.6 to 8.7%, K+ from 6 to 35%, and NO3- from 2.7 to 17%. The increases in Ca2+ and SO42- flux were proportional to the increase in water flux, the flux of HCO3- increased proportionally less than water flux, and NO3- and K+ were proportionally greater than water flux. Increased subsurface flow accounted for most of the increase in non-limiting nutrient loss. For limiting nutrients, loss of plant uptake and increased shallow subsurface flow accounted for the greater loss. Seasonal ion concentration patterns in streamwater and subsurface flow were similar.Research on the effects of vegetation manipulation on snowpack, soil water, and streamwater chemistry and flux has been underway at the Fraser Experimental Forest (FEF), CO, since 1982. Greater than 95% of FEF snowmelt passes through watersheds as subsurface flow where soil processes significantly alter meltwater chemistry. To better understand the mechanisms accounting for annual variation in watershed streamwater ion concentration and flux with snowmelt, we studied subsurface water flow, its ion concentration, and flux in conterminous forested and clear cut plots. Repetitive patterns in subsurface flow and chemistry were apparent. Control plot subsurface flow chemistry had the highest ion concentrations in late winter and fall. When shallow subsurface flow occurred, its Ca2+, SO42-, and HCO3- concentrations were lower and K+ higher than deep flow. The percentage of Ca2+, NO3-, SO42-, and HCO3- flux in shallow depths was less and K+ slightly greater than the percentage of total flow. Canopy removal increased precipitation reaching the forest floor by about 40%, increased peak snowpack water equivalent (SWE) > 35%, increased the average snowpack Ca2+, NO3-, and NH4+ content, reduced the snowpack K+ content, and increased the runoff four-fold. Clear cutting doubled the percentage of subsurface flow at shallow depths, and increased K+ concentration in shallow subsurface flow and NO3- concentrations in both shallow and deep flow. The percentage change in total Ca2+, SO42-, and HCO3- flux in shallow depths was less than the change in water flux, while that of K+ and NO3- flux was greater. Relative to the control, in the clear cut the percentage of total Ca

Need to improve SWMM's subsurface flow routing algorithm for green infrastructure modeling

EPA Science Inventory

SWMM can simulate various subsurface flows, including groundwater (GW) release from a subcatchment to a node, percolation out of storage units and low impact development (LID) controls, and rainfall derived inflow and infiltration (RDII) at a node. Originally, the subsurface flow...

Large-scale assessment of present day and future groundwater recharge and its sensitivity to climate variability in Europe's karst regions

NASA Astrophysics Data System (ADS)

Hartmann, A. J.; Gleeson, T. P.; Wagener, T.; Wada, Y.

2016-12-01

Karst aquifers in Europe are an important source of fresh water contributing up to half of the total drinking water supply in some countries. Karstic groundwater recharge is one of the most important components of the water balance of karst systems as it feeds the karst aquifers. Presently available large-scale hydrological models do not consider karst heterogeneity adequately. Projections of current and potential future groundwater recharge of Europe's karst aquifers are therefore unclear. In this study we compare simulations of present (1991-2010) and future (2080-2099) recharge using two different models to simulate groundwater recharge processes. One model includes karst processes (subsurface heterogeneity, lateral flow and concentrated recharge), while the other is based on the conceptual understanding of common hydrological systems (homogeneous subsurface, saturation excess overland flow). Both models are driven by the bias-corrected 5 GCMs of the ISI-MIP project (RCP8.5). To further assess sensitivity of groundwater recharge to climate variability, we calculate the elasticity of recharge rates to annual precipitation, temperature and average intensity of rainfall events, which is the median change of recharge that corresponds to the median change of these climate variables within the present and future time period, respectively. Our model comparison shows that karst regions over Europe have enhanced recharge rates with greater inter-annual variability compared to those with more homogenous subsurface properties. Furthermore, the heterogeneous representation shows stronger elasticity concerning climate variability than the homogeneous subsurface representation. This difference tends to increase towards the future. Our results suggest that water management in regions with heterogeneous subsurface can expect a higher water availability than estimated by most of the current large-scale simulations, while measures should be taken to prepare for increasingly variable groundwater recharge rates.

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

Qiusheng, Y., E-mail: qsyan@gdut.edu.cn; Senkai, C., E-mail: senkite@sina.com; Jisheng, P., E-mail: panjisheng@gdut.edu.cn

Different machining processes were used in the single crystal SiC wafer machining. SEM was used to observe the surface morphology and a cross-sectional cleavages microscopy method was used for subsurface cracks detection. Surface and subsurface cracks characteristics of single crystal SiC wafer in abrasive machining were analysed. The results show that the surface and subsurface cracks system of single crystal SiC wafer in abrasive machining including radial crack, lateral crack and the median crack. In lapping process, material removal is dominated by brittle removal. Lots of chipping pits were found on the lapping surface. With the particle size becomes smaller,more » the surface roughness and subsurface crack depth decreases. When the particle size was changed to 1.5µm, the surface roughness Ra was reduced to 24.0nm and the maximum subsurface crack was 1.2µm. The efficiency of grinding is higher than lapping. Plastic removal can be achieved by changing the process parameters. Material removal was mostly in brittle fracture when grinding with 325# diamond wheel. Plow scratches and chipping pits were found on the ground surface. The surface roughness Ra was 17.7nm and maximum subsurface crack depth was 5.8 µm. When grinding with 8000# diamond wheel, the material removal was in plastic flow. Plastic scratches were found on the surface. A smooth surface of roughness Ra 2.5nm without any subsurface cracks was obtained. Atomic scale removal was possible in cluster magnetorheological finishing with diamond abrasive size of 0.5 µm. A super smooth surface eventually obtained with a roughness of Ra 0.4nm without any subsurface crack.« less

Real rock-microfluidic flow cell: A test bed for real-time in situ analysis of flow, transport, and reaction in a subsurface reactive transport environment.

PubMed

Singh, Rajveer; Sivaguru, Mayandi; Fried, Glenn A; Fouke, Bruce W; Sanford, Robert A; Carrera, Martin; Werth, Charles J

2017-09-01

Physical, chemical, and biological interactions between groundwater and sedimentary rock directly control the fundamental subsurface properties such as porosity, permeability, and flow. This is true for a variety of subsurface scenarios, ranging from shallow groundwater aquifers to deeply buried hydrocarbon reservoirs. Microfluidic flow cells are now commonly being used to study these processes at the pore scale in simplified pore structures meant to mimic subsurface reservoirs. However, these micromodels are typically fabricated from glass, silicon, or polydimethylsiloxane (PDMS), and are therefore incapable of replicating the geochemical reactivity and complex three-dimensional pore networks present in subsurface lithologies. To address these limitations, we developed a new microfluidic experimental test bed, herein called the Real Rock-Microfluidic Flow Cell (RR-MFC). A porous 500μm-thick real rock sample of the Clair Group sandstone from a subsurface hydrocarbon reservoir of the North Sea was prepared and mounted inside a PDMS microfluidic channel, creating a dynamic flow-through experimental platform for real-time tracking of subsurface reactive transport. Transmitted and reflected microscopy, cathodoluminescence microscopy, Raman spectroscopy, and confocal laser microscopy techniques were used to (1) determine the mineralogy, geochemistry, and pore networks within the sandstone inserted in the RR-MFC, (2) analyze non-reactive tracer breakthrough in two- and (depth-limited) three-dimensions, and (3) characterize multiphase flow. The RR-MFC is the first microfluidic experimental platform that allows direct visualization of flow and transport in the pore space of a real subsurface reservoir rock sample, and holds potential to advance our understandings of reactive transport and other subsurface processes relevant to pollutant transport and cleanup in groundwater, as well as energy recovery. Copyright © 2017 Elsevier B.V. All rights reserved.

A field study of colloid transport in surface and subsurface flows

NASA Astrophysics Data System (ADS)

Zhang, Wei; Tang, Xiang-Yu; Xian, Qing-Song; Weisbrod, Noam; Yang, Jae E.; Wang, Hong-Lan

2016-11-01

Colloids have been recognized to enhance the migration of strongly-sorbing contaminants. However, few field investigations have examined combined colloid transport via surface runoff and subsurface flows. In a headwater catchment of the upper Yangtze River, a 6 m (L) by 4 m (W) sloping (6°) farmland plot was built by cement walls to form no-flow side boundaries. The plot was monitored in the summer of 2014 for the release and transport of natural colloids via surface runoff and subsurface flows (i.e., the interflow from the soil-mudrock interface and fracture flow from the mudrock-sandstone interface) in response to rain events. The water sources of the subsurface flows were apportioned to individual rain events using a two end-member model (i.e., mobile pre-event soil water extracted by a suction-cup sampler vs. rainwater (event water)) based on δ18O measurements. For rain events with high preceding soil moisture, mobile pre-event soil water was the main contributor (generally >60%) to the fracture flow. The colloid concentration in the surface runoff was 1-2 orders of magnitude higher than that in the subsurface flows. The lowest colloid concentration was found in the subsurface interflow, which was probably the result of pore-scale colloid straining mechanisms. The rainfall intensity and its temporal variation govern the dynamics of the colloid concentrations in both surface runoff and subsurface flows. The duration of the antecedent dry period affected not only the relative contributions of the rainwater and the mobile pre-event soil water to the subsurface flows but also the peak colloid concentration, particularly in the fracture flow. The <10 μm fine colloid size fraction accounted for more than 80% of the total suspended particles in the surface runoff, while the colloid size distributions of both the interflow and the fracture flow shifted towards larger diameters. These results highlight the need to avoid the application of strongly-sorbing agrochemicals (e.g., pesticides, phosphorus fertilizers) immediately before rainfall following a long no-rain period because their transport in association with colloids may occur rapidly over long distances via both surface runoff and subsurface flows with rainfall.

Use of sinkhole and specific capacity distributions to assess vertical gradients in a karst aquifer

USGS Publications Warehouse

McCoy, K.J.; Kozar, M.D.

2008-01-01

The carbonate-rock aquifer in the Great Valley, West Virginia, USA, was evaluated using a database of 687 sinkholes and 350 specific capacity tests to assess structural, lithologic, and topographic influences on the groundwater flow system. The enhanced permeability of the aquifer is characterized in part by the many sinkholes, springs, and solutionally enlarged fractures throughout the valley. Yet, vertical components of subsurface flow in this highly heterogeneous aquifer are currently not well understood. To address this problem, this study examines the apparent relation between geologic features of the aquifer and two spatial indices of enhanced permeability attributed to aquifer karstification: (1) the distribution of sinkholes and (2) the occurrence of wells with relatively high specific capacity. Statistical results indicate that sinkholes (funnel and collapse) occur primarily along cleavage and bedding planes parallel to subparallel to strike where lateral or downward vertical gradients are highest. Conversely, high specific capacity values are common along prominent joints perpendicular or oblique to strike. The similarity of the latter distribution to that of springs suggests these fractures are areas of upward-convergent flow. These differences between sinkhole and high specific capacity distributions suggest vertical flow components are primarily controlled by the orientation of geologic structure and associated subsurface fracturing. ?? 2007 Springer-Verlag.

Effects of spatial configuration of imperviousness and green infrastructure networks on hydrologic response in a residential sewershed

NASA Astrophysics Data System (ADS)

Lim, Theodore C.; Welty, Claire

2017-09-01

Green infrastructure (GI) is an approach to stormwater management that promotes natural processes of infiltration and evapotranspiration, reducing surface runoff to conventional stormwater drainage infrastructure. As more urban areas incorporate GI into their stormwater management plans, greater understanding is needed on the effects of spatial configuration of GI networks on hydrological performance, especially in the context of potential subsurface and lateral interactions between distributed facilities. In this research, we apply a three-dimensional, coupled surface-subsurface, land-atmosphere model, ParFlow.CLM, to a residential urban sewershed in Washington DC that was retrofitted with a network of GI installations between 2009 and 2015. The model was used to test nine additional GI and imperviousness spatial network configurations for the site and was compared with monitored pipe-flow data. Results from the simulations show that GI located in higher flow-accumulation areas of the site intercepted more surface runoff, even during wetter and multiday events. However, a comparison of the differences between scenarios and levels of variation and noise in monitored data suggests that the differences would only be detectable between the most and least optimal GI/imperviousness configurations.

Surface and subsurface cracks characteristics of single crystal SiC wafer in surface machining

NASA Astrophysics Data System (ADS)

Qiusheng, Y.; Senkai, C.; Jisheng, P.

2015-03-01

Different machining processes were used in the single crystal SiC wafer machining. SEM was used to observe the surface morphology and a cross-sectional cleavages microscopy method was used for subsurface cracks detection. Surface and subsurface cracks characteristics of single crystal SiC wafer in abrasive machining were analysed. The results show that the surface and subsurface cracks system of single crystal SiC wafer in abrasive machining including radial crack, lateral crack and the median crack. In lapping process, material removal is dominated by brittle removal. Lots of chipping pits were found on the lapping surface. With the particle size becomes smaller, the surface roughness and subsurface crack depth decreases. When the particle size was changed to 1.5µm, the surface roughness Ra was reduced to 24.0nm and the maximum subsurface crack was 1.2µm. The efficiency of grinding is higher than lapping. Plastic removal can be achieved by changing the process parameters. Material removal was mostly in brittle fracture when grinding with 325# diamond wheel. Plow scratches and chipping pits were found on the ground surface. The surface roughness Ra was 17.7nm and maximum subsurface crack depth was 5.8 µm. When grinding with 8000# diamond wheel, the material removal was in plastic flow. Plastic scratches were found on the surface. A smooth surface of roughness Ra 2.5nm without any subsurface cracks was obtained. Atomic scale removal was possible in cluster magnetorheological finishing with diamond abrasive size of 0.5 µm. A super smooth surface eventually obtained with a roughness of Ra 0.4nm without any subsurface crack.

Discrepancy between Snowmelt and Soil Infiltration

NASA Astrophysics Data System (ADS)

Fassnacht, S. R.

2017-12-01

A majority of snowmelt enters the soil and is either transmitted through or stored in the soil. Snowmelt has been estimated from the decrease in snow mass of a snow pillow and soil infiltration has been estimated from near surface TDR probes. Here, these data are from a set of Snow Telemetry (SNOTEL) stations across Colorado. While seasonal totals are similar, it is shown that there is a disconnect between the amount of water melted in a day and the increased daily volume of water measured in the near sub-surface. It is surmised that these differences are a function of the data collection methods, the infiltration rate, and possible lateral flow. An examination of daily infiltration volumes at depth shows a further disconnect, as it is likely that lateral flow complicates the measurements to a true three dimensional problem. The data are informative to illustrate the transmission of meltwater into the soil; methods for improvement are explored.

Preferential flow, diffuse flow, and perching in an interbedded fractured-rock unsaturated zone

NASA Astrophysics Data System (ADS)

Nimmo, John R.; Creasey, Kaitlyn M.; Perkins, Kim S.; Mirus, Benjamin B.

2017-03-01

Layers of strong geologic contrast within the unsaturated zone can control recharge and contaminant transport to underlying aquifers. Slow diffuse flow in certain geologic layers, and rapid preferential flow in others, complicates the prediction of vertical and lateral fluxes. A simple model is presented, designed to use limited geological site information to predict these critical subsurface processes in response to a sustained infiltration source. The model is developed and tested using site-specific information from the Idaho National Laboratory in the Eastern Snake River Plain (ESRP), USA, where there are natural and anthropogenic sources of high-volume infiltration from floods, spills, leaks, wastewater disposal, retention ponds, and hydrologic field experiments. The thick unsaturated zone overlying the ESRP aquifer is a good example of a sharply stratified unsaturated zone. Sedimentary interbeds are interspersed between massive and fractured basalt units. The combination of surficial sediments, basalts, and interbeds determines the water fluxes through the variably saturated subsurface. Interbeds are generally less conductive, sometimes causing perched water to collect above them. The model successfully predicts the volume and extent of perching and approximates vertical travel times during events that generate high fluxes from the land surface. These developments are applicable to sites having a thick, geologically complex unsaturated zone of substantial thickness in which preferential and diffuse flow, and perching of percolated water, are important to contaminant transport or aquifer recharge.

Preferential flow, diffuse flow, and perching in an interbedded fractured-rock unsaturated zone

USGS Publications Warehouse

Nimmo, John R.; Creasey, Kaitlyn M; Perkins, Kimberlie; Mirus, Benjamin B.

2017-01-01

Layers of strong geologic contrast within the unsaturated zone can control recharge and contaminant transport to underlying aquifers. Slow diffuse flow in certain geologic layers, and rapid preferential flow in others, complicates the prediction of vertical and lateral fluxes. A simple model is presented, designed to use limited geological site information to predict these critical subsurface processes in response to a sustained infiltration source. The model is developed and tested using site-specific information from the Idaho National Laboratory in the Eastern Snake River Plain (ESRP), USA, where there are natural and anthropogenic sources of high-volume infiltration from floods, spills, leaks, wastewater disposal, retention ponds, and hydrologic field experiments. The thick unsaturated zone overlying the ESRP aquifer is a good example of a sharply stratified unsaturated zone. Sedimentary interbeds are interspersed between massive and fractured basalt units. The combination of surficial sediments, basalts, and interbeds determines the water fluxes through the variably saturated subsurface. Interbeds are generally less conductive, sometimes causing perched water to collect above them. The model successfully predicts the volume and extent of perching and approximates vertical travel times during events that generate high fluxes from the land surface. These developments are applicable to sites having a thick, geologically complex unsaturated zone of substantial thickness in which preferential and diffuse flow, and perching of percolated water, are important to contaminant transport or aquifer recharge.

Cometary activity and nucleus modelling: a new approach

NASA Astrophysics Data System (ADS)

Möhlmann, D.

1996-06-01

The phenomena of comet splittings with an average frequency of about one splitting per 100 years and comet (Chen and Jewitt, Icarus108, 265-271, 1994), and the restriction of cometary activity to well-defined small areas at the almost passive and mantle covered surface (Keller et al., ESA SP-250, Vol. II, pp. 363-364, 1986) are at present driving challenges to models of structure and evolution of comet nuclei. Extending the presently discussed models by incorporating lateral subsurface transport of sublimed volatiles, there appears the possibility that the places of sublimation are different from those of activity (the so-called active areas). Then, there is no necessity to distinguish between different surface properties at active and passive areas, assuming, e.g. an uncovered icy surface at active areas. Active areas are simply the very local "source sites" where the accumulated subsurface flows from distant regions reach the surface. The pressure driven subsurface flows of volatiles may not only leave the comet at its surface, they may penetrate via cracks, etc. also deeply into the nucleus. There they can cause a further growth of cracks and also new cracks. This can be a cause for the observed regular splittings. Furthermore, actual models (Kührt and Keller, Icarus109, 121-132, 1994; Skorov and Rickman, Planet. Space Sci.43, 1587-1594, 1995) of the gas transport through porous comet surface crusts can be interpreted as to give first indications for thermodynamical parameters in heat conducting and porous cometary crusts which are appropriate for 1 AU conditions to permit the temporary existence of a layer with fluid subsurface water within these crusts. This exciting result of the possible temporary existence of subsurface warm water in comets which approach the Sun within about 1 AU makes a cometary subsurface chemistry much more efficient than expected hitherto.

Fractal topography and subsurface water flows from fluvial bedforms to the continental shield

USGS Publications Warehouse

Worman, A.; Packman, A.I.; Marklund, L.; Harvey, J.W.; Stone, S.H.

2007-01-01

Surface-subsurface flow interactions are critical to a wide range of geochemical and ecological processes and to the fate of contaminants in freshwater environments. Fractal scaling relationships have been found in distributions of both land surface topography and solute efflux from watersheds, but the linkage between those observations has not been realized. We show that the fractal nature of the land surface in fluvial and glacial systems produces fractal distributions of recharge, discharge, and associated subsurface flow patterns. Interfacial flux tends to be dominated by small-scale features while the flux through deeper subsurface flow paths tends to be controlled by larger-scale features. This scaling behavior holds at all scales, from small fluvial bedforms (tens of centimeters) to the continental landscape (hundreds of kilometers). The fractal nature of surface-subsurface water fluxes yields a single scale-independent distribution of subsurface water residence times for both near-surface fluvial systems and deeper hydrogeological flows. Copyright 2007 by the American Geophysical Union.

Virtual experiments: a new approach for improving process conceptualization in hillslope hydrology

NASA Astrophysics Data System (ADS)

Weiler, Markus; McDonnell, Jeff

2004-01-01

We present an approach for process conceptualization in hillslope hydrology. We develop and implement a series of virtual experiments, whereby the interaction between water flow pathways, source and mixing at the hillslope scale is examined within a virtual experiment framework. We define these virtual experiments as 'numerical experiments with a model driven by collective field intelligence'. The virtual experiments explore the first-order controls in hillslope hydrology, where the experimentalist and modeler work together to cooperatively develop and analyze the results. Our hillslope model for the virtual experiments (HillVi) in this paper is based on conceptualizing the water balance within the saturated and unsaturated zone in relation to soil physical properties in a spatially explicit manner at the hillslope scale. We argue that a virtual experiment model needs to be able to capture all major controls on subsurface flow processes that the experimentalist might deem important, while at the same time being simple with few 'tunable parameters'. This combination makes the approach, and the dialog between experimentalist and modeler, a useful hypothesis testing tool. HillVi simulates mass flux for different initial conditions under the same flow conditions. We analyze our results in terms of an artificial line source and isotopic hydrograph separation of water and subsurface flow. Our results for this first set of virtual experiments showed how drainable porosity and soil depth variability exert a first order control on flow and transport at the hillslope scale. We found that high drainable porosity soils resulted in a restricted water table rise, resulting in more pronounced channeling of lateral subsurface flow along the soil-bedrock interface. This in turn resulted in a more anastomosing network of tracer movement across the slope. The virtual isotope hydrograph separation showed higher proportions of event water with increasing drainable porosity. When combined with previous experimental findings and conceptualizations, virtual experiments can be an effective way to isolate certain controls and examine their influence over a range of rainfall and antecedent wetness conditions.

Dynamic coupling of subsurface and seepage flows solved within a regularized partition formulation

NASA Astrophysics Data System (ADS)

Marçais, J.; de Dreuzy, J.-R.; Erhel, J.

2017-11-01

Hillslope response to precipitations is characterized by sharp transitions from purely subsurface flow dynamics to simultaneous surface and subsurface flows. Locally, the transition between these two regimes is triggered by soil saturation. Here we develop an integrative approach to simultaneously solve the subsurface flow, locate the potential fully saturated areas and deduce the generated saturation excess overland flow. This approach combines the different dynamics and transitions in a single partition formulation using discontinuous functions. We propose to regularize the system of partial differential equations and to use classic spatial and temporal discretization schemes. We illustrate our methodology on the 1D hillslope storage Boussinesq equations (Troch et al., 2003). We first validate the numerical scheme on previous numerical experiments without saturation excess overland flow. Then we apply our model to a test case with dynamic transitions from purely subsurface flow dynamics to simultaneous surface and subsurface flows. Our results show that discretization respects mass balance both locally and globally, converges when the mesh or time step are refined. Moreover the regularization parameter can be taken small enough to ensure accuracy without suffering of numerical artefacts. Applied to some hundreds of realistic hillslope cases taken from Western side of France (Brittany), the developed method appears to be robust and efficient.

Mechanistic assessment of hillslope transpiration controls of diel subsurface flow: a steady-state irrigation approach

Treesearch

H.R. Barnard; C.B. Graham; W.J. van Verseveld; J.R. Brooks; B.J. Bond; J.J. McDonnell

2010-01-01

Mechanistic assessment of how transpiration influences subsurface flow is necessary to advance understanding of catchment hydrology. We conducted a 24-day, steady-state irrigation experiment to quantify the relationships among soil moisture, transpiration and hillslope subsurface flow. Our objectives were to: (1) examine the time lag between maximum transpiration and...

A porewater - based stable isotope approach for the investigation of subsurface hydrological processes

NASA Astrophysics Data System (ADS)

Garvelmann, J.; Külls, C.; Weiler, M.

2011-10-01

Predicting and understanding subsurface flowpaths is still a crucial issue in hydrological research. We present an experimental approach to reveal present and past subsurface flowpaths of water in the unsaturated and saturated zone. Two hillslopes in a humid moutainous catchment have been investigated. The H2O(liquid) - H2O(vapor) equilibration laser spectroscopy method was used to obtain high resolution δ2H vertical depth profiles of porewater at various points along a fall line of a pasture hillslope in the southern Black Forest, Germany. The Porewater Stable Isotope Profile (PSIP) approach was developed to use the integrated information of several vertical depth profiles of deuterium along two transects at the hillslopes. Different shapes of depth profiles were observed in relation to hillslope position. The statistical variability (inter-quartile range and standard deviation) of each profile was used to characterize different types of depth profiles. The profiles upslope or with a weak affinity for saturation as indicated by a low topographic wetness index preserve the isotopic input signal by precipitation with a distinct seasonal variability. These observations indicate mainly vertical movement of soil water in the upper part of the hillslope before sampling. The profiles downslope or at locations with a strong affinity for saturation do not show a similar seasonal isotopic signal. The input signal is erased in the foothills and a large proportion of pore water samples are close to the isotopic values of δ2H in stream water during base flow. Near the stream indications for efficient mixing of water from lateral subsurface flow paths with vertical percolation are found.

Modeling Subsurface Storm and Tile Drain Systems in GSSHA with SUPERLINK

DTIC Science & Technology

2014-09-01

side is computed as . ( )e Kq d m m L   2 0 01 2 (7) de is defined as ( van der Molen and Wesseling (1991)) ERDC/CHL TR-14-11 15...Conservation Service. Van der Molen , W.H., and J. Wesseling. 1991. A solution in closed form and a series solution to replace the tables for thickness of...effective lateral hydraulic conductivity (cm hr-1) C = 1 in the present version. Hooghoudt ( van Schilfgaarde 1974) characterized flow to cylindrical

Impacts of microtopographic snow redistribution and lateral subsurface processes on hydrologic and thermal states in an Arctic polygonal ground ecosystem: a case study using ELM-3D v1.0

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

Bisht, Gautam; Riley, William J.; Wainwright, Haruko M.

Microtopographic features, such as polygonal ground, are characteristic sources of landscape heterogeneity in the Alaskan Arctic coastal plain. In this study, we analyze the effects of snow redistribution (SR) and lateral subsurface processes on hydrologic and thermal states at a polygonal tundra site near Barrow, Alaska. We extended the land model integrated in the E3SM to redistribute incoming snow by accounting for microtopography and incorporated subsurface lateral transport of water and energy (ELM-3D v1.0). Multiple 10-year-long simulations were performed for a transect across a polygonal tundra landscape at the Barrow Environmental Observatory in Alaska to isolate the impact of SRmore » and subsurface process representation. When SR was included, model predictions better agreed (higher R 2, lower bias and RMSE) with observed differences in snow depth between polygonal rims and centers. The model was also able to accurately reproduce observed soil temperature vertical profiles in the polygon rims and centers (overall bias, RMSE, and R 2 of 0.59°C, 1.82°C, and 0.99, respectively). The spatial heterogeneity of snow depth during the winter due to SR generated surface soil temperature heterogeneity that propagated in depth and time and led to ~ 10 cm shallower and ~ 5 cm deeper maximum annual thaw depths under the polygon rims and centers, respectively. Additionally, SR led to spatial heterogeneity in surface energy fluxes and soil moisture during the summer. Excluding lateral subsurface hydrologic and thermal processes led to small effects on mean states but an overestimation of spatial variability in soil moisture and soil temperature as subsurface liquid pressure and thermal gradients were artificially prevented from spatially dissipating over time. The effect of lateral subsurface processes on maximum thaw depths was modest, with mean absolute differences of ~ 3 cm. Our integration of three-dimensional subsurface hydrologic and thermal subsurface dynamics in the E3SM land model will facilitate a wide range of analyses heretofore impossible in an ESM context.« less

Impacts of microtopographic snow-redistribution and lateral subsurface processeson hydrologic and thermal states in an Arctic polygonal ground ecosystem

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

Bisht, Gautam; Riley, William J.; Wainwright, Haruko M.

Microtopographic features, such as polygonal ground, are characteristic sources of landscape heterogeneity in the Alaskan Arctic coastal plain. We analyze the effects of snow redistribution (SR) and lateral subsurface processes on hydrologic and thermal states at a polygonal tundra site near Barrow, Alaska. We extended the land model integrated in the ACME Earth System Model (ESM) to redistribute incoming snow by accounting for microtopography and incorporated subsurface lateral transport of water and energy (ALMv0-3D). Three 10-years long simulations were performed for a transect across polygonal tundra landscape at the Barrow Environmental Observatory in Alaska to isolate the impact of SRmore » and subsurface process representation. When SR was included, model results show a better agreement (higher R 2 with lower bias and RMSE) for the observed differences in snow depth between polygonal rims and centers. The model was also able to accurately reproduce observed soil temperature vertical profiles in the polygon rims and centers (overall bias, RMSE, and R 2 of 0.59°C, 1.82°C, and 0.99, respectively). The spatial heterogeneity of snow depth during the winter due to SR generated surface soil temperature heterogeneity that propagated in depth and time and led to ~10 cm shallower and ~5 cm deeper maximum annual thaw depths under the polygon rims and centers, respectively. Additionally, SR led to spatial heterogeneity in surface energy fluxes and soil moisture during the summer. Excluding lateral subsurface hydrologic and thermal processes led to small effects on mean states but an overestimation of spatial variability in soil moisture and soil temperature as subsurface liquid pressure and thermal gradients were artificially prevented from spatially dissipating over time. The effect of lateral subsurface processes on active layer depths was modest with mean absolute difference of ~3 cm. Finally, our integration of three-dimensional subsurface hydrologic and thermal subsurface dynamics in the ACME land model will facilitate a wide range of analyses heretofore impossible in an ESM context.« less

Impacts of microtopographic snow redistribution and lateral subsurface processes on hydrologic and thermal states in an Arctic polygonal ground ecosystem: a case study using ELM-3D v1.0

DOE PAGES

Bisht, Gautam; Riley, William J.; Wainwright, Haruko M.; ...

2018-01-08

Microtopographic features, such as polygonal ground, are characteristic sources of landscape heterogeneity in the Alaskan Arctic coastal plain. In this study, we analyze the effects of snow redistribution (SR) and lateral subsurface processes on hydrologic and thermal states at a polygonal tundra site near Barrow, Alaska. We extended the land model integrated in the E3SM to redistribute incoming snow by accounting for microtopography and incorporated subsurface lateral transport of water and energy (ELM-3D v1.0). Multiple 10-year-long simulations were performed for a transect across a polygonal tundra landscape at the Barrow Environmental Observatory in Alaska to isolate the impact of SRmore » and subsurface process representation. When SR was included, model predictions better agreed (higher R 2, lower bias and RMSE) with observed differences in snow depth between polygonal rims and centers. The model was also able to accurately reproduce observed soil temperature vertical profiles in the polygon rims and centers (overall bias, RMSE, and R 2 of 0.59°C, 1.82°C, and 0.99, respectively). The spatial heterogeneity of snow depth during the winter due to SR generated surface soil temperature heterogeneity that propagated in depth and time and led to ~ 10 cm shallower and ~ 5 cm deeper maximum annual thaw depths under the polygon rims and centers, respectively. Additionally, SR led to spatial heterogeneity in surface energy fluxes and soil moisture during the summer. Excluding lateral subsurface hydrologic and thermal processes led to small effects on mean states but an overestimation of spatial variability in soil moisture and soil temperature as subsurface liquid pressure and thermal gradients were artificially prevented from spatially dissipating over time. The effect of lateral subsurface processes on maximum thaw depths was modest, with mean absolute differences of ~ 3 cm. Our integration of three-dimensional subsurface hydrologic and thermal subsurface dynamics in the E3SM land model will facilitate a wide range of analyses heretofore impossible in an ESM context.« less

Impacts of microtopographic snow-redistribution and lateral subsurface processeson hydrologic and thermal states in an Arctic polygonal ground ecosystem

DOE PAGES

Bisht, Gautam; Riley, William J.; Wainwright, Haruko M.; ...

2018-01-08

Microtopographic features, such as polygonal ground, are characteristic sources of landscape heterogeneity in the Alaskan Arctic coastal plain. We analyze the effects of snow redistribution (SR) and lateral subsurface processes on hydrologic and thermal states at a polygonal tundra site near Barrow, Alaska. We extended the land model integrated in the ACME Earth System Model (ESM) to redistribute incoming snow by accounting for microtopography and incorporated subsurface lateral transport of water and energy (ALMv0-3D). Three 10-years long simulations were performed for a transect across polygonal tundra landscape at the Barrow Environmental Observatory in Alaska to isolate the impact of SRmore » and subsurface process representation. When SR was included, model results show a better agreement (higher R 2 with lower bias and RMSE) for the observed differences in snow depth between polygonal rims and centers. The model was also able to accurately reproduce observed soil temperature vertical profiles in the polygon rims and centers (overall bias, RMSE, and R 2 of 0.59°C, 1.82°C, and 0.99, respectively). The spatial heterogeneity of snow depth during the winter due to SR generated surface soil temperature heterogeneity that propagated in depth and time and led to ~10 cm shallower and ~5 cm deeper maximum annual thaw depths under the polygon rims and centers, respectively. Additionally, SR led to spatial heterogeneity in surface energy fluxes and soil moisture during the summer. Excluding lateral subsurface hydrologic and thermal processes led to small effects on mean states but an overestimation of spatial variability in soil moisture and soil temperature as subsurface liquid pressure and thermal gradients were artificially prevented from spatially dissipating over time. The effect of lateral subsurface processes on active layer depths was modest with mean absolute difference of ~3 cm. Finally, our integration of three-dimensional subsurface hydrologic and thermal subsurface dynamics in the ACME land model will facilitate a wide range of analyses heretofore impossible in an ESM context.« less

Impacts of microtopographic snow redistribution and lateral subsurface processes on hydrologic and thermal states in an Arctic polygonal ground ecosystem: a case study using ELM-3D v1.0

NASA Astrophysics Data System (ADS)

Bisht, Gautam; Riley, William J.; Wainwright, Haruko M.; Dafflon, Baptiste; Yuan, Fengming; Romanovsky, Vladimir E.

2018-01-01

Microtopographic features, such as polygonal ground, are characteristic sources of landscape heterogeneity in the Alaskan Arctic coastal plain. Here, we analyze the effects of snow redistribution (SR) and lateral subsurface processes on hydrologic and thermal states at a polygonal tundra site near Barrow, Alaska. We extended the land model integrated in the E3SM to redistribute incoming snow by accounting for microtopography and incorporated subsurface lateral transport of water and energy (ELM-3D v1.0). Multiple 10-year-long simulations were performed for a transect across a polygonal tundra landscape at the Barrow Environmental Observatory in Alaska to isolate the impact of SR and subsurface process representation. When SR was included, model predictions better agreed (higher R2, lower bias and RMSE) with observed differences in snow depth between polygonal rims and centers. The model was also able to accurately reproduce observed soil temperature vertical profiles in the polygon rims and centers (overall bias, RMSE, and R2 of 0.59 °C, 1.82 °C, and 0.99, respectively). The spatial heterogeneity of snow depth during the winter due to SR generated surface soil temperature heterogeneity that propagated in depth and time and led to ˜ 10 cm shallower and ˜ 5 cm deeper maximum annual thaw depths under the polygon rims and centers, respectively. Additionally, SR led to spatial heterogeneity in surface energy fluxes and soil moisture during the summer. Excluding lateral subsurface hydrologic and thermal processes led to small effects on mean states but an overestimation of spatial variability in soil moisture and soil temperature as subsurface liquid pressure and thermal gradients were artificially prevented from spatially dissipating over time. The effect of lateral subsurface processes on maximum thaw depths was modest, with mean absolute differences of ˜ 3 cm. Our integration of three-dimensional subsurface hydrologic and thermal subsurface dynamics in the E3SM land model will facilitate a wide range of analyses heretofore impossible in an ESM context.

Seasonal GPR Signal Changes in Two Contrasting Soils in the Shale Hills Catchment

NASA Astrophysics Data System (ADS)

Lin, H.; Zhang, J.; Doolittle, J. A.

2011-12-01

Repeated GPR surveys in different seasons, combined with real-time soil water monitoring, provide a useful methodology to reveal subsurface hydrologic processes and their underlying mechanisms in different soils and hillslopes. This was demonstrated in the Shale Hills Critical Zone Observatory using two contrasting soils over several dry and wet seasons. Our results showed that 1) the radar reflection in the BC-C horizon interface in the deep Rushtown soil became clearer as soil became wetter, which was linked to lateral flow above this horizon interface that increased the contrast, and 2) the reflection in the soil-bedrock interface and the weathered-unweathered rock interface in the shallow Weikert soil become intermittent as soil became wetter, which was attributed to non-uniform distribution of water in bedrock fractures that created locally strong contrast, leading to point scatter of GPR reflection. This study shows the optimal time for using GPR to detect soil horizon interfaces, the value of nondestructive mapping of soil-rock moisture distribution patterns, and the possibility of identifying preferential flow pathways in the subsurface.

Vertical Subsurface Flow Mixing and Horizontal Anisotropy in Coarse Fluvial Aquifers: Structural Aspects

NASA Astrophysics Data System (ADS)

Huggenberger, P.; Huber, E.

2014-12-01

Detailed descriptions of the subsurface heterogeneities in coarse fluvial aquifer gravel often lack in concepts to distinguish between the essence and the noise of a permeability structure and the ability to extrapolate site specific hydraulic information at the tens to several hundred meters scale. At this scale the heterogeneity strongly influences the anisotropies of the flow field and the mixing processes in groundwater. However, in many hydrogeological models the complexity of natural systems is oversimplified. Understanding the link between the dynamics of the surface processes of braided-river systems and the resulting subsurface sedimentary structures is the key to characterizing the complexity of horizontal and vertical mixing processes in groundwater. From the different depositional elements of coarse braided-river systems, the largest permeability contrasts can be observed in the scour-fills. Other elements (e.g. different types of gravel sheets) show much smaller variabilities and could be considered as a kind of matrix. Field experiments on the river Tagliamento (Northeast Italy) based on morphological observation and ground-penetrating radar (GPR) surveys, as well as outcrop analyses of gravel pit exposures (Switzerland) allowed us to define the shape, sizes, spatial distribution and preservation potential of scour-fills. In vertical sections (e.g. 2D GPR data, vertical outcrop), the spatial density of remnant erosional bounding surfaces of scours is an indicator for the dynamics of the braided-river system (lateral mobility of the active floodplain, rate of sediment net deposition and spatial distribution of the confluence scours). In case of combined low aggradation rate and low lateral mobility the deposits may be dominated by a complex overprinting of scour-fills. The delineation of the erosional bounding surfaces, that are coherent over the survey area, is based on the identification of angular discontinuities of the reflectors. Fence diagrams and horizontal time-slices from GPR data are used to construct simplified 3D hydraulic properties distribution models and to derive anisotropy patterns. On the basis of this work, conceptual models could be designed and implemented into numerical models to simulate the flow field and mixing in heterogeneous braided-river deposits.

Geophysical Characterization of Groundwater-Fault Dynamics at San Andreas Oasis

NASA Astrophysics Data System (ADS)

Faherty, D.; Polet, J.; Osborn, S. G.

2017-12-01

The San Andreas Oasis has historically provided a reliable source of fresh water near the northeast margin of the Salton Sea, although since the recent completion of the Coachella Canal Lining Project and persistent drought in California, surface water at the site has begun to disappear. This may be an effect of the canal lining, however, the controls on groundwater are complicated by the presence of the Hidden Springs Fault (HSF), a northeast dipping normal fault that trends near the San Andreas Oasis. Its surface expression is apparent as a lineation against which all plant growth terminates, suggesting that it may form a partial barrier to subsurface groundwater flow. Numerous environmental studies have detailed the chemical evolution of waters resources at San Andreas Spring, although there remains a knowledge gap on the HSF and its relation to groundwater at the site. To better constrain flow paths and characterize groundwater-fault interactions, we have employed resistivity surveys near the surface trace of the HSF to generate profiles of lateral and depth-dependent variations in resistivity. The survey design is comprised of lines installed in Wenner Arrays, using an IRIS Syscal Kid, with 24 electrodes, at a maximum electrode spacing of 5 meters. In addition, we have gathered constraints on the geometry of the HSF using a combination of ground-based magnetic and gravity profiles, conducted with a GEM walking Proton Precession magnetometer and a Lacoste & Romberg gravimeter. Seventeen gravity measurements were acquired across the surface trace of the fault. Preliminary resistivity results depict a shallow conductor localized at the oasis and discontinuous across the HSF. Magnetic data reveal a large contrast in subsurface magnetic susceptibility that appears coincident with the surface trace and trend of the HSF, while gravity data suggests a shallow, relatively high density anomaly centered near the oasis. These data also hint at a second, previously undocumented fault bounding the opposite margin of the oasis and trending subparallel to the HSF. We thus speculate that the Hidden Springs Fault and this possible secondary fault act as partial barriers to lateral subsurface flow and form a structural wedge, localizing groundwater beneath the oasis.

Role of Climatic Variability on Fate and Transport of LNAPL Pollutants in Subsurface

NASA Astrophysics Data System (ADS)

Gupta, P. K.; Yadav, B. K.

2017-12-01

Climatic variability affects groundwater resources both directly through replenishment by surface recharge and indirectly via changes in groundwater extraction patterns. Remediation of polluted groundwater due to the release of mono-aromatic hydrocarbons such as light non-aqueous phase liquids (LNAPL) is of particular concern under changing climatic conditions because of their higher water solubility and wide coverage in the subsurface. Thus, the aim of this study was to investigate the impact of these shallow groundwater extremes on biodegradation and transport of toluene, the selected LNAPL, in subsurface using a series of practical and numerical experiments. An air tight three-dimensional sand tank setup (60cm-L×30cm-W×60cm-D) embedded with horizontal and vertical layers of sampling ports was fabricated using a glass sheet of thickness 7 mm. Clean sand having an average grain size of 0.5-1.0 mm was packed homogeneously for creating an unconfined aquifer. Pure phase of toluene was released from the top surface to create a pool of the LNAPL around the groundwater table which was maintained at 35 cm from the tank bottom. Initially, a constant water flux was allowed to flow to maintain a pore water velocity of 1.2 m/day in lateral direction to mimic a base groundwater flow regime. Subsequently, faster and slow groundwater velocity regimes were developed by changing the water flux through the saturated zone keeping the water table location at the same level. The observed breakthrough curves at different ports showed that lateral and transverse transport of the LNAPL was more prominent as compared to its vertical movement. The increased vertical movement of the LNAPL along with an enhanced dissolution rate under the faster groundwater flow condition shows the crucial role of mechanical dispersion and the shear force acting on the water-LNAPL interface. The rate of degradation was found high for the case of faster pore water velocities due to dependency of the degradation kinetics on substrate concentration. The observed data were compared well with the simulated curves for all the three cases of groundwater flow conditions. The results of this study are of direct use in applying bioremediation technique in the field and for planning of LNAPL polluting industrial locations under changing climatic conditions.

Shallow subsurface storm flow in a forested headwater catchment: Observations and modeling using a modified TOPMODEL

USGS Publications Warehouse

Scanlon, Todd M.; Raffensperger, Jeff P.; Hornberger, George M.; Clapp, Roger B.

2000-01-01

Transient, perched water tables in the shallow subsurface are observed at the South Fork Brokenback Run catchment in Shenandoah National Park, Virginia. Crest piezometers installed along a hillslope transect show that the development of saturated conditions in the upper 1.5 m of the subsurface is controlled by total precipitation and antecedent conditions, not precipitation intensity, although soil heterogeneities strongly influence local response. The macroporous subsurface storm flow zone provides a hydrological pathway for rapid runoff generation apart from the underlying groundwater zone, a conceptualization supported by the two‐storage system exhibited by hydrograph recession analysis. A modified version of TOPMODEL is used to simulate the observed catchment dynamics. In this model, generalized topographic index theory is applied to the subsurface storm flow zone to account for logarithmic storm flow recessions, indicative of linearly decreasing transmissivity with depth. Vertical drainage to the groundwater zone is required, and both subsurface reservoirs are considered to contribute to surface saturation.

Anatomy Of The ‘LuSi’ Mud Eruption, East Java

NASA Astrophysics Data System (ADS)

Tingay, M. R.

2009-12-01

Early in the morning of the 29th of May 2006, hot mud started erupting from the ground in the densely populated Porong District of Sidoarjo, East Java. With initial flow rates of ~5000 cubic meters per day, the mud quickly inundated neighbouring villages. Over two years later and the ‘Lusi’ eruption has increased in strength, expelling over 90 million cubic meters of mud at an average rate of approximately 100000 cubic meters per day. The mud flow has now covered over 700 hectares of land to depths of over 25 meters, engulfing a dozen villages and displacing approximately 40000 people. In addition to the inundated areas, other areas are also at risk from subsidence and distant eruptions of gas. However, efforts to stem the mud flow or monitor its evolution are hampered by our overall lack of knowledge and consensus on the subsurface anatomy of the Lusi mud volcanic system. In particular, the largest and most significant uncertainties are the source of the erupted water (shales versus deep carbonates), the fluid flow pathways (purely fractures versus mixed fracture and wellbore) and disputes over the subsurface geology (nature of deep carbonates, lithology of rocks between shale and carbonates). This study will present and overview of the anatomy of the Lusi mud volcanic system with particular emphasis on these critical uncertainties and their influence on the likely evolution of disaster.

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.

Enhanced Remedial Amendment Delivery to Subsurface Using Shear Thinning Fluid and Aqueous Foam

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

Zhong, Lirong; Szecsody, James E.; Oostrom, Martinus

2011-04-23

A major issue with in situ subsurface remediation is the ability to achieve an even spatial distribution of remedial amendments to the contamination zones in an aquifer or vadose zone. Delivery of amendment to the aquifer using shear thinning fluid and to the vadose zone using aqueous foam has the potential to enhance the amendment distribution into desired locations and improve the remediation. 2-D saturated flow cell experiments were conducted to evaluate the enhanced sweeping, contaminant removal, and amendment persistence achieved by shear thinning fluid delivery. Bio-polymer xanthan gum solution was used as the shear thinning fluid. Unsaturated 1-D columnmore » and 2-D flow cell experiments were conducted to evaluate the mitigation of contaminant mobilization, amendment uniform distribution enhancement, and lateral delivery improvement by foam delivery. Surfactant sodium lauryl ether sulfate was used as the foaming agent. It was demonstrated that the shear thinning fluid injection enhanced the fluid sweeping over a heterogeneous system and increased the delivery of remedial amendment into low-permeability zones. The persistence of the amendment distributed into the low-perm zones by the shear thinning fluid was prolonged compared to that of amendment distributed by water injection. Foam delivery of amendment was shown to mitigate the mobilization of highly mobile contaminant from sediments under vadose zone conditions. Foam delivery also achieved more uniform amendment distribution in a heterogeneous unsaturated system, and demonstrated remarkable increasing in lateral distribution of the injected liquid compared to direct liquid injection.« less

Performance of a pilot showcase of different wetland systems in an urban setting in Singapore.

PubMed

Quek, B S; He, Q H; Sim, C H

2015-01-01

The Alexandra Wetlands, part of PUB's Active, Beautiful, Clean Waters (ABC Waters) Programme, showcase a surface flow wetland, an aquatic pond and a sub-surface flow wetland on a 200 m deck built over an urban drainage canal. Water from the canal is pumped to a sedimentation basin, before flowing in parallel to the three wetlands. Water quality monitoring was carried out monthly from April 2011 to December 2012. The order of removal efficiency is sub-surface flow (81.3%) >aquatic pond (58.5%) >surface flow (50.7%) for total suspended solids (TSS); sub-surface (44.9%) >surface flow (31.9%) >aquatic pond (22.0%) for total nitrogen (TN); and surface flow (56.7%) >aquatic pond (39.8%) >sub-surface flow (5.4%) for total phosphorus (TP). All three wetlands achieved the Singapore stormwater treatment objectives (STO) for TP removal, but only the sub-surface flow wetland met the STO for TSS, and none met the STO for TN. Challenges in achieving satisfactory performance include inconsistent feed water quality, undesirable behaviour such as fishing, release of pets and feeding of animals in the wetlands, and canal dredging during part of the monitoring period. As a pilot showcase, the Alexandra Wetlands provide useful lessons for implementing multi-objective wetlands in an urban setting.

Discharge-nitrate data clustering for characterizing surface-subsurface flow interaction and calibration of a hydrologic model

NASA Astrophysics Data System (ADS)

Shrestha, R. R.; Rode, M.

2008-12-01

Concentration of reactive chemicals has different chemical signatures in baseflow and surface runoff. Previous studies on nitrate export from a catchment indicate that the transport processes are driven by subsurface flow. Therefore nitrate signature can be used for understanding the event and pre-event contributions to streamflow and surface-subsurface flow interactions. The study uses flow and nitrate concentration time series data for understanding the relationship between these two variables. Unsupervised artificial neural network based learning method called self organizing map is used for the identification of clusters in the datasets. Based on the cluster results, five different pattern in the datasets are identified which correspond to (i) baseflow, (ii) subsurface flow increase, (iii) surface runoff increase, (iv) surface runoff recession, and (v) subsurface flow decrease regions. The cluster results in combination with a hydrologic model are used for discharge separation. For this purpose, a multi-objective optimization tool NSGA-II is used, where violation of cluster results is used as one of the objective functions. The results show that the use of cluster results as supplementary information for the calibration of a hydrologic model gives a plausible simulation of subsurface flow as well total runoff at the catchment outlet. The study is undertaken using data from the Weida catchment in the North-Eastern Germany, which is a sub-catchment of the Weisse Elster river in the Elbe river basin.

Exact three-dimensional spectral solution to surface-groundwater interactions with arbitrary surface topography

USGS Publications Warehouse

Worman, A.; Packman, A.I.; Marklund, L.; Harvey, J.W.; Stone, S.H.

2006-01-01

It has been long known that land surface topography governs both groundwater flow patterns at the regional-to-continental scale and on smaller scales such as in the hyporheic zone of streams. Here we show that the surface topography can be separated in a Fourier-series spectrum that provides an exact solution of the underlying three-dimensional groundwater flows. The new spectral solution offers a practical tool for fast calculation of subsurface flows in different hydrological applications and provides a theoretical platform for advancing conceptual understanding of the effect of landscape topography on subsurface flows. We also show how the spectrum of surface topography influences the residence time distribution for subsurface flows. The study indicates that the subsurface head variation decays exponentially with depth faster than it would with equivalent two-dimensional features, resulting in a shallower flow interaction. Copyright 2006 by the American Geophysical Union.

Export of nutrients and major ionic solutes from a rain forest catchment in the Central Amazon Basin

NASA Astrophysics Data System (ADS)

Lesack, Lance F. W.

1993-03-01

The relative roles of base flow runoff versus storm flow runoff versus subsurface outflow in controlling total export of solutes from a 23.4-ha catchment of undisturbed rain forest in the central Amazon Basin were evaluated from water and solute flux measurements performed over a 1 year period. Solutes exported via 173 storms during the study were estimated from stream water samples collected during base flow conditions and during eight storms, and by utilizing a hydrograph separation technique in combination with a mixing model to partition storm flow from base flow fluxes. Solutes exported by subsurface outflow were estimated from groundwater samples from three nests of piezometers installed into the streambed, and concurrent measurements of hydraulic conductivity and hydraulic head gradients. Base flow discharge represented 92% of water outflow from the basin and was the dominant pathway of solute export. Although storm flow discharge represented only 5% of total water outflow, storm flow solute fluxes represented up to 25% of the total annual export flux, though for many solutes the portion was less. Subsurface outflow represented only 2.5% of total water outflow, and subsurface solute fluxes never represented more than 5% of the total annual export flux. Measurement errors were relatively high for storm flow and subsurface outflow fluxes, but cumulative measurement errors associated with the total solute fluxes exported from the catchment, in most cases, ranged from only ±7% to 14% because base flow fluxes were measured relatively well. The export fluxes of most solutes are substantially less than previously reported for comparable small catchments in the Amazon basin, and these differences cannot be reconciled by the fact that storm flow and subsurface outflows were not appropriately measured in previous studies.

Independent Review of Simulation of Net Infiltration for Present-Day and Potential Future Climates

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

Review Panel: Soroosh Sorooshian, Ph.D., Panel Chairperson, University of California, Irvine; Jan M. H. Hendrickx, Ph.D., New Mexico Institute of Mining and Technology; Binayak P. Mohanty, Ph.D., Texas A&M University

The DOE Office of Civilian Radioactive Waste Management (OCRWM) tasked Oak Ridge Institute for Science and Education (ORISE) with providing an independent expert review of the documented model and prediction results for net infiltration of water into the unsaturated zone at Yucca Mountain. The specific purpose of the model, as documented in the report MDL-NBS-HS-000023, Rev. 01, is “to provide a spatial representation, including epistemic and aleatory uncertainty, of the predicted mean annual net infiltration at the Yucca Mountain site ...” (p. 1-1) The expert review panel assembled by ORISE concluded that the model report does not provide a technicallymore » credible spatial representation of net infiltration at Yucca Mountain. Specifically, the ORISE Review Panel found that: • A critical lack of site-specific meteorological, surface, and subsurface information prevents verification of (i) the net infiltration estimates, (ii) the uncertainty estimates of parameters caused by their spatial variability, and (iii) the assumptions used by the modelers (ranges and distributions) for the characterization of parameters. The paucity of site-specific data used by the modeling team for model implementation and validation is a major deficiency in this effort. • The model does not incorporate at least one potentially important hydrologic process. Subsurface lateral flow is not accounted for by the model, and the assumption that the effect of subsurface lateral flow is negligible is not adequately justified. This issue is especially critical for the wetter climate periods. This omission may be one reason the model results appear to underestimate net infiltration beneath wash environments and therefore imprecisely represent the spatial variability of net infiltration. • While the model uses assumptions consistently, such as uniform soil depths and a constant vegetation rooting depth, such assumptions may not be appropriate for this net infiltration simulation because they oversimplify a complex landscape and associated hydrologic processes, especially since the model assumptions have not been adequately corroborated by field and laboratory observations at Yucca Mountain.« less

Coupled three-layer model for turbulent flow over large-scale roughness: On the hydrodynamics of boulder-bed streams

NASA Astrophysics Data System (ADS)

Pan, Wen-hao; Liu, Shi-he; Huang, Li

2018-02-01

This study developed a three-layer velocity model for turbulent flow over large-scale roughness. Through theoretical analysis, this model coupled both surface and subsurface flow. Flume experiments with flat cobble bed were conducted to examine the theoretical model. Results show that both the turbulent flow field and the total flow characteristics are quite different from that in the low gradient flow over microscale roughness. The velocity profile in a shallow stream converges to the logarithmic law away from the bed, while inflecting over the roughness layer to the non-zero subsurface flow. The velocity fluctuations close to a cobble bed are different from that of a sand bed, and it indicates no sufficiently large peak velocity. The total flow energy loss deviates significantly from the 1/7 power law equation when the relative flow depth is shallow. Both the coupled model and experiments indicate non-negligible subsurface flow that accounts for a considerable proportion of the total flow. By including the subsurface flow, the coupled model is able to predict a wider range of velocity profiles and total flow energy loss coefficients when compared with existing equations.

Interactions between bedforms, turbulence and pore flow

NASA Astrophysics Data System (ADS)

Blois, G.; Best, J.; Sambrook Smith, G.; Hardy, R. J.; Lead, J.

2010-12-01

A widespread occurrence of flow-form interaction in rivers is represented by subaqueous bedforms such as dunes. Many models have been proposed to explain how bedform generation and evolution are driven by turbulent flow structures that control the incipient motion of cohesionless sediments and later bedform development. However, most of these models have assumed such bedforms to be migrating over an impermeable bed, and that any surface-subsurface flow interaction is negligible. However, for some gravel-bed rivers the porosity can be high, up to 43%, which may result in significant flow both through the permeable bed (hyporheic flow) and across the surface-subsurface interface. The mass and momentum exchange occurring at the interface may have a strong impact on the structure of turbulent flow in the near-bed region. In the case of a dune, its topography induces a local pressure gradient that enhances flow across the interface. This results in a flow structure that may be radically different from that commonly proposed by past work. This paper presents results from a simplified laboratory model akin to a fine-grained bedform generated on top of a coarser sediment bed. Particle imaging velocimetry (PIV) measurements were conducted in order to characterise flow both over and underneath an idealised 2-dimensional dune (0.41 m long, 0.056 m high and having a leeside angle of 27°) overlaying a packed bed of uniform size spheres (D = 0.04 m diameter). Experiments were conducted in free surface flow conditions (Froude number = 0.1; Reynolds number = 25,000) for one bedform height: flow depth ratio (0.31). The flow above the dune was measured using a standard PIV technique while a novel endoscopic PIV (EPIV) system allowed collection of flow data within the pore spaces beneath the dune. The results show that topographically-induced subsurface flow significantly modifies the structure of flow in the leeside of the dune, resulting in a flow field that is radically different to traditional concepts of leeside flow. The pressure gradient across the bedform controls the direction and intensity of flow within the bed and across the interface. Specifically, the low pressure region induced by flow separation at the dune crest causes a pulsating jet flow from the bed into the free flow downstream of the dune. Fluid upwelling is particularly intense at the toe of the leeside and gradually decreases downstream. The interaction between the free-flow and hyporheic flow is significant; in the leeside, flow reattachment is entirely absent, and recirculation in the separation zone is replaced by a mechanism of asymmetric alternate vortex shedding. Hyporheic flow thus controls the dynamics of flow in the leeside and near-wake region. The paper discusses the implications of these results for the morphodynamics of coarse-sediment bedforms.

A porewater-based stable isotope approach for the investigation of subsurface hydrological processes

NASA Astrophysics Data System (ADS)

Garvelmann, J.; Külls, C.; Weiler, M.

2012-02-01

Predicting and understanding subsurface flowpaths is still a crucial issue in hydrological research. We present an experimental approach to reveal present and past subsurface flowpaths of water in the unsaturated and saturated zone. Two hillslopes in a humid mountainous catchment have been investigated. The H2O(liquid) - H2O(vapor) equilibration laser spectroscopy method was used to obtain high resolution δ2H vertical depth profiles of pore water at various points along two fall lines of a pasture hillslope in the southern Black Forest, Germany. The Porewater-based Stable Isotope Profile (PSIP) approach was developed to use the integrated information of several vertical depth profiles of deuterium along transects at the hillslope. Different shapes of depth profiles were observed in relation to hillslope position. The statistical variability (inter-quartile range and standard deviation) of each profile was used to characterize different types of depth profiles. The profiles upslope or with a weak affinity for saturation as indicated by a low topographic wetness index preserve the isotopic input signal by precipitation with a distinct seasonal variability. These observations indicate mainly vertical movement of soil water in the upper part of the hillslope before sampling. The profiles downslope or at locations with a strong affinity for saturation do not show a similar seasonal isotopic signal. The input signal is erased in the foothills and a large proportion of pore water samples are close to the isotopic values of δ2H in streamwater during base flow conditions indicating the importance of the groundwater component in the catchment. Near the stream indications for efficient mixing of water from lateral subsurface flow paths with vertical percolation are found.

Subsurface conditions in hydrothermal vents inferred from diffuse flow composition, and models of reaction and transport

NASA Astrophysics Data System (ADS)

Larson, B. I.; Houghton, J. L.; Lowell, R. P.; Farough, A.; Meile, C. D.

2015-08-01

Chemical gradients in the subsurface of mid-ocean ridge hydrothermal systems create an environment where minerals precipitate and dissolve and where chemosynthetic organisms thrive. However, owing to the lack of easy access to the subsurface, robust knowledge of the nature and extent of chemical transformations remains elusive. Here, we combine measurements of vent fluid chemistry with geochemical and transport modeling to give new insights into the under-sampled subsurface. Temperature-composition relationships from a geochemical mixing model are superimposed on the subsurface temperature distribution determined using a heat flow model to estimate the spatial distribution of fluid composition. We then estimate the distribution of Gibb's free energies of reaction beneath mid oceanic ridges and by combining flow simulations with speciation calculations estimate anhydrite deposition rates. Applied to vent endmembers observed at the fast spreading ridge at the East Pacific Rise, our results suggest that sealing times due to anhydrite formation are longer than the typical time between tectonic and magmatic events. The chemical composition of the neighboring low temperature flow indicates relatively uniform energetically favorable conditions for commonly inferred microbial processes such as methanogenesis, sulfate reduction and numerous oxidation reactions, suggesting that factors other than energy availability may control subsurface microbial biomass distribution. Thus, these model simulations complement fluid-sample datasets from surface venting and help infer the chemical distribution and transformations in subsurface flow.

Triaxial thermopile array geo-heat-flow sensor

DOEpatents

Carrigan, C.R.; Hardee, H.C.; Reynolds, G.D.; Steinfort, T.D.

1990-01-01

A triaxial thermopile array geothermal heat flow sensor is designed to measure heat flow in three dimensions in a reconstituted or unperturbed subsurface regime. Heat flow can be measured in conductive or permeable convective media. The sensor may be encased in protective pvc tubing and includes a plurality of thermistors and an array of heat flow transducers produce voltage proportional to heat flux along the subsurface regime and permit direct measurement of heat flow in the subsurface regime. The presence of the thermistor array permits a comparison to be made between the heat flow estimates obtained from the transducers and heat flow calculated using temperature differences and Fourier's Law. The device is extremely sensitive with an accuracy of less than 0.1 Heat Flow Units (HFU) and may be used for long term readings. 6 figs.

Evaluate the Relative Importance of Subsurface Lateral Energy Exchange to Ground Heat Flux and Energy Balance over the Heterogeneous Surface of a Sub-tropical Wetland

NASA Astrophysics Data System (ADS)

CUI, W.; Chui, T. F. M.

2016-12-01

Subsurface lateral water and energy exchanges are often ignored in methods involving a surface energy balance under the homogeneity assumption, which may affect the estimation of evapotranspiration over a heterogeneous surface. Wetlands, however, are heterogeneous with vegetated areas and open water, making it difficult to accurately measure and estimate evapotranspiration. This study estimated the subsurface lateral energy exchange between the reed bed and shallow open water of a wetland within Mai Po Nature Reserve in Hong Kong, and further discussed its relative importance to the ground heat flux and energy balance over the wetland surface. An array of water level and temperature sensors were installed in the reed bed and the adjacent water, together with an eddy covariance system. The results suggested that the lateral energy exchange was over 30% of ground heat flux for half of the monitoring period, and should therefore be accounted for during the measurement of ground heat flux. However, the lateral energy exchange could not explain the energy balance disclosure at the site, as the variation was in phase with the residual of energy budget during the summer but was out of phase during the winter. Furthermore, this study developed a convolution model to estimate the lateral energy exchange based on air temperature which is readily available at many sites worldwide. This study overall enhanced our understanding of the subsurface lateral energy exchange, and possibly our estimation of evapotranspiration in heterogeneous environment.

Micro Imaging Spectrometer for Subsurface Studies of Martian Soil: Ma_Miss

NASA Astrophysics Data System (ADS)

de Sanctis, M. C.; Coradini, A.; Ammannito, E.; Boccaccini, A.; di Iorio, T.; Battistelli, E.; Capanni, A.

2012-03-01

Ma_Miss (Mars Multispectral Imager for Subsurface Studies) is a spectrometer devoted to observe the lateral wall of the borehole generated by the drill installed on the ExoMars Pasteur Rover to perform in situ investigations in the Mars subsurface.

Regional variations in the stability and diffusion of water-ice in the Martian regolith

NASA Technical Reports Server (NTRS)

Mellon, Michael T.; Jakosky, Bruce M.

1992-01-01

Geologic evidence suggests subsurface water-ice has played an important role in the formation of Martian landforms. Forms of mass-wasting such as debris aprons and flow patterns on valley floors suggest creep deformation of ice-laden soil, while thermokarst and chaotic terrain suggest once extensive deposits of ground ice that were later removed. The global distribution of ice-related morphology was mapped. The mapping showed regional variation, in both latitude and longitude, in the distribution of debris aprons, concentric fill craters, and 'softened' crater profiles.

Subsurface Controls on Stream Intermittency in a Semi-Arid Landscape

NASA Astrophysics Data System (ADS)

Dohman, J.; Godsey, S.; Thackray, G. D.; Hale, R. L.; Wright, K.; Martinez, D.

2017-12-01

Intermittent streams currently constitute 30% to greater than 50% of the global river network. In addition, the number of intermittent streams is expected to increase due to changes in land use and climate. These streams provide important ecosystem services, such as water for irrigation, increased biodiversity, and high rates of nutrient cycling. Many hydrological studies have focused on mapping current intermittent flow regimes or evaluating long-term flow records, but very few have investigated the underlying causes of stream intermittency. The disconnection and reconnection of surface flow reflects the capacity of the subsurface to accommodate flow, so characterizing subsurface flow is key to understanding stream drying. We assess how subsurface flow paths control local surface flows during low-flow periods, including intermittency. Water table dynamics were monitored in an intermittent reach of Gibson Jack Creek in southeastern Idaho. Four transects were delineated with a groundwater well located in the hillslope, riparian zone, and in the stream, for a total of 12 groundwater wells. The presence or absence of surface flow was determined by frequent visual observations as well as in situ loggers every 30m along the 200m study reach. The rate of surface water drying was measured in conjunction with temperature, precipitation, subsurface hydraulic conductivity, hillslope-riparian-stream connectivity and subsurface travel time. Initial results during an unusually wet year suggest different responses in reaches that were previously observed to occasionally cease flowing. Flows in the intermittent reaches had less coherent and lower amplitude diel variations during base flow periods than reaches that had never been observed to dry out. Our findings will help contribute to our understanding of mechanisms driving expansion and contraction cycles in intermittent streams, increase our ability to predict how land use and climate change will affect flow regimes, and improve management of our critical water resources.

Removal of nutrients from septic tank effluent with baffle subsurface-flow constructed wetlands

Treesearch

Lihu Cui; Ying Ouyang; Weizhi Yang; Zhujian Huang; Qiaoling Xu; Guangwei Yu

2015-01-01

Three new baffle flow constructed wetlands (CWs), namely the baffle horizontal flow CW (Z1), baffle vertical flow CW (Z2) and baffle hybrid flow CW (Z3), along with one traditional horizontal subsurface flow CW (Z4) were designed to test the removal efficiency of nitrogen (N) and phosphorus (P) from the septic tank effluent under varying hydraulic retention times (HRTs...

Monitoring Subsurface Fluid Flow Using Perfluorocarbon Tracers: Another Tool Potentially Available for Subsurface Fluid Flow Assessments

EPA Pesticide Factsheets

Perfluorocarbon Tracers (PFTs) Complement stable Isotopes and Geochemistry for Verifying, Assessing or Modeling Fluid Flow. Geochemistry, Isotopes and PFT’s complement Geophysics to monitor and verify plume movement, leakage to shallow aquifers or surface

Lysimeter apparatus

DOEpatents

Clark, Don T.; Erickson, Eugene E.; Casper, William L.; Everett, David M.; Hubbell, Joel M.; Sisson, James B.

2005-09-06

A suction lysimeter for sampling subsurface liquids includes a lysimeter casing having a drive portion, a reservoir portion, and a tip portion, the tip portion including a membrane through which subsurface liquids may be sampled; a fluid conduit coupled in fluid flowing relation relative to the membrane, and which in operation facilitates the delivery of the sampled subsurface liquids from the membrane to the reservoir portion; and a plurality of tubes coupled in fluid flowing relation relative to the reservoir portion, the tubes in operation facilitating delivery of the sampled subsurface liquids from the reservoir portion for testing. A method of sampling subsurface liquids comprises using this lysimeter.

Joint application of Geoelectrical Resistivity and Ground Penetrating Radar techniques for the study of hyper-saturated zones. Case study in Egypt

NASA Astrophysics Data System (ADS)

Mesbah, Hany S.; Morsy, Essam A.; Soliman, Mamdouh M.; Kabeel, Khamis

2017-06-01

This paper presents the results of the application of the Geoelectrical Resistivity Sounding (GRS) and Ground Penetrating Radar (GPR) for outlining and investigating of surface springing out (flow) of groundwater to the base of an service building site, and determining the reason(s) for the zone of maximum degree of saturation; in addition to provide stratigraphic information for this site. The studied economic building is constructed lower than the ground surface by about 7 m. A Vertical Electrical Sounding (VES) survey was performed at 12 points around the studied building in order to investigate the vertical and lateral extent of the subsurface sequence, three VES's were conducted at each side of the building at discrete distances. And a total of 9 GPR profiles with 100- and 200-MHz antennae were conducted, with the objective of evaluating the depth and the degree of saturation of the subsurface layers. The qualitative and quantitative interpretation of the acquired VES's showed easily the levels of saturations close to and around the studied building. From the interpretation of GPR profiles, it was possible to locate and determine the saturated layers. The radar signals are penetrated and enabled the identification of the subsurface reflectors. The results of GPR and VES showed a good agreement and the integrated interpretations were supported by local geology. Finally, the new constructed geoelectrical resistivity cross-sections (in contoured-form), are easily clarifying the direction of groundwater flow toward the studied building.

Triaxial thermopile array geo-heat-flow sensor

DOEpatents

Carrigan, Charles R.; Hardee, Harry C.; Reynolds, Gerald D.; Steinfort, Terry D.

1992-01-01

A triaxial thermopile array geothermal heat flow sensor is designed to measure heat flow in three dimensions in a reconstituted or unperturbed subsurface regime. Heat flow can be measured in conductive or permeable convective media. The sensor may be encased in protective pvc tubing and includes a plurality of thermistors and an array of heat flow transducers arranged in a vertical string. The transducers produce voltage proportional to heat flux along the subsurface regime and permit direct measurement of heat flow in the subsurface regime. The presence of the thermistor array permits a comparison to be made between the heat flow estimates obtained from the transducers and heat flow calculated using temperature differences and Fourier's Law. The device is extremely sensitive with an accuracy of less than 0.1 Heat Flow Units (HFU) and may be used for long term readings.

Researchers Mine Information from Next-Generation Subsurface Flow Simulations

DOE PAGES

Gedenk, Eric D.

2015-12-01

A research team based at Virginia Tech University leveraged computing resources at the US Department of Energy's (DOE's) Oak Ridge National Laboratory to explore subsurface multiphase flow phenomena that can't be experimentally observed. Using the Cray XK7 Titan supercomputer at the Oak Ridge Leadership Computing Facility, the team took Micro-CT images of subsurface geologic systems and created two-phase flow simulations. The team's model development has implications for computational research pertaining to carbon sequestration, oil recovery, and contaminant transport.

Implementing a conceptual model of physical and chemical soil profile evolution

NASA Astrophysics Data System (ADS)

Kirkby, Mike

2017-04-01

When soil profile composition is generalised in terms of the proportion, p, of bedrock remaining (= 1 - depletion ratio), then other soil processes can also be expressed in terms of p, and 'soil depth' described by the integral of (1-p) down to bedrock. Soil profile evolution is expressed as the advance of a sigmoidal weathering front into the critical zone under the action of upward ionic diffusion of weathering products; downward advection of solutes in percolating waters, with loss of (cleanish) water as evapotranspiration and (solute-laden) water as a lateral sub-surface flow increment; and mechanical denudation increment at the surface. Each component responds to the degree of weathering. Percolation is limited by precipitation, evapotranspiration demand and the degree of weathering at each level in the profile which diverts subsurface flow. Mechanical removal rates are considered to broadly increase as weathering proceeds, as grain size and dilation angle decreases. The implication of these assumptions can be examined for steady state profiles, for which observed relationships between mechanical and chemical denudation rates; and between chemical denudation and critical zone depth are reproduced. For non-steady state evolution, these relationships break down, but provide a basis for linking critical zone with hillslope/ landform evolution.

Data on the solute concentration within the subsurface flows of Little Lost Man Creek in response to a transport experiment, Redwood National Park, northwest California

USGS Publications Warehouse

Zellweger, Gary W.; Kennedy, V.C.; Bencala, K.E.; Avanzino, R.J.; Jackman, A.P.; Triska, F.J.

1986-01-01

A solute transport experiment was conducted on a 327-m reach of Little Lost Man Creek, a small stream in Humboldt County, California. Solutes were injected for 20 days. Chloride was used as a conservative tracer; lithium, potassium, and strontium were used as reactive tracers. In addition, nitrate and phosphate were added as biological reactants. Eighteen shallow wells were dug along the length of the study reach, 1-10 m laterally from the edge of the stream. The wells and sites in the stream were monitored for the injected solutes during and after the injection. Solute concentrations in the wells and stream are indicative of transport properties of stream and subsurface channel flow. This report presents the results of the analyses of the well samples and chemical data relevant to the interpretation of hydrological and chemical interaction between the stream and adjacent channel subsurface flows in the streambed. Calculations of the percentage of streamwater in the wellwater were made from conservative tracer measurements. The composition of wellwater ranged from 47% to 100% streamwater with most values above 90%. The time for water to travel from the beginning of the study reach to the wells was approximately three times as great as the travel time in the stream at the same distance down the reach. The three conclusions that can be drawn are (1) water in the stream exchanges extensively with water in the rest of the channel; (2) the interstitial water in the channel gravels achieves almost the same composition as the stream; and (3) under low flow conditions the stream gravels contain a significant portion of the stream volume. Plots of normalized chloride, lithium, and strontium concentrations at three stream sites are included. (Author 's abstract)

The Influence of Plant Root Systems on Subsurface Flow: Implications for Slope Stability

EPA Science Inventory

Although research has explained how plant roots mechanically stabilize soils, in this article we explore how root systems create networks of preferential flow and thus influence water pressures in soils to trigger landslides. Root systems may alter subsurface flow: Hydrological m...

MODELING THREE-DIMENSIONAL SUBSURFACE FLOW, FATE AND TRANSPORT OF MICROBES AND CHEMICALS (3DFATMIC)

EPA Science Inventory

A three-dimensional model simulating the subsurface flow, microbial growth and degradation, microbial-chemical reaction, and transport of microbes and chemicals has been developed. he model is designed to solve the coupled flow and transport equations. asically, the saturated-uns...

Modelling rapid subsurface flow at the hillslope scale with explicit representation of preferential flow paths

NASA Astrophysics Data System (ADS)

Wienhöfer, J.; Zehe, E.

2012-04-01

Rapid lateral flow processes via preferential flow paths are widely accepted to play a key role for rainfall-runoff response in temperate humid headwater catchments. A quantitative description of these processes, however, is still a major challenge in hydrological research, not least because detailed information about the architecture of subsurface flow paths are often impossible to obtain at a natural site without disturbing the system. Our study combines physically based modelling and field observations with the objective to better understand how flow network configurations influence the hydrological response of hillslopes. The system under investigation is a forested hillslope with a small perennial spring at the study area Heumöser, a headwater catchment of the Dornbirnerach in Vorarlberg, Austria. In-situ points measurements of field-saturated hydraulic conductivity and dye staining experiments at the plot scale revealed that shrinkage cracks and biogenic macropores function as preferential flow paths in the fine-textured soils of the study area, and these preferential flow structures were active in fast subsurface transport of artificial tracers at the hillslope scale. For modelling of water and solute transport, we followed the approach of implementing preferential flow paths as spatially explicit structures of high hydraulic conductivity and low retention within the 2D process-based model CATFLOW. Many potential configurations of the flow path network were generated as realisations of a stochastic process informed by macropore characteristics derived from the plot scale observations. Together with different realisations of soil hydraulic parameters, this approach results in a Monte Carlo study. The model setups were used for short-term simulation of a sprinkling and tracer experiment, and the results were evaluated against measured discharges and tracer breakthrough curves. Although both criteria were taken for model evaluation, still several model setups produced acceptable matches to the observed behaviour. These setups were selected for long-term simulation, the results of which were compared against water level measurements at two piezometers along the hillslope and the integral discharge response of the spring to reject some non-behavioural model setups and further reduce equifinality. The results of this study indicate that process-based modelling can provide a means to distinguish preferential flow networks on the hillslope scale when complementary measurements to constrain the range of behavioural model setups are available. These models can further be employed as a virtual reality to investigate the characteristics of flow path architectures and explore effective parameterisations for larger scale applications.

Lysimeter methods and apparatus

DOEpatents

Clark, Don T.; Erickson, Eugene E.; Casper, William L.; Everett, David M.; Hubbell, Joel M.; Sisson, James B.

2004-12-07

A suction lysimeter for sampling subsurface liquids includes a lysimeter casing having a drive portion, a reservoir portion, and a tip portion, the tip portion including a membrane through which subsurface liquids may be sampled; a fluid conduit coupled in fluid flowing relation relative to the membrane, and which in operation facilitates the delivery of the sampled subsurface liquids from the membrane to the reservoir portion; and a plurality of tubes coupled in fluid flowing relation relative to the reservoir portion, the tubes in operation facilitating delivery of the sampled subsurface liquids from the reservoir portion for testing. A method of sampling subsurface liquids comprises using this lysimeter.

Using an Integrated Hydrologic Model to Assess the Ecohydrological Impacts of Change on a Mountain Headwaters Critical Zone

NASA Astrophysics Data System (ADS)

Collins, C.; Maxwell, R. M.; Visser, A.

2016-12-01

The critical zone is the region of the Earth's crust where hydrogeology, ecology, and climate interact. As many critical zone processes are fundamental, the significance of studying critical zone processes goes beyond understanding the local ecohydrological setting. Therefore studying critical zone governing processes requires an interdisciplinary approach that integrates simulation and observation. In this study, a high-resolution integrated hydrologic model, ParFlow-CLM, was developed for the Providence Creek watershed. Providence Creek is a highly instrumented critical zone observatory (CZO) located in the southern Sierra Nevada Mountains, a region currently experiencing a range of short-term responses (i.e. tree mortality) to a severe four-year drought. Sources of plant water use, pathways and residence times of water through the subsurface are identified using a suite of isotopic signatures and numerical particle tracking. Implications of using a fully coupled integrated hydrologic model accompanied by tracer analysis include better understanding of water partitioning and water storage in the regolith and vegetation water use during drought time conditions. The importance of subsurface storage, plant available water and lateral flow during the 2012-2015 drought to mitigate vegetation stress are addressed and verified against observed tree mortality. The stream flow response to tree mortality in the aftermath of the drought, analogous to the Colorado Mountain Pine Beetle case, provides insight into the potential effects of proposed forest management practices.

Subsurface energy storage and transport for solar-powered geysers on Triton

NASA Technical Reports Server (NTRS)

Kirk, Randolph L.; Soderblom, Laurence A.; Brown, Robert H.

1990-01-01

The location of active geyser-like eruptions and related features close to the current subsolar latitude on Triton suggests a solar energy source for these phenomena. Solid-state greenhouse calculations have shown that sunlight can generate substantially elevated subsurface temperatures. A variety of models for the storage of solar energy in a subgreenhouse layer and for the supply of gas and energy to a geyser are examined. 'Leaky greenhouse' models with only vertical gas transport are inconsistent with the observed upper limit on geyser radius of about 1.5 km. However, lateral transport of energy by gas flow in a porous N2 layer with a block size on the order of a meter can supply the required amount of gas to a source region about 1 km in radius. The decline of gas output to steady state may occur over a period comparable with the inferred active geyser lifetime of 5 earth years. The required subsurface permeability may be maintained by thermal fracturing of the residual N2 polar cap. A lower limit on geyser source radius of about 50 to 100 m predicted by a theory of negatively buoyant jets is not readily attained.

Modeling the time-varying interaction between surface water and groundwater bodies

NASA Astrophysics Data System (ADS)

Gliege, Steffen; Steidl, Jörg; Lischeid, Gunnar; Merz, Christoph

2016-04-01

The countless kettle holes (small lakes) in the Late Pleistocene landscapes of Northern Europe have important ecological and hydrological functions. On the one hand they act as depressions in which water and solutes of mainly agriculturally used catchments accumulate. On the other hand they operate as biochemical reactors with respect to greenhouse gas emissions, carbon sequestration, and as major sinks for nutrients and contaminants. Even small kettle holes often are hydraulically connected to the uppermost groundwater system: Groundwater discharges into the kettle hole on one side, and the aquifer is recharged from the kettle hole water body on the other side. Thus kettle hole biogeochemical processes are both affected by groundwater and vice versa. Groundwater flow direction and velocity into and out of the kettle hole often is not stable over time. Groundwater flow direction might reverse at the downstream part, resulting in repeated recycling of groundwater and corresponding solute turnover within the kettle holes. A sound understanding of this intricate interplay is a necessary prerequisite for better understanding of the biogeochemistry of this terrestrial-aquatic interface. A numerical experiment was used to quantify the lateral solute exchange between a kettle hole and the surrounding groundwater. A vertical cross section through the real existing catchment of a kettle hole was chosen. Glacial till represents the lower boundary. The heterogeneity of the subsurface was reproduced by various parameterizations of the soil hydraulic properties as well as varying the thickness of the unconfined aquifer or the lateral boundary conditions. In total 24 different parameterizations were implemented in the modeling software HydroGeoSphere (HGS). HGS is suitable to calculate the fluid exchange between surface and subsurface simultaneously and in a physically based way. The simulation runs were done for the period from November 1994 to October 2014. All results were analyzed with regard to the intensity and duration of water exchange between kettle hole and surrounding groundwater. Finally the three variants with the smallest, average and largest number of days where water flow is directed from kettle hole to groundwater were chosen to extend the scope of this study to include solute transport. Therefore a non-reactive substance was injected at the kettle hole bed. As a result, the fluid dynamics and the spread of water flowing from surface to subsurface and vice versa could be closely monitored. Finally, the choice of these variants provides information on the range of distances and duration of water exchange between kettle hole and surrounding groundwater.

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.

Let's Go Off the Grid: Subsurface Flow Modeling With Analytic Elements

NASA Astrophysics Data System (ADS)

Bakker, M.

2017-12-01

Subsurface flow modeling with analytic elements has the major advantage that no grid or time stepping are needed. Analytic element formulations exist for steady state and transient flow in layered aquifers and unsaturated flow in the vadose zone. Analytic element models are vector-based and consist of points, lines and curves that represent specific features in the subsurface. Recent advances allow for the simulation of partially penetrating wells and multi-aquifer wells, including skin effect and wellbore storage, horizontal wells of poly-line shape including skin effect, sharp changes in subsurface properties, and surface water features with leaky beds. Input files for analytic element models are simple, short and readable, and can easily be generated from, for example, GIS databases. Future plans include the incorporation of analytic element in parts of grid-based models where additional detail is needed. This presentation will give an overview of advanced flow features that can be modeled, many of which are implemented in free and open-source software.

Development of stream-subsurface flow module in sub-daily simulation of Escherichia coli using SWAT

NASA Astrophysics Data System (ADS)

Kim, Minjeong; Boithias, Laurie; Cho, Kyung Hwa; Silvera, Norbert; Thammahacksa, Chanthamousone; Latsachack, Keooudone; Rochelle-Newall, Emma; Sengtaheuanghoung, Oloth; Pierret, Alain; Pachepsky, Yakov A.; Ribolzi, Olivier

2017-04-01

Water contaminated with pathogenic bacteria poses a large threat to public health, especially in the rural areas in the tropics where sanitation and drinking water facilities are often lacking. Several studies have used the Soil and Water Assessment Tool (SWAT) to predict the export of in-stream bacteria at a watershed-scale. However, SWAT is limited to in-stream processes, such as die-off, resuspension and, deposition; and it is usually implemented on a daily time step using the SCS Curve Number method, making it difficult to explore the dynamic fate and transport of bacteria during short but intense events such as flash floods in tropical humid montane headwaters. To address these issues, this study implemented SWAT on an hourly time step using the Green-Ampt infiltration method, and tested the effects of subsurface flow (LATQ+GWQ in SWAT) on bacterial dynamics. We applied the modified SWAT model to the 60-ha Houay Pano catchment in Northern Laos, using sub-daily rainfall and discharge measurements, electric conductivity-derived fractions of overland and subsurface flows, suspended sediments concentrations, and the number of fecal indicator organism Escherichia coli monitored at the catchment outlet from 2011 to 2013. We also took into account land use change by delineating the watershed with the 3-year composite land use map. The results show that low subsurface flow of less than 1 mm recovered the underestimation of E. coli numbers during the dry season, while high subsurface flow caused an overestimation during the wet season. We also found that it is more reasonable to apply the stream-subsurface flow interaction to simulate low in-stream bacteria counts. Using fecal bacteria to identify and understand the possible interactions between overland and subsurface flows may well also provide some insight into the fate of other bacteria, such as those involved in biogeochemical fluxes both in-stream and in the adjacent soils and hyporheic zones.

Three-Dimensional Subsurface Flow, Fate and Transport of Microbes and Chemicals (3DFATMIC) Model

EPA Pesticide Factsheets

This model simulates subsurface flow, fate and transport of contaminants that are undergoing chemical or biological transformations. The model is applicable to transient conditions in both saturated and unsaturated zones.

Two-Dimensional Subsurface Flow, Fate and Transport of Microbes and Chemicals (2DFATMIC) Model

EPA Pesticide Factsheets

This model simulates subsurface flow, fate, and transport of contaminants that are undergoing chemical or biological transformations. This model is applicable to transient conditions in both saturated and unsaturated zones.

Upscaling of reaction rates in reactive transport using pore-scale reactive transport model

NASA Astrophysics Data System (ADS)

Yoon, H.; Dewers, T. A.; Arnold, B. W.; Major, J. R.; Eichhubl, P.; Srinivasan, S.

2013-12-01

Dissolved CO2 during geological CO2 storage may react with minerals in fractured rocks, confined aquifers, or faults, resulting in mineral precipitation and dissolution. The overall rate of reaction can be affected by coupled processes among hydrodynamics, transport, and reactions at the (sub) pore-scale. In this research pore-scale modeling of coupled fluid flow, reactive transport, and heterogeneous reaction at the mineral surface is applied to account for permeability alterations caused by precipitation-induced pore-blocking. This work is motivated by the observed CO2 seeps from a natural analog to geologic CO2 sequestration at Crystal Geyser, Utah. A key observation is the lateral migration of CO2 seep sites at a scale of ~ 100 meters over time. A pore-scale model provides fundamental mechanistic explanations of how calcite precipitation alters flow paths by pore plugging under different geochemical compositions and pore configurations. In addition, response function of reaction rates will be constructed from pore-scale simulations which account for a range of reaction regimes characterized by the Damkohler and Peclet numbers. Newly developed response functions will be used in a continuum scale model that may account for large-scale phenomena mimicking lateral migration of surface CO2 seeps. Comparison of field observations and simulations results will provide mechanistic explanations of the lateral migration and enhance our understanding of subsurface processes associated with the CO2 injection. This work is supported as part of the Center for Frontiers of Subsurface Energy Security, an Energy Frontier Research Center funded by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences under Award Number DE-SC0001114. Sandia National Laboratories is a multi-program laboratory managed and operated by Sandia Corporation, a wholly owned subsidiary of Lockheed Martin Corporation, for the U.S. Department of Energy's National Nuclear Security Administration under contract DE-AC04-94AL85000.

Effects Of Bedrock Shape And Hillslope Gradient On The Pore-Water Pressure Development: Implication For Slope Stability

NASA Astrophysics Data System (ADS)

Lanni, Cristiano; McDonnell, Jeff

2010-05-01

Shallow Landslides are one of the most important causes of loss of human life and socio-economic damage related to the hydro-geological risk issues. The danger of these phenomena is related to their speed of development, the diffculty of foreseeing their location, and the high density of individual phenomena, whose downhill trajectories have a relevant probability of interfering with urbanized areas. Research activity on precipitation-induced landslides has focused mainly on developing predictive understanding of where and when landslides are likely to occur. Nevertheless, some major aspects that may be related to activation of landslides have been poorly investigated. For instance, landslide susceptibility zones are generally predicted assuming constant thickness of soil over an impervious bedrock layer. Nevertheless, recent studies showed subsurface topography could be a first order control for subsurface water-flow dynamics, because of the effects of its own irregular shape. Tromp-van Meerveld and McDonnell (2006) argued that connectivity of patches of transient saturation were a necessary prerequisite for exceeding the rainfall threshold necessary to drive lateral flow. Connectivity - "how the hillslope architecture controls the filling and spilling of isolated patches of saturation" (Hopp and McDonnell, 2009) - appears to be a possible unifying concept and theoretical platform for moving hillslope and watershed hydrology forward. Connectivity could also have important implications on triggering of shallow landslides, because the particular shape of bedrock may limit the water-flow downhill. Here we present a number of virtual numerical experiments performed to investigate the role of bedrock shape and hillslope gradient on pore-water pressure development. On this purpose, our test is represented by the subsurface topography of the Panola Experiment Hillslope (PEH). That is because scientific literature on PEH provides substantial documentation about the role of bedrock layer on subsurface water-flow dynamics. We also exploit the concept of Downslope Index (DWI) (Hjerdt et al., 2004) and Upslope Contributing Area (UCA) as indicators of the areas more susceptible to landslide. The results indicate that bedrock shape influences the max pore-water pressure, even with different hillslope gradients; meanwhile, hillslope gradient affects the persistence-time of the max pore-water pressure. Moreover, results suggest DWI as an useful index to improve the capability of the very-used SHALSTAB model to assess for landslide susceptibility areas.

Experimental quantification of solute transport through the vadose zone under dynamic boundary conditions with dye tracers and optical methods.

NASA Astrophysics Data System (ADS)

Cremer, Clemens; Neuweiler, Insa

2017-04-01

Knowledge of subsurface solute transport processes is vital to investigate e.g. groundwater contamination, nutrient uptake by plant roots and to implement remediation strategies. Beside field measurements and numerical simulations, physical laboratory experiments represent a way to establish process understanding and furthermore validate numerical schemes. Atmospheric forcings, such as erratically varying infiltration and evaporation cycles, subject the shallow subsurface to local and temporal variations in water content and associated hydraulic conductivity of the prevailing porous media. Those variations in material properties can cause flow paths to differ between upward and downward flow periods. Thereby, the unsaturated subsurface presents a highly complicated, dynamic system. Following an extensive systematical numerical investigation of flow and transport through bimodal, unsaturated porous media under dynamic boundary conditions (Cremer et al., 2016), we conduct physical laboratory experiments in a 22 cm x 8 cm x 1 cm flow cell where we introduce structural heterogeneity in the form sharp material interfaces between different porous media. In all experiments, a constant pressure head is implemented at the lower boundary, while cyclic infiltration-evaporation phases are applied at the soil surface. As a reference case a stationary infiltration with a rate corresponding to the cycle-averaged infiltration rate is applied. By initial application of dye tracers, solute transport within the domain is visualized such that transport paths and redistribution processes can be observed in a qualitative manner. Solute leaching is quantified at the bottom outlet, where breakthrough curves are obtained via spectroscopy. Liquid and vapor flow in and out of the domain is obtained from multiple balances. Thereby, the interplay of material structural heterogeneity and alternating flow (transport) directions and flow (transport) paths is investigated. Results show lateral transport through the material interface which differs between the stationary (unilateral) and dynamic cases (bilateral). This qualitative observation is confirmed by breakthrough curves for dynamic experiments which generally show the trend of faster initial breakthrough and increased tailing when compared to stationary infiltration results. Literature Cremer, C.J.M., I. Neuweiler, M. Bechtold, J. Vanderborght (2016): Solute Transport in Heterogeneous Soil with Time-Dependent Boundary Conditions, Vadose Zone Journal 15 (6) DOI: 10.2136/vzj2015.11.0144

Treatment of laboratory wastewater in a tropical constructed wetland comparing surface and subsurface flow.

PubMed

Meutia, A A

2001-01-01

Wastewater treatment by constructed wetland is an appropriate technology for tropical developing countries like Indonesia because it is inexpensive, easily maintained, and has environmentally friendly and sustainable characteristics. The aim of the research is to examine the capability of constructed wetlands for treating laboratory wastewater at our Center, to investigate the suitable flow for treatment, namely vertical subsurface or horizontal surface flow, and to study the effect of the seasons. The constructed wetland is composed of three chambered unplanted sedimentation tanks followed by the first and second beds, containing gravel and sand, planted with Typha sp.; the third bed planted with floating plant Lemna sp.; and a clarifier with two chambers. The results showed that the subsurface flow in the dry season removed 95% organic carbon (COD) and total phosphorus (T-P) respectively, and 82% total nitrogen (T-N). In the transition period from the dry season to the rainy season, COD removal efficiency decreased to 73%, T-N increased to 89%, and T-P was almost the same as that in the dry season. In the rainy season COD and T-N removal efficiencies increased again to 95% respectively, while T-P remained unchanged. In the dry season, COD and T-P concentrations in the surface flow showed that the removal efficiencies were a bit lower than those in the subsurface flow. Moreover, T-N removal efficiency was only half as much as that in the subsurface flow. However, in the transition period, COD removal efficiency decreased to 29%, while T-N increased to 74% and T-P was still constant, around 93%. In the rainy season, COD and T-N removal efficiencies increased again to almost 95%. On the other hand, T-P decreased to 76%. The results show that the constructed wetland is capable of treating the laboratory wastewater. The subsurface flow is more suitable for treatment than the surface flow, and the seasonal changes have effects on the removal efficiency.

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

Aleman, S.E.

This report documents a finite element code designed to model subsurface flow and contaminant transport, named FACT. FACT is a transient three-dimensional, finite element code designed to simulate isothermal groundwater flow, moisture movement, and solute transport in variably saturated and fully saturated subsurface porous media.

Sources and processes affecting the distribution of dissolved Nd isotopes and concentrations in the West Pacific

NASA Astrophysics Data System (ADS)

Behrens, Melanie K.; Pahnke, Katharina; Schnetger, Bernhard; Brumsack, Hans-Jürgen

2018-02-01

In the Atlantic, where deep circulation is vigorous, the dissolved neodymium (Nd) isotopic composition (expressed as ɛNd) is largely controlled by water mass mixing. In contrast, the factors influencing the ɛNd distribution in the Pacific, marked by sluggish circulation, is not clear yet. Indication for regional overprints in the Pacific is given based on its bordering volcanic islands. Our study aims to clarify the impact and relative importance of different Nd sources (rivers, volcanic islands), vertical (bio)geochemical processes and lateral water mass transport in controlling dissolved ɛNd and Nd concentration ([Nd]) distributions in the West Pacific between South Korea and Fiji. We find indication for unradiogenic continental input from South Korean and Chinese rivers to the East China Sea. In the tropical West Pacific, volcanic islands supply Nd to surface and subsurface waters and modify their ɛNd to radiogenic values of up to +0.7. These radiogenic signatures allow detailed tracing of currents flowing to the east and differentiation from westward currents with open ocean Pacific ɛNd composition in the complex tropical Pacific zonal current system. Modified radiogenic ɛNd of West Pacific intermediate to bottom waters upstream or within our section also indicates non-conservative behavior of ɛNd due to boundary exchange at volcanic island margins, submarine ridges, and with hydrothermal particles. Only subsurface to deep waters (3000 m) in the open Northwest Pacific show conservative behavior of ɛNd. In contrast, we find a striking correlation of extremely low (down to 2.77 pmol/kg Nd) and laterally constant [Nd] with the high-salinity North and South Pacific Tropical Water, indicating lateral transport of preformed [Nd] from the North and South Pacific subtropical gyres into the study area. This observation also explains the previously observed low subsurface [Nd] in the tropical West Pacific. Similarly, Western South Pacific Central Water, Antarctic Intermediate Water, and Lower Circumpolar Deep Water in the southern and equatorial West Pacific are marked by vertically and laterally almost invariant [Nd] indicating a dominance of conservative behavior of [Nd]. In contrast, Central and Intermediate Water in the North West Pacific are characterized by increasing [Nd] with depth reflecting Nd release from particles. Overall, our data demonstrate a dominant lateral transport control on [Nd] distributions and clear non-conservative modification of ɛNd in the West Pacific. The latter affords tracing of surface and subsurface zonal transport in the tropical Pacific, but prevents the use of ɛNd as strictly conservative tracer of the major meridionally circulating water masses in the West Pacific between 15°S and 28°N.

Contribution of lateral terrestrial water flows to the regional hydrological cycle: A joint soil-atmospheric moisture tagging procedure with WRF-Hydro

NASA Astrophysics Data System (ADS)

Arnault, Joel; Wei, Jianhui; Zhang, Zhenyu; Wagner, Sven; Kunstmann, Harald

2017-04-01

Water resources management requires an accurate knowledge of the behavior of the regional hydrological cycle components, including precipitation, evapotranspiration, river discharge and soil water storage. Atmospheric models such as the Weather Research and Forecasting (WRF) model provide a tool to evaluate these components. The main drawback of these atmospheric models, however, is that the terrestrial segment of the hydrological cycle is reduced to vertical infiltration, and that lateral terrestrial water flows are neglected. Recent model developments have focused on coupled atmospheric-hydrological modeling systems, such as WRF-hydro, in order to take into account subsurface, overland and river flow. The aim of this study is to investigate the contribution of lateral terrestrial water flows to the regional hydrological cycle, with the help of a joint soil-atmospheric moisture tagging procedure. This procedure is the extended version of an existing atmospheric moisture tagging method developed in WRF and WRF-Hydro (Arnault et al. 2017). It is used to quantify the partitioning of precipitation into water stored in the soil, runoff, evapotranspiration, and potentially subsequent precipitation through regional recycling. An application to a high precipitation event on 23 June 2009 in the upper Danube river basin, Germany and Austria, is presented. Precipitating water during this day is tagged for the period 2009-2011. Its contribution to runoff and evapotranspiration decreases with time, but is still not negligible in the summer 2011. At the end of the study period, less than 5 % of the precipitating water on 23 June 2009 remains in the soil. The additionally resolved lateral terrestrial water flows in WRF-Hydro modify the partitioning between surface and underground runoff, in association with a slight increase of evapotranspiration and recycled precipitation. Reference: Arnault, J., R. Knoche, J. Wei, and H. Kunstmann (2016), Evaporation tagging and atmospheric water budget analysis with WRF: A regional precipitation recycling study for West Africa, Water Resour. Res., 52, 1544-1567, doi:10.1002/2015WR017704.

The thermal impact of subsurface building structures on urban groundwater resources - A paradigmatic example.

PubMed

Epting, Jannis; Scheidler, Stefan; Affolter, Annette; Borer, Paul; Mueller, Matthias H; Egli, Lukas; García-Gil, Alejandro; Huggenberger, Peter

2017-10-15

Shallow subsurface thermal regimes in urban areas are increasingly impacted by anthropogenic activities, which include infrastructure development like underground traffic lines as well as industrial and residential subsurface buildings. In combination with the progressive use of shallow geothermal energy systems, this results in the so-called subsurface urban heat island effect. This article emphasizes the importance of considering the thermal impact of subsurface structures, which commonly is underestimated due to missing information and of reliable subsurface temperature data. Based on synthetic heat-transport models different settings of the urban environment were investigated, including: (1) hydraulic gradients and conductivities, which result in different groundwater flow velocities; (2) aquifer properties like groundwater thickness to aquitard and depth to water table; and (3) constructional features, such as building depths and thermal properties of building structures. Our results demonstrate that with rising groundwater flow velocities, the heat-load from building structures increase, whereas down-gradient groundwater temperatures decrease. Thermal impacts on subsurface resources therefore have to be related to the permeability of aquifers and hydraulic boundary conditions. In regard to the urban settings of Basel, Switzerland, flow velocities of around 1 md -1 delineate a marker where either down-gradient temperature deviations or heat-loads into the subsurface are more relevant. Furthermore, no direct thermal influence on groundwater resources should be expected for aquifers with groundwater thicknesses larger 10m and when the distance of the building structure to the groundwater table is higher than around 10m. We demonstrate that measuring temperature changes down-gradient of subsurface structures is insufficient overall to assess thermal impacts, particularly in urban areas. Moreover, in areas which are densely urbanized, and where groundwater flow velocities are low, appropriate measures for assessing thermal impacts should specifically include a quantification of heat-loads into the subsurface which result in a more diffuse thermal contamination of urban groundwater resources. Copyright © 2017 Elsevier B.V. All rights reserved.

Coupling a three-dimensional subsurface flow model with a land surface model to simulate stream-aquifer-land interactions

NASA Astrophysics Data System (ADS)

Huang, M.; Bisht, G.; Zhou, T.; Chen, X.; Dai, H.; Hammond, G. E.; Riley, W. J.; Downs, J.; Liu, Y.; Zachara, J. M.

2016-12-01

A fully coupled three-dimensional surface and subsurface land model is developed and applied to a site along the Columbia River to simulate three-way interactions among river water, groundwater, and land surface processes. The model features the coupling of the Community Land Model version 4.5 (CLM4.5) and a massively-parallel multi-physics reactive tranport model (PFLOTRAN). The coupled model (CLM-PFLOTRAN) is applied to a 400m×400m study domain instrumented with groundwater monitoring wells in the Hanford 300 Area along the Columbia River. CLM-PFLOTRAN simulations are performed at three different spatial resolutions over the period 2011-2015 to evaluate the impact of spatial resolution on simulated variables. To demonstrate the difference in model simulations with and without lateral subsurface flow, a vertical-only CLM-PFLOTRAN simulation is also conducted for comparison. Results show that the coupled model is skillful in simulating stream-aquifer interactions, and the land-surface energy partitioning can be strongly modulated by groundwater-river water interactions in high water years due to increased soil moisture availability caused by elevated groundwater table. In addition, spatial resolution does not seem to impact the land surface energy flux simulations, although it is a key factor for accurately estimating the mass exchange rates at the boundaries and associated biogeochemical reactions in the aquifer. The coupled model developed in this study establishes a solid foundation for understanding co-evolution of hydrology and biogeochemistry along the river corridors under historical and future hydro-climate changes.

The effects of spatial heterogeneity and subsurface lateral transfer on evapotranspiration estimates in large scale Earth system models

NASA Astrophysics Data System (ADS)

Rouholahnejad, E.; Fan, Y.; Kirchner, J. W.; Miralles, D. G.

2017-12-01

Most Earth system models (ESM) average over considerable sub-grid heterogeneity in land surface properties, and overlook subsurface lateral flow. This could potentially bias evapotranspiration (ET) estimates and has implications for future temperature predictions, since overestimations in ET imply greater latent heat fluxes and potential underestimation of dry and warm conditions in the context of climate change. Here we quantify the bias in evaporation estimates that may arise from the fact that ESMs average over considerable heterogeneity in surface properties, and also neglect lateral transfer of water across the heterogeneous landscapes at global scale. We use a Budyko framework to express ET as a function of P and PET to derive simple sub-grid closure relations that quantify how spatial heterogeneity and lateral transfer could affect average ET as seen from the atmosphere. We show that averaging over sub-grid heterogeneity in P and PET, as typical Earth system models do, leads to overestimation of average ET. Our analysis at global scale shows that the effects of sub-grid heterogeneity will be most pronounced in steep mountainous areas where the topographic gradient is high and where P is inversely correlated with PET across the landscape. In addition, we use the Total Water Storage (TWS) anomaly estimates from the Gravity Recovery and Climate Experiment (GRACE) remote sensing product and assimilate it into the Global Land Evaporation Amsterdam Model (GLEAM) to correct for existing free drainage lower boundary condition in GLEAM and quantify whether, and how much, accounting for changes in terrestrial storage can improve the simulation of soil moisture and regional ET fluxes at global scale.

Study on removal efficiency of nitrogen and phosphorus from agricultural wastewater by subsurface flow constructed wetland

NASA Astrophysics Data System (ADS)

Ling, Zhen; Li, Jie

2018-03-01

Subsurface Flow Constructed Wetland Plant 5 kinds of perennial herbs, there are Canna, Water onion, Iris, Calamus, Reed. Foucs on Subsurface Flow Constructed Wetlands on agricultural wastewater nitrogen and phosphorus removal effect. Research results: Different plants TP removal efficiency from high to low is Iris> reed> calamus> water onion> canna.And TN removal efficiency from high to low is reed> water onion> iris> calamus> canna. Compared with the blank test land, Wetland plants improves TN removal and TP removal is higher than TN. Wetland plants can reduce the PH of experimental water.

Contribution of thermal infrared images on the understanding of the subsurface/atmosphere exchanges on Earth.

NASA Astrophysics Data System (ADS)

Lopez, Teodolina; Antoine, Raphaël; Baratoux, David; Rabinowicz, Michel

2017-04-01

High temporal resolution of space-based thermal infrared images (METEOSAT, MODIS) and the development of field thermal cameras have permitted the development of thermal remote sensing in Earth Sciences. Thermal images are influenced by many factors such as atmosphere, solar radiation, topography and physico-chemical properties of the surface. However, considering these limitations, we have discovered that thermal images can be used in order to better understand subsurface hydrology. In order to reduce as much as possible the impact of these perturbing factors, our approach combine 1) field observations and 2) numerical modelling of surface/subsurface thermal processes. Thermal images of the Piton de la Fournaise volcano (Réunion Island), acquired by hand, show that the Formica Leo inactive scoria cone and some fractures close to the Bory-Dolomieu caldera are always warmer, inducing a thermal difference with the surrounding of at least 5°C and a Self-Potential anomaly [1, 2]. Topography cannot explain this thermal behaviour, but Piton de la Fournaise is known as highly permeable. This fact allows the development of an air convection within the whole permeable structure volcanic edifice [2]. Cold air enters the base of the volcano, and exits warmer upslope, as the air is warmed by the geothermal flow [1,2]. Then, we have decided to understand the interaction between subsurface hydrogeological flows and the humidity in the atmosphere. In the Lake Chad basin, regions on both sides of Lake Chad present a different thermal behaviour during the diurnal cycle and between seasons [3]. We propose that this thermal behaviour can only be explained by lateral variations of the surface permeability that directly impact the process of evaporation/condensation cycle. These studies bring new highlights on the understanding of the exchanges between subsurface and the atmosphere, as the presence of a very permeable media and/or variations of the surface permeability may enhance or not the evaporation/condensation cycle. [1] Antoine et al. (2009). J. Volcanol. Geotherm. Res., 183(3-4), 228-1140. [2] Antoine et al. (2017). Geothermics, 65, 81-98. [3] Lopez et al. (2016). Surv. Geophys., 37 (2), 471-502.

Time-lapse 3-D seismic imaging of shallow subsurface contaminant flow.

PubMed

McKenna, J; Sherlock, D; Evans, B

2001-12-01

This paper presents a physical modelling study outlining a technique whereby buoyant contaminant flow within water-saturated unconsolidated sand was remotely monitored utilizing the time-lapse 3-D (TL3-D) seismic response. The controlled temperature and pressure conditions, along with the high level of acquisition repeatability attainable using sandbox physical models, allow the TL3-D seismic response to pore fluid movement to be distinguished from all other effects. TL3-D seismic techniques are currently being developed to monitor hydrocarbon reserves within producing reservoirs in an endeavour to improve overall recovery. However, in many ways, sandbox models under atmospheric conditions more accurately simulate the shallow subsurface than petroleum reservoirs. For this reason, perhaps the greatest application for analogue sandbox modelling is to improve our understanding of shallow groundwater and environmental flow mechanisms. Two fluid flow simulations were conducted whereby air and kerosene were injected into separate water-saturated unconsolidated sand models. In both experiments, a base 3-D seismic volume was recorded and compared with six later monitor surveys recorded while the injection program was conducted. Normal incidence amplitude and P-wave velocity information were extracted from the TL3-D seismic data to provide visualization of contaminant migration. Reflection amplitudes displayed qualitative areal distribution of fluids when a suitable impedance contrast existed between pore fluids. TL3-D seismic reflection tomography can potentially monitor the change in areal distribution of fluid contaminants over time, indicating flow patterns. However, other research and this current work have not established a quantifiable relationship between either normal reflection amplitudes and attenuation and fluid saturation. Generally, different pore fluids will have unique seismic velocities due to differences in compressibility and density. The predictable relationships that exist between P-wave velocity and fluid saturation can allow a quantitative assessment of contaminant migration.

The Effect of Projectile Density and Disruption on the Crater Excavation Flow-Field

NASA Technical Reports Server (NTRS)

Anderson, Jennifer L. B.; Schultz, P. H.

2005-01-01

The ejection parameters of material excavated by a growing crater directly relate to the subsurface excavation flow-field. The ejection angles and speeds define the end of subsurface material streamlines at the target surface. Differences in the subsurface flow-fields can be inferred by comparing observed ejection parameters of various impacts obtained using three-dimensional particle image velocimetry (3D PIV). The work presented here investigates the observed ejection speeds and angles of material ejected during vertical (90 impact angle) experimental impacts for a range of different projectile types. The subsurface flow-fields produced during vertical impacts are simple when compared with that of oblique impacts, affected primarily by the depth of the energy and momentum deposition of the projectile. This depth is highly controlled by the projectile/target density ratio and the disruption of the projectile (brittle vs. ductile deformation). Previous studies indicated that cratering efficiency and the crater diameter/depth ratio were affected by projectile disruption, velocity, and the projectile/target density ratio. The effect of these projectile properties on the excavation flow-field are examined by comparing different projectile materials.

Impacts of ditch blocking on peatland hydrology - the benefits of long-term monitoring

NASA Astrophysics Data System (ADS)

Holden, Joseph; Green, Sophie; Baird, Andy; Chapman, Pippa; Evans, Chris; Grayson, Richard

2016-04-01

A long-term field trial was conducted on a blanket peatland in North Wales. Twelve ditches were studied. After an initial monitoring period, eight of the ditches had peat dams installed a few metres apart along their entire length (dammed), four of these ditches were also partially infilled through bank reprofiling (reprofiled). Four ditches were left open with no dams or reprofiling (open). These 12 ditches and the surrounding peat were then monitored for a further 4 years. The effect of ditch blocking on local water tables was spatially highly variable but small overall (of the order of 2-3 cm) because the site, despite having ditches, already had relatively shallow water tables (medians within the upper 10 cm of the peat profile). An initial five-fold reduction in discharge occurred in ditches that had been dammed or reprofiled. However, there was evidence of a slow change over time in ditch flow at the site in subsequent years, with the overall volume of water leaving the dammed or reprofiled ditch weirs increasing per unit of rainfall to around twice that which occurred in the first year after the restoration. These changes were not observed in the open ditches. There was therefore clear evidence of the benefits of long-term monitoring as hydrological impacts in the first year after ditch blocking were very different from those in later years as the site conditions gradually changed. The additional water that flowed in later periods of the study from the blocked ditch catchments occurred in the form of a more continuously-flowing baseflow with fewer dry periods. The cause of this increase was related to changes in subsurface flow pathways in the peat in the aftermath of re-wetting. We show that these subsurface pathways mean that even in sloping blanket peatlands, the catchment areas for peatland ditches may be very different from that expressed by surface topography alone. Therefore, peatland studies that have estimated aerially-weighted water or carbon fluxes from one or two open or blocked ditches and where such data have also been used in upscaling estimates, need to be treated with caution.

Post-injection Multiphase Flow Modeling and Risk Assessments for Subsurface CO2 Storage in Naturally Fractured Reservoirs

NASA Astrophysics Data System (ADS)

Jin, G.

2015-12-01

Subsurface storage of carbon dioxide in geological formations is widely regarded as a promising tool for reducing global atmospheric CO2 emissions. Successful geologic storage for sequestrated carbon dioxides must prove to be safe by means of risk assessments including post-injection analysis of injected CO2 plumes. Because fractured reservoirs exhibit a higher degree of heterogeneity, it is imperative to conduct such simulation studies in order to reliably predict the geometric evolution of plumes and risk assessment of post CO2injection. The research has addressed the pressure footprint of CO2 plumes through the development of new techniques which combine discrete fracture network and stochastic continuum modeling of multiphase flow in fractured geologic formations. A subsequent permeability tensor map in 3-D, derived from our preciously developed method, can accurately describe the heterogeneity of fracture reservoirs. A comprehensive workflow integrating the fracture permeability characterization and multiphase flow modeling has been developed to simulate the CO2plume migration and risk assessments. A simulated fractured reservoir model based on high-priority geological carbon sinks in central Alabama has been employed for preliminary study. Discrete fracture networks were generated with an NE-oriented regional fracture set and orthogonal NW-fractures. Fracture permeability characterization revealed high permeability heterogeneity with an order of magnitude of up to three. A multiphase flow model composed of supercritical CO2 and saline water was then applied to predict CO2 plume volume, geometry, pressure footprint, and containment during and post injection. Injection simulation reveals significant permeability anisotropy that favors development of northeast-elongate CO2 plumes, which are aligned with systematic fractures. The diffusive spreading front of the CO2 plume shows strong viscous fingering effects. Post-injection simulation indicates significant upward lateral spreading of CO2 resulting in accumulation of CO2 directly under the seal unit because of its buoyancy and strata-bound vertical fractures. Risk assessment shows that lateral movement of CO2 along interconnected fractures requires widespread seals with high integrity to confine the injected CO2.

Sequential geophysical and flow inversion to characterize fracture networks in subsurface systems

DOE PAGES

Mudunuru, Maruti Kumar; Karra, Satish; Makedonska, Nataliia; ...

2017-09-05

Subsurface applications, including geothermal, geological carbon sequestration, and oil and gas, typically involve maximizing either the extraction of energy or the storage of fluids. Fractures form the main pathways for flow in these systems, and locating these fractures is critical for predicting flow. However, fracture characterization is a highly uncertain process, and data from multiple sources, such as flow and geophysical are needed to reduce this uncertainty. We present a nonintrusive, sequential inversion framework for integrating data from geophysical and flow sources to constrain fracture networks in the subsurface. In this framework, we first estimate bounds on the statistics formore » the fracture orientations using microseismic data. These bounds are estimated through a combination of a focal mechanism (physics-based approach) and clustering analysis (statistical approach) of seismic data. Then, the fracture lengths are constrained using flow data. In conclusion, the efficacy of this inversion is demonstrated through a representative example.« less

Sequential geophysical and flow inversion to characterize fracture networks in subsurface systems

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

Mudunuru, Maruti Kumar; Karra, Satish; Makedonska, Nataliia

Subsurface applications, including geothermal, geological carbon sequestration, and oil and gas, typically involve maximizing either the extraction of energy or the storage of fluids. Fractures form the main pathways for flow in these systems, and locating these fractures is critical for predicting flow. However, fracture characterization is a highly uncertain process, and data from multiple sources, such as flow and geophysical are needed to reduce this uncertainty. We present a nonintrusive, sequential inversion framework for integrating data from geophysical and flow sources to constrain fracture networks in the subsurface. In this framework, we first estimate bounds on the statistics formore » the fracture orientations using microseismic data. These bounds are estimated through a combination of a focal mechanism (physics-based approach) and clustering analysis (statistical approach) of seismic data. Then, the fracture lengths are constrained using flow data. In conclusion, the efficacy of this inversion is demonstrated through a representative example.« less

3DFATMIC: THREE DIMENSIONAL SUBSURFACE FLOW, FATE AND TRANSPORT OF MICROBES AND CHEMICALS MODEL - USER'S MANUAL VERSION 1.0

EPA Science Inventory

This document is the user's manual of 3DFATMIC, a 3-Dimensional Subsurface Flow, Fate and Transport of Microbes and Chemicals Model using a Lagrangian-Eulerian adapted zooming and peak capturing (LEZOOMPC) algorithm.

Shallow subsurface temperatures and some estimates of heat flow from the Colorado Plateau of northeastern Arizona

USGS Publications Warehouse

Sass, J.H.; Stone, Claudia; Bills, D.J.

1982-01-01

Temperature data to depths of a few hundred meters were obtained from 29 wells in northeastern Arizona; 12 in the region surrounding the San Francisco Volcanic Field, 8 in the Black Mesa area, and 9 in the south-central Colorado Plateau which includes the White Mountains. Although there was evidence for local hydrologic disturbances in many temperature profiles, most wells provided an estimate of the conductive thermal gradient at the site. A few thermal conductivities were measured and were combined with published regional averages for the north-central part of the Colorado Plateau to produce crude estimates of regional heat flux. None of the wells was accessible below the regional aquifers. To these depths, heat flow in the area of the San Francisco Volcanic Field appears to be controlled primarily by regional lateral water movement having a significant downward vertical component of velocity. The mean heat flow of 27?5 mW^-2 is only a third to a quarter of what we would expect in this tectonic setting. The heat that is being carried laterally and downward probably is being discharged at low enthalpy and low elevation in springs and streams of the Colorado Plateau and Mogollon Rim. In the vicinity of Black Mesa, heat-flow averages about 60 mW^-2, characteristic of the 'cool interior' of the Colorado Plateau. North of the White Mountain Volcanic Field, the average heat flow is about 95 mW^-2.

Effects of Surface and Subsurface Bed Material Composition on Gravel Transport and Flow Competence Relations—Possibilities for Prediction

NASA Astrophysics Data System (ADS)

Bunte, K.; Abt, S. R.; Swingle, K. W.; Cenderelli, D. A.; Gaeuman, D. A.

2014-12-01

Bedload transport and flow competence relations are difficult to predict in coarse-bedded steep streams where widely differing sediment supply, bed stability, and complex flow hydraulics greatly affect amounts and sizes of transported gravel particles. This study explains how properties of bed material surface and subsurface size distributions are directly related to gravel transport and may be used for prediction of gravel transport and flow competence relations. Gravel transport, flow competence, and bed material size were measured in step-pool and plane-bed streams. Power functions were fitted to gravel transport QB=aQb and flow competence Dmax=cQd relations; Q is water discharge. Frequency distributions of surface FDsurf and subsurface FDsub bed material were likewise described by power functions FDsurf=hD j and FDsub=kDm fitted over six 0.5-phi size classes within 4 to 22.4 mm. Those gravel sizes are typically mobile even in moderate floods. Study results show that steeper subsurface bed material size distributions lead to steeper gravel transport and flow competence relations, whereas larger amounts of sediment contained in those 6 size bedmaterial classes (larger h and k) flatten the relations. Similarly, steeper surface size distributions decrease the coefficients of the gravel transport and flow competence relations, whereas larger amounts of sediment within the six bed material classes increase the intercepts of gravel transport and flow competence relations. Those relations are likely causative in streams where bedload stems almost entirely from the channel bed as opposed to direct (unworked) contributions from hillslopes and tributaries. The exponent of the subsurface bed material distribution m predicted the gravel transport exponent b with r2 near 0.7 and flow competence exponent d with r2 near 0.5. The intercept of bed surface distributions h increased the intercept a of gravel transport and c of the flow competence relations with r2 near 0.6.

Significance of Bioturbated Layer (BTL) and Deep Groundwater Storage on Runoff in Steep Saprolitic Tropical Lowlands Catchment

NASA Astrophysics Data System (ADS)

Cheng, Y.; Ogden, F. L.; Zhu, J.

2016-12-01

Bioturbated soil layers (BTLs) play a significant role in hydrological response and provisioning of ecosystem services in steep, saprolitic, tropical lowlands catchments. In this study, a new physically-based model formulation was developed for testing of runoff generation hypotheses. A main feature in the model formulation is explicit simulation of hydrological processes in the BTL including macropores, which our field observations show are ubiquitous, and deep groundwater stores that provide streamflow during the dry season The numerical model developed includes two main flow paths in the BTL, including one-dimensional (1D) vertical infiltration and two-dimensional (2D) lateral flows in both macropores and the soil matrix. Hydrological processes incorporated along with the BTL processes include intercepted rainfall, evapotranspiration, 2D surface flow and 1D deep groundwater discharge. This model was first tested in a 6.5 ha secondary succession catchment, that is under study by the Smithsonian Tropical Research Institute, Agua Salud project in Panama, which is dominated by steep slopes. With the incorporation of lateral macropore flow mechanism in the BTL, the model performs better than only including soil matrix flow in the BTL especially in simulating baseflow dynamics, which illustrates the importance of preferential flow from the BTL to stream discharge dynamics. The increase in the BTL thickness promotes more flow through the BTL and increases storage in both the BTL and the deep groundwater reservoir, but decreases the total streamflow and overland flow. Lateral macropore diameter distribution influences flows more than the macropore number or distribution type. The model has thus far passed falsification tests during the early wet season. Complexity in subsurface storage and base flow generation offer a new challenge for this model. The overall objective is to develop a model formulation that is useful in practical applications related to land-use management, provisioning of ecosystem services, and water security in similar tropical settings with distinct dry and wet seasons or in the humid tropics during periods of drought.

Glaciation and regional ground-water flow in the Fennoscandian Shield: Site 94

USGS Publications Warehouse

Provost, Alden M.; Voss, Clifford I.; Neuzil, C.E.

1998-01-01

Results from a regional-scale ground-water flow model of the Fennoscandian shield suggest that ground-water flow is strongly affected by surface conditions associated with climatic change and glaciation. The model was used to run a series of numerical simulations of variable-density ground-water flow in a 1500-km-long and approximately 10-km-deep cross-section that passes through southern Sweden. Ground-water flow and shield brine transport in the cross-sectional model are controlled by an assumed time evolution of surface conditions over the next 140 ka. Simulations show that, under periglacial conditions, permafrost may locally or extensively impede the free recharge or discharge of ground water. Below cold-based glacial ice, no recharge or discharge of ground water occurs. Both of these conditions result in the settling of shield brine and consequent freshening of near-surface water in areas of natural discharge blocked by permafrost. The presence of warm-based ice with basal melting creates a potential for ground-water recharge rates much larger than under present, ice-free conditions. Recharging basal meltwater can reach depths of a few kilometers in a few thousand years. The vast majority of recharged water is accommodated through storage in the volume of bedrock below the local area of recharge; regional (lateral) redistribution of recharged water by subsurface flow is minor over the duration of a glacial advance (~10 ka). During glacial retreat, the weight of the ice overlying a given surface location decreases, and significant upward flow of ground water may occur below the ice sheet due to pressure release, despite the continued potential for recharge of basal meltwater. Excess meltwater must exit from below the glacier through subglacial cavities and channels. Subsurface penetration of meltwater during glacial advance and up-flow during glacial retreat are greatest if the loading efficiency of the shield rock is low. The maximum rate of ground-water discharge occurs at the receding ice margin, and some discharge occurs below incursive post-glacial seas. The simulation results suggest that vertical movement of deep shield brines induced by the next few glacial cycles should not increase the concentration of dissolved solids significantly above present-day levels. However, the concentration of dissolved solids should decrease significantly at depths of up to several kilometers during periods of glacial meltwater recharge. The meltwater may reside in the subsurface for periods exceeding 10 ka and may bring oxygenated conditions to an otherwise reducing chemical environment.

Subsurface Flow and Moisture Dynamics in Response to Swash Motions: Effects of Beach Hydraulic Conductivity and Capillarity

NASA Astrophysics Data System (ADS)

Geng, Xiaolong; Heiss, James W.; Michael, Holly A.; Boufadel, Michel C.

2017-12-01

A combined field and numerical study was conducted to investigate dynamics of subsurface flow and moisture response to waves in the swash zone of a sandy beach located on Cape Henlopen, DE. A density-dependent variably saturated flow model MARUN was used to simulate subsurface flow beneath the swash zone. Values of hydraulic conductivity (K) and characteristic pore size (α, a capillary fringe property) were varied to evaluate their effects on subsurface flow and moisture dynamics in response to swash motions in beach aquifers. The site-specific modeling results were validated against spatiotemporal measurements of moisture and pore pressure in the beach. Sensitivity analyses indicated that the hydraulic conductivity and capillary fringe thickness of the beach greatly influenced groundwater flow pathways and associated transit times in the swash zone. A higher value of K enhanced swash-induced seawater infiltration into the beach, thereby resulting in a faster expansion of a wedge of high moisture content induced by swash cycles, and a flatter water table mound beneath the swash zone. In contrast, a thicker capillary fringe retained higher moisture content near the beach surface, and thus, significantly reduced the available pore space for infiltration of seawater. This attenuated wave effects on pore water flow in the unsaturated zone of the beach. Also, a thicker capillary fringe enhanced horizontal flow driven by the larger-scale hydraulic gradient caused by tides.

Integrated geoelectrical survey for groundwater and shallow subsurface evaluation: case study at Siliyin spring, El-Fayoum, Egypt

NASA Astrophysics Data System (ADS)

Metwaly, Mohamed; El-Qady, Gad; Massoud, Usama; El-Kenawy, Abeer; Matsushima, Jun; Al-Arifi, Nasser

2010-09-01

Siliyin spring is one of the many natural fresh water springs in the Western Desert of Egypt. It is located at the central part of El-Fayoum Delta, which is a potential place for urban developments and touristic activities. Integrated geoelectrical survey was conducted to facilitate mapping the groundwater resources and the shallow subsurface structures in the area. Twenty-eight transient electromagnetic (TEM) soundings, three vertical electrical soundings (VES) and three electrical resistivity tomography (ERT) profiles were carried out around the Siliyin spring location. The dense cultivation, the rugged topography and the existence of infra structure in the area hindered acquiring more data. The TEM data were inverted jointly with the VES and ERT, and constrained by available geological information. Based on the inversion results, a set of geoelectrical cross-sections have been constructed. The shallow sand to sandy clay layer that forms the shallow aquifer has been completely mapped underneath and around the spring area. Flowing of water from the Siliyin spring is interconnected with the lateral lithological changes from clay to sand soil. Exploration of the extension of Siliyin spring zone is recommended. The interpretation emphasizes the importance of integrating the geoelectrical survey with the available geological information to obtain useful, cheap and fast lithological and structural subsurface information.

A physically-based Distributed Hydrologic Model for Tropical Catchments

NASA Astrophysics Data System (ADS)

Abebe, N. A.; Ogden, F. L.

2010-12-01

Hydrological models are mathematical formulations intended to represent observed hydrological processes in a watershed. Simulated watersheds in turn vary in their nature based on their geographic location, altitude, climatic variables and geology and soil formation. Due to these variations, available hydrologic models vary in process formulation, spatial and temporal resolution and data demand. Many tropical watersheds are characterized by extensive and persistent biological activity and a large amount of rain. The Agua Salud catchments located within the Panama Canal Watershed, Panama, are such catchments identified by steep rolling topography, deep soils derived from weathered bedrock, and limited exposed bedrock. Tropical soils are highly affected by soil cracks, decayed tree roots and earthworm burrows forming a network of preferential flow paths that drain to a perched water table, which forms at a depth where the vertical hydraulic conductivity is significantly reduced near the bottom of the bioturbation layer. We have developed a physics-based, spatially distributed, multi-layered hydrologic model to simulate the dominant processes in these tropical watersheds. The model incorporates the major flow processes including overland flow, channel flow, matrix and non-Richards film flow infiltration, lateral downslope saturated matrix and non-Darcian pipe flow in the bioturbation layer, and deep saturated groundwater flow. Emphasis is given to the modeling of subsurface unsaturated zone soil moisture dynamics and the saturated preferential lateral flow from the network of macrospores. Preliminary results indicate that the model has the capability to simulate the complex hydrological processes in the catchment and will be a useful tool in the ongoing comprehensive ecohydrological studies in tropical catchments, and help improve our understanding of the hydrological effects of deforestation and aforestation.

Sub-crop geologic map of pre-Tertiary rocks in the Yucca Flat and northern Frenchman Flat areas, Nevada Test Site, southern Nevada

USGS Publications Warehouse

Cole, James C.; Harris, Anita G.; Wahl, Ronald R.

1997-01-01

This map displays interpreted structural and stratigraphic relations among the Paleozoic and older rocks of the Nevada Test Site region beneath the Miocene volcanic rocks and younger alluvium in the Yucca Flat and northern Frenchman Flat basins. These interpretations are based on a comprehensive examination and review of data for more than 77 drillholes that penetrated part of the pre-Tertiary basement beneath these post-middle Miocene structural basins. Biostratigraphic data from conodont fossils were newly obtained for 31 of these holes, and a thorough review of all prior microfossil paleontologic data is incorporated in the analysis. Subsurface relationships are interpreted in light of a revised regional geologic framework synthesized from detailed geologic mapping in the ranges surrounding Yucca Flat, from comprehensive stratigraphic studies in the region, and from additional detailed field studies on and around the Nevada Test Site.All available data indicate the subsurface geology of Yucca Flat is considerably more complicated than previous interpretations have suggested. The western part of the basin, in particular, is underlain by relics of the eastward-vergent Belted Range thrust system that are folded back toward the west and thrust by local, west-vergent contractional structures of the CP thrust system. Field evidence from the ranges surrounding the north end of Yucca Flat indicate that two significant strike-slip faults track southward beneath the post-middle Miocene basin fill, but their subsurface traces cannot be closely defined from the available evidence. In contrast, the eastern part of the Yucca Flat basin is interpreted to be underlain by a fairly simple north-trending, broad syncline in the pre-Tertiary units. Far fewer data are available for the northern Frenchman Flat basin, but regional analysis indicates the pre- Tertiary structure there should also be relatively simple and not affected by thrusting.This new interpretation has implications for ground water flow through pre-Tertiary rocks beneath the Yucca Flat and northern Frenchman Flat areas, and has consequences for ground water modeling and model validation. Our data indicate that the Mississippian Chainman Shale is not a laterally extensive confining unit in the western part of the basin because it is folded back onto itself by the convergent structures of the Belted Range and CP thrust systems. Early and Middle Paleozoic limestone and dolomite are present beneath most of both basins and, regardless of structural complications, are interpreted to form a laterally continuous and extensive carbonate aquifer. Structural culmination that marks the French Peak accommodation zone along the topographic divide between the two basins provides a lateral pathway through highly fractured rock between the volcanic aquifers of Yucca Flat and the regional carbonate aquifer. This pathway may accelerate the migration of ground-water contaminants introduced by underground nuclear testing toward discharge areas beyond the Nevada Test Site boundaries. Predictive three-dimensional models of hydrostratigraphic units and ground-water flow in the pre-Tertiary rocks of subsurface Yucca Flat are likely to be unrealistic due to the extreme structural complexities. The interpretation of hydrologic and geochemical data obtained from monitoring wells will be difficult to extrapolate through the flow system until more is known about the continuity of hydrostratigraphic units.

Gradual caldera collapse at Bárdarbunga volcano, Iceland, regulated by lateral magma outflow.

PubMed

Gudmundsson, Magnús T; Jónsdóttir, Kristín; Hooper, Andrew; Holohan, Eoghan P; Halldórsson, Sæmundur A; Ófeigsson, Benedikt G; Cesca, Simone; Vogfjörd, Kristín S; Sigmundsson, Freysteinn; Högnadóttir, Thórdís; Einarsson, Páll; Sigmarsson, Olgeir; Jarosch, Alexander H; Jónasson, Kristján; Magnússon, Eyjólfur; Hreinsdóttir, Sigrún; Bagnardi, Marco; Parks, Michelle M; Hjörleifsdóttir, Vala; Pálsson, Finnur; Walter, Thomas R; Schöpfer, Martin P J; Heimann, Sebastian; Reynolds, Hannah I; Dumont, Stéphanie; Bali, Eniko; Gudfinnsson, Gudmundur H; Dahm, Torsten; Roberts, Matthew J; Hensch, Martin; Belart, Joaquín M C; Spaans, Karsten; Jakobsson, Sigurdur; Gudmundsson, Gunnar B; Fridriksdóttir, Hildur M; Drouin, Vincent; Dürig, Tobias; Aðalgeirsdóttir, Guðfinna; Riishuus, Morten S; Pedersen, Gro B M; van Boeckel, Tayo; Oddsson, Björn; Pfeffer, Melissa A; Barsotti, Sara; Bergsson, Baldur; Donovan, Amy; Burton, Mike R; Aiuppa, Alessandro

2016-07-15

Large volcanic eruptions on Earth commonly occur with a collapse of the roof of a crustal magma reservoir, forming a caldera. Only a few such collapses occur per century, and the lack of detailed observations has obscured insight into the mechanical interplay between collapse and eruption. We use multiparameter geophysical and geochemical data to show that the 110-square-kilometer and 65-meter-deep collapse of Bárdarbunga caldera in 2014-2015 was initiated through withdrawal of magma, and lateral migration through a 48-kilometers-long dike, from a 12-kilometers deep reservoir. Interaction between the pressure exerted by the subsiding reservoir roof and the physical properties of the subsurface flow path explain the gradual, near-exponential decline of both collapse rate and the intensity of the 180-day-long eruption. Copyright © 2016, American Association for the Advancement of Science.

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

Hammond, Glenn Edward; Bao, J; Huang, M

Hyporheic exchange is a critical mechanism shaping hydrological and biogeochemical processes along a river corridor. Recent studies on quantifying the hyporheic exchange were mostly limited to local scales due to field inaccessibility, computational demand, and complexity of geomorphology and subsurface geology. Surface flow conditions and subsurface physical properties are well known factors on modulating the hyporheic exchange, but quantitative understanding of their impacts on the strength and direction of hyporheic exchanges at reach scales is absent. In this study, a high resolution computational fluid dynamics (CFD) model that couples surface and subsurface flow and transport is employed to simulate hyporheicmore » exchanges in a 7-km long reach along the main-stem of the Columbia River. Assuming that the hyporheic exchange does not affect surface water flow conditions due to its negligible magnitude compared to the volume and velocity of river water, we developed a one-way coupled surface and subsurface water flow model using the commercial CFD software STAR-CCM+. The model integrates the Reynolds-averaged Navier-Stokes (RANS) equation solver with a realizable κ-ε two-layer turbulence model, a two-layer all y + wall treatment, and the volume of fluid (VOF) method, and is used to simulate hyporheic exchanges by tracking the free water-air interface as well as flow in the river and the subsurface porous media. The model is validated against measurements from acoustic Doppler current profiler (ADCP) in the stream water and hyporheic fluxes derived from a set of temperature profilers installed across the riverbed. The validated model is then employed to systematically investigate how hyporheic exchanges are influenced by surface water fluid dynamics strongly regulated by upstream dam operations, as well as subsurface structures (e.g. thickness of riverbed and subsurface formation layers) and hydrogeological properties (e.g. permeability). The results suggest that the thickness of riverbed alluvium layer is the dominant factor for reach-scale hyporheic exchanges, followed by the alluvium permeability, the depth of the underlying impermeable layer, and the assumption of hydrostatic pressure.« less

Reconnecting tile drainage to riparian buffer hydrology for enhanced nitrate removal.

PubMed

Jaynes, D B; Isenhart, T M

2014-03-01

Riparian buffers are a proven practice for removing NO from overland flow and shallow groundwater. However, in landscapes with artificial subsurface (tile) drainage, most of the subsurface flow leaving fields is passed through the buffers in drainage pipes, leaving little opportunity for NO removal. We investigated the feasibility of re-routing a fraction of field tile drainage as subsurface flow through a riparian buffer for increasing NO removal. We intercepted an existing field tile outlet draining a 10.1-ha area of a row-cropped field in central Iowa and re-routed a fraction of the discharge as subsurface flow along 335 m of an existing riparian buffer. Tile drainage from the field was infiltrated through a perforated pipe installed 75 cm below the surface by maintaining a constant head in the pipe at a control box installed in-line with the existing field outlet. During 2 yr, >18,000 m (55%) of the total flow from the tile outlet was redirected as infiltration within the riparian buffer. The redirected water seeped through the 60-m-wide buffer, raising the water table approximately 35 cm. The redirected tile flow contained 228 kg of NO. On the basis of the strong decrease in NO concentrations within the shallow groundwater across the buffer, we hypothesize that the NO did not enter the stream but was removed within the buffer by plant uptake, microbial immobilization, or denitrification. Redirecting tile drainage as subsurface flow through a riparian buffer increased its NO removal benefit and is a promising management practice to improve surface water quality within tile-drained landscapes. Copyright © by the American Society of Agronomy, Crop Science Society of America, and Soil Science Society of America, Inc.

Hydrology and subsurface transport of oil-field brine at the U.S. Geological Survey OSPER site "A", Osage County, Oklahoma

USGS Publications Warehouse

Herkelrath, W.N.; Kharaka, Y.K.; Thordsen, J.J.; Abbott, M.M.

2007-01-01

Spillage and improper disposal of saline produced water from oil wells has caused environmental damage at thousands of sites in the United States. In order to improve understanding of the fate and transport of contaminants at these sites, the U.S. Geological Survey carried out multidisciplinary investigations at two oil production sites near Skiatook Lake, Oklahoma. As a part of this effort, the hydrology and subsurface transport of brine at OSPER site "A", a tank battery and pit complex that was abandoned in 1973, was investigated. Based on data from 41 new boreholes that were cored and completed with monitoring wells, a large (???200 m ?? 200 m ?? 20 m) plume of saline ground water was mapped. The main dissolved species are Na and Cl, with TDS in the plume ranging as high as 30,000 mg/L. Analysis of the high barometric efficiency of the wells indicated a confined aquifer response. Well-slug tests indicated the hydraulic conductivity is low (0.3-7.0 cm/day). Simplified flow and transport modeling supports the following conceptual model: (1) prior to the produced water releases, recharge was generally low (???1 cm/a); (2) in ???60 a of oil production enough saline produced water in pits leaked into the subsurface to create the plume; (3) following abandonment of the site in 1973 and filling of Skiatook Reservoir in the mid-1980s, recharge and lateral flow of water through the plume returned to low values; (4) as a result, spreading of the brine plume caused by mixing with fresh ground water recharge, as well as natural attenuation, are very slow.

Spatial and temporal variation of residence time and storage volume of subsurface water evaluated by multi-tracers approach in mountainous headwater catchments

NASA Astrophysics Data System (ADS)

Tsujimura, Maki; Yano, Shinjiro; Abe, Yutaka; Matsumoto, Takehiro; Yoshizawa, Ayumi; Watanabe, Ysuhito; Ikeda, Koichi

2015-04-01

Headwater catchments in mountainous region are the most important recharge area for surface and subsurface waters, additionally time and stock information of the water is principal to understand hydrological processes in the catchments. However, there have been few researches to evaluate variation of residence time and storage volume of subsurface water in time and space at the mountainous headwaters especially with steep slope. We performed an investigation on age dating and estimation of storage volume using simple water budget model in subsurface water with tracing of hydrological flow processes in mountainous catchments underlain by granite, Paleozoic and Tertiary, Yamanashi and Tsukuba, central Japan. We conducted hydrometric measurements and sampling of spring, stream and ground waters in high-flow and low-flow seasons from 2008 through 2012 in the catchments, and CFCs, stable isotopic ratios of oxygen-18 and deuterium, inorganic solute constituent concentrations were determined on all water samples. Residence time of subsurface water ranged from 11 to 60 years in the granite catchments, from 17 to 32 years in the Paleozoic catchments, from 13 to 26 years in the Tertiary catchments, and showed a younger age during the high-flow season, whereas it showed an older age in the low-flow season. Storage volume of subsurface water was estimated to be ranging from 10 ^ 4 to 10 ^ 6 m3 in the granite catchments, from 10 ^ 5 to 10 ^ 7 m3 in the Paleozoic catchments, from 10 ^ 4 to 10 ^ 6 m3 in the Tertiary catchments. In addition, seasonal change of storage volume in the granite catchments was the highest as compared with those of the Paleozoic and the Tertiary catchments. The results suggest that dynamic change of hydrological process seems to cause a larger variation of the residence time and storage volume of subsurface water in time and space in the granite catchments, whereas higher groundwater recharge rate due to frequent fissures or cracks seems to cause larger storage volume of the subsurface water in the Paleozoic catchments though the variation is not so considerable. Also, numerical simulation results support these findings.

Baseline hydraulic performance of the Heathrow constructed wetlands subsurface flow system.

PubMed

Richter, K M; Margetts, J R; Saul, A J; Guymer, I; Worrall, P

2003-01-01

A constructed wetland treatment system has been commissioned by BAA (formerly the British Airports Authority) in order to attenuate airfield runoff contaminated with de-icant and other potentially polluting materials from Heathrow Airport. Airfield runoff containing de-icants has the potential to impose significant oxygen demands on water bodies. The site consists of a number of integrated treatment systems, including a 1 ha rafted reed bed canal system and a 2 ha sub-surface flow gravel reed bed. This research project is concerned with the performance of the subsurface flow reed beds, though attention will be paid in this paper to the operation of the whole system. Prior to the planting of the subsurface flow reed beds, flow-tracing experiments were carried out on the three different types of subsurface flow beds, so that the baseline performance of the system could be quantified. In association, data regarding the soil organic matter content was also collected prior to the planting of the beds. As expected, soil organic matter content is observed to be negligible within the bed, though a small amount of build up was observed in localised areas on the surface of the beds. This was attributed to the growth of algae in depressions where standing water persisted during the construction phase. Few studies exist which provide detailed measurements into the cause and effect of variations in hydraulic conductivity within an operational reed bed system. The data presented here form the baseline results for an ongoing study into the investigation of the change in hydraulic conductivity of an operational reed bed system.

Evidence for preferential flow through sandstone aquifers in Southern Wisconsin

USGS Publications Warehouse

Swanson, S.K.; Bahr, J.M.; Bradbury, K.R.; Anderson, K.M.

2006-01-01

Sandstones often escape extensive hydrogeologic characterization due to their high primary porosity and perceived homogeneity of permeability. This study provides evidence for laterally extensive, high permeability zones in the Tunnel City Group, an undeformed, Cambrian-aged sandstone unit that exists in the subsurface throughout much of central and southern Wisconsin, USA. Several discrete high-permeability zones were identified in boreholes using flow logging and slug tests, and the interconnectedness of the features was tested using a site-specific numerical model for springs in the region. Explicit incorporation of a high-permeability layer leads to improvements in the flux calibration over simulations that lack the features, thus supporting the hydraulic continuity of high-permeability zones in the sandstone aquifer over tens of kilometers. The results suggest that stratigraphically controlled heterogeneities like contrasts in lithology or bedding-plane fractures, which have been shown to strongly influence the flow of groundwater in more heterogeneous sedimentary rocks, may also deserve close examination in sandstone aquifers. ?? 2005 Elsevier B.V. All rights reserved.

THE HYDROLOGIC SYSTEM: GEOMORPHIC AND HYDROGEOLOGIC CONTROLS ON SURFACE AND SUBSURFACE FLOW REGIMES IN RIPARIAN MEADOW ECOSYSTEMS IN THE CENTRAL GREAT BASIN

EPA Science Inventory

Riparian corridors in upland watersheds in the Great Basin of central Nevada contain the majority of the region's biodiversity. Water, in both surface and subsurface flow regimes, is an important resource sustaining these sensitive ecosystems and other similar riparian ecosystem...

Subsurface temperatures and surface heat flow in the Michigan Basin and their relationships to regional subsurface fluid movement

USGS Publications Warehouse

Vugrinovich, R.

1989-01-01

Linear regression of 405 bottomhole temperature (BHT) measurements vs. associated depths from Michigan's Lower Peninsula results in the following equation relating BHT and depth: BHT(??C) = 14.5 + 0.0192 ?? depth(m) Temperature residuals, defined as (BHT measured)-(BHT calculated), were determined for each of the 405 BHT's. Areas of positive temperature residuals correspond to areas of regional groundwater discharge (determined from maps of equipotential surface) while areas of negative temperature residuals correspond to areas of regional groundwater recharge. These relationships are observed in the principal aquifers in rocks of Devonian and Ordovician age and in a portion of the principal aquifer in rocks of Silurian age. There is a similar correspondence between high surface heat flow (determined using the silica geothermometer) and regional groundwater discharge areas and low surface heat flow and regional groundwater recharge areas. Post-Jurassic depositional and tectonic histories suggest that the observed coupling of subsurface temperature and groundwater flow systems may have persisted since Jurassic time. Thus the higher subsurface palaeotemperatures (and palaeogeothermal gradients) indicated by recent studies most likely pre-date the Jurassic. ?? 1989.

Conceptual Model Evaluation using Advanced Parameter Estimation Techniques with Heat as a Tracer

NASA Astrophysics Data System (ADS)

Naranjo, R. C.; Morway, E. D.; Healy, R. W.

2016-12-01

Temperature measurements made at multiple depths beneath the sediment-water interface has proven useful for estimating seepage rates from surface-water channels and corresponding subsurface flow direction. Commonly, parsimonious zonal representations of the subsurface structure are defined a priori by interpretation of temperature envelopes, slug tests or analysis of soil cores. However, combining multiple observations into a single zone may limit the inverse model solution and does not take full advantage of the information content within the measured data. Further, simulating the correct thermal gradient, flow paths, and transient behavior of solutes may be biased by inadequacies in the spatial description of subsurface hydraulic properties. The use of pilot points in PEST offers a more sophisticated approach to estimate the structure of subsurface heterogeneity. This presentation evaluates seepage estimation in a cross-sectional model of a trapezoidal canal with intermittent flow representing four typical sedimentary environments. The recent improvements in heat as a tracer measurement techniques (i.e. multi-depth temperature probe) along with use of modern calibration techniques (i.e., pilot points) provides opportunities for improved calibration of flow models, and, subsequently, improved model predictions.

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

Noble gas and hydrocarbon tracers in multiphase unconventional hydrocarbon systems: Toward integrated advanced reservoir simulators

NASA Astrophysics Data System (ADS)

Darrah, T.; Moortgat, J.; Poreda, R. J.; Muehlenbachs, K.; Whyte, C. J.

2015-12-01

Although hydrocarbon production from unconventional energy resources has increased dramatically in the last decade, total unconventional oil and gas recovery from black shales is still less than 25% and 9% of the totals in place, respectively. Further, the majority of increased hydrocarbon production results from increasing the lengths of laterals, the number of hydraulic fracturing stages, and the volume of consumptive water usage. These strategies all reduce the economic efficiency of hydrocarbon extraction. The poor recovery statistics result from an insufficient understanding of some of the key physical processes in complex, organic-rich, low porosity formations (e.g., phase behavior, fluid-rock interactions, and flow mechanisms at nano-scale confinement and the role of natural fractures and faults as conduits for flow). Noble gases and other hydrocarbon tracers are capably of recording subsurface fluid-rock interactions on a variety of geological scales (micro-, meso-, to macro-scale) and provide analogs for the movement of hydrocarbons in the subsurface. As such geochemical data enrich the input for the numerical modeling of multi-phase (e.g., oil, gas, and brine) fluid flow in highly heterogeneous, low permeability formations Herein we will present a combination of noble gas (He, Ne, Ar, Kr, and Xe abundances and isotope ratios) and molecular and isotopic hydrocarbon data from a geographically and geologically diverse set of unconventional hydrocarbon reservoirs in North America. Specifically, we will include data from the Marcellus, Utica, Barnett, Eagle Ford, formations and the Illinois basin. Our presentation will include geochemical and geological interpretation and our perspective on the first steps toward building an advanced reservoir simulator for tracer transport in multicomponent multiphase compositional flow (presented separately, in Moortgat et al., 2015).

Hydrogeologic setting and ground water flow beneath a section of Indian River Bay, Delaware

USGS Publications Warehouse

Krantz, David E.; Manheim, Frank T.; Bratton, John F.; Phelan, Daniel J.

2004-01-01

The small bays along the Atlantic coast of the Delmarva Peninsula (Delaware, Maryland, and Virginia) are a valuable natural resource, and an asset for commerce and recreation. These coastal bays also are vulnerable to eutrophication from the input of excess nutrients derived from agriculture and other human activities in the watersheds. Ground water discharge may be an appreciable source of fresh water and a transport pathway for nutrients entering the bays. This paper presents results from an investigation of the physical properties of the surficial aquifer and the processes associated with ground water flow beneath Indian River Bay, Delaware. A key aspect of the project was the deployment of a new technology, streaming horizontal resistivity, to map the subsurface distribution of fresh and saline ground water beneath the bay. The resistivity profiles showed complex patterns of ground water flow, modes of mixing, and submarine ground water discharge. Cores, gamma and electromagnetic-induction logs, and in situ ground water samples collected during a coring operation in Indian River Bay verified the interpretation of the resistivity profiles. The shore-parallel resistivity lines show subsurface zones of fresh ground water alternating with zones dominated by the flow of salt water from the estuary down into the aquifer. Advective flow produces plumes of fresh ground water 400 to 600 m wide and 20 m thick that may extend more than 1 km beneath the estuary. Zones of dispersive mixing between fresh and saline ground water develop on the upper, lower, and lateral boundaries of the the plume. the plumes generally underlie small incised valleys that can be traced landward to stream draining the upland. The incised valleys are filled with 1 to 2 m of silt and peat that act as a semiconfining layer to restrict the downward flow of salt water from the estuary. Active circulation of both the fresh and saline ground water masses beneath the bay is inferred from the geophysical results and supported by geochemical data.

Inputs and losses by surface runoff and subsurface leaching for pastures managed by continuous or rotational stocking.

PubMed

Owens, L B; Barker, D J; Loerch, S C; Shipitalo, M J; Bonta, J V; Sulc, R M

2012-01-01

Pasture management practices can affect forage quality and production, animal health and production, and surface and groundwater quality. In a 5-yr study conducted at the North Appalachian Experimental Watershed near Coshocton, Ohio, we compared the effects of two contrasting grazing methods on surface and subsurface water quantity and quality. Four pastures, each including a small, instrumented watershed (0.51-1.09 ha) for surface runoff measurements and a developed spring for subsurface flow collection, received 112 kg N ha(-1) yr(-1) and were grazed at similar stocking rates (1.8-1.9 cows ha(-1)). Two pastures were continuously stocked; two were subdivided so that they were grazed with frequent rotational stocking (5-6 times weekly). In the preceding 5 yr, these pastures received 112 kg N ha(-1) yr(-1) after several years of 0 N fertilizer and were grazed with weekly rotational stocking. Surface runoff losses of N were minimal. During these two periods, some years had precipitation up to 50% greater than the long-term average, which increased subsurface flow and NO(3)-N transport. Average annual NO(3)-N transported in subsurface flow from the four watersheds during the two 5-yr periods ranged from 11.3 to 22.7 kg N ha(-1), which was similar to or less than the mineral-N received in precipitation. Flow and transport variations were greater among seasons than among watersheds. Flow-weighted seasonal NO(3)-N concentrations in subsurface flow did not exceed 7 mg L(-1). Variations in NO(3)-N leached from pastures were primarily due to variable precipitation rather than the effects of continuous, weekly rotational, or frequent rotational stocking practices. This suggests that there was no difference among these grazing practices in terms of NO(3)-N leaching. Copyright © by the American Society of Agronomy, Crop Science Society of America, and Soil Science Society of America, Inc.

Salinization Sources Along the Lower Jordan River Under Draught Conditions

NASA Astrophysics Data System (ADS)

Holtzman, R.; Shavit, U.; Segal, M.; Vengosh, A.; Farber, E.; Gavrieli, I.

2003-12-01

The Lower Jordan River, once a flowing freshwater river, is suffering from an ongoing reduction of discharge and water quality. The river flows between the Sea of Galilee and the Dead Sea, an aerial distance of about 105 Km. The severe reduction is caused by an excessive exploitation of its sources and diversion of sewage and agricultural drainage into the river. The extreme low flows and low water quality threaten the natural existence of the river and its potential use for agriculture. In spite of its importance, little research has been done in the river. The objectives of the study were to measure the discharge and water composition along the river and to evaluate the main sources that control its flow and chemical characteristics. The hypothesis of the study was that interaction with subsurface flows significantly affects the river flow and chemical composition. The research is based on a detailed field study, which included flow rate measurements in the river and its tributaries, water sampling and analysis and mass balance calculations of water and solutes. A portable Acoustic Doppler Velocimeter (ADV) was used to measure velocities and bathymetry at different locations across the river sections. Due to accessibility constraints, a floating traverse construction, which enables the ADV's deployment from one bank of the river, was developed. It was found that flow rate ranges between 500-1,100 L/s in northern (upstream) sections and 300-1,650 L/s in the south. This low discharge represents a significant reduction from historical values and is lower than recent published estimations. This research represents base flows only, as the measurements were done during a period of two consecutive draught years. Calculated mass balance of water flows in the northern sections shows that the subsurface source contributes to the river around 200-670 L/s (30-80% of the river flow). Calculations of solute balance show that the subsurface flows add 20-50% of the mass of solutes (e.g. Sulfate) that flows in the river. The assumption of a hydraulic gradient that points at inflows from subsurface flows is encouraged by high water levels measured in nearby piezometers. Possible natural subsurface sources include shallow groundwater or rising of water from deep formations. The existence of adjacent thermal wells strengthens the reasonability of such water rise. Possible anthropogenic sources include return flows and effluents. The results are consistent and agree with the geochemical and isotopic analyses. It is concluded that the impact of the subsurface component on the Jordan River is significant and must be taken into consideration, for future water management schemes and implementation of the Peace Treaty between Israel and Jordan.

Heat flow and subsurface temperature as evidence for basin-scale ground-water flow, North Slope of Alaska

USGS Publications Warehouse

Deming, D.; Sass, J.H.; Lachenbruch, A.H.; De Rito, R. F.

1992-01-01

Several high-resolution temperature logs were made in each of 21 drillholes and a total of 601 thermal conductivity measurements were made on drill cuttings and cores. Near-surface heat flow (??20%) is inversely correlated with elevation and ranges from a low of 27 mW/m2 in the foothills of the Brooks Range in the south, to a high of 90 mW/m2 near the north coast. Subsurface temperatures and thermal gradients estimated from corrected BHTs are similarly much higher on the coastal plain than in the foothills province to the south. Significant east-west variation in heat flow and subsurface temperature is also observed; higher heat flow and temperature coincide with higher basement topography. The observed thermal pattern is consistent with forced convection by a topographically driven ground-water flow system. Average ground-water (Darcy) velocity in the postulated flow system is estimated to be of the order of 0.1 m/yr; the effective basin-scale permeability is estimated to be of the order of 10-14 m2. -from Authors

Towards a high resolution, integrated hydrology model of North America: Diagnosis of feedbacks between groundwater and land energy fluxes at continental scales.

NASA Astrophysics Data System (ADS)

Maxwell, Reed; Condon, Laura

2016-04-01

Recent studies demonstrate feedbacks between groundwater dynamics, overland flow, land surface and vegetation processes, and atmospheric boundary layer development that significantly affect local and regional climate across a range of climatic conditions. Furthermore, the type and distribution of vegetation cover alters land-atmosphere water and energy fluxes, as well as runoff generation and overland flow processes. These interactions can result in significant feedbacks on local and regional climate. In mountainous regions, recent research has shown that spatial and temporal variability in annual evapotranspiration, and thus water budgets, is strongly dependent on lateral groundwater flow; however, the full effects of these feedbacks across varied terrain (e.g. from plains to mountains) are not well understood. Here, we present a high-resolution, integrated hydrology model that covers much of continental North America and encompasses the Mississippi and Colorado watersheds. The model is run in a fully-transient manner at hourly temporal resolution incorporating fully-coupled land energy states and fluxes with integrated surface and subsurface hydrology. Connections are seen between hydrologic variables (such as water table depth) and land energy fluxes (such as latent heat) and spatial and temporal scaling is shown to span many orders of magnitude. Model results suggest that partitioning of plant transpiration to bare soil evaporation is a function of water table depth and later groundwater flow. Using these transient simulations as a proof of concept, we present a vision for future integrated simulation capabilities.

Application of a pore-scale reactive transport model to a natural analog for reaction-induced pore alterations

DOE PAGES

Yoon, Hongkyu; Major, Jonathan; Dewers, Thomas; ...

2017-01-05

Dissolved CO 2 in the subsurface resulting from geological CO 2 storage may react with minerals in fractured rocks, confined aquifers, or faults, resulting in mineral precipitation and dissolution. The overall rate of reaction can be affected by coupled processes including hydrodynamics, transport, and reactions at the (sub) pore-scale. In this work pore-scale modeling of coupled fluid flow, reactive transport, and heterogeneous reactions at the mineral surface is applied to account for permeability alterations caused by precipitation-induced pore-blocking. This paper is motivated by observations of CO 2 seeps from a natural CO 2 sequestration analog, Crystal Geyser, Utah. Observations alongmore » the surface exposure of the Little Grand Wash fault indicate the lateral migration of CO 2 seep sites (i.e., alteration zones) of 10–50 m width with spacing on the order of ~100 m over time. Sandstone permeability in alteration zones is reduced by 3–4 orders of magnitude by carbonate cementation compared to unaltered zones. One granular porous medium and one fracture network systems are used to conceptually represent permeable porous media and locations of conduits controlled by fault-segment intersections and/or topography, respectively. Simulation cases accounted for a range of reaction regimes characterized by the Damköhler (Da) and Peclet (Pe) numbers. Pore-scale simulation results demonstrate that combinations of transport (Pe), geochemical conditions (Da), solution chemistry, and pore and fracture configurations contributed to match key patterns observed in the field of how calcite precipitation alters flow paths by pore plugging. This comparison of simulation results with field observations reveals mechanistic explanations of the lateral migration and enhances our understanding of subsurface processes associated with the CO 2 injection. In addition, permeability and porosity relations are constructed from pore-scale simulations which account for a range of reaction regimes characterized by the Da and Pe numbers. Finally, the functional relationships obtained from pore-scale simulations can be used in a continuum scale model that may account for large-scale phenomena mimicking lateral migration of surface CO 2 seeps.« less

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

Yoon, Hongkyu; Major, Jonathan; Dewers, Thomas

Dissolved CO 2 in the subsurface resulting from geological CO 2 storage may react with minerals in fractured rocks, confined aquifers, or faults, resulting in mineral precipitation and dissolution. The overall rate of reaction can be affected by coupled processes including hydrodynamics, transport, and reactions at the (sub) pore-scale. In this work pore-scale modeling of coupled fluid flow, reactive transport, and heterogeneous reactions at the mineral surface is applied to account for permeability alterations caused by precipitation-induced pore-blocking. This paper is motivated by observations of CO 2 seeps from a natural CO 2 sequestration analog, Crystal Geyser, Utah. Observations alongmore » the surface exposure of the Little Grand Wash fault indicate the lateral migration of CO 2 seep sites (i.e., alteration zones) of 10–50 m width with spacing on the order of ~100 m over time. Sandstone permeability in alteration zones is reduced by 3–4 orders of magnitude by carbonate cementation compared to unaltered zones. One granular porous medium and one fracture network systems are used to conceptually represent permeable porous media and locations of conduits controlled by fault-segment intersections and/or topography, respectively. Simulation cases accounted for a range of reaction regimes characterized by the Damköhler (Da) and Peclet (Pe) numbers. Pore-scale simulation results demonstrate that combinations of transport (Pe), geochemical conditions (Da), solution chemistry, and pore and fracture configurations contributed to match key patterns observed in the field of how calcite precipitation alters flow paths by pore plugging. This comparison of simulation results with field observations reveals mechanistic explanations of the lateral migration and enhances our understanding of subsurface processes associated with the CO 2 injection. In addition, permeability and porosity relations are constructed from pore-scale simulations which account for a range of reaction regimes characterized by the Da and Pe numbers. Finally, the functional relationships obtained from pore-scale simulations can be used in a continuum scale model that may account for large-scale phenomena mimicking lateral migration of surface CO 2 seeps.« less

A one-dimensional model of subsurface hillslope flow

Treesearch

Jason C. Fisher

1997-01-01

Abstract - A one-dimensional, finite difference model of saturated subsurface flow within a hillslope was developed. The model uses rainfall, elevation data, a hydraulic conductivity, and a storage coefficient to predict the saturated thickness in time and space. The model was tested against piezometric data collected in a swale located in the headwaters of the North...

Ecohydrologic Investigations of Shallow Lateral Subsurface Flow in Tropical Soils using Time-Lapse Surface Electrical Resistivity Tomography

NASA Astrophysics Data System (ADS)

Ogden, F. L.; Mojica, A.; Abebe, N. A.; Smithsonian Tropical Research Institute, Panama Canal Watershed Experiment, Agua Salud Project

2010-12-01

The hydrologic effects of deforestation and aforestation in the tropics remain an area of active research. Hydrologic predictions of land-use change effects remain elusive. One of the unique features of catchment hydrology in the tropics is the effect of intense, continuous biological activity by insects, shrubs, trees, and small mammals. Sapprolitic soils derived from weathered bedrock cover widespread areas. These soils have low matrix permeabilities on the order of 1 mm/h, are 10 to 20 m in thickness and have relatively low activity because they have been depleted of light cations by annual rainfall over 2000 mm. As part of the Smithsonian Tropical Research Institute, Panama Canal Watershed Experiment, Agua Salud Project, we have observed shallow subsurface flow in tropical soils in central Panama using an introduced salinity contrast and surface electrical resistivity tomography (ERT). In 2009 and 2010, experiments were conducted in a 30 year-old secondary succession forest, and in two former pasture sites that were planted with native timber species and teak, respectively, in 2008. At each site, saline water (NaCl tagged with LiBr) was introduced to the soil using two different methods: soil pits and ponded surface applications. Results showed the strongest response in the case of ponded surface applications with observed changes in resistivity between -50% and 50%. In soil pit applications, the change in electrical resistivity varied from -10% to 10%. Results suggest that in the case of surface application, a transient perched water table is created near the bottom of the bioturbation layer that activates the downslope macropore network and results in bulk flow velocities that are significantly higher than observed soil matrix permeabilities. When heavy rainfall occurred during tests, increased mobility of the salinity contrast more clearly showed the active layer where most flow occurred. Time-series ERT observations enabled measurements of downslope bulk flow velocities over 1 m/h, presumably due to the existing downslope macroporosity network. These observations are being used to estimate macroporosity network properties and constrain hydrologic model parameters in different land uses. These results show that these non-invasive tests are a useful tool to determine the distribution of downslope lateral flow generated from pit and surface-applied saline solutions. ERT experimental results from a hillslope-scale experiment in central Panama, showing change in electrical conductivity from 30-minutes to 330-minutes after continuous injection of salinity contrast at x=0.

Influence of vertical and lateral heat transfer on permafrost thaw, peatland landscape transition, and groundwater flow

USGS Publications Warehouse

Kurylyk, Barret L.; Masaki, Masaki; Quinton, William L.; McKenzie, Jeffrey M.; Voss, Clifford I.

2016-01-01

Recent climate change has reduced the spatial extent and thickness of permafrost in many discontinuous permafrost regions. Rapid permafrost thaw is producing distinct landscape changes in the Taiga Plains of the Northwest Territories, Canada. As permafrost bodies underlying forested peat plateaus shrink, the landscape slowly transitions into unforested wetlands. The expansion of wetlands has enhanced the hydrologic connectivity of many watersheds via new surface and near-surface flow paths, and increased streamflow has been observed. Furthermore, the decrease in forested peat plateaus results in a net loss of boreal forest and associated ecosystems. This study investigates fundamental processes that contribute to permafrost thaw by comparing observed and simulated thaw development and landscape transition of a peat plateau-wetland complex in the Northwest Territories, Canada from 1970 to 2012. Measured climate data are first used to drive surface energy balance simulations for the wetland and peat plateau. Near-surface soil temperatures simulated in the surface energy balance model are then applied as the upper boundary condition to a three-dimensional model of subsurface water flow and coupled energy transport with freeze-thaw. Simulation results demonstrate that lateral heat transfer, which is not considered in many permafrost models, can influence permafrost thaw rates. Furthermore, the simulations indicate that landscape evolution arising from permafrost thaw acts as a positive feedback mechanism that increases the energy absorbed at the land surface and produces additional permafrost thaw. The modeling results also demonstrate that flow rates in local groundwater flow systems may be enhanced by the degradation of isolated permafrost bodies.

Evaluation of SCS-CN method using a fully distributed physically based coupled surface-subsurface flow model

NASA Astrophysics Data System (ADS)

Shokri, Ali

2017-04-01

The hydrological cycle contains a wide range of linked surface and subsurface flow processes. In spite of natural connections between surface water and groundwater, historically, these processes have been studied separately. The current trend in hydrological distributed physically based model development is to combine distributed surface water models with distributed subsurface flow models. This combination results in a better estimation of the temporal and spatial variability of the interaction between surface and subsurface flow. On the other hand, simple lumped models such as the Soil Conservation Service Curve Number (SCS-CN) are still quite common because of their simplicity. In spite of the popularity of the SCS-CN method, there have always been concerns about the ambiguity of the SCS-CN method in explaining physical mechanism of rainfall-runoff processes. The aim of this study is to minimize these ambiguity by establishing a method to find an equivalence of the SCS-CN solution to the DrainFlow model, which is a fully distributed physically based coupled surface-subsurface flow model. In this paper, two hypothetical v-catchment tests are designed and the direct runoff from a storm event are calculated by both SCS-CN and DrainFlow models. To find a comparable solution to runoff prediction through the SCS-CN and DrainFlow, the variance between runoff predictions by the two models are minimized by changing Curve Number (CN) and initial abstraction (Ia) values. Results of this study have led to a set of lumped model parameters (CN and Ia) for each catchment that is comparable to a set of physically based parameters including hydraulic conductivity, Manning roughness coefficient, ground surface slope, and specific storage. Considering the lack of physical interpretation in CN and Ia is often argued as a weakness of SCS-CN method, the novel method in this paper gives a physical explanation to CN and Ia.

Coupling a three-dimensional subsurface flow and transport model with a land surface model to simulate stream-aquifer-land interactions (CP v1.0)

NASA Astrophysics Data System (ADS)

Bisht, Gautam; Huang, Maoyi; Zhou, Tian; Chen, Xingyuan; Dai, Heng; Hammond, Glenn E.; Riley, William J.; Downs, Janelle L.; Liu, Ying; Zachara, John M.

2017-12-01

A fully coupled three-dimensional surface and subsurface land model is developed and applied to a site along the Columbia River to simulate three-way interactions among river water, groundwater, and land surface processes. The model features the coupling of the Community Land Model version 4.5 (CLM4.5) and a massively parallel multiphysics reactive transport model (PFLOTRAN). The coupled model, named CP v1.0, is applied to a 400 m × 400 m study domain instrumented with groundwater monitoring wells along the Columbia River shoreline. CP v1.0 simulations are performed at three spatial resolutions (i.e., 2, 10, and 20 m) over a 5-year period to evaluate the impact of hydroclimatic conditions and spatial resolution on simulated variables. Results show that the coupled model is capable of simulating groundwater-river-water interactions driven by river stage variability along managed river reaches, which are of global significance as a result of over 30 000 dams constructed worldwide during the past half-century. Our numerical experiments suggest that the land-surface energy partitioning is strongly modulated by groundwater-river-water interactions through expanding the periodically inundated fraction of the riparian zone, and enhancing moisture availability in the vadose zone via capillary rise in response to the river stage change. Meanwhile, CLM4.5 fails to capture the key hydrologic process (i.e., groundwater-river-water exchange) at the site, and consequently simulates drastically different water and energy budgets. Furthermore, spatial resolution is found to significantly impact the accuracy of estimated the mass exchange rates at the boundaries of the aquifer, and it becomes critical when surface and subsurface become more tightly coupled with groundwater table within 6 to 7 meters below the surface. Inclusion of lateral subsurface flow influenced both the surface energy budget and subsurface transport processes as a result of river-water intrusion into the subsurface in response to an elevated river stage that increased soil moisture for evapotranspiration and suppressed available energy for sensible heat in the warm season. The coupled model developed in this study can be used for improving mechanistic understanding of ecosystem functioning and biogeochemical cycling along river corridors under historical and future hydroclimatic changes. The dataset presented in this study can also serve as a good benchmarking case for testing other integrated models.

Coupling a three-dimensional subsurface flow and transport model with a land surface model to simulate stream–aquifer–land interactions (CP v1.0)

DOE PAGES

Bisht, Gautam; Huang, Maoyi; Zhou, Tian; ...

2017-12-12

A fully coupled three-dimensional surface and subsurface land model is developed and applied to a site along the Columbia River to simulate three-way interactions among river water, groundwater, and land surface processes. The model features the coupling of the Community Land Model version 4.5 (CLM4.5) and a massively parallel multiphysics reactive transport model (PFLOTRAN). The coupled model, named CP v1.0, is applied to a 400 m × 400 m study domain instrumented with groundwater monitoring wells along the Columbia River shoreline. CP v1.0 simulations are performed at three spatial resolutions (i.e., 2, 10, and 20 m) over a 5-year periodmore » to evaluate the impact of hydroclimatic conditions and spatial resolution on simulated variables. Results show that the coupled model is capable of simulating groundwater–river-water interactions driven by river stage variability along managed river reaches, which are of global significance as a result of over 30 000 dams constructed worldwide during the past half-century. Our numerical experiments suggest that the land-surface energy partitioning is strongly modulated by groundwater–river-water interactions through expanding the periodically inundated fraction of the riparian zone, and enhancing moisture availability in the vadose zone via capillary rise in response to the river stage change. Meanwhile, CLM4.5 fails to capture the key hydrologic process (i.e., groundwater–river-water exchange) at the site, and consequently simulates drastically different water and energy budgets. Furthermore, spatial resolution is found to significantly impact the accuracy of estimated the mass exchange rates at the boundaries of the aquifer, and it becomes critical when surface and subsurface become more tightly coupled with groundwater table within 6 to 7 meters below the surface. Inclusion of lateral subsurface flow influenced both the surface energy budget and subsurface transport processes as a result of river-water intrusion into the subsurface in response to an elevated river stage that increased soil moisture for evapotranspiration and suppressed available energy for sensible heat in the warm season. The coupled model developed in this study can be used for improving mechanistic understanding of ecosystem functioning and biogeochemical cycling along river corridors under historical and future hydroclimatic changes. The dataset presented in this study can also serve as a good benchmarking case for testing other integrated models.« less

Coupling a three-dimensional subsurface flow and transport model with a land surface model to simulate stream–aquifer–land interactions (CP v1.0)

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

Bisht, Gautam; Huang, Maoyi; Zhou, Tian

A fully coupled three-dimensional surface and subsurface land model is developed and applied to a site along the Columbia River to simulate three-way interactions among river water, groundwater, and land surface processes. The model features the coupling of the Community Land Model version 4.5 (CLM4.5) and a massively parallel multiphysics reactive transport model (PFLOTRAN). The coupled model, named CP v1.0, is applied to a 400 m × 400 m study domain instrumented with groundwater monitoring wells along the Columbia River shoreline. CP v1.0 simulations are performed at three spatial resolutions (i.e., 2, 10, and 20 m) over a 5-year period to evaluate themore » impact of hydroclimatic conditions and spatial resolution on simulated variables. Results show that the coupled model is capable of simulating groundwater–river-water interactions driven by river stage variability along managed river reaches, which are of global significance as a result of over 30 000 dams constructed worldwide during the past half-century. Our numerical experiments suggest that the land-surface energy partitioning is strongly modulated by groundwater–river-water interactions through expanding the periodically inundated fraction of the riparian zone, and enhancing moisture availability in the vadose zone via capillary rise in response to the river stage change. Meanwhile, CLM4.5 fails to capture the key hydrologic process (i.e., groundwater–river-water exchange) at the site, and consequently simulates drastically different water and energy budgets. Furthermore, spatial resolution is found to significantly impact the accuracy of estimated the mass exchange rates at the boundaries of the aquifer, and it becomes critical when surface and subsurface become more tightly coupled with groundwater table within 6 to 7 meters below the surface. Inclusion of lateral subsurface flow influenced both the surface energy budget and subsurface transport processes as a result of river-water intrusion into the subsurface in response to an elevated river stage that increased soil moisture for evapotranspiration and suppressed available energy for sensible heat in the warm season. The coupled model developed in this study can be used for improving mechanistic understanding of ecosystem functioning and biogeochemical cycling along river corridors under historical and future hydroclimatic changes. The dataset presented in this study can also serve as a good benchmarking case for testing other integrated models.« less

Coupling a three-dimensional subsurface flow and transport model with a land surface model to simulate stream–aquifer–land interactions (CP v1.0)

DOE PAGES

Bisht, Gautam; Huang, Maoyi; Zhou, Tian; ...

2017-01-01

A fully coupled three-dimensional surface and subsurface land model is developed and applied to a site along the Columbia River to simulate three-way interactions among river water, groundwater, and land surface processes. The model features the coupling of the Community Land Model version 4.5 (CLM4.5) and a massively parallel multiphysics reactive transport model (PFLOTRAN). The coupled model, named CP v1.0, is applied to a 400 m × 400 m study domain instrumented with groundwater monitoring wells along the Columbia River shoreline. CP v1.0 simulations are performed at three spatial resolutions (i.e., 2, 10, and 20 m) over a 5-year period to evaluate themore » impact of hydroclimatic conditions and spatial resolution on simulated variables. Results show that the coupled model is capable of simulating groundwater–river-water interactions driven by river stage variability along managed river reaches, which are of global significance as a result of over 30 000 dams constructed worldwide during the past half-century. Our numerical experiments suggest that the land-surface energy partitioning is strongly modulated by groundwater–river-water interactions through expanding the periodically inundated fraction of the riparian zone, and enhancing moisture availability in the vadose zone via capillary rise in response to the river stage change. Meanwhile, CLM4.5 fails to capture the key hydrologic process (i.e., groundwater–river-water exchange) at the site, and consequently simulates drastically different water and energy budgets. Furthermore, spatial resolution is found to significantly impact the accuracy of estimated the mass exchange rates at the boundaries of the aquifer, and it becomes critical when surface and subsurface become more tightly coupled with groundwater table within 6 to 7 meters below the surface. Inclusion of lateral subsurface flow influenced both the surface energy budget and subsurface transport processes as a result of river-water intrusion into the subsurface in response to an elevated river stage that increased soil moisture for evapotranspiration and suppressed available energy for sensible heat in the warm season. The coupled model developed in this study can be used for improving mechanistic understanding of ecosystem functioning and biogeochemical cycling along river corridors under historical and future hydroclimatic changes. The dataset presented in this study can also serve as a good benchmarking case for testing other integrated models.« less

Coupling a three-dimensional subsurface flow and transport model with a land surface model to simulate stream–aquifer–land interactions (CP v1.0)

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

Bisht, Gautam; Huang, Maoyi; Zhou, Tian

A fully coupled three-dimensional surface and subsurface land model is developed and applied to a site along the Columbia River to simulate three-way interactions among river water, groundwater, and land surface processes. The model features the coupling of the Community Land Model version 4.5 (CLM4.5) and a massively parallel multiphysics reactive transport model (PFLOTRAN). The coupled model, named CP v1.0, is applied to a 400 m × 400 m study domain instrumented with groundwater monitoring wells along the Columbia River shoreline. CP v1.0 simulations are performed at three spatial resolutions (i.e., 2, 10, and 20 m) over a 5-year periodmore » to evaluate the impact of hydroclimatic conditions and spatial resolution on simulated variables. Results show that the coupled model is capable of simulating groundwater–river-water interactions driven by river stage variability along managed river reaches, which are of global significance as a result of over 30 000 dams constructed worldwide during the past half-century. Our numerical experiments suggest that the land-surface energy partitioning is strongly modulated by groundwater–river-water interactions through expanding the periodically inundated fraction of the riparian zone, and enhancing moisture availability in the vadose zone via capillary rise in response to the river stage change. Meanwhile, CLM4.5 fails to capture the key hydrologic process (i.e., groundwater–river-water exchange) at the site, and consequently simulates drastically different water and energy budgets. Furthermore, spatial resolution is found to significantly impact the accuracy of estimated the mass exchange rates at the boundaries of the aquifer, and it becomes critical when surface and subsurface become more tightly coupled with groundwater table within 6 to 7 meters below the surface. Inclusion of lateral subsurface flow influenced both the surface energy budget and subsurface transport processes as a result of river-water intrusion into the subsurface in response to an elevated river stage that increased soil moisture for evapotranspiration and suppressed available energy for sensible heat in the warm season. The coupled model developed in this study can be used for improving mechanistic understanding of ecosystem functioning and biogeochemical cycling along river corridors under historical and future hydroclimatic changes. The dataset presented in this study can also serve as a good benchmarking case for testing other integrated models.« less

Subsurface Water Flow and its Subsequent Impact on Chemical Behavior

USDA-ARS?s Scientific Manuscript database

The impact of the subsurface stratigraphy on crop growth and agrichemical behavior has been studied for several years at the OPE3 research site located at the USDA-ARS Beltsville Agricultural Research Center, in Beltsville Maryland. This site contains subsurface restricting layers that have been id...

Improving Flash Flood Prediction in Multiple Environments

NASA Astrophysics Data System (ADS)

Broxton, P. D.; Troch, P. A.; Schaffner, M.; Unkrich, C.; Goodrich, D.; Wagener, T.; Yatheendradas, S.

2009-12-01

Flash flooding is a major concern in many fast responding headwater catchments . There are many efforts to model and to predict these flood events, though it is not currently possible to adequately predict the nature of flash flood events with a single model, and furthermore, many of these efforts do not even consider snow, which can, by itself, or in combination with rainfall events, cause destructive floods. The current research is aimed at broadening the applicability of flash flood modeling. Specifically, we will take a state of the art flash flood model that is designed to work with warm season precipitation in arid environments, the KINematic runoff and EROSion model (KINEROS2), and combine it with a continuous subsurface flow model and an energy balance snow model. This should improve its predictive capacity in humid environments where lateral subsurface flow significantly contributes to streamflow, and it will make possible the prediction of flooding events that involve rain-on-snow or rapid snowmelt. By modeling changes in the hydrologic state of a catchment before a flood begins, we can also better understand the factors or combination of factors that are necessary to produce large floods. Broadening the applicability of an already state of the art flash flood model, such as KINEROS2, is logical because flash floods can occur in all types of environments, and it may lead to better predictions, which are necessary to preserve life and property.

Role of rainwater induced subsurface flow in water-level dynamics and thermoerosion of shallow thermokarst ponds on the Northeastern Qinghai-Tibet Plateau

NASA Astrophysics Data System (ADS)

Pan, X.; Yu, Q.; You, Y.

2014-12-01

Understanding hydrological and thermal regimes of thermokarst lakes is of great importance for predicting their responses to climate change. However, mechanism of water-level dynamics and associated thermal effects on thermoerosion of thermokarst lakes are still not well understood on the Qinghai-Tibet Plateau (QTP). In this study, we investigate two typical shallow thermokarst ponds (namely small lakes) in a warm permafrost region with thick active layer on the northeastern QTP through quantifying water budget. Results demonstrate that, rainfall induced subsurface lateral flow dominates pond water-level regime. Annual variation of pond water-level relies on areal water budget of surrounding active layer, particularly the high variable of precipitation. Besides, it is worth noting the extraordinary warming during the late ice-cover period, because marked air gap between upper ice-cover and underlying water, led by the upward thawing of thick ice-cover, might result in greenhouse-like condition due to the unique weather that strong solar radiation and little snowpack. This hydrological mechanism also exerts evident impacts on thermal regime and thermoerosion of the shallow thermokarst ponds, and they are closely related to retreat of thermokarst pondshore and underlying permafrost degradation. These findings imply a localized model addressing the unique hydrological and thermal regimes of thermokarst lakes would be essential to study the evolution of these shallow rainwater dominated thermokarst ponds on the QTP.

A locally conservative stabilized continuous Galerkin finite element method for two-phase flow in poroelastic subsurfaces

NASA Astrophysics Data System (ADS)

Deng, Q.; Ginting, V.; McCaskill, B.; Torsu, P.

2017-10-01

We study the application of a stabilized continuous Galerkin finite element method (CGFEM) in the simulation of multiphase flow in poroelastic subsurfaces. The system involves a nonlinear coupling between the fluid pressure, subsurface's deformation, and the fluid phase saturation, and as such, we represent this coupling through an iterative procedure. Spatial discretization of the poroelastic system employs the standard linear finite element in combination with a numerical diffusion term to maintain stability of the algebraic system. Furthermore, direct calculation of the normal velocities from pressure and deformation does not entail a locally conservative field. To alleviate this drawback, we propose an element based post-processing technique through which local conservation can be established. The performance of the method is validated through several examples illustrating the convergence of the method, the effectivity of the stabilization term, and the ability to achieve locally conservative normal velocities. Finally, the efficacy of the method is demonstrated through simulations of realistic multiphase flow in poroelastic subsurfaces.

Finding the best windows: An apparent environmental threshold determines which diffuse flows are dominated by subsurface microbes

NASA Astrophysics Data System (ADS)

Olins, H. C.; Rogers, D.; Scholin, C. A.; Preston, C. J.; Vidoudez, C.; Ussler, W.; Pargett, D.; Jensen, S.; Roman, B.; Birch, J. M.; Girguis, P. R.

2014-12-01

Hydrothermal vents are hotspots of microbial primary productivity often described as "windows into the subsurface biosphere." High temperature vents have received the majority of research attention, but cooler diffuse flows are as, if not more, important a source of heat and chemicals to the overlying ocean. We studied patterns of in situ gene expression and co-registered geochemistry in order to 1) describe the diversity and physiological poise of active microbial communities that span thermal and geochemical gradients from active diffuse flow to background vent field seawater, and 2) determine to what extent seawater or subsurface microbes were active throughout this environment. Analyses of multiple metatranscriptomes from 5 geochemically distinct sites (some from samples preserved in situ) show that proximate diffuse flows showed strikingly different transcription profiles. Specifically, caldera background and some diffuse flows were similar, both dominated by seawater-derived Gammaproteobacteria despite having distinct geochemistries. Intra-field community shows evidence of increased primary productivity throughout the entire vent field and not just at individual diffuse flows. In contrast, a more spatially limited, Epsilonproteobacteria-dominated transcription profile from the most hydrothermally-influenced diffuse flow appeared to be driven by the activity of vent-endemic microbes, likely reflecting subsurface microbial activity. We suggest that the microbial activity within many diffuse flow vents is primarily attributable to seawater derived Gammaproteobacterial sulfur oxidizers, while in certain other flows vent-endemic Epsilonproteobactiera are most active. These data reveal a diversity in microbial activity at diffuse flows that has not previously been recognized, and reshapes our thinking about the relative influence that different microbial communities may have on local processes (such as primary production) and potentially global biogeochemical cycles.

[Correlation of substrate structure and hydraulic characteristics in subsurface flow constructed wetlands].

PubMed

Bai, Shao-Yuan; Song, Zhi-Xin; Ding, Yan-Li; You, Shao-Hong; He, Shan

2014-02-01

The correlation of substrate structure and hydraulic characteristics was studied by numerical simulation combined with experimental method. The numerical simulation results showed that the permeability coefficient of matrix had a great influence on hydraulic efficiency in subsurface flow constructed wetlands. The filler with a high permeability coefficient had a worse flow field distribution in the constructed wetland with single layer structure. The layered substrate structure with the filler permeability coefficient increased from surface to bottom could avoid the short-circuited flow and dead-zones, and thus, increased the hydraulic efficiency. Two parallel pilot-scale constructed wetlands were built according to the numerical simulation results, and tracer experiments were conducted to validate the simulation results. The tracer experiment result showed that hydraulic characteristics in the layered constructed wetland were obviously better than that in the single layer system, and the substrate effective utilization rates were 0.87 and 0.49, respectively. It was appeared that numerical simulation would be favorable for substrate structure optimization in subsurface flow constructed wetlands.

Profiling USGA putting greens using GPR

USDA-ARS?s Scientific Manuscript database

All USGA-specification putting greens require a subsurface drainage system. A typical subsurface installation is a herringbone pattern of buried 100-mm dia. PVC drainage pipes, designed such that the central main line is placed along the line of maximum slope. Laterals are spaced no more than 5 m, r...

Rapid runoff via shallow throughflow and deeper preferential flow in a boreal catchment underlain by frozen silt (Alaska, USA)

USGS Publications Warehouse

Koch, Joshua C.; Ewing, Stephanie A.; Striegl, Robert G.; McKnight, Diane M.

2013-01-01

In high-latitude catchments where permafrost is present, runoff dynamics are complicated by seasonal active-layer thaw, which may cause a change in the dominant flowpaths as water increasingly contacts mineral soils of low hydraulic conductivity. A 2-year study, conducted in an upland catchment in Alaska (USA) underlain by frozen, well-sorted eolian silt, examined changes in infiltration and runoff with thaw. It was hypothesized that rapid runoff would be maintained by flow through shallow soils during the early summer and deeper preferential flow later in the summer. Seasonal changes in soil moisture, infiltration, and runoff magnitude, location, and chemistry suggest that transport is rapid, even when soils are thawed to their maximum extent. Between June and September, a shift occurred in the location of runoff, consistent with subsurface preferential flow in steep and wet areas. Uranium isotopes suggest that late summer runoff erodes permafrost, indicating that substantial rapid flow may occur along the frozen boundary. Together, throughflow and deep preferential flow may limit upland boreal catchment water and solute storage, and subsequently biogeochemical cycling on seasonal to annual timescales. Deep preferential flow may be important for stream incision, network drainage development, and the release of ancient carbon to ecosystems

Channel water balance and exchange with subsurface flow along a mountain headwater stream in Montana, United States

USGS Publications Warehouse

Payn, R.A.; Gooseff, M.N.; McGlynn, B.L.; Bencala, K.E.; Wondzell, S.M.

2009-01-01

Channel water balances of contiguous reaches along streams represent a poorly understood scale of stream-subsurface interaction. We measured reach water balances along a headwater stream in Montana, United States, during summer base flow recessions. Reach water balances were estimated from series of tracer tests in 13 consecutive reaches delineated evenly along a 2.6 km valley segment. For each reach, we estimated net change in discharge, gross hydrologic loss, and gross hydrologic gain from tracer dilution and mass recovery. Four series of tracer tests were performed during relatively high, intermediate, and low base flow conditions. The relative distribution of channel water along the stream was strongly related to a transition in valley structure, with a general increase in gross losses through the recession. During tracer tests at intermediate and low flows, there were frequent substantial losses of tracer mass (>10%) that could not be explained by net loss in flow over the reach, indicating that many of the study reaches were concurrently losing and gaining water. For example, one reach with little net change in discharge exchanged nearly 20% of upstream flow with gains and losses along the reach. These substantial bidirectional exchanges suggest that some channel interactions with subsurface flow paths were not measurable by net change in flow or transient storage of recovered tracer. Understanding bidirectional channel water balances in stream reaches along valleys is critical to an accurate assessment of stream solute fate and transport and to a full assessment of exchanges between the stream channel and surrounding subsurface.

Channel water balance and exchange with subsurface flow along a mountain headwater stream in Montana, United States

USGS Publications Warehouse

Payn, R.A.; Gooseff, M.N.; McGlynn, B.L.; Bencala, K.E.; Wondzell, S.M.

2009-01-01

Channel water balances of contiguous reaches along streams represent a poorly understood scale of stream-subsurface interaction. We measured reach water balances along a headwater stream in Montana, United States, during summer base flow recessions. Reach water balances were estimated from series of tracer tests in 13 consecutive reaches delineated evenly along a 2.6 km valley segment. For each reach, we estimated net change in discharge, gross hydrologic loss, and gross hydrologic gain from tracer dilution and mass recovery. Four series of tracer tests were performed during relatively high, intermediate, and low base flow conditions. The relative distribution of channel water along the stream was strongly related to a transition in valley structure, with a general increase in gross losses through the recession. During tracer tests at intermediate and low flows, there were frequent substantial losses of tracer mass (>10%) that could not be explained by net loss in flow over the reach, indicating that many of the study reaches were concurrently losing and gaining water. For example, one reach with little net change in discharge exchanged nearly 20% of upstream flow with gains and losses along the reach. These substantial bidirectional exchanges suggest that some channel interactions with subsurface flow paths were not measurable by net change in flow or transient storage of recovered tracer. Understanding bidirectional channel water balances in stream reaches along valleys is critical to an accurate assessment of stream solute fate and transport and to a full assessment of exchanges between the stream channel and surrounding subsurface. Copyright 2009 by the American Geophysical Union.

Enhanced Amendment Delivery to Subsurface Using Shear Thinning Fluid and Aqueous Foam for Metal, Radionuclide, and NAPL Remediation

NASA Astrophysics Data System (ADS)

Zhong, L.; Szecsody, J.; Li, X.; Oostrom, M.; Truex, M.

2010-12-01

In many contamination sites, removal of contaminants by any active remediation efforts is not practical due to the high cost and technological limitations. Alternatively, in situ remediation is expected to be the most important remediation strategy. Delivery of reactive amendment to the contamination zone is essential for the reactions between the contaminants and remedial amendments to proceed in situ. It is a challenge to effectively deliver remedial amendment to the subsurface contamination source areas in both aquifer and vadose zone. In aquifer, heterogeneity induces fluid bypassing the low-permeability zones, resulting in certain contaminated areas inaccessible to the remedial amendment delivered by water injection, thus inhibiting the success of remedial operations. In vadose zone in situ remediation, conventional solution injection and infiltration for amendment delivery have difficulties to achieve successful lateral spreading and uniform distribution of the reactive media. These approaches also tend to displace highly mobile metal and radionuclide contaminants such as hexavalent chromium [Cr(VI)] and technetium (Tc-99), causing spreading of contaminations. Shear thinning fluid and aqueous foam can be applied to enhance the amendment delivery and improve in situ subsurface remediation efficiency under aquifer and vadose zone conditions, respectively. Column and 2-D flow cell experiments were conducted to demonstrate the enhanced delivery and improved remediation achieved by the application of shear thinning fluid and foam injection at the laboratory scale. Solutions of biopolymer xanthan gum were used as the shear thinning delivering fluids. Surfactant sodium lauryl ether sulfate (STEOL CS-330) was the foaming agent. The shear thinning fluid delivery (STFD) considerably improved the sweeping efficiency over a heterogeneous system and enhanced the non-aqueous liquid phase (NAPL) removal. The delivery of amendment into low-perm zones (LPZs) by STFD also increased the persistence of amendment solution in the LPZs after injection. Immobilization of Tc-99 was improved when a reductant was delivered by foam versus by water-based solution to contaminated vadose zone sediments. Foam delivery remarkably improved the lateral distribution of fluids compared to direct liquid injection. In heterogeneous vadose zone formation, foam injection increased the liquid flow in the high permeable zones into which very limited fluid was distributed during liquid infiltration, demonstrating improved amendment distribution uniformity in the heterogeneous system by foam delivery.

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.

Analytical solution for tension-saturated and unsaturated flow from wicking porous pipes in subsurface irrigation: The Kornev-Philip legacies revisited

NASA Astrophysics Data System (ADS)

Kacimov, A. R.; Obnosov, Yu. V.

2017-03-01

The Russian engineer Kornev in his 1935 book raised perspectives of subsurface "negative pressure" irrigation, which have been overlooked in modern soil science. Kornev's autoirrigation utilizes wicking of a vacuumed water from a porous pipe into a dry adjacent soil. We link Kornev's technology with a slightly modified Philip (1984)'s analytical solutions for unsaturated flow from a 2-D cylindrical pipe in an infinite domain. Two Darcian flows are considered and connected through continuity of pressure along the pipe-soil contact. The first fragment is a thin porous pipe wall in which water seeps at tension saturation; the hydraulic head is a harmonic function varying purely radially across the wall. The Thiem solution in this fragment gives the boundary condition for azimuthally varying suction pressure in the second fragment, ambient soil, making the exterior of the pipe. The constant head, rather than Philip's isobaricity boundary condition, along the external wall slightly modifies Philip's formulae for the Kirchhoff potential and pressure head in the soil fragment. Flow characteristics (magnitudes of the Darcian velocity, total flow rate, and flow net) are explicitly expressed through series of Macdonald's functions. For a given pipe's external diameter, wall thickness, position of the pipe above a free water datum in the supply tank, saturated conductivities of the wall and soil, and soil's sorptive number, a nonlinear equation with respect to the total discharge from the pipe is obtained and solved by a computer algebra routine. Efficiency of irrigation is evaluated by computation of the moisture content within selected zones surrounding the porous pipe.