arrow-left

All pages
gitbookPowered by GitBook
1 of 3

Loading...

Loading...

Loading...

1:4.2 Climate Change

The impact of global climate change on water supply is a major area of research. Climate change can be simulated with SWAT+ by manipulating the climatic input that is read into the model (precipitation, temperature, solar radiation, relative humidity, wind speed, potential evapotranspiration and weather generator parameters). A less time-consuming method is to set adjustment factors for the various climatic inputs.

SWAT+ will allow users to adjust precipitation, temperature, solar radiation, relative humidity, and carbon dioxide levels in each subbasin. The alteration of precipitation, temperature, solar radiation and relative humidity are straightforward:

Rday=Rdayβˆ—(1+adjpcp100)R_{day}=R_{day}*(1+\frac{adj_{pcp}}{100})Rday​=Rdayβ€‹βˆ—(1+100adjpcp​​) 1:4.2.1

where RdayR_{day}Rday​ is the precipitation falling in the subbasin on a given day (mm H2_22​O), and adjpcpadj_{pcp}adjpcp​ is the % change in rainfall.

1:4.2.2

where is the daily maximum temperature (Β°C), and is the change in temperature (Β°C).

1:4.2.3

where is the daily minimum temperature (Β°C), and is the change in temperature (Β°C).

1:4.2.4

where is the daily mean temperature (Β°C), and is the change in temperature (Β°C).

1:4.2.5

where is the daily solar radiation reaching the earth’s surface (), and is the change in radiation ().

1:4.2.6

where is the relative humidity for the day expressed as a fraction, and is the change in relative humidity expressed as a fraction.

SWAT+ allows the adjustment terms to vary from month to month so that the user is able to simulate seasonal changes in climatic conditions.

Changes in carbon dioxide levels impact plant growth. As carbon dioxide levels increase, plant productivity increases and plant water requirements go down. The equations used to account for the impact of carbon dioxide levels on plant water requirements are reviewed in Chapters 2:2 and 5:2. When carbon dioxide climate change effects are being simulated, the Penman-Monteith equation must be used to calculate potential evapotranspiration. This method has been modified to account for impacts on potential evapotranspiration levels.

Table 1:4-2: SWAT+ input variables that pertain to climate change.

Description
Source Name
Input Name
Input File

Chapter 1:4 Climate Customization

SWAT+ is capable of simulating a number of climate customization options. Orographic impacts on temperature and precipitation for watersheds in mountainous regions can be simulated. The model will also modify climate inputs for simulations that are looking at the impact of climatic change in a given watershed. Finally, SWAT+ allows a weather forecast period to be incorporated into a simulation to study the effects of predicted weather in a watershed.

.sub

: increase or decrease in relative humidity for month

huminc

.sub

: carbon dioxide level at the start of the simulation (ppmv)

co2

Potential evapotranspiration method

pet

pet

Tmx=Tmx+adjtmpT_{mx}=T_{mx}+adj_{tmp}Tmx​=Tmx​+adjtmp​
TmxT_{mx}Tmx​
adjtmpadj_{tmp}adjtmp​
Tmn=Tmn+adjtmpT_{mn}=T_{mn}+adj_{tmp}Tmn​=Tmn​+adjtmp​
TmnT_{mn}Tmn​
adjtmpadj_{tmp}adjtmp​
Tβ€Ύav=Tβ€Ύav+adjtmp\overline T_{av} =\overline T_{av} +adj_{tmp}Tav​=Tav​+adjtmp​
Tβ€Ύav\overline T_{av}Tav​
adjtmpadj_{tmp}adjtmp​
Hday=Hday+adjradH_{day}= H_{day}+ adj_{rad}Hday​=Hday​+adjrad​
HdayH_{day}Hday​
MJmβˆ’2MJ m^{-2}MJmβˆ’2
adjradadj_{rad}adjrad​
MJmβˆ’2dβˆ’1MJ m^{-2} d^{-1}MJmβˆ’2dβˆ’1
Rh=Rh+adjhmdR_h=R_h +adj_{hmd}Rh​=Rh​+adjhmd​
RhR_hRh​
adjhmdadj_{hmd}adjhmd​
CO2CO_2CO2​

