Research Article  Open Access
Yang Chen, Wen Zhao, Yong Huang, Pengjiao Jia, "Investigation of Land Subsidence Based on the Column Element Settlement Model in a SoftSoil Area", Geofluids, vol. 2019, Article ID 9530135, 16 pages, 2019. https://doi.org/10.1155/2019/9530135
Investigation of Land Subsidence Based on the Column Element Settlement Model in a SoftSoil Area
Abstract
The SuXiChang area is located in the Yangtze River Delta in southeastern Jiangsu Province, China. More than 2 m cumulative land subsidence has occurred since 1980s. A range of monitoring programs, geological investigations, and numerical modeling has been implemented in order to establish a regional rehabilitation plan. In this paper, the column element settlement model (CESM) has been established to investigate the coupling effect of soil selfweight, upper loads, and groundwater withdrawal on land subsidence. The model has been verified by comparing with the land subsidence records in Changzhou and further applied to the investigation area. The results show that extensive groundwater withdrawal is the main cause of the land subsidence. Settlement caused by soil selfweight increases with the depth and decreases with upper loads, while settlement caused by upper loads increases with the intensity of upper loads and decreases with depth. The proportion of ground settlement caused by soil selfweight, upper loads, and groundwater withdrawal is also investigated.
1. Introduction
Land subsidence is an environmental geological phenomenon which may cause significant damage and heavy economic losses. Ground water exploration and highrise building construction are the main anthropogenic activities that induce notable land subsidence [1–3]. All these human activities increase both the dynamic and static loading in soils which result in the compression of soil layers [4–6]. Significant land subsidence is induced by the superimposition effect of groundwater withdrawal, soil selfweight, and upper loads [7]. Groundwater drawdown is the main culprit that could significantly induce ground settlement by increasing the effective stress of soils [8]. Massive exploration of groundwater leads to the decreasing of the water table as well as the ground elevation [9–12]. Many researchers [13, 14] proved that the effect of soil sedimentation and continuous soil consolidation cannot be neglected. Ground movement also occurs with the change of the stress state due to upper loads from buildings and foundations, which could lead to severe volumetric changes within the soil [7]. The relationship between ground settlement and consolidation time under different ground water tables was also proposed in some literatures [15–18]. Li [19] did further study on ground settlement caused by soil selfweight, indicating that soil selfweight increases nonlinearly with depth. The impact of sedimentation on ground settlement increases with upper loads, soil compressibility, and soil thickness [20–22].
Analytical expressions were derived to calculate the excess pore water pressure and settlement caused by different types of loading. Shi et al. [23] proposed a 3D groundwater mode and a 1D deformation model to simulate the groundwater level and soil deformation. Some practical methods are also available, such as the relationships among effective soil weight and total settlement of normally consolidated clay [24]. With very distinguished soil properties, the SuXiChang area has experienced the most serious land subsidence due to excessive groundwater withdrawal in the late 1980s. More than 2 m cumulative land subsidence caused significant damage and heavy economic losses [25]. The SuXiChang area is also one of the most known in the world for the presence of earth fissures which has been observed within the areas and along the boundaries of these areas by many researchers [26, 27]. A novel approach to simulate the soil deformation has been introduced by Chen et al. [28]. From the geological investigation report of this area, groundwater is mainly stored in aquifers and some of them are stored in clay and silty clay layers due to the wide distribution of lowpermeability soil layers. Previous studies mainly focus on the homogeneous isotropic properties of hydrogeological parameters in soil layers, which is not suitable for the real cases [2, 5]. The variation of soil parameters in horizontal directions may cause difference in land subsidence which is the main concern in practical engineering.
The column element settlement model (CESM) is introduced in the present work to consider different characteristics of soil layers in horizontal direction. The study area is divided into several subareas in horizontal direction. Relative soil parameters are assigned to corresponding soil column according to the geological survey data. The accuracy of this proposed method was assessed by comparing numerical results with land subsidence records in study areas. Then the developed approach is applied to simulate the coupling effect of land subsidence induced by soil selfweight, upper loads, and groundwater withdrawal. The proportion and evolution of land subsidence induced by soil selfweight, upper loads, and groundwater withdrawal were discussed.
2. Methodology
2.1. Establishment of CESM
There are many limitations to consider the soil as homogeneous isotropic media because soils are always a threedimensional medium with the physical, mechanical, and hydrogeological parameter variations. CESM was proposed to depict the differences of hydrogeological structures, parameters, and physical mechanical index of soil layers in horizontal direction with the following properties: (1)CESM is a concept based on effective stress principles, that is, the total stress, effective stress, and pore water pressure are in a dynamic equilibrium. Soil compression increases with the increase of effective stress and decrease of pore water pressure(2)The column element shape can be either a cylinder or prism according to the geographic properties in the study area, and the soil parameters among each element are mutually independent (Figure 1). The settlement at the center of each element represents the settlement of this element. The accuracy of this model can be improved by the smaller element size. Soil inhomogeneity and anisotropy can be reflected by different partitions within the study area. The settlement distribution of the whole study area can be obtained through each element settlement(3)The sizes and shapes of each column can vary randomly based on the partitions of the study area. The soil parameters in vertical direction can be obtained from the borehole survey data, while the horizontal soil parameters of each column can be obtained by collecting the borehole data in different areas
(a) Cylindrical element
(b) Prismatic element
3. Solution of CESM
3.1. Materials and Parameters
3.1.1. Regional Geology and Hydrogeology
The SuXiChang area is located in the Yangtze River Delta, Jiangsu Province, southeastern China, bordering the west of Shanghai (Figure 2). It covers over a 17000 square kilometer area with a large population of over 21 million. The average elevation is less than 6 m above the sea level [23].
Due to the lack of in situ physical and mechanical parameters, the engineering geological analogy method was adopted to obtain the soil parameters of the study site [29]. The engineering geology analogy is mainly carried out based on the following aspects: regional location, lithological characteristics, and types of aquifer and aquitard.
Physical and mechanical parameters of Changzhou Xinzha, Wuxi Meilianghu Pumping Station Project, Taicang Danghe River hub project, Runyang Yangtze River Bridge south anchorage base, Taizhou Yangtze River Bridge Jiajiang base, and Sutong Yangtze River Bridge base were collected and applied to other areas with similar regional location, lithological characteristics, or types of aquifer and aquitard. The maximum and minimum values of the physical and mechanical parameters of the same stratum are collected and analyzed. The soil parameters of the SuXiChang area can be derived by the interpolation of obtained soil parameters in the study area.
Since the land subsidence is well known in the study area, several geological surveys have been conducted by the local Geological Survey Institute. Considering the soil parameters in Table 1, they have been derived from the Groundwater Resources Evaluation and Planning Report in the study area. Sixteen drilling pumping tests and geotechnical tests have been conducted to obtained the hydrogeological and physical mechanical parameters in the study area. The in situ test indicates that the mechanical properties of the soil in different elements share almost the same parameters in the same layer except two parameters, that is, transmissibility and water storage coefficient. Material parameters used in CESM numerical simulations are depicted in Table 1.

