The main objective of this paper is to present a model for generating synthetic rainfall totals on various timescales to be applicable for a variety of uses. Many large-scale ecological and water resource models require daily, monthly and yearly rainfall data as input to the model. As historical data provides only one realisation, synthetic generated rainfall totals are needed to assess the impact of rainfall variability on water resources systems (Srikanthan, In: MODSIM2005, Melbourne, Dec. 2005, pp. 1915–1921, 2005). Thus, our preferred model should simulate rainfall for yearly, monthly and daily periods. We believe that, for water supply issues, no higher resolution is needed, although higher resolution would be useful in models designed to measure the risk of local flooding. The critical factors are daily, monthly and yearly totals and daily, monthly and yearly variation. A model for generating yearly totals will be described using traditional time series methods. This model, along with a similarly constructed daily generation model by Piantadosi et al. (A New Model for Correlated Daily Rainfall, 2008), will be cascaded to start with a synthetic yearly total, then generate a synthetic sequence of monthly totals (through selection from a large number of realisations) that match the yearly total, and subsequently perform a similar operation for sequences of daily totals to match the required monthly totals. We present a new model for the generation of synthetic monthly rainfall data, which we demonstrate for Parafield in Adelaide, South Australia. The rainfall for each month of the year is modelled as a non-negative random variable from a mixed distribution with either a zero outcome or a strictly positive outcome. We use maximum likelihood to find parameters for both the probability of a zero outcome and the gamma distribution that best matches the observed probability density for the strictly positive outcomes. We describe a new model that generates correlated monthly rainfall totals using a diagonal band copula with a single parameter to generate lag-1 correlated random numbers. Our model preserves the marginal monthly distributions and, hence, also preserves the monthly and yearly means. We show that, for Parafield, the correlation between rainfall totals for successive months is not significant, and so, it is reasonable to assume independence. This is, however, not true for daily rainfall. We describe a new model that generates correlated daily rainfall totals using a diagonal band copula with a single parameter to generate lag-1 correlated random numbers.
[1]
D. Wilks,et al.
The weather generation game: a review of stochastic weather models
,
1999
.
[2]
T. McMahon,et al.
Stochastic generation of annual, monthly and daily climate data: A review
,
2001
.
[3]
R M Cooke,et al.
Monte Carlo sampling for generalized knowledge dependence with application to human reliability.
,
1986,
Risk analysis : an official publication of the Society for Risk Analysis.
[4]
M. Parlange,et al.
Overdispersion phenomenon in stochastic modeling of precipitation
,
1998
.
[5]
William H. Press,et al.
Book-Review - Numerical Recipes in Pascal - the Art of Scientific Computing
,
1989
.
[6]
Richard Coe,et al.
A Model Fitting Analysis of Daily Rainfall Data
,
1984
.
[7]
F. A. Seiler,et al.
Numerical Recipes in C: The Art of Scientific Computing
,
1989
.
[8]
Roger David Braddock,et al.
A model for seasonal variation of rainfall at Adelaide and Turen
,
1996
.
[9]
R. Nelsen.
An Introduction to Copulas
,
1998
.
[10]
William H. Press,et al.
The Art of Scientific Computing Second Edition
,
1998
.
[11]
Phil Howlett,et al.
Matching the grade correlation coefficient using a copula with maximum disorder
,
2007
.
[12]
Phil Howlett,et al.
Analysis of a Practical Control Policy for Water Storage in Two Connected Dams
,
2005
.
[13]
Daniel Lewandowski,et al.
Generalized diagonal band copulas
,
2005
.
[14]
G. E. Berrios,et al.
. New York: Cambridge
,
2000
.
[15]
Phil Howlett,et al.
Simulation of monthly rainfall totals
,
2008
.
[16]
R Srikanthan.
Stochastic Generation of Daily Rainfall Data
,
2005
.
[17]
D. C. Bowden,et al.
MAXIMUM LIKELIHOOD ESTIMATION FOR MIXTURES OF TWO NORMAL DISTRIBUTIONS
,
1973
.