USGS Science⌕ Search

SEARCH · USGS Science

Results for “International Journal of Applied Earth Observations and Geoinformation”

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,689 records · Page 35Linked to original sources
USGS2026

New measurements indicate that natural geologic methane emissions from microseepage in the Michigan Basin are likely negligible

The magnitude of natural geologic methane (CH 4 ) emissions to the atmosphere (including emissions of fossil CH 4 from offshore and onshore gas and oil seeps, diffuse microseepage, mud volcanoes, volcanic vents, and geothermal areas) is highly uncertain. The largest component of geologic emissions is thought to be microseepage, which is the diffuse flux of CH 4 from soils across large areas of productive hydrocarbon basins. The accuracy of existing bottom-up estimates of microseepage is limited by low spatial coverage of published microseepage measurements. We present the first soil—atmosphere CH 4 flux measurements from Michigan Basin, USA. Results from 335 measurements taken during summer and winter seasons across a large portion of the basin suggest microseepage is nonexistent in the sampled region. Even areas with predictive features for microseepage (e.g., underlying mature, organically rich source rocks, proven gas accumulations, faults, and lineaments) yield null or negative fluxes, suggesting that CH 4 emissions from microseepage are negligible throughout our study region. A Monte Carlo method was used to place an upper limit on the regional-mean microseepage flux, in which synthetic patchy microseepage distributions were generated and tested against our observations to assess the impact of possible emission hot spots that were missed by our sampling strategy. Our analysis finds it is very unlikely that regional-mean emissions are as high as assumed in a previous global microseepage study. The observed lack of seepage may be explained by groundwater flow, active methanotrophy, glacial sediments, and bedded salt deposits, which could inhibit vertical gas migration and release to atmosphere.

Michigan↗
USGS2024

Extending the Boore and Abrahamson (2023) modified square-root-impedance method for the development of site amplifications consistent with the full-resonance approach to a range of VS30 values

The square-root-impedance (SRI) method is commonly used to approximate the seismic site amplifications computed using the full-resonance (FR) method for gradient shear-wave velocity ( V S ) profiles that are smoothly varying with depth. The SRI site amplifications have been observed to systematically underpredict the FR site amplifications by a ratio of FR/SRI amplifications around 1.05 to 1.3 across a wide frequency range (Boore, 2013). Recently, Boore and Abrahamson (2023; hereafter, BA23) related this difference in the SRI and FR methods to differences in the exponent η of the ratio of seismic impedances between the two methods. They proposed the implementation of a modified frequency-dependent η in the SRI method to improve its match to the FR site amplifications. This modified η was derived using only five V S profiles. We investigate the performance of the BA23 η for a wide range of realistic gradient V S profiles with V S30 ranging from 180 to 1500 m/s. These gradient V S profiles are constructed using two power-law functions of depth and are constrained by the assigned VS30 value, the depth and velocity of the half-space, and depths to shear-wave velocity horizons of 1.0 and 2.5 km/s ( Z 1.0 and Z 2.5 ) based on western United States sites. Despite observing a V S30 dependence of η, we find that the BA23 η generally works reasonably well for the range of V S profiles analyzed. Using the VS30 -dependent η derived in this study results in improvements in matching the FR site amplification compared to using the BA23 η. These improvements are more pronounced for the soft-site conditions and become modest to negligible for the stiff site conditions

Bulletin of the Seismological Society of America↗
USGS2024

Post-glacial stratigraphy and late Holocene record of great Cascadia earthquakes in Ozette Lake, Washington, USA

