Model help:TopoFlow-Evaporation-Energy Balance: Difference between revisions
No edit summary |
m (Text replacement - "http://csdms.colorado.edu/wiki/" to "https://csdms.colorado.edu/wiki/") |
||
(20 intermediate revisions by 3 users not shown) | |||
Line 3: | Line 3: | ||
1) Log in to the wiki | 1) Log in to the wiki | ||
2) Create a new page for each model, by using the following URL: | 2) Create a new page for each model, by using the following URL: | ||
* | * https://csdms.colorado.edu/wiki/Model help:<modelname> | ||
* Replace <modelname> with the name of a model | * Replace <modelname> with the name of a model | ||
3) Than follow the link "edit this page" | 3) Than follow the link "edit this page" | ||
Line 28: | Line 28: | ||
|width="20%"| Component status | |width="20%"| Component status | ||
|width="60%"| Enabled / Disabled | |width="60%"| Enabled / Disabled | ||
|width="20%"| - | |width="20%"| [-] | ||
|- | |- | ||
| Input directory | | Input directory | ||
Line 121: | Line 121: | ||
==Uses ports== | ==Uses ports== | ||
< | |||
• Meteorology <br /> | |||
• Channels (surface water flow in a network of channels with trapezoidal cross-section) <br /> | |||
• Snow (Snowmelt) <br /> | |||
• Infil (Infiltration) <br /> | |||
• Satzone (Subsurface flow in saturated zone) <br /> | |||
==Provides ports== | ==Provides ports== | ||
< | |||
• Evap (Evaporation) <br /> | |||
• Configure (tabbed dialog GUI to change settings) <br /> | |||
• Run (only if used as the Driver) <br /> | |||
==Main equations== | ==Main equations== | ||
* Evaporation rate | * Evaporation rate | ||
::::{| | ::::{| | ||
|width=500px|<math>ET=\left (1000 Q_{et}\right) / \left ( | |width=500px|<math>ET=\left (1000 \ast Q_{et}\right) / \left (\rho_{water} \ast L_{v}\right) </math> | ||
|width=50px align="right"|(1) | |width=50px align="right"|(1) | ||
|} | |} | ||
Line 139: | Line 147: | ||
* Conduction energy flux | * Conduction energy flux | ||
::::{| | ::::{| | ||
|width=500px|<math>Q_{c}= K_{soil} \left (T_{soil_x} - T_{surf} \right) \left ( 100 / x \right) </math> | |width=500px|<math>Q_{c}= K_{soil} \ast \left (T_{soil_x} - T_{surf} \right) \ast \left ( 100 / x \right) </math> | ||
|width=50px align="right"|(3) | |width=50px align="right"|(3) | ||
|} | |} | ||
* Sensible heat flux | * Sensible heat flux | ||
::::{| | ::::{| | ||
|width=500px|<math>Q_{h}= | |width=500px|<math>Q_{h}= \rho_{air} \ast c_{air} \ast D_{h} \ast \left (T_{air} - T_{surf}\right) </math> | ||
|width=50px align="right"|(4) | |width=50px align="right"|(4) | ||
|} | |} | ||
* Bulk exchange coeff. (neutrally stable conditions) | * Bulk exchange coeff. (neutrally stable conditions) | ||
::::{| | ::::{| | ||
|width=500px|<math>D_{n}= u_{z} | |width=500px|<math>D_{n}= u_{z} \ast \kappa^2 / LN [ \left ( z - h_{snow}\right) / z0_{air}] ^2 </math> | ||
|width=50px align="right"|(5) | |width=50px align="right"|(5) | ||
|} | |} | ||
* Bulk exchange coeff. for heat | * Bulk exchange coeff. for heat | ||
::::{| | ::::{| | ||
|width=500px|<math>D_{h}= D_{n} / \left (1 + \left (10 Ri \right)\right) </math> | |width=500px|<math>D_{h}=\left\{\begin{matrix} D_{n} / \left (1 + \left (10 \ast Ri \right)\right) & stable: T_{air} > T_{surf} \\ D_{n} \ast [ 1 - \left ( 10 \ast Ri \right) ] & unstable: T_{air} < T_{surf} \end{matrix}\right. </math> | ||
|width=50px align="right"|(6) | |width=50px align="right"|(6) | ||
|} | |} | ||
* | * Richardson's number | ||
::::{| | ::::{| | ||
|width=500px|<math> | |width=500px|<math>Ri= g \ast z \ast \left (T_{air} - T_{surf} \right) / [ u_{z}^2 \left ( T_{air} + 273.15 \right) ] </math> | ||
|width=50px align="right"|(7) | |width=50px align="right"|(7) | ||
|} | |} | ||
Line 218: | Line 221: | ||
| m | | m | ||
|- | |- | ||
| | | ρ<sub>air</sub> | ||
| density of the air | | density of the air | ||
| kg / m<sup>3</sup> | | kg / m<sup>3</sup> | ||
Line 226: | Line 229: | ||
| J /kg | | J /kg | ||
|- | |- | ||
| L<sub>v</sub> | | L<sub>v</sub> | ||
| latent heat of vaporization, water (2500000) | | latent heat of vaporization, water (2500000) | ||
Line 266: | Line 268: | ||
| Rechardson's number | | Rechardson's number | ||
| - | | - | ||
|- | |||
| ρ<sub>water</sub> | |||
| density of the water | |||
| kg / m<sup>3</sup> | |||
|- | |- | ||
|} | |} | ||
Line 271: | Line 277: | ||
</div> | </div> | ||
==Notes== | ==Notes== | ||
'''''Notes on Input Parameters''''' | |||
If net total radiation has been measured, it can be entered as Q<sub>SW</sub> and then Q<sub>LW</sub> can be set to zero. Any meteorological variables entered here (such as T<sub>air</sub>) are automatically shared with other other processes, such as Snowmelt and Precipitation. | If net total radiation has been measured, it can be entered as Q<sub>SW</sub> and then Q<sub>LW</sub> can be set to zero. Any meteorological variables entered here (such as T<sub>air</sub>) are automatically shared with other other processes, such as Snowmelt and Precipitation. | ||
Line 280: | Line 288: | ||
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. | 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. | ||
'''''Notes on the Equations''''' | |||
All variables and their units can be seen by expanding the Nomenclature section above. | |||
Wherever (d > 0), evaporation results in a reduction in the surface flow depth. Wherever (d = 0), water is removed from subsurface storage. If the 1D Richards' equation is used for infiltration, then the evaporation rate is applied as a surface boundary condition and alters | Wherever (d > 0), evaporation results in a reduction in the surface flow depth. Wherever (d = 0), water is removed from subsurface storage. If the 1D Richards' equation is used for infiltration, then the evaporation rate is applied as a surface boundary condition and alters | ||
In the equation for computing D<sub>n</sub>, the reference height, z, is reduced by the computed snow depth, h<sub>snow</sub>. It follows that z must be chosen so as to be greater than any possible snow depth. | In the equation for computing D<sub>n</sub>, the reference height, z, is reduced by the computed snow depth, h<sub>snow</sub>. It follows that z must be chosen so as to be greater than any possible snow depth. | ||
==Examples== | ==Examples== | ||
Line 289: | Line 300: | ||
<span class="remove_this_tag">Follow the next steps to include images / movies of simulations:</span> | <span class="remove_this_tag">Follow the next steps to include images / movies of simulations:</span> | ||
* <span class="remove_this_tag">Upload file: | * <span class="remove_this_tag">Upload file: https://csdms.colorado.edu/wiki/Special:Upload</span> | ||
* <span class="remove_this_tag">Create link to the file on your page: <nowiki>[[Image:<file name>]]</nowiki>.</span> | * <span class="remove_this_tag">Create link to the file on your page: <nowiki>[[Image:<file name>]]</nowiki>.</span> | ||
Line 295: | Line 306: | ||
==Developer(s)== | ==Developer(s)== | ||
[[User: | [[User:Peckhams|Scott Peckham]] | ||
==References== | ==References== | ||
Line 307: | Line 318: | ||
==Links== | ==Links== | ||
* [[ | |||
* [[ | '''''Related Help Pages''''' | ||
* [[Model help:TopoFlow-Evaporation-Priestley_Taylor]] | |||
* [[Model help:TopoFlow-Evaporation-Read_File]] | |||
* [[Model help:TopoFlow-Meteorology]] | |||
'''''Model Metadata''''' | |||
* [[Model:TopoFlow-Evaporation-Energy_Balance]] | |||
* [[Model:TopoFlow-Evaporation-Priestley_Taylor]] | |||
* [[Model:TopoFlow-Snowmelt-Energy Balance]] | |||
* [[Model:TopoFlow-Meteorology]] | |||
* [[Model:TopoFlow]] | |||
[[Category:Modules]] | [[Category:Modules]] |
Latest revision as of 17:18, 19 February 2018
TopoFlow-Evaporation-Energy Balance
This module is used as the evaporation process component (Energy Balance method) for a D8-based, spatial hydrologic model.
Model introduction
This process component is part of a spatially-distributed hydrologic model called TopoFlow, but it can now be used as a stand-alone model.
Model parameters
Uses ports
• Meteorology
• Channels (surface water flow in a network of channels with trapezoidal cross-section)
• Snow (Snowmelt)
• Infil (Infiltration)
• Satzone (Subsurface flow in saturated zone)
Provides ports
• Evap (Evaporation)
• Configure (tabbed dialog GUI to change settings)
• Run (only if used as the Driver)
Main equations
- Evaporation rate
<math>ET=\left (1000 \ast Q_{et}\right) / \left (\rho_{water} \ast L_{v}\right) </math> (1)
- Energy flux used to evaporate water
<math>Q_{et}=\left (Q_{SW} + Q_{LW} + Q_{c} + Q_{h}\right) </math> (2)
- Conduction energy flux
<math>Q_{c}= K_{soil} \ast \left (T_{soil_x} - T_{surf} \right) \ast \left ( 100 / x \right) </math> (3)
- Sensible heat flux
<math>Q_{h}= \rho_{air} \ast c_{air} \ast D_{h} \ast \left (T_{air} - T_{surf}\right) </math> (4)
- Bulk exchange coeff. (neutrally stable conditions)
<math>D_{n}= u_{z} \ast \kappa^2 / LN [ \left ( z - h_{snow}\right) / z0_{air}] ^2 </math> (5)
- Bulk exchange coeff. for heat
<math>D_{h}=\left\{\begin{matrix} D_{n} / \left (1 + \left (10 \ast Ri \right)\right) & stable: T_{air} > T_{surf} \\ D_{n} \ast [ 1 - \left ( 10 \ast Ri \right) ] & unstable: T_{air} < T_{surf} \end{matrix}\right. </math> (6)
- Richardson's number
<math>Ri= g \ast z \ast \left (T_{air} - T_{surf} \right) / [ u_{z}^2 \left ( T_{air} + 273.15 \right) ] </math> (7)
Symbol | Description | Unit |
---|---|---|
QSW | net shortwave radiation | [W / m2] |
QLW | net longwave radiation | [W / m2] |
Tair | air temperature | deg C |
Tsurf | surface (snow) temperature | deg C |
Tsoil_x | soil temperature at depth x | deg C |
x | reference depth in soil | m |
Ksoil | thermal conductivity of soil | W / (m deg_C ) |
uz | wind velocity at height z | m / s |
z | reference height for wind (above land surface) | m |
z0 | surface roughness height (with no snow) | m |
h0snow | initial snow depth | m |
ρair | density of the air | kg / m3 |
cair | specific heat of air | J /kg |
Lv | latent heat of vaporization, water (2500000) | m / m |
g | gravitational constant, Earth = 9.81 | m / s 2 |
κ | von Karman's constant, equals to 0.41 | - |
ET | evaporation rate | mm / sec |
Qet | energy flux used to evaporate water | W / m2 |
Qc | conduction energy flux (between surf. and subsurf.) | W / m2 |
Qh | sensible heat flux | W / m2 |
Dn | bulk exchange coeff. (neutrally stable conditions) | m / s |
Dh | bulk exchange coeff. for heat | m / s |
Ri | Rechardson's number | - |
ρwater | density of the water | kg / m3 |
Notes
Notes on Input Parameters
If net total radiation has been measured, it can be entered as QSW and then QLW can be set to zero. Any meteorological variables entered here (such as Tair) are automatically shared with other other processes, such as Snowmelt and Precipitation.
Snow depth is tracked internally, and can increase from snowfall or decrease by melting, starting from its initial value. In the current version there is no mechanism for redistribution of snow.
For each 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.
Notes on the Equations
All variables and their units can be seen by expanding the Nomenclature section above.
Wherever (d > 0), evaporation results in a reduction in the surface flow depth. Wherever (d = 0), water is removed from subsurface storage. If the 1D Richards' equation is used for infiltration, then the evaporation rate is applied as a surface boundary condition and alters
In the equation for computing Dn, the reference height, z, is reduced by the computed snow depth, hsnow. It follows that z must be chosen so as to be greater than any possible snow depth.
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: https://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
Brutsaert, W. (1975) On a derivable formula for long-wave radiation from clear skies, Water Resources Research, 11, 742-744.
Dingman, S.L (2002) Physical Hydrology, 2nd ed., Prentice Hall, New Jersey. (see Chapter 7, pp. 285-299)
Schlicting, H. (1960) Boundary Layer Theory, 4th ed., McGraw-Hill, New York, 647 pp.
Zhang, Z., D.L. Kane and L.D. Hinzman (2000) Development and application of a spatially-distributed Arctic hydrological and thermal process model (ARHYTHM), Hydrological Processes, 14, 1017-1044.
Links
Related Help Pages
- Model help:TopoFlow-Evaporation-Priestley_Taylor
- Model help:TopoFlow-Evaporation-Read_File
- Model help:TopoFlow-Meteorology
Model Metadata