Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “PLoS ONE”

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,186 records · Page 8Linked to original sources
USGS2024

Onset of aftershocks: Constraints on the Rate-and-State model

Aftershock rates typically decay with time t after the mainshock according to the Omori–Utsu law, R (t)=K(c+t) −p ⁠ , with parameters K , c , and p . The rate‐and‐state (RS) model, which is currently the most popular physics‐based seismicity model, also predicts an Omori–Utsu decay with p = 1 and a c ‐value that depends on the size of the coseismic stress change. Because the mainshock‐induced stresses strongly vary in space, the c ‐value should vary accordingly. Short‐time aftershock incompleteness (STAI) in earthquake catalogs has prevented a detailed test of this prediction so far, but the newly developed a ‐positive method for reconstructing the true earthquake rate now allows its testing. Using previously published slip models, we calculate the coseismic stress changes for the six largest mainshocks in Southern California in recent decades and estimate the maximum shear as a scalar proxy of the coseismic stress tensor. Aftershock rates reconstructed for events in different stress ranges show that the rates follow a power law with p = 1 independent of stress with no clear sign of a c ‐value. The onset of the power‐law decay is abrupt and more delayed in areas with smaller stress changes. The observations do not necessarily contradict the RS model, as STAI limits the resolution for early aftershocks, and the RS model can reproduce the observations for specific Aσ values. However, the observations lead to strong constraints, namely Aσ <10 kPa and a power‐law decay of the background rate with distance to the fault, with exponent 2.7.

Seismological Research Letters↗
USGS2025

System characterization report on the Environmental Mapping and Analysis Program (EnMAP)

This report addresses system characterization of the Environmental Mapping and Analysis Program hyperspectral sensor by the DLR (German Aerospace Center, ground segment project management), GFZ (Deutsches Geoforschungszentrum, science lead) and is part of a series of system characterization reports produced and delivered by the U.S. Geological Survey Earth Resources Observation and Science Cal/Val Center of Excellence. These reports present and detail the methodology and procedures for characterization; present technical and operational information about the EnMAP hyperspectral sensor; and provide a summary of test measurements, data retention practices, data analysis results, and conclusions. The Earth Resources Observation and Science Cal/Val Center of Excellence system characterization team completed data analyses to characterize the geometric (interior and exterior), and radiometric performances of the EnMAP hyperspectral sensor. Results of these analyses indicate that the Environmental Mapping and Analysis Program has a band-to-band geometric performance in the range of −0.135 to 0.15 pixel, geometric performance relative to the Operational Land Imager in the range of −27.716 meters (−0.92 pixel) to 32.892 meters (1.09 pixels) offset in comparison to Landsat 8 Operational Land Imager, offset of a radiometric comparison in the range of −0.012 to 0.020, slope of a radiometric comparison in the range of 0.947 to 1.031.

Open-File Report↗
USGS2026

Software to support remote sensing of river discharge based on critical flow theory

Water resource management requires accurate observations of streamflow but standard field methods for measuring river discharge ( Q ) are costly and can be hazardous for equipment and personnel. Remote sensing has become a viable alternative, but many image-based techniques require field data for calibration and depth and velocity can seldom be mapped with a single sensor. A new approach based on critical flow theory, in contrast, allows both of these attributes to be inferred from readily available image data. This technique only pertains to sites with standing waves, called undular hydraulic jumps (UHJs), but a recent investigation demonstrated its potential to provide accurate discharge estimates. This paper introduces software designed to facilitate Inferring Q from UHJs Identified in River Images (InQUIRI). The package includes modules for retrieving data from image servers, making the measurements of wavelength and width required to calculate discharge, inferring a representative wavelength from a profile digitized along a wave train, combining multiple estimates to obtain an ensemble median discharge, and assessing accuracy via comparison to gage records from the U.S. Geological Survey. By making these steps easier to implement, InQUIRI enables users to apply the workflow to a variety of UHJ-containing images. Accumulating more case studies, some successful and others less so, would help constrain the range of applicability of the critical flow approach and foster development of refined guidelines for selecting and measuring waves. The software described herein could play an important role in promoting informed use of this new technique for non-contact streamflow measurement.

