As is true of many tropical regions, northeastern Puerto Rico is an ecologically sensitive area with biological life that is highly elevation dependent on precipitation and temperature. Climate change has the potential to increase the risk of losing endemic species and habitats. Consequently, it is important to explore the pattern of trends in precipitation and temperature along an elevation gradient. Statistical derivatives of a frequently sampled dataset of precipitation and temperature at 20 sites along an elevation gradient of 1000 m in northeastern Puerto Rico were examined for trends from 2001 to 2013 with nonparametric methods accounting for annual periodic variations such as yearly weather cycles. Overall daily precipitation had an increasing trend of around 0.1 mm day 21 yr 21 . The driest months of the annual dry, early, and late rainfall seasons showed a small increasing trend in the precipitation (around day 21 yr 21 ). There was strong evidence that precipitation in the driest months of each rainfall season increased faster at higher elevations (0.02 mm day 21 more increase for 100-m elevation gain) and some evidence for the same pattern in precipitation in all months of the year but at half the rate. Temperature had a positive trend in the daily minimum (around 0.028C yr 21 ) and a negative trend in the daily maximum whose size is likely an order of magnitude larger than the size of the daily minimum trend. Physical mechanisms behind the trends may be related to climate change; longer-term studies will need to be undertaken in order to assess the future climatic trajectory of tropical forests.

This study quantitatively explores whether land cover changes have a substantive impact on simulated streamflow within the tropical island setting of Puerto Rico. The Precipitation Runoff Modeling System (PRMS) was used to compare streamflow simulations based on five static parameterizations of land cover with those based on dynamically varying parameters derived from four land cover scenes for the period 1953-2012. The PRMS simulations based on static land cover illustrated consistent differences in simulated streamflow across the island. It was determined that the scale of the analysis makes a difference: large regions with localized areas that have undergone dramatic land cover change may show negligible difference in total streamflow, but streamflow simulations using dynamic land cover parameters for a highly altered subwatershed clearly demonstrate the effects of changing land cover on simulated streamflow. Incorporating dynamic parameterization in these highly altered watersheds can reduce the predictive uncertainty in simulations of streamflow using PRMS. Hydrologic models that do not consider the projected changes in land cover may be inadequate for water resource management planning for future conditions.(KEY TERMS: land use/land cover change; urbanization; surface water hydrology; Precipitation Runoff Modeling System; geospatial analysis; Caribbean.)

The impact of Hurricane Maria on the U.S. Caribbean was used to study the causes of remotely-sensed spatial variation in the effects of (1) vegetation index loss and (2) landslide occurrence. The vegetation index is a measure of canopy ‘greenness’, a combination of leaf chlorophyll, leaf area, canopy cover and structure. A generalized linear model was made for each kind of effect, using idealized maps of the hurricane forces, along with three landscape characteristics that were significantly associated. In each model, one of these characteristics was forest fragmentation, and another was a measure of disturbance-propensity. For the greenness loss model, the hurricane force was wind, the disturbance-propensity measure was initial greenness, and the third landscape characteristic was fraction forest cover. For the landslide occurrence model, the hurricane force was rain, the disturbance-propensity measure was amount of land slope, and the third landscape characteristic was soil clay content. The model of greenness loss had a pseudo R2 of 0.73 and showed the U.S. Caribbean lost 31% of its initial greenness from the hurricane, with 51% lost from the initial in the Luquillo Experimental Forest (LEF) from Hurricane Maria along with Hurricane Irma. More greenness disturbance was seen in areas with less wind sheltering, higher elevation and topographic sides. The model of landslide occurrence had a pseudo R2 of 0.53 and showed the U.S. Caribbean had 34% of its area and 52% of the LEF area with a landslide density of at least one in 1 km2 from Hurricane Maria. Four experiments with parameters from previous storms of wind speed, storm duration, rainfall, and forest structure over the same storm path and topographic landscape were run as examples of possible future scenarios. While intensity of the storm makes by far the largest scenario difference, forest fragmentation makes a sizable difference especially in vulnerable areas of high clay content or high wind susceptibility. This study showed the utility of simple hurricane force calculations connected with landscape characteristics and remote-sensing data to determine forest susceptibility to hurricane effects.

A module to simulate glacier runoff, PRMSglacier, was added to PRMS (Precipitation Runoff Modeling System), a distributed‐parameter, physical‐process hydrological simulation code. The extension does not require extensive on‐glacier measurements or computational expense but still relies on physical principles over empirical relations as much as is feasible while maintaining model usability. PRMSglacier is validated on two basins in Alaska, Wolverine, and Gulkana Glacier basin, which have been studied since 1966 and have a substantial amount of data with which to test model performance over a long period of time covering a wide range of climatic and hydrologic conditions. When error in field measurements is considered, the Nash‐Sutcliffe efficiencies of streamflow are 0.87 and 0.86, the absolute bias fractions of the winter mass balance simulations are 0.10 and 0.08, and the absolute bias fractions of the summer mass balances are 0.01 and 0.03, all computed over 42 years for the Wolverine and Gulkana Glacier basins, respectively. Without taking into account measurement error, the values are still within the range achieved by the more computationally expensive codes tested over shorter time periods.

Abstract. The degree to which cloud immersion provides water in addition to rainfall, suppresses transpiration, and sustains tropical montane cloud forests (TMCFs) during rainless periods is not well understood. Climate and land use changes represent a threat to these forests if cloud base altitude rises as a result of regional warming or deforestation. To establish a baseline for quantifying future changes in cloud base, we installed a ceilometer at 100 m altitude in the forest upwind of the TMCF that occupies an altitude range from ∼ 600 m to the peaks at 1100 m in the Luquillo Mountains of eastern Puerto Rico. Airport Automated Surface Observing System (ASOS) ceilometer data, radiosonde data, and CloudAerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO) satellite data were obtained to investigate seasonal cloud base dynamics, altitude of the trade-wind inversion (TWI), and typical cloud thickness for the surrounding Caribbean region. Cloud base is rarely quantified near mountains, so these results represent a first look at seasonal and diurnal cloud base dynamics for the TMCF. From May 2013 to August 2016, cloud base was lowest during the midsummer dry season, and cloud bases were lower than the mountaintops as often in the winter dry season as in the wet seasons. The lowest cloud bases most frequently occurred at higher elevation than 600 m, from 740 to 964 m. The Luquillo forest low cloud base altitudes were higher than six other sites in the Caribbean by ∼ 200-600 m, highlighting the importance of site selection to measure topographic influence on cloud height. Proximity to the oceanic cloud system where shallow cumulus clouds are seasonally invariant in altitude and cover, along with local trade-wind orographic lifting and cloud formation, may explain the dry season low clouds. The results indicate that climate change threats to low-elevation TMCFs are not limited to the dry season; changes in synoptic-scale weather patterns that increase frequency of drought periods during the wet seasons (periods of higher cloud base) may also impact ecosystem health.