HEC-HMS Users Manual 4.2
HEC-HMS Users Manual 4.2
HEC-HMS Users Manual 4.2
User's Manual
Version 4.2
August 2016
The public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching
existing data sources, gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this
burden estimate or any other aspect of this collection of information, including suggestions for reducing this burden, to the Department of Defense, Executive
Services and Communications Directorate (0704-0188). Respondents should be aware that notwithstanding any other provision of law, no person shall be
subject to any penalty for failing to comply with a collection of information if it does not display a currently valid OMB control number.
PLEASE DO NOT RETURN YOUR FORM TO THE ABOVE ORGANIZATION.
1. REPORT DATE (DD-MM-YYYY) 2. REPORT TYPE 3. DATES COVERED (From - To)
August 2016 Computer Program
Document
4. TITLE AND SUBTITLE 5a. CONTRACT NUMBER
Hydrologic Modeling System HEC-HMS
5b. GRANT NUMBER
User's Manual
5c. PROGRAM ELEMENT NUMBER
14. ABSTRACT
The program features an integrated work environment with data entry utilities, computation
engine, and results reporting tools. A graphical user interface allows seamless movement
between the different parts of the program. Simulation results are stored in the Data Storage
System HEC-DSS and can be used in conjunction with other software for studies of water
availability, urban drainage, flow forecasting, future urbanization impact, reservoir spillway
design, flood damage reduction, floodplain regulation, and systems operation.
15. SUBJECT TERMS
Hydrology, watershed, precipitation runoff, river routing, flood frequency, flood control,
water supply, computer simulation, environmental restoration.
16. SECURITY CLASSIFICATION OF: 17. LIMITATION 18. NUMBER 19a. NAME OF RESPONSIBLE
a. REPORT b. ABSTRACT c. THIS PAGE OF OF PERSON
ABSTRACT PAGES
U U U UU 612 19b. TELEPHONE NUMBER
User's Manual
Version 4.2
August 2016
Phone 530.756.1104
Fax 530.756.8250
Email [email protected]
Hydrologic Modeling System HEC-HMS, User's Manual
Please consider the environment before printing this document, and recycle if you do.
ii
Table of Contents
PREFACE XIII
INTRODUCTION 1
SCOPE ................................................................................................................................................. 1
HISTORY .............................................................................................................................................. 1
CAPABILITIES........................................................................................................................................ 3
Watershed Physical Description .................................................................................................... 3
Meteorology Description ................................................................................................................ 4
Hydrologic Simulation .................................................................................................................... 6
Parameter Estimation .................................................................................................................... 6
Analyzing Simulations .................................................................................................................... 6
Forecasting Future Flows .............................................................................................................. 6
Sediment and Water Quality .......................................................................................................... 7
GIS Connection .............................................................................................................................. 7
LIMITATIONS ......................................................................................................................................... 7
Model Formulation ......................................................................................................................... 8
Flow Representation ...................................................................................................................... 8
DOCUMENTATION CONVENTIONS ........................................................................................................... 9
REFERENCES ....................................................................................................................................... 9
INSTALLING AND RUNNING THE PROGRAM 11
OPERATING SYSTEM REQUIREMENTS .................................................................................................. 11
HARDWARE REQUIREMENTS AND RECOMMENDATIONS ......................................................................... 11
INSTALLATION ..................................................................................................................................... 12
Microsoft Windows Operating System ......................................................................................... 12
Linux Operating System ............................................................................................................... 13
RUNNING THE PROGRAM ..................................................................................................................... 14
Microsoft Windows Operating System ......................................................................................... 15
Linux Operating System ............................................................................................................... 15
COMMAND LINE OPERATION ................................................................................................................ 15
MANAGING MEMORY ALLOCATION ....................................................................................................... 17
ADDITIONAL RESOURCES .................................................................................................................... 17
OVERVIEW 19
PROGRAM SCREEN ............................................................................................................................. 19
Menu System ............................................................................................................................... 19
Toolbar ......................................................................................................................................... 25
Watershed Explorer ..................................................................................................................... 25
Desktop ........................................................................................................................................ 27
Component Editor ........................................................................................................................ 28
Message Log................................................................................................................................ 29
PROGRAM SETTINGS .......................................................................................................................... 29
DATA CONVENTIONS ........................................................................................................................... 31
Saving Properties ......................................................................................................................... 31
Number Formatting ...................................................................................................................... 32
Date and Time Formatting ........................................................................................................... 32
Units Conversion.......................................................................................................................... 33
Interpolation ................................................................................................................................. 34
APPLICATION STEPS ........................................................................................................................... 34
Create a New Project ................................................................................................................... 34
Enter Shared Project Data ........................................................................................................... 35
Describe the Physical Watershed ................................................................................................ 37
Describe the Meteorology ............................................................................................................ 41
Enter Simulation Time Windows .................................................................................................. 42
Simulate and View Results .......................................................................................................... 42
iii
Table of Contents
iv
Table of Contents
v
Table of Contents
vi
Table of Contents
vii
Table of Contents
viii
Table of Contents
ix
Table of Contents
x
Table of Contents
xi
Table of Contents
xii
Preface
PREFACE
This Manual is not intended to teach you how to do hydrologic
engineering or even hydrology. It does not describe the mathematical
equations for the various models included in the program. So what
does it do? This Manual will teach you how to use the various features
and capabilities of the program. It works very well to simply read the
Manual through starting at the beginning. If you read the Manual in
front of your computer with the program up and running, it will work
even better. However, the Manual works equally well as an occasional
reference when you cannot remember exactly how to perform a certain
task or need to check the parameter definitions for a particular method.
The scope of this Manual does not mean that we think engineering
applications or mathematical analysis are unimportant. In fact, both of
those things are vital to producing good engineering plans and
designs. We feel they are so important that we have created a
separate Manual for each of them. The Technical Reference Manual
provides detailed descriptions of each of the models included in the
program. You can expect to find the mathematical derivation of the
model equations, details on the numerical schemes employed in the
program to solve the equations, and specific guidance on parameter
estimation. Consequently, it focuses less on using the program and
more on understanding the science of hydrology. The Applications
Guide provides practical suggestions for using the program to perform
engineering work. We selected a number of typical projects that
engineers often encounter and showed how the program can be used
to provide real answers. Consequently, it focuses less on using the
program and more on the engineering process.
Many engineers, computer specialists, and student interns have
contributed to the success of this project. Each one has made
valuable contributions that enhance the overall success of the
program. Nevertheless, the completion of this version of the program
was overseen by David J. Harris while Christopher N. Dunn was
director of the Hydrologic Engineering Center. Development and
testing of this release was led by William A. Scharffenberg.
xiii
Preface
xiv
Chapter 1 Introduction
CHAPTER 1
Introduction
The Hydrologic Modeling System is designed to simulate the
precipitation-runoff processes of dendritic watershed systems. It is
designed to be applicable in a wide range of geographic areas for
solving the widest possible range of problems. This includes large
river basin water supply and flood hydrology, and small urban or
natural watershed runoff. Hydrographs produced by the program are
used directly or in conjunction with other software for studies of water
availability, urban drainage, flow forecasting, future urbanization
impact, reservoir spillway design, flood damage reduction, floodplain
regulation, and systems operation.
Scope
The program is a generalized modeling system capable of
representing many different watersheds. A model of the watershed is
constructed by separating the hydrologic cycle into manageable pieces
and constructing boundaries around the watershed of interest. Any
mass or energy flux in the cycle can then be represented with a
mathematical model. In most cases, several model choices are
available for representing each flux. Each mathematical model
included in the program is suitable in different environments and under
different conditions. Making the correct choice requires knowledge of
the watershed, the goals of the hydrologic study, and engineering
judgment.
The program features a completely integrated work environment
including a database, data entry utilities, computation engine, and
results reporting tools. A graphical user interface allows the seamless
movement between the different parts of the program. Program
functionality and appearance are the same across all supported
platforms.
History
The computation engine draws on over 30 years experience with
hydrologic simulation software. Many algorithms from HEC-1 (HEC,
1998), HEC-1F (HEC, 1989), PRECIP (HEC, 1989), and HEC-IFH
(HEC, 1992) have been modernized and combined with new
algorithms to form a comprehensive library of simulation routines.
Future versions of the program will continue to modernize desirable
1
Chapter 1 Introduction
2
Chapter 1 Introduction
Capabilities
The program has an extensive array of capabilities for conducting
hydrologic simulation. Many of the most common methods in
hydrologic engineering are included in such a way that they are easy to
use. The program does the difficult work and leaves the user free to
concentrate on how best to represent the watershed environment.
3
Chapter 1 Introduction
Meteorology Description
Meteorologic data analysis is performed by the meteorologic model
and includes shortwave radiation, longwave radiation, precipitation,
evapo-transpiration, and snowmelt. Not all of these components are
required for all simulations. Simple event simulations require only
precipitation, while continuous simulation additionally requires evapo-
transpiration. Generally, snowmelt is only required when working with
watersheds in cold climates. Selection of energy-balance methods for
evapo-transpiration or snowmelt also requires radiation methods.
Four different methods for analyzing historical precipitation are
included. The user-specified hyetograph method is for precipitation
data analyzed outside the program. The gage weights method uses
an unlimited number of recording and non-recording gages. The
Thiessen technique is one possibility for determining the weights. The
inverse distance method addresses dynamic data problems. An
unlimited number of recording and non-recording gages can be used to
4
Chapter 1 Introduction
5
Chapter 1 Introduction
Hydrologic Simulation
The time span of a simulation is controlled by control specifications.
Control specifications include a starting date and time, ending date and
time, and a time interval.
A simulation run is created by combining a basin model, meteorologic
model, and control specifications. Run options include a precipitation
or flow ratio, capability to save all basin state information at a point in
time, and ability to begin a simulation run from previously saved state
information.
Simulation results can be viewed from the basin map. Global and
element summary tables include information on peak flow, total
volume, and other variables. A time-series table and graph are
available for elements. Results from multiple elements and multiple
simulation runs can also be viewed. All graphs and tables can be
printed.
Parameter Estimation
Most parameters for methods included in subbasin and reach elements
can be estimated automatically using optimization trials. Observed
discharge must be available for at least one element before
optimization can begin. Parameters at any element upstream of the
observed flow location can be estimated. Seven different objective
functions are available to estimate the goodness-of-fit between the
computed results and observed discharge. Two different search
methods can be used to minimize the objective function. Constraints
can be imposed to restrict the parameter space of the search method.
Analyzing Simulations
Analysis tools are designed to work with simulation runs to provide
additional information or processing. Currently, the only tool is the
depth-area analysis tool. It works with simulation runs that have a
meteorologic model using the frequency storm method. Given a
selection of elements, the tool automatically adjusts the storm area and
generates peak flows represented by the correct storm areas.
6
Chapter 1 Introduction
GIS Connection
The power and speed of the program make it possible to represent
watersheds with hundreds of hydrologic elements. Traditionally, these
elements would be identified by inspecting a topographic map and
manually identifying drainage boundaries. While this method is
effective, it is prohibitively time consuming when the watershed will be
represented with many elements. A geographic information system
(GIS) can use elevation data and geometric algorithms to perform the
same task much more quickly. A GIS companion product has been
developed to aid in the creation of basin models for such projects. It is
called the Geospatial Hydrologic Modeling Extension (HEC-GeoHMS)
and can be used to create basin and meteorologic models for use with
the program.
Limitations
Every simulation system has limitations due to the choices made in the
design and development of the software. The limitations that arise in
this program are due to two aspects of the design: simplified model
formulation, and simplified flow representation. Simplifying the model
formulation allows the program to complete simulations very quickly
while producing accurate and precise results. Simplifying the flow
representation aids in keeping the compute process efficient and
reduces duplication of capability in the HEC software suite.
7
Chapter 1 Introduction
Model Formulation
All of the mathematical models included in the program are
deterministic. This means that the boundary conditions, initial
conditions, and parameters of the models are assumed to be exactly
known. This guarantees that every time a simulation is computed it will
yield exactly the same results as all previous times it was computed.
Deterministic models are sometimes compared to stochastic models
where the same boundary conditions, initial conditions, and
parameters are represented with probabilistic distributions. Plans are
underway to develop a stochastic capability through the addition of a
Monte Carlo analysis tool.
All of the mathematical models included in the program use constant
parameter values, that is, they are assumed to be time stationary.
During long periods of time it is possible for parameters describing a
watershed to change as the result of human or other processes at
work in the watershed. These parameter trends cannot be included in
a simulation at this time. There is a limited capability to break a long
simulation into smaller segments and manually change parameters
between segments. Plans are underway to develop a variable
parameter capability, through an as yet undetermined means.
All of the mathematical models included in the program are uncoupled.
The program first computes evapotranspiration and then computes
infiltration. In the physical world, the amount of evapotranspiration
depends on the amount of soil water. The amount of infiltration also
depends on the amount of soil water. However, evapotranspiration
removes water from the soil at the same time infiltration adds water to
the soil. To solve the problem properly the evapotranspiration and
infiltration processes should be simulated simultaneously with the
mathematical equations for both processes numerically linked. At this
time the program computes evapo-transpiration and infiltration in a
loosely coupled arrangement and consequent errors due to the limited
coupling are minimized as much as possible by using a small time
interval for calculations. While preparations have been made to
support the inclusion of fully coupled plant-surface-soil models, none
have been added at this time.
Flow Representation
The design of the basin model only allows for dendritic stream
networks. The best way to visualize a dendritic network is to imagine a
tree. The main tree trunk, branches, and twigs correspond to the main
river, tributaries, and headwater streams in a watershed. The key idea
is that a stream does not separate into two streams. The basin model
allows each hydrologic element to have only one downstream
connection so it is not possible to split the outflow from an element into
8
Chapter 1 Introduction
Documentation Conventions
The following conventions are used throughout the manual to describe
the graphical user interface:
References
Corps of Engineers. 1952. Engineer Manual 1110-2-1411: Standard
Project Flood Determinations. U.S. Army, Washington, DC.
Hydrologic Engineering Center. June 1998. HEC-1 Flood Hydrograph
Package: User's Manual. U.S. Army Corps of Engineers, Davis, CA.
9
Chapter 1 Introduction
10
Chapter 2 Installing and Running the Program
CHAPTER 2
1 GB of memory recommended.
11
Chapter 2 Installing and Running the Program
Installation
Installation packages for the program are available from the Hydrologic
Engineering Center (HEC) website where the current version of the
program is always available. Old versions of the program are archived
and can be downloaded. However, old versions are not maintained,
contain bugs and errors, and may not function correctly with current
versions of the supported operating systems.
12
Chapter 2 Installing and Running the Program
13
Chapter 2 Installing and Running the Program
The installation program for this operating system is run from the
command line, so terminal access is required. If you would like to
install the program for all users of the system, you will need root
access during the installation. Otherwise, the program package can be
installed into any user writable directory on the system and be run from
that location. The program does not require root privileges to run. In
general, the changes required to install on this operating system
require the skills of a system administrator. Please contact your
system administrator to install the program for you if youre unsure
how. You may need to refer your administrator to this section of the
manual.
After you have obtained the install package and proper permissions,
use the following steps to install the program:
1. Download hechms42.bin to your Linux box. Keep note where
you put it. In the example below, we put it into /tmp.
2. Change directory to the location where you want to install the
HEC-HMS 4.2 program files.
cd <location where HEC-HMS 4.2 is to be installed>
4. Run hechm42.bin:
</tmp>/hechms42.bin
5. You will be presented with the Terms and Conditions for use of
HEC-HMS. Read the TCU and then type in "yes" if you agree.
6. Once you have accepted the TCU, HEC-HMS 4.2 will be
extracted to <current working directory>/hechms42.
HEC-HMS 4.2 is now installed. When future versions of the program
become available, you may have each version separately installed on
your computer. You will need to carefully organize the installation
locations so that each version can be kept separate.
14
Chapter 2 Installing and Running the Program
separately for each user, unless the users take steps to make the
projects available to all users.
If you would like to run HEC-HMS 4.2 from any location on your file
system without changing directory to the install location first, you must
make a change to the hms script. There is a variable defined in the
script called HMS_HOME which specifies where HEC-HMS is
installed. By default, it is set to the current directory ("."). Change that
variable to the fully qualified path of where HEC-HMS 4.2 installed.
After changing the script to reflect the full path to the program install
directory, you need to modify the PATH variable to include the program
directory:
set path = ($path <HEC-HMS Install>)
15
Chapter 2 Installing and Running the Program
Once you have created the script file, it can be used with the program
from the command line. The program will start and automatically
process the script. The first line is used to setup the scripting
environment and make the program data model accessible to the
script. The second line opens an existing project and the third line
computes an existing simulation run. The final line of the script exits
the program.
To use a script on the Microsoft Windows operating system, begin by
opening a command window and changing directories to the
installation folder. The installation folder is not standardized and
depends on where you chose to install the program. One possibility
would look like the following:
C:\Program Files\HEC\HEC-HMS\4.2>
At the command prompt, type the following to launch the program and
run the script, where the last argument is the complete path to a script:
hec-hms.cmd s C:\hmsproj\Tenk\compute.script
At the console, type the following to launch the program and run the
script, where the last argument is the complete path to the script:
hms s /usr/smith/hmsproj/tenk/compute.script
The program will not be visible while it is running the script. However,
the commands in the script will be carried out. Any messages
generated while computing the simulation run will be written to the log
16
Chapter 2 Installing and Running the Program
file. All results will be stored in the output Data Storage System (DSS)
file. Inspection of the log file will reveal any errors, warning, or notes
and results can be read from the DSS file.
Linux:
<HEC-HMS Install>/hms
In all three files, search for the string -Xmx512M and replace the 512
with the number of MBs you would like to allow the program to use.
Additional Resources
The program includes an online help system that is automatically
installed when the program is installed. The help system is equivalent
to the User's Manual, Technical Reference Manual, and Applications
Guide. The various documents are also available separately.
Three sample projects are included with the program. The "Castro"
project shows how the program can be used for basic hydrology. The
"Tenk" project demonstrates the capability of the program to use
gridded precipitation. Finally, the "Tifton" project exhibits continuous
simulation with the soil moisture accounting method. The projects are
17
Chapter 2 Installing and Running the Program
18
Chapter 3 Overview
CHAPTER 3
Overview
This chapter describes the basics of working with the program. It
includes descriptions of the main parts of the interface. Specific details
of when and how data are saved is also included. Conventions are
provided for the formatting of input data, the use of units, and
interpolation. An outline of the best way to use the program is also
provided.
Program Screen
The program screen contains a title bar, menu system, toolbar, and
four panes. These panes will be referred to as the Watershed
Explorer, Desktop, Component Editor, and the Message Log as shown
in Figure 1. The title bar displays the version of the program used and
the location of the currently-open project. The other parts of the
program screen are discussed in detail in this chapter.
Menu System
The menu system contains several menus to help you use the
program. Each menu contains a list of related commands. For
example, the Parameters menu contains a list of commands to open
global parameter tables for viewing and editing parameters required by
hydrologic elements in the selected basin model. Items in an individual
menu are inactive (cannot be selected) if the command can not be
carried out by the program at the current time.
Commands for managing the opened project are available from the
File menu. File menu items and the resulting actions are provided in
Table 1. The last four projects opened are shown at the bottom of the
File menu. Click on one of the project names to open the project.
The Edit menu contains commands for editing hydrologic elements in
the selected basin model. If no basin model is selected, then all
commands in this menu are inactive. Edit menu commands and the
resulting actions are provided in Table 2. Additional Edit menu
includes commands select mouse tools for use in the Basin Model
Map.
The View menu contains a list of commands for working in the basin
map. These commands are inactive if no basin model is open in the
19
Chapter 3 Overview
Desktop. A list of View menu items and the resulting actions are
provided in Table 3.
Figure 1. The main program screen with Watershed Explorer in the upper
left, Component Editor in the lower left, Message Log at the
bottom, and Desktop using the remaining area.
20
Chapter 3 Overview
21
Chapter 3 Overview
22
Chapter 3 Overview
23
Chapter 3 Overview
Table 5. Continued.
Table 6. Commands available from the Results menu. Items marked with
(1) are available for simulations runs, items with (2) are available
for optimization trials, and items with (3) are available for depth-
area analyses.
24
Chapter 3 Overview
Toolbar
Like the menu system, the toolbar groups tools with similar uses. Also,
tools are inactive if the resulting action for the tool can not be carried
out by the program. Table 8 contains a description of tools on the
toolbar.
Watershed Explorer
The Watershed Explorer was developed to provide quick access to all
components and simulations in a project. For example, the user can
easily navigate from a basin model to a precipitation gage and then to
a meteorologic model without using menu commands or opening
additional windows. The Watershed Explorer is divided into three
parts: Components, Compute, and Results.
A typical view of the "Components" tab of the Watershed Explorer is
shown in Figure 2. The hierarchal structure of model components,
such as basin models, meteorologic models, etc., is available from the
"Components" tab. The Watershed Explorer organizes model
components into individual folders. If all component types were in a
project, then the Components tab would contain six main folders, Basin
Models, Meteorologic Models, Control Specifications, Time-Series
Data, Paired Data, and Grid Data. The Watershed Explorer expands
when one of these main component folders is selected. When a
component in a main folder is selected, the Watershed Explorer
expands even more to show sub-components. For example, when the
25
Chapter 3 Overview
Tool Action
26
Chapter 3 Overview
Desktop
The Desktop holds a variety of windows including global parameter
editors, and most importantly the basin map. Result windows including
graphs, summary tables, and time-series tables can be shown in the
Desktop, or optionally, outside the Desktop. All other windows cannot
be moved outside of the Desktop area. The basin map is used to
develop a basin model. Hydrologic elements (subbasin, reach,
reservoir, etc.) are added from the toolbar and connected to represent
the physical drainage network of the study area. Background maps
can be imported to help visualize the watershed. The Castro 1 basin
map is shown in Figure 3.
27
Chapter 3 Overview
Component Editor
When a component or sub-component in the Watershed Explorer or
hydrologic element in the basin map is active (use the mouse and click
on the component name in the Watershed Explorer or select the
hydrologic element in the basin map), a specific component editor will
open. Data for model components is entered in the Component Editor.
Required data is indicated with a red asterisk. For example, loss
parameter data for a subbasin element is entered in the Component
Editor for the subbasin. The Component Editor for the Castro 1 basin
model is shown in Figure 4.
Figure 3. Basin map for a basin model named Castro 2. The basin map is
shown in the Desktop area of the program screen.
28
Chapter 3 Overview
Message Log
Notes, warnings, and errors are shown in the Message Log, as seen in
Figure 5. These messages are useful for identifying why a simulation
run failed or why a requested action, like opening a project, was not
completed.
Program Settings
Additional program settings are available. Click the Tools menu and
select the Program Settings command to access the Program
Settings window. This window allows you to change any of the
optional settings as shown in Figure 6.
The "General" tab contains a number of settings for the general
behavior of the program. This includes details such as how much
warning you want before making irreversible changes in modeling
components, where to display global editors, a display option for the
decimal separator, language for displaying calendar dates, and a
display option for latitudes and longitudes.
29
Chapter 3 Overview
The "Basin Map" tab contains settings for the basin map window. The
settings include the option to automatically recenter the map when
selecting an element. The selection highlight color and the grid line
color can also be selected.
The "Defaults" tab contains settings for the default unit system,
subbasin component methods, and reach component methods. These
are the settings that will be used when new components are created.
You can always change the setting after the component is created.
The "Compute" tab contains settings for closing the compute progress
if the compute is successful and other details of how a compute is
performed. Importantly, the buffer length controls how much memory
is required during a compute. When the length is set to a large value,
the compute is optimized for speed but may exhaust available
memory. When the length is set to a small value, very large computes
(large number of elements or very long time windows) can be
completed but the speed may be slower.
The "Results" tab contains a number of settings for managing
simulations and the results that are produced. Control is provided for
30
Chapter 3 Overview
supplemental results that give extra information but are not critical to
the simulation, output of gridded results when using the ModClark
transform method, and the location to display results. Finally, you can
specify the number of decimal digits to use when displaying different
categories of results.
The "Messages" tab contains settings for the error, warning, and note
messages. You can select whether a sound is played when a
particular type of message occurs. You can also choose the color that
is used to print the message in the Message Log.
Data Conventions
The program must manage many different kinds of data that are part of
a wide variety of components. Part of the management is the saving
and display of the data entered by the user. The remainder of the
management is the use of the data during a compute.
Saving Properties
The program uses the concept of the current component to manage
saving updates or changes to properties. The current component is
the item in the project that is currently selected in the Watershed
Explorer. For example, if you click on a subbasin icon in the
Watershed Explorer it will become highlighted, and its editor will be
shown in the Component Editor. The subbasin becomes the current
component at that moment. You may make changes to the properties
of the subbasin on any of the tabs in the Component Editor. When you
switch between tabs, any changes are automatically updated in the
subbasin. The changes are also updated when you leave the tabs and
click anywhere else in the program interface. Even though the updates
are made in the subbasin properties, the changes are not saved to
disk. Updates and changes are only saved to disk when the project is
saved. You can save the project two different ways.
The first way to save all pending changes is from the File menu. Click
on the File menu and then select the Save command. All components
currently open in the project will be saved to disk, including the current
basin model with its hydrologic elements, current meteorologic model,
current control specifications, and all other project components.
The second way to save all pending changes is from the toolbar. Click
the Save Current Project button. All components in the project will be
saved to disk.
31
Chapter 3 Overview
Number Formatting
Each country of the world has socially accepted conventions for
formatting numbers. It would be difficult for the program to support
every single convention in use around the world. However, one
important factor is the symbol used for the decimal separator. For
example, in the United States a number would typically be formatted
as follows:
12.34
32
Chapter 3 Overview
Units Conversion
Almost all initial conditions and parameter data for the various methods
included in the program require units. For example, hydraulic
conductivity has units of either millimeters per hour (mm/hr) or inches
per hour (in/hr) depending on the unit system of the basin model. The
unit system is one of the properties of the basin model and
meteorologic model. If the system international unit system is selected
for a basin model, then the hydraulic conductivity should be entered in
mm/hr. However, the conductivity should be entered in in/hr if the U.S.
customary unit system is selected. The units of an initial condition or
parameter are shown in parenthesis after the label.
The value of initial conditions and parameter data are automatically
converted when the unit system of a basin model or meteorologic
model is changed. For example, suppose a basin model used the U.S.
customary unit system and a hydraulic conductivity was entered as
0.23 in/hr. If the unit system were changed to system international, the
conductivity would be automatically converted to 5.84 mm/hr. The
conversions are performed according to standards specified by the
National Institutes of Standards and Technology (Taylor 1995).
Time-series data, paired data, and grid data components each have
their own unit system based on the units of the data. The unit system
is determined automatically from the units of the data. For example, a
discharge time-series gage with units of cubic meters per second
(M3/S) will be in the system international unit system. The units are
selected by the user for manual entry data, but are read automatically
from the record header for external DSS data. Data is automatically
converted to the correct unit system during a compute.
33
Chapter 3 Overview
Interpolation
Time-series data and gridset data are usually defined with a fixed time
interval, though some data may be defined on an irregular basis. All of
the different types of simulations happen with a fixed time interval, as
specified in the control specifications. When the time interval of the
time-series or gridset data does not match the time interval of the
compute, the data is automatically interpolated. A linear interpolation
in time is used.
Paired data components use a limited number of points to represent a
curve, such as a storage-discharge curve. However, the curve
represents continuous data. Linear interpolation is used when a
dependent value is required for an independent value between two
specified values. The interpolation is performed between the
dependent values corresponding to the closest available independent
values on each side of the requested value. Some paired data
components use irregularly spaced values to represent an annual
pattern, such as a groundmelt pattern. Linear interpolation in time is
used on these components.
Application Steps
The program is designed with reusable data sets that can be
independently developed. However, some data sets depend on others
for important definitions. For example, gages must be created before
they can be used in basin or meteorologic models. Consequently,
there is a necessary sequence to successfully obtain results. The
remainder of this chapter provides an overview of the best procedure
for obtaining computation results.
34
Chapter 3 Overview
Precipitation Altitude
Radiation Concentration
Windspeed Percent
Humidity Sunshine
All of the different types of paired data:
Elevation-area Duration-precipitation
Elevation-discharge ATI-meltrate
Inflow-diversion ATI-coldrate
Cumulative probability
35
Chapter 3 Overview
Elevation
Open a component manager to add shared data to a project. Go to
the Components menu and select Time-Series Data Manager,
Paired Data Manager, or Grid Data Manager command. Each one of
these component managers contains a menu for selecting the type of
data to create or manage. The Paired Data Manager with the
Storage-Discharge data type selected is shown in Figure 7. Once the
data type is selected, you can use the buttons on the right side of the
component manager to add a New, Copy, Rename, and Delete a data
type. There is also a button that opens an editor where you can
change the description of the selected component. In the case of time-
series data, the manager contains two extra buttons to add or delete
time windows. A time window is needed for entering or viewing time-
series data.
36
Chapter 3 Overview
37
Chapter 3 Overview
Hydrologic Description
Element
Subbasin The subbasin is used to represent the physical
watershed. Given precipitation, outflow from the
subbasin element is calculated by subtracting
precipitation losses, calculating surface runoff,
and adding baseflow.
Reach The reach is used to convey streamflow in the
basin model. Inflow to the reach can come from
one or many upstream elements. Outflow from
the reach is calculated by accounting for
translation and attenuation. Channel losses can
optionally be included in the routing.
Junction The junction is used to combine streamflow from
elements located upstream of the junction.
Inflow to the junction can come from one or many
upstream elements. Outflow is calculated by
summing all inflows.
Source The source element is used to introduce flow into
the basin model. The source element has no
inflow. Outflow from the source element is
defined by the user.
Sink The sink is used to represent the outlet of the
physical watershed. Inflow to the sink can come
from one or many upstream elements. There is
no outflow from the sink.
Reservoir The reservoir is used to model the detention and
attenuation of a hydrograph caused by a
reservoir or detention pond. Inflow to the
reservoir element can come from one or many
upstream elements. Outflow from the reservoir
can be calculated using one of three routing
methods.
Diversion The diversion is used for modeling streamflow
leaving the main channel. Inflow to the diversion
can come from one or many upstream elements.
Outflow from the diversion element consists of
diverted flow and non-diverted flow. Diverted
flow is calculated using input from the user. Both
diverted and non-diverted flows can be
connected to hydrologic elements downstream of
the diversion element.
38
Chapter 3 Overview
Figure 8. Reach component editor including data for routing and loss/gain
methods. The downstream element is also shown.
Figure 9. Global editor for the initial and constant loss method. Only the
elements selected in the basin map are shown initially. You can
switch between viewing the elements in the initial selection and
viewing all elements in the basin model that use the selected
method. You can also control how the elements are sorted.
39
Chapter 3 Overview
40
Chapter 3 Overview
41
Chapter 3 Overview
Figure 10. Control specifications component editor. The starting date and
time are required, along with the ending date and time.
42
Chapter 3 Overview
Results can be accessed from the basin map and the Watershed
Explorer, "Results" tab. Results are available as long as a simulation
run has been successfully computed and no edits have been made
after the compute to any component used by the simulation run. For
example, if the time of concentration parameter was changed for a
subbasin element after the simulation run was computed, then results
are no longer available for that hydrologic element and any elements
downstream of it in the basin model. The simulation run must be
recomputed for results to become available.
The simulation must be selected (from the Compute menu or
Watershed Explorer) before results can be accessed from the basin
map. After the simulation run is selected, select the hydrologic
element where you want to view results. While the mouse is located
on top of the element icon, click the right mouse button. In the menu
that opens, select the View Results option. Three result types are
available: Graph, Summary Table, and Time-Series Table (Figure
12). These results can also be accessed through the toolbar and the
Results menu. A hydrologic element must be selected before the
toolbar buttons and options from the Results menu are active. A global
summary table is available from the toolbar and Results menu. The
global summary table contains peak flows and time of peak flows for
each hydrologic element in the basin model.
43
Chapter 3 Overview
Figure 12. Accessing results for the current simulation run using the basin
map.
44
Chapter 3 Overview
Figure 13. Selecting simulation run results from the Watershed Explorer.
45
Chapter 3 Overview
Figure 14. Graph comparing outflow in the same element in two different
simulations runs. The correct time-series were selected in the
Watershed Explorer and then the graph toolbar button was
pressed.
46
Chapter 3 Overview
References
Taylor, B. 1995. Special Publication 811: Guide for the Use of the
International System of Units (SI). United States Department of
Commerce, National Institute of Standards and Technology.
Gaithersburg, MD.
47
Chapter 3 Overview
48
Chapter 4 Projects and Control Specifications
CHAPTER 4
Projects
A project represents all of the input data and simulation output
necessary to answer an engineering question. Possible questions
could include the following:
How will the average in-stream flow depth for the month of July
change after the operation schedule of an irrigation diversion is
modified?
49
Chapter 4 Projects and Control Specifications
Figure 15. Creating a new project. The location may be changed to suit
your preferences for storing and organizing projects.
Opening a Project
Open a project by selecting the File Open menu option. After this
option is selected, the Open an Existing Project window will open as
shown in Figure 16. This window contains all projects in the watershed
list. The watershed list contains all projects previously opened by the
program. Open a project by clicking the project name and pressing the
Open button. You may also double-click a project name to open it.
You cannot press the Open button if no project is selected. If the
project you want to open is not in the list, press the Browse button.
The Select Project File window opens that lets you navigate to the
directory containing the desired project (Figure 17). If you change your
mind and do not want to open a project, press the Cancel button or the
X button in the upper right to close the window. The Open an Existing
Project window can also be opened by pressing the Open a Project
button on the toolbar.
50
Chapter 4 Projects and Control Specifications
Figure 16. Opening a project. Previously opened projects are kept on a list
for easy reference.
Figure 17. Browsing to open a project that has not been opened previously.
Select the project file with extension HMS to open that project.
Copying a Project
Copy a project by selecting the File Save As menu option. After
this option is selected, the Save Current Project As window will open
where you can name, select the location on your computer or network
server to store the project, and describe the project as shown in Figure
51
Chapter 4 Projects and Control Specifications
You are not required to use the default directory and may copy a
project to the local computer or a network resource anywhere your
security limitations permit. All project files will be copied to a folder
(the folder name is the name of the project) in the location entered. If
the project description is long you can use the button to the right of the
description field to open an editor. Press the Copy button when you
are satisfied with the name, location, and description. You cannot
press the Copy button if no name or location is specified. If you
change your mind and do not want to create a copy, press the Cancel
button or the X button in the upper right to close the window.
There are two options in this window that need extra attention. The
first option, Copy external DSS data (time-series, paired, grid data),
will copy all external DSS records defined in the project and place the
copied records into DSS files in the new project folder. Time-series,
paired, and grid data are copied into separate DSS files. All
references to these DSS records are automatically updated. The
second option, Copy basin files (grid cell, maps), will copy all grid
cell files and background map layer files used by all basin models in
the project and paste them into the new project folder. The
background map layer files are placed in a folder named "Maps." All
references to grid cell and map layer files are automatically updated.
Renaming a Project
Rename a project by selecting the File Rename menu option. After
this option is selected, the Rename Current Project window opens
(Figure 19). The window shows the current project name and
description and contains a text box for you to enter the new project
52
Chapter 4 Projects and Control Specifications
name and description. If the project description is long you can use
the button to the right of the description field to open an editor. Press
the Rename button when you are satisfied with the name and
description. You cannot press the Rename button if no name is
specified. If you change your mind and do not want to rename the
project, press the Cancel button or the X button in the upper right to
close the window.
Deleting a Project
Delete a project by selecting the File Delete menu item. The Delete
Current Project window will open (Figure 20). This window contains
the project name and description and three choices for deleting the
project. The first option, Watershed entry in watershed list only (all
files remain intact), will only delete the project from the watershed list.
The watershed list contains all projects previously opened by the
program. No files are deleted when this option is selected. The
second option, Watershed entry and all watershed files (external
files remain intact), will delete the project from the project list and
delete all project files. Project files include *.hms, *.basin, etc. (refer to
Table 12 for a complete list of project files). No external files are
deleted when this option is selected. External files include DSS, grid
cell, and background map layer files. The third option, Watershed
entry, watershed files, and all external files, will delete the project
from the watershed list, all files created by the program for this project,
and all external files referenced by the project (DSS, grid cell, and
background map layer files). It is very important to keep in mind that
the third option will delete external files which could be used by other
projects. After the appropriate option is selected, press the Delete
button. Once a project has been deleted it cannot be retrieved or
undeleted. If you change your mind and do not want to delete the
53
Chapter 4 Projects and Control Specifications
project, press the Cancel button or the X button in the upper right to
close the window.
Project Properties
Basic properties and settings for the project are provided in the
Component Editor for the project. Access the Component Editor on
any tab of the Watershed Explorer. The highest level folder is labeled
with the project name; click on the project name to display the
Component Editor (Figure 21).
Figure 21. Project component editor accessed by clicking on the top node
of the Watershed Explorer. The top node is labeled with the
project name.
54
Chapter 4 Projects and Control Specifications
If the description is short, you can type it directly into the description
field. For long descriptions, you can press the editor button to the right
of the description field. The description editor makes it easier to enter
and edit long descriptions.
Output results are written to the project DSS file by default. However,
you can change the file where results are written. Results will not be
available after changing the output DSS file until the component is
recomputed. If you wish, you can press the folder button to the right of
the file field to open a file browser for locating the file. The file browser
is set to locate files with the DSS extension. If you change your mind
about searching for a different output file, you can press the Cancel
button or the X button in the upper right of the browser window. Press
the Select button after you locate and highlight the desired file. It is
possible to have more than one project write output results to the same
DSS file. However, the separate projects are not synchronized and
extreme care must be taken in naming components in order to avoid
conflicts.
55
Chapter 4 Projects and Control Specifications
Table 12. Files created by the program for storing component data.
Sample names have been used assuming a project named
Castro, that contains a basin model named North Branch, a
meteorologic model named Historic Calib, and a control
specifications named Oct 1977. It is also assumed that all
features of the program are used including an optimization trial
named Est Baseflow, and a saved state named Mid Oct 77.
Filename Description
Castro.access Control file; contains no data.
Castro.dsc Catalog of the project DSS file to speed finding
specific records.
Castro.dss Project DSS file containing manual entry time-
series gages and paired data. Default file for
output results.
Castro.gage List of all time-series gages.
Castro.grid List of all grid data.
Castro.hms Project definition including lists of basin models,
meteorologic models, and control specifications.
Castro.log All messages are recorded in the log file except
for messages generated during a compute. Each
compute component has its own log file.
Castro.nals List of all depth-area analyses.
Castro.out List of all data read or written from DSS files.
Castro.pdata Listing of all paired data.
Castro.run List of simulation runs with their properties and
the time of last compute.
North_Branch.basin Configuration for element properties, channel
network, zones, and map layers. One per basin.
Historic_Calib.met Configuration for precipitation, evapotranspiration,
and snowmelt data. One per meteorologic model.
Oct_1977.control Configuration for start time, end time, and time
interval. One per control specifications.
Castro.trial List of all the optimization trials.
Est_Baseflow.optim Optimization trial properties including search
method, objective function, and parameters.
Castro.fcstList list of all forecast alternatives.
DryStart.forecast Forecast alternatives properties including start
and end times, blending controls, and parameter
adjustments. One per forecast alternative
Castro.stateIndex List of all the saved states.
Mid_Oct_77.state Saved state information for all the hydrologic
elements in a specific basin model.
56
Chapter 4 Projects and Control Specifications
Computed Results
By default, all computed results are stored in the project DSS file.
However, the user has the option of changing the file where computed
results are stored. Any manually entered time-series or paired data
will continue to be stored in the project DSS file regardless of the file
used for computed results. If the computed results are stored in a
different file, that file can be safely deleted without affecting the
program. However, it is not possible to view results after deleting the
file where they are stored.
The result from each element is stored in a separate record. Some
elements compute different types of results; each result is stored in a
separate record. The record is identified with a pathname. Each
57
Chapter 4 Projects and Control Specifications
record pathname contains six parts called the A-part, B-part, C-part, D-
part, E-part, and F-part. The pathname parts are separated with a
slash and may contain spaces. The complete pathname, including
slashes and spaces, can total up to 256 uppercase characters. The
following is an example of a pathname for the computed flow at a
hydrologic element named "Sand Cr" in simulation run "Plan 3A":
//SAND CR/FLOW/01JAN1985/1HOUR/RUN:PLAN 3A/
Security Limitations
The program can create a project on the local computer or on any
accessible network device. Creating a project requires the user to
have read and write permission for the folder that will contain the new
project. Usually the system administrator determines where a user has
permission to read and write. Depending on the security settings
assigned to the user by the system administrator, the program may not
be able to create a project in some folders. As a user, be sure you
understand where you have permission to create new projects. No
other privileges beyond read and write permission are required to use
the program; it is fully compatible with operating in a so-called reduced
privilege environment.
The program automatically checks all of the project files for read and
write permission every time the project is opened. If any of the files
are read-only, then the program will not be able to open the project and
an error message will be displayed. If a project was previously
58
Chapter 4 Projects and Control Specifications
Control Specifications
Control specifications are one of the main components in a project,
even though they do not contain much parameter data. Their principle
purpose is to control when simulations start and stop, and what time
interval is used in the simulation.
59
Chapter 4 Projects and Control Specifications
the new control specifications. You cannot press the Create button if
no name is specified for the new specifications. If you change your
mind and do not want to create a new control specifications, press the
Cancel button or the X button in the upper right to return to the Control
Specifications Manager window.
60
Chapter 4 Projects and Control Specifications
61
Chapter 4 Projects and Control Specifications
Figure 25. Renaming a control specifications. This was accessed from the
Control Specifications Manager.
62
Chapter 4 Projects and Control Specifications
63
Chapter 4 Projects and Control Specifications
64
Chapter 4 Projects and Control Specifications
Time Window
Each control specifications sets the time window over which a
simulation will be performed. The window is specified using a separate
start date, start time, end date, and end time. There is no limit on the
length of a time window, or the number of simulation time steps it can
contain. The program contains simulation methods suitable for both
event and continuous simulation. Whether a particular simulation is
considered event or continuous depends on the length of the time
window set in the control specifications and the methods chosen in the
basin model.
The program is capable of processing dates from 1583 AD through
4000 AD. All dates are processed according to the rules specified in
the Gregorian calendar. The format for specifying a date is to use two
digits for the day, followed by the three-letter month abbreviation, and
finally the four digit year. Two digit years are never used for entering
or displaying dates. For example, the date February 8, 1936 should be
entered as follows:
08Feb1936
It is very important to use the correct format or the date you enter may
be incorrectly interpreted. If the program is not able to interpret a date,
the entry field will become blank. The same format is used for both
start and end dates, and for dates throughout the program.
The program processes times assuming an arbitrary local time zone
that does not observe summer time (daylight savings in the United
States). It uses 24-hour clock time instead of AM or PM notation.
Time windows can only be entered with minute resolution. Times may
range from 00:00 at the beginning of a day to 23:59 at the end. If a
time of 24:00 is entered, it is automatically converted to 00:00 on the
following day. For example, the time of 2:15:00 PM should be entered
as follows:
14:15
It is very important to use the correct format, including the colon, or the
time may be incorrectly interpreted. The same format is used for start
and end times, and for times throughout the program.
Enter or edit the time window in the Component Editor for the control
specifications. Access the editor from the Watershed Explorer on the
"Components" tab by clicking on the desired control specifications icon
(Figure 29). Type the date and time information in the appropriate
fields.
65
Chapter 4 Projects and Control Specifications
Time Interval
Each control specifications includes the time interval that will be used
to perform computations during a simulation. The same interval will be
used when viewing time-series results from the simulation. In some
cases, computations may be carried out at a shorter time interval and
results interpolated to the interval specified in the control
specifications. Some methods in the basin model have parameters
that are sensitive to time interval. Those parameters have to be
estimated with knowledge of the time interval in the control
specifications.
Time-series gage data and grid data are interpolated to the time
interval during a simulation. The original data stored for the gage or
grid is not altered; the available data is interpolated as part of the
simulation process and is not retained. The interpolation is performed
linearly in time between the available data points.
Specify the time interval in the Component Editor for the control
specifications (Figure 29). Access the editor from the Watershed
Explorer on the Components tab by clicking on the desired control
specifications icon. Select the desired time interval from the list of
available choices. The possible choices range from 1 minute to 1 day.
It is important that any minutes included in a start or end time be an
integer multiple of the time interval. For example, if the time interval
were chosen to be 10 minutes, then the start or end time could be
10:00, 10:10, 10:20 or other multiples of 10 minutes. For the same
time interval, the start or end time could not be 10:07 since that is not
an integer multiple of the time interval.
66
Chapter 4 Projects and Control Specifications
67
Chapter 4 Projects and Control Specifications
results are stored assuming the local time zone. However, gridded
precipitation data is typically published with a time zone of GMT. The
gridded precipitation method includes an option to shift the data to the
local time zone so that it may be aligned with other input data during
the simulation. The gridded output results are normally stored in the
local time zone. By using the grid time shift, the gridded output results
are shifted after the simulation is complete so that the gridded output
will align with the original gridded precipitation data.
This feature is considered an interim solution to processing input data
with different time zone definitions. Eventually, a complete solution will
be implemented where all time-series and gridded input data are
labeled with a time zone and the output time-series and gridded
simulation results will also be labeled with a time zone.
68
Chapter 4 Projects and Control Specifications
correct HEC-1 file. The browser will allow you to select HEC-1 files,
which usually end with the DAT or HC1 extensions. Once you have
located and selected the desired HEC-1 file, press the Select button.
If you change your mind about selecting a file, you can press the
Cancel button or the X button in the upper right of the browser to
return to the Import HEC-1 File window without selecting a file. If you
change your mind about importing a HEC-1 file, you can press the
Cancel button or the X button in the upper right of the Import HEC-1
File window.
Figure 30. Selecting a HEC-1 file to import. The default names could be
changed to better reflect the components.
After you specify the HEC-1 file to import, you must enter the name of
the basin model, meteorologic model, and control specifications. You
can use the default names that are provided or you can enter your own
names. You can always rename the components at a later time. If the
HEC-1 file includes ZR records for loading time-series data from a
HEC-DSS file, then you should specify the DSS file that contains the
data. Without specifying the file, you will have to manually connect to
the time-series data. Press the Import button when you have finished
specifying the HEC-1 file name to import, and the names of the three
components. You cannot import a HEC-1 file unless the file name and
all three component names are specified. When you press the Import
button, the import process will read the HEC-1 file and create the
appropriate parts of the three components.
Unsupported Features
Most of the features in the HEC-1 program are also available in this
program. Some features related to economic estimates were not
incorporated into this program because they are inconsistent with
69
Chapter 4 Projects and Control Specifications
70
Chapter 4 Projects and Control Specifications
References
Hydrologic Engineering Center. June 1998. HEC-1 Flood Hydrograph
Package: User's Manual. U.S. Army Corps of Engineers, Davis, CA.
Hydrologic Engineering Center. October 1994. HEC-DSS User's
Guide and Utility Manuals: User's Manual. U.S. Army Corps of
Engineers, Davis, CA.
Hydrologic Engineering Center. March 1984. HMR52 Probable
Maximum Precipitation (Eastern United States): User's Manual. U.S.
Army Corps of Engineers, Davis, CA.
71
Chapter 4 Projects and Control Specifications
72
Chapter 5 Shared Component Data
CHAPTER 5
Time-Series Data
Hydrologic models often require time-series of precipitation data for
estimating basin-average rainfall. A time-series of flow data, often
called observed flow or observed discharge, is helpful for calibrating a
model and is required for optimization. Other kinds of time-series data
are used as well. Time-series data is stored in a project as a gage.
The program separates different types of data with different gage
types. Gage data only has to be entered one time. The gages are part
of the project and can be shared by multiple basin or meteorologic
models.
73
Chapter 5 Shared Component Data
Figure 32. Creating a new discharge gage after pressing the New button
in the Time-Series Data Manager.
the right of the time-series data list can be used to manage existing
data or create new data. To create a new time-series gage, press the
New button. After you press the button a window will open (Figure
32) where you can name and describe the new gage. A default name
is provided for the new gage; you can use the default or replace it with
your own choice. A description can also be entered. If the description
is long, you can press the button to the right of the description field to
open an editor.
The editor makes it easy to enter and edit long descriptions. When
you are satisfied with the name and description, press the Create
button to finish the process of creating the new time-series gage. You
cannot press the Create button if no name is specified. If you change
your mind and do not want to create a new time-series gage, press the
Cancel button or the X button in the upper right to return to the time-
series data manager.
74
Chapter 5 Shared Component Data
Copying a Gage
There are two ways to copy a time-series gage. Both methods for
copying a gage create an exact duplicate with a different name. Once
the copy has been made it is independent of the original and they do
not interact.
The first way to create a copy is to use the time-series data manager,
which is accessed from the Components menu. First, select the data
type of the time-series gage you want to copy from the Data Type
menu. Then, select the time-series gage you want to copy by clicking
on it in the list of current time-series gages. The selected gage is
highlighted after you select it. After you select a gage you can press
the Copy button on the right side of the window. A new window will
open where you can name and describe the copy that will be created
(Figure 33). A default name is provided for the copy; you can use the
default or replace it with your own choice.
Figure 33. Creating a copy of a gage after pressing the Copy button in the
Time-Series Data Manager.
75
Chapter 5 Shared Component Data
A description can also be entered; if it is long you can use the button to
the right of the description field to open an editor. When you are
satisfied with the name and description, press the Copy button to finish
the process of copying the selected time-series gage. You cannot
press the Copy button if no name is specified. If you change your
mind and do not want to copy the selected gage, press the Cancel
button or the X button in the upper right to return to the time-series
data manager.
The second way to create a copy is from the Watershed Explorer, on
the "Components" tab. Move the mouse over the time-series
component you wish to copy, then press the right mouse button
(Figure 34). A context menu is displayed that contains several choices
including copy. Click the Create Copymenu option. A new window
will open where you can name and describe the copy that will be
created. A default name is provided for the copy; you can use the
default or replace it with your own choice. A description can also be
entered; if it is long you can use the button to the right of the
description field to open an editor. When you are satisfied with the
name and description, press the Copy button to finish the process of
copying the selected time-series gage. You cannot press the Copy
button if no name is specified. If you change your mind and do not
want to copy the gage, press the Cancel button or the X button in the
upper right of the window to return to the Watershed Explorer.
Renaming a Gage
There are two ways to rename a time-series gage. Both methods for
renaming a gage changes its name and then all references to the old
name are automatically updated to the new name.
The first way to perform a rename is to use the time-series data
manager, which you can access from the Components menu. First,
select the data type of the time-series gage you want to rename from
the Data Type menu. Then, select the time-series gage you want to
rename by clicking on it in the list of current time-series gages. The
selected gage is highlighted after you select it. After you select a gage
you can press the Rename button on the right side of the window. A
window will open where you can provide the new name (Figure 35).
You can also change the description at the same time. If the new
description will be long, you can
76
Chapter 5 Shared Component Data
Figure 35. Renaming a gage after pressing the Rename button in the
Time-Series Data Manager.
use the button to the right of the description field to open an editor.
When you are satisfied with the name and description, press the
Rename button to finish the process of renaming the selected time-
series gage. You cannot press the Rename button if no name is
specified. If you change your mind and do not want to rename the
selected gage, press the Cancel button or the X button in the upper
right of the window to return to the time-series data manager.
The second way to rename is from the Watershed Explorer, on the
"Components" tab. Move the mouse over the time-series component
you wish to rename, then press the right mouse button (Figure 36). A
context menu is displayed; select the Rename command and the
highlighted name will change to editing mode. You can then move the
cursor with the arrow keys on the keyboard or by clicking with the
mouse. You can also use the mouse to select some or all of the name.
Change the name by typing with the keyboard. When you have
finished changing the name, press the Enter key to finalize your
77
Chapter 5 Shared Component Data
choice. You can also finalize your choice by clicking elsewhere in the
Watershed Explorer. If you change your mind while in editing mode
and do not want to rename the selected gage, press the Escape key.
Deleting a Gage
There are two ways to delete a time-series gage. Both methods for
deleting a gage will remove it from the project and then automatically
update all references to that gage. Once a gage has been deleted it
cannot be retrieved or undeleted. Any references to the deleted gage
will switch to using no gage, which is usually not a valid choice during
a simulation. At a later time you will have to go to those components
and manually select a different gage.
The first way to perform a deletion is to use the time-series data
manager, which you can access from the Components menu. First,
select the data type of the time-series gage you want to delete from the
Data Type menu. Then, select the time-series gage you want to
delete by clicking on it in the list of current time-series gages. The
selected gage is highlighted after you select it. After you select a gage
you can press the Delete button on the right side of the window (Figure
37). A window will open where you must confirm that you want to
delete the selected gage. Press the OK button to delete the gage. If
you change your mind and do not want to delete the selected gage,
press the Cancel button or the X button in the upper right to return to
the time-series data manager.
Figure 37. Preparing to delete a gage from the Time-Series Data Manager.
78
Chapter 5 Shared Component Data
Time Windows
Time windows are used to separate the time-series data into
manageable sections. You may choose to have a separate time
window for each event. Alternately you may have several time
windows for a continuous record to break it into months or years. You
may choose to have a combination of time window types and they may
overlap. All time windows use the same data units, time interval, and
other properties discussed in the following sections.
The program is capable of processing dates from 1583 AD through
4000 AD. All dates are processed according to the rules specified in
the Gregorian calendar. The format for specifying a date is to use two
digits for the day, followed by the three-letter month abbreviation, and
finally the four digit year. Two digit years are never used for entering
or displaying dates. For example, the date October 2, 1955 should be
entered as follows:
02Oct1955
79
Chapter 5 Shared Component Data
It is very important to use the correct format or the date you enter may
be incorrectly interpreted. If the program is not able to interpret a date,
the entry field will become blank. The same format is used for both
start and end dates, and for dates throughout the program.
The program processes times assuming an arbitrary local time zone
that does not observe summer time (daylight savings in the United
States). It uses 24-hour clock time instead of AM or PM notation.
Time windows can only be entered with minute resolution. Times may
range from 00:00 at the beginning of a day to 23:59 at the end. If a
time of 24:00 is entered, it is automatically converted to 00:00 on the
following day. For example, the time of 6:20:00 PM should be entered
as follows:
18:20
It is very important to use the correct format, including the colon, or the
time may be incorrectly interpreted. The same format is used for start
and end times, and for times throughout the program.
There are two ways to create a new time window. The first way is from
the Time-Series Data Manager, accessed by clicking the Components
menu and then selecting the Time-Series Data Manager command.
Select the desired data type, then click on a time-series data
component in the list; the component will become highlighted. Press
the Add Window button to create a new time window. The Add Time-
Series Data Time Window window will open where you can enter the
start date and other information as shown in Figure 39. You can either
enter the information manually, or select a control specifications. If you
select a control specifications, the start and end time in that control
specifications will be used for the new time window. Press the Add
button to create the new time window. The window will remain open
for adding additional time windows. When you are finished, press the
Close button or the X button in the upper corner of the Add Time-
Series Data Time Window. The second way to create a new time
window is directly from the Watershed Explorer. Select a time-series
component by clicking on it or one of the existing time windows. Keep
the mouse over the gage or time window icon and click the right mouse
button. A context menu appears as shown in Figure 40; click the
Create Time Window command to create a new time window. The
same window shown in Figure 39 will open for creating a new time
window.
80
Chapter 5 Shared Component Data
Figure 39. Creating a new time window for a gage, beginning from the
Time-Series Data Manager.
Figure 40. Creating a new time window from the Watershed Explorer.
There are two ways to delete a time window. The first way is from the
Time-Series Data Manager, accessed by clicking the Components
menu and then selecting the Time-Series Data Manager command.
Select the desired data type, then click on a time-series gage in the
list; it will become highlighted. Press the Delete Window button to
delete a time window. The Delete Time-Series Data Time Window
window will open where you can select the window to delete (Figure
41). Click on the desired window and it will become highlighted. Press
the Delete button to delete the highlighted time window. If you change
your mind and do not want to delete a time window, press the X button
in the upper corner of the Delete Time-Series Data Time Window
window. The second way to delete a time window is directly from the
Watershed Explorer. Select a time window for a time-series gage; it
will become highlighted. Keep the mouse over the time window icon
and click the right mouse button. A context menu appears as shown in
81
Chapter 5 Shared Component Data
Figure 42; click the Delete Time Window command to delete the
selected time window.
Figure 41. Selecting a time window to delete from a gage after pressing the
Delete Window button in the Time-Series Data Manager.
Figure 42. Deleting a selected time window from a gage after choosing it in
the Watershed Explorer.
The user can change the start date, start time, end date, and end time
of an existing time window. Use the Watershed Explorer to select the
time window you wish to change. Click on the time window under the
correct time-series component. The component will become the
selected component and its data will be shown in the Component
Editor as seen in Figure 43; the "Time Window" tab is automatically
selected. Change the start date or other properties to the desired
values. Click on a different tab in the Component Editor or elsewhere
in the program interface to make the changes take affect.
82
Chapter 5 Shared Component Data
Figure 43. Using the time window component editor for a gage to view and
edit the start and end time for a window.
Data Source
The data source determines how the data for a time-series component
will be entered. There are three alternative methods for entering the
data:
1. The user may enter the time-series values by manually typing
each one. A variation of this method is using the clipboard to
copy and paste instead of manually typing. This first method is
known as "Manual Entry."
2. The user may load data into a Data Storage System (HEC-DSS)
file. This method can be facilitated using the program HEC-
DSSVue to create individual HEC-DSS files. The gage is linked
to the time-series values stored in the file. This second method
is known as "Single Record HEC-DSS."
3. The user may store alternative versions of the time-series
values in a HEC-DSS file. This method is used exclusively for
forecasting operations with the Corps Water Management
System (CWMS) and planning studies with the Watershed
Analysis Tool (HEC-WAT). This third method is known as
"Multiple Record HEC-DSS."
The selected data source will determine how the time-series values are
entered. The values may be manually typed or linked to a HEC-DSS
file. Compare Figure 44 and Figure 45 to see the differences in
specifyin the data sources.
83
Chapter 5 Shared Component Data
Figure 45. Component editor for a precipitation gage retrieving data from a
Data Storage System (HEC-DSS) file. Of the different data
types, only precipitation gages have latitude and longitude.
Data Units
The data units can only be selected for a manual entry time-series
gage (Figure 44); they are retrieved automatically for both HEC-DSS
options (Figure 45). Most types of time-series data have only two
options for units; one for the system international unit system and one
for the U.S. customary unit system. For example, discharge gages can
use cubic meters per second (M3/S) or cubic feet per second (CFS).
84
Chapter 5 Shared Component Data
Time Interval
The time interval can only be selected for a manual entry time-series
gage (Figure 44); it is retrieved automatically for both HEC-DSS
options (Figure 45). An interval must be selected from the available
choices that range from 1 minute to 1 day. All time windows defined
for a time-series component must use the same time interval.
Generally you should choose the time interval before entering any data
for the component. However, if you change the time interval after
entering data, the data will be adjusted to the new time interval. When
the time interval is made shorter, the data for each time window will be
adjusted so that it still begins at the start of the time window. The data
will have the new, shorter time interval and there will be missing data
from the last specified value to the end of the time window. When the
time interval is made longer, the data for each time window will be
adjusted so that it still begins at the start of the time window. The data
will have the new, longer time interval and the end of the time window
will be advanced so that no data is lost.
Select the time interval for a time-series gage using the Component
Editor. Access the editor by selecting a time-series gage in the
Watershed Explorer. The "Time-Series Gage" tab in the Component
Editor will display the time interval if the manual entry option is
selected.
85
Chapter 5 Shared Component Data
Elevation
Certain atmospheric variables show a strong trend with elevation. Air
pressure and air temperature both decrease as the elevation
increases. Relative humidity increases as the elevation increases.
The physics of water and air can be used to develop relationships for
how these variables change with elevation. Therefore, calculating
atmospheric variables at a location in the watershed requires both the
elevation where the measurement was performed and the elevation
where the variable should be calculated. The elevation must be
entered for data types that require it. In general, the ground surface
elevation (above sea level) at the measurement site is used. The
value may be entered as either meters or feet, depending on the unit
system of the gage. The requirements for entering latitude and
longitude are shown in Table 14.
Reference Height
Some atmospheric variables change dramatically near the ground
surface. Windspeed is very close to zero at the ground surface and
increases as distance above the ground increases. Air temperature
and relative humidity change significantly with distance above a snow
covered surface. Therefore, calculating these atmospheric variables at
the ground surface requires information about the height above the
ground where the measurements were performed. While it is common
for meteorologic instruments to be mounted on a tower 10 meters
above the ground surface, other distances may be used at some
installations. The instrument reference height must be entered for data
types that require it. The value may be entered as either meters or
feet, depending on the unit system of the gage. The requirements for
entering latitude and longitude are shown in Table 14.
86
Chapter 5 Shared Component Data
87
Chapter 5 Shared Component Data
Manual Entry
Manual entry is the simpliest way to get data values for a time-series
gage. The values are manually typed for each date and time. All of
the data should be reviewed after entry to verify accurate typing. If no
value was recorded for a date and time, then the value should be left
blank. The program does not automatically fill in missing data. Data
recovery should be considered for extended periods of missing data.
Manual entry can also be used in combination with the clipboard to
transfer data values from another computer program. Time-series
values may be stored in a spreadsheet or other program. The values
can be selected in the other program and copied to the clipboard. The
values copied to the clipboard can be pasted in HEC-HMS. Using the
clipboard can save a significant amount of manual typing and reduce
errors when the time-series values are available in another program.
Manual entry time-series values are automatically stored in a HEC-
DSS file in the project directory. The program uses the name of the
gage to create a HEC-DSS file for storing the manual entry data.
Manual entry values can be entered and then edited as many times as
necessary to correct typing errors. Each time the data values are
modified manually by the user, the data is automatically updated in the
HEC-DSS file. One benefit of using a HEC-DSS file as an automatic
storage solution is that the file can be reduced for other modeling
projects.
88
Chapter 5 Shared Component Data
89
Chapter 5 Shared Component Data
match that choice will be shown in the selection table. If you make
choices in several filters, only pathnames that satisfy all of the choices
will be shown in the selection table.
The program observes a very strict set of rules for data type and units
within the record pathnames. Rules governing the C-part of the
pathname are also enforced. Data cannot be used unless is follows
the rules correctly; error messages will be generated if you attempt to
use an invalid C-part, data type, or units. The acceptable data types
for the different types of time-series data are shown in Table 15. The
correct unit labels are shown in Table 16.
Table 15. Internal DSS data type label for different types of time-series
data.
90
Chapter 5 Shared Component Data
91
Chapter 5 Shared Component Data
Table 16. Internal DSS units label for different types of time-series data.
92
Chapter 5 Shared Component Data
Table
The data for the current time window is shown in tabular form on the
"Table" tab of the Component Editor (Figure 46). If you select a time-
series gage in the Watershed Explorer, only the tab for the "Time-
Series Gage" is shown in the Component Editor. If you select a time
window under a time-series gage in the Watershed Explorer, the
"Table" tab will be added to the Component Editor. Data in the table
can be edited if the gage uses manual entry. Data is not editable if the
93
Chapter 5 Shared Component Data
gage retrieves data from a DSS file. Furthermore, if the gage uses
DSS data and no time-series data is available for the specified time
window, then the table will not contain any data.
Figure 46. Manually entering data for a temperature gage. The fill
command will be used to linearly interpolate between two known
temperatures.
You can enter all of the data for each time window one value at a time
in the table. However, there are tools to help you enter the data
quickly. The table includes support for the clipboard. This means you
can copy data stored in a spreadsheet or other file and then paste it
into the table. You can also use the fill tool to enter or adjust data
values in the table. Select the cells in the table you wish to fill and click
the right mouse button. A context menu is displayed; select the Fill
command. The Fill Table Options window opens for you to control the
process of filling and adjusting cell values. Options include
interpolating the values between the first and last cell in the selection,
several choices for interpolating and replacing missing values, copying
the first selected cell value to all other selected cells, adding a constant
value to all selected cells, and multiplying the selected cell values by a
constant. Press the OK button to apply your choice, or the Cancel
button to return to the table without making any changes.
Graph
The data for the current time window is shown in graphical form on the
"Graph" tab of the Component Editor (Figure 47). If you select a time-
series gage in the Watershed Explorer, only the tab for the "Time-
Series Gage" is shown in the Component Editor. If you select a time
94
Chapter 5 Shared Component Data
Figure 47. Viewing data for a discharge gage connected to a HEC-DSS file.
Paired Data
Hydrologic models often require the use of paired data to describe
inputs that are functional in form. Functional data defines a dependant
variable in terms of an independent variable. For most cases, the
function must be monotonically increasing which means it only
increases and never decreases. Examples of paired data include unit
hydrographs and stage-discharge curves. The program separates
different types of paired data with different data types. Paired data
only has to be entered one time. The data are part of the project and
can be shared by multiple basin or meteorologic models.
95
Chapter 5 Shared Component Data
Reference source not found. for a complete list of paired data types.
When a data type is selected, the manager will show all paired data of
the same type. The buttons to the right of the paired data list can be
used to manage existing data or create new data. To create a paired
data curve, press the New button. After you press the button a
window will open (Figure 49) where you can name and describe the
new curve. A default name is provided for the new curve; you can use
the default or replace it with your own choice. A description can also
be entered. If the description is long, you can press the button to the
right of the description field to open an editor. The editor makes it
easier to enter and edit long descriptions. When you are satisfied with
the name and description, press the Create button to finish the
process of creating the new paired data curve. You cannot press the
Create button if no name is specified. If you change your mind and do
not want to create a paired data curve, press the Cancel button or the
X button in the upper right to return to the paired data manager.
Figure 48. Paired Data Manager accessed from the Components menu.
Figure 49. Creating a new cross section after pressing the New button in
the Paired Data Manager.
96
Chapter 5 Shared Component Data
Copying a Curve
There are two ways to copy a paired data curve. Both methods for
copying a curve create an exact duplicate with a different name. Once
the copy has been made it is independent of the original and they do
not interact.
The first way to create a copy is to use the paired data manager, which
is accessed from the Components menu. First, select the data type of
paired data curve you want to copy from the Data Type menu. Then,
select the paired data curve you want to copy by clicking on it in the list
of current paired data curves. The selected curve is highlighted after
you select it. After you select a curve you can press the Copy
button on the right side of the window. A new window will open where
you can name and describe the copy that will be created (Figure 50).
A default name is provided for the copy; you can use the default or
replace it with your own choice. A description
Figure 50. Creating a copy of a gage after pressing the Copy button in the
Paired Data Manager.
can also be entered; if it is long you can use the button to the right of
the description field to open an editor. When you are satisfied with the
name and description, press the Copy button to finish the process of
copying the selected paired data curve. You cannot press the Copy
button if no name is specified. If you change your mind and do not
want to copy the selected curve, press the Cancel button or the X
button in the upper right to return to the paired data manager.
The second way to create a copy is from the Watershed Explorer, on
the Components tab. Move the mouse over the paired data
component you wish to copy, then press the right mouse button
(Figure 51). A context menu is displayed that contains several choices
including copy. Click the Create Copy menu option. A new window
will open where you can name and describe the copy that will be
created. A default name is provided for the copy; you can use the
default or replace it with your own choice. A description can also be
entered; if it is long you can use the button to the right of the
description field to open an editor. When you are satisfied with the
name and description, press the Copy button to finish the process of
97
Chapter 5 Shared Component Data
copying the selected paired data curve. You cannot press the Copy
button if no name is specified. If you change your mind and do not
want to copy the curve, press the Cancel button or the X button in the
upper right of the window to return to the Watershed Explorer.
Renaming a Curve
There are two ways to rename a paired data curve. Both methods for
renaming a curve changes its name and then all references to the old
name are automatically updated to the new name.
The first way to perform a rename is to use the paired data manager,
which you can access from the Components menu. First, select the
data type of paired data curve you want to rename from the Data Type
menu. Then, select the paired data curve you want to rename by
clicking on it in the list of current curves. The selected curve is
highlighted after you select it. After you select a curve you can press
the Rename button on the right side of the window. A window will
open where you can provide the new name (Figure 52). You can also
change the description at the same time. If the new description will be
long, you can use the button to the right of the field to open an editor.
When you are satisfied with the name and description, press the
Rename button to finish the process of renaming the selected paired
data curve. You cannot press the Rename button if no name is
specified. If you change your mind and do not want to rename the
selected curve, press the Cancel button or the X button in the upper
right of the window to return to the paired data manager.
98
Chapter 5 Shared Component Data
Figure 52. Renaming a cross section after pressing the Rename button in
the Paired Data Manager.
Deleting a Curve
There are two ways to delete a paired data. Both methods for deleting
a curve will remove it from the project and then automatically update all
references to that curve. Once a curve has been deleted it cannot be
99
Chapter 5 Shared Component Data
Figure 54. Preparing to delete a cross section from the Paired Data
Manager.
100
Chapter 5 Shared Component Data
Data Source
The data source determines how the data for a paired data will be
stored. Data may be entered manually or retrieved from the Data
Storage System (HEC-DSS). Manual entry means that you must enter
each of the values in the curve or pattern. Data can be retrieved from
a DSS file without having to enter the values, but the data must be
correctly loaded into the file. The type of properties you specify for the
paired data will be determined by the data source, as discussed in the
following sections. Compare Figure 56 and Figure 57 to see the
difference in entering data for manual or HEC-DSS data sources.
101
Chapter 5 Shared Component Data
Figure 56. Component editor for a unit hydrograph paired data curve with
manually-entered data. Only unit hydrographs have a time
interval.
Data Units
The data units can only be selected for a manual entry paired data
(Figure 56); they are retrieved automatically for the DSS option (Figure
57). Paired data have only two options for units; one for the system
international unit system and one for the U.S. customary unit system.
For example, a storage-discharge curve can use 1000 cubic meters
and cubic meters per second (1000M3 - M3/S) or acre feet and cubic
feet per second (ACFT - CFS). The units available in the "Units" field
will depend on the paired data type of the selected component.
Generally you should choose the data units before entering any data
for the component. However, if you change the units after entering
data, the data will be adjusted to the new units. There is no units
conversion during the adjustment. The values are all kept the same
but the assigned units are changed. This is helpful when the data is
entered without first checking to make sure the data units are in the
desired unit system.
Select the data units for a paired data component using the
Component Editor. Access the editor by selecting a paired data in the
102
Chapter 5 Shared Component Data
Watershed Explorer. The "Paired Data " tab in the Component Editor
will display the data units if the manual entry option is selected.
Time Intervals
Unit hydrographs are the only type of paired data that can have a time
interval as seen in Figure 56. The time interval can only be selected
for a manual entry unit hydrograph; it is retrieved automatically for the
DSS option. An interval must be selected from the available choices
that range from 1 minute to 24 hours. Generally you should choose
the time interval before entering data for the unit hydrograph.
However, if you change the time interval after entering data, the data
will be adjusted to the new time interval.
Select the time interval for a unit hydrograph paired data using the
Component Editor. Access the editor by selecting a unit hydrograph
paired data in the Watershed Explorer. The "Paired Data" tab in the
Component Editor will display the time interval if the manual entry
option is selected.
103
Chapter 5 Shared Component Data
part, E-part, and F-part. The pathname parts are separated with a
slash and may contain spaces. The complete pathname, including
slashes and spaces, can total up to 256 uppercase characters. The
following is an example of a storage-discharge curve pathname:
//BEAR CREEK/STORAGE-OUTFLOW///COMPUTED/
If you do not know the full pathname of the record you wish to use, you
can use the pathname browser to specify it. You must select a DSS
file first before the browser is available. Press the Select DSS
Pathname button to the right of the "DSS Pathname" field to open the
browser. The browser initially shows all of the records in the specified
DSS file, organized by pathname in the selection table. You can scroll
through the list and select a record pathname by clicking on it. Press
the Select button at the bottom of the browser to choose that record
and return to the Component Editor. If you change your mind and do
not want to select a record pathname, press the Cancel button or the
X button in the upper right of the Select Pathname From HEC-DSS File
window. You can reduce the number of record pathnames shown in
the selection table using the "Search by Parts" filters. A separate filter
selection is shown for each of the six pathname parts. By selecting a
choice for a filter, only pathnames that match that choice will be shown
in the selection table. If you make choices in several filters, only
pathnames that satisfy all of the choices will be shown in the selection
table.
The program observes a preferred order for paired data. For example,
storage-outflow is preferred to outflow-storage. However, the program
is capable of using data specified backwards from the preferred order.
The program does observe a very strict set of rules for data type and
units within the record pathnames. Data cannot be used unless is
follows the rules correctly; error messages will be generated if you
attempt to use an invalid data type or units. The units of both halves of
the paired data must be in the same unit system. The acceptable data
types for the different types of paired data are shown in Table 17. The
correct unit labels are shown in Error! Reference source not found..
104
Chapter 5 Shared Component Data
Table 17. Preferred data order and internal DSS data type label for
different types of paired data.
105
Chapter 5 Shared Component Data
Table 18. Internal DSS units for different types of paired data. The unit
hydrograph and patterns uses dates in the year 2000.
106
Chapter 5 Shared Component Data
Table
The data for the current paired data is shown in tabular form on the
"Table" tab of the Component Editor (Figure 58). Data in the table can
be edited if the paired data uses manual entry. Data is not editable if
the paired data retrieves data from a DSS file. Furthermore, if the
paired data uses DSS data and no data is available, then the table will
not contain any data.
107
Chapter 5 Shared Component Data
You can enter all of the data one value at a time in the table. However,
there are tools to help you enter the data quickly. The table includes
support for the clipboard. This means you can copy data stored in a
spreadsheet or other file and then paste it into the table. You can also
use the fill tool to enter or adjust data values in the table. Select the
cells in the table you wish to fill and click the right mouse button. A
context menu is displayed; select the Fill command. The Fill Table
Options window opens for you to control the process of filling and
adjusting cell values. Options include interpolating the values between
the first and last cell in the selection, several choices for interpolating
and replacing missing values, copying the first selected cell value to all
other selected cells, adding a constant value to all selected cells, and
multiplying the selected cell values by a constant. Press the OK button
to apply your choice, or the Cancel button to return to the table without
making any changes.
Graph
The data for the current paired data is shown in graphical form on the
"Graph" tab of the Component Editor (Figure 59). Data in the graph
cannot be edited regardless of whether the paired data uses manual
entry or retrieves data from an external DSS file. If no data is
available, then the graph will not contain any data.
108
Chapter 5 Shared Component Data
Grid Data
Some of the methods included in the program operate on a grid cell
basis. This means that parameters must be entered for each grid cell.
It also means that boundary conditions like precipitation must be
available for each grid cell. One example of grid data is the SCS curve
number grid, which is used as a parameter grid. Another example is
the temperature gridset, which is a time-series of grids. The program
separates different types of grid data with different data types. Grid
data only has to be entered one time. The data are part of the project
and can be shared by multiple basin or meteorologic models.
109
Chapter 5 Shared Component Data
long, you can press the button to the right of the description field to
open an editor. The editor makes it easier to enter and edit long
descriptions. When you are satisfied with the name and description,
press the Create button to finish the process of creating the new grid.
You cannot press the Create button if no name is specified. If you
change your mind and do not want to create the grid, press the Cancel
button or the X button in the upper right to return to the grid data
manager.
Figure 60. Grid Data Manager accessed from the Components menu.
Figure 61. Creating a new precipitation gridset after pressing the New
button in the Grid Data Manager.
Copying a Grid
There are two ways to copy a grid. Both methods for copying a grid
create an exact duplicate with a different name. Once the copy has
been made it is independent of the original and they do not interact.
The first way to create a copy is to use the grid data manager, which is
accessed from the Components menu. First, select the grid data type
110
Chapter 5 Shared Component Data
you want to copy from the Data Type menu. Then, select the grid you
want to copy by clicking on it in the list of current grids. The selected
grid is highlighted after you select it. After you select a grid you can
press the Copy button on the right side of the window. A new
window will open where you can name and describe the copy that will
be created (Figure 62). A default name is provided for the copy; you
can use the default or replace it with your own choice. A description
can also be entered; if it is long you can use the button to the right of
the description field to open an editor. When you are satisfied with the
name and description, press the Copy button to finish the process of
copying the selected grid. You cannot press the Copy button if no
name is specified. If you change your mind and do not want to copy
the selected grid, press the Cancel button or the X button in the upper
right to return to the grid data manager.
The second way to create a copy is from the Watershed Explorer, on
the Components tab. Move the mouse over the grid data you wish to
copy, then press the right mouse button (Figure 63). A context menu
is displayed that contains several choices including copy. Click the
Create Copymenu option. A new window will open where you can
name and describe the copy that will be created. A default name is
provided for the copy; you can use the default or replace it with your
own choice. A description can also be entered; if it is long you can use
the button to the right of the description field to open an editor. When
you are satisfied with the name and description, press the Copy button
to finish the process of copying the selected grid. You cannot press
the Copy button if no name is specified. If you change your mind and
do not want to copy the grid, press the Cancel button or the X button in
the upper right of the window to return to the Watershed Explorer.
111
Chapter 5 Shared Component Data
Renaming a Grid
There are two ways to rename a grid. Both methods for renaming a
grid changes its name and then all references to the old name are
automatically updated to the new name.
The first way to perform a rename is to use the grid data manager,
which you can access from the Components menu. First, select the
grid data type you want to rename from the Data Type menu. Then,
select the grid you want to rename by clicking on it in the list of current
grids. The selected grid is highlighted after you select it. After you
select a grid you can press the Rename button on the right side of
the window. A new window will open where you can provide the new
name (Figure 64). You can also change the description at the same
time. If the new description will be long, you can use the button to the
right of the description field to open an editor. When you are satisfied
with the name and description, press the Rename button to finish the
process of renaming the selected grid. You cannot press the Rename
button if no name is specified. If you change your mind and do not
want to rename the selected grid, press the Cancel button or the X
button in the upper right of the window to return to the grid data
manager.
The second way to rename is from the Watershed Explorer, on the
Components tab. Select the grid you want to rename by clicking on it
in the Watershed Explorer; it will become highlighted. Keep the mouse
over the selected grid and click the left mouse button again (Figure
65). The highlighted name will change to editing mode. You can then
move the cursor with the arrow keys on the keyboard or by clicking
with the mouse. You can also use the mouse to select some or all of
112
Chapter 5 Shared Component Data
the name. Change the name by typing with the keyboard. When you
have finished changing the name, press the Enter key to finalize your
choice. You can also finalize your choice by clicking elsewhere in the
Watershed Explorer. If you change your mind while in editing mode
and do not want to rename the selected grid, press the Escape key.
Deleting a Grid
There are two ways to delete a grid. Both methods for deleting a grid
will remove it from the project and then automatically update all
references to that grid. Once a grid has been deleted it cannot be
retrieved or undeleted. Any references to the deleted grid will switch to
using no grid, which is usually not a valid choice during a simulation.
At a later time you will have to go to those components and manually
select a different grid.
The first way to perform a deletion is to use the grid data manager,
which you can access from the Components menu. First, select the
grid data type you want to delete from the Data Type menu. Then,
113
Chapter 5 Shared Component Data
select the grid you want to delete by clicking on it in the list of current
grids. The selected grid is highlighted after you select it. After you
select a grid you can press the Delete button on the right side of the
Figure 66. Preparing to delete a precipitation gridset from the Grid Data
Manager.
window (Figure 66). A window will open where you must confirm that
you want to delete the selected grid. Press the OK button to delete the
grid. If you change your mind and do not want to delete the selected
grid, press the Cancel button or the X button in the upper right to
return to the grid data manager.
The second way to delete a grid is from the Watershed Explorer, on
the Components tab. Select the grid you want to delete by clicking
on it in the Watershed Explorer; it will become highlighted. Keep the
mouse over the selected grid and click the right mouse button (Figure
67). A context menu is displayed that contains several choices
including delete. Click the Delete menu option. A window will open
where you must confirm that you want to delete the selected grid.
Press the OK button to delete the grid. If you change your mind and
do not want to delete the selected grid, press the Cancel button or the
X button in the upper right to return to the Watershed Explorer.
114
Chapter 5 Shared Component Data
115
Chapter 5 Shared Component Data
Figure 68. Component editor for a precipitation gridset. Grid data can only
be stored in a HEC-DSS file.
You must also specify the pathname to retrieve from the selected DSS
file (Figure 68). You may type the complete pathname if you know it.
Each pathname contains six parts called the A-part, B-part, C-part, D-
part, E-part, and F-part. The pathname parts are separated with a
slash and may contain spaces. The complete pathname, including
slashes and spaces, can total up to 256 uppercase characters. The
following is an example of a pathname for a percolation parameter
grid:
//SOLDIER CR WATERSHED/PERCOLATION///GIS ESTIMATE/
If you do not know the full pathname of the record you wish you use,
you can use the pathname browser to specify it. You must select a
DSS file first before the browser is available. Press the Select DSS
Pathname button to the right of the "DSS Pathname" field to open the
browser. The browser initially shows all of the records in the specified
DSS file, organized by pathname in the selection table. You can scroll
through the list and select a record pathname by clicking on it. Press
the Select button at the bottom of the browser to choose that record
and return to the Component Editor. If you change your mind and do
not want to select a record pathname, press the Cancel button or the
X button in the upper right of the Select Pathname From HEC-DSS File
window. You can reduce the number of record pathnames shown in
116
Chapter 5 Shared Component Data
the selection table using the "Search by Parts" filters. A separate filter
selection is shown for each of the six pathname parts. By selecting a
choice for a filter, only pathnames that match that choice will be shown
in the selection table. If you make choices in several filters, only
pathnames that satisfy all of the choices will be shown in the selection
table.
Table 19. Internal DSS data type label for different types of grid data.
117
Chapter 5 Shared Component Data
118
Chapter 5 Shared Component Data
119
Chapter 5 Shared Component Data
Table 20. Internal DSS units label for different types of grid data.
120
Chapter 5 Shared Component Data
121
Chapter 5 Shared Component Data
References
Hydrologic Engineering Center. January 2003. HEC Data Storage
System Visual Utility Engine: User's Manual. U.S. Army Corps of
Engineers, Davis, CA.
122
Chapter 6 Watershed Physical Description
CHAPTER 6
Basin Models
Basin models are one of the main components in a project. Their
principle purpose is to convert atmospheric conditions into streamflow
at specific locations in the watershed. Hydrologic elements are used to
break the watershed into manageable pieces. They are connected
together in a dendritic network to form a representation of the stream
system. Background maps can be used to aid in placing the elements
in a spatial context.
123
Chapter 6 Watershed Physical Description
Figure 69. Creating a new basin model. This was accessed by opening the
Basin Model Manager from the Components menu, and then
pressing the New button.
can press the button to the right of the description field to open an
editor. The editor makes it easier to enter and edit long descriptions.
When you are satisfied with the name and description, press the
Create button to finish the process of creating the new basin model.
You cannot press the Create button if no name is specified for the new
model. If you change your mind and do not want to create a new basin
model, press the Cancel button or the X button in the upper right to
return to the Basin Model Manager window.
124
Chapter 6 Watershed Physical Description
Figure 71. Copying a basin model from the Wateshed Explorer. The Copy
Basin Model window will appear after the Create Copy menu
command is selected.
125
Chapter 6 Watershed Physical Description
Figure 72. Renaming a basin model. This was accessed from the Basin
Model Manager.
126
Chapter 6 Watershed Physical Description
127
Chapter 6 Watershed Physical Description
Figure 74. Preparing to delete a basin model from the Basin Model
Manager. A confirmation will be required after pressing the
Delete button.
128
Chapter 6 Watershed Physical Description
mind, you can press the Cancel button or the X button in the upper
right to return to the main program window without importing a model.
After you make your selection, the basin model will be checked for
dependent data such as time-series gages. The user is given the
opportunity to also import any required data. The program
automatically copies the selected file and any additional data into the
project folder and adds the various components to the project.
Figure 76. Basin model component editor showing option properties for
gridded subbasins and flow computations.
Gridded Subbasins
Subbasins that use the ModClark gridded transform method are
considered gridded subbasins. Loss rate and surface transform
calculations will be carried out on a grid cell basis. Properties of the
grid cells are specified in a special grid cell file. The file specifies
which grid cells are in each subbasin, along with the properties of each
cell including location, area within the subbasin, and distance to the
subbasin outlet. Regardless of which transform method is selected, if
gridded meteorological methods will be used, the grid cell file is also
required. One file is used for all subbasins in the basin model and its
format is described in Appendix B. If you wish to use gridded
subbasins, you must prepare the file external to the program and
specify it as part of the basin model parameter data. Specify the
129
Chapter 6 Watershed Physical Description
complete filename of the grid cell file in the Component Editor for the
basin model (Figure 76). Access the Component Editor from the
Watershed Explorer on the "Components" tab by clicking on a basin
model icon. You can use the file browser button to the right of the
entry field to aid in finding the file.
Local Flow
Local flow can be computed at junctions. Local flow is defined as the
sum of all subbasin and source outflows entering a junction. The
subbasin and source elements must be connected directly to the
junction. Any subbasin or source outflow that moves through a routing
element is no longer considered local flow. When local flow is
disabled, a junction element computes outflow as the sum of all inflow
from any type of element. No local flow is computed. When local flow
is enabled, a junction continues to compute outflow as the sum of all
inflow. However, in addition to computing outflow, it additionally
computes local flow (sum of subbasin and source elements). Enabling
or disabling local flow is done in the Component Editor for the basin
model (Figure 76). Access the Component Editor from the Watershed
Explorer on the "Components" tab by clicking on a basin model icon.
Flow Ratio
Flow ratios can be used to increase or decrease the computed flow by
a fixed ratio; they can only be applied to subbasin and source
elements. Once the flow ratios are turned on, each subbasin and
source can have a separate ratio, or no ratio. It is not necessary to
enter a ratio for every element in order to have a ratio at one element.
The ratio 1.0 is used if no ratio is specified. The calculations for
computing outflow proceed normally according to the method choices
and parameter data for each element. At the conclusion of normal
processing, the flow ratio is applied to produce the final outflow.
Flow ratios are enabled or disabled in the Component Editor for the
basin model (Figure 76). After disabling the flow ratios, any ratios
specified for subbasin and source elements in the basin model are
removed and cannot be retrieved again or undeleted. The actual flow
ratio is specified in the Component Editor for the element on the
"Options" tab. Access the editor from the Watershed Explorer on the
"Components" tab by clicking on the element icon. The last tab in the
Component Editor is the "Options" tab. The flow ratio field is disabled
when flow ratios are disabled in the basin model. The field becomes
enabled when flow ratios are enabled in the basin model. The flow
ratio field is never enabled for elements other than subbasins and
sources.
130
Chapter 6 Watershed Physical Description
Missing Flow
Missing inflow data for an element can be set to zero. Under some
conditions it may be possible for source or subbasin elements to
produce outflow with missing values. Downstream routing elements
generally cannot process missing data. When missing flow data is not
replaced, any element that encounters missing inflow data will halt a
simulation with an error message. When missing flow data is replaced,
the missing inflow data is set to zero and a message is generated that
indicates how many values were missing. Processing in the routing
element proceeds normally after any missing inflow data is set to zero.
Setting the action to take with missing inflow data is done on the
Component Editor for the basin model (Figure 76).
Unit System
Each basin model must be in either United States customary units
(sometimes called English units) or in system international units (also
called metric units). All parameter data in a basin model must be in the
same unit system. If you change the unit system, all data will be
automatically converted to the new unit system. All time-series data,
paired data, and gridded data referenced in a basin model will be in its
own unit system. If necessary these referenced data are automatically
converted to the unit system of the basin model during a simulation.
Select the unit system using the Component Editor for the basin model
(Figure 76). Access the Component Editor by clicking the basin model
icon on the "Components" tab of the Watershed Explorer. If you
change the unit system, all data is automatically converted to the new
selection.
Sediment
The movement of sediment in the watershed can be included as part of
the hydrology of the basin model. When the sediment is disabled, no
sediment processing will take place in the basin model. When the
sediment is enabled, appropriate sediment simulation components are
added to elements in the basin model. Subbasin elements include
surface erosion and wash-off. Reach elements include erosion,
deposition, and sediment transport. Reservoir elements include
sediment settling and sediment transport. Source, junction, diversion,
and sink elements pass sediment through the channel network.
Water Quality
Nutrient water quality (nitrogen and phosphorus) can be included as
part of the hydrology of the basin model. When the water quality is
disabled, no water quality processing will take place in the basin
131
Chapter 6 Watershed Physical Description
Background Maps
Background maps provide a spatial context for the hydrologic elements
composing a basin model. The maps are not actually used in the
compute process, but they can be very helpful in showing the spatial
relationship between elements. They are commonly used for showing
the boundaries of a watershed or the location of streams. They can
also be useful for showing supplemental information such as the
location of levees. The use of background maps is always optional.
Background maps are added and removed using the Background
Maps window (Figure 77). Access the window by clicking the View
menu and selecting the Background Maps command. The menu
command is only available if the Basin Model Map window is currently
open. Any changes made to the current background maps or their
draw order take affect immediately as soon as they are made.
Add a map to the list of current background maps by pressing the
Add button. A file browser opens for you to select the map file you
wish to add. Six different map file formats can be used. The HMS
map file and geo-reference image file formats are described in
Appendix D. File formats for the other map types are defined by their
creators (Autodesk 2005, ESRI 1997, USGS 1999). Select the type of
map file you wish to add using the "Files of type" selection list at the
bottom of the file browser.
Locate the map file you wish to add and select it. Press the Select
button to add the map to the list of current background maps. If you
change your mind and do not want to add a map, press the Cancel
button or the X button in the upper right of the Select window. The
program does not include any coordinate transformation tools so you
132
Chapter 6 Watershed Physical Description
must independently make certain that all maps you wish to use are in
the same coordinate system.
Figure 77. Background maps can be added to provide a spatial context for
the basin model elements. Maps can be moved up or down in
the drawing order.
133
Chapter 6 Watershed Physical Description
Maximum Extents
The maximum extents are designed to limit the area shown in the
Basin Model Map window to a reasonable view. This is especially
helpful if background maps are used which cover a very large region;
much larger than the watershed. Four different methods are available
for specifying the maximum extents. Setting the maximum extents is
done from the Maximum Extents window (Figure 78) accessed by
selecting the View Maximum Extents command.
The Manual option will set the extents to a region specified by the
user. When this option is selected, the "Manual view extents" box will
be available. The "Minimum Northing" corresponds to the South edge
and the "Maximum Northing" corresponds to the North edge of the
extents. The "Minimum Easting" corresponds to the West edge and
the "Maximum Easting" corresponds to the East edge of the extents.
You may enter the appropriate values in the coordinate system of the
current background maps. The northing and easting values can
automatically be filled in from the current view in the basin map display
by pressing the Set button.
The Union of All Elements option will set the extents to the
combined extents of all the hydrologic elements in the basin model. As
elements are added and removed, the maximum extents will
automatically be updated. The "Element Buffer" can be used to make
sure there is some space around the elements. When the extents are
computed, the buffer is computed as the percentage of the width and
height of the elements themselves.
134
Chapter 6 Watershed Physical Description
Figure 78. The maximum extents are used to control the view size of the
basin model, and facilitate use of maps larger than the
watershed of interest.
The Union of All Maps option will set the extents to the combined
extents of all background maps. As maps are added and removed, the
maximum extents will automatically be updated.
The Union of All Maps and Elements options will set the
extents to the combined extents of all background maps and
hydrologic elements in the basin model. As elements are added and
removed, or as background maps are added and removed, the extents
will automatically be updated. The "Element Buffer" will be applied to
the elements as part of computing the extents.
Background Gridlines
Background gridlines are useful for helping to measure approximate
distance in the Basin Model Map window. They are not drawn at a
specific scale; the relative scale changes depending on the current
view. There is only one way to turn the grid lines on and off.
Turn the gridlines on and off using the View menu. Click on the View
menu and select the Draw Gridlines command. A check is shown
next to the menu command to indicate that the gridlines are turned on.
The lines are drawn first in the map using a light gray color, then all
other maps and the hydrologic elements are drawn. You can turn off
the gridlines by clicking on the same menu command again.
135
Chapter 6 Watershed Physical Description
136
Chapter 6 Watershed Physical Description
where you wish to begin and hold the left mouse button. Drag the
mouse in any direction to change the view. The new view will be
redrawn when you release the mouse button.
It can be very difficult to find the selected element or elements in a
basin model with many hydrologic elements. A special tool has been
created to zoom to the current element selection. Begin by creating an
element selection. You can click on an element in the Watershed
Explorer, or use the Edit Select Special command to create a
current selection. You can zoom to the selection by clicking the View
menu and selecting the Zoom To Selected command. The Basin
Model Map window will automatically zoom to an appropriate scale for
the current element selection and center the view on the selection.
After zooming in and moving around the basin model, you may want to
quickly return to viewing the entire model. You can immediately zoom
to the maximum extents. Click the View menu and select the Zoom
To Maximum Extents command. The Basin Model Map window will
automatically zoom out to the coordinates specified as the maximum
extents and redraw.
137
Chapter 6 Watershed Physical Description
The geometric data necessary to provide the map objects drawing for
elements is stored in external files that use the ESRI Shape File
format. Three separate files are used for subbasin, reservoir, and
reach geometric data. Currently the files must be created externally to
the program. Additionally, the files must be stored in the same
directory as the basin model, which is usually stored in the project
directory. The files use a strict naming convention as illustrated below:
Basin Model "Baseline" Baseline.basin
Subbasin polygon geometry Baseline_Subbasin.shp
Reservoir polygon geometry Baseline_Reservoir.shp
Reach polygon geometry Baseline_Reach.shp
138
Chapter 6 Watershed Physical Description
Hydrologic Elements
Hydrologic elements are the basic building blocks of a basin model.
An element represents a physical process such as a watershed
catchment, stream reach, or confluence. Each element represents
part of the total response of the watershed to atmospheric forcing.
Seven different element types have been included in the program:
subbasin, reach, reservoir, junction, diversion, source, and sink.
An element uses one or more mathematical models to describe a
physical process occurring in the watershed. Sometimes the model is
only a good approximation of the original physical process over a
limited range of environmental conditions. Data availability and the
required parameters of a model can also determine fitness. To make
the program suitable for many different conditions, most elements have
more than one model or method for approximating the physical
process. For example, there are seven different methods for
specifying overland flow for a subbasin.
139
Chapter 6 Watershed Physical Description
create the element. Click the left mouse button. A window will open
where you can name and describe the new hydrologic element that will
be created, as shown in Figure 79. A default name is provided for the
new element; you can use the default or replace it with your own
choice. A description can also be entered. If the description is long,
you can press the button to the right of the description field to open an
editor. The editor makes it easier to enter and edit long descriptions.
When you are satisfied with the name and description, press the
Create button to finish the process of creating the new hydrologic
element. You cannot press the Create button if no name is specified
for the new element. If you change your mind and do not want to
create a new hydrologic element, press the Cancel button or the X
button in the upper right to return to the Basin Model Map window.
After you finish creating one element, the element creation tool is still
selected and you can create additional new elements of the same type.
Figure 79. Creating a subbasin in a new basin model. The subbasin tool
was selected on the toolbar. The cursor changed from cross
hairs back to a pointer when the left mouse button was pressed
140
Chapter 6 Watershed Physical Description
Most of the elements are created by clicking the mouse button over the
location where you want to create the element. However, two mouse
clicks are required when creating a reach element. To create a reach,
begin by selecting the correct tool from the toolbar. Next click once
over the location you want to be the upstream end of the reach. Then
as you move the mouse you will see a line connecting the mouse to
the upstream end of the reach. Move the mouse to the location you
want to be the downstream end of the reach and click. You will finish
the process of creating a reach by selecting a name and entering the
optional description.
Copying an Element
Hydrologic elements can be copied to the clipboard from one basin
model and then exact duplicates can be pasted in either the same or a
different basin model. Begin by opening the basin model that contains
the element or elements you wish to copy. Select the hydrologic
element you wish to copy by clicking on it with the arrow tool. You may
also select more than one element in the map. The selected element
or elements become highlighted after the selection. After you make a
selection, place the mouse over a selected element and press the right
mouse button as shown in Figure 80. A context menu is displayed that
contains several choices including copy. Click the Copy Element
command. If more than one element is selected, select the Copy
Elements command. You can now paste the selection into the same
or a different basin model. The copy command is only available if
there is at least one element selected.
141
Chapter 6 Watershed Physical Description
Figure 80. Copying a group of selected elements from the basin map to the
clipboard; the elements will remain in the basin model. If the
elements were cut, they would be on the clipboard but removed
from the basin model.
Pasting an Element
An element in the clipboard can be pasted into a basin model. Pasting
an element does not remove it from the clipboard so you can paste the
same selection into multiple basin models. Begin by opening the basin
model into which you wish to paste elements from the clipboard.
Select the arrow tool and move the mouse to the location where you
wish to paste the element. If there are multiple elements on the
clipboard then you should place the mouse where you want the center
of the element grouping to be located. Press the right mouse button as
shown in Figure 81. A context menu is displayed that contains several
choices including paste. Click the Paste Element command. The
paste command is only available if there is at least one element on the
clipboard.
142
Chapter 6 Watershed Physical Description
Figure 81. Pasting an element selection into the basin map. The mouse
was placed where the pasted element should be added to the
basin map and the right mouse button was pressed. Now the
Paste Element command can be selected.
Cutting an Element
An element can be cut from a basin model and removed to the
clipboard. Cutting an element places an exact but independent copy of
the element on the clipboard and then deletes it from the basin model.
Once an element is in the clipboard it can be pasted into the same
basin model from which it was copied or it can be pasted into a
different model. Begin by opening the basin model from which you
wish to cut an element. Select the hydrologic element you wish to cut
by clicking on it with the arrow tool. You may also select more than
one element in the map. The selected element or elements become
highlighted after the selection. After you make a selection, place the
mouse over a selected element and press the right mouse button. A
143
Chapter 6 Watershed Physical Description
Renaming an Element
Rename a hydrologic element using the "Components" tab of the
Watershed Explorer. Select the element you wish to rename by
clicking on it in the Watershed Explorer; it will become highlighted.
Keep the mouse over the selected element and click the left mouse
button again. The highlighted name will change to editing mode as
shown in Figure 82. You can then move the cursor with the arrow keys
on the keyboard or by clicking with the mouse. You can also use the
mouse to select some or all of the name. Change the name by typing
with the keyboard. When you have finished changing the name, press
the Enter key to finalize your choice. You can also finalize your choice
by clicking elsewhere on the "Components" tab. If you change your
mind while in editing mode and do not want to rename the selected
hydrologic element, press the Escape key.
Figure 82. Renaming an element in the Watershed Explorer. You can right-
click on an element and select the Rename command.
Alternately, click with the left mouse button on an element that is
already selected to immediately enter editing mode and change
the name.
Deleting an Element
There are three ways to delete a hydrologic element. All methods for
deleting an element remove it from the basin model. Once an element
has been deleted it cannot be retrieved or undeleted.
The first way to perform a deletion is from the Basin Model Map
window. Select the hydrologic element you wish to delete by clicking
on it with the arrow tool. You may also select more than one element
in the map. The selected element or elements become highlighted
after the selection. After you make a selection, place the mouse over a
selected element and press the right mouse button (Figure 83). A
144
Chapter 6 Watershed Physical Description
145
Chapter 6 Watershed Physical Description
146
Chapter 6 Watershed Physical Description
147
Chapter 6 Watershed Physical Description
Element Inventory
The element inventory provides a listing of some or all of the
hydrologic elements in a basin model. This is helpful for reviewing the
various element names and descriptions. To access the inventory,
click the Parameters menu and select the Element Inventory
command. A sample inventory is shown in Figure 86. If there is a
current element selection, only the selected elements will be shown in
the inventory when the window opens; you can switch to showing all
elements by using the selection control at the top of the window. All
hydrologic elements in the basin model will be shown if no elements
are currently selected.
Figure 86. Typical element inventory for a basin model. You can switch
between viewing the initial element selection from the basin
map, or all elements.
148
Chapter 6 Watershed Physical Description
You can also find and select a hydrologic element in the Basin Model
Map window. Begin by making sure the Arrow Tool is selected on the
toolbar. Examine the element icons in the Basin Model Map window
until you find the one you wish to select. Click on it with the arrow
cursor and it will become selected. The selected element is
highlighted in the Basin Model Map window and in the Watershed
Explorer. After you make the first selection, you can add to the
selection by holding the shift key and clicking on additional elements.
You can select several elements simultaneously by using the arrow
cursor to drag a box around the desired elements. To select using a
drag box, move the mouse to a blank area of the map. Hold the left
mouse button and drag the mouse. A box is drawn to show which
elements will be selected when you release the mouse button.
It can be difficult to find a specific hydrologic element in a large basin
model with many elements. A special tool is included that can be used
to help locate a single element or a group of elements. Access the tool
by clicking the Edit menu and selecting the Select Special command
(Figure 87). The Select Special window will remain open until it is
closed by pressing the Close button or the X button at the upper
corner of the window. Elements are selected in the basin model
according to the selections on the Select Special window every time
the Apply button is pressed. There are three components to a special
selection.
The first component of a special selection is the element name. You
may enter a specific name, for example, the name of an element
appearing in an error message. You can find all elements that begin
with a specific sequence of letters by entering those letters followed by
an asterisk. For example, entering sub* will find all elements with
names that start with "sub". You can also find all elements that end
with a specific sequence of letters by entering an asterisk followed by
those letters. For example, entering *creek will find all elements with
names that end in "creek". Finally you can find all elements that
contain a specific sequence of letters by entering an asterisk, the
letters, and ending with another asterisk. For example, entering
*basin* will find all elements with names that contain "basin". By
default the "Element Name" is set to an asterisk so that the special
selection will find all element names.
The second component of a special selection takes the result of
limiting by element name and further limits it based on the element
type. Elements will only be selected if they match the name criteria
and their type is checked on. For example, to only search among the
source elements, check off all element types except source. You may
check on one, several, or all element types. The default is to search
among all element types.
149
Chapter 6 Watershed Physical Description
Figure 87. Using the special selection tool to select all junctions whose
names begins with the letters "West".
Flow Network
The flow network is the skeleton that connects hydrologic elements
together into a representation of the stream system in the watershed.
Each link in the network is a one-way connector that takes outflow from
an element and connects it as inflow to a downstream element. The
connection information of the flow network along with the drainage
area at each element is used to sort the elements in hydrologic order.
Moving Elements
Hydrologic elements are moved in the Basin Model Map window; they
are moved the same way regardless of whether they are connected to
150
Chapter 6 Watershed Physical Description
Rescaling Elements
The best practice for creating basin models is to add background maps
before creating any elements. However, sometimes a suitable
background map is not immediately available and must be added later.
At other times it is desirable to add elements from a different basin
model into the current basin model. In both situations it is rare that the
elements are at the correct locations or even in the same scale. For
example, the elements already in a basin model may have horizontal
151
Chapter 6 Watershed Physical Description
units that range from 100 to 125. A background map that must be
added may have horizontal units that range from 10,000 to 15,000. A
similar situation applies when combining elements from multiple basin
models. In these cases it is necessary to rescale the elements so they
cover the same range as a background map or a different set of
elements. This rescaling should be done before attempting to merge
all the elements into a single flow network.
The Rescale Elements command on the View menu can be used to
adapt elements to a background map or elements from another basin
model. Begin by selecting two or more elements in the Basin Model
Map window. Choose the appropriate command from the Rescale
Elements menu command. Selecting the 25% or 50% option will
reduce the area covered by the elements, while selecting the 150% or
200% will increase the area covered by the elements. The elements
will anchor to the bottom, left corner. When the scale of the elements
is reduced, the elements will shift from the top toward the bottom, and
from the right edge toward the left. Conversely, when the scale of the
elements is increased, the elements will shift up and to the right. The
scaling is applied linearly so that elements in the middle of the group
move half as much as elements at the top or right edge of the group.
152
Chapter 6 Watershed Physical Description
connected, the command will not appear on the menu. After you click
the Connect Downstream command, the cursor will change to cross
hairs as shown in Figure 88. Move the mouse to position the cross
hairs over the element to which you wish to connect. Once the mouse
is positioned, press the left mouse button. A connection link will be
shown between the upstream and downstream elements.
To disconnect an element in the Basin Model Map window, start by
selecting the Arrow Tool from the toolbar. Next click the icon of the
element you wish to disconnect from its downstream element; the
element you click will become highlighted. Click the right mouse
button to see a context menu that includes Delete Connection. If the
element is not connected, the command will not appear. After you
click the Delete Connection command, the connection link between
the upstream and downstream element will be removed.
153
Chapter 6 Watershed Physical Description
You can also create and delete connections from the Component
Editor for the element (Figure 89). To access the Component Editor,
click on the desired element in the Basin Model Map window, or on the
"Components" tab of the Watershed Explorer. The Component Editor
shows the downstream element. If the element is not connected to a
downstream element, the selection list will show None as the selected
element. To connect the element currently shown in the Component
Editor, select an element on the list. Only elements that could
potentially be downstream of the current element are shown. To
disconnect the current element, change the selection in the list to the
None selection.
Hydrologic Order
Hydrologic order means that elements from the basin model are shown
in order from the headwaters to the outlet. The key properties used to
determine hydrologic order are the flow network connections and
drainage area. From each element in the basin model it is possible for
the program to automatically determine the downstream element, and
the upstreams elements. These connections create the basic structure
of the network. The area entered for subbasins and sources is used to
determine the drainage area for any other downstream element in the
network. While area is optional in source elements for many purposes,
the hydrologic order downstream of a source cannot be properly
determined if the area is not specified. A simple basin model is shown
in Figure 90, with the corresponding Watershed Explorer shown in
Figure 91.
154
Chapter 6 Watershed Physical Description
155
Chapter 6 Watershed Physical Description
156
Chapter 6 Watershed Physical Description
Figure 92. Manually moving an element in the hydrologic order using the
right mouse menu.
157
Chapter 6 Watershed Physical Description
158
Chapter 6 Watershed Physical Description
Zones
Zones provide a way to group together elements with similar
properties. The grouping could be performed spatially so that a large
watershed is subdivided with a separate zone for each tributary.
Another way to subdivide spatially is to group all subbasins between
upstream and downstream gage locations. Alternatively, the grouping
could be performed parametrically so that a watershed is subdivided
with a separate zone for each type of land use. The properties used to
group elements into zones are defined by the user and may use any
meaningful criteria. Zones are defined and managed as part of the
basin model. Once defined, the primary use of zones is for setting
parameter values within the Forecast Alternative, described fully in a
later chapter. In brief, during a forecast it becomes possible to set
parameter values once and have them immediately applied to all
subbasins in a zone.
There are three organizational levels to the use of zones. The lowest
organizational level is individual subbasins. The subbasins contain the
actual parameter data and are located within the watershed. The
middle organizational level is a Zone. Each Zone contains one or
more subbasins. Subbasins are grouped into zones using one of the
schemes described in the previous paragraph. The highest level is a
Zone Configuration. Each Zone Configuration includes one or more
Zones. A Basin Model may include multiple configurations.
159
Chapter 6 Watershed Physical Description
160
Chapter 6 Watershed Physical Description
161
Chapter 6 Watershed Physical Description
Figure 97. Renaming a zone configuration. This was accessed from the
Zone Configuration Manager.
162
Chapter 6 Watershed Physical Description
163
Chapter 6 Watershed Physical Description
164
Chapter 6 Watershed Physical Description
165
Chapter 6 Watershed Physical Description
Figure 103. Elements are added to the selected zone "Alligator Creek" by
pressing the Add button. Multiple elements were selected by
holding the control key and clicking each desired element.
Elements can also be added by double-clicking each one in the
left panel.
There are two ways to add an element to the selected zone. The first
way is to click on an unassigned element in the left side of the
Elements window. You may hold the control key and click on
additional elements to add more than one element at a time. The
selected elements are highlighted after clicking. Press the Add button
to add the selected elements to the selected zone.
The second way to add an element to the selected zone is to double-
click it in the list of unassigned elements. When an element is double-
166
Chapter 6 Watershed Physical Description
Figure 104. Elements are returned from the selected zone "Alligator Creek"
to the unassigned elements by pressing the Remove button.
Elements can also be removed by double-clicking each one
shown in the right panel.
The right side of the Elements window shows all elements that have
been assigned to a zone. The zones selector at the upper right of the
window can be used to switch between zones in the selected zone
configuration. You must select the zone from which you wish to
167
Chapter 6 Watershed Physical Description
remove elements. The left side of the Elements window shows all
elements have not been assigned to a zone.
There are two ways to remove an element from the selected zone.
The first way is to click on an element in the right side of the Elements
window. You may hold the control key and click on additional
elements to remove more than one element at a time. The selected
elements are highlighted after clicking. Press the Remove button to
return the selected elements to the list of unassigned elements.
The second way to remove an element from the selected zone is to
double-click it in the list of zone elements. When an element is double-
clicked, it will be removed from the currently selected zone and placed
on the list of unassigned elements.
168
Chapter 6 Watershed Physical Description
Figure 105. Assigning the "AligCr-04" subbasin to the "Alligator Creek" zone
in the current zone configuration.
169
Chapter 6 Watershed Physical Description
Figure 106. Viewing the parameter data for subbasins assigned to zone
"Alligator Creek" in the global editor, and preparing to view the
subbasins assigned to zone "Hackberry Creek."
170
Chapter 6 Watershed Physical Description
Figure 107. Viewing the results for subbasins assigned to zone "Alligator
Creek" in the global summary table, and preparing to view
subbasins assigned to zone "Hackberry Creek."
Computation Points
Elements in the basin model can be designated as computation points.
Generally only elements with observed flow, observed stage, or other
observed data are designated as computation points. Computation
points are designed primarily for use with simulation runs (Chapter 12)
and forecast alternatives (Chapter 14). When used with a simulation
run, the features of a computation point are a customizable editor that
uses slider bars to adjust parameters upstream of the computation
point, and result graphs with customizable time-series selection.
Results from the simulation run are automatically recomputed as the
slider bars are adjusted and the result graphs update dynamically.
When used with forecast alternatives, the principal feature of a
computation point is to make upstream to downstream calibration more
efficient.
171
Chapter 6 Watershed Physical Description
using the right mouse button (Figure 108). Choose the Select
Computation Point command to select the element as a computation
point. It is only possible to select the element as a computation point if
it is not already selected. An element selected as a computation point
shows a small red circle added to the icon in the Basin Model Map
window and also in the Watershed Explorer.
172
Chapter 6 Watershed Physical Description
selected; the selected row in the table is highlighted. Press the Close
button when you are done selecting elements and you will return to the
Computation Point Manager window.
173
Chapter 6 Watershed Physical Description
174
Chapter 6 Watershed Physical Description
175
Chapter 6 Watershed Physical Description
Figure 113. Preparing to select parameters for the customizable editor at the
"Hackberry Creek" computation point.
Figure 114. The parameters management window shows all the parameters
for a computation point selected in the Computation Point
Manager.
176
Chapter 6 Watershed Physical Description
177
Chapter 6 Watershed Physical Description
178
Chapter 6 Watershed Physical Description
179
Chapter 6 Watershed Physical Description
Figure 117. Changing the minimum and maximum parameter values in the
customizable editor at a computation point.
180
Chapter 6 Watershed Physical Description
Figure 118. Selected results for the first graph at computation point
"Hackberry Creek." The outflow and observed flow are included
from the computation point, plus the outflow from several
upstream subbasins.
181
Chapter 6 Watershed Physical Description
Figure 119. Selecting the outflow time-series from a subbasin element for
addition to a result graph at a computation point. Multiple time-
series can be selected simultaneously by holding the control key
and selecting additional time-series results.
182
Chapter 6 Watershed Physical Description
References
Autodesk, Inc. 2005. "AutoCAD 2006 Customization Guidebook."
San Rafael, CA.
Barnes Jr., Harry H. 1967. "Roughness Characteristics of Natural
Channels." U.S. Geological Survey, Denver, CO.
Environmental Systems Research Institute, Inc. July 1998. "ESRI
Shapefile Technical Description." Redlands, CA.
U.S. Geological Survey. September 1999. "Standards for Digital Line
Graphs: Part 2 Specifications." Department of the Interior,
Washington, DC.
183
Chapter 6 Watershed Physical Description
184
Chapter 7 Subbasin Elements
CHAPTER 7
Subbasin Elements
A subbasin is an element that usually has no inflow and only one
outflow. It is one of only two ways to produce flow in the basin model.
Outflow is computed from meteorologic data by subtracting losses,
transforming excess precipitation, and adding baseflow. Optionally the
user can choose to include a canopy component to represent
interception and evapotranspiration. It is also optional to include a
surface component to represent water caught in surface depressional
storage. The subbasin can be used to model a wide range of
catchment sizes.
Figure 121. Subbasin component editor. All element editors include the
basin model and element names, description, and downstream
connection.
185
Chapter 7 Subbasin Elements
Dynamic Canopy
This method includes an interception storage capacity and a crop
coefficient that change in time. The storage capacity specifies the
amount of water that can be held in the canopy before precipitation
begins falling through to the ground surface. No through-fall will occur
until the canopy storage is full. The canopy will use all potential evapo-
transpiration until the storage is emptied. The potential
evapotranspiration is multiplied by the crop coefficient to determine the
amount of evapotranspiration from canopy storage and later the
surface and soil components. Only after the canopy storage has been
emptied will unused potential evapotranspiration be used by the
186
Chapter 7 Subbasin Elements
187
Chapter 7 Subbasin Elements
188
Chapter 7 Subbasin Elements
The canopy can be set to only evaporate water from storage and
extract water from the soil during dry time periods with no precipitation.
Version 4.1 of the program and all previous versions only performed
evapotranspiration calculations during dry periods. Alternately the
canopy can be set to evaporate water from storage and extract water
from the soil during both dry and wet periods. The choice for
simultaneous precipitation and evapotranspiration can improve results
when using a long time interval or during a snowmelt simulation.
There are two methods for extraction of water from the soil. The
Simple method extracts water at the potential evapo-transpiration
rate, and can be used with the deficit constant or soil moisture
accounting loss rate methods. The Tension Reduction method
extracts water at the potential evapo-transpiration rate from the gravity
zone but reduces the rate when extracting from the tension zone; it can
only be used with the soil moisture accouting loss rate method.
Simple Canopy
This method is a simple representation of a plant canopy. All
precipitation is intercepted until the canopy storage capacity is filled.
Once the storage is filled, all further precipitation falls to the surface, or
directly to the soil if no representation of the surface is included. All
potential evapotranspiration will be used to empty the canopy storage
until the water in storage has been eliminated. The potential
evapotranspiration is multiplied by the crop coefficient to determine the
amount of evapotranspiration from canopy storage and later the
surface and soil components. Only after the canopy storage has been
emptied will unused potential evapotranspiration be used by the
surface and soil components. The Component Editor is shown in
Figure 124.
The initial condition of the canopy should be specified as the
percentage of the canopy storage that is full of water at the beginning
of the simulation.
189
Chapter 7 Subbasin Elements
190
Chapter 7 Subbasin Elements
191
Chapter 7 Subbasin Elements
Simple Surface
This method is a simple representation of the soil surface. All
precipitation or precipitation through-fall that arrives on the soil surface
is captured in storage until the storage capacity of the surface is filled.
Water in surface storage infiltrates into the soil whenever it is present
in storage. That is, water will infiltrate even when the storage capacity
is not full. Surface runoff will be generated when the storage capacity
is filled, and the precipitation through-fall rate exceeds the infiltration
rate. The Component Editor is shown in Figure 126.
The initial condition of the surface should be specified as the
percentage of the surface storage that is full of water at the beginning
of the simulation.
Surface storage represents the maximum amount of water that can be
held on the soil surface before surface runoff begins. The amount of
storage is specified as an effective depth of water.
192
Chapter 7 Subbasin Elements
193
Chapter 7 Subbasin Elements
The initial deficit is the initial condition for the method. It is the amount
of water required to fill the soil layer to the maximum storage.
The maximum deficit specifies the amount of water the soil layer can
hold, specified as a depth. An upper bound is the depth of the active
soil layer multiplied by the porosity. However, in most cases such an
estimate must be reduced by the permanent wilting point and for other
conditions that reduce the water holding capacity.
The constant rate defines the percolation rate when the soil layer is
saturated. Saturated hydraulic conductivity is a good approximation.
The percentage of the subbasin which is directly connected impervious
area can be specified. No loss calculations are carried out on the
impervious area; all precipitation on that portion of the subbasin
becomes excess precipitation and subject to surface storage and direct
runoff.
194
Chapter 7 Subbasin Elements
Exponential Loss
The exponential loss method is empirical and generally speaking
should not be used without calibration. It represents incremental
infiltration as a logarithmically decreasing function of accumulated
infiltration. It includes the option for increased initial infiltration when
the soil is particularly dry before the arrival of a storm. It should not be
used for continuous simulation because it is a function of cumulative
infiltration and does not include any type of recovery. The Component
Editor is shown in Figure 128.
The initial range is the amount of initial accumulated infiltration during
which the loss rate is increased. This parameter is considered to be a
function primarily of antecedent soil moisture deficiency and is usually
storm-dependent.
The initial coefficient specifies the starting loss rate coefficient on the
exponential infiltration curve. It is assumed to be a function of
infiltration characteristics and consequently may be correlated with soil
type, land use, vegetation cover, and other properties of a subbasin.
195
Chapter 7 Subbasin Elements
196
Chapter 7 Subbasin Elements
A maximum deficit grid must be selected from the list of choices. The
selection list will show all moisture deficit grids available in the grid
data manager. You can use a chooser to select a grid by pressing the
grid button next to the selection list. The chooser shows all of the
moisture deficit grids in the grid data manager. Click on a grid to view
the description. This is especially helpful since initial deficit and
maximum storage grids are both stored as moisture deficit grids.
A constant rate grid must be selected from the list of choices. The
selection list will show all percolation rate grids available in the grid
data manager. You can use the chooser to select a grid by pressing
the grid button next to the selection list, or select directly from the list of
choices.
An impervious grid must be selected from the list of choices. The
selection list will show all impervious area grids available in the grid
197
Chapter 7 Subbasin Elements
data manager. You can use the chooser to select a grid by pressing
the grid button next to the selection list.
The initial deficit grid ratio must be entered. The actual initial deficit
value at each grid cell is obtained by multiplying the value obtained
from the initial deficit grid by the specified ratio. The default value is
1.0.
The maximum deficit grid ratio must be entered. The actual maximum
deficit value at each grid cell is obtained by multiplying the value
obtained from the maximum deficit grid by the specified ratio. The
default value is 1.0.
The constant rate grid ratio must be entered. The actual constant loss
rate value at each grid cell is obtained by multiplying the value
obtained from the constant rate grid by the specified ratio. The default
value is 1.0.
The impervious grid ratio must be entered. The actual percentage of
impervious area at each grid cell is obtained by multiplying the value
obtained from the impervious grid by the specified ratio. The default
value is 1.0.
198
Chapter 7 Subbasin Elements
Figure 131. Gridded Green and Ampt loss rate component editor.
The wetting front suction grid must be selected from the list of choices.
The selection list will show all of the water potential grids available in
the grid data manager. You can use the chooser to select a grid by
pressing the grid button next to the selection list.
The hydraulic conductivity grid must be selected from the list of
choices. The selection list will show all of the percolation rate grids
specified in the grid data manager, because conductivity is a type of
percolation and has the same units. You can use the chooser to select
a grid by pressing the grid button next to the selection list.
An impervious grid must be selected from the list of choices. The
selection list will show all impervious area grids available in the grid
data manager. You can use the chooser to select a grid by pressing
the grid button next to the selection list.
199
Chapter 7 Subbasin Elements
200
Chapter 7 Subbasin Elements
The maximum infiltration grid is selected from the percolation grids that
have been previously defined in the grid data manager. The grid
should specify the maximum infiltration rate at each grid cell. This is
the upper bound on infiltration; the actual infiltration at any cell in a
particular time interval is a linear function of the surface and soil
storage in the cell, if a surface method is selected. Without a selected
surface method, water will always infiltrate at the maximum rate. You
may use a chooser to select the grid. The grid selections will be
disabled unless you have previously created grids in the grid data
manager. Likewise, you must also select a soil percolation and
groundwater percolation grids. All infiltration and percolation grids use
the same type of parameter grid so descriptions for the grids are
important.
The initial condition of the soil should be specified as the percentage of
the soil storage that is full of water at the beginning of the simulation.
The same percentage will be applied to every grid cell. Likewise, you
must specify the initial storage of each groundwater layer.
The soil storage grid must be selected from the grids that have been
previously defined in the grid data manager. The grid should specify
the maximum soil storage in each grid cell. You may use a chooser to
select the grid by pressing the grid button next to the selection list.
You will not be able to select a grid if no grids have been created in the
grid data manager. Likewise, you must select a tension storage grid,
and a storage grid for each groundwater layer. Because tension
storage is contained within the total soil storage, the tension storage at
each cell must be less than the soil storage at the same cell.
201
Chapter 7 Subbasin Elements
Groundwater coefficient grids must be selected for the upper and lower
groundwater layers. The selected grid should specify the storage
coefficient for each cell in the layer. The coefficient is used as the time
lag on a linear reservoir for transforming water in storage to become
lateral outflow. Contributions from each grid cell are accumulated to
determine the total amount of flow available to become baseflow.
202
Chapter 7 Subbasin Elements
the difference in volume at the end of two adjacent time intervals. The
Component Editor is shown in Figure 135.
You may optionally enter an initial abstraction. The initial abstraction
defines the amount of precipitation that must fall before surface excess
results. However, it is not the same as an initial interception or initial
loss since changing the initial abstraction changes the infiltration
response later in the storm. If this value is left blank, it will be
automatically calculated as 0.2 times the potential retention, which is
calculated from the curve number.
You must enter a curve number. This should be a composite curve
number that represents all of the different soil group and land use
combinations in the subbasin. The composite curve number should
not include any impervious area that will be specified separately as the
percentage of impervious area.
The percentage of the subbasin which is directly connected impervious
area can be specified. Any percentage specified should not be
included in computing the composite curve number. No loss
calculations are carried out on the impervious area; all precipitation on
that portion of the subbasin becomes excess precipitation and subject
to direct runoff.
203
Chapter 7 Subbasin Elements
The initial water content gives the initial saturation of the soil at the
beginning of a simulation. It should be specified in terms of volume
ratio.
The residual water content specifies the amount of water remaining in
the soil after all drainage has ceased. It should be specified in terms of
volume ratio. It may be determined in the laboratory or estimated from
the soil texture.
The saturated water content specifies the maximum water holding
capacity in terms of volume ratio. It is often assumed to be the total
porosity of the soil.
The bubbling pressure, also known as the wetting front suction, must
be specified. It is generally assumed to be a function of the soil
texture.
The pore size distribution determines how the total pore space is
distributed in different size classes. It is typically assumed to be a
function of soil texture.
The hydraulic conductivity must also be specified, typically as the
effective saturated conductivity. It can be estimated from field tests or
approximated by knowing the soil texture.
The percentage of the subbasin which is directly connected impervious
area can be specified. No loss calculations are carried out on the
impervious area; all precipitation on that portion of the subbasin
becomes excess precipitation and subject to direct runoff.
Optionally, a temperature gage may be selected for adjusting the water
density, water viscosity, and matric potential based on temperature. If
204
Chapter 7 Subbasin Elements
205
Chapter 7 Subbasin Elements
206
Chapter 7 Subbasin Elements
207
Chapter 7 Subbasin Elements
208
Chapter 7 Subbasin Elements
Figure 139. Plane editor for the kinematic wave transform method. The
same properties are required for the second plane, if it is used.
The first parameter is the length. The typical length of the conceptual
plane should be entered. For impervious areas, this should be the
average flow length from the point where precipitation falls, to where
the runoff first enters a collection gutter or channel. For pervious
areas, this should likewise be the average flow length. When using the
two planes to represent pervious and impervious areas, it is helpful to
adopt a convention for which plane number (one or two) represents the
pervious area.
The slope is also required for each plane that will be used. The slope
should represent the average slope along the flow line from the point
where precipitation falls to where the runoff first enters a gutter or
channel.
Roughness is the principal difference between a plane meant to
represent impervious or pervious area. Roughness coefficients for
natural areas are much higher than for developed areas. However, in
both cases the surface roughness coefficients are higher than typical
Manning's roughness coefficients used for open channel flow.
209
Chapter 7 Subbasin Elements
Figure 140. Subcollector editor for the kinematic wave transform method.
The same properties are required for the collector, if it is used.
The length should be the average distance from the beginning of the
subcollector to where it will enter a collector. Usually this value will be
measured from maps of the stormwater collection system.
The slope should be the average slope along the average flow length.
This value may be estimated from maps if they provide sufficient
vertical resolution. Field survey data may be necessary to actually
determine elevations of the gutter or channel. If the slope varies
significantly throughout the system, this may become a determining
factor in how the system is broken into subcollectors, collectors, and
main channel.
210
Chapter 7 Subbasin Elements
211
Chapter 7 Subbasin Elements
Figure 141. Channel editor for the kinematic wave transform method.
212
Chapter 7 Subbasin Elements
ModClark Transform
The ModClark method is a linear, quasi-distributed transform method
that is based on the Clark conceptual unit hydrograph. It
fundamentally represents the subbasin as a collection of grid cells.
The Clark method uses a time-area curve and the time of
concentration to develop a translation hydrograph. By contrast, the
ModClark method eliminates the time-area curve and instead uses a
separate travel time index for each grid cell. The travel time index for
each cell is scaled by the overall time of concentration. Excess
precipitation falling on each grid cell is lagged by the scaled time index
and then routed through a linear reservoir. The outputs from the linear
reservoirs of the cells are combined to produce the final hydrograph.
External utility programs are used to create the grid cell information at
one of the following cell sizes: 10, 20, 50, 100, 200, 500, 1000, 2000,
5000, 10000 meters. The Component Editor is shown in Figure 142.
213
Chapter 7 Subbasin Elements
The percentage of unit runoff occurring before the peak flow is not
uniform across all watersheds because it depends on flow length,
ground slope, and other properties of the watershed. By changing the
percentage of unit runoff before the peak, alternate unit hydrographs
can be computed for watersheds with varying topography and other
conditions that effect runoff. The percentage of runoff occurring before
the peak is reflected in the peak rate factor (PRF). It has been found
that flat watersheds typically have a lower PRF that may be as small
as 100. Steeper watersheds have a larger PRF thay may range up to
600. The default unit hydrograph has a PRF of 484. Unit hydrographs
with specific peak rate factors are defined in the National Engineering
Handbook (NRCS, 2007). A very specific case of the unit hydrograph
for flat coastal watersheds is known as the Delmarva unit hydrograph
(Welle, Woodward, and Moody, 1980).
The standard lag is defined as the length of time between the centroid
of precipitation mass and the peak flow of the resulting hydrograph.
Examination of the equations used in deriving the curvilinear unit
hydrograph show that the lag time can be computed as the duration of
unit precipitation divided by two, plus 60% of the time of concentration.
214
Chapter 7 Subbasin Elements
The standard lag is defined as the length of time between the centroid
of precipitation mass and the peak flow of the resulting hydrograph.
Many relationships for estimating lag from subbasin characteristics
have been developed for different regions.
The peaking coefficient measures the steepness of the hydrograph
that results from a unit of precipitation. It ranges typically from 0.4 to
0.8 with lower values associated with steep-rising hydrographs. It is
estimated using the best judgement of the user, or possibly from
locally-developed relationships to watershed physical features.
Ft Worth District Method
The Ft Worth method is the result of regional urban studies conducted
by the U.S. Army Corps of Engineers, Ft. Worth District in the Dallas-Ft
215
Chapter 7 Subbasin Elements
Worth, Austin, and San Antonio areas. The method is applicable for
watersheds substantially similar to the watersheds included in the
regional study. The user enters values for five physical characteristics
of a subbasin. Equations from the regional study use the parameter
values to calculate the Snyder standard lag. The Component Editor is
shown in Figure 145.
The total length is measured from the outlet of the subbasin, along the
main watercourse to the most hydraulically remote point on the
subbasin boundary.
The centroid length is measured from the outlet of the subbasin, along
the main watercourse, to a point opposite the subbasin centroid.
The weighted slope is the slope of the main watercourse between
points located at 10 percent and 85 percent of the length of the main
watercourse, measured from the subbasin outlet.
216
Chapter 7 Subbasin Elements
The Fort Worth method does not include estimation of the peaking
coefficient. This parameter must be estimated by another method and
entered directly.
Tulsa District Method
The Tulsa method is the result of regional urban studies conducted by
the U.S. Army Corps of Engineers, Tulsa District using watersheds
across Oklahoma. The method is applicable for watersheds
substantially similar to the watersheds included in the regional study.
The user enters values for four physical characteristics of a subbasin.
The values are used with equations developed in the regional study for
estimating the standard lag and peaking coefficient. The Component
Editor is shown in Figure 146.
Figure 146. Snyder unit hydrograph transform method editor showing Tulsa
District parameter estimation.
The total length is measured from the outlet of the subbasin, along the
main watercourse to the most hydraulically remote point on the
subbasin boundary.
The centroid length is measured from the outlet of the subbasin, along
the main watercourse, to a point opposite the subbasin centroid.
The weighted slope is the slope of the main watercourse between
points located at 10 percent and 85 percent of the length of the main
watercourse, measured from the subbasin outlet.
Channelization percentage is estimated as the percentage of the
subbasin where the drainage system has been improved to more
efficiently carry runoff to the outlet. Such improvements may include
street gutters, straightened channels, concrete lined channels, or other
effeciencies. This parameter provides a way to account for the effects
of urbanization on the hydrologic response of the subbasin.
217
Chapter 7 Subbasin Elements
218
Chapter 7 Subbasin Elements
219
Chapter 7 Subbasin Elements
220
Chapter 7 Subbasin Elements
221
Chapter 7 Subbasin Elements
222
Chapter 7 Subbasin Elements
223
Chapter 7 Subbasin Elements
224
Chapter 7 Subbasin Elements
225
Chapter 7 Subbasin Elements
Recession Baseflow
The recession baseflow method is designed to approximate the typical
behavior observed in watersheds when channel flow recedes
exponentially after an event. This method is intended primarily for
event simulation. However, it does have the ability to automatically
reset after each storm event and consequently may be used for
continuous simulation. It does not conserve mass within the subbasin.
The Component Editor is shown in Figure 155.
The initial baseflow at the beginning of a simulation must be specified.
Two methods are available for specifying the initial condition: initial
discharge and initial discharge per area. Using the first method, you
must specify the initial baseflow as a discharge with units of volume
per time. This method is particularly good when there is observed
streamflow data at the outlet of the subbasin for determining the initial
flow in the channel. In the second method you specify the initial
baseflow as a volume per area per time. This method is better suited
when general guidelines for watershed yield must be used to estimate
the initial flow.
226
Chapter 7 Subbasin Elements
References
Natural Resources Conservation Service (NRCS). 2007. National
Engineering Handbook: Chapter 16 Hydrographs. Washington, DC.
Welle, P., D. Woodward, and H. Moody. 1980. "A Dimensionless Unit
Hydrograph for the Delmarva Peninsula." American Society of
Agricultural Engineers Paper 80-2013, St. Joseph, Michigan.
227
Chapter 7 Subbasin Elements
228
Chapter 8 Reach Elements
CHAPTER 8
Reach Elements
A reach is an element with one or more inflow and only one outflow.
Inflow comes from other elements in the basin model. If there is more
than one inflow, all inflow is added together before computing the
outflow. Outflow is computed using one of the several available
methods for simulating open channel flow. Optionally it may include a
method for representing interactions with the subsurface. The reach
element can be used to model rivers and streams.
229
Chapter 8 Reach Elements
Figure 156. Reach component editor. All element editors include the basin
model and element name, description, and downstream
connection. The reach editor also has the routing and loss/gain
method selections.
230
Chapter 8 Reach Elements
The slope should be the average slope for the whole reach. This value
may be estimated from maps if they provide sufficient vertical
resolution. Field survey data may be necessary to actually determine
elevations of the channel bed for calculating the slope. If the slope
varies significantly throughout the stream represented by the reach, it
may be necessary to use multiple reaches with different slopes.
The Manning's n roughness coefficient should be the average value for
the whole reach. This value can be estimated from pictures of streams
with known roughness coefficient (Barnes, 1967) or by calibration.
The number of subreaches is used as a hint to the program when it
determines the correct distance step to use during routing calculations.
Criteria based on steepness of the inflow hydrograph and other factors
are used to automatically determine the correct distance and time
steps for solving the kinematic wave equation. The default value is 2
but may be optionally increased.
Five options are provided for specifying the cross section shape: circle,
deep, rectangle, trapezoid, and triangle. The circle shape cannot be
used for pressure flow or pipe networks, but is suitable for representing
a free water surface inside a pipe. The deep shape should only be
used for flow conditions where the flow depth is approximately equal to
the flow width. Depending on the shape you choose, additional
information will have to be entered to describe the size of the cross
section shape. This information may include a diameter (circle) bottom
width (deep, rectangle, trapezoid) or side slope (trapezoid, triangle).
The side slope is dimensionless and entered as the units of horizontal
distance per one unit of vertical distance.
231
Chapter 8 Reach Elements
Lag Routing
The lag routing method only represents the translation of flood waves.
It does not include any representation of attenuation or diffusion
processes. Consequently, it is best suited to short stream segments
with a predicable travel time that doesn't vary with flow depth. The
Component Editor is shown in Figure 158.
The only parameter is the lag time in minutes. Inflow to the reach is
delayed in time by an amount equal to the specified lag, and then
becomes outflow.
232
Chapter 8 Reach Elements
Figure 159. Modified Puls routing method editor. The optional elevation-
discharge function is disabled because no curves are currently
available for selection. The invert can be entered after a curve is
selected.
233
Chapter 8 Reach Elements
Muskingum Routing
The Muskingum routing method uses a simple conservation of mass
approach to route flow through the stream reach. However, it does not
assume that the water surface is level. By assuming a linear, but non-
level, water surface it is possible to account for increased storage
during the rising side of a flood wave and decreased storage during the
falling side. By adding a travel time for the reach and a weighting
between the influence of inflow and outflow, it is possible to
approximate attenuation. The Component Editor is shown in Figure
160.
Muskingum-Cunge Routing
The Muskingum-Cunge routing method is based on the combination of
the conservation of mass and the diffusion representation of the
conservation of momentum. It is sometimes referred to as a variable
coefficient method because the routing parameters are recalculated
every time step based on channel properties and the flow depth. It
represents attenuation of flood waves and can be used in reaches with
a small slope. The Component Editor is shown in Figure 161.
234
Chapter 8 Reach Elements
The time interval selection provides two options. The program can
automatically select a fixed time interval that maintains numeric
stability during the most steeply rising hydrograph. Alternately, the
program can automatically vary the time interval to take as long a time
interval as possible while also maintaining numeric stability.
The length should be the total length of the reach element. Usually
this value will be measured from maps of the watershed.
The slope should be the average slope for the whole reach. This value
may be estimated from maps if they provide sufficient vertical
resolution or by field surveying. If the slope varies significantly
throughout the stream represented by the reach, it may be necessary
to use multiple reaches with different slopes.
The Manning's n roughness coefficient should be the average value for
the whole reach. This value can be estimated from pictures of streams
with known roughness coefficient (Barnes, 1967) or by calibration.
Optionally, an invert elevation may be specified. When used, the flow
depth computed during the routing is added to the invert elevation to
compute the stage.
Five options are provided for specifying the cross section shape: circle,
eight point, rectangle, trapezoid, and triangle. The circle shape cannot
be used for pressure flow or pipe networks, but is suitable for
representing a free water surface inside a pipe. Depending on the
shape you choose, additional information will have to be entered to
describe the size of the cross section shape. This information may
include a diameter (circle) bottom width (rectangle, trapezoid) or side
235
Chapter 8 Reach Elements
236
Chapter 8 Reach Elements
program provides a different level of detail and not all methods are
equally adept at representing a particular stream. Further, because of
differing data requirements, some loss/gain methods are only
compatible with certain routing methods.
The loss/gain method for a reach is selected on the Component Editor
for the reach element. Access the Component Editor by clicking the
reach element icon on the "Components" tab of the Watershed
Explorer (Figure 156). You can also access the Component Editor by
clicking on the element icon in the basin map, if the map is currently
open. You can select a loss/gain method from the list of three
available choices. If you choose the None method, the reach will
perform routing calculations without including any losses or gains to
the channel. Use the selection list to choose the method you wish to
use. Each reach may use a different method or several reaches may
use the same method.
When a new reach is created, it is automatically set to use the default
loss/gain method specified in the program settings. You may change
the loss/gain method for a reach at any time using the Component
Editor for the reach element. Since a reach can only use one loss/gain
method at a time, you will be warned when changing methods that the
old parameter data will be lost. You can turn off this warning in the
program settings. You can change the loss/gain method for several
reaches simultaneously. Click on the Parameters menu and select
the Loss/Gain Change Method command. The loss/gain method
you choose will be applied to the selected reaches in the basin model,
or to all reaches if none are currently selected.
The parameters for each loss/gain method are presented on a
separate Component Editor from the reach element editor. The
"Loss/Gain" editor is always shown next to the "Routing" editor. The
information shown on the loss/gain editor will depend on which method
is currently selected.
Constant Loss/Gain
The constant loss/gain method uses an empirical relationship to
calculate channel loss using a fixed flow rate reduction and a ratio of
the flow. It does not include any capability to representing gaining
streams. A fixed flow rate is subtracted from the routed flow and then
the remainder is multiplied by a ratio. The reduced flow becomes the
outflow for the reach. This method is compatible with all routing
methods. The Component Editor is shown in Figure 163.
The flow rate parameter specifies the amount of flow to be subtracted
from the inflow. It may be zero or greater than zero; if the specified
value is zero then no flow rate reduction will occur. If the specified
237
Chapter 8 Reach Elements
value is greater than the inflow for a time step, then no outflow will
result for that interval.
The fraction is used to reduce the inflow in linear proportion to the flow
rate. After the constant value is subtracted from the routed inflow, the
remainder is multiplied by the value one minus the fraction. The
fraction must be between zero and one.
Percolation Loss/Gain
The percolation method uses a constant infiltration rate in combination
with the inundated area in the reach to compute channel loss. It does
not include any capability to representing gaining streams. This
method is only compatible with the modified Puls and Muskingum-
Cunge routing methods. When used with the modified Puls method,
the optional elevation-discharge function and invert elevation must be
specified in the routing parameters. The inundated area is calculated
differently depending on the routing method. In combination with the
modified Puls method, the current storage and outflow are combined
with the elevation-discharge curve to compute inundated area. With
the Muskingum-Cunge method, the wetted perimeter is multiplied by
the reach length to compute area. The Component Editor is shown in
Figure 164.
The rate parameter specifies percolation in terms of a flow rate per
area. The inundated area is multiplied by the rate to determine the
channel loss for each time interval.
238
Chapter 8 Reach Elements
239
Chapter 8 Reach Elements
240
Chapter 9 Reservoir Elements
CHAPTER 9
Reservoir Elements
A reservoir is an element with one or more inflow and one computed
outflow. Inflow comes from other elements in the basin model. If there
is more than one inflow, all inflow is added together before computing
the outflow. It is assumed that the water surface in the reservoir pool
is level. Several methods are available for defining the storage
properties of the reservoir. The element can be used to model
reservoirs, lakes, and ponds. Access the Component Editor by clicking
the reservoir element icon on the "Components" tab of the Watershed
Explorer (Figure 165).
Figure 165. Reservoir component editor using the outflow curve routing
method with a storage-discharge function.
241
Chapter 9 Reservoir Elements
Storage Method
There are three different options for specifying the storage relationship.
The first option is the Storage-Discharge choice. The user must
select a storage-discharge curve from the available curves in the
Paired Data Manager. The second option is the Elevation-
Storage-Discharge choice. The user must select both a storage-
discharge curve and elevation-storage curve from the Paired Data
Manager. The final option is the Elevation-Area-Discharge
choice. In this case the user must select both an elevation-area curve
and an elevation-discharge curve from the Paired Data Manager. With
this choice, the program automatically transforms the elevation-area
curve into an elevation-storage curve using the conic formula.
Regardless of which option is selected, the routing is always performed
using only the storage-discharge curve. After the routing is complete
using the storage-discharge curve, the program will compute the
elevation and surface area for each time step, depending on the
selected storage method.
242
Chapter 9 Reservoir Elements
Initial Condition
The initial condition sets the amount of storage in the reservoir at the
beginning of a simulation. Therefore, the simplest option is to specify
the Storage as a volume of water in the reservoir. For convenience,
other options are also provided. The Inflow=Outflow method takes
the inflow to the reservoir at the beginning of the simulation, and uses
the storage-discharge curve to determine the storage required to
produce that same flowrate as the outflow from the reserovir. Some
storage methods permit the specification of Elevation as the initial
condition. In this case, the elevation provided by the user is used to
interpolate a storage value from the elevation-storage curve. Other
storage methods permit the specification of Discharge as the initial
condition. In this case, the storage is interpolated from the storage-
discharge curve. The initial condition options depend on the selected
storage method and are shown in Table 21.
243
Chapter 9 Reservoir Elements
Table 21. Available initial condition options for different storage methods
used with the outflow curve routing method.
Storage Method
There are two different options for specifying the storage relationship,
as shown in Figure 166. The first option is the Elevation-Storage
choice. The user must select an elevation-storage curve from the
available curves in the Paired Data Manager. After the routing is
complete, the program will compute the elevation and storage for each
time interval. The second option is the Elevation-Area choice,
which requires the selection of an elevation-area curve from the
available curves in the Paired Data Manager. With this choice, the
program automatically transforms the elevation-area curve into an
elevation-storage curve using the conic formula. After the routing is
complete, the program will compute the elevation, surface area, and
storage for each time interval.
You must select an appropriate function to define the selected storage
method. For example, if you select the Elevation-Storage method
you must select an appropriate elevation-storage paired data function
that defines the storage characteristics of the reservoir. The
appropriate selection list will be shown directly under the storage
method selection list. The necessary paired data functions must be
defined in the paired data manager before they can be used in the
reservoir. Choose an appropriate function in each selection list. If you
wish, you can use a chooser by clicking the paired data button next to
the selection list. A chooser will open that shows all of the paired data
functions of that type. Click on a function to view its description.
244
Chapter 9 Reservoir Elements
Initial Condition
The initial condition sets the amount of storage in the reservoir at the
beginning of a simulation. The choices depend on the method
selected for specifying the storage characteristics of the reservoir.
When the Elevation-Storage method is selected, you may choose
to specify the initial elevation or the initial storage. When the
Elevation-Area method is selected, you must specify an initial
elevation.
Figure 166. Reservoir component editor using the specified release routing
method with an elevation-storage function.
245
Chapter 9 Reservoir Elements
Storage Method
There are two different options for specifying the storage relationship.
The first option is the Elevation-Storage choice, as shown in
Figure 167. The user must select an elevation-storage curve from the
available curves in the Paired Data Manager. After the routing is
complete, the program will compute the elevation and storage for each
time interval. The second option is the Elevation-Area choice,
246
Chapter 9 Reservoir Elements
Figure 167. Reservoir component editor using the outflow structures routing
method with three outlets, one spillway, and one dam top
representing an emergency spillway.
247
Chapter 9 Reservoir Elements
the selection list. A chooser will open that shows all of the paired data
functions of that type. Click on a function to view its description.
Initial Condition
The initial condition sets the amount of storage in the reservoir at the
beginning of a simulation. Therefore, the simplest option is to specify
the Storage as a volume of water in the reservoir. For convenience,
other options are also provided. The Inflow=Outflow method takes
the inflow to the reservoir at the beginning of the simulation, and
determines the pool elevation necessary to cause that outflow through
the outlet structures. The pool elevation is then used in the elevation-
storage curve to determine the matching storage. The pool
Elevation method can also be selected for the initial condition. In
this case, the elevation provided by the user is used to interpolate a
storage value from the elevation-storage curve. The initial condition
options depend on the selected storage method and are shown in
Table 22.
Table 22. Available initial condition options for different storage methods
used with the outlet structures routing method.
Tailwater Method
The selected tailwater method determines how submergence will be
calculated for the individual structures specified as part of the
reservoir. When a structure is submerged, the discharge through the
structure will decrease in accordance with the physics of the structure
and the tailwater elevation for each time interval. Only one tailwater
method can be selected and it is applied to all structures specified as
part of the reservoir.
248
Chapter 9 Reservoir Elements
clicking the paired data button next to the selection list. The rating
curve should be specified in the same vertical datum as the function
used to describe the storage characteristics of the reservoir.
249
Chapter 9 Reservoir Elements
Outlets
Outlets can only be included in reservoirs using the Outflow
Structures routing method. Outlets typically represent structures
near the bottom of the dam that allow water to exit in a controlled
manner. They are often called gravity outlets because they can only
move water when the head in the reservoir is greater than the head in
the tailwater. Up to 10 independent outlets can be included in the
reservoir. Select the number of outlets you wish to include. An icon
for each outlet will be added to the reservoir icon in the Watershed
Explorer. You will need to click on the individual outlet icon to enter
250
Chapter 9 Reservoir Elements
parameter data for it. There are two different methods for computing
outflow through an outlet: culvert or orifice.
Culvert Outlet
The culvert outlet allows for partially full or submerged flow through a
culvert with a variety of cross-sectional shapes. It can account for inlet
controlled outflow or outlet control. A typical culvert outlet is shown in
Figure 168.
You must select a solution method for the culvert: inlet control, outlet
control, or automatic. You may select Inlet Control if it is known
that at all times during a simulation the culvert outflow will be controlled
by a high pool elevation in the reservoir. You may likewise select
Outlet Control if it is known that at all times the culvert outflow will
be controlled by a high tailwater condition. In general it is best to
select Automatic control and the program will automatically
determine the controlling inlet or outlet condition.
You must select the number of identical barrels. This can be used to
specify several culvert outlets that are identical in all parameters.
There can be up to 10 identical barrels.
251
Chapter 9 Reservoir Elements
The length of the culvert must be specified. This should be the overall
length of the culvert including any projection at the inlet or outlet.
The inlet elevation must be specified as the invert elevation at the
bottom of the culvert on the inlet side. The inlet side is always
assumed to be in the reservoir pool. This should be measured in the
same vertical datum as the paired data functions defining the storage
characteristics of the reservoir.
The entrance coefficient describes the energy loss as water moves into
the inlet of the culvert. Values may range from 0.2 up to 1.0.
The exit coefficient describes the energy loss that occurs when water
expands as it leaves the culvert outlet. Typically the value is 1.0.
The outlet elevation must be specified as the invert elevation at the
bottom of the culvert on the outlet side. The outlet side is always
assumed to be in the reservoir tailwater. This should be measured in
252
Chapter 9 Reservoir Elements
the same vertical datum as the paired data functions defining the
storage characteristics of the reservoir.
A Manning's n value should be entered that describes the roughness in
the culvert. At this time, the same n value must be used for the entire
length of the culvert, as well as the entire top, sides, and bottom.
Orifice Outlet
The orifice outlet assumes sufficient submergence on the outlet for
orifice flow conditions to dominate. It should not be used to represent
an outlet that may flow only partially full. The inlet of the structure
should be submerged at all times by a depth at least 0.2 times the
diameter. A typical orifice outlet editor is shown in Figure 169.
You must select the number of identical barrels. This can be used to
specify several culvert outlets that are identical in all parameters.
There can be up to 10 identical barrels.
The center elevation specifies the center of the cross-sectional flow
area. This should be measured in the same vertical datum as the
paired data functions defining the storage characteristics of the
reservoir. It is used to compute the head on the outlet, so no flow will
be released until the reservoir pool elevation is above this specified
elevation.
The cross-sectional flow area of the outlet must be specified. The
orifice assumptions are independent of the shape of the flow area.
The dimensionless discharge coefficient must be entered. This
parameter describes the energy loss as water exits the reservoir
through the outlet.
253
Chapter 9 Reservoir Elements
Spillways
Spillways can only be included in reservoirs using the Outflow
Structures routing method. Spillways typically represent structures
at the top of the dam that allow water to go over the dam top in a
controlled manner. Up to 10 independent spillways can be included in
the reservoir. Select the number of spillways you wish to include. An
icon for each spillway will be added to the reservoir icon in the
Watershed Explorer. You will need to click on the individual outlet icon
to enter parameter data for it. There are three different methods for
computing outflow through a spillway: broad-crested, ogee, and user
specified. The broad-crested and ogee methods may optionally
include gates. If no gates are selected, then flow over the spillway is
unrestricted. When gates are included, the flow over the spillway will
be controlled by the gates. Up to 10 independent gates may be
included on a spillway.
Broad-Crested Spillway
The broad-crested spillway allows for controlled flow over the top of the
reservoir according to the weir flow assumptions. A typical broad-
crested spillway editor is shown in Figure 170.
The crest elevation of the spillway must be specified. This should be
measured in the same vertical datum as the paired data functions
defining the storage characteristics of the reservoir.
The length of the spillway must be specified. This should be the total
width through which water passes.
The discharge coefficient accounts for energy losses as water enters
the spillway, flows through the spillway, and eventually exits the
spillway. Depending on the exact shape of the spillway, typical values
254
Chapter 9 Reservoir Elements
255
Chapter 9 Reservoir Elements
The rating curve describing flow over the spillway must be selected.
Before it can be selected it must be created in the Paired Data
Manager as an elevation-discharge function. The function must be
calculated external to the program on the basis of advanced spillway
hydraulics or experimentation.
256
Chapter 9 Reservoir Elements
Spillway Gates
Spillway gates are an optional part of specifying the configuration of a
spillway. They may be included on either broad-crested or ogee
spillways. The number of gates to use for a spillway is specified on the
spillway editor (Figure 170 and Figure 171). An icon for each gate will
be added to the spillway icon under the reservoir icon in the Watershed
Explorer. You will need to click on the individual gate icon to enter
parameter data for it. There are two different methods for computing
outflow through a gated spillway: sluice or radial. In both cases you
may specify the number of identical units; each identical unit has
exactly the same parameters, including how the gate is controlled.
Sluice Gate
A sluice gate moves up and down in a vertical plane above the spillway
in order to control flow. The water passes under the gate as it moves
over the spillway. For this reason it is also called a vertical gate or
underflow gate. The editor is shown in Figure 173.
257
Chapter 9 Reservoir Elements
Radial Gate
A radial gate rotates above the spillway with water passing under the
gate as it moves over the spillway. This type of gate is also known as
a tainter gate. The editor is shown in Figure 174.
258
Chapter 9 Reservoir Elements
Dam Tops
Dam tops can only be included in reservoirs using the Outflow
Structures routing method. These typically represent the top of the
dam, above any spillways, where water goes over the dam top in an
uncontrolled manner. In some cases it may represent an emergency
spillway. Up to 10 independent dam tops can be included in the
reservoir. Select the number of dam tops you wish to include. An icon
for each dam top will be added to the reservoir icon in the Watershed
Explorer. You will need to click on the individual dam top icon to enter
parameter data for it. There are two different methods for computing
outflow through a dam top: level or non-level.
Level Dam Top
The level dam top assumes flow over the dam can be represented as a
broad-crested weir. The calculations are essentially the same as for a
broad-crested spillway. They are included separately mostly for
conceptual representation of the reservoir structures. A typical level
dam top is shown in Figure 175.
259
Chapter 9 Reservoir Elements
Figure 175. Dam top editor with the level overflow method selected.
The crest elevation of the dam top must be specified. This should be
measured in the same vertical datum as the paired data functions
defining the storage characteristics of the reservoir.
The length of the dam top must be specified. This should be the total
width through which water passes, excluding any amount occupied by
spillways if any are included.
The discharge coefficient accounts for energy losses as water
approaches the dam top and flows over the dam. Depending on the
exact shape of the dam top, typical values range from 1.10 to 1.66 in
System International units (2.0 to 3.0 US Customary units).
Non-Level Dam Top
The non-level dam top assumes that the top of the dam can be well-
represented by a cross section with eight station-elevation pairs. A
separate flow calculation is carried out for each segment of the cross
section. The broad-crested weir assumptions are made for each
segment. A typical non-level dam top is shown in Figure 176.
A cross section must be selected which describes the shape of the top
of the dam with a simplified eight point shape. From abutment to
abutment of the dam, but should not include any spillways that may be
included. It may be necessary to use multiple dam tops to represent
the different sections of the dam top between spillways. The cross
section should extend from the dam top up to the maximum water
surface elevation that will be encountered during a simulation. The
cross section must be defined in the paired data manager before it can
be used in the reservoir element.
The discharge coefficient accounts for energy losses as water
approaches the dam top and flows over the dam. The same value is
used for all segments of the dam top. Typical values range from 2.6 to
4.0 depending on the exact shape of the dam top.
260
Chapter 9 Reservoir Elements
Figure 176. Dam top editor with the non-level overflow method selected.
Pumps
Pumps can only be included in reservoirs using the Outflow
Structures routing method. These typically represent pumps in
interior ponds or pump stations that are intended to move water out of
the reservoir and into the tailwater when gravity outlets alone cannot
move sufficient water. Up to 10 independent pumps can be included in
the reservoir. Select the number of pumps you wish to include. An
icon for each pump will be added to the reservoir icon in the
Watershed Explorer. You will need to click on the individual pump icon
to enter parameter data for it. There is only one method for computing
outflow through a pump: head-discharge pump.
Head-Discharge Pump
The head-discharge pump is designed to represent low-head, high-
flow situations. This means that the pump is designed for high flow
rates against a relatively small head. The pump can be controlled to
come on and shut off as the reservoir pool elevation changes. A
typical head-discharge pump is shown in Figure 177.
The number of identical units must be specified. This allows data to be
entered only once when there are multiple pump units with exactly the
same parameters.
261
Chapter 9 Reservoir Elements
The intake elevation defines the elevation in the reservoir pool where
the pump takes in water. This should be measured in the same
vertical datum as the paired data functions defining the storage
characteristics of the reservoir.
The line elevation defines the highest elevation in the pressure line
from the pump to the discharge point. This should be measured in the
same vertical datum as the paired data functions defining the storage
characteristics of the reservoir.
You must specify the elevation when the pump turns on. This should
be measured in the same vertical datum as the paired data functions
defining the storage characteristics of the reservoir. Once the pump
turns on, it will remain on until the reservoir pool elevation drops below
the trigger elevation to turn the pump off.
You must specify the elevation when the pump turns off. This should
be measured in the same vertical datum as the paired data functions
defining the storage characteristics of the reservoir. This elevation
must be lower than the elevation at which the pump turns on.
The specification of a minimum rest time is optional. If it is used, once
a pump shuts off it must remain off the specified minimum rest time
even if the reservoir pool elevation reaches the trigger elevation to turn
the pump on.
The specification of a minimum run time is optional. If it is used, once
a pump turns on it must remain on the specified minimum run time
262
Chapter 9 Reservoir Elements
even if the reservoir pool elevation drops below the trigger elevation to
turn the pump off. The only exception is if the pool elevation drops
below the intake elevation, then the pump will shut off even though the
minimum run time is not satisfied.
The equipment loss includes all energy losses between the intake and
discharge points, including the pump itself. This loss is added to the
head difference due to reservoir pool elevation and tailwater elevation
to determine the total energy against which the pump must operate.
The head-discharge curve describes the capacity of the pump as a
function of the total head. Total head is the head difference due to
reservoir pool elevation and tailwater elevation, plus equipment loss. A
curve must be defined as an elevation-discharge function in the paired
data manager before it can be selected for a pump in the reservoir.
You can press the paired data button next to the selection list to use a
chooser. The chooser shows all of the available elevation-discharge
functions in the project. Click on a function to view its description.
Dam Break
Dam break can only be included in reservoirs using the Outflow
Structures routing method. Only one dam break can be included in
the reservoir. Choose whether you wish to include dam break. An
icon for the dam break will be added to the reservoir icon in the
Watershed Explorer. You will need to click on the dam break icon to
enter parameter data for it. There are two different methods for
computing outflow through a dam break: overtop and piping.
Overtop Dam Break
The overtop dam break (Figure 178) is designed to represent failures
caused by overtopping of the dam. These failures are most common
in earthen dams but may also occur in concrete arch, concrete gravity,
or roller compacted dams as well. The failure begins when
appreciable amounts of water begin flowing over or around the dam
face. The flowing water will begin to erode the face of the dam. Once
erosion begins it is almost impossible to stop the dam from failing. The
method begins the failure at a point on the top of the dam and expands
it in a trapezoidal shape until it reaches the maximum size. Flow
through the expanding breach is modeled as weir flow.
The top elevation is the top of the dam face. The breach may be
initiated at a lower elevation than the top depending on the selection of
the trigger. This information is used to constrain the top of the breach
opening as it grows. It should be measured in the same vertical datum
as the paired data functions defining the storage characteristics of the
reservoir.
263
Chapter 9 Reservoir Elements
Figure 178. Dam break editor with the overtop breach method selected.
264
Chapter 9 Reservoir Elements
failure will begin. For the Specific Time method, the breach will
begin opening at the specified time regardless of the reservoir pool
elevation. When specifying an elevation, it should be measured in the
same vertical datum as the paired data functions defining the storage
characteristics of the reservoir.
The progression method determines how the breach grows from
initiation to maximum size during the development time. Select the
Linear method to have the breach grow in equal increments of depth
and width. Select the Sine Wave method to have the breach grow
quickly in the early part of breach development and more slowly as it
reaches maximum size. The speed varies according to the first quarter
cycle of a since wave. Select the User Curve method to have the
breach grow according to a specified pattern. You will need to select a
curve in the selection list, which will show all percentage curves
defined in the paired data manager. The independent variable should
range from 0 to 100 percent and define the percentage of the
development time. The dependent variable should define the
percentage opening of the maximum breach size. The function must
be monotonically increasing.
Piping Dam Break
The piping dam break is designed to represent failures caused by
piping inside the dam. These failures typically occur only in earthen
dams. The failure begins when water naturally seeping through the
dam core increases in velocity and quantity enough to begin eroding
fine sediments out of the soil matrix. If enough material erodes, a
direct piping connect may be established from the reservoir water to
the dam face. Once such a piping connect is formed it is almost
impossible to stop the dam from failing. The method begins the failure
at a point in the dam face and expands it as a circular opening. When
the opening reaches the top of the dam, it continues expanding as a
trapezoidal shape. Flow through the circular opening is modeled as
orifice flow while in the second stage it is modeled as weir flow.
The piping dam break (Figure 179) uses many of the same parameters
as the overtop dam break. The top elevation, bottom elevation, bottom
width, left slope, and right slope all are used to describe a trapezoidal
breach opening that will be the maximum opening in the dam. These
are only used once the piping opening transitions to an open breach.
The parameters for development time, trigger method, and progression
method are also the same for defining when the failure initiates, how
long it takes to attain maximum breach opening, and how the breach
develops during the development time. The remaining parameters,
unique to piping dam break, are described below.
265
Chapter 9 Reservoir Elements
Figure 179 Dam break editor with the piping breach method selected.
The piping elevation indicates the point in the dam where the piping
failure first begins to form. This should be measured in the same
vertical datum as the paired data functions defining the storage
characteristics of the reservoir.
The piping coefficient is used to model flow through the piping opening
as orifice flow. As such, the coefficient represents energy losses as
water moves through the opening.
Dam Seepage
Dam seepage can only be included in reservoirs using the Outflow
Structures routing method. Most dams have some water seeping
through the face of the dam. The amount of seepage depends on the
elevation of water in the dam, the elevation of water in the tailwater,
the integrity of the dam itself, and possibly other factors. In some
situations, seepage from the pool through the dam and into the
tailwater can be a significant source of discharge that must be
modeled. Interior ponds may discharge seepage water but in some
situations water in the main channel may seep through the levee or
dam face and enter the pool. Both of these situations can be
represented using the dam seepage structure.
There can only be one dam seepage structure in a reservoir that must
represent all sources and sinks of seepage. When water seeps out of
266
Chapter 9 Reservoir Elements
Additional Release
An additional release can only be included in reservoirs using the
Outflow Structures routing method. In most situations a dam can
be properly configured using various outlet structures such as
spillways, outlets, etc. The total outflow from the reservoir can be
calculated automatically using the physical properties entered for each
of the included structures. However, some reservoirs may have an
additional release beyond what is represented by the various physical
structures. In many cases this additional release is a schedule of
managed releases achieved by operating spillway gates. The
additional release can be used in combination with other outlet
structures to determine the total release from the reservoir.
The additional release that will be specified must be stored as a
discharge gage. The appropriate gage can be selected in the editor as
267
Chapter 9 Reservoir Elements
Evaporation
Evaporation can only be included in reservoirs using the Outflow
Structures routing method. Additionally, the reservoir must be set
to use the elevation-area storage option. Water losses due to
evaporation may be an important part of the water balance for a
reservoir, especially in dry or desert environments. The evaporation
losses are different from other structures because they do not
contribute to either main or auxiliary outflow. They are accounted
separately and available for review with the other time-series results
for the reservoir. An evaporation depth is computed for each time
interval and then multiplied by the current surface area.
Monthly Evaporation
The monthly evaporation method can be used to specify a separate
evaporation rate for each month of the year, entered as a total depth
for the month. The evaporation data must be developed through
separate, external analysis and entered as shown in Figure 182.
268
Chapter 9 Reservoir Elements
Figure 182. Evaporation editor showing the monthly evaporation method for
a reservoir.
269
Chapter 9 Reservoir Elements
270
Chapter 10 Source, Junction, Diversion, and Sink Elements
CHAPTER 10
Source
A source is an element with no inflow, one outflow, and is one of only
two ways to produce flow in the basin model. The source can be used
to represent boundary conditions to the basin model such as measured
outflow from reservoirs or unmodeled headwater regions. The
boundary condition can also represent contributing area modeled in a
separate basin model. Access the Component Editor by clicking the
source element icon on the "Components" tab of the Watershed
Explorer (Figure 183).
Figure 183. Source component editor. All element editors include the basin
model and element name, description, and downstream
connection. The source editor also has an area and other
properties depending on the method.
Representative Area
Specification of a representative area for the source is optional. If the
area is specified, then it will be possible for the program to
271
Chapter 10 Source, Junction, Diversion, and Sink Elements
Discharge Gage
When the Discharge Gage method is selected for the source, you
must select a time-series discharge gage. The gage should record the
discharge to use for each time interval during a simulation. If there is
missing data in the record and the basin model options are set to
replace missing data, a zero flow rate will be substituted for each
missing data value. If the basin model is not set to replace missing
data, any missing data will cause the simulation to stop and an error
message will be displayed.
272
Chapter 10 Source, Junction, Diversion, and Sink Elements
Constant Flow
When the "Constant Flow" method is selected for the source, you must
enter a flow rate. The same flow rate is used for all time steps during a
simulation. Specify the flow rate on the Component Editor for the
source element. Access the Component Editor by clicking the source
element icon on the "Components" tab of the Watershed Explorer
(Figure 185). You can also access the Component Editor by clicking
on the element icon in the basin map, if the map is currently open.
273
Chapter 10 Source, Junction, Diversion, and Sink Elements
that varies throughout the year. The percent annual pattern must be
defined in the paired data manager before it can be used in the source
element.
Junction
A junction is an element with one or more inflows and only one outflow.
All inflow is added together to produce the outflow by assuming zero
storage at the junction. It is usually used to represent a river or stream
confluence.
The junction element does not have any special data or properties; it
only has the standard Component Editor used by all elements. Access
the Component Editor by clicking the junction element icon on the
"Components" tab of the Watershed Explorer (Figure 186). You can
also access the Component Editor by clicking on the element icon in
the basin map, if the map is currently open.
Diversion
A diversion is an element with two outflows, main and diverted, and
one or more inflows. Inflow comes from other elements in the basin
model. If there is more than one inflow, all inflow is added together
before computing the outflows. Six methods are available for
computing the diverted flow that will be taken out of the channel. All
flow that is not diverted becomes the main outflow. Diverted flow can
be connected to an element that is computationally downstream. The
diversion can be used to represent weirs or pump stations that divert
flow into canals or off-stream storage. The diversion element includes
optional properties for limiting the amount of diverted flow. Access the
Component Editor by clicking the diversion element icon on the
"Components" tab of the Watershed Explorer (Figure 187).
274
Chapter 10 Source, Junction, Diversion, and Sink Elements
275
Chapter 10 Source, Junction, Diversion, and Sink Elements
276
Chapter 10 Source, Junction, Diversion, and Sink Elements
Figure 188. Specifying the constant diversion flowrate without using the
optional annual adjustment pattern.
277
Chapter 10 Source, Junction, Diversion, and Sink Elements
A rating curve must be selected for the channel. The curve should
give the stage for the entire range of inflows that will occur during a
simulation. The curve must be monotonically increased. It must be
defined in the paired data manager before it can be selected.
Optionally, a rating curve may be entered for tailwater. The curve
should give the tailwater stage in the area where the diverted flow is
discharged. It is used to automatically account for submergence of the
weir. The curve must be monotonically increasing. It must be defined
in the paired data manager before it can be selected. If the optional
rating curve is not specified, then flow over the weir is computed
assuming no tailwater influence.
278
Chapter 10 Source, Junction, Diversion, and Sink Elements
A rating curve must be selected for the channel. The curve should
give the stage for the entire range of inflows that will occur during a
simulation. The curve must be monotonically increased. It must be
defined in the paired data manager before it can be selected. You can
press the paired data button next to the selection list to use a chooser.
279
Chapter 10 Source, Junction, Diversion, and Sink Elements
280
Chapter 10 Source, Junction, Diversion, and Sink Elements
The equipment loss includes all energy losses between the intake and
discharge points. These losses are are sometimes called static losses
because the do not change very much even as the water surface
elevation in the channel fluctuates. Components of the loss include
entrance losses to the suction line, losses in the pump itself, pipe
friction losses, bend losses in the line, and exit losses at the end of the
discharge line. This loss is added to the dynamic head to determine
the total head against which the pump must operate.
The head-discharge curve describes the pumping capability of the
pump as a function of the total head. The total head is the sum of the
equipment loss and the dynamic head. The dynamic is head is first
estimated as the difference between the water surface elevation in the
channel and the line elevation. If the water surface elevation is above
the line elevation, then the estimated value will be zero. Secondarily,
the estimate is adjusted for tailwater submergence. This second stage
is only necessary if an optional tailwater rating curve is specified.
When specified, the tailwater water surface elevation is compared to
the line elevation. If the tailwater exceeds the line elevation, then the
depth of submergence over the line elevation is added to the initial
estimate of the dynamic head. The head-discharge curve is used to
calculate the diverted flow given the calculated total head. A curve
must be defined as an elevation-discharge function in the paired data
manager before it can be selected for a pump in the diversion element.
You can press the paired data button next to the selection list to use a
chooser. The chooser shows all of the available elevation-discharge
functions in the project. Click on a function to view its description.
281
Chapter 10 Source, Junction, Diversion, and Sink Elements
Sink
A sink is an element with one or more inflows but no outflow. Multiple
inflows are added together to determine the total amount of water
entering the element. Sinks can be used to represent the lowest point
of an interior drainage area or the outlet of the basin model.
The sink element does not have any special data or properties; it only
has the standard Component Editor used by all elements. Access the
Component Editor by clicking the sink element icon on the
"Components" tab of the Watershed Explorer (Figure 193). You can
also access the Component Editor by clicking on the element icon in
the basin map, if the map is currently open.
282
Chapter 11 Meteorology Description
CHAPTER 11
Meteorology Description
This chapter describes how meteorology information is entered into the
program using a meteorologic model. The meteorologic model is
responsible for preparing the boundary conditions that act on the
watershed during a simulation. Consequently, a meteorologic model is
prepared for use with one or more basin models. If the basin model
contains subbasin elements, then the meteorologic model must specify
how precipitation will be generated for each subbasin. Snowmelt may
optionally be included if environmental conditions require it.
Evapotranspiration should be included if the basin model is configured
for continuous simulation using any of the following loss methods:
deficit constant, gridded deficit constant, soil moisture accounting,
gridded soil moisture accounting. Finally, when you select energy
balance methods for snowmelt or evapotranspiration it will also be
necessary to select shortwave and longwave radiation methods.
Meteorologic Models
Meteorologic models are one of the main components in a project.
The principle purpose is to prepare meteorologic boundary conditions
for subbasins. Consequently, you must create at least one basin
model before creating a meteorologic model. A meteorologic model
can be used with many different basin models. However, results
computed by the meteorologic model will be matched with the
subbasins in the basin models using the name of the subbasin. If
subbasins in different basin models have the same name, they will
both receive the same boundary conditions from the meteorologic
model. Careful naming of subbasins is necessary so that the the
correct boundary conditions are computed for each one.
283
Chapter 11 Meteorology Description
meteorologic model, press the New button. After you press the
button a window (Figure 194) will open where you can name and
describe the new meteorologic model that will be created. A default
name is provided for the new model; you can use the default or replace
it with your own choice. A description can also be entered. If the
description is long, you can press the button to the right of the
description field to open an editor. The editor makes it easier to enter
and edit long descriptions. When you are satisfied with the name and
description, press the Create button to finish the process of creating
the new meteorologic model. You cannot press the Create button if no
name is specified for the new model. If you change your mind and do
not want to create a new meteorologic model, press the Cancel button
or the X button in the upper right to return to the Meteorologic Model
Manager window.
284
Chapter 11 Meteorology Description
mind and do not want to copy the selected meteorologic model, press
the Cancel button or the X button in the upper right to return to the
Meteorologic Model Manager window.
285
Chapter 11 Meteorology Description
Figure 197. Renaming a meteorologic model. This was accessed from the
Meteorologic Model Manager.
286
Chapter 11 Meteorology Description
287
Chapter 11 Meteorology Description
288
Chapter 11 Meteorology Description
Figure 201. Meteorologic model component editor for selecting the radiation,
precipitation, evapotranspiration, and snowmelt methods.
289
Chapter 11 Meteorology Description
Precipitation Method
Precipitation includes all rain and snowfall over the land surface. The
determination of liquid or frozen precipitation is performed in the
snowmelt method; the precipitation method only determines the
amount of water falling to the land surface. There are nine different
precipitation methods including an option for no precipitation. If you
plan to use the meteorologic model with basin models that contain
subbasins, you must choose a precipitation method. If any of the
subbasins use the ModClark gridded transform method you must use
the gridded precipitation method. You can use the "no precipitation"
option if the basin models do not contain subbasins.
Select the precipitation method in the Component Editor for the
meteorologic model (Figure 201). Only one precipitation method can
be selected at a time. You can always change the precipitation
method. When you change the precipitation method, all data for the
old method is deleted and cannot be retrieved or undeleted.
290
Chapter 11 Meteorology Description
Evapotranspiration Method
The meteorologic model is responsible for computing the potential
evapotranspiration over the land surface. The potential
evapotranspiration is the amount of evapotranspiration that would
occur if there was an infinite supply of water in the soil. The potential
amount can generally be computed based on atmospheric conditions.
Including evapotranspiration is usually not important for short events,
but becomes critical for longer simulations. There are seven different
evapotranspiration methods including an option for no
evapotranspiration. If you plan to use the meteorologic model with
basin models that contain subbasins, you may need to choose an
evapotranspiration method. The following loss methods require
evapotranspiration boundary conditions: deficit constant, gridded deficit
constant, soil moisture accounting, gridded soil moisture accounting. If
you use a meteorologic model with subbasins using any of those loss
methods, you should select an evapotranspiration method. If you use
those loss methods but no evapotranspiration method is selected in
the meteorologic model, then evapotranspiration is assumed to be zero
in the subbasins. A gridded evapotranspiration method is provided for
use with the ModClark gridded transform method. Using a gridded
evapotranspiration method provides separate boundary conditions for
each grid cell. Using a non-gridded evapotranspiration method means
that the same boundary condition will be applied at each grid cell in a
subbasin.
Select the evapotranspiration method in the Component Editor for the
meteorologic model (Figure 201). Only one evapotranspiration method
can be selected at a time. You can always change the
evapotranspiration method. When you change the evapotranspiration
method, all data for the old method is deleted and cannot be retrieved
or undeleted.
Some evapotranspiration methods require parameter data for each
subbasin. Other methods use the same data for all subbasins.
Depending on the method you choose, the Watershed Explorer will be
updated to indicate which data is required.
291
Chapter 11 Meteorology Description
Snowmelt Method
The snowmelt method determines if precipitation is rain or snow using
air temperature. The precipitation method is responsible for estimating
the total depth of combined rain and snow. Once the amount of
precipitation is determined by the precipitation method, the state of the
precipitation is determined by the snowmelt method. Subsequently,
the snowmelt method tracks the accumulation and melting of the
snowpack throughout the winter season. The snowmelt methods can
also operate during the summer season when precipitation is generally
purely rain and no snowpack accumulates. There are five different
snowmelt methods including an option for no snowmelt. Snowmelt is
only required if you plan to use the meteorologic model with basin
models that contain subbasins and the watershed receives
precipitation as snow, or as both rain and snow. If any of the
subbasins use the ModClark gridded transform method you must use a
gridded snowmelt method.
Select the snowmelt method in the Component Editor for the
meteorologic model (Figure 201). Only one snowmelt method can be
selected at a time. You can always change the snowmelt method.
When you change the snowmelt method, all data for the old method is
deleted and cannot be retrieved or undeleted.
Some snowmelt methods require parameter data for each subbasin,
and for elevation bands within each subbasin. Other methods use the
same data for all subbasins. Depending on the method you choose,
the Watershed Explorer will be updated to indicate which data is
required.
Unit System
Each meteorologic model must be in either United States customary
units (sometimes called English units) or in system international units
(also called metric units). All parameter data in a meteorologic model,
whether for all subbasins or individual subbasins, must be in the same
unit system. If you change the unit system, all data will be
automatically converted to the new unit system. All time-series data,
paired data, and gridded data referenced in a meteorologic model will
be in its own unit system. If necessary, these referenced data are
automatically converted to the unit system of the meteorologic model
during a simulation.
Select the unit system using the Component Editor for the
meteorologic model (Figure 201). Access the Component Editor by
clicking on the meteorologic model icon on the "Components" tab of
the Watershed Explorer. If you change the unit system, all data is
automatically converted to the new selection.
292
Chapter 11 Meteorology Description
Missing Data
Missing data is a common problem with atmospheric measurements.
The exact meaning of missing data often depends on the type of data.
Missing precipitation data often means that no precipitation was
occurring. Care must be used, but it is often possible to assume that
zero rainfall can be substituted for missing data. Other types of
atmospheric measurements are more difficult to manage when data is
missing. When air temperature is missing there is usually no
reasonable substitution that can be made.
There are two options for managing missing data in the meteorologic
model. The first option is to abort the compute with an error message.
Using this option means that when any missing data is encountered,
the simulation is stopped and an error message will report the time-
series gage or grid that contained the missing data. Using this
approach is the best because it requires all atmospheric data to be
reviewed and corrected before use in a simulation. The second option
is to assume a default value. For example, with this option selected
any missing precipitation data will be assumed to be zero. However,
allowing missing data to be replaced with a default value sometimes
leads to inconsistent results that can be difficult to intrepret.
Furthermore, some types of data such as temperature have no default
value. The compute will always abort with an error message when
missing data is encountered for data types with no defined default
value.
Select the action to take for missing data using the Component Editor
for the meteorologic model (Figure 201). Access the Component
Editor by clicking on the meteorologic model icon on the "Components"
tab of the Watershed Explorer.
293
Chapter 11 Meteorology Description
across all basin models that will be used with the same meteorologic
model.
294
Chapter 11 Meteorology Description
Atmospheric Variables
Atmospheric variables are measurable properties of the atmosphere
that change over time in response to external influences. The
measured values provide insight to short term weather phenomena as
well as long term climate trends. The measured values are also input
to a number of physically-based models used to quantify various
components of the hydrologic cycle including evapotranspiration and
snowmelt, among others. Meteorologists and climatologists have
established standard scientific methods for making measurements of
the atmospheric variables.
The most common atmospheric variable is air temperature which
measures the internal energy of the air. The average speed of moving
air is defined as the windspeed. The air pressure is defined as the
absolute pressure which is actually the weight of the atmosphere. The
amount of water vapor in the air can be measured either as relative
humidity or dew point temperature.
Atmospheric variables are required for each subbasin when certain
methods are selected for shortwave radiation, longwave radiation,
evapotranspiration, and snowmelt. The same variables may be
required for more than one method, for example, a shortwave radiation
method and evapotranspiration method may both require the same
variable. When necessary, time-series gages must be selected for
each subbasin in order to specify the atmospheric variables. Click on
the subbasin in the Watershed Explorer as shown in Figure 203 to
access the Component Editor for the atmospheric variables (Figure
204). Gages only need to be selected for atmospheric variables that
are required based on selections in the meteorologic model. Variables
are enabled when they are required. Variables are disabled when they
are not required and no selection can be made.
Observations for the atmospheric variables must be defined in the
Time-Series Data Manager before they can be selected in the
Component Editor. Chapter 5 of this manual describes how to create
and manage time-series data. The available gages are shown in the
selection list for each variable. If there are many different gages
available, you may wish to choose a gage from the selector accessed
with the gage button next to the selection list. The selector displays
the description for each gage, making it easier to select the correct
one.
295
Chapter 11 Meteorology Description
Shortwave Radiation
Shortwave radiation is a radiant energy produced by the sun with
wavelengths ranging from infared through visible to ultraviolet.
Shortwave radiation is therefore exclusively associated with daylight
hours for a particular location on the Earth's surface. The energy
arrives at the top of the Earth's atmosphere with a flux (Watts per
square meter) that varies very little during the year and between years.
Consequently the flux is usually taken as a constant for hydrologic
296
Chapter 11 Meteorology Description
Bristow Campbell
The Bristow Campbell method (Bristow and Campbell, 1984) uses a
conceptual approach to estimating the shortwave radiation at the land
surface. During the daylight hours, any clouds present in the
atmosphere will block some portion of incoming solar radiation which
reduces solar heating and results in a lower temperature. Conversely,
lack of clouds permits much more of the solar radiation to pass through
the atmosphere which allows greater heating and generally higher air
temperatures. In theory the daily temperature range should be small
on cloudy days and large on non-cloudy days. This correlation
between temperature range and incoming solar radiation is exploited
as a simple way to compute shortwave radiation using only air
temperature.
The Bristow Campbell method includes a Component Editor with
parameter data for each subbasin in the meteorologic model. The
Watershed Explorer provides access to the shortwave component
editor using a picture of solar radiation (Figure 205).
An air temperature gage must be selected in the atmospheric variables
for each subbasin.
297
Chapter 11 Meteorology Description
298
Chapter 11 Meteorology Description
Figure 206. Entering atmosphere and temperature data for a subbasin using
the Bristow Campbell shortwave radiation method.
FAO56
The FAO56 method implements the algorithm detailed by Allen,
Pereira, Raes, and Smith (1998). The algorithm calculates the solar
declination and solar angle for each time interval of the simulation,
using the coordinates of the subbasin, Julian day of the year, and time
at the middle of the interval. The solar values are used to compute the
extra-terrestrial radiation for each subbasin. Total daylight hours are
computed based on the Julian day and compared to the number of
actual sunshine hours. Shortwave radiation arriving at the ground
surface is then computed using the most common relationship
accounting for reduction in sunshine hours due to cloud cover.
The Watershed Explorer provides access to the shortwave component
editors using a picture of solar radiation (Figure 207). The FAO56
method includes a Component Editor with parameter data for all
subbasins in the meteorologic model (Figure 208); the same data is
also shared with the FAO56 longwave radiation method. A
Component Editor is also included for each subbasin (Figure 209).
An air temperature gage and a windspeed gage must be selected in
the atmospheric variables for each subbasin. Additionally, the water
299
Chapter 11 Meteorology Description
Figure 208. Entering the longitude of the central meridian of the local time
zone.
300
Chapter 11 Meteorology Description
Gridded Shortwave
The gridded shortwave method is designed to work with the ModClark
gridded transform. However, it can be used with other area-average
transform methods as well. The most common use of the method is to
utilize gridded shortwave radiation estimates produced by an external
model, for example, a dynamic atmospheric model. If it is used with a
transform method other than ModClark, an area-weighted average of
the grid cells in the subbasin is used to compute the shortwave
radiation time-series for each subbasin.
The gridded shortwave method includes a Component Editor with
parameter data for all subbasins in the meteorologic model. The
Watershed Explorer provides access to the shortwave component
editor using a picture of solar radiation (Figure 210).
301
Chapter 11 Meteorology Description
Figure 210. A meteorologic model using the gridded shortwave method with
a component editor for all subbasins in the meteorologic model.
Figure 211. Specifying the shortwave radiation data source for the gridded
shortwave method.
Specified Pyranograph
A pyranometer is an instrument that can measure incoming solar
shortwave radiation. They are not part of basic meteorological
observation stations, but may be included at first order stations. This
method may be used to import observed values from a pyranometer or
it may be used to import estimates produced by an external model.
This is the recommended choice for use with the Priestley Taylor
302
Chapter 11 Meteorology Description
Figure 213. Specifying the shortwave radiation time-series gage for each
subbasin.
303
Chapter 11 Meteorology Description
Longwave Radiation
All living and non-living bodies emit longwave radiation. The
magnitude of the radiation is proportional to the temperature
(measured in Kelvin degrees) of the body raised to the fourth power.
Significant sources of longwave radiation in hydrologic applications
include the atmosphere itself, and any clouds that may be present
locally in the atmosphere. Clouds usually have a higher heat content
and higher temperature than clear atmosphere, and therefore there is
increased downwelling longwave radiation on cloudy days. Whether
the atmosphere and clouds are a net source of longwave radiation to
the land surface depends on their temperature relative to the land
surface temperature. In most cases, the net longwave radiation is
incoming during the daylight hours, and outgoing during the night
hours.
The longwave radiation method included in the meteorologic model is
only necessary when energy balance methods are used for
evapotranspiration or snowmelt. The options available cover a range
of detail from simple to complex. Simple specified methods are also
available for input of a time-series or grid. Each option produces the
downwelling longwave radiation arriving at the land surface.
Absorption and radiation by the land surface is computed within the
subbasin element. More detail about each method is provided in the
following sections.
FAO56
The FAO56 method implements the algorithm detailed by Allen,
Pereira, Raes, and Smith (1998). The algorithm calculates the solar
declination and solar angle for each time interval of the simulation,
using the coordinates of the subbasin, Julian day of the year, and time
at the middle of the interval. The solar values are used to compute the
clear sky radiation for each subbasin. The computed shortwave
radiation is then divided by the calculated clear sky radiation to
estimate the fraction of shortwave radiation reaching the ground.
Finally, the fraction of shortwave radiation reaching the ground is
combined with the vapor pressure to calculate a reduction to the
downwelling longwave radiation based on the Stefan-Bolrzman
equation.
The FAO56 method includes a Component Editor with parameter data
for all subbasins in the meteorologic model; the same data is also
shared with the FAO56 shortwave radiation method. When FAO56 is
used for both shortwave and longwave radiation, the icon will
preferentially show the shortwave icon. A Component Editor is also
included for each subbasin. The Watershed Explorer provides access
304
Chapter 11 Meteorology Description
Figure 214. The FAO56 shortwave and FAO56 longwave methods use the
same editor for all subbasins.
305
Chapter 11 Meteorology Description
Figure 215. The FAO56 shortwave and FAO56 longwave methods both
require the time zone central meridian.
Gridded Longwave
The gridded longwave method is designed to work with the ModClark
gridded transform. However, it can be used with other area-average
transform methods as well. The most common use of the method is to
utilize gridded longwave radiation estimates produced by an external
model, for example, a dynamic atmospheric model. If it is used with a
transform method other than ModClark, an area-weighted average of
the grid cells in the subbasin is used to compute the longwave
radiation time-series for each subbasin.
The gridded longwave method includes a Component Editor with
parameter data for all subbasins in the meteorologic model. The
Watershed Explorer provides access to the longwave component
editor using a picture of downwelling radiation from clouds (Figure
217).
306
Chapter 11 Meteorology Description
Figure 217. A meteorologic model using the gridded longwave method with a
component editor for all subbasins in the meteorologic model.
Figure 218. Specifying the longwave radiation data source for the gridded
longwave method.
Satterlund
The Satterlund method (Satterlund, 1979) uses temperature and vapor
pressure to estimate incoming longwave radiation. While certain
physical limits are enforced, the approach is fundamentally empirical.
It provides good estimates under mostly clear skies and continues to
perform well even as temperature drops below the freezing point. It
307
Chapter 11 Meteorology Description
308
Chapter 11 Meteorology Description
Specified Pyrgeograph
A pyrgeometer is an instrument that can measure downwelling
longwave radiation. They are not part of basic meteorological
observation stations, but may be included at first order stations. This
method may be used to import observed values from a pyrgeometer or
it may be used to import estimates produced by an external model.
The specified pyrgeograph method includes a Component Editor with
parameter data for all subbasins in the meteorologic model. The
Watershed Explorer provides access to the longwave component
editors using a picture of heat radiation from a cloud (Figure 221).
The Component Editor for all subbasins in the meteorologic model
includes the time-series gage of longwave radiation for each subbasin
(Figure 222). A radiation gage must be selected for a subbasin. The
current gages are shown in the selection list.
309
Chapter 11 Meteorology Description
Figure 222. Specifying the longwave radiation time-series gage for each
subbasin.
Precipitation
Precipitation is water falling over the land surface. It is caused when
water vapor in the atmosphere condenses on airborne nuclei such as
dust particles. Condensation continues to add to the droplet forming
around the nuclei until the weight of the droplet exceeds the ability of
wind currents to keep it aloft. Precipitation can be caused by different
types of storms including stratiform, convection, and cyclone which
each have typical characteristics of duration, intensity, and spatial
extent. Precipitation includes the liquid form known as rain as well as
a variety of frozen forms including sleet, snow, graupel, and hail. Most
hydrologic purposes can be met by limiting consideration to rain and
snow. The determination of the rain or snow state is made separately
in the snowmelt portion of the meteorologic model.
The precipitation method included in the meteorologic model is
required whenever a basin model includes subbasin elements. The
options available include several that process gage measurements, a
statistical method that uses depth-duration data, several design storms
options, and a gridded method that can be used with radar rainfall
data. Each option produces a hyetograph of precipitation falling over
each subbasin. More detail about each method is provided in the
following sections.
Frequency Storm
The frequency storm method is designed to produce a synthetic storm
from statistical precipitation data. The most common source of
statistical data in the United States is the National Weather Service.
Typically the data is given in the form of maps, where each map shows
the expected precipitation depth for a storm of specific duration and
exceedance probability. This method is designed to use data collected
310
Chapter 11 Meteorology Description
311
Chapter 11 Meteorology Description
Figure 224. Precipitation editor for all subbasins when the frequency storm
method is selected.
The intensity position determines where in the storm the period of peak
intensity will occur. Changing the position does not change the total
depth of the storm, but does change how the total depth is distributed
in time during the storm. You may select 25%, 33%, 50%, 67%, or
75% from the list of choices. If the storm duration is selected to be 6
hours and the 25% position is selected, the peak intensity will occur
1.5 hours after the beginning of the storm. The default selection is
50%.
The storm area is used to automatically compute the depth-area
reduction factor. In most cases the specified storm area should be
equal to the watershed drainage area at the point of evaluation. The
same hyetograph is used for all subbasins. Optionally you may leave
the storm area blank. When no storm area is specified, each subbasin
will have a different hyetograph computed using the subbasin area as
the storm area.
Precipitation depth values must be entered for all durations from the
peak intensity to the total storm length. Values for durations less than
the peak intensity duration, or greater than the total storm duration are
312
Chapter 11 Meteorology Description
Gage Weights
The gage weights method is designed to work with recording and non-
recording precipitation gages. Recording gages typically measure
precipitation as it occurs and then the raw data are converted to a
regular time step, such as 1 hour. Non-recording gages usually only
provide an estimate of the total storm depth. The user can choose any
method to develop the weights applied to each gage when calculating
the hyetograph for each subbasin. For increased flexibility, the total
storm depth and the temporal pattern are developed separately for
each subbasin. Optionally, an index depth can be assigned to each
gage and subbasin. The index is used to adjust for regional bias in
annual or monthly precipitation.
313
Chapter 11 Meteorology Description
Figure 226. Selecting options for the gage weights precipitation method.
314
Chapter 11 Meteorology Description
Figure 228. Specifying indices for each subbasin using average annual
precipitation depth. Other methods may be used to estimate the
index.
315
Chapter 11 Meteorology Description
Figure 229. Entering optional total depth and index for the precipitation
gages.
The optional index is entered for each gage. If you enter an index for a
gage, it can only be used during a simulation if you also specify an
index for all gages used in a subbasin and also specify an index for the
subbasin. Turning off the "Use Indexing" option will remove all indexes
from the meteorologic model.
Total storm gages are created and managed directly from the
meteorologic manager. To create a new total storm gage, access the
"Total Storm Gages" tab and enter a gage name in the first column of
the last row (Figure 230). The last row is always kept blank for
creating new gages. You can rename a gage by typing over the name
in the first column. You can delete a gage by deleting its name from
the first column and any data from the other columns. You must
always enter a total depth for the gage that represents the total
precipitation during the simulation.
316
Chapter 11 Meteorology Description
Figure 230. Creating a total storm gage and entering the total storm depth
and optional index.
317
Chapter 11 Meteorology Description
Figure 232. Entering depth and time weights for the selected gages.
Gridded Precipitation
The gridded precipitation method is designed to work with the
ModClark gridded transform. However, it can be used with other area-
average transform methods as well. The most common use of the
method is to utilize radar-based precipitation estimates. Using
additional software, it is possible to develop a gridded representation of
gage data or to use output from atmospheric models. If it is used with
a transform method other than ModClark, an area-weighted average of
the grid cells in the subbasin is used to compute the precipitation
hyetograph for each subbasin.
The gridded precipitation method includes a Component Editor with
parameter data for all subbasins in the meteorologic model. The
Watershed Explorer provides access to the precipitation component
editor using a picture of raindrops (Figure 233).
318
Chapter 11 Meteorology Description
The time shift can be used to correct for precipitation grids stored with
a time zone offset. All calculations during a simulation are computed
assuming an arbitrary local time zone that does not observe summer
time (daylight savings in the United States). It is common for
precipitation data from radar sources to be referenced in Coordinated
Universal Time (UTC). Set the shift to zero if all the time-series and
grid data is referenced in same local time zone. If other data sources
such as observed discharge or temperature are referenced in local
time and the precipitation grid is in UTC, select the correct shift so that
the precipitation data will match the rest of the data. Local time zones
located to the West of the zero longitude line will use a positive shift
when the precipitation grid is referenced in UTC. Local time zones
located to the East of the zero longitude line will use a negative shift
when the precipitation grid is referenced in UTC.
HMR 52 Storm
The HMR 52 storm is one approach to computing the probable
maximum precipitation for a watershed as detailed in
Hydrometeorological Report No 52 (Hansen, Schreiner, and Miller,
1982). Concentric ellipses are used to construct the storm spatial
pattern where each ellipse represents an isohyet of precipitation depth.
The storm is located over the watershed by specifying the center of the
pattern and the angle of the major axis of the ellipses. Total
precipitation depth is computed using a specified storm area and area-
319
Chapter 11 Meteorology Description
320
Chapter 11 Meteorology Description
The HMR 52 storm is 72 hours long and begins at the simulation start
time. The peak intensity parameter specifies the period within the 72-
hour storm when the precipitation rate will be greatest. The 6-hour
period of peak intensity can be set to begin as early as hour 24 of the
storm or as late as hour 60 of the storm. The depth of rain falling
during the period of peak intensity is subdivided into 1-hour increments
using the parameter for the ratio of the 1-hour to 6-hour depth.
The total storm area must be specified. Additionally, a duration-
precipitation function must be selected for the range of possible storm
areas. In general, the curves are constructed using data available in
Hydrometeorological Report No 51 (Schreiner and Riedel, 1978).
Each duration-precipitation function must be defined in the paired data
manager before it can be selected. You can press the paired data
button next to the selection list to use a chooser. The chooser shows
all of the available duration-precipitation functions in the project.
The HMR 52 storm method creates a storm from input specified by the
user including storm center, orientation, and area. A new analysis tool
is under development which would automatically select these three
parameters in order to maximize precipitation volume over the
watershed, or peak runoff, or volume, or reservoir pool elevation.
Inverse Distance
The inverse distance method was originally designed for application in
real-time forecasting systems. It can use recording gages that report
on a regular interval like 15 minutes or 1 hour. It can also use gages
that only report daily precipitation totals. Because it was designed for
real-time forecasting, it has the ability to automatically switch from
using close gages to using more distant gages when the closer gages
stop reporting data. The latitude and longitude of the gages is used to
determine closeness to one or more nodes specified in each subbasin.
The distance at which gages are used for a subbasin is controlled by
the search distance. Optionally, an index depth can be assigned to
each gage and subbasin node. The index is used to adjust for regional
bias in annual or monthly precipitation.
The inverse distance method includes several options for processing
precipitation gage data (Figure 236). The indexing option is used to
adjust gage data when there are regional trends in precipitation
patterns. Using the index option requires the specification of an index
value for each precipitation gage and each node in each subbasin.
The average annual precipitation total is often used as the index at a
gage, and the estimated average annual precipitation at a node
location is often used as the index for a subbasin node. Alternately,
the monthly average values may be used at each gage and each node.
The search distance option can be used to limit the influence distance
321
Chapter 11 Meteorology Description
Figure 236. Selecting options for the inverse distance precipitation method.
322
Chapter 11 Meteorology Description
Each precipitation gage may use any of the allowable time intervals
from 1 minute to 24 hours. When a gage is selected as a daily gage, it
is assumed during calculations that the daily precipitation depth is
known but there is no timing information. Usually daily gages are
stored with a 24-hour time interval, but any interval may be used. Daily
gage data is only used during processing for days where the entire day
is within the simulation time window. When a gage is not selected as a
daily gage, the data it contains is interpolated to the simulation time
step. The appropriate setting should be made in the "Daily Gage"
column.
The optional index is entered for each gage. If you enter an index for a
gage, it can only be used during a simulation if you also specify an
index for all gages used in a subbasin and also specify an index for
each node in that subbasin. Turning off the "Use Indexing" option will
remove all indexes from the meteorologic model.
The Component Editor for each subbasin in the meteorologic model is
used to create one or more nodes and enter the coordinates for each
node (Figure 239). A minimum of one node is required. To create a
new node, access the "Node Weights" tab and enter a node name in
the first column of the last row. The last row is always kept blank for
creating new nodes. You can rename a node by typing over the name
in the first column. You can delete a node by deleting its name from
the first column and any data from the other columns. You must
always enter a weight for a node. The weight controls how the final
hyetograph is computed for the subbasin from the hyetographs
computed at each node. If you have enabled the indexing option, you
should also enter the index for each node. The index is used to adjust
for regional bias in annual or monthly precipitation.
323
Chapter 11 Meteorology Description
All nodes that you have created will be shown on the "Latitudes" and
"Longitudes" tabs. You must enter the appropriate coordinate
information for each node. Coordinates may be entered using degree-
minute-second format or alternately may be entered using decimal
degree format. The data entry format is selected in the Program
Settings.
SCS Storm
The Soil Conservation Service (SCS) storm method implements the
design storm patterns developed developed for use in rural,
agricultural watersheds (NRCS, 2007). The original methodology was
developed to provide guidance when designing safe water storage
facilities for agricultural applications. However, the storm has been
applied in a variety of other situations. The same hyetograph is
computed for all subbasins.
The SCS storm method includes a Component Editor with parameter
data for all subbasins in the meteorologic model. The Watershed
Explorer provides access to the precipitation component editor using a
picture of raindrops (Figure 240).
324
Chapter 11 Meteorology Description
Figure 240. A meteorologic model using the SCS storm precipitation method
with a component editor for all subbasins.
Figure 241. Precipitation component editor for the SCS storm precipitation
method.
Specified Hyetograph
The specified hyetograph method allows the user to specify the exact
time-series to use for the hyetograph at subbasins. This method is
useful when precipitation data will be processed externally to the
program and essentially imported without alteration. This method is
325
Chapter 11 Meteorology Description
326
Chapter 11 Meteorology Description
Figure 244. Selecting a gage for each subbasin. Total depth override is
optional.
Optionally, you may enter a total depth for each subbasin. If no total
depth is entered, the depth will be the sum of the data actually stored
in the precipitation gage. However, if a total depth is entered for a
subbasin, the exact pattern is maintained but the magnitude of
precipitation at each time step is adjusted proportionally so that the
specified depth is applied over the entire simulation. Total depth can
be specified for no subbasins, one subbasin, many subbasins, or all
subbasins. It is not required to enter the depth for all subbasins in
order to specify it for just one subbasin.
327
Chapter 11 Meteorology Description
Figure 245. A meteorologic model using the standard project storm method
with a component editor for all subbasins.
328
Chapter 11 Meteorology Description
Evapotranspiration
Evapotranspiration is the the combination of evaporation from the
ground surface and transpiration by vegetation. It includes both
evaporation of free water from the surface of vegetation and the land
surface. It also includes transpiration which is the process of
vegetation extracting it from the soil through the plant root system.
Whether by evaporation or transpiration, water is returned from the
land surface or subsurface to the atmosphere. Even though
evaporation and transpiration are taken together, transpiration is
responsible for the movement of much more water than evaporation.
Combined evapotranspiration is often responsible for returning 50 or
even 60% of precipitation back to the atmosphere. The theoretical
evapotranspiration, also called the potential evapotranspiration, serves
as the upper limit for what can happen on the land surface based on
atmospheric conditions. In all cases, the meteorologic model is
computing the potential evapotranspiration and subbasins will calculate
actual evapotranspiration based on soil water limitations.
The evapotranspiration method included in the meteorologic model is
only necessary when using continuous simulation loss rate methods in
subbasins: deficit constant, gridded deficit constant, soil moisture
329
Chapter 11 Meteorology Description
Annual Evapotranspiration
The annual evapotranspiration method is designed to work with a
maximum daily rate combined with an optional pattern of variation
throughout the year. Specifying only a daily rate can produce good
results for simulations lasting days to weeks if evapotranspiration is
fairly consistent each day. The optional pattern can be used to adjust
the applied evapotranspiration rate during simulations lasting weeks to
years.
The annual evapotranspiration method includes a Component Editor
with parameter data for all subbasins in the meteorologic model. The
Watershed Explorer provides access to the evapotranspiration
component editor using a picture of a water pan.
330
Chapter 11 Meteorology Description
the optional percentage pattern is not used, the value entered for the
daily rate should generally be the average daily potential
evapotranspiration rate over the duration of the simulation. When a
pattern will be added, the value entered for the daily rate should
generally be the largest potential evapotranspiration for any day
occurring during the simulation.
Figure 249. Entering rates and selecting a percent pattern for each subbasin.
331
Chapter 11 Meteorology Description
332
Chapter 11 Meteorology Description
333
Chapter 11 Meteorology Description
displays the description for each gridset, making it easier to select the
correct one.
Monthly Average
The monthly average method is designed to work with measured pan
evaporation data. However, it can also be used with data collected
with the eddy correlation technique or other modern methods.
Regardless of how they are collected, the data are typically presented
as the average depth of evaporated water each month. Maps or
tabular reports can be found for each month and used with this
method.
The monthly average method includes a Component Editor with
parameter data for each individual subbasin in the meteorologic model.
The Watershed Explorer provides access to the evapotranspiration
component editor using a picture of a water pan (Figure 254).
334
Chapter 11 Meteorology Description
Figure 255. Entering rate and pan coefficient data for a subbasin in a
meteorologic model using the monthly average
evapotranspiration method.
Penman Monteith
The Penman Monteith method implements the Penman Monteith
equations for computing evapotranspiration at less than a daily time
interval as detailed by Allen, Pereira, Raes, and Smith (1998). The
equations are based on a combination of an energy balance with a
mass transfer. The maximum possible evapotranspiration is
moderated by an aerodynamic resistance due to friction as air flows
over the vegetation. A bulk surface resistance is added in series with
the aerodynamic resistance to account for limitations to water vapor
flow at the leaf surfaces and at the soil. The parameterization is
entirely dependent on the atmospheric conditions.
The Penman Monteith method requires shortwave radiation and
longwave radiation be included in the meteorologic model. The
algorithm of Allen, Pereira, Raes, and Smith (1998) is followed most
closely when the FAO56 shortwave method and the FAO56 longwave
method are selected.
The Penman Monteith method includes a Component Editor with
parameter data for each individual subbasin in the meteorologic model.
335
Chapter 11 Meteorology Description
Priestley Taylor
The Priestley Taylor method (Priestley and Taylor, 1972) uses a
simplified energy balance approach where the soil water supply is
assumed to be unlimited. Simplified forms of latent and sensible
336
Chapter 11 Meteorology Description
337
Chapter 11 Meteorology Description
Specified Evapotranspiration
The specified evapotranspiration method allows the user to specify the
exact time-series to use for the potential evapotranspiration at
subbasins. This method is useful when atmospheric and vegetation
data will be processed externally to the program and essentially
imported without alteration. This method is also useful when a single
evapotranspiration observation measurement can be used to represent
what happens over a subbasin.
The specified evapotranspiration method uses a Component Editor
with parameter data for all subbasins in the meteorological model. The
Watershed Explorer provides access to the gage component editor
using a picture of a water pan (Figure 260).
338
Chapter 11 Meteorology Description
model. The data may actually be from daily pan measurements, hourly
eddy covariance measurements, or could be the result of complex
calculations exterior to the program. Regardless, the time-series must
be stored as a gage. You may use the same gage for more than one
subbasin. For each subbasin in the table, select the gage to use for
that subbasin. Only evapotranspiration gages already defined in the
Time-Series Data Manager will be shown in the selection list.
Snowmelt
Precipitation generally occurs as snow when the air temperature over
the land surface is below the freezing point. In general the snow will
accumulate on the land surface so long as the air temperature remains
below the freezing point. In some watersheds it is part of the annual
hydrologic cycle that snowfall accumulates in a snowpack during the
winter months. Snowpacks are not static but undergo continual ice
crystal transformations in response to diurnal temperature fluctuations.
Ablation or melting of the snowpack occurs when atmospheric
conditions transfer energy sufficient to raise the temperature of the
snowpack to the freezing point, and then transform the snow to liquid
water through the heat of fusion for ice. The most common way to
measure the water content of the snowpack is by the snow water
equivalent or SWE. The SWE is the depth of water that results from
melting a unit column of the snowpack.
The snowmelt method included in the meteorologic model is only
necessary when the air temperature may be below the freezing point
during a simulation, or a snowpack may exist at the beginning of a
simulation. Currently the only option available is a temperature index
approach that includes a conceptual representation of the snowpack
energy. It computes the liquid water available at the soil surface which
339
Chapter 11 Meteorology Description
340
Chapter 11 Meteorology Description
zone that does not observe summer time (daylight savings in the
United States). In some cases, the temperature data may be stored in
universal time, also known as Greenwich Mean Time (GMT). If the
simulation will happen in a local time zone because of other data
sources such as observed discharge, select the correct shift so that the
temperature data will match the rest of the data. For local time zones
West of the zero longitude line, the shift will be negative. For local time
zones East of the zero longitude line, the shift will be positive.
Figure 263. Component editor for the gridded temperature index snowmelt
method.
341
Chapter 11 Meteorology Description
342
Chapter 11 Meteorology Description
The difference between the base temperature and the air temperature
defines the temperature index used in calculating snowmelt. The
meltrate is multiplied by the difference between the air temperature
and the base temperature to estimate the snowmelt amount. If the air
temperature is less than the base temperature, then the amount of
melt is zero. Typically, the base temperature should be 0C.
The wet meltrate is used during time intervals when precipitation is
falling as rain, and the rainfall rate is greater than the rain rate limit. It
represents the rate at which the snowpack melts when it is raining on
the snowpack.
The rain rate limit discriminates between dry melt and wet melt. The
wet meltrate is applied as the meltrate when it is raining at rates
greater than the rain rate limit. If the rain rate is less than the rain rate
limit, the meltrate is computed as if there were no precipitation. The
default value of 0mm/day is used if no value is entered, meaning that
even a trace of precipitation results in the use of the wet melt rate.
A meltrate must be calculated for time intervals when the precipitation
rate is less than the rain rate limit. The calculation starts with the
meltrate antecedent temperature index. A coefficient is used to update
the antecedent meltrate index from one time interval to the next. The
default value of 0.98 is used if no value is entered.
An antecedent temperature index meltrate function is used to calculate
a meltrate from the current meltrate index. The function must be
specified separately in the Paired Data Manager before is can be used
in the snow melt method. The function should define appropriate
meltrates to use over the range of meltrate index values that will be
encountered during a simulation. The melt rate for dry conditions
typically ranges from 1 to 4 mm/degC-day.
Optionally, you may adjust the meltrate computed from the index
meltrate function. A meltrate pattern may be specified that defines the
percentage adjustment as a function of the time of year. If no meltrate
pattern is selected, the meltrate will be computed only from the
antecedent temperature index and the meltrate function. The pattern
must be specified separately in the Paired Data Manager before is can
be used in the snow melt method.
The cold limit accounts for the rapid changes in temperature that the
snowpack undergoes during high precipitation rates. When the
precipitation rate exceeds the specified cold limit, the antecedent cold
content index is set to the temperature of the precipitation. If the
temperature is above the base temperature, the cold content index is
set to the base temperature. If the temperature is below the base
temperature, the cold content index is set to the actual temperature. If
343
Chapter 11 Meteorology Description
the precipitation rate is less than the cold limit, the cold content index is
computed as an antecedent index. The default value of 0mm/day is
used if no value is entered, meaning that even a trace of snowfall will
reset the cold content.
The cold content antecedent temperature index coefficient is used to
update the antecedent cold content index from one time interval to the
next. This is a separate index from the one used to update the
meltrate index. A default value of 0.5 is used if no value is entered.
An antecedent temperature index cold content function is used to
calculate a cold content from the current cold content index. The
function must be specified separately in the Paired Data Manager
before is can be used in the snow melt method. The function should
define appropriate cold contents to use over the range of cold content
index values that will be encountered during a simulation. The cold
rate typically ranges from 1.22 to 1.32 mm/degC-day.
The maximum liquid water capacity specifies the amount of melted
water that must accumulate in the snowpack before liquid water
becomes available at the soil surface for infiltration or runoff. Typically,
the maximum liquid water held in the snowpack is on the order of 3%-
5% of the snow water equivalent, although it can be higher. Liquid
water can persist in the snow only if the snowpack temperature is at
0C; at which point the cold content is zero. The maximum is entered
as a percentage of the snow water equivalent.
Heat from the ground can cause snowmelt, especially if the snowpack
accumulates on ground that is only partially frozen or completely
unfrozen. In these cases the warm ground is insulated by the
snowpack. Heat from the warm ground will cause the bottom of the
snowpack to melt. Two methods are available for specifying the
melting of the snowpack due to contact with unfrozen ground. A fixed
value can be entered; the same amount of melt is computed for the
snowpack regardless of atmospheric conditions above the pack or the
time of year. An annual pattern can alternately be entered; the pattern
specifies the meltrate due to contact with the ground as a function of
the time of year. The pattern must be entered in the Paired Data
Manager before is can be used in the snowmelt method.
Temperature Index
The temperature index method is an extension of the degree-day
approach to modeling a snowpack. A typical approach to the degree
day is to have a fixed amount of snowmelt for each degree above
freezing. This method includes a conceptual representation of the cold
energy stored in the pack along with a limited memory of past
conditions and other factors to compute the amount of melt for each
344
Chapter 11 Meteorology Description
345
Chapter 11 Meteorology Description
The difference between the base temperature and the air temperature
defines the temperature index used in calculating snowmelt. The
meltrate is multiplied by the difference between the air temperature
and the base temperature to estimate the snowmelt amount. If the air
temperature is less than the base temperature, then the amount of
melt is zero. Typically, the base temperature should be 0C.
346
Chapter 11 Meteorology Description
the antecedent meltrate index from one time interval to the next. The
default value of 0.98 is used if no value is entered.
An antecedent temperature index meltrate function is used to calculate
a meltrate from the current meltrate index. The function must be
specified separately in the Paired Data Manager before is can be used
in the snow melt method. The function should define appropriate
meltrates to use over the range of meltrate index values that will be
encountered during a simulation. The melt rate for dry conditions
typically ranges from 1 to 4 mm/degC-day.
Optionally, you may adjust the meltrate computed from the index
meltrate function. A meltrate pattern may be specified that defines the
percentage adjustment as a function of the time of year. If no meltrate
pattern is selected, the meltrate will be computed only from the
antecedent temperature index and the meltrate function. The pattern
must be specified separately in the Paired Data Manager before is can
be used in the snow melt method.
The cold limit accounts for the rapid changes in temperature that the
snowpack undergoes during high precipitation rates. When the
precipitation rate exceeds the specified cold limit, the antecedent cold
content index is set to the temperature of the precipitation. If the
temperature is above the base temperature, the cold content index is
set to the base temperature. If the temperature is below the base
temperature, the cold content index is set to the actual temperature. If
the precipitation rate is less than the cold limit, the cold content index is
computed as an antecedent index. The default value of 0mm/day is
used if no value is entered, meaning that even a trace of snowfall will
reset the cold content.
The cold content antecedent temperature index coefficient is used to
update the antecedent cold content index from one time interval to the
next. This is a separate index from the one used to update the
meltrate index. A default value of 0.5 is used if no value is entered.
An antecedent temperature index cold content function is used to
calculate a cold content from the current cold content index. The
function must be specified separately in the Paired Data Manager
before is can be used in the snow melt method. The function should
define appropriate cold contents to use over the range of cold content
index values that will be encountered during a simulation. The cold
rate typically ranges from 1.22 to 1.32 mm/degC-day.
The maximum liquid water capacity specifies the amount of melted
water that must accumulate in the snowpack before liquid water
becomes available at the soil surface for infiltration or runoff. Typically,
the maximum liquid water held in the snowpack is on the order of 3%-
347
Chapter 11 Meteorology Description
348
Chapter 11 Meteorology Description
Each subbasin is broken into one or more elevation bands; each band
has its own parameter data. One elevation band may be used to
represent a subbasin with very little terrain variation. Subbasins with
large elevation variations should use multiple elevation bands. Create
an elevation band by clicking with the right mouse button on the
snowmelt icon under a subbasin (Figure 267). A context menu is
displayed that allows you to create a new elevation band. You can
also create an elevation band by clicking with the right mouse button
on any existing elevation band. The same context menu is displayed
that allows you to create a new elevation band. Delete an elevation
band by clicking on it with the right mouse button. A context menu is
displayed that allows you to delete the elevation band. Clicking on the
elevation band icon will access the Component Editor used for each
band (Figure 268).
349
Chapter 11 Meteorology Description
350
Chapter 11 Meteorology Description
References
Allen, R.G., L.S. Pereira, D. Raes, and M. Smith. 1998. "Crop
Evapotranspiration: Guidelines for Computing Crop Water
Requirements." FAO Irrigation and Drainage Paper No. 56, Rome,
Italy.
Bristow, K.L. and G.S. Campbell. 1984. "On the Relationship
Between Incoming Solar Radiation and Daily Maximum and Minimum
Temperature." Agricultural and Forest Meteorology, vol 31, pp 159-
166.
Engineering Manual 1110-2-1411. March 1065. Standard Project
Flood Determination. U.S. Army Corps of Engineers, Davis, CA.
FAO Irrigation and Drainage Paper 56. 1998. Crop
Evapotranspiration: Guidelines for Computing Crop Water
Requirements. Food and Agriculture Organization, Rome.
351
Chapter 11 Meteorology Description
352
Chapter 12 Hydrologic Simulation
CHAPTER 12
Hydrologic Simulation
This chapter describes how watershed and meteorology information is
combined to simulate the hydrologic response. The simulation run is
the primary mode for performing simulations. In effect the simulation
run applies the meteorology conditions to the watershed land surface
to determine the runoff. In addition to configuring and computing
simulation runs, this chapter also describes how to access the
computed results.
Simulation Runs
Simulation runs are the primary method of computing results. Each
run is composed of one meteorologic model, one basin model, and one
control specifications. Results can be visualized as graphs, summary
tables, and time-series tables either through the basin map or from the
Watershed Explorer. In addition to selecting the meteorologic, basin,
and control components, advanced features for controlling the run are
also included.
353
Chapter 12 Hydrologic Simulation
button when you are satisfied with the name you have entered and are
ready to proceed to the next step.
The second step of creating a simulation run is to select a basin model.
All of the basin models in the project are shown and you must select
one before proceeding to the next step. By default the first basin
model in the table is selected. The selected model is highlighted. You
can use your mouse to select a different basin model by clicking on it
in the table of available choices. You can also use the arrow keys on
your keyboard to select a different model. Press the Next> button
when you are satisfied with the basin model you have selected and are
ready to proceed to the next step. Press the <Back button if you wish
to return to the previous step and change the name for the new
simulation run.
Figure 269. Beginning the process of creating a new simulation run using the
Simulation Run Manager.
354
Chapter 12 Hydrologic Simulation
Figure 270. Entering a name for a new simulation run. The remaining steps
are to select a basin model, meteorologic model, and control
specifications.
Copying a Run
There are two ways to copy a simulation run. Both methods for
copying a run create an exact duplicate with a different name. Once
the copy has been made it is independent of the original and they do
not interact.
The first way to create a copy is to use the Simulation Run Manager,
which is accessed from the Compute menu. Select the simulation run
you wish to copy by clicking on it in the list of available simulation runs.
The selected run is highlighted after you select it. After you select a
run you can press the Copy button on the right side of the window.
A new Copy Simulation Run window (Figure 271) will open where you
can name and describe the copy that will be created. A default name
355
Chapter 12 Hydrologic Simulation
is provided for the copy; you can use the default or replace it with your
own choice. A description can also be entered; if it is long you can use
the button to the right of the description field to open an editor. When
you are satisfied with the name and description, press the Copy button
to finish the process of copying the selected simulation run. You
cannot press the Copy button if no name is specified. If you change
your mind and do not want to copy the selected simulation run, press
the Cancel button or the X button in the upper right to return to the
Simulation Run Manager window.
The second way to copy is from the "Compute" tab of the Watershed
Explorer. Move the mouse over the simulation run you wish to copy,
then press the right mouse button (Figure 272). A context menu is
displayed that contains several choices including copy. Click the
Create Copy command. A new Copy Simulation Run window will
open where you can name and describe the copy that will be created.
A default name is provided for the copy; you can use the default or
replace it with your own choice. A description can also be entered; if it
is long you can use the button to the right of the description field to
open an editor. When you are satisfied with the name and description,
press the Copy button to finish the process of copying the selected
simulation run. You cannot press the Copy button if no name is
specified. If you change your mind and do not want to copy the
selected simulation run, press the Cancel button or the X button in the
upper right of the Copy Simulation Run window to return to the
Watershed Explorer.
356
Chapter 12 Hydrologic Simulation
Figure 272. Copying a simulation run from the Watershed Explorer. The
Copy Simulation Run window will appear after the Create Copy
menu command is selected.
Renaming a Run
There are two ways to rename a simulation run. Both methods for
renaming a run change its name and then all references to the old run
name are automatically updated to the new name.
The first way to perform a rename is to use the Simulation Run
Manager, which you can access from the Compute menu. Select the
simulation run you wish to rename by clicking on it in the list of
available simulation runs. The selected run is highlighted after you
select it. After you select a run you can press the Rename button on
the right side of the window. A new Rename Simulation Run window
(Figure 273) will open where you can provide the new name. If you
wish you can also change the description at the same time. If the new
description will be long, you can use the button to the right of the
description field to open an editor. When you are satisfied with the
name and description, press the Rename button to finish the process
of renaming the selected simulation run. You cannot press the
Rename button if no name is specified. If you change your mind and
do not want to rename the selected simulation run, press the Cancel
button or the X button in the upper right of the Rename Simulation Run
window to return to the Simulation Run Manager window.
357
Chapter 12 Hydrologic Simulation
Figure 273. Renaming a simulation run. This was accessed from the
Simulation Run Manager.
Deleting a Run
There are two ways to delete a simulation run. Both methods for
deleting a run remove it from the project, automatically update all
references to that run, and delete previously computed results. Once a
run has been deleted it cannot be retrieved or undeleted. Any
references to the deleted run will switch to using no simulation run,
which is usually not a valid choice during a simulation. At a later time
358
Chapter 12 Hydrologic Simulation
Figure 275. Preparing to delete a simulation run from the Simulation Run
Manager. A confirmation will be required after pressing the
Delete button.
The second way to delete is from the "Compute" tab of the Watershed
Explorer. Move the mouse over the simulation run you wish to delete
and press the right mouse button (Figure 276). A context menu is
displayed that contains several choices including delete. Click the
Delete command. A window will open where you must confirm that
you wish to delete the selected run. Press the OK button to delete the
run. If you change your mind and do not want to delete the selected
simulation run, press the Cancel button or the X button in the upper
right to return to the Watershed Explorer.
359
Chapter 12 Hydrologic Simulation
Importing a Run
A simulation run can be imported into the current project. The run to
be imported is selected from a source project. Once you select the run
to be imported, all of the necessary simulation components will be
imported into the current project. This includes the basin model,
meteorologic model, control specifications, and any supporting data.
The supporting data include background maps and grid cell file for the
basin model, and shared data such as time-series gages, paired data,
and grid data used either in the basin or meteorologic model. A
simulation run will be created in the current project using the
components that have been imported.
To import a simulation run, begin by selecting the
FileImportSimulation Run menu command. Use the file chooser
to select the project file for the project that contains the simulation run
you wish to import (Figure 277). You will then be provided a selection
list of available simulation runs in that source project (Figure 278).
Select the simulation run that you wish to import and press the Import
button.
360
Chapter 12 Hydrologic Simulation
Figure 278. Selecting the simulation run to be imported from the source
project into the current project.
361
Chapter 12 Hydrologic Simulation
Selecting Components
The principal task when creating a simulation run using the wizard is
the selection of a basin model, meteorologic model, and control
specifications. However, you can change the components you wish to
use at any time using the Component Editor for a simulation run.
Access the Component Editor from the "Compute" tab of the
Watershed Explorer. If necessary, click on the "Simulation Runs"
folder to expand it and view the available simulation runs. Click on the
simulation run node you wish to edit and its Component Editor will
automatically be shown (Figure 279). The Component Editor contains
three lists for selecting the basin model, meteorologic model, and
control specifications, respectively. The first selection list shows all of
the basin models currently available in the project. Likewise, the
second list contains all of the meteorologic models and the third list
contains all of the control specifications. As the user, you are
responsible for knowing which components are intended to be used
together to produce simulation results.
362
Chapter 12 Hydrologic Simulation
Figure 279. The simulation run component editor can be used to change a
selected component after the run has been created.
363
Chapter 12 Hydrologic Simulation
Figure 280. Entering a precipitation ratio for all of the subbasins in the
simulation run. If a flow ratio were selected, it could be applied
to subbasins or sources or both subbasins and sources.
364
Chapter 12 Hydrologic Simulation
one long simulation. Suppose that a long simulation covers all of 1995
and 1996. You could configure a simulation run for 1995 and save
states at the end of the simulation. You could then configure a
simulation for 1996 that uses as start states the states saved at the
end of the 1995 run. The results obtained by breaking the simulation
into two runs and using states would be identical to the results
obtained from a single long run. This approach can be useful when
performing continuous simulation for many years or decades. In the
second case of real-time forecasting operations, the typical approach
is often to compute a simulation run daily and forecast three to five
days into the future. The simulation can be configured to run for five
days, saving the state variables at the end of the first day. When the
forecast is updated on the subsequent day, it can start from the saved
states of the previous day and continue forward. In this way,
continuously updated results can be obtained without requiring the
simulation to start at the last known watershed conditions which might
be very old.
Saving out the start variable during a simulation run, or starting a run
from saved state variables is optional. If you do not use start states,
then the initial conditions will be specified from the basin and
meteorologic models. If you do use the start states, the initial
conditions specified in the basin and meteorologic models will be
overridden by the values contained in the saved states used as the
start states. If you do not use save states, the simulation will run
normally. If you do choose to save states, the simulation will pause
momentarily to save the state variables and then continue normally.
When start states are used, they must have been saved at exactly the
same time as the start time in the selected control specifications.
Selecting optional start states for the simulation run is accessed from
the "Compute" tab of the Watershed Explorer. Click on the "Start
States" node under the simulation run to display the Component Editor
for the start states (Figure 281). The "Start States" and other optional
365
Chapter 12 Hydrologic Simulation
Figure 282. Configuring a simulation run to save states at the end of the
compute time window.
366
Chapter 12 Hydrologic Simulation
367
Chapter 12 Hydrologic Simulation
368
Chapter 12 Hydrologic Simulation
Figure 283. Computing simulation run results to an element that has been
designated a computation point.
369
Chapter 12 Hydrologic Simulation
370
Chapter 12 Hydrologic Simulation
Adjust a parameter value for an element using the slider bar for that
parameter. Each slider bar is labeled with the element name and
parameter name, including units for the parameter. Change the value
by clicking with the mouse and dragging the slider. The precise
selected value is shown on the right of each slider. You may optionally
change the value by typing in the editor field where the precise value is
displayed.
The results for the computation point and all elements upstream of it
are recomputed immediately after changing a parameter value with a
slider bar. In many watersheds the computations require only a few
seconds to complete. However, some complex watersheds require
more time to recompute the results. In these cases you can use the
checkbox in the upper right to control when results are recomputed.
When the checkbox is engaged, results are recomputed after every
slider change. When the checkbox is disengaged, results are only
recomputed when the Apply button is pressed.
The result graphs will automatically update to reflect changes made
with the customizable editor. Updated results automatically replace
the results shown in each graph as soon as a recompute is complete.
Sometimes you may wish to abandon all the parameter adjustments
you have made and start over with the calibration. Press the Reset
button to return the parameter values in the customizable editor to the
values in place when the editor first opened.
371
Chapter 12 Hydrologic Simulation
372
Chapter 12 Hydrologic Simulation
373
Chapter 12 Hydrologic Simulation
Figure 287. Viewing the global summary table for a simulation run. Volume
in millimeters or inches is computed by dividing the volume by
the drainage area.
You can choose the units for volume results. Actual volume is shown
in either thousands of cubic meters (THOU M3) or in acre-feet (AC-
FT). You can use the buttons above the table of results to change
volume units.
You can view the elements listed either in alphabetic or hydrologic
order. When hydrologic order is selected, the elements are shown in
the same order as in the Watershed Explorer. This order represents
the order determined from the flow network plus any manual
adjustments. Optionally you can switch to viewing elements in
alphabetic order using the selection box above the table of results.
Individual Elements
There are three ways to access results for individual elements in the
current simulation run. All methods show the same results. You must
select one or more elements in the basin map before you can view
results for that element.
The first method for viewing individual element results is to use the
Results menu. You may need to first open the basin model by clicking
on it in the Watershed Explorer, on the "Components" tab. The current
simulation run is shown in brackets in the basin map title bar. Select
one or more elements in the basin map by clicking with the arrow tool.
With an element selected in the basin map, click the Results menu
374
Chapter 12 Hydrologic Simulation
The information included in the summary table (Figure 289) also varies
by element type but always includes the peak flow, time of peak flow,
and outflow volume. The time-series table includes the same
information as the graph but in numerical format (Figure 290). If you
selected more than one element, then one result for each element will
open. The menu commands on the Results menu will only be enabled
if the results for the current simulation run do not need to be
recomputed, and there is an element selection in the basin map.
375
Chapter 12 Hydrologic Simulation
Figure 290. Element time-series table for a junction with observed flow.
The second method for viewing individual element results is to use the
element icon in the basin map. Again you may need to open the basin
model and select a current simulation run before attempting to view
results. Click on an element to highlight it. Keep the mouse over the
element and press the right mouse button. A context menu is
376
Chapter 12 Hydrologic Simulation
377
Chapter 12 Hydrologic Simulation
summary table will open. It is exactly the same table that can be
viewed for the current simulation run. The table includes one row for
each element in the basin model and columns for element name,
drainage area, peak flow, time of peak flow, and total outflow volume.
Figure 291. Accessing simulation run results from the Watershed Explorer.
Individual Elements
Each element in the basin model is shown in the Watershed Explorer
under the simulation run node. These elements are listed in hydrologic
order below the global summary table. The results for each element
are accessed by clicking on its node. The first item listed for each
element is the graph; click on the "Graph" node to view the result
(Figure 291). It is exactly the same graph that can be viewed for the
current simulation run. The information included in the graph varies by
element type but always includes outflow. Optional items such as
observed flow, computed stage, and observed stage are also included.
Similarily, the summary table and time-series table can also be
378
Chapter 12 Hydrologic Simulation
379
Chapter 12 Hydrologic Simulation
Figure 292. Selecting computed outflow from the same element in two
different simulation runs. Other types of time-series data could
also be selected.
380
Chapter 12 Hydrologic Simulation
After you have opened a time-series table or graph, you may add
additional time-series results. Position the mouse over the time-series
result in the Watershed Explorer that you wish to add to the graph or
table. Press and hold the left mouse button and then drag the mouse
over the top of the graph or table where you want the result to be
added. The mouse cursor will change to indicate which tables and
graphs can accept the additional time-series. Release the mouse
button while it is over the desired table or graph and it will be
automatically updated to show the additional time-series results.
381
Chapter 12 Hydrologic Simulation
changed. It is also possible to change the properties for the axis, title,
gridlines, patterns, and legend.
382
Chapter 13 Model Optimization
CHAPTER 13
Model Optimization
This chapter describes how observed streamflow can be used to
optimize model performace by automatically estimating parameters.
The optimization process begins with initial parameter estimates and
adjusts them so that the simulated results match the observed
streamflow as closely as possible. Two different search algorithms are
provided that move from the initial parameter estimates to the final best
parameter estimates. A variety of objective functions are provided to
measure the goodness of fit between the simulated and observed
streamflow in different ways. While parameter estimation using
optimization does not produce perfect results, it can be a valuable aid
when calibrating models.
Optimization Trials
Parameter estimation is the process of adapting a general model to a
specific watershed. Some parameters can be estimated directly from
field measurements. For example, the area that must be entered for a
subbasin element can be measured directly in the field using standard
surveying procedures or from maps developed through surveying.
Other parameters can be estimated indirectly from field
measurements. In this case, the field measurement does not result in
a value that can be input directly to the program. However, the field
measurement can provide a strong recommendation for a parameter in
the program based on previous experience. For example,
measurements of soil texture are correlated with parameters such as
hydraulic conductivity. Finally, there are parameters that can only be
estimated by comparing computed results to observations such as
measured streamflow. Even for parameters of the first two types, there
is often enough uncertainty in the true parameter value to require some
adjustment of the estimates in order for the model to closely follow the
observed streamflow.
The quantitative measure of the goodness-of-fit between the computed
result from the model and the observed flow is called the objective
function. An objective function measures the degree of variation
between computed and observed hydrographs. It is equal to zero if
the hydrographs are exactly identical. The key to automated
parameter estimation is a search method for adjusting parameters to
minimize the objective function value and find optimal parameter
values. A minimum objective function is obtained when the parameter
383
Chapter 13 Model Optimization
384
Chapter 13 Model Optimization
Figure 295. Beginning the process of creating a new optimization trial using
the Optimization Trial Manager.
The first step of creating an optimization trial is to provide the name for
the new trial (Figure 296). A default name is provided for the new
optimization trial; you can use the default or replace it with your own
choice. After you finish creating the trial you can add a description to
it. If you change your mind and do not want to create a new
optimization trial, you can press the Cancel button at the bottom of the
wizard or the X button in the upper right corner of the wizard. The
Cancel button can be pressed at any time you are using the wizard.
Press the Next> button when you are satisfied with the name you have
entered and are ready to proceed to the next step.
385
Chapter 13 Model Optimization
Figure 296. Entering a name for a new optimization trial. The remaining
steps are to select a basin model, the optimization location, and
a meteorologic model.
386
Chapter 13 Model Optimization
387
Chapter 13 Model Optimization
The second way to copy is from the "Compute" tab of the Watershed
Explorer. Move the mouse over the optimization trial you wish to copy
and press the right mouse button (Figure 298). A context menu is
displayed that contains several choices including copy. Click the
Create Copy command. The Copy Optimization Trial window will
open where you can name and describe the copy that will be created.
A default name is provided for the copy; you can use the default or
replace it with your own choice. A description can also be entered; if it
is long you can use the button to the right of the description field to
open an editor. When you are satisfied with the name and description,
press the Copy button to finish the process of copying the selected
optimization trial. You cannot press the Copy button if no name is
specified. If you change your mind and do not want to copy the
selected optimization trial, press the Cancel button or the X button in
the upper right of the Copy Optimization Trial window to return to the
Watershed Explorer.
Figure 298. Copying an optimization trial from the Watershed Explorer. The
Copy Optimization Trial window will appear after the Create
Copy menu command is selected.
388
Chapter 13 Model Optimization
Figure 299. Renaming an optimization trial. This was accessed from the
Optimization Trial Manager.
389
Chapter 13 Model Optimization
390
Chapter 13 Model Optimization
The second way to delete is from the "Compute" tab of the Watershed
Explorer. Select the optimization trial you wish to delete by clicking on
it in the Watershed Explorer; it will become highlighted (Figure 302).
Keep the mouse over the selected trial and click the right mouse
button. A context menu is displayed that contains several choices
including delete. Click the Delete command. A window will open
where you must confirm that you wish to delete the selected trial.
Press the OK button to delete the trial. If you change your mind and
do not want to delete the selected optimization trial, press the Cancel
button or the X button in the upper right to return to the Watershed
Explorer.
391
Chapter 13 Model Optimization
Selecting Components
One of the principal tasks when creating an optimization trial using the
wizard is the selection of the components that will be used to compute
optimization results. The components include the basin model and the
hydrologic element in the basin model where the objective function will
be computed. The components also include the meteorologic model.
These components are selected when creating a new optimization trial
with the wizard. However, you can change the basin model and
meteorologic model you wish to use at any time using the Component
Editor for the optimization trial. Access the Component Editor from the
"Compute" tab of the Watershed Explorer (Figure 303). If necessary,
click on the "Optimization Trials" folder to expand it and view the
available optimization trials. The Component Editor contains a basin
model selection list that includes all of the basin models in the project
where the basin model has at least one element with observed flow.
The Component Editor also contains a meteorologic model selection
list that includes all of the meteorologic models in the project. The
element where the objective function will be computed can also be
changed, but that function is performed in a different Component Editor
that will be described later.
392
Chapter 13 Model Optimization
Figure 303. The optimization trial component editor can be used to change
the selected components, time window, and search tolerance.
Search Method
Two search methods are available for minimizing the objective function
and finding optimal parameter values. The univariate gradient method
evaluates and adjusts one parameter at a time while holding other
parameters constant. The Nelder and Mead method uses a downhill
simplex to evaluate all parameters simultaneously and determine
which parameter to adjust. The default method is the univariate
gradient method. Also, the univariate gradient method is always used
if only one parameter is selected.
Selecting the search method for the optimization trial is accessed from
the "Compute" tab of the Watershed Explorer in the optimization trial
Component Editor. Click on the optimization trial node to display the
Component Editor for the optimization trial (Figure 303). If necessary,
393
Chapter 13 Model Optimization
Objective Function
The objective function measures the goodness-of-fit between the
computed outflow and observed streamflow at the selected element.
Eight different functions are provided that measure the goodness-of-fit
in different ways. The RMS error function calculates the classical root
mean square error. The peak-weighted RMS error function is a
modification of the standard root mean square error that gives greatly
increased weight to flows above average and less weight to flows
below average. The RMS log error uses the root-mean-square error of
the logarithm-transformed flows to emphasize small flows compared to
large flows. The sum of squared residuals function gives increased
weight to large errors and less weight to small errors. The sum of
absolute residuals function gives equal weight to large and small
errors. The percent error in peak flow function ignores the entire
hydrograph except for the single peak flow value. The percent error in
volume function ignores peak flow or timing considerations in favor of
the volume. The time-weighted function gives greater weight to error
near the end of the optimization time window and less weight to errors
early in the window.
Select the objective function method for the optimization trial on the
Component Editor for the objective function (Figure 304). Access the
Component Editor from the "Compute" tab of the Watershed Explorer.
If necessary, click on the "Optimization Trials" folder to expand it and
view the available optimization trials in the project. Click on the
optimization trial node to expand it and see the objective function node.
394
Chapter 13 Model Optimization
Click on the objective function node to view the editor. The objective
function method is the first item in the Component Editor.
Figure 304. Specifying properties for the objective function using the
component editor.
395
Chapter 13 Model Optimization
396
Chapter 13 Model Optimization
397
Chapter 13 Model Optimization
Method Parameter
Dynamic Initial Storage
Gridded Simple Initial Storage
Simple Initial Storage
Maximum Storage
Method Parameter
Gridded Simple Initial Storage
Simple Initial Storage
Maximum Storage
Method Parameter
Curve Number Initial Abstraction
Constant Rate
Deficit Constant Initial Deficit
Maximum Deficit
Constant Rate
Exponential Initial Range
Initial Coefficient
Coefficient Ratio
Exponent
Green Ampt Initial Content
Saturated Content
Suction
Conductivity
Gridded Curve Initial Abstraction Ratio
Number Potential Retention Scale Factor
398
Chapter 13 Model Optimization
Method Parameter
Gridded Deficit Initial Deficit Ratio
Constant Maximum Deficit Ratio
Constant Rate Ratio
Impervious Area Ratio
Gridded SMA Soil Initial Storage
Groundwater 1 and 2 Initial Storage
Initial Constant Initial Loss
Constant Rate
Smith Parlange Initial Content
Residual Content
Saturated Content
Bubbling Pressure
Pore Distribution
Conductivity
Beta Zero
Soil Moisture Soil Initial Storage
Accounting Soil Storage
Soil Tension Storage
Soil Percolation
Groundwater 1 and 2 Initial Storage
399
Chapter 13 Model Optimization
Method Parameter
Clark Time of Concentration
Storage Coefficient
Kinematic Wave Plane Roughness
Subcollector, Collector, and
Main Channel Manning's n
ModClark Time of Concentration
Storage Coefficient
SCS Time Lag
S-Graph Time Lag
Snyder Peaking Coefficient
Standard Lag
Method Parameter
Bounded Initial Flow Rate or Initial Flow Rate per Area 1
Recession Recession Constant
Linear Reservoir GW 1 and 2 Initial Flow Rate or
Initial Flow Rate per Area 1
Groundwater 1 and 2 Storage Coefficient
Groundwater 1 and 2 Number of Steps
Nonlinear Initial Flow Rate or Initial Flow Rate per Area 1
Boussinesq Characteristic Length
Hydraulic Conductivity
Drainable Porosity
Threshold Ratio or Threshold Flow Rate 2
Recession Initial Flow Rate or Initial Flow Rate per Area 1
Recession Constant
Threshold Ratio or Threshold Flow Rate 2
1 The available parameter depends on the method selected for specifying the initial condition.
2 The available parameter depends on the method selected for specifying the recession threshold.
400
Chapter 13 Model Optimization
Method Parameter
Kinematic Wave Manning's n
Lag Lag
Modified Puls Subreaches
Initial Flow
Muskingum K
X
Subreaches
Muskingum Cunge Manning's n
Straddle Stager Lag
Duration
401
Chapter 13 Model Optimization
The initial value is the starting point for the parameter estimation
process. The search method will begin searching from that point for
optimal parameter values. The default initial value is the parameter
value in the basin model that was selected for the optimization trial.
You may change the initial value without affecting the basin model.
It is possible to lock a parameter. When a parameter is locked, the
initial value is used and no adjustments are made during the search
process.
The minimum parameter value can be used to narrow the lower end of
the range of values that will be used by the search method. Likewise,
the maximum parameter value can be used to narrow the upper end of
the range of values that will be used by the search method. A good
source of information for narrowing the search range is preliminary
estimates from field measurements or manual calibration. Default
values for the minimum and maximum are provided based on physical
and numerical limits. The search may continue outside the specified
range. When it does so, a penalty is applied that is proportional to the
distance outside the specified range. The penalty nudges the search
for optimal parameter values back to the range between the specified
minimum and maximum.
402
Chapter 13 Model Optimization
accessed from the Watershed Explorer and the basin map as soon as
they are available. The current optimization trial is shown in the title of
the basin map.
The first way to select the current optimization trial is from the toolbar.
The compute selection box shows all of the compute components in
the project; optimization trials appear after simulations runs in the
selection box. Each optimization trial in the selection box includes the
prefix Trial followed by the name of the optimization trial. Click on
the toolbar selection box and choose an optimization trial to become
the current trial.
The second way to select the current optimization trial is from the
Watershed Explorer, on the "Compute" tab. Select the optimization
trial you wish to become the current trial by clicking on it in the
Watershed Explorer. The optimization trials are listed alphabetically in
the "Optimization Trials" folder.
403
Chapter 13 Model Optimization
404
Chapter 13 Model Optimization
405
Chapter 13 Model Optimization
time, and end time. It indicates which objective function method was
used, the start and end of the evaluation period for the function, and
the final value of the function when the search method finished. It
gives the volume, peak flow, time of peak flow, and time to the center
of mass of the computed and observed hydrographs. Finally, it
provides volume and peak flow differences between the computed and
observed hydrographs. The objective function table is accessed from
the Results menu. When the Basin Map is open and a trial is
selected, choose the Objective Function Summary command to view
the results. Sample results are shown in Figure 309.
Figure 309. Objective function table for an optimization trial. The table
presents statistical results of the computed and observed flow at
the optimization location.
406
Chapter 13 Model Optimization
407
Chapter 13 Model Optimization
408
Chapter 13 Model Optimization
Figure 312. Flow comparison graph showing the computed flow plotted
against the observed flow. Different symbols are used for points
before and after the observed peak flow.
409
Chapter 13 Model Optimization
Figure 313. Flow residuals showing the computed minus the observed flow
at the optimization location.
410
Chapter 13 Model Optimization
Figure 314. Objective function graph showing the value of the objective
function after each iteration in the search.
Individual Elements
There are three ways to access results for individual elements in the
current optimization trial. All methods show the same results. You
must select one or more elements in the basin map before you can
view results for that element. However, not all elements in the basin
map may have results. The optimization trial is based on a simulation
run that uses a specific basin model. The optimization trial works with
the elements upstream of the location where observed flow is used to
evaluate the objective function. The trial does not perform simulations
at any elements downstream of the element with observed flow.
Therefore, results will only be available at elements upstream of the
location with observed flow. Any downstream elements will have
results shown as unavailable.
The first method for viewing individual element results is to use the
Results menu. You may need to first open the basin model by clicking
on it in the Watershed Explorer, on the "Components" tab. The current
optimization trial is shown in brackets in the basin map title bar. Select
one or more elements in the basin map by clicking with the arrow tool.
With an element selected in the basin map, click the Results menu
and select the Element Graph command, the Element Summary
Table command, or the Element Time-Series Table command. The
appropriate result will automatically be shown in the Desktop area.
411
Chapter 13 Model Optimization
412
Chapter 13 Model Optimization
413
Chapter 13 Model Optimization
Figure 316. Element summary table for a subbasin after an optimization trial.
414
Chapter 13 Model Optimization
The third method for viewing individual element results is to use the
buttons on the toolbar. First you must open the basin model and select
one or more elements in the basin map by clicking with the arrow tool.
Once you have a selected element, click the graph, summary table, or
time-series table buttons. The button for selecting a graph shows a
line plot. The button for the summary table shows a plain table. The
button for the time-series table shows a table plus a clock. The
appropriate result will automatically open. If you selected more than
one element, then one result for each element will be added to the
Desktop. The toolbar buttons for viewing element results will only be
enabled if the results for the current optimization trial do not need to be
recomputed, and there is an element selection in the basin map, and
the selected element is upstream of the objective function evaluation
element with observed flow.
Trial Results
The objective function summary can be accessed by clicking on the
"Objective Function Summary" node in the Watershed Explorer (Figure
318). The table will be displayed. It is exactly the same table that can
be viewed for the current optimization trial. Similarly, the other results
for optimized parameters, hydrograph comparison, flow comparison,
flow residuals, and objective function can also be accessed in the
same way.
415
Chapter 13 Model Optimization
Individual Elements
Each element above the objective function evaluation location is
shown in the Watershed Explorer under the optimization trial node.
These elements are listed in hydrologic order after the objective
function graph. The results for each element are accessed by clicking
on its node. The first item listed for each element is the graph; click on
the "Graph" node to view the result (Figure 318). It is exactly the same
graph that can be viewed for the current optimization trial. The
information included in the graph varies by element type but always
includes outflow. Optional items such as observed flow, computed
stage, and observed stage are also included. Similarily, the summary
table and time-series table can also be accessed by clicking on the
"Summary Table" or "Time-Series Table" node, respectively.
416
Chapter 13 Model Optimization
417
Chapter 13 Model Optimization
After you have opened a time-series table or graph, you may add
additional time-series results. Position the mouse over the time-series
result you wish to add to the graph or table. Press and hold the left
mouse button and then drag the mouse over the top of the graph or
table where you want the result to be added. The mouse cursor will
change to indicate which tables and graphs can accept the additional
time-series. Release the mouse button while it is over the desired
table or graph and it will be automatically updated to show the
additional time-series results.
418
Chapter 13 Model Optimization
419
Chapter 13 Model Optimization
420
Chapter 14 Forecasting Streamflow
CHAPTER 14
Forecasting Streamflow
Streamflow forecasts are fundamental to effective operation of flood
control reservoirs and levee systems. Forecasts may also support
emergency operations by providing estimates of the timing and extent
of expected hazardous or damaging flood conditions. Forecasts are
based on recent meteorologic and hydrologic conditions in the
watershed, and may also incorporate predicted future meteorologic
conditions. Although most often used to predict anticipated flood
conditions, streamflow forecasts may also support water supply,
hydropower, environmental flow requirements, and other operational
needs.
Forecast Alternatives
Streamflow forecasting typically includes simulation of past and future
conditions. The process begins with the selection of a time of forecast.
Typically the time of forecast represents the last available time for
meteorologic observations of precipitation, temperature, and other
variables. If observations of streamflow, stage, or reservoir pool
elevation are available, the last available value is also generally near
the time of forecast. The simulation is started hours or days before the
time of forecast. Results computed between the start time and the
forecast time may be called the "look back period." When
observations of current watershed conditions are available, they may
be compared with computed results from the look back period to make
calibration adjustments that improve model performance. Meteorologic
observations are not available after the time of forecast and predictions
of future values are used. For example, a quantitative precipitation
forecast (QPF) provides a meteorologic prediction of future
precipitation depths. Similar predictions are used for other
meteorologic variables such as temperature. The future streamflow
response is simulated based on the predicted meteorologic conditions.
This period of time in the future may be called the "forecast."
Time is often of the essence when performing streamflow forecasting
and modeling techniques must consider this fact. During the
development of a forecast it is common to group similar subbasins
together and apply parameter adjustments across each group. This is
supported in forecast alternatives through the use of zones. Creating
zones is described in Chapter 6. Zones are used in a forecast
alternative to quickly set parameter adjustments when a group of
421
Chapter 14 Forecasting Streamflow
422
Chapter 14 Forecasting Streamflow
423
Chapter 14 Forecasting Streamflow
Figure 322. Entering a name for a new forecast alternative. The remaining
steps are to select a basin model, select a meteorologic model,
and select a time interval and zone configurations.
424
Chapter 14 Forecasting Streamflow
425
Chapter 14 Forecasting Streamflow
The second way to copy is from the "Compute" tab of the Watershed
Explorer. Move the mouse over the forecast alternative you wish to
copy and press the right mouse button (Figure 324). A context menu
is displayed that contains several choices including copy. Click the
Create Copy command. The Copy Forecast Alternative window will
open where you can name and describe the copy that will be created.
A default name is provided for the copy; you can use the default or
replace it with your own choice. A description can also be entered; if it
is long you can use the button to the right of the description field to
open an editor. When you are satisfied with the name and description,
press the Copy button to finish the process of copying the selected
forecast alternative. You cannot press the Copy button if no name is
specified. If you change your mind and do not want to copy the
selected forecast alternative, press the Cancel button or the X button
in the upper right of the Copy Forecast Alternative window to return to
the Watershed Explorer.
426
Chapter 14 Forecasting Streamflow
Figure 325. Renaming a forecast alternative. This was accessed from the
Forecast Alternative Manager.
427
Chapter 14 Forecasting Streamflow
change your mind while in editing mode and do not want to rename the
selected forecast alternative, press the Escape key.
428
Chapter 14 Forecasting Streamflow
The second way to delete is from the "Compute" tab of the Watershed
Explorer. Select the forecast alternative you wish to delete by clicking
on it in the Watershed Explorer; it will become highlighted (Figure 328).
Keep the mouse over the selected alternative and click the right mouse
button. A context menu is displayed that contains several choices
including delete. Click the Delete command. A window will open
where you must confirm that you wish to delete the selected
alternative. Press the OK button to delete the alternative. If you
change your mind and do not want to delete the selected forecast
alternative, press the Cancel button or the X button in the upper right
to return to the Watershed Explorer.
429
Chapter 14 Forecasting Streamflow
Selecting Components
One of the principal tasks when creating a forecast alternative using
the wizard is the selection of the components that will be used to
compute forecast results. The components include the basin model
and the meteorologic model. These components are selected when
creating a new forecast alternative with the wizard. However, you can
change the basin model and meteorologic model you wish to use at
any time using the Component Editor for the forecast alternative.
Access the Component Editor from the "Compute" tab of the
Watershed Explorer (Figure 329). If necessary, click on the "Forecast
Alternatives" folder to expand it and view the available forecast
alternatives. The Component Editor contains a basin model selection
list that includes all of the basin models in the project. The Component
Editor also contains a meteorologic model selection list that includes all
of the meteorologic models in the project.
430
Chapter 14 Forecasting Streamflow
431
Chapter 14 Forecasting Streamflow
432
Chapter 14 Forecasting Streamflow
Figure 330. Entering parameter adjustments for the Alligator Creek zone.
The initial deficit will be decreased by 5 mm at each subbasin,
the maximum deficit will be multiplied by 1.25 at each subbasin,
and the constant rate will be increased by 2 mm/hr at each
subbasin.
The bottom half of the parameter editor (Figure 330) shows details
about the parameter values for the elements in the selected zone.
Clicking on a row in the upper half of the parameter editor selects a
zone. The lower half of the editor will show one row for each element
in the selected zone. Two columns are shown for each parameter in
the loss rate, transform, baseflow, or routing method. The first column
labeled Base Value shows the default parameter value taken from the
basin model. The second column labeled Final Value shows the
adjusted parameter value that will be used during the forecast.
Adjusted parameter values are calculated using any zone adjustments
entered in the upper half of the editor. The adjusted parameter values
are updated everytime the Apply button is pressed.
433
Chapter 14 Forecasting Streamflow
Figure 331. Entering an override for the AlgCr-02 subbasin within the
Alligator Creek zone. The default final value has been changed
for the Initial Deficit parameter.
434
Chapter 14 Forecasting Streamflow
Slider Adjustments
A slider adjustment is an alternate way to view and change a zone
parameter adjustment or an element parameter override. An example
of a slider adjustment is shown in Figure 332. Instead of typing a new
value for the adjustment in a data entry table, the slider can be moved
with the mouse to quickly obtain a new factor. Slider adjustments
enhance existing zone parameter adjustments and element parameter
overrides. Configuration information is required for each adjustment or
override that will be manipulated with a slider. Preparing to use a
slider adjustment begins with selecting the adjustments or override and
entering configuration information.
Figure 332. Using the slider adjustments editor to rapidly change the zone
adjustments for loss rate parameters.
435
Chapter 14 Forecasting Streamflow
row and multiple elements in multiple zones are all shown together.
The tables will be empty if no adjustments or overrides have been
selected for slider adjustments.
436
Chapter 14 Forecasting Streamflow
zone currently selected in the upper half. Click on a zone name in the
upper half to select it; the selected zone will remain highlighted until a
different zone is selected. Move to the bottom half and select an
element by clicking on it. The available element parameter overrides
will be shown after an element is selected. Overrides are selected in
the same way as zone parameter adjustments by pressing the Select
button or by double clicking on an override. Press the Close button
when you have finished selecting element override adjustments.
437
Chapter 14 Forecasting Streamflow
Figure 335. Preparing to remove two slider adjustments. Hold the control
key and click multiple adjustments before pressing the Remove
button.
Selecting slider adjustments is the first step in preparing them for use.
Each slider adjustment must also be configured with the properties
necessary to operate properly in the slider editor (Figure 332). Access
the settings for the current forecast alternative by selecting the
Compute menu and clicking the Forecast Slider Settings
command. The table in the upper half shows the slider adjustments for
zone parameters (Figure 336). The first column displays the zone
name while the second column shows the parameter name. The next
column shows the selected adjustment type. The type may be a value,
increment, or factor. The units of each zone parameter are shown for
reference. The minimum value sets the left side of the slider range
while the maximum sets the right side of the range. The minimum and
maximum values should be specified carefully giving consideration to
the nature of each parameter.
The table in the lower half shows the slider adjustments for element
parameter overrides (Figure 336). The zone name, element name,
and parameter name is shown for each slider adjustment. Use the
same approach for entering the settings for element parameter
overrides as for zone parameters. The slider adjustment settings can
be closed when all configurations are complete.
438
Chapter 14 Forecasting Streamflow
439
Chapter 14 Forecasting Streamflow
440
Chapter 14 Forecasting Streamflow
The observed flow from the selected gage is multiplied by the specified
ratio in order to calculate the initial baseflow estimated value. The ratio
is only used if the observed flow is between the minimum and
maximum values. It is required for the ratio to be greater than zero
and it may be greater than one. However, using a ratio greater than
one should generally only be done if the observed flow is not
downstream of the subbasin; this should be done only after careful
analysis.
The estimated value is calculated from the observed flow gage using
the ratio and substituting the default value according to the minimum
and maximum values. At first, the estimated value is also copied as
the final value. The final value can be edited to override the estimated
value. Editing the final value may be necessary if current flow
conditions diverge from the assumptions used to develop the ratio.
The final value for each subbasin shown in the table will be reset to the
estimated value when the Reset button is pressed.
441
Chapter 14 Forecasting Streamflow
Replace All With Zero choice will substitute zero discharge for
all time intervals after the forecast time even if some observed or
specified values may be availalable from the time-series manger. The
Replace Missing With Zero choice will substitute zero for any
missing values between the forecast time and end time. The Replace
Missing With Last choice will substitute the last observed or
specified value prior to the missing value during the period from the
forecast time to the end time.
Figure 338. Setting reservoir initial pool elevation at the forecast start time.
The Forecast Reservoirs editor uses two tabs for separating initial
pool elevation at the start time from resetting the pool elevation at the
forecast time. Each reservoir in the basin model is shown on a
separate row in the editor. The "Start Time" tab corresponds to the
start time of the forecast alternative while the "Forecast Time" tab is for
the forecast time of the alternative.
The "Start Time" tab is where a stage gage may be selected for each
reservoir. Pool elevation is taken from a stage gage that must be
created in the time-series gage manager before it can be selected for a
reservoir. The selected gage should represent a pool elevation
observation at the reservoir. The observed elevation will be shown for
the start time of the forecast alternative.
A minimum and maximum elevation may be entered to qualify the use
of the observed stage gage for calculating the initial storage. The
minimum and maximum values are used in conjunction with the default
value. If the observed value is less than the minimum, then the default
will be used instead. The default will also be used if the observed
value is greater than the maximum. The minimum or maximum can be
left blank if the observed value should be used without a lower or
442
Chapter 14 Forecasting Streamflow
Figure 339. Reseting the reservoir pool elevation at the forecast time.
The "Forecast Time" tab allows the reservoir storage to be reset at the
forecast time of the alternative. The selection of the stage gage,
minimum value, maximum value, and default value are shown for each
reservoir. However, these properties are set on the "Start Time" tab
and are only repeated for easy reference.
The "Forecast Time" tab shows the observed elevation at the start time
of the forecast alternative. The observed value will be blank if the data
is missing at the forecast time. The estimated value is calculated from
the observed stage gage by substituting the default value according to
the minimum and maximum values. The default value is also
substituted if the observed value is missing. At first, the estimated
value is also copied as the final value. The final value can be edited to
override the estimated value. The final value for each reservoir shown
443
Chapter 14 Forecasting Streamflow
in the table will be reset to the estimated value when the Reset button
is pressed.
444
Chapter 14 Forecasting Streamflow
Figure 340. Setting the blending properties for each element with observed
flow.
445
Chapter 14 Forecasting Streamflow
The first way to select the current forecast alternative is from the
toolbar. The compute selection box shows all of the compute
components in the project; forecast alternatives appear after
optimization trials in the selection box. Each forecast alternative in the
selection box includes the prefix Forecast followed by the name of
the alternative. Click on the toolbar selection box and choose a
forecast alternative to become the current alternative.
The second way to select the current forecast alternative is from the
Watershed Explorer, on the "Compute" tab. Select the forecast
alternative you wish to become the current alternative by clicking on it
in the Watershed Explorer. The forecast alternatives are listed
alphabetically in the "Forecast Alternatives" folder.
446
Chapter 14 Forecasting Streamflow
447
Chapter 14 Forecasting Streamflow
448
Chapter 14 Forecasting Streamflow
449
Chapter 14 Forecasting Streamflow
Individual Elements
There are three ways to access results for individual elements in the
current forecast alternative. All methods show the same results. You
must select one or more elements in the basin map before you can
view results for that element.
The first method for viewing individual element results is to use the
Results menu. You may need to first open the basin model by clicking
on it in the Watershed Explorer, on the "Components" tab. The current
forecast alternative is shown in brackets in the basin map title bar.
Select one or more elements in the basin map by clicking with the
arrow tool. With an element selected in the basin map, click the
Results menu and select the Element Graph command, the Element
Summary Table command, or the Element Time-Series Table
command. The appropriate result will automatically be shown in the
Desktop area. The information included in the graph (Figure 343)
varies by element type, but always includes outflow. Optional items
such as observed flow, computed stage, and observed stage are also
included.
450
Chapter 14 Forecasting Streamflow
The information included in the summary table (Figure 344) also varies
by element type but always includes the peak flow, time of peak flow,
and outflow volume. The time-series table includes the same
information as the graph but in numerical format (Figure 345). If you
selected more than one element, then one result for each element will
open. The menu commands on the Results menu will only be enabled
if the results for the current forecast alternative do not need to be
recomputed, and there is an element selection in the basin map.
Figure 344. Element summary table for a junction after a forecast alternative.
451
Chapter 14 Forecasting Streamflow
The second method for viewing individual element results is to use the
element icon in the basin map. Again you may need to open the basin
model and select a current forecast alternative before attempting to
view results. Click on an element to highlight it. Keep the mouse over
the element and press the right mouse button. A context menu is
displayed that contains several choices including View Results. The
name of the current forecast alternative is shown in brackets after the
menu command. Move the mouse into the View Results submenu
and select any of the result commands: Graph, Summary Table, or
Time-Series Table. The appropriate result will open. If you selected
more than one element, then one result for each element will open.
The View Results menu command will only be enabled if the results
for the current forecast alternative do not need to be recomputed.
The third method for viewing individual element results is to use the
buttons on the toolbar. First you must open the basin model, select a
current forecast alternative, and select one or more elements in the
basin map by clicking with the arrow tool. Once you have a selected
element, click the graph, summary table, or time-series table buttons.
The button for selecting a graph shows a line plot. The button for the
summary table shows a plain table. The button for the time-series
table shows a table plus a clock. The appropriate result will
452
Chapter 14 Forecasting Streamflow
automatically open. If you selected more than one element, then one
result for each element will be added to the Desktop. The toolbar
buttons for viewing element results will only be enabled if the results
for the current forecast alternative do not need to be recomputed, and
there is an element selection in the basin map.
Individual Elements
Each element in the basin model used in the forecast alternative is
shown in the Watershed Explorer under the forecast alternative node.
These elements are listed in hydrologic order. The results for each
element are accessed by clicking on its node. The first item listed for
each element is the graph; click on the "Graph" node to view the result
(Figure 346). It is exactly the same graph that can be viewed for the
current forecast alternative. The information included in the graph
varies by element type but always includes outflow. Optional items
such as observed flow, computed stage, and observed stage are also
included. Similarily, the summary table and time-series table can also
be accessed by clicking on the "Summary Table" or "Time-Series
Table" node, respectively.
453
Chapter 14 Forecasting Streamflow
454
Chapter 14 Forecasting Streamflow
455
Chapter 14 Forecasting Streamflow
After you have opened a time-series table or graph, you may add
additional time-series results. Position the mouse over the time-series
result you wish to add to the graph or table. Press and hold the left
mouse button and then drag the mouse over the top of the graph or
table where you want the result to be added. The mouse cursor will
change to indicate which tables and graphs can accept the additional
time-series. Release the mouse button while it is over the desired
table or graph and it will be automatically updated to show the
additional time-series results.
456
Chapter 14 Forecasting Streamflow
457
Chapter 14 Forecasting Streamflow
458
Chapter 15 Depth-Area Reduction
CHAPTER 15
Depth-Area Reduction
One of the principal results in many water resource studies is a flow-
frequency curve. A common way to develop a frequency curve is to
use frequency precipitation events in combination with a hydrologic
model. A study of a large watershed usually includes multiple
evaluation points in the stream network. The depth-area analysis
automates the process of producing flow estimates due to frequency
precipitation at multiple points of interest within a watershed.
Depth-Area Analyses
A depth-area analysis is designed to assist in the process of
developing peak flows in a watershed due to frequency precipitation.
For example, it may be the goal to compute the streamflow expected to
result from the storm with 1% exceedance probability. The usual
assumption is that the 1% precipitation will result in the 1% flow, and
that the hydrology model can accurately simulate the conversion of
precipitation to flow. Developing the flow values throughout a
watershed for a specific exceedance probability is a key step in
developing flow frequency curves.
Storm area is a key component of a frequency storm. The area is
used to determine the appropriate depth-area reduction factor. The
depth-area reduction factor is used to account for the observation from
empirical data that average precipitation intensity decreases as the
area of a storm increases. Thus the average precipitation intensity of a
200 km2 storm is generally much less than the intensity of a 10 km2
storm. Proper determination of the flow at a given evaluation point in
the watershed requires that the frequency storm have an area equal to
the drainage area at that point. Failure to match the storm area to the
drainage area can result in significant under or over estimation of the
flow. It is simple enough to setup a simulation run with the
meteorologic model set to compute flow at a particular evaluation
point. However, it is often the case that the flow must be estimated at
many points in the watershed.
The depth-area analysis automates the process of developing a
separate frequency storm for each evaluation point. The analysis
requires the specification of an underlying simulation run that must use
a meteorologic model set to use the frequency storm precipitation
method. The analysis also includes the listing of various analysis
points where the flow should be evaluated. When the analysis
459
Chapter 15 Depth-Area Reduction
460
Chapter 15 Depth-Area Reduction
Figure 350. The first step of creating a new depth-area analysis is to enter a
name for the new analysis. The remaining step is to select a
simulation run that uses the frequency precipitation method.
461
Chapter 15 Depth-Area Reduction
The second way to copy is from the "Compute" tab of the Watershed
Explorer. Move the mouse over the depth-area analysis you wish to
copy and press the right mouse button (Figure 352). A context menu
462
Chapter 15 Depth-Area Reduction
463
Chapter 15 Depth-Area Reduction
Figure 353. Renaming a depth-area analysis. This was accessed from the
Depth-Area Analysis Manager.
464
Chapter 15 Depth-Area Reduction
The second way to delete is from the "Compute" tab of the Watershed
Explorer. Select the depth-area analysis you wish to delete by clicking
on it in the Watershed Explorer; it will become highlighted (Figure 356).
Keep the mouse over the selected analysis and click the right mouse
button. A context menu is displayed that contains several choices
including delete. Click the Delete command. A window will open
where you must confirm that you wish to delete the selected analysis.
465
Chapter 15 Depth-Area Reduction
Press the OK button to delete the analysis. If you change your mind
and do not want to delete the selected depth-area analysis, press the
Cancel button or the X button in the upper right to return to the
Watershed Explorer.
466
Chapter 15 Depth-Area Reduction
Figure 357. Selecting the simulation run used as the basis for the depth-area
analysis.
467
Chapter 15 Depth-Area Reduction
468
Chapter 15 Depth-Area Reduction
469
Chapter 15 Depth-Area Reduction
470
Chapter 15 Depth-Area Reduction
Figure 360. Peak flow summary table showing the computed peak flow for
each analysis point. Results correspond to the hypothetical
storm with an area matched to the drainage area.
Individual Elements
There are three ways to access results for individual elements in the
current depth-area analysis. All methods show the same results. You
must select one or more elements in the basin map before you can
view results for that element. However, not all elements in the basin
map may have results. The depth-area analysis is based on a
simulation run that uses a specific basin model. The depth-area
analysis works with selected analysis points. The results are separate
471
Chapter 15 Depth-Area Reduction
for each analysis point. Results for an analysis point are only
computed at the elements upstream of that analysis point; all other
elements in the basin model will have no results for that analysis point.
When an element is upstream of more than one analysis point, you will
need to indicate the analysis point for which results are desired.
The first method for viewing individual element results is to use the
Results menu. You may need to first open the basin model by clicking
on it in the Watershed Explorer, on the "Components" tab. The current
depth-area analysis is shown in brackets in the basin map title bar.
Choose the analysis point for which results will be viewed, by selecting
from the Results Select Point menu. The currently selected
analysis point is shown in brackets in the menu. Next select one or
more elements in the basin map by clicking with the arrow tool. With
an element selected in the basin map, click the Results menu and
select the Element Graph command, the Element Summary Table
command, or the Element Time-Series Table command. The
appropriate result will automatically be shown in the Desktop area.
The information included in the graph (Figure 361) varies by element
type, but always includes outflow. Optional items such as observed
flow, computed stage, and observed stage are also included.
472
Chapter 15 Depth-Area Reduction
The information included in the summary table (Figure 362) also varies
by element type but always includes the peak flow, time of peak flow,
and outflow volume. The time-series table (Figure 363) includes the
same information as the graph but in numerical format. If you selected
more than one element, then one result for each element will open.
The menu commands on the Results menu will only be enabled if the
results for the current depth-area analysis do not need to be
recomputed, and an analysis point has been selected on the Results
menu.
473
Chapter 15 Depth-Area Reduction
474
Chapter 15 Depth-Area Reduction
The second method for viewing individual element results is to use the
element icon in the basin map. Again you may need to open the basin
model before attempting to view results. You may also need to choose
the analysis point for which results will be viewed, by selecting from the
Results Select Point menu. Finally, click on an element to
highlight it. Keep the mouse over the element and press the right
mouse button. A context menu is displayed that contains several
choices including View Results. The name of the current depth-area
analysis is shown in brackets after the menu command. Move the
mouse into the View Results submenu and select any of the result
commands: Graph, Summary Table, or Time-Series Table. The
appropriate result will open. If you selected more than one element,
then one result will open for each element. The View Results menu
command will only be enabled if the results for the current depth-area
analysis do not need to be recomputed, and an analysis point has
been selected on the Results menu.
The third method for viewing individual element results is to use the
buttons on the toolbar. First you must open the basin model, select a
current analysis point on the Results menu, and finally select one or
more elements in the basin map by clicking with the arrow tool. Once
you have a selected element, click the graph, summary table, or time-
series table buttons. The button for selecting a graph shows a line
plot. The button for the summary table shows a plain table. The
button for the time-series table shows a table plus a clock. The
appropriate result will automatically open. If you selected more than
one element, then one result for each element will be added to the
Desktop. The toolbar buttons for viewing element results will only be
enabled if the results for the current depth-area analysis do not need to
be recomputed, and there is an element selection in the basin map,
and the selected element is upstream of the analysis point selected on
the Results menu.
475
Chapter 15 Depth-Area Reduction
Individual Elements
Each analysis point is shown as a folder in the Watershed Explorer
under the depth-area analysis node. The folder for each analysis point
includes all of the elements above that analysis point. The results for
each element are accessed by clicking on its node. The first item
listed for each element is the graph; click on the "Graph" node to view
the result (Figure 364). It is exactly the same graph that can be viewed
for the current depth-area analysis. The information included in the
graph varies by element type but always includes outflow. Optional
items such as observed flow, computed stage, and observed stage are
also included. Similarily, the summary table and time-series table can
476
Chapter 15 Depth-Area Reduction
477
Chapter 15 Depth-Area Reduction
478
Chapter 15 Depth-Area Reduction
mouse button and then drag the mouse over the top of the graph or
table where you want the result to be added. The mouse cursor will
change to indicate which tables and graphs can accept the additional
time-series. Release the mouse button while it is over the desired
table or graph and it will be automatically updated to show the
additional time-series results.
Figure 366. Comparing results from the same element at two different
analysis points. The depth-area reduction factor changes due to
different drainage area at the two analysis points.
479
Chapter 15 Depth-Area Reduction
480
Chapter 16 Assessing Model Uncertainty
CHAPTER 16
Uncertainty Analyses
Uncertainty assessment is the process of determining the total error in
the simulated watershed response, for example, the flow at the outlet.
The simulated flow at the outlet actually depends on many individual
components each with its own error. There is error in the meteorologic
data because it is generally impossible to accurately measure
precipitation at the same spatial and temporal scale as the land
surface processes. There is error in the models of the hydrologic
processes becauses it is generally impossible to include every possible
process at the scale it occurs, for example, animal burrows or plant
transpiration. There is error in the model parameter values because
the equations are solved at a scale ranging from meters up to whole
subbasins and area-average values must be used. The error in the
whole watershed response includes all of these individual errors and
the complex way in which they interact.
The watershed model is very complex and it may not be possible to
attribute the total error to individual components. An error in the
precipitation data may be compensated with a corresponding error in
the infiltration parameter values. An error in the mathematical
formulation of the infiltration or transpiration process may be
compensated with a corresponding error in the parameter values. An
assumption about process scale may lead to effective area-averaged
parameter values that do not match values that would be measured in
the field. It is usually not possible to determine the exact error in each
481
Chapter 16 Assessing Model Uncertainty
482
Chapter 16 Assessing Model Uncertainty
483
Chapter 16 Assessing Model Uncertainty
wizard. The Cancel button can be pressed at any time you are using
the wizard. Press the Next> button when you are satisfied with the
name you have entered and are ready to proceed to the next step.
Figure 368. Entering a name for a new uncertainty analysis. The remaining
steps are to select a basin model and a meteorologic model.
484
Chapter 16 Assessing Model Uncertainty
The second way to copy is from the "Compute" tab of the Watershed
Explorer. Move the mouse over the uncertainty analysis you wish to
copy and press the right mouse button (Figure 370) A context menu is
displayed that contains several choices including copy. Click the
Create Copy command. The Copy Uncertainty Analysis window will
open where you can name and describe the copy that will be created.
A default name is provided for the copy; you can use the default or
485
Chapter 16 Assessing Model Uncertainty
486
Chapter 16 Assessing Model Uncertainty
the description field to open an editor. When you are satisfied with the
name and description, press the Rename button to finish the process
of renaming the selected uncertainty analysis. You cannot press the
Rename button if no name is specified. If you change your mind and
do not want to rename the selected uncertainty analysis, press the
Cancel button or the X button in the upper right of the Rename
Uncertainty Analysis window to return to the Uncertainty Analysis
Manager window.
Figure 371. Renaming an uncertainty analysis. This was accessed from the
Uncertainty Analysis Manager.
487
Chapter 16 Assessing Model Uncertainty
488
Chapter 16 Assessing Model Uncertainty
The second way to delete is from the "Compute" tab of the Watershed
Explorer. Select the uncertainty analysis you wish to delete by clicking
on it in the Watershed Explorer; it will become highlighted (Figure 374).
Keep the mouse over the selected analysis and click the right mouse
button. A context menu is displayed that contains several choices
including delete. Click the Delete command. A window will open
where you must confirm that you wish to delete the selected analysis.
Press the OK button to delete the analysis. If you change your mind
and do not want to delete the selected uncertainty analysis, press the
Cancel button or the X button in the upper right to return to the
Watershed Explorer.
489
Chapter 16 Assessing Model Uncertainty
Selecting Components
One of the principal tasks when creating an uncertainty analysis using
the wizard is the selection of the components that will be used to
compute uncertainty results. The components include the basin model
and the meteorologic model. These components are selected when
creating a new uncertainty analysis with the wizard. However, you can
change the basin model and meteorologic model you wish to use at
any time using the Component Editor for the uncertainty analysis.
Access the Component Editor from the "Compute" tab of the
Watershed Explorer (Figure 375). If necessary, click on the
"Uncertainty Analyses" folder to expand it and view the available
uncertainty analyses. The Component Editor contains a basin model
selection list that includes all of the basin models in the project. The
Component Editor also contains a meteorologic model selection list
that includes all of the meteorologic models in the project.
490
Chapter 16 Assessing Model Uncertainty
491
Chapter 16 Assessing Model Uncertainty
492
Chapter 16 Assessing Model Uncertainty
493
Chapter 16 Assessing Model Uncertainty
Method Parameter
Dynamic Initial Storage
Gridded Simple Initial Storage
Simple Initial Storage
Maximum Storage
Method Parameter
Gridded Simple Initial Storage
Simple Initial Storage
Maximum Storage
Method Parameter
Curve Number Initial Abstraction
Constant Rate
Impervious Area
Deficit Constant Initial Deficit
Maximum Deficit
Constant Rate
Impervious Area
Exponential Initial Range
Initial Coefficient
Coefficient Ratio
Exponent
Impervious Area
Green Ampt Initial Content
Saturated Content
Suction
Conductivity
Impervious Area
Gridded Curve Number Initial Abstraction Ratio
Potential Retention Scale Factor
494
Chapter 16 Assessing Model Uncertainty
Method Parameter
Gridded Deficit Constant Initial Deficit Ratio
Maximum Deficit Ratio
Constant Rate Ratio
Impervious Area Ratio
Gridded SMA Soil Initial Storage
Groundwater 1 and 2 Initial Storage
Initial Constant Initial Loss
Constant Rate
Impervious Area
Smith Parlange Initial Content
Residual Content
Saturated Content
Bubbling Pressure
Pore Distribution
Conductivity
Beta Zero
Impervious Area
Soil Moisture Soil Initial Storage
Accounting
Soil Storage
Soil Tension Storage
Soil Infiltration
Soil Percolation
Groundwater 1 and 2 Initial Storage
Groundwater 1 and 2 Storage
Groundwater 1 and 2 Percolation
Groundwater 1 and 2 Coefficient
Impervious Area
495
Chapter 16 Assessing Model Uncertainty
Method Parameter
Clark Time of Concentration
Storage Coefficient
Kinematic Plane Length
Wave Plane Slope
Plane Roughness
Subcollector, Collector, and Main Channel Length
Subcollector, Collector, and Main Channel Slope
Subcollector, Collector, and Main Channel Manning's n
Subcollector, Collector, and Main Channel Bottom Width 1
Subcollector, Collector, and Main Channel Side Slope 2
ModClark Time of Concentration
Storage Coefficient
SCS Time Lag
S-Graph Time Lag
Snyder Peaking Coefficient
Standard Lag
1 Bottom width can be selected for cross section shapes: deep, rectangle, trapezoid.
2 Side slope can be selected for cross section shapes: trapezoid, triangle.
Method Parameter
Bounded Initial Flow Rate or Initial Flow Rate per Area 1
Recession Recession Constant
Linear GW 1 and 2 Initial Flow Rate or Initial Flow Rate per Area
1
496
Chapter 16 Assessing Model Uncertainty
2 The available parameter depends on the method selected for specifying the recession threshold.
Table 35. Reach routing parameters available for sampling.
Method Parameter
Kinematic Wave Energy Slope
Manning's n
Diameter 1
Bottom Width 2
Side Slope 3
Lag Lag
Modified Puls Subreaches
Initial Flow
Muskingum K
X
Subreaches
Muskingum Cunge Energy Slope
Manning's n
Diameter 1
Bottom Width 2
Side Slope 3
Straddle Stager Lag
Duration
1 Available for the circle cross section.
2 Available for deep, rectangle, and trapezoid cross section.
3 Available for trapezoid and triangle cross section.
Method Parameter
Initial Condition Initial Elevation 1
Initial Outflow 1
Initial Storage 1
Orifice Outlet Coefficient
Culvert Entrance Coefficient
Exit Coefficient
Manning's n
Level Dam Top Coefficient
Nonlevel Dam Top Coefficient
497
Chapter 16 Assessing Model Uncertainty
Method Parameter
Broadcrested Coefficient
Spillway
Overtop Breach Top Elevation
Bottom Elevation
Bottom Width
Left Side Slope
Right Side Slope
Development Time
Piping Breach Top Elevation
Bottom Elevation
Piping Elevation
Bottom Width
Left Side Slope
Right Side Slope
Piping Coefficient
Development Time
1 The available parameter depends on the method selected for specifying the initial condition.
498
Chapter 16 Assessing Model Uncertainty
There are three different methods for sampling the selected parameter.
The Component Editor will display the appropriate property information
depending on the sampling method that you choose. The Simple
Distribution sampling method allows you to choose an analytical
probability distribution and then samples the parameter independently
of all other parameters selected in the uncertainty analysis. The
Monthly Distribution sampling method allows you to choose an
analytical probability distribution and enter unique distribution
properties for each month; the parameter is sampled independently of
all other parameters in the analysis. The Regression With
Additive Error sampling method allows you to choose a previously
selected parameter and then define a linear relationship between this
dependent parameter and the previously selected parameter; an
epsilon error term can be added to the linear relationship. The
information specified for each of the three sampling methods is
explained in the following sections.
Simple Distribution
499
Chapter 16 Assessing Model Uncertainty
Figure 379. Setting sampling properties for a parameter using the Simple
Distribution method.
500
Chapter 16 Assessing Model Uncertainty
Function Formula
Beta 1 ( )1 ( )1
() =
(, ) ( )+1
Alpha=p, Beta=q, Lower=min, Upper=max
B (p,q) is the Beta function
Exponential 1 ( )
() =
Mu=, Shift=min
Gamma ( )1
() =
()
501
Chapter 16 Assessing Model Uncertainty
Figure 380. Setting sampling properties for a parameter using the Monthly
Distribution method.
All of the available PDF choices include minimum and maximum limits
that are imposed on the parameter sampling process. The parameter
will be repeatedly sampled until a value is obtained that falls between
the minimum and maximum. These constraints are separate from the
mathematical features of the selected PDF. The limits are useful for
enforcing physical process limitations or restricting the sampling to a
reasonable range for the selected parameter. The same minimum and
maximum is used for all months of the year.
Regression With Additive Error
502
Chapter 16 Assessing Model Uncertainty
503
Chapter 16 Assessing Model Uncertainty
504
Chapter 16 Assessing Model Uncertainty
Figure 382. Preparing to choose the element and time-series output results
for an uncertainty analysis.
Figure 383. Choosing the elements and time-series for output results.
The Results window displays the elements and time-series that will be
stored when the uncertainty analysis is computed. Press the Select
button to choose additional time-series. The Select Results window
(Figure 384) is used to choose time-series at specific elements. All of
the elements in the uncertainty analysis are shown on the left side of
the window. Click on an element to choose a result from among the
time-series computed at the element. Only one element can be
selected at a time. After clicking on an element, the available time-
505
Chapter 16 Assessing Model Uncertainty
series at the element will be shown on the right side of the window.
Select the time-series that should be added to the stored results.
Multiple time-series can be selected simultaneously by holding the
control key and clicking on additional time-series. Press the Select
button and the selected time-series for the selected element will be
added to the stored results. Additional elements and time-series can
be selected. Press the Close button when you have finished choosing
results.
Figure 384. Choosing two time-series at an element for the output results.
The Results window is also used to stop storing elements and time-
series data when the uncertainty analysis is computed. Time-series
that have been removed will not be stored the next time the uncertainty
analysis is computed. Select the row in the table that should be
removed from the stored results. Multiple rows can be selected
simultaneously by holding the control key and clicking on additional
rows (Figure 385). Press the Remove button and the selected
elements and time-series will no longer be stored. Press the Close
button when you have finished removing elements and time-series.
506
Chapter 16 Assessing Model Uncertainty
507
Chapter 16 Assessing Model Uncertainty
508
Chapter 16 Assessing Model Uncertainty
509
Chapter 16 Assessing Model Uncertainty
produced the result could change while the result is open for viewing.
In this case, the open result will immediately show an indication that
data has changed and the uncertainty analysis needs to be
recomputed. After the analysis is recomputed the open results are
automatically updated with the new results and the indication is
updated with the date and time of the most-recent compute.
Sampled Parameters
The uncertainty analysis includes a number of parameters that are
sampled during the simulation. A new parameter value is sampled for
each sample of the simulation. The sequence of sampled values from
the first sample to the last sample is stored separately for each
parameter and may be viewed in a table. Each parameter sample can
be accessed by clicking the "Parameter" node on the "Results" tab of
the Watershed Explorer (Figure 387). The parameter sample table will
open. The table includes one row for each sample of the simulation.
The first column of the table displays the sample number while the
second column displays the sampled parameter value. An example of
a parameter sample is shown in Figure 388.
510
Chapter 16 Assessing Model Uncertainty
511
Chapter 16 Assessing Model Uncertainty
Figure 388. A parameter sample showing the sampled parameter value for
each sample of the uncertainty analysis simulation.
Individual Elements
The uncertainty analysis may include one or more locations where
output results have been selected. Special tabular and graphical
results may be viewed when outflow or pool elevation is selected. The
element outflow is available at all seven element types including
subbasin, reach, junction, source, diversion, reservoir, and sink. The
pool elevation is only available at reservoir elements where the storage
is specified with an elevation-storage or elevation-area curve. Other
computed time-series may be selected and will be stored in the
uncertainty analysis output HEC-DSS file. These other computed
time-series may be viewed and analyzed using the HEC-DSSVue
software.
512
Chapter 16 Assessing Model Uncertainty
513
Chapter 16 Assessing Model Uncertainty
Figure 390. A statistic output showing the outflow volume for each sample of
the uncertainty analysis simulation.
514
Chapter 16 Assessing Model Uncertainty
Figure 391. Outflow result graph for an uncertainty analysis showing the
mean and other statistical measures for each time interval.
A table and a graph are available when pool elevation is selected for a
reservoir location. The table displays the maximum pool elevation for
each sample of the uncertainty analysis. The first column of the table
displays the sample number while the second column shows the
maximum pool elevation. The table can be accessed by expanding the
reservoir node in the Watershed Explorer and clicking on the result
node in the same way described above for accessing outflow results.
An example of a maximum pool elevation table is shown in Figure 392.
515
Chapter 16 Assessing Model Uncertainty
516
Chapter 16 Assessing Model Uncertainty
Figure 393. Pool elevation result graph for an uncertainty analysis showing
the mean and other statistical measures for each time interval.
517
Chapter 16 Assessing Model Uncertainty
518
Chapter 17 Erosion and Sediment Transport
CHAPTER 17
519
Chapter 17 Erosion and Sediment Transport
Figure 394. Activating the erosion and sediment transport features in the
basin model.
520
Chapter 17 Erosion and Sediment Transport
same transport potential method will be used at all reaches within the
basin model. The available choices are shown in Table 38. A
cohesive transport potential method can also be selected. When
selected, transport of cohesive sediment is computed in addition to the
non-cohesive sediment.
Table 38. Transport potential functions for calculating the amount of
sediment that can be carried by the stream flow. Type is non-
cohesive (NC) or cohesive (CO). Method is excess shear (ES),
stream power (SP), or regression (RE).
521
Chapter 17 Erosion and Sediment Transport
reservoirs that require the calculation of a fall velocity will use the
method selected in the sediment properties. The available methods for
calculating fall velocity include Report 12 (Interagency Committee,
1957), Rubey (1933), Toffaleti (1968), and Van Rijn (1993). The
default selection is Rubey.
Two different grade scale systems are available. The selected system
will determine the number of grain sizes used for calculating erosion
and sediment transport throughout the basin model. Calculations of
erosion, deposition, settling, and transport are carried out separately
and independently for each size in the selected grade scale. The
Clay Silt Sand Gravel system only recognizes four broad size
classifications. The AGU 20 system recognizes 20 different size
classifications, generally breaking each of the four broad size
classifications into sub-classifications.
Subbasin
A subbasin element represents a catchment where precipitation falls
and causes surface runoff. Erosion within the catchment will result
from a combination of several different physical processes. Rain drops
cause erosion when they impact the ground surface and break apart
the top layer of the soil, dislodging soil particles to move with the
overland flow. Overland flow also imparts an erosive energy to the
ground surface that may further break apart the top soil layer. As the
overland flow rate increases, the flow becomes concentrated in rills
which focuses erosive energy and further erodes the surface. Total
erosion is closely linked to precipitation rate, land surface slope, and
condition of the surface. In some cases, soil eroded high in the
catchment may deposit before reaching the outlet of the subbasin.
Certain simulation features are common to all of the surface erosion
methods available for the subbasin element. Each erosion method
computes the total sediment load transported out of the subbasin
during a storm. This calculation process is repeated for each storm
during the simulation time window. The sediment load must be
distributed into a time-series of sediment discharge from the subbasin.
The distribution of sediment is based on the computed hydrograph and
the power function approach of Haan et al. (1994). The power function
is given by
=
where ct is the sediment concentration at time t, ki is the proportion of
the load for the current event i to the total annual load, Qt is the
subbasin discharge (flowrate) at time t, and a is an exponent entered
by the user.
522
Chapter 17 Erosion and Sediment Transport
Figure 396. Selecting the surface erosion method for a subbasin element.
523
Chapter 17 Erosion and Sediment Transport
Build-up Wash-off
The build-up wash-off erosion method is designed for urban
environments. In these environments, sediment accumulates in street
curbs due to wind deposit and erosion from pervious areas adjacent to
the curbs. The sediment accumulates during dry periods between
storms. During a storm, the accumulated sediment is flushed from the
street curbs by stormwater runoff. The method may optionally include
street sweeping operations designed to mechanically remove
accumulated sediment. A typical example of source data for
parameterizing the method is Breault et al. (2005) though data sources
should always be recognized as highly local. The Component Editor is
shown in Figure 397.
The initial time is an initial condition for the method. It specifies the
number of days since the last sweeping operation when the simulation
begins.
The half time specifies the number of days required for half of the
maximum solids to accumulate in the street curb under continuously
dry conditions.
524
Chapter 17 Erosion and Sediment Transport
Modified USLE
The modified USLE method (Williams, 1975) was adapted from the
original Universal Soil Loss Equation. The original equation was based
on precipitation intensity, and consequently could not differentiate
between storms with low or high infiltration. With high infiltration, there
is little surface runoff and little accompanying surface erosion.
Conversely, low infiltration events have relatively more surface runoff
and consequently more surface erosion. The modifications to the
original USLE equation changed the formulation to calculate erosion
from surface runoff instead of precipitation. The other components of
the original formulation remained the same. The method works best in
agricultural environments where it was developed. However, some
525
Chapter 17 Erosion and Sediment Transport
The erodibility factor describes the difficulty of eroding the soil. The
factor is a function of the soil texture, structure, organic matter content,
and permeability. Typical values range from 0.05 for unconsolidated
loamy sand to 0.75 for silty and clayey loam soils.
The topographic factor describes the suseptability to erosion due to
length and slope. It is based on the observation that long slopes have
more erosion than short slopes, and steep slopes have more erosion
than flat slopes. Typical values range from 0.1 for short and flat slopes
to 10 for long or steep slopes.
The cover factor describes the influence of plant cover on surface
erosion. Bare ground is the most suseptable to erosion while a thick
vegetation cover significantly reduces erosion. Typical values range
from 1.0 for bare ground, to 0.1 for fully mulched or covered soils, to as
small as 0.0001 for forest soils with a well developed soil O horizon
under a dense tree canopy.
The practice factor describes the effect of specific soil conservation
practices, sometimes called best management practices. Agricultural
practices could include strip cropping, terracing, or contouring.
Construction and urban practices could include silt fences, hydro
seeding, and settling basins. It is difficult to establish general ranges
for these practices as they are usually highly specific.
Only some precipitation events will cause surface erosion. The
threshold can be used to set the lower limit for runoff events that cause
526
Chapter 17 Erosion and Sediment Transport
erosion. Events with a peak flow less than the threshold will have no
erosion or sediment yield.
The exponent is used to distribute the sediment load into a time-series
sedigraph. A small value flattens the sedigraph compared to the
hydrograph. A large value heightens the sedigraph compared to the
hydrograph.
The gradation curve defines the distribution of the total sediment load
into grain size classes and subclasses. The gradation curve is defined
as a diameter-percentage function in the Paired Data Manager. The
current functions are shown in the selection list. If there are many
different functions available, you may wish to choose a function from
the selector accessed with the paired data button next to the selection
list. The selector displays the description for each function, making it
easier to select the correct one.
Reach
Sediment processes within a reach are directly linked to the capacity of
the stream flow to carry eroded soil. The transport capacity of the flow
can be calculated from the flow parameters and sediment properties.
If the stream can transport more sediment than is contained in the
inflow, additional sediment will be eroded from the stream bed and
entrained in the flow. However, if the flow in the reach cannot
transport the sediment of the inflow, entrained sediment will settle and
be deposited to the reach bed. The sediment transport capacity is
calculated using the transport potential method.
The deposition of sediment from the water column to the stream bed
requires time. The fall velocity of reach grain size class provides a
physical basis for determing how much time is required for sediment in
excess of the transport capacity to settle from the water column to the
stream bed. The settling velocity is calculated and multiplied by the
time interval to detmine the settling distance in one time interval. This
settling distance is then compared to the flow depth calculated during
flow routing in order to determine the fraction of calculated deposition
which is actually permitted during a time interval. The approach here
is similar to the one used in HEC-RAS.
The erosion of sediment from the stream bed and entraining into the
stream flow requires time. Erosion limitation has been observed but
the theoretical basis is less clear than for deposition limitation. An
empirical rule is used wherein a characteristic flow length necessary
for erosion is calculated. Erosion is limited when the length of the
reach divided by the flow depth is less than 30. This empirical rule is
similar to one used in HEC-RAS.
527
Chapter 17 Erosion and Sediment Transport
The reach bed is represented with a two-layer model. The upper layer
represents the top of the stream bed which actively interacts with the
flow field. This layer responds relatively quickly to changes in flow
rate. The lower layer is linked to the upper layer with a simple bed
mixing algorithm. The lower layer represents the underlying substrate
of the reach. Long-term processes that lead to a reach being either a
sediment sink or sediment source within the watershed are
represented through the lower layer. The two-layer model is also
capable of representing an "armouring" condition, which reduces the
rate of erosion as the upper layer coarsens.
As a simplifying assumption, the cross section shape of the channel is
fixed. There is no feedback between sediment dynamics and flow
routing. A characteristic width and depth of the sediment bed within
the reach is used to represent the two-layer model.
528
Chapter 17 Erosion and Sediment Transport
Fisher's Dispersion
The Fisher's dispersion method is based on an analysis of the
advection and diffusion of sediment within a reach (Fisher et al., 1979).
This is the most detailed of the sediment routing methods and requires
more data than some of the other available methods. Separate
specification of the advection and diffusion parameters for each grain
size class permits the large-grained sediments to move slower than the
fine-grained sediments. For each time interval, sediment from
upstream elements is added to the sediment already in the reach.
After erosion or deposition is calculated, the remaining available
sediment is translated in the reach by a travel time and attenuated
through a diffusion process. The advection and diffusion of the
sediment is linked to the velocity of water in the reach which is
calculated during the flow routing. The Component Editor is shown in
Figure 400.
529
Chapter 17 Erosion and Sediment Transport
The dispersion coefficient must be specified for each grain class (clay,
silt, sand, gravel). The dispersion coefficient indicates the diffusion of
the particles during transit through the reach and is closely connected
to channel geometry. The dispersion coefficient can vary over several
orders of magnitude and often must be adjusted during calibration.
Some guidance is available for estimating the dispersion coefficient, for
example, Kashefipour and Falconer (2002). The travel time must also
be specified for each grain class (clay, silt, sand, gravel) and is often
close to the travel time for water in the reach. When the AGU 20 grain
size classification is used, the same dispersion and retention values
are used for all subclasses of a grain class.
The initial gradation curve defines the distribution of the bed sediment
by grain size at the beginning of the simulation. The same curve is
used in both the upper and lower layers of the two-layer bed model.
The gradation curve is defined as a diameter-percentage function in
the Paired Data Manager. The current functions are shown in the
selection list. If there are many different functions available, you may
wish to choose a function from the selector accessed with the paired
data button next to the selection list. The selector displays the
description for each function, making it easier to select the correct one.
Using an erosion limit is optional. When the erosion limit is
deactivated, erosion is limited only by the transport capacity of the
flow. When the erosion limit is activated, actual erosion is reduced
when the ratio of reach length to flow depth is less than 30. The
erosion limit is usually encountered only in very short reaches.
Using a deposition limit is optional. When the deposition limit is
deactivated, sediment in excess of the transport capacity is deposited
completely. When the deposition limit is activated, sediment is limited
by the flow depth calculated during flow routing and the fall velocity of
each grain size. The fall velocity is computed using the method
selected with the basin model properties. Using a deposition limit
requires the specification of the water temperature in the reach. You
may specify a fixed temperature or select a temperature time-series
gage. Temperature gages must be created in the Time-Series Data
Manager before they can be used for the sediment method.
The width of the sediment bed must be specified. The width should be
typical of the reach and is used in computing the volume of the upper
and lower layers of the bed model. The depth of the bed must also be
specified. The depth should be typical of the total depth of the upper
and lower layers of the bed, representing the maximum depth of mixing
over very long time periods.
The active layer factor is used to calculate the depth of the upper layer
of the bed model. At each time interval, the upper layer depth is
530
Chapter 17 Erosion and Sediment Transport
Linear Reservoir
The linear reservoir method uses a simple linear reservoir to route
each grain size through the reach. For each time interval, available
sediment is calculated from the upstream sediment and local erosion
or deposition. The available sediment in each grain size class is
routed through a linear reservoir independently of the hydrologic
routing of the flow. This allows sediment of different grain sizes to
move at different speeds through the reach. Calibration with observed
data significantly improves the results. The Component Editor is
shown in Figure 401.
The retention parameter is equivalent to the storage coefficient in the
linear reservoir used to route the sediment through the reach. The
routing is performed separately for clay, silt, sand, and gravel. When
the AGU 20 grain size classification is used, the same retention is used
for all subclasses of each class. The value of the retention parameter
is analogous to the median length of time for each sediment class to
transit the reach. It may change for each sediment size class but is
often close to the travel time for water in the reach.
The initial gradation curve defines the distribution of the bed sediment
by grain size at the beginning of the simulation. The same curve is
used in both the upper and lower layers of the two-layer bed model.
531
Chapter 17 Erosion and Sediment Transport
Uniform Equilibrium
The uniform equilibrium method assumes that sediment is translated
instantaneously through the reach. It is the simplest method because
it does not compute any temporal lag for the sediment passing through
the reach. Sediment enters the reach from upstream elements. The
transport capacity for each grain size is calculated to determine the
deposition or erosion state. The available sediment is computed
subject to the limitations on deposition and erosion. The available
sediment is then passed out of the reach regardless of the velocity
determined during flow routing. The Component Editor is shown in
Figure 402.
532
Chapter 17 Erosion and Sediment Transport
The initial gradation curve defines the distribution of the bed sediment
by grain size at the beginning of the simulation. The same curve is
used in both the upper and lower layers of the two-layer bed model.
The gradation curve is defined as a diameter-percentage function in
the Paired Data Manager. The current functions are shown in the
selection list. If there are many different functions available, you may
wish to choose a function from the selector accessed with the paired
data button next to the selection list. The selector displays the
description for each function, making it easier to select the correct one.
Using an erosion limit is optional. When the erosion limit is
deactivated, erosion is limited only by the transport capacity of the
flow. When the erosion limit is activated, actual erosion is reduced
when the ratio of reach length to flow depth is less than 30. The
erosion limit is usually encountered only in very short reaches.
Using a deposition limit is optional. When the deposition limit is
deactivated, sediment in excess of the transport capacity is deposited
completely. When the deposition limit is activated, sediment is limited
by the flow depth calculated during flow routing and the fall velocity of
each grain size. The fall velocity is computed using the method
selected with the basin model properties. Using a deposition limit
requires the specification of the water temperature in the reach. You
may specify a fixed temperature or select a temperature time-series
gage. Temperature gages must be created in the Time-Series Data
Manager before they can be used for the sediment method.
The width of the sediment bed must be specified. The width should be
typical of the reach and is used in computing the volume of the upper
and lower layers of the bed model. The depth of the bed must also be
specified. The depth should be typical of the total depth of the upper
533
Chapter 17 Erosion and Sediment Transport
and lower layers of the bed, representing the maximum depth of mixing
over very long time periods.
The active layer factor is used to calculate the depth of the upper layer
of the bed model. At each time interval, the upper layer depth is
computed as the d90 of the sediment in the upper layer, multiplied by
the active layer factor.
Volume Ratio
The volume ratio method links the transport of sediment to the
transport of flow in the reach using a conceptual approach. For each
time interval, sediment from upstream elements is added to the
sediment already in the reach. The deposition or erosion of sediment
is calculated for each grain size to determine the available sediment for
routing. The proportion of available sediment that leaves the reach in
each time interval is assumed equal to the proportion of stream flow
that leaves the reach during that same interval. This means that the all
grain sizes are transported through the reach at the same rate, even
though erosion and deposition are determined separately for each
grain size. This assumption essentially limits the advection velocity of
the sediment to the bulk water velocity. Therefore, this method is a
higher fidelity option than the Equilibrium method (without requiring
much more data) but less precise than the other available methods
which require substantially more data. This method is very similar to
the approach used in the Soil and Water Assessment Tool model
(Gassman, 2007). The Component Editor is shown in Figure 403.
The initial gradation curve defines the distribution of the bed sediment
by grain size at the beginning of the simulation. The same curve is
used in both the upper and lower layers of the two-layer bed model.
The gradation curve is defined as a diameter-percentage function in
the Paired Data Manager. The current functions are shown in the
534
Chapter 17 Erosion and Sediment Transport
selection list. If there are many different functions available, you may
wish to choose a function from the selector accessed with the paired
data button next to the selection list. The selector displays the
description for each function, making it easier to select the correct one.
Using an erosion limit is optional. When the erosion limit is
deactivated, erosion is limited only by the transport capacity of the
flow. When the erosion limit is activated, actual erosion is reduced
when the ratio of reach length to flow depth is less than 30. The
erosion limit is usually encountered only in very short reaches.
Using a deposition limit is optional. When the deposition limit is
deactivated, sediment in excess of the transport capacity is deposited
completely. When the deposition limit is activated, sediment is limited
by the flow depth calculated during flow routing and the fall velocity of
each grain size. The fall velocity is computed using the method
selected with the basin model properties. Using a deposition limit
requires the specification of the water temperature in the reach. You
may specify a fixed temperature or select a temperature time-series
gage. Temperature gages must be created in the Time-Series Data
Manager before they can be used for the sediment method.
The width of the sediment bed must be specified. The width should be
typical of the reach and is used in computing the volume of the upper
and lower layers of the bed model. The depth of the bed must also be
specified. The depth should be typical of the total depth of the upper
and lower layers of the bed, representing the maximum depth of mixing
over very long time periods.
The active layer factor is used to calculate the depth of the upper layer
of the bed model. At each time interval, the upper layer depth is
computed as the d90 of the sediment in the upper layer, multiplied by
the active layer factor.
Reservoir
A reservoir provides an ideal environment for removing sediment from
water. Quiescent water, or completely still water, has no capacity to
suspend sediment because the horizontal velocity of the water is zero.
Therefore, given enough time all sediment in the water column will
settle to the bottom of the reservoir. Time is important because
sediment does not settle instantaneously. There is a very small
amount of friction between each sediment grain and the water
molecules that surround it. This friction creates a turbulent
environment that slows the descent of each sediment grain. The
velocity at which a particular grain will settle is due to its size and to a
lesser extent, the density and shape of the grain. Large grains will
settle very quickly while small grains will settle slower. The very fine-
535
Chapter 17 Erosion and Sediment Transport
grained sediments such as silt and clay may require years to settle
from the water column. The greater the residence time of the
reservoir, the higher the percentage of the inflowing sediment will settle
to the bottom. However, the very smallest clay particles may never
settle due to electrical charges in the minerials composing the particles
and the charges in water molecules caused by hydrogen bonding.
Sediment entering a reservoir and settling to the bottom eliminates
some of the storage space of the reservoir. Over time the cumulative
storage lost to sediment settling may be significant. The program will
compute the sediment balance in the reservoir over time. However,
the volume of sediment accumulating in the reservoir is not linked to
the storage characteristics of the reservoir used for flow routing.
536
Chapter 17 Erosion and Sediment Transport
Specified Sediment
The specified sediment method uses observations of sediment at the
outlet of the reservoir to compute the sediment leaving the reservoir at
each time interval. The sediment observed at the outlet may be more
or less than the incoming sediment during a time interval. Settling is
presumed to occur when the observed sediment is less than the
incoming sediment. When the observed sediment is more than the
537
Chapter 17 Erosion and Sediment Transport
538
Chapter 17 Erosion and Sediment Transport
The only parameter for the Chen sediment trap method is the depth-
averaged temperature of the water in the reservoir. You may specify a
fixed temperature or select a temperature time-series gage.
Temperature gages must be created in the Time-Series Data Manager
before they can be used for the sediment method.
Source
A source element provides a way to introduce inflow to the channel
network. Usually the inflow represents a portion of the upstream
watershed that is not explicitly modeled in the basin model.
Sometimes the inflow is computed in a separate simulation model,
where the outflow hydrograph from that model becomes the inflow at
the source element. Other times the inflow is observed at a flow
measurement station and becomes the inflow at the source element. It
is critical in a basin model that contains a source element to introduce
both the water and sediment inflow at the source.
539
Chapter 17 Erosion and Sediment Transport
Annual Load
The annual load method applies an average annual sediment load to
the discharge hydrograph with a methodology that may be suitable for
continuous simulation applications where the land use is not changing
in time. The total load must be distributed into a time-series in a
manner that is consistent with the outflow hydrograph at the source.
The method adopted in the program is to use an exponential form
developed by Haan et al. (1994). The exponential equation is given by
=
where ct is the sediment concentration at time t, k is the ratio of the
total sediment load for the time window to the total discharge for the
540
Chapter 17 Erosion and Sediment Transport
The total load specifies the accumulated sediment over all grain size
classes for a typical year. The total load depends on the land use of
the area upstream of the source element and can be highly variable
due to local conditions. The load is fixed for the simulation time
window.
The exponent is used to distribute the sediment load into a time-series
sedigraph. A small value flattens the sedigraph compared to the
hydrograph. A large value heightens the sedigraph compared to the
hydrograph.
The gradation curve defines the distribution of the total sediment load
into grain size classes and subclasses. The gradation curve is defined
as a diameter-percentage function in the Paired Data Manager. The
current functions are shown in the selection list. If there are many
different functions available, you may wish to choose a function from
the selector accessed with the paired data button next to the selection
list. The selector displays the description for each function, making it
easier to select the correct one.
Specified Load
The specified sediment method uses observed sediment load time-
series data to determine the sediment to add to the flow at the source
element. Observations of sediment are less common than
observations of stream flow. In cases where observations are not
available at the location of the source element, this method could be
used for sediment computed from an external model of sediment
discharge. The Component Editor is shown in Figure 409.
541
Chapter 17 Erosion and Sediment Transport
Junction
A junction element is used to combine the outflow from multiple
upstream elements in order to produce one outflow. The outflow from
the junction is usually connected downstream to a reach or other
elements. All water that enters a junction must leave the junction
because no water is stored at the junction. Therefore, all of the inflows
to the junction are added together in order to compute the total inflow
hydrograph to the junction. Similarily, all sediment in each of the
inflows is added to compute the total sedigraph. The total inflow
hydrograph and total inflow sedigraph are connected to the
downstream element.
Diversion
A diversion element is used to represent a hydraulic control structure
that directs some of the stream flow out of the channel and into a
storage pond or side channel. The diversion method calculates the
amount of flow to be diverted and subtracts it from the inflow to the
element. The diverted flow may be connected to a downstream
542
Chapter 17 Erosion and Sediment Transport
543
Chapter 17 Erosion and Sediment Transport
Passage Efficiency
Diversion structures are usually designed primarily to control the
movement of water, redirecting some portion of the inflow away from
the main channel. In general the d50 of the sediment in the diverted
flow is smaller than the d50 of the sediment in the main channel. The
passage efficiency method provides a simple approach to controlling
the differential sediment movement that typically occurs at diversion
structures. The Component Editor is shown in Figure 411.
The clay percentage in the diverted flow must be specified. The
incoming clay load to the diversion element is converted to a
concentration. Subsequently the clay concentration in the diverted
flow is calculated as the incoming concentration multiplied by the
specified percentage. A mass balance is then performed in order to
calculate the load in the diverted flow and the load remaining in the
channel. When the AGU 20 grain size classification is used, the same
percentage is used for all subclasses of clay.
The silt percentage in the diverted flow must be specified. The silt
percentage may be close to the clay percentage if the velocity through
the diversion structure is high. When the AGU 20 grain size
classification is used, the same percentage is used for all subclasses
of silt.
The sand percentage in the diverted flow must be specified. Unless
the diversion structure is specifically designed to direct sediment from
the main channel, the sand percentage will be much smaller than the
percentage for clay and silt. When the AGU 20 grain size classification
is used, the same percentage is used for all subclasses of sand.
The gravel percentage in the diverted flow must be specified. The
percentage for gravel is often close to zero. When the AGU 20 grain
544
Chapter 17 Erosion and Sediment Transport
Specified Load
The specified sediment method uses observed sediment load time-
series data to determine the sediment to remove from the inflow to the
diversion element and transfer with the diverted flow. Remaining
sediment remains in the channel and continues to be routed
downstream of the diversion. If observations of sediment are not
available where the diverted flow is discharged out of the channel, then
an external model could be used to comput the sediment load. The
Component Editor is shown in Figure 412.
The total sediment discharge for the diverted flow is specified with a
time-series of sediment load. The time-series should represent the
total of all grain size classes. The time-series data must be defined in
the Time-Series Data Manager before it can be selected in the
sediment editor. The time-series button next to the selection list can
be used to assist in selecting the correct sediment gage.
The total sediment discharge in the diverted flow is apportioned into
grain size classes using a grain size distribution curve. The grain size
distribution curve is specified with a diameter-percentage function.
The function data must be defined in the Paired Data Manager before
it can be selected in the sediment editor. The paired data button next
to the selection list can be used to assist in selecting the correct
function.
Sink
A sink element is used to mark the outlet of a watershed model. No
water leaves the sink. Therefore, all of the inflows to the sink are
added together in order to compute the total inflow hydrograph to the
sink. Just as water cannot leave the sink, sediment also cannot leave
545
Chapter 17 Erosion and Sediment Transport
the sink. Similar to the water, all sediment in each of the inflows is
added to compute the total sedigraph.
546
Chapter 17 Erosion and Sediment Transport
Subbasin
The erosion results for each subbasin include the time-series of total
sediment load leaving the subbasin. The total sediment is the sum of
all grain size classes. Also included is the time-series of sediment load
in each grain size class. The total sediment concentration in the
outflow is also included among the results. The organization of the
results is shown in Figure 413.
547
Chapter 17 Erosion and Sediment Transport
Reach
The sediment results for each reach include the time-series of total
sediment load entering the reach, and the total sediment load leaving
the reach. The total sediment is the sum of all grain size classes. Also
included is the time-series of sediment load in each grain size class
both coming into the reach and leaving the reach. The organization of
the results is shown in Figure 414.
548
Chapter 17 Erosion and Sediment Transport
Reservoir
The sediment results for each reservoir include the time-series of total
sediment load entering the reservoir, and the total sediment load
leaving the reservoir. Additionally, the time-series of total sediment
load settled to the bottom of the reservoir is available, along with the
time-series of sediment load in suspension in the reservoir pool. For
each of the four categories, the total sediment is the sum of all grain
size classes. Also included is the time-series of sediment load in each
grain size class for each of the four categories of sediment at a
reservoir. The organization of the results is shown in Figure 415.
549
Chapter 17 Erosion and Sediment Transport
Source
The sediment results for each source include the time-series of total
sediment load leaving the source. The total sediment is the sum of all
grain size classes. Also included is the time-series of sediment load
within each grain size class.
Junction
The sediment results for each junction include the time-series of total
sediment load leaving the junction. The total sediment is the sum of all
sediment from each of the upstream elements. Also included is the
time-series of sediment load for each grain size class, computed
similarly to the total load but specific to each grain size.
Diversion
The sediment results for each diversion include the time-series of total
sediment load entering the diversion, and the total sediment load
leaving the diversion. The total sediment is the sum of all grain size
classes. Additionally, the total sediment load in the diverted flow is
available. Also included is the time-series of sediment load in each
grain size class for each of the three categories of sediment at a
diversion.
Sink
The sediment results for each sink include the time-series of total
sediment load entering the sink. The total sediment is the sum of all
sediment from each of the upstream elements. Also included is the
time-series of sediment load for each grain size class, computed
similarly to the total load but specific to each grain size.
References
Ackers, P. and W.R. White. 1973. "Sediment Transport: New
Approach and Analysis." Journal of the Hydraulics Division, vol 99, no
HY11, pp 2040-2060.
Breault, R.F., K.P. Smith, and J.R. Sorenson. 2005. Residential
Street-dirt Accumulation Rates and Chemical Composition, and
Removal Efficiencies by Mechanical- and Vacuum-type Sweepers,
New Bedford, Massachusetts, 2003-04. U.S. Geological Survey,
Scientific Investigations Report 2005-5184.
Chen, C.N. 1975. "Design of Sediment Retention Basins."
Proceedings of the National Symposium on Urban Hydrology and
Sediment Control, pp 285-298.
550
Chapter 17 Erosion and Sediment Transport
551
Chapter 17 Erosion and Sediment Transport
552
Chapter 18 Water Quality
CHAPTER 18
Water Quality
Water flowing over the land surface comes in contact with soil and
organic matter. Many different chemical compounds are present in the
soil and may dissolve upon contact with the water and subsequently be
carried along with the water to the catchment outlet. The organic
matter is constantly decomposing, and in the process, various
chemical compounds are released and may also be entrained in the
surface flow. Catchment outflows enter the channel and move
downstream. As water moves down the stream channel, the chemical
compounds carried along in the water undergo transformations
mediated by algae and bacteria. Nutrients of nitrogen and phosphorus
are particularly noteworthy among the various chemical compounds.
The algae and bacteria communities respond to the nutrients and key
water quality criteria such as dissolved oxygen are impacted. This and
other water quality criteria are closely linked to the presence and
movement of water and are best evaluated in conjunction with the
watershed hydrology. At this time, water quality evaluations are limited
to the nutrients nitrogen and phosphorus.
553
Chapter 18 Water Quality
554
Chapter 18 Water Quality
555
Chapter 18 Water Quality
Figure 417. Activating the water quality features in the basin model.
Figure 418. Specifying dissolved oxygen properties that will be applied to all
elements within the basin model.
556
Chapter 18 Water Quality
557
Chapter 18 Water Quality
Table 40. Algae and dissolved oxygen parameters used in water quality
modeling.
558
Chapter 18 Water Quality
559
Chapter 18 Water Quality
Many of the parameters used for water quality modeling are sensitive
to temperature. The reactions are assumed to happen at the water
temperature and provision is made for specifying the temperature for
each element. Temperature dependence requires the parameters be
adjusted as the temperature changes. The formula is used to adjust
the parameter value as a function of temperature is given as:
() = (20) 20
where k is a parameter that depends on temperature, T is the current
water temperature, k(20) is the standard parameter value at 20C, and
is a temperature correction coefficient specific to each parameter. If
a parameter shows temperature dependence, the temperature
correction coefficient is given immediately following the parameter in
the Component Editor.
Subbasin
The water quality in streams begins not in the streams, but on the land
surface adjacent to the streams. This is the portion of the watershed
that is represented by subbasin elements. Many physical processes
are continusouly happening on the land surface and determine the
quality of the water leaving the subbasins and entering the stream
channel. These processes occur naturally but may be affected by land
use alterations.
Many watersheds contain plants growing on the landscape. Natural
plant communities include forests and grasslands. Some areas may
have been converted to agricultural production. In both cases, plants
560
Chapter 18 Water Quality
extract nutrients (nitrogen and phosphorus) from the soil and convert
the nutrients to live organic matter. Agricultural production may include
the application of natural or synthetic fertilizers that contain nitrogen
and phosphorus components, increasing the amount of nutrients
above natural levels. When the plants die, the organic matter remains
on the soil surface where the decay process begins. Products of the
decay process include various forms of nitrogen and phosphorus.
These products may be entrained in infiltrating precipitation
replenishing the nutrient stores in the soil. Decay products may also
be entrained or dissolved in overland flow that leaves the subbasin as
outflow. Nutrients entrained in outflow from land surface areas are the
primary source of nutrients in the stream channel.
When performing water quality calculations, the outflow from the
subbasin must also include constituent concentrations. Several
methods are available for determining the concentrations to be linked
with the outflow hydrograph.
561
Chapter 18 Water Quality
Constant Concentration
The constant concentration method may be useful for simulations
where steady-state conditions prevail most of the time. It may also be
useful for preliminary model configuration before moving to a more
advanced method. The Component Editor is shown in Figure 420.
562
Chapter 18 Water Quality
Specified Concentration
The specified concentration method is designed for situations where
the concentration in the outflow from a catchment is measured at a
monitoring station. Alternately, it can be used if the concentration is
calculated externally from the program. In either case, the user may
enter a time-series of concentration for a particular constituent. The
Component Editor is shown in Figure 421.
563
Chapter 18 Water Quality
Reach
The reach element is used for representing streams and rivers. The
flow routing performed in the reach forms the basis of the water quality
simulation. Results from the flow routing are used as the hydronamic
field for the advection-diffusion routing of the water quality constituents.
The advection-diffusion algorithm used for routing water quality
constituents is based on the QUICKEST scheme described by Leonard
(1971). It is an explicit third-order accurate method, where the second-
order approximation polynomial is centered on the upstream control
volume. It was designed for solving problems characterized by highly
advective unsteady flows, such as can be routinely found in moderate
to steep sloped streams. Because of the third-order accuracy, it is
possible for unrealistic undershoots and overshoots to appear in the
solution, especially when steep concentration fronts are simulated.
These artifacts are significantly reduced by employing the ULTIMATE
scheme as part of the solution algorithm (Leonard, 1991). The
combined QUICKEST-ULTIMATE solution scheme performs very well
564
Chapter 18 Water Quality
565
Chapter 18 Water Quality
566
Chapter 18 Water Quality
Figure 423. The Nutrient Simulation Module water quality method for a reach
element.
The initial water quality conditions in the reach can be set using one of
two available methods. The Inflow Concentrations method will
initialize each water quality node in the reach using the inflows to the
reach element. All inflows to the reach are combined and a mass
balance is perfomed to calculate the concentration for each state
variable. The concentration at the beginning of the simulation time
window is used as the initial concentration at each node in the reach.
The Specified Values method allows you to set the concentration
of each state variable independent of the concentrations in the inflow.
You must enter a value for each state variable and the same value will
be used as the initial concentration at each node in the reach.
The air temperature must be specified for the reach. It may be
specified as an average value for the entire simulation duration. The
same average value is used for each time interval of the simulation.
Choose the Average air temperature method and enter a value.
Alternatively, the temperature may be specified using a time-series
gage. The gage must be defined in the time-series manager before it
can be used in the reach element. Choose the Time-Series air
temperature method and then select the correct gage from the list of
choices for temperature gages.
The windspeed must be specified for the reach. It may be specified as
an average value for the entire simulation duration. The same average
value is used for each time interval of the simulation. Choose the
Average windspeed method and enter a value. Alternatively, the
windspeed may be specified using a time-series gage. The gage must
be defined in the time-series manager before it can be used in the
reach element. Choose the Time-Series windspeed method and
then select the correct gage from the list of choices for windspeed
gages.
567
Chapter 18 Water Quality
= 0.0112 2 *
where K is the dispersion coefficient, is the cross section mean
velocity, W is the channel width, and du* is the derivative of the shear
velocity. The minimum and maximum limits are applied after the
calculation.
Reservoir
The reservoir element is used for representing lakes and reservoirs.
The storage and outflow routing performed in the reservoir forms the
basis of the water quality simulation. Results from the flow routing are
used as the starting point for the simulation of the nutrients.
Nutrients in the reservoir water will undergo biological and chemical
transformations while in storage. A variety of process kinetics
representations have been proposed for use in nutrient simulation and
can be found in the published literature. Currently only one approach
is included within the program for reservoirs: the Nutrient Simulation
Module. The Nutrient Simulation Module (Zhang and Johnson, 2012)
provides for simulation of nutrients and dissolved oxygen in the water
column.
Nutrients are suspended in the water that enters the reservoir. The
water within the reservoir is assumed to be well-mixed and the nutrient
concentrations are therefore assumed to be uniform throughout the
water in storage. For a given time during the simulation, the nutrients
in the reservoir inflow are added to the nutrient concentrations already
in reservoir storage. The constituent concentrations in the storage
water are updated to reflect the added nutrient mass of the inflow. The
Nutrient Simulation Module is used to compute the change in each
nutrient concentration according to biological and physical processes.
The concentrations at the end of the process calculations are applied
to the outflow determined during the flow routing. The mass of
nutrients leaving the reservoir are then calculated based on a volume
568
Chapter 18 Water Quality
Figure 424. Selecting the water quality method for a reservoir element.
569
Chapter 18 Water Quality
570
Chapter 18 Water Quality
Figure 425. The Nutrient Simulation Module water quality method for a
reservoir element.
Source
The source element is frequently used to represent a portion of the
watershed that will not be explicitly represented in the basin model.
The stream flow from the unmodeled portion of the watershed is
introduced to the basin model at the point of the source element. The
inflow at the source element may be from a stream gaging station or
an external simulation model. The hydrograph at the source element
is added into the channel network and joins the stream flow computed
within the basin model.
571
Chapter 18 Water Quality
Constant Concentration
The constant concentration method may be useful for simulations
where steady-state conditions prevail most of the time. It may also be
useful for preliminary model configuration before moving to a more
advanced method. The Component Editor is shown in Figure 427.
572
Chapter 18 Water Quality
Specified Concentration
The specified concentration method is designed for situations where
the concentration in the inflow from a source is measured at a
monitoring station. Alternately, it can be used if the concentration is
calculated externally from the program. In either case, the user may
enter a time-series of concentration for a particular constituent. The
Component Editor is shown in Figure 428.
573
Chapter 18 Water Quality
Junction
A junction element takes inflow from multiple upstream elements to
create a combined outflow. The upstream elements could be two
reach elements that join at a confluence. Alternately, the upstream
elements could be a reach element and a subbasin element that
represent a significant lateral inflow to the stream. The purpose of the
junction is to provide conservation of mass when the outflows from
multiple elements come together and continue downstream. The
conservation of mass is performed assuming there is no storage at the
junction. Therefore, all inflow to the junction during a time interval
leaves the junction immediately as outflow.
When performing water quality calculations, the junction also performs
a mass balance for each constituent. The concentration of each
constitute in the inflow is multiplied by the inflow rate to determine the
mass of each constituent entering the junction. The mass calculated in
this manner for each upstream element is totaled across all
constituents. Then the total mass of each constituent is converted
574
Chapter 18 Water Quality
back to a concentration leaving the junction using the outflow rate from
the junction.
Diversion
A diversion element begins calculations with a combined inflow
hydrograph from upstream elements. Based on the selected diversion
method, some portion of the inflow is removed from the main channel
and sent to the diversion flow. The remaining flow becomes the
outflow from the diversion and continues down the main channel. The
conservation of mass that considers the inflow, diversion flow, and
outflow, is performed assuming there is no storage at the diversion.
Therefore, all inflow to the diversion during a time interval leaves the
diversion either as diversion flow or outflow.
When performing water quality calculations, the diversion maintains
the concentration of constituents through a mass balance. The
concentration of each constituent in the combined inflow is known.
The same concentration of each constituent is passed to the diversion
flow, assuming the diversion structure is not able to affect the
concentration of dissolved constituents. Likewise, the same
concentration of each constituent in the inflow is passed to the outflow.
Sink
A sink element takes inflow from multiple upstream elements to create
a combined flow that is removed from the channel network. For each
time interval, the inflow from each upstream element is accumulated to
calculate the total inflow for the interval. The total inflow is presented
as the flow rate at the sink.
When performing water quality calculations, the sink performs a mass
balance for each constituent. The concentration of each constitute in
each upstream inflow is multiplied by the inflow rate to determine the
mass of each constituent. The mass calculated in this manner is
totaled across all upstream elements for each constituent. Then the
total mass of each constituent is converted back to a concentration
with respect to the total flow at the sink.
575
Chapter 18 Water Quality
Subbasin
The water quality results for a subbasin are the time-series of
concentration for each of the constituents. All eight of the constituents
are included. The results represent the quality of the water leaving the
subbasin and entering the downstream channel network. The
organization of the results is very similar to results for the junction
element.
Reach
The water quality results for a reach include the time-series of
concentration for each constituent, both entering the reach and leaving
the reach. All eight constituents in the nutrient water quality simulation
are available at each location of the reach. The organization of the
results is shown in Figure 429.
576
Chapter 18 Water Quality
Reservoir
The water quality results for a reservoir include the time-series of
concentration for each constituent, both entering the reservoir and
leaving the reservoir. The concentration of each constituent in the
reservoir pool is also available. All eight constituents in the nutrient
water quality simulation are available at each location of the reservoir.
The organization of the results is shown in Figure 430.
577
Chapter 18 Water Quality
Source
The water quality results for a source are the time-series of
concentration for each of the constituents. All eight of the constituents
are included. The results represent the quality of the water leaving the
source and entering the downstream channel network. The
organization of the results is very similar to results for the junction
element.
Junction
The water quality results for a junction are the time-series of
concentration for each of the constituents. All eight of the constituents
are included, and represent a mass balance of the upstream inflows
with the constituent concentrations in each inflow. The results
represent the quality of the water leaving the junction and continuing
578
Chapter 18 Water Quality
Diversion
The water quality results for a diversion include the time-series of
concentration for each constituent entering the diversion. Separate
results are also available for the concentration of each constituent in
the main outflow that continues downstream in the channel network,
and in the diverted flow. All eight constituents in the nutrient water
quality simulation are available at each location of the diversion. The
organization of the results is very similar to results for the reservoir
element.
Sink
The water quality results for a sink are the time-series of concentration
for each of the constituents. All eight of the constituents are included,
and represent a mass balance of the upstream inflows with the
constituent concentrations in each inflow. The results represent the
579
Chapter 18 Water Quality
References
Cahyono, M. 1993. Three-Dimensional Numerical Modelling of
Sediment Transport Processes in Non-Stratified Estuarine and Coastal
Waters. Ph.D. Thesis, Dept of Civil Engineering, University of
Bradford, Bradford, England.
Fisher, H.B., E.J. List, R.C.Y. Koh, J. Imberger, and N.H. Brooks.
1979. Mixing in Inland and Coastal Waters. Academic Press, New
York.
Leonard, B. 1971. "A Stable and Accurate Convective Modeling
Procedure Based on Quadratic Upstream Interpolation." Computer
Methods in Applied Mechanics and Engineering, vol 19, pp 59-98.
Leonard, B. 1991. "The ULTIMATE Conservative Difference Scheme
Applied to Unsteady One-Dimensional Advection." Computer Methods
in Applied Mechanics and Engineering, vol 88, pp71-74.
Zhang, Z. and B. Johnson. 2012. DRiverine Nutrient Simulation
Modules Developed for Hydrologic and Hydraulic Models. U.S. Army
Corps of Engineers, Environmental Laboratory, Vicksburg, MS.
Technical Report ERDC/EL TR-12-X.
580
Appendix A Data Storage in HEC-DSS
APPENDIX A
Descriptors
AREA Computed area for a reservoir element.
ATI-COLD Computed antecedent temperature index for cold content
in the snowmelt method.
ATI-COLDRATE Manually entered antecedent temperature index
versus coldrate function for a paired data function.
ATI-MELT Computed antecedent temperature index for meltrate in
snowmelt.
ATI-MELTRATE Manually entered antecedent temperature index
versus meltrate function for a paired data function.
COLD CONTENT Computed cold content in snowmelt.
CROP COEFFICIENT Manually entered time-series of crop
coefficient for a time-series gage.
DISTANCE-ELEVATION Manually entered cross section in the paired
data manager.
ELEVATION Computed pool elevation for a reservoir element.
ELEVATION-OBSERVED Observed pool elevation for a reservoir
element using the optional observed elevation time-series gage.
ELEVATION-RESIDUAL Residual elevation for a reservoir element
with observed elevation. The residual is calculated as computed flow
minus the observed elevation.
ELEVATION-AREA A manually entered elevation-area function
defined in the paired data manager.
581
Appendix A Data Storage in HEC-DSS
582
Appendix A Data Storage in HEC-DSS
583
Appendix A Data Storage in HEC-DSS
584
Appendix A Data Storage in HEC-DSS
585
Appendix A Data Storage in HEC-DSS
586
Appendix B Grid Cell File Format
APPENDIX B
File Definition
The grid-cell file begins with the keyword "Parameter Order" followed
by a colon and parameter keywords indicating the order for reading
parameters from the file (Figure 432). The keyword "End" must be on
a line by itself after the "Parameter Order" line. Valid parameter
keywords are shown in Table 42. Parameter keywords are not case
sensitive and are separated by spaces. If the parameter order is not
defined, it is assumed to be: Xcoord Ycoord TravelLength Area. The
coordinate system of Xcoord and Ycoord used in the file must match
the coordinate system used in the gridded DSS records. Typically the
coordinate system will be either hydrologic rainfall analysis project
(HRAP) or standard hydrologic grid (SHG).
Table 42. Parameter keyword definitions.
587
Appendix B Grid Cell File Format
The data for a subbasin begins with the keyword "Subbasin" followed
by a colon and the subbasin identifier. One line beginning with the
keyword "Grid Cell" follows for each cell in the subbasin. Data for the
subbasin ends with the keyword "End". Keywords are not case
sensitive and may contain spaces. Blank lines can be included and
lines beginning with "#" are ignored as comments. The same grid-cell
file can be referenced by more than one subbasin, allowing data for
many subbasins to be stored in the same file. The identifier for a
subbasin must be exactly the same in the grid-cell file as it is in the
basin model.
588
Appendix B Grid Cell File Format
589
Appendix B Grid Cell File Format
590
Appendix C Map File Format
APPENDIX C
File Definition
Watershed boundary and stream line features are both defined in the
same file, which is in plain ASCII format. Each feature type is
contained in a separate section of the file; it is not important which
section is first in the file. Each section begins with the keyword
"MapGeo" followed by a colon and either "BoundaryMap" or
"RiverMap" (Figure 433).
A map segment defines a list of map coordinates that are connected
by a line. A closed segment defines a polygon and an open segment
defines a line. Closed segments are used for watershed boundaries
and open segments are used for stream lines. Each segment begins
with the keyword "MapSegment" followed by a colon and either
"Closed" or "Open." The last coordinate in a closed segment is
automatically connected to the first coordinate.
Segment coordinates are defined with x-y pairs. Map features are
automatically scaled in the Basin Model screen. Coordinates are
therefore independent of projection, units, and offset. All segments
must be in the same coordinate system.
591
Appendix C Map File Format
MapGeo: BoundaryMap
MapSegment: closed
582242.875000, 4174922.500000
582220.875000, 4174961.500000
582205.625000, 4175013.750000
581981.000000, 4174672.750000
582025.812500, 4174696.250000
582068.812500, 4174711.000000
MapSegment: closed
582810.125000, 4174024.500000
582874.687500, 4173973.750000
582950.687500, 4173902.750000
582554.000000, 4174000.250000
582667.687500, 4174003.750000
582810.125000, 4174024.500000
MapGeo: RiverMap
MapSegment: open
582750.187500, 4176706.000000
582687.000000, 4176594.000000
582657.375000, 4176468.500000
582613.125000, 4176359.500000
582482.125000, 4174521.500000
582555.250000, 4174377.500000
582555.250000, 4174378.000000
MapSegment: open
582941.500000, 4175098.500000
582920.500000, 4175009.750000
582912.312500, 4174956.500000
582699.375000, 4174540.500000
582618.250000, 4174468.250000
Figure 433. Sample background map file.
592
Appendix D HEC-HMS and HEC-1 Differences
APPENDIX D
Recession Baseflow
The recession baseflow method includes a recession constant for
specifying the rate at which recession flow decreases with time. In
HEC-HMS the parameter is defined as the ratio of the current
recession flow to the recession flow one day earlier. The HEC-1
parameter, called RTIOR, is defined as the ratio of the current
recession flow to the flow one hour later. The following equation can
be used to convert an HEC-1 recession constant for use in HEC-HMS:
1
Recession Constant =
( RTIOR )24
593
Appendix D HEC-HMS and HEC-1 Differences
c
=
( ) (
(10 w z y ) + 16 z y 2 z 2 + 1 + (5 w 2 ) + 6 w y z 2 + 1 )
v (
3 (w + 2 z y ) w + 2 y z 2 + 1 )
where w is the bottom width, z is the side slope, and y is the flow
depth. Velocity is computed using Manning's formula and the
properties of the cross section.
In HEC-1 the cross-sectional flow area is computed as:
1
Q m
A =
594
Appendix D HEC-HMS and HEC-1 Differences
595
Appendix D HEC-HMS and HEC-1 Differences
596
Appendix E Terms and Conditions of Use
APPENDIX E
597
Appendix E Terms and Conditions of Use
598