Namoi groundwater uncertainty analysis

## **Abstract** \n\nThis 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.\n\n\n\nThe dataset contains the predictions of maximum drawdown and time to maximum drawdown at all groundwater model nodes in the Namoi subregion, constrained by the observations of groundwater level, river flux and mine water production rates. The dataset also contains the scripts required for and the results of the sensitivity analysis. The dataset contains all the scripts to generate these results from the outputs of the groundwater model (Namoi groundwater model dataset) and all the spreadsheets with the results. The methodology and results are described in Janardhanan et al. (2017)\n\n\n\nReferences\n\nJanardhanan S, Crosbie R, Pickett T, Cui T, Peeters L, Slatter E, Northey J, Merrin LE, Davies P, Miotlinski K, Schmid W and Herr A (2017) Groundwater numerical modelling for the Namoi subregion. Product 2.6.2 for the Namoi subregion from the Northern Inland Catchments Bioregional Assessment. Department of the Environment and Energy, Bureau of Meteorology, CSIRO and Geoscience Australia, Australia., http://data.bioregionalassessments.gov.au/product/NIC/NAM/2.6.2.\n\n## **Dataset History** \n\nThe workflow that underpins this dataset is captured in 'NAM_MF_UA_workflow.png'.\n\n\n\nSpreadsheet NAM_MF_dmax_Predictions_all.csv is sourced from dataset 'Namoi groundwater model' and contains the name, coordinates, Bore_ID in the model, layer number, the name of the objective function and the minimum, maximum, median, 5th percentile and 95th percentile of the design of experiment runs of maximum drawdown (dmax) for each groundwater model node. The individual results for each node for each run of the design of experiment is stored in spreadsheet 'NAM_MF_dmax_DoE_Predictions_all.csv' The equivalent files for time to maximum drawdown (tmax) are 'NAM_MF_tmax_Predictions_all.csv' and 'NAM_MF_tmax_DoE_Predictions_all.csv'.\n\n\n\nThese files are combined with the file 'NAM_MF_Observations_all.csv', which contains the observed values for groundwater levels, mine dewatering rates and river flux, and the files NAM_MF_dist_hobs.csv, NAM_MF_dist_rivers.csv, NAM_MF_dist_mines.csv, which contain the distances of the predictions to each mine, groundwater level observation and river, in python script 'NAM_MF_datawranling.csv'. This script selects only those predictions where the 95th percentile of dmax is less than 1 cm for further analysis. The subset of predictions is stored in 'NAM_MF_dmax_Predictions.csv','NAM_MF_tmax_Predictions.csv', 'NAM_MF_dmax_DoE_Predictions.csv','NAM_MF_tmax_DoE_Predictions.csv'. The output spreadsheet 'NAM_MF_Observations.csv' has the observations and the distances to the selected predictions.\n\n\n\nAs the simulated equivalents to the observations are part of the predictions dataset, these files are combined in python script NAM_MF_OFs.py to generate the objective function values for each run and each prediction. The objective function values are weighted sums of the residuals, stored in NAM_MF_DoE_hres.csv, NAM_MF_DoE_mres.csv, NAM_MF_DoE_rres.csv, according to the distance to the predictions and the results are stored in NAM_MF_DoE_OFh.csv, NAM_MF_DoE_OFm.csv, NAM_MF_DoE_OFr.csv. The threshold values for each objective function and prediction are stored in NAM_MF_OF_thresholds.csv. Python script NAM_MF_OF_wrangling.py further post-processes this information to generate the acceptance rates, saved in spreadsheet NAM_MF_dmax_Predictions_ARs.csv\n\n\n\nPython script NAM_MF_CreatePosterior.py selects the results from the design of experiment run that satisfy the acceptance criteria. The results form the posterior predictive distributions stored in NAM_MF_dmax_Posterior.csv and NAM_MF_tmax_Posterior.csv. These are further summarised in NAM_MF_Predictions_summary.csv.\n\n\n\nThe sensitivity analysis is done with script NAM_MF_SI.py, which uses the results of the design of experiment together with the parameter values, stored in NAM_MF_DoE_Parameters.csv and their description (name, range, transform) in NAM_MF_Parameters.csv. The resulting sensitivity indices for dmax, tmax and river, head and minewater flow observations are stored in NAM_MF_SI_dmax.csv, NAM_MF_SI_tmax.csv, NAM_MF_SI_river.csv, NAM_MF_SI_mine.csv and NAM_MF_SI_head.csv. The intermediate files, ending in xxxx, are the results grouped per 100 predictions. The scripts NAM_MF_SI_collate.py and NAM_MF_SI_collate.slurm collate these.\n\n## **Dataset Citation** \n\nBioregional Assessment Programme (2017) Namoi groundwater uncertainty analysis. Bioregional Assessment Derived Dataset. Viewed 11 December 2018, http://data.bioregionalassessments.gov.au/dataset/36bd27e9-58d2-4bf2-8e4a-54b22ac98cfb.\n\n## **Dataset Ancestors** \n\n* **Derived From** [NSW Office of Water GW licence extract linked to spatial locations NIC v2 (28 February 2014)](https://data.gov.au/data/dataset/505e8a67-90b4-408a-aa0d-038a9fb14f42)\n\n* **Derived From** [Namoi hydraulic conductivity measurements](https://data.gov.au/data/dataset/5f88517d-8154-411d-907f-4e2c2d12a912)\n\n* **Derived From** [Namoi NGIS Bore analysis for 2012](https://data.gov.au/data/dataset/1c828f1b-3e9c-46ae-983a-6ea2218dba48)\n\n* **Derived From** [Namoi groundwater model alluvium extent](https://data.gov.au/data/dataset/e7e79c42-89aa-4278-b14c-d7c67ebd4d89)\n\n* **Derived From** [Surface Geology of Australia, 1:1 000 000 scale, 2012 edition](https://data.gov.au/data/dataset/8284767e-b5b1-4d8b-b8e6-b334fa972611)\n\n* **Derived From** [Namoi Leapfrog geological model](https://data.gov.au/data/dataset/26747362-8e20-49df-ab2c-918fba839aa4)\n\n* **Derived From** [Historical Mining Footprints DTIRIS NAM 20150914](https://data.gov.au/data/dataset/fb0a43f1-9022-4844-b5b5-b91f72a1d9f9)\n\n* **Derived From** [Gippsland Project boundary](https://data.gov.au/data/dataset/27413de5-d13a-4231-ac79-fc77f4cbb5f7)\n\n* **Derived From** [Bioregional Assessment areas v04](https://data.gov.au/data/dataset/e8a2d577-c5c5-4e2c-b0fa-53e4b2d4a034)\n\n* **Derived From** [Natural Resource Management (NRM) Regions 2010](https://data.gov.au/data/dataset/1d54e38f-4051-4f0c-a350-c7dbd8eba65b)\n\n* **Derived From** [Soil and Landscape Grid National Soil Attribute Maps - Clay 3 resolution - Release 1](https://data.gov.au/data/dataset/f8640540-4bb7-42ee-995a-219881e67705)\n\n* **Derived From** [GEODATA TOPO 250K Series 3, File Geodatabase format (.gdb)](https://data.gov.au/data/dataset/96ebf889-f726-4967-9964-714fb57d679b)\n\n* **Derived From** [Bioregional_Assessment_Programme_Catchment Scale Land Use of Australia - 2014](https://data.gov.au/data/dataset/6f72f73c-8a61-4ae9-b8b5-3f67ec918826)\n\n* **Derived From** [GEODATA TOPO 250K Series 3](https://data.gov.au/data/dataset/a0650f18-518a-4b99-a553-44f82f28bb5f)\n\n* **Derived From** [NSW Office of Water Groundwater Licence Extract NIC- Oct 2013](https://data.gov.au/data/dataset/808fd05a-f451-4170-a1d3-8439105c0273)\n\n* **Derived From** [Geological Provinces - Full Extent](https://data.gov.au/data/dataset/0a064b3b-2a1e-4672-80d9-03333be67aad)\n\n* **Derived From** [Bioregional Assessment areas v03](https://data.gov.au/data/dataset/96dbf469-5463-4f4d-8fad-4214c97e5aac)\n\n* **Derived From** [BOM, Australian Average Rainfall Data from 1961 to 1990](https://data.gov.au/data/dataset/fd91f2d4-2cc8-4d5d-9f67-8fe8af1e2676)\n\n* **Derived From** [GIS analysis of HYDMEAS - Hydstra Groundwater Measurement Update: NSW Office of Water - Nov2013](https://data.gov.au/data/dataset/d414c703-aabd-43af-81e0-30aab4d9dfb1)\n\n* **Derived From** [NSW Catchment Management Authority Boundaries 20130917](https://data.gov.au/data/dataset/0dc4272e-c081-4ed4-9645-aa336ffacd85)\n\n* **Derived From** [Australian 0.05º gridded chloride deposition v2](https://data.gov.au/data/dataset/c1649bd7-227f-41ff-9964-b55479bef640)\n\n* **Derived From** [Hydstra Groundwater Measurement Update - NSW Office of Water, Nov2013](https://data.gov.au/data/dataset/08b34518-c947-4665-82e5-c144b5b24e01)\n\n* **Derived From** [Namoi dryland diffuse groundwater recharge](https://data.gov.au/data/dataset/858327d1-a558-4edc-9315-8779924f3632)\n\n* **Derived From** [Bioregional Assessment areas v01](https://data.gov.au/data/dataset/0fd0d820-30b3-451f-a126-a50040426999)\n\n* **Derived From** [Bioregional Assessment areas v02](https://data.gov.au/data/dataset/e414b1b3-c42e-42e9-9cc4-2093054aa35f)\n\n* **Derived From** [Namoi groundwater model](https://data.gov.au/data/dataset/5ad7be74-e091-4f95-a2bf-6033ae7b933a)\n\n* **Derived From** [Namoi bore locations, depth to water for June 2012](https://data.gov.au/data/dataset/1deed9de-9299-4092-a68d-57852acf1a0c)\n\n* **Derived From** [NSW Office of Water Groundwater Entitlements Spatial Locations](https://data.gov.au/data/dataset/0f4da0d0-9e82-42f3-84da-1148e0b730c1)\n\n* **Derived From** [Victoria - Seamless Geology 2014](https://data.gov.au/data/dataset/2872d02e-66cb-42b6-9e5a-63abc8ad871b)\n\n* **Derived From** [Namoi NSW Office of Water groundwater licence BA purpose](https://data.gov.au/data/dataset/15a48026-27f1-4239-b278-32a23e34db5f)\n\n

