Regional flow in a complex coastal aquifer system: Combining voxel geological modelling with regularized calibration

Rena Meyer; Peter Engesgaard; Anne Sophie Høyer; Flemming Jørgensen; Giulio Vignoli; Torben O. Sonnenborg

doi:10.1016/j.jhydrol.2018.05.020

Regional flow in a complex coastal aquifer system: Combining voxel geological modelling with regularized calibration

Rena Meyer^*, Peter Engesgaard, Anne Sophie Høyer, Flemming Jørgensen, Giulio Vignoli, Torben O. Sonnenborg

^*Corresponding author for this work

Geology

7 Citations (Scopus)

Abstract

Low-lying coastal regions are often highly populated, constitute sensitive habitats and are at the same time exposed to challenging hydrological environments due to surface flooding from storm events and saltwater intrusion, which both may affect drinking water supply from shallow and deeper aquifers. Near the Wadden Sea at the border of Southern Denmark and Northern Germany, the hydraulic system (connecting groundwater, river water, and the sea) was altered over centuries (until the 19th century) by e.g. the construction of dikes and drains to prevent flooding and allow agricultural use. Today, massive saltwater intrusions extend up to 20 km inland. In order to understand the regional flow, a methodological approach was developed that combined: (1) a highly-resolved voxel geological model, (2) a ∼1 million node groundwater model with 46 hydrofacies coupled to rivers, drains and the sea, (3) Tikhonov regularization calibration using hydraulic heads and average stream discharges as targets and (4) parameter uncertainty analysis. It is relatively new to use voxel models for constructing geological models that often have been simplified to stacked, pseudo-3D layer geology. The study is therefore one of the first to combine a voxel geological model with state-of-the-art flow calibration techniques. The results show that voxel geological modelling, where lithofacies information are transferred to each volumetric element, is a useful method to preserve 3D geological heterogeneity on a local scale, which is important when distinct geological features such as buried valleys are abundant. Furthermore, it is demonstrated that simpler geological models and simpler calibration methods do not perform as well. The proposed approach is applicable to many other systems, because it combines advanced and flexible geological modelling and flow calibration techniques. This has led to new insights in the regional flow patterns and especially about water cycling in the marsh area near the coast based on the ability to define six predictive scenarios from the linear analysis of parameter uncertainty. The results show that the coastal system near the Danish-German border is mainly controlled by flow in the two aquifers separated by a thick clay layer, and several deep high-permeable buried valleys that connect the sea with the interior and the two aquifers. The drained marsh area acts like a huge regional sink limiting submarine groundwater discharge. With respect to water balance, the greatest sensitivity to parameter uncertainty was observed in the drained marsh area, where some scenarios showed increased flow of sea water into the interior and increased drainage. We speculate that the massive salt water intrusion may be caused by a combination of the preferential pathways provided by the buried valleys, the marsh drainage and relatively high hydraulic conductivities in the two main aquifers as described by one of the scenarios. This is currently under investigation by using a salt water transport model.

Original language	English
Journal	Journal of Hydrology
Volume	562
Pages (from-to)	544-563
Number of pages	20
ISSN	0022-1694
DOIs	https://doi.org/10.1016/j.jhydrol.2018.05.020
Publication status	Published - 1 Jul 2018

Keywords

Highly heterogeneous coastal aquifer system
Linear uncertainty analysis
Regional groundwater flow model (MODFLOW)
Regularized inversion (PEST)
Voxel geology

Access to Document

10.1016/j.jhydrol.2018.05.020

Cite this

