View source for Energy Balance

An energy balance (EB) method is one of the options included in the GSSHA model to estimate snowfall melting.  This method is admittedly simple and other factors, soil freezing, change in overland roughness, etc, are not yet considered.  This is an area of active research and model development at ERDC.  As inputs, the EB method requires hourly hydrometeorological data values of air temperature (''T<sub>a</sub>''), relative humidity (''r<sub>h</sub>''), wind speed (''U'') barometric pressure (''P<sub>a</sub>'') and cloud cover.  

In the energy budget model the amount of heat available is applied to the snowpack and the amount of meltwater is calculated.  The simplest representation of the snowpack is used; each 80 calories of heat added to the snowpack results in the release of 1 cm<sup>3</sup> of meltwater (Linsley et al., 1982, Gray and Prowse, 1993).  This method ignores complex snowpack behavior, such as ripening of the snowpack and refreezing of meltwater.  Hourly values of hydrometeorological variables allow both seasonal and diurnal variations in climatic conditions to be included in the heat balance. 

The amount of heat, Q (cal cm<sup>-2</sup> hr<sup>-1</sup>) available is computed from the components of the energy balance.  In GSSHA the following components are accounted for:

''Q<sup>*</sup> - net radiation (in - out), 
''Q<sub>v</sub>'' - heat in precipitation, 
''Q<sub>e</sub>'' - heat transferred by sublimation and evaporation, and  
''Q<sub>h</sub>'' - sensible heat transfer due to turbulence.

For non-precipitation periods the net radiation is typically the dominant source of energy for melting of the snowpack (Gray and Prowse, 1993).  The net radiation is computed using Stephan-Boltzman’s law, with the assumptions that incoming radiation can be computed from the ambient temperature, ''T<sub>a</sub>''  (C), and outgoing radiation is computed assuming the snowpack is at 0&deg; C (Bras, 1990):

{| 
|- 
| 
: 
| width=550 | ''Q<sup>*</sup> = 49.56 x 10-10(''T<sub>a</sub>'' + 273)4 – 27 || (92)
|}

Precipitation falling on the snowpack at temperatures above 0&deg; transmits the difference in heat between the raindrop and the snowpack.  Assuming the snowpack is at 0&deg; C and the rainfall is a ambient temperature the difference in heat energy is:

{| 
|- 
| 
: 
| width=550 | ''Q<sub>v</sub>'' = ''IT<sub>a</sub>'' || (93)	
|}

where:  ''I'' is the precipitation intensity (cm/hr).  Heats transferred from evaporation, sublimation, and turbulent energy are usually much smaller parts of the heat balance and are ignored in many computations (Gray and Prowse, 1993).  However, convective exchange can be significant (Linsley et al., 1982).  If the dew point is below the temperature of the snowpack, assumed to be 0&deg; C, then condensation occurs and heat is transferred (Linsley et al., 1982; Gray and Prowse, 1993).  Estimates of turbulent and latent heat exchange are usually based on measurements of air temperature, humidity, and wind speed (Gray and Prowse, 1993).  During periods of melt, the temperature of the snowpack is 0&deg; C and the saturated vapor pressure (es) is 6.11 mb (Linsley et al., 1982).  The latent heat exchange is computed assuming the latent heat of evaporation/condensation is 600 cal g<sup>-1</sup> (Anderson, 1968) and a water density of 1 g cm<sup>-3</sup> as:

{| 
|- 
| 
: 
| width=550 | [[Image:Equation094.gif]] || (94)	
|}

where:  ''r<sub>h</sub>'' is the relative humidity (%), &fnof;(V) = 0.0002 ''U'' (km/hr) (Anderson 1978), where ''U''U is the wind speed (m s<sup>-1</sup>).  Employing the Bowen ratio (Bowen, 1926) the sensible heat transfer is computed assuming the snow pack temperature is at 0&deg; C, latent heat of evaporation is 600 cal g<sup>-1</sup>, density of water is 1 g cm<sup>-3</sup>, and the Bowen ratio coefficient is 0.61 x 10-3 C<sup>-1</sup> (Bras, 1990) as: 

{| 
|- 
| 
: 
| width=550 | [[Image:Equation095.gif]] || (95) 
|}

where:  ''P<sub>a</sub>'' is the atmospheric pressure (mb).

The total energy to melt snow is calculated using these energy fluxes as shown in Equation 96.  The amount of water melted is then calculated using the latent heat of melting (~80 cal cm-3) as shown in Equation 97.

{| 
|- 
| 
: 
| width=550 | ''Q<sub>melt</sub> = Q<sub>a</sub> - Q<sub>bs</sub> + Q<sub>e</sub> + Q<sup>h</sup> + Q<sub>p</sub> || (96)
|}

{| 
|- 
| 
: 
| width=550 | ''M<sub>EB</sub> = [(Q<sub>melt</sub> / 100) / 80.0 ] * dt || (97)
|}

	Q<sub>melt</sub> = total energy available to melt snow (cal cm<sup>-2</sup> hr<sup>-1</sup>)
	M<sub>EB</sub> = melt calculated using energy balance snow melt routine (mm SWE)


For non-precipitation periods, the energy budget is calculated at an hourly time step (same as the standard hydrometeorological data), so diurnal changes in energy inputs are included in the model formulation.   During precipitation periods the energy budget is updated each overland flow routing time step (generally less than 5 minutes).