趙先麗 蔡福 丁抗抗 王賽頔 王笑影 張慧 賈慶宇 溫日紅 吳航

（1．中國(guó)氣象局沈陽(yáng)大氣環(huán)境研究所，遼寧沈陽(yáng) 110166；2.遼寧省氣象信息中心，遼寧沈陽(yáng) 110166；3.遼寧省氣象服務(wù)中心，遼寧沈陽(yáng) 110166；4.錦州市生態(tài)與農(nóng)業(yè)氣象中心，遼寧錦州 121001）

引言

IPCC第六次評(píng)估報(bào)告表明，1850—1900年全球地表平均溫度已上升約1.00℃，2011—2020年平均溫升相比工業(yè)化前（1850—1900年）增高了1.09℃，到21世紀(jì)中期全球平均溫度預(yù)計(jì)將升高達(dá)到或超過(guò)1.50℃，溫度劇烈升高導(dǎo)致世界各地氣候變化加?。?］。氣候變暖背景下，氣候變化對(duì)農(nóng)業(yè)生產(chǎn)的影響更加顯著，特別是作物生長(zhǎng)季內(nèi)氣候條件的變化。東北地區(qū)是中國(guó)增溫最顯著且受氣候變化影響最敏感的地區(qū)之一，作為以雨養(yǎng)農(nóng)業(yè)為主的糧食主產(chǎn)區(qū)，氣候變化對(duì)作物的影響更大［2-5］。1960—2019年?yáng)|北地區(qū)年平均氣溫增加率為0.28℃／10 a，≥0℃積溫增加率為61.16℃／10 a，≥10℃積溫增加率為57.09℃／10 a［6］；1960—2010年?yáng)|北地區(qū)平均氣溫增加率為0.34℃／10 a，其中20世紀(jì)80年代末增溫顯著，日最低氣溫增加率為0.48℃／10 a，日最高氣溫增加率為0.24℃／10 a，≥0℃積溫增加率為60.13℃／10 a，≥10℃積溫增加率為62.79℃／10 a，而年降水量（-6.82 mm／10 a）呈減少的趨勢(shì)且波動(dòng)較大，特別是1981—2010年降水波動(dòng)幅度更大［7-10］。玉米是東北地區(qū)重要的糧食作物，2019年?yáng)|北地區(qū)玉米種植面積已達(dá)1.28×107hm2，產(chǎn)量達(dá)8.87×107t，約占全國(guó)玉米總產(chǎn)量的30.9%，占東北糧食總產(chǎn)量的64.2%［11-12］，可見(jiàn)玉米生產(chǎn)在糧食安全和農(nóng)業(yè)生產(chǎn)中具有舉足輕重的作用。

玉米是喜溫作物，生長(zhǎng)發(fā)育各階段都需要適宜的氣候條件才能完成，生長(zhǎng)季農(nóng)業(yè)氣候資源是玉米生長(zhǎng)和發(fā)育過(guò)程的決定因子。因此，探討氣候變暖背景下東北地區(qū)春玉米水熱資源的變化及其潛在的影響，具有重要的意義。研究表明［13-15］，近年來(lái)東北地區(qū)熱量資源增加，降水資源減少或變化不明顯。關(guān)于東北地區(qū)春玉米物候的相關(guān)研究較多，但由于選取的觀測(cè)站點(diǎn)及農(nóng)業(yè)管理措施、施肥、品種等不同，物候期的變化也不同。研究表明［16-19］，東北地區(qū)玉米41.5%的站點(diǎn)出苗期顯著提前，17.0%的站點(diǎn)顯著推遲；41.5%的站點(diǎn)抽雄期顯著提前，13.2%的站點(diǎn)顯著推遲；43.4%的站點(diǎn)成熟期顯著推遲。對(duì)于生育階段，43.4%的站點(diǎn)營(yíng)養(yǎng)生長(zhǎng)階段長(zhǎng)度縮短，56.6%的站點(diǎn)延長(zhǎng)；22.6%的站點(diǎn)生殖生長(zhǎng)階段長(zhǎng)度縮短，77.4%的站點(diǎn)延長(zhǎng)；33.9%的站點(diǎn)全生育期長(zhǎng)度顯著延長(zhǎng)［16，20-23］。

