Skip to main content

Uncertainty analysis for flood inundation modelling with a random floodplain roughness field

Abstract

Background

This study tested a first-order perturbation method based on Karhunen-Loevè expansion (FP-KLE), to analyze flood inundation modeling under uncertainty. The floodplain roughness over a 2-dimensional domain was assumed to be a statistically heterogeneous field with log-normal distributions. Firstly, we attempted to use KLE to decompose the random field of log-transferred floodplain roughness N(x), which was based on the eigenvalues and eigenfunctions of the covariance function of N(x), and a set of orthogonal normal random variables. Secondly, the maximum flow depths were expanded by the first-order perturbation method by using the same set of random variables as used in the KLE decomposition. Then, a flood inundation model, named FLO-2D, was adopted to numerically solve the corresponding perturbation expansions.

Results

To illustrate the methodology, a one-in-five-years flood event was chosen as the study case. The results indicated that the mean of the maximum flow-depth field obtained from the proposed method was fairly close to that from Monte Carlo Simulation (MCS), but the standard deviation was somewhat higher. However, the FP-KLE method was computationally more efficient than MCS.

Conclusions

The study verified the applicability of FP-KLE in handling uncertainties of flood modeling in a more efficient manner; further test with multiple inputs of random fields is desired.

Background

Flood inundation mapping is important for helping identify potential flood profiles, and assist decision makers in planning and mitigation operations (Jung and Merwade, 2012). Due to inherent complexity and large numbers of parameters, flood inundation modelling processes often involve large uncertainties, especially at topographical formation where no sufficient data are available for model parameterization, calibration and validation. As a key global parameter in flood inundation modelling, it has been recognized that the Manning’s roughness coefficient (Manning’s n) could notably affect flood inundation predictions (Beven and Binley, 1992; Aronica et al., 1998; Hall et al., 2005). Furthermore, due to different types of vegetation, surface irregularity and non-uniform/unsteady flows over the floodplain, the roughness coefficient has demonstrated a high heterogeneity (Pappenberger et al., 2005; Liu, 2009). Previously, many studies were devoted to analyse uncertainty propagation from manning’s roughness coefficients during flood modelling, such as Generalized Likelihood Uncertainty Estimation (GLUE) and traditional Monte Carlo simulations (MCS) (Aronica et al., 1998; Aronica et al., 2002; Van Vuren et al., 2005; Reza Ghanbarpour et al., 2011; Jung and Merwade, 2012). However, most of these studies tend to adopt homogenous values of roughness coefficient over the study domain, and this could lead to significant discrepancies between the observed data and simulation results.

It has been realized that the geological-related process involving spatial distribution of uncertainty, such as roughness coefficients, can be more practically considered as a random field (Ghanem and Spanos, 1991; Lu and Zhang, 2007). Subsequently, the governing equations describing one- or multi-dimensional physical processes will be transformed into stochastic. The simulation results (such as the flow depth and flow velocity) are probabilistic distributions of the outputs, instead of deterministic quantifications. In recent decades, new stochastic methodologies based on Karhunen-Loevè expansion (KLE) decomposition have been proposed, which could consider the spatial distribution of uncertain parameters in generating realizations (Huang et al., 2001; Zhang and Lu, 2004; Liu et al., 2006). KLE decomposition was firstly developed to deal with uncertain random field with spatial distribution and the corresponding uncertainties were treated with an orthogonal basis of polynomial chaos (Ghanem and Spanos, 1991). Roy and Grilli (1997) firstly applied KLE to decompose the random field of hydraulic conductivity in a 2-dimensional (2D) stationary groundwater flow modelling system. Later on, Zhang and Lu (2004) developed a new approach called KLME (which combined KLE with moment equation method) and applied it to analyze the first and second moments of the hydraulic head. Despite of quite a large number of analytical approaches been developed for assessing parameter uncertainty in flood modeling, studies on the uncertainty parameters with spatial distribution are relatively limited. Liu and Matthies (2010) attempted to combine KLE and Hermite polynomial chaos expansion and examine the uncertainty from inflow, topography and roughness coefficient over the entire flood modelling domain using stochastic 2D shallow water equations. In this study, KLE is to be tested in decomposing the random field of floodplain roughness coefficients (keeping the channel roughness coefficients fixed) within a coupled 1-dimensional (1D) (for channel flow) and 2D (for floodplain flow) physical flood inundation model (i.e. FLO-2D). The method is effective in saving computational efforts without compromising the accuracy of uncertainty assessment. Previously, no such attempt was made using FLO-2D.

