張志高， 孫梓欣， 張秀麗， 郭可欣， 李卓婭， 郝海姣， 蔡茂堂

（1.安陽(yáng)師范學(xué)院資源環(huán)境與旅游學(xué)院，河南 安陽(yáng) 455000；2.河南大學(xué)地理與環(huán)境學(xué)院，河南 開封 475004；3.河南大學(xué)黃河文明與可持續(xù)發(fā)展研究中心，河南 開封 475001；4.西南石油大學(xué)地球科學(xué)與技術(shù)學(xué)院，四川 成都 610500；5.中國(guó)地質(zhì)科學(xué)院地質(zhì)力學(xué)研究所，北京 100081）

IPCC 第六次評(píng)估報(bào)告指出，全球地表溫度正以前所未有的速度上升，2011—2020 年全球平均表面溫度相比1850—1900年上升1.09 ℃，2001—2020年較工業(yè)化前增暖0.99 ℃［1］，北半球中高緯度地區(qū)氣溫上升更為明顯［2］。隨著全球變暖的持續(xù)發(fā)展，氣溫、光照、降水等氣候要素也相應(yīng)發(fā)生變化，進(jìn)而對(duì)作物生長(zhǎng)發(fā)育以及農(nóng)業(yè)生態(tài)環(huán)境產(chǎn)生顯著影響［3-4］。生長(zhǎng)季是溫度和土壤濕度條件適宜作物生長(zhǎng)的時(shí)期，對(duì)生態(tài)系統(tǒng)功能發(fā)揮有重要控制因素［5］。在全球氣候變暖背景下，生長(zhǎng)季變化會(huì)顯著影響生態(tài)系統(tǒng)及農(nóng)業(yè)作物生產(chǎn)［6-7］。因此，生長(zhǎng)季時(shí)空演變格局及其影響因素受到了國(guó)內(nèi)外學(xué)者的廣泛關(guān)注［8-11］。眾多學(xué)者分別以物候觀測(cè)［11-12］、歸一化植被指數(shù)（NDVI）［13］和地表氣溫［14］等方面資料進(jìn)行生長(zhǎng)季相關(guān)研究。物候觀測(cè)和NDVI 數(shù)據(jù)序列相對(duì)較短，而氣象站點(diǎn)氣溫?cái)?shù)據(jù)空間覆蓋范圍廣、持續(xù)時(shí)間長(zhǎng)，可用于長(zhǎng)時(shí)間尺度上氣候生長(zhǎng)季研究［15］。相關(guān)研究表明，在全球變暖背景下北半球大部分地區(qū)生長(zhǎng)季呈延長(zhǎng)趨勢(shì)［16-18］。中國(guó)大部分地區(qū)如內(nèi)蒙古［19］、華北地區(qū)［20］、青藏高原［21］和東北地區(qū)［22］等生長(zhǎng)季開始日提前、結(jié)束日推遲、長(zhǎng)度延長(zhǎng)，但區(qū)域性特征明顯，青藏高原地區(qū)生長(zhǎng)季增加最多，北方地區(qū)次之，南方地區(qū)最少［23］。在生長(zhǎng)季變化影響因素方面，Xia等［24］發(fā)現(xiàn)，多年代際氣候變化（multidecadal climate variability,MDV）與大西洋多年代際振蕩指數(shù)貢獻(xiàn)了1980—2009 年北半球生長(zhǎng)季延長(zhǎng)的53%；Dong 等［21］認(rèn)為，1960—2009 年青藏高原生長(zhǎng)季指數(shù)變化與海拔高度密切相關(guān)；吳蓓蕾等［25］研究發(fā)現(xiàn)，中國(guó)平均生長(zhǎng)季開始日提前和長(zhǎng)度延長(zhǎng)主要源于春季升溫，生長(zhǎng)季結(jié)束日推遲則與秋季變暖有關(guān)；Jiang 等［26］研究表明，新疆地區(qū)生長(zhǎng)季延長(zhǎng)主要源于結(jié)束日的推遲；Wu 等［27］揭示了風(fēng)速變化對(duì)秋季物候的影響，發(fā)現(xiàn)風(fēng)速減弱延遲了高緯度植被秋季落葉期；還有學(xué)者分析了生長(zhǎng)季氣候資源的變化［28］。

