Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Hydrogeology Journal”

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

Occurrence of volatile organic compounds in aquifers of the United States

Samples of ambient ground water were collected during 1985-2002 from 3,498 wells in 98 aquifer studies throughout the United States. None of the sampled wells were selected because of prior knowledge of nearby contamination. Most of these samples were analyzed for 55 volatile organic compounds (VOCs) to characterize their national occurrence. Volatile organic compounds were found in samples collected from 90 of the 98 aquifer studies. Occurrence frequencies of one or more VOCs for the 98 aquifer studies ranged from 0 to about 77% at an assessment level of 0.2 microgram per liter (??g/l). The aquifer studies with the largest occurrence frequencies were in southern Florida, southern New York, southern California, New Jersey, and Nevada. Trihalomethanes and solvents were the most frequently occurring VOC groups. Of the 55 VOCs included in this assessment, 42 occurred in at least one sample at an assessment level of 0.2 ??g/l. Chloroform, perchloroethene, and methyl tert-butyl ether were the most frequently occurring VOCs. Many factors, such as the hydrogeology of the aquifer, use of VOCs, land use, and the transport and fate properties of VOCs, affect the occurrence of VOCs in ground water. ?? 2008 American Water Resources Association.

Journal of the American Water Resources Associatio↗

A trans-dimensional Bayesian Markov chain Monte Carlo algorithm for model assessment using frequency-domain electromagnetic data

A meaningful interpretation of geophysical measurements requires an assessment of the space of models that are consistent with the data, rather than just a single, ‘best’ model which does not convey information about parameter uncertainty. For this purpose, a trans-dimensional Bayesian Markov chain Monte Carlo (MCMC) algorithm is developed for assessing frequencydomain electromagnetic (FDEM) data acquired from airborne or ground-based systems. By sampling the distribution of models that are consistent with measured data and any prior knowledge, valuable inferences can be made about parameter values such as the likely depth to an interface, the distribution of possible resistivity values as a function of depth and non-unique relationships between parameters. The trans-dimensional aspect of the algorithm allows the number of layers to be a free parameter that is controlled by the data, where models with fewer layers are inherently favoured, which provides a natural measure of parsimony and a significant degree of flexibility in parametrization. The MCMC algorithm is used with synthetic examples to illustrate how the distribution of acceptable models is affected by the choice of prior information, the system geometry and configuration and the uncertainty in the measured system elevation. An airborne FDEM data set that was acquired for the purpose of hydrogeological characterization is also studied. The results compare favorably with traditional least-squares analysis, borehole resistivity and lithology logs from the site, and also provide new information about parameter uncertainty necessary for model assessment.

Geophysical Journal International↗

Hydrogeologic controls on groundwater discharge and nitrogen loads in a coastal watershed

Submarine groundwater discharge (SGD) is a small portion of the global water budget, but a potentially large contributor to coastal nutrient budgets due to high concentrations relative to stream discharge. A numerical groundwater flow model of the Inland Bays Watershed, Delaware, USA, was developed to identify the primary hydrogeologic factors that affect groundwater discharge rates and transit times to streams and bays. The distribution of groundwater discharge between streams and bays is sensitive to the depth of the water table below land surface. Higher recharge and reduced hydraulic conductivity raised the water table and increased discharge to streams relative to bays compared to the Reference case (in which 66% of recharge is discharged to streams). Increases to either factor decreased transit times for discharge to both streams and bays compared to the Reference case (in which mean transit times are 56.5 and 94.3 years, respectively), though sensitivity to recharge is greater. Groundwater-borne nitrogen loads were calculated from nitrogen concentrations measured in discharging fresh groundwater and modeled SGD rates. These loads combined with long SGD transit times suggest groundwater-borne nitrogen reductions and estuarine water quality improvements will lag decades behind implementation of efforts to manage nutrient sources. This work enhances understanding of the hydrogeologic controls on and uncertainties in absolute and relative rates and transit times of groundwater discharge to streams and bays in coastal watersheds.

Journal of Hydrology↗

Multivariate analyses with end-member mixing to characterize groundwater flow: Wind Cave and associated aquifers

