青藏高原东北部残积母质土壤发育过程研究<br/>——以青海湖北部宁夏剖面为例

青藏高原东北部残积母质土壤发育过程研究
——以青海湖北部宁夏剖面为例

2022-08-08谢丽倩鄂崇毅孙满平先巴吉

干旱区地理(汉文版) 2022年4期

谢丽倩，鄂崇毅,3，赵霞,3，李萍，张晶，孙满平，先巴吉

（1.青海师范大学青藏高原地表过程与生态保育教育部重点实验室，青海西宁 810001；2.青海师范大学地理科学学院/青海省自然地理与环境过程重点实验室，青海西宁 810008；3.青海省人民政府-北京师范大学高原科学与可持续发展研究院，青海西宁 810008）

土壤母质是土壤发育的基础，对土壤的性质与发生方向有深远影响。根据土壤母质与母岩的关系可分为残积母质和运积母质两种类型：残积母质是母岩在气候、生物等改造作用下原地风化后形成的松散碎屑物或疏松风化壳，母质成分继承了母岩的特点；运积母质是经重力、流水、风力等多种外力作用搬运而来的风化物质，通常其物质组成与母岩具有较大差异［1］。

青藏高原地势高耸，寒冷干燥，物理风化强烈，化学风化微弱，加之该区河流、冰川、风沙、湖泊等外力作用广泛，形成了大量的第四纪松散堆积物［2-4］，因此，该区土壤发育的物质基础以运积母质为主。随着土壤年代学工作的开展，近年来对高原风积母质土壤形成演变过程有了一定的认识［2-5］。然而，在野外考察中不乏残积母质发育土壤的存在，说明在一些特殊的气候时段，化学风化作用较强，原地风化也是高原土壤的重要成壤方式。对残积母质土壤剖面进行系统的年代学和土壤理化特征分析，可以有效补充高原土壤的形成演化过程，对区域土壤资源的保护与开发具有重要的指导意义。

青海湖地处青藏高原东北部，是中国最大的内陆湖，其流域是阻止柴达木沙漠东扩、共和沙地北扩的重要生态屏障，流域内的土壤安全更是维系生态环境安全的关键［6］。系统研究青海湖流域不同类型土壤的形成发育过程，对理解高原土壤的形成演化过程具有重要的参考价值。通过野外考察，在青海湖北部发现有较大面积的花岗岩残积母质土壤发育。因此，本文以青海湖北部的宁夏（NX）残积母质土壤剖面为研究对象，采用光释光（Optically stimulated luminescence，OSL）测年方法获取该剖面残积母质土壤发育的年代，基于土壤剖面理化特征，结合青海湖地区古气候环境演变记录，探讨青海湖流域土壤发育的过程和发育模式。

1 材料与方法

1.1 研究区概况

青海湖地区位于青藏高原东北部，年均气温为1.9 ℃，年均降水量为350～400 mm，降水集中在夏季，属典型的高原半干旱高寒气候［7］，土壤类型主要有高山寒漠土、高山草甸土、高山草原土、山地栗钙土、风沙土等［8］，地表植被以高原嵩草（Kobresia myosuroides）、芨芨草（Achnatherum splendens）、早熟禾属（Poaspp.）、狼毒（Stellera chamaejasme）、紫花针茅（Stipa purpurea）等为主［9］。

1.2 土壤剖面与样品采集

宁夏（NX）剖面（99°50′20″E，37°15′14″N，海拔3231.7 m）位于青海省海北州刚察县宁夏村附近（图1），剖面厚约40 cm，0～35 cm 为土壤较细颗粒，35～40 cm为花岗岩风化碎屑物，40 cm以下为花岗岩母岩，上覆植被为狼毒和芨芨草。本研究在剖面10 cm、20 cm、26 cm 和36 cm 处采集4 个OSL 测年样品，依次编号为NX1～NX4；以5 cm 土层厚度为间隔采集8个土壤散样，并在底部采集花岗岩母岩，编号为NX5，用于元素地球化学测试。

图1 采样点位置图Fig.1 Location of sampling point

1.3 光释光测年

