随着全球气候变暖, 干旱事件频发, 降水波动对荒漠生态系统的影响也日益加剧^[1-2]。据IPCC(2013)报告, 未来全球强降雨的强度和频率都会上升, 同时部分地区也将遭遇更为严重和频繁的旱灾^[3]。干旱事件将对生态系统的生物多样性^[4]、生产力^[5-7]、碳、氮循环^{[4, 8-10]}、土壤微生物群落^[11]等产生深远影响。目前, 为了揭示极端干旱事件等对不同生态系统的影响, 国际上已建立了IDE(International Drought Experiment)全球联网研究系统。而广泛分布于我国干旱和半干旱地区的荒漠生态系统由于其自身的脆弱性, 决定了其对极端干旱等气候条件的响应更为敏感^{[5, 8]}。因此, 深入研究降水波动等对荒漠区植物群落结构和功能的影响对于荒漠生态系统生物多样性保育和科学管理具有重要意义。

物种多样性和生物量是表征和评价植物群落组成和功能的重要指标^[12]。在干旱荒漠地区, 植物群落的物种多样性和生产力水平对于维持荒漠生态系统的功能稳定性, 以及改善区域生态环境具有重要作用^[13]。近年来, 国内外学者对物种多样性、生物量与降水变化的内在关系开展了深入的研究。大部分的研究侧重于探讨物种多样性、生产力与干旱、降水的年际变化之间的关系^{[1, 4, 14]}, 研究工作主要集中在草原、荒漠等地区^{[2, 7-8]}。例如, Gherardi和Sala在美国新墨西哥州荒漠通过6年的降水控制试验发现, 年际间降水变异系数的增加显著地降低了生态系统的初级生产力; 不同生活型植物对降水变异的响应截然相反, 其中草本的生物量下降了81%, 而灌木的生物量增加了67%。因此, 全球气候变化引起的降水量大幅波动会导致草本多度的下降, 而促进灌木的定居和生长^[1]。而极端干旱事件将对植物多样性产生严重的影响, 如在美国明尼苏达州经历了50年一遇的极端干旱事件后导致草地植物丰富度平均下降37%, 地上生物量下降47%^[4]。国内学者在降水变化对植物群落的影响方面也开展了研究, 叶鑫^[15]、杨秀静^[16]、徐满厚^[17]、孙羽^[18]等对高寒草地、草原、荒漠进行了降水调控对植物群落物种多样性和生物量影响的研究; 白春利等^[19]对内蒙古草原植物群落特征的研究表明, 水分添加会提高群落的丰富度指数、均匀度指数和多样性指数, 同时也发现植被群落特征的变化对水分的响应有赖于时间的积累。Wang等^[14]不同人工草地群落物种多样性和生产力在干旱和正常降水条件下的研究表明, 干旱和正常降雨条件下, 物种丰富度和生产力之间不相关; 生产力低的人工草地群落较生产力高的群落对干旱的耐受性更为显著。相对而言, 我国在干旱荒漠地区开展的相关研究较少, 亟需通过降水控制试验揭示极端干旱及降雨事件对荒漠生态系统植物多样性和生产力等的影响机制。

荒漠生态系统广泛分布于我国西北地区, 该地区年降水量稀少且时空分布不均, 蒸发强烈, 水资源相对匮乏, 属于典型的温带大陆性气候^[20]。本研究区域位于阿拉善荒漠南缘(地处宁夏自治区和内蒙古阿拉善盟交界处), 以草原化荒漠植被为主, 主要由灌木和草本层片组成, 物种多样性和植被覆盖度较低, 而且植被遭受破坏后生态恢复的进程非常缓慢或难以恢复^[21-22]。本研究按照IDE的标准布置了降水梯度控制试验平台, 以揭示不同降水处理(极端干旱、干旱以及降水增加)条件下荒漠植物多样性水平的变化以及不同生活型植物在生物量方面的响应规律, 研究工作可为全球气候变化背景下荒漠生态系统的科学管理提供一定的科学依据。

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

研究区位于中国科学院沙坡头沙漠研究站小红山荒漠生态系统野外综合观测场(37°29′N, 104°25′E), 地处宁夏中卫市境内的腾格里沙漠东南缘, 属于阿拉善高原荒漠与荒漠草原过渡带, 多年平均气温为10.6℃左右, 平均降雨量为186mm, 土壤为棕钙土, 土壤剖面分化较明显, 腐殖质层厚度为15—30cm。地下水埋深无法被植物根系直接吸收利用, 植物和土壤的主要水分来源于天然降水, 降水主要集中在7—9月^[23]。植物群落组成主要以细枝盐爪爪(Kalidium foliatum)、珍珠猪毛菜(Salsola passerina)和红砂(Reaumuria soongorica)为优势灌木, 以小画眉草(Eragrostis minor)、碱蓬(Suaeda glauca)为优势草本。