Arizona, Colorado, New Mexico, Utah↗
USGS2026

Site response models based on geometric parameters for southern California sedimentary basins

Site response in sedimentary basins is influenced by complex three-dimensional (3D) features, including trapping of seismic waves, focusing of seismic energy and basin resonance. Current ground motion models (GMMs) incorporate basin effects using one-dimensional parameters like V S30 and shear wave velocity isosurface depths, which are limited in capturing lateral and 3D effects. To address these limitations, we develop seismic site response models based on novel parameters that represent multi-dimensional properties of the Los Angeles Basin (LAB) geometry and shear wave velocity. We define a basin shape for the LAB using depth to subsurface geologic interfaces associated with the oldest sedimentary deposits (depth to a particular shear wave velocity horizon, i.e., 1.5 km/s - z 1.5 ) and the depth to the crystalline basement ( z cb ) which are determined using geologic cross sections and community seismic velocity model profiles. We explore a suite of geometric descriptors computed for the LAB and southern California, from which three parameters with the greatest predictive potential are selected and evaluated using empirical ground motion residual analyses in combination with the Boore et al. GMM. The results demonstrate that the zonal heterogeneity index ( ), standard deviation of the absolute difference between z 1.5 and z cb ( ) and standard deviation of z cb ( ) each provide a reduction in site-to-site variability ( ϕ S2S ) of empirical GMMs. The reduction in ϕ S2S is period-dependent, with average decreases of 3%, 26% and 6% for , , and , respectively. Although these reductions are modest from an engineering application perspective, they are statistically significant, underscoring the inherent difficulty in fully characterising complex basin effects. Collectively, these findings indicate that the inclusion of basin-specific geometric parameters yields measurable, albeit incremental, improvements in site response prediction and establishes a framework for the progressive refinement of seismic hazard characterisation within sedimentary basins.

California↗
USGS2026

Earthquake rupture arrest from depth‐dependent frictional stability on the Pütürge segment of the East Anatolian Fault Zone, Türkiye (Turkey)

Determining why earthquake ruptures stop is a central challenge in earthquake science and seismic hazard assessment. The Pütürge segment of the East Anatolian Fault Zone, Türkiye, exhibits shallow creep (<∼3 km depth) yet hosts M > 6.5 earthquake ruptures at greater depth. Here, we evaluate whether variations in frictional stability along this segment aided arrest of the 2020 M 6.7 Elaziğ and 2023 M 7.8 Pazarcık earthquake ruptures. Analysis of Sentinel‐1 Synthetic Aperture Radar imagery indicates the 2023 M 7.8 rupture propagated laterally into a metamorphic massif within the Pütürge segment, where slip rapidly decayed below detection limits. Creepmeters along this segment recorded no significant surface afterslip (<3 mm) in the following year. To investigate this fault‐slip behavior, we conducted triaxial friction experiments on Pütürge fault gouge sampled from an outcrop exposure. The gouge, composed primarily of muscovite, quartz, and calcite, is velocity strengthening at conditions approximating 0–2.5 km depth and velocity weakening at 4–5 km depth. This transition to velocity‐weakening friction is associated with enhanced comminution and shear localization observed microstructurally. Our results suggest that depth‐dependent frictional stability of the Pütürge fault segment facilitates rupture nucleation and propagation at depth while maintaining rupture‐arresting behavior near Earth’s surface.

Pütürge segment of the East Anatolian Fault Zone↗
USGS2025

Manipulation of the Symbiodiniaceae microbiome confers multigenerational impacts on symbioses and reproductive ecology of its Exaiptasia diaphana host