OSL测年在青海师范大学自然地理与环境过程重点实验室光释光分室完成，测试仪器为Risø TL/OSL-DA-20-C/D型热/光释光仪，辐照源为人工β源90Sr/90Y。前处理过程参考E等［11］。全部样品等效剂量率(Dose，De)采用粗颗粒（63～90 μm）石英单片再生剂量法［12-13］进行测定，环境剂量率(Dose rate)中铀（U）、钍（Th）、钾（K）元素含量由电感耦合等离子体质谱法进行测定，再通过计算转化为β、γ剂量［14］，宇宙射线的贡献采用Prescott 等［15］提出方法计算。考虑到青海湖北部降水较少，蒸发量大，冻土发育，最终计算样品的含水量为（10±5）%，具体信息见表1。

1.4 土壤理化性质测试

土壤粒度实验在青海省自然地理与环境过程重点实验室粒度分室完成，使用仪器为英国马尔文公司制造的Mastersizer 2000 激光粒度仪（测量范围：0.02～2000 μm）。参考鹿化煜等［16］的实验法进行实验前处理：称取0.3～0.5 g 过20 目标准筛的土壤，在加热条件下分别加入10%的H2O2和10%的HCl去除土壤中的有机质和碳酸盐，静置后加入分散剂（NaPO3)6，最后对土壤样品进行上机测试。

土壤中的元素地球化学分析在西安地质调查中心测定，测试样品使用全岩样品，为去除土壤中碳酸盐的影响，另取样品NX1-1～NX1-8 进行化学淋滤处理：加10%的HCl 浸泡约30 min，静置2 h 后抽取上清液，将不溶于酸的部分经去离子水清洗3次，再放入60 ℃烘箱烘干，冷却后研磨至200 目以下。常量元素含量按照《GB/T14506.28-2010》标准使用荷兰帕纳科公司生产的Axios 型波长色散X 荧光光谱仪进行测试，微量元素按照《GB/T14506.30-2010》标准用美国热电公司生产的iCAP RQ 型电感耦合等离子体质谱仪进行测定。

土壤总有机碳（Total organic carbon，TOC）使用德国Elemertar 公司生产的vario TOC cube 总有机碳分析仪进行测试。测试样品经过950 ℃高温催化，燃烧产生CO2，利用IR（Infrared Radiation）检测CO2总量，再通过标准曲线计算得到TOC 的含量。TOC测试实验在青海省自然地理与环境过程重点实验室有机碳分室完成。

2 结果与分析

2.1 年代结果

样品NX2（63～90 μm）的衰退曲线（图2a）显示出石英信号在前2 s快速衰减到背景值，表明石英颗粒以快组分为主、晒退充分，即该样品适用于OSL测年；剂量生长曲线（图2b）中该样品的释光信号指数拟合较优。从NX 剖面4 个年代结果可以看出，NX剖面的土壤年代集中在早全新世，表层10 cm年代～0.4 ka，至下部20 cm土壤年代跳跃至～9 ka，表现为明显的沉积间断（表1）。

图2 NX2衰退及生长曲线Fig.2 Decay and growth curves of the sample from NX2

表1 NX剖面OSL测年结果及相关参数Tab.1 OSL dating results and related parameters of profile NX

2.2 土壤的物理特征

由测试结果（表2）可知，NX剖面土壤的组成特点为：粉砂＞砂＞黏粒，剖面在30～40 cm 与30 cm 以上土壤的机械组成差异明显，30～40 cm主要以砂为主，而30 cm以上土壤粉砂是主要的组成部分，黏粒含量也由下往上逐渐增加。粉黏比由下往上不断减小，可见剖面底部物质颗粒明显粗于上部土壤。

表2 NX剖面的物理性质Tab.2 Physical properties of profile NX

2.3 土壤的化学特征

从表3可看出，NX剖面主要以SiO2、Al2O3、Fe2O3和K2O 等氧化物为主，其中SiO2、Al2O3含量占比最高。与上地壳［17］（Upper continental crust，UCC）相比（图3），NX 剖面上部土壤的MgO、MnO、P2O5和TiO2含量与UCC 接近，Fe2O3和Na2O 含量低于UCC；CaO含量高于UCC；底部母岩和残积母质与上覆土壤的元素含量差异较大，上覆土壤和底部物质成分明显不同。

