USGS Science⌕ Search

SEARCH · USGS Science

Search USGS Science

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.

At least 19 recordsLinked to original sources
USGS2028

Migration ecology and stopover population size of Red Knots Calidris canutus rufa at Mingan Archipelago after exiting the breeding grounds

Populations of migratory birds present unique conservation challenges given the often vast distances separating critical resources throughout the annual cycle. Migration areas close to the breeding grounds represent a link between two key stages of the annual cycle, and understanding migration ecology as birds exit the breeding grounds may be particularly informative for successful conservation. We studied migration phenology and stopover ecology of an endangered subspecies of the Red Knot Calidris canutus rufa at a migration area relatively close to its breeding range. Using mark-recapture/resight data and a Jolly-Seber model for open populations, we described the arrival and departure schedules, stopover duration, and passage population size at the Mingan Archipelago, Quebec, Canada. Red Knots arrived at the study area in two distinct waves of birds separated by approximately 22 days. Nearly 30% of the passage population arrived in the first wave of arrivals during 15–18 July, and approximately 22% arrived in a second wave during 8–11 August. The sex-ratio in the stopover population at the time of the first wave was slightly skewed toward females, whereas the second wave was heavily skewed toward males. Because males remain on the breeding grounds to care for young, this may reflect successful breeding in the year of our study. The estimated stopover duration (population mean) was 11 days (95% credible interval: 10.3–11.7 days), but stopover persistence was variable throughout the season. We estimated a passage population size of 9,450 birds (8,355–10,710), a minimum estimate for reasons related to the duration of our sampling. Mingan Archipelago is thus an important migration area for this endangered subspecies and could be a priority in conservation planning. Our results also emphasize the advantages of mark-recapture/resight approaches for estimating migration phenology and stopover persistence.

Quebec↗
USGS2027

Ground failure resulting from earthquake ground motions

The ground failure resulting from an earthquake, at all scales, is ultimately controlled by the details of near-surface earthquake ground motions, the spatially distributed mechanical characteristics of surface materials and how the two interact. This chapter introduces the aspects of ground-motion seismology that are most relevant to the triggering of different types of ground failure, such as landslides and liquefaction. In particular, the chapter focuses on different sources of local modification (e.g. amplification) of earthquake ground motions in the shallow subsurface, known as site effects. The chapter ends with an overview of the controls on, general spatial patterns of, and approaches to regional modelling of the manifestation of ground failure that may be triggered by an earthquake.

Book chapter↗
USGS2026

Distribution, abundance, and habitat characteristics of Coastal Cactus Wrens ( Campylorhynchus brunneicapillus ) in San Diego County, California—2025 data summary

Executive Summary We surveyed for coastal Cactus Wren ( Campylorhynchus brunneicapillus ) in 498 established plots in San Diego County, California, in 2025, encompassing 4 genetic clusters (Otay, Lake Jennings, Sweetwater/Encanto, and San Pasqual). Of the 498 plots, 372 have been surveyed every year starting in 2020. Two surveys were completed at each plot between March 1 and July 31. Cactus Wrens were detected in 220 plots (44 percent of plots). In plots that have been consistently surveyed since 2020, Cactus Wrens were detected in 35 percent of plots, similar to 2024 (34 percent), 2021 (35 percent), and 2020 (35 percent), and higher than in 2023 (26 percent) and 2022 (30 percent). Across all survey plots we detected 184 Cactus Wren territories in 2025. We documented 107 territories in plots that have been consistently surveyed since 2020, which is similar to the number of territories in 2024 (106) and 2020 (108), more than we documented in 2022 (93) and 2023 (85), but slightly less than the peak in 2021 (113). Across all survey plots the number of territories increased from 2024 to 2025 in the San Pasqual genetic cluster, increased slightly in the Lake Jennings cluster, decreased slightly in the Sweetwater/Encanto genetic cluster and remained the same in the Otay genetic cluster. At least 83 percent of Cactus Wren territories were occupied by pairs, and 100 fledglings were observed in 2025. We observed three banded Cactus Wrens in 2025. Adult birds ranged from 6 to 9 years old and occupied the same territories as in previous years. Banding of Cactus Wrens ceased in 2019 and as a result we have a very low number of banded wrens remaining in the population, which makes it difficult to detect movement of wrens between genetic clusters. We detected no movement of banded Cactus Wrens between genetic clusters in 2025. Vegetation at Cactus Wren survey plots was dominated by coastal sage scrub shrubs, such as California sagebrush ( Artemisia californica ), California buckwheat ( Eriogonum fasciculatum ), and lemonade berry ( Rhus integrifolia ). Cactus Wrens preferentially selected plots dominated by California buckwheat in the Otay genetic cluster and avoided plots dominated by lemonade berry and jojoba ( Simmondsia chinensis ). Cactus Wrens also appeared to prefer plots dominated by broom baccharis ( Baccharis sarothroides ) while avoiding plots dominated by California sagebrush in the Lake Jennings cluster. No definitive signs of fungal pathogens were observed on cactus within and around survey plots. Blue elderberry ( Sambucus mexicana ) was detected at 43 percent of plots, and Cactus Wrens preferentially occupied more plots with elderberry than were available. Very little dead or unhealthy cactus was observed within all survey plots, and Cactus Wrens preferentially occupied plots with less dead cactus while occupying plots with unhealthy cactus in proportion to its availability. Eighty-nine percent of plots had more than 5 percent of cactus crowded or overtopped by vines and shrubs, and Cactus Wren occupied plots in proportion to the degree of shrub and vine crowding or overtopping. Just over half of the plots had greater than 5 percent non-native annual cover, with Cactus Wrens preferentially occupying plots with less non-native cover.

