SEARCH · USGS Science

Results for “Environmental Science and Technology”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

1,659 records · Page 48Linked to original sources

USGS2025

Spatiotemporal variations in strain release and seismic rupture in multifault systems: An example from Panamint Valley, southeastern California

Geometrically complex, multifault ruptures have been observed in recent, damaging earthquakes in southeastern California, sparking renewed efforts to identify physical conditions that promote or inhibit fault discontinuity-spanning coseismic ruptures. The likelihood of ruptures propagating across fault discontinuities is thought to be partly controlled by fault geometries, rupture direction, and the history of strain release. However, these parameters vary in space and time over multiple earthquake cycles, making it difficult to forecast the likelihood that an earthquake on one fault will trigger rupture on a nearby fault. Here we use tectono-geomorphic mapping of a geometrically complex fault zone in Panamint Valley, southeastern California, to assess spatiotemporal variations of paleo-rupture patterns and geometries of fault discontinuities over multiple earthquake cycles. First, we identify ten generations of late Pleistocene to Holocene alluvium using geomorphic parameters and luminescence dating to constrain ages of alluvium and bracket late Holocene earthquake timing. Then, we quantify slip kinematics using high-resolution structure from motion digital surface models. We find the Panamint Valley transtensional relay (PVTR) hosted four late Holocene earthquakes, bracketed to ~5.8–3.4 ka, ~3.8–2.2 ka, ~2.4–0.6 ka, and ~0.64–0.16 ka, with ~0.6–1.1 m of slip per event, correlative to Mw ≈ 6.7–6.9 earthquakes. Additionally, we find similarities in earthquake timing on the Ash Hill, PVTR, and Panamint Valley faults and similarities in the slip magnitude and slip kinematics between the Ash Hill and PVTR faults, implying that the PVTR may co-rupture with nearby faults. Paleo-rupture patterns indicate that seismogenic strain transfer may occur through the PVTR, along different combinations of fault segments and jump distances, over multiple earthquake cycles. These data highlight the utility of tectono-geomorphic mapping in evaluating paleo-rupture patterns and suggest that the PVTR may act to propagate and/or arrest rupture between the Ash Hill and Panamint Valley faults.

California↗

USGS2025

Capturing the uncertainty of seismicity observations in earthquake rate estimates: Implications for probabilistic seismic hazard analysis and the USGS National Seismic Hazard Model

The rate of earthquakes in a region is a fundamental input to Probabilistic Seismic Hazard Analysis. We present a Monte Carlo method for computing that rate from seismicity catalogs while including a range of data and analysis uncertainties. This method is applied to regions for which the b value is assumed to be spatially invariant. Each region is broken down into epochs for which each epoch is estimated to have a uniform magnitude of completeness ( ⁠⁠ M c ). The distribution of earthquake rates for M ≥ M c is determined for each epoch by considering the Poisson likelihood of rates given the number of observed earthquakes with M ≥ M c ⁠ . We use a Monte Carlo process to include the uncertainty in b , ⁠ M c , and individual event magnitudes. The result for each epoch is the joint distribution of the Poisson rate of earthquakes with magnitudes larger than the minimum value used to calculate hazard ( ⁠⁠ M 1 ) and the Gutenberg–Richter b values, which control the extrapolation to other magnitudes. The rate for each region is either the duration‐weighted average over the epochs or, to better capture temporal variations, we also consider mixture models. The mixture models also provide an avenue to allow temporal variations in b values. To implement this joint distribution in a logic tree, we use the mean and 95% confidence branches, each of which is parameterized with an M ≥ M 1 rate and b value. We explore different ways of defining those branches, as well as non‐Gutenberg–Richter branches, and their impact on hazard estimates. The mean hazard, but not the fractiles, is robust with respect to these choices. To illustrate these new methods, we use synthetic data and catalogs from recent U.S. Geological Survey National Seismic Hazard Models for the Central and Eastern United States and for Puerto Rico and the U.S. Virgin Islands.

central and eastern United States, Puerto Rico, U.↗

USGS2026

Storage timescales and the crystal growth history recorded in the layered reservoir of the Unit 3 Kulanaokuaiki Tephra Member of the Uēkahuna Ash from Kīlauea volcano

