数据可用性声明:本文关联实体数据集已在《全球变化数据仓储电子杂志(中英文)》出版,可获取:https://doi.org/10.3974/geodb.2020.08.03.V1
1 前言
土壤水分一般指单位体积土壤中所含液态水的体积,是全球气候变化、陆地表面水文 过程、碳循环的关键物理量之一[1, 4]。土壤水分通过改变土壤热容量,控制地表蒸发和植被蒸腾,影响植被生长[5, 7]。因此,获取土壤水分含量数据不仅对于分析陆地生态系统演替、碳氮水循环过程具有关键科学意义,而且对于开展旱涝灾害预警、进行农作物估产具有重要实际应用价值[8, 11]。
卫星遥感技术是获取全球尺度、连续时间序列陆地表层土壤水分数据的重要手段,海量丰富的卫星土壤水分产品为开展全球气候演化分析提供了前所未有的机遇[12, 13]。但受到卫星扫描间隙、射频干扰、植被光学厚度、冰冻季节等影响,基于卫星遥感反演的土壤水分产品存在大量的空值区域[14, 16]。不同微波波段(如C波段、X波段、K波段、Ka波段、L波段)对地表土壤水分的敏感度也不尽相同[17]。为了提高卫星土壤水分产品的时空序列完整性和数据质量,欧洲太空局(European Space Agency,ESA)在2010年起致力于融合多源卫星遥感数据,形成一套迄今为止时间序列最长的多源卫星数据融合全球地表土壤水分产品(Essential Climate Variable Soil Moisture,ECV SM)[18, 20]。ECV土壤水分产品的时间序列长度、空间序列完整性以及数据精度相较于单波段卫星土壤水分产品有了显著提升,但与同化产品相比其空间覆盖率仍有较大提升空间。因此,融入新的卫星土壤水分产品是提升ECV土壤水分产品空间完整性和数据质量的有效手段。本数据集将2015年问世的L波段SMAP(Soil Moisture Active Passive)土壤水分产品[21]与ECV土壤水分产品融合[22],经过重投影、重采样和插补[23, 26],最终得到了2015年3月31日至2019年12月31日全球地表土壤水分0.25 分辨率逐日数据集。
2 数据集元数据简介
《全球土壤水分产品融合数据集(2015–2019)》[27]的名称、作者、地理区域、数据年代、时间分辨率、空间分辨率、数据集组成、数据出版与共享服务平台、数据共享政策等信息见表1。
| 条目 | 描述 | ||||
|---|---|---|---|---|---|
| 数据集名称 | 全球土壤水分产品融合数据集(2015–2019) | ||||
| 数据集短名 | Global_SM | ||||
| 作者信息 | 刘杨晓月 ABB-3889-2020,广东省科学院广州地理研究所,lyxy@lreis.ac.cn | ||||
| 地理区域 | 全球 | 数据年代 2015–2019 时间分辨率逐日 | |||
| 数据格式 | .tif | 数据量 | 6.71 GB | 空间分辨率 | 0.25×0.25 |
| 数据集组成 | 包括1,737个数据文件,自2015年3月31日至2019年12月31日期间每天有一个数据文件,以SM-yyyymmdd.tif 形式命名,如SM-20160101.tif 是2016年1月1日的全球土壤水分融合数据 | ||||
| 基金项目 | 博士后创新人才支持计划(BX20200100);国家对地观测科学数据中心(NODAOP2020002);南方海洋科学与工程广东省实验室(广州)人才团队引进重大专项(GML2019ZD0301) | ||||
| 数据计算环境 | Python 2.7 | ||||
| 出版与共享服务平台 | 全球变化科学研究数据出版系统http://www.geodoi.ac.cn | ||||
| 地址 | 北京市朝阳区大屯路甲11号100101,中国科学院地理科学与资源研究所 | ||||
| 数据共享政策 | 全球变化科学研究数据出版系统的“数据”包括元数据(中英文)、通过《全球变化数据仓储电子杂志(中英文)》发表的实体数据和通过《全球变化数据学报(中英文)》发表的数据论文。其共享政策如下:(1)“数据”以最便利的方式通过互联网系统免费向全社会开放,用户免费浏览、免费下载;(2)最终用户使用“数据”需要按照引用格式在参考文献或适当的位置标注数据来源;(3)增值服务用户或以任何形式散发和传播(包括通过计算机服务器)“数据”的用户需要与《全球变化数据学报(中英文)》编辑部签署书面协议,获得许可;(4)摘取“数据”中的部分记录创作新数据的作者需要遵循10%引用原则,即从本数据集中摘取的数据记录少于新数据集总记录量的10%,同时需要对摘取的数据记录标注数据来源[28] | ||||
| 数据和论文检索系统 | DOI,DCI,CSCD,WDS/ISC,GEOSS,China GEOSS,Crossref | ||||
3 数据研发方法
3.1 数据来源
本数据集中所采用的ECV土壤水分产品为基于多源主动微波(ERS-1、ERS-2、MetOp-A、ASCAT)和被动微波产品(SMMR、SSM/I、TMI、AMSR-E、AMSR-2、Windsat、SMOS)融合的全球逐日土壤水分产品[21]。本数据集采用的SMAP土壤水分产品为基于L波段被动微波辐射计反演的2015–2019年全球逐日土壤水分产品,空间分辨率为36 km[22]。大量已有研究表明L波段对地表土壤水分的敏感度优于其他微波波段,对地面的穿透深度与土壤水分地面监测传感器的埋深最为接近。研发机构从SMAP卫星传感器设计到反演算法构建均强化了对地面人为射频的抗干扰能力,经验证表明,SMAP土壤水分产品的数值精度和时空序列拟合度均优于构成ECV的每种卫星土壤水分产品[23, 29, 31]。
数据集评价过程中使用的土壤水分地面实测数据源于国际土壤水分网络[32],本研究选用的8个土壤水分地面实测网络的基本属性信息如表2所示。
| 名称 | 所在国家 | 站点数量 | 区域气候 | 主要土地覆被类型 |
|---|---|---|---|---|
| REMEDHUS | 西班牙 | 20 | 温带海洋气候 | 耕地、灌丛 |
| FR_Aqui | 法国 | 4 | 地中海气候 | 耕地、林地 |
| FMI | 瑞典 | 20 | 亚寒带针叶林气候 | 稀树草原 |
| HOBE | 丹麦 | 27 | 温带海洋气候 | 耕地、林地 |
| BIEBRZA_S | 波兰 | 18 | 温带大陆性气候 | 草地、湿地 |
| TERENO | 德国 | 5 | 温带海洋气候 | 耕地、林地 |
| RMSN | 罗马尼亚 | 19 | 温带大陆性气候 | 耕地、林地 |
| SMOSMANIA | 法国 | 21 | 地中海气候 | 多类型混合 |
3.2 研究方法
本研究旨在提高ECV土壤水分产品数据完整性同时促进数据精度稳中有升。首先,利用土壤水分地面实测数据对SMAP和ECV土壤水分产品进行质量评价。对实测土壤水分数据采取质量控制措施,实测值一般以小时为单位测定土壤含水量,仅当一天内有效数据超过12小时才认定为有效数据,计算加权平均值得到各网络各站点日尺度土壤水分实测值。选取拟合优度(R)、偏差(Bias)和无偏均方根误差(ubRMSE)作为误差参数来系统验证ECV、SMAP数据精度。结果表明SMAP升降轨土壤水分产品精度显著优于ECV数据,可作为ECV插补的优质数据源。第二,为提高有效数据百分比,计算SMAP升轨与降轨数据的加权平均得到日尺度SMAP土壤水分产品。对SMAP进行投影转换和空间重采样,使其与ECV空间属性一致。第三,基于Python读取并遍历ECV日尺度数据,检索空值区域。第四,采用逐日数据逐行逐列遍历的方式,使用与ECV空间属性一致的SMAP日尺度土壤水分产品对ECV进行插补,得到融合的全球土壤水分产品,并对该产品的空间完整性和数据精度进行评价验证。具体数据研制流程如图1所示:
4 数据结果与验证
4.1 数据集组成
数据集包含1,737个文件,格式为.tif,时间范围从2015年3月31日到2019年12月31日,空间范围覆盖全球,时间分辨率为日尺度,空间分辨率为0.25°×0.25°(约25 km× 25 km)。数据单位为m3m–3,即单位体积土壤的含水量,取值范围是[0, 1]。文件命名方式为SM-yyyymmdd.tif。
4.2 数据结果
图2为ECV土壤水分产品与融合后的土壤水分产品对比,由图可知,土壤水分的时空分布特征与区域季相节律周期高度吻合,主要取值范围是0–0.5 m3m–3。以2016年1月1日、4月1日、7月1日、10月1日的数据作为融合产品在冬、春、夏、秋四季的样例,比较可得融合后的土壤水分产品空间覆盖完整性较原始ECV土壤水分产品有非常显著的改善。冬春季节高纬度地区地表温度恒定在0 ℃以下,地表为冻土状态,微波仅能测定土壤中液态水的含量,故表现为空值。格陵兰岛和南极大陆终年被冰雪覆盖,同样表现为空值。
4.3 数据结果验证
4.3.1 完整性评价
由图3可知,融合后土壤水分产品的全球地表覆盖率较ECV土壤水分产品显著提升约20%,填补了ECV在亚马孙热带雨林、刚果盆地热带雨林地区的土壤水分数据空白,为进行土壤水分时空尺度连续分析提供便利。同时,可以观察到土壤水分高覆盖率区域主要分布在60°S–60°N的中低海拔、非植被密集覆盖之处。高寒气候区冬季寒冷漫长,微波信号无法测得液态水回波信号;此外,微波信号难以穿透植被密集覆盖区(>5 kgm–2)。
4.3.2 精度验证
利用取自8个土壤水分网络共计134个站点实测数据,对ECV和SMAP升(A)、降(D)轨土壤水分产品开展验证评价。结果如图4–6所示,图中自上而下的水平线分别代表最大值、上四分位、中位数、下四分位和最小值,虚线是均值,点代表异常值。SMAP土壤水分产品拟合优度和偏差的评价结果优于ECV,无偏均方根误差的评价结果与ECV相当,充分表明SMAP土壤水分产品的质量可靠性,可以作为有效的融合数据来改善ECV土壤水分产品的完整性。
利用实测数据对土壤水分融合结果进行精度验证,验证结果如表3所示,融合土壤水分产品与ECV土壤水分产品精度水平相当[33],在REMEDHUS、FR_Aqui、RSMN和SMOSMANIA的精度优于ECV,能够有效捕捉土壤水分时空序列演化趋势,精确拟合土壤水分地面实测值[34]。总体来说,土壤水分融合结果可以准确反映地面实测数据的分布情况及数值变化。
为进一步分析融合土壤水分产品取值分布与地面实测数据的拟合情况,本研究绘制了实测数据、ECV、土壤水分融合产品的概率分布函数曲线(Probability Distribution Function,PDF)。如图7所示,三类土壤水分数据的取值均符合正态分布,地面实测数据(黑线)的取值离散化较为显著,相比而言,ECV土壤水分产品(红线)取值聚集现象明显,而土壤水分融合产品(蓝线)能够有效调整ECV数据取值过于密集分布情况,与实测值曲线更为接近。
| 土壤水分实测网络 | 拟合优度R | 偏差Bias | 无偏均方根误差ubRMSE | 拟合优度R | 偏差Bias | 无偏均方根误差ubRMSE |
|---|---|---|---|---|---|---|
| REMEDHUS | 0.75 | 0.09 | 0.05 | 0.77 | 0.08 | 0.04 |
| FR_Aqui | 0.76 | 0.11 | 0.04 | 0.77 | 0.11 | 0.04 |
| FMI | 0.11 | ‒0.08 | 0.07 | ‒0.04 | 0.00 | 0.06 |
| BIEBRZA_S-1 | 0.60 | ‒0.34 | 0.14 | 0.58 | ‒0.35 | 0.13 |
| TERENO | 0.67 | 0.02 | 0.06 | 0.57 | 0.02 | 0.06 |
| RSMN | 0.56 | 0.11 | 0.05 | 0.58 | 0.10 | 0.05 |
| SMOSMANIA | 0.62 | 0.08 | 0.10 | 0.68 | 0.07 | 0.06 |
5 讨论和总结
为保证数据质量,在土壤水分融合前后利用分布在8个实测网络的134个站点实测数据对ECV、SMAP升降轨以及融合后的土壤水分产品进行了验证,结果表明本研究获得的土壤水分融合结果能够有效拟合地面实测数据的取值和时空序列演化趋势。本研究的主要贡献是在保证精度和稳定性的基础上提升了ECV土壤水分产品的空间的覆盖完整性,得到一种高精度、高地表覆盖率的土壤水分融合产品。
本研究以地面站点实测数据作为“理想真值”验证卫星土壤水分产品精度。然而,站点尺度数据与0.25°×0.25°分辨率栅格像元的空间代表性差异巨大,且在下垫面性质复杂的区域尤为显著。因此,基于地面站点实测数据的精度评价能够在一定程度上检验土壤水分融合产品的质量,但不能完全准确地等价于数据集精度。
空值图斑是卫星土壤水分产品普遍存在的现象,本研究基于高精度SMAP对ECV土壤水分产品进行插补,是提升卫星土壤水分产品空间覆盖率的有效方式之一。基于数学模型构建土壤水分与多源地表参数(如降水、气温、植被指数)的映射关系反演土壤水分也是重要的卫星土壤水分插补方法,有待进一步深入研究。
作者分工:王卓颖采集和处理了数据、验证了数据精度、撰写了数据论文;刘杨晓月对数据集的研制做了总体设计、构建了模型和算法等。






