基于RWEQ模型的典型沙区风蚀治理成效研究

张寒冰1,2 高阳1,2,† 崔艳智3

1.中国农业大学土地科学与技术学院, 北京 100193; 2.国土资源部农用地质量与监控重点实验室, 北京 100193; 3.中日友好环境保护中心, 北京 100029; †通信作者, E-mail: yanggao@cau.edu.cn

摘要 以通辽市科尔沁左翼后旗(科左后旗)为研究区, 通过修正的风蚀方程(revised wind erosion equation, RWEQ)计算科左后旗2010—2016年的风蚀强度时空变化, 评价通辽市开展科尔沁沙地“双千万亩”综合治理工程的风沙治理效果。研究结果表明: 1) 2010—2016年, 科左后旗风蚀强度呈现明显下降的趋势, 总风蚀量从2010年的3893.87万t下降到2016年的1883.56万t, 降幅为51.63%; 2) 工程实施区内平均风蚀强度高于全旗平均风蚀强度, 且相对风蚀系数在2016年明显下降, 说明项目建设成效显著, 有效地削弱了风蚀较高的项目区内的侵蚀; 3) 项目区建设具有一定的域外效应, 范围在3km左右, 项目实施后, 不同距离缓冲区的相对风蚀系数均明显下降, 且越靠近项目区下降越明显, 从而减小了缓冲区之间的差距。

关键词 风蚀; 修正的风蚀方程(RWEQ); “双千万亩”综合治理工程; 科尔沁左翼后旗

土壤风蚀是土壤表层在风力作用下, 相对较细小的颗粒被吹蚀、搬运与沉积的过程[1], 是半干旱地区土地退化的主要原因[2–3]。作为内蒙古自治区面积最大的沙地, 科尔沁沙地是我国北方七大强风蚀区之一, 不合理放牧和种植导致土壤风蚀加剧, 扬尘起沙严重, 地表土壤平均粒径增大[4–5], 又因其特殊的地理位置而成为华北平原沙尘的主要来源, 严重影响当地以及华北地区居民的生产生活[6]

为缓解风蚀造成的土地退化和环境污染问题[7], 2000年6月, 国家在科尔沁沙地启动京津风沙源治理一期工程, 2012年启动二期工程, 通过荒山造林、封山育林和退耕还草等一系列措施增加科尔沁沙地的植被覆盖, 治沙止漠。经过十多年的建设, 工程区风蚀模数整体上呈下降趋势, 相对容易治理的沙地已经得到初步改善, 林地与牧草地面积增大, 风蚀程度有所减轻[8]。但是, 目前科尔沁沙地仍有部分“远沙大沙”没有得到有效的治理, 区域环境恶化状况没有从根本上得到扭转, 沙漠化形势仍然严峻[9]。为进一步加快防沙治沙步伐, 改善沙区群众生产生活环境, 2014年, 位于科尔沁沙地腹地的通辽市政府开展“双千万亩”综合治理工程, 计划利用7年时间完成1000万亩林业生态治理工程和1000万亩草原生态治理工程, 逐步恢复疏林草原生态系统, 以期减缓土壤风蚀和土地退化, 提高当地生态环境质量。

修正的风蚀方程(revised wind erosion equation, RWEQ)是一种以较高时空分辨率对区域土壤风蚀状况进行长时间序列估算, 从而有效预测风蚀量的模型[10], 可以为土地沙化防治提供依据[11]。江凌等[12]利用RWEQ模型计算青海省的风蚀情况, 发现在国内区域尺度上运用RWEQ模型具有一定的可行性。巩国丽等[13]将RWEQ模型应用于内蒙古锡林郭勒盟, 计算1991—2010年间20年的土壤风蚀空间分布, 发现位于锡林郭勒盟西部的浑善达克沙地地区土壤侵蚀最严重。郭中领提出RWEQ模型在我国北方的参数修订方案[14], 并分析环境治理政策对区域风蚀强度的影响[7]。但是, 现有研究多侧重于探究风蚀强度的时空分布, 与区域政策管理, 尤其是工程治理成效评估相结合的研究较少[15]。此外, 现有研究在小时风速估计方面多采用模型内插模拟的方法, 在一定程度上低估了局地风蚀量[16]。在局地起沙风速估计中, 多采用系统默认值, 对当地实际情况考虑不足[17]。科学客观地评价风沙防治工程建设成效, 度量区域风蚀强度的时空变化特征, 评估工程建设对风蚀的削减作用, 可为区域风沙治理与管控提供科学依据。

