USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Hydrology”

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 235 records · Page 13Linked to original sources
USGS2024

Diel temperature signals track seasonal shifts in localized groundwater contributions to headwater streamflow generation at network scale

Groundwater contributions to streamflow sustain aquatic ecosystem resilience; streams without significant groundwater inputs often have well-coupled air and water temperatures that degrade cold-water habitat during warm low flow periods. Widespread uncertainty in stream-groundwater connectivity across space and time has created disparate predictions of energy and nutrient fluxes across headwater networks, hindering predictions of cold-water habitat resilience under climate change scenarios. Recently, annual paired air and water temperature signals have been harnessed to indicate stream water thermal sensitivity and the dominance of deep versus shallow groundwater influence, although the utility of diel air–water temperature signal metrics for hydrologic inference has remained unexplored. Here we analyzed two consecutive years of locally paired, air–water temperature data from 47 headwater stream sites in the Catskill Mountains, New York, USA, and discovered characteristic seasonal patterns in diel temperature signal sinusoid metrics (amplitude ratio, phase lag, and mean ratio) driven by shifts in streamflow generation mechanisms and stream network position. Hydrologic interpretations of observed patterns were supported by stream heat budget model scenarios and additional analysis of paired air–water temperature data from two streams in Shenandoah National Park, Virginia, USA, with well characterized stream-groundwater connectivity. We found that within smaller tributaries, streamflow generation transitions from runoff to groundwater dominance were driven by hillslope drying during seasonal periods of lower precipitation. This was evidenced by significant correlations (p < 0.01) between daily water:air temperature signal amplitudes (non-linear decreases of ∼ 50 %) and derived base-flow index at 22 of the 28 sites, indicating enhanced local groundwater influence on streamflow promotes decoupling of diel air–water temperature signals. Additionally, ratios between daily water:air temperature signal means were lower in tributaries (∼0.68) when compared to main-stem (∼0.8) sites, increasing linearly throughout the observational period. In conceptual stream heat budget models, groundwater inflow had minimal effects on daily phase lags (∼0.2 hr), but increases in fractional groundwater discharge (0–50 %) depressed daily amplitude (∼20 % to 50 %) and mean ratios (∼15 %), supporting the sensitivity of daily metrics to interpreted changes in seasonal groundwater contributions to streamflow. During observational periods (i.e., April through October 2021 and 2022), significant differences (p < 0.01) between tributary and main-stem air–water metrics occurred when base-flow contributions were highest (∼0.93 vs. ∼ 0.68), as sites lower in the network had daily temperature metrics dominated by stream channel thermal inertia, rather than local groundwater connectivity, showing enhanced air–water diel signal coupling during warmer, drier periods. Divergent air temperature coupling across the network was interpreted as being driven by distance from local groundwater source zones, additional lateral groundwater inflows do not contribute a meaningful fraction to channel discharge lower in the network. Given the growing footprint of stream temperature observations, diel air–water temperature signals can provide distributed metrics sensitive to upstream groundwater discharge. Consequently, these metrics can support ongoing efforts by resource managers and researchers seeking to forecast the resilience of cold-water habitat to climate warming and changing precipitation regimes in mountain headwater streams.

New York, Virginia↗
USGS2008

Temporal downscaling of decadal sediment load estimates to a daily interval for use in hindcast simulations

In this study we used hydrologic proxies to develop a daily sediment load time-series, which agrees with decadal sediment load estimates, when integrated. Hindcast simulations of bathymetric change in estuaries require daily sediment loads from major tributary rivers, to capture the episodic delivery of sediment during multi-day freshwater flow pulses. Two independent decadal sediment load estimates are available for the Sacramento/San Joaquin River Delta, California prior to 1959, but they must be downscaled to a daily interval for use in hindcast models. Daily flow and sediment load data to the Delta are available after 1930 and 1959, respectively, but bathymetric change simulations for San Francisco Bay prior to this require a method to generate daily sediment load estimates into the Delta. We used two historical proxies, monthly rainfall and unimpaired flow magnitudes, to generate monthly unimpaired flows to the Sacramento/San Joaquin Delta for the 1851-1929 period. This step generated the shape of the monthly hydrograph. These historical monthly flows were compared to unimpaired monthly flows from the modern era (1967-1987), and a least-squares metric selected a modern water year analogue for each historical water year. The daily hydrograph for the modern analogue was then assigned to the historical year and scaled to match the flow volume estimated by dendrochronology methods, providing the correct total flow for the year. We applied a sediment rating curve to this time-series of daily flows, to generate daily sediment loads for 1851-1958. The rating curve was calibrated with the two independent decadal sediment load estimates, over two distinct periods. This novel technique retained the timing and magnitude of freshwater flows and sediment loads, without damping variability or net sediment loads to San Francisco Bay. The time-series represents the hydraulic mining period with sustained periods of increased sediment loads, and a dramatic decrease after 1910, corresponding to a reduction in available mining debris. The analogue selection procedure also permits exploration of the morphological hydrograph concept, where a limited set of hydrographs is used to simulate the same bathymetric change as the actual set of hydrographs. The final daily sediment load time-series and morphological hydrograph concept will be applied as landward boundary conditions for hindcasting simulations of bathymetric change in San Francisco Bay.

