USGS Science⌕ Search

SEARCH · USGS Science

Results for “Science of the Total Envionrment”

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.

1,687 records · Page 30Linked to original sources
USGS2025

Hydrologic response of groundwater and streamflow to natural and anthropogenic drivers of change in headwaters of the upper Colorado River basin during recent wet (1982–1999) and drought (2000–2022) conditions

Study region: Headwaters of the upper Colorado River basin (UCOL), USA Study focus: Surface-water and groundwater numerical models incorporating water-use information were used to investigate changes in climate, water use, and simulated hydrologic responses of snow processes, evapotranspiration, groundwater, and streamflow during recent wet (1982–1999) and drought (2000–2022) periods in the headwater subregions of the upper Colorado River basin. New hydrologic insights for the region: Decreases in average streamflow between wet and drought periods ranged from 20 % in the Colorado River headwaters subregion to 23 % in the Gunnison River headwaters subregion. Like streamflow, average surface runoff was statistically less during the drought than the wet period, with decreases from 24–31 % in the headwaters. On a volume basis, runoff decreases were greater than streamflow decreases in both the Colorado River and Gunnison River headwaters. Although the amount of water-year groundwater discharge to streams remained nearly the same between the wet and drought periods, groundwater as a percentage of streamflow increased between the wet and drought periods, highlighting the importance of groundwater in sustaining streamflow during drought conditions. Multiple linear regression analyses revealed that snowmelt-only models were better than the best precipitation and temperature models at explaining streamflow variability from all headwater subregions for both the wet and drought periods.

Arizona, Colorado, New Mexico, Utah, Wyoming↗
USGS2024

Characterizing directivity in small (M 2.4-5) aftershocks of the Ridgecrest sequence

Directivity, or the focusing of energy along the direction of an earthquake rupture, is a common property of earthquakes of all sizes and can cause increased hazard due to azimuthally dependent ground‐motion amplification. For small earthquakes, the effects of directivity are generally less pronounced due to reduced rupture size, yet the directivity in small events can bias source property estimates and provide important insights into general regional faulting patterns. However, due to observational limitations, directivity is usually only measured and modeled for large events. As such, many studies of small earthquakes either ignore directivity altogether or assume a constant rupture direction for all events in a cluster. In our study, we apply a refined directivity fitting method constrained with two separate methods of source deconvolution to the dataset of aftershocks of the 2019 Ridgecrest earthquakes, which contain a large number of well‐recorded small‐to‐mid sized earthquakes occurring in close proximity to each other. The revealed directivity of 100+ small (M 2.4–5) earthquakes is highly heterogeneous and primarily oblique to and away from the main fault strike, suggesting a complex postseismic stress redistribution. In addition, the energy focusing effect of directivity appears to bias the selection of high‐quality data from stations in the direction of rupture, leading to average stress‐drop increases of 50% if directivity is not accounted for.

California↗
USGS2024

Site response in the Walnut Creek–Concord region of the San Francisco Bay, California: Ground motion amplification in a fault-bounded basin

Thirty‐seven portable accelerometers were deployed in the eastern San Francisco Bay communities of Walnut Creek and Concord to study site response in a fault‐bounded, urban, sedimentary basin. Local earthquakes were recorded for a period of two years from 2017 to 2019 resulting in 101 well‐recorded events. Site response is estimated by two methods: the reference site spectral ratio method and a source‐site spectral inversion method. The reference site spectral ratio method allows investigation of the variability of site amplification with source azimuth and frequency. The source‐site spectral inversion method yields the best least‐squares fit to site response for a database of ground‐motion records. Both methods show substantial amplification in the Walnut Creek–Concord basin below 2 Hz indicating strong surface‐wave development. Greater amplification is seen for sources aligned along the long axis of the basin. Inversion using close‐in sources at short distances yields lower amplification at longer periods than the entire data set due to reduced surface‐wave generation for steeper angles of incidence. Inversion of site response spectra for shallow shear‐wave velocity using a global search algorithm yields V S30 values consistent with generalized mapping results based on geology and topography but with greater variability due to local site variations. 3D finite‐element modeling shows greater amplification in the Walnut Creek–Concord basin with a basin‐edge effect likely contributing to higher ground motions. Topography is also seen to lead to increased scattering and shadowing effects.

