Search USGS⌕ Search

Geology topics

Research about Puerto Rico

Source-linked reports with geographic coverage including Puerto Rico.

At least 55 records · Page 3Linked to original sources

Understory plant communities show resistance to drought, hurricanes, and experimental warming in a wet tropical forest

Global climate change has led to rising temperatures and to more frequent and intense climatic events, such as storms and droughts. Changes in climate and disturbance regimes can have non-additive effects on plant communities and result in complicated legacies we have yet to understand. This is especially true for tropical forests, which play a significant role in regulating global climate. We used understory vegetation data from the Tropical Responses to Altered Climate Experiment (TRACE) in Puerto Rico to evaluate how plant communities responded to climate warming and disturbance. The TRACE understory vegetation was exposed to a severe drought (2015), 2 years of experimental warming (4°C above ambient in half of the plots, 2016–2017 and 2018–2019), and two major hurricanes (Irma and María, September 2017). Woody seedlings and saplings were censused yearly from 2015 to 2019, with an additional census in 2015 after the drought ended. We evaluated disturbance-driven changes in species richness, diversity, and composition across ontogeny. We then used Bayesian predictive trait modeling to assess how species responded to disturbance and how this might influence the functional structure of the plant community. Our results show decreased seedling richness after hurricane disturbance, as well as increased sapling richness and diversity after warming. We found a shift in species composition through time for both seedlings and saplings, yet the individual effects of each disturbance were not significant. At both ontogenetic stages, we observed about twice as many species responding to experimental warming as those responding to drought and hurricanes. Predicted changes in functional structure point to disturbance-driven functional shifts toward a mixture of fast-growing and drought-tolerant species. Our findings demonstrate that the tropical forest understory community is more resistant to climatic stressors than expected, especially at the sapling stage. However, early signs of changes in species composition suggest that, in a warming climate with frequent droughts and hurricanes, plant communities might shift over time toward fast-growing or drought-tolerant species.

Bosque experimental de Luquillo, Luquillo Experime↗

U.S. Geological Survey response to Hurricane Maria flooding in Puerto Rico and characterization of peak streamflows observed September 20–22, 2017

Hurricane Maria struck the island of Puerto Rico on September 20, 2017, as a Category 4 storm. The hurricane traversed the island from southeast to northwest and produced recorded 48-hour rainfall totals of up to 30.01 inches. Estimates of the human death toll range from 2,975 to 4,645, possibly more. The U.S. Geological Survey (USGS) hydrologic monitoring network sustained substantial wind and flood damage during the hurricane. Eighty-five of the 300 hydrologic monitoring stations operating in Puerto Rico and the U.S. Virgin Islands prior to the passage of Hurricane Maria were destroyed or damaged. During the weeks and months after the hurricane, USGS field crews in Puerto Rico prioritized repair of the hydrologic monitoring network and collected hydrologic information to characterize the magnitude of observed peak streamflows at 20 streamgage and to develop new theoretical stage-streamflow relations for 58 streamgages where stream channels were substantially altered; the theoretical stage-streamflow relations were used to estimate Hurricane Maria peak streamflows for 39 of those sites. As part of a pilot program, USGS field crews installed continuous slope-area monitoring equipment at two remote streamgages to automate the collection of high-streamflow stage data. Hurricane Maria peak streamflows and rankings were determined for 73 USGS streamgages in Puerto Rico. New rank 1 period-of-record peak streamflows occurred at 28 sites, rank 2 period-of-record peak streamflows occurred at 17 sites, and rank 3 period-of-record peak streamflows occurred at 9 sites; period-of-record peak streamflows at the remaining 19 sites either ranked from 4th to 20th or were not ranked. Annual exceedance probabilities for 53 unregulated peak streamflows ranged from greater than 50.0 percent (recurrence interval of less than 2 years) to 0.3 percent (recurrence interval of 333 years), with the majority (28 of 53) in the range of 10.0–2.1 percent (recurrence intervals of 10–48 years). A comparison of period-of-record ranks for the largest flood events that have occurred in Puerto Rico since the 1960s indicated that Hurricane Maria produced more record peak streamflows than either Hurricane Hortense in 1996 or Hurricane Georges in 1998. Limited pre-1960s hydrologic data preclude quantitative comparison with earlier storms. As part of this study, a maximum peak-streamflow envelope curve for Puerto Rico was developed using historical peak-streamflow information available through 2017. Other post-Hurricane Maria USGS activities summarized in this report include (1) Global Navigation Satellite System surveys at all stations in the USGS hydrologic monitoring network, used to tie the network to the Puerto Rico Vertical Datum of 2002; and (2) telemetered monitoring of the Lago Guajataca Dam in northwestern Puerto Rico, which was damaged and at risk of failure from October to December 2017.

