姜姍姍 占車生 王會肖 富可榮 王飛宇

摘要：利用分布式時變增益水文模型（DTVGM）的產流機制改進CLM3.5（Community Land Model version 3.5）的水文過程參數(shù)化，同時嵌入地下水開采對流域水循環(huán)影響的模塊，并考慮經(jīng)濟社會用水，建立能夠描述流域自然-人文過程的大尺度陸地水循環(huán)模型（CLM-DTVGM）。以地下水開采嚴重的海河流域為研究區(qū)，以中科院青藏高原研究所建立的中國區(qū)域高時空分辨率地面氣象要素數(shù)據(jù)集作為大氣強迫驅動，針對海河流域1980年-2010年的水循環(huán)過程進行模擬，分析地下水開采活動對流域水循環(huán)過程的影響。結果表明：地下水開采活動引起海河流域地表徑流和土壤濕度整體有所減少，蒸散發(fā)則呈增加趨勢，且各水循環(huán)要素空間分異較大。

關鍵詞：海河流域；陸面過程；地下水開采；水循環(huán)；模擬研究

中圖分類號：P641 文獻標志碼：A 文章編號：1672-1683（2016）04-0054-06

Abstract：By improving the hydrological parameterization schemes of the land surface model （CLM3.5） using the Distributed Time Variant Gain Model（DTVGM） which was developed independently，embedding a module which reflects the impact of groundwater exploitation on water cycle process，considering domestic water，industrial water and ecological water，a large-scale land water cycle model was developed which could describe the nature and human process of a basin.The Haihe River Basin where groundwater is extracted seriously was chosen as the study area，using high spatial and temporal resolution meteorological element data sets of China established by the Institute of Tibetan Plateau as atmospheric forcing drivers，this paper conducted simulation test for the water cycle process in 1980-2010 in the Haihe River Basin， and analyzed the impacts of groundwater exploitation activities on the water cycle process.The results showed that：groundwater exploitation caused the decrease in runoff and soil moisture in the overall basin；the evapotranspiration in the overall basin increased， and the spatial variation of water cycle elements was large.

Key words：the Haihe River Basin；land surface processes；groundwater exploitation；water cycle；simulation

近年來，人類對水資源的開發(fā)利用日趨廣泛，極大地擾亂了自然水循環(huán)過程。尤其在地表水資源量無法滿足用水需求或地表水污染嚴重的地區(qū)，持續(xù)開采地下水不但會造成地面沉降、濕地減少、海水入侵以及地下水污染，而且會導致陸地水資源的枯竭[1-2]，同時由于更多的水分被抽取至地表，引起局地土壤濕度和地面溫度的變化，進而影響水循環(huán)過程[3-7]。因此，考慮地下水開采活動對流域水循環(huán)影響研究是十分必要的。

已有很多學者就地下水開采活動對水循環(huán)的影響進行了研究探索：束龍倉和Xunhong Chen以內布拉斯加州普拉特河谷為例，采用Visual MODFLOW研究了地下水開采對河流流量衰減的影響，認為地下水位高出河水位的值越大，抽水對河流基流量減少的影響越明顯，對河流滲漏增加的影響微弱[8]；賈仰文等綜合考慮水文氣象、地下水開采、南水北調等因素的影響，建立了海河流域二元水循環(huán)模型，得出不同規(guī)劃水平年水資源管理推薦方案[9-10]；王中根等針對海河流域的地下水超采問題，利用較為成熟的流域地表水模型SWAT與地下水模型MODFLOW構建了海河流域地表水與地下水耦合模型，用于評估地下水開采對流域水資源的影響[11]。這些模型主要基于水量平衡，側重于對降雨徑流和水分收支的模擬與計算，缺乏對能量平衡的考慮，無法表達陸氣界面的水氣能量交換過程[12]。雖然水文模型研究在耦合社會經(jīng)濟水循環(huán)、嵌入人類活動模塊方面有了一定進展，但是為了適應全球氣候變化研究，水文模型還需要從大尺度和高精度模擬方向進行新的突破[13]。陸面過程模式作為水文模型和氣候模型聯(lián)系的橋梁，可以作為大尺度水文模型發(fā)展的基礎。耦合陸面過程模式和水文模型可以保留水文模型的水量平衡模擬方法，尤其是產匯流過程等，而能量平衡、植被作用等采用陸面模式參數(shù)化方案進行描述，以此改善陸面模式對水文過程描述的不足，提高模式模擬預報的精度。