Symbiodiniaceae-associated microbiota strongly influence cnidarian symbioses. We systematically reduced the bacterial and fungal communities associated with Symbiodiniaceae to study potential effects on the cnidarian holobiont Exaiptasia diaphana (Aiptasia). Clonal anemones were inoculated with xenic Breviolum minutum (SSB01) and microbiome-manipulated cultures after antibacterial or antifungal treatment. The asexual reproduction of pedal laceration allowed for three generations of clonal aposymbiotic Aiptasia to be utilised in this study, from the initial adult generation (G0), to the first (G1), and second (G2) generation. We inoculated small and large G1 Aiptasia with SSB01 algae and monitored onset of symbiosis, rate of algal proliferation, and holobiont characteristics. Sequencing the 16S and 18S rRNA gene regions identified significant differences in the bacterial and fungal communities of the G0 and G1 generations, alongside differences between the size classes of small and large G1 anemones. The microbiome of larger G1 individuals was distinct to the smaller G1 anemones, suggesting a microbiome maturation process. Control Breviolum minutum cultures exhibited a significantly greater proliferation rate in large G1 anemones when compared to antibacterial or antifungal treated cultures, whereas the opposite trend was documented in the small G1 anemones. Although no differences were observed between algal photochemical parameters, or the growth and polyp activity of G1 juveniles, we observed a significant influence in the production of G2 clones between treatments. Overall, we provide strong ecological implications of manipulating Symbiodiniaceae microbiome, not for the algae themselves, but for the maturation of the host Aiptasia, as well as for the cnidarian holobiont over multiple generations.

The ISME Journal↗
USGS1976

An environmental assessment of impacts of coal development on the water resources of the Yampa River basin, Colorado and Wyoming; Phase I work plan

Coal resources of the western United States are being developed at ever-increasing rates, causing concerns of the effects of mining and associated activities on the environment. The Yampa River basin in northwestern Colorado and south-central Wyoming is undergoing economic development of its coal, oil and gas, and uranium resources. The Yampa River basin assessment is a 21-year program of the U.S. Geological Survey. It is designed primarily to assess the availability and quality of the basin's water resources. The basin assessment also will evaluate potential environmental and selected socioeconomic impacts of energy-resource development plans proposed by mining and power companies. This report serves as a project work plan for the basin assessment's first-phase work activities. Economic development in the Yampa River basin will result in higher levels of residuals (noneconomic byproducts) being discharged into the environment. Residuals-management encompasses the technological and economic interrelationships that should be considered in order to comply with Federal and State environmental-control regulations regarding air quality, water quality, and disposal of solid wastes. The basin assessment, the first phase that is described by this report, is designed to provide Federal, State, and local decision-makers with the basic environmental information needed to formulate and to evaluate policies for the development of energy resources. The techniques applied and conclusions reached in the Yampa River basin assessment should aid similar studies of other energy-rich basins in the western United States.

Colorado, Wyoming↗
USGS2024

System characterization report on the Earth Surface Mineral Dust Source Investigation (EMIT) sensor

Executive Summary This report addresses system characterization of the Earth Surface Mineral Dust Source Investigation (EMIT) sensor, an imaging spectrometer developed by the National Aeronautics and Space Administration. This report is part of a series of system characterization reports produced and delivered by the U.S. Geological Survey Earth Resources Observation and Science Cal/Val Center of Excellence. These reports present and detail the methodology and procedures for characterization; present technical and operational information about the specific sensing system being evaluated; and provide a summary of test measurements, data retention practices, data analysis results, and conclusions. The Earth Resources Observation and Science Cal/Val Center of Excellence system characterization team completed data analyses to characterize the geometric (interior and exterior) and radiometric performances. Results of these analyses indicate that the EMIT sensor has a band-to-band geometric performance in the range of −0.355 to 0.210 pixel with a few exceptions of shortwave infrared channels. Geometric offset relative to the Landsat 8 Operational Land Imager ranged from −15.966 meters (−0.266 pixel) to 43.844 meters (0.731 pixel). Offset of a radiometric comparison ranged from −0.016 to 0.025, and slope of a radiometric comparison ranged from 0.837 to 0.985. EMIT agreed with Radiometric Calibration Network measurements within 5 percent across most of the spectral channels.

Open-File Report↗
USGS2026

Effects of a low-lipid diet on the gut microbiome and head kidney transcriptome of juvenile Chinook Salmon

