SEARCH · Geology Reports

Results for “Journal of Agricultural, Biological and Environmental Statistics”

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

SourceResearch connectionFrom yearThrough year

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

1,689 records · Page 85Linked to original sources

USGS2026

Importance of fish in the diet of Adélie penguins across multiple life stages

Over recent decades, the Adélie penguin ( Pygoscelis adeliae ) population has grown across most of its range, the exception being the northern coast of the western Antarctic Peninsula (WAP). In the Ross Sea, the very large Cape Crozier colony grew to reach ~9% of the global population. Demographic factors driving the Cape Crozier trend have yet to be identified. However, low chick fledging mass - a negative influence in the northern WAP where the diet is mostly less-energy-dense krill - is not showing an effect in the Cape Crozier colony. Previously, we hypothesized that Cape Crozier fledglings, and post-breeding adults, must be finding sufficient higher-quality prey once free of the prey-depleted colony foraging area, with subsequent higher survival. Here we explore penguin diet within and outside the colony foraging area using stable isotope analysis (SIA) of feathers: those grown by post-breeding adults outside the foraging area, prior to moult, compared to those of near-fledged chicks raised on meals obtained within the foraging area. Second, we explore whether a major difference occurs between breeding and post-breeding diet of adults from a small, nearby colony (Cape Royds) that exploits a smaller, little-depleted preyscape. Finally, we compare the results to SIA analyses published by others to explore whether the pre-moult diet of Ross Island adults, mostly foraging within the north-eastern Ross Sea, compares with that of adults nesting in the southern WAP (Bellingshausen Sea coast). At the latter, populations are not decreasing, with diet containing an appreciable contribution of fish. The results confirmed that fish contributed greatly to chick diets at both Ross Island colonies (δ 15 N 10.6–11.4). Once adults were released from central-place foraging, SIA levels (δ 15 N 9.7–12.9) were similar to those of the southern WAP. Foraging on energy-dense fish in coastal Antarctic waters could be one reason Adélie penguins are among the most abundant and therefore most ecologically successful penguin species.

Antarctic Science↗

USGS2026

Corundum discovered by SuperCam and the Perseverance rover at Jezero crater, Mars

Mars is primarily composed of mafic mineral assemblages and their alteration products, but small, scattered rocks strewn across the landscape offer clues to greater petrological diversity. While traversing the Jezero crater rim, the Perseverance rover encountered several plagioclase-rich light-toned float rocks. SuperCam identified the distinctive signature of corundum (α-Al 2 O 3 ) in these rocks using time-resolved luminescence spectroscopy. Two strong peaks (692.7 and 694.1 nm) with millisecond lifetimes, and additional supporting lines, are consistent with Cr 3+ substitution for Al 3+ in corundum. Corundum forms in Al-rich, Si-depleted environments through magmatic or metamorphic processes. Given the rocks' small size, association with plagioclase, and location on the crater rim, we interpret the most plausible formation scenario to be impact induced metamorphism at the interface of a felsic and a mafic/ultramafic member with the likely action of fluids at some stage, although other possibilities are not excluded.

Geophysical Research Letters↗

USGS2025

Rare earth elements on the Moon

Rare earth elements (REEs) are a scarce but vital resource for our modern economies and lifestyles. Since the late 1990s, China has supplied the vast majority of the world’s refined REEs. Increasing global demand has broadened the search for REE deposits to unconventional places, including the Moon. Although most lunar rocks have very low REE concentrations, Apollo samples showed that one type of lunar rock containing potassium (K), REEs, and phosphorus (P)—known by the acronym KREEP—has high concentrations of REEs. Data from orbiting satellites have identified locations where substantial deposits of KREEP are likely. The viability of mining these deposits depends on the evolution of REE economics, the development of the Earth-Moon infrastructure, and the findings from future lunar mineral exploration missions.

Fact Sheet↗

USGS2025

2024 Surprise Inlet landslides: Insights from a prototype landslide‐triggered tsunami monitoring system in Prince William Sound, Alaska

