农用地质量因素与过程研究进展

1、农用地质量因素与过程研究进展农用地质量因素与过程研究进展2013年春季学术报告会年春季学术报告会创新良田工程科技创新良田工程科技 重构资源环境基础重构资源环境基础电话：电话O）；）M）Email: YFHUANGCAU.EDU.CN黄黄 元元 仿仿提提 纲纲 一、研究方向与内容一、研究方向与内容二、近期研究重点与远期规划二、近期研究重点与远期规划三、部分研究进展三、部分研究进展一、研究方向与内容一、研究方向与内容通过开展农用地质量形成通过开展农用地质量形成过程、演变规律及其驱动过程、演变规律及其驱动因素的应用基础研究，为因素的应用基础研究，

2、为建立农用地质量全面、高建立农用地质量全面、高效、持续提升的长效机制，效、持续提升的长效机制，提供关键理论支撑。提供关键理论支撑。农用地质量农用地质量因素与过程因素与过程农用地多功能农用地多功能诊断与评价诊断与评价农用地质量与产能调控农用地质量与产能调控农用地质量调查与监测农用地质量调查与监测基本农田质量保护与管理基本农田质量保护与管理一、研究方向与内容一、研究方向与内容研究农用地质量关键要素的物理、化学和研究农用地质量关键要素的物理、化学和生物过程的相互作用和耦合机制，揭示农生物过程的相互作用和耦合机制，揭示农用地质量形成过程，实现过程定量化；用地质量形成过程，实现过程定量化；建立农用地质量

3、要素的数学模型和评价方建立农用地质量要素的数学模型和评价方法，分析农用地质量及其关键要素的时间法，分析农用地质量及其关键要素的时间变化趋势和空间分布规律；变化趋势和空间分布规律；研究农用地质量因素、指标和等级变化的研究农用地质量因素、指标和等级变化的多尺度转换方法和模型；多尺度转换方法和模型；研究农用地质量演变的驱动因素，揭示农研究农用地质量演变的驱动因素，揭示农用地质量的演变趋势。用地质量的演变趋势。二、近期研究重点与远期规划二、近期研究重点与远期规划（一）近期研究重点（一）近期研究重点以以耕地耕地及其及其生产力生产力为重点，兼顾多功能，研究耕地质量因素为重点，兼顾多功能，研究耕地质量因素（

4、要素）。（要素）。因素因素地学因素地学因素气候因素气候因素工程（科技）因素工程（科技）因素生态（健康）因素生态（健康）因素耕地耕地农用地功能功能生产力生产力多功能多功能质量质量认定管理认定管理评价评价模拟模拟应用应用提升技术提升技术趋势分析趋势分析时空变化规律与时空变化规律与尺度转换研究尺度转换研究二、近期研究重点与远期规划二、近期研究重点与远期规划（一）近期研究重点（一）近期研究重点关键关键地学地学要素的物理、化学和生物过程的相互作用要素的物理、化学和生物过程的相互作用及其及其与作物生产力与作物生产力耦合机制，实现过程定量化耦合机制，实现过程定量化。其中，重点。其中，重点关注水及其驱动的水肥

5、耦合过程与耕地生产力。关注水及其驱动的水肥耦合过程与耕地生产力。关键要素的时间变化趋势关键要素的时间变化趋势与与空间分布规律空间分布规律，研究研究关键关键指指标和等级变化的标和等级变化的空间空间尺度转换尺度转换。研究耕地质量要素的综合方法，包括分等定级的研究耕地质量要素的综合方法，包括分等定级的评价方评价方法法与农田生产力模型。与农田生产力模型。选择耕地质量选择耕地质量关键要素关键要素，近期重点研究以下几方面：，近期重点研究以下几方面：共性基础研究（共性基础研究（1-3）二、近期研究重点与远期规划二、近期研究重点与远期规划（一）近期研究重点（一）近期研究重点全球气候变化与土地流转新形式下的耕地

