All major rivers (54 reaches) within the Figure 14.(Geoscience Australia, ) are represented in the model. Some additional small reaches were added to ensure that change in – flux could be generated along rivers represented in the AWRA-R node-link network. The polylines used to define the locations of rivers in the groundwater model are shown in
Also shown are the Australian Water Resources Assessment river model (AWRA-R) catchments (red lines) and the surface water nodes (yellow dots) where the change in surface water – groundwater flux is aggregated.
Data: Bioregional Assessment Programme ()
The River (RIV) package for MODFLOW-USG uses a series of model cells to represent the river boundary condition. Close to the CSG and coal mining developments the model was discretised as 300 m cells. In these regions, the river was decomposed into a sequence of points, spaced at 300 m intervals. This means surface water – groundwater flux is recorded only at 300 m intervals, rather than continuously along a river reach. The parameters necessary for the RIV package are stage height, conductance and river-bed geometry.
A steady-state river stage for the first stress period of the model was implemented based on average river stages for the historical period (1983–2012). The steady-state stages were interpolated from measurements obtained at 54 gauge sites located within the model domain (Figure 14). Transient river stages at the 54 gauges were derived from rating curves based on historical records and the flow volumes from the surface water modelling.
Another key controlling parameter is the, CRIV, which is defined as follows:
where Kn and Mn represent the hydraulic conductivity and thickness of the riverbed, respectively, and Ln and Wn represent the length and width of the river reach at node n, respectively (220.127.116.11.). Due to a lack of data in the Namoi subregion, hydraulic conductivity of the riverbed was initially assigned values that range from 1e–6 to 1e0 m/day. Since the gaining and losing nature of a river depends largely on the relative elevation of the river stage with respect to the groundwater level in the , depth of the riverbed below the land surface was identified as the most important parameter that influences the surface water – groundwater in the Namoi subregion. This was further tested in the groundwater model stress-test and is reported in Section
It is widely accepted that riverbed hydraulic conductivity is primarily a function of reach geometry, streamflow velocity, composition and erodibility of catchment, and bed disturbance frequency (). The riverbed hydraulic conductivity generally increases with riverbed slope when other factors are similar ( ). The current study explored the conceptualisation of hydraulic conductivity variation in direct relationship with the riverbed slope. However, it was found that the surface water – groundwater interaction was mostly sensitive to the depth of the river bottom in relation to the groundwater level in the watertable aquifer.
Depth of the river bottom was considered to vary between 2 to 15 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 3 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 riverbed’s vertical position downwards in relation to the land surface in the model. This shift is varied in the analysis.
18.104.22.168.3.2 AWRA-R baseflows
The river model contains 54 nodes at which themodel can provide change in – groundwater flux 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 under the . The (positive) and leakages (negative) are summed for all the points in the link upstream of a node in order to calculate the total change in surface water – groundwater flux 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 change in surface water – groundwater flux. The river model runs using daily time steps, and the groundwater model results are interpolated linearly to provide this.
Product Finalisation date
- 22.214.171.124 Methods
- 126.96.36.199 Review of existing models
- 188.8.131.52 Model development
- 184.108.40.206 Boundary and initial conditions
- 220.127.116.11 Implementation of the coal resource development pathway
- 18.104.22.168 Parameterisation
- 22.214.171.124 Observations and predictions
- 126.96.36.199 Uncertainty analysis
- 188.8.131.52 Limitations
- Currency of scientific results
- Contributors to the Technical Programme
- About this technical product