本研究基于RWEQ模型, 以通辽市科尔沁左翼后旗(简称科左后旗)为例, 使用时间分辨率较高的实测小时风速参与RWEQ模型的计算, 并通过实地测量科左后旗剖面土壤含水量来修正当地起沙风速, 提高风因子计算精度, 最终得到2010, 2012, 2014和2016年4期风蚀强度的时空变化, 测度重点项目区内外风蚀变化的差异, 评估风沙防治工程成效, 为典型沙区生态系统保护和防风固沙工作提供科学依据。

1 研究区概况

通辽市科左后旗(121.5°—123.7°E, 42.7°—43.7°N)地处科尔沁沙地东南端(图1), 年降水量为350~ 450mm, 植被类型自东向西由森林草原向典型草原、荒漠草原过渡, 土地利用以耕地、草地和林地为主, 耕地集中在水利灌溉条件较好的东部及东南部地区, 草地广泛分布在西部及北部地区。旗内多风沙, 风积物覆盖范围广, 流动和半流动沙丘占全旗面积的47.2%, 沙地土壤表层颗粒结构松散, 有一定的粗化现象, 可风蚀性较强[18], 风沙干旱、水土流失、土地沙化和草场退化是该旗的主要特点, 生态环境十分脆弱, 土壤风蚀危害大。

作为典型的半干旱沙区[19]及强风蚀区[20], 科左后旗是《内蒙古自治区林业发展“十三五”规划》中东部风蚀沙化治理的重点区域, 也是“京津风沙源治理工程(一期及二期)”、“天然林保护”、“退耕还林(草)”等一系列防沙治沙工程的重点实施区域。2014年, 科左后旗开始实施科尔沁沙地“双千万亩”综合治理工程, 依照“因地制宜、适地适树”的原则, 对沙化严重地区进行林业生态综合治理, 并根据立地条件实施退耕还林还草, 逐步引导农牧民实现全部退耕, 恢复天然林原貌。至2020年, 科左后旗拟通过强化治沙造林措施, 使未治理沙化土地得到初步保护与治理, 森林覆盖率年均提高1.2%以上, 草原植被盖度年均提高2%以上[21]。重点项目区主要分布在努古斯台、阿古拉及好通高速道路沿线等地(表1)。

2 研究方法与数据来源

2.1 RWEQ模型

width=423.7,height=300.95

图1 科左后旗位置及“双千万亩”重点项目区分布

Fig. 1 Location of Horqin Left Back Banner and the key project area of “Double Ten Million Mu”

表1 科左后旗“双千万亩”重点项目区范围及实施进展

Table 1 Scope and implementation progress of the key project area of “Double Ten Million Mu”

项目区治理措施范围总实施目标截止2016年底实施情况 好通高速好通高速科左境内道路两侧200m缓冲区内栽植樟子松、五角枫等; 1000m以内区域划定为全年禁牧区域好通高速公路沿线, 涉及努古斯台、巴嘎塔拉、甘旗卡3个苏木镇规划治理面积4.5万亩, 涉及公路段长89km完成3万亩 阿古拉项目区立地条件较好的地块栽植杨树、樟子松、五角枫等; 流动沙丘、沙化严重的起伏较大的固定沙丘栽植锦鸡儿、紫穗槐、五角枫等; 全年禁牧阿古拉镇色音呼都嘎、吉力吐、玛琳秋、白兴吐、东花灯、舍根、乌旦塔拉7个嘎查规划治理面积25万亩完成25万亩 努古斯台项目区栽种樟子松、油松、五角枫、榆树、蒙桑、锦鸡儿、紫穗槐、金叶榆、红刺榆等; 全年禁牧努古斯台镇努古斯台、套海、花卜、图古日格、道乃营子、瑙木塔拉6个嘎查规划治理面积50万亩完成30万亩

说明: 资料来源于科左后旗林业局。

