苏北海涂围垦区土壤水分空间变异性及其协同克立格估值

以苏北滨海滩涂围垦区为研究区域,应用地统计学方法研究了剖面不同深度土壤含水量的空间变异特征,以0～20 cm表土层含水量作为协同变量,分析了不同深度土壤含水量的交互半方差特征,探讨了土壤含水量的普通克立格和协同克立格估值方法.结果表明:各层土壤含水量都表现为弱变异强度;受结构性因素和随机性因素的共同作用,各土层含水量均呈中等空间自相关性;协同区域化变量间均表现为正相关,且协同区域化变量的空间结构性要优于单一变量;协同克立格对土壤水分局部变异的描述较普通克立格更为详细;与普通克立格法相比,协同克立格法估值产生的均方误差减小10.1%～21.3%,平均标准误差减小11.3%～17.7%,预测值和实测值间的相关系数提高15.9%～26.4%.该研究为滨海滩涂地区土壤水分分区管理和水盐优化调控提供了一定的理论参考.     

黑土区小流域土壤速效磷空间分布模拟方法

为探索东北黑土小流域土壤AP空间分布模拟的最佳方法,综合考虑区域尺度、景观格局和采样方法等信息,选取两个典型黑土小流域(光荣和海沟河流域)作为研究对象,结合半变异分析、回归分析、主成分分析等方法,比较了不同空间插值方法(反距离权重法、径向基函数法、普通克里格、协同克里格、多元线性回归模型、地理加权回归模型、回归克里格和地理加权回归克里格)对AP空间模拟精度的影响。结果表明:(1)引入主成分分析后能够提高回归克里格和地理加权回归克里格方法的模拟精度(4.5%和2.4%);(2)地理加权回归克里格方法可以作为最优空间插值方法模拟黑土区小流域AP的空间分布格局。地理加权回归克里格方法相较于传统插值方法能在一定程度上提高黑土小流域土壤AP空间模拟精度,为小流域尺度黑土养分管理提供技术支撑。     

土壤速效钾养分含量空间插值方法比较研究

土壤养分连续空间分布数据是土壤信息系统工作的基础，土壤养分空间插值的研究因此变得尤为重要。对陕西省周至县北部猕猴桃适生区土壤进行采样，以对猕猴桃生长作用较为密切的土壤速效钾含量为研究对象，用普通克里格（OK）、样条函数（Spline）、趋势面拟合（TSA）、距离权重反比法（IDW）等常用插值方法对采样点进行插值获取土壤速效钾空间分布图，进行交叉验证，结果表明能够反映出结构性影响的克里格插值方法明显优于其它方法，其中又以球形模型为最佳，样条函数、距离权重反比法在采样点密集区也能内插出较好的效果，但其受采样点密度影响较大．在采样点稀疏的地区内插结果较差。     

基于磁感式土壤表观电导率空间变异性的插值方法比较

针对目前黄河三角洲地区存在的土壤盐渍障碍问题，以该地区典型地块为研究对象，运用磁感大地电导率仪(EM38)，结合GIS与地统计学研究了不同样点密度下土壤表观电导率的空间变异特征，确定了最佳的空间插值模式，并采用偏差指数对各分布图层的空间相似性进行了评价。结果表明：不同样点密度下的土壤表观电导率均呈中等变异强度，并服从对数正态分布；各样点密度下的土壤表观电导率均表现为强空间相关性，其空间变异主要表现在小于10 m的田间尺度上，且对于预测精度，泛克立格>普通克立格>反距离权重>局部多项式。偏差指数法分析表明，各插值方法分布图的空间相似性均随样点密度的降低而下降；对于相同的样点密度减小比例下的空间信息保留度，泛克立格>普通克立格>反距离权重>局部多项式，采用泛克立格法可以在保证预测精度的基础上合理降低样点密度。该研究为黄河三角洲地区盐渍土地磁感式田间数据采集的合理密度确定以及优化空间插值模式选取提供了理论依据与技术参考。     

考虑数据变换的泾河流域月降雨空间插值

