Geology topics

Research about Puerto Rico

Source-linked reports with geographic coverage including Puerto Rico.

At least 19 recordsLinked to original sources

USGS2026

Characterizing metabolic responses of Eleutherodactylus frogs in Puerto Rico to different thermal treatments: Implications for conservation

Amphibians are increasingly vulnerable to extinction due to physiological stress induced by projected global warming and associated drying trends. Understanding the underlying mechanisms and consequences of these responses is crucial to developing effective conservation strategies. Puerto Rico hosts 15 endemic Eleutherodactylus frog species, which are particularly vulnerable to environmental changes due to poor vagility and narrow thermotolerance. In this study, we report routine metabolic rates (RMRs) of four species in response to four temperature treatments (15, 20, 25 and 30°C) at two temporal scales (24 and 72 hours) using indirect calorimetry. Results indicated that three species exhibited increased RMR with rising temperatures. The RMR at 25°C for E. wightmanae was spurious, and it was the only species that did not survive at 30°C. The Arrhenius breakpoint analysis indicated that E. juanariveroi RMR decreased between 25 and 30°C. In contrast, E. antillensis and E. coqui RMR increased with increasing assay temperatures. We discuss how diel metabolic demands decrease the scope at both thermal extremes for E. antillensis and E. coqui , and how habitat temperature may impinge energetic challenges to those coping with documented warming trends, such as E. wightmanae .

Puerto Rico↗

USGS2026

Flood-inundation maps for Río Grande de Arecibo in and near Utuado Municipio, Puerto Rico, 2026

The U.S. Geological Survey (USGS) created digital flood-inundation maps for a 3.4-mile reach of Río Grande de Arecibo in the municipality of Utuado Municipio, Puerto Rico. Water-surface profiles were computed for the stream reach by using a one-dimensional, steady-state, step-backwater model. This model was calibrated to the current (2026) stage-streamflow relation (rating curve number 11.0) for the USGS streamgage 50024950, Río Grande de Arecibo below Utuado, Puerto Rico. The resulting hydraulic model was used to compute 11 water-surface profiles for water levels (flood stages) ranging from 14.00 to 34.00 feet at the streamgage and ranging from what are referred to as “minor flood stage” to above “major flood stage,” as reported by the National Weather Service. The 34.00-foot stage exceeds the peak stage (32.92 feet) recorded at this streamgage in 1998. The simulated water-surface profiles were used in combination with a digital elevation model derived from lidar data to map the inundated areas associated with each flood profile. The flood-inundation maps and the supporting hydraulic model produced by this study can be used by emergency managers and local officials to define flood hazard areas and assess flood mitigation strategies to help protect life and property, to coordinate flood response activities such as evacuations and road closures, and to aid post-flood recovery efforts. Information about map limitations, associated uncertainties, and how to view the maps and access the associated data online is included in this report.

Puerto Rico↗

USGS2026

Impacts of transient critical zone evolution on landslide susceptibility and hazard Across the Utuado batholith, Puerto Rico

Landscapes in humid climates often experience rapid chemical and mechanical weathering that can change groundwater fluxes, drainage pathways, and frictional strength via alterations in subsurface mineral composition and matrix porosity. These changes to hillslope critical zone structure and strength should control landslide susceptibility, size, and failure style. Here we test the impact of differential weathering on landslide failure processes by exploiting a natural experiment along the dioritic Utuado batholith in Puerto Rico's Cordillera Central, where a sharp weathering gradient follows an upstream-propagating knickzone into an uplifted shoreline platform that was affected by Hurricane Maria in 2017. Through a combination of landslide mapping, field geophysical and geotechnical surveys, topographic analysis, calibrated stream profile inversion, 3D back-analysis of hillslope strength, and debris-flow inundation assessments, we show that differential weathering induces changes to critical zone structure and consequently landslide susceptibility and hazard. Our mapping reveals that the hurricane-induced landslides preferentially cluster along less weathered hillslopes despite a marked increase in frictional strength. Jointly inverted seismic refraction and electrical resistivity profiles show that the frictionally weak, heavily weathered hillslopes are highly porous and dry, which we posit acts to stabilize slopes through limiting deeper hydrologic connectivity that can produce shallow water tables with destabilizing lateral seepage forces during storms. Debris-flow inundation mapping and analysis reveals that major contributors to debris-flow growth, channel steepness and landslide coalescence, scale linearly with hillslope friction angle. This leads to higher-mobility debris flows along knickzone reaches, which disproportionately affects communities at the base of the plateau.

