基于Sentinel-1雷达数据的青藏高原地区土壤水分反演研究
Soil Moisture Retrieval in the Tibetan Plateau based on Sentinel-1 Radar Data
通讯作者:
收稿日期: 2021-07-01 修回日期: 2022-10-29
基金资助: |
|
Received: 2021-07-01 Revised: 2022-10-29
作者简介 About authors
方西瑶(1996-),女,甘肃陇南人,硕士研究生,主要从事微波遥感土壤水分反演E⁃mail:
关键词:
Keywords:
本文引用格式
方西瑶, 蒋玲梅, 崔慧珍.
Fang Xiyao, Jiang Lingmei, Cui Huizhen.
1 引 言
青藏高原有“世界屋脊”之称,平均海拔高于4 000 m,雄踞亚洲东部,是中纬度大气环流中庞大的阻碍物 [1]。由于特殊的地理位置和地貌特征,青藏高原在东南亚季风的形成演化和全球气候系统的大气环流中有着重要影响,同时也是我国和东南亚地区生态和气候系统稳定的重要屏障[2]。而作为地表水存储的重要组成部分,土壤水分既是地表—大气能量交换的关键参数之一,也是影响生态稳定和农业生产的重要因素,高空间分辨率的土壤水分有助于加深对水文过程的理解、改进降雨和蒸散发的估算、提升灾害预警能力,在当地和区域的农业管理、土壤侵蚀建模方面也至关重要[3]。因此获取高空间分辨率的土壤水分对青藏高原地区的水文过程研究、生态环境保护和生产生活等方面都具有重要的意义。
目前青藏高原地区的土壤水分获取方法主要分为实地测量、陆面水文过程建模和遥感技术手段。陆面水文过程建模可以获得大面积的深层土壤水分信息,但模型估算土壤水分存在很大的不确定性[4]。地面测量虽然可以进行逐点的观测,获得中小尺度的高精度土壤水分信息,但观测成本高、空间上非常稀疏,在人难以到达的地区无法进行地面观测,不能满足实际应用的需求[5]。遥感技术为土壤水分区域性“面”信息的获取提供了有效途径。随着遥感技术发展,卫星的空间分辨率、重访周期不断提高,遥感手段已成为了大面积监测土壤水分的主流方法。其中微波遥感可以穿透云雾和小雨,不受天气的影响,具有较强的全天候、全天时的工作能力,且微波对地物具有穿透性,对地表土壤水分相对敏感,主/被动微波遥感是土壤水分遥感反演的有效手段。
微波遥感中,被动微波遥感利用土壤的微波发射率进行土壤水分反演[6-7],常用于土壤水分反演的卫星有多通道扫描微波辐射计(Scanning Multichannel Microwave Radiometer, SMMR)、高级微波扫描辐射计/地球观测系统(Advanced Microwave Scanning Radiometer for EOS, AMSR-E)、风云三号卫星(FengYun-3B/C/D, FY-3B/C/D)、土壤湿度与海水盐度卫星(Soil Moisture and Ocean Salinity mission, SMOS)和土壤水分主被动探测卫星(Soil Moisture Active and Passive, SMAP)等。被动微波遥感虽然对土壤水分敏感,时间分辨率高,但其空间分辨率低。目前,国际上已有的被动微波土壤水分产品如AMSR-E/2、FY-3B/C/D、SMOS、SMAP等卫星土壤水分产品空间分辨率主要为25—40 km,无法满足水文气候建模和农业、水资源管理的应用需求。
主动微波遥感向目标发射电磁波,利用回波信息分析地物的特征,土壤水分会直接影响地表土壤的介电常数,介电常数与后向散射系数间有着本质关系,由此可通过后向散射系数反演出土壤水分信息。相较之下,主动微波遥感既可以全天时、全天侯进行工作,又能较大范围地获取更精细尺度的地表信息,在区域精细尺度土壤水分监测方面具有潜力,越来越多的学者尝试利用主动微波遥感进行土壤水分反演研究。
由于青藏高原特殊的气候环境和地理位置,一些研究针对该地区利用主动微波遥感开展了土壤水分反演工作。如武胜利等[8]基于积分方程模型(Integrated Equation Model, IEM)和一阶植被散射模型,利用不同时相的TRMM/PR数据获取了土壤水分变化量信息,该工作结果表明ku波段范围内的TRMM/PR数据在一定程度上也可以反映青藏高原地区的土壤水分变化情况,不同站点反演值的均方根误差小于0.07 m³/m³。王建明[9]基于简化的地表散射模型,利用ERS-l/2散射计的前、中天线以不同入射角观测目标的特点,用单时相的不同入射角观测数据发展了估算裸露地表土壤水分绝对值的算法,反演结果与实测值的相关系数为0.6—0.7,标准偏差约为0.04 m³/m³。在此基础上,结合水云模型(Water-Cloud Model, WCM)校正了植被的影响,估算了青藏高原地表土壤水分的时空分布,土壤水分估值的季节变化与研究区的降雨分布/水气输送路线也较为一致。上述工作都是基于粗尺度的雷达和散射计估算土壤水分,空间分辨率为公里级。随着主动雷达技术的发展,一些学者在青藏高原地区开展了更精细尺度的土壤水分研究。如Xu等[10]用简单的线性表达式表征了裸土地表散射和土壤水分之间的关系,并用水云模型描述了植被的影响,基于简化的水云模型,利用中分辨率成像光谱仪(Moderate Resolution Imaging Spectroradiometer, MODIS)的反射率数据和Sentinel-1观测数据估算了青藏高原中部地区的土壤水分,空间分辨率为500 m。土壤水分估算值与实测值之间的均方根误差小于0.07 m³/m³,但该线性模型属于经验模型,算法缺乏一定的普适性。Wang等[11]基于半经验模型Oh模型,利用ALOS PALSAR-2 ScanSAR成像模式数据估算了优化后的的裸土地表粗糙度参数,并基于改进的水云模型,估算了植被的有效散射反照率和地表土壤水分,土壤水分反演值与地面测量值之间有较高的一致性。需要注意的是,虽然已有不少研究利用主动微波在青藏高原地区估算了土壤水分,但目前的研究工作中空间分辨率多在百米或公里级以上。
综上所述,获取青藏高原地区高空间分辨率的土壤水分信息在科学研究和生产生活方面都具有重要价值。因此,为了满足区域尺度研究的需求,获取青藏高原地区高空间分辨率的土壤水分,本文基于Sentinel-1 合成孔径雷达(Synthetic Aperture Radar, SAR)数据,结合微波散射模型,利用神经网络方法发展青藏高原地区高空间分辨率(20 m)土壤水分反演算法,以期为后续青藏高原地区的科学研究提供数据支撑。
2 研究区与数据预处理
2.1 研究区概况
图1
图1
研究区地表覆盖示意图
(地表覆盖数据来源:Global Map V.1 Land Cover (GLCNMO))
Fig.1
Land cover of the study area
2.2 Sentinel-1数据及预处理
Sentinel-1卫星是欧空局哥白尼计划中的地球观测卫星系统,其数据具有以下特点:①时间分辨率高,单颗卫星的重访周期为12 d,两颗卫星组成星座,其重复周期可缩短至6 d;②空间分辨率高,最高可达5 m;③成像模式多,有单极化和多极化等不同的极化方式,成像系统采用4种成像模式,可精确确定卫星的位置和姿态角;④相较于其他的雷达数据,Sentinel-1数据可以免费获取,数据更新较快[13]。因此研究选取Sentinel-1 SAR数据进行青藏高原地区高空间分辨土壤水分反演的研究。
Sentinel-1卫星雷达的主要工作方式是干涉宽幅(Interferometric Wide swath, IW)模式和波浪模式(Wave Model),陆地上的默认模式为干涉宽幅模式,波浪模式主要用于获取海洋参数,另外还有条带成像(Stripmap Model, SM)模式和超宽幅(Extra Wide swath, EW)两个附加模式[14]。Sentinel-1雷达传感器的每一种模式都可以生成4种产品:Level-0级产品、Level-1级SLC产品、Level-1级GRD(Ground Range Detected)产品和Level-2级OCN(Ocean)产品,对于波浪模式,不提供Level-0和Level-级产品。
采用Sentinel-1数据为干涉宽幅模式下的Level-1级GRD数据,空间分辨率可达20m,数据由Copernicus Open Access Hub (https:∥scihub.copernicus.eu/)提供。根据研究区气候条件下载相应地区的2017、2018 年地表非冻结期间(5月—10月)的影像数据,数据的详细信息见表1。
表1 雷达数据信息
Table1
站点 | 影像时间(年月日) | 极化方式 |
---|---|---|
曲麻莱 | 20170603、20170615、20170627、20170709、20170721、20170802、20170814、20170826、20170907、20170919、20171013、20171025、20180517、20180529、20180610、20180704、20180716、20180728、20180821、20180914、20181008、20181020 | VV,VH |
查拉坪、玛沁 | 20170505、20170610、20170622、20170704、20170716、20170728、20170809、20170821、20170902、20170914、20171008、20171020、20171101、20180512、20180524、20180629、20180711、20180723、20180804、20181015 | VV,VH |
那曲、 纳木错 | 20170501、20170513、20170525、20170606、20170618、20170630、20170712、20170724、20170805、20170817、20170829、20170910、20170922、20171004、20171016、20171028、20180508、20180601、20180613、20180625、20180707、20180824、20180905、20180917、20180929 | VV,VH |
研究主要利用ESA提供的SNAP软件对Sentinel-1雷达数据进行了预处理,主要包括多视处理、辐射定标、地形辐射校正等。多视处理中距离向与方位向的视数设置为2,经多视处理,影像像元的空间分辨率约为20 m。地形校正过程中的DEM数据来自于SRTM (Shuttle Radar Topography Mission)DEM产品。最后根据地面站点的地理位置信息提取相应的VV极化和VH极化下的Sentinel-1后向散射系数和局部入射角。
2.3 MODIS植被指数数据
由于植被会影响植被覆盖地表的雷达观测信号,土壤水分反演时需要考虑植被的影响。目前归一化植被指数(Normalized Difference Vegetation Index, NDVI)已被广泛应用于表征植被状况,并用于量化植被对土壤表面后向散射的影响[15]。尽管有研究指出NDVI在植被高覆盖区容易饱和,但本文的研究区植被稀疏,全年的NDVI最大值不超过0.6,饱和效应可以忽略。因此,利用NDVI表示植被的生长状况,以考虑土壤水分反演中植被的影响。由于高空间分辨率的光学数据Landsat 8、Sentinel-2在研究区时间分辨率较低,受云影响严重,研究时段内的完整性较差。此外,研究挑选了2018年5月25日云量较少的Sentinel-2数据,与同一天的MODIS NDVI数据做比较,发现MODIS与Sentinel-2位于站点所在像元的NDVI差值为0.012 7,差异非常小,几乎可以忽略。因此研究选用了MODIS双星的 MOD13Q1、MYD13Q1产品的 NDVI数据来校正雷达土壤水分反演中的植被影响。该数据是16 d日合成的连续、无云产品,两数据集结合后时间分辨率为8 d。虽然MODIS的NDVI产品空间分辨率为250 m,与处理后的雷达数据存在一定的空间差异,但本文的研究区地势相对平坦,地表覆盖类型均为稀疏低矮的高寒草原,均一性较好。
根据研究区站点的经纬度信息,研究下载了与雷达数据相近时间段的MOD13Q1、MYD13Q1影像,MODIS数据具体时间跨度为2017年4月23日至11月1日、2018年4月23日至11月1日。由于研究获取的雷达数据影像时间与 MODIS 数据影像时间并不完全一致,所以在选取雷达影像数据对应的NDVI值时以时间最近为原则。最终,大部分站点的每组雷达数据与MODIS数据时间差不超过1 d,最大不超过4 d。
对MODIS数据按经纬度0.002 5°间隔重采样,并提取出观测站点所在像元的NDVI。受云、大气等影响,NDVI 时序数据中有不规则的锯齿状波动,通过Savitzky-Golay滤波对NDVI时序数据的不规则波动进行了平滑。
2.4 地面数据
地面站点的实测数据来自“青藏高原积雪深度连续自动观测网络数据集”[12],该观测网络每个站点都配有SR50A 超声波测距传感器、Model 109 温度传感器和 Hydra Probe II 土壤多参数传感器,每10 min一次获取2 m空气温度、积雪深度、土壤温度、5 cm深度的土壤水分(体积含水量)、土壤电导率和土壤复介电常数。经分析,数据集中土壤体积含水量精度约为3%,土壤温度精度约为0.6 ℃。
提取研究区站点数据集中与雷达观测时间最接近的实测土壤水分和地表温度数据,时间差小于5 min。由于站点分布较为稀疏,雷达数据量较大,处理较为繁琐,因此,研究选取了雷达影像上较为集中的5个站点(曲麻莱、查拉坪、那曲、纳木错、玛沁)的观测数据进行卫星反演验证与分析。剔除地表温度为负值的数据后获得了5个观测站点2017年和2018年地表非冻结期间5月—10月的101组卫星、实测数据。由于地面测量数据和雷达卫星影像的缺失,各观测站数据量不一致,曲麻莱观测站21组(2017年11组,2018年10组),查拉坪观测站18组(2017年12组,2018年6组),那曲观测站24组(2017年15组,2018年9组),纳木错观测站25组(2017年16组,2018年9组),玛沁观测站缺失2018年的地面测量数据,只有2017年间的13组数据可用。文中利用2017年的数据集调整算法阈值和训练,2018年的数据用于测试与检验算法。
3 方 法
3.1 理论模型
水云模型中假设:①植被层由许多大小、形状相同的散射微粒组成;②植被和土壤间的多次散射项可以忽略。最终传感器接收到的后向散射系数可以分为两部分:植被的直接散射项和经植被两次衰减的土壤后向散射项,具体公式如下:
其中:
3.2 反演方法
由于机器学习方法依赖于训练数据,仅依靠野外测量和卫星观测的样本点难以满足典型地表的范围,因此本文结合微波散射模型模拟数据集来扩充训练数据集,以增强训练数据的代表性和反演过程的物理性。基于观测与模型模拟数据集,借助于人工神经网络方法来训练雷达观测后向散射系数与土壤水分之间的非线性关系,进一步反演地表土壤水分。
研究思路为:首先将裸露地表与植被覆盖地表的电磁波散射模型耦合,即将AIEM模型与Oh模型结合,用于描述裸露地表不同极化的后向散射系数,采用零阶的植被散射模型—水云模型描述植被的散射特征。由于水云模型中未知参数A、B的值取决于植被的类型和结构,因此研究利用地面观测站点的土壤水分、站点所在像元的植被含水量、雷达观测数据来校正水云模型中的A、B参数,获得水云模型中A、B参数的最优值。然后基于参数优化后的水云模型建立不同雷达入射角、土壤水分、植被覆盖度、粗糙度情况下的地表散射模拟数据。将该模拟数据和部分实测数据作为神经网络算法的训练集进行训练测试以进一步反演土壤水分。
4 结果和分析
4.1 水云模型参数优化
青藏高原的土壤主要是壤质土或砂质土[9,36]。土壤介电常数模型Dobson模型中统一取土壤中砂土含量为35%,粘土含量为20%,参考世界土壤数据库的土壤质地类型数据[37],土壤容重设置为1.61 g/cm3。根据Sentinel-1过境时的土壤温度实测均值,土壤温度设置为10 ℃,频率取Sentinel-1的中心频率5.405 GHz。设置均方根高度与相关长度的初始值分别为0.8 cm和8 cm,这两个参数迭代区间分别为0.1—3 cm,1—20 cm。利用2017年期间的67组卫星观测与实测数据,构建模型模拟后向散射系数与雷达观测后向散射系数间的代价函数,基于代价函数最小化准则,用Nelder-Mead优化方法[38]分别求解VV极化和VH极化下参数A、B的最优解,同时得到一组优化的地表粗糙度参数。最终得到研究区最优的均方根高度与相关长度均值为0.4 cm、5 cm,VV极化下A、B参数最优解分别为0.019和0.183,VH极化下A、B参数最优解分别为0.003和0.173。
利用2017年的实测数据、卫星观测数据对参数优化后模型模拟能力进行验证。结果如图2所示,优化后模型模拟的站点雷达后向散射系数与Sentinel-1卫星观测值在时间序列上的比较。经统计,在VV极化下观测值与模拟值相关系数为0.98,相比于VV极化,VH极化对土壤水分的响应相对较差,对植被、粗糙度的敏感性更高,因此拟合度稍差,但相关系数达到0.87,两种极化下的均方根误差都比较小,总体而言,参数优化后的模型能较好地表征研究区站点的地表散射特性,可进一步用于土壤水分反演。
图2
图2
参数优化后的模拟数据与Sentinel-1卫星数据的比较
Fig.2
Comparison of optimized simulated data with Sentinel-1 satellite data
4.2 建立模拟数据集
为了增强反演方法的物理机理,文中利用4.1节中参数优化后的微波散射模型在有效的地表参数范围内进行模拟,建立模拟数据集用以神经网络的训练与测试。其中土壤水分、雷达入射角、植被含水量的相应范围根据实测数据和观测数据设定,地表相关长度和均方高度的取值范围根据4.1节中的优化结果设置。
同时为了使模拟数据库具有广泛的适用性,可以表征更广泛地表、卫星参数条件下的后向散射特征,模型输入参数应取较宽的范围。文中模拟所用的具体参数及取值范围如表2所示。
表2 模拟数据集参数设置范围
Table 2
模型参数 | 最小值 | 最大值 | 步长 | 单位 |
---|---|---|---|---|
土壤水分 | 0.02 | 0.5 | 0.01 | m³/ m³ |
入射角 | 25 | 50 | 5 | degree |
均方根高度 | 0.1 | 3 | 0.1 | cm |
表面相关长度 | 1 | 20 | 1 | cm |
植被含水量 | 0 | 0.2 | 0.02 | kg/m2 |
频率 | 5.405GHz |
图3
图3
实测数据与模拟数据的关系
Fig.3
The relationship between measured data and simulated data
4.3 土壤水分反演及验证
利用BP神经网络模型反演土壤水分,模型由输入层、隐藏层和输出层组成。输入层的神经元接受外界信息,通过中间隐藏层的各神经元,将信息传递给输出层的各神经元,进行信息的正向传递,若输入层的实际输出与期望的输出不符,则进行误差反向传递,修正各层神经元的权值,直到输出的误差减小到可接受的程度或达到设定的学习次数。本文最终目的是根据卫星观测值得到土壤水分信息,因此对神经网络模型输入端参数设为:VH 极化后向散射系数、VV 极化后向散射系数、植被含水量和雷达入射角;输出参数设为土壤水分。
对已经确定结构的模型进行多次训练,用2017年的所有实测数据对每次的训练结果反演验证。以反演结果的相关性系数、均方根误差、偏差作为依据,选择反演效果最好的训练结果。根据训练结果,2017年的土壤水分实测值与模拟值间相关系数约为0.784,均方根误差较小,为0.064 m³/m³,土壤水分估算偏差约为-0.009 m³/m³。
利用2018年的数据对算法进行检验,土壤水分反演值与实测值的相关系数为0.82,均方根误差为0.052 m³/m³,该结果表明利用Sentinel-1数据在青藏高原地区反演高精度、高空间分辨率的土壤水分有较大潜力。从相关系数等指标看,2018年的验证结果要优于2017年,这是因为反演算法在玛沁站点的效果较差(2017年的验证结果中R仅为0.526),2018年的验证数据缺少玛沁站点的观测,因此验证的总体精度得到了提升。
图4
图4
2018年土壤水分反演值与实测土壤水分的比较
Fig.4
Comparison between estimated and in situ soil moisture in 2018
图5展示了各个站点反演得到的土壤水分与实测土壤水分的时间序列图。在不同站点、不同时间段,反演算法有不同的表现。在那曲观测站除2017年5月外,反演土壤水分与实测土壤水分变化趋势具有较高的相关性。纳木错站点在2017年7月至10月上旬、2018年6月至10月的反演精度较高,土壤水分反演值与实测值的变化趋势也较为一致。曲麻莱站点在整个研究期内反演误差都较小,反演精度较高,反演值能捕捉到实测土壤水分的变化。玛沁站点由于地面实测数据缺失,只有2017年时间段的验证,但在7月至10月中旬时间段内的土壤水分反演值和实测值有很好的一致性。反演算法在查拉坪站点的误差较大,但反演值与实测值的变化趋势相似,相关系数大于0.5。
图5
图5
土壤水分实测值与土壤水分反演值在时间序列上的变化
Fig 5
Time series of soil moisture measurements and soil moisture retrievals
表3统计了各观测站的反演精度(土壤水分实测值与土壤水分反演值的相关系数R、均方根误差RMSE、偏差Bias)以及实测土壤水分与卫星观测雷达后向散射系数的相关性(土壤水分分别与VV、VH极化后向散射系数的相关系数Rvv、Rvh)。我们发现对于各观测站,其反演精度在一定程度上受该站点实测土壤水分和雷达后向散射系数的相关性影响。那曲、曲麻莱站点土壤水分与雷达VV极化的观测值之间有较高的相关性,其站点土壤水分的反演精度也相对较好。均方根误差大、反演精度较差的查拉坪、纳木错站点,其Rvv、Rvh也较低。同时图5和表3的结果表明,在研究期内大部分站点的反演结果都存在低估现象,这可能与地表粗糙度参数的误差有关。此外随着土壤水分变大,雷达后向散射系数对土壤水分变化的敏感性逐渐降低[43],这也会使算法在土壤水分较大站点(如查拉坪站)的反演值偏低。土壤水分反演值低估的具体原因,还需在后续工作中收集更多的实测地面参数进一步探究。
表3 各站点土壤水分反演精度及土壤水分测量值与雷达后向散射系数相关性
Table 3
站点 | R | RMSE | BIAS | Rvv | Rvh |
---|---|---|---|---|---|
那曲 | 0.828 | 0.038 | -0.021 | 0.818 | 0.632 |
纳木错 | 0.717 | 0.079 | 0.058 | 0.655 | 0.183 |
查拉坪 | 0.662 | 0.074 | -0.056 | 0.397 | 0.519 |
玛沁 | 0.526 | 0.055 | -0.034 | 0.613 | 0.360 |
曲麻莱 | 0.726 | 0.040 | -0.016 | 0.804 | 0.704 |
4.4 区域土壤水分反演
上述的基于神经网络的土壤水分反演算法,经地面观测检验效果较好。因此研究利用该算法开展了区域土壤水分的应用。考虑到算法的适用性,选择以观测站为中心,对观测站周围1 km范围的区域进行土壤水分反演,以探究基于Sentinel-1的区域土壤水分反演合理性与可行性。图6—7分别是2017年7月12日纳木错和那曲区域、2017年7月16日玛沁、查拉坪区域、2017年6月3日曲麻莱区域空间分辨率为20 m的土壤水分、VV极化雷达后向散射系数空间分布。从结果来看,各个地区土壤水分的空间变化与VV极化的雷达后向散射系数的空间分布响应较好,土壤水分呈现出了较大的空间异质性,高空间分辨率的土壤水分空间分布图显示了较多的空间细节信息。该方法在获取区域高空间分辨率的土壤水分研究中具有一定的潜力。
图6
图6
不同地区1 km×1 km范围内的土壤水分(m³/ m³)
Fig.6
Spatial distribution of soil moisture (m³/ m³) within 1 km×1 km in different regions
图7
图7
不同地区1 km×1 km范围内的VV极化后向散射系数(dB)
Fig.7
Spatial distribution of VV polarization backscattering coefficient (dB) within 1 km×1 km in different regions
5 结 语
青藏高原地理位置、地貌特征特殊,高空间分辨率的土壤水分对该地区在全球水循环与能量循环等研究中有重要作用。本文基于Sentinel-1数据、MODIS归一化植被指数和土壤水分地面观测数据,结合水云模型,AIEM模型和Oh模型,利用BP神经网络算法反演了青藏高原部分地区的高空间分辨率土壤水分,并对其进行了验证。主要结论如下:
(1)微波散射模型模拟数据的加入增强了BP神经网络算法在土壤水分反演过程中的物理性,土壤水分反演值与地面实测值具有较高的相关性,其相关系数为0.784—0.82,总体均方根误差小于0.064 m3/m3。时间序列上,土壤水分反演值能够捕捉到实测值的变化趋势。
(2)C波段的Sentinel-1雷达数据在获取青藏高原地区的高空间分辨率土壤水分方面具有一定的潜力,高空间分辨率(20 m)的Sentinel-1土壤水分反演值能够在空间分布上能够表征更多的土壤水分细节信息。
此外,本文的研究中仍存在一些不足,如:青藏高原中东部地区有较为丰富的土壤有机碳,本文所用的土壤介电常数模型未考虑有机碳的影响;由于缺乏实测的植被含水量数据,模型中估算植被含水量方法在研究区并未得到有效的验证;地表粗糙度是影响土壤水分反演的主要因素之一,文中缺少野外的粗糙度实测数据,这对模型优化具有一定影响,可能会进一步影响土壤水分的反演精度。目前论文中使用的观测站点稀疏、数量较少,优化的水云模型参数在大区域内的普适性存在一定限制。此外,MODIS数据空间分辨率粗,在地形、植被状况变化剧烈的区域,NDVI与雷达数据存在一定的空间尺度差异。下一步工作还需考虑上述因素,进一步优化大尺度下高空间分辨率土壤水分的反演算法,以期提高高空间分辨率土壤水分反演精度。
参考文献
A discussion on the effect of Qinghai Tibet Plateau uplifting on geological and natural environment in China
[J].
青藏高原隆起及其对中国地质自然环境影响的探讨
[J]
Role of the Tibetan Plateau thermal forcing in the summer climate patterns over subtropical Asia
[J].
A roadmap for high-resolution satellite soil moisture applications - confronting product characteristics with user requirements
[J].
Validation of soil moisture simulation in landsurface parameterisation schemes with HAPEX data
[J].
A semi-empirical approach for surface soil water content estimation from radar data without a-priori information on surface roughness
[J].
Soil moisture retrieval from space: the Soil Moisture and Ocean Salinity (SMOS) mission
[J].
Recent advances on algorithms and applications of soil moisture retrieval from microwave remote sensing
[J].
The application of TRMM/PR to the study of soil moisture change in Tibetan Plateau
[J]
TRMM/PR在青藏高原土壤水分变化研究中的应用
[J].
Study on the method of estimating soil moisture in Qinghai-Tibet Plateau based on ERS scatterometer data
[D].
基于ERS散射计数据的青藏高原土壤水分估算方法研究
[D].
Monitoring surface soil moisture content over the vegetated area by integrating optical and SAR satellite observations in the permafrost region of Tibetan Plateau
[J].
Soil moisture retrievals using ALOS2-ScanSAR and MODIS synergy over Tibetan Plateau
[J].
An automatic measurement dataset of snow depth on the Tibetan Plateau (2015-2016)
[J]
2015—2016年青藏高原积雪深度自动观测数据集
[J].
Sentinel-1 Satellite Overview
[J].
Sentinel-1卫星综述
[J].
Application of Sentinel-1 satellite and data products
[J].
Sentinel-1卫星数据产品应用探讨
[J].
Assessment of different vegetation parameters for parameterizing the coupled water cloud model and advanced integral equation model for soil moisture retrieval using time series Sentinel-1A data
[J].
Emission of rough surfaces calculated by the integral equation method with comparison to three-dimensional moment method simulations
[J].
Backscattering from a randomly rough dielectric surface
[J].
An update of AIEM model with multiple scattering of rough surface
[C]∥
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].
Soil moisture inversion research in bare region based on L-band radar data
[J].
基于L波段的裸土区土壤水分微波遥感反演研究
[J].
Time-series approach to estimate the soil moisture of a subsidence area by using dual polarimetric radar data
[J].
Vegetation modeled as a water cloud
[J].
Michigan microwave canopy scattering model
[J].
Monitoring of wheat crops using the backscattering coefficient and the interferometric coherence derived from Sentinel-1 in semi-arid areas
[J].
Retrieval of high-resolution soil moisture through combination of Sentinel-1 and Sentinel-2 data
[J].
Combined use of sentinel-1 SAR and Landsat sensors products for residual soil moisture retrieval over agricultural fields in the upper Blue Nile Basin,Ethiopia
[J].
Retrieval of both soil moisture and texture using TerraSAR-X images
[J].
Dense temporal series of C- and L-band SAR data for soil moisture retrieval over agricultural crops
[J].
An improved approach for soil moisture estimation in gully fields of the loess plateau using Sentinel-1A radar images
[J].
Soil moisture evaluation using machine learning yechniques on Synthetic Aperture Radar (SAR) and land surface model
[C]∥
Biomass assessment of agricultural crops using multi-temporal dual-polarimetric terraSAR-X data
[J].
A comparison of algorithms for retrieving soil moisture from ENVISAT/ASAR images
[J].
The potential of ALOS-2 and Sentinel-1 Radar data for soil moisture retrieval with high spatial resolution over agroforestry Areas, China
[J].
A soil particle-size distribution dataset for regional land and climate modelling in China
[J].
The harmonized world soil database
[C]∥
A simplex method for function minimization
[J].
A prototype ann based algorithm for the soil moisture retrieval from l- band in view of the incoming SMAP mission
[C]∥
Sequential ensembles tolerant to Synthetic Aperture Radar (SAR) soil moisture retrieval errors
[J].
Error in radar-derived soil moisture due to roughness parameterization: An analysis based on synthetical surface profiles
[J].
Using SAR-derived vegetation descriptors in a water cloud model to improve soil moisture retrieval
[J].
/
〈 |
|
〉 |
