GLO MF Model v02

## **Abstract** 

The dataset was derived by the Bioregional Assessment Programme from multiple source datasets. The source datasets are identified in the Lineage field in this metadata statement. The processes undertaken to produce this derived dataset are described in the History field in this metadata statement.



The dataset contains all the input files and scripts to run the alluvial MODFLOW models for the Avon and Karuah alluvial systems in the Gloucester subregion as described in product 2.6.2, Groundwater numerical modelling for the Gloucester subregion (Peeters et al. 2016), as well as a single instance of the model output files.



This model is run 10.000 times to provide inputs for the uncertainty analysis and the surface water model, as described in products 2.6.2 (Peeters et al. 2016) and 2.6.1 (Zhang et al. 2016)



References

Peeters L, Dawes W, Rachakonda P, Pagendam D, Singh R, Pickett T, Frery E, Marvanek S, and McVicar T (2016) Groundwater numerical modelling for the Gloucester subregion. Product 2.6.2 for the Gloucester subregion from the Northern Sydney Basin Bioregional Assessment. Department of the Environment, Bureau of Meteorology, CSIRO and Geoscience Australia, Australia. http://data.bioregionalassessments.gov.au/product/NSB/GLO/2.6.2.



Zhang Y Q, Viney N R, Peeters L J M, Wang B, Yang A, Li L T, McVicar T R, Marvanek S P, Rachakonda P K, Shi X G, Pagendam D E and Singh R M (2016) Surface water numerical modelling for the Gloucester subregion. Product 2.6.1 for the Gloucester subregion from the Northern Sydney Basin Bioregional Assessment. Department of the Environment, Bureau of Meteorology, CSIRO and Geoscience Australia, Australia., Department of the Environment, Bureau of Meteorology, CSIRO and Geoscience Australia, Australia., http://data.bioregionalassessments.gov.au/product/NSB/GLO/2.6.1.

## **Purpose** 

The 10,000 MODFLOW outputs have been combined with 10,000AWRA-L model outputs to evaluate impacts of additional coal resources development (ACRD) on streamflow at 30 simulation nodes.

## **Dataset History** 

The dataset has all files and scripts necessary to execute the 10.000 runs on the linux platform of the CSIRO High Performance Cluster computers. Powerpoint file 'GLO_MF_flowchart.pptx' describes the workflow of how all the files link together and the order in which they need to be executed.



The main script is 'GLO_MF.py'. This scripts reads a single parameter combination from file 'GLO_DoE_Parameters.csv' in dataset GLO MF dmax v01 and the outputs of the analytical element model in dataset GLO AEM Model v02 (GUID: e76e05a1-3d3c-4014-a163-66a342278ea2) to generate the various input files for a MODFLOW run. The script subsequently runs the three variants (no development, baseline, crdp) for both the Avon and Karuah models. The script finishes with copying the essential output files to the output folders and deleting the intermediate files.



The folder 'grids' has the ascii grids of the IBOUND array for the MODFLOW models (grids/GLO_Avon_IBOUND.asc and grids/GLO_Karuah_IBOUND.asc), i.e. a binary grid where the active model cells have a value of 1. The inactive cells are 0. These grids are derived from the extent of the alluvium in dataset Avon Wards River alluvium v03 (GUID: 72858ae2-da7f-4a97-8481-62aab1e3c184). These files are copied into identical files ac-cells.csv and nf-cells.csv, the IBOUND array for Avon layer 1 and 2 respectively, and km-cells.csv and l2-cells.csv, the IBOUND array for Karuah layer 1 and 2 respectively.



The folder 'grids' also has grids with the DEM values, from dataset GLO DEM 1sec SRTM MGA56, in files grids/GLO_Avon_href.dat and grids/GLO_Karuah_href.dat. These files are used as the top of the model (topo-0.txt for Avon and km-surf-0.txt for Karuah), the base of layer 1 (topo-b.txt for Avon and km-base-0.txt for Karuah), the extinction depth surface for evapotranspiration as topography minus 2m (topo-2.txt for Avon and km-surf-2.txt) and the reference surface for constant heads in layer 2 (topo-x.txt for Avon and km-surf-x.txt for Karuah). These files are referred to in MODFLOW input files gloc.evt, gloc.bas, gloc.dis and karauah.evt, karuah.bas and karauh.dis.



The outputs of the GLO AEM model, the drawdown in the weathered zone for baseline and crdp at selected output locations is interpolated to the modflow grid in GLO_MF.py and combined with the reference surface for constant heads to create the constant head files required for the MODFLOW run. The python scripts documents the interpolation algorithm and the naming of these intermediate constant head files.



