Numerical Modeling of Coupled Water Flow and Heat Transport in Soil and Snow

NASA Astrophysics Data System (ADS)

Kelleners, T.

2015-12-01

A numerical model is developed to calculate coupled water flow and heat transport in seasonally frozen soil and snow. Both liquid water flow and water vapor flow are included. The effect of dissolved ions on soil water freezing point depression is included by combining an expression for osmotic head with the Clapeyron equation and the van Genuchten soil water retention function. The coupled water flow and heat transport equations are solved using the Thomas algorithm and Picard iteration. Ice pressure is always assumed zero and frost heave is neglected. The new model is tested using data from a high-elevation rangeland soil that is subject to significant soil freezing and a mountainous forest soil that is snow-covered for about 8 months of the year. Soil hydraulic parameters are mostly based on measurements and only vegetation parameters are fine-tuned to match measured and calculated soil water content, soil & snow temperature, and snow height. Modeling statistics for both systems show good performance for temperature, intermediate performance for snow height, and relatively low performance for soil water content, in accordance with earlier results with an older version of the model.

Modeling Surface Water Flow in the Atchafalaya Basin

NASA Astrophysics Data System (ADS)

Liu, K.; Simard, M.

2017-12-01

While most of the Mississippi River Delta is sinking due to insufficient sediment supply and subsidence, the stable wetlands and the prograding delta systems in the Atchafalaya Basin provide a unique opportunity to study the constructive interactions between riverine and marine forcings and their impacts upon coastal morphology. To better understand the hydrodynamics in this region, we developed a numerical modeling system for the water flow through the river channel - deltas - wetlands networks in the Atchafalaya Basin. Determining spatially varying model parameters for a large area composed of such diverse land cover types poses a challenge to developing an accurate numerical model. For example, the bottom friction coefficient can not be measured directly and the available elevation maps for the wetlands in the basin are inaccurate. To overcome these obstacles, we developed the modeling system in three steps. Firstly, we modeled river bathymetry based on in situ sonar transects and developed a simplified 1D model for the Wax Lake Outlet using HEC-RAS. Secondly, we used a Bayesian approach to calibrate the model automatically and infer important unknown parameters such as riverbank elevation and bottom friction coefficient through Markov Chain Monte Carlo (MCMC) simulations. We also estimated the wetland elevation based on the distribution of different vegetation species in the basin. Thirdly, with the lessons learnt from the 1D model, we developed a depth-averaged 2D model for the whole Atchafalaya Basin using Delft3D. After calibrations, the model successfully reproduced the water levels measured at five gauges in the Wax Lake Outlet and the modeled water surface profile along the channel agreed reasonably well with our LIDAR measurements. In addition, the model predicted a one-hour delay in tidal phase from the Wax Lake Delta to the upstream gauge. In summary, this project presents a procedure to initialize hydrology model parameters that integrates field

Modeling the Death Valley regional ground-water flow system

USGS Publications Warehouse

D'Agnese, F. A.; Faunt, C.C.; Hill, M.C.

2004-01-01

The development of a regional ground-water flow model of the Death Valley region in the southwestern United States is discussed in the context of the fourteen guidelines of Hill. This application of the guidelines demonstrates how they may be used for model calibration and evaluation, and to direct further model development and data collection.

Mathematical modelling of surface water-groundwater flow and salinity interactions in the coastal zone

NASA Astrophysics Data System (ADS)

Spanoudaki, Katerina; Kampanis, Nikolaos A.

2014-05-01

Coastal areas are the most densely-populated areas in the world. Consequently water demand is high, posing great pressure on fresh water resources. Climatic change and its direct impacts on meteorological variables (e.g. precipitation) and indirect impact on sea level rise, as well as anthropogenic pressures (e.g. groundwater abstraction), are strong drivers causing groundwater salinisation and subsequently affecting coastal wetlands salinity with adverse effects on the corresponding ecosystems. Coastal zones are a difficult hydrologic environment to represent with a mathematical model due to the large number of contributing hydrologic processes and variable-density flow conditions. Simulation of sea level rise and tidal effects on aquifer salinisation and accurate prediction of interactions between coastal waters, groundwater and neighbouring wetlands requires the use of integrated surface water-groundwater models. In the past few decades several computer codes have been developed to simulate coupled surface and groundwater flow. In these numerical models surface water flow is usually described by the 1-D Saint Venant equations (e.g. Swain and Wexler, 1996) or the 2D shallow water equations (e.g. Liang et al., 2007). Further simplified equations, such as the diffusion and kinematic wave approximations to the Saint Venant equations, are also employed for the description of 2D overland flow and 1D stream flow (e.g. Gunduz and Aral, 2005). However, for coastal bays, estuaries and wetlands it is often desirable to solve the 3D shallow water equations to simulate surface water flow. This is the case e.g. for wind-driven flows or density-stratified flows. Furthermore, most integrated models are based on the assumption of constant fluid density and therefore their applicability to coastal regions is questionable. Thus, most of the existing codes are not well-suited to represent surface water-groundwater interactions in coastal areas. To this end, the 3D integrated

Assessing the urban water balance: the Urban Water Flow Model and its application in Cyprus.

PubMed

Charalambous, Katerina; Bruggeman, Adriana; Lange, Manfred A

2012-01-01

Modelling the urban water balance enables the understanding of the interactions of water within an urban area and allows for better management of water resources. However, few models today provide a comprehensive overview of all water sources and uses. The objective of the current paper was to develop a user-friendly tool that quantifies and visualizes all water flows, losses and inefficiencies in urban environments. The Urban Water Flow Model was implemented in a spreadsheet and includes a water-savings application that computes the contributions of user-selected saving options to the overall water balance. The model was applied to the coastal town of Limassol, Cyprus, for the hydrologic years 2003/04-2008/09. Data were collected from the different authorities and hydrologic equations and estimations were added to complete the balance. Average precipitation was 363 mm/yr, amounting to 25.4 × 10(6)m(3)/yr, more than double the annual potable water supply to the town. Surface runoff constituted 29.6% of all outflows, while evapotranspiration from impervious areas was 21.6%. Possible potable water savings for 2008/09 were estimated at 5.3 × 10(3) m(3), which is 50% of the total potable water provided to the area. This saving would also result in a 6% reduction of surface runoff.

Using a Content Management System for Integrated Water Quantity, Quality and Instream Flows Modeling

NASA Astrophysics Data System (ADS)

Burgholzer, R.; Brogan, C. O.; Scott, D.; Keys, T.

2017-12-01

With increased population and water demand, in-stream flows can become depleted by consumptive uses and dilution of permitted discharges may be compromised. Reduced flows downstream of water withdrawals may increase the violation rate of bacterial concentrations from direct deposition by livestock and wildlife. Water storage reservoirs are constructed and operated to insure more stable supplies for consumptive demands and dilution flows, however their use comes at the cost of increased evaporative losses, potential for thermal pollution, interrupted fish migration, and reduced flooding events that are critical to maintain habitat and water quality. Due to this complex interrelationship between water quantity, quality and instream habitat comprehensive multi-disciplinary models must be developed to insure long-term sustainability of water resources and to avoid conflicts between drinking water, food and energy production, and aquatic biota. The Commonwealth of Virginia funded the expansion of the Chesapeake Bay Program Phase 5 model to cover the entire state, and has been using this model to evaluate water supply permit and planning since 2009. This integrated modeling system combines a content management system (Drupal and PHP) for model input data and leverages the modularity of HSPF with the custom segmentation and parameterization routines programmed by modelers working with the Chesapeake Bay Program. The model has been applied to over 30 Virginia Water Permits, instream flows and aquatic habitat models and a Virginias 30 year water supply demand projections. Future versions will leverage the Bay Model auto-calibration routines for adding small-scale water supply and TMDL models, utilize climate change scenarios, and integrate Virginia's reservoir management modules into the Chesapeake Bay watershed model, feeding projected demand and operational changes back up to EPA models to improve the realism of future Bay-wide simulations.

Effects of linking a soil-water-balance model with a groundwater-flow model

USGS Publications Warehouse

Stanton, Jennifer S.; Ryter, Derek W.; Peterson, Steven M.

2013-01-01

A previously published regional groundwater-flow model in north-central Nebraska was sequentially linked with the recently developed soil-water-balance (SWB) model to analyze effects to groundwater-flow model parameters and calibration results. The linked models provided a more detailed spatial and temporal distribution of simulated recharge based on hydrologic processes, improvement of simulated groundwater-level changes and base flows at specific sites in agricultural areas, and a physically based assessment of the relative magnitude of recharge for grassland, nonirrigated cropland, and irrigated cropland areas. Root-mean-squared (RMS) differences between the simulated and estimated or measured target values for the previously published model and linked models were relatively similar and did not improve for all types of calibration targets. However, without any adjustment to the SWB-generated recharge, the RMS difference between simulated and estimated base-flow target values for the groundwater-flow model was slightly smaller than for the previously published model, possibly indicating that the volume of recharge simulated by the SWB code was closer to actual hydrogeologic conditions than the previously published model provided. Groundwater-level and base-flow hydrographs showed that temporal patterns of simulated groundwater levels and base flows were more accurate for the linked models than for the previously published model at several sites, particularly in agricultural areas.

ENVIRONMENTAL RESEARCH BRIEF : ANALYTIC ELEMENT MODELING OF GROUND-WATER FLOW AND HIGH PERFORMANCE COMPUTING

EPA Science Inventory

Several advances in the analytic element method have been made to enhance its performance and facilitate three-dimensional ground-water flow modeling in a regional aquifer setting. First, a new public domain modular code (ModAEM) has been developed for modeling ground-water flow ...

Adaptive hydrological flow field modeling based on water body extraction and surface information

NASA Astrophysics Data System (ADS)

Puttinaovarat, Supattra; Horkaew, Paramate; Khaimook, Kanit; Polnigongit, Weerapong

2015-01-01

Hydrological flow characteristic is one of the prime indicators for assessing flood. It plays a major part in determining drainage capability of the affected basin and also in the subsequent simulation and rainfall-runoff prediction. Thus far, flow directions were typically derived from terrain data which for flat landscapes are obscured by other man-made structures, hence undermining the practical potential. In the absence (or diminutive) of terrain slopes, water passages have a more pronounced effect on flow directions than elevations. This paper, therefore, presents detailed analyses and implementation of hydrological flow modeling from satellite and topographic images. Herein, gradual assignment based on support vector machine was applied to modified normalized difference water index and a digital surface model, in order to ensure reliable water labeling while suppressing modality-inherited artifacts and noise. Gradient vector flow was subsequently employed to reconstruct the flow field. Experiments comparing the proposed scheme with conventional water boundary delineation and flow reconstruction were presented. Respective assessments revealed its advantage over the generic stream burning. Specifically, it could extract water body from studied areas with 98.70% precision, 99.83% recall, 98.76% accuracy, and 99.26% F-measure. The correlations between resultant flows and those obtained from the stream burning were as high as 0.80±0.04 (p≤0.01 in all resolutions).

TopoDrive and ParticleFlow--Two Computer Models for Simulation and Visualization of Ground-Water Flow and Transport of Fluid Particles in Two Dimensions

USGS Publications Warehouse

Hsieh, Paul A.

2001-01-01

This report serves as a user?s guide for two computer models: TopoDrive and ParticleFlow. These two-dimensional models are designed to simulate two ground-water processes: topography-driven flow and advective transport of fluid particles. To simulate topography-driven flow, the user may specify the shape of the water table, which bounds the top of the vertical flow section. To simulate transport of fluid particles, the model domain is a rectangle with overall flow from left to right. In both cases, the flow is under steady state, and the distribution of hydraulic conductivity may be specified by the user. The models compute hydraulic head, ground-water flow paths, and the movement of fluid particles. An interactive visual interface enables the user to easily and quickly explore model behavior, and thereby better understand ground-water flow processes. In this regard, TopoDrive and ParticleFlow are not intended to be comprehensive modeling tools, but are designed for modeling at the exploratory or conceptual level, for visual demonstration, and for educational purposes.

Comparison of ground-water flow model particle-tracking results and isotopic data in the Mojave River ground-water basin, southern California, USA

USGS Publications Warehouse

Izbicki, John A.; Stamos, Christina L.; Nishikawa, Tracy; Martin, Peter

2004-01-01

Flow-path and time-of-travel results for the Mojave River ground-water basin, southern California, calculated using the ground-water flow model MODFLOW and particle-tracking model MODPATH were similar to flow path and time-of-travel interpretations derived from delta-deuterium and carbon-14 data. Model and isotopic data both show short flow paths and young ground-water ages throughout the floodplain aquifer along most the Mojave River. Longer flow paths and older ground-water ages as great as 10,000 years before present were measured and simulated in the floodplain aquifer near the Mojave Valley. Model and isotopic data also show movement of water between the floodplain and regional aquifer and subsequent discharge of water from the river to dry lakes in some areas. It was not possible to simulate the isotopic composition of ground-water in the regional aquifer away from the front of the San Gabriel and San Bernardino Mountains - because recharge in these areas does not occur under the present-day climatic conditions used for calibration of the model.

Multiregional input-output model for the evaluation of Spanish water flows.

PubMed

Cazcarro, Ignacio; Duarte, Rosa; Sánchez Chóliz, Julio

2013-01-01

We construct a multiregional input-output model for Spain, in order to evaluate the pressures on the water resources, virtual water flows, and water footprints of the regions, and the water impact of trade relationships within Spain and abroad. The study is framed with those interregional input-output models constructed to study water flows and impacts of regions in China, Australia, Mexico, or the UK. To build our database, we reconcile regional IO tables, national and regional accountancy of Spain, trade and water data. Results show an important imbalance between origin of water resources and final destination, with significant water pressures in the South, Mediterranean, and some central regions. The most populated and dynamic regions of Madrid and Barcelona are important drivers of water consumption in Spain. Main virtual water exporters are the South and Central agrarian regions: Andalusia, Castile-La Mancha, Castile-Leon, Aragon, and Extremadura, while the main virtual water importers are the industrialized regions of Madrid, Basque country, and the Mediterranean coast. The paper shows the different location of direct and indirect consumers of water in Spain and how the economic trade and consumption pattern of certain areas has significant impacts on the availability of water resources in other different and often drier regions.

Regression modeling of ground-water flow

USGS Publications Warehouse

Cooley, R.L.; Naff, R.L.

1985-01-01

Nonlinear multiple regression methods are developed to model and analyze groundwater flow systems. Complete descriptions of regression methodology as applied to groundwater flow models allow scientists and engineers engaged in flow modeling to apply the methods to a wide range of problems. Organization of the text proceeds from an introduction that discusses the general topic of groundwater flow modeling, to a review of basic statistics necessary to properly apply regression techniques, and then to the main topic: exposition and use of linear and nonlinear regression to model groundwater flow. Statistical procedures are given to analyze and use the regression models. A number of exercises and answers are included to exercise the student on nearly all the methods that are presented for modeling and statistical analysis. Three computer programs implement the more complex methods. These three are a general two-dimensional, steady-state regression model for flow in an anisotropic, heterogeneous porous medium, a program to calculate a measure of model nonlinearity with respect to the regression parameters, and a program to analyze model errors in computed dependent variables such as hydraulic head. (USGS)

Integrated water flow model and modflow-farm process: A comparison of theory, approaches, and features of two integrated hydrologic models

USGS Publications Warehouse

Dogrul, Emin C.; Schmid, Wolfgang; Hanson, Randall T.; Kadir, Tariq; Chung, Francis

2016-01-01

Effective modeling of conjunctive use of surface and subsurface water resources requires simulation of land use-based root zone and surface flow processes as well as groundwater flows, streamflows, and their interactions. Recently, two computer models developed for this purpose, the Integrated Water Flow Model (IWFM) from the California Department of Water Resources and the MODFLOW with Farm Process (MF-FMP) from the US Geological Survey, have been applied to complex basins such as the Central Valley of California. As both IWFM and MFFMP are publicly available for download and can be applied to other basins, there is a need to objectively compare the main approaches and features used in both models. This paper compares the concepts, as well as the method and simulation features of each hydrologic model pertaining to groundwater, surface water, and landscape processes. The comparison is focused on the integrated simulation of water demand and supply, water use, and the flow between coupled hydrologic processes. The differences in the capabilities and features of these two models could affect the outcome and types of water resource problems that can be simulated.

Digital-computer model of ground-water flow in Tooele Valley, Utah

USGS Publications Warehouse

Razem, Allan C.; Bartholoma, Scott D.

1980-01-01

A two-dimensional, finite-difference digital-computer model was used to simulate the ground-water flow in the principal artesian aquifer in Tooele Valley, Utah. The parameters used in the model were obtained through field measurements and tests, from historical records, and by trial-and-error adjustments. The model was calibrated against observed water-level changes that occurred during 1941-50, 1951-60, 1961-66, 1967-73, and 1974-78. The reliability of the predictions is good in most parts of the valley, as is shown by the ability of the model to match historical water-level changes.

Documentation of a finite-element two-layer model for simulation of ground-water flow

USGS Publications Warehouse

Mallory, Michael J.

1979-01-01

This report documents a finite-element model for simulation of ground-water flow in a two-aquifer system where the two aquifers are coupled by a leakage term that represents flow through a confining layer separating the two aquifers. The model was developed by Timothy J. Durbin (U.S. Geological Survey) for use in ground-water investigations in southern California. The documentation assumes that the reader is familiar with the physics of ground-water flow, numerical methods of solving partial-differential equations, and the FORTRAN IV computer language. It was prepared as part of the investigations made by the U.S. Geological Survey in cooperation with the San Bernardino Valley Municipal Water District. (Kosco-USGS)

Three-dimensional hydrogeologic framework model for use with a steady-state numerical ground-water flow model of the Death Valley regional flow system, Nevada and California

USGS Publications Warehouse

Belcher, Wayne R.; Faunt, Claudia C.; D'Agnese, Frank A.

2002-01-01

The U.S. Geological Survey, in cooperation with the Department of Energy and other Federal, State, and local agencies, is evaluating the hydrogeologic characteristics of the Death Valley regional ground-water flow system. The ground-water flow system covers an area of about 100,000 square kilometers from latitude 35? to 38?15' North to longitude 115? to 118? West, with the flow system proper comprising about 45,000 square kilometers. The Death Valley regional ground-water flow system is one of the larger flow systems within the Southwestern United States and includes in its boundaries the Nevada Test Site, Yucca Mountain, and much of Death Valley. Part of this study includes the construction of a three-dimensional hydrogeologic framework model to serve as the foundation for the development of a steady-state regional ground-water flow model. The digital framework model provides a computer-based description of the geometry and composition of the hydrogeologic units that control regional flow. The framework model of the region was constructed by merging two previous framework models constructed for the Yucca Mountain Project and the Environmental Restoration Program Underground Test Area studies at the Nevada Test Site. The hydrologic characteristics of the region result from a currently arid climate and complex geology. Interbasinal regional ground-water flow occurs through a thick carbonate-rock sequence of Paleozoic age, a locally thick volcanic-rock sequence of Tertiary age, and basin-fill alluvium of Tertiary and Quaternary age. Throughout the system, deep and shallow ground-water flow may be controlled by extensive and pervasive regional and local faults and fractures. The framework model was constructed using data from several sources to define the geometry of the regional hydrogeologic units. These data sources include (1) a 1:250,000-scale hydrogeologic-map compilation of the region; (2) regional-scale geologic cross sections; (3) borehole information, and (4

Measurement Of Multiphase Flow Water Fraction And Water-cut

NASA Astrophysics Data System (ADS)

Xie, Cheng-gang

2007-06-01

This paper describes a microwave transmission multiphase flow water-cut meter that measures the amplitude attenuation and phase shift across a pipe diameter at multiple frequencies using cavity-backed antennas. The multiphase flow mixture permittivity and conductivity are derived from a unified microwave transmission model for both water- and oil-continuous flows over a wide water-conductivity range; this is far beyond the capability of microwave-resonance-based sensors currently on the market. The water fraction and water cut are derived from a three-component gas-oil-water mixing model using the mixture permittivity or the mixture conductivity and an independently measured mixture density. Water salinity variations caused, for example, by changing formation water or formation/injection water breakthrough can be detected and corrected using an online water-conductivity tracking technique based on the interpretation of the mixture permittivity and conductivity, simultaneously measured by a single-modality microwave sensor.

Experimental Measurements and Mathematical Modeling of Static and Dynamic Characteristics of Water Flow in a Long Pipe

NASA Astrophysics Data System (ADS)

Jablonska, J.; Kozubkova, M.

2017-08-01

Static and dynamic characteristics of flow in technical practice are very important and serious problem and can be solved by experimental measurement or mathematical modeling. Unsteady flow presents time changes of the flow and water hammer can be an example of this phenomenon. Water hammer is caused by rapid changes in the water flow by means the closure or opening of the control valve. The authors deal with by hydraulic hammer at the multiphase flow (water and air), its one-dimensional modeling (Matlab SimHydraulics) and modeling with the use of the finite volume method (Ansys Fluent) in article. The circuit elements are defined by static and dynamic characteristics. The results are verified with measurements. The article evaluates different approaches, their advantages, disadvantages and specifics in solving of water hammer.

The Development of a Real Time Surface Water Flow Model to Protect Public Water Intakes in West Virginia

NASA Astrophysics Data System (ADS)

Zegre, N.; Strager, M.

2015-12-01

In January of 2014 West Virginia experienced a chemical spill upstream of a public water intake on the Elk River near Charleston, West Virginia that made the water unusable for 300,000 people for weeks. In response to this disaster, state officials enacted legislation to protect the future public water intake locations by requiring the delineation of zones of critical concern that extend a five hour travel time above the intakes. Each zone is defined by the travel time and buffered along the river mainstem and tributary locations to identify future potential threats to the water supply. While this approach helps to identify potential problems before they occur, the need existed to be able to respond to a spill with information regarding the real travel time of a spill to an intake with consideration of actual stream flow at the time of the spill. This study developed a real time surface flow model to protect the public water intakes using both regional and seasonal variables. Bayesian statistical inference enabled confidence levels to be placed on flow estimates and used to show the probability for the time steps as water approached an public water intake. The flow model has been incorporated into both a smartphone app and web-based tool for better emergency response and management of water resources throughout the state.

Snow mass and river flows modelled using GRACE total water storage observations

NASA Astrophysics Data System (ADS)

Wang, S.

2017-12-01

Snow mass and river flow measurements are difficult and less accurate in cold regions due to the hash environment. Floods in cold regions are commonly a result of snowmelt during the spring break-up. Flooding is projected to increase with climate change in many parts of the world. Forecasting floods from snowmelt remains a challenge due to scarce and quality issues in basin-scale snow observations and lack of knowledge for cold region hydrological processes. This study developed a model for estimating basin-level snow mass (snow water equivalent SWE) and river flows using the total water storage (TWS) observations from the Gravity Recovery and Climate Experiment (GRACE) satellite mission. The SWE estimation is based on mass balance approach which is independent of in situ snow gauge observations, thus largely eliminates the limitations and uncertainties with traditional in situ or remote sensing snow estimates. The model forecasts river flows by simulating surface runoff from snowmelt and the corresponding baseflow from groundwater discharge. Snowmelt is predicted using a temperature index model. Baseflow is predicted using a modified linear reservoir model. The model also quantifies the hysteresis between the snowmelt and the streamflow rates, or the lump time for water travel in the basin. The model was applied to the Red River Basin, the Mackenzie River Basin, and the Hudson Bay Lowland Basins in Canada. The predicted river flows were compared with the observed values at downstream hydrometric stations. The results were also compared to that for the Lower Fraser River obtained in a separate study to help better understand the roles of environmental factors in determining flood and their variations with different hydroclimatic conditions. This study advances the applications of space-based time-variable gravity measurements in cold region snow mass estimation, river flow and flood forecasting. It demonstrates a relatively simple method that only needs GRACE TWS

Three-dimensional numerical model of ground-water flow in northern Utah Valley, Utah County, Utah

USGS Publications Warehouse

Gardner, Philip M.

2009-01-01

A three-dimensional, finite-difference, numerical model was developed to simulate ground-water flow in northern Utah Valley, Utah. The model includes expanded areal boundaries as compared to a previous ground-water flow model of the valley and incorporates more than 20 years of additional hydrologic data. The model boundary was generally expanded to include the bedrock in the surrounding mountain block as far as the surface-water divide. New wells have been drilled in basin-fill deposits near the consolidated-rock boundary. Simulating the hydrologic conditions within the bedrock allows for improved simulation of the effect of withdrawal from these wells. The inclusion of bedrock also allowed for the use of a recharge model that provided an alternative method for spatially distributing areal recharge over the mountains.The model was calibrated to steady- and transient-state conditions. The steady-state simulation was developed and calibrated by using hydrologic data that represented average conditions for 1947. The transient-state simulation was developed and calibrated by using hydrologic data collected from 1947 to 2004. Areally, the model grid is 79 rows by 70 columns, with variable cell size. Cells throughout most of the model domain represent 0.3 mile on each side. The largest cells are rectangular with dimensions of about 0.3 by 0.6 mile. The largest cells represent the mountain block on the eastern edge of the model domain where the least hydrologic data are available. Vertically, the aquifer system is divided into 4 layers which incorporate 11 hydrogeologic units. The model simulates recharge to the ground-water flow system as (1) infiltration of precipitation over the mountain block, (2) infiltration of precipitation over the valley floor, (3) infiltration of unconsumed irrigation water from fields, lawns, and gardens, (4) seepage from streams and canals, and (5) subsurface inflow from Cedar Valley. Discharge of ground water is simulated by the model to (1

UNSAT-H Version 2. 0: Unsaturated soil water and heat flow model

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

Fayer, M.J.; Jones, T.L.

1990-04-01

This report documents UNSAT-H Version 2.0, a model for calculating water and heat flow in unsaturated media. The documentation includes the bases for the conceptual model and its numerical implementation, benchmark test cases, example simulations involving layered soils and plant transpiration, and the code listing. Waste management practices at the Hanford Site have included disposal of low-level wastes by near-surface burial. Predicting the future long-term performance of any such burial site in terms of migration of contaminants requires a model capable of simulating water flow in the unsaturated soils above the buried waste. The model currently used to meet thismore » need is UNSAT-H. This model was developed at Pacific Northwest Laboratory to assess water dynamics of near-surface, waste-disposal sites at the Hanford Site. The code is primarily used to predict deep drainage as a function of such environmental conditions as climate, soil type, and vegetation. UNSAT-H is also used to simulate the effects of various practices to enhance isolation of wastes. 66 refs., 29 figs., 7 tabs.« less

Do Europa's Mountains Have Roots? Modeling Flow Along the Ice-Water Interface

NASA Astrophysics Data System (ADS)

Cutler, B. B.; Goodman, J. C.

2016-12-01

Are topographic features on the surface of Europa and other icy worlds isostatically compensated by variations in shell thickness (Airy isostasy)? This is only possible if variations in shell thickness can remain stable over geologic time. In this work we demonstrate that local shell thickness perturbations will relax due to viscous flow in centuries. We present a model of Europa's ice crust which includes thermal conduction, viscous flow of ice, and a mobile ice/water interface: the topography along the ice-water interface varies in response to melting, freezing, and ice flow. Temperature-dependent viscosity, conductivity, and density lead to glacier-like flow along the base of the ice shell, as well as solid-state convection in its interior. We considered both small scale processes, such as an isostatically-compensated ridge or lenticula, or heat flux from a hydrothermal plume; and a larger model focusing on melting and flow on the global scale. Our local model shows that ice-basal topographic features 5 kilometers deep and 4 kilometers wide can be filled in by glacial flow in about 200 years; even very large cavities can be infilled in 1000 years. "Hills" (locally thick areas) are removed faster than "holes". If a strong local heat flux (10x global average) is applied to the base of the ice, local melting will be prevented by rapid inflow of ice from nearby. On the large scale, global ice flow from the thick cool pole to the warmer and thinner equator removes global-scale topography in about 1 Ma; melting and freezing from this process may lead to a coupled feedback with the ocean flow. We find that glacial flow at the base of the ice shell is so rapid that Europa's ice-water interface is likely to be very flat. Local surface topography probably cannot be isostatically compensated by thickness variations: Europa's mountains may have no roots.

MODFLOW-2005 : the U.S. Geological Survey modular ground-water model--the ground-water flow process

USGS Publications Warehouse

Harbaugh, Arlen W.

2005-01-01

This report presents MODFLOW-2005, which is a new version of the finite-difference ground-water model commonly called MODFLOW. Ground-water flow is simulated using a block-centered finite-difference approach. Layers can be simulated as confined or unconfined. Flow associated with external stresses, such as wells, areal recharge, evapotranspiration, drains, and rivers, also can be simulated. The report includes detailed explanations of physical and mathematical concepts on which the model is based, an explanation of how those concepts are incorporated in the modular structure of the computer program, instructions for using the model, and details of the computer code. The modular structure consists of a MAIN Program and a series of highly independent subroutines. The subroutines are grouped into 'packages.' Each package deals with a specific feature of the hydrologic system that is to be simulated, such as flow from rivers or flow into drains, or with a specific method of solving the set of simultaneous equations resulting from the finite-difference method. Several solution methods are incorporated, including the Preconditioned Conjugate-Gradient method. The division of the program into packages permits the user to examine specific hydrologic features of the model independently. This also facilitates development of additional capabilities because new packages can be added to the program without modifying the existing packages. The input and output systems of the computer program also are designed to permit maximum flexibility. The program is designed to allow other capabilities, such as transport and optimization, to be incorporated, but this report is limited to describing the ground-water flow capability. The program is written in Fortran 90 and will run without modification on most computers that have a Fortran 90 compiler.

Integrated Water Flow Model (IWFM), A Tool For Numerically Simulating Linked Groundwater, Surface Water And Land-Surface Hydrologic Processes

NASA Astrophysics Data System (ADS)

Dogrul, E. C.; Brush, C. F.; Kadir, T. N.

2006-12-01

The Integrated Water Flow Model (IWFM) is a comprehensive input-driven application for simulating groundwater flow, surface water flow and land-surface hydrologic processes, and interactions between these processes, developed by the California Department of Water Resources (DWR). IWFM couples a 3-D finite element groundwater flow process and 1-D land surface, lake, stream flow and vertical unsaturated-zone flow processes which are solved simultaneously at each time step. The groundwater flow system is simulated as a multilayer aquifer system with a mixture of confined and unconfined aquifers separated by semiconfining layers. The groundwater flow process can simulate changing aquifer conditions (confined to unconfined and vice versa), subsidence, tile drains, injection wells and pumping wells. The land surface process calculates elemental water budgets for agricultural, urban, riparian and native vegetation classes. Crop water demands are dynamically calculated using distributed soil properties, land use and crop data, and precipitation and evapotranspiration rates. The crop mix can also be automatically modified as a function of pumping lift using logit functions. Surface water diversions and groundwater pumping can each be specified, or can be automatically adjusted at run time to balance water supply with water demand. The land-surface process also routes runoff to streams and deep percolation to the unsaturated zone. Surface water networks are specified as a series of stream nodes (coincident with groundwater nodes) with specified bed elevation, conductance and stage-flow relationships. Stream nodes are linked to form stream reaches. Stream inflows at the model boundary, surface water diversion locations, and one or more surface water deliveries per location are specified. IWFM routes stream flows through the network, calculating groundwater-surface water interactions, accumulating inflows from runoff, and allocating available stream flows to meet specified or

SIMPLE MODEL OF ICE SEGREGATION USING AN ANALYTIC FUNCTION TO MODEL HEAT AND SOIL-WATER FLOW.

USGS Publications Warehouse

Hromadka, T.V.; Guymon, G.L.

1984-01-01

This paper reports on the development of a simple two-dimensional model of coupled heat and soil-water flow in freezing or thawing soil. The model also estimates ice-segregation (frost-heave) evolution. Ice segregation in soil results from water drawn into a freezing zone by hydraulic gradients created by the freezing of soil-water. Thus, with a favorable balance between the rate of heat extraction and the rate of water transport to a freezing zone, segregated ice lenses may form.

Numerical modeling of coupled water flow and heat transport in soil and snow

Treesearch

Thijs J. Kelleners; Jeremy Koonce; Rose Shillito; Jelle Dijkema; Markus Berli; Michael H. Young; John M. Frank; William Massman

2016-01-01

A one-dimensional vertical numerical model for coupled water flow and heat transport in soil and snow was modified to include all three phases of water: vapor, liquid, and ice. The top boundary condition in the model is driven by incoming precipitation and the surface energy balance. The model was applied to three different terrestrial systems: A warm desert bare...

The Impact of Rhizosphere Processes on Water Flow and Root Water Uptake

NASA Astrophysics Data System (ADS)

Schwartz, Nimrod; Kroener, Eva; Carminati, Andrea; Javaux, Mathieu

2015-04-01

For many years, the rhizosphere, which is the zone of soil in the vicinity of the roots and which is influenced by the roots, is known as a unique soil environment with different physical, biological and chemical properties than those of the bulk soil. Indeed, in recent studies it has been shown that root exudate and especially mucilage alter the hydraulic properties of the soil, and that drying and wetting cycles of mucilage result in non-equilibrium water dynamics in the rhizosphere. While there are experimental evidences and simplified 1D model for those concepts, an integrated model that considers rhizosphere processes with a detailed model for water and roots flow is absent. Therefore, the objective of this work is to develop a 3D physical model of water flow in the soil-plant continuum that take in consideration root architecture and rhizosphere specific properties. Ultimately, this model will enhance our understanding on the impact of processes occurring in the rhizosphere on water flow and root water uptake. To achieve this objective, we coupled R-SWMS, a detailed 3D model for water flow in soil and root system (Javaux et al 2008), with the rhizosphere model developed by Kroener et al (2014). In the new Rhizo-RSWMS model the rhizosphere hydraulic properties differ from those of the bulk soil, and non-equilibrium dynamics between the rhizosphere water content and pressure head is also considered. We simulated a wetting scenario. The soil was initially dry and it was wetted from the top at a constant flow rate. The model predicts that, after infiltration the water content in the rhizosphere remained lower than in the bulk soil (non-equilibrium), but over time water infiltrated into the rhizosphere and eventually the water content in the rhizosphere became higher than in the bulk soil. These results are in qualitative agreement with the available experimental data on water dynamics in the rhizosphere. Additionally, the results show that rhizosphere processes

Use of a ground-water flow model with particle tracking to evaluate ground-water vulnerability, Clark County, Washington

USGS Publications Warehouse

Snyder, D.T.; Wilkinson, J.M.; Orzol, L.L.

1996-01-01

A ground-water flow model was used in conjunction with particle tracking to evaluate ground-water vulnerability in Clark County, Washington. Using the particle-tracking program, particles were placed in every cell of the flow model (about 60,000 particles) and tracked backwards in time and space upgradient along flow paths to their recharge points. A new computer program was developed that interfaces the results from a particle-tracking program with a geographic information system (GIS). The GIS was used to display and analyze the particle-tracking results. Ground-water vulnerability was evaluated by selecting parts of the ground-water flow system and combining the results with ancillary information stored in the GIS to determine recharge areas, characteristics of recharge areas, downgradient impact of land use at recharge areas, and age of ground water. Maps of the recharge areas for each hydrogeologic unit illustrate the presence of local, intermediate, or regional ground-water flow systems and emphasize the three-dimensional nature of the ground-water flow system in Clark County. Maps of the recharge points for each hydrogeologic unit were overlaid with maps depicting aquifer sensitivity as determined by DRASTIC (a measure of the pollution potential of ground water, based on the intrinsic characteristics of the near-surface unsaturated and saturated zones) and recharge from on-site waste-disposal systems. A large number of recharge areas were identified, particularly in southern Clark County, that have a high aquifer sensitivity, coincide with areas of recharge from on-site waste-disposal systems, or both. Using the GIS, the characteristics of the recharge areas were related to the downgradient parts of the ground-water system that will eventually receive flow that has recharged through these areas. The aquifer sensitivity, as indicated by DRASTIC, of the recharge areas for downgradient parts of the flow system was mapped for each hydrogeologic unit. A number of

Documentation of the Surface-Water Routing (SWR1) Process for modeling surface-water flow with the U.S. Geological Survey Modular Ground-Water Model (MODFLOW-2005)

USGS Publications Warehouse

Hughes, Joseph D.; Langevin, Christian D.; Chartier, Kevin L.; White, Jeremy T.

2012-01-01

A flexible Surface-Water Routing (SWR1) Process that solves the continuity equation for one-dimensional and two-dimensional surface-water flow routing has been developed for the U.S. Geological Survey three-dimensional groundwater model, MODFLOW-2005. Simple level- and tilted-pool reservoir routing and a diffusive-wave approximation of the Saint-Venant equations have been implemented. Both methods can be implemented in the same model and the solution method can be simplified to represent constant-stage elements that are functionally equivalent to the standard MODFLOW River or Drain Package boundary conditions. A generic approach has been used to represent surface-water features (reaches) and allows implementation of a variety of geometric forms. One-dimensional geometric forms include rectangular, trapezoidal, and irregular cross section reaches to simulate one-dimensional surface-water features, such as canals and streams. Two-dimensional geometric forms include reaches defined using specified stage-volume-area-perimeter (SVAP) tables and reaches covering entire finite-difference grid cells to simulate two-dimensional surface-water features, such as wetlands and lakes. Specified SVAP tables can be used to represent reaches that are smaller than the finite-difference grid cell (for example, isolated lakes), or reaches that cannot be represented accurately using the defined top of the model. Specified lateral flows (which can represent point and distributed flows) and stage-dependent rainfall and evaporation can be applied to each reach. The SWR1 Process can be used with the MODFLOW Unsaturated Zone Flow (UZF1) Package to permit dynamic simulation of runoff from the land surface to specified reaches. Surface-water/groundwater interactions in the SWR1 Process are mathematically defined to be a function of the difference between simulated stages and groundwater levels, and the specific form of the reach conductance equation used in each reach. Conductance can be

Application of nonlinear-regression methods to a ground-water flow model of the Albuquerque Basin, New Mexico

USGS Publications Warehouse

Tiedeman, C.R.; Kernodle, J.M.; McAda, D.P.

1998-01-01

This report documents the application of nonlinear-regression methods to a numerical model of ground-water flow in the Albuquerque Basin, New Mexico. In the Albuquerque Basin, ground water is the primary source for most water uses. Ground-water withdrawal has steadily increased since the 1940's, resulting in large declines in water levels in the Albuquerque area. A ground-water flow model was developed in 1994 and revised and updated in 1995 for the purpose of managing basin ground- water resources. In the work presented here, nonlinear-regression methods were applied to a modified version of the previous flow model. Goals of this work were to use regression methods to calibrate the model with each of six different configurations of the basin subsurface and to assess and compare optimal parameter estimates, model fit, and model error among the resulting calibrations. The Albuquerque Basin is one in a series of north trending structural basins within the Rio Grande Rift, a region of Cenozoic crustal extension. Mountains, uplifts, and fault zones bound the basin, and rock units within the basin include pre-Santa Fe Group deposits, Tertiary Santa Fe Group basin fill, and post-Santa Fe Group volcanics and sediments. The Santa Fe Group is greater than 14,000 feet (ft) thick in the central part of the basin. During deposition of the Santa Fe Group, crustal extension resulted in development of north trending normal faults with vertical displacements of as much as 30,000 ft. Ground-water flow in the Albuquerque Basin occurs primarily in the Santa Fe Group and post-Santa Fe Group deposits. Water flows between the ground-water system and surface-water bodies in the inner valley of the basin, where the Rio Grande, a network of interconnected canals and drains, and Cochiti Reservoir are located. Recharge to the ground-water flow system occurs as infiltration of precipitation along mountain fronts and infiltration of stream water along tributaries to the Rio Grande; subsurface

MODFLOW-2000, The U.S. Geological Survey Modular Ground-Water Model - User Guide to Modularization Concepts and the Ground-Water Flow Process

USGS Publications Warehouse

Harbaugh, Arlen W.; Banta, Edward R.; Hill, Mary C.; McDonald, Michael G.

2000-01-01

MODFLOW is a computer program that numerically solves the three-dimensional ground-water flow equation for a porous medium by using a finite-difference method. Although MODFLOW was designed to be easily enhanced, the design was oriented toward additions to the ground-water flow equation. Frequently there is a need to solve additional equations; for example, transport equations and equations for estimating parameter values that produce the closest match between model-calculated heads and flows and measured values. This report documents a new version of MODFLOW, called MODFLOW-2000, which is designed to accommodate the solution of equations in addition to the ground-water flow equation. This report is a user's manual. It contains an overview of the old and added design concepts, documents one new package, and contains input instructions for using the model to solve the ground-water flow equation.

Simulation of cylindrical flow to a well using the U.S. Geological Survey Modular Finite-Difference Ground-Water Flow Model

USGS Publications Warehouse

Reilly, Thomas E.; Harbaugh, Arlen W.

1993-01-01

Cylindrical (axisymmetric) flow to a well is an important specialized topic of ground-water hydraulics and has been applied by many investigators to determine aquifer properties and determine heads and flows in the vicinity of the well. A recent modification to the U.S. Geological Survey Modular Three-Dimensional Finite-Difference Ground-Water Flow Model provides the opportunity to simulate axisymmetric flow to a well. The theory involves the conceptualization of a system of concentric shells that are capable of reproducing the large variations in gradient in the vicinity of the well by decreasing their area in the direction of the well. The computer program presented serves as a preprocessor to the U.S. Geological Survey model by creating the input data file needed to implement the axisymmetric conceptualization. Data input requirements to this preprocessor are described, and a comparison with a known analytical solution indicates that the model functions appropriately.

A three-dimensional numerical model of predevelopment conditions in the Death Valley regional ground-water flow system, Nevada and California

USGS Publications Warehouse

D'Agnese, Frank A.; O'Brien, G. M.; Faunt, C.C.; Belcher, W.R.; San Juan, C.

2002-01-01

In the early 1990's, two numerical models of the Death Valley regional ground-water flow system were developed by the U.S. Department of Energy. In general, the two models were based on the same basic hydrogeologic data set. In 1998, the U.S. Department of Energy requested that the U.S. Geological Survey develop and maintain a ground-water flow model of the Death Valley region in support of U.S. Department of Energy programs at the Nevada Test Site. The purpose of developing this 'second-generation' regional model was to enhance the knowledge an understanding of the ground-water flow system as new information and tools are developed. The U.S. Geological Survey also was encouraged by the U.S. Department of Energy to cooperate to the fullest extent with other Federal, State, and local entities in the region to take advantage of the benefits of their knowledge and expertise. The short-term objective of the Death Valley regional ground-water flow system project was to develop a steady-state representation of the predevelopment conditions of the ground-water flow system utilizing the two geologic interpretations used to develop the previous numerical models. The long-term objective of this project was to construct and calibrate a transient model that simulates the ground-water conditions of the study area over the historical record that utilizes a newly interpreted hydrogeologic conceptual model. This report describes the result of the predevelopment steady-state model construction and calibration. The Death Valley regional ground-water flow system is situated within the southern Great Basin, a subprovince of the Basin and Range physiographic province, bounded by latitudes 35 degrees north and 38 degrees 15 minutes north and by longitudes 115 and 118 degrees west. Hydrology in the region is a result of both the arid climatic conditions and the complex geology. Ground-water flow generally can be described as dominated by interbasinal flow and may be conceptualized as

Simulation of the Regional Ground-Water-Flow System and Ground-Water/Surface-Water Interaction in the Rock River Basin, Wisconsin

USGS Publications Warehouse

Juckem, Paul F.

2009-01-01

A regional, two-dimensional, areal ground-water-flow model was developed to simulate the ground-water-flow system and ground-water/surface-water interaction in the Rock River Basin. The model was developed by the U.S. Geological Survey (USGS), in cooperation with the Rock River Coalition. The objectives of the regional model were to improve understanding of the ground-water-flow system and to develop a tool suitable for evaluating the effects of potential regional water-management programs. The computer code GFLOW was used because of the ease with which the model can simulate ground-water/surface-water interactions, provide a framework for simulating regional ground-water-flow systems, and be refined in a stepwise fashion to incorporate new data and simulate ground-water-flow patterns at multiple scales. The ground-water-flow model described in this report simulates the major hydrogeologic features of the modeled area, including bedrock and surficial aquifers, ground-water/surface-water interactions, and ground-water withdrawals from high-capacity wells. The steady-state model treats the ground-water-flow system as a single layer with hydraulic conductivity and base elevation zones that reflect the distribution of lithologic groups above the Precambrian bedrock and a regionally significant confining unit, the Maquoketa Formation. In the eastern part of the Basin where the shale-rich Maquoketa Formation is present, deep ground-water flow in the sandstone aquifer below the Maquoketa Formation was not simulated directly, but flow into this aquifer was incorporated into the GFLOW model from previous work in southeastern Wisconsin. Recharge was constrained primarily by stream base-flow estimates and was applied uniformly within zones guided by regional infiltration estimates for soils. The model includes average ground-water withdrawals from 1997 to 2006 for municipal wells and from 1997 to 2005 for high-capacity irrigation, industrial, and commercial wells. In addition

The One-Water Hydrologic Flow Model - The next generation in fully integrated hydrologic simulation software

NASA Astrophysics Data System (ADS)

Boyce, S. E.; Hanson, R. T.

2015-12-01

The One-Water Hydrologic Flow Model (MF-OWHM) is a MODFLOW-based integrated hydrologic flow model that is the most complete version, to date, of the MODFLOW family of hydrologic simulators needed for the analysis of a broad range of conjunctive-use issues. MF-OWHM fully links the movement and use of groundwater, surface water, and imported water for consumption by agriculture and natural vegetation on the landscape, and for potable and other uses within a supply-and-demand framework. MF-OWHM is based on the Farm Process for MODFLOW-2005 combined with Local Grid Refinement, Streamflow Routing, Surface-water Routing Process, Seawater Intrusion, Riparian Evapotranspiration, and the Newton-Raphson solver. MF-OWHM also includes linkages for deformation-, flow-, and head-dependent flows; additional observation and parameter options for higher-order calibrations; and redesigned code for facilitation of self-updating models and faster simulation run times. The next version of MF-OWHM, currently under development, will include a new surface-water operations module that simulates dynamic reservoir operations, the conduit flow process for karst aquifers and leaky pipe networks, a new subsidence and aquifer compaction package, and additional features and enhancements to enable more integration and cross communication between traditional MODFLOW packages. By retaining and tracking the water within the hydrosphere, MF-OWHM accounts for "all of the water everywhere and all of the time." This philosophy provides more confidence in the water accounting by the scientific community and provides the public a foundation needed to address wider classes of problems such as evaluation of conjunctive-use alternatives and sustainability analysis, including potential adaptation and mitigation strategies, and best management practices. By Scott E. Boyce and Randall T. Hanson

Multiphase flow experiments, mathematical modeling and numerical simulation of the water - gas - solute movement

NASA Astrophysics Data System (ADS)

Li, Y.; Ma, X.; Su, N.

2013-12-01

The movement of water and solute into and through the vadose zone is, in essence, an issue of immiscible displacement in pore-space network of a soil. Therefore, multiphase flow and transport in porous media, referring to three medium: air, water, and the solute, pose one of the largest unresolved challenges for porous medium fluid seepage. However, this phenomenon has always been largely neglected. It is expected that a reliable analysis model of the multi-phase flow in soil can truly reflect the process of natural movement about the infiltration, which is impossible to be observed directly. In such cases, geophysical applications of the nuclear magnetic resonance (NMR) provides the opportunity to measure the water movements into soils directly over a large scale from tiny pore to regional scale, accordingly enable it available both on the laboratory and on the field. In addition, the NMR provides useful information about the pore space properties. In this study, we proposed both laboratory and field experiments to measure the multi-phase flow parameters, together with optimize the model in computer programming based on the fractional partial differential equations (fPDE). In addition, we establish, for the first time, an infiltration model including solute flowing with water, which has huge influence on agriculture and soil environment pollution. Afterwards, with data collected from experiments, we simulate the model and analyze the spatial variability of parameters. Simulations are also conducted according to the model to evaluate the effects of airflow on water infiltration and other effects such as solute and absorption. It has significant meaning to oxygen irrigation aiming to higher crop yield, and shed more light into the dam slope stability. In summary, our framework is a first-time model added in solute to have a mathematic analysis with the fPDE and more instructive to agriculture activities.

Experimental and numerical modelling of surface water-groundwater flow and pollution interactions under tidal forcing

NASA Astrophysics Data System (ADS)

Spanoudaki, Katerina; Bockelmann-Evans, Bettina; Schaefer, Florian; Kampanis, Nikolaos; Nanou-Giannarou, Aikaterini; Stamou, Anastasios; Falconer, Roger

2015-04-01

Surface water and groundwater are integral components of the hydrologic continuum and the interaction between them affects both their quantity and quality. However, surface water and groundwater are often considered as two separate systems and are analysed independently. This separation is partly due to the different time scales, which apply in surface water and groundwater flows and partly due to the difficulties in measuring and modelling their interactions (Winter et al., 1998). Coastal areas in particular are a difficult hydrologic environment to represent with a mathematical model due to the large number of contributing hydrologic processes. Accurate prediction of interactions between coastal waters, groundwater and neighbouring wetlands, for example, requires the use of integrated surface water-groundwater models. In the past few decades a large number of mathematical models and field methods have been developed in order to quantify the interaction between groundwater and hydraulically connected surface water bodies. Field studies may provide the best data (Hughes, 1995) but are usually expensive and involve too many parameters. In addition, the interpretation of field measurements and linking with modelling tools often proves to be difficult. In contrast, experimental studies are less expensive and provide controlled data. However, experimental studies of surface water-groundwater interaction are less frequently encountered in the literature than filed studies (e.g. Ebrahimi et al., 2007; Kuan et al., 2012; Sparks et al., 2013). To this end, an experimental model has been constructed at the Hyder Hydraulics Laboratory at Cardiff University to enable measurements to be made of groundwater transport through a sand embankment between a tidal water body such as an estuary and a non-tidal water body such as a wetland. The transport behaviour of a conservative tracer was studied for a constant water level on the wetland side of the embankment, while running a

Hypothetical Modeling of Redox Conditions Within a Complex Ground-Water Flow Field in a Glacial Setting

USGS Publications Warehouse

Feinstein, Daniel T.; Thomas, Mary Ann

2009-01-01

This report describes a modeling approach for studying how redox conditions evolve under the influence of a complex ground-water flow field. The distribution of redox conditions within a flow system is of interest because of the intrinsic susceptibility of an aquifer to redox-sensitive, naturally occurring contaminants - such as arsenic - as well as anthropogenic contaminants - such as chlorinated solvents. The MODFLOW-MT3D-RT3D suite of code was applied to a glacial valley-fill aquifer to demonstrate a method for testing the interaction of flow patterns, sources of reactive organic carbon, and availability of electron acceptors in controlling redox conditions. Modeling results show how three hypothetical distributions of organic carbon influence the development of redox conditions in a water-supply aquifer. The distribution of strongly reduced water depends on the balance between the rate of redox reactions and the capability of different parts of the flow system to transmit oxygenated water. The method can take account of changes in the flow system induced by pumping that result in a new distribution of reduced water.

The value of long-term monitoring in the development of ground-water-flow models

USGS Publications Warehouse

Feinstein, Daniel T.; Hart, David J.; Krohelski, James T.

2004-01-01

As environmental issues have come to the forefront of public concern, so has the awareness of the importance of ground water in the overall water cycle and as a source of the Nation’s drinking water. Heightened interest has spawned a host of scientific enterprises (Taylor and Alley, 2001). Some activities are directed toward collection of water-level data and related information to monitor the physical and chemical state of the resource. Other activities are directed at interpretive studies undertaken, for example, to optimize the location of new water-supply wells or to protect rivers and lakes fed by ground water. An important type of interpretive study is the computer ground-water-flow model that inte- grates field data in a mathematical framework. Long-term, systematic collection of hydro- logic data is crucial to the construction and testing of ground-water models so that they can reproduce the evolution of flow systems and forecast future conditions. 

Simulation of ground-water flow, surface-water flow, and a deep sewer tunnel system in the Menomonee Valley, Milwaukee, Wisconsin

USGS Publications Warehouse

Dunning, C.P.; Feinstein, D.T.; Hunt, R.J.; Krohelski, J.T.

2004-01-01

Numerical models were constructed for simulation of ground-water flow in the Menomonee Valley Brownfield, in Milwaukee, Wisconsin. An understanding of ground-water flow is necessary to develop an efficient program to sample ground water for contaminants. Models were constructed in a stepwise fashion, beginning with a regional, single-layer, analytic-element model (GFLOW code) that provided boundary conditions for a local, eight layer, finite-difference model (MODFLOW code) centered on the Menomonee Valley Brownfield. The primary source of ground water to the models is recharge over the model domains; primary sinks for ground water within the models are surface-water features and the Milwaukee Metropolitan Sewerage District Inline Storage System (ISS). Calibration targets were hydraulic heads, surface-water fluxes, vertical gradients, and ground-water infiltration to the ISS. Simulation of ground-water flow by use of the MODFLOW model indicates that about 73 percent of recharge within the MODFLOW domain circulates to the ISS and 27 percent discharges to gaining surface-water bodies. In addition, infiltration to the ISS comes from the following sources: 36 percent from recharge within the model domain, 45 percent from lateral flow into the domain, 15 percent from Lake Michigan, and 4 percent from other surface-water bodies. Particle tracking reveals that the median traveltime from the recharge point to surface-water features is 8 years; the median time to the ISS is 255 years. The traveltimes to the ISS are least over the northern part of the valley, where dolomite is near the land surface. The distribution of traveltimes in the MODFLOW simulation is greatly influenced by the effective porosity values assigned to the various lithologies.

Determinants of virtual water flows in the Mediterranean.

PubMed

Fracasso, Andrea; Sartori, Martina; Schiavo, Stefano

2016-02-01

The aim of the paper is to investigate the main determinants of the bilateral virtual water (water used in the production of a commodity or service) flows associated with international trade in agricultural goods across the Mediterranean basin. We consider the bilateral gross flows of virtual water in the area and study what export-specific and import-specific factors are significantly associated with virtual water flows. We follow a sequential approach. Through a gravity model of trade, we obtain a "refined" version of the variable we aim to explain, one that is free of the amount of flows due to pair-specific factors affecting bilateral trade flows and that fully reflects the impact of country-specific determinants of virtual water trade. A number of country-specific potential explanatory variables, ranging from water endowments to trade barriers, from per capita GDP to irrigation prices, is presented and tested. To identify the variables that help to explain the bilateral flows of virtual water, we adopt a model selection procedure based on model averaging. Our findings confirm one of the main controversial results in the literature: larger water endowments do not necessarily lead to a larger 'export' of virtual water, as one could expect. We also find some evidence that higher water irrigation prices reduce (increase) virtual water 'exports' ('imports'). Copyright © 2015 Elsevier B.V. All rights reserved.

Evaluation of 2D shallow-water model for spillway flow with a complex geometry

USDA-ARS?s Scientific Manuscript database

Although the two-dimensional (2D) shallow water model is formulated based on several assumptions such as hydrostatic pressure distribution and vertical velocity is negligible, as a simple alternative to the complex 3D model, it has been used to compute water flows in which these assumptions may be ...

Fracture control of ground water flow and water chemistry in a rock aquitard

USGS Publications Warehouse

Eaton, T.T.; Anderson, M.P.; Bradbury, K.R.

2007-01-01

There are few studies on the hydrogeology of sedimentary rock aquitards although they are important controls in regional ground water flow systems. We formulate and test a three-dimensional (3D) conceptual model of ground water flow and hydrochemistry in a fractured sedimentary rock aquitard to show that flow dynamics within the aquitard are more complex than previously believed. Similar conceptual models, based on regional observations and recently emerging principles of mechanical stratigraphy in heterogeneous sedimentary rocks, have previously been applied only to aquifers, but we show that they are potentially applicable to aquitards. The major elements of this conceptual model, which is based on detailed information from two sites in the Maquoketa Formation in southeastern Wisconsin, include orders of magnitude contrast between hydraulic diffusivity (K/Ss) of fractured zones and relatively intact aquitard rock matrix, laterally extensive bedding-plane fracture zones extending over distances of over 10 km, very low vertical hydraulic conductivity of thick shale-rich intervals of the aquitard, and a vertical hydraulic head profile controlled by a lateral boundary at the aquitard subcrop, where numerous surface water bodies dominate the shallow aquifer system. Results from a 3D numerical flow model based on this conceptual model are consistent with field observations, which did not fit the typical conceptual model of strictly vertical flow through an aquitard. The 3D flow through an aquitard has implications for predicting ground water flow and for planning and protecting water supplies. ?? 2007 National Ground Water Association.

Ensemble kalman filtering to perform data assimilation with soil water content probes and pedotransfer functions in modeling water flow in variably saturated soils

USDA-ARS?s Scientific Manuscript database

Data from modern soil water contents probes can be used for data assimilation in soil water flow modeling, i.e. continual correction of the flow model performance based on observations. The ensemble Kalman filter appears to be an appropriate method for that. The method requires estimates of the unce...

Development of a stream-aquifer numerical flow model to assess river water management under water scarcity in a Mediterranean basin.

PubMed

Mas-Pla, Josep; Font, Eva; Astui, Oihane; Menció, Anna; Rodríguez-Florit, Agustí; Folch, Albert; Brusi, David; Pérez-Paricio, Alfredo

2012-12-01

Stream flow, as a part of a basin hydrological cycle, will be sensible to water scarcity as a result of climate change. Stream vulnerability should then be evaluated as a key component of the basin water budget. Numerical flow modeling has been applied to an alluvial formation in a small mountain basin to evaluate the stream-aquifer relationship under these future scenarios. The Arbúcies River basin (116 km(2)) is located in the Catalan Inner Basins (NE Spain) and its lower reach, which is related to an alluvial aquifer, usually becomes dry during the summer period. This study seeks to determine the origin of such discharge losses whether from natural stream leakage and/or induced capture due to groundwater withdrawal. Our goal is also investigating how discharge variations from the basin headwaters, representing potential effects of climate change, may affect stream flow, aquifer recharge, and finally environmental preservation and human supply. A numerical flow model of the alluvial aquifer, based on MODFLOW and especially in the STREAM routine, reproduced the flow system after the usual calibration. Results indicate that, in the average, stream flow provides more than 50% of the water inputs to the alluvial aquifer, being responsible for the amount of stored water resources and for satisfying groundwater exploitation for human needs. Detailed simulations using daily time-steps permit setting threshold values for the stream flow entering at the beginning of the studied area so surface discharge is maintained along the whole watercourse and ecological flow requirements are satisfied as well. The effects of predicted rainfall and temperature variations on the Arbúcies River alluvial aquifer water balance are also discussed from the outcomes of the simulations. Finally, model results indicate the relevance of headwater discharge management under future climate scenarios to preserve downstream hydrological processes. They also point out that small mountain basins

A comparison of hydrologic models for ecological flows and water availability

USGS Publications Warehouse

Caldwell, Peter V; Kennen, Jonathan G.; Sun, Ge; Kiang, Julie E.; Butcher, John B; Eddy, Michelle C; Hay, Lauren E.; LaFontaine, Jacob H.; Hain, Ernie F.; Nelson, Stacy C; McNulty, Steve G

2015-01-01

Robust hydrologic models are needed to help manage water resources for healthy aquatic ecosystems and reliable water supplies for people, but there is a lack of comprehensive model comparison studies that quantify differences in streamflow predictions among model applications developed to answer management questions. We assessed differences in daily streamflow predictions by four fine-scale models and two regional-scale monthly time step models by comparing model fit statistics and bias in ecologically relevant flow statistics (ERFSs) at five sites in the Southeastern USA. Models were calibrated to different extents, including uncalibrated (level A), calibrated to a downstream site (level B), calibrated specifically for the site (level C) and calibrated for the site with adjusted precipitation and temperature inputs (level D). All models generally captured the magnitude and variability of observed streamflows at the five study sites, and increasing level of model calibration generally improved performance. All models had at least 1 of 14 ERFSs falling outside a +/−30% range of hydrologic uncertainty at every site, and ERFSs related to low flows were frequently over-predicted. Our results do not indicate that any specific hydrologic model is superior to the others evaluated at all sites and for all measures of model performance. Instead, we provide evidence that (1) model performance is as likely to be related to calibration strategy as it is to model structure and (2) simple, regional-scale models have comparable performance to the more complex, fine-scale models at a monthly time step.

A Modular Three-Dimensional Finite-Difference Ground-Water Flow Model

USGS Publications Warehouse

McDonald, Michael G.; Harbaugh, Arlen W.; Guo, Weixing; Lu, Guoping

1988-01-01

This report presents a finite-difference model and its associated modular computer program. The model simulates flow in three dimensions. The report includes detailed explanations of physical and mathematical concepts on which the model is based and an explanation of how those concepts are incorporated in the modular structure of the computer program. The modular structure consists of a Main Program and a series of highly independent subroutines called 'modules.' The modules are grouped into 'packages.' Each package deals with a specific feature of the hydrologic system which is to be simulated, such as flow from rivers or flow into drains, or with a specific method of solving linear equations which describe the flow system, such as the Strongly Implicit Procedure or Slice-Successive Overrelaxation. The division of the program into modules permits the user to examine specific hydrologic features of the model independently. This also facilita development of additional capabilities because new packages can be added to the program without modifying the existing packages. The input and output systems of the computer program are also designed to permit maximum flexibility. Ground-water flow within the aquifer is simulated using a block-centered finite-difference approach. Layers can be simulated as confined, unconfined, or a combination of confined and unconfined. Flow associated with external stresses, such as wells, areal recharge, evapotranspiration, drains, and streams, can also be simulated. The finite-difference equations can be solved using either the Strongly Implicit Procedure or Slice-Successive Overrelaxation. The program is written in FORTRAN 77 and will run without modification on most computers that have a FORTRAN 77 compiler. For each program ,module, this report includes a narrative description, a flow chart, a list of variables, and a module listing.

Modeling and Simulation of Hydropower Station Diversion System's characteristic line method by introducing water head to flow feedback

NASA Astrophysics Data System (ADS)

Guangwen, Xu; Xi, Li; Ze, Yao

2018-06-01

To solve the damping problem of water hammer wave in the modeling method of water diversion system of hydropower station, this paper introduces the feedback regulation technology from head to flow, that is: A fixed water head is taken out for flow feedback, and the following conclusions are obtained through modeling and simulation. Adjusting the feedback coefficient F of the water head to the flow rate can change the damping characteristic of the system, which can simulate the attenuation process of the water shock wave in the true water diversion pipeline. Even if a small feedback coefficient is introduced, the damping effect of the system is very obvious, but it has little effect on the amplitude of the first water shock wave after the transition process. Therefore, it is feasible and reasonable to introduce water head to flow rate feedback coefficient F in hydraulic turbine diversion system.

Fracture control of ground water flow and water chemistry in a rock aquitard.

PubMed

Eaton, Timothy T; Anderson, Mary P; Bradbury, Kenneth R

2007-01-01

There are few studies on the hydrogeology of sedimentary rock aquitards although they are important controls in regional ground water flow systems. We formulate and test a three-dimensional (3D) conceptual model of ground water flow and hydrochemistry in a fractured sedimentary rock aquitard to show that flow dynamics within the aquitard are more complex than previously believed. Similar conceptual models, based on regional observations and recently emerging principles of mechanical stratigraphy in heterogeneous sedimentary rocks, have previously been applied only to aquifers, but we show that they are potentially applicable to aquitards. The major elements of this conceptual model, which is based on detailed information from two sites in the Maquoketa Formation in southeastern Wisconsin, include orders of magnitude contrast between hydraulic diffusivity (K/S(s)) of fractured zones and relatively intact aquitard rock matrix, laterally extensive bedding-plane fracture zones extending over distances of over 10 km, very low vertical hydraulic conductivity of thick shale-rich intervals of the aquitard, and a vertical hydraulic head profile controlled by a lateral boundary at the aquitard subcrop, where numerous surface water bodies dominate the shallow aquifer system. Results from a 3D numerical flow model based on this conceptual model are consistent with field observations, which did not fit the typical conceptual model of strictly vertical flow through an aquitard. The 3D flow through an aquitard has implications for predicting ground water flow and for planning and protecting water supplies.

A modification of the finite-difference model for simulation of two dimensional ground-water flow to include surface-ground water relationships

USGS Publications Warehouse

Ozbilgin, M.M.; Dickerman, D.C.

1984-01-01

The two-dimensional finite-difference model for simulation of groundwater flow was modified to enable simulation of surface-water/groundwater interactions during periods of low streamflow. Changes were made to the program code in order to calculate surface-water heads for, and flow either to or from, contiguous surface-water bodies; and to allow for more convenient data input. Methods of data input and output were modified and entries (RSORT and HDRIVER) were added to the COEF and CHECKI subroutines to calculate surface-water heads. A new subroutine CALC was added to the program which initiates surface-water calculations. If CALC is not specified as a simulation option, the program runs the original version. The subroutines which solve the ground-water flow equations were not changed. Recharge, evapotranspiration, surface-water inflow, number of wells, pumping rate, and pumping duration can be varied for any time period. The Manning formula was used to relate stream depth and discharge in surface-water streams. Interactions between surface water and ground water are represented by the leakage term in the ground-water flow and surface-water mass balance equations. Documentation includes a flow chart, data deck instructions, input data, output summary, and program listing. Numerical results from the modified program are in good agreement with published analytical results. (USGS)

Three Principles of Water Flow in Soils

NASA Astrophysics Data System (ADS)

Guo, L.; Lin, H.

2016-12-01

Knowledge of water flow in soils is crucial to understanding terrestrial hydrological cycle, surface energy balance, biogeochemical dynamics, ecosystem services, contaminant transport, and many other Critical Zone processes. However, due to the complex and dynamic nature of non-uniform flow, reconstruction and prediction of water flow in natural soils remain challenging. This study synthesizes three principles of water flow in soils that can improve modeling water flow in soils of various complexity. The first principle, known as the Darcy's law, came to light in the 19th century and suggested a linear relationship between water flux density and hydraulic gradient, which was modified by Buckingham for unsaturated soils. Combining mass balance and the Buckingham-Darcy's law, L.A. Richards quantitatively described soil water change with space and time, i.e., Richards equation. The second principle was proposed by L.A. Richards in the 20th century, which described the minimum pressure potential needed to overcome surface tension of fluid and initiate water flow through soil-air interface. This study extends this principle to encompass soil hydrologic phenomena related to varied interfaces and microscopic features and provides a more cohesive explanation of hysteresis, hydrophobicity, and threshold behavior when water moves through layered soils. The third principle is emerging in the 21st century, which highlights the complex and evolving flow networks embedded in heterogeneous soils. This principle is summarized as: Water moves non-uniformly in natural soils with a dual-flow regime, i.e., it follows the least-resistant or preferred paths when "pushed" (e.g., by storms) or "attracted" (e.g., by plants) or "restricted" (e.g., by bedrock), but moves diffusively into the matrix when "relaxed" (e.g., at rest) or "touched" (e.g., adsorption). The first principle is a macroscopic view of steady-state water flow, the second principle is a microscopic view of interface

Computing nonhydrostatic shallow-water flow over steep terrain

USGS Publications Warehouse

Denlinger, R.P.; O'Connell, D. R. H.

2008-01-01

Flood and dambreak hazards are not limited to moderate terrain, yet most shallow-water models assume that flow occurs over gentle slopes. Shallow-water flow over rugged or steep terrain often generates significant nonhydrostatic pressures, violating the assumption of hydrostatic pressure made in most shallow-water codes. In this paper, we adapt a previously published nonhydrostatic granular flow model to simulate shallow-water flow, and we solve conservation equations using a finite volume approach and an Harten, Lax, Van Leer, and Einfeldt approximate Riemann solver that is modified for a sloping bed and transient wetting and drying conditions. To simulate bed friction, we use the law of the wall. We test the model by comparison with an analytical solution and with results of experiments in flumes that have steep (31??) or shallow (0.3??) slopes. The law of the wall provides an accurate prediction of the effect of bed roughness on mean flow velocity over two orders of magnitude of bed roughness. Our nonhydrostatic, law-of-the-wall flow simulation accurately reproduces flume measurements of front propagation speed, flow depth, and bed-shear stress for conditions of large bed roughness. ?? 2008 ASCE.

Validation of a Ground-Water Flow Model of the Mississippi River Valley Alluvial Aquifer Using Water-Level and Water-Use Data for 1998-2005 and Evaluation of Water-Use Scenarios

USGS Publications Warehouse

Gillip, Jonathan A.; Czarnecki, John B.

2009-01-01

A ground-water flow model of the Mississippi River Valley alluvial aquifer in eastern Arkansas, developed in 2003 to simulate the period of 1918-98, was validated with the addition of water-level and water-use data that extended the observation period to 2005. The original model (2003) was calibrated using water-level observations from 1972, 1982, 1992, and 1998, and water-use data through 1997. The original model subsequently was used to simulate water levels from 1999 to 2049 and showed that simulation of continued pumping at the 1997 water-use rate could not be sustained indefinitely without causing dry cells in the model. After publication of the original ground-water flow model, a total of 3,616 water-level observations from 698 locations measured during the period of 1998 to 2005 became available. Additionally, water-use data were compiled and used for the same period, totaling 290,005 discrete water-use values from 43,440 wells with as many as 39,169 wells pumping in any one year. Total pumping (which is primarily agricultural) for this 8-year period was about 2.3 trillion cubic feet of water and was distributed over approximately 10,340 square miles within the model area. An updated version of the original ground-water flow model was used to simulate the period of 1998-2005 with the additional water-level and water-use data. Water-level observations for 1998-2005 ranged from 74 to 293 feet above National Geodetic Vertical Datum of 1929 across the model area. The maximum water-level residual (observed minus simulated water-level values) for the 3,616 water-level observations was 52 feet, the minimum water-level residual was 60 feet, the average annual root mean squared error was 8.2 feet, and the annual average absolute residual was 6.0 feet. A correlation coefficient value of 0.96 was calculated for the line of best fit for observed to simulated water levels for the combined 1998-2005 dataset, indicating a good fit to the data and an acceptable validation

Water Flow Performance of a Superscale Model of the Fastrac Liquid Oxygen Pump

NASA Technical Reports Server (NTRS)

Skelley, Stephen; Zoladz, Thomas

1999-01-01

As part of the National Aeronautics and Space Administration's ongoing effort to lower the cost of access to space, the Marshall Space Flight Center has developed a rocket engine with 60,000 pounds of thrust for use on the Reusable Launch Vehicle technology demonstrator slated for launch in 2000. This gas generator cycle engine, known as the Fastrac engine, uses liquid oxygen and RP-1 for propellants and includes single stage liquid oxygen and RP-1 pumps and a single stage supersonic turbine on a common shaft. The turbopump design effort included the first use and application of new suction capability prediction codes and three-dimensional blade generation codes in an attempt to reduce the turbomachinery design and certification costs typically associated with rocket engine development. To verify the pump's predicted cavitation performance, a water flow test of a superscale model of the Fastrac liquid oxygen pump was conducted to experimentally evaluate the liquid oxygen pump's performance at and around the design point. The water flow test article replicated the flow path of the Fastrac liquid oxygen pump in a 1.582x scale model, including scaled seal clearances for correct leakage flow at a model operating speed of 5000 revolutions per minute. Flow entered the 3-blade axial-flow inducer, transitioned to a shrouded, 6-blade radial impeller, and discharged into a vaneless radial diffuser and collection volute. The test article included approximately 50 total and static pressure measurement locations as well as flush-mounted, high frequency pressure transducers for complete mapping of the pressure environment. The primary objectives of the water flow test were to measure the steady-state and dynamic pressure environment of the liquid oxygen pump versus flow coefficient, suction specific speed, and back face leakage flow rate. Results showed excellent correlation between the predicted and experimentally measured pump head rise at low suction specific speeds. Likewise

Revised Multi-Node Well (MNW2) Package for MODFLOW Ground-Water Flow Model

USGS Publications Warehouse

Konikow, Leonard F.; Hornberger, George Z.; Halford, Keith J.; Hanson, Randall T.; Harbaugh, Arlen W.

2009-01-01

Wells that are open to multiple aquifers can provide preferential pathways to flow and solute transport that short-circuit normal fluid flowlines. Representing these features in a regional flow model can produce a more realistic and reliable simulation model. This report describes modifications to the Multi-Node Well (MNW) Package of the U.S. Geological Survey (USGS) three-dimensional ground-water flow model (MODFLOW). The modifications build on a previous version and add several new features, processes, and input and output options. The input structure of the revised MNW (MNW2) is more well-centered than the original verion of MNW (MNW1) and allows the user to easily define hydraulic characteristics of each multi-node well. MNW2 also allows calculations of additional head changes due to partial penetration effects, flow into a borehole through a seepage face, changes in well discharge related to changes in lift for a given pump, and intraborehole flows with a pump intake located at any specified depth within the well. MNW2 also offers an improved capability to simulate nonvertical wells. A new output option allows selected multi-node wells to be designated as 'observation wells' for which changes in selected variables with time will be written to separate output files to facilitate postprocessing. MNW2 is compatible with the MODFLOW-2000 and MODFLOW-2005 versions of MODFLOW and with the version of MODFLOW that includes the Ground-Water Transport process (MODFLOW-GWT).

Description and field analysis of a coupled ground-water/surface-water flow model (MODFLOW/BRANCH) with modifications for structures and wetlands in southern Dade County, Florida

USGS Publications Warehouse

Swain, E.D.; Howie, Barbara; Dixon, Joann

1996-01-01

A coupled surface-water model (BRANCH) and ground-water model (MODFLOW) model were tested to simulate the interacting wetlands/surface-water/ ground-water system of southern Dade County. Several options created for the MODFLOW ground- ground-water model were used in representing this field situation. The primary option is the MODBRANCH interfacing software, which allows leakage to be accounted for between the MODFLOW ground-water model and the BRANCH dynamic model for simulation of flow in an interconnected network of open channels. A modification to an existing software routine, which is referred to as BCF2, allows cells in MODFLOW to rewet when dry--a requirement in representing the seasonal wetlands in Dade County. A companion to BCF2 is the modified evapotranspiration routine EVT2. The EVT2 routine changes the cells where evapotranspiration occurs, depending on which cells are wet. The Streamlink package represents direct connections between the canals and wetlands at locations where canals open directly into overland flow. Within the BRANCH model, the capability to represent the numerous hydraulic structures, gated spillways, gated culverts, and pumps was added. The application of these modifications to model surface-water/ground-water interactions in southern Dade County demonstrated the usefulness of the coupled MODFLOW/BRANCH model. Ground-water and surface-water flows are both simulated with dynamic models. Flow exchange between models, intermittent wetting and drying, evapotranspiration, and hydraulic structure operations are all represented appropriately. Comparison was made with a simulation using the RIV1 package instead of MODBRANCH to represent the canals. RIV1 represents the canals by user-defined stages, and computes leakage to the aquifer. Greater accuracy in reproducing measured ground- water heads was achieved with MODBRANCH, which also computes dynamic flow conditions in the canals, unlike RIV1. The surface-water integrated flow and transport

First status report on regional ground-water flow modeling for the Paradox Basin, Utah

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

Andrews, R.W.

1984-05-01

Regional ground-water flow within the principal hydrogeologic units of the Paradox Basin is evaluated by developing a conceptual model of the flow regime in the shallow aquifers and the deep-basin brine aquifers and testing these models using a three-dimensional, finite-difference flow code. Semiquantitative sensitivity analysis (a limited parametric study) is conducted to define the system response to changes in hydrologic properties or boundary conditions. A direct method for sensitivity analysis using an adjoint form of the flow equation is applied to the conceptualized flow regime in the Leadville limestone aquifer. All steps leading to the final results and conclusions aremore » incorporated in this report. The available data utilized in this study is summarized. The specific conceptual models, defining the areal and vertical averaging of litho-logic units, aquifer properties, fluid properties, and hydrologic boundary conditions, are described in detail. Two models were evaluated in this study: a regional model encompassing the hydrogeologic units above and below the Paradox Formation/Hermosa Group and a refined scale model which incorporated only the post Paradox strata. The results are delineated by the simulated potentiometric surfaces and tables summarizing areal and vertical boundary fluxes, Darcy velocities at specific points, and ground-water travel paths. Results from the adjoint sensitivity analysis include importance functions and sensitivity coefficients, using heads or the average Darcy velocities to represent system response. The reported work is the first stage of an ongoing evaluation of the Gibson Dome area within the Paradox Basin as a potential repository for high-level radioactive wastes.« less

CityWaterBalance: Track Flows of Water Through an Urban System

EPA Science Inventory

CityWaterBalance provides a reproducible workflow for studying an urban water system. The network of urban water flows and storages can be modeled and visualized. Any city may be modeled with preassembled data, but data for US cities can be gathered via web services using this p...

Ground-water models for water resources planning

USGS Publications Warehouse

Moore, John E.

1980-01-01

In the past decade hydrologists have emphasized the development of computer-based mathematical models to aid in the understanding of flow, the transport of solutes, transport of heat, and deformation in the groundwater system. These models have been used to provide information and predictions for water managers. Too frequently, groundwater was neglected in water-resource planning because managers believed that it could not be adequately evaluated in terms of availability, quality, and effect of development on surface water supplies. Now, however, with newly developed digital groundwater models, effects of development can be predicted. Such models have been used to predict hydrologic and quality changes under different stresses. These models have grown in complexity over the last 10 years from simple one-layer flow models to three-dimensional simulations of groundwater flow which may include solute transport, heat transport, effects of land subsidence, and encroachment of salt water. This paper illustrates, through case histories, how predictive groundwater models have provided the information needed for the sound planning and management of water resources in the United States. (USGS)

Quantifying water flow and retention in an unsaturated fracture-facial domain

USGS Publications Warehouse

Nimmo, John R.; Malek-Mohammadi, Siamak

2015-01-01

Hydrologically significant flow and storage of water occur in macropores and fractures that are only partially filled. To accommodate such processes in flow models, we propose a three-domain framework. Two of the domains correspond to water flow and water storage in a fracture-facial region, in addition to the third domain of matrix water. The fracture-facial region, typically within a fraction of a millimeter of the fracture wall, includes a flowing phase whose fullness is determined by the availability and flux of preferentially flowing water, and a static storage portion whose fullness is determined by the local matric potential. The flow domain can be modeled with the source-responsive preferential flow model, and the roughness-storage domain can be modeled with capillary relations applied on the fracture-facial area. The matrix domain is treated using traditional unsaturated flow theory. We tested the model with application to the hydrology of the Chalk formation in southern England, coherently linking hydrologic information including recharge estimates, streamflow, water table fluctuation, imaging by electron microscopy, and surface roughness. The quantitative consistency of the three-domain matrix-microcavity-film model with this body of diverse data supports the hypothesized distinctions and active mechanisms of the three domains and establishes the usefulness of this framework.

Water exchange, mixing and transient storage between a saturated karstic conduit and the surrounding aquifer: Groundwater flow modeling and inputs from stable water isotopes

NASA Astrophysics Data System (ADS)

Binet, S.; Joigneaux, E.; Pauwels, H.; Albéric, P.; Fléhoc, Ch.; Bruand, A.

2017-01-01

Water exchanges between a karstic conduit and the surrounding aquifer are driven by hydraulic head gradient at the interface between these two domains. The case-study presented in this paper investigates the impact of the geometry and interface conditions around a conduit on the spatial distribution of these exchanges. Isotopic (δ18O and δD), discharge and water head measurements were conducted at the resurgences of a karst system with a strong allogenic recharge component (Val d'Orléans, France), to estimate the amounts of water exchanged and the mixings between a saturated karstic conduit and the surrounding aquifer. The spatio-temporal variability of the observed exchanges was explored using a 2D coupled continuum-conduit flow model under saturated conditions (Feflow®). The inputs from the water heads and stable water isotopes in the groundwater flow model suggest that the amounts of water flowing from the aquifer are significant if the conduit flow discharges are less than the conduit flow capacity. This condition creates a spatial distribution of exchanges from upstream where the aquifer feeds the conduit (recharge area) to downstream where the conduit reaches its maximum discharge capacity and can feed the aquifer (discharge area). In the intermediate transport zone no exchange between the two domains takes place that brings a new criterion to delineate the vulnerable zones to surface water. On average, 4% of the water comes from the local recharge, 80% is recent river water and 16% is old river water. During the November 2008 flood, both isotopic signatures and model suggest that exchanges fluctuate around this steady state, limited when the river water level increases and intensified when the river water level decreases. The existence of old water from the river suggests a transient storage at the aquifer/conduit interface that can be considered as an underground hyporheic zone.

Ground Water Flow No Longer A Mystery

ERIC Educational Resources Information Center

Lehr, Jay H.; Pettyjohn, Wayne A.

1976-01-01

Examined are the physical characteristics of ground water movement. Some potential pollution problems are identified. Models are used to explain mathematical and hydraulic principles of flow toward a pumping well and an effluent stream, flow around and through lenticular beds, and effects of pumping on the water table. (Author/MR)

Water Flow Performance of a Superscale Model of the Fastrac Liquid Oxygen Pump

NASA Technical Reports Server (NTRS)

Skelley, Stephen; Zoladz, Thomas

2001-01-01

As part of the National Aeronautics and Space Administration's ongoing effort to lower the cost of access to space, the Marshall Space Flight Center has developed a rocket engine with 60,000 pounds of thrust for use on the Reusable Launch Vehicle technology demonstrator slated for launch in 2000. This gas generator cycle engine, known as the Fastrac engine, uses liquid oxygen and RP-1 for propellants and includes single stage liquid oxygen and RP-1 pumps and a single stage supersonic turbine on a common shaft. The turbopump design effort included the first use and application of new suction capability prediction codes and three-dimensional blade generation codes in an attempt to reduce the turbomachinery design and certification costs typically associated with rocket engine development. To verify the pump's predicted cavitation performance, a water flow test of a superscale model of the Fastrac liquid oxygen pump was conducted to experimentally evaluate the liquid oxygen pump's performance at and around the design point. The water flow test article replicated the flow path of the Fastrac liquid oxygen pump in a 1.582x scale model, including scaled seal clearances for correct leakage flow at a model operating speed of 5000 revolutions per minute. Flow entered the 3-blade axial-flow inducer, transitioned to a shrouded, 6- blade radial impeller, and discharged into a vaneless radial diffuser and collection volute. The test article included approximately 50 total and static pressure measurement locations as well as flush-mounted, high frequency pressure transducers for complete mapping of the pressure environment. The primary objectives of the water flow test were to measure the steady-state and dynamic pressure environment of the liquid oxygen pump versus flow coefficient, suction specific speed, and back face leakage flow rate. Initial results showed acceptable correlation between the predicted and experimentally measured pump head rise at low suction specific speeds

Worse than imagined: Unidentified virtual water flows in China.

PubMed

Cai, Beiming; Wang, Chencheng; Zhang, Bing

2017-07-01

The impact of virtual water flows on regional water scarcity in China had been deeply discussed in previous research. However, these studies only focused on water quantity, the impact of virtual water flows on water quality has been largely neglected. In this study, we incorporate the blue water footprint related with water quantity and grey water footprint related with water quality into virtual water flow analysis based on the multiregional input-output model of 2007. The results find that the interprovincial virtual flows accounts for 23.4% of China's water footprint. The virtual grey water flows are 8.65 times greater than the virtual blue water flows; the virtual blue water and grey water flows are 91.8 and 794.6 Gm 3 /y, respectively. The use of the indicators related with water quantity to represent virtual water flows in previous studies will underestimate their impact on water resources. In addition, the virtual water flows are mainly derived from agriculture, chemical industry and petroleum processing and the coking industry, which account for 66.8%, 7.1% and 6.2% of the total virtual water flows, respectively. Virtual water flows have intensified both quantity- and quality-induced water scarcity of export regions, where low-value-added but water-intensive and high-pollution goods are produced. Our study on virtual water flows can inform effective water use policy for both water resources and water pollution in China. Our methodology about virtual water flows also can be used in global scale or other countries if data available. Copyright © 2017 Elsevier Ltd. All rights reserved.

Ground-Water Flow Model for the Spokane Valley-Rathdrum Prairie Aquifer, Spokane County, Washington, and Bonner and Kootenai Counties, Idaho

USGS Publications Warehouse

Hsieh, Paul A.; Barber, Michael E.; Contor, Bryce A.; Hossain, Md. Akram; Johnson, Gary S.; Jones, Joseph L.; Wylie, Allan H.

2007-01-01

This report presents a computer model of ground-water flow in the Spokane Valley-Rathdrum Prairie (SVRP) aquifer in Spokane County, Washington, and Bonner and Kootenai Counties, Idaho. The aquifer is the sole source of drinking water for more than 500,000 residents in the area. In response to the concerns about the impacts of increased ground-water withdrawals resulting from recent and projected urban growth, a comprehensive study was initiated by the Idaho Department of Water Resources, the Washington Department of Ecology, and the U.S. Geological Survey to improve the understanding of ground-water flow in the aquifer and of the interaction between ground water and surface water. The ground-water flow model presented in this report is one component of this comprehensive study. The primary purpose of the model is to serve as a tool for analyzing aquifer inflows and outflows, simulating the effects of future changes in ground-water withdrawals from the aquifer, and evaluating aquifer management strategies. The scale of the model and the level of detail are intended for analysis of aquifer-wide water-supply issues. The SVRP aquifer model was developed by the Modeling Team formed within the comprehensive study. The Modeling Team consisted of staff and personnel working under contract with the Idaho Department of Water Resources, personnel working under contract with the Washington Department of Ecology, and staff of the U.S. Geological Survey. To arrive at a final model that has the endorsement of all team members, decisions on modeling approach, methodology, assumptions, and interpretations were reached by consensus. The ground-water flow model MODFLOW-2000 was used to simulate ground-water flow in the SVPR aquifer. The finite-difference model grid consists of 172 rows, 256 columns, and 3 layers. Ground-water flow was simulated from September 1990 through September 2005 using 181 stress periods of 1 month each. The areal extent of the model encompasses an area of

Benchmark testing of DIII-D neutral beam modeling with water flow calorimetry

DOE PAGES

Rauch, J. M.; Crowley, B. J.; Scoville, J. T.; ...

2016-06-02

Power loading on beamline components in the DIII-D neutral beam system is measured in this paper using water flow calorimetry. The results are used to benchmark beam transport models. Finally, anomalously high heat loads in the magnet region are investigated and a speculative hypothesis as to their origin is presented.

Modeling water quality, temperature, and flow in Link River, south-central Oregon

USGS Publications Warehouse

Sullivan, Annett B.; Rounds, Stewart A.

2016-09-09

The 2.1-km (1.3-mi) Link River connects Upper Klamath Lake to the Klamath River in south-central Oregon. A CE-QUAL-W2 flow and water-quality model of Link River was developed to provide a connection between an existing model of the upper Klamath River and any existing or future models of Upper Klamath Lake. Water-quality sampling at six locations in Link River was done during 2013–15 to support model development and to provide a better understanding of instream biogeochemical processes. The short reach and high velocities in Link River resulted in fast travel times and limited water-quality transformations, except for dissolved oxygen. Reaeration through the reach, especially at the falls in Link River, was particularly important in moderating dissolved oxygen concentrations that at times entered the reach at Link River Dam with marked supersaturation or subsaturation. This reaeration resulted in concentrations closer to saturation downstream at the mouth of Link River.

Ground-water models for water resource planning

USGS Publications Warehouse

Moore, J.E.

1983-01-01

In the past decade hydrogeologists have emphasized the development of computer-based mathematical models to aid in the understanding of flow, the transport of solutes, transport of heat, and deformation in the ground-water system. These models have been used to provide information and predictions for water managers. Too frequently, ground-water was neglected in water resource planning because managers believed that it could not be adequately evaluated in terms of availability, quality, and effect of development on surface-water supplies. Now, however, with newly developed digital ground-water models, effects of development can be predicted. Such models have been used to predict hydrologic and quality changes under different stresses. These models have grown in complexity over the last ten years from simple one-layer models to three-dimensional simulations of ground-water flow, which may include solute transport, heat transport, effects of land subsidence, and encroachment of saltwater. Case histories illustrate how predictive ground-water models have provided the information needed for the sound planning and management of water resources in the USA. ?? 1983 D. Reidel Publishing Company.

Development and application of a screening model for simulating regional ground-water flow in the St. Croix River basin, Minnesota and Wisconsin

USGS Publications Warehouse

Feinstein, Daniel T.; Buchwald, Cheryl A.; Dunning, Charles P.; Hunt, Randall J.

2006-01-01

A series of databases and an accompanying screening model were constructed by the U.S. Geological Survey, in cooperation with the National Park Service, to better understand the regional ground-water-flow system and its relation to stream drainage in the St. Croix River Basin. The St. Croix River and its tributaries drain about 8,000 square miles in northeastern Minnesota and northwestern Wisconsin. The databases contain information for the entire St. Croix River Basin pertaining to well logs, lithology, thickness of lithologic groups, ground-water levels, streamflow, and well pumpage. Maps and generalized cross sections created from the compiled data show the lithologic groups, extending from the water table to the crystalline bedrock, through which ground water flows. These lithologic groups are: fine-grained unconsolidated deposits; coarse-grained unconsolidated deposits; sandstone bedrock; carbonate bedrock; and other bedrock lithologies including shale, siltstone, conglomerate, and igneous intrusions. The steady-state screening model treats the ground-water-flow system as a single layer with transmissivity zones that reflect the distribution of lithologic groups, and with recharge zones that correspond to general areas of high or low evapotranspiration. The model includes representation of second- and higher-order streams and municipal and other high-capacity production wells. The analytic-element model code GFLOW was used to simulate the regional ground-water flow, the water-table surface across the St. Croix River Basin, and base-flow contributions from ground water to streams. In addition, the model routes tributary base flow through the stream network to the St. Croix River. The parameter-estimation inverse model UCODE was linked to the GFLOW model to select the combination of parameter values best able to match over 5,000 water-level measurements and base-flow estimates at 22 streamflow-gaging stations. Results from the calibrated screening model show

Finite-difference model to simulate the areal flow of saltwater and fresh water separated by an interface

USGS Publications Warehouse

Mercer, James W.; Larson, S.P.; Faust, Charles R.

1980-01-01

Model documentation is presented for a two-dimensional (areal) model capable of simulating ground-water flow of salt water and fresh water separated by an interface. The partial differential equations are integrated over the thicknesses of fresh water and salt water resulting in two equations describing the flow characteristics in the areal domain. These equations are approximated using finite-difference techniques and the resulting algebraic equations are solved for the dependent variables, fresh water head and salt water head. An iterative solution method was found to be most appropriate. The program is designed to simulate time-dependent problems such as those associated with the development of coastal aquifers, and can treat water-table conditions or confined conditions with steady-state leakage of fresh water. The program will generally be most applicable to the analysis of regional aquifer problems in which the zone between salt water and fresh water can be considered a surface (sharp interface). Example problems and a listing of the computer code are included. (USGS).

Analytical Model of Water Flow in Coal with Active Matrix

NASA Astrophysics Data System (ADS)

Siemek, Jakub; Stopa, Jerzy

2014-12-01

This paper presents new analytical model of gas-water flow in coal seams in one dimension with emphasis on interactions between water flowing in cleats and coal matrix. Coal as a flowing system, can be viewed as a solid organic material consisting of two flow subsystems: a microporous matrix and a system of interconnected macropores and fractures. Most of gas is accumulated in the microporous matrix, where the primary flow mechanism is diffusion. Fractures and cleats existing in coal play an important role as a transportation system for macro scale flow of water and gas governed by Darcy's law. The coal matrix can imbibe water under capillary forces leading to exchange of mass between fractures and coal matrix. In this paper new partial differential equation for water saturation in fractures has been formulated, respecting mass exchange between coal matrix and fractures. Exact analytical solution has been obtained using the method of characteristics. The final solution has very simple form that may be useful for practical engineering calculations. It was observed that the rate of exchange of mass between the fractures and the coal matrix is governed by an expression which is analogous to the Newton cooling law known from theory of heat exchange, but in present case the mass transfer coefficient depends not only on coal and fluid properties but also on time and position. The constant term of mass transfer coefficient depends on relation between micro porosity and macro porosity of coal, capillary forces, and microporous structure of coal matrix. This term can be expressed theoretically or obtained experimentally. W artykule zaprezentowano nowy model matematyczny przepływu wody i gazu w jednowymiarowej warstwie węglowej z uwzględnieniem wymiany masy między systemem szczelin i matrycą węglową. Węgiel jako system przepływowy traktowany jest jako układ o podwójnej porowatości i przepuszczalności, składający się z mikroporowatej matrycy węglowej oraz z

A novel star extraction method based on modified water flow model

NASA Astrophysics Data System (ADS)

Zhang, Hao; Niu, Yanxiong; Lu, Jiazhen; Ouyang, Zibiao; Yang, Yanqiang

2017-11-01

Star extraction is the essential procedure for attitude measurement of star sensor. The great challenge for star extraction is to segment star area exactly from various noise and background. In this paper, a novel star extraction method based on Modified Water Flow Model(MWFM) is proposed. The star image is regarded as a 3D terrain. The morphology is adopted for noise elimination and Tentative Star Area(TSA) selection. Star area can be extracted through adaptive water flowing within TSAs. This method can achieve accurate star extraction with improved efficiency under complex conditions such as loud noise and uneven backgrounds. Several groups of different types of star images are processed using proposed method. Comparisons with existing methods are conducted. Experimental results show that MWFM performs excellently under different imaging conditions. The star extraction rate is better than 95%. The star centroid accuracy is better than 0.075 pixels. The time-consumption is also significantly reduced.

Quasi 3D modeling of water flow and solute transport in vadose zone and groundwater

NASA Astrophysics Data System (ADS)

Yakirevich, A.; Kuznetsov, M.; Weisbrod, N.; Pachepsky, Y. A.

2013-12-01

The complexity of subsurface flow systems calls for a variety of concepts leading to the multiplicity of simplified flow models. One commonly used simplification is based on the assumption that lateral flow and transport in unsaturated zone is insignificant unless the capillary fringe is involved. In such cases the flow and transport in the unsaturated zone above groundwater level can be simulated as a 1D phenomenon, whereas through groundwater they are viewed as 2D or 3D phenomena. A new approach for a numerical scheme for 3D variably saturated flow and transport is presented. A Quasi-3D approach allows representing flow in the 'vadose zone - aquifer' system by a series of 1D Richards' equations solved in variably-saturated zone and by 3D-saturated flow equation in groundwater (modified MODFLOW code). The 1D and 3D equations are coupled at the phreatic surface in a way that aquifer replenishment is calculated using the Richards' equation, and solving for the moving water table does not require definition of the specific yield parameter. The 3D advection-dispersion equation is solved in the entire domain by the MT3D code. Using implicit finite differences approximation to couple processes in the vadose zone and groundwater provides mass conservation and increase of computational efficiency. The above model was applied to simulate the impact of irrigation on groundwater salinity in the Alto Piura aquifer (Northern Peru). Studies on changing groundwater quality in arid and semi-arid lands show that irrigation return flow is one of the major factors contributing to aquifer salinization. Existing mathematical models do not account explicitly for the solute recycling during irrigation on a daily scale. Recycling occurs throughout the unsaturated and saturated zones, as function of the solute mass extracted from pumping wells. Salt concentration in irrigation water is calculated at each time step as a function of concentration of both surface water and groundwater

Ground water flow modeling with sensitivity analyses to guide field data collection in a mountain watershed

USGS Publications Warehouse

Johnson, Raymond H.

2007-01-01

In mountain watersheds, the increased demand for clean water resources has led to an increased need for an understanding of ground water flow in alpine settings. In Prospect Gulch, located in southwestern Colorado, understanding the ground water flow system is an important first step in addressing metal loads from acid-mine drainage and acid-rock drainage in an area with historical mining. Ground water flow modeling with sensitivity analyses are presented as a general tool to guide future field data collection, which is applicable to any ground water study, including mountain watersheds. For a series of conceptual models, the observation and sensitivity capabilities of MODFLOW-2000 are used to determine composite scaled sensitivities, dimensionless scaled sensitivities, and 1% scaled sensitivity maps of hydraulic head. These sensitivities determine the most important input parameter(s) along with the location of observation data that are most useful for future model calibration. The results are generally independent of the conceptual model and indicate recharge in a high-elevation recharge zone as the most important parameter, followed by the hydraulic conductivities in all layers and recharge in the next lower-elevation zone. The most important observation data in determining these parameters are hydraulic heads at high elevations, with a depth of less than 100 m being adequate. Evaluation of a possible geologic structure with a different hydraulic conductivity than the surrounding bedrock indicates that ground water discharge to individual stream reaches has the potential to identify some of these structures. Results of these sensitivity analyses can be used to prioritize data collection in an effort to reduce time and money spend by collecting the most relevant model calibration data.

An integrated approach to conjunctive-use analysis with the one-water hydrologic flow model, MODFLOW-OWHM

USGS Publications Warehouse

Boyce, Scott E.; Hanson, Randall T.

2015-01-01

The MODFLOW-2005 (MF) family of hydrologic simulators has diverged into multiple versions designed for specific needs, thus limiting their use to their respective designs. The One-Water Hydrologic Flow Model (MF-OWHM v1.0) is an integrated hydrologic flow model that is an enhanced fusion of multiple MF versions. While maintaining compatibility with existing MF versions, MF-OWHM includes: linkages for coupled heads, flows, and deformation; facilitation of self-updating models, additional observation and parameter options for higher-order calibrations; and redesigned code for faster simulations. This first release of MF-OWHM incorporates MODFLOW-2005 and the Farm Process (MF-FMP2), with new features (FMP3), combined with Local Grid Refinement (MF-LGR), Streamflow Routing (SFR), Surfacewater Routing Process (SWR), Seawater Intrusion (SWI), Riparian Evapotranspiration (RIP-ET), the Newton Formulation (MF-NWT), and more. MF-OWHM represents a complete integrated hydrologic model that fully links the movement and use of groundwater, surface water, and imported water for consumption by agriculture and natural vegetation on the landscape, and for potable and other uses. By retaining and keeping track of the water during simulation of the hydrosphere, MF-OWHM accounts for “all of the water everywhere and all of the time.” This provides the foundation needed to address integrated hydrologic problems such as evaluation of conjunctive-use alternatives and sustainability analysis, including potential adaptation and mitigation strategies, and best management practices.

Documentation of the Unsaturated-Zone Flow (UZF1) Package for modeling Unsaturated Flow Between the Land Surface and the Water Table with MODFLOW-2005

USGS Publications Warehouse

Niswonger, Richard G.; Prudic, David E.; Regan, R. Steven

2006-01-01

Percolation of precipitation through unsaturated zones is important for recharge of ground water. Rain and snowmelt at land surface are partitioned into different pathways including runoff, infiltration, evapotranspiration, unsaturated-zone storage, and recharge. A new package for MODFLOW-2005 called the Unsaturated-Zone Flow (UZF1) Package was developed to simulate water flow and storage in the unsaturated zone and to partition flow into evapotranspiration and recharge. The package also accounts for land surface runoff to streams and lakes. A kinematic wave approximation to Richards? equation is solved by the method of characteristics to simulate vertical unsaturated flow. The approach assumes that unsaturated flow occurs in response to gravity potential gradients only and ignores negative potential gradients; the approach further assumes uniform hydraulic properties in the unsaturated zone for each vertical column of model cells. The Brooks-Corey function is used to define the relation between unsaturated hydraulic conductivity and water content. Variables used by the UZF1 Package include initial and saturated water contents, saturated vertical hydraulic conductivity, and an exponent in the Brooks-Corey function. Residual water content is calculated internally by the UZF1 Package on the basis of the difference between saturated water content and specific yield. The UZF1 Package is a substitution for the Recharge and Evapotranspiration Packages of MODFLOW-2005. The UZF1 Package differs from the Recharge Package in that an infiltration rate is applied at land surface instead of a specified recharge rate directly to ground water. The applied infiltration rate is further limited by the saturated vertical hydraulic conductivity. The UZF1 Package differs from the Evapotranspiration Package in that evapotranspiration losses are first removed from the unsaturated zone above the evapotranspiration extinction depth, and if the demand is not met, water can be removed

Mathematical Modeling of Fluid Flow in a Water Physical Model of an Aluminum Degassing Ladle Equipped with an Impeller-Injector

NASA Astrophysics Data System (ADS)

Gómez, Eudoxio Ramos; Zenit, Roberto; Rivera, Carlos González; Trápaga, Gerardo; Ramírez-Argáez, Marco A.

2013-04-01

In this work, a 3D numerical simulation using a Euler-Euler-based model implemented into a commercial CFD code was used to simulate fluid flow and turbulence structure in a water physical model of an aluminum ladle equipped with an impeller for degassing treatment. The effect of critical process parameters such as rotor speed, gas flow rate, and the point of gas injection (conventional injection through the shaft vs a novel injection through the bottom of the ladle) on the fluid flow and vortex formation was analyzed with this model. The commercial CFD code PHOENICS 3.4 was used to solve all conservation equations governing the process for this two-phase fluid flow system. The mathematical model was reasonably well validated against experimentally measured liquid velocity and vortex sizes in a water physical model built specifically for this investigation. From the results, it was concluded that the angular speed of the impeller is the most important parameter in promoting better stirred baths and creating smaller and better distributed bubbles in the liquid. The pumping effect of the impeller is increased as the impeller rotation speed increases. Gas flow rate is detrimental to bath stirring and diminishes the pumping effect of the impeller. Finally, although the injection point was the least significant variable, it was found that the "novel" injection improves stirring in the ladle.

Cavitating flow during water hammer using a generalized interface vaporous cavitation model

NASA Astrophysics Data System (ADS)

Sadafi, Mohamadhosein; Riasi, Alireza; Nourbakhsh, Seyed Ahmad

2012-10-01

In a transient flow simulation, column separation may occur when the calculated pressure head decreases to the saturated vapor pressure head in a computational grid. Abrupt valve closure or pump failure can result in a fast transient flow with column separation, potentially causing problems such as pipe failure, hydraulic equipment damage, cavitation or corrosion. This paper reports a numerical study of water hammer with column separation in a simple reservoir-pipeline-valve system and pumping station. The governing equations for two-phase transient flow in pipes are solved based on the method of characteristics (MOC) using a generalized interface vaporous cavitating model (GIVCM). The numerical results were compared with the experimental data for validation purposes, and the comparison indicated that the GIVCM describes the experimental results more accurately than the discrete vapor cavity model (DVCM). In particular, the GIVCM correlated better with the experimental data than the DVCM in terms of timing and pressure magnitude. The effects of geometric and hydraulic parameters on flow behavior in a pumping station with column separation were also investigated in this study.

Experimental investigation of the flow in a simplified model of water lubricated axial thrust bearing

NASA Astrophysics Data System (ADS)

Kirschner, O.; Ruprecht, A.; Riedelbauch, S.

2014-03-01

In hydropower plants the axial thrust bearing takes up the hydraulic axial thrust of the runner and, in case of vertical shafts, the entire weight of all rotating masses. The use of water lubricated bearings can eliminate the oil leakage risk possibly contaminating the environment. A complex flow is generated by the smaller film thickness due to the lower viscosity of water compared with oil. Measurements on a simplified hydrostatic axial trust bearing model were accomplished for validating CFD analysis of water lubricated bearings. In this simplified model, fixed pads are implemented and the width of the gap was enlarged to create a higher resolution in space for the measurements. Most parts of the model were manufactured from acrylic glass to get optical access for measurement with PIV. The focus of these measurements is on the flow within the space between two pads. Additional to the PIV- measurement, the pressure on the wall of the rotating disk is captured by pressure transducers. The model bearing measurement results are presented for varied operating conditions.

MODOPTIM: A general optimization program for ground-water flow model calibration and ground-water management with MODFLOW

USGS Publications Warehouse

Halford, Keith J.

2006-01-01

MODOPTIM is a non-linear ground-water model calibration and management tool that simulates flow with MODFLOW-96 as a subroutine. A weighted sum-of-squares objective function defines optimal solutions for calibration and management problems. Water levels, discharges, water quality, subsidence, and pumping-lift costs are the five direct observation types that can be compared in MODOPTIM. Differences between direct observations of the same type can be compared to fit temporal changes and spatial gradients. Water levels in pumping wells, wellbore storage in the observation wells, and rotational translation of observation wells also can be compared. Negative and positive residuals can be weighted unequally so inequality constraints such as maximum chloride concentrations or minimum water levels can be incorporated in the objective function. Optimization parameters are defined with zones and parameter-weight matrices. Parameter change is estimated iteratively with a quasi-Newton algorithm and is constrained to a user-defined maximum parameter change per iteration. Parameters that are less sensitive than a user-defined threshold are not estimated. MODOPTIM facilitates testing more conceptual models by expediting calibration of each conceptual model. Examples of applying MODOPTIM to aquifer-test analysis, ground-water management, and parameter estimation problems are presented.

Spatial variability analysis of combining the water quality and groundwater flow model to plan groundwater and surface water management in the Pingtung plain

NASA Astrophysics Data System (ADS)

Chen, Ching-Fang; Chen, Jui-Sheng; Jang, Cheng-Shin

2014-05-01

As a result of rapid economic growth in the Pingtung Plain, the use of groundwater resources has changed dramatically. The groundwater is quite rich in the Pingtung plain and the most important water sources. During the several decades, a substantial amount of groundwater has been pumped for the drinking, irrigation and aquaculture water supplies. However, because the sustainable use concept of groundwater resources is lack, excessive pumping of groundwater causes the occurrence of serious land subsidence and sea water intrusion. Thus, the management and conservation of groundwater resources in the Pingtung plain are considerably critical. This study aims to assess the conjunct use effect of groundwater and surface water in the Pingtung plain on recharge by reducing the amount of groundwater extraction. The groundwater quality variability and groundwater flow models are combined to spatially analyze potential zones of groundwater used for multi-purpose in the Pingtung Plain. First, multivariate indicator kriging (MVIK) is used to analyze spatial variability of groundwater quality based on drinking, aquaculture and irrigation water quality standards, and probabilistically delineate suitable zones in the study area. Then, the groundwater flow model, Processing MODFLOW (PMWIN), is adopted to simulate groundwater flow. The groundwater flow model must be conducted by the calibration and verification processes, and the regional groundwater recovery is discussed when specified water rights are replaced by surface water in the Pingtung plain. Finally, the most suitable zones of reducing groundwater use are determined for multi-purpose according to combining groundwater quality and quantity. The study results can establish a sound and low-impact management plan of groundwater resources utilization for the multi-purpose groundwater use, and prevent decreasing ground water tables, and the occurrence of land subsidence and sea water intrusion in the Pingtung plain.

A modular finite-element model (MODFE) for areal and axisymmetric ground-water-flow problems, Part 1: Model Description and User's Manual

USGS Publications Warehouse

Torak, L.J.

1993-01-01

A MODular, Finite-Element digital-computer program (MODFE) was developed to simulate steady or unsteady-state, two-dimensional or axisymmetric ground-water flow. Geometric- and hydrologic-aquifer characteristics in two spatial dimensions are represented by triangular finite elements and linear basis functions; one-dimensional finite elements and linear basis functions represent time. Finite-element matrix equations are solved by the direct symmetric-Doolittle method or the iterative modified, incomplete-Cholesky, conjugate-gradient method. Physical processes that can be represented by the model include (1) confined flow, unconfined flow (using the Dupuit approximation), or a combination of both; (2) leakage through either rigid or elastic confining beds; (3) specified recharge or discharge at points, along lines, and over areas; (4) flow across specified-flow, specified-head, or bead-dependent boundaries; (5) decrease of aquifer thickness to zero under extreme water-table decline and increase of aquifer thickness from zero as the water table rises; and (6) head-dependent fluxes from springs, drainage wells, leakage across riverbeds or confining beds combined with aquifer dewatering, and evapotranspiration. The report describes procedures for applying MODFE to ground-water-flow problems, simulation capabilities, and data preparation. Guidelines for designing the finite-element mesh and for node numbering and determining band widths are given. Tables are given that reference simulation capabilities to specific versions of MODFE. Examples of data input and model output for different versions of MODFE are provided.

A modular finite-element model (MODFE) for areal and axisymmetric ground-water-flow problems; Part 1, Model description and user's manual

USGS Publications Warehouse

Torak, Lynn J.

1992-01-01

A MODular, Finite-Element digital-computer program (MODFE) was developed to simulate steady or unsteady-state, two-dimensional or axisymmetric ground-water flow. Geometric- and hydrologic-aquifer characteristics in two spatial dimensions are represented by triangular finite elements and linear basis functions; one-dimensional finite elements and linear basis functions represent time. Finite-element matrix equations are solved by the direct symmetric-Doolittle method or the iterative modified, incomplete-Cholesky, conjugate-gradient method. Physical processes that can be represented by the model include (1) confined flow, unconfined flow (using the Dupuit approximation), or a combination of both; (2) leakage through either rigid or elastic confining beds; (3) specified recharge or discharge at points, along lines, and over areas; (4) flow across specified-flow, specified-head, or head-dependent boundaries; (5) decrease of aquifer thickness to zero under extreme water-table decline and increase of aquifer thickness from zero as the water table rises; and (6) head-dependent fluxes from springs, drainage wells, leakage across riverbeds or confining beds combined with aquifer dewatering, and evapotranspiration.The report describes procedures for applying MODFE to ground-water-flow problems, simulation capabilities, and data preparation. Guidelines for designing the finite-element mesh and for node numbering and determining band widths are given. Tables are given that reference simulation capabilities to specific versions of MODFE. Examples of data input and model output for different versions of MODFE are provided.

Construction of estimated flow- and load-duration curves for Kentucky using the Water Availability Tool for Environmental Resources (WATER)

USGS Publications Warehouse

Unthank, Michael D.; Newson, Jeremy K.; Williamson, Tanja N.; Nelson, Hugh L.

2012-01-01

Flow- and load-duration curves were constructed from the model outputs of the U.S. Geological Survey's Water Availability Tool for Environmental Resources (WATER) application for streams in Kentucky. The WATER application was designed to access multiple geospatial datasets to generate more than 60 years of statistically based streamflow data for Kentucky. The WATER application enables a user to graphically select a site on a stream and generate an estimated hydrograph and flow-duration curve for the watershed upstream of that point. The flow-duration curves are constructed by calculating the exceedance probability of the modeled daily streamflows. User-defined water-quality criteria and (or) sampling results can be loaded into the WATER application to construct load-duration curves that are based on the modeled streamflow results. Estimates of flow and streamflow statistics were derived from TOPographically Based Hydrological MODEL (TOPMODEL) simulations in the WATER application. A modified TOPMODEL code, SDP-TOPMODEL (Sinkhole Drainage Process-TOPMODEL) was used to simulate daily mean discharges over the period of record for 5 karst and 5 non-karst watersheds in Kentucky in order to verify the calibrated model. A statistical evaluation of the model's verification simulations show that calibration criteria, established by previous WATER application reports, were met thus insuring the model's ability to provide acceptably accurate estimates of discharge at gaged and ungaged sites throughout Kentucky. Flow-duration curves are constructed in the WATER application by calculating the exceedence probability of the modeled daily flow values. The flow-duration intervals are expressed as a percentage, with zero corresponding to the highest stream discharge in the streamflow record. Load-duration curves are constructed by applying the loading equation (Load = Flow*Water-quality criterion) at each flow interval.

Numerical model of water flow and solute accumulation in vertisols using HYDRUS 2D/3D code

NASA Astrophysics Data System (ADS)

Weiss, Tomáš; Dahan, Ofer; Turkeltub, Tuvia

2015-04-01

Keywords: dessication-crack-induced-salinization, preferential flow, conceptual model, numerical model, vadose zone, vertisols, soil water retention function, HYDRUS 2D/3D Vertisols cover a hydrologically very significant area of semi-arid regions often through which water infiltrates to groundwater aquifers. Understanding of water flow and solute accumulation is thus very relevant to agricultural activity and water resources management. Previous works suggest a conceptual model of dessication-crack-induced-salinization where salinization of sediment in the deep section of the vadose zone (up to 4 m) is induced by subsurface evaporation due to convective air flow in the dessication cracks. It suggests that the salinization is induced by the hydraulic gradient between the dry sediment in the vicinity of cracks (low potential) and the relatively wet sediment further from the main cracks (high potential). This paper presents a modified previously suggested conceptual model and a numerical model. The model uses a simple uniform flow approach but unconventionally prescribes the boundary conditions and the hydraulic parameters of soil. The numerical model is bound to one location close to a dairy farm waste lagoon, but the application of the suggested conceptual model could be possibly extended to all semi-arid regions with vertisols. Simulations were conducted using several modeling approaches with an ultimate goal of fitting the simulation results to the controlling variables measured in the field: temporal variation in water content across thick layer of unsaturated clay sediment (>10 m), sediment salinity and salinity the water draining down the vadose zone to the water table. The development of the model was engineered in several steps; all computed as forward solutions by try-and-error approach. The model suggests very deep instant infiltration of fresh water up to 12 m, which is also supported by the field data. The paper suggests prescribing a special atmospheric

Ground-water flow in the New Jersey Coastal Plain

USGS Publications Warehouse

Martin, Mary

1998-01-01

Ground-water flow in 10 aquifers and 9 intervening confining units of the New Jersey Coastal Plain was simulated as part of the Regional Aquifer System Analysis. Data on aquifer and confining unit characteristics and on pumpage and water levels from 1918 through 1980 were incorporated into a multilayer finite-difference model. The report describes the conceptual hydrogeologic model of the unstressed flow systems, the methods and approach used in simulating flow, and the results of the simulations.

Modeling of Water Flow Processes in the Soil-Plant-Atmosphere System: The Soil-Tree-Atmosphere Continuum Model

NASA Astrophysics Data System (ADS)

Massoud, E. C.; Vrugt, J. A.

2015-12-01

Trees and forests play a key role in controlling the water and energy balance at the land-air surface. This study reports on the calibration of an integrated soil-tree-atmosphere continuum (STAC) model using Bayesian inference with the DREAM algorithm and temporal observations of soil moisture content, matric head, sap flux, and leaf water potential from the King's River Experimental Watershed (KREW) in the southern Sierra Nevada mountain range in California. Water flow through the coupled system is described using the Richards' equation with both the soil and tree modeled as a porous medium with nonlinear soil and tree water relationships. Most of the model parameters appear to be reasonably well defined by calibration against the observed data. The posterior mean simulation reproduces the observed soil and tree data quite accurately, but a systematic mismatch is observed between early afternoon measured and simulated sap fluxes. We will show how this points to a structural error in the STAC-model and suggest and test an alternative hypothesis for root water uptake that alleviates this problem.

Simulation of ground-water flow and evaluation of water-management alternatives in the upper Charles River basin, eastern Massachusetts

USGS Publications Warehouse

DeSimone, Leslie A.; Walter, Donald A.; Eggleston, John R.; Nimiroski, Mark T.

2002-01-01

Ground water is the primary source of drinking water for towns in the upper Charles River Basin, an area of 105 square miles in eastern Massachusetts that is undergoing rapid growth. The stratified-glacial aquifers in the basin are high yield, but also are thin, discontinuous, and in close hydraulic connection with streams, ponds, and wetlands. Water withdrawals averaged 10.1 million gallons per day in 1989?98 and are likely to increase in response to rapid growth. These withdrawals deplete streamflow and lower pond levels. A study was conducted to develop tools for evaluating water-management alternatives at the regional scale in the basin. Geologic and hydrologic data were compiled and collected to characterize the ground- and surface-water systems. Numerical flow modeling techniques were applied to evaluate the effects of increased withdrawals and altered recharge on ground-water levels, pond levels, and stream base flow. Simulation-optimization methods also were applied to test their efficacy for management of multiple water-supply and water-resource needs. Steady-state and transient ground-water-flow models were developed using the numerical modeling code MODFLOW-2000. The models were calibrated to 1989?98 average annual conditions of water withdrawals, water levels, and stream base flow. Model recharge rates were varied spatially, by land use, surficial geology, and septic-tank return flow. Recharge was changed during model calibration by means of parameter-estimation techniques to better match the estimated average annual base flow; area-weighted rates averaged 22.5 inches per year for the basin. Water withdrawals accounted for about 7 percent of total simulated flows through the stream-aquifer system and were about equal in magnitude to model-calculated rates of ground-water evapotranspiration from wetlands and ponds in aquifer areas. Water withdrawals as percentages of total flow varied spatially and temporally within an average year; maximum values were

An electrode polarization impedance based flow sensor for low water flow measurement

NASA Astrophysics Data System (ADS)

Yan, Tinghu; Sabic, Darko

2013-06-01

This note describes an electrode polarization impedance based flow sensor for low water flow measurement. It consists of two pairs of stainless steel electrodes set apart and inserted into a non-conductive flow tube with each pair of electrodes placed diametrically at the opposite sides. The flow sensor is modeled as a typical four-electrode system of which two electrodes are current-carrying and the other two serve as output pick ups. The polarization impedances of the two current carrying electrodes are affected by water flows resulting in changes of differential potential between the two pick-up electrodes which are separated by the same fluid. The interrogation of the two excitation electrodes with dc biased ac signals offers significantly higher sensor sensitivities to flow. The prototype flow sensor constructed for a 20 mm diameter pipeline was able to measure water flow rate as low as tested at 1.06 l h-1 and remained sensitive at a flow rate of 25.18 l h-1 when it was driven with a sinusoidal voltage at 1000 Hz with a peak ac amplitude of 2 V and a dc offset of +8 V. The nonlinear characteristics of the sensor response indicate that the sensor is more sensitive at low flows and will not be able to measure at very high flows. Additional experiments are needed to evaluate the influences of impurities, chemical species, ions constituents, conductivity and temperature over a practical range of residential water conditions, the effects of fluctuating ground signals, measurement uncertainty, power consumption, compensation of effects and practical operations. The flow sensor (principle) presented may be used as (in) a secondary sensor in combination with an existing electronic water meter to extend the low end of measurement range in residential water metering.

Simulation of ground-water flow and land subsidence in the Antelope Valley ground-water basin, California

USGS Publications Warehouse

Leighton, David A.; Phillips, Steven P.

2003-01-01

Antelope Valley, California, is a topographically closed basin in the western part of the Mojave Desert, about 50 miles northeast of Los Angeles. The Antelope Valley ground-water basin is about 940 square miles and is separated from the northern part of Antelope Valley by faults and low-lying hills. Prior to 1972, ground water provided more than 90 percent of the total water supply in the valley; since 1972, it has provided between 50 and 90 percent. Most ground-water pumping in the valley occurs in the Antelope Valley ground-water basin, which includes the rapidly growing cities of Lancaster and Palmdale. Ground-water-level declines of more than 200 feet in some parts of the ground-water basin have resulted in an increase in pumping lifts, reduced well efficiency, and land subsidence of more than 6 feet in some areas. Future urban growth and limits on the supply of imported water may continue to increase reliance on ground water. To better understand the ground-water flow system and to develop a tool to aid in effectively managing the water resources, a numerical model of ground-water flow and land subsidence in the Antelope Valley ground-water basin was developed using old and new geohydrologic information. The ground-water flow system consists of three aquifers: the upper, middle, and lower aquifers. The aquifers, which were identified on the basis of the hydrologic properties, age, and depth of the unconsolidated deposits, consist of gravel, sand, silt, and clay alluvial deposits and clay and silty clay lacustrine deposits. Prior to ground-water development in the valley, recharge was primarily the infiltration of runoff from the surrounding mountains. Ground water flowed from the recharge areas to discharge areas around the playas where it discharged either from the aquifer system as evapotranspiration or from springs. Partial barriers to horizontal ground-water flow, such as faults, have been identified in the ground-water basin. Water-level declines owing to

Patterns and age distribution of ground-water flow to streams

USGS Publications Warehouse

Modica, E.; Reilly, T.E.; Pollock, D.W.

1997-01-01

Simulations of ground-water flow in a generic aquifer system were made to characterize the topology of ground-water flow in the stream subsystem and to evaluate its relation to deeper ground-water flow. The flow models are patterned after hydraulic characteristics of aquifers of the Atlantic Coastal Plain and are based on numerical solutions to three-dimensional, steady-state, unconfined flow. The models were used to evaluate the effects of aquifer horizontal-to-vertical hydraulic conductivity ratios, aquifer thickness, and areal recharge rates on flow in the stream subsystem. A particle tracker was used to determine flow paths in a stream subsystem, to establish the relation between ground-water seepage to points along a simulated stream and its source area of flow, and to determine ground-water residence time in stream subsystems. In a geometrically simple aquifer system with accretion, the source area of flow to streams resembles an elongated ellipse that tapers in the downgradient direction. Increased recharge causes an expansion of the stream subsystem. The source area of flow to the stream expands predominantly toward the stream headwaters. Baseflow gain is also increased along the reach of the stream. A thin aquifer restricts ground-water flow and causes the source area of flow to expand near stream headwaters and also shifts the start-of-flow to the drainage basin divide. Increased aquifer anisotropy causes a lateral expansion of the source area of flow to streams. Ground-water seepage to the stream channel originates both from near- and far-recharge locations. The range in the lengths of flow paths that terminate at a point on a stream increase in the downstream direction. Consequently, the age distribution of ground water that seeps into the stream is skewed progressively older with distance downstream. Base flow ia an integration of ground water with varying age and potentially different water quality, depending on the source within the drainage basin

Water tunnel flow visualization using a laser

NASA Technical Reports Server (NTRS)

Beckner, C.; Curry, R. E.

1985-01-01

Laser systems for flow visualization in water tunnels (similar to the vapor screen technique used in wind tunnels) can provide two-dimensional cross-sectional views of complex flow fields. This parametric study documents the practical application of the laser-enhanced visualization (LEV) technique to water tunnel testing. Aspects of the study include laser power levels, flow seeding (using flourescent dyes and embedded particulates), model preparation, and photographic techniques. The results of this study are discussed to provide potential users with basic information to aid in the design and setup of an LEV system.

Regional groundwater-flow model of the Lake Michigan Basin in support of Great Lakes Basin water availability and use studies

USGS Publications Warehouse

Feinstein, D.T.; Hunt, R.J.; Reeves, H.W.

2010-01-01

A regional groundwater-flow model of the Lake Michigan Basin and surrounding areas has been developed in support of the Great Lakes Basin Pilot project under the U.S. Geological Survey's National Water Availability and Use Program. The transient 2-million-cell model incorporates multiple aquifers and pumping centers that create water-level drawdown that extends into deep saline waters. The 20-layer model simulates the exchange between a dense surface-water network and heterogeneous glacial deposits overlying stratified bedrock of the Wisconsin/Kankakee Arches and Michigan Basin in the Lower and Upper Peninsulas of Michigan; eastern Wisconsin; northern Indiana; and northeastern Illinois. The model is used to quantify changes in the groundwater system in response to pumping and variations in recharge from 1864 to 2005. Model results quantify the sources of water to major pumping centers, illustrate the dynamics of the groundwater system, and yield measures of water availability useful for water-resources management in the region. This report is a complete description of the methods and datasets used to develop the regional model, the underlying conceptual model, and model inputs, including specified values of material properties and the assignment of external and internal boundary conditions. The report also documents the application of the SEAWAT-2000 program for variable-density flow; it details the approach, advanced methods, and results associated with calibration through nonlinear regression using the PEST program; presents the water-level, drawdown, and groundwater flows for various geographic subregions and aquifer systems; and provides analyses of the effects of pumping from shallow and deep wells on sources of water to wells, the migration of groundwater divides, and direct and indirect groundwater discharge to Lake Michigan. The report considers the role of unconfined conditions at the regional scale as well as the influence of salinity on groundwater flow

Simulating Water Flow in Variably Saturated Soils - Exploring the Advantage of Three-dimensional Models

NASA Astrophysics Data System (ADS)

Hopp, L.; Ivanov, V. Y.

2010-12-01

There is still a debate in rainfall-runoff modeling over the advantage of using three-dimensional models based on partial differential equations describing variably saturated flow vs. models with simpler infiltration and flow routing algorithms. Fully explicit 3D models are computationally demanding but allow the representation of spatially complex domains, heterogeneous soils, conditions of ponded infiltration, and solute transport, among others. Models with simpler infiltration and flow routing algorithms provide faster run times and are likely to be more versatile in the treatment of extreme conditions such as soil drying but suffer from underlying assumptions and ad-hoc parameterizations. In this numerical study, we explore the question of whether these two model strategies are competing approaches or if they complement each other. As a 3D physics-based model we use HYDRUS-3D, a finite element model that numerically solves the Richards equation for variably-saturated water flow. As an example of a simpler model, we use tRIBS+VEGGIE that solves the 1D Richards equation for vertical flow and applies Dupuit-Forchheimer approximation for saturated lateral exchange and gravity-driven flow for unsaturated lateral exchange. The flow can be routed using either the D-8 (steepest descent) or D-infinity flow routing algorithms. We study lateral subsurface stormflow and moisture dynamics at the hillslope-scale, using a zero-order basin topography, as a function of storm size, antecedent moisture conditions and slope angle. The domain and soil characteristics are representative of a forested hillslope with conductive soils in a humid environment, where the major runoff generating process is lateral subsurface stormflow. We compare spatially integrated lateral subsurface flow at the downslope boundary as well as spatial patterns of soil moisture. We illustrate situations where both model approaches perform equally well and identify conditions under which the application of a

Estimating water flow through a hillslope using the massively parallel processor

NASA Technical Reports Server (NTRS)

Devaney, Judy E.; Camillo, P. J.; Gurney, R. J.

1988-01-01

A new two-dimensional model of water flow in a hillslope has been implemented on the Massively Parallel Processor at the Goddard Space Flight Center. Flow in the soil both in the saturated and unsaturated zones, evaporation and overland flow are all modelled, and the rainfall rates are allowed to vary spatially. Previous models of this type had always been very limited computationally. This model takes less than a minute to model all the components of the hillslope water flow for a day. The model can now be used in sensitivity studies to specify which measurements should be taken and how accurate they should be to describe such flows for environmental studies.

A framework for modeling connections between hydraulics, water surface roughness, and surface reflectance in open channel flows

USGS Publications Warehouse

Legleiter, Carl; Mobley, Curtis D.; Overstreet, Brandon

2017-01-01

This paper introduces a framework for examining connections between the flow field, the texture of the air-water interface, and the reflectance of the water surface and thus evaluating the potential to infer hydraulic information from remotely sensed observations of surface reflectance. We used a spatial correlation model describing water surface topography to illustrate the application of our framework. Nondimensional relations between model parameters and flow intensity were established based on a prior flume study. Expressing the model in the spatial frequency domain allowed us to use an efficient Fourier transform-based algorithm for simulating water surfaces. Realizations for both flume and field settings had water surface slope distributions positively correlated with velocity and water surface roughness. However, most surface facets were gently sloped and thus unlikely to yield strong specular reflections; the model exaggerated the extent of water surface features, leading to underestimation of facet slopes. A ray tracing algorithm indicated that reflectance was greatest when solar and view zenith angles were equal and the sensor scanned toward the Sun to capture specular reflections of the solar beam. Reflected energy was concentrated in a small portion of the sky, but rougher water surfaces reflected rays into a broader range of directions. Our framework facilitates flight planning to avoid surface-reflected radiance while mapping other river attributes, or to maximize this component to exploit relationships between hydraulics and surface reflectance. This initial analysis also highlighted the need for improved models of water surface topography in natural rivers.

Instability of water-ice interface under turbulent flow

NASA Astrophysics Data System (ADS)

Izumi, Norihiro; Naito, Kensuke; Yokokawa, Miwa

2015-04-01

It is known that plane water-ice interface becomes unstable to evolve into a train of waves. The underside of ice formed on the water surface of rivers are often observed to be covered with ice ripples. Relatively steep channels which discharge melting water from glaciers are characterized by beds covered with a series of steps. Though the flowing agent inducing instability is not water but gas including water vapor, a similar train of steps have been recently observed on the Polar Ice Caps on Mars (Spiral Troughs). They are expected to be caused by the instability of water-ice interface induced by flowing fluid on ice. There have been some studies on this instability in terms of linear stability analysis. Recently, Caporeale and Ridolfi (2012) have proposed a complete linear stability analysis in the case of laminar flow, and found that plane water-ice interface is unstable in the range of sufficiently large Reynolds numbers, and that the important parameters are the Reynolds number, the slope angle, and the water surface temperature. However, the flow inducing instability on water-ice interface in the field should be in the turbulent regime. Extension of the analysis to the case of fully developed turbulent flow with larger Reynolds numbers is needed. We have performed a linear stability analysis on the instability of water-ice interface under turbulent flow conditions with the use of the Reynolds-averaged Navier-Stokes equations with the mixing length turbulent model, the continuity equation of flow, the diffusion/dispersion equation of heat, and the Stefan equation. In order to reproduce the accurate velocity distribution and the heat transfer in the vicinity of smooth walls with the use of the mixing length model, it is important to take into account of the rapid decrease in the mixing length in the viscous sublayer. We employ the Driest model (1956) to the formulation. In addition, as the thermal boundary condition at the water surface, we describe the

The importance of base flow in sustaining surface water flow in the Upper Colorado River Basin

USGS Publications Warehouse

Miller, Matthew P.; Buto, Susan G.; Susong, David D.; Rumsey, Christine

2016-01-01

The Colorado River has been identified as the most overallocated river in the world. Considering predicted future imbalances between water supply and demand and the growing recognition that base flow (a proxy for groundwater discharge to streams) is critical for sustaining flow in streams and rivers, there is a need to develop methods to better quantify present-day base flow across large regions. We adapted and applied the spatially referenced regression on watershed attributes (SPARROW) water quality model to assess the spatial distribution of base flow, the fraction of streamflow supported by base flow, and estimates of and potential processes contributing to the amount of base flow that is lost during in-stream transport in the Upper Colorado River Basin (UCRB). On average, 56% of the streamflow in the UCRB originated as base flow, and precipitation was identified as the dominant driver of spatial variability in base flow at the scale of the UCRB, with the majority of base flow discharge to streams occurring in upper elevation watersheds. The model estimates an average of 1.8 × 1010 m3/yr of base flow in the UCRB; greater than 80% of which is lost during in-stream transport to the Lower Colorado River Basin via processes including evapotranspiration and water diversion for irrigation. Our results indicate that surface waters in the Colorado River Basin are dependent on base flow, and that management approaches that consider groundwater and surface water as a joint resource will be needed to effectively manage current and future water resources in the Basin.

The importance of base flow in sustaining surface water flow in the Upper Colorado River Basin

NASA Astrophysics Data System (ADS)

Miller, Matthew P.; Buto, Susan G.; Susong, David D.; Rumsey, Christine A.

2016-05-01

The Colorado River has been identified as the most overallocated river in the world. Considering predicted future imbalances between water supply and demand and the growing recognition that base flow (a proxy for groundwater discharge to streams) is critical for sustaining flow in streams and rivers, there is a need to develop methods to better quantify present-day base flow across large regions. We adapted and applied the spatially referenced regression on watershed attributes (SPARROW) water quality model to assess the spatial distribution of base flow, the fraction of streamflow supported by base flow, and estimates of and potential processes contributing to the amount of base flow that is lost during in-stream transport in the Upper Colorado River Basin (UCRB). On average, 56% of the streamflow in the UCRB originated as base flow, and precipitation was identified as the dominant driver of spatial variability in base flow at the scale of the UCRB, with the majority of base flow discharge to streams occurring in upper elevation watersheds. The model estimates an average of 1.8 × 1010 m3/yr of base flow in the UCRB; greater than 80% of which is lost during in-stream transport to the Lower Colorado River Basin via processes including evapotranspiration and water diversion for irrigation. Our results indicate that surface waters in the Colorado River Basin are dependent on base flow, and that management approaches that consider groundwater and surface water as a joint resource will be needed to effectively manage current and future water resources in the Basin.

Multiphase flow modelling of volcanic ash particle settling in water using adaptive unstructured meshes

NASA Astrophysics Data System (ADS)

Jacobs, C. T.; Collins, G. S.; Piggott, M. D.; Kramer, S. C.; Wilson, C. R. G.

2013-02-01

Small-scale experiments of volcanic ash particle settling in water have demonstrated that ash particles can either settle slowly and individually, or rapidly and collectively as a gravitationally unstable ash-laden plume. This has important implications for the emplacement of tephra deposits on the seabed. Numerical modelling has the potential to extend the results of laboratory experiments to larger scales and explore the conditions under which plumes may form and persist, but many existing models are computationally restricted by the fixed mesh approaches that they employ. In contrast, this paper presents a new multiphase flow model that uses an adaptive unstructured mesh approach. As a simulation progresses, the mesh is optimized to focus numerical resolution in areas important to the dynamics and decrease it where it is not needed, thereby potentially reducing computational requirements. Model verification is performed using the method of manufactured solutions, which shows the correct solution convergence rates. Model validation and application considers 2-D simulations of plume formation in a water tank which replicate published laboratory experiments. The numerically predicted settling velocities for both individual particles and plumes, as well as instability behaviour, agree well with experimental data and observations. Plume settling is clearly hindered by the presence of a salinity gradient, and its influence must therefore be taken into account when considering particles in bodies of saline water. Furthermore, individual particles settle in the laminar flow regime while plume settling is shown (by plume Reynolds numbers greater than unity) to be in the turbulent flow regime, which has a significant impact on entrainment and settling rates. Mesh adaptivity maintains solution accuracy while providing a substantial reduction in computational requirements when compared to the same simulation performed using a fixed mesh, highlighting the benefits of an

Modeling axisymmetric flow and transport

USGS Publications Warehouse

Langevin, C.D.

2008-01-01

Unmodified versions of common computer programs such as MODFLOW, MT3DMS, and SEAWAT that use Cartesian geometry can accurately simulate axially symmetric ground water flow and solute transport. Axisymmetric flow and transport are simulated by adjusting several input parameters to account for the increase in flow area with radial distance from the injection or extraction well. Logarithmic weighting of interblock transmissivity, a standard option in MODFLOW, can be used for axisymmetric models to represent the linear change in hydraulic conductance within a single finite-difference cell. Results from three test problems (ground water extraction, an aquifer push-pull test, and upconing of saline water into an extraction well) show good agreement with analytical solutions or with results from other numerical models designed specifically to simulate the axisymmetric geometry. Axisymmetric models are not commonly used but can offer an efficient alternative to full three-dimensional models, provided the assumption of axial symmetry can be justified. For the upconing problem, the axisymmetric model was more than 1000 times faster than an equivalent three-dimensional model. Computational gains with the axisymmetric models may be useful for quickly determining appropriate levels of grid resolution for three-dimensional models and for estimating aquifer parameters from field tests.

Modeling of ground-water flow in subsurface Austin Chalk and Taylor marl in Ellis County, Texas, near the superconducting super collider site

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

Mace, R.E.

1993-02-01

Numerical models are useful tools for developing an understanding of ground-water flow in sparsely characterized low-permeability aquifers. Finite-difference, cross-sectional models of Cretaceous chalk and marl formations near the Superconducting Super Collider (SSC) were constructed using MODFLOW to evaluate ground-water circulation paths and travel times. Weathered and fractured zones with enhanced permeability were included to assess the effect these features had on flow paths and times. Pump tests, slug tests, packer tests, core tests, and estimates were used to define hydraulic properties for model input. The model was calibrated with water-level data from monitor wells and from wire-line piezometers near amore » test shaft excavated by the SSC project. A ratio of vertical-to-horizontal permeability of 0.0085 was estimated through model calibration. A chalk-to-marl permeability ratio of 18 was needed to reproduce artesian head in a well completed in chalk beneath marl. Hydraulic head distributions and ground-water flow paths reflected local, intermediate, and regional flow systems with recharge beneath upland surface-water divides and discharge in valleys. Most of the flow (99%) occurred in the weathered zone, with average residence times of 5 to 10 years. Residence time in unweathered chalk bedrock was substantially longer, at an average of 1.7 Ma. As expected, the model demonstrated that deep and rapid ground-water circulation might occur in fracture zones. Particle paths calculated using MODPATH showed that ground-water travel times from recharge areas to the SSC subsurface facilities might be 20 to 60 years where flow is through fracture zones.« less

Ground-water flow model of the Boone formation at the Tar Creek superfund site, Oklahoma and Kansas

USGS Publications Warehouse

Reed, T.B.; Czarnecki, John B.

2006-01-01

Extensive mining activities conducted at the Tar Creek Superfund site, one of the largest Superfund sites in the United States, pose substantial health and safety risks. Mining activities removed a total of about 6,000,000 tons of lead and zinc by 1949. To evaluate the effect of this mining on the ground-water flow, a MODFLOW 2000 digital model has been developed to simulate ground-water flow in the carbonate formations of Mississippian age underlying the Tar Creek Superfund site. The model consists of three layers of variable thickness and a grid of 580 rows by 680 columns of cells 164 feet (50 meters) on a side. Model flux boundary conditions are specified for rivers and general head boundaries along the northern boundary of the Boone Formation. Selected cells in layer 1 are simulated as drain cells. Model calibration has been performed to minimize the difference between simulated and observed water levels in the Boone Formation. Hydraulic conductivity values specified during calibration range from 1.3 to 35 feet per day for the Boone Formation with the larger values occurring along the axis of the Miami Syncline where horizontal anisotropy is specified as 10 to 1. Hydraulic conductivity associated with the mine void is set at 50,000 feet per day and a specific yield of 1.0 is specified to represent that the mine void is filled completely with water. Residuals (the difference between measured and simulated ground-water altitudes) has a root-mean-squared value of 8.53 feet and an absolute mean value of 7.29 feet for 17 observed values of water levels in the Boone Formation. The utility of the model for simulating and evaluating the possible consequences of remediation activities has been demonstrated. The model was used to simulate the emplacement of chat (mine waste consisting of fines and fragments of chert) back into the mine. Scenarios using 1,800,000 and 6,500,000 tons of chat were run. Hydraulic conductivity was reduced from 50,000 feet per day to 35 feet

Marangoni flow in an evaporating water droplet

NASA Astrophysics Data System (ADS)

Xu, Xuefeng; Luo, Jianbin

2007-09-01

Marangoni effect has been observed in many liquids, but its existence in pure water is still a debated problem. In the present work, the Marangoni flow in evaporating water droplets has been observed by using fluorescent nanoparticles. Flow patterns indicate that a stagnation point where the surface flow, the surface tension gradient, and the surface temperature gradient change their directions exists at the droplet surface. The deduced nonmonotonic variation of the droplet surface temperature, which is different from that in some previous works, is explained by a heat transfer model considering the adsorbed thin film of the evaporating liquid droplet.

PHYSICAL MODELING OF CONTRACTED FLOW.

USGS Publications Warehouse

Lee, Jonathan K.

1987-01-01

Experiments on steady flow over uniform grass roughness through centered single-opening contractions were conducted in the Flood Plain Simulation Facility at the U. S. Geological Survey's Gulf Coast Hydroscience Center near Bay St. Louis, Miss. The experimental series was designed to provide data for calibrating and verifying two-dimensional, vertically averaged surface-water flow models used to simulate flow through openings in highway embankments across inundated flood plains. Water-surface elevations, point velocities, and vertical velocity profiles were obtained at selected locations for design discharges ranging from 50 to 210 cfs. Examples of observed water-surface elevations and velocity magnitudes at basin cross-sections are presented.

Simulation of ground-water/surface-water flow in the Santa Clara-Calleguas ground-water basin, Ventura County, California

USGS Publications Warehouse

Hanson, Randall T.; Martin, Peter; Koczot, Kathryn M.

2003-01-01

Ground water is the main source of water in the Santa Clara-Calleguas ground-water basin that covers about 310 square miles in Ventura County, California. A steady increase in the demand for surface- and ground-water resources since the late 1800s has resulted in streamflow depletion and ground-water overdraft. This steady increase in water use has resulted in seawater intrusion, inter-aquifer flow, land subsidence, and ground-water contamination. The Santa Clara-Calleguas Basin consists of multiple aquifers that are grouped into upper- and lower-aquifer systems. The upper-aquifer system includes the Shallow, Oxnard, and Mugu aquifers. The lower-aquifer system includes the upper and lower Hueneme, Fox Canyon, and Grimes Canyon aquifers. The layered aquifer systems are each bounded below by regional unconformities that are overlain by extensive basal coarse-grained layers that are the major pathways for ground-water production from wells and related seawater intrusion. The aquifer systems are bounded below and along mountain fronts by consolidated bedrock that forms a relatively impermeable boundary to ground-water flow. Numerous faults act as additional exterior and interior boundaries to ground-water flow. The aquifer systems extend offshore where they crop out along the edge of the submarine shelf and within the coastal submarine canyons. Submarine canyons have dissected these regional aquifers, providing a hydraulic connection to the ocean through the submarine outcrops of the aquifer systems. Coastal landward flow (seawater intrusion) occurs within both the upper- and lower-aquifer systems. A numerical ground-water flow model of the Santa Clara-Calleguas Basin was developed by the U.S. Geological Survey to better define the geohydrologic framework of the regional ground-water flow system and to help analyze the major problems affecting water-resources management of a typical coastal aquifer system. Construction of the Santa Clara-Calleguas Basin model required

Flow in water-intake pump bays: A guide for utility engineers. Final report

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

Ettema, R.

1998-09-01

This report is intended to serve as a guide for power-plant engineers facing problems with flow conditions in pump bays in water-intake structures, especially those located alongside rivers. The guide briefly introduces the typical prevailing flow field outside of a riverside water intake. That flow field often sets the inflow conditions for pump bays located within the water intake. The monograph then presents and discusses the main flow problems associated with pump bays. The problems usually revolve around the formation of troublesome vortices. A novel feature of this monograph is the use of numerical modeling to reveal diagnostically how themore » vortices form and their sensitivities to flow conditions, such as uniformity of approach flow entering the bay and water-surface elevation relative to pump-bell submergence. The modeling was carried out using a computer code developed specially for the present project. Pump-bay layouts are discussed next. The discussion begins with a summary of the main variables influencing bay flows. The numerical model is used to determine the sensitivities of the vortices to variations in the geometric parameters. The fixes include the use of flow-control vanes and suction scoops for ensuring satisfactory flow performance in severe flow conditions; notably flows with strong cross flow and shallow flows. The monograph ends with descriptions of modeling techniques. An extensive discussion is provided on the use of numerical model for illuminating bay flows. The model is used to show how fluid viscosity affects bay flow. The effect of fluid viscosity is an important consideration in hydraulic modeling of water intakes.« less

Evaluation of the ground-water flow model for northern Utah Valley, Utah, updated to conditions through 2002

USGS Publications Warehouse

Thiros, Susan A.

2006-01-01

This report evaluates the performance of a numerical model of the ground-water system in northern Utah Valley, Utah, that originally simulated ground-water conditions during 1947-1980 and was updated to include conditions estimated for 1981-2002. Estimates of annual recharge to the ground-water system and discharge from wells in the area were added to the original ground-water flow model of the area.The files used in the original transient-state model of the ground-water flow system in northern Utah Valley were imported into MODFLOW-96, an updated version of MODFLOW. The main model input files modified as part of this effort were the well and recharge files. Discharge from pumping wells in northern Utah Valley was estimated on an annual basis for 1981-2002. Although the amount of average annual withdrawals from wells has not changed much since the previous study, there have been changes in the distribution of well discharge in the area. Discharge estimates for flowing wells during 1981-2002 were assumed to be the same as those used in the last stress period of the original model because of a lack of new data. Variations in annual recharge were assumed to be proportional to changes in total surface-water inflow to northern Utah Valley. Recharge specified in the model during the additional stress periods varied from 255,000 acre-feet in 1986 to 137,000 acre-feet in 1992.The ability of the updated transient-state model to match hydrologic conditions determined for 1981-2002 was evaluated by comparing water-level changes measured in wells to those computed by the model. Water-level measurements made in February, March, or April were available for 39 wells in the modeled area during all or part of 1981-2003. In most cases, the magnitude and direction of annual water-level change from 1981 to 2002 simulated by the updated model reasonably matched the measured change. The greater-than-normal precipitation that occurred during 1982-84 resulted in period-of-record high

Accounting for intracell flow in models with emphasis on water table recharge and stream-aquifer interaction: 1. Problems and concepts

USGS Publications Warehouse

Jorgensen, Donald G.; Signor, Donald C.; Imes, Jeffrey L.

1989-01-01

Intracell flow is important in modeling cells that contain both sources and sinks. Special attention is needed if recharge through the water table is a source. One method of modeling multiple sources and sinks is to determine the net recharge per cell. For example, for a model cell containing both a sink and recharge through the water table, the amount of recharge should be reduced by the ratio of the area of influence of the sink within the cell to the area of the cell. The reduction is the intercepted portion of the recharge. In a multilayer model this amount is further reduced by a proportion factor, which is a function of the depth of the flow lines from the water table boundary to the internal sink. A gaining section of a stream is a typical sink. The aquifer contribution to a gaining stream can be conceptualized as having two parts; the first part is the intercepted lateral flow from the water table and the second is the flow across the streambed due to differences in head between the water level in the stream and the aquifer below. The amount intercepted is a function of the geometry of the cell, but the amount due to difference in head across the stream bed is largely independent of cell geometry. A discharging well can intercept recharge through the water table within a model cell. The net recharge to the cell would be reduced in proportion to the area of influence of the well within the cell. The area of influence generally changes with time. Thus the amount of intercepted recharge and net recharge may not be constant with time. During periods when the well is not discharging there will be no intercepted recharge even though the area of influence from previous pumping may still exist. The reduction of net recharge per cell due to internal interception of flow will result in a model-calculated mass balance less than the prototype. Additionally the “effective transmissivity” along the intercell flow paths may be altered when flow paths are occupied by

Patterns and rates of ground-water flow on Long Island, New York

USGS Publications Warehouse

Buxton, Herbert T.; Modica, Edward

1992-01-01

Increased ground-water contamination from human activities on Long Island has prompted studies to define the pattern and rate of ground-water movement. A two-dimensional, fine-mesh, finite-element model consisting of 11,969 nodes and 22,880 elements was constructed to represent ground-water flow along a north-south section through central Long Island. The model represents average hydrologic conditions within a corridor approximately 15 miles wide. The model solves discrete approximations of both the potential and stream functions. The resulting flownet depicts flow paths and defines the vertical distribution of flow within the section. Ground-water flow rates decrease with depth. Sixty-two percent of the water flows no deeper than the upper glacial (water-table) aquifer, 38 percent enters the underlying Magothy aquifer, and only 3.1 percent enters the Lloyd aquifer. The limiting streamlines for flow to the Magothy and Lloyd aquifers indicate that aquifer recharge areas are narrow east-west bands through the center of the island. The recharge area of the Magothy aquifer is only 5.4 miles wide; that of the Lloyd aquifer is less than 0.5 miles. The distribution of ground-water traveltime and a flownet are calculated from model results; both are useful in the investigation of contaminant transport or the chemical evolution of ground water within the flow system. A major discontinuity in traveltime occurs across the streamline which separates the flow subsystems of the two confined aquifers. Water that reaches the Lloyd aquifer attains traveltimes as high as 10,000 years, whereas water that has not penetrated deeper than the Magothy aquifer attains traveltimes of only 2,000 years. The finite-element approach used in this study is particularly suited to ground-water systems that have complex hydrostratigraphy and cross-sectional symmetry.

Hydraulic-property estimates for use with a transient ground-water flow model of the Death Valley regional ground-water flow system, Nevada and California

USGS Publications Warehouse

Belcher, Wayne R.; Elliott, Peggy E.; Geldon, Arthur L.

2001-01-01

The Death Valley regional ground-water flow system encompasses an area of about 43,500 square kilometers in southeastern California and southern Nevada, between latitudes 35? and 38?15' north and longitudes 115? and 117?45' west. The study area is underlain by Quaternary to Tertiary basin-fill sediments and mafic-lava flows; Tertiary volcanic, volcaniclastic, and sedimentary rocks; Tertiary to Jurassic granitic rocks; Triassic to Middle Proterozoic carbonate and clastic sedimentary rocks; and Early Proterozoic igneous and metamorphic rocks. The rock assemblage in the Death Valley region is extensively faulted as a result of several episodes of tectonic activity. This study is comprised of published and unpublished estimates of transmissivity, hydraulic conductivity, storage coefficient, and anisotropy ratios for hydrogeologic units within the Death Valley region study area. Hydrogeologic units previously proposed for the Death Valley regional transient ground-water flow model were recognized for the purpose of studying the distribution of hydraulic properties. Analyses of regression and covariance were used to assess if a relation existed between hydraulic conductivity and depth for most hydrogeologic units. Those analyses showed a weak, quantitatively indeterminate, relation between hydraulic conductivity and depth.

Well logging interpretation of production profile in horizontal oil-water two phase flow pipes

NASA Astrophysics Data System (ADS)

Zhai, Lu-Sheng; Jin, Ning-De; Gao, Zhong-Ke; Zheng, Xi-Ke

2012-03-01

Due to the complicated distribution of local velocity and local phase hold up along the radial direction of pipe in horizontal oil-water two phase flow, it is difficult to measure the total flow rate and phase volume fraction. In this study, we carried out dynamic experiment in horizontal oil-water two phases flow simulation well by using combination measurement system including turbine flowmeter with petal type concentrating diverter, conductance sensor and flowpassing capacitance sensor. According to the response resolution ability of the conductance and capacitance sensor in different range of total flow rate and water-cut, we use drift flux model and statistical model to predict the partial phase flow rate, respectively. The results indicate that the variable coefficient drift flux model can self-adaptively tone the model parameter according to the oil-water two phase flow characteristic, and the prediction result of partial phase flow rate of oil-water two phase flow is of high accuracy.

The Landlab v1.0 OverlandFlow component: a Python tool for computing shallow-water flow across watersheds

NASA Astrophysics Data System (ADS)

Adams, Jordan M.; Gasparini, Nicole M.; Hobley, Daniel E. J.; Tucker, Gregory E.; Hutton, Eric W. H.; Nudurupati, Sai S.; Istanbulluoglu, Erkan

2017-04-01

Representation of flowing water in landscape evolution models (LEMs) is often simplified compared to hydrodynamic models, as LEMs make assumptions reducing physical complexity in favor of computational efficiency. The Landlab modeling framework can be used to bridge the divide between complex runoff models and more traditional LEMs, creating a new type of framework not commonly used in the geomorphology or hydrology communities. Landlab is a Python-language library that includes tools and process components that can be used to create models of Earth-surface dynamics over a range of temporal and spatial scales. The Landlab OverlandFlow component is based on a simplified inertial approximation of the shallow water equations, following the solution of de Almeida et al.(2012). This explicit two-dimensional hydrodynamic algorithm simulates a flood wave across a model domain, where water discharge and flow depth are calculated at all locations within a structured (raster) grid. Here, we illustrate how the OverlandFlow component contained within Landlab can be applied as a simplified event-based runoff model and how to couple the runoff model with an incision model operating on decadal timescales. Examples of flow routing on both real and synthetic landscapes are shown. Hydrographs from a single storm at multiple locations in the Spring Creek watershed, Colorado, USA, are illustrated, along with a map of shear stress applied on the land surface by flowing water. The OverlandFlow component can also be coupled with the Landlab DetachmentLtdErosion component to illustrate how the non-steady flow routing regime impacts incision across a watershed. The hydrograph and incision results are compared to simulations driven by steady-state runoff. Results from the coupled runoff and incision model indicate that runoff dynamics can impact landscape relief and channel concavity, suggesting that, on landscape evolution timescales, the OverlandFlow model may lead to differences in

Experimental validation of a 2D overland flow model using high resolution water depth and velocity data

NASA Astrophysics Data System (ADS)

Cea, L.; Legout, C.; Darboux, F.; Esteves, M.; Nord, G.

2014-05-01

This paper presents a validation of a two-dimensional overland flow model using empirical laboratory data. Unlike previous publications in which model performance is evaluated as the ability to predict an outlet hydrograph, we use high resolution 2D water depth and velocity data to analyze to what degree the model is able to reproduce the spatial distribution of these variables. Several overland flow conditions over two impervious surfaces of the order of one square meter with different micro and macro-roughness characteristics are studied. The first surface is a simplified representation of a sinusoidal terrain with three crests and furrows, while the second one is a mould of a real agricultural seedbed terrain. We analyze four different bed friction parameterizations and we show that the performance of formulations which consider the transition between laminar, smooth turbulent and rough turbulent flow do not improve the results obtained with Manning or Keulegan formulas for rough turbulent flow. The simulations performed show that using Keulegan formula with a physically-based definition of the bed roughness coefficient, a two-dimensional shallow water model is able to reproduce satisfactorily the flow hydrodynamics. It is shown that, even if the resolution of the topography data and numerical mesh are high enough to include all the small scale features of the bed surface, the roughness coefficient must account for the macro-roughness characteristics of the terrain in order to correctly reproduce the flow hydrodynamics.

Oahu Groundwater Flow Model

DOE Data Explorer

Nicole Lautze

2015-01-01

Groundwater flow model for the island of Oahu. Data is from the following sources: Rotzoll, K., A.I. El-Kadi. 2007. Numerical Ground-Water Flow Simulation for Red Hill Fuel Storage Facilities, NAVFAC Pacific, Oahu, Hawaii - Prepared TEC, Inc. Water Resources Research Center, University of Hawaii, Honolulu.; Whittier, R.B., K. Rotzoll, S. Dhal, A.I. El-Kadi, C. Ray, G. Chen, and D. Chang. 2004. Hawaii Source Water Assessment Program Report – Volume VII – Island of Oahu Source Water Assessment Program Report. Prepared for the Hawaii Department of Health, Safe Drinking Water Branch. University of Hawaii, Water Resources Research Center. Updated 2008.; and Whittier, R. and A.I. El-Kadi. 2009. Human and Environmental Risk Ranking of Onsite Sewage Disposal Systems – Final. Prepared by the University of Hawaii, Dept. of Geology and Geophysics for the State of Hawaii Dept. of Health, Safe Drinking Water Branch. December 2009.

Evaluation of a watershed model for estimating daily flow using limited flow measurements

USDA-ARS?s Scientific Manuscript database

The Soil and Water Assessment Tool (SWAT) model was evaluated for estimation of continuous daily flow based on limited flow measurements in the Upper Oyster Creek (UOC) watershed. SWAT was calibrated against limited measured flow data and then validated. The Nash-Sutcliffe model Efficiency (NSE) and...

Hydrogeology, simulated ground-water flow, and ground-water quality, Wright-Patterson Air Force Base, Ohio

USGS Publications Warehouse

Dumouchelle, D.H.; Schalk, C.W.; Rowe, G.L.; De Roche, J.T.

1993-01-01

Ground water is the primary source of water in the Wright-Patterson Air Force Base area. The aquifer consists of glacial sands and gravels that fill a buried bedrock-valley system. Consolidated rocks in the area consist of poorly permeable Ordovician shale of the Richmondian stage, in the upland areas, the Brassfield Limestone of Silurian age. The valleys are filled with glacial sediments of Wisconsinan age consisting of clay-rich tills and coarse-grained outwash deposits. Estimates of hydraulic conductivity of the shales based on results of displacement/recovery tests range from 0.0016 to 12 feet per day; estimates for the glacial sediments range from less than 1 foot per day to more than 1,000 feet per day. Ground water flow from the uplands towards the valleys and the major rivers in the region, the Great Miami and the Mad Rivers. Hydraulic-head data indicate that ground water flows between the bedrock and unconsolidated deposits. Data from a gain/loss study of the Mad River System and hydrographs from nearby wells reveal that the reach of the river next to Wright-Patterson Air Force Base is a ground-water discharge area. A steady-state, three-dimensional ground-water-flow model was developed to simulate ground-water flow in the region. The model contains three layers and encompasses about 100 square miles centered on Wright-Patterson Air Force Base. Ground water enters the modeled area primarily by river leakage and underflow at the model boundary. Ground water exits the modeled area primarily by flow through the valleys at the model boundaries and through production wells. A model sensitivity analysis involving systematic changes in values of hydrologic parameters in the model indicates that the model is most sensitive to decreases in riverbed conductance and vertical conductance between the upper two layers. The analysis also indicates that the contribution of water to the buried-valley aquifer from the bedrock that forms the valley walls is about 2 to 4

Efficacy of bedrock erosion by subglacial water flow

NASA Astrophysics Data System (ADS)

Beaud, F.; Flowers, G. E.; Venditti, J. G.

2015-09-01

Bedrock erosion by sediment-bearing subglacial water remains little-studied, however the process is thought to contribute to bedrock erosion rates in glaciated landscapes and is implicated in the excavation of tunnel valleys and the incision of inner gorges. We adapt physics-based models of fluvial abrasion to the subglacial environment, assembling the first model designed to quantify bedrock erosion caused by transient subglacial water flow. The subglacial drainage model consists of a one-dimensional network of cavities dynamically coupled to one or several Röthlisberger channels (R-channels). The bedrock erosion model is based on the tools and cover effect, whereby particles entrained by the flow impact exposed bedrock. We explore the dependency of glacial meltwater erosion on the structure and magnitude of water input to the system, the ice geometry and the sediment supply. We find that erosion is not a function of water discharge alone, but also depends on channel size, water pressure and on sediment supply, as in fluvial systems. Modelled glacial meltwater erosion rates are one to two orders of magnitude lower than the expected rates of total glacial erosion required to produce the sediment supply rates we impose, suggesting that glacial meltwater erosion is negligible at the basin scale. Nevertheless, due to the extreme localization of glacial meltwater erosion (at the base of R-channels), this process can carve bedrock (Nye) channels. In fact, our simulations suggest that the incision of bedrock channels several centimetres deep and a few meters wide can occur in a single year. Modelled incision rates indicate that subglacial water flow can gradually carve a tunnel valley and enhance the relief or even initiate the carving of an inner gorge.

MODELING TO EVOLVE UNDERSTANDING OF THE SHALLOW GROUND WATER FLOW SYSTEM BENEATH THE LIZZIE RESEARCH SITE, NC

EPA Science Inventory

The purpose of the modeling effort presented here is to evolve a conceptual model of ground-water flow at the Lizzie, NC research site using analytic solutions and field observations. The resulting analytic element parameterization of boundary conditions, aquifer transmissivitie...

Computer program for simulation of variable recharge with the U. S. Geological Survey modular finite-difference ground-water flow model (MODFLOW)

USGS Publications Warehouse

Kontis, A.L.

2001-01-01

The Variable-Recharge Package is a computerized method designed for use with the U.S. Geological Survey three-dimensional finitedifference ground-water flow model (MODFLOW-88) to simulate areal recharge to an aquifer. It is suitable for simulations of aquifers in which the relation between ground-water levels and land surface can affect the amount and distribution of recharge. The method is based on the premise that recharge to an aquifer cannot occur where the water level is at or above land surface. Consequently, recharge will vary spatially in simulations in which the Variable- Recharge Package is applied, if the water levels are sufficiently high. The input data required by the program for each model cell that can potentially receive recharge includes the average land-surface elevation and a quantity termed ?water available for recharge,? which is equal to precipitation minus evapotranspiration. The Variable-Recharge Package also can be used to simulate recharge to a valley-fill aquifer in which the valley fill and the adjoining uplands are explicitly simulated. Valley-fill aquifers, which are the most common type of aquifer in the glaciated northeastern United States, receive much of their recharge from upland sources as channeled and(or) unchanneled surface runoff and as lateral ground-water flow. Surface runoff in the uplands is generated in the model when the applied water available for recharge is rejected because simulated water levels are at or above land surface. The surface runoff can be distributed to other parts of the model by (1) applying the amount of the surface runoff that flows to upland streams (channeled runoff) to explicitly simulated streams that flow onto the valley floor, and(or) (2) applying the amount that flows downslope toward the valley- fill aquifer (unchanneled runoff) to specified model cells, typically those near the valley wall. An example model of an idealized valley- fill aquifer is presented to demonstrate application of the

Freshwater Ecosystem Service Flow Model To Evaluate Regional Water Security: A Case Study In Beijing-Tianjin-Hebei Region, China

NASA Astrophysics Data System (ADS)

Li, D.; Li, S.

2016-12-01

Freshwater service, as the most important support ecosystem service, is essential to human survival and development. Many studies have evidenced the spatial differences in the supply and demand of ecosystem services and raised the concept of ecosystem service flow. However, rather few studies quantitatively characterize the freshwater service flow. This paper aims to quantify the effect of freshwater ecosystem service flow on downstream areas in Beijing-Tianjin-Hebei (BTH) region, China over 2000, 2005 and 2010. We computed the freshwater ecosystem service provision with InVEST model. We calculated freshwater ecosystem service consumption with water quota method. We simulated the freshwater ecosystem service flow using our simplified flow model and assessed the regional water security with the improved freshwater security index. The freshwater provision service mainly depends on climatic factors that cannot be influenced by management, while the freshwater consumption service is constrained by human activities. Furthermore, the decrease of water quota for agricultural, domestic and industrial water counteracts the impact of increasing freshwater demand. The analysis of freshwater ecosystem service flow reveals that the majority area of the BTH (69.2%) is affected by upstream freshwater. If freshwater ecosystem service flow is considered, the water safety areas of the whole BTH account for 66.9%, 66.1%, 71.3%, which increase 6.4%, 6.8% and 5.7% in 2000, 2005 and 2010, respectively. These results highlight the need to understand the teleconnections between distant freshwater ecosystem service provision and local freshwater ecosystem service use. This approach therefore helps managers choose specific management and investment strategies for critical upstream freshwater provisions across different regions.

Simulation of ground-water flow and evaluation of water-management alternatives in the Assabet River Basin, Eastern Massachusetts

USGS Publications Warehouse

DeSimone, Leslie A.

2004-01-01

Water-supply withdrawals and wastewater disposal in the Assabet River Basin in eastern Massachusetts alter the flow and water quality in the basin. Wastewater discharges and stream-flow depletion from ground-water withdrawals adversely affect water quality in the Assabet River, especially during low-flow months (late summer) and in headwater areas. Streamflow depletion also contributes to loss of aquatic habitat in tributaries to the river. In 19972001, water-supply withdrawals averaged 9.9 million gallons per day (Mgal/d). Wastewater discharges to the Assabet River averaged 11 Mgal/d and included about 5.4 Mgal/d that originated from sources outside of the basin. The effects of current (2004) and future withdrawals and discharges on water resources in the basin were investigated in this study. Steady-state and transient ground-water-flow models were developed, by using MODFLOW-2000, to simulate flow in the surficial glacial deposits and underlying crystalline bedrock in the basin. The transient model simulated the average annual cycle at dynamic equilibrium in monthly intervals. The models were calibrated to 19972001 conditions of water withdrawals, wastewater discharges, water levels, and nonstorm streamflow (base flow plus wastewater discharges). Total flow through the simulated hydrologic system averaged 195 Mgal/d annually. Recharge from precipitation and ground-water discharge to streams were the dominant inflow and outflow, respectively. Evapotranspiration of ground water from wetlands and non-wetland areas also were important losses from the hydrologic system. Water-supply withdrawals and infiltration to sewers averaged 5 and 1.3 percent, respectively, of total annual out-flows and were larger components (12 percent in September) of the hydrologic system during low-flow months. Water budgets for individual tributary and main stem subbasins identified areas, such as the Fort Meadow Brook and the Assabet Main Stem Upper subbasins, where flows resulting from

Electrokinetic instability in microchannel ferrofluid/water co-flows

PubMed Central

Song, Le; Yu, Liandong; Zhou, Yilong; Antao, Asher Reginald; Prabhakaran, Rama Aravind; Xuan, Xiangchun

2017-01-01

Electrokinetic instability refers to unstable electric field-driven disturbance to fluid flows, which can be harnessed to promote mixing for various electrokinetic microfluidic applications. This work presents a combined numerical and experimental study of electrokinetic ferrofluid/water co-flows in microchannels of various depths. Instability waves are observed at the ferrofluid and water interface when the applied DC electric field is beyond a threshold value. They are generated by the electric body force that acts on the free charge induced by the mismatch of ferrofluid and water electric conductivities. A nonlinear depth-averaged numerical model is developed to understand and simulate the interfacial electrokinetic behaviors. It considers the top and bottom channel walls’ stabilizing effects on electrokinetic flow through the depth averaging of three-dimensional transport equations in a second-order asymptotic analysis. This model is found accurate to predict both the observed electrokinetic instability patterns and the measured threshold electric fields for ferrofluids of different concentrations in shallow microchannels. PMID:28406228

Ground-Water Flow Model of the Sierra Vista Subwatershed and Sonoran Portions of the Upper San Pedro Basin, Southeastern Arizona, United States, and Northern Sonora, Mexico

USGS Publications Warehouse

Pool, D.R.; Dickinson, Jesse

2007-01-01

A numerical ground-water model was developed to simulate seasonal and long-term variations in ground-water flow in the Sierra Vista subwatershed, Arizona, United States, and Sonora, Mexico, portions of the Upper San Pedro Basin. This model includes the simulation of details of the groundwater flow system that were not simulated by previous models, such as ground-water flow in the sedimentary rocks that surround and underlie the alluvial basin deposits, withdrawals for dewatering purposes at the Tombstone mine, discharge to springs in the Huachuca Mountains, thick low-permeability intervals of silt and clay that separate the ground-water flow system into deep-confined and shallow-unconfined systems, ephemeral-channel recharge, and seasonal variations in ground-water discharge by wells and evapotranspiration. Steady-state and transient conditions during 1902-2003 were simulated by using a five-layer numerical ground- water flow model representing multiple hydrogeologic units. Hydraulic properties of model layers, streamflow, and evapotranspiration rates were estimated as part of the calibration process by using observed water levels, vertical hydraulic gradients, streamflow, and estimated evapotranspiration rates as constraints. Simulations approximate observed water-level trends throughout most of the model area and streamflow trends at the Charleston streamflow-gaging station on the San Pedro River. Differences in observed and simulated water levels, streamflow, and evapotranspiration could be reduced through simulation of climate-related variations in recharge rates and recharge from flood-flow infiltration.

Conceptual model and numerical simulation of the ground-water-flow system in the unconsolidated deposits of the Colville River Watershed, Stevens County, Washington

USGS Publications Warehouse

Ely, D. Matthew; Kahle, Sue C.

2004-01-01

Increased use of ground- and surface-water supplies in watersheds of Washington State in recent years has created concern that insufficient instream flows remain for fish and other uses. Issuance of new ground-water rights in the Colville River Watershed was halted by the Washington Department of Ecology due to possible hydraulic continuity of the ground and surface waters. A ground-water-flow model was developed to aid in the understanding of the ground-water system and the regional effects of ground-water development alternatives on the water resources of the Colville River Watershed. The Colville River Watershed is underlain by unconsolidated deposits of glacial and non-glacial origin. The surficial geologic units and the deposits at depth were differentiated into aquifers and confining units on the basis of areal extent and general water-bearing characteristics. Five principal hydrogeologic units are recognized in the study area and form the basis of the ground-water-flow model. A steady-state ground-water-flow model of the Colville River Watershed was developed to simulate September 2001 conditions. The simulation period represented a period of below-average precipitation. The model was calibrated using nonlinear regression to minimize the weighted differences or residuals between simulated and measured hydraulic head and stream discharge. Simulated inflow to the model area was 53,000 acre-feet per year (acre-ft/yr) from precipitation and secondary recharge, and 36,000 acre-ft/yr from stream and lake leakage. Simulated outflow from the model was primarily through discharge to streams and lakes (71,000 acre-ft/yr), ground-water outflow (9,000 acre-ft/yr), and ground-water withdrawals (9,000 acre-ft/yr). Because the period of simulation, September 2001, was extremely dry, all components of the ground-water budget are presumably less than average flow conditions. The calibrated model was used to simulate the possible effects of increased ground-water pumping

Tafilalet OASIS System: Water Resources Management and Investigation by GIS and Groundwater Flow Model

NASA Astrophysics Data System (ADS)

Bouaamlat, I.; Larabi, A.; Faouzi, M.

2014-12-01

The geographical location of Tafilalet oasis system (TOS) in the south of the valley of Ziz (Morocco) offers him a particular advantage on the plane of water potential. The surface water which comes from humid regions of the High Atlas and intercepted by a dam then converged through the watercourse of Ziz towards the plain of the TOS, have created the conditions for the formation of a water table relatively rich with regard to the local climatic conditions (arid climate with recurrent drought). Given the role of the water table in the economic development of the region, a hydrogeological study was conducted to understand the impact of artificial recharge and recurrent droughts on the development of the groundwater reserves of TOS. In this study, a three-dimensional model of groundwater flow was developed for the TOS, to assist the decision makers as a "management tool" in order to assess alternative schemes for development and exploitation of groundwater resources based on the variation of artificial recharge and drought. The results from this numerical investigation of the TOS aquifer shows that the commissioning of the dam to control the flows of extreme flood and good management of water releases, has avoided the losses of irrigation water and consequently the non-overexploitation of the groundwater. So that with one or two water releases per year from the dam of flow rate more than 28 million m3/year it is possible to reconstruct the volume of water abstracted by wells. The idea of lowering water table by pumping wells is not exactly true, as well the development of groundwater abstraction has not prevented the wound of water table in these last years, the pumping wells accompanied more than it triggers the lowering of water table and it is mainly the succession of dry periods causing the decreases of the piezometric level. This situation confirms the important role that groundwater plays as a "buffer" during the drought periods.

South Asia river-flow projections and their implications for water resources

NASA Astrophysics Data System (ADS)

Mathison, C.; Wiltshire, A. J.; Falloon, P.; Challinor, A. J.

2015-12-01

South Asia is a region with a large and rising population, a high dependence on water intense industries, such as agriculture and a highly variable climate. In recent years, fears over the changing Asian summer monsoon (ASM) and rapidly retreating glaciers together with increasing demands for water resources have caused concern over the reliability of water resources and the potential impact on intensely irrigated crops in this region. Despite these concerns, there is a lack of climate simulations with a high enough resolution to capture the complex orography, and water resource analysis is limited by a lack of observations of the water cycle for the region. In this paper we present the first 25 km resolution regional climate projections of river flow for the South Asia region. Two global climate models (GCMs), which represent the ASM reasonably well are downscaled (1960-2100) using a regional climate model (RCM). In the absence of robust observations, ERA-Interim reanalysis is also downscaled providing a constrained estimate of the water balance for the region for comparison against the GCMs (1990-2006). The RCM river flow is routed using a river-routing model to allow analysis of present-day and future river flows through comparison with available river gauge observations. We examine how useful these simulations are for understanding potential changes in water resources for the South Asia region. In general the downscaled GCMs capture the seasonality of the river flows but overestimate the maximum river flows compared to the observations probably due to a positive rainfall bias and a lack of abstraction in the model. The simulations suggest an increasing trend in annual mean river flows for some of the river gauges in this analysis, in some cases almost doubling by the end of the century. The future maximum river-flow rates still occur during the ASM period, with a magnitude in some cases, greater than the present-day natural variability. Increases in river flow

Analytical solutions for benchmarking cold regions subsurface water flow and energy transport models: one-dimensional soil thaw with conduction and advection

USGS Publications Warehouse

Kurylyk, Barret L.; McKenzie, Jeffrey M; MacQuarrie, Kerry T. B.; Voss, Clifford I.

2014-01-01

Numerous cold regions water flow and energy transport models have emerged in recent years. Dissimilarities often exist in their mathematical formulations and/or numerical solution techniques, but few analytical solutions exist for benchmarking flow and energy transport models that include pore water phase change. This paper presents a detailed derivation of the Lunardini solution, an approximate analytical solution for predicting soil thawing subject to conduction, advection, and phase change. Fifteen thawing scenarios are examined by considering differences in porosity, surface temperature, Darcy velocity, and initial temperature. The accuracy of the Lunardini solution is shown to be proportional to the Stefan number. The analytical solution results obtained for soil thawing scenarios with water flow and advection are compared to those obtained from the finite element model SUTRA. Three problems, two involving the Lunardini solution and one involving the classic Neumann solution, are recommended as standard benchmarks for future model development and testing.

Cavity formation and surface modeling of laser milling process under a thin-flowing water layer

NASA Astrophysics Data System (ADS)

Tangwarodomnukun, Viboon

2016-11-01

Laser milling process normally involves a number of laser scans over a workpiece to selectively remove the material and then to form cavities with shape and dimensions required. However, this process adversely causes a heat accumulation in work material, which can in turn damage the laser-milled area and vicinity in terms of recast deposition and change of material properties. Laser milling process performing in a thin-flowing water layer is a promising method that can overcome such damage. With the use of this technique, water can flush away the cut debris and at the same time cool the workpiece during the ablation. To understand the potential of this technique for milling application, the effects of process parameters on cavity dimensions and surface roughness were experimentally examined in this study. Titanium sheet was used as a workpiece to be milled by a nanosecond pulse laser under different water flow velocities. A smooth and uniform cut feature can be obtained when the metal was ablated under the high laser pulse frequency and high water flow velocity. Furthermore, a surface model based on the energy balance was developed in this study to predict the cavity profile and surface roughness. By comparing to the experiments, the predicted profiles had a good agreement with the measured ones.

Geohydrology and conceptual model of a ground-water-flow system near a Superfund site in Cheshire, Connecticut

USGS Publications Warehouse

Stone, J.R.; Barlow, P.M.; Starn, J.J.

1996-01-01

Degradation of ground-water quality has been identified in an area of the north-central part of the town of Cheshire, Connecticut. An investigation by the U.S. Geological Survey, in cooperation with the U.S. Environmental Protection Agency, was done during 1994-95 to characterize the unconsolidated glacial deposits and the sedimentary bedrock, integrate the local geohydrologic conditions with the regional geohydrologic system, and develop a conceptual understanding of ground-water flow in the study area. A regional ground-water-flow model developed for the region near the study area indicates that perennial streams, including Judd Brook and the Tenmile River, form hydrologic divides that separate the larger region into hydraulically independent flow systems. In the local study area, synoptic water-level measurements made in June 1995 indicate that ground water near the water table flows west and southwestward from the low hill on the eastern side of the area toward the pond and wetlands along Judd Brook. Water-level data indicate that there is good hydraulic connection between the unconsolidated materials and underlying fractured bedrock. Unconsolidated materials in the study area consist principally of glacial stratified deposits that are fine sand, silt, and clay of glaci- olacustrine origin; locally these overlie thin glacial till. The glacial sediments range in thickness from a few feet to about 25 ft in the eastern part of the study area and are as much as 100 ft thick in the western and southern part of the study area beneath the Judd Brook and Tenmile River valleys. Fluvial redbeds of the New Haven Arkose underlie the glacial deposits in the region; in the study area, the redbeds consist of (1) channel sandstone units, which are coarse sandstone to fine conglomerate, generally in 6- to 15-ft- thick sequences; and (2) overbank mudstone units, which are siltstone and silty sandstone with some fine sandstone, generally in 6- to 50-ft-thick sequences. Thin

Water flow through temperate glaciers

USGS Publications Warehouse

Fountain, A.G.; Walder, J.S.

1998-01-01

Understanding water movement through a glacier is fundamental to several critical issues in glaciology, including glacier dynamics, glacier-induced floods, and the prediction of runoff from glacierized drainage basins. to this end we have synthesized a conceptual model os water movement through a temperate glacier from the surface to the outlet stream. Processes that regulate the rate and distribution of water input at the glacier surface and that regulate water movement from the surface to the bed play important but commonly neglected roles in glacier hydrology. Where a glacier is covered by a layer of porous, permeable firn (the accumulation zone), the flux of water to the glacier interior varies slowly because the firn temporarily stores water and thereby smooths out variations in the supply rate. In the firn-free ablation zone, in contrast, the flux of water into the glacier depends directly on the rate of surface melt or rainfall and therefore varies greatly in time. Water moves from the surface to the bed through an upward branching arborescent network consisting of both steeply inclined conduits, formed by the enlargement of intergranular veins, and gently inclined conduits, sprqwned by water flow along the bottoms of near-surface fractures (crevasses). Englacial drainage conduits deliver water to the glacier bed at a linited number of points, probably a long distance downglacier of where water enters the glacier. Englacial conduits supplied from the accumulation zone are quasi steady state features that convey the slowly varying water flux delivered via the firn. their size adjusts so that they are usually full of water and flow is pressurized. In contrast, water flow in englacial conduits supplied from the ablation area is pressurized only near times of peak daily flow or during rainstorms; flow is otherwise in an open-channel configuration. The subglacial drainage system typically consists of several elements that are distinct both morpphologically and

Simulation of an urban ground-water-flow system in the Menomonee Valley, Milwaukee, Wisconsin using analytic element modeling

USGS Publications Warehouse

Dunning, C.P.; Feinstein, D.T.

2004-01-01

A single-layer, steady-state analytic element model was constructed to simulate shallow ground-water flow in the Menomonee Valley, an old industrial center southwest of downtown Milwaukee, Wisconsin. Project objectives were to develop an understanding of the shallow ground-water flow system and identify primary receptors of recharge to the valley. The analytic element model simulates flow in a 18.3 m (60 ft) thick layer of estuarine and alluvial sediments and man-made fill that comprises the shallow aquifer across the valley. The thin, laterally extensive nature of the shallow aquifer suggests horizontal-flow predominates, thus the system can appropriately be modeled with the Dupuit-Forchheimer approximation in an analytic element model. The model was calibrated to the measured baseflow increase between two USGS gages on the Menomonee River, 90 head measurements taken in and around the valley during December 1999, and vertical gradients measured at five locations under the river and estuary in the valley. Recent construction of the Milwaukee Metropolitan Sewer District Inline Storage System (ISS) in the Silurian dolomite under the Menomonee Valley has locally lowered heads in the dolomite appreciably, below levels caused by historic pumping. The ISS is a regional hydraulic sink which removes water from the bedrock even during dry weather. The potential effect on flow directions in the shallow aquifer of dry-weather infiltration to the ISS was evaluated by adjusting the resistance of the line-sink strings representing the ISS in the model to allow infiltration from 0 to 100% of the reported 9,500 m3/d. The best fit to calibration targets was found between 60% (5,700 m3/d) and 80% (7,600 m3/d) of the reported dry-weather infiltration. At 60% infiltration, 65% of the recharge falling on the valley terminates at the ISS and 35% at the Menomonee River and estuary. At 80% infiltration, 73% of the recharge terminates at the ISS, and 27% at the river and estuary. Model

Modelling Water Flow through Paddy Soils under Alternate Wetting and Drying Irrigation Practice

NASA Astrophysics Data System (ADS)

Shekhar, S.; Mailapalli, D. R.; Das, B. S.; Raghuwanshi, N. S.

2017-12-01

Alternate wetting and drying (AWD) irrigation practice in paddy cultivation requires an optimum soil moisture stress (OSMS) level at which irrigation water savings can be maximized without compromising the yield reduction. Determining OSMS experimentally is challenging and only possible with appropriate modeling tools. In this study, field experiments on paddy were conducted in thirty non-weighing type lysimeters during dry seasons of 2016 and 2017. Ten plots were irrigated using continuous flooding (CF) and the rest were irrigated with AWD practice at 40mb and 75mb soil moisture stress levels. Depth of ponding and soil suction at 10, 40 and 70 cm from the soil surface were measured daily from all lysimeter plots. The measured field data were used in calibration and validation of Hydrus-1D model and simulated the water flow for both AWD and CF plots. The Hydrus-1D is being used to estimate OSMS for AWD practice and compared the seasonal irrigation water input and deep percolation losses with CF practice.

Integrated Modelling on Flow and Water Quality Under the Impacts of Climate Change and Agricultural Activities

NASA Astrophysics Data System (ADS)

SHI, J.

2014-12-01

Climate change is expected to have a significant impact on flooding in the UK, inducing more intense and prolonged storms. Frequent flooding due to climate change already exacerbates catchment water quality. Land use is another contributing factor to poor water quality. For example, the move to intensive farming could cause an increase in faecal coliforms entering the water courses. In an effort to understand better the effects on water quality from land use and climate change, the hydrological and estuarine processes are being modelled using SWAT (Soil and Water Assessment Tool), linked to a 2-D hydrodynamic model DIVAST(Depth Integrated Velocity and Solute Transport). The coupled model is able to quantify how much of each pollutant from the catchment reaches the harbour and the impact on water quality within the harbour. The work is focused on the transportation and decay of faecal coliforms from agricultural runoff into the rivers Frome and Piddle in the UK. The impact from the agricultural land use and activities on the catchment river hydrology and water quality are evaluated. The coupled model calibration and validation showed the good model performance on flow and faecal coliform in the watershed and estuary.

Comparison of vertical discretization techniques in finite-difference models of ground-water flow; example from a hypothetical New England setting

USGS Publications Warehouse

Harte, Philip T.

1994-01-01

Proper discretization of a ground-water-flow field is necessary for the accurate simulation of ground-water flow by models. Although discretiza- tion guidelines are available to ensure numerical stability, current guidelines arc flexible enough (particularly in vertical discretization) to allow for some ambiguity of model results. Testing of two common types of vertical-discretization schemes (horizontal and nonhorizontal-model-layer approach) were done to simulate sloping hydrogeologic units characteristic of New England. Differences of results of model simulations using these two approaches are small. Numerical errors associated with use of nonhorizontal model layers are small (4 percent). even though this discretization technique does not adhere to the strict formulation of the finite-difference method. It was concluded that vertical discretization by means of the nonhorizontal layer approach has advantages in representing the hydrogeologic units tested and in simplicity of model-data input. In addition, vertical distortion of model cells by this approach may improve the representation of shallow flow processes.

Surface velocity divergence model of air/water interfacial gas transfer in open-channel flows

NASA Astrophysics Data System (ADS)

Sanjou, M.; Nezu, I.; Okamoto, T.

2017-04-01

Air/water interfacial gas transfer through a free surface plays a significant role in preserving and restoring water quality in creeks and rivers. However, direct measurements of the gas transfer velocity and reaeration coefficient are still difficult, and therefore a reliable prediction model needs to be developed. Varying systematically the bulk-mean velocity and water depth, laboratory flume experiments were conducted and we measured surface velocities and dissolved oxygen (DO) concentrations in open-channel flows to reveal the relationship between DO transfer velocity and surface divergence (SD). Horizontal particle image velocimetry measurements provide the time-variations of surface velocity divergence. Positive and negative regions of surface velocity divergence are transferred downstream in time, as occurs in boil phenomenon on natural river free-surfaces. The result implies that interfacial gas transfer is related to bottom-situated turbulence motion and vertical mass transfer. The original SD model focuses mainly on small-scale viscous motion, and this model strongly depends on the water depth. Therefore, we modify the SD model theoretically to accommodate the effects of the water depth on gas transfer, introducing a non-dimensional parameter that includes contributions of depth-scale large-vortex motion, such as secondary currents, to surface renewal events related to DO transport. The modified SD model proved effective and reasonable without any dependence on the bulk mean velocity and water depth, and has a larger coefficient of determination than the original SD model. Furthermore, modeling of friction velocity with the Reynolds number improves the practicality of a new formula that is expected to be used in studies of natural rivers.

Well balancing of the SWE schemes for moving-water steady flows

NASA Astrophysics Data System (ADS)

Caleffi, Valerio; Valiani, Alessandro

2017-08-01

In this work, the exact reproduction of a moving-water steady flow via the numerical solution of the one-dimensional shallow water equations is studied. A new scheme based on a modified version of the HLLEM approximate Riemann solver (Dumbser and Balsara (2016) [18]) that exactly preserves the total head and the discharge in the simulation of smooth steady flows and that correctly dissipates mechanical energy in the presence of hydraulic jumps is presented. This model is compared with a selected set of schemes from the literature, including models that exactly preserve quiescent flows and models that exactly preserve moving-water steady flows. The comparison highlights the strengths and weaknesses of the different approaches. In particular, the results show that the increase in accuracy in the steady state reproduction is counterbalanced by a reduced robustness and numerical efficiency of the models. Some solutions to reduce these drawbacks, at the cost of increased algorithm complexity, are presented.

Dynamics of drinking water biofilm in flow/non-flow conditions.

PubMed

Manuel, C M; Nunes, O C; Melo, L F

2007-02-01

Drinking water biofilm formation on polyvinyl chloride (PVC), cross-linked polyethylene (PEX), high density polyethylene (HDPE) and polypropylene (PP) was followed in three different reactors operating under stagnant or continuous flow regimes. After one week, a quasi-steady state was achieved where biofilm total cell numbers per unit surface area were not affected by fluctuations in the concentration of suspended cells. Metabolically active cells in biofilms were around 17-35% of the total cells and 6-18% were able to form colony units in R(2)A medium. Microbiological analysis showed that the adhesion material and reactor design did not affect significantly the biofilm growth. However, operating under continuous flow (0.8-1.9 Pa) or stagnant water had a significant effect on biofilm formation: in stagnant waters, biofilm grew to a less extent. By applying mass balances and an asymptotic biofilm formation model to data from biofilms grown on PVC and HDPE surfaces under turbulent flow, specific growth rates of bacteria in the biofilm were found to be similar for both materials (around 0.15 day(-1)) and much lower than the specific growth rates of suspended bacteria (around 1.8 day(-1)).

Global scale groundwater flow model

NASA Astrophysics Data System (ADS)

Sutanudjaja, Edwin; de Graaf, Inge; van Beek, Ludovicus; Bierkens, Marc

2013-04-01

As the world's largest accessible source of freshwater, groundwater plays vital role in satisfying the basic needs of human society. It serves as a primary source of drinking water and supplies water for agricultural and industrial activities. During times of drought, groundwater sustains water flows in streams, rivers, lakes and wetlands, and thus supports ecosystem habitat and biodiversity, while its large natural storage provides a buffer against water shortages. Yet, the current generation of global scale hydrological models does not include a groundwater flow component that is a crucial part of the hydrological cycle and allows the simulation of groundwater head dynamics. In this study we present a steady-state MODFLOW (McDonald and Harbaugh, 1988) groundwater model on the global scale at 5 arc-minutes resolution. Aquifer schematization and properties of this groundwater model were developed from available global lithological model (e.g. Dürr et al., 2005; Gleeson et al., 2010; Hartmann and Moorsdorff, in press). We force the groundwtaer model with the output from the large-scale hydrological model PCR-GLOBWB (van Beek et al., 2011), specifically the long term net groundwater recharge and average surface water levels derived from routed channel discharge. We validated calculated groundwater heads and depths with available head observations, from different regions, including the North and South America and Western Europe. Our results show that it is feasible to build a relatively simple global scale groundwater model using existing information, and estimate water table depths within acceptable accuracy in many parts of the world.

Estimating Monthly Water Withdrawals, Return Flow, and Consumptive Use in the Great Lakes Basin

USGS Publications Warehouse

Shaffer, Kimberly H.; Stenback, Rosemary S.

2010-01-01

Water-resource managers and planners require water-withdrawal, return-flow, and consumptive-use data to understand how anthropogenic (human) water use affects the hydrologic system. Water models like MODFLOW and GSFLOW use calculations and input values (including water-withdrawal and return flow data) to simulate and predict the effects of water use on aquifer and stream conditions. Accurate assessments of consumptive use, interbasin transfer, and areas that are on public supply or sewer are essential in estimating the withdrawal and return-flow data needed for the models. As the applicability of a model to real situations depends on accurate input data, limited or poor water-use data hampers the ability of modelers to simulate and predict hydrologic conditions. Substantial differences exist among the many agencies nationwide that are responsible for compiling water-use data including what data are collected, how the data are organized, how often the data are collected, quality assurance, required level of accuracy, and when data are released to the public. This poster presents water-use information and estimation methods summarized from recent U.S. Geological Survey (USGS) reports with the intent to assist water-resource managers and planners who need estimates of monthly water withdrawals, return flows, and consumptive use. This poster lists references used in Shaffer (2009) for water withdrawals, consumptive use, and return flows. Monthly percent of annual withdrawals and monthly consumptive-use coefficients are used to compute monthly water withdrawals, consumptive use, and return flow for the Great Lakes Basin.

Simulation of ground-water flow in the Cedar River alluvial aquifer flow system, Cedar Rapids, Iowa

USGS Publications Warehouse

Turco, Michael J.; Buchmiller, Robert C.

2004-01-01

Model results indicate that the primary sources of inflow to the modeled area are infiltration from the Cedar River (53.0 percent) and regional flow in the glacial and bedrock materials (34.1 percent). The primary sources of outflow from the modeled area are discharge to the Cedar River (45.4 percent) and pumpage (44.8 percent). Current steady-state pumping rates have increased the flow of water from the Cedar River to the alluvial aquifer by 43.8 cubic feet per second. Steady-state and transient hypothetical pumpage scenarios were used to show the relation between changes in pumpage and changes in infiltration of water from the Cedar River. Results indicate that more than 99 percent of the water discharging from municipal wells infiltrates from the Cedar River, that the time required for induced river recharge to equilibrate with municipal pumpage may be 150 days or more, and that ground-water availability in the Cedar Rapids area will not be significantly affected by doubling current pumpage as long as there is sufficient flow in the Cedar River to provide recharge.

Modeling two-phase flow in three-dimensional complex flow-fields of proton exchange membrane fuel cells

NASA Astrophysics Data System (ADS)

Kim, Jinyong; Luo, Gang; Wang, Chao-Yang

2017-10-01

3D fine-mesh flow-fields recently developed by Toyota Mirai improved water management and mass transport in proton exchange membrane (PEM) fuel cell stacks, suggesting their potential value for robust and high-power PEM fuel cell stack performance. In such complex flow-fields, Forchheimer's inertial effect is dominant at high current density. In this work, a two-phase flow model of 3D complex flow-fields of PEMFCs is developed by accounting for Forchheimer's inertial effect, for the first time, to elucidate the underlying mechanism of liquid water behavior and mass transport inside 3D complex flow-fields and their adjacent gas diffusion layers (GDL). It is found that Forchheimer's inertial effect enhances liquid water removal from flow-fields and adds additional flow resistance around baffles, which improves interfacial liquid water and mass transport. As a result, substantial improvements in high current density cell performance and operational stability are expected in PEMFCs with 3D complex flow-fields, compared to PEMFCs with conventional flow-fields. Higher current density operation required to further reduce PEMFC stack cost per kW in the future will necessitate optimizing complex flow-field designs using the present model, in order to efficiently remove a large amount of product water and hence minimize the mass transport voltage loss.

MODFLOW-2000, the U.S. Geological Survey Modular Ground-Water Model -Documentation of the Hydrogeologic-Unit Flow (HUF) Package

USGS Publications Warehouse

Anderman, E.R.; Hill, M.C.

2000-01-01

This report documents the Hydrogeologic-Unit Flow (HUF) Package for the groundwater modeling computer program MODFLOW-2000. The HUF Package is an alternative internal flow package that allows the vertical geometry of the system hydrogeology to be defined explicitly within the model using hydrogeologic units that can be different than the definition of the model layers. The HUF Package works with all the processes of MODFLOW-2000. For the Ground-Water Flow Process, the HUF Package calculates effective hydraulic properties for the model layers based on the hydraulic properties of the hydrogeologic units, which are defined by the user using parameters. The hydraulic properties are used to calculate the conductance coefficients and other terms needed to solve the ground-water flow equation. The sensitivity of the model to the parameters defined within the HUF Package input file can be calculated using the Sensitivity Process, using observations defined with the Observation Process. Optimal values of the parameters can be estimated by using the Parameter-Estimation Process. The HUF Package is nearly identical to the Layer-Property Flow (LPF) Package, the major difference being the definition of the vertical geometry of the system hydrogeology. Use of the HUF Package is illustrated in two test cases, which also serve to verify the performance of the package by showing that the Parameter-Estimation Process produces the true parameter values when exact observations are used.

Gas-liquid mass transfer and flow phenomena in the Peirce-Smith converter: a water model study

NASA Astrophysics Data System (ADS)

Zhao, Xing; Zhao, Hong-liang; Zhang, Li-feng; Yang, Li-qiang

2018-01-01

A water model with a geometric similarity ratio of 1:5 was developed to investigate the gas-liquid mass transfer and flow characteristics in a Peirce-Smith converter. A gas mixture of CO2 and Ar was injected into a NaOH solution bath. The flow field, volumetric mass transfer coefficient per unit volume ( Ak/V; where A is the contact area between phases, V is the volume, and k is the mass transfer coefficient), and gas utilization ratio ( η) were then measured at different gas flow rates and blow angles. The results showed that the flow field could be divided into five regions, i.e., injection, strong loop, weak loop, splashing, and dead zone. Whereas the Ak/V of the bath increased and then decreased with increasing gas flow rate, and η steadily increased. When the converter was rotated clockwise, both Ak/V and η increased. However, the flow condition deteriorated when the gas flow rate and blow angle were drastically increased. Therefore, these parameters must be controlled to optimal conditions. In the proposed model, the optimal gas flow rate and blow angle were 7.5 m3·h-1 and 10°, respectively.

Hydrogeologic Setting and Ground-Water Flow in the Leetown Area, West Virginia

USGS Publications Warehouse

Kozar, Mark D.; Weary, David J.; Paybins, Katherine S.; Pierce, Herbert A.

2007-01-01

The Leetown Science Center is a research facility operated by the U.S. Geological Survey that occupies approximately 455-acres near Kearneysville, Jefferson County, West Virginia. Aquatic and fish research conducted at the Center requires adequate supplies of high-quality, cold ground water. Three large springs and three production wells currently (in 2006) supply water to the Center. The recent construction of a second research facility (National Center for Cool and Cold Water Aquaculture) operated by the U.S. Department of Agriculture and co-located on Center property has placed additional demands on available water resources in the area. A three-dimensional steady-state finite-difference ground-water flow model was developed to simulate ground-water flow in the Leetown area and was used to assess the availability of ground water to sustain current and anticipated future demands. The model also was developed to test a conceptual model of ground-water flow in the complex karst aquifer system in the Leetown area. Due to the complexity of the karst aquifer system, a multidisciplinary research study was required to define the hydrogeologic setting. Geologic mapping, surface- and borehole-geophysical surveys, stream base-flow surveys, and aquifer tests were conducted to provide the hydrogeologic data necessary to develop and calibrate the model. It would not have been possible to develop a numerical model of the study area without the intensive data collection and methods developments components of the larger, more comprehensive hydrogeologic investigation. Results of geologic mapping and surface-geophysical surveys verified the presence of several prominent thrust faults and identified additional faults and other complex geologic structures (including overturned anticlines and synclines) in the area. These geologic structures are known to control ground-water flow in the region. Results of this study indicate that cross-strike faults and fracture zones are major

Scarce water resources and scarce data: Estimating recharge for a complex 3D groundwater flow model in arid regions

NASA Astrophysics Data System (ADS)

Gräbe, A. C.; Guttman, J.; Rödiger, T.; Siebert, C.; Merz, R.; Kolditz, O.

2012-12-01

Semi-arid to arid regions are usually characterized by a scarcity of precipitation and a lack of stream flow. Especially in desert environments, groundwater is one of the most important fresh water sources and its recharge is basically controlled by two main mechanisms: the direct regional infiltration of precipitation in the mountains and interdrainage areas in the first place and secondly the flood water infiltration through ephemeral channel beds (transmission loss). Due to extensive spatio-temporal data scarcity, direct quantitative estimations of groundwater recharge are often difficult to perform, and numerical models simulating the water fluxes, have to be applied to enable a quantitative approximation of the groundwater recharge. We made an assumption about the quantity of recharge for the subsurface catchment of the western Dead Sea escarpment, which is at the same time the input for the complex groundwater flow model of the Judea Group Aquifer. This can only be suggested if the hydrogeological situation in the tectonically complex region is fully understood. A number of simplified models of the Judea Group aquifer have been formulated and employed using a two-dimensional (one horizontal layered) numerical simulation of groundwater flow (Baida et al. 1978; Goldschtoff & Shachnai, 1980; Guttman, 2000; Laronne Ben-Itzhak & Gvirtzmann, 2005). However, all previous approaches focused only on a limited area of the Judea Group aquifer. We developed a high resolution regional groundwater flow model for the entire western basin of the Dead Sea. Whereas the structural model could be defined using a large geological dataset, the challenge was to generate the groundwater flow model with only limited well data. With the help of the scientific software OpenGeoSys (OGS) the challenge was reliably solved resulting in a simulation of the hydraulic characteristics (hydraulic conductivity and hydraulic head) of the cretaceous aquifer system, which was calibrated using PEST.

Moderated, Water-Based, Condensational Particle Growth in a Laminar Flow

PubMed Central

Hering, Susanne V.; Spielman, Steven R.; Lewis, Gregory S.

2014-01-01

Presented is a new approach for laminar-flow water condensation that produces saturations above 1.5 while maintaining temperatures of less than 30°C in the majority of the flow and providing an exiting dew point below 15°C. With the original laminar flow water condensation method, the particle activation and growth occurs in a region with warm, wetted walls throughout, which has the side-effect of heating the flow. The “moderated” approach presented here replaces this warm region with a two sections – a short, warm, wet-walled “initiator”, followed by a cool-walled “moderator”. The initiator provides the water vapor that creates the supersaturation, while the moderator provides the time for particle growth. The combined length of the initiator and moderator sections is the same as that of the original, warm-walled growth section. Model results show that this new approach reduces the added heat and water vapor while achieving the same peak supersaturation and similar droplet growth. Experimental measurements confirm the trends predicted by the modeling. PMID:24839342

Simulation of ground-water flow in the Vevay Township area, Ingham County, Michigan

USGS Publications Warehouse

Luukkonen, Carol L.; Simard, Andreanne

2004-01-01

Ground water is the primary source of water for domestic, public-supply, and industrial use within the Tri-County region that includes Clinton, Eaton, and Ingham Counties in Michigan. Because of the importance of this ground-water resource, numerous communities, including the city of Mason in Ingham County, have begun local Wellhead Protection Programs. In these programs, communities protect their groundwater resource by identifying the areas that contribute water to production wells and potential sources of contamination, and by developing methods to manage and minimize threats to the water supply. In addition, some communities in Michigan are concerned about water availability, particularly in areas experiencing water-level declines in the vicinity of quarry dewatering operations. In areas where Wellhead Protection Programs are implemented and there are potential threats to the water supply, residents and communities need adequate information to protect the water supply.In 1996, a regional ground-water-flow model was developed by the U.S. Geological Survey to simulate ground-water flow in Clinton, Eaton, and Ingham Counties. This model was developed primarily to simulate the bedrock ground-waterflow system; ground-water flow in the unconsolidated glacial sediments was simulated to support analysis of flow in the underlying bedrock Saginaw aquifer. Since its development in 1996, regional model simulations have been conducted to address protection concerns and water availability questions of local water-resources managers. As a result of these continuing model simulations, additional hydrogeologic data have been acquired in the Tri-County region that has improved the characterization of the simulated ground-water-flow system and improved the model calibration. A major benefit of these updates and refinements is that the regional Tri-County model continues to be a useful tool that improves the understanding of the ground-water-flow system in the Tri-County region

Two-dimensional advective transport in ground-water flow parameter estimation

USGS Publications Warehouse

Anderman, E.R.; Hill, M.C.; Poeter, E.P.

1996-01-01

Nonlinear regression is useful in ground-water flow parameter estimation, but problems of parameter insensitivity and correlation often exist given commonly available hydraulic-head and head-dependent flow (for example, stream and lake gain or loss) observations. To address this problem, advective-transport observations are added to the ground-water flow, parameter-estimation model MODFLOWP using particle-tracking methods. The resulting model is used to investigate the importance of advective-transport observations relative to head-dependent flow observations when either or both are used in conjunction with hydraulic-head observations in a simulation of the sewage-discharge plume at Otis Air Force Base, Cape Cod, Massachusetts, USA. The analysis procedure for evaluating the probable effect of new observations on the regression results consists of two steps: (1) parameter sensitivities and correlations calculated at initial parameter values are used to assess the model parameterization and expected relative contributions of different types of observations to the regression; and (2) optimal parameter values are estimated by nonlinear regression and evaluated. In the Cape Cod parameter-estimation model, advective-transport observations did not significantly increase the overall parameter sensitivity; however: (1) inclusion of advective-transport observations decreased parameter correlation enough for more unique parameter values to be estimated by the regression; (2) realistic uncertainties in advective-transport observations had a small effect on parameter estimates relative to the precision with which the parameters were estimated; and (3) the regression results and sensitivity analysis provided insight into the dynamics of the ground-water flow system, especially the importance of accurate boundary conditions. In this work, advective-transport observations improved the calibration of the model and the estimation of ground-water flow parameters, and use of

Modelling Soil Heat and Water Flow as a Coupled Process in Land Surface Models

NASA Astrophysics Data System (ADS)

García González, Raquel; Verhoef, Anne; Vidale, Pier Luigi; Braud, Isabelle

2010-05-01

To improve model estimates of soil water and heat flow by land surface models (LSMs), in particular in the first few centimetres of the near-surface soil profile, we have to consider in detail all the relevant physical processes involved (see e.g. Milly, 1982). Often, thermal and iso-thermal vapour fluxes in LSMs are neglected and the simplified Richard's equation is used as a result. Vapour transfer may affect the water fluxes and heat transfer in LSMs used for hydrometeorological and climate simulations. Processes occurring in the top 50 cm soil may be relevant for water and heat flux dynamics in the deeper layers, as well as for estimates of evapotranspiration and heterotrophic respiration, or even for climate and weather predictions. Water vapour transfer, which was not incorporated in previous versions of the MOSES/JULES model (Joint UK Land Environment Simulator; Cox et al., 1999), has now been implemented. Furthermore, we also assessed the effect of the soil vertical resolution on the simulated soil moisture and temperature profiles and the effect of the processes occurring at the upper boundary, mainly in terms of infiltration rates and evapotranspiration. SiSPAT (Simple Soil Plant Atmosphere Transfer Model; Braud et al., 1995) was initially used to quantify the changes that we expect to find when we introduce vapour transfer in JULES, involving parameters such as thermal vapour conductivity and diffusivity. Also, this approach allows us to compare JULES to a more complete and complex numerical model. Water vapour flux varied with soil texture, depth and soil moisture content, but overall our results suggested that water vapour fluxes change temperature gradients in the entire soil profile and introduce an overall surface cooling effect. Increasing the resolution smoothed and reduced temperature differences between liquid (L) and liquid/vapour (LV) simulations at all depths, and introduced a temperature increase over the entire soil profile. Thermal

Non-Darcy flow of water through woodchip media

USDA-ARS?s Scientific Manuscript database

A denitrifying bioreactor is a system where a carbon substrate (commonly woodchips) is used to reduce nitrate concentration in water flow. Knowledge of intrinsic permeability of woodchip media in different types of this system is of great importance for design and modeling. For many years, water flo...

Flow path oscillations in transient ground-water simulations of large peatland systems

USGS Publications Warehouse

Reeve, A.S.; Evensen, R.; Glaser, P.H.; Siegel, D.I.; Rosenberry, D.

2006-01-01

Transient numerical simulations of the Glacial Lake Agassiz Peatland near the Red Lakes in Northern Minnesota were constructed to evaluate observed reversals in vertical ground-water flow. Seasonal weather changes were introduced to a ground-water flow model by varying evapotranspiration and recharge over time. Vertical hydraulic reversals, driven by changes in recharge and evapotranspiration were produced in the simulated peat layer. These simulations indicate that the high specific storage associated with the peat is an important control on hydraulic reversals. Seasonally driven vertical flow is on the order of centimeters in the deep peat, suggesting that seasonal vertical advective fluxes are not significant and that ground-water flow into the deep peat likely occurs on decadal or longer time scales. Particles tracked within the ground-water flow model oscillate over time, suggesting that seasonal flow reversals will enhance vertical mixing in the peat column. The amplitude of flow path oscillations increased with increasing peat storativity, with amplitudes of about 5 cm occurring when peat specific storativity was set to about 0.05 m-1. ?? 2005 Elsevier B.V. All rights reserved.

Squirt flow due to interfacial water films in hydrate bearing sediments

NASA Astrophysics Data System (ADS)

Sell, Kathleen; Quintal, Beatriz; Kersten, Michael; Saenger, Erik H.

2018-05-01

Sediments containing gas hydrate dispersed in the pore space are known to show a characteristic seismic anomaly which is a high attenuation along with increasing seismic velocities. Currently, this observation cannot be fully explained albeit squirt-flow type mechanisms on the microscale have been speculated to be the cause. Recent major findings from in situ experiments, using the gas in excess and water in excess formation method, and coupled with high-resolution synchrotron-based X-ray micro-tomography, have revealed the systematic presence of thin water films between the quartz grains and the encrusting hydrate. The data obtained from these experiments underwent an image processing procedure to quantify the thicknesses and geometries of the aforementioned interfacial water films. Overall, the water films vary from sub-micrometer to a few micrometers in thickness. In addition, some of the water films interconnect through water bridges. This geometrical analysis is used to propose a new conceptual squirt flow model for hydrate bearing sediments. A series of numerical simulations is performed considering variations of the proposed model to study seismic attenuation caused by such thin water films. Our results support previous speculation that squirt flow can explain high attenuation at seismic frequencies in hydrate bearing sediments, but based on a conceptual squirt flow model which is geometrically different than those previously considered.

Effects of spatially distributed sectoral water management on the redistribution of water resources in an integrated water model

NASA Astrophysics Data System (ADS)

Voisin, Nathalie; Hejazi, Mohamad I.; Leung, L. Ruby; Liu, Lu; Huang, Maoyi; Li, Hong-Yi; Tesfa, Teklu

2017-05-01

Realistic representations of sectoral water withdrawals and consumptive demands and their allocation to surface and groundwater sources are important for improving modeling of the integrated water cycle. To inform future model development, we enhance the representation of water management in a regional Earth system (ES) model with a spatially distributed allocation of sectoral water demands simulated by a regional integrated assessment (IA) model to surface and groundwater systems. The integrated modeling framework (IA-ES) is evaluated by analyzing the simulated regulated flow and sectoral supply deficit in major hydrologic regions of the conterminous U.S, which differ from ES studies looking at water storage variations. Decreases in historical supply deficit are used as metrics to evaluate IA-ES model improvement in representating the complex sectoral human activities for assessing future adaptation and mitigation strategies. We also assess the spatial changes in both regulated flow and unmet demands, for irrigation and nonirrigation sectors, resulting from the individual and combined additions of groundwater and return flow modules. Results show that groundwater use has a pronounced regional and sectoral effect by reducing water supply deficit. The effects of sectoral return flow exhibit a clear east-west contrast in the hydrologic patterns, so the return flow component combined with the IA sectoral demands is a major driver for spatial redistribution of water resources and water deficits in the US. Our analysis highlights the need for spatially distributed sectoral representation of water management to capture the regional differences in interbasin redistribution of water resources and deficits.

Virtual water flows and trade liberalization.

PubMed

Ramirez-Vallejo, J; Rogers, P

2004-01-01

The linkages between agricultural trade and water resources need to be identified and analyzed to better understand the potential impacts that a full liberalization, or lack thereof, will have on water resources. This paper examines trade of virtual water embodied in agricultural products for most countries of the world. The main purpose of the paper, however, is to examine the impact of trade liberalization on virtual-water trade in the future. Based on a simulation of global agricultural trade, a scenario of full liberalization of agriculture was used to assess the net effect of virtual water flows from the relocation of meat and cereals' trade. The paper also identifies the main reasons behind the changes in the magnitude and direction of the net virtual water trade over time, and shows that virtual water trade flows are independent of water resource endowments, contrary to what the Heckscher-Ohlin Theorem states. Finally, based on a formal model, some input demand functions at the country level are estimated. The estimates of the income and agricultural support elasticities of demand for import of virtual water have the expected sign, and are statistically significant. Variables found to have some explanatory power of the variance of virtual water imports are average income; population; agriculture as value added; irrigated area, and exports of goods and services.

Enhanced water vapour flow in silica microchannels and interdiffusive water vapour flow through anodic aluminium oxide (AAO) membranes

NASA Astrophysics Data System (ADS)

Lei, Wenwen; McKenzie, David R.

2015-12-01

Enhanced liquid water flows through carbon nanotubes reinvigorated the study of moisture permeation through membranes and micro- and nano-channels. The study of water vapour through micro-and nano-channels has been neglected even though water vapour is as important as liquid water for industry, especially for encapsulation of electronic devices. Here we measure moisture flow rates in silica microchannels and interdiffusive water vapour flows in anodic aluminium oxide (AAO) membrane channels for the first time. We construct theory for the flow rates of the dominant modes of water transport through four previously defined standard configurations and benchmark it against our new measurements. The findings show that measurements of leak behaviour made using other molecules, such as helium, are not reliable. Single phase water vapour flow is overestimated by a helium measurement, while Washburn or capillary flow is underestimated or for all channels when boundary slip applies, to an extent that depends on the slip length for the liquid phase flows.

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

VISCOPLASTIC FLUID MODEL FOR DEBRIS FLOW ROUTING.

USGS Publications Warehouse

Chen, Cheng-lung

1986-01-01

This paper describes how a generalized viscoplastic fluid model, which was developed based on non-Newtonian fluid mechanics, can be successfully applied to routing a debris flow down a channel. The one-dimensional dynamic equations developed for unsteady clear-water flow can be used for debris flow routing if the flow parameters, such as the momentum (or energy) correction factor and the resistance coefficient, can be accurately evaluated. The writer's generalized viscoplastic fluid model can be used to express such flow parameters in terms of the rheological parameters for debris flow in wide channels. A preliminary analysis of the theoretical solutions reveals the importance of the flow behavior index and the so-called modified Froude number for uniformly progressive flow in snout profile modeling.

Physiological water model development

NASA Technical Reports Server (NTRS)

Doty, Susan

1993-01-01

The water of the human body can be categorized as existing in two main compartments: intracellular water and extracellular water. The intracellular water consists of all the water within the cells and constitutes over half of the total body water. Since red blood cells are surrounded by plasma, and all other cells are surrounded by interstitial fluid, the intracellular compartment has been subdivided to represent these two cell types. The extracellular water, which includes all of the fluid outside of the cells, can be further subdivided into compartments which represent the interstitial fluid, circulating blood plasma, lymph, and transcellular water. The interstitial fluid surrounds cells outside of the vascular system whereas plasma is contained within the blood vessels. Avascular tissues such as dense connective tissue and cartilage contain interstitial water which slowly equilibrates with tracers used to determine extracellular fluid volume. For this reason, additional compartments are sometimes used to represent these avascular tissues. The average size of each compartment, in terms of percent body weight, has been determined for adult males and females. These compartments and the forces which cause flow between them are presented. The kidneys, a main compartment, receive about 25 percent of the cardiac output and filters out a fluid similar to plasma. The composition of this filtered fluid changes as it flows through the kidney tubules since compounds are continually being secreted and reabsorbed. Through this mechanism, the kidneys eliminate wastes while conserving body water, electrolytes, and metabolites. Since sodium accounts for over 90 percent of the cations in the extracellular fluid, and the number of cations is balanced by the number of anions, considering the renal handling sodium and water only should sufficiently describe the relationship between the plasma compartment and kidneys. A kidney function model is presented which has been adapted from a

Slugging Flow of Water Draining from the Bottom of a Non-Vented Container

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

Charles W. Solbrig

2010-06-01

Experiments were run to observe the behavior of water exiting through an orifice at the bottom of an non-vented container. Initially, the container is nearly full of water with a small air space on top. Once the orifice was uncovered, the slugging rate and the drain rate of the water leaving the container were measured. Upon initially opening the orifice, water drains out until the air pressure above the water reduces enough that the air pressure drop from inside to outside of the container supports the water column and the water stops flowing. Air then enters the container through themore » orifice forming a bubble, which grows until it detaches and bubbles through the water to reach the air space. Once the bubble enters, this added air increases the pressure in the air space enough to allow the water to start flowing out again. This cycle of flow out, flow stoppage, air inflow, and bubble breakoff continues over and over until the hole is closed or the container empties. This is referred to as the “slugging cycle.” A mechanism is proposed to describe the slugging cycle which is modeled analytically. This paper presents the description of the experiments, data obtained, the mechanistic model, and comparison of the model to the experimental data. The model predicts outflow rates close to experimental values. Flow rates from non-vented containers are more than 10 to 20 less than vented containers. The bubbles which must enter the container periodically to increase the internal air pressure stop the water flow momentarily so are responsible for this large decrease in flow rate. Swirl induced in the non-vented container causes the flow rates to increase by a factor of two. The flow rate out of a non-vented container is independent of water height which is in direct contrast to a vented container where the flow rate is proportional to the square root of the water height. The constant rate is due to the container pressure. The higher the water level, the lower the

Application of nonlinear least-squares regression to ground-water flow modeling, west-central Florida

USGS Publications Warehouse

Yobbi, D.K.

2000-01-01

A nonlinear least-squares regression technique for estimation of ground-water flow model parameters was applied to an existing model of the regional aquifer system underlying west-central Florida. The regression technique minimizes the differences between measured and simulated water levels. Regression statistics, including parameter sensitivities and correlations, were calculated for reported parameter values in the existing model. Optimal parameter values for selected hydrologic variables of interest are estimated by nonlinear regression. Optimal estimates of parameter values are about 140 times greater than and about 0.01 times less than reported values. Independently estimating all parameters by nonlinear regression was impossible, given the existing zonation structure and number of observations, because of parameter insensitivity and correlation. Although the model yields parameter values similar to those estimated by other methods and reproduces the measured water levels reasonably accurately, a simpler parameter structure should be considered. Some possible ways of improving model calibration are to: (1) modify the defined parameter-zonation structure by omitting and/or combining parameters to be estimated; (2) carefully eliminate observation data based on evidence that they are likely to be biased; (3) collect additional water-level data; (4) assign values to insensitive parameters, and (5) estimate the most sensitive parameters first, then, using the optimized values for these parameters, estimate the entire data set.

Functional test of pedotransfer functions to predict water flow and solute transport with the dual-permeability model MACRO

NASA Astrophysics Data System (ADS)

Moeys, J.; Larsbo, M.; Bergström, L.; Brown, C. D.; Coquet, Y.; Jarvis, N. J.

2012-07-01

Estimating pesticide leaching risks at the regional scale requires the ability to completely parameterise a pesticide fate model using only survey data, such as soil and land-use maps. Such parameterisations usually rely on a set of lookup tables and (pedo)transfer functions, relating elementary soil and site properties to model parameters. The aim of this paper is to describe and test a complete set of parameter estimation algorithms developed for the pesticide fate model MACRO, which accounts for preferential flow in soil macropores. We used tracer monitoring data from 16 lysimeter studies, carried out in three European countries, to evaluate the ability of MACRO and this "blind parameterisation" scheme to reproduce measured solute leaching at the base of each lysimeter. We focused on the prediction of early tracer breakthrough due to preferential flow, because this is critical for pesticide leaching. We then calibrated a selected number of parameters in order to assess to what extent the prediction of water and solute leaching could be improved. Our results show that water flow was generally reasonably well predicted (median model efficiency, ME, of 0.42). Although the general pattern of solute leaching was reproduced well by the model, the overall model efficiency was low (median ME = -0.26) due to errors in the timing and magnitude of some peaks. Preferential solute leaching at early pore volumes was also systematically underestimated. Nonetheless, the ranking of soils according to solute loads at early pore volumes was reasonably well estimated (concordance correlation coefficient, CCC, between 0.54 and 0.72). Moreover, we also found that ignoring macropore flow leads to a significant deterioration in the ability of the model to reproduce the observed leaching pattern, and especially the early breakthrough in some soils. Finally, the calibration procedure showed that improving the estimation of solute transport parameters is probably more important than the

A Simple Model for the Viscosity of Rhyolites as a Function of Temperature, Pressure and Water Content: Implications for Obsidian Flow Emplacement

NASA Astrophysics Data System (ADS)

Whittington, A. G.; Romine, W. L.

2014-12-01

Understanding the dynamics of rhyolitic conduits and lava flows, requires precise knowledge of how viscosity (η) varies with temperature (T), pressure (P) and volatile content (X). In order to address the paucity of viscosity data for high-silica rhyolite at low water contents, which represent water saturation at near-surface conditions, we made 245 viscosity measurements on Mono Craters (California) rhyolites containing between 0.01 and 1.1 wt.% H2O, at temperatures between 796 and 1774 K using parallel plate and concentric cylinder methods at atmospheric pressure. We then developed and calibrated a new empirical model for the log of the viscosity of rhyolitic melts, where non-linear variations due to temperature and water content are nested within a linear dependence of log η on P. The model was fitted to a total of 563 data points: our 245 new data, 255 published data from rhyolites across a wide P-T-X space, and 63 data on haplogranitic and granitic melts under high P-T conditions. Statistically insignificant parameters were eliminated from the model in an effort to increase parsimony and the final model is simple enough for use in numerical models of conduit or lava flow dynamics: log η = -5.142+(13080-2982log⁡(w+0.229))/(T-(98.9-175.9 log⁡(w+0.229)))- P(0.0007-0.76/T ) where η is in Pa s, w is water content in wt.%, P is in MPa and T is in K. The root mean square deviation (rmsd) between the model predictions and the 563 data points used in calibration is 0.39 log units. Experimental constraints have led previously to spurious correlations between P, T, X and η in viscosity data sets, so that predictive models may struggle to correctly resolve the individual effects of P, T and X, and especially their cross-correlations. The increasing water solubility with depth inside a simple isothermal sheet of obsidian suggests that viscosity should decrease by ~1 order of magnitude at ~20m depth and by ~2 orders of magnitude at ~100m depth. If equilibrium water

The shallow water equations as a hybrid flow model for the numerical and experimental analysis of hydro power stations

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

Ostermann, Lars; Seidel, Christian

2015-03-10

The numerical analysis of hydro power stations is an important method of the hydraulic design and is used for the development and optimisation of hydro power stations in addition to the experiments with the physical submodel of a full model in the hydraulic laboratory. For the numerical analysis, 2D and 3D models are appropriate and commonly used.The 2D models refer mainly to the shallow water equations (SWE), since for this flow model a large experience on a wide field of applications for the flow analysis of numerous problems in hydraulic engineering already exists. Often, the flow model is verified bymore » in situ measurements. In order to consider 3D flow phenomena close to singularities like weirs, hydro power stations etc. the development of a hybrid fluid model is advantageous to improve the quality and significance of the global model. Here, an extended hybrid flow model based on the principle of the SWE is presented. The hybrid flow model directly links the numerical model with the experimental data, which may originate from physical full models, physical submodels and in-situ measurements. Hence a wide field of application of the hybrid model emerges including the improvement of numerical models and the strong coupling of numerical and experimental analysis.« less

Water Flow in the High Plains Aquifer in Northwestern Oklahoma

USGS Publications Warehouse

Luckey, Richard R.; Osborn, Noel I.; Becker, Mark F.; Andrews, William J.

2000-01-01

The High Plains is a major agricultural area, supported primarily by water from the High Plains aquifer, which is used to irrigate wheat and corn and to raise cattle and swine. The U.S. Geological Survey (USGS) and the Oklahoma Water Resources Board (OWRB) began a study of the High Plains aquifer in 1996. One purpose of the study was to develop a ground-water flow model that the OWRB could use to allocate the amount of water withdrawn from the a aquifer. The study area in Oklahoma covers all or parts of Beaver, Cimarron, Dewey, Ellis, Harper, Texas, and Woodward Counties. To provide appropriate hydrologic boundaries for the ground-water flow model, the study area was expanded to include parts of Colorado, Kansas, New Mexico, and Texas.

Ground-water flow in the shallow aquifer system at the Naval Weapons Station Yorktown, Virginia

USGS Publications Warehouse

Smith, Barry S.

2001-01-01

The Environmental Directorate of the Naval Weapons Station Yorktown, Virginia, is concerned about possible contamination of ground water at the Station. Ground water at the Station flows through a shallow system of layered aquifers and leaky confining units. The units of the shallow aquifer system are the Columbia aquifer, the Cornwallis Cave confining unit, the Cornwallis Cave aquifer, the Yorktown confining unit, and the Yorktown-Eastover aquifer. The Eastover-Calvert confining unit separates the shallow aquifer system from deeper confined aquifers beneath the Station. A three-dimensional, finite-difference, ground-water flow model was used to simulate steady-state ground-water flow of the shallow aquifer system in and around the Station. The model simulated ground-water flow from the peninsular drainage divide that runs across the Lackey Plain near the southern end of the Station north to King Creek and the York River and south to Skiffes Creek and the James River. The model was calibrated by minimizing the root mean square error between 4 7 measured and corresponding simulated water levels. The calibrated model was used to determine the ground-water budget and general directions of ground-water flow. A particle-tracking routine was used with the calibrated model to estimate groundwater flow paths, flow rates, and traveltimes from selected sites at the Station. Simulated ground-water flow velocities of the Station-area model were small beneath the interstream areas of the Lackey Plain and Croaker Flat, but increased outward toward the streams and rivers where the hydraulic gradients are larger. If contaminants from the land surface entered the water table at or near the interstream areas of the Station, where hydraulic gradients are smaller, they would migrate more slowly than if they entered closer to the streams or the shores of the rivers where gradients commonly are larger. The ground-water flow simulations indicate that some ground water leaks downward from

Description and comparison of selected models for hydrologic analysis of ground-water flow, St Joseph River basin, Indiana

USGS Publications Warehouse

Peters, J.G.

1987-01-01

The Indiana Department of Natural Resources (IDNR) is developing water-management policies designed to assess the effects of irrigation and other water uses on water supply in the basin. In support of this effort, the USGS, in cooperation with IDNR, began a study to evaluate appropriate methods for analyzing the effects of pumping on ground-water levels and streamflow in the basin 's glacial aquifer systems. Four analytical models describe drawdown for a nonleaky, confined aquifer and fully penetrating well; a leaky, confined aquifer and fully penetrating well; a leaky, confined aquifer and partially penetrating well; and an unconfined aquifer and partially penetrating well. Analytical equations, simplifying assumptions, and methods of application are described for each model. In addition to these four models, several other analytical models were used to predict the effects of ground-water pumping on water levels in the aquifer and on streamflow in local areas with up to two pumping wells. Analytical models for a variety of other hydrogeologic conditions are cited. A digital ground-water flow model was used to describe how a numerical model can be applied to a glacial aquifer system. The numerical model was used to predict the effects of six pumping plans in 46.5 sq mi area with as many as 150 wells. Water budgets for the six pumping plans were used to estimate the effect of pumping on streamflow reduction. Results of the analytical and numerical models indicate that, in general, the glacial aquifers in the basin are highly permeable. Radial hydraulic conductivity calculated by the analytical models ranged from 280 to 600 ft/day, compared to 210 and 360 ft/day used in the numerical model. Maximum seasonal pumping for irrigation produced maximum calculated drawdown of only one-fourth of available drawdown and reduced streamflow by as much as 21%. Analytical models are useful in estimating aquifer properties and predicting local effects of pumping in areas with

Numerical and Experimental Modeling of the Recirculating Melt Flow Inside an Induction Crucible Furnace

NASA Astrophysics Data System (ADS)

Asad, Amjad; Bauer, Katrin; Chattopadhyay, Kinnor; Schwarze, Rüdiger

2018-06-01

In the paper, a new water model of the turbulent recirculating flow in an induction furnace is introduced. The water model was based on the principle of the stirred vessel used in process engineering. The flow field in the water model was measured by means of particle image velocimetry in order to verify the model's performance. Here, it is indicated that the flow consists of two toroidal vortices similar to the flow in the induction crucible furnace. Furthermore, the turbulent flow in the water model is investigated numerically by adopting eddy-resolving turbulence modeling. The two toroidal vortices occur in the simulations as well. The numerical approaches provide identical time-averaged flow patterns. Moreover, a good qualitative agreement is observed on comparing the experimental and numerical results. In addition, a numerical simulation of the melt flow in a real induction crucible furnace was performed. The turbulent kinetic energy spectrum of the flow in the water model was compared to that of the melt flow in the induction crucible furnace to show the similarity in the nature of turbulence.

Analysis of internal flow characteristics of a smooth-disk water-brake dynamometer

NASA Technical Reports Server (NTRS)

Evans, D. G.

1973-01-01

The principal of absorbing power with an enclosed partially submerged rotating disk through the turbulent viscous shearing of water is discussed. Reference information is used to develop a flow model of the water brake. A method is then presented that uses vector diagrams to relate the effects of rotational flow, through flow, and secondary flow to power absorption. The method is used to describe the operating characteristics of an example 111-cm (43.7-in.) diameter water brake. Correlating performance parameters are developed in a dimensional analysis.

Hydrology and simulation of ground-water flow, Lake Point, Tooele County, Utah

USGS Publications Warehouse

Brooks, Lynette E.

2006-01-01

Water for new residential development in Lake Point, Utah may be supplied by public-supply wells completed in consolidated rock on the east side of Lake Point. Ground-water flow models were developed to help understand the effect the proposed withdrawal will have on water levels, flowing-well discharge, spring discharge, and ground-water quality in the study area. This report documents the conceptual and numerical ground-water flow models for the Lake Point area.The ground-water system in the Lake Point area receives recharge from local precipitation and irrigation, and from ground-water inflow from southwest of the area. Ground water discharges mostly to springs. Discharge also occurs to evapotranspiration, wells, and Great Salt Lake. Even though ground water discharges to Great Salt Lake, dense salt water from the lake intrudes under the less-dense ground water and forms a salt-water wedge under the valley. This salt water is responsible for some of the high dissolved-solids concentrations measured in ground water in Lake Point.A steady-state MODFLOW-2000 ground-water model of Tooele Valley adequately simulates water levels, ground-water discharge, and ground-water flow direction observed in Lake Point in 1969 and 2002. Simulating an additional 1,650 acre-feet per year withdrawal from wells causes a maximum projected drawdown of about 550 feet in consolidated rock near the simulated wells and drawdown exceeding 80 feet in an area encompassing most of the Oquirrh Mountains east of Lake Point. Drawdown in most of Lake Point ranges from 2 to 10 ft, but increases to more than 40 feet in the areas proposed for residential development. Discharge to Factory Springs, flowing wells, evapotranspiration, and Great Salt Lake is decreased by about 1,100 acre-feet per year (23 percent).The U.S. Geological Survey SUTRA variable-density ground-water-flow model generates a reasonable approximation of 2002 dissolved-solids concentration when simulating 2002 withdrawals. At most

Calculating ground water transit time of horizontal flow through leaky aquifers.

PubMed

Braunsfurth, Angelika C; Schneider, Wilfried

2008-01-01

The calculation of ground water transit times is one important factor in ground water protection. In this paper, we present an analytical solution for the transit time for a Dupuit-type flow system applicable to saturated flow through a horizontal leaky aquifer discharging to a downgradient fixed-head boundary under steady-state conditions. We investigate the influence of leakage when comparing the resulting travel times of our model based on head-dependent leakage with the commonly used model with no leakage and a simplified model with constant leakage. The results show significant differences in the position of the water divide and transit time, suggesting that leakage cannot be ignored.

What maintains the waters flowing in our rivers?

NASA Astrophysics Data System (ADS)

Vasconcelos, Vitor Vieira

2017-07-01

This article discusses how new contributions from hydrogeological science in the 20th and 21st centuries have allowed for a better understanding of the processes that affect the maintenance of river flows. Moreover, the way in which this knowledge has been conveyed beyond academia and has been gradually incorporated into public policy for natural resource management is also discussed. This article explains the development of several approaches used to understand the relationships among the management of aquifers, vegetation and river flows, including water balance, aquifer recharge, the piston effect, seasonal effects, and safe and sustainable yields. Additionally, the current challenges regarding the modeling of hydrological processes that integrate groundwater and surface waters are discussed. Examples of studies applied in Brazil that demonstrate these processes and stimulate thought regarding water management strategies are presented. In light of the case studies, it is possible to propose different strategies, each adapted for specific hydrogeological context to maximize aquifer recharge or base flow maintenance. Based on these strategies, the role of infiltration ponds and other artificial recharge techniques is re-evaluated in the context of the mitigation of environmental impacts on the maintenance of river flows. Proposals for the improvement of public policies regarding the payment of related environmental services to stimulate investment in aquifer recharge and the maintenance of base flow, for which the goal is to attain win-win-win situations for the environment, farmers and water users, while preventing land speculation, are discussed. Lastly, a conceptual model for the dissemination of hydrogeological knowledge in public policies is provided, and its challenges and possibilities are discussed.

Assessing Receiving Water Quality Impacts due to Flow Path Alteration in Residential Catchments, using the Stormwater and Wastewater Management Model

NASA Astrophysics Data System (ADS)

Wolosoff, S. E.; Duncan, J.; Endreny, T.

2001-05-01

The Croton water supply system, responsible for supplying approximately 10% of New York City's water, provides an opportunity for exploration into the impacts of significant terrestrial flow path alteration upon receiving water quality. Natural flow paths are altered during residential development in order to allow for construction at a given location, reductions in water table elevation in low lying areas and to provide drainage of increased overland flow volumes. Runoff conducted through an artificial drainage system, is prevented from being attenuated by the natural environment, thus the pollutant removal capacity inherent in most natural catchments is often limited to areas where flow paths are not altered by development. By contrasting the impacts of flow path alterations in two small catchments in the Croton system, with different densities of residential development, we can begin to identify appropriate limits to the re-routing of runoff in catchments draining into surface water supplies. The Stormwater and Wastewater Management Model (SWMM) will be used as a tool to predict the runoff quantity and quality generated from two small residential catchments and to simulate the potential benefits of changes to the existing drainage system design, which may improve water quality due to longer residence times.

FUNDAMENTALS OF GROUND-WATER MODELING

EPA Science Inventory

Ground-water flow and contaminant transport modeling has been used at many hazardous waste sites with varying degrees of success. odels may be used throughout all phases of the site investigation and remediation processes. eveloping a better understanding of ground-water modeling...

Two-dimensional, steady-state model of ground-water flow, Nevada Test Site and vicinity, Nevada-California

USGS Publications Warehouse

Waddell, R.K.

1982-01-01

A two-dimensional, steady-state model of ground-water flow beneath the Nevada Test Site and vicinity has been developed using inverse techniques. The area is underlain by clastic and carbonate rocks of Precambrian and Paleozoic age and by volcanic rocks and alluvium of Tertiary and Quaternary age that have been juxtaposed by normal and strike-slip faulting. Aquifers are composed of carbonate and volcanic rocks and alluvium. Characteristics of the flow system are determined by distribution of low-conductivity rocks (barriers); by recharge originating in the Spring Mountains, Pahranagat, Timpahute, and Sheep Ranges, and in Pahute Mesa; and by underflow beneath Pahute Mesa from Gold Flat and Kawich Valley. Discharge areas (Ash Meadows, Oasis Valley, Alkali Flat, and Furnace Creek Ranch) are upgradient from barriers. Sensitivities of simulated hydraulic heads and fluxes to variations in model parameters were calculated to guide field studies and to help estimate errors in predictions from transport modeling. Hydraulic heads and fluxes are very sensitive to variations in the greater magnitude recharge/discharge terms. Transmissivity at a location may not be the most important transmissivity for determining flux there. Transmissivities and geometries of large barriers that impede flow from Pahute Mesa have major effects on fluxes elsewhere; as their transmissivities are decreased, flux beneath western Jackass Flats and Yucca Mountains is increased as water is diverted around the barriers. Fortymile Canyon is underlain by highly transmissive rocks that cause potentiometric contours to vee upgradient; increasing their transmissivity increases flow through them, and decreases it beneath Yucca Mountain. (USGS)

Two-dimensional ground-water flow model of the Cretaceous aquifer system of Lee County and vicinity, Mississippi

USGS Publications Warehouse

Kernodle, John Michael

1981-01-01

A two-dimensional ground-water flow model of the Eutaw-McShan and Gordo aquifers in the area of Lee County, Miss., was successfully calibrated and verified using data from six long-term observation wells and two intensive studies of areal water levels. The water levels computed by the model were found to be most sensitive to changes in simulated aquifer hydraulic conductivity and to changes in head in the overlying Coffee Sand aquifer. The two-dimensional model performed reasonably well in simulating the aquifer system except possibly in southern Lee County and southward where a clay bed at the top of the Gordo Formation partially isolated the Gordo from the overlying Eutaw-McShan aquifer. The verified model was used to determine theoretical aquifer response to increased ground-water withdrawal to the year 2000. Two estimated rates of increase and five possible well field locations were examined. (USGS)

A theory for modeling ground-water flow in heterogeneous media

USGS Publications Warehouse

Cooley, Richard L.

2004-01-01

Construction of a ground-water model for a field area is not a straightforward process. Data are virtually never complete or detailed enough to allow substitution into the model equations and direct computation of the results of interest. Formal model calibration through optimization, statistical, and geostatistical methods is being applied to an increasing extent to deal with this problem and provide for quantitative evaluation and uncertainty analysis of the model. However, these approaches are hampered by two pervasive problems: 1) nonlinearity of the solution of the model equations with respect to some of the model (or hydrogeologic) input variables (termed in this report system characteristics) and 2) detailed and generally unknown spatial variability (heterogeneity) of some of the system characteristics such as log hydraulic conductivity, specific storage, recharge and discharge, and boundary conditions. A theory is developed in this report to address these problems. The theory allows construction and analysis of a ground-water model of flow (and, by extension, transport) in heterogeneous media using a small number of lumped or smoothed system characteristics (termed parameters). The theory fully addresses both nonlinearity and heterogeneity in such a way that the parameters are not assumed to be effective values. The ground-water flow system is assumed to be adequately characterized by a set of spatially and temporally distributed discrete values, ?, of the system characteristics. This set contains both small-scale variability that cannot be described in a model and large-scale variability that can. The spatial and temporal variability in ? are accounted for by imagining ? to be generated by a stochastic process wherein ? is normally distributed, although normality is not essential. Because ? has too large a dimension to be estimated using the data normally available, for modeling purposes ? is replaced by a smoothed or lumped approximation y?. (where y is a

Geohydrology and Numerical Simulation of the Ground-Water Flow System of Molokai, Hawaii

USGS Publications Warehouse

Oki, Delwyn S.

1997-01-01

A two-dimensional, steady-state, areal ground-water flow model was developed for the island of Molokai, Hawaii, to enhance the understanding of (1) the conceptual framework of the ground-water flow system, (2) the distribution of aquifer hydraulic properties, and (3) the regional effects of ground-water withdrawals on water levels and coastal discharge. The model uses the finite-element code AQUIFEM-SALT, which simulates flow of fresh ground water in systems that may have a freshwater lens floating on denser underlying saltwater. Model results are in agreement with the general conceptual model of the flow system on Molokai, where ground water flows from the interior, high-recharge areas to the coast. The model-calculated ground-water divide separating flow to the northern and southern coasts lies to either the north or the south of the topographic divide but is generally not coincident with the topographic divide. On the basis of model results, the following horizontal hydraulic conductivities were estimated: (1) 1,000 feet per day for the dike-free volcanic rocks of East and West Molokai, (2) 100 feet per day for the marginal dike zone of the East Molokai Volcano, (3) 2 feet per day for the West Molokai dike complex, (4) 0.02 feet per day for the East Molokai dike complex, and (5) 500 feet per day for the Kalaupapa Volcanics. Three simulations to determine the effects of proposed ground-water withdrawals on water levels and coastal discharge, relative to model-calculated water levels and coastal discharge for 1992-96 withdrawal rates, show that the effects are widespread. For a withdrawal rate of 0.337 million gallons per day from a proposed well about 4 miles southeast of Kualapuu and 3 miles north of Kamiloloa, the model-calculated drawdown of 0.01 foot or more extends 4 miles southeast and 6 miles northwest from the well. For a withdrawal rate of 1.326 million gallons per day from the same well, the model-calculated drawdown of 0.01 foot or more extends 6 miles

A characterisation model to address the environmental impact of green water flows for water scarcity footprints.

PubMed

Quinteiro, Paula; Rafael, Sandra; Villanueva-Rey, Pedro; Ridoutt, Bradley; Lopes, Myriam; Arroja, Luís; Dias, Ana Cláudia

2018-06-01

The development of methods to assess the potential environmental impact of green water consumption in life cycle assessment has lagged behind those for blue water use, which are now routinely applied in industrial and policy-related studies. This represents a critical gap in the assessment of land-based production systems and the ability to inform policy related to the bio-economy. Combining satellite remote sensing and meteorological data sets, this study develops two new sets of spatially-differentiated and globally applicable characterisation factors (CFs) to assess the environmental impact of green water flows in LCA. One set of CFs addresses the impact of shifts in water vapour flow by evapotranspiration on blue water availability (CFWS) and the other set of CFs addresses moisture recycling within a basin (CFWA). Furthermore, as an additional and optional step, these two indicators are combined into an aggregated green water scarcity indicator, representing the global variability of green water scarcity. The values obtained for CFWA show that there are significant changes in green water flows that were returned to the atmosphere in Alaska (covered by open shrublands) and in some central regions of China (covered by grasslands and barren or sparsely vegetated land), where precipitation levels are lower than 10 mm/yr. The results obtained for CFWS indicate that severe perturbations in surface blue water production occur, particularly in central regions of China (covered by grasslands), the southeast of Australia (covered by evergreen broadleaf forest) and in some central regions of the USA (covered by grassland and evergreen needleleaf forest). The application of the green water scarcity CFs enables the evaluation of the potential environmental impact due to green water consumption by agricultural and forestry products, informing both technical and non-technical audiences and decision-makers for the purpose of strategic planning of land use and to identify green

South Asia river flow projections and their implications for water resources

NASA Astrophysics Data System (ADS)

Mathison, C.; Wiltshire, A. J.; Falloon, P.; Challinor, A. J.

2015-06-01

South Asia is a region with a large and rising population and a high dependance on industries sensitive to water resource such as agriculture. The climate is hugely variable with the region relying on both the Asian Summer Monsoon (ASM) and glaciers for its supply of fresh water. In recent years, changes in the ASM, fears over the rapid retreat of glaciers and the increasing demand for water resources for domestic and industrial use, have caused concern over the reliability of water resources both in the present day and future for this region. The climate of South Asia means it is one of the most irrigated agricultural regions in the world, therefore pressures on water resource affecting the availability of water for irrigation could adversely affect crop yields and therefore food production. In this paper we present the first 25 km resolution regional climate projections of river flow for the South Asia region. ERA-Interim, together with two global climate models (GCMs), which represent the present day processes, particularly the monsoon, reasonably well are downscaled using a regional climate model (RCM) for the periods; 1990-2006 for ERA-Interim and 1960-2100 for the two GCMs. The RCM river flow is routed using a river-routing model to allow analysis of present day and future river flows through comparison with river gauge observations, where available. In this analysis we compare the river flow rate for 12 gauges selected to represent the largest river basins for this region; Ganges, Indus and Brahmaputra basins and characterize the changing conditions from east to west across the Himalayan arc. Observations of precipitation and runoff in this region have large or unknown uncertainties, are short in length or are outside the simulation period, hindering model development and validation designed to improve understanding of the water cycle for this region. In the absence of robust observations for South Asia, a downscaled ERA-Interim RCM simulation provides a

Numerical model of water flow in a fractured basalt vadose zone: Box Canyon Site, Idaho

NASA Astrophysics Data System (ADS)

Doughty, Christine

2000-12-01

A numerical model of a fractured basalt vadose zone has been developed on the basis of the conceptual model described by Faybishenko et al. [[his issue]. The model has been used to simulate a ponded infiltration test in order to investigate infiltration through partially saturated fractured basalt. A key question addressed is how the fracture pattern geometry and fracture connectivity within a single basalt flow of the Snake River Plain basalt affect water infiltration. The two-dimensional numerical model extends from the ground surface to a perched water body 20 m below and uses an unconventional quasi-deterministic approach with explicit but highly simplified representation of major fractures and other important hydrogeologic features. The model adequately reproduces the majority of the field observation and provides insights into the infiltration process that cannot be obtained by data collection alone, demonstrating its value as a component of field studies.

Water flow in high-speed handpieces.

PubMed

Cavalcanti, Bruno Neves; Serairdarian, Paulo Isaías; Rode, Sigmar Mello

2005-05-01

This study measured the water flow commonly used in high-speed handpieces to evaluate the water flow's influence on temperature generation. Different flow speeds were evaluated between turbines that had different numbers of cooling apertures. Two water samples were collected from each high-speed handpiece at private practices and at the School of Dentistry at São José dos Campos. The first sample was collected at the customary flow and the second was collected with the terminal opened for maximum flow. The two samples were collected into weighed glass receptacles after 15 seconds of turbine operation. The glass receptacles were reweighed and the difference between weights was recorded to calculate the water flow in mL/min and for further statistical analysis. The average water flow for 137 samples was 29.48 mL/min. The flow speeds obtained were 42.38 mL/min for turbines with one coolant aperture; 34.31 mL/min for turbines with two coolant apertures; and 30.44 mL/min for turbines with three coolant apertures. There were statistical differences between turbines with one and three coolant apertures (Tukey-Kramer multiple comparisons test with P < .05). Turbine handpieces with one cooling aperture distributed more water for the burs than high-speed handpieces with more than one aperture.

A Device to Emulate Diffusion and Thermal Conductivity Using Water Flow

ERIC Educational Resources Information Center

Blanck, Harvey F.

2005-01-01

A device designed to emulate diffusion and thermal conductivity using flowing water is reviewed. Water flowing through a series of cells connected by a small tube in each partition in this plastic model is capable of emulating diffusion and thermal conductivity that occurs in variety of systems described by several mathematical equations.

Comparison of an algebraic multigrid algorithm to two iterative solvers used for modeling ground water flow and transport

USGS Publications Warehouse

Detwiler, R.L.; Mehl, S.; Rajaram, H.; Cheung, W.W.

2002-01-01

Numerical solution of large-scale ground water flow and transport problems is often constrained by the convergence behavior of the iterative solvers used to solve the resulting systems of equations. We demonstrate the ability of an algebraic multigrid algorithm (AMG) to efficiently solve the large, sparse systems of equations that result from computational models of ground water flow and transport in large and complex domains. Unlike geometric multigrid methods, this algorithm is applicable to problems in complex flow geometries, such as those encountered in pore-scale modeling of two-phase flow and transport. We integrated AMG into MODFLOW 2000 to compare two- and three-dimensional flow simulations using AMG to simulations using PCG2, a preconditioned conjugate gradient solver that uses the modified incomplete Cholesky preconditioner and is included with MODFLOW 2000. CPU times required for convergence with AMG were up to 140 times faster than those for PCG2. The cost of this increased speed was up to a nine-fold increase in required random access memory (RAM) for the three-dimensional problems and up to a four-fold increase in required RAM for the two-dimensional problems. We also compared two-dimensional numerical simulations of steady-state transport using AMG and the generalized minimum residual method with an incomplete LU-decomposition preconditioner. For these transport simulations, AMG yielded increased speeds of up to 17 times with only a 20% increase in required RAM. The ability of AMG to solve flow and transport problems in large, complex flow systems and its ready availability make it an ideal solver for use in both field-scale and pore-scale modeling.

Documentation of a computer program to simulate aquifer-system compaction using the modular finite-difference ground-water flow model

USGS Publications Warehouse

Leake, S.A.; Prudic, David E.

1991-01-01

Removal of ground water by pumping from aquifers may result in compaction of compressible fine-grained beds that are within or adjacent to the aquifers. Compaction of the sediments and resulting land subsidence may be permanent if the head declines result in vertical stresses beyond the previous maximum stress. The process of permanent compaction is not routinely included in simulations of ground-water flow. To simulate storage changes from both elastic and inelastic compaction, a computer program was written for use with the U.S. Geological Survey modular finite-difference ground- water flow model. The new program, the Interbed-Storage Package, is designed to be incorporated into this model. In the Interbed-Storage Package, elastic compaction or expansion is assumed to be proportional to change in head. The constant of proportionality is the product of the skeletal component of elastic specific storage and the thickness of the sediments. Similarly, inelastic compaction is assumed to be proportional to decline in head. The constant of proportionality is the product of the skeletal component of inelastic specific storage and the thickness of the sediments. Storage changes are incorporated into the ground-water flow model by adding an additional term to the right-hand side of the flow equation. Within a model time step, the package appropriately apportions storage changes between elastic and inelastic components on the basis of the relation of simulated head to the previous minimum (preconsolidation) head. Two tests were performed to verify that the package works correctly. The first test compared model-calculated storage and compaction changes to hand-calculated values for a three-dimensional simulation. Model and hand-calculated values were essentially equal. The second test was performed to compare the results of the Interbed-Storage Package with results of the one-dimensional Helm compaction model. This test problem simulated compaction in doubly draining

The origin of high and low flows in the river Rhine: particle tracing and water quality calculations in a distributed hydrological model

NASA Astrophysics Data System (ADS)

Schellekens, Jaap; van Gils, Jos; Christophe, Christophe; Sperna-Weiland, Frederiek; Winsemius, Hessel

2013-04-01

The ability to quickly link a complete water quality model to any distributed hydrological model can be of great value. It provides the hydrological modeller with more information on the performance of the model by being able to add particle tracing and independent mass balance calculations to an existing distributed hydrological model. It also allows for full catchment water quality calculations forced by emissions to different hydrological compartments, taking into account the relevant processes in the different compartments of the hydrological model. A combined distributed hydrological model and hydrochemical model (Delwaq) have been combined within the modeling framework OpenStreams to model large scale hydrological processes in the Rhine basin upstream of the Dutch border at Lobith. Several models have been setup to evaluate (1) the origin of high and low flows in the Rhine basin based on subcatchment contribution and (2) the contribution of different land covers to the total flow with special reference to urban land cover. In addition (3) the relative share of fast and slow runoff components in the total river discharge has been quantified, as well as the age of these two fractions, both as a function of time. Finally (4) the transmission of a pollutant released in infiltrating water and undergoing sorption has been simulated, as a first test for implementing full water quality modelling. The results of a thirty-five year run using daily time steps for 1975 to 2010 were analysed for monthly average contribution to the total flow of each subcatchment and the different land cover types both for average flow conditions and for the top ten and bottom ten flow percentiles. Furthermore, a number of high and low flow events have been analysed in detail. They reveal the large contribution of the basin area upstream of Basel to the dry season flow, especially during the driest summers. Flood conditions in the basin have a more varied origin with the Moselle being the

Geohydrology of Storage Unit III and a combined flow model of the Santa Barbara and foothill ground-water basins, Santa Barbara County, California

USGS Publications Warehouse

Freckleton, John R.; Martin, Peter; Nishikawa, Tracy

1998-01-01

The city of Santa Barbara pumps most of its ground water from the Santa Barbara and Foothill ground-water basins. The Santa Barbara basin is subdivided into two storage units: Storage Unit I and Storage Unit III. The Foothill basin and Storage Unit I of the Santa Barbara basin have been studied extensively and ground-water flow models have been developed for them. In this report, the geohydrology of the Santa Barbara ground- water basin is described with a special emphasis on Storage Unit III in the southwestern part of the basin. The purposes of this study were to summarize and evaluate the geohydrology of Storage Unit III and to develop an areawide model of the Santa Barbara and Foothill basins that includes the previously unmodeled Storage Unit III. Storage Unit III is in the southwestern part of the city of Santa Barbara. It is approximately 3.5 miles long and varies in width from about 2,000 feet in the southeast to 4,000 feet in the north-west. Storage Unit III is composed of the Santa Barbara Formation and overlying alluvium. The Santa Barbara Formation (the principal aquifer) consists of Pleistocene and Pliocene(?) unconsolidated marine sand, silt, and clay, and it has a maximum saturated thickness of about 160 feet. The alluvium that overlies the Santa Barbara Formation has a maximum saturated thickness of about 140 feet. The storage unit is bounded areally by faults and low-permeability deposits and is underlain by rocks of Tertiary age. The main sources of recharge to Storage Unit III are seepage from Arroyo Burro and infiltration of precipitation. Most of the recharge occurs in the northwest part of the storage unit, and ground water flows toward the southeast along the unit's long axis. Lesser amounts of recharge may occur as subsurface flow from the Hope Ranch subbasin and as upwelling from the underlying Tertiary rocks. Discharge from Storage Unit III occurs as pumpage, flow to underground drains, underflow through alluvium in the vicinity of Arroyo

Simulation of ground-water flow in glaciofluvial aquifers in the Grand Rapids area, Minnesota

USGS Publications Warehouse

Jones, Perry M.

2004-01-01

A calibrated steady-state, finite-difference, ground-waterflow model was constructed to simulate ground-water flow in three glaciofluvial aquifers, defined in this report as the upper, middle, and lower aquifers, in an area of about 114 mi2 surrounding the city of Grand Rapids in north-central Minnesota. The calibrated model will be used by Minnesota Department of Health and communities in the Grand Rapids area in the development of wellhead protection plans for their water supplies. The model was calibrated through comparison of simulated ground-water levels to measured static water levels in 351 wells, and comparison of simulated base-flow rates to estimated base-flow rates for reaches of the Mississippi and Prairie Rivers. Model statistics indicate that the model tends to overestimate ground-water levels. The root mean square errors ranged from +12.83 ft in wells completed in the upper aquifer to +19.10 ft in wells completed in the middle aquifer. Mean absolute differences between simulated and measured water levels ranged from +4.43 ft for wells completed in the upper aquifer to +9.25 ft for wells completed in the middle aquifer. Mean algebraic differences ranged from +9.35 ft for wells completed in the upper aquifer to +14.44 ft for wells completed in the middle aquifer, with the positive differences indicating that the simulated water levels were higher than the measured water levels. Percentage errors between simulated and estimated base-flow rates for the three monitored reaches all were less than 10 percent, indicating good agreement. Simulated ground-water levels were most sensitive to changes in general-head boundary conductance, indicating that this characteristic is the predominant model input variable controlling steady-state water-level conditions. Simulated groundwater flow to stream reaches was most sensitive to changes in horizontal hydraulic conductivity, indicating that this characteristic is the predominant model input variable controlling

Root growth, water uptake, and sap flow of winter wheat in response to different soil water conditions

NASA Astrophysics Data System (ADS)

Cai, Gaochao; Vanderborght, Jan; Langensiepen, Matthias; Schnepf, Andrea; Hüging, Hubert; Vereecken, Harry

2018-04-01

How much water can be taken up by roots and how this depends on the root and water distributions in the root zone are important questions that need to be answered to describe water fluxes in the soil-plant-atmosphere system. Physically based root water uptake (RWU) models that relate RWU to transpiration, root density, and water potential distributions have been developed but used or tested far less. This study aims at evaluating the simulated RWU of winter wheat using the empirical Feddes-Jarvis (FJ) model and the physically based Couvreur (C) model for different soil water conditions and soil textures compared to sap flow measurements. Soil water content (SWC), water potential, and root development were monitored noninvasively at six soil depths in two rhizotron facilities that were constructed in two soil textures: stony vs. silty, with each of three water treatments: sheltered, rainfed, and irrigated. Soil and root parameters of the two models were derived from inverse modeling and simulated RWU was compared with sap flow measurements for validation. The different soil types and water treatments resulted in different crop biomass, root densities, and root distributions with depth. The two models simulated the lowest RWU in the sheltered plot of the stony soil where RWU was also lower than the potential RWU. In the silty soil, simulated RWU was equal to the potential uptake for all treatments. The variation of simulated RWU among the different plots agreed well with measured sap flow but the C model predicted the ratios of the transpiration fluxes in the two soil types slightly better than the FJ model. The root hydraulic parameters of the C model could be constrained by the field data but not the water stress parameters of the FJ model. This was attributed to differences in root densities between the different soils and treatments which are accounted for by the C model, whereas the FJ model only considers normalized root densities. The impact of differences in

Conceptual model for quantifying pre-smolt production from flow-dependent physical habitat and water temperature

USGS Publications Warehouse

Williamson, S. C.; Bartholow, J. M.; Stalnaker, C. B.

1993-01-01

A conceptual model has been developed to test river regulation concepts by linking physical habitat and water temperature with salmonid population and production in cold water streams. Work is in progress to examine numerous questions as part of flow evaluation and habitat restoration programmes in the Trinity River of California and elsewhere. For instance, how much change in pre-smolt chinook salmon (Oncorhynchus tshawytscha) production in the Trinity River would result from a different annual instream allocation (i.e. up or down from 271 × 106 m3released in the late 1980s) and how much change in pre-smolt production would result from a different release pattern (i.e. different from the 8.5 m3 s−1 year-round release). The conceptual model is being used to: design, integrate and improve young-of-year population data collection efforts; test hypotheses that physical habitat significantly influences movement, growth and mortality of salmonid fishes; and analyse the relative severity of limiting factors during each life stage. The conceptual model, in conjunction with previously developed tools in the Instream Flow Incremental Methodology, should provide the means to more effectively manage a fishery resource below a regulated reservoir and to provide positive feedback to planning of annual reservoir operations.

The physics of confined flow and its application to water leaks, water permeation and water nanoflows: a review.

PubMed

Lei, Wenwen; Rigozzi, Michelle K; McKenzie, David R

2016-02-01

This review assesses the current state of understanding of the calculation of the rate of flow of gases, vapours and liquids confined in channels, in porous media and in permeable materials with an emphasis on the flow of water and its vapour. One motivation is to investigate the relation between the permeation rate of moisture and that of a noncondensable test gas such as helium, another is to assist in unifying theory and experiment across disparate fields. Available theories of single component ideal gas flows in channels of defined geometry (cylindrical, rectangular and elliptical) are described and their predictions compared with measurement over a wide range of conditions defined by the Knudsen number. Theory for two phase flows is assembled in order to understand the behaviour of four standard water leak configurations: vapour, slug, Washburn and liquid flow, distinguished by the number and location of phase boundaries (menisci). Air may or may not be present as a background gas. Slip length is an important parameter that greatly affects leak rates. Measurements of water vapour flows confirm that water vapour shows ideal gas behaviour. Results on carbon nanotubes show that smooth walls may lead to anomalously high slip lengths arising from the properties of 'confined' water. In porous media, behaviour can be matched to the four standard leaks. Traditional membrane permeation models consider that the permeant dissolves, diffuses and evaporates at the outlet side, ideas we align with those from channel flow. Recent results on graphite oxide membranes show examples where helium which does not permeate while at the same time moisture is almost unimpeded, again a result of confined water. We conclude that while there is no a priori relation between a noncondensable gas flow and a moisture flow, measurements using helium will give results within two orders of magnitude of the moisture flow rate, except in the case where there is anomalous slip or confined water

Flow measurements in a water tunnel using a holocinematographic velocimeter

NASA Technical Reports Server (NTRS)

Weinstein, Leonard M.; Beeler, George B.

1987-01-01

Dual-view holographic movies were used to examine complex flows with full three-space and time resolution. This approach, which tracks the movement of small tracer particles in water, is termed holocinematographic velocimetry (HCV). A small prototype of a new water tunnel was used to demonstrate proof-of-concept for the HCV. After utilizing a conventional flow visualization apparatus with a laser light sheet to illuminate tracer particles to evaluate flow quality of the prototype tunnel, a simplified version of the HCV was employed to demonstrate the capabilities of the approach. Results indicate that a full-scale version of the water tunnel and a high performance version of the HCV should be able to check theoretical and numerical modeling of complex flows and examine the mechanisms operative in turbulent and vortex flow control concepts, providing an entirely unique instrument capable, for the first time, of simultaneous three-space and time measurements in turbulent flow.

Stochastic Ground Water Flow Simulation with a Fracture Zone Continuum Model

USGS Publications Warehouse

Langevin, C.D.

2003-01-01

A method is presented for incorporating the hydraulic effects of vertical fracture zones into two-dimensional cell-based continuum models of ground water flow and particle tracking. High hydraulic conductivity features are used in the model to represent fracture zones. For fracture zones that are not coincident with model rows or columns, an adjustment is required for the hydraulic conductivity value entered into the model cells to compensate for the longer flowpath through the model grid. A similar adjustment is also required for simulated travel times through model cells. A travel time error of less than 8% can occur for particles moving through fractures with certain orientations. The fracture zone continuum model uses stochastically generated fracture zone networks and Monte Carlo analysis to quantify uncertainties with simulated advective travel times. An approach is also presented for converting an equivalent continuum model into a fracture zone continuum model by establishing the contribution of matrix block transmissivity to the bulk transmissivity of the aquifer. The methods are used for a case study in west-central Florida to quantify advective travel times from a potential wetland rehydration site to a municipal supply wellfield. Uncertainties in advective travel times are assumed to result from the presence of vertical fracture zones, commonly observed on aerial photographs as photolineaments.

Effects of spatially distributed sectoral water management on the redistribution of water resources in an integrated water model: SECTORAL WATER MANAGEMENT IN IA-ESM

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

Voisin, Nathalie; Hejazi, Mohamad I.; Leung, L. Ruby

To advance understanding of the interactions between human activities and the water cycle, an integrated terrestrial water cycle component has been developed for Earth system models. This includes a land surface model fully coupled to a river routing model and a generic water management model to simulate natural and regulated flows. A global integrated assessment model and its regionalized version for the U.S. are used to simulate water demand consistent with the energy technology and socio-economics scenarios. Human influence on the hydrologic cycle includes regulation and storage from reservoirs, consumptive use and withdrawal from multiple sectors ( irrigation and non-irrigation)more » and overall redistribution of water resources in space and time. As groundwater provides an important source of water supply for irrigation and other uses, the integrated modeling framework has been extended with a simplified representation of groundwater as an additional supply source, and return flow generated from differences between withdrawals and consumptive uses from both groundwater and surface water systems. The groundwater supply and return flow modules are evaluated by analyzing the simulated regulated flow, reservoir storage and supply deficit for irrigation and non-irrigation sectors over major hydrologic regions of the conterminous U.S. The modeling framework is then used to provide insights on the reliability of water resources by isolating the reliability due to return flow and/or groundwater sources of water. Our results show that high sectoral ratio of withdrawals over consumptive demand adds significant stress on the water resources management that can be alleviated by reservoir storage capacity. The return flow representation therefore exhibits a clear east-west contrast in its hydrologic signature, as well as in its ability to help meet water demand. Groundwater use has a limited hydrologic signature but the most pronounced signature is in terms of

Prediction of unsaturated flow and water backfill during infiltration in layered soils

NASA Astrophysics Data System (ADS)

Cui, Guotao; Zhu, Jianting

2018-02-01

We develop a new analytical infiltration model to determine water flow dynamics around layer interfaces during infiltration process in layered soils. The model mainly involves the analytical solutions to quadratic equations to determine the flux rates around the interfaces. Active water content profile behind the wetting front is developed based on the solution of steady state flow to dynamically update active parameters in sharp wetting front infiltration equations and to predict unsaturated flow in coarse layers before the front reaches an impeding fine layer. The effect of water backfill to saturate the coarse layers after the wetting front encounters the impeding fine layer is analytically expressed based on the active water content profiles. Comparison to the numerical solutions of the Richards equation shows that the new model can well capture water dynamics in relation to the arrangement of soil layers. The steady state active water content profile can be used to predict the saturation state of all layers when the wetting front first passes through these layers during the unsteady infiltration process. Water backfill effect may occur when the unsaturated wetting front encounters a fine layer underlying a coarse layer. Sensitivity analysis shows that saturated hydraulic conductivity is the parameter dictating the occurrence of unsaturated flow and water backfill and can be used to represent the coarseness of soil layers. Water backfill effect occurs in coarse layers between upper and lower fine layers when the lower layer is not significantly coarser than the upper layer.

Simulation of Ground-Water Flow and Effects of Ground-Water Irrigation on Base Flow in the Elkhorn and Loup River Basins, Nebraska

USGS Publications Warehouse

Peterson, Steven M.; Stanton, Jennifer S.; Saunders, Amanda T.; Bradley, Jesse R.

2008-01-01

Irrigated agriculture is vital to the livelihood of communities in the Elkhorn and Loup River Basins in Nebraska, and ground water is used to irrigate most of the cropland. Concerns about the sustainability of ground-water and surface-water resources have prompted State and regional agencies to evaluate the cumulative effects of ground-water irrigation in this area. To facilitate understanding of the effects of ground-water irrigation, a numerical computer model was developed to simulate ground-water flow and assess the effects of ground-water irrigation (including ground-water withdrawals, hereinafter referred to as pumpage, and enhanced recharge) on stream base flow. The study area covers approximately 30,800 square miles, and includes the Elkhorn River Basin upstream from Norfolk, Nebraska, and the Loup River Basin upstream from Columbus, Nebraska. The water-table aquifer consists of Quaternary-age sands and gravels and Tertiary-age silts, sands, and gravels. The simulation was constructed using one layer with 2-mile by 2-mile cell size. Simulations were constructed to represent the ground-water system before 1940 and from 1940 through 2005, and to simulate hypothetical conditions from 2006 through 2045 or 2055. The first simulation represents steady-state conditions of the system before anthropogenic effects, and then simulates the effects of early surface-water development activities and recharge of water leaking from canals during 1895 to 1940. The first simulation ends at 1940 because before that time, very little pumpage for irrigation occurred, but after that time it became increasingly commonplace. The pre-1940 simulation was calibrated against measured water levels and estimated long-term base flow, and the 1940 through 2005 simulation was calibrated against measured water-level changes and estimated long-term base flow. The calibrated 1940 through 2005 simulation was used as the basis for analyzing hypothetical scenarios to evaluate the effects of

ANALYSIS AND REDUCTION OF LANDSAT DATA FOR USE IN A HIGH PLAINS GROUND-WATER FLOW MODEL.

USGS Publications Warehouse

Thelin, Gail; Gaydas, Leonard; Donovan, Walter; Mladinich, Carol

1984-01-01

Data obtained from 59 Landsat scenes were used to estimate the areal extent of irrigated agriculture over the High Plains region of the United States for a ground-water flow model. This model provides information on current trends in the amount and distribution of water used for irrigation. The analysis and reduction process required that each Landsat scene be ratioed, interpreted, and aggregated. Data reduction by aggregation was an efficient technique for handling the volume of data analyzed. This process bypassed problems inherent in geometrically correcting and mosaicking the data at pixel resolution and combined the individual Landsat classification into one comprehensive data set.

A water tunnel flow visualization study of the vortex flow structures on the F/A-18 aircraft

NASA Technical Reports Server (NTRS)

Sandlin, Doral R.; Ramirez, Edgar J.

1991-01-01

The vortex flow structures occurring on the F/A-18 aircraft at high angles of attack were studied. A water tunnel was used to gather flow visualization data on the forebody vortex and the wing leading edge extension vortex. The longitudinal location of breakdown of the leading edge vortex was found to be consistently dependent on the angle of attack. Other parameters such as Reynolds number, model scale, and model fidelity had little influence on the overall behavior of the flow structures studied. The lateral location of the forebody vortex system was greatly influenced by changes in the angle of sideslip. Strong interactions can occur between the leading edge extension vortex and the forebody vortex. Close attention was paid to vortex induced flows on various airframe components of the F/A-18. Reynolds number and angle of attack greatly affected the swirling intensity, and therefore the strength of the studied vortices. Water tunnel results on the F/A-18 correlated well with those obtained in similar studies at both full and sub scale levels. The water tunnel can provide, under certain conditions, good simulations of realistic flows in full scale configurations.

Laminar and Turbulent Flow in Water

ERIC Educational Resources Information Center

Riveros, H. G.; Riveros-Rosas, D.

2010-01-01

There are many ways to visualize flow, either for laminar or turbulent flows. A very convincing way to show laminar and turbulent flows is by the perturbations on the surface of a beam of water coming out of a cylindrical tube. Photographs, taken with a flash, show the nature of the flow of water in pipes. They clearly show the difference between…

Eye wash water flow direction study: an evaluation of the effectiveness of eye wash devices with opposite directional water flow.

PubMed

Fogt, Jennifer S; Jones-Jordan, Lisa A; Barr, Joseph T

2018-01-01

New designs of eye wash stations have been developed in which the direction of water flow from the fountain has been reversed, with two water streams originating nasally in both eyes and flowing toward the temporal side of each eye. No study has been done to determine the ideal direction of water flow coming from the eye wash in relation to the eye. Ophthalmic eye examinations were conducted before and after the use of two eye wash stations with opposite water flow directionality. Fluorescein was instilled in both eyes before using an eye wash to measure the effectiveness of the water flow. Subjects were surveyed upon their experiences using the eye washes. Ophthalmic examination found no significant difference in the efficacy of the eye washes with nasal-to-temporal water flow when compared to temporal-to-nasal water flow direction.

A high-resolution land model coupled with groundwater lateral flow, human water regulation and the changes in soil freeze-thaw fronts

NASA Astrophysics Data System (ADS)

Xie, Z.; Zeng, Y.; Liu, S.; Gao, J.; Jia, B.; Qin, P.

2017-12-01

Both anthropogenic water regulation and groundwater lateral flow essentially affect groundwater table patterns. Their relationship is close because lateral flow recharges the groundwater depletion cone, which is induced by over-exploitation. And the movement of frost and thaw fronts (FTFs) affects soil water and thermal characteristics, as well as energy and water exchanges between land surface and the atmosphere. In this study, schemes describing groundwater lateral flow, human water regulation and the changes in soil freeze-thaw fronts were developed and incorporated into the Community Land Model 4.5. Then the model was applied in Heihe River Basin(HRB), an arid and semiarid region, northwest China. High resolution ( 1 km) numerical simulations showed that groundwater lateral flow driven by changes in water heads can essentially change the groundwater table pattern with the deeper water table appearing in the hillslope regions and shallower water table appearing in valley bottom regions and plains. Over the last decade, anthropogenic groundwater exploitation deepened the water table by approximately 2 m in the middle reaches of the HRB and rapidly reduced the terrestrial water storage, while irrigation increased soil moisture by approximately 0.1 m3 m-3. The water stored in the mainstream of the Heihe River was also reduced by human surface water withdrawal. The latent heat flux was increased by 30 W m-2 over the irrigated region, with an identical decrease in sensible heat flux. The simulated groundwater lateral flow was shown to effectively recharge the groundwater depletion cone caused by over-exploitation. The offset rate is higher in plains than mountainous regions. In addition, the simulated FTFs depth compared well with the observed data both in D66 station (permafrost) and Hulugou station (seasonally frozen ground). Over the HRB, the upstream area is permafrost region with maximum thawed depth at 2.5 m and lower region is seasonal frozen ground region

A flow-simulation model of the tidal Potomac River

USGS Publications Warehouse

Schaffranek, Raymond W.

1987-01-01

A one-dimensional model capable of simulating flow in a network of interconnected channels has been applied to the tidal Potomac River including its major tributaries and embayments between Washington, D.C., and Indian Head, Md. The model can be used to compute water-surface elevations and flow discharges at any of 66 predetermined locations or at any alternative river cross sections definable within the network of channels. In addition, the model can be used to provide tidal-interchange flow volumes and to evaluate tidal excursions and the flushing properties of the riverine system. Comparisons of model-computed results with measured watersurface elevations and discharges demonstrate the validity and accuracy of the model. Tidal-cycle flow volumes computed by the calibrated model have been verified to be within an accuracy of ? 10 percent. Quantitative characteristics of the hydrodynamics of the tidal river are identified and discussed. The comprehensive flow data provided by the model can be used to better understand the geochemical, biological, and other processes affecting the river's water quality.

Imaging water velocity and volume fraction distributions in water continuous multiphase flows using inductive flow tomography and electrical resistance tomography

NASA Astrophysics Data System (ADS)

Meng, Yiqing; Lucas, Gary P.

2017-05-01

This paper presents the design and implementation of an inductive flow tomography (IFT) system, employing a multi-electrode electromagnetic flow meter (EMFM) and novel reconstruction techniques, for measuring the local water velocity distribution in water continuous single and multiphase flows. A series of experiments were carried out in vertical-upward and upward-inclined single phase water flows and ‘water continuous’ gas-water and oil-gas-water flows in which the velocity profiles ranged from axisymmetric (single phase and vertical-upward multiphase flows) to highly asymmetric (upward-inclined multiphase flows). Using potential difference measurements obtained from the electrode array of the EMFM, local axial velocity distributions of the continuous water phase were reconstructed using two different IFT reconstruction algorithms denoted RT#1, which assumes that the overall water velocity profile comprises the sum of a series of polynomial velocity components, and RT#2, which is similar to RT#1 but which assumes that the zero’th order velocity component may be replaced by an axisymmetric ‘power law’ velocity distribution. During each experiment, measurement of the local water volume fraction distribution was also made using the well-established technique of electrical resistance tomography (ERT). By integrating the product of the local axial water velocity and the local water volume fraction in the cross section an estimate of the water volumetric flow rate was made which was compared with a reference measurement of the water volumetric flow rate. In vertical upward flows RT#2 was found to give rise to water velocity profiles which are consistent with the previous literature although the profiles obtained in the multiphase flows had relatively higher central velocity peaks than was observed for the single phase profiles. This observation was almost certainly a result of the transfer of axial momentum from the less dense dispersed phases to the water

CALCULATION OF NONLINEAR CONFIDENCE AND PREDICTION INTERVALS FOR GROUND-WATER FLOW MODELS.

USGS Publications Warehouse

Cooley, Richard L.; Vecchia, Aldo V.

1987-01-01

A method is derived to efficiently compute nonlinear confidence and prediction intervals on any function of parameters derived as output from a mathematical model of a physical system. The method is applied to the problem of obtaining confidence and prediction intervals for manually-calibrated ground-water flow models. To obtain confidence and prediction intervals resulting from uncertainties in parameters, the calibrated model and information on extreme ranges and ordering of the model parameters within one or more independent groups are required. If random errors in the dependent variable are present in addition to uncertainties in parameters, then calculation of prediction intervals also requires information on the extreme range of error expected. A simple Monte Carlo method is used to compute the quantiles necessary to establish probability levels for the confidence and prediction intervals. Application of the method to a hypothetical example showed that inclusion of random errors in the dependent variable in addition to uncertainties in parameters can considerably widen the prediction intervals.

A water tunnel flow visualization study of the F-15

NASA Technical Reports Server (NTRS)

Lorincz, D. J.

1978-01-01

Water tunnel studies were performed to qualitatively define the flow field of the F-15 aircraft. Two lengthened forebodies, one with a modified cross-sectional shape, were tested in addition to the basic forebody. Particular emphasis was placed on defining vortex flows generated at high angles of attack. The flow visualization tests were conducted in the Northrop diagnostic water tunnel using a 1/48-scale model of the F-15. Flow visualization pictures were obtained over an angle-of-attack range to 55 deg and sideslip angles up to 10 deg. The basic aircraft configuration was investigated in detail to determine the vortex flow field development, vortex path, and vortex breakdown characteristics as a function of angle of attack and sideslip. Additional tests showed that the wing upper surface vortex flow fields were sensitive to variations in inlet mass flow ratio and inlet cowl deflection angle. Asymmetries in the vortex systems generated by each of the three forebodies were observed in the water tunnel at zero sideslip and high angles of attack.

Numerical simulation of water and sand blowouts when penetrating through shallow water flow formations in deep water drilling

NASA Astrophysics Data System (ADS)

Ren, Shaoran; Liu, Yanmin; Gong, Zhiwu; Yuan, Yujie; Yu, Lu; Wang, Yanyong; Xu, Yan; Deng, Junyu

2018-02-01

In this study, we applied a two-phase flow model to simulate water and sand blowout processes when penetrating shallow water flow (SWF) formations during deepwater drilling. We define `sand' as a pseudo-component with high density and viscosity, which can begin to flow with water when a critical pressure difference is attained. We calculated the water and sand blowout rates and analyzed the influencing factors from them, including overpressure of the SWF formation, as well as its zone size, porosity and permeability, and drilling speed (penetration rate). The obtained data can be used for the quantitative assessment of the potential severity of SWF hazards. The results indicate that overpressure of the SWF formation and its zone size have significant effects on SWF blowout. A 10% increase in the SWF formation overpressure can result in a more than 90% increase in the cumulative water blowout and a 150% increase in the sand blowout when a typical SWF sediment is drilled. Along with the conventional methods of well flow and pressure control, chemical plugging, and the application of multi-layer casing, water and sand blowouts can be effectively reduced by increasing the penetration rate. As such, increasing the penetration rate can be a useful measure for controlling SWF hazards during deepwater drilling.

Modifications of the U.S. Geological Survey modular, finite-difference, ground-water flow model to read and write geographic information system files

USGS Publications Warehouse

Orzol, Leonard L.; McGrath, Timothy S.

1992-01-01

This report documents modifications to the U.S. Geological Survey modular, three-dimensional, finite-difference, ground-water flow model, commonly called MODFLOW, so that it can read and write files used by a geographic information system (GIS). The modified model program is called MODFLOWARC. Simulation programs such as MODFLOW generally require large amounts of input data and produce large amounts of output data. Viewing data graphically, generating head contours, and creating or editing model data arrays such as hydraulic conductivity are examples of tasks that currently are performed either by the use of independent software packages or by tedious manual editing, manipulating, and transferring data. Programs such as GIS programs are commonly used to facilitate preparation of the model input data and analyze model output data; however, auxiliary programs are frequently required to translate data between programs. Data translations are required when different programs use different data formats. Thus, the user might use GIS techniques to create model input data, run a translation program to convert input data into a format compatible with the ground-water flow model, run the model, run a translation program to convert the model output into the correct format for GIS, and use GIS to display and analyze this output. MODFLOWARC, avoids the two translation steps and transfers data directly to and from the ground-water-flow model. This report documents the design and use of MODFLOWARC and includes instructions for data input/output of the Basic, Block-centered flow, River, Recharge, Well, Drain, Evapotranspiration, General-head boundary, and Streamflow-routing packages. The modification to MODFLOW and the Streamflow-Routing package was minimized. Flow charts and computer-program code describe the modifications to the original computer codes for each of these packages. Appendix A contains a discussion on the operation of MODFLOWARC using a sample problem.

Geohydrology and simulated ground-water flow in northwestern Elkhart County, Indiana

USGS Publications Warehouse

Arihood, L.D.; Cohen, D.A.

1998-01-01

In 1994, the U.S. Geological Survey, in cooperation with the U.S. Environmental Protection Agency and the City of Elkhart, developed a ground-water model of the Elkhart, Indiana, area to determine the avail-ability and source of water at potential new well fields. The modeled area covered 190 square miles of northwestern Elkhart County and a small part of southern Michigan. Three Superfund sites and several other sites in this area are undergoing environmental cleanup. The model would be used to guide the location of well fields so that Superfund sites and environmental cleanup areas would not be within recharge areas for the well fields. The City of Elkhart obtains its water supply from two aquifers separated by a generally continuous confining unit. The upper aquifer is composed primarily of sand and gravel of glacial origin. Thickness of the upper aquifer ranges from 0 to 116 feet and averages 47 feet. The lower aquifer is composed of sand and gravel with interbedded lenses of silt and clay. Thickness of the lower aquifer ranges from 1 to 335 feet and averages 35 feet. The intervening confining unit is composed of silt and clay with interbedded sand and gravel; the confining unit ranges from 0 to 177 feet, with an average thickness of 27 feet. Flow through the aquifers is generally horizontal vertically downward from the upper aquifer, through the confining unit, and into the lower aquifer, except where flow is vertically upward at the St. Joseph River and other large streams. The hydraulic characteristics of the aquifers and confining unit were estimated by analyzing aquifer-test data from well drillers? logs and by calibration of the model. The horizontal hydraulic conductivity of the upper aquifer is 170 feet per day within about 1 mile of the St. Joseph and Elkhart Rivers and 370 feet per day at distances greater than about 1 mile. The horizontal hydraulic conductivity of the lower aquifer is 370 feet per day throughout the modeled area, with the exception

Fundamentals of Ground-Water Modeling

EPA Pesticide Factsheets

This paper presents an overview of the essential components of ground-water flow and contaminant transport modeling in saturated porous media. While fractured rocks and fractured porous rocks may behave like porous media with respect to many flow and...

Hydro-dynamic damping theory in flowing water

NASA Astrophysics Data System (ADS)

Monette, C.; Nennemann, B.; Seeley, C.; Coutu, A.; Marmont, H.

2014-03-01

Fluid-structure interaction (FSI) has a major impact on the dynamic response of the structural components of hydroelectric turbines. On mid-head to high-head Francis runners, the rotor-stator interaction (RSI) phenomenon always has to be considered carefully during the design phase to avoid operational issues later on. The RSI dynamic response amplitudes are driven by three main factors: (1) pressure forcing amplitudes, (2) excitation frequencies in relation to natural frequencies and (3) damping. The prediction of the two first factors has been largely documented in the literature. However, the prediction of fluid damping has received less attention in spite of being critical when the runner is close to resonance. Experimental damping measurements in flowing water on hydrofoils were presented previously. Those results showed that the hydro-dynamic damping increased linearly with the flow. This paper presents development and validation of a mathematical model, based on momentum exchange, to predict damping due to fluid structure interaction in flowing water. The model is implemented as an analytical procedure for simple structures, such as cantilever beams, but is also implemented in more general ways using three different approaches for more complex structures such as runner blades: a finite element procedure, a CFD modal work based approach and a CFD 1DOF approach. The mathematical model and all three implementation approaches are shown to agree well with experimental results.

Modeling Water Waves with Smoothed Particle Hydrodynamics

DTIC Science & Technology

2013-09-30

SPH Model for Water Waves and Other Free Surface Flows ...Lagrangian nature of SPH allows the modeling of wave breaking, surf zones, ship waves, and wave-structure interaction, where the free surface becomes...proving to be a competent modeling scheme for free surface flows in three dimensions including the complex flows of the surf zone. As the GPU

Experimental Study on the Flow Regimes and Pressure Gradients of Air-Oil-Water Three-Phase Flow in Horizontal Pipes

PubMed Central

Al-Hadhrami, Luai M.; Shaahid, S. M.; Tunde, Lukman O.; Al-Sarkhi, A.

2014-01-01

An experimental investigation has been carried out to study the flow regimes and pressure gradients of air-oil-water three-phase flows in 2.25 ID horizontal pipe at different flow conditions. The effects of water cuts, liquid and gas velocities on flow patterns and pressure gradients have been studied. The experiments have been conducted at 20°C using low viscosity Safrasol D80 oil, tap water and air. Superficial water and oil velocities were varied from 0.3 m/s to 3 m/s and air velocity varied from 0.29 m/s to 52.5 m/s to cover wide range of flow patterns. The experiments were performed for 10% to 90% water cuts. The flow patterns were observed and recorded using high speed video camera while the pressure drops were measured using pressure transducers and U-tube manometers. The flow patterns show strong dependence on water fraction, gas velocities, and liquid velocities. The observed flow patterns are stratified (smooth and wavy), elongated bubble, slug, dispersed bubble, and annular flow patterns. The pressure gradients have been found to increase with the increase in gas flow rates. Also, for a given superficial gas velocity, the pressure gradients increased with the increase in the superficial liquid velocity. The pressure gradient first increases and then decreases with increasing water cut. In general, phase inversion was observed with increase in the water cut. The experimental results have been compared with the existing unified Model and a good agreement has been noticed. PMID:24523645

Experimental study on the flow regimes and pressure gradients of air-oil-water three-phase flow in horizontal pipes.

PubMed

Al-Hadhrami, Luai M; Shaahid, S M; Tunde, Lukman O; Al-Sarkhi, A

2014-01-01

An experimental investigation has been carried out to study the flow regimes and pressure gradients of air-oil-water three-phase flows in 2.25 ID horizontal pipe at different flow conditions. The effects of water cuts, liquid and gas velocities on flow patterns and pressure gradients have been studied. The experiments have been conducted at 20 °C using low viscosity Safrasol D80 oil, tap water and air. Superficial water and oil velocities were varied from 0.3 m/s to 3 m/s and air velocity varied from 0.29 m/s to 52.5 m/s to cover wide range of flow patterns. The experiments were performed for 10% to 90% water cuts. The flow patterns were observed and recorded using high speed video camera while the pressure drops were measured using pressure transducers and U-tube manometers. The flow patterns show strong dependence on water fraction, gas velocities, and liquid velocities. The observed flow patterns are stratified (smooth and wavy), elongated bubble, slug, dispersed bubble, and annular flow patterns. The pressure gradients have been found to increase with the increase in gas flow rates. Also, for a given superficial gas velocity, the pressure gradients increased with the increase in the superficial liquid velocity. The pressure gradient first increases and then decreases with increasing water cut. In general, phase inversion was observed with increase in the water cut. The experimental results have been compared with the existing unified Model and a good agreement has been noticed.

Groundwater Flow Model for Taos, New Mexico

NASA Astrophysics Data System (ADS)

Burck, P. W.; Barroll, P. W.; Core, A. B.; Rappuhn, D.

2003-12-01

The New Mexico Office of the State Engineer - Hydrology Bureau (OSE) has developed a regional groundwater flow model for Taos, New Mexico. The MODFLOW 2000 model will serve as a tool to evaluate alternatives in settlement negotiations in an on-going water rights adjudication. If current settlement negotiations fail, it is conceivable that the model might be used in support of litigation. OSE produced the model in cooperation with technical representatives of the various parties to the adjudication. Regional hydrogeologic data including well records, aquifer test results, stream flow measurements and seepage studies have been shared relatively freely among the parties. A recent deep drilling program conducted in conjunction with the negotiation effort has added substantially to the hydrogeologic data set. Among the hydrologic processes simulated by the model are mountain front recharge; areal recharge from precipitation; evapotranspiration; discharge from springs; river and stream flow; accretions to groundwater from irrigation return flow, seepage from acequias, canals, and ditches, and deep percolation; and pumping by municipal entities and mutual domestic water users associations. The resulting model files are available for all parties to review and evaluate. Comments are assessed and many have resulted in significant improvements to the model. At this stage, however, it is unclear whether adopting this cooperative approach will increase the likelihood of model acceptance by the parties.

COMPILATION OF GROUND WATER MODELS

EPA Science Inventory

The full report presents an overview of currently available computer-based simulation models for ground-water flow, solute and heat transport, and hydrogeochemistry in both porous media and fractured rock. Separate sections address multiphase flow and related chemical species tra...

Trees influence preferencial flow and water uptake in tropical savanna

NASA Astrophysics Data System (ADS)

Benegas, Laura; Bargues-Tobella, Aida; Hasselquist, Niles; Malmer, Anders; Ilstedt, Ulrik

2017-04-01

To address potential competition between trees and grasses for soil water, and to disentangle the main process responsible for local soil water dynamics in pasture ecosystems, we conducted a study of the soil water content and water source partitioning of grasses and trees within a pasture in the Copan River catchment, Honduras. We used differences in the 2H/1H (δD) isotopic signature of soil water (δSW) and the local meteoric water line (LMWL; δLMWL) as a relative index of evaporation, following a recent model proposed by Hasselquist et al (under review). The model uses Lc-excess calculated as the absolute value of the difference between measured δD and that predicted by the local meteoric water line (lc-excess = ¦δDM - δDP¦). Lc-excess values close to zero indicate little difference between soil water samples and local precipitation, whereas larger values indicate a greater degree of evaporation .()...(adapted from Landwehr and Coplen, 2006). From the relation between Lc-excess and SWC, we can tease apart different processes by which trees influence local soil water dynamics, where one such processes indicate that if preferential flow, i.e quick flows through macropores that by-pass the soil matrix, is the main pathway for water movement in the soil, then the Lc-excess values of soil water at deeper depths will be closer to zero than those of the surface soil, whereas relatively higher Lc-excess values would indicate increasing dominance of matrix flow. We found that soil underneath trees was wetter than underneath grasses at the dry season and we can relate this with a lack of clear relationship between Lc-excess and SWC and with the treés apparent shift to groundwater sources for root uptake especially in the dry season. Due to the positive correlation between Lc-excess and SWC under trees and due to the lower Lc-excess values found at subsoil below trees during the dry season, we can infer that preferential flow is also facilitated by the trees

49 CFR 229.111 - Water-flow indicator.

Code of Federal Regulations, 2010 CFR

2010-10-01

... 49 Transportation 4 2010-10-01 2010-10-01 false Water-flow indicator. 229.111 Section 229.111....111 Water-flow indicator. (a) Steam generators shall be equipped with an illuminated visual return water-flow indicator. (b) Steam generators shall be equipped with an operable test valve or other means...

49 CFR 229.111 - Water-flow indicator.

Code of Federal Regulations, 2011 CFR

2011-10-01

... 49 Transportation 4 2011-10-01 2011-10-01 false Water-flow indicator. 229.111 Section 229.111....111 Water-flow indicator. (a) Steam generators shall be equipped with an illuminated visual return water-flow indicator. (b) Steam generators shall be equipped with an operable test valve or other means...

49 CFR 229.111 - Water-flow indicator.

Code of Federal Regulations, 2013 CFR

2013-10-01

... 49 Transportation 4 2013-10-01 2013-10-01 false Water-flow indicator. 229.111 Section 229.111....111 Water-flow indicator. (a) Steam generators shall be equipped with an illuminated visual return water-flow indicator. (b) Steam generators shall be equipped with an operable test valve or other means...

49 CFR 229.111 - Water-flow indicator.

Code of Federal Regulations, 2014 CFR

2014-10-01

... 49 Transportation 4 2014-10-01 2014-10-01 false Water-flow indicator. 229.111 Section 229.111....111 Water-flow indicator. (a) Steam generators shall be equipped with an illuminated visual return water-flow indicator. (b) Steam generators shall be equipped with an operable test valve or other means...

49 CFR 229.111 - Water-flow indicator.

Code of Federal Regulations, 2012 CFR

2012-10-01

... 49 Transportation 4 2012-10-01 2012-10-01 false Water-flow indicator. 229.111 Section 229.111....111 Water-flow indicator. (a) Steam generators shall be equipped with an illuminated visual return water-flow indicator. (b) Steam generators shall be equipped with an operable test valve or other means...

A conceptual model of ground-water flow in the eastern Snake River Plain aquifer at the Idaho National Laboratory and vicinity with implications for contaminant transport

USGS Publications Warehouse

Ackerman, Daniel J.; Rattray, Gordon W.; Rousseau, Joseph P.; Davis, Linda C.; Orr, Brennon R.

2006-01-01

Ground-water flow in the west-central part of the eastern Snake River Plain aquifer is described in a conceptual model that will be used in numerical simulations to evaluate contaminant transport at the Idaho National Laboratory (INL) and vicinity. The model encompasses an area of 1,940 square miles (mi2) and includes most of the 890 mi2 of the INL. A 50-year history of waste disposal associated with research activities at the INL has resulted in measurable concentrations of waste contaminants in the aquifer. A thorough understanding of the fate and movement of these contaminants in the subsurface is needed by the U.S. Department of Energy to minimize the effect that contaminated ground water may have on the region and to plan effectively for remediation. Three hydrogeologic units were used to represent the complex stratigraphy of the aquifer in the model area. Collectively, these hydrogeologic units include at least 65 basalt-flow groups, 5 andesite-flow groups, and 61 sedimentary interbeds. Three rhyolite domes in the model area extend deep enough to penetrate the aquifer. The rhyolite domes are represented in the conceptual model as low permeability, vertical pluglike masses, and are not included as part of the three primary hydrogeologic units. Broad differences in lithology and large variations in hydraulic properties allowed the heterogeneous, anisotropic basalt-flow groups, andesite-flow groups, and sedimentary interbeds to be grouped into three hydrogeologic units that are conceptually homogeneous and anisotropic. Younger rocks, primarily thin, densely fractured basalt, compose hydrogeologic unit 1; younger rocks, primarily of massive, less densely fractured basalt, compose hydrogeologic unit 2; and intermediate-age rocks, primarily of slightly-to-moderately altered, fractured basalt, compose hydrogeologic unit 3. Differences in hydraulic properties among adjacent hydrogeologic units result in much of the large-scale heterogeneity and anisotropy of the

Variable thickness transient ground-water flow model. Volume 3. Program listings

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

Reisenauer, A.E.

1979-12-01

The Assessment of Effectiveness of Geologic Isolation Systems (AEGIS) Program is developing and applying the methodology for assessing the far-field, long-term post-closure safety of deep geologic nuclear waste repositories. AEGIS is being performed by Pacific Northwest Laboratory (PNL) under contract with the Office of Nuclear Waste Isolation (OWNI) for the Department of Energy (DOE). One task within AEGIS is the development of methodology for analysis of the consequences (water pathway) from loss of repository containment as defined by various release scenarios. Analysis of the long-term, far-field consequences of release scenarios requires the application of numerical codes which simulate the hydrologicmore » systems, model the transport of released radionuclides through the hydrologic systems to the biosphere, and, where applicable, assess the radiological dose to humans. Hydrologic and transport models are available at several levels of complexity or sophistication. Model selection and use are determined by the quantity and quality of input data. Model development under AEGIS and related programs provides three levels of hydrologic models, two levels of transport models, and one level of dose models (with several separate models). This is the third of 3 volumes of the description of the VTT (Variable Thickness Transient) Groundwater Hydrologic Model - second level (intermediate complexity) two-dimensional saturated groundwater flow.« less

A New Method of Studying the Flow of the Water Along the Bottom of a Model of a Flying-boat Hull

NASA Technical Reports Server (NTRS)

Ward, Kenneth E

1940-01-01

A new method of studying the flow of the water along the bottom of a model of a flying-boat hull is described. In this method, the model is fitted with a transparent bottom and is divided down the center line by a bulkhead. The flow is observed and photographed through one-half of the model by means of the diffused illumination from a battery of lamps contained in the other half of the model. Photographs of the flow, particularly of the changes that occur when the step ventilates, are shown. The results of the present investigation indicate that the method has considerable promise, chiefly in connection with motion-picture studies.

Analog Model Study of Ground-Water Flow in the Rehoboth Bay Area, Delaware.

DTIC Science & Technology

The study concerns ground- water flow in the Rehoboth Bay area, Delaware, a coastal area which depends on ground water for its fresh- water supply...Increased pumping demands may threaten to lower the water table and allow salt- water intrusion into the wells. The study was conducted using a viscous...use two different glycerin solutions to make observations and predict interactions between fresh and salt water in nature. Results indicate the

Kauai Groundwater Flow Model

DOE Data Explorer

Nicole Lautze

2015-01-01

Groundwater flow model for Kauai. Data is from the following sources: Whittier, R. and A.I. El-Kadi. 2014. Human and Environmental Risk Ranking of Onsite Sewage Disposal Systems For the Hawaiian Islands of Kauai, Molokai, Maui, and Hawaii – Final. Prepared by the University of Hawaii, Dept. of Geology and Geophysics for the State of Hawaii Dept. of Health, Safe Drinking Water Branch. September 2014.; and Whittier, R.B., K. Rotzoll, S. Dhal, A.I. El-Kadi, C. Ray, G. Chen, and D. Chang. 2004. Hawaii Source Water Assessment Program Report – Volume IV – Island of Kauai Source Water Assessment Program Report. Prepared for the Hawaii Department of Health, Safe Drinking Water Branch. University of Hawaii, Water Resources Research Center. Updated 2015.

Simulation of ground-water flow in the Intermediate and Floridan aquifer systems in Peninsular Florida

USGS Publications Warehouse

Sepúlveda, Nicasio

2002-01-01

A numerical model of the intermediate and Floridan aquifer systems in peninsular Florida was used to (1) test and refine the conceptual understanding of the regional ground-water flow system; (2) develop a data base to support subregional ground-water flow modeling; and (3) evaluate effects of projected 2020 ground-water withdrawals on ground-water levels. The four-layer model was based on the computer code MODFLOW-96, developed by the U.S. Geological Survey. The top layer consists of specified-head cells simulating the surficial aquifer system as a source-sink layer. The second layer simulates the intermediate aquifer system in southwest Florida and the intermediate confining unit where it is present. The third and fourth layers simulate the Upper and Lower Floridan aquifers, respectively. Steady-state ground-water flow conditions were approximated for time-averaged hydrologic conditions from August 1993 through July 1994 (1993-94). This period was selected based on data from Upper Floridan a quifer wells equipped with continuous water-level recorders. The grid used for the ground-water flow model was uniform and composed of square 5,000-foot cells, with 210 columns and 300 rows.

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

Shale Gas Well, Hydraulic Fracturing, and Formation Data to Support Modeling of Gas and Water Flow in Shale Formations

NASA Astrophysics Data System (ADS)

Edwards, Ryan W. J.; Celia, Michael A.

2018-04-01

The potential for shale gas development and hydraulic fracturing to cause subsurface water contamination has prompted a number of modeling studies to assess the risk. A significant impediment for conducting robust modeling is the lack of comprehensive publicly available information and data about the properties of shale formations, shale wells, the process of hydraulic fracturing, and properties of the hydraulic fractures. We have collated a substantial amount of these data that are relevant for modeling multiphase flow of water and gas in shale gas formations. We summarize these data and their sources in tabulated form.

Environmental water incentive policy and return flows

NASA Astrophysics Data System (ADS)

Qureshi, M. E.; Schwabe, K.; Connor, J.; Kirby, M.

2010-04-01

With increasing urban, industrial, and agricultural water demand and projected reduced supply under climate change, allocations to the environment are critically low in many arid and semiarid basins. Consequently, many governments are striving to augment environmental flows, often through market-oriented mechanisms that involve compensating irrigated agriculture, the largest water user in most basins, for reducing diversions. A widely documented challenge with policies to recover water for the environment arises because part of the water diversion reduction can form the basis for downstream consumptive water rights or environmental flows. This article gives an empirical comparison of two incentive policies to acquire water for environmental flows for a part of the Murray-Darling Basin (MDB), Australia. One policy consists of paying irrigators and water delivery firms to make capital and management investments that improve on-farm irrigation and water-conveyance; the other policy consists of having the government buy water from irrigators on the active MDB water market. The results show that the first option results in relatively larger return flow reduction, while the second option tends to induce significant irrigated land retirement with relatively large reductions in consumptive use and small reductions in return flow. In cases where irrigation losses result in little useful return flow (e.g., evaporative loss reduction or during drought in some instances), efficiency-improving investments may provide some cost-effective opportunities. Where a large portion of loss forms valuable return flow, it is difficult to make a case for the cost-effectiveness of policies involving payments for investments in irrigation and conveyance system upgrades.

Investigation and mitigation of condensation induced water hammer by stratified flow experiments

NASA Astrophysics Data System (ADS)

Kadakia, Hiral J.

This research primarily focuses on the possibility of using stratified flow in preventing an occurrence of condensation induced water hammer (CIWH) in horizontal pipe involving steam and subcooled water. A two-phase flow loop simulating the passive safety systems of an advanced light water reactor was constructed and a series of stratified flow experiments were carried out involving a system of subcooled water, saturated water, and steam. Special instruments were designed to measure steam flow rate and subcooled liquid velocity. These experiments showed that when flow field conditions meet certain criteria CIWH does occur. Flow conditions used in experiments were typically observed in passive safety systems of an advanced light water cooled reactor. This research summarizes a) literature research and other experimental data that signify an occurrence of CIWH, b) experiments in an effort to show an occurrence of CIWH and the ability to prevent CIWH, c) qualitative and quantitative results to underline the mechanism of CIWH, d) experiments that show CIWH can be prevented under certain conditions, and e) guidelines for the safe operating conditions. Based on initial experiment results it was observed that Bernoulli's effect can play an important role in wave formation and instability. A separate effect table top experiment was constructed with plexi-glass. A series of entrance effect tests and stratified experiments were carried out with different fluids to study wave formation and wave bridging. Special test series experiments were carried out to investigate the presence of a saturated layer. The effect of subcooled water and steam flow on wedge length and depth were recorded. These experiments helped create a model which calculates wedge and depth of wedge for a given condition of steam and subcooled water. A very good comparison between the experiment results and the model was obtained. These experiments also showed that the presence of saturated layer can mitigate

Use of a ground-water flow model to delineate contributing areas to the Puchack Well Field, Pennsauken township and vicinity, Camden county, New Jersey

USGS Publications Warehouse

Pope, Daryll A.; Watt, Martha K.

2005-01-01

The New Jersey Department of Environmental Protection (NJDEP) Well Head Protection Program, developed in response to the 1986 Federal Safe Drinking Water Act Amendments, requires delineation of Well Head Protection Areas (WHPA's), commonly called contributing areas, for all public and non-community water-supply wells in New Jersey. Typically, WHPA's for public community water-supply wells in New Jersey are delineated using a two-dimensional ground-water flow model incorporating the regional hydraulic gradient; however, NJDEP guidelines allow for the use of a three-dimensional flow model to delineate contributing areas to wells in complex hydrogeologic settings. The Puchack well field in Pennsauken Township, Camden County, N.J., is an area of strong hydraulic connection between the Lower aquifer of the Potomac-Raritan-Magothy aquifer system and the Delaware River. Interactions among and within the public-supply well fields in the area are complex. To delineate the contributing area to the Puchack well field, the U.S. Geological Survey, in cooperation with the NJDEP, developed an 11-layer ground-water flow model of the Potomac-Raritan-Magothy aquifer system in the Pennsauken Township area to simulate flow in the vicinity of the well field. The model incorporates the interaction between the aquifer system and the Delaware River, and includes boundary flows from an existing regional model of the Camden area. Recharge used in the model ranged from 4.5 to 14 inches per year, and horizontal hydraulic conductivity ranged from 50 to 250 feet per day. Values of vertical hydraulic conductivity ranging from 0.001 to 0.5 feet per day were assigned to zones created on the basis of variations in hydrogeologic conditions observed in geophysical logs from wells. A steady-state simulation was used to calibrate the model to synoptic water-level data collected in March 1998. Near the Puchack well field, simulated heads generally were within 1 foot of the measured heads in both the

Continental-scale river flow in climate models

NASA Technical Reports Server (NTRS)

Miller, James R.; Russell, Gary L.; Caliri, Guilherme

1994-01-01

The hydrologic cycle is a major part of the global climate system. There is an atmospheric flux of water from the ocean surface to the continents. The cycle is closed by return flow in rivers. In this paper a river routing model is developed to use with grid box climate models for the whole earth. The routing model needs an algorithm for the river mass flow and a river direction file, which has been compiled for 4 deg x 5 deg and 2 deg x 2.5 deg resolutions. River basins are defined by the direction files. The river flow leaving each grid box depends on river and lake mass, downstream distance, and an effective flow speed that depends on topography. As input the routing model uses monthly land source runoff from a 5-yr simulation of the NASA/GISS atmospheric climate model (Hansen et al.). The land source runoff from the 4 deg x 5 deg resolution model is quartered onto a 2 deg x 2.5 deg grid, and the effect of grid resolution is examined. Monthly flow at the mouth of the world's major rivers is compared with observations, and a global error function for river flow is used to evaluate the routing model and its sensitivity to physical parameters. Three basinwide parameters are introduced: the river length weighted by source runoff, the turnover rate, and the basinwide speed. Although the values of these parameters depend on the resolution at which the rivers are defined, the values should converge as the grid resolution becomes finer. When the routing scheme described here is coupled with a climate model's source runoff, it provides the basis for closing the hydrologic cycle in coupled atmosphere-ocean models by realistically allowing water to return to the ocean at the correct location and with the proper magnitude and timing.

Model simulations of flood and debris flow timing in steep catchments after wildfire

NASA Astrophysics Data System (ADS)

Rengers, F. K.; McGuire, L. A.; Kean, J. W.; Staley, D. M.; Hobley, D. E. J.

2016-08-01

Debris flows are a typical hazard on steep slopes after wildfire, but unlike debris flows that mobilize from landslides, most postwildfire debris flows are generated from water runoff. The majority of existing debris flow modeling has focused on landslide-triggered debris flows. In this study we explore the potential for using process-based rainfall-runoff models to simulate the timing of water flow and runoff-generated debris flows in recently burned areas. Two different spatially distributed hydrologic models with differing levels of complexity were used: the full shallow water equations and the kinematic wave approximation. Model parameter values were calibrated in two different watersheds, spanning two orders of magnitude in drainage area. These watersheds were affected by the 2009 Station Fire in the San Gabriel Mountains, CA, USA. Input data for the numerical models were constrained by time series of soil moisture, flow stage, and rainfall collected at field sites, as well as high-resolution lidar-derived digital elevation models. The calibrated parameters were used to model a third watershed in the burn area, and the results show a good match with observed timing of flow peaks. The calibrated roughness parameter (Manning's n) was generally higher when using the kinematic wave approximation relative to the shallow water equations, and decreased with increasing spatial scale. The calibrated effective watershed hydraulic conductivity was low for both models, even for storms occurring several months after the fire, suggesting that wildfire-induced changes to soil-water infiltration were retained throughout that time. Overall, the two model simulations were quite similar suggesting that a kinematic wave model, which is simpler and more computationally efficient, is a suitable approach for predicting flood and debris flow timing in steep, burned watersheds.

Model simulations of flood and debris flow timing in steep catchments after wildfire

USGS Publications Warehouse

Rengers, Francis K.; McGuire, Luke; Kean, Jason W.; Staley, Dennis M.; Hobley, D.E.J

2016-01-01

Debris flows are a typical hazard on steep slopes after wildfire, but unlike debris flows that mobilize from landslides, most post-wildfire debris flows are generated from water runoff. The majority of existing debris-flow modeling has focused on landslide-triggered debris flows. In this study we explore the potential for using process-based rainfall-runoff models to simulate the timing of water flow and runoff-generated debris flows in recently burned areas. Two different spatially distributed hydrologic models with differing levels of complexity were used: the full shallow water equations and the kinematic wave approximation. Model parameter values were calibrated in two different watersheds, spanning two orders of magnitude in drainage area. These watersheds were affected by the 2009 Station Fire in the San Gabriel Mountains, CA, USA. Input data for the numerical models were constrained by time series of soil moisture, flow stage, and rainfall collected at field sites, as well as high-resolution lidar-derived digital elevation models. The calibrated parameters were used to model a third watershed in the burn area, and the results show a good match with observed timing of flow peaks. The calibrated roughness parameter (Manning's $n$) was generally higher when using the kinematic wave approximation relative to the shallow water equations, and decreased with increasing spatial scale. The calibrated effective watershed hydraulic conductivity was low for both models, even for storms occurring several months after the fire, suggesting that wildfire-induced changes to soil-water infiltration were retained throughout that time. Overall the two model simulations were quite similar suggesting that a kinematic wave model, which is simpler and more computationally efficient, is a suitable approach for predicting flood and debris flow timing in steep, burned watersheds.

Hydrogeology and simulation of ground-water flow near the Lantana Landfill, Palm Beach County, Florida

USGS Publications Warehouse

Russell, G.M.; Wexler, E.J.

1993-01-01

The Lantana landfill in Palm Beach County has a surface that is 40 to 50 feet above original ground level and consists of about 250 acres of compacted garbage and trash. Parts of the landfill are below the water table. Surface-resistivity measurements and water-quality analyses indicate that leachate-enriched ground water along the eastern perimeter of the landfill has moved about 500 feet eastward toward an adjacent lake. Concentrations of chloride and nutrients within the leachate-enriched ground water were greater than background concentrations. The surficial aquifer system in the area of the landfill consists primarily of sand of moderate permeability, from land surface to a depth of about 68 feet deep, and consists of sand interbedded with sandstone and limestone of high permeability from a depth of about 68 feet to a depth of 200 feet. The potentiometric surface in the landfill is higher than that in adjacent areas to the east, indicating ground-water movement from the landfill toward a lake to the east. Steady-state simulation of ground-water flow was made using a telescoping-grid technique where a model covering a large area is used to determine boundaries and fluxes for a finer scale model. A regional flow model encompassing a 500-square mile area in southeastern Palm Beach County was used to calculate ground-water fluxes in a 126.5-square mile subregional area. Boundary fluxes calculated by the subregional model were then used to calculate boundary fluxes for a local model of the 3.75-square mile area representing the Lantana landfill site and vicinity. Input data required for simulating ground-water flow in the study area were obtained from the regional flow models, thus, effectively coupling the models. Additional simulations were made using the local flow model to predict effects of possible remedial actions on the movement of solutes in the ground-water system. Possible remedial actions simulated included capping the landfill with an impermeable layer

Modeling two-phase flow in PEM fuel cell channels

NASA Astrophysics Data System (ADS)

Wang, Yun; Basu, Suman; Wang, Chao-Yang

2008-05-01

This paper is concerned with the simultaneous flow of liquid water and gaseous reactants in mini-channels of a proton exchange membrane (PEM) fuel cell. Envisaging the mini-channels as structured and ordered porous media, we develop a continuum model of two-phase channel flow based on two-phase Darcy's law and the M2 formalism, which allow estimate of the parameters key to fuel cell operation such as overall pressure drop and liquid saturation profiles along the axial flow direction. Analytical solutions of liquid water saturation and species concentrations along the channel are derived to explore the dependences of these physical variables vital to cell performance on operating parameters such as flow stoichiometric ratio and relative humility. The two-phase channel model is further implemented for three-dimensional numerical simulations of two-phase, multi-component transport in a single fuel-cell channel. Three issues critical to optimizing channel design and mitigating channel flooding in PEM fuel cells are fully discussed: liquid water buildup towards the fuel cell outlet, saturation spike in the vicinity of flow cross-sectional heterogeneity, and two-phase pressure drop. Both the two-phase model and analytical solutions presented in this paper may be applicable to more general two-phase flow phenomena through mini- and micro-channels.

Estimation of daily flow rate of photovoltaic water pumping systems using solar radiation data

NASA Astrophysics Data System (ADS)

Benghanem, M.; Daffallah, K. O.; Almohammedi, A.

2018-03-01

This paper presents a simple model which allows us to contribute in the studies of photovoltaic (PV) water pumping systems sizing. The nonlinear relation between water flow rate and solar power has been obtained experimentally in a first step and then used for performance prediction. The model proposed enables us to simulate the water flow rate using solar radiation data for different heads (50 m, 60 m, 70 m and 80 m) and for 8S × 3P PV array configuration. The experimental data are obtained with our pumping test facility located at Madinah site (Saudi Arabia). The performances are calculated using the measured solar radiation data of different locations in Saudi Arabia. Knowing the solar radiation data, we have estimated with a good precision the water flow rate Q in five locations (Al-Jouf, Solar Village, AL-Ahsa, Madinah and Gizan) in Saudi Arabia. The flow rate Q increases with the increase of pump power for different heads following the nonlinear model proposed.

Modifications to the modular three-dimensional finite-difference ground-water flow model used for the Columbia Plateau Regional Aquifer-System Analysis, Washington, Oregon, and Idaho

USGS Publications Warehouse

Hansen, A.J.

1993-01-01

The report documents modifications to the U.S. Geological Survey's modular three-dimensional finite-difference ground-water flow model used for a regional aquifer-system analysis of the Columbia Plateau. The report, which describes the concepts and mathematical basis for the modifications, is intended for potential users who are familiar with the original modular model. The modifications permit flow from a layer to any adjacent layer, allow the model to retain a cell of a layer that has been cut completely through by a canyon, and allow placing ground-water flow barriers on only specified branch conductances; a special version of the modified model uses a convergent grid. The report describes the data-input items that this modified model must read.

One-dimensional flow model of the river-hyporheic zone system

NASA Astrophysics Data System (ADS)

Pokrajac, D.

2016-12-01

The hyporheic zone is a shallow layer beneath natural streams that is characterized by intense exchange of water, nutrients, pollutants and thermal energy. Understanding these exchange processes is crucial for successful modelling of the river hydrodynamics and morphodynamics at various scales from the river corridor up to the river network scale (Cardenas, 2015). Existing simulation models of hyporheic exchange processes are either idealized models of the tracer movement through the river-hyporheic zone system (e.g. TSM, Bencala and Walters, 1983) or detailed models of turbulent flow in a stream, coupled with a conventional 2D Darcian groundwater model (e.g. Cardenas and Wilson, 2007). This paper presents an alternative approach which involves a simple 1-D simulation model of the hyporheic zone system based on the classical SWE equations coupled with the newly developed porous media analogue. This allows incorporating the effects of flow unsteadiness and non-Darcian parameterization od the drag term in the hyporheic zone model. The conceptual model of the stream-hyporheic zone system consists of a 1D model of the open channel flow in the river, coupled with a 1D model of the flow in the hyporheic zone via volume flux due to the difference in the water level in the river and the hyporheic zone. The interaction with the underlying groundwater aquifer is neglected, but coupling the present model with any conventional groundwater model is straightforward. The paper presents the derivation of the 1D flow equations for flow in the hyporheic zone, the details of the numerical scheme used for solving them and the model validation by comparison with published experimental data. References Bencala, K. E., and R. A. Walters (1983) "Simulation of solute transport in a mountain pool-and-riffle stream- a transient storage model", Water Resources Reseach 19(3): 718-724. Cardenas, M. B. (2015) "Hyporheic zone hydrologic science: A historical account of its emergence and a

The iFlow modelling framework v2.4: a modular idealized process-based model for flow and transport in estuaries

NASA Astrophysics Data System (ADS)

Dijkstra, Yoeri M.; Brouwer, Ronald L.; Schuttelaars, Henk M.; Schramkowski, George P.

2017-07-01

The iFlow modelling framework is a width-averaged model for the systematic analysis of the water motion and sediment transport processes in estuaries and tidal rivers. The distinctive solution method, a mathematical perturbation method, used in the model allows for identification of the effect of individual physical processes on the water motion and sediment transport and study of the sensitivity of these processes to model parameters. This distinction between processes provides a unique tool for interpreting and explaining hydrodynamic interactions and sediment trapping. iFlow also includes a large number of options to configure the model geometry and multiple choices of turbulence and salinity models. Additionally, the model contains auxiliary components, including one that facilitates easy and fast sensitivity studies. iFlow has a modular structure, which makes it easy to include, exclude or change individual model components, called modules. Depending on the required functionality for the application at hand, modules can be selected to construct anything from very simple quasi-linear models to rather complex models involving multiple non-linear interactions. This way, the model complexity can be adjusted to the application. Once the modules containing the required functionality are selected, the underlying model structure automatically ensures modules are called in the correct order. The model inserts iteration loops over groups of modules that are mutually dependent. iFlow also ensures a smooth coupling of modules using analytical and numerical solution methods. This way the model combines the speed and accuracy of analytical solutions with the versatility of numerical solution methods. In this paper we present the modular structure, solution method and two examples of the use of iFlow. In the examples we present two case studies, of the Yangtze and Scheldt rivers, demonstrating how iFlow facilitates the analysis of model results, the understanding of the

Environmental flows in hydro-economic models

NASA Astrophysics Data System (ADS)

Pereau, Jean-Christophe; Pryet, Alexandre

2018-03-01

The protection of environmental flows, as a management objective for a regulating agency, needs to be consistent with the aquifer water balance and the degree of resource renewability. A stylized hydro-economic model is used where natural recharge, which sustains environmental flows, is considered both in the aquifer water budget and in the welfare function as ecosystem damage. Groundwater recharge and the associated natural drainage may be neglected for aquifers containing fossil water, where the groundwater is mined. However, when dealing with an aquifer that constitutes a renewable resource, for which recharge is not negligible, natural drainage should explicitly appear in the water budget. In doing so, the optimum path of net extraction rate does not necessarily converge to the recharge rate, but depends on the costs associated with ecosystem damages. The optimal paths and equilibrium values for the water volume and water extraction are analytically derived, and numerical simulations based on the Western La Mancha aquifer (southwest Spain) illustrate the theoretical results of the study.

A controlled experiment in ground water flow model calibration

USGS Publications Warehouse

Hill, M.C.; Cooley, R.L.; Pollock, D.W.

1998-01-01

Nonlinear regression was introduced to ground water modeling in the 1970s, but has been used very little to calibrate numerical models of complicated ground water systems. Apparently, nonlinear regression is thought by many to be incapable of addressing such complex problems. With what we believe to be the most complicated synthetic test case used for such a study, this work investigates using nonlinear regression in ground water model calibration. Results of the study fall into two categories. First, the study demonstrates how systematic use of a well designed nonlinear regression method can indicate the importance of different types of data and can lead to successive improvement of models and their parameterizations. Our method differs from previous methods presented in the ground water literature in that (1) weighting is more closely related to expected data errors than is usually the case; (2) defined diagnostic statistics allow for more effective evaluation of the available data, the model, and their interaction; and (3) prior information is used more cautiously. Second, our results challenge some commonly held beliefs about model calibration. For the test case considered, we show that (1) field measured values of hydraulic conductivity are not as directly applicable to models as their use in some geostatistical methods imply; (2) a unique model does not necessarily need to be identified to obtain accurate predictions; and (3) in the absence of obvious model bias, model error was normally distributed. The complexity of the test case involved implies that the methods used and conclusions drawn are likely to be powerful in practice.Nonlinear regression was introduced to ground water modeling in the 1970s, but has been used very little to calibrate numerical models of complicated ground water systems. Apparently, nonlinear regression is thought by many to be incapable of addressing such complex problems. With what we believe to be the most complicated synthetic

Simulation of ground-water flow in the Saginaw Aquifer, Clinton, Eaton, and Ingham counties, Michigan

USGS Publications Warehouse

Holtschlag, David J.; Luukkonen, Carol L.; Nicholas, J.R.

1996-01-01

A numerical model was developed to simulate ground-water flow in the Tri-County region, which consists of Clinton, Eaton, and Ingham Counties, Michigan. This region includes a nine-township area surrounding Lansing, Michigan. The model simulates the regional response of the Saginaw aquifer to major groundwater withdrawals associated with public-supply wells. The Saginaw aquifer, which is in the Grand River and Saginaw Formations of Pennsylvanian age, is the primary source of ground water for Tri-County residents. The Saginaw aquifer is overlain by glacial deposits, which also are important ground-water sources in some locations. Flow in the Saginaw aquifer and the glacial deposits is simulated by discretizing the flow system into model cells arranged in two layers. Each cell, which corresponds to a land area of 0.0625 square mile, represents the locally averaged properties of the system. The spatial variation of hydraulic properties controlling ground-water flow was estimated by geostatistical analysis of 4,947 well logs. Parameter estimation, a form of nonlinear regression, was used to calibrate the flow model. Results of steady-state ground-water-flow simulations show close agreement between water flowing into and out of the model area for 1992 pumping conditions; standard error of the difference between simulated and measured heads is 14.7 feet. Simulation results for three alternative pumping scenarios for the year 2020 show that the glacial aquifer could be dewatered in places if hypothetical increases in pumping are not distributed throughout the Tri-County region. Contributing areas to public-supply wells in the nine-township area were delineated by a particle-tracking analysis. These areas cover about 121 square miles. Contributing areas for particles having travel times of 40 years or less cover about 42 square miles. Results of tritium sampling support results of model simulations to delineate contributing areas.

Simulation of flow and water quality of the Arroyo Colorado, Texas, 1989-99

USGS Publications Warehouse

Raines, Timothy H.; Miranda, Roger M.

2002-01-01

A model parameter set for use with the Hydrological Simulation Program—FORTRAN watershed model was developed to simulate flow and water quality for selected properties and constituents for the Arroyo Colorado from the city of Mission to the Laguna Madre, Texas. The model simulates flow, selected water-quality properties, and constituent concentrations. The model can be used to estimate a total maximum daily load for selected properties and constituents in the Arroyo Colorado. The model was calibrated and tested for flow with data measured during 1989–99 at three streamflow-gaging stations. The errors for total flow volume ranged from -0.1 to 29.0 percent, and the errors for total storm volume ranged from -15.6 to 8.4 percent. The model was calibrated and tested for water quality for seven properties and constituents with 1989–99 data. The model was calibrated sequentially for suspended sediment, water temperature, biochemical oxygen demand, dissolved oxygen, nitrate nitrogen, ammonia nitrogen, and orthophosphate. The simulated concentrations of the selected properties and constituents generally matched the measured concentrations available for the calibration and testing periods. The model was used to simulate total point- and nonpoint-source loads for selected properties and constituents for 1989–99 for urban, natural, and agricultural land-use types. About one-third to one-half of the biochemical oxygen demand and nutrient loads are from urban point and nonpoint sources, although only 13 percent of the total land use in the basin is urban.

Laboratory and Physical Modelling of Building Ventilation Flows

NASA Astrophysics Data System (ADS)

Hunt, Gary

2001-11-01

Heating and ventilating buildings accounts for a significant fraction of the total energy budget of cities and an immediate challenge in building physics is for the design of sustainable, low-energy buildings. Natural ventilation provides a low-energy solution as it harness the buoyancy force associated with temperature differences between the internal and external environment, and the wind to drive a ventilating flow. Modern naturally-ventilated buildings use innovative design solutions, e.g. glazed atria and solar chimneys, to enhance the ventilation and demand for these and other designs has far outstripped our understanding of the fluid mechanics within these buildings. Developing an understanding of the thermal stratification and movement of air provides a considerable challenge as the flows involve interactions between stratification and turbulence and often in complex geometries. An approach that has provided significant new insight into these flows and which has led to the development of design guidelines for architects and ventilation engineers is laboratory modelling at small-scale in water tanks combined with physical modelling. Density differences to drive the flow in simplified plexiglass models of rooms or buildings are provided by fresh and salt water solutions, and wind flow is represented by a mean flow in a flume tank. In tandom with the experiments, theoretical models that capture the essential physics of these flows have been developed in order to generalise the experimental results to a wide range of typical building geometries and operating conditions. This paper describes the application and outcomes of these modelling techniques to the study of a variety of natural ventilation flows in buildings.

Evaluating the sources of water to wells: Three techniques for metamodeling of a groundwater flow model

USGS Publications Warehouse

Fienen, Michael N.; Nolan, Bernard T.; Feinstein, Daniel T.

2016-01-01

For decision support, the insights and predictive power of numerical process models can be hampered by insufficient expertise and computational resources required to evaluate system response to new stresses. An alternative is to emulate the process model with a statistical “metamodel.” Built on a dataset of collocated numerical model input and output, a groundwater flow model was emulated using a Bayesian Network, an Artificial neural network, and a Gradient Boosted Regression Tree. The response of interest was surface water depletion expressed as the source of water-to-wells. The results have application for managing allocation of groundwater. Each technique was tuned using cross validation and further evaluated using a held-out dataset. A numerical MODFLOW-USG model of the Lake Michigan Basin, USA, was used for the evaluation. The performance and interpretability of each technique was compared pointing to advantages of each technique. The metamodel can extend to unmodeled areas.

Measuring scarce water saving from interregional virtual water flows in China

NASA Astrophysics Data System (ADS)

Zhao, X.; Li, Y. P.; Yang, H.; Liu, W. F.; Tillotson, M. R.; Guan, D.; Yi, Y.; Wang, H.

2018-05-01

Trade of commodities can lead to virtual water flows between trading partners. When commodities flow from regions of high water productivity to regions of low water productivity, the trade has the potential to generate water saving. However, this accounting of water saving does not account for the water scarcity status in different regions. It could be that the water saving generated from this trade occurs at the expense of the intensified water scarcity in the exporting region, and exerts limited effect on water stress alleviation in importing regions. In this paper, we propose an approach to measure the scarce water saving associated with virtual water trade (measuring in water withdrawal/use). The scarce water is quantified by multiplying the water use in production with the water stress index (WSI). We assessed the scarce water saving/loss through interprovincial trade within China using a multi-region input-output table from 2010. The results show that interprovincial trade resulted in 14.2 km3 of water loss without considering water stress, but only 0.4 km3 scarce water loss using the scarce water concept. Among the 435 total connections of virtual water flows, 254 connections contributed to 20.2 km3 of scarce water saving. Most of these connections are virtual water flows from provinces with lower WSI to that with higher WSI. Conversely, 175 connections contributed to 20.6 km3 of scarce water loss. The virtual water flow connections between Xinjiang and other provinces stood out as the biggest contributors, accounting for 66% of total scarce water loss. The results show the importance of assessing water savings generated from trade with consideration of both water scarcity status and water productivity across regions. Identifying key connections of scarce water saving is useful in guiding interregional economic restructuring towards water stress alleviation, a major goal of China’s sustainable development strategy.

Simulation of the effects of ground-water withdrawals and recharge on ground-water flow in Cape Cod, Martha's Vineyard, and Nantucket Island basins, Massachusetts

USGS Publications Warehouse

Masterson, John P.; Barlow, Paul M.

1994-01-01

The effects of changing patterns of ground-water pumping and aquifer recharge on the surface-water and ground-water hydrologic systems were determined for the Cape Cod, Martha's Vineyard, and Nantucket Island Basins. Three-dimensional, transient, ground-water-flow modelS that simulate both freshwater and saltwater flow were developed for the f1ow cells of Cape Cod which currently have large-capacity public-supply wells. Only the freshwater-flow system was simulated for the Cape Cod flow cells where public-water supply demands are satisfied by small-capacity domestic wells. Two- dimensional, finite-difference, change models were developed for Martha's Vineyard and Nantucket Island to determine the projected drawdowns in response to projected in-season pumping rates for 180 days of no aquifer recharge. Results of the simulations indicate very little change in the position of the freshwater-saltwater interface from predevelopment flow conditions to projected ground-water pumping and recharge rates for Cape Cod in the year 2020. Results of change model simulations for Martha's Vineyard and Nantucket Island indicate that the greatest impact in response to projected in-season ground-water pumping occurs at the pumping centers and the magnitude of the drawdowns are minimal with respect to the total thickness of the aquifers.

Assessment of Effectiveness of Geologic Isolation Systems. Variable thickness transient ground-water flow model. Volume 2. Users' manual

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

Reisenauer, A.E.

1979-12-01

A system of computer codes to aid in the preparation and evaluation of ground-water model input, as well as in the computer codes and auxillary programs developed and adapted for use in modeling major ground-water aquifers is described. The ground-water model is interactive, rather than a batch-type model. Interactive models have been demonstrated to be superior to batch in the ground-water field. For example, looking through reams of numerical lists can be avoided with the much superior graphical output forms or summary type numerical output. The system of computer codes permits the flexibility to develop rapidly the model-required data filesmore » from engineering data and geologic maps, as well as efficiently manipulating the voluminous data generated. Central to these codes is the Ground-water Model, which given the boundary value problem, produces either the steady-state or transient time plane solutions. A sizeable part of the codes available provide rapid evaluation of the results. Besides contouring the new water potentials, the model allows graphical review of streamlines of flow, travel times, and detailed comparisons of surfaces or points at designated wells. Use of the graphics scopes provide immediate, but temporary displays which can be used for evaluation of input and output and which can be reproduced easily on hard copy devices, such as a line printer, Calcomp plotter and image photographs.« less

Estimated Water Flows in 2005: United States

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

Smith, C A; Belles, R D; Simon, A J

2011-03-16

Flow charts depicting water use in the United States have been constructed from publicly available data and estimates of water use patterns. Approximately 410,500 million gallons per day of water are managed throughout the United States for use in farming, power production, residential, commercial, and industrial applications. Water is obtained from four major resource classes: fresh surface-water, saline (ocean) surface-water, fresh groundwater and saline (brackish) groundwater. Water that is not consumed or evaporated during its use is returned to surface bodies of water. The flow patterns are represented in a compact 'visual atlas' of 52 state-level (all 50 states inmore » addition to Puerto Rico and the Virgin Islands) and one national water flow chart representing a comprehensive systems view of national water resources, use, and disposition.« less

Simulation of ground-water flow in the Mojave River basin, California

USGS Publications Warehouse

Stamos, Christina L.; Martin, Peter; Nishikawa, Tracy; Cox, Brett F.

2001-01-01

The proximity of the Mojave River ground-water basin to the highly urbanized Los Angeles region has led to rapid growth in population and, consequently, to an increase in the demand for water. The Mojave River, the primary source of surface water for the region, normally is dry-except for a small stretch of perennial flow and periods of flow after intense storms. Thus, the region relies almost entirely on ground water to meet its agricultural and municipal needs. Ground-water withdrawal since the late 1800's has resulted in discharge, primarily from pumping wells, that exceeds natural recharge. To better understand the relation between the regional and the floodplain aquifer systems and to develop a management tool that could be used to estimate the effects that future stresses may have on the ground-water system, a numerical ground-water flow model of the Mojave River ground-water basin was developed, in part, on the basis of a previously developed analog model. The ground-water flow model has two horizontal layers; the top layer (layer 1) corresponds to the floodplain aquifer and the bottom layer (layer 2) corresponds to the regional aquifer. There are 161 rows and 200 columns with a horizontal grid spacing of 2,000 by 2,000 feet. Two stress periods (wet and dry) per year are used where the duration of each stress period is a function of the occurrence, quantity of discharge, and length of stormflow from the headwaters each year. A steady-state model provided initial conditions for the transient-state simulation. The model was calibrated to transient-state conditions (1931-94) using a trial-and-error approach. The transient-state simulation results are in good agreement with measured data. Under transient-state conditions, the simulated floodplain aquifer and regional aquifer hydrographs matched the general trends observed for the measured water levels. The simulated streamflow hydrographs matched wet stress period average flow rates and times of no flow at the

Coupling Radar Rainfall to Hydrological Models for Water Abstraction Management

NASA Astrophysics Data System (ADS)

Asfaw, Alemayehu; Shucksmith, James; Smith, Andrea; MacDonald, Ken

2015-04-01

The impacts of climate change and growing water use are likely to put considerable pressure on water resources and the environment. In the UK, a reform to surface water abstraction policy has recently been proposed which aims to increase the efficiency of using available water resources whilst minimising impacts on the aquatic environment. Key aspects to this reform include the consideration of dynamic rather than static abstraction licensing as well as introducing water trading concepts. Dynamic licensing will permit varying levels of abstraction dependent on environmental conditions (i.e. river flow and quality). The practical implementation of an effective dynamic abstraction strategy requires suitable flow forecasting techniques to inform abstraction asset management. Potentially the predicted availability of water resources within a catchment can be coupled to predicted demand and current storage to inform a cost effective water resource management strategy which minimises environmental impacts. The aim of this work is to use a historical analysis of UK case study catchment to compare potential water resource availability using modelled dynamic abstraction scenario informed by a flow forecasting model, against observed abstraction under a conventional abstraction regime. The work also demonstrates the impacts of modelling uncertainties on the accuracy of predicted water availability over range of forecast lead times. The study utilised a conceptual rainfall-runoff model PDM - Probability-Distributed Model developed by Centre for Ecology & Hydrology - set up in the Dove River catchment (UK) using 1km2 resolution radar rainfall as inputs and 15 min resolution gauged flow data for calibration and validation. Data assimilation procedures are implemented to improve flow predictions using observed flow data. Uncertainties in the radar rainfall data used in the model are quantified using artificial statistical error model described by Gaussian distribution and

A Low Mach Number Model for Moist Atmospheric Flows

DOE PAGES

Duarte, Max; Almgren, Ann S.; Bell, John B.

2015-04-01

A low Mach number model for moist atmospheric flows is introduced that accurately incorporates reversible moist processes in flows whose features of interest occur on advective rather than acoustic time scales. Total water is used as a prognostic variable, so that water vapor and liquid water are diagnostically recovered as needed from an exact Clausius–Clapeyron formula for moist thermodynamics. Low Mach number models can be computationally more efficient than a fully compressible model, but the low Mach number formulation introduces additional mathematical and computational complexity because of the divergence constraint imposed on the velocity field. Here in this paper, latentmore » heat release is accounted for in the source term of the constraint by estimating the rate of phase change based on the time variation of saturated water vapor subject to the thermodynamic equilibrium constraint. Finally, the authors numerically assess the validity of the low Mach number approximation for moist atmospheric flows by contrasting the low Mach number solution to reference solutions computed with a fully compressible formulation for a variety of test problems.« less