修正的风蚀方程(Revised Wind Erosion Equa-tion, RWEQ)计算多因素影响下的区域土壤转运总量(SL), 作为该区域的土壤风蚀量[11]。首先计算受气象、土壤侵蚀、土壤结皮、地表粗糙度和植被覆盖状况等因素影响的风力最大转运容量Qmax(式(1)), 当土壤转移量达到Qmax的63.2%时, 对应的地块长度称为关键地块长度S(式(2)), 然后通过QmaxS计算得到区域土壤转运总量SL (式(3))。

Qmax=109.8×(WF×EF×SCF×K'×C), (1)

S=150.71×(WF×EF×SCF×K'×C)−0.3711, (2)

width=81.7,height=32.25 (3)

式中: Qmax为风力最大运转量(kg/m); S为关键地块长度(m); SL 为土壤损失量(kg/m2); x 表示距上风向距离(m), 取50m; WF表示气象因子(kg/m); EF表示土壤可侵蚀程度(无量纲); SCF 表示土壤结皮因子(无量纲); 表示土壤粗糙度因子(无量纲); C 表示植被覆盖因子(无量纲)。

2.1.1 气象因子(WF)

风是风蚀发生的基本动力。气象因子的计算公式为

width=163.85,height=44.6, (4)

式中, U2为2m高处风速(m/s), Ut为2m处的起沙临界风速(m/s)。需要注意的是, 起沙风速一般随沙中含水量的增加而增加, 在含水量为0.8%左右时, 起沙风速约为4.5~5.0m/s[22]。科尔沁沙地的沙土以中细沙为主, 本研究组于2016年6—9月在科左后旗乌旦塔拉村布控12个土壤剖面, 采用TDR100时域反射仪测度土壤水分含量, 发现土壤含水量均在1%以下, 此处起沙风速取5m/s。Nd 表示一次实验的观测天数(d), N为一次实验中观察总次数; ρ为空气密度(kg/m3); g表示重力加速度(m/s2), 本文取9.8m/s2; SW为土壤湿度因子; SD为雪盖因子。其中ρ由海拔高度EL(km)和绝对温度T(K)计算得到:

width=166.5,height=32.25, (6)

width=99.4,height=44.6, (7)

width=148.85,height=29.55, (8)

ETP为潜在相对蒸发量(mm)[23], R为降雨量(mm), I为灌溉量(mm), 取I=0; Rd为降雨或灌溉的次数或天数(d); Nd为观测天数(d); SR为太阳辐射总量(cal/cm2); DT为平均温度(ºC)。

SD =1- P(雪深>25.4mm), (9)

P为此次计算时段内雪盖深度大于25.4mm的概率。

2.1.2 土壤可蚀因子(EF)与土壤结皮因子(SCF)

EF值和SCF值与土壤的基本物理化学性质相关[2425]。本研究采用如下公式计算:

width=121.5,height=27.4width=10.2,height=10.2

width=156.35,height=27.4, (10)

width=156.4,height=29.55, (11)

式中: Sa为土壤砂粒含量(5.50%~93.60%); Si为土壤粉砂含量(0.50%~69.50%); Cl为黏土含量(5.00%~ 39.30%); Sa/Cl为土壤砂粒和黏土含量比(1.20%~ 53.00%); OM表示有机质含量(0.32%~4.74%); CaCO3表示碳酸钙含量(0~25.20%), 括号中百分数为RWEQ模型要求的数据含量范围, 本研究区土壤属性绝大部符合模型要求, 对少数不符合模型准入范围的土壤参数, 采用RWEQ内嵌的土壤质地资料建议的输入参数。

2.1.3 地表粗糙度因子(K')

地表糙度分为土垄糙度Kr 和随机糙度Crr, 地表糙度因子K'的计算公式如下:

width=106.95,height=17.2。 (12)

由于本研究区土地随机粗糙度的空间数据难以获得, 故 Crr 值取 0。地形粗糙长度一般通过 Smith-Carson方程计算:

width=74.7,height=27.4, (13)

式中, Δ表示距离 L 范围内的海拔高程差(cm), L为地势起伏参数(cm), 本研究中 L 取5×105cm[26]

2.1.4 植被覆盖因子(C)

植被覆盖度是影响风蚀的关键因子, 其覆盖程度直接影响近地表风速以及土地粗糙程度, 计算公式如下:

C=e−0.0438VEG, (13)

width=124.15,height=29.55, (14)

式中: VEG表示植被覆盖度(%); NDVIii栅格中的NDVI值; NDVIsoil为研究区裸地对应NDVI值, 取最小5%对应值; NDVImax为研究区年NDVI最大值。

2.2 相对风蚀系数

气候和植被覆盖是影响风蚀的主要原因。科左后旗项目区建设主要通过增加植被覆被的手段改善风蚀状况, 并不能改变或影响气候等不可控因子。不可控的天气状况导致科左后旗各年份之间风蚀总量不同。为更好地表征项目内外风蚀程度的相对变化, 本文设置相对风蚀系数α来表示项目区相对于同年份全旗风蚀的程度, 避免由于不同年份气候差异带来的不确定性, 计算公式如下:

width=29.55,height=29.55, (15)

其中, Qi表示该年项目区平均风蚀量, Qo表示该年全旗平均风蚀量。α 值越大, 相对风蚀程度越高; 当α>1时, 表明项目区平均风蚀程度高于全旗平均水平; 当α=1时, 表明项目区平均风蚀程度与全旗平均水平一致。

2.3 数据来源

气象因子中的风因子和土壤湿度因子利用中国气象科学数据共享服务网(http://cdc.cma.gov.cn)提供的科左后旗及周边14个气象台站的风速、降水和温度等数据, 经克里金差值计算得到, 其中风速采用小时风速(每小时前后的10分钟平均风速), 每天24期。此外, 本研究组于2015年6—10月和2016年6—9月赴科左后旗进行野外调查, 并在甘旗卡镇乌旦塔拉村进行野外土壤剖面观测, 获得土壤含水量等数据。日照时数来源于中国气象科学数据共享服务网的中国辐射日值数据集。雪盖因子采用中国西部环境与生态科学数据中心(http://westdc.westgis.ac. cn)的中国雪深长时间序列数据集计算得到。DEM数据来源于ASTER数据生成的全球数字高程模型Global Digital Elevation Model(https://wist.echo.nasa. gov/~wist/api/imswelcome/), 空间分辨率为30m。土壤砂土、粉砂土与黏土、CaCO3和有机碳含量数据来自世界土壤数据库(Harmonized World Soil Database, HWSD), 并结合实地测量数据矫正得到。植被覆盖度采用MODIS MOD13Q1V006每16天一期最大合成的NDVI数据计算得到。所有数据经重采样, 空间分辨率为250m, 时间分辨率为16天, 全年共23期, 计算结果累加得到年总风蚀量。

3 结果与分析

3.1 科左后旗2010—2016年风蚀强度的时空变化

width=442.05,height=340.9

图2 2010—2016年科左后旗风蚀强度的空间分布

Fig. 2 Spatial pattern of wind erosion of Horqin Left Back Banner from 2010 to 2016

依照水利部《土壤侵蚀分类分级标准(SL190—2007)》分类标准, 将科左后旗境内风蚀强度划分为微度侵蚀(0~2t/(hm2·a))、轻度侵蚀(2~25t/ (hm2·a))、中度侵蚀(25~50t/(hm2·a))、强烈侵蚀(50 ~80t/(hm2·a))和极强烈侵蚀(80~150t/(hm2·a))5 个等级。计算结果显示, 科左后旗风蚀强度以轻度侵蚀和中度侵蚀为主(图2), 4年平均风蚀量为19.49t/(hm2·a), 2010, 2012和2014年均呈现明显的西高东低的空间分布规律, 但2016年全旗风蚀量无明显空间趋势。在年内变化上, 取科左后旗4年各期风蚀量的平均值, 并按月统计(图3), 结果表明全旗年内风蚀量较大的时间段集中于春季, 3—5月风蚀尤为强烈, 占全年风蚀量的89.16%。

在时间变化上, 2010—2016年科左后旗风蚀强度明显下降(图4), 总风蚀量从2010年的3893.87万t下降为2016年的1883.56万t, 降幅为51.63%, 平均风蚀强度从2010年的33.86t/(hm2·a)下降到2016年的16.37t/(hm2·a), 年均下降3.49t/hm2, 全旗风蚀状况得到明显改善。此外, 2010—2016年不同风蚀程度面积占比也明显改变, 微度和轻度风蚀面积占比呈上升趋势, 强烈风蚀和极强烈风蚀面积占比明显下降(图5)。其中, 强烈风蚀地区面积占比从2010年的24.02%下降到2016年的8.08%; 微度风蚀从2.21%增加至9.59%, 轻度风蚀从41.71%增加至65.48%, 增幅较大。从总体上看, 自2010年以来, 作为京津风沙源治理工程项目区之一, 科左后旗风蚀量下降明显, 强烈风蚀与极强烈风蚀侵蚀地区的面积显著减少, 风蚀状况得到一定程度的缓解。

width=206.85,height=119

图3 科左后旗年内各月风蚀量分布(2010—2016)

Fig. 3 Monthly wind erosion distribution of Horqin Left Back Banner (2010—2016)

width=185.5,height=117.6

图4 科尔沁左翼后旗2010—2016年年均风蚀趋势

Fig. 4 Trend of annual average wind erosion in the Horqin Left Back Banner from 2010 to 2016

3.2 科左后旗项目区2010—2016年风蚀强度的时空变化

依照通辽市“双千万亩”项目分区治理方案, 科左后旗工程选取努古斯台、阿古拉及好通高速道路沿线等明沙集聚、沙化最严重的村镇作为重点治理区域。根据以上测算结果, 科左后旗项目区2010—2016年平均风蚀量为24.51t/(hm2·a), 明显高于科左后全旗的平均值(19.48t/hm2·a), 说明区域选择合理; 从时间上看, 如表2所示, 2010年项目区平均风蚀量为43.03t/(hm2·a), 2016年降至18.11t/(hm2·a), 下降幅度达到57.91%。

此外, 2010—2016年科左后旗项目区相对风蚀系数 α 分别为1.27, 1.37, 1.19和1.10, 从图 6 的项目区 α 变化曲线(实线)明显看出, 2010—2014年, “双千万亩工程”尚未开展时, 项目区相对风蚀系数较高, 即项目区平均风蚀状况与全旗平均状况的差距更明显; 2016年项目区相对风蚀系数大幅度下降, 说明区域内平均风蚀状况与全旗平均状况的差距减小, 仅为全旗平均风蚀量的1.10倍。由此可见, 自2014年“双千万亩”工程实施后, 与全旗其他区域相比, 原本风蚀强度较大的项目区风蚀状况出现明显改善, 到2016年, 项目区风蚀强度与全旗平均水平基本上一致, 说明2014—2016年间“双千万亩”工程建设效果较显著。

width=331.65,height=127.65

图5 科左后旗2010—2016年各类风蚀程度面积占比

Fig. 5 The proportion of different levels of wind erosion in the Horqin Left Back Banner from 2010 to 2016

表2 2010—2016年各区域年均风蚀

Table 2 Annual average wind erosion in each region from 2010 to 2016

年份年均风蚀/(t·hm−2·a−1) 项目区内1km缓冲区2km缓冲区3km缓冲区4km缓冲区5km缓冲区全旗 201043.0339.9737.0134.8333.2133.1533.8581 201229.3726.3924.0622.3821.1821.1421.3834 20147.557.196.576.225.905.946.3291 201618.1116.7215.9215.5315.2515.6816.378

3.3 科左后旗项目区周边地区风蚀强度的时空变化

科左后旗“双千万亩”工程建设以防风林栽种及禁牧等措施为主[27]。栽种防风林可在覆盖土壤表面形成保护层, 减缓近地表风速, 阻碍风湿过程中砂砾的运输过程, 进而减缓风蚀[28]。此外, 防风林还可通过吸引传播种子的动物来增加周围地区树木和灌木种子的沉积, 进一步通过生物方式加速自然的演变过程, 提升周边地域植被覆盖度[29]。可见, 防风林建设不仅可以减缓项目区内风蚀状况, 也可减缓周边地区风蚀状况, 即具有一定的域外效应。因此, 为进一步探究项目区周边地区的风蚀变化, 我们在项目区外1km处、1~2km处、2~3km处、3~4km处和4~5km处划分缓冲区(以下称为1km缓冲区、2km缓冲区、3km缓冲区、4km缓冲区和5 km缓冲区), 分别统计2010—2016年各缓冲区的风蚀变化量(表2), 并计算相对风蚀系数(图6)。

在时间变化上, 2014年项目启动后, 科左后旗项目区外不同距离缓冲区的相对风蚀系数均下降(图6), 说明项目区建设降低了项目区域内及域外一定距离缓冲区的风沙侵蚀。2010—2016年各区域相对风蚀系数的标准差逐渐下降, 分别为0.12, 0.15, 0.11和0.06, 不同区域间风蚀程度差距减小, 进一步说明项目区治理取得一定成效。

width=205.65,height=195.1

图6 不同缓冲区2010—2016年相对风蚀系数变化

Fig. 6 Changes of relative wind erosion coefficient in different buffer regions from 2010 to 2016

在空间分布上, 项目区的风蚀程度高于全旗的平均水平, 距项目区距离越远, 风蚀程度越低(表2)。3km缓冲区2010—2016年的相对风蚀系数分别为1.03, 1.04, 0.98和0.97(图6), 即项目区外2~3km处风蚀强度已与全旗平均水平相当。此外, 2016年项目区内及1~5km缓冲区的相对风蚀系数分别下降0.16, 0.15, 0.012, 0.08, 0.05和0.02, 说明越靠近项目区, 相对风蚀系数降低越明显。但是, 在3km缓冲区外, 相对风蚀系数的变化不明显, 4km缓冲区与5km缓冲区相对风蚀系数相差不大, 与全旗平均水平相近。2014—2016年 4km缓冲区相对风蚀系数几乎没有变化, 5km缓冲区的相对风蚀系数甚至出现一定程度的增加, 因此项目区建设的域外效应边界约在项目区外3km处。

4 讨论

与文献[30]相比, 本研究模拟得到的区域风蚀强度略高, 但时空分布格局基本上一致。在利用传统方法的研究中, 多通过统计型风速生成器(wind Generator)模拟风速, 以期达到模型要求, 这可能会在一定程度上低估风速[16]。本研究采用国家气象局信息中心从2005年开始收集的小时风速作为因子计算的数据来源, 风速的时间分辨率更高。因此, 模拟得到的风蚀量较大。

此外, 以风速为代表的气候因素是影响风蚀的主要原因, 风沙治理主要通过增加植被覆被的手段来改善区域风蚀状况, 不能改变或影响气候等不可控因子。本研究模拟结果中, 科左后旗2014年风蚀量较小, 2016年风蚀量略有增加, 但仍低于项目启动前2010年及2012年的两期结果。统计研究区周边各气象站2010—2016年的大风天气(2m处风速大于8m/s)的频次, 发现2010—2016年大风次数分别为289, 182, 92和192次, 2014年大风发生的次数明显少于其他三期, 这也是导致2014年风蚀量小的重要因素。

本研究存在一些不足之处。1)RWEQ风蚀模型中, 计算风蚀因子所需的风速为2m处风速, 但我国气象站数据为10m处风速, 计算中使用的2m处风速是用《中华人民共和国气象行业标准QX/T 81—2007》提供的校正方式校正得到, 因缺乏实地风速的检测与校正, 与真实值相比, 可能存在一些偏差。2)本文的主要目的是探讨科左后旗防风固沙项目的建设成效, 故研究的时间范围(2010, 2012, 2014和2016年)是根据项目建设时间确定的。由于气候具有一定的随机性, 在接下来的研究中, 应采用较长的时间序列探究科左后旗以及京津风沙源治理工程重点实施区域的风蚀强度变化。3)本文提出防风林的建设具有一定的域外效应, 但由于本研究区域仅限科左后旗, 未考虑周边地区的防风固沙建设工作对科左后旗风蚀强度的影响, 因此, 扩大研究区范围是今后工作的重点之一。

5 结论

“双千万亩”整治工程是通辽市改善当地环境、治理土壤风蚀的重要尝试。本文以通辽市科左后旗为研究区, 通过RWEQ模型计算区域内风蚀强度, 评估其实施效果, 主要结论如下。

1)2010—2016年科左后旗风蚀强度明显下降, 总风蚀量从2010年的3893.87万t下降到2016年的1883.56万t, 降幅为51.63%, 防风固沙整治工程项目在全旗取得一定的成效。在空间分布上, 科左后旗风蚀强度自东向西递增, 风蚀较强的时间主要集中在3—5月, 因此, 下一步应加强春季的风沙治理工作, 并将资金和力量集中在西部地区。