California↗
USGS2026

The National Wildlife Disease Database

The National Wildlife Disease Database (NWDD) is a national information hub that integrates wildlife disease data across jurisdictions to guide natural resource management and deliver early warnings of unusual events and emerging threats, strengthening public health preparedness and protecting the agricultural economy.

Fact Sheet↗
USGS2026

Evolution of summit magma pressure during the 2018 caldera collapse at Kılauea volcano

We utilized geodetic data from the 16 May to 4 August 2018 caldera collapse sequence to refine the geometries and pressure evolution of two summit reservoirs at Kīlauea volcano: the shallower Halema‘uma‘u (HMM) and deeper South Caldera (SC) reservoirs. We modeled magma reservoirs as pressurized spheroidal cavities and constrained reservoir geometries with Global Positioning System (GPS) and Interferometric Synthetic Aperture Radar displacements. Inverting high-rate GPS and borehole tilt time series recovered time-dependent reservoir pressure changes and caldera ring fault shear stress. During the collapse sequence, the HMM reservoir accommodated episodic, co-collapse pressure increases of 1.9 MPa (median) and a net 5.8 MPa (median) decrease. Simultaneously, the SC reservoir underwent monotonic depressurization of 6.7 MPa (median). Negligible co-collapse pressure changes in the SC reservoir imply it was not exclusively supplying magma to the eruption vent. These results provide a coherent view of magma reservoir pressure evolution during the 2018 eruption.

Hawaii↗
USGS2026

Slip tendency and dilation tendency of Quaternary faults in the Great Basin region, Arizona, California, Idaho, Nevada, Oregon, Utah, and Wyoming

This report presents the results of slip tendency and dilation tendency analysis of Quaternary faults in the Great Basin region, in the states of Arizona, California, Idaho, Nevada, Oregon, Utah, and Wyoming. This effort was undertaken as part of two studies: the U.S. Department of Energy Geothermal Technologies Office-funded Innovative Geothermal Exploration through Novel Investigations of Undiscovered Systems (INGENIOUS) project, which is a multidisciplinary, multi-institution effort to develop new methodologies and best practices to accelerate the discovery of new, commercially viable hydrothermal resources; and the U.S. Geological Survey Geothermal Resources Investigations Project. The goal of this work is to identify faults and fault segments that are appropriately oriented to be stress-loaded for slip or to dilate under the ambient stress conditions. Both conditions may make such faults more likely to host permeability and fluid flow relative to other faults in the region. Slip tendency and dilation tendency were calculated for all fault segments in the study area from the U.S. Geological Survey Quaternary Fault and Fold Database. Results demonstrate that faults with high slip tendency and dilation tendency are common in the Great Basin region. The areas containing many high slip tendency and dilation tendency fault segments are associated with high strain rates. These areas also host a relatively high proportion of the known hydrothermal systems. This indicates that faults with high slip tendency and dilation tendency in areas with high strain-rate are potentially most prospective for hosting hydrothermal systems. The database published in coordination with this report contains the slip tendency and dilation tendency results data, including the fault geometry and stress data used to calculate these quantities, for both the Great Basin study area and the smaller INGENIOUS study area.

