三北防護(hù)林工程區(qū)氣候變化分析

王鵬濤， 延軍平， 蔣 沖， 劉憲鋒, 廖光明

(1.陜西師范大學(xué) 旅游與環(huán)境學(xué)院， 陜西 西安 710119； 2.北京師范大學(xué) 全球變化與地球科學(xué)研究院 地表過程與資源生態(tài)國(guó)家重點(diǎn)實(shí)驗(yàn)室， 北京 100875； 3.北京師范大學(xué) 資源學(xué)院， 北京 100875)

1978年11月，為改善中國(guó)北方的生態(tài)環(huán)境，防止三北地區(qū)風(fēng)沙災(zāi)害和水土流失的惡化，正式啟動(dòng)三北防護(hù)林工程，工程建設(shè)35 a以來取得了一定的生態(tài)效益[1]。20世紀(jì)以來全球進(jìn)入變暖期[2-5]，而防護(hù)林的建設(shè)期，正處在近百年全球變暖加劇期[4-6]。中國(guó)對(duì)氣候變暖有著顯著的響應(yīng)[4-5]，但這種響應(yīng)存在區(qū)域差異[7]，三北防護(hù)林建設(shè)以來，工程區(qū)氣候?qū)θ蜃兓憫?yīng)的時(shí)間、幅度如何，這種變化是否會(huì)阻礙防護(hù)林生態(tài)效益的發(fā)揮都是亟需探討的問題。

多年來，針對(duì)全國(guó)氣候變化及北方地區(qū)的局部氣候變化的研究較多[4-15]，可作為研究區(qū)域氣候背景的科學(xué)依據(jù)，但未充分揭示出三北防護(hù)林工程區(qū)的氣候變化特點(diǎn)。而關(guān)于三北防護(hù)林的研究，又多集中在防護(hù)林建設(shè)給區(qū)域帶來的生態(tài)效應(yīng)[16-21]，對(duì)該區(qū)多年來的氣候變化的響應(yīng)研究較少。趙國(guó)藏[21]選用氣溫及相對(duì)濕度的指標(biāo)，分析了1951—1988年工程實(shí)施引起的區(qū)域氣候變化。但是研究序列過短，對(duì)氣候變化的空間差異也未作分析。王強(qiáng)[22]在分析植被覆被變化對(duì)氣候變化的響應(yīng)時(shí)，采用1982—2006年氣象數(shù)據(jù)，對(duì)工程區(qū)25 a的氣溫降水變化作了研究，可是該研究也存在時(shí)間序列過短的問題，對(duì)于防護(hù)林建設(shè)前后氣候變化的階段特征缺乏對(duì)比。

本文選取117個(gè)站點(diǎn)52 a的氣象數(shù)據(jù)，對(duì)三北防護(hù)林多年氣候變化進(jìn)行分析，探討該區(qū)氣候?qū)θ蜃兓捻憫?yīng)及其對(duì)三北防護(hù)林工程建設(shè)的可能影響，以期為防護(hù)林后續(xù)建設(shè)布局及改善區(qū)域生態(tài)環(huán)境等提供科學(xué)依據(jù)。

1 研究區(qū)概況

三北防護(hù)林工程區(qū)位于中國(guó)北半部，東起黑龍江省的賓縣，西至新疆維吾爾自治區(qū)的烏孜別克里山，北抵國(guó)界線，南沿天津、汾河、渭河、洮河下游、布爾汗達(dá)山、喀喇昆侖山，東西長(zhǎng)4 480 km，南北寬560～1 460 km，面積約4.069×106km2，占國(guó)土總面積的42.4%。建設(shè)范圍包括黑龍江、吉林、遼寧、北京、天津、河北、山西、內(nèi)蒙古、陜西、甘肅、寧夏、青海、新疆13個(gè)省(自治區(qū)、直轄市)的551個(gè)縣(市、旗、區(qū))[23]。

