2017年8月8日21时19分46秒在四川省北部阿坝州九寨沟县发生7.0级地震,震中位于九寨沟核心景区西部5km附近。强震往往会导致当地的土壤侵蚀量增加,造成当地水土流失和环境破坏。例如,“5·12”汶川地震发生后,在地震影响下,受灾区的水土保持功能在不同程度上受损:四川理县震后侵蚀量相比2001年增加了12万t/a左右。茶坪山区水土严重受损,震后侵蚀量相比震前1年增加了72.3%。龙门山区相比震前1年侵蚀面积增加了86.3km2。北川县震后土壤侵蚀面积相比2000年全国土壤侵蚀遥感的调查成果增加275.13km2。1999年台湾集集发生的M-W7.6级大地震,在地震初期年均侵蚀深度达到6.04mm,是震前的2.26倍。
由此可见,地震会导致当地的水土保持功能受损,在侵蚀严重的区域,人民的生产生活和环境会受到严重的威胁;因此,在灾区不同区域,特别是像九寨沟这样以旅游业为主的区域,对土壤侵蚀进行定量评价,划分出受地震影响严重的区域,并分析该区域的空间分布规律,能够有助于现在或将来将有限的资源最有效的投入到灾后水土保持工作中,这对灾后生态恢复具有现实的意义。
国内外已有大量关于土壤侵蚀的研究。目前土壤侵蚀的模型主要分2种:1)以观测数据建立的经验模型;2)基于侵蚀机理建立的物理模型。经验模型出现的时间比较早,且模型结构也比物理模型精简,计算更加方便。
在经验模型方面,美国于1985年提出了通用土壤流失方程(USLE),被美国农业部以农业手册的方式颁布执行。USLE模型为土壤侵蚀研究提供了一个简单实用的工具,在全世界范围内得到广泛的应用。早期该模型在降雨条件下,可以预测年均土壤侵蚀量,可以计算出由于细沟侵蚀、面蚀造成的水土流失。
随着空间信息技术的发展,运用USLE模型对区域尺度上的土壤侵蚀进行定量评价成为现实。如Batjes利用0.5经度乘0.5纬度网格全球数据,通过USLE模型实现全球范围内土壤侵蚀的定量评价。我国学者通过修改USLE模型参数因子的计算方法,得到适合于我国土壤状况的土壤侵蚀模型。刘宝元等在2001基于USLE模型提出了中国土壤流失方程(CSLE)。卜兆宏等在2003年改进的侵蚀模型与USLE形式一样,只是修正了各个侵蚀因子的算法。
作者结合前人研究,在空间信息技术(GIS)支持下应用USLE模型定量计算出九寨沟地震前后受灾区的土壤侵蚀强度,并对地震前后土壤侵蚀参数作对比分析,划分出受地震影响较大的区域,探讨该区域土壤侵蚀的空间分布规律与距发震断层距离、坡度、高程和坡向4方面之间的关系。
研究结果表明:
1)
从侵蚀参数变化来看,相比震前,震后侵蚀面积增加8.1%,侵蚀总量提高28%,侵蚀强度增强18.39%。在侵蚀强度类别中轻微度、轻度和中度的面积占比在震后全部下降,强度侵蚀及以上侵蚀的面积占比都有所上升,因此震后水土保持工作应重视对中度侵蚀及以上的侵蚀区域的治理。
2)
对震后易受地震影响区域的空间分布规律和距断层距离进行分析。结果表明距发震断层距离越近,易受地震影响区域的面积越大,在距发震断层2km范围内,震后易影响区的面积最大,随着距断层距离逐渐增大,易影响区的面积也逐渐减小,在距发震断层6km以外,几乎没有易影响区。且发现易影响区的侵蚀强度相比震前同区域的侵蚀强度也呈现一定的规律,在易影响区范围内,随着距断层距离逐渐增大,易受影响区的侵蚀强度增量急剧减小,因此应该加强对距发震断层2km区域范围的水土保持工作。
3)
研究区易受地震影响的区域主要发生在高程2000~4000m,
占到总面积的99%,在30°~60°坡度带内,侵蚀面积分别占到易影响区总面积的89.25%,易受地震影响的区域在北坡、东北坡和东南坡的侵蚀面积达到61.07%;因此高程在2~4km、坡度在30°~60°和坡向在北坡、东北坡和东南坡的区域为最易受地震影响的区域,加强该区域的水土保持工作是非常有意义的。
本研究可初步得出震后水土保持功能易受损区域的空间分布规律,为震后水土流失防治提供指导,对震后灾区尤其是受灾景区的水土保持工作提供科学依据。