氣候變化影響下我國(guó)冬油菜物候期時(shí)空演變分析

黃穩(wěn)清　黃洪宇　蔣范晨　葉智燕　閆超

摘要：農(nóng)作物物候是農(nóng)業(yè)生產(chǎn)決策的重要依據(jù)，其變化關(guān)系到糧食的可持續(xù)生產(chǎn)，其中，保證油料作物生產(chǎn)是保障我國(guó)糧食生產(chǎn)的重要組成部分。油菜作為我國(guó)第一大油料作物，研究其物候信息隨時(shí)間和空間的變化，可為油菜生產(chǎn)提供策略依據(jù)。本研究利用1992—2013年農(nóng)業(yè)物候觀測(cè)站及氣象站數(shù)據(jù)，通過線性回歸的方法，探索我國(guó)冬油菜物候期的時(shí)空演變及對(duì)氣候因子響應(yīng)程度。結(jié)果表明：1）空間上，緯度每增加1°，冬油菜播種期提前2.84 d，開花期、成熟期分別推遲5.3和3.74 d。2）時(shí)間上，我國(guó)冬油菜的現(xiàn)蕾期、開花期和成熟期平均推遲了3.1、0.03和0.13 d/10 a。3）氣候因子響應(yīng)程度上，溫度、降水和日照3個(gè)氣候因子中，溫度對(duì)冬油菜物候期的影響最大，生長(zhǎng)期內(nèi)溫度每增加1 ℃，現(xiàn)蕾期、開花期和成熟期分別提前6.05、1.94和0.43 d。按照1992—2021年的氣候變化趨勢(shì)，到2050年我國(guó)冬油菜現(xiàn)蕾期、開花期、成熟期將平均提前9.96、2.58和0.88 d。氣候變化將導(dǎo)致作物物候期的變化，進(jìn)而影響到作物的生產(chǎn)。為應(yīng)對(duì)氣候變化對(duì)我國(guó)冬油菜物候期的影響，建議培育和種植生長(zhǎng)期更長(zhǎng)的品種，以保證我國(guó)冬油菜的生產(chǎn)。

關(guān)鍵詞：冬油菜;物候期;物候變化;氣候變化

中圖分類號(hào)：S565.4;S162文獻(xiàn)標(biāo)志碼：A

油菜是世界四大油料作物之一，是我國(guó)食用植物油生產(chǎn)的主要原料來源，兼具油、花、蜜等多維度利用價(jià)值[1]。保證油料作物生產(chǎn)是我國(guó)糧食生產(chǎn)安全的重要組成部分，具有重要的貿(mào)易戰(zhàn)略意義。油菜是我國(guó)第一大油料作物，目前我國(guó)的油菜種植規(guī)模不斷增長(zhǎng)，但在油菜的生產(chǎn)過程中仍然面臨著諸多挑戰(zhàn)[2]。我國(guó)油菜生產(chǎn)中冬油菜種植比例最大，冬油菜物候期變化將影響冬油菜的產(chǎn)量及質(zhì)量。因此，對(duì)冬油菜物候期變化的研究能夠?yàn)槲覈?guó)冬油菜的生產(chǎn)提供決策依據(jù)。