2) “双千万亩”项目工程实施前, 项目区平均风蚀强度明显大于全旗平均风蚀; 2014年项目区开始建设之后, 项目区内相对风蚀系数明显下降。项目区2014—2016的建设成效显著, 有效地缓解了项目区内的风沙侵蚀。

3)项目区的建设具有一定的域外效应。距项目区不同距离缓冲区内, “双千万亩”项目实施后, 相对风蚀系数都明显下降, 且越靠近项目区, 下降越明显, 且各区域之间的差距减小。距离风蚀区2~3km范围内的风蚀强度与全旗风蚀强度基本上一致, 相对风蚀系数下降不显著, 项目区的域外效应范围约为3km。

参考文献

[1]丁国栋. 风沙物理学. 北京: 中国林业出版社, 2010

[2]Yan Y, Xin X, Xu X, et al. Quantitative effects of wind erosion on the soil texture and soil nutrients under different vegetation coverage in a semiarid steppe of northern China. Plant & Soil, 2013, 369 (1/2): 585–598

[3]陈渭南, 董光荣. 中国北方土壤风蚀问题研究的进展与趋势. 地球科学进展, 1994, 9(5): 6–12

[4]常学礼, 鲁春霞, 高玉葆. 人类经济活动对科尔沁沙地风沙环境的影响. 资源科学, 2003, 25(5): 78–83

