تعیین شدت بیابان‌زایی بر اساس شاخص‌های طیفی با استفاده از تصاویر سنتینل2 (منطقه مورد مطالعه: استان سیستان و بلوچستان)

پیشینه و هدف سطوح مختلف پوشش گیاهی آلبدوی متفاوتی دارند. از طرفی آلبدوی سطحی یکی از مهم‌ترین مولفه‌های تعادل تابش سطحی است که با بررسی رفتار آن می‌توان به شدت تخریب و بیابان‌زایی پی برد. پوشش گیاهی به دلیل اینکه عاملی برای پایداری سطح زمین محسوب می‌شود می‌تواند یکی از مهم‌ترین مولفه‌های کلیدی در مناطق خشک برای کاهش اثرات فرسایش و بیابان‌زایی به‌حساب آید. گسترش بیابان‌زایی و تغییر در میزان پوشش گیاهی به‌نوبه خود ازجمله عوامل تغییر در آلبدو به شمار می‌روند. هدف از این مطالعه تعیین شدت بیابان‌زایی بر اساس شاخص‌های طیفی آلبیدو (albedo)، اندازه دانه خاک سطحی (tgsi) و شاخص نرمال‌شده تفاضل پوشش گیاهی (ndvi) با استفاده از فناوری سنجش‌ازدور می‌باشد. شناسایی مناطق تخریب شده در کوتاه‌ترین زمان و با کمترین هزینه با استفاده از تصاویر سنتینل2 با قدرت تفکیک مکانی 10 متری از اهداف این پژوهش می‌باشد. در این پژوهش برای اولین بار در منطقه سیستان و بلوچستان بر اساس شاخص‌های طیفی با استفاده از تصاویر سنتینل2، بهترین شاخص برای پایش شدت بیابان‌زایی در مناطق خشک معرفی خواهد شد.مواد و روش ها برای ارزیابی شدت بیابان‌زایی و شناخت شاخص مناسب برای تهیه نقشه شدت بیابان‌زایی گام‌های زیر انجام گرفت؛ 1) انتخاب تصاویر و انجام عملیات پیش‌پردازش تصاویر با استفاده از نرم‌افزار snap. 2) محاسبه شاخص‌های tgsi، ndvi و albedo، 3) بررسی رابطه همبستگی بین شاخص‌ها با استفاده از نرم‌افزار spss®24، 4) تهیه نقشه شدت بیابان‌زایی منطقه و به دست آوردن معادله شدت بیابان‌زایی با استفاده از نرم‌افزار arcgis®10.3. در گام اول این پژوهش داده‌های ماهواره sentinel2a مربوط به سنسور msil1c 20 اوت سال 2020 انتخاب شد. تصاویر به‌گونه‌ای انتخاب گردید که فصل رویش گیاهان یک‌ساله و موقتی نباشد و همچنین روزی انتخاب شد که پوشش ابری وجود نداشته باشد تصاویر مورد نیاز از سایت http://scihub.copernicus.eu/ دانلود و مورد استفاده قرار گرفت.نتایج و بحث نتایج حاصل از مدل رگرسیون خطی بین دو شاخص ndvi و albedo نشان داد که این دو شاخص با یکدیگر دارای همبستگی منفی می‌باشند و به ترتیب میزان ضریب همبستگی در منطقه سوران و زابل برابر با 0.76 و 0.63 بود. نتایج نشان داد که با افزایش میزان شاخص ndvi از میزان شاخص آلبدو کاسته می‌شود. همچنین نتایج حاصل از مدل رگرسیون خطی بین دو شاخص tgsi و albedo نشان داد که این دو شاخص با یکدیگر رابطه قوی و مثبتی داشته و به ترتیب میزان ضریب همبستگی برای منطقه سوران و زابل برابر با 0.78 و 0.81 بود. نتایج نشان داد که با افزایش میزان شاخص tgsi بر میزان شاخص آلبدو افزوده می‌شود. شدت بیابان‌زایی در مناطق مورد مطالعه بر اساس معادله  تعیین گردید و با استفاده از روش شکست طبیعی جنکس (natural breaks jenks) در نرم‌افزار arcgis شدت بیابان‌زایی به 5 درجه؛ مناطق بدون بیابان‌زایی، مناطق با شدت کم بیابان‌زایی، مناطق با شدت متوسط بیابان‌زایی، مناطق با شدت زیاد بیابان‌زایی، و مناطق با شدت خیلی زیاد بیابان‌زایی تقسیم گردید. در این پژوهش بر اساس داده‌های ماهواره سنتینل2 شاخص‌های albedo، ndvi و tgsi استخراج گردید. نتایج بررسی رگرسیون خطی بین دو شاخص ndvi و albedo نشان داد که یک رابطه منفی و قوی بین این دو شاخص وجود دارد که با نتایج پژوهش‌های مشابه مطابقت دارد. نتایج حاصل از همبستگی بالا و منفی بر اساس این دو شاخص به این مفهوم است که هر گونه افزایش در مقدار شاخص پوشش گیاهی ndvi منجر به کاهش در میزان آلبدوی سطح خواهد شد. از طرفی مناطق با آلبدوی بالا بیانگر تخریب پوشش گیاهی و برهنه بودن خاک می‌باشد. در مناطقی که کلاس شدت زیاد بیابان‌زایی مشاهده گردید مقدار شاخص آلبدوی سطحی بالا می‌باشد و مقدار شاخص پوشش گیاهی کم است. طبقه‌بندی شدت بیابان‌زایی در منطقه سیستان بر اساس مدل  albedondvi نشان می‌دهد که 27.73 درصد منطقه در کلاس بدون شدت بیابان‌زایی، 18.03 درصد در کلاس کم شدت، 32.92 درصد از منطقه از نظر بیابان‌زایی در کلاس شدت متوسط بیابان‌زایی، 20.3 درصد در کلاس شدید و تنها 1.02 درصد از منطقه در کلاس خیلی شدید بیابان‌زایی قرار گرفته است. همچنین طبقه‌بندی شدت بیابان‌زایی در سوران بر اساس مدل albedondvi نشان می‌دهد که 4.82 درصد منطقه بدون شدت بیابان‌زایی، 8.44 درصد در کلاس کم، 50.97 درصد از منطقه از نظر بیابان‌زایی در کلاس شدت متوسط، 34.48 درصد در کلاس شدید و 1.3 درصد از منطقه در کلاس خیلی شدید بیابان‌زایی قرار گرفته است. بیشترین درصد شدت بیابان‌زایی مربوط به کلاس شدت متوسط می‌باشد. نتایج رگرسیون خطی بین دو شاخص tgsi و albedo نیز نشان داد که یک رابطه مثبت و قوی بین این دو شاخص وجود دارد. نتایج بیانگر این است که رابطه بین شاخص tgsi و albedo نسبت به رابطه بین دو شاخص ndvi و albedo قوی‌تر و از ضریب همبستگی بالاتری در هر دو منطقه برخوردار می‌باشد که از دلایل عمدۀ آن می‌توان به پراکندگی پوشش گیاهی در مناطق خشک اشاره نمود. رابطه شاخص tgsi با albedo خصوصیات فضایی مناطق عاری از پوشش گیاهی و همچنین مناطق با پوشش گیاهی خیلی کم را برای تعیین شدت بیابان‌زایی بهتر نشان می‌دهد. شاخص tgsi منعکس کننده اندازه ذرات درشت خاک سطحی می‌باشد که رابطه مثبتی با ذرات ریز محتوای ماسه خاک سطحی دارد. هر چه اندازه ذرات خاک سطحی درشت دانه‌تر باشد میزان بیابان‌زایی شدت بیشتری خواهد داشت. در مناطقی که محتوای زیاد ماسه ریز در اندازه ذرات خاک سطحی وجود داشته باشد مقادیر بالای شاخص tgsi قابل مشاهده خواهد بود.نتیجه گیری در این پژوهش ما با استفاده از تکنیک سنجش از راه دور و بر اساس تصاویر چند طیفی سنتینل2 برای اولین بار در ایران به استخراج شدت بیابان‌زایی در دو منطقه متفاوت از استان خشک سیستان و بلوچستان پرداختیم. بر اساس بازتاب طیفی اتفاق افتاده از سطح زمین و قدرت تفکیک مکانی 10 متری شدت بیابان‌زایی را در دو منطقه مورد مطالعه قرار دادیم. بر اساس نتایج این پژوهش پیشنهاد می‌گردد به‌صورت ترکیبی از دو مدل albedondvi و albedotgsi برای پایش و تهیه نقشه‌های شدت بیابان‌زایی در مناطق خشک ایران استفاده گردد. نتایج این پژوهش نشان داد که مناطق بدون بیابان‌زایی و شدت کم بیابان‌زایی بر اساس مدل albedotgsi به نحو مطلوب‌تری نشان داده می‌شوند.

