پیش‌بینی تغییرات کاربری اراضی با مدل زنجیره مارکوف و سلول خودکار (مطالعه موردی: منطقه هیرکانی)

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

نویسندگان

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

2 گروه علوم و مهندسی جنگل، دانشکده کشاورزی و منابع طبیعی و علوم دریایی، دانشگاه تربیت مدرس، نور، ایران

3 گروه مرتع و آبخیزداری، دانشکده کشاورزی و منابع طبیعی، دانشگاه محقق اردبیلی، اردبیل، ایران

چکیده

سابقه و هدف: مدل‌سازی و پیش‌بینی تغییرات کاربری اراضی جهت استفاده پایدار از زمین و آگاهی از میزان تغییرات احتمالی در آینده امری ضروری است. آشکارسازی و مدل‌سازی این تغییرات با پردازش تصاویر ماهواره‌ای ابزار سودمندی جهت درک تغییرات محیط زیستی مرتبط با فعالیت‌های انسانی به شمار می‌روند. استفاده از سامانه اطلاعات جغرافیایی و سنجش از دور، اطلاعات به‌هنگامی را پیرامون پدیده‌های سطح زمین ارائه می‌دهد. نظارت بر روند این تغییرات به منظور استفاده پایدار از زمین می‌تواند با استفاده از داده‌های سنجش از دور چند زمانی با موفقیت حاصل شود. با توجه به اهمیت زیست‌محیطی منطقه رویشی هیرکانی، تنوع ژنتیکی و کاربردهای اکولوژیکی و گردشگری آن، آگاهی از وسعت تغییرات و تبدیل کاربری‌های مختلف آن مهم‌ترین جنبه برنامه‌ریزی کاربری اراضی و توسعه پایدار است. همچنین، از آنجایی‌که تا کنون پایش تغییرات کاربری اراضی در گذشته و پیش‌بینی آن در آینده در کل منطقه هیرکانی صورت نگرفته است، این پژوهش، با تلفیق مدل زنجیره مارکوف و سلول‌های خودکار تحت عنوان "مارکوف-CA" به شبیه‌سازی و پیش‌بینی تغییرات کاربری اراضی منطقه هیرکانی از گذشته تا آینده می‌پردازد.
مواد و روش‌ها: برای آشکارسازی و مدل‌سازی تغییرات کاربری اراضی، تصاویر لندست 5 سنجنده TM سال‌های 2000، 2005 و 2010 و لندست 8 سنجنده OLI سال‌های 2015 و 2020 میلادی تجزیه و تحلیل و مورد استفاده قرار گرفتند. پس از اعمال فیلترهای مکانی، زمانی و ابری (ابرناکی کمتر از 5 درصد)، تصاویر تمام مقاطع زمانی با استفاده از فیلتر میانگین در پلتفرم گوگل ارث انجین (GEE) ترکیب شدند و سپس، با استفاده از الگوریتم یادگیری ماشین نظارت شده جنگل تصادفی در شش کلاس (جنگل، مرتع، کشاورزی، مناطق ساخته شده، اراضی بایر و پهنه‌های آبی) طبقه‌بندی شدند. دقت طبقه‌بندی بر اساس شاخص کاپا و دقت کلی محاسبه شد. شبیه‌سازی کاربری اراضی سال 2020 با استفاده از نقشه‌های کاربری سال‌های 2000 و 2010 و به کمک مدل تلفیقی سلول‌های خودکار و قوانین انتقال رگرسیون لوجستیک در نرم‎افزار MATLAB 2021a انجام گرفت. برای این منظور از متغیرهای مکانی شامل شاخص نرمال شده تفاوت پوشش گیاهی (NDVI)، شاخص تفاضلی نرمال شده آب (NDWI)، شاخص نرمال شده تفاوت مناطق ساخته شده (NDBI)، ارتفاع و شیب به‎عنوان عوامل مؤثر بر تغییرات در الگوریتم جنگل تصادفی استفاده شد. در راستای انتخاب نمونه‌ها برای کاربری‌های در نظر گرفته شده، سعی شد تا نمونه‌ها پراکنش همگن و مناسبی در سطح منطقه داشته باشند. از 70 درصد نمونه‌های برداشت شده به صورت بصری از روی تصویر ماهواره‌ای، برای اجرای مدل (نمونه‌های تعلیمی یا آموزشی) و از 30 درصد باقیمانده برای ارزیابی دقت مدل (نمونه‌های تست یا اعتبارسنجی) استفاده شد. همچنین، مدل مارکوف-CA برای پیش‌بینی تغییرات کاربری اراضی در منطقه هیرکانی برای سال‌های 2070 و 2100 تلفیق شدند.
یافته‌ها: بر اساس ارزیابی صحت داده‌های مکانی، برای سال‌های 2000 تا 2020، مقادیر دقت کلی و ضریب کاپا بالاتر از 80 درصد به‌دست آمد که حاکی از انطباق زیاد تصاویر پیش‌بینی شده مدل با واقعیت زمینی و کارایی سامانه GEE در پردازش و طبقه‌بندی تصاویر لندست می‌باشد. در دوره 2000 تا 2010، کاهش 8/3 درصدی در پوشش جنگل و تبدیل آن به اراضی مرتعی و کشاورزی رخ داد. در دوره 2010 تا 2020 بیشترین میزان تبدیل از کاربری مرتع به کشاورزی به وسعت 872,78 هکتار بود. با مدل سلول‌های خودکار و قوانین رگرسیون لجستیک، شاخص کاپا برای پیش‌بینی کاربری‌ها در سال‌های 2070 و 2100 به ترتیب 76 و 72 درصد به‌دست آمد. ضریب کاپای حاصل از مقایسه این دو نقشه با مقدار 31/87 درصد نشان‌دهنده قابلیت بالای مدل مارکوف-CA در شبیه‌سازی تغییرات کاربری اراضی برای آینده بود. نتایج پیش‌بینی تغییرات کاربری اراضی با استفاده از تخمین‌گر انتقال مارکوف مبیّن تخریب و کاهش سطح کاربری‌های جنگل، مرتع و پهنه‌های آبی در سال‌های 2070 و 2100 بود. با تغییر پوشش کاربری‌های مختلف بین سال‌های 2020 تا 2070، سطح اراضی جنگلی نسبت به دوره 2020-2000، 168,680 هکتار کاهش خواهد یافت. سطح پهنه‌های آبی نیز تا سال 2070 کاهش خواهد یافت و به میزان 8,667 هکتار به کاربری مرتع تبدیل خواهد شد. تا سال 2100 و طی یک دوره زمانی 30 ساله، سطح اراضی جنگلی 8,900 هکتار کاهش خواهد یافت و به پهنه‌های آبی و اراضی بایر تبدیل خواهد شد. بیشترین تغییرات کاربری اراضی در منطقه مورد مطالعه شامل تبدیل جنگل‌ها به مراتع و اراضی کشاورزی، تبدیل مراتع به کاربری کشاورزی و مناطق ساخته‌شده، و نیز تبدیل پهنه‌های آبی به اراضی بایر بوده است. در واقع، در صورت تداوم روند تغییر کاربری اراضی و تبدیل غیر اصولی آن‌ها به یکدیگر، عرصه‌های منابع طبیعی به‌شدت کاهش خواهند یافت.
نتیجه‌گیری: در مجموع، صحت ارزیابی بر اساس شاخص‌های مستخرج از ماتریس خطا مبیّن این بود که سامانه GEE با پردازش داده‌های لندست در سری‌های مختلف زمانی توانست ابزار مهمی جهت ارزیابی و پایش تغییرات کاربری باشد. روند تغییرات در منطقه مورد مطالعه به گونه‌ای بود که کاربری‌های کشاورزی، ساخته شده و بایر روند افزایشی داشتند که دلیل این امر، تبدیل سایر کاربری‌ها به آن‌ها بود. در مقابل، در مقایسه با جنگل و مرتع، بیشترین روند کاهشی در منطقه مربوط به پهنه‌های آبی بود که ادامه این روند تخریب، پیامدهای نامطلوب زیست‌محیطی را به دنبال خواهد داشت. افزایش سطح کاربری‌های مناطق ساخته شده و کشاورزی را می‌توان به در توجیه افزایش جمعیت و گسترش شهرنشینی برای رفع نیازهای غذایی این جمعیت و همچنین مدیریت صحیح در زمینه کشاورزی در نظر گرفت. از آنجاییکه با بررسی تغییرات کاربری اراضی، روندهای مربوط به فرآیندهای تخریب، جنگل‌زدایی، بیابان‌زایی و فقدان تنوع زیستی در یک منطقه مشخص می‌گردد، با آگاهی از این تغییرات در طول زمان، می‌توان چالش‌های آینده را پیش‌بینی کرد و با ارائه اقدامات پیشگیرانه، مدیریت بهتری را اعمال نمود.

