Astronomy:Solar zenith angle
The solar zenith angle is the zenith angle of the sun, i.e., the angle between the sun’s rays and the vertical direction. It is the complement to the solar altitude or solar elevation, which is the altitude angle or elevation angle between the sun’s rays and a horizontal plane.[1][2] At solar noon, the zenith angle is at a minimum and is equal to latitude minus solar declination angle. This is the basis by which ancient mariners navigated the oceans.[3]
Solar zenith angle is normally used in combination with the solar azimuth angle to determine the position of the Sun as observed from a given location on the surface of the Earth.
Formula
[math]\displaystyle{ \cos \theta_s = \sin \alpha_s = \sin \Phi \sin \delta + \cos \Phi \cos \delta \cos h }[/math]
where
- [math]\displaystyle{ \theta_s }[/math] is the solar zenith angle
- [math]\displaystyle{ \alpha_s }[/math] is the solar altitude angle, [math]\displaystyle{ \alpha_s = 90^\circ - \theta_s }[/math]
- [math]\displaystyle{ h }[/math] is the hour angle, in the local solar time.
- [math]\displaystyle{ \delta }[/math] is the current declination of the Sun
- [math]\displaystyle{ \Phi }[/math] is the local latitude.
Derivation of the formula using the subsolar point and vector analysis
While the formula can be derived by applying the cosine law to the zenith-pole-Sun spherical triangle, the spherical trigonometry is a relatively esoteric subject.
By introducing the coordinates of the subsolar point and using vector analysis, the formula can be obtained straightforward without incurring the use of spherical trigonometry.[4]
In the Earth-Centered Earth-Fixed (ECEF) geocentric Cartesian coordinate system, let [math]\displaystyle{ (\phi_{s}, \lambda_{s}) }[/math] and [math]\displaystyle{ (\phi_{o}, \lambda_{o}) }[/math] be the latitudes and longitudes, or coordinates, of the subsolar point and the observer's point, then the upward-pointing unit vectors at the two points, [math]\displaystyle{ \mathbf{S} }[/math] and [math]\displaystyle{ \mathbf{V}_{oz} }[/math], are
[math]\displaystyle{ \mathbf{S}=\cos\phi_{s}\cos\lambda_{s}{\mathbf i}+\cos\phi_{s}\sin\lambda_{s}{\mathbf j}+\sin\phi_{s}{\mathbf k}, }[/math] [math]\displaystyle{ \mathbf{V}_{oz}=\cos\phi_{o}\cos\lambda_{o}{\mathbf i}+\cos\phi_{o}\sin\lambda_{o}{\mathbf j}+\sin\phi_{o}{\mathbf k}. }[/math]
where [math]\displaystyle{ {\mathbf i} }[/math], [math]\displaystyle{ {\mathbf j} }[/math] and [math]\displaystyle{ {\mathbf k} }[/math] are the basis vectors in the ECEF coordinate system.
Now the cosine of the solar zenith angle, [math]\displaystyle{ \theta_{s} }[/math], is simply the dot product of the above two vectors
[math]\displaystyle{ \cos\theta_{s} = \mathbf{S}\cdot\mathbf{V}_{oz} = \sin\phi_{o}\sin\phi_{s} + \cos\phi_{o}\cos\phi_{s}\cos(\lambda_{s}-\lambda_{o}). }[/math]
Note that [math]\displaystyle{ \phi_{s} }[/math] is the same as [math]\displaystyle{ \delta }[/math], the declination of the Sun, and [math]\displaystyle{ \lambda_{s}-\lambda_{o} }[/math] is equivalent to [math]\displaystyle{ -h }[/math], where [math]\displaystyle{ h }[/math] is the hour angle defined earlier. So the above format is mathematically identical to the one given earlier.
Additionally, Ref. [4] also derived the formula for solar azimuth angle in a similar fashion without using spherical trigonometry.
Minimum and Maximum
At any given location on any given day, the solar zenith angle, [math]\displaystyle{ \theta_{s} }[/math], reaches its minimum, [math]\displaystyle{ \theta_\text{min} }[/math], at local solar noon when the hour angle [math]\displaystyle{ h = 0 }[/math], or [math]\displaystyle{ \lambda_{s}-\lambda_{o}=0 }[/math], namely, [math]\displaystyle{ \cos\theta_\text{min} = \cos(|\phi_{o}-\phi_{s}|) }[/math], or [math]\displaystyle{ \theta_\text{min} = |\phi_{o}-\phi_{s}| }[/math]. If [math]\displaystyle{ \theta_\text{min} \gt 90^{\circ} }[/math], it is polar night.
And at any given location on any given day, the solar zenith angle, [math]\displaystyle{ \theta_{s} }[/math], reaches its maximum, [math]\displaystyle{ \theta_\text{max} }[/math], at local midnight when the hour angle [math]\displaystyle{ h = -180^{\circ} }[/math], or [math]\displaystyle{ \lambda_{s}-\lambda_{o}=-180^{\circ} }[/math], namely, [math]\displaystyle{ \cos\theta_\text{max} = \cos(180^{\circ}-|\phi_{o}+\phi_{s}|) }[/math], or [math]\displaystyle{ \theta_\text{max} = 180^{\circ}-|\phi_{o}+\phi_{s}| }[/math]. If [math]\displaystyle{ \theta_\text{max} \lt 90^{\circ} }[/math], it is polar day.
Caveats
The calculated values are approximations due to the distinction between common/geodetic latitude and geocentric latitude. However, the two values differ by less than 12 minutes of arc, which is less than the apparent angular radius of the sun.
The formula also neglects the effect of atmospheric refraction.[5]
Applications
Sunrise/Sunset
Sunset and sunrise occur (approximately) when the zenith angle is 90°, where the hour angle h0 satisfies[2] [math]\displaystyle{ \cos h_0 = -\tan \Phi \tan \delta. }[/math]
Precise times of sunset and sunrise occur when the upper limb of the Sun appears, as refracted by the atmosphere, to be on the horizon.
Albedo
A weighted daily average zenith angle, used in computing the local albedo of the Earth, is given by [math]\displaystyle{ \overline{\cos \theta_s} = \frac{\displaystyle \int_{-h_0}^{h_0} Q \cos \theta_s \, \text{d}h}{\displaystyle \int_{-h_0}^{h_0} Q \, \text{d}h} }[/math] where Q is the instantaneous irradiance.[2]
Summary of special angles
For example, the solar elevation angle is :
- 90° if you are on the equator, a day of equinox, at a solar hour of twelve
- near 0° at the sunset or at the sunrise
- between −90° and 0° during the night (midnight)
An exact calculation is given in position of the Sun. Other approximations exist elsewhere.[6]
See also
- Azimuth
- Solar azimuth angle
- Horizontal coordinate system
- List of orbits
- Photovoltaic mounting system § Orientation and inclination
- Position of the Sun
- Sun path
- Sunrise
- Sunset
- Sun transit time
References
- ↑ Jacobson, Mark Z. (2005). Fundamentals of Atmospheric Modeling (2nd ed.). Cambridge University Press. p. 317. ISBN 0521548659. https://archive.org/details/fundamentalsatmo00jaco.
- ↑ 2.0 2.1 2.2 Hartmann, Dennis L. (1994). Global Physical Climatology. Academic Press. p. 30. ISBN 0080571638. https://archive.org/details/globalphysicalcl00hart.
- ↑ Bonan, Gordon (2005). Ecological climatology: concepts and applications. Cambridge University Press. p. 62. ISBN 9781316425190. https://www.cambridge.org/core/books/ecological-climatology/D146443B007985BC366B2512345692C0. Retrieved 13 November 2019.
- ↑ 4.0 4.1 Zhang, T., Stackhouse, P.W., Macpherson, B., and Mikovitz, J.C., 2021. A solar azimuth formula that renders circumstantial treatment unnecessary without compromising mathematical rigor: Mathematical setup, application and extension of a formula based on the subsolar point and atan2 function. Renewable Energy, 172, 1333-1340. DOI: https://doi.org/10.1016/j.renene.2021.03.047
- ↑ Woolf, Harold M. (1968). "On the computation of solar elevation angles and the determination of sunrise and sunset times". NASA Technical Memorandu, X-1646 (Washington, D.C.): 3.
- ↑ livioflores-ga. "Equation to know where the Sun is at a given place at a given date-time". http://answers.google.com/answers/threadview/id/782886.html. Retrieved 9 March 2013.
 |  |
