Introduction
The Exponential Decay model applies an exponential decay to the stored constituent mass in a node or link.
Scale
This model operates at a river reach scale.
Principal Developer
eWater Solutions
Authors
Rachel Blakers, Andrew Davidson
Version
Source v4.1
Availability/Conditions
The Exponential Decay model is automatically installed with Source. For Lumped constituent routing, it can be applied to constituents in:
...
The decay model is not currently available for Storage Routing Links if Marker routing is used.
Definition of Symbols
Symbol | Description |
M(t) | Constituent mass in storage (instantaneous) |
V(t) | Storage volume (instantaneous) |
Qin(t) | Water volume flowing into the storage (instantaneous) |
Cin(t) | Constituent concentration flowing into the storage (instantaneous) |
Qout(t) | Water volume flowing out of the storage (instantaneous) |
Cout(t) | Constituent concentration flowing out of the storage (instantaneous) |
k | Exponential decay constant |
t | Time |
h | Exponential decay half-life (input by user) |
Background
Exponential Decay
Exponential decay describes the decay of a quantity at a rate that is directly proportional to the amount present (Equation 1).
...
|
|
Constituent Decay in a Fully Mixed Store
In Source, the model describing the exponential decay of a constituent must consider inflows and outflows from the storage (here the term "storage" is used in a general sense and refers to the water stored in a river reach or reservoir, for example). Assuming instantaneous, fully mixed conditions, the change in constituent mass over time can be determined from the following conservation of mass relationship:
...
Equation 9 forms the basis of the Exponential Decay model in Source. The solution becomes less accurate as the half-life becomes smaller (particularly if the half-life is smaller than the time step) and as the flushing rate (Qout /Vt+1) becomes larger. It also follows from Equation 9 that if Vt+1 = 0 then Mt+1 = 0.
Example
If the half-life is equal to the time step, then the constituent mass after one time step will be equal to one half of the previous mass (ignoring losses and gains from outflows and inflows). For example, consider a half-life of h = 1 day, a daily time step, zero inflows and zero outflows. Equation 9 tells us that, after 1 day, the constituent mass will be:
...
which is one half of the initial mass.
Solution Methodology - Flow Phase
The Exponential Decay model can be applied in Storage and Reservoir nodes for both Marker and Lumped constituent routing. Currently, it can only be applied in Storage Routing Links if Lumped constituent routing is used. The model is run each time step after the flow phase calculations for the relevant node or link:
...
CASE 1: If the end-of-time-step storage volume is 0 and the outflow is 0, then all constituent mass is deposited.
CASE 2: If the end-of-time-step storage volume is 0 and the outflow is greater than 0, then all constituent mass leaves the storage with the outflow.
CASE 3: If the end-of-time-step storage volume is greater than 0, then Equation 9 applies.
|
Input Data
A single value for the half-life parameter in seconds. Smith et al. (2011) provide some estimates of pesticide half-lives, Birgand et al. (2007) and Rao et al. (2009) provide estimates of decay rates for nutrients in streams.
Parameters or Settings
Model parameters are summarised in Table 1.
| Anchor | ||||
|---|---|---|---|---|
|
Parameter | Description | Units | Default | Range |
Half-life | The half-life of the constituent | seconds | 86400 (equivalent to 1 day) | 0 to ∞ |
Outputs
A time series of stored constituent mass.
Reference List
Birgand, F., R.W. Skaggs, G.M. Chescheir, J.W. Gilliam (2007) Nitrogen Removal in Streams of Agricultural Catchments-A Literature Review. Critical Reviews in Environmental Science and Technology; 2007; 37, 5; ProQuest Agriculture Journals, p. 381
Rao, P.S.C., N.B. Basu, S. Zanardo, G. Botter, A. Rinaldo (2009) Contaminant load-discharge relationships across scales in engineered catchments: Order out of complexity. 18th World IMACS / MODSIM Congress, Cairns, Australia 13-17 July 2009, p. 1886-1892. http://mssanz.org.au/modsim09
Smith, R., R. Turner, S. Vardy, M. Warne (2011) Using a convolution integral model for assessing pesticide dissipation time at the end of catchments in the Great Barrier Reef Australia. In F. Chan, D. Marinova, R.S. Anderssen (eds) Modsim2011, 19th International Congress on Modelling and Simulation. Modelling and Simulation Society of Australia and New Zealand, December 2011, pp. 2064-2070. ISBN: 978-0-9872143-1-7. http://www.mssanz.org.au/modsim2011/E5/smith.pdf