Principal component analysis (PCA) applied to hydrochemical data has been used with end-member mixing to characterize groundwater flow to a limited extent, but aspects of this approach are unresolved. Previous similar approaches typically have assumed that the extreme-value samples identified by PCA represent end members. The method presented herein is different from previous work in that (1) end members were not assumed to have been sampled but rather were estimated and constrained by prior knowledge; (2) end-member mixing was quantified in relation to hydrogeologic domains, which focuses model results on major hydrologic processes; (3) a method to select an appropriate number of end members using a series of cluster analyses is presented; and (4) conservative tracers were weighted preferentially in model calibration, which distributed model errors of optimized values, or residuals, more appropriately than would otherwise be the case. The latter item also provides an estimate of the relative influence of geochemical evolution along flow paths in comparison to mixing. This method was applied to groundwater in Wind Cave and the associated karst aquifer in the Black Hills of South Dakota, USA. The end-member mixing model was used to test a hypothesis that five different end-member waters are mixed in the groundwater system comprising five hydrogeologic domains. The model estimated that Wind Cave received most of its groundwater inflow from local surface recharge with an additional 33% from an upgradient aquifer. Artesian springs in the vicinity of Wind Cave primarily received water from regional groundwater flow.

South Dakota↗

Multi-lithofacies alluvial characterization via airborne electromagnetic-borehole fusion using ordinary interval kriging and geologic constraints

Airborne electromagnetic (AEM) survey provides extensive spatial coverage and detailed resolution at the near surface and can be used to develop hydrogeological models. However, utilization of AEM data is not straightforward because AEM resistivity is an indirect measurement for inferring sediment types. This study develops an ordinary interval kriging (OIK) algorithm and a resistivity-to-multi-lithofacies (R2ML) data fusion workflow for multi-lithofacies alluvial characterization. OIK utilizes irregular interval data to construct three-dimensional (3D) resistivity fields from one-dimensional inverted AEM resistivity models. The R2ML workflow maps the resistivity field generated from OIK into a multi-facies lithological model, incorporating geologic constraints derived from well logs and geological observations. The numerical and real-world cases demonstrate that OIK is computationally efficient, accounts for 3D anisotropy, and minimizes the smoothing effect, thereby preserving resistivity contrasts and reducing interpolation uncertainty. The methodology is applied to lithologic characterization of the Mississippi River Valley alluvial aquifer (MRVA) in the Shellmound area, Mississippi, U.S. A frequency-domain AEM survey was conducted to support groundwater studies for the managed aquifer recharge (MAR) to the MRVA. The resulting lithological model, including four types of lithofacies—clay, very fine sand, fine-medium sands, and graveliferous sands, illustrates the geomorphological processes of the MRVA and implies potential MAR. The alignment between the lithological model and existing geological and hydrogeological investigations demonstrates that OIK and R2ML workflow effectively capture the subsurface architecture of the MRVA. The methods have broad applicability for characterizing alluvial aquifers through AEM-borehole data fusion, supporting sustainable groundwater management.

Mississippi↗

Low transient storage and uptake efficiencies in seven agricultural streams: implications for nutrient demand

We used mass load budgets, transient storage modeling, and nutrient spiraling metrics to characterize nitrate (NO 3 − ), ammonium (NH 4 + ), and inorganic phosphorus (SRP) demand in seven agricultural streams across the United States and to identify in-stream services that may control these conditions. Retention of one or all nutrients was observed in all but one stream, but demand for all nutrients was low relative to the mass in transport. Transient storage metrics ( A s /A , F med 200 , T str , and q s ) correlated with NO 3 − retention but not NH 4 + or SRP retention, suggesting in-stream services associated with transient storage and stream water residence time could influence reach-scale NO 3 − demand. However, because the fraction of median reach-scale travel time due to transient storage ( F med 200 ) was ≤1.2% across the sites, only a relatively small demand for NO 3 − could be generated by transient storage. In contrast, net uptake of nutrients from the water column calculated from nutrient spiraling metrics were not significant at any site because uptake lengths calculated from background nutrient concentrations were statistically insignificant and therefore much longer than the study reaches. These results suggest that low transient storage coupled with high surface water NO 3 − inputs have resulted in uptake efficiencies that are not sufficient to offset groundwater inputs of N. Nutrient retention has been linked to physical and hydrogeologic elements that drive flow through transient storage areas where residence time and biotic contact are maximized; however, our findings indicate that similar mechanisms are unable to generate a significant nutrient demand in these streams relative to the loads.

