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

نویسندگان

1 دانشجوی دکتری ژئومورفولوژی دانشگاه حکیم سبزواری

2 استادیار دانشگاه حکیم سبزواری

3 استاد دانشگاه حکیم سبزواری

10.22034/gp.2021.11306

چکیده

فرونشست زمین در طی سال‌های گذشته به دلایل مختلف از جمله رشد بی رویه جمعیت، برداشت بیش از حد از منابع آب زیرزمینی، عوامل تکتونیکی به همراه عوامل دیگر سبب بروز مشکلات و معضلات فراوان به زمین‌های کشاورزی، ساختمان‌های مسکونی، جاده‌ها، خطوط انتقال نیرو و ... شده است. این پدیده در چند سال گذشته آسیب‌های فراوان به دشت‌های استان فارس به خصوص ساختمان‌های مسکونی و زمین‌های کشاورزی وارد نموده است. پایش میزان نرخ فرونشست و نیز پرداختن به علت یا علل تاثیر گذار آن جهت کنترل و مدیریت خطر دارای اهمیت فراوان است. در این پژوهش به منظور بررسی وضعیت فرونشست درآبخوان دشت نورآباد ممسنی از داده های سنتینل 1A در بازه زمانی (2018-2015) و روش تداخل سنجی راداری که با مقایسه فازهای دو تصویر راداری که از یک منطقه در دو زمان مختلف اخذ شده اند، قادر به تعیین تغییرات سطح زمین در آن بازه ی زمانی می باشد در محیط نرم افراز( SNAP) انجام گرفت. جهت تهیه نقشه های افت و سطح آب زیرزمینی در سال‌های حداقل 1382 و حداکثر1394 از روش درونیابی IDW استفاده شد و سپس میان دو پارامتر افت آب زیرمینی و فرونشست سطح دشت با استفاده از مدل GWR تحلیل فضایی صورت پذیرفت و خروجی حاصل از این تحلیل نشان داد که حداکثر میزان R2 در بخش های مرکزی آبخوان 79 /0 می باشد و در بخش های جنوبی این مقدار کاهش می یابد. همچنین میزان مقادیر برآورد شده با پیش بینی شده نشان می دهد که وضعیت خطا در بخش مرکزی بین 0006 .0تا 0.0007 می باشد. نشان از اعتبار سنجی بالای مدل و میزان دقت بالای بین خروجی های حاصل از دو پارامتر می باشد.

کلیدواژه‌ها

موضوعات

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

Spatial analysis of land subsidence and groundwater loss using the GWR model (Case study: Noorabad Mamasani aquifer,)

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

  • najmeh shafiei 1
  • LAILA GOLIMKOTARI 2
  • Aboalghasem Amir Ahmadi 3
  • Rahman Zandi 2

1 PhD student of geomorphology of Hakim Sabzevari University

2 Assistant professor Hakim Sabzevari University

3 Professor Hakim Sabzevari University

چکیده [English]

