A near-global, high resolution land surface parameter dataset for the variable infiltration capacity model

Hydrologic models predict the spatial and temporal distribution of water and energy at the land surface. Currently, parameter availability limits global-scale hydrologic modelling to very coarse resolution, hindering researchers from resolving fine-scale variability. With the aim of addressing this problem, we present a set of globally consistent soil and vegetation parameters for the Variable Infiltration Capacity (VIC) model at 1/16° resolution (approximately 6 km at the equator), with spatial coverage from 60°S to 85°N. Soil parameters derived from interpolated soil profiles and vegetation parameters estimated from space-based MODIS measurements have been compiled into input files for both the Classic and Image drivers of the VIC model, version 5. Geographical subsetting codes are provided, as well. Our dataset provides all necessary land surface parameters to run the VIC model at regional to global scale. We evaluate VICGlobal's ability to simulate the water balance in the Upper Colorado River basin and 12 smaller basins in the CONUS, and their ability to simulate the radiation budget at six SURFRAD stations in the CONUS.

Hydro-Geomorphologic-Based Water Budget at Event Time-Scale in A Mediterranean Headwater Catchment (Southern Italy)

The Ciciriello catchment is a 3 km2 drainage sub-basin of the Bussento river basin, located in the southern part of the Campania Region (Southern Italy). Since 2012, this catchment has been studied using an interdisciplinary approach—geomorphological, hydrogeological, and hydrological—and a hydro-chemical monitoring system. Following previous research, the aim of this paper is to calibrate, on this catchment, the hydrologic parameters for a water budget at event time-scales using the HEC-HMS model, adopting object-based hydro-geomorphological class features. Firstly, lumped modeling was performed to calibrate the hydrologic parameters from 20 observed hydrographs at the downstream monitoring station of the Ciciriello catchment. Then, physical-based rainfall–runoff modeling was conducted using three different procedures: (1) applying the recession coefficients to each outlet with a newly defined hydro-geomorphologic index (HGmI); (2) assessing the storage coefficient for each sub-basin as a weighted mean of HGmI; and (3) using the storage coefficient associated with the largest HGmI in the sub-basin. The adopted procedures were tested using diverse goodness-of-fit indices, resulting in good performance when the object-based hydro-geomorphotypes were used for the parameter calibration. The adopted procedure can thus contribute to improvements in rainfall–runoff and water budget modeling in similar ungauged catchments in Mediterranean, hilly, and forested landscapes.

Impacts of Climate Change on Hydroclimatic Conditions of U.S. National Forests and Grasslands

The conterminous United States includes national forests and grasslands that provide ecological, social, economic, recreational, and aesthetic services. Future climate change can alter long-term hydroclimatic conditions of national forests and grasslands and lead to negative consequences. This study characterizes shifts in hydroclimatology and basin characteristics of US National Forests (NFs) and National Grasslands (NGs) in response to climate change over the 21st century under the DRY, MIDDLE, and WET climate models with the representative concentration pathway (RCP) 8.5 emission scenario. Climatic projections for three climate models ranging from the driest to wettest conditions were obtained from the Multivariate Adaptive Constructed Analogs (MACA) dataset. Then, the variable infiltration capacity (VIC) model was used to model hydrological responses of the selected future climates. Changes in regional hydroclimatic conditions of NFs and NGs were assessed by the magnitude and direction of movements in the Budyko space. The Fu’s equation was applied to estimate changes in basin characteristics. The results indicate that NFs and NGs are likely to experience larger changes in basin characteristics compared to the average of the United States. In general, across the conterminous US, the NFs in mountainous regions are likely to have larger changes in hydroclimatic variables than NFs with lower elevation and NGs. Comparing Forest Service regions, Pacific Northwest, Intermountain, and Northern regions may have a less arid climate with lower freshwater availability. The Southwestern, Northern, Intermountain, and Rocky Mountain regions are likely to experience higher shifts in their basin characteristics. This study can help environmental scientists, and land and water managers improve future land management plans.

The Applicability of SWOT’s Non-Uniform Space–Time Sampling in Hydrologic Model Calibration

