氣候變化對(duì)黑龍江省自然植被氣候生產(chǎn)潛力的影響

田寶星 宮麗娟 楊帆 陳晶 翟墨 張洋

摘要：為了解1970—2014年黑龍江省自然植被氣候生產(chǎn)潛力的分布規(guī)律和變化趨勢(shì)及其對(duì)氣候變化的響應(yīng)，利用黑龍江省80個(gè)氣象站觀測(cè)數(shù)據(jù)，結(jié)合Miami模型，通過小波分析和EOF分析等方法，對(duì)黑龍江省自然植被氣候生產(chǎn)潛力的變化特征進(jìn)行分析。結(jié)果表明：1970年以來，氣溫、降水、蒸散和標(biāo)準(zhǔn)氣候生產(chǎn)潛力呈增加趨勢(shì)，2000—2014年糧食安全存在隱患。標(biāo)準(zhǔn)氣候生產(chǎn)潛力（W）與蒸散氣候生產(chǎn)潛力年代際變化相近。4種氣候生產(chǎn)潛力主周期均為26～31年。W空間分布呈正、負(fù)相間。W與氣溫、降水量和蒸散量呈正相關(guān)，氣溫是主要限制因子。“暖濕型”氣候?qū)有利，生產(chǎn)潛力增加幅度7.41%～ 14.82%。研究結(jié)果對(duì)優(yōu)化農(nóng)業(yè)布局，提高農(nóng)業(yè)氣候資源利用率以及促進(jìn)糧食作物穩(wěn)產(chǎn)高產(chǎn)等有參考價(jià)值。

關(guān)鍵詞：植被；氣候變化；Miami模型；氣候生產(chǎn)潛力；黑龍江

中圖分類號(hào)：P49文獻(xiàn)標(biāo)志碼：A論文編號(hào)：cjas20190900204

Effect of Climate Change on Climate-induced Potential Productivity of Natural Vegetation in Heilongjiang

Tian Baoxing1，2， Gong Lijuan1，2， Yang Fan3， Chen Jing4， Zhai Mo1，2， Zhang Yang3

（1Heilongjiang Province Institute of Meteorological Science， Harbin 150030， Heilongjiang， China； 2Innovation and Opening Laboratory of Regional Eco-meteorology in Northeast of China Meteorological Administration， Harbin 150030， Heilongjiang， China；3Heilongjiang Meteorological Bureau Office Service Centre， Harbin 150030， Heilongjiang， China；4Bayan Meteorological Bureau， Harbin 151800， Heilongjiang， China）

Abstract： To understand the distribution rule and change trend of climate-induced potential productivity of natural vegetation and its response to climate change in Heilongjiang from 1970 to 2014， the authors used the Miami model， wavelet analysis and EOF to analyze the variation characteristics of the climate- induced potential productivity of natural vegetation based on observation data of 80 meteorological stations. The results showed that： the four types of climate- induced potential productivity （air temperature， precipitation， evapotranspiration and standard climate-induced potential productivity） had increased since 1970， however， there were hidden dangers in food security from 2000 to 2014. The standard climate- induced potential productivity （W） and evapotranspiration climate- induced potential productivity had similar interdecadal variation. The major periodic of the four types of climate-induced potential productivity was 26 to 31 years.The W presented alternate positive and negative variations in space， and it was positively correlated with the air temperature， precipitation， and evapotranspiration， the air temperature was the main restriction factor. Warm-wet climate is beneficial to W， and the potential productivity will increase within a range of 7.41% to 14.82% . The results could provide a basic method for optimizing agricultural distribution， improving the utilization of agricultural climate resources and promoting stable and high yield production of grain crops.

Keywords： Vegetation； Climate Change； Miami Model； Climate-induced Potential Productivity； Heilongjiang

0引言

