耕地质量下降、水土流失、土壤酸化、土壤退化、生物多样性减少等问题都对当代中国的粮食安全造成影响^[1]。土壤是农业生产的基础，土壤的肥沃程度和灌溉条件直接决定了农作物的生长发育状况。气候、生物、母质、地形、时间、内动力地质作用以及人类活动等因素都对土壤的发生产生影响，而降水因素是影响土壤发生的最主要气候因素之一^[2]。河流、湖泊等地表水体类型是水资源的重要组成部分，地表水体时空分布与大气降水之间存在着密切联系^[3-4]，因此地表水体时空分布与土壤发生间也存在直接联系。与此同时，如何对土壤和水这类自然资源类型的空间分布格局进行有效的定量化表达对粮食安全问题的相关研究也具有一定意义。土壤多样性(pedodiversity)是21世纪以来世界土壤地理学研究领域的前沿性内容^[5-8]，它可以定量化地描述土壤的空间变异性和不同土壤类型的分布格局问题，为土壤发生、土壤资源保护和可持续合理利用提供理论及数据支持。近年来，基于土壤多样性的相关研究结合了地形地貌^[9]、数学结构^[10-11]、土壤资源保护^[12]、土地利用变化^[13-15]、自然遗产^[16]等不同研究要素，使土壤多样性这一研究专题在国内外均获得了长足发展。

土壤多样性理论与研究方法在21世纪初引入中国^[17-18]，涉及该领域的报道近年来出现了系列性研究^[19-23]。在本项目的前期研究中，为深入探索水要素在土壤形成中所起的作用，使用了一种改良的土壤多样性计量方法并尝试将其引入水、土资源间的关联性分析中，这些研究发现地表水体的空间分布与区域内的土壤种类构成存在一定联系^[24]。但这些研究结论是在市县级别(面积低于5 000 km²)^[24]的研究区数据中获得，研究尺度较小，相关研究结论是否具有普适性仍存疑问，同时前期研究内容仅分析了地表水体对区域土壤构成的影响作用，并未分析土壤分布对水体空间分布的反向影响。以上述研究内容为理论和方法基础，以中国中部典型农业区为例，本研究探索了地表水体分布对区域土壤构成的影响，然后从地理空间分布格局角度分析了具体土壤类型(以研究区内代表性土类为例)对区域地表水体分布的影响，以期对水、土两种重要自然要素在区域空间分布格局上的交互关系研究以及水要素在土壤形成中的作用评价进行进一步探索，并为计量土壤学的发展、区域土地资源的合理配置和农业发展提供理论及数据支持。

1 材料与方法 1.1 研究区概况

选取河南省中部、东部和北部的典型区域为研究区，包括平顶山市(部分区域)、漯河市(部分区域)、许昌市、郑州市、安阳市、鹤壁市、濮阳市、新乡市、开封市、商丘市和周口市(部分区域)，区域总面积66 574 km²。区内地形多以平原为主，高程范围–52~1 662 m，其中86% 的区域高程值低于200 m，为区域内的农业生产提供了有利条件。研究区属暖温带季风气候，四季分明，农业基础条件较好，以小麦–玉米轮作种植为主，是河南省小麦主产区，该区的粮食产出对河南省乃至全国粮食安全具有重要意义。研究区境内河流分属海河流域、黄河流域及淮河流域，年平均降水量700 mm，土壤类型共包含14个土类、32个亚类、94个土属。

1.2 数据来源

研究选用2000—2014年美国Landsat系列卫星数据(其中2013年及以后数据使用OLI_TIRS传感器数据，其余使用TM或ETM+传感器数据)，所有数据均获取自5月份左右，云量均小于5%。其余研究数据包括基于第二次全国土壤普查的研究区土壤图矢量数据、行政区划矢量数据、DEM数据、统计年鉴数据等。数据运行环境为ENVI 4.8和ArcGIS 10.2。

1.3 研究方法

研究使用以下改良的土壤多样性计量方法^[13]对水、土资源的分布格局特征和关联性进行综合分析：

${Y_h} = \frac{{ - \sum\limits_{i = 1}^S {{P_i}\ln {p_i}} }}{{\ln S}}$

(1)