Ozette Lake is an ~100-m-deep coastal lake located along the outer coast of the Olympic Peninsula (Washington, USA); it is situated above the locked portion of the northern Cascadia megathrust but also relatively isolated from active crustal faults and intraslab earthquakes. Here we present a suite of geophysical and geological evidence for earthquake-triggered mass transport deposits (MTDs) and related turbidite deposition in Ozette Lake since ca. 14 ka. Comprehensive high-resolution bathymetry data, seismic reflection profiles, and sediment cores are used to characterize the post-glacial stratigraphic framework and examine paleoseismic evidence in the lacustrine sediments. Stacked sequences of MTDs along the steep eastern flanks of the lake appear to grade basin-ward from thick, chaotic, blocky masses to thin, parallel-bedded turbidite beds. The discrete turbidite event layers are separated by fine-grained (silt and clay) lake sedimentation. The event layers are observed throughout the lake, but the physical characteristics of the deposits vary considerably depending on proximity to primary depocenters, steep slopes, and subaqueous deltas. A total of 30–34 event deposits are observed in the post-glacial record. Radiometric dating was used to reconstruct a detailed sedimentation history over the last ~5.5 k.y., develop an age model, and estimate the recurrence (365–405 yr) for the most recent 12 event layers. Based on sedimentological characteristics, temporal overlap with other regional paleoseismic chronologies, and recurrence estimates, at least 10 of the dated event layers appear to be sourced from slope failures triggered by intense shaking during megathrust ruptures; the recurrence interval for these 10 events is 440–560 yr. Thus, Ozette Lake contains one of the longest and most robust geological records of repeated shaking along the northern Cascadia subduction zone.

Washington↗
USGS2024

A global assessment of SAOCOM-1 L-band stripmap data for InSAR characterization of volcanic, tectonic, cryospheric, and anthropogenic deformation

SAOCOM-1 is an L-band (23.5 cm) synthetic aperture radar (SAR) constellation made up of two satellites launched in 2018 and 2020 by Comisión Nacional de Actividades Espaciales (CONAE, Argentina). In this contribution, we present a global summary of interferometric SAR (InSAR) observations of ground deformation with SAOCOM-1 stripmap data for tracking volcanic, tectonic, glacier, and anthropogenic deformation. These examples include: 1) episodes of unrest at volcanoes in the Aleutian Islands, Southern Andes, and Italy, with line-of-sight (LOS) deformation from 4 cm/yr in InSAR time series to ~70 cm in interferograms; 2) dike intrusions in Hawai’i; 3) earthquakes in the Andean fold and thrust belt and the East Anatolian fault; 4) ice flow of the Southern Patagonia icefield; and 5) subsidence due to lithium brine extraction in the Salar de Atacama basin (northern Chile). Comparisons between SAOCOM-1, ALOS-2 SM3, Sentinel-1, and TerraSAR-X/ TanDEM-X/PAZ (TSX/TDX/PAZ) mean velocities from InSAR time series show a 1:1 ± 3% correlation in the LOS velocity, which highlights the high accuracy of SAOCOM-1 data. The minimum deformation that we measured in individual interferograms is 4 ± 0.6 cm. One limitation of SAOCOM-1 is the lack of a global acquisition program, which reduces its global and broader applications. Considering the repeat periods, background observation program, and lack of a controlled orbital tube, the best suited targets for SAOCOM-1 InSAR are two. First, volcanoes that deform with secular rates located in vegetated regions in mid- and high-latitudes, and/or that undergo transient episodes of fast deformation in which C-band coherence is lost quickly. Second, glaciers where coherence can be sustained during the repeat period of eight days.

IEEE Transactions on Geoscience and Remote Sensing↗
USGS2025

Temporal and spatial comparison of coal mine ventilation methane emissions and mitigation quantified using PRISMA satellite data and on-site measurements

Emission monitoring at the facility level (bottom-up, BU) is key for accurate reporting of coal mine methane (CMM) emissions. Recent advances in aerial and satellite observations (top-down, TD) indicate that these methods have the potential to support CMM emissions monitoring and reporting of, as well as track the effectiveness of, mitigation actions. That said, studies have shown discrepancies between BU and TD estimations. Performing TD monitoring with concurrent BU measurements at the same mines may help address the observed discrepancies and improve quantification methods to narrow the gap between BU and TD data. This paper presents a comparison of concurrent BU-TD methane emission fluxes monitored from two ventilation shafts at a southwestern Virginia, USA, longwall mine to complement the existing body of studies on satellite-based monitoring of coal mines by incorporating continuous ground-based monitoring with concurrent TD monitoring of methane emission fluxes to address the gap and provide valuable insights into temporal emission patterns. The shafts were monitored on multiple dates between 2020 and 2023. BU monitoring was performed at the exhaust fans, while TD quantifications used PRISMA hyperspectral satellite data and two different wind reanalysis datasets (i.e., GEOS-FP and ERA5). This mine also offered a unique opportunity for BU-TD comparisons before and after ventilation air methane (VAM) oxidizer operation, which was installed at one of the shafts in 2022. The results showed that TD-estimated mean fluxes were generally lower than BU data, which were attributed to quantification difficulties associated with the low albedo caused by heavy vegetation and the terrain of the area. However, despite the discrepancies in mean emissions, both the interquartile range and the data range of the distributions generally overlapped, and the estimates correctly showed the emission trends.

