“雙碳”目標(biāo)背景下山西省地溫時(shí)空變化分析

范秀平 戴翠賢 賈斯 蘇睿杰 方紅娟 岳改慧

摘 要：【目的】探索山西省近40 a地溫變化規(guī)律，為其地?zé)衢_(kāi)發(fā)、碳達(dá)峰路徑選擇提供參考?！痉椒ā窟x取1981—2020年山西省108個(gè)臺(tái)站地表溫度數(shù)據(jù)，利用氣候傾向率、突變檢驗(yàn)、小波分析、經(jīng)驗(yàn)正交函數(shù)分解等方法進(jìn)行分析?！窘Y(jié)果】山西省近40 a地溫呈明顯上升趨勢(shì)，上升速率為0.499 ℃/10 a。地溫呈現(xiàn)“南高北低”分布，高緯地區(qū)低于低緯地區(qū)，高海拔地區(qū)低于低海拔地區(qū)。呂梁山脈走向?qū)Φ販厣涤绊懨黠@，西北部升溫速率最快。【結(jié)論】山西省地溫與時(shí)間相關(guān)系數(shù)高值區(qū)基本與趨勢(shì)系數(shù)高值區(qū)重合。地溫同位相變化占主導(dǎo)地位，即趨于整體一致的偏高或偏低變化，并在1999年后趨向于更明顯的升溫趨勢(shì)。近40 a存在18 a、12 a、6 a左右的周期震蕩，未來(lái)一段時(shí)間山西地溫趨向夏季更高、秋季更低。地溫變化與氣溫變化密切相關(guān)，約2/3站點(diǎn)地溫變化與降水量呈負(fù)相關(guān)；約1/3站點(diǎn)地溫變化與日照時(shí)數(shù)顯著相關(guān)。

關(guān)鍵詞：山西?。坏乇頊囟?；時(shí)空變化

中圖分類號(hào)：P423???? 文獻(xiàn)標(biāo)志碼：A???? 文章編號(hào)：1003-5168（2024）08-0097-04

DOI：10.19968/j.cnki.hnkj.1003-5168.2024.08.019

Analysis of Spatiotemporal Changes of Land Surface Temperature in Shanxi for Carbon Peaking and Carbon Neutrality Goals

FAN Xiuping1 DAI Cuixian2 JIA Si3 SU Ruijie1 FANG Hongjuan2 YUE Gaihui1

（1.Shanxi Meteorological Information Center， Taiyuan 030006， China; 2.Sanmenxia Meteorological Bureau of Henan Province， Sanmenxia 472000， China; 3.Xinxiang Changyuan Meteorological Bureau of Henan Province， Changyuan 453400， China）

Abstract： [Purposes] This paper aims to explore the variation rule of land surface temperature in Shanxi Province during the past 40 years， and provide reference for geothermal development and carbon peak path selection in Shanxi Province. [Methods] This paper selected the land surface temperature data of 108 stations in Shanxi from 1981 to 2020，used the methods of climate propensity rate analysis，M～K mutation detection， wavelet analysis etc analyze land surface temperature. [Findings] The land surface temperature of nearly 40 a in Shanxi showed a significant upward trend， with an increase rate of 0.499 ℃/10 a. The land surface temperature high in the south and low in the north， high latitudes are clearly lower than those in low latitudes， high altitudes are significantly lower than low altitudes. The trend of the Lüliang Mountain has a significant impact on the rise and fall of land surface temperature. The northwest has the fastest rate of warming. [Conclusions] The high value area of correlation coefficient of land surface temperature and time in the whole province basically coincides with the high value area of trend coefficient. The phase change of the land surface temperature is dominant， the high or low change tends to be consistent as a whole， and tends to be a more obvious warming trend after 1999. there are periodic oscillations of about 18 a， 12 a and 6 a， and the land surface temperature in Shanxi has a tendency to be higher in summer and lower in autumn in the future. The land surface temperature variation is closely related to air temperature variation. About 2/3 of the stations showed negative correlation between land surface temperature and precipitation. About 1/3 of the sites have significant correlation between land surface temperature and sunshine duration.

