Versions Compared

Key

  • This line was added.
  • This line was removed.
  • Formatting was changed.

The Sacramento Model is a catchment water balance model that relates runoff to rainfall with daily data. The model contains five stores and has sixteen parameters.

Scale

Model operates at a functional unit scale and a daily time-step.

Principal developer

Cooperative Research Centre for Catchment Hydrology. The Sacramento Model was originally developed for the United States National Weather Service and State of California, Department of Water Resources by Burnash et al (1973).

Scientific provenance

The Sacramento Model has been widely applied in the USA by the National Weather Service. In Australia, the Sacramento Model has been the model of choice in water resources studies that have applied IQQM and has therefore been very widely applied in Queensland and New South Wales. There have been applications in other Australian states, although these have been fewer in number.

Version

Source v3.8.12

Availability/conditions

The Sacramento Model is provided with the Source installer. It is also available through the Rainfall Runoff Library on eWater Toolkit.

Theory

Introduction

The Sacramento Model is a continuous rainfall-runoff model used to generate daily stream flow from daily rainfall and potential evapotranspiration data. It uses soil moisture accounting to simulate the water balance within the catchment (i.e. functional unit). The conceptual layout of the model is shown in Figure 1.  At each model time step the sequence of calculations is:

 

  1. Soil moisture depletion by evapotranspiration and soil moisture redistribution.
  2. Soil moisture replenishment by rainfall and percolation, and streamflow generation.

 

The internal Sacramento Model calculations represent water quantities using units of depth (in millimetres). The model outputs are converted to volumes by multiplying by the catchment area.

 



Figure 1. Sacramento Model

Selecting stream flow data to use in a river-basin-scale simulation study needs information about the reliability of the data. It is best to use data which are most representative of the stream flow from the catchment. Observed data would normally be selected, except where the data are of poor quality or of unknown reliability.

 

There are five stores in the Sacramento Model:

  • Upper zone tension water (UZTW);
  • Upper zone free water (UZFW);
  • Lower zone tension water (LZTW);
  • Lower zone primary free water (LZFWP); and
  • Lower zone supplementary free water (LZFWS).

The tension water stores represent the volume of water that is held in the soil matrix by surface tension. Water can only be removed from the tension water stores by evapotranspiration. In the case of the free water stores, water can move through the soil vertically and laterally to other stores, and be discharged as interflow (upper zone) or baseflow (lower zone).

The Sacramento Model divides the catchment into impervious and pervious areas. The impervious area is the portion of the catchment that is covered by lakes, rivers, pavement and other impervious surfaces that are directly connected to the stream network. The impervious area produces runoff from any rainfall while the pervious area only produces runoff when rainfall is sufficiently heavy. The Sacramento model also allows the user to specify a variable impervious area, which is a portion of the catchment that can become impervious when the catchment is quite saturated.

Streamflow generated by the Sacramento Model is made up of four main forms:

  1. impervious runoff from permanent impervious areas and direct runoff from variable impervious areas,
  2. surface runoff, which occurs when Upper Zone Free Water storage is full and the precipitation intensity exceeds the rate of percolation and interflow
  3. interflow resulting from the lateral drainage of the Upper Zone Free Water storage
  4. baseflow, which is composed of primary and supplemental baseflow

Impervious runoff, direct runoff and surface runoff have no time delay in the Sacramento model, occurring in the same time interval as the rain that generated the runoff component. Interflow generally has a time delay in terms of days, supplemental baseflow a delay of weeks or months and primary baseflow a delay of months or years.The generation of surface runoff, interflow and baseflow depends on the amount of water in each soil moisture store relative to that store’s capacity, and the rate at which water moves into and out of the stores.

Further details are in the following sections.

Soil moisture depletion and redistribution

Depletion

The upper limit of evapotranspiration is the potential evapotranspiration demand, ED, specified by the modeller. Soil moisture depletion by evapotranspiration takes place as follows.