Journal of Environmental Quality↗

Regional Aquifer-System Analysis Program of the U.S. Geological Survey, 1978-1992

The major ground-water systems of the United States have been investigated by the U.S. Geological Survey (USGS) through its Regional Aquifer-System Analysis (RASA) Program. During the first 15 years of the program (1978-92), 25 regional aquifer systems, including the most heavily pumped aquifers in the Nation, were intensively studied. As of mid-1992, 18 of the regional aquifer studies are completed or nearly so; 7 of the regional aquifer studies are ongoing, and compilation of a national ground-water atlas is in progress. This report summarizes the status of each RASA study and briefly describes the hydrology of the 25 regional aquifer systems. Important study results and examples of applications of study results are presented for some of the completed RASA investigations. The major contributions of the RASA Program are (1) assembly of data from numerous local studies and long-time data collections into systematic regional data bases; (2) comprehensive descriptions of the geologic, hydrologic, and geochemical characteristics of the regional aquifer systems; and (3) an understanding of how the regional ground-water-flow systems function under natural (predevelopment) and current (developed) conditions. To provide the comprehensive system descriptions, many of the RASA studies present, for the first time, maps depicting the hydrogeologic frameworks, water chemistry, potentiometric surfaces, and other aspects of entire regional aquifer systems. To provide an understanding of how the flow systems function, several of the completed RASA studies provide, for the first time, hydrologic budgets for both predevelopment and developed conditions. The results of the RASA Program are contained in nearly 900 reports published by the USGS, as well as various State and local agencies, and in articles published in scientific journals. A series of U.S. Geological Survey Professional Papers are being published to summarize and synthesize the results of each RASA study.

Circular↗

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

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

Nebraska↗

Low resistivity and permeability in actively deforming shear zones on the San Andreas Fault at SAFOD

The San Andreas Fault Observatory at Depth (SAFOD) scientific drillhole near Parkfield, California crosses the San Andreas Fault at a depth of 2.7 km. Downhole measurements and analysis of core retrieved from Phase 3 drilling reveal two narrow, actively deforming zones of smectite-clay gouge within a roughly 200 m-wide fault damage zone of sandstones, siltstones and mudstones. Here we report electrical resistivity and permeability measurements on core samples from all of these structural units at effective confining pressures up to 120 MPa. Electrical resistivity (~10 ohm-m) and permeability (10-21 to 10-22 m2) in the actively deforming zones were one to two orders of magnitude lower than the surrounding damage zone material, consistent with broader-scale observations from the downhole resistivity and seismic velocity logs. The higher porosity of the clay gouge, 2 to 8 times greater than that in the damage zone rocks, along with surface conduction were the principal factors contributing to the observed low resistivities. The high percentage of fine-grained clay in the deforming zones also greatly reduced permeability to values low enough to create a barrier to fluid flow across the fault. Together, resistivity and permeability data can be used to assess the hydrogeologic characteristics of the fault, key to understanding fault structure and strength. The low resistivities and strength measurements of the SAFOD core are consistent with observations of low resistivity clays that are often found in the principal slip zones of other active faults making resistivity logs a valuable tool for identifying these zones.

Journal of Geophysical Research↗

Soviet books and publications on hydrology (continental) and hydrogeology: titles and some notes on obtaining Soviet monographs

