Potential for satellite remote sensing of ground water

Predicting hydrologic behavior at regional scales requires heterogeneous data that are often prohibitively expensive to acquire on the ground. As a result, satellite-based remote sensing has become a powerful tool for surface hydrology. Subsurface hydrology has yet to realize the benefits of remote sensing, even though surface expressions of ground water can be monitored from space. Remotely sensed indicators of ground water may provide important data where practical alternatives are not available. The potential for remote sensing of ground water is explored here in the context of active and planned satellite-based sensors. Satellite technology is reviewed with respect to its ability to measure ground water potential, storage, and fluxes. It is argued here that satellite data can be used if ancillary analysis is used to infer ground water behavior from surface expressions. Remotely sensed data are most useful where they are combined with numerical modeling, geographic information systems, and ground-based information.     

A patch hierarchy approach to modeling surface and subsurface hydrology in complex flood‐plain environments

Geomorphology interacts with surface‐ and ground‐water hydrology across multiple spatial scales. Nonetheless, hydrologic and hydrogeologic models are most commonly implemented at a single spatial scale. Using an existing hydrogeologic computer model, we implemented a simple hierarchical approach to modeling surface‐ and ground‐water hydrology in a complex geomorphic setting. We parameterized the model to simulate ground‐ and surface‐water ?ow patterns through a hierarchical, three‐dimensional, quantitative representation of an anabranched montane alluvial ?ood plain (the Nyack Flood Plain, Middle Fork Flathead River, Montana, USA). Comparison of model results to ?eld data showed that the model provided reasonable representations of spatial patterns of aquifer recharge and discharge, temporal patterns of ?ood‐water storage on the ?ood plain, and rates of ground‐water movement from the main river channel into a large lateral spring channel on the ?ood plain, and water table elevation in the alluvial aquifer. These results suggest that a hierarchical approach to modeling ground‐ and surface‐water hydrology can reproduce realistic patterns of surface‐ and ground‐water ?ux on alluvial ?ood plains, and therefore should provide an excellent ‘quantitative laboratory’ for studying complex interactions between geomorphology and hydrology at and across multiple spatial scales. Copyright © 2004 John Wiley & Sons, Ltd.     

An Imminent Human Resource Crisis in Ground Water Hydrology?

Anecdotal evidence, mostly from the United States, suggests that it has become increasingly difficult to find well-trained, entry-level ground water hydrologists to fill open positions in consulting firms and regulatory agencies. The future prospects for filling positions that require training in ground water hydrology are assessed by considering three factors: the market, the numbers of qualified students entering colleges and universities, and the aging of the existing workforce. The environmental and water resources consulting industry has seen continuous albeit variable growth, and demand for environmental scientists and hydrologists is expected to increase significantly. Conversely, students' interest and their enrollment in hydrology and water resources programs have waned in recent years, and the interests of students within these departments have shifted away from ground water hydrology in some schools. This decrease in the numbers of U.S. students graduating in hydrology or emphasizing ground water hydrology is coinciding with the aging of and pending retirement of ground water scientists and engineers in the baby boomer generation. We need to both trigger the imagination of students at the elementary school level so that they later want to apply science and math and communicate the career opportunities in ground water hydrology to those high school and college graduates who have acquired the appropriate technical background. Because the success of a consulting firm, research organization, or regulatory agency is derived from the skills and judgment of the employees, human resources will be an increasingly more critical strategic issue for many years.     

Seismic safety analysis of tall concrete dams,investigation and insights on critical challenges

This paper discusses critical and potentially controversial issues related to the seismic safety of tall concrete dams. These include the seismic input at a dam site, the effective treatment of the damage-rupture process, and the consideration of compressibility of reservoir water for hydrodynamic pressure. Major challenges to currently popular but questionable treatments of these critical problems are presented. Insights and additional research on these critical challenges are emphasized and explained based on prior published works of the author. More reasonable alternatives to dealing with these potentially controversial problems are provided in light of engineering practice in China. First, the design seismic input at depth as deconvoluted from an arbitrarily selected recorded accelerogram at a control point of an artificially developed free-field surface with the elevation of the dam crest is difficult for engineering projects to accept as appropriate. It may be more reasonable to use the design seismic incident motions as half of the ground surface motions from seismic safety analyses obtained from deterministic or probability approaches conducted by seismologists according to approved standards or guidelines. Second, since seismic damage to the dam must be estimated separately following uniaxial tensile and compressive experimental damage evolution rules, a simplified and realistic nonlinear elastic model is proposed as an alternative to the plastic-damage coupling model, which is very complex and includes assumptions based on a number of uncertainties. Finally, the effect of the reflection coefficient for compressibility of reservoir water on hydrodynamic pressures is very sensitive. The notion that the applied unified reflection coefficient at the reservoir bottom could be frequency-dependent and exhibit a significant variability in space as confirmed by field tests is questionable. To neglect the compressibility of reservoir water it may be closer to engineering practice at present.     

