鄱阳湖湿地土壤微生物群落结构沿地下水位梯度分异特征

张广帅,于秀波,张全军,李　雅,刘　宇,段后浪

1 中国科学院地理科学与资源研究所,中国科学院生态系统网络观测与模拟重点实验室,北京　100101 2 中国科学院大学,北京　100049

湿地是陆地生态系统和水域生态系统相互作用形成的自然综合体和独特生态系统类型,在维护区域生态平衡和生物多样性保护等方面具有重要作用[1],故而被称为“地球之肾”和“生态超市”[2]。由于湿地生态系统独特的环境条件和丰富的物种多样性,物质交换和能量流动过程十分复杂,因而其生物地球化学过程相对于陆地生态系统更为复杂[3]。

土壤微生物是使土壤具有生命力的最主要成分,与土壤肥力和土壤健康水平密切相关[4]。湿地生态系统中,土壤微生物不仅在有机质分解、腐殖质合成以及物质转化和能量传递过程中发挥关键作用[5,6],还控制着沼泽湿地类型的分异与演替和土壤碳库的动态变化[7]。在洪泛平原湿地生态系统中,土壤微生物群落结构对环境因子的变化尤其是土壤湿度、pH、质地结构和营养元素的可利用程度响应敏感[8]。水文情势的变化是湿地生态系统结构和功能的主要控制因素[9],湿地水文状态的改变不仅制约着土壤微生物的生理和功能代谢水平,还影响微生物功能群落的组成[10]。地下水位变化所引起的干湿交替环境通过其有氧和无氧环境的转变为土壤微生物的群落组成形成了生理压力[11]。微生物群落结构沿地下水位梯度的分异与土壤水分、透气性、有机养分和植被等因子的差异密切相关[12]。高地下水位梯度的湿润条件下,细菌更具有生存优势[13]。

鄱阳湖是我国第一大淡水湖泊,受流域来水和长江水位共同作用,年内洪枯水位可达12m,这种独特的水文节律变化下形成了大面积干湿交替的洲滩湿地[14]。沿不同的高程梯度,地下水位变化呈明显的空间地带性分布,进而在洲滩尺度上形成了不同的微地带性环境。相关研究表明,水文情势的变化对鄱阳湖植被群落分布、生物量积累和物种多样性造成了不同程度的影响[15-16]。目前在鄱阳湖湿地的研究主要集中在水文过程与景观变化[17-18]、土壤微生物与植被关系[19-20]等方面,这些研究对认识鄱阳湖湿地的变化规律奠定了基础。Coulis等[21]认为短期内生态环境因子的变化不会直接改变微生物的群落结构,只能改变微生物活性,枯水期内洲滩尺度上不同地下水位环境梯度是否会引起土壤微生物群落结构的变化及其与微环境因子的相关性研究尚不足。本文以鄱阳湖典型碟形子湖白沙湖洲滩湿地为研究对象,研究比较不同地下水位梯度土壤湿地微生物生物量和群落结构的分异特征,并对相关影响因素进行分析,以期为深入了解枯水期鄱阳湖湿地生态系统结构稳定性的维持机制及湿地生物地球化学循环过程,揭示土壤地下生物群系应对枯水期环境变化的响应策略提供科学依据,进而对鄱阳湖湿地生态系统的有效保护具有重要意义。

1　材料与方法

1.1　研究区概况

鄱阳湖地处江西北部,长江中下游,地理坐标为115°49′—116°46′E,28°24′—29°46′N,南北相隔170km,东西最大跨度70km。鄱阳湖属东亚季风区,为典型亚热带季风气候,夏季炎热多雨,冬季低温少雨,年平均降水量为1387—1795mm,年内、年际变化较大。年降水量最大2452.8mm(1954年),最小1082.6mm(1978年),最大6个月降水量(3—8月)占全年降水量的74.4%。年平均蒸发量800—1200mm,约有一半集中在温度最高的7—9月。多年平均气温17.6℃,最冷月(1月)平均气温5.1℃,最热月(7月)平均气温29.5℃。受鄱阳湖水文节律季节性变化的影响,呈现出河、湖、滩交错的湿地景观格局。

