关于土壤导气率流域尺度的空间分布特征研究,农业基础科学论文.docx

《关于土壤导气率流域尺度的空间分布特征研究,农业基础科学论文.docx》由会员分享，可在线阅读，更多相关《关于土壤导气率流域尺度的空间分布特征研究,农业基础科学论文.docx（10页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、关于土壤导气率流域尺度的空间分布特征研究,农业基础科学论文近年来人们对土壤与环境间互相作用以及精到准确调节农田土壤气体传导状况问题的关注度不断增加，数值模型研究土壤气体传输经过及其有效性与其对作物生长和生态环境的影响程度逐步成为现实。怎样发展简单可行的土壤气体传输动力学参数确定方式方法成为各国学者关注的热门问题。对此研究者在大量试验基础上，相继提出确定土壤气体传输动力学参数的不同方式方法，即直接测定方式方法和间接推求方式方法。两种方式方法各具优势，在生产实际中广泛应用。由于土壤气体传输经过比拟复杂，不仅与土壤质地、构造、含水率、容重、土壤有机质含量等土壤基本物理化学特征有关，而且遭到气候条件、

2、土地利用方式、生长植物群落的种类、田间耕作与管理方式等的强烈影响，有些方式方法难以准确估计田间土壤导气率，在很大程度上存在局限性，尤其难以在大范围区域研究土壤气体传输动力学特征参数。 除此之外，土壤气体传输存在严重空间变异性，其空间变异特征与诸多影响因素间存在何种关系，怎样准确获取一定范围内土壤导气率的平均值，此方面仍缺乏系统研究。综合国内外研究结果发现: 各研究者在其研究区域内作出相关研究，得出相关结论，这些研究区域各自的土壤物理化学性质不相一致，进而导致研究结果存在一定分歧。在同一个研究区域内，不同含水率、容重的条件下，土壤导气率的变异和空间格局也是不完全一样的。各个影响因素造成土壤导气率

3、变异，这些因素与变异之间是严密联络的。因而需要进一步研究土壤气体传输综合特征与各影响因素间的关系，说明土壤导气率空间变异特征，发展确定土壤导气率简单易行的估算方式方法，才能实现利用动力学模型分析不同尺度下土壤气体传输动力经过，为农业、水利和生态环境行业的生产实际提供有效的方式方法。 基于此，本文以陕西省泾惠渠灌区麦田土壤为研究对象，展开关于土壤导气率流域尺度的空间分布特征研究，并利用 Pearson 相关系数及交相互关关系分析各个影响因素(土壤物理基本参数) 与土壤导气率空间分布间的关系。 1、 材料与方式方法 1. 1 研究区大概情况 泾惠渠灌区位于陕西省关中平原，地理坐标介于东经 108

4、34 34 109 21 35 ，北纬 34 25 20 34 41 40 ，灌区北邻黄土台塬，西、南、东三面有泾河、渭河、石川河环绕，清峪河由灌区北部自西向东穿越，为一较完好的水文地质单元; 地面和地下径流排泄条件好，灌区东西长约 70 km，南北宽约 20 km，总面积 1 180 km2; 地势总趋势是自西北向东南倾斜，海拔高度350 450 m，地面坡降一般为1. 66% 3. 33% 。多年平均降水量为 512 mm，蒸发量1 212 mm，年平均气温 13. 6。 1. 2 取样与测量方式方法 本研究选取灌区麦田为测点，以 2. 5 km 的采样间距布设，除去 盲点 外，共 212

5、 个测点，图 1 为布置点示意图。采样点用手持 GPS 定位，记录采样点的经纬度，定位的样点经格式转换成 ArcGIS 能辨别的 . shp 格式。对研究区进行地图数字化，经投影转换后生成研究区样点分布图，录入采样点基本物理参数的信息后，即可用于研究区土壤物理参数及土壤导气率的变异分析。笔者曾利用 PL-300 型土壤导气率测定仪展开研究，详见文献1，5。 土壤物理基本参数的测定方式方法如下: 含水率(采集地表 0 10 cm 土层土样) 、容重通过枯燥法测量; 土壤颗粒组成利用马尔文激光分析仪测定。每个测点土样用密封袋带回，风干后研磨，过 100 目土壤筛，用于制备供试有机质含量(重铬酸钾外

6、加热法) 的土样。 1. 3 数据分析与处理 1. 3. 1 空间变异分析 利用 SPSS 19. 0 软件进行数据统计分析，以空间变异理论为基础，利用 ArcGIS 软件里的地统计模块(Geostatistical Analyst) ，采用空间克里格插值法(Ordinary Kriging) ，构建土壤物理基本参数与土壤导气率的空间变异分布图。克里格插值法详见文献。 1. 3. 2 交相互关系数确实定 变量 A 和 B，在一定滞后间距下，两个变量的交相互关系数可表示为 当 h =0 时，交相互关系数等于经典统计学计算得到的线性相关系数。通常用 t 检验验证交相互关系数的显著性，t 值计算公式

