西藏尼洋河水环境特征多元统计分析*

刘海平，钟国辉，叶少文，张良松，张惠娟，范丽卿，巴 桑，李钟杰

(1：西藏大学农牧学院，林芝 860000;2：中国科学院水生生物研究所淡水生态与生物技术国家重点实验室，武汉 430072;3：中国科学院大学，北京 1000494：福建省水产技术推广总站，福州 350003)

*取解释率绝对值最大的作为每一个主成分解释程度较高的水质理化指标，并用加粗来表示.

西藏尼洋河水环境特征多元统计分析*

刘海平1，2，3，钟国辉1，叶少文2，张良松4，张惠娟1，范丽卿1，巴 桑1，李钟杰2**

(1：西藏大学农牧学院，林芝 860000;2：中国科学院水生生物研究所淡水生态与生物技术国家重点实验室，武汉 430072;3：中国科学院大学，北京 1000494：福建省水产技术推广总站，福州 350003)

自2008年11月到2009年10月按照月份监测西藏林芝地区尼洋河水环境特征.结果显示，根据Pearson相关性矩阵，可将尼洋河26项理化指标分为3大类，分别是低度相关性理化指标、中度相关性理化指标和高度相关性理化因子.析因分析显示，可将26项理化指标分为4大类，第1类为尼洋河水体常规的理化指标、第2类为尼洋河水体硬度相关指标、第3类为海拔及第4类为影响水生生物生长的水质理化指标.PCA分析发现，平水期、丰水期以及枯水期水环境特征明显不同，以上3期各自聚为一类，而尼洋河沿程方面，水环境特征则不存在差异.构建了尼洋河海拔、底层水温与相关理化因子的一元回归方程.建议加强对尼洋河水环境的监测工作，积极推进尼洋河水域生态环境的可持续发展进程.

西藏；尼洋河；水质理化；时空动态；主成分分析；析因分析；Pearson相关性矩阵

青藏高原素有“亚洲水塔”之誉，是亚洲十大河流之源，湖泊星罗密布，河流纵横交错，由于受到强烈的地质条件的影响以及独特的气候作用，特别是高原地区，高海拔区域的气候变暖效应较低海拔区域明显，由于海拔所引起的气候变暖效应能够加速改变高山生态系统、水生态系统的生物多样性[1]，随着全球气候变暖，青藏高原冰雪融量以及降水量都在发生改变，导致湖面水位升高和河流径流量增大[2]，这些水源的补给将影响和改变水环境特征[3]，从而间接影响水生生物群落结构和水体的生态服务功能，同时，随着水温的升高，将影响河流水生生物的分布、数量以及水体的健康[4].除此之外，社会发展和人为干扰因子的存在，高原河流、湖泊正在遭受着一系列现实的问题，如：流域用水总量不断上升、污染排放总量持续增加、流域主要水体水质不断恶化[5]，丁海容[6]在对拉萨市城区段水环境进行监测发现，由于拉萨市区内大量未经处理的生活污水及工业废水直排到中干渠，再由中干渠在拉鲁湿地出口全部汇入流沙河，导致流沙河水质现状较差，部分水质超过《地表水环境质量标准》(GB 3838-2002)Ⅲ类水质标准，COD污染超标严重，最大倍数达17.5倍.可见，长期坚持开展高原河流、湖泊的水环境监测工作，不论是从科学研究价值还是社会现实意义角度，都显得尤为重要.本文通过对尼洋河水环境特征的分析，判断尼洋河各类水质理化指标之间的相关性，建立水温、海拔与相关理化指标的一元回归方程，从而为尼洋河水环境的监测提供科学的理论依据，为尼洋河水域生态的可持续发展提供理论基础，也为鱼类保护、水电站建设评估提供科学依据.

1 实验设计和数据处理

关于研究区域概况、采站设置、采样时间等内容参考同系列文章[7-9]，水体理化指标采集和分析方法具体参见表1.

表1 尼洋河水体水质理化指标及代码、单位、保存方法和分析方法

通过SPSS 18.0建立尼洋河水质理化指标Pearson相关性矩阵，用析因分析方法探析尼洋河水质理化主要解释指标，构建海拔、底层水温与相关理化因子的一元回归方程.基于R语言(版本号：2.14.1)，用PCA方法(principal component analysis)对尼洋河水质理化指标时空特征进行分析.

