USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Environmental Quality”

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,662 records · Page 50Linked to original sources
USGS1995

The distribution of bromide in water in the Floridan aquifer system, Duval County, northeastern Florida

Although Duval County, Florida, has ample ground-water resources for public supply, the potential exists for a problem with excessive disinfectant by-products. These disinfectant by-products result from the treatment of raw water containing low concentrations of bromide and naturally occurring organic compounds. Because of this potential problem, the relation of bromide concentrations to aquifer tapped, well location and depth, and chemical characteristics of water in the Floridan aquifer system underlying Duval County were studied to determine if these relations could be applied to delineate water with low-bromide concentrations for future supplies. In 1992, water samples from 106 wells that tap the Floridan aquifer system were analyzed for bromide and major dissolved constituents. A comparison of bromide concentrations from the 1992 sampling with data from earlier studies (1979-80) indicates that higher bromide concentrations were detected during the earlier studies. The difference between the old and new data is probably because of a change in analytical methodology in the analysis of samples. Bromide concentrations exceeded the detection limit (0.10 milligrams per liter) in water from 28 of the 106 wells (26 percent) sampled in 1992. The maximum concentration was 0.56 milligrams per liter. There were no relations between bromide and major dissolved constituents, well depth, or aquifer tapped that would be useful for determining bromide concentrations. Areal patterns of bromide occurrence are not clearly defined, but areas with relatively high bromide concentrations tend to be located in a triangular area near the community of Sunbeam, Florida, and along the St. Johns River throughout Duval County.

Florida↗
USGS1984

Techniques for estimating magnitude and frequency of floods on streams in Indiana

Equations are presented for estimating the magnitude and frequency of floods at ungaged sites on unregulated and nonurban streams in Indiana. The equations were developed by multiple-regression, analysis of basin characteristics and peak-flow statistical data from 242 gaged locations in Indiana, Ohio, and Illinois. The State of Indiana was divided into seven areas on the basis of the regression analysis. A set of equations for estimating peak discharges with recurrence intervals of 2, IO, 25, 50, and 100 years was developed for each area. Significant basin characteristics in the equations are drainage area, channel length, channel slope, mean annual precipitation, storage, precipitation intensity, and a runoff coefficient. Standard errors of estimate for the equations range from 24 to 45 percent. Methods are also presented for estimating flood magnitude and frequency at sites on gaged streams. Flood-frequency data based on observed peaks are given for 270 gaged locations. Twenty of these are on regulated streams, and six are on urban streams. Basin characteristics are also included car 245 of the gaged locations on unregulated and nonurban streams. No techniques are given for estimating flood magnitude and frequency at ungaged sites on regulated or urban streams. A rainfall-runoff model was tlsed to synthesize long-term peak data at 11 gaged locations on small streams. Flood-frequency curves developed from the long-term synthetic data were combined with curves based on short-term observed data to provide weighted estimates of flood magnitude and frequency at the rainfall-runoff stations.

Indiana↗
USGS2024

Predicted exposure of communities in southeastern United States to climate-related coastal hazards

A rigorous analysis of 21st Century multi-hazard exposure for U.S. Southeast Atlantic coastal communities indicates that up to 70% of residents will be exposed daily to shallow and emerging groundwater by ~2100, 15 times higher than from surficial flooding alone. This threat further exacerbates other coastal stressors, such as flooding, subsidence, and beach erosion, that impact these vulnerable coastal communities and ecosystems under expected future climate change scenarios.

Florida, Georgia, North Carolina, South Carolina, ↗
USGS2026

Chronic exposure to waterborne nickel significantly reduced growth of juvenile crayfish (Faxonius virilis)

