近年来，以流域为单元的生态学研究已经成为人与自然相互关系研究领域的热点。结构与功能完善的流域生态系统能有效释放生产潜力并放大系统的生态经济效益，对改善地方生态环境、促进当地经济建设至关重要。但在自然和社会因素的双重影响下，流域多呈现出不同程度的景观破碎化现象^[1-2]，进而导致其生态功能减弱并逐步退化。景观破碎化是描述景观格局的重要参数，直接影响着景观的结构、功能及其内部的生态功能的变化^[3-4]，对其进行研究对深入理解景观格局的形成与变化机制具有重要意义。

大尺度上景观破碎化的研究主要是基于土地利用展开，而土地利用类型的组合可以呈出现不同的景观空间格局，因此，在进行景观破碎化研究时更需强调空间格局与生态过程的关系^[5]。当前国内外对景观破碎化的研究大多还是使用多个景观指数来综合描述景观破碎化的过程^[6-8]，而对定量测度景观破碎化的程度及对景观破碎化的空间变异特征的研究相对较少^[9-11]，定量测度和空间变异结合的研究也鲜有报道。

宁夏黄河流域位于我国西北内陆，属典型的脆弱内陆流域生态系统。由于流域气候整体呈现暖干化变化，使地表水资源减少，而人口和经济增长导致了土地利用的快速转变和地下水资源持续减少，流域原有景观格局发生改变，呈现出不同程度的景观破碎化，使原本脆弱的生态环境进一步恶化。本文选取有效粒度尺寸m_eff)^[12]来定量分析宁夏黄河流域的景观破碎化水平，并以地统计学的方法确定流域景观破碎化的特征尺度，对流域景观破碎化空间异质性进行分析，以此揭示流域生态状况及空间变异特征，为流域土地利用管理、景观稳定及可持续发展提供科学依据。

1 研究区域概况

黄河自宁夏北部穿境而过，境内全长397 km，有清水河、苦水河等支流注入，形成了以黄河主干及其支流组成的北部宁夏黄河流域。流域范围介于35°50′—39°23′N，104°17′—107°12′E之间，面积约为4.16万km²，占宁夏总面积的81%。流域内地形复杂多样，南部为黄土丘陵区、中部为山间平原、北部为宁夏平原，三者形成了流域内面积最大的三级地势阶梯，其余地区山地、丘陵、平原和沙漠等交错分布。气候属于典型的大陆性气候，年平均气温11.5—14℃，日照时数2200—3000h，无霜期 110—140d，全年平均降水量为200—800mm。根据河流及用水调度的差异，将宁夏黄河流域分为引黄灌区、黄右区间、黄左区间、清水河、苦水河、红柳沟和盐池内流区7个分区(图 1)。

2 研究方法 2.1 数据来源及处理

文获取了1985年和2010年两期Landsat-5 TM影像数据(30m×30m，来源于“国际科学数据服务平台”)、 DEM(90m×90m，来源于“国际科学数据服务平台”)等数据。使用ERDAS IMAGINE 9.2对TM影像进行辐射校正、几何纠正和配准，采用监督分类提取耕地、林地、草地、水域、建设用地和未利用地6类主要景观类型。参考与研究区自然条件相似的已知区域建立对应判读标志，并结合人-机交互式对错分的地区进行修改，分类结果借助Google Earth 软件^[13]、相关航拍数据等进行验证^[14-15]，结果表明1985年和2010年两期影像分类的Kappa系数分别达到为0.80和0.78，解译精度满足研究要求(图 2)。

2.2 景观破碎化指数选取

选取有效粒度尺寸(m_eff)来定量分析流域景观破碎化水平，利用软件Fragstats 4.1及Excel进行景观破碎化指数的计算和统计分析。有效粒度尺寸融合了生态过程、景观组分与空间格局，从而可以更为综合、客观地表征景观的破碎化状况。有效粒度尺寸越小，破碎化程度就越高^[4]。计算公式如下^[12]：

(1)

式中，m_eff(j)表示为景观j的有效粒度尺寸，n为景观 j中非破碎斑块的数量，A_ij表示景观j中斑块i的面积大小，A_j为景观j的面积大小。该指数的取值范围是：最小值为栅格大小，此时相邻斑块之间的类型均不相同；最大值为景观面积，此时该景观具有唯一的类型。

2.3 景观破碎化空间异质性分析

空间异质性是产生空间景观格局的主要原因^[16-17]，空间异质性分析必须考虑景观格局的尺度效应^[18]。本研究在连续空间尺度上对景观破碎化进行空间计算，借助半变异函数，以地统计学方法^[19-20]揭示变量的空间异质性，以此来确定合适的分析尺度(特征尺度)。

(1) 半变异函数 半变异函数是地理现象分布中的空间依赖性与空间异质性的一个综合性衡量指标，主要参数包括：块金值 C₀、基台值 C₀+C 及变程 A₀。其中块金值与基台值之比 C₀ /(C₀+C) 则反映了随机部分引起的空间异质性占总空间异质性的比重。其定义为：

(2)

