邢蘭芹, 姚眾眾

(1.陜西學(xué)前師范學(xué)院 環(huán)境與資源管理系, 陜西 西安 710100; 2.南京信息工程大學(xué) 地理與遙感學(xué)院, 江蘇 南京 210044)

氣候生產(chǎn)潛力是指當(dāng)其他條件(如土壤、養(yǎng)分、二氧化碳等)處于最適狀況時，充分和合理利用當(dāng)?shù)氐墓狻帷⑺畾夂蛸Y源，單位面積土地上可能獲得的最高生物學(xué)產(chǎn)量或農(nóng)業(yè)產(chǎn)量[1-5]。由于氣候資源的差異，氣候生產(chǎn)潛力在區(qū)域間存在較大差異，其時空動態(tài)一直是國內(nèi)外全球變化研究的熱點(diǎn)之一。近年來，國內(nèi)外對氣候生產(chǎn)潛力的研究己經(jīng)取得一定成果[6-9]，并建立了一些比較成熟的分析模型[10-13]，對氣候生產(chǎn)潛力的探究提供了重要的理論基礎(chǔ)。

江蘇省是中國糧食的主要生產(chǎn)區(qū)，糧食產(chǎn)量直接影響中國糧食市場的供求關(guān)系，隨著人口的增長，耕地面積的不斷下降以及全球氣候的變化，糧食安全問題日益突出。而如何在有限的耕地面積下充分合理的利用氣候資源，提高單位面積的作物產(chǎn)量，成為解決這一問題的關(guān)鍵。所以，探究氣候生產(chǎn)潛力變化趨勢和主要影響因子，充分發(fā)揮氣候生產(chǎn)潛力，對提高農(nóng)業(yè)生產(chǎn)力，指導(dǎo)農(nóng)業(yè)生產(chǎn)，具有重要的理論和現(xiàn)實(shí)指導(dǎo)意義?；诖耍狙芯繑M利用1971—2010年江蘇省氣象數(shù)據(jù)，應(yīng)用Miami模型、Thornthwaite Memorial模型，探究江蘇省40 a的溫度生產(chǎn)潛力、降水生產(chǎn)潛力和氣候生產(chǎn)潛力的變化規(guī)律，并對氣候資源利用率及氣候生產(chǎn)潛力的影響因素進(jìn)行剖析，以期為區(qū)域農(nóng)業(yè)生產(chǎn)提供理論依據(jù)。

1 研究區(qū)概況和研究方法

1.1 研究區(qū)概況

江蘇省位于中國內(nèi)地的東部區(qū)域的中部，處于淮河、長江下游，江蘇省東部臨近黃海，北部接山東，西部連接安徽，東南與浙江省、上海市交界是長江三角洲的組成部分，江蘇省位于東經(jīng)120°57′—116°18′，北緯35°20′—30°45′。江蘇省海岸線全長為955 km，面積約為1.073×105km2，占中國總面積的1.07%，長江貫穿江蘇省東西426 km，京杭大運(yùn)河貫穿南北719 km。

在地形上，江蘇省以平原為主，主要由長江中下游、蘇北等平原組成，地勢就全國而言是最低的一個省區(qū)，平原面積大約7×104km2，約占江蘇省面積的70%，比例占據(jù)中國各省份中的首位。在氣候上，江蘇省處在由溫帶向亞熱帶的過度性氣候帶，溫暖濕潤，四季分明，并以淮河為界，在淮河以北是暖溫帶濕潤、半濕潤的季風(fēng)氣候，在淮河以南則是亞熱帶濕潤季風(fēng)氣候。

本文選取1971—2010年江蘇省域范圍內(nèi)具有空間代表性的13個氣象站點(diǎn)觀測月平均氣溫、月平均降水?dāng)?shù)據(jù)作為數(shù)據(jù)源，站點(diǎn)分別為徐州市、贛榆區(qū)、常州、東臺、高郵、淮陰、溧陽市、呂泗鎮(zhèn)、南京、南通市、射陽縣、吳縣東山鎮(zhèn)、盱眙縣。氣溫和降水?dāng)?shù)據(jù)來源于中國氣象數(shù)據(jù)網(wǎng)。糧食播種面積和糧食產(chǎn)量來源于江蘇省統(tǒng)計(jì)年鑒。

1.2 研究方法