Puerto Rico↗

Modeling and estimating co-occurrence between the invasive Shiny Cowbird and its Puerto Rican hosts

Invasive species threaten island biodiversity globally. For example, the Shiny Cowbird ( Molothrus bonariensis ) parasitizes many of Puerto Rico’s endemic species, particularly in the open forests in the island’s southwest. Less is known, however, about cowbird parasitism in the agro-ecological highlands, which contain a patchwork of forests, shaded-coffee plantations, and coffee farms without shade. In this paper, we estimated co-occurrence rates, a potential indicator of parasitism rates, between the cowbird and four host species across these three land uses, hypothesizing that cowbirds would most likely co-occur with their hosts in shaded-coffee farms. We also hypothesized that the presence of host species would increase the probability of cowbird occurrence. To investigate these hypotheses, we developed three Bayesian hierarchical occupancy models: one where the hosts and parasite occurred independently, one that used the latent host species richness as a predictor of cowbird occurrence, and one that used each latent host occurrence state as predictors. These methods addressed observation errors and appropriately propagated error to our predictions of co-occurrence rates. We selected the best performing model using WAIC, then used it to predict co-occurrence rates. While there was some evidence that host species richness increased the probability of cowbirds, the parsimonious model assumed no interaction. With this model, we found that cowbirds were more likely to overlap with certain hosts in shaded-coffee plantations. This may suggest increased parasitism at these plantations, potentially presenting challenges for managers who advocate for shade restoration to gain ecological services such as biodiversity conservation.

Biological Invasions↗

Extreme rainstorms drive exceptional organic carbon export from forested humid-tropical rivers in Puerto Rico

Extreme rainfall events in the humid-tropical Luquillo Mountains, Puerto Rico export the bulk of suspended sediment and particulate organic carbon. Using 25 years of river carbon and suspended sediment data, which targeted hurricanes and other large rainstorms, we estimated biogenic particulate organic carbon yields of 65 ± 16 tC km −2 yr −1 for the Icacos and 17.7 ± 5.1 tC km −2 yr −1 for the Mameyes rivers. These granitic and volcaniclastic catchments function as substantial atmospheric carbon-dioxide sinks, largely through export of river biogenic particulate organic carbon during extreme rainstorms. Compared to other regions, these high biogenic particulate organic carbon yields are accompanied by lower suspended sediment yields. Accordingly, particulate organic carbon export from these catchments is underpredicted by previous yield relationships, which are derived mainly from catchments with easily erodible sedimentary rocks. Therefore, rivers that drain petrogenic-carbon-poor bedrock require separate accounting to estimate their contributions to the geological carbon cycle.

Puerto Rico↗

Mature diffuse tectonic block boundary revealed by the 2020 southwestern Puerto Rico seismic sequence

Distributed faulting typically tends to coalesce into one or a few faults with repeated deformation. The progression of clustered medium-sized (≥Mw4.5) earthquakes during the 2020 seismic sequence in southwestern Puerto Rico (SWPR), modeling shoreline subsidence from InSAR, and sub-seafloor mapping by high-resolution seismic reflection profiles, suggest that the 2020 SWPR seismic sequence was distributed across several short intersecting strike-slip and normal faults beneath the insular shelf and upper slope of Guayanilla submarine canyon. Multibeam bathymetry map of the seafloor shows significant erosion and retreat of the shelf edge in the area of seismic activity as well as slope-parallel lineaments and submarine canyon meanders that typically develop over geological time. The T -axis of the moderate earthquakes further matches the extension direction previously measured on post early Pliocene (∼>3 Ma) faults. We conclude that although similar deformation has likely taken place in this area during recent geologic time, it does not appear to have coalesced during this time. The deformation may represent the southernmost part of a diffuse boundary, the Western Puerto Rico Deformation Boundary, which accommodates differential movement between the Puerto Rico and Hispaniola arc blocks. This differential movement is possibly driven by the differential seismic coupling along the Puerto Rico—Hispaniola subduction zone. We propose that the compositional heterogeneity across the island arc retards the process of focusing the deformation into a single fault. Given the evidence presented here, we should not expect a single large event in this area but similar diffuse sequences in the future.

Puerto Rico↗