Virginia↗
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↗
USGS2025

A newly identified creeping strand of the Concord fault, San Francisco Bay Area

The Concord fault constitutes a major branch of the Pacific–North America transform plate boundary in Northern California, bridging the strike‐slip Bartlett Springs ‐ Green Valley Fault system to the north with the Greenville and Calaveras Faults to the south. Like many faults in the San Francisco Bay Area its long‐term slip is partially accommodated by aseismic slip (creep). Although creep has been recognized and monitored on the northern half of the fault for decades, the precise location of the southern half of the fault and its slip rate—whether accommodated seismically or aseismically—has remained enigmatic. How slip transfers between the Concord and Greenville or Calaveras faults to the south remains an outstanding question. New field observations presented here indicate that the active trace of the fault south of downtown Concord is not where previously interpreted and is indeed actively creeping. We report observations of shallow creep continuing >7 km farther south along the Concord fault than previously reported, along a fault strand not previously recognized for most of its length. This is evident as right‐laterally deflected concrete curbs and sidewalk slabs on both sides of every street that crosses the fault at a high angle in southeast Concord and northeast Walnut Creek. We document the magnitude and location of these deflections to estimate accumulated right‐lateral aseismic slip expressed in engineered structures. Offsets of these piercing lines range from 8 to 18 cm, over widths varying from narrow breaks along centimeter‐scale concrete joints to 10‐m‐wide zones of deflection. Significantly, this active trace is ∼400 m west of where the Quaternary active trace has previously been inferred, placing it within—rather than bounding—the built area of suburban Concord. Slip along the fault has already caused infrastructure damage. These results revise our understanding of the southern Concord fault and help constrain its seismic potential.

California↗
USGS2025

Combining multisite tsunami and deformation modeling to constrain slip distributions for the 1700 C.E. Cascadia earthquake

A major earthquake ruptured the Cascadia subduction zone (CSZ) on 26 January 1700. Key paleoseismic evidence associated with this event include tsunami deposits, stratigraphic evidence of coastal coseismic subsidence, written Japanese records of a tsunami unaccompanied by earthquake shaking, and margin‐wide turbidites found offshore and in lacustrine environments. Despite this wealth of independent clues, important details about this event remain unresolved. Dating uncertainties do not conclusively establish whether the proxies are from one earthquake or a sequence of them, and we have limited knowledge of the likely slip distributions of the event or events. Here, we use a catalog of 37,500 candidate synthetic ruptures between M w 7.8 and 9.2 and simulate their resulting coseismic deformation and tsunami inundation. Each model is then compared against estimated Japan tsunami arrivals, regional coastal subsidence records, and local paleotsunami deposits mapped at six different coastal marshes and one coastal lake along the CSZ. We find that seven full‐margin ruptures with a median magnitude of M w 9.1 satisfy all three constraints. We favor one M w 9.11 model that best matches all site paleoseismic observations and suggests that the Cascadia megathrust slipped up to ∼30 m and must have shallow geodetic coupling. We also find that some sequences composed of three or four ruptures can still satisfy the observations, yet no sequences of two ruptures can. Sequences are differentiated into three groups based on whether they contain a mainshock rupture located in the south (>44° N) or further north. All sequences contain unruptured portions of the megathrust and most contain mainshocks with peak slip above 40 m. The fit of the geologic evidence from sequences is poor in comparison to single‐event models. Therefore, sequences are generally less favored compared to full‐margin events.