Crayfish are critical functional components of aquatic ecosystems. Previous research has documented adverse effects of mineral extraction on crayfish. Here, we characterize potential risks of mining-derived waterborne nickel (Ni) to crayfish by documenting the effects of dissolved Ni on growth and food consumption of juvenile virile crayfish ( Faxonius virilis) in a 28-day chronic laboratory exposure. Nominal Ni concentrations ranged from 31.25 to 500 micrograms per liter (µg/L; pH = 7.96 ± 0.20, hardness = 150 ± 1 milligrams per liter as calcium carbonate). Crayfish survival, carapace length, and wet weight were measured. After 28 days of exposure, a 24-h feeding trial was performed to determine differences in food consumption. During the growth trial, 99% of crayfish survived. Change in wet weight and final wet weight were the most sensitive endpoints, with 20% effect concentrations of 24.8 and 22.6 µg/L Ni, respectively. Crayfish exposed to an average of 438 µg/L Ni consumed 41% less, and weighed 65.1% less, than control crayfish. These results suggest chronic, sublethal exposure to waterborne Ni may have negative effects on crayfish growth. Reduced growth and consumption rates in crayfish could have wide-ranging consequences throughout aquatic ecosystems since crayfish are consumers, prey, keystone trophic regulators, and ecosystem engineers. Finally, these results could inform bioenergetics and may be coupled with population models to predict potential changes in population sizes of native and invasive crayfishes.

Ecotoxicology↗
USGS2023

Thickness of the Saudi Arabian crust

As part of a joint Saudi Geological Survey (SGS) and U.S. Geological Survey (USGS) project, we analyzed P-wave receiver functions from seismic stations covering most of the Kingdom of Saudi Arabia to map the thickness of the crust across the Arabia Plate. We present an update of crustal-thickness estimates and fill in gaps for the western Arabian Shield and the rifted margin at the Red Sea (the coastal plain), as well as the eastern Arabian Platform. We applied a conventional H-k stacking algorithm and included careful attention to stacking weights, two forms of sedimentary corrections for stations located on the Arabian Platform, and additional processing for noisy stations. We obtained useful results at 154 stations from 898 teleseismic events over a 2-year period from 1995–1997 (for non-SGS stations) and a 6-year period from 2008–2014 (for SGS stations). Average crustal thickness (that is, depth to the Mohorovičić discontinuity [Moho] below the surface) beneath the Red Sea coastal plain (the rift margin) is 29 kilometers (km), beneath the volcanic fields (known in Arabic as harra [plural] or harrat [singular]) is 35 km, beneath the Arabian Shield (excluding harrats) is 37 km, and beneath the Arabian Platform is 38 km. Crustal thinning appears not to extend east of the rift escarpment, suggesting uniform extension that is no broader at depth than at the surface. In contrast to some previous claims that the Arabian Platform crust is thicker than that of the Arabian Shield, we find no statistically significant difference between their whole crustal thicknesses. However, the average subsedimentary crustal thickness (that is, the crystalline crust) for stations on the Arabian Platform is 34 km, 3 km thinner than the crust of the Arabian Shield. Individual station P-wave (pressure) velocity and S-wave (shear) velocity ratios ( V P / V S ) are highly variable for the Arabia Plate, ranging from 1.60 to 1.97 and averaging 1.75, with a standard deviation of 0.07. There are no statistically significant differences between V P / V S ratios of the different geologic regions of Saudi Arabia. Similar V P / V S ratios, coupled with similar crustal thicknesses for harrats and the Arabian Shield, indicate that Cenozoic magmatism has contributed negligibly to crustal growth.

Professional Paper↗
USGS2024

Distribution of ancient carbon in groundwater and soil gas from degradation of petroleum near the Red Hill Bulk Fuel Storage Facility, O‘ahu, Hawai‘i

The groundwater below the Red Hill Bulk Fuel Storage Facility (the facility) in Oʻahu, Hawaiʻi, contains fuel compounds from past spills. This study used carbon-14 analyses to distinguish fuel-derived carbon from background carbon, along with other biodegradation indicators, to address two goals: (1) determine the extent and migration direction of groundwater affected by residual fuel below the facility and (2) determine if residual fuel locations in the subsurface could be identified by analyzing soil gas at the surface above the facility. Groundwater from 19 wells was sampled between September 2022 and April 2023. Nonvolatile dissolved organic carbon (NVDOC) from a well presumed to be unaffected by past spills contained 38 percent ancient carbon indicating a natural source of ancient carbon in the subsurface. The NVDOC concentrations and ancient carbon percentages indicate fuel biodegradation products are likely present on the north and south of Red Hill with the greatest effects at well RHMW02 near the 2014 spill site. The NVDOC concentrations are almost three times higher than diesel range organic (DRO) concentrations in groundwater from the same sites. Major ion data indicate that iron reduction is an important biodegradation process. Soil probe samples and soil carbon traps were used to determine the carbon-14 content of soil carbon dioxide. Ancient carbon from fuel biodegradation was not detected at any soil probe or carbon trap site in contrast to a 2017 study which reported ancient carbon detections. A reanalysis of the 2017 results using a range of local values for background carbon-14 indicates that ancient carbon from fuel biodegradation was probably only detected in lower tunnel exhaust system samples and not in any soil carbon trap samples. Measurements of carbon dioxide efflux with a dynamic closed chamber were highly variable. The soil gas results indicate that soil gas measurements at land surface were not useful for detecting residual fuel at the facility.