研究试验样地位于江西南矶湿地国家级自然保护区内的白沙湖,地处赣江三角洲前缘,属典型的内陆河口湿地区。白沙湖为典型的碟形湖泊,洲滩地势平坦,受鄱阳湖季节性周期水文变化影响,4—9月份保护区除南山岛和矶山岛外,其余地方均被淹没,白沙湖与大湖连通,10月份湖水逐渐消退,归入碟形洼地,不同高程的洲滩相继出露。白沙湖洲滩发育有大面积的草洲,优势植被为灰化苔草(Carexcinerascens)、南荻(Triarrhenalutarioriparia)、芦苇(Phragmitesaustralis)、虉草(Phalarisarundinacea)等,其中灰化苔草分布最为广泛,几乎遍及整个湿地洲滩,盖度在85%—100%。试验样地位置如图1所示,面积为200m×300m, 样地内优势植被为灰化苔草,伴生种为南荻和芦苇。从湖岸向湖心按照地下水位划分为4个梯度(GT-A,GT-B,GT-C,GT-D),梯度带之间相隔约100m。不同梯度带灰化苔草生长态势、土壤温度及地下水位分异明显(表1)。每个地下水位梯度带内根据微地形差异设置3—6个采样点,同一梯度带内采样点间隔50—80m,每个采样点设置3个重复。每个采样点处埋入细口径(内径为2cm)PVC管(埋入地下约70cm)用于测定地下水位。

1.2　样品采集与分析

2016年11月下旬,利用土钻(内径5cm)对每个采样点进行五点混合取样,采集0—20cm土壤样品约500g装入聚乙烯自封袋内,编号后置于保鲜盒运回实验室进行后续分析。土壤样品运回实验室后,挑拣出石块等杂物后,每份样品分为两份,一份置于-20℃冷藏箱中用于测定土壤微生物磷脂脂肪酸(PLFA)测定,另一部分自然风干后磨碎风干,分别过20目和100目筛,用于土壤理化性质的分析。

图1　研究区地理位置与样地设置Fig.1　Location of study area GT-A, GT-B, GT-C, GT-D 代表不同的地下水位梯度带(从湖岸向湖心方向依次为A, B, C, D)

土壤理化性质分析具体方法如下：土壤含水量采用烘干法测定,土壤颗粒组成采样马尔文激光粒度仪测定；土壤容重(BD)采用环刀法测定；土壤pH采用水土比2.5∶1 pH计测定；土壤有机碳(SOC)采用重铬酸钾外加热法测定；土壤总氮(TN)采用半微量凯氏定氮法测定；土壤总磷(TP)采用钼锑抗比色法测定。土壤微生物量碳(MBC)和微生物量氮(MBN)的测定采用氯仿熏蒸-K2SO4浸提法[22]。

表1　采样点基本信息

*采样点生物量为苔草生物量, GT-A, GT-B, GT-C, GT-D 代表不同的地下水位梯度带(从湖岸向湖心方向依次为A, B, C, D)

土壤微生物磷脂脂肪酸(PLFA)的分析测定参照Frostegård 等[23]的KOH-甲醇溶液甲酯化方法进行,以十九烷酸为内标采用GC-MS (TRACE GC Ultra ISQ) 测定,全自动进样装置,检测磷脂脂肪酸的方法为：进样口温度280℃,分流比20∶1,柱温箱程序升温为70℃起始,保持1min,以20℃/min升温至170℃,维持2min,再以5℃/min升温至280℃,维持5min,最后以40℃/min升温至3000℃,维持1.5min。

对不同地下水位条件下的土壤进行PLFA 分析,共检测出30种不同的脂肪酸。通常情况下每类特征标记脂肪酸仅存在于特定的微生物类群内,且占有优势地位[24]。因此不同的特征标记脂肪酸可以表征特定的生物类群。其中,i13:0, i14:0, i15:0, a15:0, i16:0, i17:0, a17:0, a16:0, i15:1 可以表征革兰氏阳性细菌；12:0,14:1ω5c, 14:0, 16:1ω6c, 15:0, cy17:0, 16:1 2OH, i17:0 3OH, cy19:0 ω8c, i15:0 3OH, 16:1ω9c 可以表征革兰氏阴性细菌；16:1ω5c,18:1ω9c,18:3ω6c(6,9,12)表征真菌；10Me17:0 和10Me 18:0 表征放线菌[25-27]。

1.3　数据处理