File GLO_MF_rch_ts.csv has the recharge time series from the GLO AWRA model (dataset GLO AWRA Model v01 -GUID: 15ca8f9d-84b4-4395-87db-ab4ff15b9f07). This timeseries is in GLO_MF.py combined with the recharge multiplier parameter into modflow grids that can be read by the RCH MODFLOW package.



The files GLO_MF_Avon_drn.csv and GLO_MF_Karuah_drn.csv have the model cells that have drain boundary conditions. These coincide with the river network in dataset BA SYD selected GA TOPO 250K data (GUID: 	ba5feac2-b35a-4611-82da-5b6213777069) plus added map features. These files get combined with the drain conductance parameter in GLO_MF.py to create the files for the DRN MODFLOW package.



Files gloc.hob and karuah.hob ensure output time series are generated at the receptors (from dataset GLO Receptors 20150828 - GUID: acdd94a6-5201-4026-8177-035f25077587).



Files gloc.drob and karuah.drob ensure output time series are generated for each of river reaches of the surface water model (from dataset GLO Receptors 20150828 - GUID: acdd94a6-5201-4026-8177-035f25077587).



Files gloc.oc and karuah.oc are the output control files for modflow.



Files gloc.sip and karuah.sip are the solver setting files for modflow.



For each of the model runs, an individual NAM file is generated by GLO_MF.py.