植物物候是植物應(yīng)對(duì)環(huán)境變化時(shí)表現(xiàn)出的一系列生活史對(duì)策，是隨氣候季節(jié)性變化的各生長(zhǎng)階段更替的自然現(xiàn)象[3-4]。過去的幾十年里，學(xué)者們進(jìn)行了大量關(guān)于植被物候變化的研究，發(fā)現(xiàn)大多數(shù)植物表現(xiàn)出了春季物候提前以及秋季物候推遲[5]。針對(duì)不同區(qū)域自然植被的研究表明，氣候變化背景下植被的物候變化普遍表現(xiàn)為生長(zhǎng)開始時(shí)間提前以及生長(zhǎng)結(jié)束時(shí)間推遲[6-10]，反映了自然植被的生長(zhǎng)受到自然氣候變化的影響。氣候變化背景下，我國(guó)學(xué)者進(jìn)行了許多農(nóng)作物物候的變化研究，包括水稻、玉米、冬小麥、大豆、棉花[11-15]等。其中，我國(guó)水稻的主要物候期提前[11]，玉米和大豆的主要物候期推遲[12-13];華北平原棉花吐絮之前的物候期提前，吐絮之后的物候期推遲[14];冬小麥?zhǔn)窃蕉魑?，越冬期之前的物候期推遲，越冬期后的物候期提前[15]。而國(guó)內(nèi)外學(xué)者針對(duì)冬油菜物候期變化的相關(guān)研究，卻得到了不同的結(jié)果。其中，AHMAD等[16]通過線性回歸的方法發(fā)現(xiàn)，1980—2014年間，巴基斯坦日均溫度的上升導(dǎo)致冬油菜播種、出苗、開花和生理成熟期每十年平均推遲6.02、3.14、3.31和1.89 d;REZAEI等[17]通過分段線性回歸的方法發(fā)現(xiàn)，1960—2013年，德國(guó)冬油菜的播種期提前了4 d，成熟期推遲了10 d;張佩等[18]通過多元線性回歸的方法，發(fā)現(xiàn)變暖導(dǎo)致江蘇冬油菜開花期提前;李文靜等[19]通過建立有效積溫與開花時(shí)間的線性關(guān)系，對(duì)長(zhǎng)江流域的冬油菜花期進(jìn)行預(yù)報(bào)，得到實(shí)際開花時(shí)間偏早的結(jié)果。作物的生長(zhǎng)受種植習(xí)慣和播種時(shí)間的影響，其物候變化與自然植被不同，并且農(nóng)藝措施差異也使得作物的物候變化更復(fù)雜，不同作物間物候的變化也存在差異，其中，冬油菜作為越冬作物，物候變化差異更為特殊。

由于我國(guó)幅員遼闊，在過去的二十多年里，我國(guó)冬油菜的物候期狀況在時(shí)間和空間上存在著巨大的差異。過去國(guó)內(nèi)冬油菜物候期的研究大多僅針對(duì)一兩個(gè)物候期或局限在小區(qū)域尺度上[18-20]，沒有反映出我國(guó)不同區(qū)域和氣候條件下冬油菜各個(gè)物候期間的差異。因此，有必要在大區(qū)域尺度上對(duì)冬油菜多個(gè)物候期進(jìn)行研究，探究影響冬油菜物候期變化的因素，為我國(guó)冬油菜生產(chǎn)提供策略依據(jù)。本研究將調(diào)查油菜物候格局的時(shí)空變化，探究油菜物候變化與溫度、日照、降水的關(guān)系，推測(cè)未來油菜物候期的變化并為油菜生產(chǎn)實(shí)踐提供建議。

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

1.1研究區(qū)概況

研究區(qū)包括我國(guó)一級(jí)農(nóng)業(yè)區(qū)中的黃土高原區(qū)、黃淮海區(qū)、西南區(qū)以及長(zhǎng)江中下游區(qū)[21]，地理范圍在98°21′—122°42′E、23°1′—40°43′N之間（圖1）。黃土高原區(qū)按氣候可分為東南部和西北部氣候區(qū)，年平均氣溫分別為12～15 ℃和6～10 ℃，降水量分別為600 mm以上和350～600 mm，地區(qū)海拔較高，晴天太陽輻射強(qiáng)，晝夜溫差大[22]。黃淮海區(qū)年均氣溫10～15 ℃，年降水量為500～950 mm，大部分降水集中在夏季[23]。西南區(qū)年均氣溫在10～22 ℃，熱量豐富，冬暖突出，年降水量600～1 800 mm，降水集中于下半年，該區(qū)光照資源較少，時(shí)空分布差異大[24]。長(zhǎng)江中下游區(qū)是我國(guó)重要的商品糧基地，屬亞熱帶季風(fēng)氣候，年平均氣溫14～18 ℃，年降水量1 000～1 500 mm，水熱條件好[25]。

1.2數(shù)據(jù)來源與處理

1.2.1數(shù)據(jù)來源

油菜物候期資料采用中國(guó)氣象數(shù)據(jù)網(wǎng)（data.cma.cn）提供的農(nóng)業(yè)物候觀測(cè)站點(diǎn)資料。該資料包含1992—2013年的油菜物候數(shù)據(jù)，記錄了油菜多個(gè)物候的發(fā)生日期，其中，有冬油菜記錄的站點(diǎn)103個(gè)。冬油菜物候期驗(yàn)證數(shù)據(jù)通過2021年實(shí)地調(diào)研獲得，經(jīng)過篩選，共獲得現(xiàn)蕾期、開花期和成熟期點(diǎn)位分別為19、17和43個(gè)。日均溫度、日照時(shí)數(shù)、降水等氣象數(shù)據(jù)同樣采用中國(guó)氣象數(shù)據(jù)網(wǎng)提供的中國(guó)地面資料日值數(shù)據(jù)集，選擇相應(yīng)或臨近的氣象站記錄作為物候站點(diǎn)的氣候條件。

