利用高分立体影像构建东北黑土山地丘陵区切沟体积估算模型

分析侵蚀沟形态特征并构建体积估算模型,对大空间尺度探究切沟侵蚀具有重要意义。该研究选取位于东北黑土山地丘陵区穆棱市的典型样区,基于WorldView-2高分立体像对影像(0.5m分辨率),利用ENVI5.3提取同分辨率DEM,选择45条切沟,在ArcGIS10.3中提取并计算切沟顶宽(TW)、底宽(BW)、沟深(D)、断面面积(CSA)、沟长(L)、面积(A)及体积(V)等参数,分析形态特征及相互关系,构建切沟体积估算模型。结果表明:1)切沟TW、BW、D、CSA、L、A及V均值分别为5.63 m、2.82 m、1.05 m、4.94 m~2、81.90 m、470.54 m~2、428.53 m~3。断面BW/TW均值为0.48,切沟以V～U型为主。TW/D全部大于1,均值为5.95,表明沟道横向侵蚀较下切侵蚀速率更快;2)切沟体积与沟长(V-L)、面积(V-A)间均有极显著幂函数关系,相比V-L关系模型,V-A关系模型具有更小的平均相对误差和更大的纳什系数,可更准确有效预测切沟体积,建议将其作为东北山地丘陵区切沟体积估算模型;3)与黄土高原及干热河谷地区相比,东北黑土区切沟体积与面积幂函数关系的指数更小,而区内山地丘陵区较漫川漫岗区更大,反映出山地丘陵区更加强烈的切沟发育状况;切沟形态特征影响V-A模型构建的精度,建议将狭长度(L/TW)作为切沟分类指标构建切沟体估算模型。研究结果可为东北黑土山地丘陵区大空间尺度切沟侵蚀的定量模拟提供方法和依据。     

松嫩典型黑土区耕地切沟密度分布特征及影响因子

东北黑土区是中国重要的商品粮生产基地，该地区耕地切沟侵蚀严重威胁农业生产及生态环境。为掌握松嫩典型黑土区耕地切沟密度分布特征及影响因子，基于全国第一次水利普查的998块抽样调查单元和1∶10000地形图，以及2018年Google高分遥感影像目视解译提取调查单元内所有切沟，采用地理探测器分析切沟密度空间分布的影响因子，分析松嫩典型黑土区耕地切沟密度的空间分布规律。结果表明：1)35.64%的耕地调查单元有切沟发育，切沟线密度变化在0.18～17.48 km/km²，面密度变化在425.95～91 921.34 m2/km²。2)地形和降水显著(P<0.05)影响松嫩典型黑土区耕地切沟密度，坡度是影响耕地切沟密度空间分布最重要的环境因子。3)依据切沟线密度划分，松嫩典型黑土区东部的黑土带耕地切沟侵蚀严重，轻度侵蚀以上面积占比35.87%，强烈侵蚀以上(切沟线密度≥2 km/km²)面积占比33.33%，剧烈侵蚀面积占比27.94%。黑钙土带有6.65%的面积遭受剧烈切沟侵蚀，星落分布于黑钙土带的边缘。系统抽样与高分辨...     

黄土丘陵沟壑区小流域侵蚀沟数量及形态特征

为分析黄土高原不同类型侵蚀沟形态特征及其区域差异,以黄土高原第一副区和第五副区小流域为例,基于无人机影像(0.15 m)及DEM(1 m),分析沟谷网络特征和各类型侵蚀沟的数量特征,并对比冲沟和切沟的二维、三维形态参数。结果表明：丘一区小流域主沟支沟比为1∶35.6,远小于丘五区小流域主沟支沟比(1∶22.0),且各级沟谷密度多高于丘五区,表明丘一区支沟数量多,发育程度较高,地表更破碎,沟谷系统形态更复杂。丘一区和丘五区侵蚀沟类型都以切沟为主,密度分别为9.72,14.22 km/km²,冲沟密度分别为4.27,4.23 km/km²,冲沟数量虽少,但沟谷面积较大。切沟在沟谷网络分级中主要集中在1级,冲沟主要分布于1～3级,干沟多为3～5级。丘一区切沟除纵比降较大外,其他形态参数均明显小于丘五区,说明丘五区切沟规模更大,但两区冲沟形态差异较小。切沟、冲沟的体积和面积之间都具有显著的幂函数关系,但丘一区冲沟和切沟的判定系数更高,表明其同类侵蚀沟在形态上更具有相似性。研究结果为黄土高原侵蚀沟准确统计和小流域水土保持综合治理提供参考。     