**摘要**

本数据集由生物区域评估计划（Bioregional Assessment Programme）基于多份源数据集衍生而来，源数据集的信息可在本元数据说明的“谱系（Lineage）”字段中查询，该衍生数据集的生成流程则见于元数据说明的“历史（History）”字段。

本数据集包含纳米伊次区域（Namoi subregion）所有地下水模型节点处的最大降深（maximum drawdown）及最大降深出现时间（time to maximum drawdown）预测结果，该结果基于地下水位、河流通量与矿井产水率的观测数据进行约束。数据集同时涵盖敏感性分析所需的脚本及分析结果，包含了从地下水模型（纳米伊地下水模型数据集，Namoi groundwater model dataset）输出中生成上述结果的全部脚本，以及存储所有结果的电子表格。本数据集的研究方法与结果详见Janardhanan等人2017年的研究成果。

参考文献
Janardhanan S, Crosbie R, Pickett T, Cui T, Peeters L, Slatter E, Northey J, Merrin LE, Davies P, Miotlinski K, Schmid W and Herr A (2017) 纳米伊次区域地下水数值模拟. 北内陆集水区生物区域评估纳米伊次区域成果2.6.2. 澳大利亚环境与能源部、澳大利亚气象局（Bureau of Meteorology）、联邦科学与工业研究组织（CSIRO）及澳大利亚地质科学局（Geoscience Australia）, 澳大利亚. http://data.bioregionalassessments.gov.au/product/NIC/NAM/2.6.2.