三北地區(qū)地勢(shì)西高東低，呈3級(jí)臺(tái)階分布。全區(qū)地貌類型多樣，山地，丘陵，平原，高原盆地，沙漠均有分布。根據(jù)區(qū)域自然特征差異，三北防護(hù)林工程區(qū)在地域上區(qū)劃為東北西部，蒙新，黃土高原，華北北部4個(gè)Ⅰ級(jí)區(qū)，22個(gè)Ⅱ級(jí)區(qū)，59個(gè)Ⅲ級(jí)區(qū)。該區(qū)氣候類型屬北溫帶大陸性季風(fēng)氣候，就全國(guó)范圍而言屬少雨區(qū)，氣候干旱，自然植被以森林草原，草原，荒漠草原為主，生態(tài)環(huán)境脆弱[24-26],再加上三北地區(qū)土地利用結(jié)構(gòu)的不合理，干旱、半干旱地區(qū)沙漠化逐年加劇，導(dǎo)致生態(tài)平衡嚴(yán)重失調(diào)和經(jīng)濟(jì)的惡性循環(huán)。

2 資料與方法

2.1 資料選取

為反映各站點(diǎn)在防護(hù)林建設(shè)前后的氣候變化，參考在一期工程規(guī)劃的縣域，選取其中代表性站點(diǎn)117個(gè)，并從氣象數(shù)據(jù)完整的1960年開始，統(tǒng)計(jì)到2011年，建立氣象資料的時(shí)間序列。劃分1960—1977年建設(shè)前為第一時(shí)段；1978年(一期工程開始)—1995年(二期工程結(jié)束)命名為第二時(shí)段，此時(shí)一期117個(gè)站點(diǎn)防護(hù)林生態(tài)效應(yīng)已經(jīng)發(fā)揮出來；1996—2011年為第三時(shí)段，以對(duì)比分析近期區(qū)域氣候變化。三個(gè)時(shí)段跨度分別為18,18,16 a，序列長(zhǎng)短接近，有利于進(jìn)行數(shù)據(jù)的比較。站點(diǎn)的氣溫、降水資料，來源于中國(guó)氣象資源共享網(wǎng)。

2.2 研究方法

本文采用氣候傾向率，距平分析法及空間插值等方法，分析該區(qū)的近52 a的氣候變化趨勢(shì)。統(tǒng)計(jì)各站點(diǎn)不同時(shí)段的氣象要素距平值(均以52 a均值為距平基準(zhǔn)值)及氣候傾向率，利用ArcGIS軟件繪制出各時(shí)段的距平分布圖及多年傾向率分布圖。為了直觀反映氣象數(shù)據(jù)的正負(fù)距平及正負(fù)傾向率變化，我們把零值線在相應(yīng)的分圖中標(biāo)注出來。

3 氣溫變化趨勢(shì)分析

3.1 氣溫的時(shí)間變化

該區(qū)近52 a來平均氣溫為6.81 ℃，年均溫的三個(gè)高值均出現(xiàn)在1996—2011年，分別為1998年(8.10 ℃)，2006年(7.82 ℃)，2007年(8.29 ℃)，3個(gè)低值均出現(xiàn)在1960—1977年時(shí)間段，分別為1967年(5.78 ℃)，1969年(5.50 ℃)，1976年(5.79 ℃)。其中最高溫出現(xiàn)在1998年，這與李慶偉等[9]研究結(jié)果一致，最低溫出現(xiàn)在1969年，與最高溫相差2.79 ℃。多年來區(qū)域氣溫呈明顯上升趨勢(shì)，由1960年的年均溫6.29 ℃上升到2011年的7.19 ℃，線性擬合增溫速率為0.346 ℃/10 a，符合全球氣溫變化趨勢(shì),但區(qū)域增溫速率明顯高于同期全球增溫速率 0.13±0.03 ℃/10 a[4,10]和全國(guó)增溫速率0.25 ℃/10 a[8]。

