植物生态学报  2016 , 40 (10): 1003-1014 https://doi.org/10.17521/cjpe.2016.0045

Orginal Article

内蒙古荒漠草原放牧强度对风沙通量和沉积物粒径的影响

李永强12**, 李治国1**, 董智3, 王忠武1, 屈志强1, 韩国栋1***

1内蒙古农业大学生态环境学院, 呼和浩特 010018
2土肥资源高效利用国家工程实验室/山东农业大学资源与环境学院, 山东泰安 271018
3山东农业大学林学院/泰山森林生态站, 山东泰安 271018

Effects of grazing intensity on windblown sediment mass flux and particle size distribution in the desert steppe of Nei Mongol, China

LI Yong-Qiang12*, LI Zhi-Guo1*, DONG Zhi3, WANG Zhong-Wu1, QU Zhi-Qiang1, HAN Guo-Dong1**

1College of Ecological and Environmental Sciences, Inner Mongolia Agricultural University, Hohhot 010018, China
2National Engineering Laboratory for Efficient Utilization of Soil and Fertilizer Resources
College of Resources and Environment, Shandong Agricultural University, Taian, Shandong 271018, China
and 3Forestry College of Shandong Agriculture University, Taishan Forestry Ecological Station, Taian, Shandong 271018, China

版权声明:  2016 植物生态学报编辑部 本文是遵循CCAL协议的开放存取期刊,引用请务必标明出处。

基金资助:  国家自然科学基金(31260123)、国家公益行业专项计划项目(201003019和200903060)和教育部草地资源重点实验室和创新团队支持项目。

展开

摘要

以内蒙古短花针茅(Stipa breviflora)草原为研究对象, 在放牧11年的样地布设BSNE集沙仪, 通过2013年4月到2015年4月4个采样期对短花针茅草原放牧强度与风沙通量的定量关系及风蚀物粒度特征进行了研究。结果表明: (1)放牧强度对风沙通量有显著影响(p < 0.05)。随着放牧强度增强, 风沙通量逐渐增加; 风沙通量对放牧强度的响应存在季节差异, 生长季的日平均风沙通量小于非生长季, 且生长季不同放牧强度间风沙通量差异较大, 其中重度放牧区风沙通量是对照区的5倍, 而非生长季不同放牧强度间风沙通量差异较小, 重度放牧区是对照区的1.7倍; (2)生长季不同高度风沙通量的自然对数和垂直高度之间存在较好的二项式拟合结果, 而非生长季风沙通量的垂直分布没有明显规律; 放牧强度间通量的垂直分布差异主要表现在0-50 cm高度; 在同样高度, 放牧强度大, 垂直通量也大; (3)风蚀沉积物中粒径≤250 μm颗粒富集度随垂直高度变化而变化; 风沙沉积物中, ≤250 μm的颗粒占沉积物总量的85%以上, 沉积物中粒径≤50 μm的颗粒有明显的富集作用, 且随垂直高度增加富集比逐渐增加; 粒径为125-250 μm和50-125 μm的颗粒表现为随垂直高度增加富集比降低的趋势, 且125-250 μm粒径的富集比显著小于50-125 μm颗粒的富集比; 荒漠草原放牧强度对风沙通量有不同程度的影响, 放牧强度越大, 风蚀越重, 风蚀作用强化了重度放牧对草地退化的影响作用。

关键词: 放牧强度 ; 风蚀 ; 风沙沉积物 ; 风沙通量 ; 荒漠草原

Abstract

Aims Grazing activities degrade soil aggregates, reduce vegetation coverage and affect the amount of deposited material, and make the land more vulnerable to wind erosion. Although livestock increase was considered as the main issue leading to the degradation, only very few studies have quantitatively investigated the relationship between grazing and soil erosion. The relationship between different stocking rates and sediment flux, and sediment soil particle was studied to reveal the mechanism of different grazing intensities on soil erosion process, to provide basic parameters for grazing optimization in the Stipa breviflora desert steppe. Methods In the Stipa breviflora desert steppe research area, BSNE collecting sand boxes were set in the randomly distributed paddock experiment sites for 11 year with different grazing intensities (0.15、0.30、0.45、0 sheep·hm-2·month-1, corresponding to light grazing LG, moderate grazing MG, heavy grazing HG and control CK, respectively). The quantitative relationship between grazing intensity and sediment flux, and the characteristics of sediment soil particle were conducted in four sampling periods through 2 years (April 2013 to April 2015).Important findings (1) Grazing intensity had a significant effect on the sediment flux (p< 0.05), and the sediment flux increased with the increase of grazing intensity. The response of sediment flux to grazing intensity was variable with season. The daily average sediment flux (13.12 g·m-1·d-1) during the period of April to October was smaller than that from October to April (18.74 g·m-1·d-1). The sediment flux difference of different grazing intensities was greater from April to October, with the 5 times daily average sand flux in the heavy grazing paddock that in the control. The average sediment flux difference of different grazing intensities was small from October to April. (2) The relationship between the natural logarithm of sediment flux at different height and the vertical height had a better binomial fitting from April to October, and there was no obvious regular pattern about flux vertical distribution from October to April, and the vertical flux difference of grazing intensities was mainly expressed in 0-50 cm layer. (3) Sand sediment particle ≤250 μm accounted for more than 85% of the total sediment, the sand sediment particle of ≤50 μm) size was significantly enriched, and the enrichment ratio increased with the increase of vertical height. The enrichment ratio of 125-250 μm particle and 50-125 μm particle decreased with the increase of vertical height, and the enrichment ratio of 125-250 μm particle was smaller than that of 50-125 μm particle (p< 0.05). Therefore grazing intensity had different influence on the sand flux in Stipa breviflora desert steppe, the greater the grazing intensity, the heavier the wind erosion was, and the effect of grazing intensity on grassland was enhanced by wind erosion.

Keywords: grazing intensity ; wind erosion ; windblown sediments ; sediment flux ; desert steppe