1.2.2物候數(shù)據(jù)處理

當(dāng)同一站點(diǎn)同一年份的單個(gè)物候期有多個(gè)物候日期記錄時(shí)，取多個(gè)記錄的中值作為該物候期的時(shí)間?？紤]氣候因子對(duì)冬油菜物候期的影響時(shí)，溫度、降水和日照為播種期開始到相應(yīng)物候期之間的日平均值。22F6B144-B11E-427D-814C-D6087EB6CED3

1.3研究方法

1.3.1顯著性檢驗(yàn)

在物候期時(shí)間序列趨勢(shì)分析中，采用曼-肯德爾（Mann-Kenda11）檢驗(yàn)法[26]，檢驗(yàn)冬油菜物候期變化的顯著性。MK檢驗(yàn)的Z絕對(duì)值大于1.28、1.64、1.96時(shí)，分別表示通過了置信度80%、90%和95%的顯著性檢驗(yàn)。在本文中選擇Z的絕對(duì)值大于1.64（90%）作為顯著性檢驗(yàn)的標(biāo)準(zhǔn)。當(dāng)考慮模型擬合效果以及模型驗(yàn)證的顯著性時(shí)采用t檢驗(yàn)方法[27]，P<0.05時(shí)為顯著，P<0.01時(shí)為極顯著。

1.3.2線性回歸

采用線性回歸的方法，對(duì)冬油菜的物候期及其氣候因子進(jìn)行線性擬合。將緯度作為物候期的影響因素時(shí)，通過一元線性回歸，擬合緯度與物候期的關(guān)系。對(duì)于油菜物候期變化及氣候因子趨勢(shì)同樣采用一元線性回歸，將一元回歸方程斜率作為變化率，正值代表正趨勢(shì)，負(fù)值代表負(fù)趨勢(shì)。探究溫度、降水、日照對(duì)冬油菜物候期的綜合影響時(shí)，則使用多元線性回歸的方法。

2結(jié)果與分析

2.1冬油菜物候期分布

基于冬油菜的播種期、現(xiàn)蕾期、開花期和成熟期，分析油菜物候期的時(shí)空分布變化。根據(jù)農(nóng)業(yè)分區(qū)，冬油菜的物候站點(diǎn)分布在西南區(qū)、長(zhǎng)江中下游區(qū)、黃土高原區(qū)和黃淮海區(qū)，其中大多數(shù)點(diǎn)位分布在西南區(qū)和長(zhǎng)江中下游區(qū)。利用物候站點(diǎn)記錄的1992—2013年的物候期時(shí)間，計(jì)算各站點(diǎn)冬油菜物候期時(shí)間的均值，分布如圖2所示。

根據(jù)物候期分布特點(diǎn)，將西南區(qū)和黃土高原區(qū)歸為西南-黃土區(qū)，長(zhǎng)江中下游區(qū)和黃淮海區(qū)歸為長(zhǎng)江-黃淮區(qū)。從冬油菜物候期時(shí)間來看，播種期總體上呈現(xiàn)從南往北逐漸提前的分布規(guī)律，在西南-黃土區(qū)這一規(guī)律尤為明顯。冬油菜播種期最早在第240天，最遲在第342天?，F(xiàn)蕾日期分布不規(guī)則，沒有明顯的緯度漸變特征。冬油菜現(xiàn)蕾期最早在發(fā)生在第317天，最遲發(fā)生在次年的第132天。冬油菜開花和成熟期總體上呈現(xiàn)從南到北逐漸推遲的趨勢(shì)，尤其是在長(zhǎng)江-黃淮區(qū)。冬油菜開花時(shí)間最早在第9天，最遲在第141天。冬油菜成熟期最早在第89天，最遲在第165天。油菜物候期中，現(xiàn)蕾期的早晚差異最大，達(dá)到了180 d，成熟期差異最小，僅76 d。除了現(xiàn)蕾期外，其他物候期最早和最晚的點(diǎn)位幾乎都分布在最南或最北邊。冬油菜與冬小麥同屬越冬作物，物候期的時(shí)空分布相似，緯度越高播種期越早。冬小麥抽穗期和成熟期緯度越低物候期越早，而冬油菜開花期和成熟期緯度越低物候期越早。冬油菜與冬小麥同樣表現(xiàn)出播種早的站點(diǎn)生長(zhǎng)季較長(zhǎng)，播種晚的站點(diǎn)生長(zhǎng)季較短的特點(diǎn)[15]。

