Professional Documents
Culture Documents
Bayesian Approach For Uncertainty Analysis of An Urban Storm Wate
Bayesian Approach For Uncertainty Analysis of An Urban Storm Wate
Misgana K. Muleta, M.ASCE 1; Jonathan McMillan 2; Geremew G. Amenu 3; and Steven J. Burian, M.ASCE 4
Abstract: The significance of uncertainty analysis (UA) to quantify reliability of model simulations is being recognized. Consequently,
literature on parameter and predictive uncertainty assessment of water resources models has been rising. Applications dealing with urban
drainage systems are, however, very limited. This study applies formal Bayesian approach for uncertainty analysis of a widely used storm
water management model and illustrates the methodology using a highly urbanized watershed in the Los Angeles Basin, California. A flexible
likelihood function that accommodates heteroscedasticity, non-normality, and temporal correlation of model residuals has been used for the
study along with a Markov-chain Monte Carlo-based sampling scheme. The solution of the UA model has been compared with the solution of
the conventional calibration methodology widely practiced in water resources modeling. Results indicate that the maximum likelihood sol
ution determined using the UA model produced runoff simulations that are of comparable accuracy with the solution of the traditional
calibration method while also accurately characterizing structure of the model residuals. The UA model also successfully determined both
parameter uncertainty and total predictive uncertainty for the watershed. Contribution of parameter uncertainty to total predictive uncertainty
was found insignificant for the study watershed, underlying the importance of other sources of uncertainty, including data and model struc
ture. Overall, the UA methodology proved promising for sensitivity analysis, calibration, parameter uncertainty, and total predictive uncer
tainty analysis of urban storm water management models. DOI: 10.1061/(ASCE)HE.1943-5584.0000705. © 2013 American Society of Civil
Engineers.
CE Database subject headings: Uncertainty principles; Calibration; Bayesian analysis; Stormwater Management; Parameters; Markov
process; Monte Carlo method; Urban areas.
Author keywords: Uncertainty analysis; Calibration; Bayesian approach; Storm water management; Parameter uncertainty; Predictive
uncertainty; Markov-chain Monte Carlo.
Likewise, continuous (long-term) simulation was not consid where Y = model simulated output; O = observed runoff; Omean =
ered because of computational concern. Short-term simulation mean of the observations, which the NSE uses as a benchmark
Table 1. Model Parameters and Ranges Used for Calibration and data, and calibration data uncertainty. Total predictive uncertainty
Uncertainty Analysis was calculated using the methodology described in Schoups
Name Description Minimum Maximum and Vrugt (2010) based on the error model parameters determined
by GL-DREAMðZSÞ. The ML parameter values determined by
SWMM parameters Percent adjustment
Width Subcatchment width (m) −90 150
GL-DREAMðZSÞ are benchmarked against the calibration results
Slope Subcatchment slope (%) −20 20 obtained using DDS and are also compared in terms of their ability
%Imperv Percentage of impervious −15 15 to fit different parts of the hydrograph.
area (%)
N-Imperv Manning n for impervious −50 100
area Results and Discussion
N-Perv Manning n for pervious area −50 100
Dstore-Imperv Depression storage for −95 100
impervious area (mm) Single-Objective Calibration
Dstore-Perv Depression storage for −95 300
Results of the single-objection calibration are summarized in Fig. 2
pervious area (mm)
%Zero-Imperv Percent of the impervious −100 100
and Tables 2 and 3. Fig. 2 compares the streamflow simulated using
area with no depression the optimal parameter values identified by DDS with the stream-
storage (%) flow observed at Swatelle station for both calibration [Figs. 2(a–d)]
Curve number NRCS runoff curve number −20 20 and verification [Figs. 2(e–h)] periods. Performance of the cali
Drying Time Time for a fully saturated soil −95 100 brated model was also tested using the efficiency criteria given
to completely dry (days) in Table 2. The graphical comparison and the efficiency criteria
Conduit n Manning’s roughness −50 100 indicate that the single-objection calibration performed very well
coefficient for conduit for both calibration and verification periods. However, as shown
Error parameters Parameter values in Fig. 2, a closer look into the characteristics of the residuals
σ0 Heteroscedasticity intercept 0 100 (i.e., the difference between the observed streamflow and the
(m3 =s) streamflow simulated by the calibrated model) depicts that the as
σ1 Heteroscedasticity slope 0 1 sumptions (i.e., homoscedasticity, Gaussian distribution, and tem
ϕ1 Lag one autocorrelation −1 1 poral independence) made by the objective functions almost always
coefficient used in model calibrations, including the MAE used for this study,
ϕ2 Lag two autocorrelation −1 1
are unjustified. The residuals exhibited heteroscedasticity (i.e., they
coefficient
ϕ3 Lag three autocorrelation −1 1
increase with the magnitude of streamflow) as shown in Fig. 2(b),
coefficient they do not follow Gaussian distribution [Fig. 2(c)], and they are
β Kurtosis parameter −1 1 temporally correlated [Fig. 2(d)] for both calibration and verifica
tion periods. Similar findings have been reported by other studies
including those by Schoups and Vrugt (2010) and regarding char
against which performance of the model simulations is evaluated; acteristics of the residuals generated from solutions of the tradi
and N = number of data points (observations). tional calibration methods.
In addition to relying on unrealistic assumptions on residuals,
the conventional calibration models attempt to identify a single best
Parameter Posteriors and Predictive Uncertainty solution based on the assumption that data, model structure, and
Posterior distributions of the 11 runoff parameters were estimated model parameters are all error-free. It is now a common knowledge
with DREAMðZSÞ and the GL function. Six additional error model that input and output data are subjected to substantial uncertainty;
parameters were considered for the GL. These include σ0 and σ1 in no model structure is a true representation of the watershed being
Eq. (7), to explicitly account for heteroscedasticity of model resid studied, and optimal parameter sets are not unique for a given
uals, the kurtosis parameter, β, in Eq. (5) to account for non- watershed (i.e., multiple set of parameters can be equally good for
normality of the residuals, and ϕ1 , ϕ2 , and ϕ3 in Eq. (8) to allow the watershed). As such, these calibration methods provide no in
for autocorrelation of the residuals. It was initially assumed formation on reliability of the optimal solution. The uncertainty
that the residual distribution is not skewed, setting skewness analysis model described in this study has been developed to
parameter, ξ, equal to 1, and then this assumption was checked a address these vital limitations.
posteriori. A total of 100,000 SWMM5 simulations were used to
sample the posterior distribution of the parameters. Convergence of
Parameter Uncertainty
GL-DREAMðZSÞ to a stable posterior PDF was monitored using the
R statistic of Gelman and Rubin (1992). Convergence is declared As previously described, 100,000 SWMM5 simulations were run
when Rj ≤ 1.2 for all j ¼ 1; : : : ; d, where d is the number of model for the GL-DREAMðZSÞ UA model. Fig. 3 shows progress of the R
parameters being analyzed (i.e., 17 in this study). The last 5,000 statistic of Gelman and Rubin (1992) that has been used to test
GL-DREAMðZSÞ runs that meet the convergence criteria were ex convergence of the UA runs. The plot indicates that the 100,000
tracted, and parameter posteriors were determined and reported for model runs used for the analysis were sufficient to meet the
each individual parameter. convergence criteria of R ≤ 1.2 for all the SWMM5 and error
Once the posterior distribution of the model parameters is model parameters considered for the study. Fig. 4 shows the pos
known, runoff predictive uncertainty can be estimated by propa terior histograms obtained for the parameters using the last 3,000
gating the different samples of the posterior distribution through simulations. As shown in Fig. 3, the convergence criteria was met
the SWMM5 model and reporting the respective prediction uncer after approximately 62,000 model simulations, implying that the
tainty ranges (e.g., 95% confidence interval). This prediction last 3,000 simulations used to generate the posterior histograms
interval, however, represents parameter uncertainty only; it doesn’t have met the convergence criteria, and each of these parameter sets
consider other sources of error, including model structure, forcing represents a reasonable SWMM5 model for the watershed.
500
(a)
Streamflow [m /s]
400
3
300
200
100
0
0 100 200 300 400 500 600
Time Count [15 minute interval]
40 1 1
Partial autocorrelation
(b) (c) (d)
0.8
Residuals (a)
20
0.5
Density
0.6
0
0.4
0
-20 0.2
-40 0 -0.5
0 200 400 -2 0 2 0 10 20 30 40 50
3 Residuals (a) Lag [15 minute interval]
Simulated flow [m /s]
400
(e)
Streamflow [m /s]
3
300
200
100
0
0 50 100 150 200 250 300 350 400 450 500 550
Time Count [15 minute interval]
50 1 1
Partial autocorrelation
(f) (g) (h)
0.8
Residuals (a)
25
0.5
Density
0.6
0
0.4
0
-25 0.2
-50 0 -0.5
0 100 200 300 -2 0 2 0 10 20 30 40 50
3 Residuals (a) Lag [15 minute interval]
Simulated flow [m /s]
Fig. 2. Comparison of observed and simulated streamflow and diagnostic plots of the residuals obtained using DDS for both calibration (a–d) and
verification (e–h) periods: (a and e) comparison between observed (dots) and simulated (solid line) streamflow; (b and f) illustration of heterosce
dasticity of the residuals; (c and g) comparison of observed PDF of the residuals to normal distribution; (d and h) illustration of autocorrelation of the
residuals; dashed lines in (d) and (h) show lower and upper bounds of the 95% confidence interval
Fig. 4(a) offers important information regarding the relative sets did not show significant variability as shown in Fig. 5.
importance of the SWMM5 parameters considered for the analysis. This suggests that the SWMM5 parameters that exhibited wide
Except for percentage imperviousness (Imperv), depression storage uncertainty range do not have substantial effect on rainfall-runoff
for impervious subareas (Dstore-Imperv), and percentage of the im characteristics of the Ballona Creek watershed. This has practical
pervious subarea with no depression storage (% Zero Imperv), the implication, for example, in terms of prioritizing resources on data
uncertainty bound of the other SWMM5 parameters is very wide. collection. Availability of more accurate data that characterize the
Nonetheless, the streamflow simulated by the last 3,000 parameter insensitive parameters may not help in improving accuracy of run
off simulation for the watershed. On the contrary, uncertainty of the
Table 2. Efficiency Criteria Values Obtained Using DDS and model predictions can be further reduced from having more reliable
GL-DREAMðZSÞ imperviousness data. Unlike the SWMM5 parameters, the error
Volume (mm) model parameters given in Fig. 4(b) produced narrow uncertainty
MAE Bias range indicating their identifiablity.
Method Period (m3 =s) NSE (%) Observed Simulated Tables 2 and 3 compare the solution of the GL-DREAMðZSÞ
DDS C 6.89 0.94 11.15 59.6 66.2 model with the single-objection calibration results. Performance
V 9.79 0.90 −18.82 83.6 68.9 of the maximum likelihood solution has been summarized using
GL-DREAMðZSÞ C 7.74 0.93 10.94 59.6 66.1 several efficiency criteria in Table 2. The results show that the
V 10.57 0.88 −20.40 83.6 67.6 ML solution performed very well but not as well as the DDS sol-
Note: C= calibration; V= verification. ution. This is not surprising because the objective function used by
Table 3. Percentage Adjustments Obtained by DDS, the ML, and the 95% This shows that the insensitive parameters have minimal impact on
Confidence Interval Obtained by GL-DREAMðZSÞ and Optimal Parameter runoff from the watershed and thus on the likelihood function. The
Values Determined for One of the Subcatchments in the Ballona Creek ML adjustments recommended for these insensitive parameters are
Watershed meaningless as any adjustment within the wide range of the pos
Percentage adjustments Parameter valuesa terior PDFs will produce almost identical runoff and likelihood
GL-DREAMðZSÞ function value.
One interesting observation is the recommendation by both ML
Lower Upper Initial Optimal
and DDS solutions to decrease percent imperviousness of the
Parameter DDS ML bound bound value value
watershed by approximately 15% so that model simulations closely
Width 137.1 117.0 −12.6 126.3 460 998 match observed runoff. Given that only the largest 72 storm drains
Slope 2.7 −7.4 −40.1 99.5 1.79 1.66 were considered in the study by ignoring hundreds of smaller storm
% Imperv −14.9 −15.0 −15.0 −12.7 70 59 drains and channels in the watershed, one would intuitively expect
N-Imperv 32.1 79.3 5.8 99.5 0.01 0.02
a solution that speeds up travel time (e.g., increase in percent
N-Perv −42.8 38.6 −46.5 99.0 0.02 0.03
Dstore-Imperv 99.8 98.2 84.9 99.9 0.25 0.50 imperviousness and steeper slope) to compensate for impact of
Dstore-Perv 234.7 1.8 72.2 293.6 1.27 1.29 the ignored storm drains on travel time. The suggested decrease
% Zero-Imperv −99.9 −100.0 −100.0 −94.5 100 0 in percentage imperviousness is believed to be related to the
Curve Number −3.2 8.4 −17.5 0.7 69 75 assumption made regarding connectivity of the impervious and
Drying Time −6.9 −93.1 −68.0 89.7 4 0.28 pervious subareas in each subwatershed. Both subarea types were
Conduit n 31.0 28.0 34.0 96.0 0.013 0.016 assumed to directly flow to the outlet of each subwatershed,
a
Initial values refer to the actual initial parameter values assigned to one of whereas in reality only a fraction of the impervious subareas
the subcatchments, and the optimal values are determined from the ML may be directly connected to the engineered drainage system. This
percent adjustments and the initial parameter values. See Table 1 for units. assumption might have led to overestimation of the effective per
centage imperviousness of the watershed in the model.
2
Predictive Uncertainty
1.8
Understanding the total predictive uncertainty associated with
model simulations is very essential for decision makers. Modeling
1.6 uncertainties are believed to arise from imprecise knowledge of the
R Statistic
0 0 0 0
-20 Width 150 -100 Slope 125 -15 Imperv -12 -20 120
N-Imperv
0 0 0 0
-75 100 80 100 0 320 -100 % Zero Imperv -94
Dstore-Imperv Dstore-Perv
N-Perv
0 0 0
-20 CN 10 -100 Drying Time 100 30 Manning's n 100
(a)
0 0 0
0 1 0.117 0.12 0.96 1
s0 x 10-6 s1 β
2000 1500
750
1000 750
0 0 0
0.845 0.85 -0.105 φ2 -0.1 0.185
φ3
0.205
(b) φ1
Fig. 4. Posterior PDFs determined using GL-DREAMðZSÞ for (a) SWMM5 runoff parameters; (b) GL error model parameters; ordinate of the plots
show frequency of occurrence
Rainfall [mm]
The runoff simulated using the optimal solution identified by the
single-objection calibration attempt was in good agreement with
the observed counterparts when evaluated graphically and by using
several goodness-of-fit measures. However, diagnostic analysis of
the residuals indicates that the assumptions of homoscedasticity,
temporal independence, and Gaussian distribution commonly made
in such traditional calibration models are unjustified. On the other
hand, the maximum likelihood solutions determined using the UA
model produced runoff simulations that are of comparable accuracy
Time Count [15 Minute Interval] with that of the single-objection calibration solutions while accu
rately characterizing the structure of the model residuals. The
Fig. 5. Comparison of observed runoff (dotted) to the 95% confidence
assumptions made by the error model used in the UA methodology
interval bounds (lines) determined using GL-DREAMðZSÞ , considering
were found consistent with the characteristics of the residuals
only parameter uncertainty for calibration and verification periods
generated from the ML solution.
600
(a)
Streamflow [m /s]
400
3
200
-200
0 100 200 300 400 500 600
Time Count [15 minute interval]
40 1 1
Partial autocorrelation
(b) (c) (d)
20 0.8
Residuals (a)
0.5
Density
0.6
0
0.4
0
-20 0.2
-40 0 -0.5
0 100 200 300 400 -2 0 2 0 10 20 30 40 50
Residuals (a) Lag [15 minute interval]
Simulated flow [m3/s]
600
(e)
Streamflow [m /s]
400
3
200
-200
0 50 100 150 200 250 300 350 400 450 500 550
Time Count [15 minute interval]
20 1 1
Partial autocorrelation
0.5
Density
0 0.5
0
-10
-20 0 -0.5
0 100 200 300 -2 0 2 0 10 20 30 40 50
Simulated flow [m3/s] Residuals (a) Lag [15 minute interval]
Fig. 6. Total predictive uncertainty (95% confidence interval) and diagnostic plots of the residuals obtained using GL-DREAMðZSÞ for calibration
(a–d) and verification (e–h) periods: (a and e) comparison between observed streamflow (dots) to the 95% confidence bounds (solid line); (b and f)
show that the residuals are homoscedastic; (c and g) comparison of observed PDF (dotted) of the residuals to the assumed distribution (solid line) in
GL; (d and h) show autocorrelation of the residuals; dashed lines in (d) and (h) show lower and upper bounds of the 95% confidence interval
In addition to accurately simulating runoff and properly charac Laboratory, Presented at 1st Ballona Wetlands Symp., Los Angeles,
terizing the residuals, the UA model has successfully determined CA.
parameter uncertainty and total predictive uncertainty. The param Deletic, A., et al. (2012). “Assessing uncertainties in urban drainage
eter posteriors showed that eight of the 11 SWMM5 parameters models.” Phys. Chem. Earth, Parts A/B/C, 42–44(2012), 3–10.
considered for the analysis exhibited wide uncertainty bound, Dotto, C. B. S., et al. (2012). “Comparison of different uncertainty
techniques in urban stormwater quantity and quality modelling.” Water
whereas the runoff simulated for the watershed considering param
Res., 46(8), 2545–2558.
eter uncertainty alone showed no appreciable variability. This sug
Duan, Q., Sorooshian, S., and Gupta, V. K. (1992). “Effective and efficient
gests that runoff is sensitive only to three (i.e., Imperv, Dstore global optimization for conceptual rainfall-runoff models.” Water
Imperv, % Zero Imperv) of the 11 parameters. Additionally, the ML Resour. Res., 28(4), 1015–1031.
solution identified by the UA model and the optimal solution EPA. (2007). “National water quality inventory: Report to congress; 2002
determined by DDS showed good agreement only for four (Imperv, reporting cycle.” Document No. EPA-841-R-07-001, Washington, DC.
Dstore Imperv, % Zero Imperv, Conduit n) of the 11 SWMM5 Fang, T., and Ball, J. E. (2007). “Evaluation of spatially variable parameters
parameters confirming nonidentifiablity of the insensitive parame in a complex system: An application of a genetic algorithm.” J. Hydro-
ters. Results also suggest that contribution of parameter uncertainty inform., 9(3), 163–173.
to total predictive uncertainty is insignificant for the study water Feyen, L., Khalas, M, and Vrugt, J. A. (2008). “Semi-distributed parameter
shed, underlying the importance of the other sources of predictive optimization and uncertainty assessment for large-scale streamflow sim
uncertainty for Ballona Creek watershed. ulation using global optimization.” Hydrol. Sci. J., 53(2), 293–308.
The 95% confidence interval determined for total predictive Freni, G., Mannina, G., and Viviani, G. (2008). “Uncertainty in urban
uncertainty using the UA model bracketed the majority of the stormwater quality modelling: The effect of acceptability threshold
in the GLUE methodology.” Water Res., 42(8–9), 2061–2072.
observed data, demonstrating reasonable accuracy of the UA
Freni, G., Mannina, G., and Viviani, G. (2009a). “Uncertainty in urban
result. Satisfactory total predictive uncertainty bounds were gen
stormwater quality modelling: The influence of likelihood measure
erated for both calibration and verification periods although the formulation in the GLUE methodology.” Sci. Total Environ., 408(1),
verification period results seem less adequate. Overall, the UA 138–145.
methodology proved promising for sensitivity analysis, calibra Freni, G., Mannina, G., and Viviani, G. (2009b). “Urban runoff modelling
tion, parameter uncertainty, and total predictive uncertainty uncertainty: Comparison among Bayesian and pseudo-Bayesian meth
analysis of urban storm water drainage models at least for the ods.” Environ. Modell. Software, 24(9), 1100–1111.
short-simulation durations considered in this study. Applications Gelman, A., and Rubin, D. B. (1992). “Inference from iterative simulation
to additional watersheds in other hydroclimatic regions can help using multiple sequences.” Stat. Sci., 7(4), 457–472.
further examine this potential. The subsequent study will investi Gironás, J., Niemann, J. D., Roesner, L. A., Rodriguez, F., and Andrieu, H.
gate application to continuous simulations and on ways to use the (2010). “Evaluation of methods for representing urban terrain in storm-
predictive uncertainty in decision making such as for optimal water modeling.” J. Hydrol. Eng., 15(1), 1–14.
SCM design applications. Gotzinger, J., and Bardossy, A. (2008). “Generic error model for calibration
and uncertainty estimation of hydrological models.” Water Resour.
Res., 44(12), W00B07.
Green, W. H., and Ampt, G. (1911). “Studies of soil physics. Part 1. The
References flow of air and water through soils.” J. Agr. Soc., 4, 1–24.
Holland, J. H. (1975). Adaptation in natural and artificial systems, Uni
Amenu, G. G. (2011). “A comparative study of water quality conditions versity of Michigan Press, Ann Arbor, MI.
between heavily urbanized and less urbanized watersheds of Los Horton, R. E. (1937). “Determination of infiltration-capacity for large
Angeles Basin.” Proc., World Environmental and Water Resources drainage basins.” Trans. Am. Geophys. Union, 18, 371–385.
Congress 2011: Bearing Knowledge for Sustainability, ASCE, Reston, Kavetski, D., Kuczera, G., and Franks, S. W. (2006). “Bayesian analysis of
VA, 680–690. input uncertainty in hydrological modeling: 1. Theory.” Water Resour.
Ball, J. E. (2009). “Discussion of ‘Automatic calibration of the U.S. EPA Res., 42(3), W03407.
SWMM model for a large urban catchment’ by J. Barco, K. M. Wong, Kuczera, G. (1983). “Improved parameter inference in catchment models,
and M. K. Stenstrom.” J. Hydraul. Eng., 135(12), 1108–1110. 1. Evaluating parameter uncertainty.” Water Resour. Res., 19(5),
Barco, J., Wong, K. M., and Stenstrom, M. K. (2008). “Automatic calibra 1151–1162.
tion of the U.S. EPA SWMM model for a large urban catchment.” Kuczera, G., Kavetski, D., Franks, S., and Thyer, M. (2006). “Towards a
J. Hydraul. Eng., 134(4), 466–474. Bayesian total error analysis of conceptual rainfall runoff models: Char
Bay, S., Jones, B. H., Schiff, K., and Washburn, L. (2003). “Water quality acterizing model error using storm-dependent parameters.” J. Hydrol.,
impacts of stormwater discharges to Santa Monica Bay.” Marine Envi 331(1–2), 161–177.
ron. Res., 56(1–2), 205–223. Kuczera, G., and Parent, E. (1998). “Monte Carlo assessment of parameter
Beven, K. (2006). “A manifesto for the equifinality thesis.” J. Hydrol., uncertainty in conceptual catchment models: The Metropolis algo
320(1–2), 28–36. rithm.” J. Hydrol., 211(1–4), 69–85.
Beven, K., and Freer, J. (2001). “Equifinality, data assimilation, and uncer Laloy, E., and Vrugt, J. A. (2012). “High-dimensional posterior exploration
tainty estimation in mechanistic modeling of complex environmental of hydrologic models using multiple-try DREAM(ZS) and high-perfor
systems using the GLUE methodology.” J. Hydrol., 249(1–4), 11–29. mance computing.” Water Resour. Res., 48, W01526.
Beven, K. J., and Binley, A. M. (1992). “The future of distributed models: Los Angeles County Department of Public Works (LACDPW). (2011).
Model calibration and uncertainty prediction.” Hydrol. Processes, 6(3), “Ballona Creek Watershed.” 〈http://ladpw.org/wmd/watershed/bc/〉
279–298. (Jul. 26, 2011).
Beven, K. J., and Freer, J. (2001). “Equifinality, data assimilation, and un Mannina, G., and Viviani, G. (2010). “An urban drainage stormwater
certainty estimation in mechanistic modeling of complex environmental quality model: Model development and uncertainty quantification.”
systems using the GLUE methodology.” J. Hydrol., 249(1–4), 11–29. J. Hydrol., 381(3–4), 248–265.
Box, G. E. P., and Tiao, G. C. (1992). Bayesian inference in statistical Mantovan, P., and Todini, E. (2006). “Hydrological forecasting uncertainty
analysis, Wiley, New York. assessment: Incoherence of the GLUE methodology.” J. Hydrol.,
Burian, S. J., McPherson, T. N., Brown, M. J., and Turin, H. J. (2000). 330(1–2), 368–381.
“Development of a stormwater model for the Ballona Creek watershed.” Montanari, A., Shoemaker, C. A., and van de Giesen, N. (2009).
Los Alamos Unclassified Rep. (LA-UR-00-1849), Los Alamos National “Introduction to special section on uncertainty assessment in surface
and subsurface hydrology: An overview of issues and challenges.” Stein, E. D., and Tiefenthaler, L. L. (2005). “Dry-weather metals and bac
Water Resour. Res., 45(12), W00B00. teria loading in an arid, urban watershed: Ballona Creek, California.”
Moradkhani, H., and Sorooshian, S. (2008). “General review of rainfall- Water Air Soil Pollut., 164(1–4), 367–382.
runoff modeling: Model calibration, data assimilation, and uncertainty Stenstrom, M. K., and Strecker, E. (1993). “Assessment of storm drain
analysis.” Hydrological modeling and water cycle: Coupling of the sources of contaminants to Santa Monica Bay: Annual pollutants load
atmospheric and hydrological models, Vol. 63, Springer, Water Science ings to Santa Monica Bay from stormwater runoff.” Rep. No. UCLA
and Technology Library, Berlin, Heidelberg, 1–23. ENG-93-62, Vol. 1, Univ. of California, Los Angeles, 1–248.
Muleta, M. K. (2010). “Comparison of model evaluation methods to Thiemann, M., Trosset, M., Gupta, H. V., and Sorooshian, S. (2001).
develop a comprehensive watershed simulation model.” World Environ “Bayesian recursive parameter estimation for hydrologic models.”
mental and Water Resources Congress 2010, ASCE, Reston, VA,
Water Resour. Res., 37(10), 2521–2535.
2492–2501.
Thyer, M., Renard, B., Kavetski, D., Kuczera, G., Franks, S. W., and
Muleta, M. K. (2012). “Sensitivity of model performance to the efficiency
Srikanthan, S. (2009). “Critical evaluation of parameter consistency
criteria used as objective function during automatic calibrations.”
J. Hydrol. Eng., 17(6), 756–767. and predictive uncertainty in hydrological modeling: A case study
Muleta, M. K., and Nicklow, J. W. (2005). “Sensitivity and uncertainty using Bayesian total error analysis.” Water Resour. Res., 45(12),
analysis coupled with automatic calibration for a distributed watershed W00B14.
model.” J. Hydrol., 306(1–4), 127–145. Tolson, B. A., and Shoemaker, C. A. (2007). “Dynamically dimensioned
Nash, J. E., and Sutcliffe, J. V. (1970). “River flow forecasting through search algorithm for computationally efficient watershed model calibra
conceptual models. Part I- A discussion of principles.” J. Hydrol., tion.” Water Resour. Res., 43(1), W01413.
125, 277–291. USGS. (2013). “The national map.” USGS, Washington, DC. 〈http://
National Oceanic, and Atmospheric Administration (NOAA). (2013). seamless.usgs.gov/〉 (Jun. 25, 2013).
“Atlas 14 point precipitation frequency estimates.” Silver Spring, Vrugt, J. A., Diks, C. G. H., Gupta, H. V., Bouten, W., and Verstraten, J. M.
MD. 〈http://hdsc.nws.noaa.gov/hdsc/pfds/pfds_map_pr.html〉 (Jun. 25, (2005). “Improved treatment of uncertainty in hydrologic modeling:
2013). Combining the strengths of global optimization and data assimilation.”
Natural Resources Conservation Service (NRCS). (2013). “Soil data mart.” Water Resour. Res., 41(1), W01017.
U.S. Dept. of Agriculture, Washington, DC. 〈http://soildatamart.nrcs Vrugt, J. A., Gupta, H. V., Bouten, W., and Sorooshian, S. (2003).
.usda.gov/〉 (Jun. 25, 2013). “A Shuffled Complex Evolution Metropolis algorithm for optimization
National Research Council (NRC). (2008). Urban stormwater management and uncertainty assessment of hydrologic model parameters.” Water
in the United States, The National Academies Press, Washington, DC. Resour. Res., 39(8), SWC 1-1–SWC 1-14.
Rossman, L. A. (2010). “Stormwater management model user’s manual, Vrugt, J. A., Ter Braak, C. J. F., Clark, M. P., Hyman, J. M., and Robinson,
version 5.” EPA/600/R-05/040, U.S. EPA, Cincinnati.
B. A. (2008). “Treatment of input uncertainty in hydrologic modeling:
Schoups, G., and Vrugt, J. A. (2010). “A formal likelihood function
Doing hydrology backward with Markov chain Monte Carlo simula
for parameter and predictive inference of hydrologic models with cor
tion.” Water Resour. Res., 44(12), W00B09.
related, heteroscedastic, and non-Gaussian errors.” Water Resour. Res.,
46(10), W10531. Vrugt, J. A., Ter Braak, C. J. F., Diks, C. G. H., Higdon, D., Robinson,
Soil Conservation Service. (1964). “National engineering handbook.” Sec B. A., and Hyman, J. M. (2009a). “Accelerating Markov chain Monte
tion 4, Hydrology, Department of Agriculture, Washington, DC, 450. Carlo simulation by differential evolution with self-adaptive random
Sorooshian, S., and Dracup, J. A. (1980). “Stochastic parameter estimation ized subspace sampling.” Int. J. Nonlinear Sci. Numer. Simul., 10(3),
procedures for hydrologic rainfall-runoff models: Correlated and heter 271–288.
oscedastic error cases.” Water Resour. Res., 16(2), 430–442. Vrugt, J. A., Ter Braak, C. J. F., Gupta, H. V., and Robinson, B. A. (2009b).
Stedinger, J. R., Vogel, R. M., Lee, S. U., and Batchelder, R. (2008). “Equifinality of formal (DREAM) and informal (GLUE) Bayesian ap
“Appraisal of the generalized likelihood uncertainty estimation (GLUE) proaches to hydrologic modeling?” Stoch. Environ. Res. Risk Assess.,
method.” Water Resour. Res., 44(12), W00B06. 23(7), 1011–1026.