式中：S和p_i的定义如下：①在土壤构成组分多样性里，S表示土壤分类单元个数，p_i表示第i个土壤类型在所有土壤类型总面积里所占的比例。在此前提下，多样性指数Y_h表示研究区内土壤构成组分的多样性特征，即所有分类单元在数量构成上的均匀程度。②在土壤(地表水体)空间分布多样性里，S表示空间网格的数目，p_i表示第i个空间网格里某个土壤类型(地表水体)的面积在该土壤类型(地表水体)总面积里所占的比例。在此前提下，多样性指数Y_h表示研究区内该土壤分类单元(地表水体)的空间分布多样性特征，它表征了单个土壤类型(地表水体)的空间分布离散性程度，即空间分布的多样性格局问题。两种情况下，多样性指数Y_h取值区间同为[0, 1]，当相对丰度分布极度不均匀，也就是当一个或者少数几个对象占支配地位时，例如研究区内只有一种土壤类型(情况①)或所有地表水体只分布在一个空间网格中(情况②)时，Y_h取值趋于0；当每个对象都均匀分布时，例如研究区内有多种土壤类型且各自面积相等(情况①)或地表水体在每个空间网格中均有分布且各自面积相等(情况②)时，Y_h取值等于1。在情况①下，Y_h、Y_h的分子部分和S分别等同于传统土壤多样性计量方法中的Pielou均匀度指数、Shannon多样性指数和丰富度指数，而情况②涉及到空间网格概念，因此在评价系统内(土壤分类系统或土地利用分类系统等)特定分类单元的空间分布离散性特征(即地理空间分布格局问题)时表现得比传统土壤多样性计量方法更加直观和精确。

本研究首先对所选研究区的遥感数据进行土地利用监督分类，提取其中的地表水体信息(主要由河流、湖泊、水库、鱼塘、高等级沟渠等地表水资源类型组成)，地表水体的分类结果由基于Google Earth高清卫片的目视解译方法进行精度验证。2010年开工的南水北调中线工程、近20年高速城市化发展造成的土地利用变化、气候变化等因素都对研究区的地表水体分布特征造成了显著影响，为获取更加客观的研究区地表水体常态分布，本研究对基础遥感数据进行初期筛选，选择平均云量小于5% 的5月份左右数据，同时每个县市选择2000年左右及2014年左右的两组数据进行土地利用监督分类，取两组地表水体分类结果的公共部分作为该县市的地表水体分布最终数据(图 2A)，而基于面积信息的面状矢量水体比单纯基于长度的线状矢量水体包含更多水要素信息，能够更加准确地评价区域地表水体的空间分布格局。经过图层叠置分析等处理后，发现最终的地表水体几乎都由河流、湖泊及少部分水库、鱼塘构成，而研究区内的水库、鱼塘等非自然水体类型也与其附近的河流、湖泊等自然水体类型密切相关。为更加客观及全面地评价研究区不同区域间土壤与地表水体空间分布格局的交互关系，将研究样区基于地理邻近原则、经度邻近原则、纬度邻近原则分别划分为面积相等的10个分区域，并分别计算3种不同的分区域分割方式下分区域内的地表水体空间分布多样性、土属构成组分多样性和土壤亚类构成组分多样性，并评价3种多样性之间的关系，即地理邻近分割原则下的图 2B，经度邻近分割原则下的图 2C，纬度邻近分割原则下的图 2D。选取研究区内的代表性土类为例，计算分区域土类空间分布多样性，结合分区域地表水体空间分布多样性的计算结果和代表性土类的理化性质，对代表性土类和地表水体之间的地理分布交互关系进行初步探索。

2 结果与分析 2.1 土壤分类单元数目与土壤构成组分多样性的相关性

计算并统计了3种分区域分割方式下，各分区域内1 km网格尺度下的地表水体空间分布多样性、土壤亚类个数及其构成组分多样性、土属个数及其构成组分多样性(表 1)，在土壤图数据分析中删除了居民点、河中沙丘、水系等非土壤类斑块。