Hawaii↗
USGS2026

Morphodynamics of snow barchans at Concordia Station, Antarctica

Barchans are crescent-shaped dunes that form and migrate under unimodal wind regimes in areas of low sediment availability. While sand barchan morphodynamics has been extensively studied, snow barchans remain poorly documented. Here, using hourly time-series images taken near Concordia Station, Antarctica, we measure the sizes (lengths 0.3 m to 13 m, widths 0.4 m to 7.5 m) and migration rates (0.13 m h −1 to 3.5 m h −1 ) of 44 populations of snow barchans during 2017–2022. We compare this dataset with established height–velocity and length-width relationships for sand barchans in both air and water. Our results show that snow barchans follow the same scaling as sand barchans, allowing us to estimate a characteristic sintering (bonding of snow particles) time of about one day. This sintering process limits the size of snow barchans when active, reduces their aspect ratios, and, once wind falls below the transport threshold, explains their preservation over a few years. During wind events, migrating barchans are too small to maintain their shape and orientation, making them highly responsive to fluctuations in wind direction. Thus, studying snow barchans improves understanding of an elementary bedform encountered in multiple environments and also of the surface properties of snow on different timescales, offering insights into snow accumulation and wind-blown snow in ice cores.

Geomorphica↗
USGS2024

Assessment of fresh groundwater discharge and saline surface-water intrusion at Operable Unit 2, North Chevalier Field Disposal Area (Site 11), Naval Air Station Pensacola, Florida, 2018–22

Site 11 is a former landfill at North Chevalier Field Disposal Area in Operable Unit 2 at Naval Air Station Pensacola, in northwest Florida. Site 11 is adjacent to Bayou Grande, a shallow, tidally influenced, saline estuary of the Pensacola Bay watershed. Federal and Florida regulators have expressed concern that contaminants detected in groundwater beneath the inland parts of Site 11 may discharge to Bayou Grande. In 2017, the Department of Defense, U.S. Navy, Naval Facilities Engineering Systems Command Southeast asked the U.S. Geological Survey to assess the occurrence of fresh groundwater discharge to Bayou Grande at Site 11 and to delineate to the extent practicable the location of groundwater discharge. Between 2018 and 2022, the U.S. Geological Survey used a multiple-lines-of-evidence approach that included a visual method and three physical methods based on the temperature difference between groundwater and surface water to assess groundwater discharge. One of the physically based methods also used the difference in specific conductance between fresh groundwater and brackish to saline surface water. Combined, the data indicate that fresh groundwater from across Site 11 discharges primarily along the shoreline of the northern and northeastern part of Site 11. The data also indicate that saline surface water from Bayou Grande intrudes tens of feet into the shallow aquifer beneath Site 11. The combined data indicate that the interface between fresh groundwater and saline surface water changes over space and time. Any new monitoring wells proposed for installation near the shoreline of Site 11 should include approaches to monitor the changes in the location of the freshwater/saltwater interface. Care would need to be taken to collect any groundwater samples at the correct season and tidal period to provide the highest probability of collecting a representative sample of Site 11 groundwater unaffected by saltwater.

Florida↗
USGS2025

Geologic map of the Guinevere Planitia quadrangle (V–30), Venus