The Kulanaokuaiki Tephra Unit 3 (900 C.E.) was a subplinian basaltic scoria eruption, the largest known to have occurred at Kīlauea’s summit. The initiation mechanism and cause(s) of this highly hazardous eruption remain poorly understood, particularly as this style of activity has not been observed in the modern period. We use olivine diffusion timescales from Fe-Mg and lithium zoning to assess magma storage timescales and evaluate evidence for mixing prior to eruption. Of the 200 crystals assessed, 55% are normally zoned for forsterite and 45% are non-zoned. Zoned crystals (> 1 mm) ( n = 96) have average compositions of Fo 89 cores and Fo 81 rims, while small crystals 0.5–1 mm ( n = 14) have Fo 85 cores and Fo 82 rims. Modeled Fo diffusion timescales range from < 1 to 12.6 years (median of 3.0 years). Lithium maps in 32 of the forsterite-zoned crystals reveal a variety of lithium zoning types, most of which we attribute to crystal growth. Many lithium profiles are coupled to phosphorus within the outer ~ 100 μm of the crystal ( n = 17), reflecting growth processes and precluding their use for diffusion studies. Overall, our Unit 3 olivine dataset displays no evidence of magmatic priming shortly prior to eruption, unlike prior observations for the basaltic subplinian eruption of the Keanakāko‘i Tephra Unit E (1650 C.E.). Thus, explosive eruption magmas may be stored for considerably longer than previously thought, and magmatic priming is not a necessary precursor for Kīlauea’s explosive eruptions. We hypothesize that gas-driven filter pressing from crystallization of microlites and macrocrysts in the reservoir slurry contributed to the explosivity of the Unit 3 eruption.

Hawaii↗

USGS2024

New 10Be-26Al isochron burial dating informs the Pliocene and Pleistocene evolution of the lower Colorado River, southwestern United States

Four new 10 Be- 26 Al isochron burial ages ranging from 4.4 to 2 Ma on ancestral Colorado River deposits in the lower Colorado River corridor (LCRC) help constrain the river’s evolution during the Pliocene and early Pleistocene. They help fill a gap between previous work that focused on older and younger deposits: Older dated deposits include the 5 Ma Bouse Formation, which records the integration of the Colorado River through a series of preexisting basins to the Gulf of California and the ca. 4.5 - 3.5 Ma Bullhead Alluvium, a 200 to 300 m thick aggregational package that immediately followed integration. The much younger, 100 - 70 ka, Chemehuevi Formation is another major aggradation package mapped throughout the LCRC. The new burial ages on the facies of Santa Fe Railway (4.37 ± 0.71 Ma), boulder conglomerate of Bat Cave Wash (2.12 ± 0.26 and 2.05 ± 0.31 Ma), and the Palo Verde alluvium (3.03 ± 0.26 Ma) partially fill in a 3.5 M.y. gap between the deposition of the Bullhead Alluvium and the Chemeheuvi Formation and document the timescales over which the Colorado River was able to remove the Bullhead aggradational package and initiate newer and smaller aggradational pulses.

Arizona, California, Nevada↗

USGS2024

Field observations and logs from the Rose Hip trench exposure across a north-facing scarp within the Seattle Fault Zone, southern Bainbridge Island, Washington

The Seattle Fault Zone is an approximately 70-km-long, east-west-trending zone of south-dipping blind reverse faults within the Puget lowland region in Washington. Because of the proximity, the Seattle Fault Zone poses a significant earthquake hazard to the Puget sound and Seattle metropolitan regions. We present preliminary mapping and trench-site information from a paleoseismic investigation across a newly identified active fault scarp located within the hanging wall of the Seattle Fault Zone on southern Bainbridge Island, Washington. The trench exposed monoclinally folded Miocene bedrock, fractured and faulted glacial-related deposits, and laminated lacustrine deposits capped by slope-derived colluvium. The observations from this investigation record late Pleistocene to Holocene north-vergent folding and faulting along this new fault scarp.

Washington↗

USGS2025

Geochemistry and radiogenic isotopes constrain the mantle source region of the Mountain Pass Intrusive Suite, California