由于受到地下水農業(yè)灌溉、河道取用水等高強度人為活動的影響，海河流域地表產匯流非線性特征明顯，水文過程模擬的不確定性較大。夏軍等[14]結合水文非線性系統(tǒng)方法與分布式流域水文模擬技術自主研發(fā)了一種分布式時變增益水文模型（DTVGM），其特點在于能模擬人類活動影響下的降雨、徑流等水文變量之間非線性關系[15]，對下墊面實時反應能力較強，尤其在海河、黑河、黃河等半干旱-半濕潤流域的模擬效果較好，且計算消耗小[16]。因此，本研究利用DTVGM的產流機制改進CLM3.5的產流過程，同時嵌入地下水開采對流域水循環(huán)影響的模塊，并考慮經(jīng)濟社會用水，構建能夠描述流域自然-人文過程的大尺度陸地水循環(huán)模型（CLM-DTVGM），然后以中國科學院青藏高原研究所基于Princeton 再分析數(shù)據(jù)建立的一套中國區(qū)域高時空分辨率地面氣象要素數(shù)據(jù)集作為大氣強迫驅動，設置兩種模式（無人類活動影響S1和僅考慮開采地下水S2），對海河流域過去30年的地表徑流、蒸散發(fā)和土壤濕度進行模擬，分析地下水開采活動對海河流域水循環(huán)過程的影響。

1 研究區(qū)和數(shù)據(jù)介紹

1.1 研究區(qū)

海河流域（圖1）位于東經(jīng)112°-120°，北緯35°-43°之間，總面積約31.8萬km²。按照流域西北高，東南低的地勢，大致可分高原、山地和平原三種地貌類型。流域西部和北部的高原及山地占流域總面積的60%，流域東部和南部以農田為主，也是主要的人口居住地，占流域總面積的40%。海河流域多年平均降水量為539 mm，自20世紀50年代以來，降水總體上呈現(xiàn)逐步減少的趨勢。流域年平均氣溫為10.8 ℃，年內變化較為平緩。海河流域作為嚴重的資源性缺水地區(qū)，水資源供需矛盾十分突出，用不足全國1.3%的水資源量，承擔著全國11%的耕地面積和10%人口用水任務 [17]。由于流域內地表水資源非常匱乏，主要通過開采地下水進行灌溉，區(qū)內已形成大范圍常年性淺層地下水位降落漏斗，地下水超采是地下水位下降的主導原因[18]。近50年來，海河流域累計超采地下水高達1 900 億m³，目前地下水年開采量約占流域總供水量的2/3[19]。大范圍超采活動已使海河流域地下水處于不可持續(xù)狀態(tài)，嚴重威脅該地區(qū)的可持續(xù)發(fā)展。

1.2 數(shù)據(jù)介紹

已有的可用于陸面過程模式輸入的再分析數(shù)據(jù)產品，在中國區(qū)域上往往存在著系統(tǒng)偏差，如以Princeton 數(shù)據(jù)。因此，本研究采用中國科學院青藏高原研究所開發(fā)的中國區(qū)域的長時間序列、高時空分辨率的陸面模式驅動數(shù)據(jù)，簡稱ITPCAS數(shù)據(jù)集。該套數(shù)據(jù)以國際上現(xiàn)有的Princeton 再分析資料、GEWEX-SRB 輻射資料以及TRMM 降水資料為基礎，利用中國氣象局常規(guī)氣象觀測數(shù)據(jù)進行校正以消除系統(tǒng)偏差，并通過空間降尺度的方式，得到最高時間分辨率為3 h、最高水平空間分辨率為0.1°的再分析數(shù)據(jù)（包含近地面氣溫、近地面氣壓、近地面空氣比濕、近地面全風速、地面向下短波輻射、地面向下長波輻射、地面降水率）[20]，可以滿足國內陸面過程研究的需要。

