SEARCH · USGS Science

Results for “River Research and Applications”

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

USGS2025

Constraining landslide frequency across the United States to inform county-level risk reduction

Informative landslide hazard estimates are needed to support landslide mitigation strategies to reduce landslide risk across the United States. Whereas existing national-scale landslide susceptibility products assess where landslides are likely to occur, they do not address how often , which is a critical element of landslide hazard and risk assessments. In particular, the U.S. Federal Emergency Management Agency's National Risk Index (NRI) requires landslide frequency estimates to inform expected annual loss estimates. We present county-level landslide frequency (landslides per area per year) estimates for the 50 US states. We applied Bayesian negative binomial regression to estimate both the expected (average) reported landslide frequency and full distribution of annual landslide counts for each county. We compared a suite of models that used combinations of landslide-susceptible area, probability of potentially triggering earthquakes, frequency of potentially triggering precipitation, and ecological region as predictors. We trained our models with landslide inventory data from counties with the most comprehensive records available nationwide and used zero-inflated negative binomial distributions as an incompleteness model to correct for temporal reporting gaps. We selected a preferred frequency model to inform the NRI based on information criteria and physically plausible parameter estimates. The model showed that average annual reported landslide frequencies vary by 5 orders of magnitude across US counties, ranging from 0.002 (0.00015–0.05) landslides 1000 km −2 yr −1 in Kusilvak Census Area, Alaska, to 29 (19–46) landslides 1000 km −2 yr −1 in Lake County, California, reflecting the country's strong variations in landslide susceptibility, earthquake probability, and other factors for which ecological region serves as a proxy. Counties with estimated frequencies in the top 20 % of all counties are predominately along the West Coast of the continental United States, in mountainous regions of the Pacific Northwest and Intermountain West, in locally steep or earthquake-prone regions of the Midwest and Southeast, along the Appalachians, in southern and southeastern Alaska, and on some Hawaiian islands. By examining the number of landslides predicted in 99th percentile years for each county, we identified that 26 % of US counties likely have potential for widespread landsliding with more than 10 landslides 1000 km −2 yr −1 , even when such large events have not been reported in the training data for that county. Overall, our results better represent the range of possible landslide frequencies and spatial variations than previous national-scale estimates reported in the NRI, and our approach can inform other risk-reduction and loss-mitigation efforts across the United States and globally.

Natural Hazards and Earth System Sciences↗

USGS2026

Rare earth element potential in coal and coal ash in the U.S. Gulf Coast

United States heavy reliance on imports of critical minerals (CMs), including rare earth elements (REEs), underscores the importance of development of domestic sources. The study objective was to quantify CM and REE concentrations in coal and coal ash in the US Gulf Coast region. CM and REE concentrations were measured for 118 samples from outcrops and 14 mines in the Gulf Coast. Results show that total REE + Yttrium (REY) concentrations (dry coal basis) are comparable to those of the upper continental crust (UCC) with localized hot spots, such as the Texas Gibbons Creek mine (REY ≤ ~ 2860 ppm). When normalized to UCC REY concentration (169 ppm, dry coal basis), REY to UCC ratios for Gulf Coast coal samples range from 0.1 to 17 (median ratio 0.6). REE extractability from lignites is high (median: 63%–93%) using environmentally benign weak acid. In addition to raw coal, coal ash from power plants could also serve as an REE source with a median ratio of REY in ash relative to coal of 4; however, extractability from coal ash is generally much lower (≤ 5% using the same weak acid as in coal). The median basket price for extracted REY as oxides from coal, assuming 70% extractability, is $3.2 per tonne of coal and $186 billion based on 58 billion metric tonnes of dry coal in the Gulf Coast. REEs important for magnets (Pr + Nd + Tb + Dy) account for ~ 80% of the total value. The corresponding median basket price for extracted REY as oxides from coal ash, assuming ~ 30% extractability, is ~$4.4 per tonne of ash and $1.2 billion based on 258 million tonnes of ash. REE production from coal would likely require co-products, such as activated carbon or humic acids, to attain economic viability. Production of REEs from coal ash could offset remediation costs related to potential water contamination. This reconnaissance study shows the potential for REE production from coal and coal ash in the Gulf Coast; however, carbon coproducts and/or societal benefits would likely be required for socioeconomic viability.

