USGS Science⌕ Search

SEARCH · USGS Science

Results for “Research Paper”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

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

At least 1,621 records · Page 90Linked to original sources
USGS1994

Petrophysical analysis of geophysical logs of the National Drilling Company-U.S. Geological Survey ground-water research project for Abu Dhabi Emirate, United Arab Emirates

A program of borehole-geophysical logging was implemented to supply geologic and geohydrologic information for a regional ground-water investigation of Abu Dhabi Emirate. Analysis of geophysical logs was essential to provide information on geohydrologic properties because drill cuttings were not always adequate to define lithologic boundaries. The standard suite of logs obtained at most project test holes consisted of caliper, spontaneous potential, gamma ray, dual induction, microresistivity, compensated neutron, compensated density, and compensated sonic. Ophiolitic detritus from the nearby Oman Mountains has unusual petrophysical properties that complicated the interpretation of geophysical logs. The density of coarse ophiolitic detritus is typically greater than 3.0 grams per cubic centimeter, porosity values are large, often exceeding 45 percent, and the clay fraction included unusual clays, such as lizardite. Neither the spontaneous-potential log nor the natural gamma-ray log were useable clay indicators. Because intrinsic permeability is a function of clay content, additional research in determining clay content was critical. A research program of geophysical logging was conducted to determine the petrophysical properties of the shallow subsurface formations. The logging included spectral-gamma and thermal-decay-time logs. These logs, along with the standard geophysical logs, were correlated to mineralogy and whole-rock chemistry as determined from sidewall cores. Thus, interpretation of lithology and fluids was accomplished. Permeability and specific yield were calculated from geophysical-log data and correlated to results from an aquifer test. On the basis of results from the research logging, a method of lithologic and water-resistivity interpretation was developed for the test holes at which the standard suite of logs were obtained. In addition, a computer program was developed to assist in the analysis of log data. Geohydrologic properties were estimated, including volume of clay matrix, volume of matrix other than clay, density of matrix other than clay, density of matrix, intrinsic permeability, specific yield, and specific storage. Geophysical logs were used to (1) determine lithology, (2) correlate lithologic and permeable zones, (3) calibrate seismic reprocessing, (4) calibrate transient-electromagnetic surveys, and (5) calibrate uphole-survey interpretations. Logs were used at the drill site to (1) determine permeability zones, (2) determine dissolved-solids content, which is a function of water resistivity, and (3) design wells accordingly. Data and properties derived from logs were used to determine transmissivity and specific yield of aquifer materials.

Water Supply Paper↗
USGS2011

Nine endangered taxa, one recovering ecosystem: Identifying common ground for recovery on Santa Cruz Island, California

It is not uncommon to have several rare and listed taxa occupying habitats in one landscape or management area where conservation amounts to defense against the possibility of further loss. It is uncommon and extremely exciting, however, to have several listed taxa occupying one island that is managed cooperatively for conservation and recovery. On Santa Cruz Island, the largest of the northern California island group in the Santa Barbara Channel, we have a golden opportunity to marry ecological knowledge and institutional "good will" in a field test of holistic rare plant conservation. Here, the last feral livestock have been removed, active weed control is underway, and management is focused on understanding and demonstrating system response to conservation management. Yet funding limitations still exist and we need to plan the most fiscally conservative and marketable approach to rare plant restoration. We still experience the tension between desirable quick results and the ecological pace of system recovery. Therefore, our research has focused on identifying fundamental constraints on species recovery at individual, demographic, habitat, and ecosystem levels, and then developing suites of actions that might be taken across taxa and landscapes. At the same time, we seek a performance middle ground that balances an institutional need for quick demonstration of hands-on positive results with a contrasting approach that allows ecosystem recovery to facilitate species recovery in the long term. We find that constraints vary across breeding systems, life-histories, and island locations. We take a hybrid approach in which we identify several actions that we can take now to enhance population size or habitat occupancy for some taxa by active restoration, while allowing others to recover at the pace of ecosystem change. We make our recommendations on the basis of data we have collected over the last decade, so that management is firmly grounded in ecological observation.

California↗
USGS1977

Cauldron subsidence of Oligocene age at Mount Lewis, Shoshone Range, Nevada: A reasonable interpretation

