盧玢宇，楊波，裴占江，史風(fēng)梅，王粟，高亞冰，劉杰

黑龍江省農(nóng)業(yè)科學(xué)院農(nóng)村能源研究所，黑龍江 哈爾濱 150086

黑龍江省氣候生產(chǎn)潛力時(shí)空演變特征研究

盧玢宇，楊波，裴占江，史風(fēng)梅，王粟，高亞冰，劉杰*

黑龍江省農(nóng)業(yè)科學(xué)院農(nóng)村能源研究所，黑龍江 哈爾濱 150086

氣候生產(chǎn)潛力的時(shí)空變化規(guī)律不僅反映各氣象因子與氣候生產(chǎn)潛力之間的配合協(xié)調(diào)程度，還能對(duì)糧食生產(chǎn)決策起到?jīng)Q定性的指導(dǎo)意義?；赥hornthwaite Memorial模型估算了黑龍江省1986—2015年作物生育期內(nèi)的氣候生產(chǎn)潛力，分析了氣候生產(chǎn)潛力與平均溫度和降水量之間的關(guān)系。采用 Morlet小波分析方法，對(duì)黑龍江省作物生育期內(nèi)的氣候生產(chǎn)潛力進(jìn)行了多時(shí)間尺度特征分析和預(yù)測(cè)，并通過GIS軟件擬合分析了黑龍江省作物生育期內(nèi)氣候生產(chǎn)潛力的空間變化特征。結(jié)果表明，1986—2015年作物生育期內(nèi)黑龍江省氣候生產(chǎn)潛力呈緩慢的降低趨勢(shì)，變化趨勢(shì)不顯著，在7012.143～11680.771 kg·hm-2·a-1之間變化，平均值為9390.362 kg·hm-2·a-1；1986—2015年生育期內(nèi)黑龍江省氣候生產(chǎn)潛力具有顯著的多時(shí)間尺度特征，存在7 a和10 a的主周期以及16 a的變化周期；擬合可知2016—2018年的氣候生產(chǎn)潛力逐年遞減，但均高于年平均值；黑龍江省生育期內(nèi)氣候生產(chǎn)潛力空間分布差異顯著，以巴彥縣為中心，氣候生產(chǎn)潛力向四周逐漸遞減，體現(xiàn)出較為明顯的經(jīng)度地帶性，大部分地區(qū)熱量條件相對(duì)充足，降水量是影響黑龍江省作物氣候生產(chǎn)潛力的主要限制因子。

氣候生產(chǎn)潛力；時(shí)空演變；Thornthwaite Memorial 模型；小波分析；黑龍江省

氣候生產(chǎn)潛力是指在一定的光、溫、水資源條件下，其他的環(huán)境因子（CO2、養(yǎng)分等）和作物群體因素處于最適宜狀態(tài)，作物利用當(dāng)?shù)氐墓?、溫、水資源的潛在生產(chǎn)力（逯亞杰，2016；張藝萌等，2015；李莉等，2014；羅永忠等，2011；侯西勇，2008）。氣候生產(chǎn)潛力模型不僅能夠反映某一區(qū)域氣象因子間的協(xié)同作用效果，而且還能為提高農(nóng)事操作措施提供建議（張銳等，2012；武永利等，2009；趙慧穎等，2008）。不同地區(qū)地貌特征、環(huán)境因子時(shí)空差異較大，因而其氣候生產(chǎn)潛力的時(shí)空變化也會(huì)表現(xiàn)出一定的差異。氣候生產(chǎn)潛力的時(shí)空動(dòng)態(tài)變化一直是國(guó)內(nèi)外研究的熱點(diǎn)，近年來，雖然對(duì)黑龍江省氣候生產(chǎn)潛力的時(shí)空分布已有研究（于盛楠等，2005；祖世亨等，2001），但對(duì)其周期性變化和趨勢(shì)變化的分析和預(yù)測(cè)還鮮有報(bào)道。