Results and discussion

Study case

A flood inundation case modified from Aronica et al. (2002) is chosen to demonstrate the proposed FP-KLE method. The related settings are listed as follows (Aronica et al., 2002; Hall et al., 2005; Bates et al., 2008): (i) 50-meter resolution DEM is used as the topographical data with 0.25-meter accuracy vertically, changing from 67.73 to 83.79 m; this is relatively moderate compared to those of the steeply-changing mountain areas; and (ii) channel cross-section is defined as rectangular with the size of 25 m (width) by 1.5 m (depth); (iii) the upstream inflow hydrograph is suggested as steady at 73 m3/s; this corresponds to a flood event with a 5-years return period (i.e. December 1992 flood lasting about 27.8 hours). The channel roughness coefficient n is fixed at 0.03. More detailed description of the study case can be referred to Bates et al. (2008). For uncertainty assessment, the floodplain roughness is assumed as a random field with log-normal probability distribution function (PDF). Monte Carlo simulation (MCS), as a traditional stochastic method, could be used to deal with the stochastic partial differential equations, i.e., 2D shallow water equations. Although MCS is easy to use, it is restricted by extensive computational needs, particularly for large-scale flood inundation modelling problems. It is desired to apply KLE method for improving the computational efficiency.

Simplifications of KLE for the random floodplain roughness field

In order to make the uncertainty analysis with KLE more practical, some simplifications are made according to Zhang and Lu (2004) for 1D model: (i) x ^ = x / L x , where L x is the length of domain at the x direction, (ii) λ ^ i = λ i / L x , (iii) normalized correlation length η ^ = η / L x , and (iv) normalized eigenfuctions f ^ i x = f i x L x . For simplicity, the tilde is omitted for the above four expressions in the followed equations and expressions within the parts of results and discussion, and conclusion.

The log floodplain roughness coefficient is assumed an exponential spatial covariance function (Roy and Grilli, 1997):

C N x , y = C N x 1 , x 2 , y 1 , y 2 = σ N e - x 1 - y 1 η x - x 2 - y 2 η y
(1)

where η x and η y represent normalized correlation lengths in the x- and y-directions, respectively; (x 1 , y 1 ) and (x 1 , y 1 ) are the spatial Cartesian coordinates of two points located in a 2D physical domain.

For a 2D rectangular modelling domain, we have eigenvalues for a 2D field and can integrate them as (Zhang and Lu, 2004):

n = 1 λ n = i = 1 λ i j = 1 λ j / σ N 2 = D σ N 2
(2)

where λ n are the eigenvalues for 2D modelling domain; λ i (i = 1, 2,…, ∞) and λ j (j = 1, 2,…, ∞) are the eigenvalues for the x- and y-directions separately; σ N is the standard deviation of the log-transformed floodplain roughness coefficient; D is the size of the 2D modelling domain.

In this study, the random field of floodplain roughness coefficient n(x) is suggested as log-normal PDF with u n  = 0.06 (mean) and σ n  = 0.001 (standard deviation), spatially. The N(x) = ln[n(x)] is a normal distribution with mean u N  = -2.936 and standard deviation σ N  = 0.495. The range of N(x) is approximated as [u N -3σ N , u N +3σ N ] which is (0.012, 0.234). To achieve both efficiency in operationality and accuracy in computation, the number of KLE terms with different normalized correlation lengths may vary with different scenarios (i.e., various scales of the domain size) with specific model settings (i.e., boundary condition settings) and floodplain roughness coefficients (i.e., changing from rural to urban areas) under consideration. In this case, the numbers of terms retained in KLE expansion in the x-direction (m x ) and y-direction (m y ) are set as 20 and 10, respectively; hence, the total number of KLE terms is 20 × 10 = 200.

The eigenvalues would monotonically reduce as index n increases, as shown in Equation (2). Figure 1(a) shows that, for different exponential spatial covariance functions (with different normalized correlation length η), the declining rate λ ¯ n ( λ ¯ n = λ n / D σ N 2 ) is different. When η becomes larger, λ ¯ n would reduce more significantly. Figure 1(b) illustrates that the summation of λ ¯ n based on a finite number of terms rather than on an infinite number, can be considered as a function of the index n. The value of λ ¯ n would gradually approach to 1 when n is increasing.

Figure 1
figure 1

