arrow-left

All pages
gitbookPowered by GitBook
1 of 3

Loading...

Loading...

Loading...

1:3.1 Precipitation

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.

1:3.1.2 Amount of Precipitation

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:

Rday=ΞΌmon+2βˆ—Οƒmonβˆ—([(SNDdayβˆ’gmon6)βˆ—gmon6+1]3βˆ’1gmon)R_{day}=\mu_{mon}+2*\sigma_{mon}*(\frac{[(SND_{day}-\frac{g_{mon}}{{6}})*\frac{g_{mon}}{{6}}+1]^3-1}{g_{mon}})Rday​=ΞΌmon​+2βˆ—Οƒmonβ€‹βˆ—(gmon​[(SNDdayβ€‹βˆ’6gmon​​)βˆ—6gmon​​+1]3βˆ’1​) 1:3.1.3

where RdayR_{day}Rday​ is the amount of rainfall on a given day (mmΒ H2Omm\space H_2OmmΒ H2​O), ΞΌmon\mu_{mon}ΞΌmon​ is the mean daily rainfall (mmΒ H2Omm\space H_2OmmΒ H2​O) for the month, Οƒmon\sigma_{mon}Οƒmon​ is the standard deviation of daily rainfall (mmΒ H2Omm\space H_2OmmΒ H2​O) for the month, SNDdaySND_{day}SNDday​ is the standard normal deviate calculated for the day, and gmong_{mon}gmon​ 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 (), is the mean daily rainfall () 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

1:3.1.1 Occurrence of Wet or Dry Day

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, Pi(W/W)P_i(W/W)Pi​(W/W) or Pi(W/D)P_i(W/D)Pi​(W/D). 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.

iii
iβˆ’1,Piβˆ’1(W/W)i-1,Pi-1(W/W)iβˆ’1,Piβˆ’1(W/W)
iii
iβˆ’1,Pi(W/D)i-1,P_i(W/D)iβˆ’1,Pi​(W/D)
Pi(D/W)=1βˆ’Pi(W/W)P_i(D/W)=1-P_i(W/W)Pi​(D/W)=1βˆ’Pi​(W/W)
Pi(W/W)=1βˆ’Pi(W/D)P_i(W/W)=1-P_i(W/D)Pi​(W/W)=1βˆ’Pi​(W/D)
Pi(D/W)P_i(D/W)Pi​(D/W)
iii
iβˆ’1i-1iβˆ’1
Pi(D/D)P_i(D/D)Pi​(D/D)
iii
iβˆ’1i-1iβˆ’1

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 ()

pcpmm

pcp_ave

average number of days of precipitation in month(= PCPMM / PCPD)

pcpd

pcp_days

: standard deviation for daily precipitation in month ()

pcpstd

pcp_sd

: skew coefficient for daily precipitation in month

pcpskw

pcp_skew

SNDday=cos(6.283βˆ—rnd2)βˆ—βˆ’2ln(rnd1)SND_{day}=cos(6.283*rnd_2)*\sqrt{-2ln(rnd_1)}SNDday​=cos(6.283βˆ—rnd2​)βˆ—βˆ’2ln(rnd1​)​
rnd1rnd_1rnd1​
rnd2rnd_2rnd2​
Rday=ΞΌmonβˆ—(βˆ’ln(rnd1))rexpR_{day}=\mu_{mon}*(-ln(rnd_1))^{rexp}Rday​=ΞΌmonβ€‹βˆ—(βˆ’ln(rnd1​))rexp
RdayR_{day}Rday​
mmΒ H2Omm\space H_2OmmΒ H2​O
ΞΌmon\mu_{mon}ΞΌmon​
mmΒ H2Omm\space H_2OmmΒ H2​O
rnd1rnd_1rnd1​
rexprexprexp
rexprexprexp

Precipitation input: 'sim' for simulated or gage name

pgage

pcp

weather-sta.cli

Pi(W/D)P_i(W/D)Pi​(W/D): probability of a wet day following a dry day in month

pr_wd

wet_dry

weather-wgn.cli

Pi(W/W)P_i(W/W)Pi​(W/W): probability of a wet day following a wet day in month

pr_ww

rexprexprexp
mmΒ H2Omm\space H_2OmmΒ H2​O
ΞΌmon\mu_{mon}ΞΌmon​
Οƒmon\sigma_{mon}Οƒmon​
mmΒ H2Omm\space H_2OmmΒ H2​O
gmong_{mon}gmon​
weather-wgn.cli
weather-wgn.cli
weather-wgn.cli
weather-wgn.cli
weather-wgn.cli