Prospective and retrospective evaluation of the U.S. Geological Survey public aftershock forecast for the 2019-2021 Southwest Puerto Rico Earthquake and aftershocks

The M w "> M w 6.4 Southwest Puerto Rico Earthquake of 7 January 2020 was accompanied by a robust fore‐ and aftershock sequence. The U.S. Geological Survey (USGS) has issued regular aftershock forecasts for more than a year since the mainshock, available on a public webpage. Forecasts were accompanied by interpretive and informational material, published in English and Spanish. Informational products included narrative “scenarios” for how the aftershock sequence could play out, infographics, and a report on the potential duration of the aftershock sequence through the next decade. Forecasts are based on the epidemic‐type aftershock sequence (ETAS) model and generated using the USGS AftershockForecaster software—an interactive graphical user interface built on the OpenSHA platform ( Field et al. , 2003 ). The initial forecast is based on past sequences in similar tectonic environments; subsequent forecasts are tuned to the ongoing sequence via Bayesian model updating. Probabilistic aftershock forecasts for the next day, week, month, and year were publicly released and archived at a daily to monthly tempo, allowing for a truly prospective test of the forecast. Here, we evaluate the forecast over the first year of the recorded aftershocks. The ETAS‐based forecast performed well overall, successfully capturing both the chance of having at least one earthquake of a given magnitude in a forecast interval as well as the non‐Poissonian distribution of the total number of aftershocks within an interval. A retrospective analysis shows that the ETAS model is a substantial improvement over the existing Reasenberg and Jones (1989) forecast model. The exercise also reveals some limitations of the current model, in particular, with respect to nonstationarities in the aftershock magnitude distribution and model parameters throughout the evolving sequence.

Puerto Rico↗

Ground failure triggered by the 7 January 2020 M6.4 Puerto Rico earthquake

The 7 January 2020 M 6.4 Puerto Rico earthquake, the mainshock of an extended earthquake sequence, triggered significant ground failure. In this study, we detail the ground failure that occurred based largely on a postearthquake field reconnaissance campaign that we conducted. We documented more than 300 landslides, mainly rock falls that were concentrated in areas where peak ground acceleration (PGA) exceeded 30% g ; sparse smaller landslides occurred in highly susceptible areas more than 50 km from the epicenter and at PGA values <10% g . Though some of the largest mass movements were in natural slopes, rock falls in road cuts had more impact because they caused widespread transportation disruption. Some structures also were damaged by landslides, but no landslide‐related fatalities were reported. Liquefaction and related lateral spreading were severe in some areas, causing damage to residential and commercial structures and a power plant. Most of the liquefaction occurred in coastal areas where shaking exceeded 50% g , though some of the most damaging instances were in Ponce, where shaking estimates were as low as 20% g . In this article, we summarize the most notable ground failures and detail the overall patterns that we observed. The observations and compiled inventory datasets presented here are valuable because no island‐wide hazard maps for earthquake‐triggered landslides or liquefaction have been developed for Puerto Rico, and postevent inventories are vital for such an effort. In general, the datasets presented here contribute to a global understanding of coseismic ground‐failure hazard by presenting detailed observations for a relatively moderate ground‐failure event with well‐constrained shaking estimates; current models and expectations are biased by the tendency to collect detailed datasets mainly for exceptional events.

Puerto Rico↗

Forest cover lessens the impact of drought on streamflow in Puerto Rico

Tropical regions are experiencing high rates of forest cover loss coupled with changes in the volume and timing of rainfall. These shifts can compromise streamflow and water provision, highlighting the need to identify how forest cover influences streamflow generation under variable rainfall conditions. Although rainfall is the key driver of streamflow regimes, the role of forests is less clear, particularly in tropical regions where forest loss is an ongoing risk. Forest cover loss alters evapotranspiration, rainfall infiltration and storage, and may increase stream ecosystem vulnerability to rainfall extremes. Puerto Rico, an island with spatially heterogenous forest cover and a marked geographic rainfall gradient, is projected to experience more frequent droughts and flash flooding. Using 15-minute streamflow data collected between 2005 and 2016 from 20 USGS stream gages and 3-hourly Multi-Source Weighted-Ensemble Precipitation rainfall estimates, we utilized flow-duration curves and linear mixed regression models to examine the role of forest cover in regulating the timing and volume of streamflow. The mixed model approach helps to account for differences in watershed characteristics. We determined the effects of rainfall and forest cover on low and peak flows in Puerto Rican streams, then evaluated changes in these relationships under dry and wet antecedent rainfall conditions. Watersheds with high forest cover had consistently greater low and peak streamflow than deforested ones under all rainfall conditions, although the effect was more marked during wet antecedent conditions, suggesting that peak flow is largely the result of saturation excess overland flow. During dry antecedent rainfall conditions, highly forested watersheds had higher streamflow than deforested ones, suggesting greater hillslope storage and release may also be at play. Our results demonstrate that forest cover generated a net increase in hillslope infiltration and storage and may lessen drought impacts on streamflow in Puerto Rico. Resilience to prolonged drought may be limited by finite water storage potential in this steep, mountainous setting, highlighting maintenance of forest cover as an important water management strategy to increase infiltration.

