The procedure used to generate daily values for maximum temperature, minimum temperature and solar radiation (Richardson, 1981; Richardson and Wright, 1984) is based on the weakly stationary generating process presented by Matalas (1967).

For each month, users provide the maximum half-hour rain observed over the entire period of record. These extreme values are used to calculate representative monthly maximum half-hour rainfall fractions.

Prior to calculating the representative maximum half-hour rainfall fraction for each month, the extreme half-hour rainfall values are smoothed by calculating three month average values:

$R_{0.5sm(mon)}=\frac{R_{0.5x(mon-1)}+R_{0.5x(mon)}+R_{0.5x(mon+1)}}{3}$ 1:3.2.1

where $R_{0.5sm(mon)}$ is the smoothed maximum half-hour rainfall for a given month (mm H$_2$O) and $R_{0.5x}$ is the extreme maximum half-hour rainfall for the specified month (mm H$_2$O). Once the smoothed maximum half-hour rainfall is known, the representative half-hour rainfall fraction is calculated using the equation:

1:3.2.2

where is the average half-hour rainfall fraction for the month, is an adjustment factor, is the smoothed half-hour rainfall amount for the month (mm HO), is the mean daily rainfall (mm HO) for the month, is the number of years of rainfall data used to obtain values for monthly extreme half-hour rainfalls, and are the number of wet days in the month. The adjustment factor is included to allow users to modify estimations of half-hour rainfall fractions and peak flow rates for runoff.

The rainfall intensity distribution given in equation 1:3.3.1 can be normalized to eliminate units. To do this, all time values are divided, or normalized, by the storm duration and all intensity values are normalized by the average storm intensity. For example,

1:3.3.2

1:3.3.3

where the normalized rainfall intensity at time , is the rainfall intensity at time T (mm/hr), is the average storm rainfall intensity (mm/hr), is the time during the storm expressed as a fraction of the total storm duration (0.0-1.0), is the time since the beginning of the storm (hr), and is the duration of the storm (hr).

Maximum half-hour rainfall is required by SWAT+ to calculate the peak flow rate for runoff. When daily precipitation data are used by the model, the maximum half-hour rainfall may be calculated from a triangular distribution using monthly maximum half-hour rainfall data or the user may choose to use the monthly maximum half-hour rainfall for all days in the month. The maximum half-hour rainfall is calculated only on days where surface runoff has been generated.

The normalized time to peak intensity is calculated by SWAT+ using a triangular distribution. The triangular distribution used to generate the normalized time to peak intensity requires four inputs: average time to peak intensity expressed as a fraction of total storm duration $(\hat t_{peakM})$, maximum time to peak intensity expressed as a fraction of total storm duration $(\hat t_{peakU})$, minimum time to peak intensity expressed as a fraction of total storm duration $(\hat t_{peakL})$and a random number between 0.0 and 1.0.

The maximum time to peak intensity, or upper limit of the triangular distribution, is set at 0.95. The minimum time to peak intensity, or lower limit of the triangular distribution is set at 0.05. The mean time to peak intensity is set at 0.25.

The triangular distribution uses one of two sets of equations to generate a normalized peak intensity for the day. If $rnd_1\le[\frac{\hat t_{peakM} - \hat t_{peakL}}{\hat t_{peakU} - \hat t_{peakL}}]$ then

$\hat t_{peak}=\hat t_{peakM} * \frac{\hat t_{peakL}+[rnd_1*(\hat t_{peakU}-\hat t_{peakL})*(\hat t_{peakM}-\hat t_{peakL})]^{0.5}}{\hat t_{peak,mean}}$ 1:3.3.9

If then

1:3.3.10

where is the time from the beginning of the storm till the peak intensity expressed as a fraction of the total storm duration (0.0-1.0),is the average time to peak intensity expressed as a fraction of storm duration, is a random number generated by the model each day, is the minimum time to peak intensity that can be generated,is the maximum time to peak intensity that can be generated, and is the mean of and .

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.