Puerto Rico↗

USGS2026

3D Elevation Program—Puerto Rico

Using 3D Elevation Program data, a modern illustration of Puerto Rico’s landscape showcases its elevation differences in vivid detail. This map has been consistently adjusted and colorized according to the highest and lowest elevations, accentuating differences across Puerto Rico. To further emphasize landforms, shaded relief has been added on top of the elevation colors.

Puerto Rico↗

USGS2026

Fault sources of the 2020 southern Puerto Rico seismic sequence and the regional tectonic structure

The 2020 southwestern Puerto Rico seismic sequence occurred in an area not known for seismic activity and included multiple normal and strike‐slip focal mechanisms. The largest earthquake, M w 6.4, caused significant damage and permanent shoreline subsidence. The locations of the ruptured faults have been explored extensively using seismological and geodetic evidence, but no consensus has been reached. Here, we report the results of a high‐resolution seismic reflection survey that identified a north‐verging 20 km long curved normal fault spanning the upper reach of the Guayanilla canyon, which we call the Guayanilla Canyon fault. Part of the fault likely ruptured during the M w 6.4. The concave‐to‐the‐north fault bounds a south‐dipping ∼1 km deep asymmetric basin. This long‐lived structure, identified also on sediment isopach and magnetic anomaly maps, appears to control the formation of the Guayanilla canyon system, the only submarine drainage system in southern Puerto Rico, suggesting a relationship among seismic activity, geomorphology, and reef development. Our seismic reflection data do not support the existence of previously published major seismogenic faults in the area, such as the north‐northeast–south‐southwest Guayanilla fault zone, the Northern Boqueron Bay‐Punta Montalva fault, or a northwest–southeast strike‐slip fault, suggested to have ruptured as part of the M w 6.4 earthquake. We identified smaller faults in the survey area, some of which can be correlated with seafloor expressions, and/or with changes in sediment thickness and magnetic anomalies, and/or with moderate‐sized earthquakes. Multiple seismic images of Investigator fault, a prominent mid‐slope feature parallel to the south coast of Puerto Rico, show buried low‐angle thrusts and an abrupt change in the attitude of the overlying sediments, leading us to interpret this feature as the expression of a boundary between the now‐inactive upper accretionary prism of Muertos trough and a buried back‐arc basin that overlies the arc basement of Puerto Rico.

Puerto Rico↗

USGS2026

Effects of warming on growth and leaf colonization by litter mat-forming fungi in a wet tropical forest in Puerto Rico

Wet tropical forests are experiencing rising temperatures and increased frequency and intensity of extreme climatic events, such as cyclones, which can increase rates of soil erosion and surface runoff. Fungal litter mats, formed by agaric decomposer fungi, play a crucial role in stabilizing slopes, preventing erosion, and aiding nutrient cycling; however, little is known about how warming affects litter mat growth and function. We investigated two litter mat-forming fungi, Gymnopus johnstonii and Marasmius aff. crinis-equi , in warmed (+4°C above ambient) and control plots in the Luquillo Experimental Forest, Puerto Rico. Growth and time-to-leaf colonization were monitored over 6 weeks in spring (both species) and summer ( G. johnstonii only). We hypothesized that warming would inhibit fungal mat growth and slow leaf colonization, particularly for G. johnstonii since it is drought sensitive. As expected, warming significantly reduced relative growth rates (RGR) in spring, though M. aff. crinis-equi showed slightly higher RGR than G. johnstonii . Leaf colonization was also delayed by 22% in warmed plots, with M. aff. crinis-equi colonizing leaves 4.3 times faster than G. johnstonii . There were significant seasonal differences in response to warming for G. johnstonii , with warming increasing RGR during the consistently wetter summer sampling period. Overall, warming led to significant inhibition of leaf colonization when conditions were dry, whereas there was a trend toward increased colonization in warm and wet conditions. Our findings suggest that warming, combined with drier conditions, is likely to suppress drought-sensitive fungal mat growth, reducing their ability to prevent nutrient and soil loss via erosion.

Luquillo↗

USGS2026

