晉祠泉域地表蒸散發的遙感計算與分析

蘭 榮 輝

(山西省地質環境監測中心，山西 太原 030024)

1 研究區概況

晉祠泉域位于山西省中部、太原盆地西側，面積為2 030 km2，其中裸露可溶巖面積391 km2。晉祠泉域總體上西北高、東南低，海拔標高在760 m～2 100 m之間[2]。區內多年(1956年—2012年)平均降水量為462.4 mm，降雨年內分布極不均勻。由于大氣降水入滲是晉祠泉域巖溶地下水的主要補給源之一，泉域內植被相對發育，地表蒸散蒸騰作用明顯，所以研究晉祠泉域地表蒸散發對于巖溶地下水補給的影響不容忽視。

2 模型原理及數據來源

2.1 模型原理

計算研究區的蒸散發量采用SEBS模型，該模型通過地表能量平衡原理演化而來，地表面能量平衡是指地表面向大氣空間支出熱量和從空間獲取熱量相平衡的情況，實際為入射到地表面的太陽凈輻射能在地表面的轉換和再分配過程。其公示表示為：

Rn=H+LE+F+P

(1)

其中，Rn為凈輻射；H為感熱通量；LE為潛熱通量；F為土壤熱通量；P為光合作用消耗量。由于光合作用消耗量一般不計，因此將式(1)轉化成式(2)來計算潛熱通量LE：

LE=Rn-H-F

(2)

將下載好的遙感數據影像通過ENVI平臺可以直接求得一般地表參數，包括凈輻射、感熱通量、土壤熱通量等，然后根據式(2)通過平臺求得研究區內每個像元的潛熱通量，繼而將其轉化為蒸散發量[1]。

2.2 數據來源

計算所用到遙感數據為NASA官方網站上獲取的MODIS數據，下載2008年—2012年研究區云量較少的數據產品，用MRT再投影工具重新對數據進行投影變換，再提取出計算時所需要的波段。

以文獻[2]里所涉及的泉域范圍為主要依據，利用ARCGIS平臺矢量化作圖來繪制泉域的邊界矢量數據。

高程數據為地理空間數據云網站數據資源中所獲取的90 m分辨率的DEM數字高程數據。

氣象數據由中國地面氣候資料網站上獲取，下載2008年—2012年5年四個氣象站(太原、興縣、原平、離石)的氣象資料，通過對資料進行整理，得出計算時所用到的氣象資料輸入數據集(文本文檔)。

3 蒸散發量計算及結果驗證

3.1 年度蒸散量計算思路

SEBS模型通過預處理過的MODIS數據及氣象數據只能計算出某一天的蒸散發量，由于數據云層的影響導致一年好多天都是無數據的，因此，通過計算出來蒸散發量與氣象站的同期蒸散量做對比來估算其他天無數據的蒸散發量。

假設通過SEBS模型計算出某月第i天的蒸散發量ESi，通過與第i天氣象站觀測量EMi的對比得到一個比例系數：

μ=ESi/EMi。

通過此比例系數與其余時段氣象站的蒸散發量EMj來估算其余時段的蒸散發量，即：

ESj=μ/EMj。

把估算出來的每天蒸散發量進行相加得到本月的蒸散發量，在得到年度中所有月份的蒸散發量后，進行匯總推算出晉祠泉域全年的蒸散發量[2]。

3.2 結果驗證

為了對計算結果進行驗證，需引入氣象站點獲得的實測蒸發值來進行對比，由于氣象站資料是通過蒸發皿測得的，因此需要先引入蒸發皿折算系數將其折算為大面積水面蒸發值，再引入水陸折算系數將其折算為陸地蒸發值，最后將兩者數據統一為一個標準進行對比。

由于太原氣象站離研究區較近，選用太原站的氣象數據來進行驗證。根據以前學者的研究，水陸折算系數冰期取0.6，非冰期取0.62，蒸發皿折算系數0.9[4]。

通過ENVI計算所得到的蒸散發量存在于每個像素中，選定氣象站所在位置可以將數據提取出來。選定2012年每個月無云日的計算量與氣象站數據折算結果進行對比，結果見圖1。由圖1中可以看出兩者結果大致統一，數據方面經過計算相對誤差都在20%以內，因此可以判定本次遙感計算結果是比較穩定的能夠達到精度要求。

4 蒸散發的影響因素

研究區地表蒸散發變化規律：經計算得出研究區2008年—2012年每年的蒸散發量值分別為406.43 mm，379.08 mm，393.45 mm，403.31 mm，429.75 mm。5年間蒸散發量均值為402.4 mm，最大值為2012年的429.75 mm，最小值為2009年的379.08 mm。

圖2是晉祠泉域2008年—2012年5年內月際蒸散發曲線。可以得出，晉祠泉域每年的蒸散發變化趨勢較為一致，都是凍融期較低，汛期較高。由于冬季植被相繼枯萎，淺層地表大部分被冰層覆蓋，蒸發蒸騰作用不明顯，造成地表蒸發量降低，某些地區會出現蒸散量極低的情況。夏季由于蒸散發量的來源降水量充足，植被茂盛，蒸騰作用較為強烈，此時蒸散發量達到年內峰值。

圖3為氣溫月平均值與蒸散發量月平均值的關系曲線，由圖中可以看出兩者存在指數正相關的關系，相關系數R2=0.675，因此可以得出結論，氣溫越高，蒸散發量越高，反之亦然。

為了驗證降雨量跟蒸散發量的關系，選取了研究區內的13個氣象站，用泰森多邊形對每個氣象站進行分區，在得出每個氣象站的控制面積后，通過與研究區總面積的比值決定每個氣象站降雨量所占比例，再進行疊加計算出研究區5年內每月的降雨量。

圖4為蒸散發量與降雨量的對比關系。可以看出晉祠泉域每年的降雨量遠遠大于蒸散發量，而兩者在年內的分布情況也較為相似，汛期降水量大，蒸散發強烈，冰期降水量變小，蒸散發量也降為最低。由圖4中還可以看出降水量一般年內的最大值較蒸散發量年內最大值晚1個～2個月出現，春季鋒面雨帶還沒到達，而氣溫回升快，蒸發量大，農作物用水增加，這就造成了華北平原春季干旱的情況。

圖5為晉祠泉域蒸散發量與高程對比關系圖。由圖5中看出高程在1 500 m以內時蒸散發量是隨著高程的變大而變大的，隨著高程的增加，蒸散發量也變得不穩定，原因可以歸結為隨著高程到達一定高度，氣候條件與植被的蒸騰作用變得不穩定，計算所用到高程較高的點較少，存在一定的誤差，造成了圖5中線性的不規律波動。

5 結語

1)晉祠泉域2008年—2012年地表蒸散發量年度變化較小，年平均值為402.4 mm。

2)晉祠泉域地表蒸散發量的年內規律為：汛期蒸散量最大，冬季凍融期蒸散量最小。

3)蒸散發量和降水量在年內的變化具有正相關關系，而降水量明顯大于蒸散發量，降水量一般年內的最大值較蒸散發量年內最大值晚1個～2個月出現。

4)氣溫越高，蒸散發量越高，反之亦然。高程在1 500 m以內時蒸散發量是隨著高程的變大而變大的，隨著高程的增加，蒸散發量變得不穩定。