張 清，岳青華，董曉華，楊百銀，魏 沖，喻 丹，張 特

(1.三峽大學(xué) 水利與環(huán)境學(xué)院，湖北 宜昌 443002；2.水資源安全保障湖北省協(xié)同創(chuàng)新中心，湖北 武漢 430072；3.中國電建集團華東勘測設(shè)計研究院有限公司，浙江 杭州，310014；4.水電水利規(guī)劃設(shè)計總院，北京 100120)

1 研究背景

蒸散發(fā)是水文循環(huán)中一個重要的過程，大約占降雨量的60%～70%，在干旱地區(qū)甚至能達到年平均降雨量的90%[1]，在區(qū)域乃至全球尺度的氣候、水文和能量過程上都是主要的調(diào)控要素[2]。蒸散發(fā)也是流域水分散失的主要途徑[3]，控制著地表和地氣交互作用中可利用的水資源量[4]，影響著農(nóng)業(yè)用水效率、區(qū)域生態(tài)功能及其干濕狀態(tài)[5]。降雨量與蒸散發(fā)量之間的差額是可以供人們利用和管理的水資源量，因此準(zhǔn)確地估算蒸散發(fā)量，對于氣候及土地利用變化對水資源的影響研究至關(guān)重要[6]。

同降雨、徑流等相比，通過測量得到實際蒸散發(fā)數(shù)據(jù)的過程是非常困難的[7]。目前，農(nóng)田尺度的蒸散發(fā)測量方法主要有基于能量平衡原理的波文比能量平衡法[8]及渦度協(xié)方差法[9]等；區(qū)域尺度主要分為地面觀測和遙感觀測兩種，地面觀測應(yīng)用最為廣泛的為激光閃爍儀及蒸發(fā)皿測量，但由于成本較高不能做到大面積推廣，只能通過數(shù)量較少的站點觀測甚至無觀測數(shù)據(jù)。以此取得的數(shù)據(jù)缺乏空間代表性，不能有效地考慮到氣候、土壤、土地利用以及坡度等因子的時空差異性對蒸散發(fā)的影響，因此需要在區(qū)域尺度上建立相應(yīng)的蒸散發(fā)估算方法與模型。遙感技術(shù)的發(fā)展豐富了水文模型的數(shù)據(jù)源，可以實現(xiàn)在水文模型中利用氣象、土壤、土地利用等數(shù)據(jù)與徑流數(shù)據(jù)建立多目標(biāo)率定以模擬所需的水文過程，并且已有較大程度的改進[10]。

近年來有較多的學(xué)者利用SWAT模型中幾種潛在蒸散發(fā)估算方法(Penman-Monteith、Hargreaves、Priestley-Taylor等公式)以及增加不同的水文變量進行流域的蒸散發(fā)模擬[11]，也有對不同區(qū)域尺度蒸散發(fā)估算方法差異的對比研究[12]。研究者在考慮其他過程時，可以采用多變量、多站點率定的方式，以此增加模型模擬的多樣性及可靠性[11]。Bai等[13]通過GRACE(gravity recovery and climate experiment)反演陸地水儲量變化與觀測徑流多目標(biāo)率定模型，提高了干旱地區(qū)水儲量變化與蒸散發(fā)的模擬精度。Parajka等[14]在SWAT模型中利用MODIS(moderate-resolution imaging spectroradiometer)降雪數(shù)據(jù)與觀測徑流數(shù)據(jù)建立多目標(biāo)率定，與僅使用徑流率定參數(shù)的情景相比，徑流模擬精度略有降低，但降雪的模擬精度得到了較大提高。Rajib等[15]采用遙感蒸散發(fā)數(shù)據(jù)與觀測徑流數(shù)據(jù)建立多目標(biāo)參數(shù)率定情景，徑流模擬精度無較大變化，蒸散發(fā)的模擬精度有較大的提升。在密西西比河流域，通過遙感數(shù)據(jù)進行徑流以及蒸散發(fā)的多目標(biāo)率定，得到NS系數(shù)分別在0.8左右和0.4左右，整體趨勢模擬地較好，證明遙感數(shù)據(jù)能改善模型對缺乏資料地區(qū)水文過程的模擬[16]。而蒸散發(fā)是一個歷時較長的循環(huán)過程，需要在陸面進行長期的監(jiān)測，故而連續(xù)一致的時間序列數(shù)據(jù)是蒸散發(fā)研究的基礎(chǔ)和前提，作為時空分辨率較高數(shù)據(jù)的典型代表，MOD16A2數(shù)據(jù)為較多研究者所采用。其算法在生成全球ET數(shù)據(jù)產(chǎn)品方面具有良好的性能，可在充分考慮植被生長變化下，提供有關(guān)全球陸地水和能源循環(huán)以及環(huán)境變化的關(guān)鍵信息。賀添等[17]在站點尺度以及流域尺度驗證了MOD16A2數(shù)據(jù)的精度，該數(shù)據(jù)對我國森林、農(nóng)田生態(tài)系統(tǒng)類型以及遼河、海河、黃河和淮河流域的蒸散發(fā)模擬精度較高。楊秀芹等[18]通過收集2001-2010年實測徑流和降水統(tǒng)計數(shù)據(jù)，基于水量平衡方程，得出MOD16蒸散發(fā)產(chǎn)品在淮河流域有較好的適用性與一致性的結(jié)論。