從年平均距平圖(圖1)來看，在1960—1977年，區(qū)域溫度先下降，1970年開始緩慢上升。各年氣溫，除1975年外，均處于負(fù)距平狀態(tài)，階段距平為-0.57 ℃，氣溫偏低。1978—1995年氣溫正距平年份達(dá)39%，階段距平為-0.09 ℃，已經(jīng)接近52 a平均值。 1996—2011年氣溫繼續(xù)上升，在2004年達(dá)到峰值，階段距平上升到0.74 ℃，增暖趨勢(shì)顯著。統(tǒng)計(jì)年均溫3個(gè)高值均出現(xiàn)在該段時(shí)間，也與此階段溫度升高的趨勢(shì)相符。

圖1 1960-2011年三北防護(hù)林工程區(qū)氣溫變化

3.2 氣溫的空間變化

3.2.1 氣溫距平空間變化 由圖2可見，1960—1977年工程區(qū)相對(duì)溫度普遍偏低，大多站點(diǎn)為負(fù)距平?？臻g分布上，氣溫距平的分布存在南北差異，以工程區(qū)東部的南北差異最為明顯。北部除吐魯番盆地、河西走廊外，黑龍江松嫩平原、吉林西部平原、內(nèi)蒙古高原、冀西北山地、新疆阿勒泰地區(qū)、柴達(dá)木地區(qū)為顯著低值區(qū)。遼寧北部、冀北山地、山西西部黃土高原、陜北地區(qū)、隴中地區(qū)、青海東部、新疆吐魯番盆地、天山南部地區(qū)及塔里木盆地為高值分布區(qū)，距平值在-0.56 ℃以上。3個(gè)正距平站河曲、庫(kù)車、承德均分布在該范圍內(nèi)。

1978—1995年區(qū)域大部分地區(qū)仍為負(fù)距平區(qū)，但是氣溫正距平站點(diǎn)明顯增加，各站點(diǎn)的溫度顯著提高，平均距平值上升0.48 ℃。區(qū)域氣溫正距平站點(diǎn)32個(gè)，其中30個(gè)都分布在東北、華北和北疆地區(qū)。新疆阿勒泰地區(qū)、黑龍江松嫩平原、遼西走廊、冀北山地為氣溫正距平區(qū)。距平高值中心已經(jīng)發(fā)生了遷移，相較于第一階段氣溫分布，陜西、山西、甘肅、青海東部、新疆南部由高值中心轉(zhuǎn)為低值中心，絕對(duì)溫度降低較大。

1996—2011年氣溫正距平區(qū)進(jìn)一步擴(kuò)大，僅西寧，承德，庫(kù)車3站為負(fù)距平，區(qū)域整體距平值達(dá)到正距平0.74 ℃，提高了0.85 ℃，增溫趨勢(shì)相對(duì)第二階段更為明顯，相對(duì)氣溫的高溫中心繼續(xù)遷移，內(nèi)蒙古及柴達(dá)木盆地增溫最顯著，成為新的高值中心，距平值可以達(dá)到0.78 ℃以上。塔里木盆地、華北平原的相對(duì)溫度較低。

圖2 三個(gè)時(shí)段氣溫距平及其變化趨勢(shì)分布

3.2.2 氣溫傾向率分布 統(tǒng)計(jì)各站點(diǎn)氣溫傾向值，除承德、河曲、庫(kù)車3站為負(fù)值，其他站均為正值，區(qū)域整體變暖趨勢(shì)明顯。整個(gè)工程區(qū)，除了南部的柴達(dá)木盆地為增溫高值區(qū)外，大致反映出南低北高的線性變化趨勢(shì)。內(nèi)蒙古高原、黑龍江松嫩平原、吉林西部平原、柴達(dá)木盆地均為顯著增暖區(qū)，華北北部、山西西部、陜北、隴東、塔里木盆地為增溫低值區(qū)，這與氣溫距平圖的反映出的高低值中心遷移變化過程相符。

另外，觀察第一階段距平分布圖2a與氣溫傾向圖2d，1960—1977年氣溫在空間上，呈現(xiàn)出“高值增溫慢，低值增溫快”的變化規(guī)律，即建設(shè)初期的相對(duì)氣溫低值區(qū)增溫速率最大，高值區(qū)域增溫速率最小的對(duì)應(yīng)變化趨勢(shì)，其中河北、遼寧、陜西、山西、甘肅、新疆等建設(shè)初期的高值中心，多年來增溫趨勢(shì)最為緩慢。這與郭志梅等[14]的研究結(jié)果基本一致，但是從圖2d分析，第三階段西北西部的增溫趨勢(shì)已明顯放緩。

