User Tools

Site Tools


start:hype_model_description:hype_land

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revision Previous revision
Next revision
Previous revision
start:hype_model_description:hype_land [2018/10/15 14:04]
cpers [Links to file reference]
start:hype_model_description:hype_land [2019/03/19 14:19] (current)
cpers [Diagnostic variables]
Line 38: Line 38:
 Some additional output variables are calculated from the soil state variables. ​ Some additional output variables are calculated from the soil state variables. ​
  
-=== Ground water level ===+=== Groundwater ​level ===
  
-A water table is calculated for each soil layer (//gwat//from the proportion ​of water-filled pores of effective porosity. If the ground water reaches above the surface, the water is calculated with 100% porosity.+The groundwater level is measured negative from surface ​(0mto bottom ​of the soil layers. A positive groundwater level means that the soil surface is below water. If the ground water table reaches above the surface, the water is calculated with 100% porosity.
  
-  IF(soil(k)-wp(k)-fc(k)>​0.0) +The water table is found in the lowest ​soil layer that is not completely filled with water. Soil layers above this layer may have water in its effective porosity, but that is not included in the groundwater level output variable.  
-    ​gwat(k) = (soil(k)-wp(k)-fc(k)-ep(k))/​ep(k) * soillayerthick(k) –  +The water table for a soil layer is calculated linearly from the proportion of water-filled pores of effective porosity part of the soil pore volume. If the soil moisture of a soil layer is at field capacity ​(or below), the groundwater level of that soil layer is at the bottom of the layer. If the pore volume is filled, the groundwater level of that soil layer is at the top of the layer
-               ​soillayerdepth(k-1) +
-  IF(gwat(1) > 0) gwat(1) = (soil(1)-wp(1)-fc(1)-ep(1))/​1000.+
  
-The water table measured as a negative from ground surface to bottom. A positive ground water table means that the land is under water. ​ 
-The lowest soil layer that is not completely filled with water is defined as the "​official"​ ground water table layer and is the one being printed. Soil layers above this may have water in its effective porosity. ​ 
  
 === Soil moisture deficit === === Soil moisture deficit ===
Line 168: Line 164:
  
 ==== Infiltration ==== ==== Infiltration ====
-Infiltration is calculated from the sum of rain and snowmelt (//​infilt0//​) .+Infiltration is calculated from the sum of rain and snowmelt (//​infilt0//​, mm/time step) .
  
 <m> infilt0 = rainfall + melt </​m> ​ <m> infilt0 = rainfall + melt </​m> ​
Line 195: Line 191:
 The actual infiltration is calculated by subtracting the macropore flow and surface runoff from the sum of snow melt and rain. The actual infiltration is calculated by subtracting the macropore flow and surface runoff from the sum of snow melt and rain.
  
-<m> infilt = infilt0 - macroflow ​– infoverflow </m>+<m> infilt = infilt0 - macroflow ​infoverflow </m>
  
 === Additional infiltration limitation by frozen soil === === Additional infiltration limitation by frozen soil ===
Line 268: Line 264:
 ==== Saturated surface runoff ==== ==== Saturated surface runoff ====
    
-Surface runoff due to a high ground water table (//satoverflow//) occurs when the water table in the upper soil layer reaches above the surface. It depends on a parameter //srrcs// which is dependent on land use. The recession parameter is corrected with the correction factor //​rrcscorr//​ for different parameter regions (parreg). It is defined as an increase.+Surface runoff due to a high ground water table (//q//, mm/time step) occurs when the water table in the upper soil layer reaches above the surface. It depends on a parameter //srrcs// which is dependent on land use. The recession parameter is corrected with the correction factor //​rrcscorr//​ for different parameter regions (parreg). It is defined as an increase.
  
 <m> srrcs = srrcs*(1+rrcscorr) </m> <m> srrcs = srrcs*(1+rrcscorr) </m>
  
-<​m> ​satoverflow ​= MAX(srrcs * (soil(1)-wp(1)-fc(1)-ep(1)),​0.) </m>+<​m> ​= MAX(srrcs * (soil(1)-wp(1)-fc(1)-ep(1)),​0.) </m>
  
 Runoff is removed from the uppermost soil layer. The total surface runoff (due to high ground water table and low infiltration capacity) is calculated and printed. Runoff is removed from the uppermost soil layer. The total surface runoff (due to high ground water table and low infiltration capacity) is calculated and printed.
Line 368: Line 364:
 ==== Soil temperature and snow depth ==== ==== Soil temperature and snow depth ====
  
-Soil layer temperature (//​soiltemp//​) is calculated as a balance of three temperatures;​ previous time step soil layer temperature,​ soil temperature at deep depth (//​deeptemp//​) and air temperature (//temp//). The weight of the deep soil is constant (0.001), while the weight of the air temperature (//​weightair//​) depends on snow depth (//​snowdepth//​) and parameters. The soil memory (//​soilmem//​) depends on depth and land use, with parameters //surfmem// and //​depthrel//​. The memory of deep soil temperature is a general parameter (//​deepmem//​).+Soil layer temperature (//​soiltemp//​) is calculated as a balance of three temperatures;​ previous time step soil layer temperature,​ soil temperature at deep depth (//​deeptemp//​) and air temperature (//T//). The weight of the deep soil is constant (0.001), while the weight of the air temperature (//​weightair//​) depends on snow depth (//​snowdepth//​) and parameters. The soil memory (//​soilmem//​) depends on depth and land use, with parameters //surfmem// and //​depthrel//​. The memory of deep soil temperature is a general parameter (//​deepmem//​).
  
 <m> soilmem = {lbrace}{ <m> soilmem = {lbrace}{
Line 377: Line 373:
 <m> weigth_{air}={1}/​{soilmem+10*snowdepth} </m> <m> weigth_{air}={1}/​{soilmem+10*snowdepth} </m>
  
-<m> deeptemp=weight_{air}*temp+(1-weight_{air})*deeptemp </m>+<m> deeptemp=weight_{air}*T+(1-weight_{air})*deeptemp </m>
  
-<m> soiltemp=weight_{air}*temp+(1-weight_{air}-weight_{deep} )*soiltemp+weight_{deep}*deeptemp </m>+<m> soiltemp=weight_{air}*T+(1-weight_{air}-weight_{deep} )*soiltemp+weight_{deep}*deeptemp </m>
  
 In the default snow depth model, snow density (//​snowdens//​) depends on the snow's age in days (//​snowage//​). Snow density for fresh snow (//​sdnsnew//​) and the increase of density with snow age (//​snowdensdt//​) are general parameters (~ 0.1 and ~0.002). The snow's age increases by one every time step, but are weighted with age (0) for any new snow. In the default snow depth model, snow density (//​snowdens//​) depends on the snow's age in days (//​snowage//​). Snow density for fresh snow (//​sdnsnew//​) and the increase of density with snow age (//​snowdensdt//​) are general parameters (~ 0.1 and ~0.002). The snow's age increases by one every time step, but are weighted with age (0) for any new snow.
Line 403: Line 399:
 ^Section ^Symbol ^Parameter/​Data ^File ^ ^Section ^Symbol ^Parameter/​Data ^File ^
 |Snow melt|//​cmlt,​ ttmp//​|//​cmlt,​ ttmp//​|[[start:​hype_file_reference:​par.txt|par.txt]]| |Snow melt|//​cmlt,​ ttmp//​|//​cmlt,​ ttmp//​|[[start:​hype_file_reference:​par.txt|par.txt]]|
-|:::|//temp//​|calculated from|[[start:​hype_file_reference:​tobs.txt|Tobs.txt]]|+|:::|//T//​|calculated from|[[start:​hype_file_reference:​tobs.txt|Tobs.txt]]|
 |Snow cover|//​stdelev//​|//​elev_std//​|[[start:​hype_file_reference:​geodata.txt|GeoData.txt]]| |Snow cover|//​stdelev//​|//​elev_std//​|[[start:​hype_file_reference:​geodata.txt|GeoData.txt]]|
 |:::​|//​fscmax,​ fscdist0, fscdist1, fsck1, fsckexp//​|//​fscmax,​ fscdist0, fscdist1, fsck1, fsckexp//​|[[start:​hype_file_reference:​par.txt|par.txt]]| |:::​|//​fscmax,​ fscdist0, fscdist1, fsck1, fsckexp//​|//​fscmax,​ fscdist0, fscdist1, fsck1, fsckexp//​|[[start:​hype_file_reference:​par.txt|par.txt]]|
Line 431: Line 427:
 <m> glac_{vol} = coef*{glac_{area}}^{exp} </m> <m> glac_{vol} = coef*{glac_{area}}^{exp} </m>
  
-<m> coef = glacvcoef*{e}^{c} </m>+<m> coef = coef0*{e}^{c} </m>
  
 The initial glacier volume if calculated from class area: The initial glacier volume if calculated from class area:
Line 441: Line 437:
 <m> glac_{area} = ({glac_{vol}*{1/​{coef}}})^{1/​exp} </m> <m> glac_{area} = ({glac_{vol}*{1/​{coef}}})^{1/​exp} </m>
  
-The equation coefficients //coef// and //exp// can have different values for specific glaciers. The first coefficient coef is calculated as the product of //EXP( c)//, where //c// is a glacier volume correction, and a general parameter (//​glacvcoef/​glacvcoef1//​) depending on glacier type. The second coefficient,​ //exp// is a general parameter (//​glacvexp/​glacvexp1//​) depending on glacier type. Glacier density (//​glacdens//​) is a general model parameter (m3 water / m3 ice).+The equation coefficients //coef// and //exp// can have different values for specific glaciers. The first coefficient coef is calculated as the product of //EXP( c)//, where //c// is a glacier volume correction, and a general parameter (//​glacvcoef/​glacvcoef1//​) depending on glacier type. The second coefficient,​ //exp// is a general parameter (//​glacvexp/​glacvexp1//​) depending on glacier type. Glacier density (//​glacdens//​) is a general model parameter (m<​sup>​3</​sup> ​water / m<​sup>​3</​sup> ​ice).
  
 Glaciers are divided into four types. The default type is mountain glacier, the alternatives are ice cap, ice sheet and infinite glacier. Glacier type is given as input, or determined by the glacier area (a threshold (//​glac2arlim//,​ a general parameter). The glacier area is used to determine the glacier type if it is not given as input and the threshold parameter is set. The glaciers will then be divided into mountain glaciers and ice caps.  Glaciers are divided into four types. The default type is mountain glacier, the alternatives are ice cap, ice sheet and infinite glacier. Glacier type is given as input, or determined by the glacier area (a threshold (//​glac2arlim//,​ a general parameter). The glacier area is used to determine the glacier type if it is not given as input and the threshold parameter is set. The glaciers will then be divided into mountain glaciers and ice caps. 
Line 481: Line 477:
 |//​c//​|//​logvolcorr//​|[[start:​hype_file_reference:​glacierdata.txt|GlacierData.txt]]| |//​c//​|//​logvolcorr//​|[[start:​hype_file_reference:​glacierdata.txt|GlacierData.txt]]|
 | |//​glactype//​|:::​| | |//​glactype//​|:::​|
-|//T//|see needed data in [[start:​hype_model_description:​processes_above_ground#​links_to_file_reference|temperature]]| |+|//T//|see needed data in [[start:​hype_model_description:​processes_above_ground#​links_to_file_reference|Links for temperature]]| |
 |//​swrad//​|calculated or from|[[start:​hype_file_reference:​swobs.txt|SWobs.txt]]| |//​swrad//​|calculated or from|[[start:​hype_file_reference:​swobs.txt|SWobs.txt]]|
-|//glacvcoef//​|//​glacvcoef//​ or //​glacvcoef1//​|[[start:​hype_file_reference:​par.txt|par.txt]]| +|//coef0//​|//​glacvcoef//​ or //​glacvcoef1//​|[[start:​hype_file_reference:​par.txt|par.txt]]| 
-|//​exp//​|//​glacvexpglacvexp1//​|:::​|+|//​exp//​|//​glacvexp// or //glacvexp1//​|:::​|
 |//cmlt, ttmp//​|//​glaccmlt,​ glacttmp//​|:::​| |//cmlt, ttmp//​|//​glaccmlt,​ glacttmp//​|:::​|
 |//​albedo<​sub>​snow</​sub>//​|calculated from //snalbmin, snalbmax, snalbkexp//​|:::​| |//​albedo<​sub>​snow</​sub>//​|calculated from //snalbmin, snalbmax, snalbkexp//​|:::​|
start/hype_model_description/hype_land.1539605054.txt.gz · Last modified: 2018/10/15 14:04 by cpers