Alabama, Florida, Louisiana, Mississippi, Texas↗

USGS1991

Wetland sedimentation and vegetation patterns near selected highway crossings in West Tennessee

Wetland sedimentation and vegetation patterns at 11 highway crossings in West Tennessee were studied from 1987 to 1989. The purpose of the study was to investigate potential adverse effects of highway crossings on wetlands. Sedimentation rates, determined from root-burial depths, were highly variable. Average rates offine-grained deposition ranged from 0.005 to 0.033 foot per year for stations in locally ponded areas and from -0.002 to 0.039 foot per year for stations in drained areas. Sedimentation rates upstream from highway crossings were not significantly different from downstream rates at 8 of the 11 study sites. Three study sites had significantly greater sedimentation rates downstream. Sand splays were observed downstream from bridges at most study sites. Vegetation patterns and tree growth appear most strongly related to hydropetiod, defined as the average length of time an area is covered by water each year. The influence of sedimentation on tree growth is difficult to separate from the influence of hydroperiod because areas with high sedimentation rates typically have long hydroperiods. Estimated hydroperiod increased no more than 1 percent because of backwater from the highway crossings at the 11 study sites, while the estimated average depth of flood-plain inundation increased by an average of 6 percent.

Tennessee↗

USGS2025

Critical mineral inventory of select IOA-IOCG deposits, southwestern USA

Critical minerals are necessary for modern technology and strategic purposes. Their increasing importance requires finding new and nontraditional resources. Samples of ore, altered, and unaltered host rock were collected from 26 iron mines and prospects in California, Nevada, and Utah to assess the potential of these deposits to host economic quantities of different critical minerals. Geochemical analyses were conducted by 61 element ICP-OES-MS sodium peroxide fusion and major elements determined by WDXRF. These deposits concentrated many critical minerals beyond what is found in average upper crustal abundances, such as Sb, As, Bi, Co, Ga, Mg, Mn, Ni, Nb, Pd, REE, Sc, Te, Sn, Ti, W, V, and Zn. However, most of these are not concentrated enough in the ore to be considered as economic resources. Those critical minerals that are enriched enough in some of these deposits to possibly be considered as by-product commodities are Ni, REE, V, and potentially Co and Ga. These enrichments were not uniform, with REE more likely to be enriched in IOA deposits, whereas Co, Ga, Ni, and V could be found enriched in either IOA or IOCG deposits.

California, Nevada, Utah↗

USGS1994

Hydrogeology of the Chickasaw National Recreation Area, Murray County, Oklahoma

