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:
By utilizing functional relationships based on observations at plot or field scales, water quality models first compute surface runoff and then use it as the primary governing variable to estimate sediment and nutrient transport. When these models are applied at watershed scales, this serial model structure, coupling a surface runoff sub‐model with a water quality sub‐model, may be inappropriate because dominant hydrological processes differ among scales. A parallel modeling approach is proposed to evaluate how best to combine dominant hydrological processes for predicting water quality at watershed scales. In the parallel scheme, dominant variables of water quality models are identified based entirely on their statistical significance using time series analysis. Four surface runoff models of different model complexity were assessed using both the serial and parallel approaches to quantify the uncertainty on forcing variables used to predict water quality. The eight alternative model structures were tested against a 25‐year high‐resolution data set of streamflow, suspended sediment discharge, and phosphorous discharge at weekly time steps. Models using the parallel approach consistently performed better than serial‐based models, by having less error in predictions of watershed scale streamflow, sediment and phosphorus, which suggests model structures of water quantity and quality models at watershed scales should be reformulated by incorporating the dominant variables. The implication is that hydrological models should be constructed in a way that avoids stacking one sub‐model with one set of scale assumptions onto the front end of another sub‐model with a different set of scale assumptions.
more » « less- NSF-PAR ID:
- 10458767
- Publisher / Repository:
- Wiley Blackwell (John Wiley & Sons)
- Date Published:
- Journal Name:
- Hydrological Processes
- Volume:
- 34
- Issue:
- 2
- ISSN:
- 0885-6087
- Page Range / eLocation ID:
- p. 404-419
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract 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. -
more » « less
Abstract 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
Abstract Wildfire is a natural component of sagebrush (
Artemisia spp.) steppe rangelands that induces temporal shifts in plant community physiognomy, ground surface conditions, and erosion rates. Fire alteration of the vegetation structure and ground cover in these ecosystems commonly amplifies soil losses by wind‐ and water‐driven erosion. Much of the fire‐related erosion research for sagebrush steppe has focused on either erosion by wind over gentle terrain or water‐driven erosion under high‐intensity rainfall on complex topography. However, many sagebrush rangelands are geographically positioned in snow‐dominated uplands with complex terrain in which runoff and sediment delivery occur primarily in winter months associated with cold‐season hydrology. Current understanding is limited regarding fire effects on the interaction of wind‐ and cold‐season hydrologic‐driven erosion processes for these ecosystems. In this study, we evaluated fire impacts on vegetation, ground cover, soils, and erosion across spatial scales at a snow‐dominated mountainous sagebrush site over a 2‐year period post‐fire. Vegetation, ground cover, and soil conditions were assessed at various plot scales (8 m2to 3.42 ha) through standard field measures. Erosion was quantified through a network of silt fences (n = 24) spanning hillslope and side channel or swale areas, ranging from 0.003 to 3.42 ha in size. Sediment delivery at the watershed scale (129 ha) was assessed by suspended sediment samples of streamflow through a drop‐box v‐notch weir. Wildfire consumed nearly all above‐ground live vegetation at the site and resulted in more than 60% bare ground (bare soil, ash, and rock) in the immediate post‐fire period. Widespread wind‐driven sediment loading of swales was observed over the first month post‐fire and extensive snow drifts were formed in these swales each winter season during the study. In the first year, sediment yields from north‐ and south‐facing aspects averaged 0.99–8.62 t ha−1at the short‐hillslope scale (~0.004 ha), 0.02–1.65 t ha−1at the long‐hillslope scale (0.02–0.46 ha), and 0.24–0.71 t ha−1at the swale scale (0.65–3.42 ha), and watershed scale sediment yield was 2.47 t ha−1. By the second year post fire, foliar cover exceeded 120% across the site, but bare ground remained more than 60%. Sediment yield in the second year was greatly reduced across short‐ to long‐hillslope scales (0.02–0.04 t ha−1), but was similar to first‐year measures for swale plots (0.24–0.61 t ha−1) and at the watershed scale (3.05 t ha−1). Nearly all the sediment collected across all spatial scales was delivered during runoff events associated with cold‐season hydrologic processes, including rain‐on‐snow, rain‐on‐frozen soils, and snowmelt runoff. Approximately 85–99% of annual sediment collected across all silt fence plots each year was from swales. The high levels of sediment delivered across hillslope to watershed scales in this study are attributed to observed preferential loading of fine sediments into swale channels by aeolian processes in the immediate post‐fire period and subsequent flushing of these sediments by runoff from cold‐season hydrologic processes. Our results suggest that the interaction of aeolian and cold‐season hydrologic‐driven erosion processes is an important component for consideration in post‐fire erosion assessment and prediction and can have profound implications for soil loss from these ecosystems. © 2019 John Wiley & Sons, Ltd. -
Abstract. Watersheds are the fundamental Earth surface functioning units that connect the land to aquatic systems. Many watershed-scale models represent hydrological processes but not biogeochemical reactive transport processes. This has limited our capability to understand and predict solute export, water chemistry and quality, and Earth system response to changing climate and anthropogenic conditions. Here we present a recently developed BioRT-Flux-PIHM (BioRT hereafter) v1.0, a watershed-scale biogeochemical reactive transport model. The model augments the previously developed RT-Flux-PIHM that integrates land-surface interactions, surface hydrology, and abiotic geochemical reactions. It enables the simulation of (1) shallow and deep-water partitioning to represent surface runoff, shallow soil water, and deeper groundwater and of (2) biotic processes including plant uptake, soil respiration, and nutrient transformation. The reactive transport part of the code has been verified against the widely used reactive transport code CrunchTope. BioRT-Flux-PIHM v1.0 has recently been applied in multiple watersheds under diverse climate, vegetation, and geological conditions. This paper briefly introduces the governing equations and model structure with a focus on new aspects of the model. It also showcases one hydrology example that simulates shallow and deep-water interactions and two biogeochemical examples relevant to nitrate and dissolved organic carbon (DOC). These examples are illustrated in two simulation modes of complexity. One is the spatially lumped mode (i.e., two land cells connected by one river segment) that focuses on processes and average behavior of a watershed. Another is the spatially distributed mode (i.e., hundreds of cells) that includes details of topography, land cover, and soil properties. Whereas the spatially lumped mode represents averaged properties and processes and temporal variations, the spatially distributed mode can be used to understand the impacts of spatial structure and identify hot spots of biogeochemical reactions. The model can be used to mechanistically understand coupled hydrological and biogeochemical processes under gradients of climate, vegetation, geology, and land use conditions.more » « less
-
more » « less
Abstract A vector‐river network explicitly uses realistic geometries of river reaches and catchments for spatial discretization in a river model. This enables improving the accuracy of the physical properties of the modeled river system, compared to a gridded river network that has been used in Earth System Models. With a finer‐scale river network, resolving smaller‐scale river reaches, there is a need for efficient methods to route streamflow and its constituents throughout the river network. The purpose of this study is twofold: (1) develop a new method to decompose river networks into hydrologically independent tributary domains, where routing computations can be performed in parallel; and (2) perform global river routing simulations with two global river networks, with different scales, to examine the computational efficiency and the differences in discharge simulations at various temporal scales. The new parallelization method uses a hierarchical decomposition strategy, where each decomposed tributary is further decomposed into many sub‐tributary domains, enabling hybrid parallel computing. This parallelization scheme has excellent computational scaling for the global domain where it is straightforward to distribute computations across many independent river basins. However, parallel computing for a single large basin remains challenging. The global routing experiments show that the scale of the vector‐river network has less impact on the discharge simulations than the runoff input that is generated by the combination of land surface model and meteorological forcing. The scale of vector‐river networks needs to consider the scale of local hydrologic features such as lakes that are to be resolved in the network.
