Geology Reports⌕ Search

Geology topics

Research about Puerto Rico, U.S. Virgin Islands

Source-linked reports with geographic coverage including Puerto Rico, U.S. Virgin Islands.

At least 19 recordsLinked to original sources

2025 USGS National Seismic Hazard Model for Puerto Rico and the U.S. Virgin Islands: Overview of model and hazard results

The U.S. Geological Survey recently updated the National Seismic Hazard Model (NSHM) for Puerto Rico and the U.S. Virgin Islands (PRVI). The first version of the PRVI NSHM was released in 2003, and therefore this 2025 update includes over 20 years of new geologic, geophysical, and engineering data, methods, and models. Updates follow similar efforts performed in the recent 2023 50-state NSHM. However, this is the first NSHM in which we: (1) apply an inversion methodology to subduction interface fault sources in the earthquake rupture forecast (ERF) model; (2) develop scaled backbone median ground-motion models and independent aleatory variability models that are applied in the ground-motion characterization (GMC) model; and (3) calculate epistemic uncertainty related to alternative scenarios in the ERF and GMC models for all grid points in the study region. Long-term time-independent mean hazard calculations were performed for peak ground acceleration and 5%-damped pseudospectral acceleration at 21 spectral periods from 0.01- to 10.0-s, for eight National Earthquake Hazards Reduction Program site conditions ranging from V S30 = 150 to 1500 m/s, and for 2%, 5%, and 10% in 50-year probabilities of exceedance (return periods of 2475, 975, and 475 years, respectively). Epistemic uncertainty, in the form of selected percentiles, is also provided for a suite of test sites and all grid points in the study region for limited periods, site conditions, and probabilities of exceedance. Selected results, including comparisons with the 2003 PRVI NSHM, are shown and discussed for selected periods, site conditions, and probabilities of exceedance. When comparing the 2025 PRVI NSHM with the 2003 PRVI NSHM, hazard is generally higher at shorter periods and lower at longer periods, as a result of updates in both ERF and GMC models. The 2025 PRVI NSHM is applicable for return periods greater than ∼475 or less than ∼10,000 years.

Puerto Rico, U.S. Virgin Islands↗

Ground-motion characterization for the 2025 U.S. National Seismic Hazard Model for Puerto Rico and the U.S. Virgin Islands

We develop the ground-motion characterization (GMC) for the 2025 U.S. National Seismic Hazard Model for Puerto Rico and the U.S. Virgin Islands (NSHM-PRVI) for earthquakes in active crustal, subduction interface, and subduction intraslab regimes. Using ground-motion models (GMMs) from the Next-Generation Attenuation (NGA)-West2 and NGA-Subduction projects, the GMC is parameterized by scaled-backbone models for median ground motions and by independent logic trees of aleatory variability. We introduce several novel GMC features into the U.S. NSHM: (1) use of regional ground-motion data for modeling median ground motions; (2) development of scaled-backbone models for median ground motions; and (3) development of independent logic trees of aleatory variability from variance components of GMMs and computed from multiple ground-motion datasets, and incorporating regional ground-motion data effects on variability. We compute probabilistic seismic hazard curves and maps to evaluate the contributions from the GMC components and for comparison with the 2003 GMC. Contributions of the GMC to epistemic uncertainty in seismic hazard are evaluated through spatial variations in epistemic uncertainty in hazard maps, comparisons between mean hazard curves and fractiles, as well as investigations of the contributions of logic-tree branches to hazard maps and curves. Comparisons between seismic hazard from the 2025 and 2003 GMCs allow examination of the changes in hazard curves and mapped values with 2% and 10% probabilities of exceedance in 50 years. The 2025 GMC exhibits modest changes in median predictions, relative to the 2003 GMC; overall, values of aleatory variability are higher, except at long periods ( ), resulting in changes to probabilistic ground motions at low probabilities of exceedance (<10% probability of exceedance in 50 years). Changes in hazard at 2% and 10% probabilities of exceedance in 50 years are also relatively modest (within 20%) at most sites, with the impacts from the 2025 GMC exhibiting minor reductions and increases at 0.2 and 1.0 s periods, respectively, relative to the 2003 GMC.

Puerto Rico, U.S. Virgin Islands↗

Suitability of ground-motion models for seismic hazard assessment in Puerto Rico and the U.S. Virgin Islands

