李偉光，侯美亭，張京紅，車秀芬，陳小敏

華南主產區(qū)雙季水稻物候變化及其與氣候條件的關系＊

李偉光1,2，侯美亭3，張京紅1,2，車秀芬1,2，陳小敏1,2

（1．海南省氣候中心，海口 570203；2．海南省南海氣象防災減災重點實驗室，?？?570203；3．中國氣象局干部學院，北京 100081）

為探索華南水稻生長發(fā)育對氣候變化的敏感性，估算氣候變化、品種更替對水稻物候期變化的貢獻，選取1981-2013年華南雙季水稻種植區(qū)物候和氣象觀測數據，采用Theil-Sen估計線性傾向率、Mann-Kendall趨勢顯著性檢驗等方法，分析雙季水稻物候期的時空變化趨勢。結果表明：與其它地區(qū)不同，華南早稻播種期提前、晚稻成熟期推遲趨勢較?。辉绲?、晚稻營養(yǎng)生長期存在縮短趨勢，而生殖生長期存在延長趨勢。早稻、晚稻的營養(yǎng)生長期、生殖生長期均隨著降水量和日照時數的增加而延長，隨著氣溫上升而縮短；平均氣溫每上升1℃，生長期長度縮短0.5～4.7d。氣候變化與品種對水稻物候期影響的貢獻率分析表明，品種更替對早稻和晚稻生長期長度的影響大于氣候變化。研究期內品種更替中，早稻傾向于選擇營養(yǎng)生長期和生殖生長期均長的品種；晚稻傾向于選擇營養(yǎng)生長期短、生殖生長期長的品種。

水稻；物候；發(fā)育期；氣候變化；貢獻率

植物物候是一年中植物生長階段的時間記錄，它反應了植物對周圍環(huán)境條件變化的響應，是全球氣候變化的有力證據[1]。天然植被物候的變化已經在全球各地觀測到，主要表現(xiàn)為春季提前、秋季推遲[2]。在農業(yè)生產中，人們會通過改良品種、改進栽培措施來延長生育期，以充分利用熱量資源，提升產量。因此，觀測到的物候變化是氣候變化和品種更新共同作用的結果[3]。國內外已有眾多學者對糧食作物開展物候變化的研究，并取得大量成果。不同地區(qū)的冬小麥對氣候變化的反應基本一致，播種期推遲[4]，小麥營養(yǎng)生長期和生長季都在縮短，而生殖生長期長度呈現(xiàn)延長趨勢[5]。中國春玉米物候期提前而夏玉米和春夏播玉米物候期推遲，營養(yǎng)生長期呈縮短、生殖生長期呈延長趨勢[6]。

水稻作為中國產量最大的糧食作物之一，針對其物候期變化的相關研究也已開展。Tao等發(fā)現(xiàn)中國的單季稻移栽期提前、成熟期推遲，營養(yǎng)生長期、生殖生長期及全生育期均延長；而長江中下游的雙季水稻物候期提前，營養(yǎng)生長期、全生育期縮短，生殖生長期延長[7]。Bai等發(fā)現(xiàn)早稻生育期提前、晚稻生育期延遲，單季稻和早稻的營養(yǎng)生長、生殖生長和總生育期均延長，而晚稻表現(xiàn)為縮短[8]。東北地區(qū)一季稻在2001?2012年間生育期表現(xiàn)出縮短而1992?2001年的生育期表現(xiàn)出延長現(xiàn)象[3]。Hu等研究發(fā)現(xiàn)氣候變暖使單季稻、早稻和晚稻生育期均呈現(xiàn)縮短趨勢[9]。氣象等環(huán)境因子也影響著水稻生育期。Liu等認為氣溫升高、日照時數增加使生育期縮短，而降水增加對延長生育期有幫助[10]。水稻分布廣泛，有秈稻、粳稻兩個亞種，地區(qū)間一季水稻、雙季水稻、早稻晚稻等種植制度不同，加之氣象條件的空間分布差異、季節(jié)增溫的不對稱性[11]共同導致了水稻物候變化的時空不一致。

