سنجش میزان خودهمبستگی فضایی دمای سطح زمین با کاربری اراضی نمونۀ مطالعه: شهر اصفهان

در سال‌های اخیر، آب‌وهوا و محیط‌زیست جهانی دچار تغییرات شدیدی شده‌است. گسترش صنعت و استفاده از سوخت‌های فسیلی به‌منظور تولید انرژی موردنیاز موجب افزایش گازهای گلخانه‌ای مانند CO و به‌ویژه 2N2O، CH4 و CFC شده‌است. افزایش گازهای جاذب گرما در جوّ زمین به افزایش دمای سطح زمین منجر خواهد شد. بنابراین بررسی دمای سطح شهرها و مناطق اطراف آنها ضروری است (Santamouris, 2013; Yue and Xu, 2013; Macintyre et al., 2018).

مطالعات زیادی در مناطق مختلف دنیا در این زمینه انجام شده‌است، ازجمله:

سوبرینو و همکاران در پژوهشی دریافتند که جزایر حرارتی شهری (UHI)^[1] عاملی مهم در مطالعات جهانی و تغییر و تعادل گرماست و نماینده‌ای برای تغییرات آب‌وهوایی به‌شمار می‌رود. برمبنای این پژوهش، مطالعۀ این پدیده و بررسی سازوکار آن در برنامه‌ریزی‌های شهری اهمیت بسیار زیادی دارد(Sobrino et al., 2004) .

وان و همکاران در پژوهشی به این نتیجه رسیدند که جزایر حرارتی یک متغیر کلیدی در تجزیه و تحلیل و مدل‌سازی بیلان انرژی سطح، رطوبت سطحی و تبخیر و تعرق و مطالعات تغییر اقلیم در مقیاس‌های گوناگون است(Wan et al., 2002).

برمبنای پژوهش آیرونز و همکاران دمای سطح زمین و ضریب گسیل سطحْ شاخص‌های مهمی در مطالعۀ مدل‌های تعادل انرژی در سطح زمین و بررسی فعل‌وانفعالات سطح زمین در مقیاس منطقه‌ای و جهانی‌اند(Irons et al., 2012).

وان در مطالعه‌ای بیان کرد که ضریب گسیل سطحْ پارامتری است که ویژگی‌های حرارتی یک جسم را توصیف می‌کند و از ویژگی‌های ذاتی جسم به‌شمار می‌رود. برمبنای نتایج این پژوهش، دقت تعیین دمای سطح زمین به دقت تعیین این پارامتر وابسته است. اطلاعات به‌دست‌آمده از دمای سطح زمین می‌تواند در زمینۀ مطالعۀ پدیده‌های گوناگون ازجمله تجزیه و تحلیل جزایر حرارتی شهر، تشخیص آنومالی‌های حرارتی^[2] مربوط به زلزله، شناسایی منابع زمین‌گرمایی، پایش خشکسالی و آتش‌سوزی سودمند باشد (Wan, 1999).

براساس پژوهش یان و بو نوع سازه، پوشش سطح زمین، ترافیک، ازدحام جمعیت و فعالیت‌های صنعتی به شکل‌گیری جزایر حرارتی شهری منجر می‌شود. مستقر شدن پوشش زمین با سنگفرش، سازه‌ها، آسفالت و دیگر ساخت‌وسازهای شهری آثار سردشدگی ابعاد طبیعی را از بین می‌برد. جزایر حرارتی با تغییر الگوی بادهای محلی، تقویت رشد ابرها و مه، افزایش تعداد رعدوبرق و ... به‌طور مستقیم و غیرمستقیم بر سلامت و آسایش ساکنان شهری اثر می‌گذارد. آثار نامطلوب این پدیده تخریب محیط‌زیست، افزایش مصرف، افزایش مقدار ازن سطح زمین و حتی افزایش میزان مرگ‌ومیر است(Yan & baur, 2007) .

کشکی و همکاران پژوهشی را با هدف ارزیابی نقش پارامترهای جغرافیایی بر شکل‌گیری دمای سطح زمین در شهر شیراز انجام دادند. چهار مؤلفۀ مهم و معنادار که ارزش ویژۀ بیش از یک داشتند و حدود 86٫75 درصد واریانس متغیرهای اولیه انتخاب شدند. این چهار مؤلفه به ترتیب اهمیت، شامل تابش‌ها (مستقیم، غیرمستقیم و کلی)، جهت شیب و فاصله از رودخانه، پوشش گیاهی و میزان شیب است. در گام بعدی، ارتباط بین مؤلفه‌ها و دمای سطح زمین با استفاده از مدل‌های OLS و GWR برقرار شد. نتایج نشان داد که مدل GWR کارایی بهتری در نمایش توزیع مکانی دمای سطح زمین در شیراز دارد. جزایرحرارتی شهر شیراز بر محدودۀ فرودگاه و مناطق خاکی و بایر پیرامون شهر منطبق است که اغلب مناطقی با کاربری غیرمسکونی‌اند (Kashki et al., 2021).

