The focus of the numerical modelling for the Clarence-Moreton Bioregional Assessment is the coal seam gas (CSG) development in the Walloon Coal Measures in the Richmond river basin (Figure 3). The Richmond river basin overlies a sedimentary basin with an alternating sequence of and (Figure 4) in which the alluvial deposits of the Richmond River and its tributaries are incised. The Walloon Coal Measures outcrop in the west of the study area. In the proposed CSG development area close to Casino, the Walloon Coal Measures are overlain by hundreds of metres of younger sedimentary bedrock, including several hydrostratigraphic units that are considered to be aquitards . In the northern part of the Richmond river basin, the Lamington Volcanics form a topographic high, and these high-elevation basalts are associated with very high rates.
The Richmond River system is perennial, with the exception of some smaller tributaries in the west such as Shannon Brook. Unlike the perennial streams, these smaller tributaries in the western part of the Richmond river basin are not incised into the Lamington Volcanics, and flow in Shannon Brook, for instance, is observed 94% of the year. Downstream from Casino, the Richmond River is tidally influenced.
Viewed from the south-east; the vertical extent is from –2500 to +1400 m Australian Height Datum (AHD); the north–south extent is 320 km; the maximum east–west extent is 140 km; the vertical exaggeration is 10
The main group identified in companion product 2.3 for the Clarence-Moreton bioregion considered in the numerical modelling is the ‘Subsurface depressurisation and dewatering’ causal pathway group. Dewatering applies to coal mines and, given that no coal mines are being modelled in the or , it will not be discussed further in relation to the . Subsurface depressurisation of the Walloon Coal Measures near Casino, may lead to in the Walloon Coal Measures and overlying aquifers including the alluvial aquifers. A change in groundwater pressures and levels in the alluvium may result in changes to – groundwater interactions and thus a reduction of streamflow in the Richmond River or its tributaries.
Different model types and model codes have been adopted in (BAs) depending on the specific requirements of each or . The companion submethodology M07 (as listed in Table 1) for groundwater modelling stipulates that a groundwater model must deliver spatially explicit outputs that can evaluate the / on groundwater resources as well as providing input to other modelling tasks, including surface water modelling, analysis and impact modelling. A formal assessment of the modelling approach adopted in the Clarence-Moreton bioregion was conducted according to the five fit-for-purpose criteria listed in Table 3.
Table 3 Assessment of groundwater numerical modelling approach in the Clarence-Moreton bioregion
18.104.22.168.2.1 Prediction of hydrological response variables
The objective of modelling that is undertaken as part of a is to probabilistically assess the of at associated with . The is quantified by the change in some predefined such as fluxes, hydraulic heads or variables depending on fluxes or hydraulic heads. Hydrological response variables need to be decided upon prior to the analysis, for example, at model node (x, y, z) at time t. The hydrological response variables adopted in the current assessment include additional drawdown (dmax, maximum difference in drawdown between the (CRDP) and , due to the additional coal resource development), year of maximum change (tmax) and leakage/baseflow changes at a number of model nodes.
In order to quantify the uncertainties associated with model predictions, groundwater modelling in BA is run probabilistically and not deterministically. Consequently, this means that model predictions cannot be presented as unique values but as . Mismatches in spatial and temporal scale between the regional nature of the modelling and the point-scale nature of model nodes mean that the modelling is not able to capture fine-scale complexities of impacts upon model nodes. In BA, the focus is on the impact of the additional coal resource development, which is attributed to the difference in results between two possible futures – the baseline and the . Hence, numerical groundwater modelling results are primarily presented as the difference in groundwater drawdown between the baseline and CRDP. This approach reduces the as the prediction is formulated as the difference between the output of two model results .
The linkage between the groundwater model and the Australian Water Resources Assessment (AWRA) landscape model (AWRA-L) is depicted in Figure 5. AWRA-L was used to generate the temporal trends of diffuse for the groundwater model. The river stage elevation time series for the groundwater model was derived from the outputs of AWRA-L, which were coupled with rating curves derived from measured data sourced from the NSW streamflow from NSW Office of Water. The groundwater model was used to generate the leakage/baseflow change at a number of gauges along the Richmond River due to the additional coal resource development. These flux changes can then be fed back into the AWRA-L model predictions to estimate changes of relevant hydrological response variables. Subsequently, changes of hydrological response variables can be used to assess the impacts on model nodes associated with
Although the groundwater model was able to deliver input for impact modelling, receptor impact modelling was not considered in this current BA. Analysis of the potential impacts of the additional coal resource development have shown minimal hydrological changes, particularly at the uppermost where median additional drawdowns does not exceed 0.2 m. The median predicted change in surface water – groundwater flux at simulated model nodes, including both leakage and baseflow, does not exceed 0.01 GL/year as shown in companion product 2.6.1 for the . As a result of these minor hydrological impacts, ecological impacts are expected to be below the detectable limit and hence the decision was to not carry out receptor impact modelling in the Clarence-Moreton bioregion.
AWRA-L = Australian Water Resources Assessment landscape module; baseline = baseline coal resource development; CRDP = coal resource development pathway; CRDP = baseline + additional coal resource development; green rectangles represent models; orange rectangles are input and/or output data of models
22.214.171.124.2.2 Design and construction
Following the guidelines of companion submethodology M07 (as listed in Table 1) for modelling , a numerical groundwater model was developed to assess the potential impacts of the on the of the . As reported in companion product 2.3 for the Clarence-Moreton bioregion , the additional coal resource development in this only includes one CSG development, the West Casino Gas Project. The groundwater model mainly focused on the Richmond river basin where the West Casino Gas Project is likely to occur. The conceptualisation and the likely that underpinned the model are identified in companion product 2.3 for the Clarence-Moreton bioregion .
The major assumptions and model choices related to model development are described in Section 126.96.36.199. The groundwater model is considered the first regional model for the area and is based on a three-dimensional geological model that was specifically developed for the current BA (see companion product 2.3 for the Clarence-Moreton bioregion ). MODFLOW-NWT was used to run the model in order to minimise the numerical instabilities that can be caused by cell drying and rewetting in the unconfined layer of the model. The effect of the assumptions and model choices on predictions are discussed in Section 188.8.131.52.
Groundwater modelling simulations were run under historical conditions for the 30-year period from 1 January 1983 to 31 December 2012. The historical climate time series was repeated three times to create a 90-year time series and modified to be consistent with a median future climate projection as described in companion submethodology M06 (as listed in Table 1) for modelling . This 30-year period was repeated to ensure that the effect of droughts and floods does not confound the comparison between time periods. Monthly stress periods were used throughout the simulation to provide monthly change of baseflow required by the surface water model. This led to a total number of 1441 stress periods from 1983 to 2102, including the first steady-state stress period.
Final outputs of predicted were produced with the aid of emulators as described in companion submethodology M07 (as listed in Table 1) for groundwater modelling and M09 (as listed in Table 1) for propagation of through models . The main goal of the groundwater model was to provide realisations that can be used to train the emulators for the . Emulators are computationally efficient approximation statistical tools compared to the slower process-based groundwater model. As described in companion submethodology M09 , the application of model emulation is to overcome the computational inefficiency of the original groundwater model, but still being able to capture the important groundwater flow features.
For those model nodes that are associated with groundwater-dependent assets and hence directly impacted by groundwater , such as stock and domestic groundwater , the groundwater model directly generates the impacts. If impact modelling is deemed necessary, then the groundwater model outputs will directly feed into the receptor impact model. For other model nodes that are indirectly impacted by groundwater extraction through an interaction with any of the various surface water features, a surface water model is required to simulate the impact of the additional coal resource development. Due to minimal regulation in the Richmond river basin, the rainfall- model AWRA-L was deemed adequate to model the Richmond River flows. The modelling results for these model nodes are presented in companion product 2.6.1 for the Clarence-Moreton bioregion .
184.108.40.206.2.3 Integration with sensitivity and uncertainty analysis workflow
ABC MCMC = Approximate Bayesian Computation Markov chain Monte Carlo; HRV = hydrological response variable
Submethodology M09 (as listed in Table 1) discusses in detail the propagation of through numerical models in the . Figure 6 summarises the uncertainty propagation workflow which consists of four major steps:
- Design of experiment: large number of evaluations with a wide range of parameter values
- Train emulators for:
- Create posterior parameter through Approximate Bayesian Computing Markov chain Monte Carlo
- Sample the posterior parameter probability distribution to generate the posterior probability distribution for each hydrological response variable at each model node.
The first step is to identify the parameters of the model chain to include in the uncertainty analysis and to define a wide range that represents the plausible range of the parameters. A large number of model chain evaluations are carried out, sampling extensively from this parameter range. For each evaluation the corresponding simulated change in hydrological response variables at the model nodes is stored, together with the simulated equivalents to the observations. The latter are summarised into objective functions, tailored to each hydrological response variable.
This information forms the basis for the subsequent uncertainty analysis. In the uncertainty analysis, the prior parameter distributions, the most likely range of the parameter values based on data and expert knowledge, is constrained with the available relevant data using the Approximate Bayesian Computation methodology. This results in a posterior parameter distribution, tailored to a specific hydrological response variable, which subsequently can be sampled to generate a probability distribution at each model node.
This type of uncertainty analysis requires a very large number of model evaluations, which is practically not feasible. This is the main reason that the original model chain in the uncertainty analysis is replaced by emulators, statistical functions that closely mimic the effect of parameter values on predictions. These emulators take little time to evaluate and are straightforward to integrate into the uncertainty analysis workflow.
In order for the model chain to be amenable for incorporation into this uncertainty analysis it needs to be scripted so that parameter values can be changed in an automated fashion, be able to be evaluated from a command line on high performance computers and, most importantly, be numerically stable so that the model converges for a wide range of parameter values.
The two models in the model chain for the have text files as input files and can be executed from the command line. The robustness of each model is tested through a stress-test in which a selection of extreme parameter combinations is evaluated. While this does not guarantee that all model evaluations will converge, it provides confidence that the majority of parameter combinations will.
Section 220.127.116.11 and Section 18.104.22.168 provide details of the implementation of this uncertainty propagation workflow for the MODFLOW model. The uncertainty analysis for AWRA-L is in Section 22.214.171.124 and Section 126.96.36.199 of companion product 2.6.1 for the Clarence-Moreton bioregion . These sections also have a qualitative uncertainty analysis that provides a structured discussion of the assumptions and model choices not included in the numerical uncertainty analysis and the perceived effect on the predictions.
188.8.131.52.2.4 Water balance components
The water balance was reported for a defined control volume in that includes all hydrologically connected changes predicted by the and models. The water balance components (e.g. , evapotranspiration, baseflow, licensed , upward flow from deeper groundwater and change in storage) were compared with estimates described in the regional-scale and localised groundwater models to provide confidence in model predictions, thus meeting the modelling guiding principles . Companion product 2.5 for the compares model estimates of water extraction from CSG with available local and regional estimates. Parameters associated with water balance components are discussed in Section 184.108.40.206 and Section 220.127.116.11.
18.104.22.168.2.5 Transparent and reproducible model outputs
The requirement that model results must be reproducible means that the models need to be run as part of a documented workflow that records the provenance of the input data, executables and outputs. This is achieved through the use of scripting in BAs. All pre-processing, model runs and post-processing is done using scripts that are made available along with the products at www.bioregionalassessments.gov.au; this ensures that all model inputs, parameters, executables and outputs are traceable, meeting the modelling guiding principles related to .
The current model described in this product is designed specifically for delivering a probabilistic assessment of the of the in the . In its current form, the model is neither suited to address any other water management questions nor provide deterministic predictions in this bioregion. The validity of all assumptions in this product should be re-assessed for any other application of the current model.
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 and conclusions
- Contributors to the Technical Programme
- About this technical product