• CN 62-1175/P
  • ISSN 1006-7639
  • 双月刊
  • 中国科技核心期刊
  • 中国学术期刊综合评价数据库统计源期刊
  • 中文科技期刊数据库收录期刊

干旱气象, 2024, 42(3): 367-375 DOI: 10.11755/j.issn.1006-7639(2024)-03-0367

论文

湖南省地表高温遥感评估指标构建和特征分析

韩沁哲,1,2, 刘海磊3, 范嘉智2,4, 吴浩,2,5, 陈磊士1,2, 欧小锋1,2, 韩沁真6

1.湖南省气象科学研究所,湖南 长沙 410118

2.气象防灾减灾湖南省重点实验室,湖南 长沙 410118

3.成都信息工程大学,四川 成都 610225

4.中国气象局气象干部培训学院湖南分院,湖南 长沙 410125

5.湖南省气候中心,湖南 长沙 410118

6.湖南省长沙市气象局,湖南 长沙 410205

Surface high temperature remote sensing evaluation index construction and characteristic analysis in Hunan Province

HAN Qinzhe,1,2, LIU Hailei3, FAN Jiazhi2,4, WU Hao,2,5, CHEN Leishi1,2, OU Xiaofeng1,2, HAN Qinzhen6

1. Hunan Meteorological Research Institute, Changsha 410118, China

2. Hunan Key Laboratory of Meteorological Disaster Prevention and Reduction, Changsha 410118, China

3. Chengdu University of Information Technology, Chengdu 610225, China

4. China Meteorological Administration Training Centre Hunan Branch, Changsha 410125, China

5. Hunan Climate Center, Changsha 410118, China

6. Changsha Meteorological Bureau of Hunan Province, Changsha 410205, China

通讯作者: 吴浩(1988—),男,硕士,高级工程师,主要从事气候变化研究。E-mail:wuhaophy@163.com

责任编辑: 王涓力;校对:邓祖琴

收稿日期: 2023-06-9   修回日期: 2023-12-12  

基金资助: 湖南省自然科学联合基金项目(2023JJ60542)
湖南省气象局重点项目(XQKJ18A003)
国家自然科学基金项目(42005058)
中国气象局气象能力提升联合研究专项重点项目(22NLTSZ003)

Received: 2023-06-9   Revised: 2023-12-12  

作者简介 About authors

韩沁哲(1980—),女,湖南长沙人,硕士,高级工程师,主要研究方向为生态遥感和防灾减灾。E-mail:6665599@qq.com

摘要

基于长时序遥感产品构建地表高温评估指标,分析湖南省近20 a地表高温分布特征,可为农业生产和防灾减灾提供技术参考。首先利用2002—2021年的MYD11A1、MYD13A1、MYD09GA长时序数据回归计算日地表高温数据集,然后基于百分比法、百分位法和趋势比率法构建地表高温遥感评估指标,并分析其时空分布特征。结果表明,使用日地表温度大于等于45 ℃面积占比作为日尺度的地表高温等级阈值划分指标与气象预警指标一致性最好,湖南省地表高温年最大值以2013年为转折点增加显著;基于百分位法确定的湖南省年尺度的地表高温极端阈值(40.2~64.1 ℃),可作为高温强度的评价指标;利用湖南省地表高温趋势率来评价湖南省地表高温年增减趋势具有较好的应用价值,长株潭城市群和洞庭湖区的快速增温、衡邵干旱走廊的极端高温需要重点关注。

关键词: 地表高温; 评估指标; 高温等级; 卫星遥感

Abstract

Based on long-term remote sensing products, this paper constructs evaluation indicators for surface high temperatures and analyzes the distribution characteristics of surface high temperatures in Hunan Province over the past 20 years, aiming to provide technical references for agricultural production, disaster prevention and mitigation. Firstly, the daily surface high temperature dataset are calculated using long-term data regression of MYD11A1, MYD13A1 and MYD09GA from 2002 to 2021. The percentage method, percentile method and trend ratio method are used to construct remote sensing evaluation indicators for surface high temperatures, then the spatial and temporal distribution characteristics of surface high temperatures are analyzed. The results indicate that when using the proportion of areas with daily surface temperatures greater than or equal to 45 ℃ as the threshold for classifying surface high temperature levels on a daily scale, it aligns best with meteorological warning indicators. The annual maximum surface high temperature in Hunan Province showed a significant increase with 2013 as a turning point. The extreme threshold for annual surface high temperatures in Hunan Province (ranging from 40.2 to 64.1 ℃), determined using the percentile method, can serve as an evaluation metric for high temperature intensity. Additionally, using the trend rate of surface high temperatures in Hunan Province to evaluate the increase or decrease trends of annual surface high temperatures has good application value. Rapid warming in the Changsha-Zhuzhou-Xiangtan urban agglomeration and the Dongting Lake area, as well as extreme high temperatures in the Hengshao drought corridor, require special attention in the future.

Keywords: surface high temperature; evaluation indicators; high temperature level; satellite remote sensing

