三峽水庫對宜昌地區天氣影響的數值模擬

魚 艇， 朱克云， 張 杰， 潘蔚琳， 丁潔瑩

(1.成都信息工程大學大氣科學學院高原大氣與環境四川省重點實驗室，四川成都610225;2.成都軍區空軍氣象中心，四川 成都610041;3.中國科學院大氣物理研究所，北京100029;4.荷蘭皇家氣象局，荷蘭 De Bilt 3731GA)

0 引言

三峽水庫位于湖北省宜昌市長江支流上，是目前世界上最大的水力發電站。水庫寬約2.3 km，平均寬度約1 km，長約600 km，水域面積約1040 km2。三峽水庫修建在兩山之間，形成了明顯的山地氣候區，湖泊效應與山地氣候區結合，導致周邊區域天氣變得更加復雜。

三峽水庫的建立引起國內外學者的廣泛關注。水庫是否會影響周邊氣候，甚至全球氣候一直是熱點科學問題。對影響范圍的研究有:王國慶、馬占山等［1-2］發現，三峽水庫主要影響當地氣溫，風速，實際的蒸發量等，并且影響的主要范圍是河流兩側約10 km。李艷等［3］則使用WRF模式研究水庫對當地氣候的影響，其結果表明水庫主要影響的是近地面層，并且影響效果隨高度的上升迅速下降，對850 hPa及以上高度層幾乎沒有影響。對降水的研究有:Wu等［4］利用美國宇航局TRMM衛星數據和數值模式MM5模擬降水率，結果表明三峽水庫的建立引起地表參數的變化，導致在水庫上游地區降水量增加，在水庫附近的地區降水量減少。王中等［5］使用WRF模式模擬的結果表明，三峽水庫建立后，明顯對當地的降水有影響，并使得該地區的風速增大。而對于溫度效應的研究，陳鮮艷等［6］系統地分析了三峽庫區局地氣候變化的特征，研究結果表明20世紀90年代后三峽庫區氣溫有顯著上升趨勢，蓄水后受水域擴大影響近庫地區的氣溫發生了一定變化，表現出冬季增溫效應，夏季有弱降溫效應，但總體以增溫為主。張強、毛以偉等［7-8］根據三峽水庫壩區周邊氣候考察氣象站的氣溫觀測資料客觀分析了氣溫隨時間和地形的變化規律。王浩、段德寅等［9-12］則使用數值模擬方法研究了三峽水庫的溫度效應，認為水庫水體白天有降溫效應，夜間有升溫效應，這種效應晴天大于陰天，夏季夜間的升溫效應小于白天的降溫效應，冬季則相反。Miller等［13］也利用中尺度區域氣候模式MM5，針對1990年3月2日到5月16日共44天的無雨日，對三峽水庫周邊局地氣候進行了數值模擬，研究表明三峽水庫的建設導致平均地面溫度下降了2.9℃，降水沒有發生明顯的變化;但是Xiao等［14］通過經驗正交函數(EOF)的方法研究得出相反的結果，他們的結果表明該地區的改變是自然的變化結果，而不是水庫所造成的影響，這值得對其進行更加詳細的分析研究。

以往研究大都是利用觀測資料進行客觀統計分析或者單一氣候要素的局地氣候效應影響，而沒有真正涉及到詳細的氣象要素，因此文中主要研究水庫建成后對宜昌地區溫度，濕度，風等氣象要素的影響，并通過CLASS模式對溫度、濕度平流影響作用進行驗證。

1 研究方法

使用WRF模式模擬水庫及周邊地區3天的天氣演變情況，然后用實際觀測數據與模式數據進行對比驗證。模式結果可以給出湖泊對宜昌地區單個氣象要素的影響及差異，以及湖泊形狀、面積對宜昌地區氣象要素的影響。最后通過CLASS模式對平流作用的影響進行驗證。

1．1WRF 模式

WRF模式是由美國國家大氣研究中心(NCAR)開發的一個完全可壓縮以及非靜力模式。模式具有多層嵌套，以及大量的物理參數。

文中所使用的資料為ECMWF再分析資料，且把該資料作為WRF模式的初始數據。由于美國懷俄明州立大學網站有近幾十年(北京時08點和20點)最近的探空數據文件，而且資料相對及時完善，因此使用該大學探空資料站所探測的宜昌(30.70°N，111.30°E)站的高空探測資料，以及NOAA衛星的地表逐小時資料來驗證WRF模式結果。