2 结果与分析

2.1 尼洋河水体理化指标相关性分析

选择26项水体理化指标以及采样点海拔信息，采用Pearson相关性矩阵，探析尼洋河水体理化因子之间的相关性，结果显示(表2)，在26项水体理化指标信息中，显著相关项(P<0.05)在10项以下的指标有11项，包括：铵态氮、硬度、铁、pH值、镁、海拔、钙、底层溶解氧、总氮、矿化度和表层溶解氧，可以将它们称为低度影响尼洋河水体理化性质相关因子，其中，与海拔相关的指标有硬度、钙、碳酸氢根、硅酸根、叶绿素a以及总碱度6项.显著相关项(P<0.05)在10～15项的理化因子有9项，包括氯离子、总磷、气温、硝酸根、硫化物、化学耗氧量、钠、碳酸氢根和碱度，可以将它们称为中度影响尼洋河水体理化性质相关因子.显著相关项(P<0.05)在15项以上的理化因子有6项，包括硅酸根、钾、叶绿素a、表层温度、底层温度以及硫酸根，可以将它们称为高度影响尼洋河水体理化性质相关因子，显著相关项最多的是底层温度，达22项.

2.2 尼洋河水体理化指标析因分析

用析因分析法分析尼洋河水质理化指标特征，经最大四次方值法旋转平方和之后，前8个主成分累计贡献率达84.89%，可将26项水质理化指标归纳为4类：第1类为尼洋河水体常规监测的理化指标，如主成分1和主成分2所展示，包括：水体的底层温度、表层温度、气温、钠离子、钾离子、硝酸根、硫酸根、碳酸氢根、氯离子、化学耗氧量、硅酸根、表层溶解氧、底层溶解氧、叶绿素a、总氮、硫化物、碱度和铵态氮等18项理化指标；第2类如主成分3所展示，为尼洋河水体的硬度相关指标，包括：硬度和钙；第3类为海拔，即主成分4；第4类为影响水生生物生长的水质理化指标，如主成分5、主成分6和主成分7所示，包括pH值、总磷、矿化度、镁和铁5项理化指标(表3、表4).

2.3 尼洋河水体理化性质时空特征分析

通过PCA解释尼洋河水体理化性质时空特征，前两个主成分解释率之和达48.7%.尼洋河沿程方面，水体的理化性质较为相似，河道底质方面，流经砂石和粘土的尼洋河水体理化性质也较为相似，但是在季节方面，出现较大的变动，主要表现在冬季和春季(枯水期)的尼洋河水体理化性质归为一类，夏季(丰水期)和秋季(平水期)的尼洋河水体理化性质各为一类，尼洋河水体理化性质与季节存在极大的关联(图1).

2.4 尼洋河水体海拔、底层温度与相关理化因子一元回归方程构建

海拔决定一个地区的温度和光照等环境因素的变化，属于宏观尺度的环境因子，通常情况下，由于海拔的升高，水域温度降低，冰冻期延长[11-12]，物种的丰富度也随之降低[13]，可见，海拔和温度对水环境的影响较大，但尼洋河水体的海拔与温度不存在显著相关性(表2)，同时，为了探究尼洋河海拔与其他相关理化因子的具体关系，构建了尼洋河海拔与相关理化因子的一元回归方程，尼洋河可与海拔构建一元回归方程的理化因子有水体硬度、钙、碳酸氢根、硅酸根和叶绿素a，且均为正相关，所构建的一元回归方程关系显著(P<0.05)(图2).选择底层温度与相关理化因子构建一元回归方程，这些相关理化因子有表层温度、气温、pH值、钠离子、钾离子、硝酸根、硫酸根、碳酸氢根、氯离子、硅酸根、矿化度、化学耗氧量、表层溶解氧、底层溶解氧、叶绿素a、硫、总碱度和铵态氮等18项理化指标，所构建的一元回归方程关系显著(P<0.05)，其中，与底层水温呈负相关的理化因子有碳酸氢根、矿化度、表层溶解氧、底层溶解氧、总碱度5项理化因子，其余13项理化因子与底层水温均呈正相关(图3).

3 讨论

3.1 尼洋河水生生物及环境季节变化基本特征及原因探析

