تحلیل احتمالاتی فرونشست زمین ناشی از برداشت آب زیرزمینی تحت شرایط غیرماندگار در خاک‌های ریزدانه

نوع مقاله : مقاله پژوهشی

نویسندگان

1 گروه مهندسی آب، دانشکده علوم کشاورزی، دانشگاه گیلان، رشت، استان گیلان، ایران

2 عضو هیات علمی گروه مهندسی آب دانشگاه گیلان

چکیده

فرونشست زمین به‌عنوان یکی از پیامدهای برداشت بی‌رویه آب‌های زیرزمینی، تهدیدی جدی برای زیرساخت‌ها و محیط‌زیست به شمار می‌رود. در این تحقیق، مدل‌سازی فرونشست ناشی از برداشت آب زیرزمینی با استفاده از یک رویکرد تحلیل احتمالاتی انجام شده است. بدین منظور، یک برنامه رایانه‌ای در محیط نرم‌افزار MATLAB بر پایه تابع تاثیر اصلاح‌شده بودریک–کنات توسعه داده شد و برای مدل‌سازی همبستگی میان پارامترهای ورودی، از چندین تابع کاپولا استفاده گردید. در تحقیق حاضر، با توجه به اهمیت نقش عدم‌قطعیت در فرآیند مدل‌سازی فرونشست، دو رویکرد شامل اعمال عدم‌قطعیت موجود در متغیرهای زمانی و متغیرهای مرتبط با خصوصیات مکانیکی خاک مورد بررسی و مقایسه قرار گرفت. نتایج مطالعه موردی انجام‌شده در منطقه فومنات استان گیلان با نوع خاک رسی سیلتی نشان داد که در شرایط اعمال عدم‌قطعیت متغیرهای زمانی، با گذشت زمان دامنه تغییرات فرونشست محاسبه‌شده کاهش یافته و مدل رفتار پایدارتری نشان می‌دهد. در مقابل، با اعمال عدم‌قطعیت موجود در خصوصیات مکانیکی خاک، افزایش پراکندگی مقادیر فرونشست در طول زمان مشاهده می‌شود. بر این اساس، یافته‌ها بیانگر این است که در ارزیابی‌های کوتاه‌مدت فرونشست باید بیشتر به خصوصیات مکانیکی خاک تاکید گردد، در صورتی که عدم‌قطعیت‌های مبتنی بر زمان در بلندمدت تاثیر غالب ایفا می‌کنند. به طور کلی، نتایج بیانگر آن است که برای خاک رسی سیلتی، تابع کاپولا فرانک مناسب‌ترین توزیع احتمال مشترک در شرایط اعمال عدم‌قطعیت متغیرهای زمانی و خصوصیات مکانیکی خاک می‌باشد.

کلیدواژه‌ها

موضوعات


عنوان مقاله [English]

Probabilistic Analysis of Land Subsidence Caused by Groundwater Extraction under Unsteady Conditions in Fine-Grained Soils

نویسندگان [English]

  • Maryam Nejati Kalashami 1
  • Amir Malekpour 2
  • Behnam Shafiei Sabet 1
1 Dept. of Water Engineering, Faculty of Agricultural Sciences, University of Guilan, Rasht, Guilan Province, Iran
2 Academic staff, Dept. Water Engineering, University of Guilan
چکیده [English]

Land subsidence, as a consequence of the excessive groundwater extraction, poses a serious threat to infrastructure and environment. In the current research, the land subsidence caused by groundwater withdrawal was investigated and modeled using a probabilistic analysis approach. A computer program was developed in MATLAB based on the modified Budryk–Knothe influence function and several copula functions were employed to model the dependence among input parameters. Given the significance of incorporating uncertainty in subsidence modeling, two approaches were considered and compared including the application of temporal random variables and the random variables related to mechanical soil properties. A case study was conducted in the Foomanat region of Guilan Province, characterized by silty clay soil. The results demonstrated that when the temporal random variables were considered, the variations in the range of calculated subsidence gradually decreased over time, indicating a more stable subsidence model. In contrast, incorporating the uncertainties associated with mechanical soil properties led to an increased scattering in subsidence values with time. Accordingly, the findings suggest that for short-term subsidence assessments, greater emphasis should be placed on soil mechanical properties, whereas temporal uncertainties become more dominant in long-term analyses. Generally, for the silty clay soil of the study area, the Frank copula was identified as the most appropriate function for modeling the joint probability distributions when the temporal random variables and random soil properties are considered.

