莫申萍 蒙小寒 梁虹 劉軍

氣候變化背景下貴港市晚稻氣候適宜度的變化特征

莫申萍1蒙小寒1梁虹1劉軍2

（1. 貴港市氣象局 廣西貴港 537100；2. 貴港市農(nóng)業(yè)農(nóng)村局 廣西貴港 537100）

為氣候變化背景下貴港市水稻生產(chǎn)提供科學(xué)參考，利用貴港市1960—2020年氣象資料和2010—2020年晚稻觀測(cè)資料，基于模糊數(shù)學(xué)計(jì)算晚稻生育期內(nèi)溫度、降水、光照和氣候適宜度，采用線性趨勢(shì)分析方法和M-K突變方法分析各適宜度的變化趨勢(shì)和突變特征，并統(tǒng)計(jì)1960—2020年各年代豐歉年比例。結(jié)果表明：貴港市晚稻溫度適宜度線性趨勢(shì)不確定，突變時(shí)間為2000年，降水、光照和氣候適宜度以下降趨勢(shì)為主，突變時(shí)間多在1990年附近；1960—1989年期間豐年占比多，1990—2020年平年占比多。氣候變化背景下，晚稻生育期內(nèi)氣候適宜度下降、氣象產(chǎn)量偏低是當(dāng)前貴港市水稻生產(chǎn)現(xiàn)狀，有關(guān)部門(mén)應(yīng)積極開(kāi)展水稻品種改善、種植布局優(yōu)化等，以應(yīng)對(duì)氣候變化帶來(lái)的不利影響。

氣候適宜度；貴港；水稻；變化特征

氣候變化已是全球不爭(zhēng)的事實(shí)。氣候變化對(duì)我國(guó)水稻生產(chǎn)產(chǎn)生了深遠(yuǎn)影響，主要表現(xiàn)為各稻區(qū)種植制度和熟性搭配的改變，水稻種植界限北移高擴(kuò)，安全生產(chǎn)期提前或推遲、安全生長(zhǎng)季延長(zhǎng)，以及水稻生育期內(nèi)氣候適宜度的變化，最終導(dǎo)致水稻產(chǎn)量、品質(zhì)發(fā)生改變[1-6]。

貴港市位于我國(guó)廣西東南部、珠江流域干線西江中游，氣候資源豐富，有利于各類(lèi)農(nóng)作物生長(zhǎng)。2021年廣西統(tǒng)計(jì)年鑒[7]顯示，貴港市水稻種植面積2.1×105hm2，占該市總農(nóng)作物種植面積的45%，單產(chǎn)達(dá)5 879 kg/hm2，因此保障水稻的正常生產(chǎn)對(duì)貴港市的糧食安全至關(guān)重要。大量研究表明，在以氣溫上升、降水和日照時(shí)數(shù)分布不均為顯著特征的氣候變化背景下，水稻生育期內(nèi)氣候適宜度已發(fā)生顯著變化，如俞芬等[8]、賴(lài)純佳等[9-10]、宋秋波等[11]對(duì)淮河流域水稻氣候適宜度的研究發(fā)現(xiàn)，水稻生育期內(nèi)光照、降水適宜度呈顯著下降趨勢(shì)；田俊[12]對(duì)江西雙季稻氣候適宜度變化進(jìn)行分析發(fā)現(xiàn)，光照適宜度下降最為明顯，且光照、降水氣候適宜度將長(zhǎng)期維持在較低水平；譚孟祥等[13]對(duì)江蘇一季稻氣候適宜度變化趨勢(shì)進(jìn)行研究得出，溫度、光照、降水和氣候適宜度都呈顯著下降趨勢(shì)。就目前研究來(lái)看，大部分學(xué)者認(rèn)為，氣候變化背景下水稻氣候適宜度降低，當(dāng)前氣候條件對(duì)水稻生產(chǎn)的弊大于利，但不同地區(qū)可能存在差異。在當(dāng)前氣候變化背景下，廣西年平均氣溫顯著上升，且顯著區(qū)域主要分布在高海拔、高緯度的桂北地區(qū)[14]，日照時(shí)數(shù)顯著下降，但貴港地區(qū)日照時(shí)數(shù)變化不明顯[15-16]，且廣西降水仍以非均勻性多尺度的周期變化為主[17]，線性變化趨勢(shì)仍不明確；由于區(qū)域氣候變化的復(fù)雜性，氣候變化對(duì)當(dāng)前貴港市水稻生產(chǎn)可能存在差別性影響。因此，本研究以貴港市晚稻觀測(cè)資料和逐日氣象資料為數(shù)據(jù)源，基于模糊數(shù)學(xué)，采用線性趨勢(shì)方法和Mann-Kendall非參數(shù)統(tǒng)計(jì)檢驗(yàn)法（簡(jiǎn)稱(chēng)M-K方法）分析1960—2020年晚稻氣候適宜度的變化特征，以期在氣候變化背景下掌握水稻氣候適宜度變化規(guī)律，為有關(guān)部門(mén)制定未來(lái)水稻生產(chǎn)計(jì)劃提供科學(xué)參考。