The whole study area is divided into 18 subareas in horizontal direction according to the soil characteristics from the Groundwater Resources Evaluation and Planning Report in Suzhou, Wuxi, and Changzhou (Figure 3). The transmissibility and water storage coefficient of different aquifers are calculated according to the drilling pumping test data. It can be seen from the table that the maximum transmissibility is 20 times larger than the minimum transmissibility and the maximum water storage coefficient is about 67 times the minimum storage coefficient (Table 2). It is also indicated from the table that the physical and mechanical properties of the aquifer as well as the hydrogeological parameters show a large discrepancy in different regions.

There is a wellestablished land subsidence station with multilayered extensometers to monitor the strata compaction in Changzhou. Thus, Changzhou is selected as the representative area to investigate the land subsidence caused by groundwater withdrawal, selfweight, and upper loading to lower the computing efforts. As it is well known, different hydrostratigraphic units have different kinds of deformation and the same unit may also exhibit different deformation characteristics, such as elasticity and elastoplasticity at different locations in different periods. Shi et al. [23] found that the deformation of both the sandy aquifer and the clay aquifer mainly possesses elastoplastic behavior. In this work, a simulation is performed based on the MohrCoulomb constitutive model and the soil performs elastoplastic behavior. According to the local geological soil properties of the study area, a conceptual model in Changzhou can be proposed by assuming as a plane strain problem [30]. Based on stratigraphic data and soil properties, the multilayer aquifer system consists of four aquifer units including four confined aquifers and four aquitard units. The left and right boundaries are set as a fixed water level with no horizontal displacement. The bottom boundary of the study area is assumed to be incompressible and impermeable.
The depth of the column element model is 156.5 m according to the average thickness of soil layers. A, B, C, D, E, F, G, H, and I located at the middle of each layer are selected as the monitoring points to indicate the settlement responses of each layer. The buried depths of the monitoring point from A to I are 3.75 m, 8.25 m, 15.25 m, 26.5 m, 40.25 m, 60.25 m, 85.25 m, 109 m, and 137.75 m, respectively. Five soil layers are considered through the soil profile of the study area. More specifically, the soil type is subclay, silty sand, mediumfine sand, gravel sand, and subsand from the top to the bottom soil layers. The crosssection of the hydrostratigraphy of the Changzhou area is presented in Figure 4.
3.2. Verification of CESM
Monitored data in Changzhou was adopted to verify the established CESM in the study area. Hierarchical signs for the monitored settlement in the Changzhou Qingliang primary school are shown in Figure 5.
(a)
(b)
3.2.1. Measured Settlement in Changzhou
Monitoring of ground settlement and piezometric head has been carried out since the end of the 20^{th} century. The depth of each monitoring point is shown in Table 3. Figure 6 depicted the cumulative ground settlement of the Changzhou Qingliang primary school from year 1984 to year 2000. Settlement at the end of each year was adopted to conduct comparison. Land subsidence increases with the investigation time and depth. Layer 3 and layer 4 account for most of the total settlement through the whole monitoring period.