黃河流域橫跨我國(guó)北方東、中、西三大地理階梯，跨越干旱、半干旱和半濕潤(rùn)區(qū)，流域上、中、下游地形地貌、植被及氣候類型分異顯著，為我國(guó)重要的生態(tài)功能區(qū)。黃河流域現(xiàn)有耕地1.3×107hm2，糧食產(chǎn)量約占全國(guó)1/3，在我國(guó)農(nóng)業(yè)生產(chǎn)中占有重要地位。因此，本文基于1960—2020年逐日氣溫?cái)?shù)據(jù)資料，運(yùn)用線性趨勢(shì)估計(jì)、Mann-Kendall突變檢驗(yàn)和相關(guān)分析等方法對(duì)近61 a來黃河流域氣候生長(zhǎng)季指標(biāo)的時(shí)空演變進(jìn)行研究，并探討其影響因素，以期為準(zhǔn)確把握氣候變化對(duì)流域生態(tài)環(huán)境的影響、充分利用氣候資源以及保障糧食安全等方面提供科學(xué)依據(jù)和參考。

1 數(shù)據(jù)來源及研究方法

1.1 研究區(qū)概況

黃河流域（32°6′～41°48′N，95°50′～119°6′E）橫貫青藏高原、內(nèi)蒙古高原、黃土高原和華北平原4個(gè)地貌單元，包含青海、四川等9 省區(qū)，流域面積約7.52×105km2。黃河流域海拔在0～6241 m，西接青藏高原，中部為黃土高原區(qū)，東鄰華北平原，東西跨度約5464 km，整體地勢(shì)西高東低，起伏顯著。流域位于中緯度地帶，受大氣環(huán)流和季風(fēng)環(huán)流的影響，跨越干旱、半干旱、半濕潤(rùn)氣候區(qū)。流域多數(shù)地區(qū)降水在200～650 mm，且年內(nèi)分布不均，多集中在6—10 月，年均溫在-4～14 ℃，全年日照時(shí)數(shù)在2000～3300 h。為了研究黃河流域生長(zhǎng)季指標(biāo)的區(qū)域分布特征，將黃河流域進(jìn)一步劃分為上、中、下游地區(qū)，其中，上游地區(qū)包括久治、景泰和包頭等40 個(gè)氣象站點(diǎn)，中游地區(qū)包括武功、靖邊和神木等42 個(gè)氣象站點(diǎn)，下游地區(qū)包括新鄉(xiāng)、鄭州和沂源等7 個(gè)氣象站點(diǎn)。

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

1960—2020 年黃河流域89 個(gè)氣象站點(diǎn)逐日氣象數(shù)據(jù)來源于中國(guó)氣象科學(xué)數(shù)據(jù)共享服務(wù)網(wǎng)（http://data.cma.cn）中國(guó)地面氣候資料日值數(shù)據(jù)集，該數(shù)據(jù)集經(jīng)過極值檢驗(yàn)和RHtest 均一性檢驗(yàn)等方面的嚴(yán)格質(zhì)量控制。研究區(qū)域及代表站點(diǎn)分布如圖1所示。

圖1 黃河流域氣象站點(diǎn)分布Fig.1 Distribution of meteorology stations in the Yellow River Basin

1.3 生長(zhǎng)季指標(biāo)計(jì)算

生長(zhǎng)季的界定尚未統(tǒng)一，參考相關(guān)資料，本文將生長(zhǎng)季長(zhǎng)度（GSL）定義為一年中5 d 滑動(dòng)平均溫度大于5 ℃到最后一次連續(xù)5 d 滑動(dòng)平均溫度大于5 ℃之間的天數(shù)［29］。生長(zhǎng)季開始（GSS）為一年中第一次連續(xù)5 d滑動(dòng)平均氣溫大于5 ℃的第一天，生長(zhǎng)季結(jié)束（GSE）為最后一次連續(xù)5 d 滑動(dòng)平均氣溫大于5 ℃的最后一天。生長(zhǎng)季內(nèi)≥10 ℃活動(dòng)積溫（AT10）為生長(zhǎng)季內(nèi)≥10 ℃日平均氣溫的累積之和，對(duì)應(yīng)的天數(shù)之和為生長(zhǎng)季內(nèi)≥10 ℃活動(dòng)積溫天數(shù)（DT10）。GSS、GSE、GSL 和DT10 單位為d，AT10 單位為℃。

