Planning Urban Growth in Ground Water Recharge Areas: Central Valley, Costa Rica

Nitrate levels in the ground water of the Central Valley of Costa Rica have increased in relation to the past. Previous studies determined that the unseweved sanitation systems in the recharge areas are the main source of nitrogen. Calculations are made in this study to estimate the maximum population density allowable without improved sewage systems in order to keep the nitrogen levels in ground water below the World Health Organization (WHO) criteria. Results were achieved employing a mass balance that involved the concentration and rate of domestic effluents and the flow rate in the aquifer, as well as an estimation of the effects caused by the agricultural activity. It was concluded that, in general terms, the population density must not exceed 45 inhabitants per hectare. Otherwise, sewage systems and treatment plants are necessary. These conclusions provide a basis for urban growth planning, which will protect ground water quality. The method used in this case should apply to similar situations.     

Importance of river water recharge to the San Joaquin Valley groundwater system

Groundwater is not a sustainable resource, unless abstraction is balanced by recharge. Identifying the sources of recharge in a groundwater basin is critical for sustainable groundwater management. We studied the importance of river water recharge to groundwater in the south‐eastern San Joaquin Valley (24,000 km², population 4 million). We combined dissolved noble gas concentrations, stable isotopes, tritium, and carbon‐14 analyses to analyse the sources, mechanisms, and timescales of groundwater recharge. Area‐representative groundwater sampling and numerical model input data enabled a stable isotope mass balance and quantitative estimates of river and local recharge. River recharge, identified by a lighter stable isotope signature, represents 47 ± 4% of modern groundwater in the San Joaquin Valley (recharged after 1950) but only 26 ± 4% of premodern groundwater (recharged before 1950). This implies that the importance of river water recharge in the San Joaquin valley has nearly doubled and is likely the result of a 40% increase in total recharge, caused by river water irrigation return flows and increased stream depletion and river recharge due to groundwater pumping. Compared with the large and long‐duration capacity for water storage in the subsurface, storage of water in rivers is limited in time and volume, as evidenced by cold river recharge temperatures resulting from fast infiltration and recharge. Groundwater banking of seasonal surface water flows and expansion of managed aquifer recharge practices therefore appear to be a natural and promising method for increasing the resilience of the San Joaquin Valley water supply system.     

Modelling lake stage and water balance of Lake Tana,Ethiopia

The level of Lake Tana, Ethiopia, fluctuates annually and seasonally following the patterns of changes in precipitation. In this study, a mass balance approach is used to estimate the hydrological balance of the lake. Water influx from four major rivers, subsurface inflow from the floodplains, precipitation, outflow from the lake constituting river discharge and evapotranspiration from the lake are analysed on monthly and annual bases. Spatial interpolation of precipitation using rain gauge data was conducted using kriging. Outflow from the lake was identified as the evaporation from the lake's surface as well as discharge at the outlet where the Blue Nile commences. Groundwater inflow is estimated using MODular three‐dimensional finite‐difference ground‐water FLOW model software that showed an aligned flow pattern to the river channels. The groundwater outflow is considered negligible based on the secondary sources that confirmed the absence of lake water geochemical mixing outside of the basin. Evaporation is estimated using Penman's, Meyer's and Thornwaite's methods to compare the mass balance and energy balance approaches. Meteorological data, satellite images and temperature perturbation simulations from Global Historical Climate Network of National Oceanographic and Atmospheric Administration are employed for estimation of evaporation input parameters. The difference of the inflow and outflow was taken as storage in depth and compared with the measured water level fluctuations. The study has shown that the monthly and annually calculated lake level replicates the observed values with root mean square error value of 0·17 and 0·15 m, respectively. Copyright © 2009 John Wiley & Sons, Ltd.     

Coupled flow and salinity transport modelling in semi-arid environments: The Shashe River Valley,Botswana

Numerical groundwater modelling is used as the base for sound aquifer system analysis and water resources assessment. In many cases, particularly in semi-arid and arid regions, groundwater flow is intricately linked to salinity transport. A case in point is the Shashe River Valley in Botswana. A freshwater aquifer located around an ephemeral stream is depleted by the combined effect of transpiration and pumping. Quantitative system analysis reveals that the amount of water taken by transpiration is far more than the quantities pumped for water supply. Furthermore, the salinity distribution in and around Shashe River Valley as well as its temporal dynamics can be satisfactorily reproduced if the transpiration is modelled as a function of groundwater salinity. The location and dynamics of the saltwater–freshwater interface are highly sensitive to the parameterization of evaporative and transpirative salt enrichment. An existing numerical code for coupled flow/transport simulations (SEAWAT) was adapted to this situation. Model results were checked against a large set of field data including water levels, water chemistry, isotope data and ground and airborne geophysical data. The resulting groundwater model was able to reproduce the long-term development of the freshwater lens located in Shashe River Valley as well as the decline in piezometric heads observed over the last decade. Furthermore, the old age of the saline water surrounding the central freshwater lens could be explained.     