British Columbia, California, Oregon, Washington↗
USGS2026

A ground-motion model derived using a generalized mean rupture distance for large slab interface earthquakes

Source–station distance is a central input to ground‐motion models (GMMs) for predicting seismic shaking. GMM development uses distance metrics including the Joyner–Boore distance, which is the shortest distance from an observation point to the surface projection of the earthquake rupture, and R rup the shortest distance to the rupture in three dimensions. Thompson and Baltay (2018) proposed the generalized mean rupture distance R p to address observed near‐fault ground‐motion saturation. R p accounts for the contribution to the shaking of all parts of the rupture and provides a simple method for incorporating spatially variable slip. They used R p to develop a GMM for shallow crustal earthquakes, assuming uniform slip. Here, we investigate the improvement offered by an R p ‐based GMM for large subduction interface earthquakes by recalibrating the path term for a published GMM ( Parker et al. , 2022 ) using R p derived from distributed slip models (DSMs). Inspection of within‐event and total residuals indicates that the recalibrated model fits the data at least as well as an alternative in which the path term was recalibrated using R rup and the same dataset. Incorporating slip information in its entirety or trimming the DSM geometry is preferable to assuming uniform moment release on a prescribed model fault that extends beyond the actual ruptured area. R p tuned to minimize model uncertainty is closer to the maximum distance to rupture for peak ground velocity (PGV) than acceleration, possibly reflecting the contribution to high‐frequency shaking of slip on local asperities. However, when the moment is concentrated far from stations, PGV is fit adequately with R p closer to R rup ⁠ . Our results suggest incorporating slip‐derived moment release through R p could improve GMMs for slab interface events, especially if the implementation of R p ‐based models is refined using a larger dataset of earthquakes with greater geographic diversity to account for regional ground‐motion variations.

Seismological Research Letters↗
USGS1989

An investigation of shallow ground-water quality near East Fork Poplar Creek, Oak Ridge, Tennessee

Alluvial soils and fill materials in and near the floodplain of East Fork Poplar Creek in Oak Ridge, Tennessee, are contaminated with various trace metals (primarily mercury), organic compounds, and radionuclides that were lost to the stream as a result of past operations at the Y-12 Plant, a nuclear-processing facility located within the U.S. Department of Energy’s Oak Ridge Reservation. Observation wells were installed in the shallow (above bedrock) aquifer in and near the floodplain of East Fork Poplar Creek and water-quality samples were collected to determine if contaminants found in the soils and fill are also present in the shallow ground water. The shallow aquifer in the East Fork Poplar Creek flood plain consists primarily of alluvial silt and clay with lesser amounts of sand and gravel. Thickness of the shallow aquifer ranges from essentially zero to as much as 20 feet. A silty-clay glei horizon is present between the base of the alluvium and the top of bedrock at most flood-plain locations and, where present, likely impedes downward ground-water movement. Water in the shallow aquifer near East Fork Poplar Creek occurs under water-table conditions. Recharge to the shallow aquifer is principally from precipitation, and discharge is through springs and seeps to East Fork Poplar Creek and its tributaries. During spring, summer; and fall, evapotranspiration also accounts for the removal of water in storage in the shallow aquifer. Water levels in the shallow aquifer fluctuate seasonally in response to variations in recharge and evapotranspiration. Generally, the depth to water in the observation wells ranged from about 1 to 4 feet below land surface in late winter, and from about 2 to 7 feet below land surface in late fall. During extremely dry periods, the water table recedes below the top of bedrock in some flood-plain areas, possibly causing East Fork Poplar Creek to lose water to the shallow aquifer along some reaches. Contaminants found in water samples collected from several of the observation wells in concentrations (total and (or) total-recoverable) which equaled or exceeded drinking-water standards established by the U.S. Environmental Protection Agency are antimony, chromium, lead, mercury, selenium, total phenols, and strontium-90. Total and dissolved uranium concentrations exceeded 1.0 microgram per liter in samples from nearly 70 percent of the wells in the East Fork Poplar Creek flood plain. Organic compounds that were identified in low concentrations in samples from a few wells in the flood plain are: Arochlor 1260, benzo(a)anthracene, benzo(b)fluoranthene, chrysene, 3,3-dichlorobenzidine, di-n-butylphthalate, N-nitrosodiphenylamine, and pyrene. Water from one well in the East Fork Poplar Creek flood plain at a contaminated fill site contained 37and 8 micrograms per liter of trichloroethene and trans-1,2-dichloroethene, respectively. Comparison of the results of total and (or) total-recoverable trace-metal determinations with those from dissolved determinations demonstrates that elevated concentrations of these substances in water collected from several of the wells in the East Fork Poplar Creek flood plain resulted from sorption of trace metals (and possibly organic compounds and radionuclides) by fine sediment suspended in the samples. The occurrence of contaminated sediment in these samples is suspected to be the result of borehole contamination during well installation.