1.4 研究方法

采用Mann-Kendall 方法［30］對(duì)1960—2020 年黃河流域氣候生長(zhǎng)季指標(biāo)（GSS、GSE、GSL、AT10 和DT10）的突變特征進(jìn)行分析，利用信噪比對(duì)其突變點(diǎn)進(jìn)行檢驗(yàn)［31］。運(yùn)用Morlet小波分析法［32］分析黃河流域生長(zhǎng)季指標(biāo)的周期變化特征。

2 結(jié)果與分析

2.1 氣候生長(zhǎng)季時(shí)間變化

由圖2 可知，1960—2020 年黃河流域GSS 在63.36～88.43 d，平均為75.90 d（3月17日），2022年最小，1976 年最大。近61 a 來GSE 平均值為310.25 d（11 月6 日），在301.10～318.06 d 之間波動(dòng)，1981 年最小，2016年最大。近61 a GSL平均值為235.35 d，其中1976 年最小為219.74 d，2020 年最大為250.84 d。1960—2020 年AT10 在2925.51～3592.65 ℃，平均值為3265.91 ℃，1976 年最小，2016 年最大。近61 a 來DT10 平均為187.54 d，1976 年最小為172.92 d，2009年最大為206.73 d。從變化趨勢(shì)來看，近61 a來GSS呈顯著的提前趨勢(shì)（P<0.01），傾向率為-2.04 d·(10a)-1，GSE 以0.85 d·(10a)-1的傾向率呈明顯的推遲趨勢(shì)（P<0.05），GSL 則呈顯著的延長(zhǎng)趨勢(shì)（P<0.01），傾向率為2.88 d·(10a)-1。AT10 和DT10 均呈顯著增加趨勢(shì)，且通過了0.01 的顯著性檢驗(yàn)，傾向率分別為70.62 ℃·(10a)-1和3.26 d·(10a)-1。

圖2 1960—2020年黃河流域生長(zhǎng)季指標(biāo)變化Fig.2 Changes of growing season indicators in the Yellow River Basin from 1960 to 2020

年代際變化方面（表1），1960—1980 年代GSS偏大，距平為正，生長(zhǎng)季開始較晚，GSE距平為負(fù)，生長(zhǎng)季結(jié)束較早，GSL 距平為負(fù)，生長(zhǎng)季長(zhǎng)度相對(duì)較短，AT10 距平為負(fù)，相對(duì)較小。1990 年代后GSS 逐年代提前，距平轉(zhuǎn)為負(fù)值，GSE 距平轉(zhuǎn)正，生長(zhǎng)季末期推遲，GSL 距平轉(zhuǎn)為正值，生長(zhǎng)季長(zhǎng)度延長(zhǎng)，AT10距平轉(zhuǎn)正，生長(zhǎng)季內(nèi)≥10 ℃活動(dòng)積溫增加。1970 年代DT10 最小，2000 年代和2010 年代DT10 最大，生長(zhǎng)季內(nèi)≥10 ℃活動(dòng)積溫天數(shù)最多。

表1 黃河流域生長(zhǎng)季指標(biāo)年代際變化Tab.1 Interdecadal variation of growing season index in the Yellow River Basin

2.2 上、中、下游地區(qū)生長(zhǎng)季指標(biāo)時(shí)間變化

