Global-scale regionalization of hydrologic model parameters

Current state-of-the-art models typically applied at continental to global scales (hereafter called macroscale) tend to use a priori parameters, resulting in suboptimal streamflow (Q) simulation. For the first time, a scheme for regionalization of model parameters at the global scale was developed. We used data from a diverse set of 1787 small-to-medium sized catchments ( 10-10,000 km(2)) and the simple conceptual HBV model to set up and test the scheme. Each catchment was calibrated against observed daily Q, after which 674 catchments with high calibration and validation scores, and thus presumably good-quality observed Q and forcing data, were selected to serve as donor catchments. The calibrated parameter sets for the donors were subsequently transferred to 0.5 degrees grid cells with similar climatic and physiographic characteristics, resulting in parameter maps for HBV with global coverage. For each grid cell, we used the 10 most similar donor catchments, rather than the single most similar donor, and averaged the resulting simulated Q, which enhanced model performance. The 1113 catchments not used as donors were used to independently evaluate the scheme. The regionalized parameters outperformed spatially uniform (i.e., averaged calibrated) parameters for 79% of the evaluation catchments. Substantial improvements were evident for all major Koppen-Geiger climate types and even for evaluation catchments>5000 km distant from the donors. The median improvement was about half of the performance increase achieved through calibration. HBV with regionalized parameters outperformed nine state-of-the-art macroscale models, suggesting these might also benefit from the new regionalization scheme. The produced HBV parameter maps including ancillary data are available via

A coupled terrestrial and aquatic biogeophysical model of the Upper Merrimack River watershed, New Hampshire, to inform ecosystem services evaluation and management under climate and land-cover change

Accurate quantification of ecosystem services (ES) at regional scales is increasingly important for making informed decisions in the face of environmental change. We linked terrestrial and aquatic ecosystem process models to simulate the spatial and temporal distribution of hydrological and water quality characteristics related to ecosystem services. The linked model integrates two existing models (a forest ecosystem model and a river network model) to establish consistent responses to changing drivers across climate, terrestrial, and aquatic domains. The linked model is spatially distributed, accounts for terrestrial–aquatic and upstream–downstream linkages, and operates on a daily time-step, all characteristics needed to understand regional responses. The model was applied to the diverse landscapes of the Upper Merrimack River watershed, New Hampshire, USA. Potential changes in future environmental functions were evaluated using statistically downscaled global climate model simulations (both a high and low emission scenario) coupled with scenarios of changing land cover (centralized vs. dispersed land development) for the time period of 1980–2099. Projections of climate, land cover, and water quality were translated into a suite of environmental indicators that represent conditions relevant to important ecosystem services and were designed to be readily understood by the public. Model projections show that climate will have a greater influence on future aquatic ecosystem services (flooding, drinking water, fish habitat, and nitrogen export) than plausible changes in land cover. Minimal changes in aquatic environmental indicators are predicted through 2050, after which the high emissions scenarios show intensifying impacts. The spatially distributed modeling approach indicates that heavily populated portions of the watershed will show the strongest responses. Management of land cover could attenuate some of the changes associated with climate change and should be considered in future planning for the region

Sensitivity analysis and parameter estimation for distributed hydrological modeling: potential of variational methods

Variational methods are widely used for the analysis and control of computationally intensive spatially distributed systems. In particular, the adjoint state method enables a very efficient calculation of the derivatives of an objective function (response function to be analysed or cost function to be optimised) with respect to model inputs. In this contribution, it is shown that the potential of variational methods for distributed catchment scale hydrology should be considered. A distributed flash flood model, coupling kinematic wave overland flow and Green Ampt infiltration, is applied to a small catchment of the Thoré basin and used as a relatively simple (synthetic observations) but didactic application case. It is shown that forward and adjoint sensitivity analysis provide a local but extensive insight on the relation between the assigned model parameters and the simulated hydrological response. Spatially distributed parameter sensitivities can be obtained for a very modest calculation effort (~6 times the computing time of a single model run) and the singular value decomposition (SVD) of the Jacobian matrix provides an interesting perspective for the analysis of the rainfall-runoff relation. For the estimation of model parameters, adjoint-based derivatives were found exceedingly efficient in driving a bound-constrained quasi-Newton algorithm. The reference parameter set is retrieved independently from the optimization initial condition when the very common dimension reduction strategy (i.e. scalar multipliers) is adopted. Furthermore, the sensitivity analysis results suggest that most of the variability in this high-dimensional parameter space can be captured with a few orthogonal directions. A parametrization based on the SVD leading singular vectors was found very promising but should be combined with another regularization strategy in order to prevent overfitting