Long‐Term Groundwater Depletion in the United States

The volume of groundwater stored in the subsurface in the United States decreased by almost 1000 km³ during 1900–2008. The aquifer systems with the three largest volumes of storage depletion include the High Plains aquifer, the Mississippi Embayment section of the Gulf Coastal Plain aquifer system, and the Central Valley of California. Depletion rates accelerated during 1945–1960, averaging 13.6 km³/year during the last half of the century, and after 2000 increased again to about 24 km³/year. Depletion intensity is a new parameter, introduced here, to provide a more consistent basis for comparing storage depletion problems among various aquifers by factoring in time and areal extent of the aquifer. During 2001–2008, the Central Valley of California had the largest depletion intensity. Groundwater depletion in the United States can explain 1.4% of observed sea‐level rise during the 108‐year study period and 2.1% during 2001–2008. Groundwater depletion must be confronted on local and regional scales to help reduce demand (primarily in irrigated agriculture) and/or increase supply.     

Significance of stemflow in groundwater recharge. 1: Evaluation of the stemflow contribution to recharge using a mass balance approach

Stemflow was evaluated in a water balance and its contribution to groundwater recharge determined. Gross precipitation, throughfall and stemflow were measured for one year in a pine forest (Tsukuba, Japan) to determine each component of the water balance in the forest. Groundwater recharge rates by stemflow and throughfall were calculated from a mass balance method using chloride in subsurface waters. The stemflow in the water balance was relatively small when estimated as a value per canopy projected area of the tree in the forest. However, the results for the mass balance of chloride in subsurface waters indicated that it was impossible to disregard the stemflow in determining groundwater recharge. Although the ratio of stemflow to the net precipitation was small in the water balance, the effect of stemflow on groundwater recharge was relatively large.     

Tracer hydrology of the data-scarce and heterogeneous Central American Isthmus

Numerous socio-economic activities depend on the seasonal rainfall and groundwater recharge cycle across the Central American Isthmus. Population growth and unregulated land use changes resulted in extensive surface water pollution and a large dependency on groundwater resources. This work combines stable isotope variations in rainfall, surface water, and groundwater of Costa Rica, Nicaragua, El Salvador, and Honduras to develop a regionalized rainfall isoscape, isotopic lapse rates, spatial–temporal isotopic variations, and air mass back trajectories determining potential mean recharge elevations, moisture circulation patterns, and surface water–groundwater interactions. Intra-seasonal rainfall modes resulted in two isotopically depleted incursions (W-shaped isotopic pattern) during the wet season and two enriched pulses during the mid-summer drought and the months of the strongest trade winds. Notable isotopic sub-cloud fractionation and near-surface secondary evaporation were identified as common denominators within the Central American Dry Corridor. Groundwater and surface water isotope ratios depicted the strong orographic separation into the Caribbean and Pacific domains, mainly induced by the governing moisture transport from the Caribbean Sea, complex rainfall producing systems across the N-S mountain range, and the subsequent mixing with local evapotranspiration, and, to a lesser degree, the eastern Pacific Ocean fluxes. Groundwater recharge was characterized by (a) depleted recharge in highland areas (72.3%), (b) rapid recharge via preferential flow paths (13.1%), and enriched recharge due to near-surface secondary fractionation (14.6%). Median recharge elevation ranged from 1,104 to 1,979 m a.s.l. These results are intended to enhance forest conservation practices, inform water protection regulations, and facilitate water security and sustainability planning in the Central American Isthmus.     

Tectonic influences on ground water quality: insight from complementary methods

A study using multiple techniques provided insight into tectonic influences on ground water systems; the results can help to understand ground water systems in the tectonically active western United States and other parts of the world. Ground water in the San Bernardino Valley (Arizona, United States and Sonora, Mexico) is the main source of water for domestic use, cattle ranching (the primary industry), and the preservation of threatened and endangered species. To improve the understanding of ground water occurrence, movement, and sustainability, an investigation was conducted using a number of complementary methods, including major ion geochemistry, isotope hydrology, analysis of gases dissolved in ground water, aquifer testing, geophysics, and an examination of surface and subsurface geology. By combining information from multiple lines of investigation, a more complete picture of the basin hydrogeology was assembled than would have been possible using fewer methods. The results show that the hydrogeology of the San Bernardino Valley is markedly different than that of its four neighboring basins in the United States. The differences include water quality, chemical evolution, storage, and residence time. The differences result from the locally unique geology of the San Bernardino Valley, which is due to the presence of a magmatically active accommodation zone (a zone separating two regions of normal faults with opposite dips). The geological differences and the resultant hydrological differences between the San Bernardino Valley and its neighboring basins may serve as a model for the distinctive nature of chemical evolution of ground water in other basins with locally distinct tectonic histories.     