黑龍江省位于中國(guó)東北最北部地區(qū)，面積約為45.4×104km2，東西分布有三江平原和松嫩平原，土地肥沃，物產(chǎn)豐富，盛產(chǎn)各種農(nóng)副產(chǎn)品，全省耕地面積約為1.19×107hm2，是中國(guó)最大的商品糧生產(chǎn)基地。然而，由于黑龍江省四季氣候變化差異顯著，其農(nóng)業(yè)生態(tài)系統(tǒng)受其影響而波動(dòng)性強(qiáng)（李天霄等，2017；袁彬，2012；王鳳玲等，2011），故研究其氣候生產(chǎn)潛力對(duì)農(nóng)業(yè)生產(chǎn)意義重大。本研究以黑龍江省為研究區(qū)域，采用Thornthwaite Memorial模型估算了黑龍江省作物生育期（5—9月）的氣候生產(chǎn)潛力，揭示了 1986—2015年黑龍江省生育期氣候生產(chǎn)潛力的時(shí)空演變特征，并預(yù)測(cè)了生育期氣候生產(chǎn)潛力的變化趨勢(shì)，以期為黑龍江省的農(nóng)事操作提供基礎(chǔ)數(shù)據(jù)支持。

1 材料與研究方法

1.1 數(shù)據(jù)來源

選取黑龍江省內(nèi) 31個(gè)氣象站的月平均氣溫和降水資料，時(shí)間序列為 1986—2015年，氣象數(shù)據(jù)來源于黑龍江省氣象局地面氣候資料。站點(diǎn)選取考慮了時(shí)間序列的連續(xù)性、完整性，地域的代表性和空間分布情況等因素。所選氣象站點(diǎn)見圖1。

圖1 1986—2015年黑龍江省氣候生產(chǎn)潛力空間分布特征Fig. 1 Spatial distribution of climatic productivity potential in Heilongjiang Province from 1986 to 2015

1.2 氣候生產(chǎn)潛力模型

Thornthwaite Memorial模型（趙慧穎等，2012；于盛楠等，2005；Lieth，1986；Lieth，1975）可以定量表征作物產(chǎn)量與溫度、降水量和蒸發(fā)量之間的關(guān)系。其計(jì)算公式如下：

式中，Wv為計(jì)算獲得的作物氣候生產(chǎn)潛力（kg·hm-2·a-1）；V 為生育期平均蒸散量（mm）；e為自然對(duì)數(shù)底數(shù)。

式中，R為生育期降水量（mm）；L為生育期平均最大蒸散量。

式中，T為生育期平均溫度℃。僅當(dāng)R＞0.316 L時(shí)，式（3）才適用，而當(dāng)R≤0.316 L時(shí)，V=R（姚玉璧等，2006；梁瑞龍等，1998；高素華等，1994）。

1.3 小波分析

小波分析（Wavelet Analysis）方法是基于傅立葉變換，在時(shí)域和頻域方面都能夠表征信號(hào)的局部特征，并且對(duì)信號(hào)進(jìn)行多尺度細(xì)化分析，得到所有頻率隨時(shí)間的變化特征和不同頻率之間的關(guān)系（Ramírez et al.，2016；Yang et al.，2016；Zhang et al.，2007）。

本文采用Morlet小波（Morlet wavelet）函數(shù)對(duì)黑龍江省的氣候生產(chǎn)潛力進(jìn)行研究。Morlet小波具有很好的時(shí)頻局部性，多用于復(fù)制信號(hào)的分解和時(shí)頻分析（孫春媛等，2016；Huang et al.，2015；任艷林，2012）。Morlet小波函數(shù)為非正交小波，定義為：

式中，Ψ(t)為小波函數(shù)；ω0為小波中心頻率；t為時(shí)間（a）。

小波方差圖能夠反映時(shí)間序列中的周期波動(dòng)及其強(qiáng)弱隨尺度變化的特性，可以表征在一個(gè)時(shí)間序列中起主要作用的時(shí)間尺度，即主周期（張佩等，2009；李春強(qiáng)等，2010），其計(jì)算公式如下：

1.4 數(shù)據(jù)處理