The Travertine District (Park) of the Chickasaw National Recreation Area, operated and maintained by the National Park Service, is near the City of Sulphur in south-central Oklahoma. The Park was established in 1902 because of its unique hydrologic setting, which includes Rock Creek, Travertine Creek, numerous mineralized and freshwater springs, and a dense cover of riparian vegetation. Since the turn of the century several flowing artesian wells have been drilled within and adjacent to the Park. Discharge from many of these springs and the numbers of flowing wells have declined substantially during the past 86 years. To determine the cause of these declines, a better understanding of the hydrologic system must be obtained. The U.S. Geological Survey, in cooperation with the National Park Service, has appraised hydrologic information obtained for the Park from several studies conducted during 1902-87. The principal geologic units referred to in this report are the Arbuckle Group and the overlying Simpson Group. These rocks are of Upper Cambrian to Middle Ordovician age and are composed of dolomitic limestone, with some sandstones and shales in the Simpson Group. Surface geologic maps give a general understanding of the regional subsurface geology, but information about the subsurface geology within the Park is poor. The Simpson and Arbuckle aquifers are the principal aquifers in the study area. The two aquifers are not differentiated readily in some parts of the study area because of the similarity of the Simpson and Arbuckle rocks; thus, both water-bearing units are referred to frequently as the Arbuckle-Simpson aquifer. The aquifers are confined under the Park, but are unconfined east and south of the Park. Precipitation on the outcrop area of the Arbuckle aquifer northeast and east of the Park recharges the freshwater springs (Antelope and Buffalo Springs) near the east boundary of the Park. The source of water from mineralized springs located in the central part of the Park, and flowing wells within and north of the Park, is believed to be a mix of waters from rocks of the Arbuckle and Simpson Groups. The source of water from two highly mineralized springs, Bromide and Medicine, that ceased to flow in the early 1970’s is believed to be from the Simpson Group. Water-quality characteristics reflect the sources of ground water in the study area. The highly mineralized springs near the western end of the Park are a sodium chloride type with dissolved solids greater than 4,500 mg/L. The freshwater springs near the eastern end of the Park are a calcium bicarbonate type with total dissolved solids of less than 400 mg/L. Flow from the artesian wells has declined substantially during the past 86 years and the wells are estimated to currently discharge only about 10 percent of the total flow reported in 1939. The depletion is believed to be caused by a gradual lowering of the hydraulic head within the aquifer. The influence on the hydrologic system of local municipal and industrial pumping from the Arbuckle-Simpson aquifer is difficult to discern because the system is much more sensitive to precipitation than to pumpage. Ground-water levels and spring flows in this region respond rapidly to precipitation. The effects of withdrawals from the City of Sulphur and Oklahoma Gas and Electric Company power-plant water-well fields are not discernible at wells and springs. The hydrologic system may be influenced by pumping, particularly during extended dry periods of several years, but the impact of pumping on the system cannot be determined without further investigation.

Oklahoma↗

USGS2025

Probing lower-crustal fault properties with frequency-dependent tidal tremor triggering

The way seismicity responds to periodic stress perturbations offers crucial insights into the processes that can trigger an earthquake. Laboratory and theoretical analyses have shown that the period of imposed forcing and source properties affect the sensitivity to triggering, but frequency-dependent triggering of tectonic faults is poorly understood. The rate of low-frequency earthquakes (LFEs) near Parkfield, California has been found to be strongly correlated with solid earth tides. Tidal forcing acts over multiple frequencies, and the sensitivity to tidal triggering of LFEs therefore provides a unique opportunity to probe the physics of earthquake triggering and underlying fault properties. Here, we find that the response of LFEs to solid earth tides at diurnal and semi-diurnal frequencies is highly variable but spatially coherent along the San Andreas Fault. Using rate-state friction modeling, we find that the variation of the amplitude of tidal modulation is mainly affected by the spatial variation of the background effective normal stress, whereas the spatially varied tidal modulation at diurnal and semi-diurnal frequencies is mainly affected by the fault frictional property. The spatial variations of the LFE response reveal a heterogeneous lower-crustal geologic structure and complex physical faulting processes below the rupture zone of eventual great earthquakes.

California↗

USGS2026

Simulation of the impacts of spring diversions on streamflow in the Strawberry Creek watershed, San Bernardino County, California, using an integrated hydrological model

The Strawberry Creek watershed, in the San Bernardino Mountains of southern California, features a group of natural springs known as Arrowhead Springs that have been augmented with diversions in the form of sub-horizontal borings and tunnels. Understanding the effects of these structures on streamflow through groundwater capture is crucial for managing surface-water resources in this watershed. In this study, we constructed the Strawberry Creek integrated hydrological model to increase this understanding. The Strawberry Creek integrated hydrological model is an integrated surface runoff and groundwater model, GSFLOW to simulate surface runoff and infiltration and groundwater flow. The model has 263 rows and 176 columns in each of its three layers. The model area includes the Strawberry Creek and four adjacent watersheds. The precipitation-runoff modeling system was calibrated using two streamflow gaging stations, and the groundwater and surface-water flow model was calibrated to the reported spring diversion discharge and a sparse number of groundwater-level measurements. The Strawberry Creek integrated hydrological model was run with and without diversions active, and simulated streamflow was compared, indicating that in the headwaters of Strawberry Creek, about 35 percent of the diversion flow was captured from base flow.

