華北水資源年代際變化及其與全球變暖之間的關聯

張書萍 祝從文 周秀驥

?

華北水資源年代際變化及其與全球變暖之間的關聯

張書萍 祝從文 周秀驥

中國氣象科學研究院，北京100081

本文利用1951～2011年期間中國臺站資料、東亞地區的探空資料、NCEP/NCAR和ERA40等大氣再分析資料，通過對水分平衡方程診斷探討了華北地區過去60年中水資源和可利用降水量的變化特征及其與大尺度環流變化之間的關系。結果發現，華北地區69%的降水量被蒸發，可利用降水僅僅為降水量的31%。夏季可利用降水是華北水資源的主要來源，華北夏季可利用降水量在80年代初發生突變減少，進入21世紀初，伴隨蒸發量的增加該地區可利用降水量進一步減少。西風帶水汽與東亞夏季風水汽是華北可利用降水的主要來源。NCAR/NCEP 和EAR40再分析資料的結果均顯示貝加爾湖一帶的位勢高度偏低和西太平洋高壓的偏強有利于該地區降水增多。利用探空資料進一步證明，蒙古以及貝加爾湖地區的溫度在對流層低層變暖和位勢高度場的加強導致了過去幾十年華北可利用降水量減少。由于貝加爾湖地區溫度變化與全球變暖存在密切關聯，監測該地區溫度的變化對預測華北水資源和東亞夏季風的長期變化具有重要的意義。

華北水資源 可利用降水 年代際變化 全球變暖

1 引言

中國水資源的變化與東亞季風的變化密切相關。大量的證據表明，伴隨全球變暖，東亞季風發生了顯著的年代際調整，與此同時中國水資源在空間分布上也發生了明顯改變，這一改變不僅加劇了中國南澇北旱的空間格局，也引起了一系列社會和科學問題并受到各領域的科學家的廣泛關注。我國水資源系統對氣候變化的承受能力十分脆弱，多數河流的徑流對大氣降水變化非常敏感（任國玉等，2008）。水資源的重新分配將對我國農業、牧業、漁業、航運、水力發電等其他部門產生進一步的影響。事實表明，全球變化加劇我國水資源南北不平衡的現狀，華北大部、華東北部和東北東部降水均呈減少趨勢，干旱化嚴重，觀測顯示海河流域入海流量逐年減少，黃河在20 世紀90 年代后發生持續的嚴重斷流（葉柏生，2004）。而長江中下游地區及東南沿海一帶洪澇災害出現頻繁，長江流域出現1998年特大洪澇災害等（丁一匯和胡國權，2003）。這些災害給人民的生命財產都帶來巨大威脅。中國水資源主要來自東亞夏季風降水，而東亞季風變異所引起降水的空間分布改變直接導致了水資源空間配置狀態的變化。華北地區是我國水資源十分貧乏的地區之一，夏季降水異常引起該地干旱化的問題一直是氣象學者們關注的焦點。研究表明，華北地區夏季降水在1960年之后發生了持續性的減少，并且具有顯著的年際和年代際變率，其轉折性的干濕變化發生在上世紀70年代末，與 東亞夏季風的持續減弱存在緊密聯系（黃榮輝 等，1999；周連童和黃榮輝，2003；徐桂玉等，2005；馬柱國和符淙斌，2006；徐康等，2011a；Zhu et al., 2012）。

天然降水中的一部分會通過蒸發的形式返回到大氣中，剩余的部分通過地表徑流的形式補充地下水。因此，天然降水雖然是水資源主要的補給來源，卻并不能等同于當地的可利用水資源。由于蒸發的水分進入大氣之后難以加以利用，因此降水中只有減去蒸發的部分才是可利用的水資源，也稱可利用降水量。宋正山等（1999）研究了華北的可利用降水量變化，發現華北地區年平均的80%左右的降水被蒸發掉，僅有20% 左右的降水滯留在 土壤、湖泊中或者可以成為河流的徑流中加以利用。研究發現，近幾十年來華北地區一方面降水急劇減少，另一方面由于氣溫不斷上升，潛在蒸發量明顯加大。特別是20世紀80年代以來，華北地區持續干旱導致該地區的水資源供需矛盾更加嚴重。