Alaska's coastal communities face growing landslide hazards owing to glacier retreat and extreme weather intensified by the warming climate, yet hazard monitoring remains challenging. As part of ongoing experimental monitoring in Prince William Sound, we detected three large landslides (0.5–2.3 M m 3 ) at Surprise Inlet on 20 September 2024, within the span of an hour. These events were identified in near real-time through seismic data and later confirmed using satellite imagery, tidal records, and infrasound. The landslides generated a modest tsunami, and a 4 cm wave was recorded by a tide gauge 18 km away, marking the first recorded landslide to reach water since monitoring began in this region in 2021. Here, we examine the detection and interpretation of these landslides using multiple data sources and modeling. We demonstrate the effectiveness of this regional seismic monitoring system and show how complementary instrumentation, where available, can enhance detection capabilities.

Alaska↗

USGS2024

Groundwater flow model for the Des Moines River alluvial aquifer near Des Moines, Iowa

Des Moines Water Works (DMWW) is a regional municipal water utility that provides residential and commercial water resources to about 600,000 customers in Des Moines, Iowa, and surrounding municipalities in central Iowa. DMWW has identified a need for increased water supply and is exploring the potential for expanding groundwater production capabilities in the Des Moines River alluvial aquifer, where it operates two radial collector wells (RCWs). The U.S. Geological Survey, in cooperation with DMWW, completed a study of the Des Moines River alluvial aquifer and interactions of the RCWs with the aquifer; no previously published model has included the existing well locations, which is the focus of this model. A conceptual and numerical groundwater flow model have been developed to characterize the Des Moines River alluvial aquifer under existing conditions, to simulate water levels observed in the RCWs, and to provide publicly accessible hydrologic data and research that advance understanding of the regional hydrologic system and can potentially be used in the future to evaluate groundwater production scenarios. Model performance was assessed by comparing observed and simulated groundwater levels that included water level elevations, water level changes, water level inequality observations, surface water streamflow, and change in surface water volume from upstream to downstream. Water table elevation in the aquifer layers is on average slightly overestimated with average absolute value error less than 1.5 meters at both RCWs and less than 2.5 meters for all observation wells in the alluvial aquifer layers. The model also accurately simulated water tables greater than the RCW design minimum (a water level threshold at which RCW pumping is reduced) in all timesteps for which water level observation data existed. Water table elevation error was higher in other model layers that were not the focus of the study, and the model did not accurately match streamflow targets.

Iowa↗

USGS2025

Favorability mapping for hydrothermal power resource assessments of the Great Basin, USA

The U.S. Geological Survey (USGS) is updating the 2008 assessment of conventional hydrothermal resources for the Great Basin in the western United States. As part of this work, the workflow for hydrothermal resource favorability maps is being modified to integrate modern data-driven machine learning (ML) methods. Improvements include: [1] using new and refined evidence layers (features); [2] using an order of magnitude more training sites (labeled examples); [3] utilizing simple but non-linear supervised ML algorithms; [4] representing positive training sites (wells with measured heat flow) with their ordinal value proportional to the magnitude of convective upflow (i.e., low, high, or very high convective signals instead of past strategies using positive-negative labels); [5] supplementing training sites with additional sites with low convective signals to represent diverse under-sampled areas where hydrothermal systems are unlikely to exist; [6] comparing with competing approaches; and [7] utilizing Monte Carlo cross-validation to estimate and evaluate prediction uncertainty. For the new favorability map, over half of the power-producing systems (i.e., 15 of 28) are predicted in the 99th percentile of most favorable locations (i.e., the highest 1 % of favorability, corresponding to 1 % of the map area), exceeding the performance of past models that have explicitly used power plants as training sites. Previous favorability maps predicted approximately half of the power-producing hydrothermal systems above the 80th percentile (i.e., 20 % of the map area). For the new favorability map, 93 % of power-producing systems (i.e., 26 of 28) are above the 80th percentile. The power-producing systems for which the new model does not perform well are either comparatively small, low-temperature systems or systems also not predicted well by prior modeling approaches, suggesting that these few systems are unusual when compared with most power-producing systems. Focusing research on these known, seemingly different systems may yield new insights and subsequent discovery of new prospects.