Calculating ground water transit time of horizontal flow through leaky aquifers

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.     

Modeled Ground Water Age Distributions

The age of ground water in any given sample is a distributed quantity representing distributed provenance (in space and time) of the water. Conventional analysis of tracers such as unstable isotopes or anthropogenic chemical species gives discrete or binary measures of the presence of water of a given age. Modeled ground water age distributions provide a continuous measure of contributions from different recharge sources to aquifers. A numerical solution of the ground water age equation of Ginn (1999) was tested both on a hypothetical simplified one‐dimensional flow system and under real world conditions. Results from these simulations yield the first continuous distributions of ground water age using this model. Complete age distributions as a function of one and two space dimensions were obtained from both numerical experiments. Simulations in the test problem produced mean ages that were consistent with the expected value at the end of the model domain for all dispersivity values tested, although the mean ages for the two highest dispersivity values deviated slightly from the expected value. Mean ages in the dispersionless case also were consistent with the expected mean ages throughout the physical model domain. Simulations under real world conditions for three dispersivity values resulted in decreasing mean age with increasing dispersivity. This likely is a consequence of an edge effect. However, simulations for all three dispersivity values tested were mass balanced and stable demonstrating that the solution of the ground water age equation can provide estimates of water mass density distributions over age under real world conditions.     

Design and Implementation of Evapotranspiration Measuring Equipment for Owens Valley, California

As part of a plant survivability and ground water study in Owens Valley, California, semipermanent installations are used to measure continuous range-land evapotranspiration in the valley's phreatophyte community. A proposed mobile installation also has been designed. The semipermanent micrometeoro-logical station collects continuous data for solution of the Bowen ratio/energy budget equation and the Penman combination equation. Three sites were chosen for this type of installation to provide a representative sampling of Owens Valley. The proposed mobile aerodynamic installation should be capable of calculating evapotranspiration by the eddy correlation method. This instrumentation will be used throughout the valley for short periods of time (up to five days). Many problems with equipment operation, calibration and design have been identified and resolved by means of improved calibration techniques, systematic error-removal techniques, reduced cycle times, modified equipment design and proper observer training. The collected evapotranspiration data will be instrumental in developing a one-dimensional evapotranspiration flux algorithm for a model of valleywide ground water flow.     

Influence of fracture anisotropy on ground water ages and chemistry,Valley and Ridge province,Pennsylvania

Model ground water ages based on chlorofluorocarbons (CFCs) and tritium/helium-3 (³H/³He) data were obtained from two arrays of nested piezometers located on the north limb of an anticline in fractured sedimentary rocks in the Valley and Ridge geologic province of Pennsylvania. The fracture geometry of the gently east plunging fold is very regular and consists predominately of south dipping to subhorizontal to north dipping bedding-plane parting and east striking, steeply dipping axial-plane spaced cleavage. In the area of the piezometer arrays, which trend north-south on the north limb of the fold, north dipping bedding-plane parting is a more dominant fracture set than is steeply south dipping axial-plane cleavage. The dating of ground water from the piezometer arrays reveals that ground water traveling along paths parallel to the dip direction of bedding-plane parting has younger ³H/³He and CFC model ages, or a greater component of young water, than does ground water traveling along paths opposite to the dip direction. In predominantly unmixed samples there is a strong positive correlation between age of the young fraction of water and dissolved sodium concentration. The travel times inferred from the model ages are significantly longer than those previously calculated by a ground water flow model, which assumed isotropically fractured layers parallel to topography. A revised model factors in the directional anisotropy to produce longer travel times. Ground water travel times in the watershed therefore appear to be more influenced by anisotropic fracture geometry than previously realized. This could have significant implications for ground water models in other areas underlain by similarly tilted or folded sedimentary rock, such as elsewhere in the Valley and Ridge or the early Mesozoic basins.     

Simulating a lake as a high-conductivity variably saturated porous medium