Tennessee↗
USGS2025

Characterizing Meteor Crater impact melts through geochemistry and textural analysis

The U.S. Geological Survey Astrogeology Science Center houses the Meteor Crater sample collection, an assemblage of over 2,500 meters of cuttings from 161 drill holes into Meteor Crater’s rim, flanks, and ejecta blanket. We have utilized this unique collection to study the composition and spatial distribution of impact-generated materials from within the ejecta blanket. Meteor Crater has historically been known to have generated only a relatively small amount of impact melt compared to other terrestrial craters of similar size. A detailed compositional and textural dataset of impact-derived melts from this impact can therefore be a useful asset in improving our understanding of crater formation, and in particular impact melt formation. We have characterized 42 impact-melt particles from Meteor Crater using a scanning electron microscope and an electron microprobe for textural and compositional analysis. We analyzed samples from six drill holes in the ejecta blanket, situated to the northwest, southeast, south, and southwest of the crater (ejecta northeast of the crater is devoid of impact melts). Impact melts were collected from drill cuttings at various depths within the ejecta blanket, ranging from a few centimeters below the surface down to ~6.5 meters. Backscattered electron (BSE) images were acquired for each analyzed impact-melt particle. To characterize the various textures and phases present in each impact melt, we also took many detailed BSE images. Our geochemical analyses include full spectral profiles using energy dispersive X-ray spectrometry and well-calibrated wavelength dispersive spectrometry for a number of phases, including minerals (olivine, pyroxene, and so on), pristine glass, and metallic inclusions. The full dataset is available in ScienceBase as a data release (Gullikson and others, 2024), accessible at https://doi.org/10.5066/P9OGAJ8P . Our goal for this Open-File Report is to provide a summary of this immense dataset, details on data collection, descriptions of the different phases observed within impact-melt particles (both geochemically and texturally), and observable trends.

Arizona↗
USGS2025

3D Converted wave reverse time migration imaging

We describe a newly developed method for recovering high-resolution images of seismic discontinuities, such as subducting slabs, in 3D. Our method makes use of converted S → P or P → S waves observed by dense arrays of seismometers to infer the locations and relative strengths of seismic discontinuities at depth in a target region. Observed direct and converted waves are backpropagated to their times of origin. The time-reversed wavefield is then separated into its constituent P and S components via the Helmholtz decomposition, and those separated wavefields are used to compute imaging functions that characterize the locations and relative strengths of seismic discontinuities. Imaging functions may be designed to use either S → P or P → S waves, so that users can target those arrivals expected to be most dominant in a given dataset. We have previously demonstrated the efficacy of our method in two dimensions, and we now present a 3D implementation of our technique which addresses the significant computational challenges posed by the size of volumetric wavefield data in three dimensions. Through a series of synthetic examples, we demonstrate that our method is capable of recovering the fine scale structure of a subducting slab given realistic station coverage and earthquake sources. We investigate optimal seismic station geometries for our technique and explore image interpretability in regions with poor data coverage. We find that linear station geometries yield more optimal, interpretable imaging functions than collections of small arrays can. We also show that our method can successfully recover bothS → P or P → S images when realistic shear earthquake sources are used, and we explore the additional computational challenges presented by the high frequency content of S waves. Our results demonstrate the potential for our technique to recover high-resolution information about subducting slabs in real-world regions, given that relatively sparse seismic arrays with only approximately 100 stations are capable of recovering interpretable imaging functions from just a few realistic earthquake sources for multiple discontinuities at significant depth in an area of approximately 400~sq~km.