6、全球气候变化与土地流转新形式下的耕地质量演变趋势质量演变趋势（宏观）（宏观）全球气候变化对耕地生产能力影响过程，全球气候变化全球气候变化对耕地生产能力影响过程，全球气候变化对中国耕地质量综合影响评估。对中国耕地质量综合影响评估。今年中央一号文件就提出加大土地流转力度，组建家庭今年中央一号文件就提出加大土地流转力度，组建家庭农场、农民合作社，把土地集中起来，两会多有关注（农场、农民合作社，把土地集中起来，两会多有关注（谢德体、秦博勇、黄志明等等）。规模化、集约化、现谢德体、秦博勇、黄志明等等）。规模化、集约化、现代化等新形式下的耕地质量演变？代化等新形式下的耕地质量演变？应用基础研究（应用基础研

7、究（4-6）二、近期研究重点与远期规划二、近期研究重点与远期规划（一）近期研究重点（一）近期研究重点中低等耕地质量形成过程与关键技术中低等耕地质量形成过程与关键技术低等别耕地的障碍因素与质量提升技术，特别是新开垦（低等别耕地的障碍因素与质量提升技术，特别是新开垦（含复垦）耕地的质量控制与快速提升技术。含复垦）耕地的质量控制与快速提升技术。中低等别（产）田占比近中低等别（产）田占比近70% ，是我国生产力的潜力所，是我国生产力的潜力所在。两会多有关注（谢德体、秦博勇等）。在。两会多有关注（谢德体、秦博勇等）。高标准基本农田创建中的质量控制。高标准基本农田创建中的质量控制。应用基础研究（应用基础研

8、究（4-6）（摘自胡存智，2012）二、近期研究重点与远期规划二、近期研究重点与远期规划（一）近期研究重点（一）近期研究重点高等别耕地的高产高效与质量定向培育关键技术高等别耕地的高产高效与质量定向培育关键技术作物高产高效的耕地质量（如土壤条件）是什么？作物高产高效的耕地质量（如土壤条件）是什么？如何定向培育耕地质量（如土壤结构与过程），为高产如何定向培育耕地质量（如土壤结构与过程），为高产高效创造土壤条件。高效创造土壤条件。高标准基本农田创建中的质量控制。高标准基本农田创建中的质量控制。应用基础研究（应用基础研究（4-6）二、近期研究重点与远期规划二、近期研究重点与远期规划（二）中远期研究规划

9、（二）中远期研究规划从耕地扩展到农用地。从耕地扩展到农用地。全面体现农用地多功能的质量要素与过程。全面体现农用地多功能的质量要素与过程。建设高标准基本农田，不能唯产能提升论，要把握农用地（耕地）质量的全面内涵，特别是耕地健康。实现农用地多功能的关键技术与集成。实现农用地多功能的关键技术与集成。把农用地视为完整的生态系统，维护农田生态系统循环过程及其与周边环境的联系，协调人地关系，实现生态文明。三、部分研究进展三、部分研究进展（一）土壤属性的时空分析（一）土壤属性的时空分析（二）土壤（二）土壤-作物系统过程模型与应用作物系统过程模型与应用（三）农田用水与水分生产力（李保国）（三）农田用水与水分生

10、产力（李保国）（四）高产高效的土壤条件与过程（任图生）（四）高产高效的土壤条件与过程（任图生）三、部分研究进展三、部分研究进展（一）土壤属性的时空分析（一）土壤属性的时空分析1、高程在土壤属性空间预测中的应用、高程在土壤属性空间预测中的应用2、不同外部趋势变量（组合）与有机质空间预测、不同外部趋势变量（组合）与有机质空间预测3、土壤质地空间预测方法的比较、土壤质地空间预测方法的比较4、土壤质地的空间三维模拟土壤质地的空间三维模拟5、 空间分析的尺度效应空间分析的尺度效应6、 土壤属性的时间演替土壤属性的时间演替.预测方法近五年共发表论文约近五年共发表论文约30篇篇（其中（其中SCI、EI收录收