California, Idaho, Nevada, Oregon, Utah↗

USGS2025

U.S. Geological Survey Groundwater Climate Response Network, 2024

As of October 2024, the U.S. Geological Survey (USGS) operated 588 sites across the United States and its territories as part of the Groundwater Climate Response Network (CRN). The CRN is comprised of wells selected to monitor the effects of climate variability, such as droughts, on groundwater levels nationwide. The CRN includes nearly 500 locations with real-time data and more than 100 sites with non-real-time data available to the public on the CRN web mapper and the USGS National Water Dashboard.

General Information Product↗

USGS2025

An unexplained tsunami: Was there megathrust slip during the 2020 Mw7.6 Sand Point, Alaska, earthquake?

On October 19, 2020, the M w 7.6 Sand Point earthquake struck south of the Shumagin Islands in Alaska. Moment tensors indicate the earthquake was primarily strike-slip, yet the event produced an enigmatic tsunami that was larger and more widespread than expected for an earthquake of that magnitude and mechanism. Using a suite of hydrodynamic, seismic, and geodetic modeling techniques, we explore plausible causes of the tsunami. We find that strike-slip models consistent with the moment tensor orientation cannot produce the observed tsunami. Hydrodynamic inversion of sea surface deformation from deep ocean and tide gauge data suggest seafloor deformation more closely matches a megathrust, rather than a strike-slip, source. Static slip inversions, using sea level and Global Navigation Satellite System data, allow for a portion of co-seismic megathrust slip that can explain tsunamigenesis. Combining all available geophysical datasets to model the kinematic rupture, we show that considerable, relatively slow, megathrust slip is allowable in the Shumagin segment, concurrent with strike-slip faulting. We hypothesize that the slow megathrust rupture does not contribute much seismic radiation allowing it to previously go unnoticed with traditional seismic monitoring.

Alaska↗

USGS2026

Cnidarian–algal partnerships structure bacterial communities during strobilation in Cassiopea xamachana

Cnidarian–algal (Symbiodiniaceae) symbioses rely on complex interactions among the cnidarian host, algal symbionts, and associated bacterial communities. In the upside-down jellyfish Cassiopea xamachana , the polyp-to-medusa transition (strobilation) requires the establishment of symbiosis with Symbiodiniaceae algal partners, yet bacterial community dynamics during this developmental process remain unknown. Here, we experimentally induced symbiosis in aposymbiotic polyps using four algal treatments: xenic Symbiodinium microadriaticum (native symbiont), xenic Breviolum minutum , antibiotic-treated B. minutum , and a photosynthetically impaired B. minutum mutant. We combined 16S rRNA gene sequencing with measurements of photosynthetic efficiency, asexual budding, and algal surface N-glycan profiles to characterize holobiont assembly during symbiosis onset and strobilation. Algal treatment structured bacterial communities in both algal cultures and polyp tissues. Our analyses identified a set of amplicon sequence variants that consistently distinguished strobilating polyps from non-strobilating aposymbiotic and mutant polyps, in addition to potential bacterial biomarkers associated with successful metamorphosis. Strobilation was associated with the enrichment of bacterial communities putatively involved in sulfur and nitrogen cycling, whereas non-strobilating aposymbiotic and mutant polyps were characterized by opportunistic bacteria and increased community variability. Together, these results reveal coordinated changes in algal physiology, surface glycan profiles, and bacterial community structure associated with successful strobilation in C. xamachana and support a model in which tripartite host–alga–bacteria interactions influence cnidarian life stage transitions.

ISME Communications↗

USGS2025

The Mammoth magnetic anomaly, Pinal County, Arizona