A common method of publication for Soviet scientists, which partly supplants periodicals, is the publication of a collection of articles on a general area of research, frequently by members of a given institution. An extensive sampling of world geologic literature for 1961 (Hawkes, 1966) showed that 33 percent of Soviet titles appeared in periodicals whereas 55 percent of North American and 70 percent of Western European literature appeared in this form. The Soviet predilection for symposia and collections of papers makes searching for information on a given subject more difficult for Westerners because the monographs in question are often not included in exchange agreements (except informal personal ones) with Western libraries and institutions, because they may be primed in small editions, and because such publications frequently escape the notice of Western abstract journals. Unless one is fortunate enough to have many personal contacts in the Soviet Union, there seems to be little alternative to at least a rudimentary knowledge of Russian in order to stay abreast of work published as monographs and in collections.

Union Of Soviet Socialist Republics↗

Interactions between ground water and wetlands, southern shore of Lake Michigan, USA

Wetlands between, and within, dune-beach complexes along the south shore of Lake Michigan are strongly affected by ground water. The hydrogeology of the glacial drift aquifer system in a 26 km 2 area was investigated to determine the effects of ground water on the hydrology and hydrochemistry of Cowles Bog and its adjacent wetlands. The investigation showed that ground water from intermediate- and regional-scale flow systems discharges to Cowles Bog from confined aquifers that underlie the wetland. These flow systems are recharged in moraines south of the dune-beach complexes. Water from the confined aquifers discharges into the surficial aquifer mainly by upward leakage through a buried till sheet that serves as the confining layer. However, the till sheet is breached below a raised peat mound in Cowles Bog, allowing direct upward discharge from the confinef aquifer into the surficial sand, marl, and peat. The shallow ground and wetland water in the area influenced by this leakage is a calcium magnesium bicarbonate type, with low tritium concentrations consistent with mixing of older ground water and more recent precipitation. Ground water and wetland water from surrounding areas are less mineralized and have higher tritium concentrations characteristic of precipitation in the late 1970s. The results of this study suggest that wetlands in complex hydrogeologic settings may be influenced by multiple ground-water flow systems that are affected by geomorphic features, stratigraphic discontinuities, and changes in sediment types. Discharge and recharge zones may both occur in the same wetland. Multidisciplinary studies incorporating hydrological, hydrochemical, geophysical, and sedimentological data are necessary to identify such complexities in wetland hydrology.

Journal of Hydrology↗

Generalized analytical model for benthic water flux forced by surface gravity waves

A generalized analytical model for benthic water flux forced by linear surface gravity waves over a series of layered hydrogeologic units is developed by adapting a previous solution for a hydrogeologic unit with an infinite thickness (Case I) to a unit with a finite thickness (Case II) and to a dual-unit system (Case III). The model compares favorably with laboratory observations. The amplitude of wave-forced benthic water flux is shown to be directly proportional to the amplitude of the wave, the permeability of the hydrogeologic unit, and the wave number and inversely proportional to the kinematic viscosity of water. A dimensionless amplitude parameter is introduced and shown to reach a maximum where the product of water depth and the wave number is 1.2. Submarine groundwater discharge (SGD) is a benthic water discharge flux to a marine water body. The Case I model estimates an 11.5-cm/d SGD forced by a wave with a 1 s period and 5-cm amplitude in water that is 0.5-m deep. As this wave propagates into a region with a 0.3-m-thick hydrogeologic unit, with a no-flow bottom boundary, the Case II model estimates a 9.7-cm/d wave-forced SGD. As this wave propagates into a region with a 0.2-m-thick hydrogeologic unit over an infinitely thick, more permeable unit, the Case III quasi-confined model estimates a 15.7-cm/d wave-forced SGD. The quasi-confined model has benthic constituent flux implications in coral reef, karst, and clastic regions. Waves may undermine tracer and seepage meter estimates of SGD at some locations. Copyright 2009 by the American Geophysical Union.

Journal of Geophysical Research C: Oceans↗

Environmental tracer evidence for connection between shallow and bedrock aquifers and high intrinsic susceptibility to contamination of the conterminous U.S. glacial aquifer