式中，h为两样本点的空间分隔距离，Z(x_i)与 Z(x_i+h)分别是区域化变量Z(x)在空间位置x_i和 x_i +h上的观测值(i = 1,2,…,N(h))，N(h)是分隔距离为h时的样本对总数。

(2) 特征尺度确定 利用Fragstats 4.1的“滑窗”功能，获取不同尺度下景观破碎化的空间变异特征，再运用GME(Geospatial Modelling Environment)生成20000个随机点，并在Arcgis10.2中对样点进行赋值，最后利用GS+7.0软件对连续尺度序列上景观破碎化指数的空间变异特征值进行计算，进而确定流域景观空间异质性的特征尺度。

3 结果与分析 3.1 流域景观破碎化的数量变化

宁夏黄河流域7个分区的面积和尺度不同，各分区的m_eff直接作为破碎化程度的比较依据将带来较大误差。因此，本文以流域分区的m_eff与其分区面积的比值来表征该流域分区的破碎化程度(百分比)，比值越大破碎化程度越低，比值越小破碎化程度越高。计算结果如表 1所示。

流域整体景观的m_eff在1985年为6326.62 km²，2010年大幅减小为2974.32 km²，减小比例高达52.99%，说明25年间景观整体的破碎化程度大幅加剧。人口和经济快速发展，具有较大斑块面积的草地和未利用地发生转变，斑块面积进一步减小，破碎化程度加剧，景观异质性加大。

从流域分区来看，各流域分区的破碎化程度不一，破碎化最剧烈的是黄左区间，其m_eff仅占该区面积的6.96%，其后依次是黄右区间、引黄灌区、盐池内流区、清水河、苦水河和红柳沟。在20世纪80年代，流域人口少、经济发展水平低，人为因素对流域景观破碎化影响较大的区域主要在人口聚集区，即引黄灌区；此外，自然地理条件在一定程度上决定景观破碎化水平，在例如黄左区间和黄右区间等地形地貌复杂程度较大的区域，人为干扰较小，但自然景观复杂程度较大，也具有较高的景观破碎化水平；而苦水河和红柳沟地形地貌较为简单，景观类型单一，人为干扰也较小，破碎化程度最小。

由于各流域分区自然地理环境的不同，人口与经济的发展也表现出差异，对流域分区的景观破碎化的影响也不相同。从1985—2010年的m_eff变化可知，清水河、苦水河、盐池内流区、黄右区间、红柳沟和黄左区间的破碎化程度均有所加剧，其中清水河和盐池内流区的变化最为剧烈，m_eff减小幅度高达52.27%和56.48%；7个分区中仅有引黄灌区的m_eff有所增大，景观破碎化程度有所减弱。清水河流域大部分属于黄土高原，地形地貌复杂、建设用地和耕地的增加、草地的减少等引起景观复杂程度加大，破碎化程度加剧；盐池内流区的降水在研究时段有所增加^[21]，植被水分条件得到改善，并且“三北”防护林工程和防沙治沙工程等引起草地、林地面积增加，未利用地大幅减少，再加之耕地和建设用地的增长，景观破碎化程度增加；在人口最为密集的引黄灌区，对耕地和建设用地的需求加大，并连接成片，而草地和未利用地等面积大幅减少，使得景观类型单一化，导致破碎化程度有所降低。

3.2 流域景观破碎化的空间变异特征 3.2.1 特征尺度

以流域1985年有效粒度尺寸(m_eff)的空间尺度变异分析结果作为景观破碎化分析尺度确定的依据，结果如图 3所示。由图可知空间变异特征值(块金值 C₀ 与基台值(C₀ + C) 的比值)随幅度的增加而增大，在4500m左右时开始相对平稳，表明在该尺度能够反映景观破碎化空间变异特征，因此在Fragstats软件中设置4500m作为空间分析的“滑窗”大小。

3.2.2 景观破碎化指数的空间变异特征

流域m_eff大小的空间分布在4500m的空间尺度下如图 4所示，高值代表景观破碎化程度低，景观类型单一，且分布较连续，低值代表景观破碎化程度高，景观类型间的组合复杂。1985年和2010年流域有效粒度尺寸的高值区和低值区空间分布状况大体相同，高值区和低值区交错分布，且低值区范围明显多于高值区。破碎化程度较大的区域主要位于流域南部、中东部、黄河沿岸及贺兰山山麓地带，在流域分区上与清水河流域、盐池内流区、引黄灌区、黄左区间、黄右区间等相对应，景观破碎化的空间异质性与流域分区景观破碎化的程度保持一致。

就1985年和2010年流域m_eff的空间变化而言，m_eff减小的区域的范围明显大于m_eff增大的区域，即景观破碎化程度加大的区域多于下降的区域，景观破碎化空间异质性表现出明显的上升趋势，其中以流域南部和中东部的变化最为明显，结果与对应的清水河流域和盐池内流区景观破碎化的数量变化也保持一致。

3.3 不同背景下的景观破碎化异质性