adjpcpadj_{pcp}adjpcp​: % change in rainfall for month

rfinc

.sub

adjtmpadj_{tmp}adjtmp​: increase or decrease in temperature for month (Β°C)

tmpinc

.sub

adjradadj_{rad}adjrad​: increase or decrease in solar radiation reaching earth’s surface for month (MJ mβˆ’2^{-2}βˆ’2)

radinc

adjhmdadj_{hmd}adjhmd​
CO2CO_2CO2​
parameters.bsn
codes.bsn

1:4.1 Elevation Bands

Orographic precipitation is a significant phenomenon in certain areas of the world. To account for orographic effects on both precipitation and temperature, SWAT+ allows up to 10 elevation bands to be defined in each subbasin. Precipitation and maximum and minimum temperatures are calculated for each band as a function of the respective lapse rate and the difference between the gage elevation and the average elevation specified for the band. For precipitation,

Rband=Rday+(ELbandβˆ’ELgage)βˆ—plapsdayspcp,yrβˆ—1000R_{band}=R_{day}+(EL_{band}-EL_{gage})*\frac{plaps}{days_{pcp,yr}*1000}Rband​=Rday​+(ELbandβ€‹βˆ’ELgage​)βˆ—dayspcp,yrβ€‹βˆ—1000plaps​ when Rday>0.01R_{day}>0.01Rday​>0.01 1:4.1.1

where RbandR_{band}Rband​ is the precipitation falling in the elevation band (mm H2_22​O), RdayR_{day}Rday​ is the precipitation recorded at the gage or generated from gage data (mm H2_22​O), ELbandEL_{band}ELband​ is the mean elevation in the elevation band (m), ELgageEL_{gage}ELgage​ is the elevation at the recording gage (m), plapsplapsplaps is the precipitation lapse rate (mm H2_22​O/km), dayspcp,yrdays_{pcp,yr}dayspcp,yr​ is the average number of days of precipitation in the subbasin in a year, and 1000 is a factor needed to convert meters to kilometers. For temperature,

1:4.1.2

1:4.1.3

1:4.1.4

where is the maximum daily temperature in the elevation band (Β°C), is the minimum daily temperature in the elevation band (Β°C), is the mean daily temperature in the elevation band (Β°C), is the maximum daily temperature recorded at the gage or generated from gage data (Β°C), is the minimum daily temperature recorded at the gage or generated from gage data (Β°C), is the mean daily temperature recorded at the gage or generated from gage data (Β°C), is the mean elevation in the elevation band (m), is the elevation at the recording gage (m), is the temperature lapse rate (Β°C/km), and 1000 is a factor needed to convert meters to kilometers.

Once the precipitation and temperature values have been calculated for each elevation band in the subbasin, new average subbasin precipitation and temperature values are calculated:

1:4.1.5

1:4.1.6

1:4.1.7

1:4.1.8

where is the daily average precipitation adjusted for orographic effects (mm HO), is the daily maximum temperature adjusted for orographic effects (Β°C), is the daily minimum temperature adjusted for orographic effects (Β°C), is the daily mean temperature adjusted for orographic effects (Β°C), is the precipitation falling in elevation band (mm H2O), is the maximum daily temperature in elevation band (Β°C), is the minimum daily temperature in elevation band (Β°C), is the mean daily temperature in elevation band (Β°C), is the fraction of subbasin area within the elevation band, and is the total number of elevation bands in the subbasin.

The only processes modeled separately for each individual elevation band are the accumulation, sublimation and melting of snow. As with the initial precipitation and temperature data, after amounts of sublimation and snow melt are determined for each elevation band, subbasin average values are calculated. These average values are the values that are used in the remainder of the simulation and reported in the output files.