1.2 实验方法

2015年在宁夏中卫小红山地区, 布设了15个3m×3m的植被样地, 通过遮雨棚和滴灌系统来设计不同的降水控制条件(发明专利号:ZL201510186556.9), 包括5个降水梯度, 即极端干旱处理(-50%, D1)、干旱处理(-25%, D2)、对照(CK)、增雨处理Ⅰ(+50%, W1)、增雨处理Ⅱ(+25%, W2), 每个处理3个重复。同期, 在试验区内架设小型气象站, 主要对降水量、空气温度、有效光合辐射(PAR)、土壤10cm和40cm深度的含水量、温度和电导率进行连续数据采集, 采用CR1000数采仪记录数据。经过一年观测研究, 于2016年9月末, 进行物种组成及生物量等群落特征的调查。在每个样方内对所有灌木进行调查, 测量并记录每一株灌木的冠幅和株高; 草本调查是在每个水分处理样方中随机选择3个50cm×50cm的小样方进行调查, 共计45个样方, 草本地上生物量采用刈割法取样, 然后在烘箱中75℃恒温烘干至恒重。

1.3 数据处理与分析

植物群落物种多样性指数计算公式如下^[24]:

丰富度指数:

Margalef指数:

Simpson多样性指数:

Shannon-Wiener多样性指数:

Pielou均匀度指数:

式中, S为样方中的物种数, p_i为第i中植物的株数占总株数的比例, N植物总株数。

为保证试验样地的长期定位连续监测, 减少人为干扰和破环, 研究前期通过大样本采样的方式建立了优势灌木冠幅、高度和生物量之间的优化估算模型^[25]。具体的模型及参数见表 1。

利用SPSS 19.0软件进行单因素方差分析, 比较不同降水处理对植物多样性和地上生物量的影响, 差异显著条件下采用LSD法进行多重比较; 利用回归曲线拟合不同降水处理和地上生物量的相关性。采用Origin 8.5软件进行绘图。

2 结果分析 2.1 不同降水控制处理的土壤水分动态变化

降水控制对土壤含水量有着一定的影响作用, 其中D1、D2、CK、W1、W2处理的10cm和40cm平均土壤体积含水量(v/v%)分别为11.38%±0.89%、8.96%±0.99%、9.57%±0.71%、14.36%±1.06%、11.97%±0.89%和9.08%±0.46%、12.39%±0.67%、14.52%±0.89%、21.39%±1.34%、18.15%±1.11%。通过土壤深度、降水控制条件、月份三因素对土壤含水量影响的方差分析发现, 土壤深度、降水控制条件、月份, 以及月份和降水控制条件的交互作用对土壤含水量的影响显著(P＜0.01), 土壤深度与月份的交互作用、降水控制条件与月份的交互作用对土壤含水量的影响不显著(P>0.05)(表 2)。不同深度土壤含水量与降水量、温度的变化趋势一致, 但10cm土壤含水量显著低于对40cm土壤含水量(图 1, 图 2)。

2.2 不同降水控制对群落物种多样性的影响 2.2.1 灌木层片物种多样性

选取不同多样性指数对灌木层片物种多样性进行计算, 分析结果表明不同实验处理下灌木多样性指数差异不显著(P>0.05)(表 3)。Margalef指数、Simpson指数、Shannon-Wiener指数最大值均出现在极端干旱处理(D1), 分别为0.86±0.14、0.47±0.09、0.88±0.16;Pielou均匀度指数随水分的添加呈现递增的趋势, 但差异性不显著(P>0.05), 最大值出现在增雨处理Ⅰ(W1), 为0.69±0.10。

2.2.2 草本层片物种多样性

对草本层片物种多样性指数进行单因素方差分析(表 4)发现, 在不同降水控制条件下, Margalef指数和Pielou均匀度指数差异性显著(P＜0.05), 但并未随着水分的添加出现增加的趋势, 最大值均出现在干旱处理(D2), 分别为0.48±0.04、0.82±0.04;Simpson指数和Shannon-Wiener指数差异性极显著(P＜0.01), 最大值出现在干旱处理(D2), 分别为0.50±0.04、0.79±0.06。

