民勤绿洲外围新月形沙丘宽高比与移动速率变化特证

2023-03-16何晨晨吴盈盈田永胜王振亭

干旱区研究 2023年2期

何晨晨，吴盈盈，田永胜，马瑞，王振亭

（1.中国科学院西北生态环境资源研究院，沙漠与沙漠化重点实验室，甘肃兰州 730000；2.中国科学院大学，北京 100049；3.甘肃农业大学林学院，甘肃兰州 730070）

河西走廊毗邻腾格里、巴丹吉林和库姆塔格三大沙漠，饱受风沙危害，对中国北方生态安全屏障建设至关重要。走廊东端的民勤地区已有六十余年治沙历史，虽然取得了举世瞩目的成绩［1］，但是土地荒漠化形势依然严峻。绿洲外围水资源匮乏，风力强劲，极易发生风蚀荒漠化。以呈半环状包围着民勤绿洲的连古城国家级自然保护区为例，其总面积的98.41%属于风蚀荒漠化土地［2］。了解风沙地貌形成与运动规律是科学开展风沙工程设计与施工的重要前提。民勤绿洲外围上风向区域内的风沙地貌以固定和半固定的灌丛沙丘为主，新月形沙丘与沙丘链次之，小型风蚀坑偶有分布。流动沙丘在自身不断移动的同时，亦为其附近的风沙活动提供物质来源，对绿洲存在与发展的现实威胁最大。

在众多沙丘类型中，新月形沙丘的研究程度较高［3］。早期野外观测主要以文字描述为主［4］，其工作量大、耗费时间长、形态特征描述不够准确。后来，插钎法、重复测量和航片分析等传统监测方法得到广泛应用［5-6］。目前，遥感影像、GIS、全站仪、实时动态差分GPS、三维激光扫描仪、无人机倾斜摄影测量等各种先进技术被应用于沙丘形态特征测量［7-9］，为深入探索沙丘形态动力学过程提供了宝贵资料。新月形沙丘的形态参数包括沙丘高度、宽度、迎风坡长度、落沙坡长度、迎风坡坡度、两翼长度、开展度和对称度等［10］。各参数之间的关系包括沙丘宽度与高度、迎风坡坡长与宽度、两翼长度与宽度、落沙坡倾角与高度、底面积与周长或体积等。其中，沙丘宽度和高度之间的关系研究最为广泛。由大量实测数据得出的一个经验关系是沙丘宽度与高度正相关［11］，可用线性关系描述［10-12］。按照新月形沙丘的迎风坡长与宽度的比值，将其形态分为肥的、矮胖的、标准的、苗条的［13-14］。按兽角形态将其分为线性的、串珠状的、打结的［15］。除了形态特征，沙丘的演变过程和移动速率是另一个研究重点。沙丘是风力、沙粒和下垫面三者相互作用的产物［16］，其形态演变和移动速率受沙源丰富程度、风速、风向、地形、植被等多种因素影响，不同地区的风沙运移特征有着不同的表现形式，看似普遍的经验和理论关系往往存在很大的区域差异。

河西走廊绿洲边围2006 年与2014—2015 年两个时段的遥感影像显示［2］，民勤地区的新月形沙丘移动最快。本研究拟在前人工作的基础上，利用野外实地测量和高分辨率遥感影像分析沙丘形态特征和移动速率，给出民勤绿洲上风向边缘区域内新月形沙丘形态参数之间及其与移动速率的定量关系，探讨该地区多年防风固沙工程的实际效果。本研究不仅直接服务于当地的防风固沙与生态建设，而且对河西走廊其他地区相关工作具有重要参考价值与借鉴意义。

1 数据与方法

1.1 研究区概况

研究区位于甘肃省民勤县城西北约30 km范围内，主体南起红崖山水库北缘，北至莱菔山脚下（102°45′～103°50′E，38°25′～38°50′N），如图1 所示。该地区属温带干旱荒漠性气候，常年干燥、降水稀少、蒸发强烈、风大沙多。年均降水量116.2 mm［17］，气温7.8 ℃，年内风沙日数139 d，8 级以上大风29 d。区域地形平坦，无地表径流，风向单一，植被稀疏，零星分布有梭梭、沙拐枣、灌丛白刺等植物。土地利用类型以沙地为主，林地占一定比例，耕地和草地零星分布。土壤类型以荒漠土为主。风蚀荒漠化土地既有历史时期自然形成［18］，也有最近几十年开垦后又弃耕退化。较为开阔的前后滩地是新月形沙丘及沙丘链存在所必须的环境条件，同时，该区位于巴丹吉林沙漠盛行风的下风向，沙源供应较为充足。因此，新月形沙丘及沙丘链是本区的主要沙丘类型［19］。该区域历来是当地防风固沙的重中之重。目前，治沙措施以大面积植树造林和围栏封育为主，已在乡道两侧3～5 km内及绿洲边缘部分地区形成了面积不等的防风固沙人工梭梭林。