晋西黄土区切沟断面特征及体积估算模型

断面特征不仅能反映切沟发育环境因子,亦有助于揭示切沟侵蚀规律;体积是表征切沟侵蚀量的关键参数,建立体积估算模型对大空间尺度切沟侵蚀定量研究有重要意义。该文选取晋西黄土区蔡家川流域为研究区,利用三维激光扫描全站仪获取31条切沟地形数据,构建切沟高精度地面高程模型(DEM),分析切沟断面特征,建立体积估算模型。研究结果表明:(1)切沟顶宽、底宽、沟深和断面面积的均值分别为9.0 m、3.8 m、6.2 m和41.8 m2,切沟断面参数值从沟头、沟中到沟口均逐渐增大,仅沟深差异不显著;断面顶宽与沟深之比的均值为1.55,说明切沟横向侵蚀速率大于沟底下切速率;与沟中、沟口相比,沟头宽深比最小。(2)回归分析显示,晋西黄土区切沟体积(V)与沟长(L)及面积(A)之间均具有显著的幂函数关系,方程的决定系数分别为0.68和0.85;平均相对误差和纳什系数显示,与V–L模型相比,V–A在预测切沟体积上更具优势。因此,建议将回归方程V=1.7097A1.1356作为晋西黄土区切沟体积估算的模型。该研究结果可为大空间尺度切沟侵蚀定量模拟、空间分布制图等提供方法和依据。     

晋西黄土区蔡家川小流域切沟的空间分布及形态特征

为了掌握切沟的空间分布及发育特点,实现小流域综合治理中切沟治理措施的精准对位配置,该研究以晋西黄土区吉县蔡家川小流域为对象,利用Quickbird影像和数字高程模型（DEM）对流域内370条切沟的形态特征和空间分布进行了提取和解译,并对其中40条切沟进行实地调查验证,采用核密度指标探讨切沟在小流域内的空间分布特征,分析坡度、坡向、高程等地形因子对切沟空间分异特征的影响。结果表明：1）蔡家川小流域内农地、果园的主要分布区（东北部和东南部）是切沟的高聚集区,而在人工林为主的中上游区属于中密集区;2）切沟数量、密度、频度和规模均随坡度增加呈现先增大后减小的趋势,坡度阈值为25°。 5°～25°的坡面切沟广泛发育,数量占研究区的71.5%,密度、频度也远高于其他坡度的坡面;3）切沟密度和频度的坡向分异与流域走向有关。切沟密度北坡最大（3.03 km/km2）,东坡最小（0.86 km/km2）,切沟频度西北坡最大（12.7条/km2）,东坡最小（6.1条/ km2）。南北向的切沟规模明显大于东西向;4）切沟数量、密度、频度随海拔增加呈现先增大后减小的趋势,在海拔为1 000～1 200 m的坡面上切沟数量占77%,海拔为1 100～1 200 m的坡面上切沟密度最大（2.17 km/km2）,海拔为1 000～1 100 m的坡面上切沟频度最高（15.5条/ km2）。随海拔升高切沟逐渐向大型化发展,沟道越发平缓,下切侵蚀逐渐减弱。研究结果可揭示流域切沟发育的重点地形以及在地形上的共性特点,为黄土区小流域切沟治理提供依据,可为科学认知切沟空间分布提供参考。     

基于Google Earth影像的横断山区沟蚀及侵蚀沟类型调查研究

[目的]探讨横断山区侵蚀沟的类型、分布特征及影响因素,为研究沟蚀对横断山区土壤侵蚀及土地退化的影响提供依据。[方法]基于Google Earth高清影像结合GIS的方法,基于可见性的基本原则,对横断山区开展侵蚀沟的抽样调查。共布设调查单元2 242个,每个单元尺寸1km×1km。[结果]共有571个调查单元发现有侵蚀沟分布,占调查单元总数的25.5%,平均沟壑密度为2.20km/km2。沟蚀在海拔低于1 500m的区域发育显著,沟壑密度随坡度增加而增长,与年降雨相关性较弱。草地是沟蚀发育的主要土地利用类型,燥红土的沟蚀分布率及平均沟壑密度显著高于其他土壤类型。通过总结国内外现有的侵蚀沟分类体系及其异同,结合调查结果,判定该次调查的横断山区侵蚀沟类型主要为切沟和冲沟(包括有泥石流事件的冲沟)。[结论]横断山区切沟和冲沟分布广泛,发育强烈,其分布受海拔、坡度、土地利用及土壤类型等因素的影响较为显著。     