主要通過改變地表參數創建湖泊來模擬對宜昌地區天氣變化的影響。所選用的時間段為2011年5月12日00:00UTC至15日00:00UTC，即2011年5月12日08:00(北京時)至15日08:00(北京時)，時間步長為6 h。該時段該地區受高空穩定的大氣環流和地面高壓系統作用，天氣晴朗，天氣系統穩定。

使用WRF模式建立3個試驗模擬當地的天氣變化情況:(1)初始試驗，假設該地區沒有湖泊。(2)試驗1，用正方形湖泊近似模擬大壩建立之后的河流湖泊，湖泊面積為1089 km2，長寬均為33 km。(3)試驗2，用矩形湖泊近似模擬未建立大壩之前的河流湖泊，湖泊面積為1215 km2，全長90 km，東半部分寬為9 km，西半部分寬為18 km。宜昌市經緯度分別為30.70°N，111.30°E，三峽水庫正好位于城市的西部。試驗以宜昌為中心共設置了3個區域，水平分辨率分別為27、9、3 km，垂直方向分為31層。所有的區域有31×31個網格，整個模式模擬的區域約656100平方公里。試驗1和試驗2分別是改變地表參數為森林或者草地的湖泊區域(圖1)。

模式中所使用的物理方案分別為:WSM 3-class單冰方案，長波輻射rrtmg方案，短波輻射Dudhia方案，Monin-Obukhov近地面層方案，Unified Noah陸面模式層選方案，YSU邊界層層選方案，Kain Fritsch積云選擇方案。

圖1 區域中陸面情況

1．2 CLASS 模式

CLASS模式是一個通過測試系統中各種物理和化學過程來檢驗系統參數的盒子模式。包含7個部分:基礎條件，風條件，輻射/地形條件，地表條件，化學條件，反射條件，高級曲面條件。由于受觀測條件的限制，因此在文中只利用CLASS模式的基本條件來分析平流在系統中的作用。計算溫度和比濕平流的公式為，φ表示環境中的溫度和比濕。

模式利用2個試驗分析宜昌地區邊界層高度變化和三峽大壩附近地區相對濕度變化情況。因為只考慮基本條件，所以將該地區這一時期的平均值改為初始的溫度和濕度。模式共假設4個試驗:沒有任何形式的平流;只含有溫度平流;只含有濕度平流;包含有溫度平流和濕度平流。

2 結果與討論

2．1 觀測值對比結果

2.1.1 觀測值與模擬值對比結果

圖2 宜昌市從5月12日01:00UTC到5月15日00:00UTC 3種模式試驗與觀測值的對比結果

如圖2(a)所示，模式模擬的溫度值的整體變化趨勢與觀測溫度值幾乎一致。而且5月12日01:00UTC至5月13日00:00UTC，3個模式模擬的溫度值與觀測溫度值近似一致。但是從5月13日00:00UTC開始出現變化，起初，觀測的溫度值比模式的溫度值高約1.5℃，但當太陽升起之后，兩者的數值又近似相同。太陽落山后，模式的溫度卻又比觀測溫度高3℃。5月14日00:00UTC至5月15日00:00UTC，白天的溫度模擬值比觀測溫度值高1℃，夜間的模擬值也比觀測溫度值高2℃。溫度的回歸系數為0.87。

如圖2(b)所示，模式中比濕值與觀測值的對比結果卻出現很大的變化差異。在白天比濕值的變化趨勢與觀測值能夠近似保持一致，時間段分別為0～12 h，27～33 h和60～69 h。但是在夜間它們卻完全不同，模式的結果明顯低于觀測值且比濕的回歸系數僅為0.1。

如圖2(c)和(d)所示，模式中風向的角度值與觀測值存在大小差異，而且風速模擬值明顯高于觀測值，但是它們整體的變化趨勢與觀測值卻近似一致。宜昌位于山谷中，白天風以偏西風為主，從庫區吹向城市。而夜間卻以東南風為主，由城市吹向庫區。這也與Stivari等［15］研究伊泰普水利工程對局地環流的影響所發現的結果一致，即湖面風白天向四周發散，而晚上則向湖面集中。

3.1.2 溫度和比濕的垂直分布情況

圖3為5月13日12:00UTC的溫度和比濕隨高度的變化情況，從圖中可發現模式中近地面層的溫度隨高度變化趨勢與觀測值趨于一致，且明顯高于觀測值，但模式中比濕的結果與觀測值完全不一致，其主要原因可能是比濕受許多參數的影響，如土壤中水分含量和邊界層高度，如果模式中這些參數是錯誤的，比濕的誤差也就可能被放大。在WRF模式中所產生的比濕結果不足以分析對城市的影響，所以將使用CLASS模式來驗證溫度平流和濕度平流對城市影響。