图3 NX剖面氧化物含量UCC标准化图Fig.3 UCC standardization diagram of oxide content in profile NX

表3 NX剖面的化学特征Tab.3 Chemical characteristics of of profile NX

沉积物的化学风化程度可以通过化学蚀变指数（Chemical index of alteration，CIA）表示［18］，表达式为：CIA=［Al2O3/(Al2O3+CaO*+Na2O+K2O)］×100。式中均为氧化物的摩尔分数，CaO*为硅酸盐矿物中的CaO［19］。为排除碳酸盐的影响，本研究通过加入10%的HCl 去除非硅酸盐矿物中的碳酸盐，对CIA值进行校正。NX 剖面中CIA 平均值为58.58，最高为61.88，底部未经风化母岩的CIA 本底值为48.34。Rb/Sr 比值在黄土高原被用作衡量风化成壤程度的指标，比值越高指示剖面化学风化强度越强［20］。NX剖面Rb/Sr 比值在35 cm 以上变化较小，介于0.34～0.55，从下往上不断增大；底部母岩的Rb/Sr 比值为4.22，残积母质为0.69，与上部土壤有明显区别。

NX 剖面TOC 含量为0.06%～2.14%（平均值为0.66%），剖面表层与底部物质的有机碳含量差异较大。样品NX1-8 接近基岩，颗粒大，有机碳含量最低，表层（0～15 cm）受到植被影响，有机碳含量大。

3 讨论

3.1 青海湖北部NX剖面物源分析

沉积物发育的土壤序列成分取决于物质来源的组成成分和风化程度。土壤中稳定元素在风化、搬运、沉积过程中能保留母岩的特征，可用于研究沉积物的物质来源［21］。在化学风化过程中，由于K常吸附在粘土矿物中不易发生淋溶，K2O/Al2O3比值受粒度影响较小［22］。而Ti、Al 的溶解度小，在化学分化过程中淋溶作用较弱，TiO2/Al2O3比值可以反映风成沉积的物质贡献，可作为判别物源变化的指标［23-24］。Zr、Nb 是赋存于碎屑岩中的稳定元素，迁移率低，可指示母岩的特征，由于其稳定性特征常用来指示物源［25-26］。

基于K2O/Al2O3和Zr/Nb比值、TiO2/Al2O3和Zr/Nb比值的图解显示：NX 剖面的母岩和残积母质与上部土壤明显不同，表明残积母质土壤剖面上部35 cm 以上土壤与残积母质和下覆母岩的物质来源不同。通过与青海湖地区的近源沉积物（河流沉积和湖相沉积［10］）对比（图4），NX 剖面具有相对较高的K2O/Al2O3比值和TiO2/Al2O3比值和较低的Zr/Nb 比值；NX 剖面与青海湖地区黄土［10］相比，上部土壤与青海湖地区黄土有相似的范围：K2O/Al2O3比值在0.2～0.3之间，TiO2/Al2O3比值在0.05～0.07之间，说明上部35 cm 土壤母质主要是风尘加积形成。下部35～40 cm 风化碎屑物的K2O/Al2O3比值和Zr/Nb 比值、TiO2/Al2O3比值和Zr/Nb 比值接近下部花岗岩母岩，说明对母岩的化学风化作用主要发生在此段。

图4 青海湖地区样品的K2O/Al2O3-Zr/Nb、TiO2/Al2O3-Zr/Nb比值Fig.4 K2O/Al2O3 vs.Zr/Nb,TiO2/Al2O3 vs.Zr/Nb ratios of samples in Qinghai Lake area

3.2 青海湖北部NX剖面土壤的发育程度

