Bayesian Analysis of the Impact of Rainfall Data Product on Simulated Slope Failure for North Carolina Locations

In the past decades, many different approaches have been developed in the literature to quantify the load-carrying capacity and geotechnical stability (or the factor of safety, Fs) of variably saturated hillslopes. Much of this work has focused on a deterministic characterization of hillslope stability. Yet, simulated Fs values are subject to considerable uncertainty due to our inability to characterize accurately the soil mantles properties (hydraulic, geotechnical, and geomorphologic) and spatiotemporal variability of the moisture content of the hillslope interior. This is particularly true at larger spatial scales. Thus, uncertainty-incorporating analyses of physically based models of rain-induced landslides are rare in the literature. Such landslide modeling is typically conducted at the hillslope scale using gauge-based rainfall forcing data with rather poor spatiotemporal coverage. For regional landslide modeling, the specific advantages and/or disadvantages of gauge-only, radar-merged and satellite-based rainfall products are not clearly established. Here, we compare and evaluate the performance of the Transient Rainfall Infiltration and Grid-based Regional Slope-stability analysis (TRIGRS) model for three different rainfall products using 112 observed landslides in the period between 2004 and 2011 from the North Carolina Geological Survey database. Our study includes the Tropical Rainfall Measuring Mission (TRMM) Multi-satellite Precipitation Analysis Version 7 (TMPA V7), the North American Land Data Assimilation System Phase 2 (NLDAS-2) analysis, and the reference truth Stage IV precipitation. TRIGRS model performance was rather inferior with the use of literature values of the geotechnical parameters and soil hydraulic properties from ROSETTA using soil textural and bulk density data from SSURGO (Soil Survey Geographic database). The performance of TRIGRS improved considerably after Bayesian estimation of the parameters with the DiffeRential Evolution Adaptive Metropolis (DREAM) algorithm using Stage IV precipitation data. Hereto, we use a likelihood function that combines binary slope failure information from landslide event and null periods using multivariate frequency distribution-based metrics such as the false discovery and false omission rates. Our results demonstrate that the Stage IV-inferred TRIGRS parameter distributions generalize well to TMPA and NLDAS-2 precipitation data, particularly at sites with considerably larger TMPA and NLDAS-2 rainfall amounts during landslide events than null periods. TRIGRS model performance is then rather similar for all three rainfall products. At higher elevations, however, the TMPA and NLDAS-2 precipitation volumes are insufficient and their performance with the Stage IV-derived parameter distributions indicates their inability to accurately characterize hillslope stability

Quantifying Parameter Sensitivity, Interaction and Transferability in Hydrologically Enhanced Versions of Noah-LSM over Transition Zones

We use sensitivity analysis to identify the parameters that are most responsible for shaping land surface model (LSM) simulations and to understand the complex interactions in three versions of the Noah LSM: the standard version (STD), a version enhanced with a simple groundwater module (GW), and version augmented by a dynamic phenology module (DV). We use warm season, high-frequency, near-surface states and turbulent fluxes collected over nine sites in the US Southern Great Plains. We quantify changes in the pattern of sensitive parameters, the amount and nature of the interaction between parameters, and the covariance structure of the distribution of behavioral parameter sets. Using Sobol s total and first-order sensitivity indexes, we show that very few parameters directly control the variance of the model output. Significant parameter interaction occurs so that not only the optimal parameter values differ between models, but the relationships between parameters change. GW decreases parameter interaction and appears to improve model realism, especially at wetter sites. DV increases parameter interaction and decreases identifiability, implying it is overparameterized and/or underconstrained. A case study at a wet site shows GW has two functional modes: one that mimics STD and a second in which GW improves model function by decoupling direct evaporation and baseflow. Unsupervised classification of the posterior distributions of behavioral parameter sets cannot group similar sites based solely on soil or vegetation type, helping to explain why transferability between sites and models is not straightforward. This evidence suggests a priori assignment of parameters should also consider climatic differences

Comment on 'Shang S. 2012. Calculating actual crop evapotranspiration under soil water stress conditions with appropriate numerical methods and time step. Hydrological Processes 26: 3338-3343. DOI: 10.1002/hyp.8405'

A previous study analyzed errors in the numerical calculation of actual crop evapotranspiration (ET(sub a)) under soil water stress. Assuming no irrigation or precipitation, it constructed equations for ET(sub a) over limited soil-water ranges in a root zone drying out due to evapotranspiration. It then used a single crop-soil composite to provide recommendations about the appropriate usage of numerical methods under different values of the time step and the maximum crop evapotranspiration (ET(sub c)). This comment reformulates those ET(sub a) equations for applicability over the full range of soil water values, revealing a dependence of the relative error in numerical ET(sub a) on the initial soil water that was not seen in the previous study. It is shown that the recommendations based on a single crop-soil composite can be invalid for other crop-soil composites. Finally, a consideration of the numerical error in the time-cumulative value of ET(sub a) is discussed besides the existing consideration of that error over individual time steps as done in the previous study. This cumulative ET(sub a) is more relevant to the final crop yield