California↗

USGS2025

Spatio-temporal evolution of distributed volcanic fields, case studies—Sierra Chichinautzin and Michoacán-Guanajuato, México

An analysis of 1,375 volcanoes in the Michoacán-Guanajuato (1,148 volcanoes in a 26,200 square-kilometer area) and Sierra Chichinautzin (227 volcanoes in a 3,500 square-kilometer area) volcanic fields in central Mexico identified patterns in the spatial and temporal distribution of past eruptions. A cluster agglomerative hierarchical method and kernel analysis confirmed that the Michoacán-Guanajuato volcanic field comprises four volcanic fields (Valle de Santiago, Uruapan, Apatzingán, and Pátzcuaro volcanic fields) controlled by different fault systems, indicating that it is not a single volcanic field but rather a group of volcanic fields (a “superfield”), each of which has distinct characteristics. In the Sierra Chichinautzin volcanic field, well-constrained isotopic ages were used to build a model of how the spatial distribution of the eruptions has changed over time. Two new 40Ar/39Ar ages from a locally recognized volcanic feature near the town of El Cantil, herein called El Cantil volcano (1,537±17 kilo-annum [ka]) and the volcanic feature at Cerro el Elefante (herein called El Elefante dome) (1,485±92 ka) belong to the oldest volcanic group identified in the Sierra Chichinautzin volcanic field, confirming the timing of the beginning of monogenetic volcanism in the region. Based on the volcanic groups identified in the Sierra Chichinautzin volcanic field, the youngest volcanism (less than 35 ka) is found only in the central-western sector of the field. Principal component analysis determined the directional trends of feeder dikes only for vents <10 ka in the Sierra Chichinautzin volcanic field. Possible magma migration paths through the crust were identified using seismic data from both volcanic fields using an earthquake catalog from 1973 to 2023, which includes 9,016 earthquakes in the Michoacán-Guanajuato volcanic field and 841 in the Sierra Chichinautzin volcanic field. The spatial distribution of the hypocenters does not highlight any trend that could be associated with superficial movement of magma in the Sierra Chichinautzin volcanic field. In the Michoacán-Guanajuato volcanic field, however, eight seismic swarms since 1997 have been detected. These swarms are interpreted to result from ascending magma. Strengthening monitoring systems and reinforcing mitigation measures to address volcanic hazards and risk are important means of preparing for future eruptions in both regions. Analysis such as those herein provide insights as to where an eruption might occur and may help mitigate volcanic hazards.

Michoacán-Guanajuato, Sierra Chichinautzin↗

USGS1998

Technique for simulating peak-flow hydrographs in Maryland

