周博揚(yáng)，王君軍，吳和俐

（貴陽(yáng)市氣象局，貴州貴陽(yáng) 550000）

0 引言

作物氣候生產(chǎn)潛力是評(píng)價(jià)一個(gè)地區(qū)在環(huán)境因素適宜情況下，結(jié)合當(dāng)?shù)氐臍夂驐l件，如光、溫、水等在單位面積上能產(chǎn)生的經(jīng)濟(jì)產(chǎn)量最大值[1-4]。通過(guò)數(shù)學(xué)方法對(duì)農(nóng)業(yè)氣候資源來(lái)進(jìn)行分析，建立數(shù)學(xué)模型綜合分析和評(píng)價(jià)當(dāng)?shù)氐臍夂蛸Y源[5-6]。在這方面研究中，李斯（H.Lieth）根據(jù)全球的年平均氣溫、年降水量與作物之間的聯(lián)系，用實(shí)際的蒸散量來(lái)計(jì)算作物氣候生產(chǎn)潛力，建立了Miami模型[7-8]。聯(lián)合國(guó)糧農(nóng)組織推薦了AEZ法，所需資料少，算法簡(jiǎn)單，從而降低了估算的精準(zhǔn)性。R·S·盧米斯[11]引進(jìn)了量子效率研究作物生產(chǎn)潛力。孫忠富等[12]用基本驅(qū)動(dòng)因子的方法，采用太陽(yáng)輻射因子建立了番茄發(fā)育生長(zhǎng)模型。諸多學(xué)者[14-15]考慮氣候變化對(duì)糧食產(chǎn)量的影響，提出建立糧食作物生育期全過(guò)程的模型，能更好地了解本地的農(nóng)業(yè)氣候資源，進(jìn)行開(kāi)發(fā)利用。

位于黔中地區(qū)的貴陽(yáng)市地處云貴高原的東坡，屬亞熱帶季風(fēng)性濕潤(rùn)氣候，地勢(shì)西南高、東北低，是一座具有高原特色的城市。近年來(lái)，全市大力發(fā)展糧食和果蔬產(chǎn)業(yè)，但在氣候變暖的大背景下，貴陽(yáng)市地區(qū)的氣候也隨之變化。該地區(qū)的氣候變化對(duì)作物生產(chǎn)潛力有哪些影響，目前還沒(méi)有研究。針對(duì)貴陽(yáng)市的氣候變化及演變規(guī)律，利用Thornthwaite模型和Miami模型作為基礎(chǔ)進(jìn)行研究，以期更合理地利用農(nóng)業(yè)氣候資源，開(kāi)發(fā)作物生產(chǎn)潛力，為“三農(nóng)”服務(wù)提供更科學(xué)更高效的保障。

1 資料與方法

1.1 資料來(lái)源

本研究選取貴陽(yáng)市有完整的氣象觀測(cè)資料的8個(gè)氣象站1971—2020年（白云為1982—2020年）的氣溫、降水量等基本氣象數(shù)據(jù)。所有數(shù)據(jù)進(jìn)行了嚴(yán)格的質(zhì)量控制，估算農(nóng)業(yè)氣候生產(chǎn)潛力。

1.2 研究方法

選取簡(jiǎn)潔又實(shí)用的Thornthwaite模型。李斯（H.Lieth）和博克斯（E.Box）在第22屆國(guó)際地理學(xué)大會(huì)上提出，他們根據(jù)氣溫和降水與作物產(chǎn)量的關(guān)系，結(jié)合干物質(zhì)積累與光溫水之間密切的聯(lián)系，得出著名的Thornthwaite模型[7-8]，得到了廣泛的應(yīng)用。見(jiàn)公式（1）、（2）、（3）。

式中，Wv為作物氣候生產(chǎn)潛力[g/（m2·a）]；

e為自然對(duì)數(shù)的底；

V為年平均蒸散量（mm）；

R為年降水量（mm）；

L為年平均最大蒸發(fā)量（mm）；

T為年平均氣溫（℃）。

邁阿密模型是李斯用最小二乘法建立的模型，選取年降水量和年平均氣溫為因子[7]?？紤]該模型是降水和溫度兩個(gè)因子，根據(jù)限制因子定律[8]，選擇二者中最小值作為樣本點(diǎn)的作物氣候生產(chǎn)潛力。見(jiàn)公式（4）、（5）。

式中，R為年降水量（mm）；

T為年平均氣溫（℃）；

e為自然對(duì)數(shù)的底；

Wt和Wr分別為以單因子的溫度和降水量計(jì)算的作物氣候生產(chǎn)潛力[g/（m2·a）]。

本文研究采用Thornthwaite模型分析作物氣候生產(chǎn)潛力年際和時(shí)空變化，采用邁阿密模型分析限制因素。

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

采用趨勢(shì)分析、方差分析等統(tǒng)計(jì)方法進(jìn)行統(tǒng)計(jì)分析?？臻g繪圖采用Surfer13.0及克里金插值法，統(tǒng)計(jì)軟件采用SPSS17.0。

2 結(jié)果與分析

2.1 年際變化特征

統(tǒng)計(jì)分析8個(gè)研究對(duì)象的作物氣候生產(chǎn)潛力的逐年數(shù)據(jù)，將其進(jìn)行區(qū)域平均，以作物氣候生產(chǎn)潛力建立時(shí)間序列。分析可知（圖1），貴陽(yáng)市多年平均作物氣候生產(chǎn)潛力為1 399.7 g/（m2·a），2015年最高，作物氣候生產(chǎn)潛力為1 503.9 g/（m2·a），2011年最低，作物氣候生產(chǎn)潛力為1 240 g/（m2·a），最高年與最低年差值為263.9 g/（m2·a）。1971—2020年以來(lái)，貴陽(yáng)地區(qū)的作物氣候生產(chǎn)潛力值呈現(xiàn)出上升趨勢(shì)，趨勢(shì)變化量6.3 g/（m2·a），并通過(guò)了0.05的顯著性檢驗(yàn)，表明均顯著增加。極值出現(xiàn)在21世紀(jì)20年代。作物氣候生產(chǎn)潛力經(jīng)歷了“增-減-增-減-增”的周期性階段