Objective Pacific salmon Oncorhynchus spp. reared in production hatcheries are typically fed high-lipid, energy-dense diets to achieve large size and high body condition prior to release. In contrast, juveniles in natural environments tend to consume low-lipid, high-protein diets, and fish reared for research or conservation purposes are sometimes fed diets that are formulated to mimic natural diets and promote wild-like phenotypes. Understanding how these alternative diets affect fish health beyond growth and body condition could ultimately contribute to improving hatchery fish fitness. Methods In this work, we evaluated changes in the fecal microbiome and gene expression of juvenile Chinook Salmon O. tshawytscha on a standard high-lipid hatchery diet versus a low-lipid diet formulated to mimic the nutrition profile of natural-origin fish. To evaluate the time scale at which diet alters the fecal microbiome, we collected longitudinal samples over a 12-week period and switched the diets of a subset of fish twice during the experiment. We used 16S ribosomal RNA gene amplicon sequencing to characterize fecal microbiome differences between fish on the two diets as well as hatchery-reared fish at a production hatchery, hatchery fish that had been captured after release into a stream, and natural-origin, stream-reared fish of similar ages. Additionally, we conducted RNA sequencing on head kidney samples from laboratory-reared fish to evaluate changes in gene expression in this important immune organ. Results We found that the low-lipid diet and the hatchery diet resulted in microbiomes that differed from the microbiome of natural-origin fish and from each other and that diet-driven changes to the microbiome could occur in under 14 d. The low-lipid diet did not result in a microbiome that resembled the microbiome of naturally produced fish. Instead, the low-lipid diet resulted in a microbiome community that was distinct from those of fish reared on the hatchery diet and fish sampled from the wild. The RNA sequencing results indicated differential enrichment of pathways related to immunity, metabolism, and hormone synthesis between fish that were fed the two experimental diets. Conclusions The results suggest that additional environmental factors influence the microbiome more strongly than diet formulation or that the low-lipid diet has a smaller effect on the microbiome than a natural, ­invertebrate-based diet. Given that the gut microbiome and systemic immune function contribute significantly to disease resistance, our findings highlight the importance of understanding how diets fed to fish in captivity may affect fish health beyond growth and body condition metrics.

Journal of Aquatic Animal Health↗
USGS2026

Groundwater dependency and hydroclimatic influences on riparian and upland vegetation productivity, Upper San Pedro, Arizona, United States

In arid and semi-arid regions, groundwater sustains vegetation through subsurface water access, yet the responses of groundwater-dependent ecosystems (GDEs) to changing hydroclimate and groundwater availability are relatively understudied. This study investigates seasonal and spatial patterns in vegetation greenness using Landsat Enhanced Vegetation Index (EVI) values across riparian and upland zones in the semi-arid Upper San Pedro (USP) watershed, southern Arizona, which experiences a bimodal precipitation regime. We paired 25 years (2000–2024) of EVI and depth to groundwater (DTG) data from 89 wells and climate metrics (precipitation and vapour pressure deficit) to quantify the sensitivity of vegetation to subsurface moisture as well as atmospheric moisture supply and demand. Vegetation at wells near the USP riparian area showed strong associations between EVI and DTG anomalies during the monsoon season, indicating sustained groundwater use even during this wet period when summer precipitation is abundant. In contrast, upland vegetation that lacked access to groundwater showed minimal sensitivity in EVI to DTG and was generally less responsive to vapour pressure deficit. Interestingly, the riparian GDEs were not decoupled from precipitation and climate variability. These results underscore the importance of groundwater for maintaining riparian productivity and highlight the utility of remote sensing in identifying vegetation-climate-groundwater linkages across heterogeneous dryland landscapes.

Arizona↗
USGS2025

Effect of land cover type on 3D deformation recovery from synthetically deformed high-resolution satellite optical imagery