1 數(shù)據(jù)與方法

1.1 數(shù)據(jù)

從貴港市氣象局獲取1960—2020年逐日氣象資料和2010—2020年晚稻觀測(cè)資料，氣象資料主要包括日平均氣溫、降水量、日照時(shí)數(shù)等，晚稻觀測(cè)資料主要包括發(fā)育期、產(chǎn)量。2010—2020年晚稻發(fā)育期的分布箱葉圖如圖1-a所示，以近11年的發(fā)育期平均值作為本研究的各發(fā)育期數(shù)據(jù)；產(chǎn)量情況如圖1-b所示，參考文獻(xiàn)[18]，對(duì)產(chǎn)量進(jìn)行趨勢(shì)產(chǎn)量和氣象產(chǎn)量分離，并求出相對(duì)氣象產(chǎn)量。

圖1 貴港市2010—2020年晚稻發(fā)育期和產(chǎn)量概況

1.2 方法

1.2.1 晚稻氣候適宜度計(jì)算方法[19-20]水稻氣候適宜度由溫度適宜度、降水適宜度和光照適宜度組成，計(jì)算公式如下：

式中，為氣候適宜度，為溫度適宜度，為降水適宜度，為光照適宜度，由beta函數(shù)構(gòu)成，計(jì)算公式如下：

式中，為溫度適宜度，值域[0，1]；T、T和T分別為某一生育階段下限溫度、最適溫度和上限溫度；為該生育階段逐日平均溫度；參數(shù)由下限溫度、最適溫度和上限溫度決定。

是實(shí)際降水量和所需降水量比值的分段函數(shù)，計(jì)算公式如下：

式中，為降水適宜度，值域[0，1]；為生育期內(nèi)的總降水量；R為生育期內(nèi)需水總量。

是底數(shù)為e的指數(shù)函數(shù)，計(jì)算公式如下：

式中，為光照適宜度，值域?yàn)?0，1]；為生育期內(nèi)平均日照時(shí)數(shù)；s為該生育階段最適日照時(shí)數(shù)，當(dāng)實(shí)際日照時(shí)數(shù)超過(guò)該值后，光照適宜度為1，表明光照條件已經(jīng)滿(mǎn)足作物生長(zhǎng)的需要，否則光照適宜度隨日照時(shí)數(shù)的增加呈指數(shù)增加。公式（2）（3）和（4）中各參數(shù)取值參考文獻(xiàn)[21]，如表1所示。

表1 晚稻氣候適宜度模型參數(shù)

1.2.2 趨勢(shì)分析和突變分析 采用線性趨勢(shì)方法和世界氣象組織推薦的M-K法對(duì)溫度適宜度、降水適宜度、光照適宜度和氣候適宜度趨勢(shì)變化和突變情況進(jìn)行分析。

2 結(jié)果與分析

2.1 氣候適宜度模型驗(yàn)證

利用2010—2020年晚稻全生育期氣候適宜度與相對(duì)氣象產(chǎn)量進(jìn)行相關(guān)分析，結(jié)果如圖2所示。得出的回歸方程=67.918–47.019通過(guò)了0.05顯著性水平檢驗(yàn)，說(shuō)明氣候適宜度與相對(duì)氣象產(chǎn)量關(guān)系顯著，建立的回歸模型較好。根據(jù)相對(duì)氣象產(chǎn)量大于5%為豐年、在–5%～5%為平年、小于–5%為歉年的劃分指標(biāo)[20]，利用回歸方程得出豐年與平年、平年與歉年的氣候適宜度閾值分別為0.76、0.62，即氣候適宜度（）大于0.76為豐年，處于0.62至0.76之間為平年，小于0.62為歉年。