California↗
USGS2013

Fens as whole-ecosystem gauges of groundwater recharge under climate change

Currently, little is known about the impact of climate change on groundwater recharge in the Sierra Nevada and southern Cascade Range of California or other mountainous regions of the world. The purpose of this study was to determine whether small alpine peat lands called fens can be used as whole-ecosystem gauges of groundwater recharge through time. Fens are sustained by groundwater discharge and are highly sensitive to changes in groundwater flow due to hydrologic disturbance including climate change. Seven fens in the Sierra Nevada and southern Cascade Range were studied over a 50-80 year period using historic aerial photography. In each aerial photograph, fen areas were identified as open lawn and partially treed areas that exhibited (1) dark brownish-green coloring or various shades of gray and black in black and white imagery and (2) mottling of colors and clustering of vegetation, which signified a distinct moss canopy with overlying clumped sedge vegetation. In addition to the aerial photography study, a climate analysis for the study sites was carried out using both measured data (U.S. Department of Agriculture Natural Resources Conservation Service SNOwpack TELemetry system) and modeled data (a downscaled version of the Parameter-elevation Regressions on Independent Slopes Model) for the period from 1951 to 2010. Over the study period, the five fens in the Sierra Nevada were found to be decreasing between 10% and 16% in delineated area. The climate analysis revealed significant increases through time in annual mean minimum temperature (T min ) between 1951-1980 and 1981-2010. In addition, April 1 snow water equivalent and snowpack longevity also decreased between 1951-1980 and 1981-2010. For the fens in the Cascade Range, there were no discernible changes in delineated area. At these sites, increases in T min occurred only within the past 20-25 years and decreases in snowpack longevity were more subtle. A conceptual model is presented, which illustrates that basic differences in hydrogeology of the Sierra Nevada vs. the Cascade Range may control the threshold at which changes in delineated fen areas are discernible. Overall, the results from this study show that fens in the Sierra Nevada have strong potential as whole ecosystem gauges for determining long-term changes in groundwater recharge under climate change. Due to either more moderate climate change and/or hydrogeological differences, fens in the southern Cascade Range currently do not appear to have the same utility. A greater sample size of fens in the Sierra Nevada is needed to confirm the general applicability of this method. In addition, future work needs to focus on integrating fen monitoring with geochemical and/or isotopic process-level studies in order to quantify changes in groundwater recharge identified using this new approach.

Journal of Hydrology↗
USGS1984

Use of a digital model to evaluate hydrogeologic controls on groundwater flow in a fractured rock aquifer at Niagara Falls, New York, U.S.A.

