HEC-HMS Users Manual 3.5
HEC-HMS Users Manual 3.5
HEC-HMS Users Manual 3.5
User's Manual
Version 3.5
August 2010
CPD-74A
Form Approved
OMB No. 0704-0188
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 Washington Headquarters Services, Directorate for Information Operations and Reports, 1215 Jefferson Davis Highway, Suite 1204, Arlington,
VA 22202-4302, and to the Office of Management and Budget, Paperwork Reduction Project (0704-0188), Washington DC 20503.
2. REPORT DATE
August 2010
8. PERFORMING ORGANIZATION
REPORT NUMBER
Distribution is unlimited.
The Hydrologic Modeling System (HEC-HMS) is designed to simulate the precipitation-runoff processes of dendritic
watershed systems. It supersedes HEC-1 and provides a similar variety of options but represents a significant advancement in
terms of both computer science and hydrologic engineering. In addition to unit hydrograph and hydrologic routing options,
capabilities include a linear quasi-distributed runoff transform (ModClark) for use with gridded precipitation, continuous
simulation with either a one-layer or more complex five-layer soil moisture method, and a versatile parameter estimation
option.
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 user 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.
Program functionality and appearance are the same across all supported platforms. It is available for Microsoft Windows,
Sun Microsystems Solaris, and Linux operating systems.
14. SUBJECT TERMS
Hydrology, watershed, precipitation runoff, river routing, flood frequency, flood control, water
supply, computer simulation, environmental restoration.
20. LIMITATION OF
ABSTRACT
Unclassified
Unclassified
NSN 7540-01-280-5500
Unclassified
Standard Form 298 (Rev. 2-89)
Prescribed by ANSI Std. Z39-18 298-102
316
Unlimited
USAPPC V1.00
User's Manual
Version 3.5
August 2010
Please recycle this document when you are finished using it.
ii
Table of Contents
PREFACE
INTRODUCTION
SCOPE ................................................................................................................................................. 1
HISTORY .............................................................................................................................................. 1
CAPABILITIES........................................................................................................................................ 2
Watershed Physical Description .................................................................................................... 2
Meteorology Description ................................................................................................................ 3
Hydrologic Simulation .................................................................................................................... 4
Parameter Estimation .................................................................................................................... 4
Analyzing Simulations .................................................................................................................... 4
GIS Connection .............................................................................................................................. 4
LIMITATIONS ......................................................................................................................................... 5
Model Formulation ......................................................................................................................... 5
Flow Representation ...................................................................................................................... 5
DOCUMENTATION CONVENTIONS ........................................................................................................... 6
REFERENCES ....................................................................................................................................... 6
OVERVIEW
17
iii
Table of Contents
39
PROJECTS .......................................................................................................................................... 39
Creating a New Project................................................................................................................. 39
Opening a Project ......................................................................................................................... 40
Copying a Project ......................................................................................................................... 41
Renaming a Project ...................................................................................................................... 42
Deleting a Project ......................................................................................................................... 43
Project Properties ......................................................................................................................... 43
DIRECTORIES AND FILES...................................................................................................................... 44
Files Generated by the Program .................................................................................................. 44
Files Specified by the User ........................................................................................................... 45
Manually Entered Time-Series and Paired Data .......................................................................... 46
Computed Results ........................................................................................................................ 46
External Time-Series, Paired, and Grid Data ............................................................................... 46
Security Limitations ...................................................................................................................... 47
CONTROL SPECIFICATIONS .................................................................................................................. 47
Creating a New Control Specifications ......................................................................................... 47
Copying a Control Specifications ................................................................................................. 48
Renaming a Control Specifications .............................................................................................. 49
Deleting a Control Specifications ................................................................................................. 50
Time Window ................................................................................................................................ 52
Time Interval ................................................................................................................................. 53
Gridded Data Output Interval........................................................................................................ 53
IMPORTING HEC-1 FILES..................................................................................................................... 54
Selecting and Processing a File ................................................................................................... 54
Unsupported Features .................................................................................................................. 55
REFERENCES ...................................................................................................................................... 57
59
iv
Table of Contents
95
SUBBASIN ELEMENTS
125
Table of Contents
REACH ELEMENTS
155
RESERVOIR ELEMENTS
165
vi
Table of Contents
189
METEOROLOGY DESCRIPTION
199
vii
Table of Contents
HYDROLOGIC SIMULATION
231
PARAMETER ESTIMATION
253
viii
Table of Contents
279
293
297
299
301
305
ix
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.
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 algorithms from
legacy software. The current research program is designed to produce new
algorithms and analysis techniques for addressing emerging problems.
The initial program release was called Version 1.0 and included most of the eventsimulation capabilities of the HEC-1 program. It did introduce several notable
improvements over the legacy software including an unlimited number of hydrograph
ordinates and gridded runoff representation. The tools for parameter estimation with
optimization were much more flexible than in previous programs. The maiden
release also included a number of "firsts" for HEC including object-oriented
development in the C++ language and multiplatform support in a program with a
graphical user interface.
The second major release was called Version 2.0 and focused on continuous
simulation. The addition of the soil moisture accounting method extended the
Chapter 1 Introduction
program from an event-simulation package to one that could work equally well with
event or continuous simulation applications. The reservoir element was also
expanded to include physical descriptions for an outlet, spillway, and overflow. An
overtopping dam failure option was also added.
The third-party graphics libraries used to implement the multi-platform interface used
in Version 1.0 and Version 2.0 were sold and soon became unavailable. Faced with
the prospect of using unsupported graphics tools, the design team evaluated
alternatives and chose to move the program to the Java language. The simulation
engine was converted to Java and a completely new interface was designed.
During the process, careful attention was paid to lessons learned from earlier
versions. The result was Version 3.0 with the new interface plus new simulation
capabilities for infiltration, reservoir outlet structures, piping dam failure, and analysis
tools.
Enhancement of the program is ongoing. HEC has a strong commitment to
continued research in emerging needs for hydrologic simulation, both in terms of
simulation techniques and representation of physical processes. Future needs are
identified by conducting our own application projects, speaking with program users,
and monitoring academic journals. HEC also has a strong commitment to continued
development of the program interface. Plans are already underway to add new
features in a future version that will make the program easier to use by providing
more flexible ways to accomplish work. New visualization concepts are also being
developed. Look for future versions to continue the tradition.
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.
Chapter 1 Introduction
Meteorology Description
Meteorologic data analysis is performed by the meteorologic model and includes
precipitation, evapotranspiration, and snowmelt. Six different historical and synthetic
precipitation methods are included. Three evapotranspiration methods are included
at this time. Currently, only two snowmelt methods are available.
Four different methods for analyzing historical precipitation are included. The userspecified 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 automatically proceed when
missing data is encountered. The gridded precipitation method uses radar rainfall
data.
Four different methods for producing synthetic precipitation are included. The
frequency storm method uses statistical data to produce balanced storms with a
specific exceedance probability. Sources of supporting statistical data include
Technical Paper 40 (National Weather Service, 1961) and NOAA Atlas 2 (National
Weather Service, 1973). While it was not specifically designed to do so, data can
also be used from NOAA Atlas 14 (National Weather Service, 2004ab). The
standard project storm method implements the regulations for precipitation when
estimating the standard project flood (Corps of Engineers, 1952). The SCS
hypothetical storm method implements the primary precipitation distributions for
design analysis using Natural Resources Conservation Service (NRCS) criteria (Soil
Conservation Service, 1986). The user-specified hyetograph method can be used
with a synthetic hyetograph resulting from analysis outside the program.
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 and total volume. 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.
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
Chapter 1 Introduction
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.
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. This program does not currently include such
coupling of the process models. Errors due to the use of uncoupled models are
minimized as much as possible by using a small time interval for calculations. While
preparations have been made to support the inclusion of 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
Chapter 1 Introduction
connection so it is not possible to split the outflow from an element into two different
downstream elements. The diversion element provides a limited capability to remove
some of the flow from a stream and divert it to a different location downstream in the
network. Likewise, a reservoir element may have an auxiliary outlet. However, in
general, branching or looping stream networks cannot be simulated with the program
and will require a separate hydraulic model which can represent such networks.
The design of the process for computing a simulation does not allow for backwater in
the stream network. The compute process begins at headwater subbasins and
proceeds down through the network. Each element is computed for the entire
simulation time window before proceeding to the next element. There is no iteration
or looping between elements. Therefore, it is not possible for an upstream element
to have knowledge of downstream flow conditions, which is the essence of backwater
effects. There is a limited capability to represent backwater if it is fully contained
within a reach element. However, in general, the presence of backwater within the
stream network will require a separate hydraulic model.
Documentation Conventions
The following conventions are used throughout the manual to describe the graphical
user interface:
Menu names, menu items, and button names are shown in bold.
Data typed into an input field or selected from a list is shown using the
courier font.
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.
Hydrologic Engineering Center. April 1992. HEC-IFH Interior Flood Hydrology
Package: User's Manual. U.S. Army Corps of Engineers, Davis, CA.
Hydrologic Engineering Center. November 1989. Water Control Software: Forecast
and Operations. U.S. Army Corps of Engineers, Davis, CA.
National Weather Service. 1961. Technical Paper 40: Rainfall Frequency Atlas for
the United States for Durations from 30 Minutes to 24 Hours and Return Periods from
1 to 100 Years. U.S. Department of Commerce, Washington, DC.
National Weather Service. 2004. NOAA Atlas 14 Precipitation-Frequency Atlas of
the United States: Volume 1 Semi Arid Southwest (Arizona, Southeast California,
Nevada, New Mexico, Utah). U.S. Department of Commerce, Silver Spring, MD.
Chapter 1 Introduction
CHAPTER 2
Solaris 10 UltraSPARC.
The program has been extensively tested on Windows XP, Solaris 10 Update 6
UltraSPARC, and RED HAT Enterprise Linux 4.
1 GB of memory recommended.
The typical hardware equipment for the Sun Microsystems Solaris installation
includes:
1 GB memory recommended.
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.
10
5. The next window is used to select the location where the program will be
installed on the local disk. It is recommended that the default location in the
C:\Program Files folder be used. Press the Next button when you are
satisfied with the installation location.
6. The next window allows you to choose if a shortcut to the program will be
placed on the desktop. Program shortcuts will automatically be created in
the Start menu under the All Programs HEC folder, so having a desktop
shortcut is optional.
7. The next and final window allows you to confirm that you are ready for
installation. If you would like to change any of the previously configured
settings, you may use the Back button. Press the Install button to install the
program with all of the configuration information specified in the previous
windows.
The installer will copy all necessary files and make additional configuration changes
to the operating system. You do not need to restart the computer after the
installation completes. At any time you can uninstall the program through the Control
Panel. When future versions of the program become available, you may have each
version separately installed on your computer.
11
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 3.4 will be extracted to
<current working directory>/hechms34.
HEC-HMS 3.4 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.
12
13
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>)
where <HEC-HMS Install> is replaced with the installation directory name. Once
your path is set to include the HEC-HMS installation directory, then you can simply
type hms from anywhere on the system to run program.
14
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 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.
computer using Microsoft Windows or Linux can usually use up to 1,350 Mb. A
computer using the 32 bit version of Sun Microsystems Solaris can often use 3,000
MB while the 64 bit version can use hundreds of gigabytes. These are general
guidelines and your situation will depend on the specifics of your hardware and other
processes that may be executing at the same time as the program. In no case
should you attempt to use more than half of the physical memory in the machine
since other applications and system processes also require memory resources.
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 kept in a space-saving, compressed form. They can be extracted for use
at any time by going to the Help menu and selecting the Install Sample Projects
command.
15
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 way to use the program is also
provided.
Program Screen
The program screen contains a title bar, menu s-ystem, 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.
Desktop
Watershed
Explorer
Component
Editor
Message Log
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.
17
Chapter 3 Overview
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.
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 Desktop. A list of View menu
items and the resulting actions are provided in Table 3.
Table 1.
Open
Open a project.
Import
Save
Save As
Delete
Rename
Backup
Exit
Table 2.
18
Action
Action
Cut
Copy
Paste
Select All
Clear Selection
Select Special
Chapter 3 Overview
Table 3.
Action
Maximum Extents
Background Maps
Draw Gridlines
Zoom In
Zoom Out
Zoom To Selected
Rescale Elements
Clear Messages
Table 4.
Action
19
Chapter 3 Overview
The Compute menu contains a list of commands for creating, computing, and
managing simulation runs, optimization trials, and analyses. A list of Compute menu
items and the resulting actions are provided in Table 5. Menu commands for creating
simulation runs, optimization trials, and analyses are only enabled with the necessary
underlying components have been created. The Check Parameters and Compute
commands are only enabled when a compute item is selected.
The Results menu contains a list of commands for viewing results from hydrologic
elements in a basin model. Three cases must be met before commands from the
Results menu are active:
1. A simulation run, optimization trial, or analysis must be selected.
2. A hydrologic element or elements must be selected.
3. The results for the simulation must be available.
Results for a simulation are available as long as the simulation computed
successfully and no changes were made by the user to components used by the
simulation after it was computed. If you change parameter data, like values in a
paired data table used by a reach element for modified-Puls routing, then the
simulation must be recomputed. A list of Results menu items and the resulting
actions are provided in Table 6.
The Tools menu offers a list of specialized commands. A list of Tools menu items
and the resulting actions are provided in Table 7. Additional tools will be added to
the menu in a future program version.
Table 5.
20
Action
Select Run
Run Manager
Select Trial
Trial Manager
Create Analysis
Select Analysis
Analysis Manager
Check Parameters
Compute
Chapter 3 Overview
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.
Hydrograph Comparison
Flow Comparison
Flow Residuals
Action
Objective Function
Element Graph
Graph Properties
Table 7.
Action
Reports
Delete Results
Project Options
Program Settings
The Help menu provides links to the online help system. The sample projects can be
installed from this menu. You may also access the terms and conditions for using the
program, and information about the program and development team; this information
is helpful when reporting program problems.
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.
21
Chapter 3 Overview
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 Basin Models folder is selected, the Watershed Explorer expands to show
Table 8.
Action
Create a new project.
Open an existing project.
Save the current project.
Print the selected item in the Desktop (basin map or result window).
Select hydrologic elements in the basin map.
Pan in the basin map.
Zoom in or out in the basin map.
Add a subbasin element to the basin map.
Add a reach element to the basin map.
Add a reservoir element to the basin map.
Add a junction element to the basin map.
Add a diversion element to the basin map.
Add a source element to the basin map.
Add a sink element to the basin map.
Compute the currently selected simulation run.
Compute the currently selected optimization trial.
Compute the currently selected analysis.
Open global summary table.
Open graph for the current element selection.
Open summary table for the current element selection.
Open time-series table for the current element selection.
22
Chapter 3 Overview
Figure 2.
all basin models in the project. The Watershed Explorer will expand again to show all
hydrologic elements in a basin model if a basin model is selected. Notice in Figure 2
that the Castro 1 basin model is selected and the Watershed Explorer is expanded to
show all hydrologic elements in the basin model. The plus/minus sign beside model
components and sub-components can be used to expand or collapse the Watershed
Explorer. The Compute tab of the Watershed Explorer contains all project
simulation runs, optimization trials, and analyses. Model results are available from
the Results tab of Watershed Explorer.
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, river 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.
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
23
Chapter 3 Overview
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.
24
Figure 3.
Basin map for a basin model named Castro 2. The basin map is
shown in the Desktop area of the program screen.
Figure 4.
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.
Figure 5.
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, taking certain other automatic actions,
where to display global editors, and a display option for the decimal separator.
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 "Results" tab contains a number of settings for managing simulations and the
results that are produced. Control is provided for supplemental results that provide
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.
25
Chapter 3 Overview
Figure 6.
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
26
Chapter 3 Overview
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.
Number Formatting
Each country of the world has socially accepted conventions for formatting numbers.
The combination of spoken language and country are combined and called a locale.
For example, there is a locale for English conventions in Canada and also a locale for
French conventions in Canada, since both languages are commonly used in that
country. Locales have been defined for almost every language and country
combination in the world. Each locale includes the language and conventions for
formatting numbers. For example, in the United States a number would be displayed
as follows:
12.34
In Austria the same number would usually be displayed as follows:
12,34
Either format can be selected using the Program Settings, as shown in Figure 6. The
setting is used to interpret all user input and to configure all displayed output data.
27
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.
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.
28
Chapter 3 Overview
Table 9.
Time-Series Data
Paired Data
Grid Data
Precipitation
Storage-discharge
Precipitation
Discharge
Elevation-storage
Temperature
Stage
Elevation-area
Solar radiation
Temperature
Elevation-discharge
Crop coefficient
Solar radiation
Inflow-diversion
Storage capacity
Windspeed
Diameter-percentage
Percolation rate
Air pressure
Cross sections
Storage coefficient
Humidity
Moisture deficit
Altitude
Percentage curves
Impervious area
Crop coefficient
ATI-meltrate functions
ATI-coldrate functions
Elevation
Sediment Load
Groundmelt patterns
Cold content
Concentration
Meltrate patterns
Percent
Meltrate ATI
Liquid water content
Snow water equivalent
Water content
Water potential
29
Chapter 3 Overview
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 timeseries data.
Figure 7.
30
Chapter 3 Overview
Table 10.
Hydrologic
Element
Description
Subbasin
Reach
Junction
Source
Sink
Reservoir
Diversion
Figure 8.
31
Chapter 3 Overview
Figure 9.
Global editor for the initial and constant loss method. Only the
elements selected in the basin map are shown, or all elements if none
were initially selected. 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.
Most hydrologic elements require parameter data so that the program can model the
hydrologic processes represented by the element. In the case of the subbasin
element, many mathematical models are available for determining precipitation
losses, transforming excess precipitation to streamflow at the subbasin outlet, and
adding baseflow. In this document the different mathematical models will be referred
to as methods. The available methods for subbasin and reach elements are shown
in Table 11. Parameter data is entered in the Component Editor. Select a hydrologic
element in the basin map or Watershed Explorer to open the correct Component
Editor as shown in Figure 8. Global parameter editors can also be used to enter or
view parameter data for many hydrologic elements as shown in Figure 9. Global
parameter editors are opened using the Parameters menu.
32
Chapter 3 Overview
Table 11.
Hydrologic Element
Subbasin
Calculation Type
Method
Canopy
Surface
Loss Rate
Transform
Baseflow
Bounded recession
Constant monthly
Linear reservoir
Nonlinear Boussinesq
Recession
Reach
Routing
Kinematic wave
Lag
Modified Puls
Muskingum
Muskingum-Cunge
Straddle stagger
Gain/Loss
Constant
Percolation
33
Chapter 3 Overview
Table 12.
Precipitation Methods
Description
Frequency storm
Gage weights
Gridded precipitation
Inverse distance
SCS storm
Applies a user specified SCS time distribution to a 24hour total storm depth.
Specified hyetograph
Figure 10.
34
Chapter 3 Overview
Figure 11.
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 any hydrologic element 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.
35
Chapter 3 Overview
Results can also be viewed from the Watershed Explorer, Results tab. Select the
simulation run and the Watershed Explorer will expand to show all hydrologic
elements in the basin model. If you select one of the hydrologic elements, the
Watershed Explorer expands again to show all result types as shown in Figure 13.
For a subbasin element, you might see outflow, incremental precipitation, excess
precipitation, precipitation losses, direct runoff, and baseflow as the output results.
Select one of these results to open a preview graph. Multiple results can be selected
and viewed by holding down the Control or Shift buttons. Results from multiple
hydrologic elements can be viewed together. Also, results from different simulation
runs can be selected and viewed. Once output types are selected in the Watershed
Explorer, a larger graph or time-series table can be opened in the Desktop by
selecting the Graph and Time-Series buttons on the toolbar.
Figure 12.
36
Accessing results for the current simulation run using the basin map.
Chapter 3 Overview
Figure 13.
37
Chapter 3 Overview
Figure 14.
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.
38
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?
What is the impact of changing the land use on 7% of the watershed from
cow pasture to home sites?
39
Figure 15.
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.
Figure 16.
40
Figure 17.
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 computer to store the project, and describe
the project as shown in Figure 18. A default location is provided. The default
location for creating new projects is specified in the program settings. Access the
program settings by clicking the Tools menu and selecting the Program Settings
command.
Figure 18.
41
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 paste them 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 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.
Figure 19.
42
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 13 for a complete list of project files). No external files are
deleted when this option is selected. External files include DSS, grid cell, and 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
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
project, press the Cancel button or the X button in the upper right to close the
window.
Figure 20.
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).
The description can be used to provide internal documentation. It is helpful to
include notes and comments within the project to remind yourself of details at a later
time. It also helps with providing information to other users who may work with the
project in the future. 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.
43
Figure 21.
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.
44
Table 13.
Description
Castro.access
Castro.dsc
Castro.dss
Castro.gage
Castro.grid
Castro.hms
Castro.log
Castro.nals
Castro.out
Castro.pdata
Castro.run
North_Branch.basin
Historic_Calib.met
Oct_1977.control
Castro.trial
Est_Baseflow.optim
Castro.stateIndex
Mid_Oct_77.state
45
can be stored at any location on the computer file system, but it is often convenient to
store them in the project directory. Additional supplementary files, related to the
project but not used by the program, can also be placed in the project directory and
will be ignored. However, when a project is copied, only files used by the program
will be copied.
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 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/
A consistent naming convention is used for assigning the different pathname parts of
the computation results (HEC, 1994). The B-part is assigned the name of the
element in the basin model. The C-part is assigned a data descriptor as described in
Appendix B. The D-part is assigned the simulation start date. The E-part is assigned
the simulation time interval. The F-part begins with a three-letter code that identifies
the type of computed result, followed by a colon and the name of the compute
component.
46
external files can store regular or irregular interval time-series data, paired data, or
grid data. They can be located anywhere on the computer or network and shared
with other programs. This program automatically determines the data type, units,
and interval from the record header.
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 accessible but becomes inaccessible, it is possible that the file
permissions were changed external to the program. It is also possible for the files to
have read-only permission if they are copied from a CD-ROM or other removable
storage media. If the files have been set to read-only permission for any reason, you
will need to manually change the permissions on the files before the program can
open the project. If you are using the Microsoft Windows operating system, you can
change file permissions using Windows Explorer. On the Sun Microsystems
Solaris and Linux operating systems, you can open a command window and use
the chmod command.
The program is designed to work with projects that may be shared by several users.
Usually shared projects will be stored on a server or network storage device, but
could be stored in a shared folder on a local computer. All users who will share the
project must have read and write permission for that folder. Even though several
users may share the project, only one user can access the project at a time. The
program automatically tracks how many users are accessing a project and limits
access as necessary.
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.
47
name is provided for the new specifications; 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 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.
Figure 22.
Figure 23.
48
The second way to copy is from the "Components" tab of the Watershed Explorer.
Move the mouse over the control specifications you wish to copy, then press the right
mouse button (Figure 24). A context menu is displayed that contains several choices
including copy. Click the Create Copycommand. A new Copy Control
Specifications 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 control specifications. You cannot press the Copy
button if no name is specified. If you change your mind and do not want to copy the
selected control specifications, press the Cancel button or the X button in the upper
right of the Copy Control Specifications window to return to the Watershed Explorer.
Figure 24.
Figure 25.
The first way to perform a rename is to use the Control Specifications Manager,
which you can access from the Components menu. Select the control specifications
you wish to rename by clicking on it in the list of current control specifications. The
selected specifications is highlighted after you select it. After you select a
specifications you can press the Rename button on the right side of the window. A
new Rename Control Specifications window (Figure 25) will open where you can
provide the new name. If you wish you can also change the description at the same
49
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
control specifications. You cannot press the Rename button if no name is specified.
If you change your mind and do not want to rename the selected control
specifications, press the Cancel button or the X button in the upper right of the
Rename Control Specifications window to return to the Control Specifications
Manager window.
The second way to rename is from the "Components" tab of the Watershed Explorer.
Select the control specifications you wish to rename by clicking on it in the
Watershed Explorer; it will become highlighted. Keep the mouse over the selected
specifications and click the right mouse button. Select the Rename command from
the menu and the highlighted name will change to editing mode as shown in Figure
26. 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 control specifications,
press the Escape key.
Figure 26.
50
control specifications, press the Cancel button or the X button in the upper right to
return to the Control Specifications Manager window.
Figure 27.
The second way to delete is from the "Components" tab of the Watershed Explorer.
Select the control specifications you wish to delete by clicking on it in the Watershed
Explorer; it will become highlighted. Keep the mouse over the selected specifications
and click the right mouse button (Figure 28). 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 specifications.
Press the OK button to delete the specifications. If you change your mind and do not
want to delete the selected control specifications, press the Cancel button or the X
button in the upper right to return to the Watershed Explorer.
Figure 28.
51
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 1 AD through 4000 AD. 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.
Figure 29.
52
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 timeseries 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 24
hours.
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.
53
The gridded state variables will be stored at the gridded data output interval.
However, the internal calculations will be performed at the time interval, and the timeseries results will use the time interval. For example, a time interval of 1 hour might
be specified for a continuous simulation that spans several months. A gridded data
output interval of 1 day could be selected in order to look at the gridded results at the
end of each day.
54
Figure 30.
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 modern engineering analysis methods for riskbased design. Other features have not been added to this program because they
use old numerical algorithms that have been replaced by superior methods in this
program and no direct translation is available. Other features may be added to this
program in the future. The so-called cards from HEC-1 that are not supported during
import are shown in Table 14. An error message will be displayed if any of those
cards are encountered.
55
Table 14.
Record Identifier
56
Alternative
IO, VS, VV
JP, KP
JR
JD
OU, OR
HB
HL
HS
KF
All computed results are stored in the project DSS file using double
precision. There is no need to specify output format.
LM, LH
PM
QP
RC
Normal depth routing to determine the modified Puls storageoutflow curve is not supported. Channel data are converted to the
Muskingum-Cunge method.
RL
UA
WP, WR
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.
57
CHAPTER 5
Time-Series Data
Hydrologic models often require time-series of precipitation data for estimating basinaverage 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.
Figure 31.
59
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 timeseries data manager.
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.
60
Creating a copy of a gage after pressing the Copy button in the TimeSeries Data Manager.
Figure 34.
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 timeseries 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 timeseries gage you want to rename from the Data Type menu. Then, select the timeseries 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
61
Figure 35.
Renaming a gage after pressing the Rename button in the TimeSeries Data Manager.
Figure 36.
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 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.
62
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 timeseries gage you want to delete from the Data Type menu. Then, select the timeseries 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.
The second way to delete a gage is from the Watershed Explorer, on the
Components tab (Figure 38). Select the time-series gage you want to delete by
clicking on it in the Watershed Explorer; it will become highlighted. Keep the mouse
over the selected gage and click the right mouse button. 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 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 Watershed Explorer.
63
Figure 38.
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 1 AD through 4000 AD. 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
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 TimeSeries Data Manager, accessed by clicking the Components menu and then
64
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
to 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.
Figure 39.
Creating a new time window for a gage, beginning from the TimeSeries Data Manager.
Figure 40.
65
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 Figure 42; click the Delete Time
Window command to delete the selected time window.
66
Figure 41.
Figure 42.
Figure 43.
Using the time window component editor for a gage to view and edit
the start and end time for a window.
You 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.
Data Source
The data source determines how the data for a time-series component will be stored.
Data may be entered manually or retrieved from the Data Storage System (HECDSS). Manual entry means that you must enter values for all of the time windows.
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 timeseries gage will be determined by the data source, as discussed in the following
sections. Compare Figure 44 and Figure 45 to see the difference in entering data for
manual entry or HEC-DSS data sources.
Figure 44.
67
Figure 45.
Data Units
The data units can only be selected for a manual entry time-series gage (Figure 44);
they are retrieved automatically for the DSS option (Figure 45). Most types of timeseries 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). The
precipitation time-series type has additional options for specifying incremental or
cumulative data. The units available in the "Units" field will depend on the time-series
type of the selected component. All time windows defined for a time-series
component must use the same time interval.
Generally you should choose the data units before entering any data for the gage.
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 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 data units if the manual entry
option is selected.
Time Interval
The time interval can only be selected for a manual entry time-series gage (Figure
44); it is retrieved automatically for the DSS option (Figure 45). An interval must be
selected from the available choices that range from 1 minute to 24 hours. 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
68
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 "TimeSeries Gage" tab in the Component Editor will display the time interval if the manual
entry option is selected.
69
Table 15.
Internal DSS data type label for different types of time-series data.
Time-Series Type
Precipitation
Discharge
Stage
Temperature
Solar Radiation
70
Label
Description
PER-CUM
INST-CUM
PER-AVER
INST-VAL
PER-AVER
INST-VAL
PER-AVER
INST-VAL
PER-AVER
Windspeed
INST-VAL
Air pressure
INST-VAL
Humidity
INST-VAL
Altitude
INST-VAL
Crop Coefficient
INST-VAL
Snow Water
Equivalent
INST-VAL
Sediment Load
PER-CUM
Concentration
INST-VAL
Percent
INST-VAL
Table 16.
Time-Series Type
Precipitation
Label
MM
IN
Discharge
Stage
Temperature
Solar Radiation
Windspeed
Air Pressure
Description
Millimeters
Inches
M3/S
CFS
Meters
FT
Feet
DEG C
Degrees centigrade
DEG F
Degrees Fahrenheit
WATT/M2
LANG/MIN
KPH
MPH
KPA
Kilo pascals
IN HG
Inches of mercury
Humidity
Relative humidity
Altitude
KM
MILE
Crop Coefficient
Snow Water Equivalent
UNSPECIF
MM
IN
Sediment Load
TONS
TONNES
Concentration
Percent
MG/L
%
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 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.
71
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 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 linearly interpolating the values between the
first and last cell in the selection, 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 window under a time-series gage in the Watershed Explorer, the
"Graph" tab will be added to the Component Editor. Data in the graph cannot be
edited regardless of whether the gage uses manual entry or retrieves data from a
72
DSS file. If no time-series data is available for the specified time window, then the
graph will not contain any data.
Figure 47.
Elevation
Some calculations cannot be performed unless the elevation of the gage is specified.
The only gages where this may need to be entered are temperature gages, air
73
pressure gages, and humidity gages. The elevation must be specified when one of
these gages is used in snowmelt calculations; use in other calculations does not
require the elevation to be specified. When required, the elevation should be the
ground elevation above sea level at the location where the measurement is collected.
The value may be entered as either meters or feet, depending on the unit system of
the gage. An example of entering the elevation of a gage can be seen in Figure 44.
Reference Height
Most meteorologic monitoring equipment is installed on a tower or some other
elevated structure. The reference height specifies the distance between the
instrument and the ground surface. The value can be entered in meters or feet,
depending on the unit system of the gage. The reference height is only required for
windspeed gages and humidity gages.
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.
74
Figure 48.
Figure 49.
Creating a new cross section after pressing the New button in the
Paired Data Manager.
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
75
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 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.
Figure 51.
76
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.
Figure 52.
Figure 53.
77
The second way to rename is from the Watershed Explorer, on the Components
tab. Move the mouse over the paired data you wish to copy, then press the right
mouse button (Figure 53). 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 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 curve, press the Escape key.
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 retrieved or undeleted. Any references
to the deleted curve will switch to using no curve, 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 curve.
The first way to perform a deletion is to use the paired data manager, which you can
access from the Components menu. First, select the data type for the paired data
curve you want to delete from the Data Type menu. Then, select the curve you want
to delete 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 Delete
button on the right side of the window (Figure 54). A window will open where you
must confirm that you want to delete the selected curve. Press the OK button to
delete the curve. If you change your mind and do not want to delete the selected
curve, press the Cancel button or the X button in the upper right to return to the
paired data manager.
Figure 54.
78
Figure 55.
The second way to delete a paired data is from the Watershed Explorer, on the
Components tab (Figure 55). Select the curve you want to delete by clicking on it in
the Watershed Explorer; it will become highlighted. Keep the mouse over the
selected curve and click the right mouse button. 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 curve. Press
the OK button to delete the curve. If you change your mind and do not want to delete
the selected curve, press the Cancel button or the X button in the upper right to
return to the Watershed Explorer.
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.
Figure 56.
79
Figure 57.
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 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.
80
81
Table 17.
Preferred data order and internal DSS data type label for different
types of paired data.
Description
Storage-Discharge
UNT
Elevation-Storage
UNT
Elevation-Area
UNT
Elevation-Discharge
UNT
Inflow-Diversion
UNT
Diameter-Percentage
UNT
Cross Section
UNT
INST-VAL
Unit Hydrograph
82
Label
Percentage Curve
UNT
ATI-Meltrate
UNT
ATI-Coldrate
UNT
Groundmelt Pattern
INST-VAL
Meltrate pattern
INST-VAL
Table 18.
Internal DSS units for different types of paired data when the preferred
order is followed. The unit hydrograph and patterns use dates in the
year 2000 to specify information for the x-axis.
Elevation-Storage
Elevation-Area
Elevation-Discharge
Inflow-Diversion
Diameter-Percentage
Cross Section
X-Label
Y-Label
THOU M3
M3/S
ACRE-FT
CFS
THOU M3
FT
ACRE-FT
THOU M2
FT
ACRE
M3/S
FT
CFS
M3/S
M3/S
CFS
CFS
MM
IN
FT
FT
Unit Hydrograph
Percentage Curve
ATI-Meltrate
ATI-Coldrate
Groundmelt Pattern
M3/S
CFS
DEGC-D
MM/DEG-D
DEGF-D
IN/DEG-D
DEG C
MM/DEG-D
DEG F
IN/DEG-D
MM/DAY
IN/DAY
Meltrate Pattern
Description
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.
83
Figure 58.
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
linearly interpolating the values between the first and last cell in the selection,
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.
Figure 59.
84
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.
Figure 60.
85
Figure 61.
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 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
Figure 62.
86
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.
Figure 63.
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
87
Figure 64.
Figure 65.
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 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, 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
88
Figure 66.
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.
Figure 67.
89
Figure 68.
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/
The convention for storing gridset data is different because a separate grid is
required for each time interval. The different records will each have all the same
pathname parts except for the D-part and E-part. The D-part indicates the starting
time and the e-part the ending time of the individual record. Any of the record
pathnames can be selected and the program will automatically retrieve the correct
data depending on the time window during a compute. The following is an example
of a temperature gridset:
//LOUISVILLE/TEMP/12JUL2003:1230/12JUL2003:1300/OBS/
90
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 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 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 grid data are shown in Table 19. The correct unit
labels are shown in
Table 20.
Table 19.
Internal DSS data type label for different types of grid data.
Time-Series Type
Label
Description
Precipitation
PER-CUM
Temperature
PER-AVER
INST-VAL
PER-AVER
INST-VAL
Solar Radiation
91
Label
Description
PER-AVER
INST-VAL
Storage Capacity
INST-VAL
Percolation Rate
INST-VAL
Storage Coefficient
INST-VAL
Moisture Deficit
INST-VAL
Impervious Area
INST-VAL
INST-VAL
Elevation
INST-VAL
Cold Content
INST-VAL
INST-VAL
Meltrate ATI
INST-VAL
Liquid Water
INST-VAL
Crop Coefficient
92
Label
Description
Snow Water
Equivalent
INST-VAL
Water Content
INST-VAL
Water Potential
INST-VAL
Table 20.
Label
MM
IN
Temperature
Solar Radiation
Description
Millimeters
Inches
DEG C
Degrees Celsius
DEG F
Degrees Fahrenheit
WATT/M2
LANG/MIN
Crop Coefficient
UNDEF
Storage Capacity
MM
Percolation Rate
MM/HR
IN
IN/HR
Dimensionless coefficient
Millimeters
Inches
Millimeters per hour
Inches per hour
Storage Coefficient
HR
Coefficient in hours
Moisture Deficit
MM
Millimeters
Impervious Area
SCS Curve Number
Elevation
Cold Content
IN
Inches
Percent of area
UNDEF
M
Meters
FT
Feet
MM
Millimeters
IN
Cold Content ATI
Meltrate ATI
Liquid Water Content
Inches
DEG C
Degrees Celsius
DEG F
Degrees Fahrenheit
DEGC-D
DEGF-D
MM
IN
Dimensionless parameter
MM
IN
Water Content
UNDEF
Water Potential
MM
IN
Millimeters
Inches
Millimeters
Inches
Fraction of saturation
Millimeters
Inches
93
References
Hydrologic Engineering Center. January 2003. HEC Data Storage System Visual
Utility Engine: User's Manual. U.S. Army Corps of Engineers, Davis, CA.
94
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.
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.
95
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.
Figure 70.
The second way to copy is from the "Components" tab of the Watershed Explorer.
Move the mouse over the basin model you wish to copy, then press the right mouse
button (Figure 71). A context menu is displayed that contains several choices
including copy. Click the Create Copy command. A new Copy Basin Model
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 basin model. You cannot press the Copy button if no name is specified. If
you change your mind and do not want to copy the selected basin model, press the
Cancel button or the X button in the upper right of the Copy Basin Model window to
return to the Watershed Explorer.
96
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.
Figure 72.
Renaming a basin model. This was accessed from the Basin Model
Manager.
97
Figure 73.
The second way to rename is from the "Components" tab of the Watershed Explorer.
Select the basin model you wish to rename by clicking on it in the Watershed
Explorer; it will become highlighted. Keep the mouse over the selected model and
click the left mouse button again. The highlighted name will change to editing mode
as shown in Figure 73. 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
basin model, press the Escape key.
98
Figure 74.
Figure 75.
99
Figure 76.
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. 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
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 downstream 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 some of all inflow. However,
in addition to computing outflow, it also computes the additional local flow. 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.
100
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 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.
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 timeseries 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
A future version of the program will provide the option of simulating the erosion of
sediment from the land surface. It will also provide the option of simulating the
transport of sediment through reaches, along with possible erosion or deposition in
the channel network. The simulation of sediment settling in reservoirs will also be
101
included. This selection option will be available when the capability is fully tested and
ready for use.
Water Quality
A future version of the program will provide the option of simulating nutrient loadings
from subbasin and source elements. It will provide for the simulation of nutrient
transport and fate in reaches as water moves through the channel network. The
simulation of nutrient transformations in reservoirs will also be included. This
selection option will be available when the capability is fully tested and ready for use.
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.
Adding, removing, and managing background maps is done from 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 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 must independently make sure all maps you wish to use are in the same
coordinate system.
102
Figure 77.
The current background maps are drawn in a specified order so that you can control
overlapping and other drawing issues. The draw order is based on the order they are
shown in the Background Maps window. The map at the bottom of the list is drawn
first. After it is drawn, the map shown above it in the list is drawn second. The map
shown at the top of the list is drawn last. The hydrologic elements are drawn after all
of the maps have been drawn. You may temporarily turn drawing for a map "on" or
"off" using the check box in the map list.
The draw order of the maps can be changed once all the desired background maps
are shown in the list of current maps. To change the order of a map you must first
select it in the list of current maps. Place the mouse over the desired map and click
the left mouse button; the map name is highlighted to show that it is selected. Press
the Move Up button to move the map up in the draw order so that other maps are
drawn first. Press the Move Down button to move the map down in the draw order
so that other maps are drawn after it. Each time you press the Move Up or Move
Down button the map is moved one position up or down. You can press the buttons
more than once to get the map to the desired position in the drawing order.
A background map can be removed. Place the mouse over the map you wish to
remove and select it by clicking the left mouse button; the map name is highlighted to
show that it is selected. Press the Remove button. The background map is not
deleted when it is removed from the list of current background maps. The files that
contain the actual background map data cannot be deleted from the program; they
must be deleted manually.
The drawing properties of a background map using the Shape File format can be
controlled. The Draw Properties button will be enabled when ever a map is selected
that uses the Shape File format. Pressing the button will open an editor that allows
control of drawing properties including colors, line style and width, shading,
transparency, and other properties. Changes may be previewed before they become
final.
103
Maximum Extents
The maximum extents are designed to limit the area shown in the Basin Model
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 by pressing the Set button.
The Union of All Elements option will set the extents to the combined extends
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.
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
104
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.
105
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 window will automatically zoom out to the coordinates specified as the
maximum extents and redraw.
Background Gridlines
Background gridlines are useful to helping to measure approximate distance in the
Basin Model 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.
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
106
atmospheric forcing. Seven different element types have been included in the
program: subbasin, reach, reservoir, junction, diversion, source, and sink.
An element uses a mathematical model to describe the physical process.
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 five different methods
for specifying the input data for a reservoir.
Figure 79.
107
subbasin, reach, reservoir, junction, source, diversion, or sink. After selecting the
appropriate tool, move the mouse over the Basin Model map; the mouse cursor
changes to cross hairs. Move the mouse until the cross hairs are over the location
where you wish to 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.
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.
108
Figure 80.
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 will 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.
109
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 context menu is displayed that contains several choices including cut.
Click the Cut Element command. You can now paste the selection into the same or
a different basin model. The cut command is only available if there is at least one
element selected.
110
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.
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. 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 context
menu is displayed that contains several choices including delete. Click the Delete
Element command. A window will open where you must confirm that you wish to
delete the selected hydrologic element. Press the OK button to delete the element or
elements. If you change your mind and do not want to delete the element selection,
press the Cancel button or the X button in the upper right to return to the Basin
Model map window.
The second way to delete an element also uses the selection in the map. Select one
or more elements that you wish to delete. Press the Delete key on the keyboard.
You will have to confirm your choice to delete.
The third way to delete is from the "Components" tab of the Watershed Explorer.
Move the mouse over the element you wish to delete and press the right mouse
button as shown in Figure 84. A context menu is displayed that contains several
111
choices including delete. Click the Delete command. You will have to confirm your
choice to delete.
112
Figure 83.
Figure 84.
Figure 85.
The program fundamentally computes flow for each of the hydrologic elements
included in a basin model. An elevation-discharge curve can be specified so that
stage can be computed as well. If a curve is specified, the stage for each time step is
113
determined by taking the computed flow and using it to interpolate elevation from the
curve. The computed stage is added to time-series results for the element and
appears in summary tables, time-series tables, and graphs.
The flow ratio is only shown for subbasin and source elements, and only if the basin
model is set to allow flow ratios. To enable or disable flow ratios, access the
Component Editor for the basin model and make the desired selection. When the
ratios are enabled in the basin model, the flow ratio can be entered for subbasin and
source elements. The ratio is applied to the computed flow from the subbasin or
source in order to compute the final outflow.
A reference flow can be specified to assist in interpreting computed flow results. The
flow is added to the element results graph as a horizontal marker line at the specified
flow value. The marker line is labeled with the specified label. The reference flow
can represent any significant flow value such as bank-full discharge, flood watch, or
levee overtopping.
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.
114
order under the basin model icon. Click on an element icon in the Watershed
Explorer and it will become selected. The selected element is highlighted in the
Watershed Explorer and in the Basin Model window. More than one element may be
selected at a time but only the Component Editor for the first selected element is
shown.
You can also find and select a hydrologic element in the Basin Model window. Begin
by making sure the Arrow Tool is selected on the toolbar. Examine the element icons
in the Basin Model 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 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 basin model 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 to 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.
115
Figure 87.
Using the special selection tool to select all junctions whose names
begins with the letters "West".
The third component of a special selection is the selection type to create. You may
replace the current selection. This means that what ever elements are found to meet
the name criteria and the element type criteria will become the selected elements.
Any previously selected elements will no longer be selected. You may alternately
add to the current selection. This means the new selection will include any elements
that are currently selected plus the elements that meet the new name and element
type criteria. Finally, you can select from the current selection. In this case the name
and element type criteria are further limited by only selecting from among the
elements that are already selected.
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 window; they are moved the
same way regardless of whether they are connected to a downstream element. To
move an element, start by selecting the Arrow Tool from the toolbar. Next click the
icon of the element you wish to move; it will become highlighted. Keep the mouse
over the element and hold the left mouse button. Drag the mouse until the element
icon is in the desired location. Release the left mouse button to finalize the move.
Reach elements are often connected between two junctions or possibly between
other element types. Reach elements must connect to the upstream and
116
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 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 applied 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. 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
117
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.
118
Figure 88.
Figure 89.
119
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 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 optionally for sources, is used to determine the
drainage area for any other element in the network. A simple basin model is shown
in Figure 90, with the corresponding Watershed Explorer shown in Figure 91.
Figure 90.
120
Figure 91.
121
position. You may also select a contiguous block of elements and move all of then
up or down position in unison.
Any selected element in the Watershed Explorer can be moved several positions at
once using the mouse. Start by clicking on an element to select it. Next, carefully
position the mouse over the selected element and hold the left mouse button, and
then drag the mouse. The shape of the cursor will change while you drag the
element and you will also see an insertion, as shown in Figure 93. The insertion line
shows where the element will be placed when you release the mouse button. Drag
the mouse up or down over the other elements until you position the mouse where
you want the element, and then release the button. The element will be moved to the
new position in the hydrologic order.
122
Figure 92.
Figure 93.
Global editors and simulation results shown throughout the program will display
elements in hydrologic order, as shown in the Watershed Explorer. Hydrologic order
will be used when the editor or result is first opened. However, each result also
includes a selection choice for switching between hydrologic order and alphabetic
order.
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.
123
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 94.
The canopy method for a subbasin is selected on the Component Editor for the
subbasin element (Figure 94). Access the Component Editor by clicking the
subbasin element icon on the "Components" tab of the Watershed Explorer. 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 canopy method from the list of three
available choices. If you choose the None method, the subbasin will not compute
any interception and all precipitation will be assumed as direct precipitation on the
ground surface, and subject to interception by the surface, and infiltration into the
125
soil. Use the selection list to choose the method you wish to use. Each subbasin
may use a different method or several subbasins may use the same method.
When a new subbasin is created, it is automatically set to use the default canopy
method specified in the project settings. You may change the canopy method for a
subbasin at any time using the Component Editor for the subbasin element. Since a
subbasin can only use one canopy 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 canopy method for several
subbasins simultaneously. Click on the Parameters menu and select the Subbasin
Methods Canopy command. The canopy method you choose will be applied to
the selected subbasins in the basin model, or to all subbasins if none are currently
selected.
The parameters for each canopy method are presented on a separate Component
Editor from the subbasin element editor. The "Canopy" editor is always shown next
to the "Subbasin" editor. If the kinematic wave transform method is selected, there
may be two canopy editors, one for each runoff plane. The information shown on the
canopy editor will depend on which method is currently selected.
Figure 95.
Simple Canopy
This method is a simple representation of a plant canopy. All precipitation is
intercepted until the 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. The canopy will consume all potential evapotranspiration until
the water in storage has been eliminated. Unused potential evapotranspiration after
126
the storage is emptied can be used by the surface and soil components. The canopy
either fills storage during precipitation or empties storage when it is not precipitating.
There is no evapotranspiration if there is precipitation. The Component Editor is
shown in Figure 96.
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.
Canopy storage represents the maximum amount of water that can be held on leaves
before through-fall to the surface begins. The amount of storage is specified as an
effective depth of water.
Figure 96.
127
warning in the program settings. You can change the surface method for several
subbasins simultaneously. Click on the Parameters menu and select the Subbasin
Methods Surface command. The surface method you choose will be applied to
the selected subbasins in the basin model, or to all subbasins if none are currently
selected.
The parameters for each surface method are presented on a separate Component
Editor from the subbasin element editor. The "Surface" editor is always shown next
to the "Canopy" editor. If the kinematic wave transform method is selected, there
may be two surface editors, one for each runoff plane. The information shown on the
surface editor will depend on which method is currently selected.
Figure 97.
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 when ever 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 98.
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.
128
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.
Figure 98.
129
Figure 99.
The maximum storage specifies the amount of water the soil layer can hold, specified
as a depth. An upper bound would be the depth of the active soil layer multiplied by
the porosity. However, in most cases such an estimate will have to be reduced by
the permanent wilting point and for other conditions that reduce the holding capacity
of the soil.
The constant rate defines the infiltration rate when the soil layer is saturated. A good
approximation is to use the saturated hydraulic conductivity.
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.
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 of accumulated infiltration. It includes the option for increased initial
infiltration when the soil is particularly dry before the arrival of a storm. Because it is
a function of cumulative infiltration and does not include any type of recovery, it
should not be used for continuous simulation. The Component Editor is shown in
Figure 100.
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.
130
Figure 100.
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.
The coefficient ratio indicates the rate at which the exponential decrease in infiltration
capability proceeds. It may be considered a function of the ability of the surface of a
subbasin to absorb precipitation and should be reasonable constant for large,
homogeneous areas.
The precipitation exponent reflects the influence of precipitation rate on subbasinaverage loss characteristics. It reflects the manner in which storms occur within an
area and may be considered a characteristic of a particular region. It varies from 0.0
up to 1.0.
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.
131
Figure 101.
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.
132
Figure 102.
The maximum storage grid ratio must be entered. The actual maximum storage
value at each grid cell is obtained by multiplying the value obtained from the
maximum storage 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.
133
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.
Figure 103.
134
Figure 104.
135
Figure 105.
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.
136
Figure 106.
Figure 107.
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.
137
Figure 108.
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 no temperature gage is
selected then a temperature of 25C (75F) is assumed to prevail. The gage must be
defined in the time-series manager before it can be selected in the component editor.
138
The beta zero parameter is used to correct the matric potential based on
temperature. It has been found to be a function of soil texture. It will only be shown
for input if a temperature gage has been selected.
Figure 109.
The initial condition of the soil should be specified as the percentage of the soil that is
full of water at the beginning of the simulation. The initial condition of the upper and
lower groundwater layers must also be specified.
The maximum infiltration rate sets the upper bound on infiltration from the surface
storage into the soil. This is the upper bound on infiltration; the actual infiltration 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.
The percentage of the subbasin which is directly connected impervious area can be
specified. All precipitation on that portion of the subbasin becomes excess
precipitation and subject to direct runoff.
139
Soil storage represents the total storage available in the soil layer. It may be zero if
you wish to eliminate soil calculations and pass infiltrated water directly to
groundwater.
Tension storage specifies the amount of water storage in the soil that does not drain
under the affects of gravity. Percolation from the soil layer to the upper groundwater
layer will occur whenever the current soil storage exceeds the tension storage.
Water in tension storage is only removed by evapotranspiration. By definition,
tension storage must be less that soil storage.
The soil percolation sets the upper bound on percolation from the soil storage into
the upper groundwater. The actual percolation rate is a linear function of the current
storage in the soil and the current storage in the upper groundwater.
Groundwater 1 storage represents the total storage in the upper groundwater layer.
It may be zero if you wish to eliminate the upper groundwater layer and pass water
percolated from the soil directly to the lower groundwater layer.
The groundwater 1 percolation rate sets the upper bound on percolation from the
upper groundwater into the lower groundwater. The actual percolation rate is a
linear function of the current storage in the upper and lower groundwater layers.
The groundwater 1 coefficient is used as the time lag on a linear reservoir for
transforming water in storage to become lateral outflow. The lateral outflow is
available to become baseflow.
Groundwater 2 storage represents the total storage in the lower groundwater layer. It
may be zero if you wish to eliminate the lower groundwater layer and pass water
percolated from the upper groundwater layer directly to deep percolation.
The groundwater 2 percolation rate sets the upper bound on deep percolation out of
the system. The actual percolation rate is a linear function of the current storage in
the lower groundwater layer.
The groundwater 2 coefficient is used as the time lag on a linear reservoir for
transforming water in storage to become lateral outflow. It is usually a larger value
that the groundwater 1 coefficient. The lateral outflow is likewise available to become
baseflow.
140
When a new subbasin is created, it is automatically set to use the default transform
method specified in the project settings. You may change the transform method for a
subbasin at any time using the Component Editor for the subbasin element. Since a
subbasin can only use one transform 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 transform method for several
subbasins simultaneously. Click on the Parameters menu and select the Subbasin
Methods Transform command. The transform method you choose will be
applied to the selected subbasins in the basin model, or to all subbasins if none are
currently selected.
The parameters for each transform method are presented on a separate Component
Editor from the subbasin element editor. The "Transform" editor is always shown
near to the "Subbasin" editor. The kinematic wave method uses multiple tabs in the
Component Editor. The information shown on the transform editor will depend on
which method is currently selected.
Figure 110.
141
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.
Figure 111.
Plane editor for the kinematic wave transform method. The same
properties are required for the second plane, if it is used.
The percentage of the subbasin area occupied by each plane must be entered. If
you only want to use one plane, then enter 100 for the first plane and 0 for the
second. You do not need to enter length, slope, or roughness if the area will be set
to zero.
The number of routing steps is used as a hint to the program when it determines the
correct distance step to use during runoff calculations. Criteria based on rate of
precipitation and other factors are used to automatically determine the correct
distance and time steps for solving the kinematic wave equation. The default value is
five.
The outflow from each plane is scaled up to the subbasin according the percentage
of the subbasin occupied by each plane type. The resulting composite outflow is
linked to a subcollector channel (Figure 112). The subcollector is intended to
represent primary collection systems such as street gutters but cross section choices
allow alternate representations. The outflow from the planes is applied to the
subcollector as lateral inflow.
142
Figure 112.
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.
The Manning's n roughness coefficient should be the average value for the whole
subcollector. It is important to remember that the parameter data entered for the
subcollector should be typical of all similar subcollectors in the subbasin.
The number of subreaches is used as a hint to the program when it determines the
correct distance step to use during routing calculations. The default value is five.
The typical area of each subcollector must be entered. This is used to determine
how to apply the composite outflow from the subcollector to the collector channel. It
essentially determines the number of subcollector channels in the subbasin.
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.
The outflow from the subcollector enters a collector as lateral inflow. The collector
channel is intended to represent small ditches or open channels that are part of an
engineered stormwater management system. The parameters for the collector are
143
exactly the same as for the subcollector. Length, slope, Manning's n, and the
number of subreaches are all the same. The area served by a typical collector must
be entered; this is used to apply lateral inflow from the collector to the main channel.
Options for the cross section shape are exactly the same as for the subcollector.
Finally, the outflow from the collector enters a main channel (Figure 113). The
channel is intended to represent the main stream in the subbasin. Outflow from the
collector is applied to the channel as lateral inflow. Optionally, upstream inflow can
be connected to the channel.
Normally it is not possible to connect upstream inflow to a subbasin element.
However, because of the conceptual representation of the kinematic wave method, it
is possible to have an upstream inflow on the main channel. In order to connect
upstream elements to the subbasin, you must select the route upstream option.
The surface runoff on the two planes is always routed using the kinematic wave
method. However, you have the choice of using Muskingum-Cunge routing in the
subcollector, collector, and main channel. All three channels use the same method.
Select the method you wish to use.
Figure 113.
If the route upstream option is turned on, the length should be the distance from the
upstream boundary to the subbasin to the outlet. If the route upstream option is
turned off, the length should be from the identifiable concentration point of the
channel to the subbasin outlet. There is only one channel is each subbasin.
The slope should be the average slope for the whole channel. 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.
The number of subreaches is used as a hint to the program when it determines the
correct distance step to use during routing calculations. The default value is five.
This parameter is only entered when the kinematic wave routing method is selected.
A cross section shape must be selected. If the kinematic wave routing method is
selected, the options will include circle, deep, rectangle, trapezoid, and triangle.
When the Muskingum-Cunge routing method is selected, the options will include
144
circular, rectangle, trapezoid, triangle, and eight point. You will have to enter a
Manning's n roughness value for the channel. If the eight point cross section shape
is used, you will also have to enter a roughness value for the left and right overbank
areas. You will have to enter other appropriate parameter data based on the cross
section choice.
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. The
Component Editor is shown in Figure 114.
The time of concentration defines the maximum travel time in the subbasin. The grid
cell in the subbasin with the largest travel time index will have exactly this specified
time of concentration. All other grid cells will have a scaled time of concentration
based on the ratio of the cell's travel time index to the maximum travel time index.
The storage coefficient is used in the linear reservoir for each grid cell. The same
coefficient is used for all cells in the subbasin.
Figure 114.
145
Natural Resources Conservation Service) found that in general the lag time can be
approximated by taking 60% of the time of concentration.
Figure 115.
Figure 116.
146
Figure 117.
147
Figure 118.
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, etc. This parameter provides a way to account for the effects of
urbanization on the hydrologic response of the subbasin.
148
The time lag is defined as the length of time between the centroid of precipitation
mass and the peak flow of the resulting hydrograph.
Figure 119.
Figure 120.
149
access the Component Editor by clicking on the element icon in the basin map, if the
map is currently open. You can select a baseflow method from the list of five
available choices. If you choose the None method, the subbasin will not compute
baseflow and the outflow will only include direct runoff from the transform method.
Use the selection list to choose the method you wish to use. Each subbasin may use
a different method or several subbasins may use the same method.
When a new subbasin is created, it is automatically set to use the default baseflow
method specified in the project settings. You may change the baseflow method for a
subbasin at any time using the Component Editor for the subbasin element. Since a
subbasin can only use one baseflow 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 baseflow method for several
subbasins simultaneously. Click on the Parameters menu and select the Subbasin
Methods Baseflow command. The baseflow method you choose will be applied
to the selected subbasins in the basin model, or to all subbasins if none are currently
selected.
The parameters for each baseflow method are presented on a separate Component
Editor from the subbasin element editor. The "Baseflow" editor is always shown near
the "Subbasin" editor. The information shown on the baseflow editor will depend on
which method is currently selected.
Figure 121.
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
150
per area per time. This method is better suited when general guidelines for
watershed yield must be used to estimate the initial flow.
The recession constant describes the rate at which baseflow recedes between storm
events. It is defined as the ratio of baseflow at the current time, to the baseflow one
day earlier.
A baseflow value must be entered for the month of January. Likewise, a value must
also be entered for each of the remaining months from February to December. This
value is used to limit the computed baseflow.
Figure 122.
151
Figure 123.
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. The same method must be
used for specifying the initial condition for both layers.
The groundwater storage coefficient is the time constant for the linear reservoir in
each layer. Since it is measured in hours, it gives a sense of the response time of
the subbasin.
The number of groundwater reservoirs can be used so that baseflow is routing
through several sequential reservoirs. Minimum attenuation is achieved when only
one routing step is selected. Attenuation of the baseflow increases as the number of
steps increases.
152
Figure 124.
There are two different methods for determining how to reset the baseflow during a
storm event: ratio to peak and threshold flow. When using the ratio to peak method,
you must specify the flow ratio to the peak. The baseflow is reset when the current
flow divided by the peak flow is falls to the specified value. For example, if a ratio of
0.2 is selected, the baseflow will be reset on the receding limb of an event
hydrograph when the flow has decreased to 20% of the event peak flow. With the
threshold flow method, the baseflow is always reset when the receding limb of the
hydrograph falls to a specified flow value, regardless of the peak flow during the
previous storm event.
The characteristic subsurface flow length must be specified. This could be estimated
as the mean distance from the subbasin boundary to the stream.
The conductivity of the soil must be specified. This could be estimated from field
tests of from the soil texture.
The drainable porosity must be specified in terms of volume ratio. The upper limit
would be the total porosity minus the residual porosity. The actual drainable porosity
depends on local conditions.
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. The Component Editor is shown in Figure 125.
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.
153
Figure 125.
The recession constant describes the rate at which baseflow recedes between storm
events. It is defined as the ratio of baseflow at the current time, to the baseflow one
day earlier.
There are two different methods for determining how to reset the baseflow during a
storm event: ratio to peak and threshold flow. When using the ratio to peak method,
you must specify the flow ratio to the peak. The baseflow is reset when the current
flow divided by the peak flow is falls to the specified value. For example, if a ratio of
0.2 is selected, the baseflow will be reset on the receding limb of an event
hydrograph when the flow has decreased to 20% of the event peak flow. With the
threshold flow method, the baseflow is always reset when the receding limb of the
hydrograph falls to a specified flow value, regardless of the peak flow during the
previous storm event.
154
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.
155
Figure 126.
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.
Figure 127.
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.
156
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.
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 128.
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.
Figure 128.
157
Figure 129.
There are two options for initial condition: specified discharge, and inflow equals
outflow. If you use the first option, you will also have to enter a discharge value. The
initial storage in the reach will be calculated from the specified discharge and the
storage-discharge function. If you use the second option, it will be assumed that the
initial outflow is the same as the initial inflow to the reach from upstream elements.
This is essentially the same as assuming a steady-state initial condition. The initial
storage will be computed from the first inflow to the reach and storage-discharge
function.
Optionally an elevation-discharge function can be selected. If used, it should
represent the depth of water for any given outflow from the reach. The function must
be monotonically increasing. The elevation-discharge function must be defined in the
paired data manager before it can be used in the reach element.
If the optional elevation-discharge function is selected, then an invert elevation
should also be specified. The flow depth is added to the invert elevation to compute
the stage.
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 130.
158
The Muskingum K is essentially the travel time through the reach. It can be
estimated from knowledge of the cross section properties and flow properties. It may
be a calibration parameter in some cases.
The Muskingum X is the weighting between inflow and outflow influence; it ranges
from 0.0 up to 0.5. In practical application, a value of 0.0 results in maximum
attenuation and 0.5 results in no attenuation. Most stream reaches require an
intermediate value found through calibration.
The number of subreaches must be entered. It affects attenuation where one
subreach gives more attenuation and increasing the number of subreaches
decreases the attenuation. A good approximation is to divide the reach length by the
product of the wave celerity and the simulation time step.
Figure 130.
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 131.
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.
159
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 slope (trapezoid,
triangle). The side slope is dimensionless and entered as the units of horizontal
distance per one unit of vertical distance. The eight point shape requires a cross
section simplified with only eight station-elevation values. The cross section is
usually configured to represent the main channel plus left and right overbank areas.
A separate Manning's n value is entered for each overbank. The cross section
should extend from the channel invert up to the maximum water surface elevation
that will be encountered during a simulation. The cross section must be created in
the Paired Data Manager before it can be used in the reach.
Figure 131.
160
Figure 132.
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
161
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 133.
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 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. If no
value is specified, then a value of zero is assumed.
Figure 133.
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 134.
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.
162
Figure 134.
163
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 135).
Figure 135.
165
method is selected. Some parameters are common to more than one routing
method. However, all possible required and optional data for each routing method is
described in the following sections.
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 storagedischarge 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.
Interpolation is used when the Elevation-Storage-Discharge or ElevationArea-Discharge storage methods are used. This means that it is not necessary
for the storage-discharge and elevation-storage curves used in the ElevationStorage-Discharge method to contain matching independent variables. The two
curves do not need to have the same storage values in each curve, or even have the
same number of rows. At compute time, the two curves selected by the user are
combined into a single routing table with three rows: elevation, storage, and
discharge. The table is initially configured using the curve selected by the user as
the Primary curve. The remaining column is interpolated from the curve not
selected as the primary curve. Finally the storage routing is completed from the
combined table using the storage and outflow columns, and then elevation and area
is calculated from the computed storages where possible. A similar procedure is also
used with the Elevation-Area-Discharge storage method.
You must select appropriate functions 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. If you select the Elevation-Storage-Discharge
method, you must select an elevation-storage and a storage-discharge paired data
function. Appropriate selection lists will be shown directly under the storage method
selection list. Any 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
166
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.
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 do depend on the selected
storage method and are shown in Table 21.
Table 21.
Storage Method
Storage-Discharge
Elevation-Storage-Discharge
Elevation-Area-Discharge
Storage Method
There are two different options for specifying the storage relationship, as shown in
Figure 136. 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.
167
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.
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 136.
168
Storage Method
There are two different options for specifying the storage relationship. The first
option is the Elevation-Storage choice, as shown in Figure 137. 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.
169
Figure 137.
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. Any 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.
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 elevationstorage 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
170
user is used to interpolate a storage value from the elevation-storage curve. The
initial condition options do depend on the selected storage method and are shown in
Table 22.
Table 22.
Storage Method
Elevation-Storage
Elevation-Area
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.
The Assume None method is used in cases where reservoir tailwater has no affect
on the reservoir outflow.
The Reservoir Main Discharge method is typically used with reservoirs that
span the stream channel and are not influenced by backwater from downstream
sources. For such cases, the tailwater below the reservoir only comes from the
reservoir releases. A rating curve defined by an elevation-discharge paired data
function must be selected to convert reservoir outflow to stage. The elevationdischarge function must be defined before it can be used in the reservoir. Choose an
appropriate function in the selection list or use a chooser by 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.
The Downstream Of Main Discharge method is typically used with reservoirs
that represent an interior pond or pump station, and the outflow from the reservoir will
be a significant impact on the downstream stage. In this case, the outflow from the
reservoir is combined with all other inflows to the element downstream of the
reservoir. That combined inflow is used in combination with a rating curve to
determine the stage for the reservoir tailwater. The elevation-discharge function for
the rating curve must be defined in the paired data manager before it can be used in
the reservoir. Choose an appropriate function in the selection list or use a chooser
by 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.
The Specified Stage method is typically used with reservoirs that represent an
interior pond or pump station, and the outflow from the reservoir will have minimal
affect on the downstream stage. In this case, the outflow from the reservoir is
adjusted for submergence based on the stage specified in a stage time-series gage.
The gage must be defined in the time-series data manager before it can be used in
the reservoir. Choose an appropriate gage in the selection list or use a chooser by
clicking the gage data button next to the selection list. The stage should be specified
in the same vertical datum as the function used to describe the storage
characteristics of the reservoir.
171
The Fixed Stage method is typically used with reservoirs that represent an interior
pond or pump station. The same stage is used for all time intervals in a simulation.
In this case, the outflow from the reservoir is adjusted for submergence based on the
specified stage.
Outlets
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 the 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 parameter data for it. There are two different methods for
computing outflow through an outlet: culvert or orifice.
172
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 138.
Figure 138.
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 select 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.
The shape specifies the cross-sectional shape of the culvert: circular, semi circular,
elliptical, arch, high-profile arch, low-profile arch, pipe arch, box, or con span. The
shape you choose will determine some of the remaining parameters in the
Component Editor. The parameters you will need to enter are shown in Figure 138.
The chart specifies the FHWA chart identification number. Only the charts that apply
to the selected shape will be shown in the selection list (Figure 138).
The scale specifies the FHWA scale identification number. Only the scales that apply
to the selected chart number will be shown in the selection list.
173
Table 23.
Listing of which parameters are required for each cross section shape.
Diameter
Rise
Span
Elliptical
Arch
Circular
Semi Circular
X
X
High-Profile Arch
Low-Profile Arch
Pipe Arch
Box
Con Span
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 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 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
139.
174
Figure 139.
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.
Spillways
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: broadcrested, 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 140.
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.
175
Figure 140.
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. Typical values range
from 2.6 to 4.0 depending on the exact shape of the spillway.
Ogee Spillway
The ogee spillway allows for controlled flow over the top of the reservoir according to
the weir flow assumptions. However, the discharge coefficient in the weir flow
equation is automatically adjusted when the upstream energy head is above or below
the design head. A typical ogee spillway editor is shown in Figure 141.
The ogee spillway may be specified with concrete or earthen abutments. These
abutments should be the dominant material at the sides of the spillway above the
crest. The selected material is used to adjust energy loss as water passes through
the spillway. The spillway may have one, two, or no abutments depending on how
the spillway or spillways in a reservoir are conceptually represented.
The ogee spillway is assumed to have an approach channel that moves water from
the main reservoir to the spillway. If there is such an approach channel, you must
specify the depth of the channel, and the energy loss that occurs between the main
reservoir and the spillway. If there is no approach channel, the depth should be the
different between the spillway crest and the bottom of the reservoir, and the loss
should be zero.
The crest elevation of the spillway must be entered. This should be measured in the
same vertical datum as the paired data functions defining the storage characteristics
of the reservoir.
The crest length of the spillway must be specified. This should be the total width
through which water passes.
The apron elevation is the elevation at the bottom of the ogee spillway structure.
This should be measured in the same vertical datum as the paired data functions
defining the storage characteristics of the reservoir.
176
Figure 141.
The apron width must be specified. This should be the total width of the spillway
bottom.
The design head is the total energy head for which the spillway is designed. The
discharge coefficient will be automatically calculated when the head on the spillway
departs from the design head.
Specified Spillway
The user-specified spillway can be used to represent spillways with flow
characteristics that cannot be represented by the broad-crested or ogee weir
assumptions. The user must create an elevation-discharge curve that represents the
spillway discharge as a function of reservoir pool elevation. At this time there is no
ability to include submergence effects on the spillway discharge. Therefore the userspecified spillway method should only be used for reservoirs with the downstream
tailwater stage cannot effect the discharge over the spillway. A typical user-specified
spillway editor is shown in Figure 142.
Figure 142.
177
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.
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 140 and Figure 141). 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 143.
The width of the sluice gate must be specified. It should be specified as the total
width of an individual gate.
The gate coefficient describes the energy losses as water passes under the gate.
Typical values are between 0.5 and 0.7 depending on the exact geometry and
configuration of the gate.
The orifice coefficient describes the energy losses as water passes under the gate
and the tailwater of the gate is sufficiently submerged. A typical value for the
coefficient is 0.8.
Figure 143.
178
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 144.
Figure 144.
The width of the radial gate must be specified. It should be specified as the total
width of an individual gate.
The gate coefficient describes the energy losses as water passes under the gate.
Typical values are between 0.5 and 0.7 depending on the exact geometry and
configuration of the gate.
The orifice coefficient describes the energy losses as water passes under the gate
and the tailwater of the gate is sufficiently submerged. A typical value for the
coefficient is 0.8.
The pivot point for the radial gate is known as the trunnion. The height of the
trunnion above the spillway must be entered.
The trunnion exponent is part of the specification of the geometry of the radial gate.
A typical value is 0.16.
The gate opening exponent is used in the calculation of flow under the gate. A
typical value is 0.72.
The head exponent is used in computing the total head on the radial gate. A typical
value is 0.62.
179
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 broadcrested 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 145.
Figure 145.
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.
180
The discharge coefficient accounts for energy losses as water approaches the dam
top and flows over the dam. Typical values range from 2.6 to 4.0 depending on the
exact shape of the dam top.
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 146.
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.
Figure 146.
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 for representation of 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
147.
181
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.
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 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.
Figure 147.
182
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 pumping capability 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 elevationdischarge 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 148) 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.
The bottom elevation defines the elevation of the bottom of the trapezoidal opening in
the dam face when the breach is fully developed. This should be measured in the
same vertical datum as the paired data functions defining the storage characteristics
of the reservoir.
The bottom width defines the width of the bottom of the trapezoidal opening in the
dam face when the breach is fully developed.
The left side slope is dimensionless and entered as the units of horizontal distance
per one unit of vertical distance. The right side slop is likewise dimensionless and
entered as the units of horizontal distance per one unit of vertical distance.
183
Figure 148.
The development time defines the total time for the breach to form, from initiation to
reaching the maximum breach size. It should be specified in hours.
There are three methods for triggering the initiation of the failure: elevation, duration
at elevation, and specific time. Depending on the method you choose, additional
parameters will be required. For the Elevation method, you will enter an elevation
when the failure should start. The breach will begin forming as soon as the reservoir
pool elevation reaches that specified elevation. For the Duration at Elevation
method, you will enter an elevation and duration to define when the failure should
start. The reservoir pool will have to remain at or above the specified elevation for
the specified duration before the 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.
184
Figure 149
185
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 potential 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 the reservoir, the seepage
is automatically taken from the reservoir storage and added to the main tailwater
discharge location. This is the mode of seepage when the pool elevation is greater
than the tailwater elevation. Seepage into the reservoir happens when the tailwater
elevation is higher than the pool elevation. In this mode the appropriate amount of
seepage is added to reservoir storage, but it is not subtracted from the tailwater.
Tabular Seepage
The tabular seepage method uses an elevation-discharge curve to represent
seepage as shown in Figure 150. Usually the elevation-discharge data will be
developed through a geotechnical investigation separate from the hydrologic study.
A curve may be specified for inflow seepage from the tailwater toward the pool, and a
separate curve selected for outflow seepage from the pool to the tailwater. The same
curve may be selected for both directions if appropriate. Any curve used for dam
seepage must first be created in the Paired Data Manager. If a curve is not selected
for one of the seepage directions, then no seepage will be calculated in that direction.
Figure 150.
Evaporation
Evaporation can only be included in reservoirs using the Outflow Structures
routing method, and that specify storage using the elevation-area 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.
186
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 151.
Figure 151.
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. 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 shown in Figure 152. The gage
must be defined in the Time-Series Data Manager before it can be selected.
187
Figure 152.
188
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 153).
Figure 153.
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 automatically compute drainage
area at downstream elements. If the area is not specified, then downstream
elements will not have a drainage area. When the source represents a region
modeled separately, then the area should be equal to the drainage area at the
source location in the watershed. When the source is used to represent a spring,
inter-basin transfer, or other movement of water, then the area should usually be set
to zero.
189
If you choose to include a representative area, specify the area on the Component
Editor for the source element (Figure 153). Access the Component Editor by clicking
the source element icon on the "Components" tab of the Watershed Explorer. You
can also access the Component Editor by clicking on the element icon in the basin
map, if the map is currently open.
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.
The time-series discharge gage must be defined in the time-series manager before it
can be used in the source element, as shown in Figure 154. Selecting the correct
gage is performed 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. You can also access the Component Editor by clicking on
the element icon in the basin map, if the map is currently open.
Figure 154.
190
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 155). You can also access the Component Editor by clicking on the
element icon in the basin map, if the map is currently open.
Figure 155.
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 156). You can also access the Component Editor by clicking on the
element icon in the basin map, if the map is currently open.
Figure 156.
191
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. Only one
method is available for computing the diverted flow: a user-specified, monotonically
increasing inflow-diversion relationship. All flow that is not diverted becomes main
outflow. Diverted outflow can be connected to an element that is computationally
downstream. The diversion can be used to represent weirs that divert flow into
canals, flumes, or off-stream storage. The 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
157).
Figure 157.
192
193
Select the inflow-diversion function from the list of available choices. The inflowdiversion function must be specified in the paired data manager before it can be
selected in the component editor.
Figure 158.
194
Figure 159.
195
Figure 160.
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 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 pumping capability 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 elevationdischarge functions in the project. Click on a function to view its description.
196
Figure 161.
The time-series discharge gage representing the diversion flows must be defined in
the time-series manager before it can be used in the diversion element, as shown in
Figure 161. Selecting the correct gage is performed on the Component Editor for the
diversion element. Access the Component Editor by clicking the diversion element
icon on the "Components" tab of the Watershed Explorer. You can also access the
Component Editor by clicking on the element icon in the basin map, if the map is
currently open.
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 162). You can also access the Component Editor by clicking on the element
icon in the basin map, if the map is currently open.
Figure 162.
197
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.
Meteorologic Models
Meteorologic models are one of the main components in a project. Their principle
purpose is to prepare meteorologic boundary conditions for subbasins.
Consequently, you must usually 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. To avoid this, you can rename the
subbasins so that separate boundary conditions are computed for each one.
199
Figure 163.
Figure 164.
The second way to copy is from the "Components" tab of the Watershed Explorer.
Move the mouse over the meteorologic model you wish to copy, then press the right
mouse button (Figure 165). A context menu is displayed that contains several
choices including copy. Click the Create Copy command. A new Copy
Meteorologic Model 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 meteorologic model. You cannot press the Copy
200
button if no name is specified. If you change your mind and do not want to copy the
selected meteorologic model, press the Cancel button or the X button in the upper
right of the Copy Meteorologic Model window to return to the Watershed Explorer.
Figure 165.
Figure 166.
201
The second way to rename is from the "Components" tab of the Watershed Explorer.
Select the meteorologic model you wish to rename by clicking on it in the Watershed
Explorer; it will become highlighted. Keep the mouse over the selected model and
click the right mouse button. Select the Rename command from the menu and the
highlighted name will change to editing mode as shown in Figure 167. 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 meteorologic model, press the Escape key.
Figure 167.
202
Figure 168.
Figure 169.
203
Precipitation Method
There are seven different precipitation methods or you can choose to have 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 170). 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.
Figure 170.
Some precipitation 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.
Also, some precipitation methods have optional settings. If the method you select
has optional settings then an additional tab may be added to the Component Editor.
Evapotranspiration Method
There are three different evapotranspiration methods or you can choose to have 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 170). 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.
204
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.
Snowmelt Method
There are two different snowmelt methods or you can choose to have 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. The precipitation method is responsible for estimating the
total depth of combined rain and snow. The snowmelt method determines if
precipitation is rain or snow, and tracks the accumulation and melt of the snowpack.
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 170). 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 timeseries 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 170). 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.
205
Figure 171.
206
Precipitation
Precipitation is one of the three components of a meteorologic model, and the one
required in almost all meteorologic models. The only time precipitation is not
required is if a basin model contains no subbasins. Any basin model that contains at
least one subbasin can only be used in a simulation with a meteorologic model that
includes a precipitation method. Some of the method choices are designed to work
with statistical data. Others are designed to work with gage data. Regardless of the
supporting data source, the precipitation method produces a hyetograph at the
correct time step for each subbasin.
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 from the
maps along with other information to compute a hyetograph for each subbasin.
This method uses the same parameter data for all subbasins in the meteorologic
model. The Watershed Explorer will show the precipitation icon one level under the
meteorologic model (Figure 172). Click on the icon to access the Component Editor
for the frequency storm (Figure 173).
Each storm has a single exceedance probability which must be selected from the list
of available choices. The choices range from 0.2 to 50 percent and generally match
the precipitation maps that are commonly available. Please note that return interval
is not used in the program.
The frequency storm method is designed to accept partial duration precipitation data,
as this is the most common type. However, annual duration data can also be entered
by selection the correct input option. Either annual or partial-duration output can be
calculated. You must select the type of input and output you wish the program to
compute. The difference between annual and partial-duration output is small for
exceedance probabilities 4% and smaller; you must select the type when the 10%,
20% or 50% exceedance probability is selected. The default selection is partial
duration input and annual duration output.
Figure 172.
207
Figure 173.
Precipitation editor for all subbasins when the frequency storm method
is selected.
The intensity duration specifies the shortest time period of the storm. Usually the
duration should be set equal to the time step of the simulation. It must be less than
the total storm duration.
Storm duration determines how long the precipitation will last. It must be longer than
the intensity duration. If the simulation duration is longer than the storm duration, all
time periods after the storm duration will have zero precipitation.
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, it only
changes 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 cannot be entered. The values must be entered
as partial-duration data. Values should be entered as the total precipitation depth
expected for the specified duration.
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
208
then the raw data are converted to a regular time step, such as 1 hour. Nonrecording 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.
This method uses separate parameter data for each gage used to compute
precipitation and also uses separate parameter data for each subbasin in the
meteorologic model. Because of the separate data and special options (Figure 174),
there are multiple tabs in the Component Editor for this method. The Watershed
Explorer may show a precipitation icon one level under the meteorologic model,
depending on options. It will always show a precipitation gages icon one level under
the meteorologic model. It will also always show a separate precipitation icon under
each subbasin icon, as shown in Figure 175. Click on any of the icons to access the
Component Editor for that portion of the gage weights method.
The gage weights precipitation icon is only shown one level under the meteorologic
model icon when the "Use Indexing" option is enabled. You can enable and disable
the indexing option from the "Options" tab of the Component Editor for the
meteorologic model (Figure 174). To access the options editor, click on the
meteorologic model icon in the Watershed Explorer and then click on the "Options"
tab in the Component Editor. If you decide to use indexing and turn on the option,
Figure 174.
Figure 175.
209
the gage weight precipitation icon will appear. Click on the icon to access the
Component Editor. The editor contains a table of all subbasins in the meteorologic
model with a separate row for each subbasin (Figure 176). For each subbasin you
can enter an index depth. The index depth is usually used to account for regional
trends in annual or monthly average precipitation. If you enter an index for a
subbasin, all precipitation gages used with that subbasin must also have a specified
index or else indexing is ignored for that subbasin.
Figure 176.
The precipitation gages icon is always shown one level under the meteorologic model
icon. When the Use Indexing option and Total Override option are both
turned off, the Component Editor for the precipitation gages only shows information
for total storm gages. If one or both of the options are turned on, the Component
Editor shows information for time-series and total storm gages. The "Time-Series
Gages" tab (Figure 177) shows all currently defined precipitation gages in the TimeSeries Data Manager. The "Total Storm Gages" tab (Figure 178) shows the gages
that have been defined in the meteorologic model.
Figure 177.
210
Entering optional total depth and index for the precipitation gages.
Figure 178.
Creating a total storm gage and entering the total storm depth and
optional index.
Total depth can be optionally entered for time-series gages. 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, the exact pattern is maintained but the
magnitude of precipitation at each time step is adjusted so that the specified depth is
applied over the entire simulation. Total depth can be specified for no time-series
gages, one gage, many gages, or all gages. Total depth is always required for a total
storm gage. Turning off the "Total Override" option will delete all total depth
information from the meteorologic model.
An index can be optionally entered for time-series gages or total storm gages. If you
enter an index for a gage, it can only be used during a simulation if all gages used in
a subbasin and the subbasin itself specify an index. Turning off the Use Indexing
option will delete 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. 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.
The parameter data describing which gages to use and what weights to apply is
specified separately for each subbasin. The Watershed Explorer will show all
subbasins in the meteorologic model and one level under each subbasin will be a
precipitation icon. Click on the icon to access the Component Editor for specifying
which gages to use and the weights for each gage.
The "Gage Selections" tab is where the gages are specified (Figure 179). All of the
available gages are shown in a table. The available gages are the precipitation timeseries gages defined in the Time-Series Data Manager plus any total storm gages
defined in the meteorologic model. For each subbasin you must separately select
which gages will be used for that subbasin.
211
Figure 179.
Figure 180.
The "Gage Weights" tab is where the weights are specified for each gage selected
for a subbasin (Figure 180). The gages are shown in a table with a separate row for
each gage. Only the gages selected previously on the "Gage Selections" tab are
included in the table. For each gage you can enter a depth weight. You can enter a
time weight for time-series gages. The values entered for the depth or time weights
are automatically normalized during the simulation. The value of the weights must be
computed separate from the program. Possible methods for computing the weights
include Thiessen polygons, inverse distance, inverse distance squared, isohyetal
mapping, or any method deemed appropriate.
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 for simulation. 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.
This method uses the same parameter data for all subbasins in the meteorologic
model. The Watershed Explorer will show the precipitation icon one level under the
meteorologic model. Click on the icon to access the Component Editor for gridded
precipitation data (Figure 181).
212
Figure 181.
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. Optionally, an index
depth can be assigned to each gage. The index is used to adjust for regional bias in
annual or monthly precipitation.
This method uses separate parameter data for each gage used to compute
precipitation and also uses separate parameter data for each subbasin in the
meteorologic model. Because of the separate data and special options (Figure 182)
there are multiple tabs in the Component Editor for this method. The Watershed
Explorer will show a precipitation gages icon one level under the meteorologic model.
It will also show a separate precipitation icon under each subbasin icon (Figure 183).
Click on any of the icons to access the Component Editor for that portion of the
inverse distance method.
213
Figure 182.
Figure 183.
The precipitation gages icon is always shown one level under the meteorologic model
icon. The Component Editor shows a table of all precipitation gages from the TimeSeries Data Manager, with a separate row for each gage. In the first column of the
table you must select if you wish to use each gage in this meteorologic model as
shown in Figure 184. The latitude and longitude location must be specified for each
gage. Enter the location information using the Component Editor for precipitation
gages in the time-series data section of the Watershed Explorer.
The second column of the Component Editor is used to select how the gage data is
processed. All gages used in the inverse distance method must be setup in the
Time-Series Data Manager before they can be used. Each 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 24hour 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.
You can enable and disable the indexing option from the "Options" tab of the
Component Editor for the meteorologic model. To access the options editor, click on
the meteorologic model icon in the Watershed Explorer and then click on the
"Options" tab in the Component Editor. If you decide to use indexing and turn on the
option, the extra column will be added to the Component Editor. For each gage you
can enter an index depth. The index depth is usually used to account for regional
214
Figure 184.
trends in annual or monthly average precipitation. If you enter an index for a gage,
all precipitation gages and subbasin nodes must also have a specified index or else
indexing is ignored.
You must create one or more nodes for each subbasin using the inverse distance
precipitation icon under each subbasin in the meteorologic model (Figure 185). The
Component Editor contains three tabs. The first "Node Weights" tab is used to
specify nodes for a subbasin, with a minimum of one required node. The remaining
tabs are used to enter the latitude and longitude of the nodes.
Nodes are created and managed directly from the Component Editor. 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. 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.
You can control how far the program searches from each subbasin node to find
precipitation gages. Specify the maximum distance to search from the "Options" tab
of the Component Editor for the meteorologic model. To access the options editor,
click on the meteorologic model icon in the Watershed Explorer and then click on the
"Options" tab in the Component Editor.
Figure 185.
215
SCS Storm
The SCS storm method implements the design storm developed by the Natural
Resources Conservation Service (formerly the Soil Conservation Service). 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.
This method uses the same parameter data for all subbasins in the meteorologic
model. The Watershed Explorer will show the precipitation icon one level under the
meteorologic model. Click on the icon to access the Component Editor for the SCS
storm (Figure 186).
Each storm has only one time distribution type which must be selected from the list of
available choices. The available types are Type 1, Type 1A, Type 2, and Type 3.
Each storm is 24 hours long. The simulation must have a duration at least 24 hours
long. All precipitation values after the first 24 hours will be zero.
Each storm has the same total precipitation depth for each subbasin in the
meteorologic model. There is no depth-area reduction or any other kind of reduction
or modification performed on the specified depth. The specified depth is applied to
the selected time distribution type to generate the hyetograph for each subbasin.
Figure 186.
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 also useful when a single precipitation gage can be used
to represent what happens over a subbasin. Several options are available to
increase control over how the data is processed (Figure 187).
Figure 187.
216
This method uses separate parameter data for each subbasin in the meteorologic
model, but all subbasins are shown in a single table. Each subbasin is shown on a
separate row of the table. The Watershed Explorer will show the precipitation icon
one level under the meteorologic model. Click on the icon to access the Component
Editor for the specified hyetograph (Figure 188).
A hyetograph must be stored as a precipitation gage before it can be used in the
meteorologic model. The data may actually be from a recording gage or could be the
result of complex calculations exterior to the program. Regardless, the hyetograph
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
precipitation gages already defined will be shown in the selection list.
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 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.
Figure 188.
217
The distribution determines how the adjusted precipitation depth for each subbasin is
shaped into a hyetograph; it must be selected from the list of available choices. The
Standard option uses the procedure specified in the Engineering Manual to
distribute the adjusted storm depth. The Southwest Division option uses a
different procedure that may be more applicable in some watersheds.
The transposition factor accounts for the location of a subbasin in the watershed,
relative to the center of the storm (Figure 190). The Engineering Manual contains an
isohyetal map with concentric rings labeled in percent. If a subbasin is generally
covered by the 120% isohyetal line, then a factor of 1.20 should be entered in the
program.
Figure 189.
Figure 190.
Evapotranspiration
Evapotranspiration is the second of the three components of a meteorologic model.
It is the combination of evaporation from the ground surface and transpiration by
vegetation. It is only required if the meteorologic model will be used with subbasins
that use continuous simulation loss rate methods: deficit constant, gridded deficit
constant, soil moisture accounting, and gridded soil moisture accounting. Even if
those methods are used in the subbasins, an evapotranspiration method is not
required. If a continuous simulation loss rate method is used and no
evapotranspiration is specified in the meteorologic model, then zero
evapotranspiration is used in the subbasins. However, evapotranspiration is often
responsible for returning 50 or even 60% of precipitation back to the atmosphere. In
general there should be a selected evapotranspiration method for continuous
simulation. In all cases, the meteorologic model is computing the potential
evapotranspiration and subbasins will calculate actual evapotranspiration based on
soil water limitations.
218
The gridded Priestley Taylor method is designed to work with the gridded ModClark
transform method. It is the same as the regular Priestley Taylor (described in a later
section) method except the Priestley-Taylor equation is applied to each grid cell using
separate boundary conditions instead of area-averaged values over the whole
subbasin. If it is used with a transform method other than ModClark, an areaweighted average of the grid cells in the subbasin is used to compute the boundary
condition of the subbasin.
This method uses the same parameter data for all subbasins in the meteorologic
model. The Watershed Explorer will show the evapotranspiration icon one level
under the meteorologic model. Click on the icon to access the Component Editor for
the gridded Priestley Taylor (Figure 191).
A crop coefficient gridset must be selected for all subbasins. Only crop coefficient
gridsets already defined will be shown in the selection list. The crop coefficient is
applied on a grid cell basis to the evapotranspiration computed with the PriestleyTaylor method to give the final potential evapotranspiration. The coefficient can be
used to adjust the potential demand for different plants, growing season, and other
physiological factors.
A dryness coefficient must be entered for all subbasins. The same coefficient is
applied to all grid cells in all subbasins. The coefficient is used to make small
corrections based on soil moisture state. A coefficient should be specified that
represents typical soil water conditions during the simulation. A value of 1.2 can be
used in humid conditions while a value of 1.3 represents an arid environment.
A solar radiation gridset must be selected for all subbasins. The available gridsets
from the Grid Data Manager will be shown in the selection list. If no gridsets are
currently available, you will not be able to make a selection. Gridsets must be
created in the manager before they can be used in the meteorologic model.
A temperature gridset must be selected for all subbasins. The current gridsets are
shown in the selection list. If there are many different gridsets available, you may
wish to choose a gridset from the selector accessed with the grid button next to the
selection list. The selector displays the description for each gridset, making it easier
to select the correct one.
Figure 191.
219
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.
This method uses separate parameter data for each subbasin in the meteorologic
model. The Watershed Explorer will show an evapotranspiration icon under each
subbasin icon in the meteorologic model. Click on the evapotranspiration icon to
access the Component Editor for monthly average (Figure 192). The editor contains
a table with a separate row for each month.
The evapotranspiration rate must be entered for each month. It is entered as the
total amount of evapotranspiration for the month. Every time step within the month
will have the same evapotranspiration rate.
The pan coefficient must also be entered for each month. The specified rate is
multiplied by the coefficient to determine the final potential rate for each month. The
coefficient is usually used to correct actual evaporation pan data to more closely
reflect plant water use.
Figure 192.
Priestley Taylor
The Priestley Taylor method implements the Priestley-Taylor equation for computing
evapotranspiration. It is capable of capturing diurnal variation in potential
evapotranspiration through the use of a solar radiation gage, so long as the
simulation time step is less than 24 hours. A crop coefficient is used to adjust the
calculated Priestley-Taylor amount.
This method uses separate parameter data for each subbasin in the meteorologic
model. The Watershed Explorer will show an evapotranspiration icon under each
subbasin icon in the meteorologic model. Click on the evapotranspiration icon to
access the Component Editor for Priestley Taylor (Figure 193).
A crop coefficient gage must be selected for a subbasin. You may use the same
gage for more than one subbasin. Only crop coefficient gages already defined will be
220
shown in the selection list. The crop coefficient is applied to the evapotranspiration
computed with the Priestley-Taylor method to give the final potential
evapotranspiration. The coefficient can be used to adjust the potential demand for
different plants, growing season, and other physiological factors.
A dryness coefficient must be entered for a subbasin. The coefficient is used to
make small corrections based on soil moisture state. A coefficient should be
specified that represents typical soil water conditions during the simulation.
A solar radiation gage must be selected for a subbasin. The available gages from
the Time-Series Data Manager will be shown in the selection list. If no gages are
currently available, you will not be able to make a selection. Gages must be created
in the manager before they can be used in the meteorologic model.
A temperature gage must be selected for a subbasin. The current gages are shown
in the selection list. 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.
Figure 193.
Snowmelt
Snowmelt is the third of the three components of a meteorologic model. It uses the
precipitation computed by the precipitation method. Using temperature, it determines
whether the precipitation previously computed was liquid rain or frozen snow. The
accumulation and melt of the snowpack is simulated in response to atmospheric
conditions. The output of the method is the liquid water available at the soil surface,
which becomes the hyetograph for the subbasin.
221
Figure 194.
This method uses the same parameter data for all subbasins in the meteorologic
model. The Watershed Explorer will show the snowmelt icon one level under the
meteorologic model. Click on the icon to access the Component Editor for gridded
temperature index (Figure 194). A temperature gridset must be selected as one of
the atmospheric boundary conditions on the model. The current gridsets are shown
in the selection list. If there are many different gridsets available, you may wish to
choose a gridset from the selector accessed with the grid button next to the selection
list. The selector displays the description for each gridset, making it easier to select
the correct one.
The initial values can be specified with gridsets, or can be initialized to default values.
All cells in each subbasin are set to the same default value when the Default
Values option is choosen. The default values are: initial SWE of 0mm, initial cold
content of 0mm, initial liquid of 0%, initial cold content ATI of 0C, and initial meltrate
ATI of 0DegC-day. Five parameter grids must be selected when the Parameter
Grids option is selected; the grids are described in the following paragraphs.
The initial snow water equivalent that exists at the beginning of the simulation must
be entered. This information is usually determined by interpolating from actual
measurements of snow water equivalent. This value can be set to zero if there is no
snow.
222
The initial cold content that exists at the beginning of the simulation must be entered.
It represents the heat required to raise the temperature of the snow pack to 0C (32F)
and is expressed as a number equivalent to mm (inches) of frozen water. Generally,
this value is not known at the start of simulation unless there is no snow, in which
case it can be set to zero. If the value is not known, it can be set to zero. The error
in doing this may be small for relatively shallow ephemeral snow covers but may
cause errors for deep, seasonal snowpacks.
The liquid water held within the snowpack at the beginning of the simulation must be
entered. For any melt or precipitation to get through the snowpack, the liquid water
holding capacity of the snow must first be satisfied. Liquid water can persist in the
snow only if the snowpack temperature is at 0C (32F); at which point the cold content
is zero. A snowpack with liquid water is said to be ripe. Generally this value is not
known at the start of the simulation unless there is no snow, in which case it can be
set to zero. If the snow is known to be colder than 0C (32F) at the start of the
simulation, this value can be set to zero.
The initial cold content antecedent temperature index is an index to the snow
temperature near the surface of the snowpack. It is calculated assuming an
approximation to the transient heat flow equations. This value is used to estimate the
cold content of the snow. It should be set to the approximate snowpack temperature
at the beginning of the simulation. If the initial temperature is not known, it can be set
to 0C (32F).
The seasonal variation of meltrate is indexed by an antecedent temperature function.
The initial meltrate ATI should be thought of as similar to the accumulated thawing
degree days. This antecedent temperature function allows the melt rate to change
as the snowpack matures and ages. If there is no snow on the ground at the start of
the simulation this term can be set to zero. It can also be set to zero if the simulation
is starting during or at the end of a cold period when air temperatures were
continually below the base temperature.
The PX temperature is used to discriminate between precipitation falling as rain or
snow. When the air temperature is less than the specified temperature, any
precipitation is assumed to be snow. When the air temperature is above the
specified temperature, any precipitation is assumed to be rain. This discrimination
temperature is usually one to two degrees above freezing.
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 assumed to be zero. Typically, the base temperature should be 0C
(32F) or close to it.
The wet meltrate is used during time periods of precipitation when the precipitation is
falling as rain, at rates greater than the rain rate limit. It represents the rate at which
the snowpack melts when it is raining on the pack.
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.
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
223
from one time interval to the next. This index generally ranges from 0.015 to 0.550
and is separate from the cold content index. A typical value for the coefficient is 0.98.
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.
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. A typical value is 20 mm/day (0.8 in/day).
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 typical value for the
coefficient is 0.84.
An antecedent temperature index cold content function is used to calculate a cold
content from the current cold content index. The index typically ranges from 0.010 to
0.025. 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 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 (32F);
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.
224
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
degree above freezing. As the snowpack internal conditions and atmospheric
conditions change, the melt coefficient also changes.
Some of the parameter data required by this method is the same for all subbasins in
the meteorologic model. The Watershed Explorer will show the snowmelt icon one
level under the meteorologic model. Click on the icon to access the Component
Editor for component data that is the same for all subbasins (Figure 195).
Figure 195.
225
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.
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. This index generally ranges from 0.015 to 0.550
and is separate from the cold content index. A typical value for the coefficient is 0.98.
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 it 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.
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 it 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, cold content index is computed as an
antecedent index. A typical value is 20 mm/day (0.8 in/day).
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 typical value for the
coefficient is 0.84.
An antecedent temperature index cold content function is used to calculate a cold
content from the current cold content index. The index typically ranges from 0.010 to
0.025. 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 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 (32F);
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
226
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.
Some of the parameter data required by this method is different for each subbasin
included in the meteorologic model. The Watershed Explorer will also show an icon
for each subbasin; the snowmelt icon is shown under each subbasin. Click on the
snowmelt icon under the subbasin to access the Component Editor for component
data that is unique for each subbasin (Figure 196).
Figure 196.
227
Figure 197.
Figure 198.
You must specify the percentage of the subbasin that each elevation band
composes. An elevation band is not required to be contiguous. The percentage
specified for each elevation band will automatically be normalized if the sum of the
percentages across all subbasins does not equal 100. There is no limit to the
number of elevation bands that can be used, but at least one is required. Typically
only one band is used in watersheds with small elevation differences. Mountainous
watersheds usually require several bands for each subbasin.
Enter the average elevation for each elevation band. Typically the specified
elevation will be either the area-weighted elevation of the band, or the average of the
highest and lowest points in the band.
The initial snow water equivalent that exists at the beginning of the simulation must
be entered. This information is usually determined by interpolating from actual
measurements of snow water equivalent. This value can be set to zero if there is no
snow.
228
The initial cold content that exists at the beginning of the simulation must be entered.
It represents the heat required to raise the temperature of the snow pack to 0C (32F)
and is expressed as a number equivalent to mm (inches) of frozen water. Generally,
this value is not known at the start of simulation unless there is no snow, in which
case it can be set to zero. If the value is not known, it can be set to zero. The error
in doing this may be small for relatively shallow ephemeral snow covers but may
cause errors for deep, seasonal snowpacks.
The liquid Water held within the snowpack at the beginning of the simulation must be
entered. For any melt or precipitation to get through the snowpack, the liquid water
holding capacity of the snow must first be satisfied. Liquid water can persist in the
snow only if the snowpack temperature is at 0C (32F); at which point the cold content
is zero. A snowpack with liquid water is said to be ripe. Generally this value is not
known at the start of the simulation unless there is no snow, in which case it can be
set to zero. If the snow is known to be cold at the start of the simulation, this value
can be set to zero.
The initial cold content antecedent temperature index is an index to the snow
temperature near the surface of the snowpack. It is calculated assuming an
approximation to the transient heat flow equations. This value is used to estimate the
cold content of the snow. It should be set to the approximate snowpack temperature
at the beginning of the simulation. If the initial temperature is not known, it can be set
to 0C (32F).
The seasonal variation of meltrate is indexed by an antecedent temperature function.
The initial meltrate ATI should be thought of as similar to the accumulated thawing
degree days. This antecedent temperature function allows the melt rate to change
as the snowpack matures and ages. If there is no snow on the ground at the start of
the simulation this term can be set to 0 DEG C-DAY (0 DEG F-DAY). It can also be
set to 0 if the simulation is starting during or at the end of a cold period when air
temperatures were continually below the base temperature.
References
Engineering Manual 1110-2-1411. March 1065. Standard Project Flood
Determination. U.S. Army Corps of Engineers, Davis, CA.
229
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. A run may form the basis for additional analysis using
optimization trials or analyses. In addition to configuring and computing simulation
runs, this chapter also describes how to access the computed results.
Simulation Runs
Simulation runs are one of the three different components that can compute results:
simulation runs, optimization trials, and analyses. 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.
231
Figure 199.
The third step of creating a simulation run is to select a meteorologic model. All of
the meteorologic models in the project are shown and you must select one before
proceeding to the next step. By default the first meteorologic model in the table is
selected. The selected model is highlighted. You can use your mouse to select a
different meteorologic 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 meteorologic 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 select a different basin model.
Figure 200.
232
Entering a name for a new simulation run. The remaining steps are to
select a basin model, meteorologic model, and control specifications.
The fourth and final step of creating a simulation run is to select a control
specifications. All of the control specifications in the project are shown and you
must select one before proceeding to the next step. By default the first control
specifications in the table is selected. The selected specifications is highlighted. You
can use your mouse to select a different control specifications by clicking on it in the
table of available choices. You can also use the arrow keys on your keyboard to
select a different specifications. Press the Finish button when you are satisfied with
the name you have entered and the components you have selected, and are ready to
create the simulation run. Press the <Back button if you wish to return to the
previous step and select a different meteorologic model.
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 current 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 201) 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 to return to the Simulation Run
Manager window.
Figure 201.
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 202). 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.
233
Figure 202.
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 current 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 203) 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.
Figure 203.
The second way to rename is from the "Compute" tab of the Watershed Explorer.
Select the simulation run you wish to rename by clicking on it in the Watershed
Explorer; it will become highlighted. Keep the mouse over the selected run and click
the left mouse button again. The highlighted name will change to editing mode
(Figure 204). 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
234
changing the name, press the Enter key to finalize your choice. You can also finalize
your choice by clicking elsewhere on the "Compute" tab. If you change your mind
while in editing mode and do not want to rename the selected simulation run, press
the Escape key.
Figure 204.
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 you will have to
go to those components and manually select a different simulation run.
The first way to perform a deletion is to use the Simulation Run Manager, which you
can access from the Compute menu. Select the simulation run you wish to delete by
clicking on it in the list of current simulation runs. The selected run is highlighted
after you select it. After you select a run you can press the Delete button on the right
side of the window. A window will open where you must confirm that you wish to
delete the selected run as shown in Figure 205. 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 Simulation
Run Manager window.
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 206). 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.
235
Figure 205.
Figure 206.
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 207). You will then be
provided a selection list of available simulation runs in that source project (Figure
208). Select the simulation run that you wish to import and press the Import button.
236
Figure 207.
Figure 208.
Selecting the simulation run to be imported from the source project into
the current project.
Depending on the data used in the basin model and meteorologic model, you may be
prompted to choose how to treat time-series, paired, and gridded data. If you select
to import the data, a copy of the data will be made in the current project directory and
appropriate gage, paired data, and grid components will be automatically created.
Any references in the basin model or meteorologic model to the time-series, paired,
and gridded data will be automatically updated to use the new shared component
data in the current project. If you choose not to import the data, the basin model and
237
meteorologic model will be imported without any of the data. In this second case, the
current project should already include the necessary time-series, paired, and gridded
data necessary for the proper simulation of the basin model and meteorologic model.
Depending on the data used in the basin model and meteorologic model, you may be
prompted to choose how to treat external data. If you choose to localize external
data, then copies of the files will be made in the current project directory. Any
references in the basin model or meteorologic model to external data will be
automatically updated to use the new copies in the current project. If you choose not
to localize the external data, the components in the current project will be updated to
reference the files in their original locations.
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 209). 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.
Figure 209.
238
specified in the basin model. Optionally, you can apply a ratio to the precipitation
computed by the meteorologic model before the precipitation is applied to the basin
model. Alternately, you may instead apply a ratio to the outflow computed by
subbasin and source elements in the basin model before routing the outflow
downstream through the element network. You must choose between applying no
ratio, a precipitation ratio, or a flow ratio. The same ratio is applied to all elements.
The ratio for the simulation run is accessed from the "Compute" tab of the Watershed
Explorer. Click on the "Ratio" node under the simulation run to display the
Component Editor for the ratio (Figure 210). The "Ratio" and other optional tabs are
always shown when ever the Component Editor is shown. Select the type of ratio
you wish to apply from the available choices. If you select the precipitation method,
you can only choose to apply the ratio to subbasins. If you select the flow method,
you can choose to apply the ratio to subbasins, sources, or both element types.
Finally enter the ratio value. When working with a precipitation ratio, the precipitation
value for each time step computed by the meteorologic model is multiplied by the
specified ratio before the precipitation is applied to the subbasin elements. When
working with a flow ratio, the outflow computed for each time step by the basin model
is multiplied by the specified ratio before the outflow is routed to the next downstream
element. You can only select whether to apply a ratio to subbasins and sources, and
enter a ratio, if a ratio method is actually selected.
Figure 210.
239
States have two main purposes: breaking long simulations into smaller time periods,
and real-time forecasting operations. In the first case of long simulations, using
states achieves exactly the same results as 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.
Figure 211.
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 times
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 211). The "Start
States" and other optional tabs are always shown when ever the Component Editor is
shown. You can choose the states you wish to use from the selection list. Only
states saved by a different simulation run using the same basin model are shown.
The only way to create states for use as saved states is by allowing the program to
save them during a separate simulation run. You will not be able to select start
states unless there are valid states available.
Selecting optional save states for the simulation run is accessed from the "Compute"
tab of the Watershed Explorer. Click on the "Save States" node under the simulation
run to display the Component Editor for save states (Figure 212). The "Save States"
and other optional tabs are always shown when ever the Component Editor is shown.
You must specify a name for the saves and optionally may enter a description. There
240
are two ways to specify when the states will be saved: at the end of the simulation
run, or at a given time during the simulation run. If you select to save states at the
end of the run, the program will automatically determine the correct date and time to
save states. If you select to save states at a specific time, you must enter the date
and time. The specified date and time must be after the start and either at or before
the end of the simulation run.
Figure 212.
Checking Parameters
The program is equipped with the capability to check parameter values in the basin
and meteorologic model as part of performing a simulation run. The check is
performed in two stages. The first stage of checking is performed before the actual
simulation run begins. During the first stage, required parameters are checked to
make sure they have been specified. They are also checked to make sure they are
between the minimum and maximum allowed values. Checking also makes certain
that the specified values are physically possible. In some cases, one required
parameter may have relationships with other required parameters. When possible,
the dependencies are checked to make sure they are valid. For example, the Clark
unit hydrograph is required to have a time of concentration. If one is not entered, the
first stage of parameter checking will detect that it is missing and issue an error
message. During the second stage, checking is performed that cannot happen until
the actual simulation begins. For example, all inflow to a reach must be greater than
or equal to zero. There is no way to determine if this condition is satisfied until all
operations necessary to compute inflow to that reach have finished. Second stage
parameter checking may result in additional messages.
241
All messages issued during first or second stage parameter checking are classified
as notes, warnings, or errors. Notes are used to communicate general information to
the user. They may also be used to indicate actions taken where the result was
within established guidelines. Warnings are used when the program must take action
in order to continue the simulation, and the actions taken may conflict with the
intentions of the user. These messages should be checked to make sure the
program made appropriate assumptions. Errors are used for problems that stop the
simulation from proceeding. These messages can be used to determine what
required parameters are missing, or what specified parameters have invalid values.
All messages are shown in the message window at the bottom of the program
screen, regardless of when or how they are produced.
Parameter checking is performed automatically as part of computing a simulation run.
However, you can manually perform the first stage parameter checking to aid in
preparing basin and meteorologic models for use. First be certain that a simulation
run is selected. To perform the check, click on the Compute menu and select the
Check Parameters command. Any messages generated during the check will be
displayed in the message window.
Computing a Run
There are three ways to compute the currently selected simulation run. All methods
automatically perform parameter checking, and if no errors are generated, proceed to
the actual computing of simulation results. Additionally, the program is designed to
be computationally efficient. Only components with data changes since the last
compute will be recomputed; the remainder of the components have not changed so
previous simulation results are still applicable. If you wish, you may force all
components to recomputed regardless of data changes since the last compute.
The first way to compute a simulation run is from the Compute menu. Click on the
Compute menu and then select the Compute Run command. The name of the
current simulation run is shown in brackets as part of the menu command. If the
command is not available, it is because there is no current simulation run; you must
first select a run. A window will automatically open that shows the progress of the
compute. You will need to manually close the window when the compute is done,
whether it failed or was successful. If you wish to force all components to be
recomputed instead of just those with data changes, hold the control key while
selecting the menu command.
The second way to compute a simulation run is from the Watershed Explorer, on the
"Compute" tab. Select the simulation run you wish to compute by clicking on it in the
Watershed Explorer with the right mouse button. A context menu is displayed that
contains several choices including compute. Click the Compute command. The
progress window will automatically open. If you wish to force all components to be
recomputed instead of just those with data changes, hold the control key while
selecting the right mouse menu command.
The third way to compute a simulation run is from the toolbar. The compute button is
enabled when ever there is a current run that can be computed. If the button is not
available or shows a different type of compute, you must first select a simulation run.
Press the button to compute the current simulation run. The progress window will
automatically open. If you wish to force all components to be recomputed instead of
just those with data changes, hold the control key while selecting the toolbar button.
242
Figure 213.
243
while the result is open for viewing. In this case, the open result will immediately
show an indication that data has changed and the simulation run needs to be
recomputed. After the run is recomputed the open results are automatically updated
with the new results and the indication is updated with the date and time of the mostrecent compute.
244
Figure 214.
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 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 215) varies by element type, but always includes outflow. Optional items
such as observed flow, computed stage, and observed stage are also included.
245
Figure 215.
The information included in the summary table (Figure 216) 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 217). 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.
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. 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 simulation run 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 simulation run do not need to be recomputed.
246
Figure 216.
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 simulation run do not need to
be recomputed, and there is an element selection in the basin map.
Figure 217.
247
Figure 218.
Individual Elements
Each element in the basin model is shown in the Watershed Explorer under the
simulation run node. These elements are listed in alphabetical order after 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 218). 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 accessed by clicking on the "Summary Table" or "Time-Series
Table" node, respectively.
248
Figure 219.
249
Figure 220.
After you have opened a time-series table or graph, you may add additional timeseries 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 timeseries. 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.
250
Figure 221.
251
CHAPTER 13
Parameter Estimation
This chapter describes how observed streamflow can be used to automatically
estimate parameters. A process called optimization is used that begins from 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 estimates to the final best 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 highly correlated with
parameters such as hydraulic conductivity. Finally, there are parameters that can
only be estimated by comparing computed results to observed results such as
observed 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 values best able to reproduce the
observed hydrograph are found. Constraints are set to insure that unreasonable
parameter values are not used.
Optimization trials are one of the three different components that can compute
results: simulation runs, optimization trials, and analyses. Each trial is based on a
simulation run. The run provides the basic framework of a basin model, meteorologic
model, and control specifications within which parameters are estimated. A variety of
result graphs and tables are available from the Watershed Explorer for evaluating the
quality of the estimation.
The iterative parameter estimation procedure used by the program is often called
optimization. Initial values for all parameters are required at the start of the
optimization trial window. A hydrograph is computed at a target element by
computing all of the upstream elements. The target must have an observed
253
hydrograph for the time period over which the objective function will be evaluated.
Only parameters for upstream elements can be estimated. The value of the objective
function is computed at the target element using the computed and observed
hydrographs. Parameter values are adjusted by the search method and the
hydrograph and objective function for the target element are recomputed. This
process is repeated until the value of the objective function is sufficiently small, or the
maximum number of iterations is exceeded. Results can be viewed after the
optimization trial is complete.
Figure 222.
The first step of creating an optimization trial is to provide the name for the new trial
(Figure 223). 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.
254
Figure 223.
Entering a name for a new optimization trial. The remaining steps are
to select a simulation run and the optimization location.
The second step of creating an optimization trial is to select a simulation run that will
form the basis of the trial. All of the simulation runs containing at least one element
with observed flow are shown and you must select one before proceeding to the next
step. By default the first simulation run in the table is selected. The selected run is
highlighted. You can use your mouse to select a different simulation run by clicking
on it in the table of available choices. You can also use the arrow keys on your
keyboard to select a different run. Press the Next> button when you are satisfied
with the simulation run 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 optimization trial.
The third and final step of creating an optimization trial is to select an element where
there is observed flow. The element is where the objective function will be evaluated.
You will only be able to perform parameter estimation at, or upstream of the selected
element. The elements in the selection list come from the basin model used in the
simulation run that was selected in the previous step. By default the first element in
the table is selected. The selected element is always highlighted. You can use your
mouse to select a different element by clicking on it in the table of available choices.
You can also use the arrow keys on your keyboard to select a different element.
Press the Finish button when you are satisfied with the name you have entered, the
simulation run and element you selected, and are ready to create the optimization
trial. Press the <Back button if you wish to return to the previous step and select a
simulation run.
Copying a Trial
There are two ways to copy an optimization trial. Both methods for copying a trial
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 Optimization Trial Manager, which is
accessed from the Compute menu. Select the optimization trial you wish to copy by
clicking on it in the list of current optimization trials. The selected trial is highlighted
255
after you select it. After you select a trial you can press the Copy button on the
right side of the window. A new Copy Optimization Trial window (Figure 224) 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 to return to the Optimization Trial
Manager window.
Figure 224.
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 225). A context menu is displayed that contains several choices
including copy. Click the Create Copy command. A new 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 225.
256
Renaming a Trial
There are two ways to rename an optimization trial. Both methods for renaming a
trial change its name and perform other necessary operations.
The first way to perform a rename is to use the Optimization Trial Manager, which
you can access from the Compute menu. Select the optimization trial you wish to
rename by clicking on it in the list of current optimization trials. The selected trial is
highlighted after you select it. After you select a trial you can press the Rename
button on the right side of the window. A new Rename Optimization Trial window
(Figure 226) 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 optimization trial. 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 trial, press the Cancel button or the X button in the upper
right of the Rename Optimization Trial window to return to the Optimization Trial
Manager window.
Figure 226.
The second way to rename is from the "Compute" tab of the Watershed Explorer.
Select the optimization trial you wish to rename by clicking on it in the Watershed
Explorer; it will become highlighted. Keep the mouse over the selected trial and click
the left mouse button again. The highlighted name will change to editing mode
(Figure 227). 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 "Compute" tab. If you change your mind
while in editing mode and do not want to rename the selected optimization trial, press
the Escape key.
Figure 227.
257
Deleting a Trial
There are two ways to delete an optimization trial. Both methods for deleting a trial
remove it from the project and automatically delete previously computed results.
Once a trial has been deleted it cannot be retrieved or undeleted.
The first way to perform a deletion is to use the Optimization Trial Manager, which
you can access from the Compute menu. Select the optimization trial you wish to
delete by clicking on it in the list of current optimization trials. The selected trial is
highlighted after you select it. After you select a trial you can press the Delete button
on the right side of the window. A window will open where you must confirm that you
wish to delete the selected trial as shown in Figure 228. 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 Optimization Trial Manager window.
258
Figure 228.
Figure 229.
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 229). 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.
Figure 230.
The optimization trial component editor can be used to change the run
used in a trial, to select the search method, and to control the 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 (Figure 230). Click on the optimization trial node to
display the Component Editor for the optimization trial. If necessary, click on the
"Optimization Trials" folder to expand it and view the available optimization trials in
the project.
259
Objective Function
The objective function measures the goodness-of-fit between the computed outflow
and observed streamflow at the selected element. Seven different functions are
provided that measure the goodness-of-fit in different ways. 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 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 RMS log error uses the root-mean-square error of the logarithmtransformed flows to emphasize small flows compared to large flows. The timeweighted function gives greater weight to error near the end of the optimization time
window and less weight to errors early in the window.
Figure 231.
260
Select the objective function method for the optimization trial on the Component
Editor for the objective function (Figure 231). 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.
Click on the objective function node to view the editor. The objective function method
is the first item in the Component Editor.
An element with observed flow was selected at the time the optimization trial was
created. This is the location where the objective function will be evaluated. You can
change the element location at any time from the Component Editor for the objective
function. The selection list shows all of the elements with observed flow in the
simulation run selected as the basis for the optimization trial. Select a different
element in the list to change where the objective function will be evaluated.
Recognize that parameters can only be estimated at locations upstream of the
selected element. Changing the location will change which parameters can be
estimated.
In general, the observed flow record at the selected element location should not
contain any missing data. However, you can specify the amount of missing data in
the observed flow that will be accepted. Specify the amount that may be missing as
a percentage. The program will check the observed flow record and only perform the
optimization if the percentage of missing flow in the record is less than the specified
amount. Time steps with missing data are ignored when computing the value of the
objective function. The default amount of missing flow is 0.0 percent.
The objective function is evaluated over a specified time window. The time window
cannot begin before the start time in the control specifications used in the underlying
simulation run. Also, the time window cannot end after the end time in the control
specifications. However, you have the option of changing the time window to be
narrower than the one given by the control specifications. By default, the start and
end time for the objective function will default to the time window in the control
specifications of the underlying simulation run.
261
262
Figure 232.
Figure 233.
Table 24.
Method
Parameter
Deficit Constant
Initial Deficit
Constant Loss Rate
Maximum Storage
Recovery Factor
Initial Constant
Initial Loss
Constant Loss Rate
Exponential
Initial Range
Exponent
Loss Coefficient Ratio
Initial Loss Coefficient
Green Ampt
Initial Loss
Hydraulic Conductivity
Wetting Front Suction
Moisture Deficit
Curve Number
Initial Abstraction
Curve Number
Grid CN
SMA
Canopy Capacity
Canopy Initial Percent
Surface Capacity
Surface Initial Percent
Soil Capacity
Tension Zone Capacity
Soil Initial Percent
Groundwater 1 and 2 Capacity
Groundwater 1 and 2 Initial Percent
Maximum Infiltration Rate
Maximum Soil Percolation Rate
Maximum Groundwater 1 and 2 Percolation Rate
Groundwater 1 and 2 Storage Coefficient
Smith Parlange
Initial Content
Residual Content
Saturated Content
Bubbling Pressure
Pore Distribution
Hydraulic Conductivity
Beta Zero
263
Table 25.
Method
Parameter
Clark
Time of Concentration
Kinematic Wave
Channel Manning's n
Storage Coefficient
Collector Manning's n
Subcollector Manning's n
ModClark
Time of Concentration
Storage Coefficient
SCS
Time Lag
Snyder
Peaking Coefficient
Standard Lag
Table 26.
Method
Bounded Recession
Parameter
Initial Flow Rate
Initial Flow Rate per Area
Recession Constant
Linear Reservoir
Nonlinear Boussinesq
Recession
264
Table 27.
Method
Lag
Parameter
Time Lag
Modified Puls
Number of Steps
Muskingum
Muskingum K
Muskingum X
Number of Steps
Muskingum Cunge
Manning's n
Figure 234.
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 corresponding basin model in the
underlying simulation run. You may change the initial value without affecting the
basin model.
265
Checking Parameters
The program is equipped with the capability to check parameter values in the basin
and meteorologic model as part of performing an optimization trial. The check is
performed in two stages, similar to the way parameters are checked as part of a
simulation run. The first stage of checking is performed before the actual
optimization begins. During the first stage, required parameters are checked to make
sure they have been specified and are within the minimum and maximum values.
During the second stage, checking is performed that cannot happen until the actual
optimization begins.
Parameter checking is performed automatically as part of computing an optimization
trial. However, you can manually perform the first stage parameter checking. First
be certain that an optimization trial is selected. To perform the check, click on the
Compute menu and select the Check Parameters command. Any messages
generated during the check will be displayed in the message window.
Computing a Trial
There are three ways to compute the currently selected optimization trial. All
methods automatically perform parameter checking, and if no errors are generated,
proceed to the actual computing of optimization results. Additionally, the program is
designed to be computationally efficient. Only components with parameters selected
for optimization will be computed for each iteration, along with the components
downstream of those components; the remainder of the components have not
changed so previous optimization results are still applicable. If you wish, you may
force all components to be recomputed for each iteration regardless of data changes.
266
The first way to compute an optimization trial is from the Compute menu. Click on
the Compute menu and then select the Compute Trial command at the bottom.
The name of the current optimization trial is shown in brackets as part of the menu
command. If the command is not available, it is because there is no current
optimization trial; you must first select a trial. A window will automatically open that
shows the progress of the compute. You will need to manually close the window
when the compute is done, whether it failed or was successful. If you wish to force
all components to be recomputed instead of just those with data changes, hold the
control key while selecting the menu command.
The second way to compute an optimization trial is from the Watershed Explorer, on
the "Compute" tab. Select the optimization trial you wish to compute by clicking on it
in the Watershed Explorer with the right mouse button. A context menu is displayed
that contains several choices including compute. Click the Compute command. The
progress window will automatically open. If you wish to force all components to be
recomputed instead of just those with data changes, hold the control key while
selecting the right mouse menu command.
The third way to compute an optimization trial is from the toolbar. The compute
button is enabled when ever there is a current trial that can be computed. If the
button is not available or shows a different type of compute, you must first select an
optimization trial. Press the button to compute the current optimization trial. The
progress window will automatically open. If you wish to force all components to be
recomputed instead of just those with data changes, hold the control key while
selecting the toolbar button.
267
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 235.
Figure 235.
268
Figure 236.
Figure 237.
269
Figure 238.
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.
270
Figure 239.
Flow residuals showing the computed minus the observed flow at the
optimization location.
Figure 240.
271
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. The information included in
the graph (Figure 241) varies by element type, but always includes outflow. Optional
items such as observed flow, computed stage, and observed stage are also included.
The information included in the summary table (Figure 242) also varies by element
type but always includes the peak flow, time of peak flow, and outflow volume. The
time-series table (Figure 243) 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 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.
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. 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 simulation run 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 simulation run do not need to be recomputed,
and the selected element is upstream of the objective function evaluation element
with observed flow.
272
Figure 241.
Figure 242.
273
Figure 243.
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 simulation run 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.
274
To begin viewing results, go to the "Results" tab of the Watershed Explorer and click
on the desired optimization trial icon. If necessary, click on the "Optimization Trials"
folder to expand it and view the optimization trials in the project. The optimization
trial icon will be disabled if any data used in the trial has changed and results need to
be recomputed. If any result is open at the time data changes, the affected results
will be automatically updated.
Trial Results
The objective function summary can be accessed by clicking on the "Objective
Function Summary" node in the Watershed Explorer (Figure 244). 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.
Figure 244.
Individual Elements
Each element above the objective function evaluation location with observed flow is
shown in the Watershed Explorer under the optimization trial node. These elements
are listed in alphabetical 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 244). 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.
275
Figure 245.
276
Figure 246.
After you have opened a time-series table or graph, you may add additional timeseries 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 timeseries. 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.
277
CHAPTER 14
Analyses
Analyses are one of the three different components that can compute results:
simulation runs, optimization trials, and analyses. Each analysis is based on one or
more simulation runs. The run provides the basic framework of a basin model,
meteorologic model, and control specifications. Different types of analyses require
different types and amounts of additional information in order to provide value-added
processing to the simulation results. The results available from each analysis
depend on the type of analysis. Specific information on the required parameters and
available results for each type of analysis is given later in this chapter in the section
describing each type of analysis.
279
you have chosen to create. Your progress will be shown in the title of the wizard,
where it indicates the current step and the total number of steps. At any step you can
press the <Back button to return to a previous step and change you choices. When
you arrive at the last step, press the Finish button to create the new analysis. You
will subsequently need to use the Component Editor for the analysis to finish
specifying properties before it can be computed.
280
Figure 247.
Figure 248.
The first step of creating a new analysis is to select the type of analysis
to create. The remaining steps depend on the type of analysis but
always including giving it a name.
Copying an Analysis
There are two ways to copy an analysis. Both methods for copying an analysis
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 Analysis Manager, which is accessed
from the Compute menu. First select the type of analysis you wish to copy using the
selection list at the top of the manager. With the correct analysis type showing,
select the analysis you wish to copy by clicking on it in the list of current analyses.
The selected analysis is highlighted after you select it. After you select a analysis
you can press the Copy button on the right side of the window. A new Copy
Analysis window (Figure 249) 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 analysis. You cannot press the Copy button if no
name is specified. If you change your mind and do not want to copy the selected
analysis, press the Cancel button or the X button in the upper right to return to the
Analysis Manager window.
Figure 249.
The second way to copy is from the "Compute" tab of the Watershed Explorer. You
may need to first click on the "Analyses" folder to expand it, then click on the folder
for the correct type of analysis. Move the mouse over the analysis you wish to copy
and press the right mouse button (Figure 250). A context menu is displayed that
contains several choices including copy. Click the Create Copy command. A new
Copy 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
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 analysis. You cannot press the Copy button if no
name is specified. If you change your mind and do not want to copy the selected
analysis, press the Cancel button or the X button in the upper right of the Copy
Analysis window to return to the Watershed Explorer.
281
Figure 250.
Renaming an Analysis
There are two ways to rename an analysis. Both methods for renaming an analysis
change its name and perform other necessary operations.
The first way to perform a rename is to use the Analysis Manager, which you can
access from the Compute menu. First select the type of analysis you wish to
rename using the selection list at the top of the manager. With the correct analysis
type showing, select the analysis you wish to rename by clicking on it in the list of
current analyses. The selected analysis is highlighted after you select it. After you
select an analysis you can press the Rename button on the right side of the
window. A new Rename Analysis window (Figure 251) 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
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 analysis, press the
Cancel button or the X button in the upper right of the Rename Analysis window to
return to the Analysis Manager window.
Figure 251.
The second way to rename is from the "Compute" tab of the Watershed Explorer.
Select the analysis you wish to rename by clicking on it in the Watershed Explorer; it
will become highlighted. You may need to first click on the "Analyses" folder to
expand it, then click on the folder for the correct type of analysis. Keep the mouse
over the selected analysis and click the left mouse button again. The highlighted
name will change to editing mode (Figure 252). You can then move the cursor with
282
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
"Compute" tab. If you change your mind while in editing mode and do not want to
rename the selected analysis, press the Escape key.
Figure 252.
Deleting an Analysis
There are two ways to delete an analysis. Both methods for deleting an analysis
remove it from the project and automatically delete previously computed results.
Once an analysis has been deleted it cannot be retrieved or undeleted.
The first way to perform a deletion is to use the Analysis Manager, which you can
access from the Compute menu. First select the type of analysis you wish to delete
using the selection list at the top of the manager. With the correct analysis type
showing, select the analysis you wish to delete by clicking on it in the list of current
analyses. The selected analysis is highlighted after you select it. After you select a
analysis you can press the Delete button on the right side of the window (Figure
253). 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 analysis, press the Cancel button or the
X button in the upper right to return to the Analysis Manager window.
The second way to delete is from the "Compute" tab of the Watershed Explorer.
Select the analysis you wish to delete by clicking on it in the Watershed Explorer; it
will become highlighted. You may need to first click on the "Analyses" folder to
expand it, then click on the folder for the correct type of analysis. Keep the mouse
over the selected analysis and click the right mouse button. A context menu is
displayed that contains several choices including delete (Figure 254). 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 analysis, press the Cancel
button or the X button in the upper right to return to the Watershed Explorer.
283
Figure 253.
Figure 254.
284
Checking Parameters
The program is equipped with the capability to check parameter values in the basin
and meteorologic model as part of performing an analysis. The check is performed in
two stages, similar to the way parameters are checked as part of a simulation run.
The first stage of checking is performed before the actual analysis begins. During the
first stage, required parameters are checked to make sure they have been specified
and are within the minimum and maximum values. During the second stage,
checking is performed that cannot happen until the actual analysis begins.
Parameter checking is performed automatically as part of computing an analysis.
However, you can manually perform the first stage parameter checking. First be
certain that an analysis is selected. To perform the check, click on the Compute
menu and select the Check Parameters command. Any messages generated
during the check will be displayed in the message window.
Computing an Analysis
There are three ways to compute the currently selected analysis. All methods
automatically perform parameter checking, and if no errors are generated, proceed to
the actual computing of analysis results. Additionally, the program is designed to be
computationally efficient. Only components with data changes since the last time the
analysis was computed will be recomputed, along with the components downstream
of those components; the remainder of the components have not changed so
previous results are still applicable. If you wish, you may force all components to be
recomputed regardless of data changes.
The first way to compute an analysis is from the Compute menu. Click on the
Compute menu and then select the compute command at the bottom of the menu.
The exact wording of the compute command will depend on the type of analysis that
is selected. When a depth-area analysis is selected, the command will read
Compute Depth-Area. The name of the current analysis is shown in brackets as
part of the menu command. If the command is not available, it is because there is no
current analysis; you must first select an analysis. A window will automatically open
that shows the progress of the compute. You will need to manually close the window
when the compute is done, whether it failed or was successful. If you wish to force
all components to be recomputed instead of just those with data changes, hold the
control key while selecting the menu command.
The second way to compute an analysis is from the Watershed Explorer, on the
"Compute" tab. Select the analysis you wish to compute by clicking on it in the
Watershed Explorer with the right mouse button. You may need to first click on the
"Analyses" folder to expand it, then click on the folder for the correct type of analysis.
A context menu is displayed that contains several choices including compute. Click
the Compute command. The progress window will automatically open. If you wish
to force all components to be recomputed instead of just those with data changes,
hold the control key while selecting the right mouse menu command.
The third way to compute an analysis is from the toolbar. The compute button is
enabled when ever there is a current analysis that can be computed. If the button is
not available or shows a different type of compute, you must first select an analysis.
Press the button to compute the current analysis. The progress window will
automatically open. If you wish to force all components to be recomputed instead of
just those with data changes, hold the control key while selecting the toolbar button.
285
Depth-Area Analysis
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
2
2
of a 200 km storm is generally much less than the intensity of a 10 km 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
computes, it automatically generates frequency storms based on the one specified in
the meteorologic model but with appropriate storm area for each analysis point. The
instantaneous peak flows at each analysis point are automatically collected and
available in a tabular display.
286
node to display the Component Editor for the depth-area analysis. You may need to
first click on the "Analyses" folder to expand it, then click on "Depth-Area Analyses"
folder. The "Depth-Area Analysis" tab is shown along with the "Analysis Points" tab,
as shown in Figure 255. The selection list shows all of the simulation runs that have
an appropriate meteorologic model using the frequency precipitation method. You
can also press the simulation run button next to the selection list to bring up a
chooser. The chooser helps to make a selection by showing the description of each
run.
Figure 255.
Figure 256.
Viewing Results
The results available for a depth-area analysis are similar to the results for
simulations runs and optimization trials. When a depth-area analysis is selected,
287
certain results are available for it on the Results menu and from the right-mouse
menu for each element. These results can be quickly accessed while navigating the
analysis though the basin model map window. Results are also available on the
"Results" tab of the Watershed Explorer. These results provide additional flexibility
since they can be accessed even when a depth-area analysis is not selected.
Results accessed here can also be combined with results from simulation runs or
optimization trials for comparison.
Figure 257.
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.
Access to the peak flow summary table is provided in multiple locations. The
summary for the current selected depth-area analysis can be accessed by clicking on
the Results menu and selecting the Peak Flow Summary command. The menu
command is only available when a depth-area analysis is selected, any the results
are available. The results are only available if the analysis has been computed and
none of the data required for the analysis has changed since it was last computed.
The menu command is disabled if results need to be recomputed.
Access to the peak flow summary table for any depth-area analysis can be accessed
through the "Results" tab of the Watershed Explorer. Click on the analysis for which
you wish to view results. You may need to first click on the "Analyses" folder to
expand it, then click on the "Depth-Area Analyses" folder. Clicking on the analysis
will cause it to expand, as shown in Figure 258. The peak flow summary table is the
first item listed under the analysis name. Click on the summary table node to view
the results.
288
Figure 258.
Individual Elements
Typical results are available for the individual elements including graphs, summary
tables, and time-series tables. What is different for depth-area analyses is that
multiple sets of results may be available for an element. One set of results will be
produced at an element for each analysis point that is downstream of the element.
The multiple results correspond to the separate hypothetical storm generated for
each analysis point. Recall that the drainage area at an analysis point is used to
configure a hypothetical storm. The storm is applied to the elements upstream of the
analysis point in order to produce results. The results from each storm are stored
separately. For example, if a particular element is upstream from three different
analysis points, then the element will have three sets of results corresponding to the
three points. Since the results are stored separately, they can all be viewed. The
standard graph, summary table, and time-series table can be viewed that resulted
from the storm computed for each downstream analysis point.
Access to the results for an element is provided in multiple locations. The results can
be accessed through the Results menu. To begin you must select the analysis point
for which subsequent results will be viewed. On the Results menu, move to the
Select Point menu and choose an analysis point. You will be able to view a graph,
summary table, or time-series table for the element once a point is selected. The
results can be accessed by clicking on the appropriate menu command on the
Results menu, or my right-clicking on the element in the basin map. In order to view
results for the same element above a different analysis point you will need to select a
different point on the Results menu. Once a point is selected it will be used for all
graphs, summary tables, and time-series tables selected through the Results menu
or basin map until a new point is selected.
289
All of the results for the analysis points with their associated elements are available
on the "Results" tab of the Watershed Explorer. Click on the analysis for which you
wish to view results. You may need to first click on the "Analyses" folder to expand it,
then click on the "Depth-Area Analyses" folder. Clicking on the analysis will cause it
to expand. Below the peak flow summary table is a folder is for each analysis point,
as shown in Figure 258. Within each folder are icons for each of the elements
upstream of the analysis point. At each element you can access the standard graph,
summary table, and time-series table. It is also possible to access the individual
time-series computed in the element.
Combining Results
All of the time-series data computd by an individual element are available for viewing.
The time-series data are listed under each element node on the "Results" tab of the
Watershed Explorer. Results can be combined by clicking on a time-series node,
and then holding the control key while clicking on additional nodes (Figure 259). The
selected time-series may come from different elements above the same analysis
point, the same element for different analysis points, or different elements in different
analyses. The selected time-series may even come from simulation runs,
optimization trials, and analyses. The selected time-series data will automatically be
partitioned into groups based on data type.
Preview graphs of selected time-series data can be opened as graphs or time-series
tables within the Desktop area. Begin by selecting the time-series you wish to
include in the graph or table. Once you have selected the desired time-series, you
can press the graph or time-series table buttons on the toolbar. The chosen timeseries will be graphed or tabulated.
After you have opened a time-series table or graph, you may add additional timeseries 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 fchange 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. Depth-area
analysis results can be added to tables or graphs that hold only analysis results, or
can be mixed with results from simulations runs and optimization trials.
All of the graphs that can be accessed open with default properties for line color, line
style, data symbols, etc. These default properties have been selected to be
appropriate for most situations. However, it is possible to customize the properties in
a graph. To change the properties, first click on the graph to select it. Next go to the
Results menu and select the Graph Properties command. An editor will open that
can be used to change the properties of the selected graph. The properties for each
time-series curve can be changed. It is also possible to change the properties for the
axis, title, gridlines, patterns, and legend. The editor used to change properties is the
same as used for simulation run graphs shown in Figure 221.
290
Figure 259.
Reviewing the outflow for the same element above two different
analysis points in the same depth-area analysis.
291
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 timeseries 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.
ELEVATION-FLOW An elevation-discharge curve defined in the reservoir element.
No longer used in the current program version.
ELEVATION-STORAGE A manually entered elevation-storage function defined in
the paired data manager.
ET-CANOPY Computed actual evapotranspiration from the canopy layer in the soil
moisture accounting loss method.
ET-POTENTIAL Potential evapotranspiration for a subbasin element as computed
by the meteorologic model.
293
ET-SOIL Computed actual evapotranspiration from the soil layer in the soil moisture
accounting loss method.
ET-SURFACE Computed actual evapotranspiration from the surface depression
layer in the soil moisture accounting loss method.
EXCESS-CANOPY Computed incremental precipitation minus canopy interception
in the soil moisture accounting loss method.
FLOW Final computed flow for an element.
FLOW-AUX Outflow from a reservoir through the auxiliary discharge.
FLOW-BASE Computed baseflow for a subbasin element.
FLOW-COMBINE Computed total inflow to an element.
FLOW-DIRECT Surface flow computed by transforming excess precipitation for a
subbasin element.
FLOW-DIVERSION Computed diversion flow for a diversion element.
FLOW-IN Total inflow to a reach element.
FLOW-LOCAL Local flow at a junction when the basin model is set to compute local
flow at junctions.
FLOW-LOCAL-SIM Computed local flow without blending for a junction when
observed flow is present, blending is used, and the basin model is set to compute
local flow at junctions.
FLOW-LOSS At a reach, the computed channel losses.
FLOW-OBSERVED Observed flow for an element using the optional observed flow
time-series gage.
FLOW-RESIDUAL Residual flow for an element with observed flow. The residual is
calculated as computed flow minus the observed flow.
FLOW-SIM Computed flow without blending at an element when observed flow is
present and blending is used.
FLOW-UNIT GRAPH Manually entered unit hydrograph in the paired data manager.
GROUNDMELT Manually entered annual groundmelt pattern in the paired data
manager.
INFILTRATION Computed infiltration from the surface layer to the soil in the soil
moisture accounting loss method. Also, the amount of loss attributed to infiltration in
the deficit constant loss method.
LIQUID WATER Computed liquid water in the snowpack in snowmelt.
MELTRATE Manually entered annual meltrate pattern in the paired data manager.
MOISTURE DEFICIT Computed moisture deficit in the deficit constant and gridded
deficit constant loss methods.
294
295
STAGE-SIM
Stage at a reach using the modified Puls or Muskingum Cunge
routing methods computed using the simulated flow depth.
STAGE-TAILWATER The computed tailwater stage at a reservoir element.
STAGE-TW-AUX The computed tailwater stage at the auxiliary discharge from a
reservoir element.
STORAGE Computed storage for a reservoir element.
STORAGE-CANOPY Computed storage depth of the canopy layer in the soil
moisture accounting loss method.
STORAGE-ELEVATION A storage-elevation curve defined in the reservoir element.
No longer used in the current program version.
STORAGE-FLOW A manually entered storage-discharge function defined in the
paired data manager.
STORAGE-GW Computed storage depth of a groundwater layer in the soil moisture
accounting loss method.
STORAGE-SOIL Computed storage depth of the soil layer in the soil moisture
accounting loss method.
STORAGE-SURFACE Computed storage depth of the surface layer in the soil
moisture accounting loss method.
SWE Computed snow water equivalent in snowmelt.
SWE-OBSERVED Observed snow water equivalent for a subbasin element using
the optional observed SWE time-series gage.
SWE-RESIDUAL Residual snow water equivalent for a subbasin element with
observed SWE. The residual is calculated as computed flow minus the observed
SWE.
TEMPERATURE Computed temperature for a subbasin elevation band in
snowmelt. Also, manually entered time-series of temperature for a time-series gage.
296
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
260). The keyword "End" must be on a line by itself after the "Parameter Order" line.
Valid parameter keywords are shown in Table 28. 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 28.
Keyword
Definition
Units
XCoord
integer value
YCoord
integer value
TravelLength
Area
kilometers
square kilometers
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.
297
Figure 260.
298
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 261).
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.
299
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 261.
300
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:
Recession Constant =
( RTIOR )24
A = y(w+ z y )
301
(10 w z y ) + 16 z y 2 z 2 + 1 + (5 w 2 ) + 6 w y z 2 + 1
c
=
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 =
and the wave speed, c, is computed as:
c = m A m 1
where
302
303
APPENDIX E
305
306