Stable-isotope and solute-chemistry approaches to flow characterization in a forested tropical watershed, Luquillo Mountains, Puerto Rico

- Stream dissolved solutes increased with decreasing elevation in watershed.
- Hydrograph separation showed brief groundwater input then dilution of solutes.
- Stable isotopes indicated quickflow was 36% of the storm hydrograph.
- Deuterium excess is a good tracer of seasonal watershed response and transit time.
- Streamflow typically comprises ∼25% quickflow, ∼75% groundwater baseflow.

The prospect of changing climate has led to uncertainty about the resilience of forested mountain watersheds in the tropics. In watersheds where frequent, high rainfall provides ample runoff, we often lack understanding of how the system will respond under conditions of decreased rainfall or drought. Factors that govern water supply, such as recharge rates and groundwater storage capacity, may be poorly quantified. This paper describes 8-year data sets of water stable isotope composition (δ2H and δ18O) of precipitation (4 sites) and a stream (1 site), and four contemporaneous stream sample sets of solute chemistry and isotopes, used to investigate watershed response to precipitation inputs in the 1780-ha Río Mameyes basin in the Luquillo Mountains of northeastern Puerto Rico. Extreme δ2H and δ18O values from low-pressure storm systems and the deuterium excess (d-excess) were useful tracers of watershed response in this tropical system. A hydrograph separation experiment performed in June 2011 yielded different but complementary information from stable isotope and solute chemistry data. The hydrograph separation results indicated that 36% of the storm rain that reached the soil surface left the watershed in a very short time as runoff. Weathering-derived solutes indicated near-stream groundwater was displaced into the stream at the beginning of the event, followed by significant dilution. The more biologically active solutes exhibited a net flushing behavior. The d-excess analysis suggested that streamflow typically has a recent rainfall component (∼25%) with transit time less than the sampling resolution of 7 days, and a more well-mixed groundwater component (∼75%). The contemporaneous stream sample sets showed an overall increase in dissolved solute concentrations with decreasing elevation that may be related to groundwater inputs, different geology, and slope position. A considerable amount of water from rain events runs off as quickflow and bypasses subsurface watershed flowpaths, and better understanding of shallow hillslope and deeper groundwater processes in the watershed will require sub-weekly data and detailed transit time modeling. A combined isotopic and solute chemistry approach can guide further studies to a more comprehensive model of the hydrology, and inform decisions for managing water supply with future changes in climate and land use.