區域降水的長期變化與大尺度環流背景下的水汽輸送特征和局地不穩定度密切相連。大氣水汽含量的異常和水汽輸送的多寡是造成降水異常的重要因素之一。東亞季風區是全球大氣含水量最高的地區之一，東亞季風在全球大氣水分循環中起到了重要的作用。研究表明，中國東部地區夏季降水不僅受到來自低緯熱帶地區的印度西南季風及南海季風的影響，還受到西太平洋副熱帶西南季風的影響。并且研究發現降水異常年對應的水汽來源和正常年對應的水汽來源明顯不同（Zhou and Yu，2005），其中華北降水異常偏多年對應中緯度西風帶異常輸送與副高西北側的西南氣流異常輸送在華北上空的輻合（田紅等，2002）。大尺度環流場在降水異常年中也有明顯的特征，有證據表明，旱澇過程主要發生在異常環流具有明顯持續性和穩定性的時期，特別是干旱，常常是某種狀態的異常環流持續發展和長期穩定的結果。張慶云等（1999）研究表明，華北降水的年際和年代際變化始于夏季東亞—西太平洋地區上空500 hPa位勢高度場以及西太平洋副熱帶高壓脊線位置的年際和年代際變化。Wang（2001）揭示了1970年代末對流層溫度和風場發生的轉變，在這次轉變之后亞洲和非洲的季風環流減弱，降水也發生了相應的變化。郭其蘊等（2003）研究了東亞夏季風的年代際變率對中國氣候的影響，指出了東亞夏季風在1970年代中期發生了年代際減弱，華北的降水也隨之由多變少。宇如聰（2008）等指出，華北干旱和長江流域多雨的趨勢是東亞對流層高層變冷造成的，變冷與對流層和平流層相互作用相關。最近的研究結果發現，中國夏季降水南澇北旱格局的長期變化與溫室氣體導致的貝加爾湖地區地表增溫存在密切聯系（Zhu et al., 2012）。

華北地區地處東亞夏季風的北緣地帶，東亞季風水汽對其水資源的變化起著重要的作用。在最近幾十年全球變暖長期變化背景下，華北地區的可利用水資源發生了怎樣的變化？導致這種變化的環流特征是什么？為此，本文首先利用中國臺站觀測的降水和溫度資料，結合現有的大氣再分析資料和臺站探空資料，通過分析華北地區空中水資源的垂直分量即降水與蒸發的變化特征，揭示出華北地區可利用水資源的變化特征及其與全球變暖之間的關聯。

2 資料與方法

本文用到的資料有1951～2011年的NCEP/ NCAR再分析月平均資料（Kalnay et al., 1996）、歐洲數值預報中心的大氣再分析資料（Uppala et al., 2005）。兩套資料的水平分辨率為2.5°×2.5°。除此之外，本文還用到了國家信息中心提供的1951～2011年臺站觀測的逐日降水和溫度資料以及亞洲地區1951～2007年月平均探空資料集。本文利用1971～2000年30年的平均值代表氣候態特征。顯著性檢驗采用檢驗方法。

我們利用探空資料計算了大氣的水汽通量。水汽平衡方程表示如下：

緯向和經向水汽輸送矢量Q與Q的計算公式分別表示為

, （3）

其中，為該單位氣柱各層大氣的風速矢量，和分別為東西風和南北風分量，是該單位氣柱各層大氣的比濕，為地面氣壓。

整層大氣的水汽輸送通量的散度計算公式表示為

將水汽通量輻合項進行離散化，則局地的水汽收支可以通過如下公式計算得出：

（6）

其中，σ為單位面積，分別代表南邊界、北邊界、東邊界、西邊界的格距長度。由于水汽在300 hPa以上含量很小可以忽略，因此在本文中計算和水汽相關的物理量整層累加值，=300 hPa。

對于較長時間（月、年），較大范圍的平均狀況而言，大氣的局地水汽儲存率非常小，較之其他項要小一個量級，可以忽略不計，因此水汽方程可以簡化為

上式說明了某個地區某一個季節或者某一個月內的降水主要來源為同季節時段中流經該區水汽的輻合以及在此地區的蒸發而來的水汽的總和。

蒸發是一個很難估算的分量，由于其在實際觀測中很難得到，很多學者都用經驗公式來估算實際蒸發量。研究發現，高橋浩一郎經驗公式在干旱與半 干旱地區對蒸發的估算較為準確（宋正山等，1999）。因此，本文采用高橋浩一郎的經驗公式來計算實際蒸散發（高橋浩一郎，1979），公式表示如下：

3 華北氣候和空中水資源的季節變化特征

為了對華北地區全區的水資源的氣候演變特征有一個整體的認識，我們首先對華北地區全區平均的年平均氣溫、年降水量、年蒸發量以及年可利用降水量的變化特征進行分析討論。本文利用全國722個臺站降水、溫度資料數據集，選取的華北地區范圍及降水站點（見圖1）。其中，華北地區的范圍定義為（35°～42.5°N，110°～120°E）區域，利用該區域平均降水代表華北降水量，計算華北降水量同全國其他站點降水量的相關，發現華北降水量的變化同內蒙古中東部、東北西部及華北地區的降水變化為正相關，相關性超過顯著性檢驗，表明所選范圍降水量變化的一致。

由于可利用水資源為降水量與蒸發量之差，而蒸發量又受到溫度變化的影響，因此，對可利用水資源的各分量的研究也至關重要。圖2表示可利用水資源及其各分量的季節變化，表1為各個分量季節平均的統計量。華北地區全區氣候的平均狀況是：氣候平均溫度為10.3°C，并具有顯著的夏高冬低的季節分布特征。平均降水量為552 mm，夏季占全年降水量的65%，為358 mm；隨后依次為秋季101 mm，占全年18%；春季77 mm，占全年14%；冬季降水量最小，為16 mm，僅占全年的3%。顯而易見，華北地區夏季降水是全年降水量的主要來源。