The Hyde Park landfill is a 15-acre (6.1 ha) chemical waste disposal site located north of Niagara Falls, New York. Underlying the site in descending order are: (1) low-permeability glacial till and lacustrine deposits; (2) a moderately permeable fractured rock aquifer - the Lockport Dolomite; and (3) a low-permeability unit - the Rochester Shale. The site is bounded on three sides by groundwater drains; the Niagara River gorge, the Niagara Power Project canal, and the Niagara Power Project buried conduits. The mechanism by which groundwater moves through fractured rocks underlying a hazardous waste site was investigated using a digital simulation approach. Three hypotheses were tested related to flow in the fractured rocks underlying Hyde Park landfill. For this purpose we used a Galerkin finite-element approximation to solve a saturated-unsaturated flow equation. A primary focus was to investigate anisotropy in the Lockport Dolomite, that is the effectiveness of horizontal (bedding) joints vs. vertical joints as water-transmitting openings. Three hydrogeologic scenarios were set up — each with prescribed limits on the hydrologic parameters. Scenario 1 specified strongly anisotropic conditions in the Lockport Dolomite (horizontal hydraulic conductivity along bedding joints exceeds vertical conductivity by 2–3 orders of magnitude), uniform areal recharge (5 in. yr. −1 or 12.7 cm yr. −1 ) except at the landfill where there is no recharge, and no flow through the base of the Rochester Shale. Scenario 2 also specified strongly anisotropic conditions in the Lockport; however, areal recharge was 6 in. yr. −1 (15.2 cm yr. −1 ) except at the landfill where the recharge was 2 in. yr. −1 (5.1 cm yr. −1 ), and outflow from the Rochester occurred. Scenario 3 specified isotropic conditions (that is, permeability along horizontal and vertical joints is the same in the Lockport Dolomite), recharge rates were the same as in scenario 2 and outflow through Rochester occurred. Scenario 2 provided the closest agreement between the simulated and measured heads while scenario 3 provided the poorest agreement. Among the three scenarios tested, scenario 2 (with strongly anisotropic conditions in the Lockport Dolomite with added recharge through the landfill cap and limited flow through the Rocherster Shale) is considered the most realistic hydrogeologic model. Based on simulation with the hydrogeologic parameters of scenario 2, groundwater flow near the Hyde Park site can be summarized as follows: 1. (1) Specific discharge (Darcy velocity) ranges from ≈0.01 to 0.1 ft. day −1 (0.003 to 0.03 m day −1 ) in the upper unit of the Lockport Dolomite to slightly more than 0.0001 ft. day −1 (0.00003 m day −1 ) in the Rochester Shale. Real velocities are highest in the upper unit of the Lockport, ranging from ≈1 to 5 ft. day −1 (0.3 to 1.5 m day −1 ) if the average effective porosity is assumed to be 0.02. 2. (2) A groundwater divide exists east of the landfill, indicating that all groundwater originating near or flowing beneath the landfill will flow toward and discharge in the gorge. 3. (3) Highest flow velocities (and presumably greatest potential for transporting chemical contaminants) occur in the upper unit of the Lockport and part of the lower unit of the Lockport Dolomite between the landfill and the gorge. The average time required for groundwater to move from the landfill to the discharge points at the gorge along selected flow paths in the Lockport Dolomite is estimated to be 5-6 yr.

New York↗
USGS2014

Analysis of projected water availability with current basin management plan, Pajaro Valley, California

The projection and analysis of the Pajaro Valley Hydrologic Model (PVHM) 34 years into the future using MODFLOW with the Farm Process (MF-FMP) facilitates assessment of potential future water availability. The projection is facilitated by the integrated hydrologic model, MF-FMP that fully couples the simulation of the use and movement of water from precipitation, streamflow, runoff, groundwater flow, and consumption by natural and agricultural vegetation throughout the hydrologic system at all times. MF-FMP allows for more complete analysis of conjunctive-use water-resource systems than previously possible with MODFLOW by combining relevant aspects of the landscape with the groundwater and surface-water components. This analysis is accomplished using distributed cell-by-cell supply-constrained and demand-driven components across the landscape within &ldquo;water-balance subregions&rdquo; (WBS) comprised of one or more model cells that can represent a single farm, a group of farms, watersheds, or other hydrologic or geopolitical entities. Analysis of conjunctive use would be difficult without embedding the fully coupled supply-and-demand into a fully coupled simulation, and are difficult to estimate a priori. The analysis of projected supply and demand for the Pajaro Valley indicate that the current water supply facilities constructed to provide alternative local sources of supplemental water to replace coastal groundwater pumpage, but may not completely eliminate additional overdraft. The simulation of the coastal distribution system (CDS) replicates: 20 miles of conveyance pipeline, managed aquifer recharge and recovery (MARR) system that captures local runoff, and recycled-water treatment facility (RWF) from urban wastewater, along with the use of other blend water supplies, provide partial relief and substitution for coastal pumpage (aka in-lieu recharge). The effects of these Basin Management Plan (BMP) projects were analyzed subject to historical climate variations and assumptions of 2009 urban water demand and land use. Water supplied directly from precipitation, and indirectly from reuse, captured local runoff, and groundwater is necessary but inadequate to satisfy agricultural demand without coastal and regional storage depletion that facilitates seawater intrusion. These facilities reduce potential seawater intrusion by about 45% with groundwater levels in the four regions served by the CDS projected to recover to levels a few feet above sea level. The projected recoveries are not high enough to prevent additional seawater intrusion during dry-year periods or in the deeper aquifers where pumpage is greater. While these facilities could reduce coastal pumpage by about 55% of the historical 2000&ndash;2009 pumpage for these regions, and some of the water is delivered in excess of demand, other coastal regions continue to create demands on coastal pumpage that will need to be replaced to reduce seawater intrusion. In addition, inland urban and agricultural demands continue to sustain water levels below sea level causing regional landward gradients that also drive seawater intrusion. Seawater intrusion is reduced by about 45% but it supplies about 55% of the recovery of groundwater levels in the coastal regions served by the CDS. If economically feasible, water from summer agricultural runoff and tile-drain returnflows could be another potential local source of water that, if captured and reused, could offset the imbalance between supply and demand as well as reducing discharge of agricultural runoff into the National Marine Sanctuary of Monterey Bay. A BMP update (2012) identifies projects and programs that will fund a conservation program and will provide additional, alternative water sources to reduce or replace coastal and inland pumpage, and to replenish the aquifers with managed aquifer recharge in an inland portion of the Pajaro Valley.

