一种基于星-地观测的多年冻土活动层厚度估算方法

本发明属于遥感领域，涉及一种活动层厚度计算方法，尤其涉及一种基于星-地观测的多年冻土活动层厚度估算方法，基于星-地观测数据，利用stefan公式对多年冻土区活动层厚度进行估算。

背景技术：

1、中分辨率成像光谱仪（modis）是搭载在美国宇航局发射的 terra 和 aqua 卫星上的一个重要的传感器。其获取的数据在全球免费的，是一种不可多得的、廉价并且实用的数据资源。modis数据产品具有光谱范围广、数据接收简单和更新频率高等特点，已广泛用于陆地、大气和海洋的综合分析以及森林、草原火灾监测和救灾等环境遥感分析。地表温度（land surface temperature, lst）是地球表面与大气之间界面的重要参量之一。它不仅是气候变化的敏感指示因子和掌握气候变化规律的重要前提，还是众多模型的直接输入参数，在气象气候、环境生态、冻土领域有广泛的应用。如今，随着地学及相关领域研究的深入和精细化，学术界对卫星遥感的全天候地表温度具有迫切的需求。

2、stefan模型是1891年奥地利科学家 josef stefan研究北极海冰厚度首次提出，之后广泛用于计算湖冰、海冰等冰体的厚度。该模型假设冰体内部的热传导速度足够快且内部温度呈线性变化。因此，当冰面温度低于相变温度时，冰体下部与湖水接触面的温度就与相变温度相同。这里的相变温度指的是物质在不同相之间发生转变时所需达到的临界温度。1943年，berggren提出冰水相变过程中的潜热释放或潜热吸收远大于干土本身热容变化。这说明，当土层内含水量较大时，将 stefan 方程应用到土壤中与应用到冰体内造成的差别并不大。因此可以用土壤中冰体的热容代替对纯冰体的热容，将其应用到多年冻土活动层厚度估算或季节性冻土的季节冻深。目前，stefan 模型是国内外用于计算多年冻土冻结和融化深度最常用的经验公式。它充分考虑了气候条件、土壤热属性和水分条件，形式简单，驱动参数少，模拟效果较好，既可以用于模拟单点的冻结融化深度，也可以较方便地模拟大尺度的活动层厚度空间分布。

3、这一模型局限于参数选取的精确性与广泛性。modis卫星可以提供全天候地表温度产品，为模拟偏远地区大范围活动层厚度提供了可能。但是，modis地表温度产品易受云或其他大气干扰，造成数据污染，对环境恶劣、难以抵达的野外多年冻土区域提供的地表温度不准确。因此，不能直接使用modis地表温度产品进行活动层厚度估算，需要对其进行校正才能需用于估算活动层厚度。在研究中提供了估算活动层厚度可能的同时，更能在实际上减少人材损失。

4、由于stefan模型参数较少，在多年冻土活动层冻融过程的研究中得到了广泛的应用。但是，单一、传统的stefan模型难以反演大范围区域的活动层厚度。随着遥感卫星技术的发展，将卫星遥感数据同地面观测资料结合在大尺度多年冻土的研究取得了很好的效果。活动层厚度(alt)的变化会严重影响寒区地表的水文循环、生态系统和碳通量的改变。多年冻土区大多地处偏远，缺乏多年冻土的分布和活动层厚度的长期观测数据，使得多年冻土层和活动层的温度和厚度的长期变化大多无法准确评估。基于modis卫星提供的地表温度数据产品，通过stefan公式计算多年冻土活动层厚度，对研究全球气候变化预测、基础设施建设、填补多年冻土区域活动层厚度变化领域研究不足具有重要意义。

技术实现思路

1、本发明为了弥补多年冻土区域活动层厚度变化研究的不足，提供一种基于星-地观测的多年冻土活动层厚度估算方法，其模型能够计算得到多年冻土区活动层厚度。

2、本发明的上述目的通过以下技术手段实现：

3、一种基于星-地观测的多年冻土活动层厚度估算方法，包括如下步骤：

4、步骤1、选定多年冻土区域并选定年份区间，依据选定的多年冻土区域和选定的年份区间分别获得对应的modis日均地表温度数据，以及土壤气候站点日平均地表温度数据、活动层厚度和设定深度处的土壤年均含水/冰量；