Influence of time‐varying frequency content in earthquake ground motions on seismic response of linear elastic systems

Earthquake ground motion records are nonstationary in both amplitude and frequency content. However, the latter nonstationarity is typically neglected mainly for the sake of mathematical simplicity. To study the stochastic effects of the time‐varying frequency content of earthquake ground motions on the seismic response of structural systems, a pair of closely related stochastic ground motion models is adopted here. The first model (referred to as ground motion model I) corresponds to a fully nonstationary stochastic earthquake ground motion model previously developed by the authors. The second model (referred to as ground motion model II) is nonstationary in amplitude only and is derived from the first model. Ground motion models I and II have the same mean‐square function and global frequency content but different features of time variation in the frequency content, in that no time variation of the frequency content exists in ground motion model II. New explicit closed‐form solutions are derived for the response of linear elastic SDOF and MDOF systems subjected to stochastic ground motion model II. New analytical solutions for the evolutionary cross‐correlation and cross‐PSD functions between the ground motion input and the structural response are also derived for linear systems subjected to ground motion model I. Comparative analytical results are presented to quantify the effects of the time‐varying frequency content of earthquake ground motions on the structural response of linear elastic systems. It is found that the time‐varying frequency content in the seismic input can have significant effects on the stochastic properties of system response. Copyright © 2016 John Wiley & Sons, Ltd.     

Numerical implementation of a backward probabilistic model of ground water contamination

Backward location and travel time probabilities can be used to characterize known and unknown sources or prior positions of ground water contamination. Backward location probability describes the position of the observed contamination at some time in the past; backward travel time probability describes the amount of time prior to observation that the contamination was released from its source or was at a particular upgradient location. The governing equation for backward probabilities is the adjoint of the governing equation for contaminant transport, but with new load terms. Numerical codes that have been written to solve the forward equations of contaminant transport, e.g., the advection-dispersion equation, can also be used to solve the adjoint equation for location and travel time probabilities; however, the interpretation of the results is different and some new approximations must be made for the load terms. We present the governing equations for backward location and travel time probabilities, and provide appropriate numerical approximations for these load terms using the cell-centered finite difference method, one of the most popular numerical methods in ground water hydrology. We discuss some additional numerical considerations for the backward model including boundary conditions, reversal of the flow field, and interpretation of the results. We illustrate the implementation of the backward probability model using hypothetical examples in one- and two-dimensional domains. We also present a three-dimensional application of a pump-and-treat remediation capture zone delineation at the Massachusetts Military Reservation. The illustrations are performed using MODFLOW-96 for flow simulations and MT3DMS for transport simulations.     

Regional gravity variations in Europe from superconducting gravimeters

Recent satellite missions (CHAMP, GRACE) are now returning data on the time variation of the gravity field with harmonic coefficients computed every 4 weeks. The promise is to achieve a sub-microgal accuracy that will define continental mass variations involving large-scale hydrology. With this in mind, we examine the time varying gravity field over central Europe using a limited number of high quality ground-based superconducting gravimeter stations within the Global Geodynamics Project (GGP). Our purpose is to see whether there are coherent signals between the individual stations and to compare the regional component with that predicted from models of continental hydrology. The results are encouraging. We have found, using empirical orthogonal eigenfunctions of the gravity data that a clear annual signal is present that is consistent in phase (low amplitudes in summer) and amplitude (1–3 microgal) with that determined from a large-scale model of land water in connection with global climate modeling. More work is required to define how the gravity field is related to large-scale soil moisture and other mass variations, and we have yet to compare our results to the latest satellite-derived data.     