水环境的稳定性随着水体中离子组成及比例的变化而变化[14]，水体中离子的来源(C)有：流域内岩石风化产物(Cw)、人为因素贡献(Canth)、大气降水(Cdry)、大气干沉降(Cwet)、生物圈贡献(Cbio)和物质再循环过程中的净迁移量(Cexch)[15]，可表示为：

C=Cw+Canth+Cdry+Cwet+Cbio+Cexch

表3 尼洋河水质理化指标总方差分解

Tab.3 Total variance explanations for physico-chemical factors of water quality index for Niyang River

理化指标初始特征值提取平方和载入旋转平方和载入合计方差百分比/%累积百分比/%合计方差百分比/%累积百分比/%合计方差百分比/%累积百分比/%Elevation8.25831.76131.7618.25831.76131.7616.14623.63823.638WTS3.93615.14046.9023.93615.14046.9025.45720.98744.625WTB2.5629.85356.7552.5629.85356.7552.4229.31453.939T1.7756.82863.5831.7756.82863.5831.8737.20461.143pH1.6106.19469.7771.6106.19469.7771.7386.68467.828TH1.4255.48175.2581.4255.48175.2581.6236.24474.072Mg2+1.2384.76380.0211.2384.76380.0211.5475.94980.021Ca0.9753.75083.771Na+0.8503.27087.040K+0.5902.27089.310Fe0.4961.90791.217NO-30.4511.73392.950SO2-40.3881.49294.442HCO-30.3341.28395.725Cl-0.3211.23696.961SiO2-30.1880.72597.686TDS0.1450.55898.244COD0.1170.45198.695DOS0.0990.38199.076DOB0.0790.30499.381Chl.a0.0590.22699.606TN0.0380.14699.752TP0.0270.10599.857S2-0.0230.09099.947TA0.0140.053100.000NH+4-N00100.000

表4 尼洋河水质理化指标旋转主成分矩阵*

*取解释率绝对值最大的作为每一个主成分解释程度较高的水质理化指标，并用加粗来表示.

图1 基于PCA分析尼洋河水体理化性质时空特征(图中数字为样点标记，其中1、2、3、4分别为采样点Ⅰ的春、夏、秋、冬季；5、6、7、8分别为采样点Ⅱ的春、夏、秋、冬季；9、10、11、12分别为采样点Ⅲ的春、夏、秋、冬季；13、14、15、16分别为采样点Ⅳ的春、夏、秋、冬季.S1指采样点Ⅰ、S2为采样点Ⅱ、S3为采样点Ⅲ、S4为采样点Ⅳ.第一主成分解释率为33.0%(图中左上图第一条黑色柱所示)，第二主成分解释率为15.7% (图中左上图第二条黑色柱所示)，前两个主成分解释率和达48.7%)Fig.1 Spatio-temporal character for physico-chemical characteristics for Niyang River based on PCA

青藏高原大部分河流是以钙、碳酸氢根为主的河流，青藏高原长江源头则以钠离子、氯离子和硫酸根为主[16]，在这样的水环境基本框架下，降水和融雪则是影响青藏高原河流径流量的主要因素[2].蒲焘等[17]对丽江盆地水环境特征研究中指出，季风期丰沛的大气降水输入对河水离子特征有较为显著的影响，Zhang等[3]指出降水对高原湖泊(纳木错)的离子和总磷的补给有着重要的作用.但是，随着青藏高原气温的逐年上升，大量冰雪融水汇入到湖泊或者河流，截止到2010年，色林错湖面水位由于冰雪融水的输入上升了近8m[18]，自1998-2005年青藏高原部分湖区(兹格塘错和措那湖)面积超过了历来最大面积的25%[19].由于冰雪中含有丰富的离子和矿物质元素[20]，也将随着涓涓细流汇入江河，由于雪水的补给，流域(内华达山脉(Sierra Nevada)附近流域)中离子(硝酸根和硫酸根)的变化呈现出较为显著的时空特征[21].

尼洋河是以雨水和冰雪融水混合补给的河流[22]，水环境特征呈现出明显的季节特征.

图2 尼洋河水体海拔与相关理化因子一元回归方程Fig.2 Regression equation between elevation and relevant physico-chemical factors for Niyang River

通过PCA分析发现，尼洋河丰水期(夏季)、平水期(秋季)以及枯水期(冬季和春季)水环境特征明显不同，以上3期各自聚为一类，而尼洋河沿程方面，水环境特征则无差异.