Flood-inundation maps for Río Grande De Loíza in and near Caguas, Puerto Rico, 2026

Digital flood-inundation maps for a 2.7-mile reach of Río Grande De Loíza in Caguas, Puerto Rico, were created by the U.S. Geological Survey. Water-surface profiles were computed for the stream reach by using a one-dimensional, steady-state, step-backwater model. The model was calibrated to the current (2025) stage-streamflow relation (rating curve) for the U.S. Geological Survey streamgage 50055000, Río Grande De Loíza, Puerto Rico. The resulting hydraulic model was then used to compute 16 water-surface profiles for water levels (flood stages) ranging from 19.00 to 34.00 feet at the streamgage; these flood stages range from “moderate flood stage” to above “major flood stage” as defined by the National Weather Service. The 34.00-foot stage exceeds the historical maximum peak stage of 33.20 feet, recorded at the streamgage in 1945. The simulated water-surface profiles were used in combination with a digital elevation model derived from light detection and ranging (lidar) data to map the inundated areas associated with each flood profile. The flood-inundation maps and the supporting hydraulic model produced by this study can be used by emergency managers and local officials to assess flood-mitigation strategies and to define flood-hazard areas to help protect life and property, to coordinate flood-response activities such as evacuations and road closures, and to aid post-flood recovery efforts.

Puerto Rico↗

USGS2026

Erosion and transport of sediments in watersheds of southwest Puerto Rico determined from strontium isotopes and mixing models

Sediment outflow from streams to coastal regions can adversely affect water quality and the health of ecosystems including Puerto Rico's southwestern carbonate platform with shallow marine communities. Knowledge of sediment provenance and erosional areas would aid effective restoration to reduce excessive sedimentation. Strontium (Sr) isotopic ratios ( 87 Sr/ 86 Sr) can be used to differentiate between sediment eroded into streams from different lithologies, or rock types, and used in mixing models to estimate proportions of source materials in downstream mixtures. The goal of this study was to use Sr isotopes to differentiate between river sediments in drainages with distinct geologies and present a case study where Sr isotope mixing models are used to identify sediment sourcing to Río Loco and Bahía de Guánica, a priority watershed of the U.S. Coral Reefs Task Force. For the Río Loco watershed, model results indicate that serpentinites in the uplands contributed approximately 71 % of sediment to Bahía de Guánica. About 21 % of sediment was from trachy-andesitic basalts in the foothills, and the remaining 8 % originated from sources on the coastal plain. Neither the uppermost headwater region of Río Loco, the Canal Principal De Riego Valle De Lajas, nor the tributary Quebrada Susúa appeared to be significant sources of sediment to Río Loco's mainstem or Bahía de Guánica at the time of the study. Our results demonstrate that sediment from the serpentinite-bearing areas of Río Loco watershed was effectively transported to the lowlands, which then discharge into Guánica Bay, and highlight the potential use of Sr isotopes as tracers of sediment provenance. This will, in turn, enable better science-based restoration efforts aimed at reducing terrestrial sediment erosion and impacts to marine ecosystems.

Puerto Rico↗

USGS2025

Trophic assessment of potential competition between invasive cichlids and sport fish in Puerto Rico reservoirs

Objective Several species of New World cichlids have recently invaded reservoirs in Puerto Rico, potentially jeopardizing established recreationally important, albeit nonnative, sport fish populations. Interactions between invasive species and important sport fish must be understood so that they can be mitigated when feasible. This study compared monthly prey consumption between three invasive cichlids (Jaguar Guapote Parachromis managuensis and two species of Amphilophus ) and two principal sport fishes (Largemouth Bass Micropterus nigricans and Butterfly Peacock Bass Cichla ocellaris ) to understand potential for competitive interactions. Methods Stomach contents were evaluated quarterly using gastric lavage for sport fish and destructive sampling techniques for invasive species. Prey proportion by mass was calculated for each individual, and mean prey proportion by mass was calculated for each species by reservoir, month, and overall. Multivariate hypothesis testing using permutational multivariate analysis of variance and visualization of potential diet overlap were performed via analysis of variance and nonmetric multidimensional scaling (NMDS), respectively. Results Fish were the most common prey type, and Threadfin Shad Dorosoma petenense was the most prevalent prey species. Although there were significant differences in diets of invasive cichlids and sport fishes ( P < 0.01), there also was considerable overlap based on low percentages of variation in diets explained by species (3–19% < pseudo- R 2 ) and NMDS visualizations. Diets varied little across reservoirs and months (≤4% pseudo- R 2 ; high overlap in NMDS space). Conclusion These data indicate that the diet of Jaguar Guapote and Amphilophus spp. may overlap with that of sport fish populations in Puerto Rico. Jaguar Guapote relied on Threadfin Shad more than other species, with Amphilophus spp. relying more on detritus and noncrayfish invertebrates. Thus, if prey are limited, Jaguar Guapote may have a more direct effect via competition for prey. Finally, we report evidence of Largemouth Bass recruitment failure and local extirpations that support a hypothesis that invasive cichlids are having negative effects on at least Largemouth Bass, although the mechanism is unclear.

