Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Environmental Geosciences”

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 37 records · Page 2Linked to original sources
USGS2020

Modeling geologic sequestration of carbon dioxide in a deep saline carbonate reservoir with TOUGH2–ChemPlugin, a new tool for reactive transport modeling

This paper outlines the development and demonstration of a new tool, TOUGH2–ChemPlugin (T2CPI) for predicting rock–water–CO 2 interaction following injection of supercritical CO 2 into a heterogeneous carbonate system. Specifically, modeling capabilities of TOUGH2, which examines multiphase flow and supercritical CO 2 behavior, were combined with the geochemical modeling capabilities of The Geochemist’s Workbench ® (GWB), using ChemPlugin TM . ChemPlugin is a self-linking re-entrant software object that, when coupled to a transport simulator, retains the flow and transport capabilities of the simulator but enables incorporation of reactive chemistry via GWB. To test and assess the capabilities of T2CPI, results from T2CPI simulations were compared to those of TOUGHREACT, using the same carbonate reservoir parameters (based on the Dollar Bay Formation of the South Florida Basin). Overall, results of simulations from TOUGHREACT and T2CPI were very similar for nearly all evaluated parameters. Dissimilarities between the two programs included qualitative differences in how TOUGHREACT and T2CPI predicted calcite dissolution and the subsequent spatial pattern of the porosity gain caused by how each handles evaporation of water near the injection point. The TOUGHREACT program is a proven, widely used tool for evaluating CO 2 –brine–rock interaction following supercritical CO 2 injection. The T2CPI tool offers similar capabilities and strengths of TOUGHREACT, with the ability to read in and use databases for a wide range of activity coefficient types. This program also has abilities to use a wide range of kinetic constraints, define those kinetic constraints with scripts or compiled libraries, account for colloidal transport, and/or account for a wide range of surface sorption models.

Florida↗
USGS2014

Reconstructing suspended sediment mercury contamination of a steep, gravel-bed river using reservoir theory

We use sediment ages and mercury (Hg) concentrations to estimate past and future concentrations in the South River, Virginia, where Hg was released between 1930 and 1950 from a manufacturing process related to nylon production. In a previous study, along a 40 km (25 mi) reach, samples were collected from 26 of 54 fine-grained deposits that formed in the lee of large wood obstructions in the channel and analyzed for grain size, Hg concentration, and organic content. We also obtained radiometric dates from six deposits. To create a history that reflects the full concentration distribution (which contains concentrations as high as 900 mg/kg [900 ppm]), here, we treat the deposits as a single reservoir exchanging contaminated sediments with the overlying water column, and assume that the total sediment mass in storage and the distribution of sediment ages are time invariant. We use reservoir theory to reconstruct the annual history of Hg concentration on suspended sediment using data from our previous study and new results presented here. Many different reconstructed histories fit our data. To constrain results, we use information from a well-preserved core (and our estimate of the total mass of Hg stored in 2007) to specify the years associated with the peak concentration of 900 mg/kg. Our results indicate that around 850 kg (1874 lb) of Hg was stored in the deposits between 1955 and 1961, compared to only 80 kg (176 lb) today. Simulations of future Hg remediation suggest that 100-yr timescales will be needed for the South River to remove Hg-contaminated sediments from the channel perimeter through natural processes.

Environmental Geosciences↗
USGS2005

Environmental impacts of oil production on soil, bedrock, and vegetation at the U.S. Geological Survey Osage-Skiatook Petroleum Environmental Research site A, Osage County, Oklahoma

