User Tools

Site Tools


start:hype_model_description:hype_np_riv_lake

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_np_riv_lake [2018/09/10 12:06]
cpers [Nitrogen and phosphorus processes in rivers and lakes]
start:hype_model_description:hype_np_riv_lake [2018/11/16 15:46] (current)
cpers [Basic assumptions]
Line 1: Line 1:
 ====== Nitrogen and phosphorus processes in rivers and lakes ====== ====== Nitrogen and phosphorus processes in rivers and lakes ======
  
 +HYPE simulates concentration of inorganic nitrogen (IN), organic nitrogen (ON), soluble reactive phosphorus (SP) and particulate phosphorus (PP) in discharge and other surface waters. As output also the total nitrogen (TN) and phosphorus (TP) concentration can be had. In addition total suspended sediments (TS) can be simulated. It consists of the sum of suspended sediments (SS) and algae (AE) simulated concentrations.
 +
 +===== Basic assumptions =====
  
 Transformations of nutrients take place in lakes and rivers. For lakes, which are divided into fast (FLP) and slow (SLP) lake parts, the process is performed only in SLP (Fig. 1). For rivers, which hold delayed water in a queue and in the damping box, the processes is performed only in the damping box.  Transformations of nutrients take place in lakes and rivers. For lakes, which are divided into fast (FLP) and slow (SLP) lake parts, the process is performed only in SLP (Fig. 1). For rivers, which hold delayed water in a queue and in the damping box, the processes is performed only in the damping box. 
Line 9: Line 12:
 The processes of denitrification,​ primary production and mineralization have been implemented for both rivers and lakes. For particulate phosphorus (PP) there is an exchange with the river sediments. The rivers dimensions are used in the calculation of these processes. The width and depth of the watercourse are calculated from a number of empirical equations (for more information on these equations see "​Modelling phosphorus transport and retention in river networks"​ by Jörgen Rosberg). The processes of denitrification,​ primary production and mineralization have been implemented for both rivers and lakes. For particulate phosphorus (PP) there is an exchange with the river sediments. The rivers dimensions are used in the calculation of these processes. The width and depth of the watercourse are calculated from a number of empirical equations (for more information on these equations see "​Modelling phosphorus transport and retention in river networks"​ by Jörgen Rosberg).
  
-<m> velocity = 10^{rivvel1} * meanflow^{rivvel2} * (flow/​meanflow)^{rivvel3} </m>+<m> velocity = 10^{vel1} * meanflow^{vel2} * (flow/​meanflow)^{vel3} </m>
  
-<m> width = 10^{rivwidth1} * (flow/​velocity)^{rivwidth2+rivwidth3*log10(flow/​velocity)} </m>+<m> width = 10^{width1} * (flow/​velocity)^{width2+width3*log10(flow/​velocity)} </m>
  
 <m> depth = (flow / velocity) / width </m> <m> depth = (flow / velocity) / width </m>
                        
-where //rivvel1, rivvel2, rivvel3, rivwidth1, rivwidth2// and //rivwidth3// are lake region dependent parameters ​located in the file par.txt, //flow// is the flow of water in the watercourse (//m3/s//) and //​meanflow//​ is a 365-day rolling average flow. The river’s width is limited by the dead volume width and a parameter //​maxwidth//​. The rivers length is estimated at square root of catchment area or is given in GeoData.txt for the main river. The watercourses **surface ​area** is calculated as the length times the width, where the maximum of the above-calculated width and dead volume width is used.+where //vel1-vel3// and //width1-width3// are lake region dependent parameters, //flow// is the flow of water in the watercourse (//m3/s//) and //​meanflow//​ is a 365-day rolling average flow. The river’s width is limited by the dead volume width and a parameter //​maxwidth//​. The rivers length is estimated at square root of catchment area or is given as input. The watercourses **bottom ​area** is calculated as the length times the width, where the maximum of the above-calculated width and dead volume width is used. The bottom area is used for the nutrient processes. If the river is a class and thus has a defined class area, the class area is used to estimate the bottom area instead.
  