Arizona, California, Idaho, Nevada, Oregon, Utah, ↗
USGS2026

Water and endangered fish in the Klamath River basin revisited: Does old age or Upper Klamath Lake surface elevation and water quality affect adult Lost River Sucker and Shortnose Sucker survival? A response to Kann and Walker

Our update responds to recent scientific critiques and again finds no evidence that lake levels or water conditions influence adult sucker survival. This reaffirms our prior work (Krause et al., 2022) and suggests that recent declines in adult survival may reflect an aging population.

Oregon↗
USGS2026

Characterization of subsurface groundwater pathways from treated wastewater effluent infiltration along the South Rim of Grand Canyon National Park, Arizona, USA

Recharge of treated wastewater to the environment is increasingly becoming a component of water resources management, especially in water-limited areas such as the Grand Canyon. Treated wastewater from the Grand Canyon South Rim Wastewater Treatment Plant has been discharged to a dry wash located along the surface expression of a fault for decades, but the fate of the wastewater is unclear. Previous research suggested possible connections between the treated wastewater discharge sites and springs along the South Rim of Grand Canyon, particularly Monument spring. Additional investigations presented in this report add further evidence to the connection between the wastewater and Monument spring, utilizing per- and polyfluoroalkyl substances, artificial sweeteners, and stable isotopes of water. A dye trace established the first direct connection between wastewater and Monument spring where the dye added to treated wastewater took about 4–6 months to move through the subsurface and was detectable for 5 months. The persistence of wastewater tracers over time and the timing of dye movement suggests a fast subsurface flow pathway. On the basis of geologic observations of faults and jointing, groundwater flow follows a possible combination of karst and other structural pathways. The dye and wastewater tracers were monitored at other springs along the South Rim of Grand Canyon and the combination of repeated dye detections and multiple wastewater tracers were only found at Monument spring. Multiple lines of evidence provide a robust tracer suite to support evidence of structural/karst controls on preferential flowpaths; faults alone are unlikely to explain the flow pathway.

Arizona↗
USGS2026

Groundwater nutrient contributions to surface water in the upper Maumee River basin

Harmful algal blooms triggered by nutrient loading have negatively affected Lake Erie water quality for decades, but the role groundwater plays in transporting nutrients to tributary streams is unclear. Here, we apply a multi-scale approach to investigate groundwater nitrogen (N) and phosphorus (P) dynamics at the field scale, along a groundwater flow path, and at the Maumee River basin scale. Unsaturated zone monitoring indicates that macropore flow is the primary mechanism for N and P leaching with nitrate and orthophosphate (OP) concentrations from macropore samplers more than double those from soil-matrix samplers. Shallow (< 4.5 m) groundwater nitrate concentrations indicate the timing of saturated-zone nitrate fluxes is crop dependent with more nitrate mobilized during the maize growing season whereas concentrations below soybeans can remain elevated (> 3 mg/L) through winter months. A three-end-member hydrograph analysis for the upper Maumee River basin indicates that tile drains transport most subsurface N and P with greater than 15,000 kg per day of nitrate transported between January and May. Orthophosphate concentrations in tile drainage at an edge-of-field site correlate well with concentrations in a nearby tributary of the Maumee River. While runoff is the dominant source of OP loading, tile drainage and groundwater discharge contributed an estimated 36% of OP load during a 2018–2022 analysis period with most subsurface loading occurring between January and March. These findings highlight how best management practices (such as cover crops and tile drainage control structures) can be used to help limit leaching and tile effluent during the nongrowing season.

Indiana↗
USGS2026

Blind spots in our understanding of the tropical nitrogen cycle