圖2 氣候適宜度與相對(duì)氣象產(chǎn)量的相關(guān)性分析

2.2 溫度適宜度趨勢(shì)變化和突變情況

貴港市晚稻各生育期及全生育期溫度適宜度的趨勢(shì)變化和突變情況見(jiàn)表2。從表2可看出，返青期和分蘗期溫度適宜度變化趨勢(shì)顯著，返青期溫度適宜度呈顯著下降趨勢(shì)，分蘗期溫度適宜度呈顯著上升趨勢(shì)；返青期和抽穗開(kāi)花期的溫度適宜度存在突變點(diǎn)，突變時(shí)間均為2000年，突變后返青期溫度適宜度較突變前減少0.05，突變后抽穗開(kāi)花期溫度適宜度較突變前增加0.03。

表2 貴港市晚稻生育期內(nèi)溫度適宜度變化趨勢(shì)及突變情況

注：**為通過(guò)0.01顯著水平檢驗(yàn)，*為通過(guò)0.05顯著水平檢驗(yàn)，“×”表示未出現(xiàn)突變特征，下同。

2.3 降水適宜度趨勢(shì)變化和突變情況

貴港市晚稻各生育期及全生育期降水適宜度的趨勢(shì)變化和突變情況如表3所示。從表3可以看出，乳熟期降水適宜度呈顯著下降趨勢(shì)，苗期、乳熟期和全生育期的降水適宜度存在突變點(diǎn)，突變時(shí)間分別為1989、1989和1990年，突變后苗期、乳熟期和全生育期降水適宜度較突變前均減少，減少量分別為0.08、0.25和0.08。

2.4 光照適宜度趨勢(shì)變化和突變情況

貴港市晚稻各生育期及全生育期光照適宜度的趨勢(shì)變化和突變情況如表4所示。從表4可以看出，分蘗期光照適宜度呈顯著下降趨勢(shì)；苗期、拔節(jié)孕穗期和全生育期的光照適宜度存在突變點(diǎn)，突變時(shí)間分別為2015、1989和2017年，突變后三者光照適宜度較突變前均減少，苗期、拔節(jié)孕穗期和全生育期減少量分別為0.10、0.09和0.12。

表3 貴港市晚稻生育期內(nèi)降水適宜度變化趨勢(shì)及突變情況

表4 貴港市晚稻生育期內(nèi)光照適宜度變化趨勢(shì)及突變情況

2.5 氣候適宜度趨勢(shì)變化和突變情況

貴港市晚稻各生育期及全生育期氣候適宜度的趨勢(shì)變化和突變情況如表5所示。從表5可以看出，乳熟期和全生育期氣候適宜度變化趨勢(shì)顯著，都呈顯著下降趨勢(shì)；苗期、分蘗期、拔節(jié)孕穗期、乳熟期和全生育期的氣候適宜度存在突變點(diǎn)，突變時(shí)間分別為1992、2012、1987、1990和1990年，突變后五者氣候適宜度較突變前都減少，減少量分別為0.06、0.10、0.06、0.06和0.04。

表5 貴港市晚稻生育期內(nèi)氣候適宜度變化趨勢(shì)及突變情況

2.6 1960—2020年豐、平、歉年數(shù)年代變化

根據(jù)2.1建立的回歸方程，分別對(duì)貴港市1960—2020年晚稻氣象條件進(jìn)行評(píng)價(jià)，并以年代為時(shí)間尺度對(duì)豐、平、歉年份數(shù)進(jìn)行統(tǒng)計(jì)，結(jié)果如圖3所示。從圖3可以看出，1960s、1970s和1980s豐年數(shù)量大于或等于平年數(shù)量，不存在歉年；1990s、2000s和2010s平年數(shù)量最多，豐年次之，歉年最少。總體而言，豐年數(shù)量從1990s開(kāi)始明顯下降，平年和歉年比例上升。

圖3 1960—2020年豐、平、歉年份數(shù)年代變化情況

