一种基于标准化陆地水储量指数的干旱遥感监测方法技术

本发明专利技术公开了一种基于标准化陆地水储量指数的干旱遥感监测方法，包括以下流程：S1、获取基于重力卫星遥感反演的陆地水储量的原始时间序列数据集；S2、依据某一时间尺度的干旱监测应用需求，计算该时间尺度上陆地水储量变化的滑动平均时间序列；S3、由不同时间尺度上陆地水储量变化的滑动平均时间序列计算出对应时间尺度的标准化陆地水储量指数STWSI；S4、依据计算所得的标准化陆地水储量指数，进行干旱监测情况判断。本发明专利技术利用全新定义的标准化陆地水储量指数实现水文、地下水、农业或生态系统的干旱监测，基于卫星遥感陆地水储量变化表征干旱状况，算法简单，参与计算的观测数据易获取，能更准确地描述干旱严重程度。能更准确地描述干旱严重程度。能更准确地描述干旱严重程度。

【技术实现步骤摘要】
一种基于标准化陆地水储量指数的干旱遥感监测方法


[0001]本专利技术涉及一种干旱监测方法，尤其涉及一种基于标准化陆地水储量指数的干旱遥感监测方法。

技术介绍

[0002]干旱是一种多发的极端气候事件，同时也是最具破坏性的自然灾害之一。频繁发生的干旱灾害会给生态环境、农业生产、经济活动、社会民生等带来严重影响。现有用于监测干旱强度的干旱指数包括多种，如（1）标准化降水指数（SPI），仅能反映大气降水供应的短缺引起的干旱变化；（2）标准化降水蒸散发指数（SPEI），需要估算潜在蒸散发（PET），但是潜在蒸散发存在多种算法，不同算法差异很大，导致潜在蒸散发具有非常大的不确定性；（3）标准化土壤湿度指数（SSMI），仅需要土壤湿度观测数据，但是目前观测稀少、时间和空间连续性差，此外，观测仪器安装和维护成本高。
[0003]上述干旱指数的具体算法方案如下：（1）标准化降水指数（SPI）：，其中，P为月累积降水时间序列，为各个月份P的多年平均值，为各个月份P的多年标准差。SPI指数仅能反映大气降水供应的短缺引起的干旱变化。需要指出的是，实际应用中由于降水数据只要服从伽玛（Gamma）分布，SPI是由服从伽玛分布的函数拟合降水的时间序列并进行标准化后得到。
[0004]（2）标准化降水蒸散发指数（SPEI）：或其中，，为降水P和潜在蒸散发PET的差值。而F为D时间序列服从log
‑
logistic分布的拟合函数：，其中，为参数；该指数算法需要估算潜在蒸散发（PET），但是潜在蒸散发存在多种算法，不同算法差异很大，导致潜在蒸散发具有非常大的不确定性。
[0005]（3）标准化土壤湿度指数（SSMI）：，其中，SM为土壤湿度时间序列；和分别为SM的均值和标准差；该指数算法仅用土壤湿度一个变量作为输入，计算简单，但土壤湿度的观测稀少，因此，该指数在时间和空间的适用性都受到限制。

技术实现思路