11、录20篇）篇）研究区：北京市平谷区研究区：北京市平谷区目标变量：土壤有机质（目标变量：土壤有机质（SOM） 土壤速效钾（土壤速效钾（K） 土壤有效磷（土壤有效磷（P） 土壤有效锰（土壤有效锰（Mn） 土壤有效铁（土壤有效铁（Fe） 土壤有效锌（土壤有效锌（Zn） 土壤土壤pH值值样本点数：样本点数：887个训练点个训练点 222个验证点个验证点辅助变量：高程辅助变量：高程预测方法：预测方法：OK法、法、RK法、法、CK法法1、 高程在土壤属性空间预测中的应用高程在土壤属性空间预测中的应用三、部分研究进展三、部分研究进展p 结论结论高程作为辅助变量，可以提高土壤属性的空间预测精度，高程作为辅助变

12、量，可以提高土壤属性的空间预测精度，但并不适合所有土壤属性。但并不适合所有土壤属性。 相关因素分析表明，能否利用高程来辅助提高土壤属性的相关因素分析表明，能否利用高程来辅助提高土壤属性的空间预测精度与土壤属性的空间自相关度具有明显的相关空间预测精度与土壤属性的空间自相关度具有明显的相关关系。关系。 CK法和法和RK法在预测能力上法在预测能力上没有表现出明显的差别。没有表现出明显的差别。 三、部分研究进展三、部分研究进展高程辅助变量在不确定性分析中的应用高程辅助变量在不确定性分析中的应用 探讨在随机模拟计算环境下，是否可以利用高程探讨在随机模拟计算环境下，是否可以利用高程辅助变量来提高土壤属性空

13、间预测的准确性、局辅助变量来提高土壤属性空间预测的准确性、局部不确定性模拟的准确性和空间不确定模拟的准部不确定性模拟的准确性和空间不确定模拟的准确性。确性。局部不确定性局部不确定性：单一点上预测值大于或小于阈值的概：单一点上预测值大于或小于阈值的概率。率。空间不确定性空间不确定性：所有预测点上的预测值同时大于或小：所有预测点上的预测值同时大于或小于阈值的概率。于阈值的概率。三、部分研究进展三、部分研究进展研究样区：密云样区、平谷样区研究样区：密云样区、平谷样区密云样区：密云样区：284个样本点个样本点平谷样区：平谷样区：547个样本点个样本点目标变量：土壤有机质目标变量：土壤有机质辅助变量：高

14、程辅助变量：高程方法：序贯高斯模拟法（方法：序贯高斯模拟法（SGS） 序贯高斯协模拟法（序贯高斯协模拟法（SGCS）密云样区（a）和平谷样区（b）高程图研究区位置与采样点分布图结论结论通过序贯高斯协模拟方法利用高程辅助变量通过序贯高斯协模拟方法利用高程辅助变量： （1）可以提高土壤有机质预测结果的准确性可以提高土壤有机质预测结果的准确性 （2）可以提高土壤有机质空间不确定模拟的准确性）可以提高土壤有机质空间不确定模拟的准确性 （3）没能提高土壤有机质局部不确定性模拟的准确）没能提高土壤有机质局部不确定性模拟的准确性性三、部分研究进展三、部分研究进展2、 利用不同外部趋势变量（组合）利用不同外部

15、趋势变量（组合） 预测土壤有机质空间分布预测土壤有机质空间分布 探讨通过高程、坡度和地形湿度指数等地形指标之间探讨通过高程、坡度和地形湿度指数等地形指标之间的组合，是否可以进一步提高土壤有机质空间预测精的组合，是否可以进一步提高土壤有机质空间预测精度。度。 比较回归克里格法（比较回归克里格法（RK）和有限最大似然法为基础的）和有限最大似然法为基础的线性无偏最优估计法（线性无偏最优估计法（BLUP-REML）在不同外部趋）在不同外部趋势度情况下预测能力的表现。势度情况下预测能力的表现。三、部分研究进展三、部分研究进展研究区位置及采样点分布图研究区位置及采样点分布图研究区：北京市平谷研究区：北京市