由冬油菜物候期分布（圖2）可知，冬油菜的播種期、開花期、成熟期有較明顯的緯度漸變趨勢(shì)，而現(xiàn)蕾期沒有表現(xiàn)出緯度漸變的特點(diǎn)。其中，在黃土高原區(qū)播種最早，越冬后的物候期到來最晚。緯度與播種期、開花期、成熟期的關(guān)系如表1所示。

由表1可知，冬油菜播種期、開花期和成熟期都體現(xiàn)了明顯的隨緯度漸變特點(diǎn)。研究區(qū)整體上表現(xiàn)為緯度每增加1°，冬油菜播種期提前2.84 d，開花期、成熟期分別推遲5.3和3.74 d;西南-黃土區(qū)緯度每增加1°，冬油菜播種期提前3.67 d，開花期、成熟期分別推遲4.74和2.77 d;長(zhǎng)江-黃淮區(qū)緯度每增加1°，冬油菜播種期提前2.89 d，開花期、成熟期分別推遲6.99和5.04 d。播種期在西南-黃土區(qū)與緯度的相關(guān)性更強(qiáng)，而開花期和成熟期在長(zhǎng)江-黃淮區(qū)與緯度的相關(guān)性更強(qiáng)。研究表明在西南-黃土區(qū)，播種期向南推遲的趨勢(shì)比長(zhǎng)江-黃淮區(qū)明顯，在長(zhǎng)江-黃淮區(qū)開花期和成熟期向北推遲的趨勢(shì)比西南-黃土區(qū)明顯，冬油菜物候期變化存在空間差異。

2.2冬油菜物候歷史趨勢(shì)及變化

探究冬油菜歷史物候期的變化，分析過去我國(guó)冬油菜物候的變化情況。冬油菜播種期受人為決策影響大，在此及下文中僅探究現(xiàn)蕾期、開花期、成熟期的趨勢(shì)及變化。通過一元回歸方程擬合多年物候期，僅對(duì)有十年記錄以上的站點(diǎn)物候進(jìn)行趨勢(shì)計(jì)算，并且通過MK檢驗(yàn)計(jì)算變化的顯著性。有效記錄的站點(diǎn)中，現(xiàn)蕾期、開花期和成熟期的點(diǎn)位分別有39、44和43個(gè)。其中27個(gè)（69.2%）站點(diǎn)現(xiàn)蕾期推遲;開花期推遲和提前的站點(diǎn)數(shù)量相同，但推遲顯著的站點(diǎn)更多;27個(gè)（62.8%）站點(diǎn)成熟期推遲。因此冬油菜現(xiàn)蕾期、開花期和成熟期趨勢(shì)都以推遲為主。1992—2013年，我國(guó)不同區(qū)域冬油菜物候期變化存在差異（圖3）。研究時(shí)間內(nèi)，我國(guó)冬油菜現(xiàn)蕾期、開花期和成熟期平均推遲了3.10、0.03和0.13 d/10 a。表明過去我國(guó)冬油菜物候期的變化存在時(shí)間和空間上的差異，而氣候特征是影響作物物候期的重要因素[28]，因此有必要進(jìn)一步研究氣候因子與冬油菜物候期的關(guān)系及影響程度。

2.3氣候因子與冬油菜物候期的關(guān)系

考慮溫度、日照和降水3個(gè)氣候因子對(duì)冬油菜物候期的影響?；诓シN開始至現(xiàn)蕾期、開花期和成熟期期間的氣候均值，研究現(xiàn)蕾期、開花期和成熟期與3個(gè)氣候因子的關(guān)系。將播種開始至現(xiàn)蕾期、開花期和成熟期的單個(gè)氣候因子與冬油菜物候期時(shí)間進(jìn)行擬合，結(jié)果如表2所示。

