استخراج شاخص رطوبت سطحی خاک (Tvdi) با استفاده از نمودار پراکندگی دما/ پوشش گیاهی و تصاویر مودیس

پیشینه و هدف رطوبت خاک به عنوان یک پارامتر مهم در کنترل بسیاری از فرآیندهای سیستم آب و هوایی، یکی از پارامترهای بنیادی محیط زیست و تاثیر مستقیم آن بر زندگی گیاهی، جانوری و میکرواورگانیسم‌ها، اهمیت آن در چرخه جهانی آب، انرژی و کربن، تبادلات انرژی بین هوا و خاک، چرخه آبی طبیعت (به ویژه در توزیع باران بین رواناب سطحی و نفوذ) و مدیریت منابع آب و خاک شناخته می شود. رطوبت خاک در فرآیندهای تعاملی بین جو و زمین و تغییرات جهانی اقلیم نقش مهمی ایفا می کند. روش های مثلثی و ذوزنقه ای با ترکیب داده های حرارتی و مرئی، از جمله روش‌های پرکاربرد در تعیین میزان رطوبت سطحی خاک می‌باشند. هدف از انجام این تحقیق برآورد رطوبت سطحی خاک، با شاخص خشکی دما و پوشش گیاهی (tvdi)، به روش مثلثی جنوب استان آذربایجان غربی با استفاده از شاخص دمای سطح زمین (lst) و شاخص پوشش گیاهی (ndvi)، طی سال های 2010 ، 2014 و 2018، است. مواد و روش ها مطالعه حاضر با استفاده از تصاویر سری زمانی سنجنده مودیس، شاخص نرمال‌شده تفاضل پوشش گیاهی (ndvi) و شاخص تغییرات دمای سطح زمین (lst)، برای برآورد شاخص رطوبت سطحی (شاخص خشکی دما و پوشش گیاهی، tvdi)، در سه بازه زمانی شامل؛ بازه زمانی اول 11 دی 1388 تا 9 دی 1389، بازه زمانی دوم 11 دی 1392 تا 9 دی 1393 و بازه زمانی سوم  11 دی 1396 تا 9 دی 1397 انجام شد. در هر بازه زمانی 12 تصویر در روز 15 ام هر ماه مورد استفاده قرار گرفت. همچنین برآورد رطوبت سطحی به دو روش یکی برقرار کردن رابطه رگرسیونی و برداشتن حداقل دما و روش دوم برقرار کردن رابطه رگرسیونی و پایین پیکسل ها انجام شد. برای ارزیابی دقت این دو روش از همبستگی رگرسیونی بین نتایج حاصل از این روش‌ها با نمونه زمینی رطوبت سطحی خاک مربوط به جهاد کشاورزی (30 نقطه) در عمق 5 تا 15 سانتی متر استفاده شد. علت انتخاب کردن این سه سال به دلیل تفاوت بارندگی زیاد در بعضی ماه های سال های مورد مطالعه است. این تحقیق که با استفاده از تصاویر مودیس برای بررسی میزان رطوبت سطحی خاک، در جنوب استان آذربایجان غربی در بخشی از منطقه غرب ایران با استفاده از روش مثلثی به دو صورت متفاوت انجام شد.  نتایج و بحث نمودار مثلث تبخیری متشکل از شاخص پوشش گیاهی و دمای سطح زمین از 11دی ماه 1388 تا 9 دی ماه 1389 تغییرات دمای زیادی رخ داده است. همین تغییراتی دمای سطح زمین باعث شده است، که نمودارها تغییرات زیادی داشته باشند. در سال 1388 طبق نمودار، حداکثر دما مردادماه و حداقل آن دی ماه و حداکثر پوشش گیاهی در اردیبهشت ماه و حداقل آن آذرماه بوده است. در سال 1392، حداکثر دما در مردادماه و حداقل آن دی ماه و حداکثر پوشش گیاهی در اردیبهشت ماه و حداقل آن دی ماه و همچنین این سال نسبت به سال 1388 نسبتاً گرمتر و خشک تر بوده است. نمودار مثلث تبخیری در سال 1398 نسبت به دو سال دیگر مورد مطالعه هم میزان بارش بیشتر بوده و هم میزان پوشش گیاهی و با توجه به گراف ها در این سال حداکثر دما در تیرماه و حداقل آن دی ماه و حداکثر پوشش گیاهی در اردیبهشت ماه و حداقل آن در دی ماه است. میزان رطوبت سطحی خاک در سال 1388 برای منطقه غرب کشور ایران که حداکثر میزان رطوبت در اردیبهشت ماه و حداقل آن در مردادماه است. در بیشتر نقشه های شاخص رطوبتی در سال 1388، حداکثر رطوبت در غرب و حداقل آن در جنوب منطقه است. نتایج نقشه های شاخص رطوبتی در سال 1392 در این سال نسبت به سال های دیگر مورد مطالعه نسبتاً خشک‌تر بوده است. سال 1392 از بارش و پوشش گیاهی کمی برخوردار است. تغییرات رطوبت در این سال نسبت به سال 1388 کمتر می باشد. حداکثر و حداقل رطوبت در سال 1392 بین 0 تا 0.6 است. حداکثر رطوبت خردادماه و حداقل آن مردادماه می باشد. نقشه های شاخص رطوبتی tvdi سال 1396 در این سال نسبت به دو سال دیگر مورد مطالعه میزان شاخص رطوبتی بیشتر داشته است. در سال 1396 بارش زیاد باعث شده، که پوشش گیاهی افزایش و دمای سطح زمین کاهش یافته و این مسئله باعث افزایش میزان شاخص رطوبتی نسبت به سال 1388 و 1392 شده است. در سال 1396 میزان پوشش گیاهی تا حد 0.89 می رسد. اما در سال های دیگر مورد مطالعه تا 0.7 بوده است. در این سال، بیشترین میزان رطوبت مربوط به اردیبهشت ماه و کمترین آن در مردادماه قرار دارد. حداکثر رطوبت در طول این سال در غرب و کمترین مقدار آن جنوب است. نتایج شاخص خشکی دما و پوشش گیاهی (tvdi) برای سال های 1388، 1392 و 1396 با استفاده از روش دوم نتایج کلی این روش شبیه به روش اول می باشد. با توجه به نتایج بدست آمده از میزان دقت دو تا روش، دقت روش اول بهتر و در کل ساده تر از روش دوم می باشد. در سال 1396 در اردیبهشت ماه با توجه به روش اول مقدار 0.66=r2 و با توجه به روش دوم مقدار 0.41=r2است. نتیجه گیری برآورد رطوبت خاک سطحی برای مدیریت بهینه منابع آب و خاک ضروری است. رطوبت خاک سطحی، متغیری مهم در چرخه آبی طبیعت است که نقش مهمی در تعادل جهانی آب و انرژی به واسطه تاثیر بر فرآیندهای هیدرولوژیک، اکولوژیک و هواشناسی دارد. بررسی دو تا روش مورد استفاده حاکی از آن می باشد که روش اول که به طور کلی در این تحقیق نیز مورد استفاده قرار گرفته بود با توجه به نتایج دقت تصاویر نسبت به داده های زمینی دقت بالاتری برخوردار است. در سال 1388 اردیبهشت ماه و مردادماه با توجه به روش اول مقدار r2 به ترتیب 0.61 و 0.57 است. در سال 1388 مقدار r2 با توجه به داده های زمینی و استفاده از روش دوم در اردیبهشت ماه و مردادماه به ترتیب 0.43 و 0.47 می باشد. همچنین، در سال 1396 مقدارr2 با استفاده از روش اول در اردیبهشت ماه 0.66 است. در سال 1396 مقدارr2 با استفاده از روش دوم در اردیبهشت ماه 0.41 است. نتایج به دست آمده از مدل برآورد شاخص رطوبتی سطحی خاک در این پژوهش نشان داد که این مدل قادر به برآورد مقدار رطوبت خاک در پهنه های وسیع جغرافیایی با دقت قابل قبول است.