Table 1:4-1: SWAT+ input variables that pertain to orographic effects.

Definition
Source Name
Input Name
Input File

pcp_days

: Temperature lapse rate (Β°C/km)

tlaps

tlaps

: Daily precipitation (mm HO)

pcp

pcp

: Daily maximum temperature (Β°C)

tmpmx

tmpmax

: Daily minimum temperature (Β°C)

tmpmn

tmpmin

Tmx,band=Tmx+(ELbandβˆ’ELgage)βˆ—tlaps1000T_{mx,band}=T_{mx}+(EL_{band}-EL_{gage})*\frac{tlaps}{1000}Tmx,band​=Tmx​+(ELbandβ€‹βˆ’ELgage​)βˆ—1000tlaps​
Tmn,band=Tmn+(ELbandβˆ’ELgage)βˆ—tlaps1000T_{mn,band}=T_{mn}+(EL_{band}-EL_{gage})*\frac{tlaps}{1000}Tmn,band​=Tmn​+(ELbandβ€‹βˆ’ELgage​)βˆ—1000tlaps​
Tβ€Ύav,band=Tβ€Ύav+(ELbandβˆ’ELgage)βˆ—tlaps1000\overline T_{av,band} =\overline T_{av} +(EL_{band}-EL_{gage})*\frac{tlaps}{1000}Tav,band​=Tav​+(ELbandβ€‹βˆ’ELgage​)βˆ—1000tlaps​
Tmx,bandT_{mx,band}Tmx,band​
Tmn,bandT_{mn,band}Tmn,band​
TmxT_{mx}Tmx​
TmnT_{mn}Tmn​
ELbandEL_{band}ELband​
ELgageEL_{gage}ELgage​
tlapstlapstlaps
Rday=βˆ‘bnd=1bRbandβˆ—frbndR_{day}= \sum_{bnd=1}^{b}R_{band}*fr_{bnd}Rday​=βˆ‘bnd=1b​Rbandβ€‹βˆ—frbnd​
Tmx=βˆ‘bnd=1bTmx,bandβˆ—frbndT_{mx}=\sum_{bnd=1}^b T_{mx,band}*fr_{bnd}Tmx​=βˆ‘bnd=1b​Tmx,bandβ€‹βˆ—frbnd​
Tmn=βˆ‘bnd=1bTmn,bandβˆ—frbndT_{mn}=\sum_{bnd=1}^b T_{mn,band} *fr_{bnd}Tmn​=βˆ‘bnd=1b​Tmn,bandβ€‹βˆ—frbnd​
Tβ€Ύav=βˆ‘bnd=1bTβ€Ύav,bandβˆ—frbnd\overline T_{av}=\sum_{bnd=1}^b \overline T_{av,band}*fr_{bnd}Tav​=βˆ‘bnd=1b​Tav,bandβ€‹βˆ—frbnd​
RdayR_{day}Rday​
2_22​
TmxT_{mx}Tmx​
TmnT_{mn}Tmn​
RbandR_{band}Rband​
bndbndbnd
Tmx,bandT_{mx,band}Tmx,band​
bndbndbnd
Tmn,bandT_{mn,band}Tmn,band​
bndbndbnd
bndbndbnd
frbndfr_{bnd}frbnd​
bbb

ELgageEL_{gage}ELgage​: Elevation of temperature station (m)

elev

elev

.tmp

plapsplapsplaps: Precipitation lapse rate (mm H2_22​O/km)

plaps

plaps

parameters.bsn

average number of days of precipitation in month , dayspcp,yr=βˆ‘days_{pcp,yr}=\sumdayspcp,yr​=βˆ‘(PCPD(mon)) for a subbasin

pcpd

tlapstlapstlaps
RdayR_{day}Rday​
2_22​
TmxT_{mx}Tmx​
TmnT_{mn}Tmn​
weather-wgn.cli
parameters.bsn
.pcp
.tmp
.tmp