The Mountain Pass carbonatite stock is the largest rare earth element (REE) deposit and only active REE mine in the United States. The carbonatite intrusion and spatially associated alkaline silicate intrusions constitute the Mountain Pass Intrusive Suite, which is located within the Mojave Province in California. Both the carbonatite and the alkaline silicate rocks are enriched in large ion lithophile elements and light REEs, and less enriched to depleted in high field strength elements, indicating the mantle source region was metasomatically enriched in incompatible trace elements. The cause of this metasomatic mantle enrichment and the genetic relationship between the carbonatite and the alkaline silicate stocks are poorly understood. In this study, major and trace element geochemical data and isotopic (Rb-Sr, Sm-Nd, and Lu-Hf) data are presented to constrain genesis of the Mountain Pass Intrusive Suite, from mantle source region to the intrusion of the stocks. Our geochemical data are consistent with derivation of the alkaline silicate and carbonatite melts through partial melting from a shared mantle source region rather than through liquid immiscibility or fractional crystallization and separation of a carbothermal fluid. Although the Rb-Sr isotopic system in the Mountain Pass Intrusive Suite is disturbed at the whole-rock scale, the isotopic systems for whole-rock Sm-Nd (εNd i = ‐2.2 ± 0.8) and zircon Lu-Hf (εHf i = 0.1 ± 1.1) are robust and support mantle derivation of the magmas. Geochemical modeling using experimentally derived partition coefficients was used to identify possible causes of enrichment in incompatible elements through metasomatism in the mantle source region. Modeling of metasomatism by melts derived by partial melting of deeply subducted carbonated sediments approximates observed Mountain Pass Intrusive Suite trace element chemistry. Scattered εHf i in inherited zircon (2.8 ± 2.6) is consistent with derivation from an arc-related environment with substantial crustal contamination. Paleotectonic studies in the Mojave Province indicate that regional subduction preceded emplacement of the Mountain Pass Intrusive Suite by ∼300 Ma. Melting of the Mountain Pass source region may have been caused by post-collisional thermal relaxation and extension.

California↗

USGS2025

Conceptual and numerical groundwater flow model of the Iowa River alluvial aquifer near Tama County, Iowa, 1980 through 2022

The Iowa River alluvial aquifer is an important source of water on the Meskwaki Settlement in Tama County, Iowa, which is land owned by the Sac & Fox Tribe of the Mississippi in Iowa (commonly known as the Meskwaki Nation). The U.S. Geological Survey constructed a groundwater flow model, including a conceptual and numerical model, of the Iowa River alluvial aquifer and underlying hydrogeologic units near the Meskwaki Settlement in Tama County, Iowa, for the period of January 1980–August 2022 to estimate the fraction of water pumped from the Iowa River alluvial aquifer by Meskwaki Settlement wells that is derived from streamflow depletion in the Iowa River and its tributaries. Streamflow depletion is a reduction in streamflow caused by groundwater pumping and includes the interception by groundwater production wells of water that otherwise would have been discharged to streams (called “captured groundwater discharge”) and induced infiltration of streamflow to the production wells. Calibrated model runs were performed with no simulated pumping and simulated pumping only at Meskwaki Settlement wells, and the change in simulated flow rates between the groundwater system and streams for the two model runs represents the amount of streamflow depletion in the Iowa River and tributary streams resulting from pumping at the Meskwaki Settlement wells. Streamflow depletion in the Iowa River and its tributaries as a percentage of simulated pumping at the Meskwaki Settlement wells was calculated by dividing this difference by the total simulated pumping rate for the Meskwaki Settlement wells. The model results demonstrate that the mean monthly streamflow depletion, including induced infiltration and captured discharge, in the Iowa River and its tributary streams as a percentage of mean monthly pumping at the Meskwaki Settlement wells was 97.4 percent and ranged from 65.4 to 112 percent. Of the total streamflow depletion, mean monthly induced recharge was 20.9 percent and ranged from 4.9 to 37.2 percent. Mean monthly captured discharge was 76.5 percent and ranged from 57.1 to 97.1 percent. These results indicate that most of the water pumped from the Meskwaki Settlement wells is the result of streamflow depletion, in the form of both induced infiltration and captured discharge.

Iowa↗

USGS2025

Mechanics and statistics of postseismic shaking

Analysis of two weeks of continuous post-seismic shaking after the 2019 M7.1 Ridgecrest, CA earthquake sequence using 4 nearby borehole seismometers reveals that continuous ground motions decay as Omori’s law in time and follow the Gutenberg-Richter distribution in logarithmic amplitude. The measured temporal decay in amplitudes agrees with predictions of the rate-and-state framework and indicates shaking amplitudes are proportional to the velocity of afterslip. Our ground motion-based statistical framework provides a basis to forecast shaking intensity in the minutes to hours after a large earthquake.

