2008年5月12日的汶川大地震主要发生在龙门山断裂带，地震灾区南北跨度为400 km，东西宽约70 km，该区域是青藏高原高寒植被区与东部湿润森林区的结合部，土地利用以林业用地为主，森林覆盖率达44.51%(Wu et al.，2012)。据调查，地震造成重灾区受损林地面积达32.87万hm²，灾害类型包括山体崩塌、滑坡、坠石、堰塞湖和泥石流，导致该区域森林覆盖率下降1.87%(骆建国等，2011)。由于山地地质构造复杂、生态系统本身较为脆弱，震后往往会发生滑坡、泥石流、塌方甚至山崩等次生灾害，进一步造成森林植被的破坏及生境破碎化程度的加剧，严重影响到该地区生态环境和动植物的生存(马文宝等，2012; Cheng et al.，2012)，因此，加速地震灾区的林业生态修复，对保障该区域生态安全、构筑长江上游生态屏障具有重要现实意义。近5年来，关于汶川地震在土壤侵蚀、环境次生灾害、重金属污染等方面的影响已有不少调查与报道。杨渺等(2013)采用GIS/RS技术监测了地震灾区石亭江上游474.41 km²的典型流域水土保持功能恢复效应，结果表明2011年震区土壤侵蚀总量比2008年减少11.92%，但土壤侵蚀强度的面积分布格局并未改变，灾区生态安全仍令人担忧。Xu等(2012)通过研究震后土壤重金属Cd，Pb空间分布特征及影响因素，认为地震导致土壤重金属污染程度加剧。地震后发生的次生灾害对土壤动物，尤其是大型土壤动物的群落结构及多样性有强烈的干扰作用(吴鹏飞等，2010)。研究表明地震灾区土壤有机质含量能较好地表征土壤抗蚀性，有机质含量的减少是导致土壤水土流失的主要因素之一(宋旭等，2010; 蔡艳等，2010)。目前，对地震引起的森林土壤系统的变化研究大多数都是从宏观层面上对水土流失和植被恢复的过程和变化开展比较和监测的(郭兵等，2012; Wu et al.，2012)，而对地震后森林土壤系统的组成、结构和功能的变化过程研究还很缺乏，导致对灾后森林土壤演变的途径、机制以及特点尚不清楚。本研究通过对地震灾区四川理县典型山地岷江柏(Cupressus chengiana)人工林土壤理化性质变化的调查和比较，进一步了解地震对森林土壤环境带来的影响及可能存在的生态后果，为震后森林土壤保护和生态修复提供理论依据。

理县地处岷江上游杂谷脑河流域(102°33'—103°30'E，30°54'—31°12'N)，地质结构属龙门山断裂带中段，地壳活动频繁，地质构造复杂，海拔1 442~3 400 m，平均坡度30°~40°，属典型的高山峡谷地貌(郭兵等，2012)。年均气温6~9 ℃，年均降水量400~700 mm，年均蒸发量1 400~1 700 mm，属典型的干旱河谷气候。

本研究选取理县的熊尔山和蒲溪沟2个典型受灾区的岷江柏人工林，土壤类型分别为山地褐土和山地棕壤，成土母质分别为碳酸岩和砂页岩。选择海拔高于2 000 m、坡度大于30°、东南/西南坡向的岷江柏人工林作为受灾调查林分。岷江柏是当地的乡土树种，也是我国特有的珍稀树种，耐寒、抗旱、根深，具有良好的水土保持和水源涵养功能。调查林分林龄为8~12年生，种植密度为1.5 m×1.5 m。岷江柏群落结构比较简单，层次明显。熊尔山调查点的主要林下植被有中华胡枝子(Lespedeza chinensis)、西南苔草(Carex austro-occidentalis)、虎榛子(Ostryopsis davidiana)和阿坝蒿(Artemisia abaensis)等。蒲溪沟调查地的主要林下植被有秦岭槲蕨(Drynaria sinic)和川甘亚菊(Ajania potaninii)等。调查区域约有近半数的森林在汶川地震中受到严重影响，土体出现翻动、崩塌、滑坡、泥石流等现象。同时，分别在熊尔山和蒲溪沟分别选取1个与调查林分立地条件相近、位置条件相对一致但没有明显的断裂、颠覆和滑坡等土地破坏现象，土体及土表植被相对完整地段作为对照地，编号分别为CK1，CK2。在熊尔山研究区选取的3个受害调查地编号为S1，S2，S3，在蒲溪沟研究区选取的2个受害调查地编号为P1，P2。样地具体情况见表 1。