圖1 華北全區平均降水量與周邊臺站觀測的降水量的相關系數分布。陰影部分為通過95%顯著性檢驗

圖2 華北全區1971～2000年氣候平均（a）溫度、降水量、蒸發量以及（b）蒸發系數和可利用降水系數逐月變化

表1 1951～2010年華北地區季節和年平均降水量（P）、蒸發量（E）、可利用降水量（P?E）、蒸發系數（α）和可利用降水系數（β）

本文利用蒸發和可利用降水系數分別表示大氣蒸發量和可利用降水相對于大氣降水量的大小。其中，蒸發系數定義為：；可利用降水系數定義為：（宋正山等，1999）。因此，蒸發愈多，可利用的降水資源愈少。

分析可以發現：華北地區區域平均的蒸發量總量為379 mm，冬夏季節差異較大，7月為最大蒸發月份。從季節平均來看，其夏季蒸發量為231 mm，占全年蒸發總量的61%，其次為秋季73 mm，占全年19%；春季62 mm，占全年16%；冬季蒸發最小，為13 mm，占全年3%。夏季蒸發最多，為主要蒸發季節。

多年氣候平均來看，華北地區區域平均的蒸發系數為0.69，可利用降水系數為0.31，表明華北地區全年降水量中有69%被蒸發返還到大氣中，可利用降水量僅為全年降水的31%。因此華北地區當地降水量的一多半用于蒸發消耗，而無法被人們使用。而從蒸發系數的季節變化中可以發現，華北地區最小的蒸發系數發生在夏季，僅為0.65，冬季、春季反而較高，最高為0.81。這一數據表明，夏季降水較多，伴隨的蒸發量也較高，但是對可利用降水的影響卻是最小。相反，在春、秋季雖然是蒸發總量不高的季節，但是蒸發的影響反而更大。也就是說：降水少的季節，蒸發的損失也就大。華北地區地面蒸發這一特征恰好加劇水資源的供需矛盾，特別是當地春季正是當地牧草返青、作物生長的時節，春灌等農事活動又非常需要水資源，較大的蒸發量引發可利用降水減少，引起更加嚴重的干旱災害（宋正山等，1999）。

由上可知，蒸發量損失占當地降水量中的比例極大，單憑當地降水量來描述華北水資源的整體變化并不十分恰當。本文用可利用降水量表示華北水資源的變化，可利用降水量為降水與蒸發之差計算得到，在水分平衡公式（公式7）中，可利用降水量代表水資源的垂直分量，其變化特征具有很大的研究價值。分析發現，華北地區的全區年平均可利用降水量為173 mm，其中夏季可利用降水量占全年的73%；依次為秋季、春季分別占全年的16%與9%；冬季最小，可利用降水量僅有3 mm，占全年的2%。由此可見，受到蒸發的影響，華北地區的可利用降水進一步集中在夏季，夏季可利用降水量已成為華北地區水資源的主要來源。然而，伴隨全球變暖和東亞夏季風的持續減弱，華北夏季潛在蒸發量不斷增加，與之相反的是華北夏季降水量持續減少，這一反向變化直接導致華北地區可利用降水量進一步減少，也使得水資源不平衡狀況更加明顯。

我們注意到，華北全區降水量、蒸發量以及可利用降水量的區域分布并不一致。由于夏季可利用降水是華北地區可利用降水量的主要來源，因此本文重點分析夏季的可利用水資源分布特征。圖3表示的是華北地區多年平均夏季蒸發、降水、可利用降水量以及蒸發系數等分量的區域分布特征。如圖所示，華北地區夏季降水明顯具有西北—東南的階梯分布特征，其分布的水平梯度并不均勻，東南部梯度大，西北部梯度小，但是在京津唐地區卻有一個小的高值中心。華北地區夏季蒸發的分布特征同降水相似，但緯向分布性比降水量區域分布均勻，這同蒸發量是選取經驗公式計算值有關，因此同實際情況必然存在一定誤差。華北可利用降水量也顯示出東南—西北遞減的分布形式，同降水量的分布比較一致，表明華北東南地區的蒸發雖然較多，但是其降水量也同樣偏多，可利用降水量相對于西北地區仍然充足。華北西北地區，特別是內 蒙古中部的可利用降水量不足50 mm，蒸發系數卻能達到0.8以上，可利用降水量稀少，水資源十分匱乏。

圖3 華北地區1971～2010年夏季平均：（a）降水量；（b）蒸發量；（c）可利用降水量；（d）蒸發系數的區域分布

4 華北氣候和水資源年代際變化特征