華南地區(qū)水熱資源充沛，限制及影響水稻物候變化的環(huán)境條件與其它地區(qū)差異較大。因此，本文對華南雙季水稻主產區(qū)早稻、晚稻的物候變化及其與氣象條件的關系展開研究，以期發(fā)現(xiàn)雙季水稻適應氣候變化的規(guī)律，提出品種適應當前氣候變化的針對性措施，進而指導該地區(qū)農業(yè)研究和生產。

1 數據與方法

1.1 農業(yè)氣象站和數據

水稻物候觀測數據來源于中國農作物生長發(fā)育數據集(http://data.cma.cn/data/cdcdetail/dataCode/AGME_ AB2_CHN_TEN.html)，數據時間范圍為1981-2013年。數據的觀測記錄依據《農業(yè)氣象觀測規(guī)范》[12]，水稻物候期包括播種、出苗、三葉、移栽、返青、分蘗、拔節(jié)、孕穗、抽穗、乳熟、成熟共11個。營養(yǎng)生長期(Vegetative growth period，VGP)為播種-抽穗階段，生殖生長期（Reproductive growth period，RGP）為抽穗-成熟階段，生長季（Whole growth period，WGP）為播種-成熟階段。

氣候觀測數據來自中國氣象科學數據共享服務網提供的中國地面氣候資料日值數據集（http://data. cma.cn/data/cdcdetail/dataCode/SURF_CLI_CHN_MUL_DAY.html），選擇日平均氣溫、降水和日照時數進行計算。

數據預處理過程分三步。第一步，選取有連續(xù)物候觀測的農業(yè)氣象觀測站點，對數據進行質量控制，僅選取數據有效率90%以上站點，有23個站點入選，分布于廣東、廣西及海南省的海南島，詳細分布見圖1。第二步，將物候期日期轉換為年日序(Day of year，DOY)。第三步，按站點和年份篩選各物候期的DOY值，分別計算生長季、營養(yǎng)生長期、生殖生長期內的光、溫、水和積溫等氣候要素數值。

圖1 研究區(qū)域位置及站點分布

注：基于國家測繪地理信息局標準地圖服務網站下載的審圖號為GS( 2019) 1698號的標準地圖制作，底圖無修改。

Note:Based on the standard map No. GS (2019) 1698, downloaded from National Bureau of surveying, mapping and geographic information, and the base map was not modified.

1.2 方法

1.2.1 積溫計算

積溫（Growing degree days, GDD）是指某個時間段內日平均氣溫大于基準溫度部分的累積值，具有重要的生物學意義，可以用來估計作物在生長期內的生長發(fā)育情況[13]。計算式為

式中，Dstart為某發(fā)育階段的起始日期，Dend是該發(fā)育階段結束的日期，TEMd是d日內的日平均氣溫，TEMbase為水稻發(fā)育的基準溫度，即發(fā)育溫度下限，設置為10℃[14]。

1.2.2傾向率及顯著性檢驗

早、晚稻物候期和生長期長度的線性傾向率采用Theil-Sen方法計算。Theil-Sen估計較普通最小二乘(OLS)對數據中異常值的容錯性更高[15]。計算式為

式中，n為分析時段的年數；Pi,k和Pj,k分別表示觀測站點在第i年和第j年的k物候期（用年日序表示，Day of year, DOY）或生長期長度（d）；Sk表示k物候期或生長期長度的傾向率（d·a-1），median為對斜率序列取中位數。Sk大于0，表示物候期或生長期長度推遲或延長，反之表示提前或縮短。

根據歷年生長期初始日和終止日觀測值，統(tǒng)計作物生長發(fā)育過程中的平均氣溫、降水量、日照時數和有效積溫。以上氣候要素的傾向率也用式（2）計算。

為檢驗變化趨勢的顯著性，進一步結合Mann-Kendall法[16]進行分析。Mann-Kendall法是一種非參數統(tǒng)計檢驗方法，其優(yōu)點是不需要測量值服從正態(tài)分布，也不要求趨勢是線性的，并且不受缺失值和異常值的影響，廣泛應用于長時間序列數據的趨勢顯著檢驗。計算方法為對于時間序列Xi，i=1,2,3,…，n定義標準化檢驗統(tǒng)計量Z。

1.2.3 水稻生長期長度對各氣候因子變化的敏感度分析

水稻生長期長度對生長期內平均氣溫、降水和日照時數變化的敏感度計算式為[17]

1.2.4 氣候變化及品種對發(fā)育期長度變化的貢獻率計算

觀測到的作物發(fā)育期變化是由作物品種更替、氣候變化以及其它農事活動共同作用的結果。作物發(fā)育期長度的趨勢可以通過品種與氣候因素引發(fā)的發(fā)育期趨勢表示[9，18-19]，即

氣候變化導致的水稻生長期長度變化趨勢[18，20]為

式中，Tphe-cli表示由氣候變化導致的作物發(fā)育期長度變化趨勢。Ttem、Tpre、Tssd分別為氣溫、降水、日照時數的變化趨勢。a、b和c為研究區(qū)域內的所有站點水稻生長期長度對平均氣溫、降水和日照時數敏感度的平均值。

作物品種引發(fā)的水稻生長期長度變化趨勢為

式中，Tphe-GDD表示由作物品種變化導致的作物發(fā)育期長度變化趨勢，TGDD表示水稻積溫隨時間變化的趨勢，TEMavg表示發(fā)育期內的平均氣溫，TEMbase表示水稻發(fā)育的基準溫度（10℃）。

氣候變化和品種更替引發(fā)的發(fā)育期長度趨勢變化的貢獻率計算式[18]分別為

2 結果與分析

2.1 雙季水稻播種期和成熟期變化趨勢

利用1981?2013年華南雙季水稻播種及成熟期數據，經過傾向率估計及顯著性檢驗，雙季水稻播種期和成熟期的變化趨勢見圖2。由圖可見，1981?2013年華南雙季水稻種植區(qū)早稻播種期提前、晚稻成熟期推遲的趨勢較小，而早稻成熟期提前、晚稻播種期推遲的趨勢明顯。早稻播種期平均每10年提前0.3d，成熟期提前1.2d。23個站點中早稻播種期呈現(xiàn)提前趨勢的站點有9個，其中7個通過0.05水平的顯著性檢驗；推遲的站點有10個，其中5個通過0.05水平的顯著性檢驗（圖2a）。早稻成熟期呈現(xiàn)提前趨勢的站點有14個，其中6個站點通過0.05水平的顯著性檢驗；推遲的站點有5個，沒有站點通過顯著性檢驗（圖2b）。晚稻播種、成熟期變化幅度更大，播種期平均每10年推遲4.8d，收獲期推遲0.6d。晚稻播種期表現(xiàn)為推遲的站點有21個，其中15個站點通過0.05水平的顯著性檢驗；僅有1個站點表現(xiàn)為播種期提前（圖2c）。晚稻收獲期有8個站點表現(xiàn)出提前，其中3個站點通過0.05水平的顯著性檢驗；14個站點推遲，其中7個站點通過0.05水平的顯著性檢驗（圖2d）。

圖2 1981?2013年雙季水稻播種期和成熟期線性傾向率的空間分布

注：帶〇的站點表示通過0.05水平的顯著性檢驗。下同。

Note: The site with 〇 passed 0.05 level significance test. The same as below.

早稻成熟、晚稻播種期變化趨勢空間分布上較為一致（圖2b、圖2c），從南到北普遍出現(xiàn)提前、推遲的趨勢，相反表現(xiàn)的站點呈零星分布。早稻播種、晚稻成熟變化趨勢的空間一致性較差（圖2a、圖2d），表現(xiàn)出推遲、提前的站點零星分布，沒有規(guī)律性的聚集。

2.2 雙季水稻生育期長度的變化趨勢

利用雙季水稻營養(yǎng)生長期、生殖生長期及生長季的長度數據，經過線性傾向率估計及顯著性檢驗，得到1981?2013年華南雙季水稻生育期長度的變化趨勢（圖3）。由圖3可見，從早稻、晚稻各生長期長度變化趨勢的空間分布來看，早晚稻生長季普遍存在縮短的趨勢；營養(yǎng)生長期、生殖生長期分開來看，前者呈現(xiàn)出縮短而后者呈現(xiàn)出延長的趨勢。早稻全生育期平均以0.7d·10a-1的速率縮短，23個站點中早稻全生育期縮短的有12個站點，其中5個通過0.05水平的顯著性檢驗；全生育期延長的有8個站點，僅2個通過0.05水平的顯著性檢驗。早稻營養(yǎng)生長期縮短趨勢更大，約以1.1d·10a-1的速率縮短。23個站點中早稻營養(yǎng)生長期縮短的有15個站點，其中5個通過0.05水平的顯著性檢驗；延長的站點僅4個。早稻生殖生長期平均以0.6d·10a-1的速率延長，半數的站點呈現(xiàn)出延長趨勢，其中4個通過0.05水平的顯著性檢驗；9個站點未表現(xiàn)出上升或下降趨勢，3個站點縮短。

圖3 雙季水稻生長期長度線性傾向率的空間分布

Note: VGP is vegetative growth period, RGP is reproductive growth period, WGP is whole growth period. The same as below.

晚稻全生育期、營養(yǎng)生長期長度變化的趨勢更為一致，普遍表現(xiàn)出縮短的趨勢，分別以4.5d·10a-1、4.4d·10a-1的速率縮短。生殖生長期表現(xiàn)出小幅度的延長趨勢，平均以0.3d·10a-1的速率延長。從變化趨勢的空間一致性來看，趨勢也較為統(tǒng)一。全生育期、營養(yǎng)生長期縮短的站點分別有16、19個，分別有12個、11個站點通過0.05水平的顯著性檢驗，但表現(xiàn)為上升的站點均未通過顯著性檢驗。超過2/3的站點表現(xiàn)出生殖生長期延長，但延長幅度較小。

2.3 雙季水稻生長季內主要氣候因子的變化趨勢

利用1981?2013年雙季水稻生長季內的平均氣溫、降水量、日照時數及積溫數據，經過線性傾向率及顯著性檢驗，得到華南雙季水稻生長季內主要氣候因子的變化趨勢。由圖4a可見，華南地區(qū)早稻生長季內降水量變化趨勢不統(tǒng)一，表現(xiàn)出上升和下降趨勢的站點分別為12個和11個。晚稻生長季內大部分站點表現(xiàn)出降水減少的趨勢（圖4b），23個站點中19個站點下降，總體平均以31.6mm·10a-1的速率下降。降水量減少的原因可能是因為晚稻生長季縮短顯著。早稻生長季內多數站點（圖4c）日照時數表現(xiàn)出增加的趨勢，其中3個站點通過0.05水平的顯著性檢驗；總體以16.1h·10a-1的速率增加。而晚稻生長季內有14個站點表現(xiàn)出日照時數下降（圖4d），其中4個站點通過0.05水平的顯著性檢驗，平均以13.7h·10a-1的速率減少。

早稻、晚稻生育期內的平均氣溫上升趨勢均較明顯（圖4e、圖4f），上升速率分別為0.2℃·10a-1和0.1℃·10a-1。早稻生長季內有15個站點平均氣溫上升，上升速率超過0.3℃的站點有11個，9個站點上升趨勢通過顯著性檢驗；呈現(xiàn)下降的站點有8個，呈零星分布。晚稻生長季內平均氣溫上升速率略小，有14個站點平均氣溫上升，其中6個站點上升趨勢通過顯著性檢驗。早稻生長季的有效積溫普遍表現(xiàn)出增加趨勢，所有站點平均以42℃·10a-1的速率增加，12個站點的增加趨勢通過顯著性檢驗，僅有3個站點表現(xiàn)出減少趨勢。與早稻不同，晚稻生長季內多半站點有效積溫表現(xiàn)出輕微下降趨勢（22℃·10a-1），23個站點中有15個站點表現(xiàn)為下降，7個站點下降趨勢通過顯著性檢驗。下降站點在兩廣地區(qū)分布相對集中，而有效積溫增加站點零星分布。有效積溫減少的原因可能也是晚稻生長季縮短引起的。

2.4 雙季水稻生長期長度對氣候因子變化的敏感性

利用雙季水稻營養(yǎng)生長期、生殖生長期長度與相應時段內的降水量、日照時數、平均氣溫數據，建立線性回歸模型，得到1981?2013年華南雙季水稻種植區(qū)生長期長度對相應時段內平均氣溫、降水量和日照時數變化的平均敏感度，結果見表1。由表可見，構建的水稻生長期長度對平均氣溫、降水量和日照時數變化的敏感度方程平均決定系數R2較大，達到0.380～0.762（P＜0.05），說明構建的方程具有較高的可靠性。各生長期長度對平均氣溫的敏感度為負值，說明早稻、晚稻的營養(yǎng)生長和生殖生長期長度均隨氣溫升高而縮短。具體數值顯示，平均氣溫每升高1℃，早稻營養(yǎng)生長期縮短4.7d，生殖生長期縮短2.1d；晚稻營養(yǎng)生長期縮短4.3d，生殖生長期縮短0.5d。與晚稻相比，早稻對平均氣溫更加敏感。早稻和晚稻營養(yǎng)生長期對平均氣溫的敏感度顯著高于生殖生長期。

各生長期長度對日照時數變化的敏感度均為正值，說明各生長期長度均隨日照時數增加而延長。早稻、晚稻不同生長期長度對日照的敏感性相近，日照時數每增加100h，生育階段長度延長4.7～7.8d。

各生長期長度對降水量變化的敏感度數值均為正值，說明早稻、晚稻的營養(yǎng)生長和生殖生長期長度均隨降水增多而延長。營養(yǎng)生長期內降水量每增加100mm生長期延長1.1～1.7d，生殖生長期內降水量每增加100mm生長期延長0.9～2.1d。

2.5 品種更替及氣候變化對生長期長度變化的貢獻率

積溫是植物內在生命活動周期的度量，積溫的變化趨勢能夠在一定程度上反應出品種的變遷。統(tǒng)計1981?2013年水稻各生育期內積溫的平均傾向率表明（表2）,華南雙季水稻主產區(qū)的早稻營養(yǎng)生長期和生殖生長期以及晚稻生殖生長期的積溫均

圖4 雙季水稻生長季內氣象要素線性傾向率的空間分布

表1 雙季水稻生長期長度對各氣候因子變化的平均敏感度及顯著性

表2 雙季水稻積溫變化趨勢折算生長期長度趨勢(CT)

表現(xiàn)出上升趨勢。早稻營養(yǎng)生長期、生殖生長期積溫分別以25.2、21.8℃·d·10a-1的速率增加。而晚稻營養(yǎng)生長期積溫表現(xiàn)為明顯減少的趨勢，速率為43.5d·10a-1；生殖生育期呈現(xiàn)增加趨勢，速率為26.6℃·d·10a-1，增加幅度與早稻生殖生長期相近。根據積溫趨勢與平均氣溫折算生長期趨勢，早稻營養(yǎng)生長期、生殖生長期、晚稻生殖生長期均表現(xiàn)延長趨勢，延長速率分別為1.9、1.2、2.0℃·d·10a-1，晚稻營養(yǎng)生長期表現(xiàn)為以2.5℃·d·10a-1速率縮短。

對1981?2013年23個站點氣候變化及品種更替對水稻生長期長度變化的平均貢獻率分析發(fā)現(xiàn)（圖5），氣候變化對華南雙季水稻生長期長度變化的貢獻較小；品種對生長期長度變化的貢獻較大且一般為正效應，即品種更替延長水稻生長期長度；僅在晚稻營養(yǎng)生長期兩者表現(xiàn)出負效應。在早稻和晚稻的生殖生長期，品種因素均表現(xiàn)出超過50%的延長貢獻率，而氣候變化貢獻率均不足5%。營養(yǎng)生長期，品種更替對早稻與晚稻生長期長度的作用相反，對早稻表現(xiàn)為延長，對晚稻表現(xiàn)為縮短。氣候變化對生長期影響最大的時期為晚稻營養(yǎng)生長期，貢獻率達到45%，表現(xiàn)為縮短生長期長度。

圖5 氣候變化及品種因素對水稻生長期長度變化的貢獻率

3 結論與討論

3.1 結論

（1）1981?2013年華南雙季水稻種植區(qū)早稻播種期、成熟期分別以平均0.3、1.2 d·10a-1的速率提前，晚稻播種期、成熟期分別以4.8、0.6 d·10a-1的速率推遲。早稻成熟期、晚稻播種期變化趨勢的空間分布相對一致，區(qū)域內大部分站點整體提前。早稻和晚稻的營養(yǎng)生長期均存在縮短，生殖生長期均存在延長的趨勢，總生育期也存在縮短的趨勢。

（2）早稻生長季內降水量變化趨勢不明顯，而晚稻生長季內降水量普遍表現(xiàn)出減少的趨勢。早稻生長季內多數地區(qū)日照時數增加而晚稻生長季內大部分區(qū)域日照時數下降。早稻、晚稻生長季內平均氣溫普遍呈上升趨勢。早稻生長季有效積溫表現(xiàn)出增加趨勢，而晚稻生長季有效積溫下降。

（3）早稻、晚稻的營養(yǎng)生長期和生殖生長期長度隨降水量和日照時數增加而延長，隨平均氣溫的上升而縮短，平均氣溫每上升1℃，生長期縮短0.5～4.7d。通過積溫反應的品種變化引起的生長期長度變化來看，早稻營養(yǎng)生長期、生殖生長期及晚稻生殖生長期存在較大的延長趨勢（1.2～2.0d·10a-1），晚稻營養(yǎng)生長期長度趨勢為以2.5d·10a-1速率縮短。貢獻率分析表明品種對生長期長度變化趨勢的影響大于氣候變化。品種更替過程中，早晚稻均傾向于選擇生殖生長期長的品種。但在晚稻營養(yǎng)生長期，品種與氣候變化均表現(xiàn)為負效應，且貢獻了90%以上的作用。

3.2 討論

1981?2013年華南雙季水稻種植區(qū)，早稻播種期提前、晚稻成熟期推遲的趨勢較小，這與該研究區(qū)域以北地區(qū)水稻物候變化趨勢存在明顯不同[7]。其原因可能是華南地區(qū)水熱資源較北方水稻種植區(qū)豐富，水稻界限溫度上的日數已經足夠雙季水稻發(fā)育。水稻種植上利用氣候變化帶來的增溫需求不強烈，所以早稻播種期提前、晚稻成熟期推遲趨勢較小。

華南雙季水稻生育期長度變化趨勢與其它地區(qū)[7-8，10]基本一致，表現(xiàn)為營養(yǎng)生長期和生長季普遍縮短，而生殖生長期延長。從作物產量形成的角度看，生殖生長期越長，接收的光合輻射越多，積累的干物質越多，產量越高[21-22]。所以延長作物生殖生長期是育種家追求的重要目標，這一點在華南早稻和晚稻中均有確實證據。從驅動水稻物候變化因素的貢獻率看，氣候變化對華南雙季水稻的影響遠小于品種更替因素，這可能是因為華南水熱充足，氣候條件的波動仍然在水稻適宜范圍內。晚稻營養(yǎng)生長期，品種更替及氣候變化一致表現(xiàn)為縮短貢獻，產生這個現(xiàn)象的原因可能是晚稻營養(yǎng)生長期是一年中氣溫最高的時段，環(huán)境氣溫可能已經接近水稻生長的最適溫度[23]。在氣溫升高的背景下，農戶主動選擇營養(yǎng)生長期短的品種，作物被動縮短生育期來共同應對高溫脅迫。

農作物的物候變化是氣候變化和農事管理共同作用的結果，區(qū)分二者貢獻對客觀評價氣候變化帶來的影響具有重要意義。通常采用的方法有統(tǒng)計[24]和作物模型[25]，但兩種方法還沒有壓倒性的優(yōu)勢。作物模型的方法可以模擬作物、環(huán)境條件及管理措施之間的相互作用，它需要更多的信息進行模型校準和驗證，而它的跨區(qū)域推廣性較差[26]。統(tǒng)計方法只是建立了作物表現(xiàn)與氣象因素之間的數理關系，而這種關系并不一定是因果關系。本研究早稻、晚稻不同發(fā)育階段對氣象因子的敏感系數大小相近，這也在一定程度上肯定了統(tǒng)計方法結果的可靠性。未來可以結合田間試驗，進一步探索其中的生理基礎，分析水稻對氣候變化的響應機制，增強水稻適應氣候變化的能力。

[1] Ge Q,Wang H,Rutishauser T,et al.Phenological response to climate change in China:a meta-analysis[J].Glob Chang Biology,2015,21(1):265-274.

[2] Chen X,Wang L,Niu Z,et al.The effects of projected climate change and extreme climate on maize and rice in the Yangtze River Basin,China[J].Agricultural and Forest Meteorology,2020:282-283.

[3] Pachauri R K,Allen M R,Barros V R,et al.Climate change 2014:Synthesis report.Contribution of Working Groups I,II and III to the Fifth assessment report of the intergovernmental panel on climate change[M]. Geneva,Switzerland:IPCC, 2014.

[4] 劉玉潔,陳巧敏,葛全勝,等.氣候變化背景下1981-2010中國小麥物候變化時空分異[J].中國科學:地球科學,2018, 48(7):888-898.

Liu Y J,Chen Q M,Ge Q S,et al.Spatiotemporal differentiation of changes in wheat phenology in China under climate change from 1981 to 2010[J].Science China:Earth Sciences,2018,48(7):888-898.(in Chinese)

[5] 張媛鈴,郭曉磊,王娜,等.基于決策系統(tǒng)模擬不同降水年型旱作冬小麥的最佳播期[J].中國農業(yè)氣象,2021,42(6): 475-485.

Zhang Y L,Guo X L,Wang N,et al.Analysis on the optimal sowing date of dry-land winter wheat under different precipitation pattern based on wheat decision system[J]. Chinese Journal of Agrometeorology,2021,42(6):475-485. (in Chinese)

[6] 秦雅,劉玉潔,葛全勝.氣候變化背景下1981-2010年中國玉米物候變化時空分異[J].地理學報,2018,73(5):906-916.

Qin Y, Liu Y J, Ge Q S. Spatiotemporal variations in maize phenology of China under climate change from 1981 to 2010[J].Acta Geographica Sinica,2018,73(5):906-916.(in Chinese)

[7] Tao F L,Zhao Z,Shi W J, et al.Single rice growth period was prolonged by cultivars shifts,but yield was damaged by climate change during 1981-2009 in China, and late rice was just opposite[J].Global change biology,2013,19(10): 3200-3209.

[8] Bai H,Xiao D.Spatiotemporal changes of rice phenology in China during 1981-2010[J].Theoretical and Applied Climatology,2020,140(3):1483-1494.

[9] Hu X,Huang Y,Sun W,et al.Shifts in cultivar and planting date have regulated rice growth duration under climate warming in China since the early 1980s[J].Agricultural and Forest Meteorology,2017,247:34-41.

[10] Liu Y,Zhou W,Ge Q.Spatiotemporal changes of rice phenology in China under climate change from 1981 to 2010[J].Climatic Change,2019,157(2):261-277.

[11] Cai Q,Liu Y,Wang Y,et al.Recent warming evidence inferred from a tree-ring-based winter-half year minimum temperature reconstruction in northwestern Yichang,South Central China,and its relation to the large-scale circulation anomalies[J].International Journal of Biometeorology,2016, 60(12):1885-1896.

[12] 國家氣象局.農業(yè)氣象觀測規(guī)范(水稻)[M].北京:氣象出版社,1993.

National Meteorological Administration.Standard for agrometeorological observation(rice)[M].Beijing:China Meteorology Press,1993.(in Chinese)

[13] Xiao L,Liu L,Asseng S,et al.Estimating spring frost and its impact on yield across winter wheat in China[J]. Agricultural and Forest Meteorology,2018,260-261:154- 164.

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

Huang W,Wu X K,Liu Y Y,et al.Impact of climate change on layout of double cropping rice in Guangxi[J].Chinese Journal of Agrometeorology,2020,41(9):539-551.(in Chinese)

[15] Fernandes R,Leblanc S G.Parametric(modified least squares) and non-parametric(Theil–Sen) linear regressions for predicting biophysical parameters in the presence of measurement errors[J].Remote Sensing of Environment, 2005,95(3):149-153.

[16] 黃家俊,張強,張生,等.基于信息熵的新疆降水時空變異特征研究[J].生態(tài)學報,2017,37(13):4444-4455.

Huang J J,Zhang Q,Zhang S,et al.Information entropy- based analysis of spatial and temporal variation in precipitation in Xinjiang[J].Acta Ecologica Sinica,2017, 37(13):4444-4455.(in Chinese)

[17] Liu Y J, Qin Y,Ge Q S.Spatiotemporal differentiation of changes in maize phenology in China from 1981 to 2010[J].Journal of Geographical Sciences,2019,29(3):351- 362.

[18] He L,Jin N,Yu Q.Impacts of climate change and crop management practices on soybean phenology changes in China[J].Science of the Total Environment,2020,707: 135638.

[19] Verón S R,De A D,Lobell D B.Impacts of precipitation and temperature on crop yields in the Pampas[J].Climatic Change,2015,130(2):235-245.

[20] Liu Y,Chen Q,Ge Q,et al.Modelling the impacts of climate change and crop management on phenological trends of spring and winter wheat in China[J].Agricultural and Forest Meteorology,2018,248:518-526.

[21] Gautam P,Lal B,Nayak A K,et al.Inter-relationship between intercepted radiation and rice yield influenced by transplanting time, method, and variety[J].International Journal of Biometeorology,2019,63(3):337-349.

[22] 吳俊,鄧啟云,袁定陽,等.超級雜交稻研究進展[J].科學通報,2016,61(35):3787-3796.

Wu J,Deng Q Y,Yuan D Y,et al.Progress of super hybrid rice research in China[J].Chinese Science Bulletin,2016, 61(35):3787-3796.(in Chinese)

[23] 宋有金,吳超,李子煜,等.水稻產量對生殖生長階段不同時期高溫的響應差異[J].中國水稻科學,2021,35(2):177-186.

Song Y J,Wu C,Li Z Y,et al.Differential responses of grain yields to high temperature in different stages of reproductive growth in rice[J].Chinese Journal of Rice Science,2021,35(2):177-186.(in Chinese)

[24] Zhao H,Fu Y H,Wang X,et al.Timing of rice maturity in China is affected more by transplanting date than by climate change[J].Agricultural and Forest Meteorology, 2016,216:215-220.

[25] He L,Asseng S,Zhao G,et al.Impacts of recent climate warming,cultivar changes,and crop management on winter wheat phenology across the Loess Plateau of China[J]. Agricultural and Forest Meteorology,2015,200:135-143.

[26] Zhang Q,Zhang W,Li T,et al.Accuracy and uncertainty analysis of staple food crop modelling by the process-based Agro-C model[J].International Journal of Biometeorology, 2021,65:587-599.

Spatiotemporal Changes of Double Cropping Rice Phenology in South China and Relationship with Climatic Conditions

LI Wei-guang1,2, HOU Mei-ting3, ZHANG Jing-hong1,2, CHE Xiu-fen1,2, CHEN Xiao-min1,2

(1.Climate Center of Hainan Province, Haikou 570203, China; 2. Key Laboratory of Meteorological Disaster Prevention and Reduction in South China Sea, Haikou 570203; 3. Cadre College of China Meteorological Administration, Beijing 100081)

In order to explore the sensitivity of rice growth and development to climate change, and to estimate the contribution of climate change and varieties change to the change of rice phenology in South China, the phenological and meteorological observation data of double cropping rice planting in South China from 1981 to 2013 were selected for analysis. Theil-Sen estimation linear tendency rate and Mann-Kendall trend significance test were used to analyze the change trend and spatial distribution characteristics of double cropping rice phenology. Different from other regions, the sowing date of early rice in South China is ahead of schedule and the mature period of late rice is delayed slightly. The vegetative growth period of early and late rice was shortened, while the reproductive growth period was prolonged. The length of vegetative and reproductive growth periods of early rice and late rice increased with the increase of precipitation and sunshine hours, and decreased with the increase of temperature; the length of growth period was shortened by 0.5?4.7 days when the average temperature increased by 1℃. The analysis of the contribution rate of climate change and varieties to the impact of rice phenology showed that the influence of varieties change on the long-term length of early and late rice was greater than that of climate change. In the past 30 years, early rice tended to choose the varieties with long vegetative and reproductive growth period, late rice tends to choose varieties with short vegetative growth period and long reproductive growth period.

Rice；Phenology；Growth phases；Climate change；Contribution

10.3969/j.issn.1000-6362.2021.12.004

李偉光,侯美亭,張京紅,等.華南主產區(qū)雙季水稻物候變化及其與氣候條件的關系[J].中國農業(yè)氣象,2021,42(12):1020-1030

2021?04?27

中國氣象局創(chuàng)新發(fā)展專項（CXFZ2021J070）；氣象軟科學項目（2021ZZXM16）

李偉光，高級工程師，研究方向為氣候變化及氣象災害對農業(yè)的影響,E-mail：163great@163.com