Effects of Shading and Aspect

From Gsshawiki
Jump to: navigation, search

Topographic Shading
For every cell within the GSSHA model a search is ran every 2 weeks to determine which hours of the day the cell is shaded from surrounding topography. The hourly Solar Azimuth Angle (SAA) and Solar Elevation Angle (SEA) in combination with basic geometry on the structured grid are used to determine if any cell within the GSSHA domain blocks the direct line between the solar location and the cell. If the direct line is blocked during an hour, Kt,x=0.0 for that hour (simulating complete shading), otherwise Kt,x=1.0 for that hour (simulating no shading). Kt,x is a reduction factor that is used to reduce the amount of shortwave radiation affecting each cell. Having Kt,x be a value of 0.0 or 1.0 is the same assumption that is made in the GEOTOP model (Zanotti, Endrizzi et al. 2004).

Aspect Angle
Due to several factors including location of the sun and the slope of land surface, the solar radiation emitted from the sun often does not hit the land surface on a perpendicular plane, resulting in less radiation per unit area (Bras 1990). To account for this, a reduction in shortwave radiation based on the aspect angle of the land surface to the location of the sun (Ks,x) is calculated hourly for each cell within the model. In order to complete this, the location of the sun is derived in three-dimensional space (sx, sy, sz) on a unit-sphere. For each cell a three-dimensional vector (cx, cy, cz) on a unit-sphere is then created that is perpendicular to the cell surface. A dot product can be used between the two vectors to calculate the solar-surface incident angle (SSIA).

An adjusted solar elevation angle (SEAadj) can then be computed, which takes into consideration the slope of the land surface within the cell. Ks,x is then a function of SEAadj.