黄土高原沟壑丘陵区沟道侵蚀与洞穴侵蚀特征

[目的]通过对黄土高原沟壑丘陵区内的沟道侵蚀和洞穴侵蚀进行特征分析,为该地区开展水土保持工作提供数据支持和决策依据。[方法]以王家沟流域为研究对象,通过资料收集、卫星影像识别和现场调查等方法对流域内沟道侵蚀和洞穴侵蚀状况进行分析。[结果](1)流域内共分布718条沟道和948个洞穴。(2)根据空间位置和地形特点将所有沟道分为冲沟、切沟、干沟、河沟。这4类沟道由于受土壤岩性组成、坡度及人为活动等因素影响,空间分布存在巨大差异。(3)43%的洞穴分布在冲沟中,51%分布在切沟中,6%分布在干沟中。河沟中没有洞穴存在。冲沟和切沟是洞穴形成和发育的主要区域。(4)洞穴沟道的沟长、面积、切深、产流及产沙量都明显大于非洞穴沟道。[结论]在黄土高原丘陵沟壑区,洞穴侵蚀是土壤侵蚀的重要组成部分,并对沟道的发育和扩展起着至关重要的作用。     

近50年来王茂沟流域侵蚀沟变化及其影响因素

基于1968年、2004年和2018年3期高分遥感影像,以王茂沟流域为典型研究区,提取了3个时间点的沟沿线、土地利用、水土保持措施、植被覆盖度和LS因子,对近50年侵蚀沟时空变化特征及其影响因素进行了分析。结果表明:(1)用遥感方法采集的信息,可基本满足侵蚀沟中长期变化分析;(2)50年来王茂沟流域侵蚀沟一直处于变化中,时间上变化速度逐渐减缓,前期和后期沟头年均前进速率分别为0.30,0.27m/a,沟壁年均扩张速率分别为0.009,0.004 m/a;空间上2个时期年均侵蚀沟变化密度大于200 m/(km^2·a)的值分别占比26.89%,7.07%,且多分布在流域中下游;(3)1968—2004年侵蚀沟变化主要受到土地利用变化和水土保持措施增加的影响,2005—2018年,侵蚀沟变化主要受到植被覆盖度进一步提高和土地利用多样性的影响。研究有助于量化分析多年治理过程中侵蚀沟的变化及其影响因素,能对退耕还林还草工程实施前后侵蚀沟治理的生态效益评价提供重要参考。     

晋西黄土区切沟对其上方集水区内土壤水分的影响

为掌握黄土高原切沟对其上方集水区内土壤水分的影响,选取晋西黄土区典型切沟上方的集水区为研究对象,采用土钻法定期测量0—10 m土层的土壤水分,探讨集水区内不同地貌部位的土壤水分状况。结果表明:(1)集水区内土壤水分状况受距离切沟远近和地貌部位的共同影响。(2)集水区内距离切沟1 m处的蓄水量显著低于距离切沟3,6 m处的蓄水量(P＜0.05),集水区内洼地的蓄水量大于坡面。(3)集水区内6个采样点浅层含水量随时间变化剧烈,深层含水量随时间变化小。(4)切沟对其上方集水区内土壤水分的水平影响距离为1~3 m,距切沟1 m处,影响深度最深可达7 m,距切沟3 m处,影响深度为2 m。在未来开展切沟治理中,建议在距离切沟1 m范围内密植低耗水的耐旱灌木或草本,利用密植的植物增加地面糙率,降低水流流速,并利用植物根系网络土体,增加抗冲性,从而起到防治切沟溯源侵蚀。     

陕北黄土区植被恢复对切沟发育速率的影响

