scholarly journals From maps to movies: high-resolution time-varying sensitivity analysis for spatially distributed watershed models

2013 ◽  
Vol 17 (12) ◽  
pp. 5109-5125 ◽  
Author(s):  
J. D. Herman ◽  
J. B. Kollat ◽  
P. M. Reed ◽  
T. Wagener

Abstract. Distributed watershed models are now widely used in practice to simulate runoff responses at high spatial and temporal resolutions. Counter to this purpose, diagnostic analyses of distributed models currently aggregate performance measures in space and/or time and are thus disconnected from the models' operational and scientific goals. To address this disconnect, this study contributes a novel approach for computing and visualizing time-varying global sensitivity indices for spatially distributed model parameters. The high-resolution model diagnostics employ the method of Morris to identify evolving patterns in dominant model processes at sub-daily timescales over a six-month period. The method is demonstrated on the United States National Weather Service's Hydrology Laboratory Research Distributed Hydrologic Model (HL-RDHM) in the Blue River watershed, Oklahoma, USA. Three hydrologic events are selected from within the six-month period to investigate the patterns in spatiotemporal sensitivities that emerge as a function of forcing patterns as well as wet-to-dry transitions. Events with similar magnitudes and durations exhibit significantly different performance controls in space and time, indicating that the diagnostic inferences drawn from representative events will be heavily biased by the a priori selection of those events. By contrast, this study demonstrates high-resolution time-varying sensitivity analysis, requiring no assumptions regarding representative events and allowing modelers to identify transitions between sets of dominant parameters or processes a posteriori. The proposed approach details the dynamics of parameter sensitivity in nearly continuous time, providing critical diagnostic insights into the underlying model processes driving predictions. Furthermore, the approach offers the potential to identify transition points between dominant parameters and processes in the absence of observations, such as under nonstationarity.

2013 ◽  
Vol 10 (8) ◽  
pp. 10775-10808
Author(s):  
J. D. Herman ◽  
J. B. Kollat ◽  
P. M. Reed ◽  
T. Wagener

Abstract. Distributed watershed models are now widely used in practice to simulate runoff responses at high spatial and temporal resolutions. Counter to this purpose, diagnostic analyses of distributed models currently aggregate performance measures in space and/or time and are thus disconnected from the models' operational and scientific goals. To address this disconnect, this study contributes a novel approach for computing and visualizing time-varying global sensitivity indices for spatially distributed model parameters. The high-resolution model diagnostics employ the method of Morris to identify evolving patterns in dominant model processes at sub-daily timescales over a six-month period. The method is demonstrated on the United States National Weather Service's Hydrology Laboratory Research Distributed Hydrologic Model (HL-RDHM) in the Blue River watershed, Oklahoma, USA. Three hydrologic events are selected from within the six-month period to investigate the patterns in spatiotemporal sensitivities that emerge as a function of forcing patterns as well as wet-to-dry transitions. Surprisingly, events with similar magnitudes and durations exhibit significantly different performance controls in space and time, indicating that the diagnostic inferences drawn from representative events will be heavily biased by the a priori selection of those events. By contrast, this study demonstrates high-resolution time-varying sensitivity analysis, requiring no assumptions regarding representative events and allowing modelers to identify transitions between modeled hydrologic regimes a posteriori. The proposed approach details the dynamics of parameter sensitivity in nearly continuous time, providing critical diagnostic insights into the underlying model processes driving predictions. Furthermore, the approach offers the potential to identify transition points between hydrologic regimes under nonstationarity.


2013 ◽  
Vol 10 (4) ◽  
pp. 4275-4299 ◽  
Author(s):  
J. D. Herman ◽  
J. B. Kollat ◽  
P. M. Reed ◽  
T. Wagener

Abstract. The increase in spatially distributed hydrologic modeling warrants a corresponding increase in diagnostic methods capable of analyzing complex models with large numbers of parameters. Sobol' sensitivity analysis has proven to be a valuable tool for diagnostic analyses of hydrologic models. However, for many spatially distributed models, the Sobol' method requires a prohibitive number of model evaluations to reliably decompose output variance across the full set of parameters. We investigate the potential of the method of Morris, a screening-based sensitivity approach, to provide results sufficiently similar to those of the Sobol' method at a greatly reduced computational expense. The methods are benchmarked on the Hydrology Laboratory Research Distributed Hydrologic Model (HL-RDHM) model over a six-month period in the Blue River Watershed, Oklahoma, USA. The Sobol' method required over six million model evaluations to ensure reliable sensitivity indices, corresponding to more than 30 000 computing hours and roughly 180 gigabytes of storage space. We find that the method of Morris is able to correctly identify sensitive and insensitive parameters with 300 times fewer model evaluations, requiring only 100 computing hours and 1 gigabyte of storage space. Method of Morris proves to be a promising diagnostic approach for global sensitivity analysis of highly parameterized, spatially distributed hydrologic models.