California↗

USGS2024

Near-real-time earthquake-induced fatality estimation using crowdsourced data and few-shot large-language models

When a damaging earthquake occurs, immediate information about casualties (e.g., fatalities and injuries) is critical for time-sensitive decision-making by emergency response and aid agencies in the first hours and days. Systems such as the Prompt Assessment of Global Earthquakes for Response (PAGER) by the U.S. Geological Survey (USGS) were developed to provide a forecast of such impacts within about 30 min of any significant earthquake globally. However, existing disaster-induced human loss estimation systems often rely on early casualty reports manually retrieved from global traditional media, which are labor-intensive, time-consuming, and have significant time latencies. Recent approaches use keyword matching and topic modeling to identify human casualty-relevant information from social media but tend to be error-prone when dealing with complex semantics in multi-lingual text data and parsing dynamically changing and conflicting human death and injury numbers shared by various unvetted sources in social media platforms. In this work, we introduce an end-to-end framework to significantly improve the timeliness and accuracy of global earthquake-induced human loss forecasting using multi-lingual, crowdsourced social media. Our framework integrates (i) a hierarchical casualty extraction model built upon large language models, prompt design, and few-shot learning to retrieve quantitative human loss claims from social media, (ii) a physical constraint-aware, dynamic-truth discovery model that discovers the truthful human loss from massive noisy and potentially conflicting human loss claims, and (iii) a Bayesian updating loss projection model that dynamically updates the final loss estimation using discovered truths. We test the framework in real-time on a series of global earthquake events in 2021 and 2022 and show that our framework effectively automates the retrieval of casualty information faster but with comparable accuracy to those now retrieved manually by the USGS. The code associated with this work is made available at: https://github.com/SusuXu-s-Lab/Hierarchical-Earthquake-Casualty-Information-Retrieval

International Journal of Disaster Risk Reduction↗

USGS1976

Uranium and trace elements in stream sediments as an exploration tool

More than 45 trace elements have been reported in anomalous amounts in the uranium ores. The specific suite of elements associated with any one uranium deposit varies according to deposit type and geologic province. The primary geochemical halo of uranium and associated trace elements in the host rock, together with secondary dispersion halo in soils and alluvium, offers a potential geochemical exploration target. Sediment from streams near low-grade uranium occurrences in arkosic sandstone of the Denver basin, Colorado, and sediment from streams in the igneous and metamorphic terrace near the Midnite mine of eastern Washington have been analyzed for acid-extractable U, Cu,Pb, Zn, Ni, Co, Fe, Mn, Mo, Cd, Cr, Ag, V, and Se, and for total Fe, Mn, Hg, As, and organic C. The -80 mesh fraction of the alluvium downstream from a small uraniferous limonite occurrence in the Denver basin shows anomalous concentrations of U, Pb, As, Ni, Mo, V, Zn, and Se. Anomalous concentrations of U persist farther downstream (as much as 1,000 m) than any of the trace elements. In another stream system southeast of Denver, above-background levels of U, Pb, As, Zn, V, Fe, and Mn occur near a uranium anomaly (defined by other methods) which has no surface mineralization. In the Midnite mine area, sediments of streams draining prospect areas showed anomalous concentrations of uranium (as much as 10 times background), but no accompanying anomalous trace-element concentrations, although analyses of Midnite mine ores show that a broad suite of trace elements are associated with uranium.

Open-File Report↗

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↗

USGS2000

Use of a ground-penetrating radar system to detect pre- and post-flood scour at selected bridge sites in New Hampshire, 1996-98