3 討論與結(jié)論

本研究得出在氣候變化背景下，水稻全生育期溫度適宜度線性變化趨勢(shì)不確定，但降水適宜度、光照適宜度和氣候適宜度呈下降趨勢(shì)，這與俞芬等[8]、賴(lài)純佳等[9-10]、宋秋波等[11]、田俊[12]、譚孟祥等[13]關(guān)于水稻氣候適宜度變化趨勢(shì)的研究結(jié)果基本一致。貴港市年總熱量資源顯著增加，但光照和降水的變化趨勢(shì)不明顯，年際變化總體向較好方向發(fā)展，但年際內(nèi)的季節(jié)性變化差異較大，特別是秋冬降水減少而春夏增加，以及夏冬不對(duì)稱(chēng)增溫等[16]，給處于下半年的晚稻生產(chǎn)帶來(lái)諸多風(fēng)險(xiǎn)，因此溫度適宜度變化趨勢(shì)仍不確定。水稻全生育期氣候適宜度以1990年為突變點(diǎn)，比廣西氣溫突變點(diǎn)早10年[14]，比廣西日照突變點(diǎn)晚15年[15]，而與廣西降水突變點(diǎn)相差不大[22]。因此，氣候適宜度的突變時(shí)間點(diǎn)與氣溫和日照的突變時(shí)間點(diǎn)存在一定的提前或滯后效應(yīng)。隨著種植技術(shù)水平的提高和水稻品種的改良，近年來(lái)水稻產(chǎn)量已大幅度提高，但受不同年份氣象條件和農(nóng)業(yè)氣象災(zāi)害影響，氣象產(chǎn)量年際間波動(dòng)性大。本研究得出1960s、1970s和1980s以豐年為主，1990s及后以平年為主，氣象條件對(duì)近年來(lái)水稻產(chǎn)量影響較大，這與黃維等[23]研究一致。

基于前人研究基礎(chǔ)，采用模糊數(shù)學(xué)定量描述晚稻生育期內(nèi)溫度、降水、光照和氣候適宜度，在一定程度上為作物氣象條件鑒定評(píng)價(jià)工作提供參考，但各指標(biāo)和參數(shù)均仍采用前人研究成果，由于水稻不同品種的耐寒性及感光性不同，因此，今后研究中應(yīng)根據(jù)水稻品種劃分不同的溫、光、水指標(biāo)，以進(jìn)一步優(yōu)化水稻氣候適宜度模型。此外，由于不同發(fā)育階段水稻對(duì)光、溫的敏感性不同，全生育階段的適宜度應(yīng)由各生育階段適宜度和權(quán)重系數(shù)的乘積累加而得[24]，但在本研究中采用平均值法進(jìn)行計(jì)算，忽略不同生育階段對(duì)全生育期的貢獻(xiàn)率；因此，在今后研究中應(yīng)采用專(zhuān)家打分法或?qū)哟畏治龇ǖ葘?duì)水稻6個(gè)生育階段的權(quán)重系數(shù)進(jìn)行確認(rèn)，使得全生育期氣候適宜度更能體現(xiàn)整個(gè)晚稻的實(shí)際生長(zhǎng)過(guò)程。

為積極應(yīng)對(duì)氣候變化對(duì)水稻生產(chǎn)的不利影響，氣象部門(mén)應(yīng)進(jìn)一步提高農(nóng)業(yè)氣象中長(zhǎng)期和短期的預(yù)報(bào)準(zhǔn)確率，及時(shí)采取人工影響天氣等防旱消雹措施，減輕氣象災(zāi)害對(duì)水稻及整個(gè)農(nóng)業(yè)生產(chǎn)的影響。農(nóng)業(yè)科研部門(mén)應(yīng)加強(qiáng)科研投入，主攻研發(fā)、選育抗倒伏、抗病力強(qiáng)的水稻品種，從水稻株形、葉面積大小等方面著手提高光能利用率，提高光溫生產(chǎn)潛力，進(jìn)一步提高水稻產(chǎn)量。