کلیدواژه‌ها

موضوعات


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

Prediction of Land Use Changes Using the Markov Chain and Cellular Automata Model (Case Study: Hyrcanian Region)

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

  • Sara Hedayati Kaliji 1
  • Seyed Mohsen Hosseini 1
  • Seyed Jalil Alavi 2
  • محدثه امیری 3
1 Department of Forest Sciences and Engineering, Faculty of Natural Resources and Marine Sciences, Tarbiat Modares University, Noor, I.R. Iran
2 Department of Forest Sciences and Engineering, Faculty of Natural Resources and Marine Sciences, Tarbiat Modares University, Noor, I.R. Iran
3 Department of Range and Watershed Management, Faculty of Agriculture and Natural Resources, University of Mohaghegh Ardabili, Ardabil, I.R. Iran
چکیده [English]

Background and Objectives: Modeling and predicting land use changes is essential for sustainable land use and awareness of potential future changes. Detection and modeling of such changes using satellite image processing serve as effective tools for understanding environmental transformations associated with human activities. The use of Geographic Information Systems (GIS) and Remote Sensing (RS) provides systematic and accurate information about land surface phenomena. Monitoring the trend of these changes for sustainable land use can be successfully achieved using multi-temporal remote sensing data. Given the ecological significance of the Hyrcanian region- due to its genetic diversity and ecological and touristic functions- understanding the extent of land use changes and conversions is a key aspect of land use planning and sustainable development. Moreover, since land use change monitoring in the past and its future prediction have not yet been conducted across the entire Hyrcanian region, this study aims to simulate and predict land use/land cover (LULC) changes from past to future by integrating the Markov chain model with Cellular Automata (CA), known as the "Markov-CA" approach. Another novelty of this study lies in the implementation of image processing and modeling workflows within the cloud-based Google Earth Engine (GEE) environment.
Materials and Methods: To detect and model land use changes, Landsat 5 TM images from the 2000, 2005 and 2010, and Landsat 8 OLI images from 2015 and 2020 were analyzed and used. After applying spatial, temporal, and cloud filters (with cloud cover less than 5%), all time-series images were combined using a mean filter in the Google Earth Engine platform, and then classified into six classes (i.e. forest, rangeland, agriculture, built-up areas, bare land, and water bodies) using the supervised Random Forest machine learning algorithm. Classification accuracy was calculated based on the Kappa index and overall accuracy. To simulate land use for the year 2020, land use maps from 2000 and 2010 were utilized in conjunction with a combined Cellular Automata model and logistic regression-based transition rules, implemented in MATLAB 2021a. Spatial variables including Normalized Difference Vegetation Index (NDVI), Normalized Difference Water Index (NDWI) and Normalized Difference Built-up Index (NDBI) indices, as well as elevation and slope, were used as influencing factors on changes in the Random Forest algorithm. To ensure the representativeness of the selected land use classes, efforts were made to achieve a homogeneous and well-distributed sampling across the study area. Seventy percent of the visually interpreted samples from satellite imagery were used for model implementation (training samples), while the remaining 30 percent were employed for accuracy assessment (test or validation samples). Additionally, the Markov-CA model was integrated to predict land use changes in the Hyrcanian region for the years 2070 and 2100.
Results: Based on the spatial accuracy assessment, for the period 2000 to 2020, the overall accuracy and Kappa coefficient were higher than 80%, indicating a high agreement between the model-predicted images and ground reality and also the efficiency of the Google Earth Engine platform in processing and classifying Landsat imagery. During the 2000-2010 period, a 3.8% decrease in forest cover was observed, primarily due to its conversion into rangelands and agriculture. In the subsequent period (2010–2020), the most significant land use conversion was from rangeland to agriculture, covering an area of 78,872 ha. Using the cellular automata model and logistic regression rules, the Kappa coefficient for land use predictions in 2070 and 2100 was estimated at 76% and 72%, respectively. The Kappa coefficient obtained from comparing these two maps was 87.31%, demonstrating the high capability of the Markov-CA model in simulating future land use changes. The predicted land use changes using the Markov transition estimator showed the degradation and reduction in forest, rangeland, and water bodies in 2070 and 2100. As a result of changes in land use/land cover between 2020 and 2070, forested areas are expected to decrease by 1,686,80 ha compared to the 2000-2020 period. Water bodies are also projected to decline by 2070, with approximately 8,667 ha expected to be converted into rangelands. By 2100, over a 30-year period, forest areas are projected to decrease by 8,900 ha, primarily transitioning into water bodies and barren lands. The most significant land-use changes in the study area include the conversion of forests to rangelands and agriculture, the transformation of rangelands into agriculture and built-up areas, and the conversion of water bodies into barren lands. In fact, if the current trend of land use changes and improper conversions continue, natural resource areas will be drastically reduced.
Conclusion: Overall, the accuracy assessment based on indices derived from the error matrix indicated that the Google Earth Engine platform, through processing multi-temporal Landsat data, proved to be an effective tool for evaluating and monitoring land use changes. The observed trend in the study area showed that agricultural lands, built-up areas, and bare lands had an increasing trend, primarily due to the conversion of other land use types into these classes. In contrast, among all land-cover types, water bodies experienced the most significant decline compared to forests and rangelands. The continuation of this degradation trend is expected to lead to adverse environmental consequences. The increase in built-up areas and agricultural lands can be justified by population growth and urbanization to meet the food needs of this population, as well as proper management in agriculture. Since by examining land use changes, trends related to degradation, deforestation, desertification, and loss of biodiversity in a region are identified, awareness of these changes over time, it is possible future challenges and implement preventive measures for better management.

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

  • Change monitoring
  • Markov transition estimator
  • Image classification
  • Markov cellular automata model
  • GEE