تاثیر تغییر کاربری/پوشش زمین بر دمای سطح زمین در منطقه ساحلی بوشهر

پیشینه و هدف افزایش تغییر کاربری و پوشش زمین در اثر توسعه شهرنشینی باعث افزایش دمای شهرها در مقیاس محلی می شود و تا حد زیادی موجب افزایش استرس اکولوژیکی می شود. در حال حاضر، بسیاری از مناطق شهری با تبدیل گسترده کاربری اراضی و ایجاد مناطق گرمایی جدید مواجهه هستند. تکنیک های سنجش از دور به طور قابل توجهی در تشخیص تغییر کاربری/پوشش زمین (lulc) و پیامدهای آن موثر هستند. سنجنده های مختلف ماهواره ای قادر به شناسایی این مناطق تغییر با استفاده از باندهای مرئی و مادون قرمز نزدیک (vnir) و امواج مادون قرمز کوتاه (swir) هستند. علاوه بر الگوریتم‌های طبقه‌بندی مرسوم کاربری/پوشش اراضی، برخی از شاخص‌های طیفی در تشخیص ویژگی‌های اراضی استفاده می‌شوند. شاخص نرمال شده تفاوت پوشش گیاهی (ndvi) را می توان به عنوان کاربردی ترین شاخص طیفی در این سناریو در نظر گرفت. شاخص نرمال شده تفاوت پوشش گیاهی یک فاکتور مهم در فرآیندهای بررسی در پایش دمای سطح زمین (lst) است و در هر مطالعۀ مرتبط با lst استفاده می شود. ndvi به طور مستقیم در تعیین گسیل مندی سطح زمین استفاده می شود و بنابراین یک عامل مهم برای تخمین lst است. همچنین طبقات کاربری/پوشش اراضی lulc با محدودیت‌های آستانه بهینه در محیط‌های فیزیکی مختلف تعیین می شود. به عنوان یک شاخص پوشش گیاهی، ndvi تا حد زیادی به تغییرات زمانی و مکانی بستگی دارد. از این رو، lst نیز با تغییرات زمانی و مکانی، تغییر می نماید. بنابراین، ارزیابی زمانی و مکانی lst و ndvi یک فاکتور مهم در تهیۀ نقشه و پایش دمای سطح زمین، به ویژه در محیط های شهری است. در این تحقیق lst و ndvi را در مردادماه در اراضی ساحلی بوشهر با استفاده از تصویر ماهواره‌ای لندست برای سال‌های 1990، 2005 و 2020 بررسی می شود. نقشه lulc با مقادیر آستانه مناسب ndvi به دست آمد. از اهداف این مطالعه؛ 1) تجزیه و تحلیل تغییرات زمانی الگوی توزیع مکانی lst در منطقه مورد مطالعه، 2) تعیین تغییرات مکانیزمانی رابطه lstndvi برای کل اراضی مورد مطالعه، و 3) بررسی تغییرات مکانیزمانی رابطۀ lstndvi در انواع مختلف کاربری/ پوشش اراضی.مواد و روش ها منطقه مطالعه اراضی شهر بوشهر که در ساحل شمالی خلیج فارس، با ابعاد 20 در 8 کیلومتر با مساحت 1011.5 کیلومترمربع و با متوسط حداقل دما 18.1 درجه سانتی گراد و متوسط حداکثر دمای 33 درجه سانتی گراد، میزان رطوبت نسبی بین 58 - 75 در صد و متوسط بارندگی سالیانه 272 میلی متر است. داده‌های مورد استفاده در این تحقیق شامل؛ داده سنجنده لندست 8 (oli) و سنجنده مادون قرمز حرارتی (tirs) در سال 2020؛ داده etm+ سال 2005، و داده tm در سال 1990 که از مرکز داده های سازمان زمین شناسی ایالات متحده (usgs) دانلود شده است. سنجنده لندست 8 tirs دارای دو باند tir (باندهای 10 و 11) است که در آنها باند 11 دارای عدم قطعیت در کالیبراسیون است. بنابراین، تنها باند 10tir  (رزولوشن 100 متر) برای مطالعه حاضر توصیه شده است. باند 10tir  به اندازه پیکسل 30×30 متر با روش کانولوشن مکعبی توسط usgs تبدیل گردید. داده های landsat 5 tm تنها دارای یک باند مادون قرمز حرارتی tir (باند 6) با وضوح 120 متر است که همچنین توسط usgs به اندازه پیکسل 30×30 متر با روش کانولوشن مکعبی تبدیل گردید. برای داده لندست tm و +etm وضوح مکانی 30 متر باندهای مرئی به مادون قرمز نزدیک (vnir) استفاده شد. روش طبقه بندی حداکثر احتمال برای تعیین صحت نقشه طبقه بندی کاربری/پوشش اراضی lulc بدست آمده از شاخص آستانه ndvi استفاده شد. در این مطالعه جهت استخراج دمای سطح زمین از الگوریتم تک پنجره ای برای بازیابی lst از سنجنده های ماهواره لندست چند زمانی استفاده شد. شاخص تفاوت پوشش گیاهی (ndvi) برای استخراج انواع مختلف کاربری/پوشش اراضی با استفاده از مقادیر آستانه استفاده شد. این مقادیر آستانه با توجه به تفاوت در محیط فیزیکی متفاوت است. محدودیت‌های آستانه ndvi بر روی تصاویر اعمال شد تا انواع مختلف مختلف کاربری/پوشش اراضی استخراج شود.نتایج و بحث مقادیر دقت کلی طبقه‌بندی lulc در سال‌های 1990، 2005 و 2020 به ترتیب 73.6%، 83.9% و 84.6% است. ضرایب کاپا برای طبقات lulc به ترتیب 0.77، 0.80 و 0.84 برای سال های 1990، 2005 و 2020 بود. در مطالعه حاضر میانگین دقت کلی و میانگین ضریب کاپا به ترتیب 80.7 درصد و 0.80 بود. بنابراین، طبقه بندی نقشه طبقه بندی کاربری/پوشش اراضی بر اساس روش آستانه ndvi قابل قبول می باشد. نتایج این تحقیق، افزایش تدریجی (1.4 درجه سانتی گراد در طول سال های 1990 - 2005 و 2 درجه سانتی گراد در طول 2005 -2020) lst در کل دوره مطالعه را نشان داد. میانگین مقدار lst در سه سال مطالعه کمترین (30.86 درجه سانتی گراد) در پوشش گیاهی و بالاترین (49.07 درجه سانتی گراد) در زمین های بایر و مناطق ساخته شده/مسکونی بود. توزیع مکانی ndvi و lst نشان دهنده یک رابطۀ معکوس است. بهترین (0/97) و کمترین همبستگی (0.80)، در حالی که همبستگی متوسط (0.89) مشاهده شد. همبستگی lst - ndvi بر روی سطح پوشش گیاهی منفی قوی (0.80) بود. lst تا حد زیادی توسط ویژگی های کاربری اراضی کنترل می شود.نتیجه گیری مطالعه حاضر به تحلیل رابطه مکانی و زمانی lst و ndvi در اراضی ساحلی بوشهر با استفاده از 3 مجموعه داده لندست برای سال‌های 1990، 2005 و 2020 پرداخت. الگوریتم تک پنجره ای در استخراج lst استفاده شد. به طور کلی، نتایج نشان داد که lst در سال های مطالعه با ndvi رابطه معکوس دارد. وجود پوشش گیاهی از عوامل اصلی منفی بودن زیاد این همبستگی است. رابطه lst - ndvi برای انواع کاربری/پوشش اراضی lulc متفاوت است. در مناطقی با پوشش گیاهی یک رابطۀ رگرسیون منفی قوی (0.80) بین lst و ndvi برقرار است، میانگین lst منطقه مورد مطالعه 3.4 درجه سانتی گراد طی سال های 1990 - 2020 افزایش یافت. تبدیل سایر کاربری اراضی به مناطق ساخته شده و اراضی بایر تاثیر زیادی بر میانگین lst در مناطق شهری دارد. روند تغییرات lst در هر دو منطقه ساخته شده تغییر یافته و بدون تغییر و اراضی بایر افزایشی بود. این مطالعه می تواند به عنوان مرجعی برای کاربری اراضی و برنامه ریزی زیست محیطی در اراضی ساحلی مورد استفاده قرار گیرد.

