• ISSN 1008-505X
  • CN 11-3996/S

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

黑土区县域土壤养分空间分布特征及其影响因子

杜可 王乐 张淑香 雷秋良

引用本文:
Citation:

黑土区县域土壤养分空间分布特征及其影响因子

    作者简介: 杜可 E-mail: 429474807@qq.com;
    通讯作者: 张淑香, E-mail:zhangshuxiang@caas.cn

Spatial distribution characteristics and influence factors of soil nutrients in black soil region counties

    Corresponding author: ZHANG Shu-xiang, E-mail:zhangshuxiang@caas.cn ;
  • 摘要: 【目的】 土壤养分是农业生产的必要条件,了解土壤养分的空间变化,对指导精准施肥有一定意义。 【方法】 本文以地统计学为基础,利用ArcGIS软件,对克山县土壤养分数据进行析取克里格以及泛克里格插值,得到土壤养分的空间分布图,并通过DEM数据提取坡度、坡向、曲率等相关地形因子,在SPSS软件中进行地形因子与土壤养分的相关分析。 【结果】 克山县土壤有机质、有效磷、速效钾和全氮含量平均值分别为43.89 g/kg、28.35 mg/kg、184.75 mg/kg和1.81 g/kg,有机质与全氮的变异系数均为0.09,属于弱变异强度,而有效磷与速效钾的变异系数分别为0.18与0.11,属于中等变异强度。有机质、速效钾及全氮的含量与海拔均呈显著的正相关关系(P<0.01),其中,有机质与全氮受海拔的影响程度最大。有机质以及全氮的含量还受到坡度的影响,与坡度呈显著的正相关性。不同土地利用方式间,旱地和菜地有机质和速效钾含量显著高于水田,而有效磷显著低于水田,全氮含量差异不显著。 【结论】 克山黑土区土壤有机质和全氮含量变异小,而有效磷和速效钾含量变异大。土壤养分的空间分布受到坡度与海拔的影响,尤其海拔对有机质与全氮的影响程度最大。土地利用类型、土壤类型以及成土母质也影响养分的分布。
  • 图 1  采样点分布图

    Figure 1.  Sampling point distribution map

    图 2  养分含量趋势分析

    Figure 2.  Trend analysis of nutrient contents

    图 3  泛克里格法 (左) 和析取克里格法 (右) 表征的土壤养分空间分布特征

    Figure 3.  Spatial distribution characteristics of soil nutrients described by universal (Left) and disjunctive kriging (Right)

    图 4  地形因子空间分布特征

    Figure 4.  Terrain factor characteristics of soil nutrients

    表 1  土壤养分统计特征值

    Table 1.  Statistical characteristics of soil nutrients

    土壤养分
    Soilnutrients
    统计值Statistics 变异系数
    CV
    偏斜度
    Skewness
    峰度
    Kurtosis
    K-S检验
    K-S test
    样本数Sample No. 最小值Min. 最大值Max. 平均数Mean 中值Median
    有机质 SOM(g/kg) 3520 35.00 53.70 43.89 43.703 0.09 0.159 –0.431 4.56E-08
    有效磷 Avail. P(mg/kg) 3545 16.40 41.50 28.35 28.10 0.18 0.178 –0.469 2.40E-07
    速效钾 Avail. K(mg/kg) 3455 131.00 239.00 184.75 186.00 0.11 –0.157 0.316 1.40E-24
    全氮 Total N(g/kg) 3524 1.44 2.20 1.81 1.81 0.09 0.086 –0.578 2.20E-05
    下载: 导出CSV

    表 2  泛克里格和析取克里格插值法相关数值

    Table 2.  Accuracy of parameters used in universal and disjunctive krigings

    项目
    Item
    插值方法
    Kriging
    标准平均值
    Standardized mean
    均方根预测误差
    Root-mean-square
    平均标准误差
    Average mean error
    标准均方根预测误差
    Root-mean-square standardized
    有机质SOM 泛克里格法Universal kriging –0.00262 2.40323 2.31138 1.03715
    析取克里格法Disjuctive kriging 0.00331 2.73138 2.91306 0.93739
    有效磷
    Avail. P
    泛克里格法Universal kriging –0.00536 3.44905 3.54562 0.95637
    析取克里格法Disjuctive kriging –0.00514 3.46199 3.94540 0.87078
    速效钾
    Avail. K
    泛克里格Universal kriging 0.00235 12.83784 13.86553 0.92381
    析取克里格Disjuctive kriging 0.00493 12.92830 15.63747 0.82311
    全氮
    Total N
    泛克里格Universal kriging 0.00186 0.09891 0.10055 0.97619
    析取克里格Disjuctive kriging 0.00211 0.11484 0.12190 0.94208
    下载: 导出CSV

    表 3  土壤养分与环境因子相关性分析

    Table 3.  Correlation analysis between soil nutrients and environmental factors

    项目
    Item
    有机质SOM 有效磷Avail. P 速效钾Avail. K 全氮Total N
    相关性Correlation 显著性Significance 相关性Correlation 显著性Significance 相关性Correlation 显著性Significance 相关性Correlation 显著性Significance
    海拔Latitude 0.391** 0.000 0.029 0.088 0.141** 0.000 0.388** 0.000
    坡度Gradient 0.054** 0.001 –0.020 0.233 0.028 0.104 0.051** 0.002
    坡向Aspect 0.029 0.081 0.001 0.935 0.024 0.157 0.032 0.057
    曲率Curvature 0.008 0.633 –0.017 0.324 –0.015 0.376 0.003 0.881
    注(Note):**表示在 0.01 水平上显著相关 Indicates a significant correlation at the 0.01 level.
    下载: 导出CSV

    表 4  不同土地利用类型的土壤养分样本数 (n) 和含量

    Table 4.  Sample number (n) and soil nutrient contents of different land uses

    土地利用方式
    Land use
    有机质SOM (g/kg) 有效磷Avail. P (mg/kg) 速效钾Avail. K (mg/kg) 全氮Total N (g/kg)
    n 含量Content n 含量Content n 含量Content n 含量Content
    旱地Dry land 3417 43.90 ± 3.82 3439 28.39 ± 5.24 3349 184.60 ± 21.22 3421 1.81 ± 0.16
    水田Paddy 97 43.61 ± 2.91 100 26.88 ± 4.71 100 189.14 ± 12.58 97 1.83 ± 0.12
    菜地Vegetable 6 43.95 ± 0.54 6 26.67 ± 0.67 6 197.50 ± 7.29 6 1.78 ± 0.02
    下载: 导出CSV

    表 5  不同土壤类型的土壤样本数(n)和养分含量

    Table 5.  Sample number (n) and soil nutrient contents of different soil types

    土类
    Soil type
    有机质 SOM (g/kg) 有效磷 Avail. P (mg/kg) 速效钾 Avail. K (mg/kg) 全氮 Total N (g/kg)
    n 含量Content n 含量Content n 含量Content n 含量Content
    黑土Black soil 2173 44.25 ± 3.89 2178 28.49 ± 5.30 2142 184.79 ± 21.37 2181 1.82 ± 0.16
    棕壤Brown soil 298 44.61 ± 3.77 310 28.04 ± 5.02 285 187.48 ± 21.19 303 1.84 ± 0.15
    草甸土Meadow soil 967 43.14 ± 3.43 997 27.96 ± 5.12 960 184.55 ± 19.39 961 1.78 ± 0.15
    黑钙土Chernozems 61 40.80 ± 3.56 59 30.89 ± 4.58 54 171.33 ± 27.36 58 1.66 ± 0.15
    沼泽土Wetland 21 43.35 ± 2.41 21 29.21 ± 4.58 14 189.64 ± 26.16 21 1.78 ± 0.11
    下载: 导出CSV

    表 6  不同成土母质的土壤样本数(n)和养分含量

    Table 6.  Sample number (n) and soil nutrient contents of different soil parent materials

    成土母质
    Parent material
    有机质 SOM (g/kg) 有效磷 Avail. P (mg/kg) 速效钾 Avail. K (mg/kg) 全氮 Total N (g/kg)
    n 含量Content n 含量Content n 含量Content n 含量Content
    黄土母质Loess 2244 44.14 ± 3.92 2243 28.49 ± 5.31 2192 184.23 ± 21.90 2246 1.82 ± 0.16
    冲积物Alluvium 1163 43.36 ± 3.56 1182 28.08 ± 5.12 1158 185.19 ± 19.22 1165 1.79 ± 0.15
    坡积物Slope deposit 75 44.94 ± 3.46 82 28.24 ± 4.97 75 190.63 ± 19.98 76 1.86 ± 0.13
    玄武岩Basalt 14 43.72 ± 2.09 14 27.11 ± 2.11 13 190.46 ± 18.04 13 1.81 ± 0.09
    火山岩Lava 3 43.70 ± 0.40 3 28.87 ± 0.75 3 205.00 ± 3.00 3 1.82 ± 0.02
    淤积物Sediment 21 43.35 ± 2.41 21 29.21 ± 4.58 14 189.64 ± 26.16 21 1.78 ± 0.11
    下载: 导出CSV
  • [1] 董敏. 基于GARBF神经网络的土壤属性信息空间插值方法研究[D]. 成都: 四川农业大学硕士学位论文, 2009.
    Dong M. Research on spatial interpolation method of soil attribute information based on GARBF Neural Network[D]. Chengdu: MS Thesis of Sichuan Agricultural University, 2009
    [2] 王辛芝, 张甘霖, 余元春, 等. 南京城市土壤pH和养分的空间分布[J]. 南京林业大学学报 (自然科学版), 2006, 30(4): 69–72
    Wang X Z, Zhang G L, Yu Y C, et al. Spatial distribution of soil pH and nutrients in urban Nanjing[J]. Journal of Nanjing Forestry University (Natural Science Edition), 2006, 30(4): 69–72
    [3] 王政权. 地统计学及在生态学中的应用[M]. 北京: 科学出版社, 1999.
    Wang Z Q. Geostatistics and its application in ecology[M]. Beijing: Science Press, 1999.
    [4] 张淑娟, 何勇, 方慧. 基于GPS和GIS的田间土壤特性的空间变异性的研究[J]. 农业工程学报, 2003, 19(2): 39–44 doi: 10.3321/j.issn:1007-4333.2003.02.010
    Zhang S J, He Y, Fang H. Spatial variability of soil properties in the field based on GPS and GIS[J]. Transactions of the Chinese Society of Agricultural Engineering, 2003, 19(2): 39–44 doi: 10.3321/j.issn:1007-4333.2003.02.010
    [5] 杨玉玲, 文启凯, 田长彦, 等. 土壤空间变异研究现状及展望[J]. 干旱区研究, 2001, 18(2): 50–55
    Yang Y L, Wen Q K, Tian C Y, et al. Research status and prospect of soil spatial variability[J]. Arid Zone Research, 2001, 18(2): 50–55
    [6] 张文博, 张福平, 苏玉波, 等. 渭河干流沿岸土壤有机质空间分布特征及其影响因素[J]. 水土保持通报, 2014, 34(1): 138–143
    Zhang W B, Zhang F P, Su Y B, et al. Spatial distribution of soil organic matter along the main stream of Weihe River and its influencing factors[J]. Bulletin of Soil and Water Conservation, 2014, 34(1): 138–143
    [7] Smith J L, Halvorson J J, Jr H B. Soil properties and microbial activity across a 500 m elevation gradient in a semi-arid environment[J]. Soil Biology & Biochemistry, 2002, 34(11): 1749–1757
    [8] 张珍明, 林绍霞, 张清海, 等. 不同土地利用方式下草海高原湿地土壤碳、氮、磷分布特征[J]. 水土保持学报, 2013, 27(6): 199–204 doi: 10.3969/j.issn.1009-2242.2013.06.039
    Zhang Z M, Lin S X, Zhang Q H, et al. Distribution characteristics of soil carbon, nitrogen and phosphorus in Caohai Plateau Wetland under different land use patterns[J]. Journal of soll and , 2013, 27(6): 199–204 doi: 10.3969/j.issn.1009-2242.2013.06.039
    [9] Shahbeik S, Afzal P, Moarefvand P, et al. Comparison between ordinary kriging (OK) and inverse distance weighted (IDW) based on estimation error. Case study: Dardevey iron ore deposit, NE Iran[J]. Arabian Journal of Geosciences, 2014, 7(9): 3693–3704 doi: 10.1007/s12517-013-0978-2
    [10] Vano M, Gennai S, Karapetsa D, et al. Comparación de los interpoladores IDW Y Kriging en la variación espacial de pH, Ca, CICE y P del suelo[J]. Agronomía Costarricense, 2008, 32(218): 95–105
    [11] Nedbalek J. Rbf neural networks for function approximation in dynamic modeling[J]. Journal of Konbin, 2008, 8(1): 223–232 doi: 10.2478/v10040-008-0115-6
    [12] 李艳, 史舟, 程街亮, 等. 辅助时序数据用于土壤盐分空间预测及采样研究[J]. 农业工程学报, 2006, 22(6): 19–55 doi: 10.3321/j.issn:1007-4333.2006.06.004
    Li Y, Shi Z, Cheng J L, et al. Spatial forecast and sampling of soil salinity by Kriging with temporally auxiliary data[J]. Transactions of the Chinese Society of Agricultural Engineering, 2006, 22(6): 19–55 doi: 10.3321/j.issn:1007-4333.2006.06.004
    [13] 于洋, 卫伟, 陈利顶, 等. 黄土高原年均降水量空间插值及其方法比较[J]. 应用生态学报, 2015, 26(4): 999–1006
    Yu Y, Wei W, Chen L D, et al. Spatial interpolation of annual precipitation in Loess Plateau and comparison of its methods[J]. Chinese Journal of Applied Ecology, 2015, 26(4): 999–1006
    [14] 付传城, 王文勇, 潘剑君, 等. 南京市溧水区土壤重金属污染不同插值方法的对比研究[J]. 土壤通报, 2014, 45(6): 1325–1333
    Fu C C, Wang W Y, Pan J J, et al. Comparative study on different interpolation methods of soil heavy metal pollution in Lishui District of Nanjing City[J]. Chinese Journal of Soil Science, 2014, 45(6): 1325–1333
    [15] Guan W H, Chamberlain R H, Sabol B M, Doering P H. Mapping submerged aquatic vegetation with GIS in the Caloosahatchee Estuary: evaluation of different interpolation methods[J]. Marine Geodesy, 1999, 22(2): 69–91 doi: 10.1080/014904199273506
    [16] 张铁婵, 常庆瑞, 刘京. 土壤养分元素空间分布不同插值方法研究——以榆林市榆阳区为例[J]. 干旱地区农业研究, 2010, 28(2): 177–182
    Zhang T C, Chang Q R, Liu J. Study on different interpolation methods for spatial distribution of soil nutrient elements: a case study of Yuyang District, Yulin City[J]. Agricultural Research in the Arid Areas, 2010, 28(2): 177–182
    [17] 秦玮, 何方, 范贵娟, 等. 不同插值方法对土壤养分空间插值结果的影响——以寿西湖农场为例[J]. 安徽农学通报, 2017, 23(19): 51–53 doi: 10.3969/j.issn.1007-7731.2017.19.027
    Qin W, He F, Fan G J, et al. Effect of different interpolation methods on spatial interpolation results of soil nutrients: a case study of Shouxi Lake Farm[J]. Anhui Agricultural Science Bulletin, 2017, 23(19): 51–53 doi: 10.3969/j.issn.1007-7731.2017.19.027
    [18] 张慧智, 史学正, 于东升, 等. 中国土壤温度的空间插值方法比较[J]. 地理研究, 2008, 27(6): 1299–1307 doi: 10.3321/j.issn:1000-0585.2008.06.009
    Zhang H Z, Shi X Z, Yu D S, et al. Comparison of spatial interpolation methods for soil temperature in China[J]. Geographical Research, 2008, 27(6): 1299–1307 doi: 10.3321/j.issn:1000-0585.2008.06.009
    [19] 秦耀东, 李保国. 应用析取克里格方法估计区域地下水埋深分布[J]. 水利学报, 1998, 29(8): 28–33 doi: 10.3321/j.issn:0559-9350.1998.08.006
    Qin Y D, Li B G. Estimation of regional groundwater depth distribution by disjunctive Kriging method[J]. Journal of Hydraulic Engineering, 1998, 29(8): 28–33 doi: 10.3321/j.issn:0559-9350.1998.08.006
    [20] 郭奇峰, 刘金辉. 析取克里格法在可回采资源量的总体估算的应用[J]. 中国矿业, 2012, 21(8): 20–22 doi: 10.3969/j.issn.1004-4051.2012.08.006
    Guo Q F, Liu J H. Application of disjunctive Kriging method in the total estimation of recoverable resources[J]. China Mining Magazine, 2012, 21(8): 20–22 doi: 10.3969/j.issn.1004-4051.2012.08.006
    [21] 鲁如坤. 土壤农业化学分析方法[M]. 北京: 中国农业科技出版社, 2000.
    Lu R K. Analytical methods of soil and agronomic chemistry[M]. Beijing: China Agricultural Science and Technology Press, 2000.
    [22] Liu Y, Fan Z P, You T H, et al. Large group decision-making (LGDM) with the participators from multiple subgroups of stakeholders: A method considering both the collective evaluation and the fairness of the alternative[J]. Computers & Industrial Engineering, 2018, 122(1): 262–272
    [23] Shen C, Bao X, Tan J, et al. Two noise-robust axial scanning multi-image phase retrieval algorithms based on Pauta criterion and smoothness constraint[J]. Optics Express, 2017, 25(14): 16235–16249 doi: 10.1364/OE.25.016235
    [24] 姜春雷. 克里格插值的加速和参数优化及其应用[D]. 哈尔滨: 中国科学院研究生院 (东北地理与农业生态研究所) 博士学位论文, 2016.
    Jiang C L. Acceleration and parameter optimization of kriging interpolation and its application[D]. Haerbin: PhD Dissertation of Chinese Academy of Sciences (Northeast Institute of Geography and Agroecology), 2016.
    [25] Li L, Dong J Dong J, et al. Prediction of the spatial distribution of bovine endemic fluorosis using ordinary Kriging[J]. Bulletin of the Veterinary Institute in Pulawy, 2015, 59(1): 161–164 doi: 10.1515/bvip-2015-0024
    [26] Viswanathan R, Jagan J, Samui P, et al. Spatial variability of rock depth using simple Kriging, ordinary Kriging, RVM and MPMR[J]. Geotechnical & Geological Engineering, 2015, 33(1): 69–78
    [27] Carr J R, Roberts K P. Application of universal kriging for estimation of earthquake ground motion: Statistical significance of results[J]. Mathematical Geology, 1989, 21(2): 255–265 doi: 10.1007/BF00893218
    [28] Elsayed K. Optimization of the cyclone separator geometry for minimum pressure drop using Co-Kriging[J]. Powder Technology, 2015, 269(1): 409–424
    [29] Daya A A. Application of disjunctive kriging for estimating economic grade distribution in an iron ore deposit: A case study of the Choghart North Anomaly, Iran[J]. Journal of the Geological Society of India, 2014, 83(5): 567–576 doi: 10.1007/s12594-014-0085-3
    [30] Brown J R, Harvey M E. Rational arithmetic mathematica functions to evaluate the one-sided one-sample K-S cumulative sample distribution[J]. Journal of Statistical Software, 2007, 26(2): 1–40
    [31] Jahufer A. Identifying unusual observations in ridge regression linear model using box-cox power transformation technique[J]. Open Journal of Statistics, 2014, 04(1): 19–26 doi: 10.4236/ojs.2014.41003
    [32] 曾春阳. 基于GIS的小班立地因子提取与立地评价研究[D]. 湖南株洲: 中南林业科技大学硕士学位论文, 2008.
    Zeng C Y. Study on extraction and site evaluation of small class site factors based on GIS[D]. Zhuzhou, Hunan: MS Thesis of Central South University of Forestry and Technology, 2008.
    [33] 杜朝正. 基于ArcGIS的坡度分析[J]. 资源开发与市场, 2009, 25(1): 17–18 doi: 10.3969/j.issn.1005-8141.2009.01.007
    Du C Z. Slope analysis based on ArcGIS[J]. Resource Development & Market, 2009, 25(1): 17–18 doi: 10.3969/j.issn.1005-8141.2009.01.007
    [34] 朱雷, 秦富仓, 苏江. 基于ArcGIS9. 3的等高线生成DEM及坡度坡向分析[J]. 内蒙古林业调查设计, 2014, 37(2): 125–128 doi: 10.3969/j.issn.1006-6993.2014.02.055
    Zhu L, Qin F C, Su J. Contour generating DEM and slope direction analysis based on ArcGIS9. 3[J]. Investigation and Design of Forestry in Inner Mongolia, 2014, 37(2): 125–128 doi: 10.3969/j.issn.1006-6993.2014.02.055
    [35] 高凤杰, 马泉来, 韩文文, 等. 黑土丘陵区小流域土壤有机质空间变异及分布格局[J]. 环境科学, 2016, 37(5): 1915–1922
    Gao F J, Ma Q L, Han W W, et al. Spatial variation and distribution pattern of soil organic matter in small watershed of black soil hilly region[J]. Environmental Science, 2016, 37(5): 1915–1922
    [36] 曹祥会, 龙怀玉, 周脚根, 等. 河北省表层土壤有机碳和全氮空间变异特征性及影响因子分析[J]. 植物营养与肥料学报, 2016, 22(4): 937–948
    Cao X H, Long H Y, Zhou J G, et al. Spatial variability of organic carbon and total nitrogen in surface soil of Hebei Province and its influencing factors[J]. Journal of Plant Nutrition and Fertilizer, 2016, 22(4): 937–948
    [37] 赵艳艳, 周华坤, 姚步青, 等. 长期增温对高寒草甸植物群落和土壤养分的影响[J]. 草地学报, 2015, 23(4): 665–671
    Zhao Y Y, Zhou H K, Yao B Q, et al. Effects of long-term warming on plant communities and soil nutrients in alpine meadow[J]. Acta Agrectir Sinica, 2015, 23(4): 665–671
    [38] 崔贝, 王纪华, 杨武德, 等. 冬小麦-夏玉米轮作区土壤养分时空变化特征[J]. 中国农业科学, 2013, 46(12): 2471–2482 doi: 10.3864/j.issn.0578-1752.2013.12.008
    Cui B, Wang J H, Yang W D, et al. Temporal and spatial variation characteristics of soil nutrients in winter wheat summer maize rotation area[J]. Scientia Agricultura Sinica, 2013, 46(12): 2471–2482 doi: 10.3864/j.issn.0578-1752.2013.12.008
    [39] 孙大生. 水分与温度对土壤碳氮磷转化的影响及其机理[D]. 杭州: 浙江大学博士学位论文, 2016.
    Sun D S. Effect of water and temperature on soil carbon, nitrogen and phosphorus transformation and its mechanism[D]. Hanghou:PhD Dissertation of Zhejiang University, 2016.
    [40] 许中旗, 闵庆文, 王英舜, 等. 人为干扰对典型草原生态系统土壤养分状况的影响[J]. 水土保持学报, 2006, 20(5): 38–42 doi: 10.3321/j.issn:1009-2242.2006.05.010
    Xu Z Q, Min Q W, Wang Y S, et al. Effects of anthropogenic disturbance on soil nutrient status in typical grassland ecosystem[J]. Journal of Soil and Water Conservation, 2006, 20(5): 38–42 doi: 10.3321/j.issn:1009-2242.2006.05.010
    [41] 庞夙, 李廷轩, 王永东, 等. 土壤速效氮、磷、钾含量空间变异特征及其影响因子[J]. 植物营养与肥料学报, 2009, 15(1): 114–120 doi: 10.3321/j.issn:1008-505X.2009.01.016
    Pang S, Li Y X, Wang Y D, et al. Spatial variation characteristics of available nitrogen, phosphorus and potassium in soil and their influencing factors[J]. Plant Nutrition and Fertilizer Science , 2009, 15(1): 114–120 doi: 10.3321/j.issn:1008-505X.2009.01.016
  • [1] 徐新朋赵士诚张云贵何萍高强 . 吉林省玉米种植区土壤养分空间变异特征研究. 植物营养与肥料学报, 2011, 19(6): 1342-1350. doi: 10.11674/zwyf.2011.1044
    [2] 曹祥会龙怀玉周脚根邱卫文雷秋良刘颖李军穆真 . 河北省表层土壤有机碳和全氮空间变异特征性及影响因子分析. 植物营养与肥料学报, 2016, 24(4): 937-948. doi: 10.11674/zwyf.15135
    [3] 庞夙李廷轩王永东余海英吴德勇 . 土壤速效氮、磷、钾含量空间变异特征及其影响因子 . 植物营养与肥料学报, 2009, 17(1): 114-120. doi: 10.11674/zwyf.2009.0116
    [4] 刘永红倪中应谢国雄徐立军钟林炳马立强 . 浙西北丘陵区农田土壤微量元素空间变异特征及影响因子. 植物营养与肥料学报, 2016, 24(6): 1710-1718. doi: 10.11674/zwyf.15343
    [5] 赵业婷李志鹏常庆瑞*齐雁冰刘京 . 西安市粮食主产区耕层土壤速效养分空间特征. 植物营养与肥料学报, 2013, 21(6): 1376-1385. doi: 10.11674/zwyf.2013.0611
    [6] 熊汉锋王运华 . 梁子湖湿地土壤养分的空间异质性. 植物营养与肥料学报, 2005, 13(5): 584-589. doi: 10.11674/zwyf.2005.0503
    [7] 姜城杨俐苹金继运张维理 . 土壤养分变异与合理取样数量. 植物营养与肥料学报, 2001, 9(3): 262-270. doi: 10.11674/zwyf.2001.0304
    [8] 许红卫高克异王珂周斌 . 稻田土壤养分空间变异与合理取样数研究. 植物营养与肥料学报, 2006, 14(1): 37-40. doi: 10.11674/zwyf.2006.0107
    [9] 王景云乔鹏炜杨军雷梅 . 基于空间结构特征的土壤Cu、Pb来源解析—以北京市东南郊污灌区为例. 植物营养与肥料学报, 2018, 24(1): 195-202. doi: 10.11674/zwyf.17195
    [10] 毛伟李文西高晖陈欣姜义杭天文龚鑫鑫陈明张月平 . 扬州市耕地土壤pH值30年演变及其驱动因子. 植物营养与肥料学报, 2017, 25(4): 883-893. doi: 10.11674/zwyf.16442
    [11] 王坷许红卫史舟 , . 土壤钾素空间变异性和空间插值方法的比较研究. 植物营养与肥料学报, 2000, 8(3): 318-322. doi: 10.11674/zwyf.2000.0312
    [12] 尹睿张华勇黄锦法林先贵王俊华曹志洪 . 保护地菜田与稻麦轮作田土壤微生物学特征的比较. 植物营养与肥料学报, 2004, 12(1): 57-62. doi: 10.11674/zwyf.2004.0111
    [13] 周涛胡华惠开基 . 宁南山区不同降水年型与地形的施肥模式研究. 植物营养与肥料学报, 2000, 8(4): 404-408. doi: 10.11674/zwyf.2000.0407
    [14] 赵亮张贺翠廉小平陆广涛朱利泉 . 喀斯特地形区水稻测土配方施肥指标体系研究. 植物营养与肥料学报, 2015, 23(4): 1056-1065. doi: 10.11674/zwyf.2015.0426
    [15] 南雅芳郭胜利李娜娜张彦军陈爱霞荣丽媛刘京 . 不同地形条件下青藏高原农田土壤有机碳的分布特征. 植物营养与肥料学报, 2013, 21(4): 946-954. doi: 10.11674/zwyf.2013.0422
    [16] 李小宇李勇于寒青张迎珍过治军 . 退耕还林坡地土壤CO2排放的空间变化: 地形的控制作用. 植物营养与肥料学报, 2015, 23(5): 1217-1224. doi: 10.11674/zwyf.2015.0515
    [17] 曹槐张晓林刘世熙和丽忠陈锦玉樊永言 . 烤烟矿质营养分布的因子分析. 植物营养与肥料学报, 2001, 9(3): 318-324. doi: 10.11674/zwyf.2001.0312
    [18] 李世清李生秀 . 影响土壤尿素水解速率的一些因子. 植物营养与肥料学报, 1999, 7(2): 156-162. doi: 10.11674/zwyf.1999.0209
    [19] 魏孝荣邵明安张兴昌高建伦 . 黄土沟壑区小流域不同地形条件下土壤锰的形态分布及其有效性 . 植物营养与肥料学报, 2008, 16(3): 439-444. doi: 10.11674/zwyf.2008.0305
    [20] 杨靖民杨靖一HoogenboomGerrit . DSSAT作物模型的统计和图形校验方法. 植物营养与肥料学报, 2012, 20(5): 1065-1074. doi: 10.11674/zwyf.2012.12006
  • 加载中