Evapotranspiration occurs firstly from UZTW:

 

Equation 1 

Where:

UZTWC = current contents of upper zone tension water store (state variable)

UZTWM = capacity of upper zone tension water store (parameter)

Then, if there is potential demand unsatisfied, evapotranspiration occurs from UZFW

Equation 2

Where:

UZFWC = current contents of upper zone free water store (state variable)

 

If there is still potential demand unsatisfied then the remainder is satisfied as much as possible from LZTW:

Equation 3

Where:

LZTWC = current contents of lower zone tension water store (state variable)

LZTWM = capacity of lower zone tension water store (parameter)

 

There is also an additional store, with current contents of ADIMC, related to the parameter ADIMP, the additional area that behaves as impervious when UZTW is full. This store can be depleted by evapotranspiration as well:

Equation 4

 

Evaporation can also occur directly from the stream as governed by the modeller specified parameter, SARVA:

Equation 5Image Removed 

Redistribution

Soil moisture is redistributed from the free water stores to the tension water stores as follows:

  • If the ratio of UZFWC to UZFWM is greater than the ratio of UZTWC to UZTWM, water is transferred from UZFW to UZTW until the ratios are equal (but the reverse does not occur).
  • If the ratio of the combined current contents of the three lower zone stores (less the modeller specified RSERV component) to the combined capacities of these stores (also less the modeller specified RSERV component) is greater than the ratio of LZTWC to LZTWM, alone, then the ratios are equalised by firstly transferring water from LZFWS to LZTW, and secondly, if there is insufficient water in LZFWS to achieve equalisation, by transferring water from LZFWP to LZTW.

 

Soil moisture replenishment and percolation, and streamflow generation

Overview

The sequence of calculations in this phase, for each model time step, is as follows:

  1. Replenish the upper zone tension water store.
  2. Compute supplemental and primary baseflow from the two lower zone free water stores.
  3. Compute total percolation from the UZFW store to the three lower zone stores.
  4. Compute interflow (from the UZFW store).
  5. Distribute percolation water between the lower zone free water stores, collectively, and the LZTW store.
  6. Distribute percolation water available to lower zone free water stores between primary and supplemental stores.
  7. Replenish upper zone free water store (UZFW) with any overflow from upper zone tension water store and, if UZFW overflows, compute surface runoff.
  8. Calculate instantaneous runoff, being the sum of surface runoff, impervious runoff ROIMP, and any contribution due to ADIMP.

Impervious areas

All rain falling on the permanently impervious area (defined by PCTIM) runs off; that is, it is not subject to depletion by evaporation. Hence:

 

 

Water movement through the stores is determined by rules, where the UZTW store receives the rain first, and when this is filled water will go to the UZFW store. The UZFW store then supplies water to the lower stores simultaneously, with a user determined split between the free water and tension water stores. When the LZFWS is filled water will go to the tension water stores.

Equation 2
Equation 3
Equation 4

Stream flow generated with the Sacramento model is made up of three flow components:

  • Surface runoff;
  • Interflow; and
  • Base flow

The generation of these components depends on the amount of water in each store relative to that store’s capacity, and the rate at which water moves into and out of the stores.

Surface runoff is either direct or occurs when UZTWS is full and the rainfall exceeds the sum of the percolation rate and the maximum interflow drainage capacity.

Interflow is generated from the UZFWS as the product of the volume of water in the store, and a drainage rate parameter, UZK.

Equation 5

Base flow is calculated in a similar manner to interflow, using the volume of water in the lower zone free stores with their corresponding drainage rate parameters, LZPK and LZSK. The base flow is then reduced by channel loss parameters, SIDE and SSOUT.

Equation 6
Equation 7
Equation 8
Equation 9

 

Evapotranspiration can only take place from upper and lower tension water stores and upper free water stores, and directly from the streams. The upper limit of evaporation is the evaporative demand, and is the product of the pan evaporation modified by the (user- specified) pan factor. Evaporation occurs firstly from the UZTWS, then from the UZFWS, and lastly from the LZTWS. Evaporation can also occur directly from the stream as set by SARVA.