Keywords： Shanxi Province; the land surface temperature; spatiotemporal changes

0 引言

氣候變化對(duì)人類社會(huì)可持續(xù)發(fā)展的嚴(yán)峻挑戰(zhàn)是各個(gè)國(guó)家面臨的重要議題。實(shí)現(xiàn)“雙碳”目標(biāo)和應(yīng)對(duì)氣候變化，需要多部門(mén)、多學(xué)科的廣泛合作與交融。地表溫度的變化趨勢(shì)、變化頻率對(duì)區(qū)域氣候［1-2］、農(nóng)作物布局、農(nóng)業(yè)生產(chǎn)、土壤養(yǎng)分變化、氣候產(chǎn)量［3］、海綿城市建設(shè)、區(qū)域地質(zhì)災(zāi)害［4-5］等有著重要影響，甚至與輸電線路運(yùn)轉(zhuǎn)［6］等密切相關(guān)。無(wú)論是基于全國(guó)地表溫度數(shù)據(jù)變化研究，還是基于區(qū)域地表溫度分析，其氣候空間特征和突變特征都可為助力經(jīng)濟(jì)轉(zhuǎn)型升級(jí)、糧食生產(chǎn)安全等提供有效參考。本研究利用山西省1981—2020年108個(gè)氣象站地表溫度資料進(jìn)行分析，探究山西省近40 a地溫變化規(guī)律，為其地?zé)衢_(kāi)發(fā)、碳達(dá)峰路徑選擇提供參考。

1 資料與方法

1.1 氣候傾向率

假設(shè)：yt = a+bx，（t = 1，2，3，4[…]n） （1）

式中：yt為對(duì)應(yīng)年份地表溫度；t為時(shí)間序列；x為年份值；b為回歸系數(shù)，b>0說(shuō)明y隨時(shí)間增加呈上升趨勢(shì)，b<0說(shuō)明y隨時(shí)間增加呈下降趨勢(shì)。

1.2 小波分析方法

對(duì)于時(shí)間序列函數(shù)f（x），小波變換定義為式（2）。

wf（a，b）= [a12（x）ψ*.x-badx] （2）

式中：wf（a，b）為小波系數(shù)；a為伸縮因子，決定小波寬度；b為平移因子。

1.3 突變檢驗(yàn)方法

采用Mann-Kendall法進(jìn)行突變點(diǎn)檢驗(yàn)。對(duì)于具有n個(gè)樣本量的時(shí)間序列x，構(gòu)造一個(gè)序列，如式（3）。

Sk = [i=1kri k=]2，3，…，n （3）

其中：ri = [0，當(dāng)xi≤xj1，當(dāng)xi>xj?? ] j=1，2，…，n

序列Sk是第i時(shí)刻數(shù)值大于j時(shí)刻數(shù)值個(gè)數(shù)的累計(jì)數(shù)。在時(shí)間序列隨機(jī)獨(dú)立的假設(shè)下，定義UFk為式（4）。

UFk = [Sk-E（Sk）VAR（Sk）]k=1，2，…，n （4）

2 分析結(jié)果

2.1 時(shí)間分布特征

山西省近40 a平均地溫為10.7（1984年）～13.6 ℃（2019年），平均為12.2 ℃，40 a地溫平均值比氣候平均值（1991—2020年）的12.5 ℃偏低0.3 ℃（如圖1所示），最大正距平為1.1 ℃，最大負(fù)距平為-1.7 ℃，近40 a上升速率為0.499 ℃/10 a。

極端最高地溫平均值為60.2（2003年）～ 67.5 ℃（2010年），上升速率為0.939 ℃/10 a。最大值為2019年7月27日的76.1 ℃，最小值為2003年7月30日的66.1 ℃。年極端最高值19站次出現(xiàn)在運(yùn)城、臨汾地區(qū)。極端最低地溫平均值為-27.1（2002年）～-16.9 ℃（2017年），上升速率為1.298 ℃/10 a。最大值為2020年12月30日的-26.2 ℃，最小值為1998年1月18日的-40.3 ℃，年極端最低值多出現(xiàn)在大同、朔州、忻州一帶。