0

PDF (991KB) 元数据 多维度评价 相关文章 收藏文章

本文引用格式 导出 EndNote Ris Bibtex

李永强, 李治国, 董智, 王忠武, 屈志强, 韩国栋. 内蒙古荒漠草原放牧强度对风沙通量和沉积物粒径的影响. 植物生态学报, 2016, 40(10): 1003-1014 https://doi.org/10.17521/cjpe.2016.0045

LI Yong-Qiang, LI Zhi-Guo, DONG Zhi, WANG Zhong-Wu, QU Zhi-Qiang, HAN Guo-Dong. Effects of grazing intensity on windblown sediment mass flux and particle size distribution in the desert steppe of Nei Mongol, China. Chinese Journal of Plant Ecology, 2016, 40(10): 1003-1014 https://doi.org/10.17521/cjpe.2016.0045

土壤风蚀(wind erosion)普遍存在于干旱半干旱区, 影响1/3以上的陆地面积(Pimentel & Kounang, 1998; Li et al., 2009), 为空气沙尘提供物质来源(李少昆等, 2008), 是下游风沙沉降区空气污染的重要来源。在北美的圣地亚哥, 引起呼吸道疾病的PM10污染物中, 有50%来自土壤颗粒, 空气污染程度与风暴频率有直接的关系(Garciachevesich et al., 2014)。在我国, 沙尘是北方许多城市的重要季节性污染物(李万元等, 2012)。风蚀造成异地污染的同时, 在风沙源区带走大量的土壤营养物质, 是荒漠化地区土壤养分变异的重要驱动力(Li et al., 2007)。沙漠化过程中的地表粗化、营养物质丧失以及植被退化等都与土壤风蚀有直接或间接的关系。1952-2001年的统计结果(李万元等, 2012)表明, 内蒙古中西部的强尘沙暴几乎占我国北方沙尘暴总频数的2/3。于国茂等(2011)按照土壤风蚀危险度模型也将内蒙古中西部地区的草原区域评定为风蚀危险型和风蚀强险型。内蒙古中西部地区的草原类型主要是荒漠草原, 是旱生性最强的草原生态系统, 承受着比其他草原类型更为严重的风蚀和放牧双重胁迫(Aubault et al., 2015)。

风蚀率通常用风蚀模数, 即单位面积、单位时间的土壤风蚀量来表示。普遍认为风沙流水平通量(Q)和风蚀率存在定量关系(Li et al., 2007), 单位时间、单位面积土壤沙尘释放Fa = K × Q (Fa为垂直通量, K为常数项)的关系被广泛认可(Shao et al., 1993), Q是风蚀研究中的关键环节(Basaran et al., 2011)。完整的植被有助于近地表土壤结构的稳定, 可以有效地减少风沙通量。在干旱、半干旱地区保持一定临界植被盖度对控制土壤风蚀的发生与发展具有重要的意义。在东北地区, 裸露的农田是重要的沙尘来源。保护草原能有效地抑制沙尘排放(Zhang et al., 2015)。植被抑制风蚀的效果可能存在一个盖度阈值。由于风蚀对土壤、植被和气候变化空间格局的敏感性, 量化影响风蚀速率的不同管理措施非常困难(Aubault et al., 2015)。关于土壤风蚀的研究结论大多通过理论推导或风洞模拟实验而获得, 研究对象大多是沙漠分布区的风沙土(董智, 2004; 杨兴华等, 2013), 或者关注于农牧交错区耕作措施对土壤风蚀的影响作用(李晓丽和申向东, 2006)。仅有少量研究关注放牧对草地风蚀的影响, 但侧重于单一放牧强度的水土流失强度和风沙流变化等。不同放牧梯度下草原风蚀过程的野外系统观测研究尚不多见。

短花针茅(Stipa breviflora)草原是内蒙古荒漠草原的重要组成部分。放牧导致短花针茅草原退化、沙化一直是内蒙古荒漠草原面临的主要环境问题。针对荒漠草原生态系统, 学者们开展了较多放牧影响的研究工作(韩国栋等, 2007; 王忠武, 2009; Kemp et al., 2013), 但是迄今我们对放牧强度对荒漠草原风蚀特征的影响并不清楚。为此, 本文作者依托内蒙古荒漠草原野外科学观测试验站建立的长期放牧实验平台, 通过2年野外BSNE集沙仪取样, 开展了不同放牧强度对荒漠草原风蚀过程的影响研究, 主要想回答两个方面的科学问题: (1)不同季节, 荒漠草原风沙通量对放牧强度如何响应?(2)荒漠草原风沙沉积物颗粒粒径的主要分布特征及其对不同放牧强度的响应有何不同?

1 材料和方法

1.1 试验地概况

试验地位于内蒙古自治区乌兰察布市四子王旗王府一队(41.79° N, 111.90° E, 海拔1450 m), 这里春季干旱多风, 夏季炎热, 年降水量311.4 mm, 降水主要集中在6-9月, 6-9月降水量占全年降水量的70%以上, 年平均气温3.4 ℃, 土壤类型为栗钙土。群落类型为短花针茅+冷蒿(Artemisia frigida)+无芒隐子草(Cleistogenes songorica)。植物群落主要由20多种植物组成, 建群种为短花针茅, 优势种为冷蒿和无芒隐子草。试验区草地在2004年随机区组放牧试验前一直沿用简单的连续放牧制度, 并且放牧牲畜超过了草地载畜能力, 多年的过度利用使这里的草地受到严重破坏, 草地生产力低下。由于气候干旱和强烈的人为活动, 在围栏放牧之初试验区草地已出现了一定程度的退化。

1.2 试验设计

