鄭 紅，魏 磊，高 峰，潘華盛

(1.黑龍江省氣候中心，黑龍江哈爾濱 150030；2.寶清縣生產(chǎn)力促進(jìn)中心，黑龍江寶清 155600；3.黑龍江省氣象臺(tái)，黑龍江哈爾濱 150030)

七星河是黑龍江省三江平原撓力河左岸的重要支流，分為內(nèi)、外七星河。內(nèi)七星河發(fā)源于雙鴨山市七星砬子山，河水經(jīng)三環(huán)泡狼豁子泄洪道向東至炮臺(tái)亮子入撓力河，全長(zhǎng)241 km，流域面積4 001 km2，其中山丘區(qū)1 600 km2。外七星河發(fā)源于完達(dá)山北麓的雙鴨山，進(jìn)入平原后，河身漫散于沼澤濕地區(qū)，流經(jīng)黑魚(yú)泡滯洪區(qū)折向東，漫行于沼澤區(qū)、名漂筏河，于解放亮子開(kāi)始有明顯河槽，轉(zhuǎn)向東北在菜嘴子以上4 km處匯入撓力河，長(zhǎng)約142 km，流域面積6 520 km2，絕大部分是平原，約占91%，山丘區(qū)面積僅為713 km2。如按平原計(jì)，內(nèi)七星河流域面積為2 401 km2，外七星河流域面積平原部分為5 933 km2，則七星河流域平原部分總面積8 334 km2(圖1)。

近二三十年來(lái)，我國(guó)水利、氣象、地理等部門(mén)科技工作者對(duì)三江平原沼澤濕地下墊面的改變和土地利用變化及其對(duì)濕地的蓄水調(diào)洪功能及徑流變化的影響[1-5]、氣候變化對(duì)水資源和農(nóng)業(yè)氣候資源及空間分布的影響[6-9]、氣候變化和水利工程建設(shè)對(duì)徑流變化的影響[10-15]進(jìn)行了研究，但針對(duì)氣候要素本身降水、蒸發(fā)、氣溫的變化特別氣候變暖情景對(duì)徑流影響的研究卻較少。筆者利用1961—2013年保安和寶清氣象、水文監(jiān)測(cè)資料，基于水量平衡方程，采用相關(guān)、趨勢(shì)和突變檢測(cè)等方法，研究七星河流域中游區(qū)段降水、氣溫、蒸發(fā)及其對(duì)徑流驅(qū)動(dòng)的影響，以期為科學(xué)地合理調(diào)配水資源、綜合利用整個(gè)三江平原水資源和農(nóng)業(yè)發(fā)展提供科學(xué)依據(jù)。

1 資料與方法

1.1資料來(lái)源選取的七星河保安站位于寶清縣西北約30 km處的七星河鄉(xiāng)內(nèi)，距五九七農(nóng)場(chǎng)23 km，是七星河流域中游區(qū)段的唯一代表性水文監(jiān)測(cè)站。所選資料包括1961—2013年逐年七星河流域中游段保安站的降水、徑流資料，撓力河中游寶清水文站的降水、蒸發(fā)、徑流資料，均來(lái)自黑龍江省水文局；寶清氣象站年平均氣溫、五九七農(nóng)場(chǎng)降水資料，來(lái)自黑龍江省氣象信息中心。其中五九七農(nóng)場(chǎng)的降水資料年限為1967—2009年。

1.2研究方法基于水量平衡方程[16]，采用氣候變化趨勢(shì)分析法、相關(guān)系數(shù)分析法、Mann-Kendall突變檢測(cè)法[17]等，分析了七星河流域中游區(qū)段降水、氣溫、蒸發(fā)及其對(duì)徑流驅(qū)動(dòng)的影響。

2 結(jié)果與分析

2.1七星河中游降水

2.1.1保安、寶清與五九七站降水關(guān)系。經(jīng)互相關(guān)計(jì)算，保安與五九七降水的相關(guān)系數(shù)為0.88，保安與寶清降水的相關(guān)系數(shù)為0.80，寶清與五九七降水的相關(guān)系數(shù)為0.69，相關(guān)系數(shù)均達(dá)到了0.01信度檢驗(yàn)，表明保安、寶清、五九七的面雨量關(guān)系是存在的。

圖1 七星河流域周邊地理形勢(shì)Fig.1 Geographic situation around Qixing River Basin