کلیدواژه‌ها [English]

  • Consolidation settlement
  • copula
  • influence function
  • pumping
  • uncertainty

Introduction

Land subsidence, as one of the most critical consequences of excessive groundwater extraction, poses significant risks to infrastructure, agriculture, and the environment. Given the growing demand for groundwater in recent years, studying and predicting this phenomenon—particularly in vulnerable regions—has become increasingly important. Time plays a pivotal role in subsidence processes and the long-term behavior of soil, and overlooking its impact can lead to less accurate predictions. A thorough analysis of this phenomenon requires accounting for the inherent uncertainties of the influential parameters. A deeper understanding of these uncertainties, along with the role of time, can improve predictive models and contribute to a more effective groundwater resource management in at-risk areas.

Materials and Methods

The objective of this research is to investigate and model the land subsidence resulting from the excessive groundwater extraction using a probabilistic analysis approach. To achieve this goal, a computer program was developed in the MATLAB environment based on a modified influence function. Copula functions were also employed to model the correlations among the input parameters. In the Fumanat region of Guilan province characterized by silty clay soil, given the significance of incorporating uncertainties in subsidence modeling, two distinct sets of variables were examined including temporal random variables and random mechanical soil properties. Additionally, a sensitivity analysis was conducted to identify the most influential soil parameters affecting the subsidence.

Results and Discussion

The results indicated that in the first approach, incorporating the temporal parameters, both the range of variability and the coefficient of variation of subsidence decreased over time. The peak of the probability density function was lowest during the early and mid-periods, which gradually increased in the long term. Analyzing the performance of copula functions revealed that the Gumbel copula exhibited the highest coefficient of variation during short and medium terms, whereas the Clayton and Frank copula provided a lower coefficient of variation. In the second approach, where the uncertainties of soil mechanical properties were investigated, the range of data variability increased over time, while the coefficient of variation remained relatively constant. Moreover, the maximum values of probability density functions showed a descending trend. Comparing the performances of copulas indicated that the Frank copula yielded the lowest coefficients of variation. These results underscore the importance of selecting an efficient copula function in the complex analyses related to land subsidence.

Conclusion

This study provides compelling evidence regarding the significance of considering uncertainties in subsidence models. Since the probabilistic analysis is applied in contexts with a high level of uncertainty, it is recommended that the short-term subsidence investigation is better to focus on mechanical properties of soil, while in the long-term analyses the temporal variables should be more considered. Furthermore, the findings indicate that the Frank copula provides the most appropriate joint probability distributions for both the temporal variables and soil properties.

Funding

This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors.

Authorship contribution

All authors contributed equally to the conceptualization of the article and writing of the original and subsequent drafts.

Data availability statement

Data will be available on request from the authors.

Acknowledgements

The authors would like to thank all the people assisted in the present study including experts from the Regional Water Company of Guilan province and deputy for research and technology of the University of Guilan.

Ethical considerations

The authors avoided data fabrication, falsification, and plagiarism, and any form of misconduct.

Conflict of interest

The authors declare no conflict of interest.

 