The Surface Water and Ocean Topography (SWOT) satellite mission, expected to launch in 2022, will enable near global river discharge estimation from surface water extents and elevations. However, SWOT’s orbit specifications provide non-uniform space–time sampling. Previous studies have demonstrated that SWOT’s unique spatiotemporal sampling has a minimal impact on derived discharge frequency distributions, baseflow magnitudes, and annual discharge characteristics. In this study, we aim to extend the analysis of SWOT’s added value in the context of hydrologic model calibration. We calibrate a hydrologic model using previously derived synthetic SWOT discharges across 39 gauges in the Ohio River Basin. Three discharge timeseries are used for calibration: daily observations, SWOT temporally sampled, and SWOT temporally sampled including estimated uncertainty. Using 10,000 model iterations to explore predefined parameter ranges, each discharge timeseries results in similar optimal model parameters. We find that the annual mean and peak flow values at each gauge location from the optimal parameter sets derived from each discharge timeseries differ by less than 10% percent on average. Our findings suggest that hydrologic models calibrated using discharges derived from SWOT’s non-uniform space–time sampling are likely to achieve results similar to those based on calibrating with in situ daily observations.

PEMIP: Post-fire erosion model inter-comparison project

Land managers often need to predict watershed-scale erosion rates after disturbance or other land cover changes. This study compared commonly used hillslope erosion models to simulate post-fire sediment yields (SY) at both hillslope and watershed scales within the High Park Fire, Colorado, U.S.A. At hillslope scale, simulated SY from four models- RUSLE, AGWA/KINEROS2, WEPP, and a site-specific regression model-were compared to observed SY at 29 hillslopes. At the watershed scale, RUSLE, AGWA/KINEROS2, and WEPP were applied to simulate spatial patterns of SY for two 14-16 km² watersheds using different scales (0.5-25 ha) of hillslope discretization. Simulated spatial patterns were compared between models and to densities of channel heads across the watersheds. Three additional erosion algorithms were implemented within a land surface model to evaluate effects of parameter uncertainty. At the hillslope scale, SY was only significantly correlated to observed SY for the empirical model, but at the watershed scale, sediment loads were significantly correlated to observed channel head densities for all models. Watershed sediment load increased with the size of the hillslope sub-units due to the nonlinear effects of hillslope length on simulated erosion. SY's were closest in magnitude to expected watershed-scale SY when models were divided into the smallest hillslopes. These findings demonstrate that current erosion models are fairly consistent at identifying areas with low and high erosion potential, but the wide range of predicted SY and poor fit to observed SY highlight the need for better field observations and model calibration to obtain more accurate simulations.

Global Reconstruction of Naturalized River Flows at 2.94 Million Reaches

Spatiotemporally continuous global river discharge estimates across the full spectrum of stream orders are vital to a range of hydrologic applications, yet they remain poorly constrained. Here we present a carefully designed modeling effort (Variable Infiltration Capacity land surface model and Routing Application for Parallel computatIon of Discharge river routing model) to estimate global river discharge at very high resolutions. The precipitation forcing is from a recently published 0.1° global product that optimally merged gauge-, reanalysis-, and satellite-based data. To constrain runoff simulations, we use a set of machine learning-derived, global runoff characteristics maps (i.e., runoff at various exceedance probability percentiles) for grid-by-grid model calibration and bias correction. To support spaceborne discharge studies, the river flowlines are defined at their true geometry and location as much as possible-approximately 2.94 million vector flowlines (median length 6.8 km) and unit catchments are derived from a high-accuracy global digital elevation model at 3-arcsec resolution (~90 m), which serves as the underlying hydrography for river routing. Our 35-year daily and monthly model simulations are evaluated against over 14,000 gauges globally. Among them, 35% (64%) have a percentage bias within ±20% (±50%), and 29% (62%) have a monthly Kling-Gupta Efficiency ≥0.6 (0.2), showing data robustness at the scale the model is assessed. This reconstructed discharge record can be used as a priori information for the Surface Water and Ocean Topography satellite mission's discharge product, thus named "Global Reach-level A priori Discharge Estimates for Surface Water and Ocean Topography". It can also be used in other hydrologic applications requiring spatially explicit estimates of global river flows.

A long-term, temporally consistent, gridded daily meteorological dataset for northwestern North America