圖1 1971—2020年貴陽(yáng)市作物氣候生產(chǎn)潛力趨勢(shì)變化

2.3 限制因素

2.2 空間變化特征

通過(guò)對(duì)貴陽(yáng)市地區(qū)8個(gè)區(qū)縣的作物氣候生產(chǎn)潛力的歷史值進(jìn)行差值處理分析（圖2a），發(fā)現(xiàn)貴陽(yáng)市的西南部一帶為作物氣候生產(chǎn)潛力高值區(qū)；貴陽(yáng)市東北部，主要在開(kāi)陽(yáng)地區(qū)為低值區(qū)，這與開(kāi)陽(yáng)的氣溫是全市最低關(guān)系密切。由作物氣候生產(chǎn)潛力距平變化分析可知（圖2b），作物氣候生產(chǎn)潛力的增加趨勢(shì)是由東北向西南遞增的。

圖2 1971—2020年貴陽(yáng)市作物氣候生產(chǎn)潛力多年平均值（a）及距平變化（b）空間分布

根據(jù)限制因素定律[8]，計(jì)算Wt和Wr中最小值作為各個(gè)區(qū)縣的作物氣候生產(chǎn)潛力，進(jìn)而找到限制因素。通過(guò)計(jì)算分析，全地區(qū)的Wt多年平均值為1 793.8 g/（m2·a），Wr多年平均值為1 595.8 g/（m2·a）；全市的Wt值均大于Wr值，表明降水因素是作物氣候生產(chǎn)潛力的主要限制因子。分析可知，開(kāi)陽(yáng)和清鎮(zhèn)存在兩個(gè)高值中心，說(shuō)明作物氣候生產(chǎn)潛力對(duì)降水比較敏感；息烽和貴陽(yáng)有兩個(gè)低值中心，說(shuō)明兩地的資源利用率較低，農(nóng)業(yè)灌溉用水量需要更多（圖3）。

圖3 1971—2020年降水量估算的作物氣候生產(chǎn)潛力多年平均值空間分布

2.4 實(shí)際與作物氣候生產(chǎn)潛力關(guān)系

黔中地區(qū)，四季分明，雨熱同季，熱量資源充足，利于農(nóng)業(yè)生產(chǎn)發(fā)展。利用貴陽(yáng)市2011—2020年實(shí)際糧食產(chǎn)量與作物氣候生產(chǎn)潛力的比值，計(jì)算出利用率。分析發(fā)現(xiàn)，10年來(lái)，貴陽(yáng)市糧食總產(chǎn)平均值為419.6 g/（m2·a），作物氣候生產(chǎn)潛力1 425.8 g/（m2·a），生產(chǎn)利用率為30%，說(shuō)明實(shí)際糧食產(chǎn)量還有較大的挖掘潛力（表1）。

表1 貴陽(yáng)市實(shí)際糧食產(chǎn)量與作物氣候生產(chǎn)潛力對(duì)比

3 結(jié)論

1）貴陽(yáng)市多年平均作物氣候生產(chǎn)潛力為1399.7g/（m2·a），最高年與最低年差值為 263.9 g/（m2·a）。1971～2020年以來(lái)，作物氣候生產(chǎn)潛力的曲線周期為“增-減-增-減-增”，目前位于增加趨勢(shì)中。

2）貴陽(yáng)市的西南部一帶，為作物氣候生產(chǎn)潛力高值區(qū)；貴陽(yáng)市東北部地區(qū)為低值區(qū)，這與開(kāi)陽(yáng)的氣溫是全市最低關(guān)系密切。作物氣候生產(chǎn)潛力是由東北向西南增加的空間變化趨勢(shì)。

3）全地區(qū)的Wt多年平均值為1 793.8 g/（m2·a），Wr多年平均值為1 595.8 g/（m2·a）；表明降水因素是作物氣候生產(chǎn)潛力的限制因子。從利用率來(lái)看，貴陽(yáng)市糧食總產(chǎn)平均值為424.4 g/（m2·a），作物氣候生產(chǎn)潛力1 418.7 g/（m2·a），生產(chǎn)利用率為30%，說(shuō)明實(shí)際糧食產(chǎn)量還有較大的挖掘潛力。因此，重視氣候變化產(chǎn)生的影響，提高水分利用率，科學(xué)規(guī)避不利的氣象因素，能有效減緩對(duì)農(nóng)業(yè)生產(chǎn)的影響，進(jìn)一步開(kāi)發(fā)作物氣候生產(chǎn)潛力。

4）在全球氣候變暖的大背景下，黔中地區(qū)的氣候也隨之變化，各種氣象要素的改變必然影響作物生長(zhǎng)發(fā)育，他們之間的聯(lián)系是相互制約的、綜合的、整體的。本研究有齊全的氣象資料作支撐，應(yīng)用了兩種數(shù)學(xué)模型，討論了貴陽(yáng)地區(qū)作物氣候生產(chǎn)潛力對(duì)氣候的響應(yīng)，但這種擬合計(jì)算與實(shí)際生產(chǎn)情況還存在一定的差距。本模型沒(méi)有考慮土壤肥力、光因子、種子優(yōu)化等因素，如何更好地高效地提高作物氣候生產(chǎn)潛力有待進(jìn)一步研究。