4 降水變化趨勢(shì)分析

4.1 降水的時(shí)間變化

區(qū)域多年平均降水309.0 mm，降水的高值在3個(gè)時(shí)段均有分布(圖3)， 3個(gè)年份分別為1964年(387.8 mm)，1990年(367.7 mm)，1998年(392.2 mm)，這與氣溫高值的分布有所差異。降水3個(gè)低值則分布在前兩個(gè)時(shí)段內(nèi)，分別為1965年(242.2 mm)，1980年(253.4 mm)，1982年(246.2 mm)。最大降水量與最小降水量相差達(dá)149.9 mm。52 a來工程區(qū)年平均降水量從1960年的317.6 mm下降到2011年的280.8 mm，下降擬合速率為3.554 mm/10 a。

從降水距平的擬合曲線來看，區(qū)域平均降水總體呈現(xiàn)波動(dòng)下降趨勢(shì)。第一時(shí)段，區(qū)域降水先波動(dòng)下降，在后半段1975年附近降水量開始緩慢增加，整個(gè)時(shí)段降水距平3.42 mm。第二時(shí)段，降水量先呈升高趨勢(shì)，在90年代初期又逐漸減少，此階段降水量平均距平增加為7.84 mm，平均降水量有所增加。1996—2011年，降水變化也是呈現(xiàn)先減少，在2006年附近又波動(dòng)增加的趨勢(shì)，平均距平為-12.7 mm，此階段降水量較常年平均值偏少。

圖3 1960-2011年三北防護(hù)林工程區(qū)降水量變化

4.2 降水的空間變化

4.2.1 降水量距平空間變化 圖4為各站點(diǎn)不同時(shí)段降水量距平的空間分布。在第一時(shí)段(圖4a)，區(qū)域降水距平3.42 mm，其空間分布基本呈現(xiàn)出東高西低的變化特征，西部地區(qū)相對(duì)降水量普遍偏少。雨量負(fù)值區(qū)分布在東部的黑龍江松嫩平原和大興安嶺及內(nèi)蒙古陰山山地。東北西部其他地區(qū)、華北北部、黃土高原均為降水距平高值區(qū)，58個(gè)正距平站均分布在此區(qū)。

圖4 3個(gè)時(shí)段降水距平及其變化趨勢(shì)分布

第二時(shí)段(圖4b)，降水正距平面積進(jìn)一步擴(kuò)大，有77個(gè)站為高值站，平均距平增加為7.84 mm?？梢钥闯鲞@一時(shí)段，工程區(qū)普遍降水增加，降水空間分布仍然表現(xiàn)出東西差異，低值區(qū)主要分布在山西西部、陜西、青海及甘肅祁連山區(qū)、塔里木盆地及天山山地、其他地區(qū)降水量普遍增加，其中東部的東北西部、華北北部地區(qū)為降水高值中心，與前一時(shí)段相比，陜西、山西、甘肅祁連山地降水量明顯下降。新疆、柴達(dá)木盆地、河西走廊地區(qū)降水有所增加。

第三時(shí)段(圖4c)，降水平均距平為-12.7 mm，整體降水呈下降趨勢(shì)。降水分布的東西差異明顯，大致以河西走廊及青藏高原為界，界線以西為降正距平區(qū)。此階段37個(gè)降水正距平站，有28個(gè)都分布在界限以西。對(duì)比前一時(shí)段，可以看出，東部地區(qū)除黑龍江松嫩平原外，其他區(qū)域普遍降水有所減少，而西部區(qū)域仍然為正距平態(tài)勢(shì)，降水量進(jìn)一步增加。