图1 研究区示意图Fig.1 Map of the research area

1.2 研究方法

1.2.1 沙丘形态参数野外实地测量和高分遥感影像分析是获取沙丘形态参数的主要方法。新月形沙丘高度和宽度等形态参数利用NTS-352 型全站仪直接测得。测量对象主要为形态完整、轮廓显著、相对独立的新月形沙丘。野外测量采用局部坐标系，落沙坡最高点与地面的垂直距离作为沙丘高度；在两翼之间任选一条平行于主风向的直线，两个翼角与其水平距离之和即为沙丘宽度。先后两次测量沙丘形态的时间为2006 年3 月与2022 年7月，沙丘个数分别为27 和25 个。时隔16 a，第一次测量的沙丘［20］大多难以确定其具体位置。相较于上次测量，第二次的空间范围更大，基本覆盖了整个研究区。

1.2.2 遥感影像来源与分析高分遥感影像数据购自高分辨率对地观测系统甘肃数据与应用中心。其中，2013年11月9日和2014年4月21日影像为高分一号卫星数据，其空间分辨率为2 m；2015年8月21日、2016年5月23日、2017年11月21日、2018年7 月21 日和2019 年10 月7 日的空间分辨率为0.81 m 的高分二号影像；2020 年8 月5 日、2021 年8 月9日和2022年5月23日为高分六号影像，其空间分辨率2 m。各时段影像预处理包括：遥感影像经正射校正、融合、镶嵌、几何校正。镶嵌处理过程需要注意选择合适的重采样方法，几何校正利用了地面控制点和几何校正数学模型来矫正非系统因素产生的误差。完成以上步骤后，利用ArcGIS软件从中提取流动新月形沙丘底面的轮廓线，最后计算沙丘宽度和底面质心位置。该地区新月形沙丘形态处于不断变化之中，难以选定轮廓线上的某个或某几个特征点代表沙丘位置。因此，以底面质心位置计算其年际移动距离。以遥感影像测量的沙丘主要位于薛百镇和大坝镇以西，仅有3 个沙丘位于绿洲北缘的莱菔山下，具体位置见图1。

影像分析采用整体坐标系，具体做法如下：（1）按以直代曲方法量测出每个沙丘底面轮廓处若干点在WGS-84 坐标系下的坐标值，计算其在三维直角坐标系下的坐标；（2）利用最小二乘法拟合出包含所有沙丘底面轮廓点的最优平面；（3）在最优平面上任选一点为坐标原点，以主风向为x方向，建立平面直角坐标系，并将所有测量点向该坐标系投影。以投影后的轮廓线计算沙丘宽度与质心位置。其中，确定质心坐标时的底面面积和一阶矩由梯形法计算。因为研究区范围不大且地形非常平坦，这种影像分析方法能够比较准确地提取到沙丘底面轮廓。对于沙丘迎风坡长度和两翼宽度等形态参数，野外实测与遥感影像量测结果之间的相对误差小于5%。

类似于前期探索［20］，仍然利用全站仪数据拟合新月形沙丘的形态参数关系。沙丘移动的具体规律则在理论推导基础上，由遥感影像数据最终确定。野外测量的新月形沙丘既有已被人工林固定的，也有完全流动的，遥感影像上量测的则均为流动沙丘。

1.2.3 风况特征民勤气象站2013—2021 年逐小时风速数据来源于中国气象数据网（http://data.cma.cn/），用于风况特征分析。输沙势常用于评估风况特征与潜在输沙强度，其计算公式为［21］：

DP=U2(U-Ut)t