We describe a spatially contiguous, temporally consistent high-resolution gridded daily meteorological dataset for northwestern North America. This >4 million km² region has high topographic relief, seasonal snowpack, permafrost and glaciers, crosses multiple jurisdictional boundaries and contains the entire Yukon, Mackenzie, Saskatchewan, Fraser and Columbia drainages. We interpolate daily station data to 1/16° spatial resolution using a high-resolution monthly 1971–2000 climatology as a predictor in a thin-plate spline interpolating algorithm. Only temporally consistent climate stations with at least 40 years of record are included. Our approach is designed to produce a dataset well suited for driving hydrological models and training statistical downscaling schemes. We compare our results to two commonly used datasets and show improved performance for climate means, extremes and variability. When used to drive a hydrologic model, our dataset also outperforms these datasets for runoff ratios and streamflow trends in several, high elevation, sub-basins of the Fraser River.

Sensitivity of Evapotranspiration Components in Remote Sensing-Based Models

Accurately estimating evapotranspiration (ET) at large spatial scales is essential to our understanding of land-atmosphere coupling and the surface balance of water and energy. Comparisons between remote sensing-based ET models are difficult due to diversity in model formulation, parametrization and data requirements. The constituent components of ET have been shown to deviate substantially among models as well as between models and field estimates. This study analyses the sensitivity of three global ET remote sensing models in an attempt to isolate the error associated with forcing uncertainty and reveal the underlying variables driving the model components. We examine the transpiration, soil evaporation, interception and total ET estimates of the Penman-Monteith model from the Moderate Resolution Imaging Spectroradiometer (PM-MOD), the Priestley-Taylor Jet Propulsion Laboratory model (PT-JPL) and the Global Land Evaporation Amsterdam Model (GLEAM) at 42 sites where ET components have been measured using field techniques. We analyse the sensitivity of the models based on the uncertainty of the input variables and as a function of the raw value of the variables themselves. We find that, at 10% added uncertainty levels, the total ET estimates from PT-JPL, PM-MOD and GLEAM are most sensitive to Normalized Difference Vegetation Index (NDVI) (%RMSD = 100.0), relative humidity (%RMSD = 122.3) and net radiation (%RMSD = 7.49), respectively. Consistently, systemic bias introduced by forcing uncertainty in the component estimates is mitigated when components are aggregated to a total ET estimate. These results suggest that slight changes to forcing may result in outsized variation in ET partitioning and relatively smaller changes to the total ET estimates. Our results help to explain why model estimates of total ET perform relatively well despite large inter-model divergence in the individual ET component estimates.

A Multialgorithm Approach to Land Surface Modeling of Suspended Sediment in the Colorado Front Range

A new paradigm of simulating suspended sediment load (SSL) with a Land Surface Model (LSM) is presented here. Five erosion and SSL algorithms were applied within a common LSM framework to quantify uncertainties and evaluate predictability in two steep, forested catchments (>1,000 km²). The algorithms were chosen from among widely used sediment models, including empirically based: monovariate rating curve (MRC) and the Modified Universal Soil Loss Equation (MUSLE); stochastically based: the Load Estimator (LOADEST); conceptually based: the Hydrologic Simulation Program-Fortran (HSPF); and physically based: the Distributed Hydrology Soil Vegetation Model (DHSVM). The algorithms were driven by the hydrologic fluxes and meteorological inputs generated from the Variable Infiltration Capacity (VIC) LSM. A multiobjective calibration was applied to each algorithm and optimized parameter sets were validated over an excluded period, as well as in a transfer experiment to a nearby catchment to explore parameter robustness. Algorithm performance showed consistent decreases when parameter sets were applied to periods with greatly differing SSL variability relative to the calibration period. Of interest was a joint calibration of all sediment algorithm and streamflow parameters simultaneously, from which trade-offs between streamflow performance and partitioning of runoff and base flow to optimize SSL timing were noted, decreasing the flexibility and robustness of the streamflow to adapt to different time periods. Parameter transferability to another catchment was most successful in more process-oriented algorithms, the HSPF and the DHSVM. This first-of-its-kind multialgorithm sediment scheme offers a unique capability to portray acute episodic loading while quantifying trade-offs and uncertainties across a range of algorithm structures.