放牧控制实验平台于2004年6月开始运行。该平台占地面积约50 hm2。实验处理采用随机区组放牧试验, 放牧区共分成3个区组, 每个区组内设有3个载畜率水平和一个不放牧的对照, 分别为轻度放牧(LG)、中度放牧(MG)、重度放牧(HG)和对照(CK) 4个处理, 随机排列。每个处理设为3次重复(LG1、LG2、LG3; MG1、MG2、MG3; HG1、HG2、HG3; CK1、CK2、CK3), 每个放牧小区面积为4.4 hm2, 共设12个放牧小区。家畜类型为绵羊, 放牧只数分别为0、4、8、12只, 对应的载畜率分别为: 0 (CK)、0.15 (LG)、0.30 (MG)、0.45 (HG) sheep·hm-2·month-1。每年6-12月为放牧期。期间每天将家畜赶入放牧区让其自由采食, 夜间赶回畜圈休息。

1.3 取样和测定方法

在每个放牧小区中间, 选择能代表对应放牧强度的典型地段, 布置5 m × 20 m的小样区, 用于安装BSNE沙尘采样器。每套BSNE沙尘采样器由不同高度的取样器构成(Fryrear & Saleh, 1993), 取样器配有风叶和可旋转轴心, 安装在支撑杆的不同高度; 采样器的高度可以调整。每个采样器集沙口宽2 cm, 高5 cm。含沙风进入集沙口后, 沙粒被捕获, BSNE沙尘采样器的集沙效率在90%以上(Shao et al., 1993)。BSNE沙尘采样器按照4 m间隔沿线安置在每个5 m × 20 m小样区, 每个小区布置4套BSNE沙尘采样器, 排列与盛行风垂直。每套系统(BSNE沙尘采样器+中心支撑金属管)的中心支撑金属管被垂直钉入地面以下至少50 cm, 保证在实验期间不发生沉降。4套BSNE沙尘采样器, 除1根1.5 m高的BSNE沙尘采样器支撑杆安装6个高度集沙盒(0.1 m、0.3 m、0.5 m、0.7 m、1.0 m、1.2 m)外, 其他3个1.2 m高的支撑杆安装3个不同高度集沙盒(0.1 m、0.5 m、1.0 m), 每个小区15个集沙盒, 12个小区, 共180个集沙盒。

BSNE沙尘采样器+中心支撑金属管在2013年春季风季开始之前安装完毕, 全部集沙时间包含两个完整年度, 分别是2013年4月-2014年4月和2014年4月-2015年4月(表示为2013-2014年和2014-2015年, 下同)。每年两次取样, 时间在10月中旬(采样期4-10月, 对应于生长季, 下同)和次年4月中旬(采样期上年10月至当年4月, 对应于非生长季, 下同), 每次取样时准确记录日期, 两次间隔日期用于计算日平均通量(4次取样, 间隔日期分别是191 d、198 d、187 d、176 d); 由于第一个采集期集沙盒掉落地面(可能是风或家畜触动导致), 在2013年10月共收集了178个沙尘样品; 重新固定集沙盒后, 后期每次取样获得180个样品, 4次取样, 总共采集710个沙尘样品。每次取样时, 在5 m × 20 m小样区内采集土壤表层和结皮样品。每次取20个土壤表层样品(0-5 cm), 12个小区共240个样品; 4次取样共480个土壤样品。安装BSNE沙尘采样器的同时, 在CK2区安装野外小型气象站(AG1000农田草地小气候站, 北京天诺基业科技有限公司, 北京), 用于测定试验期间温度、降水量、空气湿度、土壤温/湿度、2.5 m高的风速风向、太阳辐射参数等指标, 所有气象指标每10 min测定一次。

1.4 计算和统计分析

风沙流水平通量(Q)可以用单位时间内垂直于风向的固定高度单位宽度的风沙质量来表示, Q由不同高度风沙流q(z)构成, q(z)有不同的拟合方程(Fryrear & Saleh, 1993)。本试验按照Shao和Raupach (1992)的方法进行拟合, 不同高度单位面积水平通量q(z)满足下面的方程(1)关系:

qz=ce(az2+bz)(1)

这里a、bc是拟合参数, z是集沙口高度(m)。先将方程两边取自然对数, 即 lnq(z) = az2 + bz + lnc, 上述公式可以通过SPSS进行多项式拟合。通过二次多项式函数建立不同垂直高度水平通量q与集沙口高度z的关系; 通过拟合的多项式方程求得常数项abc, 完成q(z)方程的建立。

需要进一步对q(z)方程进行从地表到1 m高度的积分计算求得Q。根据积分方程(2)可以求得1 m高度1 m宽度截面范围内的水平通量Q

Q=0lmq(z)dz(2)

用Excel 2003和Origin 9.0软件构建图表, 采用Matlab进行积分计算。采用IBM SPSS V21.0软件对不同载畜率草原的风沙通量差异进行单因素方差分析, 多重比较采用最小显著差异法(LSD)进行, p < 0.05差异显著, 参数以平均值±标准偏差(mean ± SD)表示。

2 结果和分析

2.1 两个采样年度气象要素特征

研究区2013-2014年和2014-2015年两个采样年度总降水量分别为289 mm和210 mm, 降水量年度差异较大, 分别是多年平均降水量的91%和69% (图1); 按照生长季和非生长季进行划分, 2013-2014年生长季和非生长季降水量分别为266.9 mm和22.1 mm, 分别比2014-2015年对应时期高71.4 mm和 7.6 mm; 生长季降水量差异主要表现在6-8月, 2013年6-8月降水总量比2014年同期高120 mm。由于大气相对湿度和降水量的相关性, 大气湿度变化与降水量有类似的变化规律。

两个采样年度主导风向季节变化特征基本类似(图2), 4-10月以西南风向为主, 10月至次年4月以西风为主; 2013-2014年大于7 m·s-1的总时间为3803 h, 2014-2015年为3798 h, 分别占到对应时段的42%和44%。但大于7 m·s-1风速覆盖时段在季节间存在差异(图2), 2013-2014年生长季平均风速大于7 m·s-1的时段达到2292 h, 占到总时段的50%, 而非生长季大于7 m·s-1的时段共计1611 h, 占统计时段的34%; 2014-2015年的生长季和非生长季, 大于7 m·s-1风速段大约占对应时段的比例都是43%。

