Model help:HydroTrend: Difference between revisions

From CSDMS
No edit summary
 
(35 intermediate revisions by 3 users not shown)
Line 3: Line 3:
1) Log in to the wiki
1) Log in to the wiki
2) Create a new page for each model, by using the following URL:
2) Create a new page for each model, by using the following URL:
   * http://csdms.colorado.edu/wiki/Model help:<modelname>
   * https://csdms.colorado.edu/wiki/Model help:<modelname>
   * Replace <modelname> with the name of a model
   * Replace <modelname> with the name of a model
3) Than follow the link "edit this page"
3) Than follow the link "edit this page"  
   * Now you will see preloaded text. Hit the button "Show Preview" at the bottom, below the edit window
   * Now you will see preloaded text. Hit the button "Show Preview" at the bottom, below the edit window
   * You will see gray text with tag: <span class="remove_this_tag">; Remove the gray text with supporting text and add your own, Hit the button "save" and your help document is all done.
   * You will see gray text with tag: <span class="remove_this_tag">; Remove the gray text with supporting text and add your own, Hit the button "save" and your help document is all done.
Line 18: Line 18:


==Model introduction==
==Model introduction==
HydroTrend is an ANSI-standard C numerical model that creates  synthetic river discharge and sediment load time series as a function of  climate trends and basin morphology and has been used to study the  sediment flux to a basin for basin filling models. As a drainage basin  simulator, the model provides time series of daily discharge hydraulics  at a river mouth, including the sediment load properties. HydroTrend was  designed to provide input to lake or shelf circulation and  sedimentation models (Steckler et al., 1996; Syvitski and Alcott, 1995b), and study the impact of land-sea fluxes given climatic change  scenarios (Moore, 1992; Syvitski and Andrews, 1994).
HydroTrend is an ANSI-standard C numerical model that creates  synthetic river discharge and sediment load time series as a function of  climate trends and basin morphology and has been used to study the  sediment flux to a basin for basin filling models. As a drainage basin  simulator, the model provides time series of daily discharge hydraulics  at a river mouth, including the sediment load properties. HydroTrend was  designed to provide input to lake or shelf circulation and  sedimentation models <ref>Steckler, M., Swift, D., Syvitski, J., Goff, J., and Niedoroda, A., 1996. Modeling the sedimentology and stratigraphy of continental margins. ''Oceanography'', '''9''', 183-188. [[https://www.jstor.org/stable/43924771?seq=1#metadata_info_tab_contents web-pdf]]</ref> <ref>Syvitski, J.P.M., and Alcott, J.M., 1995. DELTA6: Numerical simulation of basin sedimentation affected by slope failure and debris flow runout. In ''Proceedings of the Pierre Beghin International Workshop on Rapid Gravitational Mass Movements'', pp. 180-195. 6-10 December, 1993, Grenoble, France.</ref>, and study the impact of land-sea fluxes given climatic change  scenarios <ref>Moore, R.D., 1992. Hydrological responses to climatic variations in a glacierized watershed: inferences from a conceptual streamflow model. In ''Using Hydrometric Data to Detect and Monitor Climate Change'', Proceedings NHRI Symposium, No. 8, April (1991), pp. 9-20, NHRI Saskatoon.</ref> <ref>Syvitski, J.P.M., and Andrews, J.T., 1994. Climate change: numerical modelling of sedimentation and coastal processes, eastern Canadian Arctic. ''Arctic and Alpine research'', '''26''', 199-212.</ref>.
HydroTrend  simulates the major processes that occur in a river basin, including:  
HydroTrend  simulates the major processes that occur in a river basin, including:  
* Glacierized areas with advances and retreats depending on the climate scenario,  
* Glacierized areas with advances and retreats depending on the climate scenario,  
Line 24: Line 24:
* Rainfall over the remaining portions  of the basin with canopy evaporation,
* Rainfall over the remaining portions  of the basin with canopy evaporation,
* Groundwater recharging and  discharging,
* Groundwater recharging and  discharging,
* The impact of reservoirs.
* The impact of lakes and reservoirs on the stream flow as well as on the sediment load due to sediment retention.
<br>
The Hydrotrend model typically runs at daily timesteps, for a user-defined number of years.


==Model parameters==
==Model parameters==
= Input Files and directories =
{|{{Prettytable}} class = "wikitable unsortable"  cellspacing="0" cellpadding="0" style="margin:0em 0em 0em 0;"
|-
!Parameter!!Description!!Unit
|-valign="top"
|width="20%"| Input directory
|width="60%"| Determine if you want to use the "GUI" interface to provide input parameter values or use a text file with the input parameters by providing the location of the file on the server.
|width="20%"| [-]
|-valign="top"
| Site prefix
| Part of the input and output file name e.g. the name of the geographic location, or project
| [-]
|-valign="top"
| Case prefix
| Part of the input and output file name that provides you the opportunity to do different scenario simulations for e.g. the same location, or project
| [-]
|}


= Run Parameters =
= Run Parameters =
Line 51: Line 36:
|-valign="top"
|-valign="top"
|width="20%"| Run duration
|width="20%"| Run duration
|width="60%"| Number of simulation time steps
|width="60%"| Number of simulation time steps (it is most consistent to use a multiple of 365 days)
|width="20%"| [years]
|width="20%"| [days]
|-
| Hypsometry file
| Used to input hypsometry files for the model
| -
|-
|}
|}


