User Tools

Site Tools


start:hype_model_description:hype_human_water

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
Next revision Both sides next revision
start:hype_model_description:hype_human_water [2018/10/12 13:59]
cpers [Irrigation water demand]
start:hype_model_description:hype_human_water [2018/11/15 09:39]
cpers
Line 24: Line 24:
  
 ^Parameter/​Data ^File ^ ^Parameter/​Data ^File ^
-|//​lrwet_area,​ lrwet_dep, lrwet_part//​|[[start:​hype_file_reference:​geodata.txt|GeoData.txt]]| +|//​lrwet_area,​ lrwet_dep, lrwet_part, mrwet_area, mrwet_dep, mrwet_part//​|[[start:​hype_file_reference:​geodata.txt|GeoData.txt]]| 
-|//​mrwet_area,​ mrwet_dep, mrwet_part//​|:::​|+
  
  
Line 83: Line 83:
 (3) Up to a defined fraction of <​m>​S_{SW}</​m>​ (<​m>​P_{iwdfrac}</​m>,​ //iwdfrac// in par.txt): ​ (3) Up to a defined fraction of <​m>​S_{SW}</​m>​ (<​m>​P_{iwdfrac}</​m>,​ //iwdfrac// in par.txt): ​
  
-<​m>​W_{I,​D,​j}=min[(S_{SW}*P_{SSWCORR}*AWC-H)×P_{iwdfrac},​(AWC-H)]</​m>​+<​m>​W_{I,​D,​j}=min[(S_{SW}*P_{SSWCORR}*AWC-H)*P_{iwdfrac},​(AWC-H)]</​m>​
  
 The fraction can be larger than 1. For example, to irrigate to a level 10% above <​m>​S_{SW}</​m>,​ <​m>​P_{iwdfrac}</​m>​ =1.1. <​m>​W_{I,​D,​j}</​m>​ is, however, limited to <​m>​AWC</​m>​. The fraction can be larger than 1. For example, to irrigate to a level 10% above <​m>​S_{SW}</​m>,​ <​m>​P_{iwdfrac}</​m>​ =1.1. <​m>​W_{I,​D,​j}</​m>​ is, however, limited to <​m>​AWC</​m>​.
Line 128: Line 128:
 If the <​m>​W_{L,​D,​s}</​m>​ volume is available at the source, the demanded water is withdrawn. Otherwise, only the available volume (<​m>​V_s</​m>​) is withdrawn. The withdrawal can also be scaled with the user-defined parameter <​m>​P_{I,​S}</​m>​ (//pirrs// in par.txt): ​ If the <​m>​W_{L,​D,​s}</​m>​ volume is available at the source, the demanded water is withdrawn. Otherwise, only the available volume (<​m>​V_s</​m>​) is withdrawn. The withdrawal can also be scaled with the user-defined parameter <​m>​P_{I,​S}</​m>​ (//pirrs// in par.txt): ​
  
-<m> W_{L,​A,​s(1)}=min(W_{L,​D,​s},​V_{s1})×P_{I,S} </m>+<m> W_{L,​A,​s(1)}=min(W_{L,​D,​s},​V_{s1})*P_{I,S} </m>
  
 <m> W_{L,​D,​s2}=W_{L,​D,​s}-{W_{L,​A,​s(1)}}/​{P_{I,​S}} </​m> ​ <m> W_{L,​D,​s2}=W_{L,​D,​s}-{W_{L,​A,​s(1)}}/​{P_{I,​S}} </​m> ​
Line 144: Line 144:
 If any groundwater demand exists, the model withdraws a user-specified fraction of <​m>​W_{L,​D,​g}</​m>​ from an unlimited source outside the model domain or if an aquifer is simulated connected to the subbasin from this aquifer. The first case conceptually represents a large deep aquifer source, which is currently outside of the model domain. <​m>​W_{L,​a,​g}</​m>​ is the abstracted groundwater and <​m>​P_{I,​G}</​m>​ (//pirrg// in par.txt) is the groundwater withdrawal fraction: If any groundwater demand exists, the model withdraws a user-specified fraction of <​m>​W_{L,​D,​g}</​m>​ from an unlimited source outside the model domain or if an aquifer is simulated connected to the subbasin from this aquifer. The first case conceptually represents a large deep aquifer source, which is currently outside of the model domain. <​m>​W_{L,​a,​g}</​m>​ is the abstracted groundwater and <​m>​P_{I,​G}</​m>​ (//pirrg// in par.txt) is the groundwater withdrawal fraction:
  
-<m> W_{L,​A,​g}=W_{L,​D,​g}×P_{I,G} </m>+<m> W_{L,​A,​g}=W_{L,​D,​g}*P_{I,G} </m>
  
 To simulate a more dynamic conjunctive use of groundwater and surface water sources, the model allows for compensation of remaining surface water demands from the groundwater source. This compensation is only allowed if both groundwater and surface water sources are used (//​0<​gw_part<​1//​),​ and if the //irrcomp// parameter is >0. The //irrcomp// parameter defines the degree of compensation allowed, i.e. the fraction of the residual surface water demand which can be met through source compensation. The compensation algorithm is as follows: if any surface water demand remains (<​m>​W_{L,​D,​s,​1}</​m>​ > 0) and the groundwater is not depleted, the groundwater withdrawal cycle is calculated once more using the scaled residual surface water demand. Finally, after possible source compensation,​ the remaining (surface) water demand at the sub-basin scale (<​m>​W_{L,​D,​1}</​m>​) is calculated. To simulate a more dynamic conjunctive use of groundwater and surface water sources, the model allows for compensation of remaining surface water demands from the groundwater source. This compensation is only allowed if both groundwater and surface water sources are used (//​0<​gw_part<​1//​),​ and if the //irrcomp// parameter is >0. The //irrcomp// parameter defines the degree of compensation allowed, i.e. the fraction of the residual surface water demand which can be met through source compensation. The compensation algorithm is as follows: if any surface water demand remains (<​m>​W_{L,​D,​s,​1}</​m>​ > 0) and the groundwater is not depleted, the groundwater withdrawal cycle is calculated once more using the scaled residual surface water demand. Finally, after possible source compensation,​ the remaining (surface) water demand at the sub-basin scale (<​m>​W_{L,​D,​1}</​m>​) is calculated.
Line 159: Line 159:
 The total water demand from the regional source (<​m>​W_{R,​D}</​m>​) is then calculated as the sum of the demand from each connected sub-basin, scaled by a parameter controlling the strength of the regional connection (<​m>​P_regirr</​m>,​ //regirr// in par.txt): The total water demand from the regional source (<​m>​W_{R,​D}</​m>​) is then calculated as the sum of the demand from each connected sub-basin, scaled by a parameter controlling the strength of the regional connection (<​m>​P_regirr</​m>,​ //regirr// in par.txt):
  
-<m> W_{R,D}= (sum{i=1}{N}{W_{R,​D,​i}} )×P_{regirr} </m>+<m> W_{R,D}= (sum{i=1}{N}{W_{R,​D,​i}} )*P_{regirr} </m>
  
 The regional demand can be met from two sources in sub-basin <​m>​D_R</​m>:​ the olake and the main river. If the regional source sub-basin has an olake, and if the //irrdam// input variable is set to 1 for that sub-basin, the model attempts to withdraw <​m>​W_{R,​D}</​m>​ first from the olake and then the residual from the main river. If not, the model only attempts to withdraw <​m>​W_{R,​D}</​m>​ from the main river. The regional abstraction (<​m>​W_{R,​a}</​m>​) is limited by the volume available at the source (<​m>​V_r</​m>​) and the scaling parameter <​m>​P_{I,​S}</​m>:​ The regional demand can be met from two sources in sub-basin <​m>​D_R</​m>:​ the olake and the main river. If the regional source sub-basin has an olake, and if the //irrdam// input variable is set to 1 for that sub-basin, the model attempts to withdraw <​m>​W_{R,​D}</​m>​ first from the olake and then the residual from the main river. If not, the model only attempts to withdraw <​m>​W_{R,​D}</​m>​ from the main river. The regional abstraction (<​m>​W_{R,​a}</​m>​) is limited by the volume available at the source (<​m>​V_r</​m>​) and the scaling parameter <​m>​P_{I,​S}</​m>:​
Line 167: Line 167:
 <m> W_{R,​D2}=W_{R,​D}-{W_{R,​A1}}/​{P_{I,​S}} </​m> ​ <m> W_{R,​D2}=W_{R,​D}-{W_{R,​A1}}/​{P_{I,​S}} </​m> ​
  
-<m> W_{R,​A2}=min(W_{R,​D2},​V_{r2} )×P_{I,S} </m>+<m> W_{R,​A2}=min(W_{R,​D2},​V_{r2} )*P_{I,S} </m>
  
 <m> W_{R,​A}=W_{R,​A1}+ W_{R,A2} </m> <m> W_{R,​A}=W_{R,​A1}+ W_{R,A2} </m>
Line 177: Line 177:
 The concentrations of the withdrawn water are the same as that of the irrigation water source. If the water originates from several sources, the volume-weighted concentration is calculated. If desired, the model can simulate sedimentation tanks, in which a defined fraction of the particulate phosphorous (//pp//) and organic nitrogen (//on//) settles: The concentrations of the withdrawn water are the same as that of the irrigation water source. If the water originates from several sources, the volume-weighted concentration is calculated. If desired, the model can simulate sedimentation tanks, in which a defined fraction of the particulate phosphorous (//pp//) and organic nitrogen (//on//) settles:
  
-<m> C_{a,​pp}=C_{src,​pp}×(1-P_{cirrsink} ) </m>+<m> C_{a,​pp}=C_{src,​pp}*(1-P_{cirrsink} ) </m>
  
-<m> C_{a,​on}=C_{src,​on}×(1-P_{cirrsink} ) </m>+<m> C_{a,​on}=C_{src,​on}*(1-P_{cirrsink} ) </m>
  
 where <​m>​C_{,​a}</​m>​ is the concentration of the abstracted water after settling, <​m>​C_{src}</​m>​ is the concentration of the source, and <​m>​P_{cirrsink}</​m>​ is the concentration reduction fraction (//​cirrsink//​ parameter in par.txt). To use sedimentation tanks in a region, the concentration reduction fraction needs to be set in par.txt (//​0<​cirrsink≤1//​). where <​m>​C_{,​a}</​m>​ is the concentration of the abstracted water after settling, <​m>​C_{src}</​m>​ is the concentration of the source, and <​m>​P_{cirrsink}</​m>​ is the concentration reduction fraction (//​cirrsink//​ parameter in par.txt). To use sedimentation tanks in a region, the concentration reduction fraction needs to be set in par.txt (//​0<​cirrsink≤1//​).
Line 189: Line 189:
 The regionally abstracted water (<​m>​W_{R,​A}</​m>​) is first distributed to each connected sub-basin (i) according to their proportional demand, and then scaled to the local scale using the respective regional efficiency: The regionally abstracted water (<​m>​W_{R,​A}</​m>​) is first distributed to each connected sub-basin (i) according to their proportional demand, and then scaled to the local scale using the respective regional efficiency:
  
-<m> W_{R,​A,​i}=W_{R,​A}× {W_{R,D,i}×P_{regirr}}/​{W_{R,​D}}×E_{R,i} </m>+<m> W_{R,​A,​i}=W_{R,​A}{W_{R,D,i}*P_{regirr}}/​{W_{R,​D}}*E_{R,i} </m>
  
 For a given sub-basin, the total amount of abstracted water available at the local scale (<​m>​W_{L,​A,​i,​tot}</​m>​) is calculated and then scaled, using the local efficiency, to represent the water applied to the soil (<​m>​W_{I,​A,​i}</​m>​):​ For a given sub-basin, the total amount of abstracted water available at the local scale (<​m>​W_{L,​A,​i,​tot}</​m>​) is calculated and then scaled, using the local efficiency, to represent the water applied to the soil (<​m>​W_{I,​A,​i}</​m>​):​
Line 195: Line 195:
 <m> W_{L,​A,​i,​tot}=W_{L,​A,​i}+W_{R,​A,​i} </m> <m> W_{L,​A,​i,​tot}=W_{L,​A,​i}+W_{R,​A,​i} </m>
  
-<m> W_{I,​A,​i}=W_{L,​A,​i,​tot}×E_{L,i} </m>+<m> W_{I,​A,​i}=W_{L,​A,​i,​tot}*E_{L,i} </m>
  
 <​m>​W_{I,​A,​i}</​m>​ is then distributed onto each irrigated class in proportion to its water demand: <​m>​W_{I,​A,​i}</​m>​ is then distributed onto each irrigated class in proportion to its water demand:
  
-<m> W_{I,​A,​j}=W_{I,​A,​i}×{W_{I,​D,​j}}/​{W_{I,​D}} </m>+<m> W_{I,​A,​j}=W_{I,​A,​i}*{W_{I,​D,​j}}/​{W_{I,​D}} </m>
  
 <​m>​W_{I,​A,​j}</​m>​ is added to the soil water of class j as additional infiltration. <​m>​W_{I,​A,​j}</​m>​ is divided between the top two soil layers according to the epotdist function, beginning with the second layer. For unlimited irrigation: <​m>​W_{I,​A,​j} =  W_{I,​D,​j}</​m>​ , <​m>​P_{I,​S} =1</​m>​ and <​m>​P_{I,​G} = 1</​m>​. <​m>​W_{I,​A,​j}</​m>​ is added to the soil water of class j as additional infiltration. <​m>​W_{I,​A,​j}</​m>​ is divided between the top two soil layers according to the epotdist function, beginning with the second layer. For unlimited irrigation: <​m>​W_{I,​A,​j} =  W_{I,​D,​j}</​m>​ , <​m>​P_{I,​S} =1</​m>​ and <​m>​P_{I,​G} = 1</​m>​.
Line 335: Line 335:
 |:::​|//​subid,​ receiver, flow//|:::| |:::​|//​subid,​ receiver, flow//|:::|
 |:::​|//​dwtr//​|[[start:​hype_file_reference:​xobs.txt|Xobs.txt]]| |:::​|//​dwtr//​|[[start:​hype_file_reference:​xobs.txt|Xobs.txt]]|
 +
 ==== Links to relevant procedures in the code ==== ==== Links to relevant procedures in the code ====
  
Line 348: Line 349:
 |:::​|water_transfer_from_outlet_lake| ::: | |:::​|water_transfer_from_outlet_lake| ::: |
  
 +===== Dams =====
  
 +Regulation of flow through dams is described in the Chapter about [[start:​hype_model_description:​hype_routing| Rivers and lakes]]. Dams of different purposes and regulation management can be simulated. See details in the sections on [[start:​hype_model_description:​hype_routing#​simple_outlet_lake_or_dam_olake| Simple outlet lake or dam]] and the special case of an [[start:​hype_model_description:​hype_routing#​outlet_lake_with_two_outlets| Outlet lake with two outlets]].
start/hype_model_description/hype_human_water.txt · Last modified: 2024/02/21 09:14 by cpers