式中：DP为输沙势，矢量单位为VU；U和Ut分别为10 m 高度的风速和起动风速，Ut取值为6 m·s-1；t为起沙风作用的时间，一般用频率表示。风能环境评价包括总输沙势（DP）、合成输沙势（RDP）、合成输沙方向（RDD）和方向变率（RDP/DP）。根据矢量合成法则，对16个方位的输沙势合成得到合成输沙势和合成输沙方向，而合成输沙势与输沙势之比为方向变率。根据输沙势的大小，将风能环境［21］划分为低能（DP＜200）、中能（200≤DP≤400）和高能环境（DP＞400）。根据方向变率，风能环境划分为高变率（RDP/DP＜0.3）、中等变率（0.3≤RDP/DP≤0.8）和低变率（RDP/DP＞0.8）。

2 结果与分析

2.1 风况特征

输沙势是反映一个区域风沙活动强度的重要指标之一［22］。民勤地区2013—2021 年输沙玫瑰图如图2 所示。年起沙风主要来自于西西北和西北风，历年合成输沙方向介于120°～135°，总体搬运方向为ESE～SE，风况分布为窄单峰型。2013—2021年输沙势依次为40 VU、34 VU、34 VU、26 VU、25 VU、20 VU、20 VU、15 VU和35 VU，年际间输沙势虽有变化，但差异较小，风能环境基本稳定。多年平均输沙势为27.7 VU。根据风能环境划分标准［21］，该地区属于低风能环境（DP＜200）。另外，该区域合成输沙势与输沙势呈现出同样的年际变化，变化趋势基本一致。在2013—2020年间呈不断下降趋势，2021 年出现上升（图3a）。合成输沙方向和方向变率处于不断波动状态（图3b），多年平均方向变率为0.84，属于低变率，表明该地区风向较为单一。但是，2013年、2014年、2019年方向变率小于0.8，属于中等变率。因此，研究区总体输沙势较小，部分年份风向变化较大，新月形沙丘的移动主要受西西北和西北风的影响。

图2 2013—2021年民勤站输沙势玫瑰图Fig.2 Rose of drift potential at Minqin station from 2013 to 2021

图3 2013—2021年民勤站风能环境变化Fig.3 The variations in the wind energy at Minqin station from 2013 to 2021

2.2 高度-宽度关系

新月形沙丘高度与宽度数据可用正比例函数［20］或幂函数［23］拟合。为了方便，进行正比例拟合，结果见图4。总体上看，虽然两次测量时隔16 a，但是形态参数关系基本没有变化，沙丘宽高比约为16。前人大量工作表明其值一般在8～18之间［24］，本研究符合这一规律。多年保持常数值，说明宽高比是新月形沙丘的特征参数。沙丘形态特征受风况、沙源和地形等因素影响。对于同一区域，上述风况特征表明多年输沙势变化不大，属于低风能环境。多年平均方向变率属于低变率，风向较为单一。在无人为干扰的情况下，新月形沙丘的参数关系也应无明显变化。

图4 新月形沙丘高度与宽度呈正比例关系Fig.4 Proportional relationship between height and width of barchans

2022 年实测数据略显分散，主要原因有二：第一，部分沙丘迎风坡上的人工梭梭林影响了其形态-动力学过程。一般而言，气流在沙丘迎风坡加速，绕过沙丘顶部在背风坡形成回流，迎风坡气流加速和背风坡回流会影响沙丘的侵蚀速率，从而影响沙丘的宽度、高度等形态参数。对于不同宽高比的沙丘，这一过程也会有所变化。沙丘宽高比越小，气流在迎风坡加速越快，背风坡回流区长度和强度均会增大［25］。地表植被覆盖会引起下垫面粗糙度改变，地表抗风蚀能力增强，从而影响沙丘近地面流场分布。在2022年测量的部分沙丘中，迎风坡存在的人工植被会减弱气流在迎风坡的加速作用，同时沙丘背风坡回流区范围更大。流场的这种变化会使得沙丘高度快速增长，以及背风坡沙尘堆积［26］。植被冠层拦截悬浮在空中的沙尘，使得沙丘上风向沙源供应不足，导致其形态发生变化。同时，迎风坡植被阻挡会使气流向沙丘的两翼绕流，风沙流会将部分沙粒带向沙丘两翼，导致沙丘宽度发生变化。若沙丘两翼附近存在其他沙丘，气流流过沙丘会产生加速效应，这种加速效应会加快两翼的侵蚀速率。