山西省近40 a冬季均值為-3.7 ℃，春季均值為14.8 ℃，夏季均值為26.6 ℃，秋季均值為11.3 ℃。各季均值隨年代際變化均呈持續(xù)上升趨勢(shì)，冬、夏季上升趨勢(shì)一致，春、秋季變化幅度相當(dāng)。

2.2 空間分布特征

1981—2020年山西省地表溫度平均值分布及趨勢(shì)系數(shù)與時(shí)間系數(shù)（陰影）疊加如圖2所示。由圖2（a）可知，由北向南均值在7.1～16.6 ℃之間，高緯地區(qū)低于低緯地區(qū)，高海拔地區(qū)低于低海拔地區(qū)。呂梁山脈走向?qū)Φ販胤植加绊懨黠@，西北部低于東南部。極端最高地溫均值在58.1～67.0 ℃之間，最小值位于交口，最大值位于永濟(jì)。極端最低在-31.2～13.3 ℃之間，最小值位于右玉，最大值位于平陸。由圖2（b）可知，各地上升速率在0.185～0.914 ℃/10 a，西北部速率相對(duì)偏快，臨汾西部山區(qū)及太行山一帶升溫較慢。地溫均值與時(shí)間的相關(guān)系數(shù)在0.33～0.90之間，時(shí)間相關(guān)系數(shù)高值區(qū)與趨勢(shì)系數(shù)高值區(qū)基本重合。

2.3 小波分析

各季節(jié)小波系數(shù)實(shí)部變換分布如圖3所示。由圖3可知，山西省近40 a準(zhǔn)18 a周期震蕩非常顯著。2019年后等值線未閉合，預(yù)示未來(lái)地溫有望回落，同時(shí)包含有12 a、6 a的小周期震蕩。冬季呈明顯的12 a周期變化且地溫有望下降。春季準(zhǔn)18 a周期震蕩與6 a周期并存。夏季呈現(xiàn)24 a的高—低周期震蕩，2008年后又反彈上升，提示未來(lái)有繼續(xù)上升傾向。秋季呈現(xiàn)6 a、12 a、18 a周期變化震蕩，2018年后等值線未閉合且信號(hào)增強(qiáng)，預(yù)示秋季地溫更加偏低。

2.4 突變檢驗(yàn)

各季Mann-Kendall統(tǒng)計(jì)量曲線如圖4所示。由圖4可知，全部UF、UB交點(diǎn)均在信度線之外，無(wú)法判斷突變點(diǎn)的存在。但從1987年開(kāi)始UF>0，并在1999年超過(guò)了顯著性水平線，即1999年后升溫加劇。冬季突變點(diǎn)在1997年；春季在2004年；夏季UF、UB在信度線內(nèi)交于1997和2004年，結(jié)合夏季小波變換分析，可認(rèn)為夏季突變點(diǎn)出現(xiàn)在1997年；秋季突變點(diǎn)在2009年。

2.5 經(jīng)驗(yàn)正交函數(shù)分解

對(duì)地溫標(biāo)準(zhǔn)化資料進(jìn)行經(jīng)驗(yàn)正交函數(shù)分解（EOF），第1模態(tài)方差貢獻(xiàn)率達(dá)79.4%，前3模態(tài)累積方差貢獻(xiàn)率為86.6%。

第一特征向量相對(duì)空間函數(shù)均為正值，屬同位相，變化范圍0.07～0.1之間，方差貢獻(xiàn)率為79.4%。即同位相變化占主導(dǎo)地位，反映山西省地溫在大多數(shù)年份整體要么偏高、要么偏低。第一特征向量所對(duì)應(yīng)的時(shí)間系數(shù)正值2019年，負(fù)值2021年，正相位與負(fù)相位事件相當(dāng)，即山西省地溫一致偏高或偏低的概率相當(dāng)，基本反映了近40 a地溫年際變化特征。第二特征向量大致表現(xiàn)為山西省南部大部、西北部與山西省中北部的反相位，變化范圍-0.19～0.25之間，方差貢獻(xiàn)率3.8%。第三特征向量大致表現(xiàn)為山西省東南大部與山西省西北大部的反相位，變化范圍-0.24～0.23之間，方差貢獻(xiàn)率3.4%，與第二特征向量貢獻(xiàn)率相當(dāng)。