@article{631183c2bfb4456587e51f9354c6b562,

title = "Regional flow in a complex coastal aquifer system: Combining voxel geological modelling with regularized calibration",

abstract = "Low-lying coastal regions are often highly populated, constitute sensitive habitats and are at the same time exposed to challenging hydrological environments due to surface flooding from storm events and saltwater intrusion, which both may affect drinking water supply from shallow and deeper aquifers. Near the Wadden Sea at the border of Southern Denmark and Northern Germany, the hydraulic system (connecting groundwater, river water, and the sea) was altered over centuries (until the 19th century) by e.g. the construction of dikes and drains to prevent flooding and allow agricultural use. Today, massive saltwater intrusions extend up to 20 km inland. In order to understand the regional flow, a methodological approach was developed that combined: (1) a highly-resolved voxel geological model, (2) a ∼1 million node groundwater model with 46 hydrofacies coupled to rivers, drains and the sea, (3) Tikhonov regularization calibration using hydraulic heads and average stream discharges as targets and (4) parameter uncertainty analysis. It is relatively new to use voxel models for constructing geological models that often have been simplified to stacked, pseudo-3D layer geology. The study is therefore one of the first to combine a voxel geological model with state-of-the-art flow calibration techniques. The results show that voxel geological modelling, where lithofacies information are transferred to each volumetric element, is a useful method to preserve 3D geological heterogeneity on a local scale, which is important when distinct geological features such as buried valleys are abundant. Furthermore, it is demonstrated that simpler geological models and simpler calibration methods do not perform as well. The proposed approach is applicable to many other systems, because it combines advanced and flexible geological modelling and flow calibration techniques. This has led to new insights in the regional flow patterns and especially about water cycling in the marsh area near the coast based on the ability to define six predictive scenarios from the linear analysis of parameter uncertainty. The results show that the coastal system near the Danish-German border is mainly controlled by flow in the two aquifers separated by a thick clay layer, and several deep high-permeable buried valleys that connect the sea with the interior and the two aquifers. The drained marsh area acts like a huge regional sink limiting submarine groundwater discharge. With respect to water balance, the greatest sensitivity to parameter uncertainty was observed in the drained marsh area, where some scenarios showed increased flow of sea water into the interior and increased drainage. We speculate that the massive salt water intrusion may be caused by a combination of the preferential pathways provided by the buried valleys, the marsh drainage and relatively high hydraulic conductivities in the two main aquifers as described by one of the scenarios. This is currently under investigation by using a salt water transport model.",

keywords = "Highly heterogeneous coastal aquifer system, Linear uncertainty analysis, Regional groundwater flow model (MODFLOW), Regularized inversion (PEST), Voxel geology",

author = "Rena Meyer and Peter Engesgaard and H{\o}yer, {Anne Sophie} and Flemming J{\o}rgensen and Giulio Vignoli and Sonnenborg, {Torben O.}",

year = "2018",

month = jul,

day = "1",

doi = "10.1016/j.jhydrol.2018.05.020",

language = "English",

volume = "562",

pages = "544--563",

journal = "Journal of Hydrology",

issn = "0022-1694",

publisher = "Elsevier",

}

TY - JOUR

T1 - Regional flow in a complex coastal aquifer system

T2 - Combining voxel geological modelling with regularized calibration

AU - Meyer, Rena

AU - Engesgaard, Peter

AU - Høyer, Anne Sophie

AU - Jørgensen, Flemming

AU - Vignoli, Giulio

AU - Sonnenborg, Torben O.

PY - 2018/7/1

Y1 - 2018/7/1

N2 - Low-lying coastal regions are often highly populated, constitute sensitive habitats and are at the same time exposed to challenging hydrological environments due to surface flooding from storm events and saltwater intrusion, which both may affect drinking water supply from shallow and deeper aquifers. Near the Wadden Sea at the border of Southern Denmark and Northern Germany, the hydraulic system (connecting groundwater, river water, and the sea) was altered over centuries (until the 19th century) by e.g. the construction of dikes and drains to prevent flooding and allow agricultural use. Today, massive saltwater intrusions extend up to 20 km inland. In order to understand the regional flow, a methodological approach was developed that combined: (1) a highly-resolved voxel geological model, (2) a ∼1 million node groundwater model with 46 hydrofacies coupled to rivers, drains and the sea, (3) Tikhonov regularization calibration using hydraulic heads and average stream discharges as targets and (4) parameter uncertainty analysis. It is relatively new to use voxel models for constructing geological models that often have been simplified to stacked, pseudo-3D layer geology. The study is therefore one of the first to combine a voxel geological model with state-of-the-art flow calibration techniques. The results show that voxel geological modelling, where lithofacies information are transferred to each volumetric element, is a useful method to preserve 3D geological heterogeneity on a local scale, which is important when distinct geological features such as buried valleys are abundant. Furthermore, it is demonstrated that simpler geological models and simpler calibration methods do not perform as well. The proposed approach is applicable to many other systems, because it combines advanced and flexible geological modelling and flow calibration techniques. This has led to new insights in the regional flow patterns and especially about water cycling in the marsh area near the coast based on the ability to define six predictive scenarios from the linear analysis of parameter uncertainty. The results show that the coastal system near the Danish-German border is mainly controlled by flow in the two aquifers separated by a thick clay layer, and several deep high-permeable buried valleys that connect the sea with the interior and the two aquifers. The drained marsh area acts like a huge regional sink limiting submarine groundwater discharge. With respect to water balance, the greatest sensitivity to parameter uncertainty was observed in the drained marsh area, where some scenarios showed increased flow of sea water into the interior and increased drainage. We speculate that the massive salt water intrusion may be caused by a combination of the preferential pathways provided by the buried valleys, the marsh drainage and relatively high hydraulic conductivities in the two main aquifers as described by one of the scenarios. This is currently under investigation by using a salt water transport model.

