第1章 绪论
　　1.1 本书创作背景
　　1. 新时期自然资源数量-质量-生态“三位一体”的管理，对耕地和土壤资源的监测评价提出了新要求
　　自然资源是人类赖以生存和发展的基础，也是地球和生态系统的重要组成部分。随着经济、社会的快速发展，自然资源的数量、质量、结构及空间布局等要素深受人类活动影响，工业化、城镇化造成的自然资源数量减少、质量下降、生态破坏等问题逐渐成为目前昀突出的矛盾之一。中国是人口大国，粮食安全是保障社会稳定的重要基础，然而经济建设过程中非农建设造成了自然资源尤其是耕地资源的大规模减少。 1990年以来，中国耕地面积累计减少超过 1000万 hm2以上，每年减少约 69万 hm2。同时，无序开发建设、自然资源的过度开发利用等导致了土地荒漠化及沙化等一系列生态问题，自然生态保护已迫在眉睫，如何合理利用自然资源、促进生态环境可持续发展已引起全球学者的普遍关注。资源利用与整个社会的宏观经济系统相联系，随着区域经济发展，人们的生态环境保护意识随之增强，生态建设技术不断创新，以牺牲资源环境来实现经济增长的方式将导致区域自然资源枯竭、生态系统崩溃。为此，单纯的资源数量管理已无法满足当前经济社会的发展需求，迫切要求自然资源管理向数量-质量-生态“三位一体”管理理念转变，即以自然资源数量管理为基础，重视自然资源的质量提升，显化自然资源的生态功能。
　　目前，在自然资源数量管控方面，我国已经建立了比较完整的技术支撑体系，通过多轮次的土地调查、自然资源调查以及年度更新调查等，已全面掌握了耕地、林地、草地、水面等自然资源的数量及空间分布，建立耕地、林地、草地、水面等自然资源占补实时台账，基本实现了逐图斑及时和准确掌握自然资源数量及其变化情况，达到了自然资源数量精细化管理的要求。自然资源的质量管理方面，经过国土、农业等部门多轮次以县级为基本行政单位的样点调查、土壤采样、测试和分析评价，耕地资源质量等级的本底数据已基本掌握。自然资源的生态功能管护方面，目前尚未形成成熟、系统的调查、评价及管理体系，是当前自然资源管理短板，需进一步加强。
　　为此，应积极探索耕地和土壤资源监测评价的新方法、新途径，尤其是要加强耕地和土壤生态环境的监测评价，以适应新时期自然资源数量-质量-生态“三位一体”管理的新要求。
　　2. 当前耕地资源质量评价的效率、精度和精细化水平，与自然资源管理的要求仍有很大差距，评价方法仍需不断改进
　　为加强耕地资源的质量管理，在国内外相关研究的基础上，我国对耕地资源质量评价的指标体系、评价方法、评价程序及评价成果确定等相继制定了相应的标准规范，也建立了比较完整的耕地资源质量等级本底数据，但评价的效率、准确性以及精细化水平等，与当前自然资源管理的要求还有很大差距，仍需不断地研究探索。其中以下 4个方面在耕地资源质量评价方法改进研究中应受到特别重视。
　　(1)土壤属性指标是耕地资源质量评价指标体系的关键内容。指标体系的构建是耕地资源质量评价的关键工作之一，目前尚未形成统一的评价指标体系，但土壤一直是耕地资源质量的关键和核心要素(傅伯杰， 1990；周生路等， 2004；陈印军等， 2011)。随着指标信息获取手段的多样化，自然、社会、经济、生态、环境等因素不断进入监测指标体系，耕地资源质量评价指标越来越全面，体系结构也越来越庞大，可以多方面、全方位、系统性地反映耕地资源质量(张凤荣等， 2001；孔祥斌等， 2007)。目前，自然资源和农业部门从不同的应用目的出发，建立了不同的耕地资源质量评价指标体系。其中，自然资源部门的《农用地质量分等规程》(GB/T 28407—2012)确定的全国耕地资源质量评价指标共 3大类 12个，立地条件(3个)：地形坡度、地表岩石露头度、障碍层距地表深度；土壤条件(6个)：土壤有机质含量、土壤酸碱度、有效土层厚度、表层土壤质地、剖面构型、盐渍化程度；利用条件(3个)：排水条件、灌溉保证率、灌溉水源。农业部门的《耕地地力调查与质量评价技术规程》(NY/T 1634—2008)确定的全国耕地地力评价指标共 6大类 66个：气象(7个)、立地条件(18个)、剖面性状(5个)、土壤理化性状 (21个)、障碍因素(7个)、土壤管理(8个)。耕地地力质量评价侧重于直接指导农民种植施肥，指导农业生产，从而提高耕地粮食产量；农用地质量分等侧重于国家对土地的宏观管理，提高资源生产能力，从而保护耕地。虽然两套评价体系出发点不同，但无论是国土资源部门还是农业部门的耕地资源质量评价体系，土壤属性指标都是评价指标体系的关键组成部分，关系到耕地资源质量评价结果的准确性及应用价值，因此土壤数据的应用在耕地资源质量评价中具有重要的地位。
　　(2)土壤数据的合理使用是耕地资源质量评价精度的重要保障。目前我国耕地资源质量评价可整合利用的抽样调查数据主要有三类：一是按土壤类型布点采样调查的土壤数据，如耕地地力调查数据、测土配方调查数据等，这是当前耕地资源质量评价土壤采样调查数据的主要来源；二是按网格法布点调查的土壤数据，如多目标区域地球化学调查的土壤数据；三是其他如按乡、村行政区布点采样调查耕地资源质量动态监测的土壤数据等。对耕地资源质量评价而言，三种抽样调查各有优劣。按土壤类型布点充分考虑了同一土壤类型在空间上相对均质的这一特征，但土种单元空间分布较为复杂且单元面积大小不均，一些面积较大的土种单元布点往往较多，而面积较小的土种单元布点较少甚至没有布点，过量布点会造成资源浪费，抽样点过少则造成代表性不足。按网格布点调查点位空间分布均匀，但没有考虑区域的空间异质性，造成在空间差异小的地区样点过密，而在空间差异较大的地区则样点不足，降低采样数据的代表性。按行政区布点抽样便于管理与数据调查，但样点分布过少且布点随意性过大。因此，在整合利用各类土壤调查数据进行耕地资源质量评价时，对布点过多区域如何进行样点及数据的优化利用，对布点过少区域如何进行样点及数据的预测补充，对优化整合后的数据如何进行精度检验评估等均有待进一步研究。
　　(3)优化的空间抽样是提高评价效率与精度的重要手段。目前耕地资源质量评价中，对土壤样点的使用还没有考虑优化选择的问题，无论精度如何均是全部使用，还未考虑通过对样点数量及布局进行优化以提高评价精度。空间抽样方法的运用为这一问题的解决提供了有效的手段。空间抽样方法是探索如何实现在一定的成本约束下(既定的样本数量)，通过合理布设样点和科学推断估计，获得尽可能高精度的抽样结果；或者，如何通过合理布设样点和科学推断估计，利用尽可能少的样点数量，来得到满足精度要求的抽样结果。目前建立了遗传算法、模拟退火等相应的计算机演化算法(Haining，2003；Wang et al.，2013；Diaz-Avalos and Mateu，2014；Goodchild，2013)，并在国内外耕地资源质量或土壤质量调查采样点空间布局研究领域受到广泛重视。但是，现有空间抽样理论的应用研究主要针对单一指标或属性的空间抽样推断，而耕地资源质量涉及土壤的众多属性及土地利用等各种因素，样点的布设需要考虑多指标的空间关联性，如何在耕地资源质量评价中通过空间抽样方法对样点数量及布局进行优化的提高评价效率与精度，还有待进一步研究探讨。
　　(4)数据处理的标准化、精细化是耕地资源质量评价的发展趋势。随着土地管理、农业集约化水平的发展，耕地资源质量成为土地评价的重点内容之一。经过多年的深入探索，我国耕地资源质量评价已经全面开展，形成了丰富的成果，建立了较为完善的耕地资源质量本底数据库。随着信息技术的发展， 3S技术、计算机技术等被广泛应用于耕地资源质量评价中，耕地资源质量评价的信息化程度逐渐提高，耕地资源质量评价体系、方法与程序被进一步完善与改进，向标准化、精确化、智能化方向发展(郧文聚和黄元仿， 2014)。2000年，全国第一轮农用地分等成果比例尺为 1∶10万或更小，主要以计算机辅助的人工处理为主；而到 2011年第二轮农用地分等成果比例尺已达到 1∶5000，开发了各类软件系统，主要以人工辅助的计算机处理为主。目前，对耕地资源质量评价方法的研究主要集中在耕地资源质量的内涵界定、评价体系的建立以及耕地资源质量等级划分方法的研究等方面，自然资源部门的《农用地质量分等规程》(GB/T 28407—2012)以及农业部门的《耕地地力调查与质量评价技术规程》(NY/T 1634—2008)对耕地资源质量评价中的基础参数、评价指标、等级划分进行了标准化规定，是目前为止昀为权威的耕地资源质量评价标准，但两套体系对于耕地资源质量评价中参评因素的数据处理仍没有提出标准的、优化的技术方法。
　　3. 大数据分析从数据本身出发，通过海量数据的深入挖掘和非线性关系的解析，为耕地资源质量评价提供了新的方法手段
　　针对当前耕地资源质量评价存在的效率低、耗时长、精细化程度不够的问题，国内外学者一直在对耕地资源质量调查评价普遍采用的抽样调查估计推测全体的方法进行研究改进(Cressie and Wikle， 2015)，并通过整合利用多源土壤数据(孙亚彬等， 2013)、空间插值优化(Wei et al.， 2015)以及运用遗传算法(Haining， 2003；Goodchild，2013)、模拟退火算法(张淑杰等， 2013)进行样点空间优化等方法来提高评价精度。
　　但就作者在江苏省东海县的前期相关研究来看，全县耕地总面积 122 482hm2，1∶5000土地利用图上耕地图斑 56 965个，耕地质量等级评价土壤属性的各类样点共 1442个，约占耕地图斑的 2.5%，每个样点覆盖约 84.9hm2。显然，以如此密度的样点来估计推断评价全县耕地质量，很难保证结果的准确度，更谈不上精细化的要求。抽样调查估计推测全体的方法属于传统的样本分析方法，旨在寻求以昀少的数据获得昀多的信息，并特别强调数据的精确性(Shin and Choi， 2015)。而区域耕地质量评价所整合利用的多源土壤数据，由于调查时间、处理测试方法等的差异，数据存在可比性问题并影响评价的精度。另外，耕地质量等级评价省级成果按照分区评价逐级汇总方法在县市成果基础上逐级汇总而成，这势必造成评价结果精度的进一步减损。
　　近年发展起来的大数据分析方法，强调使用所有数据分析，允许数据优劣掺杂、不全面和不精确(Roth and Almeida， 2015；Wei et al.， 2015)。大数据分析通过寻找事物之间的相关关系而非因果关系，比以往方法更容易、更快捷、更清楚地分析事物，通过找到现象的良好关联物，捕捉现在和预测未来(李德仁等， 2014；Suthaharan， 2014)。大数据分析无须在还没有收集数据之前就把分析建立在早已设立的少量假设的基础之上，而是直接由数据说话，用数据驱动的大数据相关关系分析法，取代基于假设的易出错的方法，不易受偏见的影响(李国杰和程学旗 ， 2012)。大数据分析能看到一些以往分析无法发现的细节，更清楚地看到样本无法揭示的细节信息，产生额外价值，获得新的发现、新的观点(宋杰等， 2014；Li et al.，2016)。目前，大数据技术的应用主要在商业领域(Rabl et al.， 2013；隋殿志等，2014)，工程和科学领域也逐渐开展大数据的探索研究(Cheng et al.， 2014)。
　　近年来，我国相关部门及科研工作者开展了大量的土地资源调查评价。国土资源部 1999年在全国部署开展了第一轮农用地质量调查评价， 2002年在全国部署开展了多目标区域地球化学调查， 2009年完成了第二次全国土地调查， 2011年在全国部署开展了第二轮耕地质量调查评价， 2011年在全国试点开展了耕地质量等级监测。农业部 2006年在全国试点开展耕地地力评价，地方农业部门开展了长期的测土配方土壤调查，等等。这些工作按照不同的目的要求，进行了大量的实地调查、采样