5、步骤2、利用土壤气候站点日平均地表温度数据对整个选定的多年冻土区域modis日均地表温度数据进行异常值去除和校正，得到校正后的modis日均地表温度数据；然后对于整个选定的多年冻土区域中各个空间位置，分别累加计算每一年中温度0摄氏度以上的校正后的modis日均地表温度数据得到对应空间位置的累积正积温；

6、步骤3、对于选定的多年冻土区域中各个土壤气候站点对应的空间位置，利用选定的年份区间内各年份中土壤气候站点的活动层厚度、设定深度处的土壤年均含水/冰量和对应年份累积正积温，推算对应年份的土壤导热系数，然后计算选定的年份区间所有土壤气候站点的土壤导热系数均值；

7、步骤4、将上述步骤得到的土壤导热系数均值、选择某一年份的整个选定的多年冻土区域的累积正积温，和所有土壤气候站点设定深度处的土壤年均含水/冰量的均值数据输入stefan模型中，计算得到选定的多年冻土区域中各个空间位置在选择的年份中活动层厚度计算值。

8、如上所述步骤1中依据选定的多年冻土区域和选定的年份区间分别获得对应的modis日均地表温度数据，以及土壤气候站点日平均地表温度数据、活动层厚度和设定深度处的土壤年均含水/冰量的具体方法为：

9、步骤1.1利用gee平台下载选定的多年冻土区域在选定的年份区间内的modis日均地表温度数据，另外在选定的多年冻土区域中土壤气候站点下载相同时间段内的土壤气候站点数据，其中土壤气候站点数据包括：地表温度，不同深度的土壤含水量及监测的活动层厚度；

10、步骤1.2对土壤气候站点数据进行统计得到土壤气候站点日平均地表温度数据、活动层厚度和设定深度处的土壤年均含水/冰量；

11、步骤1.3将步骤1.1下载的modis日均地表温度数据和步骤1.2统计得到的土壤气候站点日平均地表温度数据、活动层厚度和设定深度处的土壤年均含水/冰量按照相同年份归类，建立以年为单位分组的多年冻土活动层厚度数据集。

12、如上所述步骤2中利用土壤气候站点日平均地表温度数据对整个选定的多年冻土区域modis日均地表温度数据进行异常值去除和校正得到校正后的modis日均地表温度数据的具体方法为：

13、步骤2.1在选定的多年冻土区域中选择所有监测地表温度的土壤气候站点，依据土壤气候站点的空间位置坐标提取相同空间位置的modis日均地表温度数据；

14、步骤2.2将统计的土壤气候站点日平均地表温度数据和相同时间以及相同空间位置的modis的日均地表温度数据两倍标准差之外的数值视为异常值去除，然后利用去除异常值之后的相同时间和相同空间位置的土壤气候站点日平均地表温度数据和modis日均地表温度数据进行线性拟合，通过拟合得到拟合参数对整个选定的多年冻土区域的modis日均地表温度数据进行校正，得到的校正后的modis日均地表温度数据。

15、如上所述步骤3中利用选定的年份区间内各年份中土壤气候站点的活动层厚度、设定深度处的土壤年均含水/冰量和对应年份累积正积温，推算对应年份的土壤导热系数的具体方法为：

16、根据土壤气候站点的活动层厚度、设定深度处的土壤年均含水/冰量和对应年份累积正积温，带入如下公式中分别推算对应年份的土壤导热系数，，

17、其中，z表示土壤气候站点的活动层厚度；表示土壤导热系数；表示累积正积温；为温度0摄氏度以上的校正后的modis日均地表温度数据；表示水或冰的密度920；为设定深度处的土壤年均含水/冰量； l表示冰融化为水的潜热334000。

18、如上所述步骤4中通过arcmap的栅格计算器完成计算过程。

19、本发明的优点及有益效果如下：

20、（1）普适性。本发明提出的一种基于星-地观测的多年冻土活动层厚度估算方法，stefan模型能计算得到多年冻土区域活动层厚度，广泛且灵活运用于寒冷区域的活动层厚度的估算工作。

21、（2）高精度性。一种基于星-地观测的多年冻土活动层厚度估算方法，相比传统多年冻土活动层厚度计算方法，具有更大的空间范围和更长的时间序列，能够捕捉到更细节的多年冻土活动层厚度变化趋势和空间分布。

22、（3）实用性。本发明提出的一种基于星-地观测的多年冻土活动层厚度估算方法，能够弥补多年冻土区域活动层厚度变化研究方面的不足。