A high-resolution Earth Mapping Resources Initiative airborne geophysical survey was flown in the southwest North American porphyry copper province to improve bedrock geologic maps and to identify areas that have unrecognized critical mineral resource potential. During the review of the aeromagnetic data, a distinctly monopolar-shaped, negative magnetic anomaly was observed at a flight elevation of 200 m above the ground with a maximum amplitude of –9500 nT. We have named this the Mammoth magnetic anomaly (MMA) because it is centered 12 km northeast of the town of Mammoth, Arizona, USA. The total field anomaly (TFA) contour of –500 nT enclosing the MMA defines an elongate shape measuring 2.5 km long by 1 km wide that trends northwest–southeast. Given the striking nature of this negative, monopolar-shaped magnetic anomaly, we conducted a ground campaign in May 2025 to determine its authenticity and potential relationship to critical mineral endowment. The MMA was confirmed on the ground with a TFA approaching –46,000 nT. Total magnetic intensity (TMI) observations routinely fell below the 18,000 nT operating floor of an industry-standard cesium-vapor total field magnetometer, and extremely low TMI measurements were corroborated along coincident traverse lines using two high dynamic range, but lower sensitivity, smartphone vector magnetometers. The lowest TMI values recorded by both smartphone magnetometers were 1000 nT and confirmed with multiple adjacent and crossing lines. Field observations suggest that this magnetic feature is caused by strong remanent magnetization within fine-grained magnetite hosted within locally altered Pinal Schist.

Arizona↗

USGS2025

Improved prediction of postfire debris flows through rainfall anomaly maps

Predicting where runoff-generated debris flows might occur during rainfall on steep, recently burned terrain is challenging. Studies of mass-movement processes in unburned areas indicate that event locations are well-predicted by rainfall anomaly, R* , in which peak observed rainfall is normalized by local rainfall climatology. Here, we use remote and field methods to map debris flows triggered within the 2020 Dolan Fire burn area in coastal California, demonstrate that a short-duration R* metric predicts debris-flow occurrence more effectively than absolute peak intensity or longer-duration rainfall metrics, and show that incorporating an R* criterion into an existing debris-flow likelihood model can reduce false positive predictions and improve accuracy. We test R * at three other climatically distinct fires in California, demonstrating its utility for mapping likely debris-flow locations in different climates. We also consider how R* can benefit postfire debris-flow prediction given recent increases in climatological variability within individual burn perimeters.

Callifornia↗

USGS2025

Gulf Coast Basin CORE-CM initiative final report

The Bureau of Economic Geology at the University of Texas at Austin (UT-BEG) is leading the Gulf Coast Carbon Ore, Rare Earth, and Critical Minerals (CORE-CM) Initiative to assess the potential to produce critical minerals (CMs), including rare earth elements (REEs) from coal, coal ash, and produced water related to oil and gas production, and related materials (alumina processing waste [red mud], heavy mineral sands, graphite, and zeolite) within the Gulf Coast Basin. This project represents the first phase in a long-term program and provides reconnaissance data that will be foundational for future work by assessing resources and suggesting plans to be conducted in future work and expanding stakeholder engagement. The project includes several tasks designed to identify, characterize, and assess several necessary aspects for development of CMs and REEs in the Gulf Coast Basin.

Gulf Coast basin↗

USGS2026

Quantifying methane emissions from a rich fen with uncrewed aircraft systems in boreal Alaska

Thawing of permafrost in northern latitudes is accelerating, potentially releasing substantial amounts of methane (CH 4 ) as forested permafrost plateaus transition into wetlands. This ecosystem shift alters the carbon exchange between the soil and atmosphere, influencing the permafrost-carbon feedback. Monitoring these changes may require measurement platforms operating across varied spatial and temporal scales. Recent advancements in small uncrewed aircraft systems (sUAS) enable high resolution CH 4 flux quantification in remote, complex terrains; however, comparisons with established methods such as eddy covariance flux towers remain limited. We used a hexacopter sUAS to quantify CH 4 emissions from the Alaska Peatland Experiment, a wetland within the Bonanza Creek Experimental Forest. Using an ensemble of methods to define the background CH 4 concentration, along with near surface emissions from soil chambers, helped constrain our flux estimates. The sUAS method yielded an average flux of 0.0077 ± 0.0019 mol s −1 CH 4 , within a factor of two concurrent tower-derived total source flux estimates (0.0036 ± 0.00042 mol s −1 CH 4 ). To assess spatial drivers of observed fluxes, we conducted a 2D footprint analysis and overlaid the results with high-resolution hyperspectral land cover classification, quantifying vegetative contributions within each footprint. This revealed higher fen representation in sUAS measurements (73.8%) than in tower footprints (58.8%), and lower tussock meadow representation (15.6% and 30.3%, respectively). These differences were consistent with known variation in vegetation-specific CH 4 emissions. Our results highlight that combining footprint modeling with land cover characterization can enhance interpretations of CH 4 fluxes and guide cross-platform comparisons.