1.2.1 Miami模型 溫度、降水是制約作物生長的重要?dú)夂蛞蜃樱瑴囟群徒邓畣为?dú)或復(fù)合會影響到作物生產(chǎn)，分析溫度和降水分別對生產(chǎn)力的影響，采用Miami模型來分別計(jì)算氣溫和降水的生產(chǎn)潛力公式為：

Yt=3 000/(1+e1.315-0.119t)

(1)

Yr=3 000/(1-e0.000 664r)

(2)

式中：t——年平均溫度(℃)；r——年降水量(mm)；Yt，Yr——年平均溫度、年降水量所決定的生產(chǎn)潛力〔kg/(hm2·a)〕； e=2.718 3。下同。

1.2.2 Thornthwaite Memorial模型 該模型依據(jù)蒸散量與氣溫、降水和作物的關(guān)系，計(jì)算氣候生產(chǎn)潛力。公式為:

Ye=3 000〔1-e-0.0 009 695(V-20)〕

(3)

V=1.05R/〔1+(1.05R/L)2〕1/2

(4)

L=300+25T+0.05T3

(5)

式中：V——年平均蒸散量(mm)；R——年均降水(mm)；L——年平均最大蒸發(fā)量(mm)；Ye——利用蒸散量計(jì)算的生產(chǎn)潛力〔kg/(hm2·a)〕。

1.2.3 氣候資源利用率公式 參照吳宜進(jìn)[14]的方法，定義氣候資源利用率為實(shí)際糧食產(chǎn)量與生產(chǎn)力之比，即：

(6)

2 結(jié)果與分析

2.1 江蘇省年均氣溫和降水變化分析

由圖1可見，1971—2010年江蘇省年均氣溫呈波動性上升趨勢。累積年平均氣溫為14.86 ℃，年均最高溫出現(xiàn)在2007年，年均氣溫為16.16 ℃，年均最低溫出現(xiàn)在1978年，為13.95 ℃，二者差值為2.21 ℃。由此可知，在全球氣候變暖大背景下，江蘇省的年均氣溫產(chǎn)生了升溫響應(yīng)。

圖1 1971-2010年江蘇省年均氣溫和年均降水量變化

1971—2010年江蘇省年均降水波動較大，但并未表現(xiàn)出明顯的增減趨勢，但在2000年后出現(xiàn)了3次降水高值。累年平均降水為1 021 mm，年均最高降水出現(xiàn)在2000年，年均降水為1 513 mm，年均最低降水出現(xiàn)在1987年，年均降水為563 mm，二者差值高達(dá)950 m。

2.2 江蘇省氣候生產(chǎn)潛力變化分析

2.2.1 江蘇省溫度、降水生產(chǎn)潛力變化分析 根據(jù)公式(1)，計(jì)算1971—2010年江蘇省全省溫度生產(chǎn)潛力(Yt)，結(jié)果見圖2。由圖2可以得出： ①40 a間，全省溫度生產(chǎn)潛力波動性較大(Yt=2.194 6x+2 534.6)，1971—1989年并未表現(xiàn)出規(guī)律性變化，1989年后呈現(xiàn)出在波動中增長趨勢； ②40 a間，溫度生產(chǎn)潛力的平均值為1 833.8 kg/(hm2·a)，最大值出現(xiàn)在2007年，為1 942.8 kg/(hm2·a)，最小值出現(xiàn)在1978年，為1 756.4 kg/(hm2·a)，二者差值為186.5 kg/(hm2·a)； ③溫度生產(chǎn)潛力與溫度表現(xiàn)出一致的變化規(guī)律。

根據(jù)公式(2)，計(jì)算1971—2010年江蘇省全省降水生產(chǎn)潛力(Yr)，結(jié)果見圖2。由圖2可以得出：①全省40 a間降水生產(chǎn)潛力變化規(guī)律不明顯，波動性較大(y=-3.516 6+13 397)； ②40 a間，降水生產(chǎn)潛力平均值為6 397.4 kg/(hm2·a)，最大值出現(xiàn)在1987年，為9 920.8 kg/(hm2·a)，最小值出現(xiàn)在2000年，為4 836.8 kg/(hm2·a)，二者差值高達(dá)5 084.0 kg/(hm2·a)，說明江蘇省年均降水量差異巨大； ③降水生產(chǎn)潛力與降水量呈現(xiàn)出相反的變化規(guī)律。

圖2 1971-2010年江蘇省氣溫生產(chǎn)潛力和降水生產(chǎn)潛力變化

