Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Constraining the response of continental-scale groundwater flow to climate change #23

Open
2 tasks done
dansand opened this issue May 15, 2024 · 12 comments
Open
2 tasks done
Labels
model approved Model approved model created Model created model published Model published new model New model review requested Review requested

Comments

@dansand
Copy link
Collaborator

dansand commented May 15, 2024

-> submitter ORCID (or name)

0000-0002-2207-6837

-> slug

mather-2022-groundwater

-> license

CC-BY-4.0

-> alternative license URL

No response

-> model category

model published in study, inverse model

-> model status

completed

-> associated publication DOI

http://dx.doi.org/10.1038/s41598-022-08384-w

-> model creators

0000-0003-3566-1557
0000-0002-3334-5764
0000-0002-6034-1881
0000-0002-7182-1864
0000-0002-6557-0237
0000-0003-3685-174X

-> title

No response

-> description

This model was developed in order to study groundwater flow on a continental scale, focusing on the Sydney–Gunnedah–Bowen Basin in Australia. Using data such as hydraulic head measurements and borehole temperatures, it predicts how water moves through deep aquifers to the surface. Coastal aquifers show fast water flow, while inland aquifers have much slower flow. The study shows that increased water extraction from inland areas could permanently change water flow patterns. This open-source model can be used for other regions and aims to support sustainable groundwater management policies

-> abstract

Numerical models of groundwater flow play a critical role for water management scenarios under climate extremes. Large-scale models play a key role in determining long range flow pathways from continental interiors to the oceans, yet struggle to simulate the local flow patterns offered by small-scale models. We have developed a highly scalable numerical framework to model continental groundwater flow which capture the intricate flow pathways between deep aquifers and the near-surface. The coupled thermal-hydraulic basin structure is inferred from hydraulic head measurements, recharge estimates from geochemical proxies, and borehole temperature data using a Bayesian framework. We use it to model the deep groundwater flow beneath the Sydney–Gunnedah–Bowen Basin, part of Australia’s largest aquifer system. Coastal aquifers have flow rates of up to 0.3 m/day, and a corresponding groundwater residence time of just 2,000 years. In contrast, our model predicts slow flow rates of 0.005 m/day for inland aquifers, resulting in a groundwater residence time of 400,000 years. Perturbing the model to account for a drop in borehole water levels since 2000, we find that lengthened inland flow pathways depart significantly from pre-2000 streamlines as groundwater is drawn further from recharge zones in a drying climate. Our results illustrate that progressively increasing water extraction from inland aquifers may permanently alter long-range flow pathways. Our open-source modelling approach can be extended to any basin and may help inform policies on the sustainable management of groundwater.

-> scientific keywords

groundwater, thermal-hydraulic, Bayesian, water-management

-> funder

NSW Department of Industry
https://ror.org/04s1m4564

-> model embargo?

No response

-> include model code ?

  • yes

-> model code/inputs DOI

https://github.com/brmather/Sydney_Basin/tree/master

-> model code/inputs notes

In the Scripts folder, HL05 was used to run the optimisation problem and HL06 was used to take the maximum a posteriori model and run it at high resolution.

-> include model output data?

  • yes

-> data creators

No response

-> model output data DOI

No response

-> model output data notes

model_output_data contains the following file types:

.h5 - Underworld2 data files
.xdmf- Underworld2 xdmf header files
.csv - Various data in csv format
.npz - data on numpy binary format
.png - image files
.pvsm - Paraview state files
.txt - data in .txt format

-> model output data size

15 Gb

-> software framework DOI/URI

https://doi.org/10.5281/zenodo.7455999

-> software framework source repository

https://github.com/underworldcode/underworld2

-> name of primary software framework (e.g. Underworld, ASPECT, Badlands, OpenFOAM)

No response

-> software framework authors

No response

-> software & algorithm keywords

Python, C, finite element, heat equation, advection-diffusion

-> computer URI/DOI

https://ror.org/04yx6dh41

-> add landing page image and caption

No response

-> add an animation (if relevant)

No response

-> add a graphic abstract figure (if relevant)