(a) Series of λ ¯ n ( λ ¯ n = λ n / D σ N 2 ) and (b) their finite cumulative sums for the 2D rectangular domain with different exponential spatial covariance functions (i.e. η=η x =η y = 0.15, 0.3, 1.0 and 4.0 respectively).

For this study case, the normalized correlation lengths are set as η x  = 0.15 and η y  = 0.3, and the total number of KLE terms is m = 200. Figure 2 shows the decreasing rate of eigenvalues and how much energy of KLE approximation is obtained. For example, if 200 KLE terms of N(x) expansion are used in KLE decomposition and the total energy of the approximation would save by 86.56%, as shown in Figure 2(b). Figures 3 shows two representations of the random fields of floodplain roughness coefficients over the 2D flood modelling domain with η x  = 0.15 and η y  = 0.3 and the number of KLE terms = 200. These figures show that the KLE decomposition of the uncertain random field is different from the Monte Carlo sampling, in which the heterogeneous profile of random field can be represented by smoother eigenpairs as expressed in Equation (2).

Figure 2
figure 2

(a) Series of λ ¯ n , and (b) its finite cumulative sum for 2D rectangular domain with defined exponential spatial covariance function.

Figure 3
figure 3

The random fields of floodplain roughness coefficients over the modelling domain for the (a) 5th and (b) 151th realizations. Domain size is divided into 76 (in x axis) × 48 (in y axis) grid elements. fn represent floodplain manning’s roughness coefficients.

Comparison with MCS

In order to verify the accuracy of the FP-KLE, the modelling results from 5,000 realizations of Monte Carlo simulations are also presented. Figure 4 shows the distribution statistics of the maximum flow depths h (x) using KLE and MCS, respectively. From Figures 4a) and 4(b), it can be seen that the mean of h (x) distribution (varying from 0 to 2.5 m) simulated by KLE reproduces well the result from MCS. From Figures 4(c) and 4(d), different distributions of the standard deviation of h (x) are found. The standard deviation of h(x) simulated by FP-KLE is somewhat higher than that calculated by MCS. This may because FP-KLE is in lower order (i.e. first-order) and less capable of achieving a high accuracy, comparing with MCS.

Figure 4
figure 4

Comparison of statistics of maximum flow depth field simulated by the first-order perturbation based on the KLE and MCS: (a) and (b) are the mean maximum depth distributions calculated by FP-KLE and MCS, respectively; (c) and (d) are the standard deviation of the maximum depth distributions calculated by FP-KLE and MCS, respectively. Domain size is divided into 76 (in x axis) × 48 (in y axis) grid elements.

Figure 5 shows a comparison of the statistics of the h (x) field along the cross-section x = 43 between FP-KLE and MCS. It seems that the mean of the h(x) along the concerned cross section simulated by FP-KLE fits very well with that simulated by MCS. However, the standard deviation from the perturbation method is higher than that from MCS. For example, at the location (x, y) = (43, 30), the standard deviation of the h (x) (i.e. 0.01262 m) is about 66.42% higher than that (i.e. 0.007583 m) calculated by MCS; the range of the standard deviation by FP-KLE is from 0 to 0.01262 m, and that by MCS is from 0 to 0.00883 m. It indicates that the FP-KLE with 200 terms may not sufficiently capture the simulated standard deviation results by MCS.

Figure 5
figure 5

Comparison of statistics of max flow depth field simulated by FP-KLE and MCS along the profile x = 43/76: (a) mean max flow depth, and (b) standard deviation of max flow depth.

Generally, the FP-KLE is capable of quantifying uncertainty propagation in a highly heterogeneous flood modelling system. By comparison, FP-KLE is proved to be more efficient than traditional MCS in terms of computational efforts. The presented approach can be used for large-scale flood domains with high spatial-variability of input parameters, and it could provide reliable predictions to the decision-makers in flood risk assessment with relatively a small number of model runs.

Conclusions

This study attempted to use a first-order perturbation called FP-KLE to investigate the impact of uncertainty associated with floodplain roughness coefficients on a 2D flooding modelling process. Firstly, the KLE decomposition for the log-transformed floodplain random field was made within a 2D rectangular flood domain represented by pairs of eigenvalue and eigenfunctions. Secondly, the first-order expansion of h (x) perturbation was applied to the maximum flow depth distribution. Thirdly, the flood inundation model, i.e. FLO-2D, was used to solve each term of the perturbation based on the FP-KLE approach. Finally, the results were compared with those obtained from traditional Monte Carlo simulation (MCS).