2.3 降水控制对群落生物量的影响 2.3.1 灌木层片生物量

通过对灌木地上生物量进行单因素方差分析表明(表 5), 不同降水控制处理对灌木地上生物量影响不显著。灌木地上生物量受水分影响差异不显著, 其总体变化是随着土壤水分含量的增加呈现先升高后降低的趋势, 最大值出现在干旱处理(D2), 为(369.39±12.41)g/m²; 珍珠猪毛菜地上生物量受水分影响差异不显著(P>0.05), 变化趋势与总体生物量一致, 最大值出现在极端干旱处理(D1), 为(146.83±67.49)g/m²; 盐爪爪地上生物量受水分影响差异显著(P＜0.05), 但生物量并未随水分的添加呈现增多的趋势, 最大值也出现在干旱处理(D2), 为(224.12±19.97)g/m²。

2.3.2 草本层片生物量

不同降水控制条件对草本植物生物量的影响显著(表 6), 草本总体生物量随水分的增加呈现递增的变化趋势, 最大值出现在增雨处理Ⅰ(W1), 为(105.69±28.60)g/m², 总体差异极显著(P＜0.01)。对于草本优势种碱蓬和小画眉草, 在不同水分条件控制下, 碱蓬的生物量变化与总体生物量变化趋势一致, 碱蓬最大值出现在增雨处理Ⅰ(W1), 为(62.45±12.15)g/m², 生物量与水分含量呈现正相关性, 且差异极显著(P＜0.01)。而小画眉草生物量的变化趋势虽然随水分的增加而递增, 但差异性不显著(P>0.05), 最大值出现在增雨处理Ⅱ(W2), 为(0.98±0.15)g/m²(图 3)。总体而言, 优势种碱蓬对总体生物量的变化有着关键作用, 而小画眉草对总体生物量的变化影响作用较小。

2.4 群落特征与土壤含水量之间的关系 2.4.1 物种多样性与土壤含水量的关系

物种多样性与土壤含水量的相关性分析表明(表 7), 灌木物种多样性和10、40cm深度的土壤含水量没有显著相关性; 草本物种多样性和40cm土壤含水量不相关, 仅和表层10cm处土壤含水量呈弱负相关性(P=0.07)。

2.4.2 生物量与土壤含水量的关系

通过对生物量与土壤含水量相关性分析表明, 4种优势植物的生物量与10cm深土壤含水量之间的相关性均不显著(P>0.05), 40cm土壤含水量主要和草本生物量、碱蓬生物量相关性显著(P＜0.05)(表 8)。

3 讨论 3.1 水分控制对灌木和草本层片物种多样性的影响

对于我国西北内陆干旱荒漠生态系统, 水分是植物生长的主要限制性因素, 极端干旱和降水的大幅波动会通过土壤水分的变化在不同程度上影响荒漠植物群落多样性水平以及养分循环过程^{[2, 4, 26]}。本研究结果表明, 在降水逐渐增加的条件下, 灌木和草本层片的物种多样性指数呈现出先升高后降低的变化形式, 但差异性均不显著(P>0.05)。这与赵新风等^[25]、曾勇等^[26]的研究结果相似, 即灌木的物种多样性对不同水分控制条件表现出的差异不显著。而白春利等^[19]在内蒙古草原开展的水分添加试验表明, 植物丰富度和物种多样性水平会随着水分增加呈上升趋势。可见, 植物群落对于干旱的响应很大程度上取决于群落结构、优势种的生活史和生物地化循环过程等^[27-28]。究其原因, 研究认为荒漠植物多样性水平的显著变化是降水等环境因素长期变化和作用的结果, 短期(例如1年)的水分变化对灌木和草本层片的物种多样性影响不显著。荒漠灌木在长期的进化适应过程中, 对于水分的胁迫已经具有一定的生理耐受性, 灌木的丰富度和多度的响应可能存在“滞后效应”, 还需要长期定位研究才能给出确切的答案。本文仅是总结了干旱处理一年的试验研究结果, 目前的结论也仅是反映短期干旱处理对荒漠植被物种多样性水平的影响。同样, 草本的物种多样性对不同的降水控制处理也不存在显著差异, 短期内发生显著变化的可能性较小。

3.2 水分控制对灌木和草本层片生物量的影响