Introduction
Today, the phenomenon of land subsidence is one of the most important geomorphological hazards on a global scale, which causes great damage to urban and rural structures, development facilities. Most of the reports from around the world related to the phenomenon of land subsidence have been related to arid and low rainfall areas. According to the US Geological Survey, a subsidence phenomenon involves a collapse or subsidence that can have a small displacement vector. Land subsidence is a geological phenomenon that causes the earth's surface to descend slowly and horizontally.
In recent decades, most joints and fissures in agricultural areas have been reported due to over-abstraction of groundwater. Recently, urban areas have also been affected by this phenomenon. These seams and cracks cause a lot of annual damage to ground-level installations such as roads and bridges, buildings, power transmission lines, oil and gas transmission pipes, water pipes and sewage systems, wells wall pipes of exploitation wells. In addition, they provide a path for surface pollutants to move to groundwater sources and contaminate groundwater by infiltrating aquifers.
Among the innovations of this study, for the first time, radar interferometry technique was used to investigate the risk of subsidence of Nurabad plain aquifer subsidence. , Identification of vulnerable areas in the aquifer area.
Methodology
The method used in this research is applied analysis. Groundwater data including piezometer and exploitation wells were used to investigate groundwater level changes and Sentinel 1 images were used to calculate radar interferometry. Geographic weighted regression model (GWR) was used to investigate the relationship between subsidence and groundwater loss parameters.
By examining the changes in groundwater level of piezometer wells in the aquifer during a period of 17 years (1380-1397) showed that an average of 15 meters of groundwater drop has occurred in the plain. This decrease in agricultural and residential areas is a priority compared to other uses. By preparing the output of radar images during 4 statistical years, the results show that in 2015 the maximum amount of leakage was 7 cm and in 2016 it increased to 8 cm. In 2017 it is equal to 9 cm and in 2018 it is equal to 10 cm. In the end, by calculating the average of these four years, the amount of subsidence at the aquifer level increases significantly, with the regions in the central and eastern parts of the aquifer having the highest subsidence at the aquifer level. Which has been developed in residential and agricultural areas of the region so that the effects and evidence of these meetings are visible in residential houses and agricultural lands according to Table (2) in the analysis section to examine the status of subsidence rates in the study area using Radar images were taken during the mentioned years and show that in 2015 it is equal to 8 square kilometers, which is compared to 2014, in 2016 it is equal to 34 square kilometers, and in 2017 and 2018 it is equal to 40 and 86 square kilometers. According to the maps, we see in these 4 consecutive years that every year b The amount of subsidence at the plain level has increased and thus subsidence at the plain level has become a kind of hazard.
Results and discussion
 High values ​​of R2 have occurred in the central parts of the model aquifer, which shows a good estimate of the model in estimating the dependent variable and the predictive explanatory variable. Also, relatively lower values ​​are located in the southern parts of the aquifer by examining the subsidence status in the plain. By creating a relationship with the groundwater drop layer, it showed that the observed values ​​with the predicted values ​​indicate a strong relationship. The layer shows the coefficient of determination R2, which is the highest coefficient of R2 in the central part, which is about -0.79 Is 0.53 Figure (7). The distribution of the remaining space of the GWR model errors shows that the model outputs are closer to the actual values. Second, the resulting map shows the low error values ​​in the range.
Occlusion
The average subsidence rate is about 4 cm per year and the average 4-year average is 16 cm. The average area of ​​the aquifer where the subsidence has taken place during 4 years is about 36 square kilometers, which includes about 5% of its area. The high concentration of wells and the greatest drop in groundwater and the concentration of existing cities and villages. In these areas, due to the water needs of the residential community, including the city of Nurabad and surrounding areas, as well as the water needs of the industrial sector, including factories and agricultural activities, access to water resources is more and groundwater extraction is very high. There has been a high decline that these factors with the intervention of the arid climate of this region has caused the high talent of these areas in the discussion of subsidence and has been identified as one of the critical areas of the city in the subsidence map. The results of radar interferometry show that uncontrolled extraction from the surface of underground reserves has caused the subsidence of the central part of the aquifer and also the presence of several aqueducts at the aquifer has caused damage to residential centers. One of the areas with the highest risk of subsidence in the plain, in the eastern, central and southern parts, the rate of groundwater loss in the aquifer within 14 statistical years is about 22 meters in the same well and the results of GWR method show local R2 in the central parts of the aquifer which is about 79% maximum and shows a high correlation between subsidence and groundwater loss Residual rate in the central part which is equal to 0.007 / 0007 which is a high coefficient of significance.

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

  • land subsidence
  • Groundwater drop
  • Interfemteri
  • GWR
  • Noorabad Mamasani Aquifer
