1:1.2.5 Daily Net Radiation
Net radiation requires the determination of both incoming and reflected short-wave radiation and net long-wave or thermal radiation. Expressing net radiation in terms of the net short-wave and long-wave components gives:
Hnet=Hday↓−α∗Hday↑+HL↓−HL↑ 1:1.2.11
or
Hnet=(1−α)Hday+Hb 1:1.2.12
where Hnet is the net radiation (MJm−2d−1), Hday is the short-wave solar radiation reaching the ground (MJm−2d−1), is the short-wave reflectance or albedo, HL is the long-wave radiation (MJm−2d−1), Hb is the net incoming long-wave radiation (MJm−2d−1) and the arrows indicate the direction of the radiation flux.
1:1.2.5.1 Net Short-Wave Radiation
Net short-wave radiation is defined as (1−α)Hday. SWAT+ calculates a daily value for albedo as a function of the soil type, plant cover, and snow cover. When the snow water equivalent is greater than 0.5 mm,
α=0.8 1:1.2.13
When the snow water equivalent is less than 0.5 mm and no plants are growing in the HRU,
α=αsoil 1:1.2.14
where αsoil is the soil albedo. When plants are growing and the snow water equivalent is less than 0.5 mm,
α=αplant(1−covsol)+αsoilcovsol 1:1.2.15
where αplant is the plant albedo (set at 0.23), and covsol is the soil cover index. The soil cover index is calculated
covsol=exp(−5.0X10−5∗CV) 1:1.2.16
where CV is the aboveground biomass and residue (kgha−1).
1:1.2.5.2 Net Long-Wave Radiation
Long-wave radiation is emitted from an object according to the radiation law:
HR=εσTK4 1:1.2.17
where HR is the radiant energy (MJm−2d−1), ε is the emissivity, σ is the Stefan-Boltzmann constant (4.90310−9MJm−2K−4d−1), and TK is the mean air temperature in Kelvin (273.15 + °C). Net long-wave radiation is calculated using a modified form of equation 1:1.2.17 (Jensen et al., 1990):
Hb=fcld(εa−εvs)σTK4 1:1.2.18
where Hb is the net long-wave radiation (MJm−2d−1), fcld is a factor to adjust for cloud cover, εa is the atmospheric emittance, and εvs is the vegetative or soil emittance.
Wright and Jensen (1972) developed the following expression for the cloud cover adjustment factor, fcld:
fcld=aHMXHday−b 1:1.2.19
where a and b are constants, Hday is the solar radiation reaching the ground surface on a given day (MJm−2d−1), and HMX is the maximum possible solar radiation to reach the ground surface on a given day (MJm−2d−1).
The two emittances in equation 1:1.2.18 may be combined into a single term, the net emittance ε′. The net emittance is calculated using an equation developed by Brunt (1932):
ε′=εa−εvs=−(a1+b1(e)) 1:1.2.20
where a1 and b1 are constants and e is the vapor pressure on a given day (kPa). The calculation of e is given in Chapter 1:2. Combining equations 1:1.2.18, 1:1.2.19, and 1:1.2.20 results in a general equation for net long-wave radiation:
Hb=−[aHMXHday−b][a1+b1(e)]σTk4 1:1.2.21
Experimental values for the coefficients a,b,a1, and b1 are presented in Table 1:1.3. The default equation in SWAT+ uses coefficient values proposed by Doorenbos and Pruitt (1977):
Hb=−[0.9HMXHday+0.1][0.34−0.139(e)]σTk4 1:1.2.22
Table 1:1-3: Experimental coefficients for net long-wave radiation equations (from Jensen et al., 1990).
Davis, California
(1.35,
-0.35)
(0.35,
-0.145)
Southern Idaho
(1.22,
-0.18)
(0.325,
-0.139)
England
not available
not available
(0.47,
-0.206)
England
not available
not available
(0.44,
-0.253)
Australia
not available
not available
(0.35,
-0.133)
General
(1.2
-0.2)
(0.39,
-0.158)
General-humid areas
(1.0
0.0)
General-semihumid areas
(1.1
-0.1)
Table 1:1-4: SWAT+ input variables used in net radiation calculations.
Last updated