对泾河流域各站点1951-2005年多年平均月降雨量数据进行偏斜度分析.并分别做平方根、立方根、对数变换;根据泾河流域自然地理特征选择普通克里格法和梯度距离平方反比法对变换前后的月均降雨量数据进行空间插值,通过交叉验证的检验方法得到以下结论:偏斜度分析对于空间插值意义不大;两种插值方法都是对数变换后的数据插值结果较优,同时考虑高程的梯度距离平方反比法无论数据变换与否其插值精度都没有普通克里格法的精度高.选择对数变换和普通克里格法进行泾河流域多年平均月降雨量插值效果较好,为日降水插值的获得提供了可靠的本底文件.     

基于混合地理加权回归与克里格的区域降水量空间插值方法

基于四川省区域范围内144个气象站点的实测降水数据,在综合考虑空间位置、地形等影响因素的基础上,采用改进的回归克里格模型,即混合地理加权回归克里格模型（MGWRK）对四川省年降水量的空间分布进行空间插值,并与普通克里格（OK）、全局回归克里格（GRK）和地理加权回归克里格（GWRK）等模型的插值效果进行对比分析。结果表明：（1）应用逐步回归法筛选确定的用于回归分析的影响因子组合为经度、纬度和坡度,可有效消除解释变量间的多重共线性,为后续的空间插值奠定基础;（2）同一回归变量在地理加权回归（GWR）与全局回归（GR）两种回归模型中的AICc（修正的赤池信息量准则,Corrected Akaike Information Criterion）值之差（ΔAICc）可用于定量判定各回归变量的空间非平稳性类型,据此将变量坡度设为全局变量,经度和纬度设为局部变量进行处理。在此基础上,通过MGWRK模型对四川省年降水量进行空间插值;（3）MGWRK插值模型综合考虑了空间位置、地形等多个影响因素及其与降水相互关系的空间非平稳性特征,相对于传统的OK和GRK法具有更高的插值精度。     

电磁感应仪用于土壤盐分空间分布的协同克立格估值研究

以黄河三角洲地区典型地块为研究对象,应用地统计学的半方差函数理论,分析了0～40 cm和40～80 cm土壤盐分的空间变异特征.在对协同区域化变量进行交互半方差分析的基础上,采用协同克立格法,以电磁感应仪EM38测得的土壤表观电导率作为协同变量,对各层次的土壤盐分进行估值.结果表明:受结构性因素和随机性因素共同作用,各层次土壤盐分均表现中等强度的变异和空间相关性;随机性因素是引起土壤盐分空间异质性的关键原因;协同区域化变量间均表现为正相关,且协同区域化变量的空间结构优于单一变量.同普通克立格法相比,协同克立格法估值产生的均方误差减小13.1%～17.8%,平均标准误差减小5.83%～17.6%,预测值和实测值间的决定系数提高37.6%～42.6%.该研究结果为黄河三角洲地区土壤盐渍化的精准定量化、科学管理与合理改良提供了一定的理论基础和实践依据.     

县域耕地土壤锌含量的协同克里格插值及采样数量优化

耕地土壤锌含量的空间分布和适宜的采样数量研究,为农产品产地评价和认证、土壤环境质量调查提供帮助。以土壤有机质、pH、全氮、碱解氮、有效铁、有效铜、有效锰、高程数据为辅助变量,采用协同克里格插值方法,对吉林省舒兰市耕地土壤锌含量进行了插值精度分析,并对不同样点数量下协同克里格法在县域尺度耕地土壤锌含量空间分布研究中的适用性进行了评价。结果表明,相同采样数量下,土壤锌协同克里格法的均方根误差比普通克里格法低,拟合效果更好;辅助因子对土壤锌协同克里格插值的影响效果大小为:全氮>有效锰>高程>有效铁>碱解氮>pH>有机质>有效铜;当高程和有效铜作为辅助因子时,土壤采样数量可以减少20%,土壤锌协同克里格插值仍能满足精度要求。协同克里格作为一种更为精确和经济的空间插值方法,可为县域尺度耕地土壤重金属含量的空间分布研究提供更多的信息和帮助     

克里格插值在区域土壤有机碳空间预测中的应用