The U.S. Geological Survey is investigating the impacts of oil and gas production on soils, groundwater, surface water, and ecosystems in the United States. Two sites in northeastern Oklahoma (sites A and B) are presently being investigated under the Osage–Skiatook Petroleum Environmental Research project. Oil wells on the lease surrounding site A in Osage County, Oklahoma, produced about 100,000 bbl of oil between 1913 and 1981. Prominent production features on the 1.5-ha (3.7-ac) site A include a tank battery, an oil-filled trench, pipelines, storage pits for both produced water and oil, and an old power unit. Site activities and historic releases have left open areas in the local oak forest adjacent to these features and a deeply eroded salt scar downslope from the pits that extends to nearby Skiatook Lake. The site is underlain by surficial sediments comprised of very fine-grained eolian sand and colluvium as much as 1.4 m (4.6 ft) thick, which, in turn, overlie flat-lying, fractured bedrock comprised of sandstone, clayey sandstone, mudstone, and shale. A geophysical survey of ground conductance and concentration measurements of aqueous extracts (1:1 by weight) of core samples taken in the salt scar and adjacent areas indicate that unusual concentrations of NaCl-rich salt are present at depths to at least 8 m (26 ft) in the bedrock; however, little salt occurs in the eolian sand. Historic aerial photographs, anecdotal reports from oil-lease operators, and tree-ring records indicate that the surrounding oak forest was largely established after 1935 and thus postdates the majority of surface damage at the site. Blackjack oaks adjacent to the salt scar have anomalously elevated chloride ( 400 ppm) in their leaves and record the presence of NaCl-rich salt or salty water in the shallow subsurface. The geophysical measurements also indicate moderately elevated conductance beneath the oak forest adjoining the salt scar.

Oklahoma↗
USGS2005

Progression of methanogenic degradation of crude oil in the subsurface

Our results show that subsurface crude-oil degradation rates at a long-term research site were strongly influenced by small-scale variations in hydrologic conditions. The site is a shallow glacial outwash aquifer located near Bemidji in northern Minnesota that became contaminated when oil spilled from a broken pipeline in August 1979. In the study area, separate-phase oil forms a subsurface oil body extending from land surface to about 1 m (3.3 ft) below the 6–8-m (20–26 ft)-deep water table. Oil saturation in the sediments ranges from 10–20% in the vadose zone to 30–70% near the water table. At depths below 2 m (6.6 ft), degradation of the separate-phase crude oil occurs under methanogenic conditions. The sequence of methanogenic alkane degradation depletes the longer chain n-alkanes before the shorter chain n-alkanes, which is opposite to the better known aerobic sequence. The rates of degradation vary significantly with location in the subsurface. Oil-coated soils within 1.5 m (5 ft) of land surface have experienced little degradation where soil water saturation is less than 20%. Oil located 2–8 m (6.6–26 ft) below land surface in areas of higher recharge has been substantially degraded. The best explanation for the association between recharge and enhanced degradation seems to be increased downward transport of microbial growth nutrients to the oil body. This is supported by observations of greater microbial numbers at higher elevations in the oil body and significant decreases with depth in nutrient concentrations, especially phosphorus. Our results suggest that environmental effects may cause widely diverging degradation rates in the same spill, calling into question dating methods based on degradation state.

Environmental Geosciences↗
USGS2005

Impacts of petroleum production on ground and surface waters: Results from the Osage-Skiatook Petroleum Environmental Research A site, Osage County Oklahoma

Kharaka has been a research hydrogeochemist with the U.S. Geological Survey since1975. His current research covers the fields of water-rock-gas interactions, CO 2 sequestration, contamination from agricultural drainage and petroleum produced water, and fluid-fault interactions. Kharaka has authored more than 100 scientific papers and book chapters and has delivered about 200 presentations. He received his Ph.D. from the University of California, Berkeley, in 1971. Thordsen is a geologist, employed with the U.S. Geological Survey since 1981 and presently working in the project team Hydrogeochemistry of Fine Grained Sediments, in the Water Resources Division, National Research Program. His recent research covers water-rock interactions and geochemical modeling in geothermal systems, fault-fluid systems, and petroleum produced water brine impacts. He received his M.S. in geology from the Ohio State University in 1988. Kakouros has been a hydrogeologist in the U.S. Geological Survey Water Resources Division in Menlo Park, California, since 1999. His current research covers the fields of Se contamination from agricultural drainage, CO 2 sequestration, and contamination from petroleum and petroleum-produced water. Kakouros received his M.Sc. from the San Jose State University in 2001. Herkelrath is a physicist in the U.S. Geological Survey's Water Resources Division National Research Program in Menlo Park, California. Herkelrath received his Ph.D. in physics from the University of Wisconsin at Madison. His recent research interests include multiphase flow and solute transport modeling in shallow ground-water systems where oil or brine is present.

Oklahoma↗
USGS2006

Measuring permanence of CO 2 storage in saline formations: The Frio experiment