Puerto Rico↗

Puerto Rico’s state of the climate 2014-2021

The climate of Puerto Rico is influenced by the changing global climate. The following chapters present the current knowledge of the geophysical and chemical drivers and signals of global climate change as they affect the climate of Puerto Rico and influence the climate-dependent services, risks, and vulnerabilities that govern human well-being. These include sustainable economic development, delivery of ecosystem services, the conservation of natural and cultural resources, resiliency in built and natural systems, and food security. The chapters draw on global expertise of land, atmosphere, and ocean geophysical interactions associated with increasing greenhouse gases that drive global warming and on local scientific expertise, data, observations, and modeled projections. They present the global warming scenario (section 1), the contribution of Puerto Rico to global climate change as GHG emissions and aerosols (section 2), the context of natural climate variability (section 3), observed and projected trends in temperature (section 4), rainfall (section 5), sea level rise (section 6), ocean acidification and sea surface warming (section 7), and the expected implications of warming climate on tropical cyclones affecting Puerto Rico (section 8).

Puerto Rico↗

Seismotectonic analysis of the 2019–2020 Puerto Rico sequence: The value of absolute earthquake relocations in improved interpretations of active tectonics

We present a new catalog of calibrated earthquake relocations from the 2019–2020 Puerto Rico earthquake sequence related to the 7 January 2020 M w 6.4 earthquake that occurred offshore of southwest Puerto Rico at a depth of 15.9 km. Utilizing these relocated earthquakes and associated moment tensor solutions, we can delineate several distinct fault systems that were activated during the sequence and show that the M w 6.4 mainshock may have resulted from positive changes in Coulomb stress from earlier events. Seismicity and mechanisms define (1) a west–southwest (∼260°) zone of seismicity comprised of largely sinistral strike‐slip and oblique‐slip earthquakes that mostly occurs later in the sequence and to the west of the mainshock, (2) an area of extensional faulting that includes the mainshock and occurs largely within the mainshock’s rupture area, and (3) an north–northeast (∼30°)‐striking zone of seismicity, consisting primarily of dextral strike‐slip events that occurs before and following the mainshock and generally above (shallower than) the normal‐faulting events. These linear features intersect within the M w 6.4 mainshock’s fault plane in southwest Puerto Rico. In addition, we show that earthquake relocations for M 4+ normal‐faulting events, when traced along their fault planes, daylight along east–west‐trending bathymetric features offshore of southwest Puerto Rico. Correlation of these normal‐faulting events with bathymetric features suggests an active fault system that may be a contributor to previously uncharacterized seismic hazards in southwest Puerto Rico.

Puerto Rico↗

Anomalous noble gas solubility in liquid cloud water: Possible implications for noble gas temperatures and cloud physics

The noble gas temperature climate proxy is an established tool that has previously been applied to determine the source of groundwater recharge, however, unanswered questions remain. In fractured media (e.g., volcanic islands) recharge can be so rapid that groundwater is significantly depleted in heavy noble gases, indicating that the water has retained noble gas concentrations from higher elevations. Previous studies of rain samples have confirmed a match to patterns seen in fractured-rock groundwater for heavy noble gases along with a significant helium excess. Snow has been shown to be a credible source for both the helium excess and the observed heavy noble gas pattern. Here, liquid cloud water samples were collected at two mountainous sites and analyzed for noble gas concentrations. A pattern like that of rainwater was found. However, an analysis of diffusive uptake of noble gases into cloud water demonstrates that droplets of 1 mm diameter and smaller should be in constant solubility equilibrium with the atmosphere. To explain this, we present a novel hypothesis that relies on the assumption that liquid water consists of two types of water molecule clusters bounded by hydrogen bonds: a low-density ice-like structure and a high-density condensed structure. In this model, the pressure gradient near the surface of a droplet resulting from surface tension could allow for the formation of a surface layer that is rich in ice-like low density clusters. This can explain both the helium excess and the heavy noble gas depletion seen in the samples.