Seismica↗
USGS2024

Evidence for low effective stress within the crust of the subducted Gorda plate from the 2022 December Mw 6.4 Ferndale earthquake sequence

Stress levels on and adjacent to megathrust faults at seismogenic depths remain a key but difficult to constrain parameter for assessing seismic hazard in subduction zones. Although strong ground motions have been observed to be generated from distinct, high-stress regions on the downdip end of the megathrust rupture areas in many great earthquakes, we lack direct constraints on the stress level in the lower seismogenic portion of the Cascadia megathrust. On 2022 December 20, a M w 6.4 strike-slip earthquake occurred near Ferndale, California in southern Cascadia and likely ruptured the Gorda slab crust in the lower seismogenic portion, providing an opportunity to assess the stress level in this region. Here, we relocate the Ferndale mainshock and the first two weeks of aftershocks using a high-resolution 3-D velocity model and estimate rupture dimensions, directivity, and stress drop for several M w 4-5 aftershocks and recent earthquakes. The aftershocks define a strike-slip fault in the slab crust striking ENE, consistent with the mainshock focal mechanism. The orientation of this fault is about 45 ° off the ideally oriented fault plane given the stress state in the slab. The aftershock zone is extensive and broad in the forward direction of the mainshock rupture but still constrained within the volume of high Vp/Vs within the slab crust. Our stress drop estimates are generally lower for M w 4-5 earthquakes located within the slab crust compared to those a few km deeper in the slab mantle. Combined, our results support a relatively low effective stress level in the vicinity of the megathrust in the lower portion of the seismogenic zone in southern Cascadia, likely due to elevated fluid pressures. Consequently, the ground motion in the onshore region above this low-stress seismogenic portion in southern Cascadia may not be as intense as that observed during great earthquakes in other subduction zones.

California↗
USGS2026

Arsenic and isotope concentrations in the lower Platte River valley of eastern Nebraska, early 1970s to 2023

The City of Lincoln, Nebraska, has been monitoring concentrations of arsenic in their source water and evaluating their options for treatment and removal since at least 2002. In 2022, the City of Lincoln, Nebr., with funding assistance from the Nebraska Water Sustainability Fund, began cooperating with the U.S. Geological Survey to examine arsenic concentrations in surface water and groundwater in the lower Platte River valley and the area around City of Lincoln Water System (LWS) well field. Arsenic data collected from the Platte River since 1974 were examined using the “weighted regression on time, discharge, and season” model, which compared the streamflow (also referred to as “discharge”), time of year, and season to estimate concentrations of arsenic. Annual mean arsenic concentrations modeled for more than 49 years at the Platte River at Louisville, Nebr., U.S. Geological Survey streamgage (station 06805500), indicated a significant increasing trend. Arsenic concentrations in the Platte River were seasonal, with the highest concentrations being observed during mid- to late summer. When seasonal patterns and streamflow were combined with arsenic concentrations in the Platte River during low streamflow conditions, groundwater contributions, which can have higher arsenic concentrations, make up a larger portion of the streamflow. Arsenic samples were collected from upstream rivers in 2022 and 2023 and were paired to analyze the arsenic contributions at the U.S. Geological Survey streamgage on the Platte River near Ashland, Nebr. (station 06801000), near the City of Lincoln well field. The arsenic concentrations from the streamgage on the Platte River near Ashland, Nebr., location, were higher than the U.S. Geological Survey streamgage on the Elkhorn River at Waterloo, Nebr. (station 06800500), and significantly lower than at the U.S. Geological Survey streamgage on the Platte River near Leshara, Nebr.(station 06796500), indicating that the Platte River usually contributes a higher concentration of arsenic than does the Elkhorn River as they join near Ashland, Nebr. During 1991–2023, six groundwater monitoring wells were analyzed to identify trends in arsenic concentrations. Two of the six wells had a positive trend during the 33-year period. One monitoring well did not reveal a long-term trend during this period but showed a trend during 2019–23, correlating to a period when the island in the middle of the Platte River was connected to the east bank of the river when manganese reducing conditions were present and groundwater levels were declining in the well. Across all wells the oxidation and reduction (redox) condition during the time of sampling was assessed. Mixed anoxic and (or) oxic redox condition was the most common redox process and the highest sampled arsenic concentrations in monitoring wells were observed in anoxic conditions driven by manganese reduction. Groundwater arsenic concentrations had seasonal variation around the City of Lincoln well field, with higher arsenic concentrations tending to be further south in comparison to samples collected further north. Isotope samples were collected and analyzed in surface water and groundwater around the LWS well field. The samples indicate that the proportion of surface water present in the LWS production wells can be higher in the spring and lower in the summer. With higher arsenic concentrations observed in the stream water during the summer period, the LWS source water can be affected by these elevated arsenic concentrations even though the proportion of surface water is lower.