Puerto Rico↗

USGS2025

Comparisons of shoreline positions from satellite-derived and traditional field- and remote-sensing techniques

Satellite-derived shorelines (SDS) have the potential to help researchers answer critical coastal science questions and support work to predict coastal change by filling in the spatial and temporal gaps present in current field-based and remote-sensing data collection methods. The U.S. Geological Survey conducted comparison analyses of traditionally sourced shorelines and SDS in diverse coastal landscapes to determine how SDS could be used in ongoing and future work across varied coastal environments and provided some initial findings that could be used for implementation. Using CoastSeg, a browser-based program for SDS detection and mapping, SDS for the period 1984–2023 for multiple locations across the United States were compared to shoreline positions from traditionally sourced shoreline data. In this report, the authors present these comparisons alongside lessons learned and challenges encountered when building SDS workflows in different coastal locations. Results show that individual SDS have larger uncertainty and yet produced similar linear trends to sparser, traditionally sourced shoreline data; because SDS methods provide orders of magnitude more data than traditional shoreline-detection methods, they can be used to evaluate shoreline behaviors. Refining average scalar slopes used in tidal corrections did not result in substantial decreases in uncertainty. Using lessons from this work to outline needs for regional implementation, initial setup time would be considerable, being on the order of weeks. However, once complete, shoreline detections and analyses are fast (on the order of minutes to hours) and achievable using a desktop computer.

Alaska, Florida, Massachusetts, Washington↗

USGS2025

Fiber-optic sensing for earthquake hazards research, monitoring and early warning

The use of fiber‐optic sensing systems in seismology has exploded in the past decade. Despite an ever‐growing library of ground‐breaking studies, questions remain about the potential of fiber‐optic sensing technologies as tools for advancing if not revolutionizing earthquake‐hazards‐related research, monitoring, and early warning systems. A working group convened to explore these topics; we comprehensively examined the application of fiber optics in various aspects of earthquake hazards, encompassing earthquake source processes, crustal imaging, data archiving, and technological challenges. There is great potential for fiber‐optic systems to advance earthquake monitoring and understanding, but to fully unlock their capabilities requires continued progress in key areas of research and development, including instrument testing and validation, increased dynamic range for applications focused on larger earthquakes, and continued improvement in subsurface and source imaging methods. A key current stumbling block results from the lack of clear data archiving requirements, and we propose an initial strategy that balances data volume requirements with preserving key data for a broad range of future studies. In addition, we demonstrate the potential for fiber‐optic sensing to impact monitoring efforts by documenting the data completeness in a number of long‐term experiments. Finally, we outline the features of a instrument testing facility that would enable progress toward reliable and standardized distributed acoustic sensing data. Overcoming these current obstacles would facilitate progress in fiber‐optic sensing and unlock its potential application to a broad range of earthquake hazard problems.

California, Oregon, Washington↗

USGS2025

Late Pleistocene kinematics of the Great Southern Puerto Rico Fault Zone, Puerto Rico