**数据集历史**

支撑本数据集的工作流程已记录于文件'NAM_MF_UA_workflow.png'中。

电子表格'NAM_MF_dmax_Predictions_all.csv'源自'纳米伊地下水模型'数据集，其中包含每个地下水模型节点的名称、坐标、模型中钻孔编号（Bore_ID）、层号、目标函数名称，以及该节点最大降深（dmax）的试验设计（design of experiment, DoE）运行结果的最小值、最大值、中位数、5%分位数与95%分位数。每次试验设计运行中各节点的单独结果存储于电子表格'NAM_MF_dmax_DoE_Predictions_all.csv'。对应最大降深出现时间（tmax）的等效文件为'NAM_MF_tmax_Predictions_all.csv'与'NAM_MF_tmax_DoE_Predictions_all.csv'。

在Python脚本'NAM_MF_datawranling.csv'中，上述文件与存储地下水位、矿井排水速率及河流通量观测值的文件'NAM_MF_Observations_all.csv'，以及存储各预测结果与各矿井、地下水位观测点及河流之间距离的文件'NAM_MF_dist_hobs.csv'、'NAM_MF_dist_rivers.csv'、'NAM_MF_dist_mines.csv'进行合并。该脚本仅筛选出最大降深95%分位数小于1厘米的预测结果用于后续分析，筛选后的预测子集存储于'NAM_MF_dmax_Predictions.csv'、'NAM_MF_tmax_Predictions.csv'、'NAM_MF_dmax_DoE_Predictions.csv'与'NAM_MF_tmax_DoE_Predictions.csv'中。输出电子表格'NAM_MF_Observations.csv'则包含观测数据及与筛选后预测结果的距离信息。