表 1 调查地基本情况 Tab.1 Basic situation of study sites

表 1 调查地基本情况 Tab.1 Basic situation of study sites 调查地 Study site 土壤类型 Soil type 面积 Area/hm2 海拔 Altitude/m 坡向及坡度 Aspect, slope 林分郁闭度 Canopy density 平均树高 Mean tree height/m DBH/cm 地表覆盖度 Coverage (%) 受灾类型 Type of disaster-affected CK1 山地褐土 Mountain cinnamon soil 3.50 2 350 东南 Southeast,30° 0.50 5.05 8.65 25 无明显受灾现象 No obvious disaster- affected phenomena S1 山地褐土 Mountain cinnamon soil 25.00 2 460 东南 Southeast,35° 0.17 3.68 5.21 14 山体崩塌 Mountain collapse S2 山地褐土 Mountain cinnamon soil 2.75 2 420 东南 Southeast,30° 0.16 3.56 5.04 13 土壤滑移 Soil sliding S3 山地褐土 Mountain cinnamon soil 14.79 2 410 东南 Southeast,30° 0.18 3.92 5.56 15 土壤侵蚀 Soil erosion CK2 山地棕壤 Mountain brown soil 5.50 2 370 西南 Southwest,30° 0.47 5.13 8.70 23 无明显受灾现象 No obvious disaster- affected phenomena P1 山地棕壤 Mountain brown soil 4.82 2 300 西南 Southwest,34° 0.14 3.14 4.15 12 土体崩塌、掩埋 Soil collapse and burial P2 山地棕壤 Mountain brown soil 9.53 2 480 西南 Southwest,35° 0.13 3.06 3.89 11 滑坡、泥石流 Landslide and debris flow

于2013年11月份进行土壤采样和调查。在每个调查地分别选择3块30 m×30 m的典型林分作为调查的标准样地，采用梅花形的布点方法设置5个取样点。表层土壤(0~20 cm)用100 cm³环刀取土，保持环刀土原状带回室内用于土壤密度、总孔隙度和颗粒组成的测定。每块样地挖5个土壤剖面，按0~20，20~40和40~60 cm土层分别取土，各土层分别采1 kg土样，装入自封袋带回室内，将土样铺平到牛皮纸上进行自然风干，拣出石砾和植物残体，然后将风干土研碎，充分混匀后用四分法依次过2，1和0.149 mm筛，分装入新的自封袋编号，进行土壤酸碱度、阳离子交换量、有机质含量、全氮含量、碱解氮含量、有效磷含量和速效钾含量等化学性质的分析测定。各林地土壤的理化性质分析数值用5个采样点的数学平均值±标准差表示。

土壤密度采用环刀法测定(NY/T 1121.4-2006); 土壤总孔隙度采用水浸泡-环刀法测定(LY/T 1215-1999); 土壤pH值采用电极电位法(水土比为2.5∶1)测定(NY/T 1377-2007); 土壤机械颗粒组成采用激光粒度仪(济南微纳)分析; 土壤全氮含量采用凯氏定氮仪法测定(LY/T 1228-1999); 碱解氮含量采用碱解扩散法测定(鲁如坤，2000); 有效磷含量采用0.5 mol·L^-1 NaHCO₃浸提—分光光度计比色法测定(鲁如坤，2000); 速效钾含量采用NH₄OAc浸提—火焰光度计法测定(鲁如坤，2000); 土壤有机质含量采用重铬酸钾氧化—外加热法测定(鲁如坤，2000); 土壤阳离子交换量采用乙酸钠-火焰光度计法测定(GB 7857-1987); 土壤黏粒的提取采用虹吸法(LY/T 1252-1999); 黏土矿物组成含量采用X射线衍射分析(SY/T 5163-1995): 称取40 mg干样制成定向片，在X射线衍射仪上分析，测试条件为Cu-Kα辐射，扫描速度采用4°·min^-1，工作电压40 kV，电流40 mA。

采用Microsoft Excel 2010软件进行数据计算处理; 采用SPSS16.0软件进行统计分析，对不同立地或不同理化性质的分析用LSD法进行多重比较(P＜0.05)和Pearson相关分析。