Equation 10

Equation 11
Equation 12
Equation 13
Equation 13A

 

The percolation to the lower stores is a key process of the Sacramento model. The driving force for percolation is the relative wetness of the UZFWS as moderated by the relative wetness of the lower zone stores.

Percolation increases when either the storage in the UZFW store increases or the storage in the lower zone stores decrease. This is equivalent to supply increasing and demand increasing respectively. Conversely, percolation decreases when the lower stores start becoming full.

The lower limit of percolation, Pbase, occurs when the lower zones are saturated, and is determined by the rate at which the lower zones drain (Equation 14). The maximum rate of percolation occurs when the lower zones are dry, and Pbase is factored up using the ZPERC parameter (Equation 15).

Equation 14
Equation 15

The actual percolation is moderated by the relative saturation of the lower and upper zones, which is the ratio of actual storage to maximum storage in these stores, to give an estimate of percolation (Equation 16).

Equation 16

Normally, the lower zone tension store would fill before water goes to the lower zone free water store. However, variations in soil types cause deviations from average conditions and therefore in the Sacramento model a fraction of the percolation (PFREE) is made available for lower zone free water stores.

Equation 17
Equation 18
Equation 19
Equation 20
Equation 21

 

Time delay tools

The direct runoff and interflow can be delayed to better represent stream flow hydrographs. The Sacramento model uses a unit hydrograph for this purpose. Each ordinate (UH1...5) represents the proportion of flow which will reach the channel outlet in successive time periods. Sacramento also uses a layered routing process.

Other factors

Other inputs that can also be altered in the calibration process include the parameter RSERV (Fraction of lower zone free water unavailable for transpiration).

Input data

The model requires daily rainfall and potential evapotranspiration data. The rainfall and evaporation data sets need to be continuous (no gaps) and overlapping.

Daily rainfall data may be obtained from rain gauges or rainfall surfaces but will need to be converted to a time series record that is spatially representative of the whole catchment. Note that the time that rainfall data are collected may be important. Very often rainfall data are collected in the morning, the usual time is 9 am, and may be more representative of the previous day’s rainfall.

Daily evaporation is an estimate of the spatially averaged evaporation rate of the catchment being modelled. This estimate is subject to the types of land uses that are in the catchment. This may be estimated by applying a crop/land use factor to daily pan or potential evapotranspiration surface data.

Daily flow data in ML/day, m3/s or mm/day may be required to calibrate the model.

Parameters or settings

The Sacramento model uses a total of sixteen parameters to simulate the water balance. Of these:

  • five define the size of soil moisture stores;
  • three calculate the rate of lateral outflows;
  • three calculate the percolation water from the upper to the lower soil moisture stores;
  • two calculate direct runoff; and
  • three calculate losses in the system
Table 1. Model Parameters
ParameterDescriptionUnitsDefaultMinMax

UZTWM

Upper Zone Tension Water Maximum. The maximum volume of water held by the upper zone between field capacity and the wilting point which can be lost by direct evaporation and evapotranspiration from soil surface. This storage is filled before any water in the upper zone is transferred to other storages.

mm

50

0

100

UZFWM

Upper Zone Free Water Maximum, this storage is the source of water for interflow and the driving force for transferring water to deeper depths.

mm

40

0

80

LZTWM

Lower Zone Tension Water Maximum, the maximum capacity of lower zone tension water. Water from this store can only be removed through evapotranspiration.

mm

130

0

400

LZFSM

Lower Zone Free Water Supplemental Maximum, the maximum volume from which supplemental base flow can be drawn.

m

23

0

50

LZFPM

Lower Zone Free Water Primary Maximum, the maximum capacity from which primary base flow can be drawn.

mm

40

0

50

UZK

The ratio of water in UZFWM, which drains as interflow each day.

 

0.245

0.000

