USGS ScienceSearch

SEARCH · USGS Science

Results for “Annual Review of Environment and Resources”

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,696 records · Page 67Linked to original sources

Helium detection as a guide for uranium exploration

Helium, a byproduct of radioactive decay, may prove to be a valuable indicator of the presence and distribution of uranium deposits. Recent technological advances permit the development of instrumentation not previously adapted for this purpose. Commercially available equipment can provide high sensitivity at low cost and allow reasonable mobility for field use. A truck-mounted mass-spectrometer, tuned for helium-4, permits immediate adjustment or modification of sampling patterns in response to accumulating data. The inlet system of the spectrometer has been designed to allow flexibility in gas analyses from various sample types--soil gas, atmosphere, or gases in water. Sensitivity of the instrument is better than 50 parts of helium per 10 9 parts of gas. Replicate samples and standards can be analyzed in only 3 minutes. Helium in soil gas is being used for the initial evaluation of the technique. A hollow probe, as much as 2 m long, is driven into the ground; and a 10-cc syringe is used to purge the probe and extract the gas. This sample is then injected into the inlet system and introduced into the spectrometer at constant pressure. The output signal is displayed on a chart recorder. Investigations can be performed on a qualitative, relative basis, or on a quantitative basis by comparison to calibrated helium standards. Preliminary field testing includes studies of the responses to variations in wind speed, temperature, barometric pressure, moisture, and sampling depth over extended time periods, as well as studies of geologic controls on the helium content in soil gas. Surveys over known uranium occurrences reveal some anomalous helium distributions. This report describes an updated technological approach to an old idea: that of using helium, a byproduct of uranium radioactive decay, as an exploration tool for uranium. Helium is the sixth most abundant gas found in the earth's atmosphere (Table 1). It is an inert gas and very mobile, mixing rapidly. When the alpha particles from the radioactive decay of uranium and thorium pick up two electrons, atoms of the isotope helium-4 are formed. Several helium-4 atoms are produced from each decay series of U-238, U-235, and Th-232 (fig. 1). The half-lives of the parent isotopes are5also shown in figure 1. One gram of uranium will produce ~10 5 atoms of helium-4 per second, and one gram of thorium produces ~2.5 x 10 4 atoms of helium-4 per second (fig. 2). Calculations of the helium produced from the crust and mantle reveal that 1,125 x 10 30 atoms per year are produced, but only 7 x 10 30 atoms per year escape from the earth (Damon and Kulp, 1958). More is being produced than is being lost; in fact, the total atmospheric content of helium could be produced in only 2 million years (MacDonald, 1963)! However, all the crustal and mantle helium-4 does not degas into the atmosphere as it is formed. It is trapped in crystal lattices and in pore spaces within the earth. What this means is that there is excess helium-4 in the earth--excess, that is, compared to the atmospheric concentration in equilibrium with the helium-4 escape rate into space and the helium-4 flux from the crust and mantle (Nicolet, 1957). Near-surface pockets of high helium-4 concentrations are known, and some natural gas fields are so enriched that they are the source for commercial production of helium. Work by Roberts and others (1975) has shown high helium concentration in the soil gas associated with geothermal areas; and work by Clark and Kugler (1973), Dyck (1975), and Goldak (1974), for example, has noted high helium concentrations in soils and waters in the vicinity of uranium deposits. These latter studies certainly indicate the potential of helium detection for use as an exploration tool for uranium.

Open-File Report

Complex carbonate ore mineralogy in the Mountain Pass carbonatite rare earth element deposit, USA