Understanding environmental controls on tropical nitrogen (N) cycling is critical for projecting flux rates and modeling ecosystem responses to rising CO 2 . Here, we ask whether geographic bias in field sampling may shape our understanding of N cycling by failing to capture the breadth of environmental heterogeneity within the tropics. We mapped empirical sampling locations of two major ecosystem N inputs and outputs—biological nitrogen fixation and denitrification—paired with a wide range of variables thought to influence N cycling. Both single and multi-variate analyses revealed that sampling density is greater in productive, wet, fertile, warmer tropical ecosystems with low seasonality—areas expected to have higher process rates—but biased away from spatially extensive biomes which may have high process rate variability and thus higher uncertainty. Even within tropical broadleaf forests, sampling is biased toward sites with unusually high soil fertility (N and phosphorus), potentially skewing biome-wide averages. Although some geographically undersampled regions share environmental characteristics with well-studied sites, others are unique, compromising inference from existing data. Finally, we find that while sampled field locations capture the range of N cycle process means/variances predicted by an ensemble of terrestrial biosphere models, this may be the result of circularity in model parameters and driving data rather than true coverage, including overrepresentation of forest plant functional types that may be parameterized from a small subset of sites. Data coverage expansion into areas with novel environmental conditions and higher uncertainty could be used to assess and improve model parameterization, as well as test general hypotheses.

Global Biogeochemical Cycles↗
USGS2026

EMIT reveals elevated calcite signatures in shorelines of Lake Powell and Lake Mead exposed by drought and Colorado River management

Nearly half of Earth’s seasonal freshwater storage variability occurs in human-managed reservoirs, resulting in large areas of reservoir-bed that are episodically inundated or subaerially exposed. Here we use hyperspectral data from NASA’s Earth Surface Mineral Dust Source Investigation (EMIT) mission to map calcite mineral assemblages across exposed surfaces in two large U.S. reservoirs, Lake Mead and Lake Powell. Results show distinct patterns in the occurrence of carbonate, clay, and mica minerals. Reservoir-deposited sediments exhibited significantly higher calcite proportions than adjacent, lithologically similar, non-inundated terrestrial surfaces. These differences likely reflect hydrologically-driven calcite precipitation, wherein warm, clear, productive reservoir waters drive calcium carbonate out of solution. We also detect reductions in calcite proportions in areas of the shoreline disturbed by human recreation. EMIT's imaging capabilities enable rapid, high-resolution mineralogical inventories across vast areas where field sampling is infeasible. We find EMIT is a scalable tool for monitoring geochemical changes associated with large water level fluctuations.

Arizona, Nevada, Utah↗
USGS2026

Drivers of ecosystem metabolism in an anthropogenically modified delta

Estuaries are under threat globally from a wide range of stressors, with rates of ecosystem metabolism – respiration (R), gross primary production (GPP), and net ecosystem metabolism (NEM) – representing a way for resource managers to quantify their impacts. In this study, we assessed how natural and anthropogenic stressors have altered metabolism within the extensively modified Sacramento-San Joaquin Delta (“Delta”) in California, USA. Metabolism was estimated using a diel open-water technique utilizing dissolved oxygen measurements from a high frequency (15 min) water quality monitoring array, and spanned 6 sites across 11 years (2014–2025). Our analyses were driven by a conceptual model of Delta carbon (C) cycling focused on 1) Invasive aquatic vegetation (IAV), 2) Phytoplankton, 3) Wastewater treatment plant (WWTP) upgrades, 4) Large-scale flow alterations, and 5) Wet vs. dry years. We found that allochthonous C fuels net heterotrophy in the outer riverine sections of the Delta, while the heavily vegetated central Delta produces autochthonous C remineralized elsewhere in the system. IAV were the dominant primary producers, with low phytoplankton production except during a historic bloom. WWTP upgrades imparted negligible changes to metabolism, with low chlorophyll a (< 5 µg L −1 ) both before and after. Introduction of a flow barrier – a management response to drought – hydrodynamically isolated deep-water channels from nearby shallow areas, resulting in decreased metabolism indicative of changes to C transport. Contrary to past research, low rainfall years were associated with increased metabolism. This study demonstrates how multi-sensor arrays can be used to quantify ecosystem responses to major environmental stressors.

California↗
USGS2026

Microgravity measurements spanning 2019–2022 indicate no change in magma supply or storage during Kilauea's reawakening

After more than 2 years of quiescence following its major 2018 flank lava effusion and summit caldera collapse, Kīlauea, Hawaiʻi, returned to frequent eruptive activity in late 2020. Microgravity measurements from the summit region in 2019 and 2022 span this reawakening and indicate a gravity increase of up to hundreds of microgals near the center of the summit caldera. Much of this increase is caused by lava that partially filled the 2018 collapse pit. We accounted for this mass by differencing digital elevation models of the summit area to determine the geometry and amount of accumulated lava and its gravity contribution at each station in the network. The remaining gravity change is best approximated by a mass increase of 6.9 × 10 10 kg (5.1–23.0 × 10 10 kg at 95% credible interval) at a depth of 480 m (121–795 m at 95% credible interval) beneath the locus of the 2018 collapse. The modeled source depth reflects the combined impacts of magma accumulation in a reservoir known from previous geophysical studies to exist at about 1–2 km depth, intrusions between that reservoir and the surface, and mass accumulation in the heavily fractured caldera floor. The spatial and temporal patterns of subsurface mass accumulation during 2019–2022 imply that there were no significant changes in magma reservoir geometry or magma supply to the volcano in the years following the 2018 collapse.