氣候變化是人類面臨的最為嚴(yán)重的環(huán)境問題之一，氣候變暖日益明顯，使得地表溫度及氣溫明顯高于以往，同時(shí)極端天氣也給農(nóng)業(yè)生產(chǎn)增加不確定性[1-2]。氣候變化已經(jīng)對(duì)人類及自然植被產(chǎn)生了廣泛和重要的影響[3]，特別是以農(nóng)業(yè)為主體的糧食安全[4]。氣候生產(chǎn)潛力指某地區(qū)光溫水等因子相互協(xié)調(diào)，作物合理利用當(dāng)?shù)貧夂蛸Y源和土壤養(yǎng)分，單位面積產(chǎn)量將達(dá)到最大值[5]。探討自然植被氣候生產(chǎn)潛力，體現(xiàn)了氣象因子的不同及相互協(xié)調(diào)程度，不僅可以預(yù)測(cè)植被生產(chǎn)潛力，可以依據(jù)氣候變化趨勢(shì)來預(yù)測(cè)未來的發(fā)展情況，是評(píng)估農(nóng)業(yè)氣候資源的依據(jù)之一[6]。

黑龍江省土壤肥沃，自然條件適宜，是中國(guó)重要的糧食主產(chǎn)區(qū)和最大的商品糧基地[7]。大豆是中國(guó)重要的糧油作物之一，黑龍江省是中國(guó)大豆種植面積和產(chǎn)量較高的地區(qū)之一，在中國(guó)大豆生產(chǎn)中占有非常重要的作用[8]。玉米是中國(guó)第一大糧食作物，黑龍江省玉米播種面積占全國(guó)10%，產(chǎn)量占全國(guó)9%左右，與國(guó)家糧食安全息息相關(guān)[9]。黑龍江水稻品質(zhì)佳，2016年水稻產(chǎn)量占全國(guó)稻谷總產(chǎn)量10.9%，因此黑龍江省水稻安全生產(chǎn)對(duì)國(guó)家糧食安全意義重大[10]。研究黑龍江省的自然植被氣候生產(chǎn)潛力，基于統(tǒng)計(jì)學(xué)原理分析不同時(shí)間的植被生產(chǎn)潛力的變化，可為實(shí)現(xiàn)優(yōu)質(zhì)高產(chǎn)高效的農(nóng)業(yè)生產(chǎn)模式提供可靠的依據(jù)和參考。

氣候生產(chǎn)潛力是指作物不受養(yǎng)分限制時(shí)，由光照、氣溫和降水共同決定的生產(chǎn)力，是一個(gè)地區(qū)作物產(chǎn)量的上限[11]。國(guó)內(nèi)外相關(guān)研究已有很多，如Miami模型、Thornthwaite Memorial模型、AEZ模型和逐級(jí)訂正模型等[12-15]。研究表明[16-18]，Miami模型體現(xiàn)了氣溫和降水因素，Thornthwaite Memorial模型則通過蒸散量來模擬生物量。

眾多學(xué)者已從不同層面、角度對(duì)黑龍江省氣候生產(chǎn)潛力進(jìn)行研究。如張興林等[19]為充分發(fā)揮黑龍江省的氣候資源優(yōu)勢(shì)，探求玉米種植新模式，探討富裕縣移栽玉米氣候生產(chǎn)潛力。王秀芬等[14]基于AEZ模型估算了黑龍江省光溫和氣候生產(chǎn)潛力。盧玢宇等[20]基于Thornthwaite Memorial模型估算了黑龍江省氣候生產(chǎn)潛力并分析氣候生產(chǎn)潛力與平均溫度和降水量之間的關(guān)系。杜國(guó)明等[21]基于GAEZ模型分析氣候變化導(dǎo)致的東北玉米生產(chǎn)潛力時(shí)空格局演變特征。上述研究所用站點(diǎn)資料較少，不足以代表黑龍江省氣候生產(chǎn)潛力狀況。本研究根據(jù)黑龍江省80個(gè)氣象站1970—2014年逐日氣象資料，采用Thornthwaite Memorial模型和Miami模型相結(jié)合估算黑龍江省自然植被標(biāo)準(zhǔn)氣候生產(chǎn)潛力及其變化，分析氣候變化背景下黑龍江省自然植被的溫度氣候生產(chǎn)潛力、降水氣候生產(chǎn)潛力、蒸散氣候生產(chǎn)潛力和標(biāo)準(zhǔn)氣候生產(chǎn)潛力的分布特征，以及熱量和降水資源對(duì)該區(qū)域氣候生產(chǎn)潛力的影響，以期為該區(qū)提高農(nóng)業(yè)生產(chǎn)力及可持續(xù)發(fā)展，合理利用氣候資源及糧食生產(chǎn)安全提供參考。