本研究在SWAT模型中利用徑流觀測數(shù)據(jù)與遙感蒸散發(fā)數(shù)據(jù)建立不同參數(shù)率定情景，填補SWAT模型僅采用徑流為變量率定參數(shù)從而模擬區(qū)域內(nèi)水文過程的缺陷，模擬蒸散發(fā)和徑流等水文過程，并對各情景模擬結(jié)果在時間與空間尺度上進行分析。本研究所用方法可推廣到其他水文模型和具有不同下墊面特征的地區(qū)。

2 研究區(qū)域概況

淮河流域地處我國東部，介于長江與黃河之間，位于東經(jīng)115°55′～121°25′，北緯30°55′～36°36′。本研究區(qū)為淮河流域上游息縣以上控制流域，位于我國河南省南部，面積約為10 200 km2，西部、西南部及東北部分布著山區(qū)和丘陵，約占總面積的30%，其他區(qū)域大多為平原，耕地為本研究區(qū)主要的土地利用方式。研究區(qū)域概況見圖1。

圖1 研究區(qū)流域概況

該流域多年年平均降雨量為920 mm，根據(jù)楊秀芹等[18]利用彭曼公式估算研究區(qū)2000-2014年潛在蒸散發(fā)量，得出該時段的平均實際蒸散發(fā)量約為623.4 mm，空間上由東部、西南向西北遞增，但逐年平均值呈下降趨勢。本研究區(qū)屬于半干旱氣候與濕潤氣候的過渡帶，冬春季節(jié)干旱少雨，夏秋季節(jié)悶熱多雨，冷暖旱澇轉(zhuǎn)變急劇。

3 數(shù)據(jù)來源與研究方法

3.1 數(shù)據(jù)來源