[5]张正偲, 董治宝. 土壤风蚀对表层土壤粒度特征的影响. 干旱区资源与环境, 2012, 26(12): 86–89

[6]何文清, 赵彩霞, 高旺盛, 等. 不同土地利用方式下土壤风蚀主要影响因子研究——以内蒙古武川县为例. 应用生态学报, 2005, 16(11): 88–92

[7]Guo Z, Huang N, Dong Z, et al. Wind erosion induced soil degradation in Northern China: status, measures and perspective. Sustainability, 2014, 6(12): 8951–8966

[8]于宝勒, 吴文俊, 赵学军, 等. 内蒙古京津风沙源治理工程土壤风蚀控制效益研究. 干旱区研究, 2016, 33(6): 1278–1286

[9]韩已文, 刘忠, 于学谦, 等. 沙化土地整治背景下的沙漠化监测预警. 农业工程学报, 2017, 33(10): 271–277

[10]Fryrcar D W, Chen W N, Lester C. Revised Wind Erosion Equation. Annals of Arid Zone, 2001, 40(3): 265–279

[11]Fryrear D W, Bilbro J D, Saleh A, et al. RWEQ: improved wind erosion technology. Journal of Soil & Water Conservation, 2000, 55(2): 183–189

[12]江凌, 肖燚, 欧阳志云, 等. 基于RWEQ模型的青海省土壤风蚀模数估算. 水土保持研究, 2015, 22 (1): 21–32