Economic concentrations of rare earth element (REE) minerals are uncommon in the Earth’s crust, with most occurring in carbonatites. Unlike most igneous rocks composed of silicate minerals, carbonatites are dominated by carbonate minerals, some of which can incorporate significant light REEs (LREEs; La, Ce, Pr, Nd). Technological applications of REEs are numerous and they have been identified as some of the most critical mineral commodities to the global economy. The Mountain Pass carbonatite stock in the Mojave Desert of California is the most economically significant REE deposit in the USA and contains a few to tens of percent (by volume) of the carbonate REE ore mineral bastnäsite. Despite the economic significance of the Mountain Pass deposit, studies of its ore mineralogy are limited. Here we present new carbonate ore mineralogy data for a compositionally diverse suite of carbonatitic rocks from the Mountain Pass stock and related dikes. Whole-rock geochemical data are integrated with mineral-scale textural and chemical data obtained by scanning electron microscopy (SEM), electron probe microanalysis (EPMA), and microRaman spectroscopy. Our results document a complex spectrum of REE-bearing carbonate minerals and intermediate mixed-layer structures. Mineral species include bastnäsite [REE(CO 3 )F], hydroxylbastnäsite [REE(CO 3 )OH], parisite [Ca(REE) 2 (CO 3 ) 3 F 2 ], synchysite [Ca(REE)(CO 3 ) 2 F], röntgenite [Ca 2 (Ce,La) 3 (CO 3 ) 5 F 3 ], and sahamalite [(Mg,Fe 2+ )(REE) 2 (CO 3 ) 4 ]. Carbonate ore mineralogy is heterogeneous within and between samples, including at the intracrystal scale. Complexly zoned crystals exhibit as many as five to six different compositional domains and syntaxial intergrowths, commonly with the more Ca-rich varieties (parisite, synchysite) forming crystal rims that surround relict bastnäsite cores. We attribute the phenocryst variability to changes in the chemistry and temperature of primary carbonatite magmas and evolved/exsolved fluids. Cross-cutting vein textures of calcite, celestine and various REE carbonate minerals, interstitial bastnäsite crystallization, breccia blocks lined by fine-grained bastnäsite, and the presence of hydroxylbastnäsite and partially hydroxylated bastnäsite point to the role of secondary hydrothermal processes in REE mineralization. Fluorcarbonate mineral compositions demonstrate that La and Ce are more structurally abundant in bastnäsite, whereas the more Ca-rich species (parisite, synchysite) contain a greater proportion of REE heavier than Pr (Nd, Sm, Eu, Gd) and Y. Atomic ratios of Pr/(Nd + Pr) are likewise variable, with the highest average value for bastnäsite (0.25) compared to parisite (0.22) and sychysite (0.21). This finding has geometallurgical implications, given that current mining operations are focused on recovery of Nd and Pr for high field strength permanent magnets and the Nd/Pr ratios are a critical factor in ore processing and magnet manufacture.

California

Geologic map of the Buckner 7.5-Minute quadrangle, Louisa County, Virginia

The Buckner 7.5-minute quadrangle straddles three terrane boundaries in the Piedmont Physiographic Province in central Virginia: the Chopawamsic terrane, the Elk Hill Complex, and the Goochland terrane. In much of the map area, the Elk Hill Complex separates the Chopawamsic and Goochland terranes. Rocks of the Chopawamsic terrane include Ordovician metavolcanic and metasedimentary rocks of the Chopawamsic Formation, Ordovician to Silurian granodiorite sheet intrusions, and Paleozoic mafic intrusions. Silurian to Devonian rocks of the Quantico Formation, mostly garnet-mica schist, crop out in the northwesternmost part of the map area, and are in unconformable contact with rocks of the Chopawamsic Formation on the southeastern limb of the Quantico synclinorium. The main map unit in the Elk Hill Complex is Neoproterozoic mica gneiss, which is in pre-metamorphic fault contact with rocks of the Chopawamsic Formation to the west. The main map unit of the Goochland terrane is the Maidens Gneiss. Except for Jurassic diabase dikes, all rocks on the Buckner 7.5-minute quadrangle were metamorphosed to amphibolite facies during the Alleghanian orogeny and preserve multiple compositional and phyllosilicate penetrative foliations. Evidence of amphibolite-facies metamorphism during the Taconic orogeny is preserved in rocks of the Elk Hill Complex. The entire width of the Maidens Gneiss on the Buckner 7.5-minute quadrangle is within the Spotsylvania high-strain zone and amphibolite-facies mylonitic textures are pervasive. Quartz veins and Jurassic diabase dikes crosscut all older rocks of the quadrangle. Multiple levels of terrace deposits are present along and near the major streams of the quadrangle. The lower terrace deposits are likely remnants of former positions of the Little River on the landscape, whereas higher deposits may be remnants of former deposits of the Atlantic Coastal Plain that covered this portion of the Piedmont Province. A linear cluster of aftershocks from the magnitude 5.8 earthquake that occurred near Mineral, Virginia, in 2011 defines the Fredericks Hall fault, which is at depth on the Buckner quadrangle. Most of the aftershocks occurred in the core of the Elk Creek antiform and have no relation to faults mapped at the surface. Several abandoned crushed stone and building stone quarries, as well as a mica prospect, exist in the quadrangle.