表 2为调查地0~20 cm土层土壤的物理性质，熊尔山研究区的土壤密度为1.15~1.18 g·m^-3，蒲溪沟研究区土壤密度为1.35~1.52 g·m^-3，熊尔山土壤密度小于蒲溪沟。受地震次生灾害影响较严重的蒲溪沟林地P1和P2的土壤总孔隙度显著减少，分别低于对照21.6%和9.8%。总孔隙度减少将导致土壤下渗能力降低，容易造成水土流失。0~20 cm土层土壤机械组成分析表明，S1，S2，S3的土壤黏粒成分相对于CK1分别降低52.0%，55.1%和57.6%，粉粒分别升高12.5%，17.2%和16.2%，砂粒分别升高13.6%，10.3%和12.2%; P1，P2的土壤黏粒成分相对于CK2分别降低54.8%和46.5%，粉粒分别升高6.5%和8.5%，砂粒分别升高31.7%和17.5%。地震引发的一系列次生灾害，如山体崩塌、泥石流、土层断裂、土壤侵蚀、滑坡和掩埋等，造成土壤中砂粉粒含量显著增多，黏粒含量显著减少。

表 2 土壤物理性质^① Tab.2 Physical properties of soil

表 2 土壤物理性质① Tab.2 Physical properties of soil 调查地 Study site 土壤密度 Soil density/(g·cm-3) 土壤总孔隙度 Total porosity(%) 土壤机械组成 Soil mechanical composition 黏粒Clay＜0.002 mm(%) 粉粒Silt 0.002~0.05 mm(%) 砂粒sand＞0.05 mm(%) CK1 1.10±0.04a 59.0±1.5a 19.8±1.2a 39.5±1.0a 40.7±0.5a   S1 1.15±0.11b 56.5±4.0ab 9.5±0.7b 44.4±1.0b 46.1±1.7b   S2 1.18±0.11b 55.5±4.0b 8.9±0.6b 46.2±0.5b 44.9±0.1b   S3 1.18±0.09b 55.6±3.4b 8.4±0.7b 45.9±0.6b 45.7±1.3b CK2 1.21±0.03a 54.2±1.3a 21.3±1.4a 41.4±0.5a 37.2±0.9a   P1 1.52±0.06b 42.5±2.2c 6.9±0.4b 44.1±0.4b 49.0±0.2bc   P2 1.35±0.09b 48.9±3.4b 11.4±1.0b 44.9±0.9b 43.7±1.6b ①同列不同小写字母表示差异显著(P＜0.05)。Different lowercase in the same column represents significant difference (P<0.05).

土壤黏粒的XRD分析表明，2个调查点所有林地土壤主要黏土矿物组成接近(表 3)，各调查地的伊利石含量高达69.8%~79.9%，其次为绿泥石和高岭石，伊蒙混层含量很少，无蒙脱石，黏土矿物的组成未受到地震次生灾害明显影响。

表 3 黏土矿物组成 Tab.3 Composition of clay minerals

表 3 黏土矿物组成 Tab.3 Composition of clay minerals 调查地Study site 伊利石Illite (%) 绿泥石Chlorite (%) 高岭石Kaolinite (%) 伊蒙混层Illite/smectite formation (%) CK1 71.1 14.7 11.0 3.2   S1 73.0 13.1 11.2 2.7   S2 79.9 19.7 14.1 6.3   S3 79.4 11.4 8.7 0.5 CK2 72.2 22.8 5.9 0.9   P1 69.8 21.6 7.7 0.9   P2 77.7 20.2 1.8 0.3

由表 4可知，受地震次生灾害影响土壤全氮含量为0.5~3.3 g·kg^-1，碱解氮含量为11.5~242.4 mg·kg^-1。全氮和碱解氮含量在不同调查地不同土层的变化幅度较大。与相应对照相比较，S1，S2和S3的全氮含量分别降低35.5%，71.0%和41.9%，碱解氮含量分别降低了40.8%，78.7%和53.8%; P1和P2的全氮含量分别降低85.7%和28.6%，碱解氮含量分别降低93.7%和46.1%。S1，S2和S3调查地的有效磷含量均随土层深度增加而降低，P1和P2有效磷含量并没有出现这一趋势，P2的有效磷含量是所有调查地中最高的，在40~60 cm土层达到最大值10.8 mg·kg^-1。P1和P2的无规则变化和异常值可能与地震次生灾害影响较大导致土层变动有关。与对照相比，震后土壤的3个土层中的有效磷含量均显著降低。调查地土壤速效钾含量总体较低，对照地为17.5~33.5 mg·kg^-1，而调查地速效钾含量最低为1.7 mg·kg^-1。S1，S2和S3的速效钾含量分别比CK1降低64.4%，80.9%和78.5%; P1和P2的速效钾含量分别比CK2降低了80.3%和81.6%。