Determining the desertification intensity based on spectral indices using Sentinel2 images (Case study: Sistan and Baluchestan province)

Background and Objective Different vegetation covers have different albedo levels. On the other hand, surface albedo is one of the most important components of surface radiation balance, which can be used to identify severely degraded and desertified regions. Vegetation can be considered as one of the most important key components in arid regions to reduce the effects of erosion and desertification due to the effects of vegetation for land surface stability. Expansion of desertification and also changes in vegetation cover, could be change the surface Albedo. The purpose of this study is to determine the desertification intensity based on spectral indices, Albedo, Topsoil Grain Size Index (TGSI) and Normalized Difference Vegetation Index (NDVI) using remote sensing technology. Identification the damaged areas with the lowest cost in the shortest time, using Sentinel2 images with a spatial resolution of 10 meters is one of the objectives of this study. Also, this study will introduce the best indicator for monitoring desertification intensity in arid regions for the first time in the Sistan and Baluchestan region based on spectral indices using Sentinel2 images.Materials and Methods The following steps were performed to evaluate the intensity of desertification and identify the appropriate indicator in order to mapping the desertification intensity: 1) Selection the images and perform image preprocessing operations using SNAP software; 2) Calculation of TGSI, NDVI and Albedo indices; 3) Investigation the correlation between indices using SPSS®24 software. 4) Preparation of desertification intensity map of the region and obtaining the equation of desertification intensity using ArcGIS®10.3 software. In the first step of this research, Sentinel2A satellite data related to MSIL1C sensor was selected on August 20, 2020. The images were selected in such a way that the growing season of the plants is not annual and temporary, and also the day was selected when there is no cloud cover. The required images were downloaded and used from the URL address: http://scihub.copernicus.eu/. Results and Discussion The results of linear regression between NDVI and Albedo indices showed that, these two indices had negative correlation, and the correlation coefficient in Souran and Zabol was 0.76 and 0.63, respectively. The results showed that with increasing NDVI, decreased of the albedo index occurred. Also, the results of linear regression model showed strong and positive relationship between TGSI and Albedo indices, as the correlation coefficient of Souran and Zabol was 0.78 and 0.81, respectively. The results showed that the TGSI and the albedo simultaneously decreased or increased. Desertification intensity in the study areas was determined based on the equation I= a × Index ± Albedo and also by using Natural Breaks (Jenks) method in ArcGIS software, desertification intensity of study areas classified to 5 degrees, 1. Without desertification, 2. Low desertification, 3. Moderate desertification, 4. Severe desertification, and 5. Extremely desertification. In this study Albedo, NDVI and TGSI indices were extracted based on Sentinel2 satellite data. The results of linear regression between NDVI and Albedo showed that there is strong negative relationship between these indices that was consistent with the results of similar studies. The high and negative correlation, means that any increase in the vegetation cover will lead to decrease the Albedo. On the other hands the areas with high Albedo, indicate degradation of vegetation cover and bare soil. In the regions with sever desertification intensity, the value of surface Albedo was high and the vegetation cover was low. Classification of desertification intensity in Sistan region based on AlbedoNDVI model showed that 27.73% of the area were in the class of without desertification intensity, 18.03% in the low class, 32.92% in the moderate class, 20.3% were in the severe class and only 1.02% of the area were in the very severe desertification intensity class. Also, the classification of desertification intensity in Souran based on AlbedoNDVI model showed 4.82% of the area without desertification, 8.44% in low class, 50.97% in moderate class, 34.48% in severe class and 1.3% of the area were in very severe desertification class. The highest percentage of desertification intensity of the area were in the moderate class. The results of linear regression between TGSI and Albedo indices also showed that there is a positive and strong relationship between these indices. The results showed that the relationship between TGSI and Albedo indices was stronger than the relationship between NDVI and Albedo indices and in both regions the correlation coefficient was higher. One of the main reasons for this is the dispersion of vegetation cover in arid areas. The relationship between TGSI and Albedo better shows the spatial characteristics of vegetationfree areas as well as areas with very low vegetation cover to determine the intensity of desertification. The TGSI index reflects the coarse particle size of the topsoil, which has a positive relationship with the fine sand content of the topsoil. Whatever the larger particle size of the topsoil, will have the greater desertification intensity. In the areas where the content of fine sand in the topsoil is high, the high range of TGSI index will be seen.Conclusion In this study, using Sentinel2 multispectral images and remote sensing technique, we extracted the intensity of desertification in different arid regions of the Sistan and Baluchestan province, for the first time in Iran. Based on the spectral reflection that occurred from the ground and the spatial resolution of 10 meters, we studied the intensity of desertification in two areas. Based on the results of this research, we suggest to use the combination of AlbedoNDVI and AlbedoTGSI models in order to monitoring the desertification intensity in arid regions of Iran. The results of this study showed that areas without desertification and low intensity of desertification are better identified based on AlbedoTGSI model.