Virginia

Slope Unit Maker (SUMak): An efficient and parameter-free algorithm for delineating slope units to improve landslide modeling

Slope units are terrain partitions bounded by drainage and divide lines. In landslide modeling, including susceptibility modeling and event-specific modeling of landslide occurrence, slope units provide several advantages over gridded units, such as better capturing terrain geometry, improved incorporation of geospatial landslide-occurrence data in different formats (e.g., point and polygon), and better accommodating the varying data accuracy and precision in landslide inventories. However, the use of slope units in regional ( > 100 km 2 ) landslide studies remains limited due, in part, to the large computational costs and/or poor reproducibility with current delineation methods. We introduce a computationally efficient algorithm for the parameter-free delineation of slope units that leverages tools from within TauDEM and GRASS, using an R interface. The algorithm uses geomorphic laws to define the appropriate scaling of the slope units representative of hillslope processes, avoiding the often ambiguous determination of slope unit size. We then demonstrate how slope units enable more robust regional-scale landslide susceptibility and event-specific landslide occurrence maps.

Natural Hazards and Earth Systems Sciences (NHESS)

Using public participatory geographic information systems (PPGIS) to explore uses and values for Mojave Trails National Monument, California

Many people ascribe a variety of values to public lands and waters, but some values are more difficult to assess and quantify than others. Public participatory geographic information systems (PPGIS) are tools that have been used to help quantify and map the public’s diverse values for a landscape. This work describes the first known Office of Management and Budget–approved use of PPGIS by a Department of the Interior bureau. The U.S. Geological Survey developed an internet-based application to aid in gathering PPGIS data, called Values Mapping for Planning in Regional Ecosystems (VaMPIRE). Further, this work describes the first pilot of the VaMPIRE application in coordination with the Bureau of Land Management to collect spatial data and other survey data regarding the public’s uses of and values for locations within Mojave Trails National Monument. We emailed the link to the VaMPIRE application to an interested party email list in 2024 with 207 valid emails and received 74 responses; we also received 47 responses from members of an off-roading social media group. Of the list of 16 value options, recreation was the most popular value for the monument, followed by wilderness and inspirational. Over 1,000 points were placed throughout the monument, indicating locations people use or value, with the locations spread throughout the entire monument. Additionally, most survey respondents stated their ability to receive benefits in locations they mapped would not change in response to a hypothetical scenario related to recreational facility development. This report describes exploratory results from the first use of the VaMPIRE tool in Mojave Trails National Monument and includes reflections on how the process went and considerations for future use of VaMPIRE.

California

3-D geological modeling for numerical flow simulation studies of gas hydrate reservoirs at the Kuparuk State 7-11-12 Pad in the Prudhoe Bay Unit on the Alaska North Slope

