關(guān)鍵詞暖冬;冬季氣候變化;Mann-Kendall檢驗(yàn);相關(guān)性分析
中圖分類號(hào)P467 文獻(xiàn)標(biāo)識(shí)碼A 文章編號(hào) 1007-7731(2025)17-0109-06
DOI號(hào) 10.16377/j.cnki.issn1007-7731.2025.17.027
Analysis of warm winter climate characteristics in Fengyang County
YUAN Xuesuo WU JinxiangWANG Gang GUO Yang (Fengyang Meteorological Bureau,F(xiàn)engyang 2331Oo, China)
AbstractThe characteristics of warm winters over the past 67 years in the study area were analyzed based on meteorologicalobservationdatafromtheFengyang National Basic Meteorological StationinAnhui,China,from1957to 2023,including winter average temperature (T) ,average maximum temperature (MTh),and average minimum temperature (MTe).Methodssuchascorrelationcoeffientanalysis,Mann-Kendalltest,and linear trend estmation were employed.Theresults indicatethatthe winteraverage temperaturesforstrong warm,warm,normal,andcold winter types in the study area were 4.3-4.7°C, 3.2-4.2°C ,
,and
,respectively. Correlation analysis revealedthatsunshine duration (S)wasnegativelycorelatedwith MTeandprecipitation (R),with correlatin coeffcients of-0.471and-O.647,respectively,while most other climatic elements showed positive correlations.The Mann-Kendalltrend and abrupt change tests demonstrated that from 1957 to 2O23,winter T,MTh,MTe,and R in the studyarea exhibited significant increasing trends,whereas S showed a significant decreasing trend. Except for Th abruptchanges were detected in all other climatic elements.Lineartrend estimation indicated thatthe climatic tendency rates for winter T , R ,and S were 0.166°C/10 a, 8.106mm/10 a, and -16.939 h/10 a, respectively. The findings provide a meteorological reference for agricultural production.
Keywords warm winter; winter climate change; Mann-Kendall test; correlation analysis
暖冬是指冬季平均氣溫高于氣候平均值的氣候現(xiàn)象,對(duì)農(nóng)業(yè)、生態(tài)以及經(jīng)濟(jì)等有直接或間接影響,不同區(qū)域暖冬氣候特征的研究逐漸成為研究人員的關(guān)注熱點(diǎn)。如阮翠冰等利用線性傾向估計(jì)、Mann-Kendall突變檢驗(yàn)和小波分析等方法分析了寧德市暖冬現(xiàn)象特征;王秀萍等采用線性傾向估計(jì)、Mann-Kendall突變檢驗(yàn)、百分位法定義極端事件閾值和合成分析等方法,分析了大連地區(qū)冬季冷暖變化趨勢(shì)及異常特征,并探討了冷冬、暖冬年大氣環(huán)流異常成因;普布桑姆等4采用線性傾向估計(jì)、相關(guān)系數(shù)分析、Mann-Kendall檢驗(yàn)等方法,分析了近32a林芝暖冬指數(shù)變化特征及其對(duì)冬小麥生育期的影響;謝睿恒等5利用1961—2020年湖南省的冬季逐日氣溫資料,分析了該地區(qū)暖冬事件的氣候特征及成因;紀(jì)詩(shī)璇等%分析了廊坊地區(qū)1971—2020年冬季平均氣溫、平均最高氣溫、年平均最低氣溫變化趨勢(shì)、冷暖冬事件變化特點(diǎn),并探究了冬季氣溫變化對(duì)流感病傳播的影響。本研究基于鳳陽(yáng)國(guó)家基本氣象站長(zhǎng)序列氣象觀測(cè)資料,采用相關(guān)系數(shù)分析、線性傾向估計(jì)、Mann-Kendall趨勢(shì)檢驗(yàn)和突變檢驗(yàn)等方法,分析了1957—2023年研究區(qū)暖冬氣候特征,對(duì)冬季平均氣溫、降水量、日照時(shí)數(shù)等主要?dú)夂蛞剡M(jìn)行時(shí)變特征分析,為當(dāng)?shù)剞r(nóng)業(yè)、生態(tài)等產(chǎn)業(yè)發(fā)展提供參考。
1材料與方法
1.1 數(shù)據(jù)來(lái)源
氣象資料來(lái)源于1957—2023年研究區(qū)國(guó)家基本氣象站,統(tǒng)計(jì)每年冬季主要?dú)夂蛞兀唧w包括:冬季平均氣溫T(12月至次年2月每日平均氣溫的算術(shù)平均值)冬季平均最高氣溫MT(同時(shí)段每日最高氣溫的算術(shù)平均值)冬季平均最低氣溫MTe(同時(shí)段每日最低氣溫的算術(shù)平均值)冬季極端最高氣溫 Th (同時(shí)段每日最高氣溫的最大值)冬季極端最低氣溫 Te (同時(shí)段每日最低氣溫的最小值)冬季降水量R(同時(shí)段每日降水量的累計(jì)值)冬季日照時(shí)數(shù)S(同時(shí)段每日日照時(shí)數(shù)的累計(jì)值)。
1.2試驗(yàn)方法
1.2.1冬季氣候概況及相關(guān)性分析 參照GB/T21983一2020《暖冬等級(jí)》,將研究區(qū)冬季劃分為“強(qiáng)暖、偏暖、正常、偏冷\"4種年型(表1),對(duì)冬季平均氣溫等7個(gè)氣候要素,利用SPSS軟件計(jì)算相關(guān)系數(shù)分析氣候要素之間的相關(guān)性,研究其時(shí)變特征。
表11957—2023年研區(qū)冬季年型劃分