氣候生產(chǎn)潛力和氣象因子的區(qū)域平均值采用算數(shù)平均法計(jì)算，運(yùn)用典型相關(guān)分析方法，對(duì)氣候生產(chǎn)潛力與平均溫度和降雨量進(jìn)行相關(guān)性分析，并在α=0.01水平上進(jìn)行顯著性檢驗(yàn)（陳廣洲等，2017；于雷等，2015）；利用GIS空間插值法對(duì)研究區(qū)域各因子進(jìn)行插值分析，探討其空間分布特征（臧星華等，2015；張藝萌等，2015）；通過小波分析識(shí)別和對(duì)比氣候生產(chǎn)潛力的多時(shí)間尺度特征，分析不同時(shí)間尺度下氣候生產(chǎn)潛力的演變規(guī)律和發(fā)展趨勢(shì)（侯西勇，2008）；采用Origin 2015軟件進(jìn)行繪圖。

2 結(jié)果與分析

2.1 黑龍江省氣候生產(chǎn)潛力空間分布特征

如圖1所示，黑龍江省作物生育期內(nèi)氣候生產(chǎn)潛力空間分布地理差異顯著，規(guī)律明顯。以巴彥縣為中心，氣候生產(chǎn)潛力向四周逐漸遞減，體現(xiàn)出較明顯的經(jīng)度地帶性，不同經(jīng)度的氣候生產(chǎn)潛力水平差別較大。黑龍江省西部和東部部分地區(qū)氣候生產(chǎn)潛力水平比較低，嫩江市、富錦市、安達(dá)縣、富裕縣、寶清縣、呼瑪縣、肇源縣、泰來縣等地氣候生產(chǎn)潛力均低于 9000 kg·hm-2·a-1；中部地區(qū)氣候生產(chǎn)潛力的水平比較高，其中，巴彥縣、尚志市、五常市、海倫市、方正縣等地氣候生產(chǎn)潛力都多年高于10000 kg·hm-2·a-1。

近 30年來，泰來縣作物生育期的平均氣候生產(chǎn)潛力值最低，為 8045.801 kg·hm-2·a-1，其氣候生產(chǎn)潛力與降雨量之間呈線性正相關(guān)，r值為0.994，通過了0.01水平上的顯著性檢驗(yàn)（n=30），氣候生產(chǎn)潛力與平均溫度之間呈負(fù)相關(guān)，r值為-0.422，在0.01水平上無顯著性相關(guān)（n=30），說明泰來縣的降雨量對(duì)氣候生產(chǎn)潛力的影響大于平均溫度。泰來縣近 30年來的年平均溫度高于其他地區(qū)，熱量條件充足，但由于該地區(qū)處于干旱半干旱地區(qū)，沙地面積大，平均降雨量處于較低水平，土地退化嚴(yán)重，生態(tài)環(huán)境脆弱，嚴(yán)重限制了該區(qū)氣候生產(chǎn)潛力的發(fā)展（張雪萍等，2005）。

巴彥縣氣候生產(chǎn)潛力水平最高，為 10596.743 kg·hm-2·a-1，在近30年時(shí)間里，巴彥縣在作物生育期內(nèi)的降雨量居于全省首位，平均溫度處于中等水平（見圖2），雨量充沛，有利于作物生長(zhǎng)，促進(jìn)了該區(qū)氣候生產(chǎn)潛力的發(fā)展。

尚志市的氣候生產(chǎn)潛力為 10496.112 kg·hm-2·a-1，位居第二，降雨量緊隨巴彥縣降雨量之后。李慶等（2008）對(duì)尚志市的農(nóng)業(yè)可持續(xù)性進(jìn)行了評(píng)價(jià)，表明尚志市作物生長(zhǎng)季降水量占全年降水量的87%，雨熱同季，水熱協(xié)同作用較好，利于作物生長(zhǎng)，促進(jìn)了該區(qū)氣候生產(chǎn)潛力的發(fā)展。

綜上，氣候生產(chǎn)潛力較高的地區(qū)是降水量較多的巴彥縣、尚志市等地區(qū)，肇源縣、泰來縣以及呼瑪縣等地區(qū)降水較少，生產(chǎn)潛力普遍偏低。