表 4 土壤氮、磷、钾养分含量^① Tab.4 Nitrogen, phosphorus, and potassium contents of soil

表 4 土壤氮、磷、钾养分含量① Tab.4 Nitrogen, phosphorus, and potassium contents of soil 土层 Soil layer 调查地 Study site 全氮 Total nitrogen/(g·kg-1) 碱解氮 Alkali-hydrolysable nitrogen/(mg·kg-1) 有效磷 Available phosphorus/(mg·kg-1) 速效钾 Available potassium/(mg·kg-1) ①同一土层不同调查地不同小写字母表示差异显著(P＜0.05)。下同。Different lowercase of different study site in the same soil layer represents significant difference (P<0.05). The same below. 0~20 cm CK1 3.6±0.10a 347.6±1.24a 15.6±0.64a 33.5±0.95a   S1 2.4±0.10b 242.4±0.61b 6.4±0.57b 10.7±1.01b   S2 1.1±0.09c 98.4±0.99c 2.5±0.99c 7.7±0.33c   S3 2.2±0.10b 200.9±1.48b 6.5±0.41b 9.5±0.98bc CK2 5.4±0.26a 414.9±11.62a 15.2±0.70a 30.4±0.94a   P1 0.5±0.10c 24.1±1.02c 1.5±0.31c 8.0±1.02b   P2 3.3±0.11b 139.6±0.54b 10.8±0.62b 5.4±0.39c 20~40 cm CK1 3.2±0.10a 293.7±1.02a 13.9±0.99a 26.1±1.02a   S1 2.4±0.10b 192.2±0.54b 3.7±1.04b 9.6±1.27b   S2 0.8±0.07c 51.6±1.02c 1.6±0.17b 4.0±0.54c   S3 1.6±0.13bc 117.4±1.11b 5.4±0.43b 5.1±1.07c CK2 3.4±0.18a 351.8±6.88a 13.2±0.60a 25.1±0.65a   P1 0.5±0.12c 11.5±1.01c 1.4±0.43c 4.7±0.23b   P2 2.2±0.06b 195.9±0.98b 7.7±0.24b 4.5±0.43b 40~60 cm CK1 2.6±0.06a 250.8±0.53a 11.4±0.63a 18.9±0.98a   S1 1.3±0.10b 93.5±0.98b 2.3±0.82bc 7.2±1.00b   S2 0.7±0.10c 40.0±1.07c 1.2±0.27c 3.6±0.99c   S3 1.5±0.09b 93.6±0.98b 2.7±0.22b 3.2±1.03c CK2 2.4±0.23a 228.8±10.84a 10.0±0.54a 17.5±0.67a   P1 0.6±0.11b 26.8±1.03c 1.5±0.30b 1.7±0.41c   P2 2.5±0.12a 200.7±1.01b 10.7±0.99a 3.5±0.39b

由表 5可知，受地震次生灾害影响，土壤pH值呈现上升趋势。熊尔山研究区受害样地的不同土层pH值均大于8.0，属于碱性土壤; 蒲溪沟研究区的土壤碱性与对照相比也明显增强，土壤平均pH值提高了11.5%。调查地土壤有机质含量低至4.4 g·kg^-1，与对照相比，S1，S2和S3的有机质含量分别降低48.6%，88.6%和50.0%，P1和P2分别降低89.2%和11.1%。震后土壤的阳离子交换量最大值为5.88 cmol·kg^-1，最小值为1.74 cmol·kg^-1，平均值为3.72 cmol·kg^-1。与对照相比，S1、S2和S3的阳离子交换量分别降低87.5%，89.3%和85.9%，P1和P2分别降低90.3%和79.7%。一般认为土壤阳离子交换量低于10 cmol·kg^-1为保肥力弱，由此看来，地震后土壤的保肥能力都很弱。2个调查地中受灾林分与对照林分相比，土壤碳氮比虽有差异但并没有呈现规律性变化; 且同一样地不同土层的土壤碳氮比也未发现规律性变化。但是可以看出所有样地的碳氮比值均小于12。

表 5 土壤pH值、阳离子交换量、有机质含量及C/N比 Tab.5 Soil pH value, cation exchange capacity, soil organic matter and C/N ratio