2.2 风沙流通量变化

根据方程(2)计算的不同放牧强度水平输沙通量见图3。年度、季节、载畜率的多因素单变量方差分析(年度和季节间风沙通量差异采用t检验方法)表明载畜率和季节变化对平均风沙通量影响达到极显著水平(p < 0.05), 并且二者之间的交互效应对通量也有显著影响(p < 0.05); 随放牧强度增加, 风沙通量呈显著增加趋势, LG、MG、HG的风沙通量分别是对照区的1.4、2.2和2.5倍。多重比较表明, 两年的HG和MG显著高于CK和LG (p < 0.05), 其他处理间无显著差异。非生长季平均风沙通量为18.74 g·m-1·

d-1, 是生长季(13.12 g·m-1·d-1)的1.4倍。风沙通量对放牧强度的响应存在季节差异; 在生长季(4-10月)风沙通量随放牧强度增强而增加的趋势比非生长季(10月至次年4月)更明显(图3), HG区风沙通量是对照区的5倍, 而非生长季不同放牧强度间风沙通量差异缩小, HG区仅是对照区的1.7倍; 2014-2015年日平均风沙通量16.46 g·m-1·d-1, 稍高于2013-2014年的15.53 g·m-1·d-1, 但年度间差异不显著(p = 0.378)。

图1   试验站2013年4月至2015年4月降水量和相对湿度的月变化。

Fig. 1   Monthly variation of precipitation and relative humidity from April 2013 to April 2015 at the experimental site.

图2   不同采样期(2013年4月至2015年4月)风向风速分布。

Fig. 2   Distribution of frequency (%) of wind speed from different directions in different sampling dates (April 2013 to April 2015).

图3   不同放牧强度在生长季和非生长季日平均风沙通量变化(平均值±标准误差, n = 3)。A, 生长季。B, 非生长季。CK, 对照; LG, 轻度放牧; MG, 中度放牧; HG, 重度放牧。不同小写字母表示同一时期不同放牧处理间差异显著(p < 0.05)。

Fig. 3   Mean horizontal mass flux versus different levels of stocking rate in growing-season and non-growing-season in two sampling years (mean ± SE, n = 3). A, Growing-season. B, Non-growing season. CK, control; LG, lightly grazed; MG, moderately grazed; HG, heavily grazed. Different lowercase letters indicate significant differences among the differently stocking rate during the same period (p < 0.05).

2.3 不同放牧强度风沙流结构

风沙通量自然对数与集沙口高度多项式拟合曲线、方程和参数见图4。生长季的拟合效果(R2平均值0.830)好于非生长季(R2平均值0.525)。拟合方程参数统计表明, 特别是2014-2015年的生长季, 各个处理的回归方程均达到显著水平(p < 0.05), 在MG区达到极显著水平(p < 0.001); 而非生长季, 风沙通量对数曲线随垂直高度变化变得更为平缓, 绝大多数处理拟合方程未达到显著水平(p > 0.05)。在不同放牧强度处理下, 拟合曲线的相对位置不同, 在0-50 cm 高度, 从下至上曲线依次为CK、LG、MG、HG。总体来看, 风沙通量自然对数和高度的二次拟合关系在生长季的效果更好, 在植被较好的季节, 风沙流结构有更显著的垂直梯度变化(图4), 而在非生长季, Shao和Raupach (1992)的风沙通量指数递减方程拟合效果有限。

图4   不同采样期风沙流水平通量对数随高度的变化。A, 2013-2014年生长季。B, 2014-2015年生长季。C, 2013-2014年非生长季。D, 2014-2015年非生长季。CK, 对照; LG, 轻度放牧; MG 中度放牧; HG, 重度放牧。

Fig. 4   Natural logarithm of sediment flux as a function of height during different sampling seasons. A, 2013-2014 growing-season. B, 2014-2015 growing-season. C, 2013-1014 non-growing-season. D, 2014-2015 non-growing-season. CK, control; LG, lightly grazed; MG, moderately grazed; HG, heavily grazed.

2.4 不同放牧强度风沙沉积物粒径特征

为了便于分析, 选取两年度生长季对照区和HG区不同高度(10 cm、30 cm、50 cm、100 cm)沉积物粒径特征进行比较(图5)。对照区和HG区风沙沉积物所有样品大于250 μm的粗砂粒含量均未超过15%; 对照区10 cm高度集沙口捕获的风沙沉积物中大于250 μm的粗砂粒平均含量为10.91%, 而在100 cm高度集沙口处, 大于250 μm的粗砂粒含量下降为7.57%, 垂直高度对大于250 μm的粗砂粒含量影响显著(p < 0.05); 在HG区, 不同高度的大于250 μm的粗砂粒含量差异不显著(p > 0.05), 平均含量为14.62%。这表明HG区的风沙沉积物中粗颗粒含量更多。小于50 μm的粉粒和黏粒含量在不同放牧强度和垂直高度也有明显的变化: 对照区10 cm高度的平均含量为44.9%, 而在100 cm高度平均含量为61.0%, 但随着高度增加, 小于50 μm的粉粒和黏粒含量并非线性增加; 对照区含量峰值出现在50 cm高度, 达到70.20%, 显著高于其他垂直高度的细颗粒含量; HG区的细颗粒含量从近地表10 cm高度的40.02%增加到100 cm高度的49.92%, 细颗粒峰值出现在70 cm高度(59.89%)。这表明HG区风沙沉积物中细颗粒含量普遍较低, 且随垂直高度的增加而变化的幅度较小。

图5   对照(CK)和重度放牧(HG)区生长季(10月采样期)不同集沙高度的风沙沉积物粒径分布曲线。A, 2013年10月CK区。B, 2014年10月CK区。C, 2013年10月HG区。D, 2014年10月HG区。10 cm、30 cm、50 cm、100 cm为集沙高度。

