Radiosity Approach for the Shortwave Surface Radiation Balance in Complex Terrain

The influence of topography on the radiation balance in complex terrain has so far been investigated either with very simple or very sophisticated approaches that are limited, respectively, by an uncontrolled spatial representation of radiative fluxes or heavy computational efforts. To bridge this gap in complexity, this paper proposes the radiosity approach, well known in computer graphics, to study anisotropic reflections of radiation in complex terrain. To this end the radiosity equation is rederived in the context of three-dimensional radiative transfer. The discretized equation is solved by means of an adapted version of progressive refinement iteration. To systematically study terrain effects, the geometrical disorder provided by the topography is considered in its simplest approximation by Gaussian random fields. These model topographies capture the most important length scales of complex terrain, namely a typical elevation and a typical valley width via the variance and the correlation length of the field, respectively. The mean reflected radiation is computed as a function of these length scales and sun elevation, thereby explicitly addressing finite system sizes and grid resolutions. A comparison with an isotropic parameterization of terrain reflections reveals that mean values are similar whereas spatial distributions vary remarkably. It is also shown that the mean reflected radiation in real topography is reasonably well characterized by the Gaussian approximation. As a final application of the method, the effective albedo of a topography is shown to vary with sun elevation and domain-averaged albedo, leading to albedo differences up to 0.025.