目录
序
第1章 绪论 1
1.1 本书创作背景 1
1.2 本书主要内容 8
参考文献 11
第2章 土壤样点、空间预测和赋值方法综合优化下的县域耕地资源质量等级评价 14
2.1 耕地质量评价方法完善研究概况 14
2.1.1 评价内容与指标体系的完善 14
2.1.2 耕地土壤属性预测方法的完善 16
2.1.3 评价样点空间抽样方法的完善 18
2.2 研究思路与基础数据 19
2.2.1 研究内容 19
2.2.2 研究技术路线 20
2.2.3 研究区概况 21
2.2.4 研究基础数据 22
2.3 样点疏密差异对土壤属性空间预测精度的影响及分区 23
2.3.1 基于不同来源样点的土壤属性预测总体精度 23
2.3.2 基于土壤属性均质单元的样点密度空间差异 34
2.3.3 土壤属性预测精度与样点密度的关系 40
2.3.4 研究区土壤样点的疏密性分区 43
2.4 样点稀缺区土壤属性空间分布的优化预测 45
2.4.1 优化预测的土壤-景观模型建立 45
2.4.2 研究区样点密集区土壤属性特征与景观因子的关系 49
2.4.3 基于土壤-景观模型的土壤属性空间分布预测 51
2.4.4 土壤-景观预测与地统计插值预测结果比较 56
2.4.5 样点稀缺区土壤属性提取及其分布特征 60
2.5 样点密集区参评土壤样点的数量及布局优化 61
2.5.1 样点数量及布局优化的空间抽样原理与智能优化算法 61
2.5.2 基于模拟退火算法的参评土壤样点数量及布局优化 63
2.5.3 考虑土壤属性空间异质性的模拟退火算法改进及运用 70
2.5.4 优选土壤样点的精度评价与来源分析 77
2.6 评价单元土壤属性数据赋值方法的优化选择 80
2.6.1 评价单元土壤属性数据赋值过程 80
2.6.2 不同土壤属性*优空间插值方法选择 81
2.6.3 评价单元土壤属性赋值方法优化选择 87
2.7 各优化方法综合运用的县域耕地质量评价 92
2.7.1 评价的方法体系 92
2.7.2 评价单元数据获取 93
2.7.3 耕地质量评价结果 95
2.7.4 评价精度分析 96
2.8 本章小结 98
参考文献 100
第3章 运用大数据分析方法的省域耕地资源质量等级评价 103
3.1 大数据研究概况 103
3.1.1 大数据研究的兴起 103
3.1.2 大数据信息来源 104
3.1.3 大数据管理与云计算 105
3.1.4 大数据分析与机器学习研究 106
3.2 省域耕地资源质量等级评价内容及技术路线 107
3.2.1 省域耕地资源质量等级评价内容 107
3.2.2 省域耕地资源质量等级评价技术路线 108
3.2.3 江苏省概况 110
3.3 耕地资源质量基础大数据集和云计算环境构建 110
3.3.1 耕地资源质量的驱动因素与被驱动因素 111
3.3.2 研究基础大数据集的构建 112
3.3.3 大数据云计算环境的构建与数据集成 117
3.4 大数据深度学习的省域耕地资源等级评价 121
3.4.1 基本流程与数据处理 122
3.4.2 大数据耕地资源质量评价研究结果与精度 127
3.4.3 大数据耕地资源质量评价空间差异与影响因素 132
3.4.4 大数据深度学习方法在耕地资源质量评价中的应用效果评价 134
3.5 方法融合下的大数据耕地资源质量评价结果的精度效应 136
3.5.1 面向对象方法 136
3.5.2 深度学习与面向对象方法融合的耕地资源质量评价结果 140
3.5.3 深度学习与面向对象方法融合的耕地资源质量评价精度 140
3.6 本章小结 143
参考文献 145
第4章 基于农田大气-土壤-作物多介质系统监测的耕地污染和风险评价 149
4.1 各介质重金属污染特征与风险概况 149
4.1.1 大气重金属污染特征与风险概况 149
4.1.2 土壤重金属污染特征与风险概况 150
4.1.3 作物重金属污染特征与风险概况 151
4.2 多介质系统监测的耕地污染和风险评价 152
4.2.1 耕地污染和风险评价对象与技术路线 152
4.2.2 蠡河流域概况 153
4.2.3 监测点布设与样品采集 156
4.2.4 样品的测试 160
4.2.5 数据分析和评价方法 163
4.3 大气降尘重金属污染特征和风险评价 169
4.3.1 大气降尘重金属时空分布特征 169
4.3.2 大气降尘重金属生态风险评价 175
4.3.3 大气降尘重金属健康风险评价 178
4.4 土壤重金属污染特征和风险评价 180
4.4.1 土壤基本理化性质 180
4.4.2 土壤重金属分布特征 182
4.4.3 土壤重金属生态风险评价 190
4.4.4 土壤重金属健康风险评价 192
4.5 作物重金属污染特征与风险评价 195
4.5.1 小麦重金属污染特征与风险评价 195
4.5.2 水稻重金属污染特征和风险评价 202
4.5.3 两种作物籽粒重金属污染特征和风险的比较 207
4.6 农田系统多介质重金属综合健康风险 211
4.6.1 饮用水重金属的健康风险 211
4.6.2 不同输入途径下重金属污染产生的健康风险比较 213
4.6.3 综合健康风险评价及空间分布 214
4.7 本章小结 215
参考文献 217
第5章 多方法联用的农田大气-土壤-作物系统重金属污染源解析 221
5.1 农田系统重金属污染源解析研究概况 221
5.1.1 大气降尘重金属源解析 221
5.1.2 土壤重金属污染源解析 221
5.1.3 作物重金属污染源解析 222
5.2 农田各介质重金属污染源解析 222
5.2.1 拟解决的关键问题 222
5.2.2 研究内容与技术路线 223
5.2.3 研究材料 223
5.2.4 重金属污染源解析方法 224
5.3 大气降尘重金属污染源解析 228
5.3.1 相关分析与富集因子分析 229
5.3.2 PMF分析和Pb同位素分析 230
5.3.3 大气降尘重金属源解析方法联用探讨 234
5.4 土壤重金属污染源解析 235
5.4.1 定性源识别 235
5.4.2 定量源解析 237
5.4.3 土壤源解析方法联用探讨 249
5.5 作物重金属污染源解析 251
5.5.1 小麦籽粒元素Pb来源解析 252
5.5.2 水稻籽粒元素Pb来源解析 254
5.5.3 两种作物籽粒Pb污染源解析结果比较 256
5.6 本章小结 257
参考文献 258
第6章 基于监测-模拟的城市地区土壤-大气多环芳烃分布变化、来源和风险预警 262
6.1 多环芳烃的来源、性质和分布 262
6.1.1 多环芳烃的来源 262
6.1.2 多环芳烃的性质 262
6.1.3 土壤多环芳烃的空间分布 263
6.2 多环芳烃分布、来源与预警的研究内容与方法 264
6.2.1 研究内容与技术路线 264
6.2.2 研究区概况 266
6.2.3 样品采集与测定 266
6.3 城市土壤-大气多环芳烃的空间分布与来源 272
6.3.1 城市土壤多环芳烃的污染水平与空间差异 272
6.3.2 城市大气PAHs时空分布特征 278
6.3.3 城市土壤-大气多环芳烃的来源解析 282
6.4 城市环境多环芳烃的多介质空间分布模拟 286
6.4.1 多环芳烃排放量估算 287
6.4.2 多介质逸度Ⅲ级模型及其参数 288
6.4.3 OC-Model与BC-Model模拟结果对比与选择 295
6.4.4 BC-Model模拟多环芳烃的空间分布 299
6.5 城市土壤多环芳烃累积变化的模拟重建与预测 303
6.5.1 多介质逸度Ⅳ级模型及其参数 304
6.5.2 城市土壤多环芳烃累积变化模拟重建与预测的方法 307
6.5.3 1978～2013年城市土壤多环芳烃的累积变化与未来预测 310
6.5.4 城市土壤多环芳烃累积变化的时空特征 313
6.6 城市土壤多环芳烃累积的风险评估与预警 314
6.6.1 健康风险评估和预警的方法 315
6.6.2 城市土壤多环芳烃累积的健康风险评估 318
6.6.3 城市土壤多环芳烃累积的健康风险预警 321
6.7 本章小结 322
参考文献 323
第7章 基于溪流溶解态有机质监测的土地利用、土壤环境变化的效应评价 327
7.1 溪流溶解态有机质性状及土地利用、土壤环境变化的影响 327
7.1.1 溪流有机质的生物地球化学功能 327
7.1.2 溪流溶解态有机质的来源 327
7.1.3 溪流溶解态有机质的组成 328
7.1.4 土地利用、土壤环境变化对溪流溶解态有机质性状的影响 329
7.2 溪流溶解态有机质检测效应评价研究内容与方法 331
7.2.1 研究内容与技术路线 331
7.2.2 研究材料与方法 333
7.2.3 样品测试分析 335
7.2.4 数据分析 339
7.3 流域溪流有机质陆源输入特征对土地利用、土壤环境变化的响应 339
7.3.1 陆源有机质输入估算模型的耦合 339
7.3.2 流域水体有机碳、氮陆源输入特征的空间变异 347
7.3.3 流域水体有机质陆源输入特征对土地利用、土壤环境变化的响应 352
7.3.4 典型溪流有机质陆源输入特征对土地利用、土壤环境变化的响应 356
7.4 典型溪流溶解态有机质性状对土地利用、土壤环境变化的响应 359
7.4.1 不同土地利用-土壤类型影响溪流中溶解态有机质输出含量的特征差异 359
7.4.2 不同土地利用-土壤类型影响溪流中溶解态有机质组成的特征差异 367
7.4.3 溪流溶解态有机质组成对土地利用、土壤环境变化的响应 376
7.5 本章小结 381
参考文献 382
第8章 基于污染物湖泊沉积记录的资源环境保护的效应评价 387
8.1 湖泊沉积与污染物的变化记录 387
8.1.1 湖泊沉积物定年 387
8.1.2 湖泊沉积物中多环芳烃的累积变化 388
8.1.3 湖泊沉积物中重金属的累积变化 389
8.2 湖泊沉积物研究内容与方法 390
8.2.1 研究内容与技术路线 390
8.2.2 研究区概况 391
8.2.3 样品的采集与测定 392
8.2.4 数据分析和评价模拟的方法 396
8.3 西太湖沉积物中多环芳烃的历史变化与物源分析 400
8.3.1 西太湖沉积环境与沉积速率测算 400
8.3.2 西太湖多环芳烃的历史变化 402
8.3.3 多环芳烃吸附沉积的影响因素 411
8.3.4 西太湖多环芳烃的