Identification of physical processes inherent in artificial neural network rainfall runoff models

The emergence of artificial neural network (ANN) technology has provided many promising results in the field of hydrology and water resources simulation. However, one of the major criticisms of ANN hydrologic models is that they do not consider/explain the underlying physical processes in a watershed, resulting in them being labelled as black‐box models. This paper discusses a research study conducted in order to examine whether or not the physical processes in a watershed are inherent in a trained ANN rainfall‐runoff model. The investigation is based on analysing definite statistical measures of strength of relationship between the disintegrated hidden neuron responses of an ANN model and its input variables, as well as various deterministic components of a conceptual rainfall‐runoff model. The approach is illustrated by presenting a case study for the Kentucky River watershed. The results suggest that the distributed structure of the ANN is able to capture certain physical behaviour of the rainfall‐runoff process. The results demonstrate that the hidden neurons in the ANN rainfall‐runoff model approximate various components of the hydrologic system, such as infiltration, base flow, and delayed and quick surface flow, etc., and represent the rising limb and different portions of the falling limb of a flow hydrograph. Copyright © 2004 John Wiley & Sons, Ltd.     

Modeling the Usefulness of Spatial Correlation Analysis on Karst Systems

Cross-correlation analyses on field data collected in karst aquifer systems can be used to develop a conceptual understanding of the aquifer. This includes the use of many data sets from the same aquifer to develop an understanding of how properties vary spatially. We focus on a method for characterizing the distribution of recharge, which is becoming increasingly important in regions where urban development encroaches on these important sources of water. Spatially varying precipitation data and cross-correlation analysis provide a means of spatially characterizing recharge locations on a karst aquifer. Our work expands on the numerical experiments conducted by Padilla and Pulido-Bosch (1995) using the numerical ground water model MODFLOW to introduce spatially varying parameters. The numerical experiments include conduit-controlled, matrix-controlled, and mixed karst systems with more than one precipitation time series input. The results show that spatially varying parameters can be inferred based on the cross-correlation of precipitation data and spring discharge. Simulations were completed using aquifer parameters derived from studies of the Barton Springs segment of the Edwards Aquifer. The simulations indicate that spatial variability within an aquifer can be inferred using cross-correlation analysis. A field study using these methods is summarized for Barton Springs near Austin, Texas.     

Interaction of Aquifer and River‐Canal Network near Well Field

The article presents semi‐analytical mathematical models to asses (1) enhancements of seepage from a canal and (2) induced flow from a partially penetrating river in an unconfined aquifer consequent to groundwater withdrawal in a well field in the vicinity of the river and canal. The nonlinear exponential relation between seepage from a canal reach and hydraulic head in the aquifer beneath the canal reach is used for quantifying seepage from the canal reach. Hantush's (1967) basic solution for water table rise due to recharge from a rectangular spreading basin in absence of pumping well is used for generating unit pulse response function coefficients for water table rise in the aquifer. Duhamel's convolution theory and method of superposition are applied to obtain water table position due to pumping and recharge from different canal reaches. Hunt's (1999) basic solution for river depletion due to constant pumping from a well in the vicinity of a partially penetrating river is used to generate unit pulse response function coefficients. Applying convolution technique and superposition, treating the recharge from canal reaches as recharge through conceptual injection wells, river depletion consequent to variable pumping and recharge is quantified. The integrated model is applied to a case study in Haridwar (India). The well field consists of 22 pumping wells located in the vicinity of a perennial river and a canal network. The river bank filtrate portion consequent to pumping is quantified.     

Seismic response of multi‐span simply supported bridges to a spatially varying earthquake ground motion

This paper carries out a parametrical study of the pounding phenomenon associated with the seismic response of multi‐span simply supported bridges with base isolation devices. In particular, the analyses focus on the causal relationship between pounding and the properties of a spatially varying earthquake ground motion. In order to include the effect of the torsional component of pounding forces on the seismic response of the whole structure, a three‐dimensional (3D) finite element model has been defined and 3D non‐linear time‐history analyses have been performed. A parametrical study on the size of the gaps between adjacent bridge decks has highlighted that the pounding effects are amplified when the spatially varying ground motion time histories at each support are considered. Because of a spatially varying input, the pounding forces can assume values 3–4 times larger than those derived by a conventional seismic analysis with uniform input or with spatial input but considering ground motion wave passage effect only. The numerical results show that in order to achieve an acceptably safe structural performance during seismic events, a correct design of the isolation devices should take into account the relative displacements calculated by means of a non‐linear time‐history analysis with multi‐support excitation. Copyright © 2002 John Wiley & Sons, Ltd.     

