Note
Need help? Please let us know in the UMEP Community.
4.5. Solar Radiation: Solar Energy on Building Envelopes (SEBE)
- Contributor:
Name
Institution
Fredrik Lindberg
Gothenburg
Dag Wästberg
Ramböll
- Introduction
The SEBE plugin (Solar Energy on Building Envelopes) can be used to calculate pixel wise potential solar energy using ground and building digital surface models (DSM). SEBE is also able to estimate irradiance on building walls. Optionally, vegetation DSMs could also be used.
The total irradiance for a roof pixel (R) on a DSM is calculated by summing the direct, diffuse and reflected radiation such as:
\[R={\sum_{i=0}^p}[(I{\omega}S + DS + G(1-S){\alpha})]\]where p is the number of patches on the hemisphere. I is the incidence direct radiation, D is diffuse radiation and G is the global radiation originating from the ith patch. α is the surface albedo and S is the shadow calculated to each pixel. ω is the Sun incidence angle given by:
\[\begin{split}{\omega}=sin(DSM_{slope}) * cos(Sun_{altitude}) * cos(Sun_{azimuth}-DSM_{aspect}) + \\ cos(DSM_{slope}) * sin(Sun_{altitude})\end{split}\]The methodology and evaluation is described in detail in Lindberg et al. (2015).
- Dialog box
- Consists of:
top section where input data is specified
bottom section for specifying the output and for running the calculations
Fig. 4.11 The dialog for the SEBE model
- Building and Ground DSM
A DSM consisting of ground and building heights. This dataset also decides the latitude and longitude used for the calculation of the Sun position.
- Vegetation Canopy DSM
A DSM consisting of pixels with vegetation heights above ground. Pixels where no vegetation is present should be set to zero.
- Vegetation Trunk Zone DSM
A DSM (geoTIFF) consisting of pixels with vegetation trunk zone heights above ground. Pixels where no vegetation is present should be set to zero.
- Use Vegetation DSMs
Tick this box if you want to include vegetation (trees and bushes) into the analysis.
- Trunk Zone DSM Exist
Tick this box if a trunk zone DSM already exist.
- Transmissivity of Light Through Vegetation (%)
Percentage of light that is penetrating through vegetation. Default value is set to 3 % according to Konarska et al. (2013).
- Percent of Canopy Height
If a trunk zone vegetation DSM is absent, this can be generated based on the height of the Canopy DSM. The default percentage is set to 25%.
- Wall Height Raster
A raster of the same size and extent as the ground and building DSM including information of the wall pixels and its height in meters above ground should be specified here. Non wall pixels should be set to zero. This raster is used to estimate irradiance on building walls and can be generated using the Wall Height and Aspect plugin located at UMEP -> Pre-processing -> Urban Geometry -> Wall Height and Aspect.
- Wall Aspect Raster
A raster of the same size and extent as the ground and building DSM including information of the wall pixels and its aspect, i.e. angle, should be specified here. For example a wall facing towards the south has a value of 180°. Non wall pixels should be set to zero. This raster is used to estimate irradiance on building walls and can be generated using the Wall Height and Aspect plugin located at UMEP -> Pre-processing -> Urban Geometry -> Wall Height and Aspect.
- Albedo
This parameter specifies the reflectivity of shortwave radiation of all surfaces (ground, roofs, walls and vegetation). It should be a value between 0 and 1. The default value is set to 0.15.
- UTC Offset (Hours)
Time zone needs to be specified. Positive numbers increase when moving east (e.g. Stockholm UTC +1). This is related to the meteorological forcing data so if ERA5 data is used, UTC should be equal to zero.
- Estimate Diffuse and Direct Shortwave Components from Global Radiation
Tick this if only global radiation is present. Diffuse and direct shortwave components will then be estimated from global radiation based on the statistical model presented by Reindl et al. (1990). If air temperature and relative humidity is present, the statistical model will perform better but it is able to estimate the components using only global shortwave radiation.
- Input Meteorological File
Input meteorological data specifically formatted to be used in UMEP. This specific format can be created using UMEP -> Pre-processing -> Meteorological data -> Prepare existing data. A dataset with hourly time resolution should be used for SEBE, preferably at least one year in length. The time should preferably be in LST for the specific location to be modelled (ERA5 data, UTC=0). Multiple years can also be used to improve the model outcome. Model output is dependent on the meteorological input data so if a short dataset is used, potential solar energy would be valid for that particular time period only.
Mandatory data is global shortwave radiation, but the model will perform best if also diffuse and direct components are available.
The direct radiation component used as input in the SOLWEIG model is not the direct shortwave radiation on a horizontal surface but the beam/direct irradiance on a plane always normal to sun rays. Hence, the relationship between global radiation and the two separate components are:
Global radiation = direct radiation * sin(h) + diffuse radiation
where h is the sun altitude. Since diffuse and direct components of short wave radiation is not common data, it is also possible to calculate diffuse and direct shortwave radiation (see above).
- Save Sky Irradiance Distribution
When the box is ticked, it is possible to save the radiation distribution from the sky vault calculated from the meteorological file. SEBE first distributes the radiation on 145 sky patches on the sky vault and then generates shadows on the DSMs based on these patches, i.e. the core loop in the model iterates 145 times. For more detailed information on this, see Lindberg et al. (2015).
- Output Folder
A specified folder where result will be saved should be specified here. One raster showing irradiance on ground and building roofs named Energyyearroof.tif is saved as well as a text file of wall irradiance (Energyyearwall.txt). Also, the ground and building DSM is saved in the output folder to be used later in a SEBE visualization plugin (UMEP -> Post-processing -> Solar Energy -> SEBE (Visualisation)).
- Run
This starts the calculations.
- Add Roof and Ground Irradiance Result Raster to Project
If this is ticked, Energyyearroof.tif will be loaded into to the map canvas.
- Close
This button closes the plugin.
- Output
As mentioned earlier, three mandatory datasets are saved if the model runs successfully. The geoTIFF Energyyearroof.tif show pixel wise total irradiance in kWh. Energyyearwall.txt show total wall irradiance for each wall column. The Energyyearwall.txt is formatted in the following way: first and second column is the relative row and column number from upper left corner of the modelled grid. The following columns are irradiance for each wall voxel starting from the ground and moving upwards as going right in each row. If zero values are found (especially at the end of the row) that means that the wall column has reached its maximum height. The column voxel is decided based on the pixel resolution of the input data. Also, the ground and building DSM is saved in the output folder for later use. If the vegetation DSMs were added, one additional file (Vegetationdata.txt) including information of vegetation height and location, is also saved. This file is also used in the SEBE visualization plugin.
Example of input data and result
Fig. 4.12 Input DSM (left) and irradiance image (right) in Gothenburg using data from 1977.
- Remarks
All DSMs need to have the same extent and pixel resolution.
This plugin is computationally intensive i.e. large grids will take a lot of time and very large grids will not be possible to use. Large grids e.g. larger than 4000000 pixels should be tiled before.
- References
Konarska J, Lindberg F, Larsson A, Thorsson S, Holmer B 2013. Transmissivity of solar radiation through crowns of single urban trees—application for outdoor thermal comfort modelling. Theoret. Appl. Climatol., 1–14 Link to Paper
Lindberg, F., Jonsson, P. & Honjo, T. and Wästberg, D. (2015) Solar energy on building envelopes - 3D modelling in a 2D environment. Solar Energy. 115 (2015) 369–378 Link to Paper
Reindl DT, Beckman WA, Duffie JA (1990) Diffuse fraction correlation. Sol Energy 45:1–7. Link to paper