James Gilluly has rejected the interpretation of Wrucke and Silberman (U.S. Geol. Survey Prof. Paper 876, 1975) that a thrust fault and tear fault mapped by Gilluly and Gates (U.S. Geol. Survey Prof. Paper 465, 1965) as structures bounding the upper plate of the Roberts Mountains thrust at Mount Lewis are parts of a ring fracture around an area that underwent volcanic collapse. In his discussion (this volume) of our paper, Gilluly fails to consider important questions that we presented in support of the subsidence hypothesis. Instead of answering these critical questions, Gilluly merely recapitulates the interpretations that he and Gates gave in Professional Paper 465. We presented new information, including a map of one critical area along the cauldron boundary where, among other significant differences in geologic interpretation, we found the ring fault where previously no steep fault was shown. Gilluly believes that the paucity of dikes along the ring fracture at Mount Lewis is highly anomalous for cauldrons. However, the amount of dike rock is comparable to that in known cauldrons (some cauldrons have none) and is what might be expected at high levels in subsidence structures that have undergone relatively little resurgent igneous activity after collapse. Gilluly concludes that in formulating our interpretation of volcanic collapse, we have ignored much evidence that he and Gates have presented on thrust faulting; in the Shoshone Range. On the contrary, we have considered their ideas and have reinterpreted them using new evidence that strongly supports the concept of cauldron subsidence at Mount Lewis.

Nevada↗
USGS1996

Methods and techniques for rigorous calibration of a cellular automaton model of urban growth

Several lessons about the process of calibration were learned during development of a self-modifying cellular automaton model to predict urban growth. This model, part of a global change research project on human-induced land transformations, was used to predict the spatial extent of urban growth 100 years into the future. The context of the prediction was to evaluate urban environmental disturbances such as land use conversion, urban heat island intensification, and greenhouse gas generation. Using data for the San Francisco Bay area as a test case, methods were developed, including interactive and statistical versions of the model, animation and visualization tools, automated testing methods, and Monte Carlo simulations. This presentation will enumerate, analyze, and discuss the lessons learned during the extensive process of model calibration. Experience with the methods developed may have broader use in assisting the rigorous calibration for other CA models, and perhaps those coupled environmental models with an extensive spatial data component. These methods are now under test as the project moves to a new data set for the Washington, D.C.-Baltimore area.

Conference Paper↗
USGS2003

Stochastic analysis of transverse dispersion in density‐coupled transport in aquifers

Spectral perturbation techniques have been used previously to derive integral expressions for dispersive mixing in concentration‐dependent transport in three‐dimensional, heterogeneous porous media, where fluid density and viscosity are functions of solute concentration. Whereas earlier work focused on evaluating longitudinal dispersivity in isotropic media and incorporating the result in a mean one‐dimensional transport model, the emphasis of this paper is on evaluation of the complete dispersion tensor, including the more general case of anisotropic media. Approximate analytic expressions for all components of the macroscopic dispersivity tensor are derived, and the tensor is shown to be asymmetric. The tensor is separated into its symmetric and antisymmetric parts, where the symmetric part is used to calculate the principal components and principal directions of dispersivity, and the antisymmetric part of the tensor is shown to modify the velocity of the solute body compared to that of the background fluid. An example set of numerical simulations incorporating the tensor illustrates the effect of density‐coupled dispersivity on a sinking plume in an aquifer. The simulations show that the effective transverse vertical spreading in a sinking plume to be significantly greater than would be predicted by a standard density‐coupled transport model that does not incorporate the coupling in the dispersivity tensor.

Water Resources Research↗
USGS2015

Beyond annual streamflow reconstructions for the Upper Colorado River Basin: a paleo-water-balance approach

In this paper, we present a methodology to use annual tree-ring chronologies and a monthly water balance model to generate annual reconstructions of water balance variables (e.g., potential evapotrans- piration ( PET ), actual evapotranspiration ( AET ), snow water equivalent ( SWE ), soil moisture storage ( SMS ), and runoff ( R )). The method involves resampling monthly temperature and precipitation from the instrumental record directed by variability indicated by the paleoclimate record. The generated time series of monthly temperature and precipitation are subsequently used as inputs to a monthly water balance model. The methodology is applied to the Upper Colorado River Basin, and results indicate that the methodology reliably simulates water-year runoff, maximum snow water equivalent, and seasonal soil moisture storage for the instrumental period. As a final application, the methodology is used to produce time series of PET , AET , SWE , SMS , and R for the 1404–1905 period for the Upper Colorado River Basin.