采用SPSS 20.0 进行数据的统计分析,采用单因素方差分析不同地下水位梯度环境下土壤理化性质、微生物群落结构的变化,多重比较方法选用Dunnett′s T3比较检验,采用双变量皮尔逊相关分析研究区微生物群落结构与土壤理化性质的相关性；采用CANOCO(Canoco for Windows 5.0)软件对不同地下水位环境梯度下不同采样点土壤微生物特征PLFAs进行主成分分析(PCA),对土壤理化性质、土壤微生物群落结构之间的关系进行冗余分析(RDA),RDA分析通过蒙特卡洛显著性检验(Monte Carlo permutation text)(P=0.02)。

根据脂肪酸含量和种类计算微生物群落Shannon 多样性指数(H′)、Pielou 均匀度指数(J)和 Simpson优势度指数(D)来评价微生物多样性特征[28-30],计算公式如下：

H′=-∑PilnPi

式中,Pi=Ni/N,Ni为 第i类脂肪酸量,N为试验检测出的总脂肪酸量,S为微生物群落中总脂肪酸的种类。

2　结果与分析

2.1　土壤属性沿地下水位梯度的分异

地下水位对TN以外的土壤理化性质影响显著(表2；P< 0.05)。随着地下水位的升高,土壤pH从4.62(GT-A)升至5.05(GT-D),增长幅度达9.3%(P< 0.05)；土壤含水量从30.8%(GT-A)增长到52.6%(GT-D),变幅达70.7%(P< 0.05)。GT-D的粘粒含量比GT-A、GT-B和GT-C分别降低了34.6%、31.2%、22.7%,粉粒含量分别降低了7.0%、4.6%和1.5%,而沙粒含量则增加了4.28、2.18和1.28倍(P< 0.05)。地下水位的升高促进了土壤TOC 和TP的流失,相对GT-A,GT-D 的TOC 和TP含量显著降低了53.5%和26.1%(P< 0.05)。不同地下水位梯度间,土壤TN差异性不显著(P> 0.05),土壤碳氮比依次表现为GT-A(10.78)> GT-B(10.20)> GT-C(9.07)> GT-D(7.63)。

表2　不同地下水位梯度土壤理化性质

不同小写字母表示同行在不同地下水位梯度带间在0.05水平上具有显著差异

2.2　地下水位对土壤微生物量及其分配比例的影响

图2　不同地下水位梯度土壤微生物MBC、MBN及其分配比例Fig.2　MBC, MBN, allocation ratio of MBC to SOC and allocation ration of MBN to TN of different ground water level gradient zonesMBC: 微生物量碳, Microbial biomass carbon；MBN：微生物量氮,Microbial biomass nitrogen； SOC：土壤有机碳,Soil organic carbon; TN: 土壤总氮, Total nitrogen

土壤生物量碳氮在不同地下水位梯度间表现为在GT-A、GT-B和GT-C之间差异不显著(图2),GT-D显著高于其他3个梯度(P< 0.05)。GT-D的微生物量碳氮分别比GT-A增长了2.8倍和4.3倍。微生物商是指土壤微生物量碳和土壤有机碳含量的百分比(MBC/SOC)。随着地下水位梯度的升高,微生物商逐渐升高,在GT-A、GT-B和GT-C之间变化不显著,相对GT-C,GT-D的微生物商升高了2.6倍。土壤微生物量氮的分配比例(MBN/TN)在不同地下水位梯度间的变化趋势与微生物商一致,GT-D 土壤微生物量氮的分配比例是GT-C的3.95倍。

2.3　地下水位对土壤微生物群落结构的影响

不同地下水位土壤微生物类群差异性显著(图3)。GT-D显著增强了微生物总PLFAs,细菌PLFA和放线菌PLFA的含量(P< 0.05),与GT-A相比,增长幅度分别为106.8%、117.2%、74.9%,其中革兰氏阳性细菌和革兰氏阴性细菌分别增长了107.9%和207.2%。真菌PLFA随着地下水位的升高呈现出低-高-低-高的波动性变化。由GT-A到GT-D,细菌PLFA占总PLFAs的比重分别为79.5%、80.9%、78.3%、83.5%,真菌PLFA占总PLFAs的比重分别为12.9%、12.9%、11.3%、8.0%,放线菌PLFA占总PLFA的比重分别为5.8%、6.2%、7.8%、4.9%。不同地下水位梯度样点土壤微生物PLFA数据PCA分析表明,沿着第1、第2主排序轴(PC1),采样点产生了明显的空间分离,GT-D和GT-A、GT-B、GT-C沿第1主排序轴能明显分开；而GT-A、GT-B、GTC之间则沿第2主排序轴明显分开(图4)。PC1和PC2分别解释了空间分异的76.7%和13.6%,累积贡献率达90.3%。在所有特征PLFA中,对第1主成分起主要作用的是i13:0、i15:0、i17:0 3OH、i15:0 3OH、a16:0; 而18:3ω6c、10Me 17:0 对第2主成分贡献较大。