Accurate reservoir evaluation requires reliable three-dimensional (3-D) geological models. This study conducted 3-D geological modeling for numerical flow simulation of the B1 sand gas hydrate reservoir at the Kuparuk State 7-11-12 pad, Prudhoe Bay Unit, Alaska North Slope. The model integrates well logs, core, and seismic data to address spatial heterogeneity in geological structures and reservoir properties. Two modeling types were performed: structural framework modeling and petrophysical property modeling. For structural framework modeling, seismic data and well log markers were used to reproduce subsurface structures characterized by a normal fault system. A volume-based modeling algorithm and stair-step gridding were applied. The resulting 3-D model comprised 2,640,000 grid cells across 264 layers, including seven fault grids. For petrophysical property modeling, total porosity was initially modeled using sequential Gaussian simulation with collocated cokriging. To reproduce the upward coarsening of the B1 sand, upscaled log-derived total porosity and a 3-D trend depicting total porosity variation were used as primary and secondary data, respectively. Gas hydrate saturation distribution was modeled similarly, with secondary data from estimated porosity distribution and seismic-derived acoustic impedance map enhancing accuracy. Results indicate higher gas hydrate saturation in the upper part of the B1 sand and areas with higher acoustic impedance. Intrinsic permeability was modeled from the total porosity and clay-bound water volume, and effective permeability was derived from the gas hydrate saturation and intrinsic permeability distributions based on the “Tokyo model”. Effective permeability distributions were influenced by the total porosity, gas hydrate saturation, and intrinsic permeability. Within the same layer, higher gas hydrate saturation leads to decreased effective permeability. In total, 100 sets of multiple scenarios were prepared, providing input data for dynamic flow simulations to evaluate the effects of lateral heterogeneity in reservoir properties and the hydraulic characteristics of faults on production behavior for preassessment before the long-term production test.

Alaska

Vulnerability assessment of groundwater influenced ecosystems in the Northeastern United States

Groundwater-influenced ecosystems (GIEs) are increasingly vulnerable due to groundwater extraction, land-use practices, and climate change. These ecosystems receive groundwater inflow as a portion of their baseflow or water budget, which can maintain water levels, water temperature, and chemistry necessary to sustain the biodiversity that they support. In some systems (e.g., springs, seeps, fens), this connection with groundwater is central to the system’s integrity and persistence. Groundwater management decisions for human use often do not consider the ecological effects of those actions on GIEs. This disparity can be attributed, in part, to a lack of information regarding the physical relationships these systems have with the surrounding landscape and climate, which may influence the environmental conditions and associated biodiversity. We estimate the vulnerability of areas predicted to be highly suitable for the presence of GIEs based on watershed (U.S. Geological Survey Hydrologic Unit Code 12 watersheds: 24–100 km 2 ) and pixel (30 m × 30 m pixels) resolution in the Atlantic Highlands and Mixed Wood Plains EPA Level II Ecoregions in the northeastern United States. We represent vulnerability with variables describing adaptive capacity (topographic wetness index, hydric soil, physiographic diversity), exposure (climatic niche), and sensitivity (aquatic barriers, proportion urbanized or agriculture). Vulnerability scores indicate that ~26% of GIEs were within 30 m of areas with moderate vulnerability. Within these GIEs, climate exposure is an important contributor to vulnerability of 40% of the areas, followed by land use (19%, agriculture or urbanized). There are few areas predicted to be suitable for GIEs that are also predicted to be highly vulnerable, and of those, climate exposure is the most important contributor to their vulnerability. Persistence of GIEs in the northeastern United States may be challenged as changes in the amount and timing of precipitation and increasing air temperatures attributed to climate change affect the groundwater that sustains these systems.

Connecticut, Delaware, District of Columbia, Maine

Potential corrosivity of untreated groundwater in Louisiana

Corrosive groundwater can cause lead, copper, and other metals to leach from pipes and plumbing fixtures in water distribution systems. Metals, if ingested, could lead to serious health implications to the nearly 2.9 million people in Louisiana who obtain their drinking water from groundwater sources. Four indices—the Langelier Saturation Index (LSI), Ryznar Stability Index (RSI), Puckorius Scaling Index (PSI), and the Potential to Promote Galvanic Corrosion (PPGC)—in addition to an analysis which normalized the results from the existing indices, the Combined Index (CI), were used to assess the corrosivity of groundwater in Louisiana and identify areas within eight major aquifers and aquifer systems with moderate to high corrosivity potential. The purpose of this study is to provide State and local governments, public water system managers, and the nearly 500,000 private well owners in Louisiana with information needed to manage drinking-water supplies and mitigate potential health risks related to leaching of metals from water pipes and fixtures. The average scores of untreated groundwater samples from approximately 375 wells by index are as follows: LSI, −1.28; RSI, 9.78; PSI, 9.34; and CI, 4.14. The PPGC does not produce a numerical score, but the total percentage of class counts can be used to assign a classification; overall, samples in Louisiana were classified as significant concern. The percentages of groundwater samples from wells classified as potentially corrosive, by index, are as follows: LSI, 53 percent; RSI, 94 percent; PSI, 81 percent; PPGC, 98 percent; and CI, 81 percent. The percentages of samples classified as indeterminate, by index, are as follows: LSI, 46 percent; RSI, 5 percent; PSI, 12 percent; PPGC, 0 percent; and CI, 18 percent.