分別計(jì)算1986—2015年黑龍江省31個(gè)氣象站在生育期內(nèi)同一水平條件下氣候生產(chǎn)潛力與同期、平均氣溫和降水量的相關(guān)系數(shù)（r），結(jié)果發(fā)現(xiàn)，氣候生產(chǎn)潛力與溫度的 r值在-0.467～0.330之間，顯著性檢驗(yàn)表明，除龍江縣（P=-0.467，n=30）、呼瑪縣（P=-0.403，n=30）、嫩江縣（P=-0.457，n=30）、泰來縣（P=-0.422，n=30）的氣候生產(chǎn)潛力與溫度之間的P值均在0.05水平上顯著相關(guān)外，其余地區(qū)間均不存在顯著相關(guān)性；降水量與氣候生產(chǎn)潛力的r值在0.980～0.996之間，均在0.01水平上顯著相關(guān)（n=30）。由此表明，黑龍江省大部分地區(qū)生育期內(nèi)熱量條件相對(duì)充足，但溫度對(duì)作物氣候生產(chǎn)潛力的影響不顯著，降水量是限制該區(qū)域作物氣候生產(chǎn)潛力的主要?dú)庀笠蜃印?/p>

2.2 黑龍江省氣候生產(chǎn)潛力時(shí)間變化特征

圖2 1986—2015年不同地區(qū)平均溫度和平均降雨量Fig. 2 Average temperature and rainfall from 1986 to 2015 for each monitoring region

黑龍江省氣候生產(chǎn)潛力在 1986—2015年的時(shí)間變化特征如圖3所示，生育期內(nèi)總氣候生產(chǎn)潛力在 7012.143～11680.771 kg·hm-2·a-1之間變化，平均值為 9390.362 kg·hm-2·a-1，年際間波動(dòng)較大。1986—2015年的氣候生產(chǎn)潛力變化呈緩慢降低趨勢(shì)，降低趨勢(shì)變化不顯著。其中，2001年黑龍江省的氣候生產(chǎn)潛力值最低，2001年的氣候特點(diǎn)為降雨偏少，作物生育期內(nèi)出現(xiàn)歷史上最為嚴(yán)重的春夏連旱，6月全省大部地區(qū)降水偏少，尤其上旬中西部大部、下旬大部農(nóng)區(qū)幾乎無降水，西部部分地區(qū)最高氣溫突破歷史極值。截至6月末，除牡丹江外其他大部地區(qū)均呈現(xiàn)不同程度旱象。雖然生育期內(nèi)全省大部分地區(qū)平均溫度比常年偏高1～2 ℃，但作物生長(zhǎng)期的降水條件較差，限制了氣候生產(chǎn)潛力的發(fā)展。1994年黑龍江省氣候生產(chǎn)潛力值最高，1994年的氣候特點(diǎn)為高溫多降水，氣溫、降水等要素均出現(xiàn)歷史極值。生育期內(nèi)全省平均氣溫比歷年偏高1.5 ℃左右，水熱同步，資源充足，極大地促進(jìn)了氣候生產(chǎn)潛力的發(fā)展。

圖3 1986—2015年黑龍江省氣候生產(chǎn)潛力年際變化Fig. 3 Interannual variability of climatic productivity potential in Heilongjiang Province from 1986 to 2015 n=31

2.3 黑龍江省氣候生產(chǎn)潛力周期變化特征

通過對(duì)近30年來黑龍江省氣候生產(chǎn)潛力的統(tǒng)計(jì)分析，發(fā)現(xiàn)生育期內(nèi)氣候生產(chǎn)潛力具有較為明顯的多時(shí)間尺度特征，如圖4所示。其中在9～10 a尺度周期表現(xiàn)出正負(fù)相位交替出現(xiàn)的現(xiàn)象，較強(qiáng)周期振蕩幾乎存在于整個(gè)研究時(shí)域內(nèi)，周期性表現(xiàn)十分顯著；在6～7 a尺度周期表現(xiàn)出正負(fù)相位交替出現(xiàn)的現(xiàn)象，較強(qiáng)周期振蕩幾乎存在于整個(gè)研究時(shí)域內(nèi)，周期性表現(xiàn)十分顯著；在11～16 a尺度周期上表現(xiàn)出負(fù)正相位交替出現(xiàn)的現(xiàn)象，周期性振蕩表現(xiàn)得比較明顯。