從圖2可以看出，近53年寶清與保安站降水年際變化較一致，保安降水呈稍上升趨勢(shì)，上升幅度為6.5 mm/10 a；寶清降水基本變化不大，上升幅度為0.7 mm/10 a。保安降水多年(1971—2000年)平均為519.7 mm，寶清降水多年(1971—2000年)平均為514.9 mm，保安降水基本持平寶清降水。年降水最大值寶清、五九七站均出現(xiàn)在1981年，分別為826.5、879.5 mm；保安卻出現(xiàn)在1994年，為849.0 mm，可見(jiàn)位于兩站中間的五九七降水較多。1981年，三江平原出現(xiàn)有記錄以來(lái)最嚴(yán)重的洪澇，糧食減產(chǎn)5成左右，個(gè)別低洼地塊造成絕產(chǎn)[18]。年降水最少值保安和寶清出現(xiàn)在1967年，分別為312.7和305.3 mm，五九七出現(xiàn)在1976年，為375.2 mm，三江平原發(fā)生嚴(yán)重干旱，糧食損失慘重。3站最大與最小降水量相差近500.0 mm。

2.1.2降水變化。從圖2可以看出，近53年來(lái)保安與寶清降水變化基本一致，選取代表站保安進(jìn)行降水變化多項(xiàng)式擬合。明顯所見(jiàn)，降水存在多與少時(shí)段的變化；1961—1979年為少水段，平均478.0 mm，較多年平均值少42.0 mm；1980—2002年為多水段，平均為540.0 mm，較多年平均值多20.0 mm；2003—2011年為少水段，平均為474.0 mm，較多年平均值少46.0 mm。2012年后保安降水將由少水轉(zhuǎn)為多水時(shí)段，這一結(jié)論是與全省降水由少水轉(zhuǎn)多水時(shí)段觀點(diǎn)相一致的[19]。

圖2 1961—2013年寶清與保安降水量變化Fig.2 Variation of precipitation in Baoqing and Bao’an during 1961-2013

2.2七星河與撓力河徑流變化關(guān)系選取撓力河中游站寶清和相距20多公里七星河中游站保安，兩站地貌特征生態(tài)環(huán)境一致，而且同在一個(gè)面雨量控制區(qū)域內(nèi)，因此兩站年徑流變化趨勢(shì)也是非常一致的(圖3)。經(jīng)計(jì)算，兩站年徑流量相關(guān)系數(shù)為0.89，遠(yuǎn)遠(yuǎn)大于0.01信度檢驗(yàn)。寶清徑流增加則保安徑流就增加，寶清年徑流以5 994.8 萬(wàn)m3/10 a在減少，保安年徑流以2 180.6 萬(wàn)m3/10 a同樣在減少。寶清多年平均徑流量43 561萬(wàn)m3，保安多年平均徑流量15 762萬(wàn)m3，寶清年徑流為保安的2.8倍。從寶清年徑流變化多項(xiàng)式擬合曲線來(lái)看，1961—1981年為徑流增多期，平均流量為54 959萬(wàn)m3，較歷年平均流量(43 561萬(wàn)m3)多11 398萬(wàn)m3；同期保安徑流較歷年平均多4 335萬(wàn)m3；但這一時(shí)期徑流量波動(dòng)也較大，近53年4次大的徑流1964、1973、1974和1981年都發(fā)生在此期，尤其1981年發(fā)生特大流量，寶清為132 600萬(wàn)m3，是歷年平均流量的3倍，保安為55 500萬(wàn)m3，是歷年的3.5倍，造成三江平原嚴(yán)重洪澇；而1976—1979年也是近53年來(lái)徑流量最少的4年，寶清平均不到10 000萬(wàn)m3，而保安為4 475萬(wàn)m3。1982—1997年為徑流平穩(wěn)期，平均流量為40 467萬(wàn)m3，較歷年平均稍多3 094萬(wàn)m3，基本變化不大；同期保安徑流量為15 756萬(wàn)m3，與歷年徑流持平。1998—2009年為徑流減少期，寶清平均流量為28 009萬(wàn)m3，較歷年平均減少17 552萬(wàn)m3；同期保安徑流量為11 162萬(wàn)m3，較歷年減少71%。2010年后特別是2013年后徑流量將開(kāi)始增多。

圖3 1961—2013年寶清與保安年徑流變化Fig.3 Annual runoff changes of Baoqing and Bao’an during 1961-2013

圖4 1961—2013年保安年降水與徑流量關(guān)系Fig.4 Relationship between precipitation and runoff in Bao’an during 1961-2013

