Solar radiation analysis equations
Release 9.3 Last modified January 13, 2012	Print all topics in : "Solar Radiation analysis"

Note: This topic was updated for 9.3.1.

The solar radiation analysis tools calculate insolation across a landscape or for specific locations, based on methods from the hemispherical viewshed algorithm developed by Rich et al. (Rich 1990, Rich et al. 1994), as further developed by Fu and Rich (2000, 2002).

The total amount of radiation calculated for a particular location or area is given as global radiation. The calculation of direct, diffuse, and global insolation are repeated for each feature location or every location on the topographic surface producing insolation maps for an entire geographic area.

Global radiation calculation

Global radiation (Global_tot) is calculated as the sum of direct (Dir_tot) and diffuse (Dif_tot) radiation of all sunmap and skymap sectors, respectively.

Globaltot = Dirtot + Diftot

Direct radiation calculation

Total direct insolation (Dir_tot) for a given location is the sum of the direct insolation (Dir_θ,α) from all sunmap sectors:

Dirtot =  ΣDirθ,α			(1)

The direct insolation from the sunmap sector (Dir_θ,α) with a centroid at zenith angle (θ) and azimuth angle (α) is calculated using the following equation:

Dirθ,α = SConst * βm(θ) * SunDurθ,α * SunGapθ,α * cos(AngInθ,α)	(2)

where:

S_Const is the solar flux outside the atmosphere at the mean earth-sun distance, known as solar constant. The solar constant used in the analysis is 1367 W/m². This is consistent with the World Radiation Center (WRC) solar constant.

β is transmisivity of the atmosphere (averaged over all wavelengths) for the shortest path (in the direction of the zenith);

m(θ) is the relative optical path length, measured as a proportion relative to the zenith path length (see equation 3 below).

SunDur_θ,α is the time duration represented by the sky sector. For most sectors, it is equal to the day interval (for example, a month) multiplied by the hour interval (for example, a half hour). For partial sectors (near the horizon), the duration is calculated using spherical geometry;

SunGap_θ,α is the gap fraction for the sunmap sector;

AngIn_θ,α is the angle of incidence between the centroid of the sky sector and the axis normal to the surface (see equation 4 below).

Relative optical length (m(θ)) is determined by the solar zenith angle and elevation above sea level. For zenith angles less than 80°, it can be calculated using the following equation:

m(θ)  = EXP(-0. 000118 * Elev - 1. 638 * 10-9 * Elev2) / cos(θ)	(3)

where:

θ is the solar zenith angle;

Elev is elevation above sea level in meters.

The effect of surface orientation is taken into account by multiplying by the cosine of the angle of incidence. Angle of incidence (AngInSky_θ,α) between the intercepting surface and a given sky sector with a centroid at zenith angle and azimuth angle is calculated using the following equation:

AngInθ,α = acos[Cos(θ)*Cos(Gz)+Sin(θ)*Sin(Gz)*Cos(α-Ga)] 	(4)

where:

G_z is the surface zenith angle;

G_a is the surface azimuth angle.

(For zenith angles greater than 80° refraction is important).

Diffuse radiation calculation

For each sky sector, the diffuse radiation at its centroid (Dif) is calculated, integrated over the time interval, and corrected by the gap fraction and angle of incidence using the following equation:

Difnθ,α  = Rglb * Pdif * Dur * SkyGapθ,α * Weightθ,α * cos(AngInθ,α)	(5)

where:

R_glb is the global normal radiation (see equation 6 below);

P_dif is the proportion of global normal radiation flux that is diffused. Typically it is approximately 0.2 for very clear sky conditions and 0.7 for very cloudy sky conditions;

Dur is the time interval for analysis;

SkyGap_θ,α is the gap fraction (proportion of visible sky) for the sky sector;

Weight_θ,α is the proportion of diffuse radiation originating in a given sky sector relative to all sectors (see equation 7 and 8 below);

AngIn_θ,α is the angle of incidence between the centroid of the sky sector and the intercepting surface.

The global normal radiation (R_glb) can be calculated by summing the direct radiation from every sector (including obstructed sectors) without correction for angle of incidence, then correcting for proportion of direct radiation, which equals to 1- P_dif:

Rglb = (SConst Σ(βm(θ)))/ (1 - Pdif)    	(6)

For the uniform sky diffuse model, Weight_θ,α is calculated as follows:

Weightθ,α = (cosθ2- cos θ1) / Divazi      	(7)

where:

θ₁ and θ₂ are the bounding zenith angles of the sky sector;

Div_azi is the number of azimuthal divisions in the skymap.

For the standard overcast sky model, Weight_θ,α is calculated as follows:

Weightθ,α  = (2cosθ2+ cos2θ2- 2cosθ1- cos2θ1) / 4 * Divazi 	(8)

Total diffuse solar radiation for the location (Dif_tot) is calculated as the sum of the diffuse solar radiation (Dif) from all the skymap sectors:

Diftot =  ΣDifθ,α 	(9)

Please visit the Feedback page to comment or give suggestions on ArcGIS Desktop Help.
Copyright © Environmental Systems Research Institute, Inc.