景观破碎化在不同自然条件和人为干扰背景下的空间异质性不同。本文以海拔和降水作为自然背景条件，来探究景观破碎化在不同自然条件下的空间异质性(图 5)；由于受限于研究区人口和经济数据，本文仅以距城市和居民住宅地等人为活动最强区域的距离大小来表征流域内人为干扰的强弱^[22-23]，分析景观破碎化在不同人为干扰强度(图 5)下的空间异质性，即距建设用地距离越小的区段人为干扰越强。不同区段的景观破碎化程度以该区段所包含上文获取到随机点来计算平均有效粒度尺寸(avg-m_eff)表示，值越小，破碎化程度越高。

3.3.1 海拔变化下的景观破碎化异质性

不同海拔区段的avg-m_eff空间分异结果显示，1985年和2010年研究区内景观破碎化异质性特征基本一致。2400 m以下区段的景观破碎化程度均较高，且变化较为平稳；在2400 m以上区段的破碎化程度随海拔升高而明显降低(图 6)。表明2400m以下区段景观破碎化受海拔影响较小，以上区段景观受海拔因子的影响强烈。

就1985和2010年流域景观破碎化的变化来看，除1064—1200 m区段外，其他海拔区段的景观破碎化均有不同程度的加剧，在1200—2400 m区段破碎化程度加剧幅度随海拔的上升而增大，在2400 m以上区段景观破碎化加剧的幅度随海拔的上升而降低。在海拔背景条件一致情况下，流域2010年景观破碎化程度较1985年发生加剧变化，表明变化主要是由于其他自然条件变化和人为因素干扰引起的。

3.3.2 降水变化下的景观破碎化异质性

通过对流域不同降水量区域的有效粒度尺寸进行计算可知(图 6)，降水对流域景观破碎化变化有明显作用。在空间变化上，景观破碎化程度在1985年和2010年均表现出随降水的增加而加剧的变化，在多年平均降水量最大的400—425mm区段，破碎化程度达到最大。呈现这种变化是由流域降水的不均匀分布、流域自北向南的三级地形地貌分布以及人口经济等人为因素共同作用造成的。

景观破碎化程度在1985年到2010年流域各区段的破碎化程度均有不同程度的加剧。宁夏黄河流域近30年降水呈下降变化，气温逐渐升高，总体向暖干化发展变化^[24-25]，并会进一步影响流域的能量平衡，从而加速流域景观破碎化进程^[7]。此外，研究时段内流域人口和经济快速发展，不同降水区段人为影响均有所加剧，景观类型发生变化，也是导致流域破碎化程度加大的原因。

3.3.3 人为干扰下的景观破碎化异质性

由不同人为干扰强度区段的avg-m_eff可知(图 5，图 7)，1985年和2010年研究区内景观破碎化程度与人为干扰强度呈现出相同的变化趋势。在距建设用地1500 m以内的强人为干扰区段，斑块类型由较为单一建设用地逐渐向斑块类型多样的城市郊区变化，因此景观破碎化程度随人为干扰强度的减小而增大；在1500—10000 m中强度人为干扰区段，景观破碎化随随人为干扰强度的降低而逐渐减弱；在10000 m以外区段，破碎化程度又有明显的加剧变化。

就1985和2010年流域景观破碎化的变化来看，人为干扰较强的1500m以内和人为影响较弱的10000 m以外区段的景观破碎化程度均有所加剧。结合图 5可知，1500m以内区域的变化主要由受较强人为干扰引起，而10000 m以外区段可能受生态修复措施和降水、气温等自然条件变化的影响更为显著。而在1500—10000m之间的中度人为干扰区段的景观破碎化程度在1985年到2010年间基本未发生改变，表明在人为干扰强度相同的情况下，由自然条件改变引起景观破碎化发生的变化有限。

4 结论

对宁夏黄河流域景观破碎化进行研究，选取有效粒度尺寸(m_eff)作为景观破碎化程度的度量参数，对流域在1985—2010年间景观破碎化程度进行分析，并利用“滑窗法”对流域景观破碎化的空间变异特征进行分析，得到以下几点结论。

在研究时段内，宁夏黄河流域内黄左区间的景观破碎化程度最高，红柳沟和苦水河的破碎化程度最低，流域整体的景观破碎化在25年间急剧增加。

研究区内景观破碎化程度大的区域面积偏大，25年间景观破碎化程度加剧的区域增长显著，以流域南部和中东部最为明显，且与相应流域分区的数量变化保持一致。

景观破碎化在2400 m以下区段受海拔影响较小，以上区段具有随海拔升高而下降的效应；降水在空间上的不均匀分布及降水量的变化均对流域景观破碎化的空间异质性及变化有重要影响；在建设用地及近郊人为干扰较强区段的景观破碎化程度随人为干扰强度的增加而增大。

研究区在25年间景观破碎化受人为干扰影响产生的变化最为显著，由自然条件改变产生的影响有限。

在未来的流域整体的管理上，应结合流域的不同流域分区、自然以及不同人为干扰等因素对流域景观破碎化的影响来做出合理的利用规划，使其生态负面影响最小化。