-       حافظی، مقدس ناصر، غفوری، محمد ،قزی، اعظم،(1389)، مشکلات قنات­های قدیمی در شهر مشهد، همایش زمین شناسی مهندسی ساخت گاه‌ها و ابنیه‌های فنی و تاریخی، صص 536-531.
-       رکنی، جعفر، حسین زاده، رضا، لشکری پور غلامرضا، ولایتی سعد الله، (1395)، بررسی فرو نشست زمین، چشم اندازها و تحولات ژئومورفولوژی ناشی از آن در دشت­های ترا کمی مطالعه موردی: دشت نیشابور،  مطالعات جغرافیایی مناطق خشک، (24)6، 38 -21.
-       زندی، رحمان، فرزین کیا ربابه، شفیعی، نجمه، (1398)، فرونشست و تداخل سنجی راداری، انتشارات ماهواره.
-       سازمان آب منطقه­ای استان فارس،گزارش بهنگام سازی اطلس منابع آب حوضه آبریز رودخانه­های هندیجان و جراحی، (1398)، جلد4، صص216-1.
-       شایان، سیاوش، یمانی، مجتبی، یادگاری، منیژه(1395)، پهنه‌بندی فرونشست زمین در حوضه آبخیز قره‌چای همدان، هیدروژئومورفولوژی، شماره 9، صص158-139.
-       شریفی کیا، محمد،(۱۳۹۱)، تعیین میزان و دامنه فرونشست به کمک تداخل سنجی راداری در دشت نوق - بهرمان، مجله برنامه ریزی و آمایش فضا، دوره شانزدهم، شماره ۳ .
-       شفیعی، نجمه، گلی مختاری، لیلا، امیراحمدی، ابوالقاسم، زندی، رحمان،(1398)، بررسی فرونشست آبخوان نورآباد با استفاده از روش تداخل سنجی راداری، نشریه پژوهش­های کمی ژئومورفولوژی،
-       عطایی، محمد، (1390)، معدن کاری زیرزمینی، جلد 3، انتشارات دانشگاه صنعتی شاهرود.
-       قضاوی، رضا، رمضانی سربندی، مجید(1396)، بررسی تأثیر تغییرات میزان بارش و برداشت از آب‌های زیرزمینی بر تغییرات کمی و کیفی آب آبخوان (مطالعه‌ی موردی: دشت رفسنجان)، فصلنامه هیدروژئومورفولوژی، شماره 12 صص 111-129.
-        Bates, R. L, and Jackson, J. A., 1980, “Glossary of Geology.” American Geological Institute. Second edition, Falls Church, Virginia, 749 p.
-        Gallowey, D.L., Burbey; T.J., 2011, Regional land subsidence accompanying groundwater extraction, Hydrogeol, J., Vol. 19, pp1459-1486.
-        Gabriel, A.K., R.M. Goldstein and H.A. Zebker, 1989. Mapping small elevation changes over large areas - Differential radar interferometry. J. Geophys. Res., 94(B7), pp. 9183-9191.
-        Helm, D.C., Looney, B.B. and Falta, R.W., 2000, Fissures in Yucca Flat dry lake bed, Nevada test site. Vadose Zone, Science and Technology Solutions. Battelle Press, Columbus, OH [CD-ROM.].
-        Faunt, C., Sneed, M., Traum, J., Brandt, R., 2015, Water availability and Land Subsidance in the Central Valley, California USA, Hydrogeol journal, spring 15(3),pp.452-470.
-        Khorrami M, Alizadeh B, Ghasemi Tousi E, Shakerian M, Maghsoudi M, Rahgozar P. 2019. How Groundwater Level Fluctuations and Geotechnical Properties Lead to Asymmetric Subsidence: A PSInSAR Analysis of Land Deformation over a Transit Corridor in the LosAngeles Metropolitan Area. Remote Sensing journal, 11(4):257-270.
-        Larson,K.J., Basaglu,H., Marino,M.A, 2001, Predication  of optimal safe ground water yield and Land subsidence in Los Banos Kettleman citty area, using a calibrated numerical simulation model, journal Hydrology, Vol 242, pp79-102
-        Lin, B., Liming, J., Hansheng, W., Qishi, S. 2016, Spatiotemporal Characterization of Land Subsidence and Uplift (2009–2010) over Wuhan in Central China Revealed by TerraSAR-X InSAR Analysis.
-        Putra, D.P.E. Setianto, A. Keokhampui, K. andFukuoka, H, 2011,Land Subsidence Risk Assessment CaseStudy: Rongkop, Gunung Kidul, Yogyakarta-Indonesia.The 4th AUN/SEED‐Net Regional Conference on Geo‐Disaster Mitigation in ASEAN“Past tragedies are the lessonsfor future mitigation”. The Royal Paradise Hotel & Spa,Phuket, Thailand.
-        Raspinia, F, Loupasakis., Constantinos; Rozosb, Dimitrios; Adamc, Nico; Moretti, Sandro, 2014, Ground subsidence phenomena in the Delta municipality region (Northern Greece): Geotechnical modeling and validation withPersistent Scatterer Interferometry, ELSEVIER International Journal, vol28, pp78-89.
-        Smith R G, Knight R , Chen J. Reeves J A , Zebker H A , Farr T , Liu  Z.2017. Estimating the permanent loss of groundwater storage in the southern San Joaquin Valley, California, Water Resources Research journal,)53(2): 147-163.
-        Sun.H, Grandstaff.D, Shagam.R, (1999). Land subsidence due to groundwater Withdrawal: potential damage of subsidence and sea level ruse in southern New Jersy, USA, Environmental Geology 37(4) April 1999 7 Q Springer-Verlag. pp 290-296.
-        Tobler, W.R, 1970, A Computer movie simulating urban growth in the Detroit region, Economic Geography, 46(2): 234-24.
-        Zebker, H. A., )2000(. “Studying the earth with interferometric radar.” Computing in Science and Engineering, Vol. 2(3), pp. 52-60.