All major rivers (31 named reaches) within the (Geoscience Australia, ) are represented in the model. Some additional small reaches were added to ensure that change in baseflows could be generated along rivers represented in the AWRA-R (river) node-link network. The polylines used to define the locations of rivers in the model are shown in Figure 16.
The Multiphysics Object-Oriented Simulation Environment (MOOSE) software used to implement and run the model cannot represent line sources and sinks, only point sources and sinks, so each river was decomposed into a sequence of points, spaced at 1 km intervals. This means baseflow is recorded only at 1 km intervals, rather than continuously along a river reach.
These points are placed at a depth of 5 m below the model’s topography. This is for two reasons. Firstly, rivers tend to incise channels below the land surface at a scale that may be too fine to be represented in the 3-second digital elevation model (Geoscience Australia, ) used to define the surface topography. At some points on the river network, surveyed information has been recorded that could be used instead, but this is not available for most of the river network. Secondly, the model further discretises this digital elevation model with resolution as coarse as 15 km, so many points at their true elevation according to the digital elevation model would lie outside the model domain (above the model topography). For these two reasons it is appropriate to shift the points’ vertical position downwards. This shift is varied in the analysis.
184.108.40.206.3.2 Baseflow and leakage
At each point, representing a 1 km section of river, along the river network, the model can lose water to represent fluxes to the river (i.e. baseflow) or gain water, representing leakage from the river, according to:
- F is the flux out of the (in kg/year) and is positive for baseflow and negative for leakage
- h is the level at the point
- h0 is the river stage height
- C is the , which has a base value of 320 ML/m/year, but which is a parameter varied in the analysis
- T is a parameter that limits the leakage. For perennial rivers T = –1 m, while for ephemeral rivers T = 0. This means there is no leakage for ephemeral rivers, while the leakage from perennial rivers does not increase once the water head is 1 m below h0.
The base value of C is motivated by assuming the riverbed has a width of order 10 m and thickness of order 1 m and conductivity of 0.1 m/day. The river stage height defaults to h0 = 3 m, for all points. However, because the depth of the points is varied in the uncertainty analysis, this effectively also varies h0.
220.127.116.11.3.3 AWRA-R baseflows
The model contains 63 nodes at which the model can provide change in baseflow estimates. Nodes mostly correspond to streamflow gauging stations, but some nodes have been included specifically for assessing hydrological changes in response to coal resource development under the and (CRDP). The baseflows (positive) and leakages (negative) are summed for all the points in the link upstream of a node in order to predict the total baseflow at the node. The numerical implementation of the groundwater model runs using nominal monthly time steps, so can only provide a monthly time series of gauge baseflows. The surface water model runs using daily time steps, and the groundwater results are interpolated linearly to provide this.
Product Finalisation date
- 18.104.22.168 Methods
- 22.214.171.124 Review of existing models
- 126.96.36.199 Model development
- 188.8.131.52 Boundary and initial conditions
- 184.108.40.206 Implementation of the coal resource development pathway
- 220.127.116.11 Parameterisation
- 18.104.22.168 Observations and predictions
- 22.214.171.124 Uncertainty analysis
- 126.96.36.199 Limitations and conclusions
- Currency of scientific results
- Contributors to the Technical Programme
- About this technical product