گوو و همکاران مطالعه‌ای با هدف کشف ویژگی‌های مکانی و زمانی عوامل تأثیرگذار بر LST چندزمانه با داده‌های چندزمانه و روش تجزیه و تحلیل مکانی در منطقۀ دالیان چین انجام دادند. نتایج نشان داد که دمای سطح زمین در این منطقه به‌تدریج با میانگین 0٫55 درجۀ سانتی‌گراد افزایش می‌یابد. همچنین همبستگی قوی و مثبتی میان توزیع مکانی LSTs وجود دارد و درجۀ برازش میان LST و عوامل گوناگون در مدل رگرسیون وزن جغرافیایی (GWR) بهتر از رویکرد حداقل مربعات معمولی است. به‌علاوه، از توزیع مکانی Local-R2 می‌توان دریافت که با گسترش فضای شهری، تأثیر مناطق شهری بر LST افزایش می‌یابد (Guo et al., 2020).

سیکی و یوهانگ در پژوهشی به‌منظور مطالعۀ الگوهای دمای سطح زمین و کاربری اراضی در هنگ‌کنگ و تغییرات فصلی، از تصاویر ماهواره‌ای و داده‌های جغرافیایی و اقلیمی استفاده کردند. نتایج نشان داد که LST تحت تأثیر انواع کاربری اراضی و پوشش گیاهی است. مدل حداقل مربعات (OLS) نشان داد که LST با شاخص NDBI و NDBal رابطۀ مثبت دارد، اما با شاخص NDVI و NDWI، ارتفاع و سرعت باد رابطۀ منفی دارد. همچنین مدل رگرسیون وزن جغرافیایی (GWR) نشان داد که LST با NDWI و ارتفاع رابطۀ منفی دارد. اندازۀ تأثیر شاخص‌ها با فصل متفاوت است(Siqi and Yuhong, 2020) .

تان و همکاران به‌منظور بررسی رابطۀ دمای سطح زمین با انواع کاربری اراضی در منطقۀ دریاچۀ Donating، از تصاویر لندست 7 استفاده کردند. نتایج نشان داد که LST از نوع سطح تأثیر می‌پذیرد. زمین‌های ساحلی، جنگلی و مزارع شالیزار نسبت‌به مناطق آبی، سواحل گلی، مرداب‌ها و جنگل‌های ساحلی مقدار LST بیشتری داشتند. LST زمستانه در منطقه تقریباً 5/3 درجۀ سانتی‌گراد افزایش داشت که ممکن است به کاهش سطح آب دریاچۀ Donating مربوط باشد. همبستگی منفی بین NDMI و LST قوی و پایدار است و بیشترین تأثیر را بر LST دارد(Tan et al., 2020).

