Quantitative estimations of ecohydrological water partitioning into evaporation and transpiration remains mostly based on plot‐scale investigations that use well‐instrumented, small‐scale experimental catchments in temperate regions. Here, we attempted to upscale and adapt the conceptual tracer‐aided ecohydrology model STARRtropics to simulate water partitioning, tracer, and storage dynamics over daily time steps and a 1‐km grid larger‐scale (2565 km2) in a sparsely instrumented tropical catchment in Costa Rica. The model was driven by bias‐corrected regional climate model outputs and was simultaneously calibrated against daily discharge observations from 2 to 30 years at four discharge gauging stations and a 1‐year, monthly streamwater isotope record of 46 streams. The overall model performance for the best discharge simulations ranged in KGE values from 0.4 to 0.6 and correlation coefficients for streamflow isotopes from 0.3 to 0.45. More importantly, independent model‐derived transpiration estimates, point‐scale residence time estimates, and measured groundwater isotopes showed reasonable model performance and simulated spatial and temporal patterns pointing towards an overall model realism at the catchment scale over reduced performance in the headwaters. The simulated catchment system was dominated by low‐seasonality and high precipitation inputs and a marked topographical gradient. Climatic drivers overrode smaller, landcover‐dependent transpiration fluxes giving a seemingly homogeneous rainfall‐runoff dominance likely related to model input bias of rainfall isotopes, oversimplistic Potential Evapotranspiration (PET) estimates and averaged Leaf Area Index (LAI). Topographic influences resulted in more dynamic water and tracer fluxes in the headwaters that averaged further downstream at aggregated catchment scales. Modelled headwaters showed greater storage capacity by nearly an order of magnitude compared to the lowlands, which also favoured slightly longer residence times (>250 days) compared to superficially well‐connected groundwater contributing to shorter streamflow residence times (<150 days) in the lowlands. Our findings confirm that tracer‐aided ecohydrological modelling, even in the data‐scarce Tropics, can help gain a first, but crucial approximation of spatio‐temporal dynamics of how water is partitioned, stored and transported beyond the experimental catchment scale of only a few km2.
Recent studies have demonstrated that compartmentalized pools of water preferentially supply either plant transpiration (poorly mobile water) or streamflow and groundwater (highly mobile water) in some catchments, a phenomenon referred to as ecohydrologic separation. The omission of processes accounting for ecohydrologic separation in standard applications of hydrological models is expected to influence estimates of water residence times and plant water availability. However, few studies have tested this expectation or investigated how ecohydrologic separation alters interpretations of stores and fluxes of water within a catchment. In this study, we compare two rainfall‐runoff models that integrate catchment‐scale representations of transport, one that incorporates ecohydrologic separation and one that does not. The models were developed for a second‐order watershed at the H.J. Andrews Experimental Forest (Oregon, USA), the site where ecohydrologic separation was first observed, and calibrated against multiple years of stream discharge and chloride concentration. Model structural variations caused mixed results for differences in calibrated parameters and differences in storage between reservoirs. However, large differences in catchment storage volumes and fluxes arise when considering only mobile water. These changes influence interpreted residence times for streamflow‐generating water, demonstrating the importance of ecohydrologic separation in catchment‐scale water and solute transport.
more » « less- Award ID(s):
- 1652293
- NSF-PAR ID:
- 10460399
- Publisher / Repository:
- Wiley Blackwell (John Wiley & Sons)
- Date Published:
- Journal Name:
- Hydrological Processes
- Volume:
- 33
- Issue:
- 20
- ISSN:
- 0885-6087
- Page Range / eLocation ID:
- p. 2658-2675
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
more » « less
Abstract Large‐scale models often use a single grid to represent an entire catchment assuming homogeneity; the impacts of such an assumption on simulating evapotranspiration (ET) and streamflow remain poorly understood. Here, we compare hydrological dynamics at Shale Hills (PA, USA) using a complex model (spatially explicit, >500 grids) and a simple model (spatially implicit, two grids using “effective” parameters). We asked two questions:
What hydrological dynamics can a simple model reproduce at the catchment scale? What processes does it miss by ignoring spatial details? Results show the simple model can reproduce annual runoff ratios and ET, daily discharge peaks (e.g., storms, floods) but not discharge minima (e.g., droughts) under dry conditions. Neither can it reproduce different streamflow from the two sides of the catchment with distinct land surface characteristics. The similar annual runoff ratios between the two models indicate spatial details are not as important as climate in reproducing annual scale ET and discharge partitioning. Most of the calibrated parameters in the simple model are within the ranges in the complex model, except that effective porosity has to be reduced to 40% of the average porosity from the complex model. The form of the storage‐discharge relationship is similar. The effective porosity in the simple model however represents the dynamic and mobile water storage in the effective drainage area of the complex model that connects to the stream and contributes to high streamflow; it does not represent the passive, immobile water storage in the often disconnected uphill areas. This indicates that an additional uphill functioning unit is needed in the simple model to simulate the full spectrum of high‐low streamflow dynamics in natural catchments. -
null (Ed.)Abstract. Quantifying how vegetation mediates water partitioning at different spatialand temporal scales in complex, managed catchments is fundamental forlong-term sustainable land and water management. Estimations fromecohydrological models conceptualising how vegetation regulates theinterrelationships between evapotranspiration losses, catchment water storage dynamics, and recharge and runoff fluxes are needed to assess water availability for a range of ecosystem services and evaluate how these might change under increasing extreme events, such as droughts. Currently, the feedback mechanisms between water and mosaics of different vegetation and land cover are not well understood across spatial scales, and the effects of different scaleson the skill of ecohydrological models needs to be clarified. We used thetracer-aided ecohydrological model EcH2O-iso in an intensively monitored 66 km2 mixed land use catchment in northeastern Germany to quantify water flux–storage–age interactions at four model grid resolutions (250, 500, 750, and 1000 m). This used a fusion of field (including precipitation, soil water, groundwater, and stream isotopes) and remote sensing data in the calibration. Multicriteria calibration across the catchment at each resolution revealed some differences in the estimation of fluxes, storages, and water ages. In general, model sensitivity decreased and uncertainty increased with coarser model resolutions. Larger grids were unable to replicate observed streamflow and distributed isotope dynamics in the way smaller pixels could. However, using isotope data in the calibration still helped constrain the estimation of fluxes, storage, and water ages at coarserresolutions. Despite using the same data and parameterisation for calibration at different grid resolutions, the modelled proportion of fluxes differed slightly at each resolution, with coarse models simulating higher evapotranspiration, lower relative transpiration, increased overland flow, and slower groundwater movement. Although the coarser resolutions also revealed higher uncertainty and lower overall model performance, the overall results were broadly similar. The study shows that tracers provide effective calibration constraints on larger resolution ecohydrological modelling and help us understand the influence of grid resolution on the simulation of vegetation–soil interactions. This is essential in interpreting associated uncertainty in estimating land use influence on large-scale “blue” (ground and surface water) and “green” (vegetation and evaporated water) fluxes, particularly for future environmental change.more » « less
-
more » « less
Abstract Considering heterogeneity in porous media pore size and connectivity is essential to predicting reactive solute transport across interfaces. However, exchange with less‐mobile porosity is rarely considered in surface water/groundwater recharge studies. Previous research indicates that a combination of pore‐fluid sampling and geoelectrical measurements can be used to quantify less‐mobile porosity exchange dynamics using the time‐varying relation between fluid and bulk electrical conductivity. For this study, we use macro‐scale (10 s of cm) advection–dispersion solute transport models linked with electrical conduction in COMSOL Multiphysics to explore less‐mobile porosity dynamics in two different types of observed sediment water interface porous media. Modeled sediment textures contrast from strongly layered streambed deposits to poorly sorted lakebed sands and cobbles. During simulated ionic tracer perturbations, a lag between fluid and bulk electrical conductivity, and the resultant hysteresis, is observed for all simulations indicating differential loading of pore spaces with tracer. Less‐mobile exchange parameters are determined graphically from these tracer time series data without the need for inverse numerical model simulation. In both sediment types, effective less‐mobile porosity exchange parameters are variable in response to changes in flow direction and fluid flux. These observed flow‐dependent effects directly impact local less‐mobile residence times and associated contact time for biogeochemical reaction. The simulations indicate that for the sediment textures explored here, less‐mobile porosity exchange is dominated by variable rates of advection through the domain, rather than diffusion of solute, for typical low‐to‐moderate rate (approximately 3–40 cm/day) hyporheic fluid fluxes. Overall, our model‐based results show that less‐mobile porosity may be expected in a range of natural hyporheic sediments and that changes in flowpath orientation and magnitude will impact less‐mobile exchange parameters. These temporal dynamics can be assessed with the geoelectrical experimental tracer method applied at laboratory and field scales.
-
more » « less
Abstract This study integrated spatially distributed field observations and soil thermal models to constrain the impact of frozen ground on snowmelt partitioning and streamflow generation in an alpine catchment within the Niwot Ridge Long‐Term Ecological Research site, Colorado, USA. The study area was comprised of two contrasting hillslopes with notable differences in topography, snow depth and plant community composition. Time‐lapse electrical resistivity surveys and soil thermal models enabled extension of discrete soil moisture and temperature measurements to incorporate landscape variability at scales and depths not possible with point measurements alone. Specifically, heterogenous snowpack thickness (~0–4 m) and soil volumetric water content between hillslopes (~0.1–0.45) strongly influenced the depths of seasonal frost, and the antecedent soil moisture available to form pore ice prior to freezing. Variable frost depths and antecedent soil moisture conditions were expected to create a patchwork of differing snowmelt infiltration rates and flowpaths. However, spikes in soil temperature and volumetric water content, as well as decreases in subsurface electrical resistivity revealed snowmelt infiltration across both hillslopes that coincided with initial decreases in snow water equivalent and early increases in streamflow. Soil temperature, soil moisture and electrical resistivity data from both wet and dry hillslopes showed that initial increases in streamflow occurred prior to deep soil water flux. Temporal lags between snowmelt infiltration and deeper percolation suggested that the lateral movement of water through the unsaturated zone was an important driver of early streamflow generation. These findings provide the type of process‐based information needed to bridge gaps in scale and populate physically based cryohydrologic models to investigate subsurface hydrology and biogeochemical transport in soils that freeze seasonally.