Geographic Information Systems (GIS)-based spatially distributed model for runoff routing

A method is proposed for routing spatially distributed excess precipitation over a watershed to produce runoff at its outlet. The land surface is represented by a (raster) digital elevation model from which the stream network is derived. A routing response function is defined for each digital elevation model cell so that water movement from cell to cell can be convolved to give a response function along a flow path and responses from all cells can be summed to give the outlet hydrograph. An example application of analysis of runoff on Waller Creek in Austin, Texas, is presented.Waller Creek Working Grou

Using the soil and water assessment tool to simulate the pesticide dynamics in the data scarce Guayas River Basin, Ecuador

Agricultural intensification has stimulated the economy in the Guayas River basin in Ecuador, but also affected several ecosystems. The increased use of pesticides poses a serious threat to the freshwater ecosystem, which urgently calls for an improved knowledge about the impact of pesticide practices in this study area. Several studies have shown that models can be appropriate tools to simulate pesticide dynamics in order to obtain this knowledge. This study tested the suitability of the Soil and Water Assessment Tool (SWAT) to simulate the dynamics of two different pesticides in the data scarce Guayas River basin. First, we set up, calibrated and validated the model using the streamflow data. Subsequently, we set up the model for the simulation of the selected pesticides (i.e., pendimethalin and fenpropimorph). While the hydrology was represented soundly by the model considering the data scare conditions, the simulation of the pesticides should be taken with care due to uncertainties behind essential drivers, e.g., application rates. Among the insights obtained from the pesticide simulations are the identification of critical zones for prioritisation, the dominant areas of pesticide sources and the impact of the different land uses. SWAT has been evaluated to be a suitable tool to investigate the impact of pesticide use under data scarcity in the Guayas River basin. The strengths of SWAT are its semi-distributed structure, availability of extensive online documentation, internal pesticide databases and user support while the limitations are high data requirements, time-intensive model development and challenging streamflow calibration. The results can also be helpful to design future water quality monitoring strategies. However, for future studies, we highly recommend extended monitoring of pesticide concentrations and sediment loads. Moreover, to substantially improve the model performance, the availability of better input data is needed such as higher resolution soil maps, more accurate pesticide application rate and actual land management programs. Provided that key suggestions for further improvement are considered, the model is valuable for applications in river ecosystem management of the Guayas River basin

Assessment of different modelling studies on the spatial hydrological processes in an arid alpine catchment

To assess the model description of spatial hydrological processes in the arid alpine catchment, SWAT and MIKE SHE were jointly applied in Yarkant River basin located in northwest China. Not only the simulated daily discharges at the outlet station but also spatiotemporal distributions of runoff, snowmelt and evapotranspiration were analyzed contrastively regarding modules' structure and algorithm. The simulation results suggested both models have their own strengths for particular hydrological processes. For the stream runoff simulation, the significant contributions of lateral interflow flow were only reflected in SWAT with a proportion of 41.4 %, while MIKE SHE simulated a more realistic distribution of base flow from groundwater with a proportion of 21.3 %. In snowmelt calculation, SWAT takes account of more available factors and got better correlations between snowmelt and runoff in temporal distribution, however, MIKE SHE presented clearer spatial distribution of snowpack because of fully distributed structure. In the aspect of water balance, less water was evaporated because of limitation of soil evaporation and less spatially distributed approach in SWAT, on another hand, the spatial distribution of actual evapotranspiration (ETa) in MIKE SHE clearly expressed influence of land use. Whether SWAT or MIKE SHE, without multiple calibrations, the model's limitation might bring in some biased opinions of hydrological processes in a catchment scale. The complementary study of combined results from multiple models could have a better understanding of overall hydrological processes in arid and scarce gauges alpine region