California↗
USGS2025

Geologic map of the Greater Antilles and Virgin Islands

Introduction The geologic map of the Greater Antilles and Virgin Islands is a compilation of information from the literature, integrated to provide a seamless geologic map of the region. This map was prepared to serve as a base map for a mineral resource assessment of the region. Several small-scale regional geologic maps of the region have been prepared in the past. This report supersedes an earlier version of the geologic map of the Greater Antilles and the Virgin Islands that was released as U.S. Geological Survey (USGS) Open-File Report 2019–1036. For this report, the regional geologic overview shown on sheet 1 covers Cuba; the island of Hispaniola, which includes Haiti and Dominican Republic; Jamaica; the Cayman Islands; Puerto Rico; and the U.S. and British Virgin Islands. Sheet 2 shows the geology of Cuba and Cayman Islands, sheet 3 the geology of Jamaica, sheet 4 the geology of Hispaniola, and sheet 5 the geology of Puerto Rico and Virgin Islands. Accompanying the maps are three tables: table 1 lists the map units (in alphabetical order by map-unit label), their ages, and the countries or territories in which they can be found. These are listed in alphabetical order for ease in locating unit names from the symbols on the maps. Table 2 is a list of sources for the geologic map compilation of Puerto Rico. Table 3 shows formations in the Devil’s Race Course Group of Jamaica.

Greater Antilles, Virgin Islands↗
USGS2025

Hydrogeologic framework and considerations for drilling and grouting of closed-loop geothermal bores in the Erie-Ontario Lowlands and Allegheny Plateau of New York State

The hydrogeologic framework at closed-loop geothermal sites in the Erie-Ontario Lowlands and Allegheny Plateau of central and western New York is the result of the complex interaction of bedrock geology, glacial geology, and groundwater hydrology, and the occurrence of petroleum and gas. Considerations for closed-loop geothermal bore installation include the thickness and character of glacial deposits, bedrock solubility and depth to competent rock, karst development, the distribution of highly permeable zones and their hydraulic heads, and the presence of saline water, gas, and oil. The hydrogeology of the Erie-Ontario Lowlands and Allegheny Plateau poses challenges to closed-loop geothermal bore drilling and casing; managing drill cuttings, discharge water, and gas; and grouting. The potential to encounter severe challenges typically increases with bore depth. This report highlights hydrogeologic considerations for closed-loop geothermal bore installation in New York’s Erie-Ontario Lowlands and Allegheny Plateau to help guide the efficient and safe development of geothermal resources in the regions.

New York↗
USGS1999

Characterization and simulation of the quantity and quality of water in the Highland Lakes, Texas, 1983-92

