USGS Science⌕ Search

SEARCH · USGS Science

Results for “Reviews of Geophysics”

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

To heal or not to heal?: 2. The moment-recurrence time behavior of repeating earthquakes in the 2011 Prague, Oklahoma aftershock sequence is consistent with laboratory healing rates

The timing and failure conditions of an earthquake are governed by the interplay between fault reloading and restrengthening. The moment-recurrence time behavior of repeating earthquakes can give observational estimates of fault healing rates; however, it is difficult to link these observed healing rates to laboratory studies of frictional healing in part because of uncertainty in lithology. Here, we study the 2011 Prague earthquake sequence, which includes repeating earthquakes in the Arbuckle group and the granitic basement, and compare them to laboratory experiments on samples of the Arbuckle and Troy granite (representative of the basement rock) (Okamoto et al., 2025, https://doi.org/10.1029/2024JB030573 ). We find three spatially distinct groups of repeating earthquakes with different moment-recurrence behavior: (a) constant moment-recurrence time in the Arbuckle group, (b) scattered moment-recurrence time at the intersection of the foreshock-mainshock fault in the granitic basement, and (c) moment-predictable behavior outside of the foreshock-mainshock fault intersection also in the granitic basement. Our observation of stagnant healing for repeating sequences in the Arbuckle group is consistent with laboratory observations of low healing rates for moderately high pore fluid pressures in Arbuckle samples. For the moment-predictable group, the source radius that is required in order to match healing rates is consistent with source radius estimations when taking into account reasonable attenuation of the P- pulse width. Overall, we observe diverse healing behaviors in the seismic families that are consistent with laboratory healing rates, providing seismic evidence that contact-scale frictional mechanisms are relevant to large-scale earthquake dynamics.

Oklahoma↗
USGS2026

Legacy of the fumigant 1,2-dibromo-3-chloropropane (DBCP) in California groundwater

The fumigant pesticide 1,2-dibromo-3-chloropropane (DBCP) was widely used in California agriculture during the 1960s and 1970s before being banned in 1979. Despite this ban, DBCP continues to contaminate groundwater due to its persistence and mobility. This study evaluates the distribution, historical trends, and projected persistence of DBCP in California using data from over 13,000 public supply wells and additional domestic, irrigation, and observation wells (1980-2022). Since 2010, DBCP has been detected in 9% of public supply wells statewide, with higher frequencies in the San Joaquin Valley (21%) and upper Santa Ana River watershed (13%), where DBCP use was most prevalent. Approximately 70% of wells had decreasing concentration trends, whereas increases were more common in deeper wells, indicating downward vertical migration of the DBCP front. Groundwater age estimates show that recharge timing aligns with the 1960s–1970s loading period, enabling reconstruction of peak inputs and providing a basis for age based modeling. To estimate future persistence, we applied a one dimensional advection–dispersion model that simulates long term declines in peak concentrations based on groundwater age, historical loading, and a 38 year degradation half life. Model projections suggest that concentrations above the maximum contaminant level may persist in a declining number of wells until approximately 2080 (range: 2048–2109), with longer persistence in the San Joaquin Valley. The simplified modeling framework, based on age distributions typical of wells capturing peak concentrations, can provide practical regional scale assessment of non-point source contaminants where long-term monitoring exists. This study highlights how the legacy of DBCP contamination will likely affect California's groundwater resources throughout the 21st century.

California↗
USGS2025

Cenozoic stratigraphy of Colorado

As a successor to previous Colorado stratigraphy charts (MS-53 Colorado Stratigraphic Chart and MS-54 Cretaceous Stratigraphy of Colorado), this Colorado Geological Survey (CGS) publication resulted from a collaboration between the CGS, USGS, and the Denver Museum of Nature and Science (DMNS). The chart was designed to illustrate Cenozoic stratigraphy spanning the state’s many sedimentary basins. It builds upon the work of dozens of colleagues and updates Richard Pearl’s seminal 1974 stratigraphy chart. The chart leverages the community’s stratigraphic work in both the subsurface and outcrop, and depicts new geochronologic constraints for many units. To facilitate comparison of strata to external forcing factors, the chart employs a linear timescale. Each unit’s dominant depositional environment is depicted, as are major mountain building events, erosional events, and regional unconformities.

Colorado↗
USGS2026

Surface rupture and slip distribution of the 2025 Mw7.7 Mandalay earthquake and updated length scaling of supershear earthquakes

The 2025 M w 7.7 Mandalay, Burma (Myanmar), earthquake ruptured 475 km of the central Sagaing fault and is the longest continental strike-slip rupture on record. The observed rupture length is 1.6–4.7 times the value expected (100–300 km) from existing length-magnitude scaling relations for strike-slip earthquakes. The earthquake resulted from shallow dextral faulting and ruptured bilaterally with supershear speeds south of the epicenter, rupturing close to three major cities in Myanmar and exposing over six million people to violent or extreme shaking. We report on the surface rupture character, length, and slip distribution based on sub-pixel correlation of Sentinel-2 (10 m) and Planet Dove (3 m) optical images and visual analysis of SkySat and WorldView (0.3–0.5 m) optical images. The earthquake had moderate surface slip (average = 3.3 m, maximum = 5.6 m, 25–75% range = 3.0–4.0 m), narrow deformation zone width (1–10 pixels in sub-pixel correlation and up to 190 meters for the detailed surface rupture mapping), and simple fault geometry (no stepovers or large changes in strike, 87% of the rupture that was mapped in detail is single-stranded). We attribute the extreme length of the Mandalay earthquake to supershear rupture speed, simple fault geometry, narrow down-dip width, and moderate surface slip. Based on a compilation of 25 supershear strike-slip earthquakes ( M w 6.5–8.6; 1979–2025), we find that the rupture length of supershear earthquakes does not fit empirical scaling relationships for strike-slip earthquakes that predict length from magnitude. A length-magnitude scaling relationship based on supershear earthquakes has a best fit of log 10 (surface rupture length) = 0.89 M w – 4.44, indicating that supershear earthquakes tend to be longer than their subshear counterparts for any given magnitude and thus may expose a greater population to shaking.

Mandalay↗
USGS2024

New U-Pb geochronology and geochemistry of Paleozoic metaigneous rocks from western Yukon and eastern Alaska, cross-border synthesis, and implications for tectonic models

The tectonic evolution of and relation between the Yukon-Tanana terrane and the Lake George assemblage, as well as other associated tectonic assemblages in western Yukon and eastern Alaska, have been debated for decades. The Yukon-Tanana terrane is widely considered to be an allochthonous rifted fragment derived from the Laurentian continental margin, whereas the Lake George assemblage and associated assemblages are currently interpreted to be part of the parautochthonous continental margin of western North America (Laurentia). To address these topics, we present 40 new U-Pb zircon ages and 20 new whole-rock geochemical analyses. We incorporate these data into a new compilation of available geological mapping for a large area that straddles the Alaska-Yukon border, together with 34 previously published U-Pb age determinations and an extensive geochemical database of metaigneous rocks from Late Devonian to Early Mississippian and middle to late Permian assemblages in this area. Magmatism in the Lake George assemblage and related assemblages occurred in two pulses from about 371 to 360 and from about 358 to 347 million years ago (Ma); geochemical discrimination diagrams indicate a large crustal component, possibly indicative of arc magmatism, for felsic metaigneous rocks and a range of tectonic environments for mafic rocks. Magmatism in the Fortymile River and related assemblages, and parts of the Nasina assemblage—all parts of the Yukon-Tanana terrane—are mainly Early Mississippian and span a crystallization age range from about 361 to 343 Ma; geochemical discrimination diagrams for these rocks indicate primarily arc geochemical signatures for both mafic and felsic rocks. Middle to late Permian crystallization ages (about 261–253 Ma) are indicated for felsic metaigneous rocks in the Klondike assemblage and some of the felsic metaigneous rocks in the Nasina assemblage. Based on our mapping, we propose the existence of a possible unconformity between the Mississippian and Permian felsic metavolcanic rocks within the Nasina assemblage that is marked by sporadic occurrences of stretched-pebble conglomerate. Our combined database supports the well-established model of a magmatic arc comprising the Fortymile River and Finlayson assemblages of the rifted Yukon-Tanana terrane continental fragment on which a middle to late Permian arc (Klondike assemblage) was later built. The assemblages of the Yukon-Tanana terrane were subsequently intruded by Late Triassic to Early Jurassic granitoids, presumably during reaccretion of the Yukon-Tanana terrane to the continental margin. Permian and Late Triassic to Early Jurassic intrusions have not been mapped in the now structurally lower plate Lake George assemblage; their absence is one of the lines of evidence that have been used to support the parautochthonous, rather than allochthonous, origin of the Lake George assemblage and related assemblages. Our new data, together with previously published ranges of igneous crystallization ages and geochemical tectonic signatures of the Late Devonian to Early Mississippian magmatic rocks in the Lake George assemblage and associated assemblages and in the Fortymile River, Nasina, and correlated assemblages of the Yukon-Tanana terrane, indicate that the currently accepted interpretation of the Lake George assemblage and associated rocks being part of parauthochthonous North America is not the only possible interpretation of this tectonic entity. Approximately half of the dated intrusive rocks in the Lake George assemblage are contemporaneous with the metaigneous rocks of the Yukon-Tanana terrane arc (<361 Ma). We speculate that our approximately 361 Ma U-Pb age for quartz syenite in part of the North American continental margin in south-central Yukon defines the beginning of rifting of the Laurentian margin. Although the currently favored model of prolonged middle Paleozoic subduction and extension in both the Yukon-Tanana terrane and parautochthonous North America allows for simultaneous middle Paleozoic magmatism on both sides of the Slide Mountain Ocean, we now propose an alternative hypothesis in which the Lake George assemblage represents a deeper part of the rifted Yukon-Tanana terrane arc. If this is the case, the absence of Permian and Late Triassic to Early Jurassic arc rocks in the Lake George assemblage could be explained either by the arcs of these ages not being wide enough to have affected the Lake George assemblage or by tectonic displacement of these arc rocks away from the Lake George assemblage. Our approximately 259 Ma U-Pb zircon age and geochemical analyses of metarhyolite in the Seventymile terrane in Alaska, which comprises remnants of the back-arc basin that separated the Yukon-Tanana terrane from the Laurentian continental margin, confirm the presence of a late middle Permian volcanic arc component to the terrane. Our approximately 319 Ma U-Pb zircon age from the Chicken assemblage (as redefined in this study) in eastern Alaska, combined with previously reported fossil ages and a U-Pb zircon age from this assemblage, indicate that it is a Late Mississippian to Early Pennsylvanian arc assemblage. We propose several other relatively young, locally developed arc assemblages outboard of the ancient continental margin of Laurentia that may correlate with the Chicken assemblage, but we consider its origin to remain an enigma.

Alaska↗
USGS2024

Assessment of long-term changes in surface-water extent within Klamath Marsh, south-central Oregon, 1985–2021

The annual maximum extent of surface water in Klamath Marsh has naturally fluctuated in response to periods of wet and dry conditions in the surrounding basin. Field observations during the 2010s indicate that the annual maximum extent of surface water has been declining and the marsh is not responding to hydrologic inputs as it had historically. This report describes the results of a hydrologic evaluation of Klamath Marsh to characterize and understand multi-year declines in the surface-water extent and increased intermittency of streamflow exiting the marsh. Landsat imagery collected during 1985–2021 was processed to create a time series of annual maximum surface-water extent to assess changes in surface-water inundation within the marsh. A 50-percent decrease in the mean surface area of annual total open-water extent (TOWE) during the latter half of the study period (2003–21) compared to the first half (1985–2003) was observed in this 37-year time-series dataset. The change in open-water extent was offset by a corresponding increase in dry land in the marsh. Time series of streamflow, groundwater level, total annual precipitation, annual mean temperature, and anthropogenic water use and water management were compiled and evaluated to improve understanding of the factors affecting TOWE. Statistically significant downward trends in the regional groundwater table and streamflow into and out of the marsh were identified as well as statistically significant upward trends in annual mean temperature. Statistically significant correlations among TOWE, streamflow, and groundwater level also were identified. The decreasing trends could not be attributed to changes in total annual precipitation or changing anthropogenic groundwater use within the study area. Declines in the open-water extent of Klamath Marsh since 2000 principally are due to a decoupling of the groundwater and surface-water system beneath the marsh because of regional declines in groundwater level. Regional increases in air temperature and the reestablishment of more than 55,000 acres of forested land within the study area have likely contributed to increasing evapotranspiration, leaving less water available for groundwater recharge and stream base flow and resulting in basin-wide declines in streamflow and groundwater levels.

Oregon↗
USGS2024

High resolution identification and quantification of diffuse deep groundwater discharge in mountain rivers using continuous boat-mounted helium measurements

Discharge of deeply sourced groundwater to streams is difficult to locate and quantify, particularly where both discrete and diffuse discharge points exist, but diffuse discharge is one of the primary controls on solute budgets in mountainous watersheds. The noble gas helium is a unique identifier of deep groundwater discharge because groundwater with long residence times is commonly enriched in helium. In this study, a portable mass spectrometer was used to measure longitudinal variation in dissolved helium concentrations in two mountainous rivers at high spatial resolution not feasible with traditional sampling techniques. Helium profiles were then simulated using a mass-balance model to quantify longitudinal variation in groundwater discharge to the receiving rivers. Results indicate helium concentrations were enriched by multiple orders of magnitude above atmospheric equilibrium in both rivers and that this persisted for up to 18 km below observed pulse inputs in the Colorado River. Helium mass-balance models match observed longitudinal patterns with the exception of sharp initial increases in helium observed in the rivers. Increased longitudinal groundwater discharge rates correspond to mapped geologic structures in both watersheds that likely transport deep geothermal water. Models show variable sensitivity to spatial assignment of input variables representing the groundwater source, illustrating the importance of collecting data from discrete groundwater discharges where possible. The methodology shows promise for field experiments designed to assess air–water exchange rates and to quantify total groundwater discharge from a combination of discrete and diffuse sources.

Colorado, Utah↗
USGS2024

Sensitivity of Alabama freshwater gastropod species to nickel exposure

Snails are effective bioindicators due to their prolific distribution, high level of endemism, and capacity to accumulate contaminants. Freshwater snails have unique ecological niches which are imperiled by land‐use change and the introduction of hazardous chemicals. To assess how environmental alterations affect gastropods, lab‐based studies are needed to characterize the toxicity of specific stressors. This can help guide policy decisions and remediation efforts. The aim of this research was to characterize acute toxicity of nickel (Ni) on endemic snails ( Somatogyrus georgianus [Walker, 1904], Elimia cahawbensis [Lea, 1861], and Elimia spp.) and measure the accumulation of Ni and mineral elements including calcium (Ca), magnesium, potassium, and sodium (Na). Snails were exposed to six concentrations (25–800 µg/L) of Ni for 96 h. Among the studied snail species, E. cahawbensis was the most sensitive to Ni, with the lowest lethal concentration where 50% of the organisms died (LC50) at 88.88 µg/L Ni after 96 h. The LC50 at 96 h for S. georgianus was 167.78 µg/L Ni, and 393.13 μg/L Ni for Elimia spp. Except for Elimia spp., mortality of the other two snail species corresponded to the whole‐body uptake of Ni. Nickel exposure also influenced Ca and Na uptake for Elimia spp. All three endemic species are potential candidate species for evaluating localized effects of human activities, and the present study provides a first step in characterizing how snails would be affected by environmental alterations. More research could further characterize potential effects of other human stressors on these endemic snail species. Future research into subindividual responses and routes of exposure can further elucidate variations in species sensitivity.

Alabama↗
USGS2025

Nanometer-scale relationships between sedimentary organic matter molecular composition, fluorescence, cathodoluminescence, and reflectance: The importance of oxygen content at low thermal maturities

Molecular characterization of sedimentary organic matter (SOM), termed macerals, is a common goal when seeking to understand petroleum generation as well as other geologic processes in deep time. However, unambiguous measurement of discrete macerals is challenging due to the small size of organic particles in sedimentary rocks, the proximity of different organic matter types to one another, mineral-organic matter interactions, and maceral mixing that occurs during SOM isolation prior to ex situ analysis. The recent advent of infrared spectrometers capable of nanometer-scale resolution and the application of these technologies to geologic samples has enabled advances in rapid, in situ molecular characterization of SOM allowing for insights into paleoenvironmental processes, such as organic matter productivity and preservation, among others. Here we employ one such technology, optical photothermal infrared (OPTIR) spectroscopy, to map SOM functional group distributions at 500-nm resolution in a sample from the Lower Cretaceous Sunniland Limestone of the South Florida Basin. Examined fields of view include occurrences of amorphous organic matter (AOM), inertinite, micrinite, solid bitumen, telalginite, and vitrinite. OPTIR data from these macerals are compared against traditional organic petrographic data from the same organic grains including fluorescence intensity and white light reflectance as well as against cathodoluminescence response, an emerging organic petrographic approach. Maceral oxygen content (using carbonyl functional group abundance as a proxy) is observed to vary widely between maceral types but correlates strongly with fluorescence and cathodoluminescence intensity as well as against reflectance. These findings highlight the important role that oxygen content plays in determining the optical properties of SOM and further demonstrate the ability of OPTIR to discriminate subtle molecular differences between SOM types.

Organic Geochemistry↗
USGS2025

High frequency and region-scale simulations of large (Mw7+) earthquakes on the southern Whidbey Island fault, Washington, USA

We simulate ground shaking in western Washington State from hypothetical M w 7.0–7.5 earthquakes on the southern Whidbey Island fault (SWIF). Ground motions are modeled considering kinematic source distributions on a complex fault plane, a 3D seismic velocity model, and region‐specific soil velocity models. We run simulations with varying model resolutions, including regional‐scale simulations with a maximum‐modeled frequency of ∼1 Hz and local‐scale simulations with a maximum‐modeled frequency of ∼2.5 Hz. Additional local‐scale simulations are run considering high‐resolution surface topography. We explore how source parameters (i.e., magnitude, hypocenter location, and dip direction) and 3D velocity structure impact peak shaking intensity and its variability. In particular, we find that earthquakes on the SWIF would likely produce strong shaking throughout the populated Puget Lowland, including in the cities of Everett, Seattle, Bellevue, and Tacoma, Washington. Simulated short‐period ( T ≤ 2 s) spectral accelerations are strong throughout the Puget Lowland, and long‐period shaking ( T ≥ 5 s) is strong in the deep regional sedimentary basins, especially the Everett and Seattle basins. Source parameters strongly influence intra‐ and interevent variability in response, primarily through changes in source and site geometry, as well as rupture directivity. We also note a potential coupling between rupture directivity and basin effects, wherein directivity pulses are seemingly guided through the region’s deep, interconnected sedimentary basins. Overall, this work highlights the impacts of 3D source, path, and site effects on seismic hazard in the U.S. Pacific Northwest and substantially expands the catalog of simulated ground motions for Puget Sound area crustal faults.

Washington↗
USGS2024

Hypothetical CO2 leakage into, and hydrological plume management within, an underground source of drinking water at a proposed CO2 storage facility, Kemper County, Mississippi, USA

A large Geologic Carbon Sequestration (GCS) hub has been proposed in Kemper County, Mississippi. The target injection interval consists of numerous Cretaceous-aged deep saline aquifers overlain by a competent and extensive regional sealing layer. Above the seal, the deepest Underground Source of Drinking Water (USDW) at the site is the Eutaw aquifer of the Eutaw Group and McShan Formation, undifferentiated. To assess potential risks of leakage from the deep sequestration reservoir, a model of a portion of the Cretaceous Eutaw Group was constructed in this study. Simulations tested various permeabilities, hypothetical leakage rates, and plume mitigation strategies utilizing existing wells. Results suggest that, under the influence of regional groundwater flow fields, leaking CO 2 would effectively bypass the existing wells, and to influence this migration would require very large water extraction rates. Therefore, to ensure plume detection, monitoring for leakage at the injection wells themselves is very important.

Mississippi↗
USGS2026

The use of fluorite geochemistry and machine learning to identify critical mineral systems

Fluorite (CaF 2 ) is a potential pathfinder to critical mineral and rare earth element (REE) deposits but its application has been limited to a narrow range of mineralization types. I show that fluorite is a robust recorder of mineralization fertility by applying statistical and machine-learning methods to a new global fluorite geochemical database. Distinct median rare earth and trace element patterns are observed among deposit types and genetic environments. Fluorite associated with carbonatites and REE deposits are relatively enriched in Sr and have minimal Eu anomalies. These characteristics define new bivariate discrimination diagrams that correctly identify 78% of carbonatite-related fluorite and 88% of fluorite from REE deposits. Random forest classifiers were developed for a wide range of mineralization types and genetic settings. Trained solely on rare earth element patterns, these models achieve accuracies of 77–79%. Higher classification accuracies (up to 88–96%) are obtained when including elements such as Sr, highlighting the significance of trace elements for optimal fluorite classification. The recognition of diagnostic fluorite compositional fingerprints, particularly in REE-fertile systems, underscores its potential as a pathfinder and indicator for critical mineral exploration in F-bearing environments.

Mineralium Deposita↗
USGS2024

How does the onset of offset influence geologic slip rates?

Geologic slip rates are typically based on the displacement accrued by a geomorphic or stratigraphic feature and the age of the offset feature. Because slip rates are commonly calculated by dividing the displacement of a faulted marker by its age, they contain two open time intervals: the elapsed time between the age of an offset feature and the age of the earthquake that displaced the feature, and the time between the present‐day and the most recent earthquake. Here, we explore the influence of including unconstrained open intervals in geologic slip rate calculations. We test the degree to which these open intervals affect geologic slip rates and their uncertainties, and we find that their influence depends primarily on mean earthquake recurrence intervals (RIs). Slip rates on faults with longer RIs, such as the Wasatch fault, can be greatly influenced by an increase of up to 20% when accounting for open intervals. In contrast, slip rates on faults with shorter RIs, such as the San Andreas fault, are only slightly influenced by the assumption that slip rates calculated over open intervals approximate those calculated over closed intervals. Our analyses indicate that faults with moderate slip rates (∼0.2–5 mm/yr) are sensitive to both open interval effects themselves, as well as methods to quantify and account for these effects. We re‐evaluate how slip rates are calculated and defined in displacement–time space using published deformation records. We explore the utility of assigning a probability distribution to the initiation of offset of the oldest faulted feature and the timing of the most recent earthquake (MRE). We find that calculating geologic slip rates without using probability distributions that capture the timing of the MRE and the onset of offset of the oldest faulted feature, especially on slow‐to‐moderate slip rate faults, can lead to systematic underestimation of average geologic slip rates.

Seismological Research Letters↗
USGS2025

Integrated stratigraphic and geochemical analysis of organic-rich intervals of the Lewis Shale in the eastern Washakie Basin, Wyoming

Geological studies in the Cretaceous Western Interior Seaway (KWIS) in North America reveal highly variable sedimentological conditions on its western shore caused by rapidly changing sea level and detrital input during the seaway’s closure. Here we examine a 601-foot (183 meters) continuous core through the lower part of the Maastrichtian Lewis Shale in the eastern Washakie Basin, Wyoming, through integrating stratigraphic and geochemical analyses to better understand organic matter deposition and preservation during the final marine transgression within the seaway. The core penetrates eight organic-rich flooding intervals (F1-F8) and a regional condensed section, informally known as the Asquith marker. The lower portions of the core record sediment input from a southern source, likely the Sierra Madre/Park Ranges, while the upper part records sediment input from a northern source, likely the Granite Mountains. This provenance transition is supported by regional stratigraphic analysis and changes in bulk mineral and trace metal composition. The Asquith marker and early flooding surfaces are enriched in oil-prone, marine organic matter deposited under dysoxic to euxinic conditions, whereas younger flooding surfaces show increased terrigenous input and poorer preservation conditions. The Asquith marker is identified as a prime oil-prone source rock. Although younger flooding surfaces also exhibit favorable source-rock properties, their generative potential is reduced due to increased clastic and terrigenous organic matter deposition caused by regression. These results emphasize the importance of sediment source variability, organic matter preservation, and changing redox conditions to provide insight into the sediment provenance and petroleum potential of the Lewis Shale in the context of the final transgression within the KWIS.

Colorado, Wyoming↗
USGS2026

A process-based model for forecasting wave runup along the coast of Georgia

Wave runup is an important nearshore process that impacts total water level, sediment transport, and coastal design. Current methods for forecasting wave runup implement an empirical model that considers offshore wave height, wave period, and generalized beach slope. In this study, the authors generated wave runup forecasts from offshore wave conditions and a system of polynomial equations derived from numerical simulations at three different still water datums for each beach profile. They developed a process-based methodology that incorporated site-specific cross-shore topobathy into the phase-resolving numerical model. A comparison between the system of equations, deterministic hydrodynamic simulations, and observed high-water marks was made using Hurricanes Matthew (2016) and Irma (2017) for 12 cases, and it showed that the polynomials were capable of being consistent with the results from full simulation runs, while not requiring hours of runtime when a forecast was needed—the differences between the polynomial and the observed high water marks ranged from 3 to 32 cm for the Irma hindcast and 9–70 cm for Matthew. Then, using forcings from Hurricanes Ian and Nicole (2022), the model predicted the occurrence of dune collision, overwash, and inundation for the coast of Georgia and suggested that wave runup was impacted by the still water level and local topobathy.

Georgia↗
USGS2025

Depth of magma crystallization and fluid exsolution beneath the porphyry-skarn Cu deposits at Santa Rita and Hanover-Fierro, New Mexico, USA

The depth level at which porphyry Cu–forming magmas fractionated and exsolved mineralizing fluids is actively debated. In the classic model, extensive magma fractionation occurs in large, upper crustal magma chambers, and concomitant fluid exsolution leads to forceful expulsion of residual magmas in the form of porphyry dikes, stocks, and breccia pipes, which subsequently serve as pathways for the mineralizing fluids. In contrast, some recent studies highlighting the role of deep crustal magma fractionation in the production of fertile magmas essentially deny the existence of upper crustal magma chambers at the time of mineralization. To address this, we conducted a detailed thermobarometric investigation of 13 intermediate to felsic, porphyritic intrusive rocks related to porphyry-skarn Cu mineralization at Santa Rita and Hanover-Fierro, New Mexico, United States, representing two premineralization magmas (61–60 Ma), seven synmineralization magmas (60–58 Ma), and four late- to postmineralization magmas (58–57 Ma). For each sample, the pressure of last magma crystallization before final magma ascent to the current exposure level was reconstructed based on Al-in-hornblende barometry of small hornblende inclusions trapped within quartz phenocrysts and through titanium-in-quartz (TitaniQ) thermobarometry of the host quartz phenocrysts themselves. Since quartz is one of the last crystallizing magmatic minerals, and no significant phenocryst growth could have occurred in small dikes and stocks after final magma emplacement, quartz phenocrysts and their contained hornblende inclusions record the depth of last magma crystallization before final magma ascent. When present, hornblende phenocrysts and hornblende inclusions within other major phenocrysts were also analyzed. Both quartz and hornblende barometers return consistent average pressures of 3.2 ± 0.4 kbar for the entire suite of pre- to postmineralization magmas, corresponding to depths of 11 to 14 km. The synmineralization magmas return even more consistent average pressures of 3.1 ± 0.2 kbar, corresponding to a depth of 12 ± 1 km. The volume of the mineralizing porphyry dikes and stocks at the emplacement level is far too small to have provided all the fluids and metals required to form the observed ore deposits. Therefore, the majority of the ore-forming fluids must have originated from the magmas that crystallized at 12 ± 1 km depth. The ore deposits, conversely, formed at ~5-km paleodepth. This implies that most of the mineralizing fluids traveled an average vertical distance of ~7 km from their magmatic source to the eventual site of ore precipitation. The relatively unaltered nature and low veining degree of deeper parts of mineralized porphyry dikes and stocks suggest that the fluid transport through these intrusive bodies occurred mostly at near-solidus conditions by means of fluid percolation along grain boundaries. In summary, our results suggest that (1) a large, upper crustal pluton exists ~7 km beneath the Santa Rita and Hanover-Fierro deposits; (2) abundant phenocryst crystallization occurred at this depth level; and (3) this pluton was the main source for the exsolution of ore-forming fluids. However, the investigated rocks have elevated whole-rock Sr/Y ratios, indicating magma fractionation at deep crustal levels. As a result, our preferred model is a combination of the two end-member models introduced above, with most magma fractionation having occurred in the deep crust and with residual, intermediate to felsic melts having ascended and accumulated at 11 to 14 km paleodepth, where they continued to crystallize with comparatively little crystal-liquid separation, before some of these magmas ascended further to shallow levels and quenched to porphyries.

New Mexico↗
USGS2024

Characterization of change in tree cover state and condition over the conterminous United States

Variability in the effects of disturbances and extreme climate events can lead to changes in tree cover over time, including partial or complete loss, with diverse ecological consequences. It is therefore critical to identify in space and time the change processes that lead to tree cover change. Studies of change are often hampered by the lack of data capable of consistently detecting different types of change. Using the Landsat satellite record to create a long time-series of land cover and land cover change, the U.S. Geological Survey Land Change Monitoring Assessment and Projection (LCMAP) project has made an annual time series of land cover across the conterminous United States for the period 1985 to 2018. Multiple LCMAP products analyzed together with map validation reference plots provide a robust basis for understanding tree cover change. In LCMAP (Collection 1.2), annual change detection is based on harmonic model breaks calculated at each Landsat pixel from the Continuous Change Detection and Classification (CCDC) algorithm. The results showed that the majority of CCDC harmonic model breaks (signifying change) indicated partial tree cover loss (associated with management practices such as tree cover thinning) as compared to complete tree cover loss (associated with practices like clearcut harvest or fire disturbance). Substantially fewer occurrences of complete tree cover loss were associated with change in land cover state. The area of annual tree cover change increased after the late 1990s and stayed high for the rest of the study period. The reference data showed that tree harvest dominated across the conterminous United States. The majority of tree cover change occurred in evergreen forests. Large estimates of disturbance-related tree cover change indicated that tree cover loss may have previously been underreported due to omission of partial tree cover loss in prior studies. This has considerable implications for forest carbon accounting along with tracking ecosystem goods and services.

Forests↗