Fig. 5   Particle size distributions for sediment samples of the BSNE different heights on control (CK) and heavily grazed (HG) treatments during growing seasons (October sampling period). A, October 2013 CK. B, October 2014 CK. C, October 2013 HG. D, October 2014 HG. 10 cm, 30 cm, 50 cm, and 100 cm indicate heights of the traps.

以对应样地表层土壤为基础, 为详细比较直径(D) ≤ 250 μm颗粒在不同高度的富集比(E), 把D ≤ 250 μm颗粒分为≤ 50 μm、50 μm < D ≤ 125 μm和125 < D ≤ 250 μm三级。图6表明, D ≤ 50 μm的颗粒E值在所有处理和垂直高度上都大于1.0, 位于基准线(E = 1.0)右侧; D为50-125 μm和125-250 μm的颗粒E值都小于1.0, 位于基准线左侧; 随着垂直高度增加, 各处理D ≤ 50 μm的颗粒E值均值逐渐增加, 从10 cm高度的1.76增加到120 cm高度的3.14。富集比最高的处理是HG区的120 cm高度, 达到3.36; 就不同放牧强度而言, D ≤ 50 μm的颗粒E值随垂直高度增加的程度(斜率)在MG和HG区比较稳定, 而CK和LG中, 50 cm高度以上的E值随高度变化较小, CK和LG的平均E值分别为3.02和2.77。在D ≤ 250 μm颗粒的3个分级中, E值最小的是125-250 μm颗粒, 平均值为0.47, 50-125 μm颗粒为0.70; 125-250 μm和50-125 μm颗粒的E值呈现基本与D ≤ 50 μm颗粒相反的变化趋势, 随高度增加, E值逐渐减小, 但变化幅度较小。

图6   不同放牧处理风沙沉积物中直径(D) ≤ 250 μm颗粒富集比随高度的变化(富集比=沉积物中某粒径含量/土壤表土中对应粒径含量)。CK, 对照。LG, 轻度放牧。MG, 中度放牧。HG, 重度放牧。

Fig. 6   Enrichment factor (E) of particle size ≤ 250 μm in sediments versus different heights of different treatments (enrichment factor = content of particle size in the sediments / content of the relative particle size in the surface soil layer). CK, Control. LG, Lightly grazed. MG, Moderately grazed. HG, Heavily grazed. D250, particle size ≤ 250 μm.

3 讨论

3.1 放牧强度对输沙通量的影响

风蚀率是风蚀研究中最受关注的问题(杜鹏飞和刘孝盈, 2012), 包括利用WEPS (Wagner, 2013)、WEAM (Shao et al., 1996)等不同风蚀预报模型的风蚀量估测(邹学勇等, 2014)、遥感和GIS技术的风蚀动态监测技术(Santini et al., 2010)、元素示踪技术等。由于尺度效应和因子的复杂性(Webb et al., 2014), 在判断具体土地利用模式或者小尺度风蚀强度时, 各种野外实测技术如风蚀盘法和集沙仪等(Fryrear, 1986)仍是风蚀率测定常用的方法。由于风蚀盘法容器材料和土壤性质的差异, 不能完全反映野外土壤复杂的表面状况, 只适合于单次风蚀事件的风蚀量估测。武川县试验站采用风蚀盘法测得10月至次年4月的裸地风蚀量(赵彩霞等, 2005)达到3700 t·km-2, 明显高于相同区域的其他研究结果(孙悦超等, 2007); 集沙仪法测定输沙通量是风沙流结构研究中最常用的方法。许多研究者开发了不同类型的集沙仪, BSNE是近年来用于野外测定风蚀过程常用的集沙仪类型。风蚀作为一个自然过程, 涉及风蚀因子的许多方面, 但土地利用模式对各种风蚀因子的影响是人类管理生态系统风蚀变化的主要原因(邹学勇等, 2014)。放牧导致草地植被盖度降低是影响风沙通量的主要原因之一。在非洲的Kalahari荒漠(Wiggs et al., 1994), 由于放牧等原因降低植被盖度60%时, 近地表风速增加200%; 放牧通过减少植被和破坏土壤表皮, 加快风力侵蚀速度(Belnap et al., 2009)。

本试验所在的长期放牧样地, 从2004年开始实施随机区组放牧试验, 最初不同载畜率间的植被盖度没有明显不同, 但从2007年开始, MG和HG区的盖度已经发生明显变化, 显著低于CK和LG区, 放牧对盖度的累积效应在第三年开始得到体现(王忠武, 2009)。两个采样年度的非生长季平均风沙通量比生长季高41%, 反映了植被盖度季节变化导致的风蚀强度的变化。与CK相比, 生长季的HG区平均风沙通量增加5.1倍; 非生长季, 荒漠草原植物盖度普遍比较低, 在风速条件大致相同的条件下(图2), 各放牧处理的风沙通量都比生长季有明显的增加, 但CK区风沙通量和HG区风沙通量差距缩小。这表明放牧对风沙通量的影响强度与不同季节的降水、盖度等因子有关。整个试验期HG区的平均风沙通量为22.76 g·m-1·d-1。与CK相比, HG导致风沙通量增加2.5倍(图3)。植被盖度和输沙通量之间不是简单的线性关系(Liu & Westphal, 2001), 植被盖度作为重要的风蚀影响因子, 也影响着近地表风速、气流和土壤的接触面积、土壤水分等其他风蚀影响因子(李万元等, 2007), 植被盖度增加有效地降低了土壤对风蚀的敏感性(Shabani et al., 2014)。荒漠草原非生长季不同放牧强度间的盖度差异比生长季间的盖度差异要小。

3.2 风沙流结构和风沙沉积物粒度特征