16、平谷区区样本点数：样本点数：70个训练个训练点，点，131个验证点个验证点目标变量：土壤有机目标变量：土壤有机质质辅助变量：高程、坡辅助变量：高程、坡度、地形湿度指数度、地形湿度指数三、部分研究进展三、部分研究进展结论结论高程、坡度和地形湿度指数等地形指标可以作为辅助变量高程、坡度和地形湿度指数等地形指标可以作为辅助变量用来提高土壤有机质的空间预测精度，尤其是通过不同的用来提高土壤有机质的空间预测精度，尤其是通过不同的地形指标之间的组合，能进一步提高预测精度。地形指标之间的组合，能进一步提高预测精度。在样本点数比较小的情况下，在预测具有空间趋势的土壤在样本点数比较小的情况下，在预测具有空间趋势

17、的土壤属性方面，有限最大似然法为基础的线性无偏最优估计法属性方面，有限最大似然法为基础的线性无偏最优估计法（BLUP-REML）是一种比较可靠的预测方法）是一种比较可靠的预测方法 。三、部分研究进展三、部分研究进展以地形因素和分类变量组合辅助变量的土壤有机质空间预测以地形因素和分类变量组合辅助变量的土壤有机质空间预测地理位置与样点分布图地理位置与样点分布图坡长因子图坡长因子图土壤质地图土壤质地图三、部分研究进展三、部分研究进展地形因素地形因素基本属性基本属性复合属性复合属性分类变量分类变量土壤质地土壤质地土地利用土地利用土壤有机质土壤有机质辅助变量组合辅助变量组合1辅助变量组合辅助变量组合2土

18、壤有机质空间分布土壤有机质空间分布空间格局的环境相关分析空间格局的环境相关分析基于环境相关分析的空间预测基于环境相关分析的空间预测方差分方差分析析皮尔森皮尔森相关分相关分析析多元线多元线性逐步性逐步回归回归普通克普通克里格法里格法回归克回归克里格法里格法分类变量以及何种分类变量对土壤有机质空间预测精度的影响分类变量以及何种分类变量对土壤有机质空间预测精度的影响三、部分研究进展三、部分研究进展结论结论u地形因素、土壤质地和母质类型等自然、内在性质决定着土壤地形因素、土壤质地和母质类型等自然、内在性质决定着土壤有机质的空间分布有机质的空间分布。母质类型样本数（个）显著性检验结果（g kg-1）长石

19、岩类风化物12214.72bcd冲积物468.59e非碳酸盐3512.33d钙质岩类风化物3314.47bcd耕种型8416.25abc硅质岩类风化物3418.04a 洪积冲积物5614.32bcd黄土性母质2112.32d人工堆垫物2213.42dc铁镁质岩类风化物1616.94ab不同高程组有机质含量盒形图不同高程组有机质含量盒形图不同母质类型各颗粒组成均值两两比较结果不同母质类型各颗粒组成均值两两比较结果三、部分研究进展三、部分研究进展高程组（高程组（m）土壤有机质含量（土壤有机质含量（g kg-1）不同不同母质母质类型类型间土间土壤有壤有机质机质含量含量差异差异显著显著u将分类变量和地

20、形因素等将分类变量和地形因素等考虑到空间预测中考虑到空间预测中较大幅度较大幅度地提高土壤有机质空间预测地提高土壤有机质空间预测精度。精度。不同辅助变量组合下的实测、预测对比图不同辅助变量组合下的实测、预测对比图三、部分研究进展三、部分研究进展（样点编号）（样点编号）不同预测方法的不同预测方法的均方根误差均方根误差相对提高值（相对提高值（RI）M2变量组合下的回归克里变量组合下的回归克里格法格法RI值达值达37.48%辅助变量组合辅助变量组合1(M1):高程高程辅助变量组合辅助变量组合2(M2):高程、母质高程、母质回归预测和普通克里格模型平滑效应，回归预测和普通克里格模型平滑效应，趋趋中趋势中

