1. 引言
灌丛沙堆是流动碎屑颗粒(砂、粉砂、粘土和盐晶粒)受植物丛的阻挡,在植物丛风影区发生堆积并随植物生长而不断扩大的生物–风积地貌类型 [1] [2] [3] ,其高度从几厘米至几米不等,发育大小主要取决于灌丛植物的冠幅和密度。常见的灌丛沙堆有柽柳沙堆、白刺沙堆、锦鸡儿沙堆和油蒿沙堆等 [4] [5] [6] [7] 。在荒漠地区含水分的沙质土地上出现或保留适应性较强的灌丛植物是灌丛沙堆形成的首要条件。我国北方荒漠、半荒漠地区,灌丛沙堆一般发育在洪积扇前缘、沙漠边缘和沙漠腹地中河床两岸、湖盆洼地等水分条件相对较好的地方,是一种广泛存在的生物–风积地貌,自20世纪80年代以来受到关注。迄今,对灌丛沙堆的研究主要涉及其空间格局 [8] 、形成演变 [9] 、沉积特征 [10] 、形态特征 [11] 和流场特征 [12] 以及阻沙效益 [13] [14] [15] 等方面。一般认为,植被、风力和物源等因素与灌丛沙堆形成演化密切相关 [16] [17] [18] 。在不同的生物气候带,其植被、风力状况存在差异,进而影响到灌丛沙堆的形成与发育。极端干旱地区的库姆塔格沙漠以流动沙丘为主,在沙漠腹地的干涸河床和阿奇克谷地也有灌丛沙堆的发育,对其研究鲜有报道。在对库姆塔格沙漠考察时,受恶劣的自然环境限制,仅对个别灌丛沙堆做了观察描述和采样,未开展该沙漠范围内灌丛沙堆的系统调查,但本文旨在通过测定研究区典型灌丛沙堆及周边可能物源沉积物粒度特征,比较灌丛沙堆周边沉积物与灌丛沙堆沉积物的异同性,探讨灌丛沙堆沙物质的来源,为该区灌丛沙堆的发育和演化提供科学数据。
2. 材料与方法
2.1. 研究区概况
研究区位于库姆塔格沙漠西北(图1),罗布泊东南,在地质构造上属塔里木盆地的东南坳陷(500~1000 m)东延部分 [19] ;在库姆塔格沙漠风沙地貌区划分中 [20] 属罗布泊南–阿奇克谷地湖积平原区的罗布泊以南灌丛沙丘亚区。宽广的丘间地主要分布泡泡刺(Nitraria sphaerocarpa)沙堆(图1)和雅丹体,靠北地区分布有柽柳(Tamaricaceae)沙堆 [21] 。因丘间地砂土含盐量较高,几乎无植被覆盖。由于本区处于小泉沟沟口的洪积(倾斜)平原上,受季节性洪水影响,在丘间地表层可见大小不等的砂砾石覆盖。据位于库姆塔格沙漠西部的测风站(39˚46'07''N, 90˚57'19''E)观测数据显示,该地年平均风速为3.7 m∙s−1,年起沙风占总观测结果的26.8%,年风向和起沙风向均为东北和东南两组风向,为钝双峰型风况,两组风向占全年风向的73.73%,其中东北风占36.46%,东南风占37.27% [22] 。
(A) 研究区位置(B) 研究区范围(C) 灌丛沙堆分布区(D) 典型灌丛沙堆
Figure 1. The nebkhas in the Kumtagh Desert
图1. 库姆塔格沙漠的灌丛沙堆
2.2. 研究方法
2015年10月考察期间,在库姆塔格沙漠选取典型泡泡刺沙堆(图1(D)),地理坐标为39˚55'12''N、91˚06'29''E,沙堆高约4 m,上覆植被株高10~40 cm,盖度70%左右。在该沙堆表层0~5 cm的深度采集代表样品。另外,为探究灌丛沙丘沉积物来源,分别在灌丛沙堆周围的丘间地、位于灌丛沙堆东北部雅丹体中自然出露的埋藏古风成砂、西部的流动沙丘和丘间地、东部和东南部丘间地等不同地貌单元进行样品采集,除雅丹体中的埋藏古风成砂外,其余样品均在20 × 20 cm范围内采集0~10 cm的表层沙样。
粒度分析采用1/3Φ间隔筛析法,粒级划分采用乌登–温特沃斯粒级标准 [23] ,即砾石(<−1Φ)、极粗沙(−1~0Φ)、粗沙(0~1Φ)、中沙(1~2Φ)、细沙(2~3Φ)、极细沙(3~4Φ)、粉沙(4~9Φ)、黏土(>9Φ)。粒度参数计算使用矩法 [24] ,并采用矩法粒度参数进行分级,参数有平均粒径(Mz)、分选系数(σ)、偏度(SK)和峰度(Kg),最后将数据结果在SPSS中进行统计分析,使用Origin软件给予图示化。
3. 结果与分析
3.1. 粒度组成与参数特征
粒度是沉积物特征的重要表征之一,其组成及分布受物质来源、搬运方式、沉积环境的综合影响。为了分析和对比库姆塔格沙漠灌丛沙堆与潜在沙源的粒度组成和参数特征,按照乌登–温特沃斯粒级标准统计现有粒度数据,在表1中给出粒度组成及其参数值。
从表1中可以看出,库姆塔格沙漠灌丛沙堆表层沉积物的粒级以极细沙为主,约占47.52%,表明极细沙是灌丛沙堆沉积物的主要组成部分,细沙和中沙含量次之,约占28.45%和14.34%;其他粒级含量极少,在0.04%~7.02%之间。灌丛丘间地的粒度组成以中沙和极细沙为主,含量分别为22.85%和22.97%。同灌丛沙堆沙粒级配相似的有西部丘间地和东南部丘间地,其粒度组成同样以极细沙为主,含量分别为35.16%和29.73%,且二者的沙粒级配在中粒组分和细粒组分上同灌丛沙堆均有一定的相似性。其余样点地表沉积物粒度含量差异明显。东部丘间地以粗沙含量为主,西部流动沙丘以中沙为主,而雅丹体中埋藏古风成砂以细沙为主。据此,通过各样点粒度组成可以发现,灌丛沙堆中的细粒组分所占比重最高。表明在该地风能条件下,各丘间地沙物质中的细沙和极细沙物质最易吹蚀,并被灌丛沙堆所拦截,小于极细沙粒级颗粒物质则易被远距离输送,大于细沙粒级颗粒物质则不易被吹蚀与搬运。
根据矩法公式对所有样品的粒度特征参数进行计算(表1)。结果表明:平均粒径(Mz)中,灌丛丘间地、西部流动沙丘、西部丘间地、东部丘间地和东南部丘间地的优势粒级为中沙,范围介于1.59~2.07Φ,其中灌丛丘间地的平均粒径最粗;灌丛沙堆和雅丹体中埋藏古风成砂均以细沙为主,但灌丛沙堆的平均粒径细于雅丹体中埋藏古风成砂;各样点的分选系数(σ)根据Collias矩法分类标准 [24] 可知,除灌丛沙堆为中等分选外(0.70~1.00),其他样点的分选均较差(1.00~2.00),说明灌丛沙堆沉积物颗粒组成的均一程度好于其他样点;偏度(SK)除灌丛丘间地、东部丘间地和东南部丘间地为“近对称”分布外,其余样点均为负偏,且其分布趋势与平均粒径(Mz)相关;灌丛沙堆和雅丹体中埋藏古风成砂的峰态(Kg)为中峰态,其余样点沙粒的峰态均为宽峰态,说明灌丛沙堆和雅丹体中埋藏古风成砂粒度分布相对集中。从粒度组成和粒度参数特征可以看出,灌丛沙堆的沉积物平均粒径细于其他样点,且分选性中等,但具有负偏和中峰态的特点。其余样点地表沉积物中,西部丘间地、东南部丘间地表层沉积物的粒度组成上同灌丛沙堆有一定的相似性,但雅丹体中埋藏古风成砂的粒度参数与灌丛沙堆有一定的相似性。
3.2. 粒度频率曲线和概率累积曲线
为进一步探讨库姆塔格沙漠灌丛沙堆的沙源及粒度特征,对各样点沉积物颗粒粒度频率曲线与灌丛沙堆粒度频率曲线进行对比分析(图2),以反映灌丛沙堆的粒度频率特征。
![](Images/Table_Tmp.jpg)
Table 1. The grain size composition and parameter of surface sediments in different geomorphological locations
表1. 研究区采样点地表沉积物粒度组成与参数特征
![](//html.hanspub.org/file/16-1770563x11_hanspub.png)
Figure 2. The grain size frequency curves of the nebkhas and other correlative points (different inter-dune lowlands and dunes)
图2. 灌丛沙堆及相关样点(各丘间地与各沙丘)粒度频率曲线
粒度频率曲线可以显示粒度分布的峰态和趋势。从图2中可以看出,灌丛沙堆的粒度分布以4个众数、负偏态、1个明显主峰为特征,且与4个丘间地样点的沉积物在粒级>1Φ范围内的粒度分布趋势表现出明显的一致性,这5个样点沉积物粒度的主峰值大多出现在3.7~4Φ之间,属极细沙组分。灌丛沙堆与雅丹体中埋藏古风成砂沉积物在粒级>1Φ范围内的粒度分布趋势较为相似,而与西部流动沙丘差异明显;在粗沙(0~1Φ)至细沙(2~3Φ)范围(1~3Φ),雅丹体中埋藏古风成砂所对应的粒度频率含量高于灌丛沙堆,但在细沙(2~3Φ)至极细沙(3~4Φ)范围,灌丛沙堆所对应的粒度频率含量高于雅丹体中埋藏古风成砂;灌丛沙堆与雅丹体中埋藏古风成砂样品的主峰值主要出现在3.7~4Φ之间,属极细沙组分。这说明在>1Φ粒级范围内的沙物质为灌丛沙堆和周围样点地表沙物质的主要粒径组分,灌丛沙堆中极细沙含量所占比重高于其他样点。据此,反映出灌丛沙堆的沙物质来源很可能是其丘间地沙物质和雅丹体中埋藏古风成砂。
概率累积曲线可以揭示沉积物与搬运营力间的关系,也可以将碎屑组分中含量较少的粗、细尾部特征放大,便于沉积成因分析 [25] 。灌丛沙堆的概率累积曲线表现出“滚动–多跳跃–悬浮”式组合的粒度特征(图3)。其中,滚动组分(0~0.67Φ)含量极少所占比重约1.10%,可能与其形成过程中短时或瞬时强风将粗粒物质吹蚀沉降于灌丛沙堆有关;跳跃组分(1~3.33Φ)含量占58.40%,悬浮组分(>3.33Φ)含量占40.40%,二者合计所占比重98.80%。由此可知,灌丛沙堆沉积物主要以跳跃组分和悬浮组分为主。
概率累积曲线的斜率可以反映特定粒度组分样品的分选性。图3中显示,滚动组分的斜率为1.38,
![](//html.hanspub.org/file/16-1770563x12_hanspub.png)
Figure 3. The grain size probability cumulative curves of the nebkhas in the Kumtagh desert
图3. 库姆塔格灌丛沙堆沉积物粒度概率累积曲线
跳跃组分的斜率为0.81,悬浮组分的斜率为2.70。据此说明悬浮组分的分选性最好,即灌丛沙堆中细沙(2~3Φ)和极细沙(3~4Φ)的分选性较其他粒级组分好。粒度概率累积曲线表现出的这种特征,与研究区风力 [22] [29] 适中有关,即在中风能环境下以跃移搬运和悬移搬运方式为主。
3.3. 灰色关联度分析
灰色关联度分析是一种定量化的比较分析方法,适用于处理多因素非典型数据 [26] 。将研究区内各样点的粒度组成与参数特征同灌丛沙堆各粒度组成与参数特征作为相应的因子参数序列进行关联度分析。结合表1可知砾石组分和极粗沙组分在灌丛沙堆中所占比例很低,不宜做为参考因子与比较因子进行比较,因此予以剔除。将其余各样点粒度组成与参数数据作为关联因子,运用无量纲化方法进行规范标准化处理,将其规范为极差化数据并建立矩阵(表2)。
表2中X0至X6依次为灌丛沙堆、灌丛丘间地、雅丹体中埋藏古风成砂、西部流动沙丘、西部丘间地、东部丘间地和东南部丘间地。其中X0为参考数列,Xi为比较数列。根据所得数据将参考数列X0与各比较数列Xi进行求差序列,并结合公式(1)求出关联系数(表3),而后按公式(2)计算求出其关联度。
(1)
式中,r为分辨系数,一般取r = 0.5。
(2)
由关联系数经运算可得各样点对灌丛沙堆的关联度,即各样点对于灌丛沙堆的等权关联度。由于各影响因子对于整体贡献各不相同,为客观评价各样点同灌丛沙堆的相关性,需运用加权关联度。为此,将各因子根据不同的重要性赋予权重WK,依据公式(3)运算得到加权关联度。
![](Images/Table_Tmp.jpg)
Table 2. Dimensionless data of correlative factors
表2. 关联因子的无量纲化值
![](Images/Table_Tmp.jpg)
Table 3. Grey correlative coefficients of the different sample points
表3. 各样点灰色关联系数
(3)
运算结果为:γ1 = 0.47,γ2 = 0.73,γ3 = 0.52,γ4 = 0.63,γ5 = 0.45,γ6 = 0.54,γ2 > γ4 > γ6 > γ3 > γ1 > γ5,即各样点与灌丛沙堆的关联度大小依次为雅丹体中埋藏古风成砂>西部丘间地>东南部丘间地>西部流动沙丘>灌丛丘间地>东部丘间地。不难看出,雅丹体中埋藏古风成砂与灌丛沙堆关联度最高,其次为西部丘间地和东南部丘间地,这与粒度分析所体现出来的相似性较为一致。据此,说明雅丹体中埋藏古风成砂、西部丘间地及东南部丘间地沉积物都可能是灌丛沙堆的物质来源。
4. 讨论
灌丛沙堆发育的影响因素主要有沙源、风况、植被,其中灌丛植被为灌丛沙堆形成的首要条件。地处极端干环境条件下的库姆塔格沙漠也不例外,灌丛沙堆的形成首先与其表面覆盖泡泡刺(Nitraria schoberi)植物且盖度较大有关。灌丛植物的生长增加了沙堆表面的粗糙度,有利于拦截通过风力运移的细颗粒物质,并使其沉积于灌丛沙堆表面。Leenders J K等 [27] 在萨赫勒地区通过野外观测灌丛对风沙拦截的研究结果证实了这种现象,李志忠等 [28] 对灌丛沙堆风洞模拟实验的结果也与之类似。王训明等 [12] 的研究认为较低的风能环境有利于灌丛沙堆的发育,对库姆塔格沙漠进行野外风沙活动观测结果显示 [22] [29] ,风况特征是决定风沙地貌的动力因子,研究区属中风能环境,风向变率为中等变率,促使灌丛沙堆分布区丘间地及其周围地表沙物质颗粒在中等风力条件下通过跃移和悬移运移方式堆积发育灌丛沙堆。
荒漠半荒漠地区能否发育为灌丛沙堆其沙源是关键因素之一,且灌丛沙堆沙物质来源多为短距离风力搬运产物 [5] 。因此,对灌丛沙堆邻近各样点进行了粒度分析,并在进一步对粒度组分及参数特征分析的基础上,得出灌丛沙堆沙物质组成与西部丘间地、东南部丘间地的沙物质颗粒组成上具有相似性,其粒度参数特征与雅丹体中埋藏古风成砂有相似性。这似乎暗示了灌丛沙堆沙物质除部分来源于其周围丘间地外,仍有部分来源于雅丹体中埋藏古风成砂。有学者对新疆策勒绿洲沙漠过渡带、南非Molopo盆地和非洲博茨瓦纳地区灌丛沙堆的研究表明 [30] ,灌丛沙堆的沙物质来源以周围丘间地为主。由于库姆塔格沙漠以流动沙丘为主,丘间地的细颗粒物质经长期吹蚀在表层残留厚度极薄(3~8 mm不等)的中粗沙粒级的粗颗粒,并覆盖下层的细沙和极细砂物质免遭进一步吹蚀起到保护作用,这种“保护作用”也称之为“保护理论”在董治宝、钱广强等 [22] 的研究中予以证实。但应用灰色关联相关性分析表明,在加权灰色关联度中与灌丛沙堆关联度最高的为雅丹体中埋藏古风成砂,不难理解,这是由于灌丛沙丘所处的上风向有崩塌垮落的破碎雅丹体,使得埋藏古风成砂出露被吹蚀运移至灌丛阻拦沉降原因所致。显然,库姆塔格沙漠发育的灌丛沙堆,沙物质来源主要是其周围丘间地及雅丹体中的埋藏古风成砂。
5. 结论
通过对库姆塔格沙漠灌丛沙堆及其丘间地沉积物和雅丹体中埋藏古风成砂粒度组成与参数特征的比较,并运用灰色关联度进行定量分析研究,对灌丛沙堆的可能物质来源获得以下认识:
1) 库姆塔格沙漠灌丛沙堆的优势粒级为极细沙,平均粒径为2.92Φ,分选性中等,偏度为负偏,峰态为中等峰态;西部丘间地和东南部丘间地在粒度组成上同灌丛沙堆有相似性;雅丹体中埋藏古风成砂的粒度参数与灌丛沙堆也有一定的相似性。这反映出灌丛沙堆的沙物质来源除丘间地沙物质外,雅丹体中埋藏古风成砂也是沙源之一。
2) 库姆塔格沙漠灌丛沙堆周围丘间地、流动沙丘和雅丹体各采样点同灌丛沙堆的加权灰色关联度大小依次为雅丹体中埋藏古风成砂 > 西部丘间地 > 东南部丘间地 > 西部流动沙丘 > 灌丛丘间地 > 东部丘间地,表明丘间地沙物质、流动沙丘及其雅丹体中埋藏古风成砂都可能是灌丛沙堆的物质来源。
3) 灌丛沙堆的物质来源是影响其发育的因素之一,雅丹体中埋藏古风成砂与灌丛沙堆的相关性较大,表明埋藏古风成砂是灌丛沙堆的物源之一,这是特定环境所造成的,即与上风向崩塌垮落的破碎雅丹体出露埋藏古风成砂被吹蚀运移至灌丛阻拦沉降有关。
基金项目
国家自然科学基金项目(编号:41271030);国家科技基础性工作专项(编号:2012FY111700)资助。