1.000

LZSK

The ratio of water in LZFSM which drains as base flow each day.

 

0.043

0

1

LZPK

The ratio of water in LZFPM, which drains as base flow each day.

 

0.009

0

1

PFREE

The minimum proportion of percolation from the upper zone to the lower zone directly available for recharging the lower zone free water stores.

 

0.063

0.000

1.000

REXP

An exponent determining the rate of change of the percolation rate with changing lower zone water storage.

 

1

0

3

ZPERC

The factor applied to Pbase to define maximum percolation rate.

 

40

0

80

SIDE

The decimal fraction of observed base flow, which leaves the basin, as groundwater flow.

 

0

0

1

SSOUT

The volume of the flow which can be conveyed by porous material in the bed of stream.

 

0.001

0.000

1.000

PCTIM

The impervious fraction of the basin, and contributes to direct runoff.

 

0.01

0.00

1.00

ADIMP

The additional fraction of pervious area, which develops impervious characteristics under soil saturation, conditions. ADIMP spill only occurs when both tension stores are full.

 

0.01

0.00

1.00

SARVA

A decimal fraction representing that portion of the basin normally covered by streams, lakes and vegetation that can deplete stream flow by evapotranspiration.

 

0.01

0.00

1.00

RSERV

Fraction of lower zone free water unavailable for transpiration

 

0.3

0.0

1.0

UH1

The first component of the unit hydrograph, ie. the proportion of runoff not lagged

 

0.9

0

1

UH2

The second component of the unit hydrograph, ie. the proportion of runoff lagged by one time-step

 

0.1

0

1

UH3

The third component of the unit hydrograph, ie. the proportion of runoff lagged by two time-steps

 

0

0

1

UH4

The fourth component of the unit hydrograph, ie. the proportion of runoff lagged by three time-steps

 

0

0

1

UH5

The fifth component of the unit hydrograph, ie. the proportion of runoff lagged by four time-steps

 

0

0

1

In Source, up to five unit hydrograph terms (UH1...UH5) can be set to lag the runoff over time. When using RRL v1.0.5 the unit hydrograph term is fixed at 1 for the first time increment, and 0 for the subsequent time increments, and therefore no unit hydrograph routing is applied in RRL. Equivalent behaviour would be achieved with Sacramento in Source by setting UH1 = 1 and UH2 through UH5 to 0.

The sum of the unit hydrograph terms should always be 1.

The Sacramento model in Source is configured with a set of default values for each model parameter. There are also upper and lower bounds for each parameter.

As with any modelling, the accuracy and reliability of the results from the Sacramento model are determined by how representative the model is of the catchment (particularly as the Sacramento model is lumped) and also by the quality of the rainfall, evaporation and stream flow data used. The accuracy and reliability of the model can be assessed using the results of comparisons with observed data.

As a rule, the calibrated parameter values of a specific catchment should not be transposed to other catchments, unless the reliability of this transposition can be assessed. The parameter set is unique to the climate, topography, size, geology, soil and vegetation type of the catchment on which it was calibrated. There is no proven methodology to adjust these parameters to other catchments, including subcatchments, of the calibrated catchment.

Output data

The model outputs daily surface and base flow. This may be saved in ML/day, m3/s or mm/day.

Reference list

Burnash, RJC, Ferral, RL & McGuire, RA 1973, A generalized streamflow simulation system: conceptual modeling for digital computers, Technical Report, Joint Federal and State River Forecast Center, US National Weather Service and California Department of Water Resources, Sacramento, CA.

Bibliography

Rainfall Runoff Library v1.0.5, June 25, 2004 (http://www.toolkit.net.au/Tools/RRL).

NOAA Document explaining the conceptualisation of the Sacramento Model.

NOAA version of the soil moisture accounting model (in Fortran).

Info
iconfalse
Note: The RRL v1.0.5 Sacramento model does not have the multiple unit hydrograph terms that are included in the Sacramento model implemented in Source.