續(xù)表11957—2023年研究區(qū)冬季年型劃分

續(xù)表11957—2023年研究區(qū)冬季年型劃分

1.2.2趨勢(shì)分析 利用Mann-Kendall檢驗(yàn)法對(duì)冬季氣候要素 T,MTh,MTe,Th,Te 進(jìn)行趨勢(shì)分析 [7-8] 。檢驗(yàn)統(tǒng)計(jì)量S計(jì)算如式(1),當(dāng)時(shí)間序列 Δngt; 10時(shí),統(tǒng)計(jì)量 Z 計(jì)算如式(2)。


式中, sign 為符號(hào)函數(shù)。當(dāng) Di-Djlt;0?=0?gt; 0時(shí), sign(Di-Dj) 分別取 -1.0,1;Zgt;0 表示增加趨勢(shì), Zlt;0 表示減少趨勢(shì)。 Z 的絕對(duì)值 ?1.28,?1.64,? 2.32時(shí)表示分別通過(guò)了信度 90%95%.99% 顯著性檢驗(yàn)。
1.2.3突變檢驗(yàn) 利用Mann-Kendall檢驗(yàn)法對(duì)冬季氣候要素 T,MTh,MTe,Th,Te,R,S 進(jìn)行突變檢驗(yàn)[7-8],設(shè)有氣候要素時(shí)間序列: ΦD2,D3,…,Dn ,構(gòu)造一秩序列 ri,ri 表示 Digt;Dj(1?j?i) 的樣本累積數(shù)Sk ,計(jì)算如式(3)。

式中, ri 表示 Digt;Dj(I?i?j) 的累加數(shù),當(dāng) Digt;Dj 時(shí), ri=1 ;當(dāng) Di?Dj 時(shí) ξ,ri=0(j=1,2,…,i) 。
在時(shí)間序列隨即獨(dú)立假定下,定義統(tǒng)計(jì)量列UFκ,E(Sk),Var(Sk) ,計(jì)算如式(4)~(6)。