4.2.2 降水傾向率分布 從降水傾向率統(tǒng)計(jì)來看(圖4d)，降水線性趨勢(shì)的東西差異明顯，全區(qū)降水呈增加趨勢(shì)的站點(diǎn)達(dá)到了59%。西部地區(qū)除庫(kù)爾勒，額濟(jì)納旗，民和站外，降水傾向率均大于0，傾向率的最大值伊寧站就位于西部，傾向率達(dá)到了21.4 mm/10 a。說明近52 a來，中國(guó)西部地區(qū)降水量顯著增加。這與葉柏生等[11]，姚惠明等[12]研究結(jié)果一致。區(qū)域東部的黑龍江松嫩平原、內(nèi)蒙古陰山山地也是顯著的增加區(qū)。東部其余地區(qū)，降水則呈現(xiàn)減少趨勢(shì)，線性趨勢(shì)圖與工程區(qū)降水高值區(qū)的空間變化結(jié)果相符，只是總的變化趨勢(shì)反映不出東部的東北西部、華北北部地區(qū)在第二時(shí)段降水增加的過程。

5 結(jié) 論

(1)近52 a來，工程區(qū)氣溫呈顯著增加趨勢(shì)，多年增溫速率為0.346 ℃/10 a。防護(hù)林建設(shè)前為偏冷期，第二時(shí)段，區(qū)域的氣溫逐步升高，第三時(shí)段為顯著偏暖期，其中1998年達(dá)到多年最高值。區(qū)域氣溫距平在三個(gè)時(shí)段表現(xiàn)出普遍低溫—高溫區(qū)增多—普遍高溫的變化規(guī)律。另外區(qū)域的高低溫區(qū)發(fā)生明顯的遷移，高溫區(qū)由黃土高原、塔里木盆地、華北平原北部遷移到第二時(shí)段的黑龍江、吉林、北疆阿勒泰地區(qū)，再到第三時(shí)段的內(nèi)蒙古及柴達(dá)木盆地。線性增溫趨勢(shì)大致表現(xiàn)出南北差異，顯著區(qū)主要集中在柴達(dá)木盆地、內(nèi)蒙古高原、黑龍江松嫩平原、吉林西部平原。另外第一時(shí)段距平分布圖與氣溫傾向圖呈現(xiàn)出低值區(qū)增溫速率最大，高值區(qū)域增溫速率最小的對(duì)應(yīng)變化趨勢(shì)。綜合氣溫各時(shí)段距平及傾向率變化，可以看出三北防護(hù)林區(qū)在全球變化的背景下，區(qū)域響應(yīng)仍是以變暖趨勢(shì)為主，但是三北防護(hù)林的工程的建設(shè)實(shí)施，可能在局部地區(qū)起到了減緩這種趨勢(shì)的作用。其中河北、吉林、遼寧、內(nèi)蒙、陜西、山西、黑龍江省，均由建設(shè)初期的升溫趨勢(shì)變?yōu)榻禍刳厔?shì)。河北、遼寧、陜西、山西、甘肅、新疆等建設(shè)初期的高溫中心增溫趨勢(shì)較為緩慢的現(xiàn)象，說明這幾個(gè)省區(qū)的三北防護(hù)林工程對(duì)其氣溫增暖起到了一定的抑制作用。

(2)多年來區(qū)域降水量呈下降趨勢(shì)，下降速率為-3.5 mm/10 a。其時(shí)段性特征建設(shè)前偏少，建設(shè)初期有所增加，但是1996年后又有明顯減少趨勢(shì)。 區(qū)域降水量距平，在三個(gè)時(shí)段發(fā)生了顯著的空間變化。第一時(shí)段東部地區(qū)普遍降水較常年平均值高。第二時(shí)段，工程區(qū)大部分地區(qū)降水成正距平，高值區(qū)仍然位于東部。第三時(shí)段工程區(qū)西部降水較東部豐沛，降水增距平區(qū)均位于西部。降水的高值區(qū)，由工程區(qū)中部的黃土高原、華北平原北部向上轉(zhuǎn)移到東北地區(qū)北部，再轉(zhuǎn)移到西部地區(qū)。各時(shí)段降水量均有正負(fù)距平站點(diǎn)，這與氣溫距平的顯著變化不同，表現(xiàn)出了降水變化的復(fù)雜性，可能是因?yàn)榕c降水受更多的因素影響有關(guān)。降水的傾向率分布，表現(xiàn)為明顯的東西差異，東部地區(qū)降水普遍減少，西部大部分地區(qū)降水均在增加，黑龍江松嫩平原、內(nèi)蒙古陰山山地也是顯著的增加區(qū)。綜合降水距平及傾向率分布，降水整體減少，而三北防護(hù)林的建設(shè)只是在局部地區(qū)起到了減緩這種趨勢(shì)的作用。