File mf2005 is the modflow 2005 v1.11.00 executable compiled for linux, based on the source files available from the USGS website (http://water.usgs.gov/ogw/modflow/MODFLOW.html).



After the modflow run is completed the relevant outputs are copied into the outputs folder. The .lst files are the list-files of modflow, which contain the water balance of the entire model at each time step. The .hout files have the timeseries of groundwater levels at each receptor. The .drout files have the timeseries of drainage flux at each of the river reaches used in AWRA-L.



The batch file GLO_MF.q runs the script GLO_MF.py for each parameter combination in GLO_DoE_Parameters.csv on the Linux high performance cluster computers at CSIRO.



Python script GLO_MF_postprocessing.py, which is run through batch file GLO_MF_postprocessing.q, creates the maximum drawdown (dmax) and time to maximum drawdown (tmax) at the receptor locations for Avon and Karuah (files outputs/dmax/GLO_MF_dmax_g_DoE_Predictions.csv, outputs/dmax/GLO_MF_dmax_k_DoE_Predictions.csv, outputs/dmax/GLO_MF_tmax_g_DoE_Predictions.csv, outputs/dmax/GLO_MF_tmax_k_DoE_Predictions.csv) as well as the estimates of the water balance used to constrain the MODFLOW models (outputs/budget/GLO_MF_gloc_budget_m3d.csv and outputs/budget/GLO_MF_gloc_budget_m3d.csv). The last set of outputs generated by the postprocessing script are the drainage fluxes that are used in the AWRA-L model: outputs/baseflow/GLO_MF_drup_base_m3d.csv, outputs/baseflow/GLO_MF_drup_crdp_m3d.csv and outputs/baseflow/GLO_MF_drup_ndev_m3d.csv

## **Dataset Citation** 

Bioregional Assessment Programme (2016) GLO MF Model v02. Bioregional Assessment Derived Dataset. Viewed 18 July 2018, http://data.bioregionalassessments.gov.au/dataset/9d732408-003e-4901-86a7-cea95b585640.

## **Dataset Ancestors** 

* **Derived From** [Standard Instrument Local Environmental Plan (LEP) - Heritage (HER) (NSW)](https://data.gov.au/data/dataset/fd1ff106-2f48-4bda-bfd1-a7206747c483)

* **Derived From** [NSW Office of Water GW licence extract linked to spatial locations - GLO v5 UID elements 27032014](https://data.gov.au/data/dataset/0115c2ba-73c6-4c98-b539-9f67594980cf)

* **Derived From** [Gloucester digitised coal mine boundaries](https://data.gov.au/data/dataset/be85b60f-49dd-485c-bb85-11d45de9792d)

* **Derived From** [Groundwater Dependent Ecosystems supplied by the NSW Office of Water on 13/05/2014](https://data.gov.au/data/dataset/1bdc6089-c114-4485-b962-35a8b3d2357f)

* **Derived From** [Geology 100K NSW selected mapsheets for Sydney Bioregion](https://data.gov.au/data/dataset/19a28cab-2acb-4c9d-b265-1535158508f2)

* **Derived From** [GLO Geological Model Extracted Horizons Final Grid XYZ V01](https://data.gov.au/data/dataset/18402b28-677e-45af-9860-067c258eea33)

* **Derived From** [NSW Office of Water GW licence extract linked to spatial locations GLOv4 UID 14032014](https://data.gov.au/data/dataset/98341388-5705-4c42-9f7a-186e16647614)

* **Derived From** [Communities of National Environmental Significance Database - RESTRICTED - Metadata only](https://data.gov.au/data/dataset/c01c4693-0a51-4dbc-bbbd-7a07952aa5f6)

* **Derived From** [National Groundwater Dependent Ecosystems (GDE) Atlas](https://data.gov.au/data/dataset/e358e0c8-7b83-4179-b321-3b4b70df857d)

* **Derived From** [Asset database for the Gloucester subregion on 12 September 2014](https://data.gov.au/data/dataset/d0c21c5c-1b30-4aae-933a-71145b935cc6)

* **Derived From** [GEODATA 9 second DEM and D8: Digital Elevation Model Version 3 and Flow Direction Grid 2008](https://data.gov.au/data/dataset/ebcf6ca2-513a-4ec7-9323-73508c5d7b93)

* **Derived From** [BA SYD GLO -  AGL Avon Wards River alluvium extent v1](https://data.gov.au/data/dataset/4484f125-8401-4324-a8a4-42fd1828e485)

* **Derived From** [National Groundwater Information System (NGIS) v1.1](https://data.gov.au/data/dataset/e157b1cd-3cb7-4cf9-a13c-a6bed576b8c7)

* **Derived From** [Groundwater Entitlement Data GLO NSW Office of Water 20150320 PersRemoved](https://data.gov.au/data/dataset/5afdddd5-f5b5-42d1-8ca6-8db19ecf1fd3)

* **Derived From** [GLO MFdmax v01](https://data.gov.au/data/dataset/df0ff3b6-b131-4cf7-9d86-33d1de01634d)

* **Derived From** [R-scripts for uncertainty analysis v01](https://data.gov.au/data/dataset/322c38ef-272f-4e77-964c-a14259abe9cf)

* **Derived From** [New South Wales 2 kilometers Residential Exclusions Zone](https://data.gov.au/data/dataset/8e9b4d01-bba2-4741-9ffd-aed0484eb14a)

* **Derived From** [Geofabric Surface Cartography - V2.1](https://data.gov.au/data/dataset/5342c4ba-f094-4ac5-a65d-071ff5c642bc)

* **Derived From** [BA SYD selected GA TOPO 250K data plus added map features](https://data.gov.au/data/dataset/ba5feac2-b35a-4611-82da-5b6213777069)

* **Derived From** [Groundwater Entitlement Data Gloucester - NSW Office of Water 20150320](https://data.gov.au/data/dataset/bf435e81-60d4-40c7-ac00-76585fc1fe77)

* **Derived From** [Collaborative Australian Protected Areas Database (CAPAD) 2010 - External Restricted](https://data.gov.au/data/dataset/47312aee-722e-4c6e-bef8-9e439480503e)

* **Derived From** [National Groundwater Dependent Ecosystems (GDE) Atlas (including WA)](https://data.gov.au/data/dataset/6dbaee0d-8813-46b1-9c13-1b796e7ed3bf)

* **Derived From** [EIS Gloucester Coal 2010](https://data.gov.au/data/dataset/04c22c59-2f7c-4ba6-a192-7bfb89f6b98d)

* **Derived From** [Report for Director Generals Requirement Rocky Hill Project 2012](https://data.gov.au/data/dataset/d8837d3d-91bb-4496-814e-e707cf65fbb2)

* **Derived From** [Species Profile and Threats Database (SPRAT) - Australia - Species of National Environmental Significance Database (BA subset - RESTRICTED - Metadata only)](https://data.gov.au/data/dataset/7276dd93-cc8c-4c01-8df0-cef743c72112)

* **Derived From** [Geological Maps Combined for NSW](https://data.gov.au/data/dataset/f507bb28-8095-43f8-901e-565699a290b5)

* **Derived From** [GEODATA TOPO 250K Series 3](https://data.gov.au/data/dataset/a0650f18-518a-4b99-a553-44f82f28bb5f)

* **Derived From** [Asset database for the Gloucester subregion on 28 May 2015](https://data.gov.au/data/dataset/5815842e-d271-4f73-9d1a-d15c90571330)

* **Derived From** [Gloucester Deep Wells Completion Reports - Geology](https://data.gov.au/data/dataset/0529ae9a-4d40-460d-b52a-fb0e5f646746)

* **Derived From** [NSW Office of Water GW licence extract linked to spatial locations GLOv3 12032014](https://data.gov.au/data/dataset/587c6834-09aa-4cd8-8972-5632880f5ae5)

* **Derived From** [EIS for Rocky Hill Coal Project 2013](https://data.gov.au/data/dataset/ebc1abb5-8b0e-43fa-95e3-2c49dd3343d7)

* **Derived From** [GLO AEM dmax v01](https://data.gov.au/data/dataset/c54640db-ca88-4ed2-8e58-6b1a198293c5)

* **Derived From** [National Heritage List Spatial Database (NHL) (v2.1)](https://data.gov.au/data/dataset/26daa8d7-a90e-47f3-982b-0df362414e65)

* **Derived From** [Asset database for the Gloucester subregion on 8 April 2015](https://data.gov.au/data/dataset/d559e708-089a-46bc-bad2-dc5068993574)

* **Derived From** [Gloucester - Additional assets from local councils](https://data.gov.au/data/dataset/47e9540a-dd75-4e04-bddb-574de1b89ecb)

* **Derived From** [NSW Office of Water combined geodatabase of regulated rivers and water sharing plan regions](https://data.gov.au/data/dataset/24157c41-c42f-4e1f-a791-a1ad18c8215d)

* **Derived From** [GLO RMS Model Depth Structure Eroded v01](https://data.gov.au/data/dataset/d2d676b1-26e5-4b5f-bcf3-60ae7fe41cec)

* **Derived From** [Asset database for the Gloucester subregion on 29 August 2014](https://data.gov.au/data/dataset/7000a82c-85f2-461a-ba35-9d7d1ba6d339)

* **Derived From** [New South Wales NSW Regional CMA Water Asset Information WAIT tool databases, RESTRICTED Includes ALL Reports](https://data.gov.au/data/dataset/a2a929ea-2cb3-42ab-86b9-87deac68acae)

* **Derived From** [Groundwater Modelling Report for Stratford Coal Mine](https://data.gov.au/data/dataset/4a9d9e6c-8d12-4d9c-837f-5c63fa03472a)

* **Derived From** [AGL Gloucester Gas Project AECOM report location map features](https://data.gov.au/data/dataset/7fca4f35-6e4e-4d37-bbae-1c1029a93a40)

* **Derived From** [Groundwater Economic Assets GLO 20150326](https://data.gov.au/data/dataset/2e314212-0677-40b8-86ff-c5166c6906bd)

* **Derived From** [NSW Office of Water Groundwater Licence Extract Gloucester - Oct 2013](https://data.gov.au/data/dataset/56a49393-abaa-43f2-8b7d-97c5d414c1cf)

* **Derived From** [New South Wales NSW - Regional - CMA - Water Asset Information Tool - WAIT - databases](https://data.gov.au/data/dataset/330532aa-66ba-44f5-984b-8a21a99661a0)

* **Derived From** [GLO Deep Well Locations and Depths of Formations V01](https://data.gov.au/data/dataset/e549473c-5b78-4279-8eeb-af86ca6522b3)

* **Derived From** [Freshwater Fish Biodiversity Hotspots](https://data.gov.au/data/dataset/6093bd55-a978-4fcc-a8fa-6881fa015329)

* **Derived From** [NSW Office of Water Groundwater licence extract linked to spatial locations GLOv2 19022014](https://data.gov.au/data/dataset/38c9e2c9-02ac-45a8-a188-936f91a58391)

* **Derived From** [GLO AEM Model v02](https://data.gov.au/data/dataset/e76e05a1-3d3c-4014-a163-66a342278ea2)

* **Derived From** [Australia - Species of National Environmental Significance Database](https://data.gov.au/data/dataset/4eff885b-0d68-45e8-b4b1-e3e5448354eb)

* **Derived From** [GLO DEM 1sec SRTM MGA56](https://data.gov.au/data/dataset/ca38ed31-e15d-4bb5-a7ef-0aeba3dad3f4)

* **Derived From** [Australia, Register of the National Estate (RNE) - Spatial Database (RNESDB) Internal](https://data.gov.au/data/dataset/878f6780-be97-469b-8517-54bd12a407d0)

* **Derived From** [Avon Wards River alluvium v03](https://data.gov.au/data/dataset/72858ae2-da7f-4a97-8481-62aab1e3c184)

* **Derived From** [NSW Office of Water Groundwater Entitlements Spatial Locations](https://data.gov.au/data/dataset/0f4da0d0-9e82-42f3-84da-1148e0b730c1)

* **Derived From** [GLO Receptors 20150828](https://data.gov.au/data/dataset/acdd94a6-5201-4026-8177-035f25077587)

* **Derived From** [Directory of Important Wetlands in Australia (DIWA) Spatial Database (Public)](https://data.gov.au/data/dataset/6636846e-e330-4110-afbb-7b89491fe567)

* **Derived From** [Geoscience Australia, 1 second SRTM Digital Elevation Model (DEM)](https://data.gov.au/data/dataset/9a9284b6-eb45-4a13-97d0-91bf25f1187b)

* **Derived From** [Collaborative Australian Protected Areas Database (CAPAD) 2010 (Not current release)](https://data.gov.au/data/dataset/7b649c6d-fdbf-40e3-b002-db521665af53)



Tags: Gloucester subregion, environment, geoscientificInformation, inlandWaters
Organization: Bioregional Assessment Program