We perform linear mixed-effects analyses with a ground-motion dataset to evaluate how well ground-motion models (GMMs) fit active crustal, subduction interface, and subduction intraslab earthquakes in Puerto Rico and the US Virgin Islands. Most of the GMMs reproduce the spatial variation in peak ground motions with earthquake magnitude and rupture distance but predict ground motions 0.3–1.0 natural log units (35%–270%) greater than observations. Two GMMs developed for Puerto Rico that are based on ground-motion records from mostly magnitude 4–5 earthquakes do not perform as well. We attribute the period-dependent overprediction in observed ground motions to differences between the observed site response and the linear site response in the GMMs. Consequently, we developed region-specific GMMs by adjusting the period-dependent linear site response coefficients and period-dependent constant coefficients to remove most of the bias between predicted ground motions and observations. For the analysis, we compile ground-motion records and process waveforms to build a dataset with 10,127 records at 72 stations from 849 magnitude 4.0–6.4 earthquakes between 1 January 2006 and 31 March 2024. The earthquakes include active crustal, subduction interface, and subduction intraslab events. We evaluate the GMMs using the time-averaged shear wave speed in the top 30 m ( ), which we compile from site surveys and proxy values computed from horizontal to vertical spectral ratios. Site terms exhibit strong consistency across GMMs and crustal and subduction earthquakes, indicating that the linear mixed-effects analysis successfully isolates the effects of local site response. The event terms show little spatial correlation and more substantial variability than in other regions, which we attribute to uncertainties in the earthquake magnitudes. This analysis guides the selection of the GMMs for the 2025 update of the National Seismic Hazard Model for Puerto Rico and the US Virgin Islands.

Puerto Rico, U.S. Virgin Islands↗

The U.S. Geological Survey 2025 Puerto Rico and U.S. Virgin Islands time-independent earthquake rupture forecast

We present the 2025 U.S. Geological Survey Puerto Rico and U.S. Virgin Islands (PRVI) time‐independent earthquake rupture forecast (ERF), developed for the 2025 update to the National Seismic Hazard Model (NSHM) for PRVI. The updated ERF improves upon a prior model from 2003, including an expanded fault inventory with slip‐rate estimates, updated seismicity catalogs, and refined subduction zone geometries and deformation models. It applies the fault‐system inversion methodology to solve for rates of ruptures on modeled faults, adapted from the 2023 NSHM (NSHM23) for the western United States, including the first application of the inversion to model rates on a U.S. subduction interface. Off‐fault and intraslab seismicity are constrained by observed seismicity and use updated methods developed for NSHM23. Uncertainties in model components are substantial, and the ERF represents epistemic uncertainties through a comprehensive logic tree consisting of 1.7 billion logic‐tree branches combined across all sources.

Puerto Rico, U.S. Virgin Islands↗

Ground-motion aleatory-variability models for Puerto Rico and the U.S. Virgin Islands

I develop independent logic trees for aleatory variability for crustal and subduction-zone (interface and intraslab) earthquakes for seismic hazards analyses in Puerto Rico and the U.S. Virgin Islands (PRVI) from existing suites of ground-motion models (GMMs) and from ground-motion datasets, including a regional PRVI dataset. The aleatory variability models are parameterized using a partially nonergodic partitioning of standard deviation that consists of independently developed between-event ( ), site-to-site ( ), and event-corrected single-station ( ) standard deviation components. The effects of nonlinear site response on aleatory variability are incorporated through additional terms that modify the standard deviation components. Because one goal of this work is to develop independent logic trees for aleatory variability that synthesize the aleatory variability models from GMMs, I make use of the functional forms of the input GMMs. The PRVI dataset contains a limited number of stations with high-quality site metadata and does not contain records from earthquakes with magnitudes greater than 6.1, so I choose not to develop the aleatory variability models from the regional dataset alone. Instead, the standard deviation components from regional ground-motion data are evaluated against the components derived from GMMs and from available global datasets, and regionalized standard deviation components are incorporated where there is evidence that regional effects exhibit substantial differences. The resulting logic trees for aleatory variability consist of models of and that are consistent with semiempirical GMMs for active crustal and subduction-zone regimes, and two alternative models of , including one model that exhibits site-to-site variability informed by PRVI data, with values that exceed global models. The aleatory variability models may be considered in future hazards assessments in PRVI to simplify the hazard calculations, to incorporate regional ground-motion variability effects, and to enable direct logic-tree weighs of aleatory variability.

Puerto Rico, U.S. Virgin Islands↗

Using gridded seismicity to forecast the long-term spatial distribution of earthquakes for the 2025 Puerto Rico and U.S. Virgin Islands National Seismic Hazard Model