眾多研究表明，華北地區年平均溫度持續升高，同全球氣候變暖一致，降水量持續減少，并具有顯著的年代際變化。為了更清楚理解華北地區的可利用水資源在全球變暖背景下的年代際變化特征，本文分別計算了每10年間各分量的統計平均值，并結合氣溫、降水、蒸發以及可利用降水的逐年變化討論華北氣候和空中水資源年代際變化特征。圖4a表示的是華北地區不同季節區域平均氣溫的年代際變化。如圖所示，華北年平均氣溫波動升高，1950～1980年代均低于氣候平均值。然而1990年代，華北氣溫驟然升高，10年平均氣溫升高0.76°C，高于氣候平均值，轉為正距平。21世紀初華北氣溫繼續攀升，達到1960年來歷史最高值。M-K突變檢驗顯示，華北地區溫度出現突變升高，突變發生在1995年左右（圖略）。我們注意到，華北地區溫度變化在各個季節表現不盡相同。春季增溫發生在1980年代，持續到21世紀初達到最大，其升溫幅度為四個季節中最大。夏季增溫發生在1990年代，1980年代反而出現氣溫降低趨勢同其他季節表現不同，似乎有先抑后揚的現象。秋季增溫自1960年代起出現，升溫平穩緩慢，并無大起大落的現象。冬季華北溫度自1960年代后也持續平穩上升，增幅高于秋季，而在1990年代起增溫明顯。綜上所述，華北地區全區年平均溫度逐年升高由春、秋、冬三季貢獻，其中1950～1960年代溫度除春季外各季節均有一定程度的下降趨勢，溫度變化最大的轉折發生在1980年代。1960年代的增溫，冬季增暖為主要的貢獻，1980～1990年代的增溫為四個季節同時溫度變暖，冬季增溫幅度最大，其次為夏季、春季、秋季。這次四季同時變化的形勢很可能是一次氣候的突變情況。而在1990年代到21世紀初的增暖主要是由于春季增暖而引起的，這進一步加劇了春季水資源的矛盾需求。

伴隨華北溫度變暖，該地區的年平均降水量持續減少，與溫度的變化趨勢恰好相反（見圖4b）。如圖所示，全年平均降水量1950～1970年代均高于氣候平均值，1980年代突然減少轉為負距平，10年平均減少了41.7 mm。伴隨降水量持續減少，21世紀初達到降水量的氣候最低值。M-K突變檢驗顯示，華北地區降水量持續減少，突變發生在上世紀1970年代末和1980年代初。同樣，華北降水變化各個季節表現不盡相同。春季降水波動變化，并未有較大減少現象，1970年代為降水量最低時期。夏季降水量變化同全年基本一致，1980年代突變減少為負距平，此后持續降低，并在21世紀初進一步減少。秋季降水波動減少，趨勢并不明顯。冬季降水量值最小，相對其他季節變化并不明顯。由此可見，華北地區全區夏季降水量逐年減少對年降水的影響至關重要。從全年降水變化來看，1970年代降水量的減少主要是由于春季降水大幅度減少導致，全年降水在21世紀初減少主要由于夏季降水大幅度減少的貢獻，而春、秋季降水反而有所增加。

圖4 華北區域平均的氣候要素每10年距平變化：（a）溫度；（b）降水量；（c）蒸發量；（d）可利用降水量

蒸發受到降水量和溫度的共同影響，其年代際變化表現出獨有的特征（見圖4c）。從10年平均來看，1950年代蒸發量最高，1960年代發生急劇減少，1970年代略增，1980年代蒸發再一次劇減至負距平，而21世紀初蒸發反而增加，同降水量的變化趨勢相反，同溫度變化趨勢相同。除此之外，我們發現，華北蒸發變化各個季節表現不盡相同。春季蒸發量自1970年代劇減后一直波動增加。夏季蒸發量波動減少，而在21世紀初不降反增。秋季蒸發在1950～1970年代高于氣候平均值，1980年代劇減9 mm，低于氣候平均值，隨后又持續增加，21世紀初達到正距平。冬季蒸發波動變化，21世紀略有增加轉為正距平。由此可見，蒸發量的量值小于降水量，其變化趨勢同降水量有相似之處，但是并非完全一致。由于夏、秋、冬季蒸發的共同減少導致了蒸發量在1980年代的大幅度減少，但是同時期的春季蒸發卻有增多的趨勢。在21世紀初，蒸發量具有顯著增多的變化，同降水量的變化趨勢相反，這一增多趨勢在四個季節均有體現，最大的增多貢獻是春季的蒸發量的增加。蒸發量的增加會導致可利用降水量進一步減少。分析可利用降水的變化可以發現（圖4d），華北地區 1950～1960年代增加了14.1 mm，從1960年代開始，可利用降水量持續減少，分別在1960～1970年代減少了29.4 mm，1970～1980年代減少了33 mm，1990年代至21世紀初減少了20.7 mm。可利用降水量的各個季節的年代際變化趨勢特征同降水量非常一致，并且夏季的可利用降水量的變化幾乎決定了全年的整體變化特征，其他季節的可利用降水量的貢獻基本可以忽略不計。夏季的可利用降水量從1950年代開始持續減少，其減少量最大的兩個階段分別是1970～1980年代，減少33.6 mm，以及1990年代至21世紀初，減少27.6 mm。其他季節的可利用降水量的年代際變化特征并不明顯。可以看出，華北水資源的嚴重減少主要是由于其夏季可利用水資源的年代際減少而引起的。