表 5 土壤pH值、阳离子交换量、有机质含量及C/N比 Tab.5 Soil pH value, cation exchange capacity, soil organic matter and C/N ratio 土层Soil layer 调查地Study site pH 土壤有机质含量Soil organic matter/(g·kg-1) 阳离子交换量Cation exchange capacity/(cmol·kg-1) 碳氮比C/N ratio 0~20 cm CK1 7.4±0.21a 76.3±5.16a 34.9±3.04a 12.3±0.91a   S1 8.2±0.21b 46.3±2.62b 4.2±0.92b 11.2±0.39b   S2 8.3±0.33b 11.5±1.74c 4.6±0.38b 6.0±1.00c   S3 8.3±0.33b 42.6±2.66b 4.5±0.39b 11.8±0.74ab CK2 5.4±0.07a 83.3±4.35a 35.6±0.43a 13.5±0.75a   P1 6.2±0.34b 4.4±0.40c 2.6±0.38c 5.3±0.58b   P2 6.3±0.32b 75.3±5.36b 5.9±1.04b 13.1±0.88a 20~40 cm CK1 7.5±0.09a 64.9±4.38a 28.2±2.08a 7.5±0.09a   S1 8.4±0.12b 31.2±3.60c 3.7±0.27b 7.6±0.87c   S2 8.3±0.28b 5.3±1.61d 2.6±0.48c 3.8±1.19d   S3 8.4±0.36b 27.8±1.73c 3.6±0.43b 10.3±0.65b CK2 5.6±0.04a 68.6±2.04a 27.1±0.69a 13.0±0.35a   P1 6.6±0.26b 10.1±1.44c 2.4±0.33c 12.8±2.56a   P2 6.5±0.33b 56.0±4.38b 5.1±0.99b 14.2±1.13a 40~60 cm CK1 7.7±0.10a 50.1±6.98a 20.8±3.05a 10.9±1.71a   S1 8.6±0.12b 21.0±2.66b 2.6±0.45c 9.7±1.10a   S2 8.3±0.15bc 5.0±1.47c 1.7±0.43c 4.3±1.29b   S3 8.4±0.39bc 25.3±1.80b 3.7±0.23bc 10.0±0.69a CK2 5.7±0.02a 48.2±1.93a 19.5±0.92a 10.5±0.49a   P1 6.6±0.27b 7.2±1.68c 3.0±1.02c 8.5±1.98b   P2 6.8±0.11b 46.6±5.62b 5.7±0.93b 10.9±1.37a

相关分析得到: 土壤有机质与全氮含量(n=105，r=0.956，P＜0.01)与碱解氮含量(n=105，r=0.784，P＜0.01)、有效磷含量(n=105，r=0.947，P＜0.01)和阳离子交换量(n=105，r=0.850，P＜0.01)均极显著正相关，与速效钾含量的相关性不显著。

本研究调查发现，研究区下层土体中碎石、砾石含量高达10%~30%，因此，只能采集0~20 cm土层土壤开展物理性质分析，未对20 cm以下土层用环刀精确地采集土壤样品。庞学勇等(2004)曾对该地27个岷江柏林下土壤密度开展调查，0~20 cm层的变化范围为0.87~1.13 g·cm^-3，而本研究调查得到的土壤密度较高，尤其是在蒲溪沟发育于页岩的土壤，该点对照土壤的密度高达1.21 g·cm^-3。地震次生灾害导致土壤密度显著提高，这一结果同李志华等(2014)的无林地中调查结果不一致，他们发现地震造成无林地的土壤密度显著降低。这一现象可能与土地利用方式有关，林地土壤地下部分根系深且密度大，再加上地上部分生物量高，对土体有一定的挤压作用(Nancy et al.，2012)，因此本研究中受灾林地的密度增大。有些森林土壤研究中发现，土壤密度常常同土壤有机质含量、土壤粗颗粒(粉粒、沙粒)含量有关(Nancy et al.，2012)，但在本调查中没有得到类似关系，这可能同调查地土壤质地粗、样地数量少有关。本研究发现土壤密度增大直接导致土壤孔隙度下降、通透性降低，因此产生地表径流，引起水土流失和肥力下降，并发生不同程度的泥石流(王东升等，2011)。

受地震次生灾害严重破坏土壤的砂粒和粉粒含量明显增多，黏粒含量急剧减少。该地区土壤黏粒的主要成分是伊利石和绿泥石，主要与成土母质类型有关(刘学录，2011)。伊利石具有一定的比表面积和大量负电荷，能通过与土壤有机质的作用形成团聚体结构(袁仁茂等，2013)。但是地震造成土壤团聚结构体破坏，导致土壤颗粒缺乏有机质的胶结，进而土壤抗侵蚀能力减弱，发生黏粒淋失(薛立等，2011)。土壤黏粒含量与有机质、全氮等养分含量呈正相关，土壤黏粒的减少直接导致土壤养分减少，肥力下降(王文艳，2013)。