結(jié)果顯示，播種期到相應(yīng)物候期的日照、降水與物候期時(shí)間的R²較小，溫度的R²最大。表明3個(gè)氣候因子中，溫度對(duì)物候期的影響最大，但是這種影響在現(xiàn)蕾期往后變?nèi)酰▓D4）。溫度對(duì)現(xiàn)蕾期、開花期和成熟期的影響呈現(xiàn)負(fù)趨勢(shì)，日均溫度越高，物候期到來越早。日照對(duì)現(xiàn)蕾期、開花期和成熟期的影響表現(xiàn)為正趨勢(shì)。降水對(duì)現(xiàn)蕾期和成熟期表現(xiàn)為負(fù)趨勢(shì)，對(duì)開花期則為正趨勢(shì)。冬油菜播種期到現(xiàn)蕾期的發(fā)育經(jīng)歷了晚秋、冬季和早春，并且晚秋和早春的熱量條件也不足。熱量條件不足，導(dǎo)致溫度對(duì)現(xiàn)蕾期影響程度最大。發(fā)育階段到達(dá)開花期和成熟期后，溫度條件逐漸改善，溫度對(duì)物候期的作用降低。相關(guān)研究表明溫度、日照和降水這3個(gè)氣候因子中，溫度是影響油菜物候日期的主要[16，29]因素，與本研究得到的結(jié)果一致。22F6B144-B11E-427D-814C-D6087EB6CED3

2.4氣候因子對(duì)冬油菜物候期的影響程度

冬油菜物候期受溫度、日照和降水的綜合影響。因此，通過播種期到物候期的日平均溫度、日照時(shí)數(shù)和降水量與冬油菜的物候期時(shí)間進(jìn)行多元線性擬合，探究氣候因子對(duì)冬油菜物候期的影響程度，擬合結(jié)果見表3。從表3可知，將溫度、日照和降水作為物候期的影響因素時(shí)，發(fā)育期間內(nèi)溫度每增加1 ℃，冬油菜的現(xiàn)蕾期、開花期和成熟期將分別提前6.05、1.94和0.43 d;日照每增加1 h，現(xiàn)蕾期、開花期和成熟期分別推遲0.86、2.62和2.47 d;降水每增加1 mm，現(xiàn)蕾期和成熟期分別提前0.62和0.82 d，開花期推遲0.59 d。該結(jié)果體現(xiàn)了氣候條件對(duì)冬油菜物候期的影響：增溫使冬油菜物候期提前，日照增多使物候期推遲，降水增加使現(xiàn)蕾期和成熟期提前，開花期推遲。利用2021年實(shí)地調(diào)研的冬油菜物候數(shù)據(jù)與相應(yīng)物候觀測(cè)站模擬的2021年冬油菜物候期相應(yīng)時(shí)間分別進(jìn)行驗(yàn)證。

驗(yàn)證結(jié)果顯示（圖5），相對(duì)于模擬結(jié)果，實(shí)際物候時(shí)間提前和推遲比例分別為43.04%和49.37%，與模擬結(jié)果相一致的比例為6.33%。相較于模擬結(jié)果，實(shí)際現(xiàn)蕾期更傾向于提前，占現(xiàn)蕾期驗(yàn)證數(shù)據(jù)的68.42%，開花期的物候時(shí)間更傾向于推遲（52.94%），成熟期推遲和提前的點(diǎn)位數(shù)量相同。3個(gè)物候期總體驗(yàn)證累計(jì)精度在±7 d、±14 d內(nèi)的比例分別為68.35%、88.61%。模型的驗(yàn)證R²和均方根誤差（RMSE）如表3所示，表明模型模擬的冬油菜現(xiàn)蕾期和開花期的平均誤差在8 d以內(nèi)，成熟期平均誤差在10 d以內(nèi)，模型的模擬效果良好。

2.5冬油菜物候趨勢(shì)及未來變化

過去的幾十年里，由于自然及人為因素導(dǎo)致氣候發(fā)生變化。受地理?xiàng)l件差異影響，不同區(qū)域氣候變化的幅度有所差異，氣候變化也會(huì)影響到冬油菜的物候期。根據(jù)表3中冬油菜物候期與氣候因子的關(guān)系，結(jié)合各物候站點(diǎn)1992—2021年的氣候變化趨勢(shì)，計(jì)算我國(guó)冬油菜物候期的變化率。例如，現(xiàn)蕾期的變化便依據(jù)1992—2021年播種期到現(xiàn)蕾期的氣候變化率計(jì)算，依此類推。

