柴达木盆地降水梯度下荒漠植被分布特征及影响因素
1.
2.
3.
4.
Vegetation distribution characteristics and influencing factors along a precipitation gradient in the Qaidam Basin
1.
2.
3.
4.
通讯作者:
收稿日期: 2025-10-30 修回日期: 2025-12-22
| 基金资助: |
|
Received: 2025-10-30 Revised: 2025-12-22
作者简介 About authors
康文蓉(1997—),女,山东莱芜人,博士研究生,主要从事生态水文研究.E-mail:
关键词:
Keywords:
本文引用格式
康文蓉, 赵文智, 张勇勇, 马登科, 马圆圆, 岳奕帆, 周国英.
Kang Wenrong, Zhao Wenzhi, Zhang Yongyong, Ma Denke, Ma Yuanyuan, Yue Yifan, Zhou Guoying.
0 引言
干旱半干旱区约占地球陆表面积的40%以上,是陆地生态系统中分布最广、对气候变化响应最敏感的区域之一[1]。其中,荒漠生态系统作为干旱区水分限制最突出的组成部分,其植被分布与结构特征对水分变化具有高度敏感性[2]。水分供给,尤其是降水的时空分布特征,是制约荒漠植被分布格局、群落结构及生态功能的首要环境因子[3]。沿降水梯度,植被覆盖度、群落组成、多样性及空间格局通常呈现出系统性变化,是认识干旱区生态系统结构与功能演变的重要切入点[4]。已有研究表明,降水通常通过调节土壤水分状况影响植物生长、生物量分配及物种组成,从而对群落多样性和空间格局产生重要影响。在多数情况下,降水增加能够促进物种丰富度和均匀度的提升,但这一响应关系往往受生态系统类型、土壤条件及物种组成等因素的共同调节,通常具有一定的时间滞后性[5-7]。同时,在不同气候背景和土壤条件下,植被对降水梯度的响应形式及其阈值特征仍存在差异,尤其是在极端干旱环境中,现有研究多依赖遥感反演或控制试验,对基于样地调查的自然降水梯度响应研究不足[6,8-9]。
随着遥感技术的发展,基于中高分辨率卫星影像的植被指数分析已成为研究干旱区植被动态及其对降水响应的重要手段,在区域尺度揭示植被覆盖变化及其时空格局方面具有显著优势[15-16]。已有研究表明,柴达木盆地植被绿度总体呈周缘升高、腹地偏低的空间格局,并在近几十年暖湿化背景下对降水变化表现出显著响应[13,17];植被物候、蒸散发及水分利用效率等生态特征亦随降水梯度和环境条件发生变化[18-20],如植被物候总体表现为生长季开始和结束时间推迟、生长期缩短,降水对物候的调控效应整体强于气温[18]。然而,受空间分辨率和混合像元效应限制,传统卫星遥感难以准确刻画荒漠植被斑块的精细结构特征,尤其是在植被稀疏、斑块尺度较小的极端干旱区,其对群落尺度空间异质性和斑块格局的表征能力仍然有限[21-22]。近年来,无人机遥感凭借厘米级空间分辨率、数据获取灵活等优势,在荒漠植被覆盖提取、斑块识别及空间格局量化等方面具有显著优势,为实现样方尺度与区域尺度之间的有效衔接提供了新的技术途径[23-24]。而目前针对柴达木盆地自然降水梯度下的群落多样性与斑块格局研究,仍存在空间与尺度覆盖不足的问题。既往研究多局限于局地生境或单一群落类型,缺少覆盖完整降水梯度的代表性样区或样带[13,25-27]。在柴达木盆地自然降水梯度背景下,将无人机获取的精细斑块格局指标与地面群落调查及土壤理化属性在相同空间尺度下进行综合分析的研究仍相对缺乏,一定程度上制约了对降水梯度下柴达木盆地荒漠植被与土壤及环境因素耦合机制的认识。
针对上述问题,本研究沿柴达木盆地年降水量200~300、100~200、50~100 mm梯度开展样地调查、无人机航测与土壤样品采集。研究目标为:①获取不同降水带的物种组成与多样性特征;②基于无人机影像量化斑块数量、聚集性与连通性等指标,揭示植被空间格局随降水梯度的响应规律;③测定并分析土壤理化属性,解析植被与环境因子的耦合关系与关键控制因子。研究旨在为高原荒漠区分区化治理和生态恢复配置提供科学依据与技术参考。
1 研究区概况与方法
1.1 研究区概况
研究区位于青藏高原东北部的柴达木盆地(图1),属典型内陆高原大陆性干旱区,水热时空分配不均、蒸发强烈,区域景观以荒漠为主。植被总体覆盖度较低,群落以耐旱灌木与多年生草本为主体,表现出对水分短缺的功能性状适应与生活型分化。研究按照年降水量将采样点划分为3个水分梯度(200~300、100~200、50~100 mm),以便在同一气候背景下比较不同水分约束强度下的群落与格局差异。
图1
图1
柴达木盆地不同降水梯度下样地分布及无人机影像示例
Fig.1
Distribution of sampling sites across different precipitation gradients in the Qaidam Basin and representative UAV images
1.2 样地选择与数据获取
样方植被调查:在3个降水梯度内(200~300、100~200、50~100 mm)各布设3个代表性样地。每个样地设置3个10 m×10 m的样方进行植被调查,主要记录群落内全部物种的种类与数量,测定各物种盖度、平均株高和冠幅,并按照生活型类别(乔木、灌木、多年生草本和一年生草本)进行划分,用于后续群落组成与多样性分析。
无人机样地植被格局:在每个样地选择3个50 m×50 m无人机航拍范围,利用DJ Phantom 4多光谱无人机,在2025年6月,选择晴朗无云的天气,于每日中午11:00—14:00光照充足时段进行低空摄影测量,飞行高度设置为20 m,并确保航线设计满足80%以上的旁向重叠率和航向重叠率,以获取高质量影像数据。所获取的多光谱影像具有5 cm×5 cm的空间分辨率,可满足植被斑块特征识别和提取需求。
1.3 数据分析与处理
群落多样性指标:基于样方数据计算多样性指数(Shannon-Wiener)、丰富度指数(Margalef)与均匀度指数(Evenness)。计算公式如下:
式中:
式中:DMg为丰富度指数;S为样方内记录到的物种数;N为总个体数。
式中:
景观格局指标:使用大疆智图对无人机影像进行前期处理,获得研究区正射影像图。首先通过目视解释对植被和裸地进行二元划分。将结果导入R语言[32]并转为分类栅格,使用terra包进行栅格读取、类型转换,按固定编码设置植被像元值为0、非植被为1。并以像元计数法计算植被盖度(%)。借助landscapemetrics包提取景观指数,包括斑块数量(NP)、最大斑块指数(LPI)和聚合度指数(AI)。
显著性与相关性检验:植被指标在不同降水区间的差异性采用Kruskal-Wallis检验;多重比较采用Bonferroni校正下的Dunn检验。植被指标(盖度、多样性指数、丰富度指数、均匀度指数)与环境变量和土壤理化性质之间的整体相关性采用Mantel检验(R语言vegan包实现)和Pearson相关性检验,并进一步在单因子层面检验环境因子与植被指标的相关性;显著性阈值均按P<0.05判定。
2 结果与分析
2.1 群落组成与多样性沿降水梯度的变化
表1 柴达木盆地荒漠植被分类及其在降水梯度下的分布
Table 1
| 生活型 | 物种名称 | 拉丁名 | 科 | 属 | 降水梯度划分/mm | ||
|---|---|---|---|---|---|---|---|
| 200~300 | 100~200 | 50~100 | |||||
一年生 草本 | 青蒿 | Artemisia caruifolia | 菊科 | 蒿属 | √ | ||
| 野葵 | Malva verticillata | 锦葵科 | 锦葵属 | √ | |||
| 碱蓬 | Suaeda glauca | 苋科 | 碱蓬属 | √ | |||
| 米蒿 | Artemisia dalai-lamae | 菊科 | 蒿属 | √ | |||
| 滨藜 | Atriplex patens | 苋科 | 滨藜属 | √ | |||
多年生 草本 | 披针叶野决明 | Thermopsis lanceolata | 豆科 | 野决明属 | √ | ||
| 芨芨草 | Neotrinia splendens | 禾本科 | 芨芨草属 | √ | |||
| 冰草 | Agropyron cristatum | 禾本科 | 冰草属 | √ | √ | ||
| 鸡冠茶 | Sibbaldianthe bifurca | 蔷薇科 | 毛莓草属 | √ | |||
| 苣荬菜 | Sonchus wightianus | 菊科 | 苦苣菜属 | √ | |||
| 蒙古韭 | Allium mongolicum | 石蒜科 | 葱属 | √ | |||
| 拐轴鸦葱 | Lipschitzia divaricata | 菊科 | 拐轴鸦葱属 | √ | |||
| 芦苇 | Phragmites australis | 禾本科 | 芦苇属 | √ | |||
| 针茅 | Stipa capillata | 禾本科 | 针茅属 | √ | √ | ||
| 灌木 | 黑沙蒿 | Artemisia ordosica | 菊科 | 蒿属 | √ | ||
| 珍珠柴 | Salsola passerina Bunge | 苋科 | 珍珠柴属 | √ | √ | ||
| 白刺 | Nitraria tangutorum | 白刺科 | 白刺属 | √ | √ | √ | |
| 红砂 | Reaumuria songarica | 柽柳科 | 红砂属 | √ | √ | ||
| 驼绒藜 | Krascheninnikovia ceratoides | 苋科 | 驼绒藜属 | √ | √ | √ | |
| 膜果麻黄 | Ephedra przewalskii | 麻黄科 | 麻黄属 | √ | |||
| 胶黄芪状棘豆 | Oxytropis tragacanthoides | 豆科 | 棘豆属 | √ | |||
| 柽柳 | Tamarix chinensis | 柽柳科 | 柽柳属 | √ | √ | ||
| 沙拐枣 | Calligonum mongolicum | 蓼科 | 沙拐枣属 | √ | |||
| 中亚紫菀木 | Asterothamnus centraliasiaticus | 菊科 | 紫菀木属 | √ | |||
| 梭梭 | Haloxylon ammodendron | 苋科 | 梭梭属 | √ | |||
图2
图2
基于盖度估算的不同降水梯度下植物科组成占比
Fig.2
Proportion of plant family composition based on coverage estimation under different precipitation gradients
图3
图3
基于盖度估算的不同降水梯度下植物生活型组成
Fig.3
Composition of plant life forms under different precipitation gradients based on coverage estimation
多样性指数进一步量化了上述差异(表2)。Shannon-Wiener指数在200~300 mm降水组显著高于100~200 mm降水组与50~100 mm降水组;物种数同样在200~300 mm降水组达到最大。均匀度指数与丰富度指数在3个降水区之间差异不显著。
表2 柴达木盆地不同降水梯度下物种多样性
Table 2
降水梯度 /mm | Shannon-Wiener 多样性指数 | Evenness 均匀度指数 | Marglef 丰富度指数 | 物种数 /个 |
|---|---|---|---|---|
| 200~300 | 1.14±0.32a | 0.72±0.12a | 0.05±0.04a | 5.11±1.54a |
| 100~200 | 0.93±0.74b | 0.51±0.38a | 0.08±0.06a | 4.56±2.96b |
| 50~100 | 0.36±0.40b | 0.38±0.40a | 0.08±0.11a | 1.92±1.16b |
2.2 植被空间格局沿降水梯度的分异特征
植被盖度在3个降水区之间两两差异显著(图4),表现为200~300 mm降水组最高、100~200 mm降水组次之、50~100 mm降水组最低,且随降水递减数据离散度增大。斑块数量在高降水组显著高于其余两组,而100~200、50~100 mm降水组之间差异不显著。聚合度指数与最大斑块指数在3个降水区之间无显著差异。
图4
图4
不同降水梯度下植被空间格局指标的变化
Fig.4
Changes in vegetation spatial pattern indices under different precipitation gradients
2.3 土壤与环境因子沿降水梯度的变化特征
表3 不同降水梯度下土壤理化性质的描述性统计
Table 3
| 降水梯度/mm | 全氮/(g·kg-1) | 标准差 | 变异系数 | ||
|---|---|---|---|---|---|
| 最小值 | 最大值 | 平均值 | |||
| 200~300 | 0.03 | 0.07 | 0.05 | 0.01 | 0.23 |
| 100~200 | 0.01 | 0.04 | 0.02 | 0.01 | 0.36 |
| 50~100 | 0.01 | 0.06 | 0.02 | 0.01 | 0.64 |
| 降水梯度/mm | 全磷/(g·kg-1) | 标准差 | 变异系数 | ||
| 最小值 | 最大值 | 平均值 | |||
| 200~300 | 0.06 | 0.68 | 0.47 | 0.21 | 0.44 |
| 100~200 | 0.37 | 1.23 | 0.60 | 0.19 | 0.31 |
| 50~100 | 0.04 | 0.60 | 0.37 | 0.16 | 0.45 |
| 降水梯度/mm | 全钾/(g·kg-1) | 标准差 | 变异系数 | ||
| 最小值 | 最大值 | 平均值 | |||
| 200~300 | 16.12 | 20.76 | 17.90 | 1.54 | 0.09 |
| 100~200 | 12.99 | 21.27 | 17.97 | 2.39 | 0.13 |
| 50~100 | 9.14 | 18.45 | 15.08 | 2.64 | 0.17 |
| 降水梯度/mm | pH | 标准差 | 变异系数 | ||
| 最小值 | 最大值 | 平均值 | |||
| 200~300 | 7.46 | 8.26 | 7.78 | 0.23 | 0.03 |
| 100~200 | 7.25 | 8.98 | 7.85 | 0.46 | 0.06 |
| 50~100 | 7.05 | 8.59 | 7.81 | 0.45 | 0.06 |
图5
图5
不同降水梯度对土壤全氮、全磷、全钾和pH的影响
Fig.5
Effects of different precipitation gradients on total nitrogen, total phosphorus, total potassium and pH of soil
图6
图6
植被-土壤-环境因子的相关性热图
Fig.6
Correlation heatmap of vegetation-soil-environmental factors
表4 植被指标与环境因子及土壤理化性质的Mantel相关分析
Table 4
| 因素 | r | P |
|---|---|---|
| 环境因素 | 0.19 | 0.006 |
| 土壤性质 | 0.17 | 0.051 |
图7
图7
单因素环境因子与主要植被指标的Mantel相关分析
Fig.7
Mantel correlation analysis of single-factor environmental factors and main vegetation indicators
3 讨论
3.1 群落多样性与植被空间格局对降水梯度的响应
在3个降水梯度上,灌木均为群落的优势生活型,但群落结构对降水变化表现出明显的响应,在200~300 mm降水组禾本科与一年生草本的相对多度显著增加,与灌木共同构成群落;而在50~100 mm降水组的极端干旱区间,群落则几乎完全由梭梭、白刺、驼绒藜等耐旱灌木所主导(表1,图2和图3)。这一变化趋势与低海拔干旱区的研究结果总体一致,如在宁夏荒漠草原的降水控制试验表明,随着降水减少,草本植物地上生物量显著下降,群落结构趋于简单,耐旱物种的重要性上升[32]。相关研究进一步表明,降水变化主要通过直接影响豆科和禾本科植物地上生物量进而间接影响植物群落地上生物量[8]。上述群落结构变化与荒漠植物的生存策略相一致,草本植物对水分变化更为敏感,在极度干旱条件下仅有少数耐旱性高的灌木能够维持稳定种群[20,33]。植物多样性随降水量的减少显著降低。200~300 mm降水组的Shannon-Wiener指数和物种丰富度明显高于其余降水组,但群落均匀度指数组间差异不显著,这可能是由于沿降水梯度的多样性下降主要源于物种丰富度的减少,而非优势度结构的显著改变[34]。
基于无人机影像的斑块指数分析表明,高降水区植被斑块面积较大、聚集度高,而低降水区植被斑块数量和聚集度降低,整体呈现出由连片化向破碎化转变的空间格局。这一特征符合荒漠区肥力岛效应,即在水分和养分双重限制下,植被倾向于在微地形和局地有利生境中形成聚集斑块,以集中利用有限资源;随着降水条件改善,这些斑块得以扩展并相互连接,从而提高植被盖度并维持较高的空间聚集性[35]。在50~100 mm降水组观测到盖度高变异性,反映了极端干旱条件下水分供给在空间上的高度不均一性。在该降水水平下,降水多以小尺度、脉冲式事件出现,通常仅能湿润表层土壤,且入渗深度受土壤质地限制[36],导致不同样地之间实际可利用水分差异显著。在此基础上,微地形对水分再分配的影响被显著放大,局地汇水条件较好的样地能够维持相对较高的植被盖度,而水分条件不好的样地则表现为极低盖度,从而造成样地间盖度差异的显著扩大。如已有研究表明,荒漠植被高覆盖区域往往沿局地低洼沟道或短暂性汇水区分布[37]。
3.2 土壤理化性质及环境因子对植被分布的影响
Mantel分析表明,植被指标与土壤理化性质在整体层面相关性不显著(表4)。这一结果表明,研究区所处的极端干旱背景下,土壤养分条件并非主导植被空间分布与群落结构变化的关键限制因子。研究区土壤养分含量整体处于较低水平,且不同样地间变异幅度相对有限,这在一定程度上削弱了土壤理化性质对植被格局的区分能力。与本研究结果不同,在年降水量为180、280、410 mm的较高降水梯度区研究表明,土壤理化性质与植被物种数、高度等指标呈现一致变化趋势[38]。这种差异可能是由于在水分相对充足的条件下,土壤养分更容易转化为限制植被生长和群落结构的重要因子;而在本研究区的极端干旱环境中,水分胁迫强度远高于养分限制,从而掩盖了土壤理化性质对植被的潜在影响。
4 结论
柴达木盆地降水梯度明显影响植物群落组成。随着降水减少,草本植物的比例逐渐降低,灌木始终为优势生活型,50~100 mm降水组群落结构变为单一灌木型。植被多样性随降水量降低而减弱,Shannon-Wiener指数在200~300 mm降水组显著高于50~100、100~200 mm降水组,其值分别为1.14±0.32、0.36±0.40、0.93±0.74。Evenness均匀度与Margalef丰富度指数在3组间差异不显著(P>0.05)。这表明在该降水梯度下的水分条件对现有物种间相对优势结构的影响相对较弱。空间格局则呈现出随降水递减的破碎趋势。低土壤养分条件下,土壤理化性质对植被格局的影响有限。全氮在200~300 mm降水组为0.05 g·kg-1;全磷与全钾在100~200 mm降水组最高,分别为0.60、17.97 g·kg-1;3组pH为7.05~8.26,组间差异均不显著(P>0.05)。综合分析表明植被指标与环境因子整体相关显著(Mantel,P<0.01),相关性强弱为地下水埋深>海拔>距离河道距离;与年均气温无显著相关(P>0.05)。
参考文献
Global desertification:building a science for dryland development
[J].
Desert ecosystems:environment and producers
[J].
Precipitation pulses and carbon fluxes in semiarid and arid ecosystems
[J].
Productivity is a poor predictor of plant species richness
[J].
Water-mediated responses of ecosystem carbon fluxes to climatic change in a temperate steppe
[J].
Red and photographic infrared linear combinations for monitoring vegetation
[J].
A continuous satellite-derived measure of global terrestrial primary production
[J].
Spatial and temporal variation of vegetation phenology and its response to climate changes in qaidam basin from 2000 to 2015
[J].
Distribution of actual evapotranspiration over qaidam basin, an arid area in China
[J].
Spatial variation in land cover and choice of spatial resolution for remote sensing
[J].
Remote sensing for biodiversity science and conservation
[J].
Lightweight unmanned aerial vehicles will revolutionize spatial ecology
[J].
High spatial resolution three-dimensional mapping of vegetation spectral dynamics using computer vision
[J].
柴达木盆地东部荒漠植物生态位特征
[J].
Total nitrogen analysis of soil and plant tissues
[J].
Part 3-chemical methods
[M].
Underground well water level observation grid dataset from 2005 to 2022
[J].
Precipitation changes alter the structure, species composition and interspecific relationships of desert steppe plant communities
[J].
Soil moisture and soluble salt content dominate changes in foliar δ 13C and δ 15N of desert communities in the Qaidam Basin,Qinghai-Tibetan Plateau
[J].
Precipitation and temperature regulate species diversity, plant coverage and aboveground biomass through opposing mechanisms in large-scale grasslands
[J].
The response of soil moisture to rainfall event size in subalpine grassland and meadows in a semi-arid mountain range:a case study in Northwestern China's Qilian Mountains
[J].
Dryland ecohydrology and climate change:critical issues and technical advances
[J].
Vegetation response to groundwater variation in arid environments:visualization of research evolution,synthesis of response types,and estimation of groundwater threshold
[J].
/
| 〈 |
|
〉 |

甘公网安备 62010202000688号