许多学者在风洞或野外都进行过风沙流结构特征的研究。Shao和Raupach (1992)、Zobeck和Stout (1996)、Namikas (2003)的拟合结果表明, 指数函数能够较好地模拟近地表风沙流结构特征。但风沙流垂直结构变化受沙粒粒径、沙源、下垫面特征、风速等因素的影响, 野外过程比风洞复杂得多, 即使是风洞试验, 结果也有较大差异。Dong等(2006)开展的风洞实验发现, 近床面附近某一高度以下的输沙量小于指数分布预测值; 而Ni等(2003)的风洞实验结果是近床面输沙量相对于指数函数的计算发生正偏差。由此可见, 不同研究者对输沙量垂向分布规律的认识存在争议, 输沙量垂直分布规律的研究结果大致可以分为三种类型: 完全符合指数递减分布、上部分符合指数递减而下部分偏小、上部分符合指数递减而下部分偏大。本研究中, 不同季节风沙流垂向分布的指数拟合方程拟合效果不同: 生长季(4-10月)由于有较好的植被覆盖, 不同垂直高度风速存在明显不同, 风沙流通量指数递减趋势比较明显; 而非生长季(10月至次年4月), 荒漠草原处于枯黄期, 植被盖度普遍偏低, 垂直高度间风速差异较小, 而不同气流强度携带的风蚀颗粒达到的高度也有很大不同, 导致不同垂直高度间风沙流通量累计量差异减小, 垂直分布曲线表现更为平缓。根据拟合曲线的相对位置可以看出, 在0-50 cm高度, 拟合曲线从下到上依次为CK、LG、MG、HG, 表明载畜率增加, 主要增强了0-50 cm层次的通量分配。

Lyles和Tararko (1986)在美国Kansas草原研究发现, 36年的风蚀作用使得表层土壤砂粒含量上升6.5%, 黏粒含量下降7.2%。在澳大利亚南部的牧场也发现, 在20周内风蚀可以使土壤中大于250 μm颗粒的含量增加, 75-210 μm和小于2 μm颗粒的含量减少(Leys & Mctainsh, 1996)。土壤中大于250 μm颗粒的含量逐渐增加, 这与本研究风沙沉积物中以小于250 μm颗粒为主的特征一致; 但对于更精细范围的50-125 μm颗粒, 因为风蚀物中该粒径范围的富集比小于1, 因此也属于风蚀不敏感颗粒。在本研究中各粒径范围, 只是小于50 μm的颗粒富集比大于1, 与澳大利亚南部牧场的敏感颗粒明显不同。也有学者认为, 细颗粒损失量实际是风蚀和大气沉降的净结果, 沉降本身会抵消一部分侵蚀造成的损失。由于细颗粒含有更多的营养物质, 细颗粒损失可能是导致草原养分枯竭、植被退化的一个重要机制(苏永中和赵文智, 2005)。在本研究区附近阴山旱作耕地风蚀颗粒研究(李晓丽和申向东, 2006)表明, 在裸露耕地上, 大于2 mm的颗粒粒径的百分含量占整个颗粒粒径的32%, 认为长期受风蚀影响, 地表已经粗化。本研究各集沙盒的风沙沉积物中未收集到大于2 mm的颗粒, 小于250 μm的颗粒占到85%, 放牧草原与耕地有完全不同的风沙粒径结构。风沙沉积物与土壤0-5 cm粒径对比结果表明, 随着高度增加, 风蚀物颗粒组成逐渐变细。小于250 μm风蚀颗粒富集图(图5)表明, 风蚀作用下, 土壤表土中小于50 μm颗粒最为活跃, 而50-125 μm和125-250 μm颗粒在表土中有相应的聚集现象。

4 结论

放牧强度对内蒙古荒漠草原风蚀过程具有重要影响, 风沙通量随放牧强度的增加呈逐渐增加的趋势, 且风沙通量对放牧强度的响应受季节变化的影响; 过度放牧导致草地植被高度和盖度降低、地表裸露、地表破碎化。季节变化涉及植被、气候、地表状况等的风蚀因子变化。荒漠草原非生长季平均风沙通量显著高于生长季; 放牧强度间的风沙通量差异主要表现在近地表层次。生长季的风沙流结构垂向分布符合指数递减函数, 而非生长季, 指数递减函数拟合效果较差; 内蒙古荒漠草原风沙沉积物以直径小于250 μm的颗粒为主, 小于50 μm的粒径在沉积物中存在明显的富集作用, 随高度增加富集比逐渐增加; 放牧导致细颗粒随垂直高度增加而累计。本文在4种载畜率条件下分析了荒漠草原不同放牧利用强度对风沙通量和沉积物粒径的影响, 为不同放牧情景下的荒漠草原风蚀模型建立提供了基础数据和参考。放牧对植被盖度和侧影盖度(lateral coverage)、地表破碎度、土壤含水量等风蚀因子的影响是导致风沙通量变化的直接原因。特别是放牧作用下株丛结构和群落结构决定的侧影盖度的变化(Li et al., 2014)可能对风蚀过程的影响更大。在未来研究中应更多地考虑放牧和株丛结构的耦合效应, 以便深入揭示荒漠草原风蚀过程对放牧管理的响应。

The authors have declared that no competing interests exist.

作者声明没有竞争性利益冲突.


参考文献

[1] Aubault H, Webb NP, Strong C, Mctainsh GH, Leys JF, Scanlan JC (2015).

Grazing impacts on the susceptibility of rangelands to wind erosion: The effects of stocking rate, stocking strategy and land condition.

Aeolian Research, 17, 89-99.

[本文引用: 2]     

[2] Basaran M, Erpul G, Uzun O, Gabriels D (2011).

Comparative efficiency testing for a newly designed cyclone type sediment trap for wind erosion measurements.

Geomorphology, 130, 343-351.

[本文引用: 1]     

[3] Belnap J, Reynolds RL, Reheis MC, Phillips SL, Urban FE, Goldstein HL (2009).