到了冬季，较雨季而言，降雨量或者融雪量极大的减少，河流处于枯水期，流速降低，流速会直接影响到着生藻类的群落结构，由于水体的冲刷作用，流速快的地方着生藻类生物量较低，另外流速与河床的底质存在着相关性，流速快的地方底质为砂石，流速慢的地方底质则为黏土或者细砂，往往粗糙的基质上着生的着生藻类物种丰富度较光滑的基质上大，尼洋河底质为黏土河段的着生藻类物种丰富度和总丰度较底质为砂石的河段大[8].

在夏季，开始有大量的雪融水和天然降水源源不断的输入，春、夏季比秋、冬季的大型底栖动物总丰度高[9].虽然有营养盐的补给，但是此时水流较急，不适宜大量浮游动物生长和繁殖，导致尼洋河夏季浮游动物生物量、物种丰富度、总丰度较低，其中夏季浮游动物物种丰富度、香农指数和均匀度指数最低，而浮游动物生物量和总丰度则仅高于冬季(归咎于低温).河水的交汇导致水体不稳定，尼洋河浮游动物香农指数和均匀度指数在2个交汇处处于低谷位置.

3.2 尼洋河水生生物及环境空间变化基本特征及原因探析

资料显示，在高海拔山地地带，海拔升高会放大气候变暖效应，加速改变高山生态系统、水生态系统的生物多样性进程[1].本研究区域，尼洋河河源米拉山口海拔逾5000m，汇入雅鲁藏布江处海拔在2900m左右，作为高原河流，与海拔呈显著相关的理化因子有硬度、钙、碳酸氢根、硅酸根、叶绿素a以及总碱度6项，且均为正相关，一元线性回归方程显示线性关系显著(P<0.05)，研究结果未显示海拔与水温之间存在显著相关.

水温受到很多因素的综合作用.海拔和河道坡度对水温产生负面影响[23]，Null等[4]研究加州内华达山脉冷水性河流水温与气候变暖的关系，发现气温每升高2℃，水温将升高1.6℃，同时在春季水温升高幅度最大，气温每升高2℃水温将升高5℃.河流形态和融雪因素控制了水温对气温的依赖程度[24]，融雪通过影响水体物理和化学性质，从而影响水生生物的分布和多样性[25].

图3 尼洋河水体底层水温与相关理化因子一元回归方程Fig.3 Regression equation between bottom dissolved oxygen and relevant physico-chemical factors for Niyang River

在本研究中，与底层水温呈显著相关的理化因子有22项(共26项检测的理化因子)，相关率高达84.6%，与之呈正相关的理化因子有表层水温、气温、pH值、钠离子、钾离子、硝酸根、硫酸根、氯离子、硅酸根、化学耗氧量、叶绿素a、硫和铵态氮13项，相关率高达50%，与之呈负相关的理化因子有碳酸氢根、矿化度、表层溶解氧、底层溶解氧、总碱度5项，以上18项指标与底层水温拟合的一元线性回归方程拟合度好(P<0.05)，其他与之显著相关的理化因子还有钙、镁、铁和总磷，仅总磷与底层温度呈显著负相关，其他3项均呈显著正相关.在设备有限的情况下，可以通过获取尼洋河底层水温，计算得知其他理化性质，有较为广泛的应用前景，也可为尼洋河水生态长期监测提供理论依据.

尼洋河沿程浮游植物和着生藻类多样性(香农指数)特点基本一致，即：尼洋河中游藻类香农指数最大，中上游河段和中下游河段藻类香农指数呈下降趋势[7-8].这与Altormatt等[26]提出的观点一致，即：在河流源头以及高海拔区域，水生昆虫的α多样性最低，水流汇集区域多的区域，水生昆虫的物种丰富度也高，中间海拔区域群落结构较为复杂.

尼洋河下游较其它河段平缓，形成了独特的“尼洋河河谷”风景区，其河道底质以泥沙和黏土为主，具有较大的比表面积，更容易为浮游植物的生长和繁殖提供各类营养盐，同时由于河谷河段水流较其他河段平缓，减少了因水流给浮游植物带来的生存威胁，与海拔因素相比，尼洋河浮游植物物种丰富度或者总丰度对底质的响应更多一些，即以泥沙和黏土为底质的河道浮游植物物种丰富度以及总丰度较以砂石为底质的河道更大[7].

