چکیده:
روشهای زیادی برای محاسبه تبخیر-تعرق وجود دارد که به دادههای زیادی نیاز دارد، اما تعدادی از روشهای فقط به دمای هوا نیاز دارند. یکی از این روشها ترونتویت است. منطقه سیستان در جنوب شرقی ایران یکی از مناطقی است که با توجه به بادهای 120 روز در ایران منحصر به فرد است. هدف از این تحقیق: 1) ارزیابی 6 روش مختلف ترونتویت موجود در منابع در مقایسه با روش فائو پنمن مانتیث و 2) اصلاح معادله برای منطقه بادخیز سیستان است. نتایج نشان داد معادله اصلی ترونتوایت مقدار تبخیر-تعرق را کم-برآورد محاسبه میکند. در بین روشهای موجود استفاده از ضریب k=0.72 دارای بهترین نتایج بود. برای واسنجی معادله ترونت-وایت ضریب دمای موثر معادله (k) باید اصلاح گردد. نتایج نشان داد مقدار بهینه k بین 755/0 تا 04/1 متغییر هست. میانگین سالانه مقدار جذر میانگین مربعات خطا (RMSE)، با توجه به مقادیر k متغییر برابر با 14/0 میلیمتر در روز بدست آمد. همچنین با استفاده از حداقلکردن مربعات خطا مقدار k را به طور ثابت 802/0 در نظر گرفتیم که مقدار RMSE برابر با 19/1 میلیمتر در روز بدست آمد. میتوان نتیجهگیری کرد که پس از اصلاح معادله ترونتوایت میتوان آن را در منطقه سیستان استفاده کرد.
IntroductionEvapotranspiration is one of the most important components of the hydrologic cycle in the nature and its exact determination is essential for water balance studies, irrigation and water resources management. One of the most accurate methods of estimating ETo is in different climates of FAO Penman-Monteith equation (PMF-56) the accuracy of this method for calculating evapotranspiration in all over the world has been proven successfully. PMF-56 method as a standard method requires a large meteorological data such as air temperature, relative humidity, wind speed at a height of 2 meters and solar radiation. Provide exact data on all areas is not and also if it is possible it is not very reliable. Some of the evapotranspiration methods require only air temperature for ETo estimation. One of these methods is Thornthwaite. The Sistan region in the southeast of Iran is one of the regions that is unique in Iran due to the 120-day winds and high day-night temperature changes. The purpose of this research is to evaluate the different existing methods of the Thornthwaite equation (6 methods) and adjustment the Thornthwaite equation for the Sistan region.Methodology (Materials and methods or Data and methods)Six Thornthwaite approaches were used in this study are described below:1) The Thornthwaite method (Thornthwaite 1948) is a temperature-based method for the estimation of ET0 as a function of the average monthly temperature.2) Camargo et al. (1999) improved the performance of the Thornthwaite method using an effective temperature (T_eff) instead of the mean temperature (T_ave) and of the daily temperature amplitude:(1) T_eff=k(T_ave+A)=1/2 k(3T_max-T_min)where k = calibration coefficient. Camargo et al. (1999) found that k = 0.72 is the best value for estimating monthly ET0.3) Pereira and Pruitt (2004) recommended k = 0.69 for estimating daily ET0.4) Trajkovic (2005a) expressed the equation reference evapotranspiration where Thornthwaite equation depend on maximum possible duration of sunshine and mean air temperature in the i-th month.5) The data from Serbian stations Palic, Belgrade, and Nis were used to calibrate the Thornthwaite equation (Trajkovic 2005b):(2) 〖ET〗_(C-TW)=0.88〖ET〗_TW+0.5656) Bautista et al. (2009) calibrated the Thornthwaite equation by changing the value of the corresponding constant p = 16. In this research, in addition to evaluating different methods and determining the best method, the value of coefficient k for Sistan region is recalibrated and its best value is stated.Results and Discussion The results based on FAO-56 PM equation show that amount of daily evapotranspiration varied from about 21 mm per day in July to 0.7 mm day-1 in the beginning of December.The annual average value of evapotranspiration according to the FAO-56 PM equation is 8.21 mm day-1.Six Thornthwaite approaches were used in this study were compared to full FAO-56 PM equation. The statistical summary including RMSE (mm day-1), MBE (mm day-1) and r (〖ET〗_eq/〖ET〗_PM ) for Sistan location is presented in Table 1. The value of RMSE in original Thornthwaite equation was equal to 4.15 mm day-1, the value of MBE was equal to -3.92 mm day-1, and the value of r was equal to 0.52, which indicates Thornthwaite equation was very poor in estimating ET0 and greatly underestimated PM values. The method of Camargo et al. (1999) with RMSE value of 3.14 mm day-1 and MBE value of -2.49 mm day-1 and r value of 0.7 had the best results among all methods. Generally, the accuracy of this method is low compared to the PMF-56 method. Among the studied approaches, only the method of Bautista et al (2009) overestimates the amount of evapotranspiration. The average amount of annual evaporation in this approache is equal to 11.84, which has been overestimated by about 44%.Table 1- The value of statistical indicators using different Thornthwaite approaches annual average r MBE (mm day-1) RMSE (mm day-1) Method 4.31 0.52 -3.92 4.15 Thornthwaite (1948) 5.75 0.70 -2.49 3.14 Camargoet al. (1999)5.06 0.61 -3.17 3.85 Pereira and Pruitt (2004)3.41 0.41 -4.80 5.25 Trajkovic (2005)4.34 0.53 -3.84 4.24 Eq.211.84 1.44 3.63 6.68 Bautista et al. (2009)The best monthly values of k (optimal value) was obtained using the trial and error method and minimizing the value. The results show that the value of coefficient k varies from 0.755 in March to 1.04 in October and its average value is 0.825. By using the solver option and based on the lowest value of the root mean square error between Thornthwaite evapotranspiration and the FAO-56 PM method, the value of k coefficient was also obtained, which was equal to 0.802. In the following, by using three k coefficient values, including optimal variable coefficient (k_opt), average variable coefficient (k_0.825) and obtained from the solver option (k_0.802), the evapotranspiration value were calculated and evaluated. Table (2) shows the statistical indices of the evapotranspiration value calculated with the adjustment Thornthwaite equation based on the coefficient k.The results showed that the use of variable k_opt coefficient with RMSE equal to 0.14 mm day-1, MBE -0.04 mm day-1 and r equal to 0.99 had the best results. It can be concluded that using of k_opt and k_0.802 we will reach satisfactory results in calculating evapotranspiration in Sistan region.Table 2- The value of statistical indices of evapotranspiration calculated by adjustment Thornthwaite based on coefficient kannual average r MBE (mm day-1) RMSE (mm day-1) values k 8.23 0.99 -0.04 0.14 k_opt8.69 1.06 0.48 1.59 0.8258.17 0.99 -0.03 1.19 0.802Conclusion The results of this study showed that the 60% of the amount of evapotranspiration in the Sistan region occurs in the four months of June to September when the winds blow for 120 days. The total evapotranspiration in the three months of December, January and February is equal to 8% of the annual evapotranspiration value. The results showed that the six existing methods of estimating evapotranspiration with Thornthwaite method have low accuracy, so its value should be recalibrated.Hence, based on this study, reference evapotranspiration can be easily calculated for the windy region of Sistan with the available metorological data and the calibrated Thornthwaite equation. In developing countries where good quality data are relatively scarce, using such simple methods may be beneficial for the farmers and local water organizations.