土壤是一个自然连续体,而空间变异性是其最重要的属性之一[1]。土壤作为作物生长养分的来源,其养分的变异必然会引起作物生长的变异[2–4]。土壤养分空间变异性研究对于科学合理地制定农田施肥方案,提高养分资源利用率,促进变量施肥技术的发展,实现精确施肥都具有重要意义[5]。
近年来,对土壤养分,例如有机质、有效磷、速效钾及全氮的空间分布的研究一直是一个热点问题。研究结果[6–7]表明,土壤养分的空间分布均不同程度地受到地形因子如坡度、坡向、曲率以及气候因子如降水、温度等的影响。如张文博等[6]发现有机质的含量随着海拔的升高而逐渐增加,阳坡的有机质含量较低,阴坡与半阴坡有机质含量较高;在国外也有类似的研究成果,如Smith等[7]的研究结果表明,土壤中全氮的含量与有机质相同,伴随海拔的升高而呈现增加的趋势。同时,土壤养分的分布情况也受到了其他因素的影响,例如土壤质地、土地利用类型等,张珍明等[8]的研究指出,不同的土地利用类型显著影响到全氮以及有效磷的含量,全氮的含量:沼泽草地 > 农用地 > 林地,而有效磷的含量:农用地 > 沼泽草地 > 林地。为了有效地了解土壤中养分的空间分布情况,近年来除了反距离权重法 (IDW) [9–10]以及径向基函数法 (RBF)[11],结合地统计学的克里格插值法[12]也被广泛应用。于洋等[13]在对研究区内平均降水量采用不同方法进行空间插值发现,在三种方法中,径向基函数法和反距离权重法产生的误差均高于普通克里格插值法。付传城等[14]对八种重金属采用不同的插值方法,对其插值精度、插值结果进行对比分析后发现,研究区内不同的重金属所对应的最优的插值方法不同,表明在实际的操作中,应根据具体的数据统计类型去选择最优的插值方法。Robert等[15]在对水生植物的研究中,加入了另一种插值方法不规则三角网 (TIN) 的插值方法,然而其插值结果并不理想,插值精度要低于普通克里格插值法 (Ordinary kriging)。在关于土壤养分的空间分布研究中,张铁婵等[16]在对该研究区内有机质、水解性氮、速效钾、有效磷以及全氮等因子为例,采用不同的加权方法发现,插值的结果受到采样密度的影响较大,反距离权重法与径向基函数法具体实施方法简单,参数易于选择,但其插值的效果要差于普通克里格插值法。秦玮等[17]的研究中却发现,对于有效磷和全氮的空间插值,局部多项式法和克里格插值法都具有较好的效果,对于有机质的空间估计而言,三种方法并未产生明显的区别,估值结果相似,速效钾则应用克里格插值法所产生的误差最小。然而在地形、气候等因子的影响下,普通的克里格插值法预测精度不高,克里格插值法中其他的插值方法如泛克里格法与析取克里格法的预测结果要更加精确[18–20]。
从以往研究来看,大多在县域尺度上的研究中,采样点的数量一般为200左右,同时未考虑影响土壤养分空间变异的影响因子。克山县作为黑龙江省重点商品粮基地县,其土壤养分的含量、分布状况对其县内的产业结构的影响极其重要,对克山县土壤养分含量分布状况及其影响因子的研究,也可为其他同类地区的土壤养分管理提供参考价值。本研究在已有研究的基础上,应用经典统计学和地统计学相结合的方法,通过泛克里格插值与析取克里格插值法,探讨了东北黑土区克山县表层土壤养分的空间变异特征,并分析影响土壤养分的空间变异的因子,分析了克山县土壤养分的特征与分布规律及影响因子,这将为克山县针对当地土壤养分含量的监测与治理提供科学依据。
1 材料与方法 1.1 研究区概况研究区域为克山县,面积约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%。土壤类型主要以黑土为主,主要种植作物为玉米和大豆。
1.2 数据采集与测试分析研究中布点以及采样的方法是按照全国耕地地力调查与质量评价技术规程进行,为了尽可能地避免施肥对于样品的影响,统一在作物收获之后进行取样。2010年 5月在研究区进行采样 (合计样点3706个),记录样点经纬度。土壤有机质测定采用重铬酸钾容量法、全氮测定采用半微量开氏法、速效钾测定采用乙酸铵提取火焰光度法、有效磷测定采用碳酸氢钠提取—钼锑抗比色法[21]。点位图和养分插值图编制采用ArcGIS地理信息系统软件 (图1)。
运用IBM Statistics SPSS 20.0软件对研究区3706个采样点的有机质、有效磷、速效钾、全氮含量进行异常值处理,最常用的方法为拉依达准则法[22–23],将标准化后的有机质、有效磷、速效钾、全氮含量绝对值大于2的点进行删除。处理后的有机质、有效磷、速效钾、全氮分别有3520、3545、3455、以及3524个样本点数据符合要求。
1.4 克里格插值克里格插值 (Kriging) 方法,是目前地统计学中最为常用的插值方法[24],插值的方法一般包括普通克里格 (Ordinary kriging)[25]、简单克里格 (Simple Kriging)[26]、泛克里格 (Universal kriging)[27]、协同克里格 (Co-kriging)[28]以及析取克里格 (Disjunctive kriging)[29]等。考虑到原始数据或者经数据转换[30–31]后均可能不符合正态分布,本研究采取泛克里格插值以及析取克里格插值方法。
2 结果与分析 2.1 土壤中主要养分描述性统计特征及趋势分析对土壤中有机质、有效磷、速效钾以及全氮的有效样点进行经典统计学分析,得到了克山县土壤有效样点主要养分的描述性统计特征值 (表1)。表中结果表明:研究区内有机质、速效钾及全氮的含量与海拔均呈显著的正相关关系(P<0.01),其中,有机质与全氮受海拔的影响程度最大。有机质以及全氮的含量还受到坡度的影响,与坡度呈显著的正相关性。
为进行泛克里格插值,对研究区内的有机质、有效磷、速效钾以及全氮含量进行趋势分析 (图2)。由图2可知,土壤养分中有机质与全氮的分布从南到北均呈现倒U形的分布趋势,从东到西则呈现缓慢增长的趋势;研究区内有效磷的含量从从西南到东北方向逐渐降低,而从东南到西北含量有增大的趋势;速效钾的含量从东到西,从南到北均呈现为倒U形的变化趋势。
因为本文中土壤养分的数据以及经过转换之后都不符合正态分布,所以采用泛克里格方法以及析取克里格法分别对区域内的区域变化量进行估值,分别得到有机质、有效磷、速效钾以及全氮的空间分布图 (图3)。
图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之间。
2.3 不同的插值方法对预测结果的影响在对土壤养分的空间分布研究过程中,克里格插值的方法要明显优于其他几种空间插值方法,在本研究中主要采用了克里格插值法中的析取克里格与泛克里格法对研究区中的有机质、有效磷、速效钾和全氮进行估值。其中,泛克里格插值法是残差分析以及趋势分析两者相结合的结果,其插值结果的精度要远远优于趋势分析以及普通的克里格插值方法,从效果图来看,泛克里格插值结果图的空间连续性要更好,纹理更加清晰,对于土壤养分的预测更加细致、精确,插值的效果也要远远优于析取克里格插值的效果图。
针对有机质、有效磷、速效钾以及全氮四种土壤养分在探究其空间分布特征时所采用的泛克里格插值与析取克里格插值两种方法,在已知插值结果的基础上,选取两种插值方法的标准平均值 (Standardized mean)、均方根预测误差 (Root-mean-square)、平均标准误差 (Average mean error)、标准均方根预测误差 (Root-mean-square standardized) 四个指标做预测结果的分析,其相关指标数据如表2。
根据统计学的原理,有机质、有效磷、速效钾、全氮的空间分布所采用的两种模型中,相关指标符合以下几种标准的则为最优模型:1) 标准平均值最接近于0;2) 均方根预测误差最小;3) 平均标准误差最接近于均方根预测误差;4) 标准均方根预测误差最接近于1。如表2所示,研究区有机质、有效磷、速效钾及全氮的空间插值方法均为泛克里格法最优[32]。
在理论上,在进行普通克里格插值的时候要求数据具有正态分布,然而在实际的工作中,调查的数据很难符合这一要求,所以在计算前,通常会对数据进行数据转换,如最常见的对数转换以及Box-Cox转换,正态分布的转换方法有很多种,其中有一些方法甚至是不可逆的,这严重影响了对插值结果的分析。即使面对一些可以转换的数据,在对其进行转换之后,进行克里格插值,然后对克里格插值的结果进行逆转换,克里格的估值必然发生错误,所以,在实际分析的时候,可以选择避开数据分布的问题,采用析取克里格方法进行空间估值。
3 讨论 3.1 地形因子与土壤养分空间分布的相关性分析土壤养分的空间分布特征普遍受到地形因子,如坡度、坡向、海拔、曲率等因素的影响。通过下载研究区内克山县30 m × 30 m精度的DEM数据(地理空间数据云http://www.gscloud.cn/),通过ArcGIS软件中的空间分析功能得到克山县坡度[33]、坡向[34]、曲率[35]等空间分布图 (图4)。
利用ArcGIS软件提取各个样本点位的地形因子数据,通过SPSS软件,得到研究区内主要土壤养分与海拔、坡度、坡向、曲率等多种影响因子的相关性分析结果如表3。
表3结果表明,研究区内有机质、速效钾及全氮的含量与海拔在0.01水平上均呈显著的正相关关系,这表明,随着海拔的升高,研究区内有机质、速效钾及全氮的含量均呈增加的趋势,其中,有机质与全氮受海拔的影响程度最大。有机质以及全氮的含量还受到坡度的影响,与坡度呈显著的正相关性,表明随着坡度的增大,研究区内有机质与全氮的含量随之增加。
3.2 土壤质地、利用类型与土壤养分含量的关系除了坡度、坡向、海拔、曲率等因素影响土壤中主要养分的含量以外,不同的土地利用方式、土壤类型、成土母质等因素均可在较大范围内影响土壤中养分的含量。从土地利用方式来看 (表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,明显高于其他几种成土母质的含量。
土壤养分的含量除受到地形因子、土壤类型、土地利用类型等因素的影响外,气候因素如降水、温度等对土壤养分含量的影响也起到关键作用。曹祥会等[36]的研究表明,气候因素主要通过影响土壤中微生物的活性、养分的淋溶迁移过程以及植被本身的生产力,从而达到影响土壤中养分的分解与输入,进一步影响其含量的作用。赵艳艳等[37]发现,长期的增温会使灌丛和草甸土壤中速效磷和全磷的含量增加,而全钾的含量降低。而崔贝等[38]的研究指出,对于有机质而言,低温的环境,通过抑制微生物的分解作用,有利于土壤中有机质的积累,较大的降水量则会增大土壤中养分的淋溶,进而影响土壤表层养分的含量,这与曹祥会等的研究结果类似。孙大生[39]的试验表明,土壤可溶性有机碳、单糖、矿化氮、活性无机磷及有机磷含量随干旱强度的增加而升高,土壤微生物生物量碳氮磷及溶解态有机氮含量随干旱强度的增加而降低。
人类的活动对自然因素的影响以及对土壤养分含量的直接影响也越来越大。许中旗等[40]研究了人类干扰活动对土壤养分含量的影响,结果表明,土壤养分含量受人类活动影响的大小随着土层深度的增加而逐渐降低,人类的活动如放牧、开垦等行为使土壤表层养分含量明显下降。同时,化肥的施用量更是直接影响到土壤中养分的含量,从而使得土壤养分的空间变异性产生较大的变化[41],人类的活动对土壤养分分布的干扰程度不容忽视,在今后的分析中,应该将人类活动这一因素考虑进去,建立一个更加科学、完善的指标体系。
4 结论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) 土壤有机质、速效钾及全氮的含量与海拔均呈极显著的正相关关系,其中,有机质与全氮受海拔的影响程度最大,也与坡度呈显著的正相关性。
[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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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 |