Coupled heat-groundwater flow model of the Sydney–Gunnedah–Bowen Basin based on the MAP estimate of material properties and boundary conditions. (A) Groundwater velocity field with coal seams outlined in grey overlain with temperature gradients measured in boreholes. This visualisation of the velocity field obtained from our model was rendered in 3D using Paraview 5.9 (https://www.paraview.org/). (B) temperature field overlain with heat flux vectors. The 2D slice was generated from our models using Matplotlib 3.4 (https://matplotlib.org/).

fig1

-> add a model setup figure (if relevant)

3D stratigraphy of the Sydney–Gunnedah–Bowen Basin. The vertical spacing of layers has been exaggerated for visual clarity. The model of the basin was rendered in 3D using Underworld.
figure_2

-> add a description of your model setup

In this paper, we apply our numerical framework to the Sydney–Gunnedah–Bowen (SGB) Basin in eastern Australia. The SGB Basin covers about 1.5 million square kilometers, and we model it in high-resolution 3D, using over 10 million cells (or 6 x 6 x 0.6 km, in the x, y, z directions, respectively) to detail flow patterns down to 12 km beneath the crust. By adjusting the model to match real-world data, it provides accurate insights into water and heat movement through deep aquifers in large areas. Temperature advection due to groundwater flow is described by the advection-diffusion equation. Darcy flux is calculated from the groundwater flow equation. Groundwater recharge and discharge are driven by changes in hydraulic head, which is set to the height of the water table at the top boundary surface. The thermal boundary conditions include a constant temperature set to the top boundary, which corresponds to the annual mean surface temperature. The side walls are assigned zero flux, and the bottom temperature boundary is an unknown variable that we invert from borehole temperature data within our Bayesian optimization scheme.

Please provide any feedback on the model submission process?

No response

@dansand dansand added new model New model review requested Review requested labels May 15, 2024
@m-te-bot
Copy link

m-te-bot bot commented May 15, 2024

Review Requested

A review of this submission has been requested from @ModelAtlasofTheEarth/model_reviewers

@dansand dansand added the model approved Model approved label May 15, 2024
@m-te-bot
Copy link

m-te-bot bot commented May 15, 2024

Model repository created at https://github.com/ModelAtlasofTheEarth/mather-2022-groundwater-1

@m-te-bot m-te-bot bot added the model created Model created label May 15, 2024
@dansand dansand removed the model approved Model approved label May 15, 2024
@ModelAtlasofTheEarth ModelAtlasofTheEarth deleted a comment from m-te-bot bot May 15, 2024
@m-te-bot
Copy link

m-te-bot bot commented May 15, 2024

Model Report

Thank you for submitting.

Using Github actions, we have regenerated a report summarising information about your model

  • Please check the report below, including the Errors and Warnings section

  • You can update any information, by editing the markdown file at the top of the issue

  • these edits will trigger the report will be regenerated

  • once you are satisfied with the results, please add a review requested Review requested label

Parsed data

Section 1: Summary of your model

Model Submitter:

Dan Sandiford (0000-0002-2207-6837)

Model Creator(s):

Model name:

mather-2022-groundwater

(this will be the name of the model repository when created)

Model long name:

Constraining the response of continental-scale groundwater flow to climate change

License:

Creative Commons Attribution 4.0 International

Model Category:

  • model published in study
  • inverse model

Model Status:

  • completed

Associated Publication title:

Constraining the response of continental-scale groundwater flow to climate change

Abstract:

Numerical models of groundwater flow play a critical role for water management scenarios under climate extremes. Large-scale models play a key role in determining long range flow pathways from continental interiors to the oceans, yet struggle to simulate the local flow patterns offered by small-scale models. We have developed a highly scalable numerical framework to model continental groundwater flow which capture the intricate flow pathways between deep aquifers and the near-surface. The coupled thermal-hydraulic basin structure is inferred from hydraulic head measurements, recharge estimates from geochemical proxies, and borehole temperature data using a Bayesian framework. We use it to model the deep groundwater flow beneath the Sydney–Gunnedah–Bowen Basin, part of Australia’s largest aquifer system. Coastal aquifers have flow rates of up to 0.3 m/day, and a corresponding groundwater residence time of just 2,000 years. In contrast, our model predicts slow flow rates of 0.005 m/day for inland aquifers, resulting in a groundwater residence time of 400,000 years. Perturbing the model to account for a drop in borehole water levels since 2000, we find that lengthened inland flow pathways depart significantly from pre-2000 streamlines as groundwater is drawn further from recharge zones in a drying climate. Our results illustrate that progressively increasing water extraction from inland aquifers may permanently alter long-range flow pathways. Our open-source modelling approach can be extended to any basin and may help inform policies on the sustainable management of groundwater.

Scientific Keywords:

  • groundwater
  • thermal-hydraulic
  • Bayesian
  • water-management

Funder(s):

Section 2: your model code, output data

** No embargo on model contents requested****Include model code:**

True

Model code existing URL/DOI:

https://github.com/brmather/Sydney_Basin/tree/master

Model code notes:

In the Scripts folder, HL05 was used to run the optimisation problem and HL06 was used to take the maximum a posteriori model and run it at high resolution.

Include model output data:

True

Model output data notes:

model_output_data contains the following file types:

.h5 - Underworld2 data files
.xdmf- Underworld2 xdmf header files
.csv - Various data in csv format
.npz - data on numpy binary format
.png - image files
.pvsm - Paraview state files
.txt - data in .txt format

Section 3: software framework and compute details

Software Framework DOI/URL:

Found software: Underworld2: Python Geodynamics Modelling for Desktop, HPC and Cloud

Software Repository:

https://github.com/underworldcode/underworld2

Name of primary software framework:

Underworld2: Python Geodynamics Modelling for Desktop, HPC and Cloud

Software & algorithm keywords:

  • Python
  • C
  • finite element
  • heat equation
  • advection-diffusion

Section 4: web material (for mate.science)

Landing page image:

Filename:

Animation:

Filename:

Graphic abstract:

Filename: fig1.png
Caption: Coupled heat-groundwater flow model of the Sydney–Gunnedah–Bowen Basin based on the MAP estimate of material properties and boundary conditions. (A) Groundwater velocity field with coal seams outlined in grey overlain with temperature gradients measured in boreholes. This visualisation of the velocity field obtained from our model was rendered in 3D using Paraview 5.9 (https://www.paraview.org/). (B) temperature field overlain with heat flux vectors. The 2D slice was generated from our models using Matplotlib 3.4 (https://matplotlib.org/).

Model setup figure:

Filename: figure_2.png
Caption: 3D stratigraphy of the Sydney–Gunnedah–Bowen Basin. The vertical spacing of layers has been exaggerated for visual clarity. The model of the basin was rendered in 3D using Underworld.
Description: In this paper, we apply our numerical framework to the Sydney–Gunnedah–Bowen (SGB) Basin in eastern Australia. The SGB Basin covers about 1.5 million square kilometers, and we model it in high-resolution 3D, using over 10 million cells (or 6 x 6 x 0.6 km, in the x, y, z directions, respectively) to detail flow patterns down to 12 km beneath the crust. By adjusting the model to match real-world data, it provides accurate insights into water and heat movement through deep aquifers in large areas. Temperature advection due to groundwater flow is described by the advection-diffusion equation. Darcy flux is calculated from the groundwater flow equation. Groundwater recharge and discharge are driven by changes in hydraulic head, which is set to the height of the water table at the top boundary surface. The thermal boundary conditions include a constant temperature set to the top boundary, which corresponds to the annual mean surface temperature. The side walls are assigned zero flux, and the bottom temperature boundary is an unknown variable that we invert from borehole temperature data within our Bayesian optimization scheme.

Errors and Warnings

Associated Publication
Error fetching metadata with application/ld+json from https://api.crossref.org/works/http://dx.doi.org/10.1038/s41598-022-08384-w: 406 Client Error: Not Acceptable for url: https://api.crossref.org/works/http://dx.doi.org/10.1038/s41598-022-08384-w
Software Framework DOI/URI
doi.org metadata record succesfully extracted in json-ld format
Submitter
ORCID metadata record succesfully extracted in json-ld format

Model creators
ORCID metadata record succesfully extracted in json-ld format
ORCID metadata record succesfully extracted in json-ld format
ORCID metadata record succesfully extracted in json-ld format
ORCID metadata record succesfully extracted in json-ld format
ORCID metadata record succesfully extracted in json-ld format
ORCID metadata record succesfully extracted in json-ld format
Could not parse Embargo date. Check format is
Model creators
Error: no data creators found
Model output DOI
Warning: No DOI/URI provided.
Landing page image
Error: No image uploaded.

Animation
Warning: No animation uploaded.

Next steps

  • once the model_reviewers team has approved the model, we will create a repository for your model

@dansand dansand added the model approved Model approved label May 15, 2024
@m-te-bot
Copy link

m-te-bot bot commented May 15, 2024

Model repository created at https://github.com/ModelAtlasofTheEarth/mather-2022-groundwater

@dansand dansand added model approved Model approved and removed model created Model created model approved Model approved labels May 15, 2024
@m-te-bot
Copy link

m-te-bot bot commented May 15, 2024

Model repository created at https://github.com/ModelAtlasofTheEarth/mather-2022-groundwater

@m-te-bot m-te-bot bot added the model created Model created label May 15, 2024
@dansand dansand added model approved Model approved and removed model approved Model approved labels May 16, 2024
@m-te-bot
Copy link

m-te-bot bot commented May 16, 2024

Model repository created at https://github.com/ModelAtlasofTheEarth/mather-2022-groundwater

@dansand dansand added model approved Model approved and removed model approved Model approved labels May 16, 2024
@m-te-bot
Copy link

m-te-bot bot commented May 16, 2024

Model repository created at https://github.com/ModelAtlasofTheEarth/mather-2022-groundwater

@dansand dansand added model approved Model approved and removed model approved Model approved labels May 16, 2024
@m-te-bot
Copy link

m-te-bot bot commented May 16, 2024

Workflow Failure ❌

Unfortunately, the workflow to create the model repo has failed. View the logs here for more information:

🔗 Workflow logs

You may need to delete the following repo that was created via this workflow run since the run was not fully successful: ModelAtlasofTheEarth/mather-2022-groundwater

@dansand dansand added new model New model and removed workflow_error model approved Model approved new model New model labels May 16, 2024
@m-te-bot
Copy link

m-te-bot bot commented May 16, 2024

Model Report

Thank you for submitting.

Using Github actions, we have regenerated a report summarising information about your model

  • Please check the report below, including the Errors and Warnings section

  • You can update any information, by editing the markdown file at the top of the issue

  • these edits will trigger the report will be regenerated

  • once you are satisfied with the results, please add a review requested Review requested label

Parsed data

Section 1: Summary of your model

Model Submitter:

Dan Sandiford (0000-0002-2207-6837)

Model Creator(s):

Model name:

mather-2022-groundwater-1

(this will be the name of the model repository when created)

Model long name:

Constraining the response of continental-scale groundwater flow to climate change

License:

Creative Commons Attribution 4.0 International

Model Category:

  • model published in study
  • inverse model

Model Status:

  • completed

Associated Publication title:

Constraining the response of continental-scale groundwater flow to climate change

Abstract:

Numerical models of groundwater flow play a critical role for water management scenarios under climate extremes. Large-scale models play a key role in determining long range flow pathways from continental interiors to the oceans, yet struggle to simulate the local flow patterns offered by small-scale models. We have developed a highly scalable numerical framework to model continental groundwater flow which capture the intricate flow pathways between deep aquifers and the near-surface. The coupled thermal-hydraulic basin structure is inferred from hydraulic head measurements, recharge estimates from geochemical proxies, and borehole temperature data using a Bayesian framework. We use it to model the deep groundwater flow beneath the Sydney–Gunnedah–Bowen Basin, part of Australia’s largest aquifer system. Coastal aquifers have flow rates of up to 0.3 m/day, and a corresponding groundwater residence time of just 2,000 years. In contrast, our model predicts slow flow rates of 0.005 m/day for inland aquifers, resulting in a groundwater residence time of 400,000 years. Perturbing the model to account for a drop in borehole water levels since 2000, we find that lengthened inland flow pathways depart significantly from pre-2000 streamlines as groundwater is drawn further from recharge zones in a drying climate. Our results illustrate that progressively increasing water extraction from inland aquifers may permanently alter long-range flow pathways. Our open-source modelling approach can be extended to any basin and may help inform policies on the sustainable management of groundwater.

Scientific Keywords:

  • groundwater
  • thermal-hydraulic
  • Bayesian
  • water-management

Funder(s):

Section 2: your model code, output data

** No embargo on model contents requested****Include model code:**

True

Model code existing URL/DOI:

https://github.com/brmather/Sydney_Basin/tree/master

Model code notes:

In the Scripts folder, HL05 was used to run the optimisation problem and HL06 was used to take the maximum a posteriori model and run it at high resolution.

Include model output data:

True

Model output data notes:

model_output_data contains the following file types:

.h5 - Underworld2 data files
.xdmf- Underworld2 xdmf header files
.csv - Various data in csv format
.npz - data on numpy binary format
.png - image files
.pvsm - Paraview state files
.txt - data in .txt format

Section 3: software framework and compute details

Software Framework DOI/URL:

Found software: Underworld2: Python Geodynamics Modelling for Desktop, HPC and Cloud

Software Repository:

https://github.com/underworldcode/underworld2

Name of primary software framework:

Underworld2: Python Geodynamics Modelling for Desktop, HPC and Cloud

Software & algorithm keywords:

  • Python
  • C
  • finite element
  • heat equation
  • advection-diffusion

Section 4: web material (for mate.science)

Landing page image:

Filename:

Animation:

Filename:

Graphic abstract:

Filename: fig1.png
Caption: Coupled heat-groundwater flow model of the Sydney–Gunnedah–Bowen Basin based on the MAP estimate of material properties and boundary conditions. (A) Groundwater velocity field with coal seams outlined in grey overlain with temperature gradients measured in boreholes. This visualisation of the velocity field obtained from our model was rendered in 3D using Paraview 5.9 (https://www.paraview.org/). (B) temperature field overlain with heat flux vectors. The 2D slice was generated from our models using Matplotlib 3.4 (https://matplotlib.org/).

Model setup figure:

Filename: figure_2.png
Caption: 3D stratigraphy of the Sydney–Gunnedah–Bowen Basin. The vertical spacing of layers has been exaggerated for visual clarity. The model of the basin was rendered in 3D using Underworld.
Description: In this paper, we apply our numerical framework to the Sydney–Gunnedah–Bowen (SGB) Basin in eastern Australia. The SGB Basin covers about 1.5 million square kilometers, and we model it in high-resolution 3D, using over 10 million cells (or 6 x 6 x 0.6 km, in the x, y, z directions, respectively) to detail flow patterns down to 12 km beneath the crust. By adjusting the model to match real-world data, it provides accurate insights into water and heat movement through deep aquifers in large areas. Temperature advection due to groundwater flow is described by the advection-diffusion equation. Darcy flux is calculated from the groundwater flow equation. Groundwater recharge and discharge are driven by changes in hydraulic head, which is set to the height of the water table at the top boundary surface. The thermal boundary conditions include a constant temperature set to the top boundary, which corresponds to the annual mean surface temperature. The side walls are assigned zero flux, and the bottom temperature boundary is an unknown variable that we invert from borehole temperature data within our Bayesian optimization scheme.

Errors and Warnings

Associated Publication
Error fetching metadata with application/ld+json from https://api.crossref.org/works/http://dx.doi.org/10.1038/s41598-022-08384-w: 406 Client Error: Not Acceptable for url: https://api.crossref.org/works/http://dx.doi.org/10.1038/s41598-022-08384-w
Software Framework DOI/URI
doi.org metadata record succesfully extracted in json-ld format
Submitter
ORCID metadata record succesfully extracted in json-ld format

Model creators
ORCID metadata record succesfully extracted in json-ld format
ORCID metadata record succesfully extracted in json-ld format
ORCID metadata record succesfully extracted in json-ld format
ORCID metadata record succesfully extracted in json-ld format
ORCID metadata record succesfully extracted in json-ld format
ORCID metadata record succesfully extracted in json-ld format
Model Repository Slug
Warning: Model repo cannot be created with proposed slug mather-2022-groundwater.
Either propose a new slug or repo will be created with name mather-2022-groundwater-1.

Could not parse Embargo date. Check format is
Model creators
Error: no data creators found
Model output DOI
Warning: No DOI/URI provided.
Landing page image
Error: No image uploaded.

Animation
Warning: No animation uploaded.

Next steps

  • once the model_reviewers team has approved the model, we will create a repository for your model

@dansand dansand added the model approved Model approved label May 16, 2024
@m-te-bot
Copy link

m-te-bot bot commented May 16, 2024

Model repository created at https://github.com/ModelAtlasofTheEarth/mather-2022-groundwater

@dansand dansand added model approved Model approved and removed model approved Model approved model created Model created labels May 17, 2024
@m-te-bot
Copy link

m-te-bot bot commented May 17, 2024

Model repository created at https://github.com/ModelAtlasofTheEarth/mather-2022-groundwater

@m-te-bot m-te-bot bot added the model created Model created label May 17, 2024
@dansand dansand added model approved Model approved and removed model approved Model approved model created Model created labels May 20, 2024
@m-te-bot
Copy link

m-te-bot bot commented May 20, 2024

Model repository created at https://github.com/ModelAtlasofTheEarth/mather-2022-groundwater

@m-te-bot m-te-bot bot added model created Model created model published Model published labels May 20, 2024
@dansand dansand changed the title M@TE submission request from [@dansand] Constraining the response of continental-scale groundwater flow to climate change May 22, 2024
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
model approved Model approved model created Model created model published Model published new model New model review requested Review requested
Projects
None yet
Development

No branches or pull requests

1 participant