If CO 2 released from fossil fuel during energy production is returned to the subsurface, will it be retained for periods of time significant enough to benefit the atmosphere? Can trapping be assured in saline formations where there is no history of hydrocarbon accumulation? The Frio experiment in Texas was undertaken to provide answers to these questions. One thousand six hundred metric tons of CO 2 were injected into the Frio Formation, which underlies large areas of the United States Gulf Coast. Reservoir characterization and numerical modeling were used to design the experiment, as well as to interpret the results through history matching. Closely spaced measurements in space and time were collected to observe the evolution of immiscible and dissolved CO 2 during and after injection. The high-permeability, steeply dipping sandstone allowed updip flow of supercritical CO 2 as a result of the density contrast with formation brine and absence of a local structural trap. The front of the CO 2 plume moved more quickly than had been modeled. By the end of the 10-day injection, however, the plume geometry in the plane of the observation and injection wells had thickened to a distribution similar to the modeled distribution. As expected, CO 2 dissolved rapidly into brine, causing pH to fall and calcite and metals to be dissolved. Postinjection measurements, including time-lapse vertical seismic profiling transects along selected azimuths, cross-well seismic topography, and saturation logs, show that CO 2 migration under gravity slowed greatly 2 months after injection, matching model predictions that significant CO 2 is trapped as relative permeability decreases.

Environmental Geosciences↗
USGS2003

Causes of hot-spot wetland loss in the Mississippi delta plain

Field surveys and sediment cores were used to estimate marsh erosion and land subsidence at Madison Bay, a well-known wetland loss hot spot in coastal Louisiana. Former marshes of Madison Bay are under about 1 m of water. Nearly two-thirds of the permanent flooding was caused by rapid subsidence in the late 1960s, whereas the other third was caused by subsequent erosion. Subsidence rates near Madison Bay since the 1960s ( 20 mm/yr) are an order of magnitude greater than deltaic subsidence rates averaged for the past 400 4000 yr ( 2 mm/yr). The rapid acceleration and unexpected decline in wetland losses in the Mississippi delta plain are difficult to explain on the basis of most physical and biogeochemical processes. There are, however, close temporal and spatial correlations among regional wetland loss, high subsidence rates, and large-volume fluid production from nearby hydrocarbon fields. The decreased rates of wetland loss since the 1970s may be related to decreased rates of subsidence caused by significantly decreased rates of subsurface fluid withdrawal. Annual fluid production from the Lapeyrouse, Lirette, and Bay Baptiste fields that encompass Madison Bay accelerated in the 1960s, peaked about 1970, and then declined abruptly. Large decreases in pore pressure in the Lapeyrouse field have likely altered subsurface stresses and reactivated a major fault that coincides with the wetland loss hot spot. Therefore, wetland losses at Madison Bay can be closely linked to rapid subsidence and possible fault reactivation induced by long-term, large-volume hydrocarbon production.

Environmental Geosciences↗
USGS1999

Coal-fired power generaion, new air quality regulations, and future U.S. coal production

Tighter new regulation of stack gas emissions and competition in power generation are driving electrical utilities to demand cleaner, lower sulfur coal. Historical data on sulfur content of produced coals shows little variability in coal quality for individual mines and individual coal-producing counties over relatively long periods of time. If coal-using power generators follow the compliance patterns established in Phase I of the 1990 Clean Air Act Amendments, then the industry's response to the tighter Phase II emissions standards will result in large amounts of coal production shifting from higher sulfur areas to areas with lower cost low sulfur coal. One reason this shift will likely occur is that currently only 30% of U.S. coal-fired electrical generating capacity is equipped with flue-gas scrubbers. In 1995, coal mines in the higher sulfur areas of the Illinois Basin and Northern and Central Appalachia employed 78% of all coal miners (>70,000 miners). A substantial geographical redistribution of the nation's coal supplies will likely lead to economic dislocations that will reach beyond local coal-producing areas.

Environmental Geosciences↗
USGS2004

Probabilistic assessment of precipitation-triggered landslides using historical records of landslide occurrence, Seattle, Washington