Parameter Sensitivity of the Noah-MP Land Surface Model with Dynamic Vegetation

The Noah land surface model with multiple parameterization options (Noah-MP) includes a routine for dynamic simulation of vegetation carbon assimilation and soil carbon decomposition processes. To use remote sensing observations of vegetation to constrain simulations from this model, it is necessary first to understand the sensitivity of the model to its parameters. This is required for efficient parameter estimation, which is both a valuable way to use observations and also a first or concurrent step in many state-updating data assimilation procedures. We use variance decomposition to assess the sensitivity of estimates of sensible heat, latent heat, soil moisture, and net ecosystem exchange made by certain standard Noah-MP configurations that include dynamic simulation of vegetation and carbon to forty-three primary user-specified parameters. This is done using thirty-two years' worth of data from ten international FluxNet sites. Findings indicate that there are five soil parameters and six (or more) vegetation parameters (depending on the model configuration) that act as primary controls on these states and fluxes

A Comparison of Methods for a Priori Bias Correction in Soil Moisture Data Assimilation

Data assimilation is being increasingly used to merge remotely sensed land surface variables such as soil moisture, snow and skin temperature with estimates from land models. Its success, however, depends on unbiased model predictions and unbiased observations. Here, a suite of continental-scale, synthetic soil moisture assimilation experiments is used to compare two approaches that address typical biases in soil moisture prior to data assimilation: (i) parameter estimation to calibrate the land model to the climatology of the soil moisture observations, and (ii) scaling of the observations to the model s soil moisture climatology. To enable this research, an optimization infrastructure was added to the NASA Land Information System (LIS) that includes gradient-based optimization methods and global, heuristic search algorithms. The land model calibration eliminates the bias but does not necessarily result in more realistic model parameters. Nevertheless, the experiments confirm that model calibration yields assimilation estimates of surface and root zone soil moisture that are as skillful as those obtained through scaling of the observations to the model s climatology. Analysis of innovation diagnostics underlines the importance of addressing bias in soil moisture assimilation and confirms that both approaches adequately address the issue

KINEROS2-AGWA: Model Use, Calibration, and Validation

KINEROS (KINematic runoff and EROSion) originated in the 1960s as a distributed event-based model that conceptualizes a watershed as a cascade of overland flow model elements that flow into trapezoidal channel model elements. KINEROS was one of the first widely available watershed models that interactively coupled a finite difference approximation of the kinematic overland flow equations to a physically based infiltration model. Development and improvement of KINEROS continued from the 1960s on a variety of projects for a range of purposes, which has resulted in a suite of KINEROS-based modeling tools. This article focuses on KINEROS2 (K2), a spatially distributed, event-based watershed rainfall-runoff and erosion model, and the companion ArcGIS-based Automated Geospatial Watershed Assessment (AGWA) tool. AGWA automates the time-consuming tasks of watershed delineation into distributed model elements and initial parameterization of these elements using commonly available, national GIS data layers. A variety of approaches have been used to calibrate and validate K2 successfully across a relatively broad range of applications (e.g., urbanization, pre- and post-fire, hillslope erosion, erosion from roads, runoff and recharge, and manure transport). The case studies presented in this article (1) compare lumped to stepwise calibration and validation of runoff and sediment at plot, hillslope, and small watershed scales; and (2) demonstrate an uncalibrated application to address relative change in watershed response to wildfire

Modelling the impacts of agricultural management practices on river water quality in Eastern England

