西南丘陵区坡地土壤水分时空变异的地统计学分析.docxVIP

西南丘陵区坡地土壤水分时空变异的地统计学分析 地球物理的差异。经典的土地球理统计的显著特征是强调变量因子的空间过程，考虑其空间分布特征和空间自我。此外，还证明，这是分析土壤水分空间特征及其变化规律最有效的方法之一。国外研究者利用遥感和野外实测的土壤水分资料,从不同尺度上分析了土壤水分的地统计特性,但他们的研究多集中在表层土壤水分的空间结构分析,对剖面土壤水分的空间结构特性研究不够,缺乏系统分析实测土壤水分空间结构的季节变化规律。国内尤其是在西南丘陵地区多是利用经典统计学方法或描述性方法研究土壤水分,用地统计方法对其时空变异特征的研究鲜有报道。利用地统计学研究西南丘陵地区平均和剖面土壤水分空间结构及其季节变化特征,揭示该尺度下土壤水分的时空变异规律,研究结果可以有效地指导土壤水分的取样设计及取样密度的确定,或进行空间内插和制图;增强对土壤水分空间格局在景观尺度上如何影响生态过程的理解;为土壤水分的有效利用与管理提供科学依据。 1 土温和降雨特征 本研究在重庆市北碚区西南农业大学试验场进行。该场地处东经106°26′,北纬30°26′,海拔230m,属亚热带湿润气候。多年平均温度18.3℃,年际间平均气温变化不大。最热月(7月)平均28.7℃,最冷月(1月)平均7.4℃,年≥0℃的总积温为6679.5℃,≥10℃的总积温为5960℃。土温仅在1月表层才低于9℃。其余均高于10℃。雨量在季节分配上极不均匀,多集中在5~9月,占年降雨量的70.0%,其中5、6、9月雨量最大,7、8月相对减少,常出现伏旱。多年平均降水1105.4mm,且具有雨日多且多秋雨和夜雨的特点,历年平均雨日161天。一年中以夏秋季雨日最多,其中:中雨(10.0~24.9mm)与大雨(25.0~49.9mm)的降雨日数以5~9月最多,达31.9天。此期常带有暴雨出现。夜雨占全年雨量的73.3%,空气湿度大,多年平均相对湿度为80%,夏季绝对湿度大。云雾多,辐射和光照不足,多年平均87108.0卡/cm2年。试验研究地为丘陵旱坡地,土壤为侏罗纪沙溪庙组紫色沙岩发育的灰棕紫泥土。 2 学习方法 2.1 土壤含水量测定 采用网格取样方式(图1),网格大小为55m,共设置30个取样点。土壤含水量用烘干法测定。在采样点从2000年1月到12月,每月15日测定一次,每个样点分6层(5、10、20、30、40和50cm)测定,每层取3个数值,取其平均数作为该层的含水量。 2.2 半方差变异函数模型 地统计分析用GS+软件,关于地统计学方法的基本原理和方法在许多文章中都有很好的描述,在此不再叙述。变异函数的计算是通过GS+进行的,一般要求数据符合正态分布,否则可能存在比例效应。若不成正态分布则须进行对数转换。半方差图描述了随机函数Z空间依赖性的组成部分,将给定偏离间距为h的半方差估计为所有由间距h隔开的观察值之间的方差平均值,下式为半方差变异函数模型: r(h)=12N(h)∑i=1n[Z(xi)?Z(xi+h)]2r(h)=12Ν(h)∑i=1n[Ζ(xi)-Ζ(xi+h)]2 式中r(h)为变异函数;h为样点空间间隔距离,称为步长;N(h)为间隔距离为h观察值的对数;Z(xi)和Z(xi+h)分别是区域化变量Z(x)在空间位置xi和xi+h的实测值。 本试验主要用球状、指数和纯块金效应模型进行拟合。一些地统计学家认为最好避免用高斯模型,因为使用高斯模型常常会导致插值不准确并且会加大插值图的人为影响。目前对套合模型(图2)的选择和其参数(变程、基台值)的计算仍具有争议。 3 结果与分析 3.10 土壤水分随气候的变化 表1是 0~50cm平均土壤水分变异函数理论模型得出的相应参数。理论变异函数可用球状模型、指数模型和纯块金效应模型拟合,决定系数变化在0.41~0.99之间,F检验为显著水平(表1),说明理论变异函数模型很好地反映了土壤水分的空间结构特性。由于土壤水分在整个观测期的大部分时间表现出基台值,反映出土壤水分在研究区域内具有平稳特性或近平稳特性。平均土壤水分空间结构的季节变化如图3所示,1~6月土壤水分含量保持在24%上下,7月下降到14.77%,7~9月变化不大,10月上升后在11月稍有下降,说明0~50cm平均土壤水分主要受气候影响。 研究结果表明土壤水分具有强烈和中等程度的空间自相关(表1),这意味着土壤水分在此尺度上具有明显的空间自相关和格局。基台值即全部样点的总变异,其变化幅度0.83~4.98,表现出明显的季节变化趋势。这种变化趋势与控制土壤水分格局的因子如土地利用、地形和降雨量等有关。1~6月期间,土壤水分的空间变异程度低且变化不大,这可能是因为土壤水分的空间格局由土地利用和地形控制。进一步观测图3可以发现,基台值的季节变化趋势与平均土壤含水量相反。此外