The following facts were found from this study: (i) for the 2D flood case, with parameter setting of η x  = 0.15, η y  = 0.3, and the number of KLE terms = 200, about 86.56% energy have been saved; this was considered sufficient for reproduction of statistical characteristics; (ii) the mean of h (x) field from FP-KLE reproduced well the results from MCS, but the standard deviation was somewhat higher; (iii) the first-order KLE-based perturbation method was computationally more efficient than MCS with comparable accuracy. Some limitations need further discussions in future studies: (i) compared with the first-order KLE-based perturbation approach, the second-order (or higher orders) perturbation may lead to more accurate result, but the required computational effort would increase dramatically; further test of the method on higher orders is desired; (ii) for this study, the simulation is in a steady-state condition; the KLE-based perturbation method for unsteady state could be further explored; (iii) the input random field in this study was assumed in normal distribution, non-normal distributions of the input random fields could be explored in the future.

Methods

Stochastic flood inundation model

FLO-2D is used in this study as the flood inundation model. It is based on the full 2D shallow water equations (also called dynamic wave momentum equation) (FLO-2D Software, 2012). In the FLO-2D modeling system, channel flow is 1D with the channel geometry represented by either rectangular or trapezoidal cross sections; meanwhile, the overland flow is modelled 2D as either sheet flow or flow in multiple channels (FLO-2D Software, 2012). Overbank flow in the channel is computed when the channel capacity is exceeded. Besides, an interface routine calculates the channel to floodplain flow exchange including return flow to the channel. More technical details of FLO-2D can be referred to Obrien et al. (1993), O'Brien et al. (1999), and D'Agostino and Tecca (2006).

To describe a 2D flood inundation process, shallow water equations can be used (Chow et al., 1988; Obrien et al., 1993; FLO-2D Software, 2012):

h x t + h x V = I
(3a)
S f = S o - h x - 1 g V V - 1 g V t
(3b)
S f = n x 2 R 4 3 V V
(3c)

where, h(x) is the flow depth; V represents the averaged-in-depth velocity in each direction x; x represents Cartesian coordinate spatially, such as x = (x, y) represents 2D Cartesian coordinate; t is the time; S o is the bed slope, and S f is the friction slope; and I is lateral flow into the channel from other sources. Equation (3a) is the continuity equation or mass conservation equation, and Equation (3b) is the momentum equation; both of them are the fundamental equations in the flood modelling. In Equation (3c), n is the Manning coefficient (Manning’s n), which is the most commonly applied friction parameter in flooding modelling. R is the hydraulic radius. Equations (3a)-(3c) are solved mathematically in eight directions by FLO-2D. In this study, n(x) is assumed as a random function spatially, and Equations (3a)-(3c) are transformed into stochastic partial differential equations with random floodplain roughness coefficients and other items within the model are considered to be deterministic. Our purpose is to solve the mean and standard deviation of the flow depth h(x), which are used to assess the uncertainty propagation during the flood inundation modelling.

Karhunen-Loevè expansion (KLE) of log floodplain manning’s roughness coefficients

Let the n(x) be a random field of floodplain manning’s roughness coefficient, and N(x) = ln[n(x, ω)], where x ϵ D and ω ϵ Ω (a probabilistic space). This means N( x ) changes spatially. Its covariance function C N (x , y) = <N’(x, ω) N’(y, ω) > is assumed as bounded, symmetric, and positive, which can be expressed as (Ghanem and Spanos, 1991):

C N x , y = n = 1 λ n f n x f n y
(4)

where λ n represent eigenvalues; f n (x) are eigenfunctions, which are orthogonal and determined by dealing with the Fredholm equation analytically or numerically as (Courant and Hilbert, 1989):

D C N x , y f n x d x = λ n f n y
(5)

where λ n and f n (x) for some specific covariance functions could be solved analytically (Zhang and Lu, 2004). The KLE representation of floodplain roughness coefficients can be expressed as:

N x , ω = μ x + n = 1 ξ n ω λ n f n x
(6)

where μ (x) represents the mean of N(x, ω); and ξ n (ω) are the orthogonal random variables with standard Gaussian distribution. For conciseness, ω is omitted in following related expressions and equations. The calculated λ n can be ranked as monotonically decreasing series with index n. In order to make the approximation effective in both computational effort and accuracy, the KLE representation of N(x) can be written in a finite form as follows (Liu et al., 2006; Li and Zhang, 2007):