5 影響夏季華北水資源的水汽和環流特征

從大氣水分平衡公式中可以發現，空中水資源的垂直分量需要其水平分量的補償。因此，華北地區夏季可利用降水量與同期水汽輸送有著密切的關系。為了描述這種關系，本文利用華北地區夏季可利用降水量與同期的水汽通量合成矢量進行了相關計算（見圖5）。由于NCEP/NCAR資料顯示上世紀70年代以前水汽輸送存在異常偏大情況（圖略），本文的資料選取1971～2010年，此處限于ERA-40的資料的長度，下圖相關選取了兩套資料的共同長度1971～2000年。

圖5a表示的是利用NCEP/NCAR大氣再分析資料計算結果，圖中發現同華北地區夏季可利用降水量相關的水汽輸送矢量顯著相關有兩部分：一部分為以蒙古高原為中心的逆時針水汽輸送通量相關系數，表示當華北夏季可利用降水量增加時，蒙古地區出現一個逆時針水汽輸送環流，顯著的相關系數分布表明西風帶水汽從新疆等西北地區向華北地區輸送的水汽通量增加；另一部分顯著相關的地區是由孟加拉灣北上的水汽，匯合了西太平洋副熱帶高壓西北側的水汽輸送通量，共同由南向北輸送到華北，同樣對應了華北地區的夏季可利用降水量偏多時期。圖5b表示的是基于ERA-40大氣再分析資料計算結果，對比NCEP/NCAR資料所得到的結果，我們發現：ERA- 40結果同樣出現了西風帶上從新疆經內蒙古到達華北的西風帶水汽輸送的高相關系數區，表明華北夏季可利用降水量增大同這一帶狀區域中的水汽輸送通量增大相關，但是其在蒙古以北地區的相關并不顯著。與NCEP資料比較一致，華北南部出現了由南向北的高水汽輸送相關區，但是其孟加拉灣水汽的相關也沒有NCEP資料計算結果顯著。

圖5 1971～2000年華北全區夏季可利用降水量與同期水汽輸送通量相關系數矢量：（a）NCEP/NCAR計算結果；（b）EAR-40計算結果。陰影為通過顯著性檢驗的區域

圖6表示的是1971～2010華北夏季可利用降水量與同期大氣環流之間的相關系數。從500 hPa位勢高度的相關系數場發現：顯著的負相關中心位于貝加爾湖西南部的蒙古地區，與之對應的為北極地區以泰梅爾半島為中心的正相關中心。同時超過顯著性檢驗的負相關系數區覆蓋低緯度的大部分地區，而在日本南部出現了正相關中心，但是未通過顯著性檢驗。這一相關分布表明，影響華北地區最主要的是以蒙古為中心的位勢高度異常，當華北地區水資源增加時，以蒙古地區為中心的位勢高度偏低。反之，這一帶的位勢高度偏高對應同期華北可利用降水量偏少。

圖6 1971～2010華北夏季可利用降水量與同期（a）500hPa位勢高度場、（b）200hPa緯向風場、（c）700hPa風矢量場、（d）110°E～120°E平均的垂直速度場的相關系數。陰影部分為通過顯著性檢驗

分析200 hPa緯向風場的相關系數場發現，緯向風場表示高空西風急流對華北地區降水的影響，從相關系數場中看到，環貝加爾湖地區的風場表現為顯著的負相關，在里海到新疆地區風場表現出顯著的正相關，這種特征還表現在東北地區。而在日本南部海域則表現出顯著的負相關區。從相關系數分布可以推斷，當華北降水偏多時，位于中國北方地區的西風帶增強，其南邊與北邊的西風帶均減弱，對應于500 hPa位勢高度場出現的低位勢高度，華北降水偏少時對應的緯向風場特征與之相反。分析700 hPa風矢量場的相關系數場發現：風速相關矢量場顯示出同位勢高度場一致的特征。即在以貝加爾湖西南為中心的區域出現逆時針風場相關，表示當以貝加爾湖西南為中心的區域出現逆時針風場時候，其攜帶的西風帶水汽也向華北地區輸送增多，夏季華北降水偏多。同時顯著相關風場還包括由南海向北輸送的分量，表明當華北可利用降水偏多時，東亞夏季風增強，夏季風水汽輸送也偏強。華北地區夏季可利用降水偏少時，風場反之。分析華北夏季可利用降水同110°E～120°E平均的垂直速度的相關系數場發現：對應華北地區可利用降水偏多，華北地區垂直速度場出現整層的垂直上升的增加，同時其南部出現垂直上升減弱。