This 1:5,000,000-scale geologic map of the Guinevere Planitia quadrangle divides the region into 15 geologic material units, defined using Magellan synthetic aperture radar (12.6-centimeter-wavelength radar system; 75 meters per pixel) datasets and including upland terrain units (2.4 percent of the surface area), plains materials units (59 percent), flow materials associated with named and unnamed eruptive centers (37.2 percent), small volcanic edifices, and impact crater materials (1.4 percent). Upland terrain units consist of tessera and lineated upland material, plains materials consist of Guinevere regional plains and Guinevere lineated and mottled plains, and flow materials consist of lobate flow material and plains-forming flow material. Specific lobate flows associated with Atanua Mons, Tuli Mons, Var Mons, and Uilata Fluctus are mapped separately. Other mapped units are impact crater material and small volcanic edifice. In addition to geologic units, we mapped linear features that show patterns of deformation or flow across the quadrangle. These consist of faults, wrinkle ridges, broad arches, channels, troughs, and flow direction indicators. The map region also contains several small volcanic features: shields, depressions, and craters. These, in combination with the plains, large volcanoes, and coronae, show the pervasive influence of volcanism across Venusian lowlands. The rims of nine identified impact features are delineated; large bright and dark haloes, which in some cases are associated with individual impact craters, are mapped as surficial mantling deposits. We documented spatial relationships using the stratigraphic and cross-cutting relationships of the quadrangle’s geologic units and features to provide a synthesis of the region’s geologic history. The upland terrain of the quadrangle indicates intense tectonic deformation and uplift. It is exposed as embayed remnants, typically within the plains, and represents the oldest geologic materials locally and across the region. Guinevere plains and the plains-forming flow unit appear to be assemblages of volcanic flows from multiple sources, including distinct coronae and corona-like structures. The temporal evolution of Guinevere lineated and mottled plains was likely protracted, with continued formation of small volcanic edifices over a long period. The morphologic and radar brightness characteristics of volcanoes in the region indicate their growth may have involved (1) multiple large-scale eruptive centers with recognizable spatial and temporal sequences, (2) extensive lava flow fields with a multitude of flows producing complex, overlapping patterns, and (3) numerous small volcanic edifices, including shields, domes, and cones. Although geologic patterns common to other regions of Venus are evident in the Guinevere Planitia quadrangle, local relative age relationships are inconsistent or unclear, preventing robust stratigraphic correlation. The mapping results do, however, indicate complicated local sequences of volcanic and tectonic activity.

Scientific Investigations Map↗
USGS2024

Bedrock geologic map of the Woodstock quadrangle, Grafton County, New Hampshire

The bedrock geology of the Woodstock 7.5-minute quadrangle consists of highly deformed metasedimentary rocks of the Central Maine trough, including the Silurian Rangeley and Perry Mountain Formations and the Devonian Littleton Formation. The central, northern, and eastern parts of the quadrangle are underlain by the oldest rocks in the area, the Rangeley Formation. In the southwest and south-central part of the quadrangle, metaturbidites of the Perry Mountain Formation and subsequent Littleton Formation overly the Rangeley Formation in a deformed F 1 synform, herein informally called the Bagley Brook basin. The metasedimentary rocks were intruded by widespread syn- to post-tectonic granitoids of the Devonian New Hampshire Plutonic Suite and minor post-metamorphic Jurassic-Cretaceous mafic dikes of the White Mountain Plutonic-Volcanic Suite. The metasedimentary rocks were affected by at least two episodes of deformation in the Devonian Acadian orogeny. The dominant regional foliation is second-generation (S 2 /D 2 ) and formed during the development of sillimanite-muscovite mineral assemblages. Large bodies of the Early Devonian Kinsman Granodiorite intruded the metasedimentary rocks semi-concordantly during D 2 deformation. Dikes of the Late Devonian Concord Granite cut the Kinsman Granodiorite and the metasedimentary rocks and were emplaced either syn- or post-D 2 . The map pattern in the Rangeley Formation is dominated by northeast to northwest trending, moderately to steeply north-dipping F 2 and F 3 folds. Map-scale F 1 folds are defined by the Bagley Brook basin. Previous division of Rangeley Formation stratigraphy in this region into “upper” and “lower” parts was not corroborated by 1:24,000-scale mapping of lithodemic units, and rocks previously mapped as part of the Smalls Falls and Madrid Formations are here reassigned to the Rangeley Formation. Some rocks previously mapped as the lower part of the Littleton Formation are now assigned to the Perry Mountain Formation. The Littleton Formation on this map is approximately equivalent to rocks previously mapped as the upper part of the same formation. Steeply dipping fractures in the quadrangle show a preferred northeast orientation, consistent with subsurface fracture orientations in the well fields near Mirror Lake. Jurassic-Cretaceous mafic dikes and normal faults show preferred northeast orientations, similar to the fractures, suggesting that the extensional stress field that controlled dike orientation during the Mesozoic also produced the dominant brittle fabrics in the area.

