Comprehensive Geochemical Evaluation and Influencing Factors of Topsoil Nutrients in a Farming Area of Dabie Mountain in Western Anhui, China
-
摘要:
土壤养分是影响农作物产量的重要指标，掌握中国主要农业生产区的土壤养分时空演变特征对合理使用农业土地资源具有重要的意义。本文以中国重要农耕区皖西大别山六安地区为研究区，采集浅层土壤（0~20cm）样品1295组，采用硫碳仪、X射线荧光光谱法、电感耦合等离子体发射光谱法等技术进行测试。应用统计学分析方法，在分析pH值的基础上，对总氮（TN）、总磷（TP）和钾（K）和土壤有机质（SOM）四个地球化学养分指标进行单指标养分等级划分和土壤养分综合等级评价，并探讨了各项指标的影响因素。结果表明：研究区98.88%的土壤为酸性土壤（其中8.03%为强酸性土壤），受地貌类型影响较大；TN分布以中等丰富为主，受土壤类型影响较大；TP以中等含量和缺乏为主，受地貌类型的影响较大；K是最为丰富的元素，其含量在中等及以上的比例达97.81%，主要受地貌类型的影响；SOM以中等和较缺乏为主，受土壤类型、土地利用方式和地貌类型的共同影响，且SOM含量与TN含量呈高度正相关。全区土壤养分等级以中等为主，占全区总面积的60.10%，其次是较缺乏土壤，占25.27%，南部山区养分状况优于北部平原区。整体上，研究区土壤具有养分较差、养分分布不均匀、普遍酸化的特点，部分地区应按需调整肥料的使用平衡土壤养分。可通过农田养分的管理，提高农业产量。
-
关键词:
- 表层土壤 /
- 电感耦合等离子体发射光谱法 /
- 土壤养分地球化学综合等级评价 /
- 土地利用方式 /
- 土壤类型 /
- 地貌类型
要点（1）单指标养分等级划分和土壤养分综合等级评价可用于分析土壤养分状况。
（2）皖西大别山区土壤养分地球化学综合等级以中等为主，较缺乏次之。 pH 值分布以酸性为主，TN以中等丰富为主，TP中等偏缺乏，K以丰富为主，SOM 中等偏缺乏。
（3）土地利用方式、土壤类型和地貌类型对pH 值、养分元素和有机质的空间分布有不同程度的影响。
HIGHLIGHTS(1) The single index nutrient grade and comprehensive evaluation of soil nutrient grade are often used to analyze soil nutrient status.
(2) The comprehensive grade of soil nutrient geochemistry in the Dabie Mountain area of Western Anhui Province is mainly medium, followed by relatively deficient. The soil in the area is acidic, TN content is medium abundant, TP content is medium to lacking, K content is the most abundant, and organic matter content is medium to lacking.
(3) Land use, soil type and geomorphic type have different effects on the spatial distribution of pH, nutrient elements and organic matter.
Abstract:Western Anhui is an important precision agricultural area in China, characterized by diverse landforms and multiple soil types. To determine the spatial distribution of soil nutrient elements and their controlling factors, we analyzed 5 nutrient indicators by collecting 1295 groups of shallow soil (0-20cm) samples and conducted a comprehensive soil nutrient geochemical grade evaluation on the basis of single indicator analysis. The results show that soil nutrients are poor in the cultivated lands of the northern region. The use of fertilizer should be adjusted as needed to balance the soil nutrients in some areas, and the management of farmland nutrients is suggested for improving the agricultural yield. The BRIEF REPORT is available for this paper at http://www.ykcs.ac.cn/en/article/doi/10.15898/j.ykcs.202206180117.
BRIEF REPORTSignificance: Precision agriculture, as an emerging farm management strategy, is partly tasked with ensuring optimal plant nutrition and proper fertilizer application, including accurate dosage, an appropriate nitrogen (N)-phosphorus (P)-potassium (K) ratio, and efficient storage and distribution of nutrition, while also preventing environmental problems such as soil acidification, the “chemical time bomb” effect, phosphorus load of surface waters, eutrophication, rapid siltation of canals and reservoirs, nitrate pollution of drinking water resources, and the accumulation and mobilization of potential toxic compounds, etc. Western Anhui is an important precision agricultural area in China, characterized by diverse landforms and multiple soil types. We find that the Dabie Mountain area of Western Anhui exhibits poor soil nutrient levels, particularly within the cultivated lands of the northern region. To ameliorate this, it is advisable to judiciously augment the application of organic fertilizers, nitrogen, and phosphorus. The implementation of a rotational system between paddy and dry fields is suggested for enhanced soil nutrient content. It is also recommended to regulate soil pH through the utilization of matured lime to facilitate increased grain production. Additionally, exploiting potassium-enriched land resources is recommended. This study provides crucial scientific evidence for establishing large-scale, concentrated, and contiguous national commodity grain production bases.
Methods: A total of 1295 topsoil samples (0-20cm) were collected in the Dabie Mountain area of Western Anhui. Soil pH was determined in a 1∶2.5 ratio of soil to water using a glass electrode. TN and TP were determined using the Kjeldahl method and X-ray fluorescence spectrometry, respectively. K and soil organic matter (SOM) were analyzed by inductively coupled plasma-optical emission spectrometry (ICP-OES) and ammonium ferrous sulfate volumetric method, respectively. The comprehensive geochemical grade evaluation of soil nutrients was conducted. The impact factors of these five geochemistry indexes were studied based mainly on the correlation analysis.
Data and Results: The comprehensive geochemical grade evaluation of soil nutrients in the region showed that the overall soil nutrient classifications in the region were medium, accounting for 60.10% of the total area, mainly distributed in the southern parts of Huoshan County and the northeastern forest area of Jinzhai County. The second prevalent category was moderately deficient soil, constituting 25.27%, scattered sporadically in a non-contiguous manner primarily in the northern part of the surveyed area in Yu’an district. Results from the analysis of 5 various indexes indicated that 98.88% of the soils in the study area were acidic, with a significant influence from geomorphic types. TN distribution was mainly medium to rich, significantly affected by soil types. TP content was primarily medium and deficient, with a significant influence from geomorphic types. K content was the most abundant element, with 97.81% having medium or higher content, mainly influenced by geomorphic types. SOM content was primarily medium to relatively deficient, influenced by soil types, land use patterns, and geomorphic types, and shows a highly positive correlation with TN content.
-
重金属在土壤中具有潜伏性强、迁移速率慢、难降解的特性[1]，通过食物链累积可危害人类健康，成为各界关注的热点。由于平均海拔高、特殊的地理位置和气候环境作用，青藏高原目前仍然是受人类活动影响最少的地区之一[2]，自然条件下青藏高原生态系统中重金属元素含量相对较低[3]，但是随着自然资源的开发利用、相邻地区的工业活动传输、宗教活动频繁、旅游活动等的影响，青藏高原土壤生态系统也遭受到了一定程度的重金属污染[4]，存在局部土壤重金属污染现象。根据Wu等[5]对于青藏高原重金属生态环境风险评价研究结果显示，其采集的青藏高原土壤样品中接近50%的采样区域受到中度或重度重金属污染，其中青海地区土壤污染较严重，西藏地区污染相对较弱。而引起青藏高原土壤重金属累积的源解析结果主要认为人为活动和地质成土作用是最重要来源[6]。目前青藏高原土壤重金属的研究整体上有待补充，现有的大量研究成果多集中于青海湖[7]、黑河[8]等流域的周边土壤和道路交通运输附近的[9-10]土壤重金属污染，受区域地貌、范围广泛等影响，研究人员对其他地区的关注度有限，总体上存在研究覆盖范围相对较窄或采样数据较早等问题[11-12]。
青藏高原的东部区域，在中国东南季风的边缘地区，是气候和生态环境敏感、脆弱和过渡地区。与青藏高原其他区域相比，人类活动、农业生产、工业活动对此区域的影响相对显著，但是整体对于该区域的研究热度较低，目前少有的部分研究显示青藏高原东北部分地区土壤重金属污染严重[13]，污染原因主要有工业、矿业、交通和农业影响等。事实上，土壤重金属污染往往是多种成因来源和作用途径叠加综合的结果，就目前的研究成果而言，仅仅简单地判别重金属元素来源类型已经不足以为区域土壤重金属污染治理提供足够的信息，需定量计算各类排放源对元素累积的相对贡献率，确定主要污染源。全面了解青藏高原土壤中重金属的污染模式、分布及来源，是土壤污染防治的有效前提。因此，本文针对青藏高原东部区域的典型高山农业县—泽库县土壤重金属的潜在生态风险评价及来源解析开展研究，采用电感耦合等离子体质谱法和发射光谱法（ICP-MS/OES）、原子荧光光谱法（AFS）等技术手段，结合多元统计分析方法和地理信息系统技术研究了土壤重金属元素（As、Cd、Co、Cr、Cu、Hg、Mn、Ni、Pb、Zn）分布特征规律，运用富集因子法（EF）、地累积指数（Igeo）、潜在生态危害指数（PERI）等进行土壤污染程度分析，通过主成分分析法（PCA）和绝对因子分析/多元线性回归分析法（ACPS-MLR）进行重金属源解析并定量阐释主要影响因素贡献率，确定主要污染源，为制定该地区环境保护、资源有效利用政策、措施落实提供科学支撑，以充实青藏高原东部区域土壤重金属的研究内容，补充区域土壤重金属评价的理论基础。
1. 研究区概况
青海省泽库县为典型的高山农业区，位于青海省东南部、青藏高原东北部，隶属于青海省黄南藏族自治州南部，地处西倾山北麓，麦秀山西南部，西秦岭西部与黄土高原的过渡地带。区域海拔在3500m以上，高差200~800m，全境东西较狭长，地势由东向西倾斜，经历了长期的地质演化，构造特点极其复杂，形成了全县不同海拔高度的植被、土壤、气候的地域差异和垂直变化，有高山带、亚高山带、滩地、河谷阶地、低山丘陵五大地带。气候为高原大陆性季风气候，属高原亚寒带湿润气候区。
研究区隶属秦、祁、昆复合造山链东段秦岭复合造山带西部，地处泛华夏大陆早古生代秦祁昆构造区和泛华夏大陆晚古生代藏北（羌塘）三江弧—盆系两大构造单元的结合部位[14]，位于青海省二级大地构造分区——宗隆务—青海南山晚古生代—早中生代裂陷槽内。出露地层，从老到新有中生界的三叠系、侏罗系和白垩系、新生界的古近系、新近系和第四系。从分布情况看，三叠系为研究区分布最广的地层，分布于区中、南、西北部，岩性为一套浅海滨海相的碎屑岩，以砂岩、板岩为主的浅变质岩系及少量泥灰岩与生物碎屑灰岩。其他区域地层岩性为以砂岩、砂砾岩为主的碎屑岩。第四系主要分布在山间洼地和山前洪积扇裙，河湖相黄土、灰黄色亚砂土、亚黏土层、冰碛沉积物均有广泛分布。复杂的地质构造和地层分布造成了研究区内土壤元素地球化学特征的复杂性。由于海拔偏高，研究区域的土壤发育程度整体普遍偏低，以物理风化为主[15]。
2. 实验部分
2.1 样品采集与分析
在充分考虑了土壤发育厚度、受人为影响少、采样难易度等情况下，采样点位尽量覆盖不同地质单元(图1a)、地貌(图1b)及土壤类型（图1c），对青海省典型高山农业县—泽库县内的表层土壤样品（0~20cm）进行采集。整体受青藏高原地质构成和气候条件限制，泽库县的土壤风化物剖面发育厚度有限，关键带样品尽量找寻发育良好的剖面进行采集，采集表层土壤每件样品不少于1kg，共采集到表层土壤43件。样品参照《生态地质调查技术要求（1∶50000）（试行）》（DD2019-09）中的土样保存要求运送到实验室进行测试分析，样品取样后用自封袋排空空气后密封保存，带回实验室后自然风干、研磨，过2mm筛去杂质后进行化学测试。分析过程中按照规范每15件样品中插入国家一级标准物质（GBW07447、GBW07385和GBW07387）和重复样1件用于控制分析数据的准确度和精密度。
![]()
图 1 研究区地质地貌和土壤类型分布示意图及采样点位图Figure 1. Map of geology (a), geomorphology (b) and soil types (c) in Zeku County. Soils sampled in different geological units (a), with different geomorphic types (b) and soil types (c) as far as possible.样品分析和质量控制按照《区域地球化学样品分析方法》（DZ/T 0279—2016）和《碳酸盐岩石化学分析方法 第30部分：44个元素量测定》（GB/T 14506.30—2010）执行。样品测试分析在有色金属西北矿产地质测试中心完成，其中：pH值使用玻璃电极法测试（水土比例为2.5∶1）；土壤有机碳含量使用重铬酸钾容量法测试；重金属元素As、Hg采用原子荧光光谱法（AFS）测试，测试仪器为AFS2202E原子荧光光谱仪；Cr、Cu、Mn、Ni、Zn采用iCAP 6300电感耦合等离子体发射光谱仪(ICP-OES)测试；Cd、Co、Pb采用7700X电感耦合等离子体质谱仪（ICP-MS）测试。
对重金属As、Cr、Cu、Mn、Ni、Zn、Hg含量分析，样品前处理方法主要是盐酸+硝酸+氢氟酸+高氯酸消解法；Cd、Co、Pb含量分析，样品前处理方法主要为硝酸+氢氟酸消解法。分析方法的测试检出限：Hg为0.005μg/g，As、Cr、Ni、Mn为0.2μg/g，Cu、Pb为0.5μg/g，Zn为0.03μg/g，Cd为0.021μg/g，Co为0.2μg/g。元素分析的准确度和精密度均超过95%以上，达到或优于DZ/T 0279—2016和GB/T 14506.30—2010的要求，分析数据可靠。
2.2 数据处理及评价方法
2.2.1 数据处理方法
数据采用Origin进行处理与分析，采用IBM SPSS statistics 19软件对数据进行平均值、最小值、变异系数、相关性分析及主成分分析等，采用富集因子（EF）、地累积指数（Igeo）、重金属潜在生态危害指数法（PERI）来评价重金属污染状况，并用绝对因子分析/多元线性回归（PCA-APCS-MLR）分析重金属元素的主要来源及相关贡献率。
2.2.2 评价方法
（1）富集因子法（EF）
富集因子法（Enrichment Factor, EF）是一种简便、易行的重金属污染评价方法，能较为准确地判断人类活动对土壤污染的贡献，进而从自然异常中剖离出人为异常[16]。计算公式为：
$$ EF=\dfrac{{\left(\dfrac{{C}_{{\rm{i}}}}{{C}_{\mathrm{r}\mathrm{e}\mathrm{f}}}\right)}_{{\rm{s}}}}{\left(\dfrac{{B}_{{\rm{i}}}}{{B}_{\mathrm{r}\mathrm{e}\mathrm{f}}}\right)}_{{\rm{B}}} $$ (1) 式中：（Ci/Cref）s是元素i的实测浓度与参比元素实测浓度的比值，（Bi/Bref）B是元素i背景值与参比元素背景值的比值。参比元素通常选择地壳中普遍大量存在的、人为污染源较少、化学稳定性好、挥发性低的元素，多数选用Th、Sc、Fe、Al等元素，因Th稳定性好、挥发性低、人为污染源少，本文将Th作为参比元素。
选取青海省土壤元素背景值作为元素背景值。元素富集程度受人类影响的大小分为6级（表1）：EF≤1表示土壤元素较参考背景而言并未发生富集。由于在自然成土状态下，EF均一般小于2，较高的EF可能指向较为显著的人为输入。
表 1 富集因子（EF）和地累积指数(Igeo)与污染程度关系Table 1. Relationship between enrichment factor (EF), geoaccumulation (Igeo) and contamination grades.EF 富集程度等级 Igeo 污染程度判定 EF≤1 未富集 Igeo≤0 未污染 1<EF≤2 轻微富集 0<Igeo≤1 未污染至中度污染 2<EF≤5 中度富集 1<Igeo≤2 中度污染 5<EF≤20 显著富集 2<Igeo≤3 中度至高度污染 20<EF≤40 强烈富集 3<Igeo≤4 重度污染 40< EF 极强富集 4<Igeo≤5 重度至极度污染 Igeo>5 极度污染 （2）地累积指数（Igeo）
地累积指数（Geoaccumulation Index, Igeo）最早由Müller[17]提出，用于定量评价沉积物等中重金属污染程度，计算公式如下：
$$ I_{{\rm{geo}}}={\rm{log}}_{2}[C_{{\rm{i}}}/(k\times B_{{\rm{i}}})] $$ (2) 式中：Ci为元素i的测定值（mg/kg）；Bi为元素i的背景值（mg/kg）；k为修正系数1.5。
本文选择青海省土壤元素背景值作为地球化学背景值。Igeo与污染程度之间的对应关系具体见表1。
（3） 重金属潜在生态风险指数（PERI）
潜在生态风险指数（Potential Ecological Risk Index, PERI）由瑞典著名地球化学家Hakanson[18]提出，计算公式为：
$$ {E}_{\mathrm{r}}^{\mathrm{i}} = {T}_{\mathrm{r}}^{\mathrm{i}} ×（ {C}_{\mathrm{d}}^{\mathrm{i}}/{C}_{\mathrm{l}}^{\mathrm{i}}) $$ (3) $$ {\rm{PERI}}= \sum _{i=1}^{m}{E}_{\mathrm{r}}^{\mathrm{i}} $$ (4) 式中：$ {E}_{\mathrm{r}}^{\mathrm{i}} $为重金属i的潜在生态危害系数；$ {C}_{\mathrm{d}}^{\mathrm{i}} $为重金属i的实测量（mg/kg）；$ {C}_{\mathrm{l}}^{\mathrm{i}} $为重金属i的土壤背景值（mg/kg）；$ {T}_{\mathrm{r}}^{\mathrm{i}} $为i的毒性响应系数。Hakanson[18]提出的重金属毒性水平顺序为：Hg>Cd>As>Pb=Cu>Cr>Zn，给出的毒性响应系数为：Hg=40，Cd=30，As=10，Pb=Cu=5，Cr=2，Zn=1。徐争启等[19]按照Hakanson的计算原则，将Co、Ni、Ti、Mn、V等元素进行了毒性响应系数计算，得出的结果是Co=Ni=5，V=2，Ti=Mn=1。
潜在生态危害系数（$ {E}_{\mathrm{r}}^{\mathrm{i}} $）由低到高分5个等级（表2），用来描述某一元素的生态危害程度。潜在生态危害指数（PERI）用来表示某一点多个元素的潜在生态危害综合值，分值由低到高分为4个等级（表2）。
表 2 潜在生态危害系数（$ {E}_{r}^{i} $）及危害指数（PERI）与风险程度Table 2. Potential ecological risk efficient ($ {E}_{\mathrm{r}}^{\mathrm{i}} $), potential ecological risk index (PERI) and risk level.$ {E}_{\mathrm{r}}^{\mathrm{i}} $ 危害程度 PERI 危害程度 $ {E}_{\mathrm{r}}^{\mathrm{i}} $≤40 轻微 PERI<150 轻微 40≤$ {E}_{\mathrm{r}}^{\mathrm{i}} $<80 中等 150≤PERI<300 中等 80≤$ {E}_{\mathrm{r}}^{\mathrm{i}} $<160 强 300≤PERI<600 强 160≤$ {E}_{\mathrm{r}}^{\mathrm{i}} $<320 很强 PERI≥600 很强 $ {E}_{\mathrm{r}}^{\mathrm{i}} $>320 极强 （4）绝对因子分析/多元线性回归受体模型（ACPS-MLR）
绝对因子分析/多元线性回归受体模型（APCS-MLR）[20-21]的基本原理是将因子分析的主因子得分转化为绝对主因子得分（APCS），各重金属含量再分别对所有的APCS进行多元线性回归，回归系数用于计算各个主因子对应的污染源对土壤中每件样本点位某重金属含量的贡献量。主要步骤为：首先对所有金属元素含量进行标准化；然后引入一件浓度为0的人为样本，得到0浓度样本的因子分数后，将因子分析得到的元素主因子得分减去0浓度样本的主因子分数后得到主成分因数APCS；以APCS为自变量，重金属含量为因变量进行多元线性回归分析，利用得到的回归系数计算每个重金属的污染源贡献量。具体计算方法详见Guo等[22]的描述。
3. 结果与讨论
3.1 表层土壤重金属含量变化特点
3.1.1 研究区农业土壤重金属含量特征
泽库县表层土壤中重金属含量的统计情况列于表3，其中土壤环境质量标准采用《土壤环境质量 农用地土壤污染风向管控标准（试行）》（GB 15618—2018）中其他土地利用类型在pH>7.5（研究区的土壤pH均在8以上）情况下的农用地土壤污染风险筛选值，中国土壤背景值采用《中国土壤地球化学参数》中的第二环境土壤元素地球化学值[23-24]，青海省表层土壤重金属背景值采用《中国土壤地球化学参数》中青海省表层土壤地球化学参数中的几何平均值[23]。变异系数分为三个等级，分别为：低度变异（CV<0.16）、中等变异（0.16<CV<0.36）和高度变异（CV>0.36）[25]。
表 3 研究区土壤重金属含量统计特征Table 3. Statistical characteristics of soil heavy metal concentrations.重金属含量统计值 As Cd Co Cr Cu Hg Mn Ni Pb Zn 样本数（件） 43 43 43 43 43 43 43 43 43 43 平均值(mg/kg) 28.95 0.14 13.85 63.15 23.84 0.017 774.14 30.65 23.20 74.96 标准差(mg/kg) 39.35 0.05 3.77 21.49 6.95 0.008 150.26 12.56 3.32 10.16 方差 1548.65 0.00 14.21 461.94 48.32 0.00 22579.12 157.80 11.00 103.13 峰度 26.73 6.94 7.75 19.77 12.99 6.016 2.55 20.88 0.47 2.11 偏度 4.90 2.09 2.02 3.56 2.76 1.975 0.73 3.87 0.29 0.26 最小值(mg/kg) 8.66 0.08 7.45 25.20 11.80 0.009 448.00 11.10 16.40 46.60 最大值(mg/kg) 254.00 0.36 30.20 178.00 57.70 0.049 1286.00 98.70 31.30 105.00 变异系数（CV） 1.36 0.35 0.27 0.34 0.29 0.450 0.19 0.41 0.14 0.14 土壤环境质量标准 (mg/kg) 25 0.6 - 250 100 3.4 - 190 170 300 中国土壤背景值(mg/kg) 9.1 0.15 11.7 63 23 0.050 552 26 25 67 青海省表层土壤背景值（mg/kg） 13.0 0.184 12.7 73 24 0.021 654 28 22 69 As含量的变化范围为8.66~254.00mg/kg，均值28.95mg/kg，超过了全国、青海省土壤背景值，从单件样品含量判断，42件表层样品均超过了中国土壤背景值，也明显高于青海省土壤表层背景值，只有1件样品As含量为8.66mg/kg，接近中国土壤背景值9.1mg/kg。与土壤环境质量标准相比，As含量均值超过了环境质量标准25mg/kg。变异系数可以反映区域重金属元素的分布差异。As的变异系数1.36超过1，属于高度变异，偏度和峰度较高，受外界影响相对较大，分布不均匀。Cd、Cr、Cu含量均值接近全国背景值，低于青海省表层土壤背景值含量，远低于土壤环境质量标准；变异系数均小于0.36，时空分散程度较小，说明含量分布空间变异性不大。Co、Ni、Zn含量均值均大于全国背景值和青海省背景值，远小于土壤环境质量标准要求。Co、Ni、Zn的变异系数分别为0.27、0.41、0.14，时空分散程度存在元素差异。Pb含量均值未超过中国土壤和青海省表层土壤背景值，从单件样品判断，25%左右的样品超过了中国背景值。根据土壤环境质量判断，所有样点值均未超过农用地土壤质量标准含量。从变异系数判断，变异系数在0.14以下，属于低度变异，Pb峰度和偏度变化较小，时空分散有差异但是不大。样品中Hg含量均较低，最大值0.05mg/kg，均值0.02mg/kg，小于全国土壤背景值，接近青海省背景值。其变异系数、峰度和偏度都较高，有受到外界影响的可能。Mn含量的变化范围为448.00~1286.00mg/kg，平均含量（774.14mg/kg）明显超过了全国和青海省背景值。但Mn的变异系数较小，分散度不大，分布整体比较均匀，偏度和峰度较小。综上，除As外，其他各重金属均值含量均未超过且远远小于土壤环境质量值，与中国土壤环境背景值和青海省表层土壤背景值相比，Cd、Hg含量均小于背景值，As、Mn均值含量远大于背景值，Co、Cr、Ni、Pb、Zn含量略大于背景值，Cu基本接近于背景值。
3.1.2 表层土壤重金属的时空分布特征
因受不同的自然环境和人类活动影响，泽库县不同区域土壤重金属元素地球化学含量在时空分布上差异显著（图2）。As（图2a）、Cu（图2b）在泽库县北部的浓度最高，随着逐渐南移，含量逐渐降低，中部和西南区域内的含量较南部和西部高。Cd（图2c）、Cr（图2d）浓度空间分布特征具有相似性，在偏北部区域含量最高，从北到南、从东到西的含量有逐渐减少的趋势，空间差异特征明显。Co（图2e）、Zn（图2f）、Ni（图2g）在东北部的麦秀镇周边的浓度含量最高，以圆形扩散方式向四周逐渐降低，西南区域含量略低于麦秀镇，为第二高区域，西北到西南方向带含量整体较低。泽库县整体的Hg含量（图2h）偏低，但区域含量变化较为明显，西部区域Hg含量明显低于东部地带，在西部和东部区域零星分布个别高含量点位，东北部麦秀镇附近含量最高。Pb含量明显表现出在泽库县正北部、西南角、东北角区域偏高，在西北、南部区域含量偏低（图2i）。在东部区域内整体的Mn含量（图2j）高于西部，由东向西逐渐降低，但在西南方向的县域附近含量也较高，接近东部最高值，北部含量最低。整个研究区内，北部、东北部区域呈现出明显的多种重金属含量高于其他区域的特点。根据区域土地利用方式调研显示，北部、东北部采样点位是泽库县两个重要矿业公司所在地附近，初步判断北部、东北部区域多种重金属含量高与矿业开发有较大关系。
![]()
图 2 表层土壤重金属元素（a） As、（b） Cu、（c） Cd、（d） Cr 、（e） Co、（f） Zn、（g） Ni 、（h） Hg、（i） Pb、（j） Mn地球化学空间分布图Figure 2. Spatial distribution of geochemical content with (a) As, (b) Cu, (c) Cd, (d) Cr, (e) Co, (f) Zn, (g) Ni, (h) Hg, (i) Pb and (j) Mn in surface soils for sampling sites. The larger the diameter of the circle, the greater the content. The different colors represented different levels of element content.3.1.3 与其他区域农业表层土壤重金属含量对比
将本研究与青藏高原其他区域、其他平原地区农业土壤重金属含量进行对比分析（表4）。与青海湖流域表土相比，本研究区域内的As含量明显较高，与玉树县、高速公路附近的土壤相比，As含量也明显偏高，与当地岩石本底As含量偏高有关[6]；Cd和Hg含量与其他区域对比发现明显偏低；Co、Cr、Cu、Ni、Pb、Zn含量除比青海湖地区及其流域表层土壤稍高外，均低于青藏高原其他地区；Mn含量处于高水平，由于Mn的亲石特性，其元素特征一般较为稳定，因此泽库县Mn含量高于其他青藏高原地区很可能与泽库县当地的岩石本底Mn含量高有很大关系。湖泊流域内的土壤元素地球化学迁移特性较非流域内土壤影响因素复杂，青海湖流域内的表层土壤重金属含量变化与本研究及青藏高原其他地区有明显区别。对于其他青藏高原地区的表层土壤来说，泽库县的旅游业和工业活动发展有限，与表4中列出的其他地区如高速公路附近、旅游业发展地区等相比[15,28-31]，人为活动影响有限，因此整体地除As、Mn外的其他重金属含量明显偏低，高山农业区域土壤的重金属含量有其特殊性。
表 4 研究区与其他地区表层土壤的重金属含量对比Table 4. Comparison of heavy metals content in topsoil between the study area and other areas.重金属
元素本研究区含量
(mg/kg)青藏高原其他地区重金属含量（mg/kg） 青海省
玉树县[26]青藏高原
东北部地区[27]青藏高原东北—
西南方向[15]环青海湖
地区[28]青海湖流域
表层土壤[29]青海省重要
交通沿线[30]一江两河流域
中部地区[31]As 28.95 20.95 - - 11.73 11.66 21.60 - Cd 0.14 0.24 0.68 0.17 0.62 - 0.19 0.21 Co 13.85 - 11.59 11.39 12.38 12.73 10.50 - Cr 63.15 72.84 83.10 70.84 41.35 54.17 74.60 82.95 Cu 23.84 24.18 40.74 23.92 19.33 19.72 22.00 34.67 Hg 0.017 0.046 0.280 - - - 0.050 - Mn 774.14 735.06 - 639.64 546.96 626.28 - 697.39 Ni 30.65 - 54.73 31.64 21.18 24.96 39.40 49.99 Pb 23.20 26.93 72.49 28.65 21.86 20.47 32.90 35.81 Zn 74.96 85.10 145.64 73.31 63.51 - 100.30 75.31 重金属
元素本研究区含量
(mg/kg)其他平原地区农田土壤含量（mg/kg） 三江平原[32] 江苏省[33] 浙江省[33] 淮北平原[34] 长三角地区[35] 江汉平原[36] As 28.95 16.87 10.24 7.25 12.1 8.14 - Cd 0.14 0.18 0.18 0.23 0.48 0.25 0.48 Co 13.85 - - - - - - Cr 63.15 69.83 71.49 47.84 72.24 68.84 - Cu 23.84 35.28 26.56 23.96 23.73 32.58 48.2 Hg 0.017 0.072 0.070 0.120 0.046 0.140 0.120 Mn 774.14 - - - - - - Ni 30.65 22.29 29.68 21.31 33.23 33.02 48.80 Pb 23.20 18.26 28.80 36.79 24.65 32.32 36.50 Zn 74.96 68.21 75.87 91.39 131.79 92.35 96.80 与中国其他平原地区农业土壤相比，本研究的As含量明显要高于其他平原地区，极大可能与青藏高原地区本底As含量高有较大关系。其他元素，在江汉平原、长三角地区、淮北平原的含量均大于本研究。浙江省和江苏省除了Ni元素略低外，Co、Cr、Cu、Cd、Hg、Pb、Zn等重金属含量均超过本研究。三江平原的Ni、Pb、Zn略低于本研究外，其余重金属含量明显高于本研究。综上，与人类活动频繁、经济相对发达地区相比，典型高山区域的农业土壤中重金属含量绝大部分要低于平原地区的典型农田土壤。泽库县内的农业方式是以牧业为主，种植次之，相对于其他平原地区高度发达的以种植型农业为主的农业方式，对于土壤中重金属富集产生的影响较小。
3.2 重金属污染评价
3.2.1 富集因子（EF）
计算EF时，土壤背景值采用表1中的青海省表层土壤背景值。泽库县表层土壤重金属富集程度为：As （2.33）>Mn （1.27）>Ni （1.17）>Co （1.16）>Zn （1.15）>Pb （1.11）>Cr （0.92）> Cd（0.82）（图3a）。As表现为中度富集，近29%的样点为中度及以上富集，提高了整体的富集等级，但是中、高富集样点较为分散，值得注意的是本次采集的样品只有一个点位的富集程度大于10以上，是由As引起的（图4），位于泽库县北部的多隆滩—越秀镇附近。Mn、Ni、Co、Zn、Pb表现为轻微富集，其中只有近2%~3%的点位为中度富集，未有点位出现高度富集。Cr、Cd的富集因子均小于1,与青海省表层土壤相比，没有发生明显富集。根据EF值结果，泽库县的土壤重金属元素除As外的富集程度均属于轻微富集，表明表层土壤受到的人为活动较小，As的富集程度最为突出也仅表现为中度富集，后期需持续关注。
![]()
图 3 泽库县表层土壤重金属（a）EF和（b）Igeo箱线图Figure 3. Box-line diagrams of (a) EF and (b) Igeo for heavy metals in topsoil from Zeku County.![]()
图 4 泽库县表层土壤重金属As的富集系数（EF）空间分布图Figure 4. Spatial distribution of enrichment factor (EF) of As in topsoil from Zeku County.
3.2.2 地累积指数（Igeo）
地累积指数分析结果显示：泽库县表土重金属均处于未受污染范围中度污染范围（−2.12<Igeo<3.70），污染程度为Cd（−1.02）<Hg（−1.00）<Cr（−0.86）<Cu（−0.64）<Ni（−0.53）<Pb（−0.52）<Co（−0.50）<Zn （−0.48）<Mn （−0.37）<As （0.13）（图3b）。与EF分析结果基本一致，泽库县表土重金属污染以As为主（−1.17~3.7），Cd污染分布范围广，但是污染程度普遍不高，分别有21%、11%、2%的样点处于未受污染-中度污染状态、中度污染状态、中度-高度污染状态。Pb、Zn表现为所有样点均未受污染。其他重金属结果显示大部分样点处于未受污染状态，仅有个别1~2个点位表示有轻微污染-中度污染的趋势。需要注意的是，部分地区出现复合污染的趋势，如北部的多隆滩—麦秀镇即点位（D073T)附近出现了As重度污染、Cu轻微污染的状态（图5），东北部分样点也出现Co、Cr、Ni三种元素轻微-中度复合污染的状态。
![]()
图 5 泽库县表层土壤重金属（a）As和（b）Cu地累积指数（Igeo）空间分布特征图Figure 5. Spatial distribution of Geoaccumulation index of heavy metals (a) As and (b) Cu in topsoil in Zeku County. Heavy As pollution combined with moderate Cu pollution occurred in the northern sampling sites.3.3 潜在生态风险指数（PERI）
单元素的潜在生态风险系数结果（表5）显示，除了As和Hg的个别点位的风险系数大于40外，其他8个元素的生态风险指数均小于40，属于轻度危险范围内。对于As来说，具有重度及中度危险的点分别位于东北部（D073T）和中北部（D029T）。Hg元素的生态风险指数最大值62.22，属于中度风险，但只有位于东北部的一个样点风险系数达到了该值，其他42件样点均处在低风险水平，可以看出Hg的整体生态风险趋势不高。
表 5 土壤中重金属元素的潜在生态风险系数统计Table 5. Comprehensive potential ecological risk coefficients of heavy metals in the soils.潜在生态风险系数（$ {E}_{\mathrm{r}}^{\mathrm{i}} $）统计值 As Cd Co Cr Cu Hg Mn Ni Pb Zn 平均值 14.85 15.59 3.64 1.15 3.31 21.59 0.79 3.65 3.52 0.72 中位数 9.08 14.13 3.52 1.11 3.24 19.05 0.78 3.45 3.53 0.72 众数 8.46 16.30 3.54 1.16 3.64 15.24 0.70 3.58 3.58 0.68 标准差 20.18 5.52 0.99 0.39 0.97 9.71 0.15 1.50 0.50 0.10 最小值 4.44 8.59 1.96 0.46 1.64 11.43 0.46 1.32 2.48 0.45 最大值 130.26 39.13 7.93 3.25 8.01 62.22 1.31 11.75 4.74 1.01 从多金属生态风险指数（PERI）空间分布（图6）来看，表层土壤的重金属（PERI变化范围为40~200）处于轻微-中度危害的范围内：42件样品处于轻度危险内，仅有1件样品的生态危害指数处于中度危害范围内，与As的最高生态风险系数点位对应。整体判断，与富集系数、地累积指数的结果类似，表层土壤的10个重金属元素的生态危害指数处于较小危害程度，个别样品有潜在生态风险。
![]()
图 6 潜在生态危害指数（PERI）的空间分布图Figure 6. Spatial distribution of potential ecological risk index (PERI) for the topsoil in Zeku County. In the northern site, the value of PERI showed higher than the other sites.3.4 土壤重金属来源解析
相关性分析和主成分分析是重金属源解析最为成熟的定性分析方法，两种方法的分析结果可以互为验证、互为补充，共同使用能保证源解析的定性分析更准确，不漏掉相似的源成分。PCA-APCS-MLR模型可有效地应用主成分分析结果，是目前重金属源解析重要的定量分析方法之一。本研究选用三种方法相结合的方式进行重金属源解析，充分利用各个方法的分析结果，定性与定量相结合，保证了解析结果的科学性与利用性。
3.4.1 各重金属元素间的相关性分析
相关性分析是研究变量之间相关关系，可以对诸多因子影响下某个因子的显著性进行判断的一种数理统计方法。土壤重金属间的相关性见表6，其中As-Cu、Cd-Cu、Cd-Pb、Co-Cr、Co-Cu、Co-Mn、Co-Ni、Co-Zn、Cr-Ni、Mn-Pb、Mn-Zn、Pb-Zn呈现极显著正相关（p<0.01）, Cd-Hg、Cr-Cu、Cd-Zn、Cr-Mn、Cu-Ni、Cu-Zn、Hg-Mn、Mn-Ni呈显著正相关关系（p<0.05）。显著相关的重金属，相关系数越大说明这些元素之间关系越强，可能具有共同或相似来源。元素Co、Ni、Cu、Zn由于在地球化学行为相似、赋存状态相似或相同、亲和力强，其相关性较好。Cr-Ni、Co-Ni相关系数最大分别达到了0.97、0.93，地球化学行为相似，极大可能有相同来源。As只与Cu表现出极显著的正相关关系，相关系数达到0.82，从地球化学特性来说两者都具有亲硫性，具有较大的可能性具有相同的来源或者污染特性。As-Hg、As-Ni/Cr/Cd、Hg-Cr/Cu等没有显著相关性，表明As与Hg和Ni及其相似元素等来源有差异。因10种重金属间不完全显著相关，则可能具有多种来源。
表 6 土壤重金属元素间相关性分析Table 6. Correlation analysis of heavy metals in the soils.元素 As Cd Co Cr Cu Hg Mn Ni Pb Zn As 1 Cd 0.24 1 Co 0.03 0.15 1 Cr −0.00 −0.10 0.90** 1 Cu 0.82** 0.21 0.44** 0.36* 1 Hg −0.09 0.38* 0.06 0.03 −0.18 1 Mn −0.05 0.44** 0.53** 0.32* 0.20 0.31* 1 Ni −0.02 0.01 0.93** 0.97** 0.34* 0.02 0.33* 1 Pb 0.21 0.52** 0.04 −0.17 0.22 0.17 0.42** −0.17 1 Zn −0.06 0.38* 0.48** 0.28 0.32* 0.12 0.62** 0.27 0.43** 1 注：“*”表示p<0.05, “**”表示p<0.01。 3.4.2 基于主成分分析（PCA）的重金属来源分析
为了进一步明确各土壤中重金属的来源，采用主成分分析法对研究区内的重金属含量进行解析。主成分分析是在保留原始数据大部分信息的情况下，把一些信息重叠、具有错综复杂关系的可观测变量归结为少数具有代表性且不具有相关性的潜在因子，对重金属含量进行提取和汇总，从而可以消除评价指标间相关影响并减少指标工作量的选择。环境学领域中经常利用这些潜在因子来研究土壤重金属的来源[37]。研究区重金属含量的巴特利特球形度（Bartlett球形度）检验结果（p=0.00<0.05）与KMO度量值检验结果（0.66>0.5）表明各个元素原始变量间存在相关性且相关程度无太大差异，数据适合进行主成分分析（表7）。对Kaiser标准化的因子进行Varimax正交旋转后，得到4个特征值大于1的主成分，贡献率分别为33.91%、22.85%、16.05%、12.18%，累积贡献率为84.98%，可解释所有重金属的大部分信息。
表 7 土壤重金属主成分分析矩阵Table 7. Principal component analysis matrix of heavy metals in the soils.元素 主成分 F1 F2 F3 F4 As 0.14 0.28 0.88 0.28 Cd 0.26 0.67 −0.06 0.44 Co 0.95 −0.18 −0.03 −0.02 Cr 0.88 −0.45 0.00 0.01 Cu 0.53 0.24 0.78 0.00 Hg 0.15 0.29 −0.48 0.68 Mn 0.64 0.46 −0.33 −0.03 Ni 0.88 −0.44 −0.02 0.07 Pb 0.16 0.79 0.01 −0.02 Zn 0.60 0.52 −0.20 −0.39 初始特征值 4.07 2.74 1.93 1.46 方差贡献率（%） 33.91 22.85 16.05 12.18 累积贡献率（%） 33.91 56.76 72.80 84.98 第一主成分（F1）贡献率（33.91%）远高于其他主成分，载荷较高的重金属为Co、Cr、Cu、Mn、Ni、Zn，分别为0.95、0.88、0.53、0.64、0.88、0.60。由表1可知，这几种元素的平均含量与青海省的背景值、国家土壤背景值较接近，变异系数较小，受人为影响较小。相关性分析中也显示几种元素间有较好的相关性（表7）。大量研究表明，Ni、Cr、Cu、Mn等在土壤中主要源于地壳，在风化活动较小的区域内，含量与其在成土母质中含量相近，与成岩成分有关，主要受自然影响较大[38]，人为影响较小，因此F1主要受以岩石风化、成图母质等为主的自然因素的影响。
第二主成分（F2）贡献率为22.85%，其中Cd、Pb、Zn具有较大载荷。由相关分析（表7）可知，Cd、Pb、Zn两两之间有相关关系，从空间分布关系来看，三者空间具有相关性。其中，Pb、Zn、Cd在交通运输中的燃料不完全燃烧和轮胎的磨损都将会释放出来[39]。Pb常作为来自机动车污染的标识性元素，Zn常作为交通污染源的标识元素，被广泛用作清洁剂和抗氧化剂，轮胎与地面的摩擦会产生含Zn的颗粒[40]。而车辆润滑油氧化后生成的有机化合物能够腐蚀含Cd、Zn等元素的金属部件及油泵，导致其向环境中释放重金属元素[41]。在青藏高原公路沿线表层土壤中的Cd、Pb、Zn含量较高，且相关性良好，根据研究区内三者之间的相关性判断主要是受交通源的影响。因此F2主要表示交通影响因素。Zn在第一和第二主成分均有一定载荷，表明其可能存在双重来源。
第三主成分（F3）贡献率为16.05%，其中As和Cu的载荷较大。As在研究区内表现出的污染趋势最为明显，且变异系数最大，As受到人为活动的影响最大。研究表明，As进入土壤主要是通过农药、灌溉、煤炭燃烧、有色金属冶炼等方式[28]，Cu的富集来源主要是金属冶炼、交通尾气排放等[6]。根据相关性分析可知，Cu与As之间有正相关性，空间分布特征具有相似性，高值区集中于泽库县北部。经调查，在泽库县北部区域，分布着该县有名的金矿开发地——瓦勒根金矿，瓦勒根金矿位置与As和Cu含量高值区高度重合，表明金矿的开发利用引起了该地土壤中As和Cu的富集。因此，认为F3主要为本地金属冶炼的工业源。Cu在第一和第三主成分均有一定载荷，表明其可能受双重来源的影响。
第四主成分（F4）贡献率为12.18%，载荷较大的重金属只有Hg。从相关性分析和载荷距离来看，Hg与其他元素距离较远，表明Hg来源有别于其他元素。从空间分布看，Hg的平均含量基本上与当地背景值相近，个别样点表现出轻微的富集现象。有研究显示，Hg主要通过大气干沉降和湿沉降进入土壤中[42]；同时煤炭燃烧和钢铁冶炼也是Hg的主要来源[43]。从泽库县的土地利用方式、产业结构等调研来看，在泽库东北部与同仁县扎毛乡接壤处是老藏沟铅锌矿，其位置与Hg最高值位置一致。另外，泽库的冬季取暖等基本上还是以燃煤为主，交通尾气、燃煤燃烧等释放的Hg通过大气沉降等进入到区域表层土壤中，但是由于大气传输的具有局里稀释作用，对于土壤中Hg含量增幅的影响有限，所以整体上从研究区来看，Hg平均含量与当地背景值基本接近。综上，F4归为“远源大气传输”。
3.4.3 基于APCS-MLR模型的重金属来源解析
利用主成分分析识别污染源后，在获得4个污染源因子数量后，对各示踪元素的浓度进行多元线性回归（MLR），得到重金属浓度与各元素浓度的回归式，回归系数用于计算各源的相对贡献率。结算结果列于表8，其中R2为拟合的相关性，ET为重金属的估计值，OT为重金属的实测值，E/O接近1说明拟合度越好，即APCS-MLR分析的可靠性高。
表 8 泽库县表层土壤重金属污染源贡献分析结果Table 8. Analysis of contribution of heavy metal pollution sources in the soils.重金属
元素贡献率（%） ET OT E/O
（%）R2 F1 F2 F3 F4 其他源 As 9.38 4.08 43.52 0.18 42.85 28.95 28.95 99.997 0.98 Cd 3.79 75.46 2.00 1.10 17.65 0.14 0.14 100.418 0.85 Co 48.35 36.91 0.57 0.16 14.01 13.85 13.85 99.986 0.97 Cr 64.49 2.36 0.36 0.09 32.70 63.14 63.15 99.983 0.99 Cu 17.51 5.51 37.29 0.22 39.47 23.84 23.84 99.996 0.97 Hg 1.53 1.34 2.06 49.39 45.69 0.02 0.02 100.00 0.90 Mn 67.68 13.70 0.36 0.20 18.07 774.13 774.15 99.998 0.86 Ni 77.99 0.94 0.43 0.05 20.59 30.63 30.65 99.931 0.99 Pb 5.51 50.75 0.62 0.15 42.97 23.17 23.21 99.828 0.81 Zn 18.07 55.54 0.14 0.03 26.22 74.95 74.96 99.993 0.91 重金属As等的多元回归拟合的误差均较小，从模型的拟合参数的相关系数来看，相关性均较好（模型相关性R2均大于0.8），说明多元线性回归的拟合度好且准确性高。从贡献率分析结果（表8），F1对Cr、Co、Mn、Ni的影响最大，贡献率分别为64.49%、48.35%、67.68%、77.99%，结合主成分分析结果以及元素间的地球化学特征，自然因素是这几种元素的主要来源；F2作为交通源因素对Cd和Pb、Zn的影响最大，贡献率分别为75.46%、50.75%、55.54%；F3为本地冶炼因素对As和Cu的影响最大，贡献率分别为43.53%和37.29%；而F4的大气远源因素对Hg的影响贡献率达到了49.39%。其他源对As、Cr、Hg、Cu、Pb的贡献也较大，分别为42.85%、32.7%、45.69%、39.47%、42.97%，表明这几种重金属还有其他来源影响到土壤中的含量，还需要结合重金属同位素示踪等方法进一步研究分析。从APCS-MLR分析结果来看，与主成分分析的结果基本保持一致，但是根据受体模型的结果显示其他源对As、Cr、Hg、Cu、Pb的贡献也较大，在主成分分析时通常选特征值大于1的因子作为主成分因子，忽略了其他重要但是特征值未大于1的因子，造成了其他源的忽略而未能识别出来。同时需要注意的是，表层土壤中的元素是多种来源的集合，针对单一某种金属通常受多元影响[20-21,44]。
4. 结论
选择中国典型高山农业区域——青海省泽库县采集表层土壤样品，通过测定10种重金属元素含量，根据数据统计分析结果，结合富集因子、地累积指数、潜在生态危害指数对土壤重金属污染进行评价，运用相关性系数、主成分分析、APCS-MLR模型，定性定量解析泽库县重金属来源。结果表明：研究区表层土壤重金属As含量均值超过了国家背景值和青海省的表层土壤背景值，变异系数最大，受到人为影响较大。Cd、Pb、Co、Cr、Cu、Ni、Zn、Hg则表现为低于青海省背景值或与之持平，时空分散程度均较小，受人为影响比较小。Mn平均含量明显超过了背景值，空间分散度不大, 整体分布均匀。重金属含量变化在空间分布上有差异，整体上北部、东北部区域偏高，其他区域偏低。从生态危害程度来看，EF和Igeo结果显示，As的富集程度最高，污染最严重，后续要持续关注As的污染程度变化；其他重金属在个别点位上存在污染，全域基本上处于中度污染状态及以下。从潜在生态危害（PESI）分析来看，单个重金属元素潜在生态风险系数，As的危害程度最大，有重度及中度危害点位，其他重金属则均处于中度及以下危害，其结果与富集系数、地累积指数结果基本一致；重金属整体生态危害系数显示，42件样品点位处于轻度危害内，1件样品点位处于中度危害，与As的最高生态风险系数点位重合。与人为活动频繁、经济发达等区域相比，泽库县这种典型的牧业为主的高山农业区域，大部分重金属含量在表层土壤中偏低。泽库县农业区域土壤重金属来源主要有4类，自然源、交通源、矿业等工业源和大气沉降的远源影响等，同时还有其他源对于As、Cr、Hg、Cu、Pb的影响也较大，后续要持续加强对于其他源的解析及确定。
针对研究结果，后期有必要持续关注As元素的生态危害程度，有效地预防其污染程度加重，还需进一步研究重金属的其他来源及成土母质、母岩的影响，以及关注人类活动对于重金属地球化学特征的影响，更科学合理地解释高山农业土壤的重金属来源，保证土地的安全利用。
-
![]()
图 1 （a）研究区地理位置图；（b）区域点位分布图
Figure 1. （a）Geographic map of the study area；（b）Distribution map of the regional points.
![]()
图 3 研究区土壤（a）TN、（b）TP、（c）K、（d）SOM空间分布图
Figure 3. Spatial distribution maps of soil （a）TN, （b）TP, （c）K and （d） organic matter in the study area.
![]()
图 4 研究区土壤养分地球化学等级图
Figure 4. Geochemical grade map of soil nutrients in the study area.
表 1 样品测试指标与对应测试方法
Table 1 Analysis parameters and methods.
元素/指标 测试仪器/测试方法 要求检出限 本实验检出限 测定范围 pH pH计（ISE） 0.1 0.1 0.1~14.0 N 凯氏氮蒸馏酸碱滴定容量法（VOL） 20mg/kg 15mg/kg 15~5000mg/kg P X射线荧光光谱法 10mg/kg 8mg/kg 8~4500mg/kg K 电感耦合等离子体发射光谱法（ICP-OES） 5.0×10−4 5.0×10−4 2.0×10−4~7.0×10−2 SOM 硫酸亚铁铵容量法（VOL） 1.0×10−3 5.0×10−4 5.0×10−4~1.0 
下载: 导出CSV
表 2 土壤养分和pH值的分级划分标准
Table 2 Standards for classification of the soil nutrient and pH.
养分等级 丰缺程度 TN含量（mg/kg） TP含量（mg/kg） K含量（mg/kg） SOM含量（%） f综合 酸碱等级 pH Ⅰ级 丰富 ≥2000 ≥1000 ≥25 ≥4 ≥4.5 强酸性 <5.0 Ⅱ级 较丰富 1500~2000 800~1000 20~25 3~4 4.5~3.5 酸性 5.0~6.5 Ⅲ级 中等 1000~1500 600~800 15~20 2~3 3.5~2.5 中性 6.5~7.5 Ⅳ级 较缺乏 750~1000 400~600 10~15 1~2 2.5~1.5 碱性 7.5~8.5 Ⅴ级 缺乏 ≤750 ≤400 ≤10 ≤1 ≤1.5 强碱性 >8.5
下载: 导出CSV
表 3 研究区土壤养分指标统计特征
Table 3 Statistics of soil nutrients in the study area.
养分 样本量
（件）含量范围 平均值 中位数 标准差 变异系数
（%）样品分布范围占比（%） 丰富 较丰富 中等 较缺乏 缺乏 TN 1295 222.60~3195.42 1182.94 1149.23 335.89 28.39 2.03 14.00 49.6 28.42 5.94 TP 1295 267.10~2249.70 633.72 584.20 272.87 43.06 10.40 14.78 23.23 26.60 24.98 K 1295 10.24~42.86 20.41 19.66 4.61 22.59 12.83 35.13 49.85 2.19 0 SOM 1295 0.35~3.57 1.19 1.16 0.38 31.47 1.79 3.05 45.34 49.11 0.70 注：TN和TP单位为mg/kg；K单位为mg/g；SOM单位为%。
下载: 导出CSV
表 4 不同土壤类型养分含量（平均值）对照
Table 4 Comparison of nutrient contents (average value) of different soil types.
土壤类型 样本数量（件） TN TP K SOM pH 平均值（mg/kg） 变异系数（%） 平均值（mg/kg） 变异系数（%） 平均值（mg/kg） 变异系数（%） 平均值（mg/kg） 变异系数（%） 平均值（mg/kg） 变异系数（%） 粗骨土 480 1039 24.90 697 34.86 21.50 21.01 1.05 26.89 5.35 5.48 水稻土 387 1344 19.22 432 37.62 17.63 16.45 1.32 20.93 5.37 8.18 棕壤 199 1176 34.77 789 24.66 23.70 23.06 1.24 39.02 5.38 5.67 黄棕壤 111 1183 31.57 891 30.08 21.80 21.83 1.36 38.29 5.36 5.39 紫色土 65 1390 14.04 419 34.44 17.26 14.49 1.37 14.65 5.31 3.00 潮土 53 1060 29.77 676 41.05 19.42 15.76 1.02 32.41 5.41 10.43
下载: 导出CSV
表 5 研究区不同土地利用方式下养分含量（平均值）对照
Table 5 Soil factor data statistics under different land utilization types in the study area.
土地利用方式 样本数量（件） TN TP K SOM pH 平均值（mg/kg） 变异系数（%） 平均值（mg/kg） 变异系数（%） 平均值（mg/kg） 变异系数（%） 平均值（mg/kg） 变异系数（%） 平均值（mg/kg） 变异系数（%） 旱地 262 972.19 0.26 804.51 0.30 22.07 0.21 0.98 0.30 5.34 0.07 林地 133 1340.41 0.41 745.75 0.29 23.80 0.22 1.43 0.47 5.48 0.11 水田（水浇地） 684 1286.48 0.21 530.93 0.50 18.84 0.20 1.28 0.23 5.41 0.06 园地 216 1007.29 0.23 677.65 0.29 21.43 0.23 1.03 0.24 5.28 0.06
下载: 导出CSV
表 6 不同地貌类型养分含量（平均值）对照
Table 6 Soil factor data statistics under different geomorphic types in the study area.
地貌类型 样本数量（件） TN TP K SOM pH 平均值 变异系数（%） 平均值（mg/kg） 变异系数（%） 平均值（mg/kg） 变异系数（%） 平均值（mg/kg） 变异系数（%） 平均值（mg/kg） 变异系数（%） 盆地 249 972.19 0.26 804.51 0.30 22.07 0.21 0.98 0.30 5.34 0.07 平原 216 1007.29 0.23 677.65 0.29 21.43 0.23 1.03 0.24 5.28 0.06 丘陵 684 1286.48 0.21 530.93 0.50 18.84 0.20 1.28 0.23 5.41 0.06 山地 133 1340.41 0.41 745.75 0.29 23.80 0.22 1.43 0.47 5.48 0.11
下载: 导出CSV
表 7 研究区各研究指标及影响因子的皮尔森相关性分析
Table 7 Pearson correlation analysis of each research index and influence factor in the study area.
指标 相关性 土壤类型 土地利用 地貌 TN TP K SOM TN 皮尔森相关 0.13** 0.10** −0.28 显著性（双尾） 0.00 0.00 0.00 TP 皮尔森相关 −0.01 −0.01 0.54** −0.3** 显著性（双尾） 0.70 0.63 0.00 0.00 K 皮尔森相关 −0.15 0.11** 0.40** −0.15** 0.25** 显著性（双尾） 0.00 0.00 0.00 0.00 0.00 SOM 皮尔森相关 0.13** 0.14** −0.16** 0.95** −0.21** −0.09** 显著性（双尾） 0.00 0.00 0.00 0.00 0.00 0.00 pH 皮尔森相关 −0.02 0.02 −0.16** −0.11** −0.06* −0.17** −0.14** 显著性（双尾） 0.44 0.39 0.00 0.00 0.04 0.00 0.00 注：“**”表示相关性在 0.01 层上显著（双尾）；“*”表示相关性在 0.05层上显著（双尾）。
下载: 导出CSV
-
[1] 邓小华, 张瑶, 田峰, 等. 湘西植烟土壤pH和主要养分特征及其相互关系[J]. 土壤, 2017, 49(1): 49−56. Deng X H, Zhang Y, Tian F, et al. pH and main nutrients of tobacco-growing soils and their relations in Western Hunan[J]. Soils, 2017, 49(1): 49−56.
[2] 张子璐, 刘峰, 侯庭钰. 我国稻田氮磷流失现状及影响因素研究进展[J]. 应用生态学报, 2019, 30(10): 3292−3302. doi: 10.13287/j.1001-9332.201910.029 Zhang Z L, Liu F, Hou T Y. Current status of nitrogen and phosphorus losses and related factors in Chinese paddy fields[J]. Chinese Journal of Applied Ecology, 2019, 30(10): 3292−3302. doi: 10.13287/j.1001-9332.201910.029
[3] 陆保国. 霍山县耕地土壤养分现状及施肥对策[J]. 安徽农学通报, 2013, 19(17): 73−75. doi: 10.16377/j.cnki.issn1007-7731.2013.17.030 Lu B G. Present situation of soil nutrients and fertilization countermeasures in Huoshan County[J]. Anhui Agricultural Science Bulletin, 2013, 19(17): 73−75. doi: 10.16377/j.cnki.issn1007-7731.2013.17.030
[4] 唐贤, 梁丰, 徐明岗, 等. 长期施用化肥对农田土壤pH影响的整合分析[J]. 吉林农业大学学报, 2020, 42(3): 316−321. doi: 10.13327/j.jjlau.2020.4506 Tang X, Liang F, Xu M G, et al. A meta-analysis of effects of long-term application of chemical fertilizer on pH of farmland soil[J]. Journal of Jilin Agricultural University, 2020, 42(3): 316−321. doi: 10.13327/j.jjlau.2020.4506
[5] 解庆锋, 周小果, 马振波, 等. 南阳盆地农耕区表层土壤酸碱度现状研究[J]. 安徽农业科学, 2021, 49(5): 69−71, 75. doi: 10.3969/j.issn.0517-6611.2021.05.019 Xie Q F, Zhou X G, Ma Z B, et al. Study on the status of surface soil acidity and alkalinity in agricultural areas of Nanyang Basin[J]. Anhui Agriculture Science, 2021, 49(5): 69−71, 75. doi: 10.3969/j.issn.0517-6611.2021.05.019
[6] Eich-Greatorex S, Sogn T A, Øgaard A F, et al. Plant availability of inorganic and organic selenium fertiliser as influenced by soil organic matter content and pH[J]. Nutrient Cycling in Agroecosystems, 2007, 79(3): 221−231. doi: 10.1007/s10705-007-9109-3
[7] 梁俊捷.气候、地形和土壤性质对我国农田土壤有效钾分布的影响[D].成都: 四川农业大学, 2018. Liang J J. Effects of climate, topography and soil properties on soil available potassium across China[D]. Chengdu: Sichuan Agricultural University, 2018.
[8] 蒋文惠.地形和土地利用对山区土壤养分空间变异的影响[D].泰安: 山东农业大学, 2014. Jiang W H. Effects of land use and topographic factors on soil nutrients variability in mountain area[D]. Tai’an: Shandong Agricultural University, 2014.
[9] 王剑, 周琳, 李志伟, 等. 凉山主要植烟土壤类型氮、磷、钾分布特征差异分析[J]. 山西农业科学, 2012, 40(11): 1191−1194. doi: 10.3969/j.issn.1002-2481.2012.11.19 Wang J, Zhou L, Li Z W, et al. Analysis of nitrogen and phosphorus potassium distribution characteristics variance for the main tobacco soil type profile in Liangshan[J]. Journal of Shanxi Agricultural Sciences, 2012, 40(11): 1191−1194. doi: 10.3969/j.issn.1002-2481.2012.11.19
[10] 张晗, 赵小敏, 欧阳真程, 等. 江西省不同农田利用方式对土壤养分状况的影响[J]. 水土保持研究, 2018, 25(6): 53−60. doi: 10.13869/j.cnki.rswc.2018.06.008 Zhang H, Zhao X M, Ouyang Z C, et al. Effectors of different farmland use types on soil nutrients in Jiangxi Province[J]. Research of Soil and Water Conservation, 2018, 25(6): 53−60. doi: 10.13869/j.cnki.rswc.2018.06.008
[11] 朱永青, 崔云霞, 李伟迪, 等. 太滆运河流域不同用地方式下土壤pH值、有机质及氮磷含量特征分析[J]. 生态与农村环境学报, 2020, 36(2): 171−178. doi: 10.19741/j.issn.1673-4831.2019.0721 Zhu Y Q, Cui Y X, Li W D, et al. Analysis of soil pH, organic matter, nitrogen and phosphate characteristics under different land use types in Taige Canal Valley[J]. Journal of Ecology and Rural Environment, 2020, 36(2): 171−178. doi: 10.19741/j.issn.1673-4831.2019.0721
[12] 刁二龙, 曹广超, 曹生奎, 等. 祁连山南坡不同土地利用方式下土壤理化性质及空间变异性分析[J]. 西南农业学报, 2019, 32(8): 1864−1871. doi: 10.16213/j.cnki.scjas.2019.8.028 Diao E L, Cao G C, Cao S K, et al. Analysis of soil physical and chemical properties and spatial variability under different land use patterns in southern slope of Qilian Mountains[J]. Southwest China Journal of Agricultural Sciences, 2019, 32(8): 1864−1871. doi: 10.16213/j.cnki.scjas.2019.8.028
[13] 王亚男, 徐梦洁, 代圆凤, 等. 毕节市耕地土壤pH的空间变异特征与影响因素[J]. 土壤, 2018, 50(2): 385−390. doi: 10.13758/j.cnki.tr.2018.02.023 Wang Y N, Xu M J, Dai Y F, et al. Spatial characteristics and influential factors of arable soil pH in Bijie, Guizhou[J]. Soils, 2018, 50(2): 385−390. doi: 10.13758/j.cnki.tr.2018.02.023
[14] 王祚林, 陈莉. 六安市农业现代化发展水平评价[J]. 合作经济与科技, 2019(2): 36−38. doi: 10.3969/j.issn.1672-190X.2019.02.013 Wang Z L, Chen L. Evaluation of agricultural modernization development level in Lu’an City[J]. Cooperative Economy and Technology, 2019(2): 36−38. doi: 10.3969/j.issn.1672-190X.2019.02.013
[15] 张孝存, 崔格冲, 杨丽, 等. 商洛金凤山地区人工林土壤pH值及养分特征分析[J]. 辽宁农业科学, 2020(3): 9−12. doi: 10.3969/j.issn.1002-1728.2020.03.003 Zhang X C, Cui G C, Yang L, et al. Analysis on soil nutrient properties and pH of different forest land in Jinfeng Mountain in Shangluo[J]. Liaoning Agricultural Sciences, 2020(3): 9−12. doi: 10.3969/j.issn.1002-1728.2020.03.003
[16] 周子方, 解燕, 易克, 等. 马龙植烟土壤pH值分布特征及其主控因素研究[J]. 中国土壤与肥料, 2019(5): 8−13. Zhou Z F, Xie Y, Yi K, et al. Distribution characteristics and influencing factors of pH value in tobacco growing soil in Malong[J]. Soil and Fertilizer Sciences in China, 2019(5): 8−13.
[17] 张驭航, 李玲, 王秀丽, 等. 河南省土壤pH值时空变化特征分析[J]. 土壤通报, 2019, 50(5): 1091−1100. doi: 10.19336/j.cnki.trtb.2019.05.12 Zhang Y H, Li L, Wang X L, et al. Temporal and spatial variation of soil pH in Henan[J]. Chinese Journal of Soil Science, 2019, 50(5): 1091−1100. doi: 10.19336/j.cnki.trtb.2019.05.12
[18] 向娇, 宋超, 石迎春, 等. 安徽省六安地区土壤pH空间变异特征及其影响因素[J]. 土壤通报, 2021, 52(1): 34−41. Xiang J, Song C, Shi Y C, et al. Spatial variation characteristics and influencing factors of soil pH in the Lu ’an area of Anhui Province[J]. Chinese Journal of Soil Science, 2021, 52(1): 34−41.
[19] 闫志浩, 胡志华, 王士超, 等. 石灰用量对水稻油菜轮作区土壤酸度、土壤养分及作物生长的影响[J]. 中国农业科学, 2019, 52(23): 4285−4295. doi: 10.3864/j.issn.0578-1752.2019.23.009 Yan Z H, Hu Z H, Wang S C, et al. Effects of lime content on soil acidity, soil nutrients and crop growth in rice-rape rotation system[J]. Scientia Agricultura Sinica, 2019, 52(23): 4285−4295. doi: 10.3864/j.issn.0578-1752.2019.23.009
[20] 姬王佳.陕北黄土区深剖面土壤水和氮素来源及其对土地利用变化的响应[D].杨凌: 西北农林科技大学, 2021. Ji W J. Sources of soil water and nitrogen and their response to land use change in deep deposits of Northern Shaanxi, China[D]. Yangling: Northwest A&F University, 2021.
[21] 张洋洋.白龙江流域土壤碳氮磷含量及空间分布特征研究[D].兰州: 兰州大学, 2017. Zhang Y Y. Study on contents and spatial distribution characteristics of soil carbon, nitrogen and phosphorus in Bailing River Basin[D]. Lanzhou: Lanzhou University, 2017.
[22] 李龙, 姚云峰, 秦富仓. 黄花甸子流域土壤全氮、速效磷、速效钾的空间变异[J]. 生态学杂志, 2015, 34(2): 373−379. doi: 10.13292/j.1000-4890.2015.0052 Li L, Yao Y F, Qin F C. Spatial variability of soil total nitrogen, available phosphorus and available potassium in Huanghuadianzi Watershed[J]. Chinese Journal of Ecology, 2015, 34(2): 373−379. doi: 10.13292/j.1000-4890.2015.0052
[23] 顾涛, 朱晓华, 赵信文, 等. 广州新垦莲藕产区莲藕品质与地球化学条件的关系[J]. 岩矿测试, 2021, 40(6): 833−845. doi: 10.3969/j.issn.0254-5357.2021.6.ykcs202106004 Gu T, Zhu X H, Zhao X W, et al. Relationship between lotus root quality and geochemical conditions in the Xinken lotus root producing area of Guangzhou[J]. Rock and Mineral Analysis, 2021, 40(6): 833−845. doi: 10.3969/j.issn.0254-5357.2021.6.ykcs202106004
[24] 许春雪, 袁建, 王亚平, 等. 沉积物中磷的赋存形态及磷形态顺序提取分析方法[J]. 岩矿测试, 2011, 30(6): 785−794. doi: 10.3969/j.issn.0254-5357.2011.06.024 Xu C X, Yuan J, Wang Y P, et al. Speciation and release mechanism of phosphorus in sediments and analysis methods for sequential extraction[J]. Rock and Mineral Analysis, 2011, 30(6): 785−794. doi: 10.3969/j.issn.0254-5357.2011.06.024
[25] 董秋瑶, 赖书雅, 宋超, 等. 南阳盆地东部山区土壤锗分布特征及其影响因素分析[J]. 环境科学, 2022, 43(6): 442−452. Dong Q Y, Lai S Y, Song C, et al. Distribution characteristics and influencing factors of germanium in soil in the eastern mountainous area of the Nanyang Basin[J]. Environmental Science, 2022, 43(6): 442−452.
[26] 王辉.皖西大别山区土壤磷、硒分布特征及影响因素[D].北京: 中国地质大学（北京）, 2018. Wang H. The distribution of phosphorus and selenium and their controlling factors of soils in Dabie Mountain, West Anhui Province[D]. Beijing: China University of Geosciences （Beijing）, 2018.
[27] 夏炎, 宋延斌, 侯进凯, 等. 河南洛阳市土壤和农作物中钼分布规律与影响因素研究[J]. 岩矿测试, 2021, 40(6): 820−832. doi: 10.3969/j.issn.0254-5357.2021.6.ykcs202106003 Xia Y, Song Y B, Hou J K, et al. Distribution law and influencing factors of molybdenum in soils and crops in Luoyang, Henan Province[J]. Rock and Mineral Analysis, 2021, 40(6): 820−832. doi: 10.3969/j.issn.0254-5357.2021.6.ykcs202106003
[28] 白娜, 王立, 孔东升. 黑河自然保护区沼泽湿地土壤化学性质的空间分布特征研究[J]. 草业学报, 2017, 26(5): 15−28. doi: 10.11686/cyxb2016269 Bai N, Wang L, Kong D S. Spatial distribution and chemical properties of marsh wetland soil in the Heihe Nature Reserve[J]. Acta Prataculture Sinica, 2017, 26(5): 15−28. doi: 10.11686/cyxb2016269
[29] 郭安廷, 崔锦霞, 许鑫, 等. 基于GIS与地统计的土壤养分空间变异研究[J]. 中国农学通报, 2018, 34(23): 72−79. doi: 10.11924/j.issn.1000-6850.casb17070049 Guo A T, Cui J X, Xu X, et al. Spatial distribution of soil nutrients based on GIS and geostatistics[J]. Chinese Agricultural Science Bulletin, 2018, 34(23): 72−79. doi: 10.11924/j.issn.1000-6850.casb17070049
[30] 梁斌, 齐实. 北京山区土壤养分空间变化特征研究[J]. 土壤, 2018, 50(4): 769−777. doi: 10.13758/j.cnki.tr.2018.04.017 Liang B, Qi S. Spatial distributions of soil nutrients in Beijing mountainous area[J]. Soils, 2018, 50(4): 769−777. doi: 10.13758/j.cnki.tr.2018.04.017
[31] 杨东伟.水田改旱作后土壤性态与土壤类型演变研究[D].杭州: 浙江大学, 2014. Yang D W. Evolution of soil characteristics and soil types after land-use conversion from paddy field to upland[D]. Hangzhou: Zhejiang University, 2014.
[32] Dong Q Y, Xiang J, Song C, et al. Spatial distribution characteristics and main controlling factors of germanium in soil of Northern Dabie Mountains, China[J]. China. Journal of Groundwater Science and Engineering, 2022, 10(4): 381−392.
[33] 周殷竹, 王彪, 刘义, 等. 青海囊谦县城周边农耕区土壤质量地球化学评价及富硒土地利用分区[J]. 干旱区资源与环境, 2020, 34(10): 93−101. doi: 10.13448/j.cnki.jalre.2020.274 Zhou Y Z, Wang B, Liu Y, et al. Geochemical evaluation of soil quality and land-use regionalization of selenium-rich soils in cultivated area around Nangqian County, Qinghai[J]. Journal of Arid Land Resources and Environment, 2020, 34(10): 93−101. doi: 10.13448/j.cnki.jalre.2020.274
[34] 杨红, 徐唱唱, 赛曼, 等. 不同土地利用方式对土壤含水量、pH值及电导率的影响[J]. 浙江农业学报, 2016, 28(11): 1922−1927. doi: 10.3969/j.issn.1004-1524.2016.11.18 Yang H, Xu C C, Sai M, et al. Effects of land use on soil moisture, pH and electrical conductivity[J]. Acta Agriculturae Zhejiangensis, 2016, 28(11): 1922−1927. doi: 10.3969/j.issn.1004-1524.2016.11.18
[35] 董秋瑶, 温皓天, 宋超, 等. 河南南阳盆地东南耕区表层土壤养分地球化学综合评价及影响因素[J]. 现代地质, 2022, 36(2): 449−461. Dong Q Y, Wen H T, Song C, et al. Comprehensive evaluation and influencing factors of surface soil nutrient chemistry in southeastern cultivated area of the Nanyang, Henan Province[J]. Geoscience, 2022, 36(2): 449−461.
[36] Wang Y, Zhang X, Huang C. Spatial variability of soil total nitrogen and soil total phosphorus under different land uses in a small watershed on the Loess Plateau, China[J]. Geoderma, 2009, 150(1): 141−149.
[37] 张阳, 屠乃美, 谢会雅, 等. 茶陵县烟稻复种区土壤养分时空变异特征[J]. 土壤通报, 2021, 52(2): 297−305. doi: 10.19336/j.cnki.trtb.2020081201 Zhang Y, Tu N M, Xie H Y, et al. Spatial and temporal variation of soil nutrients in tobacco-rice multiple-cropping area of Chaling[J]. Chinese Journal of Soil Science, 2021, 52(2): 297−305. doi: 10.19336/j.cnki.trtb.2020081201
-
期刊类型引用(9)
1. 黄春霞，黄立幸. 广西废弃矿区土壤中Cd、Pb、As污染特征及环境质量评价研究. 矿产与地质. 2024(01): 168-173 . 
百度学术
2. 周李蕾. 土壤重金属污染的环境影响评价及管理分析. 四川有色金属. 2024(01): 54-57 . 百度学术
3. 张妍，赵新雷，冯雪珍，郭亚娇. 河南荥阳市耕地土壤重金属分布特征及来源解析. 岩矿测试. 2024(02): 330-343 . 本站查看
4. 南维鸽，董治宝，薛亮，张智，哈延庆，邵天杰. 青藏高原重要交通国道路侧土壤重金属分布特征及生态风险评价. 环境科学. 2024(08): 4825-4836 . 百度学术
5. 谢勇，黄秋雨. 地质样品中重金属污染的来源与分布特征. 世界有色金属. 2024(12): 229-231 . 百度学术
6. 韩善博，丑亚玲，曹伟，刘长东，郭焕明，彭尔兴. 寒区土壤重金属迁移转化与修复治理研究进展. 冰川冻土. 2024(04): 1369-1385 . 百度学术
7. 王立美，张雷雷. 土壤重金属污染对环境的影响与防治措施. 农村科学实验. 2024(19): 48-50 . 百度学术
8. 徐苏霞. 土壤重金属污染对农业种植的危害及治理策略. 黑龙江粮食. 2024(09): 119-121 . 百度学术
9. 李斯林，王贺封，刘佳，张安兵，魏凯濠，李斌. 冀南非露天煤矿土壤重金属风险评价与影响因素分析. 岩矿测试. 2024(05): 769-782 . 本站查看
其他类型引用(2)
计量
- 文章访问数: 144
- HTML全文浏览量: 25
- PDF下载量: 60
- 被引次数: 11
京公网安备 11010202008159号