-The water temperature (//​watertemp//​) is used in some of the process calculations. It is calculated through weighting the air temperature (//​airtemp//​) and yesterday'​s water temperature. The weighting constant (similar to moving average period) //​watertemp_days//​ is set to 20 days for rivers and is by default 5 days for lakes. For lakes a depth dependent //​watertemp_days//​ can also be used. Then the weighting constant is set to the equivialent days of the lake's depths up to maximum the general parameter //​laketemp//​ plus 5 days. Note that the water temperature can fall below <​m>​0^{circ}C</​m>​. 
  
-  IF(laketemp>​0) ​watertemp_days ​= MIN(MAX(lake_depth,​5),​5+laketemp)  +The water temperature ​(//​T<​sub>​w</​sub>//​) is used in some of the process calculations. It is calculated through weighting the air temperature (//T//) and yesterday'​s water temperature. The weighting constant (similar to moving average period) //​T<​sub>​wdays</​sub>//​ is set to 20 days for rivers and is by default 5 days for lakes. For lakes a depth dependent //​T<​sub>​wdays</​sub>//​ can also be used. Then the weighting constant is set to the equivialent days of the lake's depths up to a maximum (//​maxpar//,​ a parameter) plus 5 days. Note that the water temperature can fall below <m>0^{circ}C</​m>​. 
-  ​weightair = 1 watertemp_days ​ + 
-  ​watertemp = (1-weightair) * watertemp + weightair * airtemp+<m> T_w = (1-{1/​{T_wdays}})*T_w + {1/​{T_wdays}}*T </​m>​ 
 + 
 +<m> T_wdays ​delim{lbrace}{  
 +   ​matrix{2}{2}{ 
 +    {MIN(MAX(lake_depth,​5),​5+maxpar)} {lake} ​ 
 +    20 river 
 +   }}{} </m> 
   ​   ​
 Alternatively the [[start:​hype_model_description:​hype_tracer#​water_temperature_t2|water temperature calculated as a tracer]] (called T2) may be used. This model option is set in info.txt. Alternatively the [[start:​hype_model_description:​hype_tracer#​water_temperature_t2|water temperature calculated as a tracer]] (called T2) may be used. This model option is set in info.txt.
 +
 +
 +
 +==== Links to file reference ====
 +
 +^Symbol ^Parameter/​Data ^File ^
 +| |//rivlen, loc_rivlen//​|[[start:​hype_file_reference:​geodata.txt|GeoData.txt]]|
 +|//vel1, vel2, vel3//​|//​rivvel1,​ rivvel2, rivvel3//​|[[start:​hype_file_reference:​par.txt|par.txt]]|
 +|//width1, width2, width3//​|//​rivwidth1,​ rivwidth2, rivwidth3//​|:::​|
 +| |//deadl, deadm, maxwidth//​|:::​|
 +|//​maxpar//​|//​laketemp//​|:::​|
 +|//​lake<​sub>​depth</​sub>//​|//​lake_depth//​|[[start:​hype_file_reference:​geodata.txt|GeoData.txt]] or [[start:​hype_file_reference:​lakedata.txt|LakeData.txt]] or [[start:​hype_file_reference:​damdata.txt|DamData.txt]]|
 +|:::​|//​gldepo,​ olldepth, gldepi, illdepth//​|[[start:​hype_file_reference:​par.txt|par.txt]]|
 +
 +
 ==== Links to relevant procedures in the code ====  ​ ==== Links to relevant procedures in the code ====  ​
  
Line 35: Line 59:
 ===== Denitrification ===== ===== Denitrification =====
  
-Denitrification,​ a sink for inorganic nitrogen in lakes and rivers, is a function of the bottom area, the IN concentration (//concfcn//) in water volume, water temperature (//tmpfcn//) and the general parameters //​denitwrl//​ (local river), //​denitwrm//​ (main river) and //denitwl// (lakes). In the concentration function, the half saturation parameter (//hsatINw//) is a general parameter, but it was in earlier HYPE versions a constant equal to 1.5 mg/L. Denitrification (//denitri_water//, //kg/day//) is limited to a maximum of 50% (//​maxdenitriwater//​) ​of the available IN pool (i.e. in SLP).+Denitrification,​ a sink for inorganic nitrogen in lakes and rivers, is a function of the bottom area, the IN concentration (//conc//) in water volume, water temperature (<​m>​T_w<​/m>) and a rate parameter. The dependence on concentration is formulated as a half saturation equation. In the concentration function, the half saturation parameter (//par<​sub>​half</​sub>​//) is a general parameter, but it was in earlier HYPE versions a constant equal to 1.5 mg/L. Denitrification (//denitr//, //kg/day//) is limited to a maximum of 50% of the available IN pool (i.e. in SLP).
  
-  ​tmpfcn = 2**((watertemp ​- 20/ 10+<​m> ​tmpfcn=delim{lbrace}{  
-  IF(watertemp ​< 5.) tmpfcn = tmpfcn * (watertemp / 5) +matrix{3}{2}{ 
-  ​IF(watertemp ​0.) tmpfcn = 0. +    0 {T_w<​0} 
-  concfcn = conc / (conc + hsatINw) +    {T_w/5}*{2^{{T_w-20}/{10}}} {0<T_w<5} 
-  ​denitri_water ​= MIN(maxdenitriwater ​* INpool, ​denitw ​* area * concfcn * tmpfcn)+    ​2^{{T_w-20}/​{10}} ​ else 
 +   ​}}{} ​</m> 
 + 
 +<​m> ​concfcn= conc/{conc+par_half} </m> 
 + 
 +<m> denitr ​= MIN(0.5* INpool,par_den*area*concfcn*tmpfcn) ​</​m>​ 
 + 
 +The bottom area (//area//) is equal to the lake area for lakes, and watercourse width multiplied by the length for rivers (see above). The rate parameter (//​par<​sub>​den</​sub>//​) depend on the water body. Local and main rivers have general parameters (//​denitwrl,​ denitwrm//​),​ while lakes in addition to use a general parameter (//​denitwl//​) can have lake specific values. 
 + 
 +==== Links to file reference ==== 
 + 
 +^Symbol ^Parameter/​Data ^File ^ 
 +|<​m>​par_half</​m>​|//​hsatINw//​|[[start:​hype_file_reference:​par.txt|par.txt]]| 
 +|<​m>​par_den</​m>​|//​denitwrl,​ denitwrm//​|:::​| 
 +|:::​|//​denitwl//​|[[start:​hype_file_reference:​par.txt|par.txt]] or [[start:​hype_file_reference:​lakedata.txt|LakeData.txt]]|
  
-The bottom area (//area//) is equal to the lake area, for lakes, and watercourse width multiplied by the length, for rivers. 
  
 ==== Links to relevant procedures in the code ====  ​ ==== Links to relevant procedures in the code ====  ​
Line 55: Line 92:
 ===== Primary production and mineralization ===== ===== Primary production and mineralization =====
  
-Primary production in lakes and rivers is a source of organic nitrogen and particulate phosphorus and a sink for inorganic nitrogen and soluble reactive phosphorus in the model. The reverse is true for mineralization. Primary production and mineralization ​is controlled by two temperature functions. The first (//​tmpfcn1//​) is solely dependent on the water temperature. It simulates the increased activity at warmer temperatures. The second (//​tmpfcn2//​) governs the relationship between primary production and mineralization and determines which one dominates. Net primary production is highest in spring (northern hemisphere) and changes into net mineralization when the temperature //T10// is less than the temperature //T20// in autumn. These two temperatures are calculated as the average water temperature of 10 and 20 days. +Primary production in lakes and rivers is a source of organic nitrogen and particulate phosphorus and a sink for inorganic nitrogen and soluble reactive phosphorus in the model. The reverse is true for mineralization. The processes is modelled together and only one is active at the time. Primary production and mineralization ​are controlled by two temperature functions. The first (//​tmpfcn1//​) is solely dependent on the water temperature ​(//​T<​sub>​w</​sub>//​). It simulates the increased activity at warmer temperatures. The second (//​tmpfcn2//​) governs the relationship between primary production and mineralization and determines which one dominates. Net primary production is highest in spring (northern hemisphere) and changes into net mineralization when the temperature //T10// is less than the temperature //T20// in autumn. These two temperatures are calculated as the average water temperature of 10 and 20 days. 
  
  
-<m> tmpfcn1 = watertemp ​/ 20. </m>+<m> tmpfcn1 = T_w / 20. </m>
  
 <m> tmpfcn2 = (T10 - T20) / 5. </m> <m> tmpfcn2 = (T10 - T20) / 5. </m>
Line 65: Line 102:
  
  
-The primary production and mineralization processes are also governed by long-term average concentrations of modelled total phosphorus ​in the lake (//TPconc//) or watercourse through a concentration function (//TPfcn//) with the general half saturation parameter //hsatTP//. The half saturation parameter was in earlier HYPE versions a constant equal to 0.05 mg/L. If phosphorus is not modelled a long-term average total phosphorus concentration as a lake region dependent parameter (//​tpmean//​) is used. If set, the long-term average concentration is reduced ​by the general parameter //​limsedPP// ​before using it in the concentration function. ​+The primary production and mineralization processes are also governed by long-term average concentrations of modelled total phosphorus (//TP//) in the lake or watercourse through a concentration function (//TPfcn//) with the general half saturation parameter //hsatTP//. The half saturation parameter was in earlier HYPE versions a constant equal to 0.05 mg/L. If phosphorus is not modelled a long-term average total phosphorus concentration as a lake region dependent parameter (//​tpmean//​) is used. A limiting concentration (//lim//) can be set, and then the long-term average concentration is reduced before using it in the concentration function. ​
  
-<m> TPfcn = (TPconc-limsedPP) ​/ (TPconc-limsedPP ​+ hsatTP</m>+<m> TPfcn = {TP-lim} {(TP-lim) + hsatTP</m>
  
-For lakes, the process is acting only in lake part SLP, while the processes are active throughout the watercourse volume. The estimated production/​mineralization (//​minprodNpot//,​ kg / day) is potential and may be limited by the availability of nutrients ​by only 50% of the available IN the pool (at the primary production) ​and ON-pool (for mineralization) can be transformed. The potential phosphorus conversion (//​minprodPpot//​) is calculated in the same way, but with its own parameter (//​wprodp//​) and a factor for phosphorus/​nitrogen ratio (//​NPratio//​ = 1/​7.2). ​Additionally, there is a restriction against SP and PP pools instead. The parameters //wprodn// and //wprodp// is generic or can be specified for each lake. The area is equal to lake area for lakes and river width multiplied by the length of the watercourse.+For lakes, the process is acting only in lake part SLP, while the processes are active throughout the watercourse volume. The estimated production/​mineralization (//​minprodNpot//,​ kg / day) is the potential ​transformation, ​and may be limited by the availability of nutrients. Only 50% of the available IN the pool (at the primary production) ​or 50% of the ON-pool (for mineralization) can be transformed. The potential phosphorus conversion (//​minprodPpot//​) is calculated in the same way, but with its own parameter (//​wprodp//​) and a factor for phosphorus/​nitrogen ratio (//​NPratio//​ = 1/​7.2). ​Similarly, there is a restriction against ​transforming maximum 50% of the SP and PP pools. The parameters //wprodn// and //wprodp// is generic or can be specified for each lake. The area is equal to lake area for lakes and bottom area for rivers (width multiplied by the length of the watercourse).
  
 <m> minprodNpot = wprodn * TPfcn * tmpfcn * area </m> <m> minprodNpot = wprodn * TPfcn * tmpfcn * area </m>
Line 75: Line 112:
 <m> minprodPpot = wprodp * TPfcn * tmpfcn * area * NPratio </m> <m> minprodPpot = wprodp * TPfcn * tmpfcn * area * NPratio </m>
   ​   ​
-For simulations of total suspended sediments, algae is simulated as a component. Nitrogen in algae is assumed to grow and decline with the same function as production and mineralisation of organic nitrogen. If nitrogen is simulated it uses the actual estimated production/​mineralisation,​ but otherwise the potential production/​mineralisation (//​minprodNpot//​) is used. The mineralisation of algae is limited to available amount, but the production is unlimited.+For simulations of total suspended sediments, algae is simulated as a component ​of the suspendend material. Nitrogen in algae is assumed to grow and decline with the same function as production and mineralisation of organic nitrogen. If nitrogen is simulated it uses the actual estimated production/​mineralisation,​ but otherwise the potential production/​mineralisation (//​minprodNpot//​) is used. The mineralisation of algae is limited to available amount, but the production is then unlimited. 
 + 
 +==== Links to file reference ==== 
 + 
 +^Symbol ^Parameter/​Data ^File ^ 
 +| |//hsatTP, tpmean//​|[[start:​hype_file_reference:​par.txt|par.txt]]| 
 +|//​lim//​|//​limsedPP//​|:::​| 
 +| |//wprodn, wprodp//​|[[start:​hype_file_reference:​par.txt|par.txt]] or [[start:​hype_file_reference:​lakedata.txt|LakeData.txt]]|
  
  
Line 90: Line 134:
 ===== Sedimentation/​Resuspension ===== ===== Sedimentation/​Resuspension =====
  
-Sedimentation in lakes is a sink for particulate phosphorus (PP) and organic nitrogen (ON), as well as suspended sediments (SS) and algae (AE). Sedimentation (//sedON, sedPP, sedSS, sedAE//, //m/day//) is calculated as a function of water concentration and lake area (//area//). The parameters //sedon//, //sedpp//, //​sedss//, ​and //sedae// are generic, but ON and PP sedimentation can be specified for each lake. The water concentration may be limited by the general parameters //​limsedON//,​ //​limsedPP//, ​and //​limsedSS//​.+Sedimentation in lakes is a sink for particulate phosphorus (PP) and organic nitrogen (ON), as well as for suspended sediments (SS) and algae (AE). Sedimentation (//sed//, //m/day//) is calculated as a function of concentration ​(//conc//) in the lake and lake area (//area//). The sedimentation rate (//​par<​sub>​sed</​sub>//​) is given by parameters ​(//sedon//, //sedpp//, //sedss//, //sedae//) which are generic, but ON and PP sedimentation can be specified for each lake. The concentration ​used in the equation ​may be limited ​(//​lim//​) ​by general parameters ​(//​limsedON//,​ //​limsedPP//,​ //​limsedSS//​), but not for AE (//lim//=0).
  
-<​m> ​sedON sedon * (waterconcON-limsedON) * area </m>+<​m> ​sed par_sed ​* (conc-lim) * area </m>
  
-<m> sedPP = sedpp * (waterconcPP-limsedPP) * area </m> 
  
-<m> sedSS = sedss * (waterconcSS-limsedSS) * area </m> 
  
-<m> sedAE = sedae * waterconcAE * area </​m>​ +In the river no particles ​are removed, but they are redistributed over time. This process is used for particulate phosphorus ​(PPand suspended sediments (SS). The process is a combination of sedimentation (//sed//, //m/day//) and resuspension (//resusp//, //m/day//). Particles in the sediments is collected in a pool (//sedimentpool//) which will increase as particles from the water volume (//waterpool//) settle ​at low water flows. The higher the water flow the lower the sedimentation ​and more particles returns to the water. This combined process is governed by the general parameter //sedexp//. The net effect is determined by the sign of the variable ​(-1 <//​sedresp//​ <1).
- +
-No particles (PP and SS) disappears in river but only redistributed over time through ​sedimentation (//sedPP, sedSS//, //m/day//) and resuspension (//resuspPP, resuspSS//, //m/day//). Particles in the sediments is simulated as a pool (//PPsedimentpool,​ SSsedimentpool//) which may be given particles from the water volume (//PPpool, SSpool//) at low water flows. The higher the water flow the less sedimentation ​while more particles returns to the water (//​resuspPP,​ resuspSS//). This combined process is governed by the general parameter //sedexp//. The net effect is determined by the sign of the flow (-1 <//​sedresp//​ <1).+
  
 <m> sedresp=max(-1.,​min(1.,​{{qbank-flow}/​qbank}^{sedexp}-{flow/​qbank}^{sedexp})) </m> <m> sedresp=max(-1.,​min(1.,​{{qbank-flow}/​qbank}^{sedexp}-{flow/​qbank}^{sedexp})) </m>
  
-<​m> ​sedPP = {lbrace}{matrix{2}{2}{{sedresp * PPpool}{sedresp > 0}{0}{sedresp < 0}}} </m>+<​m> ​sed = {lbrace}{matrix{2}{2}{{sedresp * waterpool}{sedresp > 0}{0}{sedresp < 0}}} </m>
  
-<​m> ​resuspPP ​= {lbrace}{matrix{2}{2}{{- sedresp * PPsedimentpool}{sedresp < 0}{0}{sedresp > 0}}} </m>+<​m> ​resusp ​= {lbrace}{matrix{2}{2}{{- sedresp * sedimentpool}{sedresp < 0}{0}{sedresp > 0}}} </m>
  
-<m> sedSS = {lbrace}{matrix{2}{2}{{sedresp * SSpool}{sedresp > 0}{0}{sedresp < 0}}} </m> 
  
-<m> resuspSS = {lbrace}{matrix{2}{2}{{- sedresp * SSsedimentpool}{sedresp < 0}{0}{sedresp ​> 0}}} </m>+where //flow// is the current river flow (//m3/s//) and //qbank// is the flow when river is filled to the brim. This flow is calculated as the second largest simulated flow in the last year. It is adjusted with a correction factor of 0.7 before use in the //sedresp// equation.
  
  
-where //flow// is the current river flow (//m3/s//) and //qbank// is the flow when river is filled to the brimThis flow is calculated as the second largest simulated flow in the last yearIt is adjusted with a correction factor of 0.7 before use in the //sedresp// equation.+ 
 +==== Links to file reference ==== 
 + 
 +^Symbol ^Parameter/Data ^File ^ 
 +|//area//|//slc_nn, area//|[[start:​hype_file_reference:​geodata.txt|GeoData.txt]]| 
 +|<​m>​par_sed<​/m>|//sedpp, sedon//​|[[start:​hype_file_reference:​par.txt|par.txt]] or [[start:​hype_file_reference:​lakedata.txt|LakeData.txt]]| 
 +|:::|//sedss, sedae//|[[start:​hype_file_reference:​par.txt|par.txt]]| 
 +|//​lim//​|//​limsedON,​ limsedPP, limsedSS//​|:::​| 
 +| |//​sedexp//​|:::​|
  
 ==== Links to relevant procedures in the code ====  ​ ==== Links to relevant procedures in the code ====  ​
Line 128: Line 175:
  
  
-Lakes can leak phosphorus from the bottom e.g. the release of phosphorus during oxygen deficiency or the mixing of the previously sedimented emissions. With two parameters the release of PP or SP is obtained for a specific lake. The parameters, ​//prodPP// and //prodSP// (//m/d//), is specified in LakeData.txt. The load (//loadPP, loadSP, kg/day//) is calculated as a function of lake water temperature (//T20//), average concentration of phosphorus in the lake (//TPconc//) and lake area (//area//). //TPconc// is currently ​hardcoded ​to 0.1 mg/l as an estimated avarage value.+Lakes can leak phosphorus from the bottom e.g. release of phosphorus during oxygen deficiency or the mixing of the previously sedimented emissions. Internal load can be simulated for SP or PP or both fractions of phosphorus. With two parameters the release of PP or SP is obtained for a specific lake (//par<​sub>​prod<​/sub>////m/d//). The load to the lake (//load, kg/day//) is calculated as a function of lake water temperature (//T<​sub>​w</​sub>​//), average concentration of phosphorus in the lake (//TPfcn//) and lake area (//area//). //TPfcn// is currently ​not depending on simulated concentration as it is set to constant ​0.1 mg/L. 
 + 
 +<m> load = {par_prod} * TPfcn * tmpfcn * area / 1000 </​m>​ 
 + 
 +<m> tmpfcn = 0.86^{delim{|}{T_w-15}{|}} ​ </​m>​ 
 + 
 +<m> TPfcn = 0.1 </​m>​ 
  
-<m> loadPP ​prodPP * TPconc * tmpfcn * area / 1000. </m>+==== Links to file reference ====
  
-<​m> ​loadSP = prodSP * TPconc * tmpfcn * area / 1000. </m>+^Symbol ^Parameter/​Data ^File ^ 
 +|<m>par_prod</m>|//prodPP, prodSP//​|[[start:​hype_file_reference:​par.txt|par.txt]]|
  
-<m> tmpfcn = 0.86^{delim{|}{T20-15}{|}} ​ </m> 
  
-<m> TPconc = 0.1 </m> 
  
 ==== Links to relevant procedures in the code ====  ​ ==== Links to relevant procedures in the code ====  ​
start/hype_model_description/hype_np_riv_lake.1536574017.txt.gz · Last modified: 2018/09/10 12:06 by cpers