For simulations where the timing of rainfall within the day is required, the daily rainfall value must be partitioned into shorter time increments. The method used in SWAT+ to disaggregate storm data was taken from CLIGEN (Nicks et al., 1995).

A double exponential function is used to represent the intensity patterns within a storm. With the double exponential distribution, rainfall intensity exponentially increases with time to a maximum, or peak, intensity. Once the peak intensity is reached, the rainfall intensity exponentially decreases with time until the end of the storm.

The exponential equations governing rainfall intensity during a storm event are:

1:3.3.1

,

The user has the option of using the monthly maximum half-hour rainfall for all days in the month or generating a daily value. The variable ISED_DET in the basin input file (.bsn) defines which option the user prefers. The randomness of the triangular distribution used to generated daily values can cause the maximum half-hour rainfall value to jump around. For small plots or microwatersheds in particular, the variability of the triangular distribution is unrealistic.

The triangular distribution used to generate the maximum half-hour rainfall fraction requires four inputs: average monthly half-hour rainfall fraction, maximum value for half-hour rainfall fraction allowed in month, minimum value for half-hour rainfall fraction allowed in month, and a random number between 0.0 and 1.0.

The maximum half-hour rainfall fraction, or upper limit of the triangular distribution, is calculated from the daily amount of rainfall with the equation:

1:3.2.3

Residuals for maximum temperature, minimum temperature and solar radiation are required for calculation of daily values. The residuals must be serially correlated and cross-correlated with the correlations being constant at all locations. The equation used to calculate residuals is:

1:3.4.1

where is a 3 × 1 matrix for day whose elements are residuals of maximum temperature (), minimum temperature () and solar radiation (), ) is a 3 × 1 matrix of the previous day’s residuals, is a 3 × 1 matrix of independent random components, and and are 3 × 3 matrices whose elements are defined such that the new sequences have the desired serial-correlation and cross-correlation coefficients. The and matrices are given by

The volume of rain is related to rainfall intensity by:

1:3.3.11

where is the amount of rain that has fallen at time (mm HO) and is the rainfall intensity at time (mm/hr).

Using the definition for rainfall intensity given in equation 1:3.3.1, equation 1:3.3.11 can be integrated to get:

Maximum temperature and solar radiation will be lower on overcast days than on clear days. To incorporate the influence of wet/dry days on generated values of maximum temperature and solar radiation, the average daily maximum temperature, , and average daily solar radiation, , in equations 1:3.4.10 and 1:3.4.12 are adjusted for wet or dry conditions.

The continuity equation relates average daily maximum temperature adjusted for wet or dry conditions to the average daily maximum temperature for the month:

1:3.4.14

where is the average daily maximum temperature for the month (°C), are the total number of days in the month,

Relative humidity is required by SWAT+ when the Penman-Monteith equation is used to calculate potential evapotranspiration. It is also used to calculate the vapor pressure deficit on plant growth. Daily average relative humidity values are calculated from a triangular distribution using average monthly relative humidity. This method was developed by J.R. Williams for the EPIC model (Sharpley and Williams, 1990).

The continuity equation relates average daily solar radiation adjusted for wet or dry conditions to the average daily solar radiation for the month:

$\mu rad_{mon}*days_{tot}=\mu Wrad_{mon}*days_{wet}+\mu Drad_{mon}*days_{dry}$ 1:3.4.19

where $\mu rad_{mon}$ is the average daily solar radiation for the month (MJ m$^{-2}$), $days_{tot}$ are the total number of days in the month, $\mu Wrad_{mon}$ is the average daily solar radiation of the month on wet days (MJ m$^{-2}$), $days_{wet}$ are the number of wet days in the month, $\mu Drad_{mon}$is the average daily solar radiation of the month on dry days (MJ m$^{-2}$), and $days_{dry}$ are the number of dry days in the month.

The wet day average solar radiation is assumed to be less than the dry day average solar radiation by some fraction:

1:3.4.20