1992—2021年冬油菜播種到成熟期物候站點(diǎn)的氣候變化情況如圖6（a）、（b）、（c）所示。冬油菜生長(zhǎng)期內(nèi)，溫度增加的站點(diǎn)占42.7%（44個(gè)），主要在長(zhǎng)江中下游區(qū)的北部。日照增加的站點(diǎn)占57.3%（59個(gè)），主要在長(zhǎng)江中下游區(qū)中部。絕大部分站點(diǎn)的降水都減少，占93.2%（97個(gè)），降水減少幅度大的區(qū)域在長(zhǎng)江中下游區(qū)中部。物候期變化結(jié)果如圖6（d）、（e）、（f）所示。云貴高原地區(qū)相鄰站點(diǎn)冬油菜開花期的變化率差異較大，而成熟期的差異較小。表明在冬油菜開花期內(nèi)，云貴高原山地垂直的氣候變化差異更大。隨著夏季到來日均溫升高，這樣的氣候變化差異被彌補(bǔ)，因此，云貴高原成熟期變化率差異較小。長(zhǎng)江中下游區(qū)北部以及黃淮海區(qū)南部開花期提前的程度較為明顯，表明該區(qū)域過去增溫明顯，很可能是城市化所致。日照增加和降水減少，使長(zhǎng)江中下游區(qū)中部的冬油菜成熟期發(fā)生推遲。按照1992—2021年的氣候變化趨勢(shì)，冬油菜現(xiàn)蕾期的最大提前速率為-3.25 d/a，最大推遲速率為0.64 d/a;開花期最大提前速率為-0.99 d/a，最大推遲速率為0.2 d/a;成熟期最大提前速率為-0.25 d/a，最大推遲速率為0.15 d/a。按照1992—2021年的氣候變化趨勢(shì)推測(cè)，相較于2021年，2050年我國(guó)冬油菜現(xiàn)蕾期、開花期和成熟期將平均提前9.96、2.58和0.88 d。

3結(jié)論與展望

文中研究了氣候變化對(duì)我國(guó)冬油菜物候期時(shí)空變化的影響，分析了冬油菜物候期對(duì)氣候變化的響應(yīng)程度，并對(duì)未來冬油菜物候期變化進(jìn)行推測(cè)。得到以下結(jié)論：

我國(guó)冬油菜播種期由北向南逐漸推遲，開花期和成熟期由南向北逐漸推遲。溫度的緯度梯度影響了冬油菜的發(fā)育，導(dǎo)致冬油菜播種期、開花期和成熟期呈現(xiàn)南北梯度分布。在播種時(shí)間和氣候的綜合影響下，冬油菜現(xiàn)蕾期沒有呈現(xiàn)出明顯的緯度特征。溫度仍是影響冬油菜物候期的主要因素，溫度對(duì)冬油菜物候期影響最大，氣候變暖使冬油菜的相應(yīng)物候期提前，而日照和降水對(duì)冬油菜物候期的變化趨勢(shì)影響并不明顯。在氣候變化背景下，1992—2013年我國(guó)大多數(shù)地區(qū)冬油菜的現(xiàn)蕾期、開花期、成熟期發(fā)生了推遲。按照1992—2021年氣候趨勢(shì)，到2050年我國(guó)冬油菜現(xiàn)蕾期、開花期、成熟期將平均提前9.96、2.58和0.88 d。冬油菜物候期提前很可能導(dǎo)致產(chǎn)量降低，因此，為保證我國(guó)冬油菜的生產(chǎn)，應(yīng)培育生育期更長(zhǎng)的品種以適應(yīng)氣候的變化。

選取了溫度、降水和日照作為冬油菜物候期變化的影響因素，結(jié)果顯示日照和降水對(duì)冬油菜物候期變化趨勢(shì)的影響并不明顯，但日照和降水作為重要的氣候因子，其對(duì)冬油菜物候期的影響值得進(jìn)一步探討。因此，在不同氣候條件下進(jìn)行地塊試驗(yàn)，量化日照和降水因素，是進(jìn)一步理解日照和降水變化對(duì)冬油菜物候期影響的重要方法。相關(guān)研究表明，氣候是影響作物物候的主要因素，但是農(nóng)作物物候也受農(nóng)藝技術(shù)、管理措施以及作物品種更替等影響。文中僅研究了溫度、日照和降水對(duì)冬油菜物候期的影響，可能忽略了其他影響因素。在以后的研究中，可以納入更多的相關(guān)影響因素，深入揭示冬油菜物候變化的原因及特點(diǎn)。參考文獻(xiàn)：

[1]何微， 李俊， 王曉梅， 等. 全球油菜供需現(xiàn)狀與我國(guó)油菜產(chǎn)業(yè)問題、對(duì)策[J/OL]. 中國(guó)油脂. https： //kns. cnki. net/kcms/detail/61. 1099. ts. 20210929. 1245. 006. html.

[2] 王漢中. 我國(guó)油菜產(chǎn)業(yè)發(fā)展的歷史回顧與展望[J]. 中國(guó)油料作物學(xué)報(bào)， 2010， 32（2）： 300-302.