第二，不同区域的风力状况和地表条件略有差别。民勤绿洲外围的新月形沙丘主要分布在前后平坦的砂砾质滩地上。其位于巴丹吉林沙漠的下风向，沙源供应充足。但由于不同区域地表起伏程度有所差异，导致沙丘丘间地不同位置风沙流存在明显差异。不同风力作用下，距离沙源越近，输沙率越大，随着远离沙源，输沙率减小［27］。新月形沙丘近地表风速和输沙率研究显示，由于气流在迎风坡的加速作用，迎风坡风速比来流风速增加25%，沙丘翼角到顶部输沙率增加了1 倍，不断增加的输沙率反映风力侵蚀的增加［28］。2022 年测量沙丘宽度和高度总体上相较于2006 年更大，即2022 年多为体积更大的新月形沙丘。沙丘迎风坡气流的加速作用更强，丘顶的风力侵蚀作用更大，沙丘的沙源丰富度和地表条件也更加复杂多样。

2.3 移动速率

2.3.1 单个沙丘移动轨迹与形态变化沙丘的移动速率可以估算沙物质的输移量，为防风固沙实践提供指导。在风力作用下，沙丘往往通过迎风坡侵蚀与背风坡堆积而发生移动。虽然宽高比保持不变，但单个沙丘的形态在移动过程中却是不断变化的。在传统的风沙地貌研究中，往往以轮廓线上某个点或某几个特征点的平均移动距离表示整个沙丘的移动。当形态变化较大时，误差很大。因此，以底面质心标识沙丘位置是更合理的做法。

图5给出了位于102.99°E，38.79°N和102.83°E，38.68°N两个沙丘的底面轮廓在2013—2022年间的变化情况。图5a中沙丘高度为9.3 m，两翼宽度分别为89 m和76 m。以2013年为参考，底面轮廓变化能够反映沙丘详细的年际移动过程。总体上，沙丘形态维持较为完整，侧向位移微小，新月形形状明显。底面质心变化反映出沙丘整体移动缓慢。沙丘年际间的形态变化主要在迎风坡基线、脊线和翼角处比较明显。沙丘迎风坡基线和脊线所代表的沙丘整体上沿主风向移动，但2020 年和2021 年出现反向移动。沙丘移动方式包含前进式、摇摆式和摆动前进式［29］。当沙源供应不足时，风沙流变为不饱和风沙流，以及常年多个风向的出现会导致新月形沙丘的反向发展。图5b 中沙丘的两翼宽度约为60 m 和40 m。沙丘的底面轮廓反映出其形态变化很剧烈，各形态参数变化大，沙丘整体沿主风向快速前移。在2013—2022年间，沙丘沿主风向主要做前进式移动，少数年份存在反向移动。在反向风风蚀的作用下，沙丘落沙坡顶部会出现反向堆积和脊线反向移动的现象，从而促进了沙丘的反向发育。根据移动速率，流动沙丘可分为四类［30］：慢速（年平均移动距离＜1 m）、中等速度（年平均移动距离＜1～5 m）、快速（年平均移动距离＜5～10 m）、极快速（年平均移动距离＞10 m 以上）。上述高大沙丘年均移动距离小于3 m，属于中等速度；低矮沙丘，沙丘年均移动距离大于10 m 以上，属于极快速移动。沙丘移动速率与沙丘的体型有关，体型越大，移动速率越小。

图5 单个沙丘的移动轨迹与形态变化Fig.5 Trajectory and morphological variation of individual barchans