Agricultural diffuse water pollution remains a notable global pressure on water quality, posing risks to aquatic ecosystems, human health and water resources and as a result legislation has been introduced in many parts of the world to protect water bodies. Due to their efficiency and cost-effectiveness, water quality models have been increasingly applied to catchments as Decision Support Tools (DSTs) to identify mitigation options that can be introduced to reduce agricultural diffuse water pollution and improve water quality. In this study, the Soil and Water Assessment Tool (SWAT) was applied to the River Wensum catchment in eastern England with the aim of quantifying the long-term impacts of potential changes to agricultural management practices on river water quality. Calibration and validation were successfully performed at a daily time-step against observations of discharge, nitrate and total phosphorus obtained from high-frequency water quality monitoring within the Blackwater sub-catchment, covering an area of 19.6 km2. A variety of mitigation options were identified and modelled, both singly and in combination, and their long-term effects on nitrate and total phosphorus losses were quantified together with the 95% uncertainty range of model predictions. Results showed that introducing a red clover cover crop to the crop rotation scheme applied within the catchment reduced nitrate losses by 19.6%. Buffer strips of 2 m and 6 m width represented the most effective options to reduce total phosphorus losses, achieving reductions of 12.2% and 16.9%, respectively. This is one of the first studies to quantify the impacts of agricultural mitigation options on long-term water quality for nitrate and total phosphorus at a daily resolution, in addition to providing an estimate of the uncertainties of those impacts. The results highlighted the need to consider multiple pollutants, the degree of uncertainty associated with model predictions and the risk of unintended pollutant impacts when evaluating the effectiveness of mitigation options, and showed that high-frequency water quality datasets can be applied to robustly calibrate water quality models, creating DSTs that are more effective and reliable

Bias adjustment of satellite rainfall data through stochastic modeling: Methods development and application to Nepal

Estimating precipitation over large spatial areas remains a challenging problem for hydrologists. Sparse ground-based gauge networks do not provide a robust basis for interpolation, and the reliability of remote sensing products, although improving, is still imperfect. Current techniques to estimate precipitation rely on combining these different kinds of measurements to correct the bias in the satellite observations. We propose a novel procedure that, unlike existing techniques, (i) allows correcting the possibly confounding effects of different sources of errors in satellite estimates, (ii) explicitly accounts for the spatial heterogeneity of the biases and (iii) allows the use of non overlapping historical observations. The proposed method spatially aggregates and interpolates gauge data at the satellite grid resolution by focusing on parameters that describe the frequency and intensity of the rainfall observed at the gauges. The resulting gridded parameters can then be used to adjust the probability density function of satellite rainfall observations at each grid cell, accounting for spatial heterogeneity. Unlike alternate methods, we explicitly adjust biases on rainfall frequency in addition to its intensity. Adjusted rainfall distributions can then readily be applied as input in stochastic rainfall generators or frequency domain hydrological models. Finally, we also provide a procedure to use them to correct remotely sensed rainfall time series.We apply the method to adjust the distributions of daily rainfall observed by the TRMM satellite in Nepal, which exemplifies the challenges associated with a sparse gauge network and large biases due to complex topography. In a cross-validation analysis on daily rainfall from TRMM 3B42 v6, we find that using a small subset of the available gauges, the proposed method outperforms local rainfall estimations using the complete network of available gauges to directly interpolate local rainfall or correct TRMM by adjusting monthly means. We conclude that the proposed frequency-domain bias correction approach is robust and reliable compared to other bias correction approaches. © 2013 Elsevier Ltd

Flash Flood Forecasting for the Semi-Arid Southwestern United States

Flash flooding in the semi-arid United States poses a significant danger to life and property. One effective way to mitigate flood risk is by implementing a rainfall-runoff model in a real-time forecast and warning system. This study investigated the feasibility of using the mechanistic, distributed semi-arid rainfall-runoff model KINEROS2 driven by high resolution radar rainfall input estimates obtained from the NEXRAD WSR-88D DHR reflectivity measurements in such a system. The original procedural paradigm-based KINEROS2 Fortran 77 code with space-time looping was recoded into an object-oriented Fortran 90 code with time-space looping for this purpose. The recoded form is now applicable to large basins, is easily future-extensible, and individual modules can be incorporated into other models.Sources of operational uncertainty in the above system were investigated for their influence over several events within a sub-basin of the USDA-ARS Walnut Gulch Experimental Watershed. Uncertainties considered were in the rainfall estimates, the model parameters, and the initial conditions. The variance-based Sobol' method of global sensitivity analysis conditioned on the observed streamflow showed that the uncertainty in the modeled response was heavily dominated by the operational variability of biases in the radar rainfall depth estimates. Sensitivities to KINEROS2 parameters indicates the need for improved representation of semi-arid hillslope hydrology in small basins, while pointing to specific influential, but poorly identified model parameters towards which field investigations should be directed. The significant influence of initial hillslope soil moisture showed the requirement of a sophisticated inter-storm model component for a continuous forecasting model.A synthetic study data was used to further explore the phenomena seen in the above real data study, of behavioral modifier set inconsistency across all events and of irreducibility in the spatial modifier ranges. The former was found to be attributable to wide uncertainty ranges in the sources of uncertainty, and the latter to the high distributed model non-linearity with associated interactions. These contribute towards a high predictive uncertainty in operational forecasting.Overall, the GLUE-based predictive uncertainty method with behavioral classification and accommodation of wide operational source uncertainty ranges is recommended as a simple and effective setup for operational flash flood forecasting