Several onshore faults in southern Puerto Rico have recently been recognized as Quaternary active. However, the kinematics of these faults, particularly any lateral component, remain largely unconstrained. It is difficult to characterize low strain‐rate faults, partially due to extensive erosional and anthropogenic landscape modification, steep relief, and frequent landsliding, limiting the preservation of geomorphic features that could serve as recorders of fault motion. Here, we constrain the kinematics along sections of the Great Southern Puerto Rico Fault Zone (GSPRFZ) on the southern coastal plain of Puerto Rico. We integrate ∼1‐m‐resolution light detection and ranging (lidar)‐derived topography, historical air photos, and field mapping to identify a series of ∼50–1200‐m‐long fault scarps and lineaments that trend northwest–southeast and extend for ≥25 km across the southern coastal plain. Fault scarps are primarily south facing, cut across topography, and displace Quaternary deposits and landforms. We document multiple offset geomorphic markers, including channel thalwegs and interfluves formed in deposits previously mapped as Quaternary piedmont alluvial plain. We observe both vertical (south‐side‐down) and right‐lateral meter‐scale displacements, which indicate that the GSPRFZ accommodates right‐lateral oblique motion in the late Pleistocene, consistent with northeast motion of the Puerto Rico and the Virgin Islands microplate away from the Hispaniola block.

Puerto Rico↗

USGS2025

Warming induces unexpectedly high soil respiration in a wet tropical forest

Tropical forests are a dominant regulator of the global carbon cycle, exchanging more carbon dioxide with the atmosphere than any other terrestrial biome. Climate models predict unprecedented climatic warming in tropical regions in the coming decades; however, in situ field warming studies are severely lacking in tropical forests. Here we present results from an in situ warming experiment in Puerto Rico, where soil respiration responses to + 4 o C warming were assessed half-hourly for a year. Soil respiration rates were 42-204% higher in warmed relative to ambient plots, representing some of the highest soil respiration rates reported for any terrestrial ecosystem. While respiration rates were significantly higher in the warmed plots, the temperature sensitivity (Q 10 ) was 71.7% lower, pointing to a mechanistic shift. Even with reduced Q 10 , if observed soil respiration rates persist in a warmer world, the feedback to future climate could be considerably greater than previously predicted or observed.

Puerto Rico↗

USGS2025

Complex sound scattering layer and water-column dynamics over a mesophotic coral ecosystem: Southwest Puerto Rico, U.S.A.

A nearly 5-month record of high-resolution temperature and acoustic backscatter profiles from the upper insular slope off southwest Puerto Rico reveals complex sound scattering layer (SSL) dynamics over a mesophotic coral ecosystem (MCE). The SSLs exhibited both diel and reverse diel vertical migration, thin layer (< 5 m) and multiple layer formations, depth modulation due to internal waves, and vertical layering in the absence of water column stratification. The long-term observations also capture SSL and water column dynamics across changing seasons and two category five hurricanes, Irma and María. The SSLs, likely comprosed of zooplankton, represent an important food source for both the sessile ( e.g. , corals and sponges) and mobile ( e.g. , fish) MCE taxa, and their effective vertical mobility underscores their importance to trophic connectivity between the upper and lower slope MCEs, as well as the shelf. Our results also underscore the challenges in adequately resolving zooplankton aggregations using conventional sampling techniques.

Puerto Rico↗

USGS2025

Slow rupture, long rise times, and multi-fault geometry: The 2020 M6.4 southwestern Puerto Rico mainshock

The M 6.4 mainshock of the southwestern Puerto Rico seismic sequence on 7 January 2020, was one of the most impactful modern earthquakes in the northeastern Caribbean. Due to its offshore location and complex aftershock distribution, its source kinematics remain poorly constrained. This active sequence illuminated a complex set of previously unrecognized structures that indicate multiple causative faults may have slipped during its rupture. Here, we utilize seismic and geodetic observations to enhance model resolution, estimate the finite slip of the mainshock, and test a multi-segment, geologically realistic fault geometry. Our refined model finds a lower rupture velocity and longer rise times than typical for an event of this magnitude. This indicates a slow-evolving rupture process that resembles characteristics of a tsunami earthquake. Although this normal/strike-slip faulting event was not tsunamigenic, these qualities, if pervasive for this region, have important implications for future seismic monitoring and hazards in southwestern Puerto Rico.

Puerto Rico↗

USGS2025

Topographic controls on landslide mobility: Modeling hurricane-induced landslide runout and debris-flow inundation in Puerto Rico