[3] 莫非， 趙鴻， 王建永， 等. 全球變化下植物物候研究的關(guān)鍵問題[J]. 生態(tài)學(xué)報(bào)， 2011， 31（9）： 2593-2601.22F6B144-B11E-427D-814C-D6087EB6CED3

[4] 趙彥茜， 肖登攀， 柏會(huì)子， 等. 中國(guó)作物物候?qū)夂蜃兓捻憫?yīng)與適應(yīng)研究進(jìn)展[J]. 地理科學(xué)進(jìn)展， 2019， 38（2）： 224-235.

[5] GE Q S， WANG H J， RUTISHAUSER T， et al. Phenological response to climate change in China： a meta-analysis[J]. Global Change Biology， 2015， 21（1）： 265-274.

[6] 吉珍霞， 裴婷婷， 陳英， 等. 黃土高原植被物候變化及其對(duì)季節(jié)性氣候變化的響應(yīng)[J]. 生態(tài)學(xué)報(bào)， 2021， 41（16）： 6600-6612.

[7] 鄧晨暉， 白紅英， 翟丹平， 等. 氣候變化背景下1964—2015年秦嶺植物物候變化[J]. 生態(tài)學(xué)報(bào)， 2017， 37（23）： 7882-7893.

[8] 黃文琳， 張強(qiáng)， 孔冬冬， 等. 1982—2013年內(nèi)蒙古地區(qū)植被物候?qū)Ω珊底兓捻憫?yīng)[J]. 生態(tài)學(xué)報(bào)， 2019， 39（13）： 4953-4965.

[9] 孔冬冬， 張強(qiáng)， 黃文琳， 等. 1982—2013年青藏高原植被物候變化及氣象因素影響[J]. 地理學(xué)報(bào)， 2017， 72（1）： 39-52.

[10]任榮儀， 賀中華， 夏傳花， 等. 貴州省不同地貌類型下植被物候時(shí)空變化及其對(duì)氣候的響應(yīng)[J/OL]. 長(zhǎng)江流域資源與環(huán)境. https： //kns. cnki. net/kcms/detail/42. 1320. X. 20201208. 1143. 002. html.

[11]ZHANG S， TAO F L， ZHANG Z. Rice reproductive growth duration increased despite of negative impacts of climate warming across China during 1981—2009[J]. European Journal of Agronomy， 2014， 54： 70-83.

[12]TAO F L， ZHANG S， ZHANG Z， et al. Maize growing duration was prolonged across China in the past three decades under the combined effects of temperature， agronomic management， and cultivar shift[J]. Global Change Biology， 2015， 20（12）： 3686-3699.

[13]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.

[14]WANG Z B， CHEN J， XING F F， et al. Response of cotton phenology to climate change on the North China Plain from 1981 to 2012[J]. Scientific Reports， 2017， 7（1）： 6628.

[15]王倩. 我國(guó)冬小麥物候變化時(shí)空特征及基于APSIM模型的播種期適宜性評(píng)價(jià)[D]. 南京：南京師范大學(xué)， 2020.

[16]AHMAD S， ABBAS G， FATIMA Z， et al. Quantification of the impacts of climate warming and crop management on canola phenology in Punjab， Pakistan[J]. Journal of Agronomy and Crop Science， 2017， 203（5）： 442-452.

[17]REZAEI E E， SIEBERT S， EWERT F. Climate and management interaction cause diverse crop phenology trends[J]. Agricultural & Forest Meteorology， 2017， 233： 55-70.

[18]張佩， 高蘋， 錢忠海， 等. 油菜花期物候主要限制因子分析及預(yù)報(bào)模型的構(gòu)建[J]. 氣象， 2020， 46（2）： 234-244.

[19]李文靜， 黃蔚薇， 李倩， 等. 長(zhǎng)江流域油菜花期預(yù)報(bào)方法研究[J]. 中國(guó)農(nóng)業(yè)資源與區(qū)劃， 2020， 41（2）： 101-108.

[20]屈柳燕， 王靖， 劉澤華， 等. 青海省共和縣小麥和油菜物候?qū)夂蜃兓捻憫?yīng)[J]. 青海草業(yè)， 2017， 26（3）： 31-36， 46.

[21]全國(guó)農(nóng)業(yè)區(qū)劃委員會(huì)《中國(guó)綜合農(nóng)業(yè)區(qū)劃》編寫組. 中國(guó)綜合農(nóng)業(yè)區(qū)劃[M]. 北京： 農(nóng)業(yè)出版社， 1981.