圖3 5月13日12:00UTC溫度和比濕的垂直分布圖

2．2 數值模式對比結果

2.2.1 3種試驗對比結果

如圖2(a)所示，通過比較初始試驗、試驗1和試驗2的溫度值發現，在有湖泊的試驗1和試驗2中，宜昌地區的夜間溫度值比初始試驗要高一點，白天的溫度則相對低一點。試驗2的模擬結果介于初始試驗和試驗1之間，因此可以得出:(1)建成水庫后的湖泊對宜昌地區溫度的影響作用比沒有湖泊(沒有建造水庫)的影響要大;(2)湖泊總面積的大小對宜昌地區溫度的影響不明顯，而與臨近城市的區域湖泊面積大小有較大關系。如圖2(b)所示，在所模擬的3天中，比濕結果具有不確定性，試驗1和試驗2的結果都顯示:在白天的比濕結果比初始試驗的值大，而在夜間卻比初始試驗的值小。如圖2(c)和(d)所示，試驗1對風向、風速的影響都大于其它兩個試驗，而初始試驗和試驗2的模擬結果卻比較相近，這說明臨近城市的區域湖泊的面積大小對宜昌地區的影響作用更大。

2.2.2 試驗區域溫度，風速，氣壓的分布情況

圖4 夜間和白天3個試驗區域中溫度，風速，氣壓的分布情況

圖4顯示了模式模擬區域中的溫度，風速，氣壓的分布情況。分別比較6個圖發現，在白天試驗1和試驗2湖面區域的溫度明顯低于初始試驗，而在夜間試驗1和試驗2中溫度卻比初始試驗高，并且試驗1湖面區域的溫度最高，試驗2次之，初始試驗溫度最低。這表明由于湖泊的作用，該區域的夜間溫度值比初始試驗要高一點，白天的溫度卻相對較低。從圖中還可以發現，大氣環流中的風場也發生了變化。湖面風不再是簡單白天向四周發散，晚上向湖面集中。而是在湖面上空形成了一個小的環流系統。在夜間，湖面的溫度相對周邊地區較高，引起湖面上空的氣壓場發生變化，湖面氣壓降低，導致在湖面上空形成一個小的局地氣旋環流系統，氣流向氣旋中心輻合。而在白天，湖面的溫度比周邊地區低，使得湖面上空的氣壓升高，導致在湖面上空形成了一個小的局地反氣旋環流系統，氣流向外輻散。這也就有利的驗證了宜昌地區風向在白天和夜間不同的原因。

2.2.3 有湖試驗溫度、相對濕度的分布情況

如圖5所示，試驗1在白天的溫度明顯低于試驗2，而在夜間試驗1的溫度卻高于試驗2。這也就充分說明了湖泊總面積的大小對宜昌地區溫度的影響不明顯，而與臨近城市區域湖泊面積大小有較大關系。

圖5 白天和夜間有湖泊試驗溫度分布圖

在白天，受太陽輻射加熱蒸發作用，湖面上空的相對濕度逐漸升高，由于風的作用，周邊區域相對濕度也隨之升高。夜間相對濕度的變化情況，如圖6所示，試驗中水汽是由東北向西南輸送，夜間湖面溫度高，高的熱通量的暖濕水汽使得湖面西南區域地區的相對濕度有明顯的降低。而在宜昌地區水汽未經過湖面區域，溫度較低，該區域的相對濕度的變化不大。因此，湖面周圍地區相對濕度的變化可能原因是溫度、濕度平流作用。

圖6 夜間有湖泊試驗相對濕度的分布圖

2．3 CLASS 模式結果

2.3.1 溫度、濕度平流作用情況

圖7 CLASS模式中平流作用的影響

利用CLASS模式分別計算試驗1中溫度平流和濕度平流，其結果是溫度平流值為0.0003 k/s，濕度平流值為0.0001 g/(kg.s)。如圖7所示，單一濕度平流對相對濕度的影響效果最大，沒有平流作用的影響效果次之，兩者共同作用的影響效果的再次之，單一溫度平流的影響效果最小。因此，湖面周圍地區相對濕度的變化，主要取決于濕度平流作用。

2.3.2 邊界層高度變化情況

溫度和比濕對邊界層高度的影響非常重要，圖8(a)顯示了在這3天中3個模擬試驗中邊界層高度隨時間的變化情況，從圖中可以看出3個試驗中邊界層高度在第一天和第三天近似是相同的，然而試驗1的邊界層高度在第二天明顯有所降低，可能原因是濕冷空氣從湖面地區輸送到城市。