Gridded (or background) seismicity models are a critical component of probabilistic seismic hazard assessments, accounting for off‐fault and smaller‐magnitude earthquakes. They are typically developed by declustering and spatially smoothing an earthquake catalog to estimate a long‐term seismicity rate that can be used to forecast future earthquakes. Here, we present new gridded seismicity models for use in the 2025 National Seismic Hazard Model (NSHM) for Puerto Rico and the U.S. Virgin Islands (PRVI). The previous PRVI NSHM was released in 2003, and our new models incorporate updates to both data and methodology. We utilize an updated earthquake catalog based on improved Puerto Rico Seismic Network data with newly characterized completeness epochs. The catalog is divided into crustal, subduction interface, and intraslab seismicity using new methods and Slab2 subduction zone geometries. To forecast the long‐term spatial distribution of earthquakes, we use an updated methodology developed for the 2023 U.S. 50‐state NSHM, considering three declustering methods and two spatial smoothing methods based on 2D Gaussian kernels. To adapt it for the complex seismotectonics of the region, we also adopt probabilistic methods to account for events with unknown depths and uncertainties in tectonic classification, and develop a new method for spatial scaling to counteract the effects of spatial variability in network coverage while maintaining the use of smaller events. Finally, we test the performance of these spatial models in forecasting the location of M w ≥ 5earthquakes in the region. Our updated methodology improves the representation of epistemic uncertainty relative to the 2003 model, and our results demonstrate the effectiveness of the new measures we have introduced to address heterogeneities in network detection and systematically evaluate forecast performance.

Puerto Rico, U.S. Virgin Islands↗

Subduction zone earthquake catalog separation tool: Implementation in the USGS 2025 Puerto Rico and U.S. Virgin Islands National Seismic Hazard Model

The U.S. Geological Survey (USGS) periodically releases updates to National Seismic Hazard Model (NSHM) for the United States and its territories leveraging current scientific knowledge and methodologies to guide public policy, building codes, and risk assessments regarding potential ground shaking due to earthquakes that may result in infrastructure damage. In subduction zones, there is a need to separate the earthquake catalog into tectonic regimes to create specific seismicity models for which the most appropriate ground‐motion models are then applied. Here, we describe newly developed methods and software, called CatSep, that classifies subduction zone events into three primary tectonic regimes: crustal, interface, and intraslab. This method incorporates information about the location of the earthquake relative to the subducting slab, the depth of the Mohorovičić discontinuity, and the earthquake’s moment tensor. Applying this method is a first step in the NSHM workflow for regions covering U.S. subduction zones. Results using this subduction zone earthquake catalog separation tool for the 2025 Puerto Rico and U.S. Virgin Islands NSHM earthquake catalog are presented and analyzed.

Puerto Rico, U.S. Virgin Islands↗

Geologic input databases for the 2025 Puerto Rico – U.S. Virgin Islands National Seismic Hazard Model update: Crustal faults component

The last National Seismic Hazard Model (NSHM) for Puerto Rico and the U.S. Virgin Islands (PRVI) was published in 2003. In advance of the 2025 PRVI NSHM update, we created three geologic input databases to summarize new onshore and offshore fault source information in the northern Caribbean region between 62°–70° W and 16°–21° N. These databases, of fault sections, fault‐zone polygons, and geologic estimates of fault activity (fault‐slip rate and earthquake recurrence intervals) at specific sites, document updates to fault parameters used in prior seismic hazard models in PRVI. Fault sources were reviewed from published studies since 2003, which document substantial changes to the understanding of fault location, geometry, or activity. New fault section sources were added for features that meet the criteria of (1) length ≥7 km, (2) unequivocal evidence of recurrent tectonic Quaternary activity, and (3) documentation that is publicly available in a peer‐reviewed source. In addition, we revised several broad areal sources, such as the Mona and Anegada extensional zones. The 2003 model included three fault sections and two fault‐zone polygons (areal sources). These databases include 35 fault sections, 6 fault‐zone polygons, and 51 earthquake geology sites. To characterize fault activity rates, slip‐rate bins were assigned based on landscape expression and paleoseismic trench observations for faults without published slip‐rate sites. Additional fault sources were evaluated but not included in these databases due to a lack of published information about fault location, geometry, or recurrent Quaternary activity. The PRVI NSHM 2025 geologic input databases describe crustal faulting; the geometries and coupling of Puerto Rico subduction zone and Muertos Trough models are considered in a separate database. Updates to the fault sections, fault‐zone polygons, and earthquake geology databases can help inform the location and recurrence rate of damaging earthquakes in the PRVI NSHM implementation.

Puerto Rico, U.S. Virgin Islands↗

A biological condition gradient for coral reefs in the US Caribbean Territories: Part I. Coral narrative rules

