With the first-order Markov-chain model, the probability of rain on a given day is conditioned on the wet or dry status of the previous day. A wet day is defined as a day with 0.1 mm of rain or more.
The user is required to input the probability of a wet day on day given a wet day on day , and the probability of a wet day on day given a dry day on day , for each month of the year. From these inputs the remaining transition probabilities can be derived:
1:3.1.1
1:3.1.2
where is the probability of a dry day on day given a wet day on day and is the probability of a dry day on day given a dry day on day .
To define a day as wet or dry, SWAT+ generates a random number between 0.0 and 1.0. This random number is compared to the appropriate wet-dry probability, or . If the random number is equal to or less than the wet-dry probability, the day is defined as wet. If the random number is greater than the wet-dry probability, the day is defined as dry.
The daily precipitation generator is a Markov chain-skewed (Nicks, 1974) or Markov chain-exponential model (Williams, 1995). A first-order Markov chain is used to define the day as wet or dry. When a wet day is generated, a skewed distribution or exponential distribution is used to generate the precipitation amount. Table 1:3-1 lists SWAT+ input variables that are used in the precipitation generator.
Numerous probability distribution functions have been used to describe the distribution of rainfall amounts. SWAT+ provides the user with two options: a skewed distribution and an exponential distribution.
The skewed distribution was proposed by Nicks (1974) and is based on a skewed distribution used by Fiering (1967) to generate representative streamflow. The equation used to calculate the amount of precipitation on a wet day is:
1:3.1.3
where is the amount of rainfall on a given day (mm HO), is the mean daily rainfall (mm HO) for the month, is the standard deviation of daily rainfall (mm HO) for the month, is the standard normal deviate calculated for the day, and is the skew coefficient for daily precipitation in the month.
The standard normal deviate for the day is calculated:
1:3.1.4
where and are random numbers between 0.0 and 1.0.
The exponential distribution is provided as an alternative to the skewed distribution. This distribution requires fewer inputs and is most commonly used in areas where limited data on precipitation events is available. Daily precipitation is calculated with the exponential distribution using the equation:
1:3.1.5
where is the amount of rainfall on a given day (mm HO), is the mean daily rainfall (mm HO) for the month, is a random number between 0.0 and 1.0, and is an exponent that should be set between 1.0 and 2.0. As the value of is increased, the number of extreme rainfall events during the year will increase. Testing of this equation at locations across the U.S. have shown that a value of 1.3 gives satisfactory results.
Table 1:3-1: SWAT+ input variables that pertain to generation of precipitation.
Definition | Source Name | Input Name | Input File |
---|---|---|---|
Precipitation input: 'sim' for simulated or gage name
pgage
pcp
: probability of a wet day following a dry day in month
pr_wd
wet_dry
: probability of a wet day following a wet day in month
pr_ww
wet_wet
Rainfall distribution code: 0-skewed, 1-exponential
IDIST
file.cio
: value of exponent (required if IDIST = 1)
REXP
file.cio
average amount of precipitation falling in month (mm H2O)
pcpmm
pcp_ave
average number of days of precipitation in month(= PCPMM / PCPD)
pcpd
pcp_days
: standard deviation for daily precipitation in month (mm HO)
pcpstd
pcp_sd
: skew coefficient for daily precipitation in month
pcpskw
pcp_skew