Model help:TopoFlow-Meteorology: Difference between revisions
No edit summary |
No edit summary |
||
Line 520: | Line 520: | ||
| surface temperature | | surface temperature | ||
| deg C | | deg C | ||
|- | |- | ||
| p<sub>0</sub> | | p<sub>0</sub> | ||
Line 680: | Line 676: | ||
| number of hours before (-) or after (+) true solar noon (TSN) | | number of hours before (-) or after (+) true solar noon (TSN) | ||
| radians | | radians | ||
|- | |||
| T | |||
| temperature, measured | |||
| degrees Celsius | |||
|- | |||
| RH | |||
| temperature, measured | |||
| degrees Celsius | |||
|- | |||
| γ<sub>dust</sub> | |||
| dust attenuation, equals to 0.08 (this is the default value) | |||
| - | |||
|- | |||
| α | |||
| aspect angle, measured clockwise from north | |||
| radians | |||
|- | |||
| β | |||
| slope angle, equal to atan(slope) | |||
| radians | |||
|- | |||
| DPY | |||
| earth days per year, equals to 365.2425 | |||
| days | |||
|- | |||
| DPY<sub>sidereal</sub> | |||
| earth days per year (sidereal), equals to 366.2425 | |||
| days | |||
|- | |||
| angle | |||
| earth tilt angle, equals to 23.4397 | |||
| degrees | |||
|- | |||
| e (in shortwave radiation equations) | |||
| earth orbit eccentricity, equal to 0.016713 | |||
| - | |||
|- | |||
|} | |||
'''Output''' | |||
{| {{Prettytable}} class="wikitable sortable" | |||
!Symbol!!Description!!Unit | |||
|- | |||
| e<sub>air</sub> | |||
| vapor pressure of air | |||
| mbar | |||
|- | |||
| e<sub>surf</sub> | |||
| vapor pressure at the surface | |||
| mbar | |||
|- | |||
| Qn<sub>SW</sub> | |||
| net shortwave radiation | |||
| W / m^2 | |||
|- | |||
| Q<sub>LW</sub> | |||
| net longwave radiation, equals to LW<sub>in</sub> - LW<sub>out</sub> | |||
| W / m^2 | |||
|- | |||
| em<sub>air</sub> | |||
| air emissivity in [0,1] | |||
| - | |||
|- | |||
| LW<sub>in</sub> | |||
| incoming longwave radiation | |||
| - | |||
|- | |||
| LW<sub>out</sub> | |||
| outgoing longwave radiation | |||
| - | |||
|- | |||
| e<sub>sat</sub> | |||
| saturation vapor pressure | |||
| mbar | |||
|- | |- | ||
| T<sub>hr</sub> | | T<sub>hr</sub> | ||
Line 692: | Line 762: | ||
| instantaneous extraterrestrial radiation flux on a horizontal plane | | instantaneous extraterrestrial radiation flux on a horizontal plane | ||
| watts / m<sup>2</sup> | | watts / m<sup>2</sup> | ||
|- | |- | ||
| e<sub>sat</sub> | | e<sub>sat</sub> | ||
Line 723: | Line 785: | ||
| M<sub>opt</sub> | | M<sub>opt</sub> | ||
| optical air mass (from Kasten and Young, 1989) | | optical air mass (from Kasten and Young, 1989) | ||
| - | | - | ||
|- | |- | ||
Line 752: | Line 810: | ||
| backscattered radiation flux | | backscattered radiation flux | ||
| watts/m<sup>2</sup> | | watts/m<sup>2</sup> | ||
|- | |- | ||
| ΔΩ | | ΔΩ | ||
Line 792: | Line 842: | ||
| instantaneous clear sky radiation | | instantaneous clear sky radiation | ||
| watts/m<sup>2</sup> | | watts/m<sup>2</sup> | ||
|- | |- | ||
| Ω | | Ω | ||
| earth rotation rate | | earth rotation rate | ||
| radians/year | | radians/year | ||
|- | |- | ||
|} | |} |
Revision as of 17:51, 31 August 2011
TopoFlow-Meteorology
The module is the meteorology process component for a D8-based, spatial hydrologic model
Model introduction
This component reads a variety of variables for the atmosphere and for the land surface from input files or as simple scalars. It then computes many additional variables, such as vapor pressure, eair, and net shortwave (solar) radiation, QnSW, using built-in shortwave and longwave radiation calculators that are based on celestial mechanics and widely-used empirical relationships. These additional variables are needed by the Snowmelt → Energy Balance and Evaporation → Energy Balance components. Direct, diffuse and back-scattered radiation fluxes are all modeled. Properties of the atmosphere (e.g. precipitation rate, P, air temperature, Tair, relative humidity, RH, and dust attenuation, γ), are used as well as surface/topographic properties (e.g. slope angle, aspect angle and surface albedo. The approach used here closely follows the one outlined in Appendix E of Dingman (2002). However, instantaneous vs. day-integrated radiation fluxes are used and the optical air mass is modeled using the widely used method of Kasten and Young.
Model parameters
Uses ports
• Snow (Snowmelt)
Provides ports
• Meteorology
• Configure (tabbed dialog GUI to change settings)
• Run (only if used as the Driver)
Main equations
Equations Used to Compute Longwave Radiation
- Net longwave radiation
[math]\displaystyle{ Q_{LW} = \left (LW_{in} - LW_{out}\right ) }[/math] (1)
- Incoming longwave radiation (using Stefan-Bolzman law)
[math]\displaystyle{ LW_{in} = em_{air} * \sigma * \left (T_{air} + 273.15 \right )^4 }[/math] (2)
- Outgoing longwave radiation (using Stefan-Bolzman law)
[math]\displaystyle{ LW_{out} = em_{surf} * \sigma * \left (T_{surf} + 273.15 \right )^4 + \left ( 1 - em_{surf} \right ) * LW_{in} }[/math] (3)
- Saturation vapor pressure
[math]\displaystyle{ e_{sat} = 6.11 * exp [\left ( 17.3 * T \right ) / \left ( T + 237.3 \right )] }[/math] (4)
- Vapor pressure of air
[math]\displaystyle{ e_{air} = e_{sat} \left (T_{air}\right ) * RH }[/math] (5)
- Emissivity of air
[math]\displaystyle{ em_{air} = [\left (1 - F \right ) * 1.72 * [e_{air}/ \left (T_{air} + 273.15 \right )]^ \left ({\frac{1}{7}} \right ) * \left (1 + 0.22 * C^2 \right )] + F }[/math] (6)
Equations Used to Compute Shortwave Radiation
- Day angle
[math]\displaystyle{ \Gamma = 2 * \pi * \left ( J - 1 \right ) / 365 }[/math] (1)
- Earth orbit eccentricity correlation
[math]\displaystyle{ E_{0} = \left ( r_{0} / r \right )^2 = 1.000110 + 0.034221 * cos \left ( \Gamma \right ) + 0.001280 * sin \left ( \Gamma \right ) + 0.000719 * cos \left ( 2 * \Gamma \right ) + 0.000077 * sin \left ( 2 * \Gamma \right ) }[/math] (2)
- Declination of the sun at a given latitude
[math]\displaystyle{ \delta = \left ( 180 / \pi \right ) * [0.399912 * cos \left ( \Gamma \right ) + 0.070257 * sin \left ( \Gamma \right ) - 0.006758 * cos \left ( 2 * \Gamma \right ) + 0.000907 * sin \left ( 2 * \Gamma \right ) - 0.002697 * cos \left ( 3 * \Gamma \right ) + 0.00148 * sin \left ( 3 * \Gamma \right )] }[/math] (3)
- Zenith angle
[math]\displaystyle{ \theta = cos ^ \left ( -1 \right ) [sin \left ( \Lambda \right ) * sin \left ( \delta \right ) + cos \left ( \Lambda \right ) + cos \left ( \Lambda \right ) * cos \left (\delta \right ) * cos \left ( \omega * t \right )] }[/math] (4)
- Sunrise offset (hours before true solar noon)
[math]\displaystyle{ T_{hr} = - {\frac{cos^ \left ( -1 \right ) [- tan \left ( \delta \right ) * tan \left ( \Lambda \right )]}{\omega}} }[/math] (5)
- Sunset offset (hours after true solar noon)
[math]\displaystyle{ T_{hs} = + {\frac{cos^ \left ( -1 \right ) [- tan \left ( \delta \right ) * tan \left ( \Lambda \right )]}{\omega}} }[/math] (6)
- Instantaneous extraterrestrial radiation flux on a horizontal plane
[math]\displaystyle{ k^' = I_{sc} * E_{0} * [cos \left ( \delta \right ) * cos \left ( \Lambda \right ) * cos \left ( \omega * t \right ) + sin \left ( \delta \right ) * sin \left ( \Lambda \right )] }[/math] (7)
- Dew point
[math]\displaystyle{ T_{dew} = {\frac{ln \left ( e\right ) + 0.4926}{0.0708 - 0.00421 * ln \left ( e \right )}} }[/math] (8)
- Precipitable water content
[math]\displaystyle{ W_{p} = 1.12 * exp \left ( 0.0614 * T_{dew} \right ) }[/math] (9)
- Optical air mass (from Kasten and Young, 1989)
[math]\displaystyle{ M_{opt} = 1 / [cos \left ( \theta \right ) + 0.50572 * \left ( 96.07995 - \theta \right ) ^ \left (-1.6364 \right )] }[/math] (10)
- Total atmospheric transmissivity
[math]\displaystyle{ \tau = exp \left ( -0.124 - 0.0207 * W_{p} - \left (0.0682 + 0.0248 * W_{p} \right ) * M_{opt} \right ) }[/math] (11)
- Instantaneous direct radiation flux
[math]\displaystyle{ k_{direct}^ ' = \tau \, k_{ET}^ ' }[/math] (12)
- Scattering attenuation
[math]\displaystyle{ \gamma_{s} = 1 - exp \left ( -0.0363 - 0.0084 * W_{p} - \left (0.0572 + 0.0173 * W_{p} \right ) * M_{opt} \right ) + \gamma_{dust} }[/math] (13)
- Instantaneous diffuse radiation flux
[math]\displaystyle{ k_{diffuse}^ ' = 0.5 * \gamma_{s} * K_{ET} ^ ' }[/math] (13)
- Instantaneous global radiation flux
[math]\displaystyle{ k_{global}^' = k_{direct}^' + k_{diffuse}^' }[/math] (14)
- Backscattered radiation flux
[math]\displaystyle{ k_{bs}^' = 0.5 * \gamma_{s} * a * K_{g}^' }[/math] (15)
- Longitude offset
[math]\displaystyle{ \Delta \Omega = tan^ \left (-1\right ) [{\frac{sin \left ( \beta \right ) * sin \left ( \alpha \right )}{cos \left ( \beta \right ) * cos \left ( \Lambda \right ) - sin \left ( \beta \right ) * sin \left ( \Lambda \right ) * cos \left ( \alpha \right )}}] }[/math] (16)
- Equivalent latitude
[math]\displaystyle{ \Lambda_{eq} = sin^ \left (-1\right ) [sin \left (\beta \right ) * cos \left ( \alpha \right ) * cos \left ( \Lambda \right ) + cos \left ( \beta \right ) * sin \left ( \Lambda \right )] }[/math] (17)
- Sunrise offset for sloped topography
[math]\displaystyle{ T_{sr} = - {\frac{cos^ \left (-1 \right ) [- tan \left ( \Lambda_{eq} \right ) * tan \left ( \delta \right ) ] - \Delta \Omega}{\omega}} }[/math] (18)
- Sunset offset for sloped topography
[math]\displaystyle{ T_{ss} = + {\frac{cos^ \left (-1 \right ) [- tan \left ( \Lambda_{eq} \right ) * tan \left ( \delta \right ) ] - \Delta \Omega}{\omega}} }[/math] (19)
- Instantaneous extra-terrestrial radiation flux for sloped topography
[math]\displaystyle{ k_{ET} ^' = I_{sc} E_{\theta} [cos \left ( \delta \right ) cos \left ( \Lambda_{eg} \right ) cos \left ( \omega t_{h} + \Delta \Omega \right ) + sin \left ( \delta \right ) sin \left ( \Lambda_{eg} \right ) ] }[/math] (20)
- Instantaneous clear sky radiation
[math]\displaystyle{ k_{cs} = \tau * k_{ET}^' + k_{diffuse}^' + k_{bs} ^' }[/math] (21)
- Earth rotation rate
[math]\displaystyle{ \Omega = 2 * \pi * DPY_{sidereal} }[/math] (22)
Symbol | Description | Unit |
---|---|---|
ρH2O | density of water | kg / m^3 |
Cpair | heat capacity of air | J / kg / K |
ρair | density of air | kg / m^3 |
precip. rate | precipitation rate | mm / hr |
Tair | air temperature | deg C |
Tsurf | surface temperature | deg C |
p0 | atmospheric pressure | mbar |
z | reference height for uz | m |
uz | wind velocity at reference height | m / s |
z0 | surface roughness length scale for wind | m |
albedo | surface albedo in [0,1] | - |
emair | surface emissivity in [0,1] | - |
dust atten. in [0, 0.3] | dust attenuation factor | - |
cloud factor | cloud factor in [0,1], 0 for no clouds | - |
canopy factor | forest canopy factor in [0,1], 0 for no canopy | - |
slope | topographic slope in [0, infinity] | - |
slope grid file | as flat binary, row-major file with 4-byte floats | - |
aspect angle | aspect angle [radians] in [0,1] | - |
aspect grid file | as flat binary, row-major file with 4-byte floats | - |
time zone offset | offset, in hours, from Greenwich Mean Time (GMT), negative for east of prime meridian, positive otherwise | - |
start month | start month for solar radiation calculations | - |
start day | start day for solar radiation calculations | - |
start hour | start hour for solar radiation calculations [decimal, 24-hour clock] | - |
h | height above mean sea level | m |
ρ | equals to ρ(h), air density at height h | kg / m3 |
ρ0 | air density at h = 0 | kg / m3 |
δ0 | equals to n0 - 1 | - |
n0 | refractive index for air at 0.7μm wavelength at h = 0 | - |
R | mean earth radius | m |
m(γ) | the relative optical air mass at solar elevation γ | - |
mabs(γ) | absolute optical air mass | - |
emsurf(γ) | emissivity of the surface (e.g. snow) | - |
σ | Stefan-Boltzman constant, equals to 5.67e-8 | W/(m2 K4 |
C | the fraction of sky covered by clouds | - |
F | the fraction of sky covered by forest canopy | - |
sc | solar constant, 1 sc equals to 1367 | watts/m2 |
J | Julian day, decimal days (J =1 on January 1, and J = 365 on December 31) | - |
Γ | day angle or fractional year, radians | - |
E0 | earth orbit eccentricity correlation | - |
Λ | latitude (for a given point on the earth's surface), in radians | - |
Λdegree | latitude (for a given point on the earth's surface), in degree | - |
δ | declination of the sun | degree |
ω | earth angular velocity, equals to 0.2618 (15 degree per hour) | radians/hour |
θ | zenith angle | radians |
t | number of hours before (-) or after (+) true solar noon (TSN) | radians |
T | temperature, measured | degrees Celsius |
RH | temperature, measured | degrees Celsius |
γdust | dust attenuation, equals to 0.08 (this is the default value) | - |
α | aspect angle, measured clockwise from north | radians |
β | slope angle, equal to atan(slope) | radians |
DPY | earth days per year, equals to 365.2425 | days |
DPYsidereal | earth days per year (sidereal), equals to 366.2425 | days |
angle | earth tilt angle, equals to 23.4397 | degrees |
e (in shortwave radiation equations) | earth orbit eccentricity, equal to 0.016713 | - |
Output
Symbol | Description | Unit |
---|---|---|
eair | vapor pressure of air | mbar |
esurf | vapor pressure at the surface | mbar |
QnSW | net shortwave radiation | W / m^2 |
QLW | net longwave radiation, equals to LWin - LWout | W / m^2 |
emair | air emissivity in [0,1] | - |
LWin | incoming longwave radiation | - |
LWout | outgoing longwave radiation | - |
esat | saturation vapor pressure | mbar |
Thr | sunrise offset, hours before true solar noon (TSN) | - |
Ths | sunset offset, hours after true solar noon (TSN) | - |
k' ET | instantaneous extraterrestrial radiation flux on a horizontal plane | watts / m2 |
esat | saturated vapor pressure | kPa |
e | vapor pressure | kpa |
Tdew | dew point | degree Celsius |
Wp | precipitation water content | cm |
Mopt | optical air mass (from Kasten and Young, 1989) | - |
Mopt | optical air mass (from Kasten and Young, 1989) | - |
τ | total atmospheric transmissivity, in [0,1] | - |
kdirect ' | instantaneous direct radiation flux | watts / m2 |
γs | scattering attenuation | - |
kdiffuse' | instantaneous diffuse radiation flux | watts/m2 |
kglobal' | instantaneous global radiation flux | watts/m2 |
kbs' | backscattered radiation flux | watts/m2 |
ΔΩ | longitude offset | radians |
Λeq | equivalent latitude | radians |
noon_offset | noon offset for sloped topography, equals to - ΔΩ / ω | hours |
Tsr | sunrise offset for sloped topography | hours |
Tss | sunset offset for sloped topography | hours |
day_length | day length for sloped topography, equals to Tss - Tsr | hours |
kET' | instantaneous extra-radiation flux for sloped topography | watts/m2 |
kcs | instantaneous clear sky radiation | watts/m2 |
Ω | earth rotation rate | radians/year |
Notes
Notes on Input Parameters
For each input variable, you may choose from the droplist of data types. For the "Scalar" data type, enter a numeric value with the units indicated in the dialog. For the other data types, enter a filename. Values in files must also use the indicated units.
Single grids and grid sequences are assumed to be stored as RTG and RTS files, respectively. Time series are assumed to be stored as text files, with one value per line. For a time series or grid sequence, the time between values must coincide with the timestep provided.
For DEMs with pixel geometry and bounding box given in terms of Geographic coordinates, the latitude and longitude of each pixel is used in the calculations. For DEMs with a "fixed-length" pixel geometry (e.g. UTM coordinates), which tend to span smaller areas, the dialog prompts for a single lat/lon pair to be used in the calculations.
In the calculation, the timestep between frames in the new grid sequence (RTS file) should typically be about one hour and should match the timestep that will be used to model the Snowmelt or Evaporation processes. The number of frames in the RTS file will depend on the start and stop times, as well as the timestep. The start time, stop time and timestep should match those used to create the new shortwave radiation file with extension "*.Qn-SW".
Notes on the Equations
All variables and their units can be seen by expanding the Nomenclature section above.
time zone: Boundaries of time zones can be very irregular and a time zone map should be consulted if you are unsure. The time zone is not simply a function of the longitude. You can select an adjacent time zone to include the effect of Daylight Savings Time. Time zones with non-integer offsets from GMT are not yet supported.
slope: Topographic slopes (not slope angles) are specified as dimensionless numbers [m/m]. A RiverTools grid (RTG file) with extension "_slope.rtg", "_mf-slope.rtg" or "_dinf-slope.rtg" can be used.
aspect: Aspect is specified as an angle measured in radians counter-clockwise from due east (the standard convention). A RiverTools grid (RTG file) with extension "_mf-angle.rtg" or "_dinf-angle.rtg" can be used for the (continuous-angle) aspect grid.
QSW is set to zero between the times of local sunset and local sunrise, so frames in the RTS file that correspond to nighttime hours will contain only zeros.
Notes on the Equations to Compute Longwave Radiation
The value 237.3 in Brutsaert's equation for emair is not a misprint.
Notes on the Equations to Compute Shortwave Radiation
The declination of the sun reaches its lowest value of -23.5 degrees on the Winter Solstice (Dec. 21/22) and reaches its highest value of 23.5 degrees on the Summer Solstice(June 21/22). It is zero for both the Vernal Equinox (Mar. 20/21) and the Autumnal Equinox (Sept. 22/23). The value of 23.4397 degrees is the fixed tilt angle of the Earth's axis from from the plane of the ecliptic.
The equation for declination of the sun is a Fourier series expansion that is supposed to have a maximum error of 0.0006 radians (less than 3 arcminutes). This should be double-checked. The equation shown here can also be found at: http://en.wikipedia.org/wiki/Declination, and is derived in the paper called "Fourier series representation of the position of the sun", by J.W. Spencer (1971), of CSIRO.
If (abs(latdeg) gt 66.5) we are above Arctic circle or below Antarctic circle. This can also happen if latdeg is an "equivalent latitude" for sloped topography. In this case, the absolute value of the argument to ACOS can exceed one and there is either no sunrise or no sunset at the given location. If the argument to ACOS is less than -1, it is set to -1. If it is greater than 1 it is set to 1.
True solar noon (TSN) for a given location is the time when the sun reaches its highest point above the horizon. At this time, the zenith angle attains its minimum value for the day and the sun is said to be "on the meridian." In order to compute the clock time at which true solar noon occurs, we start with the number 12, then add the following three corrections:
(1) ΔTc is the time difference between true solar noon and local clock noon, without accounting for any arbitrary time zone adjustments. It is computed from the so-called Equation of Time and is related to the figure-8-shaped "analemma".
(2)ΔTtz is computed as the (signed) difference between the local longitude and the longitude at the center of the local time zone, divided by 15. The Earth turns through 15 degrees of longitude in one hour. The longitude at the center of a time zone is given by its GMT offset (and integer between -12 and 12) times 15.
(3) ΔTdst is a correction for Daylight Savings Time, generally an integer. (Reference: Whitman, A.M. (2003) A Simple Expression for the Equation of Time (online document: http://www.sunspot.noao.edu/sunspot/pr/ answerbook/expl-5.html)).
Sunrise and sunset occur when zenith angle, θ, is equal toπ/2, so cos(θ) = 0. Using the equation for the zenith angle we can then solve for time offsets (e.g. sunrise and sunset offsets) as th, in hours.
Many of these equations are from Dingman, Appendix E. In his notation, a lower-case "k" indicates an instantaneous radiation flux. In expressions for "k", the argument "th" is the number of hours before (th < 0) or after (th > 0) true solar noon. Upper-case "K" then indicates the integral of "k" when th is allowed to vary from the sunrise offset to the sunset offset. Primes on "k" or "K" indicate values for a horizontal plane which must be corrected (as explained by Dingman) for sloped topography. For example, the sunrise and sunset times for a given location depend on the slope and aspect of the local topography.
The equation for kET' is similar to Dingman, App. E, E-25, but his was integrated over one day and this one is instantaneous. Need to double check the integrating the equation here over one day gives E-25. See solar_funcs.py, which is imported by met_base.py. Note that ω has units of radians/hour and ΔΩ (longitude offset) has units of radians.
TopoFlow provides two methods for computing the saturated vapor pressure, esat = e*, as a function of temperature, T. One is from Brutsaert (1975) (and used in Dingman) and the other is from Satterlund (1979). When plotted, the two curves look very similar.
By definition, the relative humidity is the ratio of the actual vapor pressure to the saturated vapor pressure. That is, RH = ea / ea *. However, relative humidity is generally measured and used together with an equation for saturated vapor pressure to compute the vapor pressure as a function of temperature.
The equation for optical air mass is from Kasten and Young (1989) and depends on latitude, declination and time of day. Instead of this equation, Dingman plots a family of curves in his Figure E-4 (p. 605) for *average daily* optical air mass. It would be interesting to integrate the equation given for Mopt over one day and compare the two by creating a contour plot.
Typical clear-sky values for γdust are between 0 and 0.2. See Dingman, p. 604-605.
A table of typical surface albedos is given by Dingman, Table D-2 on page 584.
The Earth orbit eccentricy *correction*, E0, depends on the day angle, Γ, while the Earth orbit eccentricity is computed as: e = (b-a)/a, where a and b are the semi-major and semi-minor axes of the elliptical orbit.
Examples
An example run with input parameters, BLD files, as well as a figure / movie of the output
Follow the next steps to include images / movies of simulations:
- Upload file: http://csdms.colorado.edu/wiki/Special:Upload
- Create link to the file on your page: [[Image:<file name>]].
See also: Help:Images or Help:Movies
Developer(s)
References
Dingman, S.L (2002) Physical Hydrology, 2nd ed., Prentice Hall, New Jersey. (see Appendix E)
Brutsaert, W. (1975) On a derivable formula for long-wave radiation from clear skies, Water Resources Research, 11, 742-744.
Kasten and Young (1989) Revised optical air mass tables and approximation formula. Applied Optics, 28 (22): 4735~4738. (for the optical air mass equation)
Liston, G. *******
Marks and Dozier (1992) Climate and Energy exchange at the Snow Surface in the Alpine Region of the Sierra Nevada 1. Meteorological Measurements and Monitoring. Water Resources Research, 28(11), 3029~3042.
Whitman, A.M. (2003) A simple expression for the equation of time, online document, http://www.sunspot.noao.edu/sunspot/pr/answerbook/expl-5.html