基于河北省第二次全国土壤普查数据,对比了常用土壤有机碳相关因子土地利用和土壤类型与普通克里格插值结合前后对土壤有机碳密度空间预测精度的差异,探讨了普通克里格插值在区域土壤有机碳空间预测中的应用。研究结果表明,土地利用能够独立解释土壤有机碳密度总方差的19.0%,与普通克里格插值结合以后能够将对土壤有机碳密度总方差的解释程度显著提高到30.2%。低级土壤分类土属能够独立解释土壤有机碳密度总方差的45.0%,但与普通克里格插值结合以后对土壤有机碳密度总方差的解释程度为44.8%,两者相差不大。因此区域空间上能否进一步应用普通克里格插值优化土壤有机碳的空间预测与所选用的土壤有机碳相关因子有关。     

海涂土壤剖面电导率的协同克立格法估值及不同取样数目的比较研究

应用地统计学的半方差函数和普通克立格法 ,分析了 2 0m× 2 0m的采样间隔上 ,同一剖面0～ 2 0、2 0～ 4 0、4 0～ 6 0、6 0～ 10 0cm四个层次上海涂土壤电导率的空间相关性和空间分布。并在对协同区域化变量进行交互半方差分析的基础上 ,用协同克立格法 ,将 0～ 2 0cm层上的土壤样品作为辅助变量 ,其他三层上初始及减少后的土壤样品作为主变量 ,来对 2 0～ 4 0、4 0～ 6 0、6 0～ 10 0cm土层上的电导率进行估值。同普通克立格法相比 ,协同克立格法在这三层估值所产生的均方根误差 (RMSE)减少了 9%到 2 5 1% ,平均标准误差(ASE)减少了 9 4 %到 2 1 9% ,预测值与实测值之间的相关系数协同克立格法比普通克立格法提高了 4 6 %到85 %。数据分析的结果表明 ,利用协同克立格法 ,深层的样品数目在被减少至 6 0 %的情况下仍能满足精度要求 ,没有明显地导致空间信息的丢失。因此 ,协同克立格作为一种更为精确和经济的方法来对样品少的或难以测试的变量进行估值不但可以大大提高估测精度 ,而且可以减少采样数目 ,节省采样的成本。     

基于多元回归的黄土沟壑区小流域土壤水分空间模拟

以延河流域羊圈沟为研究对象,基于土壤水分空间自相关理论,在土壤水分环境影响因子多元回归基础上,建立了土壤水分空间化模型,并探讨了土壤水分空间分布特征。结果表明:研究区内土壤水分空间分布具有正相关性,且相关性显著,在空间上表现为聚集状态;表层土壤水分与植被覆盖度、高程、坡度、地表粗糙度、地形起伏度和地表切割程度等因素密切相关,且与地表切割程度呈正相关关系,而与植被覆盖度、高程、坡度、地表粗糙度以及地形起伏度因子呈负相关;与普通克里格(OK)、反距离加权(IDW)、径向基函数(RBF)和全局多项式(GPI)的土壤水分空间内插方法相比,多元回归+普通克里格方法(MROK)的拟合程度最高,是表层土壤水分空间化较好模型;土壤水分空间分布格局明显,主要表现为西部和东北部土壤水分含量低,东南部土壤水分含量高。该研究揭示了土壤水分环境因子相关性及空间分异规律,对区域生态文明建设具有一定的促进意义。     

采样数量与空间插值方法对华北平原典型区土壤质量评价空间预测精度的影响

空间插值方法的差异以及采样数量不同对土壤质量评价空间分布的预测精度会产生影响。本文通过随机抽样的方法,以山东省禹城市的土壤质量指数为例,从359个土壤采样点中抽取了340,170,90,50,30五个样本子集,通过普通克里格、简单克里格、反距离加权法和样条函数法4种插值方法,分别对其空间变异和布局进行了解析和预测。结果表明：不同插值方法对预测精度影响不显著,而采样点数量则显著影响了土壤质量指数空间分布预测的精度。本文提出在华北平原县域尺度上,以土壤质量评价作为调查目的的土壤采样中,90个样点是比较适宜的采样数量。同时,将我们的结果与Cochran最佳采样数量计算公式获取的最佳采样量比较后发现,Cochran方法获取的最佳采样数量明显偏低,若不考虑实际的空间变异情况,仅仅使用Cochran公式可能会导致土壤质量空间预测不准确。     