Ninety years of historical landslide records were used as input to the Poisson and binomial probability models. Results from these models show that, for precipitation-triggered landslides, approximately 9 percent of the area of Seattle has annual exceedance probabilities of 1 percent or greater. Application of the Poisson model for estimating the future occurrence of individual landslides results in a worst-case scenario map, with a maximum annual exceedance probability of 25 percent on a hillslope near Duwamish Head in West Seattle. Application of the binomial model for estimating the future occurrence of a year with one or more landslides results in a map with a maximum annual exceedance probability of 17 percent (also near Duwamish Head). Slope and geology both play a role in localizing the occurrence of landslides in Seattle. A positive correlation exists between slope and mean exceedance probability, with probability tending to increase as slope increases. Sixty-four percent of all historical landslide locations are within 150 m (500 ft, horizontal distance) of the Esperance Sand/Lawton Clay contact, but within this zone, no positive or negative correlation exists between exceedance probability and distance to the contact.

Washington↗
USGS2004

Transient hazard model using radar data for predicting debris flows in Madison County, Virginia

During the rainstorm of June 27, 1995, roughly 330-750 mm of rain fell within a 16-hour period, initiating floods and over 600 debris flows in a small area (130 km2) of Madison County, VA. We developed a distributed version of Iverson's transient response model for regional slope stability analysis for the Madison County debris flows. This version of the model evaluates pore-pressure head response and factor of safety on a regional scale in areas prone to rainfall-induced shallow (<2-3 m) landslides. These calculations used soil properties of shear strength and hydraulic conductivity from laboratory measurements of soil samples collected from field sites where debris flows initiated. Rainfall data collected by radar every 6 minutes provided a basis for calculating the temporal variation of slope stability during the storm. The results demonstrate that the spatial and temporal variation of the factor of safety correlates with the movement of the storm cell. When the rainstorm was treated as two separate rainfall events and a larger hydraulic conductivity and friction angle than the laboratory values were used, the timing and location of landslides predicted by the model were in closer agreement with eyewitness observations of debris flows. Application of spatially variable initial pre-storm water table depth and soil properties may improve both the spatial and temporal prediction of instability.

Environmental & Engineering Geoscience↗
USGS2005

Colorado Late Cenozoic Fault and Fold Database and Internet Map Server: User-friendly technology for complex information

Internet mapping applications for geologic data allow simultaneous data delivery and collection, enabling quick data modification while efficiently supplying the end user with information. Utilizing Web-based technologies, the Colorado Geological Survey's Colorado Late Cenozoic Fault and Fold Database was transformed from a monothematic, nonspatial Microsoft Access database into a complex information set incorporating multiple data sources. The resulting user-friendly format supports easy analysis and browsing. The core of the application is the Microsoft Access database, which contains information compiled from available literature about faults and folds that are known or suspected to have moved during the late Cenozoic. The database contains nonspatial fields such as structure type, age, and rate of movement. Geographic locations of the fault and fold traces were compiled from previous studies at 1:250,000 scale to form a spatial database containing information such as length and strike. Integration of the two databases allowed both spatial and nonspatial information to be presented on the Internet as a single dataset (http://geosurvey.state.co.us/pubs/ceno/). The user-friendly interface enables users to view and query the data in an integrated manner, thus providing multiple ways to locate desired information. Retaining the digital data format also allows continuous data updating and quick delivery of newly acquired information. This dataset is a valuable resource to anyone interested in earthquake hazards and the activity of faults and folds in Colorado. Additional geologic hazard layers and imagery may aid in decision support and hazard evaluation. The up-to-date and customizable maps are invaluable tools for researchers or the public.

Environmental & Engineering Geoscience↗
USGS2005

Assessing the vulnerability of a municipal well field to contamination in a karst aquifer