Louisiana

Syn-magmatic subsidence during the early stages of continental rifting in the Mesoproterozoic—A reanalysis of legacy data for the Midcontinent Rift, western Lake Superior

The Midcontinent Rift system (ca. 1.1 Ga) is a 2000-km-long series of elongated volcanic and sedimentary troughs and associated intrusive centers exposed chiefly in the Lake Superior region of North America. The rift system represents a long history of intense magmatism and subsequent sedimentation that was arrested by far-field tectonic events before sea-floor spreading was established. The premature cessation preserved a record of processes related to the beginning of continental rifting. The rift system under Lake Superior has been long studied using seismic-reflection data collected as part of the Great Lakes International Multidisciplinary Program on Crustal Evolution (GLIMPCE). We reexamine GLIMPCE Line C by developing a detailed velocity model for time to depth conversion constrained by other legacy data. We corroborate the model and develop a geologic interpretation using gravity and magnetic modeling and ties to geology mapped onshore. We recognize superposed subsiding sedimentary and volcanic basins for the southern half of the Line C depth section. This interpretation differs from previous paradigms that show major crustal faults that bound half-grabens or full grabens. We conclude that high-velocity (6.9 km/s) intrusive zones rather than major crustal faults border the sides of the basins. We speculate that the volcanic basin represents the initiation of seaward dipping reflectors. The syn-magmatic subsidence can be explained by dike injection and volcanic loading. Discrete lava basins throughout the region likely subsided at different times in a disorganized manner along the rift trend, raising questions about the long-term role of lithospheric thinning and melt generation.

Michigan, Minnesota, Ontario, Wisconsin

Evolution of arsenic speciation during thermal treatment simulating wildfire heating in arsenic-rich sediments

Understanding arsenic transformations during wildfires can help better constrain the environmental impacts of increased wildfire intensity, frequency, and burned area. Previous studies have monitored the evolution of metal(loid) speciation at specific time/temperature endpoints, hampering the comprehension of the dynamic transformations of metal(loid)s during wildfires. Here, a novel approach, namely in situ time-resolved X-ray absorption near edge structure (TR-XANES) spectroscopy, was used for the first time to monitor the heat-induced redox transformations of arsenic (As) in two As-rich sediments with different characteristics. No arsenic transformations were observed in either sediment at room temperature or when the sediments were heated to a maximum temperature of 400 ºC. When heated to a maximum temperature of 600 or 700 ºC, As underwent complex, dynamic, and partially reversible redox transformations, the extent of which varied with sediment properties, initial As speciation, and heating temperature and duration. In the case of sediments initially containing As(V) and a low sulfide/sulfate ratio, the As(V) was reduced to As(III) over a short period of time, followed by immediate reoxidation of As(III) to As(V). The extent of the transient As reduction increased with increasing maximum temperature. The final As speciation at the end of the heating experiment was not substantially different from that at the beginning of the experiment. In contrast, in sediments that were best described by a more complex combination of reference compounds (As(V) coprecipitated with or adsorbed on iron oxide, arseniosiderite, arsenopyrite) with a high sulfide/sulfate ratio, As(V) phases were reduced to As(III) and As(III) was further reduced to As(-I) (at 700 ºC, only), followed by delayed reoxidation of As(-I) to As(III) (at 700 ºC, only) and reoxidation of As(III) to As(V). Approximately 9% and 26% of all As remained as As(III) ( i.e., was not reoxidized to As(V)) at the end of the experiment when heated to a maximum temperature of 600 and 700 ºC, respectively, which is higher than the relative abundance of As(III) in the native sediment. These fire-induced As transformations can potentially enhance As mobilization to surface water and groundwater, posing risks to environmental and human health.