1材料與方法

1.1資料來源和區(qū)域概況

黑龍江省介于121°11’—135°05’E、43°26’—53°33’N之間，屬大陸季風(fēng)性氣候，全省氣候的主要特點(diǎn)是春季干旱，夏季多雨，秋季早霜，冬季寒冷。當(dāng)?shù)赜隉嵬厩彝恋刭Y源對(duì)農(nóng)業(yè)發(fā)展有得天獨(dú)厚的優(yōu)勢(shì)，盛產(chǎn)大豆、玉米、水稻等糧食作物[22-23]。

選取黑龍江省80個(gè)氣象站1970—2014年逐日氣溫和降水資料，材料來自黑龍江省氣象數(shù)據(jù)中心。

1.2研究方法

1.2.1氣候生產(chǎn)潛力模型在植被生產(chǎn)潛力估算模型中，Miami模型從植被的生理生態(tài)角度出發(fā)，通過計(jì)算該區(qū)年降水量和平均氣溫等要素來確定植被生產(chǎn)潛力，中國(guó)學(xué)者通過氣象資料將模型進(jìn)行驗(yàn)證，認(rèn)為該模型可以應(yīng)用在自然植被生產(chǎn)潛力研究[24]。Thornthwaite Memorial模型是Miami模型的演變和延伸，展現(xiàn)植被生產(chǎn)潛力受氣候因素變化的影響機(jī)理[25-26]。

1.3統(tǒng)計(jì)分析

基于研究區(qū)域逐日氣溫降水資料，采用Microsoft Office Excel工作表進(jìn)行數(shù)據(jù)整理及分析，繪制折線圖和柱狀圖；采用Matlab軟件，根據(jù)小波分析中的Morlet變換，編寫數(shù)據(jù)處理程序，繪制小波方差圖；采用Surfer軟件，將EOF函數(shù)特征向量空間變化特征展現(xiàn)出來。

2結(jié)果與分析

2.1黑龍江省自然植被氣候生產(chǎn)潛力變化特征

2.1.1年代際變化特征由圖1和表1可見，黑龍江省自然植被氣候生產(chǎn)潛力在1970—2014年的序列變化，氣溫、降水和蒸散生產(chǎn)潛力變化均呈現(xiàn)增加趨勢(shì)，氣候傾向率分別為21.82、13.84、13.33 kg/（hm2·a）。在年代際尺度中，1970s（1970—1979年）、1980s（1980—1989年）和1990s（1990—1999年）3種生產(chǎn)潛力變化均呈增加趨勢(shì)，說明氣象條件變化有利于氣候生產(chǎn)潛力的發(fā)揮；2000s（2000—2009年）、2010s（2010—2014年）3種生產(chǎn)潛力變化不穩(wěn)定，受氣溫和降水等因素限制，糧食安全存在安全隱患。標(biāo)準(zhǔn)氣候生產(chǎn)潛力（W）呈增加趨勢(shì)，平均值為7015 kg/（hm2·a），最大值為7778 kg/（hm2·a），最小值為6145 kg/（hm2·a），氣候傾向率為13.02 kg/（hm2·a），從整個(gè)研究時(shí)域看，標(biāo)準(zhǔn)氣候生產(chǎn)潛力（W）與蒸散氣候生產(chǎn)潛力（WV）由大到小排列，均為1990s>2010s>2000s>1980s>1970s。水熱配比情況，1970s、1980s、1990s和2010s這4個(gè)時(shí)期WT與WR的比值均大于0.85，可見水熱配比狀態(tài)較好，氣候生產(chǎn)潛力的發(fā)揮，熱量起主要限制作用。而2000s時(shí)期WT與WR的比值大于1，可見這個(gè)時(shí)期熱量資源充足，降水對(duì)氣候生產(chǎn)潛力發(fā)揮也起到限制作用。從整體看，研究時(shí)域黑龍江省水熱配比情況較好（WT與WR的比值≥0.85），熱量資源對(duì)氣候生產(chǎn)潛力的發(fā)揮影響更大。