图5中的质心轨迹显示两个沙丘均存在显著的反向移动。而当主风向和次风向大致相反时，才会出现这种运动。常兆丰等［31］对甘肃河西绿洲边缘新月形沙丘分析发现，单一主风向是维持新月形沙丘形态的关键因素，促使沙丘顶点和沙脊线重合；而反向风促使沙丘顶点与沙脊线分离。对民勤沙区风速风向数据的统计分析显示［32］，起沙风主要为春季的西北风，其次为夏季的东南风。反向风除了导致新月形沙丘的这种整体反向移动外，还改变了脊线处的风沙流结构并造成丘顶与脊线相对位置的变化［33］。沙丘移动主要受风速和风向等因素的影响。结合风况分析，研究时段输沙势总体较小，属于低风能环境。单个沙丘移动轨迹反映出沙丘总体向前移动，但2014年沙丘出现明显的反向移动（图5b），2014年方向变率为0.75，属于中等变率，说明风向变化较大，风向多变是沙丘反向移动的重要因素。由于民勤地区常年受西北风控制，东风、东南风较弱，难以完全抵偿西北风产生的作用。即使年内风向变率较大，但年输沙势较小不足以使沙丘发生反向移动。因此，沙丘年际间的反向移动可能是强沙尘暴事件所致。以2007 年3 月27 日一场强沙尘暴为例［34］，在强沙尘暴事件前后，对图5a 沙丘监测发现：脊线高度减小（15±5）cm，长度增加（7.00±2.00）m；移动距离（21±5）cm，沙丘体积减小（5.1±2.0）×103m3。这说明强沙尘暴事件对沙丘的移动速率和形态特征有重要的影响，也是导致沙丘反向移动的重要原因。此外，低矮沙丘底面轮廓的面积随着时间推移不断缩小，这表明在移动过程中体积在减小，沙丘向着缩小的趋势发展。而高大沙丘保持相对稳定，造成这一结果的原因与沙丘自身的体积大小有关。

2.3.2 移动速率与形态参数关系当沙丘形态和来流情况保持不变时，可由质量守恒定律导出移动速度与沙丘高度成反比的理论关系［35-36］，有时，新月形沙丘在几个月内的总体移动几乎完全由一两场强沙尘暴造成。同时，反向风的存在及遥感影像成像时间的年度差异也为确定沙丘沿主风向的移动速率带来不小困难。为了减少逐年计算带来的不确定性，本文以2013 年质心位置为参考点，由2014—2022年间质心移动的最大距离计算9 a间的平均移动速率，沙丘宽度取为10次影像分析的平均值。图6 给出了年移动速率随沙丘宽度的变化情况，沙丘宽度范围在45～300 m，移动速率在1～16 m·a-1。其中，快速移动沙丘占52.63%，极快速移动沙丘占15.79%，中等移动速度沙丘占31.58%。因此，该区域沙丘主要以快速移动沙丘为主，移动速率大于5 m·a-1。

沙丘移动速率与宽度之间建立关系之后，对研究区风况和沙丘宽度的测量就可以快速估计沙丘的移动速率。通过宽度反映沙丘的移动速率大大降低了研究的难度。同一风力条件下，沙丘的移动速率不同。这表明风况对沙丘移动速率的影响被其他因素削弱，导致移动速率差异［37］。前文已证实沙丘宽高比保持常数的经验关系。因此，移动速率应该反比于沙丘宽度。由最小二乘法得到的拟合优度为0.83，比例常数为718.52 m2·a-1。对沙丘的移动方向［38］、移动距离［39］、移动速率［40］研究发现，新月形沙丘的移动速率与高度、宽度均成反比，与本文结果一致。同时，沙丘移动速率和宽度的比值变化还表明，随着宽度增大，沙丘移动速率的变化趋于平缓。沙丘宽度由50 m 增加到150 m，移动速率的变化范围在3～16 m·a-1；而沙丘宽度由200 m增加到300 m，沙丘移动速率的变化范围在2～7 m·a-1。这说明沙丘达到一定的宽度后，移动速率变化的波动范围更小，同一风况下，沙丘更趋于稳定。对独立的新月形沙丘，沙丘宽度与移动速率之间存在的反比关系，可以用不同宽度沙丘之间的沙粒质量差异来解释。沙丘的移动是沙粒从迎风坡侵蚀重新分布到背风坡的过程。宽度更大的沙丘需要更长的时间来重建和平衡这一过程，因此移动速率更慢［41］。这一过程可能导致沙丘在移动过程中形态特征发生变化。沙丘动力学模型研究表明，恒定沙源条件下，在风沙流中不同大小的新月形沙丘可能处于不同的不平衡状态［42］。结合单个沙丘多年移动轨迹的变化以及沙丘移动速率与宽度的关系可知，沙丘宽度不但反映移动速率的快慢，而且由于移动速率之间的差异，不同宽度的沙丘形态向不同方向发育。随着时间推移，宽大沙丘形态基本保持不变或者变大，而窄小沙丘倾向于缩小或者消失。

3 讨论

3.1 风况