PDF (13410KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

韩沁哲, 刘海磊, 范嘉智, 吴浩, 陈磊士, 欧小锋, 韩沁真. 湖南省地表高温遥感评估指标构建和特征分析[J]. 干旱气象, 2024, 42(3): 367-375 DOI:10.11755/j.issn.1006-7639(2024)-03-0367

HAN Qinzhe, LIU Hailei, FAN Jiazhi, WU Hao, CHEN Leishi, OU Xiaofeng, HAN Qinzhen. Surface high temperature remote sensing evaluation index construction and characteristic analysis in Hunan Province[J]. Arid Meteorology, 2024, 42(3): 367-375 DOI:10.11755/j.issn.1006-7639(2024)-03-0367

0 引言

20世纪90年代以来,在全球气候变暖背景下,极端高温事件致灾强度大、发生频次高、影响范围广(封国林等,2012;Sun et al.,2014;董思言等,2014),对农业生产、城市建设及社会经济发展造成极大的不确定性(刘琳和徐宗学,2014)。湖南是盛夏持续高温干旱发生频率较高的地区之一(罗伯良等,2010),且自2000年以后,湖南高温热害发生次数明显增多(张超等,2012),2003、2009、2011、2013、2022年是湖南典型高温年份。其中,2013年6月29日至8月19日,湖南省84县市最高气温超过39.0 ℃,57县市最高气温超过40.0 ℃,高温过程造成农作物受灾面积207.6万hm2,绝收面积42.5万hm2,直接经济损失170.2亿元人民币(罗伯良和李易芝,2014);2022年7—11月出现持续近5个月的历史罕见夏秋高温天气,平均气温创1961年以来历史新高,同时期的降水异常偏少造成湖南历史最强夏秋连旱,全省流域面积50、100 km2以上河流断流的分别有88、49条(①数据来源于《2022年湖南气候影响评价》。),灌溉条件差的水稻减产或绝收,油茶、柑橘、葡萄落果严重或植株干枯死亡,油菜无法播种或移栽(②数据来源于《湖南省2022年9月(10月、11月)农业气象月报》。)。构建区域高温遥感评估指标可从卫星遥感监测角度揭示湖南省地表极端高温分布规律,推动高温遥感监测规范化,对于区域农业生产结构调整、能源水资源储备和森林防火调度等都具有重要指导意义(邹德全等,2021)。

国内外学者对极端高温监测评估方法开展了一系列研究。如侯威等(2012)基于随机重排去趋势波动分析法寻找收敛点确定极端高温阈值;陈敏等(2013)利用年、月累计高温有效积温和中位数、四分位数统计高温变化特征;王钊等(2016)使用极端气候指标计算软件RClimdex计算5种极端温度指数分析秦岭南北极端温度变化特征及其与区域增暖的关系;程肖侠等(2020)建立了单站高温过程综合强度指标;纪玲玲等(2020)采用气候倾向率、距平累积、灰色关联度等方法建立了吉林省高温过程综合评估指数;林爱兰等(2021)利用区域站点夏季日最高气温第80百分位值来评估区域性的高温强度;邓甜甜等(2023)用趋势分析方法结合站点数据来分析渭河流域高温热浪时空变化特征。上述极端高温评估方法的研究,多基于台站观测数据,但传统的地面气象观测站点密度低且离散(陈岳等,2022),只能采用空间插值法实现网格温度的估算(马诺等,2020),或引入气象数值模式再分析产品,但因为土地利用资料更新缓慢(温晓培等,2022),数值模式仍会造成城市和地形复杂区域的模拟偏差(董春卿等,2021)。遥感是唯一可以提供二维栅格陆面温度信息的手段,可以快速获取大面积区域地表温度(朱怀松等,2007),因此,利用遥感数据进行地表高温遥感评估相对更为客观。

本文以湖南省为研究区域,基于2002—2021年空间分辨率为1 km的MYD11A1长时序卫星遥感数据,以及空间分辨率为500 m的归一化植被指数16 d反演产品MYD13A1、逐日地表反射率产品MYD09GA,参考国外学者基于遥感时序数据的日最高气温估算方法(Sheng et al.,2012;Georgy et al.,2014;Elizabeth,2015),通过多元线性回归方程实现空间分辨率百米级的逐日地表高温监测,并进一步利用百分比法、百分位法和趋势比率法等气候统计方法构建湖南省地表高温卫星遥感监测评估指标,可作为地表高温遥感监测定量定性评估的依据,进而研究极端高温的区域分布状况,以期提升极端高温事件的卫星遥感监测业务能力。

1 数据与方法

1.1 数据来源

指标建立使用美国地质勘探局提取和分析样本的应用工具(Application for Extracting and Exploring Analysis Ready Samples,AppEEARS)(https://appeears.earthdatacloud.nasa.gov/)下载的4组经美国国家航空航天局(National Aeronautics and Space Administration,NASA)官方发布的数据:湖南省2002—2021年MODIS(Moderate-resolution Imaging Spectroradiometer)地表温度(Land Surface Temperature,LST)逐日反演产品MYD11A1、归一化植被指数(Normalized Difference Vegetation Index,NDVI)16 d反演产品MYD13A1,以及逐日地表反射率产品MYD09GA(Tang et al.,2022)。湖南省2002—2021年97个国家级气象站的逐小时气温、日最高气温和裸土的地表温度、日最高地温数据来自国家气象信息中心。

辅助数据包括:“地理空间数据云”(https://www.gscloud.cn/)公布的ASTER GDEM v2 30 m空间分辨率数字高程(李振林和王晶,2013);中国科学院资源环境科学与数据中心(https://www.resdc.cn/)公布的中国1:10万土地利用现状遥感监测数据(刘纪远等,2020)和中国1:100万地貌类型空间分布数据(周成虎等,2009;程维明和周成虎,2014);美国国家物理数据中心(https://www.ngdc.noaa.gov/dmsp)发布的2017年稳定夜间灯光数据集(白彩全等,2017)。

1.2 研究区概况

湖南省地处云贵高原向江南丘陵和南岭山脉过渡的地带,呈东南西三面环山、东北面开口的马蹄形地貌,从西到东分别为武陵雪峰山、南岭和幕阜山等中低海拔山脉,气温垂直变化较大;东北面开口处为地势最低的洞庭湖平原,为湘、资、沅、澧四水汇聚地,是长江流域重要的调蓄湖泊(图1)。从洞庭湖顺着省内最大的河流湘江南下,有湖南省经济发展的核心区长株潭城市群,是近年来省内城镇化建设进程最快、规模最大的地区;湖南省中部的衡邵盆地,分布着面积30.38万 hm2的岩溶地貌,形成风化剥蚀台地,岩面温度高于其他地貌的地表温度(莫跃爽和周秋文,2019),且易发生春夏之交的迎风坡地形雨和盛夏初秋的背风坡焚风增温效应。

图1

图1   湖南省地形和区域概况

Fig.1   The distribution of terrain and regional overview of Hunan Province


1.3 计算地表反射率

参考Jones等(2004)利用LST产品与夜间最低气温建立回归方程的方法,以及Kawashima等(2000)用LST、NDVI与气温建立二元回归模型减少误差的方法,徐伟燕等(2017)考虑了LST、NDVI、地表反照率、经度、纬度、高程6个影响因子建立回归方程估算最高气温,发现最高温模型中地表温度、地表反照率、NDVI的偏相关系数分别为0.97、-0.37、0.45。考虑到地表反照率越大,地表吸收的短波辐射越小,地表温度也越低,因此本文引入地表反射率,与LST、NDVI共同建立日最高气温和最高地表温度的多元线性回归方程。其中,地表反射率计算公式(Liang,2001)如下:

α=0.16α1+0.291α2+0.243α3+0.116α4+0.112α5+0.081α7-0.001 5

式中:α为地表反射率,α1~α5α7为MYD09GA数据第1~5以及第7波段地表反射率。

1.4 估算日最高温度

首先把LST、NDVI和地表反射率3种卫星产品处理成长时序数据矩阵,空间匹配参考曹广真等(2010)的重采样方法,使所有卫星产品的时空分辨率一致(表1),然后去除缺测值和异常值,得到符合正态分布特征的实验样本45 015个,随机挑选70%样本做回归,30%样本做检验。日最高气温和最高地表温度的多元线性回归方程如下:

Tairmax=0.9×LST-1.11×NDVI-1.84×α+5.1
Tgroundmax=1.3×LST+7.0×NDVI-2.9×α+10.3

式中:Tairmax(℃)为日最高气温;Tgroundmax(℃)为日最高地表温度。

表1   数据预处理方法

Tab.1  Data preprocessing methods

数据类型数据来源预处理过程
LSTMYD11A1通过质量控制字段把云区和异常值赋为NaN空值,由于MYD11A1的1×1格对应的MYD13A1和MYD09GA的2×2格,因此通过1格复制为4格的方法,得到同空间尺度数据
NDVIMYD13A1MYD13A1 时间分辨率为16 d,因此一个数据对应当天至16 d内的MYD11A1和MYD09GA数据
地表反射率MYD09GA利用Liang的方法(Liang,2001)计算地表反射率,时间分辨率为1 d

新窗口打开| 下载CSV


1.5 效果对比

表2可见,MYD11A1反演产品卫星过境时间与同时次气象站观测的地面温度(0 cm)数据均方根误差(Root Mean Square Error,RMSE)较大。主要是由于遥感反演的地表温度是某一像元的平均地表温度,可能包含多个下垫面类型,而同一时空条件下不同下垫面类型的地表发射率不同,其向外辐射能量不同(余卫国等,2019);而气象站逐小时观测的地面温度(0 cm)是一个点上裸土表层的真实分子动力学温度;因此遥感监测到的地表温度易受地表反射率和植被冠层的影响,造成与站点观测到的裸土地表温度偏差较大。MYD11A1与气象站同时次观测到的气温、地温形成的长时序数据序列相关性均较好,说明温度变化趋势一致(张富荣等,2015;刘强军等,2020)。

表2   MYD11A1及推导值与国家气象站观测值误差和相关性对比

Tab.2  Comparison of error and correlation between MYD11A1 and observations of national meteorological stations

产品类型RRMSE/℃RRMSE/℃
过境同时次
气温
过境同时次
地温
过境同时次
气温
过境同时次
地温
日最高气温日最高地表温度日最高气温日最高地表温度
MYD11A10.950.934.333.430.950.9313.8716.37
回归方程计算结果0.940.941.992.51

新窗口打开| 下载CSV


计算MYD11A1及多元线性回归方程推导出的日最高气温、日最高地表温度与气象站当日观测到的最高气温和最高地表温度的RMSE,发现其值分别从13.87、16.37 ℃减小到1.99、2.51 ℃。这主要是LST产品监测到的是瞬时值,该时次不一定是当日最高温度出现时间,通过多元线性回归方程将植被指数和地表反射率带入推算最高温度,可有效降低误差。

1.6 地表高温趋势率

参考管延龙等(2015)通过线性回归分析方法对最高、最低地温的线性倾向率进行分析,利用趋势比率法计算地表高温趋势率,评估地表高温强度的年增减趋势。趋势率负值为下降趋势,正值为增加趋势。

Ci=a0+a1ti

式中:Ci (℃)为第i年年最高地表温度;a0为常数;a1(℃·a-1)为年地表高温趋势率;ti为第i年减去基础年份得到的年份数。

2 指标建立与结果分析

2.1 日尺度地表高温评估指标建立

基于百分比法,建立日尺度的地表高温评估指数:

D=NTmaxN×100%

式中:D(%)为像元数百分比,是地表高温等级评估指数;NTmax为最高地表温度大于某个阈值的像元数,根据长时序日最高地表温度数据计算;N为有效值区域(非云、非异常值、非空值)的像元数。

依据湖南省气象灾害预警技术标准(表3),利用2002—2021年的地面站观测的逐日日最高气温数据,分等级进行全省高温预警日期统计,计算出湖南省近20 a的所有高温预警事件。

表3   湖南省高温预警指标

Tab.3  High air temperature warning index in Hunan Province

新窗口打开| 下载CSV


参考郑嘉雯等(2019)在集合预报极端低温预警信号阈值确定中所使用的TS评分法,测试地表高温评估指数D在40~50 ℃不同温度阈值、5%~95%不同地表高温面积百分比,与高温预警等级的判识结果进行对比,指标适用性定量评分标准公式如下:

TS评分=正确判识个数/(正确判识个数+误判个数)
样本数归一化指数=(当前样本数-最少样本数)/(最多样本数-最少样本数)
正确率评分=TS评分×样本数归一化指数

式中:正确判识个数为该条件下正确判识的总数;误判个数为地表高温指标判识为高温,但当日却没有高温预警的个数;当前样本数为满足该条件下当前判识出的样本数;最少样本数为5%~95%不同面积占比条件下判识出的最少样本数;最多样本数指5%~95%不同面积占比条件下判识出的最多样本数。云干扰导致的数据缺失过多使漏判率不具有代表性,因此在TS评分计算中不考虑。

通过正确率评分指标计算发现,黄色预警在湖南省内地表温度≥45 ℃面积占比大于35%时匹配性最好,正确率评分最高[图2(a)],而湖南省内地表温度≥45 ℃面积占比大于40%后评分下降,原因是可判识出的高温样本数减少,这样就可以保证正确判识率的情况下,尽量多的筛选出高温样本。以此类推计算出橙色预警在湖南省内地表温度≥45 ℃面积占比大于55%时匹配性最好,正确率评分最高;红色预警在湖南省内地表温度≥45 ℃面积占比大于80%时匹配性最好,正确率评分最高。通过图2(a)正确率评分的比较,得到日尺度的地表高温等级划分指标(表4)。

图2

图2   湖南省内地表温度≥45 ℃面积占比的高温预警正确率评分(a)和2002—2021年地表高温发生次数和最大值年际变化(b)

Fig.2   Accuracy score of high temperature early warning for the area proportion with surface temperature ≥45 ℃ in Hunan Province (a) and interannual variation of frequency and maximum value of high LST from 2002 to 2021 (b)


表4   日尺度地表高温等级划分指标

Tab.4  Daily high LST classification index

高温等级像元数百分比
黄色[35%,55%)
橙色[55%,80%)
红色≥80%

新窗口打开| 下载CSV


利用日尺度的卫星遥感地表高温等级划分指标,统计湖南省2002—2021年地表高温发生次数,可以得到地表高温发生频率的年变化曲线[图2(b)]。可见,自2013年来高温发生频率有所增加,且2013年后每隔3 a出现一次高温多发年份,这与贾艳青等(2017)研究发现的长江三角洲地区极端气温存在3~5 a振荡周期,以及黄卓禹等(2016)通过小波方差图计算出的湖南省5 a特征时间尺度上高温日数变化平均周期为3.5 a的结论一致。地表高温发生次数较多的年份与高温站次出现较多的年份基本一致(2003、2011、2013年),年地表高温最大值较高的年份与夏季高温距平偏高的年份基本一致(2003、2009、2011、2013年)(戴泽军等,2014)。2013年为该气候期的高温次数极多年,且从2012—2015年呈陡升陡降的变化特征。

2004—2012年地表高温年最大值为59.0~63.6 ℃;2013年为转折点,出现近20 a年最大值(66.8 ℃),2014年出现近20 a年最小值(58.7 ℃),2015—2021年年最大值在62.8~66.6 ℃内波动,与2004—2012年比较,平均值增加2.2 ℃。

2.2 年尺度地表高温评估指标建立

21世纪初,政府间气候变化专门委员会(Intergovernmental Panel on Climate Change,IPCC)第三与第四次评估报告对极端天气事件给出明确定义:极端天气事件是指发生概率小于观测记录概率密度函数第10个百分位数或超过第90个百分位数的天气事件。在气象行业标准《极端高温监测指标》(QX/T 280—2015)(国家气候中心,2015)中,利用百分位法确定年尺度的极端高温阈值,即取气候标准期(连续30 a)内任一高温指标每年的极值和次极值,构建一个包含60个样本的序列;对序列从小到大排序,取第95百分位数为偏大的极端阈值,大于或等于该阈值的事件为极端偏大事件。郑嘉雯等(2019)在极端温度预报中也把第95个百分位及以上的温度值定义为极端高温来计算检验评分。参考这一计算方法,鉴于MODIS长时序遥感产品的累计时间有限,本文以2002—2021年为一个气候期(连续20 a)计算年最高地表温度和年次高地表温度,从小到大排序,取第95百分位数定义年尺度的地表高温极端阈值,作为年尺度的地表高温定量评估依据。

图3为基于年尺度地表高温评估指标计算出的2002—2021年湖南省年地表高温极端阈值分布(地表高温极端阈值为40.2~64.1 ℃)以及基于湖南省1991—2020年一个标准气候期(30 a)的国家气象站观测的年最高气温和年次高气温排序得到的年气温极端高温阈值分布(极端气温阈值为36.0~41.8 ℃),由于站点的地表温度观测数据时间跨度不够,因此使用气温进行计算。对比两者分布特征大体一致,但卫星数据估算的地表高温阈值区域代表性更精细,与侯浩然等(2018)的研究结论一致。

图3

图3   2002—2021年遥感监测年地表高温极端阈值(a)和1991—2020年气象站极端气温阈值(b)分布

Fig.3   The distribution of the extreme thresholds of surface high temperature from remote sensing monitoring from 2002 to 2021 (a) and the extreme air temperature threshold from meteorological stations from 1991 to 2020 (b)


表5可以看出,洞庭湖周边区域、长株潭城市群、衡邵盆地等0~200 m的低海拔地区(图1),年地表高温极端阈值为64.0 ℃左右;200~500 m的中海拔地区,主要为山地丘陵之间的盆地,年地表高温极端阈值为62.8 ℃左右;而海拔500 m以上的山地年地表高温极端阈值随海拔高度增加而下降,但降速并不均匀,平均为海拔高度每上升100 m地表高温极端阈值下降1.0 ℃。

表5   不同海拔地表高温极端阈值范围

Tab.5  Threshold range of extreme high LST at different altitudes

海拔高度/m极端高温阈值/℃海拔高度/m极端高温阈值/℃海拔高度/m极端高温阈值/℃海拔高度/m极端高温阈值/℃
0~10042.6~64.0>500~60044.4~60.7>1 000~1 10040.4~52.1>1 500~1 60040.8~50.7
>100~20046.1~64.1>600~70042.1~58.2>1 100~1 20041.7~52.0>1 600~1 70040.5~50.3
>200~30044.0~62.8>700~80042.9~57.6>1 200~1 30041.1~51.2>1 700~1 80040.4~50.9
>300~40045.3~62.4>800~90042.0~54.8>1 300~1 40040.5~51.4>1 800~1 90040.2~47.5
>400~50045.6~62.4>900~1 00041.9~52.9>1 400~1 50040.4~50.5>1 900~2 00040.6~47.1

新窗口打开| 下载CSV


2.3 年尺度地表高温趋势率

根据长时序日最高地表温度数据,计算2002—2021年湖南省年最高地表温度分布,再进一步逐像元计算地表高温年趋势率(图4),分析趋势率的空间变化特征。可以看出,地表高温增加趋势最为明显的是长株潭城市群(图4),其次为洞庭湖区、衡阳市和永州北部。

图4

图4   2002—2021年湖南省年最高地表温度趋势率(单位:℃·a-1

Fig.4   The trend rate of annual maximum surface temperature in Hunan Province from 2002 to 2021 (Unit: ℃·a-1


城镇和建设用地地表高温上升趋势明显(表6),其次为河湖湿地水田、农村居民点和平原旱地,而林地、草地和山丘旱地趋势持平。黄鹤楼等(2020)研究认为城市区域增温很大部分由城市化进程带来的热岛效应造成,而其他受人类活动影响较小的区域,增温更能体现气候变暖特征。从地貌类型看,洞庭湖区的湖泊、平原增温明显,位于衡邵盆地的剥蚀台地次之,山地趋势持平或略减小(表7)。

表6   不同地表类型的地表高温年趋势率

Tab.6  Annual trend rates of high LST for different land surface types 单位:℃·a-1

类型趋势率类型趋势率类型趋势率类型趋势率
林地0.002 6中覆盖草地0.017 6城镇0.121 8平原水田0.076 4
灌木林0.012 20.055 1农村居民点0.072 9山地旱地0.007 4
疏林地0.022 60.088 6其他建设用地0.128 8丘陵旱地0.032 3
其他林地0.021 3水塘0.064 1山地水田0.014 8平原旱地0.076 4
高覆盖草地0.003 9湿地0.091 7丘陵水田0.039 6

新窗口打开| 下载CSV


表7   不同地貌类型的地表高温年趋势率

Tab.7  Annual trend rates of high LST for different geomorphic types 单位:℃·a-1

类型趋势率类型趋势率类型趋势率
冲积平原0.079 9洪积平原0.097 5湖泊0.096 9
低海拔山地0.019 1剥蚀台地0.036 7中海拔山地-0.006 8

新窗口打开| 下载CSV


3 结论与讨论

本文基于长时序遥感产品开展地表高温时空特征分析,并与预警信号进行比对,以百分位法和趋势比率法,针对湖南省年尺度地表高温极端阈值、地表高温趋势率开展研究。

基于百分比法建立湖南省日尺度的地表高温等级划分指标,并依据《湖南省气象灾害预警技术标准》计算历史高温样本进行正确率评分,发现黄色预警在湖南省内地表温度≥45 ℃面积占比大于35%时匹配性最好,橙色预警在湖南省内地表温度≥45 ℃面积占比大于55%时匹配性最好,正确率评分最高;红色预警在湖南省内地表温度≥45 ℃面积占比大于80%时匹配性最好,可作为日尺度高温评估和高温频次的评价指标。利用该指标统计湖南省2002—2021年高温发生次数,发现自2013年来高温发生频率有所增加,且2013年后每隔3 a出现一次高温多发年,且以2013年为转折点,年地表高温最大值波动区间从59.0~63.6 ℃提升到62.8~66.6 ℃。

基于百分位法确定湖南省年尺度的地表高温极端阈值,可作为高温强度的评价指标。利用该指标计算出湖南省年地表高温极端阈值(40.2~64.1 ℃)分布,与国家气象站观测值计算出的年气温极端高温阈值(36.0~41.8 ℃)分布特征大体一致,但台站区域代表性弱于卫星监测结果。另外,海拔高度每上升100 m,地表高温极端阈值下降1.0 ℃。

基于趋势比率法计算湖南省地表高温趋势率,可作为地表高温年增减趋势的评价指标。从行政区域看,趋势率增加最明显的是长株潭城市群和洞庭湖区、衡阳市和永州北部;从地表类型看,趋势率增加最明显的是城镇和建设用地;从地貌类型看,趋势率增加最明显的是洞庭湖区。

综上所述,从气候分析角度来说,应关注2013年以后的极端高温频发事件,今后应重点针对长株潭城市群和洞庭湖区的快速增温、衡邵剥蚀台地的极端高温进一步开展研究。随着静止气象卫星的发展,分钟级地表温度监测产品已在逐步应用,下一步工作将针对静止气象卫星的极端高温估算开展研究,并探索该高温评估指标的应用效果。

参考文献

白彩全, 易行, 何晨, 2017.

夜间灯光遥感数据应用研究的文献计量分析——以美国国防气象卫星计划运行线扫描系统为例

[J]. 测绘科学, 42(8): 116-123.

[本文引用: 1]

曹广真, 侯鹏, 毛显强, 2010.

北京市城市化对地表温度时空特征的影响

[J]. 气象, 36(3): 19-26.

[本文引用: 1]

陈敏, 耿福海, 马雷鸣, , 2013.

近138年上海地区高温热浪事件分析

[J]. 高原气象, 32(2): 597-607.

DOI      [本文引用: 1]

深入了解特大城市高温热浪的变化特征和发生条件对于科学防御其引发的灾害具有重要的实际意义。根据中国气象局高温日和高温预警信号阈值定义了上海地区高温热浪三级标准, 并引入高温有效积温, 基于1873-2010年上海市徐家汇站日最高气温资料分析了上海地区高温热浪的多尺度时频特征; 基于环流指数和NCEP再分析资料, 探讨了上海地区高温热浪异常偏强年的大气环流异常特征。结果表明: (1)将高温有效积温与高温热浪发生频次结合, 可更合理地表征高温热浪的炎热程度。(2)近138年来上海共发生214次高温热浪事件, 其平均高温有效积温为8.3 ℃, 7月高温热浪发生频次和高温有效积温均多于8月。(3)近138年来上海地区高温热浪有三段持续偏多与偏强期: 19世纪90年代初至19世纪末、 20世纪20年代末至50年代初、 80年代末尤其是21世纪初以来; 最强的高温热浪事件出现在1934年; 而21世纪初的10年, 其炎热程度正呈明显上升趋势。(4)上海地区高温热浪异常偏强年的主要环流特征是夏季北半球副热带高压异常偏强、 西太平洋副热带高压异常偏西和印缅槽异常偏弱, 如高温热浪异常偏强年夏季平均的西太平洋副热带高压西伸脊点伸展至我国东南沿海约122°E附近, 而

陈岳, 寇卫利, 李莹, , 2022.

农作物遥感灾损评估研究进展

[J]. 江苏农业科学, 50(16): 12-20.

[本文引用: 1]

程维明, 周成虎, 2014.

多尺度数字地貌等级分类方法

[J]. 地理科学进展, 33(1): 23-33.

DOI      [本文引用: 1]

参考已出版的全国各级各类比例尺的地貌类型图的分类方案及图例表达,探讨了中国1:100 万数字地貌的等级分类方法,采用基于形态、成因、物质和年龄等地貌要素,综合反映地貌特征的等级分类指标和分类体系,初步构建了中国多个国家基本比例尺(即1:400 万、1:100 万、1:50 万、1:25 万、1:5 万)数字地貌等级分类方法,发展了由连续分布的多边形图斑反映形态成因类型,以及由离散的点、线和面图斑共同反映形态结构类型的数字地貌类型数据组织方式,构建了多尺度数字地貌类型的编码方法。该研究可为发展多尺度地貌类型图的编制提供方法基础,也可为当前正在进行的地理国情监测工程的大比例尺地貌类型信息普查提供分类规范和技术支持。

程肖侠, 蔡新玲, 李明, , 2020.

陕西省高温时空变化特征及其定量化评估方法

[J]. 干旱区研究, 37(4): 956-965.

[本文引用: 1]

戴泽军, 李易芝, 刘志雄, , 2014.

1961—2013年湖南夏季高温气候特征

[J]. 干旱气象, 32(5): 706-711.

DOI      [本文引用: 1]

基于1961~2013年湖南82 个台站夏季逐日最高气温数据,采取线性趋势分析等方法分析了近53 a湖南高温(最高气温35 ℃)气候特征。结果表明,湖南省夏季多年平均高温日数高温均值及极端最高气温的大值区主要分布在湘东偏南及湘中偏北的部分地区。全省区域平均的多年平均高温站次及高温值总体呈线性增加趋势,7月中旬到8月上旬高温站次出现相对较多,7月8日起高温均值达到36 ℃以上并一直维持到8月末高温最大持续日数总体呈现从湘中衡阳株洲等地市向四周递减的分布趋势,极端持续高温日数阈值与高温最大持续日数空间分布较相似,但极端持续高温事件发生次数较多的台站主要分布在湘北湘东南及湘西南。

邓甜甜, 耿广坡, 杨睿, , 2023.

1980—2020年渭河流域高温热浪时空变化特征

[J]. 干旱区地理, 46(2): 211-221.

[本文引用: 1]

董春卿, 郭媛媛, 张磊, , 2021.

基于CLDAS的格点温度预报偏差订正方法

[J]. 干旱气象, 39(5): 847-856.

DOI      [本文引用: 1]

利用国家气象信息中心CLDAS格点温度实况、中央气象台SCMOC格点温度预报以及山西省站点观测温度,采用非独立性检验综合评估CLDAS在山西区域的适用性。在此基础上,采用滑动训练期订正方案,基于格点实况开展SCMOC温度预报场的客观订正。结果表明:(1)复杂地形对山西CLDAS格点温度实况的精度有一定影响,但最高气温的分析精度优于最低气温,表明地形对最低气温的偏差影响更显著,高海拔地区CLDAS最低气温一般对应为负偏差,低海拔地区一般对应为正偏差。(2)CLDAS格点温度实况的偏差空间分布具有时间延续性,进行简单的系统偏差订正后,最高、最低气温格点实况的精度分别提升1.1%、9.7%,与站点观测更为吻合。(3)基于改进后的CLDAS格点温度实况,采用滑动偏差订正方案,显著改善了山西省SCMOC温度预报的准确率。2019年,滑动偏差订正后的24 h时效最高、最低气温预报准确率较SCMOC温度预报分别提升2.7%、4.7%,订正后的短期温度预报质量有较大提高,优于预报员主观预报。

董思言, 徐影, 周波涛, , 2014.

基于CMIP5模式的中国地区未来高温灾害风险预估

[J]. 气候变化研究进展, 10(5): 365-369.

[本文引用: 1]

封国林, 侯威, 支蓉, , 2012. 极端气候事件的检测、诊断与可预测性研究[M]. 北京: 科学出版社.

[本文引用: 1]

管延龙, 王让会, 李成, , 2015.

天山北麓1963—2010年0 cm最高与最低地表温度变化特征

[J]. 干旱气象, 33(4): 587-594.

DOI      [本文引用: 1]

基于天山北麓7个气象站1963~2010年逐月0 cm最高、最低地温资料,采用线性趋势分析、Mann-Kendall检验、Morlet小波等方法,分析了天山北麓地温变化特征。结果表明:(1)近48 a来,天山北麓0 cm最高地温以精河为高值中心,总体上呈西高东低的特征,最低地温以乌苏、乌鲁木齐为高值中心的中西部地区普遍较高;(2)0 cm地温呈显著上升趋势,最低地温增幅尤为显著,达0.87 ℃/10 a;(3)0 cm最低地温在2002年发生突变,而最高地温未出现突变;(4)0 cm地温异常年份主要发生在2006年之后,以偏暖为主。

国家气候中心, 2015. 极端高温监测指标(GX/T280—2015)[S]. 北京: 国家气候中心.

[本文引用: 1]

侯浩然, 丁凤, 黎勤生, 2018.

近20年来福州城市热环境变化遥感分析

[J]. 地球信息科学学报, 20 (3): 385-395.

DOI      [本文引用: 1]

城市热环境是城市微气候的重要组成部分,已成为近年来的研究热点。受制于卫星传感器较低的热红外波段空间分辨率,此类数据反演得到的地表温度难以反映城市热环境的实际情况。为解决这一困境,本文利用空间降尺度HUTS算法反演得到30 m空间分辨率的福州市中心城区1994年5月12日、2003年5月29日和2016年7月27日3个时相的地表温度影像。在此基础上,结合土地利用等数据对热环境的时空变化做定量分析,并进一步引入景观指数,分析近20年间福州市中心城区高温度等级斑块的形态变化。结果表明:① 近20年间随着城市拓展,福州市建成区的高温区域面积从35.75 km<sup>2</sup>增加到184.11 km<sup>2</sup>,高温度等级斑块不断从市中心向四周扩散;② 市中心的特高温斑块和高温斑块趋向破裂、分散,聚集程度下降,次高温斑块的面积与占比均大幅提升,成为建成区内高温区域的主要组成部分;③ 城市热岛比例指数URI由0.39上升到0.52,热岛效应明显加强。总体上,近20 a间福州市建成区的热环境变化较大,其中鼓楼区南部、台江区和晋安区南部的高温区域聚集现象有所改善,而仓山区、马尾区和闽侯县的大部分区域在经历快速城市化过程后温度等级明显升高。

侯威, 章大全, 钱忠华, , 2012.

基于随机重排去趋势波动分析的极端高温事件研究及其综合指标的建立

[J]. 高原气象, 31(2): 329-341.

[本文引用: 1]

基于中国气象局公布的1961-2006年中国165个国际交换站无缺测的逐日平均气温资料, 利用随机重排去趋势波动分析(S-DFA)方法, 计算并分析了中国极端高温事件阈值的空间分布特征, 并对S-DFA方法在实际资料中的应用进行了检验。从可预报性的角度给出了极端高温事件强度综合指标的定义。该综合指标将极端高温事件的发生频次和强度综合起来, 兼顾不同地区各自特有的区域气候背景, 进一步说明综合指标定义的合理性。基于极端高温事件综合指标的空间分布规律, 将1961-2006年间中国极端高温事件分为4个不同等级的地区。极端高温事件综合指标在20世纪90年代初期之前一直保持平稳的波动变化, 之后则一直处于上升之中, 尤其是在90年代中后期开始迅速上升。

黄鹤楼, 丁烨毅, 姚日升, , 2020.

宁波城市化进程对夏季极端气温和高温热浪的影响

[J]. 气象与环境科学, 43(2): 56-62.

[本文引用: 1]

黄卓禹, 王式功, 李晓坤, , 2016.

近35 a湖南夏季高温日数时空变化特征分析

[J]. 贵州气象, 40(1): 35-39.

[本文引用: 1]

纪玲玲, 袭祝香, 陈卓, , 2020.

吉林省高温天气气候特征及其过程等级评估

[J]. 气象与环境学报, 36(3): 49-54.

[本文引用: 1]

贾艳青, 张勃, 张耀宗, , 2017.

长江三角洲地区极端气温事件变化特征及其与ENSO的关系

[J]. 生态学报, 37(19): 6 402-6 414.

[本文引用: 1]

李振林, 王晶, 2013.

ASTER GDEM与实测数据精度对比及其影响因素分析

[J]. 测绘与空间地理信息, 36(11): 150-153.

[本文引用: 1]

林爱兰, 谷德军, 彭冬冬, , 2021.

近60年我国东部区域性持续高温过程变化特征

[J]. 应用气象学报, 32(3): 302-314.

[本文引用: 1]

刘琳, 徐宗学, 2014.

西南5省市极端气候指数时空分布规律研究

[J]. 长江流域资源与环境, 23(2): 294-301.

[本文引用: 1]

刘强军, 梁赟, 宋军芳, , 2020.

近42年山西省晋城市地温与气温时空特征研究

[J]. 江苏农业科学, 48(6): 240-245.

[本文引用: 1]

罗伯良, 李易芝, 2014.

2013 年夏季湖南严重高温干旱及其大气环流异常

[J]. 干旱气象, 32(4): 593-598.

DOI      [本文引用: 1]

利用湖南省97个台站降水气温资料和NECP/NCAR再分析资料以及NOAA卫星观测的OLR场资料,系统分析了2013年6月30日至8月14日湖南持续高温干旱的时空分布及同期大气环流异常特征。结果表明:2013年夏季是湖南1951年以来降水最少高温干旱程度最严重的一年西太平洋副热带高压强度偏强西伸脊点偏西,湖南处在副高控制下,盛行下沉气流是发生持续高温干旱的直接原因;西太平洋副高内增强的下沉气流,致使多数县市高温日数和极端高温突破历史同期最高记录位于菲律宾附近的OLR低值带中心偏西,致使副热带地区下沉区偏西,进而造成西太平洋副高偏强偏西此外,湖南地区上空为水汽输送异常辐散中心,这在一定程度上促使了干旱的发生发展。

罗伯良, 彭莉莉, 张超, 2010.

2009年盛夏湖南持续高温干旱及同期大气环流异常分析

[J]. 干旱气象, 28(1): 20-25.

[本文引用: 1]

&nbsp;利用湖南省97个台站降水、气温资料、NCEP/NCAR再分析资料,分析了2009年8月湖南持续高温干旱的时空分布及同期大气环流的异常特征。结果表明:高温时段从8月中旬初开始一直持续到下旬后期,其中8月中旬到下旬高温范围和持续时间大部分地区为1959年以来最严重的时段。西太平洋副热带高压面积偏大、强度偏强、西伸脊点偏西,湖南地区在副高控制下,盛行下沉气流是引起持续高温干旱的直接原因。从南海到湖南地区存在经向水汽输送的负异常中心,这种水汽输送形势有利于持续高温干旱的发展。

马诺, 唐冶, 常存, , 2020.

基于DEM的气温空间插值方法比较

[J]. 干旱气象, 38(3): 465-471.

[本文引用: 1]

莫跃爽, 周秋文, 2019.

贵州省1960—2014年不同地貌类型的气温变化特征

[J]. 水土保持研究, 26(4): 166-170.

[本文引用: 1]

王钊, 彭艳, 魏娜, 2016.

近52 a秦岭南北极端温度变化及其与区域增暖的关系

[J]. 干旱气象, 34(2): 269-275.

DOI      [本文引用: 1]

利用中国气象局提供的地面气象站基本气象要素日值数据集(V3.0)中均一化温度数据,分析了1961~2012年陕西境内秦岭山脉南北两侧4个地貌单元平均最高温度和最低温度的趋势分布特征,同时采用极端气候指标计算软件RClimdex计算了5种极端温度指数,并分析其变化特征及其与区域增暖的关系。结果表明:不同季节秦岭地区极端温度变化存在较大的差异,平均最高温度春季增暖信号最明显,而平均最低温度冬季升温明显,不同的增暖趋势导致了秦岭地区春季、秋季气温日较差变大,冬季、夏季气温日较差变小。陕北黄土高原、关中盆地、秦岭南坡和汉水流域平均最高、最低温度变化趋势基本一致,但变化幅度存在一定差异,其中秦岭北部黄土高原和关中盆地平均最低、最高温度的变化幅度均大于南部的秦岭南坡和汉水流域,尤其平均最低温度关中盆地增幅更明显。秦岭北部2区域极端最低温度相关指数的变化幅度大于极端最高温度指数,而南部2区域前者的变化幅度小于后者。秦岭山脉区域增暖与平均最高、最低温度变化密切相关,还受极端温度变化的影响,其北部地区增暖主要是暖夜增加的贡献,而南部地区增暖主要与暖昼增加有关。

温晓培, 吴炜, 李昌义, , 2022.

土地利用资料的更新对四川盆地高温天气数值模拟的影响

[J]. 干旱气象, 40(5): 868-878.

DOI      [本文引用: 1]

土地利用情况影响区域天气和气候,采用与实况更为接近的下垫面资料能够为数值预报提供有力支撑。使用中尺度数值模式WRF对2019年8月四川盆地高温过程进行数值模拟试验,对比基于2019年更新的MODIS土地利用数据与WRF自带土地利用资料的差异,分析下垫面资料的更新对数值模拟效果的影响。结果表明:两组土地利用资料存在明显差异,更新后的MODIS资料刻画更为精细,模拟区域的土地利用类型更加丰富。2 m气温对土地利用资料十分敏感,更新土地利用资料后两次个例平均的2 m气温2 ℃以内的模拟准确率提高6.2%,日最高气温模拟准确率提高31.3%;而日最低气温模拟准确率下降2.1%。土地利用资料变化明显的四川盆地内部较更新前增温明显,局部增温超过4 ℃,最高气温模拟偏低的状况明显改善。选取四川盆地内部地表类型由农田转为城市和建筑的合川站作为典型站点进行分析,发现土地利用类型更新后发射率减小,植被气孔阻抗增大,粗糙度增大,陆面参数的变化导致地表对边界层大气向上的感热输送增加,地表蒸发及大气水分供应减少,气温升高,边界层高度升高。模拟初期土地利用资料的更新造成模拟差异的局地性较强,随着模拟时间增加,逐渐影响中高层及周边大气。土地利用资料的改变通过影响陆面参数的选择,导致模拟结果的差异。更新土地利用资料后可有效改善2 m气温及最高气温的模拟效果。

徐伟燕, 孙睿, 周爽, , 2017.

基于遥感和GIS的日最高最低气温估算

[J]. 北京师范大学学报, 53(3): 344-350.

[本文引用: 1]

余卫国, 房世波, 齐月, , 2019.

ASTER卫星数据地表温度产品精度评价

[J]. 干旱气象, 37(6): 987-992.

[本文引用: 1]

张超, 彭莉莉, 黄晚华, , 2012.

1961—2010年湖南气候变化特征及其对烟草种植的影响

[J]. 湖南农业大学学报, 38(5): 482-486.

[本文引用: 1]

张富荣, 刘志鹏, 冯淑霞, , 2015.

朝阳地区近30年地表温度的变化特征

[J]. 贵州农业科学, 43(5): 220-222.

[本文引用: 1]

郑嘉雯, 髙丽, 任宏利, , 2019.

基于T639集合预报的我国2016年极端温度预报检验

[J]. 气象, 45(4): 469-482.

[本文引用: 2]

周成虎, 程维明, 钱金凯, , 2009.

中国陆地1:100万数字地貌分类体系研究

[J]. 地球信息科学学报, 11(6):707-724.

[本文引用: 1]

地貌分类体系是地貌图研制的关键之一,本文在总结国内外地貌及分类研究的基础上,借鉴20世纪80年代的中国1∶100万地貌图制图规范,基于遥感影像、数字高程模型和计算机自动制图等技术条件,归纳总结了数字地貌分类过程中应遵循的几大原则,分析了它们之间的相互关系,讨论了数字地貌分类的各种指标:包括形态、成因、物质组成和年龄等,提出了中国陆地1∶100万数字地貌三等六级七层的数值分类方法,扩展了以多边形图斑反映形态成因类型,以点、线、面图斑共同反映形态结构类型的数字地貌数据组织方式,并详细划分了各成因类型的不同层次、不同级别的地貌类型。中国1∶100万数字地貌分类体系的研究,为遥感等多源数据的陆地地貌解析和制图提供了规范,也为《中华人民共和国地貌图集》的编制奠定了基础,同时为全国大、中比例尺地貌图的分类和编制研究提供了借鉴。

朱怀松, 刘晓锰, 裴欢, 2007.

热红外遥感反演地表温度研究现状

[J]. 干旱气象, 25(2): 17-21.

[本文引用: 1]

地球表面温度是一个重要的水文、气象参数,它影响着大气、海、陆之间的感热和潜热交换,是诸多研究领域不可或缺的基础资料,精确定量反演陆面温度的成果将推动旱灾预报和作物缺水研究、农作物产量估算、数值天气预报、全球气候变化和全球碳平衡等领域研究的进展。因此,利用卫星遥感资料进行地表温度的反演已成为目前遥感定量研究中的重要任务之一。本文阐述了遥感反演地表温度的原理及各种方法,对各种方法所要解决的关键问题及优缺点做了评述,最后展望了遥感反演地表温度的发展趋势。

邹德全, 邹承立, 田洪进, , 2021.

遵义市极端气温指数的时空变化特征

[J]. 贵州农业科学, 49(2): 146-153.

[本文引用: 1]

GOOD E, 2015.

Daily minimum and maximum surface air temperatures from geostationary satellite data

[J]. Journal of Geophysical Research: Atmospheres, 120(6): 2 306-2 324.

GEORGY V M, ROGER L K, RAJA R, et al, 2014.

Statistical estimation of daily maximum and minimum air temperatures from MODIS LST data over the State of Mississippi

[J]. GIScience Remote Sensing, 43(1): 78-110.

[本文引用: 1]

JONES P, JEDLOVEC C, SUGGS R, et al, 2004.

Using MODIS LST to estimate minimum air temperatures at night

[C]// the 13th Conference on Satellite Meteorology and Oceanography, American Meteorological Society, Norfolk: 13-18.

[本文引用: 1]

KAWASHIMA S, ISHIDA T, MINOMURA M, et al, 2000.

Relations between surface temperature and air temperature on a local scale during winter nights

[J]. Journal of Applied Meteorology, 39(9): 1 570-1 579.

[本文引用: 1]

LIANG S L, 2001.

Narrowband to broadband conversions of land surface albedo: I.algorithms

[J]. Remote Sensing of Environment, 76(2): 213-238.

[本文引用: 2]

SHENG P L, NATHAN J M, JOSEPH P M, et al, 2012.

Evaluation of estimating daily maximum air temperature with MODIS data in east Africa

[J]. International Journal of Applied Earth Observation and Geoinformation, 18: 128-140.

[本文引用: 1]

SUN Y, ZHANG X, ZWIERS F W, et al, 2014.

Rapid increase in the risk of extreme summer heat in Eastern China

[J]. Nature Climate Change, 4(12): 1 082-1 085.

[本文引用: 1]

TANG R, LI Z L, LIU M, et al, 2022.

A moisture-based triangle approach for estimating surface evaporative fraction with time-series of remotely sensed data

[J]. Remote Sensing of Environment:An Interdisciplinary Journal, 280, 113212. DOI:10.1016/j.rse.2022.113212.

[本文引用: 1]

/