Origin and stability of landslides in Chaya County, Lancang River Basin, Tibet
-
摘要:
地球内外动力耦合叠加人类工程活动下的大型滑坡灾害是青藏高原工程建设面临的工程地质问题之一。针对在澜沧江流域察雅县城南侧斜坡新发现的9处大型—特大型滑坡，采取野外调查、试验测试和数值计算，分析了滑坡的特征、成因和现状稳定性。结果表明：①察雅县城滑坡具有规模大、切割深、地形阻止系数大等特征，滑坡堆积体现状具有多期蠕滑变形，表现为牵引式滑移的变形破坏特点，再次活动主要与降雨、河流侵蚀、人类工程活动等因素有关；②综合滑坡的发育特征、运动特征、物质结构与组成、断层与历史地震4个方面，初步判断察雅县城滑坡的诱发因素为内动力地质作用，为历史地震滑坡；③古妥尔滑坡堆积体在天然状态下处于基本稳定-稳定状态，在强降雨状态下滑坡前缘局部或整体处于欠稳定状态，较易-易发生滑动。新发现的大型—特大型滑坡拓宽了青藏高原"三江"地区地震滑坡的分布范围，对研究昌都及周边地区的地震历史、断裂活动性、地貌形成演化等方面提供了重要的佐证资料，滑坡堆积体现状稳定性分析结果可为县城的扩建及地质灾害风险管控提供一些参考。
Abstract:Large-scale landslide disasters under the coupling of internal and external dynamics of the earth and human engineering activities are one of the engineering geological problems faced by major engineering constructions on the Qinghai-Tibet Plateau. Based on the field investigation, experimental test and numerical calculation of the newly discovered 9 large-extra large landslides on the southern slope of Chaya County in the Lancang River Region, the characteristics, causes and current stability of the landslides were analyzed. The results show that the Chaya landslide has the characteristics of large scale, deep cutting, and large topographic blocking coefficient; the landslide accumulation is represented by multi-stage creep deformation characterized by traction slip; and the deformation is mainly related to factors such as rainfall, river erosion and human engineering activities. The integration of the development characteristics, movement features, material structure and composition, faults and historical earthquakes of the landslides indicates that the inducing factor of the landslide in Chaya County is internal dynamic earthquake action, and it is a historical earthquake landslide. The Guttor landslide accumulation body is basically stable or stable in its natural state, and the front edge of the landslide is in a partially or overall under-stable state under heavy rainfall conditions, which is relatively easy to slip. The newly discovered large-extra-large ancient landslide provides important supporting information for the study of the earthquake history, fault activity and geomorphological evolution of Changdu and surrounding areas. The evaluation results of current stability of landslide accumulations can provide some scientific basis for the expansion of the county and the risk management as well as control of geological disasters.
-
中亚造山带是显生宙以来全球最大的陆壳增生与改造地带，其主体是由古亚洲洋岩石圈板片从中元古代晚期开始，不断向两侧古陆之下俯冲而形成岛弧和增生杂岩，并最终闭合使得两侧陆缘碰撞而形成的[1-7]。位于内蒙古东南部的中亚造山带南缘，通常被认为是中朝古板块与西伯利亚古板块之间的古亚洲洋海洋盆地最后消失的地区[3, 8-10] (图 1-a)。近年，综合古地磁、古生物、岩浆岩等关键性地质记录资料，夹于中朝古板块和西伯利亚古板块之间的古亚洲洋最终沿索伦-西拉木伦河缝合带闭合的观点得到了许多学者的认同[4-5, 9-23]。然而，由于对地质记录不同的解释及残余洋壳有限的保留和出露，使古亚洲洋最终闭合的时间及该缝合带的东延位置仍存在很大不确定性。研究区位于中亚造山带南缘，紧邻松辽盆地西缘，靠近索伦-西拉木伦河缝合带向东延伸的大致部位，是研究该缝合带东延进入盆地的关键部位。项目组在区域地质调查工作期间，为了对古亚洲洋最终闭合的位置和时间给予进一步的约束，对西拉木伦河一线的岩石和构造进行了较系统的观察和研究，获得了一些新的发现。在研究区原定分布较广泛的燕山期花岗岩中识别出早二叠世岩浆作用的记录，通过岩石地球化学和锆石Hf同位素研究，结合区域地质资料及前人研究成果，对该岩体成因进行了探讨，以期对该区晚古生代构造演化及缝合带的位置提供约束。
![]()
1. 地质背景
研究区位于内蒙古赤峰市阿鲁科尔沁旗南部地区，隶属大兴安岭南段东坡，毗邻松辽盆地西缘，大地构造位置上处于传统划定的二连浩特-贺根山缝合带以南、索伦-西拉木伦河缝合带偏北，一定程度上保存了中国东北地区古生代洋盆最终消失的地质记录(图 1-a)。研究区以南为中朝古板块及其北缘中段的古生代活动陆缘增生区，两者以赤峰-巴彦敖包断裂为界，以发育近东西向分布的早古生代俯冲增生杂岩和弧岩浆岩，以及具规律性成因的晚古生代活动陆缘型侵入岩和火山-沉积岩系为特征[3, 5, 10-11, 24-27]；而中朝古板块北缘中段通过索伦-西拉木伦河缝合带与西伯利亚古板块南缘相接，该区在显生宙期间经历了地壳显著地同造山水平生长和后造山垂向增生，以发育岛弧增生杂岩、多类型蛇绿岩及巨量新生地壳熔融而成的岩浆岩为特征，可能是西伯利亚南缘面向类似现今西太平洋多岛洋盆的格局而长期演化至洋盆消亡的结果[28-31]。
研究区出露的最老地层单元为零星分布于宝日乌苏镇—天山口镇一带的石炭系，主要由片岩、石英岩和结晶灰岩构成，发育中、小型石灰岩工业矿床；而在大兴安岭南缘广泛发育的二叠系在研究区也有出露，为大石寨组(原定青凤山组)海相中基性—酸性火山岩，其与下伏石炭系为断层接触；侏罗系在区内出露最广泛，呈面状分布，主要由陆相中酸性火山岩和砂页岩等碎屑岩构成，其中中侏罗世碎屑岩中产小型煤矿；白垩系在研究区分布有限，出露于阿鲁科尔沁旗北部福兴屯地区，以灰黑色辉石安山岩为主；而新生界在区内以第四纪冲积砂砾层、风成沙等松散堆积物为主，含少量河谷橄榄玄武岩[32]。区域地质调查报告①显示，研究区没有出露早于中生代的岩浆作用产物，而是主要记录了燕山中期以来的岩浆活动，以小规模的侏罗纪—白垩纪中酸性侵入体和火山作用为主，但从本文的年代学研究看，研究区至少经历过海西期岩浆活动。此外，研究区遭受了多期构造活动，断裂构造较发育。
2. 岩体野外地质和岩石学特征
安乐屯岩体呈小岩株出露，面积约5 km2(图 1-b)。野外露头观察表明，岩体经历了不同程度的风化作用和地质构造活动，节理较发育，可见后期酸性或基性岩脉侵入，同时岩体整体发生了变形，可见片理化现象。岩体与晚侏罗世下兴安岭组中酸性火山岩呈断层接触，并被燕山中期花岗岩侵入，但围岩不明①。岩体样品采样位置位于安乐屯以东露头良好的小山顶部附近，露头尺度可见组成岩体的岩石类型发生变化(图版Ⅰ-a、b)，经岩石薄片显微镜下观察，该岩体主要由灰白色中细粒石英二长闪长岩(BL03-4和BL03-7)-花岗闪长岩(BL03-2、BL03-3和BL03-5)-黑云母二长花岗岩(BL03-1)组成(图版Ⅰ-c~f)。岩石样品整体上呈中细粒半自形结晶结构，片麻状构造，但主要造岩矿物组成及含量略有不同：黑云母二长花岗岩主要由石英(约28%)、斜长石(约35%)、碱性长石(约30%)和黑云母(约5%)组成(图版Ⅰ-c)；花岗闪长岩主要由石英(约25%)、斜长石(约45%)、碱性长石(约20%)和黑云母(约8%)组成(图版Ⅰ-d)；石英二长闪长岩主要由碱性长石(约35%)、斜长石(约40%)、石英(约10%)、黑云母(约10%)和角闪石(约2%)组成(图版Ⅰ-e、f)。此外，岩石中包含少量锆石、磷灰石及磁铁矿等副矿物(2%~3%)，局部发生了绿帘石化和绿泥石化，而且石英二长闪长岩中的角闪石普遍存在暗化现象。
![]()
图版Ⅰa.安乐屯岩体花岗岩野外露头；b.安乐屯岩体石英二长闪长岩野外露头；c、d.花岗岩样品显微照片(正交偏光)；e.石英二长闪长岩样品显微照片(单偏光)；f.石英二长闪长岩样品显微照片(正交偏光)。Af—碱性长石；Bi—黑云母；Pl—斜长石；Q—石英；Mt—磁铁矿图版Ⅰ.3. 分析方法
3.1 锆石U-Pb定年
样品先后采用浮选和电磁选完成单矿物筛选，锆石的挑选由河北省区域地质调查大队地质实验室完成。锆石的阴极发光(CL)图像采集和U-Pb同位素分析均在中国地质大学(武汉)地质过程与矿产资源国家重点实验室完成。实验室应用Agilent 7500a ICP-MS仪器，按照标准测定程序进行锆石U-Pb测定。实验中采用高纯氦为剥蚀物质载气，用美国国家标准技术研究院研制的人工合成硅酸盐玻璃标准参考物质NIST SRM610进行仪器最佳化，样品测定时用哈佛大学标准锆石91500为外部校正标样，以监测标样和样品的仪器条件是否一致。实验中采用8 Hz激光频率、52 mJ激光强度和32 μm激光束斑直径。样品同位素数据处理采用ICPMSDataCal[33-34]，运用Andersen方法进行同位素比值校正，以扣除普通铅对定年结果的影响[35]。年龄计算和谐和图的绘制均采用国际标准程序Isoplot(ver3.0)完成[36]。详细的实验步骤和数据处理方法见参考文献[34]。锆石LA-ICP-MS U-Pb分析结果见表 1。
表 1 内蒙古阿鲁科尔沁旗安乐屯岩体锆石U-Th-Pb分析结果Table 1. Zircon U-Th-Pb data of the Anletun pluton in Ar Horqin Banner, Inner Mongolia样品号 含量/10-6 Th/U 比值 年龄/Ma Th U 207Pb/235U 1σ 206Pb/238U 1σ 208Pb/232Th 1σ 207Pb/235U 1σ 206Pb/238U 1σ 208Pb/232Th 1σ BL03-1-01 89 179 0.50 0.32348 0.01240 0.04348 0.00043 0.01316 0.00036 285 10 274 3 264 7 BL03-1-02 64 137 0.46 0.29881 0.01647 0.04360 0.00054 0.01315 0.00044 265 13 275 3 264 9 BL03-1-03 48 111 0.43 0.31618 0.01866 0.04401 0.00061 0.01487 0.00054 279 14 278 4 298 11 BL03-1-04 76 145 0.53 0.30907 0.01331 0.04296 0.00056 0.01097 0.00038 273 10 271 3 220 8 BL03-1-05 244 352 0.69 0.33950 0.01024 0.04340 0.00037 0.01465 0.00035 297 8 274 2 294 7 BL03-1-06 72 152 0.47 0.32433 0.01200 0.04344 0.00049 0.01332 0.00043 285 9 274 3 268 9 BL03-1-07 70 152 0.46 0.34798 0.01548 0.04371 0.00051 0.01530 0.00050 303 12 276 3 307 10 BL03-1-08 191 292 0.65 0.31868 0.01265 0.04388 0.00041 0.01375 0.00034 281 10 277 3 276 7 BL03-1-09 65 149 0.44 0.30551 0.01408 0.04278 0.00049 0.01536 0.00055 271 11 270 3 308 11 BL03-1-10 122 190 0.64 0.30211 0.01094 0.04401 0.00041 0.01398 0.00033 268 9 278 3 281 7 BL03-1-11 93 164 0.57 0.30695 0.01392 0.04334 0.00043 0.01423 0.00038 272 11 274 3 286 8 BL03-1-12 90 187 0.48 0.31149 0.01194 0.04361 0.00043 0.01402 0.00035 275 9 275 3 281 7 BL03-1-13 104 176 0.59 0.31045 0.01361 0.04418 0.00049 0.01312 0.00033 275 11 279 3 263 7 BL03-1-14 135 288 0.47 0.35014 0.01089 0.04836 0.00050 0.01546 0.00035 305 8 304 3 310 7 BL03-1-15 139 178 0.78 0.36058 0.01481 0.04334 0.00053 0.01598 0.00052 313 11 274 3 320 10 BL03-1-16 61 125 0.49 0.31625 0.01697 0.04389 0.00052 0.01437 0.00055 279 13 277 3 288 11 BL03-1-17 127 268 0.48 0.31873 0.01430 0.04333 0.00042 0.01397 0.00039 281 11 273 3 280 8 BL03-1-18 52 132 0.39 0.32111 0.01713 0.04356 0.00053 0.01525 0.00052 283 13 275 3 306 10 BL03-1-19 66 134 0.49 0.32898 0.01442 0.04342 0.00046 0.01419 0.00048 289 11 274 3 285 10 BL03-1-20 84 159 0.53 0.33255 0.01287 0.04385 0.00045 0.01556 0.00138 292 10 277 3 312 27 BL03-1-21 102 203 0.50 0.31970 0.01225 0.04420 0.00041 0.01409 0.00034 282 9 279 3 283 7 BL03-1-22 142 279 0.51 0.30964 0.00977 0.04344 0.00037 0.01335 0.00027 274 8 274 2 268 5 BL03-1-23 85 186 0.46 0.32663 0.01261 0.04346 0.00042 0.01430 0.00038 287 10 274 3 287 7 3.2 主量和微量元素
分析样品经显微镜下鉴定后，选择新鲜样品在无污染设备中进行加工粉碎，样品的粗碎和研磨工作在吉林大学地球科学院实验室完成。主量元素分析测试在中国科学院地质与地球物理研究所采用玻璃熔片大型X射线荧光光谱法(XRF)完成。微量元素在中国地质大学(武汉)地质过程与矿产资源国家重点实验室采用电感耦合等离子质谱(ICP-MS)法分析测试完成。对照国际标准参考物质BHVO-1(玄武岩)、BCR-2(玄武岩)和AGV-1(安山岩)的分析结果表明，主量元素分析精度和准确度优于5%，微量元素的分析精度和准确度一般优于10%。样品的主量和微量元素分析结果见表 2。
表 2 内蒙古阿鲁科尔沁旗安乐屯岩体主量、微量和稀土元素分析结果Table 2. Major, trace and rare earth elements compositions of the Anletun pluton in Ar Horqin Banner, Inner Mongolia样品 BL03-1 BL03-2 BL03-3 BL03-5 BL03-5r BL03-4 BL03-4r BL03-7 SiO2 71.00 70.61 70.92 70.41 54.70 54.41 TiO2 0.30 0.31 0.30 0.35 0.99 0.99 Al2O3 14.35 14.63 14.74 14.74 19.53 19.48 TFe2O3 3.42 3.49 3.56 3.67 6.70 6.67 MnO 0.04 0.04 0.06 0.04 0.13 0.13 MgO 0.86 0.94 1.01 0.92 3.39 3.37 CaO 2.36 2.95 2.96 2.42 7.10 7.06 Na2O 3.93 4.26 4.13 4.48 4.43 4.39 K2O 2.95 1.42 1.14 1.73 1.01 1.01 P2O5 0.07 0.08 0.08 0.06 0.31 0.31 烧失量 0.92 1.18 1.18 1.12 1.64 1.60 总计 100.20 99.91 100.08 99.94 99.93 99.43 Mg# 0.33 0.35 0.36 0.33 0.50 0.50 A/CNK 1.03 1.05 1.10 1.08 0.92 0.92 Na2O/K2O 1.33 3.00 3.62 2.59 4.39 4.35 Li 6.87 8.39 7.06 8.85 8.66 17.4 30.84 Be 0.82 0.82 0.83 0.86 0.85 1.13 2.26 Sc 7.10 5.52 7.89 2.23 2.25 10.7 16.05 V 40.56 45.31 45.17 43.54 44.15 124 160 Cr 1.71 1.85 2.13 1.89 1.83 47.06 85.8 Co 5.51 6.33 5.92 6.01 6.13 16.4 18.2 Ni 1.60 1.34 1.81 1.51 1.52 32.5 58.1 Cu 5.27 3.78 5.42 7.38 7.54 17.6 3.65 Zn 25.3 47.9 270.0 27.4 28.5 399 6486 Ga 12.92 13.66 12.87 14.78 14.71 22.77 22.60 Rb 60.56 50.12 44.95 42.34 41.96 45.09 151 Sr 391 497 421 684 695 1134 692 Y 7.59 1.78 10.84 1.41 1.54 13.44 14.16 Zr 90.9 95.1 97.1 106 117 93.1 105 Nb 2.78 1.56 3.40 1.83 1.86 4.70 4.55 Cs 1.87 2.12 2.68 1.84 1.85 3.14 3.29 Ba 1215 786 700 603 615 287 1048 La 11.14 9.70 13.9 9.63 9.70 15.9 20.03 Ce 17.99 14.32 23.16 13.28 13.53 37.87 46.07 Pr 2.13 1.50 2.80 1.38 1.37 4.95 5.92 Nd 7.52 4.84 10.0 4.26 4.39 21.2 24.14 Sm 1.40 0.59 1.87 0.52 0.51 4.31 4.66 Eu 0.72 0.71 0.66 0.80 0.80 1.38 1.25 Gd 1.42 0.46 1.85 0.38 0.40 3.73 3.84 Tb 0.22 0.05 0.30 0.05 0.04 0.53 0.53 Dy 1.36 0.27 1.87 0.22 0.23 2.84 2.82 Ho 0.29 0.06 0.40 0.04 0.04 0.53 0.54 Er 0.80 0.18 1.14 0.14 0.16 1.43 1.46 Tm 0.13 0.03 0.18 0.02 0.02 0.19 0.20 Yb 0.85 0.24 1.32 0.18 0.20 1.17 1.25 Lu 0.14 0.05 0.23 0.04 0.04 0.18 0.18 Hf 2.35 2.27 2.58 2.49 2.80 2.27 2.73 Ta 0.21 0.10 0.30 0.12 0.11 0.27 0.32 Pb 7.90 13.75 26.75 5.92 6.07 78.0 340 Th 2.45 0.81 4.32 1.30 1.32 2.19 3.82 U 0.49 0.30 0.83 0.36 0.37 0.57 1.04 Eu/Eu* 1.56 4.04 1.07 5.25 5.22 1.03 0.88 LREE 40.90 31.66 52.44 29.86 30.30 85.69 102 HREE 5.21 1.35 7.30 1.08 1.15 10.60 10.82 LREE/HREE 7.85 23.45 7.18 27.65 26.35 8.08 9.43 ∑REE 46.11 33.01 59.74 30.94 31.45 96.29 113 (La/Yb)N 9.40 28.45 7.54 38.08 34.90 9.78 11.47 注: Mg#= Mg2+/(Mg2++TFe2+); A/CNK = mole[Al2O3/(CaO+Na2O+K2O)]; Eu/Eu* =(Eu/0.0735)/[(Gd/0.259)+ (Sm/0.195)]/2;LREE = La+Ce+Pr+Nd+Sm+Eu; HREE = Gd+Tb+Dy+Ho+Er+Tm+Yb+Lu; ∑REE = LREE+HREE; (La/Yb)N =(La/0.687)/(Yb/0.493)；主量元素含量单位为%, 微量和稀土元素含量单位为10-6 3.3 锆石Hf同位素
锆石原位微区Hf同位素分析在中国地质科学院矿产资源研究所自然资源部成矿作用与资源评价重点实验室利用激光剥蚀多接收杯等离子体质谱(LA-MC-ICP-MS)完成。采用Neptune型多接受器等离子体质谱仪(MC-ICP-MS)，以氦气为载气，激光束斑直径为44 μm，激光脉冲速率和能量分别为6 Hz和100 mJ。详细仪器操作条件和分析方法可参照文献[37-38]。εHf(0)和εHf(t)值及模式年龄计算中，现今球粒陨石的176Lu/177Hf和176Hf/177Hf值分别采用0.0332和0.282772[39]，而亏损地幔的176Lu/177Hf和176Hf/177Hf值分别采用0.0384和0.28325[40]，二阶段模式年龄(TDM2)采用大陆地壳的fBCC(-0.65)进行计算[41]。分析数据的离线处理(样品信号区间的选择和同位素质量分馏校正)采用软件ICPMSDataCal完成[33]。样品中锆石的Hf同位素分析结果见表 3。
表 3 内蒙古阿鲁科尔沁旗安乐屯岩体锆石Hf同位素分析结果Table 3. Zircon Hf isotopic compositions of the Anletun pluton in Ar Horqin Banner, Inner Mongolia样品号 t/Ma 176Yb/177Hf 176Lu/177Hf 176Hf/177Hf 2σ εHf(0) εHf(t) 2σ TDM1/Ma TDM2/Ma fLu/Hf BL03-1-01 275 0.062493 0.001648 0.282485 0.000021 -10.2 -4.4 0.8 1104 1577 -0.95 BL03-1-02 275 0.061973 0.001681 0.282474 0.000019 -10.5 -4.8 0.7 1120 1601 -0.95 BL03-1-03 275 0.050284 0.001513 0.282497 0.000019 -9.7 -4.0 0.7 1082 1548 -0.95 BL03-1-04 275 0.049935 0.001621 0.282525 0.000020 -8.7 -3.0 0.7 1045 1485 -0.95 BL03-1-05 275 0.050878 0.001443 0.282556 0.000020 -7.7 -1.9 0.7 997 1415 -0.96 BL03-1-06 275 0.055263 0.001663 0.282471 0.000019 -10.6 -4.9 0.7 1123 1607 -0.95 BL03-1-07 275 0.055092 0.001847 0.282506 0.000018 -9.4 -3.7 0.6 1079 1530 -0.94 BL03-1-08 275 0.062363 0.001615 0.28245 0.000023 -11.4 -5.7 0.8 1153 1655 -0.95 BL03-1-09 275 0.064859 0.002239 0.282434 0.000022 -12.0 -6.3 0.8 1195 1697 -0.93 BL03-1-10 275 0.059236 0.001538 0.282517 0.000019 -9.0 -3.2 0.7 1054 1502 -0.95 BL03-1-11 275 0.05794 0.001578 0.282479 0.000021 -10.4 -4.6 0.8 1110 1589 -0.95 BL03-1-12 275 0.064227 0.001708 0.282469 0.000019 -10.7 -5.0 0.7 1129 1613 -0.95 BL03-1-13 275 0.048154 0.001763 0.282467 0.000019 -10.8 -5.1 0.7 1132 1617 -0.95 BL03-1-14 304 0.046523 0.001176 0.282487 0.000017 -10.1 -3.6 0.6 1087 1550 -0.96 BL03-1-15 275 0.077408 0.001962 0.282463 0.000021 -10.9 -5.2 0.8 1144 1629 -0.94 BL03-1-16 275 0.051141 0.001378 0.282452 0.000020 -11.3 -5.5 0.7 1143 1648 -0.96 BL03-1-17 275 0.044533 0.001812 0.282422 0.000024 -12.4 -6.7 0.9 1199 1720 -0.95 BL03-1-18 275 0.052647 0.001446 0.28241 0.000021 -12.8 -7.0 0.7 1203 1741 -0.96 BL03-1-19 275 0.053443 0.001465 0.282503 0.000020 -9.5 -3.7 0.7 1072 1533 -0.96 BL03-1-20 275 0.054552 0.001556 0.282461 0.000018 -11.0 -5.2 0.6 1135 1628 -0.95 BL03-1-21 275 0.043382 0.00141 0.282467 0.000018 -10.8 -5.0 0.6 1122 1614 -0.96 BL03-1-22 275 0.040575 0.001235 0.282481 0.000017 -10.3 -4.5 0.6 1097 1581 -0.96 BL03-1-23 275 0.042734 0.001302 0.282466 0.000017 -10.8 -5.0 0.6 1120 1614 -0.96 4. 分析结果
4.1 锆石U-Pb年龄
二长花岗岩样品BL03-1采自阿鲁科尔沁旗南部安乐屯东近山顶处(GPS坐标为北纬43°48′46.10″、东经120°08′30.10″)(图 1-b)。样品锆石晶体颗粒呈无色或浅褐色，柱状或六方双锥状，长宽比多介于2∶1~3∶1之间。阴极发光图像显示，大部分锆石呈自形晶或半自形晶，内部结构清晰，发育典型的岩浆生长环带，显示岩浆成因锆石的特征(图 2)。其较高的Th/U值(0.39~0.78；表 1)，也暗示岩浆成因[42-43]。锆石U-Pb年龄谐和图显示，所有分析数据均分布在谐和线上及其附近(图 3-a)。23个测点的206Pb/238U年龄值介于270~304 Ma之间(表 1)，除1个测点给出了304 ± 3 Ma的206Pb/238U年龄外，其余22个测点的年龄值相对集中(270 ± 3~279 ± 3 Ma)，其206Pb/238U年龄加权平均值为275 ± 1 Ma(MSWD = 0.57)(图 3-b)，代表了安乐屯岩体的形成时代。
![]()
图 2 安乐屯二长花岗岩部分锆石阴极发光(CL)图像(图中实线圆圈为锆石U-Pb测年位置，虚线圆圈为Hf同位素测定位置，括号内数字为εHf(t)值)Figure 2. CL images of some analyzed zircons from Anletun monzogranite![]()
图 3 安乐屯二长花岗岩锆石U-Pb年龄谐和图Figure 3. Concordia U-Pb diagram of zircons from Anletun monzogranite4.2 地球化学特征
4.2.1 主量元素
岩石地球化学分析测试结果表明，安乐屯岩体SiO2含量介于54.41%~71.00%之间，整体上属于中酸性侵入岩体，K2O含量为1.01%~2.95%，Na2O/K2O值为1.33~4.39，显示相对富钠、低钾的特征(表 2)。此外，安乐屯岩体TFe2O3为3.42%~6.70%、MgO为0.86%~3.39%、CaO为2.36%~7.10%，总体上Mg#值介于0.33~0.50之间。主量元素投图显示，安乐屯岩体总体上落入亚碱性系列和中钾钙碱性系列(图 4)。而铝饱和指数A/CNK≤1.1(0.92~1.10)，A/NK = 1.49~2.34，显示准铝质岩浆的特征(图 5)。哈克图解显示，安乐屯岩体主量元素与SiO2含量之间显示明显的规律性变化：除K2O与SiO2为正相关性外，其余主量元素TiO2、Al2O3、TFe2O3、MgO、CaO、Na2O和P2O5均随SiO2含量升高而降低(图 6)。安乐屯岩体主量元素与SiO2含量较好的相关关系暗示2个端元可能存在一定的成因联系。
![]()
图 4 安乐屯岩体TAS(a)和SiO2-K2O图解Figure 4. Plots of total alkalis versus SiO2(TAS)(a)and SiO2 versus K2O(b) for the Anletun pluton![]()
图 6 安乐屯岩体主量元素哈克图解Figure 6. Harker variation diagrams for major elements of the Anletun pluton4.2.2 微量元素
安乐屯岩体整体上稀土元素总量偏低(∑REE= 30.94×10-6~113×10-6)，轻、重稀土元素分馏较明显，LREE/HREE值为7.18~27.65，(La/Yb)N值介于7.54~38.08之间；石英二长闪长岩端元几乎不存在Eu异常(Eu/Eu* = 0.88~1.03)，而花岗岩端元以正Eu异常为特征，尤其样品BL03-2和BL03-5显示明显的正Eu异常(Eu/Eu* = 4.04~5.25)，存在重稀土元素分馏且显著亏损的特征(图 7-a)。在微量元素蛛网图(图 7-b)中，2种端元大体上显示类似的配分形式，均表现为富集Rb、Ba、Sr等大离子亲石元素(LILEs)，而相对亏损Nb、Ta、Ti等高场强元素(HFSEs)，类似于俯冲带岩浆岩的特征[46]；所不同的是，石英二长闪长岩端元的配分曲线整体较高，且具有较明显的负Ba异常，而花岗岩端元具有显著的Sm、Nd负异常，尤其样品BL03-2和BL03-5亏损的更明显。
![]()
4.3 锆石Hf同位素特征
安乐屯岩体二长花岗岩测年样品原位锆石Hf同位素分析测试结果显示，其锆石Hf同位素组成较均一，除1颗捕获锆石的176Hf/177Hf值为0.282487外，其余代表岩体形成年龄的22颗锆石的176Hf/177Hf值为0.282410~0.282556，以岩体形成年龄(275 Ma)计算得出，εHf(t)值均为负值，介于-7.0 ~ -1.9之间，投影到兴蒙造山带和燕山褶皱带交汇范围内及其附近区域，并与古—中元古代平均地壳演化范围一致(图 8)。岩体两阶段模式年龄TDM2为1415~1741 Ma，与中朝古板块北缘增生区内早二叠世岩浆岩给出的同位素模式年龄范围较一致[26-27]。
![]()
5. 讨论
5.1 安乐屯岩体的形成时代
安乐屯岩体最早被确定为燕山早期第一侵入期岩浆活动的产物，时代相当于早侏罗世①。前人根据研究区广泛发育燕山期岩浆活动，笼统地将安乐屯岩体划分为同一时期岩浆作用的产物，但该岩体无论从岩貌还是岩石类型都明显区别于研究区其他侏罗纪侵入岩。那么，安乐屯岩体的形成时代到底如何？由于一直缺乏精确的年代学资料的约束，本文对安乐屯岩体进行了锆石U-Pb年代学研究。定年结果表明，二长花岗岩样品(BL03-1)中22个测点给出了较一致且谐和的206Pb/238U年龄，其年龄加权平均值275 ± 1 Ma(MSWD = 0.57)应代表了该岩体的形成时代(图 3)，即早二叠世，而非前人认为的早侏罗世。近年来，在中朝古板块北缘达茂旗、镶黄旗、克什克腾旗、翁牛特旗等地区，以及西伯利亚古板块南缘苏左旗—锡林浩特—林西一带均有早二叠世岩浆事件的报道[20, 25-26, 31, 55]，与本文确定的安乐屯岩体形成年龄一致或相近，表明区域上存在早二叠世岩浆活动，进一步佐证了本文定年结果的可靠性。
5.2 岩石成因
从岩石学和地球化学特征看，安乐屯岩体由石英二长闪长岩-花岗闪长岩-二长花岗岩组成，出现暗色矿物黑云母和角闪石，SiO2含量为54.41% ~ 71.00%，Na2O含量为3.93% ~ 4.48%，Na2O/K2O为1.33 ~ 4.39，Mg#为0.33 ~ 0.50，A/CNK≤1.1，整体上具Eu的正异常特征，属于钙碱性岩系、准铝质岩石(图 4、图 5)，这些特征与I-型花岗岩相似。随着SiO2含量的增加，安乐屯岩体的其他主要氧化物呈现有规律的变化，其中，除K2O呈线性增加趋势外，其他氧化物TiO2，Al2O3，TFe2O3，MgO，CaO，Na2O和P2O5均随SiO2含量的增加而减少(图 6)，这表明岩浆上升侵位过程中发生了一定程度的分离结晶作用。由于在研究区及邻区没有发现同时代偏基性的岩石类型(SiO2＜54.41%)，笔者推测中性端元分析样品BL03-4和BL03-7的地球化学属性一定程度上应该代表了母岩浆成分特征，其较高的Nb/Ta(14~17)、Zr/Hf(39~41)值也表明并非岩浆高度分异的产物。基于此推断，酸性端元TFe2O3，MgO和CaO的减少可能与较早结晶的角闪石矿物的分离有关，而Ti和P的亏损一般是榍石和磷灰石等副矿物的分离结晶导致的，而就其稀土元素特征差异而言，可能是相关矿物分离结晶不均一性导致的。在长英质岩浆中Eu异常一般主要受控于斜长石，岩体整体上Eu的正异常特征可能与斜长石的堆晶作用有关，岩体较低的Rb/Sr值(0.04~0.22)也表明熔体可能包含较高程度的斜长石堆晶[56]。那么，显然石英二长闪长岩中Ba元素的相对亏损不大可能由斜长石的分离结晶导致，考虑到Ba替代K亦可赋存于黑云母中，所以其Ba元素的相对亏损可能与黑云母的分离结晶有关。花岗闪长岩样品BL03-2和BL03-5表现为除Eu元素外，中-重稀土(Sm~Lu)的极度亏损，研究表明中稀土元素含量主要受控于普通角闪石[57]，尤其在Dy和Er之间其具有最高的分配系数，所以普通角闪石的行为可以对REE(稀土元素)型式的形状起决定性的作用，而Eu元素随角闪石分离结晶应表现的亏损可能已经被斜长石的堆晶效应掩盖，所以普通角闪石的不均一性分离结晶作用一定程度上可以解释酸性端元的稀土元素特征。
安乐屯岩体具有低的Rb/Sr值、较低的初始176Hf/177Hf值(0.282410~0.282556)、负的εHf(t)值(-7.0 ~ -1.9)及老的模式年龄(1415~1741 Ma)，结合其I-型花岗岩成因特征及Nb、Ta亏损而Th相对富集的壳源属性，推断安乐屯岩体主要起源于古老下地壳物质的部分熔融。安乐屯岩体岩石类型包括石英二长闪长岩、花岗闪长岩和二长花岗岩，整体上相对富钠、贫钾，具有中钾钙碱性系列特征，属于准铝质系列岩石，加之富水矿物角闪石的出现，暗示其俯冲带成因特征(图 4、图 5)。岩体相对富集Rb、Ba、K等大离子亲石元素，亏损Nb、Ta、Ti等高场强元素的微量元素特征也反映了岩浆形成于与俯冲带有关的陆缘环境或岛弧构造环境(图 7-b)。此外，在Y-Nb和(Y+Nb)-Rb判别图解(图 9-a、b)中，岩体样品均投影到火山弧花岗岩区，同样暗示了类似的构造环境；而在Ta/Yb-Th/Yb构造环境判别图解(图 9-c)中，岩石样品均投影到了活动大陆边缘范围内[58]。综上推断，安乐屯岩体应是古亚洲洋俯冲消减背景下引起古老下地壳物质部分熔融的产物，并在岩浆上升侵位过程中发生了一定程度的角闪石、黑云母等矿物的分离结晶作用。
![]()
图 9 安乐屯岩体Y-Nb(a)、(Y+Nb)-Rb(b)和Ta/Yb-Th/Yb(c)相关图解MORB—大洋中脊玄武岩; WPB—板内玄武岩; WPVZ—板内火山岩带; ACM—活动大陆边缘; 0A—大洋岛弧Figure 9. Y-Nb(a),(Y+Nb)-Rb(b) and Ta/Yb-Th/Yb(c) diagrams of the Aletun pluton5.3 构造意义
研究区位于的内蒙古中部地区以发育多条蛇绿岩带和弧岩浆岩带为特点[59-62]，基于对这些关键性地质记录的认识不同，产生了针对该地区地质历史不同的构造演化模型。就研究区早二叠世的构造背景而言，有学者依据面状分布的早二叠世岩浆活动同时波及中朝古板块北缘和内蒙古中部造山带地区，并结合岩石组合特征推断该时期研究区处于块体碰撞之后的伸展环境[27]，也有学者认为古亚洲洋至少在晚泥盆世之前已经闭合，随后区域上发育广泛的裂谷活动，早二叠世岩浆岩应是大陆裂谷环境下的地质记录[63-68]。但本文的研究结果表明，研究区早二叠世安乐屯岩体构造判别图解显示其明显不同于后碰撞或裂谷环境花岗岩，而与活动大陆边缘弧花岗岩类似，暗示该时期古亚洲洋应该还未关闭。前人通过对晚古生代岩浆活动的研究也证实，研究区以南的中朝古板块北部存在晚石炭世—早二叠世连续的岩浆活动，形成了近东西向带状分布的岩浆弧，可以与现今典型的安第斯型活动大陆边缘弧相对比，其形成应与古亚洲洋板块向中朝古板块的俯冲有关[25-26, 51]，同样，前人通过对研究区以北苏左旗—西乌旗一线晚石炭世大陆边缘弧岩浆带的识别，结合该带以北地区贺根山蛇绿岩及早二叠世碱性花岗岩的存在，以及区内广泛分布的与俯冲体系相关的早二叠世晚期形成的大石寨组火山岩，认为它们共同构成了区域上晚古生代存在的沟-弧-盆体系，并推断其应与古亚洲洋板块的北向俯冲有关[28, 69-74]。综上可知，古亚洲洋最终关闭位置索伦-西拉木伦河缝合带的南北两侧在早二叠世仍均处于古洋壳俯冲的构造背景，两者都与安乐屯岩体成因背景相吻合。那么，安乐屯岩体的形成到底与哪个俯冲体系有关？这在一定程度上也制约着古亚洲洋最终闭合的位置。
前人大量的研究显示，中朝古板块北缘晚古生代岩浆岩的同位素组成与西伯利亚古板块南缘明显不同，中朝古板块北缘的岩浆岩一般具有负的εNd(t)值和εHf(t)值，暗示它们来自富集岩石圈地幔或古老地壳物质的熔融[24, 26, 47, 52-54, 75-79]，而西伯利亚古板块南缘的岩浆岩以正的εNd(t)值和εHf(t)值为特征，表明其起源于亏损地幔或新生地壳物质[14, 50, 70, 73]。基于此，通过对比缝合带两侧与俯冲相关的晚石炭世—早二叠世岩浆岩的锆石Hf同位素数据发现，安乐屯岩体锆石Hf同位素特征与中朝古板块北缘岩浆岩给出的εHf(t)值范围较一致，而明显不同于西伯利亚古板块南缘(图 8-b)。以此推断，安乐屯岩体的源区物质组成更亲缘于中朝古板块北缘的物质特性，其应与中朝古板块北缘同时代岩浆岩一样，形成于古亚洲洋板块向南俯冲于中朝古板块之下的构造环境。基于此，笔者推断安乐屯岩体很可能隶属于中朝古板块北缘古生代增生带，而古亚洲洋最终沿索伦-西拉木伦河关闭，向东延伸至研究区附近时，大洋的关闭位置很可能位于安乐屯岩体的北侧(图 1-a)。本文对索伦-西拉木伦河缝合带向东延伸位置提供了一定的约束，但确定板块缝合最终位置的关键地质记录的挖掘还需要进一步的工作。
6. 结论
(1) 内蒙古阿鲁科尔沁旗南部安乐屯岩体形成于275 ± 1 Ma，即早二叠世，而非前人认为的早侏罗世。
(2) 该岩体由石英二长闪长岩-花岗闪长岩-二长花岗岩组成，主要起源于古老下地壳物质的部分熔融，通过区域资料对比分析推断，岩体的形成与古亚洲洋板块向南俯冲于中朝古板块之下有关，指示其应隶属于中朝古板块北缘古生代增生带的一部分。
致谢: 中国地质调查局应用地质研究中心熊德清高级工程师、西华大学刘建康副教授、中国地质科学院探矿工艺研究所李洪梁工程师在研究工作中给予了帮助，西南石油大学李旭德硕士生、西南交通大学肖克锋硕士生协助完成了野外调查工作，审稿专家在文章修改过程中提出了宝贵意见，在此一并表示感谢。 -
![]()
图 1 西藏昌都地区地质地震简图(a，据参考文献[20]修改)及研究区滑坡分布图(b) (历史地震数据据参考文献[21-22])
KL—阿尼玛卿-昆仑结合带；GL—甘孜－理塘结合带；JS—金沙江结合带；AL—哀牢山结合带；NL—北澜沧江结合带；SL—南澜沧江结合带；LS—龙木错-双湖结合带；CM—昌宁-孟连结合带；BN—班公湖-怒江结合带；YL—印度河-雅鲁藏布江结合带
Figure 1. Geological and history earthquake distribution of Changdu (a) and landslide distribution of the study area (b)
![]()
图 2 察雅县城滑坡群滑坡后壁特征
Figure 2. Landslide cliff characteristics of landslides in Chaya county
![]()
图 3 察雅县城滑坡全貌及前缘变形
Figure 3. Panorama and deformation of the landslides in Chaya County
![]()
图 4 中铝新村南西侧3#滑坡滑坡体蠕滑变形
Figure 4. Creep deformation of the No.3 landslide on the southwest side of Zhongluxincun
![]()
图 7 古妥尔滑坡剖面
1—含粘土碎石土；2—页岩；3—断层破碎带；4—泥岩；5—粉砂岩；6—泥质粉砂岩；7—砂岩；8—滑面；9—逆冲断层；10—钻孔及编号; Q4del—第四系滑坡堆积；J1w—下侏罗统汪布组；T3d—上三叠统夺盖拉组
Figure 7. Profile of Gutuoer landslide
![]()
图 9 地震与非地震滑坡体积与地形阻止系数关系[18]及与察雅滑坡的对比
Figure 9. The volume and topographic blocking coefficient relation between seismic and non-seismic landslides, and the comparison with Chaya landslides
![]()
图 10 察雅县城滑坡区斜坡结构特征
1—河谷；2—顺向飘倾坡；3—顺向层面坡；4—顺向伏倾坡；5—斜顺向坡；6—横向坡；7—斜逆向坡；8—逆向坡；9—滑坡范围
Figure 10. Slope structure of landslides area in Chaya County
![]()
图 11 计算模型图
a—基岩；b—滑体；c—前部已变形滑体；d—最前部已变形滑体
Figure 11. Calculation model of landslide stability
表 1 察雅县城滑坡群各滑坡基本特征
Table 1 Basic characteristics of the landslides in Chaya County
编号 滑坡名称 滑动方向 最大滑移距离/m 前后缘最大高差/m 质心高度/m 长×宽/m 长宽比 堆积体平均厚度/m 体积/m3 平均坡度 后缘距分水岭距离/m 地形阻止系数 H1 古妥尔滑坡 0° 2200 880 490 2200×500 3.7 30 2300×104 9° 280 0.22 H2 司布日学滑坡 0° 2300 960 310 1200×650 1.8 35 2200×104 8° 310 0.23 H3 德日东滑坡 9° 2400 1040 550 2100×500 4.2 40 3360×104 25° 250 0.30 H4 德日西滑坡 4° 1600 820 430 1070×850 1.3 35 2400×104 31° 610 0.39 H5 察雅县中学南侧滑坡 28° 1340 760 370 1000×500 2.0 35 1130×104 30° 150 0.42 H6 中铝新村南东侧滑坡 55° 1170 680 380 630×690 0.9 20 600×104 28° 370 0.42 H7 中铝新村南西侧1#滑坡 42° 1070 560 320 900×325 2.8 40 880×104 29° 630 0.36 H8 中铝新村南西侧2#滑坡 22° 1200 660 390 850×375 2.3 30 620×104 23° 120 0.39 H9 中铝新村南西侧3#滑坡 26° 1300 720 330 1100×510 2.2 40 1900×104 29° 250 0.42 
下载: 导出CSV
表 2 研究区不同降雨频率下的年最大日降雨量估算结果
Table 2 Calculated results of maximum daily rainfall in a year at different frequencies
P/% KP H24P/mm 1 1.436 55.14 2 1.374 52.76 10 1.212 46.54
下载: 导出CSV
表 3 模型物理力学参数
Table 3 Physical and mechanical parameters of the model
材料 天然容重/(kN·m-3) C/kPa φ/° 后部滑坡体 26.2 9.8 25.6 中前部滑坡体 25.8 10.3 25.2 前部滑坡体 25.5 10.5 24.6 基岩层(泥质粉砂岩) 24 16.2 8.5
下载: 导出CSV
表 4 古妥尔滑坡堆积体稳定性计算结果
Table 4 Calculation results of stability for Gutour landslide
滑坡体编号 L1 L2 L3 L4 L5 工况一 天然状态 1.939 1.722 1.109 1.18 1.101 工况二 10年一遇降雨 1.933 1.592 1.113 1.147 1.077 工况三 50年一遇降雨 1.915 1.575 1.104 1.066 1.032 工况四 100年一遇降雨 1.903 1.487 1.099 1.059 1.006
下载: 导出CSV
-
王思敬. 地球内外动力耦合作用与重大地质灾害的成因初探[J]. 工程地质学报, 2002, 10(2): 115-117. doi: 10.3969/j.issn.1004-9665.2002.02.001 Scheidegger A E. Tectonic Predesign of Mass Movements, with Examples from the Chinese Himalaya[J]. Geomorphology, 1998, 26(1/2/3): 37-46. http://www.sciencedirect.com/science/article/pii/S0169555X98000506
Martel S J. Mechanics of Landslide Initiation as a Shear F racture Phenomenon[J]. Marine Geology, 2004, 203(3/4): 319-339. http://www.soest.hawaii.edu/martel/Martel.pubs.pdf/Martel_2004_Marine_Geo.pdf
居恢扬. 断裂构造对滑坡的控制意义[C]//滑坡文集: 第三集. 北京: 中国铁道出版社, 1982: 47-55. 吴玮江, 王念秦. 甘肃滑坡灾害[M]. 兰州: 兰州大学出版社, 2006: 82-87. Scheidegger A E, 艾南山. 武都地区的滑坡和泥石流[J]. 水土保持学报, 1987, 1(2): 19-27. https://www.cnki.com.cn/Article/CJFDTOTAL-TRQS198702002.htm 张永双, 苏生瑞, 吴树仁, 等. 强震区断裂活动与大型滑坡关系研究[J]. 岩石力学与工程学报, 2011, 30(增刊2): 3503-3513. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX2011S2020.htm 许强, 黄润秋. 5·12汶川大地震诱发大型崩滑灾害动力特征初探[J]. 工程地质学报, 2008, 16(6): 721-729. doi: 10.3969/j.issn.1004-9665.2008.06.001 黄润秋, 李为乐. 5·12汶川大地震触发地质灾害的发育分布规律研究[J]. 岩石力学与工程学报, 2008, 27(12): 2585-2592. doi: 10.3321/j.issn:1000-6915.2008.12.028 杨为民, 黄晓, 张春山, 等. 白龙江流域坪定-化马断裂带滑坡特征及其形成演化[J]. 吉林大学学报: 地球科学版, 2014, 44(2): 574-583. https://www.cnki.com.cn/Article/CJFDTOTAL-CCDZ201402017.htm 李晓, 李守定, 陈剑, 等. 地质灾害形成的内外动力耦合作用机制[J]. 岩石力学与工程学报, 2008, 27(9): 1792-1806. doi: 10.3321/j.issn:1000-6915.2008.09.006 吴树仁, 王涛, 石玲, 等. 2008汶川大地震极端滑坡事件初步研究[J]. 工程地质学报, 2010, 18(2): 145-159. doi: 10.3969/j.issn.1004-9665.2010.02.001 邹谨敞, 邵顺妹. 海原地震滑坡及其分布特征探讨[J]. 内陆地震, 1996, (1): 1-6. https://www.cnki.com.cn/Article/CJFDTOTAL-LLDZ601.000.htm 乔彦肖, 马中社, 吕凤军. 汶川地震地质灾害发育特点及动因机制分析[J]. 中国地质, 2009, 36(3): 736-741. doi: 10.3969/j.issn.1000-3657.2009.03.020 邓龙胜, 范文. 宁夏海原8.5级地震诱发黄土滑坡的变形破坏特征及发育机理[J]. 灾害学, 2013, 28(3): 30-37. doi: 10.3969/j.issn.1000-811X.2013.03.007 蒋瑶, 吴中海, 李家存, 等. 2010年玉树7.1级地震诱发滑坡特征及其地震地质意义[J]. 地质学报, 2014, 88(6): 1157-1176. https://www.cnki.com.cn/Article/CJFDTOTAL-DZXE201406016.htm 殷跃平. 汶川八级地震滑坡特征分析[J]. 工程地质学报, 2009, 17(1): 29-38. doi: 10.3969/j.issn.1004-9665.2009.01.004 樊晓一. 地震与非地震诱发滑坡的运动特征对比研究[J]. 岩土力学, 2010, 31(S2): 31-37. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX2010S2008.htm 潘桂棠, 李兴振, 王立全, 等. 青藏高原及邻区大地构造单元初步划分[J]. 地质通报, 2002, 21(11): 701-707. doi: 10.3969/j.issn.1671-2552.2002.11.002 安显银, 王启宇, 李勇, 等. 西藏昌都市察雅地区侏罗纪恐龙化石新发现[J]. 地质通报, 2021, 40(1): 189-193. http://dzhtb.cgs.cn/gbc/ch/reader/view_abstract.aspx?file_no=20210115&flag=1 秦四清, 李国梁, 薛雷, 等. 中国西南地区某些地震区未来震情研判[J]. 地球物理学进展, 2013, 28(5): 2407-2432. https://www.cnki.com.cn/Article/CJFDTOTAL-DQWJ201305022.htm 秦四清, 薛雷, 黄鑫, 等. 西藏地区未来强震预测[J]. 地球物理学进展, 2010, 25(6): 1879-1886. doi: 10.3969/j.issn.1004-2903.2010.06.001 钟康惠, 刘肇昌, 舒良树, 等. 澜沧江断裂带的新生代走滑运动学特点[J]. 地质论评, 2004, (1): 1-8. doi: 10.3321/j.issn:0371-5736.2004.01.001 郭长宝, 吴瑞安, 蒋良文, 等. 川藏铁路雅安-林芝段典型地质灾害与工程地质问题[J]. 现代地质, 2021, 35(1): 1-17. https://www.cnki.com.cn/Article/CJFDTOTAL-XDDZ202101002.htm 王腾, 孙晓光, 李白萍. 昌都市近36a暴雨气候特征分析[J]. 暴雨灾害, 2017, 36(1): 75-80. doi: 10.3969/j.issn.1004-9045.2017.01.010 张树轩, 杨为民, 程小杰, 等. 甘肃天水红旗山黄土滑坡群成因及稳定性分析[J]. 中国地质, 2017, 44(5): 924-937. https://www.cnki.com.cn/Article/CJFDTOTAL-DIZI201705008.htm 陈永明, 石玉成. 中国西北黄土地区地震滑坡基本特征[J]. 地震研究, 2006, 29(3): 276-280. doi: 10.3969/j.issn.1000-0666.2006.03.012 孟庆华, 孙炜锋, 张春山, 等. 陕西宝鸡地区胡家山滑坡风险性评价[J]. 地质通报, 2011, 30(7): 1155-1165. doi: 10.3969/j.issn.1671-2552.2011.07.018 -
期刊类型引用(11)
1. 刘娟娟，林一杰，张新英. 农田土水稻系统中硒与镉的相互关系研究. 农业与技术. 2025(01): 1-5 . 
百度学术
2. 谭卓贤，杜建军，孙星，易琼，徐培智，张木. 石灰、磷酸盐及硅酸盐对土壤硒有效性及水稻累积硒的影响. 江苏农业学报. 2024(03): 450-456 . 百度学术
3. 刘瑞，周卫军，彭素华，商贵铎，周雨舟，李敏. 湖南省石门县耕地土壤硒含量特征及影响因素解析. 土壤通报. 2023(04): 840-847 . 百度学术
4. Xiu-jin Liu，Ke Yang，Fei Guo，Shi-qi Tang，Ying-han Liu，Li Zhang，Hang-xin Cheng，Fei Liu. Effects and mechanism of igneous rock on selenium in the tropical soil-rice system in Hainan Province, South China. China Geology. 2022(01): 1-11 . 必应学术
5. 王蕊，陈楠，张二喜. 基于总量与形态的矿区周边土壤重金属生态风险与健康风险评估. 环境科学. 2022(03): 1546-1557 . 百度学术
6. 吴超，孙彬彬，成晓梦，周国华，贺灵，曾道明，梁倍源. 丘陵山区多目标区域地球化学调查不同成因表层土壤代表性研究——以浙江绍兴地区为例. 地质通报. 2022(09): 1539-1549 . 本站查看
7. 袁宏伟，陈江均，郭腾达，吴艳君，杨敏. 巴彦淖尔市临河区狼山镇和新华镇一带富硒土壤地球化学特征及影响因素. 地质与勘探. 2022(05): 1027-1041 . 百度学术
8. 刘冰权，沙珉，谢长瑜，周强强，魏星星，周梵. 江西赣县清溪地区土壤硒地球化学特征和水稻根系土硒生物有效性影响因素. 岩矿测试. 2021(05): 740-750 . 百度学术
9. 杨泽，刘国栋，戴慧敏，张一鹤，肖红叶，阴雨超. 黑龙江兴凯湖平原土壤硒地球化学特征及富硒土地开发潜力. 地质通报. 2021(10): 1773-1782 . 本站查看
10. 郭军，刘明，汤恒佳，廖琦，邢新丽. 湖南汨罗市范家园地区富硒土壤特征及其影响因素研究. 华南地质. 2021(04): 387-397 . 百度学术
11. 于炎炎. 安徽砀山县土壤锌地球化学特征及影响因素. 矿产与地质. 2021(06): 1147-1155+1170 . 百度学术
其他类型引用(2)
计量
- 文章访问数: 2592
- HTML全文浏览量: 579
- PDF下载量: 1560
- 被引次数: 13