Water Resources Research↗
USGS2010

Kinetically influenced terms for solute transport affected by heterogeneous and homogeneous classical reactions

Simulation of transport affected by heterogeneous or homogeneous reversible reactions requires a choice between local equilibrium-based and kinetics-based models. The error associated with the use of equilibrium-based models is equivalent to the error of neglecting certain mathematical terms in the governing kinetics-based transport equations. Identification and evaluation of these kinetically influenced terms can therefore aid in the development of criteria for applicability of local equilibrium- based transport models. This paper extends a four-step derivation procedure, previously presented for cases of transport affected by surface reactions, to transport problems involving homogeneous reactions (solution phase complex formation or oxidation-reduction) and/or precipitation-dissolution reactions. Derivations for these classes of reactions are used to illustrate the manner in which mathematical differences between reaction classes are reflected in the mathematical derivation procedures required to identify kinetically influenced terms. Simulation results for a case of transport affected by a single solution phase complexation reaction and for a case of transport affected by a precipitation-dissolution reaction are used to demonstrate the nature of departures from equilibrium-controlled transport as well as the use of kinetically influenced terms in determining criteria for the applicability of the local equilibrium assumption. A final derivation for a multireaction problem demonstrates the application of the generalized procedure to a case of transport affected by reactions of several classes.

Water Resources Research↗
USGS2010

Apparent dispersion in transient groundwater flow

This paper investigates the effects of large-scale temporal velocity fluctuations, particularly changes in the direction of flow, on solute spreading in a two-dimensional aquifer. Relations for apparent longitudinal and transverse dispersivity are developed through an analytical solution for dispersion in a fluctuating, quasi-steady uniform flow field, in which storativity is zero. For transient flow, spatial moments are evaluated from numerical solutions. Ignored or unknown transients in the direction of flow primarily act to increase the apparent transverse dispersivity because the longitudinal dispersivity is acting in a direction that is not the assumed flow direction. This increase is a function of the angle between the transient flow vector and the assumed steady state flow direction and the ratio of transverse to longitudinal dispersivity. The maximum effect on transverse dispersivity occurs if storativity is assumed to be zero, such that the flow field responds instantly to boundary condition changes.

Water Resources Research↗
USGS1965

Identification and measurement of chlorinated organic pesticides in water by electron-capture gas chromatography

Pesticides, in minute quantities, may affect the regimen of streams, and because they may concentrate in sediments, aquatic organisms, and edible aquatic foods, their detection and their measurement in the parts-per-trillion range are considered essential. In 1964 the U.S. Geological Survey at Menlo Park, Calif., began research on methods for monitoring pesticides in water. Two systems were selected--electron-capture gas chromatography and microcoulometric-titration gas chromatography. Studies on these systems are now in progress. This report provides current information on the development and application of an electron-capture gas chromatographic procedure. This method is a convenient and extremely sensitive procedure for the detection and measurement of organic pesticides having high electron affinities, notably the chlorinated organic pesticides. The electron-affinity detector is extremely sensitive to these substances but it is not as sensitive to many other compounds. By this method, the chlorinated organic pesticide may be determined on a sample of convenient size in concentrations as low as the parts-per-trillion range. To insure greater accuracy in the identifications, the pesticides reported were separated and identified by their retention times on two different types of gas chromatographic columns.

Water Supply Paper↗
USGS2022

Observations on Whooping Crane parental provisioning of chicks

Crane chicks are dependent on parent birds for provisioning during the first few months of life, but no study has examined this provisioning in detail. In 2014 research staff at the U.S. Geological Survey, Eastern Ecological Science Center (formerly Patuxent Wildlife Research Center), in Laurel, Maryland, made multiple observations of parent whooping cranes ( Grus americana ) feeding or interacting with their chick during the 3 months from hatching to fledging. Both parents participated in the feeding of the chick and only 1 chick was raised by each pair of parent whooping cranes. Initially feeding frequency was low (0-20 times per hr), but as the chick absorbed its yolk sac and required food, feeding frequency increased to a high of 105 times per hour. Whooping crane parents fed their chick from 0 to 105 times per hour. Feeding frequency peaked around day 19, then decreased after the chicks reached 40 days of age but continued at a low level during the entire 3 months from hatch to fledging. Because feeding frequency observed for this study was very low at fledging, the use of feeding by alloparents as a measure of chick-alloparent bond may not be practical.