2.2.2 江蘇省氣候生產(chǎn)潛力變化分析 根據(jù)公式(3)，(4)，(5)，計(jì)算1971—2010年江蘇省全省氣候生產(chǎn)潛力，結(jié)果見圖3。由圖3可以看出：①40 a間，江蘇省氣候生產(chǎn)潛力波動性較大，并未表現(xiàn)出明顯的上升趨勢；②40 a間，氣候生產(chǎn)潛力平均值為1 372.2 kg/(hm2·a)，最大值出現(xiàn)在2007年，為1 512.8 kg/(hm2·a)，最小值出現(xiàn)在1987年，為1 094.1 kg/(hm2·a)，二者差值為418.7 kg/(hm2·a)；③從5 a滑動平均曲線來看，1971—1990年氣候生產(chǎn)潛力在波動中降低，1991—2010年又呈現(xiàn)出波動中增加趨勢。

圖3 江蘇省氣候生產(chǎn)潛力和糧食單產(chǎn)變化

對比1977—2010年江蘇省糧食單產(chǎn)發(fā)現(xiàn)，34 a間，江蘇省糧食單產(chǎn)呈現(xiàn)出在波動中逐步上升趨勢，從1977年的299.6 kg/hm2增長為2010年的612.4 kg/hm2，增長量達(dá)312.8 kg/hm2，年均增長率為2.98%；從5 a滑動平均曲線看，1977—1986年糧食單產(chǎn)增長速率較快，年均增長率達(dá)7.17%，1987—1993年糧食單產(chǎn)較為平穩(wěn)，1993—1998年出現(xiàn)新的增長高峰，1993年后糧食單產(chǎn)較為穩(wěn)重有升。

3 江蘇省氣候資源利用率分析

3.1 江蘇省氣候資源利用率變化

根據(jù)氣候資源利用率公式(6)，計(jì)算江蘇省1977—2010年氣候資源利用率見圖4。

注：1977年之前糧食單產(chǎn)數(shù)據(jù)缺失，導(dǎo)致1977之前的氣候資源利用率的數(shù)據(jù)無法獲得。

圖4江蘇省氣候資源利用率變化

由圖4可見，1997—2010年江蘇省氣候資源利用率逐步提高，最高氣候資源利用率出現(xiàn)在1997年，為47%，最低值出現(xiàn)在1977年，為23%；5 a滑動平均曲線同樣表明江蘇省氣候資源利用率在34 a間有了較大提升，特別是在1977—1993年氣候資源利用率大幅提升，表明該階段江蘇省農(nóng)業(yè)生產(chǎn)條件有了較大改善，大大提升了農(nóng)業(yè)生產(chǎn)效率；1993年后氣候資源利用率波動中有小幅上升。

3.2 江蘇省氣候生產(chǎn)潛力影響因素分析

對氣候生產(chǎn)潛力及相關(guān)參數(shù)進(jìn)行相關(guān)性分析結(jié)果見表1。由表1可知，年均溫度與氣候生產(chǎn)潛力并不存在顯著相關(guān)性，而與糧食單產(chǎn)存在顯著相關(guān)性，表明溫度對糧食產(chǎn)量影響較大，而對氣候生產(chǎn)潛力影響不明顯；年均降水與氣候生產(chǎn)潛力存在顯著相關(guān)性，對氣候生產(chǎn)潛力影響較大，而與糧食單產(chǎn)不存在顯著相關(guān)性，表明降水對糧食單產(chǎn)的影響并不明顯，這也與江蘇省農(nóng)田水利設(shè)施高度發(fā)達(dá)、農(nóng)田灌溉可以得到有效保障有關(guān)；氣候生產(chǎn)潛力和糧食單產(chǎn)的相關(guān)性并不顯著，但糧食產(chǎn)量與氣候生產(chǎn)潛力表現(xiàn)了趨同的增長趨勢，特別是氣候資源利用率在研究時段內(nèi)表現(xiàn)出增長趨勢，說明氣候生產(chǎn)潛力對糧食單產(chǎn)產(chǎn)生了一定的影響，但并不是糧食單產(chǎn)的決定因素，進(jìn)一步表明了隨著糧食作物品種的優(yōu)化改良、耕作制度更趨合理、大田管理水平日益科學(xué)等生產(chǎn)力水平的提高，糧食產(chǎn)量逐年提高，且受傳統(tǒng)氣候資源的限制逐漸降低。