研究發現，以貝加爾湖南部為中心的區域位勢高度異常對華北地區的降水存在較強的影響（徐康等，2011a，2011b）。因此，我們利用華北可利用降水同垂直的溫度以及環流特征相關來研究垂直環流與華北可利用降水之間的關系。圖7表示的是基于NCEP/NCAR大氣再分析資料計算的華北可利用降水指數與110°E～120°E平均的大氣環流場的相關系數，陰影為通過顯著性檢驗的部分。分析可以發現，華北可利用降水同華北及其以北地區的低層溫度為負相關，而在同樣緯度的700 hPa以上高層為正相關（圖7a），并且這一正相關區域隨高度的降低而南偏，在低層覆蓋我國江南地區。這一特征表明，高緯度溫度變化具有高低層反向變化的特征，高緯度蒙古以及貝加爾湖地區低層溫度升高，同華北可利用降水變化反向，而在高層溫度降低，同華北可利用降水變化一致。圖7b為同一區域的位勢高度的相關關系，從圖中看到，華北可利用降水同高度場呈現負相關，負相關區域在華北北部地區從低層到300 hPa層，而在華北當地上空300 hPa為中心出現正相關關系，表明這一地區的位勢高度同降水變化相一致，為降低趨勢。圖7c為同緯向風場的相關，從圖中看到，緯向風場出現了和位勢高度場相應一致的變化，在位勢高度場的最大相關中心的南、北兩側，分別出現了西風的正相關與負相關區域，表明當華北降水減少，位勢高度升高，在變化的最大中心南部出現緯向風的加速，在其北部出現緯向風的減弱。圖7d為同經向風場的相關關系，圖中看到，華北可利用降水同經向風的相關關系從21°N以北全部為負相關，并且通過顯著性檢驗，表明經向風由南向北均為減弱趨勢。

圖7 基于NCEP/NCAR大氣再分析資料計算的華北夏季可利用降水與110°E～120°E平均的（a）溫度場、（b）位勢高度場、（c）緯向風場、（d）經向風場相關的緯度—高度剖面。陰影為通過95%顯著性檢驗的區域

我們注意到，大氣再分析資料的可靠性一直以來存在爭議（李建等，2010）。為了檢驗NCEP/ NCAR資料的準確性，我們利用中國臺站探空資料進行了類似的分析（見圖8）。由于探空資料對風速的觀測是在1980年之后，因此風場相關的時段為1980～2007年。如圖8所示，探空資料出現了同NCEP資料幾乎一致的變化，因此，采用NCEP/ NCAR大氣再分析資料所描述的影響華北的環流和水汽特征具有較高的可信度。

圖8 探空資料計算華北夏季可利用降水與110°E～120°E平均的（a）溫度場、（b）位勢高度場、（c）緯向風場、（d）經向風場相關的緯度—高度剖面。陰影為通過95%顯著性檢驗的區域

6 結論和討論

本文利用1951～2010年中國華北地區臺站降水、溫度資料、NCEP/NCAR 和ERA-40再分析格點資料等，根據高橋浩一郎（1979）給出蒸發計算的經驗公式，計算和分析了近幾十年華北地區氣候和水資源年代際變化，揭示了影響華北可利用降水的水汽輸送和環流變化。主要結論如下：

（1）華北地區降水量69%被蒸發損耗，可利用降水僅占全年降水的31%。近60年以來華北空中水資源表現出明顯的減少趨勢。雖然最近10年該地區春季以及秋季降水增加，但全年降水依然表現出減少趨勢。華北地區蒸發量與降水量變化趨勢基本一致，蒸發在夏、秋、冬季均表現出顯著 減少趨勢，但1980年代之后春季的蒸發卻有所增多。進入21世紀，該地區的大氣蒸發量顯著增多。其中，春季蒸發量的增加是對這一趨勢的最大貢獻。

（2）華北夏季可利用降水量基本代表全年可利用降水量的變化特征。夏季及年平均的可利用降水量自1950年代起便持續減少，減少量最大的兩個階段是70～80年代（減少33.6 mm）以及90年代到21世紀初（減少27.6 mm）。華北水資源的嚴重減少主要是由于夏季可利用水資源的年代際減少而引起的。

（3）西風帶水汽與東亞夏季風水汽是華北可利用降水的主要來源，其中蒙古地區的東亞高緯度環流異常對華北夏季降水具有重要的影響。當華北地區可利用降水減少時，以蒙古地區為中心的位勢高度場異常偏強，西風急流在華北地區減弱，低緯度風場出現順時針環流特征，垂直速度場也出現下沉氣流。反之，當該地區位勢高度異常偏弱時，西風急流在華北地區增強，低緯度風場出現逆時針環流，垂直速度場表現為上升運動增強，對應同期華北可利用降水量增多。

已有研究表明，蒙古以及貝加爾湖地區的溫度在低層變暖與全球變暖存在密切關聯。我們注意到，該地區700 hPa以上出現變冷趨勢，且變冷趨勢一致延伸到我國江南地區。同樣，位勢高度場在蒙古以及貝加爾湖地區增高，但是在華北上空300 hPa為中心出現位勢高度場的降低。緯向風場出現同位勢高度一致的變化，在位勢高度升高中心的南北兩側分別出現了東風異常和西風異常，而經向風表現出南北一致的減弱趨勢。因此，監測該地區溫度的變化對預測華北水資源和東亞夏季風的長期趨勢具有重要的意義。