Ground-penetrating radar was used to measure the depth and extent of existing and infilled scour holes and previous scour surfaces at seven bridges in New Hampshire from April 1996 to November 1998. Ground-penetrating-radar survey techniques initially were used by the U.S. Geological Survey to study streambed scour at 30 bridges. Sixteen of the 30 bridges were re-surveyed where floods exceeded a 2-year recurrence interval. A 300-megahertz signal was used in the ground-penetrating radar system that penetrated through depths as great as 20 feet of water and as great as 32 feet of streambed materials. Existing scour-hole dimensions, infilled thickness, previous scour surfaces, and streambed materials were detected using ground-penetrating radar. Depths to riprap materials and pier footings were identified and verified with bridge plans. Post data-collection-processing techniques were applied to assist in the interpretation of the data, and the processed data were displayed and printed as line plots. Processing included distance normalization, migration, and filtering but processing was kept to a minimum and some interference from multiple reflections was left in the record. Of the 16 post-flood bridges, 22 ground-penetrating-radar cross sections at 7 bridges were compared and presented in this report. Existing scour holes were detected during 1996 (pre-flood) data collection in nine cross sections where scour depths ranged from 1 to 3 feet. New scour holes were detected during 1998 (post-flood) data collection in four cross sections where scour depths were as great as 4 feet deep. Infilled scour holes were detected in seven cross sections, where depths of infilling ranged from less than 1 to 4 feet. Depth of infilling by means of steel rod and hammer was difficult to verify in the field because of cobble and boulder streambeds or deep water. Previous scour surfaces in streambed materials were identified in 15 cross sections and the depths to these surfaces ranged from 1 to 10 feet below the streambed. Riprap materials or pier footings were identified in all cross sections. Calculated record depths generally agree with bridge plans. Pier footings were exposed at two bridges and steel pile was exposed at one bridge. Exposures were verified by field observations.

New Hampshire↗

USGS2026

Metallurgical coal—Deposits, production, resources, market dynamics, and supply chain risks

Plain Language Summary Metallurgical coal (met coal; consumed to produce coke for steelmaking) must meet specific chemical and physical specifications. In 2023, the conterminous United States produced 66 million short tons (mst) of met coal, consumed 15.85 mst domestically, exported 51.1 mst, and imported 0.7 mst. Most met coal was produced in the Appalachian Basin, but there are also resources that meet the specifications for met coal in the Western United States.

conterminous United States↗

USGS2025

Subduction zone earthquake catalog separation tool: Implementation in the USGS 2025 Puerto Rico and U.S. Virgin Islands National Seismic Hazard Model

The U.S. Geological Survey (USGS) periodically releases updates to National Seismic Hazard Model (NSHM) for the United States and its territories leveraging current scientific knowledge and methodologies to guide public policy, building codes, and risk assessments regarding potential ground shaking due to earthquakes that may result in infrastructure damage. In subduction zones, there is a need to separate the earthquake catalog into tectonic regimes to create specific seismicity models for which the most appropriate ground‐motion models are then applied. Here, we describe newly developed methods and software, called CatSep, that classifies subduction zone events into three primary tectonic regimes: crustal, interface, and intraslab. This method incorporates information about the location of the earthquake relative to the subducting slab, the depth of the Mohorovičić discontinuity, and the earthquake’s moment tensor. Applying this method is a first step in the NSHM workflow for regions covering U.S. subduction zones. Results using this subduction zone earthquake catalog separation tool for the 2025 Puerto Rico and U.S. Virgin Islands NSHM earthquake catalog are presented and analyzed.

Puerto Rico, U.S. Virgin Islands↗

USGS2026

Interplay between tectonics and submarine mass transport deposits in Cortes Basin: New high-resolution geophysics in the Outer California Borderland

The Outer California Borderland (OCB) is an active transform plate boundary offshore Southern California, where the relationship between faulting and submarine mass transport deposits (MTDs) remains poorly understood. Onshore paleoseismic data provide high-resolution earthquake records, whereas marine geophysical data capture longer-term histories. Offshore fault systems pose hazards to infrastructure and dense coastal populations, particularly when linked to submarine landslides. We present new high-resolution geophysical data set (cruise SR2303), including bathymetric and CHIRP sub-bottom data integrated with legacy seismic reflection data and chronostratigraphic constraints from ODP Site 1012 to examine Quaternary MTD recurrence and tectonic controls in the Cortes Basin, OCB. Bathymetry shows deformational features, including slide scarps and previously unmapped fault segments with evidence of Holocene activity. CHIRP profiles reveal 10 stacked MTDs in the East Cortes Basin and 8 in the West Cortes Basin, spanning ∼752 ka with an average recurrence of ∼83.6 ± 1 ka. Acoustic imaging shows 7 MTD intervals coinciding with fault offset increments and fault growth suggesting earthquake-triggered mass wasting. A strong association between MTD occurrences and sea-level extremes also supports glacio-eustatic contribution to slope failure. Stratigraphic correlations suggest quasi-synchronous MTDs across the eastern and western areas, likely triggered by larger eathquakes in the Quaternary. Although the identified MTDs occur relatively far from the Southern California coast, they still pose a potential tsunamigenic hazard requiring further assessment. Moreover, if linked to earthquakes along major strike-slip faults, for example, the Ferrelo fault, the MTDs may provide valuable proxies to constrain rupture scenarios and fault connectivity within the understudied OCB.