式中, E(Sk),Var(Sk) 是累積數(shù) Sk 的均值和方差。當(dāng) k=1 時(shí), UF1=0 。 UFk 為標(biāo)準(zhǔn)正態(tài)分布,給定一顯著水平 α ,得到臨界值 Uα ,如 α=0.05 ,臨界值 Uα= ±1.96 。當(dāng) |UFk|gt;Uα ,表明序列存在一個(gè)明顯的增長(zhǎng)或減少趨勢(shì),通過(guò)信度檢驗(yàn)可知其是否具有趨勢(shì)。將此方法引用到反序列中,再重復(fù)上述計(jì)算過(guò)程,并使計(jì)算值乘以-1,得到 UBk 。分析 UFk 和 UBk 曲線圖,當(dāng) UFkgt;0 時(shí),序列呈上升趨勢(shì); UFklt;0 時(shí),呈下降趨勢(shì),若 UFk 超過(guò)臨界值則說(shuō)明上升或下降趨勢(shì)顯著。當(dāng) UFk 曲線與 UBk 曲線出現(xiàn)交點(diǎn)時(shí),且交點(diǎn)在臨界值之間,則突變開始的時(shí)間就是交點(diǎn)所對(duì)應(yīng)的時(shí)間。
1.2.4線性趨向分析 選取對(duì)農(nóng)業(yè)生產(chǎn)影響較大的氣候要素 (T,R,S) 作為因變量, (y) 年份記為自變量 (x) 。建立一元線性回歸方程如式(7)。
y=ax+b
式中, Ωa 為回歸系數(shù), b 為常數(shù), a×10 為氣候傾向率。
2 結(jié)果與分析
2.1冬季氣候概況
由表2可知,偏暖年型的 T 在
,強(qiáng)暖年型的 T 在 4.3~4.7'C ;正常年型的 T 在 2.4~ 3.1°C ,偏冷年型的 T 在 0.4~2.3°C 。強(qiáng)暖年型中,MTh最高可達(dá) 11.0°C ;偏暖年型中,MTh最高可達(dá)10.0° ;正常年型中,MTh在
;偏冷年型中,MTh最低為 5.0°C 。強(qiáng)暖年型中,MTe最高為1.1°C ;偏暖年型中,MTe最高為 0.1°C ;正常年型中,MTe在 -2.5~0°C ;偏冷年型中,MTe最低值為-4.1°C 。強(qiáng)暖年型中, Th 最高為 27.1°C ;偏暖年型中, Th 最高為 25.9°C ;正常年型中, Th 在 15.4~ 24.6°C ;偏冷年型中, Th 最低為 12.8°C ,最高為24.6°C ,與正常年型最高值持平。強(qiáng)暖年型中, Te 最高為 -5.7° ;偏暖年型中, Te 最高值 -7.3°C ;正常年型中, Te 在
;偏冷年型中, Te 最低為
,此值高于正常年型的最低值
,這是強(qiáng)冷空氣活動(dòng)的年際差異所致,因此正常年型也需注意短時(shí)間的極端低溫。強(qiáng)暖年型中, R 最高為175.8mm ;偏暖年型中, R 最高為 155.3mm ;正常年型中, R 在 26.4~195.2mm ;偏冷年型中, R 最低為18.5mm ;不同年型 R 的變幅較大,其中正常年型 R 的變幅最大。強(qiáng)暖年型中, s 最大為 529.3h ;偏暖年型中, s 最大為 578.9h ;正常年型中, s 在 212.4~ 643.3h ;偏冷年型中, s 最小為 324.0h ;不同年型 s 的變幅較大,其中正常年型 s 的變幅最大。
表2研究區(qū)冬季氣候概況