Conference Paper↗
USGS2010

Conjunctive use of groundwater and surface water for irrigated agriculture: Risk aversion

In examining the South Platte system in Colorado where surface water and groundwater are used conjunctively for irrigation, we find the actual installed well capacity is approximately sufficient to irrigate the entire area. This would appear to be an overinvestment in well capacity. In this paper we examine to what extent groundwater is being developed as insurance against periods of low streamflow. Using a simulation model which couples the hydrology of a conjunctive stream aquifer system to a behavioral-economic model which incorporates farmer behavior in such a system, we have investigated the economics of an area patterned after a reach of the South Platte Valley in Colorado. The results suggest that under current economic conditions the most reasonable groundwater pumping capacity is a total capacity capable of irrigating the available acreage with groundwater. Installing sufficient well capacity to irrigate all available acreage has two benefits: (1) this capacity maximizes the expected net benefits and (2) this capacity also minimizes the variation in annual income: it reduces the variance to essentially zero. As pumping capacity is installed in a conjunctive use system, the value of flow forecasts is diminished. Poor forecasts are compensated for by pumping groundwater.

Colorado↗
USGS2008

The Role of Eolian Sediment in the Preservation of Archeologic Sites Along the Colorado River Corridor in Grand Canyon National Park, Arizona

Since the closure of Glen Canyon Dam in 1963, the natural hydrologic and sedimentary systems along the Colorado River in the Grand Canyon reach have changed substantially (see, for example, Andrews, 1986; Johnson and Carothers, 1987; Webb and others, 1999b; Rubin and others, 2002; Topping and others, 2003; Wright and others, 2005; Hazel and others, 2006b). The dam has reduced the fluvial sediment supply at the upstream boundary of Grand Canyon National Park by about 95 percent. Regulation of river discharge by dam operations has important implications for the storage and redistribution of sediment in the Colorado River corridor. In the absence of floods, sediment is not deposited at elevations that regularly received sediment before dam closure. Riparian vegetation has colonized areas at lower elevations than in predam time when annual floods removed young vegetation (Turner and Karpiscak, 1980). Together, these factors have caused a systemwide decrease in the size and number of subaerially exposed fluvial sand deposits since the 1960s, punctuated by episodic aggradation during the exceptional high-flow intervals in 1983-84, 1996, and 2004 and by sediment input from occasional tributary floods (Beus and others, 1985; Schmidt and Graf, 1987; Kearsley and others, 1994; Hazel and others, 1999; Schmidt and others, 2004; Wright and others, 2005). When the Bureau of Reclamation sponsored the creation of the Glen Canyon Environmental Studies (GCES) research initiative in 1982, research objectives included physical and biologic resources, whereas the effects of dam operations on cultural resources were not addressed (Fairley and others, 1994; Fairley, 2003). In the early 1980s, it was widely believed that because few archeologic sites were preserved within the river's annual-flood zone, cultural features would not be greatly affected by dam operations. Recent studies, however, indicate that alterations in the flow and sediment load of the Colorado River by Glen Canyon Dam operations may affect archeologic sites within the river corridor, even above the annual flood limit (Hereford and others, 1993, Yeatts, 1996, 1997; Thompson and Potochnik, 2000; Draut and others, 2005). (The annual-flood zone is defined here by the mean annual predam flood of 2,410 m3/s; the 'predam flood limit', the highest elevation at which fluvial deposits are present locally, was approximately equivalent to a rare, major flood of 8,500 m3/s; Topping and others, 2003.) Of about 500 archeologic sites documented between Glen Canyon Dam and Separation Canyon (255 river miles), more than 330 are considered to be within the area of potential effect (APE) of dam operations (Fairley and others, 1994; Neal and others, 2000; Fairley, 2005). The APE was designated by the National Park Service (NPS) to include the area below the peak stage of the 1884 flood; though previously believed to have reached 8,490 m3/s, this flood was shown by Topping and others (2003) to have peaked at 5,940 m3/s. Archeologic research and monitoring in Grand Canyon National Park focus increasingly on the potential effects of Glen Canyon Dam operations on the landscape in which these sites exist. Many archeologic sites in or on sedimentary deposits are being eroded, owing to eolian deflation and gully incision (Leap and others, 2000; Neal and others, 2000; Fairley, 2003, 2005). Hereford and others (1993) first suggested that gully incision of sedimentary deposits, and the base level to which small drainage systems respond, were linked to dam operations; they hypothesized that pronounced arroyo incision, which occurs during rainfall runoff, was caused by lowering of the effective base level at the mouths of ephemeral drainages to meet the new postdam elevation of high-flow sediment deposition, about 3 to 4 m below the lowest predam alluvial terraces. Thompson and Potochnik (2000) modified that hypothesis to include the restorative effects of fluvial deposition in the mouths of gullies and ar