由表2 可知，近61 a 來流域上、中、下游地區(qū)GSS 分別為88.81 d、66.73 d 和53.63 d，上游地區(qū)生長(zhǎng)季開始最晚（3 月30），下游地區(qū)開始最早（2 月23），相差35 d；下游地區(qū)GSE 為334.03 d，上游為297.33 d，相差36.70 d，上游地區(qū)生長(zhǎng)季結(jié)束最早（10月24），下游地區(qū)結(jié)束最晚（11月30），相差37 d；下游地區(qū)GSL 為281.40 d，上游地區(qū)為208.02 d，上下游地區(qū)生長(zhǎng)季長(zhǎng)度相差73.38 d；下游地區(qū)AT10為4619.85 ℃，上游地區(qū)為2472.16 ℃，相差2147.69 ℃；DT10 在下游地區(qū)最大為232.88 d，上游地區(qū)最小為156.31 d，相差76.57 d。從變化趨勢(shì)來看，黃河流域上、中、下游地區(qū)各指標(biāo)變化趨勢(shì)一致，從變化幅度來看，下游地區(qū)GSS、GSE、GSL 和AT10 指標(biāo)變化幅度最大，上游地區(qū)次之，中游地區(qū)變幅最小，上游地區(qū)DT10 變化幅度最大，下游次之，中游地區(qū)最小。

表2 1960—2020年黃河流域生長(zhǎng)季指標(biāo)變化Tab.2 Changes of growing season index in the Yellow River Basin from 1960 to 2020

2.3 氣候生長(zhǎng)季空間分布特征

由圖3 可知，黃河流域GSS 的變化范圍為37～164 d，近61 a平均值為75.90 d，即3月17日，整體呈由東到西、由南向北逐漸推遲趨勢(shì)，流域東南部生長(zhǎng)季開始最早，日序在37～64 d，即2 月6 日—3 月5日，流域西南部達(dá)日、瑪多和清水河等站生長(zhǎng)季開始最晚，日序在129～164 d，即5 月9 日—6 月14 日。黃河流域GSE 在255～346 d，平均為310.25 d，即11月6 日，分布趨勢(shì)與始期相反，整體呈由東到西、由南向北逐漸提前的趨勢(shì)，流域東南部生長(zhǎng)季結(jié)束最晚，日序在327～346 d，即11月23日—12月12日，流域西南部瑪沁、達(dá)日、瑪多和清水河等站生長(zhǎng)季結(jié)束最早，日序在255～273 d，即9 月12 日—9 月30日。黃河流域GSL 在92～310 d，平均為235.35 d，整體呈由東到西，由南向北逐漸縮短趨勢(shì)，流域東南部GSL 最長(zhǎng)，在265 d 以上，西南部瑪多和清水河生長(zhǎng)季最短，在105 d 以下。黃河流域AT10 在33～5153 ℃，流域東南部積溫最高，在4142 ℃以上，西南部瑪多和清水河最低，在188.42 ℃以下。DT10 在3～251 d，流域東南部積溫天數(shù)最多，在203 d 以上，西南部瑪多和清水河最低，在20 d 以下，在空間分布上與AT10 的空間分布大致相似，即積溫天數(shù)越多，積溫越高，天數(shù)越少，積溫越低。

圖3 黃河流域氣候生長(zhǎng)季指標(biāo)的空間分布Fig.3 Spatial distribution of climate growing season index in the Yellow River Basin

2.4 氣候生長(zhǎng)季變化趨勢(shì)空間分布特征