2.2 相關(guān)性分析
由表3可知,T與 MTh,MTe,Th,Te 的相關(guān)性在0.01水平具有統(tǒng)計(jì)學(xué)意義,其中,T與MTe相關(guān)系數(shù)最大,為 0.850;MTh 與 MTe,Th,Te 的相關(guān)性在0.01水平具有統(tǒng)計(jì)學(xué)意義,與 s 的相關(guān)性在0.05水平具有統(tǒng)計(jì)學(xué)意義,其中,MTh與Th相關(guān)系數(shù)最大,為0.485;MTe與 Te,R,S 的相關(guān)性在0.01水平具有統(tǒng)計(jì)學(xué)意義,其中,MTe與 R 相關(guān)系數(shù)最大,為 0.508:R 與S的相關(guān)系數(shù)為-0.647,在0.01水平具有統(tǒng)計(jì)學(xué)意義。說(shuō)明研究區(qū)冬季主要?dú)夂蛞刂g存在一定的相關(guān)性。
表3研究區(qū)冬季氣候要素之間的相關(guān)性

注:**和*分別表示相關(guān)性在0.01和0.05水平具有統(tǒng)計(jì)學(xué)意義。
2.3 趨勢(shì)分析
由表4可知,1957—2023年,研究區(qū)冬季氣候要素 T,MTh,MTe,Th,Te,R 呈上升趨勢(shì),S呈下降趨勢(shì)。其中 T,MTh,S 的變化趨勢(shì)在0.01水平具有統(tǒng)計(jì)學(xué)意義, MTe,R 上升趨勢(shì)在0.05水平具有統(tǒng)計(jì)學(xué)意義。這種氣候變化可能與溫室氣體濃度增加有關(guān)。
表4研究區(qū)冬季氣候要素Mann-Kendall趨勢(shì)分析

注:**和*分別表示變化趨勢(shì)在0.01和0.05水平具有統(tǒng)計(jì)學(xué)意義。
2.4 突變檢驗(yàn)
由表5可知,1957—2023年,研究區(qū)冬季氣候要素 T,MTh,MTe,Te,R,S 的變化趨勢(shì)均存在一定突變現(xiàn)象。其中,S始終為下降趨勢(shì),其余要素突變前多為下降或持平趨勢(shì);突變后各要素的上升趨勢(shì)逐漸明顯; Th 未發(fā)生突變。其中,光、溫、水通常被認(rèn)為是影響農(nóng)業(yè)生產(chǎn)的基本氣候要素,故以圖1更直觀、全面地展示 T,R,S 的突變檢驗(yàn)過(guò)程,其余要素的檢驗(yàn)圖省略。
表5研究區(qū)冬季氣候要素Mann-Kendall突變檢驗(yàn)

圖11957—2023年研究區(qū)氣候要素Mann-Kendall突變檢驗(yàn)

2.5 線性傾向分析
由圖2可知,1957—2023年,研究區(qū)冬季平均氣溫氣候傾向率為 0.166°C/10a,T 與時(shí)間的相關(guān)系數(shù)為0.342(由圖中 R2 值開平方求得,下同),通過(guò)信度為 99% 的顯著性檢驗(yàn);與表4中Mann-Kendall趨勢(shì)檢驗(yàn)結(jié)果相一致;研究區(qū)冬季降水量氣候傾向率為 8.106mm/10a,R 與時(shí)間的相關(guān)系數(shù)為0.365,通過(guò)信度為 99% 的顯著性檢驗(yàn),與表4中Mann-Kendall趨勢(shì)檢驗(yàn)結(jié)果基本一致,信度更高;研究區(qū)冬季日照時(shí)數(shù)氣候傾向率為 -16.939h/10a,S 與時(shí)間的相關(guān)系數(shù)為-0.417,通過(guò)信度為 99% 的顯著性檢驗(yàn),與表4中Mann-Kendall趨勢(shì)檢驗(yàn)結(jié)果完全一致。
圖21957一2023年研究區(qū)氣候變化趨勢(shì)