where is the average daily solar radiation of the month on wet days (MJ m), is the average daily solar radiation of the month on dry days (MJ m), and is a scaling factor that controls the degree of deviation in solar radiation caused by the presence or absence of precipitation. The scaling factor, , is set to 0.5 in SWAT+.

To calculate the dry day average solar radiation, equations 1:3.4.19 and 1:3.4.20 are combined and solved for :

1:3.4.21

Incorporating the modified values into equation 1:3.4.12, SWAT+ calculated the solar radiation on a wet day using the equation:

1:3.4.22

and the solar radiation on a dry day using the equation:

1:3.4.23

Table 1:3-4: SWAT+ input variables that pertain to generation of temperature and solar radiation.

The triangular distribution used to generate daily relative humidity values requires four inputs: mean monthly relative humidity, maximum relative humidity value allowed in month, minimum relative humidity value allowed in month, and a random number between 0.0 and 1.0.

The maximum relative humidity value, or upper limit of the triangular distribution, is calculated from the mean monthly relative humidity with the equation:

$R_{hUmon}=R_{hmon}+(1-R_{hmon})*exp(R_{hmon}-1)$ 1:3.5.4

where $R_{hUmon}$ is the largest relative humidity value that can be generated on a given day in the month, and $R_{hmon}$ is the average relative humidity for the month.

The minimum relative humidity value, or lower limit of the triangular distribution, is calculated from the mean monthly relative humidity with the equation:

1:3.5.5

where is the smallest relative humidity value that can be generated on a given day in the month, and is the average relative humidity for the month.

The triangular distribution uses one of two sets of equations to generate a relative humidity value for the day. If then

1:3.5.6

If then

1:3.5.7

where is the average relative humidity calculated for the day, is a random number generated by the model each day, is the average relative humidity for the month, is the smallest relative humidity value that can be generated on a given day in the month, is the largest relative humidity value that can be generated on a given day in the month, and is the mean of and .

The daily generated values are determined by multiplying the residual elements generated with equation 1:3.4.1 by the monthly standard deviation and adding the monthly average value.

$T_{mx}=\mu mx_{mon} + \chi_i(1)*\sigma mx_{mon}$ 1:3.4.10

$T_{mn}=\mu mn_{mon} + \chi_i(2)*\sigma mn_{mon}$ 1:3.4.11

$H_{day}=\mu rad_{mon} + \chi_i(3)*\sigma rad_{mon}$ 1:3.4.12

where $T_{mx}$ is the maximum temperature for the day (°C), $\mu mx_{mon}$ is the average daily maximum temperature for the month (°C), $\chi_i(1)$ is the residual for maximum temperature on the given day, $\sigma mx_{mon}$ is the standard deviation for daily maximum temperature during the month (°C), $T_{mn}$ is the minimum temperature for the day (°C), $\mu mn_{mon}$ is the average daily minimum temperature for the month (°C), is the residual for minimum temperature on the given day, is the standard deviation for daily minimum temperature during the month (°C), is the solar radiation for the day (MJ m), is the average daily solar radiation for the month (MJ m), is the residual for solar radiation on the given day, and is the standard deviation for daily solar radiation during the month (MJ m).

The user is required to input standard deviation for maximum and minimum temperature. For solar radiation the standard deviation is estimated as ¼ of the difference between the extreme and mean value for each month.

1:3.4.13

where is the standard deviation for daily solar radiation during the month (MJ m), is the maximum solar radiation that can reach the earth’s surface on a given day (MJ m), and is the average daily solar radiation for the month (MJ m).

Relative humidity is defined as the ratio of the actual vapor pressure to the saturation vapor pressure at a given temperature:

$R_{hmon}=\frac{e_{mon}}{e^o_{mon}}$ 1:3.5.1

where $R_{hmon}$ is the average relative humidity for the month, $e_{mon}$ is the actual vapor pressure at the mean monthly temperature (kPa), and $e^o_{mon}$ is the saturation vapor pressure at the mean monthly temperature (kPa). The saturation vapor pressure,$e^o_{mon}$ , is related to the mean monthly air temperature with the equation:

$e^o_{mon}=exp[\frac{16.78*\mu tmp_{mon}-116.9}{\mu tmp_{mon}+237.3}]$ 1:3.5.2

where $e^o_{mon}$ is the saturation vapor pressure at the mean monthly temperature (kPa), and $\mu tmp_{mon}$ is the mean air temperature for the month (°C). The mean air temperature for the month is calculated by averaging the mean maximum monthly temperature, , and the mean minimum monthly temperature, .

The dew point temperature is the temperature at which the actual vapor pressure present in the atmosphere is equal to the saturation vapor pressure. Therefore, by substituting the dew point temperature in place of the average monthly temperature in equation 1:3.5.2, the actual vapor pressure may be calculated:

1:3.5.3

where is the actual vapor pressure at the mean month temperature (kPa), and is the average dew point temperature for the month (°C).

Wind speed is required by SWAT+ when the Penman-Monteith equation is used to calculate potential evapotranspiration. Mean daily wind speed is generated in SWAT+ using a modified exponential equation:

$\mu _{10m}=\mu wnd_{mon}*(-ln(rnd_1))^{0.3}$ 1.3.6.1

where $\mu _{10m}$ is the mean wind speed for the day (m s$^{-1}$), $\mu wnd_{mon}$ is the average wind speed for the month (m s$^{-1}$), and $rnd_1$ is a random number between 0.0 and 1.0.

Table 1:3-6: SWAT+ input variables that pertain to generation of wind speed.

To incorporate the effect of clear and overcast weather on generated values of relative humidity, monthly average relative humidity values can be adjusted for wet or dry conditions.

The continuity equation relates average relative humidity adjusted for wet or dry conditions to the average relative humidity for the month:

$R_{hmon}*days_{tot}=R_{hWmon}*days_{wet}+R_{hDmon}*days_{dry}$ 1:3.5.8

where $R_{hmon}$ is the average relative humidity for the month, $days_{tot}$ are the total number of days in the month, $R_{hWmon}$ is the average relative humidity for the month on wet days, $days_{wet}$ are the number of wet days in the month, $R_{hDmon}$ is the average relative humidity of the month on dry days, and $days_{dry}$ are the number of dry days in the month.

The wet day average relative humidity is assumed to be greater than the dry day average relative humidity by some fraction:

1:3.5.9

where is the average relative humidity of the month on wet days, is the average relative humidity of the month on dry days, and is a scaling factor that controls the degree of deviation in relative humidity caused by the presence or absence of precipitation. The scaling factor, , is set to 0.9 in SWAT+.

To calculate the dry day relative humidity, equations 1:3.5.8 and 1:3.5.9 are combined and solved for :

1:3.5.10

To reflect the impact of wet or dry conditions, SWAT+ will replace with on wet days or on dry days in equations 1:3.5.4 through 1:3.5.7.

Table 1:3-5: SWAT+ input variables that pertain to generation of relative humidity.

SWAT+ requires daily values of precipitation, maximum and minimum temperature, solar radiation, relative humidity and wind speed. The user may choose to read these inputs from a file or generate the values using monthly average data summarized over a number of years.

SWAT+ includes the WXGEN weather generator model (Sharpley and Williams, 1990) to generate climatic data or to fill in gaps in measured records. This weather generator was developed for the contiguous U.S. If the user prefers a different weather generator, daily input values for the different weather parameters may be generated with an alternative model and formatted for input to SWAT+.

The occurrence of rain on a given day has a major impact on relative humidity, temperature and solar radiation for the day. The weather generator first independently generates precipitation for the day. Once the total amount of rainfall for the day is generated, the distribution of rainfall within the day is computed if the Green & Ampt method is used for infiltration. Maximum temperature, minimum temperature, solar radiation and relative humidity are then generated based on the presence or absence of rain for the day. Finally, wind speed is generated independently.

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.