The efficient design and management of many bridges, culverts, embankments, and flood-protection structures may require the estimation of time-of-inundation and (or) storage of floodwater relating to such structures. These estimates can be made on the basis of information derived from the peak-flow hydrograph. Average peak-flow hydrographs corresponding to a peak discharge of specific recurrence interval can be simulated for drainage basins having drainage areas less than 500 square miles in Maryland, using a direct technique of known accuracy. The technique uses dimensionless hydrographs in conjunction with estimates of basin lagtime and instantaneous peak flow. Ordinary least-squares regression analysis was used to develop an equation for estimating basin lagtime in Maryland. Drainage area, main channel slope, forest cover, and impervious area were determined to be the significant explanatory variables necessary to estimate average basin lagtime at the 95-percent confidence interval. Qualitative variables included in the equation adequately correct for geographic bias across the State. The average standard error of prediction associated with the equation is approximated as plus or minus (+/-) 37.6 percent. Volume correction factors may be applied to the basin lagtime on the basis of a comparison between actual and estimated hydrograph volumes prior to hydrograph simulation. Three dimensionless hydrographs were developed and tested using data collected during 278 significant rainfall-runoff events at 81 stream-gaging stations distributed throughout Maryland and Delaware. The data represent a range of drainage area sizes and basin conditions. The technique was verified by applying it to the simulation of 20 peak-flow events and comparing actual and simulated hydrograph widths at 50 and 75 percent of the observed peak-flow levels. The events chosen are considered extreme in that the average recurrence interval of the selected peak flows is 130 years. The average standard errors of prediction were +/- 61 and +/- 56 percent at the 50 and 75 percent of peak-flow hydrograph widths, respectively.

Maryland↗

USGS2024

Earthquake-triggered ground-failure inventory associated with the M7.1 2018 Southcentral Alaska earthquake

The 30 November 2018, magnitude (Mw) 7.1 earthquake in Southcentral Alaska triggered substantial landslides, liquefaction, and ground cracking throughout the region, resulting in widespread geotechnical damage to buildings and infrastructure. Despite a challenging reconnaissance and remote-sensing environment, we constructed a detailed digital inventory of ground failure associated with the event from several sources. Sources included information derived from remotely sensed data, and data compiled from literature, social media postings, and earthquake damage information compiled by local, state, and federal agencies. Each instance of ground failure within the inventory contains information on the location and type of observed ground failure, and the methods and data used to document the occurrence. Where high-quality data, such as LIDAR or satellite imagery, were available and showed the ground-failure instance clearly, the extent is mapped as a polygon or polyline. All other locations are mapped as points. There are a total of 886 ground-failure instances documented within the inventory (400 landslides, 286 liquefaction features, and 200 features unattributed to specific processes). A semi-quantitative confidence scheme is used to describe mapping certainty associated with each ground-failure feature. This inventory represents a relatively moderate ground-failure-triggering event that occurred in a subarctic environment. This data paper describes the content within the inventory, the inventory data collection procedures, and limitations of the data. Events of this type are not often documented in detail; thus, adding the inventory data to the US Geological Survey Open Repository of Earthquake-Triggered Ground-Failure Inventories further diversifies the datasets available to the scientific community to be used to better understand and model earthquake-triggered ground failure.

Alaska↗

USGS2026

Evaluating the U.S. Geological Survey’s earthquake shaking hazard forecasts and their implications on seismic risk

We analyze the last six update cycles of the long-term probabilistic earthquake shaking hazard forecast from the U.S. Geological Survey (USGS) and discuss the changes in hazard estimates from the 1996 to the latest 2023 update of the National Seismic Hazard Model (NSHM) for the conterminous U.S. We summarize how our understanding of earthquake hazards has evolved over the last two decades and quantify the implications of changing hazards on people, buildings, and lifeline infrastructure. The net effect of changes in hazard estimates, with reference to the mean hazard from the 2023 NSHM update that has 2% probability of exceedance in 50 years, suggests an overall increase in total geographic areas facing very strong shaking levels (modified Mercalli Intensity (MMI) of VII or more) when compared with any of the previous cycles. This increase puts ~159.2 million people (an increase from 111.8 million when using the 2018 NSHM), ~47.3 million residential buildings (an increase from 32.5 million with the 2018 NSHM), and ~$25.6 trillion of economic exposure of all buildings (an increase from $17.4 trillion when using the 2018 NSHM) at risk from earthquakes. We demonstrate that small changes in hazard estimates do not necessarily imply small changes in risk estimates. Therefore, changes in risk estimates can be used to highlight key USGS NSHM updates for improving risk mitigation efforts.

Earthquake Spectra Journal↗

USGS2026

Ringtail (Bassariscus astutus) survival in southwestern Oregon