經(jīng)濟社會數(shù)據(jù)為2000年的數(shù)據(jù)。其中人口、GDP、農田面積數(shù)據(jù)是基于國家統(tǒng)計局城市社會經(jīng)濟調查總隊編纂出版的《中國城市統(tǒng)計年鑒2000》，并由中國科學院資源環(huán)境科學數(shù)據(jù)中心（www.resdc.cn/english/default.asp）處理為分辨率1 km×1 km的柵格數(shù)據(jù)。這組數(shù)據(jù)在本研究中被處理為0.25°×0.25°的分辨率，以適應陸面模式的要求。

工業(yè)、生活和農業(yè)的單位用水數(shù)據(jù)來自《中國水資源公報2000》（表1），并且在流域內保持固定不變。

2 研究方法

2.1 模型介紹

CLM3.5模式是NCAR發(fā)布的新一代陸面過程模式，是在CLM3.0的基礎上對陸面參數(shù)和水文過程加以改進，引進并完善了徑流、地下水、碳循環(huán)和凍土過程，其具體物理過程在文獻中有較為詳盡的描述[21-22]。CLM3.5原有的地表產流采用TOPMODEL模型中的SIMTOP參數(shù)化方案，即根據(jù)蓄滿產流和超滲產流的機制來計算，其中關鍵參數(shù)是計算單元的飽和因子，其依賴表層土壤的不透水面積，計算較為復雜。

本研究采用DTVGM模型的時變增益因子來改進CLM3.5中的產流模型，即考慮降雨徑流的非線性關系，以及產流過程中土壤濕度不同引起的產流量變化，通過時變增益因子簡化了水文循環(huán)系統(tǒng)的輸入輸出之間復雜的非線性關系，達到與一般Volterra泛函級數(shù)相同的模擬效果。

在對CLM3.5產流模塊改進的基礎上，構建人類活動影響概化模型（圖2），并成功嵌入CLM-DTVGM中。為了滿足每個時間步長內的總需水量Dt，人類需要從附近的河流和含水層汲取水源，從河流中和含水層汲取的水量分別記為Qs和Qg，而開采的水資源量主要用于人類生活Dd、工業(yè)生產Di和農業(yè)灌溉Da三個方面。其中，生活和工業(yè)用水主要消耗于蒸發(fā)，而剩余的水量作為廢水（Dg）返回河道；農業(yè)灌溉用水則作為有效降水降落到土壤表面，并繼續(xù)參加隨后的產流等計算過程。

在圖2中的水資源開采部分，從河流汲取的地表水供水量Qs在CLM3.5中主要從每個格點的總徑流（地表徑流與地下徑流之和）中扣除；而從含水層中汲取的地下水供水量Qg是在計算陸地水儲量時扣除，因此有：

在圖2中的水資源利用部分，工業(yè)和生活產生的廢水量Dg視為α（Di+Dd），且被直接從模式格點柱內移除，不再參與格點柱內的計算（α為工業(yè)和生活用水中返回河道的廢水比例），而模式中的蒸發(fā)量相應地增加（1-α）×（Di+Dd），到達地表的有效降水量也因灌溉而增加Da。

基于這種人類活動影響概化模型，模擬設置了兩種情景：無人類活動影響（S1）、考慮開采地下水（S2）。通過比較兩種情景的模擬結果，可以探討人類對水資源開采利用對水循環(huán)過程產生的影響。本次模擬的區(qū)域范圍設為112°-120°E，35°-43°N，時間范圍是1980年-2010年，輸入數(shù)據(jù)空間分辨率為0.1°×0.1°，輸出數(shù)據(jù)的空間分辨率為0.5°×0.5°。其中，需水量采用2000年統(tǒng)計值（表1）。

2.2 參數(shù)率定與模型驗證