1960—2020 年黃河流域生長(zhǎng)季指標(biāo)變化趨勢(shì)空間分布如圖4 和表3 所示，由圖4 可知，近61 a 來黃河流域GSS 整體呈提前趨勢(shì)，全部89 個(gè)氣象站點(diǎn)傾向率＜0，傾向率在-4.1～-0.2 d·(10a)-1之間，其中有71.9%站點(diǎn)通過了0.05顯著性檢驗(yàn)（表3），在空間分布上呈一定的區(qū)域集聚特征，流域東南部河南、山東地區(qū)與流域中部?jī)?nèi)蒙古西部、寧夏地區(qū)提前幅度較大，西部青海地區(qū)提前幅度較小，在-1.6 d·(10a)-1以內(nèi)。GSE 整體上呈延遲趨勢(shì)，傾向率在-0.6～2.6 d·(10a)-1之間，在89 個(gè)氣象站點(diǎn)中，有82 個(gè)站點(diǎn)呈遞增趨勢(shì)，且80.9%的站點(diǎn)達(dá)到了0.05 顯著性水平，流域中西部地區(qū)和東南部山東、河南地區(qū)推遲幅度較大，中東部山西地區(qū)推遲幅度較小。GSL 整體呈現(xiàn)出延長(zhǎng)趨勢(shì)，在0～5.5 d·(10a)-1之間，有88 個(gè)站點(diǎn)呈延長(zhǎng)趨勢(shì)，其中，80.9%的站點(diǎn)呈顯著延長(zhǎng)趨勢(shì)，流域東南部河南、山東地區(qū)和中西部?jī)?nèi)蒙古西部、寧夏地區(qū)形成兩個(gè)高值中心，傾向率在3.3 d·(10a)-1以上，流域中東部山西地區(qū)傾向率較小，在2.4 d·(10a)-1以下。AT10 整體呈現(xiàn)出上升趨勢(shì)，傾向率在2.4～148.6 ℃·(10a)-1，全部89 個(gè)站點(diǎn)呈上升趨勢(shì)，其中，93.3%站點(diǎn)通過了0.05顯著性檢驗(yàn)，流域中西部?jī)?nèi)蒙古西部、寧夏地區(qū)增幅較大，在88.4 ℃·(10a)-1以上，流域南部陜西地區(qū)、流域西部青海、四川地區(qū)增幅較小。DT10 整體呈上升趨勢(shì)，傾向率在-0.2～7.2 d·(10a)-1之間，88 個(gè)站點(diǎn)呈上升趨勢(shì)，其中，79.8%的站點(diǎn)通過了0.05 顯著性檢驗(yàn)，流域南部陜西地區(qū)增幅較小，在2.6 d·(10a)-1以下，流域西南部四川地區(qū)增幅較大，在4.3 d·(10a)-1以上。

表3 黃河流域生長(zhǎng)季指標(biāo)空間趨勢(shì)檢驗(yàn)Tab.3 Spatial trend test of growing season index in the Yellow River Basin

2.5 氣候生長(zhǎng)季突變特征

由圖5 可知，近61 a 來黃河流域GSS 正序列UF曲線呈波動(dòng)下降趨勢(shì)，反序列UB曲線呈上升趨勢(shì)，二者相交于1998 年，突變點(diǎn)在信度線內(nèi)，且通過了信噪比檢驗(yàn)，確定其突變點(diǎn)為1998 年，突變前生長(zhǎng)季日序平均為79.19 d（3 月20 日），突變后為70.45 d（3 月11 日），提前了9 d（圖5a）。GSE 的UF曲線在1994年后在0值之上，表明該時(shí)期UF曲線上升趨勢(shì)明顯，UB曲線呈波動(dòng)下降趨勢(shì)，UF和UB曲線存在多個(gè)交點(diǎn)，但均未通過信噪比檢驗(yàn)，故GSE 不存在突變點(diǎn)（圖5b）。GSL 的正序列UF和反序列UB曲線相交于2002 年，在信度線以內(nèi)，且通過了信噪比檢驗(yàn)，其突變點(diǎn)為2002 年，突變后GSL 平均為242.98 d，較突變前的231.90 d 顯著延長(zhǎng)了11.08 d（圖5c）。AT10和DT10的正序列UF和反序列UB曲線均相交于1998 年，均在信度線內(nèi)，且通過了信噪比檢驗(yàn)，故AT10 和DT10 的突變點(diǎn)為1998 年，突變后AT10 和DT10 分別較突變前增加了288.92 ℃和13.49 d。

2.6 生長(zhǎng)季周期變化特征

1960—2020年黃河流域生長(zhǎng)季指標(biāo)Morlet小波分析結(jié)果如表4 所示（圖略），由表4 可知，GSS 存在28 a、19 a與13 a左右的周期變化，其中，28 a左右的周期為第一主周期，19 a 左右周期為第二主周期，13 a 左右周期為第三主周期。近61 a 來黃河流域GSE、GSL、AT10 和DT10 均存在28 a 左右的主周期變化，而其他周期不明顯。