[22]趙艷霞， 王馥棠， 劉文泉. 黃土高原的氣候生態(tài)環(huán)境、氣候變化與農(nóng)業(yè)氣候生產(chǎn)潛力[J]. 干旱地區(qū)農(nóng)業(yè)研究， 2003（4）： 142-146.

[23]劉志娟， 楊曉光， 王文峰. 氣候變化背景下中國(guó)農(nóng)業(yè)氣候資源變化Ⅳ. 黃淮海平原半濕潤(rùn)暖溫麥-玉兩熟灌溉農(nóng)區(qū)農(nóng)業(yè)氣候資源時(shí)空變化特征[J]. 應(yīng)用生態(tài)學(xué)報(bào)， 2011， 22（4）： 905-912.22F6B144-B11E-427D-814C-D6087EB6CED3

[24]代姝瑋， 楊曉光， 趙孟， 等. 氣候變化背景下中國(guó)農(nóng)業(yè)氣候資源變化Ⅱ. 西南地區(qū)農(nóng)業(yè)氣候資源時(shí)空變化特征[J]. 應(yīng)用生態(tài)學(xué)報(bào)， 2011， 22（2）： 442-452.

[25]尹國(guó)應(yīng)， 張洪艷， 張良培. 2001—2019年長(zhǎng)江中下游農(nóng)業(yè)干旱遙感監(jiān)測(cè)及植被敏感性分析[J/OL]. 武漢大學(xué)學(xué)報(bào)（信息科學(xué)版）. https： //doi. org/10. 13203/j. whugis20210172.

[26]KENDALL M G. Rank correlation methods[J]. British Journal of Psychology， 1990， 25（1）： 86-91.

[27]GAUVREAU K， PAGANO M. Students t test[J]. Nutrition， 1993， 9（4）： 386.

[28]WANG S， WANG E L， WANG F， et al. Phenological development and grain yield of canola as affected by sowing date and climate variation in the Yangtze River Basin of China[J]. Crop and Pasture Science， 2012， 63（5）： 478-488.

[29]廖桂平， 官春云. 不同播期對(duì)不同基因型油菜產(chǎn)量特性的影響[J]. 應(yīng)用生態(tài)學(xué)報(bào)， 2001， 12（6）： 853-858.

（責(zé)任編輯：曾晶）

Spatial-temporal Evolution of Winter Rapeseed Phenology

under Climate Change in China

HUANG Wenqing HUANG Hongyu JIANG Fanchen YE Zhiyan YAN Chao

（1.Key Laboratory of Spatial Data Mining & Information Sharing of Ministry of Education， Fuzhou University， Fuzhou 350116， China;

2.The Academy of Digital China（Fujian）， Fuzhou 350116， China）Abstract： Crop phenology is an important basis for agricultural production decision-making， and its change is related to the sustainable production of food， among which ensuring oil crop production is an important part of guaranteeing food production in China. Rapeseed （Brassica napus L.） is the largest oil crop in China. The changes of its phenological information with time and space were studied to provide a strategic basis for the production of Rapeseed. Based on the data of agricultural phenological stations and meteorological stations from 1992 to 2013， this paper explored the spatial-temporal evolution of winter rapeseed phenology in China and its response to climate factors through linear regression method. The results show as follows： 1） Spatially， the sowing date of winter rapeseed was advanced by 2.84 days， flowering and maturity were delayed by 5.3 days and 3.74 days， respectively， with each increase of latitude of 1°. 2） In terms of time， the budding， flowering and maturity stages of winter rapeseed in China were delayed by 3.1， 0.03 and 0.13 days per decades on average. 3） Among the three climate factors， temperature had the greatest influence on the phenological stage of winter rapeseed， and the budding stage， flowering stage and maturity stage were advanced by 6.05， 1.94 and 0.43 days for each 1 ℃ increase in temperature during the growth period. According to the climate change trend from 1992 to 2021， the bud， flowering and maturity stages of winter rapeseed in China will advance by 9.96， 2.58 and 0.88 days on average by 2050. Climate change will lead to the change of crop phenology， which will affect crop production. In order to cope with the effects of climate change on winter rapeseed phenology in China， it is suggested to cultivate and plant varieties with longer growing period to ensure the production of winter rapeseed in China.

Key words： winter rapeseed; phenology; phenological change; climate change22F6B144-B11E-427D-814C-D6087EB6CED3