Covering a large portion of the northern conterminous United States (1.87 x 106 km2), the glacial aquifer serves as the primary water supply for 39 million public and domestic water users. Mean groundwater age, groundwater age distribution, and susceptibility to land surface contamination, using a new metric (Susceptibility Index; SI) based on the full age distribution and less prone to bias than estimated mean age, is reported for 168 public and domestic wells across the aquifer. Comparison of groundwater age metrics between well networks of varying spatial scale suggest an extensive sample network of equally spaced, long screened interval wells can be used to characterize aquifer wide groundwater age. Estimated mean age ranges from 1 to 50,000 years and, according to the composite age distribution, approximately 63 percent of all sampled water recharged after 1950 (i.e., modern) and 18 percent of the sampled water was recharged greater than 10,000 years ago. The later finding strongly suggests a connection between the glacial aquifer and underlying bedrock aquifers. Statistical analysis of glacial aquifer hydrogeology and age metrics show groundwater ages are young (less than few 100 years) and more susceptible to land surface contamination (larger SI) in unconfined and shallow portions of the aquifer. Old groundwater (greater than 1000 years) is more often associated with thicker sequences of fine grain sediments and/or shallow bedrock. Calculated SI is shown to be more strongly related to the number of land surface contaminants detected than mean age or fraction modern. Statistical analysis of SI and hydrogeology indicates SI is largely dictated by well depth and confinement. This study demonstrates how sample network design can be used to characterize groundwater age of large aquifers with a limited number of samples and how interpretation of environmental tracers can be used to improve conceptual models of groundwater aquifers and identify groundwater susceptible to contamination.

Glacial aquifer↗

Transport and fate of nitrate at the ground-water/surface-water interface

Although numerous studies of hyporheic exchange and denitrification have been conducted in pristine, high-gradient streams, few studies of this type have been conducted in nutrient-rich, low-gradient streams. This is a particularly important subject given the interest in nitrogen (N) inputs to the Gulf of Mexico and other eutrophic aquatic systems. A combination of hydrologic, mineralogical, chemical, dissolved gas, and isotopic data were used to determine the processes controlling transport and fate of NO 3 − in streambeds at five sites across the USA. Water samples were collected from streambeds at depths ranging from 0.3 to 3 m at three to five points across the stream and in two to five separate transects. Residence times of water ranging from 0.28 to 34.7 d m −1 in the streambeds of N-rich watersheds played an important role in allowing denitrification to decrease NO 3 − concentrations. Where potential electron donors were limited and residence times were short, denitrification was limited. Consequently, in spite of reducing conditions at some sites, NO 3 − was transported into the stream. At two of the five study sites, NO 3 − in surface water infiltrated the streambeds and concentrations decreased, supporting current models that NO 3 − would be retained in N-rich streams. At the other three study sites, hydrogeologic controls limited or prevented infiltration of surface water into the streambed, and ground-water discharge contributed to NO 3 − loads. Our results also show that in these low hydrologic-gradient systems, storm and other high-flow events can be important factors for increasing surface-water movement into streambeds.

Journal of Environmental Quality↗

Groundwater inflow toward a preheated volcanic conduit: Application to the 2018 eruption at Kīlauea Volcano, Hawai’i

The many successes in volcano forecasting over the past several decades owe mainly to pattern recognition, both in monitoring data and the geologic record. During the early stages of the 2018 Kīlauea eruption, the conceptual model of Stearns (1925), based on the explosive 1924 Kīlauea eruption, was highly influential. This model postulates that explosions are triggered by liquid-water inflow into a recently vacated magma conduit. Modern quantitative modeling approaches, supplemented by hydrogeologic data unavailable in 1925, yield a more nuanced view. Results demonstrate that liquid-water inflow would likely be delayed by months to years, owing to the inability of liquid water to transit a zone of very hot rock surrounding the conduit. The exercise demonstrates the use of physically based modeling to supplement traditional volcano-forecasting approaches during an ongoing event.

Hawaii↗

Groundwater flow into Lake Michigan from Wisconsin