尼洋河作为雅鲁藏布江5大一级支流之一，在工布地区社会和经济发展过程中发挥着重要的作用，目前，多布电站已在建，林拉高速即将竣工，林拉铁路工程正在有条不紊地推进，由此衍生的采砂场和水泥加工基地遍地开花，如何在社会发展过程中保持水生生态系统服务功能，因此对尼洋河水域生态可持续发展提出了严峻的挑战.加强对尼洋河水域生态持续性监测，及时向有关部门反馈水生态变化信息，保证尼洋河流域生态可持续发展，从而推动社会和经济的区域有序发展，这也是本系列论文的出发点和落脚点.

致谢：感谢西藏农牧学院动物科学学院强巴央宗院长，动物科学教研室刘锁珠主任对本项研究工作的全力支持，感谢西藏林芝地区农牧局李耿民书记、嘎路副局长对本项研究工作积极配合，感谢闫宏伟、益西曲云、刘金凤、旦增平措同学协助完成水质理化样品的采集工作.

[1] Pepin N， Bradley RS， Diaz HFetal. Elevation-dependent warming in mountain regions of the world.NatureClimateChange， 2015， 5(5)： 424-430.

[2] Bookhagen B， Burbank DW. Toward a complete Himalayan hydrological budget： Spatiotemporal distribution of snowmelt and rainfall and their impact on river discharge.JournalofGeophysicalResearch， 2010， 115：1-25.

[3] Zhang YL， Kang SC， Li CLetal. Wet deposition of precipitation chemistry during 2005-2009 at a remote site (Nam Co Station) in central Tibetan Plateau.JournalofAtmosphericChemistry， 2012， 69(3)： 187-200.

[4] Null SE， Viers JH， Deas MLetal. Stream temperature sensitivity to climate warming in California’s Sierra Nevada： impacts to coldwater habitat.ClimaticChange， 2013， 116(S1)： 149-170.

[5] 李中杰，郑一新，张大为等.滇池流域近20年社会经济发展对水环境的影响.湖泊科学，2012，24(6)：875-882.DOI 10.18307/2012.0610.

[6] 丁海容.拉萨市城区段水环境污染总量控制研究[学位论文].成都：四川大学，2005.

[7] 刘海平，叶少文，杨学峰等.西藏尼洋河水生生物群落时空动态及与环境因子关系：1.浮游植物.湖泊科学，2013，25(5)：695-706.DOI 10.18307/2013.0511.

[8] 刘海平，叶少文，杨学峰等.西藏尼洋河水生生物群落时空动态及与环境因子关系：2.着生藻类.湖泊科学，2013，25(6)：907-915.DOI 10.18307/2013.0615.

[9] 刘海平，叶少文，杨学峰等.西藏尼洋河水生生物群落时空动态及与环境因子关系：3.大型底栖动物.湖泊科学，2014，26(1)：154-160.DOI 10.18307/2014.0119.

[10] 国家环境保护总局《水和废水监测分析方法》编委会.水和废水监测分析方法：第4版.北京：中国环境科学出版社，2002.

[11] Rundle SD， Jenkins A， Ormerod SJ. Macroinvertebrate communities in streams in Himalaya， Nepal.FreshwaterBiology， 1993， 30(1)： 169-180.

[12] Ao M， Alfred JRB， Gupta A. Studies on some lotic systems in the north-eastern hill regions of India.Limnologica， 1984， 15(1)： 135-141.

[13] Suren AM. Macroinvertebrate communities of streams in western Nepal：effects of altitude and land use.FreshwaterBiology， 1994， 32(2)： 323-336.

[14] 王鼎臣.关于碳酸盐型水质稳定性指数判定法的讨论.水处理技术，1994，20(4)：219-229.

[15] 何 敏.小流域风化剥蚀作用及碳侵蚀通量的初步研究[学位论文].昆明：昆明理工大学，2009.

[16] Huang X， Sillanpaa M， Gjessing ETetal. Water quality in the Tibetan Plateau： Major ions and trace elements in the headwaters of four major Asian rivers.ScienceoftheTotalEnvironment， 2009， 407(24)： 6242-6254.

