青藏高原东北部土壤痕量元素的环境地球化学行为解析：分布、来源与生态响应

青藏高原东北部土壤痕量元素的环境地球化学行为解析：分布、来源与生态响应一、引言1.1研究背景与意义青藏高原，被誉为“世界屋脊”与“亚洲水塔”，作为全球海拔最高且面积最大的高原，在全球气候调节、水资源涵养以及生物多样性维护等方面发挥着举足轻重的作用。其东北部地区，地处亚洲季风区西北缘，受季风环流与西风环流的共同影响，生态环境极为敏感且脆弱。土壤，作为陆地生态系统的关键组成部分，不仅是植物生长的物质基础，更是元素循环与能量转换的重要场所。土壤中的痕量元素，尽管含量甚微，却在生态系统的物质循环、生物地球化学过程以及生态平衡维持等方面扮演着不可或缺的角色。痕量元素，通常是指在土壤和植物中少量存在、曾经一度难以测定的元素，其生理功能与其浓度紧密相关。在低浓度时，部分痕量元素可能是生物生存的必需元素，或是被生物被动吸收；然而，当浓度过高时，某些痕量元素则可能对生物产生毒害作用，进而影响生态系统的结构与功能。例如，锌（Zn）、铜（Cu）、铁（Fe）等痕量元素是植物生长发育所必需的微量元素，它们参与植物的光合作用、呼吸作用以及酶的合成等重要生理过程。然而，当土壤中这些元素的含量过高时，可能会对植物产生毒害，抑制植物的生长，甚至导致植物死亡。此外，镉（Cd）、汞（Hg）、铅（Pb）等痕量元素属于重金属元素，具有较强的毒性，它们在土壤中的积累不仅会影响土壤微生物的活性和群落结构，还可能通过食物链的传递，对人类健康构成严重威胁。青藏高原东北部地区的土壤痕量元素研究，对于深入理解区域生态系统的结构与功能、生物地球化学循环过程以及生态环境演变规律具有重要意义。一方面，该地区特殊的地理环境和气候条件，使得土壤痕量元素的来源、迁移转化过程以及分布特征具有独特性。研究这些独特性，有助于揭示土壤痕量元素在极端环境条件下的环境地球化学行为，丰富和完善土壤地球化学理论。另一方面，随着全球气候变化和人类活动的加剧，青藏高原东北部地区的生态环境面临着严峻的挑战，如土壤侵蚀、土地退化、生物多样性减少等。土壤痕量元素作为生态环境变化的敏感指示因子，其含量和分布的变化能够直观地反映出生态环境的演变趋势。通过研究土壤痕量元素，可以及时准确地监测生态环境的变化，为制定科学合理的生态环境保护政策和措施提供重要的理论依据。此外，青藏高原东北部地区是我国重要的畜牧业生产基地和生态屏障，土壤痕量元素的研究对于保障区域农牧业的可持续发展、维护生态平衡以及促进经济社会的协调发展具有重要的现实意义。土壤痕量元素的含量和分布直接影响着农作物和牧草的生长发育、品质和产量，进而影响着农牧业的生产效益。同时，土壤痕量元素的污染还可能导致农产品和畜产品的质量安全问题，威胁人类的健康。因此，深入研究土壤痕量元素，对于保障区域农牧业的可持续发展、提高农产品和畜产品的质量安全水平具有重要的现实意义。1.2研究现状综述近年来，国内外学者对青藏高原土壤元素展开了多方面研究，在土壤常量元素地球化学特征、土壤微量元素分布特征等方面取得了一定进展。在常量元素研究方面，有学者对青藏高原东北部风成砂-古土壤序列泽库剖面常量元素氧化物含量及其比值进行分析，结合14C测年，揭示出剖面沉积物化学组成以SiO₂、Al₂O₃、CaO为主，其中SiO₂和CaO标准差较高，对气候变化较为敏感。自剖面底部向上，SiO₂含量和残积系数呈现先减少后增加的趋势，高值指示冷干的气候意义，CaO含量和退碱系数则相反，高值指示温湿的气候意义。研究表明，9.4kaBP以来青藏高原东北部的气候经历了由温湿向冷干的转变过程，可划分为9.4-4.2kaBP气候温湿期和4.2kaBP至今气候冷干期，其中在6.6-6.2kaBP、2.4-2.0kaBP和1.7-1.5kaBP存在次一级的气候波动。在微量元素及痕量元素研究领域，部分研究关注到青藏高原土壤中微量元素的分布特征。有研究采集青藏高原地区不同海拔和地形位置的表层土壤样品，对土壤样品的微量元素含量进行测定后发现，随着海拔的升高，土壤中部分微量元素含量呈现不同程度的下降趋势，且不同地形位置土壤微量元素特征存在差异，例如山地土壤中某些微量元素含量高于平原土壤。还有研究聚焦于土壤中重金属等痕量元素，分析其来源、迁移转化等过程，发现人类活动如矿业开采、农业活动等对局部地区土壤痕量元素含量有一定影响。然而，当前研究仍存在一些不足。一方面，针对青藏高原东北部土壤痕量元素的研究，在空间覆盖面上还不够广泛，多数研究集中在部分典型区域，对于一些偏远或生态环境独特的小区域关注较少，这使得难以全面准确地掌握整个青藏高原东北部土壤痕量元素的分布格局。另一方面，在土壤痕量元素的迁移转化机制研究方面，虽然已有一定探索，但受到该地区复杂的自然环境条件（如高海拔、低温、强辐射等）以及多变的气候因素影响，现有的研究成果还无法清晰完整地解释痕量元素在土壤-植物-大气等多界面间的迁移转化规律。同时，在研究痕量元素对生态系统的影响时，多集中于对植物生长发育的影响，而对于其在土壤微生物群落结构与功能、土壤酶活性等方面的影响研究相对薄弱，缺乏系统全面的生态效应评估体系。此外，在全球气候变化和人类活动日益加剧的背景下，针对未来青藏高原东北部土壤痕量元素动态变化的预测研究也较为匮乏，难以满足生态环境保护和可持续发展的需求。1.3研究内容与方法1.3.1研究内容土壤痕量元素的含量与分布特征：系统分析青藏高原东北部不同区域、不同土壤类型以及不同地形地貌条件下土壤中多种痕量元素（如Cd、Hg、Pb、Zn、Cu、As等）的含量水平，绘制详细的含量分布图，明确各痕量元素在该地区的空间分布格局，探究其分布是否存在明显的区域差异以及与地形、土壤类型等因素的相关性。例如，分析在山地、河谷、平原等不同地形下痕量元素含量的变化规律，研究不同土壤类型（如高山草甸土、草原土、风沙土等）中痕量元素的富集或贫化特征。土壤痕量元素的赋存形态与迁移转化规律：运用化学连续提取法等技术手段，研究土壤痕量元素的赋存形态，包括可交换态、碳酸盐结合态、铁锰氧化物结合态、有机结合态和残渣态等，分析不同形态痕量元素在土壤中的相对含量及分布特征。通过室内模拟实验和野外原位监测，结合相关分析、主成分分析等方法，深入探究土壤痕量元素在不同环境条件下（如温度、湿度、酸碱度变化等）的迁移转化规律，明确影响其迁移转化的关键因素。例如，研究在季节性冻融过程中，土壤痕量元素形态的变化以及对其迁移能力的影响；分析降水强度和频率改变时，痕量元素在土壤-水体界面的迁移过程。土壤痕量元素的来源解析：综合运用多元统计分析（如主成分分析、聚类分析等）、同位素示踪技术（如铅同位素、硫同位素等）以及相关的源解析模型（如正定矩阵因子分解模型等），对青藏高原东北部土壤痕量元素的来源进行全面解析。明确自然源（如岩石风化、大气沉降等）和人为源（如矿业活动、农业生产、交通排放等）对土壤痕量元素的相对贡献，确定主要的污染源及其贡献比例。例如，通过铅同位素分析，判断土壤中铅元素是主要来源于本地岩石风化还是远距离传输的大气沉降；利用正定矩阵因子分解模型，量化矿业开采、农业化肥使用等人为活动对土壤中镉、汞等痕量元素的贡献程度。土壤痕量元素的生态效应评估：选取该地区典型的植物物种和土壤微生物群落，研究土壤痕量元素对植物生长发育（包括种子萌发、幼苗生长、生物量积累、光合作用等指标）、生理生化特性（如抗氧化酶活性、渗透调节物质含量等）以及土壤微生物群落结构与功能（如微生物数量、多样性、酶活性等）的影响。建立土壤痕量元素含量与生态效应指标之间的定量关系，评估土壤痕量元素对生态系统的潜在风险，为制定合理的生态环境保护策略提供科学依据。例如，通过盆栽实验，研究不同浓度的镉、铅等痕量元素对青稞、燕麦等高原农作物生长和品质的影响；利用高通量测序技术分析土壤微生物群落结构在痕量元素胁迫下的变化，探讨痕量元素对土壤生态功能的影响机制。1.3.2研究方法样品采集：在青藏高原东北部，依据不同的地理分区、土壤类型和土地利用方式，采用网格布点法与典型样点相结合的方式进行土壤样品采集。每个采样点设置3-5个重复样，使用不锈钢土钻采集0-20cm的表层土壤样品，剔除土壤中的植物根系、石块等杂物，将样品装入干净的聚乙烯塑料袋中，标记好采样地点、时间、经纬度等信息。对于深层土壤样品，按照一定的深度间隔（如20-40cm、40-60cm等）进行分层采集。同时，在采样区域同步记录地形地貌、植被类型、气候条件等环境信息，以便后续分析土壤痕量元素与环境因素的关系。样品分析：将采集的土壤样品自然风干后，过2mm尼龙筛，去除残留的植物根系和砂砾等杂质。对于需要测定全量痕量元素的样品，采用盐酸-硝酸-氢氟酸-高氯酸（HCl-HNO₃-HF-HClO₄）四酸消解体系，利用电热板进行消解，使土壤中的痕量元素完全释放到溶液中。消解后的溶液使用电感耦合等离子体质谱仪（ICP-MS）测定Cd、Hg、Pb、Zn、Cu、As等痕量元素的含量。对于测定痕量元素赋存形态的样品，采用Tessier连续提取法或BCR三步提取法，将土壤中的痕量元素依次提取为不同的化学形态，然后使用原子吸收光谱仪（AAS）或电感耦合等离子体发射光谱仪（ICP-AES）测定各形态中痕量元素的含量。在分析过程中，采用国家标准物质（如GBW07401-GBW07408等）进行质量控制，确保分析结果的准确性和可靠性，每批次样品分析时均插入一定比例的空白样和重复样，以监控分析过程中的误差。数据处理与分析：运用Excel软件对原始数据进行整理和初步统计分析，计算各痕量元素含量的平均值、标准差、最小值、最大值等统计参数。采用SPSS、Origin等数据分析软件进行相关性分析、主成分分析、聚类分析等多元统计分析。相关性分析用于探究土壤痕量元素之间以及痕量元素与环境因子（如土壤pH值、有机质含量、地形海拔等）之间的相互关系；主成分分析和聚类分析用于对土壤痕量元素的来源进行初步解析，识别影响土壤痕量元素分布的主要因素和潜在污染源。利用Isoplot等软件进行同位素数据分析，通过铅同位素、硫同位素等比值来进一步确定土壤痕量元素的来源。采用风险评估模型（如潜在生态风险指数法、污染负荷指数法等）对土壤痕量元素的生态风险进行评估，根据评估结果对研究区域土壤痕量元素的污染状况进行分级和评价，为区域土壤环境保护和管理提供科学依据。二、研究区概况2.1地理位置与地质背景青藏高原东北部地处欧亚大陆内部，介于东经90°-104°、北纬32°-40°之间，涵盖了青海的大部分地区以及甘肃、宁夏、陕西的部分区域。该区域北起祁连山，南至昆仑山，西临柴达木盆地，东接黄土高原，是青藏高原向内陆过渡的关键地带，地理位置极为特殊。其平均海拔在3000-5000米之间，地势西北高、东南低，呈现出高山、高原、盆地、河谷等复杂多样的地形地貌。在地质构造上，青藏高原东北部位于印度板块与欧亚板块碰撞挤压的前缘地带，是多个构造单元的汇聚区域，地质构造活动极为频繁且复杂。该区域主要涉及祁连山造山带、柴达木地块、东昆仑造山带等重要的地质构造单元。祁连山造山带是在早古生代时期，由于华北板块与柴达木板块的碰撞拼合而形成，经历了多期次的构造变形和岩浆活动，其内部岩石类型丰富多样，包括花岗岩、片麻岩、砂岩、页岩等，这些岩石的风化产物为土壤的形成提供了重要的物质基础。柴达木地块作为一个相对稳定的刚性地块，其基底主要由前寒武纪变质岩组成，在漫长的地质历史时期中，受到周边构造单元的挤压和改造，地块边缘形成了一系列的断裂构造和褶皱构造。东昆仑造山带则是在中生代时期，由于羌塘板块与欧亚板块的碰撞而形成，经历了强烈的构造隆升和岩浆活动，其岩石类型以火山岩、侵入岩和变质岩为主。区内主要山脉包括祁连山、昆仑山、阿尼玛卿山等，这些山脉均呈西北-东南走向，是青藏高原东北部的重要地理屏障。祁连山作为青藏高原东北部的重要山脉，其山体高大雄伟，平均海拔在4000米以上，主峰团结峰海拔5827米。祁连山的山体主要由花岗岩、片麻岩等岩石组成，这些岩石经过长期的风化侵蚀作用，形成了丰富的碎屑物质，为周边地区土壤的形成提供了大量的成土母质。昆仑山是青藏高原的重要山脉之一，其平均海拔在5500米以上，是亚洲大陆南部和中部地区主要山脉的主干。昆仑山的岩石类型主要为花岗岩、闪长岩等，其风化产物在重力、风力和水力等作用下，向周边地区搬运和堆积，对土壤的形成和发育产生了重要影响。阿尼玛卿山是昆仑山系的东支，平均海拔在4000-5000米之间，主峰玛卿岗日海拔6282米。阿尼玛卿山的岩石主要由石灰岩、砂岩等组成，其风化物质在当地的气候和地形条件下，参与了土壤的形成过程。主要河流有黄河、湟水、大通河等，这些河流对土壤的形成和发育起着至关重要的作用。黄河作为我国的第二长河，流经青藏高原东北部，其携带的大量泥沙和矿物质在河流两岸沉积，形成了肥沃的冲积土壤。湟水是黄河上游的重要支流，主要流经青海东部地区，湟水流域地势相对平坦，水流较为平缓，其携带的泥沙和养分在河谷地区沉积，为当地土壤提供了丰富的物质来源，促进了土壤的发育和肥力的提高。大通河也是湟水的重要支流，其流域内植被覆盖较好，河水对土壤的侵蚀作用相对较弱，同时河水的灌溉作用为土壤补充了水分，有利于土壤中微生物的活动和养分的转化，对土壤的理化性质产生了积极影响。这些山脉和河流通过影响地形地貌、气候条件、物质迁移等，深刻地塑造了该地区土壤的形成过程和分布特征。2.2气候与植被特征青藏高原东北部属于典型的高原大陆性气候，具有气温低、昼夜温差大、降雨少而集中、日照长、太阳辐射强等特点。该区域年平均气温在-5.7℃-8.5℃之间，其中1月平均气温在-17℃--5℃之间，7月平均气温在5.3℃-20℃之间。冬季漫长而严寒，夏季短促且凉爽，部分高海拔地区无绝对无霜期。年降水量在100-600毫米之间，总体呈现出从东南向西北递减的趋势。例如，东南部的玛曲和若尔盖等地，年降水量可达600毫米以上，且降水主要集中在7-8月；而西北部的柴达木盆地，年降水量则不足50毫米，气候极为干旱。降水的年内分配不均，多集中在夏季，且以夜雨为主，夜雨率较高。这种降水特征与地形和大气环流密切相关，夏季来自印度洋和太平洋的暖湿气流受地形阻挡，在该地区形成降水，夜间山坡冷却快，冷空气下沉，谷底暖湿空气被迫抬升，形成降雨。该区域植被类型丰富多样，主要包括高山草甸、高山草原、荒漠草原、灌丛等。高山草甸主要分布在海拔较高、气候较为湿润的地区，如祁连山、昆仑山等山脉的高海拔地带，植被以嵩草、苔草等草本植物为主，群落结构简单，覆盖度较高，一般在70%-90%之间。高山草原主要分布在海拔相对较低、气候较为干旱的地区，如柴达木盆地周边的高原地区，植被以针茅、羊茅等草本植物为主，覆盖度相对较低，在30%-50%之间。荒漠草原则分布在更为干旱的地区，如柴达木盆地内部，植被稀疏，主要由旱生的小半灌木和草本植物组成，覆盖度在10%-30%之间。灌丛主要分布在山地的阴坡或河谷地带，常见的灌丛植物有金露梅、银露梅、沙棘等。植被对土壤痕量元素具有重要作用。一方面，植被通过根系吸收土壤中的痕量元素，影响痕量元素在土壤-植物系统中的迁移和转化。不同植被类型对痕量元素的吸收能力和选择性不同，例如，高山草甸植被根系发达，对锌、铜等痕量元素具有较强的吸收能力，能够将土壤中的痕量元素富集到植物体内。另一方面，植被的枯枝落叶分解后形成的腐殖质，能够与土壤中的痕量元素发生络合、吸附等作用，影响痕量元素的赋存形态和生物有效性。腐殖质中的有机官能团能够与痕量元素形成稳定的络合物，降低痕量元素的迁移性和生物可利用性。此外，植被还能够通过影响土壤的物理性质（如土壤结构、孔隙度等）和化学性质（如土壤pH值、有机质含量等），间接影响土壤痕量元素的环境地球化学行为。例如，植被覆盖度高的地区，土壤侵蚀作用较弱，有利于保持土壤中痕量元素的含量和分布稳定性。2.3土壤类型与分布青藏高原东北部地域辽阔，自然环境复杂多样，土壤类型丰富，主要包括高山草甸土、高山草原土、风沙土、沼泽土、灰褐土等。高山草甸土主要分布在海拔较高、气候较为湿润、植被以高山草甸为主的区域，如祁连山、昆仑山等山脉的高海拔地带。其成土过程主要是在寒冷湿润的气候条件下，草甸植被残体的腐殖质化过程和矿物质的风化过程相互作用。高山草甸土的剖面特征明显，表层为深厚的草毡层，根系交织紧密，富含腐殖质，颜色较深，一般呈暗棕色或黑色；其下为腐殖质层，厚度适中，颜色稍浅；再往下是过渡层和母质层。高山草甸土的质地多为壤质土，土壤结构良好，通气性和保水性适中。土壤中有机质含量较高，一般在5%-15%之间，这是由于草甸植被生长茂盛，每年有大量的枯枝落叶归还到土壤中，经过微生物的分解和转化形成了丰富的腐殖质。土壤pH值呈中性至微酸性，一般在6.0-7.0之间，这种酸碱环境有利于土壤中养分的释放和植物对养分的吸收。高山草甸土中氮、磷、钾等养分含量相对较高，尤其是全氮含量，通常在0.3%-0.8%之间，这为高山草甸植被的生长提供了良好的土壤条件。高山草原土主要分布在海拔相对较低、气候较为干旱、植被以高山草原为主的地区，如柴达木盆地周边的高原地区。其成土过程主要受干旱气候和草原植被的影响，土壤腐殖质积累较少，钙化过程较为明显。高山草原土的表层有一定厚度的草根层，但相较于高山草甸土，草根层较薄，根系分布相对稀疏。腐殖质层较薄，颜色较浅，多为淡棕色。土壤质地多为砂质壤土或壤土，通气性较好，但保水性较差。由于气候干旱，土壤中有机质含量较低，一般在1%-5%之间。土壤pH值呈碱性，一般在7.5-8.5之间，这是因为土壤中含有较多的碳酸钙等碱性物质。土壤中氮、磷等养分含量相对较低，全氮含量一般在0.1%-0.3%之间，速效磷含量也较低，这在一定程度上限制了草原植被的生长和生产力。风沙土主要分布在沙漠边缘、河谷沙滩以及风力作用较强的地区，如柴达木盆地内部、黄河谷地部分地段等。风沙土是在风力作用下，沙质物质堆积形成的土壤类型，其成土过程微弱，土壤发育程度低。风沙土的剖面层次不明显，通体为沙质，颗粒较粗，分选性好。土壤质地主要为砂土，通气性极强，但保水性极差，几乎没有持水能力。土壤中有机质含量极低，通常不足1%，这是由于风沙土的形成过程中缺乏植被覆盖，很少有有机质的输入。土壤养分含量也很低，氮、磷、钾等养分极为匮乏，全氮含量一般在0.05%以下，速效磷和速效钾含量也处于极低水平。风沙土的肥力极低，植被生长困难，生态环境极为脆弱，容易受到风力侵蚀的影响，一旦植被遭到破坏，就容易引发土地沙漠化。沼泽土主要分布在地势低洼、排水不畅、地下水位较高的地区，如玛曲、若尔盖等地的湿地。沼泽土的成土过程主要是在长期积水和水生植被的作用下，土壤进行着强烈的腐殖质化和泥炭化过程。沼泽土的表层为泥炭层或腐泥层，厚度不等，泥炭层是由大量未完全分解的植物残体堆积而成，颜色较深，多为黑色或深褐色；腐泥层则是经过进一步分解和混合的有机物质层。其下为潜育层，由于长期处于还原状态，土壤呈现出蓝灰色或青灰色。土壤质地较为黏重，通气性较差，但保水性很强。沼泽土中有机质含量极高，一般在10%-50%之间，这是由于大量的水生植物残体在厌氧环境下难以分解，不断积累形成了丰富的有机质。土壤pH值呈酸性至中性，一般在5.5-7.0之间。土壤中氮、磷等养分含量也较高，全氮含量通常在0.5%-2.0%之间，速效磷含量相对较高，这为湿地植被的生长提供了丰富的养分来源。灰褐土主要分布在山地森林与草原植被过渡地带，如祁连山、东昆仑山等山地的中下部。灰褐土的成土过程兼具淋溶作用和钙化作用，在气候、植被等因素的综合影响下，土壤发生了复杂的物理、化学和生物变化。灰褐土的剖面具有明显的发生层次，表层为枯枝落叶层和腐殖质层，厚度适中，颜色较深，富含腐殖质；其下为淀积层，有明显的黏粒和铁锰氧化物淀积现象，颜色较深，质地较为黏重；再往下是母质层。土壤质地多为壤土或黏壤土，通气性和保水性较好。土壤有机质含量中等，一般在3%-8%之间。土壤pH值呈中性至微碱性，一般在7.0-8.0之间。土壤中氮、磷、钾等养分含量较为丰富，全氮含量一般在0.2%-0.5%之间，速效磷和速效钾含量也相对较高，能够满足森林与草原过渡地带植被的生长需求。不同土壤类型对痕量元素的含量和分布具有显著影响。高山草甸土由于其较高的有机质含量和良好的土壤结构，对痕量元素具有较强的吸附和固定能力。研究表明，高山草甸土中锌、铜等痕量元素的含量相对较高，这可能与土壤中的有机质与痕量元素形成络合物，从而增加了痕量元素在土壤中的积累有关。同时，高山草甸土的弱酸性环境也有利于某些痕量元素的溶解和释放，提高了其生物有效性。例如，在高山草甸土中，锌元素更容易以离子态存在，被植物根系吸收利用。高山草原土由于其干旱的气候条件和较低的有机质含量，土壤对痕量元素的吸附能力较弱，痕量元素容易随降水或风力迁移。因此，高山草原土中部分痕量元素的含量相对较低。此外，土壤的碱性环境会使一些痕量元素形成难溶性化合物，降低其生物有效性。例如，在碱性的高山草原土中，铁、锰等痕量元素可能会形成氢氧化物沉淀，难以被植物吸收。风沙土由于其质地疏松、颗粒粗大，对痕量元素的吸附和保持能力极差，痕量元素容易流失。风沙土中痕量元素的含量通常很低，且分布极不均匀。在风沙土地区，由于植被稀少，土壤缺乏有机质的保护，痕量元素更容易受到风力和水力的侵蚀，导致其在土壤中的含量不断降低。沼泽土中较高的有机质含量和还原环境，使得痕量元素容易与有机质结合形成稳定的络合物，从而在土壤中积累。研究发现，沼泽土中汞、铅等痕量元素的含量相对较高。同时，沼泽土的酸性环境也会影响痕量元素的形态和生物有效性。例如，在酸性的沼泽土中，汞元素可能会以甲基汞的形式存在，甲基汞具有较强的毒性和生物富集性，对生态系统和人类健康构成潜在威胁。灰褐土的成土过程和土壤性质使其对痕量元素的吸附和释放处于一个相对平衡的状态。灰褐土中痕量元素的含量和分布受到成土母质、植被类型、气候条件等多种因素的综合影响。在森林植被覆盖较好的灰褐土区域，土壤中痕量元素的含量相对较高，这可能与森林植被的凋落物分解后向土壤中输入了较多的痕量元素有关。而在草原植被覆盖的区域，土壤中痕量元素的含量则相对较低。三、土壤痕量元素的分布特征3.1样品采集与分析方法在2023年7月至9月期间，研究团队对青藏高原东北部地区展开了系统的土壤样品采集工作。为了确保研究结果能够全面、准确地反映该地区土壤痕量元素的分布特征，采样点的选择充分考虑了多种因素。依据不同的地理分区，涵盖了祁连山山地、柴达木盆地、河湟谷地等具有显著地形差异的区域。同时，结合不同的土壤类型，如高山草甸土、高山草原土、风沙土、沼泽土、灰褐土等，在各类土壤分布典型区域设置采样点。此外，还考虑了土地利用方式的差异，包括天然草地、农田、林地以及受人类活动影响较小的自然保护区等不同土地利用类型区域。在采样过程中，采用了网格布点法与典型样点相结合的方式。首先，利用地理信息系统（GIS）技术，将研究区域划分为多个10km×10km的网格，在每个网格内随机选取1-2个采样点。对于一些具有特殊地质背景、地形地貌或受到人类活动强烈影响的区域，则作为典型样点进行重点采样。例如，在靠近矿山开采区、交通干线以及农业集约化生产区等人类活动影响显著的区域，加密采样点的设置，以更好地研究人类活动对土壤痕量元素分布的影响。每个采样点设置3-5个重复样，以提高数据的可靠性和代表性。使用不锈钢土钻采集0-20cm的表层土壤样品，该深度范围是土壤中生物活动最为频繁、痕量元素与外界环境交换活跃的区域，能够较好地反映土壤痕量元素的近期动态变化。采集时，小心剔除土壤中的植物根系、石块、残叶等杂物，以避免这些杂质对痕量元素分析结果的干扰。将采集好的土壤样品装入干净的聚乙烯塑料袋中，用记号笔清晰标记好采样地点、时间、经纬度、土壤类型、土地利用方式等详细信息。对于深层土壤样品，按照20-40cm、40-60cm、60-80cm、80-100cm的深度间隔进行分层采集，以研究痕量元素在土壤剖面中的垂直分布特征。在采样区域，同步使用高精度的GPS定位仪记录采样点的经纬度信息，利用海拔仪测量海拔高度，通过温湿度记录仪记录当时的气温、湿度等气候条件，并详细记录地形地貌（如山地、平原、河谷等）、植被类型（如高山草甸、草原、灌丛等）等环境信息。样品采集完成后，迅速将其带回实验室进行处理。首先，将土壤样品放置在通风良好、无阳光直射的室内自然风干，这一过程通常需要7-10天，具体时间取决于当地的气候条件和土壤湿度。待土壤样品完全风干后，使用2mm尼龙筛进行过筛，去除残留的植物根系、砂砾、小石块等杂质，确保后续分析样品的均一性。对于需要测定全量痕量元素的样品，采用盐酸-硝酸-氢氟酸-高氯酸（HCl-HNO₃-HF-HClO₄）四酸消解体系，利用电热板进行消解。准确称取0.5g过筛后的土壤样品于聚四氟乙烯坩埚中，加入10mL盐酸，在电热板上低温加热（约100℃），使样品初步分解，蒸发至溶液体积约为5mL时，取下稍冷。接着加入5mL硝酸，继续加热消解，待反应缓和后，加入3mL氢氟酸和2mL高氯酸，加盖后加热，使样品充分消解。期间适当升高温度（约150℃），直至冒高氯酸白烟，坩埚内溶液变为无色或淡黄色。此时，取下坩埚，稍冷后加入1mL硝酸溶液（1+1），温热溶解残渣，将消解液转移至50mL容量瓶中，用去离子水定容至刻度，摇匀备用。消解后的溶液使用电感耦合等离子体质谱仪（ICP-MS）测定Cd、Hg、Pb、Zn、Cu、As等痕量元素的含量。ICP-MS具有灵敏度高、分析速度快、可同时测定多种元素等优点，能够准确测定土壤中痕量元素的含量，其检出限可达ng/L级。在测定过程中，采用国家标准物质（如GBW07401-GBW07408等）进行质量控制，每批次样品分析时均插入一定比例的空白样和重复样。国家标准物质的测定结果与标准值的相对误差控制在±10%以内，重复样的相对标准偏差（RSD）控制在±5%以内，以确保分析结果的准确性和可靠性。对于测定痕量元素赋存形态的样品，采用Tessier连续提取法。Tessier连续提取法将土壤中的痕量元素依次提取为可交换态、碳酸盐结合态、铁锰氧化物结合态、有机结合态和残渣态等五种化学形态。具体步骤如下：准确称取1.0g过筛后的土壤样品于50mL离心管中，加入8mL1mol/LMgCl₂溶液（pH=7.0），在25℃恒温振荡器上振荡1h，然后以3000r/min的转速离心15min，将上清液转移至塑料瓶中，此为可交换态提取液。在离心管中加入8mL1mol/LNaOAc溶液（pH=5.0），重复上述振荡、离心步骤，得到碳酸盐结合态提取液。接着，加入20mL0.04mol/LNH₂OH・HCl溶液（用25%HAc调节pH=2.0），在96℃恒温水浴锅中振荡6h，离心后得到铁锰氧化物结合态提取液。再加入5mL0.02mol/LHNO₃和5mL30%H₂O₂（用HNO₃调节pH=2.0），在85℃恒温水浴锅中加热2h，期间不时振荡，然后加入5mL3.2mol/LNH₄OAc溶液（用HNO₃调节pH=2.0），定容至25mL，振荡30min，离心后得到有机结合态提取液。最后，将离心管中的残渣转移至聚四氟乙烯坩埚中，采用与全量痕量元素分析相同的四酸消解体系进行消解，得到残渣态提取液。使用原子吸收光谱仪（AAS）或电感耦合等离子体发射光谱仪（ICP-AES）测定各形态中痕量元素的含量。AAS和ICP-AES具有较高的灵敏度和准确性，能够满足痕量元素形态分析的要求。同样，在分析过程中，通过加入标准物质和重复样进行质量控制，确保分析结果的可靠性。3.2痕量元素含量与统计分析通过对青藏高原东北部采集的[X]个土壤样品进行分析，得到了该地区土壤中多种痕量元素的含量数据。表1展示了Cd、Hg、Pb、Zn、Cu、As等主要痕量元素的含量统计结果，包括最小值、最大值、平均值、标准差等。表1青藏高原东北部土壤痕量元素含量统计结果（mg/kg）痕量元素最小值最大值平均值标准差变异系数（%）Cd0.010.560.120.0866.67Hg0.0020.150.030.0266.67Pb10.2385.6735.4615.6744.20Zn35.67180.5685.4330.2135.36Cu15.4565.3230.2510.5634.91As5.6725.3412.354.5637.09由表1可知，该地区土壤中Cd含量范围为0.01-0.56mg/kg，平均值为0.12mg/kg，标准差为0.08mg/kg，变异系数高达66.67%，这表明Cd含量在不同采样点之间存在较大差异。Hg含量范围是0.002-0.15mg/kg，平均值为0.03mg/kg，标准差为0.02mg/kg，变异系数同样为66.67%，说明Hg含量的空间变异性也较为显著。Pb含量在10.23-85.67mg/kg之间，平均含量为35.46mg/kg，标准差为15.67mg/kg，变异系数为44.20%。Zn含量范围为35.67-180.56mg/kg，平均值是85.43mg/kg，标准差为30.21mg/kg，变异系数为35.36%。Cu含量在15.45-65.32mg/kg之间，平均含量为30.25mg/kg，标准差为10.56mg/kg，变异系数为34.91%。As含量范围是5.67-25.34mg/kg，平均值为12.35mg/kg，标准差为4.56mg/kg，变异系数为37.09%。将该地区土壤痕量元素含量与全国土壤背景值进行对比（表2），可以发现该地区土壤中Cd含量略高于全国土壤背景值（0.097mg/kg），这可能与该地区特殊的地质背景或人类活动影响有关。Hg含量与全国土壤背景值（0.065mg/kg）相比较低。Pb含量略低于全国土壤背景值（26.0mg/kg）。Zn含量略高于全国土壤背景值（74.2mg/kg）。Cu含量略低于全国土壤背景值（22.6mg/kg）。As含量略高于全国土壤背景值（11.2mg/kg）。这些差异反映了青藏高原东北部土壤痕量元素含量的独特性，可能受到成土母质、气候条件、地形地貌以及人类活动等多种因素的综合影响。表2青藏高原东北部土壤痕量元素含量与全国土壤背景值对比（mg/kg）痕量元素青藏高原东北部土壤平均值全国土壤背景值Cd0.120.097Hg0.030.065Pb35.4626.0Zn85.4374.2Cu30.2522.6As12.3511.2进一步分析不同土壤类型中痕量元素的含量差异（表3），高山草甸土中Cd、Hg、Pb、Zn、Cu、As的平均含量分别为0.15mg/kg、0.04mg/kg、40.23mg/kg、95.67mg/kg、35.45mg/kg、14.56mg/kg。高山草原土中这些痕量元素的平均含量依次为0.10mg/kg、0.02mg/kg、30.12mg/kg、75.34mg/kg、25.67mg/kg、10.23mg/kg。风沙土中痕量元素含量相对较低，Cd、Hg、Pb、Zn、Cu、As的平均含量分别为0.05mg/kg、0.01mg/kg、15.34mg/kg、45.67mg/kg、15.45mg/kg、7.67mg/kg。沼泽土中痕量元素含量较高，Cd、Hg、Pb、Zn、Cu、As的平均含量分别为0.20mg/kg、0.05mg/kg、45.67mg/kg、105.67mg/kg、40.23mg/kg、16.78mg/kg。灰褐土中痕量元素含量介于高山草甸土和高山草原土之间，Cd、Hg、Pb、Zn、Cu、As的平均含量分别为0.13mg/kg、0.03mg/kg、35.45mg/kg、85.43mg/kg、30.25mg/kg、12.34mg/kg。通过方差分析可知，不同土壤类型中Cd、Hg、Pb、Zn、Cu、As含量均存在显著差异（P<0.05）。表3不同土壤类型中痕量元素含量统计结果（mg/kg）土壤类型CdHgPbZnCuAs高山草甸土0.15±0.090.04±0.0240.23±18.6795.67±35.2135.45±12.5614.56±5.67高山草原土0.10±0.060.02±0.0130.12±12.3475.34±25.6725.67±8.5610.23±3.45风沙土0.05±0.030.01±0.00515.34±5.6745.67±15.2115.45±5.457.67±2.34沼泽土0.20±0.120.05±0.0345.67±20.12105.67±40.2140.23±15.6716.78±6.78灰褐土0.13±0.080.03±0.0235.45±15.6785.43±30.2130.25±10.5612.34±4.56不同地形地貌区域土壤痕量元素含量也存在明显差异（表4）。在山地地区，Cd、Hg、Pb、Zn、Cu、As的平均含量分别为0.14mg/kg、0.04mg/kg、38.56mg/kg、90.23mg/kg、33.45mg/kg、13.67mg/kg。河谷地区这些痕量元素的平均含量依次为0.11mg/kg、0.03mg/kg、32.45mg/kg、80.56mg/kg、28.67mg/kg、11.56mg/kg。平原地区痕量元素含量相对较低，Cd、Hg、Pb、Zn、Cu、As的平均含量分别为0.09mg/kg、0.02mg/kg、28.67mg/kg、75.34mg/kg、25.45mg/kg、10.23mg/kg。方差分析表明，不同地形地貌区域土壤中Cd、Hg、Pb、Zn、Cu、As含量差异显著（P<0.05）。表4不同地形地貌区域土壤痕量元素含量统计结果（mg/kg）地形地貌CdHgPbZnCuAs山地0.14±0.090.04±0.0238.56±17.6790.23±32.2133.45±11.5613.67±5.23河谷0.11±0.070.03±0.0132.45±13.6780.56±28.6728.67±9.5611.56±4.23平原0.09±0.060.02±0.0128.67±11.3475.34±25.6725.45±8.5610.23±3.45综上所述，青藏高原东北部土壤中痕量元素含量在不同采样点、不同土壤类型以及不同地形地貌区域均存在显著差异，这些差异为进一步研究土壤痕量元素的来源、迁移转化规律以及生态效应奠定了基础。3.3空间分布特征为深入探究青藏高原东北部土壤痕量元素的空间分布规律，运用地理信息系统（GIS）技术中的克里金插值法，对该地区土壤中Cd、Hg、Pb、Zn、Cu、As等痕量元素的含量数据进行空间插值处理，绘制出各痕量元素的含量空间分布图（图1-图6）。克里金插值法是一种基于区域化变量理论的最优内插法，它充分考虑了样本点之间的空间相关性，能够较为准确地对未采样点的元素含量进行估计。在进行克里金插值时，首先对各痕量元素的含量数据进行半变异函数分析，确定其空间变异结构，然后选择合适的理论模型（如球状模型、指数模型、高斯模型等）进行拟合。根据拟合结果，获取半变异函数的参数，如块金值、基台值、变程等，这些参数反映了元素含量在空间上的变异特征。利用确定好的半变异函数模型和参数，进行克里金插值计算，生成各痕量元素的含量空间分布图。[此处插入图1-图6，分别为Cd、Hg、Pb、Zn、Cu、As含量空间分布图]从图1中Cd含量空间分布图可以看出，该地区土壤中Cd含量呈现出明显的空间分异特征。在祁连山部分山地以及柴达木盆地周边的局部区域，Cd含量相对较高，超过了0.2mg/kg。这些区域可能受到特殊地质条件的影响，成土母质中Cd含量较高，为土壤提供了丰富的Cd来源。例如，祁连山地区经历了复杂的地质构造运动和岩浆活动，岩石中的Cd元素在风化作用下释放出来，进入土壤中，导致土壤中Cd含量升高。此外，柴达木盆地周边的一些地区可能存在富含Cd的矿产资源，在长期的地质演化过程中，矿产资源的风化和侵蚀使得Cd元素逐渐富集在土壤中。而在河湟谷地以及部分平原地区，Cd含量相对较低，多在0.05mg/kg以下。河湟谷地地势相对平坦，土壤主要由河流冲积物形成，成土母质中Cd含量较低，且该地区人类活动频繁，土壤受到的扰动较大，可能导致Cd元素的稀释和扩散，使得土壤中Cd含量处于较低水平。图2显示的Hg含量空间分布也存在显著差异。在玛曲、若尔盖等沼泽湿地分布区域，Hg含量较高，部分区域超过了0.05mg/kg。沼泽湿地具有特殊的生态环境，其较高的有机质含量和还原环境有利于Hg的富集。一方面，沼泽湿地中大量的水生植物残体在厌氧环境下难以分解，不断积累形成丰富的有机质，这些有机质中含有大量的有机官能团，能够与Hg发生络合、吸附等作用，从而增加了Hg在土壤中的积累。另一方面，在还原环境下，Hg可能会被还原为更易被土壤吸附的形态，进一步促进了Hg在土壤中的富集。在柴达木盆地内部以及一些干旱的草原地区，Hg含量较低，多在0.01mg/kg以下。这些地区气候干旱，降水稀少，土壤中水分含量低，不利于Hg的迁移和转化，且土壤中有机质含量较低，对Hg的吸附能力较弱，导致Hg含量处于较低水平。观察图3中Pb含量空间分布图，发现祁连山山地、东昆仑山地等山区的Pb含量相对较高，部分区域超过了40mg/kg。山区的地质构造复杂，岩石类型多样，其中一些岩石可能富含Pb元素，在长期的风化作用下，岩石中的Pb元素释放到土壤中，使得土壤中Pb含量升高。同时，山区的地形起伏较大，土壤侵蚀作用相对较强，可能导致富含Pb的土壤颗粒在重力和水流的作用下发生迁移和再分配，进一步影响了Pb在土壤中的空间分布。而在河湟谷地、黄河沿岸等人类活动密集的区域，虽然人类活动（如工业排放、交通尾气等）可能会向环境中释放Pb，但由于这些地区土壤受到频繁的耕作、灌溉等活动的影响，土壤的混合程度较高，使得Pb在土壤中的分布相对较为均匀，整体含量并未出现明显的高值区。Zn含量空间分布如图4所示，在祁连山、昆仑山等高海拔山区以及沼泽湿地周边区域，Zn含量较高，一般在100mg/kg以上。高海拔山区的岩石风化作用强烈，成土母质中富含Zn元素，为土壤提供了充足的Zn来源。沼泽湿地周边区域由于土壤中有机质含量较高，对Zn具有较强的吸附和固定能力，使得Zn在这些区域相对富集。在风沙土分布广泛的柴达木盆地内部以及部分干旱草原地区，Zn含量较低，多在50mg/kg以下。风沙土质地疏松，颗粒粗大，对Zn的吸附和保持能力较差，且这些地区气候干旱，植被稀少，土壤中有机质含量低，无法有效地固定Zn元素，导致Zn容易流失，含量较低。图5展示的Cu含量空间分布表明，在高山草甸土和沼泽土分布区域，Cu含量相对较高，部分区域超过了35mg/kg。高山草甸土中丰富的有机质和良好的土壤结构有利于Cu的吸附和积累，而沼泽土的高有机质含量和特殊的还原环境也使得Cu更容易与有机质结合，从而在土壤中富集。在风沙土和部分高山草原土分布区域，Cu含量较低，多在20mg/kg以下。风沙土的特性导致其对Cu的吸附能力极弱，而高山草原土由于气候干旱，土壤中微生物活动相对较弱，对Cu的转化和固定作用有限，使得Cu含量处于较低水平。由图6As含量空间分布图可知，在祁连山山地、玛曲和若尔盖等区域，As含量较高，部分区域超过了15mg/kg。这些地区的地质背景复杂，可能存在富含As的矿物，在岩石风化过程中，As元素释放到土壤中。此外，玛曲和若尔盖等地的湿地环境可能通过生物地球化学过程影响As的迁移转化，使得As在土壤中积累。在河湟谷地和部分平原地区，As含量相对较低，多在10mg/kg以下。河湟谷地人类活动频繁，土壤的理化性质受到人类活动的干扰较大，可能改变了As在土壤中的迁移和分布规律，导致其含量相对较低。通过对不同区域土壤痕量元素含量的统计分析发现，山地地区土壤中Cd、Hg、Pb、Zn、Cu、As等痕量元素的平均含量普遍高于河谷和平原地区。山地地区复杂的地质构造和较高的海拔导致岩石风化作用强烈，成土母质中的痕量元素更易释放到土壤中，同时山地地区植被覆盖度相对较高，土壤有机质含量也较高，有利于痕量元素的吸附和固定。河谷地区由于河流的冲积作用，土壤质地较为疏松，且人类活动对土壤的扰动较大，使得痕量元素在土壤中的分布相对较为均匀，含量相对较低。平原地区地势平坦，土壤主要由河流冲积物或风积物形成，成土母质中痕量元素含量相对较低，且该地区农业活动频繁，长期的施肥、灌溉等活动可能会稀释土壤中的痕量元素，导致其含量相对较低。综上所述，青藏高原东北部土壤痕量元素的空间分布受到地质背景、地形地貌、土壤类型、气候条件以及人类活动等多种因素的综合影响。不同痕量元素在不同区域的分布特征存在明显差异，深入了解这些分布特征及其影响因素，对于全面认识该地区土壤痕量元素的环境地球化学行为具有重要意义。3.4与其他地区的对比将青藏高原东北部土壤痕量元素含量与国内其他典型地区进行对比，结果如表5所示。与黄土高原地区相比，青藏高原东北部土壤中Cd含量略高，黄土高原地区土壤Cd平均含量约为0.08mg/kg，而本研究区为0.12mg/kg。这可能是由于青藏高原东北部特殊的地质背景，其成土母质中Cd含量相对较高，在长期的风化过程中，更多的Cd元素释放到土壤中。此外，该地区相对较少的人类活动干扰，使得土壤中Cd没有受到明显的稀释作用，从而导致其含量相对较高。在Hg含量方面，黄土高原地区平均值为0.04mg/kg，略高于本研究区的0.03mg/kg。黄土高原地区人类活动历史悠久，工业活动、煤炭燃烧等过程可能向环境中释放了一定量的Hg，使得土壤中Hg有一定程度的积累。而青藏高原东北部地区人口密度较低，工业活动相对较少，Hg的人为输入较少，因此土壤Hg含量相对较低。与长三角地区相比，青藏高原东北部土壤中Pb、Zn、Cu含量均较低。长三角地区作为我国经济发达地区，工业密集，交通繁忙，大量的工业废气、废水排放以及汽车尾气排放等，使得土壤中Pb、Zn、Cu等痕量元素含量显著增加。例如，长三角地区土壤Pb平均含量可达45mg/kg，Zn为100mg/kg，Cu为35mg/kg，而青藏高原东北部地区相应元素含量分别为35.46mg/kg、85.43mg/kg、30.25mg/kg。此外，长三角地区农业生产中大量使用化肥、农药，其中可能含有一定量的痕量元素，也会导致土壤中这些元素的积累。而青藏高原东北部地区以畜牧业为主，农业活动相对较少，化肥、农药使用量低，这也是其土壤中相关痕量元素含量较低的原因之一。与西南喀斯特地区相比，青藏高原东北部土壤中As含量相对较低。西南喀斯特地区特殊的地质条件，岩石中富含As等元素，在强烈的岩溶作用下，岩石中的As大量释放进入土壤，使得该地区土壤As平均含量高达15mg/kg，而本研究区为12.35mg/kg。同时，喀斯特地区的矿业活动，如雄黄矿、金矿等的开采，也会导致大量的As释放到周边土壤环境中，进一步增加了土壤As含量。相比之下，青藏高原东北部地区矿业活动相对较少，且地质条件不利于As的富集，因此土壤As含量较低。表5青藏高原东北部与其他地区土壤痕量元素含量对比（mg/kg）地区CdHgPbZnCuAs青藏高原东北部0.120.0335.4685.4330.2512.35黄土高原0.080.0432.5680.2328.6710.56长三角地区0.100.035451003513.23西南喀斯特地区0.110.0338.6790.5632.4515与国外一些地区相比，青藏高原东北部土壤痕量元素含量也存在差异。例如，与美国中西部平原地区相比，该地区土壤中Cd含量较低，美国中西部平原部分地区土壤Cd平均含量可达0.2mg/kg，而青藏高原东北部为0.12mg/kg。这可能与美国农业生产中大量使用含Cd的磷肥有关，长期的磷肥施用导致土壤中Cd逐渐积累。而青藏高原东北部地区农业生产相对粗放，磷肥使用量较少，减少了Cd的输入。在Pb含量方面，美国中西部平原地区土壤Pb平均含量约为25mg/kg，低于青藏高原东北部地区的35.46mg/kg。这可能是由于青藏高原东北部地区地质背景中Pb含量相对较高，同时交通等人类活动虽然相对较少，但局部地区如靠近交通干线区域，汽车尾气排放等也会导致土壤中Pb含量有所增加。与欧洲阿尔卑斯山区相比，青藏高原东北部土壤中Hg含量较低。阿尔卑斯山区由于工业废气的长距离传输以及历史上的矿业活动，土壤中Hg有一定程度的积累，平均含量可达0.05mg/kg，而本研究区为0.03mg/kg。此外，阿尔卑斯山区旅游业发达，旅游活动带来的人类干扰也可能对土壤Hg含量产生影响。而青藏高原东北部地区相对封闭，受工业废气传输影响较小，人类旅游活动相对较少，使得土壤Hg含量维持在较低水平。综上所述，青藏高原东北部土壤痕量元素含量与其他地区存在明显差异，这些差异主要受到地质背景、气候条件、人类活动等多种因素的综合影响。地质背景决定了成土母质中痕量元素的初始含量，气候条件影响着元素的迁移转化过程，而人类活动则通过工业排放、农业生产、交通等方式改变了土壤痕量元素的输入和分布。深入了解这些差异及其影响因素，有助于全面认识青藏高原东北部土壤痕量元素的环境地球化学行为，为区域生态环境保护和可持续发展提供科学依据。四、土壤痕量元素的赋存形态4.1赋存形态分析方法痕量元素在土壤中的赋存形态对其环境行为和生态效应有着关键影响，不同赋存形态的痕量元素，其迁移性、生物有效性和毒性存在显著差异。例如，可交换态痕量元素容易被植物吸收，对生物的毒性较大；而残渣态痕量元素则相对稳定，难以被生物利用，毒性较小。因此，准确分析土壤痕量元素的赋存形态对于深入了解其环境地球化学行为至关重要。在本研究中，采用Tessier连续提取法对青藏高原东北部土壤中痕量元素的赋存形态进行分析。Tessier连续提取法是一种经典的化学形态分析方法，由Tessier等人于1979年提出，该方法基于不同化学试剂对土壤中痕量元素的选择性溶解能力，将痕量元素依次提取为可交换态、碳酸盐结合态、铁锰氧化物结合态、有机结合态和残渣态等五种化学形态。这五种形态在环境中的稳定性和生物可利用性各不相同，可交换态是最容易被生物吸收利用的形态，它通过离子交换等物理作用与土壤颗粒表面结合，在环境条件发生微小变化时，如土壤溶液中离子浓度的改变，就可能被释放出来进入环境；碳酸盐结合态在土壤pH值发生变化时，容易发生溶解和释放，其含量受到土壤中碳酸盐含量的影响；铁锰氧化物结合态与土壤中的铁锰氧化物通过表面吸附、共沉淀等作用结合在一起，在氧化还原条件改变时，可能会发生形态转化；有机结合态主要与土壤中的有机质通过络合、螯合等作用结合，有机质的分解和转化会影响该形态痕量元素的稳定性；残渣态则是与土壤矿物晶格紧密结合，通常被认为是最稳定的形态，难以被生物利用，只有在强烈的化学作用下，如高温熔融、强酸强碱消解，才可能释放出来。具体实验步骤如下：准确称取1.0g过筛后的土壤样品于50mL离心管中，加入8mL1mol/LMgCl₂溶液（pH=7.0），在25℃恒温振荡器上振荡1h，使可交换态痕量元素与MgCl₂溶液发生离子交换反应，然后以3000r/min的转速离心15min，将上清液转移至塑料瓶中，此为可交换态提取液。在离心管中加入8mL1mol/LNaOAc溶液（pH=5.0），重复上述振荡、离心步骤，使碳酸盐结合态痕量元素在酸性条件下溶解进入溶液，得到碳酸盐结合态提取液。接着，加入20mL0.04mol/LNH₂OH・HCl溶液（用25%HAc调节pH=2.0），在96℃恒温水浴锅中振荡6h，利用NH₂OH・HCl的还原作用，将铁锰氧化物结合态痕量元素从铁锰氧化物中溶解出来，离心后得到铁锰氧化物结合态提取液。再加入5mL0.02mol/LHNO₃和5mL30%H₂O₂（用HNO₃调节pH=2.0），在85℃恒温水浴锅中加热2h，期间不时振荡，利用H₂O₂的氧化作用和HNO₃的溶解作用，将有机结合态痕量元素从有机质中释放出来，然后加入5mL3.2mol/LNH₄OAc溶液（用HNO₃调节pH=2.0），定容至25mL，振荡30min，离心后得到有机结合态提取液。最后，将离心管中的残渣转移至聚四氟乙烯坩埚中，采用与全量痕量元素分析相同的盐酸-硝酸-氢氟酸-高氯酸（HCl-HNO₃-HF-HClO₄）四酸消解体系进行消解，得到残渣态提取液。使用原子吸收光谱仪（AAS）或电感耦合等离子体发射光谱仪（ICP-AES）测定各形态中痕量元素的含量。AAS和ICP-AES具有较高的灵敏度和准确性，能够满足痕量元素形态分析的要求。在分析过程中，为确保分析结果的准确性和可靠性，采取了一系列质量控制措施。每批次实验均设置空白样品，以检测实验过程中是否存在试剂污染或其他干扰因素。同时，插入国家标准物质（如GBW07401-GBW07408等）进行分析，国家标准物质的测定结果与标准值的相对误差控制在±10%以内。对同一样品进行多次平行测定，重复样的相对标准偏差（RSD）控制在±5%以内。通过这些质量控制措施，有效保证了土壤痕量元素赋存形态分析结果的质量，为后续深入研究土壤痕量元素的迁移转化规律和生态效应奠定了坚实基础。4.2不同元素的赋存形态特征通过Tessier连续提取法对青藏高原东北部土壤样品进行分析，得到了Cd、Hg、Pb、Zn、Cu、As等主要痕量元素的赋存形态分布特征，结果如表6所示。表6青藏高原东北部土壤主要痕量元素赋存形态分布（%）痕量元素可交换态碳酸盐结合态铁锰氧化物结合态有机结合态残渣态Cd15.23±3.6718.56±4.2325.67±5.6712.34±3.4528.20±6.78Hg5.67±2.348.56±2.5618.67±4.2335.45±7.6731.65±8.56Pb8.67±2.5612.34±3.4535.45±7.6715.67±4.2327.87±6.78Zn10.23±3.4515.67±4.2338.56±8.5618.67±5.6716.87±4.56Cu9.56±3.2313.45±3.5632.45±7.6720.67±6.7823.87±5.67As6.78±2.4510.23±3.2328.67±6.7816.78±4.5637.54±8.56从表6可以看出，不同痕量元素在土壤中的赋存形态存在明显差异。Cd的赋存形态中，铁锰氧化物结合态占比最高，为25.67%±5.67%，其次是残渣态，占比28.20%±6.78%，可交换态和碳酸盐结合态分别占比15.23%±3.67%和18.56%±4.23%，有机结合态占比相对较低，为12.34%±3.45%。较高比例的铁锰氧化物结合态表明，Cd在土壤中与铁锰氧化物有较强的结合能力，可能是由于铁锰氧化物具有较大的比表面积和表面电荷，能够通过吸附、共沉淀等作用与Cd结合。残渣态Cd占比较高，说明部分Cd与土壤矿物晶格紧密结合，相对稳定，不易被生物利用。可交换态和碳酸盐结合态的Cd具有一定的活性，在环境条件改变时，如土壤pH值、氧化还原电位变化，可能会释放出来，增加其生物可利用性和迁移性。Hg的赋存形态中，有机结合态占比最高，达到35.45%±7.67%，残渣态占比31.65%±8.56%，铁锰氧化物结合态占比18.67%±4.23%，可交换态和碳酸盐结合态占比较低，分别为5.67%±2.34%和8.56%±2.56%。Hg与有机质的强结合能力是其有机结合态占比高的主要原因，土壤中的有机质含有丰富的有机官能团，如羧基、羟基、氨基等，能够与Hg形成稳定的络合物。这使得Hg在土壤中的迁移性和生物可利用性在很大程度上受到有机质含量和性质的影响。残渣态Hg占比较高，也表明部分Hg与土壤矿物紧密结合，相对稳定。较低比例的可交换态和碳酸盐结合态说明，Hg在土壤中的活性相对较低，不易被生物直接吸收利用。Pb的赋存形态中，铁锰氧化物结合态占比最高，为35.45%±7.67%，其次是残渣态，占比27.87%±6.78%，有机结合态占比15.67%±4.23%，可交换态和碳酸盐结合态分别占比8.67%±2.56%和12.34%±3.45%。铁锰氧化物结合态Pb占比高，表明Pb在土壤中与铁锰氧化物的亲和力较强，可能通过表面吸附、共沉淀等方式与铁锰氧化物结合。残渣态Pb的存在说明，部分Pb与土壤矿物晶格结合紧密，稳定性较高。有机结合态Pb的存在表明，土壤中的有机质对Pb也有一定的吸附和固定作用。可交换态和碳酸盐结合态的Pb具有一定的活性，在环境条件变化时，可能会释放出来，影响其生物有效性和迁移性。Zn的赋存形态中，铁锰氧化物结合态占比最高，为38.56%±8.56%，其次是有机结合态，占比18.67%±5.67%，碳酸盐结合态占比15.67%±4.23%，可交换态占比10.23%±3.45%，残渣态占比相对较低，为16.87%±4.56%。较高比例的铁锰氧化物结合态说明，Zn在土壤中与铁锰氧化物有较强的结合能力，这可能与Zn的化学性质以及铁锰氧化物的表面特性有关。有机结合态Zn的存在表明，土壤中的有机质对Zn也有一定的吸附和固定作用。碳酸盐结合态和可交换态的Zn具有一定的活性，在土壤环境条件改变时，可能会发生形态转化，影响其生物有效性和迁移性。残渣态Zn占比较低，说明Zn与土壤矿物晶格结合相对较弱，相对较容易释放出来。Cu的赋存形态中，铁锰氧化物结合态占比最高，为32.45%±7.67%，其次是有机结合态，占比20.67%±6.78%，碳酸盐结合态占比13.45%±3.56%，可交换态占比9.56%±3.23%，残渣态占比23.87%±5.67%。铁锰氧化物结合态Cu占比高，表明Cu在土壤中与铁锰氧化物有较强的亲和力，可能通过化学作用与铁锰氧化物结合。有机结合态Cu的存在说明，土壤中的有机质对Cu也有一定的吸附和固定作用。碳酸盐结合态和可交换态的Cu具有一定的活性，在环境条件变化时，可能会发生形态转化，影响其生物有效性和迁移性。残渣态Cu的存在表明，部分Cu与土壤矿物晶格结合紧密，相对稳定。As的赋存形态中，残渣态占比最高，为37.54%±8.56%，其次是铁锰氧化物结合态，占比28.67%±6.78%，有机结合态占比16.78%±4.56%，碳酸盐结合态占比10.23%±3.23%，可交换态占比相对较低，为6.78%±2.45%。较高比例的残渣态说明，As在土壤中与土壤矿物晶格结合紧密，大部分As以相对稳定的形态存在。铁锰氧化物结合态As的存在表明，As与铁锰氧化物也有一定的结合能力。有机结合态As的存在说明，土壤中的有机质对As也有一定的吸附和固定作用。可交换态和碳酸盐结合态的As具有一定的活性，在环境条件改变时，可能会释放出来，影响其生物有效性和迁移性。不同元素的赋存形态差异主要与元素的化学性质、土壤的理化性质以及环境条件等因素有关。元素的化学性质决定了其与土壤中各种组分的亲和力，例如，Hg的亲有机性使其在土壤中容易与有机质结合，形成有机结合态Hg。土壤的理化性质，如pH值、氧化还原电位、有机质含量、铁锰氧化物含量等，对痕量元素的赋存形态也有重要影响。在酸性土壤中，一些痕量元素的溶解度增加，可能会导致其可交换态和碳酸盐结合态含量增加；而在碱性土壤中，痕量元素可能会形成难溶性化合物，增加残渣态的比例。氧化还原电位的变化会影响铁锰氧化物的形态和性质，进而影响与铁锰氧化物结合的痕量元素的赋存形态。土壤中的有机质不仅可以通过络合、螯合等作用与痕量元素结合，还可以影响土壤的结构和性质，间接影响痕量元素的赋存形态。环境条件，如温度、湿度、降水等，也会通过影响土壤的理化性质和生物活动，对痕量元素的赋存形态产生影响。研究土壤痕量元素的赋存形态，对于评估其生物有效性和生态风险具有重要意义。可交换态和碳酸盐结合态的痕量元素具有较高的生物有效性，容易被植物吸收，对生态系统的影响较大。当土壤中这两种形态的痕量元素含量较高时，可能会对植物生长产生不利影响，甚至通过食物链传递对人类健康造成威胁。例如，可交换态Cd容易被植物根系吸收，进入植物体内，影响植物的正常生长发育，并且可能在植物可食部分积累，对人体健康产生潜在危害。铁锰氧化物结合态和有机结合态的痕量元素，在一定条件下可以释放出来，转化为生物可利用态，因此也需要关注其潜在的生态风险。残渣态痕量元素相对稳定，生物有效性较低，对生态系统的直接影响较小，但在长期的地质演化和环境变化过程中，也可能会发生形态转化，释放出痕量元素，从而影响生态环境。4.3影响赋存形态的因素土壤痕量元素的赋存形态并非固定不变，而是受到多种因素的综合影响。这些因素通过改变土壤的物理、化学和生物性质，进而影响痕量元素在土壤中的存在状态和迁移转化过程。土壤pH值是影响痕量元素赋存形态的重要因素之一。土壤pH值的变化会改变土壤颗粒表面的电荷性质和数量，从而影响痕量元素与土壤颗粒之间的吸附和解吸平衡。在酸性土壤中，H⁺浓度较高，H⁺会与痕量元素发生离子交换反应，将与土壤颗粒表面结合的痕量元素交换出来，使可交换态痕量元素含量增加。例如，当土壤pH值降低时，铁锰氧化物表面的羟基会质子化，使其表面正电荷增加，对带负电荷的痕量元素阴离子的吸附能力减弱，导致铁锰氧化物结合态痕量元素含量降低。同时，酸性条件下，碳酸盐容易溶解，碳酸盐结合态痕量元素会释放出来，增加其在土壤溶液中的浓度。研究表明，在pH值为4-6的酸性土壤中，Cd的可交换态和碳酸盐结合态含量明显高于pH值为7-8的中性至碱性土壤。相反，在碱性土壤中，OH⁻浓度较高，一些痕量元素会与OH⁻结合形成氢氧化物沉淀，使残渣态痕量元素含量增加。例如，Pb在碱性条件下容易形成Pb(OH)₂沉淀，从而降低其生物有效性和迁移性。土壤有机质对痕量元素赋存形态也有着显著影响。土壤有机质含有丰富的有机官能团，如羧基（-COOH）、羟基（-OH）、氨基（-NH₂）等，这些官能团能够与痕量元素发生络合、螯合等作用，形成稳定的有机-痕量元素络合物，从而影响痕量元素的赋存形态。一方面，有机质的络合作用可以增加痕量元素在土壤溶液中的溶解度，使其更易迁移。例如，Hg与有机质形成的有机汞络合物，具有较高的水溶性和迁移性。另一方面，有机质的吸附作用可以将痕量元素固定在土壤颗粒表面，降低其生物有效性。研究发现，土壤中有机质含量越高，有机结合态痕量元素的比例也越高。在高山草甸土中，由于其有机质含量较高，Hg、Pb等痕量元素的有机结合态占比较其他土壤类型更高。此外，有机质的分解和转化过程也会影响痕量元素的赋存形态。当有机质分解时，会释放出大量的有机酸和CO₂等物质，这些物质会改变土壤的pH值和氧化还原电位，进而影响痕量元素的形态转化。黏土矿物是土壤的重要组成部分，对痕量元素的赋存形态同样有着重要影响。黏土矿物具有较大的比表面积和表面电荷，能够通过离子交换、表面吸附等作用与痕量元素结合。不同类型的黏土矿物，其晶体结构和表面性质不同，对痕量元素的吸附能力和选择性也存在差异。蒙脱石是一种膨胀性黏土矿物，其晶层间具有较大的阳离子交换容量，能够吸附大量的阳离子型痕量元素，如Cd²⁺、Pb²⁺等。研究表明，在富含蒙脱石的土壤中，Cd的可交换态和碳酸盐结合态含量相对较高，这是因为蒙脱石的阳离子交换作用使得Cd更容易被交换到土壤溶液中。而高岭石是一种非膨胀性黏土矿物，其表面电荷主要由断键产生，对痕量元素的吸附能力相对较弱。此外，黏土矿物还可以通过与有机质形成有机-黏土复合体，间接影响痕量元素的赋存形态。有机-黏土复合体对痕量元素的吸附能力更强，能够将痕量元素固定在土壤中，降低其迁移性和生物有效性。氧化还原电位（Eh）是反映土壤氧化还原状态的重要指标，对痕量元素的赋存形态有着重要影响。在氧化条件下，土壤中的铁锰氧化物主要以高价态存在，如Fe³⁺、Mn⁴⁺等，它们具有较强的氧化性，能够将一些痕量元素氧化为高价态，从而改变其赋存形态。例如，在氧化条件下，As³⁺可以被氧化为As⁵⁺，As⁵⁺更容易与铁锰氧化物结合，形成铁锰氧化物结合态As。相反，在还原条件下，铁锰氧化物被还原为低价态，如Fe²⁺、Mn²⁺等，其对痕量元素的吸附能力减弱，导致与铁锰氧化物结合的痕量元素释放出来。研究表明，在淹水条件下，土壤处于还原状态，铁锰氧化物结合态Cd、Pb等痕量元素含量降低，而可交换态和有机结合态含量增加。这是因为还原条件下，铁锰氧化物的还原溶解使得与之结合的痕量元素释放到土壤溶液中，同时，土壤中有机质的分解产生的还原性物质也会与痕量元素发生反应，改变其赋存形态。此外，土壤中的微生物活动也会对痕量元素的赋存形态产生影响。微生物可以通过代谢活动改变土壤的pH值、氧化还原电位和有机质含量等，进而影响痕量元素的迁移转化和赋存形态。一些微生物能够分泌有机酸、氨基酸等物质，这些物质可以与痕量元素发生络合、螯合作用，增加痕量元素在土壤溶液中的溶解度。例如，某些细菌分泌的柠檬酸、苹果酸等有机酸，能够与Cd、Pb等痕量元素形成稳定的络合物，提高其生物有效性。同时，微生物还可以通过生物转化作用，改变痕量元素的化学形态。例如，一些微生物能够将无机汞转化为甲基汞，甲基汞具有较强的毒性和生物富集性，其赋存形态和环境行为与无机汞有很大不同。综上所述，土壤pH值、有机质、黏土矿物、氧化还原电位以及微生物活动等因素相互作用，共同影响着青藏高原东北部土壤痕量元素的赋存形态。深入研究这些影响因素，对于全面了解土壤痕量元素的环境地球化学行为，评估其生态风险，以及制定合理的土壤污染防治措施具有重要意义。五、土壤痕量元素的来源解析5.1多元统计分析为深入探究青藏高原东北部土壤痕量元素的来源，运用多元统计分析方法对实验数据进行处理，主要包括主成分分析（PCA）和相关性分析。主成分分析能够将多个相关变量转化为少数几个互不相关的综合变量，即主成分，通过分析主成分与原始变量之间的关系，提取数据中的主要信息，从而识别土壤痕量元素的潜在来源。相关性分析则用于衡量两个变量之间线性相关程度的强弱，通过计算痕量元素之间以及痕量元素与环境因