2.3降水變化對(duì)七星河中游徑流驅(qū)動(dòng)的影響從圖4可以看出，近53年來(lái)七星河中游的保安降水呈微升趨勢(shì)，徑流呈微降趨勢(shì)，保安降水與徑流均具有10年左右的重現(xiàn)期。經(jīng)相關(guān)計(jì)算，降水與徑流之間的相關(guān)系數(shù)達(dá)0.72，遠(yuǎn)大于0.01信度檢驗(yàn)，表明保安降水多則七星河中游徑流量就多，降水少則徑流少；但從撓力河與七星河之間的面雨量代表站寶清來(lái)說(shuō)，與七星河徑流量相關(guān)系數(shù)為0.74，高于保安與徑流的相關(guān)系數(shù)，尤其寶清氣象站的年降水與保安徑流的相關(guān)系數(shù)為0.79，相關(guān)關(guān)系更大一些。經(jīng)統(tǒng)計(jì)，保安降水≥600 mm的年份共有13年，則對(duì)應(yīng)保安徑流除1988、2000、2012年偏少外，其余10年徑流偏多；寶清1988和2012年降水都是偏少的(表1)。如保安降水距平百分率取≥120%為降水特多年，則共有9年，僅2000年徑流偏少，其他8年則對(duì)應(yīng)徑流偏多2.5成。表明降水量是驅(qū)動(dòng)徑流變化的主要原因之一。另外，從七星河徑流的減少趨勢(shì)來(lái)看，主要是在1982年后，土地利用發(fā)生較大變化，它中游區(qū)域有大量水稻田種植面積的增加，如五九七農(nóng)場(chǎng)水稻1984年僅有13 hm2，到2009年已增至16 000 hm2，擴(kuò)大了1 231倍。友誼農(nóng)場(chǎng)由1984年的377 hm2已增加至2006年的21 347 hm2，擴(kuò)大了57倍。由于水稻面積大幅增加，抽取大量地下水和地表水灌溉，使地下水減少，水位下降，這樣增加了地表徑流補(bǔ)充土壤包氣帶的水分，影響了地表徑流的流出量，導(dǎo)致了七星河徑流量的減少。這也表明土地利用的變化和水利工程發(fā)展也是驅(qū)動(dòng)地表徑流減少的另一主要原因之一。

表1保安、寶清氣象站年降水量≥600mm和對(duì)應(yīng)保安年徑流量

Table1Annualprecipitation(≥600mm)andcorrespondingannualrunoffofBao’anandBaoqingweatherstation

年份Year降水量Precipitation∥mm保安Bao’an寶清Baoqing徑流Runoff∥萬(wàn)m319626286062517119716255162486319726297312505119746006963728119817998285550019886054631232019917486223429019948497104549020006335699207200268164219800200961959818030201261650012790201364566424830

2.4蒸發(fā)變化對(duì)七星河徑流的影響選取與保安距離23 km的唯一有水面蒸發(fā)監(jiān)測(cè)的寶清水文站，它的蒸發(fā)過(guò)程完全可以代表七星河中游段的蒸發(fā)。從圖5可以看出，近53年來(lái)寶清蒸發(fā)呈稍下降趨勢(shì)，以14.6 mm/10 a的速率下降，這與黑龍江省風(fēng)速減少有關(guān)[20]；蒸發(fā)與徑流相關(guān)系數(shù)為-0.25，很接近信度0.05檢驗(yàn)，表明蒸發(fā)量大則徑流就減少的反相關(guān)關(guān)系。但從蒸發(fā)對(duì)降水時(shí)段變化來(lái)看表現(xiàn)是不同的，多年(1971—2000年)平均蒸發(fā)量為753 mm，1961—1974年呈基本持平趨勢(shì)，1975—1987年呈增加趨勢(shì)，1988—2007年呈減少趨勢(shì)，2008年后呈增加趨勢(shì)。1961—1987年蒸發(fā)與對(duì)應(yīng)徑流變化的相關(guān)系數(shù)達(dá)-0.55，通過(guò)信度0.01檢驗(yàn)；而1988—2013年相關(guān)系數(shù)為-0.14，相關(guān)很低，沒(méi)有關(guān)系可尋；也可認(rèn)為寶清的水面蒸發(fā)由1988年開(kāi)始發(fā)生了轉(zhuǎn)折性變化。1961—1987年蒸發(fā)與徑流呈反位相變化，其中1976—1979年蒸發(fā)量較歷年平均多2成以上，對(duì)應(yīng)徑流卻出現(xiàn)了53年來(lái)極低的4年。1988—2013年反相關(guān)關(guān)系的年份僅有10年，其中1995—2008年中僅2002年為反位相變化，其他13年為正相關(guān)關(guān)系；可見(jiàn)1988年后蒸發(fā)與徑流的反相關(guān)關(guān)系已不復(fù)存在，這也充分表明此時(shí)期也正是大力擴(kuò)大水稻面積，引水灌溉，開(kāi)渠打井，地下水大大減少，造成徑流減少。因此，土地利用和水利工程建設(shè)改變了蒸發(fā)對(duì)徑流的驅(qū)動(dòng)反相關(guān)影響。