玉米不同生育期對(duì)水熱因子的敏感性不同，研究春玉米不同生育期水熱因子的變化特征非常重要。因此，利用遼寧省錦州地區(qū)春玉米物侯和氣象觀測(cè)資料，分析氣候變化背景下春玉米生育期的變化特征及其與水熱條件的關(guān)系，探討影響春玉米生育期的水熱因子，以期為應(yīng)對(duì)未來(lái)氣候變化措施提供參考，為合理的開(kāi)發(fā)利用農(nóng)業(yè)氣候資源提供科學(xué)依據(jù)。

1 資料與方法

1.1 研究區(qū)概況

錦州農(nóng)業(yè)氣象試驗(yàn)站（121°10′E，41°09′N）位于遼寧省西南部地區(qū)，海拔高度為25 m，屬于中國(guó)七大玉米帶中的東北玉米帶，典型的溫帶季風(fēng)氣候區(qū)，年平均氣溫為9.50℃，1月平均氣溫為-8.00℃，7月平均氣溫為24.40℃，年降水量為540～640 mm。玉米生育期為5—9月，土壤為典型棕壤。

1.2 資料來(lái)源

所用資料來(lái)源于錦州市生態(tài)與農(nóng)業(yè)氣象中心，包括1980—1984年和2004—2020年春玉米資料和氣象觀測(cè)資料，春玉米觀測(cè)資料為發(fā)育期日期（包括播種、出苗、三葉、七葉、拔節(jié)、抽雄、乳熟和成熟日期），主要為丹玉6號(hào)、丹玉39、丹玉42、丹玉47、丹玉405、農(nóng)大95、中科10號(hào)、錦玉118、海禾69、先玉335、沈海69等中晚熟品種（生育期長(zhǎng)度為120～140 d）。氣象資料包括日平均氣溫、日最高氣溫、日最低氣溫和日降水量。

1.3 研究方法

采用相關(guān)分析法、通徑分析法、趨勢(shì)系數(shù)和線性傾向率法［23-25］分析春玉米發(fā)育期的變化特征及其與水（降水量、有效降水、累積降水和累積有效降水）熱（日最高氣溫、日最低氣溫、日平均氣溫、≥10℃活動(dòng)積溫和生長(zhǎng)度日）條件的關(guān)系；其中，苗期是指播種日期—出苗普遍期，三葉期是指出苗普遍期—三葉普遍期，七葉期是指三葉普遍期—七葉普遍期，拔節(jié)期是指七葉普遍期—拔節(jié)普遍期，抽雄期是指拔節(jié)普遍期—抽雄普遍期，乳熟期是指抽雄普遍期—乳熟普遍期，成熟期是指乳熟普遍期—成熟普遍期，營(yíng)養(yǎng)生長(zhǎng)期是指播種日期—抽雄普遍期，生殖生長(zhǎng)期是指抽雄普遍期—成熟普遍期，全生育期是指播種日期—成熟普遍期。

1.3.1 通徑分析

通徑分析是分析多個(gè)自變量和因變量之間線性關(guān)系的一種重要統(tǒng)計(jì)方法。根據(jù)通徑分析，可以獲得每個(gè)自變量對(duì)因變量影響的直接作用和間接作用，具有直觀、精確的優(yōu)點(diǎn)，可以全面反映每個(gè)自變量對(duì)因變量的作用效應(yīng)，揭示每個(gè)變量對(duì)研究結(jié)果的相對(duì)重要性。

對(duì)于多個(gè)自變量和因變量相互關(guān)聯(lián)的系統(tǒng)，n個(gè)自變量x（i＝1，2，…，n）和1個(gè)因變量y之間的回歸方程為

直接通徑系數(shù)（Pi）表示自變量xi對(duì)因變量y的直接作用，公式為

式（2）中，bxi為自變量的回歸系數(shù)；σxi為自變量的標(biāo)準(zhǔn)差；σy為因變量標(biāo)準(zhǔn)差。