در پژوهش کاراکوش، به‌منظور بررسی رابطۀ میان NDVI، LULC و LST در مرکز شهر سیواس و اطراف آن، تصاویر ماهواره‌ای لندست در یک دورۀ 26ساله (2015‑1989) بررسی شد و شدت UHI نیز به‌دست آمد. نتایج نشان داد که طی این دوره، مناطق شهری و زمین‌های کشاورزی افزایش یافته‌اند، در حالی که زمین‌های بایر کاهش یافته‌است. اراضی شهری و اراضی لخت بیشترین مقدار LST را داشته‌اند و دمای سطح ساختمان‌های شهری روند نوسانی داشته‌است، در حالی که دمای مناطق روستایی روند کاهشی داشته‌است. شدت مثبت UHI مشاهده، و همچنین تشکیل جزایر گرمایی شهری مشخص شد (2019 (Karakus,.

در پژوهش زولو و همکاران، با هدف بررسی رابطۀ میان تغییرات دمای سطح زمین (LST) و افزایش مناطق شهری طی یک دورۀ آماری 2001‑2011 در درۀ پو (Po)، از الگوهای مختلف رشد شهری (UGSP) استفاده شد. نتایج نشان داد که یک دامنۀ کاهش LST (در هر شرایطی) با افزایش تدریجی مناطق شهری همراه است. علاوه بر این، الگوهای خوشه‌ای رشد شهری ازنظر آماری رابطۀ معناداری با روز، شب و شرایط روزانه دارند، در حالی که الگوی پراکندۀ رشد شهری فقط با شب شباهت دارد (Zullo et al., 2019).

هی و همکاران در پژوهشی تأثیر عامل توپوگرافی ارتفاع بر LST را در مناطق کوهستانی بررسی کردند. نتایج نشان داد که میان پارامترهای توپوگرافی و دمای سطح ارتباطی با همبستگی زیاد وجود دارد. برای منطقۀ مدنظر در این پژوهش، ارتباط میان LST و ارتفاعْ ارتباط خطی معکوس بود(He et al., 2019) .

تران و همکاران به‌منظور مشخص کردن رابطۀ میان تغییر پوشش زمین و دمای سطح زمین منطقۀ درون‌شهری هانوی (منطقۀ دشتی و مسطح) از لندست 5، 7 و 8 استفاده کردند. نتایج نشان داد که LST به روش غیرخطی انواع LULC بستگی دارد، تجزیه و تحلیل کانون با استفاده از آمار Getis Ord Gi * امکان تجزیه و تحلیل تغییر الگوی LST را در طول زمان فراهم می‌کند، UHI هم از منظر شهری و هم از نوع توسعۀ شهری تأثیر پذیرفته‌است و پیش‌بینی الگوی LST و بررسی اثر UHI می‌تواند با مدل پیشنهادی با استفاده از رگرسیون غیرخطی و سناریوهای تغییر LULC شبیه‌سازی‌شده انجام شود (Tran et al., 2017)

سورش و همکاران برای برآورد دمای سطح زمین در منطقه‌ای در هند از تصاویر سنجندۀ TM ماهوارۀ لندست 5 و سنجندۀ ETM ماهوارۀ لندست 7 سال‌های 1990، 2001 و 2010 استفاده کردند. آنها دریافتند که متوسط دما با روند ثابتی در حال افزایش است(Suresh et al., 2016) .

ژگاناتان و همکاران طی پژوهشی در شهر چین به این نتیجه دست یافتند که شاخص دمای سطح زمین همبستگی مناسبی با متغیرهای اقلیمی نظیر دما، بارش، تعرق و تبخیر در گسترۀ وسیعی از شرایط محیطی دارد. همچنین ارتباط مستقیمی میان دمای سطح زمین و پوشش گیاهی موجود در یک محیط شهری وجود دارد (Jeganathan et al., 2016).

ژو و همکاران در پژوهش دیگری از طریق تحلیل‌های همبستگی و رگرسیون خطی چندگانه، تأثیر پوشش اراضی بر دمای سطح زمین شهر بالتیمور ایالات متحده را بررسی کردند. در این پژوهش، به مثبت بودن ضرایب رگرسیون مربوط به درصد پوشش ساختمان‌ها و پیاده‌روها در  LSTبه معنی تأثیر مثبت بر افزایش دمای سطح زمین و منفی بودن ضرایب پوشش گیاهی و آبی به معنی تأثیر کاهشی بر دمای سطح زمین اشاره شده‌است(Zhou et al., 2011) .

جیانگ و تاین در پژوهشی در شهر چین از طریق شاخص‌های گیاهی، آثار زیان‌بار تغییرات پوشش و کاربری سطح زمین بر دمای سطح زمین را بر مبنای سه تصویر سنجندۀ TM و ETM+ بررسی کردند. برای این کار، شاخص پوشش گیاهی دما (TVX) را از تصاویر به‌دست آوردند. نتایج نشان داد که تغییر کاربری زمین عاملی مهم در افزایش دمای سطح زمین است. همچنین نشان‌دهندۀ درجه‌حرارت زیاد در مناطق با پوشش گیاهی پراکنده و درجه‌حرارت کم در مناطق با پوشش گیاهی متراکم بود (Jiang & Tian, 2010).

 وانگ و همکاران با استفاده از روش شبکۀ عصبیMultilayer perceptron ، دمای سطح زمین را از داده‌های فراطیفی حرارتی به‌دست آوردند. در این روش، ابتدا روی داده‌های فراطیفی حرارتی یک تبدیل PCA زده می‌شود تا ابعاد داده کاهش یابد. سپس خروجی این تبدیل به شبکۀ عصبی وارد می‌شود و خروجی این شبکه دمای سطح زمین است. دمای به‌دست‌آمده با این روش با محصول حرارتی سنجندۀ infrared Atmospheric Sounding Interferometer ارزیابی شد (Wang et al., 2013).

احمدی و داداشی‌رودباری (1396) در پژوهشی با عنوان «شناسایی جزایر حرارتی شهری مبتنی بر رویکرد زیست‌محیطی با استفاده از تصاویر لندست 8 در دورۀ آماری 2013 تا 2015 در شهر اصفهان» دریافتند که به‌دلیل وجود جزایر حرارتی سرد (UCI)، میان مرکز و حومه شیب حرارتی تندی وجود دارد و بزرگ‌ترین جزیرۀ حرارتی شهری در منطقۀ 6 شناسایی شد.

در پژوهش حاضر، شهر اصفهان با توجه به جمعیت زیاد، صنعتی بودن و وجود آثار گردشگری به‌عنوان منطقۀ موردمطالعه انتخاب شد.

روش‌شناسی پژوهش

منطقۀ مطالعه

شهرستان اصفهان کلان‌شهری در مرکز ایران و مرکز استان اصفهان است. این شهرستان با مساحت 550 کیلومترمربع بین عرض جغرافیایی "24٫22 '32 °32 تا "18٫49 '48 °32 و طول جغرافیایی "27٫82 '33 °51 تا "50٫93 '52 °51 قرار دارد. اصفهان از غرب به خمینی‌شهر و نجف‌آباد، از جنوب به کوه صفه، از جنوب شرق به شهر جدید بهارستان در دامنۀ شمالی رشته‌کوه کم‌ارتفاع لاشتر، از شمال به شاهین‌شهر و از شرق به دشت سجزی محصور می‌شود.

دلیل احداث اصفهان وجود رودخانۀ زاینده‌رود است که از کوه‌های زاگرس به‌ویژه زردکوه سرچشمه می‌گیرد. این شهر حدود 1570 متر از سطح دریا ارتفاع دارد. اصفهان 15 منطقه با مساحت 55072 هکتار دارد که بزرگ‌ترین آن منطقۀ 12 با 8223 هکتار و کوچک‌ترین آن منطقۀ 1 با 810 هکتار است. در سرشماری عمومی نفوس و مسکن سال 1395، جمعیت شهر اصفهان 1٫961٫260 نفر بوده‌است. از میان 15 منطقۀ شهر اصفهان، منطقۀ 8 با 239756 نفر جمعیت بیشترین جمعیت و منطقۀ 11 با 58841 نفر جمعیت کمترین جمعیت را دارد.

شکل 1- نقشۀ منطقۀ مطالعه (نویسندگان، 1400)

Figure 1- Map of the study area (Authors, 1400)

روش پژوهش

شش تصویر انتخاب‌شده تصاویر ماهوارۀ لندست 8 از شهر اصفهان در تابستان سال 2019 است. این تصاویر به زمانی با کمترین میزان پوشش ابر مربوط است. به‌منظور آماده‌سازی داده‌ها، تصحیحات اتمسفری و ژئومتریک روی داده‌های ماهواره‌ای انجام، و از داده‌های رقومی ارتفاع استفاده شد. لندست دو سنسور دارد، یکی OLI (مشاهدۀ زمینی) و دیگری TIRS (مشاهدۀ حرارتی). این دو سنسور با هم 11 باند را تشکیل می‌دهند. تصاویر انتخاب‌شده طول موج 11٫30-10٫30 میکرومتر و تفکیک مکانی 100 متر و باند حرارتی 10 دارند. دمای سطح زمین (LST) و NDVI در گوگل‌ارث انجین برای تابستان سال 2019 محاسبه، و تصحیحات رادیومتریکی روی دیتاها اجرا شد. الگوهای پراکنش فضایی سه حالت دارند: خوشه‌ای، پراکنده و تصادفی.

خودهمبستگی فضایی

در دسته‌بندی الگوهای فضایی، چه خوشه‌ای و چه پراکنده و چه تصادفی، می‌توان شباهت‌ها و تفاوت‌های هر جفت از واحدهای ناحیه‌ای مجاور را اندازه گرفت. وقتی که این شباهت‌ها و تفاوت‌ها برای الگوی فضایی مختصر شوند، خودهمبستگی فضایی اتفاق می‌افتد. خودهمبستگی فضایی به معنای این است که صفت‌های بررسی‌شده خودهمبسته‌اند و همبستگی آنها قابل گواهی به نظم جغرافیایی پدیده‌هاست (Clark & Hosking, 1986). خودهمبستگی فضایی پایدار به معنای این است که صفات پدیده‌های جغرافیایی با هم ارتباطی قوی دارند (خواه قوی یا ضعیف). اگر ضرایب ویژگی پدیده‌های جغرافیایی همسایه ارتباط و نظم ظاهری مشخصی نداشته باشد، ارتباط فضایی ضعیف یا الگوی تصادفی دارد. خودهمبستگی فضایی وسیلۀ تحلیلی معتبری برای نحوۀ تحول الگوی فضایی در طول زمان است. استنباط پژوهش به فهم بیشتر چگونگی تحول الگوهای فضایی از گذشته تا کنون کمک می‌کند. برای سنجش همبستگی فضایی، آماره‌هایی وجود دارد که تصدیق می‌دهند با نقاط یا پلی‌گون‌ها (سطوح ناحیه‌ای) کار شده‌است (رهنما و ذبیحی، 1390).

گونه‌های ملاک خودهمبستگی فضایی

الگوهای متفاوتی برای سنجش آماره‌های (خودهمبستگی) فضایی وجود دارد، شاخص موران و ضریب گری و انتخاب محتمل دیگر، آمارۀ عمومی‌‌ G (Lee & Wong, 2001). شاخص موران و گری مشخصه‌های مشترکی دارند، اما مشخصات (خواص) آماری آنها متفاوت است. بیشتر تحلیلگران با شاخص موران موافق‌ترند، زیرا توزیع مشخصاتش خوشایندتر است (Clark & Hosking, 1986). هنوز هر دو روش بر مقایسۀ ارزش‌های واحد ناحیه‌ای همسایگی استوارند. در صورتی که واحدهای ناحیه‌ای همسایگی در طول ناحیه ارزش‌های مشترکی داشته باشند، سرانجام آماره‌ها بر یک خودهمبستگی فضایی قوی دلالت دارند. چنانچه واحدهای ناحیه‌ای همسایگی ارزش‌های بسیار ناهمسانی داشته باشند، آماره‌ها خودهمبستگی فضایی منفی بسیار قوی را نشان می‌دهند. به هر روی، دو مدل روش‌های متفاوتی را برای مقایسۀ ضرایب همسایگی‌ها به‌کار می‌گیرند.

شاخص میانگین نزدیک‌ترین همسایگی (Average Nearest Neighbor)

شاخص میانگین نزدیک‌ترین همسایگی بر اندازه‌گیری فاصلۀ تک‌تک نقطه‌ها تا نزدیک‌ترین همسایۀ آنها مبتنی است و در تعیین همگرایی و واگرایی انواع به‌کار می‌رود. با این شاخص می‌توان فهمید که آیا توزیع نقاط تصادفی است یا خیر و نوع الگوی پراکنش چگونه است. ابزار میانگین نزدیک‌ترین فاصلۀ همسایگی ابتدا فاصلۀ میان نقطۀ مرکزی هر عارضه با نقطۀ مرکزی نزدیک‌ترین همسایگی‌اش را اندازه‌گیری، و سپس میانگین نزدیک‌ترین همسایگی‌ها را محاسبه می‌کند. اگر میانگین فاصلۀ محاسبه‌شده از میانگین توزیع تصادفی فرضی کمتر باشد، آنگاه می‌توان نتیجه گرفت که توزیع پدیدۀ بررسی‌شده در فضا به‌صورت خوشه‌ای است. اگر میانگین فاصلۀ محاسبه‌شده بیش از میانگین توزیع تصادفی فرضی باشد، آنگاه می‌توان نتیجه گرفت که عوارض به‌صورت پراکنده در فضا توزیع شده‌اند (عسگری، 1390).

آمارۀ عمومی G

در شناخت شکل‌های گوناگونی از طبقه‌بندی، الگوهای فضایی کارآمد نیستند. این الگوها گاه نقاط داغ و گاه نقاط سرد تمرکز نامیده می‌شوند. شاخص موران به‌خوبی برای خواص آماری و توصیف همبستگی فضایی جهانی ساخته شده‌است. برای نمونه اگر ارزش‌های بالا نزدیک یکدیگر باشند، شاخص موران و ضریب گری بر برخورد همبستگی فضایی مثبت نسبتاً زیادی دلالت دارند. این طبقه (خوشه) از ارزش بالا ممکن است به‌عنوان نقطۀ تمرکز (داغ) نامیده شود، اما خودهمبستگی فضایی مثبت زیاد نشان داده‌شده با شاخص موران و ضریب گری ممکن است به‌وسیلۀ ارزش‌های پایین مجاور با یکدیگر به‌وجود آمده باشد. این نوع از خوشه می‌تواند نقطۀ سرد نامیده می‌شود. شاخص موران نمی‌تواند این دو نوع خودهمبستگی فضایی را متمایز کند.

شکل 2- مراحل بررسی تغییرات با تصاویر ماهوارۀ لندست

Figure 2- Steps to examine changes with Landsat satellite imagery

یافته‌های پژوهش

دمای سطح زمین تحت تأثیر سازه‌های متفاوتی ازجمله شرایط توپوگرافی، عوامل اقلیمی و جوّی و پوشش گیاهی قرار دارد. با تغییر هریک از این عوامل، دمای سطح زمین نیز تغییر می‌کند. شکل 3 نقشۀ دمای سطح زمین در تابستان سال 2019 را نشان می‌دهد. همان‌طور که در این شکل نیز نشان داده شده‌است، منطقۀ 1 (محله‌های کوشک، خرم، جامی، علیقلی‌آقا، جوزدان الیادران، صائب و ...) و 3 (محله‌های چرخاب، باغ کاران، خواجو، ملک، گلزار، قلعه طبره، احمدآباد، شهشهان، امامزاده اسماعیل، سرتاوه، سرچشمه، سنبلستان و جویباره) و غرب و جنوب غرب منطقۀ 10 (محله‌های عطار نیشابوری، شیخ طوسی، لاله، فجر، پروین، دشتستان، مصلی، سروستان و عسگریه) و شمال منطقۀ 15 و همچنین شمال شرق منطقۀ 13 کمترین مقدار دمای سطح زمین را دارند. در مناطق 2، 4، 5، 6، 7، 12، 14 و 15 نیز دمای سطح زمین بیشترین مقدار را دارد که دلیل اصلی آن، کاهش پوشش گیاهی و شیب و جهت این مناطق است که باعث تابش مستقیم نور خورشید به این مناطق می‌شود. شکل 4 نشان می‌دهد که دمای سطح زمین روند افزایشی داشته‌است.

شکل 3- نقشۀ دمای سطح زمین در تابستان 2019 در شهر اصفهان (نویسندگان، 1400)

Figure 3- Surface temperature map in the summer of 2019 in the city of Isfahan

(Authors, 2021)

شکل 4- نمودار دمای سطح زمین در تابستان 2019 در شهر اصفهان (منبع: نویسندگان، 1400)

Figure 4- Surface temperature chart in the summer of 2019 in the city of Isfahan

(Source: Authors, 2021)

کاربری اراضی

روش استفاده‌شده برای طبقه‌بندی اطلاعاتْ روش نظارت‌شده است. در این روش، برای رده‌بندی یاخته‌ها از نمونه‌های آموزشی استفاده می‌شود. بدین معنا که با تعریف یاخته‌های مشخص از تصویر برای هریک از کلاس‌ها، عمل طبقه‌بندی در قالب کلاس‌های در نظر گرفته‌شده انجام می‌شود. همچنین لازم به توضیح است که از الگوریتم حداکثر مشابهت برای طبقه‌بندی نظارت استفاده شده‌است. در این روش، فرض بر این است که توزیع داده‌های هر طبقه براساس توزیع نرمال در اطراف یاختۀ میانگین آن طبقه قرار گرفته‌است.

بیشترین کاربری در سطح همۀ شهرهای ایران به کاربری مسکونی مربوط است. بیشترین درصد کاربری مسکونی در سطح مناطق شهر اصفهان به منطقۀ یک مربوط است. با توجه به شکل 5 مربوط به سال 2019 جهت گسترش به سمت شرق و جنوب شرق است و بیشترین مساحت کاربری متعلق به کاربری ساخته‌شده است. مساحت این کاربری 293٫63 کیلومترمربع است که 49٫86 درصد از کل مساحت را دربرگرفته‌است. کمترین مساحت مربوط به کاربری آب است که فقط 1٫05 درصد از کل مساحت را دربرگرفته‌است. کاربری باغات به 95٫69 کیلومترمربع (16٫25 درصد) و کاربری پوشش گیاهی به 106٫90 کیلومترمربع رسیده‌است. کاربری فاقد پوشش‌گیاهی 86٫42 کیلومترمربع است که 14٫67 درصد از کل مساحت را دربرگرفته‌است. کاربری باغات در مناطق 9 و 13 (باغ فردوس و باغ بانوان ناژوان)، منطقۀ 15 (در خیابان خاتون‌آبادی و خیابان پی‌وی‌سی) و منطقۀ 12 دیده می‌شود. کاربری پوشش گیاهی در منطقۀ 6 متمرکز شده‌است. کاربری ساخته‌شده در شمال نقشه در منطقۀ 12 گسترش یافته و شهرک صنعتی محمودآباد اصفهان تأسیس شده‌است. در منطقۀ 14، شهرک امام حسین تأسیس شده و گسترش به سمت شمال شرق بوده‌است.

شکل 5- نقشۀ کاربری اراضی شهر اصفهان در سال 2019 (نویسندگان، 1400)

Figure 5- Land use map of Isfahan city in 2019 (Source: Writers, 2021)

جدول 1- مساحت کاربری‌ها به کیلومترمربع سال 2019 (نویسندگان، 1400)

Table 1- Area of uses per square kilometer in 2019 (Authors, 2021)

هریک از طبقات کاربری اراضی تغییرات دمای سطح زمین ازنظر میانگین، حداقل و حداکثر مشخص شده‌است. در جدول 2، نتایج نشان می‌دهد که بیشترین میانگین به پوشش گیاهی با مقدار 38٫25 درجۀ سانتی‌گراد و کمترین میانگین به کاربری آب با مقدار 29٫32 درجۀ سانتی‌گراد مربوط است.

جدول 2- تغییرات دمای سطح زمین در کاربری اراضی ازنظر میانگین، حداکثر و حداقل (نویسندگان، 1400)

Table 2- Ground surface temperature changes in land use in terms of average and maximum and minimum (Authors, 2021)

خودهمبستگی فضایی دمای سطح زمین در منطقۀ مطالعاتی

مقادیر Z-score زیاد است (710٫18)، در صورتی که P-value نزدیک صفر است. نتایج نشان می‌دهد که خوشه‌بندی با سطح اطمینان زیاد وجود دارد و پراکندگی در حالت خوشه‌بندی بالا قرار گرفته‌است که در شکل 6 دیده می‌شود. استان اصفهان از الگوی خوشه‌بندی بالا تبعیت می‌کند. این یافته با نتیجۀ پژوهش اسدی و باعقیده (1398) مطابقت دارد.

شکل 6- خودهمبستگی فضایی دمای سطح زمین (نویسندگان، 1400)

Figure 6- Spatial autocorrelation of surface temperature (Authors, 2021)

تحلیل لکه‌های داغ

در این تحلیل، یک امتیاز z محاسبه می‌شود. این امتیاز نشان می‌دهد که در کدام محله‌ها، داده‌ها با مقادیر زیاد یا کم خوشه‌بندی شده‌اند. به هر عارضه در چهارچوب عوارضی که در همسایگی‌اش قرار دارد، توجه می‌شود و عارضه لکۀ داغ به‌شمار می‌آید که هم خودش و هم همسایگانش ازنظر آماری معنادارند. این نقشه (شکل 7) نشان می‌دهد که در کدام محله‌ها ازنظر آماری معنادار نیستند و چه محله‌هایی لکه‌های داغ و لکه‌های سرد را نمایش می‌دهند. اگر مقادیر z منفی باشند، نشان‌دهندۀ لکه‌های سرد است و مناطقی که z مثبت دارند، لکه‌های داغ را نمایش می‌دهند. برای لکه‌های داغ و سرد، سه سطح اطمینان 99%، 95% و 90% قابل مشاهده است. با توجه به شکل 7، محله‌های واقع در مناطق 1، 7، 3 و 8، محله‌های غرب منطقۀ 10 و محله‌های واقع در شمال مناطق 5 و 6 دارای لکه‌های سرد با سطح اطمینان 99% و محله‌های واقع در شرق مناطق 2 و 9 با سطح اطمینان 90% و محله‌های واقع در جنوب مناطق 7 و 14 با سطح اطمینان 95% مشاهده می‌شود. در صورتی که محله‌های واقع در مناطق 5 و 6 و غرب منطقۀ 4 و به‌صورت لکه‌های داغ در منطقۀ 15 با سطح اطمینان 99% با رنگ قرمز پررنگ، محله‌های واقع در غرب منطقۀ 12 و لکۀ داغ در مرکز منطقۀ 15 با سطح اطمینان 95% و محله‌های واقع در منطقۀ 2 و همچنین در منطقۀ 15 با رنگ قرمز کم‌رنگ با سطح اطمینان 90% قابل مشاهده است. مناطقی که با رنگ زرد کم‌رنگ دیده می‌شوند، ازنظر آماری معنادار نیستند.

شکل 7- نقشۀ خوشه‌بندی لکه‌های داغ و سرد (نویسندگان، 1400)

Figure 7- Clustering of hot and cold spots (Authors, 2021)

تحلیل لکه‌های داغ و سرد

شکل 8 تصویر الف نمایش‌دهندۀ نواحی کاربری اراضی با خوشه‌های سرد 99% است. در این شکل، کاربری ساخته‌شده بیشترین مساحت با مقدار 293 کیلومترمربع و کاربری مناطق فاقد پوشش گیاهی کمترین مساحت با مقدار 333900 کیلومترمربع را شامل می‌شود.

در نموداری که نواحی کاربری اراضی با خوشه‌های سرد 95% را نشان می‌دهد (شکل 8 تصویر ب)، بیشترین مساحت متعلق به کاربری ساخته‌شده با مقدار 17627400 مترمربع و کمترین مساحت متعلق به کاربری پوشش گیاهی با مقدار 110700 مترمربع است.

در نواحی با خوشه‌های سرد 90% (شکل 8 تصویر پ)، بیشترین مساحت متعلق به کاربری ساخته‌شده با مساحت 7771500 مترمربع و کمترین مساحت متعلق به کاربری پوشش گیاهی با مساحت 34200 مترمربع است.

در نموداری که ازنظر آماری معنادار نیست (شکل 8 تصویر ت)، باغات با مساحت 76596300 مترمربع، کاربری پوشش گیاهی با مساحت 45270000 مترمربع، کاربری ساخته‌شده با مساحت 170870400 مترمربع، کاربری آب با مساحت 3038400 مترمربع (کمترین مساحت) و کاربری فاقد پوشش گیاهی با مساحت 79249500 مترمربع نشان داده شده‌است.

در نمودار لکه‌های داغ 90% (شکل 8 تصویر ث) بیشترین مساحت متعلق به کاربری ساخته‌شده با مساحت 10936800 مترمربع و کمترین مساحت متعلق به کاربری آب با مساحت 99000 مترمربع است.

در نمودار لکه‌های داغ 95% (شکل 8 تصویر ج) کاربری آب کمترین سهم را با مساحت 4500 مترمربع و کاربری ساخته‌شده بیشترین سهم را با مساحت 990000 مترمربع دارد.

نمودار لکه‌های داغ 99% (خوشه‌های گرم) (شکل 8 تصویر چ) نشان می‌دهد که کاربری پوشش گیاهی بیشترین سهم را با مساحت 57249000 مترمربع و کاربری آب کمترین سهم را با مساحت 479700 مترمربع دارد.

شکل 8- نمودارهای سهم کاربری اراضی در لکه‌های داغ و سرد (نویسندگان، 1400)

Figure 8- Charts of land use share in hot and cold spots (Source: Authors, 2021)

نتیجه‌گیری

پژوهش حاضر تغییرات الگوی خودهمبستگی فضایی دمای سطح زمین و ارتباط آن را با کاربری اراضی بررسی کرده‌است. به همین دلیل از شاخص لکه‌های داغ و میانگین نزدیک‌ترین همسایگی استفاده شده‌است. برای اعمال روش‌ها، ArcGIS به‌کار رفته‌است.

نتایج حاصل از شاخص میانگین نزدیک‌ترین همسایگی نشان داد که در اصفهان تغییرات همبستگی فضایی دمای سطح زمین با سطح اطمینان بالا وجود دارد و پراکندگی در سطح اطمینان 99% حالت خوشه‌ای قرار گرفته‌است. براساس الگوی لکه‌های داغ، مناطق 5 و 6 که در جنوب اصفهان قرار گرفته‌اند نقش زیادی در شکل‌گیری و ایجاد نواحی جزایر حرارتی گرمایی با خوشه‌های گرم 99% داشته‌اند. در حالی که نواحی خودهمبستگی فضایی منفی (نواحی خوشه‌بندی سرد 99%) در مرکز شهر اصفهان، یعنی مناطق 1 و 3 به‌دلیل وجود کاربری آب (سی‌وسه پل)، فاقد الگوی معناداری به‌لحاظ آماری بوده‌اند. نتایج نشان داد که مناطق دارای پوشش گیاهی و کاربری آب اثر جزیرۀ گرمایی را تضعیف می‌کنند و مناطقی که مناطق ساخته‌شده در آنها گسترش یافته‌است، اثر جزایر گرمایی را سرعت می‌بخشند. بر این مبنا، در برنامه‌ریزی و توسعۀ شهر، باید به فضای سبز شهری توجه بیشتری کرد. مطالعات آمار فضایی می‌تواند الگوی مناسب و جدیدی برای استفادۀ مسئولان، سیاست‌گذاران و برنامه‌ریزان شهری باشد.

[1]. Urban heat islands

[2]. Thermal anomalies