In 2017, Hurricane Maria triggered more than 70 000 landslides in Puerto Rico. After initiation, these predominantly shallow landslides were mobilized to varying extents – some landslides only traveled partway downslope, whereas others reached drainage channels and were mobilized into long-traveled debris flows that could severely impact roads and infrastructure. Thus, forecasting potential landslide runout and inundation zones is critical for estimating landslide and debris-flow hazards. Here we conduct an in-depth topographic analysis of landslide-affected areas from nine study areas and apply a linked modeling technique to estimate locations susceptible to varying degrees of landslide runout in the Lares, Utuado, and Naranjito municipalities. We find that the longest runout lengths are observed on high-relief escarpments, although highly mobile long-runout debris flows also occurred in lower-relief dissected uplands. These topographic differences indicate that landslides that are initiated under similar conditions and possess equal potential to be mobilized as debris flows may not travel the same distances or affect the same areal extent. Our modeling approach allows the local topography to automatically control the implementation of two runout methods: (1) H/L runout zones are assigned directly downslope of landslide source zones, and (2) debris-flow inundation zones are estimated in the presence of a channel network. Debris-flow volumes are calculated as a function of area-integrated growth factors, estimated as a function of the upstream areas susceptible to shallow landslides. Applying our empirical modeling scheme over an area of 560 km 2 , our results highlight the efficacy of our methods for the assessment of the potential for landslide runout and debris-flow inundation over diverse terrains with varied susceptibility.

Puerto Rico↗

USGS2025

Smectite-rich horizons in Inceptisols trigger shallow landslides in tropical granitic terranes

Puerto Rico was affected by >70,000 landslides in the wake of 2017 Hurricane Maria, and landslide prevalence was especially high in the Utuado region in the Cordillera Central. Landslide density was highest where soil parent material is granodiorite; landslide slip surfaces tended to be shallow (<60 cm), and often were mobilized rapidly and with long runout distances. This study combines field observations with soil mineralogy (bulk and clay fractions), soil geochemistry (bulk fraction), and soil strength as determined by field cone penetrometer testing (CPT) to test the hypothesis that clay-rich subsoil horizons function as slip planes when water-logged. Soil pits were excavated to depths of ∼200 cm in Ultisols on an undulating plateau and to ∼100 cm in Inceptisols on steep slopes (36-43 o ) that flank the plateau and cone penetrometer tests (CPT) were done within 2 m of the soil pit. Six pits were located adjacent to scarps from previous landslides, enabling analysis of soil profiles downward through extrapolated slip surfaces. Results from X-ray diffraction (XRD), X-ray fluorescence (XRF) and thermogravimetric analysis (TGA) indicate that soils are heterogeneous, often with subsoil horizons enriched in clay minerals and immobile elements (Al, Fe, Ti). Inceptisols on steep slopes often contain smectite-rich horizons at 30–60 cm depth that appear to function as slip surfaces; in other Inceptisols, such horizons are not present and landslide susceptibility is potentially lower. In Ultisols, soil mineralogy is dominated by kaolinite with minor halloysite, and increased kaolinite content at soil depths ≥80 cm at some sites suggests potential slip surfaces enhancing probability of landslides. The origin of clay-rich horizons appears to be (1) fractures in granodiorite that facilitate water flow and leaching, accelerating mineral dissolution during early weathering stages, and (2) smectite-rich buried soils under permeable colluvium likely deposited by a prior mass wasting event. Where clay-rich layers occur beneath more-permeable horizons, rapid infiltration then absorption of water in clay-rich subsoil horizons causes decreased shear strength and increased landslide susceptibility.

Puerto Rico↗

USGS2025

Hybrid coral reef restoration can be a cost-effective nature-based solution to provide protection to vulnerable coastal populations

Coral reefs can mitigate flood damages by providing protection to tropical coastal communities whose populations are dense, growing fast, and have predominantly lower-middle income. This study provides the first fine-scale, regionally modeled valuations of how flood risk reductions associated with hybrid coral reef restoration could benefit people, property, and economic activity along Florida and Puerto Rico’s 1005 kilometers of reef-lined coasts. Restoration of up to 20% of the regions’ coral reefs could provide flood reduction benefits greater than costs. Reef habitats with the greatest benefits are shallow, nearshore, and fronting low-lying, vulnerable communities, which are often where reef impacts and loss are the greatest. Minorities, children, the elderly, and those below the poverty line could receive more than double the hazard risk reduction benefits of the overall population, demonstrating that reef restoration as a nature-based solution can have positive returns on investment economically and socially by providing protection to the most vulnerable people.

Florida↗