2.1.2周期變化小波分析見圖2。研究區(qū)域1970—2014年4種氣候生產(chǎn)潛力的周期變化，主要體現(xiàn)在26～31、18～24、11～16年等3種時(shí)間尺度，體現(xiàn)研究區(qū)域4種氣候生產(chǎn)潛力周期變化的一致性。由圖2可見，W、WT、WR、WV4種氣候生產(chǎn)潛力小波方差的主要峰值依次對(duì)應(yīng)28、28、27、28年左右尺度，為第一主周期且周期振蕩最強(qiáng)；第二主周期為15、22、5、15年；第三主周期為5、13、15、5年。這3個(gè)周期尺度掌控研究時(shí)域4種氣候生產(chǎn)潛力的變化特征。

2.1.3空間變化特征在EOF分析中，研究區(qū)域自然植被標(biāo)準(zhǔn)氣候生產(chǎn)潛力（W）的前3個(gè)特征向量的方差貢獻(xiàn)率分別是45.02%、12.57%和7.26%，累計(jì)方差貢獻(xiàn)率是64.85%，因此可以體現(xiàn)研究區(qū)域自然植被標(biāo)準(zhǔn)氣候生產(chǎn)潛力（W）空間變化特征。

第1特征向量方差貢獻(xiàn)率是45.02%，體現(xiàn)研究區(qū)域自然植被標(biāo)準(zhǔn)氣候生產(chǎn)潛力變化的顯著空間分布狀態(tài)。由圖3a可見，黑龍江省自然植被標(biāo)準(zhǔn)氣候生產(chǎn)潛力（W）呈現(xiàn)正變化，區(qū)內(nèi)高值中心在泰來附近且等值線密集，次高值中心在雙鴨山附近且等值線較密集，低值中心在呼中附近且等值線相對(duì)稀疏，其振幅以西部和東部為中心向周圍逐漸遞減，標(biāo)準(zhǔn)氣候生產(chǎn)潛力受地理因素影響較大，區(qū)域特征明顯，水熱條件匹配較好的區(qū)域標(biāo)準(zhǔn)氣候生產(chǎn)潛力相對(duì)偏高，同時(shí)反映了植被生長(zhǎng)狀況地域性顯著，偏北（偏南）地區(qū)氣溫較低（高），均抑制植被的生長(zhǎng)。由圖3b可見，第1特征向量所對(duì)應(yīng)的時(shí)間系數(shù)絕對(duì)值最大，這一時(shí)刻空間分布狀態(tài)最典型，但時(shí)間系數(shù)正、負(fù)相間，標(biāo)準(zhǔn)氣候生產(chǎn)潛力年際變化相反。