Tritium balance in macro‐scale river basins analysed through distributed hydrological modelling

The recession of bomb tritium in river discharge of large basins indicates a contribution of slowly moving water. For an appropriate interpretation it is necessary to consider different runoff components (e.g. direct runoff and ground water components) and varying residence times of tritium in these components. The spatially distributed catchment model (tracer aided catchment model, distributed; TAC^D) and a tritium balance model (TRIBIL) were combined to model process‐based tritium balances in a large German river basin (Weser 46 240 km²) and seven embedded sub‐basins. The hydrological model (monthly time step, 2 × 2 km²) estimated the three major runoff components: direct runoff, fast‐moving and slow‐moving ground water for the period of 1950 to 1999. The model incorporated topography, land use, geomorphology, geology and hydro‐meteorological data. The results for the different basins indicated a contribution of direct runoff of 30–50% and varying amounts for fast and slow ground water components. Combining these results with the TRIBIL model allowed us to estimate the residence time of the components. Mean residence times of 8 to 14 years were found for the fast ground water component, 21 to 93 years for the slow ground water component and 14 to 50 years for an overall mean residence time within these basins. Balance calculations for the Weser basin indicate an over‐estimation of loss of tritium through evapotranspiration (more than 60%) and decay (10%). About 28% were carried in stream‐flow where direct runoff contributed about 12% and ground water runoff 13% in relation to precipitation input over the studied 50‐year period. Neighbouring basins and nuclear power plants contributed about 1% each over this time period. Copyright © 2007 John Wiley & Sons, Ltd.     

Response of multiply supported rigid plate to spatially correlated seismic excitations

The ground motions produced by an earthquake at the multiple support points of large structures can differ considerably. In this study, the quasi-static, the dynamic and the total structural responses of a multiply supported rigid plate to spatially varying multiple ground motion inputs are analysed. The results, together with those obtained by neglecting ground motion phase shifts and those by neglecting ground motion coherency losses, are compared with the corresponding structural responses to a single input. The comparisons show that the single input method overestimates translational responses and underestimates rotational responses. They also show that the responses are sometimes overestimated and sometimes underestimated by neglecting the coherency loss effects between the multiple ground motion excitations.     

Using semivariogram parameter uncertainty in hydrogeological applications

Geostatistical estimation (kriging) and geostatistical simulation are routinely used in ground water hydrology for optimal spatial interpolation and Monte Carlo risk assessment, respectively. Both techniques are based on a model of spatial variability (semivariogram or covariance) that generally is not known but must be inferred from the experimental data. Where the number of experimental data is small (say, several tens), as is not unusual in ground water hydrology, the model fitted to the empirical semivariogram entails considerable uncertainty. If all the practical results are based on this unique fitted model, the final results will be biased. We propose that, instead of using a unique semivariogram model, the full range of models that are inside a given confidence region should be used, and the weight that each semivariogram model has on the final result should depend on its plausibility. The first task, then, is to evaluate the uncertainty of the model, which can be efficiently done by using maximum likelihood inference. The second task is to use the range of plausible models in applications and to show the effect observed on the final results. This procedure is put forth here with kriging and simulation applications, where the uncertainty in semivariogram parameters is propagated into the final results (e.g., the prediction of ground water head). A case study using log-transmissivity data from the Vega de Granada aquifer, in southern Spain, is given to illustrate the methodology.     

Influence of injection conditions on field tracer experiments