表4 1960—2020 年黃河流域生長(zhǎng)季指標(biāo)周期變化特征Tab.4 Characteristics of periodic changes of index in the Yellow River Basin during the growing season from 1960 to 2020

2.7 生長(zhǎng)季與溫度變化的關(guān)系趨勢(shì)分析

由表5 可知，黃河流域上游地區(qū)生長(zhǎng)季指標(biāo)與溫度變化相關(guān)性最強(qiáng)，中游次之，下游最弱。黃河流域平均GSS 與溫度變化顯著負(fù)相關(guān)，均通過了0.01的顯著性檢驗(yàn)，與春季平均溫度負(fù)相關(guān)性最強(qiáng)，相關(guān)系數(shù)為-0.65，說明GSS 提前與春季變暖有關(guān)。分流域來看，黃河流域上、中、下游地區(qū)均與溫度變化呈負(fù)相關(guān)關(guān)系，其中，與春季溫度負(fù)相關(guān)性最強(qiáng)，相關(guān)系數(shù)分別為-0.83、-0.55 和-0.34，相關(guān)性依次減弱。黃河流域平均GSE 與秋季平均溫度相關(guān)性最強(qiáng)，說明GSE 延遲與秋季變暖有關(guān)，上、中、下游地區(qū)GSE 均與秋季溫度相關(guān)性最強(qiáng)，相關(guān)系數(shù)分別為0.82、0.66 和0.42。黃河流域平均GSL 與秋季平均溫度相關(guān)性最強(qiáng)，說明GSL 延長(zhǎng)與秋季變暖有關(guān)，上游地區(qū)與春季溫度相關(guān)性最強(qiáng)，相關(guān)系數(shù)達(dá)0.80，中游地區(qū)GSL 延長(zhǎng)主要與秋季增溫有關(guān)，相關(guān)系數(shù)為0.62，下游地區(qū)主要與春季增溫有關(guān)，相關(guān)系數(shù)為0.41。

表5 黃河流域生長(zhǎng)季指標(biāo)與各季節(jié)平均溫度的相關(guān)關(guān)系Tab.5 Correlation between growing season index and seasonal average temperature in the Yellow River Basin

3 討論

1960—2020 年黃河流域GSS 呈提前趨勢(shì)，GSE呈推遲趨勢(shì)，GSL呈延長(zhǎng)趨勢(shì)，這與全國(guó)生長(zhǎng)季變化趨勢(shì)一致［25］。從變化幅度來看，黃河流域GSS、GSE和GSL 變化趨勢(shì)分別為-2.04 d·(10a)-1、0.85 d·(10a)-1和2.88 d·(10a)-1，黃河流域GSS 提前幅度明顯大于全國(guó)平均水平及其他多數(shù)地區(qū)，小于華北地區(qū)［20］，而GSE 延遲幅度要小于全國(guó)平均水平，GSL 延長(zhǎng)幅度大于全國(guó)及其他大部分區(qū)域，小于華北［20］和青藏高原地區(qū)［21］，這可能與黃河流域1960—2020年增溫趨勢(shì)0.36 ℃·(10a)-1顯著高于全國(guó)增溫速率0.25 ℃·(10a)-1［33］有關(guān)。近61 a 來黃河流域GSS 提前了12.4 d，GSE 推遲了5.2 d，GSL 延長(zhǎng)了17.6 d，GSS 貢獻(xiàn)了生長(zhǎng)季延長(zhǎng)的70.5%，GSE 貢獻(xiàn)了29.5%，GSS 的提前對(duì)生長(zhǎng)季延長(zhǎng)的影響更為重要。

黃河流域GSS 提前主要與春季增溫有關(guān)，GSE延遲主要與秋季增溫有關(guān)，這與吳蓓蕾等［25］對(duì)全國(guó)生長(zhǎng)季的研究結(jié)果一致，黃河流域GSL 的延長(zhǎng)主要與秋季增溫有關(guān)，這與全國(guó)平均GSL 延長(zhǎng)主要源于春季升溫的結(jié)論不同，這可能與各區(qū)域氣溫變化不同有關(guān)。