第2特征向量方差貢獻(xiàn)率是12.57%，是研究區(qū)域自然植被標(biāo)準(zhǔn)氣候生產(chǎn)潛力一種顯著的空間分布。由圖4a可見，黑龍江省自然植被標(biāo)準(zhǔn)氣候生產(chǎn)潛力（W）呈現(xiàn)正、負(fù)變化趨勢(shì)，高值中心在泰來附近，中心絕對(duì)值在0.25以上，低值中心在克東附近。西部部分地區(qū)第2特征向量為正值，其他地區(qū)為負(fù)值，說明其標(biāo)準(zhǔn)氣候生產(chǎn)潛力變化型式相反。這種分布類型，突顯了西部部分地區(qū)地理分布狀況差異，水熱匹配條件影響標(biāo)準(zhǔn)氣候生產(chǎn)潛力的高低。圖4b可見第2特征向量對(duì)應(yīng)的時(shí)間系數(shù)絕對(duì)值次之，時(shí)間系數(shù)正、負(fù)相間，說明在時(shí)間系數(shù)正值（負(fù)值）的年份里，西部部分地區(qū)標(biāo)準(zhǔn)氣候生產(chǎn)潛力偏高（低），其他地區(qū)標(biāo)準(zhǔn)氣候生產(chǎn)潛力偏低（高），氣候生產(chǎn)潛力分布型式相反特征明顯。

第3特征向量方差貢獻(xiàn)率是7.26%，對(duì)研究區(qū)域自然植被標(biāo)準(zhǔn)氣候生產(chǎn)潛力空間分布有一定的意義。由圖5a可見，西南部和東南部部分地區(qū)第3特征向量為正值，其他地區(qū)為負(fù)值，區(qū)內(nèi)高值中心在寶清附近，低值中心在哈爾濱附近。圖5b可以看出時(shí)間系數(shù)絕對(duì)值較小，自然植被標(biāo)準(zhǔn)氣候生產(chǎn)潛力分布狀態(tài)較典型，第3特征向量和時(shí)間系數(shù)特征與第2特征向量相似。 2.2黑龍江省自然植被氣候生產(chǎn)潛力對(duì)氣候變化的響應(yīng)

2.2.1氣象因子和自然植被氣候生產(chǎn)潛力相互之間關(guān)系由圖6可見，黑龍江省自然植被W與氣溫、降水和蒸散量的變化一致，相關(guān)系數(shù)分別為0.371、0.832和0.995，均通過0.05的信度檢驗(yàn)。W與年氣溫、年降水量和年蒸散量定量正相關(guān)，即年平均氣溫每上升1℃、年降水量和年蒸散量每增加10 mm，W依次變化為199.91、44.46、233.46 kg/（hm2·a）。由此可見，黑龍江省地區(qū)蒸散量對(duì)W的影響較大，根據(jù)前文公式可知，氣溫對(duì)W的影響占主導(dǎo)地位，降水量對(duì)W的影響相對(duì)較小。

2.2.2標(biāo)準(zhǔn)氣候生產(chǎn)潛力對(duì)氣溫和降水量的敏感性分析為分析氣象因子對(duì)標(biāo)準(zhǔn)氣候生產(chǎn)潛力的影響，探討W與年平均氣溫、年降水量的定量關(guān)系，通過檢驗(yàn)并建立擬合回歸方程，如式（8）所示。

方差分析得出，F(xiàn)統(tǒng)計(jì)量值為319.470，顯著性概率P<0.001，回歸方程非常顯著，且正相關(guān)（R=0.969）關(guān)系。