DTVGM 模型的參數(shù)見表2，在進行數(shù)值模擬前需根據(jù)實測流量資料進行參數(shù)率定。本文采用海河流域灤縣和觀臺水文站的實測流量資料率定上述參數(shù)。圖3為灤縣站（1997年-2006）和觀臺站（1997年-2006）實測及模擬年徑流過程，可以看出實測和模擬年徑流過程總體上一致，率定的參數(shù)能夠有效地模擬研究區(qū)域的徑流過程。

3 結果分析

3.1 兩種情景下海河流域年平均徑流深時空變化

基于S1和S2兩種情景模式，由CLM-DTVGM對海河流域1980年-2010年徑流進行模擬，得到地下水開采活動影響下的海河流域多年平均徑流深空間分布變化情況）如圖4、圖5）。兩種情景下海河流域多年平均徑流深整體分布趨勢大致相同：從西北至東南逐漸增大，其中平原區(qū)和入?？谔帍搅鬏^大，山區(qū)徑流較??；S2比S1情景下徑流深有所減少，尤其是在東部平原區(qū)變化較為明顯，這與東部平原區(qū)地下水開采較為嚴重有關。

3.2 兩種情景下海河流域年平均蒸散發(fā)時空特征

基于S1和S2兩種情景模式，由CLM-DTVGM對海河流域1980年-2010年蒸散發(fā)進行模擬，得到多年平均蒸散發(fā)空間分布情況（圖6、圖7）可以看出從西北到東南，兩種情景下海河流域多年平均蒸散發(fā)逐漸增大，這與海河流域降雨的空間分布呈現(xiàn)出較好的一致性?？紤]地下水開采活動后，局部地區(qū)尤其在平原區(qū)的取用水高值區(qū)和入?？诘貐^(qū)，蒸發(fā)變異更明顯。

3.3 兩種情景下海河流域土壤濕度時空特征

S1和S2兩種情景下，海河流域1980年—2010年多年平均土壤濕度（地表以下3.43 m）空間分布模擬結果如圖8、圖9所示，可以看出兩種情景下多年平均土壤濕度空間特征基本一致，其中山區(qū)偏低，東北、西部及西南較高，總體表現(xiàn)為東北、西南高于中部地區(qū)，這主要受氣溫和降水空間分布的影響。考慮地下水開采活動后，研究區(qū)局部的土壤濕度有較為明顯的變化，整個流域呈降低現(xiàn)象，在平原取用水高值區(qū)降低程度較為明顯。

4 結論和討論

本研究采用構建的CLM-DTVGM和中國科學院青藏高原所氣候驅動數(shù)據(jù)集，選用無人類活動影響和考慮開采地下水兩種情景，對過去30年海河流域地表徑流、蒸散發(fā)和土壤濕度變化進行模擬分析，得到如下結論：地下水開采導致海河流域地表徑流整體上呈減少趨勢，但在不同空間分布上影響程度不同，其中東部平原區(qū)地表徑流減少較為明顯；地下水開采導致流域蒸散發(fā)整體上呈增加趨勢，但平原區(qū)的取用水高值區(qū)和入海口地區(qū)蒸散發(fā)呈減少趨勢；地下水開采導致海河流域土濕整體上呈下降趨勢，考慮地下水開采活動后，研究區(qū)局部的土壤濕度有較為明顯的變化，其中在研究區(qū)東北部和西部山區(qū)土壤濕度有一定增加，平原取用水高值區(qū)土壤濕度降低。

Zou等[23]用采用CLM3.5陸面參數(shù)化方案的區(qū)域氣候模式RegCM4對海河流域1970年-2000年的陸面過程進行了模擬，得到的多年平均地表徑流約為48.8 mm。本研究不考慮人類活動過程的情景下模擬的多年平均地表徑流約為35 mm，略小于Zou等模擬結果。不過，本研究模擬時期為1980年-2010年，屬于氣候變暖和人類活動加劇時期，海河流域地表徑流已經(jīng)呈現(xiàn)了減少的趨勢[24]，所以本文模擬結果也相對合理的。Zou等模擬的多年平均蒸散發(fā)為534.7 mm，本研究不考慮人類活動過程的情景下模擬的多年平均蒸散發(fā)為440 mm，蒸散發(fā)模擬偏小可能與近年來海河流域地表徑流減少有關，與徑流的模擬結果保持一致，并且2006年研究資料顯示海河流域多年平均蒸散發(fā)為425.4 mm[25]，進一步驗證了模擬結果的可靠性。