在突變特征方面，黃河流域生長(zhǎng)季指標(biāo)分別于1998年和2002年發(fā)生突變，這與華北地區(qū)［20］研究結(jié)果一致，研究表明，全球氣候狀態(tài)轉(zhuǎn)變發(fā)生于20 世紀(jì)80年代中后期［34］，氣候狀態(tài)的轉(zhuǎn)變導(dǎo)致氣候變率增大及極端氣候事件增加，導(dǎo)致生長(zhǎng)季始期提前，生長(zhǎng)季末期延遲和生長(zhǎng)季延長(zhǎng)。

本文利用黃河流域89 個(gè)氣象站點(diǎn)氣溫?cái)?shù)據(jù)來研究生長(zhǎng)季指標(biāo)變化，還不能較好地反應(yīng)區(qū)域內(nèi)植被生長(zhǎng)動(dòng)態(tài)，后續(xù)的研究中應(yīng)結(jié)合物候觀測(cè)和歸一化植被指數(shù)（NDVI）等進(jìn)一步探究流域生長(zhǎng)季指標(biāo)變化。此外，本研究探討了黃河流域生長(zhǎng)季指標(biāo)與各季溫度的關(guān)系，城市化進(jìn)程、土地利用變化及環(huán)流因子對(duì)生長(zhǎng)季的影響有待進(jìn)一步研究。

4 結(jié)論

（1）1960—2020 年 黃 河 流 域GSS 顯 著 提 前[-2.04 d·(10a)-1]，GSE 呈推遲趨勢(shì)[0.85 d·(10a)-1]，GSL 顯著延長(zhǎng)[2.88 d·(10a)-1]，AT10 和DT10 顯著增加，傾向率分別為70.62 ℃·(10a)-1和3.26 d·(10a)-1，但區(qū)域差異較大，下游GSS 開始最早（2 月23），上游最晚（3 月30），上游GSE 結(jié)束最早（10 月24），下游最晚（11月30），下游GSL最長(zhǎng)為334.03 d，上游最短為297.33 d，下 游AT10 為4619.85 ℃，上 游 為2472.16 ℃，下游DT10 最大為232.88 d，上游最小為156.31 d。年代際變化來看，1990 年代以來，黃河流域GSS 逐年代提前，GSE 推遲，GSL 顯著延長(zhǎng)，AT10和DT10明顯上升。

（2）近61 a 黃河流域整體上由東到西，由南向北GSS 逐 漸 推 遲，GSE 逐 漸 提 前，GSL 逐 漸 縮 短，AT10和DT10逐漸減少。黃河流域上、中、下游地區(qū)生長(zhǎng)季指標(biāo)變化趨勢(shì)一致，下游地區(qū)變化幅度最大，上游地區(qū)次之，中游地區(qū)變幅最小。

（3）M-K 突變分析表明，近61 a 來黃河流域GSS、AT10 和DT10 于1998 年發(fā)生突變，突變后GSS提前了9 d、AT10 和DT10 分別增加了288.92 ℃和13.49 d。GSE 不存在明顯突變點(diǎn)。GSL 于2002 年發(fā)生突變，突變后GSL 延長(zhǎng)了11.08 d。Morlet 小波分析表明，近61 a 來黃河流域生長(zhǎng)季指標(biāo)存在28 a左右的周期變化。

（4）近61 a 來，黃河流域生長(zhǎng)季長(zhǎng)度延長(zhǎng)了17.6 d，GSS貢獻(xiàn)了生長(zhǎng)季延長(zhǎng)的70.5%，GSE貢獻(xiàn)了29.5%。相關(guān)分析表明，黃河流域GSS 提前主要與春季升溫有關(guān)，GSE延遲主要源于秋季增溫，上游和下游地區(qū)GSL 延長(zhǎng)主要源于春季增溫，中游地區(qū)GSL延長(zhǎng)主要與秋季變暖有關(guān)。