綜上所述，表明蒸發(fā)對(duì)徑流的驅(qū)動(dòng)影響基本呈反相變化關(guān)系，即蒸發(fā)量大則徑流量減少；反之亦然。尤其1961—1987年蒸發(fā)對(duì)徑流的驅(qū)動(dòng)影響明顯；而1988—2013年蒸發(fā)對(duì)徑流的驅(qū)動(dòng)影響的關(guān)系已不明顯，因此水利工程和降水變化成為此時(shí)段對(duì)徑流驅(qū)動(dòng)的主要影響因子。

圖5 1961—2013年寶清水面蒸發(fā)與保安徑流Fig.5 Water surface evaporation of Baoqing and runoff of Bao’an during 1961-2013

2.5氣候變化對(duì)七星河徑流驅(qū)動(dòng)的影響

2.5.1氣溫變化對(duì)徑流驅(qū)動(dòng)的影響。選取七星河與撓力河之間唯一代表站寶清氣象站的年氣溫和七星河保安徑流，從圖6可以看出，近53年寶清年平均氣溫以0.37 ℃/10 a的速率上升，而徑流以2 180萬(wàn)m3/10 a的速率下降，寶清多年(1971—2000年)年平均氣溫為4.0 ℃。經(jīng)計(jì)算，氣溫與徑流兩序列相關(guān)系數(shù)達(dá)-0.31，通過(guò)信度0.01檢驗(yàn)，反相關(guān)關(guān)系明顯，表明氣溫升高則徑流減少，反之氣溫降低則徑流就大。近53年來(lái)寶清年均氣溫高低時(shí)段變化明顯，為準(zhǔn)確地劃分氣溫時(shí)段的變化，常采用Mann-Kendall突變趨勢(shì)法[17]計(jì)算，結(jié)果表明，1988年為突變開(kāi)始時(shí)間，突變前(1961—1987年)為低溫時(shí)段，平均氣溫為3.4 ℃，低于歷年平均0.6 ℃；其中假設(shè)低于年平均氣溫1.0 ℃為特低年，則有1964、1965、1966、1969、1972、1974、1977、1984年共8年，平均3年就有1次特低年發(fā)生；1969年發(fā)生了53年以來(lái)的極值低溫年，較歷年平均低2.1 ℃。突變后(1988—2013年)為氣溫較高時(shí)段，平均氣溫為4.7 ℃，高于歷年平均0.6 ℃，2個(gè)時(shí)段均差為1.3 ℃；在此氣溫偏高時(shí)段里，1次低溫年也未發(fā)生，僅有2年稍低年(2000、2009年)低于歷年平均0.3 ℃，表明氣候變暖后低溫年數(shù)極大減少，這與黑龍江省及全國(guó)年氣溫變化具有同步一致的變化特征[21]。