切沟是黄土高原常见的地貌形态,是流域产沙的主要来源,而植被常被看作抑制切沟发育的关键因子,为探究退耕还林(草)工程以来切沟发育速率的特征,本研究使用吴起研究区间隔6年的2期同时相QuickBird影像和分辨率为5m的DEM,应用GIS技术分析研究区内的切沟发育速率及其影响因子.结果表明:2007-2013年研究区内小流域的年均切沟面积变化比例均＜0.5％,其中无淤地坝小流域的年均切沟面积变化比例为0.13％,淤地坝小流域的年均切沟面积变化比例为0.12％;相关分析表明:在该研究区内沟间地的植被覆盖度是影响切沟发育的最主要影响因子;通过对30个无淤地坝小流域的年均切沟面积变化比例和植被覆盖度做非线性回归(R2 =0.713)可得,当植被覆盖度＞0.6时,其抑制切沟发育的效果最明显;方差分析表明:其他影响因子值(植被覆盖度和地貌因子)相近的条件下,有无淤地坝的6组小流域切沟发育速率无显著差异.总之,沟间地植被覆盖度是影响切沟发育的主要影响因子,而淤地坝对发育在沟缘线上的切沟无显著影响,促进植被恢复有利于抑制切沟发育.     

特大暴雨下不同土地利用类型坡面切沟发育特征

切沟侵蚀是黄土高原丘陵沟壑区水土流失的重要形式之一，然而极端暴雨条件下不同土地利用类型坡面切沟侵蚀研究还鲜见报道。该研究以陕北2017年"7·26"特大暴雨为例，研究了岔巴沟流域3种土地利用类型（农地、休闲地和撂荒地）坡面切沟发育形态特征及体积估算模型。结果表明：1）农地、休闲地和撂荒地切沟长度分布在20 m内的占比分别为55.6%、34.8%和44.8%；农地切沟平均深度为110 cm，分别比休闲地和撂荒地高18.3%、19.2%；农地和休闲地切沟平均宽深比分别为0.87和0.84，横断面呈"宽-浅型"，而撂荒地切沟呈"方型"（宽深比1.01）。2）撂荒地切沟侵蚀体积分别比农地和休闲地减少47.8%和28.3%，表明植被恢复有效地削弱了极端暴雨作用下的切沟侵蚀。3）农地切沟不同坡段侵蚀体积由高到低为下坡、上坡、中坡，而休闲地和撂荒地切沟侵蚀体积沿坡长方向呈递增趋势；3种土地利用类型切沟在上坡段的沟岸拓宽速率大于下切速率，中下坡则相反。4）农地、休闲地和撂荒地切沟侵蚀体积均与切沟长度、横断面面积呈极显著幂函数关系（P<0.001），横断面面积是切沟体积估算更为有效的参数。研究结果可为黄土高原丘陵沟壑区不同土地利用类型坡面切沟侵蚀体积估算及其防治提供重要依据。     

黄土高原王东沟小流域土壤表层全磷空间变异性

[目的]分析黄土高原小流域尺度地貌单元和土地利用方式对全磷空间变异的影响。确定全磷空间变异规律,为土地生产管理提供依据。[方法]2017年在黄土高原王东沟小流域采集了267个土壤样品进行了分析研究。用ArcGIS 10.2绘制样点全磷含量分布图,并分析不同地貌单元和土地利用方式下全磷含量的空间变异特征。[结果]王东沟小流域内不同地貌单元下全磷含量呈现为:塬面塬坡沟道的分布规律,不同的土地利用类型下全磷含量呈现为:果园废弃果园农地林地荒草地的分布规律。土地利用方式和地貌单元对全磷空间变异的影响极为显著。[结论]王东沟小流域土地利用经过14a不断的调整,全磷的含量有了大幅度的提升,土地利用方式显著改变了全磷的空间分布格局。     

基于MIKE模型的不同淤地坝型组合情景对小流域侵蚀动力和输沙量的影响

为科学认识黄土高原淤地坝建设对小流域侵蚀动力过程的影响,通过分布式水文模型MIKE SHE和一维水动力模型MIKE 11耦合模拟了王茂沟流域的洪水过程,并计算了流域主沟不同断面的侵蚀动力参数和不同坝型组合的减沙效益。结果表明:(1)淤地坝建设在不同程度上改变了小流域沟道的侵蚀动力分布,坝系建成后沟道的侵蚀动力参数减幅最大。(2)小流域径流侵蚀功率表现为上中游剧烈变化,下游趋于稳定,且下游径流侵蚀功率明显小于上游。(3)淤地坝建设可以有效减小流域的输沙量,其中单独建设骨干、中型、小型坝相比流域未建坝时,输沙模数分别减少24.74%,47.11%,64.11%;流域坝系建成后减沙效益最明显,流域输沙量减少83.92%。研究成果可为黄土高原淤地坝减沙效益评估提供科学参考。     