Hawaii↗
USGS2026

How Indigenous socio-cultural systems resist, accept, or direct ecological transformation: A research agenda

Anthropogenic ecological transformations are challenging conventional environmental management paradigms, as maintaining historical ecological structures and functions becomes increasingly untenable. The resist–accept–direct (RAD) framework offers a strategic tool to navigate these shifts. However, as a product of the dominant socio-cultural system, the RAD framework risks perpetuating the same patterns of thinking that have contributed to unsustainable social-ecological dynamics. Indigenous socio-cultural systems, considered here through the lens of long-term coevolution with their environments, offer valuable empirical models for adaptive responses to environmental change. Such Indigenous socio-cultural responses, grounded in narratives, governance institutions, value systems, and other cultural elements collectively referred to as Indigenous Knowledge (IK), hold relevance for addressing today’s intertwined environmental and social challenges. We provide a targeted review of existing literature at the intersection of IK and the RAD framework, along with a broader conceptual overview of relevant insights across related fields. These reveal contrasting assumptions about Indigenous responses to environmental change, such as consistently resisting to preserve cultural fidelity or uniformly accepting, reflecting deeper cultural, temporal, and ecological complexity. We propose a structured research agenda to investigate Indigenous socio-cultural responses to environmental change through empirical research, theoretical synthesis, hypothesis development, and practical application. This agenda is grounded in comprehensive collaboration, emphasizing shared power and integrating both emic and etic perspectives throughout all stages of research and implementation. By bridging the adaptive dynamics of coevolved social-ecological systems and mainstream management strategies, an IK–RAD framework could enhance the efficacy of responses to global environmental change.

Human Ecology↗
USGS2026

Phase evolution and elemental partitioning in impact melts from Meteor Crater, Arizona

The impact of the Canyon Diablo iron meteorite in Arizona ∼61,000 years ago resulted in partial vaporization of the meteorite and target rocks, producing a vapor plume and melt particles with taenite, olivine and clinopyroxene in glassy matrix. The resulting Meteor Crater is ∼1.2 km in diameter and ∼200 m deep. Although the crater has been the natural laboratory for impact cratering, the microstructure and composition of its impact-melt particles remain poorly studied. Here microstructure and composition were investigated using optical and scanning electron microscopy, energy-dispersive X-ray spectroscopy, electron backscatter diffraction, and synchrotron microdiffraction on particles collected from the ejecta blanket to reconstruct their temporal and mineralogical evolution during ultrafast quenching and to explore how mafic minerals form from sedimentary-derived melts. Elemental mapping reveals complex chemical partitioning among glass and mafic phases. Formation of Fe-Ni globules with face-centered cubic (fcc) taenite indicates oxygen fugacity well below atmospheric in the vapor cloud. During rapid cooling, skeletal olivine crystallized, followed by nucleation of clinopyroxenes to form spherical clusters. The glassy matrix contains two glasses separated by sharp boundaries, likely representing immiscible melts derived from different lithologies. The hot particles were deposited on the surface among quartz- and dolomite-rich ejecta, followed by post-depositional precipitation of calcite and barite within glassy nodules. The results demonstrate that mafic minerals can crystallize from melts derived largely from sedimentary target rocks, with Fe supplied by the impactor, and that elemental partitioning and mineral textures record strong disequilibrium crystallization during rapid cooling.

Arizona↗
USGS2026

Map depicting locations susceptible to landslide hazards triggered by prolonged intense rainfall, Utuado Municipio, Puerto Rico