2.6 與氣溫、降水量、日照時(shí)數(shù)的相關(guān)分析

地溫與氣溫相關(guān)系數(shù)在0.59～0.96之間，在0.01水平上均呈顯著正相關(guān)；與降水量相關(guān)系數(shù)在-0.47～0.25之間，約2/3站點(diǎn)與降水量呈負(fù)相關(guān)；與日照時(shí)數(shù)相關(guān)系數(shù)在-0.62～0.48之間，僅有32站在0.1水平上顯著相關(guān)，其中顯著正相關(guān)15站，顯著負(fù)相關(guān)17站。

3 結(jié)論

山西省近40 a地溫介于10.7～13.6 ℃之間，上升速率為0.499 ℃/10 a。極端最高、最低地溫也呈上升趨勢(shì)，氣候傾向率分別為0.939 ℃/10 a、1.298 ℃/10 a。地溫分布由南向北次第減小，各地均值在7.1～16.6 ℃之間，高緯地區(qū)明顯低于低緯地區(qū)，高海拔地區(qū)低于低海拔地區(qū)。呂梁山脈走向?qū)Φ販胤植加绊懨黠@，西北部升溫速率快于東南部。

近40 a地溫存在3類尺度周期變化規(guī)律。準(zhǔn)18 a周期震蕩非常顯著，12 a、6 a左右的小周期震蕩包含在大震蕩周期中。山西省1999年之后上升趨勢(shì)更加明顯，冬、春、夏、秋季突變點(diǎn)分別發(fā)生在1997年、2004年、1997年、2009年。地溫同位相變化占主導(dǎo)地位，最大概率是趨于整體一致的偏高或偏低變化。地溫變化與氣溫變化密切相關(guān)。約2/3站點(diǎn)與降水量呈負(fù)相關(guān)，約1/3站點(diǎn)與日照時(shí)數(shù)顯著相關(guān)。

參考文獻(xiàn)：

［1］王佳琳，潘志華，韓國(guó)琳，等.1961—2010年中國(guó)0 cm地溫變化特征及其與氣溫變化的關(guān)系［J］.資源科學(xué)，2016，38（9）：1733-1741.

［2］劉強(qiáng)軍，梁赟，宋軍芳，等.近42年山西省晉城市地溫與氣溫時(shí)空特征研究［J］.江蘇農(nóng)業(yè)科學(xué)，2020，48（6）：240-245.

［3］吉春容，穆彩蕓，胡啟瑞，等.干旱區(qū)哈密瓜生長(zhǎng)和產(chǎn)量與氣候要素的關(guān)系［J］.生態(tài)環(huán)境學(xué)報(bào)，2021，30（1）：10-18.

［4］趙美亮，曹廣超，曹生奎，等.1980—2017年青海省地表溫度時(shí)空變化特征［J］.干旱區(qū)研究，2021，38（1）：178-187.

［5］蔡霞，李亞軍，李智才，等.山西省春季地溫的特征對(duì)崩塌滑坡的影響研究［J］.干旱區(qū)資源與環(huán)境，2022，36（9）：158-165.

［6］張翠華，張文煜，張艷品.輸電線路雷擊前近距離氣象要素變化特征分析：以石家莊市為例［J］.氣象與減災(zāi)研究，2019，42（4）：295-300.

收稿日期：2023-10-17

基金項(xiàng)目：山西省氣象局面上項(xiàng)目（SXKMSQH20226328）。

作者簡(jiǎn)介：范秀平（1974—），女，本科，高級(jí)工程師，研究方向：氣象數(shù)據(jù)質(zhì)控評(píng)估變化分析。