近20年黄土高原土地利用/植被覆盖变化特征及其成因

近20年来黄土高原土地利用结构发生了显著变化,生态环境取得举世瞩目的成效,揭示黄土高原土地利用/植被覆盖的变化特征及成因是促进区域土地利用空间格局优化的关键环节。运用GIS技术,对黄土高原2000年、2010年、2017年的遥感数据进行转移矩阵,系统研究了黄土高原地区、各省份和生态分区的变化特征并定性分析了造成土地利用变化的主要驱动力。结果表明:(1)黄土高原的土地利用类型主要以耕地、林地和草地为主。2000—2017年,耕地、草地和水域都呈减少趋势,分别减少了19 860 km²,42 410 km²和5 630 km²,而林地、建设用地和未利用地面积有所增加,分别增加30 650 km²,9 840 km²和27 410 km²; 草地为主要的转入与转出地类。(2)黄土高原7个省域主要土地类型与区域变化规律基本一致,草地退化现象在甘肃与宁夏的不同时间段较为严重,具体表现为甘肃省在2000—2010年有7 830 km²草地转移为未利用地,而宁夏在2010—2017年有8 900 km²草地转入未利用地。(3)黄土高塬沟壑区、丘陵沟壑区、风沙区和土石山区及河谷平原区水域面积持续减少,且主要向耕地发生转移,同时生态区土地利用主要发生耕地与草地的动态转移。(4)黄土高原土地利用/植被覆盖变化的主要驱动力是社会因子,其中国家政策的实施起到了主要的推动和促进作用,而人类活动等的影响相对小一些。研究结果可为黄土高原土地利用规划提供科学支撑,也为该区生态保护和高质量发展服务。     

陕北黄土区切沟内陷穴发育特征研究

通过对陕北黄土区六道沟流域切沟内陷穴进行调查和统计分析,探讨陷穴的发育形态、分布特征及其影响因素。结果表明:陷穴主要形态类型有浅碟形、漏斗形、竖井状以及葫芦状,其断面形态接近于圆形或椭圆形,陷穴形态规模相对偏小,处于成熟期的比例达72.8%;陷穴长轴、短轴和深度之间均具有显著的正相关关系,其中长轴和短轴之间具有较好的线性关系;陷穴主要发育位置为沟道中上游及支沟沟头处,多与黄土桥和跌水相伴而生;切沟为陷穴的形成和发育提供了条件,陷穴亦可促进切沟进一步发展。上述结果完善了切沟内陷穴的发生发育特点和空间分布特征,为深入认识陷穴侵蚀机理,进一步研究陷穴侵蚀状况、陷穴与切沟发育关系以及区域水土流失治理等提供了基础。     

暴雨条件下黄土高原长陡坡耕地细沟侵蚀特征

细沟侵蚀是黄土高原坡耕地侵蚀的重要形式,已有的坡耕地细沟研究成果主要集中于裸土细沟侵蚀特征及其发生机制的小区或水槽试验,极端暴雨条件下细沟侵蚀观测很少,对于田间细沟侵蚀的预报仍然缺少实用的方法。该研究旨在基于实地测量数据,分析极端暴雨条件下黄土高原长陡坡耕地细沟侵蚀特征。2017年无定河流域"7.26"暴雨后,在暴雨中心附近选择15个样地进行细沟侵蚀测量,结合Google影像和无人机航拍影像,调查分析不同类型坡耕地细沟侵蚀特征。研究结果表明:1)暴雨条件下坡耕地是坡面细沟侵蚀的主要地类。长陡坡裸露坡耕地侵蚀模数为22 478 t/km~2,坡长20 m时,约为种植作物的平作坡耕地的1.5倍,且坡长越长,差异越大。等高垄作能够有效减少细沟侵蚀。撂荒坡耕地细沟侵蚀模数仅为裸露坡耕地的12%,草地未发现细沟侵蚀。2)幂函数可以很好地拟合细沟侵蚀模数与坡长的关系,裸露坡耕地坡长指数为0.831。暴雨强度越大,坡长指数也越大。3)上坡来水在坡耕地上造成了严重的细沟侵蚀,与地块上承接的降水相比,上坡来水对细沟侵蚀的影响更大。坡面浅沟汇水明显减小细沟侵蚀强度,浅沟发育程度越高,细沟侵蚀强度就越小。该研究可为估算黄土高原丘陵区不同类型坡耕地的细沟侵蚀提供重要参考。     