The limits of detection for earthquake surface deformation in the spatial domain have improved with advances in remote sensing imagery data availability, resolution, and analysis. Sub-pixel correlation and digital elevation model (DEM) differencing from sub-meter, earthquake-spanning satellite optical imagery has enhanced surface rupture mapping and deformation measurements. However, knowledge of measurement accuracy and uncertainty is limited. To address this, we construct orthophotos and digital elevation models (DEMs) from repeat high resolution (∼0.5 m) satellite optical imagery along two sections of the Garlock fault, California with clear fault geomorphology and differing land cover. We deform later sets of DEMs and images with synthetic earthquakes containing both diffuse and discrete horizontal and vertical displacements. Sub-pixel image correlation and DEM differencing demonstrate how vegetation degrades recovered displacement accuracy. In barren land cover, horizontal displacements are detectable to an expected ∼1/10th-pixel size. With shrubs, trees, and grass, detectable displacements increase to >1/2-pixel size, and filtering results by correlation score and using elevation values as input rather than image values improves accuracy. Vertical displacement detection thresholds remain lower in vegetation, at >1-pixel size. Higher slope angles degrade displacement recovery, worsened by vegetation. Diminishing seasonal separation improves accuracy over vegetated regions, though not to the level achieved in barren environments. These results will inform research and operational efforts on the utility of high resolution satellite optical imagery for detecting deformation in varied land cover. Furthermore, they reveal where alternative measurements, such as from LiDAR or radar interferometry, are required to mitigate the effects of vegetation and capture fine-scale crustal deformation.

California↗
USGS2026

Making many out of one: Synthetic geologic deformation model distributions for use in USGS NSHM25‐PRVI Puerto Rico-U.S. Virgin Island update

A key use‐case of geologic slip rates is within deformation models used in probabilistic seismic hazard analyses. Field‐derived geologic slip rates have formed the cornerstone of deformation models in such applications for decades. Recent advancements in seismic hazard analyses have expanded the use of faults for which geologic slip rates are not well constrained using categorical slip rate estimates. Because of these advancements, application of a geologic deformation model for use in 2025 U.S. Geological Survey National Seismic Hazard Model Puerto Rico‐U.S. Virgin Islands (NSHM25‐PRVI) proved challenging due to: (1) a lack of field‐based geologic slip rates, and (2) a lack of epistemic uncertainty distributions within a broad range of estimated slip rates. Preliminary versions of the NSHM25‐PRVI model sampled these slip rate bins in a coincident manner along preferred and extreme value branches, which yielded untenable correlations in mean hazard results. To minimize the influence of correlated uncertainties amid these challenges, we develop a synthetic epistemic uncertainty distribution for deformation rate on each crustal fault. Each fault has a weighting schema across four possible distribution shapes: uniform, normal, triangular favoring local minima, and triangular favoring local maxima. The synthetic distributions are then sampled several times for each logic tree branch. The results provide a more realistic distribution of rates across the study region as compared with using correlated extrema sampling. This exploration of our method in a small region like PRVI can pave the way for larger‐scale, more complicated applications (e.g., western United States).

Puerto Rico, US Virgin Islands↗
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↗
USGS2026

Simulation of groundwater flow to evaluate hydrogeologic controls on a PFAS plume, Coakley Landfill Superfund Site, Rockingham County, New Hampshire

Per- and polyfluoroalkyl substances (PFAS), including perfluorooctanoic acid (PFOA) and perfluorooctanesulfonic acid (PFOS), have been detected at combined concentrations above 2,000 nanograms per liter (ng/L) at groundwater seep locations near the Coakley Landfill Superfund site, in North Hampton, New Hampshire. The landfill was active from 1972 to 1985. An impermeable cap was placed on the landfill in 1998. The adjacent area to the Coakley Landfill has many water supply wells, and transport of PFAS compounds to the wells is a concern. Fracture anisotropy in the underlying bedrock aquifer complicates the understanding of PFAS transport because groundwater preferentially travels along fractures that may not align with the prevailing groundwater flow direction. In 2018, the U.S Environmental Protection Agency and the U.S. Geological Survey began an investigation of the groundwater flow from the Coakley Landfill site. This report describes the modification of a numerical groundwater-flow model for the local area around the Coakley Landfill and summarizes findings of the investigation. In addition, this report includes a brief description of PFOA and PFOS occurrence, a discussion of model construction, evaluation of model performance through calibration, and discussion of simulation results for two periods (before and after capping). Limitations are also discussed. Results show that simulated groundwater flow moves from the Coakley Landfill to the west and north. Advective transport modeling using particle tracking shows that groundwater from the landfill discharges primarily to streams to the west and north, and a small amount is transported to distal wells. Dilution of contaminants through advection and dispersion likely plays a role in whether PFAS compounds from the landfill will be detected above laboratory reporting levels at distal wells.