如圖5所示，小波方差圖有2個(gè)相對(duì)明顯的峰值，分別對(duì)應(yīng)7 a和10 a的時(shí)間尺度，均有一定的波動(dòng)能量，而在7 a以下時(shí)間尺度、8～9 a以及11 a的時(shí)間尺度上，波動(dòng)能量都比較微弱，16 a時(shí)間尺度上存在的周期性變化，仍需要更長(zhǎng)時(shí)間的驗(yàn)證，因此認(rèn)為，松花江流域氣候生產(chǎn)潛力在中小尺度上的變化有1個(gè)7 a和10 a的主周期，另外還有1個(gè)潛在的16 a中長(zhǎng)尺度變化周期。

圖4 1986—2015年黑龍江省氣候生產(chǎn)潛力小波系數(shù)變化Fig. 4 Variation of Morlet wavelet coefficients of climatic productivity potential in Heilongjiang Province from 1986 to 2015

圖5 198—2015年氣候生產(chǎn)潛力Morlet小波方差圖Fig. 5 Morlet wavelet variance of climatic productivity potential from 1986 to 2015

2.4 氣候生產(chǎn)潛力典型尺度的趨勢(shì)預(yù)測(cè)

圖6 氣候生產(chǎn)潛力主周期尺度小波系數(shù)圖Fig. 6 Morlet wavelet coefficients on real component of major cycle scale of climatic productivity potential

典型尺度（主周期）的氣候生產(chǎn)潛力過程線可以揭示小波分析的多尺度變化，也能體現(xiàn)其預(yù)測(cè)趨勢(shì)。如圖6所示，黑龍江氣候生產(chǎn)潛力在7 a時(shí)間尺度上經(jīng)歷了3個(gè)正負(fù)相位的波動(dòng)，平均變化周期為9 a；在10 a時(shí)間尺度上經(jīng)歷了2個(gè)正負(fù)相位的波動(dòng)，平均變化周期為12 a；在16 a時(shí)間尺度上經(jīng)歷了1個(gè)正負(fù)相位的波動(dòng)，平均變化周期為20 a。3個(gè)主周期7 a、10 a、16 a疊加的總過程線，代表著 1986—2015年黑龍江省松花江流域氣候生產(chǎn)潛力在總時(shí)間尺度的變化情況，依次經(jīng)歷了 1986—1996年10年的正向位，1997—2008年11年的負(fù)向位，以及2009—2015年6年的正向位。由此可知，1986—2015年，黑龍江省松花江流域氣候生產(chǎn)潛力的波動(dòng)主要受到7 a、10 a和16 a 3個(gè)主周期時(shí)間尺度的影響，其擬合公式如下所示：

式（6）可以預(yù)測(cè)未來幾年黑龍江省作物生育期內(nèi)氣候生產(chǎn)潛力變化情況。計(jì)算可知，2016—2018年的氣候生產(chǎn)潛力呈現(xiàn)逐年遞減趨勢(shì)，但高于年平均值。

3 討論

黑龍江省作為中國(guó)重要的商品糧生產(chǎn)基地之一，分析并預(yù)測(cè)其氣候變化規(guī)律具有重要的現(xiàn)實(shí)意義。黑龍江省大部分地區(qū)在生育期內(nèi)雖然熱量條件相對(duì)充足，但溫度并不是影響作物氣候生產(chǎn)潛力的主要因素，降水量才是限制該區(qū)域作物氣候生產(chǎn)潛力的主要?dú)庀笠蜃?。干旱缺水使豐富的土地資源不能得到充分利用，農(nóng)作物生長(zhǎng)發(fā)育受到限制，從而導(dǎo)致其氣候生產(chǎn)潛力水平偏低，這與祖世亨等（2001）的研究結(jié)果一致。若積極實(shí)施人工增雨和引流灌溉等工程，使其光熱資源與水資源得到最優(yōu)配置，預(yù)計(jì)未來農(nóng)作物氣候產(chǎn)量會(huì)有很大的提升空間。對(duì)黑龍江省作物生育期內(nèi)的氣候生產(chǎn)潛力進(jìn)行時(shí)間尺度分析，發(fā)現(xiàn)在16 a時(shí)間尺度上存在1個(gè)正負(fù)相位波動(dòng)的周期性變化。如果選擇的時(shí)間尺度再大一些，16 a時(shí)間尺度上氣候生產(chǎn)潛力的周期性變化就能得到更好的體現(xiàn)。本文對(duì)黑龍江省作物生育期內(nèi)的氣候生產(chǎn)潛力進(jìn)行了周期變化和趨勢(shì)變化的分析和預(yù)測(cè)，這對(duì)于干旱和半干旱地區(qū)作物氣候生產(chǎn)潛力的研究具有一定的參考意義。由于現(xiàn)有數(shù)據(jù)的時(shí)間尺度稍小，一定程度上限制了小波分析功能和優(yōu)勢(shì)的發(fā)揮。