Puerto Rico↗

Analysis of body condition indices reveals different ecotypes of the Antillean manatee

Assessing the body condition of wild animals is necessary to monitor the health of the population and is critical to defining a framework for conservation actions. Body condition indices (BCIs) are a non-invasive and relatively simple means to assess the health of individual animals, useful for addressing a wide variety of ecological, behavioral, and management questions. The Antillean manatee ( Trichechus manatus manatus ) is an endangered subspecies of the West Indian manatee, facing a wide variety of threats from mostly human-related origins. Our objective was to define specific BCIs for the subspecies that, coupled with additional health, genetic and demographic information, can be valuable to guide management decisions. Biometric measurements of 380 wild Antillean manatees captured in seven different locations within their range of distribution were obtained. From this information, we developed three BCIs (BCI 1 = UG/SL, BCI 2 = W/SL 3 , BCI 3 = W/(SL*UG 2 )). Linear models and two-way ANCOVA tests showed significant differences of the BCIs among sexes and locations. Although our three BCIs are suitable for Antillean manatees, BCI 1 is more practical as it does not require information about weight, which can be a metric logistically difficult to collect under particular circumstances. BCI 1 was significantly different among environments, revealing that the phenotypic plasticity of the subspecies have originated at least two ecotypes—coastal marine and riverine—of Antillean manatees.

Puerto Rico↗

Experimental warming and its legacy effects on root dynamics following two hurricane disturbances in a wet tropical forest

Tropical forests are expected to experience unprecedented warming and increases in hurricane disturbances in the coming decades; yet, our understanding of how these productive systems, especially their belowground component, will respond to the combined effects of varied environmental changes remains empirically limited. Here we evaluated the responses of root dynamics (production, mortality, and biomass) to soil and understory warming (+4°C) and after two consecutive tropical hurricanes in our in situ warming experiment in a tropical forest of Puerto Rico: Tropical Responses to Altered Climate Experiment (TRACE). We collected minirhizotron images from three warmed plots and three control plots of 12 m 2 . Following Hurricanes Irma and María in September 2017, the infrared heater warming treatment was suspended for repairs, which allowed us to explore potential legacy effects of prior warming on forest recovery. We found that warming significantly reduced root production and root biomass over time. Following hurricane disturbance, both root biomass and production increased substantially across all plots; the root biomass increased 2.8-fold in controls but only 1.6-fold in previously warmed plots. This pattern held true for both herbaceous and woody roots, suggesting that the consistent antecedent warming conditions reduced root capacity to recover following hurricane disturbance. Root production and mortality were both related to soil ammonium nitrogen and microbial biomass nitrogen before and after the hurricanes. This experiment has provided an unprecedented look at the complex interactive effects of disturbance and climate change on the root component of a tropical forested ecosystem. A decrease in root production in a warmer world and slower root recovery after a major hurricane disturbance, as observed here, are likely to have longer-term consequences for tropical forest responses to future global change.

Puerto Rico↗

Development of regression equations for the estimation of the magnitude and frequency of floods at rural, unregulated gaged and ungaged streams in Puerto Rico through water year 2017