New Hampshire↗
USGS2024

Static and dynamic strain in the 1886 Charleston, South Carolina, earthquake

During the 1886 Mw 7.3 Charleston, South Carolina, earthquake, three railroads emanating from the city were exposed to severe shaking. Expansion joints in segmented railroad tracks are designed to allow railroad infrastructure to withstand a few parts in 10,000 of thermoelastic strain. We show that, in 1886, transient contractions exceeding this limiting value buckled rails, and transient extensions pulled rails apart. Calculated values for dynamic strain in the meizoseismal region are in reasonable agreement with those anticipated from the relation between strain and moment magnitude proposed by Barbour et al. (2021) and exceed estimated tectonic strain released by the earthquake by an order of magnitude. Almost all of the documented disturbances of railroad lines, including evidence for shortening of the rails, can thus be ascribed to the effects of dynamic strain changes, not static strain. Little or no damage to railroads was reported outside the estimated 10 −4 dynamic strain contour. The correspondence between 10 −3 and 2×10 −4 contours of dynamic strain and Mercalli intensity 9 and 8, anticipated from the dependence of each quantity on peak ground velocity, suggests it may be possible to use railroad damage to quantitatively estimate shaking intensity. At one location, near Rantowles, ≈20 km west of Charleston, a photograph of buckled track taken one day after the earthquake has been cited as evidence for shallow dextral slip and has long focused a search for a causal fault in this region. Photogrammetric analysis reveals that the buckle was caused by transient contraction of <10 cm with no dextral offset. Our results further weaken the evidence for faulting in the swamps and forests south of the Ashley River in 1886, hitherto motivated by the photograph and limited macroseismic evidence for high‐intensity shaking.

South Carolina↗
USGS2026

Comparative assessment of STIC sensors, streamflow and rain gauges for quantifying river connectivity in intermittent systems

In intermittent stream systems, including those occurring in Texas, USA, the severity of low-flow conditions, duration of seasonal disconnection, and frequency of no-flow events have been amplified by drought. Documentation of these no-flow events is necessary to evaluate ecosystem health. However, many intermittent reaches remain un-gauged given that perennial river sec-tions are often prioritized for gauge placement. Our objectives were to 1) document stream flow using Stream Temperature, Intermittency, and Conductivity (STIC) loggers to determine the frequency and duration of no-flow events in intermittent tributaries of the Colorado River, Texas and 2) compare logger data to publicly available data from streamflow discharge and precipitation gauge networks to understand differences among these data types for drying event characterization. We use these comparisons to summarize benefits and limitations of the application of in-stream data loggers. STIC loggers were deployed at 19 sites, one in each pool and riffle habitat of a stream reach. STIC loggers recorded a measurement of relative conductance every six hours from June 2022 to March 2024, which was used to determine the presence or absence of flow connectivity in a reach. No-flow duration among intermittent reaches varied between 37 and 270 days across tributaries during an ongoing drought in the study area. Overall, logger data was more precise than discharge data for characterizing no-flow events or precipitation data when documenting presence of water in the stream channel due to runoff. Lack of discharge gauges in intermittent tributaries left large sections of stream reaches undocumented and resulted in mischaracterization of flow patterns. Drought severity across the tributaries did not follow longitudinal patterns that would be expected by the climatic precipitation gradient of the study area. More research is needed to determine if factors such as population size affect severity. Likewise, precipitation data did not correlate well with logger water presence data, lacking consideration for groundwater recharge, soil hydrophobicity, and surface compaction. This study shows that to monitor no-flow events, detailed spatial datasets are necessary and that STIC loggers are useful tools that provide data to fill spatial information gaps and facilitate more accurate flow characterization and water presence data in intermittent systems.

Texas↗
USGS2024

Peri-Gondwanan sediment in the Arkoma Basin derived from the north: The detrital zircon record of a uniquely concentrated non-Laurentian source signal in the late Paleozoic