4 結(jié)論

（1）黑龍江省作物生育期內(nèi)的氣候生產(chǎn)潛力以巴彥縣為中心，氣候生產(chǎn)潛力向四周逐漸遞減，體現(xiàn)出較為明顯的經(jīng)度地帶性，不同經(jīng)度的氣候生產(chǎn)潛力水平差別較大。黑龍江省大部分地區(qū)熱量條件相對(duì)充足，降水量是限制該區(qū)域作物氣候生產(chǎn)潛力的主要?dú)庀笠蜃印?/p>

（2）1986—2015年黑龍江省作物生育期內(nèi)的氣候生產(chǎn)潛力在 7012.143～11680.771 kg·hm-2·a-1之間變化，年際間波動(dòng)較大，總體呈緩慢降低趨勢(shì)，但降低趨勢(shì)變化不顯著。

（3）通過對(duì)近 30年黑龍江省作物生育期內(nèi)的氣候生產(chǎn)潛力進(jìn)行分析，發(fā)現(xiàn)氣候生產(chǎn)潛力在中小尺度上的變化有1個(gè)7 a和10 a的主周期，還有1個(gè)潛在 16 a的中長(zhǎng)尺度變化周期。通過擬合預(yù)測(cè)2016—2018年的氣候生產(chǎn)潛力呈現(xiàn)逐年遞減趨勢(shì)，但高于年平均值。

HUANG P, ZHANG J, TANG Y, et al. 2015. Spatial and temporal distribution of PM2.5pollution in Xi’an City, China [J]. International journal of environmental research and public health, 12(6): 6608-6625.LIETH H. 1975. Modeling the primary productivity of the world [M]. New York: Springer: 237-263.

LIETH H. 1986. 生物圈的第一性生產(chǎn)力[M]. 王業(yè)蔻等譯. 北京: 科學(xué)出版社.

RAMíREZ I J, GRADY S C. 2016. El Ni?o, climate, and cholera associations in Piura, Peru, 1991—2001: a wavelet analysis [J].Ecohealth, 13(1): 83-99.

YANG L, CAI X J, ZHANG H, et al. 2016. Interdependence of foreign exchange markets: A wavelet coherence analysis [J]. Economic Modelling, 55: 6-14.

ZHANG Z, ZHANG Q, JIANG T. 2007. Changing features of extreme precipitation in the Yangtze River basin during 1961—2002 [J]. Journal of Geographical Sciences, 17(1): 33-42.

陳廣洲, 李鑫海, 焦利鋒, 等. 2017. 2000—2012年淮南煤礦區(qū)植被凈初級(jí)生產(chǎn)力的時(shí)空變化特征[J]. 生態(tài)環(huán)境學(xué)報(bào), 26(2): 196-203.

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

侯西勇. 2008. 1951—2000年中國(guó)氣候生產(chǎn)潛力時(shí)空動(dòng)態(tài)特征[J]. 干旱區(qū)地理, 31(5): 723-730.

李春強(qiáng), 杜毅光, 李保國(guó). 2010. 1965—2005 年河北省降水量變化的小波分析[J]. 地理科學(xué)進(jìn)展, 29(11): 1340-1344.

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