对上述统计数据进行一元线性回归分析发现，在地理邻近和经度邻近分区域分割方式下，土壤亚类个数与土壤亚类构成组分多样性、土属个数与土属构成组分多样性之间都存在明显的正相关关系(图 3A、图 3B、图 4A、图 4B)，但在纬度邻近分区域分割方式下，两组数据却并不存在明显联系(图 3C、图 4C)。对研究结果进行F检验，在95% 的置信水平下，图 3A、图 3B、图 4A、图 4B中的F值分别为0.001 391、0.000 074、0.031 923、0.003 967，均小于0.05，回归分析结果有效；图 3C、图 4C中的F值分别为0.883 052、0.749 711，均大于0.05，回归分析结果无效。地理邻近和经度邻近分区域分割方式下的此种数据表现与前期研究中的相关结论^[24]相一致，而纬度邻近分区下的数据表现不一致。在此做出如下结论：随着区域内土壤分类单元数量的增加，区域内的土壤构成一般也将会变得更加均匀，但这一现象在不同纬度地区的横向对比中无明显表现。

2.2 地表水体空间分布多样性与土壤分类单元数目的相关性

研究发现，在地理邻近和经度邻近分区域分割方式下，地表水体空间分布多样性与土壤亚类个数、土属个数之间都存在一定的负相关关系(图 5A、图 5B、图 6A、图 6B)，但在纬度邻近分区域分割方式下，两组数据却并不存在明显联系(图 5C、图 6C)。对研究结果进行F检验，在95% 的置信水平下，图 5A、图 5B、图 6A、图 6B中的F值分别为0.162 385、0.032 202、0.086 470、0.053 044，少部分回归分析结果有效；图 5C、图 6C的F值分别为0.929 129、0.442 140，回归分析结果无效。上述研究发现与前期研究中的相关结论^[24]存在明显区别，前期研究在小尺度(5 000 km²左右)研究区数据支持下发现地表水体的离散性分布特征与区域土壤分类单元数目之间存在更加明显的负相关关系。在此做出如下结论：区域地表水体更加离散的空间分布特征会造成区域内土壤分类单元数目的减少，但这种现象随着研究区面积的增加将会逐渐弱化，也即在大尺度研究区内，除水要素以外的其他成土因素将对区域土壤构成产生主要影响，但在小尺度研究区内，生物、母质、地形、时间、内动力地质作用以及人类活动等成土因素相对固定，因此小尺度研究区内的水要素分布特征将对区域土壤构成产生主要影响。

2.3 代表性土类分布对区域地表水体空间分布格局的影响

在以往常规的成土因素相关分析中，一般均着眼于气候、生物、母质、地形、时间、内动力地质作用以及人类活动等因素对土壤形成的作用分析，而从土壤形成角度探讨土壤分布现状与这些成土因素交互影响的相关报道较少出现。上一章节的研究内容主要分析区域地表水体的空间分布特征对区域内土壤构成混乱程度的整体影响，本节尝试分析单一土壤类型与区域地表水体空间分布的交互影响。

当土壤分类级别由土属变为土类时，不同土壤类别间的理化性质具有更大的区别，为深入探索区域土壤分布与地表水体空间分布格局之间存在的客观联系，基于前期研究的相关方法及理论，对研究区内的土壤分类单元(以土类为例)进行筛选，选择潮土和褐土这两种土类作为研究区的代表性土壤类型进行后续研究，这两种土类同时满足面积最大(分别占研究区土壤总面积的65.08% 和19.93%)和空间分布最广泛(1 km网格尺度下的土壤空间分布多样性均大于0.8)的代表性土壤类型选取原则(表 2、图 7)。

潮土是一种受地下潜水影响和作用的半水成土壤，是由潴育化过程和旱耕熟化影响的腐殖质积累过程发育而形成的。潮土分布地区地形平坦，地下水埋深较浅，水热资源较为丰富，一般多具有比较悠久的地区农业历史。潮土是河南省面积最大、分布最广的一个土类^[25]。计算并统计了3种分区域分割方式下，各分区域内1 km网格尺度下的潮土、褐土空间分布多样性(表 3)，结合各分区域分割方式下的地表水体空间分布多样性数据(表 1)，研究发现，随着区域内潮土空间分布多样性的升高，地表水体的空间分布离散性也有升高的趋势，这一现象在地理邻近和经度邻近分区域分割方式下表现明显(图 8A、图 8B)，但在纬度邻近分区方式下无表现(图 8C)。对上述研究结果进行F检验，在95% 的置信水平下，图 8A、图 8B、图 8C中的F值分别为0.017 086、0.000 727、0.486 264，地理邻近和经度邻近分区域分割方式下的回归分析结果有效，纬度邻近分区方式下无效。研究认为形成该数据现象的原因有两个：一是区域地表水体分布对土壤形成的影响，即潮土的成土条件之一是区域性的潮湿水分状况，因此分区域内高分布离散性的地表水体为潮土的形成提供了有利条件；二是发挥潮土生产潜力的前提是发展灌溉，即研究区内潮土区广阔的适种性要求加强农田水利设施建设，而潮土区多为平原地形也为农田水利设施的建设提供了客观上的有利条件，因此潮土分布离散性的增强往往给区域内地表水体分布离散性的增强提供了反作用力。