Line 71: Line 61:
| The standard deviation about the trend line that the annual temperatures will have.
| The standard deviation about the trend line that the annual temperatures will have.
| [°C]
| [°C]
|-valign="top"
| Monthly average temperature
| Average temperature for each month
| [°C]
|-valign="top"
| Standard deviation of monthly annual temperature
| The standard deviation about the trend line that the monthly temperatures will have, based on daily data analysis.
| [°C]
|-
|}
|}


Line 78: Line 77:
!Parameter!!Description!!Unit
!Parameter!!Description!!Unit
|-valign="top"
|-valign="top"
|width="20%"| Starting mean annual precipitation
|width="20%"| Total annual precipitation
|width="60%"| Annual total spatial river basin average precipitation rates for the beginning of the simulation
|width="60%"| Annual total spatial river basin average precipitation rates for the beginning of the simulation
|width="20%"| [m/year]
|width="20%"| [m/year]
Line 89: Line 88:
| The standard deviation about the trend line that the annual precipitation will have.
| The standard deviation about the trend line that the annual precipitation will have.
| [m/year]
| [m/year]
|-valign="top"
|width="20%"| Total monthly precipitation for each of the months of the year
|width="60%"| Total spatially-averaged river basin precipitation for each month
|width="20%"| [mm]
|-valign="top"
| Standard deviation of mean monthly precipitation
| The standard deviation for the monthly precipitation (as calculated from daily data).
| [mm]
|-
|}
|}