图(4)表(6)
计量
  • 文章访问数:  49
  • HTML全文浏览量:  30
  • PDF下载量:  11
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-08-21
  • 刊出日期:  2018-11-01

黑土区县域土壤养分空间分布特征及其影响因子

    作者简介:杜可 E-mail: 429474807@qq.com
    通讯作者: 张淑香, zhangshuxiang@caas.cn
  • 中国农业科学院农业资源与农业区划研究所,北京 100081

摘要:  目的 土壤养分是农业生产的必要条件,了解土壤养分的空间变化,对指导精准施肥有一定意义。 方法 本文以地统计学为基础,利用ArcGIS软件,对克山县土壤养分数据进行析取克里格以及泛克里格插值,得到土壤养分的空间分布图,并通过DEM数据提取坡度、坡向、曲率等相关地形因子,在SPSS软件中进行地形因子与土壤养分的相关分析。 结果 克山县土壤有机质、有效磷、速效钾和全氮含量平均值分别为43.89 g/kg、28.35 mg/kg、184.75 mg/kg和1.81 g/kg,有机质与全氮的变异系数均为0.09,属于弱变异强度,而有效磷与速效钾的变异系数分别为0.18与0.11,属于中等变异强度。有机质、速效钾及全氮的含量与海拔均呈显著的正相关关系(P<0.01),其中,有机质与全氮受海拔的影响程度最大。有机质以及全氮的含量还受到坡度的影响,与坡度呈显著的正相关性。不同土地利用方式间,旱地和菜地有机质和速效钾含量显著高于水田,而有效磷显著低于水田,全氮含量差异不显著。 结论 克山黑土区土壤有机质和全氮含量变异小,而有效磷和速效钾含量变异大。土壤养分的空间分布受到坡度与海拔的影响,尤其海拔对有机质与全氮的影响程度最大。土地利用类型、土壤类型以及成土母质也影响养分的分布。

