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

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

نویسندگان

دانشگاه تهران

چکیده

تعیین توزیع منطقه‌ای دمای بیشینة روزانه نقش مهمی در اجرای منطقه‌ای مدل‌های اقلیم-محصول و به‌خصوص ارزیابی نواحی تحت تنش دمایی پوشش‌های گیاهی دارد. هدف تحقیق حاضر ارزیابی دقت روش های درون یابی در برآورد منطقه ای دمای بیشینة روزانه در استان اصفهان و نیز مقایسه دقت این روش ها در دوسال مختلف با تعداد ایستگاه های متفاوت و تأثیرآن در میزان دقت روش های درون یابی می باشد. در این پژوهش برای درون یابی دمای بیشینة روزانه استان اصفهان ازداده های ایستگاه های سازمان هواشناسی کشور در سال های 1992 و 2007 استفاده شد. سال 1992 با تعداد30 ایستگاه (19 ایستگاه همدیدی و11 ایستگاه اقلیم شناسی) به‌عنوان سال با تعداد ایستگاه کم و سال 2007 با تعداد 54 ایستگاه (31 ایستگاه همدیدی و23 ایستگاه اقلیم شناسی) به عنوان سال با تعداد ایستگاه های تقریباً دو برابر سال 1992 انتخاب شدند.جهت منطقه ای نمودن داده های نقطه ای دمای بیشینة روزانه روش های درون یابی شامل: وزنی عکس فاصله، کریجینگ، کوکریجینگ، کریجینگ-رگرسیون، روش-رگرسیون چندگانه خطی واسپیلاین بکار گرفته شد. ارزیابی نتایج درون یابی و انتخاب بهترین روش به کمک چهار آمارۀ ریشه میانگین مربعات خطا، میانگین خطای مطلق و ضریب همبستگی انجام شد. بر اساس نتایج بدست آمده روش رگرسیون چندگانه خطی برای درون یابی داده های دمای بیشینه روزانة سال 1992 (با ریشه میانگین مربعات خطا از 41/1 تا 03/4، میانگین خطای مطلق از 98/0 تا55/2 و ضریب همبستگی از61/0 تا 95/0) و روش رگرسیون چندگانه خطی برای درون یابی داده های دمای بیشینه روزانة سال 2007 (با ریشه میانگین مربعات خطا از 99/0 تا 98/3، میانگین خطای مطلق از 77/0 تا 92/2 و ضریب همبستگی از 32/0 تا 97/0) با داشتن کمترین خطا به‌عنوان بهترین روش بوده وروش های کریجینگ-رگرسیون، کریجینگ،کوکریجینگ، وزنی عکس فاصله واسپیلاین در مراتب بعدی قرار گرفتند. با توجه به نتایج به‌دست آمده از این مطالعه برای افزایش دقت تخمین پیشنهاد می‌شود نسبت به احداث ایستگاه های جدید به خصوص در مناطق با تنوع ارتفاعی، اقدام به‌عمل آید.

کلیدواژه‌ها


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

Comparative Study of Regional Estimation Methods for Daily Maximum Temperature (A Case Study of the Isfahan Province)

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

  • Ghamar Fadavi
  • Javad Bazarafshan
University of Tehran
چکیده [English]