圖8 邊界層高度變化

利用CLASS模式對第二天的平流作用進行驗證，分別計算了第二天的溫度平流和濕度平流，其結果是溫度平流為 -0.0002 K/s，濕度平流為0.00007 g/(kg.s)。從圖8(b)中可以看出，溫度平流和濕度平流都使得邊界層的高度低于初始試驗，而且二者共同作用使得邊界層高度最低。單一濕度平流影響使得邊界層高度比初始試驗略有降低，而單一溫度平流影響卻使得邊界層高度降低了很多。這說明溫度平流對邊界層高度的變化起主導作用，影響效果最明顯，濕度平流對邊界層高度的變化影響較小，二者共同作用的影響效果是最強的。

3 結論

主要利用WRF模式來模擬湖泊區域3天的天氣演變情況，并利用CLASS模式對其平流作用進行了驗證。得到以下結論:

(1)宜昌的溫度在夜間有所增加，白天則相反。對宜昌地區溫度影響最大的是正方形湖面，這說明宜昌地區的溫度主要受臨近城市區域湖泊面積大小的影響而非湖泊總面積大小。

(2)由于三峽水庫的建立，局地的大氣環流也發生了變化。湖面風并非是簡單的白天向四周發散，晚上向湖面集中，而是在湖面上空形成了一個小的環流系統。

(3)湖面周邊區域相對濕度的變化主要取決于濕度平流作用。

(4)溫度平流對邊界層高度的變化起主導作用，影響效果最明顯，濕度平流對邊界層高度的變化影響較小，二者共同作用的影響效果是最強的。

主要通過模式創建湖泊模型模擬水庫對宜昌地區天氣變化的影響，對宜昌溫度變化做出一定的分析，并提出了局地環流系統和對溫濕平流的作用。但是由于三峽地區地形地貌的復雜，以及考慮到模式本身分辨率、敏感性的限制和物理過程參數化方案的不同，所以還需要更多的研究進行進一步的完善和驗證。

［1］ 王國慶，張建云，賀瑞敏，等.三峽水利工程對區域氣候影響的初步分析，自主創新與持續增長第十一屆中國科協年會(2009).

［2］ 馬占山，張強，秦琰琰.三峽水庫對區域氣候影響的數值模擬分析［J］.長江流域資源和環境，2010，19(9):1045-1053.

［3］ 李艷，高陽華，陳鮮艷，等.三峽下墊面變化對區域氣候效應的影響研究［J］.南京大學學報:自然科學版，2011，47(3):330-38.

［4］ Wu L G，Q Zhang.Three Gorges Dam affects regional precipitation［J］.Geophysical Research Letters，2006，33(13).

［5］ 王中，杜欽，白瑩瑩.三峽下墊面變化對重慶氣象要素影響的數值模擬［J］.西南大學學報:自然科學版，2012，34，(3):102-109.

［6］ 陳鮮艷，張強，葉殿秀，等.三峽庫區局地氣候變化.長江流域資源與環境［J］.2009，18(1):47-51.

［7］ 張強，萬素琴，毛以偉，等.三峽庫區復雜地形下的氣溫變化特征［J］.氣候研究進展，2005，1(4):164-167.

［8］ 毛以偉，陳正洪，王玨，等.三峽水庫壩區蓄水前水體對水庫周邊氣溫的影響［J］.氣象科技，2008，33(4):334-339.

［9］ 王浩，傅抱璞.水體的溫度效應［J］.氣象科學，1991，11(3):233-243.

［10］ 王浩.深淺水體不同氣候效應的初步研究［J］.南京大學學報:自然科學版，1993，29(3):517-522.

［11］ 王浩.陸地水體對氣候影響的數值研究［J］.海洋與湖沼，1991，22(5):467-473.

［12］ 段德寅，傅抱璞，王浩，等.三峽工程對氣候的影響及其對策［J］.湖南師范大學學報:自然科學版，1996，19(1):87-92.

［13］ Miller N L，Jin J，T sang C F.Local climate sensitivity of the Three Gorges Dam［J］.Geophysical Research Letters，2005，32.

［14］ Xiao C，R C Yu.Precipitation characteristics in the Three Gorges Dam vicinity［J］.International Journal of Climatology，2010，30(13):2021-2024.

［15］ Stivari，S?nia M S，Amauri P de Oliveira，Hugo A Karam，et al.2003:Patterns of Local Circulation in the Itaipu Lake Area:Numerical Simulations of Lake Breeze［J］.J.Appl.Meteor.，2003，42:37-50.