3.2.2. Land Subsidence using CESM
The monitoring period in the Changzhou Qingliang primary school is from year 2000 to year 2005 for both the groundwater level and land subsidence, and thus, the effective comparison between monitored and calculated ground settlement is from year 2000 to 2005. The numerical model was established to achieve the ground settlement in Changzhou. The obtained soil parameters are used to predict the settlement of each layer.
The measured groundwater level was adopted to simulate the real cases. The upper load intensity is assumed to be 1 kPa. The ground settlement of each layer under different depths in partition #4 was selected to conduct comparison between the monitored settlement and calculated settlement. Figure 7 shows that the relative difference between the calculated settlement and measured settlement data can be neglected by given material parameters. The total settlement increases with the groundwater level. Monitored and calculated settlements from the ground surface to layer 3 show good agreement in years 2000, 2003, and 2004 (Figure 7). Settlements calculated in years 2001 and 2002 are lower than monitored settlement data. The settlement calculated in year 2005 is higher than the monitored settlement because of the upper load changes in the monitoring site at that year. This discrepancy may result from the uncertainties of soil parameters and the variation of groundwater levels in real cases. Calculated and monitored settlements from layer 4 to layer 7 exhibit the same trend which means that the numerical model confirms with the monitored data (Figure 8). Thus, the CESM can be adopted to obtain the land subsidence in the whole study area.
(a) Surface point
(b) Layer 1
(c) Layer 3
4. Results and Discussion
4.1. Ground Settlement Induced by SelfWeight
The ground settlement caused by soil selfweight is calculated by the layerwise summation method [31]. The aquifer in the study area can be divided into several layers of different thickness , and the deformations of each layer are , and the total ground settlement is where is the total ground settlement, is the ground settlement of each layer, is the number of soil layer, and it can be calculated as where is the compression coefficient of layer , is the void ratio before compression of layer , is the additional stress in layer (kPa), and is the thickness of layer .
The study area can be divided into 9 layers and 18 partitions according to the geological investigation. Unconfined aquifer, confined aquifer, and aquitard are distributed from top to the bottom. Settlement is different from layer to layer because of the variation of soil properties. The groundwater levels are 0 m, 7.5 m, and 156.5 m below the ground surface. The calculated ground settlement under different groundwater levels is shown in Table 4.