21、趋势，须对残差再进行预测以降低预测，须对残差再进行预测以降低预测残差。在残差。在M2辅助下的回归克里格法预测结辅助下的回归克里格法预测结果与实测值最为接近。果与实测值最为接近。研究区：北京市延庆县南部北京市延庆县南部采样间隔：300m 300m采样点：210个（训练点186个，验证点24个）插值方法：基于基于对称对数比方法转换对称对数比方法转换的的普通克里格法普通克里格法基于非对称对数比方法转换基于非对称对数比方法转换的普通克里格法的普通克里格法成分克里格成分克里格法法3、土壤质地空间预测方法的比较、土壤质地空间预测方法的比较三、部分研究进展三、部分研究进展结果结果表明表明 成分克里格方法预测

22、精度最高，RMSE相对提高将近17%。成分克里格是直接对土壤质地进行插值，能保证其是最优无偏估计。 非对称对数比转换的非对称对数比转换的普通克里格普通克里格基于成分数据特征和环境变量辅助的土壤质地空间插值基于成分数据特征和环境变量辅助的土壤质地空间插值地理位置、样点和高程图地理位置、样点和高程图土地利用与土壤图土地利用与土壤图三、部分研究进展三、部分研究进展土壤有机质土壤有机质土壤类型土壤类型土壤容重土壤容重高程高程土壤质地土壤质地土土 壤壤 质质 地地 空空 间间 分分 布布已考虑成分数据特征的方法已考虑成分数据特征的方法方差分方差分析析成分克里格成分克里格法（法（CK）基于对称对基于对称对