Aghayi Gil Chalan, M. (2023). Sensitivity analysis and application of copula functions in the probabilistic analysis of subsidence induced by groundwater extraction using influence function method. (Masters Thesis, Guilan University, Rasht).(In Persian)
Aghayar, M. S., Malekpour, A., & Shafiei Sabet, B. (2020). Varied groundwater pumping-induced subsidence in stratified soil condition using Influence Function and GIS-derived data. Iranian Water Researches Journal, 14(4), 153-162. (In Persian)
Alibeikloo, M., Khabbaz, H., & Fatahi, B. (2022). Random field reliability analysis for time-dependent behaviour of soft soils considering spatial variability of elastic visco-plastic parameters. Reliability Engineering & System Safety, 219, 108254. https://doi.org/10.1016/j.ress.2021.108254 .
Bagheri-Gavkosh, M., Hosseini, S. M., Ataie-Ashtiani, B., Sohani, Y., Ebrahimian, H., Morovat, F., & Ashrafi, S. (2021). Land subsidence: A global challenge. Science of The Total Environment, 778, 146193. https://doi.org/10.1016/j.scitotenv.2021.146193.
Bai, Y., Wang, Y., Chen, Y., & Zhang, L. (2020). Probabilistic analysis of the controls on groundwater depth using Copula Functions. Hydrology Research, 51(3), 406-422. https://doi.org/10.2166/nh.2020.147.
Budhu, M., & Adiyaman, I. B. (2010). Mechanics of land subsidence due to groundwater pumping. International journal for numerical and analytical methods in geomechanics, 34(14), 1459-1478. https://doi.org/10.1002/nag.863.
Cai, Y., Verdel, T., & Deck, O. (2014). On the topography influence on subsidence due to horizontal underground mining using the influence function method. Computers and Geotechnics, 61, 328-340. https://doi.org/10.1016/j.compgeo.2014.06.003.
Chen, M., Tomás, R., Li, Z., Motagh, M., Li, T., Hu, L., ... & Gong, X. (2016). Imaging land subsidence induced by groundwater extraction in Beijing (China) using satellite radar interferometry. Remote Sensing, 8(6), 468. https://doi.org/10.3390/rs8060468.
Chu, S. T. (1994). Transient Radius of influence model. Journal of irrigation and drainage engineering, 120(5), 964-969. https://doi.org/10.1061/(ASCE)0733-9437(1994)120:5(964).
Choubari, R. Y., Estalaki, A. M., & Varaki, M. E. (2025). Bivariate probabilistic slope stability analysis using copulas and back analysis by Markov chain Monte Carlo. Earth Science Informatics, 18(3), 291. https://doi.org/10.1007/s12145-025-01789-x .
Davydzenka, T., Tahmasebi, P., & Shokri, N. (2024). Unveiling the global extent of land subsidence: The sinking crisis. Geophysical Research Letters, 51(4), e2023GL104497. https://doi.org/10.1029/2023GL104497.
Don, N. C., Hang, N. T. M., Araki, H., Yamanishi, H., & Koga, K. (2006). Groundwater resources management under environmental constraints in Shiroishi of Saga plain, Japan. Environmental geology, 49, 601-609. https://doi.org/10.1007/s00254-005-0109-9.
Fallahchai, M. M. (2016). Natural Regeneration of (Populus caspica Bornm.) and its relationship with soil physical and chemical properties (Case Study: Safrabaste Region in East of Guilan province). Journal of Plant Research (Iranian Journal of Biology), 29(1), 118-129. (In Persian)
Galloway, Devin L., & Burbey, Thomas J. (2011), Review: regional land subsidence accompanying groundwater extraction, Hydrogeology Journal, 19( 8), 1459–1486. https://doi.org/10.1007/ s10040-011-0775-5.
Ghabraie, B., Ren, G., Barbato, J., & Smith, J. V. (2017). A predictive methodology for multi-seam mining induced subsidence. International Journal of Rock Mechanics and Mining Sciences, 93, 280-294. https://doi.org/10.1016/j.ijrmms.2017.02.003
Gholami Estalkhi Kohi, T., Navabian, M., & Eemaeili Varaki, M. (2022). Investigation of the Dimensions of Open Drains in Paddy Fields of Land Consolidation Project. Iranian Journal of Water Research in Agriculture (Formerly Soil and Water Sciences), 36(3), 313-328. https://sid.ir/paper/1058041/en. (In Persian)
Hung, W. C., Hwang, C., Liou, J. C., Lin, Y. S., & Yang, H. L. (2012). Modeling aquifer-system compaction and predicting land subsidence in central Taiwan. Engineering Geology, 147, 78-90. https://doi.org/10.1016/j.enggeo.2012.07.018.
keyvan behjou, F., Masarat, F., Ghanbari, S., & Sasanifar, S. (2024). The effect of logging activities on the soil characteristics of Islam forests in Guilan. Journal of Environmental Science Studies, 9(3), 9158-9147. http://dx.doi.org/10.22034/JESS.2024.423944.2164 . (In Persian)
Khademy Kamal, M., Naghdi, R., Salehi, A., & Yousefpour, R. (2023). Changes evaluation of some soil properties under effects of management characteristics (Case study: forests of Gisom region of Guilan province). Forest and Wood Products, 75(4), 309-319. https://dor.isc.ac/dor/20.1001.1.50522008.1401.75.4.1.1 (In Persian)
Kazempour Larsari, Z., Esmaeili Varaki, M., & Malekpour, A. (2019). Laboratory study of scour downstream of stepped-labyrinth weirs. Iranian Journal of Soil and Water Research 49 (6), 1227-1241. (In Persian) https://doi.org/10.22059/ijswr.2018.246313.667796
Malekpour, A., Farookhroo, S., & Hosseini, M. (2019). Effect of Uncertainty of Soil Compressibility on Temporal Variations of Pore-Water Pressure at the Foundation of Structure. Water and Soil Science, 29(2), 199-212. (In Persian)
Malekpour, A., Niknam Yousefi, H., & Esmaili Varaki, M. (2018). Slope Stability Analysis after Decline of Water Level at the River Meander Considering Strain and Monte-Carlo Probabilistic Method. Water and Soil Science, 28 (2), 1-13. (In Persian)
Malekpour, A., Sadeghian, N., & Farrokhi, M. J. (2023). Bivariate probabilistic analysis of temporal variations of pore water pressure during consolidation process in structural foundation. Irrigation and Drainage Structures Engineering Research, 24(90), 91-113. https://doi.org/10.22092/idser.2023.363197.1551. (In Persian)
Malinowska, A., Hejmanowski, R., & Dai, H. (2020). Ground movements modeling applying adjusted influence function. International journal of mining science and technology, 30(2), 243-249. https://doi.org/10.1016/j.ijmst.2020.01.007.
Mehdidoost Roudbaneh, I., Janatrostami, S., Ashrafzadeh, A., & Javadi, S. (2023). Accuracy assessment of groundwater recharge estimation using SWAT and MODFLOW in paddy fields (Case study: Astane-Kouchsefahan aquifer). Iranian Journal of Soil and Water Research, 54(10), 1541-1564. https://doi.org/10.22059/ijswr.2023.363264.669547 . (In Persian)
Perżyło, D. (2023). A new model for forecasting of land surface subsidence caused by underground extraction of deposits. International Journal of Environmental Science and Technology, 20(4), 3733-3748. https://doi.org/10.1007/s13762-022-04225-5.
Ren, G., & Li, J. (2008). A study of angle of draw in mining subsidence using numerical modeling techniques. Electron. J. Geotech. Eng, 13, 1-14.
Ren, G., Buckeridge, J., & Li, J. (2015). Estimating land subsidence induced by groundwater extraction in unconfined aquifers using an influence function method. Journal of Water Resources Planning and Management, 141(7), 04014084. https://doi.org/10.1061/(ASCE)WR.1943-5452.0000479.
Rezai Poor Jolandan, S., Pourbabaei, H., Salehi, A. (2022). Effect of Conservation on Physical and Chemical Properties of Gisum Forest Park Soil and Dr. Dorostkar Warehouse in province Guilan. Natural Ecosystems of Iran ,12(4), 1-12. (In Persian)
Sadegh, M., Ragno, E., & AghaKouchak, A. (2017). Multivariate C opula A nalysis T oolbox (MvCAT): describing dependence and underlying uncertainty using a B ayesian framework. Water Resources Research, 53(6), 5166-5183. https://doi.org/10.1002/2016WR020242.
Sadeghian, N., Malekpour Estalaki, A., & Calamak, M. (2024). Probabilistic internal erosion analysis in stratified and unstratified foundations of embankment dams using copulas. Natural Hazards, 120, 12989–13007. https://doi.org/10.1007/s11069-024-06722-x .
Saeidi, A., Deck, O., Seifaddini, M., Heib, M. A., & Verdel, T. (2022). An improved methodology for applying the influence function for subsidence hazard prediction. Georisk: Assessment and Management of Risk for Engineered Systems and Geohazards, 16(2), 347-359. https://doi.org/10.1080/17499518.2021.1875247.
Sepúlveda-García, J. J., & Alvarez, D. A. (2022). On the use of copulas in geotechnical engineering: A tutorial and state-of-the-art-review. Archives of computational methods in engineering, 29(7), 4683-4733. https://doi.org/10.1007/s11831-022-09760-5.
Wang, T., Zhou, G., Wang, J., & Wang, D. (2020). Impact of spatial variability of geotechnical properties on uncertain settlement of frozen soil foundation around an oil pipeline. Geomechanics and Engineering, 20(1), 19-28. https://doi.org/10.12989/gae.2020.20.1.019 .
Zahmatkesh Biazar, F. (2018). Effect of uncertainty of local geotechnical properties on probability of land subsidence under groundwater pumping condition. (Masters Thesis, Guilan University, Rasht).(In Persian)
Zarinibahador, M., Nabiollahi, K., & Norouzi, M. (2015). Influence of different slope aspects on some soil properties and forest soils evolution (Case study: Rostam abad region, Guilan province). Water and Soil, 29(3), 648-662. https://doi.org/10.22067/jsw.v0i0.31149. (In Persian)