Professional Paper↗
USGS2003

Geothermal GIS coverage of the Great Basin, USA: Defining regional controls and favorable exploration terrains

A geographic information system (GIS) of geothermal resources, built last year for the state of Nevada, is being expanded to cover the Great Basin, USA. Data from that GIS is being made available to industry, other researchers, and the public via a web site at the Great Basin Center for Geothermal Energy, Reno, Nevada. That web site features a search engine, supports ArcExplorer?? for on-line map construction, and provides downloadable data layers in several formats. Though data collection continues, preliminary analysis has begun. Contour maps of geothermal temperatures, constructed using geothermometer temperatures calculated from a Great Basin geochemical database compiled by the Geo-Heat Center, reveal distinctive trends and patterns. As expected, magmatic-type and extensional-type geothermal systems have profoundly different associations, with magmatic-type systems following major tectonic boundaries, and extensional-type systems associating with regionally high heat flow, thin crust, active faulting, and high extensional strain rates. As described by earlier researchers, including Rowen and Wetlaufer (1981) and Koenig and McNitt (1983), high-temperature (> 100??C) geothermal systems appear to follow regional northeast trends, most conspicuously including the Humboldt structural zone in Nevada, the "Black Rock-Alvord Desert" trend in Oregon and Nevada, and the "Newcastle-Roosevelt" trend in Utah and Nevada. Weights-of-evidence analyses confirm a preference of high-temperature geothermal systems for young northeast-trending faults, but the distribution of geothermal systems correlates even better with high rates of crustal extension, as measured from global positioning system (GPS) stations in Nevada. A predictive map of geothermal potential based only on areas of high extensional strain rates and high heat flux does an excellent job of regionally predicting the location of most known geothermal systems in Nevada, and may prove useful in identifying blind systems.

Conference Paper↗
USGS2021

Learning augmented methods for matching: Improving invasive species management and urban mobility

With the success of machine learning, integrating learned models into real-world systems has become a critical chal- lenge. Naively applying predictions to combinatorial opti- mization problems can incur high costs, which has motivated researchers to consider learning augmented algorithms that can make use of faulty or incomplete predictions. Inspired by two matching problems in computational sustainability where data are abundant, we consider the learning augmented min weight matching problem where some nodes are revealed online while others are known a priori, e.g., by being pre- dicted by machine learning. We develop an algorithm that is able to make use of this extra information and provably im- proves upon pessimistic online algorithms. We evaluate our algorithm on two settings from computational sustainability – the coordination of opportunistic citizen scientists for inva- sive species management and the matching between taxis and riders under uncertain trip duration predictions. In both cases, we perform extensive experiments on real-world datasets and find that our method outperforms baselines, showing how learning augmented algorithms can reliably improve solu- tions for problems in computational sustainability

New York↗
USGS1997

Temporal gravity and height changes of the Yellowstone caldera, 1977 - 1994

This paper describes the longest record of gravity measurements in the area of the Yellowstone caldera, Wyoming. The temporal gravity changes, at the ±12 µGal (10 −8 ms −2 ) precision level, are compared with changes in heights from leveling and GPS. The gravity field decreased across the caldera from 1977 to 1983 during the uplift and attained a maximum decrease of up to −60± 2 µGal along the Caldera axis. The gravity field then reversed polarity to increasing values, of up to 60±12 µGal between 1986 and 1993. The ratio between height and gravity changes varied during the entire time, but converged over the latter period following the free-air gravity gradient. General ground deformation deduced from leveling showed caldera-wide uplift of ∼15 mm/a during the period of gravity decrease, then from leveling and GPS, subsidence of ∼25 mm/a during the gravity increase.

