Flood prediction in NSW (a comment on Micevski et al., 2014)

Micevski, T., A. Hackelbusch, K. Haddad, G. Kuczera and A. Rahman (2014). “Regionalisation of the parameters of the log-Pearson 3 distribution: a case study for New South Wales, Australia.” Hydrological Processes. DOI: 10.1002/hyp.10147 DOI: 10.1002/hyp.10147

Micevski et al develop a regional flood frequency analysis approach for the eastern part of NSW. In this post I present an example based on the simplest outputs and interpretation of their paper. Please refer to their paper for all the details and if you are interested in the sophisticated method they propose.

Micevski et al suggest that, when considering this area of NSW as a whole, it is possible to estimate the three parameters of the Log Pearson 3 (LPIII) distribution of flood magnitudes i.e. the mean, standard deviation and skewness.  In fact, the standard deviation and skewness are constant. Standard deviation is 1.2318 and skewness is -0.3567 (see table 3 of the paper).

From table 3 in the paper, an equation can be developed for the mean of the LPIII distribution.

\text{mean} = 0.3326 + 1.322 \log_{10}A +4.7182 \log_{10}(I_{12,50})

Where, A = area in km2 and I12,50 is the rainfall intensity for a 12 hour, 50 year ARI rain storm.

We can use this relationship to estimate flood magnitudes anywhere in eastern NSW.


Estimating the 100-year flood for Gillamatong Creek near Braidwood.

A flood study for creeks near Braidwood was undertaken by Cardo Willing in 2005 (link to report). They estimated the 100-year flood for Gillamatong Creek as 518 m3/s.

Lets make a similar estimate using the approach suggested in Micevski et al. (2014).

area = 77 km2

Lat = -35.45o

Lon = 149.8o

Looking at the 1987 IFD values from the BoM website I12,50= 12.4 mm/h

We will undertake the analysis for a 100-year flood. Parameters can be defined and the mean flood estimated as follows.

area <- 77
I12.50 <- 12.4
flood.skew <- -0.3567
flood.sd <- 1.2318
flood.aep <- 1/100

flood.mean <- 0.3326 + 1.322 * log10(area) + 4.7182 * log10(I12.50)

## [1] 7.986

To calculate the 100 year flood, we use the standard frequency factor approach for the LPIII distribution – however this is the tricky part of the Micevski et al.’s paper; they use natural logs for flood magnitudes and base 10 logs for catchment attributes. Also the flow units are ML/d not cumec. Therefore, the equation we need is:

\log_e (Q) = \mu + \sigma K(\gamma)

Where K is the frequency factor.

The frequency factor for the 100 year flood can be estimated in R using using the quape3 function from the lmom package. Frequency factors can also be looked up in Australian Rainfall and Runoff (Book IV, Table 2.2). Calculations for the 100 year flood are:

flood.K <- quape3(1 - flood.aep, para = c(0, 1, flood.skew))

## [1] 2.062

So our loge(Q100) will be:

flood.mean + flood.sd * flood.K

## [1] 10.53

Therefore Q100 = 3.72 × 104 ML/d or 431 cumec.

So the estimate is smaller than that of Cardo and Willing but in the same ballpark and reasonable for a quick approximate method.

A single snippet of code to calculate a flood estimate is as follows:

area <- 77
I12.50 <- 12.4
flood.skew <- -0.3567
flood.sd <-  1.2318
flood.aep <- 0.01
B0 <- 0.3326
B1 <- 1.3222
B2 <- 4.7182 
Q100 <- exp(B0 + B1*log10(area) + B2*log10(I12.50) + flood.sd * quape3(1-flood.aep, para = c(0, 1, flood.skew)))

It’s generally good practice, when providing the estimate of some statistic, to also estimate the uncertainty.  Micevski et al., list means and standard deviations of all their coefficients although not the covariances. Figure 5 in Micevski et al’s paper show uncertainties (posterior distributions) for a few stream gauges.

Leave a Reply

Fill in your details below or click an icon to log in:

WordPress.com Logo

You are commenting using your WordPress.com account. Log Out / Change )

Twitter picture

You are commenting using your Twitter account. Log Out / Change )

Facebook photo

You are commenting using your Facebook account. Log Out / Change )

Google+ photo

You are commenting using your Google+ account. Log Out / Change )

Connecting to %s