As coral reef condition and sustainability continue to decline worldwide, losses of critical habitat and their ecosystem services have generated an urgency to understand and communicate reef response to management actions, environmental contamination, and natural disasters. Increasingly, coral reef protection and restoration programs emphasize the need for robust assessment tools for protecting high-quality waters and establishing conservation goals. Of equal importance is the need to communicate assessment results to stakeholders, beneficiaries, and the public so that environmental consequences of decisions are understood. The Biological Condition (BCG) model provides a structure to evaluate the condition of a coral reef in increments of change along a gradient of human disturbance. Communication of incremental change, regardless of direction, is important for decision makers and the public to better understand what is gained or lost depending on what actions are taken. We developed a narrative (qualitative) Biological Condition Gradient (BCG) from the consensus of a diverse expert panel to provide a framework for coral reefs in US Caribbean Territories. The model uses narrative descriptions of biological attributes for benthic organisms to evaluate reefs relative to undisturbed or minimally disturbed conditions. Using expert elicitation, narrative decision rules were proposed and deliberated to discriminate among six levels of change along a gradient of increasing anthropogenic stress. Narrative rules for each of the BCG levels are presented to facilitate the evaluation of benthic communities in coral reefs and provide specific narrative features to detect changes in coral reef condition and biological integrity. The BCG model can be used in the absence of numeric, or quantitative metrics, to evaluate actions that may encroach on coral reef ecosystems, manage endangered species habitat, and develop and implement management plans for marine protected areas, watersheds, and coastal zones. The narrative BCG model is a defensible model and communication tool that translates scientific results so the nontechnical person can understand and support both regulatory and non-regulatory water quality and natural resource programs.

Puerto Rico, U.S. Virgin Islands↗

Monitoring storm tide, flooding, and precipitation from Hurricane Maria in Puerto Rico and the U.S. Virgin Islands, September 2017