EarthArXiv↗
USGS2026

Effect of mineral deposit data on predictions from the three-part approach to quantitative mineral resource assessment—A study of 16 previous U.S. Geological Survey assessments

The three-part approach to quantitative mineral resource assessment requires information about the properties of undiscovered mineral deposits in an assessment area. These properties are unknown, so the properties of discovered mineral deposits of the same mineral deposit type are used instead. In the three-part approach, these discovered mineral deposits come from around the world, and their properties constitute the pooled data for that mineral deposit type. Alternatively, these discovered mineral deposits could come from the assessment area, and their properties constitute the tract data for that mineral deposit type. Tract data may be more representative of the undiscovered mineral deposits in the assessment area than the pooled data. The goal of this study was to determine whether resource predictions using pooled data are equivalent to resource predictions using tract data. To this end, 16 previous U.S Geological Survey assessments were studied. For each assessment, resources were predicted for one undiscovered mineral deposit in the assessment area. One set of predictions used pooled data, and another used tract data. The two sets of predictions were compared with an equivalence test, using the six assessment statistics that are commonly reported for mineral resource assessments. Practical equivalence is the condition that two corresponding assessment statistics are within a factor of 1.5 of one another. For each of 2 assessments, all 6 assessment statistics were practically equivalent. For both assessments, the assessment statistics from the pooled data, relative to the corresponding assessment statistics from the tract data, ranged from 1.30 times smaller to 1.03 times larger. For each of 14 assessments, 1 or more of the 6 assessment statistics were not practically equivalent. The assessment statistics from the pooled data, relative to the corresponding assessment statistics from the tract data, ranged from 26.6 times smaller to 5.53 times larger. The use of pooled data has been a standard procedure in the three-part approach since at least 1986. The 16 assessments in this study are not a representative sample of those prior assessments that used pooled data. So, it is inappropriate to use the study results to infer whether pooled data affected the resource predictions for those prior assessments.

Scientific Investigations Report↗
USGS2000

Quality assurance and analysis of water levels in wells on Pahute Mesa and vicinity, Nevada Test Site, Nye County, Nevada

Periodic and continual water-level data from 1963 to 1998 were compiled and quality assured for 65 observation wells on Pahute Mesa and vicinity, Nye County, Nevada. As part of the quality assurance of all water levels, ancillary data pertinent to computing hydraulic heads in wells were compiled and analyzed. Quality-assured water levels that were not necessarily in error but which did not represent static heads in the regional aquifer system, or required some other qualification, were flagged. Water levels flagged include those recovering from recent pumping or well construction, water levels affected by nuclear tests, and measurements affected by borehole deviations. A cursory examination of about 30 wells with available water-level and down-hole temperature data indicate that water levels in most wells on Pahute Mesa would not be significantly affected by temperature if corrected to 95 degrees Fahrenheit. Wells with large corrections (greater than 10 feet) are those with long water columns (greater than 1,500 feet of water above the assumed point of inflow) in combination with mean water-column temperatures exceeding 105 degrees Fahrenheit. Water-level fluctuations in wells on Pahute Mesa are caused by several factors including infiltration of precipitation, barometric pressure, Earth tides, ground-water pumpage, and seismic events caused by tectonic activity and underground nuclear testing. No observed water-level fluctuations were attributed to a naturally occurring earthquake. The magnitude and duration of changes in water levels caused by nuclear tests are affected by the test size and the distance from a well to the test. Identifying water levels that might be affected by past nuclear tests is difficult because pre-testing water-level data are sparse. Hydrologically significant trends were found in 13 of 25 wells with multiple years of water-level record. The largest change in water levels (1,029 feet in 25 years) occurred in well U-19v PS 1D as a result of the Almendro nuclear test. Likely explanations for trends in most of the wells are either changes in precipitation patterns that affect recharge rates to the ground-water system, pumping effects from water-supply well U-20 WW, or a combination of these two factors.