California↗
USGS2012

A remote sensing approach for estimating the location and rate of urban irrigation in semi-arid climates

Urban irrigation is an important component of the hydrologic cycle in many areas of the arid and semiarid western United States. This paper describes a new approach that uses readily available datasets to estimate the location and rate of urban irrigation. The approach provides a repeatable methodology at 1/3 km 2 resolution across a large urbanized area (500 km 2 ). For this study, Landsat Thematic Mapper satellite imagery, air photos, climatic records, and a land-use map were used to: (1) identify the fraction of irrigated landscaping in urban areas, and (2) estimate the monthly rate of irrigation being applied to those areas. The area chosen for this study was the San Fernando Valley in Southern California. Identifying irrigated areas involved the use of 29 satellite images, air photos, and a land-use map. The fraction of a pixel that consists of irrigated landscaping (F irr ) was estimated using a linear-mixture model of two land-cover endmembers (selected pixels within a satellite image that represent a targeted land-cover). The two endmembers were impervious and fully-irrigated landscaping. In the San Fernando Valley, we used airport buildings, runways, and pavement to represent the impervious endmember; golf courses and parks were used to represent the fully irrigated endmember. The average F irr using all 29 satellite scenes was 44%. F irr calculated from hand-digitizing using air photos for 13 randomly selected single-family-residential neighborhoods showed similar results (42%). Estimating the rate of irrigation required identification of a third endmember: areas that consisted of urban vegetation but were not irrigated. This "nonirrigated" endmember was used to compute a Normalized Difference Vegetation Index (NDVI) surplus, defined as the difference between the NDVI signals of the irrigated and nonirrigated endmembers. The NDVI signals from irrigated areas remains relatively constant throughout the year, whereas the signal from nonirrigated areas rises and falls seasonally due to precipitation. The areas between airport runways were chosen to represent the nonirrigated endmember. Water-delivery records from 65 spatially-distributed single-family neighborhoods, consisting of nearly 1800 homes, were correlated with the NDVI surplus. The results show a strong exponential correlation ( r 2 = 0.94). In the absence of water-delivery records, which can be difficult to obtain, a surrogate was identified: the landscape evapotranspiration rate (ET L ). ET L was used to scale NDVI surplus (which is dimensionless) to irrigation rates using an exponential scaling function. The monthly irrigation rates calculated from satellite and climatic data compared well with irrigation rates calculated from actual water-delivery data using a paired Wilcoxan signed-rank test ( p = 0.0063). Identification of F irr at the pixel scale, along with identification of the irrigation rate for a fully-irrigated pixel, allows for mapping of urban irrigation over large areas. Maps showing the location and rate of monthly irrigation for the San Fernando study area were computed for January and August 1997.

California↗
USGS1994

Hydrogeologic analysis of the saturated-zone ground-water system, under Yucca Mountain, Nevada

