基于小波变换与神经网络的石羊河流域夏季地温预测模型研究
Prediction model of summer land surface temperature in the Shiyang River basin based on the wavelet transform and neural network
通讯作者:
编委: 周成林
收稿日期: 2018-10-11 修回日期: 2020-02-15
基金资助: |
|
Received: 2018-10-11 Revised: 2020-02-15
作者简介 About authors
贾东于(1989-),男,河北邢台人,讲师,2018年在中国科学院西北生态环境资源研究院获博士学位,从事陆面过程与数值模拟研究.E-mail:
关键词:
Keywords:
本文引用格式
贾东于, 李开明, 聂晓英, 袁春霞, 李清峰, 高福元.
JIA Dongyu, LI Kaiming, NIE Xiaoying, YUAN Chunxia, LI Qingfeng, GAO Fuyuan.
0 引言
地表温度(LST, Land Surface Temperature)对研究水平衡物理过程和地表能量平衡至关重要, 其空间分布和时间演化对水文、 气象及气候应用的研究是至关重要的[1-3]。地温已被广泛应用于森林火险评估[4]、 城市热岛(UHI, Urban Heat Island)[5-6]和多年冻土[7]的监测研究中。同时, 地表温度的升高或降低可能改变环境气候的影响因子(例如加速融雪、 沙漠化及增加降水), 而反过来又可以增强(正反馈)或减弱(负反馈)这种变化。因此, 地温变化在气候反馈效应中起着重要作用。另外, 在气候变化研究中, 理解地温及其与影响因素之间的时空关系对预测全球温度变化至关重要[8]。李帅等[9]和吉春荣等[10]指出, 浅层地温会显著影响到农作物的耕种, 由于浅层地温会受到土壤性质、 土壤含水量以及辐射等诸多因素的影响, 成为预报中的难点。同时, 日平均、 日最高和日最低地温也被广泛应用在各种生态模型和陆面模型中, 是模型中的重要参数之一[11]。
目前, 对地温的研究主要集中在三个方面: 影响气候变化[12-14]、 预报预测和遥感反演[15-16]。其中, 地温预报主要是通过气象因子构建地温预报模型来进行研究[12]。宋晓巍等[17]利用气温建立了不同深度地温的预测模型, 研究发现在气温变化较大时, 模拟值与实测值误差较大, 这有可能是由于地温的延迟效应。Maduako等[16]首先利用遥感图像(Landsat 7 ETM)得到地温, 之后采用神经网络算法利用过去的地温值预测尼日利亚伊科姆市未来的地温, 其研究结果证实了神经网络在学习、 理解非线性复杂数据集和准确预测方面的有效性。吴春晓等[18]利用神经网络算法建立了浅层、 深层地温的推演模型, 研究发现70%以上的平均反演误差在0.5 ℃以下。同时, 在对日最低地表温度和日最高地表温度的模拟中发现神经网络的模拟效果明显好于逐步回归。Aliihsand等[19]利用人工神经网络(ANN)和Landsat-8卫星影像的时间序列分析, 在干旱环境的局部尺度上模拟地表温度日变化(LST), 结果表明神经网络方法将是预测和模拟干旱环境中昼夜地表温度变化的有力工具。Zhan等[20]利用三个输入节点和三个输出节点构造了一个两层BP神经网络模型, 成功地利用气象观测数据得到了青海省青海湖附近的地温数据。Bozorgi等[21]以大伊斯法罕市为例, 利用人工神经网络算法, 探讨城市地温与绿地覆盖空间格局之间的非线性关系, 研究发现神经网络算法能比较好的估算出城市地温。综上可知, 神经网络方法能够比较准确的模拟地温, 但目前研究主要关注浅层日平均地温, 并未对日最高、 日最低地温的预测进行探讨, 尤其缺乏对地气温差显著的干旱、 半干旱地区(例如石羊河流域)的研究。
本文利用小波变换结合BP(Back Propagation)神经网络的方法, 通过常规气象观测资料来预测石羊河流域地温(日平均地温, 日最低地温和日最高地温)的变化趋势。同时, 对预测效果及误差也进行了分析, 希望能为有效弥补干旱、 半干旱区地温观测缺失及其与局地气候的关系提供一些参考。
1 资料与研究方法
1.1 研究区域概况
石羊河主要发源于祁连山, 消失于腾格里和巴丹吉林沙漠之间的民勤北部荒漠地区, 平均海拔1 300 ~ 5 000 m, 其是河西走廊内陆河中重要的水系之一。该流域地势南高北低, 南部为祁连山褶皱, 中部为走廊凹陷, 北部为阿拉善台地及北山断块(图1)。由于地处中国内陆地区, 该流域主要受温带大陆性气候作用, 干旱少雨[22]。其中, 流域上游年均气温低于6 ℃, 年降水量400 ~ 600 mm, 年蒸发量700 ~ 1 200 mm; 中游年均气温低于6 ~ 8 ℃, 年降水量150 ~ 300 mm, 年蒸发量1 300 ~ 2 000 mm; 下游年均气温高于8 ℃, 年降水量小于150 mm, 年蒸发量2 000 ~ 2 600 mm[22-24]。
图1
图1
石羊河流域及气象站点
Fig.1
Map showing the Shiyang River basin and the meteorological observation stations
1.2 数据与方法
1.2.1 数据介绍
气象数据来自国家气象信息中心制作的中国地面气候资料日值数据集(V3.0), 选取石羊河流域的民勤、 永昌、 武威三个国家基准、 基本气象站的观测数据。日平均气温、 相对湿度、 0 cm地温均为4次定时(02:00、 08:00、 14:00、 20:00, 北京时间, 下同)观测值的平均值, 其数据均进行了严格的质量控制。本文所选取的气象资料主要来自于石羊河流域的肃南、 永昌、 武威、 民勤、 古浪气象站1998 - 2017年6 - 8月的日观测数据。
图2
图2
石羊河流域地表下垫面特征
Fig.2
The underlying surface in the Shiyang River basin; Minqin (a), Yongchang (b) and Wuwei (c)
图3
图3
石羊河流域的日平均气温(a)、 日最高气温(b)、 日最低气温(c)和日降水概率(d)
Fig.3
Variation ranges of daily mean temperature (a), daily maximum temperature (b) and daily minimum temperature (c), and daily precipitation probability distribution (d) in the Shiyang River basin
1.2.2 研究方法
图4为小波神经网络预测地温的流程图。在本研究中通过对地温使用三层小波分解, 即得到逼近信号A3和细节信号D3、 D2和D1。经过相关性分析(表1), 本文选取与地温相关系数大于0.5的变量(空气温度、 气压、 相对湿度和蒸发量)作为训练样本的输入参数。训练样本为1998 - 2012年石羊河流域三个站点6 - 8月的日观测资料, 测试样本为2013 - 2017年6 - 8月地温的日观测资料。本文对神经网络模型选择的隐层节点数为7个。采用Matlab 2018a进行预测程序的建模, 将预测数值通过小波重构, 最终分别获得2013 - 2017年三个台站日平均地温、 日最高地温和日最低地温的预测值。其中, 小波变换和小波神经网络方法在下文详细阐述。
图4
图4
小波神经网络预测地温流程
Fig.4
The forecasting process of LST by wavelet transform and neural network
表1 石羊河流域地温与气象要素的相关系数
Table 1
气象要素 | 日平均地温 | 日最高地温 | 日最低地温 |
---|---|---|---|
蒸发量/mm | 0.62 | 0.54 | 0.09 |
气压/hPa | 0.62 | 0.46 | 0.34 |
气温/℃ | 0.88 | 0.71 | 0.41 |
相对湿度/% | 0.67 | 0.70 | 0.18 |
降水/mm | 0.35 | 0.40 | 0.03 |
风速/(m·s-1) | 0.02 | 0.12 | 0.22 |
水平辐射 | 0.15 | -0.12 | 0.06 |
(1)小波变换
式中: c是母小波的频率, i表示虚数。它的傅里叶变换如下:
对于任意函数
式中:
图5
图5
民勤日平均地温的小波分解
Fig.5
Wavelet decomposition of the daily average land surface temperature in Minqin
(2)小波BP神经网络
本文采用小波分析结合BP(Back Propagation)神经网络模型。人工神经网络可视为能够模拟生物神经系统某些行为的数学模型[25]。这个术语的现代定义是一种人工构造, 其行为是基于简单的神经元网络。神经元与权重连接在一起, 以便它们可以协作地处理信息并将这些信息存储在这些权重上。神经网络与其他模型相比有很多优点, 主要包括: (1)能学习高度复杂的模型并对噪声进行稳健的训练。同时, 少量的观测误差对模型影响不大。(2)具有较强适应能力[28]。这些特性使得神经网络成为一种重要的人工智能技术。BP神经网络是一种前馈神经网络, 其中神经元是分层的, 每个神经元只能与下一层神经元连接。因此, 在这样的网络中不存在定向循环。
图6
输出层的激活函数
BP网络可以从一组样本数据(包括输入和相应的预期输出)中学习复杂的非线性输入 - 输出函数。图6显示了BP神经网络的结构。
2 结果与分析
石羊河流域的武威、 永昌和民勤三地分布在石羊河流域的三个主要区域, 其下垫面地表类型和人口数量也有所差异, 因此, 选取这三个观测点能够较好反映该方案在整个石羊河流域的适用性。
其中, 图7 ~ 9为武威、 民勤、 永昌三地地温预测值和观测值的比较, 图9为日平均、 日最高和日最低地温预测的误差概率分布图。由图7 ~ 9可知, 三个气象站日平均地温的预测均为最佳, 其中预测值和观测值的相关系数均大于0.87, 其预测效果优于日最高和日最低地温。同时, 结合图10(a)可知, 三个站点90%以上的平均地温预测误差分布在-4 ~ 4 ℃之间, 3 ℃以内的预报概率均大于84%, 误差区间较小并呈现标准正态分布。其中, 人口密度较少和地表植被覆盖度较少的民勤地区预测效果最好, 这主要是因为该地区人为活动对地温的影响较小, 地温和气象因子关系密切。民勤地区日平均地温的预测值和模拟值的相关系数达到0.91, 3 ℃以内的预测概率达到86%, 而误差平方和、 均方根误差为三个站点中的最小值。永昌地区日平均地温预测值和观测值的相关系数虽然只有0.87, 但是预报误差在3 ℃以内的概率达到96%, 说明对该区域地温数值大小的预测效果较好。而武威地区日平均地温预测值和观测值的相关系数达到0.91, 说明预测值与观测值的变化趋势能够基本吻合。
图7
图7
小波神经网络预测石羊河流域日平均地温效果
Fig.7
The simulated and measured summer daily average land surface temperatures in Minqin (a), Yongchang (b), Wuwei (c) and the simulated error of land surface temperatures (d) in the Shiyang River basin
图8
图8
小波神经网络预测石羊河流域日最高地温效果
Fig.8
The simulated and measured summer daily maximum land surface temperatures in Minqin (a), Yongchang (b), Wuwei (c), and the simulated error of summer daily maximum temperatures (d) in the Shiyang River basin
图9
图9
小波神经网络预测石羊河流域日最低地温效果
Fig.9
The simulated and measured summer daily minimum land surface temperatures in Minqin (a), Yongchang (b), Wuwei (c), and the simulated error of summer daily minimum temperatures (d) in the Shiyang River basin
图10
图10
石羊河流域地温预测误差正态分布
Fig.10
Normal distribution of land surface temperature prediction error in the Shiyang River basin: daily average land surface temperature (a), daily maximum land surface temperature (b), daily minimum land surface temperature (c)
日最高地温的预测结果显示, 模型对日最高地温的变化趋势也有较好的模拟效果。其中, 民勤、 永昌和武威地区预测值与观测值的相关系数分别为0.82、 0.83和0.84。结合表2可知, 永昌地区的模拟效果最好, 3 ℃以内的预测概率达到83%。但是由图8(b)可知, 在2016年后模型对永昌地区日最高地温高估的概率较大, 其中的原因还值得更深入的探究。而民勤、 武威两地的日最高地温的变化区间较大, 变化范围为15.2 ~ 74.3 ℃, 因此预测值与实际观测数值的误差平方和、 标准差也较大。同时观察图8(c)可知, 对武威地区的日平均地温出现较多低估, 参考城市和人口规模可知, 这是由于武威地区的人口规模远大于永昌和民勤地区, 城市热岛效应会对地表温度也产生一定影响, 因此对武威地区的地温预估仅利用气象因子会出现一定偏差, 未来应当加入人为活动因子来修正该区域的地温预测。
表2 石羊河流域地温预报的准确性评估
Table 2
预报变量 | 相关系数R | 误差平方和SSE | 标准差RMSE | 预测误差3 ℃以内 |
---|---|---|---|---|
民勤日平均地温 | 0.91 | 2 294 | 2.238 | 86% |
民勤日最高地温 | 0.82 | 19 270 | 6.486 | 61% |
民勤日最低地温 | 0.71 | 2 629 | 2.396 | 87% |
永昌日平均地温 | 0.87 | 2 872 | 2.504 | 96% |
永昌日最高地温 | 0.83 | 21 750 | 6.891 | 83% |
永昌日最低地温 | 0.66 | 2 372 | 2.276 | 83% |
武威日平均地温 | 0.91 | 2 655 | 2.408 | 84% |
武威日最高地温 武威日最低地温 | 0.84 0.72 | 22 930 2 055 | 7.075 2.118 | 52% 94% |
3 讨论与结论
3.1 讨论
由于本研究为了实现预测的普适性并未剔除复杂天气条件下的观测数据, 而地温与气温等气象要素之间的关系在晴天、 多云天和阴雨天时存在较大差异, 这就会在一定程度上影响神经网络的学习和预测效果。另外, 由于石羊河流域国家基本气象站点分布密度有限, 符合本研究的站点较少, 故研究结论可能有一定偏差。下一步的研究中, 拟将不同的天气类型进行区分, 分别进行地温的精细化预报, 并将结合遥感数据和野外台站数据, 弥补观测台站较少对神经网络训练造成的影响。
3.2 结论
本文利用小波分析结合神经网络通过武威、 永昌、 民勤三个气象台站的4种常规气象观测资料来预测石羊河流域夏季地温的变化, 分析了石羊河流域日平均、 日最高和日最低地温, 实现了利用常规气象数据对地温的预测, 该方法简明有效, 对农业生产十分重要。主要结论如下:
(1) 利用小波变换结合BP神经网络的方法, 通过常规气象观测资料预测石羊河流域地温(日平均地温, 日最低地温和日最高地温)的整体效果较好, 能够反映该研究区域5年夏季地温的变化趋势。其中, 3 ℃以内的预报概率均大于80%, 这说明本方法能够准确的预测夏季地温的变化。
(2) 日平均地温的预测效果明显优于日最高和日最低地温。这与输入层中气象要素的时间分辨率密切相关, 由于输入的气象数据均为日平均数据, 并不能从时间尺度有效对应日最高、 日最低地温的变化, 因此会造成一些误差。然而, 本研究预测结果与实际观测的相关系数均大于0.66, 说明本方法仍可以适用于日最高和日最低地温的预测, 但在未来需要进行改进。
(3) 石羊河流域地温的预测效果在温度较高区域(例如民勤和武威)优于温度较低区域(例如永昌)。同时, 本方法对日最高地温的预测效果优于日最低地温。通过文中分析发现, 该方案对石羊河流域的日最低地温存在较大概率的低估, 应当在未来研究中考虑添加其他相关变量, 提高对日最低地温预测的准确性。
参考文献
An analysis of GOES and NOAA derived land surface temperatures estimated over a boreal forest
[J].
Generation of a time-consistent land surface temperature product from MODIS data
[J].
Responses of soil moisture and thermal conductivity to precipitation in the mesa of the Loess Plateau
[J].
Assessment of seasonal forest fire risk using NOAA-AVHRR: a case study in central Mexico
[J].
Remote sensing image-based analysis of the relationship between urban heat island and land use/cover changes
[J].
Estimation of land surface temperature-vegetation abundance relationship for urban heat island studies
[J].
Using the MODIS land surface temperature product for mapping permafrost: an application to northern Québec and Labrador, Canada
[J].
Spatio-temporal variability in remotely sensed land surface temperature and its relationship with physiographic variables in the Russian Altay Mountains
[J].
Variation features and prediction of the topsoil (0 ~ 20 cm) ground temperature of spring in Heilongjiang Province
[J].
黑龙江省春季浅层(0 ~ 20 cm)地温变化特征及预报
[J].
The relationship of ground temperature at 5 cm soil depth of plastic mulching with air temperature and its prediction in cotton area of Shihezi Oasis
[J].
石河子绿洲棉区覆膜内5 cm地温变化与气温的关系及预报
[J].
Empirical models for estimating daily maximum, minimum and mean air temperatures with MODIS land surface temperatures
[J].
Estimating evapotranspiration using improved fractional vegetation cover and land surface temperature space
[J].
Prediction method of land surface temperature based on numerical forecast products
[J].
基于数值预报产品的地表温度预报方法研究
[J].
Managing land use transformation and land surface temperature change in Anyigba Town, Kogi State, Nigeria
[J].
A preliminary view on the estimation of land surface temperature under cloud cover from thermal remote sensing data
[J].
热红外遥感图像中云覆盖像元地表温度估算初论
[J].
Prediction of land surface temperature (LST) changes within Ikom City in Nigeria using artificial neural network (ANN)
[J].
Change characteristics of shallow ground temperature and its prediction during the sowing period in Shenyang area
[J].
沈阳地区浅层地温变化特征及播种期地温预报研究
[J].
Ground temperature deduction model based on BP neural network
[J].
基于BP神经网络的地温推演模型
[J].
Modeling diurnal land surface temperature on a local scale of an arid environment using artificial neural network (ANN) and time series of Landsat-8 derived spectral indexes
[J].
The study of the land surface temperature retrieval with the proposed neural network model
[C]//
Land surface temperature estimating in urbanized landscapes using artificial neural networks
[J].
Long-term trends of major climatic variables and runoff in the Shiyang River basin
[J].
石羊河流域主要气象要素及径流变化趋势分析
[J].
Drought characteristics in the Shiyang River basin during the recent 50 years based on a composite index
[J].
基于综合气象干旱指数的石羊河流域近50年气象干旱特征分析
[J].
Temporal and spatial variation of vegetation index and remote sensing monitoring of desertification in Shiyang River basin
[J].
石羊河流域植被指数时空变化及荒漠化遥感监测
[J].
Wave parameters and functions in wavelet analysis with filtering
[J].
A power prediction method for photovoltaic power plant based on wavelet decomposition and artificial neural networks
[J].
Fundamentals of neural networks: architectures, algorithms, and applications
[M].
BP neural network-based effective fault localization
[J].
Distribution characteristics of 80 cm ground temperature and prediction method based on solar radiation
[J].
全国80 cm地温分布特征和基于太阳辐射的预测方法
[J].
Statistical analysis of the relationship between ground temperature and temperature prediction model in Wuhan
[C]//
武汉市地气温关系的统计分析与地温预测模型
[C]//
/
〈 |
|
〉 |