Sediment losses and gains across a gradient of livestock grazing and plant invasion in a cool, semi-arid grassland, Colorado Plateau, USA.

Aeolian Research, 1, 27-43.

[本文引用: 1]     

[4] Dong Z (2004).

Research on Farmland Wind-Sand Disaster of Oasis and Its Control Mechanism in Ulan Buh Desert. PhD dissertation,

Beijing Forestry University, Beijing. 59-68.(in Chinese with English abstract).
[董智 (2004).

乌兰布和沙漠绿洲农田沙害及其控制机理研究

. 博士学位论文, 北京林业大学, 北京. 59-68.]

[本文引用: 1]     

[5] Dong ZB, Qian GQ, Luo WY, Wang HT (2006).

Analysis of the mass flux profiles of an aeolian saltating cloud.

Journal of Geophysical Research, 111, 1-11.

[6] Du PF, Liu XY (2012).

A review of measurement of wind erosion rate in China research of soil and water conservation.

Research of Soil and Water Conservation, 19, 275-281.(in Chinese with English abstract)
[杜鹏飞, 刘孝盈 (2012).

中国土壤风蚀速率实测研究述评. 水土保持研究,

19, 275-281.]

[本文引用: 1]     

[7] Fryrear DW (1986).

A field dust sampler.

Journal of Soil and Water Conservation, 41, 117-120.

[本文引用: 1]     

[8] Fryrear DW, Saleh A (1993).

Field wind erosion: Vertical distribution.

Soil Science, 155, 230-294.

[本文引用: 2]     

[9] Garciachevesich P, Alvarado S, Neary DG, Valdes R, Valdes J, Aguirre JJ, Mena M, Pizarro R, Jofre P, Vera M, Olivares C (2014).

Respiratory disease and particulate air pollution in Santiago Chile: Contribution of erosion particles from fine sediments.

Environmental Pollution, 187, 202-205.

[本文引用: 1]     

[10] Han GD, Jiao SY, Biligetu, Aodenggaowa (2007).

Effects of plant species diversity and productivity under different stocking rates in the Stipa breviflora Griseb. desert steppe.

Acta Ecologica Sinica, 27, 182-188.(in Chinese with English abstract)
[韩国栋, 焦树英, 毕力格图, 敖登高娃 (2007).

短花针茅草原不同载畜率对植物多样性和草地生产力的影响.

生态学报,27, 182-188.]

[本文引用: 1]     

[11] Kemp DR, Han GD, Hou XY, Michalk DL, Hou FJ, Wu JP, Zhang YJ (2013).

Innovative grassland management systems for environmental and livelihood benefits.

Proceedings of the National Academy of Sciences of the United States of America, 110, 8369-8374.

[本文引用: 1]     

[12] Leys J, Mctainsh GH (1996).

Sediment fluxes and particle grain-size characteristics of wind-eroded sediments in southeastern Australia.

Earth Surface Processes and Landforms, 21, 661-671.

[本文引用: 1]     

[13] Li J, Okin GS, Alvarez L, Epstein HE (2007).

Quantitative effects of vegetation cover on wind erosion and soil nutrient loss in a desert grassland of southern New Mexico, USA.

Biogeochemistry, 85, 317-332.

[本文引用: 2]     

[14] Li J, Okin GS, Epstein HE (2009).

Effects of enhanced wind erosion on surface soil texture and characteristics of windblown sediments.

Journal of Geophysical Research, 114, 1-8.

[本文引用: 1]     

[15] Li J, Okin GS, Tatarko J, Webb NP, Herrick J (2014).

Consistency of wind erosion assessments across land use and land cover types: A critical analysis.

Aeolian Research, 15, 253-260.

[本文引用: 1]     

[16] Li SK, Lu M, Wang KR, Wang X (2008).

Soil erosion of the main ground surface types influences on formation of dust storm in south Xinjiang.

Scientia Agricultura Sinica, 41, 3158-3167.(in Chinese with English abstract)
[李少昆, 路明, 王克如, 王旭 (2008).

南疆主要地表类型土壤风蚀对形成沙尘暴天气的影响.

中国农业科学,41, 3158-3167.]

[本文引用: 1]     

[17] Li WY, Dong ZB, SH, Wang SG, Ao YH (2012).

Spatial distribution and temporal tendency of the dust storm frequency in north and northwest China for the recent nearly half a century.

Climate Change Research Letters, (1), 1-12.(in Chinese with English abstract)
[李万元, 董治宝, 吕世华, 王式功, 奥银焕 (2012).

近半个世纪中国北方沙尘暴的空间分布和时间变化规律回顾

. 气候变化研究快报, (1), 1-12.]

[本文引用: 2]     

[18] Li WY, Sheng ZB, SH, Li YH (2007).

Sensitivity tests of factors influencing wind erosion.

Journal of Desert Research, 27, 78-87.(in Chinese with English abstract)
[李万元, 沈志宝, 吕世华, 李耀辉 (2007).

风蚀影响因子的敏感性试验.

中国沙漠,27, 78-87.]

[本文引用: 1]     

[19] Li XL, Shen XD (2006).

Experimental study on the distribution characteristics of the saltation particle of aeolian sediment in bare tillage.

Transactions of the CSAE, 22, 74-77.(in Chinese with English abstract)
[李晓丽, 申向东 (2006).

裸露耕地土壤风蚀跃移颗粒分布特征的试验研究.

农业工程学报,22, 74-77.]

[本文引用: 2]     

[20] Liu M, Westphal DL (2001).

A study of the sensitivity of simulated mineral dust production to model resolution.

Journal of Geophysical Research, 106, 18099-18112.

[本文引用: 1]     

[21] Lyles L, Tararko J (1986).

Wind erosion effects on soil texture and organic matter.

Journal of Soil and Water Conservation, 41, 191-193.

[22] Namikas S (2003).

Field measurement and numerical modelling of aeolian mass flux distributions on a sandy beach.