(3)綜合氣溫及降水變化，52 a來工程區(qū)總體暖干化趨勢(shì)增加，其中氣溫變化以南北分布差異為主，而降水變化的東西差異變化更為顯著。增暖顯著區(qū)主要位于工程區(qū)北部高緯區(qū)，降水增加區(qū)則主要集中在工程區(qū)西部?？梢耘袛喑觯嗄陙碓谖鞑恳员苯⒗仗┑貐^(qū)、柴達(dá)木盆地暖濕化趨勢(shì)最為明顯，西部其余地區(qū)則表現(xiàn)為微弱的暖濕化。在東部，內(nèi)蒙古陰山地區(qū)及河套平原地區(qū)、黑龍江松嫩平原有暖濕化趨勢(shì)，內(nèi)蒙古呼倫貝爾草原、大興安嶺山區(qū)、遼嫩平原均呈明顯暖干化趨勢(shì)。黃土高原、華北北部地區(qū)、吉林，遼寧則為微弱的暖干化趨勢(shì)。

(4)一個(gè)地區(qū)的氣候變化，是由自然氣候波動(dòng)和人類活動(dòng)的影響兩類因素造成的[27]。就工程區(qū)的氣候變化而言，三北防護(hù)林的建設(shè)對(duì)工程區(qū)小區(qū)域的氣候變化起到了一定的改善作用[17-19]，但并沒有改變區(qū)域長(zhǎng)期的氣候變化趨勢(shì)。多年來,東部大部呈暖干化趨勢(shì)，西部地區(qū)雖然為暖濕化趨勢(shì)，但由于西部地區(qū)年降水總量較少[11-12]，這種變化并不能改善工程區(qū)西部干旱的氣候特征。工程區(qū)整體氣候的發(fā)展趨勢(shì)，將會(huì)限制防護(hù)林工程的實(shí)施。因此，繼續(xù)深入研究工程區(qū)氣候變化規(guī)律，制定適應(yīng)氣候變化的工程建設(shè)策略，對(duì)改善區(qū)域生態(tài)環(huán)境顯得很有必要。

[參考文獻(xiàn)]

[1]姜鳳岐,于占源,曾德慧,等.氣候變化對(duì)三北防護(hù)林的影響與應(yīng)對(duì)策略[J].生態(tài)學(xué)雜志,2009,28(9):1702-1705.

[2]IPCC. Summary for Policymakers of Climate Change 2007: The Physical Science Basis Contribution of Working Group 1 to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change[M].Cambridge: Cambridge University Press， 2007．

[3]Jones P D. The influence of ENSO on global temperatures[J]. Climate Monitor, 1988,17(3)：80-89.

[4]施雅風(fēng).全球和中國(guó)變暖特征及未來趨勢(shì)[J].自然災(zāi)害學(xué)報(bào),1996,5(2):5-14.

[5]王紹武.近百年氣候變化與變率的診斷研究[J].氣象學(xué)報(bào),1994,52(3):261-273.

[6]沙萬英,邵雪梅,黃玫. 20世紀(jì)80年代以來中國(guó)的氣候變暖及其對(duì)自然區(qū)域界線的影響[J].中國(guó)科學(xué)(D輯):地球科學(xué),2002,32(4):317-326.

[7]李宗省,何元慶,辛惠娟,等.我國(guó)橫斷山區(qū)1960—2008年氣溫和降水時(shí)空變化特征[J].地理學(xué)報(bào),2010,65(5):563-579.