間接通徑系數(shù)Pij，如xj對(duì)y的間接通徑系數(shù)公式為

式（3）中，Pij為自變量xi通過(guò)自變量xj對(duì)因變量y的間接作用；rij為兩個(gè)性狀間的相關(guān)系數(shù)；Pjy為xj對(duì)y的直接通徑系數(shù)。相關(guān)系數(shù)的平方即為決定系數(shù)，在多元回歸分析中，通徑系數(shù)的平方即為決定系數(shù)，變量間相關(guān)的密切程度通過(guò)決定系數(shù)大小決定，某個(gè)自變量x對(duì)因變量y的決定系數(shù)為

式（4）中，Pi為自變量xi對(duì)因變量y的通徑系數(shù)。

兩個(gè)不同自變量對(duì)因變量的決定系數(shù)dyij為

式（5）中，Pyi、Pyj分別表示自變量xi和xj與因變量y的直接通徑系數(shù)；rij為自變量xi與xj的相關(guān)系數(shù)［25］。

1.3.2 ≥10℃活動(dòng)積溫（DT10）

活動(dòng)溫度指高于或等于生物學(xué)下限溫度的日平均溫度，≥10℃活動(dòng)積溫（DT10，℃·d）是玉米生長(zhǎng)的重要熱量指標(biāo)之一，用玉米生長(zhǎng)季內(nèi)日平均氣溫≥10℃持續(xù)期間日平均氣溫之和表示［26］，即

式（6）中，a、b分別為玉米生長(zhǎng)季內(nèi)逐日平均氣溫≥10℃的起始日期和終止日期，采用5日滑動(dòng)平均法確定；Ti為日平均氣溫≥10℃持續(xù)期間日平均氣溫（℃）。

1.3.3 生長(zhǎng)度日（GDD）

生長(zhǎng)度日（Growing Degree Days，GDD）是分析作物生長(zhǎng)的整個(gè)熱量指標(biāo)，主要是指在某一特定環(huán)境因子影響下，植物能夠完成整個(gè)生育階段所積累的有效積溫，能夠較好地表征作物生育期的熱量資源。GDD公式為

式（7）中，Tmax為日最高氣溫；Tmin為日最低氣溫；Tb為玉米生長(zhǎng)季生物學(xué)下限溫度10℃，玉米生長(zhǎng)的溫度范圍為10～30℃。本文采用兩種計(jì)算方法，第一種（GDD1），（Tmax＋Tmin）／2＝Tave，Tave為日平均氣溫，當(dāng)Tave＜10℃時(shí)，記為10℃，當(dāng)Tave＞30℃時(shí)，記為30℃；第二種（GDD2），當(dāng)Tmin和Tmax＜10℃時(shí)，記為10℃，當(dāng)Tmin和Tmax＞30℃時(shí)，記為30℃［13，27-28］。

1.3.4 有效降水（Pe）

有效降水（Effective Precipitation，Pe）指干旱地區(qū)滿足作物蒸發(fā)消耗的降水量，能夠滲到地下，被植物根系所吸收的降水。本文采用有效降水［29-30］經(jīng)驗(yàn)公式為

式（8）中，Pe為有效降水量（mm）；a為有效降水系數(shù)；P為日降水量（mm）。通常認(rèn)為5 mm以下降水為無(wú)效降水，因此a的取值采用以下經(jīng)驗(yàn)系數(shù)

1.3.5 累積降水和累積有效降水

累積降水指某一生育階段日降水量的累積求和；累積有效降水指某一生育階段日有效降水的累積求和。

2 結(jié)果分析

2.1 春玉米生育期變化趨勢(shì)

2.1.1 播種日期和成熟日期變化

錦州地區(qū)春玉米播種日期變化較大（圖1a），距平變化為-13.00～19.00 d，整體呈推后的趨勢(shì)，傾向率為0.92 d·a-1；春玉米成熟日期整體呈推后的趨勢(shì)，傾向率為0.90 d·a-1，其中1980—1983年均為提前趨勢(shì)，距平變化為-23.00～-13.00 d，2004年后基本為正距平（圖1b），距平變化為-4.00～13.00 d，2004年以后成熟日期呈微弱的提前趨勢(shì)，傾向率為-0.10 d·a-1。