岷江柏生长的土壤环境为碱性或强碱性，林下土壤pH值为6.00~9.14(庞学勇等，2005)，本研究选取的调查林地的pH在此范围中，因此所选调查地具有较强的代表性。本调查点2种土壤类型(山地褐土和山地棕壤)的土壤pH值受地震次生灾害影响与对照相比出现升高趋势，并随土层深度增加而升高。这一现象同土壤中有机质和黏粒流失有关，森林土壤中有机质常在盐基离子含量少的土壤中呈现酸性或偏酸性(Nancy et al.，2012)，但是调查地土壤瘠薄，有机质大量淋失造成了土壤pH值改变。岷江柏分布区0~20，20~40和40~60 cm土层有机质含量的平均值分别为(90.06±8.72)，(38.31±4.25)和(32.58±6.99)g·kg^-1(庞学勇等，2005)，与本研究选取的对照地的有机质含量相近，而地震破坏林地的不同土层有机质含量差异显著，土壤有机质含量显著降低且随土层深度增加无变化规律，主要是地震扰动造成土壤剖面结构改变(Wu et al.，2012)。该地区土壤土层较薄，泥石流及滑坡等造成岩石裸露，地表裸露后又没有足够的凋落物和植物残体补充，也是造成有机质缺乏的原因(刘旦旦等，2013; 文炯等，2009)。调查地土壤有机质含量与全氮含量、碱解氮含量、有效磷含量和阳离子交换量之间均显著正相关，说明有机质含量对土壤理化特性影响很大(文炯等，2009)。调查地的CEC均小于10 cmol·kg^-1，说明保肥能力很弱。除速效钾外，研究林地的土壤养分含量在土壤剖面均表现为0~20 cm>20~40 cm＞40~60 cm土层，这与庞学勇等(2005)的研究结果相同。C/N比是衡量土壤C，N营养平衡状况的指标，土壤微生物获得平衡营养的C/N比约为25(林青等，2011)。但是调查地的土壤C/N比大部分在12以下，说明地震造成土壤的碳损失率大于氮损失率，这与前人的研究一致(Wu et al.，2012)。

生态修复的目标是通过人工设计和恢复措施，在受地震干扰破坏的生态系统的基础上，恢复和重新建立一个具有自我恢复能力的健康的生态系统(杨立东，2011)。调查地土壤存在明显的障碍因子: 土壤有机质含量减少，碳酸岩母质上发育的土壤碱性严重升高，表层土壤质地较粗，黏粒含量少。围绕上述障碍因子提出以下修复建议: 1)增加土壤有机质。根据研究地属高坡度林地的实际情况建议种植适于当地干冷环境的绿肥品种如鼠茅草(Vulpia myuros)、毛苕子(Vicia villosa)和苜蓿(Medicago sativa)等，达到有机质长期有效输入的目的; 2)改良黏土矿物。补充土壤黏粒，优化土壤颗粒组成，适当添加缓冲性能强的黏土矿物如蛭石和凹凸棒土等，促进土壤有机无机复合体的形成，从而提高土壤阳离子交换量;3)对碳酸岩发育土壤，应改善土壤pH值。施用酸性有机材料如褐煤(pH5.6)等修复该地土壤，可以在一定程度上降低土壤pH值，其巨大的比表面积和丰富的孔状结构也能提高土壤通气性; 4)降低人为干扰。在进行土壤修复的林地设置警示牌和栅栏，避免伐木、放牧和施工等人为活动的干扰。

综上所述，汶川地震对理县地区人工岷江柏林下土壤理化性质具有很强的破坏影响，主要受灾类型为山体崩塌、滑坡、泥石流和土体断裂及掩埋。岷江柏林下的山地褐土和山地棕壤均出现土壤碱性增强、密度升高的现象，并且土壤黏粒含量大幅度降低，砂粒和粉粒明显增多。地震在很大程度上引起了土壤剖面扰动，导致不同采样地不同土层土壤养分含量出现较大波动和明显降低，受灾林地土壤全氮含量、碱解氮含量、有效磷含量、速效钾含量、有机质含量和阳离子交换量均出现显著下降。因此，该地区应该加强对地震次生灾害破坏的森林土壤生态修复的重视，尽快改善岷江柏林下土壤环境、提高林地生产力。