California

The petrogenesis of Þingmúli volcano, East Fjords, Iceland

In this work we revisit Þingmúli volcano (Þ = Th), a classic locality known as an example of a complete tholeiitic differentiation. Þingmúli is a ~ 9.5 Ma extinct central volcano located in the East Fjords of Iceland, in which the whole compositional spectrum from basalt to rhyolites have erupted. These volcanic products have been previously considered as petrogenetically related by an ideal fractionation trend, regardless any temporal relationship or volumetric considerations. Here we report new whole-rock geochemistry, mineral chemistry, isotope analyses, estimation of residence times of the different eruptive deposits, and an update of the original petrogenetic model. Our results highlight that an enriched source, likely spinel lherzolites, generated transitional-alkaline basaltic melts after 15–20% of partial melting at depths of 40–45 km. Many of these basaltic melts erupted at various stages of the volcano's history, while others remained longer in the volcanic plumbing system. These evolved by fractional crystallisation into basaltic andesite magmas with a residence time of ~5 years based on the crystal size distribution of the plagioclase population. Isotopic differences between the basalts/basaltic andesites ( 87 Sr/ 86 Sr ~ 0.7034; 143 Nd/ 144 Nd ~ 0.51315) and the erupted rhyolites ( 87 Sr/ 86 Sr ~ 0.7037; 143 Nd/ 144 Nd ~ 0.51304) indicate that the latter are not petrogenetically related to the former. Therefore, instead of a fractional crystallisation mechanism to generate the rhyolites, we propose the partial melting of ignimbrite layers located beneath the volcano. The broad range of trace element concentrations in andesites and dacites and their different isotopic values compared to the basalts strongly suggest that these magmas have been generated by magma mixing between basaltic and rhyolitic melts, similar to modern day Icelandic volcanoes such as Hekla. These results highlight the need to revisit previously studied Icelandic classic localities and reassess their traditionally proposed petrogenetic models.

Þingmúli volcano, East Fjords

Estimating groundwater level records using MOVE.1 and computing monthly percentiles from estimated groundwater records in Massachusetts

The U.S. Geological Survey, in cooperation with the Massachusetts Department of Environmental Protection, performed record extensions on groundwater levels at select wells using the Maintenance of Variance Extension type 1 (MOVE.1) method. The groundwater levels estimated from these record extensions were used to compute monthly percentiles to improve future determinations of a groundwater index. In Massachusetts, 27 of 29 short-record study wells with continuous groundwater levels between 0.8 and 8.1 years were suitable for record extensions; 37 long-record index wells were used to extend the groundwater level records at the study wells. The index well selected to pair with a study well was chosen based on Pearson correlation coefficient values; cross-correlation between the two wells; geologic and topographic similarity; and smallest distance spanning the wells. Each study well and its corresponding index well have 1 or more years of concurrent, overlapping data; a Pearson correlation coefficient that exceeded a threshold value of 0.8; and a similar aquifer type and hydrologic characteristics. Of the 29 study wells, 2 showed poor correlations with all index wells and were not considered for record extensions. Performance metrics used to assess the accuracy of the MOVE.1 models indicated that most models provided reasonable estimates of groundwater levels. Root mean square error values ranged from 0.097 to 2.292 feet, with a median of 0.536 foot. Nash-Sutcliffe efficiency coefficient values ranged from 0.623 to 0.996, with a median value of 0.759. Generally, study wells in close geographical proximity to their index well resulted in stronger model performance. The average length of groundwater level records was extended by 14.1 years to a new average of 18.1 years. The estimated groundwater level records from the MOVE.1 models resulted in an increase in the range of highest and lowest groundwater levels at 23 of 27 wells. The increase in range of groundwater levels was between 0.08 to 7.95 feet. Monthly percentiles for State drought indices were computed from the estimated MOVE.1 records and observed records through December 31, 2021. Percentiles computed from estimated records show an average groundwater level about 1.0 foot lower than observed data at the 2d percentile and 0.1 foot lower at the 30th percentile.

Massachusetts

Potentially toxic elements in wild Agassiz’s desert tortoises: Tissue concentrations and association with disease