圖1 1980—1984年和2004—2020年錦州地區(qū)春玉米播種日期（a）和成熟日期（b）距平變化Fig.1 Variation of an anomaly of sowing date（a）and maturity date（b）of spring maize in Jinzhou from 1980 to 1984 and 2004 to 2020

2.1.2 生育期長(zhǎng)度變化

錦州地區(qū)春玉米出苗期長(zhǎng)度總體呈顯著（P＜0.05，下同）縮短的趨勢(shì)（表1），距平變化為-4.82～4.18 d，傾向率為-0.28 d·a-1（圖2a）；三葉期長(zhǎng)度總體呈縮短的趨勢(shì)，距平變化為-3.14～3.86 d，傾向率為-0.18 d·a-1（圖2b）；七葉期長(zhǎng)度總體呈縮短的趨勢(shì)，距平變化為-7.68～7.32 d，傾向率為-0.35 d·a-1（圖2c）；拔節(jié)期長(zhǎng)度總體呈極顯著（P＜0.01，下同）縮短的趨勢(shì)，距平變化為-5.95～9.05，傾向率為-0.26 d·a-1（圖2d）；抽雄期長(zhǎng)度總體呈顯著延長(zhǎng)的趨勢(shì)，距平變化為-5.59～7.41 d，傾向率為0.33 d·a-1（圖2e）；乳熟期長(zhǎng)度總體呈弱的縮短趨勢(shì)，距平變化為-8.00～8.00 d，傾向率為-0.08 d·a-1（圖2f）；成熟期長(zhǎng)度總體呈極顯著延長(zhǎng)的趨勢(shì)，距平變化為-20.27～12.73 d，傾向率為0.80 d·a-1（圖2g）。

圖2 1980—1984年和2004—2020年錦州地區(qū)春玉米出苗期（a）、三葉期（b）、七葉期（c）、拔節(jié)期（d）、抽雄期（e）、乳熟期（f）和成熟期（g）長(zhǎng)度距平變化Fig.2 Variation of length anomaly of seedling（a），three-leaf（b），seven-leaf（C），jointing（d），tasseling（e），milk maturity（F），and maturity（g）of spring maize in Jinzhou from 1980-1984 and 2004-2020

表1 1980—1984年和2004—2020年錦州地區(qū)春玉米生育期長(zhǎng)度變化趨勢(shì)Table 1 Changing trend of growth period length of spring maize in Jinzhou from 1980 to 1984 and 2004 to 2020

錦州地區(qū)春玉米營(yíng)養(yǎng)生長(zhǎng)期長(zhǎng)度總體呈極顯著縮短的趨勢(shì)（表1），距平變化為-13.18～10.82 d，傾向率為-0.74 d·a-1（圖3a）；生殖生長(zhǎng)期長(zhǎng)度總體呈極顯著延長(zhǎng)的趨勢(shì)，距平變化為-17.27～11.73 d，傾向率為0.72 d·a-1（圖3b）；全生育期長(zhǎng)度總體呈弱的延長(zhǎng)趨勢(shì)，距平變化為-14.45～20.55 d，傾向率為-0.02 d·a-1（圖3c）。

圖3 1980—1984年和2004—2020年錦州地區(qū)春玉米營(yíng)養(yǎng)生長(zhǎng)期（a）、生殖生長(zhǎng)期（b）和全生育期（c）長(zhǎng)度距平變化Fig.3 Variation of length anomaly of the vegetative period（a），reproductive period（b），and whole growth period（c）of spring maize in Jinzhou from 1980 to 1984 and from 2004 to 2020

2.2 水熱條件變化趨勢(shì)