Hurricane Maria made landfall south of Yabucoa Harbor, Puerto Rico, as a category 4 hurricane with maximum sustained winds of 155 miles per hour on September 20, 2017. The hurricane devastated much of the U.S. Virgin Islands and Puerto Rico. The U.S. Geological Survey (USGS), in cooperation with Federal Emergency Management Agency, installed a temporary monitoring network of 13 water-level and barometric pressure sensors along the coast of Puerto Rico prior to the storm. In addition to the temporary sensors, the USGS maintains 99 permanent real-time streamgages and 36 real-time precipitation stations in Puerto Rico. The real-time data, updated hourly, during and after the hurricane are displayed in the USGS Flood Event Viewer ( https://stn.wim.usgs.gov/FEV/#MariaSeptember2017 ) and in the USGS National Water Information System. The USGS measured 181 coastal and riverine high-water marks throughout Puerto Rico after the storm. Water elevations are referenced to the Puerto Rico Vertical Datum of 2002 (PRVD02) and local datums in Puerto Rico and to the U.S. Virgin Islands Vertical Datum of 2009 (VIVD09) in the U.S. Virgin Islands. Data from the Hurricane Maria storm-tide network are available in tab-delimited, American Standard Code for Information Interchange (ASCII) format and Network Common Data Form (NetCDF) format by site for each sensor through the USGS Flood Event Viewer.

Puerto Rico, U.S. Virgin Islands↗

Monitoring storm tide and flooding from Hurricane Irma along the U.S. Virgin Islands, Puerto Rico, and the Southeastern United States, September 2017

Hurricane Irma skirted the northern coasts of the U.S. Virgin Islands and Puerto Rico, with maximum sustained winds of 185 miles per hour (mi/h) on September 6, 2017. The hurricane first made landfall in Florida near Cudjoe Key, in the lower Florida Keys, with maximum sustained winds of 130 mi/h on September 10, 2017. The hurricane made a second Florida landfall on Marco Island, Florida, with maximum sustained winds of 115 mi/h on September 10, 2017. The U.S. Geological Survey (USGS), in cooperation with Federal Emergency Management Agency, deployed a temporary monitoring network of water-level and barometric pressure sensors at 249 locations along the Puerto Rico, Florida, Georgia, and South Carolina coasts to record the timing, areal extent, and magnitude of hurricane storm tide and coastal flooding generated by the hurricane. Immediately following the passage of Hurricane Irma, the sensors were retrieved, and the data were disseminated on the USGS Flood Event Viewer ( https://stn.wim.usgs.gov/FEV/#IrmaSeptember2017 ). The storm-tide peak data values were verified by comparing data from hydrologic recorders and nearby high-water marks (HWMs). Following the hurricane, 508 independent HWM locations were flagged and surveyed relative to the North American Vertical Datum of 1988, National Geodetic Vertical Datum of 1929, or a local datum along the southeastern U.S. coast, and to Puerto Rico Vertical Datum of 2002 in Puerto Rico. Most HWMs were in Florida because of the path of the hurricane. The data from the Hurricane Irma storm-tide network are available on a provisional basis in tab-delimited, American Standard Code for Information Interchange (ASCII) format and Network Common Data Form (NetCDF) format by site for each sensor by using the USGS Flood Event Viewer.

Puerto Rico, U.S. Virgin Islands↗

Drought in the U.S. Caribbean: Impacts to freshwater ecosystems

Healthy and functioning freshwater ecosystems are needed for successful conservation and management of native fish and invertebrate species, and the services they provide to human communities, across the U.S. Caribbean. Yet streams, rivers, and reservoirs are vulnerable to the effects of extreme weather events, urbanization, energy and water development, and other environmental and human-caused disturbances (Neal et al., 2009). One major management concern is the impact of prolonged drought on freshwater ecosystems. Drought impacts streamflow, dissolved oxygen content, water quality, stream connectivity, available habitat, and other important freshwater habitat characteristics necessary for sustaining fish and invertebrate populations (Covich et al., 2006). These changes can impact species interactions, abundance, life history events, and the presence of native and non-native species (Larsen, 2000; Covich et al., 2006; Ramírez et al., 2018). Drought impacts aquatic ecosystems and species both in the short-term and long-term, depending on the severity and duration of the event (e.g. Covich et al., 2006 ). In Puerto Rico, all native freshwater fish, shrimp, and snail species spend part of their lives in estuarine and marine ecosystems and depend on being able to move between these habitats to survive, so maintaining connectivity is key (e.g., Engman et al., 2017 ). Freshwater ecosystems also provide recreational, cultural, and ecological value to humans ( Kwak et al., 2007 ; Neal et al., 2009 ). For example, some communities in Puerto Rico engage in artisanal shrimp and freshwater crab fishing ( Neal et al., 2009 ). Artisanal fishing for postlarvae gobioids, known colloquially as “cetí” also occurs at the river mouths of large drainages and has strong cultural significance in parts of Puerto Rico, such as Arecibo (Kwak et al., 2016). The U.S. Virgin Islands (USVI) is particularly sensitive to drought, because almost all streams are ephemeral and typically only flow after rainfall. These intermittent channels, known locally as “ghuts”, run down the surface of steep slopes, rather than through the ground, and are important sources of freshwater. Natural springs are often located in ghuts and can form pools of freshwater that serve as habitat for wetland and migratory birds, freshwater shrimp and fish, and amphibians ( Nemeth and Platenburg, 2007 ; Gardner, 2008 ).

Puerto Rico, U.S. Virgin Islands↗

Earthquakes and tsunamis in Puerto Rico and the U.S. Virgin Islands

Many earthquakes and tsunamis have occurred in the northeastern Caribbean, where the movements of the Earth's surface plates are rapid and complicated. Future such events pose serious hazards to the 3.7 million people who live in Puerto Rico and the U.S. Virgin Islands.

Puerto Rico, U.S. Virgin Islands↗

United States National seismograph network

The USGS National Earthquake Information Center (NEIC) has planned and is developing a broadband digital seismograph network for the United States. The network will consist of approximately 150 seismograph stations distributed across the contiguous 48 states and across Alaska, Hawaii, Puerto Rico and the Virgin Islands. Data transmission will be via two-way satellite telemetry from the network sites to a central recording facility at the NEIC in Golden, Colorado. The design goal for the network is the on-scale recording by at least five well-distributed stations of any seismic event of magnitude 2.5 or greater in all areas of the United States except possibly part of Alaska. All event data from the network will be distributed to the scientific community on compact disc with read-only memory (CD-ROM).

Puerto Rico, U.S. Virgin Islands↗

Water resources data, Puerto Rico and the U.S. Virgin Islands, water year 1987

Water-resources data for surface-water, quality-of-water, and ground-water records for the 1987 water year for Puerto Rico and the U.S Virgin Islands, consists of records of discharge, water quality of streams, and water levels of wells. This report contains discharge records for 49 streamflow-gaging stations, and 1 crest-stage, partial-record streamflow stations; water quality records-for 16 streamflow-gaging stations, 42 ungaged streamsites, 11 lake sites, 1 lagoon, and 1 bay; and water-level records for 56 observation wells. These data represent that part of the National Water Data System collected by the U.S. Geological Survey and cooperating local and federal agencies in Puerto Rico and the U.S. Virgin Islands.

Puerto Rico, U.S. Virgin Islands↗