HYPE Model Documentation

This page is read only. You can view the source, but not change it. Ask your administrator if you think this is wrong.
====== Processes above ground ======

===== Temperature and precipitation =====

HYPE has many possibilities to adjust the precipitation and temperature given in the Pobs.txt and Tobs.txt files. Adjustment is done because the supplied data can be of different reference height or (known to be) wrong in some areas, or because data vary over the subbasin but the input is given as one value, and or to compensate for non modelled processes.

==== Temperature adjustments ====
Subbasin air temperatures (<m>T_i</m>) in the input file are normally assumed to be for the subbasin’s average elevation (//elev//). The same adjustments are made for maximum and minimum temperatures if available.
 
(1) If data with other reference height is used it is possible to adjust for that. The subbasin temperature can be adjusted depending on subbasin elevation with general parameter //tcelevadd// (//° / 100m//). This adjustment assumes observations are located at sea level. 

(2) Alternatively it can be adjusted based on the difference between subbasin elevation and temperature observation elevation by general parameter //tcobselev// (//° / 100m//). The temperature observation elevation (<m>elev_obs</m>) then need to be supplied to the model. 

(3) It can be adjusted based on month with the monthly dependent parameter //monthlapse// (//° / 100m//)). These three subbasin elevation adjustments should not be used together. 
Subbasin temperature can be adjusted equally over all subbasins within a region with the parameter region dependent parameter //tempcorr//.

<m> T_gc = T_i + tempcorr - tcelevadd * elev/100 - tcobselev * (elev - elev_obs)/100 - monthlapse*elev/100 </m>

The temperature can also be adjusted for each class depending on their deviation from the subbasin average elevation (<m>{Delta}h</m>). The class-dependent temperature (//T//) is calculated using the parameter //tcalt//. The temperature lapse rate often has a value of 0.6 (//° / 100m//). 

<m> T = T_gc - tcalt*{{Delta} h}/100 </m>

==== Precipitation adjustments ====
Subbasin input precipitation (<m>P_i</m>) can be adjusted equally over all subbasins with the general parameter //pcaddg// or for some subbasins with the parameter region dependent parameter //preccorr//. Additionally it is possible to adjust precipitation for undercatch with different parameters (//pcurain,pcusnow//) depending on if the precipitation falls as snow or rain. 

Subbasin precipitation (<m>P_gc</m>) is for subbasin average elevation (//elev//), but can be adjusted for elevation variations within the subbasin. The precipitation of a class (//P//) is adjusted for classes where the class average elevation is greater than a threshold (general model parameter //pcelevth//). The adjustment is determined by a general parameter (//pcelevadd//) that is the correction per 100m. The class elevation adjustment can alternatively be determined from the basin standard deviation of elevation (<m>elev_std</m>) and a parameter //pcelevstd//. The class height adjustment is limited by a general parameter //pcelevmax//. The precipitation of a class can additionally be adjusted with land-use dependent parameter //pcluse//, e.g. for interception evaporation.  

<m> P_gc=P_i*(1+pcaddg)*(1+preccorr)*(1+(pcurain*(1-snfrac)+pcusnow*snfrac)) </m>

<m> pc_{height}=delim{lbrace}{ 
matrix{2}{2}{
    0 {elev+{Delta}h<pcelevth}
    {MIN({{elev+{Delta}h-pcelevth}/{100}}*pcelevadd+{{elev_std}/{100}}*pcelevstd,pcelevmax)}  else
   }}{} </m>

<m> P=P_gc*(1+pc_{height} )*(1-pcluse) </m>

Where <m>{Delta}h</m> is a class's elevation deviation from the subbasin average elevation and //snfrac// is the average fraction of precipitation that falls as snow calculated from subbasin temperature (<m>T_i</m>) and temperature thresholds (see equation below), or from input.
==== Rainfall and snowfall separation ====
The rain/snow fraction of precipitation is calculated based on temperature or given as an input time series. Different temperatures can be used in the equation, i.e. basin average or class temperature. When the air temperature (//T//) is around the threshold temperature for mixed precipitation (land-use dependent parameter //ttmp// plus general parameter //ttpd//) both rain and snow. The interval for mixed precipitation is given by the parameter //ttpi//. For temperature below threshold minus //ttpi//, the precipitation is assumed to be in solid form only and is added to the snowpack. If the air temperature is greater than the threshold temperature plus //ttpi//, the precipitation is assumed to be solely in liquid form. For intermediate temperatures, the precipitation is assumed to be a mixture of liquid and solid forms i.e. as both rain and snow. The proportion (<m>a_rain</m>) of precipitation that falls as rain depends linearly on the temperature. 


<m> a_{rain}=delim{lbrace}{ 
matrix{3}{2}{
    0 {T<ttmp+ttpd-ttpi}
    1 {T>ttmp+ttpd+ttpi}
    {{(T-(ttmp+ttpd-ttpi))}/{(2*ttpi)}}  else
   }}{} </m>

Alternatively snowfall fraction (//sffrac//) may be read from input file and <m>a_rain</m> calculated as: 

<m> a_{rain}=1-sffrac </m>

The amount of rainfall and snowfall is calculated from the corrected precipitation (//P//).

<m> rainfall=P*a_{rain} </m>
<m> snowfall=P*(1-a_{rain} ) </m>

==== Snowfall distribution ====

The optional snowfall distribution models redistribute the snowfall of a catchment among the classes depending on wind. The model uses Winstral coefficients. For each subbasin and class and wind direction a Winstral coefficient (//WSF//) is given as input data to the model. Wind direction is calculated from forcingdata of westerly (U) and southerly (V) wind. The wind direction is translated to the closest quadrant or octant, and the Winstral coefficient of that direction is used.

The Winstral coefficient is, with some scaling parameters, used to calculate the relative snowfall for each class this time step. There are two models; one that use a linear function (ref Gustafsson et al, 2015) and one that is log-linear (ref Clemenzi et al., 2020). 

<m> relative_snowfall = 1. + WSFscale *WSFluse * (WSF + WSFbias) </m> 

<m> relative_snowfall = 10 ^ (WSFscale * WSFluse * (WSF + WSFbias)) </m> 

//WSFscale// and //WSFbias// are general parameter, //WSFluse// is a scaling parameter depending on land-use. The relative snowfall is limited to the interval [1 - //sfdmax//, 1 + //sfdmax//] for the linear model and [0,//sfdmax//] for the log-linear. //sfdmax// is a general parameter. The relative snowfall is normalized within each subbasin so that the subbasin mean snowfall is preserved. The snowfall distribution is applied only if snow is falling for a fraction of the subbasin area that is larger than a limit (//limit//).

NOTE: The snowfall distribution model only works for simulation of the winter period.


==== Links to file reference ====

^Section	^Symbol ^Parameter/Data ^File ^
|Temperature adjustments|<m>T_i</m>| |[[start:hype_file_reference:tobs.txt|Tobs.txt]]|
|:::|//elev//|//elev_mean//|[[start:hype_file_reference:geodata.txt|GeoData.txt]]|
|:::|<m>{Delta}h</m>|//dhslc_nn//|:::|
|:::|<m>elev_obs</m>|//tobselev//|[[start:hype_file_reference:forckey.txt|ForcKey.txt]]|
|:::| |//tcelevadd, tcobselev, monthlapse, tempcorr, tcalt//|[[start:hype_file_reference:par.txt|par.txt]]|
|Precipitation adjustments|<m>P_i</m>| |[[start:hype_file_reference:pobs.txt|Pobs.txt]]|
|:::| |//pcaddg, preccorr, pcurain, pcusnow, pcelevth, pcelevadd, pcelevstd, pcelevmax, pcluse//|[[start:hype_file_reference:par.txt|par.txt]]|
|:::|<m>{Delta}h</m>|//dhslc_nn//|[[start:hype_file_reference:geodata.txt|GeoData.txt]]|
|:::|<m>elev_std</m>|//elev_std//|:::|
|:::|//snfrac//|calculated based on <m>T_i</m> or from |[[start:hype_file_reference:sfobs.txt|SFobs.txt]]|
|Rainfall and snowfall calculations| |//ttmp, ttpd, ttpi//|[[start:hype_file_reference:par.txt|par.txt]]|
|:::|//sffrac//| |[[start:hype_file_reference:sfobs.txt|SFobs.txt]]|
|Snowfall distribution|U, V | |[[start:hype_file_reference:uwobs.txt|UWobs.txt]], [[start:hype_file_reference:vwobs.txt|VWobs.txt]]|
|:::|//WSF//|//wsf_nn_d// |[[start:hype_file_reference:geodata.txt|GeoData.txt]]|
|:::| |//WSFscale,WSFluse,WSFbias,sfdmax// |[[start:hype_file_reference:par.txt|par.txt]]|
|:::|//limit// |//sfdlim// |:::|






==== Links to relevant procedures in the code ====  

^ Modules (file)      ^ Procedures      ^
| [[http://hype.sourceforge.net/doxy-html/namespacemodelmodule.html|modelmodule (model_hype.f90)]]    | model     |
| [[http://hype.sourceforge.net/doxy-html/namespacesoil__processes.html|soil_processes (soil_proc.f90)]]| calculate_snow |
| ::: | infiltration |
|[[http://hype.sourceforge.net/doxy-html/namespaceatmospheric__processes.html|atmospheric_processes (atm_proc.f90)]]	| calculate_class_atmospheric_forcing |
| ::: | calculate_rain_snow_from_precipitation |
| ::: | calculate_subbasin_temperature |
| ::: | calculate_subbasin_precipitation |
| ::: | set_atmospheric_parameters_corrections |
| ::: | calculate_winddirspeed |
| ::: | calculate_snowfall_distribution |

===== Evaporation =====

Potential evaporation (//epot// in mm) is calculated based on the temperature if it is not read in from file (Xobs.txt). Alternative PET models exist, and is described [[start:hype_model_description:processes_above_ground#alternative_potential_evaporation_models|below]]. When the air temperature (//T//) is greater than the threshold temperature //ttmp// evaporation is assumed to occur. Snow melting, snow density and evaporation use the same threshold temperature. The basic potential evapotranspiration (<m>epot_base</m>) depends on the land use dependent rate parameter //cevp//.

<m> cseason = 1 + cevpam*sin(2*pi*(dayno-cevpph)/365) </m>

<m> epot_{base} = (cevp * cseason) * (T-ttmp) </m>

A seasonal factor //cseason// adjusts the potential evaporation rate (//cevp//) e.g. making it higher in the spring when the air is often dry, and lower in autumn when the air is often more humid than in spring. The factor is sinusoidal with two parameters //cevpam// and //cevpph//. It is not used if //cevpam// is zero. A //cevpph// around 45 days give a maximum correction in mid May (dayno=45+91=136). The minimum correction will then be a half year later in September (dayno=136+182). For an earlier maximum, reduce //cevpph//. 

The basic potential evapotranspiration may be adjusted with a regional correction factor (//cevpcorr//) equally over the year depending on parameter region.

<m> epot = epot_{base} * (1 + cevpcorr) </m>

Evaporation from soil is assumed to occur from the two upper layers. The potential evaporation is assumed to decrease exponentially with depth (depending on the parameter //epotdist//). The potential evaporation is divided between the two layers (//epotfrac//) with the distribution depending on the potential evaporation in the midpoint of each soil layer (figure 1). This is then used by approximating to a rectangle. Since soil layers differ between classes, the evaporation distribution do to.

<m> epot1 = EXP(- epotdist*soillayerdepth(1){/}2) </m>

<m> epot2 = EXP(- epotdist*(soillayerdepth(1)+{soillayerdepth(2)-
          soillayerdepth(1)}/2)) </m>
          
<m>  area1 = soillayerdepth(1)*epot1 </m>

<m>  area2 = (soillayerdepth(2)-soillayerdepth(1))*epot2 </m>

<m> epotfrac1 = area1 / {area1 + area2} </m>

<m>  epotfrac2 = area2 / {area1 + area2} </m>



|{{:start:hype_model_description:potentialevaporation.png?400|}}|
|Figure 1 The distribution of potential evaporation between the top two soil layers.|


The actual evaporation from a soil layer (//evap//) is limited by the availability of water in the soil (//soil//) above the wilting point (//wp//, mm). Evaporation is at potential rate only if the water exceeds field capacity (//fc//, mm) or a (large) proportion (general parameter //lp//) of field capacity. In between these limits evaporation increase linearly.

<m> evapp= {lbrace}{
 \matrix{3}{2}
    {0 {soil-wp<0}
    epot*epotfrac {soil-wp>lp*fc}
    epot*epotfrac*{{soil-wp}/{lp*fc}} {else} }} </m>

<m> evap = MAX(evapp,soil-wp) </m>


|{{:start:hype_model_description:evap.png?400|}}|
|Figure 2 Evaporation as a function of soil water.|

The actual evaporation may also depend on soil temperature (//soiltemp//). It is then reduced for temperatures above land use parameter //ttrig// and depend on two other land use parameters (//tredA//, //tredB//) as well (Figure 3). 

|{{:start:hype_model_description:pet_soiltemp.png?400|}}|
|Figure 3 Soil temperature factor for reduction of soil evapotranspiration. Parameter values: //ttrip//=1, //tredA//=0.5, //tredB//=1.|

The soil temperature evapotranspiration reduction is calculated as:

<m> factor = 1-e^( - tredA*(soiltemp-ttrig)^tredB) </m>

<m> evapp = evapp*factor </m>

The actual soil evaporation is set to zero for temperatures below the threshold temperature and for negative potential evaporation estimates (condensation). It may also be affected by frozen soil model, which then limit evaporation to the liquid part of soil water.

A river with an area (i.e. has a class), flooded floodplains and lakes are assumed to evaporate at the potential rate, when the air temperature is above the threshold temperature (//ttmp//). Evaporation is limited by the water body's volume.


==== Alternative potential evaporation models ====

HYPE give the option to exchange the default potential evaporation model for another model. Only the basic potential evaporation (<m>epot_base</m>) differs between models. Models 0-2 only use air temperature forcing. Models 3-5 want shortwave radiation and minimum and maximum daily air temperature, although if lacking approximations are made (see section [[start:hype_model_description:processes_above_ground#input_to_alternative_potential_evaporation_models|Input to alternative potential evaporation models]] below). In addition, model 5 wants relative humidity and wind speed if available.
Note, regardless of potential evaporation model, the actual evaporation is limited to temperatures above parameter //ttmp//.

=== Model 0 (default) === 
As described above; evapotranspiration depends on the rate parameter //cevp// and air temperature (//T//) above a threshold //ttmp//. If //epot// is given in Xobs.txt those values are used.

<m> epot_{base} = (cevp * cseason) * (T-ttmp) </m>

=== Model 1 ===
Model 1 is the same as model 0, but it will not be using //epot// from input data, even if it is present.

=== Model 2 - Modified Jensen-Haise/McGuiness ===
The modified Jensen-Haise/McGuiness model follow Oudin et al. (2005). The potential evaporation depends on extraterrestrial radiation (//radext//), latent heat of vaporization (<m>lambda</m>) and temperature (//T//). Two general parameters (//jhtadd// and //jhtscale//) are used and one land use dependent (crop coefficient //kc2// or //kc//).

<m> epot_{base} = {kc/jhtscale} * MAX(0,{radext/lambda}*(T+jhtadd)) </m>

=== Model 3 - Modified Hargreaves-Samani ===
The Hargreaves-Samani evaporation is modified to limit the "turbidity-factor". The potential evaporation depends on extraterrestrial radiation (//radext//), latent heat of vaporization (<m>lambda</m>), temperature (//T//) and turbidity (//turbidity//). One general parameter (//krs//) is used and one land use dependent (crop coefficient //kc3// or //kc//).

<m> epot_{base} = MAX(0,kc*0.0023 * {radext/lambda}*{turbidity/krs}*(T+17.8)) </m>

=== Model 4 - Priestly-Taylor ===
The Priestly-Taylor potential evaporation depends on net downward radiation (//netrad//), slope of saturated vapour pressure curve (//dsatvap//), latent heat of vaporization (<m>lambda</m>) and a psychrometric constant (<m>gamma</m>). One general parameter (//alfapt//) are used and one land use dependent (crop coefficient //kc4// or //kc//).

<m> epot_{base} = MAX(0,kc*alfapt * {dsatvap*netrad/lambda*(dsatvap+gamma)}) </m>

=== Model 5 - FAO Penman-Monteith ===
The FAO Penman-Monteith potential evaporation depends on net downward radiation (//netrad//), slope of saturated vapour pressure curve (//dsatvap//), saturated and actual vapour pressure (//satvap// and //actvap//), temperature (//T//), wind speed (//wind//) and a psychrometric constant (<m>gamma</m>). One land use dependent parameter (crop coefficient //kc5// or //kc//) is used.

<m> epot_{base} = MAX(0,kc*{0.408 * dsatvap*netrad + gamma*{{900}/{T+273}}*wind*(satvap-actvap)}/{dsatvap+gamma*(1+0.34*wind)}) </m>

==== Input to alternative potential evaporation models ====

Summary of alternative input to PET models, and link to file reference.
^ Model ^ Parameters ([[start:hype_file_reference:par.txt|par.txt]])^ Static data ([[start:hype_file_reference:geodata.txt|GeoData.txt]])^ Forcing data ([[start:hype_file_reference#observation_data_files|files]])^
| 0 | cevp, ttmp | | Xobs:repo, Tobs |
| 1 | cevp, ttmp | | Tobs |
| 2: modified Jensen-Haise/McGuiness | jhtadd, jhtscale, kc | latitude | Tobs |
| 3: modified Hargreaves-Samani | kc, krs | elevation, latitude | SWobs, Tobs |
| 3: modified Hargreaves-Samani | kc, krs | elevation, latitude | TMINobs, TMAXobs, Tobs |
| 4: Priestly-Taylor | alb, alfapt, kc | elevation, latitude | RHobs, SWobs, TMINobs, TMAXobs, Tobs |
| 4: Priestly-Taylor | alb, alfapt, kc, krs | elevation, latitude | SWobs, TMINobs, TMAXobs, Tobs |
| 4: Priestly-Taylor | alb, alfapt, kc, krs | elevation, latitude | TMINobs, TMAXobs, Tobs |
| 5: FAO Penman-Monteith | alb, kc, roughness, zphd, zwind, zwish | elevation, latitude | RHobs, SWobs, TMINobs, TMAXobs, Tobs, Uobs |
| 5: FAO Penman-Monteith | alb, kc, mwind | elevation, latitude | RHobs, SWobs, TMINobs, TMAXobs, Tobs |
| 5: FAO Penman-Monteith | alb, kc, krs, roughness, zphd, zwind, zwish | elevation, latitude | SWobs, TMINobs, TMAXobs, Tobs, Uobs |
| 5: FAO Penman-Monteith | alb, kc, krs, mwind | elevation, latitude | SWobs, TMINobs, TMAXobs, Tobs |
| 5: FAO Penman-Monteith | alb, kc, krs, roughness, zphd, zwind, zwish | elevation, latitude | TMINobs, TMAXobs, Tobs, Uobs |
| 5: FAO Penman-Monteith | alb, kc, krs, mwind | elevation, latitude | TMINobs, TMAXobs, Tobs |


=== Actual vapour pressure ===
Actual vapour pressure (//actvap//) is calculated following FAO recommended procedure and function/data priority. Depending on availability of minimum, mean and maximum relative humidity (//rh//) and minimum, mean and maximum air temperature (<m>T_min, T, T_max</m>) equations with different combinations of saturated vapour pressure (calculated from temperature) times relative humidity is used. For example:

<m> actvap = {satvap(T_max)*rh_min+satvap(T_min)*rh_max}/2 </m>

<m> actvap = satvap(T_min)*rh_max </m>

<m> actvap = satvap(T)*rh_mean </m>

<m> actvap = satvap(T_min)*1 </m>

In case not enough data is available, minimum temperature is calculated from turbidity and a general parameter (//krs//) and the last equation of those above is used.

<m> T_min = T - 0.5*(turbidity/krs)^2 </m>

Eventually actual vapour pressure is limited by the calculated saturated vapour pressure. 

=== Air pressure ===
Air pressure (//pa//) is calculated as a function of elevation for the class (//elev//).

<m> pa = 101.3*({293-0.0065*elev}/{293})^5.26 </m>

=== Shortwave radiation ===
Shortwave radiation (//swrad//) is supplied as input forcing time series or otherwise calculated from extraterrestrial radiation. 

<m> swrad=radext*turbidity </m>


=== Extraterrestrial radiation ===

Extraterrestrial solar radiation (//radext//) is estimated from day of year and latitude. The equations used comes from FAO. The day of the year are used to estimate distance to the sun and declination. The sunset hour angle are calculated from latitude, but with special care for high latitudes (polar night and midnight sun). These variables together with latitude and the solar constant eventually give the current extraterrestrial radiation for each subbasin.

=== Latent heat of vaporization ===
Latent heat of vaporization (<m>lambda</m>) is a function of temperature (//T//).

<m> lambda = 2.501-0.002361*T </m>

=== Net downward radiation ===
The net downward radiation (//netrad//) is used explicitly for PET model 4 and 5 (Priestly-Taylor and FAO Penman-Monteith). The net radiation is calculated following FAO recommended procedure. It is calculated as net shortwave radiation minus net longwave radiation. Net shortwave radiation (<m>net_short</m>) is calculated from the shortwave radiation (//swrad//) and the land use dependent albedo parameter (//alb//). Net longwave radiation (<m>net_long</m>) is calculated using temperature (<m>T_min, T_max, T</m>), actual vapour pressure (//actvap//) and relative shortwave radiation (//relsh//) if those are available, otherwise it is set to zero. The relative shortwave radiation is shortwave radiation in relation to clear sky shortwave radiation.

<m> net_short = swrad*(1-alb) </m>

<m> net_long = 4.903*10^9*{{(T_max+273.15)^4+(T_min+273.15)^4}/2}*(0.34-0.14*actvap^0.5)*(1.35*relsh-0.35) </m>

<m> relsh = turbidity/clearturb </m>

=== Psychrometric constant ===
The psychrometric constant (<m>gamma</m>) is a function of air pressure (//pa//) and latent heat of vaporization (<m>lambda</m>) following FAO.

<m> gamma = 0.001013*pa/{0.622*lambda} </m>

=== Saturated vapour pressure ===
Saturated vapour pressure (//satvap//, kPa) is calculated from temperature (//temp//) following FAO. If daily minimum and maximum air temperature is available //satvap// is set to the average of the saturated vapour pressure calculated for each of those two temperatures, otherwise it is calculated from daily average air temperature.

<m> satvap=0.6108*{EXP({17.27*temp}/{temp+237.3})} </m>

=== Slope of saturated vapour pressure curve ===
The slope of saturated vapour pressure temperature function (//dsatvap//) is used explicitly for PET model 4 and 5 (Priestly-Taylor and FAO Penman-Monteith). The slope is calculated from daily mean air temperature (//T//).

<m> dsatvap=4098*0.6108*{EXP({17.27*T}/{T+237.3})}*{{1}/{(T+237.3)^2}} </m>

=== Turbidity ===

The turbidity factor is used explicitly for the PET model 3 - modified Hargreaves-Samani. If shortwave radiation (//swrad//) has been given as forcing data time serie, the turbidity factor is calculated as 

<m> turbidity = swrad/radext </m>

but limited by a minimum turbidity value (0.25) and an estimated clearsky turbidity. The clearsky turbidity (//clearturb//) is estimated by the Ångström formula (FAO): 

<m> clearturb = 0.75+elev*0.00002 </m>

where //elev// is subbasin elevation in meter above sea level. If no shortwave radiaton time series are given, but time series of daily minimum and maximum temperature are, the turbidity is calculated as in the "ordinary" Hargreaves-Samani turbidity function:

<m> turbidity = krs*SQRT(T_max - T_min) </m>

but still limited by the minimum turbidity value and the calculated clearsky turbidity.

=== Wind speed ===
Wind speed is used in PET model 5 - FAO Penman-Monteith. Wind speed (//wind//) may be given as a constant general parameter (//mwind//) or as a forcing data time series. The time serie wind is given for each subbasin (<m>U_i</m>). It is possible to adjust the time serie wind speed to different height than observations. If the general parameters //zwind//, //zwish//, //roughness//, and //zpdh// is set, wind speed is adjusted with the transformation factor //windtrans//.

<m> wind=U_i*windtrans </m>

<m> windtrans = {ln(zwind-zpdh)-ln(roughness)}/{ln(zwish-zpdh)-ln(roughness)} </m>

==== Links to file reference ====

^Section	^Symbol ^Parameter/Data ^File ^
|Evaporation| |//cevp, ttmp, cevmam, cevpph, cevpcorr, epotdist, lp, ttrig, tredA, tredB// |[[start:hype_file_reference:par.txt|par.txt]]|
|:::|//wp//|calculated from //wcwp, wcwp1, wcwp2, wcwp3// and soillayerdepth|:::|
|:::|//fc//|calculated from //wcfc, wcfc1, wcfc2, wcfc3// and soillayerdepth|:::|
|:::|//soillayerdepth//| |[[start:hype_file_reference:geoclass.txt|GeoClass.txt]]|
|Alternative potential evaporation models|//epot//| |[[start:hype_file_reference:xobs.txt|Xobs.txt]]|
|:::| |//jhtadd, jhtscale, kc, kc2, kc3, kc4, kc5, krs, alfapt//|[[start:hype_file_reference:par.txt|par.txt]]|
|Input to alternative potential evaporation models| |//alb, mwind, zwind, zwish, roughness, zpdh//|[[start:hype_file_reference:par.txt|par.txt]]|
|:::|<m>U_i</m>| |[[start:hype_file_reference:uobs.txt|Uobs.txt]]|
|:::|//swrad//|calculated or from |[[start:hype_file_reference:swobs.txt|SWobs.txt]]|
|:::|//elev//|//elev_mean//|[[start:hype_file_reference:geodata.txt|GeoData.txt]]|
|:::| |//latitude//|:::|


==== Links to relevant procedures in the code ====

^ Modules (file) ^ Procedures ^
| [[http://hype.sourceforge.net/doxy-html/namespacesoil__processes.html|soil_processes (soil_proc.f90)]] | initiate_soil_water |
| ::: |calculate_potential_evaporation |
| ::: |calculate_actual_soil_evapotranspiration |
| [[http://hype.sourceforge.net/doxy-html/namespacesurfacewater__processes.html|surfacewater_processes (sw_proc.f90)]] | calculate_river_evaporation |
| ::: |calculate_actual_lake_evaporation |

==== References ====

Richard G. Allen, Luis S. Pereira, Dirk Raes & Martin Smith, 1998. Crop Evapotranspiration – Guidelines for Computing Crop Water Requirements. FAO Irrigation and drainage paper 56. Rome, Italy: Food and Agriculture Organization of the United Nations. ISBN 92-5-104219-5

Oudin, L., F. Hervieu, C. Michel, C. Perrin, V. Andreassian, F. Anctil and C. Loumagne 2005. Which potential evapotranspiration input for a lumped rainfall–runoff model?: Part 2—Towards a simple and efficient potential evapotranspiration model for rainfall–runoff modelling, Journal of Hydrology 303(1-4):290-306.


===== Atmospheric deposition of nitrogen and phosphorus =====

==== Wet deposition ====
Atmospheric deposition in the form of wet deposition of IN and SP is added as a concentration of rainfall. You can specify a time series of the concentration in precipitation in Xobs.txt (implemented for IN, SP, and T1). Otherwise nitrogen wet deposition (as a concentration) is specified in GeoData.txt for each subbasin, while the wet deposition of phosphorus is specified by a general model parameters (//wetdepsp//). 

If parameter (//aloadconst//) is set, the wet deposition load in considered constant and not dependent on precipitation corrections. Thus the concentration of precipitation is changed to keep the load constant, when precipitation correction is applied. 



==== Atmospheric deposition to the soil ====

Nitrogen dry deposition is specified in GeoData.txt for each subbasin and different vegetation groups, while dry deposition of phosphorus is specified by a model parameter //drydepPP// (land use dependent). 

Dry deposition of nitrogen and phosphorus is added to the snow or, if there is no snow, directly to the ground. Phosphorus deposition is added to PartP-pool in the upper soil layer and nitrogen deposits to the IN dissolved in soil water if the parameter //ponatm// is not set. A new concentration of IN in the soil water is then calculated. The parameter //ponatm// indicates that some of the nitrogen deposition will be added to organic nitrogen pool (//fastN//) instead of to the IN pool.

For both nitrogen and phosphorus wet deposition is added through the concentration of precipitation. Precipitation may fall as snow and be added to the snow pack or if it is rain added to the potential infiltration. Depending on the fate of the infiltration, the wet atmospheric deposition will mix with the recieving water (e.g. surface runoff, soil layer water, macropore flow).
==== Atmospheric deposition to rivers and lakes ====

Nitrogen dry deposition on water surfaces is specified in GeoData.txt for each subbasin, while dry deposition of phosphorus is specified by a model parameter //drydepPP// (land use dependent). 

For both nitrogen and phosphorus wet deposition is added through the concentration of precipitation (as described above), while for dry deposition an amount of the nutrient is added to the river (if it has a class-area) or lake water.

Alternatively wet deposition of phosphorus on water surfaces can be specified by a general model parameter (//wetdepspl//) as a load. Monthly load of IN atmosperic depostion on water surfaces can be specified in GeoData.txt for each subbasin.

==== Links to file reference ====

^Section ^Parameter/Data ^File ^
|Wet deposition|//cpT1, cpIN, cpSP//|[[start:hype_file_reference:xobs.txt|Xobs.txt]]|
|:::|//wetdep_n//|[[start:hype_file_reference:geodata.txt|GeoData.txt]]|
|:::|//wetdepsp, aloadconst//|[[start:hype_file_reference:par.txt|par.txt]]|
|Atmospheric deposition to the soil|//drydep_n1,drydep_n2//|[[start:hype_file_reference:geodata.txt|GeoData.txt]]|
|:::|//drydepPP, ponatm//|[[start:hype_file_reference:par.txt|par.txt]]|
|Atmospheric deposition to rivers and lakes|//drydep_n3, deploadn1-deploadn12//|[[start:hype_file_reference:geodata.txt|GeoData.txt]]|
|:::|//drydepPP, wetdepspl//|[[start:hype_file_reference:par.txt|par.txt]]|

==== Links to relevant procedures in the code ====

^Modules (file)	^Procedures ^
|[[http://hype.sourceforge.net/doxy-html/namespacenpc__soil__processes.html|npc_soil_processes (npc_soil_proc.f90)]]|set_class_precipitation_concentration_and_load|
|:::|add_dry_deposition_to_landclass|
|:::|atmdep_in_loss|
|[[http://hype.sourceforge.net/doxy-html/namespaceatmospheric__processes.html|atmospheric_processes (atm_proc.f90)]]|set_precipitation_concentration|
|[[http://hype.sourceforge.net/doxy-html/namespacenpc__surfacewater__processes.html|npc_surfacewater_processes (npc_sw_proc.f90)]]|add_dry_deposition_to_river|
|:::|add_dry_deposition_to_lake|
HYPE Model Documentation

User Tools

Site Tools

Page Tools