The effective conservation and management of small carnivore populations requires understanding species’ life-history traits and identifying important vital rates that drive population trajectories. However, many of these species are rare or elusive and of state or federal conservation concern, and demographic information is often lacking and difficult to obtain. At the northern limit of their range in Oregon, USA, ringtail ( Bassariscus astutus ) occupy mid-elevation forests and are a species of conservation concern because of their limited distribution and suspected low density. We initiated a radio-telemetry study in 2020 to estimate monthly and annual survival of ringtail in southwest Oregon. We monitored 26 ringtail from November 2020 to October 2022 and estimated survival rates using a known-fate framework and Program MARK. Model-averaged monthly survival estimates ranged from lows of 0.963 (SE = 0.022, 95% CI = 0.887–0.988) in spring (Feb–May) of 2021 to highs of 0.980 (SE = 0.018, 95% CI = 0.889–0.997) in summer (Jun–Oct) of 2022. Model-averaged estimates of annual survival were 0.695 (SE = 0.176, 95% CI = 0.310–0.920) during 2020–2021 and 0.728 (SE = 0.168, 95% CI = 0.336–0.934) during 2021–2022. Predation was the leading cause of mortality, but notably, no mortalities could be attributed to avian predation. Survival rates were much higher within forested landscapes in Oregon than for populations in the southwestern United States, which may have implications for regional conservation and management strategies.

Oregon↗

USGS2025

Anatectic origin of Mississippian spodumene-bearing pegmatites in western Maine during orogenic plateau collapse

Spodumene pegmatites are an important lithium source, but the processes and tectonic settings in which they form are poorly understood. The Rumford pegmatite district surrounding Plumbago Mountain, western Maine, is host to numerous spodumene pegmatites, including the Plumbago North pegmatite (a world-class spodumene resource). Competing petrogenetic models for these spodumene pegmatites include (1) highly fractionated melts of the Mooselookmeguntic igneous complex and (2) anatexis. We tested these hypotheses by constraining the geologic, magmatic, metamorphic, and tectonic history of the Plumbago Mountain area with detailed geologic mapping and U-(Th)-Pb geochronology. The Silurian Rangeley Formation records initial isoclinal folding prior to, and contact-related metamorphism synchronous with, the intrusion of the 417 ± 4 Ma Plumbago Mountain pluton. Peak amphibolite facies metamorphism and crustal melting occurred during the ca. 410 to 400 Ma Acadian orogeny. Pulsed emplacement of the Mooselookmeguntic igneous complex occurred between ca. 389 and 356 Ma. Cassiterite U-Pb dates of spodumene pegmatites (333–327 Ma) are ≥23 m.y. younger than nearby granitic plutons, strongly arguing against the fractional crystallization model. Metamorphic monazite and xenotime (346–328 Ma) and 330 to 308 Ma 40 Ar/ 39 Ar hornblende dates indicate metamorphism coeval with spodumene pegmatite emplacement, supporting anatectic models. Reheating, anatexis, and spodumene pegmatite emplacement occurred during collapse of the 380 to 330 Ma Acadian orogenic plateau. Lithium enrichment may be linked to one or more stages of partial melting of metasedimentary and plutonic rocks during the formation, tenure, and collapse of the Acadian altiplano and emphasizes the role of anatexis in producing spodumene pegmatites of economic significance.

Maine↗

USGS2000

Assessment of potential for natural attenuation of chlorinated ethenes and ethanes in ground water at a petrochemical reclamation site, Harris County, Texas

