溧阳近55年来气温变化特征分析 溧阳近55年来气温变化特征分析 （南京信息工程大学应用气象学院，南京210044） 摘要：基于溧阳基本站1953—2022年气象数据，分析了该站点1956—1998年气温的月平均、季节平均、极值；1953-2022年任意一年1月、7月气温的变率；1953-2022年任意1年稳定通过0℃、5℃、10℃、15℃和20℃等天数及起止日期结果表明，每年1月平均气温较低，7、8月平均气温较高；每年夏季平均气温最高，冬季平均气温最低；每年温度最大值约在31—32℃，最小值约在-2—-4℃；溧阳稳定通过0℃、5℃、10℃、15℃和20℃的初日呈提前趋势，终日呈推迟趋势，初终间持续时间呈增加趋势 关键词：平均气温界限温度初终日气温变化 引言 近年来，地球气候正经历着一场以变暖为主要特征的显著变化气候变化已经对许多自然系统造成了影响全球气候变化不仅影响自然系统和人类生存环境,也将影响世界经济发展和社会进步[1]因此,气候变化已经不再仅仅是一个科学问题,而日渐成为人们共同关心的重大社会问题[2] 农业界限温度是标志某些重要物候现象或农事活动的开始、终止或转折点的日平均温度,对确定地区布局、耕作制度、结构调整等都具有十分重要的意义[3-4]。
溧阳地处苏浙皖三省交界的群山竹海之畔，山区海拔400—800m，属中纬度亚热带丘陵山区气候，气候温和、四季分明、光照充足、雨量充沛 对溧阳的气候变化，已有一些研究成果，如姜瑞等[5]对溧阳白茶种植的气候适应性分析鉴于此，本文通过分析溧阳气象站的主要气象要素资料，揭示溧阳近55年来气温的变化特征 1.资料与方法 本文选取的气象观测资料为溧阳基本站（58345）1953—2022年的月平均、季节平均、极值；任意一年1月、7月气温的变率；任意1年稳定通过0℃、5℃、10℃、15℃和20℃等天数及起止日期 在Excel、Matlab中进行数据处理及作图,采用相关方法进行线性拟合 2.气温变化分析 2.1 １９５６－１９９８年月平均气温变化 （1）程序： ata=z(37:552,:) save data load data tem=data(:,4:34)/10; season=zeros(172,3); aver_mon=zeros(516,1); for i=1:516 aaa=find(tem(i,:)==3274.4); aver_mon(i,1)=(sum(tem(i,:))-length(aaa)*3274.4)/(31-length(aaa)); end （2）结果： （3）结果分析：从表格可以看出，每年1月平均气温较低，最低温度约在3℃左右，其中1977年1月平均气温为-1.0℃;7、8月平均气温较高,最高温度在28℃左右；并且每年的月平均温度起伏不大，上升趋势不明显。
2.2 １９５６－１９９８年季节平均气温变化 （1）程序： for i=1:172 season(i,1:3)=aver_mon((3*i-2):(3*i),1)'; end aver_sea=mean(season,2); （3）结果分析：从表格可以看出，每年夏季平均气温最高，冬季平均气温最低由于溧阳雨量充沛，对温度有一定的调节作用，所以夏季温度和周围城市相比比较正常，冬季温度也不至于太低其中冬季增暖幅度最大，但夏季与之相反呈变凉趋势 2.3 １９５６－１９９８年温度极值 （1）程序： yearmax=zeros(43,1); yearmin=zeros(43,1); [x,y]=find(tem==3274.4); tem(x,y)=NaN; for i=1:43 yearmax(i)=max(max(tem(12*(i-1)+1:12*(i-1)+12,:))); yearmin(i)=min(min(tem(12*(i-1)+1:12*(i-1)+12,:))); end （3）结果分析：从表格可以看出，每年温度最大值约在31—32℃，出现在夏季，最小值约在-2—-4℃，出现在冬季；在1988年和1992年温度最高，最大值同为33.6℃，在1977年温度最低，最小值为-8.9℃；历年极端温度的起伏不大，在1℃左右摆动。
２．４任意一年1月、7月气温的变率 ２．４．1 １９５６年1月、7月气温变率 绝对变率： 相对变率： （１月） l=length(find(tem(1,:)==NaN)); onesum= (31-l)*aver_mon(1); sum=0; for i=1:31 if(tem(1,i)~=NaN) sum=sum+abs(tem(1,i)-aver_mon(1)); i=1+i; end end oneabs=sum; val_1=oneabs/(31-l); vrl_1=oneabs/onesum*100; （７月） l=length(find(tem(7,:)==NaN)); onesum= (31-l)*aver_mon(7); sum=0; for i=1:31 if(tem(7,i)~=NaN) sum=sum+abs(tem(7,i)-aver_mon(7)); i=1+i; end end oneabs=sum; val_7=oneabs/(31-l); vrl_7=oneabs/onesum*100; （2）结果： 1月绝对变率：2.88 1月相对变率：117.16 7月绝对变率：1.32 7月相对变率：4.40 ２．５稳定通过界限温度的天数及起止日期 ２．５．1 １９５６年稳定通过0℃、5℃、10℃、15℃和20℃等天数及起止日期。
采用五日滑动平均法：由春(秋)季开始出现高于(低于)某界限温度之日起，在依次算出的每个连续5日平均气温中，挑选出第1个其后不再出现低于(高于)该界限温度的五天，并在这五天中挑选出第1个(最后1个)≥该界限温度的日期，作为稳定通过该界限温度的初(终)日；（1）程序： info=zeros(10,373); tmp=[tem(1,1:31)]; tmp=[tmp tem(2,1:29)]; tmp=[tmp tem(3,1:31)]; tmp=[tmp tem(4,1:30)]; tmp=[tmp tem(5,1:31)]; tmp=[tmp tem(6,1:30)]; tmp=[tmp tem(7,1:31)]; tmp=[tmp tem(8,1:31)]; tmp=[tmp tem(9,1:30)]; tmp=[tmp tem(10,1:31)]; tmp=[tmp tem(11,1:30)]; tmp=[tmp tem(12,1:31)]; for i=2:365 if tmp(i)==NaN if(tmp(i+1)==NaN) tmp(i)=tmp(i-1); end tmp(i)=tmp(i-1)+tmp(i+1); end end tmp2=smooth(tmp,5,'moving');%5′????ˉ???ù for k=1:5 aa5=find(tmp2>(k*5-5)); info(2*k-1,1)=length(aa5); info(2*k-1,2)=aa5(1); j=2; for i=2:length(aa5) if((aa5(i)-aa5(i-1))>1) info(2*k,j)=aa5(i); info(2*k-1,j+1)=aa5(i+1); j=j+1; end end info(2*k,j)=aa5(i); end （2）结果： ２．５．2 溧阳站稳定通过10、15、20、℃初日年代际变化表1显示，溧阳稳定通过10℃、15℃、20℃初日均在 21世纪以来有显著地提早，稳定通过10℃初日自 21 世纪显著提前12d; 稳定通过15℃初日在1953-1979 年有提早趋势，但在80年代又推后了几日，在90年后又进入提早阶段; 稳定通过20℃初日在70、80年代为最迟，比50、60年代和之后的90年代迟了5-6d，到21世纪以来则有了显著提早。
表1 1953-2022年溧阳站稳定通过10、15、20、℃初日年代际变化 时间段稳定通过 10℃初日（月/日） 稳定通过 15℃初日（月/日） 稳定通过 20℃初日（月/日） 50年代04-03 04-03 03-30 03-31 03-30 03-18 05-04 04-30 04-22 04-27 04-20 04-18 05-20 05-20 05-26 05-25 05-20 05-12 60年代70年代80年代90年代近几年 数据显示，10℃、15℃、20℃初日均有提前趋势，其中以稳定通过15℃的初日变化最大综上所述，溧阳站稳定通过10℃、15℃、20℃初日在90年代后均已进入提前阶段 ２．６作出1953-2022年年气温变化曲线及距平曲线图２．６.1 1953-2022年年气温变化曲线 （1）程序： M=660;N=34; m=1; for i=1:12:M k=0; a1(m)=0; for n=i:(i+11) for j=1:N if a(n,j)<=500 k=k+1; a1(m)=a1(m)+a(n,j); end end end a2(m)=a1(m)/(k*10); m=m+1; end a3=a2'; for i=1:55 a3(i,2)=1952+i; end plot(a3(:,2),a3(:,1),'b-o') ti。