根據(jù)上述結(jié)果，分析氣溫突變前后時(shí)段所對(duì)應(yīng)的保安徑流時(shí)段變化。突變前(1961—1987年)氣溫與徑流量相關(guān)系數(shù)為-0.23(未通過(guò)檢驗(yàn))，27年平均徑流量為18 751萬(wàn)m3，高于多年平均2成。突變后(1988—2013年)氣溫與徑流量相關(guān)系數(shù)為-0.15(未通過(guò)檢驗(yàn))，26年平均徑流量13 432萬(wàn)m3，少于多年平均1.5成。表明突變前氣溫偏低減少了地表土壤和水面蒸發(fā)量，造成徑流偏多，反之突變后氣溫偏高、蒸發(fā)量加大則徑流偏少。從個(gè)例分析來(lái)看，突變前低溫年取年平均氣溫距平≤0.5 ℃的年份，有1962、1964、1965、1966、1969、1972、1974、1976、1977、1980、1981、1984、1987年共13年，其中1976、1977、1980、1984年對(duì)應(yīng)徑流偏少，1976和1977年徑流偏少完全是由于降水特少造成；其他年份對(duì)應(yīng)徑流偏多，尤其1981年是近53年來(lái)徑流出現(xiàn)最高值，不但降水最多，且年均氣溫低于歷年0.5 ℃。突變后高溫年取年平均氣溫距平≥0.5 ℃的年份，有1975、1982、1989、1990、1992、1993、1994、1995、1997、1998、2001、2002、2003、2004、2007、2008、2011、2013年共18年，其中1994、2002、2013年對(duì)應(yīng)徑流偏多或特多；其他15年則對(duì)應(yīng)徑流偏少。尤其2008年氣溫較歷年偏高1.8 ℃，為近53年內(nèi)最高值，則對(duì)應(yīng)徑流偏少53%。表明七星河流域氣候暖化后氣溫偏高只要達(dá)到≥0.5 ℃，將加速?gòu)搅鞯臏p少。由此可見(jiàn)，氣溫變化雖不能對(duì)徑流變化直接產(chǎn)生驅(qū)動(dòng)影響，但能起到加速或減少蒸發(fā)的驅(qū)動(dòng)作用。

圖6 1961—2013年寶清年均氣溫與保安徑流Fig.6 Average annual temperature of Baoqing and runoff of Bao’an during 1961-2013

圖7 1988—2013年保安年降水與徑流Fig.7 Annual precipitation and runoff of Bao’an during 1988-2013

2.5.2氣候變暖大背景下降水對(duì)徑流驅(qū)動(dòng)影響。氣候變暖導(dǎo)致了極端天氣頻發(fā)，降水減少，干旱增加[20]，從保安降水對(duì)徑流變化的影響也可揭示這一點(diǎn)。前面已給出氣候突變時(shí)間在1988年，經(jīng)計(jì)算，突變前(1961—1987年)保安降水與徑流的相關(guān)系數(shù)為0.78，而突變后降水與徑流的相關(guān)系數(shù)為0.82，表明突變后降水在減少而徑流量也在減少的年份大大增加(圖7)。從降水變化趨勢(shì)看，突變后降水以19.0 mm/10 a的速率在減少，多于突變前降水減少速率(7.5 mm/10 a)，同時(shí)徑流也以1 825.8萬(wàn)m3/10 a的速率在減少。統(tǒng)計(jì)表明，在突變后的26年中，出現(xiàn)降水偏多12年則對(duì)應(yīng)徑流偏多有6年，出現(xiàn)降水偏少12年則對(duì)應(yīng)徑流偏少也是12年(其中有2年降水正常不再統(tǒng)計(jì)之內(nèi))，而且其中有8年降水特少(≤85%)，則對(duì)應(yīng)有8年徑流也特少。綜上表明，氣候變暖，徑流減少，降水成為驅(qū)動(dòng)徑流變化影響的主要因子。

3 結(jié)論與討論

通過(guò)對(duì)近53年來(lái)七星河流域中游段降水、氣溫、蒸發(fā)及對(duì)徑流驅(qū)動(dòng)變化關(guān)系的分析，特別是突變前后(即變暖前后)2個(gè)時(shí)段變化特征的對(duì)比分析，結(jié)果表明，大氣降水在水循環(huán)過(guò)程中仍然是驅(qū)動(dòng)徑流變化的主要因子；蒸發(fā)對(duì)徑流具有明顯的反位相變化特征，尤其是在氣候突變前1961—1987年，蒸發(fā)起到了加速和減少?gòu)搅鞯尿?qū)動(dòng)影響。而氣溫變化雖不能對(duì)徑流變化直接產(chǎn)生驅(qū)動(dòng)影響，但卻能起到增強(qiáng)或減弱蒸發(fā)的驅(qū)動(dòng)作用和對(duì)降水起到不同冷暖氣候條件下的降水對(duì)徑流變化的驅(qū)動(dòng)影響。另外還分析了氣候變化特別是土地利用變化及水利工程建設(shè)對(duì)徑流變化起著的驅(qū)動(dòng)作用。該研究對(duì)于在氣候變化背景下規(guī)劃七星河、撓力河流域水利工程建設(shè)和水資源綜合利用及管理，實(shí)現(xiàn)農(nóng)業(yè)可持續(xù)發(fā)展具有現(xiàn)實(shí)和長(zhǎng)遠(yuǎn)意義。