N ˜ x = μ x + n = 1 m ξ n λ n f n x
(7)

Since λ n and f n (x) always appear in pairs, let g n x = λ n f n x . Then, Equation (7) can be simplified into:

N ˜ x = μ x + n = 1 m ξ n g n x
(8)

where N ˜ x is the approximation for the random floodplain roughness field. Generally, the more the terms used in g n (x), the more accurate the random floodplain roughness filed is represented; this implies that more random field energy is saved, but with an increasing computational efforts. For simplicity of expression, the tilde over the symbol N is omitted in the followed expressions and functions.

Perturbation method

In this study, the fluctuation of the max flow depths, as one of the important indicators of the flood inundation simulation, is affected by the spatial variability of the floodplain roughness values N(x). The maximum flow depths h(x) can be expressed with a perturbation expansion in an infinite series as follows (Phoon et al., 2002; Liu et al., 2006; Li and Zhang, 2007; Lu and Zhang, 2007):

h x = i = 0 h i x
(9)

where, h(i) (∙) is the ith order perturbation term based on the standard deviation of N(x) (denoted as σ N ).

Substituting Equations (8) and (9) into Equations (3a-3c), we can obtain the ith order term of the expansion h(i) (x) and each order of perturbation is calculated based on σ N . For example, the first-order perturbation expansion for h (x) can be expressed as h(x) = h(0) (x) + h(1) (x). It can be seen that the higher the order of the term h(i) (∙) kept in the expansion of h(x), the more energy could be retained in the expansion; hence more corrections are provided for the statistical moments (i.e. mean and variation) of the simulation results (Roy and Grilli, 1997). However, in this study, considering the computational requirements of the flood modelling, only the first-order perturbation expansion based on KLE is investigated.