4.2. Evolution of Aquifer Settlement with Time
The differential equation governing onedimensional consolidation theory is as follows: where is the coefficient of consolidation of the soil, is the excess pore water pressure, is the distance measured downward from the surface of the consolidating layer, and is time [32].
The average degree of consolidation can be assumed as , 20%, 40%, 60%, 80%, and 90%. where is the dimensionless time factor and can be derived through and . , where is additional stress, is total stress, is the pore water pressure. is the average void ratio of soil during consolidation, and , , and are the thickness, permeability, and weight of soil, respectively. The coefficient of compressibility and equivalent coefficient of compressibility can be derived from the following equation [9, 15]: where is the equivalent coefficient of compressibility, which is calculated as 0.19 MPa^{1}. is the elastic modulus, is the thickness of each soil layer, and is the total thickness of the aquifer. Consolidation coefficient can be derived by substituting the parameters in Table 1 into the following equation: and the consolidation time can be expressed as follows:
The evolution of ground settlement under different groundwater depths (GWD) is depicted in Table 5. Ground settlement under soil selfweight increases with the increase of settling time and groundwater depth.

4.3. Ground Settlement under Different Upper Loads
Numerical analysis was conducted to investigate land subsidence induced by different upper loads. Two scenarios were considered as follows: (1)The groundwater table is assumed to be 0 m and load intensity is assumed to be 10 kPa, 25 kPa, 35 kPa, 50 kPa, and 80 kPa. Land subsidence of the study area under different load intensities is shown in Figures 9 and 10. It can be noted that the ground settlement increases gradually with the load intensity and time. Maximum land subsidence occurs at the surface of soil layers which indicates the depth effect of load intensity. Land subsidence decreases with the buried depth under the same load because the additional stress decreases with the soil depth, while it increases with the load intensity at the same monitoring point. The land subsidence increases linearly with the increase of loading time
Relationships between the upper loading and settlement of point A, point E, and point H in the unconfined aquifer, I_{lower} confined aquifer, and fourth aquitard were plotted in Figure 11. The loadsettlement relationship of point A, point E, and point H, i.e., equation (8), was also plotted in Figure 11. This relationship can be used to predict the land subsidence under different load intensities in the study area. The total land subsidence under different load intensities increases linearly with the load intensity and the depth, while the discrepancy between the real and predicted settlements may become larger under higher loading intensity due to plastic deformation of soil. where represents the magnitude of load intensity and represents land subsidence. (2)The groundwater table is assumed to be 4 m with a series variation of load intensity. Three scenarios were assumed, that is, soil selfweight changes linearly and load intensity increases first and then decreases to simulate the rising and falling of the groundwater table (Table 6)

As depicted in Figures 12–14, for case Π01, land subsidence increases with loading time and intensity and decreases with depth. For case Π02, settlement of soil layers increases with load and then decreases due to the change of load intensity. Settlement changes in a hysteresis way compared to loading. For point A, the ground settlement at the 4^{th} year is 4.2 mm less than the 7^{th} year under the same loading 40 kPa and the ground settlement at the 1^{st} year is 16.7 mm less than the 10^{th} year under the same loading 10 kPa. For case Π03, the land subsidence changes with loading in each soil layer.
4.4. Proportion of Land Subsidence Caused by Different Reasons
Three scenarios were considered including different groundwater levels and upper loadings. Land subsidence induced by different reasons in partition #4 was investigated to emphasize the proportion of land subsidence induced by selfweight, upper loads, and groundwater withdrawal. The land subsidence induced by different reasons are normalized by dividing total ground settlement. For example, the percentage of ground settlement caused by pumping is calculated by , where , , , and are the land subsidence induced by groundwater withdrawal, total ground settlement, the ground settlement induced by soil selfweight, and upper loads, respectively. (1)The groundwater level directly decreases from 0 m to 4 m below the ground surface; the loading intensity is 1 kPa
Comparison of settlement caused by selfweight, upper loading, and groundwater withdrawal is shown in Figure 15. Settlement caused by groundwater withdrawal and upper loading dominant from 0 to 200 years almost accounting for 97% of the total settlement. The settlement caused by selfweight gradually increases and reaches to 6.5% of the total settlement. However, settlement caused by drainage is higher than selfweight all the time. Figures 15(a) and 15(b) are the settlement proportions for different reasons at point A (depth 3.75 m) and point E (depth 40.25 m), respectively. Settlement caused by groundwater withdrawal and upper loading accounts for 98% of the total settlement from year 0 to 100. Settlement caused by selfweight is larger than that caused by pumping and loading after 150 years. The settlement caused by soil selfweight accounts for 70% of the total settlement. Analysis shows that settlement caused by groundwater drainage and upper loading dominates within 200 years and settlement caused by selfweight gradually prevails and becomes larger than that of groundwater withdrawal and loading with the increasing depth of aquifers. (2)The groundwater level gradually decreases from 0 m to 5 m with a velocity of 0.5 m per year; the loading intensity is 1 kPa
(a) Point A (depth 3.75 m)
(b) Point E (depth 40.25 m)
The comparison of settlement caused by selfweight, loading, and groundwater withdrawal is shown in Figure 16. As the groundwater table decreases, the settlement caused by groundwater withdrawal and soil selfweight increases, while loadinduced settlement decreases, resulting in an increasing percentage of ground settlement caused by pumping and selfweight (Figure 16(a)). Settlement caused by pumping and loading almost accounts for 95% of the total settlement at the first hundred years. Soil selfweight increases with the increase of depth. Settlement caused by soil selfweight reaches to 60% of the total settlement in a short time which exceeds the settlement caused by groundwater withdrawal and upper loading.
(a) Point A (depth 3.75 m)
(b) Point E (depth 40.25 m)
Settlement caused by groundwater withdrawal and selfweight increases with groundwater depth which results in the percentage of settlement caused by loading decreasing at the same time. Settlement caused by groundwater withdrawal has the highest percentage, settlement caused by loading is smaller, and settlement caused by selfweight is the smallest when the aquifer is buried shallow. Settlement caused by selfweight has the highest percentage and settlement caused by drainage is smaller while loading accounts for the smallest when the aquifer is deeply buried. (3)The initial upper loading is assumed to be 1 kPa and increases to the assumed maximum upper loading of 500 kPa (Table 7)