Introduction Landslides commonly occur in Puerto Rico during or soon after rainstorms with prolonged intense rainfall. Rainfall-triggered landslides present significant hazards to the built environment and human safety (for example, Jibson, 1986; Larsen and Torres Sanchez, 1992; Bessette-Kirton and others, 2019a). These hazards were illustrated during September 2017 (Hughes and others, 2019) when Hurricane María triggered more than 70,000 landslides across Puerto Rico (fig. 1), causing loss of life and widespread destruction (figs. 2–4). In Utuado Municipio, about 18,000 landslides were observed after Hurricane María (Hughes and others, 2019). Most intense rainfall-triggered landslides in Puerto Rico occur in soil overlying bedrock and travel downslope from where they initiate (for example, Jibson, 1986; Bessette-Kirton and others, 2019a, 2020). Hazards to human life and infrastructure can differ substantially between where landslides initiate and where they travel. At landslide initiation locations, the ground breaks up (fails) as it moves downslope, removing support for buildings, roadways, and other structures located upon the landslides (fig. 3A, B). This ground failure often causes structures to collapse and suffer severe damage or complete destruction, rendering building occupants at extreme risk. Roadways and other vital infrastructure, such as water, power, and communication lines, are often severed. During intense rainfall events, landslide initiation and movement downslope can occur in seconds to hours. Downslope from where they initiate, landslides in Puerto Rico often move very rapidly (faster than a person can walk or run) and can impact and destroy or fill structures with mud and debris (fig. 4) as well as bury roadways (fig. 3A, C). Moving landslides present extreme hazards to human safety; some fatalities during Hurricane María resulted from mud and debris that flowed through windows and doors of a home and buried occupants. Many landslides in Puerto Rico become fast-moving liquefied masses (debris flows), which are especially common in channels where they can quickly travel kilometers downstream. Multiple debris flows may occur in a channel during a single rainstorm (some channels experienced more than 50 debris flows during Hurricane María; Bessette-Kirton and others [2020] ). This hazard map separately depicts locations and the relative likelihood of where landslides may initiate, where they may move downslope, and where channelized debris flows may travel during or soon after intense rainfall like that during Hurricane María. An estimated average of 399 millimeters of rainfall fell across Utuado during September 19–21, 2017 (Ramos-Scharrón and Arima, 2019). Lower rainfall amounts, durations, and intensities also trigger landsliding in Puerto Rico, although landslides and debris flows become more numerous as rainfall intensities, durations, and cumulative amounts increase. The map does not show potential landslide locations in the karst terrain of northern Utuado (Bawiec, 1999) because Hurricane María-induced landslides were very sparse in this region (Hughes and others, 2019; Baxstrom and others, 2021a), even though rainfall amounts were similar to those in other parts of Utuado. Based on data from Hughes and others (2019), landslides were about 10 times more likely to occur outside of the karst terrane than within it, and our preliminary analyses indicated relatively low susceptibility to rainfall-induced landslide hazards within the karst. Landslide-prone soils appear to be poorly developed on the very steep slopes present in the karst region, where weathering by dissolution appears to dominate.

Puerto Rico, Utuado Municipio↗
USGS2026

Map depicting locations susceptible to landslide hazards triggered by prolonged intense rainfall, Naranjito Municipio, Puerto Rico

Introduction Landslides commonly occur in Puerto Rico during or soon after rainstorms with prolonged intense rainfall. Rainfall-triggered landslides present significant hazards to the built environment and human safety (for example, Jibson, 1986; Larsen and Torres Sanchez, 1992; Bessette-Kirton and others, 2019a). These hazards were illustrated during September 2017 (Hughes and others, 2019), when Hurricane María triggered more than 70,000 landslides across Puerto Rico (fig. 1), causing loss of life and widespread destruction (figs. 2–4). In Naranjito Municipio, about 2,000 landslides were observed after Hurricane María (Hughes and others, 2019). Most intense-rainfall-triggered landslides in Puerto Rico occur in soil overlying bedrock and travel downslope from where they initiated (for example, Jibson, 1986; Bessette-Kirton and others, 2019a; 2020). Hazards to human life and infrastructure can differ substantially between where landslides initiate and where they travel. At landslide initiation locations, the ground breaks up (fails) as it moves downslope, removing support for buildings, roadways, and other structures located upon the landslides (figs. 3B, 4A, B). This ground failure often causes structures to collapse and suffer severe damage or complete destruction, rendering building occupants at extreme risk. Roadways and other vital infrastructure, such as water, power, and communication lines are often severed. During intense rainfall events, landslide initiation and movement downslope can occur in seconds to hours. Downslope from where they initiate, landslides in Puerto Rico often move, or runout, very rapidly (faster than a person can walk or run) and can impact and destroy or fill structures with mud and debris (figs. 3A, 4A, C) as well as bury roadways (fig. 4A). Moving landslides present extreme hazards to human safety; some fatalities during Hurricane María resulted from mud and debris that flowed through windows and doors of a home and buried occupants. Many landslides in Puerto Rico become fast-moving debris flows, which are especially common in channels where they can quickly travel kilometers downstream. Multiple debris flows may occur in a channel during a single rainstorm (some channels experienced more than 50 debris flows during Hurricane María; Bessette-Kirton and others, 2020). This map separately depicts locations and the relative likelihood of where landslides may initiate, where they may move downslope, and where channelized debris flows may travel during or soon after intense rainfall like that during Hurricane María. An estimated average of 529 millimeters of rainfall fell across Naranjito during September 19–21, 2017 (Ramos-Scharrón and Arima, 2019). Lower rainfall amounts, durations, and intensities also trigger landsliding in Puerto Rico, although landslides and debris flows become more numerous as rainfall intensities, durations, and cumulative amounts increase.