Proposed expansion of extractive lime-rock mines near the Miami-Dade County Northwest well field and Everglades wetland areas has garnered intense scrutiny by government, public, environmental stakeholders, and the media because of concern that mining will increase the risk of pathogen contamination. Rock mines are excavated to the same depth as the well field's primary producing zone. The underlying karst Biscayne aquifer is a triple-porosity system characterized by (1) a matrix of interparticle porosity and separate vug porosity; (2) touching-vug porosity that forms preferred, stratiform passageways; and, less commonly, (3) conduit porosity formed by thin solution pipes, bedding-plane vugs, and cavernous vugs. Existing ground-water flow and particle tracking models do not provide adequate information regarding the ability the aquifer to limit the advective movement of pathogens and other contaminants. Chemical transport and colloidal mobility properties have been delineated using conservative and microsphere-surrogate tracers for Cryptosporidium parvum. Forced-gradient tests were executed by introducing conservative tracers into injection wells located 100 m (328 ft) from a municipal-supply well. Apparent mean advective velocity between the wells is one to two orders of magnitude greater than previously measured. Touching-vug, stratiform flow zones are efficient pathways for tracer movement at the well field. The effective porosity for a continuum model between the point of injection and tracer recovery ranges from 2 to 4 percent and is an order of magnitude smaller than previously assumed. Existing well-field protection zones were established using porosity estimates based on specific yield. The effective, or kinematic, porosity of a Biscayne aquifer continuum model is lower than the total porosity, because high velocities occur along preferential flow paths that result in faster times of travel than can be represented with the ground-water flow equation. Tracer tests indicate that the relative ease of contaminant movement to municipal supply wells is much greater than previously considered.

Environmental & Engineering Geoscience↗
USGS2005

Using multiple chemical indicators to assess sources of nitrate and age of groundwater in a karstic spring basin

Human health and ecological concerns have arisen due to a steady increase in nitrate-N concentrations during the past 40 years in Fannin Springs (0.3-4.7 mg/L), a regional discharge point with an average flow of >2.8 m3/second (>100 ft3/second) for water from the karstic Upper Floridan aquifer (UFA). Multiple chemical indicators (major dissolved species, 15N and 18O of nitrate, dissolved gases, 78 pesticides and degradates, and 67 organic compounds typically found in domestic and industrial wastewater) and transient tracers (3H/3He, chlorofluorocarbons [CFCs], sulfur hexafluoride [SF6]) were analyzed in water samples from nine wells along three transects and in spring water to assess groundwater age and potential contaminant sources. Land use is predominantly agricultural (52 percent) and forest (31 percent) in the 320 km2 (124 mi2) spring basin, which was delineated from a potentiometric-surface map of the UFA using high-resolution water-level data. Nitrate-N concentrations were highly variable in the oxic UFA and ranged from <0.02 to 4.7 mg/L. ?? 15N-NO3 values (3.4-9.9 per mil) indicated that nitrate contamination originated from inorganic sources (synthetic fertilizer) and organic sources (manure spreading or waste disposal). Higher nitrate concentrations and the younger age of spring water relative to water from upgradient wells indicate better communication with N sources at the surface. Apparent ages of groundwater correlated positively with well depth (P < 0.05) and were younger in water from wells nearer to the spring (<8 years) compared with other wells (10-50 years). Most transient tracer concentrations were consistent with binary mixing curves representing mixtures of water recharged during the past 10 years and older water (recharged before 1940). Young water mixing fractions ranged from 0.07 to 0.90. Trace levels of herbicides found in groundwater and spring water were indicative of applications for vegetative control in agricultural and other land-use types.

Environmental & Engineering Geoscience↗
USGS2002

Gravitational failure of sea cliffs in weakly lithified sediment

Gravitational failure of sea cliffs eroded into weakly lithified sediment at several sites in California involves episodic stress-release fracturing and cantilevered block falls. The principal variables that influence the gravitational stability are tensional stresses generated during the release of horizontal confining stress and weakening of the sediment with increased saturation levels. Individual failures typically comprise less than a cubic meter of sediment, but large areas of a cliff face can be affected by sustained instability over a period of several days. Typically, only the outer meter or so of sediment is removed during a failure episode. In-place sediment saturation levels vary over time and space, generally being higher during the rainy season but moderate to high year-round. Laboratory direct-shear tests show that sediment cohesion decreases abruptly with increasing saturation level; the decrease is similar for all tested sediment if the cohesion is normalized by the maximum, dry-sediment cohesion. Large failures that extend over most or all of the height of the sea cliff are uncommon, but a few large wedge-shaped failures sometimes occur, as does separation of large blocks at sea cliff-gully intersections.

Environmental & Engineering Geoscience↗
USGS2002

Comparing landslide inventories: The map depends on the method