3結(jié)論
本研究利用研究區(qū)1957—2023年的冬季平均氣溫 T, 冬季平均最高氣溫MTh等氣候資料,參照GB/T21983—2020《暖冬等級(jí)》,利用Mann-Kendall法進(jìn)行趨勢(shì)分析和突變檢驗(yàn),對(duì)氣候變化特征進(jìn)行分析,得出以下結(jié)論。
(1)研究區(qū)冬季強(qiáng)暖年型的 T 在 4.3~4.7'C ;偏暖年型的 T 在 3.2~4.1'C ;正常年型的 T 在 2.4~ 3.1°C ;偏冷年型的 T 在 0.4~2.3°C 。在強(qiáng)暖、偏暖年型中MTh最高分別為11.0和 10.0°C;MTe 最高分別為1.1和 0.1°C 。研究區(qū)冬季冷空氣活動(dòng)頻繁,不同年型均需關(guān)注 -10°C 以下的極端最低氣溫,做好防寒保暖工作。
(2)研究區(qū)冬季 T 與 MTh,MTe,Th,Te 呈正相關(guān),與 R,S 相關(guān)性無(wú)統(tǒng)計(jì)學(xué)意義;MTh與 s 呈正相關(guān);MTe與 R 呈正相關(guān),與 s 呈負(fù)相關(guān); R 與 s 呈負(fù)相關(guān)。說(shuō)明各氣候要素之間有一定的相關(guān)性,可以互相影響。
(3)1953—2023年,研究區(qū)冬季 T,MTh,MTe 、Th,Te,R 總體均呈現(xiàn)上升趨勢(shì), s 呈現(xiàn)下降趨勢(shì);其中 T,R 和 s 的傾向率分別為 0.166°C/10: 小8.106mm/10 a和 -16.939h/10a;T.MTh.MTe.Te.R 的時(shí)間序列均存在突變點(diǎn),突變點(diǎn)之后的上升趨勢(shì)逐漸明顯; s 也存在突變,但始終是下降趨勢(shì)。
本研究揭示了研究區(qū)近67年冬季4種年型主要?dú)夂蛞氐幕厩闆r及其相互關(guān)系,有助于對(duì)冬季光、溫、水等氣候要素配置情況的總體把握,為農(nóng)業(yè)等相關(guān)領(lǐng)域研究提供氣候參考。
參考文獻(xiàn)
[1]馬鋒敏,陳麗娟,李想,等.中國(guó)冬季氣溫不同年代際的季節(jié)內(nèi)變化特征及成因分析[J].大氣科學(xué),2022,46(3):573-589.
[2]阮翠冰,陳小英,黃文霖.1971—2022年寧德市暖冬現(xiàn)象特征分析[J]寧德師范學(xué)院學(xué)報(bào)(自然科學(xué)版),2023,35(1):86-91.
[3]王秀萍,金巍,侯彥澤,等.1971—2021年大連地區(qū)冬季冷暖變化特征及大氣環(huán)流異常成因分析[J]氣象與環(huán)境學(xué)報(bào),2024,40(5):30-39.
[4]普布桑姆,次仁旺姆,杜軍,等.西藏林芝暖冬指數(shù)對(duì)冬小麥生長(zhǎng)發(fā)育的影響[J].高原山地氣象研究,2024,44(3):145-150.
[5]謝睿恒,涂瀟琴,羅瀟,等.1961—2020年湖南省暖冬事件氣候特征分析[J].農(nóng)業(yè)災(zāi)害研究,2023,13(3):67-70.
[6]紀(jì)詩(shī)璇,沈芳,白千川,等.河北廊坊1971—2020年冷暖冬氣候變化特征分析[J].農(nóng)業(yè)災(zāi)害研究,2022,12(7):127-129.
[7]袁學(xué)所,夏志芬,葛慶云,等.鳳陽(yáng)近61年干燥指數(shù)變化的Mann-Kendall檢驗(yàn)[J].科技通報(bào),2019,35(8):151-154.
[8]陳中平,徐強(qiáng).Mann-Kendall檢驗(yàn)法分析降水量時(shí)程變化特征[J].科技通報(bào),2016,32(6):47-50.
(責(zé)任編輯:胡立萍)