2. 海原县自然资源局,755200,宁夏海原;
3. 内蒙古自治区水利科学研究院,010050,呼和浩特;
4. 宁夏灵武白芨滩国家级自然保护区管理局,751400,宁夏灵武
中国水土保持科学 ![]() ![]() |
沙丘移动是风沙活动的主要表现形式,也是判断当地生态环境改善与否的重要指标[1]。有关沙丘移动[2-7]、沙丘形态变化[8-11]和沙丘演变的影响因素[1-3, 11-13]研究,均揭示沙丘动态变化是一种多因素主导下的复杂地貌演变过程。影响沙丘移动的因素有沙丘自身的内在因素和外部环境的外在因素。内在因素包括沙丘高度、宽度、底面积、体积、沙粒机械组成、沙粒含水量等,特别是沙丘移动速率和沙丘高度呈明显负相关关系[5, 14]。沙丘外在因素主要指风、气温、降水量、植被等。众多学者认为沙丘移动速率密切依赖于风速大小[5-6, 14-16],通过野外观测、遥感影像监测等手段证实风在沙丘移动和沙丘形态演变过程中均扮演至关重要的角色[6, 15-17];植被覆盖度(fractional vegetation coverage, FVC)是影响沙丘移动的又一关键因素[1, 12-13, 18],相关研究[12]均证明FVC越高,沙丘移动速度越小;移动速率还受气温影响,气温升高会导致沙粒干燥度增加,从而导致沙丘移动速率增大;降水量通过影响沙粒的黏聚性以及团聚能力从而控制沙丘移动[14];此外,地表温度也是一个季节性影响因子,在冬季地温降至0 ℃以下,沙丘上部湿沙层冻结,可影响地表流沙的可蚀性从而降低风蚀[19],由此影响沙丘的移动和形态演变。
综上,沙丘移动是多因素作用结果,某一时段内各因子的同时作用最终决定该时段沙丘的移动速率,但以往欠缺系统分析各影响因子对特定时段内沙丘移动的贡献率,从而难以判断沙丘移动的关键控制因子。鉴于此,笔者选择灵武市白芨滩国家级自然保护区流动和半固定沙丘,应用高分辨率遥感影像计算多年间沙丘的移动速率,利用软件Canoco 5.0定量分析该区域各环境因素对沙丘移动速度的贡献率,以期为今后深入认识外界环境因子对沙丘移动的作用以及荒漠化防治工作提供理论依据和科学指导。
1 研究区概况白芨滩国家级自然保护区地处毛乌素沙地西南部的宁夏灵武市境内(图 1),2000年4月被国务院确定为国家级自然保护区(E 106°20′22″~106°37′19″, N 37°49′05″~38°20′54″),平均海拔1 200 m以上,总面积748.43 km2[20]。
![]() |
图 1 研究区示意图 Fig. 1 Location of the study area |
白芨滩自然保护区属大陆性季风气候,主导风向冬季为西北风,夏季为东南风,气候特点四季分明、春迟秋早、干旱少雨、蒸发量大、温差较大、日照充足、风沙天气较多,多年平均气温10.40 ℃,多年年均降水量200 mm左右,全年日照时间3 011 h[20],多年平均蒸发量近2 000 mm,多年年均风速2.5 m/s。
研究区位于自然保护区中部,面积270 km2。区域内沙丘广泛分布,沙丘类型主要为新月形沙丘链或格状沙丘,地貌景观以自然的流动沙丘和人工种植柠条等固沙灌木后形成的固定半固定沙丘交错构成。生长的植被为多年人工治沙直接或间接成果,常见的有柠条(Caragana korshinskii)、沙芥(Pugionium cornutum)、沙米(Agriophyllum squarrosum)、花棒(Corethrodendron scoparium)、沙拐枣(Calligonum mongolicum)等。
2 数据与方法笔者测量沙丘移动速率的遥感影像来自于软件91卫图助手,8期影像的拍摄时间见表 1,所有影像的空间分辨率均为0.47 m。各时段的植被盖度、平均气温、平均地温、累积降水量数据均来自谷歌地球引擎(Google Earth Engine,GEE),其中FVC利用像元二分模型计算获取[20]。风况数据来源于美国国家气候数据中心,下载灵武气象站的风速和风向数据用于计算输沙势。
![]() |
表 1 沙丘地貌格局参数的变化 Tab. 1 Variation of landscape pattern parameters of dunes |
选取白芨滩自然保护区内3个呈倒三角分布的沙丘地为样地(图 1),每个样地面积300 m×300 m。沙丘地貌格局特征参数有沙丘个数、脊线总长、脊线平均长度、沙丘间距以及缺陷密度[21]。笔者利用ArcGIS对8期影像上的沙脊线等进行矢量化,再统计计算出各时段沙丘地貌格局的特征参数值。
单个沙丘的形态特征参数用长度、宽度和高度表示[17]。同样用ArcGIS提取8期影像上沙丘的长度和宽度。沙丘高度用沙丘背风坡长度与休止角取32°时的正切值乘积来计算[22]。
2.2 合成输沙势及沙丘移动距离(速率)的计算采用下载的风况数据计算合成输沙势(resultant drift potential,RDP)和合成输沙方向(resultant drift direction,RDD),计算输沙势的临界起沙风速设为6 m/s[23]。
每个样地内用五点法取5个沙丘,3个样地共选取15个沙丘。同一沙丘在2期影像的沙脊线发生偏移,偏移距离代表沙丘移动距离。在2期影像各自的沙脊线上分别选取10个测点并连线形成测距线,测距线与影像间隔时段内的RDD平行,每个样地内50条测距线的平均长度代表该样地沙丘的平均移动距离d。沙丘移动速率:
$ v=\frac{\sum\limits_{i=1}^{10} d_i}{50 t} 。$ |
式中:v为沙丘移动速率,m/a;di为测距线长度,m;i=1,2,…,10;t为影像间隔时间,a。
2.3 冗余分析计算贡献率采用软件Canoco 5.0,对选取的自变量RDP、平均气温、平均LST、累积降水量、FVC,和对应的因变量沙丘移动速率,通过RDA获取各环境因素对沙丘移动速率的贡献率。
3 结果与分析 3.1 沙丘地貌格局的变化特征如表 1所示,自2008年1月至2021年1月,研究区内沙丘脊线数量由2008年的192条减少到2021年的101条,减少47.4%;沙脊线平均长度呈递增趋势,2021年是2008年的1.7倍,13 a间增长25.7 m。即研究区2008—2021年的沙丘地貌格局变化表现为沙丘合并,沙丘数量减少但单个沙丘体量增大。
3.2 沙丘移动速率的变化特征沙丘移动速率整体上呈减小趋势,自2008至2021年各时段的移动速率呈极显著的线性递减(P<0.01)(表 2)。依据移动速率大小相近原则,可划分为5个阶段,每个相邻阶段的移动速率相差0.22~0.28 m/a,呈阶梯式稳步下降趋势。第Ⅰ至第Ⅴ阶段的沙丘移动速率依次为:5.62、5.34、5.12、4.86和4.64 m/a。
![]() |
表 2 各时段沙丘移动速率 Tab. 2 Moving rates of dunes in each period |
沙丘自身形态影响沙丘移动,笔者通过测量2020年1月影像上15个沙丘的宽度、背风坡长度、高度及2020-1-11—2021-1-11间沙丘的移动速率并进行拟合分析,发现沙丘移动速率与沙丘宽度、高度均有极显著的线性负相关关系(P<0.01),表明高度和宽度是控制沙丘移动速率的重要内因(图 2)。
![]() |
图 2 沙丘移动速率与宽度、高度的关系 Fig. 2 Correlation between dune moving rate and width, height |
对沙丘缺陷密度和移动速率进行拟合分析(图 3),发现两者间存在极显著的线性正相关关系(P<0.01),即沙丘缺陷密度越大,沙丘移动速率越快;而各时段的RDP与沙丘移动速率也存在极显著的线性正相关关系(P<0.01),即RDP是控制沙丘移动的主要外因(图 4)。
![]() |
图 3 沙丘移动速率与缺陷密度的关系 Fig. 3 Correlation between defect density and moving rate of dunes |
![]() |
RDP: Resultant drift potential. The same below. 图 4 沙丘移动速率与合成输沙势的关系 Fig. 4 Correlation between RDP and moving rate of a dune |
各环境因子对沙丘移动速率的解释度由大到小依次为RDP(78.7%)、FVC、累积降水量、平均地温、平均气温。其中,RDP的P值在0.006水平上,表明其对沙丘移动速率的影响极显著,其余4个因素均未达到显著水平(表 3)。
![]() |
表 3 各环境因子对沙丘移动速率的贡献 Tab. 3 Contributions of various environmental factors to the moving rate of dunes |
本研究印证风是决定沙丘移动速度的主要因子[5-6, 14-16],而其他几个环境因子的贡献率不高。以往研究表明植被通过增加地表粗糙度,削弱近地表层风速,迫使沙丘移动减慢或终止[1],但本研究的植被贡献率仅为6.7%,此低贡献率可能与FVC整体偏低有关。沙丘地可依据FVC大小划分为:流动沙丘地(0<FVC<5%)、半流动沙丘地(6%<FVC<20%)、半固定沙丘地(21%<FVC<50%)、固定沙丘地(51%<FVC<1)[13]。根据以上划分方法,本研究区整体处于半流动沙丘类型(仅有2个时段的FVC略高于20%),且样地内植被分布不均,<5%的流动沙丘面积所占比例>40%,因此FVC的贡献率较低。
沙丘表层土壤水分含量越高,其抗风蚀能力越强[14]。降水后沙丘沙的渗透性极强,且沙丘表面蒸发量大,因此降水对沙丘表面可蚀性的增强作用主要发挥在降水后的几天内;降雨也通过增加土壤含水量促进植物生长间接降低沙丘的移动,但由于本研究区样地内FVC偏低甚至无植被,所以降水通过植被影响沙丘移动的作用可忽略,由此解释2008—2021年间研究区降水量对沙丘移动速率的影响极小。
4.2 沙丘的地温—水分—抗蚀—移动四者间的联动关系本研究中2015-7-31—2019-10-30和2019-10-31—2020-1-11这2个时段,后一时段的RDP远大于前一时间段(21.25 UV>9.29 UV),则沙丘移动速率应大于前一时段。但因白芨滩自然保护区沙丘地的冻结期为11月中旬至翌年3月上旬[19],后一时段的地温大部分时间都在0 ℃以下,致使沙丘地表的土壤水变为固态,且地温可通过影响土壤水分含量和土壤水分状态(液态或固态)控制沙丘表面的可蚀性[24],沙丘表面冻结后其易蚀沙源量大大减少[25],因此即使后一时段RDP较高,其沙丘移动速率仍与前一时段的基本相同(相差5 cm/a)。综上,沙丘的地温—水分—抗蚀—移动四者间为层层递进式的控制,地温长时段在0 ℃以下时,即季节性冻结期内地温可能是沙丘移动速率的另一关键决定因子。
5 结论1) 2008—2021年间,宁夏灵武市白芨滩自然保护区流动沙地沙丘脊线总长度逐步减少,沙丘脊线平均长度随脊线数的减少而增大,表明沙丘在多年移动过程中发生了合并。
2) 沙丘平均移动速率随沙丘的合并逐渐减慢,由2008年的5.62 m/a降至2021年的4.64 m/a。
3) 沙丘高度、宽度和沙丘缺陷密度是影响沙丘移动速率的内因,沙丘高度宽度与移动速率均呈极显著的线性负相关关系,而缺陷密度与之呈极显著的线性正相关关系(P<0.01)。
4) 基于Canoco的分析表明RDP是唯一对白芨滩流动沙丘移动有显著控制作用的外因,贡献率达78.7%。但在季节性冻结期内,地温可能是影响沙丘移动速率的另一关键外因。
[1] |
许明静, 吕萍, 肖南, 等. 毛乌素沙地西北部植被覆盖对沙丘移动的影响[J]. 中国沙漠, 2020, 40(4): 71. XU Mingjing, LU Ping, XIAO Nan, et al. Effect of vegetation cover on dune migration in northwest Mu Us Sandy Land[J]. Journal of Desert Research, 2020, 40(4): 71. |
[2] |
MAHMOUD A, NOVELLINO A, HUSSAIN E, et al. The use of SAR offset tracking for detecting sand dune movement in Sudan[J]. Remote Sensing, 2020, 12(20): 3410. DOI:10.3390/rs12203410 |
[3] |
SHI Xuegang, LIU Shizeng, CHANG Zhaofeng. The barchan dunes movement regularity and its influence factors: A case study in the Hexi Desert area of Gansu[J]. International Journal of Ecology, 2017, 6(2): 93. DOI:10.12677/IJE.2017.62011 |
[4] |
YANG Zhuanling, QIAN Guangqiang, DONG Zhibao, et al. Migration of barchan dunes and factors that influence migrationin the Sanlongsha dune field of the northern Kumtagh Sand Sea, China[J]. Geomorphology, 2021, 378(3): 107615. |
[5] |
王静璞, 刘连友, 沈玲玲. 基于Google Earth的毛乌素沙地新月形沙丘移动规律研究[J]. 遥感技术与应用, 2013, 28(6): 1094. WANG Jingpu, LIU Lianyou, SHEN Lingling. Research of the barchan dunes movement in the Mu Us Sandy Land on Google Earth software[J]. Remote Sensing Technology and Application, 2013, 28(6): 1094. |
[6] |
杨军怀, 董治宝, 刘铮瑶, 等. 库鲁克沙漠风沙地貌与沙丘移动[J]. 中国沙漠, 2019, 39(4): 1. YANG Junhuai, DONG Zhibao, LIU Zhengyao, et al. Aeolian geomorphology and dune migration in the Quruq Desert, China[J]. Journal of Desert Research, 2019, 39(4): 1. |
[7] |
刘宇胜. 阿拉善北部戈壁地区新月形沙丘移动规律研究[D]. 呼和浩特: 内蒙古农业大学, 2018: 19. LIU Yusheng. The Study on the Law of Barchan dunes' Movement in the Gobi Area in northern Alxa[D]. Hohhot: Inner Mongolia Agricultural University, 2018: 19. |
[8] |
宁文晓, 王振亭. 巴丹吉林沙漠东南部典型高大沙丘形态参数分析[J]. 福建农林大学学报(自然科学版), 2018, 47(6): 755. NING Wenxiao, WANG Zhenting. Analysis of morphological parameters of mega-dunes in the southeast of Badain Jaran Desert[J]. Journal of Fujian Agriculture and Forestry University (Natural Science Edition), 2018, 47(6): 755. |
[9] |
李爱敏, 韩致文. 新月形沙丘形态参数与移动速度的关系[J]. 中国沙漠, 2020, 40(1): 29. LI Aimin, HAN Zhiwen. Relationship between moving speed and morphological parameters of barchan dunes[J]. Journal of Desert Research, 2020, 40(1): 29. |
[10] |
李建军, 焦菊英, 曹雪, 等. 柴达木盆地沙丘移动的空间分异及对形态参数的响应[J]. 农业工程学报, 2021, 37(7): 309. LI Jianjun, JIAO Juying, CAO Xue, et al. Spatial regionalization and response to morphological parameters of dune migration in the Qaidam Basin of China[J]. Transactions of the CSAE, 2021, 37(7): 309. |
[11] |
DONG Zhibao, WANG Xunming, CHEN Guangting. Monitoring sand dune advance in the Taklimakan Desert[J]. Geomorphology, 2000, 35(3): 219. |
[12] |
石学刚. 民勤绿洲边缘新月形沙丘移动规律与环境因子的关系研究[D]. 兰州: 甘肃农业大学, 2017: 20. SHI Xuegang. Study on movement regularity of barchan dunes and its relationship with environment factors along the oasis fringe in Minqin[D]. Lanzhou: Gansu Agricultural University, 2017: 20. |
[13] |
王静璞, 王光镇, 韩柳, 等. 毛乌素沙地不同固沙措施下沙丘的移动特征[J]. 甘肃农业大学学报, 2017, 52(2): 54. WANG Jingpu, WANG Guangzhen, HAN Liu, et al. Barchan dunes movement of different sand fixation measures in Mu Us Sandy Land[J]. Journal of Gansu Agricultural University, 2017, 52(2): 54. |
[14] |
严应存, 高贵生, 刘宝康, 等. 青海湖周边地区沙丘移动监测及其气候驱动力分析[J]. 中国沙漠, 2009, 29(4): 617. YAN Yingcun, GAO Guisheng, LIU Baokang, et al. Monitoring of aeolian dunes movement in surrounding Qinghai Lake Areas and climate-driven forces analysis[J]. Journal of Desert Research, 2009, 29(4): 617. |
[15] |
潘凯佳, 张正偲, 梁爱民. 反向沙丘近地层气流变化及其对沙丘形态的影响[J]. 中国沙漠, 2021, 41(2): 1. PAN Kaijia, ZHANG Zhengsi, LIANG Aimin. The characteristics of inner boundary layer airflow over reversing dunes and its effect on dune morphological evolution[J]. Journal of Desert Research, 2021, 41(2): 1. |
[16] |
ZHANG Huiwen, WU Zhen, HU Jing, et al. Numerical simulation of wind field and sand flux in crescentic sand dunes[J]. Scientific Reports, 2021, 11(1). |
[17] |
王宁波, 李生宇, 王海峰, 等. 塔克拉玛干沙漠腹地垄间地上覆沙丘形态的空间变化特征及其成因[J]. 干旱区地理, 2014, 37(1): 89. WANG Ningbo, LI Shengyu, WANG Haifeng, et al. Spatial variation and formation mechanism of dune morphology in the dominant wind direction on interdune corridors of complex ridges in central Taklimakan Desert[J]. Arid Land Geography, 2014, 37(1): 89. |
[18] |
余沛东, 陈银萍, 李玉强, 等. 植被盖度对沙丘风沙流结构及风蚀量的影响[J]. 中国沙漠, 2019, 39(5): 29. YU Peidong, CHEN Yinping, LI Yuqiang, et al. Influence of vegetation coverage on sand flow structure and wind erosion yield with wind tunnel experiment as a case[J]. Journal of Desert Research, 2019, 39(5): 29. |
[19] |
秦海琴, 张萍, 展秀丽, 等. 宁夏河东沙地沙丘冻融过程的时空差异[J]. 冰川冻土, 2020, 42(3): 909. QIN Haiqin, ZHANG Ping, ZHAN Xiuli, et al. Spatial and temporal differences of freeze-thaw processes on dunes in the east desert of the Yellow River in Ningxia[J]. Journal of Glaciology and Geocryology, 2020, 42(3): 909. |
[20] |
宋超, 余琦殷, 王瑞霞, 等. 基于植被覆盖度的宁夏灵武白芨滩自然保护区防风固沙功能时空变化研究[J]. 生态学报, 2021, 41(8): 3131. SONG Chao, YU Qiyin, WANG Ruixia, et al. Spatio-temporal variation of windbreak and sand fixation functions based on vegetation coverage in Baijitan National Nature Reserve, Ningxia[J]. Acta Ecologica Sinica, 2021, 41(8): 3131. |
[21] |
肖涛. 基于高分辨率遥感影像的沙丘形态及演变特征研究: 以响水河中游右岸沙丘群为例[D]. 长沙: 湖南师范大学, 2017: 24. XIAO Tao. A study of dune morphology and its evolution based on high resolution remote sensing data: An example from dunes in the right flank of Xiangshui River middle reaches[D]. Changsha: Hunan Normal University, 2017: 24. |
[22] |
刘建宝, 王乃昂, 程弘毅, 等. 沙丘沙休止角影响因素实验研究[J]. 中国沙漠, 2010, 30(4): 758. LIU Jianbao, WANG Naiang, CHENG Hongyi, et al. Influencing factors of repose angle of sand dunes: An experiment[J]. Journal of Desert Research, 2010, 30(4): 758. |
[23] |
庞营军, 吴波, 贾晓红, 等. 毛乌素沙地风况及输沙势特征[J]. 中国沙漠, 2019, 39(1): 62. PANG Yingjun, WU Bo, JIA Xiaohong, et al. Characteristics of wind regime and drift potential in Mu Us Sandy Land[J]. Journal of Desert Research, 2019, 39(1): 62. |
[24] |
赵纳祺. 乌兰布和沙漠沿黄段冻结对沙丘坡面风积沙运移特征影响[D]. 呼和浩特: 内蒙古农业大学, 2018: 36. ZHAO Naqi. Effects of freezing on migration characteristics of aeolian sand in sand dune along bank of Yellow River in Ulan Buh Desert[D]. Hohhot: Inner Mongolia Agricultural University, 2018: 36. |
[25] |
秦海琴, 展秀丽, 张萍. 宁夏白芨滩固定沙丘浅层土壤水分时空变异[J]. 云南农业大学学报(自然科学), 2020, 35(1): 156. QIN Haiqin, ZHAN Xiuli, ZHANG Ping. Temporal and spatial variation of soil moisture in the shallow dune fixed on Baijitan in Ningxia[J]. Journal of Yunnan Agricultural University (Natural Science), 2020, 35(1): 156. |