7、 式中 n 计算交相互关系数所用数据的对数当式(2) 计算得到的 t 值大于查表得到的值时，则以为该交相互关系数显著。 2、 结果与分析 2. 1 土壤物理基本参数及导气率描绘叙述性统计分析 土壤物理基本参数及导气率的描绘叙述性统计分析见表 1。用变异系数判定其变异程度: Cv 0. 1 属于弱变异性; 0. 1 Cv 1 为中等变异性; Cv 1 属于强变异性。由此看出研究区域容重为弱变异，有机质、饱和含水率、饱和度与土壤导气率属于中等变异。 采用地统计学分析的前提要求就是所有变量服从正态分布，即变量的大小分布符合自然状态。而土壤导气率服从对数正态分布，因而将其进行取自然对数的处理。通过表

8、1 中偏度、峰度以及 K-S 法进行正态检验(p = 0. 05) ，亦能看出土壤物理基本参数以及经过对数处理过的土壤导气率属于近似正态分布。 2. 2 土壤物理基本参数与导气率的空间变异及其分布格局 将与各测点土样对应的地理坐标经纬度，以及土壤物理基本参数、土壤导气率等相关数据录入ArcGIS 软件，在地统计模块的运行环境下，拟合得到各参数的变异函数值，详见表 2。采用变异函数理论模型拟合的方式方法，推荐其拟合最优值，即可得到结果: 土壤颗粒组成(黏粒、粉粒、砂粒) 和土壤有机质的最佳拟合模型为指数型; 容重、饱和含水率与饱和度的最佳拟合模型为球形模型; 土壤导气率的最佳拟合模型为指数型。表

9、 2 可见，各个变异函数的残差平方和均很小，接近于零，且决定系数 R2均大于 0. 8，讲明上述理论关系模型能够描绘叙述试验数组的变异函数，并且能够较好地反映出理论变异函数与试验数组的变异函数间的关系。 表 2 显示各个土壤物理基本参数与土壤导气率的分维数(D) 均大于 1. 9，讲明各土壤物理基本参数与土壤导气率的空间异质性较强; 容重、饱和含水率、饱和度与土壤导气率空间自相关(C/(C0+ C) ) 均在 0. 9 以上，讲明由自相关部分引起的空间异质性占总空间异质性的程度较大。也有学者用空间自相关的大小作为断定系统内变量空间相关性程度的根据，指出当 C/(C + C0) 为 0 0. 2

10、5、0. 25 0. 75 和0. 75 1 时，分别表示清楚变量具有较弱、中等和较强的空间依靠性。本文研究结果表示清楚: 土壤有机质的空间自相关为 0. 5，其空间依靠程度属于中等水平;容重、饱和含水率、饱和度与土壤导气率，具有较强程度的空间依靠性。 以最大相关距离作为确定采样间距的参考。由表 2 可见，土壤颗粒组成(黏粒、粉粒、砂粒) 最优拟合模型为指数型，且三者变程 A0的数值很接近，空间变异分析具备一致性，因而推荐颗粒组成的采样间距为变程的 3 倍，大约为 6 8 km。 容重、饱和含水率的最优拟合模型为球形模型，推荐变程作为采样间距的参考值，大约为 3 km。土壤导气率的最优拟合模型

11、为指数型模型，推荐 3 倍变程作为采样间距的参考值，大约为 7. 5 km。同时各参数最大相关距离都没有超出滞后距离(30. 23 km)的范围，符合逻辑。而试验设计的 2. 5 km 采样间距符合空间变异分析的采样要求，试验设计是合理的。 2. 3 Kringing 插值结果 以空间变异理论为基础，利用 ArcGIS 软件里的地统计模块，建立半方差函数模型(图 2) ，应用普通克里格法进行最优内插，构建土壤物理基本参数与土壤导气率的空间变异分布图(图 3) 。 总体上，各土壤物理基本参数的空间差异较大，在空间布局上呈现出不规律的斑块状分布。从图 3中能够看出: 颗粒组成黏粒、粉粒及砂粒含量分