Naranjito Municipio, Puerto Rico↗
USGS2026

Map depicting locations susceptible to landslide hazards triggered by prolonged intense rainfall, Lares Municipio, Puerto Rico

Introduction Landslides commonly occur in Puerto Rico during or soon after rainstorms with prolonged intense rainfall. Rainfall-triggered landslides present significant hazards to the built environment and human safety (for example, Jibson, 1986; Larsen and Torres Sanchez, 1992; Bessette-Kirton and others, 2019a). These hazards were illustrated during September 2017 (Hughes and others, 2019), when Hurricane María triggered more than 70,000 landslides across Puerto Rico (fig. 1), causing loss of life and widespread destruction (figs. 2 and 3). In Lares Municipio, about 3,500 landslides were observed after Hurricane María (Hughes and others, 2019). Most intense-rainfall-triggered landslides in Puerto Rico occur in soil overlying bedrock and travel downslope from where they initiated (for example, Jibson, 1986; Bessette-Kirton and others, 2019a; 2020). Hazards to human life and infrastructure can differ substantially between where landslides initiate and where they travel. At landslide initiation locations, the ground breaks up (fails) as it moves downslope, removing support for buildings, roadways, and other structures located upon the landslides (figs. 2, 3A, B). This ground failure often causes structures to collapse and suffer severe damage or complete destruction, rendering building occupants at extreme risk. Roadways and other vital infrastructure, such as water, power, and communication lines are often severed. During intense rainfall events, landslide initiation and movement downslope can occur in seconds to hours. Downslope from where they initiate, landslides in Puerto Rico often move, or runout, very rapidly (faster than a person can walk or run) and can impact and destroy or fill structures with mud and debris (figs. 3A, C) as well as bury roadways (fig. 3A). Moving landslides present extreme hazards to human safety; some fatalities during Hurricane María resulted from mud and debris that flowed through windows and doors of a home and buried occupants. Many landslides in Puerto Rico become fast-moving debris flows, which are especially common in channels where they can quickly travel kilometers downstream. Multiple debris flows may occur in a channel during a single rainstorm (some channels experienced more than 50 debris flows during Hurricane María; Bessette-Kirton and others, 2020). This map separately depicts locations and the relative likelihood of where landslides may initiate, where they may move downslope, and where channelized debris flows may travel during or soon after intense rainfall like that during Hurricane María. An estimated average of 256 millimeters of rainfall fell across Lares during September 19–21, 2017 (Ramos-Scharrón and Arima, 2019). Lower rainfall amounts, durations, and intensities also trigger landsliding in Puerto Rico, although landslides and debris flows become more numerous as rainfall intensities, durations, and cumulative amounts increase. The map does not show potential landslide locations in the karst terrain of northern Lares (Bawiec, 1999) because Hurricane María-induced landslides were very sparse in this region (Hughes and others, 2019; Baxstrom and others, 2021a), even though rainfall amounts were similar to those in other parts of Lares. Based on data from Hughes and others (2019), landslides were about ten times more likely to occur outside of the karst terrane than within it, and our preliminary analyses indicate relatively low susceptibility to rainfall-induced landslide hazards within the karst. Landslide-prone soils appear to be poorly developed on the very steep slopes present in the karst region, where weathering by dissolution appears to dominate.

Lares Municipio, Puerto Rico↗