Calibration of ground water transport models is often performed using results of field tracer experiments. However, little attention is usually paid to the influence, on resulting breakthrough curves, of injection conditions and well-aquifer interactions, more particularly of the influence of the possible trapping of the tracer in the injection wellbore. Recently, a new mathematical and numerical approach has been developed to model injection conditions and well-aquifer interactions in a very accurate way. Using an analytical solution derived from this model, a detailed analysis is made of the evolution of the tracer input function in the aquifer. By varying injection conditions from one simulation to another, synthetic breakthrough curves are generated with the SUFT3D ground water flow and transport finite-element simulator. These tests show clearly that the shape of the breakthrough curves can be dramatically affected by injection conditions. Using generated breakthrough curves as "actual" field results, a calibration of hydrodispersive parameters is performed, neglecting the influence of injection conditions. This shows that neglecting the influence of actual injection conditions can lead to (1) errors on fitted parameters and (2) misleading identification of the active transport processes. Conclusions and guidelines are drawn in terms of proposed methodologies for better controlling the tracer injection in the field, in order to minimize risk of misinterpretation of results.     

中国大陆地电场日变化特征及影响因素

作为地电场最重要的背景变化组成成分之一,地电场日变化普遍存在。本文选取中国大陆10个地电场台站观测数据,应用时序迭加和谐波分析方法,研究地电场日变化的时间域和空间域基本特征。地电场日变化大多呈现"两峰一谷"或"两谷一峰"的变化特征,日变化特征一年四季基本相同,但是幅度、相位略有差异。此外,利用地磁、地貌、电性结构等资料,对地电场日变化的复杂性及与空间电磁活动、地下电性结构及潮汐相关性进行分析,认为地电场日变化特征本质是受空间电磁活动和地下介质电性结构、水文、地质等因素共同作用的结果。     

In Situ Measurement of Electro-Osmotic Fluxes and Conductivity Using Single Wellbore Tracer Tests

Electro-osmosis (EO), the movement of water through porous media in response to an electric field, offers a means for extracting contaminated ground water from fine-grained sediments, such as clays, that are not easily amenable to conventional pump-and-treat approaches. The EO-induced water flux is proportional to the voltage gradient in a manner analogous to the flux dependence on the hydraulic gradient under Darcy's law. The proportionality constant, the soil electro-osmotic conductivity or k_eo, is most easily measured in soil cores using bench-top tests, where flow is one-dimensional and interfering effects attributable to Darcy's law can be directly accounted for. In contrast, quantification of EO fluxes and k_eo in the field under deployment conditions can be difficult because electrodes are placed in ground water wells that may be screened across a heterogeneous mixture of lithologies. As a result, EO-induced water fluxes constitute an approximate radial flow system that is superimposed upon a Darcy flow regime through permeable pathways that may or may not be coupled with hydraulic head differences created by the EO-induced water fluxes. A single well comparative tracer test, which indirectly measures EO fluxes by comparing wellbore tracer dilution rates between background and EO-induced water fluxes, may provide a means for routinely quantifying the efficacy of EO systems in such settings. EO fluxes measured in field tests through this technique at a ground water contamination site were used to estimate a mean k_eo value through a semianalytic line source model of the electric field. The resulting estimate agrees well with values reported in the literature and with values obtained with bench-top tests conducted on a soil core collected in the test area.     

Using Orthogonal Array Sampling to Cope with Uncertainty in Ground Water Problems

Uncertainty in ground water hydrology originates from different sources. Neglecting uncertainty in ground water problems can lead to incorrect results and misleading output. Several approaches have been developed to cope with uncertainty in ground water problems. The most widely used methods in uncertainty analysis are Monte Carlo simulation (MCS) and Latin hypercube sampling (LHS), developed from MCS. Despite the simplicity of MCS, many runs are required to achieve a reliable result. This paper presents orthogonal array (OA) sampling as a means to cope with uncertainty in ground water problems. The method was applied to an analytical stream depletion problem. To examine the convergence rate of the OA sampling, the results were compared to MCS and LHS. This study shows that OA can be applied to ground water problems. Results reveal that the convergence rate of the OA sampling is faster than MCS and LHS, with a smaller error of estimate when applied to a stream depletion problem.     

Kinematic wave modelling in water resources: a historical perspective

The history of the kinematic wave theory and its applications in water resources are traced. It is shown that the theory has found its niche in water resources and its applications are so widespread that they may well constitute what may be termed ‘kinematic wave hydrology’. Few theories have been applied in hydrology and water resources as extensively as the kinematic wave theory. This theory, however, is not without limitations and when it is applied they must be so recognized. Copyright © 2001 John Wiley & Sons, Ltd.