褐土是河南省仅次于潮土的第二大土类，是在暖温带半湿润季风气候、干旱森林和灌木草原植被下，经过黏化过程和钙积过程发育而成的半淋溶性土壤，褐土区地下水位多在10 m以下，有的达数十米，因此水要素一般不参加褐土的成土过程^[25]。由于褐土在研究区内的分布特征所限，在某些分区内无褐土分布(分区内褐土面积为0)或仅有极少量分布(分区内褐土面积占研究区褐土总面积的比例低于0.1%)，为客观分析褐土的空间分布对地表水体空间分布造成的影响，在回归分析中删除了地理邻近分区域分割方式下的4、7、8、10和经度邻近分区域分割方式下的7、8、9、10号分区数据点。研究发现，在地理邻近和经度邻近分区域分割方式下，随着区域内褐土空间分布多样性的升高，地表水体空间分布离散性有降低趋势(图 9A、图 9B)，但在纬度邻近分区方式下该数据现象无表现(图 9C)。对上述研究结果进行F检验，在95% 的置信水平下，图 9A、图 9B、图 9C中的F值分别为0.074 677、0.072 122、0.446 920，3种分区方式下的回归分析结果均无效，但地理邻近和经度邻近分区域分割方式下的F值很接近0.05，这主要是由于删除无效数据点以后造成数据点总数过少引起的，但图 9A、图 9B的回归方程仍具有一定分析价值。后续研究中，将单独对褐土的主要分布区域进行分区域分析，应该能得到更加准确的拟合方程。本研究认为，由于水要素一般不参与褐土的成土过程，因此褐土与地表水体之间在空间分布格局上存在一定的负相关变化关系。

本节研究结果表明，改良的土壤多样性计量方法可以应用于某一具体土壤类型的空间分布格局与水要素空间分布格局的交互关系研究，而且这种定量化的分析方法不仅适用于潮土这类与水要素存在密切联系的土壤类型，同时还适用于褐土这类与水要素没有太大联系的土壤类型。

3 结论

1) 河南省最具代表性的土类——潮土的空间分布格局与区域地表水体的空间分布格局存在密切联系，两者的空间分布离散性之间存在较为明显的正相关关系，一元线性回归分析中的判定系数R²达到0.5以上；而另一代表性土类——褐土的空间分布格局与区域地表水体的空间分布格局之间却表现出较为明显的负相关变化关系，判定系数R²也达到0.5以上。上述数据现象与水要素在两种土壤类型成土过程中所起的重要程度密切相关。

2) 在分区域研究中，地理邻近和经度邻近分区域分割方式下的数据结果表现相近，但其均与纬度邻近分区的数据结果不一致。鉴于地理邻近分区理念与“相邻相近”的地理学第一定律更加贴合，在后续研究中将主要使用地理邻近分区方式探索不同区域水、土资源间的交互关系。

3) 土壤和地表水体空间分布格局之间的交互关系在大尺度和小尺度研究区中存在不同表现，从水、土资源的地理空间分布角度来看，水要素作为土壤形成过程中的必要环节，在小尺度区域中，由于生物、母质等其他成土因素发生突变的概率较小，因此水要素对土壤形成的影响占比较大；但在大尺度区域中，相比于水要素，其他成土因素一般存在较大的变化，因此其他成土因素对土壤形成的影响更加重要。

4) 基于资源空间分布离散性分析的土壤多样性计量方法及理论除了能够定量化描述区域地表水体空间分布对区域内土壤丰富度指数的影响，还能有效地应用于某个具体土壤类型与地表水要素空间分布格局的交互关系研究，这在本研究中得到有效证明。在未来相关研究中，将首先考虑研究区尺度，即研究区面积大小对研究结论的影响，另外考虑尝试分析不同地形、母质等条件下水、土资源空间分布格局交互性的不同体现。同时，本研究的相关研究理念及方法也为其他资源类型的空间分布格局评价及交互关系研究提供了新的思路。