Introduction: As the statistical time series are in short period and the meteorological station are not distributed well in mountainous area determining of climatic criteria are complex. Therefore, in recent years interpolation methods for establishment of continuous climatic data have been considered. Continuous daily maximum temperature data are a key factor for climate-crop modeling which is fundamental for water resources management, drought, and optimal use from climatic potentials of different regions. The main objective of this study is to evaluate different interpolation methods for estimation of regional maximum temperature in the Isfahan province.
Materials and Methods: Isfahan province has about 937,105 square kilometers, between 30 degree and 43 minutes to 34 degree and 27 minutes North latitude equator line and 49 degree and 36 minutes to 55 degree and 31 minutes east longitude Greenwich. It is located in the center of Iran and it's western part extend to eastern footage of the Zagros mountain range. It should be mentioned that elevation range of meteorological stations are between 845 to 2490 in the study area. This study was done using daily maximum temperature data of 1992 and 2007 years of synoptic and climatology stations of I.R. of Iran meteorological organization (IRIMO). In order to interpolate temperature data, two years including 1992 and 2007 with different number of meteorological stations have been selected the temperature data of thirty meteorological stations (17 synoptic and 13 climatologically stations) for 1992 year and fifty four meteorological stations (31 synoptic and 23 climatologically stations) for 2007 year were used from Isfahan province and neighboring provinces. In order to regionalize the point data of daily maximum temperature, the interpolation methods, including inverse distance weighted (IDW), Kriging, Co-Kriging, Kriging-Regression, multiple regression and Spline were used. Therefore, for this allocated data (24 days for each year and 2 days for each month) were used for different interpolation methods. Using difference measures viz. Root Mean Square Error (RMSE), Mean Bias Error (MBE), Mean Absolute Error (MAE) and Correlation Coefficient (r), the performance and accuracy of each model were tested to select the best method.
Results and Discussion: The assessment of normalizing condition of data was done using Kolmogrov-Smirnov test at ninety five percent (95%) level of significance in Mini Tab software. The results show that distribution of daily maximum temperature data had no significant difference with normal distribution for both years. Weighed inverse distance method used for estimation daily maximum temperature, for this purpose, root mean square error (RMSE) for different status of power (1 to 5) and number of station (5,10,15 and20) was calculated. According to the minimum RMSE, power for 2 and number of station for 15 in 2007 and power for 1 and number of station for 5 in 1992 were obtained as optimum power and number of station. The results also show that in regression equation the correlation coefficient were more than 0.8 for the most of the days. The regression coefficient of elevation (h) and latitude (y) were almost negative for all the month and the regression coefficient of longitude (x) was positive, showing that decreasing temperature with increasing elevation and increasing temperature with increasing longitude. The results revealed that for Kriging method the Gussian model had the best semivariogram and after that spherical and exponential were in the next order, respectively for 2007 year. In the year 1992, spherical and Gussian models had better semivariogram among others. Elevation was the best variable to improve Co-kriging method as auxiliary data. such that The correlation coefficient between temperature and elevation was more than 0.5 for all days. The results also show that for Co-Kriging method the spherical model had the best semivariogram and after that the exponential and Gussian were in the next order, respectively for 2007 year. In the year 1992, the best model of semivariogram was the linear model and after that the spherical and Gussian models had better semivariogram in the next order.
Conclusion: The results revealed that the application of multiple regression method for interpolation produced less errors between observed and estimated maximum temperature in 1992 (RMSE ranges from 1.41 to 4.03, MAE ranges from 0.98 to 2.55, and r ranges from 0.61 to 0.95). For 2007 year, the best estimation was performed by multiple regression and Kriging-Regression (RMSE=ranges from 0.99 to 3.98, MAE ranges from0.77 to 2.92, and r ranges from 0.32 to 0.97). Kriging, Co-Kriging, IDW, Spline methods were also reasonably performed well and could be used as the second order of priority .In addition, with increasing number of stations in 2007 as compared to 1992, the overall accuracy of model performance in estimation of daily maximum temperature have been improved.

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

  • Daily maximum temperature
  • Interpolation
  • Isfahan