1980—1984年和2004—2020年錦州地區(qū)春玉米各生育期水熱條件的變化趨勢(shì)不同，其中熱量條件變化顯著，而水分條件變化不顯著。對(duì)于熱量條件，春玉米多數(shù)生育期氣溫、GDD和DT10呈增加趨勢(shì)，其中出苗期的日平氣溫、日最高氣溫和日最低氣溫均極顯著增加，抽雄期、乳熟期、營(yíng)養(yǎng)生長(zhǎng)期、全生育期的日平均氣溫和日最低氣溫均極顯著增加；拔節(jié)期的GDD2和DT10均極顯著減少；抽雄期、成熟期及生殖生長(zhǎng)期、全生育期的GDD1、GDD2和DT10均極顯著增加。對(duì)于水分條件，成熟期的累積降水量和累積有效降水顯著增加（表2）。

表2 1980—1984年和2004—2020年錦州地區(qū)春玉米生育期水熱條件變化趨勢(shì)Table 2 Changing trend of hydrothermal conditions during the growth period of spring maize in Jinzhou from 1980 to 1984 and from 2004 to 2020

2.3 水熱條件與生育期長(zhǎng)度關(guān)系

2.3.1 相關(guān)分析

玉米生育期的影響因子中，品種影響約占1／4，而氣象條件的影響則占3／4，可見(jiàn)玉米生育期長(zhǎng)度的關(guān)鍵影響因子是氣候條件［18］。錦州地區(qū)春玉米各生育期長(zhǎng)度基本與DT10和GDD呈極顯著正相關(guān)，相關(guān)系數(shù)在0.60以上；出苗期和營(yíng)養(yǎng)生長(zhǎng)期長(zhǎng)度與日平均氣溫、日最高氣溫、日最低氣溫均呈極顯著負(fù)相關(guān)；成熟期和全生育期長(zhǎng)度與累積降水和累積有效降水呈極顯著正相關(guān)，其次為出苗期，再次為三葉期、七葉期和營(yíng)養(yǎng)生長(zhǎng)期、生殖生長(zhǎng)期（表3）。

表3 1980—1984年和2004—2020年錦州地區(qū)春玉米生育期水熱條件與生育期長(zhǎng)度的相關(guān)分析Table 3 Correlation analysis between hydrothermal conditions and growth period length of spring maize in Jinzhou from 1980 to 1984 and 2004 to 2020

2.3.2 通徑分析

為進(jìn)一步確定水熱條件對(duì)春玉米生育期長(zhǎng)度的影響效應(yīng)，對(duì)其進(jìn)行通徑分析。DT10對(duì)春玉米各生育期長(zhǎng)度的直接作用最大，直接通徑系數(shù)為0.80以上；其次為GDD，直接通徑系數(shù)為0.50以上；而水分條件對(duì)春玉米生育期長(zhǎng)度的影響較小，主要對(duì)全生育期長(zhǎng)度影響稍大（表4）。

表4 1980—1984年和2004—2020年錦州地區(qū)春玉米生育期水熱條件與生育期長(zhǎng)度的通徑系數(shù)Table 4 Path analysis between hydrothermal conditions and growth period length of spring maize in Jinzhou from 1980 to 1984 and 2004 to 2020

間接作用中（表略），對(duì)于出苗期、三葉期、七葉期、拔節(jié)期、抽雄期、乳熟期、成熟期及營(yíng)養(yǎng)生長(zhǎng)期、生殖生長(zhǎng)期，GDD和DT10通過(guò)互相影響進(jìn)而影響各生育期長(zhǎng)度；對(duì)于全生育期，日平均降水和累積降水通過(guò)互相影響對(duì)其長(zhǎng)度貢獻(xiàn)較大，其次為日平均降水和累積降水通過(guò)對(duì)DT10影響對(duì)其長(zhǎng)度貢獻(xiàn)其次。

3 結(jié)論與討論

