A Coupled One-Dimensional Numerical Simulation of the Land Subsidence Process in a Multilayer Aquifer System due to Hydraulic Head Variation in the Pumped Layer

Table 1

Material properties of the multilayer aquifer system used in the numerical simulation.

Property	Layered aquifer
	Phreatic aquifer	Aquitard I	Confined aquifer I	Aquitard II
Solid density, [kg/m^3]	2.00 × 10^3	1.85 × 10^3	2.00 × 10^3	1.90 × 10^3
Elastic modulus, [MPa]	200	100	200	150
Poisson’s ratio, [dimensionless]	0.25	0.35	0.25	0.35
Internal friction angle, [°]	28	23	28	16
Cohesion, [MPa]	0	0.3	0	0.6
Hydraulic conductivity, [10^−4 cm/s]	7200	368.4	5000	53.58
Porosity, [dimensionless]	0.25	0.4	0.25	0.4
Property	Confined aquifer II	Aquitard III	Confined aquifer III
Solid density, [kg/m^3]	2.00 × 10^3	1.90 × 10^3	2.00 × 10^3
Elastic modulus, [MPa]	250	200	300
Poisson’s ratio, [dimensionless]	0.25	0.3	0.25
Internal friction angle, φ [°]	30	19	34
Cohesion, [MPa]	0	0.8	0
Hydraulic conductivity, [10^−4 cm/s]	3200	9.91	620
Porosity, [dimensionless]	0.25	0.38	0.3