生态系统服务功能模块之防风固沙（一）

前言：（碎碎念：周末偷了个小懒，今天原本打算不写的，晚上打开csdn发现增粉了！！！不行，为了粉丝，今天必须再写一篇）今天分享关于防风固沙的计算，这部分内容需要收集的数据巨多，我分成几个文章进行这部分的说明，下面让我们开始吧！

1、防风固沙的定义

生态系统通过增加土壤抗风能力，降低风力侵蚀和风沙危害的功能。

2、防风固沙的计算公式（以下公式我用的是中华人民共和国国家生态环境标准HJ 1173—2021）

3、防风固沙参数计算

我们今天计算S潜中的气象因子WF（碎碎念：兄弟们顶住，孩子搞这部分崩溃了无数次，哭）

，根据公式，我们分别计算WF、EF、SCF、K‘。

WF：气象因子

（1）风力因子wf

①U1为2米处的起沙风速，根据文献，一般取5m/s。

②u2为气象站2m处月均风速值（m/s）。

③Nd为各月2m处风速大于5m/s的天数。

②与③这两个参数可使用GEE平台ERA5数据集，空间分辨率为0.1°（即10km）。（需要代码的私信我哦）

（2）

:空气密度，g:重力加速度。对于这两个参数，我看好多文献中采用的定值，即，。我在实验时没有采用这种方法，而是通过纬度、海拔等计算的。感兴趣的话关注我，下期主要分享这两个参数的计算。

（3）土壤湿度因子SW

①：潜在蒸散发量

我采用的是中国1km逐月潜在蒸散发数据集，空间分辨率1km，时间范围为1901-2022，可从青藏高原数据中心免费获取，获取链接：国家青藏高原科学数据中心。格式为nc格式，需要转化为Tif。采用arcgis或者编程实现。

②R：月降雨量

中国逐月降水数据集，空间分辨率为1km.时间范围1901—2022, 可从青藏高原数据中心免费获取。获取链接为：https://data.tpdc.ac.cn/zh-hans/data/faae7605-a0f2-4d18-b28f-5cee413766a2

格式为nc格式，需要转化为Tif。采用arcgis或者编程实现。

或者通过GEE的ERA5数据集下载数据。

③I灌溉量

全球灌溉农田灌溉用水量遥感估算数据集，空间分辨率为1km。时间范围2011—2018, 可从青藏高原数据中心免费获取，获取链接国家青藏高原科学数据中心。有文献给出I=0，看大家的研究时间范围吧，如果研究时间不在2011-2018之间，可以另I=0。我由于研究2010、2015、2020年的，所以我没考虑灌溉量了。 

④：每月降雨或灌溉次数或天数

我是通过GEE的ERA5数据集实现的。（需要代码的点个关注私信我哦）

⑤N：实验天数，即30

求出每个月的SW，然后12月求和得最终需要得SW结果。

（4）雪盖因子SD

SD=积雪厚度在25.4mm以下的天数与一年中总天数的比值。

我下载的是中国雪深长时间序列数据集，空间分辨率为2.5km，时间尺度：日，时间范围：1979-2023，数据链接：国家青藏高原科学数据中心。

OK，气象因子WF计算到这里就结束了，我们今天就分享到这里，我们下期再见！同样的以一句话结束今天的分享：幸福是我们奋斗出来的，加油干！！！