Background: Desert tortoise ( Gopherus agassizii ) populations have continued to decline due to infectious and other diseases, predation, and habitat alteration. The potential contribution of minerals and heavy metals to tortoise health and susceptibility to disease remains uncertain. Objective: The objective of this study was to evaluate the results of elemental analysis of trace minerals and macrominerals in scute keratin, kidney, and liver tissue from ill and dying desert tortoises salvaged for necropsy between 1993 and 2000. Methods: Salvaged tortoises were categorized by age (adult, juvenile), geographic location, and primary disease based on necropsy findings. A subset of tortoises that were injured or killed by vehicular trauma or predation but with no notable pathologic abnormalities were used as controls. A panel of 21 trace minerals and 6 macrominerals was analyzed in scute keratin, kidney, and liver tissue samples by inductively-coupled plasma spectrometry and atomic absorption spectrophotometry. Results: Necropsies were done on 46 tortoises, including 9 juveniles salvaged from 5 regions in the Colorado and Mojave deserts of California. Primary diseases were cutaneous dyskeratosis (n=9), infection/ inflammation (n=8), malnutrition (n=7), mycoplasmosis (n=5), and urolithiasis (n=3); 14 tortoises were classified as controls. Concentrations of elements differed significantly by tissue, age, desert region, and disease (P < 0.05). Tortoises with cutaneous dyskeratosis had significantly higher Se concentrations, primarily in keratin and liver, than tortoises with other diseases (P < 0.001). Juveniles were more likely than adults to have high Pb, Sn, and Zn levels (P < 0.05). All tortoises had detectable levels of more than one potentially toxic heavy metal, including As, Cd, Cr, Hg, Ni, Pb, Sn and V. Conclusions: Potentially toxic elements are frequently found in tissue from tortoises in desert regions of California, with significantly higher concentrations in diseased tortoises. Metal exposure from soils, mining, historic and ongoing military activities, and other human activities could increase susceptibility to disease in desert tortoises.

California

The 17 January 1994 Northridge, California, earthquake: A retrospective analysis

The 17 January 1994 Northridge, California, earthquake was a watershed event, with far-reaching societal and scientific impacts. The earthquake, which occurred in the early days of both broadband seismic networks and the Internet, spurred advances in seismic monitoring, real-time systems, and development of data products. Motivated by the 30th anniversary of the earthquake, we present a brief retrospective of the earthquake and its impact, and reconsider both ground motions and the aftershock distribution using modern tools and best-available data. With improvements in instrumentation and analysis methodology, recent earthquakes continue to reveal increasing complexity of ground motions, fault systems, and earthquake ruptures. Even in the absence of data from state-of-the art instrumentation, a retrospective consideration of ground motion data from the Northridge earthquake reveals complexities beyond what could be characterized (and modeled) thirty years ago. Aftershock relocations for both the 1971 Sylmar and 1994 Northridge earthquakes also reveal an updated view of fault complexity. Our study does provide a cautionary tale regarding legacy data sets and research results that are not easily accessible, which can result in discrepancies between catalog data and products from best-available science. We also briefly describe outreach products produced as part of the anniversary commemoration.

California

Ground motions from finite faults with unknown geometry

Earthquake emergency response and loss modeling require accurate models of ground‐motion intensity metrics (IMs), such as those produced by ShakeMap. In the critical minutes and hours after an earthquake, however, knowledge of the earthquake source geometry limits the accuracy of these models, especially in areas where seismic instrumentation is insufficient to constrain the shaking intensity pattern. Previous attempts to overcome this limitation have centered on computing the mean distance metrics from a suite of assumed source geometries. That approach is inadequate with modern ground‐motion models (GMMs) because they typically include terms for the hanging wall effect, depth to top of rupture, and other nonlinear terms. Such approaches also propagate the distance uncertainty through the GMMs and may not properly allocate the additional uncertainty into the between‐ and within‐event components. Here we present Finite‐Fault SIMulation with Mixed Effects Regression (FFSIMMER), a method to address these limitations by computing IM distributions from a suite of rupture realizations that fully accounts for all geometry‐related GMM terms and allows for the separation of the between‐ and within‐event components of the total uncertainty. The accuracy of this method is limited by the available constraints on the source parameters. To address this, FFSIMMER can incorporate constraints based on knowledge of the source, such as the approximate strike and dip based on knowledge of the plate tectonics and focal mechanism of the earthquake. We demonstrate the improvement in accuracy that can be achieved with the 2025 M 7.7 Mandalay, Burma (Myanmar), earthquake by incorporating constraints on the range of plausible source parameters.