References

  • Aronica G, Hankin B, Beven K: Uncertainty and equifinality in calibrating distributed roughness coefficients in a flood propagation model with limited data. Adv Water Resour 1998,22(4):349–365. 10.1016/S0309-1708(98)00017-7

    Article  Google Scholar 

  • Aronica G, Bates PD, Horritt MS: Assessing the uncertainty in distributed model predictions using observed binary pattern information within GLUE. Hydrol Process 2002,16(10):2001–2016. 10.1002/hyp.398

    Article  Google Scholar 

  • Bates P, Fewtrel T, Trigg M, Neal J: LISFLOOD-FP User manual and technical note, code release 4.3.6. University of Bristol; 2008.

    Google Scholar 

  • Beven K, Binley A: The future of distributed models: model calibration and uncertainty prediction. Hydrol Process 1992,6(3):279–298. 10.1002/hyp.3360060305

    Article  Google Scholar 

  • Chow VT, Maidment DR, Mays LW: Applied hydrology. McGraw-Hill International Editions, New York; 1988.

    Google Scholar 

  • Courant R, Hilbert D: Methods of Mathematical Physics. Interscience Publishers, New York; 1989.

    Book  Google Scholar 

  • D'Agostino V, Tecca PR: Some considerations on the application of the FLO-2D model for debris flow hazard assessment. In Monitoring, Simulation, Prevention and Remediation of Dense and Debris Flows. Edited by: Lorenzini G, Brebbia CA, Emmanouloudis DE. WIT Transactions on Ecology and the Envionment, 90; 2006:159–170.

    Chapter  Google Scholar 

  • FLO-2D Software, Inc: FLO-2D Reference Manual 2009. 2012.https://www.flo-2d.com/download

    Google Scholar 

  • Ghanem RG, Spanos PD: Stochastic Finite Elements: A Spectral Approach. Springer, New York; 1991.

    Book  Google Scholar 

  • Hall J, Tarantola S, Bates PD, Horritt MS: Distributed sensitivity analysis of flood inundation model calibration. Journal of Hydraulic Engineering-ASCE 2005,131(2):117–126. 10.1061/(ASCE)0733-9429(2005)131:2(117)

    Article  Google Scholar 

  • Huang SP, Quek ST, Phoon KK: Convergence study of the truncated Karhunen–Loeve expansion for simulation of stochastic processes. Int J Numer Methods Eng 2001,52(9):1029–1043. 10.1002/nme.255

    Article  Google Scholar 

  • Jung Y, Merwade V: Uncertainty Quantification in Flood Inundation Mapping Using Generalized Likelihood Uncertainty Estimate and Sensitivity Analysis. J Hydrol Eng 2012,17(4):507–520. 10.1061/(ASCE)HE.1943-5584.0000476

    Article  Google Scholar 

  • Li H, Zhang D: Probabilistic collocation method for flow in porous media: Comparisons with other stochastic methods. Water Resour Res 2007,43(9):W09409. doi:10.1029/2006WR005673 doi:10.1029/2006WR005673

    Google Scholar 

  • Liu D: Uncertainty Quantification with Shallow Water Equations. Thesis, Institute of Scientific Computing, Technical University of Braunschweig, Germany; 2009.

    Google Scholar 

  • Liu DS, Matthies HG: Uncertainty quantification with spectral approximations of a flood model. 9th World Congress on Computational Mechanics and 4th Asian Pacific Congress on Computational Mechanics, WCCM/APCOM 2010, July 19, 2010 - July 23, 2010. Institute of Physics Publishing, Sydney, NSW, Australia; 2010.

    Google Scholar 

  • Liu G, Zhang D, Lu Z: Stochastic uncertainty analysis for unconfined flow systems. Water Resour Res 2006,42(9):W09412.

    Google Scholar 

  • Lu Z, Zhang D: Stochastic Simulations for Flow in Nonstationary Randomly Heterogeneous Porous Media Using a KL‒Based Moment‒Equation Approach. Multiscale Modeling & Simulation 2007,6(1):228–245. 10.1137/060665282

    Article  Google Scholar 

  • O’Brien JS, Julien PY, Fullerton WT: Two-dimensional water flood and mudflow simulation. J of Hydraulic Engineering-ASCE 1993,119(2):244–261. 10.1061/(ASCE)0733-9429(1993)119:2(244)

    Article  Google Scholar 

  • O’Brien JS, Fullerton WT, Finch DM, Whitney JC, Kelly JF, Loftin SR: Simulation of Rio Grande floodplain inundation Using FLO-2D. In Rio Grande ecosystems: linking land, water, and people: toward a sustainable future for the Middle Rio Grande Basin. 1998 June 2–5. U.S. Department of Agriculture, Forest Service, Rocky Mountain Research Station, Albuquerque, NM. Proc. RMRS-P-7. Ogden, UT; 1999:52–60.

    Google Scholar 

  • Pappenberger F, Beven K, Horritt M, Blazkova S: Uncertainty in the calibration of effective roughness parameters in HEC-RAS using inundation and downstream level observations. J Hydrol 2005,302(1–4):46–69.

    Article  Google Scholar 

  • Phoon KK, Huang SP, Quek ST: Simulation of second-order processes using Karhunen-Loeve expansion. Comput Struct 2002,80(12):1049–1060. 10.1016/S0045-7949(02)00064-0

    Article  Google Scholar 

  • Reza Ghanbarpour M, Salimi S, Saravi MM, Zarei M: Calibration of river hydraulic model combined with GIS analysis using ground-based observation data. Res J Appl Sci Eng Technol 2011,3(5):456–463.

    Google Scholar 

  • Roy RV, Grilli ST: Probabilistic analysis of flow in random porous media by stochastic boundary elements. Engineering Analysis with Boundary Elements 1997,19(3):239–255. 10.1016/S0955-7997(97)00009-X

    Article  Google Scholar 

  • Van Vuren S, De Vriend H, Ouwerkerk S, Kok M: Stochastic modelling of the impact of flood protection measures along the river waal in the Netherlands. Nat Hazards 2005,36(1–2):81–102.

    Article  Google Scholar 

  • Zhang D, Lu Z: An efficient, high-order perturbation approach for flow in random porous media via Karhunen-Loève and polynomial expansions. J Comput Phys 2004,194(2):773–794. 10.1016/j.jcp.2003.09.015

    Article  Google Scholar 

Download references

Acknowledgement

This research was supported by Singapore’s Ministry of Education (MOM) AcRF Tier 1 Project (M4010973.030) and Tier 2 Project (M4020182.030).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Xiaosheng Qin.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

YH carried out the modelling studies and drafted the manuscript; XSQ helped out in the result analysis, and participated in the coordination and refinement of the manuscript. Both authors read and approved the final manuscript.

Authors’ original submitted files for images

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 2.0 International License (https://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Huang, Y., Qin, X. Uncertainty analysis for flood inundation modelling with a random floodplain roughness field. Environ Syst Res 3, 9 (2014). https://doi.org/10.1186/2193-2697-3-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/2193-2697-3-9

Keywords