[13]巩国丽, 刘纪远, 邵全琴. 基于RWEQ的20世纪90年代以来内蒙古锡林郭勒盟土壤风蚀研究. 地理科学进展, 2014, 33(6): 825–834

[14]郭中领. RWEQ模型参数修订及其在中国北方应用研究[D]. 北京: 北京师范大学, 2012

[15]南岭, 杜灵通, 展秀丽. 土壤风蚀可蚀性研究进展. 土壤, 2014, 46(2): 204–211

[16]Van Donk S J, Wagner L E,Skidmore E L, et al. Comparison of the Weibull model with measured wind speed distributions for stochastic wind genera-tion. 2005, 48(2): 503–510

[17]朱好, 张宏升. 沙尘天气过程临界起沙因子的研究进展. 地球科学进展, 2011, 26(1): 30–38

[18]马玉凤, 介冬梅, 严平, 等. 科尔沁东部沙地土壤可风蚀性研究. 世界地质, 2007, 26(3): 338–344

[19]杨振常, 李玉霖, 崔夺, 等. 半干旱典型沙区1951 —2005年主要气象要素及潜在蒸散量的变化趋势研究. 中国沙漠, 2012, 32(5): 1384–1392

[20]韩小红, 杨联安, 周欢水, 等. 我国北方强风蚀区空间分布格局及特征分析. 干旱区资源与环境, 2008, 22(3): 113–116