丁一匯, 胡國權. 2003. 1998年中國大洪水時期的水汽收支研究[J]. 氣象學報, 61 (2): 129–145. Ding Yihui, Hu Guoquan. 2003. A study on water vapor budget over China during the 1998 severe flood periods [J]. Acta Meteorologica Sinica (in Chinese), 61 (2): 129–145.

高橋浩一郎. 1979. 月平均溫度月降水量以及蒸發量的推定方式[J]. 天氣(日本), 26 (12): 759–763. Gaoqiaohaoyilang. 1979. Evaporation estimating method based on monthly average temperature and precipitation [J]. Weather (Japan), 26 (12): 759–763

郭其蘊, 蔡靜寧, 邵雪梅, 等. 2003. 東亞夏季風的年代際變率對中國氣候的影響[J]. 地理學報, 58 (4): 569–576. Guo Qiyun, Cai Jingning, Shao Xuemei, et al. 2003. Interdecadal variability of East?Asian summer monsoon and its impact on the climate of China [J]. Acta Geographica Sinica, 58(4): 569–576.

黃榮輝, 徐予紅, 周連童. 1999. 我國夏季降水的年代際變化及華北干旱化趨勢[J]. 高原氣象, 18 (4): 465–476. Huang Ronghui, Xu Yuhong, Zhou Liantong. 1999. The interdecadal variation of summer precipitations in China and the drought trend in North China [J]. Plateau Meteorology (in Chinese), 18 (4): 465–476.

Kalnay E, Kanamitsu M, Kistler R, et al. 1996. The NCEP/NCAR 40-year reanalysis project [J]. Bull. Amer. Meteor. Soc., 77 (3): 437–472.

李建, 宇如聰, 陳昊明, 等. 2010. 對三套再分析資料中國大陸地區夏季降水量的評估分析[J]. 氣象, 36 (12): 1–9. Li Jian, Yu Rucong, Chen Haoming, et al. 2010. Evaluation and analyses of summer rainfall over mainland China in three reanalysis datasets [J]. Meteorological Monthly (in Chinese), 36 (12): 1–9.

馬柱國, 符淙斌. 2006. 1951～2004年中國北方干旱化的基本事實[J]. 科學通報, 51 (20): 2429–2439. Ma Zhuguo, Fu Congbin. 2006. 1951–2004 basic fact of drought in northern China [J]. Chinese Science Bulletin (in Chinese), 51 (20): 2429–2439.

任國玉, 姜彤, 李維京, 等. 2008. 氣候變化對中國水資源情勢影響綜合分析[J]. 水科學進展, 19 (6): 772–779. Ren Guoyu, Jiang Tong, Li Weijing, et al. 2008. An integrated assessment of climate change impacts on China’s water resources [J]. Advanced in Water Science (in Chinese), 19 (6): 772–779.

宋正山, 楊輝, 張慶云. 1999. 華北地區水資源各分量的時空變化特征[J]. 高原氣象, 18 (4): 552–566. Song Zhengshan, Yang Hui, Zhang Qingyun. 1999. The characteristics of temporal and spatial variation of water resource components over North China [J]. Plateau Meteorology (in Chinese), 18 (4): 552–566.

田紅, 郭品文, 陸維松. 2002. 夏季水汽輸送特征及其與中國降水異常的關系[J]. 南京氣象學院學報, 25 (4): 498–502. Tian Hong, Guo Pingwen, Lu Weisong. 2002. Features of water vapor transfer by summer monsoon and their relations to rainfall anomalies over China [J]. Journal of Nanjing Institute of Meteorology (in Chinese), 25 (4): 498–502.

Uppala S M, K?llberg P W, Simmons A J, et al. 2005. The ERA-40 re-analysis [J]. Quart. J. Roy. Meteor. Soc., 131 (612): 2961–3012.

Wang H J. 2001. The weakening of the Asian monsoon circulation after the end of 1970’s [J]. Adv. Atmos. Sci., 18: 376–386.

徐桂玉, 楊修群, 孫旭光. 2005. 華北降水年代際、年際變化特征與北半球大氣環流的聯系[J]. 地球物理學報, 48 (3): 511–518. Xu Guiyu, Yang Xiuqun, Sun Xuguang. 2005. Interdecadal and interannual variation characteristics of rainfall in North China and its relation with the Northern Hemisphere atmospheric circulations [J]. Chinese J. Geophys. (in Chinese), 48 (3): 511–518.

徐康, 何金海, 祝從文. 2011a. 近50年中國東部夏季降水與貝加爾湖地表氣溫年代際變化的關系[J]. 氣象學報, 69 (4): 570–580. Xu Kang, He Jinhai, Zhu Congwen. 2011a. The interdecadallinkageof the summer precipitation in eastern China with the surface air temperature over Lake Baikal in the past 50 years [J]. Acta Meteor. Sinica (in Chinese), 69 (4): 570–580.