李慶, 尚杰, 于法穩(wěn). 2008. 黑龍江省尚志市農(nóng)業(yè)可持續(xù)性評(píng)價(jià)[J]. 綠色經(jīng)濟(jì), (4): 91-94.

李天霄, 付強(qiáng), 孟凡香, 等. 2017. 黑龍江省降水變化趨勢(shì)及其對(duì)農(nóng)業(yè)生產(chǎn)的影響研究[J]. 灌溉排水學(xué)報(bào), 36(5): 103-108.

梁瑞龍, 蔣炎炎, 朱積余. 1998. 廣西植被潛在生產(chǎn)力的估算[J]. 廣西林業(yè)科學(xué), 27(2): 68-74.

逯亞杰. 2016. 黃土高原地區(qū)氣候生產(chǎn)潛力估算及其對(duì)氣候變化的響應(yīng)[D]. 咸陽: 西北農(nóng)林科技大學(xué): 1-2.

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

任艷林. 2012. 1965—2011 年河北塞罕壩地區(qū)降水量變化規(guī)律的小波分析[J]. 北京大學(xué)學(xué)報(bào)(自然科學(xué)版), 48(6): 918-924.

孫春媛, 李令軍, 趙文吉, 等. 2016. 基于小波變換的北京市 PM2.5時(shí)空分布特征及成因分析[J]. 生態(tài)環(huán)境學(xué)報(bào), 25(8): 1343-1350.

王鳳玲, 王曉明. 2011. 氣候變化對(duì)黑龍江省農(nóng)業(yè)的影響[J]. 黑龍江氣象, 28(4): 24.

武永利, 盧淑賢, 王云峰, 等. 2009. 近 45年山西省氣候生產(chǎn)潛力時(shí)空變化特征分析[J]. 生態(tài)環(huán)境學(xué)報(bào), 18(2): 567-571.

楊重一, 龐士力, 孫彥坤. 2010. 黑龍江省作物氣候生產(chǎn)潛力估算[J]. 東北農(nóng)業(yè)大學(xué)學(xué)報(bào), 41(3): 75-78.

姚玉璧, 王毅榮, 張存杰, 等. 2006. 黃土高原作物氣候生產(chǎn)力對(duì)氣候變化的響應(yīng)[J]. 南京氣象學(xué)院學(xué)報(bào), 29(1): 101-106.

于雷, 洪永勝, 耿雷, 等. 2015. 基于偏最小二乘回歸的土壤有機(jī)質(zhì)含量高光譜估算[J]. 農(nóng)業(yè)工程學(xué)報(bào), 31(14): 103-109.

于盛楠, 鄒秋菊. 2005. 用李斯公式分析氣候變化對(duì)作物氣候生產(chǎn)潛力的影響[J]. 黑龍江氣象, (2): 11-13.

袁彬. 2012. 氣候變化下東北春玉米氣候生產(chǎn)潛力及農(nóng)業(yè)氣候資源利用率[D]. 北京: 中國(guó)氣象科學(xué)研究院: 1-12.

臧星華, 魯垠濤, 姚宏, 等. 2015. 中國(guó)主要大氣污染物的時(shí)空分布特征研究[J]. 生態(tài)環(huán)境學(xué)報(bào), 24(8): 1322-1329.

張佩, 聶義波. 2009. 小波方差分析伊河流域降雨徑流周期特征[J]. 水資源與水工程學(xué)報(bào), 20(2): 156-158.

張銳, 劉普幸, 張克新. 2012. 新疆草地氣候生產(chǎn)潛力變化特征及對(duì)氣候響應(yīng)的預(yù)測(cè)研究[J]. 中國(guó)沙漠, 32(1): 181-187.

張憲洲. 1993. 我國(guó)自然植被第一性生產(chǎn)力的估算與分布[J]. 自然資源,(1): 15-21.

張雪萍, 郭艷清, 高梅香, 等. 2005. 黑龍江省西部沙地生態(tài)經(jīng)濟(jì)系統(tǒng)能值分析——以泰來縣為例[J]. 經(jīng)濟(jì)地理, 25(5): 651-654.

張藝萌, 張雪松, 郭婷婷, 等. 2015. 遼西北地區(qū)氣溫和降水變化對(duì)氣候生產(chǎn)潛力的影響[J]. 中國(guó)農(nóng)業(yè)氣象, 36(2): 203-211.

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