Line 97: Line 105:
|-valign="top"
|-valign="top"
|width="20%"| Lithology factor
|width="20%"| Lithology factor
|width="60%"| Sediment production varies with lithology, hard versus weak lithology: (0.5 - 3):
|width="60%"| Sediment production varies with lithology, hard versus weak lithology: (0.5 - 3)<ref name="Syvitski_et_al_2007">Syvitski, J.P.M. and Milliman, J.D., 2007, Geology, geography and humans battle for dominance over the delivery of sediment to the coastal ocean. ''J. Geology'' '''115''', 1–19. doi:[http://dx.doi.org/10.1086/509246 10.1086/509246]</ref>:
* L=0.5 for basins comprised principally of hard, acid plutonic and/or high-grade metamorphic rocks.
* L=0.5 for basins comprised principally of hard, acid plutonic and/or high-grade metamorphic rocks.
* L=0.75 for basins comprised of mixed, mostly hard lithology, sometimes including shield material.
* L=0.75 for basins comprised of mixed, mostly hard lithology, sometimes including shield material.
Line 107: Line 115:
|-valign="top"
|-valign="top"
| Anthropogenic facor
| Anthropogenic facor
| Anthropogenic factor (0.3 - 8.0), disturbance to landscape:
| Anthropogenic factor (0.3 - 8.0), disturbance to landscape <ref name="Syvitski_et_al_2007"/>:
* Eh = 0.3: for basins with a high density population (PD > 200 km2) and GNP/capita > $15K/y  
* E<sub>h</sub> = 0.3: for basins with a high density population (PD > 200 km2) and GNP/capita > $15K/y  
* Eh = 1.0: for basins with low human footprint (PD smaller than 50 km2)
* E<sub>h</sub> = 1.0: for basins with low human footprint (PD smaller than 50 km2)
* Eh = 2.0 - 8.0: for basins with the high density population (PD > 200 km2) and GNP/capita is lower than $1K/y  
* E<sub>h</sub> = 2.0 - 8.0: for basins with the high density population (PD > 200 km2) and GNP/capita is lower than $1K/y  
| [-]
| [-]
|-valign="top"
|-valign="top"
| Lapse rate
| Lapse rate
| The change in temperature per change in elevation. Parameter used to determine the snow - rain transition zone
| The change in temperature per change in elevation. Parameter used to determine the snow - rain transition zone.
| [°C/km]
| [°C/km]
|-valign="top"
|-valign="top"
Line 122: Line 130:
|-valign="top"
|-valign="top"
| Change in ELA
| Change in ELA
| Is the linear change per year of the equilibrium line altitude due to a shift in the climate
| Is the linear change per year of the equilibrium line altitude due to a shift in the climate.
| [m/year]
| [m/year]
|-valign="top"
|-valign="top"
Line 134: Line 142:
|-valign="top"
|-valign="top"
| Mean volume of reservoir
| Mean volume of reservoir
| If the reservoir capacity is more than 0.5km<sup>3</sup>, Trap Efficiency (TEbasin) will be calculated based by: (TEbasin = 1.0 - (0.05 / exp(Rvol/RQbar)<sup>0.5</sup>, where Rvol = Reservoir volume and RQbar = the mean inflow discharge)<ref>Vörösmarty, C.J., Meybeck, M., Fekete, B., and Sharma, K. (1997) The potential of neo-Castorization on sediment transport by the global network of rivers. ''Human Impact on Erosion and Sedimentation'', IAHS, '''245''', 261-273.</ref>.<br>If the reservoir capacity is less than 0.5km<sup>3</sup>, (TEbasin = ( 1.0 - (1.0 / (1 + 0.0021 *D * ((Rvol * 1e9) / Rarea)))) where Rvol = Reservoir volume and Rarea (km<sup>2</sup>) = drainage area above the Reservoir) and D, set to 0.1, represents the reservoir characteristics <ref>Verstraeten G., and Poesen, J. (2000) Estimating trap efficiency of small reservoirs and ponds: methods and implications for the assessment of sediment yield. ''Progress in Physical Geography'', '''24''', 219-251</ref>.
| If the reservoir capacity is more than 0.5km<sup>3</sup>, Trap Efficiency (TEbasin) will be calculated based by: (TEbasin = 1.0 - (0.05 / exp(Rvol/RQbar)<sup>0.5</sup>, where Rvol = Reservoir volume and RQbar = the mean inflow discharge)<ref>Vörösmarty, C.J., Meybeck, M., Fekete, B., and Sharma, K. (1997) The potential of neo-Castorization on sediment transport by the global network of rivers. ''Human Impact on Erosion and Sedimentation'', IAHS, '''245''', 261-273.</ref>.<br>If the reservoir capacity is less than 0.5km<sup>3</sup>, (TEbasin = ( 1.0 - (1.0 / (1 + 0.0021 *D * ((Rvol * 1e9) / Rarea)))) where Rvol = Reservoir volume and Rarea (km<sup>2</sup>) = drainage area above the Reservoir) and D, set to 0.1, represents the reservoir characteristics <ref>Verstraeten G., and Poesen, J. (2000) Estimating trap efficiency of small reservoirs and ponds: methods and implications for the assessment of sediment yield. ''Progress in Physical Geography'', '''24''', 219-251. doi:[http://dx.doi.org/10.1177/030913330002400204 10.1177/030913330002400204]</ref>.
| [km<sup>3</sup>]
| [km<sup>3</sup>]
|-valign="top"
|-valign="top"
Line 148: Line 156:
|-valign="top"
|-valign="top"
|width="20%"| k
|width="20%"| k
|width="60%"| River mouth velocity coefficient (Where: v=k*Q<sup>m</sup>, w=a*Q<sup>b</sup>, d=c*Q<sup>f</sup> and the discharge: Q=w*v*d.<br>Therefore: a*c*k = 1 and b+m+f = 1).
|width="60%"| River mouth velocity coefficient (Where: v=k*Q<sup>m</sup>, w=a*Q<sup>b</sup>, d=c*Q<sup>f</sup> and the discharge: Q=w*v*d.<br>Therefore: a*c*k = 1 and b+m+f = 1) <ref name="Leopold_Maddock_1953">Leopold, L.B. and T. Maddock, 1953. “The Hydraulic Geometry of Stream Channels and Some Physiographic Implications”, ''U.S. Geological Survey Professional Paper'' '''252'''.</ref>.
|width="20%"| [m/s]
|width="20%"| [m/s]
|-valign="top"
|-valign="top"
| m
| m
| River mouth velocity exponent (Where: v=k*Q<sup>m</sup>, w=a*Q<sup>b</sup>, d=c*Q<sup>f</sup> and the discharge: Q=w*v*d.<br>Therefore: a*c*k = 1 and b+m+f = 1).
| River mouth velocity exponent (Where: v=k*Q<sup>m</sup>, w=a*Q<sup>b</sup>, d=c*Q<sup>f</sup> and the discharge: Q=w*v*d.<br>Therefore: a*c*k = 1 and b+m+f = 1) <ref name="Leopold_Maddock_1953"/>.
| [-]
| [-]
|-valign="top"
|-valign="top"
| a
| a
| River mouth width coefficient  (Where: v=k*Q<sup>m</sup>, w=a*Q<sup>b</sup>, d=c*Q<sup>f</sup> and the discharge: Q=w*v*d.<br>Therefore: a*c*k = 1 and b+m+f = 1).
| River mouth width coefficient  (Where: v=k*Q<sup>m</sup>, w=a*Q<sup>b</sup>, d=c*Q<sup>f</sup> and the discharge: Q=w*v*d.<br>Therefore: a*c*k = 1 and b+m+f = 1) <ref name="Leopold_Maddock_1953"/>.
| [m]
| [m]
|-valign="top"
|-valign="top"
| b
| b
| River mouth width exponent  (Where: v=k*Q<sup>m</sup>, w=a*Q<sup>b</sup>, d=c*Q<sup>f</sup> and the discharge: Q=w*v*d.<br>Therefore: a*c*k = 1 and b+m+f = 1).
| River mouth width exponent  (Where: v=k*Q<sup>m</sup>, w=a*Q<sup>b</sup>, d=c*Q<sup>f</sup> and the discharge: Q=w*v*d.<br>Therefore: a*c*k = 1 and b+m+f = 1) <ref name="Leopold_Maddock_1953"/>.
| [-]
| [-]
|-valign="top"
|-valign="top"
Line 180: Line 188:
|}
|}


= Groudwater =
= Groundwater =


{|{{Prettytable}} class = "wikitable unsortable"  cellspacing="0"  cellpadding="0" style="margin:0em 0em 0em 0;"
{|{{Prettytable}} class = "wikitable unsortable"  cellspacing="0"  cellpadding="0" style="margin:0em 0em 0em 0;"
Line 199: Line 207:
|-valign="top"
|-valign="top"
| Groundwater coefficient
| Groundwater coefficient
|
| The shallow groundwater (sub-surface storm flow) coefficient is used in draining the groundwater pool to the river.
| [m<sup>3</sup>]
| [m<sup>3</sup>]
|-valign="top"
|-valign="top"
Line 217: Line 225:
!Parameter!!Description!!Unit
!Parameter!!Description!!Unit
|-valign="top"
|-valign="top"
|width="20%"|Output directory
|width="20%"| Output directory
|width="60%"|
|width="60%"| The location where the output files will be stored.
|width="20%"|[-]
|width="20%"| [-]
|-valign="top"
|-valign="top"
|Interval between output files
| Interval between output files
|
| Specify how often you want output to be written. ''Notice: files will be really large if you choose to save every timestep''.
|[-]
| [-]
|-valign="top"
|-valign="top"
|Mean Velocity file
| Mean Velocity file
|
| Output file prefix for variable, velocity, in NetCdf format
|[-]
| [-]
|-valign="top"
|-valign="top"
|Mean Width file
| Mean Width file
|
| Output file prefix for variable, width, in NetCdf format
|[-]
| [-]
|-valign="top"
|-valign="top"
|Mean Depth file
| Mean Depth file
|
| Output file prefix for variable, depth, in NetCdf format
|[-]
| [-]
|-valign="top"
|-valign="top"
|Mean Water Discharge file
| Mean Water Discharge file
|
| Output file prefix for variable, water discharge, in NetCdf format
|[-]
| [-]
|-valign="top"
|-valign="top"
|Mean Sediment Discharge file
| Mean Sediment Discharge file
|
| Output file prefix for variable, sediment discharge, in NetCdf format
|[-]
| [-]
|-valign="top"
|-valign="top"
|Mean Bedload Flux file
| Mean Bedload Flux file
|
| Output file prefix for variable, bedload flux, in NetCdf format
|[-]
| [-]
|}
|}
<headertabs/>
<headertabs/>
Line 269: Line 277:
| = ω B Q<sup>0.31</sup> A<sup>0.5</sup> R T
| = ω B Q<sup>0.31</sup> A<sup>0.5</sup> R T
| for T ≥ 2 °C
| for T ≥ 2 °C
| = Long-term suspended sediment load at the river mouth [kg/s]
| = Long-term suspended sediment load at the river mouth [kg/s] <ref name="Syvitski_et_al_2007"/>
| (2a) or
| (2a) or
|-
|-
| Qs
| Qs
| = 2ω B Q<sup>0.31</sup> A<sup>0.5</sup> R
| = 2ω B Q<sup>0.31</sup> A<sup>0.5</sup> R T
| for T < 2 °C
| for T < 2 °C
| = Long-term suspended sediment load at the river mouth [kg/s]
| = Long-term suspended sediment load at the river mouth [kg/s] <ref name="Syvitski_et_al_2007"/>
| (2b)
| (2b)
|-
|-
Line 281: Line 289:
| = L (1 - T<sub>E</sub>) E<sub>h</sub>
| = L (1 - T<sub>E</sub>) E<sub>h</sub>
|  
|  
| = importance of geology and human factors
| = Expression that capture the importance of geology and human factors. <ref name="Syvitski_et_al_2007"/>
| (3)
| (3)
|-
|-
Line 287: Line 295:
| =  ψ<sub>[''i'']</sub> (Q<sub>[''i'']</sub> / Q)<sup>C<sub>a</sub></sup>
| =  ψ<sub>[''i'']</sub> (Q<sub>[''i'']</sub> / Q)<sup>C<sub>a</sub></sup>
|  
|  
| = Daily suspended sediment load at the river mouth [kg/s]
| = Daily suspended sediment load at the river mouth [kg/s] <ref>Morehead, M.D., Syvitski, J.P.M., Hutton, E.W.H., Peckham, S.D., 2003. Modeling the temporal variability in the flux of sediment from ungauged river basins. ''Global and Planetary Change'', '''39''', 95-110. doi:[http://dx.doi.org/10.1016/S0921-8181(03)00019-5 10.1016/S0921-8181(03)00019-5]</ref>
| (4)
| (4)
|-
|-
Line 293: Line 301:
| = (ρ<sub>s</sub> / ρ<sub>s</sub> - ρ) * (ρ g Q<sub>[''i'']</sub><sup>β</sup> S e<sub>b</sub>) / (g tan λ)
| = (ρ<sub>s</sub> / ρ<sub>s</sub> - ρ) * (ρ g Q<sub>[''i'']</sub><sup>β</sup> S e<sub>b</sub>) / (g tan λ)
| when u ≥ u<sub>cr</sub>  
| when u ≥ u<sub>cr</sub>  
| = Daily bedload at the river mouth [kg/s]
| = Daily bedload at the river mouth [kg/s] <ref>Bagnold, R.A., 1966. An approach to the sediment transport problem from general physics. ''US Geological Survey Professional Paper'', '''422''', 1-37.</ref>
| (5)
| (5)
|-
|-
|}
|}
 
<br>
 
<div class="NavFrame collapsed" style="text-align:left">
'''Nomenclature'''
  <div class="NavHead">Nomenclature</div>
{| {{Prettytable}} class = "wikitable sortable"  
  <div class="NavContent">
{| {{Prettytable}} class="wikitable sortable"
!Symbol!!Description!!Unit
!Symbol!!Description!!Unit
|-
|-
| Q
| Q
| Long-term water discharge
| Long-term water discharge
| [m<sup>3</sup>/s]
| m<sup>3</sup>/s
|-
|-
| Q<sub>''r''</sub>
| Q<sub>''r''</sub>
| Water discharge generated by rainfall
| Water discharge generated by rainfall
| [m<sup>3</sup>/s]
| m<sup>3</sup>/s
|-
|-
| Q<sub>''n''</sub>
| Q<sub>''n''</sub>
| Water discharge generated by nival melt
| Water discharge generated by nival melt
| [m<sup>3</sup>/s]
| m<sup>3</sup>/s
|-
|-
| Q<sub>''ice''</sub>  
| Q<sub>''ice''</sub>  
| Water discharge generated by glacier melt
| Water discharge generated by glacier melt
| [m<sup>3</sup>/s]
| m<sup>3</sup>/s
|-
|-
| Q<sub>''Ev''</sub>
| Q<sub>''Ev''</sub>
| Water discharge loss by evapo-transpiration processes
| Water discharge loss by evapo-transpiration processes
| [m<sup>3</sup>/s]
| m<sup>3</sup>/s
|-
|-
| Q<sub>''g''</sub>
| Q<sub>''g''</sub>
| Water discharge loss or generated by ground water
| Water discharge loss or generated by ground water
| [m<sup>3</sup>/s]
| m<sup>3</sup>/s
|-
|-
| Qs
| Qs
| Long-term suspended sediment load (30yrs or longer)  
| Long-term suspended sediment load (30yrs or longer)  
| [kg/s]
| kg/s
|-
|-
| ω
| ω
Line 337: Line 346:
| A
| A
| Drainage basin area
| Drainage basin area
| [km<sup>2</sup>]
| km<sup>2</sup>
|-
|-
| R
| R
| Drainage basin relief
| Drainage basin relief
| [km]
| km
|-
|-
| T
| T
| Drainage basin temporal and spatial mean temperature
| Drainage basin temporal and spatial mean temperature
| [°C]
| °C
|-
|-
| L
| L
| Lithology factor
| Lithology factor
|[-]
| -
|-
|-
| T<sub>E</sub>
| T<sub>E</sub>
| Trapping efficiency of reservoirs / lakes
| Trapping efficiency of reservoirs / lakes
| [-]
| -
|-
|-
| E<sub>h</sub>
| E<sub>h</sub>
| Anthropogenic factor
| Anthropogenic factor
| [-]
| -
|-
|-
| Qs[''i'']
| Qs[''i'']
| Daily suspended sediment load
| Daily suspended sediment load
| [kg/s]
| kg/s
|-
|-
| Ψ<sub>[''i'']</sub>
| Ψ<sub>[''i'']</sub>
| Daily random variable with a log-normal distribution
| Daily random variable with a log-normal distribution
| [-]
| -
|-
|-
| Q<sub>[''i'']</sub>
| Q<sub>[''i'']</sub>
| Daily water discharge
| Daily water discharge
| [m<sup>3</sup>/s]
| m<sup>3</sup>/s
|-
|-
| C<sub>a</sub>
| C<sub>a</sub>
| Annual rating term coefficient with a normal distrbution
| Annual rating term coefficient with a normal distrbution
| [-]
| -
|-
|-
| Qb<sub>[''i'']</sub>
| Qb<sub>[''i'']</sub>
| Daily bedload
| Daily bedload
| [kg/s]
| kg/s
|-
|-
| ρ<sub>s</sub>
| ρ<sub>s</sub>
| Sand density
| Sand density
| [kg/m<sup>3</sup>]
| kg/m<sup>3</sup>
|-
|-
| ρ
| ρ
| Fluid density
| Fluid density
| [kg/m<sup>3</sup>]
| kg/m<sup>3</sup>
|-
|-
| g
| g
| Acceleration due to gravity
| Acceleration due to gravity
| [m/s<sup>2</sup>]
| m/s<sup>2</sup>
|-
|-
| β
| β
| Bedload rating term
| Bedload rating term
| [-]
| -
|-
|-
| S
| S
| Slope of the riverbed
| Slope of the riverbed
| [m/m]
| m/m
|-
|-
| e<sub>b</sub>
| e<sub>b</sub>
| Bedload efficiency
| Bedload efficiency
| [-]
| -
|-
|-
| λ
| λ
| Limiting angle of response of sediment grains lying on the river bed
| Limiting angle of response of sediment grains lying on the river bed
| [-]
| -
|-
|-
| u
| u
| Stream velocity
| Stream velocity
| [m/s]
| m/s
|-
|-
| u<sub>cr</sub>
| u<sub>cr</sub>
| Critical stream velocity
| Critical stream velocity
| [m/s]
| m/s
|-
|-
|}
|}
  </div>
</div>


==Notes==
==Notes==
<span class="remove_this_tag">Any notes, comments, you want to share with the user</span>  
<span class="remove_this_tag">


'''Note 1:''' Hydrotrend Run Time in WMT is user-defined. All relevant other components in WMT run at daily timestep, thus HydroTrend runs for a number of years but it is listed in days. Example: I want to run HydroTrend for 100 years, I specify the simulation duration to be 36,500. In WMT you can not adjust for leap year (in the stand-alone ANSII C-code for HydroTrend you can run with 'real climate' input values).
Note though that the simulation run time specified by the user is now rounded up to the nearest year or, if it's already an exact multiple of a year, is kept at that year. So,
    366 days -> 2 year
    365 days -> 1 year</span>
<br>
'''Note 2:''' If you want to upload your own hypsometry file, you will need to follow instructions on the file format precisely. More information can be found here:
[[https://csdms.colorado.edu/wiki/Model:HydroTrend#HYDRO0.HYPS| Hypsometry file format nder metadata for this model]]
<span class="remove_this_tag">Numerical scheme</span>
<span class="remove_this_tag">Numerical scheme</span>


==Examples==
==Examples==
Line 427: Line 447:


<span class="remove_this_tag">Follow the next steps to include images / movies of simulations:</span>
<span class="remove_this_tag">Follow the next steps to include images / movies of simulations:</span>
* <span class="remove_this_tag">Upload file: http://csdms.colorado.edu/wiki/Special:Upload</span>
* <span class="remove_this_tag">Upload file: https://csdms.colorado.edu/wiki/Special:Upload</span>
* <span class="remove_this_tag">Create link to the file on your page: <nowiki>[[Image:<file name>]]</nowiki>.</span>
* <span class="remove_this_tag">Create link to the file on your page: <nowiki>[[Image:<file name>]]</nowiki>.</span>



Latest revision as of 10:53, 16 October 2019

The CSDMS Help System

HydroTrend

Climate driven hydrological transport model

Model introduction

HydroTrend is an ANSI-standard C numerical model that creates synthetic river discharge and sediment load time series as a function of climate trends and basin morphology and has been used to study the sediment flux to a basin for basin filling models. As a drainage basin simulator, the model provides time series of daily discharge hydraulics at a river mouth, including the sediment load properties. HydroTrend was designed to provide input to lake or shelf circulation and sedimentation models [1] [2], and study the impact of land-sea fluxes given climatic change scenarios [3] [4]. HydroTrend simulates the major processes that occur in a river basin, including:

  • Glacierized areas with advances and retreats depending on the climate scenario,
  • Snow accumulation in the winter and melt in the subsequent spring/summer,
  • Rainfall over the remaining portions of the basin with canopy evaporation,
  • Groundwater recharging and discharging,
  • The impact of lakes and reservoirs on the stream flow as well as on the sediment load due to sediment retention.


The Hydrotrend model typically runs at daily timesteps, for a user-defined number of years.

Model parameters

Parameter Description Unit
Run duration Number of simulation time steps (it is most consistent to use a multiple of 365 days) [days]
Hypsometry file Used to input hypsometry files for the model -
Parameter Description Unit
Starting mean annual temperature Mean annual temperature at the river mouth at the start of the simulation. (Notice: this is not a spatial average basin width parameter!) [°C]
Change in mean annual temperature The trend or change per year in the annual temperature [°C/year]
Standard deviation of mean annual temperature The standard deviation about the trend line that the annual temperatures will have. [°C]
Monthly average temperature Average temperature for each month [°C]
Standard deviation of monthly annual temperature The standard deviation about the trend line that the monthly temperatures will have, based on daily data analysis. [°C]
Parameter Description Unit
Total annual precipitation Annual total spatial river basin average precipitation rates for the beginning of the simulation [m/year]
change in mean annual precipitation The trend or change in the total annual precipitation [m/year/year]
Standard deviation of mean annual precipitation The standard deviation about the trend line that the annual precipitation will have. [m/year]
Total monthly precipitation for each of the months of the year Total spatially-averaged river basin precipitation for each month [mm]
Standard deviation of mean monthly precipitation The standard deviation for the monthly precipitation (as calculated from daily data). [mm]
Parameter Description Unit
Lithology factor Sediment production varies with lithology, hard versus weak lithology: (0.5 - 3)[5]:
  • L=0.5 for basins comprised principally of hard, acid plutonic and/or high-grade metamorphic rocks.
  • L=0.75 for basins comprised of mixed, mostly hard lithology, sometimes including shield material.
  • L=1.0 for basins comprised of volcanic, mostly basaltic rocks, or carbonate outcrops, or mixture of hard and soft lithologies.
  • L=1.5 for basins characterized by a predominance of softer lithologies, but having a significant area of harder lithologies.
  • L=2 for fluvial systems draining a high proportion of sedimentary rocks, unconsolidated sedimentary cover, or alluvial deposits.
  • L=3 for basins having an abundance of exceptionally weak material, such as crushed rock, or loess deposits, or shifting sand dunes.
[-]
Anthropogenic facor Anthropogenic factor (0.3 - 8.0), disturbance to landscape [5]:
  • Eh = 0.3: for basins with a high density population (PD > 200 km2) and GNP/capita > $15K/y
  • Eh = 1.0: for basins with low human footprint (PD smaller than 50 km2)
  • Eh = 2.0 - 8.0: for basins with the high density population (PD > 200 km2) and GNP/capita is lower than $1K/y
[-]
Lapse rate The change in temperature per change in elevation. Parameter used to determine the snow - rain transition zone. [°C/km]
Starting ELA Glacier equilibrium line altitude is the long-term balance point along a glacier. The ELA is where the amount of accumulated snow and ablated water are equal. [m]
Change in ELA Is the linear change per year of the equilibrium line altitude due to a shift in the climate. [m/year]
Dry precipitation (nival and ice) evaporation fraction The percentage of the dry precipitation (nival&ice) which will be evaporated. [-]
River length The length of the main stem of the river to determine the lag time before water reaches the river outlet. [km]
Mean volume of reservoir If the reservoir capacity is more than 0.5km3, Trap Efficiency (TEbasin) will be calculated based by: (TEbasin = 1.0 - (0.05 / exp(Rvol/RQbar)0.5, where Rvol = Reservoir volume and RQbar = the mean inflow discharge)[6].
If the reservoir capacity is less than 0.5km3, (TEbasin = ( 1.0 - (1.0 / (1 + 0.0021 *D * ((Rvol * 1e9) / Rarea)))) where Rvol = Reservoir volume and Rarea (km2) = drainage area above the Reservoir) and D, set to 0.1, represents the reservoir characteristics [7].
[km3]
Drainage area of reservoir The upstream area of the river basin that drains into the reservoir. [km2]
Parameter Description Unit
k River mouth velocity coefficient (Where: v=k*Qm, w=a*Qb, d=c*Qf and the discharge: Q=w*v*d.
Therefore: a*c*k = 1 and b+m+f = 1) [8].
[m/s]
m River mouth velocity exponent (Where: v=k*Qm, w=a*Qb, d=c*Qf and the discharge: Q=w*v*d.
Therefore: a*c*k = 1 and b+m+f = 1) [8].
[-]
a River mouth width coefficient (Where: v=k*Qm, w=a*Qb, d=c*Qf and the discharge: Q=w*v*d.
Therefore: a*c*k = 1 and b+m+f = 1) [8].
[m]
b River mouth width exponent (Where: v=k*Qm, w=a*Qb, d=c*Qf and the discharge: Q=w*v*d.
Therefore: a*c*k = 1 and b+m+f = 1) [8].
[-]
Average river mouth velocity The average stream velocity, used in the model to determine the flow routing [m/s]
Constant annual baseflow Defines the constant annual baseflow which occurs in the basin. This is analogous to a deep groundwater pool. Baseflow should be bigger than 0 and smaller than the total precipitation discharge. [m3/s]
Trapping efficiency Deltaic area typically trap sediment e.g. due to the low slope angle of the area. The fraction of deltaic trapping (0 - 1) can be provided here. [-]
Delta gradient Delta gradient determines the slope of the riverbed that determines the amount of bedload reaching the river mouth. [m/m]
Parameter Description Unit
Maximum groundwater storage Maximum capacity spacial averaged shallow ground water pole of the drainage basin [m3]
Minimum groundwater storage Minimum capacity spacial averaged shallow ground water pole of the drainage basin [m3]
Initial groundwater storage Spatial averaged shallow ground water pole of the drainage basin at the start of the simulation. [m3]
Groundwater coefficient The shallow groundwater (sub-surface storm flow) coefficient is used in draining the groundwater pool to the river. [m3]
Groundwater exponent The shallow groundwater (sub-surface storm flow) exponent is used in draining the groundwater pool to the river. [-]
Ksat Saturated hydraulic conductivity describes water movement through saturated media. See table for saturated hydraulic conductivity values in relation to texture. [mm/day]
Parameter Description Unit
Output directory The location where the output files will be stored. [-]
Interval between output files Specify how often you want output to be written. Notice: files will be really large if you choose to save every timestep. [-]
Mean Velocity file Output file prefix for variable, velocity, in NetCdf format [-]
Mean Width file Output file prefix for variable, width, in NetCdf format [-]
Mean Depth file Output file prefix for variable, depth, in NetCdf format [-]
Mean Water Discharge file Output file prefix for variable, water discharge, in NetCdf format [-]
Mean Sediment Discharge file Output file prefix for variable, sediment discharge, in NetCdf format [-]
Mean Bedload Flux file Output file prefix for variable, bedload flux, in NetCdf format [-]

Uses ports

This will be something that the CSDMS facility will add

Provides ports

This will be something that the CSDMS facility will add

Main equations

Q = Qr + Qn + Qice - QEv ± Qg = Water discharge at the river mouth [m3/s] (1)
Qs = ω B Q0.31 A0.5 R T for T ≥ 2 °C = Long-term suspended sediment load at the river mouth [kg/s] [5] (2a) or
Qs = 2ω B Q0.31 A0.5 R T for T < 2 °C = Long-term suspended sediment load at the river mouth [kg/s] [5] (2b)
B = L (1 - TE) Eh = Expression that capture the importance of geology and human factors. [5] (3)
(Qs[i] / Qs) = ψ[i] (Q[i] / Q)Ca = Daily suspended sediment load at the river mouth [kg/s] [9] (4)
Qb[i] = (ρs / ρs - ρ) * (ρ g Q[i]β S eb) / (g tan λ) when u ≥ ucr = Daily bedload at the river mouth [kg/s] [10] (5)


Notes

Note 1: Hydrotrend Run Time in WMT is user-defined. All relevant other components in WMT run at daily timestep, thus HydroTrend runs for a number of years but it is listed in days. Example: I want to run HydroTrend for 100 years, I specify the simulation duration to be 36,500. In WMT you can not adjust for leap year (in the stand-alone ANSII C-code for HydroTrend you can run with 'real climate' input values).

Note though that the simulation run time specified by the user is now rounded up to the nearest year or, if it's already an exact multiple of a year, is kept at that year. So,

   366 days -> 2 year
   365 days -> 1 year 


Note 2: If you want to upload your own hypsometry file, you will need to follow instructions on the file format precisely. More information can be found here: [Hypsometry file format nder metadata for this model] Numerical scheme

Examples

An example run with input parameters, BLD files, as well as a figure / movie of the output

Follow the next steps to include images / movies of simulations:

See also: Help:Images or Help:Movies

Developer(s)

Albert Kettner

Key HydroTrend Papers

  • Kettner, A.J., and Syvitski, J.P.M., 2008. HydroTrend version 3.0: a Climate-Driven Hydrological Transport Model that Simulates Discharge and Sediment Load leaving a River System. Computers & Geosciences, 34(10), 1170-1183, doi:10.1016/j.cageo.2008.02.008.
  • Syvitski, J.P.M., Morehead, M.D. and Nicholson, M., 1998. HYDROTREND: A Climate-driven Hydrologic-Transport Model for Predicting Discharge and Sediment Loads to Lakes or Oceans. Computers & Geosciences, 24(1), 51-68, doi:10.1016/S0098-3004(97)00083-6.
  • Syvitski, J.P.M., and J.M. Alcott, 1995. RIVER3: Simulation of River Discharge and Sediment Transport. Computers and Geosciences, 21(1), 89-101, doi:10.1016/0098-3004(94)00062-Y.

References

  1. Steckler, M., Swift, D., Syvitski, J., Goff, J., and Niedoroda, A., 1996. Modeling the sedimentology and stratigraphy of continental margins. Oceanography, 9, 183-188. [web-pdf]
  2. Syvitski, J.P.M., and Alcott, J.M., 1995. DELTA6: Numerical simulation of basin sedimentation affected by slope failure and debris flow runout. In Proceedings of the Pierre Beghin International Workshop on Rapid Gravitational Mass Movements, pp. 180-195. 6-10 December, 1993, Grenoble, France.
  3. Moore, R.D., 1992. Hydrological responses to climatic variations in a glacierized watershed: inferences from a conceptual streamflow model. In Using Hydrometric Data to Detect and Monitor Climate Change, Proceedings NHRI Symposium, No. 8, April (1991), pp. 9-20, NHRI Saskatoon.
  4. Syvitski, J.P.M., and Andrews, J.T., 1994. Climate change: numerical modelling of sedimentation and coastal processes, eastern Canadian Arctic. Arctic and Alpine research, 26, 199-212.
  5. 5.0 5.1 5.2 5.3 5.4 Syvitski, J.P.M. and Milliman, J.D., 2007, Geology, geography and humans battle for dominance over the delivery of sediment to the coastal ocean. J. Geology 115, 1–19. doi:10.1086/509246
  6. Vörösmarty, C.J., Meybeck, M., Fekete, B., and Sharma, K. (1997) The potential of neo-Castorization on sediment transport by the global network of rivers. Human Impact on Erosion and Sedimentation, IAHS, 245, 261-273.
  7. Verstraeten G., and Poesen, J. (2000) Estimating trap efficiency of small reservoirs and ponds: methods and implications for the assessment of sediment yield. Progress in Physical Geography, 24, 219-251. doi:10.1177/030913330002400204
  8. 8.0 8.1 8.2 8.3 Leopold, L.B. and T. Maddock, 1953. “The Hydraulic Geometry of Stream Channels and Some Physiographic Implications”, U.S. Geological Survey Professional Paper 252.
  9. Morehead, M.D., Syvitski, J.P.M., Hutton, E.W.H., Peckham, S.D., 2003. Modeling the temporal variability in the flux of sediment from ungauged river basins. Global and Planetary Change, 39, 95-110. doi:10.1016/S0921-8181(03)00019-5
  10. Bagnold, R.A., 1966. An approach to the sediment transport problem from general physics. US Geological Survey Professional Paper, 422, 1-37.

Links

Model:HydroTrend