[1] 陳剛起，張文芬.三江平原沼澤對(duì)河川徑流影響的初步探討[J].地理科學(xué)，1982，2(3)：254-263.

[2] 劉興土.三江平原沼澤濕地的蓄水與調(diào)洪功能 [J].濕地科學(xué)，2007，5(1)：64-68.

[3] 王薇，李傳奇.蓄滯洪區(qū)的功能、價(jià)值與多目標(biāo)利用 [J].水利發(fā)展研究，2004，4(9)：26-28.

[4] 陳剛起.三江平原沼澤徑流的實(shí)驗(yàn)研究 [C]//黃錫疇.中國(guó)沼澤研究.北京：科學(xué)出版社，1989：120-125.

[5] 于宏敏，劉赫男，潘華盛，等.黑龍江七星河濕地面積銳減原因的研究[J].自然災(zāi)害學(xué)報(bào)，2013，23(6)：234-239.

[6] 潘華盛，張桂華，李向東.氣候變暖對(duì)黑龍江省水資源的影響及對(duì)策[J].黑龍江水專學(xué)報(bào).2004，31(1)：1-4.

[7] 李彩俠，李俏，王雅珍，等.近52年黑龍江省農(nóng)業(yè)氣候資源變化特征分析[J].水土保持研究，2014，21(6)：187-192.

[8] 王英，曹明奎，陶波，等.全球氣候變化背景下中國(guó)降水量空間格局的變化特征[J].地理研究，2006，25(6)：1031-1040.

[9] 王國(guó)慶，張建云，劉九夫，等.氣候變化和人類活動(dòng)對(duì)河川徑流影響的定量分析[J].中國(guó)水利，2008(2)：55-58.

[10] 陳晶，楊帆，田寶興，等．黑龍江省黑土區(qū)近50年來(lái)四季氣溫時(shí)空變化特征分析[J].水土保持研究，2013，20(3)：150-154.

[11] 林凱榮，何艷虎，陳曉宏.氣候變化及人類活動(dòng)對(duì)東江流域徑流影響的貢獻(xiàn)分解研究[J].水利學(xué)報(bào)，2012，43(11)：1312-1321.

[12] 劉春蓁，夏軍.氣候變暖條件下水文循環(huán)變化檢測(cè)與歸因研究的幾點(diǎn)認(rèn)識(shí)[J].氣候變化研究進(jìn)展，2010，6(5)：313-318.

[13] 李峰平，章光新，董李勤.氣候變化對(duì)水循環(huán)與水資源的影響研究綜述[J].地理科學(xué)，2013，33(4)：457-464.

[14] 陳亞寧，李稚，范煜婷，等.西北干旱區(qū)氣候變化對(duì)水文水資源影響研究進(jìn)展[J].地理學(xué)報(bào)，2014，69(9)：1295-1304.

[15] 王國(guó)慶，張建云，賀瑞敏.環(huán)境變化對(duì)黃河中游汾河徑流情勢(shì)的影響研究[J].水科學(xué)進(jìn)展，2006，17(6)： 853-858.

[16] 徐冬梅，劉曉民.水文水利計(jì)算[M].鄭州：黃河水利出版社，2013：11-12.

[17] 魏鳳英.現(xiàn)代氣候統(tǒng)計(jì)診斷與預(yù)測(cè)技術(shù)[M].2版.北京：氣象出版社，2007：43-47，20-26，69-71.

[18] 潘華盛，劉育生.三江平原天氣氣候與農(nóng)業(yè)[M].北京：氣象出版社，1991：88-106.

[19] 潘華盛，王勇，鄒濱，等.黑龍江省2012年后由暖干轉(zhuǎn)暖濕型氣候變化及預(yù)測(cè)分析[J].黑龍江大學(xué)工程學(xué)報(bào)，2013，4(3)：45-50.

[20] 周秀杰，張桂華，鄭紅，等.黑龍江省氣候變暖對(duì)極端天氣氣事件的影響[J].氣象，2004，30(11)：47-50.

[21] 周秀杰，王風(fēng)玲，潘華盛，等.近60年來(lái)黑龍江省與東北及全國(guó)氣溫變化特點(diǎn)分析[J].自然災(zāi)害學(xué)報(bào)，2013，22(2)：124-129.