氣候變化背景下,晚稻生育期內(nèi)溫度適宜度變化趨勢(shì)不確定，光照、降水和氣候適宜度下降，突變后較突變前適宜度顯著降低，平、歉年比例逐步升高，氣候適宜度總體偏低和氣象產(chǎn)量不高且波動(dòng)性大是貴港市水稻生產(chǎn)現(xiàn)狀，有關(guān)部門(mén)應(yīng)積極開(kāi)展水稻品種改良、種植布局優(yōu)化等研究，以應(yīng)對(duì)氣候變化帶來(lái)的不利影響。

[1] 潘根興, 高民, 胡國(guó)華, 等. 氣候變化對(duì)中國(guó)農(nóng)業(yè)生產(chǎn)的影響[J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào), 2011, 30(9): 1 698-1 706.

[2] 廖雪萍, 黃梅麗, 雍陽(yáng)陽(yáng), 等. 氣候變化對(duì)廣西農(nóng)業(yè)影響的研究進(jìn)展與展望[J]. 氣象研究與應(yīng)用, 2020, 41(4): 72-80.

[3] 黃維, 吳炫柯, 劉永裕, 等. 氣候變化對(duì)廣西雙季稻種植布局的影響[J]. 中國(guó)農(nóng)業(yè)氣象, 2020, 41(9): 539-551.

[4] 黃桂珍, 李莉, 匡昭敏. 氣候變化背景下廣西水稻發(fā)育期變化特征分析[J]. 現(xiàn)代農(nóng)業(yè)科技, 2018(11): 5-13.

[5] 尹朝靜. 氣候變化對(duì)中國(guó)水稻生產(chǎn)的影響研究[D]. 武漢:華中農(nóng)業(yè)大學(xué), 2017.

[6] 凌霄霞, 張作林, 翟景秋, 等. 氣候變化對(duì)中國(guó)水稻生產(chǎn)的影響研究進(jìn)展[J]. 作物學(xué)報(bào), 2019, 45(3): 323-334.

[7] 廣西壯族自治區(qū)統(tǒng)計(jì)局. 廣西統(tǒng)計(jì)年鑒2021[M]. 北京: 中國(guó)統(tǒng)計(jì)出版社, 2021: 321-341.

[8] 俞芬, 千懷遂, 段海來(lái). 淮河流域水稻的氣候適宜度及其變化趨勢(shì)分析[J]. 地理科學(xué), 2008, 28(4): 537-542.

[9] 賴(lài)純佳, 千懷遂, 段海來(lái), 等. 淮河流域雙季稻氣候適宜度及其變化趨勢(shì)[J]. 生態(tài)學(xué)雜志, 2009, 28(11): 2 339- 2 346.

[10] 賴(lài)純佳, 千懷遂, 段海來(lái), 等. 淮河流域小麥-水稻種植制度的氣候適宜性[J]. 中國(guó)農(nóng)業(yè)科學(xué), 2011, 44(14): 2 868-2 875.

[11] 宋秋洪, 千懷遂, 俞芬, 等. 全球氣候變化下淮河流域冬小麥氣候適宜性評(píng)價(jià)[J]. 自然資源學(xué)報(bào), 2009, 24(5): 890-897.

[12] 田俊. 江西省雙季晚稻氣候適宜度變化特征和趨勢(shì)預(yù)測(cè)[J]. 安徽農(nóng)業(yè)科學(xué), 2011, 39(28): 17 485-17 488.

[13] 譚孟祥, 景元書(shū), 曹海寧. 江蘇省一季稻生長(zhǎng)季氣候適宜度及其變化趨勢(shì)分析[J]. 江蘇農(nóng)業(yè)科學(xué), 2016, 44(1): 349-353.

[14] 王瑩, 蘇永秀, 李政. 1961—2010年廣西氣溫變化對(duì)全球變暖的響應(yīng)[J]. 自然資源學(xué)報(bào), 2013, 28(10): 1 707-1 717.

[15] 葉瑜, 蘇永秀, 李政, 等. 近50年廣西日照時(shí)數(shù)時(shí)空變化特征分析[J]. 中國(guó)農(nóng)學(xué)通報(bào), 2013, 29(14): 196-201.

[16] 蒙小寒, 黃維, 莫申萍, 等. 1961—2020年貴港市農(nóng)業(yè)氣候資源變化特征[J]. 熱帶農(nóng)業(yè)科學(xué), 2021, 41(11): 93-99.

