氣候變化和人類活動對基流的影響
——以密云水庫上游潮河流域?yàn)槔?/h1>
2019-01-16 10:48:24王曼玉蔡永茂李春征張志強(qiáng)
中國水土保持科學(xué) 2018年6期
王曼玉，蔡永茂，李春征，張志強(qiáng)?
(1.北京林業(yè)大學(xué)水土保持學(xué)院，100083，北京；2.北京市八達(dá)嶺林場，102112，北京；3.北京市共青林場，101300，北京)
基流作為河川徑流的重要組成部分，是水量平衡分析、水資源優(yōu)化配置和水資源管理的重要內(nèi)容，也是維護(hù)河流生態(tài)健康的重要基礎(chǔ)[1]?；魍ǔＳ傻叵滤推渌舆t的水資源補(bǔ)給河川徑流的部分水量組成，氣候變化和人類活動對基流的影響已成為國內(nèi)外學(xué)者研究的熱點(diǎn)。白樂等[2]發(fā)現(xiàn)人類活動是影響禿尾河流域基流變化的主要因素；王敬哲等[3]發(fā)現(xiàn)呼圖壁河基流主要受溫度影響；郭軍庭等[4]計(jì)算蔡家川流域基流指數(shù)，發(fā)現(xiàn)基流指數(shù)隨降水的增加而減少，且農(nóng)地、灌叢、人工林對基流的影響為負(fù)效應(yīng)；L. Ahiablame等[5]計(jì)算了降水和農(nóng)業(yè)用地影響密蘇里河流域基流變化的彈性系數(shù)，結(jié)果發(fā)現(xiàn)降水的增加導(dǎo)致基流增加，而農(nóng)業(yè)用地增加則會導(dǎo)致基流減少。由于流域基流量不能通過直接測定的方法獲得，目前廣泛應(yīng)用的流域基流分割方法主要包括4大類：直接分割法、水量平衡法、同位素法和時間序列分割法[6]。其中，數(shù)字濾波法操作簡單，可以較好地模擬人工分割過程線，減少分割基流的主觀因素[7-8]。黃國如[9]通過對不同數(shù)字濾波法分割基流的結(jié)果進(jìn)行對比，認(rèn)為Chapman-Maxwell濾波法具有精確方便的優(yōu)點(diǎn)。
潮河是北京市主要飲用水源密云水庫的主要支流之一，近50年來該流域徑流在氣候變化與土地利用變化的影響下呈現(xiàn)顯著的下降趨勢[10]。郭軍庭等[11]和Wang等[12]定量研究了氣候變化和土地利用變化對潮河流域徑流量減少的貢獻(xiàn)率。然而，人類活動和氣候變化對基流的影響還缺乏深入的研究，筆者基于Chapman-Maxwell濾波法分割基流，并采用統(tǒng)計(jì)學(xué)方法，研究潮河流域密云水庫上游人類活動和氣候變化對基流的影響，旨在為認(rèn)識和評價水庫水資源變化、制定科學(xué)合理的水資源和流域管理方案提供科學(xué)依據(jù)。
1 研究區(qū)概況
潮河流域位于北京市與河北省的東部，E 116°7′～117°35′，N 40°19′～41°38′，發(fā)源于河北省豐寧縣，匯入北京市密云水庫。密云水庫上游流域總面積4 855.9 km2(圖1)。研究區(qū)氣候類型屬于溫帶季風(fēng)大陸性半濕潤半干旱氣候，四季分明，干旱和冷暖變化明顯，多年平均氣溫為9.16 ℃，年平均降水量490 mm，汛期6—9月之間，汛期降水量占總降水量75%以上，且多以暴雨形式出現(xiàn)。
圖1 潮河流域地形圖及雨量站、氣象站、水文站分布圖Fig.1 Terrain, precipitation station, meteorological station and hydrological station in Chaohe watershed
2 數(shù)據(jù)來源及研究方法
2.1 數(shù)據(jù)收集整理
筆者選取密云水庫上游潮河流域內(nèi)的4個雨量站、1個標(biāo)準(zhǔn)氣象站、1個水文站以及該流域周邊2個氣象站(圖1)的實(shí)測氣象數(shù)據(jù)、降水?dāng)?shù)據(jù)以及日徑流數(shù)據(jù)進(jìn)行分析。這些數(shù)據(jù)由河北水文總站、北京市水文總站和海河流域水利委員會提供。流域降水量采用4個雨量站實(shí)測資料由泰森多邊形法計(jì)算得到，潛在蒸散發(fā)量(potential evapotranspiration，PET)采用Penman-Monteith法計(jì)算得到。由于潛在蒸發(fā)散空間變異不大，其值為根據(jù)3個氣象站資料分別計(jì)算得到的PET的均值。流域土地利用分別解譯自1987年、1999年、2009年3期分辨率為30 m的Landsat TM影像遙感數(shù)據(jù)。將經(jīng)過地理坐標(biāo)配準(zhǔn)和幾何精度矯正的遙感數(shù)據(jù)與野外實(shí)測調(diào)查相結(jié)合，經(jīng)過人機(jī)交互判讀和監(jiān)督分類將土地利用類型分為有林地、灌木林地、草地、耕地、水域、居民區(qū)與建筑用地和未利用地等7類[13]。
2.2 研究方法
2.2.1 基流分割Chapman-Maxwell濾波法 濾波法是分割高頻信號和低頻信號的有效工具，流域快速響應(yīng)的直接徑流與高頻信號相似，慢速響應(yīng)的基流與低頻信號相似，該法被廣泛應(yīng)用于基流分割中[14]。其濾波方程為
(1)
式中：Qi為i時刻的總徑流，mm；Qb,i為i時刻的基流，mm；k為退水系數(shù)，一般情況下k取0.95時，基流分割結(jié)果最為合理[9，15]。
2.2.2 趨勢和突變點(diǎn)的Mann-Kendall檢驗(yàn)
1)趨勢檢驗(yàn)。
Mann-Kendall(以下簡稱MK)檢驗(yàn)被廣泛用于檢驗(yàn)水文、氣象等時間序列的趨勢[5,10]，MK檢驗(yàn)公式如下：
(2)
式中：e為樣本個數(shù)，Xn、Xm是連續(xù)時間序列上的值。當(dāng)Xn-Xm>0時，sgn=1; 當(dāng)Xn-Xm=0時，sgn=0；當(dāng)Xn-Xm<0時，sgn=-1。假設(shè)時間序列沒有變化，則S均值為0，計(jì)算方差V(S)的公式為
(3)
MK統(tǒng)計(jì)量公式如下：
(4)
通常用標(biāo)準(zhǔn)Z檢驗(yàn)表示數(shù)據(jù)變化的趨勢，對于給定的置信水平α，若|Z|≥Z1-α/2,即時間序列在置信水平α上存在明顯的上升或下降趨勢。當(dāng)Z值為正時，表示該時間序列上的值為增加趨勢，反之則為減小趨勢。
2)突變點(diǎn)分析。
設(shè)序列值為X1，X2…Xk，St表示第n個樣本Xn>Xm(1(5)
在原時間序列獨(dú)立的假定下，St的均值和方差分別為：
(6)
(7)
將St標(biāo)準(zhǔn)化：
(8)
2.2.3 雙累積曲線法 雙累積曲線法可以建立降水-徑流間的年累積值回歸關(guān)系，因此被廣泛應(yīng)用于計(jì)算降水與人類活動對年徑流量變化的相對貢獻(xiàn)率，筆者采用該方法分析將降水與人類活動對基流的影響程度[16]。步驟如下：
1)建立基準(zhǔn)期年累積降水量與基準(zhǔn)期年累積基流量的回歸方程
∑Qb=a∑P+b。
(9)
式中：Qb為基準(zhǔn)期基流量，mm；P為基準(zhǔn)期降水量，mm；a為回歸方程的相關(guān)系數(shù)，b為回歸方程截距。
2)根據(jù)基準(zhǔn)期所得回歸方程，帶入研究期累積降水量求研究期累積基流值，再反推研究期的年基流值，則
r=Qb12-Qb1，
(10)
l=Qb2-Qb12。
(11)
式中：r為降水引起基流變化部分，mm；Qb12為利用基準(zhǔn)期所得到的回歸方程推算得到的研究期計(jì)算年基流值，mm；Qb1為基準(zhǔn)期分離基流，mm；l為人類活動引起基流變化部分，mm；Qb2為研究期分離基流，mm。
3 結(jié)果與分析
3.1 降水、PET、徑流和基流變化趨勢和突變點(diǎn)
潮河流域1963—2015年降水、PET、徑流和基流呈現(xiàn)不同的變化趨勢(表1)。流域多年平均降水量為486.8 mm，以平均每年0.5 mm的幅度下降，但這一趨勢并不顯著(P=0.24)(圖2a)，并于2008年發(fā)生突變(圖3a)。PET以平均每年0.9 mm的速度顯著上升(P=0.02)(圖2b)，年平均PET為1 185.8 mm，突變點(diǎn)發(fā)生在2007年(圖3b)。流域徑流以平均每年0.6 mm的速度極顯著下降(P<0.001)(圖2c)，多年平均徑流量為49.5 mm，突變點(diǎn)位于1999年(圖3c)。流域基流呈極顯著下降趨勢(P<0.001)(圖2d)，平均基流量為10.9 mm，平均每年減少0.1 mm，且于1999年發(fā)生突變(圖3d)。由此可見，潮河流域基流量與徑流量的變化趨勢相同，但徑流量變化幅度大于基流量的變化幅度。降水和PET均未在1999年發(fā)生突變，說明影響基流量突變的因素除了氣候因素外還有人類活動的影響。
表1潮河流域年降水、PET、徑流和基流量的MK趨勢檢驗(yàn)及突變點(diǎn)分析
Tab.1MK trend analysis and change point of annual precipitation, PET, streamflow and baseflow in Chaohe watershed mm
圖2 潮河流域年降水量、年P(guān)ET、年徑流量、年基流量線性變化趨勢Fig.2 Linear trend of annual precipitation, annual PET, annual streamflow and annual baseflow in Chaohe watershed
UFt為序列值；為反序列值；a為0.05顯著水平臨界值。UFt: Sequence value; Inverse sequence value; a significant horizontal threshold of 0.05. 圖3 潮河流域年降水、年P(guān)ET、年徑流量和年基流量突變點(diǎn)分析結(jié)果Fig.3 Change points of annual precipitation, annual PET, annual streamflow and annual baseflow in Chaohe watershed
3.2 降水和PET對基流的影響
降水和PET是氣候變化中影響基流的2個主要因素。潮河流域基流與降水呈顯著正相關(guān)關(guān)系(P<0.01)，相關(guān)系數(shù)為0.64(圖4a)；基流與PET呈顯著負(fù)相關(guān)關(guān)系(P<0.01)，相關(guān)系數(shù)為0.27(圖4b)。1963—2015年，降水量的下降和PET的上升共同導(dǎo)致了基流量的減少，但降水對基流的影響比PET對基流的影響更大?；魍蛔凕c(diǎn)前后，年平均降水量減少16.4%(84.4 mm),年平均PET增加7.3%(84.4 mm)，年平均基流減少66.8%(9.2 mm)。從上述結(jié)果可知，PET在突變點(diǎn)前后的變化微小，可以忽略不計(jì)。從水文循環(huán)的角度看，基流的顯著變化不僅是由于氣候變化引起的，人類活動在其中起到了重要作用。
圖4 潮河流域降水-基流、PET-基流關(guān)系圖Fig.4 Precipitation- baseflow and PET- baseflow in Chaohe watershed
3.3 降水和人類活動對基流的影響
以突變點(diǎn)1999年為界，將潮河流域年降水與基流系列分為變化前的1963—1998年(基準(zhǔn)期)和變化后的1999—2015年(變化期)進(jìn)行對比分析?；鶞?zhǔn)期年平均基流量為13.8 mm，變化期年平均基流量為4.6 mm，突變點(diǎn)前后基流量下降了66.8%(表2)。建立1963—2015年累積基流量與累積降水量回歸關(guān)系(P<0.001)(圖5)，1963—1998年2者相關(guān)性較高(R2=0.991 8)，因此可認(rèn)為1963—1998年基流變化是受人類影響較小的時期。將此時期作為基準(zhǔn)期，用來還原1999—2015年受人類影響較小時期的情況，若基流只受降水影響，雙累積曲線中的趨勢應(yīng)與基準(zhǔn)期的走勢保持一致，然而曲線于1999年發(fā)生向下偏轉(zhuǎn)，說明人類活動在很大程度上導(dǎo)致基流量減少。在不考慮其他氣候因素變化的情況下，人類活動對基流在突變點(diǎn)前后變化的最大可能貢獻(xiàn)率高達(dá)81.10%, 降水導(dǎo)致基流變化的貢獻(xiàn)率為18.90%(表2)?？梢?，人類活動是潮河流域突變點(diǎn)后基流減少的主導(dǎo)因素。
表2 降水和人類活動對潮河流域基流的影響Tab.2 Impacts of precipitation and human activities on baseflow in Chaohe watershed
圖5 潮河流域降水量-基流量雙累積曲線Fig.5 Double accumulative curve of precipitation-baseflow in Chaohe watershed
4 討論
4.1 氣候因素對流域基流的影響
在全球變暖的總體趨勢下，潮河流域PET呈現(xiàn)顯著上升趨勢，降水量呈現(xiàn)不顯著下降趨勢，從而導(dǎo)致徑流量和基流量呈現(xiàn)極顯著下降的趨勢。盡管PET上升的趨勢較降水下降的趨勢更為顯著，但降水對基流的影響更為劇烈，這與D. L. Ficklin等[17]的結(jié)論一致。王云琦等[18]也指出流域內(nèi)降水量對產(chǎn)流量的影響大于蒸散量對產(chǎn)流量的影響。潮河流域基流量隨降水量的減小而減小，由于基流具有穩(wěn)定性和滯后性，在長時間尺度下，降水會通過影響產(chǎn)流下墊面特征，進(jìn)而對基流產(chǎn)生影響[19]。1999年以來，北京遭遇連續(xù)7年的干旱，降水量下降明顯，降水落地后先與土壤作用，待土壤飽和后轉(zhuǎn)換成徑流或被蒸發(fā)[20]；因此在降水量持續(xù)減少的時期，大多數(shù)降水都用來補(bǔ)給土壤，只有少部分形成徑流，徑流量下降，導(dǎo)致基流量減少，這是1999年基流發(fā)生突變的原因之一。當(dāng)然，降水并不是影響基流突變的最主要因素(圖5)，降水對潮河流域1999—2015年基流下降的貢獻(xiàn)率只占18.90%(表2)。李子君等[21]通過分析潮河流域降水和徑流的變化關(guān)系得到，降水并不是影響徑流發(fā)生劇烈變化的主要原因。郭軍庭等[13]也證明降水不是影響潮河流域產(chǎn)流發(fā)生劇烈變化的最主要原因。然而 Li等[16]通過繪制基流與降水量之間的雙累積曲線計(jì)算氣候變化與人類活動對基流變化的貢獻(xiàn)率時，發(fā)現(xiàn)氣候變化對基流突變的影響占70.8%，降水是影響加拿大Similkameen河上游基流突變的主要原因；因此有必要對不同流域影響基流變化的因素進(jìn)行研究，這有利于科學(xué)評價和利用水資源。
4.2 人類活動對基流的影響
影響流域基流動態(tài)變化的主要人為活動包括土地利用變化、大規(guī)模的造林綠化等生態(tài)修復(fù)工程、水利工程措施建設(shè)、流域取用水等。該流域蓄水工程大部分建于20世紀(jì)70年代末至80年代初，小型水庫蓄水減少了水庫下游的河川徑流量，基流也相應(yīng)地減少。從20世紀(jì)80年代后，該流域迎來了大規(guī)模的水土保持綜合治理，1989年至2000年水利部將潮河流域列為國家級水土保持重點(diǎn)治理區(qū)，1998年國家林業(yè)局實(shí)施退耕還林工程，2001年后又開展一系列生態(tài)和治理項(xiàng)目[22]，這些工程項(xiàng)目的實(shí)施導(dǎo)致了潮河流域土地利用發(fā)生了重大變化。潮河流域主要土地利用類型為耕地、有林地、灌木林地和草地，其面積之和占研究區(qū)總面積97%以上。1987—2015年期間，潮河流域基流呈極顯著下降趨勢(P<0.001)，郭軍庭等[13]的研究顯示，同期流域有林地和灌木林面積分別增加了58.11%和26.90%，耕地和草地面積分別減少了27.75%和60.24%，比較突變點(diǎn)前后不同土地利用類型面積的變化率(圖6)，突變點(diǎn)后耕地和草地面積變化率分別比突變點(diǎn)前減少38.02%和19.41%，有林地和灌木林地面積變化率分別比突變點(diǎn)前增加4.40%、10.07%，水域、居民地和未利用地面積占比小，因此對流域基流的影響不做討論。一般而言，與草地和農(nóng)田相比，林地水分蒸騰作用更強(qiáng)，單位面積的林地比農(nóng)田和草地蒸發(fā)更多的水分，從而使地表徑流和基流減少[23]。Sun等[24]基于概念模型研究流域內(nèi)草地轉(zhuǎn)化為林地后的變化，結(jié)果發(fā)現(xiàn)草地轉(zhuǎn)換為林地后流域的平均徑流量減少。劉柏君等[25]基于SWAT模型模擬流域不同土地利用條件下基流的百分比，發(fā)現(xiàn)全裸、全草、全耕、全林地4種條件下全耕地的基流量最大，這是由于耕地的下層土壤更為緊實(shí)，入滲率低，因此更容易產(chǎn)生基流。
潮河流域的突變點(diǎn)大約發(fā)生在造林后的20年左右，在這期間，大量草地和耕地轉(zhuǎn)換為林地，潮河流域植被覆蓋率提高，植被的截留量和蒸發(fā)量增加，植被的截留可以減少降水對土壤水和地下水的補(bǔ)給，蒸發(fā)量的提升減少了土壤向地下水的補(bǔ)給，使潮河流域基流減小[26]。密云水庫上游建有26座小型水庫，在一定程度上增加了蒸發(fā)損失和滲漏量。1999年以來，北京經(jīng)歷7年大旱，密云水庫入庫量急劇減少，加劇了北京市人均供需水量的矛盾，潮河流域農(nóng)業(yè)用水、工業(yè)用水以及生活用水約37.6%來自地下水[23]，這也是造成流域內(nèi)基流減少的原因之一。雙累積曲線法求得的人類活動對基流的貢獻(xiàn)率為人類活動最大可能貢獻(xiàn)率，該法分離的是降水和非降水的水文要素，因而最終結(jié)果可能會低估氣候變化的貢獻(xiàn)率而高估人類活動的貢獻(xiàn)率。
圖6 潮河流域突變點(diǎn)前后土地利用面積變化率Fig.6 Change rate of landuse area before and after the abrupt change point of Chaohe watershed
5 結(jié)論
筆者通過潮河流域上游1963—2015年觀測數(shù)據(jù)，分析了降水、PET與基流的年際變化趨勢，從氣候變化和人類活動的角度，探究二者與基流變化的關(guān)系，并進(jìn)一步發(fā)掘了影響基流變化的主導(dǎo)因素。研究結(jié)果表明：1963—2015年潮河流域上游降水呈現(xiàn)減少趨勢、PET顯著增加，基流顯著減少，并且在1999年發(fā)生突變；在氣候變化方面，降水對基流的影響程度大于PET，降水的減少會導(dǎo)致基流的減少，相反，PET的增加會導(dǎo)致基流的減少；與氣候變化相比，人類活動是影響基流突變的主導(dǎo)因素，其貢獻(xiàn)率為81.10%。
猜你喜歡
基于濃度與流量突變的河流總磷通量估算農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(2024年3期)2024-11-06 00:00:00 潮關(guān)村：作別傷痕累累的過往北京紀(jì)事(2024年4期)2024-03-28 12:03:07 美國長期氣象干旱對河道基流的影響研究水利水電快報(2023年11期)2023-11-23 21:51:43 繪制和閱讀降水量柱狀圖《學(xué)習(xí)方法報》地理商務(wù)星球七年級(2023年13期)2023-10-11 02:44:46 降水量是怎么算出來的啟蒙(3-7歲)(2019年8期)2019-09-10 03:09:08 1988—2017年呼和浩特市降水演變特征分析江西農(nóng)業(yè)(2018年23期)2018-02-11 07:26:59 大潮河散文選刊·下半月(2017年10期)2017-10-23 22:05:53 基于改進(jìn)Tennant法的洋河流域生態(tài)基流估算研究水利科技與經(jīng)濟(jì)(2017年5期)2017-04-22 02:39:30 基于小波變換的三江平原旬降水量主周期識別水利科技與經(jīng)濟(jì)(2016年6期)2016-04-22 05:08:06 不同生態(tài)基流對水利樞紐發(fā)電損益的影響水利規(guī)劃與設(shè)計(jì)(2016年7期)2016-02-28 15:06:25
王曼玉，蔡永茂，李春征，張志強(qiáng)?
(1.北京林業(yè)大學(xué)水土保持學(xué)院，100083，北京；2.北京市八達(dá)嶺林場，102112，北京；3.北京市共青林場，101300，北京)
基流作為河川徑流的重要組成部分，是水量平衡分析、水資源優(yōu)化配置和水資源管理的重要內(nèi)容，也是維護(hù)河流生態(tài)健康的重要基礎(chǔ)[1]?；魍ǔＳ傻叵滤推渌舆t的水資源補(bǔ)給河川徑流的部分水量組成，氣候變化和人類活動對基流的影響已成為國內(nèi)外學(xué)者研究的熱點(diǎn)。白樂等[2]發(fā)現(xiàn)人類活動是影響禿尾河流域基流變化的主要因素；王敬哲等[3]發(fā)現(xiàn)呼圖壁河基流主要受溫度影響；郭軍庭等[4]計(jì)算蔡家川流域基流指數(shù)，發(fā)現(xiàn)基流指數(shù)隨降水的增加而減少，且農(nóng)地、灌叢、人工林對基流的影響為負(fù)效應(yīng)；L. Ahiablame等[5]計(jì)算了降水和農(nóng)業(yè)用地影響密蘇里河流域基流變化的彈性系數(shù)，結(jié)果發(fā)現(xiàn)降水的增加導(dǎo)致基流增加，而農(nóng)業(yè)用地增加則會導(dǎo)致基流減少。由于流域基流量不能通過直接測定的方法獲得，目前廣泛應(yīng)用的流域基流分割方法主要包括4大類：直接分割法、水量平衡法、同位素法和時間序列分割法[6]。其中，數(shù)字濾波法操作簡單，可以較好地模擬人工分割過程線，減少分割基流的主觀因素[7-8]。黃國如[9]通過對不同數(shù)字濾波法分割基流的結(jié)果進(jìn)行對比，認(rèn)為Chapman-Maxwell濾波法具有精確方便的優(yōu)點(diǎn)。
潮河是北京市主要飲用水源密云水庫的主要支流之一，近50年來該流域徑流在氣候變化與土地利用變化的影響下呈現(xiàn)顯著的下降趨勢[10]。郭軍庭等[11]和Wang等[12]定量研究了氣候變化和土地利用變化對潮河流域徑流量減少的貢獻(xiàn)率。然而，人類活動和氣候變化對基流的影響還缺乏深入的研究，筆者基于Chapman-Maxwell濾波法分割基流，并采用統(tǒng)計(jì)學(xué)方法，研究潮河流域密云水庫上游人類活動和氣候變化對基流的影響，旨在為認(rèn)識和評價水庫水資源變化、制定科學(xué)合理的水資源和流域管理方案提供科學(xué)依據(jù)。
1 研究區(qū)概況
潮河流域位于北京市與河北省的東部，E 116°7′～117°35′，N 40°19′～41°38′，發(fā)源于河北省豐寧縣，匯入北京市密云水庫。密云水庫上游流域總面積4 855.9 km2(圖1)。研究區(qū)氣候類型屬于溫帶季風(fēng)大陸性半濕潤半干旱氣候，四季分明，干旱和冷暖變化明顯，多年平均氣溫為9.16 ℃，年平均降水量490 mm，汛期6—9月之間，汛期降水量占總降水量75%以上，且多以暴雨形式出現(xiàn)。
圖1 潮河流域地形圖及雨量站、氣象站、水文站分布圖Fig.1 Terrain, precipitation station, meteorological station and hydrological station in Chaohe watershed
2 數(shù)據(jù)來源及研究方法
2.1 數(shù)據(jù)收集整理
筆者選取密云水庫上游潮河流域內(nèi)的4個雨量站、1個標(biāo)準(zhǔn)氣象站、1個水文站以及該流域周邊2個氣象站(圖1)的實(shí)測氣象數(shù)據(jù)、降水?dāng)?shù)據(jù)以及日徑流數(shù)據(jù)進(jìn)行分析。這些數(shù)據(jù)由河北水文總站、北京市水文總站和海河流域水利委員會提供。流域降水量采用4個雨量站實(shí)測資料由泰森多邊形法計(jì)算得到，潛在蒸散發(fā)量(potential evapotranspiration，PET)采用Penman-Monteith法計(jì)算得到。由于潛在蒸發(fā)散空間變異不大，其值為根據(jù)3個氣象站資料分別計(jì)算得到的PET的均值。流域土地利用分別解譯自1987年、1999年、2009年3期分辨率為30 m的Landsat TM影像遙感數(shù)據(jù)。將經(jīng)過地理坐標(biāo)配準(zhǔn)和幾何精度矯正的遙感數(shù)據(jù)與野外實(shí)測調(diào)查相結(jié)合，經(jīng)過人機(jī)交互判讀和監(jiān)督分類將土地利用類型分為有林地、灌木林地、草地、耕地、水域、居民區(qū)與建筑用地和未利用地等7類[13]。
2.2 研究方法
2.2.1 基流分割Chapman-Maxwell濾波法 濾波法是分割高頻信號和低頻信號的有效工具，流域快速響應(yīng)的直接徑流與高頻信號相似，慢速響應(yīng)的基流與低頻信號相似，該法被廣泛應(yīng)用于基流分割中[14]。其濾波方程為
(1)
式中：Qi為i時刻的總徑流，mm；Qb,i為i時刻的基流，mm；k為退水系數(shù)，一般情況下k取0.95時，基流分割結(jié)果最為合理[9，15]。
2.2.2 趨勢和突變點(diǎn)的Mann-Kendall檢驗(yàn)
1)趨勢檢驗(yàn)。
Mann-Kendall(以下簡稱MK)檢驗(yàn)被廣泛用于檢驗(yàn)水文、氣象等時間序列的趨勢[5,10]，MK檢驗(yàn)公式如下：
(2)
式中：e為樣本個數(shù)，Xn、Xm是連續(xù)時間序列上的值。當(dāng)Xn-Xm>0時，sgn=1; 當(dāng)Xn-Xm=0時，sgn=0；當(dāng)Xn-Xm<0時，sgn=-1。假設(shè)時間序列沒有變化，則S均值為0，計(jì)算方差V(S)的公式為
(3)
MK統(tǒng)計(jì)量公式如下：
(4)
通常用標(biāo)準(zhǔn)Z檢驗(yàn)表示數(shù)據(jù)變化的趨勢，對于給定的置信水平α，若|Z|≥Z1-α/2,即時間序列在置信水平α上存在明顯的上升或下降趨勢。當(dāng)Z值為正時，表示該時間序列上的值為增加趨勢，反之則為減小趨勢。
2)突變點(diǎn)分析。
設(shè)序列值為X1，X2…Xk，St表示第n個樣本Xn>Xm(1 (5) 在原時間序列獨(dú)立的假定下，St的均值和方差分別為： (6) (7) 將St標(biāo)準(zhǔn)化： (8) 2.2.3 雙累積曲線法 雙累積曲線法可以建立降水-徑流間的年累積值回歸關(guān)系，因此被廣泛應(yīng)用于計(jì)算降水與人類活動對年徑流量變化的相對貢獻(xiàn)率，筆者采用該方法分析將降水與人類活動對基流的影響程度[16]。步驟如下： 1)建立基準(zhǔn)期年累積降水量與基準(zhǔn)期年累積基流量的回歸方程 ∑Qb=a∑P+b。 (9) 式中：Qb為基準(zhǔn)期基流量，mm；P為基準(zhǔn)期降水量，mm；a為回歸方程的相關(guān)系數(shù)，b為回歸方程截距。 2)根據(jù)基準(zhǔn)期所得回歸方程，帶入研究期累積降水量求研究期累積基流值，再反推研究期的年基流值，則 r=Qb12-Qb1， (10) l=Qb2-Qb12。 (11) 式中：r為降水引起基流變化部分，mm；Qb12為利用基準(zhǔn)期所得到的回歸方程推算得到的研究期計(jì)算年基流值，mm；Qb1為基準(zhǔn)期分離基流，mm；l為人類活動引起基流變化部分，mm；Qb2為研究期分離基流，mm。 潮河流域1963—2015年降水、PET、徑流和基流呈現(xiàn)不同的變化趨勢(表1)。流域多年平均降水量為486.8 mm，以平均每年0.5 mm的幅度下降，但這一趨勢并不顯著(P=0.24)(圖2a)，并于2008年發(fā)生突變(圖3a)。PET以平均每年0.9 mm的速度顯著上升(P=0.02)(圖2b)，年平均PET為1 185.8 mm，突變點(diǎn)發(fā)生在2007年(圖3b)。流域徑流以平均每年0.6 mm的速度極顯著下降(P<0.001)(圖2c)，多年平均徑流量為49.5 mm，突變點(diǎn)位于1999年(圖3c)。流域基流呈極顯著下降趨勢(P<0.001)(圖2d)，平均基流量為10.9 mm，平均每年減少0.1 mm，且于1999年發(fā)生突變(圖3d)。由此可見，潮河流域基流量與徑流量的變化趨勢相同，但徑流量變化幅度大于基流量的變化幅度。降水和PET均未在1999年發(fā)生突變，說明影響基流量突變的因素除了氣候因素外還有人類活動的影響。 表1潮河流域年降水、PET、徑流和基流量的MK趨勢檢驗(yàn)及突變點(diǎn)分析 Tab.1MK trend analysis and change point of annual precipitation, PET, streamflow and baseflow in Chaohe watershed mm 圖2 潮河流域年降水量、年P(guān)ET、年徑流量、年基流量線性變化趨勢Fig.2 Linear trend of annual precipitation, annual PET, annual streamflow and annual baseflow in Chaohe watershed UFt為序列值；為反序列值；a為0.05顯著水平臨界值。UFt: Sequence value; Inverse sequence value; a significant horizontal threshold of 0.05. 圖3 潮河流域年降水、年P(guān)ET、年徑流量和年基流量突變點(diǎn)分析結(jié)果Fig.3 Change points of annual precipitation, annual PET, annual streamflow and annual baseflow in Chaohe watershed 降水和PET是氣候變化中影響基流的2個主要因素。潮河流域基流與降水呈顯著正相關(guān)關(guān)系(P<0.01)，相關(guān)系數(shù)為0.64(圖4a)；基流與PET呈顯著負(fù)相關(guān)關(guān)系(P<0.01)，相關(guān)系數(shù)為0.27(圖4b)。1963—2015年，降水量的下降和PET的上升共同導(dǎo)致了基流量的減少，但降水對基流的影響比PET對基流的影響更大?；魍蛔凕c(diǎn)前后，年平均降水量減少16.4%(84.4 mm),年平均PET增加7.3%(84.4 mm)，年平均基流減少66.8%(9.2 mm)。從上述結(jié)果可知，PET在突變點(diǎn)前后的變化微小，可以忽略不計(jì)。從水文循環(huán)的角度看，基流的顯著變化不僅是由于氣候變化引起的，人類活動在其中起到了重要作用。 圖4 潮河流域降水-基流、PET-基流關(guān)系圖Fig.4 Precipitation- baseflow and PET- baseflow in Chaohe watershed 以突變點(diǎn)1999年為界，將潮河流域年降水與基流系列分為變化前的1963—1998年(基準(zhǔn)期)和變化后的1999—2015年(變化期)進(jìn)行對比分析?；鶞?zhǔn)期年平均基流量為13.8 mm，變化期年平均基流量為4.6 mm，突變點(diǎn)前后基流量下降了66.8%(表2)。建立1963—2015年累積基流量與累積降水量回歸關(guān)系(P<0.001)(圖5)，1963—1998年2者相關(guān)性較高(R2=0.991 8)，因此可認(rèn)為1963—1998年基流變化是受人類影響較小的時期。將此時期作為基準(zhǔn)期，用來還原1999—2015年受人類影響較小時期的情況，若基流只受降水影響，雙累積曲線中的趨勢應(yīng)與基準(zhǔn)期的走勢保持一致，然而曲線于1999年發(fā)生向下偏轉(zhuǎn)，說明人類活動在很大程度上導(dǎo)致基流量減少。在不考慮其他氣候因素變化的情況下，人類活動對基流在突變點(diǎn)前后變化的最大可能貢獻(xiàn)率高達(dá)81.10%, 降水導(dǎo)致基流變化的貢獻(xiàn)率為18.90%(表2)?？梢?，人類活動是潮河流域突變點(diǎn)后基流減少的主導(dǎo)因素。 表2 降水和人類活動對潮河流域基流的影響Tab.2 Impacts of precipitation and human activities on baseflow in Chaohe watershed 圖5 潮河流域降水量-基流量雙累積曲線Fig.5 Double accumulative curve of precipitation-baseflow in Chaohe watershed 在全球變暖的總體趨勢下，潮河流域PET呈現(xiàn)顯著上升趨勢，降水量呈現(xiàn)不顯著下降趨勢，從而導(dǎo)致徑流量和基流量呈現(xiàn)極顯著下降的趨勢。盡管PET上升的趨勢較降水下降的趨勢更為顯著，但降水對基流的影響更為劇烈，這與D. L. Ficklin等[17]的結(jié)論一致。王云琦等[18]也指出流域內(nèi)降水量對產(chǎn)流量的影響大于蒸散量對產(chǎn)流量的影響。潮河流域基流量隨降水量的減小而減小，由于基流具有穩(wěn)定性和滯后性，在長時間尺度下，降水會通過影響產(chǎn)流下墊面特征，進(jìn)而對基流產(chǎn)生影響[19]。1999年以來，北京遭遇連續(xù)7年的干旱，降水量下降明顯，降水落地后先與土壤作用，待土壤飽和后轉(zhuǎn)換成徑流或被蒸發(fā)[20]；因此在降水量持續(xù)減少的時期，大多數(shù)降水都用來補(bǔ)給土壤，只有少部分形成徑流，徑流量下降，導(dǎo)致基流量減少，這是1999年基流發(fā)生突變的原因之一。當(dāng)然，降水并不是影響基流突變的最主要因素(圖5)，降水對潮河流域1999—2015年基流下降的貢獻(xiàn)率只占18.90%(表2)。李子君等[21]通過分析潮河流域降水和徑流的變化關(guān)系得到，降水并不是影響徑流發(fā)生劇烈變化的主要原因。郭軍庭等[13]也證明降水不是影響潮河流域產(chǎn)流發(fā)生劇烈變化的最主要原因。然而 Li等[16]通過繪制基流與降水量之間的雙累積曲線計(jì)算氣候變化與人類活動對基流變化的貢獻(xiàn)率時，發(fā)現(xiàn)氣候變化對基流突變的影響占70.8%，降水是影響加拿大Similkameen河上游基流突變的主要原因；因此有必要對不同流域影響基流變化的因素進(jìn)行研究，這有利于科學(xué)評價和利用水資源。 影響流域基流動態(tài)變化的主要人為活動包括土地利用變化、大規(guī)模的造林綠化等生態(tài)修復(fù)工程、水利工程措施建設(shè)、流域取用水等。該流域蓄水工程大部分建于20世紀(jì)70年代末至80年代初，小型水庫蓄水減少了水庫下游的河川徑流量，基流也相應(yīng)地減少。從20世紀(jì)80年代后，該流域迎來了大規(guī)模的水土保持綜合治理，1989年至2000年水利部將潮河流域列為國家級水土保持重點(diǎn)治理區(qū)，1998年國家林業(yè)局實(shí)施退耕還林工程，2001年后又開展一系列生態(tài)和治理項(xiàng)目[22]，這些工程項(xiàng)目的實(shí)施導(dǎo)致了潮河流域土地利用發(fā)生了重大變化。潮河流域主要土地利用類型為耕地、有林地、灌木林地和草地，其面積之和占研究區(qū)總面積97%以上。1987—2015年期間，潮河流域基流呈極顯著下降趨勢(P<0.001)，郭軍庭等[13]的研究顯示，同期流域有林地和灌木林面積分別增加了58.11%和26.90%，耕地和草地面積分別減少了27.75%和60.24%，比較突變點(diǎn)前后不同土地利用類型面積的變化率(圖6)，突變點(diǎn)后耕地和草地面積變化率分別比突變點(diǎn)前減少38.02%和19.41%，有林地和灌木林地面積變化率分別比突變點(diǎn)前增加4.40%、10.07%，水域、居民地和未利用地面積占比小，因此對流域基流的影響不做討論。一般而言，與草地和農(nóng)田相比，林地水分蒸騰作用更強(qiáng)，單位面積的林地比農(nóng)田和草地蒸發(fā)更多的水分，從而使地表徑流和基流減少[23]。Sun等[24]基于概念模型研究流域內(nèi)草地轉(zhuǎn)化為林地后的變化，結(jié)果發(fā)現(xiàn)草地轉(zhuǎn)換為林地后流域的平均徑流量減少。劉柏君等[25]基于SWAT模型模擬流域不同土地利用條件下基流的百分比，發(fā)現(xiàn)全裸、全草、全耕、全林地4種條件下全耕地的基流量最大，這是由于耕地的下層土壤更為緊實(shí)，入滲率低，因此更容易產(chǎn)生基流。 潮河流域的突變點(diǎn)大約發(fā)生在造林后的20年左右，在這期間，大量草地和耕地轉(zhuǎn)換為林地，潮河流域植被覆蓋率提高，植被的截留量和蒸發(fā)量增加，植被的截留可以減少降水對土壤水和地下水的補(bǔ)給，蒸發(fā)量的提升減少了土壤向地下水的補(bǔ)給，使潮河流域基流減小[26]。密云水庫上游建有26座小型水庫，在一定程度上增加了蒸發(fā)損失和滲漏量。1999年以來，北京經(jīng)歷7年大旱，密云水庫入庫量急劇減少，加劇了北京市人均供需水量的矛盾，潮河流域農(nóng)業(yè)用水、工業(yè)用水以及生活用水約37.6%來自地下水[23]，這也是造成流域內(nèi)基流減少的原因之一。雙累積曲線法求得的人類活動對基流的貢獻(xiàn)率為人類活動最大可能貢獻(xiàn)率，該法分離的是降水和非降水的水文要素，因而最終結(jié)果可能會低估氣候變化的貢獻(xiàn)率而高估人類活動的貢獻(xiàn)率。 圖6 潮河流域突變點(diǎn)前后土地利用面積變化率Fig.6 Change rate of landuse area before and after the abrupt change point of Chaohe watershed 筆者通過潮河流域上游1963—2015年觀測數(shù)據(jù)，分析了降水、PET與基流的年際變化趨勢，從氣候變化和人類活動的角度，探究二者與基流變化的關(guān)系，并進(jìn)一步發(fā)掘了影響基流變化的主導(dǎo)因素。研究結(jié)果表明：1963—2015年潮河流域上游降水呈現(xiàn)減少趨勢、PET顯著增加，基流顯著減少，并且在1999年發(fā)生突變；在氣候變化方面，降水對基流的影響程度大于PET，降水的減少會導(dǎo)致基流的減少，相反，PET的增加會導(dǎo)致基流的減少；與氣候變化相比，人類活動是影響基流突變的主導(dǎo)因素，其貢獻(xiàn)率為81.10%。3 結(jié)果與分析
3.1 降水、PET、徑流和基流變化趨勢和突變點(diǎn)
3.2 降水和PET對基流的影響
3.3 降水和人類活動對基流的影響
4 討論
4.1 氣候因素對流域基流的影響
4.2 人類活動對基流的影響
5 結(jié)論
猜你喜歡
農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(2024年3期)2024-11-06 00:00:00
北京紀(jì)事(2024年4期)2024-03-28 12:03:07
水利水電快報(2023年11期)2023-11-23 21:51:43
《學(xué)習(xí)方法報》地理商務(wù)星球七年級(2023年13期)2023-10-11 02:44:46
啟蒙(3-7歲)(2019年8期)2019-09-10 03:09:08
江西農(nóng)業(yè)(2018年23期)2018-02-11 07:26:59
散文選刊·下半月(2017年10期)2017-10-23 22:05:53
水利科技與經(jīng)濟(jì)(2017年5期)2017-04-22 02:39:30
水利科技與經(jīng)濟(jì)(2016年6期)2016-04-22 05:08:06
水利規(guī)劃與設(shè)計(jì)(2016年7期)2016-02-28 15:06:25