During the assembly of Pangea, peri-Gondwanan terranes collided with the eastern and southern margins of Laurentia and brought with them unique detrital zircon U-Pb signatures. Discriminating between individual peri-Gondwanan terranes in the detrital record is difficult due to their similar geologic histories. However, characterization of this provenance is critical for understanding late Paleozoic sediment routing during development of Pangea. Along southeastern Laurentia, in the Arkoma Basin (present-day Arkansas and eastern Oklahoma, southeastern United States), we identified Middle Pennsylvanian (Desmoinesian) strata that exhibit a concentrated peri-Gondwanan detrital zircon signature (e.g., ca. 800–550 Ma). Although several southern peri-Gondwanan terranes (e.g., Maya, Suwannee) are closer to the Arkoma Basin, geologic data, such as predominantly north-to-south paleocurrents and proximal-to-distal facies relationships in these Desmoinesian strata, support a northern source (e.g., Ganderia, Avalonia, Meguma). Further evidence of a northern source comes from detrital zircon source mapping, which reveals the persistence of this peri-Gondwanan signal in depocenters to the north of the basin after the signal had diminished in the Arkoma Basin. To this end, bottom-up detrital zircon source modeling, source mapping, regional stratigraphy, paleocurrent data, and sandstone petrography allow us to reconstruct the evolution of this Middle Pennsylvanian (Desmoinesian) sediment pathway in the context of intraplate and plate-margin tectonic activity. This reconstruction documents processes affecting Earth’s surface (e.g., tectonics, climate) during the assembly of Pangea and describes in detail part of a dynamic continental-scale drainage system.

Geosphere↗
USGS2025

Comparison of creek and bay influences on salt marsh sediment budget and deposition patterns

The resilience of salt marshes with low organic production depends on their effective capture and retention of mineral sediment from adjacent waters. Little prior work has directly compared mechanisms of sediment import from wave-influenced marsh boundaries against those of tidal creeks. We used simultaneous deployment of net-deposition tiles and oceanographic sensors to identify the timing and magnitude of sediment import/export to, and redistribution within, a marsh in south San Francisco Bay. As the marsh has both an eroding bay-exposed scarp and a prominent tidal creek, we investigated the mechanisms and magnitudes of sediment import from the marsh-bay versus the marsh-creek interface. The strong daily sea breezes of the summer season produced most of the wave-driven erosion of the marsh scarp and controlled suspended sediment concentrations; the winter season had weaker winds punctuated by a few storms. A large seasonal difference in suspended sediment concentrations influenced both flood and ebb sediment fluxes to the marsh and led to much higher rates of import in the summer. Both bay-side and creek-side processes were important to total marsh sediment budget. Bay-side sediment contributions were more variable in time due to the bay-influenced environment, and creek-side contributions were overall larger, reflecting the large proportion of the marsh fed by creek water. Sediment was redistributed throughout the system, with erosion near the bay-edge, accretion near the creek-edge and slow import to the marsh interior. The marsh was net importing sediment in the summer and exporting in the winter from different rates of these processes; on an annual scale, the marsh was net importing despite rapid lateral marsh loss. These findings emphasize that a positive sediment budget does not imply a stable marsh and that both creek- and edge-side dynamics are important for marsh sedimentation and geomorphic trajectories. Further, we expand understandings of non-storm and seasonal controls on marsh sedimentation.

California↗
USGS2025

Altitude of the potentiometric surface and depth to water in the Mississippi River Valley alluvial aquifer, spring 2022

Potentiometric-surface and depth-to-water maps for spring 2022 were created for the Mississippi River Valley alluvial aquifer (MRVA) using groundwater-altitude data from 1,136 wells completed in the MRVA and from the altitude of the top of the water surface in area rivers from 160 streamgages. The potentiometric-surface and depth-to-water maps for 2022 were created to support investigations to characterize the MRVA as part of the U.S. Geological Survey Water Availability and Use Science Program. Sufficient data were available to map the potentiometric surface and depth to water of the MRVA for spring 2022 for about 83 percent of the aquifer area. The potentiometric contours ranged from 0 to 340 feet (ft) above the North American Vertical Datum of 1988. The regional direction of groundwater gradient was generally to the south-southwest, except in areas of groundwater-altitude depressions, where the groundwater gradient was into the depression, and near rivers, where the groundwater gradient can be from aquifer to the river or from the river into the aquifer. There are large depressions in the potentiometric-surface map in the lower one-half of the Cache region and in much of the Grand Prairie and Delta regions. Depth to water in the MRVA, spring 2022, by well ranged from 5.00 ft above land surface to 145.66 ft below land surface.