The configuration of the southward-sloping water table under Yucca Mountain is dominated by an abrupt decline of 300 m over a distance of less than 2 km. This northeast-striking zone of large hydraulic gradient (of 0.15 or more) separates an area of moderate gradient (of about 0.015) to the north from an area of very small gradient (0.0001) to the south. The position of the large gradient does not correlate well with any evident geologic feature in the upper 0.5 km of the mountain, but we suggest that buried geologic features are present that can explain all the geohydrologic observations. The three areas of differing hydraulic gradient under Yucca Mountain are parts of hydrogeologic domains that extend more than 70 km to the northeast. On a regional basis, the moderate and very small gradients generally correspond to areas underlain by exceptionally thick Tertiary volcanic sections and a highly transmissive Paleozoic carbonate aquifer, respectively. The regional large gradient and water-table decline are spatially associated with a contact in the Paleozoic rocks between clastic rocks and carbonates. This contact marks a large abrupt drop in the effective base of the hydrologic system because it is the upgradient boundary of the deep carbonate aquifer, which has a thickness of 5 km. An aeromagnetic high follows the regional-scale domain of large gradient under northern Yucca Mountain from outcrops of a magnetite-bearing clastic confining unit to the east, indicating that the regional correlation of the steep water-table decline with the upgradient boundary of the deep carbonate aquifer may extend to Yucca Mountain. Five additional features may be related to an explanation for the large hydraulic gradient: (1) anomalously low heat flow has been measured deep in the volcanic section south of the water-table decline, suggesting underflow of cool water in the deep carbonate aquifer; (2) the lower tuff sequence, of 0.5–1 km in thickness, which underlies most of Yucca Mountain, is largely replaced in the volcanic section by lavas in the area of the large gradient; (3) an analysis of the hydrogeology of the tuff section under Yucca Mountain indicates that transmissivity in the tuffs increases to the south; (4) a northeast-trending gravity low is present immediately south of the water-table decline; (5) units in the lower part of the volcanic section are 50–100% thicker in the area of the gravity low than to the north and south. The abrupt stratigraphic thickening into the area of the gravity low indicates that the low represents a buried graben with its northern bounding fault centered beneath the abrupt water-table decline. These geologic features of the zone of large gradient under Yucca Mountain suggest two possible hydrogeologic models. First, the northern bounding fault of the buried graben may provide a highly permeable pathway (a drain) through the brittle lavas in the lower part of the volcanic section under northern Yucca Mountain. The drain would allow flow from the tuff aquifer north of the decline to be captured by the deep carbonate aquifer, resulting in the heat-flow low, the abrupt water-table decline, and the transition to a very small hydraulic gradient. Alternatively, the northern bounding fault of the buried graben may be the effective northern limit of the tuff aquifer under Yucca Mountain because the permeability in the tuffs north of the fault may have been diminished by hydrothermal alteration. In this second model, the large gradient marks the point where the small southward flow of water through the altered volcanic rocks to the north abruptly drops into the tuff aquifer. In either case, heads in the tuff aquifer in the area of very small gradient may be regulated partly by upward flow from the deep carbonate aquifer. This upward flow under southern Yucca Mountain is indicated by linear thermal highs along fault zones, by ground-water isotopic data suggesting inmixing of waters from the deep carbonate aquifer into the tuff aquifer, and by the upward hydraulic gradient found in a drill hole that penetrates the deep carbonate aquifer under southern Yucca Mountain.

Journal of Hydrology↗
USGS1979

Effects of karst and geologic structure on the circulation of water and permeability in carbonate aquifers