Bulletin of the Seismological Society of America

REDPy: A Python tool for automated repeating earthquake detection and visualization

Detecting and cataloging seismic events are among the most fundamental tasks in seismology. Many standardized tools for these tasks exist, including the open‐source package repeating earthquake detector in Python (REDPy). REDPy generates an organized catalog of seismic events from continuous waveform data, in which events are automatically separated into groups (“families”) by their waveform similarity through cross‐correlation. REDPy also automatically generates various outputs that allow a user to visualize important trends in the catalog, which may be used in real time or in retrospective analyses to allow rapid identification of interesting features. The code was designed for near‐real‐time volcano monitoring but is applicable across a broad range of use cases in seismology and seismoacoustics. In this article, the utility and performance of REDPy are demonstrated on two highly seismogenic volcanic eruption sequences: the onset of the dome‐building eruption of Mount St. Helens, Washington, from 2004 to 2005, and the entirety of the summit caldera collapse sequence of Kīlauea, Hawai‘i, in 2018. This article is meant to be a companion to the documentation of the code; in addition to detailing the basic required inputs, script functionality, and resulting outputs, the reasonings behind several important design decisions are also discussed.

Seismological Research Letters

A glimpse into the future of tectonic tremor monitoring

Tectonic tremor is a weak, long-duration seismic signal often observed in subduction zones and on some other plate-bounding faults. Because of tremor's characteristically low amplitude (and low signal-to-noise) and lack of clear phase arrivals, detecting and locating tremor usually requires techniques distinct from those applied to typical earthquakes. Major advances in detection and understanding of tremor have derived in the past from a powerful combination of new data and new analysis techniques. In a recent study, Sagae et al. (2025, https://doi.org/10.1029/2025jb031348 ) exploit that combination again, developing a new machine-learning based workflow and applying it to the S-net cabled seismic network in the Japan trench offshore northern Honshu. Their approach, although complex, succeeds in detecting several times more tremor activity than earlier studies, resulting in new insights and providing a blueprint for similar approaches that could be applied elsewhere. As real-time earthquake monitoring adopts similar tools, it may present an opportunity to bring tremor monitoring into operational workflows. In turn, this could solidify tremor monitoring as a component of future operational earthquake forecasting.

JGR Solid Earth

Characterization and validation of tidally calibrated strains from the Alto Tiberina Near Fault Observatory Strainmeter Array (TABOO-NFO-STAR)

Six horizontal borehole tensor strainmeters (TSM1-6) installed from Fall 2021 to Spring 2022 comprise the Alto Tiberina Near Fault Observatory Strainmeter Array (STAR), providing an unprecedented opportunity to investigate seismic and aseismic deformation from hazardous high- and low-angle normal faults in Italy. Prior to use in tectonic applications, they require in-situ calibration and correction for non-tectonic signals. We tidally calibrate the instruments, characterize the calibration uncertainty, and test the results against environmental and earthquake signals originating from local to teleseismic distances. The STAR sites demonstrably deviate from assumptions common to the standard manufacturer's calibrations, including negative areal coupling at TSM3-6. While the tidally calibrated strains have ~3-56% uncertainty, the calibrated dynamic strains show interstation precision and accuracy to nanostrain levels, and static coseismic offsets in the array footprint are within uncertainty. TSM3 records a complex series of strains that may arise from dynamically triggered near-borehole fracture slip and fluid flow that does not appear to affect its sensitivity to lower strain rate deformation. Future calibration improvement may be afforded with longer stable timeseries, particularly for TSM4. Overall, our analyses demonstrate expanded geodetic capability for detecting deformation in the Alto Tiberina Near Fault Observatory.

Alto Tiberina Near Fault Observatory