（1）1980—1984年和2004—2020年錦州地區(qū)春玉米播種日期呈延后的趨勢(shì)，受人為因素和氣候因素影響，波動(dòng)變化較大，成熟日期2004年以后波動(dòng)變化相對(duì)較小。播種日期的主觀決定因素較多，播種時(shí)間不僅受當(dāng)時(shí)的溫度降水和土壤墑情等環(huán)境因素影響，還受到耕作方式、主觀意愿和個(gè)人經(jīng)驗(yàn)等因素的影響，而成熟日期受播種日期和水熱條件影響相對(duì)較大。春玉米各生育期長(zhǎng)度變化趨勢(shì)不同，除抽雄期、成熟期和生殖生長(zhǎng)期，其他生育期長(zhǎng)度均呈縮短的趨勢(shì)，其中出苗期長(zhǎng)度顯著縮短，拔節(jié)期和營(yíng)養(yǎng)生長(zhǎng)期長(zhǎng)度極顯著縮短，抽雄期長(zhǎng)度顯著延長(zhǎng)，成熟期和生殖生長(zhǎng)期長(zhǎng)度極顯著延長(zhǎng)，與李正國(guó)等［17］、秦雅等［18］和崔耀平等［31］研究結(jié)論一致。綜上可見(jiàn)，由于營(yíng)養(yǎng)生長(zhǎng)期長(zhǎng)度呈極顯著縮短的趨勢(shì)，而生殖生長(zhǎng)期長(zhǎng)度呈極顯著延長(zhǎng)的趨勢(shì)，造成整個(gè)生育期長(zhǎng)度變化不明顯。

（2）春玉米各生育期水熱條件變化趨勢(shì)不同，熱量條件變化顯著，而水分條件變化不顯著；GDD和DT10多數(shù)生育期呈增加的趨勢(shì)，與部分學(xué)者［13，18，23，32］的研究結(jié)果一致，其中抽雄期、成熟期和生殖生長(zhǎng)期、全生育期的GDD和DT10極顯著增加；出苗期的日平氣溫、日最高氣溫和日最低氣溫均極顯著增加，抽雄期、乳熟期、營(yíng)養(yǎng)生長(zhǎng)期和全生育期的日平均氣溫和日最低氣溫極顯著增加；拔節(jié)期的GDD和DT10極顯著減少；抽雄期、成熟期及生殖生長(zhǎng)期、全生育期的GDD和DT10極顯著增加，成熟期的累積降水量和累積有效降水顯著增加。綜上可見(jiàn)，錦州地區(qū)熱量條件對(duì)春玉米的生長(zhǎng)發(fā)育速度影響較大，進(jìn)而影響發(fā)育期長(zhǎng)度。

（3）通過(guò)對(duì)錦州地區(qū)春玉米水熱條件與生育期長(zhǎng)度的相關(guān)分析和通徑分析表明，對(duì)春玉米生育期長(zhǎng)度影響較大的因素為DT10和GDD，水分條件影響較小，其中出苗期、成熟期及全生育期長(zhǎng)度受水分條件影響稍大，水熱因子之間相互制約、相互影響，共同影響作物的整個(gè)生育期。基于日平均氣溫計(jì)算的生長(zhǎng)度日GDD1對(duì)錦州地區(qū)春玉米生育期長(zhǎng)度負(fù)向效應(yīng)較大，表明隨著日平均氣溫升高，生育期長(zhǎng)度將縮短，可見(jiàn)基于兩種算法計(jì)算的GDD對(duì)生育期影響也不同。

本文僅研究了水熱條件對(duì)錦州地區(qū)春玉米生育期長(zhǎng)度影響及其與生育期長(zhǎng)度的關(guān)系，未考慮品種、栽培措施和人為主觀因素等其他因素對(duì)生育期長(zhǎng)度的影響。作物的生育期是對(duì)整個(gè)水熱條件變化的一個(gè)適應(yīng)過(guò)程，除了水熱條件變化的不確定性和波動(dòng)性，還有人為的調(diào)整和控制。未來(lái)東北地區(qū)氣候變化將更大，極端天氣也將更頻繁出現(xiàn)，復(fù)合型災(zāi)害天氣也將頻繁出現(xiàn)，進(jìn)而造成水熱條件波動(dòng)變化也將變大，對(duì)春玉米生產(chǎn)的影響更大，需采用新方法和新技術(shù)進(jìn)行光、溫、水等氣候條件及農(nóng)藝措施、災(zāi)害和人為因素等方面的深入研究，合理利用作物自然資源，提升玉米對(duì)氣候極端變化的適應(yīng)性，確保遼寧地區(qū)春玉米的穩(wěn)產(chǎn)和高產(chǎn)。