图3　不同地下水位梯度土壤的总PLFA,细菌,真菌,放线菌,革兰氏阳性细菌,革兰氏阴性细菌Fig.3　The amount of soil PLFAs (total PLFAs, total bacteria, gram-positive bacteria, gram-negative bacteria, Actinomycetes and Fungi of different ground water level gradient zones不同小写字母表示在0.05水平上具有显著差异

图4　PLFA 指纹图谱主成分分析(A)土壤微生物群落,(B)磷脂脂肪酸Fig.4　Principle component analysis of PLFA profiles from soil microbial communities,(A)soil microbial community structure, (B) phospholipid fatty acids contributing to soil microbial community ordination pattern

在地下水位梯度升高的过程中,真菌/细菌(F/B)由0.26(GT-A)降至0.11(GT-D),压力响应指数即支链反异构饱和脂肪酸与支链异构饱和脂肪酸的比值(SA/SI)由0.47(GT-A)升至0.79(GT-D)。4个地下水位梯度带中,GT-D均具有最低的革兰氏阳性细菌/革兰氏阴性细菌(G(+)/G(-))和单一不饱和脂肪酸/支链脂肪酸(M/B)(表3；P< 0.05)。

表3不同地下水位梯度土壤微生物群落对环境的响应特征指数

Table3Fungal∶BacterialPLFAS(F/B),Gram+∶Gram-PLFAs(G+/G-),Themonounsaturated∶BranchedPLFAs(M/B),Sumofanteiso∶SumofisobranchedPLFAs(SA/SI)ofdifferentgroundwaterlevelgradients

梯度带GradientzonesG(+)/G(-)F/BM/BSA/SIGT-A2．12±0．22b0．17±0．02b0．26±0．02b0．47±0．05a GT-B2．06±0．26b0．16±0．02b0．30±0．06b0．53±0．08a GT-C2．36±0．69b0．14±0．03b0．18±0．06a0．53±0．08a GT-D1．45±0．11a0．10±0．02a0．11±0．02a0．79±0．11b

不同小写字母表示同列不同行在0.05水平上具有显著差异;F/B：真菌与细菌的比值,the ratio of fungi to bacterial; M/B：单一不饱和脂肪酸与支链脂肪酸的比值,the ratio of monounsaturated fatty acid to branched-chain fatty acid； SA/SI：支链反异构饱和脂肪酸与支链异构饱和脂肪酸的比值,the ratio of branched anti-isomeric saturated fatty acids to branched chain isomeric saturated fatty acids

不同地下水位梯度湿地土壤微生物群落的多样性分析表明, Shannon-Wiener多样性指数(H′)和Pielou 均匀度指数(J)的最低值出现在GT-D,而最高值出现在GT-B。Simpson优势度指数(D)则与之相反(图5；P< 0.05),表明地下水位梯度的升高降低了土壤微生物结构多样性,而促进了优势种群的分布。

图5　不同地下水位梯度土壤微生物群落Shannon-Wiener多样性指数(H′),Pielou 均匀度指数(J)和Simpson优势度指数(D)Fig.5　Shannon-Wiener diversity index, Pielou evenness index and Simpson dominance index of soil microbial communities under different ground water level gradient zones不同小写字母表示在0.05水平上具有显著差异

2.4　土壤微生物群落结构与土壤属性的关系

相关分析表明(表4),湿地土壤化学性质中,pH与G(+)、G(-)、细菌、放线菌、总PLFAs和优势度指数呈极显著正相关(P<0.01),而与多样性指数和优势度指数呈极显著负相关(P<0.01)。TOC与多样性指数和均匀度指数呈极显著正相关(P<0.01)而与优势度指数呈极显著负相关(P<0.01)。C∶N与G(-)、放线菌和优势度指数呈极显著负相关(P<0.01),而与多样性指数和均匀度指数呈极显著正相关(P<0.01)。湿地土壤物理结构中,容重粘粒和粉粒与微生物群落结构表现出相同的关系,即与真菌以外的微生物群落和优势度指数呈极显著负相关(P<0.01),而与多样性指数和均匀度指数呈极显著正相关(P<0.01),土壤沙粒含量则表现出相反的趋势(P<0.01)。MBC 与G(+)、G(-)、 细菌、放线菌和总PLFA均呈极显著正相关(P<0.01),与优势度指数呈显著正相关(P<0.05),而与多样性指数呈极显著负相关(P<0.01),MBN则与G(+)、G(-)、细菌、总PLFA和优势度指数呈极显著正相关(P<0.01),与放线菌呈显著正相关(P<0.05),与多样性指数呈极显著负相关(P< 0.01)。

表4　微生物群落结构与土壤环境因子间的相关性分析

** 表示在0.01水平上具有极显著相关性；*表示在0.05水平上具有显著相关性

土壤微生物商和土壤环境因子间具有显著的线性关系(图6),MBC/SOC 在研究区范围内随着pH和土壤含水量的增大而增大,随着容重和粘粒含量的增大而减小。但是与pH、容重和粘粒含量相比,MBC/SOC与含水量间的回归拟合效果明显较差(R2=0.20,P=0.03)。MBC/SOC对pH变化的响应最敏感(R2=0.72,P<0.001),其次为粘粒含量。

图6　土壤微生物商与土壤环境因子的关系Fig.6　Relationship between microbial quotient and soil environmental factors

对所有地下水位带湿地土壤微生物群落结构进行冗余分析发现,前2个排序轴累积贡献值为94.8%,其中第1排序轴解释了82.9%的变异(图7)。在所有土壤环境因子中,Monte Carlo 检验表明,土壤微生物群落结构与 pH(F=17.3,P=0.002),土壤含水量(F=3.8,P=0.04),和 C∶N(F=2.7,P=0.05)显著相关。pH和土壤含水量与第1轴呈显著正相关,而C∶N与第1轴呈显著负相关。通过单个PLFA与环境因子相关性排序可知,革兰氏阴性细菌PLFA标记(i15:0 3OH、i17:0 3OH、15:1ω6c、15:0、cy17:0)和革兰氏阳性细菌PLFA标记(i13:0、i14:0、a16:0)与土壤含水量呈显著正相关,即这些PLFAs的相对丰度值随着土壤含水量的升高而升高。真菌特征PLFA(18:3ω6c和18:1ω9c)与C∶N呈显著正相关。革兰氏阳性细菌PLFA标记(i15:0、14:0)、革兰氏阴性细菌PLFA标记(14:1ω5c)和真菌PLFA标记(16:1ω5c)与pH呈显著正相关。革兰氏阳性细菌PLFA标记(a17:0)、革兰氏阴性细菌PLFA标记(12:0)与沙粒含量呈显著正相关。放线菌PLFA标记(10Me18:0、10Me16:0)与各土壤环境因子相关性均不显著。

图7　土壤微生物优势群落与环境因子的RDA排序Fig.7　Redundancy analysis (RDA) of the microbial distribution patterns constrained by soil physicochemical properties along ground water level gradients

3　讨论

水文条件是湿地形成和发育的驱动因子。当水淹没洲滩土壤时,进入土壤的氧气减少,起分解作用的好氧微生物活动减弱,有机质分解速率低,土壤有机碳以积累为主[31]。而在湿地土壤干湿交替地条件下,会引起土壤收缩和膨胀,使土壤团聚体中有机质暴露,对土壤微生物代谢具有激发效应,促进有机质的分解[32]。本研究发现不同地下水位梯度间土壤理化性质差异明显,湿地土壤机械组成的分异主要是因为水体波动过程中的冲刷和筛选作用导致细粒物质被水流带走而使得比重较大的粗颗粒物积累下来[33]。此外,地下水位波动引起的干湿交替过程会引起土壤胶体的形态变化,破坏土壤团聚体结构。土壤pH, 颗粒组成和容重的变化使原来不能被分解的有机质因团聚体的崩溃而加速分解,所以土壤有机碳的含量显著降低[34]。张全军等[35]研究发现高盖度、长势旺盛的苔草群落有利于养分保留和有机碳的固持,由于靠近水体的梯度带苔草盖度、株高、生物量等生物性状显著低于其他梯度,故而不利于土壤有机碳的积累。土壤全磷含量主要受成土母质、气候环境条件影响,受植物群落影响不大[36],高地下水位梯度土壤磷含量显著降低主要是因为可溶性磷酸盐被淋溶分解[37],而研究表明可溶性有机碳中的芳香类物质能够保护土壤粘粒抵抗外界的侵蚀[38],进而又促进了粘粒含量的减少。

土壤微生物量可以反映土壤同化和矿化能力的大小,是土壤活性大小的标志[21]。本研究中靠近水体的高地下水位梯度具有最高的生物量碳氮,可能是因为pH和土壤含水量综合作用的结果。枯水季节期间,洲滩土壤含水量在30%—55%之间,土壤含水量的升高有利于微生物活性的增强。由于枯水期洲滩土壤pH在4.5—5.5之间,高地下水位梯度pH和微生物量碳氮均高于其他梯度,说明研究区土壤微生物更喜欢弱酸性的环境。

湿地土壤地下水位改变了土壤微生物群落组成。随着地下水位梯度空间位置的升高,PLFA总量、细菌群落、放线菌群落的PLFA含量显著增加,最高值均出现在地下水位梯度最高的GT-D。Mentzer等研究发现湿地水文情势的变化能够促进土壤微生物群落结构的持久性变化[39]。湿地水文波动是影响土壤氧气浓度的驱动型因素,但是在枯水季节,好氧环境和厌氧环境会根据降水或湖泊水位变化同时出现在高程梯度较低的区域[40],因此在高程低,靠近水体的地带,水位会在地面表层上下波动,进而导致干湿交替过程。本研究中的结果表明高地下水位区域频繁的干湿交替能够促进细菌和放线菌群落的丰度,与Foulquier等的研究结果一致[41]。真菌群落对酸性的耐受程度高于细菌[42],因此本研究中真菌群落并没有随pH产生明显变化。

由于土壤微生物量碳氮的周转快,所以它们的分配比例越大,土壤有机碳等的周转速率就越大[43]。靠近水体高地下水位梯度带土壤微生物商和微生物量氮分配比例最高,说明土壤有机碳氮周转最快,这也解释了土壤有机碳在此梯度带含量低的原因。微生物商与土壤含水量之间的相关性远远弱于pH、容重和粘粒含量。研究表明,短时期的淹水对土壤有机碳的周转速率影响不明显[43],说明地下水位梯度对微生物商的影响主要是通过土壤pH和土壤质地的作用。

相关分析和RDA分析都表明,在地下水位梯度空间分异过程中,pH、土壤含水量和土壤机械组成是影响土壤微生物群落结构的主要环境因素。因为土壤pH会对微生物吸附、胞外酶的分泌以及养分利用产生不同的作用,故而其对微生物群落结构的影响机制比较复杂。一般来说,细菌和放线菌适合生长在微碱性或弱酸性环境下,而真菌则适合生长在酸性条件下[44]。因此细菌和放线菌的PLFA含量与土壤pH呈极显著正相关,而随着pH的升高,真菌所占的比重逐渐降低。土壤机械组成对微生物群落结构的直接影响主要是因为土壤微生物主要存在于土壤颗粒的空隙或表面。研究发现,细菌主要存在于0.5—5μm的空隙中,而放线菌也适合生存在相对较大的空隙中[45],土壤容重随着土壤孔隙度的增大而减小,而本研究中细菌和放线菌的含量与土壤容重具有显著的负相关性,可能是因为在湿地环境中土壤颗粒空隙的增大有利于细菌和放线菌的存活。Cook等发现,土壤微生物群落中细菌在高水位条件下活性高而真菌则在低水位条件下活性较高[46]。无论是在陆地生态系统还是在湿地生态系统土壤水分都是土壤微生物群落结构的重要控制因素。湿地土壤地下水位的抬升,最直接的影响就是使土壤含水量升高。土壤碳氮比可以反映微生物生长是碳限制还是氮限制,同时也是影响土壤微生物群落结构的重要因子。碳氮比越小说明碳的可利用程度越低,由于真菌需要有机残体作为生命活动的载体,与外源有机碳输入和土壤有机质积累关系密切[47],所以本研究中细菌和放线菌群落对碳限制所带来的环境压力耐受性显著高于真菌群落,有机碳可利用程度的降低所带来的负面效应并没有影响pH升高所带来的正面效应。在湿地土壤环境下,地下水位的空间变化不仅关系到土壤水分条件和氧气浓度的变化,而且还反映了土壤结构和养分的可利用性,因此对土壤微生物群落结构影响复杂。土壤微生物群落的多样性指数和均匀度指数均随着地下水位梯度升高而降低,说明地下水位梯度所综合表征的土壤环境梯度的变化对土壤微生物具有明显的筛选作用,不同的湿地微生态环境塑造了不同的微生物群落,这种微生物群落结构的变化可能是湿地枯水季节地下水位短期分异的结果[21]。

真菌/细菌(F/B)常被用来衡量微生物群落对环境变化的响应程度,进而能够反映土壤微生物的生理状态,比如有机质的分解[47]。与其他地下水位梯度带相比,GT-D具有最低的F/B,说明高地下水位梯度碳的积累效率较慢[48],与微生物商所表现的结果一致。Gordon等[49]研究发现在干湿交替环境中F/B并未发生明显变化,间接指示pH可能是影响F/B变化的主要环境因子。单一不饱和脂肪酸/支链脂肪酸(M/B)可以用来指示厌氧与好氧微生物的比例[50]。本研究中GT-C和GT-D的M/B显著低于GT-A和GT-B,说明随着地下水位的升高,厌氧微生物的比重在逐渐增多,而好氧微生物比重的增长速度低于厌氧微生物,这与土壤含水量的变化是一致的。支链反异构饱和脂肪酸/支链异构饱和脂肪酸(SA/SI)能够反映微生物群落环境压力的变化,一般随着温度的降低和有机碳可利用程度的增加而升高[51-52]。由于pH的升高和土壤颗粒的变化促进了土壤有机质的分解和微生物的代谢强度,并且随着地下水位的升高,土壤温度降低,所以本研究中SA/SI 随地下水位呈显著增大趋势。

鄱阳湖湿地不同地下水位环境梯度间土壤微生物群落结构存在较大差异,由于本研究实验样地区域为单一灰化苔草群落,所以土壤微生物群落结构的分异主要是由于地下水位变化及其所驱动的土壤理化性质的改变所引起的。本研究则发现,在枯水期湿地洲滩出露阶段,鄱阳湖湿地短时期内所表现出的不同的地下水位环境梯度空间异质性既改变了微生物活性又同时显著改变了土壤微生物群落结构,进而有可能对有机质分解、土壤呼吸等湿地生态过程产生不同影响,同时靠近水体、地下水位相对较高的环境梯度带土壤碳氮周转速率快,是湿地生物地球化学循环的热点区域。

4　结论

(1)地下水位是湿地洲滩微地形、高程和水文情势综合作用的结果。地下水位的升高显著改变了土壤化学性质和物理结构。随着地下水位梯度带从最低到最高,土壤pH从4.62升至5.05,沙粒含量增加了4.28倍,而有机碳、粘粒和粉粒含量分别降低了53.5%、34.6% 和 7.0%。

(2)湿地地下水位改变驱动下土壤环境因子的变化进一步改变了湿地土壤微生物量碳氮和土壤碳氮周转速率。与地下水位最低的梯度带(GT-A)相比,当平均地下水位在地表上下波动时(GT-D),微生物量碳氮分别增长了2.82和4.30倍,微生物商提高了5.77倍,微生物量氮分配比例提高了7.15倍。

(3)湿地地下水位梯度驱动下湿地微环境的变化短期内改变了土壤微生物群落结构,与地下水位最低的梯度(GT-A)相比,干湿交替环境(GT-D)下,微生物总PLFAs,细菌PLFA和放线菌PLFA的含量分别增长了106.79%、117.24%、74.89%,随着地下水位升高,真菌PLFA比重逐渐降低而厌氧微生物比重增多。

(4)湿地地下水位梯度的升高增加了土壤微生物群落的环境压力,进而降低了土壤微生物群落多样性。高地下水位梯度(GT-D)具有最低的F/B(0.11)和最高的SA/SI(0.79)。不同地下水位环境梯度土壤微生物的多样性指数和均匀度指数表现为GT(A)> GT(B)> GT(C)> GT(D),而优势度指数则与之相反。

(5)在湿地地下水位变化过程中,影响土壤有机碳周转的主要土壤环境因子为pH和土壤质地,而影响土壤微生物群落结构的关键限制性因子为pH,土壤含水量,沙粒含量和C∶N。