Page 100 - 2023年第54卷第11期
P. 100
大 NDVI值。NDVI计算公式如下:
NDVI = (NIR - Red)?(NIR + Red) (2)
式中:NIR为 Landsat影像近红外波段(0.76~0.96μ m)的反射率;Red为 Landsat影像红色波段(0.62~
0.69μ m)的反射率;NDVI的数值范围为 - 1~1 ,数值越高表示地物为植被的可能性越大,通常 NDVI
大于 0被认为是植被。
植被覆盖度 FVC运用以下公式计算得到:
FVC = (NDVI - NDVI )?(NDVI - NDVI ) (3)
soil veg soil
式中:NDVI 为地表完全覆盖植被时的最大 NDVI值;NDVI 为地表纯裸地(或无植被覆盖)时的最小 NDVI值。
soil
veg
景观指数的计算以 1990—2020年土地利用数据为基础。土地利用数据来源于逐年土地覆被产品
( Landsat - derivedannuallandcoverproductofChina(CLCD),https:??zenodo.org?record?5816591) [28] ,空
间分辨率为 30m,时间序列为 1990—2020年。水系功能连通选取的景观多样性指标包括香浓多样性
指数 SHDI和辛普森多样性指数 SIDI,计算公式如下:
n
SHDI =- ∑ plnp i (4)
i
i =1
n
∑
SIDI =1 - p 2 (5)
i
i =1
式中:n为景观中斑块类型总数;p为斑块类型 i在景观中出现的频率。
i
实测数据包括径流、水库与生态闸相关数据,均来源于塔里木河流域管理局。数据清单如下:径
流数据包括 1990—2020年阿拉尔、新其满、英巴扎、乌斯满、恰拉水文站月径流量、大西海子水库
日下泄流量等;水库相关数据包括水库分布、数量、库容、水库建成时间等;生态闸数据包括生态闸
数量、分布、建成时间、设计流量、灌溉面积等。
2.2 研究方法
2.2.1 基于拓扑空间网络的形态连通 根据水系连通特点,本文从形态、结构和功能三方面分析水系
连通对输水工程的响应。其中,形态响应主要运用图论法和景观生态学理论,将河网水系抽象成节点、
廊道、斑块,以此构建拓扑空间网络表征河流外观连通情况 [18] ,来反映河网发育程度与形态特征,这是
流域物质循环、能量流动、生物迁移的物理基础。具体而言,水系形态连通响应指标有节点数量、廊道
数量、斑块数量、水系环度、节点连接率、河网连接度、生态网络连通度等 [14,18] ,见表 1。
表 1 水系形态连通指标与计算公式
准则层 指标层 权重 公式 指标说明
节点数量 V 0.13 V V为节点数量,河流交汇处即为节点
廊道数量 N 0.13 N N为河流廊道数量
斑块数量 PN 0.14 PN为景观中包含的斑块数量
N - V + 1
水系环度 α 0.15 α = α为水系环度,区域内水系连接成环的水平
2 V - 5
N β 为节点连接率,0 ≤ β≤ 3; β≤ 1时,水 系 为 树 枝 状 结 构; β >1
形态连通指标 节点连接率 β 0.15 β =
V 时,水系为格状结构或回路结构
N γ为河网连接度,水系之间形成连通通道的能力,0 ≤γ≤1,0表
河网连接度 γ 0.15 γ =
3(V - 2 ) 示节点间没有连线,1表示每个节点间都相互连通
e C为景观生态网络的连接度,反映区域内各节点依靠河流廊道相
生态网络连接度 C 0.15 ξ 互连通的程度;e为河流廊道总长度;A为区域面积;ξ 为水系网
C =
槡 VA 络变形系数,是各节点间实际河流廊道总长度与直线总长度之比
2.2.2 基于河流连续体的时空维结构连通 水系结构连通响应主要从空间和时间多个维度表征水系连
通时空演变动态,包括空间上的横向、纵向与时间变化过程 [18] 。横向连通主要指河道- 洪泛区在垂直
于河流流向上的连通性,纵向连通主要指河源- 尾闾湖?河口沿河流流向上的连通性;时间上,水系连
通往往会随丰枯水期呈周期性节律变化。本文使用的结构连通指标与计算公式见表 2。
6
— 1 3 2 —