The upper loading increases gradually to investigate the settlement caused by loading intensity. The groundwater level in the II confined aquifer decreases from 0 m to 40 m below the ground surface in 100 years. Settlement caused by each factor and total settlement are shown in Figures 17–19.
(a) Point A (depth 3.75 m)
(b) Point E (depth 40.25 m)
(a) Point A (depth 3.75 m)
(b) Point E (depth 40.25 m)
(a) Point A (depth 3.75 m)
(b) Point E (depth 40.25 m)
Figure 17 shows that settlement caused by groundwater withdrawal can be much more significant than the other factors, i.e., the selfweight and loading intensity. The percentage of loadinginduced settlement decreases with the soil depth, as compared to those between point A and point E. Settlement caused by soil selfweight accounts for almost 30% of the total settlement when the soil is deeply buried. Figure 18 shows that settlement caused by loading and soil selfweight increases with the loading intensity and buried depth. Settlement caused by loading accounts for the largest percentage of total settlement when the loading magnitude is 500 kPa while settlement caused by groundwater withdrawal and selfweight is relatively small (Figure 19).
To sum up, land subsidence caused by loading increases with the loading intensity and decreases with the buried depth and the contribution of groundwater withdrawal and selfweight decreases with the loading intensity. Settlement caused by groundwater withdrawal accounts for most of the total settlement and increases with time at the depth of 40.25 m. And settlement caused by loading and soil selfweight decreases with the loading intensity at the same depth. The percentage of settlement caused by soil selfweight decreases with the loading intensity which is smaller than settlement caused by groundwater withdrawal. Settlement caused by soil selfweight accounts for the smallest percentage of total settlement at the depth of 3.75 m at any load intensities.
5. Conclusions
The column element settlement model is established to consider the soil anisotropic properties in horizontal direction. The numerical model was calibrated against the monitored data in the study area and then further applied to investigate the percentage of land subsidence induced by soil selfweight, upper loading, and groundwater withdrawal. The main conclusions are summarized below: (1)The column element settlement model is proposed, and soil parameters from field tests were adopted to the numerical model. The new model can be used to investigate the coupling effect of soil selfweight, upper loads, and groundwater withdrawal on land subsidence(2)The predicted land subsidence shows good agreement with the monitored results. Land subsidence caused by soil selfweight increases with time and groundwater depth. Land subsidence increases gradually with the load intensity and time(3)Normalized land subsidence was adopted to evaluate the percentage of land subsidence caused by different reasons. Land subsidence caused by groundwater withdrawal prevails through the whole monitoring period except for upper load intensities over 500 kPa. Settlement induced by soil selfweight and upper loads varies with different loadings and depths. The proportion of land subsidence caused by soil selfweight, upper loads, and groundwater withdrawal varies with examining time
The simulation procedure proposed in this study is quite generic and can be easily extended to simulate more sophisticated soil layer distributions. While a more appropriate soil constitutive model should be developed to satisfy different soil properties due to the complex deformation characteristics, which could be elastic, elastoplastic, and viscoelastic. The percentage of land subsidence induced by groundwater withdrawal is the main culprit for large ground deformation which deserves more attention in practical engineering.
Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that there is no conflict of interest regarding the publication of this paper.
Acknowledgments
This study was supported by the National Natural Science Foundation of China (Grant no. 51878127).
References
 Z. Cui and Y. Tang, “Domestic and international recent situation and research of land subsidence disasters,” Northwestern Seismological Journal, vol. 29, no. 3, pp. 275–278, 2007. View at: Google Scholar
 J. Hu, Z. Sui, and J. Chen, “Study on geological environmental effect after banning groundwater drainage in SuXiChang area,” Jiangsu Geology, vol. 30, no. 4, pp. 261–264, 2006. View at: Google Scholar
 K. S. Li, “Predicting the density profile and settlement of marine deposits,” in Proc. of International Conference on Computational Methods in Structural and Geotechnical Engineering Hong Kong, pp. 1562–1567, China Translation & Printing Services Ltd, 1994. View at: Google Scholar
 A. K. Gabriel, R. M. Goldstein, and H. A. Zebker, “Mapping small elevation changes over large areas: differential radar interferometry,” Journal of Geophysical Research, vol. 94, no. B7, pp. 9183–9191, 1989. View at: Publisher Site  Google Scholar
 H. Jiang, “Problems and discussion in the study of land subsidence in the SuzhouWuxiChangzhou area,” Quaternary Sciences, vol. 25, no. 1, pp. 29–33, 2005. View at: Google Scholar
 J. Qian and Z. Yin, Principle and Calculation of Geotechnical Engineering, China Water Conservancy and Hydropower Press, Beijing, 1996.
 P. Teatini, M. Ferronato, G. Gambolati, and M. Gonella, “Groundwater pumping and land subsidence in the EmiliaRomagna coastland, Italy: modeling the past occurrence and the future trend,” Water Resources Research, vol. 42, no. 1, 2006. View at: Publisher Site  Google Scholar
 J. Sundell, E. Haaf, T. Norberg, C. Alén, M. Karlsson, and L. Rosén, “Risk mapping of groundwaterdrawdowninduced land subsidence in heterogeneous soils on large areas,” Risk Analysis, vol. 39, no. 1, pp. 105–124, 2019. View at: Publisher Site  Google Scholar
 D. L. Galloway, The Frequency Dependence of Aquifer System Elastic Storage CoefficientsImplications for Estimates of Aquifer Hydraulic Properties and Aquifer System Compaction, U. S. Geological Survey, Sacramento, 1995.
 G. G. Meyerhof, “Shallow foundations,” Journal of Soil Mechanics & Foundations Division ASCE, vol. 19, no. SM2, pp. 21–31, 1965. View at: Google Scholar
 J. E. T. Moncur and R. L. Pollock, “Scarcity rents for water: a valuation and pricing model,” Land Economics, vol. 64, no. 1, pp. 62–72, 1988. View at: Publisher Site  Google Scholar
 X. Wang, Finite Layer Analysis of Ground Settlement Caused by Drainage, Hohai University, Nanjing, 2007.
 T. J. Burbey, “Stressstrain analyses for aquifersystem characterization,” Ground Water, vol. 39, no. 1, pp. 128–136, 2001. View at: Publisher Site  Google Scholar
 A. A. Fungaroli and S. R. Pranger, “Unit weights of a normally consolidated soil,” in Proc. of 7th international conference on soils Mechanics and Foundation Engineering, pp. 131–134, Mexico, 1969. View at: Google Scholar
 R. Bravo, J. R. Rogers, and T. G. Cleveland, Analysis of Groundwater Level Fluctuations and Borehole Extensometer Data from the Baytown Area, Houston, TX. Land Subsidence, International Association of Hydrological Sciences Publication 1218, Washington, 1991.
 J. B. Burland, M. Burbidge, E. Wilson, and TERZAGHI, “Settlement of foundations on sand and gravel,” Proceedings of the Institution of Civil Engineers, vol. 78, no. 6, pp. 1325–1381, 1985. View at: Publisher Site  Google Scholar
 R. W. Hexem and E. Heady, Water Production Function for Irrigated Agriculture, Iowa State University Press, 1978.
 R. L. Johnson and G. V. Johnson, Economic Valuation of Natural Resources: Issues, Theory and Applications, Westiview Press, Boulder. Colorado, 1990.
 G. Li, Advanced Soil Mechanics, Tsinghua University Press, Beijing, 2006.
 B. Li, Study of Large Strain Consolidation Theory in Layered Soft Clay Soil Foundation, Zhejiang University, Hangzhou, 1999.
 M. W. A. Van der Kooij, D. V. Halsema, W. Groenewoud et al., “Satellite Radar Measurements for Land Subsidence Detection,” in Proceedings of Land Subsidence, Barends and Schroder, Balkema, Rotterdam, 1995. View at: Google Scholar
 S. Ye and S. Dai, “Comparing of results of two dimensional, quasi three dimensional and three dimensional models for groundwater,” Hydrogeology & Engineering Geology, vol. 30, no. 5, pp. 23–27, 2003. View at: Google Scholar
 X. Shi, R. Fang, J. Wu, H. Xu, Y. Sun, and J. Yu, “Sustainable development and utilization of groundwater resources considering land subsidence in Suzhou, China,” Engineering Geology, vol. 124, pp. 77–89, 2012. View at: Publisher Site  Google Scholar
 V. K. Barwell and D. R. Lee, “Determination of horizontaltovertical hydraulic conductivity ratios from seepage measurements on lake beds,” Water Resources Research, vol. 17, no. 3, pp. 565–570, 1981. View at: Publisher Site  Google Scholar
 S. Ye, Y. Wang, J. Wu et al., “Characterization of earth fissures in South Jiangsu, China,” Proceedings of the International Association of Hydrological Sciences, vol. 372, pp. 249–253, 2015. View at: Publisher Site  Google Scholar
 G. Y. Wang, G. You, B. Shi, Z. L. Qiu, H. Y. Li, and M. Tuck, “Earth fissures in Jiangsu Province, China and geological investigation of Hetang earth fissure,” Environmental Earth Sciences, vol. 60, no. 1, pp. 35–43, 2010. View at: Publisher Site  Google Scholar
 Y. S. Xu, S. L. Shen, Z. Y. Cai, and G. Y. Zhou, “The state of land subsidence and prediction approaches due to groundwater withdrawal in China,” Natural Hazards, vol. 45, no. 1, pp. 123–135, 2008. View at: Publisher Site  Google Scholar
 Y. Chen, W. Zhao, J. Y. Han, and P. J. Jia, “A CEL study of bearing capacity and failure mechanism of strip footing resting on cφ soils,” Computers and Geotechnics, vol. 111, pp. 126–136, 2019. View at: Publisher Site  Google Scholar
 S. Qi, F. Wu, Z. Ding, and T. Liu, “Study on mechanism of deep fractures of Jinping first stage hydropower station by engineering geology analogy,” Chinese Journal of Rock Mechanics and Engineering, vol. 8, p. 033, 2004. View at: Google Scholar
 Z. Wang, Y. Zhang, J. Wu, J. Yu, and X. Gong, “Numerical simulation of earth fissures due to groundwater withdrawal,” Proceedings of the International Association of Hydrological Sciences, vol. 372, pp. 395–398, 2015. View at: Publisher Site  Google Scholar
 C. Peng and G. Yang, “A simplified method for determining ep curve of soft soil and Its application to analyzing nonlinear settlement of foundation,” Rock and Soil Mechanics, vol. 29, no. 6, pp. 1706–1710, 2008. View at: Google Scholar
 J. Lovisa, W. Read, and N. Sivakugan, “Time factor in consolidation: critical review,” International Journal of Geomechanics, vol. 13, no. 1, pp. 83–86, 2013. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2019 Yang Chen et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.