Mississippi River Valley alluvial aquifer↗
USGS2024

Effects of stochastically-simulated near-fault ground motions on soil liquefaction

The scarcity of historically recorded near-fault ground motions poses a challenge to systematically understanding the influence of near-fault effects on various types of seismic demands for engineering purposes. In particular, the current state of knowledge of the influence of ground-shaking intensity on soil liquefaction and its consequences does not specifically account for the effects of near-fault ground motion characteristics. In this study, the influence of near-fault ground motions on liquefaction triggering and lateral spreading are investigated using non-linear modeling of a hypothetical liquefiable soil column in the finite-element computational platform OpenSees subjected to simulated ground motion time series that represent strong earthquake shaking in the near field. The simulated ground motion time series and resulting datasets are based on a parametric stochastic model and are developed for a range of source and path parameters to represent a realistic variability of ground motion characteristics. Dependencies between ground motion intensity measures (IMs) and liquefaction demand parameters are investigated for near-fault pulse and nonpulse-like ground motion sets. Evolutionary IMs, such as cumulative absolute velocity (CAV) and the time-varying magnitude-adjusted peak ground acceleration (PGAM), are considered in developing liquefaction triggering probability density functions. Post-liquefaction triggering responses such as lateral spreading displacements are examined in relation to PGA M and CAV. The ground motion simulations are validated by comparing their liquefaction-capacity PGA M fragilities and post-triggering CAV vulnerability relationships to historical records from the 1994 Northridge earthquake in California, USA. Finally, a path forward for future studies that includes finding systematic differences in the IM-liquefaction demand relationships between near-fault and far-field stochastic ground motion sets is outlined.

Conference Paper↗
USGS2026

Spectral characteristics and mapping of lithium-rich playas in the Basin and Range Province, western United States

This report presents a comprehensive set of prospectivity maps for lithium in brines located within the playas of the Basin and Range Province of western United States. The prospectivity assessment is based on criteria established in the descriptive deposit model outlined in a previous U.S. Geological Survey report, with a focus on regional remote sensing methodologies. The playas within the study area have been systematically ranked according to several key factors, including size, the presence and abundance of source rocks, vegetation as an indicator of water availability, documented prospects, and data products derived from the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) sensor. Notably, the deposits in the playa of Clayton Valley make it stand out as the sole lithium brine-producing playa in the United States. The ASTER sensor provides optimal regional coverage, characterized by a 60-kilometer swath width and unique multispectral capabilities that are not available with other remote sensing instruments. This report and the accompanying maps are intended to serve as resources for further multidisciplinary assessments of lithium or other elements or commodities where prospectivity can be inferred from the presence of proxies such as alteration minerals.

Arizona, Baja California, California, Chihuahua, I↗
USGS2025

Long-term dynamics of earthquake swarms in the Yellowstone caldera

The factors controlling the spatial distribution and temporal evolution of earthquake swarms in volcanic systems remain unclear. We leverage leading-edge deep learning algorithms and a detailed three-dimensional velocity model to construct a 15-year high-resolution earthquake catalog of the Yellowstone caldera region. More than half of the region’s earthquakes are clustered into swarm-like families characterized by episodes of hypocenter expansion and migration. Adjacent earthquake swarms, separated by long quiescent periods, are found to be a dominant feature. We suggest that these swarms are controlled by the interplay between slowly diffusing aqueous fluids and rapid episodic fluid injections, which may result from the breaking of permeability seals. Our analyses also indicate that clustered seismicity beneath the caldera occurs on relatively immature, rougher fault structures, compared to more planar faults outside. Our results provide additional context for understanding seismicity in hydrothermal systems, highlighting the key role played by long-term fluid diffusion processes in driving the occurrence of earthquake swarms.

Idaho, Montana, Wyoming↗