實際上，海河流域地表徑流、蒸散發(fā)和土濕變化的原因很復雜，包括降雨量減少，氣溫增加等氣候變化以及水資源過度開發(fā)利用等，在今后的研究中還應設置不同的模擬情景，分析其他要素對水循環(huán)過程的影響。

參考文獻（References）：

[1] Konikow L F，Kendy E.Groundwater depletion：a global problem [J].Hydrogeology Journal，2005（13）：317-320.

[2] 郭孟卓，趙輝.世界地下水資源利用與管理現(xiàn)狀[J].中國水利，2005（3）：59-62.（GUO Meng-zhuo，ZHAO Hui.Groundwater utilization and present management condition in the whole world [J].China Water Resources，2005（3）：59-62.（in Chinese））

[3] Vorosmarty C J，Sahagian D.Anthropogenic disturbance of the terrestrial water cycle [J].Bioscience，2000，50（9）：752-765.

[4] Pielke R A.Influence of the spatial distribution of vegetation and soils on the prediction of cumulus convective rainfall [J].Reviews of Geophysics，2001，39（2）：151-177.

[5] Kanamitsu M，Mo K C.Dynamical effect of land surface processes on summer precipitation over the southwestern United States [J].Journal of Climate，2003，16（3）：496-509.

[6] Betts A K.Understanding hydrometeorology using global models [J].Bulletin of the American Meteorological Society，2004，85（11）：1673-1688.

[7] Haddeland I，Skaugen T，Lettenmaier D P.Anthropogenic impacts on continental surface water fluxes [J].Geophysical Research Letters，2006，33（8）.doi：10.1029/2006GL026047.

[8] 束龍倉，CHEN Xun-hong，地下水開采對河流流量衰減的影響分析[J].水利學報，2003， 2：112-116.（SHU Long-cang，CHEN Xun-hong.Impact of groundwater pumping on stream flow depletion[J].Journal of Hydraulic Engineering，2003（2）：112-116.（in Chinese））

[9] 賈仰文，王浩，甘泓，等.海河流域二元水循環(huán)模型開發(fā)及其應用—II 水資源管理戰(zhàn)略研究應用[J].水科學進展， 2010， 21（1）：9-15.（JIA Yang-wen， WANG Hao， GAN Hong et al.Development and application of dualistic water cycle model in Haihe River Basin：II.Strategic research and application for water resources management[J].Advances in Water Science，2010，21（1）：9-15.（in Chinese））

[10] 賈仰文，王浩，周祖昊，等.海河流域二元水循環(huán)模型開發(fā)及其應用—I 模型開發(fā)與驗證[J].水科學進展， 2010， 21（1）：1-8.（JIA Yang-wen， WANG Hao， ZHOU Zu-hao ，et al.Development and application of dualistic water cycle model in Haihe River Basin：I.Model development and validation[J].Advances in Water Science， 2010， 21（1）：1-8.（in Chinese） ）

[11] 王中根， 朱新軍， 李尉， 等.海河流域地表水與地下水耦合模擬[J].地理科學進展， 2011， 30（11）：1345-1353.（ WANG Zhong-gen， ZHU Xin-jun， LI Wei， et al.A coupled surface-water/groundwater model for Haihe River Basin[J].Process in Geography，2011， 30（11）：1345-1353.（in Chinese））

[12] 蘇鳳閣，郝振純.陸面水文過程研究綜述[J].地球科學進展，2001，16（6）：795-800.（SU Feng-ge，HAO Zhen-chun.Review of land-surface hydrologocal process parameterization[J].Advances in Earth Sciences，2001，16（6）：795-800.（in Chinese））

