The Macrobenthos Community Structure and Its Relationships with Water Quality Factors in the Liujiang River Basin
-
摘要: 以柳江流域23个点位为例，研究了不同水文季节(3月和6月)水质因子(DO、CODMn、NH3-N、As、Hg、Cd、Pb)对底栖动物群落结构的影响，2次调查共鉴定出底栖动物4门7纲16目54科68属(种)，种类相似系数为83.76%，物种丰富度表现为枯水期略高于丰水期，底栖动物优势种在枯、丰水期分别有4种和2种，具有明显的季节演替. DCA和TWINSPAN分析结果显示：柳江流域底栖动物群落结构在时空分布上存在一定差异，主要表现为枯、丰水期均有3个不同的底栖动物类群大致依次分布在该流域干流及支流的源头、中上游及下游，群落内指示物种呈现出由清洁种逐渐过渡到耐污种的变化. 指示物种分析结果(ISA)显示：在枯水期，共有8种底栖动物能作为柳江源头和大环江上游群落的指示种，分别为四节蜉属(Baetis sp.)、锯形蜉(Serratella sp.)、溪泥甲一属(Zaitzevia sp.)、短脉纹石娥属(Cheumatopsyche sp.)、蜗虫(Dugesia sp.)、朝大蚊属(Antocha)、摇蚊属(Chironomus sp.)和细蜉属(Caenis sp.)，其他类群无显著指示种；在丰水期，共有5种底栖动物能作为群落划分的指示种，其中，方格短沟蜷(Semisulcospira cancellata)、钉螺(Oncomelania sp.)是柳江中上游的指示种，汉森安春蜓(Amphigomphus hansoni)、河蚬(Corbicula fluminea)、卷扁螺(Gyraulus compressus)是柳江下游的指示种. CCA分析结果显示：CODMn和DO是枯、丰水期均对柳江流域底栖动物群落分布有重要影响的水质驱动因子，而重金属因子(Cd、Pb、Hg)在枯水期对底栖动物分布影响较大, NH3-N在丰水期对底栖动物分布影响较大.Abstract: In order to study the influence of water quality factors (DO, CODMn, NH3-N, As, Hg, Cd, Pb) on macrobenthos community structure in different hydrological seasons in the Liujiang River basin, investigations were conducted in the Liujiang River basin and 4 phyla, 7 classes, 16 orders, 54 families and 68 genera (species) of macrobenthos were identified. The Czekanowsiki is 83.76% and the species richness is slightly higher in the dry season than in the wet season. There were 4 and 2 dominant species of macrobenthos in dry and wet seasons respectively, showing obvious seasonal succession. The results of TWINSPAN and DCA analysis show that there are certain differences in the temporal and spatial distribution of the macrobenthos community structure in the Liujiang River basin. In both the dry and the wet season, there are 3 different macrobenthos communities distributed respectively around the source, the upper and middle reaches and the lower reaches of the main stream and tributaries of the river. The indicator species in the community show a gradual transition from clean species to pollution-tolerant species. The results of indicator species analysis shows that, in the dry season, there are a total of 8 macrobenthos species that can be used as indicator species for the communities at the source of the Liujiang River and the upper reaches of the Dahuan River, namely Baetis sp., Serratella sp., Zaitzevia sp., Cheumatopsyche sp., Dugesia sp., Antocha, Chironomus sp. and Caenis sp., and the other communities have no significant indicator species. In the wet season, there are 5 macrobenthos species that can be used as indicator species for community division, among which Semisulcospira cancellata and Oncomelania sp. are indicator species in the middle and upper reaches of the Liujiang River and Amphigomphus hansoni, Corbicula fluminea and Gyraulus compressus are indicator species in the lower reaches of the Liujiang River. The CCA analysis results show that CODMn and DO are the water quality factors that have important effects on the macrobenthos community structure in the Liujiang River basin in the dry and wet seasons, the heavy metal factors (Cd, Pb, Hg) have a greater impact on the distribution of macrobenthos in the dry season, and NH3-N has a greater impact on the distribution of macrobenthos in the wet season.
-
密度泛函理论(DFT)[1-4]从理论计算的角度为材料设计和微观机理分析提供了一种有效的方法。通过DFT计算可以获得物质的电子结构、光学性质、磁矩以及弹性常数等性质，这为理解材料以及加速材料设计提供了理论指导[5-7]。最近，随着半导体技术的飞速发展以及计算机工业的不断革新，半导体集成电路越做越小。但是当集成电路芯片接近原子尺度时，发热问题越来越突出。为了克服这一障碍，使微电子器件更加强大，材料上的突破是关键。目前半金属材料是研究热点，半金属铁磁体在自旋电子学方面具有非常大的潜在应用。半金属性是指材料在一个自旋方向上表现出非金属性而在另一个方向上表现出金属性。自旋电子学的发展解决了由传统的半导体电子技术组成设备的诸多缺点，如：数据处理速度慢、能耗高、存储容量小等[8]。目前已经研究出许多可用于自旋电子学的半金属材料，例如钙钛矿、Heusler合金、硫族化合物等[9-10]。Heusler合金因其具有高自旋极化率、高居里温度等优良的物理特性而得到广泛研究和关注[11]，而且Heusler合金表现出的高热电效率使其能有效地将热量转化成电能，有利于环保以及节能技术的开发。BAAZIZ等[6]通过DFT计算发现ZrCoFeX(X=Si，Ge)表现出半金属铁磁性，而且对这2种化合物的力学、光学、电磁特性也进行了深入的研究。ZOULIKHA等[12]对CoFeYGe(Y=Ti，Cr)的磁性、弹性常数、电磁以及热电性能进行DFT模拟后发现材料具有0.489 eV和0.640 eV的带隙宽度，表现出优良的热电性能。LIU等[5]发现了一系列满足半金属性质的新的四元Heuler化合物CrLaCoZ(Z=Al，Ga，In，Ge，Sn，Pb)，在这些化合物中除了CrLaCoSn外其他化合物的居里温度明显高于室温，有利于作为磁隧道结中的电极材料。
Heusler化合物由于其优良的物理性质而得到了广泛关注，例如各种奇异的磁现象、高居里温度和高自旋极化率等[13-16]。理论计算中正分的Heusler化合物按照元素种类可分为二元、三元和四元Heusler化合物，其中常见的三元Heusler化合物可进一步分为三元半Heusler化合物和三元全Heusler化合物，其化学式分别表示为XYZ和X2YZ(X和Y代表过渡金属原子，而Z代表主族原子)。当三元全Heusler化合物中的一个过渡金属原子被另一个不同的过渡金属原子取代时可得到化学式为XX′YZ的四元Heusler化合物，这种新型的四元Heusler化合物也称为LiMgPdSn型Heusler化合物。此类四元结构具有1∶1∶1∶1的化学计量比，相对于三元Heusler化合物而言具有可忽略不计的无序散射性质而表现出低功耗的特点[11]，所以四元LiMgPdSn型的新结构类型成为了新的研究热点。另外，已经证明Heusler化合物半金属的总磁矩与价电子总数之间满足Slater-Pauling规则[17-20]，而四元Heusler化合物半金属也满足这一规则，综合以上优势，四元Heusler化合物中具有100%自旋极化率的半金属铁磁体毫无疑问地成为了自旋电子学器件的理想候选材料[11, 21-22]。
根据Slater-Pauling规则：Heusler半金属的分子自旋磁矩(Mt)和价电子数(Zt)的差值等于18、24或28，本文前期研究已对满足Mt=Zt-28的Slater-Pauling规则的四元Heusler化合物进行了较为系统的研究[23-25]，但发现满足差值为18和24的四元Heusler化合物的研究较少，鲜有计算半金属铁磁体性质的研究[26]。本文对大量四元Heusler化合物进行了详细的DFT计算，并对其力学性质、态密度和能带进行了分析。
1. 研究方法
采用基于密度泛函理论的全势线性缀加平面波方法(FP-LAPW)进行计算[27-28]，使用的软件为WIEN2k代码包，所用的交换关联势为在合金计算中使用广泛的Perdew-Burke-Ernzerhof的广义梯度近似(GGA)。在具体计算中，截断参数
Rmt×Kmax=9, 其中，Rmt是最小的muffin-tin球半径，Kmax是指平面波扩展中的最大K矢量。在muffin-tin球内的球谐函数中扩展到lmax=10。四面体网格积分中选取5 000个k点，在自洽循环中选择能量和电荷作为收敛标准，收敛标准为默认值。
四元Heusler化合物的晶体结构由沿主对角线的4个独立坐标A(0, 0, 0)、B(1/4, 1/4, 1/4)、C(1/2, 1/2, 1/2)和D(3/4, 3/4, 3/4)组成[29]，可以将4个等效晶位视为4个面心立方子晶格，用晶胞对角线连接起来，空间群为F-43m(216)。A和B位原子选取Fe、Co、Ni、V、Cr、Mn等元素；C位原子为Ti、V、Cr、Mn等元素，C位原子不再选用符合Mt = Zt -28规则的Cu和Zn元素；D位原子选为Al、As或者Si。经验证此类四元化合物的铁磁态能量均低于非磁态能量，故默认计算材料的铁磁态结构。对于每种组合，先以能量和电荷作为收敛标准进行体积优化，再使用最低能量对应的晶格常数作为平衡晶格常数计算性质。
2. 结果与讨论
2.1 力学性质
通过计算弹性常数Cij探究物质的组成稳定性，根据Born-Huang标准[13]，立方晶系的力学稳定性标准如下：
C11>C12,C44>0,(C11+2C12)>0。 (1) 表 1列出了计算得到的四元Heusler化合物经过体积优化后的平衡晶格常数(a0)、弹性常数(Cij)、剪切模量(G)、杨氏模量(E)、体积模量(B)、泊松比(ν)、剪切各向异性常数(A)和Pugh比(B/G)[30-32]，计算公式如下：
B=C11+2C123, (2) G=C11−C12+3C445, (3) E=9BG3B+G, (4) ν=3B−2G2(3B+G), (5) A=2C44C11−C12。 (6) 表 1 全有全无分配法结果Table 1. The result of all-or-nothing assignment methodID Material C11 C12 C44 ν B/G A a0 B E G 1 FeVTiAs 176.4 157.2 20.4 0.452 10.169 2.1 5.96 163.4 46.7 16.1 2 CoCrTiSi 314.1 148.6 58.4 0.350 2.990 0.7 5.80 203.8 183.9 68.1 3 FeCrTiSi 257.7 166.7 93.3 0.333 2.657 2.1 5.81 197 197.7 74.2 4 CoMnVAs 285.5 168.6 89.0 0.335 2.702 1.5 5.77 207.4 205 76.8 5 FeMnTiAl 278.9 136.7 82.8 0.314 2.355 1.2 5.89 184 205.3 78.1 6 FeVTiAl 192.3 119.7 114.0 0.258 1.734 3.1 6.04 143.8 208.7 82.9 7 FeCrTiAs 280.4 140.1 93.0 0.305 2.226 1.3 5.86 186.7 218.8 83.9 8 MnCrTiAs 286.6 132.8 96.4 0.292 2.073 1.3 5.90 183.6 228.9 88.6 9 FeCrVAl 264.5 145.6 109.1 0.292 2.073 1.8 5.84 185.1 230.7 89.3 10 MnCrVAs 319.0 164.8 101.0 0.315 2.364 1.3 5.79 216.1 240.3 91.4 11 FeCrVSi 334.5 176.1 129.5 0.294 2.092 1.6 5.70 228.8 283 109.4 12 FeMnVAl 369.2 127.5 117.4 0.260 1.750 0.97 5.76 207.9 299.3 118.8 13 CoFeTiSi 285.6 178.1 169.0 0.259 1.740 3.1 5.73 213.9 309.4 122.9 14 MnCrVSi 393.2 119.0 114.3 0.255 1.704 0.8 5.75 210.3 309.6 123.4 15 CoVTiAs 219.1 130.4 203.3 0.161 1.144 4.6 5.97 159.8 324.6 139.7 16 CoFeCrSi 284.1 212.1 194.0 0.266 1.804 5.4 5.62 236 331.1 130.8 17 CoCrVAl 261.2 152.6 197.3 0.202 1.347 3.6 5.82 188.7 337 140.1 18 CoFeMnSi 391.2 152.7 151.4 0.251 1.675 1.3 5.61 232.1 346.6 138.5 19 CoFeMnAs 344.9 121.2 177.2 0.193 1.295 1.6 5.74 195.6 360.4 151.1 20 FeMnTiSi 402.5 125.4 166.7 0.212 1.401 1.2 5.74 217.8 376.8 155.5 21 CoFeVAl 269.4 184.4 283.1 0.16 1.138 6.7 5.73 212.7 433.6 186.9 22 CoFeVAs 682.8 -49.7 78.4 0.126 1.004 0.2 5.78 194.3 435.9 193.6 23 CoFeVSi 525.2 81.5 169.9 0.174 1.202 0.8 5.66 229.3 447.9 190.7 24 CoMnTiAs 193.0 192.6 364.5 0.088 0.880 1 687.5 5.83 192.6 476.1 218.8 25 CoVTiSi 213.9 154.0 473.5 -0.043 0.587 15.8 5.90 173.8 566.5 296.1 26 FeMnVAs 1 990.8 -602.8 89.9 -0.133 0.457 0.1 5.74 261.8 993.6 572.7 27 CoMnCrAl 2 880.1 -1 175.1 3 000.9 -0.75 0.067 1.5 5.72 174.3 1 307.1 2 611.6 28 CoMnTiAl 415.1 134.6 6 562.5 -0.781 0.057 46.8 5.86 227.9 1 751.3 3 993.6 29 CrVTiSi -49.6 245.7 17.5 0.687 -3.002 -0.1 6.00 145.8 -163.9 -48.6 30 CrVTiAs -400.6 449.6 -408.6 8.002 -0.400 1 6.03 166.1 -7 475.7 -415.2 31 MnVTiSi 177.7 210.8 65.9 0.422 6.049 -4 5.93 199.3 93.7 32.9 32 MnCrTiAl 234.1 100.1 -32.8 0.476 20.303 -0.5 6.04 144.7 21 7.1 33 CoCrTiAl 95.2 182.6 110.4 0.356 3.142 -2.5 5.95 153.2 132.2 48.7 34 FeMnCrSi 618.9 101.3 -44.3 0.372 3.561 -0.2 5.59 273.9 211 76.9 35 CoMnVSi 639.8 21.1 -1 921.4 -3.934 -0.221 -6.2 5.66 227 6 039.1 -1 029.1 36 FeMnCrAs -20 332.5 1 069.8 -20.4 -1.489 -0.082 0.002 5.73 351.8 4 198.5 -4 292.7 37 CoFeMnAl -6.8 395.0 126.6 0.508 -59.45 -0.6 5.69 261.1 -13.2 -4.4 38 CoMnCrAs -230.2 377.5 186.6 0.528 -18.26 -0.6 5.75 174.7 -29.2 -9.6 39 CoFeCrAs 315.7 142.2 -3 903.1 -1.527 -0.087 -45 5.75 199.8 2 429.5 -2307.2 W(L) 28(11) 28(11) 28(11) 24(15) 15(24) 38(1) — — — — 注：加粗的数据表明化合物满足相应的性质，如力学稳定性、金属键贡献大、韧性和剪切各向异性。W(L)中，W和L分别表示满足、不满足上述性质的化合物个数。 通过Born-Huang标准发现共有28种化合物满足力学稳定性。Pugh[33]给出了判断材料脆性和韧性的简单规则：通过计算B/G的值来分析材料的脆性和韧性，若B/G > 1.75，材料表现出韧性，否则材料表现出脆性。通过比较发现有15种化合物表现出韧性，其他化合物则表现为脆性。泊松比ν>0.25表示金属键对化合物成分的贡献较大；否则，共价键的贡献较大。有24种化合物的金属键对化合物成分的贡献较大。对于各向同性的晶体，A=1[34]，很明显，表中除了CrVTiAs为各向同性晶体外，其他化合物均为各向异性，而且满足力学稳定性的28种化合物均为各向异性。为了验证数据的可靠性，图 1为FeVTiAl、CoCrTiAl、FeCrVAl、CoMnTiAs、FeCrTiSi、CoFeMnAs、CoCrTiSi、FeCrVSi以及CoFeMnSi这9种化合物在晶格常数以及体积模量上与现有研究数据[35-37]的对比。在晶格常数上的最大差值体现在化合物FeVTiAl上(Δa0=0.004 nm)，而在体积模量的最大差值体现在化合物CoCrTiAl上(ΔB=13.4)。对比结果表明：本文的结果与现有研究结果相近(图 1B)，说明本文所用方法可靠。
![]()
图 1 不同化合物的晶格常数、体积模量比较Figure 1. Comparison of different compounds in terms of lattice constant and bulk modulus力学稳定性对于Heusler化合物来说具有重要作用，所以进一步通过图 2对表 1中满足力学稳定性的28种化合物进行分析。杨氏模量表示化合物抵抗弹性形变的能力，而剪切模量表示化合物抵抗切应变的能力，两者均反映化合物的刚度性质。
![]()
图 2 杨氏模量、剪切模量的比较注：□表示化合物满足韧性且金属键贡献大。Figure 2. Comparison of Young's modulus and shear modulus从图 2可以看出，绝大多数化合物的杨氏模量大于剪切模量，表明化合物抵抗弹性形变的能力强于抵抗切应变的能力。继续对这些满足力学稳定性的化合物的韧性以及金属性进行分析，金属性越强表明材料的熔点以及硬度越大，而韧性又在一定程度上反应材料抗断裂的能力。对同时满足韧性以及金属键贡献大的化合物进行标注，发现此类化合物主要分布在低E和G值附近，共有11种化合物，即图 2中□对应的化合物。
2.2 电磁性质
根据Slater-Pauling规则，可以知道半金属化合物的总自旋磁矩是一个整数。表 2列出了总磁矩为整数的化合物的各原子磁矩和间隙磁矩、总自旋磁矩、半金属带隙和价电子数。FeVTiAl和MnVTiSi遵循Mt= Zt-18的关系，而CoMnCrAl、CoFeTiSi、CoMnTiAs、CoFeCrSi、CoFeMnSi、CoFeMnAs、FeMnCrAs、CoFeMnAl、CoMnCrAs以及CoFeCrAs遵循Mt=Zt-24的关系。
表 2 符合Slater-Pauling标准的化合物的电磁性质Table 2. Electromagnetic properties of compounds meeting Slater-Pauling criteriaID 化合物 Mi MX M′X MY MZ Mt GAP Zt 1 CoFeMnSi 0.01 0.90 0.59 2.54 -0.04 4 0.695 28 2 CoFeMnAs 0.11 1.04 1.02 2.81 0.01 5 0.224 29 3 CoFeMnAl -0.03 0.82 -0.10 2.35 -0.03 3 0.632 27 4 CoMnCrAs 0.20 1.14 -0.34 1.98 0.03 3 0.477 27 5 CoMnCrAl 0.02 0.70 -0.98 1.28 -0.02 1 0.285 25 6 CoFeCrSi 0.06 1.09 0.36 1.53 -0.03 3 0.665 27 7 CoFeCrAs 0.19 1.01 0.82 1.97 0.01 4 0.514 28 8 FeVTiAl 0.63 -0.49 1.61 0.25 0.000 3 2 0.601 20 9 MnVTiSi 0.52 -0.09 1.45 0.13 -0.004 2 0.457 20 10 FeMnCrAs 0.22 0.43 -0.58 1.91 0.02 2 0.528 26 11 CoFeTiSi -0.11 0.61 0.61 -0.10 -0.01 1 0.676 25 12 CoMnTiAs -0.07 0.67 0.46 -0.07 0.004 1 0.392 25 注：原子间隙磁矩(Mi，单位为μB)和各原子磁矩(MX、M′X、MY、MZ，单位为μB)；总自旋磁矩(Mt，单位为μB)；半金属带隙GAP(eV)和价电子数(Zt)；加粗的数据表明在局部磁矩中其对总磁矩的贡献最大。 在原子磁矩中所有化合物的主族原子的磁矩(MZ)和间隙磁矩(Mi)对总磁矩的贡献微乎其微。一方面，总磁矩主要由过渡金属原子的磁矩(MX，M′X，MY)贡献而来，而它们的贡献度与化合物的类型有关，但总的来说MY对总磁矩的贡献最大。另一方面，在这12个化合物中CoFeMnSi拥有最大的半金属带隙(0.695 eV)；而CoFeMnAs的带隙最小(0.224 eV)。由于这2种化合物只是主族原子不同，说明虽然主族原子对总磁矩的贡献可以忽略，但它能够有效的调整化合物带隙的大小，相似的结论也能从两组化合物中看出：(1)CoMnCrAl和CoMnCrAs；(2)CoFeCrSi和CoFeCrAs。既然主族元素的变化导致了这样的结果，不同主族元素的电子排布影响了半金属带隙，因为在费米能级附近的态密度变化的影响，甚至导致直接影响半金属性，对于半金属带隙出现在自旋向下方向的化合物，主族元素变化影响了自旋向下费米能级处的态密度之后就会出现这样的现象。
为了进一步研究这些磁矩为整数的化合物，计算得到自旋向上和自旋向下的总态密度(TDOS，图 3)。这12种化合物在平衡晶格常数下都显示出半金属的态密度特性。费米能级为0 eV，带隙出现在自旋向下的方向，由自旋向上和自旋向下态密度的不对称可知，这12个化合物都具有磁性，磁性强弱也可由表 2的总磁矩数值大小体现。从图 3可以看出，当主族原子改变时，即使化合物其他过渡金属原子保持不变，其态密度的形状也会发生巨大改变，例如：(1)CoFeMnSi、CoFeMnAs和CoFeMnAl；(2)CoMnCrAs和CoMnCrAl；(3)CoFeCrSi和CoFeCrAs。由此正好印证了前述对于半金属带隙能否出现的分析，所以总态密度的变化正是由于主族元素的分态密度(PDOS)变化所致。由于研究的样品数量较多，限于篇幅，本文未给出各化合物中每个元素的分态密度，依据前期研究[24]所给出的分态密度也能证明。另外，当化合物使用其他基底的晶格常数时，其态密度会发生移动，但总体形态保持不变[23]，所以在此类薄膜样品的实验制备过程中，基底材料的晶向带来的轻微应力可以起到调制费米能级的作用。
自旋极化率使用费米能级处的总态密度计算[38]：
p=ρ↑(Ef)−ρ↓(Ef)ρ↑(Ef)+ρ↓(Ef)×100%。 (7) 通过自旋极化率的计算可以得到这12种化合物的自旋极化率均为100%。图 4为这12种化合物对应的能带结构图，从自旋向下的实线可以直观读出半金属带隙的大小。CoFeCrAs、CoMnCrAs、CoMnTiAs和CoFeMnAs具有直接间隙，带隙分别为0.514、0.477、0.392、0.224 eV，并且都位于Γ点。而MnVTiSi、CoFeCrSi、FeMnCrAs、CoFeMnSi、CoFeTiSi、CoMnCrAl、FeVTiAl以及CoFeMnAl为间接带隙，其中CoFeMnSi的半金属间隙最大(约0.695 eV)。通过分析能带图发现，具有直接带隙的化合物均含有Co和As原子。当缺少其中一个原子时，化合物则表现出间接带隙。例如不含Co原子的FeMnCrAs以及不含As原子的CoFeCrSi、CoFeMnSi、CoFeTiSi、CoMnCrAl和CoFeMnAl。在本研究体系中，这是一个非常有意思的现象，可见这里同时含有Co和As原子的四元半金属化合物自旋向下方向的电子填充在Γ点处导带底的能量最低，同时在Γ点处的价带顶的能量最高。进一步根据原子的投影态密度以及能带图的原子特征分析可知，对于4个含有Co和As的直接半金属带隙化合物，其自旋向下方向能带的导带底主要由Co和As的s电子提供形成类s能带，价带顶是由Co原子的d电子提供的，而且都出现在了Γ点。在费米能级附近过渡金属Co原子s和d轨道电子和As的s轨道电子对于半金属形成直接带隙起到了至关重要的作用。
![]()
图 4 12种半金属化合物在自旋向上(虚线)和自旋向下(实线)方向的能带图Figure 4. The band of 12 half-metallic compounds in the spin-up (dashed line) and spin-down (solid line) directions3. 结论
对四元Heusler化合物进行第一性原理计算，通过弹性常数分析表明有28种四元Heusler化合物表现出力学稳定性，通过态密度图和能带图发现有12种化合物表现出半金属性，可作为自旋电子学的候选材料。CoFeCrAs、CoMnCrAs、CoMnTiAs和CoFeMnAs具有直接间隙，带隙分别为0.514 eV、0.477 eV、0.392 eV和0.224 eV。而MnVTiSi、CoFeCrSi、FeMnCrAs、CoFeMnSi、CoFeTiSi、CoMnCrAl、FeVTiAl和CoFeMnAl为间接带隙。化合物CoFeMnSi的半金属带隙最大，约为0.695 eV。磁矩分析表明12种化合物均满足Slater-Pauling规则。
-
![]()
图 2 柳江各样点底栖动物个体总数分布图
Figure 2. The amounts of macrobenthos in the Liujiang River basin
![]()
图 3 柳江各样点底栖动物物种丰富度分布图
Figure 3. The species richness of macrobenthos in the Liujiang River basin
![]()
图 4 柳江流域大型底栖动物采样点在DCA轴上的分布
Figure 4. The distribution of sampling sites of Liujiang River basin in the DCA
![]()
图 5 柳江流域环境因子与物种组成的CCA排序图(枯水期)
Figure 5. The CCA biplot of environmental variables and species composition in the Liujiang River basin (dry season)
1: Cheumatopsyche sp.; 2: Macromidae; 3: Antocha; 4: Simulium sp.; 5: Chironomus sp.; 6: Epeorus sp.; 7: Choroterpes sp.; 8: Baetis sp.; 9: Afronurus rubromaculata; 10: Labiobaetis ancoralis; 11: Labiobaetis atrebatinus orientalis ; 12: Procloeon sp.; 13: Caenis sp.; 14: Serratella sp.; 15: Micronecta; 16: Zaitzevia sp.; 17: Protohermes sp.; 18: Sulcospira cf hainanensis; 19: Sinotaia sp.; 20: Stenothyra sp.; 21: Radix sp.; 22: Semisulcospira sp.; 23: Tarebia granifera ; 24: Physa sp.; 25: Galba sp.; 26: Limnoperna lacustris; 27: Macrobrachium sp.; 28: caridina sp.; 29: Corbicula fluminea; 30: Glossiphonia sp.; 31: Dugesia sp.
![]()
图 6 柳江流域环境因子与物种组成的CCA排序图(丰水期)
Figure 6. The CCA biplot of environmental variables and species composition in the Liujiang River basin (wet season)
1: Cheumatopsyche sp2; 2: Amphigomphus hansoni; 3: Leptogomphus intermedius; 4: Baetis sp.; 5: Labiobaetis ancoralis; 6: Labiobaetis atrebatinus orientalis; 7: Habrophlebiodes gilliesi; 8: Ordobrevia sp.; 9: Bezzia sp.; 10: Culex sp.; 11: Chironomussp.; 12: Micronecta sp.; 13: Micronecta quadristrigata; 14: Neocaridina serrata; 15: Corbicula fluminea; 16: Galba pervia; 17: Gyraulus compressus; 18: Sinotaia sp.; 19: Bellamya sp.; 20: Physa acuta; 21: Semisulcospira cancellata; 22: Melanodies sp.; 23: Tarebia granifera; 24: Oncomelania sp.; 25: Radix auricularia; 26: Radix swinhoei.
表 1 柳江流域采样点地理位置分布表
Table 1 The locations of sampling sites in the Liujiang River basin
采样点 河流 干/支流 经度 纬度 采样点 河流 干/支流 经度 纬度 S1 都柳江 干流 107°35.29′ 25°51.49′ S13 大环江 支流 108°18.14′ 24°53.46′ S2 都柳江 干流 108°30.46′ 25°55.13′ S14 大环江 支流 108°09.10′ 24°43.34′ S3 都柳江 干流 108°54.16′ 25°44.37′ S15 大环江 支流 108°08.49′ 24°42.51′ S4 都柳江 干流 109°24.23′ 25°42.21′ S16 小环江 支流 108°28.04′ 24°34.24′ S5 融江 干流 109°23.51′ 25°13.12′ S17 龙江 支流 108°28.12′ 24°34.04′ S6 融江 干流 109°15.12′ 25°03.15′ S18 龙江 支流 108°39.13′ 24°29.60′ S7 打狗河 支流 107°46.57′ 25°34.49′ S19 柳江 干流 109°10.23′ 24°40.15′ S8 打狗河 支流 107°53.22′ 25°25.04′ S20 柳江 干流 109°15.03′ 24°31.45′ S9 龙江 支流 107°51.33′ 24°50.05′ S21 柳江 干流 109°23.53′ 24°19.47′ S10 龙江 支流 107°55.02′ 24°44.56′ S22 柳江 干流 109°40.19′ 23°57.21′ S11 大环江 支流 108°15.19′ 25°23.26′ S23 柳江 干流 109°31.52′ 23°50.30′ S12 大环江 支流 108°12.21′ 50°00.54′ 
下载: 导出CSV
表 2 柳江流域大型底栖无脊椎动物优势种组成(Y≥0.02)
Table 2 The composition of dominant macrobenthos species in the Liujiang River basin(Y≥0.02)
时期 种类 占总个数的比例Pi/% 出现频率fi/% 优势度Y 枯水期 萝卜螺属Radix 7.64 69.57 0.053 摇蚊属Chironomus 5.13 43.48 0.022 小划蝽属Micronecta 16.80 39.13 0.066 囊螺属Physa 6.60 39.13 0.026 丰水期 新米虾属Neocaridina 10.15 47.83 0.049 小划蝽属Micronecta 14.75 21.74 0.032
下载: 导出CSV
表 3 柳江流域大型底栖动物主要群落分布类型
Table 3 The distribution types of main macrobenthic communities in the Liujiang River basin
时间 类群 采样点位 枯水期 A1 S7、S8、S12、S13、S14、S15、S19、S22 A2 S2、S3、S4、S5、S6、S9、S10、S16、S17、S18、S20、S21、S23 A3 S1、S11 丰水期 B1 S7、S11、S12、S13、S14、S19、S21、S22、S23 B2 S2、S3、S4、S5、S6、S16、S17 B3 S1、S8、S15
下载: 导出CSV
表 4 柳江大型底栖动物群落划分及指示物种分析
Table 4 The community division of macrobenthos in the Liujiang River basin and the indicator species analysis
时期 类群 指示种 指示值 P 枯水期 A3 Baetis sp. 1.000 0.005** Serratella sp. 1.000 0.005** Zaitzeviasp. 1.000 0.005** Cheumatopsyche sp. 0.999 0.005** Dugesia sp. 0.967 0.005** Antocha 0.994 0.001** Chironomus sp. 0.918 0.019* Caenis sp. 0.705 0.048* 丰水期 B2 Semisulcospira cancellata 0.845 0.010** Oncomelania sp. 0.756 0.032* B3 Amphigomphus hansoni 0.816 0.018* Corbicula fluminea 0.816 0.012* Gyraulus compressus 0.816 0.018* 注: * *表示P＜0.01，*表示P＜0.05.
下载: 导出CSV
表 5 柳江流域水质理化因子特征(平均值±标准差)
Table 5 The characteristics of water quality factors in the Liujiang River basin(Mean±SD)
水质因子的质量浓度 枯水期 丰水期 P A1 A2 A3 总 B1 B2 B3 总 ρ(DO)/(mg·L-1) 7.67±1.58 8.29±0.89 9.21±0.21 8.15±1.19 10.37±0.85 10.46±0.80 11.71±1.19 10.61±0.96 0.000** ρ(CODMn)/(mg·L-1) 1.65±0.59 1.37±0.42 1.05±0.21 1.44±0.49 1.96±0.42 1.51±0.29 1.90±0.44 1.78±0.42 0.003** ρ(NH3-N)/(mg·L-1) 0.15±0.14 0.27±0.62 0.02±0.02 0.21±0.47 0.29±0.21 0.18±0.14 0.24±0.20 0.24±0.18 0.009** ρ(As)/(μg·L-1) 8.630±8.314 6.460±7.881 11.000±4.243 7.609±7.674 1.956±1.528 1.437±0.506 1.057±0.501 1.623±1.128 0.228 ρ(Hg)/(μg·L-1) 0.025±0.055 0.007±0.017 — 0.013±0.035 0.032±0.023 0.033±0.024 0.020±0.020 0.031±0.022 0.000** ρ(Cd)/(μg·L-1) 0.109±0.308 0.02±0.038 — 0.049±0.181 0.314±0.323 0.181±0.077 0.207±0.139 0.248±0.234 0.000** ρ(Pb)/(μg·L-1) 2.611±2.639 1.509±1.323 1.26±0.438 1.870±1.868 15.732±36.913 3.126±0.935 6.060±4.031 9.561±25.395 0.001** 注: * *表示P＜0.01，*表示P＜0.05.
下载: 导出CSV
表 6 CCA排序轴对应的特征值、物种-水质因子相关性及累积变化率
Table 6 The Eigenvalues, species-environment correlation and changing percentages for CCA axis
时期 排序轴 特征值 物种-水质因子相关性 物种数据累积变化率/% 物种-水质因子关系累积变化率/% 枯水期 1 0.579 0.894 8.5 23.6 2 0.549 0.901 16.5 45.9 3 0.411 0.854 22.5 62.7 4 0.370 0.891 28.0 77.7 丰水期 1 0.582 0.909 11.8 32.7 2 0.408 0.851 20.1 55.6 3 0.310 0.867 26.4 73.0 4 0.216 0.805 30.8 85.1
下载: 导出CSV
表 7 CCA前2轴与水质因子相关系数
Table 7 The correlations between water quality factors and the first two axes of CCA
水质因子 枯水期 丰水期 轴1 轴2 轴1 轴2 DO -0.813 8 -0.032 9 -0.520 9 -0.607 0 CODMn 0.104 4 0.852 0 -0.818 1 0.013 7 NH3-N 0.233 6 0.097 5 0.160 7 0.346 2 As 0.271 0 -0.054 8 0.259 9 0.217 4 Hg 0.057 7 0.433 3 — — Cd 0.046 2 0.572 0 0.115 3 0.287 9 Pb 0.246 0 0.461 6 -0.008 7 0.007 7
下载: 导出CSV
-
[1] 戴纪翠, 倪晋仁. 底栖动物在水生生态系统健康评价中的作用分析[J]. 生态环境, 2008, 17(5): 2107-2111. doi: 10.3969/j.issn.1674-5906.2008.05.075 DAI J C, NI J R. Roles of benthos in the aquatic ecosystem health assessment, China[J]. Ecology and Environment, 2008, 17(5): 2107-2111. doi: 10.3969/j.issn.1674-5906.2008.05.075
[2] BENKE A C. Concepts and patterns of invertebrate production in running waters[J]. Sil Proceedings, 1993, 25(1): 15-38.
[3] 段学花, 王兆印, 徐梦珍. 底栖动物与河流生态评价[M]. 北京: 清华大学出版社, 2010. [4] 迟国梁, 赵颖, 官昭瑛, 等. 广东横石水河大型底栖动物群落与环境因子的关系[J]. 生态学报, 2010, 30(11): 2836-2845. https://www.cnki.com.cn/Article/CJFDTOTAL-STXB201011005.htm CHI G L, ZHAO Y, GUAN Z Y, et al. Relationship between macroinvertebrate and environmental variables in Hengshishui River, Guangdong, China[J]. Acta Ecologica Sinica, 2010, 30(11): 2836-2845. https://www.cnki.com.cn/Article/CJFDTOTAL-STXB201011005.htm
[5] 刘祥, 陈凯, 陈求稳, 等. 淮河流域典型河流夏秋季底栖动物群落特征及其与环境因子的关系[J]. 环境科学学报, 2016, 36(6): 1928-1938. https://www.cnki.com.cn/Article/CJFDTOTAL-HJXX201606005.htm LIU X, CHEN K, CHEN Q W, et al. The community structure of macroinvertebrate and its relationship to the environmental factors in summer andautumn within typical reaches of Huai River Basin[J]. Acta Scientiae Circumstantiae, 2016, 36(6): 1928-1938. https://www.cnki.com.cn/Article/CJFDTOTAL-HJXX201606005.htm
[6] 吴召仕, 蔡永久, 陈宇炜, 等. 太湖流域主要河流大型底栖动物群落结构及水质生物学评价[J]. 湖泊科学, 2011, 23(5): 686-694. https://www.cnki.com.cn/Article/CJFDTOTAL-FLKX201105005.htm WU S S, CAI Y J, CHEN Y W, et al. Assemblage structure investigation of macrozoobenthos and water quality bioassessmentof the main river systems in Taihu Basin[J]. Lake Sciences, 2011, 23(5): 686-694. https://www.cnki.com.cn/Article/CJFDTOTAL-FLKX201105005.htm
[7] 曹然, 黎征武, 毛建忠, 等. 北江大型底栖无脊椎动物群落结构及水质的生物评价[J]. 水资源保护, 2017, 33(4): 80-87. https://www.cnki.com.cn/Article/CJFDTOTAL-SZYB201704013.htm CAO R, LI Z W, MAO J Z, et al. Benthic macroinvertebrate community structure andbioassessment of water quality of Beijiang River[J]. Water Resources Protection, 2017, 33(4): 80-87. https://www.cnki.com.cn/Article/CJFDTOTAL-SZYB201704013.htm
[8] 渠晓东, 曹明, 邵美玲, 等. 雅砻江(锦屏段)及其主要支流的大型底栖动物[J]. 应用生态学报, 2007(1): 158-162. doi: 10.3321/j.issn:1001-9332.2007.01.027 QU X D, CAO M, SHAO M L, et al. Macrobenthos in Jinping reach of Yalong River and its main tributaries[J]. Chinese Journal of Applied Ecology, 2007(1): 158-162. doi: 10.3321/j.issn:1001-9332.2007.01.027
[9] DUFRENE M, LEGENDRE P. Species assemblages and indicator species: the need for a flexible asymmetrical approach[J]. Ecological Monographs, 1997, 67(3): 345-366. http://orbi.ulg.ac.be/bitstream/2268/145383/1/Dufrene Legendre 1997 Ecol Monographs.pdf
[10] 张清华, 韦永著, 曹建华, 等. 柳江流域饮用水源地重金属污染与健康风险评价[J]. 环境科学, 2018, 39(4): 1598-1607. https://www.cnki.com.cn/Article/CJFDTOTAL-HJKZ201804019.htm ZHANG Q H, WEI Y Z, CAO J H, et al. Heavy metal pollution of the drinking water sources in the Liujiang River basin, and related health risk assessments[J]. Environmental Science, 2018, 39(4): 1598-1607. https://www.cnki.com.cn/Article/CJFDTOTAL-HJKZ201804019.htm
[11] 孟伟, 张远, 渠晓东. 河流生态调查技术方法[M]. 北京: 科学出版社, 2011. [12] RICHOUX P, BOURNAUD M, USSEGLIO-POLATERA P, et al. 淡水无脊椎动物系统分类、生物及生态学[M]. 刘威, 王旭涛, 黄少峰, 译. 北京: 中国水利水电出版社, 2015. [13] 国家环保局本书编委会. 水和废水监测分析方法[M]. 北京: 中国环境科学出版社, 1989. [14] 李佳静, 刘威, 邓培雁. 柳江流域不同水文季节水质对附生硅藻群落的影响研究[J]. 生态科学, 2020, 39(1): 10-19. https://www.cnki.com.cn/Article/CJFDTOTAL-STKX202001002.htm LI J J, LIU W, DENG P Y. The effects of water quality on epilithic diatom community in different hydrological seasons in Liujiang River basin[J]. Ecological Science, 2020, 39(1): 10-19. https://www.cnki.com.cn/Article/CJFDTOTAL-STKX202001002.htm
[15] 王备新. 大型底栖无脊椎动物水质生物评价研究[D]. 南京: 南京农业大学, 2003. [16] 殷旭旺, 徐宗学, 高欣, 等. 渭河流域大型底栖动物群落结构及其与环境因子的关系[J]. 应用生态学报, 2013, 24(1): 218-226. https://www.cnki.com.cn/Article/CJFDTOTAL-YYSB201301035.htm YIN X W, XU Z X, GAO X, et al. Macrobenthos community structure and its relationships with environmental factors in WeiheRiver basin, Northwest China[J]. Chinese Journal of Applied Ecology, 2013, 24(1): 218-226. https://www.cnki.com.cn/Article/CJFDTOTAL-YYSB201301035.htm
[17] 李巧. 指示值方法及其在昆虫中的应用[J]. 应用昆虫学报, 2011, 48(2): 457-462. https://www.cnki.com.cn/Article/CJFDTOTAL-KCZS201102042.htm LI Q. Indicator Value (IndVal) method and its application[J]. Chinese Journal of Applied Entomology, 2011, 48(2): 457-462. https://www.cnki.com.cn/Article/CJFDTOTAL-KCZS201102042.htm
[18] 程谦伟, 孟陆丽, 易弋, 等. 柳江市区段底栖动物的生态研究[J]. 安徽农业科学, 2010, 38(17): 9048-9050. doi: 10.3969/j.issn.0517-6611.2010.17.085 CHENG Q W, MENG L L, YI Y, et al. Ecosystem survey of zoobenthos in the cantonal of Liujiang[J]. Anhui Agricultural Sciences, 2010, 38(17): 9048-9050. doi: 10.3969/j.issn.0517-6611.2010.17.085
[19] BUSS D F, BAPTISTA D F, SILVEIRA M P, et al. Influence of water chemistry and environmental degradation on macroinvertebrate assemblages in a river basin in south-east Brazil[J]. Hydrobiologia, 2002, 481(1/2/3): 125-136. http://www.onacademic.com/detail/journal_1000034817910310_a52a.html
[20] 姚凯胜, 和雅静, 赵永晶, 等. 流溪河大型底栖动物群落结构及其影响因素[J]. 中国环境监测, 2020, 36(3): 83-93. https://www.cnki.com.cn/Article/CJFDTOTAL-IAOB202003013.htm YAO K S, HE Y J, ZHAO Y J, et al. Macrozoobenthos community structure and environmental determinants of the Liuxi River[J]. Environmental Monitoring in China, 2020, 36(3): 83-93. https://www.cnki.com.cn/Article/CJFDTOTAL-IAOB202003013.htm
[21] GRAA M A S, PINTO P, CORTES R, et al. Factors affecting macroinvertebrate richness and diversity in portuguese streams: a two-scale analysis[J]. International Review of Hydrobiology, 2004, 89(2): 151-164. doi: 10.1002/iroh.200310705
[22] SLOANE P I W, NORRIS R H. Relationships of AUSRIVAS-Basedmacroinvertebrate predictive model outputs to a metal pollution gradient[J]. Journal of the North American Benthological Society, 2003, 22(3): 457-471. doi: 10.2307/1468274
[23] 周玉春, 支崇远, 王路, 等. 柳江水体重金属污染健康风险评价[J]. 安徽农业科学, 2013, 41(10): 4370-4372, 4375. doi: 10.3969/j.issn.0517-6611.2013.10.050 ZHOU Y C, ZHI C Y, WANG L, et al. Health risk assessment of heavy metals pollution in Liujiang River[J]. Anhui Agricultural Sciences, 2013, 41(10): 4370-4372, 4375. doi: 10.3969/j.issn.0517-6611.2013.10.050
[24] 吴东浩, 于海燕, 吴海燕, 等. 基于大型底栖无脊椎动物确定河流营养盐浓度阈值——以西苕溪上游流域为例[J]. 应用生态学报, 2010, 21(2): 483-488. https://www.cnki.com.cn/Article/CJFDTOTAL-YYSB201002035.htm WU D H, YU H Y, WU H Y, et al. Estimation of river nutrients the resholds based on benthic macroinvertebrate assemblages: a case study in the upper reaches of Xitiao Stream in Zhejiang, China[J]. Chinese Journal of Applied Ecology, 2010, 21(2): 483-488. https://www.cnki.com.cn/Article/CJFDTOTAL-YYSB201002035.htm
[25] 曹然, 黎征武, 刘威, 等. 广东北江大型底栖动物群落与水质相关性研究[J]. 华南师范大学学报(自然科学版), 2017, 49(2): 115-121. doi: 10.6054/j.jscnun.2017084 CAO R, LI Z W, LIU W, et al. Relationship between macrobenthos and water quality in Beijiang River Basin[J]. Journal of South China Normal University (Nature Science Edition), 2017, 49(2): 115-121. doi: 10.6054/j.jscnun.2017084
[26] DURAN M. Monitoring water quality using benthic macroinvertebrates and physicochemical parameters of Behzat Stream in Turkey[J]. Polish Journal of Environmental Studies, 2006, 15(5): 709-717. http://www.pjoes.com/pdf-87923-21782?filename=Monitoring%20Water%20Quality.pdf
[27] 张莹, 吕振波, 徐宗法, 等. 山东半岛南部海湾底栖动物群落生态特征及其与水环境的关系[J]. 生态学报, 2011, 31(15): 4455-4467. https://www.cnki.com.cn/Article/CJFDTOTAL-STXB201115032.htm ZHANG Y, LV Z B, XU Z F, et al. Ecological characteristics of macrobenthic communities and their relation to water environmental factors in fourbays of southern Shandong Peninsula[J]. Acta Ecologica Sinica, 2011, 31(15): 4455-4467. https://www.cnki.com.cn/Article/CJFDTOTAL-STXB201115032.htm
计量
- 文章访问数: 489
- HTML全文浏览量: 142
- PDF下载量: 55