Nevada↗
USGS2024

Intense alteration on early Mars revealed by high-aluminum rocks at Jezero Crater

The NASA Perseverance rover discovered light-toned float rocks scattered across the surface of Jezero crater that are particularly rich in alumina ( ~ 35 wt% Al 2 O 3 ) and depleted in other major elements (except silica). These unique float rocks have heterogeneous mineralogy ranging from kaolinite/halloysite-bearing in hydrated samples, to spinel-bearing in dehydrated samples also containing a dehydrated Al-rich phase. Here we describe SuperCam and Mastcam-Z observations of the float rocks, including the first in situ identification of kaolinite or halloysite on another planet, and dehydrated phases including spinel and apparent partially dehydroxylated kaolinite. The presence of spinel in these samples is likely detrital in origin, surviving kaolinitization, pointing to an ultramafic origin. However, the association of low hydration with increased Al 2 O 3 abundances suggests heating-induced dehydration which could have occurred during the lithification or impact excavation of these rocks. Given the orbital context of kaolinite-bearing megabreccia in the Jezero crater rim, we propose an origin for these rocks involving intense aqueous alteration of the parent material, followed by dehydration/lithification potentially through impact processes, and dispersion into Jezero crater through flood or impact-related processes.

Communications Earth & Environment↗
USGS2025

Fluid inclusion constraints on the geometry of the magmatic plumbing system beneath Mauna Loa – Part 2: Xenoliths

Mauna Loa volcano erupts crystal-poor material at its summit and more crystal-rich material on its rift zones. Some of the more olivine-rich lava flows contain xenoliths with diverse mineralogy, including cumulate harzburgites with high-Mg# orthopyroxenes and high-Fo olivines (both > 84). Previous experimental work and thermodynamic modelling has proposed that high-Mg# orthopyroxenes only crystallize from Mauna Loa melts at high pressures (> 6 kbar, > 20 km), leading to suggestions that there is a region of sub-Moho magma storage at Mauna Loa in addition to the geophysically imaged magma reservoir at 2–5 km depth below the summit. We use melt and fluid inclusion barometry combined with thermodynamic models to further investigate this suggestion. Fluid inclusion data from harzburgites and dunitic xenoliths yield storage depths remarkably similar to those found in non-xenolithic crystals from lavas and tephras, with a clear peak at ~ 2–3 km (below the summit). Depths from melt inclusions in these xenoliths overlap with fluid inclusion pressures, ruling out the possibility of fluid inclusion re-equilibration during a period of stalling in a shallower reservoir. We examine five different thermodynamic models and find that the minimum pressure of olivine-orthopyroxene co-saturation varies by ~ 4 kbar (~ 12 km). These models also fail to predict that orthopyroxene is stable in ~ 15–80% of compositionally relevant experimental charges which grew orthopyroxene. Overall, this shows that phase stability modelling is an unreliable method of determining magma storage depth at Mauna Loa. We suggest that model discrepancies reflect a lack of experimental constraints on orthopyroxene stability at > 1200 ℃ and 0.01–5 kbar. Based on the presence of large oikocrystic orthopyroxenes completely enclosing rounded olivine chadacrysts, we suggest that these harzburgitic xenoliths formed through the reaction of intruding melts with olivine mush piles within the Mauna Loa edifice at ~ 3 km depth below the summit, with no need for a deeper storage reservoir. The predominance of pre-eruptive shallow storage means that there is more chance of detecting reservoir destabilization with geophysical monitoring techniques compared to a scenario where melts are supplied from sub-Moho reservoirs.

Hawaii↗