甘肃省小陇山林区降雨型浅层滑坡汇水面积-坡度关系及其影响因素

沟头发生侵蚀的地形临界模型可有效预测侵蚀沟的形成条件,浅层滑坡失稳形成的洼地也是沟头形成的方式之一。为探究浅层滑坡临界起动模型的特点,以甘肃省天水市小陇山林区的降雨型浅层滑坡为研究对象,运用汇水面积-坡度关系,构建临界起动模型,与黄土高原典型侵蚀沟(浅沟、切沟)的临界起动模型进行对比分析,并探讨土地利用类型、植被类型和土壤质地对该模型的影响。结果表明:(1)浅层滑坡临界起动模型为S=3.50A_s^-0.34,其侵蚀阈值为3.50,大于黄土高原典型浅沟(0.96)和切沟(1.54)的侵蚀阈值。研究区浅层滑坡一般发生于土层较薄的陡坡地带,其平均坡度(S=1.26)大于浅沟(S=0.35)与切沟(S=0.46),单位汇水面积(A=89.08 m²/m)小于浅沟(A=920.93 m²/m)和切沟(A=1 129.82 m²/m)。(2)汇水面积与坡度平方的乘积(AS²)代表了沟头产生侵蚀的能量指标值。研究区浅层滑坡AS²值在269.1~5 703.2 m²,平均值为1 772.97 m²,黄土高原浅沟AS²值在4.74~892.66 m²,切沟在41~814 m²,启动能量值方面,浅沟<切沟<浅层滑坡。(3)土地利用类型、植被类型和土壤质地通过影响土壤的抗冲力、渗透性和黏粒含量,从而对浅层滑坡的起动难易程度产生影响。在不同的土地利用方式中,农地最易发生侵蚀,其次是林地。油松林附近浅层滑坡的抗侵蚀能力高于日本落叶松林。研究结果为探究浅层滑坡的起动条件提供理论依据。     

Ephemeral gully recognition and accuracy evaluation using deep learning in the hilly and gully region of the Loess Plateau in China

Ephemeral gullies are widely distributed in the hilly and gully region of the Loess Plateau and play a unique role in the slope gully erosion system. Rapid and accurate identification of ephemeral gullies impacts the distribution law and development trend of soil erosion on the Loess Plateau. Deep learning algorithms can quickly and accurately process large data samples that recognize ephemeral gullies from remote sensing images. Here, we investigated ephemeral gullies in the Zhoutungou watershed in the hilly and gully region of the Loess Plateau in China using satellite and unmanned aerial vehicle images and combined a deep learning image semantic segmentation model to realize automatic recognition and feature extraction. Using Accuracy, Precision, Recall, F1value, and AUC, we compared the ephemeral gully recognition results and accuracy evaluation of U-Net, R2U-Net, and SegNet image semantic segmentation models. The SegNet model was ranked first, followed by the R2U-Net and U-Net models, for ephemeral gully recognition in the hilly and gully region of the Loess Plateau. The ephemeral gully length and width between predicted and measured values had RMSE values of 6.78 m and 0.50 m, respectively, indicating that the model has an excellent recognition effect. This study identified a fast and accurate method for ephemeral gully recognition in the hilly and gully region of the Loess Plateau based on remote sensing images to provide an academic reference and practical guidance for soil erosion monitoring and slope and gully management in the Loess Plateau region.     

基于最大似然法与矩法的黄土高原小流域土壤碳氮空间变异分析

以黄土高原寺底沟小流域为研究对象,根据不同土地利用方式采集46个样点的土壤样品,通过地统计方法对土壤有机碳和全氮的空间变异特征进行了分析。采用受限最大似然法(REML)和矩法(MOM)两种方法分别对变异函数进行了估计,通过交叉检验选择克里金预测效果较好的变异函数进行地统计插值。(1)与矩法(MOM)相比,在多数情况下受限最大似然法(REML)估计的变异函数进行克里金插值更加准确。(2)土层深度对土壤全氮空间变异影响较小,对土壤有机碳影响较大,表层土壤有机碳含量及变异程度明显高于下层土壤。(3)土地利用方式对土壤有机碳和全氮的空间分布有重要影响,灌木林和天然草地土壤有机碳和全氮水平最高,弃耕地其次,梯田、果园、人工草地最低,表明退耕还林对提高土壤碳氮水平有重要贡献。