徐康, 祝從文, 何金海. 2011b. 近50年環貝加爾湖區變暖對中國華北夏季降水的影響機理[J]. 高原氣象, 30 (2): 309–317. Xu Kang, Zhu Congwen, He Jinhai. 2011b. Impact of the surface air temperature warming around Lake Baikal on trend of summer precipitation in North China in the past 50 years [J]. Plateau Meteorology (in Chinese), 30 (2): 309–317.

葉柏生, 楊大慶, 丁永建, 等. 2004. 我國過去50 a來降水變化趨勢及其對水資源的影響(Ⅰ): 年系列[J]. 冰川凍土, 26 (5): 587–594. Ye Bosheng, Yang Daqing, Ding Yongjian, et al. 2004. Variation trend of precipitation and its impact on water resources in China during last 50 years (Ⅰ): Annual variation [J]. Journal of Glaciology and Geocryology (in Chinese), 26 (5): 587–594.

宇如聰, 周天軍, 李建, 等. 2008. 中國東部氣候年代際變化三維特征的研究進展[J]. 大氣科學, 32 (4): 893–905. Yu Rucong, Zhou Tianjun, Li Jian, et al. 2008. Progress in the studies of three-dimensional structure of interdecadal climate change over eastern China [J]. Chinese Journal of Atmospheric Sciences (in Chinese), 32 (4): 893–905.

張慶云. 1999. 1880年以來華北降水及水資源的變化[J]. 高原氣象, 18 (4): 486–495. Zhang Qingyun. 1999. The variations of the precipitation and water resources in North China since 1880 [J]. Plateau Meteorology (in Chinese), 18 (4): 486–495.

Zhou Tianjun, Yu Rucong. 2005. Atmospheric water vapor transport associated with typical anomalous summer rainfall patterns in China [J]. Journal of Geophysical Research, 110: D08104, doi: 10.1029/ 2004JD005413.

周連童, 黃榮輝. 2003. 關于我國夏季氣候年代際變化特征及其可能成因的研究[J]. 氣候與環境研究, 8 (3): 274–290. Zhou Liantong, Huang Ronghui. 2003. Research on the characteristics of interdecadal variability of summer climate in China and its possible cause [J]. Climate and Environmental Research (in Chinese), 8 (3): 274–290.

Zhu C W, Wang B, Qian W H, et al. 2012. Recent weakening of northern East Asian summer monsoon: A possible response to global warming [J]. Geophys. Res. Lett., 39: L09701, doi: 10.1029/2012GL051155.

張書萍, 祝從文, 周秀驥. 2014. 華北水資源年代際變化及其與全球變暖之間的關聯[J]. 大氣科學, 38 (5): 1005?1016, doi:10.3878/j.issn.1006-9895. 1401.13231. Zhang Shuping, Zhu Congwen, Zhou Xiuji. 2014. Decadal variability of water resources in North China and its linkage to the global warming [J]. Chinese Journal of Atmospheric Sciences (in Chinese), 38 (5): 1005?1016.

Decadal Variability of Water Resources in North China and Its Linkage to Global Warming

ZHANG Shuping, ZHU Congwen, and ZHOU Xiuji

,100081

On the basis of the water budget equation, we discuss the decadal variability of water resources and available precipitation in North China, and address its link to boreal large-scale circulation and global warming between 1951 and 2011. The data sets utilized include monthly station-observed surface air temperature and precipitation, sounding-observed air temperature and geopotential height, and NCEP/NCAR and ERA40 reanalyzed atmospheric components. Our results suggest that annual average evaporation accounts for 69% of the total water resource in contrast to 31% by the available precipitation in North China. Most of North China’s water resources come from the East Asian summer monsoon (EASM) season. The water resources in this region exhibit an obvious decadal decrease around the 1980s, and drops further in the first decade of 21st century with the increase of local evaporation. Analysis suggests that westerly winds and the EASM jointly supply the water resources in North China, and the lower and higher geopotential height (GPH) around Lake Baikal and the western North Pacific may enhance water resources in this region. This is clearly supported by the NCEP/NCAR and ERA40 reanalyzed data sets, but also by the atmospheric sounding over East Asia. Evidence indicates that the global warming-related low-level warming that caused an anomalous anti-cyclone around Lake Baikal resulted in a decrease in water resources in North China. Therefore monitoring the surface air temperature around Lake Baikal may provide a good indicator for long-term forecasting of water resources and the EASM.

Water resource in North China, Available precipitation, Decadal variability, Global warming

1006?9895(2014)05?1005?12

P467

A

10.3878/j.issn.1006-9895.1401.13231

2013?07?30，2014?01?24 收修定稿

國家自然科學基金項目41221064，中國氣象科學研究院基本科研業務項目Z2010Z001、Z2010Z003、2013Z002，中國科學院碳專項XDA05090408

張書萍，女，1986年出生，碩士研究生，主要從事氣候變化及東亞季風研究。E-mail: zsp_1986@126.com

祝從文，E-mail: tomzhu@cams.cma.gov.cn