[17] 蒲 焘，何元庆，朱国锋等.丽江盆地地表-地下水的水化学特征及其控制因素.环境科学，2012，33(1)：48-54.

[18] Meng K， Shi XH， Wang Eetal. High-altitude salt lake elevation changes and glacial ablation in Central Tibet， 2000-2010.ChineseScienceBulletin， 2012， 57(5)： 525-534.

[19] Liu JS， Wang SY， Yu SMetal. Climate warming and growth of high-elevation inland lakes on the Tibetan Plateau.GlobalPlanetChange， 2009， 67(3/4)： 209-217.

[20] 王 平，皇翠兰，刘子东.西藏希夏邦马峰地区雪冰化学特征.环境科学，1988，9(1)：23-26.

[21] Sickman JO， Melack JM. Nitrogen and sulfate export from high elevation catchments of the Sierra Nevada， California.Water，Air，andSoilPollution， 1998， 105(1/2)： 217-226.

[22] 关志华，陈传友.西藏河流与湖泊.北京：科学出版社，1984.

[23] Segura C， Caldwell P， Sun Getal. A model to predict stream water temperature across the conterminous USA.HydrologicalProcesses， 2015， 29(9)： 2178-2195.

[24] Lisi PJ， Schindler DE， Cline TJetal. Watershed geomorphology and snowmelt control stream thermal sensitivity to air temperature.GeophysicalResearchLetters， 2015， 42(9)： 3380-3388.

[25] Slemmons KEH， Saros JE， Simon K. The influence of glacial meltwater on alpine aquatic ecosystems： a review.EnvironmentalScience-Processes&Impacts， 2013， 15(10)： 1794-1806.

[26] Altermatt F， Seymour M， Martinez N. River network properties shape α-diversity and community similarity patterns of aquatic insect communities across major drainage basins.JournalofBiogeography， 2013，40(12)： 2249-2260.

Multivariate statistical analysis of water environment for Niyang River， the branch of the Yarlung Zangbo River， Tibet

LIU Haiping1，2，3， ZHONG Guohui1， YE Shaowen2， ZHANG Liangsong4， ZHANG Huijuan1， FAN Liqing1， BA Sang1& LI Zhongjie2

(1：AgriculturalandAnimalHusbandryCollegeofTibetUniversity，Linzhi860000，P.R.China2：StateKeyLaboratoryofFreshwaterEcologyandBiotechnology，InstituteofHydrobiology，ChineseAcademyofSciences，Wuhan430072，P.R.China3：UniversityofChineseAcademyofSciences，Beijing100049，P.R.China4：FujianMarineProductsTechnicalPromotionStation，Fuzhou350003，P.R.China)

We have monitored and analyzed water environment for Niyang River， the branch of the Yarlung Zangbo River， Tibet， based on seasons from November 2008 to October 2009. The results are as follows： Firstly， based on Pearson correlation matrix， physico-chemical factors for Niyang River can be divided into three categories， which are low correlation physicochemical factors， mid correlation physico-chemical factors and high correlation physico-chemical factors. Secondly， based on factorial analysis， physicochemical factors for Niyang River can be divided into four categories， which are normal physico-chemical factors， hardness relevant physico-chemical factors， elevation， and physicochemical factors affecting growth of aquatic organism. Thirdly， based on PCA， the characteristics of level period， wet season and dry season are significantly different， there is no difference along Niyang River from upstream to downstream. Last but not least， we have structured regression equation between elevation and relevant physicochemical factors， regression equation between bottom water temperature and relevant physicochemical factors. The paper also suggests paying close attention to water environment monitoring， in order to warrantee the sustainable development of Niyang River aquatorium ecology.Keywords： Tibet； Niyang River； physico-chemical characteristics； spatio-temporal dynamic； PCA； factorial analysis； Pearson correlation matrix

J.LakeSci.(湖泊科学)， 2015， 27(6)： 1187-1196

DOI 10.18307/2015.0625

©2015 byJournalofLakeSciences

*西藏科技厅地区基金重点项目、西藏财政厅农业技术推广项目、农业部公益性行业专项(201403012)和国家自然科学基金项目(31160140)联合资助.

2015-05-26收稿；2015-07-20收修改稿.

刘海平(1981～)，男，博士，副教授；E-mail：luihappying@163.com.

**通信作者；E-mail： zhongjie@ihb.ac.cn.