由于观测值的模拟等效结果属于预测数据集的一部分，因此在Python脚本'NAM_MF_OFs.py'中，将这些文件合并以生成每次试验运行及各预测结果对应的目标函数值。目标函数值为残差的加权和，根据各预测结果的距离分别存储于'NAM_MF_DoE_hres.csv'、'NAM_MF_DoE_mres.csv'与'NAM_MF_DoE_rres.csv'，最终结果存储于'NAM_MF_DoE_OFh.csv'、'NAM_MF_DoE_OFm.csv'与'NAM_MF_DoE_OFr.csv'。各目标函数与各预测结果的阈值存储于'NAM_MF_OF_thresholds.csv'。Python脚本'NAM_MF_OF_wrangling.py'会对上述信息进行进一步后处理，生成接受率并存储于电子表格'NAM_MF_dmax_Predictions_ARs.csv'中。

Python脚本'NAM_MF_CreatePosterior.py'会筛选出满足接受准则的试验设计运行结果，这些结果构成后验预测分布，存储于'NAM_MF_dmax_Posterior.csv'与'NAM_MF_tmax_Posterior.csv'中，后续进一步汇总为'NAM_MF_Predictions_summary.csv'。

敏感性分析通过脚本'NAM_MF_SI.py'完成，该脚本使用试验设计的运行结果与存储于'NAM_MF_DoE_Parameters.csv'中的参数值，以及存储于'NAM_MF_Parameters.csv'中的参数说明（名称、范围、变换方式）。针对最大降深、最大降深出现时间、河流、水头及矿井水流观测值的敏感性指数结果分别存储于'NAM_MF_SI_dmax.csv'、'NAM_MF_SI_tmax.csv'、'NAM_MF_SI_river.csv'、'NAM_MF_SI_mine.csv'与'NAM_MF_SI_head.csv'。以'xxxx'结尾的中间文件为每100个预测结果分组后的结果，脚本'NAM_MF_SI_collate.py'与'NAM_MF_SI_collate.slurm'用于合并这些中间文件。

**数据集引用**

生物区域评估计划（2017） 纳米伊地下水不确定性分析. 生物区域评估衍生数据集. 2018年12月11日查阅，http://data.bioregionalassessments.gov.au/dataset/36bd27e9-58d2-4bf2-8e4a-54b22ac98cfb.

**数据集祖先**

