Sentinel-1和Sentinel-2协同反演地表土壤水分
1.
2.
3.
4.
5.
6.
Sentinel-1 and Sentinel-2 Synergistic Retrieval of Surface Soil Moisture
1.
2.
3.
4.
5.
6.
通讯作者:
收稿日期: 2022-05-22 修回日期: 2022-10-26
基金资助: |
|
Received: 2022-05-22 Revised: 2022-10-26
作者简介 About authors
杜绍杰(1992-),男,河南开封人,硕士研究生,主要从事主动微波土壤水分反演研究E⁃mail:
关键词:
Keywords:
本文引用格式
杜绍杰, 赵天杰, 施建成, 马春锋, 邹德富, 王振, 姚盼盼, 彭志晴, 郑景耀.
Du Shaojie, Zhao Tianjie, Shi Jiancheng, Ma Chunfeng, Zou Defu, Wang Zhen, Yao Panpan, Peng Zhiqing, Zheng Jingyao.
1 引 言
土壤水分是影响陆地表层水循环、能量循环、碳循环中的重要参数[1-2],是研究陆表—大气之间能量和物质交换的关键变量[3-4]。表层土壤水分既受到气候变化的影响[5],又对气候变化具有反馈作用[6]。土壤水分在空间上具有高度异质性[7],目前主流的基于被动微波遥感影像获取的土壤水分产品,如SMOS(Soil Moisture Ocean Salinity)、SMAP(Soil Moisture Active Passive)等在小区域尺度不足以表征土壤水分的空间异质性,无法满足生态环境监测[8]、农作物评估[9]、洪水预警[10,11]等领域的应用需求,因此发展高分辨率的土壤水分反演方法具有重要意义。
针对裸地,研究者们已经建立了理论模型、经验模型以及半经验模型描述微波信号和土壤(随机粗糙面)的相互作用关系。传统的理论模型包括基尔霍夫近似模型(物理光学模型和几何光学模型是基尔霍夫近似在不同假设条件下的特例)、小扰动模型等,这些模型在实际应用中因适用粗糙度范围较窄而受到限制。以积分方程模型(Integral Equation Model, IEM)为代表的随机粗糙面散射模型因具有更大的粗糙度适用范围而获得了广泛的应用[14],其改进模型AIEM[15](Advanced Integral Equation Model, AIEM)更是在密度谱函数和光谱反射系数泛化方面具有独特的优势[16]。常见的经验模型有Oh模型[17-19]、Dubois模型[20]以及Shi等[21]利用IEM模型发展的半经验模型。
由于植被的影响,以上模型不能直接应用于植被覆盖地表的土壤水分反演。为描述植被与微波信号的相互作用,Ulaby等[22]提出的密歇根微波植被散射模型(Michigan Microwave Canopy Scattering, MIMICS),较完整地描述了森林冠层的散射特性,但该模型对植被冠层刻画较为复杂,模型参数较多难以在实际反演中进行应用。Attema和Ulaby提出考虑植被单次散射的水云模型(Water Cloud Model, WCM),由于模型较为简单,且能保证模拟精度,该模型广泛适用于农作物等低矮植被覆盖下土壤水分的反演[23]。随着植被增加,其对后向散射的贡献逐渐占据主导地位,土壤水分反演的难度与误差增大。
消除地表粗糙度和植被对雷达后向散射的影响是雷达土壤水分反演算法的关键[24-25]。学者们发展了查找表方法、多极化数据反演方法、迭代方法、基于神经网络的方法等一系列土壤水分反演方法[26],查找表法借助于微波散射模型建立不同地表参数下后向散射系数的表格,反演过程中根据实测后向散射系数查找对应的土壤水分信息。查找表法简单高效[27],但反演精度受限于所使用模型对地表后向散射描述的精度。基于多极化数据的方法可以利用不同极化雷达数据分离植被和粗糙度对后向散射的贡献[28-30],赵天杰等[31]将同极化与交叉极化结合进行地表粗糙度的估计以分离粗糙度对后向散射的影响,提升了土壤水分的反演精度。迭代方法构建代价函数并利用遗传算法、支持向量机、模拟退火等优化方法反演土壤水分[32]。近年来,学者利用神经网络在土壤水分反演中也取得了较好的效果,神经网络算法首先利用物理模型生成训练数据和验证数据,并建立后向散射系数与土壤水分之间的复杂非线性关系来反演土壤水分[33-34]。利用物理模型和神经网络算法大都需要精确的地表参数作为输入数据,但地表参数的空间异质性以及模型的欠精确性使得神经网络方法在获取大范围、高鲁棒性的反演结果仍然较为困难[35]。
以多时相观测数据为基础的变化检测等方法为解决土壤水分反演中的欠定问题提供了方案。地表粗糙度的变化相对于土壤水分的变化是缓慢的,因此假定在短观测周期内粗糙度不变,利用雷达多次观测数据,可以认为后向散射系数的变化是由土壤水分的变化引起[36]。Wagner等[37]基于ERS散射计建立后向散射系数和土壤水分之间的关系反演土壤水分,证明了变化检测法在土壤水分反演中的巨大潜力。对于植被覆盖地表,需要将植被对后向散射的影响从总的后向散射中剥离[38-39]。姜灵海等[40]利用水云模型校正植被覆盖下的后向散射系数进而利用变化检测法反演土壤水分,Zribi等[41]基于散射计数据通过建立归一化植被指数(Normalized Differential Vegetation Index, NDVI)和后向散射系数差之间的经验关系,有效消除了植被带来的影响。Gao等[42]提出两种消除植被影响的方法,一种基于观测周期内后向散射减去最干燥土壤对应的后向散射的差,另一种基于连续两个日期观测到的后向散射系数之间的差,并发现后者具有更好的反演效果。
为了建立简单、有效且合理的函数表征植被对后向散射影响的关系,本文将后向散射差表示为任一时刻观测与该观测周期内最小后向散射的差,同时最小后向散射系数对应于该周期内最小土壤水分。因植被覆盖地表下得到的后向散射最大差值不总是对应着裸土条件下的后向散射最大差值,故利用NDVI对后向散射差值的植被影响进行校正,得到裸土表面的后向散射最大差值。在五道梁地区与土壤水分地面观测数据对比评估该模型的适用性。
2 研究区与数据源
2.1 研究区概况
图1
研究区内部地形地貌较为复杂,如图1(b),研究区域海拔在4 523—5 119 m。北部呈线状分布的低海拔条带为楚玛尔河河谷,河谷两岸地形较为平坦,主要为裸露和稀疏植被地表,土壤质地主要以含有砂粒较多的砂土为主。中南部为山丘,夏季山丘土壤融化,为山谷以及两侧提供了大量水源,南北两侧平坦地带分布有大量的热熔湖河,山丘两侧较为平坦,海拔相近,但地表覆盖差异明显,山丘北部植被发育程度高,主要地表覆盖类型为高原草甸,在靠近山丘的数公里内为高原沼泽草甸,土壤类型主要为粘土。南部植被发育较低,主要以裸地及高原荒漠为主。
研究区为长宽36 km×36 km的方形区域,如图1(c),红色三角形标志为五道梁温湿度观测网络下的10个土壤水分观测站点,10个站点涵盖了研究区内典型的地表覆盖类型。
2.2 数据源
2.2.1 地面测量数据
五道梁土壤温湿度观测网络建立于2019年9月,观测站涵盖10个土壤湿度监测站点,观测仪器每30 min获取不同土壤深度的土壤水分数据。利用基于烘干法测量的土壤水分数据对观测仪器获取的数据进行校准。
土壤水分地面观测数据的时间覆盖范围为2019年9月至2020年10月。考虑Sentinel-1雷达卫星的过境时间以及C波段雷达的穿透特性,本文选取上午7:30和8:00两次5 cm深度的土壤水分观测平均值作为验证数据。
2.2.2 Sentinel-1
Sentinel系列卫星由欧空局(European Space Agency, ESA)为实施“哥白尼”计划而开发,目前已发展了涵盖雷达和光学卫星的4个系列卫星(Sentinel-1/2/3/5),在陆地和海洋监测、气候变化以及空气质量等研究中具有重要作用。
Sentinel-1是基于C波段的雷达成像系统,由两颗极轨卫星A星和B星组成。其双星系统在部分地区的重访周期可由12 d缩短到6 d。Sentinel-1提供4种极化产品且具有4种工作模式,分别为条带模式(Stripmap Model, SM)、干涉宽幅模式(Interferometric Wide Swath, IW)、超宽幅模式(Extra-Wide Swath,EW)和波模式(Wave Mode,WM)。本研究使用IW模式下的GRD产品。
为了避免不同传感器和不同轨道带来的误差,本文选取覆盖研究区的A星共60景降轨数据,极化方式为VV极化。数据起始范围为2019年1月至2020年12月,空间分辨率为5 m × 20 m、重访周期为12 d。GEE已经对来源于欧空局的原始数据进行如下处理:①热噪声去除;②辐射校正;③基于SRTM(Shuttle Radar Topography Mission, SRTM)30 m分辨率数据的地形校正。在GEE处理的基础上,本文对研究区60景影像依次进行:
(1)基于RefinedLee滤波去除相干斑噪声;
(2)局地入射角归一化:
地表的雷达后向散射对入射角度敏感。考虑到研究区地形地貌较为复杂,基于Sentinel-1的轨道参数结合SRTM 30 m分辨率数字高程模型计算了每一个像元的局地入射角。为了后续在同一个角度下与SMAP产品进行分析,根据
基中:θ为局地入射角;
(3)雷达图像升尺度至100 m分辨率
将100 m分辨率下的每个像元的值设定为对应的100个10 m分辨率像元的平均值,后向散射平均到100 m分辨率可以减弱像元内的散斑效应,并降低异常值的影响,增加数据的稳定性。最后将处理后的后向散射系数转换为dB单位。
2.2.3 Sentinel-2
Sentinel-2高分辨率多光谱成像卫星系统由运行在近极地轨道上的A、B两颗卫星组成。双星联合观测可以达到5 d的重访周期。其携带的多光谱成像仪(MSI)覆盖了13个波段,提供10 m分辨率的可见光和10/20 m分辨率的近红外数据。研究采用GEE云平台提供的Sentinel2 L2A产品,L2A产品由L1C数据经大气校正转换为大气底部反射率生成。数据时间从2019年1月至2020年12月共288景数据。
未经去云的像元使后续计算出现偏差,波段“QA60”是一个包含云掩码信息的波段,用于去除被云覆盖的区域,因此对每一景预处理后的数据去除云像元并计算NDVI。NDVI的计算公式如下:
其中:Red(Band-4)和NIR(Band-8)分别为红光和近红外波段的反射率。
去除云像元后的数据计算出来的NDVI仍然含有异常值,表现为某像元的NDVI值在时间序列上相比前一次观测明显升高随后在下次观测又明显降低,经过分析这些异常NDVI值对应于云层在地表投射的阴影,为了降低云阴影带来的误差,根据像元的NDVI时间序列信息,检索在时间序列上异常变化的像元并将其判定为云阴影去除。
将NDVI小于0的像元认为主要为水体,通过与光学数据验证对比,该方法可以较好的去除楚玛尔河河床和热熔湖对应的像元。为了与Sentinel-1雷达数据对应,研究选取了与Sentinel-1获取时间相差3 d内的NDVI数据,并将10 m分辨率的NDVI数据重采样到100 m分辨率。
原始NDVI数据在去除云及云阴影之后含有大量空值像元。NDVI代表植被的生长状况与覆盖度,因此在短时间内变化较小,本文选取NDVI空缺的日期前后两次观测的平均值进行空值填充。
3 反演方法
3.1 裸土条件下的变化检测法
裸露土壤的后向散射主要受土壤水分和土壤表面粗糙度的影响,当裸土粗糙度变化较小的时候,可以认为其对后向散射的贡献保持稳定。因此利用雷达对同一区域进行连续观测,可以忽略粗糙度对后向散射变化的影响,此时认为观测到的地表后向散射变化主要由土壤水分变化引起。因此建立后向散射变化对土壤水分变化的响应关系可以实现对土壤水分的反演。
假设裸土条件下雷达极化后向散射与土壤水分之间存在线性响应关系,有如下公式:
其中:
在重复观测周期内,当土壤湿度达到最大值
同样当土壤湿度达到最小值
观测周期内土壤水分分最大变化量
其中:
对于观测周期内第t天的土壤水分
其中:
根据后向散射的变化量可以获得土壤水分的变化量。
3.2 考虑植被对后向散射的影响
对于植被覆盖地表,后向散射的变化不仅受到土壤水分变化的影响,也受到地表植被覆盖的影响。假设将总的后向散射变化量表示为裸土下后向散射的变化量与植被影响的和:
其中:
如图2,结合观测到的植被覆盖地表后向散射系数变化量及对应的植被指数,即可求得相对于裸土条件下的后向散射系数变化量。观测到的植被覆盖地表后向散射系数变化量可由Sentinel-1观测数据计算得出,NDVI由Sentinel-2计算得出。
图2
图2
NDVI对后向散射系数变化量的影响关系
Fig.2
The influence of NDVI on the variation of backscattering coefficient
因此对于给定的NDVI下的后向散射变化量,通过
当
其中:
图3
图3
研究区域不同NDVI下后向散射变化量的趋势
Fig.3
The trend of backscattering variation under different NDVI in the study area
4 结果与分析
4.1 与地面站点观测数据对比
利用以上发展的变化检测方法进行土壤水分反演,为了得到土壤水分,从五道梁地面温湿度观测网络的土壤水分实测数据中得到需要的土壤水分最小值和最大值。并将10个观测站点对应的100 m分辨率土壤水分反演结果和地面观测数据进行对比。如图4所示,本文中提出的变化检测法表现出了良好的性能,反演值与实测值相关系数达到0.904,反演结果相对实测含水量较大,均方根误差为0.062 m3/m3。
图4
图4
土壤水分反演值与地面观测值对比
Fig.4
Comparison of soil moisture inversion values and ground observations
4.2 站点时间序列对比分析
图5为研究区内10个站点土壤水分观测值与反演值在观测周期内的变化趋势。从图中可知,研究区内观测站点的土壤水分在一年中季节性波动明显,同一时期不同的站点之间土壤水分差异较大,体现了土壤水分空间分布的异质性,这种差异与植被覆盖等相关[46]。同时说明发展植被校正下的高空间分辨率土壤水分算法和产品具有重要的意义。冬季研究区表层土壤未冻水分变化较为平稳,本文反演的为冻土中的未冻水,10个站点的最大土壤水分均不超过0.1 m3/m3。5月份随着气温上升,冻土融化[47],10个观测站点的土壤水分均开始增加,其中站点7和站点8在观测周期内土壤水分最大可达到0.5 m3/m3以上。其他站点夏季最大土壤水分在0.25 m3/m3附近。
图5
图5
10个站点土壤水分反演值与地面观测值的时间序列验证
Fig.5
Time series verification of soil moisture inversion values and ground observations at 10 sites
变化检测法反演的土壤水分在10个站点均能表达出土壤水分在一年中的变化趋势,表1展示了10个点的反演值与地面观测值之间的验证指标。10个站点反演的土壤水分值均与地面观测值相关性较好。2号站点在10个站点中的相关系数最小(R=0.692),其他9个地面观测站点相关系数都达到或接近0.8,特别是站点7和站点8相关系数分别为R=0.915和R=0.944,站点7和站点8位于研究区中部的高原草甸,这两个站点在观测周期内土壤水分最大分别达到0.60 m3/m3和0.50 m3/m3,表明变化检测算法在高土壤水分区间仍具有良好的适用性。
表1 10个站点土壤水分反演值与地面观测值的验证指标
Table 1
站点 | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 |
---|---|---|---|---|---|---|---|---|---|---|
R | 0.806 | 0.692 | 0.849 | 0.785 | 0.811 | 0.772 | 0.915 | 0.944 | 0.856 | 0.876 |
RMSE/(m3/m3) | 0.032 | 0.049 | 0.020 | 0.082 | 0.032 | 0.078 | 0.110 | 0.082 | 0.027 | 0.038 |
ubRMSE/(m3/m3) | 0.031 | 0.049 | 0.017 | 0.047 | 0.029 | 0.052 | 0.072 | 0.050 | 0.027 | 0.036 |
夏季,土壤水分对后向散射的影响占据主导因素,不同的站点的最大土壤水分值差异较大,站点1最大土壤水分为0.18 m3/m3,站点7的土壤水分最大值则可以达到0.60 m3/m3,不同站点的反演值均能很好的表达土壤水分的变化趋势,根据地面实测数据,站点4属于砂粒土质,站点5位于研究区的一处沙丘旁,两个站点土壤持水能力较弱,其土壤水分的变化主要由降雨引起,且其夏季土壤水分表现出明显的波动,其中站点4土壤水分观测数据在几天之内从0.3 m3/m3下降到0.1 m3/m3,本文的反演算法可以很好的捕捉到地面土壤水分的变化。整个观测周期内两个站点的相关系数分别达到R=0.785和R=0.811,均方根误差RMSE=0.082 m3/m3和RMSE=0.032 m3/m3。
4.3 研究区平均土壤水分变化分析
为了进一步评估变化检测法的反演结果,将2019年10月至2020年10月研究区域内10个站点的土壤水分反演数据和地面观测数据求平均值进行分析。如图6蓝色点为10个站点土壤水分反演值的平均值,黑色曲线为对应时期10个地面观测站点记录的土壤水分平均值,从时间序列上分析,反演数据在冬季表现出对土壤水分的高估,反演数据在夏季与地面观测数据相关性较好,整个观测周期内二者相关系数达到0.939,均方根误差为0.044 m3/m3,进一步验证了变化检测法进行土壤水分估算的潜力。图中绿色线段为国家气象局五道梁气象站的降雨量数据,数据来自于国家气象科学数据中心(http:∥data.cma.cn/)。在整个研究区,土壤水分的变化明显受降雨驱动,冬季降雨较少,土壤水分整体平稳且在0.05 m3/m3附近,4—5月份随着降雨增加,土壤逐渐变湿,降雨集中在6月至8月份,降雨量最大接近20 mm,土壤水分在此期间达到最大值0.287 m3/m3。这与土壤水分地面观测数据和雷达反演数据的趋势吻合。
图6
5 结 论
雷达遥感是获取高分辨率土壤水分的重要方式,本研究针对传统变化检测方法对植被影响考虑的不足之处进行改进,将雷达后向散射的变化表示为植被变化与裸土后向散射变化的函数关系,并认为年内观测到的后向散射变化的最大差值并不对应与裸土条件下后向散射变化的最大差值,提出了一种改进的土壤水分变化检测方法。
在青藏高原五道梁地面观测站对改进后的变化检测法进行验证,该方法在10个地面观测站点均具有较高的相关性,可以敏锐的捕捉到短时间内降雨带来的土壤水分的变化。在冬季土壤水分反演值出现高估,可能与冬季土壤冻结程度差异导致土壤介电常数变化以及粗糙度受冻结形变的影响发生改变有关。
研究提出的变化检测法假设了在观测周期内土壤粗糙度不变,但是在实际的应用中,土壤粗糙度往往不是稳定不变的。未来的研究中若可以将土壤粗糙度对后向散射系数的改变纳入考虑,将会进一步改善变化检测算法的反演精度。同时,本文在植被校正中基于NDVI建立经验函数、不同植被指数的选择影响植被校正的结果[48],同时该经验函数依赖于特定研究区,在其他地区是否适用需要进一步的分析。基于实测站点数据进行干湿参考值的确定,未来可考虑利用被动土壤水分产品获得区域干湿参考值以实现区域土壤水分的反演。
参考文献
Investigating soil moisture–climate interactions in a changing climate: A review
[J].
Recent advances of L-band application in the passive microwave remote sensing of soil moisture and its prospects
[J].
被动微波反演土壤水分的L波段新发展及未来展望
[J].
Soil moisture retrievals using L-band radiometry from variable angular ground-based and airborne observations
[J].
Ground observation-based analysis of soil moisture spatiotemporal variability across a humid to semi-humid transitional zone in China
[J].
Soil moisture: A central and unifying theme in physical geography
[J].
Impact of soil moisture‐climate feedbacks on CMIP5 projections: First results from the GLACE‐CMIP5 experiment
[J].
On the spatio-temporal dynamics of soil moisture at the field scale
[J].
Soil moisture measurement for ecological and hydrological watershed-scale observatories: A review
[J].
Assimilation of remotely sensed soil moisture and vegetation with a crop simulation model for maize yield prediction
[J].
Building a flood-warning framework for ungauged locations using low resolution, open-access remotely sensed surface soil moisture, precipitation, soil, and topographic information
[J].
Impact of different satellite soil moisture products on the predictions of a continuous distributed hydrological model
[J].
Retrieval of high-resolution soil moisture through combination of Sentinel-1 and Sentinel-2 data
[J].
Progress in inversion of surface parameters from microwave remote sensing
[J].
微波遥感地表参数反演进展
[J].
Progress and development trend of soil moisture microwave remote sensing retrieval me-thod
[J].
土壤水分微波反演方法进展和发展趋势
[J].
Backscattering from a randomly rough dielectric surface
[J].
The emissionof rough surfaces calculated by the integral equation method with a comparison to a Three-dimensional moment method simulations
[J].
A transition model for the reflection coefficient in surface scattering
[J].
An empirical model and an inversion technique for radar scattering from bare soil surfaces
[J].
Semi-empirical model of the ensemble-averaged differential Mueller matrix for microwave backscattering from bare soil surfaces
[J].
Quantitative retrieval of soil moisture content and surface roughness from multipolarized radar observations of bare soil surfaces
[J].
Measuring soil moisture with imaging radars
[J].
Estimation of bare surface soil moisture and surface roughness parameter using L-band SAR image data
[J].
Michigan microwave canopy scattering model
[J].
Vegetation modeled as a water cloud
[J].
Soil moisture experiment in the Luan River supporting new satellite mission opportunities
[J].
On current limits of soil moisture retrieval from ERS-SAR data
[J].
Mapping soil moisture using Radarsat-2 data and local autocorrelation statistics
[J].
Estimation of bare surface soil moisture with L-band multi-polarization radar measurements
[C]∥
Potential use of Radarsat-2 polarimetric parameters for estimating soil moisture in prairie areas
[C]∥
Soil moisture retrieval from Sentinel-1 time-series data over croplands of Northeastern Thailand
[J].
Join inversion of soil moisture using active and passive microwave data
[J].
利用主被动微波数据联合反演土壤水分
[J].
Potential of L-band radar for retrieval of canopy and subcanopy parameters of boreal forests
[J].
Soil moisture estimation with ERS-1 SAR data in the East-German Loess soil area
[J].
Soil moisture mapping using Sentinel-1 images: Algorithm and preliminary validation
[J].
Microwave-based vegetation descriptors in the parameterization of water cloud model at L-band for soil moisture retrieval over croplands
[J].
Soil moisture retrieval from active spaceborne microwave observations: An evaluation of current techniques
[J].
A method for estimating soil moisture from ERS scatterometer and soil data
[J].
Vegetation modeled as a water cloud
[J].
Soil moisture retrieval over croplands using dual-pol L-band GRD SAR data
[J].
An improved change detection method for soil moisture retrieval using Sentinel-1 and SAMP data
[C]∥
A method for soil moisture estimation in Western Africa based on the ERS scatterometer
[J].
Synergetic use of Sentinel-1 and Sentinel-2 data for soil moisture mapping at 100 m resolution
[J].
The impact of thermo⁃karst lake formation on soil environment of alpine meadow in permafrost regions in the Beiluhe basin of the Tibetan Plateau
[J].
Operational mapping of soil moisture using synthetic aperture radar data:Application to the Touch basin (France)
[J].
Drying-Wetting changes of surface soil moisture and the influencing factors in permafrost regions of the Qinghai-Tibet Plateau,China
[J].
Soil hydrological characteristics of the degrading meadow in permafrost regions in the Beilu River basin
[J].
北麓河流域多年冻土区退化草甸的土壤水文特征分析
[J].
Integration study on oasis soil moisture inversion using ALOS-2 and Landsat-8
[J].DOI:10.13203/j.whugis20220008 [本文引用: 1]
联合ALOS-2和Landsat-8的绿洲土壤水分反演模型研究
[J].
/
〈 |
|
〉 |