Redox conditions in the Numerous Sand Channels Zone beneath a petrochemical reclamation site in Harris County, Texas, range from sulfate reducing to methanogenic as indicated by the presence of methane in ground water and the range of molecular hydrogen concentrations. Assessment of the potential for reductive dechlorination using BIOCHLOR as a screening tool indicated conditions favoring anaerobic degradation of chlorinated organic compounds in the Numerous Sand Channels Zone. Evidence supporting reductive dechlorination includes apparently biogenic cis -1,2-dichloroethene; an increased ratio of 1,2-dichloroethane to 1,1,2-trichloroethane downgradient from the assumed contaminant source area; ethene and methane concentrations greater than background concentrations within the area of the contaminant plume; and a positive correlation of the ratio of ethene to vinyl chloride as a function of methane concentrations. The body of evidence presented in this report argues for hydrogenolysis of trichloroethene to cis -1,2-dichloroethene; of 1,1,2-trichloroethane to 1,2-dichloroethane; and of vinyl chloride to ethene within the Numerous Sand Channels Zone. Simulations using BIOCHLOR yielded apparent first-order decay constants for reductive dechlorination in the sequence tetrachloroethene trichloroethene cis -1,2-dichloroethene vinyl chloride ethene within the range of literature values reported for each compound and apparent first-order decay constants for reductive dechlorination in the sequence 1,1,2-trichloroethane 1,2-dichloroethane slightly greater than literature values reported for each compound along the upgradient segment of a simulated ground-water flowpath. Except for vinyl chloride, apparent rates of reductive dechlorination for all simulated species show a marked decrease along the downgradient segment of the simulated ground-water flowpath. Evidence for reductive dechlorination of chlorinated ethenes within the Numerous Sand Channels Zone indicates potential for natural attenuation of chlorinated ethenes. Reductive dechlorination of chlorinated ethanes apparently occurs to a lesser extent, indicating relatively less potential for natural attenuation of chlorinated ethanes. Additional data are needed on the concentrations and distribution of chlorinated ethenes and ethanes in individual fine sand intervals of the Numerous Sand Channels Zone. This information, combined with lower minimum reporting levels for future chloroethane analyses, might enable a more complete and quantitative assessment of the potential for natural attenuation at the site.

Texas↗

USGS2025

Earthquake stress-drop values delineate spatial variations in maximum shear stress in the Japanese forearc lithosphere

Earthquake stress drop (Δσ) may increase with depth and stress in the brittle lithosphere. However, the range of uncertainty in Δσ and the lack of constraints on absolute stress make it difficult to establish whether they are correlated. Here, we investigate Δσ dependence on depth and maximum shear stress ( τ max ) based on ~11 years of seismicity in the northeastern Japanese forearc following the 2011 Tohoku-Oki megathrust earthquake. We interpret Δσ estimates computed using both individual spectra and spectral-ratio methods and find that Δσ exhibits a clear depth dependence within the seismically active upper ~60 km of the forearc lithosphere ( ~ 0.8 MPa per 10 km). We further compare Δσ values with quantitative τ max estimates from finite-element models of force balance. We find that median Δσ values increase with τ max in the brittle forearc lithosphere and that earthquake stress release is proportional to τ max . The dependence of Δσ on τ max explains the apparent depth dependence of Δσ and suggests that average Δσ values provide a relative measure of the stress at failure. In the northeastern Japanese forearc, Δσ values remained roughly constant in the decade following the Tohoku-Oki earthquake, suggesting negligible changes in failure stress in the forearc since the mainshock.

Communications Earth and Environment↗

USGS2026

Climate-adaptive urban planning: Quantitative assessment of drought impacts and practical strategies for climate-resilient urban green spaces