Sedimentology, 50, 303-326.

[23] Ni JR, Li ZS, Mendoza C (2003).

Vertical profiles of aeolian sand mass flux.

Geomorphology, 49, 205-218.

[本文引用: 1]     

[24] Pimentel D, Kounang N (1998).

Ecology of soil erosion in ecosystems.

Ecosystems, 1, 416-426.

[本文引用: 1]     

[25] Santini M, Caccamo G, Laurenti A, Noce S, Valentini R (2010).

A multi-component GIS framework for desertification risk assessment by an integrated index.

Applied Geography, 30, 394-415.

[本文引用: 1]     

[26] Shabani F, Kumar L, Esmaeili A (2014).

Improvement to the prediction of the USLE K factor.

Geomorphology, 204, 229-234.

[本文引用: 1]     

[27] Shao Y, Raupach M (1992).

The overshoot and equilibration of saltation.

Journal of Geophysical Research, 97, 20559-20564.

[本文引用: 1]     

[28] Shao Y, Raupach M, Findlater P (1993).

Effect of saltation bombardment on the entrainment of dust by wind.

Journal of Geophysical Research, 98, 12719-12726.

[本文引用: 2]     

[29] Shao Y, Raupach M, Leys J (1996).

A model for predicting aeolian sand drift and dust entrainment on scales from paddock to region.

Soil Research, 34, 309-342.

[本文引用: 1]     

[30] Su YZ, Zhao WZ (2005).

Soil organic carbon dynamics: Wind erosion effect.

Acta Ecologica Sinica, 25, 2049-2054.
(in Chinese with English abstract) [苏永中, 赵文智 (2005).

土壤有机碳动态: 风蚀效应. 生态学报,

25, 2049-2054.]

[本文引用: 1]     

[31] Sun YC, Ma SS, Chen Z, Zhao YL, Sun YR (2007).

Test and analysis of wind erosion of land surface soil in arid and semi-arid regions in north areas of Yinshan Mountain.

Transactions of the CSAE, 23, 1-5.(in Chinese with English abstract)
[孙悦超, 麻硕士, 陈智, 赵永来, 孙宇瑞 (2007).

阴山北麓干旱半干旱区地表土壤风蚀测试与分析.

农业工程学报,23, 1-5.]

[本文引用: 1]     

[32] Wagner LE (2013).

A history of wind erosion prediction models in the united states department of agriculture: The wind erosion prediction system (WEPS).

Aeolian Research, 10, 9-24.

[本文引用: 1]     

[33] Wang ZW (2009).

Effect of Stocking Rate on Ecosystem stability of Stipa brevilora Desert Steppe. PhD dissertation,

Inner Mongolia Agricultural University, Hohhot. 5-35.(in Chinese)
[王忠武 (2009).

载畜率对短花针茅荒漠草原生态系统稳定性的影响

. 博士学位论文, 内蒙古农业大学, 呼和浩特. 5-35.]

[本文引用: 2]     

[34] Webb NP, Okin GS, Brown S (2014).

The effect of roughness elements on wind erosion: The importance of surface shear stress distribution.

Journal of Geophysical Research, 119, 6066-6084.

[本文引用: 1]     

[35] Wiggs GFS, Livingstone I, Thomas DSG, Bullard JE (1994).

Effect of vegetation removal on airflow patterns and dune dynamics in the southwest Kalahari desert.

Land Degradation and Development, 5, 13-24.

[本文引用: 1]     

[36] Yang XH, He Q, Ali M. Huo W, Liu XC (2013).

Observational study on near-surface horizontal san-dust flux of sandstorms in the southeastern fringe of the Taklimakan desert.

Journal of Desert Research, 33, 1299-1304.(in Chinese with English abstract)
[杨兴华, 何清, 艾力·买买提依明, 霍文, 刘新春 (2013).

塔克拉玛干沙漠东南缘沙尘暴过程中近地表沙尘水平通量观测研究. 中国沙漠,

33, 1299-1304.]

[本文引用: 1]     

[37] Yu GM, Liu Y, Yan Y, Hu YF (2011).

Soil wind erosion risk assessment in the middle part of Inner Mongolia Plateau during 2000 to 2008.

Scientia Geographica Sinica, 31, 1493-1499.
(in Chinese with English abstract) [于国茂, 刘越, 艳燕, 胡云锋 (2011).

2000-2008年内蒙古中部地区土壤风蚀危险度评价.

地理科学,31, 1493-1499.]

[38] Zhang XL, Zhou QQ, Chen WW, Wang YY, Tong DQ (2015).

Observation and modeling of black soil wind-blown erosion from cropland in northeastern China.

Aeolian Research, 19, 153-162.

[本文引用: 1]     

[39] Zhao CX, Zheng DW, He WQ (2005).

Vegetation cover changes over time and its effects on resistance to wind erosion.

Acta Phytoecologica Sinica, 29, 68-73.(in Chinese with English abstract)
[赵彩霞, 郑大玮, 何文清 (2005).

植被覆盖度的时间变化及其防风蚀效应.

植物生态学报,29, 68-73.]

[本文引用: 1]     

[40] Zobeck TM, Stout JE (1996).

The wolfforth field experiment: A wind erosion study.

Soil Science, 161, 616-632.

[41] Zou XY, Zhang CL, Cheng H, Kang LQ, Wu XX, Chang CP, Wang ZL, Zhang F, Li JF, Liu CC, Liu B, Tian JL (2014).

Classification and representation of factors affecting soil wind erosion in a model.

Advances in Earth Science, 29, 875-889.(in Chinese with English abstract)
[邹学勇, 张春来, 程宏, 亢力强, 吴晓旭, 常春平, 王周龙, 张峰, 李继峰, 刘辰琛, 刘博, 田金鹭 (2014).

土壤风蚀模型中的影响因子分类与表达.

地球科学进展,29, 875-889.]

[本文引用: 2]     

/