* **源自** [新南威尔士州水务局（NSW Office of Water）关联空间位置的地下水许可证提取数据（NIC v2，2014年2月28日）](https://data.gov.au/data/dataset/505e8a67-90b4-408a-aa0d-038a9fb14f42)
* **源自** [纳米伊水力传导率测量数据](https://data.gov.au/data/dataset/5f88517d-8154-411d-907f-4e2c2d12a912)
* **源自** [2012年纳米伊NGIS钻孔分析数据](https://data.gov.au/data/dataset/1c828f1b-3e9c-46ae-983a-6ea2218dba48)
* **源自** [纳米伊地下水模型冲积层范围数据](https://data.gov.au/data/dataset/e7e79c42-89aa-4278-b14c-d7c67ebd4d89)
* **源自** [澳大利亚地表地质图（1:1 000 000比例尺，2012版）](https://data.gov.au/data/dataset/8284767e-b5b1-4d8b-b8e6-b334fa972611)
* **源自** [纳米伊Leapfrog地质模型](https://data.gov.au/data/dataset/26747362-8e20-49df-ab2c-918fba839aa4)
* **源自** [历史采矿足迹（DTIRIS NAM 20150914）](https://data.gov.au/data/dataset/fb0a43f1-9022-4844-b5b5-b91f72a1d9f9)
* **源自** [吉普斯兰项目边界数据](https://data.gov.au/data/dataset/27413de5-d13a-4231-ac79-fc77f4cbb5f7)
* **源自** [生物区域评估区域v04](https://data.gov.au/data/dataset/e8a2d577-c5c5-4e2c-b0fa-53e4b2d4a034)
* **源自** [2010年自然资源管理（NRM）区域数据](https://data.gov.au/data/dataset/1d54e38f-4051-4f0c-a350-c7dbd8eba65b)
* **源自** [土壤与景观网格国家土壤属性图——3分辨率黏土数据（版本1）](https://data.gov.au/data/dataset/f8640540-4bb7-42ee-995a-219881e67705)
* **源自** [GEODATA TOPO 250K系列3 文件地理数据库格式（.gdb）](https://data.gov.au/data/dataset/96ebf889-f726-4967-9964-714fb57d679b)
* **源自** [生物区域评估计划——澳大利亚集水区尺度土地利用数据（2014）](https://data.gov.au/data/dataset/6f72f73c-8a61-4ae9-b8b5-3f67ec918826)
* **源自** [GEODATA TOPO 250K系列3](https://data.gov.au/data/dataset/a0650f18-518a-4b99-a553-44f82f28bb5f)
* **源自** [新南威尔士州水务局地下水许可证提取数据（NIC-2013年10月）](https://data.gov.au/data/dataset/808fd05a-f451-4170-a1d3-8439105c0273)
* **源自** [地质省——全范围数据](https://data.gov.au/data/dataset/0a064b3b-2a1e-4672-80d9-03333be67aad)
* **源自** [生物区域评估区域v03](https://data.gov.au/data/dataset/96dbf469-5463-4f4d-8fad-4214c97e5aac)
* **源自** [澳大利亚气象局1961-1990年平均降雨量数据](https://data.gov.au/data/dataset/fd91f2d4-2cc8-4d5d-9f67-8fe8af1e2676)
* **源自** [HYDMEAS水文测量数据分析——新南威尔士州水务局地下水测量更新：2013年11月](https://data.gov.au/data/dataset/d414c703-aabd-43af-81e0-30aab4d9dfb1)
* **源自** [新南威尔士州集水区管理局边界数据（20130917）](https://data.gov.au/data/dataset/0dc4272e-c081-4ed4-9645-aa336ffacd85)
* **源自** [澳大利亚0.05°网格氯化物沉降数据v2](https://data.gov.au/data/dataset/c1649bd7-227f-41ff-9964-b55479bef640)
* **源自** [Hydstra地下水测量更新——新南威尔士州水务局，2013年11月](https://data.gov.au/data/dataset/08b34518-c947-4665-82e5-c144b5b24e01)
* **源自** [纳米伊旱地弥散地下水补给数据](https://data.gov.au/data/dataset/858327d1-a558-4edc-9315-8779924f3632)
* **源自** [生物区域评估区域v01](https://data.gov.au/data/dataset/0fd0d820-30b3-451f-a126-a50040426999)
* **源自** [生物区域评估区域v02](https://data.gov.au/data/dataset/e414b1b3-c42e-42e9-9cc4-2093054aa35f)
* **源自** [纳米伊地下水模型](https://data.gov.au/data/dataset/5ad7be74-e091-4f95-a2bf-6033ae7b933a)
* **源自** [2012年6月纳米伊钻孔位置及水深数据](https://data.gov.au/data/dataset/1deed9de-9299-4092-a68d-57852acf1a0c)
* **源自** [新南威尔士州水务局地下水权利空间位置数据](https://data.gov.au/data/dataset/0f4da0d0-9e82-42f3-84da-1148e0b730c1)
* **源自** [维多利亚州——无缝地质数据2014](https://data.gov.au/data/dataset/2872d02e-66cb-42b6-9e5a-63abc8ad871b)
* **源自** [纳米伊新南威尔士州水务局地下水许可证BA用途数据](https://data.gov.au/data/dataset/15a48026-27f1-4239-b278-32a23e34db5f)