[21]科左后旗林业局. 科尔沁沙地“双千万亩”综合治理工程努古斯台项目区简介[EB/OL]. (2015–02–28) [2017–11–18]. http://124.205.185.15/xlyj/subsite/ 33286/33301/181961.html

[22]殷海军. RWEQ模型的改进及其在风沙区的应用研究[D]. 北京: 中国农业大学, 2007

[23]Samani Z A, Pessarakli M. Estimating potential crop evapotranspiration with minimum data in Arizona. American Society of Agricultural Engineers, 1986, 29 (2): 522–524

[24]Fryrear D W, Bilbro J D. Wind erosion control with residues and related practices // Unger P W. Managing agricultural residues. Boca Raton, FL, 1994: 7–18

[25]Li Y, Hagen G, Guilfoyle T J. Altered morphology in transgenic tobacco plants that overproduce cytokinins in specific tissues and organs. Developmental Bio-logy, 1992, 153(2): 386–395

[26]李军, 游松财, 黄敬峰. 基于GIS的中国陆地表面粗糙度长度的空间分布. 上海交通大学学报(农业科学版), 2006, 24(2): 185–189

[27]内蒙古自治区人民政府. 内蒙古自治区林业发展“十三五”规划[EB/OL]. (2017–04–02)[2017–11–20]. http://www.baotou.gov.cn/info/2057/127635.htm

