Berechnung der Sonneneinstrahlung
Mit Spatial Analyst-Lizenz verfügbar.
Mit den Werkzeugen für die Analyse der Sonneneinstrahlung wird die Sonneneinstrahlung in einer Landschaft oder an bestimmten Orten anhand von Methoden berechnet, die aus dem von Rich et al. (Rich 1990, Rich et al. 1994) entwickelten und von Fu und Rich (2000, 2002) weiterentwickelten hemisphärischen Sichtfeldalgorithmus stammen.
Die Gesamtmenge der für einen bestimmten Ort oder eine bestimmte Fläche berechneten Strahlung wird als globale Strahlung angegeben. Die Berechnung der direkten, diffusen und globalen Sonneneinstrahlung wird für jede Feature-Position oder jeden Ort auf der topografischen Oberfläche wiederholt, womit Sonneneinstrahlungskarten für einen vollständigen geographischen Bereich erzeugt werden.
Die folgenden Abschnitte enthalten Informationen zu den zur Berechnung der Sonneneinstrahlung verwendeten Gleichungen.
Berechnung der globalen Strahlung
Die globale Strahlung (\(Global_{tot}\)) wird als Summe der direkten (\(Dir_{tot}\)) und diffusen (\(Dif_{tot}\)) Strahlung aller Sonnenkarten- bzw. Himmelskartensektoren berechnet.
\(Global_{tot} = Dir_{tot} + Dif_{tot}\)
Direkte Sonneneinstrahlung
Die gesamte direkte Sonneneinstrahlung (\(Dir_{tot}\)) für einen bestimmten Ort wird als Summe der direkten Sonneneinstrahlung (\(Dir_{θ,α}\)) aus allen Sonnenkartensektoren berechnet:
\(\large Dir_{tot} = Σ \space Dir_{θ,α} \space \space \space \normalsize (1)\)
Die direkte Sonneneinstrahlung aus dem Sonnenkartensektor (\(Dir_{θ,α}\)) mit einem durch Zenitwinkel (\(θ\)) und Azimutwinkel (\(α\)) bestimmten Schwerpunkt wird anhand der folgenden Gleichung berechnet:
\(\large Dir_{θ,α} = S_{Const} * β ^{m(θ)} * SunDur_{θ,α} * SunGap_{θ,α} * cos(AngIn_{θ,α} ) \space \space \space \normalsize (2)\)
Dabei gilt:
\(S_{Const}\) ist der solare Strahlungsfluss außerhalb der Atmosphäre bei mittlerem Abstand zwischen Erde und Sonne, als Solarkonstante bezeichnet. Die in der Analyse verwendete Solarkonstante beträgt 1367 W/m2. Dies entspricht der vom World Radiation Center (WRC) verwendeten Solarkonstante.
β ist die Transmissivität der Atmosphäre (Durchschnitt über alle Wellenlängen) auf dem kürzesten Weg (in Richtung des Zenits).
\(m(θ)\) ist die relative Länge des optischen Pfades, gemessen als Verhältnis zur Pfadlänge zum Zenit (siehe Gleichung 3 unten).
\(SunDur_{θ,α}\) ist die vom Himmelssektor dargestellte Zeitdauer. In den meisten Sektoren ist dies gleich dem Tagesintervall (z. B. ein Monat), multipliziert mit dem Stundenintervall (z. B. eine halbe Stunde). Bei Teilsektoren (nahe dem Horizont) wird die Dauer über die sphäroidische Geometrie berechnet.
\(SunGap_{θ,α}\) ist der Lückenanteil des Sonnenkartensektors.
\(AngIn_{θ,α}\) ist der Einfallswinkel zwischen dem Schwerpunkt des Himmelssektors und der Oberflächennormalen (siehe Gleichung 4 unten).
Die relative optische Länge, \(m(θ)\), wird anhand des Zenitwinkels der Sonne und der Höhe über dem Meeresspiegel ermittelt. Bei Zenitwinkeln unter 80° kann sie über die folgende Gleichung berechnet werden:
\(m(θ) = \Large \frac {EXP (-0.000118 \space * \space Elev \space - 1.638 * 10 ^{-9} \space * \space Elev ^{2} )}{ cos(θ)} \space \space \space \normalsize (3)\)
Dabei gilt:
\(θ\) ist der Zenitwinkel der Sonne.
\(Elev\) ist die Höhe über dem Meeresspiegel in Metern.
Die Auswirkung der Oberflächenausrichtung wird durch die Multiplikation mit dem Cosinus des Einfallswinkels berücksichtigt. Der Einfallswinkel (\(AngInSky_{θ,α}\)) zwischen der abfangenden Fläche und einem gegebenen Himmelssektor mit einem Schwerpunkt im Zenit- und Azimutwinkel wird anhand der folgenden Gleichung berechnet:
\(\large AngIn_{θ,α} = acos( Cos(θ) * Cos(G_{z}) + Sin(θ) * Sin(G_{z}) * Cos(α-G_{a})) \space \space \space \normalsize (4)\)
Dabei gilt:
\(G_z\) ist der Zenitwinkel der Oberfläche. (Beachten Sie, dass bei Zenitwinkeln über 80° die Brechung eine wichtige Rolle spielt.)
\(G_a\) ist der Azimutwinkel der Oberfläche.
Berechnung der diffusen Strahlung
Für jeden Himmelssektor wird die diffuse Strahlung an dessen Schwerpunkt (\(Dif\)) berechnet, über das Zeitintervall integriert und für den Lückenanteil und den Einfallswinkel korrigiert. Dabei wird die folgende Gleichung verwendet:
\(\large Dif_{θ,α} = R_{glb} * P_{dif} * Dur * SkyGap_{θ,α} * Weight_{θ,α} * cos(AngIn_{θ,α}) \space \space \space \normalsize (5)\)
Dabei gilt:
\(\large R_{glb}\) ist die globale Normalstrahlung (siehe Gleichung 6 unten).
\(\large P_{dif}\) ist der diffuse Anteil des globalen normalen Strahlungsflusses. Dieser liegt meist bei etwa 0,2 bei sehr klarem Himmel bzw. bei 0,7 bei stark bewölktem Himmel.
\(\large Dur\) ist das Zeitintervall für die Analyse.
\(\large SkyGap_{θ,α}\) ist der Lückenanteil (Anteil des sichtbaren Himmels) im Himmelssektor.
\(\large Weight_{θ,α}\) ist das Verhältnis zwischen dem Anteil diffuser Strahlung aus einem bestimmten Himmelssektor und allen Sektoren (siehe Gleichungen 7 und 8 unten).
\(\large AngIn_{θ,α}\) ist der Einfallswinkel zwischen dem Schwerpunkt des Himmelssektors und der abfangenden Fläche.
Die globale Normalstrahlung (\(R_{glb}\)) kann durch Summierung der direkten Strahlung aus allen Sektoren (auch den verdeckten Sektoren) ohne Korrektur für den Einfallswinkel, aber mit anschließender Korrektur für den Anteil der direkten Strahlung (\(1-P_{dif}\)) berechnet werden:
\(\large R_{glb} = (S_{Const} \space Σ(β ^{m(θ)})) / (1 - P_{dif}) \space \space \space \normalsize (6)\)
Im einheitlichen diffusen Himmelsmodell wird \(Weight_{θ,α}\) wie folgt berechnet:
\(\large Weight_{θ,α} = (cosθ_{2} - cosθ_{1}) / Div_{azi} \space \space \space \normalsize (7)\)
Dabei gilt:
\(θ_1\) und \(θ_2\) sind die begrenzenden Zenitwinkel des Himmelssektors.
\(Div_{azi}\) ist die Anzahl der Azimutabschnitte auf der Himmelskarte.
Im Standardstrahlungsmodell bei Bewölkung wird \(Weight_{θ,α}\) wie folgt berechnet:
\(\large Weight_{θ,α} = (2 cosθ_{2} + cos2θ_{2} - 2 cosθ_{1} - cos 2θ_{1}) / 4 \space Div_{azi} \space \space \space \normalsize (8)\)
Diffuse Gesamtsonnenstrahlung für den Ort ($Dif_) is calculated as the sum of the diffuse solar radiation (\(Dif\)) aus allen Sektoren der Himmelskarte:
\(\large Dif_{tot} = Σ \space Dif_{θ,α} \space \space \space \normalsize (9)\)
Referenzen
Fu, P. 2000. A Geometric Solar Radiation Model with Applications in Landscape Ecology. Ph.D. Thesis, Department of Geography, University of Kansas, Lawrence, Kansas, USA.
Fu, P. und P. M. Rich. 2000. The Solar Analyst 1.0 Manual. Helios Environmental Modeling Institute (HEMI), USA.
Fu, P. und P. M. Rich. 2002. "A Geometric Solar Radiation Model with Applications in Agriculture and Forestry." Computers and Electronics in Agriculture 37:25–35.
Rich, P. M., R. Dubayah, W. A. Hetrick und S. C. Saving. 1994. "Using Viewshed Models to Calculate Intercepted Solar Radiation: Applications in Ecology. American Society for Photogrammetry and Remote Sensing Technical Papers, 524–529.
Rich, P. M. und P. Fu. 2000. "Topoclimatic Habitat Models." Proceedings of the Fourth International Conference on Integrating GIS and Environmental Modeling.