[13] Arnell N W.A simple water balance model for the simulation of streamflow over a large geophysical domain[J].Hydrol，1999，217，314-355.

[14] 夏軍，王綱勝，呂愛峰，等.分布式時變增益流域水循環(huán)模擬.地理學報，2003，58（5）：789-796.（XIA Jun， WANG Gang-sheng，LV Ai-feng，et al.A research on distributed time variant gain modeling[J].Journal of Geography Sciences，2003，58（5）：789-796.（in Chinese））

[15] LI L，XIA J，XU C，et al.Evaluation of the subjective factors of the GLUE method and comparison with the formal Bayesian method in uncertainty assessment of hydrological models[J].Hydrol，2010，390，210-221.

[16] 夏軍，葉愛中，王綱勝.黃河流域時變增益分布式水文模型（I）-模型的原理與結構[J].武漢大學學報：工學版.2005（6）：10-15.（XIA Jun，YE Ai-zhong，WANG Gang-sheng.A distributed time-variant gain model applied to Yellow River[J].Engineering Journey of Wuhan University，2005（6）：10-15.（in Chinese））

[17] [JP5]陳雷.在海河水利委員會干部大會上的講話.水利部水利信息中心，2009，http：//www.mwr.cn：81/ldxw/cl/zyjh/20090508082850ea5049.aspx.（CHEN Lei.Speech at the meeting of cadres in the Haihe River Commission.Misitry of Water Resources， Ministry of Water Resources Information Center，2009， http：//www.mwr.cn：81/ldxw/cl/zyjh/20090508082850ea5049.aspx.（in Chinese））

[18] 張光輝，費宇紅，李惠娣，等.海河流域平原淺層地下水位持續(xù)下降動因與效應[J].干旱區(qū)資源與環(huán)境，2002，16（2）：32-36.（ ZHANG Guang-hui， FEI Yu-hong， LI Hui-di，et al.Mechanism and consequence of continual decline of shallow ground water level in Haihe River Basin[J].Journal of Arid Land Resources and Environment，2002，16（2）：32-36.（in Chinese））

[19] 水利部.2008 年中國水資源公報.2010，http：//www.mwr.gov.cn/zwzc/hygb/szygb/.（ Ministry of Water Resources， China Water Resources Bulletin in 2008.2010， http：//www.mwr.gov.cn/zwzc/hygb/szygb/.（in Chinese））

[20] 何杰.中國區(qū)域高時空分辨率地面氣象要素數(shù)據(jù)集的建立[D].北京：中國科學院青藏高原研究所，2010.（HE Jie.Development of surface meteorological dataset of China with high temporal and spatial resolution[D].Beijing：Institute of Qinghai Tibet Plateau，Chinese Academy of Sciences，2010.（in Chinese））

[21] Oleson K W， Dai Y， Bonan G， et al.Technical description of the community land model（CLM）[J].Tech.Note NCAR/TN-461+ STR，2004.

[22] Oleson K W， Niu G Y， Yang Z L， et al.CLm3.5 documentation[J].National Center for Atmospheric Research， Boulder， 2007，34.

[23] ZOU J，XIE Z，YU Y，et al.Climatic responses to anthropogenic groundwater exploitation：a case study of the Haihe River Basin，Northern China[J].Climate Dynamics，2014，42（7-8）：2125-2145.

[24] 張建云，王國慶.氣候變化對水文水資源影響研究[M].北京：科學出版社.2007.（ZHANG Jian-yun，WANG Guo-qing.Research on climate change impact on hydrology and water resources[M].Beijing：Science Press，2007（in Chinese））

[25] 中華人民共和國水利部.中華人民共和國發(fā)展和改革委員會.2006.中國水資源及其開發(fā)利用調查評價.（Ministry of Water Resources，People′s Republic of China.Development and Reform Commission，People′s Republic of China .2006.Development and utilization of water resources in China and its investigation and evaluation.（in Chinese））