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

莫申萍 蒙小寒 梁虹 劉軍

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

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

（1. 貴港市氣象局 廣西貴港 537100；2. 貴港市農業農村局 廣西貴港 537100）

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

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

氣候變化已是全球不爭的事實。氣候變化對我國水稻生產產生了深遠影響，主要表現為各稻區種植制度和熟性搭配的改變，水稻種植界限北移高擴，安全生產期提前或推遲、安全生長季延長，以及水稻生育期內氣候適宜度的變化，最終導致水稻產量、品質發生改變[1-6]。

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

1 數據與方法

1.1 數據

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

圖1 貴港市2010—2020年晚稻發育期和產量概況

1.2 方法

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

式中，為氣候適宜度，為溫度適宜度，為降水適宜度，為光照適宜度，由beta函數構成，計算公式如下：

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

是實際降水量和所需降水量比值的分段函數，計算公式如下：

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

是底數為e的指數函數，計算公式如下：

式中，為光照適宜度，值域為(0，1]；為生育期內平均日照時數；s為該生育階段最適日照時數，當實際日照時數超過該值后，光照適宜度為1，表明光照條件已經滿足作物生長的需要，否則光照適宜度隨日照時數的增加呈指數增加。公式（2）（3）和（4）中各參數取值參考文獻[21]，如表1所示。

表1 晚稻氣候適宜度模型參數

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

2 結果與分析

2.1 氣候適宜度模型驗證

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

圖2 氣候適宜度與相對氣象產量的相關性分析

2.2 溫度適宜度趨勢變化和突變情況

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

表2 貴港市晚稻生育期內溫度適宜度變化趨勢及突變情況

注：**為通過0.01顯著水平檢驗，*為通過0.05顯著水平檢驗，“×”表示未出現突變特征，下同。

2.3 降水適宜度趨勢變化和突變情況

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

2.4 光照適宜度趨勢變化和突變情況

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

表3 貴港市晚稻生育期內降水適宜度變化趨勢及突變情況

表4 貴港市晚稻生育期內光照適宜度變化趨勢及突變情況

2.5 氣候適宜度趨勢變化和突變情況

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

表5 貴港市晚稻生育期內氣候適宜度變化趨勢及突變情況

2.6 1960—2020年豐、平、歉年數年代變化

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

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

3 討論與結論

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

基于前人研究基礎，采用模糊數學定量描述晚稻生育期內溫度、降水、光照和氣候適宜度，在一定程度上為作物氣象條件鑒定評價工作提供參考，但各指標和參數均仍采用前人研究成果，由于水稻不同品種的耐寒性及感光性不同，因此，今后研究中應根據水稻品種劃分不同的溫、光、水指標，以進一步優化水稻氣候適宜度模型。此外，由于不同發育階段水稻對光、溫的敏感性不同，全生育階段的適宜度應由各生育階段適宜度和權重系數的乘積累加而得[24]，但在本研究中采用平均值法進行計算，忽略不同生育階段對全生育期的貢獻率；因此，在今后研究中應采用專家打分法或層次分析法等對水稻6個生育階段的權重系數進行確認，使得全生育期氣候適宜度更能體現整個晚稻的實際生長過程。

為積極應對氣候變化對水稻生產的不利影響，氣象部門應進一步提高農業氣象中長期和短期的預報準確率，及時采取人工影響天氣等防旱消雹措施，減輕氣象災害對水稻及整個農業生產的影響。農業科研部門應加強科研投入，主攻研發、選育抗倒伏、抗病力強的水稻品種，從水稻株形、葉面積大小等方面著手提高光能利用率，提高光溫生產潛力，進一步提高水稻產量。

氣候變化背景下,晚稻生育期內溫度適宜度變化趨勢不確定，光照、降水和氣候適宜度下降，突變后較突變前適宜度顯著降低，平、歉年比例逐步升高，氣候適宜度總體偏低和氣象產量不高且波動性大是貴港市水稻生產現狀，有關部門應積極開展水稻品種改良、種植布局優化等研究，以應對氣候變化帶來的不利影響。

[1] 潘根興, 高民, 胡國華, 等. 氣候變化對中國農業生產的影響[J]. 農業環境科學學報, 2011, 30(9): 1 698-1 706.

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

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

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

[5] 尹朝靜. 氣候變化對中國水稻生產的影響研究[D]. 武漢:華中農業大學, 2017.

[6] 凌霄霞, 張作林, 翟景秋, 等. 氣候變化對中國水稻生產的影響研究進展[J]. 作物學報, 2019, 45(3): 323-334.

[7] 廣西壯族自治區統計局. 廣西統計年鑒2021[M]. 北京: 中國統計出版社, 2021: 321-341.

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

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

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

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

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

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

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

[15] 葉瑜, 蘇永秀, 李政, 等. 近50年廣西日照時數時空變化特征分析[J]. 中國農學通報, 2013, 29(14): 196-201.

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

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

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

[19] 黃璜. 中國紅黃壤地區作物生產的氣候生態適應性研究[J]. 自然資源學報, 1996, 11(4): 340-346.

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

[21] 黃淑娥, 田俊, 吳慧峻. 江西省雙季水稻生長季氣候適宜度評價分析[J]. 中國農業氣象, 2012, 33(4): 527-533.

[22] 宋苗, 楊東. 近51年來廣西不同等級降水的變化特征[J]. 中國農學通報, 2014, 30(35): 201-209.

[23] 黃維, 鄭慧琳, 劉永裕, 等. 南充市一季稻高溫熱害的發生規律及其對產量的影響[J]. 中國農學通報, 2018, 34(30): 101-109.

[24] 黃維, 楊沈斌, 陳德, 等. 蘇皖鄂地區一季稻氣候適宜度模型的構建[J]. 江蘇農業科學, 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年廣西科技重點研發計劃項目《基于“3S”技術的水稻高溫熱害預測監測和災害區劃技術研究與示范》（No.桂科AB17195037）；中國氣象局預報員專項（No.CXFZ2021Z034）；貴氣科（No.2022001）。

莫申萍（1980—），女，學士，農藝師，研究方向為農業氣象應用，E-mail：mspping@163.com。

蒙小寒（1980—），男，學士，高級工程師，主要研究方向為農業氣象，E-mail：42101971@qq.com。

（責任編輯 林海妹）