12、布图具有相反趋势，砂粒含量相对较低的地方，黏粒、粉粒含量相对较高。容重的大小决定了土壤孔隙率的大小，容重越大孔隙率越小。因而容重与饱和度在空间分布上具有相反的趋势。土壤导气率与容重空间布局具有一定程度的一致性，与饱和度的空间布局具有相反趋势。容重越小，砂粒含量越高的地方，孔隙率越高，孔隙构造越发达，导气率数值也就越高。饱和度反映土壤孔隙含水的状况，饱和度小讲明奉献于气体传输的有效孔隙多，土壤导气率数值越高。这归因于土壤孔隙中气相与液相并存，土壤构造及其孔隙的几何构造特征决定了土壤气体传输的特性。 2. 4 土壤导气率与影响因素的相关性分析 利用 SPSS 软件对土壤导气率与主要影响因素之间作相

13、关性分析，主要影响因素包括土壤颗粒组成、有机质质量比、容重、饱和度，将 计 算 出 的Pearson 相关系数列于表 3。 由表 3 的分析结果来看，土壤导气率与有机质的 Pearson 相关系数绝对值介于 0 0. 5 之间，讲明两者之间微弱相关; 土壤导气率与容重的 Pearson 相关系数绝对值介于 0. 5 0. 8 之间，两者显著相关;土壤导气率与饱和度的 Pearson 相关系数绝对值介于 0. 8 0. 1 之间，两者高度相关。由于在地统计分析时，对于选取的尺度较小的研究区域来讲，DEM地图的精度也是导致土壤导气率与某些影响因素之间相关性较差的原因之一，也不能完全讲土壤颗粒组成与

14、导气率之间微弱相关。 以上利用 Pearson 相关系数分析各个影响因素与土壤导气率之间的相关性的结果也有一定道理。首先，土壤孔隙中气相与液相并存，饱和度高了，含水率相对就高，导气率自然就随之降低; 其次，容重决定了土壤孔隙率的大小，奉献于土壤通气能力的孔隙受制于容重的大小; 再次，在质地一样的情况下，土壤有机质含量高，土壤颗粒团圆体构造好，孔隙发达，导气率大; 最后是决定土壤质地划分的土壤颗粒组成。基于以上分析，在利用土壤物理基本参数推求土壤导气率的研究中优先考虑饱和度，其次依次考虑容重、有机质、土壤颗粒组成。 2. 5 土壤导气率与影响因素的交相互关性分析 在经典统计学中，Pearson

15、相关系数没有考虑空间位置，它仅代表两个变量在同一采样点上的相关性，交相互关图就能够反映出两个相关的变量在多大的滞后距离上存在相关性。 根据泾惠渠灌区土壤样本的实测数据，分析土壤导气率与各个土壤物理基本参数的交相互关关系，推求这些影响因素对土壤导气率的影响范围，各个影响因素与土壤导气率的交相互关图见图 4。由图 4 可见，各因素对土壤导气率的影响格局均不同;对于同一因从来讲，在不同方向上，交相互关系数也不同。对于容重和饱和度这两个重要的影响因从来讲，无论是正相关还是负相关，随着滞后间距的增大，交相互关系数的绝对值减小，讲明随着容重、饱和度与土壤导气率之间变量的间距增大，其相关程度随之降低。在 9

16、5% 的置信水平上，容重、饱和度对土壤导气率的影响范围均为 20 20 km。图 4也能讲明容重与饱和度是影响土壤导气率的主要影响因素。 3、 结论 (1) 说明了该灌区麦田土壤导气率与土壤物理基本参数的空间变异特征，土壤导气率分维数大于1. 9，反映出该研究区域土壤气体动力学参数具有较强的空间变异特征。 (2) 容重、饱和含水率、饱和度与土壤导气率的空间自相关(C/(C0+ C) ) 均在 0. 9 以上，讲明由自相关部分引起的空间异质性占总空间异质性的程度较大。推荐土壤导气率的采样间距的参考值大约为7. 5 km，颗粒组成采样间距为 6 8 km，容重、饱和含水率采样间距的参考值大约为 3 km。 (3) 土壤导气率与容重的 Pearson 相关系数绝对值为 0. 595，两者显著相关，其与饱和度的 Pearson相关系数绝对值为 0. 959，两者高度相关; 交相互关关系讲明容重、饱和度是影响土壤导气率的主成因素，在 95%的置信水平上，容重、饱和度对土壤导气率的影响范围均为 20 20 km。 以下为参考文献： 1. 王卫华，王全九，武向博，等 黑河中游绿洲麦田土壤水气热参数田间尺度空间分布特征J 农业工程学报，2020，29(9) :94 102