降水量空间插值方法在小样本区域的比较研究

针对降水量观测点数量少的小样本区域,寻求适合本区域的降水量空间插值方法,采用已知观测点降水量数据对待估点降水量进行预测,得到可信度较高的降水量预测结果,对指导当地生产生活具有十分重要的意义。在分析美国Upper Sangamon流域9个站点的72个月月降水量数据的基础上,比较了基于ArcGIS的普通克里金法、反距离权重法、样条函数法、趋势面法4种空间插值方法的插值结果与观测数据,分析了降水量空间插值方法对内插点和外推点的插值精度及应用不同插值方法时观测点数量对插值精度的影响。结果表明:反距离权重法较优于其他3种空间插值方法。空间插值精度随观测点数量变化而变化,对于按照距待估点距离由远及近减少观测点的选点方法,反距离权重法插值精度较高,受观测点数量影响较小,样条函数法与趋势面法受观测点数量影响较大,普通克里金法处于两者之间。     

基于数字高程模型的复杂地形下的黑龙江平均气温空间插值

空间化的气象资源被广泛应用于生态系统模拟与农业区划中,但常规插值模型对地形复杂、气象站点稀少的地区气象要素空间分布模拟的精度不高。本文在定量分析海拔高度、经纬度等因子对气温空间分布影响的基础上,选取黑龙江1957-2004年1月、4月、7月和10月的平均气温数据,利用DEM模型和相关辅助信息进行复杂地形条件下的气象资源空间分布模拟。精度交叉验证的结果表明,基于DEM辅助相关信息的插值精度明显高于常规的逆距离权重法、样条函数法、普通克里格法,插值精度提高了0.5-0.9℃。气温的空间分布趋势也更加合理,更好的体现出气温分布的空间异质性。     

非线性克里格方法在日降水量分布估算中的应用

为研究非正态分布条件下的日降水量空间分布,选择湖北省75个气象观测站2010年7月的降水过程资料,通过正态变换和变差函数分析,引入多元高斯克里格方法进行逐日降水量空间分布估算,并与其它方法进行对比分析.结果表明：（1）通过合理的各向异性参数设置及理论变差函数选择,能够充分凸显降水量分布特征.（2）对比发现,反距离加权法（IDW）、普通克里格法（OK）和多元高斯克里格法（MGK）3种估算方法中,MGK估算的平均误差最小且误差曲线最平稳,是该地区非正态分布日降水量较适合的空间分布估算方法.（3）使用MGK方法,能有效改善IDW带来的“牛眼效应”以及OK法造成的过度平滑,得到具有较高精度的日降水量数据集.     

新疆地区年降水量的空间插值方法研究

为提高新疆地区年降水量的空间插值效果,通过增设模拟气象站点、补充实测站点数据空间变异性信息,采用确定性函数和地统计学的综合插值方法进行新疆地区年降水量的空间插值与比较验证分析。结果表明：（1）利用地统计方法分析多年平均降水量的最佳抽样间距,确定气象站点影响区,建立降水量与经纬度、高程、坡度、坡向、NDVI间的多元回归方程,增设模拟气象站点,可以有效增加降水量空间变异信息,显著提高降水量空间插值效果。（2）通过对比发现,反距离加权法（IDW）、复二次径向基函数法（MRBF）和普通克里格法（OK）3种空间插值方法中,MRBF的插值效果最好,是该地区降水空间分布的最适插值方法。（3）在新疆地区采用逐步综合的方法,即以多元回归方法增设模拟气象站点和MRBF进行年降水量空间插值,全面考虑影响因子与空间变异性,去除了多元回归带来的＂边界效应＂,弥补了实测降水量空间信息,得到了具有较高精度的年降水量数据集。     

基于无人机影像快速估算矿山排土场边坡土壤侵蚀速率的方法