风况是影响沙丘形态发育的重要动力因素，风速和风向决定沙丘的高度、宽度、体积和移动速率等。沙丘形态演化和移动不仅受大气环流的控制，同时受局地气流的影响［43］。宏观上，区域风况受大气环流的控制，已有研究对风况和不同类型的沙丘建立联系。窄单峰型与新月形沙丘相联系，双峰型通常与线形沙丘相联系，而复合型则与星状沙丘相联系［44］。在沙源和风向变率均低的地区会形成新月形沙丘，沙源的增加导致新月形沙丘合并成横向沙丘，而风向变率的增加导致双向风环境下线状沙丘的形成。研究区风况特征显示年起沙风主要来自于西西北和西北风，风况分布为窄单峰型，多年风向变率为低变率，风向较为单一，有利于新月形沙丘的形成和演化。因此，沙丘宽高比常年保持常数值（图4）。沙丘不仅对区域风况做出响应，沙丘局部也会产生复杂的二次流动。沙丘背风坡可形成分离流、重附未偏向流、重附偏向流和反向流［45］。沙粒从迎风坡被侵蚀并在背风坡沉积时，沙丘向下风向移动。沙丘移动是风沙量从迎风坡重新分布到背风坡的结果。有研究认为这一过程处于平衡状态，也有认为是非平衡状态。对两个典型沙丘的移动轨迹分析，高大沙丘多年形态维持稳定，移动速率缓慢（图5a）。沙粒被侵蚀后重新分布可以认为是平衡状态。相对低矮沙丘形态变化剧烈，沙丘底面轮廓不断减小反映出沙丘体积在减少（图5b），这表明沙丘在移动过程中侵蚀量大于堆积量，使得沙丘不断缩小，移动过程可以认为是非平衡状态，沙丘最终可能会缩小或者消失。因此，区域大气环流和沙丘局部产生的二次流动共同塑造沙丘的形态，影响沙丘移动。

3.2 植被

随着民勤绿洲防风固沙工程的实施［1］，沙面固定限制了风沙流的沙源供应，改变了地表下垫面的性质。防沙固沙工程主要是通过植物（或生物）措施、机械措施、化学措施等降低沙源供应、改变下垫面性质达到阻沙、固沙和输沙的目的。主要影响因素包括沙粒粒度变化、水分含量、粘合剂、植被种类、植物结构、形态和密度，以及其他表面粗糙度。植被作为沙源供给和风沙运输的主要限制因素，本研究区人工植被是影响沙丘形态和移动的主要因素，是沙丘宽高比分散的原因之一。沙生植物具有发达的根系，能固结周围的沙粒，植被腐烂后有利于有机质的聚集，促进沙的成土作用和固定，同时改变土壤化学性质［46-47］和粒度分布［48］；同时植物冠层能降低近地表风沙运输的有效性。对裸露沙丘和植被覆盖沙丘监测发现［49］，植被覆盖沙丘的移动速率与高度的相关性比裸露沙丘更差。与裸露沙丘相比，固定沙丘的两翼、迎风坡、背风坡、丘顶的移动速度存在差异。植被盖度增加，地表输沙能力和沙源供应受到限制，新月形沙丘会逐渐向抛物线形沙丘转变［50］。地表植被覆盖率仅为15%时，大部分沙粒运输得到抑制［51］。植物冠层对风沙流的影响取决于植物密度、形状或形态、分布和高度［52］，群落高大、枝叶密集的植物种固沙能力更强。对不同类型植被覆盖下沙丘表面沉积物粒度分析［53］，地表沉积物细组分含量与植物冠幅、枝高和分枝数成正相关。在低密度的情况下，单个灌木形成的湍流和涡旋脱落又加强了地表的风沙侵蚀［54］。植物高度增加，导致沙丘高度增加，沙丘长度减少［55］。研究区流动沙丘各参数间表现出良好的相关性，沙丘体积越大，移动速率越慢。同时移动速率反比于沙丘宽度（图6），沙丘总体上向东南方向移动，这些结论与前人研究基本一致［56］。不同之处在于部分新月形沙丘存在人工梭梭林，导致参数间的相关性有所减弱，沙丘宽高比分散。植被对沙丘形态特征和移动速率［57］、近地面流场［58-59］、沙丘表层粒度分布［60］等方面影响深远。植被可以增加地表粗糙度，拦截气流中携带的沙物质，影响沙粒在植物附近沉积［61］。同时，植物削弱了近地面风速，减少被吹扬搬运沙粒的数量，这都会导致沙源供应不足，从而影响沙丘的形态。