Geophysical Research Letters↗
USGS2016

Migratory bird habitat in relation to tile drainage and poorly drained hydrologic soil groups

The Prairie Pothole Region (PPR) is home to more than 50% of the migratory waterfowl in North America. Although the PPR provides an abundance of temporary and permanent wetlands for nesting and feeding, increases in commodity prices and agricultural drainage practices have led to a trend of wetland drainage. The Northern Shoveler is a migratory dabbling duck species that uses wetland habitats and cultivated croplands in the PPR. Richland County in North Dakota and Roberts County in South Dakota have an abundance of wetlands and croplands and were chosen as the study areas for this research to assess the wetland size and cultivated cropland in relation to hydrologic soil groups for the Northern Shoveler habitat. This study used geographic information system data to analyze Northern Shoveler habitats in association with Natural Resource Conservation Service soil data. Habitats, which are spatially associated with certain hydrologic soil groups, may be at risk of artificial drainage installations because of their proximity to cultivated croplands and soil lacking in natural drainage that may become wet or inundated. Findings indicate that most wetlands that are part of Northern Shoveler habitats were within or adjacent to cultivated croplands. The results also revealed soil hydrologic groups with high runoff potential and low water transmission rates account for most of the soil within the Northern Shoveler‘s wetland and cropland habitats. Habitats near agriculture with high runoff potential are likely to be drained and this has the potential of reducing Northern Shoveler habitat.

Conference Paper↗
USGS2018

Streams do work: Measuring the work of low-order streams on the landscape using point clouds

The mutable nature of low-order streams makes regular updating of surface water maps necessary for accurate representation. Low-order streams make up roughly half the streams in the conterminous United States by length, and small inaccuracies in stream head location can result in significant error in stream reach, order, and density. Reliable maps of stream features are vital for hydrologic modeling, ecosystem research, and boundary monitoring. High resolution digital elevation models derived from lidar data have shown promise in low order stream modeling yet forested high relief landscapes and low relief agricultural areas remain challenging. Here we present early results from research analyzing lidar point clouds to identify features and patterns that may be used in low-order stream identification and classification in challenging geographic conditions. This work has identified characteristics derived from point clouds that correlate with the presence of streams and stream heads and show promise for mapping small streams. In low topographic relief agricultural areas, cross sections collected at regular intervals along drainage channels extracted as 3D lines show a significant jump in value and variance of profile curvature standard deviation at stream heads. In high relief areas, observations show potential for stream mapping by identifying trends in riparian zone structure. Lidar return point density from riparian vegetation under 30 feet tall dips in the vicinity of intermittent stream heads. Also seen is an increase in point density above 60 feet downstream of stream heads. The trends found here likely reflect a change in vegetation structure relative to the presence of streams.

Illinois, Iowa, North Carolina↗
USGS2007

A method to estimate groundwater depletion from confining layers

Although depletion of storage in low‐permeability confining layers is the source of much of the groundwater produced from many confined aquifer systems, it is all too frequently overlooked or ignored. This makes effective management of groundwater resources difficult by masking how much water has been derived from storage and, in some cases, the total amount of water that has been extracted from an aquifer system. Analyzing confining layer storage is viewed as troublesome because of the additional computational burden and because the hydraulic properties of confining layers are poorly known. In this paper we propose a simplified method for computing estimates of confining layer depletion, as well as procedures for approximating confining layer hydraulic conductivity ( K ) and specific storage ( S s ) using geologic information. The latter makes the technique useful in developing countries and other settings where minimal data are available or when scoping calculations are needed. As such, our approach may be helpful for estimating the global transfer of groundwater to surface water. A test of the method on a synthetic system suggests that the computational errors will generally be small. Larger errors will probably result from inaccuracy in confining layer property estimates, but these may be no greater than errors in more sophisticated analyses. The technique is demonstrated by application to two aquifer systems: the Dakota artesian aquifer system in South Dakota and the coastal plain aquifer system in Virginia. In both cases, depletion from confining layers was substantially larger than depletion from the aquifers.

Water Resources Research↗