為了模擬氣候生產(chǎn)潛力對(duì)氣溫和降水變化的敏感程度，本研究設(shè)置5種氣候情景，將年平均氣溫降低2℃和1℃、氣溫不變、升高1℃和2℃；年降水量減少 20%和10%、降水量不變、增加10%和20%后計(jì)算標(biāo)準(zhǔn)氣候生產(chǎn)潛力，由于研究區(qū)域地處高緯度，太陽(yáng)能豐富，日照充足，基本滿足大豆生長(zhǎng)發(fā)育需要，以多年平均氣溫和降水量計(jì)算氣候生產(chǎn)潛力，比較氣候變化前后黑龍江省自然植被標(biāo)準(zhǔn)氣候生產(chǎn)潛力變化。僅考慮單因子變化，氣溫升高或降低時(shí)，W的變化幅度要比降水大。但實(shí)際上氣溫和降水同時(shí)變化，對(duì)W產(chǎn)生不同影響。由表2可見，當(dāng)氣溫不變，降水增加（或減少）10%、20%，W分別增加（或減少）3.58%、7.14%；當(dāng)降水不變，氣溫升高（或降低）1、2℃，W分別增加（或減少）3.84%、7.67%。因氣溫變化引起的W變化幅度大于降水引起的W變化幅度，說明氣溫主導(dǎo)研究區(qū)域W的大小。當(dāng)氣溫降低1、2℃，降水量增加（或減少）10%、20%時(shí)，W將減小0.26%～14.82%；然而氣溫降低1℃，降水量增加20%時(shí)，W將增加3.31%。當(dāng)氣溫升高1、2℃，降水量增加（或減少）10%、20%時(shí)，W將增加0.26%～14.82%；然而氣溫升高1℃，降水量減少20%時(shí)，W將減少3.31%。與上述分析，氣溫對(duì)W的影響占主導(dǎo)地位一致。研究區(qū)域“冷干型”氣候?qū)?biāo)準(zhǔn)氣候生產(chǎn)潛力不利，“暖濕型”氣候?qū)?biāo)準(zhǔn)氣候生產(chǎn)潛力有利。假設(shè)研究區(qū)域氣候變化的特征趨勢(shì)與上文分析一致，未來將向“暖濕型”發(fā)展，將有利于標(biāo)準(zhǔn)氣候生產(chǎn)潛力發(fā)揮。

3結(jié)論與討論

（1）黑龍江省自然植被氣候生產(chǎn)潛力（WT、WR、WV、W）均呈現(xiàn)增加趨勢(shì)，年代際變化差異明顯，整個(gè)研究時(shí)域水熱配比情況較好，大部分時(shí)期熱量起主要限制作用，降水也是主要限制因子之一。

（2）黑龍江省自然植被氣候生產(chǎn)潛力（WT、WR、WV、W）存在26～31、18～24、11～16年等3種時(shí)間尺度的周期變化，第一主周期分別為28、28、27、28年。

（3）黑龍江省自然植被標(biāo)準(zhǔn)氣候生產(chǎn)潛力空間分布存在3種分布類型，主要表現(xiàn)型式呈現(xiàn)一致的正變化趨勢(shì)，次要表現(xiàn)型式呈現(xiàn)正、負(fù)相間變化；高值中心主要集中在西南部地區(qū)，低值中心主要集中在西北部地區(qū)，自然植被標(biāo)準(zhǔn)氣候生產(chǎn)潛力受地理環(huán)境影響變化較明顯，區(qū)域特征差異顯著；水熱條件匹配較好（較差）區(qū)域的標(biāo)準(zhǔn)氣候生產(chǎn)潛力偏高（偏低）。

（4）黑龍江省自然植被標(biāo)準(zhǔn)氣候生產(chǎn)潛力對(duì)氣象要素變化響應(yīng)敏感，在年際尺度中，氣溫每上升1℃、降水量和蒸散量每增加10 mm，W依次變化為199.91、44.46、233.46 kg/（hm2·a），氣溫是主要限制因子。“冷干型”氣候?qū)?biāo)準(zhǔn)氣候生產(chǎn)潛力不利，“暖濕型”氣候?qū)?biāo)準(zhǔn)氣候生產(chǎn)潛力有利；未來將向“暖濕型”發(fā)展，研究區(qū)域?qū)⒂欣跇?biāo)準(zhǔn)氣候生產(chǎn)潛力發(fā)揮。

（5）Miami模型和Thornthwaite Memorial模型估計(jì)的凈第一性生產(chǎn)力稍大及利用的年平均實(shí)際蒸散量綜合第考慮了氣溫和降水以及下墊面性質(zhì)的影響作用[32]，通過光合作用與植物冠層蒸散伴隨的植物體對(duì)植物生長(zhǎng)機(jī)理以及能量的內(nèi)在轉(zhuǎn)換機(jī)制描述，展現(xiàn)光合作用與植被產(chǎn)量關(guān)系，因此理論上更具可靠性。該模型僅為統(tǒng)計(jì)模型，虧缺緊密的生理、生態(tài)特征及機(jī)理支撐依據(jù)，還需進(jìn)一步完善，增加其他氣象因子如日照、風(fēng)速和相對(duì)濕度等[33]，這將是未來時(shí)期重要任務(wù)之一。