Urban green spaces (UGSs) are vital for enhancing a city’s resilience and livability; however, their functionality is increasingly jeopardized by drought, particularly in water-scarce regions. This study evaluates drought impact on UGSs in Metropolitan Adelaide, Australia, a representative semi-arid urban system, using satellite-derived Normalized Difference Vegetation Index (NDVI) time-series data spanning 2000–2020. Vegetation dynamics were analyzed through Seasonal-Trend decomposition using Loess (STL), standardized anomaly assessment, lagged Pearson correlation, Ordinary Least Squares (OLS) regression, and Mann–Kendall trend analysis. To isolate climatically sensitive signals, 29 urban lawn patches were examined separately from mixed urban canopy, given their shallow root systems and direct dependence on surface moisture. NDVI declined by approximately 0.09 units during the Millennium Drought (2001–2009), with summer greenness deficits reaching 24% below the 20-year benchmark. Temperature was the dominant driver of lawn NDVI variability (r = −0.863, R 2 = 74.5%), substantially exceeding the effect of rainfall (r = 0.156, R 2 = 2.4%). El Niño–Southern Oscillation (ENSO) cycles modulated vegetation responses, with La Niña years supporting recovery and El Niño years amplifying decline. Post-drought recovery remained incomplete, with NDVI deficits of 8–20% persisting through 2020; full recovery was observed only in 2017, coinciding with the highest recorded summer rainfall. No significant directional trend was detected over the full study period (Mann–Kendall τ = 0.005, p = 0.908). These findings demonstrate that heat, rather than water limitation alone, is the primary driver of vegetation stress in urban systems, highlighting the benefits of integrated management strategies that address both warming and moisture deficits to sustain urban green infrastructure under future climate conditions. We introduce the concept of “urban greenery drought,” referring to a form of vegetation stress in managed urban landscapes where greenness is reduced primarily by elevated temperature and atmospheric demand despite water availability.

Adelaide↗

USGS2025

Macroseismology

In this chapter I discuss the use of so-called macroseismic data, i.e., reports of damage and other effects of shaking on humans and the built environment, to improve the characterization of earthquakes and the ground motions they produce. Macroseismic data are critical not only to investigate earthquakes that occurred before the start of the instrumental era in seismology, but are also, by virtue of the spatially rich nature of the data, of tremendous potential importance to investigate early instrumental and even recent earthquakes. Fueled by both the advent of on-line systems that now produce unprecedented volumes of macroseismic intensity data and by recognition of the importance of understanding key historical earthquakes, there has been a growing appreciation for the value of macroseismic data. It is, however, critical to understand both the unique challenges and the unique opportunities associated with these data. I review the evolution of thought and practices regarding analysis of macroseismic data, and the use of such data to improve the characterization of historic and early instrumental earthquakes.

Book chapter↗

USGS2026

Timescales of cumulate mobilization and mixing for the 1868 A.D. eruption of Mauna Loa, Island of Hawai‘i

The deadly 1868 A.D. eruption of Mauna Loa’s lower Southwest Rift Zone (Island of Hawai‘i) included a M7.9 earthquake and associated tsunami and landslides, demonstrating the severe hazards posed by Earth’s largest active subaerial volcano. To better understand the relationship between intense seismic activity, dike emplacement, magma storage, transport histories, and mobilization of olivine cumulates at Mauna Loa, we examine compositional zoning of olivine in the 1868 lava flows. Samples range from basalt (< 10% olivine) to picrite (30–40% olivine). The olivine cargo is heterogeneous (Fo 78.2–89.2 ; forsterite = [Mg/(Mg + Fe) × 100]) but dominated by ~ Fo 89 cores that lie above the Fe-Mg equilibrium field of host glasses. Crystal rims < Fo 80 are due to post-eruptive modification in slow cooling lava flows. Minor element compositions fall within the range of other Mauna Loa olivine erupted in the past 200 years. Olivine crystals exhibit both normal and complex Fo zoning patterns that yield timescales of diffusive re-equilibration that range from 3 to 258 days, with 72% of crystals recording 71 days or less. These timescales correspond to magmatic priming of the summit reservoir system ~ 2 months prior to the eruption and the M7.9 earthquake likely facilitated the transport of the crystal-rich summit-derived magmas downrift shortly prior to eruption. If the recently proposed faster Fe-Mg diffusion coefficient is used, timescales instead range from < 1 day to 25 days, with most recording 1 week or less. In this scenario, most of the olivine zoning would have to have been generated after the M7.9 earthquake perturbed the system.

Hawaii↗