2013 ◽  
Vol 17 (7) ◽  
pp. 2893-2903 ◽  
Author(s):  
J. D. Herman ◽  
J. B. Kollat ◽  
P. M. Reed ◽  
T. Wagener

Abstract. The increase in spatially distributed hydrologic modeling warrants a corresponding increase in diagnostic methods capable of analyzing complex models with large numbers of parameters. Sobol' sensitivity analysis has proven to be a valuable tool for diagnostic analyses of hydrologic models. However, for many spatially distributed models, the Sobol' method requires a prohibitive number of model evaluations to reliably decompose output variance across the full set of parameters. We investigate the potential of the method of Morris, a screening-based sensitivity approach, to provide results sufficiently similar to those of the Sobol' method at a greatly reduced computational expense. The methods are benchmarked on the Hydrology Laboratory Research Distributed Hydrologic Model (HL-RDHM) over a six-month period in the Blue River watershed, Oklahoma, USA. The Sobol' method required over six million model evaluations to ensure reliable sensitivity indices, corresponding to more than 30 000 computing hours and roughly 180 gigabytes of storage space. We find that the method of Morris is able to correctly screen the most and least sensitive parameters with 300 times fewer model evaluations, requiring only 100 computing hours and 1 gigabyte of storage space. The method of Morris proves to be a promising diagnostic approach for global sensitivity analysis of highly parameterized, spatially distributed hydrologic models.


2010 ◽  
Vol 11 (3) ◽  
pp. 781-796 ◽  
Author(s):  
Jonathan J. Gourley ◽  
Scott E. Giangrande ◽  
Yang Hong ◽  
Zachary L. Flamig ◽  
Terry Schuur ◽  
...  

Abstract Rainfall estimated from the polarimetric prototype of the Weather Surveillance Radar-1988 Doppler [WSR-88D (KOUN)] was evaluated using a dense Micronet rain gauge network for nine events on the Ft. Cobb research watershed in Oklahoma. The operation of KOUN and its upgrade to dual polarization was completed by the National Severe Storms Laboratory. Storm events included an extreme rainfall case from Tropical Storm Erin that had a 100-yr return interval. Comparisons with collocated Micronet rain gauge measurements indicated all six rainfall algorithms that used polarimetric observations had lower root-mean-squared errors and higher Pearson correlation coefficients than the conventional algorithm that used reflectivity factor alone when considering all events combined. The reflectivity based relation R(Z) was the least biased with an event-combined normalized bias of −9%. The bias for R(Z), however, was found to vary significantly from case to case and as a function of rainfall intensity. This variability was attributed to different drop size distributions (DSDs) and the presence of hail. The synthetic polarimetric algorithm R(syn) had a large normalized bias of −31%, but this bias was found to be stationary. To evaluate whether polarimetric radar observations improve discharge simulation, recent advances in Markov Chain Monte Carlo simulation using the Hydrology Laboratory Research Distributed Hydrologic Model (HL-RDHM) were used. This Bayesian approach infers the posterior probability density function of model parameters and output predictions, which allows us to quantify HL-RDHM uncertainty. Hydrologic simulations were compared to observed streamflow and also to simulations forced by rain gauge inputs. The hydrologic evaluation indicated that all polarimetric rainfall estimators outperformed the conventional R(Z) algorithm, but only after their long-term biases were identified and corrected.


Author(s):  
Mehmet Cüneyd Demirel ◽  
Julian Koch ◽  
Gorka Mendiguren ◽  
Simon Stisen