表1 江蘇省氣候生產(chǎn)潛力與相關(guān)參數(shù)相關(guān)性分析

為進(jìn)一步說明氣候生產(chǎn)潛力與年均溫度和年均降水的關(guān)系，本文選取1971—2010年年均溫度、降水和氣候生產(chǎn)潛力數(shù)據(jù)，建立三者的數(shù)學(xué)關(guān)系，得到線性回歸結(jié)果見表2。由表2可見，年均溫度和年均降水對氣候生產(chǎn)潛力的影響系數(shù)均達(dá)到極顯著水平，且與氣候生產(chǎn)潛力呈正相關(guān)關(guān)系，但年均降水對氣候生產(chǎn)潛力的影響(系數(shù))小于年均溫度對氣候生產(chǎn)潛力的影響(系數(shù))，方程整體顯著性較好。

表2 江蘇省氣候生產(chǎn)潛力的影響因素分析

4 結(jié) 論

(1) 1971—2010年江蘇省年均氣溫呈現(xiàn)出波動性上升趨勢，特別是1993年以后氣溫的增長趨勢更為明顯；年均降水的年度差異性較大，并未表現(xiàn)出規(guī)律性變化；溫度生產(chǎn)潛力和降水生產(chǎn)潛力的變化規(guī)律與年均氣溫和年均降水相似。

(2) 1971—1990年氣候生產(chǎn)潛力在波動中降低，1991—2010年又呈現(xiàn)出波動中增加趨勢；糧食單產(chǎn)和氣候資源利用率在1977—2010年均表現(xiàn)出波動性上升趨勢，表明二者受氣候的影響相對較小。

(3) 對氣候生產(chǎn)潛力與年均氣溫和年均降水的數(shù)理關(guān)系分析表明，盡管二者均與氣候生產(chǎn)潛力具有較為顯著的相關(guān)性，但年均降水的影響相對微弱，溫度的影響相對更大。

[1] 黃秉維.中國農(nóng)業(yè)生產(chǎn)潛力—光合潛力—農(nóng)業(yè)生態(tài)環(huán)境研究[M].北京:氣象出版社,1989:17-24.

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

[3] 肖國強(qiáng),張強(qiáng),工靜.全球氣候變化對農(nóng)業(yè)生態(tài)系統(tǒng)的影響研究進(jìn)展[J].應(yīng)用生態(tài)學(xué)報,2007,18(8):1877-1885.

[4] 張強(qiáng),王靜.全球氣候變化對農(nóng)業(yè)生態(tài)系統(tǒng)的影響研究進(jìn)展[J]. 應(yīng)用生態(tài)學(xué)報,2007,18(8):1877-1885.

[5] 郭岐峰,傅碩齡.我國農(nóng)業(yè)生產(chǎn)潛力的研究進(jìn)展[J].地理研究,1992,11(4):105-115.

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

[7] 沈思淵,席承藩.淮北渦河流域農(nóng)業(yè)自然生產(chǎn)潛力模型與分析[J].自然資源學(xué)報,1991,6(1):22-33.

[8] 鐘章奇,王錚,夏海濱,等.全球氣候變化下中國農(nóng)業(yè)生產(chǎn)潛力的空間演變[J].自然資源學(xué)報,2015,30(12):2018-2031.

[9] Kaminski J, Kan I, Fleischer A. Structural land-use analysis of agricultural adaptation to climate change: A proactive approach[J]. American Journal of Agricultural Economics, 2013,95(1):70-93.

[10] 陳國南.用邁阿密模型測算我國生物生產(chǎn)力的初步嘗試[J].自然資源學(xué)報,1987,2(3):270-278.

[11] 趙安,趙小敏. FAO-AEZ法計(jì)算氣候生產(chǎn)潛力的模型及應(yīng)用分析[J].江西農(nóng)業(yè)大學(xué)學(xué)報,1998,20(4):528-533.

[12] 林忠輝,莫興國,項(xiàng)月琴.作物生長模型研究綜述[J].作物學(xué)報,2003,29(5):750-758.

[13] Leith H. Modeling the primary productivity of the world[J]. Nature and Resources, 1972,8(2):5-10.

[14] 吳宜進(jìn),熊安元,楊荊安,等.湖北的氣候生產(chǎn)力與農(nóng)業(yè)持續(xù)發(fā)展[J].長江流域資源與環(huán)境,1999,8(4):405-410.