1 Department of Physical Geography, Utrecht University, The Netherlands
2 Department of Subsurface and Groundwater Systems, Deltares Research Institute, Utrecht, The Netherlands
3 Division of Water Resources Planning and Investigation for the South of Vietnam (DWRPIS), Ho Chi Minh city, Vietnam
Interbed
Subsidence Confined aquifer
Unconfined aquifer
water table
Creep
Ripening Oxidation
Consolidation
Creep
Isostasy Consolidation
Creep
Consolidation
T otal
Subsidence
S hallow
D eep
DRIVERS OF SUBSIDENCE
Colourcode:
Loading Artificial lowering
of groundwater table
Fluid extraction Tectonics &
Isostasy
Process Natural driver Antropogenic driver
Natural loading Infrastruc tur e Buildings Lar ge construc tions Drainage of sur fac e wat er Gr oundwat er Ear th crust dynamics
H ydr ocarbons
Aquitard
Aquitard
Bedrock Consolidation
Seismicity Autocompaction
Boundary shallow/deep
w w w
waaaatttteeerrrr tttttaaaabbbblllleeee
C C C Cr r r r re e e ee e e e e
R R R R
Ri i i ip p p pe e e e en n n ni i i i in n n ng g g g g O
O O O
Ox x x x xi i i id d d d da a a at t t t ti i io o o on n n n
C C C C
Co o o o on n n n ns s s so o o o ol l l l li i i id d d d da a a at t t t ti i i io o o o on n n n C
C C C
Cr r r re e e e ee e e e ep p p p C
C C
Co o o on n n ns s s so o o ol l l li id d d da a a at t ti i i io o o o on n n n A
A A A
Au u u u ut t t t to o o oc c c co o o o om m m m mp p p p pa a a ac c c c ct t t ti i i io o o on n n n n
U U U U
Un n n n nc c c co o o o on n n nfi fi fi fi fin n n ne e e ed d d d d C
C C C
Cr r r r re e e e ee e e ep p p p
A A A A A A A A A A
Aq q q q q q q q q qu u u u u u u u u u u u ui i i i i i i i i it t t t t t t t t ta a a a a a a a a ar r r r r r r r r rd d d d d d d d d d
B B B B B B B B B
Booooooooooooouuuuuuuuunnnnnnnnnndddddddddddaaaaaaaaaarrrrrrrrryyyyyyyyyy yyy
A A A A A A
Aq q q q q q qu u u u u u u ui i i i i i i i it t t t t ta a a a a ar r r r rd d d d d d d C
C C Cr r r re e e ee e e e d
d d
d a a a aq q q q qu u u u ui i i if f f f fe e e e e
Subsidence Colourcode: Process Natural driver Antropogenic driver
y y y y C
C C C
Co o o o on n n n nfi fi fi fi fin n n n ne e e e ed d d d d a a a a aq q q q qu u u u ui i i i if f f f fe e e e er r r r C C C Cr r r r re e e e ee e e e ep p p p
s s s
shhhhhaaaallllllllloooowwww/////dddddeeeeeeeepppp
C C C C
Co o o o on n n ns s s so o o ol l l li i i id d d da a a a at t t ti i i i i
A A A A A A A
Aq q q q q q q q q q qu u u u u u u u ui i i i i i i i i it t t t t t t t t ta a a a a a a a a ar r r r r r r rd d d d d d d d d p
p p p p p p p e e e er r r r C
C C C
Co o o o on n n n ns s s s so o o ol l l li i i id d d d da a a a at t t ti i i io o o o on n n n n
I I I I I I I I I
In n n n n n n n nt t t t t t t t t te e e e e e e er r r r r r r rb b b b b b b b b be e e e e e e e ed d d d d d d d d d
i i io o o o on n n n n
Interbed
Subsidence Confined aquifer
Unconfined aquifer
water table
S S
Se e e e eis s s sm m m m mi i i ic c c ci it t t t ty y y y C
C C C
Co o o o on n n ns s s s so o o o oli i i id d d d da a a at t t ti i i io o o o on n n n n S
S S S
Se e e e ei i i is s s s sm m m m mi i i i ic c c ci i it t t t ty y y y y Creep
Ripening Oxidation
I I
Is s s so o o os s s s st t t ta a a a as s s sy y y y Consolidation
Creep
Isostasy Consolidation
Creep
Consolidation
T otal
Subsidence
S hallow
D eep
PROCESSES OF SUBSIDENCE
Colourcode:
Loading Artificial lowering
of groundwater table
Fluid extraction Tectonics &
Isostasy
Process Natural driver Antropogenic driver
Natural loading Infrastruc tur e Buildings Lar ge construc tions Drainage of sur fac e wat er Gr oundwat er Ear th crust dynamics
H ydr ocarbons
Aquitard
Aquitard
Consolidation Seismicity Autocompaction
Boundary shallow/deep
Bedrock
1: Standard penetration test 2: Logging while drilling 3: Vertical electronical sounding
Available data Additions to existing models Outcomes
Lithological borehole descriptions
Geological and geohydrological
cross-sections
Hydraulic head and extraction well data
Physical and chemical sediment properties
Lithostratigraphical analysis Conceptual models of delta
evolution
Palaeogeographical analysis (Depositional environments
& sediment preservation)
3D lithological subsurface model
3D lithological interpolation
Better understanding aquifer - aquitard architecture / properties
Improved distribution of sediment properties
Improved geo-hydrological model
Gain
2. Evaluate future groundwater management scenario’s
Identify subsidence drivers and develop mitigation strategies
Supporting decision-making towards sustainable groundwater management
1. Unravel the subsidence balance Hydrological model with
subsidence module:
Results
C
Ca an n T Th ho o C
Ci it ty y
Subsidence in the Mekong Delta, Vietnam:
Impact of groundwater extraction
Introduction
Land subsidence rates of ~1-4 cm yr
-1are measured in the low-lying Vietnamese Mekong Delta (Fig. 1).
These relatively high subsidence rates are attributed to groundwater extraction. On daily basis over two million m
3of groundwater is extracted from the upper 500 m of the multi-aquifer subsurface. As a result, hydraulic heads in aquifers are dropping, on average 0.3-0.7 m yr
-1.
P.S.J. Minderhoud 1,2 , G. Erkens 2,1 , V.H. Pham 1,2,3 , B.T. Vuong 3 , E. Stouthamer 1
Land subsidence increases flood risk, and, on the longer term, threatens the delta with drowning.
To evaluate the impact of future land subsidence, we need to go from measurements to predictions.
Here we present our approach to assess the subsidence potential of the multi-aquifer subsurface of the Mekong delta due to groundwater extraction under different groundwater management scenarios.
Figure 1. Satellite based (InSAR) subsidence rates measured between 2006-2010 for the Mekong Delta. Data © JAXA, METI 2011 (Erban et al., 2014).
Figure 2. Schematization of the main subsidence drivers and processes within the upper (phreatic) aquifer and deeper (confined) aquifer(s). Both natural and anthropogenic drivers are distinguished. The subsidence balance equation is given on the left side, being the total sum of all shallow and deep subsidence rates.
Figure 3. Workflow of the approach to develop the 3D lithological subsurface model and construct an improved geo-hydrological model to enable subsidence modeling.
From monitoring to predicting
Total measured subsidence at the earth surface is the sum of subsidence resulting from all natural and human-induced drivers, the subsidence balance (Fig. 2). We distinguish between shallow and deep drivers and processes of subsidence.
To determine the contribution of groundwater extraction, and to go from measuring to predicting subsidence, the cumulative signal needs to be unraveled. We develop a 3D geohydrogical model to model hydrology and calculate subsidence to evaluate the impact of groundwater extraction.
References
Erban, L. E., Gorelick, S. M., & Zebker, H. A. (2014). Groundwater extraction, land subsidence, and sea-level rise in the Mekong Delta, Vietnam. Environmental Research Letters, 9(8), 1–6.
* The results depicted are preliminary model outputs before model calibration.
Acknowledgements
This poster is part of a PhD research carried out by P.S.J. Minderhoud at the Dept. of Physical Geography, Utrecht University, The Netherlands. The PhD project is funded by NWO-WOTRO (W 07.69.105), Deltares and TNO-Geological Survey of the Netherlands.
The Division of Water Resources Planning and Investigation for the South of Vietnam (DWRPIS), Ho Chi Minh city, Vietnam is thanked for providing subsurface and hydrological data for this research.
Figure 4. Measured hydraulic head time series from monitoring wells near Can Tho city, central Mekong delta. Filter depths between brackets.
Figure 5. Visualisation of the 3D geo-hydrological model in iMOD (MODFLOW shell by Deltares) showing the DEM and subsurface architecture.
Figure 6. Spatial variability of hydraulic head decline of the Middle Pleistocene aquifer after a 25-year model run (1990-2015)*.
Figure 7. Total calculated subsidence for all layers (1990-2015) modeled using the coupled SUB-Cr module in iMOD (NEN-Bjerrum method)*.
-6.00 -5.00 -4.00 -3.00 -2.00 -1.00 0.00 1.00 2.00
1995 1997 1999 2001 2003 2005 2007 2009 2011 2013
Holocene (+2 to -69 m)
Late Pleistocene (-69 to -130 m) Middle Pleistocene (-130 to -204 m) Upper Pliocene (-277 to -326)
Middle Pliocene (-326 to -398) -6.00
-5.00 -4.00 -3.00 -2.00 -1.00 0.00 1.00 2.00
1995 1997 1999 2001 2003 2005 2007 2009 2011 2013
Holocene (+2 to -69 m) Late Pleistocene (-69 to -130 m) Middle Pleistocene (-130 to -204 m) Upper Pliocene (-277 to -326 m) Middle Pliocene (-326 to -398 m) Hydraulic heads in aquifers:
Faculty of Geosciences Dep. Physical Geography
rlands.