The results of the natural processes caused by solution and leaching of limestone, dolomite, gypsum, salt and other soluble rocks, is known as karst. Development of karst is commonly known as karstification, which may have a pronounced effect on the topography, hydrology and environment, especially where such karst features as sinkholes and vertical solution shafts extend below the land surface and intersect lateral solution passages, cavities, caverns and other karst features in carbonate rocks. Karst features may be divided into two groups: (1) surficial features that do not extend far below the surface; and (2) karst features such as sinkholes that extend below the surface and affect the circulation of water below. The permeability of the most productive carbonate aquifers is due chiefly to enlargement of fractures and other openings by circulation of water. Important controlling factors responsible for the development of karst and permeability in carbonate aquifers include: (1) climate, topography, and presence of soluble rocks; (2) geologic structure; (3) nature of underground circulation; and (4) base level. Another important factor is the condition of the surface of the carbonate rocks at the time they are exposed to meteoric water. A carbonate rock surface, with soil or relatively permeable, less soluble cover, is more favorable for initiation of karstification and solution than bare rocks. Water percolates downward through the cover to the underlying carbonate rocks instead of running off on the surface. Also, the water becomes more corrosive as it percolates through the permeable cover to the underlying carbonate rocks. Where there is no cover or the cover has been removed, the carbonate rocks become case hardened and resistant to erosion. However, in regions underlain not only by carbonate rocks but also by beds of anhydrite, gypsum and salt, such as the Hueco Plateau in southeastern New Mexico, subsurface solution may occur where water without natural acids moves down from bare rock surfaces through cracks to the beds that are more soluble than carbonate rocks. For example, in the area of Carlsbad Caverns in southeastern New Mexico, much of the water responsible for solution that formed the caverns apparently entered the groundwater system through large open fractures and did not form sinkhole topography. East of the Carlsbad Caverns, however, in the Pecos River Valley where the carbonate rocks are overlain by the less soluble Ogallala Formation of Late Tertiary age, solution began along escarpments as the Pecos River and its tributaries cut through the less soluble cover. As these escarpments retreated, sinkholes and other karst features developed. Joints or fractures are essential for initiation of downward percolation of water in compact carbonate rocks such as some Paleozoic limestone in which there is no intergranular permeability. Also joints or fractures and bedding planes may be essential in the initiation of lateral movement of water in the zone of saturation. Where conditions of recharge and discharge are favorable, groundwater may move parallel to the dip. However, the direction of movement of water in most carbonate rocks is not necessarily down dip or parallel to the dip. The general direction of movement of both surface and groundwater may be parallel to the strike in a breached anticline. Faults may restrict the lateral movement of water, especially if water-bearing beds are faulted against relatively impervious beds. Conversely, some fault may serve as avenues through which water may move as, for example, in the Cretaceous Edwards aquifer in the San Antonio area, Texas. Karst aquifers, chiefly carbonate rocks, may be placed in three groups according to water-bearing capacity. Water in aquifers of group 1 occurs chiefly in joints, fractures, and other openings that have not been enlarged by solution. The yield of wells is small. Aquifers in group 2 , with low to intermediate yields, are those in which water occurs in joints and fractures with some cavities and channels enlarged by solution. Aquifers in group 3 are those in which the yield of wells and springs range from intermediate to very large. This group includes five of the most productive aquifers in the U.S.A. The water-bearing beds of all of these productive aquifers, except the Biscayne aquifer in southeastern Florida, contain buried paleokarst in which the permeability has been reactivated and enlarged by the present circulation system.

Journal of Hydrology↗
USGS1963

Applications of geohydrologic concepts in geology

Subsurface water, an active agent in many geologic proceses, must be considered in interpreting geologic phenomena. Principles of the occurrence, distribution, and movement of subsurface waters are well established and readily applicable. In many interpretations in geologic literature, geohydrologic principles have been employed realistically, but in many others these principles have been either ignored or violated. Explanations of genesis of underclays and associated deposits afford some examples wherein principles of movement and activity of vadose and ground water have been ignored and others in which they have been used advantageously. Postulates stating that waters percolate downward from swamp areas do not allow for the usual movement of subsurface water in such environments. The idea that sediments were leached by vadose water after uplift satisfies the geohydrologic requirements. Weathering and solution form porous and permeable zones subjacent to unconformities in dense rocks such as carbonates and granites; this illustrates the geohydrologic and economic significance of unconformities. Examples are Mohawkian carbonate aquifers of northern Illinois and oil-bearing limestones of Mississippian age of eastern Montana. The flushing effects of meteoric water and other hydrodynamic factors active during erosion periods are important elements in the genesis and concentration of brines. Explanation of the origin and occurrence of brines must include consideration of the geohydrologic environments throughout their geologic history.

Journal of Hydrology↗
USGS1965

Natural recharge and localization of fresh ground water in Kuwait

Fresh ground water (200 parts per million total dissolved solids and upwards) occurs in portions of Pleistocene sandstone aquifers beneath basins and wadis in north Kuwait where the mean rainfall is about five inches per year. The fresh water is surrounded and underlain by brackish water (> 4000 ppm TDS). Drilling and testing show that fresh water saturation is restricted to wadis and basin areas; in Rawdatain basin it attains a maximum thickness of about 110 feet and a lateral extent of about seven miles. The fresh ground water represents recharge localized, during infrequent, torrential rain storms, in areas of concentrated runoff where sediments in the vadose zone are moderately permeable and depth to the water table is generally less than a hundred feet. Concentration of runoff appears to be the primary control in the localization of recharge. The fresh water percolates downward to the ground-water reservoir following rare storms, then flows in the direction of hydraulic gradient and gradually becomes brackish. Theoretical delineation of the recharge area and ground-water flow pattern in Rawdatain was confirmed by tritium and C 14 dating of the water. Brackish ground-water conditions prevail from water table downward in areas where rainfall infiltrates essentially where it falls, permeability of sediments in the vadose zone is low, or the water table is several hundred feet below land surface. In these areas, rainfall is retained and lost within the soil zone or becomes mineralized during deep percolation.