Information content of spatially distributed ground-based measurements for hydrologic-parameter calibration in mixed rain-snow mountain headwaters

Parameters in hydrologic models used in mixed rain-snow regions are often uncertain to calibrate and overfitted on streamflow. To contribute addressing these challenges, we used an algorithm that assesses modeling performances through time (Dynamic Identifiability Analysis) to quantify the information content of spatially distributed ground-based measurements for identifying optimal parameter values in the Precipitation Runoff Modeling System (PRMS) model. Including spatially distributed ground-based measurements in Identifiability Analysis allowed us to unambiguously estimate more parameter values than only using streamflow (seven parameters instead of two out of a pool of thirty-three). Peaks in information gain were obtained when using dew-point temperature to identify precipitation phase-partitioning parameters. Multi-attribute identifiability analysis also yielded optimal parameter values that were temporally less variable than those estimated using streamflow alone. Overall, identifying parameter values using ground-based measurements improved the simulation of key drivers of the surface-water budget, such as air temperature and precipitation-phase partitioning. However, parameters simulating surface-to-subsurface mass fluxes like snow accumulation and melt or evapotranspiration were poorly identified by any attribute and so emerged as key sources of predictive uncertainty for this distributed-parameter hydrologic model. This work demonstrates the value of expanded ground-based measurements for identifying parameters in distributed-parameter hydrologic models and so diagnosing their conceptual uncertainty across the water budget

Predictability of seasonal runoff in the Mississippi River basin

Recent advances in climate prediction and remote sensing offer the potential to improve long-lead streamflow forecasts and to provide better land surface state estimates at the time of forecast. We characterize predictability of runoff at seasonal timescales in the Mississippi River basin due to climatic persistence (represented by El Niño-Southern Oscillation and the Arctic Oscillation) and persistence related to the initial land surface state (soil moisture and snow). These climate and land surface state indicators, at varying lead times, are then used in a multiple linear regression to explain the variance of seasonal average runoff. Soil moisture dominates runoff predictability for lead times of 1 1/2 months, except in summer in the western part of the basin, where snow dominates. For the western part of the basin, the land surface state has a stronger predictive capability than climate indicators through leads of two seasons; climate indicators are more important in the east at lead times of one season or greater. Modest winter runoff predictability exists at a lead time of 3 seasons due to both climate and soil moisture, but this is in areas producing little runoff and is therefore of lessened importance. Local summer runoff predictability is limited to the western mountainous areas (generating high runoff) through a lead of 2 seasons. This could be useful to water managers in the western portion of the Mississippi River basin, because it suggests the potential to provide skillful forecast information earlier in the water year than currently used in operational forecasts

Tools for Assessing Climate Impacts on Fish and Wildlife

Climate change is already affecting many fish and wildlife populations. Managing these populations requires an understanding of the nature, magnitude, and distribution of current and future climate impacts. Scientists and managers have at their disposal a wide array of models for projecting climate impacts that can be used to build such an understanding. Here, we provide a broad overview of the types of models available for forecasting the effects of climate change on key processes that affect fish and wildlife habitat (hydrology, fire, and vegetation), as well as on individual species distributions and populations. We present a framework for how climate-impacts modeling can be used to address management concerns, providing examples of model-based assessments of climate impacts on salmon populations in the Pacific Northwest, fire regimes in the boreal region of Canada, prairies and savannas in the Willamette Valley-Puget Sound Trough-Georgia Basin ecoregion, and marten Martes americana populations in the northeastern United States and southeastern Canada. We also highlight some key limitations of these models and discuss how such limitations should be managed. We conclude with a general discussion of how these models can be integrated into fish and wildlife management