土壤的化学风化程度是评价土壤发育程度的重要标准之一，CIA 值作为衡量化学风化程度的重要指标，可表征土壤的发育程度。当CIA 值介于50～65之间，反映寒冷、干燥的气候条件下较低的化学风化程度［27］。NX 剖面的CIA 平均值为58.58，基本处于弱化学风化阶段，符合青藏高原东北部寒冷干燥的气候特征。Rb 主要存在于云母、粘土矿物等难风化的矿物中，Sr主要存在于易风化的碳酸盐矿物中，在风化成壤过程中，Sr 相对于Rb 更容易流失［28-29］。因此Rb/Sr 比值越大，风化和淋溶作用越强［20］。NX 剖面的Rb/Sr 比值在残积母质层以上变化较小（0.34～0.55），底部由于母岩地球化学性质的影响，与上部风尘物质的Rb/Sr 比值差异较大。母岩的CIA 为48.34、Rb/Sr 比值为4.22，代表了未风化状态的本底值，随着化学风化程度加强，CIA 值增大。此外，由于上部风尘沉积物受淋溶和风化作用影响，Sr迁移至残积母质层［30］，使得残积母质的Rb/Sr 比值为0.69。残积母质以上的风尘沉积物质的CIA和Rb/Sr比值变化较小，呈现从下往上增大的趋势，表明土壤剖面上部的风化程度大于下部。在图5 中，TOC 和粉黏比呈相反趋势，表层（0～15 cm）的TOC 含量较高，受到植被和人类活动的改造作用，碳输入量增加，对土壤的改造作用强。而表层的粉黏比最小，底部最大，说明表层土壤在植物根系的改造作用下颗粒变细，土壤的发育程度高于底部土壤。

图5 年代、CIA、Rb/Sr、TOC、粉黏比随深度变化特征图Fig.5 Characteristic diagram of age,CIA,Rb/Sr,TOC and silt/clay ratio components with depth

3.3 成土母质发育时代及气候背景

前人对青海湖流域的风成砂、黄土、青海湖钻孔以及湖面变化等古环境记录体的研究，为该区气候演变提供了大量的证据，从而为了解流域内土壤发育的气候背景提供支撑［6,31-33］。这些古环境记录表明早全新世是青海湖流域气候干暖，风沙活动强烈的时期。

NX 剖面底部残积母质发育的年代应不早于10.02±0.79 ka，更早的末次盛冰期由于气候过于寒冷，化学风化非常微弱，末次冰消期16.0～11.5 ka 期间气候虽然略微转好，但在已发表的风成沉积记录中古土壤发育在此时期的记录较少［6,34］，因此NX残积母质发育的年代应在早全新世期间。这说明尽管早全新世整体风沙活动强烈，但气候较为温暖［29,33,35-36］，在适当的降水条件下，其化学风化作用也较强，裸露的花岗岩可能在地衣苔藓等先锋植物和气候作用下发生化学风化过程［37-38］并形成薄层的残积母质；后期随着残积母质上的植被演替，草本植被迁入，风尘物质被快速捕获和堆积，此时气候与植被直接作用于风尘母质，青藏高原干冷的气候和低矮的草本植被对下覆母岩的风化作用显著减弱甚至停滞。考虑青藏高原整体的气候以冷干为主，变化并不显著，此种早期残积母质发育与后期风尘母质加积的混合母质发育过程在整个全新世和现代时段都应在进行，但需更多剖面进行验证。

35 cm 以上风尘母质的年代在9.71～8.67 ka 之间，与前人在青藏高原东北部土壤剖面的风沙活动主要集中在末次冰消期和早全新世时期的结论相呼应［6,34,39］。这表明在早全新世干暖的气候背景下，青海湖湖面下降，湖滩裸露［31］，为流域内风沙活动提供更多的风尘物质，为区域土壤的发育提供丰富的风尘母质。因此NX剖面风尘母质的发育时间集中在早全新世可能与该区丰富的风尘供应有关。NX 剖面在20 cm 处8.67±0.62 ka 跳跃至10 cm 处的0.41±0.04 ka，出现～9 ka 的沉积间断。该间断的形成一方面与气候干湿状况相关，据青海湖湖面记录，湖面在7.5～2.5 ka 呈上升趋势［31］，这指示了剖面在间断的9 ka 中气候相对暖湿，限制了风尘物质的释放，风尘沉积速率较低［40］；另一方面，沉积间断也可能与剖面所处位置相关，由于地形风强盛，不利于风尘保存。前人基于高密度OSL 方法揭示了青藏高原东北部风尘沉积物广泛存在沉积间断［2,6,39,41］。因此，本剖面的间断并非个例。