The Highland Lakes, located in central Texas, are a series of seven reservoirs on the Colorado River (Lake Buchanan, Inks Lake, Lake Lyndon B. Johnson, Lake Marble Falls, Lake Travis, Lake Austin, and Town Lake). The reservoirs provide hydroelectric power for the area. In addition, Lake Austin and Town Lake also provide the public water supply for the Austin metropolitan area. Saline water released from Natural Dam Salt Lake during 1987–89 caused increased concern among water managers that high-salinity water entering the Highland Lakes could result in waterquality problems, necessitating additional treatment of the water. The maximum dissolved solids concentrations for the reservoirs after the saline inflow were about two to three times the average concentrations before the inflow. The maximum concentrations of chloride and sulfate after the inflow were about three to five times the average concentrations before the inflow. The concentrations of dissolved solids, chloride, and sulfate in Lake Buchanan, Inks Lake, Lake Lyndon B. Johnson, and Lake Marble Falls were less than the concentrations of the applicable water-quality standards by the end of 1990. Concentrations of these constituents in Lake Travis, Lake Austin, and Town Lake did not decrease to previous levels, which were less than the concentrations of the applicable waterquality standards, until the end of 1991. Constituent concentrations for Lake Buchanan and Inks Lake; for Lake Lyndon B. Johnson and Lake Marble Falls; and for Lake Travis, Lake Austin, and Town Lake were similar because of the relative storage capacities and location of tributary inflows. From the initial increase in constituent concentrations in Lake Buchanan (summer 1987) in response to the saline inflow, the high-salinity water passed through the entire Highland Lakes in about 3.5 years. A mathematical mass-balance model was used to simulate the input and movement of highsalinity water through the Highland Lakes and to estimate monthly mean concentrations of dissolved solids, chloride, and sulfate for wet, average, and dry hydrologic conditions. The simulated median monthly concentrations during the 10-year simulation period for each reservoir generally are larger for the average condition than for the wet condition and generally are larger for the dry condition than for the average condition. The simulated concentrations of dissolved solids, chloride, and sulfate decreased to levels less than the concentrations of the applicable water-quality standards in about 2 to 5 years after the saline water inflow of 1987–89 was simulated for the three hydrologic conditions. Results from the simulations indicate that saline inflows to the Highland Lakes similar to those of the releases from Natural Dam Salt Lake during 1987–89 are unlikely to cause large increases in future concentrations of dissolved solids, chloride, and sulfate in the Highland Lakes. The results also indicate that high-salinity water will continue to be diluted as it is transported downstream through the Highland Lakes, even during extended dry periods.

Texas↗
USGS1976

Solute transport and modeling of water quality in a small stream

An injection of chloride, sodium, and stable strontium was made at a constant rate for 3 hours into Uvas Creek, Santa Clara County, Calif., to determine the mass transport processes in a small stream. Five observation points were selected within a 610-metre reach of the stream below the injection site. Water samples were collected at the observation points during and immediately after the injection. A mathematical model of the stream was obtained by solving analytically and optimally the one-dimensional mass transport equation of the solutes in the stream. Comparison of field results with a simplified mathematical model indicates the dominance of convection in the behavior of sodium and chloride. The concentration of chloride and sodium can be closely simulated by the model. However, strontium cannot be well represented by the simplified model, which contains a first-order decay-type sink.

California↗
USGS2026

Geologic framework, hydrostratigraphy, and ichnology of the Blanco, Payton, and Rough Hollow 7.5-minute quadrangles, Blanco, Comal, Hays, and Kendall Counties, Texas

During 2023, the U.S. Geological Survey, in cooperation with the Edwards Aquifer Authority, revised an initial characterization completed during 2015–16 of the geologic framework, hydrostratigraphy, and ichnology of the Edwards and Trinity aquifers in the Blanco, Payton, and Rough Hollow 7.5-minute quadrangles in Blanco, Comal, Hays, and Kendall Counties, Texas. The purpose of this report is to present the updated geologic framework, hydrostratigraphy, and ichnology of the Trinity and Edwards Groups in those quadrangles. Rocks exposed in the study area are of the Lower Cretaceous Trinity Group and lower part of the Fort Terrett Formation of the Lower Cretaceous Edwards Group. The faulting and fracturing in the study area are part of the Balcones Fault Zone, an extensional system of faults active during the Paleocene to middle Eocene that generally trends southwest to northeast in south-central Texas. Hydrostratigraphically, the rocks exposed in the study area contain a section of the Edwards (Balcones Fault Zone) aquifer, upper zone of the Trinity aquifer, and middle zone of the Trinity aquifer. In the study area, the only hydrostratigraphic units of the Edwards (Balcones Fault Zone) aquifer remaining are (from top to bottom) VII and VIII, which cap several hills. The mapped hydrostratigraphic units of the upper zone of the Trinity aquifer are (from top to bottom) the cavernous, Camp Bullis, upper evaporite, fossiliferous, and lower evaporite. The mapped hydrostratigraphic units of the middle zone of the Trinity aquifer are (from top to bottom) the Bulverde, Little Blanco, Twin Sisters, Doeppenschmidt, Herff Falls (where present), Rust, Honey Creek, Hensell, and Cow Creek hydrostratigraphic units.

