عنوان به زبان ديگر :
Multipurpose Calibration of SWAT Model in Estimating Runoff, Evapotranspiration, and Crop Yield (A Case Study: Mahabad Plain)
پديد آورندگان :
رجا، اميد دانشگاه تهران - پرديس كشاورزي و منابع طبيعي - گروه مهندسي آبياري و آباداني , پارسي نژاد، مسعود دانشگاه تهران - پرديس كشاورزي و منابع طبيعي - گروه مهندسي آبياري و آباداني , تجريشي، مسعود دانشگاه صنعتي شريف - دانشكده مهندسي عمران
كليدواژه :
ارزيابي , جريان سطحي , مصرف واقعي , تركيب كشت , WaPOR
چكيده فارسي :
مدل SWAT به طور گسترده در مطالعات هيدرولوژيك بزرگ مقياس با توجه به كاربري زمين و شيوه هاي مديريت كشاورزي استفاده مي شود. در اغلب مطالعات عملكرد اين مدل براساس داده هاي ميداني رواناب سطحي واسنجي مي گردد و به ديگر مولفه هاي بيلان آب از جمله تبخيروتعرق كمتر پرداخته شده است. هدف از اين مطالعه واسنجي چندمنظوره مدل SWAT با شاخص هاي تكميلي از جمله تبخيروتعرق و عملكرد محصولات غالب در دشت مهاباد در استان آذربايجان غربي بود. واسنجي و اعتبارسنجي مدل براساس داده هاي ايستگاه هيدرومتري گرديعقوب واقع در خروجي محدوده مطالعه و داده هاي تبخيروتعرق و عملكرد محصولات غالب منطقه شامل گندم، جو، ذرت، چغندرقند، يونجه، سيب و انگور انجام شد. دوره هاي واسنجي و اعتبارسنجي به ترتيب 6 و 4 سال بود. با توجه به دقت داده ها در فرآيند واسنجي چندمنظوره، ضرايب 1، 0/9 و 0/8 به ترتيب براي داده هاي هيدرومتري، عملكرد و تبخيروتعرق در تابع هدف نرم افزار SWAT-CUP و روش SUFI2 درنظر گرفته شد. به طور كلي تحليل حساسيت مدل منجر به انتخاب 43 مورد از حساس ترين پارامترها با توجه به حساسيت آنها به دبي رودخانه، تبخيروتعرق و عملكرد گياه شد به طوري كه پتانسيل واحد گرمايي لازم براي رشد گياه و شاخص LAI از مهمترين پارامترها براي شبيه سازي عملكرد و تبخيروتعرق محصولات در فرآيند واسنجي مدل بهدست آمد. نتايج نشان داد مدل SWAT به خوبي توانسته است جريان سطحي را در دوره هاي واسنجي (0/85=NSE؛ 0/87 =R2) و اعتبارسنجي (0/92=NSE؛ 0/89 =R2) شبيه سازي كند. برآورد جريان سطحي داراي عدم قطعيت بالايي به خصوص در جريان كم در دوره واسنجي بود كه منجر به بيش برآورد (%4/5 =PBIAS) جريان سطحي شده است. در ارزيابي عملكرد مدل SWAT در برآورد تبخيروتعرق واقعي، اين مقادير با خروجي مدل سنجش از دور WaPOR مورد مقايسه قرار گرفت. نتايج نشان داد معيار RMSE در دوره واسنجي بين 13/8 تا 66/4 ميليمتر و در دوره اعتبارسنجي بين 23/9 تا 53/2 ميليمتر متغير است. در عين حال، مقادير تبخيروتعرق واقعي برآورد شده مدل SWAT با برآورد نياز آبي (پتانسيل) مدل CROPWAT مورد مقايسه قرار گرفت. در اين بررسي معيار RMSE در دوره واسنجي بين 24/4 (%3/6) تا 49/4 (%6/5) ميليمتر و در دوره اعتبارسنجي بين 21/7 (%6/2) تا 53/4 (%7/2) ميليمتر براي محصولات مختلف به دست آمد. مقادير معيار PBIAS نشان داد برآورد تبخيروتعرق مدل SWAT در هر دو دوره واسنجي و اعتبارسنجي كمتر از برآورد نياز آبي CROPWAT بوده است. به طور كلي براي اكثر محصولات، مقادير تبخيروتعرق حاصل از مدل SWAT و WAPOR به مراتب كمتر از مقادير تبخيروتعرق برآورد شده توسط مدل CROPWAT است كه مي تواند نشان دهنده كم آبياري باشد. برآورد عملكرد محصولات توسط مدل SWAT با مقادير عملكرد محصولات گزارش شده در منطقه مقايسه شد. نتايج حاكي از برآورد نسبي مدل براي عملكرد محصول بود، در عين اين كه مدل SWAT تا حدي عملكرد محصولات را بيش برآورد كرده است. بررسي حاضر نشان داد واسنجي چندمنظوره مدل هيدرولوژيكي بزرگ مقياس SWAT در يك منطقه تحت آبياري علاوه بر استفاده از شاخص متداول جريان رودخانه، با درنظر گرفتن پارامترهاي تكميلي از قبيل تبخيروتعرق و عملكرد محصولات مي تواند منجر به افزايش دقت نتايج مدل شود.
چكيده لاتين :
In most studies, the SWAT model is calibrated based on surface runoff in hydrometric stations, and the evapotranspiration component, which has a great impact on the estimation of balance components, is not considered or less considered in the calibration process. Therefore, the aim of this study was to estimate the yield and evapotranspiration values of major crops in Mahabad plain in West Azerbaijan province, through multivariate calibration of SWAT model. Calibration and validation of the model based on the data of a hydrometric station (Gerdyaghoub station located at the exit of the plain), and crop yield and evapotranspiration of the major crops of the region including wheat, barley, corn, sugar beet, alfalfa, apple, and grape. A six year calibration period followed by a four-year validation period was used. In the multi-purpose calibration process, the coefficients of 1, 0.9, and 0.8, respectively, for hydrometric data, yield, and evapotranspiration were considered in the objective function of SWAT-CUP software and SUFI2 method. In general, the sensitivity analysis of the model led to the selection of 43 of the most sensitive parameters according to their sensitivity to river discharge, crop yield and evapotranspiration. Also, the heat unit potential required for full plant growth and Leaf area index (LAI) were the most important parameters to simulate the yield and evapotranspiration crops in the model calibration process. The results showed that the SWAT model was able to simulate the surface flow in the periods of calibration (NSE = 0.85; R2 = 0.87), and validation (NSE = 0.92; R2 = 0.89), as well. The model has high uncertainty in surface flow simulation, especially over-estimation in low flow during calibration period (PBIAS = 4.5%). In evaluating the performance of the SWAT model in estimating actual evapotranspiration, these values were compared with the output of the WaPOR remote sensing model. The results showed that the RMSE criteria varies between 13.8 to 66.4 mm in the calibration period, and between 23.9 to 53.2 mm in the validation period. At the same time, the estimated actual evapotranspiration values of the SWAT model were compared with the estimated crop evapotranspiration of the CROPWAT model. In this study, RMSE criteria in the calibration period between 24.4 (3.6%) to 49.4 (6.5%) mm, and in the validation period between 21.7 (6.2%) to 4.4 53 (7.2%) mm were obtained for different crops. PBIAS values showed that the estimate of SWAT model evapotranspiration in both calibration and validation periods was lower than the estimate of CROPWAT crop evapotranspiration. In general, for most crops, the evapotranspiration values obtained from the SWAT and WAPOR models are far less than the evapotranspiration values estimated by the CROPWAT model that can indicate deficit-irrigation. Crops yield performance estimates were compared by SWAT model with crop yield performance values reported in the region. The results showed a relative estimate of the model for crop performance, while the SWAT model somewhat over-estimated crop performance. The present study showed that multifunctional calibration of large-scale hydrological model SWAT in an irrigated area in addition to using the usual river flow index, taking into account additional parameters such as evapotranspiration and crop yield can increase the accuracy of model results.