Extraction of soil moisture index (TVDI) using a scatter diagram temperature/vegetation and MODIS images

Background and Objective Soil moisture is an important parameter in controlling many processes of the climate system, one of the basic parameters of the environment and its direct impact on the plant, animal and microorganisms, its importance in the global cycle of water, energy and carbon, the energy exchange between air and soil is known for its natural water cycle (especially in the distribution of rain between surface runoff and infiltration) and the management of water and soil resources. Soil moisture plays an important role in the interactive processes between the atmosphere and the earth and global climate change. Triangular and trapezoidal methods combining thermal and visible data are the most commonly used methods for determining the amount of soil surface moisture. The aim of this study is to estimate the surface moisture of the soil (TVDI), by the triangular method in the south of West Azerbaijan province using land temperature index (LST) and vegetation index (NDVI), during 2010, 2014 and 2018. Materials and Methods The present study using MODIS timing series images, NDVI index and LST index, to estimate the surface moisture index (Temperature–Vegetation Dryness Index, TVDI), in three time periods including; the first time period from 1 January 2010 to 30 December 2010 and the second period is from 1 January 2014 to 30 December 2014 and the third period is from 1 January 2018 to 30 December 2018. During each period, 12 images were used on the 15th day of each month. Also, surface moisture was estimated by two methods, one was to establish a high regression relationship and remove the minimum temperature, and the second method was to establish a high and low regression relationship of the pixels. To evaluate the accuracy of these two methods, a regression correlation between the results of these methods with the soil surface moisture content of the Agricultural Jihad (30 points) at a depth of 5 to 15 cm was used. The reason for choosing these three years is due to the difference in high rainfall in some months of the studied years. This study was conducted in the south of the province of West Azerbaijan, which is part of the western region of Iran. Results and Discussion The evaporative triangle diagram consisting of the vegetation index and the surface temperature of the earth in 2010 from January to December month has seen many temperature changes. These same changes in the Earth’s surface temperature have caused that the graphs have many changes. During the 2010 year, according to the chart, the maximum temperature was August and the minimum was January, and the maximum vegetation was May and the minimum was December. In 2014, the maximum temperature in August and the minimum in January and the maximum vegetation in May and the minimum in January and also this year were relatively warmer and drier than in 2010. The evaporative triangle chart in 2018 is rainier than the other two years studied, and the amount of vegetation and according to the graphs in this year, the maximum temperature in July and the minimum is January and the maximum vegetation is May and the minimum in January. The surface moisture level of the soil in 2010 for the western region of Iran, which is the maximum moisture level in May and the minimum in August. In most of the 2010 moisture index maps, the maximum humidity in the west and the lowest in the South of this region. The results of the moisture index maps in 2014 this year have been relatively drier than in other years studied. In 2014 has little rainfall and vegetation. Humidity changes this year are lower than in 2010. The maximum and minimum humidity in 2014 was between 0 and 0.6. The maximum humidity is June and the minimum is August. The TVDI moisture index maps for 2018 have had more moisture indicators this year than in the other two years. In 2018, heavy rains caused the vegetation to increase and the ground temperature to decrease, and this has led to an increase in the moisture index compared to 2010 and 2014. In 2018, the vegetation reached 0.89. But in other years it has been studied up to 0.7. This year, the high humidity is in May and the lowest in August. The maximum humidity during this year is in the west and the lower is in the south. The results of the TVDI index for 2010, 2014 and 2018, using the second method, the general results of this method are similar to the first method. Based on the results obtained from the accuracy of both methods, we conclude that the accuracy of the first method is better and generally simpler than the second method. In 2018, in May, according to the first method, the amount of R2 = 0.67, and also according to the second method, the amount of R2 = 0.41. Conclusion Estimation of surface soil moisture is essential for optimal management of water and soil resources. Surface soil moisture is an important variable in the water cycle of nature, which plays an important role in the global balance of water and energy through its impact on hydrological, ecological and meteorological processes. Examination of the two methods used indicates that the first method, which was also used in this research in general, has higher accuracy than the terrestrial fields due to the results of image accuracy.  In 2010, the months of May and August, according to the first method are R2 = 0.61 and 0.57. In 2010, the amount of R2 according to terrestrial data and the use of the second method in May and August are R2 = 0.43 and 0.47. Also, in 2018, the value of R2 using the first method in May is 0.66. In 2018, the value of R2 using the second method in May is 0.41. The results of the soil surface moisture index, in this study, showed that this model is able to estimate the amount of soil moisture in large geographical areas with acceptable accuracy.