Texas↗
USGS2023

Preliminary geologic map of the Southern Santa Rosa Mountains and Borrego Badlands, San Diego County, Southern California

This investigation delineates the geologic framework of an area of 75 square kilometers (km 2 ) located west of the Salton Sea in southern California (fig. 1, on sheet 1). The study area encompasses the south flank of the Santa Rosa Mountains and the eastern part of the Borrego Badlands (sheet 1). In this study area, regionally important stratigraphic and structural elements collectively inform the late Cenozoic geologic evolution of the Anza-Borrego sector of the Salton Trough province. Critical stratigraphic and structural elements in the map area include the following: The well exposed sequence of late Cenozoic, nonmarine sedimentary rocks that filled the Anza-Borrego subbasin (fig. 1) of the Salton Trough; A tectonic boundary that—in the southern Santa Rosa Mountains—separates the sedimentary strata from underlying crystalline rocks of Peninsular Ranges type. This tectonic boundary, named the West Salton Detachment Fault System by Axen and Fletcher (1998), is projected to underlie all late Cenozoic sedimentary strata in the Anza-Borrego subbasin of the Salton Trough; A variety of transpressional, transtensional, and strike-slip structures that have deformed the late Cenozoic sedimentary strata and collectively guided syntectonic and posttectonic depositional events within the Anza-Borrego subbasin of the Salton Trough; and The southeasternmost surface expression of the Clark Fault, a major strand of the dextral San Jacinto Fault Zone. Geologic mapping and analysis for this investigation focused on clarifying geologic relations among these four stratigraphic and structural aspects in the map area.

California↗
USGS2024

Use of a numerical groundwater-flow model and projected climate scenarios to simulate the effects of future climate conditions on base flow for reach 1 of the Washita River alluvial aquifer and Foss Reservoir storage, western Oklahoma

To better understand the relation between climate variability and future groundwater resources in reach 1 of the Washita River alluvial aquifer and Foss Reservoir in western Oklahoma, the U.S. Geological Survey, in cooperation with the Bureau of Reclamation, used a previously published numerical groundwater-flow model and climate-model data to investigate changes in base flow and reservoir storage by evaluating three scenarios. The three projected climate scenarios were (1) a central-tendency scenario, (2) a warmer/drier scenario, and (3) a less-warm/wetter scenario. To estimate future base flow and groundwater availability in western Oklahoma, specifically in reach 1 of the Washita River alluvial aquifer, downscaled climate-model data from 231 Coupled Model Intercomparison Project phase 5 (CMIP5) projections coupled with a previously published numerical groundwater-flow model were used to compare the effects of different climate scenarios on the aquifer. Changes in base flow and groundwater-level elevations during a 30-year baseline scenario (1985–2014) and the three 30-year projected climate scenarios (2050–79) under central-tendency, warmer/drier, and less-warm/wetter climatic conditions were assessed by using the calibrated model. In the simulations, the amount of base flow and reservoir storage declined in the central-tendency and warmer/drier scenarios compared to the amount of base flow and reservoir storage under historical climatic conditions (baseline scenario). Mean annual change in reservoir storage decreased from the baseline scenario the most in the warmer/drier scenario, followed by the central-tendency scenario, but increased in the less-warm/wetter scenario compared to the baseline scenario. At the end of the simulation period (2079), the largest magnitude differences in groundwater-level elevations in all three projected climate scenarios relative to the baseline scenario occurred upstream from Foss Reservoir. Results from incorporating downscaled climate projections into localized numerical groundwater-flow models can highlight potential future changes in and implications for groundwater resources and availability.

Oklahoma↗
USGS2025

New paleomagnetic constraints on the eruption timing, stratigraphy, and post-emplacement deformation of the Picture Gorge Basalt within the Columbia River Basalt Group