English Abstract

  • 土壤是一个自然连续体,而空间变异性是其最重要的属性之一[1]。土壤作为作物生长养分的来源,其养分的变异必然会引起作物生长的变异[24]。土壤养分空间变异性研究对于科学合理地制定农田施肥方案,提高养分资源利用率,促进变量施肥技术的发展,实现精确施肥都具有重要意义[5]

    近年来,对土壤养分,例如有机质、有效磷、速效钾及全氮的空间分布的研究一直是一个热点问题。研究结果[67]表明,土壤养分的空间分布均不同程度地受到地形因子如坡度、坡向、曲率以及气候因子如降水、温度等的影响。如张文博等[6]发现有机质的含量随着海拔的升高而逐渐增加,阳坡的有机质含量较低,阴坡与半阴坡有机质含量较高;在国外也有类似的研究成果,如Smith等[7]的研究结果表明,土壤中全氮的含量与有机质相同,伴随海拔的升高而呈现增加的趋势。同时,土壤养分的分布情况也受到了其他因素的影响,例如土壤质地、土地利用类型等,张珍明等[8]的研究指出,不同的土地利用类型显著影响到全氮以及有效磷的含量,全氮的含量:沼泽草地 > 农用地 > 林地,而有效磷的含量:农用地 > 沼泽草地 > 林地。为了有效地了解土壤中养分的空间分布情况,近年来除了反距离权重法 (IDW) [910]以及径向基函数法 (RBF)[11],结合地统计学的克里格插值法[12]也被广泛应用。于洋等[13]在对研究区内平均降水量采用不同方法进行空间插值发现,在三种方法中,径向基函数法和反距离权重法产生的误差均高于普通克里格插值法。付传城等[14]对八种重金属采用不同的插值方法,对其插值精度、插值结果进行对比分析后发现,研究区内不同的重金属所对应的最优的插值方法不同,表明在实际的操作中,应根据具体的数据统计类型去选择最优的插值方法。Robert等[15]在对水生植物的研究中,加入了另一种插值方法不规则三角网 (TIN) 的插值方法,然而其插值结果并不理想,插值精度要低于普通克里格插值法 (Ordinary kriging)。在关于土壤养分的空间分布研究中,张铁婵等[16]在对该研究区内有机质、水解性氮、速效钾、有效磷以及全氮等因子为例,采用不同的加权方法发现,插值的结果受到采样密度的影响较大,反距离权重法与径向基函数法具体实施方法简单,参数易于选择,但其插值的效果要差于普通克里格插值法。秦玮等[17]的研究中却发现,对于有效磷和全氮的空间插值,局部多项式法和克里格插值法都具有较好的效果,对于有机质的空间估计而言,三种方法并未产生明显的区别,估值结果相似,速效钾则应用克里格插值法所产生的误差最小。然而在地形、气候等因子的影响下,普通的克里格插值法预测精度不高,克里格插值法中其他的插值方法如泛克里格法与析取克里格法的预测结果要更加精确[1820]

    从以往研究来看,大多在县域尺度上的研究中,采样点的数量一般为200左右,同时未考虑影响土壤养分空间变异的影响因子。克山县作为黑龙江省重点商品粮基地县,其土壤养分的含量、分布状况对其县内的产业结构的影响极其重要,对克山县土壤养分含量分布状况及其影响因子的研究,也可为其他同类地区的土壤养分管理提供参考价值。本研究在已有研究的基础上,应用经典统计学和地统计学相结合的方法,通过泛克里格插值与析取克里格插值法,探讨了东北黑土区克山县表层土壤养分的空间变异特征,并分析影响土壤养分的空间变异的因子,分析了克山县土壤养分的特征与分布规律及影响因子,这将为克山县针对当地土壤养分含量的监测与治理提供科学依据。

    • 研究区域为克山县,面积约3186.2平方公里,位于松嫩平原北部,属于黑龙江省齐齐哈尔市,介于东经125°10′57″~126°8′18″,北纬47°50′51″~48°33′47″。克山县海拔最高为381.7 m,最低为198.7 m,全境平均海拔236.9 m,属寒温带大陆季风气候,年平均气温2.4℃,有效积温2400℃,年降水量499 mm。克山县地处小兴安岭南麓与松嫩平原的过渡地带,北部、中部为丘陵漫岗区,地形起伏变化大,南部是冲积平原,地势东北高西南低,分为低山丘陵、丘陵边缘岗地、漫岗平原、冲积平原、河滩地等地貌类型,丘陵漫岗地占80%,平原区占14%,洼地占6%。土壤类型主要以黑土为主,主要种植作物为玉米和大豆。

    • 研究中布点以及采样的方法是按照全国耕地地力调查与质量评价技术规程进行,为了尽可能地避免施肥对于样品的影响,统一在作物收获之后进行取样。2010年 5月在研究区进行采样 (合计样点3706个),记录样点经纬度。土壤有机质测定采用重铬酸钾容量法、全氮测定采用半微量开氏法、速效钾测定采用乙酸铵提取火焰光度法、有效磷测定采用碳酸氢钠提取—钼锑抗比色法[21]。点位图和养分插值图编制采用ArcGIS地理信息系统软件 (图1)。

      图  1  采样点分布图

      Figure 1.  Sampling point distribution map

    • 运用IBM Statistics SPSS 20.0软件对研究区3706个采样点的有机质、有效磷、速效钾、全氮含量进行异常值处理,最常用的方法为拉依达准则法[2223],将标准化后的有机质、有效磷、速效钾、全氮含量绝对值大于2的点进行删除。处理后的有机质、有效磷、速效钾、全氮分别有3520、3545、3455、以及3524个样本点数据符合要求。

    • 克里格插值 (Kriging) 方法,是目前地统计学中最为常用的插值方法[24],插值的方法一般包括普通克里格 (Ordinary kriging)[25]、简单克里格 (Simple Kriging)[26]、泛克里格 (Universal kriging)[27]、协同克里格 (Co-kriging)[28]以及析取克里格 (Disjunctive kriging)[29]等。考虑到原始数据或者经数据转换[3031]后均可能不符合正态分布,本研究采取泛克里格插值以及析取克里格插值方法。

    • 对土壤中有机质、有效磷、速效钾以及全氮的有效样点进行经典统计学分析,得到了克山县土壤有效样点主要养分的描述性统计特征值 (表1)。表中结果表明:研究区内有机质、速效钾及全氮的含量与海拔均呈显著的正相关关系(P<0.01),其中,有机质与全氮受海拔的影响程度最大。有机质以及全氮的含量还受到坡度的影响,与坡度呈显著的正相关性。

      土壤养分
      Soilnutrients
      统计值Statistics 变异系数
      CV
      偏斜度
      Skewness
      峰度
      Kurtosis
      K-S检验
      K-S test
      样本数Sample No. 最小值Min. 最大值Max. 平均数Mean 中值Median
      有机质 SOM(g/kg) 3520 35.00 53.70 43.89 43.703 0.09 0.159 –0.431 4.56E-08
      有效磷 Avail. P(mg/kg) 3545 16.40 41.50 28.35 28.10 0.18 0.178 –0.469 2.40E-07
      速效钾 Avail. K(mg/kg) 3455 131.00 239.00 184.75 186.00 0.11 –0.157 0.316 1.40E-24
      全氮 Total N(g/kg) 3524 1.44 2.20 1.81 1.81 0.09 0.086 –0.578 2.20E-05

      表 1  土壤养分统计特征值

      Table 1.  Statistical characteristics of soil nutrients

      为进行泛克里格插值,对研究区内的有机质、有效磷、速效钾以及全氮含量进行趋势分析 (图2)。由图2可知,土壤养分中有机质与全氮的分布从南到北均呈现倒U形的分布趋势,从东到西则呈现缓慢增长的趋势;研究区内有效磷的含量从从西南到东北方向逐渐降低,而从东南到西北含量有增大的趋势;速效钾的含量从东到西,从南到北均呈现为倒U形的变化趋势。

      图  2  养分含量趋势分析

      Figure 2.  Trend analysis of nutrient contents

    • 因为本文中土壤养分的数据以及经过转换之后都不符合正态分布,所以采用泛克里格方法以及析取克里格法分别对区域内的区域变化量进行估值,分别得到有机质、有效磷、速效钾以及全氮的空间分布图 (图3)。

      图  3  泛克里格法 (左) 和析取克里格法 (右) 表征的土壤养分空间分布特征

      Figure 3.  Spatial distribution characteristics of soil nutrients described by universal (Left) and disjunctive kriging (Right)

      图3中泛克里格法与析取克里格法空间模拟插值结果显示,研究区内有机质与全氮的分布情况相似,高值集中分布在研究区中部地区的发展乡、西城镇以及北联镇等地以及最北面的北兴镇一带,有机质的含量处于44.50~54.00 g/kg之间,全氮含量处于1.81~2.20 g/kg之间,南部及北部部分地区有机质与全氮的含量较低,处于35.00~44.50 g/kg之间和1.43~1.81 g/kg之间。研究区内有效磷的含量表现为中部东边的发展乡、西联乡与西建乡与南部的双河乡含量较高,处于28.95~41.50 mg/kg之间,其他区域含量较低,处于16.40~28.95 mg/kg之间。速效钾的含量高值集中分布于研究区的中部与北部地区的西河镇、北联镇、北兴镇一带,主要处在185.00~239.00 mg/kg之间,南部与东部的双河乡与西村镇的含量较低,处在131.00~185.00 mg/kg之间。

    • 在对土壤养分的空间分布研究过程中,克里格插值的方法要明显优于其他几种空间插值方法,在本研究中主要采用了克里格插值法中的析取克里格与泛克里格法对研究区中的有机质、有效磷、速效钾和全氮进行估值。其中,泛克里格插值法是残差分析以及趋势分析两者相结合的结果,其插值结果的精度要远远优于趋势分析以及普通的克里格插值方法,从效果图来看,泛克里格插值结果图的空间连续性要更好,纹理更加清晰,对于土壤养分的预测更加细致、精确,插值的效果也要远远优于析取克里格插值的效果图。

      针对有机质、有效磷、速效钾以及全氮四种土壤养分在探究其空间分布特征时所采用的泛克里格插值与析取克里格插值两种方法,在已知插值结果的基础上,选取两种插值方法的标准平均值 (Standardized mean)、均方根预测误差 (Root-mean-square)、平均标准误差 (Average mean error)、标准均方根预测误差 (Root-mean-square standardized) 四个指标做预测结果的分析,其相关指标数据如表2

      项目
      Item
      插值方法
      Kriging
      标准平均值
      Standardized mean
      均方根预测误差
      Root-mean-square
      平均标准误差
      Average mean error
      标准均方根预测误差
      Root-mean-square standardized
      有机质SOM 泛克里格法Universal kriging –0.00262 2.40323 2.31138 1.03715
      析取克里格法Disjuctive kriging 0.00331 2.73138 2.91306 0.93739
      有效磷
      Avail. P
      泛克里格法Universal kriging –0.00536 3.44905 3.54562 0.95637
      析取克里格法Disjuctive kriging –0.00514 3.46199 3.94540 0.87078
      速效钾
      Avail. K
      泛克里格Universal kriging 0.00235 12.83784 13.86553 0.92381
      析取克里格Disjuctive kriging 0.00493 12.92830 15.63747 0.82311
      全氮
      Total N
      泛克里格Universal kriging 0.00186 0.09891 0.10055 0.97619
      析取克里格Disjuctive kriging 0.00211 0.11484 0.12190 0.94208

      表 2  泛克里格和析取克里格插值法相关数值

      Table 2.  Accuracy of parameters used in universal and disjunctive krigings

      根据统计学的原理,有机质、有效磷、速效钾、全氮的空间分布所采用的两种模型中,相关指标符合以下几种标准的则为最优模型:1) 标准平均值最接近于0;2) 均方根预测误差最小;3) 平均标准误差最接近于均方根预测误差;4) 标准均方根预测误差最接近于1。如表2所示,研究区有机质、有效磷、速效钾及全氮的空间插值方法均为泛克里格法最优[32]

      在理论上,在进行普通克里格插值的时候要求数据具有正态分布,然而在实际的工作中,调查的数据很难符合这一要求,所以在计算前,通常会对数据进行数据转换,如最常见的对数转换以及Box-Cox转换,正态分布的转换方法有很多种,其中有一些方法甚至是不可逆的,这严重影响了对插值结果的分析。即使面对一些可以转换的数据,在对其进行转换之后,进行克里格插值,然后对克里格插值的结果进行逆转换,克里格的估值必然发生错误,所以,在实际分析的时候,可以选择避开数据分布的问题,采用析取克里格方法进行空间估值。

    • 土壤养分的空间分布特征普遍受到地形因子,如坡度、坡向、海拔、曲率等因素的影响。通过下载研究区内克山县30 m × 30 m精度的DEM数据(地理空间数据云http://www.gscloud.cn/),通过ArcGIS软件中的空间分析功能得到克山县坡度[33]、坡向[34]、曲率[35]等空间分布图 (图4)。

      图  4  地形因子空间分布特征

      Figure 4.  Terrain factor characteristics of soil nutrients

      利用ArcGIS软件提取各个样本点位的地形因子数据,通过SPSS软件,得到研究区内主要土壤养分与海拔、坡度、坡向、曲率等多种影响因子的相关性分析结果如表3

      项目
      Item
      有机质SOM 有效磷Avail. P 速效钾Avail. K 全氮Total N
      相关性Correlation 显著性Significance 相关性Correlation 显著性Significance 相关性Correlation 显著性Significance 相关性Correlation 显著性Significance
      海拔Latitude 0.391** 0.000 0.029 0.088 0.141** 0.000 0.388** 0.000
      坡度Gradient 0.054** 0.001 –0.020 0.233 0.028 0.104 0.051** 0.002
      坡向Aspect 0.029 0.081 0.001 0.935 0.024 0.157 0.032 0.057
      曲率Curvature 0.008 0.633 –0.017 0.324 –0.015 0.376 0.003 0.881
      注(Note):**表示在 0.01 水平上显著相关 Indicates a significant correlation at the 0.01 level.

      表 3  土壤养分与环境因子相关性分析

      Table 3.  Correlation analysis between soil nutrients and environmental factors

      表3结果表明,研究区内有机质、速效钾及全氮的含量与海拔在0.01水平上均呈显著的正相关关系,这表明,随着海拔的升高,研究区内有机质、速效钾及全氮的含量均呈增加的趋势,其中,有机质与全氮受海拔的影响程度最大。有机质以及全氮的含量还受到坡度的影响,与坡度呈显著的正相关性,表明随着坡度的增大,研究区内有机质与全氮的含量随之增加。

    • 除了坡度、坡向、海拔、曲率等因素影响土壤中主要养分的含量以外,不同的土地利用方式、土壤类型、成土母质等因素均可在较大范围内影响土壤中养分的含量。从土地利用方式来看 (表4),菜地的有机质含量最高为43.95 g/kg,与其他两种类型并无显著差异,旱地的有效磷含量最高为28.39 mg/kg,含量明显高于菜地与灌溉水田,菜地中速效钾的含量最高为197.50 mg/kg,明显高于其他两种类型,菜地的全氮含量最低为1.78 g/kg,与其他两种用地差异显著。从土壤类型来看 (表5),棕壤的有机质含量最高为44.61 g/kg,但是与除黑钙土外其他土类差距并不明显,黑钙土中有机质含量最低,明显低于其他几种土壤;黑钙土的有效磷含量最高为30.89 mg/kg,草甸土的有效磷含量最低为27.96 mg/kg;黑钙土的速效钾与全氮的含量要明显低于其他几种土壤,其含量分别为171.33 mg/kg和1.66 g/kg。从成土母质来看 (表6)。坡积物的有机质与全氮含量最高,分别为44.94 g/kg和1.86 g/kg,淤积物的有效磷含量最高为29.21 mg/kg,几种母质的有机质、有效磷与全氮含量并无显著差异。火山岩的速效钾含量最高,为205.00 mg/kg,明显高于其他几种成土母质的含量。

      土地利用方式
      Land use
      有机质SOM (g/kg) 有效磷Avail. P (mg/kg) 速效钾Avail. K (mg/kg) 全氮Total N (g/kg)
      n 含量Content n 含量Content n 含量Content n 含量Content
      旱地Dry land 3417 43.90 ± 3.82 3439 28.39 ± 5.24 3349 184.60 ± 21.22 3421 1.81 ± 0.16
      水田Paddy 97 43.61 ± 2.91 100 26.88 ± 4.71 100 189.14 ± 12.58 97 1.83 ± 0.12
      菜地Vegetable 6 43.95 ± 0.54 6 26.67 ± 0.67 6 197.50 ± 7.29 6 1.78 ± 0.02

      表 4  不同土地利用类型的土壤养分样本数 (n) 和含量

      Table 4.  Sample number (n) and soil nutrient contents of different land uses

      土类
      Soil type
      有机质 SOM (g/kg) 有效磷 Avail. P (mg/kg) 速效钾 Avail. K (mg/kg) 全氮 Total N (g/kg)
      n 含量Content n 含量Content n 含量Content n 含量Content
      黑土Black soil 2173 44.25 ± 3.89 2178 28.49 ± 5.30 2142 184.79 ± 21.37 2181 1.82 ± 0.16
      棕壤Brown soil 298 44.61 ± 3.77 310 28.04 ± 5.02 285 187.48 ± 21.19 303 1.84 ± 0.15
      草甸土Meadow soil 967 43.14 ± 3.43 997 27.96 ± 5.12 960 184.55 ± 19.39 961 1.78 ± 0.15
      黑钙土Chernozems 61 40.80 ± 3.56 59 30.89 ± 4.58 54 171.33 ± 27.36 58 1.66 ± 0.15
      沼泽土Wetland 21 43.35 ± 2.41 21 29.21 ± 4.58 14 189.64 ± 26.16 21 1.78 ± 0.11

      表 5  不同土壤类型的土壤样本数(n)和养分含量

      Table 5.  Sample number (n) and soil nutrient contents of different soil types

      成土母质
      Parent material
      有机质 SOM (g/kg) 有效磷 Avail. P (mg/kg) 速效钾 Avail. K (mg/kg) 全氮 Total N (g/kg)
      n 含量Content n 含量Content n 含量Content n 含量Content
      黄土母质Loess 2244 44.14 ± 3.92 2243 28.49 ± 5.31 2192 184.23 ± 21.90 2246 1.82 ± 0.16
      冲积物Alluvium 1163 43.36 ± 3.56 1182 28.08 ± 5.12 1158 185.19 ± 19.22 1165 1.79 ± 0.15
      坡积物Slope deposit 75 44.94 ± 3.46 82 28.24 ± 4.97 75 190.63 ± 19.98 76 1.86 ± 0.13
      玄武岩Basalt 14 43.72 ± 2.09 14 27.11 ± 2.11 13 190.46 ± 18.04 13 1.81 ± 0.09
      火山岩Lava 3 43.70 ± 0.40 3 28.87 ± 0.75 3 205.00 ± 3.00 3 1.82 ± 0.02
      淤积物Sediment 21 43.35 ± 2.41 21 29.21 ± 4.58 14 189.64 ± 26.16 21 1.78 ± 0.11

      表 6  不同成土母质的土壤样本数(n)和养分含量

      Table 6.  Sample number (n) and soil nutrient contents of different soil parent materials

    • 土壤养分的含量除受到地形因子、土壤类型、土地利用类型等因素的影响外,气候因素如降水、温度等对土壤养分含量的影响也起到关键作用。曹祥会等[36]的研究表明,气候因素主要通过影响土壤中微生物的活性、养分的淋溶迁移过程以及植被本身的生产力,从而达到影响土壤中养分的分解与输入,进一步影响其含量的作用。赵艳艳等[37]发现,长期的增温会使灌丛和草甸土壤中速效磷和全磷的含量增加,而全钾的含量降低。而崔贝等[38]的研究指出,对于有机质而言,低温的环境,通过抑制微生物的分解作用,有利于土壤中有机质的积累,较大的降水量则会增大土壤中养分的淋溶,进而影响土壤表层养分的含量,这与曹祥会等的研究结果类似。孙大生[39]的试验表明,土壤可溶性有机碳、单糖、矿化氮、活性无机磷及有机磷含量随干旱强度的增加而升高,土壤微生物生物量碳氮磷及溶解态有机氮含量随干旱强度的增加而降低。

      人类的活动对自然因素的影响以及对土壤养分含量的直接影响也越来越大。许中旗等[40]研究了人类干扰活动对土壤养分含量的影响,结果表明,土壤养分含量受人类活动影响的大小随着土层深度的增加而逐渐降低,人类的活动如放牧、开垦等行为使土壤表层养分含量明显下降。同时,化肥的施用量更是直接影响到土壤中养分的含量,从而使得土壤养分的空间变异性产生较大的变化[41],人类的活动对土壤养分分布的干扰程度不容忽视,在今后的分析中,应该将人类活动这一因素考虑进去,建立一个更加科学、完善的指标体系。

    • 1) 研究区内有机质、有效磷、速效钾和全氮含量的平均值分别为43.89 g/kg、28.35 mg/kg、184.75 mg/kg和1.81 g/kg,有机质与全氮的变异系数均为0.09,属于弱变异强度,而有效磷与速效钾的变异系数分别为0.18与0.11,属于中等变异强度。

      2) 研究区内有机质、有效磷、速效钾和全氮均不符合正态分布,经过对数转换和Box-Cox转换之后仍不符合正态分布,经过数据的趋势分析发现土壤养分的数据存在明显的主导趋势,故采用泛克里格方法以及避开数据分布问题的析取克里格方法对其进行空间插值。且通过对比空间插值效果图以及相关指标,泛克里格插值要优于析取克里格插值。

      3) 土壤有机质、速效钾及全氮的含量与海拔均呈极显著的正相关关系,其中,有机质与全氮受海拔的影响程度最大,也与坡度呈显著的正相关性。

参考文献 (41)

目录

    /

    返回文章
    返回