Journal of Hydrology↗
USGS1971

Effects of karst features on circulation of water in carbonate rocks in coastal areas

The normal balance between fresh water in coastal aquifers and sea water applies also to carbonate-rock aquifers that have been karstified, but there are local modifications in the balance that need to be considered. Uneven distribution of permeability, expressed by a network of solution channels bounded by relatively impermeabler rock, causes an uneven distribution of head of the environmental water along the seacoast. Where sinkholes and (or) vertical solution shafts below sea level penetrate the aquifer, the fresh ground water may discharge through these karst features if the fresh-water head is greater than that of the salt water. However, under some conditions the salt-water head may exceed that of the fresh water, and the direction of movement is reversed as sea water flows into the aquifer. This sea-water flow into the aquifer occurs (1) where sinkholes, acting as “cased wells,” penetrate less permeable rock before reaching a lateral solution channel and (2) where (or when) the fresh-water head is less than that required to balance the salt water. On Andros Island, Bahamas, the range in tide (as much as 5 feet) from low tide to high tide is sufficient to cause such a reversal locally. During low tide the salt-water head becomes sufficiently low that the ground-water head exceeds that of the sea water, and the ground water flows through the sinkholes to the ocean floor; during high tide sea water flows in the sinkholes. In the Adriatic Sea along the coast of Yugoslavia, apparently the fresh-water head is sufficient to produce perennial springs in some localities, but in other areas, as in the Bay of Kastela near Split, the fresh-water head becomes low enough during some seasons that the flow is reversed and salt water enters the aquifer through the karst features. The development of sinkholes and other karst features near the present coast and extending below sea level occurred generally during a low stand of the Pleistocene sea when the top of the saturated zone stood lower than the bottom of the deepest sinkholes or natural wells. Integrated evaluations of (1) the distribution of permeability in coastal karst regions and (2) the principles relating to the dynamic balance between fresh aquifer water and sea water are leading to better knowledge of methods that may salvage much karst water which is lost to the sea.

Split↗
USGS1974

Water quality effects of seepage from earthen dams

Analyses of surface and seepage waters from selected floodwater retarding structures in west-central Oklahoma, U.S.A. show the salinity of seepage to be several times greater than that of the impounded waters. The increases in concentration of the several chemical components are not proportional. This phenomenon appears to be caused largely by simple solution, which is closely related to the chemical character of the geologic formation that provides the reservoir site and the earth fill of the dam. Concentration of certain chemical components in the seepage water progressively decreases as the structure ages. This is probably a function of depletion, which is related to the amounts and solubility of the parent materials subjected to solution. In contrast, the concentration of certain other components, such as iron and calcium, increases with time. The chemical activities within the accumulating mud on the bottom of the reservoir apparently cause these increases.

Oklahoma↗
USGS1976

Bicarbonate content of groundwater in carbonate rock in eastern North America

In carbonate-rock terrane the most effective solution occurs where soil and vegetative cover facilitate biogenic production and storage of CO 2 in the soil until part of it is carried downward in percolating water. Bicarbonate data for groundwater in eastern North America are examined in the light of these conditions, of the CO 2 content of soil gas, and of the timing of groundwater recharge relative to seasonal changes in soil temperature. There appears to be no well-marked relation between latitude and bicarbonate content of groundwater in this region. Interplay of all the factors listed above, and of lithology and soil types, is evidently such that under optimum conditions the amount of solution of carbonate rock is roughly comparable, per unit of groundwater recharge, over the region from Ohio to Florida (and, perhaps, to Puerto Rico and Yucatan). Relatively low HCO 3 concentrations observed in much of the southern United States are attributed to low production and storage of CO 2 in sandy soil that is poor in organic matter. On the other hand, concentrations observed in Ontario, New York and Michigan are markedly higher than is to be expected from interrelations of the factors considered. These high values, apparently anomalous, are attributed in part to solution of granular glacial drift derived largely from carbonate rock.

Journal of Hydrology↗
USGS1977

The age of groundwater in the Lincolnshire Limestone, England and its relevance to the flow mechanism