[目的]探究无人机高清影像快速估算矿山排土场边坡土壤侵蚀速率的方法,为无长期地形监测地区的土壤侵蚀强度估算提供了新方法。[方法]通过无人机和RTK-GPS分别获取研究区影像及地面控制点坐标,利用Agisoft PhotoScan Professional 1.1.2软件生成研究区空间点云,利用ArcGIS 10.3生成高精度DEM,通过目视遥感解译和现场调查,快速估算研究区土壤侵蚀速率。[结果]本研究区排土场边坡土壤流失面积11 052.53m2,土壤侵蚀总体积835.41m3,土壤侵蚀速率为4 043.82t/(km2·a)。使用普通克里金法插值生成DEM数据精度为0.26m。[结论]通过本研究方法计算出的排土场土壤侵蚀强度属于中度侵蚀,土壤侵蚀速率较为符合实际情况,精度也满足了本研究的需要,同时,在高密度点存在的情况下,各种插值方法生成DEM数据误差相差不大,使用普通克里金法能够使结果最优。     

土壤盐分和含水量的空间变异性及其CoKriging估值——以黄河三角洲地区典型地块为例

以黄河三角洲地区典型地块为研究区，运用经典统计学和地统计学相结合的方法研究了不同深度土层盐分和含水量的空间变异特征。结果表明：研究区土壤盐分和含水量普遍较高，土壤表层积盐作用明显；各层土壤盐分均属于中等变异强度，土壤含水量除表层为中等变异强度外，其余各层为弱变异强度；受结构性因素和随机性因素的共同作用，各土层盐分和含水量均具有中等的空间相关性。对Kriging插值结果分析表明，研究区各土层盐分和含水量在一定的范围内均存在着空间上的相关性和变异性，表层土壤盐分和含水量的空间分布主要受微地形和气候条件的影响，地下水性质是主导深层土壤盐分及含水量空间分布的主要因素。分别采用表层土壤盐分及相应土层含水量对深层土壤盐分进行CoKriging估值，均能提高估值精度，采用表层土壤盐分对深层土壤盐分进行CoKriging估值的精度要明显高于采用土壤含水量，其估计方差减少百分数最大达到167．36％。     

站点密度及插值方法对ET0空间插值精度的影响

流域参考作物蒸发蒸腾量（ET0）插值方法的研究对流域尺度作物耗水时空变化规律有重要意义。该文通过海河流域162个国家农业气象站3a（2003－2005年）旬值气象资料,利用Penman-Monteith公式计算了这些站点ET0,采用ArcGIS软件中常用的Spline、IDW和Ordinary Kriging（OK）法,以及近些年研究较多的线性回归Regression插值法,对不同站点密度条件下的ET0进行空间插值。分析了各空间插值方法在不同站点密度条件下的优劣性,并且给出了本流域内各种站点密度范围条件下计算ET0最适宜的插值方法。结果表明以站点密度1.3个/万km2为界,当站点密度低于此密度时,推荐使用Regression插值法;当站点密度大于1.3个/万km2时,推荐使用IDW和OK插值法;当站点密度大于4.3个/万km2,以上三种插值法并无显著差别;不推荐使用Spline插值法。     

典型黑土区村级尺度侵蚀沟演变

侵蚀沟的演变是东北水土流失过程的重要方面.利用1968和2009年亚米级遥感影像,获取典型黑土区村级尺度的侵蚀沟分布、耕地垄向和土地利用数据,基于1∶1万地形图获取等高线、坡度、坡长和垄向坡度数据,在地理信息系统和遥感技术的支持下,分析典型黑土区41 a村级尺度侵蚀沟的动态变化及其影响因素.结果表明:1)41 a间,典型黑土区侵蚀沟密度有明显增加的趋势,基于分辨率为亚米级遥感影像解译的侵蚀沟密度较以往研究结果大幅增加;2)自然坡度0.25°～1.5°范围内为该区域侵蚀沟密度变化最大的区域,8°以上区域由于实施退耕还林措施,侵蚀沟密度降低;3)300～500 m坡长区域侵蚀沟密度出现极值;4)由于实施改垄措施,41 a后垄向坡度在0～0.25°范围内耕地面积增加,大于0.25°范围内耕地面积均减少.