Nebraska↗
USGS2024

Limited evidence of late Quaternary tectonic surface deformation in the eastern Tennessee seismic zone, USA

The ~300-km-long eastern Tennessee seismic zone (ETSZ), USA, is the second-most seismically active region east of the Rocky Mountains. Seismicity generally occurs below the Paleozoic fold-and-thrust belt within the Mesoproterozoic basement, at depths of 5–26 km, and earthquake magnitudes during the instrumental record have been moment magnitude ( M w )≤4.8. Evidence of surface deformation may not exist or be difficult to detect because of the vegetated and soil-mantled landscape, landslides, locally steep topography, anthropogenic landscape modification, or long, irregular recurrence intervals between surface-rupturing earthquakes. Despite the deep seismicity, analog models indicate that accumulation of strike-slip or oblique-slip displacement at depth could be expected to propagate upward through the Paleozoic section, producing a detectable surficial signal of distributed faulting. To identify potential surface deformation, we interrogated the landscape at different spatial scales. We evaluated morphotectonic and channel metrics, such as channel sinuosity and catchment-scale hypsometry. Additionally, we mapped possible fault-related topographic features on 1-m lidar. Finally, we integrated our observations with available bedrock and Quaternary surficial mapping and subsurface geophysical data. At a regional scale, most morphotectonic and channel metrics have a strong lithologic control. Within smaller regions of similar lithology, we observe changes in landscape metrics like channel sinuosity and catchment-scale hypsometry that spatially correlate with new lineaments identified in this study and previously mapped east–west Cenozoic faults. These faults have apparent left-lateral offsets, are optimally oriented to slip in the current stress field, and match kinematics from recent focal mechanisms, but do not clearly preserve evidence of late Pleistocene or Holocene tectonic surface deformation. Most newly mapped lineaments might be explained by either tectonic or non-tectonic origins, such as fluvial or karst processes. We also re-evaluated a previously described paleoseismic site and interpret that the exposure does not record evidence of late Pleistocene faulting but instead is explained by fluvial stratigraphy.

Tennessee↗
USGS2026

Lunar analog study using portable gamma-ray neutron detector: Radiochemical mapping of silicic-to-basaltic volcanic terrains in the San Francisco Volcanic Field, Arizona