23、数比转换的数比转换的普通克里格普通克里格法（法（O K -SLR）基于对称基于对称对数比转对数比转换的回归换的回归克里格法克里格法(RK-SLR）将成分数据特征和环境变量考虑到土壤质地空间插值中将成分数据特征和环境变量考虑到土壤质地空间插值中最优辅最优辅助变量助变量组合组合数据转换数据转换三、部分研究进展三、部分研究进展对称对数比转换对称对数比转换（SLR）结论结论u综合考虑综合考虑均方根误差均方根误差（RMSE）、标准化克里格方差标准化克里格方差（MSDR） 、RMSE相对提高值和相对提高值和Aitchisons distance (DA) ，成，成分克里格（分克里格（CK）预测精度最高。）

24、预测精度最高。RIRI相对提高值，以对称对数比转换的普通克里格法为参考方法相对提高值，以对称对数比转换的普通克里格法为参考方法（成分克里格法）（成分克里格法）（基于对称对数比转（基于对称对数比转换的回归克里格法）换的回归克里格法）（土壤颗粒）（土壤颗粒）均方根误差提高值均方根误差提高值三、部分研究进展三、部分研究进展47%13 %46%46%8%6%不同预测方法的不同预测方法的Aitchisons distance图图t检验显示，成分克里格法与基于对称对数比转换的检验显示，成分克里格法与基于对称对数比转换的回归克里格法与基于对称对数比转换的普通克里格法回归克里格法与基于对称对数比转换的普通克里

25、格法DA存在显著差异，表明前两种方法明显优于后者。存在显著差异，表明前两种方法明显优于后者。CK较为集中较为集中三、部分研究进展三、部分研究进展( )iz x（1）顺序指示模拟算法在土壤质地空间分布三维模顺序指示模拟算法在土壤质地空间分布三维模拟中的应用拟中的应用15 km2区域，139个土壤剖面，运用指示变异函数和顺序指示模拟算法分析和模拟土壤质地的空间分布。4、土壤质地的空间分布三维模拟算法研究土壤质地的空间分布三维模拟算法研究三、部分研究进展三、部分研究进展三维顺序指示模拟实现三维顺序指示模拟实现结果表明，土壤质地在垂直方向上具有很强的空间自相关性；而在水平方向上空间自相关性较弱，表明土

26、壤质地有可能在水平方向变异较大，或者本研究所采样的水平采样间隔稍大。顺序指示模拟算法在这种尺度下能够反映土壤质地层次的空间变异。三、部分研究进展三、部分研究进展（2）顺序指示和转移概率指示模拟方法对小尺度土壤顺序指示和转移概率指示模拟方法对小尺度土壤粘粒含量空间模拟的对比研究粘粒含量空间模拟的对比研究以1 m3土体，运用激光粒径分析仪对粘粒含量进行了测定，并将测定数据转化为指示数据。分别运用指示变异函数和转移概率变异函数进行空间变异分析。 三、部分研究进展三、部分研究进展顺序指示模拟（SIS）方法低估了主变量，同时高估了辅变量，而转移概率指示模拟（TPROGS）能将这两种变量很好的再现。此外，

27、与实测数据相比，在SIS结果中粘粒含量分布的空间连贯性显著地减少，而在TPROGS中却能很好的再现。因此，TPROGS更能够反映粘粒含量的三维空间分布特征。顺序指示模拟结果三维实现顺序指示模拟结果三维实现转移概率指示模拟结果三维实现转移概率指示模拟结果三维实现三、部分研究进展三、部分研究进展5、 空间分析的尺度效应空间分析的尺度效应对比分析不同空间尺度对比分析不同空间尺度下矩量法（下矩量法（MoM）和有）和有限最大似然法（限最大似然法（REML）在计算变异函数准确性在计算变异函数准确性方面的表现。方面的表现。三、部分研究进展三、部分研究进展空间尺度对变异函数的影响空间尺度对变异函数的影响研究区

28、位置：北京市大兴区研究区位置：北京市大兴区尺度样区选择：尺度样区选择：40km、30km、20km、10km目标变量：土壤有机质、土壤有效锰目标变量：土壤有机质、土壤有效锰 有限最大似然法在提高空间预测精度方面效果是否明显与有限最大似然法在提高空间预测精度方面效果是否明显与研究区的空间尺度紧紧相关。研究区的空间尺度紧紧相关。 土壤有机质和土壤速效锰的分析结果表明，有限最大似然土壤有机质和土壤速效锰的分析结果表明，有限最大似然法在法在20km和和30km尺度具有比较好的表现尺度具有比较好的表现 。三、部分研究进展三、部分研究进展三、部分研究进展三、部分研究进展6、 土壤属性的时间演替土壤属性的时

29、间演替-6.000.003.006.00OM2000-19800km5km10km15km20km(a) 1980(b) 2000(c) 2000-19800.006.0010.0020.0030.000km5km10km15km20kmOM20000.006.0010.0020.0030.000KM5KM10KM15KM20KMOM1980(a)1980(b)2000(c)2000-1980土壤有机质含量变化土壤有机质含量变化平均增长平均增长46.81%(从从 8.78 到到 12.89g/kg )三、部分研究进展三、部分研究进展大华山镇峪口镇靠山集乡镇罗营乡王辛庄乡黄松峪乡夏各庄镇熊儿寨乡

30、东高村镇南独乐河镇马坊镇韩庄地区山东庄镇平谷镇刘家店乡马昌营镇大兴庄乡海子水库1980-2005年平谷区农用地有机质含量变化图048km单位：g/kg-14 - -8-8 - -4-4 - 00 - 22 - 55 - 77 - 101980-2008土壤有机质含量下降？土壤有机质含量下降？过量施用化肥，综合影响，过量施用化肥，综合影响，综合因素！综合因素！土壤有机质也有质量问题土壤有机质也有质量问题秸秆还田、绿肥作物！秸秆还田、绿肥作物！三、部分研究进展三、部分研究进展（二）土壤（二）土壤-作物系统过程模型与应用作物系统过程模型与应用 土壤水分胁迫修正 氮素转化运移参数修正 土壤氮素胁迫修正 土壤水分模拟子模块 土壤热运动参数修正 氮素转化运移参数修正 作物根系吸水 水氮优化管理措施 土壤温度模拟子模块 作物生长发育模拟子模块 土壤氮素模拟子模块 作物根系吸水 不同的水氮管理措施 图图 5.1 联合模型总体关系图联合模型总体关系图 Fig 5.1 Framework of main model 联合模型总体关系图联合模型总体关系图土壤水土壤水分模拟分模拟子模块子模块土壤氮土壤氮素模拟素模拟子模块子模块综合了胡克林、李保国、陈研等人的部分研究成果综合了胡克林、李保国、陈研等人的部分研究成果气象模块