The Picture Gorge Basalt (PGB) is part of the Miocene Columbia River Basalt Group (CRBG). The PGB, which outcrops in eastern Oregon, is considered coincident in time with the voluminous Grande Ronde Basalt. New radiometric ages have expanded the age‐range of the PGB, suggesting it began erupting prior to the Steens Basalt to its south and continued until after cessation of the Grande Ronde Basalt eruptions, an interval of 1.5 Ma. However, the existing paleomagnetism of the PGB implies this eruption timeline is an overestimate. To reconcile the radiometric and paleomagnetic timescales for the PGB, we conducted a paleomagnetic study on sections of the PGB to construct a detailed, high‐quality magnetostratigraphy. Our data indicate the stratigraphically lowest lava flows in the PGB are of reversed polarity, revealing a new paleomagnetic transition with the PGB and a reversed (R)–normal (N)–reversed (R) sequence. This suggests one of two timeline possibilities for PGB volcanism: (a) eruptions began and during through CRBG polarity chrons R0–N0–R1, penecontemporaneous with Steens Basalt, or (b) eruptions began and persisted during CRBG polarity chrons R1–N1–R2. Our work supports a longer interval of PGB volcanism than was suggested by previous paleomagnetic data but is at odds with the suggestion that PGB eruptions lasted through the entire main CRBG. We favor a scenario wherein PGB eruptions begin with R0 and continue into the R1 paleomagnetic interval. The paleomagnetic results also record a ∼18° vertical‐axis rotation of east‐central Oregon after ∼16 Ma with respect to stable North America.

Oregon↗
USGS1991

Reconnaissance of the occurrence of agricultural chemicals in ground water in Haywood, Lake, Obion and Shelby Counties, Tennessee

Data on the occurrence of agricultural chemicals in ground wafer in Tennessee are sparse. The surficial alluvial aquifer is an important source of domestic water supply in West Tennessee, and potentially is subject fo contamination from the application of agricultural chemicals in the area. Nineteen shallow wells completed in the alluvial aquifers in areas of high density agricultural use were sampled in the winter and again in the summer of 1988 to ascertain the occurrence of agricultural chemical in ground water. Although no triazine herbicides or organophosphorus insecticides were detected in any of the wells sampled, elevated nitrite plus nitrate (as nitrogen) concentrations were detected. Results from the winter sampling period indicate a range of nitrite plus nitrate (as nitrogen) concentrations of less than 0.1 to 7.8 milligrams per liter with a median concentration of 2.6 milligrams per liter. Results from the summer sampling period indicate a range of nitrite plus nitrate (as nitrogen) concentrations of less than 0.1 to 8.9 milligrams per liter, median, 2.5 milligrams per liter. The highest concentrations occurred in the shallowest wells, and, in one instance, in a shallow well near a heavily irrigated field.

Tennessee↗
USGS2025

Utilizing downhole datasets for modelling the aeromagnetic signature of the Iron Creek Co-Cu deposit in the Idaho Cobalt Belt

The Idaho Cobalt Belt in east-central Idaho is host to some of the largest domestic Co resources, including the informal Iron Creek deposit. The two main ore zones of this deposit, the Iron Creek and the Ruby, are hosted in greenschist-grade interbedded argillite/siltstone and quartz-rich units of the Mesoproterozoic Apple Creek Formation of Lemhi Group. The primary ore mineral is cobalt-bearing pyrite, which occurs with pyrrhotite, chalcopyrite, and magnetite. This study integrates mineralogical and core-scale geophysical measurements with regional-scale aeromagnetic surveys. The high magnetite content within the Ruby zone produces elevated magnetic susceptibility, but the relatively limited spatial footprint of the ore zone would produce a small-scale anomaly that may be overlooked in regional surveys. The low magnetite content in the Iron Creek zone results in low magnetic susceptibility, creating a relatively low amplitude geophysical response. By characterizing the magnetic properties and mineralogy of these ore zones, this study enhances the interpretation of aeromagnetic data, enabling the identification of small or faint anomalies as potential Co targets. These findings improve can improve exploration strategies, both within the Idaho Cobalt Belt as well as for similar deposit types globally.