[28]van de Ven T A M, Fryrear D W, Spaan W P. Vegeta-tion characteristics and soil loss by wind. Journal of Soil & Water Conservation, 1989, 44(4): 347–349

[29]Harvey C A. Windbreaks enhance seed dispersal into agricultural landscapes in Monteverde, Costa Rica. Ecological Applications, 2000, 10(1): 155–173

[30]Ling J, Yi X, Hua Z, et al. Spatio-temporal variation of wind erosion in Inner Mongolia of China between 2001 and 2010. Chinese Geographical Science, 2016, 26(2): 155–164

Evaluation the Effects of Wind Erosion Control Projects in Typical Sandy Area Based on RWEQ Model

ZHANG Hanbing1,2, GAO Yang1,2,†, CUI Yanzhi1

1. College of Land Science and Technology, Beijing 100193; 2. Key Laboratory of Agricultural Land Quality, Ministry of Land and Resources, Beijing 100193; 3. Sino-Japan Friendship Centre for Environmental Protection, Beijing 100029; † Corresponding author, E-mail: yanggao@cau.edu.cn

Abstract The revised wind erosion equation (RWEQ) is used to calculate the spatial and temporal changes of the wind erosion in Horqin Left Back Banner from 2010 to 2016, to evaluate the implementation effect of “Double Ten Million Mu” comprehensive management project in Horqin sandy land carried out by Tongliao Government. The results show that the wind erosion decreases significantly in Horqin Left Back Banner from 2010 to 2016, and the total wind erosion decreases from 38.93 million tons in 2010 to 18.83 million tons in 2016, down by 51.63%. The average wind erosion in the project implementation area is higher than the average wind erosion in Horqin Left Back Banner, which shows that the project construction achieves remarkable achievements and reduces the erosion effectively in the project area with high wind erosion. Project area construction has certain extraterritorial effect, and the positive effect of the project construction works not only within the project area, but also in the surrounding areas. The effective range is about 3km far from the project area. After two years of project implementation, the relative wind erosion coefficient of different distance buffers decreases, and the closer to the project area the buffer is, the more positive effects it behaves.

Key words wind erosion; RWEQ; “Double Ten Million Mu” comprehensive management project; Horqin Left Back Banner

doi: 10.13209/j.0479-8023.2019.010

收稿日期: 2018-03-26;

修回日期: 2018-04-17;

网络出版日期: 2019-02-28

国家自然科学基金(41501087)资助