[0006]为了解决上述技术所存在的不足之处，本专利技术提供了一种基于标准化陆地水储量指数的干旱遥感监测方法。
[0007]为了解决以上技术问题，本专利技术采用的技术方案是：一种基于标准化陆地水储量指数的干旱遥感监测方法，包括以下流程：S1、获取基于重力卫星遥感反演的陆地水储量的原始时间序列数据集；S2、依据某一时间尺度的干旱监测应用需求，计算该时间尺度上陆地水储量变化的滑动平均时间序列；S3、由不同时间尺度上陆地水储量变化的滑动平均时间序列计算出对应时间尺度的标准化陆地水储量指数STWSI；S4、依据计算所得的标准化陆地水储量指数，进行干旱监测情况判断。
[0008]进一步地，步骤S1中，所获取的重力卫星遥感反演数据包括但不限于GRACE或GRACE
‑
FO TWS时间序列数据。
[0009]进一步地，步骤S2中，某一时间尺度n的陆地水储量变化的滑动平均时间序列TWSA由原始陆地水储量的时间序列TWS在该时间尺度上的滑动平均计算得到，如公式所示：,其中，TWSA为某一陆地区域上某一时间尺度的陆地水储量变化的滑动平均时间序列；TWS为原始陆地水储量的时间序列。
[0010]进一步地，步骤S3中，标准化陆地水储量指数STWSI的指数表征公式为：,其中，或为在某一时间尺度的时间序列服从广义极值分布的拟合函数；为该时间尺度上多年的G时间序列的平均值；是该时间尺度的G时间序列的标准差。
[0011]进一步地，步骤S4中，STWSI小于0，表示陆地地表及地下水含量低于正常值，呈现水分亏缺状态；反之，STWSI大于0，表示陆地地表及地下水含量大于正常值，呈现水分盈余的状态。
[0012]进一步地，干旱等级分为轻度干旱、中度干旱、重度干旱和极端干旱，不同干旱等级间的分隔点为数据的1σ、1.5σ和2σ，σ为标准差；
‑
1<STWSI≤0，表示轻度干旱；
‑
1.5<STWSI≤
‑
1，表示中度干旱；
‑
2<STWSI≤
‑
1.5，表示重度干旱；STWSI≤
‑
2，表示极端干旱。
[0013]本专利技术公开了一种基于标准化陆地水储量指数的干旱遥感监测方法，利用全新定义的标准化陆地水储量指数实现水文、地下水、农业或生态系统的干旱监测；本专利技术基于卫星遥感陆地水储量变化表征干旱状况，算法简单，参与计算的观测数据易获取，能更准确地描述湖泊、农田、森林、草地等陆地生态系统中发生的干旱严重程度。
附图说明
[0014]图1为本专利技术标准化陆地水储量指数（STWSI）的计算流程图。
[0015]图2为本专利技术实施例中鄱阳湖地区2002
‑
2023年在不同时间尺度的标准化陆地水储量指数的时间序列图。
[0016]图3为本专利技术实施例中长江流域在2022年7
‑
9月的标准化陆地水储量指数空间分布图。
具体实施方式
[0017]下面结合附图和具体实施方式对本专利技术作进一步详细的说明。
[0018]本专利技术提出一种基于标准化陆地水储量指数的干旱遥感监测方法，是一种新的基于卫星遥感监测的干旱指数算法，基于卫星遥感陆地水储量负异常变化表征干旱状况，能够准确地描述水文、地下水、农田、森林、草地等陆地生态系统中发生的干旱程度。
[0019]首先，本专利技术定义了一个新的干旱指数，为标准化陆地水储量指数（英文名：Standardized Terrestrial Water Storage Index （STWSI）），该标准化陆地水储量指数在不同时间尺度的计算流程如图1所示，计算步骤如下：S1、获取基于重力卫星遥感反演的陆地水储量的原始时间序列数据集（GRACE TWS）；S2、依据某一时间尺度（如1、3、6、12等月尺度）的干旱监测应用需求，计算该时间尺度上陆地水储量变化的滑动平均时间序列（TWSA）；S3、由不同时间尺度上陆地水储量变化的滑动平均时间序列（TWSA）计算出对应时间尺度的标准化陆地水储量指数；S4、依据计算所得的标准化陆地水储量指数，进行干旱监测情况判断。
[0020]其中，所获取的重力卫星遥感反演数据包括但不限于GRACE或GRACE
‑
FO TWS时间序列数据；某一时间尺度（n）的陆地水储量变化的滑动平均时间序列（TWSA）由原始陆地水储量的时间序列（TWS）在该时间尺度上的滑动平均计算得到，如公式所示：，标准化陆地水储量指数的指数表征公式为：，其中，TWSA为某一陆地区域上某一日或月等时间尺度的陆地水储量变化的时间序列，单位是cm；为该时间尺度上多年的陆地水储量变化的平均值，是该时间尺度上多年的陆地水储量变化的标准差；或,其中，或为TWSA在某一日或月等时间尺度的时间序列服从广义极值分布（GEV）分布的拟合函数；为该时间尺度上多年的G时间序列的平均值；是该时间尺度的G时间序列的标准差。
[0021]计算得到标准化陆地水储量指数本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于标准化陆地水储量指数的干旱遥感监测方法，其特征在于：包括以下流程：S1、获取基于重力卫星遥感反演的陆地水储量的原始时间序列数据集；S2、依据某一时间尺度的干旱监测应用需求，计算该时间尺度上陆地水储量变化的滑动平均时间序列；S3、由不同时间尺度上陆地水储量变化的滑动平均时间序列计算出对应时间尺度的标准化陆地水储量指数STWSI；S4、依据计算所得的标准化陆地水储量指数，进行干旱监测情况判断。2.根据权利要求1所述的基于标准化陆地水储量指数的干旱遥感监测方法，其特征在于：步骤S1中，所获取的重力卫星遥感反演数据包括但不限于GRACE或GRACE
‑
FO TWS时间序列数据。3.根据权利要求1所述的基于标准化陆地水储量指数的干旱遥感监测方法，其特征在于：步骤S2中，某一时间尺度n的陆地水储量变化的滑动平均时间序列TWSA由原始陆地水储量的时间序列TWS在该时间尺度上的滑动平均计算得到，如公式所示：,其中，TWSA为某一陆地区域上某一时间尺度的陆地水储量变化的滑动平均时间序列；TWS为原始陆地水储量的时间序列。4.根据权利要求3所述的基于标准化陆地水储量指数的干旱遥感监测方法，其特征在于：步骤S3中，标准化陆地水储量指数STWSI的指数表征公式为：,其中，TWSA为某一陆地区域上某一时间尺...

【专利技术属性】
技术研发人员：张选泽，任婵月，田静，郑玉峰，苗平，田小强，张永强，
申请(专利权)人：鄂尔多斯市气象局鄂尔多斯市水旱灾害防御技术中心，
类型：发明
国别省市：