Groundwater samples from the Lincolnshire Limestone have been analysed for tritium, radiocarbon, and the stable-isotope ratios 13 C 12 C "> 13 C 12 C , 18 O 16 O "> 18 O 16 O and D/H. The age of the water increases in a downgradient direction below overlying confining deposits and reaches a maximum age greater than 25,000 years within 15 km of the outcrop. The δ 13 C ratio ultimately attains the unusually low negative values of −2‰ in a downgradient direction; this is approaching that of the aquifer matrix which is +2.35‰. The reason is believed to be exchange of carbon between the groundwater and the matrix by a continuous process of precipitation and further solution of calcium carbonate. The δ 18 O and δ D ratios imply that recharge of the aquifer during the late Pleistocene took place in the spring and autumn rather than the winter as at present. The data are interpreted by assuming that movement of water through the saturated zone is a combination of flow in fissures and “piston flow” through the micro-fissures and pores of the rock. The mechanism of water movement in the saturated zone is dominated by relatively rapid flow in fissures, but the fissure flow includes a contribution of much older water from “intergranular” storage which enters the fissures from the rock matrix by pressure differentials in the fissure distribution system and by diffusion. The distribution of water of different ages in the aquifer is closely related to recent groundwater abstraction patterns.

Journal of Hydrology↗
USGS1977

Solution of three-dimensional groundwater flow equations using the strongly implicit procedure

A three-dimensional numerical model has been coded to use the strongly implicit procedure for solving the finite-difference approximations to the ground-water flow equation. The model allows for: (1) the representation of each aquifer and each confining bed by several layers; and (2) the use of an anisotropic hydraulic conductivity at each finite-difference block. The model is compared with a previously developed quasi-three-dimensional model by simulating the steady-state flow in an aquifer system in the Piceance Creek Basin, Colorado. The aquifer system consists of two aquifers separated by a leaky confining bed. The upper aquifer receives recharge from precipitation and is hydraulically connected to streams. For this problem, in order to make a valid comparison of results, a single layer was used to represent each aquifer. Furthermore, the need for a layer to represent the confining bed was eliminated by incorporating the effects of vertical leakage into the vertical component of the anisotropic hydraulic conductivity of the adjacent aquifers. Thus, the problem was represented by only two layers in each model with a total of about 2,100 equations. This restricted the effects of flow in the confining layer to the vertical component, but simulations with a third layer in the three-dimensional model permitting horizontal flow in the confining bed show that the two-layer approach is reasonable. Convergence to a solution of this problem takes about one minute of computer time on the IBM/155. This is about 30 times faster than the time required using the quasi-three-dimensional model.

Colorado↗
USGS1978

Uranium-isotope variations in groundwaters of the Floridan aquifer and Boulder Zone of south Florida

Water samples from four wells from the main Floridan aquifer (300-400 m below mean sea level) in southeast Florida exhibit 234 U 233 U activity ratios that are significantly lower than the secular equilibrium value of 1.00. Such anomalous values have been observed previously only in waters from sedimentary aquifers in the near-surface oxidizing environments. These four wells differ from six others, all producing from the same general horizon, in being located in cavernous highly transmissive zones. We hypothesize that the low activity ratios are indicative of a relic circulation pattern whereby water from the surface aquifer was channelled to lower levels when sea level was much lower. At a deeper cavernous level, known as the Boulder Zone (800-1,000 m below mean sea level), the U isotopes, along with other chemical constituents, show progressive changes with increasing distance from an inferred flow source in the Straits of Florida. This tends to support the hypothesized landward flow (though with a more northerly component) of cold seawater in the extensively transmissive Boulder Zone.

Florida↗
USGS1978

A field evaluation of subsurface and surface runoff. II. Runoff processes

Combined use of radioisotope tracer, flow rate, specific conductance and suspended-sediment measurements on a large field plot near Stanford, California, has provided more detailed information on surface and subsurface storm runoff processes than would be possible from any single approach used in isolation. Although the plot was surficially uniform, the runoff processes were shown to be grossly nonuniform, both spatially over the plot, and laterally and vertically within the soil. The three types of processes that have been suggested as sources of storm runoff (Horton-type surface runoff, saturated overland flow, and rapid subsurface throughflow) all occurred on the plot. The nonuniformity of the processes supports the partial- and variable-source area concepts. Subsurface storm runoff occurred in a saturated layer above the subsoil horizon, and short travel times resulted from flow through macropores rather than the soil matrix. Consideration of these observations would be necessary for physically realistic modeling of the storm runoff process.

California↗