Landslide inventory maps are generally prepared by interpreting the geomorphic expression of landsliding on aerial photos, topographic maps, or on the ground. Distinctive landslide geomorphology allows the recognition and mapping of landslides, although there are always landslides that have very subtle expression and are not identified. The difficulties of mapping landslides based on their geomorphic expression are amplified in heavily forested terrain. The ground surface is obscured by tree cover on aerial photographs, and landslide-related features are often hidden. This limitation affects not only aerial photo interpretation, but also interpretation of topographic maps, which are based on aerial photographs. We compared five maps showing landslides in the Laurel Quadrangle in the Santa Cruz Mountains, California. These include a geologic map, a map prepared for the county based on interpretation of aerial photographs, a map prepared by us based on aerial photographs and compilation of previous work, a map of features interpreted from the U.S. Geological Survey 7.5-minute topographic map, and a detailed field-based landslide map. Comparison of these maps shows that the geologic map identifies few landslides, but most landslides on the geologic map are also shown on the other maps. The two maps based mainly on aerial photo interpretation tend to show the larger slides, but there is only about 60 percent correspondence of landslide areas between the two. Comparing the reconnaissance techniques with the much more detailed field mapping shows that the reconnaisance maps emphasize the large slides of bedrock and identify a lower percentage of shallow debris slides and debris flows.

Environmental & Engineering Geoscience↗
USGS2003

Geomorphic and hydrologic assessment of erosion hazards at the Norman municipal landfill, Canadian River floodplain, central Oklahoma

The Norman, Oklahoma, municipal landfill closed in 1985 after 63 years of operation, because it was identified as a point source of hazardous leachate composed of organic and inorganic compounds. The landfill is located on the floodplain of the Canadian River, a sand-bed river characterized by erodible channel boundaries and by large variation in mean monthly discharges. In 1986, floodwaters eroded riprap protection at the southern end of the landfill and penetrated the landfill's clay cap, thereby exposing the landfill contents. The impact of this moderate-magnitude flood event (Q 12 ) was the catalyst to investigate erosion hazards at the Norman landfill. This geomorphic investigation analyzed floodplain geomorphology and historical channel changes, flood-frequency distributions, an erosion threshold, the geomorphic effectiveness of discharge events, and other factors that influence erosion hazards at the landfill site. The erosion hazard at the Norman landfill is a function of the location of the landfill with respect to the channel thalweg, erosional resistance of the channel margins, magnitude and duration of discrete discharge events, channel form and hydraulic geometry, and cumulative effects related to a series of discharge events. Based on current climatic conditions and historical channel changes, a minimum erosion threshold is set at bankfull discharge (Q = 572 m 3 /s). The annual probability of exceeding this threshold is 0.53. In addition, this analysis indicates that peak stream power is less informative than total energy expenditures when estimating the erosion potential or geomorphic effectiveness of discrete discharge events. On the Canadian River, long-duration, moderate-magnitude floods can have larger total energy expenditures than shorter-duration, high-magnitude floods and therefore represent the most serious erosion hazard to floodplain structures.

Environmental & Engineering Geoscience↗
USGS2015

Timing of susceptibility to post-fire debris flows in the western USA

Watersheds recently burned by wildfires can have an increased susceptibility to debris flow, although little is known about how long this susceptibility persists, and how it changes over time. We here use a compilation of 75 debris-flow response and fire-ignition dates, vegetation and bedrock class, rainfall regime, and initiation process from throughout the western U.S. to address these issues. The great majority (85 percent) of debris flows occurred within the first 12 months following wildfire, with 71 percent within the first six months. Seven percent of the debris flows occurred between 1 and 1.5 years after a fire, or during the second rainy season to impact an area. Within the first 1.5 years following fires, all but one of the debris flows initiated through runoff-dominated processes, and debris flows occurred in similar proportions in forested and non-forested landscapes. Geologic materials affected how long debris-flow activity persisted, and the timing of debris flows varied within different rainfall regimes. A second, later period of increased debris flow susceptibility between 2.2 and 10 years after fires is indicated by the remaining 8 percent of events, which occurred primarily in forested terrains and initiated largely through landslide processes. The short time period between fire and debris-flow response within the first 1.5 years after ignition, and the longer-term response between 2.2 and 10 years after fire, demonstrate the necessity of both rapid and long-term reactions by land managers and emergency-response agencies to mitigate hazards from debris flows from recently burned areas in the western U.S.

Environmental & Engineering Geoscience↗