在荒漠地区, 生物量作为表征植被长势重要指标对降水变化的响应尤为敏感。本研究发现降水变化对珍珠猪毛菜地上生物量和灌木总体生物量的影响不显著(P>0.05)。相关研究也得到相同的结果, 例如在半干旱区对灌木群落进行长期控水处理后, 刘峻杉等^[29]等也发现灌木生物量并未发生明显的变化。可见, 为适应荒漠地区降水的年际波动, 灌木在长期的进化进程中不断适应, 已形成了适应干旱等胁迫环境生理生态特征。此外, 对不同人工草地群落物种多样性和生产力在干旱和正常降水条件下的研究表明, 生产力低的人工草地群落较生产力高的群落对干旱的耐受性更为显著^[14]。降水控制对草本总体生物量和碱蓬生物量影响极显著(P＜0.01), 而对于小画眉草的影响不显著(P>0.05), 但均随水分添加呈现出增加的趋势。这与李文娇^[30], 白春利^[19]等人的研究相似, 单一的水分添加可以增加草本植物的生物量。荒漠生态系统的生产力相对较低, 在受到干旱胁迫时生物量不会发生大幅下降。但是, 从荒漠植物生活型及其利用水分来源的角度分析, 深根系的灌木主要利用深层土壤水分, 因此, 生物量不会发生大幅波动; 而浅根系一年生或多年生草本生物量的积累主要和浅表层的土壤水分密切相关。因此, 对降水波动的响应更为敏感。此外, 草本植物生物量的积累也和水分利用效率等生理特征密切相关, 例如本研究中碱蓬作为藜科C₄植物在水分竞争方面较小画眉草更具优势, 因此对水分的响应更为敏感。

3.3 物种多样性和生物量变化与土壤含水量的关系

植物多样性是衡量群落内物种分布的均匀程度和数量的一个指标, 它与植被生长的环境有密切的关系, 是群落特征的综合体现。本研究结果表明, 灌木和草本物种多样性与土壤含水量之间的相关性均未达到显著水平(P>0.05)。特定群落内植物多样性水平是长期进化和对环境不断适应的体现, 在一定的土壤水分范围内, 群落内植物的丰富度水平会保持在一定范围, 受到干旱和降水波动的影响, 植物的多度水平首先会作出一定的响应。当然, 极端气候事件也会导致植物丰富度下降, 部分稀有种也可能会面临灭绝的危险^[4]。例如在美国明尼苏达州经历了50年一遇的极端干旱事件后导致草地植物丰富度平均下降37%, 地上生物量下降47%^[4]。国内学者对黄土高原天然草地群落结果与土壤水分关系的研究表明, 物种多样性与土壤含水量之间为正相关关系^[31]。本研究区域为草原化荒漠地区, 荒漠植物普遍具有一定的抗旱能力, 土壤水分的波动在植物耐受范围内时不会对植物多样性水平产生影响。在经历极端气候事件后, 长期的监测研究才能对植物多样性和水分之间的关系给出科学合理的回答。

本研究区优势灌木的生物量与土壤含水量(10cm和40cm深度)的相关性不显著。荒漠灌木根系发达, 主根可达到80—100cm, 后期可通过对更深层土壤含水量的测定明确二者之间关系。而草本植物的根系主要分布在土壤浅表层, 因此浅表层土壤含水量的变化直接影响草本层片的生物量。本研究表明, 降水与土壤含水量的年变化一致, 表层(10cm)土壤含水量对草本生物量的影响不显著, 浅层(40cm)土壤含水量对草本生物量的影响显著(表 7)。究其原因, 荒漠土壤土质疏松、毛管作用弱、通透性好, 导致土壤水分垂直分布差异较大, 其中表层的水分不稳定, 容易蒸发, 不易被植物根系利用; 深层土壤对有效降水具有“累积效应”, 水分相对稳定。在美国新墨西哥州荒漠的降水控制试验发现, 草本对于降雨的年际波动表现出“饱和效应”, 即干旱年份对草本层片产生的负效应要明显高于湿润年份的降水增加引起的正效应, 因此, 年际水平上降水的大幅波动会降低草本的生产力。当全球变化引起降水量大幅波动时会导致草本多度的下降, 而促进灌木的定居和生长^[1]。

4 结论

短期(1年)干旱和降水增加处理对灌木和草本层片的物种多样性以及灌木生物量的影响不显著; 但对草本生物量有显著影响。草本生物量主要取决于浅层(40cm)的土壤水分状况。