The effect of land use change/land cover on land surface temperature in the coastal area of Bushehr

Background and Objective Urbanization accelerates the ecological stress by warming the local or global cities for a large extent. Many urban areas are suffering from huge land conversion and resultant new heat zones. Remote sensing techniques are significantly effective in detecting the land use/land cover (LULC) change and its consequences. Several satellite sensors are capable to identify these change zones by using their visible and nearinfrared (VNIR) and shortwave infrared (SWIR) bands. Apart from the conventional LULC classification algorithms, some spectral indices are used in detecting specific land features. Normalized difference vegetation index (NDVI) can be considered the most applied spectral index in this scenario. NDVI is a dominant factor in LST derivation processes and is used invariably in any LSTrelated study. NDVI is directly used in the determination of land surface emissivity and thus is a significant factor for LST estimation. It also determines the LULC categories by its optimum threshold limits in the different physical environments. Being a vegetation index, NDVI depends largely on seasonal variation. Hence, LST is also regulated by the change of seasons. Thus, seasonal evaluation of LST and NDVI is an important task in LST mapping and monitoring, especially in an urban landscape. In this research, LST and NDVI in August in the coastal lands of Bushehr are investigated using Landsat satellite images for the years 1990, 2005 and 2020. The LULC map was obtained with suitable threshold values of NDVI. The objectives of this study are; 1) to analyze the temporal changes of the LST spatial distribution pattern in the study area, 2) to determine the spatialtemporal changes of the LSTNDVI relationship for the whole studied land, and 3) to investigate the spatialtemporal changes of the LST relationship NDVI in different types of land use/cover.Materials and Methods The land study area of Bushehr city, which is on the northern coast of the Persian Gulf, with dimensions of 20 × 8 km2 an area of 1011.5 km2 and with an average minimum temperature of 18.1 oC and an average maximum temperature of 33 oC, relative humidity between 5875% and the average annual rainfall is 272 mm.  The data used in this research include; Landsat 8 (OLI) and Thermal Infrared Sounder (TIRS) data in 2020; 2005 ETM+ data, and 1990 TM data downloaded from the United States Geological Survey (USGS) (https://earth explorer.usgs.gov). The Landsat 8 TIRS instrument has two TIR bands (bands 10 and 11), in which band 11 has calibration uncertainty. Therefore, only TIR band 10 (100 m resolution) is recommended for the present study. The 10 TIR band was converted to a pixel size of 30 × 30 meters by the USGS cubic convolution method. Landsat 5 TM data has only one TIR thermal infrared band (band 6) with 120 m resolution, which was also converted by USGS to 30 × 30 m pixel size by cubic convolution method. For Landsat TM and ETM+ data, the spatial resolution of 30 m visible to nearinfrared (VNIR) bands was used. The maximum likelihood classification method was applied to validate NDVI thresholdbased LULC classification. In this study, the monowindow algorithm was applied to retrieve LST from multitemporal Landsat satellite sensors. NDVI can extract different types of LULC by using the optimum threshold values. These threshold values can differ with respect to the differences in the physical environment. The NDVI threshold limits were applied to the images to extract the different LULC types.Results and Discussion The overall accuracy values of the LULC classification were 73.6%, 83.9%, and 84.6% in 1990, 2005, and 2020, respectively. The kappa coefficients for the LULC classification were 0.77, 0.80, and 0.84 in 1990, 2005, and 2020, respectively. In the present study, the average overall accuracy and average kappa coefficient were 80.7% and 0.80, respectively. Thus, the NDVI threshold methodbased LULC classification was significantly validated. The results of this research showed a gradual rising (1.4 °C during 1990–2005 and 2 °C during 2005–2020) of LST during the whole period of study. The mean LST value for three study years was the lowest (30.86 °C) on green vegetation and the highest (49.07 °C) on bare land and builtup areas. The spatial distribution of NDVI and LST reflects an inverse relationship. The best (0.97) and the least (0.80) correlation, respectively, whereas a moderate (0.89) correlation was noticed. This LSTNDVI correlation was strong negative (0.80) on the vegetation surface. The LST is greatly controlled by landuse characteristics.Conclusion The present study analyzes the spatial, and temporal relationship of LST and NDVI in Bushehr coastal lands using 3 Landsat data sets for 1990, 2005, and 2020. The monowindow algorithm was applied in deriving LST. In general, the results showed that LST is inversely related to NDVI, irrespective of any year. The presence of vegetation is the main responsible factor for high negativity. The LSTNDVI relationship varies for specific LULC types. The green area presents a strong negative (0.80) regression. The mean LST of the study area was increased by 3.4 °C during 19902020. The conversion of other lands into the builtup area and bare land influences a lot on the mean LST of the city. Both the changed and unchanged builtup area and bare land suffer from the increasing trend of LST. This study can be used as a reference for land use and environmental planning on coastal land.