Hydrologic models are conventionally constrained and evaluated using point measurements of streamflow, which represents an aggregated catchment measure. As a consequence of this single objective focus, model parametrization and model parameter sensitivity are typically not reflecting other aspects of catchment behavior. Specifically for distributed models, the spatial pattern aspect is often overlooked. Our paper examines the utility of multiple performance measures in a spatial sensitivity analysis framework to determine the key parameters governing the spatial variability of predicted actual evapotranspiration (AET). Latin hypercube one-at-a-time (LHS-OAT) sampling strategy with multiple initial parameter sets was applied using the mesoscale hydrologic model (mHM) and a total of 17 model parameters were identified as sensitive. The results indicate different parameter sensitivities for different performance measures focusing on temporal hydrograph dynamics and spatial variability of actual evapotranspiration. While spatial patterns were found to be sensitive to vegetation parameters, streamflow dynamics were sensitive to pedo-transfer function (PTF) parameters. Above all, our results show that behavioral model definition based only on streamflow metrics in the generalized likelihood uncertainty estimation (GLUE) type methods require reformulation by incorporating spatial patterns into the definition of threshold values to reveal robust hydrologic behavior in the analysis.


2011 ◽  
Vol 8 (6) ◽  
pp. 10739-10780
Author(s):  
V. Ruiz-Villanueva ◽  
M. Borga ◽  
D. Zoccatelli ◽  
L. Marchi ◽  
E. Gaume ◽  
...  

Abstract. The 2 June 2008 flood-producing storm on the Starzel river basin in South-West Germany is examined as a prototype for organized convective systems that dominate the upper tail of the precipitation frequency distribution and are likely responsible for the flash flood peaks in this region. The availability of high-resolution rainfall estimates from radar observations and a rain gauge network, together with indirect peak discharge estimates from a detailed post-event survey, provides the opportunity to study the hydrometeorological and hydrological mechanisms associated with this extreme storm and the ensuing flood. Radar-derived rainfall, streamgauge data and indirect estimates of peak discharges are used along with a distributed hydrologic model to reconstruct hydrographs at multiple locations. The influence of storm structure, evolution and motion on the modeled flood hydrograph is examined by using the "spatial moments of catchment rainfall" (Zoccatelli et al., 2011). It is shown that downbasin storm motion had a noticeable impact on flood peak magnitude. Small runoff ratios (less than 20%) characterized the runoff response. The flood response can be reasonably well reproduced with the distributed hydrological model, using high resolution rainfall observations and model parameters calibrated at a river section which includes most of the area impacted by the storm.


2020 ◽  
Vol 102 (4) ◽  
pp. 690-704 ◽  
Author(s):  
Pascal Paul

This paper studies how monetary policy jointly affects asset prices and the real economy in the United States. I develop an estimator that uses high-frequency surprises as a proxy for the structural monetary policy shocks. This is achieved by integrating the surprises into a vector autoregressive model as an exogenous variable. I use current short-term rate surprises because these are least affected by an information effect. When allowing for time-varying model parameters, I find that compared to the response of output, the reaction of stock and house prices to monetary policy shocks was particularly low before the 2007–2009 financial crisis.


2020 ◽  
Author(s):  
Bibi S Naz ◽  
Wendy Sharples ◽  
Klaus Goergen ◽  
Stefan Kollet

<p> <span>High-resolution large-scale predictions of hydrologic states and fluxes are important for many regional-scale applications and water resource management. However, because of uncertainties related to forcing data, model structural errors arising from simplified representations of hydrological processes or uncertain model parameters, model simulations remain uncertain. To quantify this uncertainty, multi-model simulations were performed at 3km resolution over the European continent using the Community Land Model (CLM3.5) and the ParFlow hydrologic model. While Parflow uses a similar approach as CLM in simulating the snow, vegetation and land-atmosphere exchange processes, it simulates three-dimensional variably saturated groundwater flow solving Richards equation and overland flow with a two-dimensional kinematic wave approximation. </span><span>The </span><span>CLM</span><span>3.5</span><span> uses a simple groundwater model to account for groundwater recharge and discharge processes. Both models were driven with the COSMO-REA6 reanalysis dataset at 6km resolution for the time period from 2000 to 2006 at an hourly time step, and both used the same datasets for the static input variables (such as topography, vegetation and soil properties). The performance of both models was analyzed through comparisons with independent observations including satellite-derived and in-situ soil moisture, evapotranspiration, river discharge, water table depth and total water storage datasets. Overall, both models capture the interannual variability in the hydrologic states and fluxes well, however differences in performance between models showed the uncertainty associated with the representation of hydrological processes, such as groundwater flow and soil moisture and its control on latent and sensible heat fluxes at the surface.</span></p>


Sign in / Sign up

Export Citation Format

Share Document