本研究的主要數(shù)據(jù)有DEM、土地利用、土壤、水文氣象以及2000-2010年MODIS數(shù)據(jù)集中的遙感蒸散發(fā)MOD16A2逐月數(shù)據(jù)(https://terra.nasa.gov/data/modis-data)。其中，DEM為于中國科學(xué)院計算機網(wǎng)絡(luò)信息中心地理空間數(shù)據(jù)云平臺(http://www.gscloud.cn/)下載的GRID數(shù)據(jù)，分辨率為90 m；土地利用數(shù)據(jù)為于中科院資源環(huán)境數(shù)據(jù)云平臺(http://www.resdc.cn/)下載的GRID數(shù)據(jù)，分辨率為1 km；土壤數(shù)據(jù)采用FAO組織構(gòu)建的Harmonized World Soil Database(http://external-world-soil-database/)下載的GRID數(shù)據(jù)；驅(qū)動因子為來自中國科學(xué)院資源科學(xué)數(shù)據(jù)中心(http://www.resdc.cn/)的GRID數(shù)據(jù)；氣象數(shù)據(jù)為下載自SWAT官網(wǎng)(https://swat.tamu.edu/software/arcswat/)的1980-2015年日降水量、氣溫、太陽輻射、風(fēng)速和相對濕度等DBF格式氣象數(shù)據(jù)；由淮河水利委員會提供的2000-2013年徑流數(shù)據(jù)以及2000-2010年蒸散發(fā)站點(息縣、南灣、石山口)蒸發(fā)皿實測水面蒸發(fā)數(shù)據(jù)；本研究利用MRT工具對MOD16A2遙感蒸發(fā)數(shù)據(jù)進行投影轉(zhuǎn)換，對數(shù)據(jù)圖像進行裁剪，并統(tǒng)計分析蒸散發(fā)數(shù)據(jù)。

3.2 研究方法

SWAT(soil and water assessment tool)模型是一個被廣泛使用且最具代表性的分布式水文模型[19-20]，20世紀(jì)90年代由美國農(nóng)業(yè)部農(nóng)業(yè)研究中心研發(fā)。SWAT模型主要基于SWRRB模型，并吸收了CREAMS、GLEAMS、EPIC、ROTO等模型的優(yōu)點，通過進一步整合，形成區(qū)域尺度乃至子流域尺度的流域水文模型，可用于預(yù)測不同土壤類型、土地利用方式、農(nóng)業(yè)管理措施對流域內(nèi)水文、泥沙、水質(zhì)運移過程的影響。該模型的核心思想在于合理地表征流域的空間變異性，通過將流域內(nèi)具有相同土地利用、土壤以及坡度等性質(zhì)的網(wǎng)格離散為不同的水文相應(yīng)單元(hydrological response units，HRUs)，使其能夠模擬氣候因素和下墊面因素的時空變化以及農(nóng)業(yè)耕作措施對流域水文循環(huán)過程的影響[21]。目前有較多使用SWAT模型的研究主要通過徑流率定模型參數(shù)，探究不同未來氣候變化情景下的水文響應(yīng)[22]，模擬調(diào)整農(nóng)田灌溉量以期維持可持續(xù)生產(chǎn)與水資源的保護[23]，模擬重金屬、營養(yǎng)物質(zhì)、泥沙的運移過程以評估減輕水資源污染的措施。

本研究應(yīng)用SWAT模型對研究區(qū)數(shù)字高程模型(digital elevation model，DEM)圖像劃分子流域(見圖1)，并提取各子流域遙感蒸散發(fā)數(shù)據(jù)MOD16A2，對比研究區(qū)內(nèi)3個蒸散發(fā)觀測站點(南灣、石山口、息縣)的實測蒸散發(fā)值與遙感蒸散發(fā)值，驗證遙感蒸散發(fā)數(shù)據(jù)在本研究區(qū)的適用性。然后建立不同的情景探究多目標(biāo)率定對SWAT模型徑流與蒸散發(fā)時空尺度模擬精度的影響。設(shè)置情景方式如下：

(1)情景1(S1)：僅使用實測徑流數(shù)據(jù)率定模型參數(shù)；

(2)情景2(S2)：僅使用遙感蒸散發(fā)月數(shù)據(jù)率定模型參數(shù)；

(3)情景3(S3)：使用徑流與遙感蒸散發(fā)月數(shù)據(jù)同時率定模型參數(shù)。

基于各情景輸出的模擬結(jié)果，對研究區(qū)多年徑流、面平均蒸散發(fā)與子流域尺度上的蒸散發(fā)過程進行模擬精度的判斷，定性及定量分析多目標(biāo)率定對水文過程模擬的影響。

3.3 模型參數(shù)選擇

在SWAT模型中，使用到DEM以及息縣的經(jīng)緯度位置來提取研究區(qū)流域以及研究區(qū)河網(wǎng)，并使用淮河水利委員會提供的河網(wǎng)shp圖進行校正。導(dǎo)入研究區(qū)內(nèi)雨量站點的坐標(biāo)，設(shè)定合適的面積閾值(約360 km2)，劃分得到39個子流域(見圖1，研究中未涉及子流域編號，故不加以敘述表明)、542個水文響應(yīng)單元(HRUs)。在建立研究區(qū)的SWAT模型后，選擇的研究時段為2000-2010年，其中2000年為預(yù)熱期，2001-2005年作為率定期，2006-2010年作為驗證期，S3在S1、S2用觀測徑流數(shù)據(jù)率定模型的基礎(chǔ)上，增加多年的遙感蒸散發(fā)數(shù)據(jù)作為額外變量。

各情景中均采用徑流曲線數(shù)(CN2)、基流系數(shù)(ALPHA_BF)等水文模型徑流率定參數(shù)并加上最大潛在葉面積指數(shù)(BLAI)、輻射利用效率(BIO_E)等生物物理相關(guān)參數(shù)[15]。其中徑流情景(S1)僅率定徑流；蒸發(fā)情景(S2)僅率定蒸散發(fā)；徑流-蒸發(fā)情景(S3)同時率定徑流與蒸散發(fā)。保持模型其他設(shè)定不變，各情景對相同時間段進行水文過程的模擬。對3種情景率定與驗證期的模擬結(jié)果采用納什效率系數(shù)(NS)、確定性系數(shù)(R2)以及標(biāo)準(zhǔn)化均方根誤差(RSR)等目標(biāo)方程進行精度的判斷。本研究中，各情景的初始參數(shù)范圍一致，率定方法為SWAT-CUP所提供SUFI-2算法[24]，對SWAT-CUP中的參數(shù)采用LH-OAT法進行敏感性分析，以降低模型的異參同效性。

4 結(jié)果與分析

4.1 MOD16A2數(shù)據(jù)的驗證

研究區(qū)內(nèi)僅有3個蒸發(fā)站點(見圖1)，分布于研究區(qū)的中部(南灣)，所在子流域多為水域與森林、中部偏下(石山口)，所在子流域多為水域與耕地、東部出口區(qū)域(息縣)，所在子流域多為耕地。圖2為2001-2010年3個蒸發(fā)站的實測蒸散發(fā)數(shù)據(jù)與相應(yīng)站點所在子流域的遙感MOD16A2蒸散發(fā)數(shù)據(jù)的對比。此時蒸發(fā)站觀測實際蒸散發(fā)折算系數(shù)取自《淮河流域淮河水系實用水文預(yù)報方案》[25]中息縣斷面以上控制流域部分的K值，4、5、10月為0.93，6-9月為0.85，其中未給出的11-3月的K值因冬季蒸散發(fā)能力較小取為0.80。

圖2 2001-2010年各蒸散發(fā)站點實測蒸散發(fā)量與MOD16A2蒸散發(fā)的對比

經(jīng)計算，息縣多年月平均ET值為58.41 mm，較MOD16A2值偏高14%；南灣多年月平均ET值為58.48 mm，相應(yīng)偏高3.36%；石山口多年月平均ET值為49.08 mm，相應(yīng)偏低13.1%。3個站點實測值與遙感MOD16A2蒸散發(fā)值線性擬合R2均在0.83以上，說明遙感蒸散發(fā)值與蒸發(fā)皿實測值顯著相關(guān)，驗證了MOD16A2數(shù)據(jù)在本研究區(qū)的適用性。

由圖2可知，2001年的測量值偏高，可能是由于2001年全國氣象表現(xiàn)出大面積降水偏少，四季均出現(xiàn)了不同程度的干旱，土壤濕度降低而蒸發(fā)皿卻為充分供水條件，且根據(jù)實測降水量數(shù)據(jù)顯示，該年降水量僅為461 mm，較多年平均降水量減少近50%。其他年份的實測蒸散發(fā)數(shù)據(jù)與MOD16A2數(shù)據(jù)的趨勢一致，圖3為蒸散發(fā)站實測數(shù)據(jù)與遙感蒸散發(fā)數(shù)據(jù)之間的相對誤差。

圖3 2001-2010年各蒸散發(fā)站實測數(shù)據(jù)與遙感蒸散發(fā)數(shù)據(jù)的相對誤差

由圖3可知, 蒸散發(fā)站實測數(shù)據(jù)與遙感蒸散發(fā)數(shù)據(jù)的相對誤差范圍基本在-50%～150%，個別點出現(xiàn)了異常偏高的情況，如息縣的2001年2月、石山口的2006年1月、息縣的2008年1月，實際觀測值均遠(yuǎn)小于MOD16觀測值。南灣的相對誤差較其他兩個站點小，多表現(xiàn)為實測蒸散發(fā)值較遙感值偏大，這與類似研究結(jié)果一致[26]。息縣與石山口兩個站點的相對誤差較大，可能與該站點土地利用類型多為水域(水庫)相關(guān)。整體來說，MOD16A2數(shù)據(jù)與蒸散發(fā)站點實測數(shù)據(jù)存在著較好的相關(guān)性，可以用于在空間上彌補研究區(qū)內(nèi)的蒸散發(fā)數(shù)據(jù)。

4.2 敏感性分析

對設(shè)定的不同情景，分別采用息縣水文站2000-2010年的月徑流數(shù)據(jù)以及對應(yīng)年份的MOD16A2遙感蒸散發(fā)數(shù)據(jù)對SWAT模型進行率定和驗證。采用SWAT模型全局敏感性分析LH-OAT法[27]，選擇對模型較為敏感的參數(shù)并增加遙感蒸散發(fā)數(shù)據(jù)參與率定[28]。

在大部分研究中，生物物理參數(shù)會采用默認(rèn)的參數(shù)范圍參與模型水文過程的模擬，這些參數(shù)能反映能量平衡機制、地氣交換以及生物生長等過程，也在很大程度上影響著植被對土壤水分的作用[29]。SWAT模型率定過程中進行全局參數(shù)敏感性分析，表1為模型參數(shù)敏感性排序。

表1 模型參數(shù)敏感性排序

由表1可知，徑流參數(shù)如CN2等依舊為模型最敏感的參數(shù)，SOL_Z、ALPHA_BNK、CANMX、DLAI、BIO_E等土壤及植物冠層、葉面積類相關(guān)參數(shù)較為敏感，相較徑流模擬情景而言，在蒸散發(fā)的模擬過程中的貢獻尤為突出，證明在模型水文過程的模擬中，植被等生物物理參數(shù)是不可忽視的。較為突出的葉面積指數(shù)類的參數(shù)主要與植物的生長狀態(tài)有關(guān)，如植物的氣孔開閉狀態(tài)、葉面積大小以及根系水分等[30]。其次，在參數(shù)的敏感性分析過程中，氣溫、降雨降雪等相關(guān)參數(shù)也表現(xiàn)得較為突出，這與Yan等[31]在淮河流域針對大網(wǎng)格季尺度蒸散發(fā)影響因素的研究一致。

4.3 模擬結(jié)果與分析

本研究區(qū)中，水文站觀測多年月平均徑流量為118.05 m3/s。在設(shè)定的不同情景下，對水文站實測徑流量與模型模擬值使用納什效率系數(shù)(NS)、確定性系數(shù)(R2)與標(biāo)準(zhǔn)化均方根誤差(RSR)指標(biāo)判斷模型對徑流量的模擬精度，結(jié)果見表2。

由表2可知，S1情景模擬多年月平均徑流量為123.08 m3/s，略高于觀測值6.99%。S2、S3情景較前兩個情景對徑流的模擬略差，分別比觀測值高13.60%與11.10%；各模擬情景在月尺度模擬徑流時，各情景率定期NS系數(shù)達到了0.71～0.86，RSR達到了0.38～0.54，驗證期NS系數(shù)達到了0.70～0.80，RSR達到了0.29～0.52?？梢奡1的模擬精度要高于其他兩個情景，且有徑流參與率定的情景精度高于僅用蒸散發(fā)率定的情景(S2)。

表2 各情景模型率定和驗證期徑流量精度

根據(jù)指標(biāo)綜合判斷，模型在模擬徑流時，僅使用蒸散發(fā)對徑流過程進行率定(S2)的效果不佳，但蒸散發(fā)與徑流共同率定(S3)的結(jié)果雖較僅率定徑流(S1)較差，但較S2略有提升，率定期NS系數(shù)從0.71上升到0.80，驗證從0.70上升到0.73。根據(jù)NS、RSR等指標(biāo)判斷，模型均達到了良好的模擬狀態(tài)，證明了SWAT模型在本研究區(qū)的適用性。3種情景的徑流過程的模擬結(jié)果如圖4所示，在徑流過程的模擬上，S2模擬的徑流峰值較S1、S3都有所增加，但也在某些時段存在著較大的偏差，如2004、2005年夏季峰值模擬偏差較大，但經(jīng)過徑流數(shù)據(jù)的補充后，峰值模擬較好；2006、2009年夏季徑流峰值較低，可能的原因是該年份較2001-2010年10年平均降水量分別降低17%、14%。

圖4 各情景月徑流量觀測值與模擬值

在本次研究中選擇了分幅號為h27v05的MOD16A2數(shù)據(jù)資料，采用ArcGIS工具箱從SWAT模型劃分的子流域中提取逐月蒸散發(fā)值，作為蒸散發(fā)觀測數(shù)據(jù)參與SWAT模型參數(shù)率定，與各情景模型模擬值對比。圖5為研究區(qū)MOD16A2與模型模擬的逐月面平均蒸散發(fā)值對比，從圖5中可以看出，在各情景下，模型模擬值與遙感觀測值都是較符合的。在研究區(qū)總流域尺度上，將SWAT模型各情景月蒸散發(fā)模擬結(jié)果與MOD16A2蒸散發(fā)數(shù)據(jù)進行線性擬合，結(jié)果如圖6所示，3個情景模型模擬蒸散發(fā)值與遙感觀測值的相關(guān)系數(shù)逐漸升高，從0.86增加至0.90；在增加了遙感蒸散發(fā)數(shù)據(jù)率定模型參數(shù)后，率定期對峰值的修正效果較為明顯，如2002、2003年夏季，但在部分年份仍出現(xiàn)了峰值模擬較差的情況，如2001、2009年夏季，可能為該年份降水減少，故而蒸散發(fā)量降低，但遙感觀測值則始終偏高。

圖5 各情景月蒸散發(fā)量遙感觀測值與模擬值對比(面平均)

圖6 各情景月蒸散發(fā)模擬值與遙感觀測值擬合

就子流域而言，增加蒸散發(fā)變量率定參數(shù)后，S2與S1情景對比，在研究區(qū)的北部和東部的NS系數(shù)有明顯的提高，圖7為各情景在子流域尺度上蒸散發(fā)的模擬精度(NS系數(shù))對比。

圖7 各情景子流域蒸散發(fā)模擬精度NS系數(shù)對比

由圖7可知，S3與S2情景對比，蒸散發(fā)與徑流共同率定模型參數(shù)的情況下，研究區(qū)的東部與中部的模擬精度明顯上升，NS系數(shù)上升的子流域個數(shù)占總個數(shù)的46%，證明增加蒸散發(fā)率定參數(shù)對模型在空間上模擬蒸散發(fā)過程有明顯的改善作用。S3與S1情景對比，NS系數(shù)上升的子流域個數(shù)占總個數(shù)的64%，表明與僅用徑流單目標(biāo)率定模型參數(shù)的情景相比，利用徑流與蒸散發(fā)多目標(biāo)率定模型參數(shù)能有效地改善模型對徑流與子流域尺度上蒸散發(fā)的模擬精度。

5 結(jié) 論

本研究通過設(shè)置3種不同的情景并對其模擬的徑流與蒸散發(fā)對比，證明在模型模擬水文過程中，植被等生物物理參數(shù)是不可忽視的，增加蒸散發(fā)變量也能有效提高SWAT模型對蒸散發(fā)過程的模擬精度，具體結(jié)論如下。

(1)通過與研究區(qū)內(nèi)3個蒸散發(fā)站點觀測數(shù)據(jù)的驗證，MOD16A2遙感蒸散發(fā)數(shù)據(jù)能有效補充本研究區(qū)缺乏的蒸散發(fā)資料，且能反演到子流域尺度上，對于研究蒸散發(fā)在空間上的變化有較大的幫助。

(2)通過SWAT模型敏感性分析，徑流參數(shù)如CN2等為最敏感的參數(shù)，SOL_Z、ALPHA_BNK、CANMX、DLAI、BIO_E等土壤及植物冠層、葉面積類相關(guān)參數(shù)也較為敏感，相較徑流模擬情景而言，在蒸散發(fā)模擬過程中的貢獻尤為突出。

(3)根據(jù)徑流與蒸散發(fā)變量設(shè)置3種參數(shù)率定情景，與徑流率定參數(shù)相比，蒸散發(fā)率定模型參數(shù)模擬徑流的精度雖有所下降，但大部分子流域蒸散發(fā)的模擬精度均呈上升趨勢；同時采用徑流與蒸散發(fā)率定參數(shù)，模擬徑流的NS系數(shù)有明顯上升，空間尺度蒸散發(fā)模擬精度有進一步提高。