JGR Atmospheres↗

USGS2026

The digital archivist: Automating legacy macroseismic data processing using large language models

Macroseismic data are a key resource to investigate shaking and damage from preinstrumental and early instrumental eras. However, data are often stored as inconsistently formatted reports describing observed shaking and damage, making manually parsing and interpreting accounts labor‐intensive. We introduce a novel workflow using Google’s Gemini 2.5 Pro large language model (LLM) to automate the extraction and structuring of macroseismic observations from summary reports. We apply this workflow to the 22 March 1957 M 5.3 Daly City, California, earthquake as a case study. We used Gemini to extract addresses, originally assigned modified Mercalli intensity values, and descriptions from each report. To address coordinate precision limits, addresses were geocoded via Google’s Geocoding application programming interface. This workflow yielded over 2300 geocoded intensity reports for the Daly City earthquake. We use the geocoded accounts, with the original report intensity assignments, to develop a shaking intensity map that in some respects rivals modern Did You Feel It? Maps. We also extract and present data for the 9 February 1971 M L 6.7 Sylmar, California, earthquake. Our results demonstrate the potential of LLMs for reliably extracting and analyzing large, unstructured macroseismic datasets. LLMs offer a scalable solution for rapidly digitizing macroseismic archives, enabling their broader use to constrain ground‐motion models in modern seismic hazard analysis and to improve our understanding of site effects in urban areas. The concepts explored here may also be applied to the handling of other legacy seismological and earth science data.

Seismological Research Letters↗

USGS2024

Separating signals in elevation data improves supervised machine learning predictions for hydrothermal favorability

A recent study identified topography (land surface elevation above sea level) as an important input dataset (feature) for predicting the location of hydrothermal systems in the Great Basin in Nevada. Yet, topography is generally a result of more than one geological process and may consequently contain multiple distinct signals. For example, the geologic evolution of the Great Basin has produced both crustal thickening (i.e., regional-scale trends in elevation) and thinning via Basin and Range extensional faulting (i.e., valley-scale topographic relief). We postulate that these geologic processes may affect the occurrence of hydrothermal systems differently. Therefore, we separate the regional trend from the valley-scale signal in the Great Basin, and then use them separately to evaluate the importance of each as predictors for hydrothermal favorability. Our prior work applying supervised machine learning (ML) using the data from the Nevada Machine Learning Project demonstrated that employing a training strategy that randomly selects negative training sites produces better performing models for predicting hydrothermal favorability than a training strategy that uses expert-selected negatives. The models created using both training strategies exhibited a west-east geographic trend in the predictions for the favorability of hydrothermal resources. These models generally predicted higher favorability in western Nevada and lower favorability in eastern Nevada. This west-east trend in predicted favorability correlates with elevation across the Great Basin, which trends higher from west to east. By separating the original elevation feature into distinct features for elevation trend (i.e., regional-scale topography) and detrended elevation (i.e., valley-scale or local relative topography), we find that models using the separated topographic signals consistently outperform competing models that use the original elevation feature. Although western Nevada still exhibits higher favorability than eastern Nevada, using separated signals for regional elevation and local structure reduces the west-east prediction trend in the region and emphasizes structures associated with hydrothermal upflow. This work emphasizes how carefully engineering features to represent geological conditions relevant to hydrothermal systems allows ML algorithms to detect important patterns for predicting hydrothermal resource favorability and leads to better model performance.

Conference Paper↗

USGS2025

Drone-based radiometric surveys provide high-resolution mine waste characterization