趙慧穎, 魏學(xué)占, 烏秋力, 等. 2008. 呼倫貝爾典型草原區(qū)牧草氣候生產(chǎn)潛力評(píng)估[J]. 干旱地區(qū)農(nóng)業(yè)研究, 26(1): 137-140.

祖世亨, 魏松林. 2001. 黑龍江省氣候生產(chǎn)潛力研究[J]. 黑龍江氣象, (2):35-39.

Spatial and Temporal Evolution of Climatic Productivity Potential in Heilongjiang Province

LU Binyu, YANG Bo, PEI Zhanjiang, SHI Fengmei, WANG Su, GAO Yabing, LIU Jie
Rural Energy Institute of Heilongjiang Academy of Agricultural Sciences, Harbin 150086, China

Heilongjiang region is one of the important commodity grain production bases in China. The rules of temporal and spatial variation of climate productivity potential not only reflect the coordination between meteorological factors and climatic potential productivity, but also play a decisive guiding role in the decision-making of grain production in this region. Climatic productivity potential during the growth period in Heilongjiang Province from 1986 to 2015 was estimated based on the Thornthwaite Memorial model and relationship between climatic productivity potential and temperature and precipitation was analyzed. Morlet wavelet analysis method was used to analyze and forecast the climatic productivity potential during the growth period in multi-time scales,and GIS software was applied to fit the spatial variation characteristics of climatic productivity potential in life circle of crops.Climatic productivity potential during the growth period in Heilongjiang Province decreased slowly from 1986 to 2015, and varied between 7012.143～11680.771 kg·hm-2·a-1, with an average of 9390.362 kg·hm-2·a-1. Precipitation was the main limiting factor that affected the climatic productivity potential at growth stage of crops in Heilongjiang Province. Climatic productivity potential at growth stage of crops in Heilongjiang Province from 1986 to 2015 had characteristics of significant multi-time scale. The main cycle scales were 7 years and 10 years, and with a long period of 16 years. It is predicted that climatic productivity potential from 2016 to 2018 will be declining year by year, but it will be higher than the annual average. Spatial distribution of climatic productivity potential in Heilongjiang Province was significantly different. The climatic productivity potential at growth stage of crops was centered on Bayan County, and gradually diminished outward, showing distinct characteristics of longitude. The heat condition in most parts of Heilongjiang Province was relative sufficient, and precipitation was the main limiting meteorological factor affecting the climatic productivity potential.

climatic productivity potential; spatial and temporal evolution; thornthwaite memorial model; wavelet analysis;Heilongjiang Province

10.16258/j.cnki.1674-5906.2017.10.003

S162.5; X16

A

1674-5906（2017）10-1659-06

盧玢宇, 楊波, 裴占江, 史風(fēng)梅, 王粟, 高亞冰, 劉杰. 2017. 黑龍江省氣候生產(chǎn)潛力時(shí)空演變特征研究[J]. 生態(tài)環(huán)境學(xué)報(bào), 26(10): 1659-1664.

LU Binyu, YANG Bo, PEI Zhanjiang, SHI Fengmei, WANG Su, GAO Yabing, LIU Jie. 2017. Spatial and temporal evolution of climatic productivity potential in Heilongjiang Province [J]. Ecology and Environmental Sciences, 26(10): 1659-1664.

中國(guó)清潔發(fā)展機(jī)制基金項(xiàng)目（2014101）；國(guó)家“十三五”重點(diǎn)研發(fā)計(jì)劃課題（2016YFD0501403）；國(guó)家科技支撐計(jì)劃課題（2015BAD21B004）；黑龍江省農(nóng)業(yè)科學(xué)院院級(jí)課題（2017ZC07）

盧玢宇（1983年生），女，助理研究員，博士，從事氣候變化及環(huán)境污染方面的研究。E-mail: lubinyu2004@163.com*通信作者：劉杰（1974年生），男，研究員，博士，從事生態(tài)農(nóng)業(yè)和氣候變化方面的研究。E-mail: Liujie1677@126.com

2017-08-04