[17] 謝鑫昌, 楊云川, 田憶, 等. 廣西降水非均勻性多尺度特征與綜合評(píng)價(jià)[J]. 地球科學(xué)進(jìn)展, 2019, 34(11): 1 152-1 164.

[18] QX/T 81-2007, 小麥干旱災(zāi)害等級(jí)[S].北京: 氣象出版社, 2007.

[19] 黃璜. 中國(guó)紅黃壤地區(qū)作物生產(chǎn)的氣候生態(tài)適應(yīng)性研究[J]. 自然資源學(xué)報(bào), 1996, 11(4): 340-346.

[20] 張建軍, 馬曉群, 許瑩. 安徽省一季稻生長(zhǎng)氣候適宜性評(píng)價(jià)指標(biāo)的建立與試用[J]. 氣象, 2013, 39(1): 88-93.

[21] 黃淑娥, 田俊, 吳慧峻. 江西省雙季水稻生長(zhǎng)季氣候適宜度評(píng)價(jià)分析[J]. 中國(guó)農(nóng)業(yè)氣象, 2012, 33(4): 527-533.

[22] 宋苗, 楊東. 近51年來(lái)廣西不同等級(jí)降水的變化特征[J]. 中國(guó)農(nóng)學(xué)通報(bào), 2014, 30(35): 201-209.

[23] 黃維, 鄭慧琳, 劉永裕, 等. 南充市一季稻高溫?zé)岷Φ陌l(fā)生規(guī)律及其對(duì)產(chǎn)量的影響[J]. 中國(guó)農(nóng)學(xué)通報(bào), 2018, 34(30): 101-109.

[24] 黃維, 楊沈斌, 陳德, 等. 蘇皖鄂地區(qū)一季稻氣候適宜度模型的構(gòu)建[J]. 江蘇農(nóng)業(yè)科學(xué), 2017, 45(2): 232-238.

Variation Characteristics of Climate Suitability of Late Rice under the Background of Climate Change in Guigang

MO Shenping1MENG Xiaohan1LIANG Hong1LIU Jun2

(1. Guigang Meteorological Bureau, Guigang, Guangxi 537100, China; 2. Agricultural and Rural Bureau of Guigang,Guigang, Guangxi, 537100, China)

In order to provide scientific reference for rice production under the background of climate change in Guigang,based on the meteorological data from 1960 to 2020 and the late rice observation data from 2010 to 2020 of Guigang,the temperature, precipitation, light and climate suitability during the growth period of late rice were calculated by fuzzy mathematics. Linear trend analysis method and M-K mutation method were used to analyze the change trend and mutation characteristics of each suitability of late rice in Guigang. The proportion of rich and poor years from 1960 to 2020 was also counted. The results showed that the linear trend of temperature suitability of late rice in Guigang was uncertain. The mutation time of temperature suitability was 2000. Precipitation, light and climate suitability mainly showed a downward trend. The mutation time were mostly near 1990. Rich years from 1960 to 1989 was high proportion while normal years from 1990 to 2020 was high proportion. The current situation of rice production in Guigang is the decline of climate suitability and low meteorological yield during the growth period of late rice under climate change. Relevant departments should actively carry out measures such as rice variety improvement and planting layout optimization to deal with the adverse effects of climate change.

climate suitability; Guigang; rice; change characteristics

S511

A

10.12008/j.issn.1009-2196.2022.06.020

2022-01-27；

2022-03-10

2017年廣西科技重點(diǎn)研發(fā)計(jì)劃項(xiàng)目《基于“3S”技術(shù)的水稻高溫?zé)岷︻A(yù)測(cè)監(jiān)測(cè)和災(zāi)害區(qū)劃技術(shù)研究與示范》（No.桂科AB17195037）；中國(guó)氣象局預(yù)報(bào)員專(zhuān)項(xiàng)（No.CXFZ2021Z034）；貴氣科（No.2022001）。

莫申萍（1980—），女，學(xué)士，農(nóng)藝師，研究方向?yàn)檗r(nóng)業(yè)氣象應(yīng)用，E-mail：mspping@163.com。

蒙小寒（1980—），男，學(xué)士，高級(jí)工程師，主要研究方向?yàn)檗r(nóng)業(yè)氣象，E-mail：42101971@qq.com。

（責(zé)任編輯 林海妹）