Airborne radiometric surveys use passive geophysical techniques to characterize geochemical variations at or near earth’s surface. These methods have been used for a variety of mapping applications, including mineral resource evaluation. However, detailed characterization of smaller geologic targets, including mine waste features, requires flying at lower altitudes and with tighter line spacing than is feasible with traditional aircraft. Here, a small uncrewed aircraft system (sUAS) equipped with a radiometric sensor was used to acquire high-resolution gamma-spectrometry over small mine waste features and a low-grade stockpile in southwestern New Mexico. The sUAS radiometric system mapped local variability within each survey area and revealed ~2–10 m wide zones where radioelements K, Th, and U may be elevated 2–10× the surrounding material. Additionally, the sUAS radiometric data revealed radioelement variability across survey sites, which correlated reasonably well with variability seen in geochemical samples at each survey site, even though samples collected from individual sites showed high local variability. The sUAS data characterized local heterogeneity within mine waste and other small geologic targets at scales of a few meters to tens of meters, which is not possible with traditional crewed aircraft, and with continuity of coverage that is not possible with ground surveys, thus filling a key gap in geophysical survey spatial resolution.

New Mexico↗

USGS2024

Declining groundwater storage expected to amplify mountain streamflow reductions in a warmer world

Groundwater interactions with mountain streams are often simplified in model projections, potentially leading to inaccurate estimates of streamflow response to climate change. Here, using a high-resolution, integrated hydrological model extending 400 m into the subsurface, we find groundwater an important and stable source of historical streamflow in a mountainous watershed of the Colorado River. In a warmer climate, increased forest water use is predicted to reduce groundwater recharge resulting in groundwater storage loss. Losses are expected to be most severe during dry years and cannot recover to historical levels even during simulated wet periods. Groundwater depletion substantially reduces annual streamflow with intermittent conditions predicted when precipitation is low. Expanding results across the region suggests groundwater declines will be highest in the Colorado Headwater and Gunnison basins. Our research highlights the tight coupling of vegetation and groundwater dynamics and that excluding explicit groundwater response to warming may underestimate future reductions in mountain streamflow.

Nature Water↗

USGS2026

Geochemical disequilibrium at the brittle-ductile transition

We investigate the microtextural, microchemical, and isotopic effects of late-stage ductile deformation in quartzite mylonites and kyanite–muscovite–quartz veins from the Raft River shear zone (Utah). Quartz microstructures record pervasive disequilibrium, expressed by unannealed features including undulatory extinction, deformation lamellae, and poorly defined fabrics, typical of waning deformation in shear zones. Microchemical disequilibrium is best preserved in kyanite quartzites, where CL-zoned kyanite records repeated fractures, overgrowth, and mineral precipitation, and in muscovite from muscovite-poor quartzite mylonites that shows minor-element zoning consistent with syn-deformational overgrowth on detrital cores. In contrast, muscovite from muscovite-rich kyanite quartzites exhibits minimal chemical zoning. These microchemical variations correlate with 40 Ar/ 39 Ar age systematics. Chemically zoned muscovite preserves variable single-step ages, including ∼150 Ma ages, reflecting retention of detrital cores. In contrast, syndeformational muscovite consistently yields Miocene ages, indicating recrystallization and new growth below argon closure temperatures that reset inherited isotopic signatures. Similar trends are observed in quartzite mylonites, where increasing quartz recrystallization and stronger crystallographic preferred orientations occur toward deeper structural levels. Together, these observations indicate increasing retrograde deformation and recrystallization with depth in the Raft River shear zone and demonstrate that strain-driven recrystallization exerts a first-order control on muscovite ages. We suggest that apparent thermochronologic gradients in retrograde shear zones may reflect recrystallization gradients rather than temperature gradients. Where cooling limits the thermal driving force for recrystallization, isotopic relics are preserved, and local deformation and fluid availability control re-equilibration. Consequently, isotopic disequilibrium—particularly in the 40 Ar/ 39 Ar system—may be the rule rather than the exception of retrograde tectonic environments.

Idaho, Utah↗