Idaho↗
USGS2025

Simulation of the impacts of projected climate change on groundwater resources in the urban, semiarid Yucaipa Valley watershed, southern California using an integrated hydrologic model

Managing water resources in semiarid watersheds is challenging due to limited supply and uncertain future climate conditions. This paper examines the impact of future climate changes on an urban watershed in southern California using an integrated hydrologic model. GSFLOW modeling software is used to simulate the nonlinear relationships between climate trends and precipitation partitioning into ET, runoff, and subsurface storage. Four global circulation models (GCMs), each with two greenhouse-gas scenarios, RCP45 and RCP85 are used to project future climate conditions. GCMs include the CanESM2, CNRM-CM5, HadGEM2-ES, and MIROC5 models. The model's simulated hydrologic conditions are compared with historical data to assess changes in water budgets and groundwater supply. Results indicate decreased groundwater storage in most scenarios due to increased natural evapotranspiration, vegetation consumptive use, and streamflow out of the watershed. Only scenarios with substantially increased future precipitation show increased groundwater storage. The study also highlights increased future aridity despite the rise in precipitation and large precipitation events forecast by GCMs, which increase the risk of urban floods and decrease stream leakage and water available to vegetation.

California↗
USGS2026

Subaerially exposed Iceberg Lake sediments: An exceptional record of historical subaqueous earthquake disturbance at the eastern edge of the Alaskan-Aleutian subduction zone

Paleoseismic records are limited in the Yakutat Terrane (eastern edge of the Alaskan-Aleutian subduction zone) due to the extensive ice cover that hinders traditional methods such as trenching of the faults, but lacustrine sediments offer an alternative archive. We investigated lakebed sediments exposed after recent outburst floods (1999 CE) at Iceberg Lake, a glacier-dammed basin whose stratigraphy was revealed after the lake drained and partially eroded. We logged outcrops across the basin and sampled their sediments. Between annually laminated background deposits, we identified sediment gravity flow beds and in-situ soft-sediment deformation structures (convolutions, sand blows, and fractures) interpreted to be earthquake-induced. Our age model links some of the uppermost depositional and deformational events to the 1958, 1964 and 1979 CE earthquakes. These results demonstrate that Iceberg Lake was a sensitive recorder of seismic shaking and its sediments hold strong potential for producing a paleoseismic record for the northern Yakutat Terrane.

Alaska↗
USGS2025

3D viscoelastic models of slip-deficit rate along the Cascadia subduction zone

Interseismic deformation in the Pacific Northwest is constrained by the horizontal crustal velocity field derived from the Global Positioning System (GPS) in addition to vertical rates derived from GPS, leveling, and tide gauge measurements. Such measurements were folded in to deformation models of fault slip rates as part of the 2023 National Seismic Hazard Model (NSHM) update. Here I build upon one of the contributing models, the viscoelastic earthquake-cycle model of Pollitz [2022]. This model permits inclusion of effects of time-dependent viscoelastic relaxation within earthquake cycles (i.e., ‘ghost transients’) and laterally variable elastic and/or ductile material properties. I lever-age these capabilities to incorporate the Cascadia megathrust into Western U.S.-wide deformation models in which crustal fault slip rates are estimated simultaneously with slip deficit rates along the interplate boundary between the descending Juan de Fuca plate and North American plate. This effort includes construction of a margin-wide model of viscoelastic structure founded on the Slab 2.0 model and probes different models of the ductile properties of the surrounding oceanic asthenosphere, continental lower crust, and mantle asthenosphere. This results in new estimates of the distribution of slip deficit rate along the ∼ 1000 km long margin, highlights the importance of correcting for glacial-isostatic adjustment effects, and permits assessment of sensitivity of results to assumed ductile properties.

California, Oregon, Washington↗