AB - Low-lying coastal regions are often highly populated, constitute sensitive habitats and are at the same time exposed to challenging hydrological environments due to surface flooding from storm events and saltwater intrusion, which both may affect drinking water supply from shallow and deeper aquifers. Near the Wadden Sea at the border of Southern Denmark and Northern Germany, the hydraulic system (connecting groundwater, river water, and the sea) was altered over centuries (until the 19th century) by e.g. the construction of dikes and drains to prevent flooding and allow agricultural use. Today, massive saltwater intrusions extend up to 20 km inland. In order to understand the regional flow, a methodological approach was developed that combined: (1) a highly-resolved voxel geological model, (2) a ∼1 million node groundwater model with 46 hydrofacies coupled to rivers, drains and the sea, (3) Tikhonov regularization calibration using hydraulic heads and average stream discharges as targets and (4) parameter uncertainty analysis. It is relatively new to use voxel models for constructing geological models that often have been simplified to stacked, pseudo-3D layer geology. The study is therefore one of the first to combine a voxel geological model with state-of-the-art flow calibration techniques. The results show that voxel geological modelling, where lithofacies information are transferred to each volumetric element, is a useful method to preserve 3D geological heterogeneity on a local scale, which is important when distinct geological features such as buried valleys are abundant. Furthermore, it is demonstrated that simpler geological models and simpler calibration methods do not perform as well. The proposed approach is applicable to many other systems, because it combines advanced and flexible geological modelling and flow calibration techniques. This has led to new insights in the regional flow patterns and especially about water cycling in the marsh area near the coast based on the ability to define six predictive scenarios from the linear analysis of parameter uncertainty. The results show that the coastal system near the Danish-German border is mainly controlled by flow in the two aquifers separated by a thick clay layer, and several deep high-permeable buried valleys that connect the sea with the interior and the two aquifers. The drained marsh area acts like a huge regional sink limiting submarine groundwater discharge. With respect to water balance, the greatest sensitivity to parameter uncertainty was observed in the drained marsh area, where some scenarios showed increased flow of sea water into the interior and increased drainage. We speculate that the massive salt water intrusion may be caused by a combination of the preferential pathways provided by the buried valleys, the marsh drainage and relatively high hydraulic conductivities in the two main aquifers as described by one of the scenarios. This is currently under investigation by using a salt water transport model.

KW - Highly heterogeneous coastal aquifer system

KW - Linear uncertainty analysis

KW - Regional groundwater flow model (MODFLOW)

KW - Regularized inversion (PEST)

KW - Voxel geology

UR - http://www.scopus.com/inward/record.url?scp=85047074383&partnerID=8YFLogxK

U2 - 10.1016/j.jhydrol.2018.05.020

DO - 10.1016/j.jhydrol.2018.05.020

M3 - Journal article

AN - SCOPUS:85047074383

SN - 0022-1694

VL - 562

SP - 544

EP - 563

JO - Journal of Hydrology

JF - Journal of Hydrology

ER -

Regional flow in a complex coastal aquifer system: Combining voxel geological modelling with regularized calibration

Abstract

Keywords

Access to Document

Other files and links

Fingerprint

Cite this