Detailed hydrogeological study has been done at six sites along the Lake Michigan shoreline in Wisconsin. At each site a flux of groundwater to the lake has been calculated for both natural conditions and the existing conditions created by pumping. The values from each site have then been extrapolated to the entire portion of the total shoreline having similar hydrogeology in order to calculate a total flow of groundwater to the lake. Sensitivity analysis with a digital model was used to define limits on the similarity of hydrogeologic conditions. The net flow calculated is 580–880 m 3 day −1 km −1 of shoreline, which falls within the previously published range of 110–8200 m 3 day −1 km −1 . Human activity may have reduced the natural flow as much as 15%. The estimated natural flow is between 7 and 11% of the surface water contribution to the lake from the study area.

Wisconsin↗

A comparison of the coupled fresh water-salt water flow and the Ghyben-Herzberg sharp interface approaches to modeling of transient behavior in coastal aquifer systems

A quasi-three dimensional finite difference model which simulates coupled, fresh water and salt water flow, separated by a sharp interface, is used to investigate the effects of storage characteristics, transmissivity, boundary conditions and anisotropy on the transient responses of such flow systems. The magnitude and duration of the departure of aquifer response from the behavior predicted using the Ghyben-Herzberg, one-fluid approach is a function of the ease with which flow can be induced in the salt water region. In many common hydrogeologic settings short-term fresh water head responses, and transitional responses between short-term and long-term, can only be realistically reproduced by including the effects of salt water flow on the dynamics of coastal flow systems. The coupled fresh water-salt water flow modeling approach is able to reproduce the observed annual fresh water head response of the Waialae aquifer of southeastern Oahu, Hawaii.

Hawaii↗

Chemical hydrogeology in natural and contaminated environments

Chemical hydrogeology, including organic and inorganic aspects, has contributed to an increased understanding of groundwater flow systems, geologic processes, and stressed environments. Most of the basic principles of inorganic-chemical hydrogeology were first established by investigations of organic-free, regional-scale systems for which simplifying assumptions could be made. The problems of groundwater contamination are causing a shift of emphasis to microscale systems that are dominated by organic-chemical reactions and that are providing an impetus for the study of naturally occurring and manmade organic material. Along with the decrease in scale, physical and chemical heterogeneity become major controls. Current investigations and those selected from the literature demonstrate that heterogeneity increases in importance as the study site decreases from regional-scale to macroscale to microscale. Increased understanding of regional-scale flow systems is demonstrated by selection of investigations of carbonate and volcanic aquifers to show how applications of present-day concepts and techniques can identify controlling chemical reactions and determine their rates; identify groundwater flow paths and determine flow velocity; and determine aquifer characteristics. The role of chemical hydrogeology in understanding geologic processes of macroscale systems is exemplified by selection of investigations in coastal aquifers. Phenomena associated with the mixing zone generated by encroaching sea water include an increase in heterogeneity of permeability, diagenesis of minerals, and formation of geomorphic features, such as caves, lagoons, and bays. Ore deposits of manganese and uranium, along with a simulation model of ore-forming fluids, demonstrate the influence of heterogeneity and of organic compounds on geochemical reactions associated with genesis of mineral deposits. In microscale environments, importance of heterogeneity and consequences of organic reactions in determining the distribution and concentrations cf. constituents are provided by several studies, including infiltration of sewage effluent and migration of creosote in coastal plain aquifers. These studies show that heterogeneity and the dominance of organically controlled reactions greatly increase the complexity of investigations.Current investigations and those selected from the literature demonstrate that heterogeneity increases in importance as the study site decreases from regional-scale to macroscale to microscale. Increased understanding of regional-scale flow systems is demonstrated by selection of investigations of carbonate and volcanic aquifers to show how application of present-day concepts and techniques can identify controlling chemical reactions and determine their rates; identify groundwater flow paths and determine flow velocity; and determine aquifer characteristics. The role of chemical hydrogeology in understanding geologic processes of macroscale systems is exemplified by selection of investigations in coastal aquifers. Ore deposits of manganese and uranium, along with a simulation model of ore-forming fluids, demonstrate the influence of heterogeneity and of organic compounds on geochemical reactions associated with genesis of mineral deposits. In microscale environments, importance of heterogeneity is illustrated by studies of infiltration of sewage effluent and migration of creosote in coastal plain aquifers.

Journal of Hydrology↗