California↗

USGS2025

Petrogenesis and mineralization potential of the Bradley Peak komatiitic basalts, Wyoming Province

The ca. 2.72 Ga Bradley Peak greenstone terrane in the Wyoming Province contains spinifex- and cumulate-textured komatiitic volcanic rocks that may host Ni-Cu-PGE sulfide mineralization. New whole rock geochemistry classifies these rocks as Al-undepleted komatiitic basalts derived from a parental melt with ~17 wt. % MgO. Isotopic data (Sm-Nd, Re-Os) and REE profiles suggest a plume source and possible assimilation of juvenile footwall volcanic rocks. Calculated sulfur concentrations at sulfide saturation and low-detection PGE geochemistry indicate these magmas did not saturate a sulfide liquid resulting in sulfide mineralization. Although the Bradley Peak komatiitic basalt flows are unlikely to host magmatic sulfide mineralization, they are likely the source of Au for regional orogenic gold deposits.

Wyoming↗

USGS2023

Paleomagnetism of the Harrat Rahat volcanic field, Kingdom of Saudi Arabia—Geologic unit correlations and geomagnetic cryptochron identifications

Paleomagnetic rock samples were collected from 173 drill sites in the Quaternary alkali basaltic volcanic field of northern Harrat Rahat, Kingdom of Saudi Arabia. Laboratory measurements on these samples established that lava flows and vent complexes—identified and mapped from field characteristics, rock types, and compositions as products of single or temporally close eruptions—typically record single, or very similar, directions of remanent magnetization. Correlations defined through geologic mapping, spatial association, geochemistry, geochronology, and identical mean remanent directions indicate at least 16 brief episodes of temporally clustered eruptions. These episodes had durations of a few centuries or less. Anomalous remanent magnetic directions were found for at least 13 mapped lavas of northern Harrat Rahat, which demonstrate that they were acquired during brief geomagnetic cryptochrons during the Brunhes Normal Polarity Chron. These uncharacteristic directions enhance the opportunity to identify common eruptive episodes, and to better understand and evaluate assessments of eruption ages based on 40 Ar/ 39 Ar geochronology. Combining paleomagnetic and regional archaeomagnetic results for the youngest eruptions allows us to evaluate their historical age assignments and, in one case, refute a previously assigned provisional age.

Harrat Rahat volcanic field↗

USGS2024

Beyond the wedge: Impact of tidal streams on salinization of groundwater in a coastal aquifer stressed by pumping and sea-level rise

Saltwater intrusion (SWI) is a well-studied phenomenon that threatens the freshwater supplies of coastal communities around the world. The development and advancement of numerical models has led to improved assessment of the risk of salinization. However, these studies often fail to include the impact of surface waters as potential sources of aquifer salinity and how they may impact SWI. Based on field-collected data, we developed a regional, variable-density groundwater model using SEAWAT for east Dover, Delaware. In this location, major users of groundwater from the surficial aquifer are the City of Dover and irrigation for agriculture. Our model includes salinized marshland and tidal streams, along with irrigation and municipal pumping wells. Model scenarios were run for 100 years and included changes in pumping rates and sea-level rise (SLR). We examined how these drivers of SWI affect the extent and location of salinization in the surficial aquifer by evaluating differences in chloride concentration near surface waters and the subsurface freshwater-saltwater interface. We found the presence of the marsh inverts the typical freshwater-saltwater wedge interface and that the edge of the interface did not migrate farther inland. Additionally, we found that tidal streams are the dominant pathways of SWI at our site with salinization from streams being exacerbated by SLR. Our results also show that spatial distribution of pumping affects both the magnitude and extent of salinization, with an increase in concentrated pumping leading to more intensive salinization than a more widely distributed increase of the same total pumping volume.

Delaware↗