Studying lunar silicic volcanism provides key insights into the Moon’s crustal evolution, magmatic processes, and volcanic history. The Lunar-VISE (Lunar Vulkan Imaging and Spectroscopy Explorer) mission will investigate the Gruithuisen domes, a unique lunar region hypothesized to have formed through silicic volcanism. Using instruments on a Firefly Aerospace lander and a Honeybee Robotics rover, Lunar-VISE will analyze mineralogy, geochemistry, and surface properties to determine the origin and evolution of the domes, with a gamma-ray and neutron spectrometer (LV-GRNS) among its payload instruments. In preparation for this mission, we conducted preliminary fieldwork using a handheld gamma-ray neutron detector with NaI(Tl) and Cs LiYCl :Ce (CLYC) scintillators. This study focuses on various rhyolitic and basaltic volcanic centers in the San Francisco Volcanic Field (SFVF) near Flagstaff, Arizona—specifically Sugarloaf Peak, Bonito Lava Flow, and Robinson Mountain, a region containing several well-characterized lunar analog sites. The SFVF was selected for this study due to its broad compositional diversity spanning basaltic to rhyolitic compositions, its well-preserved volcanic morphologies, and its extensive use in previous NASA field campaigns and astronaut training exercises, making it an ideal terrestrial laboratory for testing planetary exploration techniques. We measured natural radioactivity, specifically from potassium (K), thorium (Th), and uranium (U) and other elements in their decay chains, which serve as diagnostic tracers of magmatic differentiation and crustal evolution processes, to assess geochemical variability across compositionally diverse terrains. The detector’s sensitivity was assessed across varying concentration levels. Regionally averaged concentrations of radioisotopes were determined by gamma-ray spectroscopy at selected sites. Our measurements reveal significant compositional variations between sites, with Sugarloaf Peak (rhyolitic, silicic dome) exhibiting the highest average radioisotope concentrations (K: 3.29 ± 0.24 wt%, U: 14.35 ± 1.81 ppm, Th: 27.14 ± 2.43 ppm), while Bonito Lava Flow (K: 1.13 ± 0.08 wt%, U: 3.86 ± 0.52 ppm, Th: 7.52 ± 1.02 ppm), and Robinson Mountain (K: 1.46 ± 0.16 wt%, U: 5.59 ± 1.06 ppm, Th: 11.75 ± 1.98 ppm) show lower concentrations aligned with basaltic compositions. Gamma-ray fluxes were elevated by a factor of approximately three to four at Sugarloaf Peak relative to nearby basaltic terrains, consistent with expected geochemical differentiation patterns. Furthermore, at meter scales, proximity to geological features significantly affects measurements. Th concentrations adjacent to a cliff face at Sugarloaf’s base were 28% higher than values measured 3-5 meters away from the same feature, demonstrating localized compositional heterogeneity. Analysis of station-by-station measurements reveals that K, U, and Th concentrations show an elevation trend at Robinson Mountain and generally higher concentrations near Sugarloaf Peak’s summit, suggesting progressive magmatic differentiation and/or the presence of more evolved lithologies at higher elevations. By mapping these elements on Earth using GRNS, we aim to optimize measurement techniques for Lunar-VISE and similar rover-borne missions by demonstrating the applicability of portable GRNS instruments for planetary surface exploration in an analog environment.

Arizona↗
USGS2024

Parameter ESTimation with the Gauss–Levenberg–Marquardt algorithm: An intuitive guide

In this paper, we review the derivation of the Gauss–Levenberg–Marquardt (GLM) algorithm and its extension to ensemble parameter estimation. We explore the use of graphical methods to provide insights into how the algorithm works in practice and discuss the implications of both algorithm tuning parameters and objective function construction in performance. Some insights include understanding the control of both parameter trajectory and step size for GLM as a function of tuning parameters. Furthermore, for the iterative Ensemble Smoother (iES), we discuss the importance of noise on observations and show how iES can cope with non-unique outcomes based on objective function construction. These insights are valuable for modelers using PEST, PEST++, or similar parameter estimation tools.

Groundwater↗
USGS2026

Petrography and mineralogy of selected pre-Middle Jurassic basement rocks beneath the Atlantic and Gulf Coastal Plains in Florida

Florida is covered by flat-lying sedimentary strata of the Atlantic and Gulf Coastal Plains. These strata have accumulated since Middle Jurassic time. The pre-Middle Jurassic, or basement, rocks that underlie the Coastal Plain in Florida are known only from drill cores and cuttings recovered from a relatively small number of boreholes. This data report presents petrographic observations and the results of X-ray diffraction analyses of basement rocks from 18 boreholes across Florida in support of the identification, discrimination, and correlation of units for a subcrop geologic map of pre-Middle Jurassic rocks composing Florida’s sub-Coastal Plain geology.

Florida↗