One approach for simulating ground water–lake interactions is to incorporate the lake into the ground water solution domain as a high-conductivity region. Previous studies have developed this approach using fully saturated models. This study extends this approach to variably saturated models, so that ground water–lake interactions may be more easily simulated with commonly used or public domain variably saturated codes that do not explicitly support coupled lake–water balance modeling. General guidelines are developed for the choices of saturated hydraulic conductivity and moisture retention and relative permeability curves for the lake region. When applied to an example ground water–lake system, model results are very similar to those from a model in which the lake is represented as a specified head boundary continuously updated by a lake mass balance. The high-conductivity region approach is most suitable for relatively simple geometries and lakes with slower and smaller fluctuations when the overall flow pattern and system fluxes, rather than the detailed flow pattern around the intersection of the lake and land surfaces, are of interest.     

History of the Sole Source Aquifer Program: A Community-Based Approach for Protecting Aquifers Used for Drinking Water Supply

The Sole Source Aquifer Program has helped prevent contamination of many community drinking water supplies. If an aquifer supplies the sole or principal source of a community's drinking water, a local ground water user may petition the Environmental Protection Agency (EPA) under the Safe Drinking Water Act for its designation and protection as a "sole source aquifer." Since 1974, residents and officials of 65 communities and multi-community areas have petitioned and received assistance from the EPA to prevent contamination of their local ground water source of drinking water. This designation means that EPA may review federal financially assisted projects to determine if they would contaminate the aquifer and cause a public health hazard. If they could cause contamination, EPA can request that the project be modified or stopped. The significance of this program in terms of population served and funds affected has been substantial, indicating the Sole Source Aquifer Program has been an important local tool for protecting ground water used as a source of drinking water. Information is given on three different examples of sole source aquifer designations protected under this program: the New Jersey Coastal Plain Aquifer System, the Great Miami River Buried Valley Aquifer System (Ohio), and the Eastern Snake River Plain Aquifer (Idaho), serving populations of 543,000, 921,000, and 275,000, respectively. In all three examples, preventing ground water contamination through the Sole Source Aquifer Program has protected the community drinking water supply.     

Impact of precipitation seasonality changes on isotopic signals in polar ice cores: a multi-model analysis

For Central Greenland, water isotope analysis indicates a temperature difference of about 10°C since the Last Glacial Maximum (LGM). However, borehole thermometry and gas diffusion thermometry indicate that LGM surface temperatures were about 20°C colder than today. Two general circulation model studies have shown that changes in the seasonal precipitation timing in Central Greenland might have caused a warm bias in the LGM water isotope proxy temperatures, and that this bias could explain the difference in the estimated paleotemperatures. Here we present an analysis of a number of atmospheric general circulation model simulations mostly done within the framework of the Paleoclimate Modeling Intercomparison Project. The models suggest that the seasonal cycle of precipitation and surface mass balance over Central Greenland at the LGM might have been very different from today. This supports the idea that the accuracy of the water isotope thermometry at the LGM in Greenland might be compromised as a result of a modified surface mass balance seasonality. However, the models disagree on the amplitude and sign of the bias. For Central East Antarctica, a strong seasonality effect on the LGM isotopic signal is not simulated by any of the analyzed models. For the mid-Holocene (6 kyr BP) the models suggest relatively weak isotope paleothermometry biases linked to changes in the surface mass balance seasonality over both ice sheets.     

Seasonal and interannual behaviour of groundwater catchment boundaries in a Chalk aquifer

Groundwater catchment boundaries and their associated groundwater catchment areas are typically assumed to be fixed on a seasonal basis. We investigated whether this was true for a highly permeable carbonate aquifer in England, the Berkshire and Marlborough Downs Chalk aquifer, using both borehole hydrograph data and a physics‐based distributed regional groundwater model. Borehole hydrograph data time series were used to construct a monthly interpolated water table surface, from which was then derived a monthly groundwater catchment boundary. Results from field data showed that the mean annual variation in groundwater catchment area was about 20% of the mean groundwater catchment area, but interannual variation can be very large, with the largest estimated catchment size being approximately 80% greater than the smallest. The flow in the river was also dependent on the groundwater catchment area. Model results corroborated those based on field data. These findings have significant implications for issues such as definition of source protection zones, recharge estimates based on water balance calculations and integrated conceptual modelling of surface water and groundwater systems. Copyright © 2015 John Wiley & Sons, Ltd.     

THE STUDY OF THE GENESIS OF SURFACE WATERS BY OBSERVING THEIR ELECTRICAL RESISTIVITIES*