The methods of computation and estimates of the magnitude of flood flows were updated for the 50-, 20-, 10-, 4-, 2-, 1-, 0.5-, and 0.2-percent chance exceedance levels for 91 streamgages on the main island of Puerto Rico by using annual peak-flow data through 2017. Since the previous flood frequency study in 1994, the U.S. Geological Survey has collected additional peak flows at additional streamgages, and Puerto Rico has experienced numerous flood events. This updated study was performed using longer annual peak-flow datasets from more stations to provide more representative equations to predict flood flows. Screening criteria for these streamgages included 10 or more years of annual peak-flow data, unregulated flow, and less than 10 percent impervious drainage area. The magnitude and frequency of floods at selected streamgages in Puerto Rico were estimated using updated methods outlined in Bulletin 17C. The new procedures include a regional skew analysis that incorporates Bayesian regression techniques, the Expected Moments Algorithm to better represent missing record and estimate parameters of the log-Pearson Type III distribution, and the Multiple Grubbs-Beck test for low outlier detection. Regional regression equations were developed to estimate peak-flow statistics at ungaged locations by using selected basin and climatic characteristics as explanatory variables. These variables were determined from digital spatial datasets and geographic information systems by using the most recent data available. Ordinary least-squares regression techniques were used to filter the basin characteristics and determine two separate regions, region 1 (west) and region 2 (east), based on residuals. A generalized least-squares procedure was used to account for cross-correlation of sites and develop the final set of equations that have drainage area as the only explanatory variable. The average standard errors of prediction ranged from 18.7 to 46.7 percent in region 1 and 33.4 to 57.6 percent in region 2 for all annual exceedance probabilities (AEPs) examined. The updated statistics showed a greater accuracy of prediction when compared to those from the previous study using drainage area as the only explanatory variable for all AEPs examined in region 1 and the 0.01 and 0.002 AEP flows for region 2. When compared to equations developed in the previous study that have drainage area, mean annual rainfall, and (or) depth-to-rock as explanatory variables, the updated statistics show a greater accuracy of prediction in region 1 at AEP flows of 0.02 and lower (that is, higher flows). Those developed for region 2 do not show a greater accuracy of prediction for any AEP flows when compared to the equations having multiple explanatory variables in the previous study. The calculated regression equations, basin characteristics, and at-site statistics will be incorporated into the U.S. Geological Survey web application, StreamStats ( https://streamstats.usgs.gov/ss/ ). This application allows users to select a location on a stream, whether gaged or ungaged, to obtain estimates of basin characteristics and flow statistics.

Puerto Rico↗

Rigorously valuing the coastal hazard risks reduction provided by potential coral reef restoration in Florida and Puerto Rico

The restoration of coastal habitats, particularly coral reefs, can reduce risks by decreasing the exposure of coastal communities to flooding hazards. In the United States, the protective services provided by coral reefs were recently assessed in social and economic terms, with the annual protection provided by U.S. coral reefs off the coasts of the State of Florida and the Commonwealth of Puerto Rico estimated to be more than 9,800 people and $859 million (2010 U.S. dollars). Hurricanes Irma and Maria in 2017 caused widespread damage to coral reefs in the State of Florida and the Commonwealth of Puerto Rico. Here we combine engineering, ecologic, geospatial, social, and economic data and tools to provide a rigorous valuation of where potential coral reef restoration could decrease the hazard faced by Florida and Puerto Rico’s reef-fronted coastal communities. The three restoration scenarios considered: (1) Ecological restoration, ‘E25’, which assumes planting 0.25-meter (m)-high corals on a (cross-shore) 25-m-wide reef; (2) Structural plus ecological, ‘S25’, which assumes emplacing a 1.00-m high structure with 0.25-m high corals on top on a 25 m wide reef; and (3) structural plus ecological, ‘S05’, which assumes emplacing a 1.00-m high structure with 0.25-m high corals on top on a 5 m wide reef. Planted corals are assumed to increase hydrodynamic roughness, thereby dissipating incident wave energy and decreasing flooding potential. We used a standardized approach to ‘place’ potential restoration projects throughout the whole (linear) extent of reefs bordering Florida and Puerto Rico to identify where coral reef restoration could be useful for meeting flood reduction benefits. We always sited potential restoration projects within the existing distribution of reefs even though many sites were far (kilometers [km]) offshore and some sites were relatively deep (up to 7 m depth). We followed risk-based valuation approaches to map flood zones at 10-square-meter resolution along all 980 km of Florida and Puerto’s Rico reef-lined shorelines for the three potential coral reef restoration scenarios and compare them to the flood zones without coral reef restoration. We quantified the potential coastal flood risk reduction provided by coral reef restoration using the latest information from the U.S. Census Bureau, Federal Emergency Management Agency, and Bureau of Economic Analysis for return-interval storm events. Using the damages associated with each storm probability, we also calculate the change in annual expected damages, a measure of the annual protection gained because of coral reef restoration. We found that the benefits of reef restoration off Florida and Puerto Rico are spatially highly variable. In most areas, we found little or no benefit from reef restoration (for example, restoration sites were far offshore or deep). However, there were a number of key areas where reef restoration could have substantial benefits for flood risk reduction. In particular, we estimated the protection gained by Florida and Puerto Rico’s coral reefs from coral reef restoration to result in: Avoided flooding to more than 5.6 square kilometers (2.16 square miles) of land annually; Avoided flooding affecting more than 3,100 people annually; Avoided direct damages of more than $124.2 million to more than 890 buildings annually; and Avoided indirect damages to more $148.7 million in economic activity owing to housing and business damage annually. Thus, the annual value of flood risk reduction provided by potential coral reef restoration in Florida and Puerto Rico is more than 3,100 people and $272.9 million (2010 U.S. dollars) in economic activity. These data provide stakeholders and decision makers with a spatially explicit, rigorous valuation of how, where, and when potential coral reef restoration in Florida and Puerto Rico can increase critical coastal storm flood reduction benefits. These results help identify areas where reef management, recovery, and restoration could potentially help reduce the risk to, and increase the resiliency of, Florida and Puerto Rico’s coastal communities.

Florida↗

Rigorously valuing the impact of Hurricanes Irma and Maria on coastal hazard risks in Florida and Puerto Rico

The degradation of coastal habitats, particularly coral reefs, raises risks by increasing the exposure of coastal communities to flooding hazards. In the United States, the physical protective services provided by coral reefs were recently assessed in social and economic terms, with the annual protection provided by U.S. coral reefs off the coasts of the State of Florida and the Commonwealth of Puerto Rico estimated to be more than 9,800 people and $859 million (2010 U.S. dollars). Hurricanes Irma and Maria in 2017 caused widespread damage to coral reefs in the State of Florida and the Commonwealth of Puerto Rico. These damages were measured in post-storm surveys of reefs and assessed in terms of their impact on reef condition and height, which are critical parameters for evaluating the coastal defense benefits of reefs. We combined engineering, ecologic, geospatial, social, and economic data and tools to value the increased risks in Florida and Puerto Rico from hurricane-induced damages to their adjacent coral reefs. We followed risk-based valuation approaches to map flooding at 10-square-meter resolution along all 980 kilometers of Florida and Puerto Rico’s reef-lined shorelines considering reef condition before (undamaged) and after (damaged) the 2017 hurricanes. We quantified the coastal flood risk increase caused by the hurricane-induced damage to the coral reefs using the latest information from the U.S. Census Bureau, Federal Emergency Management Agency, and Bureau of Economic Analysis for return-interval storm events. Using the damages associated with each storm probability, we also calculated the change in annual expected damages, a measure of the annual protection lost because of the reef damage caused by the 2017 hurricanes. We found that the damages to the coral reefs off Florida and Puerto Rico from Hurricanes Irma and Maria increased future risks significantly. In particular, we estimated the protection lost by Florida and Puerto Rico’s coral reefs from the 2017 hurricanes to result in: Increased flooding to more than 10.72 square kilometers (4.14 square miles) of land annually; Increased flooding affecting more than 4,300 people annually; Increased direct damages of more than $57.2 million to more than 1,800 buildings annually; and Increased indirect damages to more $124.3 million in economic activity owing to housing and business damage annually. Thus, the annual value of increased flood risk caused by the damage to Florida and Puerto Rico’s coral reefs from hurricanes in 2017 is more than 4,300 people and $181.5 mil-lion (2010 U.S. dollars) in economic impacts. These data provide stakeholders and decision makers with a spatially explicit, rigorous valuation of how, where, and when the damage from the 2017 hurricanes decreased critical coastal storm flood reduction benefits to Florida and Puerto Rico’s coral reefs. These results help identify areas where reef management, recovery, and restoration could potentially help reduce the risk to, and increase the resiliency of, Florida and Puerto Rico’s coastal communities.

Florida↗

Climate change is creating a mismatch between protected areas and suitable habitats for frogs and birds in Puerto Rico

Climate change is altering the spatial distribution of many species around the world. In response, we need to identify and protect suitable areas for a large proportion of the fauna so that they persist through time. This exercise must also evaluate the ability of existing protected areas to provide safe havens for species in the context of climate change. Here, we combined passive acoustic monitoring, semi-automatic species identification models, and species distribution models of 21 bird and frog species based on past (1980–1989), present (2005–2014), and future (2040–2060) climate scenarios to determine how species distributions relate to the current distribution of protected areas in Puerto Rico. Species detection/non-detection data were acquired across ~ 700 sampling sites. We developed always-suitable maps that characterized suitable habitats in all three time periods for each species and overlaid these maps to identify regions with high species co-occurrence. These distributions were then compared with the distribution of existing protected areas. We show that Puerto Rico is projected to become dryer by 2040–2060, and precipitation in the warmest quarter was among the most important variables affecting bird and frog distributions. A large portion of always-suitable areas (ASA) is outside of protected areas (> 80%), and the percent of protected areas that overlaps with always-suitable areas is larger for bird (75%) than frog (39%) species. Our results indicate that present protected areas will not suffice to safeguard bird and frog species under climate change; however, the establishment of larger protected areas, buffer zones, and connectivity between protected areas may allow species to find suitable niches to withstand environmental changes.

Puerto Rico↗

Estimated water withdrawals and use in Puerto Rico, 2015

Water withdrawals and use in Puerto Rico for 2015 were estimated at 2,372 million gallons per day (Mgal/d), which was 21 percent less than withdrawals and use for 2010. The 2015 total water withdrawal and use estimates were the lowest since 1990 and coincided with a substantial decline of 25 percent in saline-water withdrawals for thermoelectric-power cooling processes from 2010 to 2015. Freshwater withdrawals were 671 Mgal/d, or 28 percent of total water withdrawals, and saline-water withdrawals were 1,701 Mgal/d, or 72 percent of total withdrawals. Fresh surface-water withdrawals were estimated at 548 Mgal/d, 10 percent less than in 2010, whereas fresh groundwater withdrawals were estimated at 122 Mgal/d, 2 percent less than in 2010. Saline surface-water withdrawals were 25 percent less than in 2010. Freshwater withdrawals were greatest for public-supply water and irrigation in 2015 and, combined, accounted for 98 percent of Puerto Rico’s total freshwater withdrawals. Withdrawals in 2015 for public-supply water (576 Mgal/d) were 14 percent lower and withdrawals for irrigation (78 Mgal/d) were 104 percent greater than in 2010, possibly because of drought conditions in agricultural counties along the south and southeast coasts in 2015. The sources for public-supply water withdrawals in 2015 included surface water (88 percent) and groundwater (12 percent). Withdrawals for other uses, which account for the remaining 2 percent of Puerto Rico’s total freshwater withdrawals, were lower in 2015 than in 2010; specifically, withdrawals for domestic self-supplied use decreased by 78 percent, industrial withdrawals decreased by 15 percent, and withdrawals for livestock decreased by 25 percent. Freshwater withdrawals for thermoelectric power and mining were greater in 2015 than in 2010, increasing by 23 percent and 5 percent, respectively. The total population of Puerto Rico decreased by 7 percent from 2010 to 2015, from 3.73 million people in 2010 to 3.47 million people in 2015. The number of people who obtained potable water from public-supply water facilities in 2015 was about 3.47 million, or about 100 percent of the population of Puerto Rico. Public-supply water deliveries for domestic use accounted for 338 Mgal/d in 2015, which is 47 percent greater than in 2010, indicating an increase in domestic per capita use from 62 to 98 gallons per person per day from 2010 to 2015. Domestic self-supplied withdrawals were estimated at 0.52 Mgal/d in 2015, for an estimated 4,708 people (less than 1 percent of Puerto Rico’s population). All domestic self-supplied withdrawals were assumed to be from groundwater sources. Irrigation freshwater withdrawals were 78 Mgal/d in 2015 and accounted for 12 percent of the total freshwater withdrawals for all uses. Surface-water deliveries from irrigation districts accounted for 44 percent of total irrigation withdrawals, whereas groundwater withdrawals accounted for 56 percent. About 37,000 acres were irrigated in 2015, a decrease of 11 percent or about 4,000 acres compared to 2010. About 99 percent of the acreage was irrigated by micro-irrigation and sprinkler systems in 2015. About 65 percent of the irrigation withdrawals were accounted for by four municipalities: Santa Isabel, Salinas, Lajas, and Juana Díaz. Altogether, freshwater withdrawals for livestock, industrial, mining, and thermoelectric power accounted for 2 percent (16.2 Mgal/d) of freshwater withdrawals for all uses, 9 percent less than in 2010. About 71 percent of the freshwater withdrawn for these categories was from groundwater sources. In 2015, 50 percent of the total freshwater withdrawn in Puerto Rico was apportioned to six municipalities: Arecibo, Trujillo Alto, Toa Alta, Villalba, Aguada, and Mayagüez. Arecibo accounted for about 18 percent of the total freshwater withdrawals, predominantly for public-supply water use. Trujillo Alto, Toa Alta, Villalba, Aguada, and Mayagüez accounted for about 32 percent (213 Mgal/d) of the total freshwater withdrawals, which were predominantly for public-supply water uses. Withdrawals in some of these municipalities are subsequently distributed to other municipalities such as those in the San Juan metro area. The Puerto Rico Aqueduct and Sewer Authority water service area for the San Juan metro area (referred to as W–102) accounted for about 28 percent of the total water delivered from public-supply water facilities to domestic users, which includes about 34 percent of the total population of Puerto Rico.

Puerto Rico↗