參考文獻(xiàn)

[1]王曉煜，楊曉光，孫爽，等.氣候變化背景下東北三省主要糧食作物產(chǎn)量潛力及資源利用率比較[J].應(yīng)用生態(tài)學(xué)報(bào)，2015，26（10）：3091-3102.

[2]馬雅麗，郭建平，趙俊芳.晉北農(nóng)牧交錯(cuò)帶作物氣候生產(chǎn)潛力分布特征及其對(duì)氣候變化的響應(yīng)[J].生態(tài)學(xué)雜志，2019，38（3）：818-827.

[3]周廣勝.氣候變化對(duì)中國(guó)農(nóng)業(yè)生產(chǎn)影響研究展望[J].氣象與環(huán)境科學(xué)，2015，38（1）：80-94.

[4]郭建平.氣候變化對(duì)中國(guó)農(nóng)業(yè)生產(chǎn)的影響研究進(jìn)展[[J].應(yīng)用氣象學(xué)報(bào)，2015，26（1）：1-11.

[5]高軍波，楚冰洋，閆軍輝，等.1960年以來河南省玉米氣候橫財(cái)潛力估算與種植空間優(yōu)化[J].農(nóng)業(yè)機(jī)械學(xué)報(bào)，2019，50（1）：245-254.

[6]郭建平.氣候變化對(duì)農(nóng)業(yè)氣候資源有效的影響評(píng)估[M].北京：氣象出版社，2016：2-49.

[7]王鵬，王婷，王瑩，等.東北地區(qū)玉米氣候生產(chǎn)潛力評(píng)估及區(qū)劃[J].大麥與谷類科學(xué)，2019，36（1）：43-48，58.

[8]胡慧杰，王猛，尹小剛，等.氣候變化下東北農(nóng)作物大豆需水量的時(shí)空變化特征[J].中國(guó)農(nóng)業(yè)大學(xué)學(xué)報(bào)，2017，22（2）：21-31.

[9]李秀芬，馬樹慶，宮麗娟，等.基于WOFOST的東北地區(qū)玉米生育期氣象條件適宜度評(píng)價(jià)[J].中國(guó)農(nóng)業(yè)氣象，2013，34（1）：43-49.

[10]朱海霞，姜麗霞，曲輝輝，等.黑龍江省水稻生育階段極端降水事件時(shí)空特征[J].氣象，2019，45（4）：522-532.

[11]龍斯玉.江蘇省農(nóng)業(yè)氣候生產(chǎn)潛力及區(qū)劃研究[J].地理科學(xué)，1985， 5（3）：218-226.

[12]陳國(guó)南.用邁阿密模型測(cè)算我國(guó)生物生產(chǎn)量的初步嘗試[J].自然資源學(xué)報(bào)，1987，2（3）：270-278.

[13]高素華，潘亞茹，郭建平.氣候變化對(duì)植物氣候生產(chǎn)力的影響[J].氣象，1994，20（1）：30-33.

[14]王秀芬，楊艷昭，尤飛.黑龍江省氣候變化及其對(duì)玉米生產(chǎn)潛力的影響[J].干旱地區(qū)農(nóng)業(yè)研究，2015，30（5）：25-29.

[15]王宗明，張柏，張樹清，等.松嫩平原農(nóng)業(yè)氣候生產(chǎn)潛力及自然資源利用率研究[J].中國(guó)農(nóng)業(yè)氣象，2005，26（1）：1-6.

[16]李莉，周宏飛，包安明.中亞地區(qū)氣候生產(chǎn)潛力時(shí)空變化特征[J].自然資源學(xué)報(bào)，2014，29（2）：285-294.

[17]于穎，范文義，楊曦光.1901-2018年小興安嶺森林NPP估算[J].林業(yè)科學(xué)，2014，50（10）：16-23.

[18]趙慧穎，田寶星，宮麗娟，等.近308年來大興安嶺北部森林植被氣候生產(chǎn)潛力及其對(duì)氣候變化的響應(yīng)[J].生態(tài)學(xué)報(bào)，2017，37（6）： 1900-1911.

[19]張興林，叢志宇，趙慧穎，等.黑龍江省富裕縣移栽玉米氣候生產(chǎn)潛力分析[J].中國(guó)農(nóng)學(xué)通報(bào)，2016，32（21）：147-152.

[20]盧玢宇，楊波，裴占江，等.黑龍江省氣候生產(chǎn)潛力時(shí)空演變特征研究[J].生態(tài)環(huán)境學(xué)報(bào)，2017，26（10）：1659-1664.

[21]杜國(guó)明，張露洋，徐新良，等.近50年氣候驅(qū)動(dòng)下東北地區(qū)玉米生產(chǎn)潛力時(shí)空演變分析[J].地理研究，2016，35（5）：864-874.

[22]田寶星，于敏，李浩然，等.黑龍江省春季低溫指數(shù)及其對(duì)作物產(chǎn)量的影響[J].中國(guó)農(nóng)學(xué)通報(bào)，2018，34（25）：97-103.

[23]田寶星，宮麗娟，楊帆，等.黑龍江省春季低溫指數(shù)小波分析[J].水土保持研究，2017，24（1）：342-345，350.

[24]趙俊芳，孔祥娜，姜月清，等.基于高時(shí)空分辨率的氣候變化對(duì)全球主要農(nóng)區(qū)氣候生產(chǎn)潛力的影響評(píng)估[J].生態(tài)環(huán)境學(xué)報(bào)，2019，28（1）： 1-6.

[25]羅永忠，成自勇，郭小芹.近40a甘肅省氣候生產(chǎn)潛力時(shí)空變化特征[J].生態(tài)學(xué)報(bào)，2011，31（1）：221-229.

[26]趙慧穎，顧潤(rùn)源，白曉慧，等.氣象水文因子對(duì)呼倫湖濕地區(qū)域植物氣候生產(chǎn)潛力的影響[J].生態(tài)環(huán)境學(xué)報(bào)，2012，21（6）：997-1003.

[27]李夏子，郭春燕，楊晶，等.氣候變化對(duì)內(nèi)蒙古牧草生產(chǎn)潛力的影響[J].中國(guó)農(nóng)學(xué)通報(bào)，2018，34（16）：106-114.

[28]周國(guó)華，羅小莉，王盤興，等.中國(guó)冬季氣溫異常EOF分析的改進(jìn)[J].大氣科學(xué)學(xué)報(bào)，2012，35（3）：295-303.

[29]符傳博，丹利，馮錦明，等.1960-2012年中國(guó)地區(qū)總運(yùn)量時(shí)空變化及其與氣溫和降水的關(guān)系[J].大氣科學(xué)，2019，43（1）：87-98.

[30]宗天韻，周瑋瑩，周平.南嶺山地1968到2015年降水的時(shí)空變化特征研究[J].生態(tài)科學(xué)，2019，38（2）：182-190.

[31]符淙斌，王強(qiáng).氣候突變的定義和檢測(cè)方法[J].大氣科學(xué)，1992，16（4）：482-493.

[32]曹宗龍.安徽省植被潛在生產(chǎn)力研究[J].綿陽(yáng)師范學(xué)院學(xué)報(bào)，2011， 30（2）：77-83.

[33]葛全勝，劉健，方修琦，等.過去200年冷暖變化的基本特征與主要暖期[J].地理學(xué)報(bào)，2013，68（5）：579-592.