Specific electrical resistivity of natural waters contains information on their genesis. The authors propose to conduct mass and regime observations of this parameter in river and stream beds. The electrical resistivities in streams flowing from under a glacier reveal details formed at the same time as the glacier. Observations in the beds of big rivers show a gradual increase in water salinity overlain by reductions by inflowing glacial waters. The diurnal and annual trend of changes in the electrical conductivity of water associated with the change in the balance of glacial and ground waters has been established near to glaciers. Resistivity observations help to locate discharge sites of sub-permafrost waters, for water.     

Energy considerations in groundwater‐ridging mechanism of streamflow generation

Groundwater ridging is the rapid rise of a shallow water table during a rainfall event, in an environment where, in the pre‐event period, the capillary fringe extends to the ground surface. Groundwater ridging is widely cited to account for the observed significant appearance of pre‐event water in a stream stormflow hydrograph. Various hypotheses have been advanced to explain the groundwater‐ridging mechanism; and most recently, from a field study site in South Africa, an energy hypothesis was proposed, which explains that groundwater‐ridging water‐table rise is a result of rapid introduction and transmission of additional pressure head into the capillary fringe from an intense rainfall at the ground surface. However, there is a need for further analysis and evidence from other field study sites to confirm and support this newly proposed energy hypothesis. The objectives of this paper are, therefore, as follows: to review previous observations on groundwater ridging, from other study sites, in order to deduce evidence of the newly proposed energy hypothesis; to present and evaluate a one‐dimensional diffusion mathematical model that can simulate groundwater‐ridging water‐table rise, based on the newly proposed energy hypothesis; and to evaluate the importance of a capillary fringe in streamflow generation. Analysis of previous observations from other study sites generally indicated that the rate of groundwater‐ridging water‐table rise is directly related to the rainfall intensity, hence confirming and agreeing with the newly proposed energy hypothesis. Additionally, theoretical results by the mathematical model agreed fairly well with the field results observed under natural rainfall, confirming that the rapidly rainfall‐induced energy is diffusively transmitted downwards through pore water, elevating the pressure head at every depth. The results in this study also support the concept of a three‐end‐member stream stormflow hydrograph and contribute to the explanation of how catchments can store water for long periods but then release it rapidly during storm events. Copyright © 2015 John Wiley & Sons, Ltd.     

Deduction of groundwater flow regime in a basaltic aquifer using geochemical and isotopic data: The Golan Heights, Israel case study

Groundwater flow-paths through shallow-perch and deep-regional basaltic aquifers at the Golan Heights, Israel, are reconstructed by using groundwater chemical and isotopic compositions. Groundwater chemical composition, which changes gradually along flow-paths due to mineral dissolution and water–rock interaction, is used to distinguish between shallow-perched and deep-regional aquifers. Groundwater replenishment areas of several springs are identified based on the regional depletion in rainwater δ¹⁸O values as a function of elevation (−0.25‰ per 100 m). Tritium concentrations assist in distinguishing between pre-bomb and post-bomb recharged rainwater.

It was found that waters emerging through the larger springs are lower in δ¹⁸O than surrounding meteoric water and poor in tritium; thus, they are inferred to originate in high-elevation regions up to 20 km away from their discharge points and at least several decades ago. These results verify the numerically simulated groundwater flow field proposed in a previous study, which considered the geological configuration, water mass balance and hydraulic head spatial distribution.     

Long‐term final void salinity prediction for a post‐mining landscape in the Hunter Valley,New South Wales,Australia

Opencast mining alters surface and subsurface hydrology of a landscape both during and post‐mining. At mine closure, following opencast mining in mines with low overburden to coal ratios, a void is left in the final landform. This final void is the location of the active mine pit at closure. Voids are generally not infilled within the mines' lifetime, because of the prohibitive cost of earthwork operations, and they become post‐mining water bodies or pit lakes. Water quality is a significant issue for pit lakes. Groundwater within coal seams and associated rocks can be saline, depending on the nature of the strata and groundwater circulation patterns. This groundwater may be preferentially drawn to and collected in the final void. Surface runoff to the void will not only collect salts from rainfall and atmospheric fallout, but also from the ground surface and the weathering of fresh rock. As the void water level rises, its evaporative surface area increases, concentrating salts that are held in solution. This paper presents a study of the long term, water quality trends in a post‐mining final void in the Hunter Valley, New South Wales, Australia. This process is complex and occurs long term, and modelling offers the only method of evaluating water quality. Using available geochemical, climate and hydrogeological data as inputs into a mass‐balance model, water quality in the final void was found to increase rapidly in salinity through time (2452 to 8909 mg l⁻¹ over 500 years) as evaporation concentrates the salt in the void and regional groundwater containing high loads of salt continues to flow into the void. Copyright © 2004 John Wiley & Sons, Ltd.