1- Anderson S. 2002. An evaluation of spatial interpolation methods on air temperature in Phoenix, AZ. http://www.cobblestoneconcepts.com/ucgis2summer/anderson/anderson.htm
2- Baghebanbashi M.V. and Sepaskhah A.R. 2001. Assessing the geostatistical methods for estimation of annual rainfall in Fars, Boushehr, Kohkilouieh va Bouirahmad provinces. Journal of Shiraz university,25(4B):432.
3- Bazgeer S., Oskuee E.A., Hagigat M., and Darban Astane A.R. 2012. Assessing the Performance of Spatial Interpolation Methods for Mapping Precipitation Data: A Case Study in Fars Province, Iran.Trends in Applied Sciences Research, 7: 467-475
4- Burrough P.A., and Mc Donnell R.A. 1998. Creating Continuous Surfaces from Point Data. In: Principles of Geographic Information Systems, Burrough, P.A. and M.F. Goodchild, R.A. McDonnell, P. Switzer, M. Worboys (Eds.), Oxford University Press, Oxford, UK., pp: 535-550.
5- Chen D., Gong L., Xu C.-Y., and Halldin S. 2007. A high-resolution gridded dataset for monthly temperaturenormals (1971–2000) in Sweden. Geogr. Ann., 89 A (4): 249–261.
6- Chuanyan Z., Zhongren N., and Guodong Ch. 2005 . Methods for modelling of temporal and spatial distribution of air temperature at landscape scale in the southern Qilian mountains, China, Ecological Modelling, 189: 209–220
7- Daly CH .2006. Guidelines for assessing the suitability of spatial climate data sets. International Journal of Climatology, 26:707-721.
8- Goula Bi, T.A., FadikaV., and Soro G.E. 2011. Improved estimation of the mean rainfall and rainfall-runoff modeling to a station with high rainfall (Tabou) in South-Western Cote D'ivoire. J. Applied Sci., 11: 512-519.
9- Hong Y. Nix H.A., Hutchinson M.F., and Booth T. 2005 .Spatial interpolation of monthly mean climate data for China. International J. of climatology, 25:1369-1379.
10- Hassani pak A.A. 1998. Geostatistics. Tehran University, Tehran. (in Persian)
11- Jarvis C.H., and Stuart N. 2001. A comparison among strategies for interpolating maximum and minimum daily air temperatures. Part II: the interaction between number of guiding variables and the type of interpolation method. J. Appl. Meteorol, 40, 1075–1084.
12- Joly D., Brossard T., Cardot H., Cavailhes J., Hilal M. and Wavresky M. 2010. Temperature interpolation based on local information: the example ofFrance. International Journal of Climatology. (wileyonlinelibrary.com) DOI: 10.1002/joc.2220
13- Mehdizadeh M. 2002. Evaluation of geostatistical methods to estimate temperature and rainfall in Ourmieh lake basin. M.Sc. Thesis, University of Tehran. (in Persian with English abstract)
14- Nadi M. 2010. Application of various interpolation techniques of climatic data for determining the most important factors affecting the tree growth at the elevated areas of Chaharbagh, Gorgan. M.Sc.Thesis, University of Tehran.(in Persian with English abstract)
15- Skirvin S.M., Marsh S.E., McClaranw M.P., and Mekoz D.M. 2003. Climatepatial variability and data resolution in asemi-arid watershed, south-eastern Arizona. Journal of Arid Environments, 54: 667–686.
16- Stahl K., Moore R.D., Floyer J.A., Asplin M.G., and McKendry,I.G. 2006. Comparison of pproaches for spatial interpolation of daily airtemperature in a large region with complex topographyand highly variable station density. Agricultural and Forest Meteorology, 139: 224–236
17- Taskinena. Sirvio H., and Vehvilainen B. 2003.Interpolation of daily temperature in finland. Nordic hydrology, 34(5):413-426.
18- Vicente-Serrano S.M., Saz-Sanchez M.A., and Cuadrat J.M. 2003. Comparative analysis of interpolation methods in the middle Ebro Valley(Spain): application to annual precipitation and temperature. Climate Research, 24: 161–180.
19- Wang Y., and Hou SH. 2009. A new interpolation method for Antarctic surface temperature. Progress in Natural Science, 19:1843–1849.
20- Willmott C.J.1982. Some comments on the evaluation of model performance. Bulletin American Meteorological Society, 63(11):1309-1313
21- Yang J.S., Wang Y.Q., and August P.V. 2004. Estimation of Land Surface Temperature Using Spatial Interpolation andSatellite-Derived Surface Emissivity.Journal ofEnvironmental Informatics, 4 (1): 37-44.
CAPTCHA Image