[8]任國(guó)玉,初子瑩,周雅清,等.中國(guó)氣溫變化研究最新進(jìn)展[J].氣候與環(huán)境研究,2005,10(4):701-716.

[9]李慶祥,李偉,鞠曉慧.1998年以來中國(guó)氣溫持續(xù)極端偏暖的事實(shí)[J].科技導(dǎo)報(bào),2006,24(4):37-40.

[10]秦大河,羅勇,陳振林,等.氣候變化科學(xué)的最新進(jìn)展:IPCC第四次評(píng)估綜合報(bào)告解析[J].氣候變化研究進(jìn)展,2007,3(6):311-314.

[11]葉柏生,李翀,楊大慶,等.我國(guó)過去50 a來降水變化趨勢(shì)及其對(duì)水資源的影響(Ⅰ):年系列[J].冰川凍土,2004,26(5):587-594.

[12]姚惠明,吳永祥,關(guān)鐵生.中國(guó)降水演變趨勢(shì)診斷及其新事實(shí)[J].水科學(xué)進(jìn)展,2013,24(1):1-10.

[13]丁一匯,戴曉蘇.中國(guó)近百年來的溫度變化[J].氣象,1994,20(12):19-26.

[14]郭志梅,繆啟龍,李雄.中國(guó)北方地區(qū)近50年來氣溫變化特征的研究[J].地理科學(xué),2005,25(4):66-72.

[15]王菱,謝賢群,李運(yùn)生,等.中國(guó)北方地區(qū)40年來濕潤(rùn)指數(shù)和氣候干濕帶界線的變化[J].地理研究,2004,23(1):45-54.

[16]張雪芹,孫楊,鄭度,等.中國(guó)干旱區(qū)溫度帶界線對(duì)氣候變暖的響應(yīng)[J].地理學(xué)報(bào),2011,66(9):1166-1178.

[17]代力民,王憲禮,王金錫.三北防護(hù)林生態(tài)效益評(píng)價(jià)要素分析[J].世界林業(yè)研究,2000,13(2):47-51.

[18]鐘中,王漢杰.三北防護(hù)林體系建設(shè)前后華北春季生態(tài)環(huán)境對(duì)比[J].氣象,2001,27(9):13-17.

[19]陳炳浩,郝玉光,陳永富.烏蘭布和沙區(qū)區(qū)域性防護(hù)林體系氣候生態(tài)效益評(píng)價(jià)的研究[J].林業(yè)科學(xué)研究,2003,16(1):63-68.

[20]王曉慧,陳永富,陳爾學(xué),等.基于遙感和GIS的三北防護(hù)林工程生態(tài)效益評(píng)價(jià)研究:以山西省中陽縣為例[J].水土保持通報(bào),2011,31(5):171-175.

[21]趙國(guó)藏.近40年來三北地區(qū)的溫濕變化[J].氣象學(xué)報(bào),1992,50(3):383-384.

[22]王強(qiáng).三北防護(hù)林工程區(qū)植被覆蓋變化特征及其對(duì)氣候變化的響應(yīng)研究[D].蘭州：西北師范大學(xué),2012.

[23]林業(yè)部三北防護(hù)林建設(shè)局.中國(guó)三北防護(hù)林體系建設(shè)[M].北京:中國(guó)林業(yè)出版社,1992：1-3.

[24]朱俊鳳.“三北”防護(hù)林地區(qū)自然資源與綜合農(nóng)業(yè)區(qū)劃[M].北京:中國(guó)林業(yè)出版社，1985:3-6.

[25]申元村,張妙弟.我國(guó)“三北”防護(hù)林地區(qū)土地資源特點(diǎn)及防護(hù)林體系合理布局的探討[J].自然資源,1993,15(5):21-27.

[26]姜鳳岐,于占源,曾德慧,等.氣候變化對(duì)三北防護(hù)林的影響與應(yīng)對(duì)策略[J].生態(tài)學(xué)雜志,2009,28(9):1702-1705.

[27]白愛娟,翟盤茂.中國(guó)近百年氣候變化的自然原因討論[J].氣象科學(xué),2007,27(5):584-590.