...
Source includes a number of optimisation techniques and statistical measures for automated model calibration and to assist modellers with the evaluation of the quality of calibration. These are mainly intended for application when calibrating catchment rainfall-runoff models in Source, but are also applicable when calibrating river system models (e.g. see Lerat et al., 2013). Optimisation techniques available The available automatic optimisation algorithms are:
- Shuffled complex evolution
- Genetic algorithms
- Uniform random sampling
- Rosenbrock method
Modellers have the option of selecting one optimisation technique , multiple optimisation techniques (in parallel), or combinations two optimisation techniques (in series).
Automated calibration requires the use of an objective function to direct the optimisation process. The Source calibration tool implements single objective function optimisation, which reduces the comparison between the observed and modelled data during the calibration period to a single number to be optimised (Multiple multiple objective optimisation is also available: , see Multi-objective optimisation /trade - off analysis - Insight - SRG for information). The following nine forms
Source implements five different basic types of objective function are available in Source:
- Match to Nash Sutcliffe Coefficient of Efficiency (NSE) of Daily Flows
- Minimise Absolute Bias between Observed and Modelled Flows (calculated using daily flows)
- Match to NSE of Daily Flows but Penalise Biased Solutions
- Match to NSE of Monthly Flows
- Match to NSE of Monthly Flows but Penalise Biased Solutions
- Combined Match to NSE and Match to Flow Duration Curve (Daily)
- Combined Match to NSE and Match to Logarithm of Flow Duration Curve (Daily)
- Combined Match to NSE of Logarithms of Daily Flows with Bias Penalty
- Combined Bias, Daily Flows and Daily Exceedance (Flow Duration) Curve (SDEB)
Further information on the first seven of these objective functions is available in Vaze et al (2011), Section 6. Guidance on model calibration is available in many publications, including various eWater Best Modelling Practice Guidelines (Black et al, 2011; Vaze et al, 2011; Black and Podger, 2012; and Lerat 2012).
Scale
...
- Flow duration (specifically, the NSE of the flow duration)
- Absolute bias
- Bias penalty
- Square-root daily, exceedance and bias
The NSE can be applied to daily or monthly data, and the NSE and flow duration objectives can be applied to data that has been transformed by taking the logarithm. Source also allows the user to create some composite objective functions, of which there are two types:
- Combinations of the individual objective functions listed above. For example, the objective for calibrating streamflow at a gauging site could be a combination of the NSE and bias penalty.
- Combinations of the objectives for different model outputs. For example, a model could be calibrated using a weighted combination of the objective functions values at two or more different gauging sites.
Scale
Typically, the optimisation techniques and statistical measures are used to compare observed and estimated data at a point, such as streamflow data at a gauging station. Both the optimisation techniques and statistical measures can be applied on a daily or monthly basis.
...
Overview information on the four optimisation techniques in Source is available in Vaze et al. (2011). Further information is in textbooks and papers, particularly for the genetic algorithm and uniform random sampling[DB1] . Publications on the shuffled complex evolution method include papers by Duan et al. (1992) and Sorooshian et al. (1993). Publications on the Rosenbrock method include the paper by Rosenbrock (1960).
...
Availability
Provided with Source.
Implementation
Structure & processes
...
Background
The optimisation techniques and statistical measures of calibration performance used in Source are well established , they and are not re- described in detail here. However, as the objective functions used in the optimisation techniques have been customised for Source, further information on these follows and as many of them rely on the Nash Sutcliffe Coefficient of Efficiency (NSE), its formulation is restated below.
The traditional formula for NSE is:
| Equation 1 |
|---|
where:
Qobsi is the observed flow on day i,
Qsimi is the modelled flow on day i,
N is the number of days
Alternatively,
| Equation 2 |
|---|
This formulation obviates the necessity to calculate the average of the observed flows before evaluating the denominator in the traditional version.
The choice of any particular objective function will depend on the intended application. Each of the pre-defined objective functions are formulated to put emphasis (reproduce as closely as possible) on different flow characteristics (Vaze et al, 2011).
- Match to Nash Sutcliffe Coefficient of Efficiency (NSE) of Daily Flows
Application of this objective function involves maximising the NSE (i.e. getting it as close to 1.0 as possible). The calculation of the NSE is in accordance with Nash and Sutcliffe (1970) and uses observed and modelled daily flow data for all days within the calibration period for which observed daily flow data, including zero flow values (i.e. cease to flow), is available.
The NSE tends to produce solutions that match high and moderate flows very well but often will produce poor fits to low flows. It will also tend to favour solutions that provide a good match to the timing and shape of runoff events (Vaze et al, 2011).
2. Minimise Absolute Bias between Observed and Modelled Flows
This objective function will produce a match on the overall volume of flow generated but often will produce a poor fit to the timing of flows (Vaze et al, 2011). It has the following form:
| Equation 3 |
|---|
The evaluation of this objective function uses observed and modelled daily flow data for all days within the calibration period for which observed daily flow data, including zero flow values, is available.
3. Match to NSE of Daily Flows but Penalise Biased Solutions
This objective function is a weighted combination of the daily NSE and a logarithmic function of bias based on Viney et al (2009), and the aim is to find its maximum value.
| Equation 4 |
|---|
where:
B is the bias; and
| Equation 5 |
|---|
The evaluation of this objective function uses observed and modelled daily flow data for all days within the calibration period for which observed daily flow data, including zero flow values, is available.
This formulation makes sure that the models are calibrated predominantly to optimise NSE while ensuring a low bias in the total streamflow. It avoids solutions that produce biased estimates of overall runoff, which can produce marginal improvements in low flow performance over the NSE objective function. However, NSE-Bias will still be strongly influenced by moderate and high flows and by the timing of runoff events, which can still often result in poor fits to low flows (Vaze et al, 2011).
4. Match to NSE of Monthly Flows
This objective function works in the same way as for the case “Match to NSE of Daily Flows” except that monthly flows are used to evaluate the NSE instead of daily flows. The Guidance on model calibration is available in many publications, including various eWater Best Modelling Practice Guidelines (Black et al., 2011; Vaze et al., 2011; Black and Podger, 2012; and Lerat, 2012).
The choice of an appropriate objective function for calibration depends on the intended application of the model. Different objective functions are designed with the intention of emphasizing the fit of modelled flow to different aspects of the observed hydrograph (Vaze et al., 2011). The objective functions available in the Source Calibration Wizard are listed in Table 1, including useful references for further information.
Table 1. List of Source calibration objective functions.Anchor Table 1 Table 1
| Objective Function Name | Description | Reference |
|---|---|---|
| NSE Daily | Maximise the NSE of daily flows | Vaze et al. (2011), Section 6 |
| NSE Monthly | Maximise the NSE of monthly flows | Vaze et al. (2011), Section 6 |
| NSE Log Daily | Maximise the NSE of the logarithm of daily flows | |
| Minimise Absolute Bias | Minimise the Absolute value of the relative bias | Vaze et al. (2011), Section 6 |
| NSE Daily & Bias Penalty | Maximise the NSE of daily flows and bias penalty | Vaze et al. (2011), Section 6 |
| NSE Log Daily & Bias Penalty | Maximise the NSE of the logarithm of daily flows and bias penalty | |
| NSE Monthly & Bias Penalty | Maximise the NSE of monthly flows and bias penalty | Vaze et al. (2011), Section 6 |
| NSE Daily & Flow Duration | Maximise the NSE of daily flows and the NSE of the flow duration | Vaze et al. (2011), Section 6 |
| NSE Daily & Log Flow Duration | Maximise the NSE of daily flows and the NSE of the flow duration of log flows | Vaze et al. (2011), Section 6 |
| Square-root Daily, Exceedance and Bias | Minimise a combination of the bias, daily Flows and daily exceedance (flow duration) curve | Lerat et al., 2013 |
Implementation Details
The Bivariate Statistics SRG User Guide entry provides general information on the objective function equations and their interpretation. Implementation details that are specific to the Source Calibration Wizard are described below.
Missing Data
It is common for observed time series of hydrological processes to contain missing values. Also, the observed and modelled time series may have different start and end dates. The Source calibration tool calculates the objective function values using only data from those time steps for which both observed and modelled data is available. In other words, the calibration objective function are calculated using observed and modelled data has been filtered to include only:
- data from within the calibration period, and
- data for time steps with complete data pairs.
Monthly Flows
The NSE Monthly objective function uses monthly streamflow values. These are calculated as follows:
- If the model is run on a daily time step, monthly flows are calculated by summing the observed and modelled daily flow values. The NSE calculation ignores observed and modelled data for all months where there are one or more days of missing data in the observed flow series.
- If the model is run on a monthly time step, then the monthly values are unchanged.
Data
Input data
Details on data to be input by the modeller are provided in the Source User Guide. Requirements for data series inputs to the various objective functions are included in the descriptions of each objective function, above.
Parameters or settings
Modellers have the option of selecting one optimisation technique, two optimisation techniques (in series), or manual optimisation. Modellers can also select which objective function they wish to use. The other parameters the modeller can input are described in Table 2.
Table 2. Objective function parameters. Refer to the Bivariate Statistics SRG User Guide entry for further information on their implementation.Anchor Table 2 Table 2
Objective Function | Parameter | Parameter Description | Units | Default | Range |
|---|---|---|---|---|---|
| NSE Daily & Flow Duration | a | Weight on NSE in the combined objective | Dimensionless | 0.5 | 0 ≤ α ≤ 1 |
| NSE Daily & Log Flow Duration | a | Weight on NSE |
...
| in the combined objective | Dimensionless | 0.5 | 0 ≤ α ≤ 1 |
Output data
Outputs include results of the evaluation of the selected objective function and other calibration performance statistics.
Reference list
Aitken, A.P. (1973). Assessing systematic errors in rainfall-runoff models. J. Hydrol, 20, 121–136.
Black, D.C. and Podger, G.M. (2012). Guidelines for modelling water sharing rules in eWater Source: towards best practice model application. eWater Cooperative Research Centre, Canberra, Australia. July. ISBN: 978-1-921543-74-6. Available via: www.ewater.com.au.
Black, D.C., Wallbrink, P.J., Jordan, P.W., Waters, D., Carroll, C., and Blackmore, J.M. (2011). Guidelines for water management modelling: towards best practice model application. eWater Cooperative Research Centre, Canberra, Australia. September. ISBN: 978-1-921543-46-3. Available via: www.ewater.com.au.
Duan, Q., Sorooshian, S. and Gupta, V. (1992). Effective and Efficient global optimization for conceptual rainfall-runoff models. Water Resources Research, 28(4), 1015-1031.
Lerat, J. (2012). Towards the adoption of uncertainty assessment in water resources models: the eWater Source uncertainty guideline. Proceedings of the 34th Hydrology and Water Resources Symposium, 19-22 November 2012, Sydney, NSW.
Lerat, J., Egan, C. A., Kim, S., Gooda, M., Loy, A., Shao, Q., and Petheram, C. (2013). Calibration of river models for the Flinders and Gilbert catchments. A technical report to the Australian Government from the CSIRO Flinders and Gilbert Agricultural Resource Assessment, part of the North Queensland Irrigated Agriculture Strategy. CSIRO Water for a Healthy Country and Sustainable Agriculture flagships, Australia.
Nash, J.E. and Sutcliffe, J.V. (1970). River flow forecasting through conceptual models, I, A discussion of principles. J. Hydrol, 10, 282–290.
Rosenbrock, H.H. (1960). An automated method of finding the greatest of least value of a function. The Computer Journal, 3, 303-307.
Sorooshian, S., Duan, Q. and Gupta, V. (1993). Calibration of rainfall-runoff models: application of global optimization to the Sacramento Soil Moisture Accounting Model. Water Resources Research, 29(4), 1185-1194.
Vaze, J., Jordan, P., Beecham, R., Frost, A., Summerell, G. (2011). Guidelines for rainfall-runoff modelling: Towards best practice model application. eWater Cooperative Research Centre, Canberra, ACT. ISBN 978-1-921543-51-7. Available via www.ewater.com.au.