生物完整性在时间尺度趋势分析中的偏差和校正研究
——以松花江流域为例
2020-09-09李中宇王业耀金小伟
阴 琨,李中宇,赵 然,王业耀,金小伟
(1:中国环境监测总站,北京 100012)(2:黑龙江省生态环境监测中心,哈尔滨 150056)
生物完整性指数(index of biotic integrity,IBI)评价方法,最初于1980s被Karr[1]提出并应用于评价鱼类完整性中. 由于IBI法适用性广、开放性强,很快在更多的生物类群和更大的尺度上得到应用,更被美国和欧盟分别用于支撑清洁水法案(Clean Water Act,CWA)和水框架指令(Water Framework Direction,WFD)对于水环境的管控[2]. 生物类群由最初的鱼类[1]拓展到多种不同生物类群[3],如大型底栖动物[4]、藻类[5]、浮游生物[6]、水生维管植物[7-8];并由单类群发展到多类群组合,如鱼类和藻类[2]. 研究区域也由最初小尺度的溪流拓展到水系[9]、流域[10]、区域[11-13]及国家范围[14]等更大的空间尺度[15-16].
IBI应用至今,也暴露出方法上的缺陷和限制因素,IBI区域特性明显,限制了评价结果在更大空间尺度的可比性[17-18];更重要的是整个方法体系依赖参照状态的确立,而如何定义和确定参照状态还不能完全标准化,由此影响到评价指标建立和评价等级确立. 目前,参照位点法仍是确定参照状态的主流方法之一[9,19]. 考虑到参照点状态有可能随时间不断改善,则需要对参照位点水生态状况随时间发生的变化进行评估分析和校正. 因此,进行不同时间尺度的分析时,参照状态的不一致会导致完整性评价结果不具直接可比性. 对水生态环境稳定、变化潜力极小的区域,较短时间变化的影响可能微乎其微,并不显著. 但现阶段,我国各大流域水体正经历着大规模的水污染治理过程,在全国各大流域处于全面污染防控和治理的作用下,水生态环境发生恢复性变化的速度和潜力都会非常显著,因此水体在相对更短的时间就会发生较为明显的变化. 以“十二五”期间太湖和松花江流域水质变化为例,太湖流域总磷、总氮浓度分别降低了19.2%和32.5%;松花江流域Ⅲ类及以上水质断面比例上升了19.9%,水环境质量发生改善性变化的速度可观[20]. 在水环境快速改善的作用下,水生生物的结构和组成也将随之发生相应的变化. 同时基于水生生物构建的生物完整性指数最初建立的参照状态(参照位点法)在时间尺度上也将随之变化. 在时间尺度的趋势分析中,每年基于变化的参照状态(参照位点)开展评价就必然会导致偏差的产生,这种偏差的大小,对趋势和年度间比较分析的影响程度,及是否可以修正,都会影响时间尺度数据的可比性和趋势分析的科学性.
本文以松花江“十二五”期间(2012-2015年)水生生物调查数据及生物完整性评价结果为例,通过年度间IBI评价结果的差异分析,定量阐释了基于参照位点法的IBI评价在趋势分析中产生的偏差和影响,并提出时间尺度趋势分析中IBI的校正方法.
1 材料与方法
1.1 研究区域概况
松花江流域位于我国东北地区,流域面积约为56万km2,全长达2328 km. 流域内气候整体上属于大陆性季风气候. 流域年平均气温一般在3~5℃之间,有长达5个多月的冰封期[21]. 从“十一五”至“十二五”流域松花江流域水环境呈现持续改善,Ⅰ~Ⅲ类水质断面比例逐年上升,2006-2010年占比在23.8%~47.6%之间,至2015年上升至65.1%[20].
1.2 研究方法
于2012-2015年6-7月连续调查了松花江流域内黑龙江(S1和S2),松花江干流(S3、S4、S11、S12、S13、S14、S15、S28和S30),嫩江(S8、S10、S19、S20和S21),第二松花江(S9、S22、S24、S26和S27)以及支流牡丹江(S16和S29)、梧桐河(S5)、雅鲁河(S17和S18)、诺敏河(S6和S7)、饮马河(S25)、伊通河(S23)(S23位点2012年设在伊通河,由于未采集到底栖生物,2013-2015年调整到汤旺河),共30个采样点,调查河流和点位分布图见图1. 研究是基于同批采样人员在6-7月,以同样的采样方法在相同位点获取的生物调查数据进行,确保数据本身的可比性和一致性. 以研究区域内可获得的最少干扰状态确定参照位点,以不存在明显的人类活动干扰迹象,研究区域上游无污染源,河岸带植被覆盖率高;及美国的快速生物评价方案(RBPs)[22]生境评分中人类活动和土地利用一项得分大于13分,同时有襀翅目昆虫存在为选择原则[23],确定S1~S10为参照位点. 采集流域内底栖动物样品,基于各年份当年生物调查数据对松花江流域进行连续4年的完整性评价. 采用复合生境采样法,底栖动物采集方法和完整性评价方法参见作者前期研究[23].
图1 松花江流域调查点位分布Fig.1 Distribution of the sampling sites in Songhua River Basin
IBI评价的基本步骤如下:初步筛选出具有代表性的大型底栖动物候选参数(25项);以箱线图(box-plot)及IQ值记分法[24]进行参数敏感性分析;以 Pearson 相关性分析进行参数冗余度剔除,筛选核心参数;对核心参数记分统一量纲;计算IBI总分值,以所有位点IBI的95th建立评价标准,完成各位点IBI评价.IBI评价核心参数为总分类单元数、EPT分类单元数、EPT密度、敏感种分类单元%、敏感种密度和Hilsenhoff 生物指数(HBI)6项.
1.3 数据处理分析
研究采用箱线图分析松花江流域IBI指数6个核心参数中正向参数(即干扰越大分值越低的参数:总分类单元数、EPT分类单元数、EPT密度、敏感种分类单元%和敏感种密度)参照组数据在2012-2015年的变化情况(2012年为初始年,2013、2014、2015年为比较年),以及反向参数(即干扰越大分值越高的参数:HBI)受损组数据的变化情况,箱线图数据用Sigma Plot version 10软件统计处理. 采用0~10连续赋分法对各年度核心参数记分[25-26],正向参数和反向参数分别按照公式(1)和(2)赋分计算:
V′i=10(Vi/VR95th)
(1)
V′i=10(1-Vi/VI95th)
(2)
式中,V′i为赋分后的参数值;Vi为赋分前参数值;VR95th为参照位点参数的95th分位数;VI95th为受损位点参数的95th分位数.
参照Mebane等[27]的评价标准划分方法以所有位点指数值分布的95th法建立IBI评价标准,以95th为IBI最优值,低于该值的分布范围进行5等分,以分值从高到低依次定义为:优、良好、一般、较差、很差.
参数赋分中采用校正系数(r)进行参数值校正,并以r值量化显示IBI分值以当年基础值进行赋分在趋势分析中产生的偏差.r为生物参数赋分计算中比较年基础值相对初始年基础值(95th)的倍数,r按照以下公式计算:
r=V(C)95th/V(B)95th
(3)
式中,V(C)95th为比较年参照点(或受损点)参数的95th,V(B)95th为初始年参照点(或受损点)参数的95th.
2 结果
2.1 参数赋分中的偏差分析和校正
2013-2015年参照样点的正向参数总分类单元数、EPT分类单元数、EPT密度、敏感种分类单元%和敏感种密度的值均不同程度地高于初始年(2012年);反向参数HBI则呈下降趋势(表1,图2),表明流域水生生物处于恢复性变化中. 以25th~75th的数值分布状况分析,总分类单元数、EPT分类单元数、EPT密度和敏感种密度4项参数值增加更显著,仅敏感种分类单元%在2013和2014年略有降低;反向参数HBI除2014年数值变化不明显,其余年份均有下降. 表1中校正系数r的结果也指明,除敏感种分类单元%2014年r为0.8,5项正向参数在各年度的r均高于1.0,分布于1.1~4.9之间;反向参数r均小于1.0,分布于0.7~0.8之间,相较初始年,各参数数值在年度间的变化明显,表明IBI以当年基础值进行赋分在趋势分析中会产生显著的偏差. 参数基础值升高,直接引起赋分后参数记分值的相对下降,导致IBI总分的相对偏低,这种偏差使得评价难以捕捉到生物状态改善的变化和程度. 参数赋分计算中,公式应修正为:V′i=10(Vi/VR95th)r(正向参数);V′i=10(1-(Vi/VI95th)r)(反向参数).
表1 2013-2015各年度6项参数的校正系数
图2 核心参数在2012-2015年的变化情况(图中箱体为参数值25th和75th分值;误差棒为参数值10th和90th分值)Fig.2 Changes of core elements from 2012 to 2015
表2 2013年评价标准校正前后IBI评价结果
表3 2013年校正前后IBI评价结果比较
2.2 评价标准划分中的偏差分析和校正
以松花江流域2013年数据为例分析了IBI评价标准建立和划分环节中产生的偏差. 据表2数据可见,初始年(2012年)和2013年经参数赋分校正后所有位点IBI分值的95th分别为46.73、64.83,分值差异明显. 2013年松花江流域IBI评价结果,以初始年和当年标准分别评价的生物完整性状况也呈现较大差异. 表3数据显示,以初始年标准评价,2013年有50%的位点生物完整性达到一般及以上,仅3.3%为很差;以2013年当年标准评价,完整性达到一般及以上的位点占比仅为40.0%,同时很差位点占比升高至16.7%. 以校正后结果进行趋势分析,2013年生物完整性出现较明显的变化,呈现状态很差的区域向较差转变的趋势;以校正前结果判断则未显示一致趋势,而呈现完整性良好和一般的区域向相对更差的状态波动的迹象(如表2中S2、S6和S7位点). 流域水生态环境评估最关注的问题集中在:状况最差区域的范围,最差区域是否出现恢复迹象,良好区域是否出现恶化. 但根据研究数据显示,以上区域的偏差均达到10%以上. 虽然获取更准确的趋势变化需要基于更大时间跨度的持续性数据分析,但以上结果仍足以表明评价环节产生的偏差会影响趋势的判定,需要对评价环节进行校正,以初始年数据建立统一的评价标准.
图3 2012-2015年各年度校正前后IBI为良好及优的结果比较Fig.3 Comparisons of results of IBI before and after calibration in each year from 2012 to 2015
2.3 校正前后IBI评价结果比较
上文利用松花江2012-2015年IBI数据,对参数赋分和评价标准划分两个环节产生的偏差进行了量化分析. 在此,进一步比较了经参数赋分和评价标准校正前后IBI评价结果在趋势分析中的整体差异.
图3呈现的各年度IBI比较结果显示,其中评价达良好及优的数据显示校正前后区域占比的差异达到4.3%~13.3%. 校正前2013年和2015年良好及优的区域占比仅为20.0%,校正后分别为26.7%和33.3%. 以校正前结果进行趋势分析,流域“十二五”期间的生物完整性并未出现可见的改善,而呈现较为明显的下降,与2012年相比,2015年完整性达到良好及优的区域占比降低了10.0%;以校正后结果进行分析,明确显示流域“十二五”期间的生物完整性呈现一定程度的恢复迹象,2014年和2015年完整性达到良好及优的区域占比分别增加了9.3%和3.3%. 完整性很差的数据比较也呈现同样的结论,校正前后区域占比差异为 6.7%~30.0%. 以2013年和2015年数据为例,校正前很差区域占比分别为10.0%和30.0%,经校正后占比分别降至3.3%和0%. 2015年的值校正前后变化较大,分析为2015年S3、S5、S11和S13几个位点生物参数变化幅度较大,EPT和敏感物种密度都明显增加,拉升了当年参数赋分中参照位点95th和评价划分中所有位点IBI95th分值所致. 根据2012-2015年所有位点的EPT密度和敏感种密度两项生物特征参数数据显示,各年度EPT密度分别为46.9、101.9、70.2和574.3个/笼;敏感种密度分别为29.7、67.0、75.4和317.9个/笼,生物恢复趋势明显. 以校正前结果判定2015年水生态质量显著变差的结果与图2核心参数的实际变化冲突. 以上经校正后得出的结论与生物特征参数显示出的实际恢复性变化结果相符,未校正所得结论与生物类群实际变化情况相悖. 综上所述,以上比较结果指明赋分和评价两个环节会产生明显的偏差,偏差很大程度上会直接影响趋势的分析判定,以初始年数据为基础进行校正在完整性趋势分析中是必要的.
3 讨论
3.1 完整性评价在趋势分析中的应用
完整性评价早期阶段,利用IBI表征生态系统状态和健康的研究多是集中在建立评价方法并开展单次评价,随着方法应用发展至今,利用IBI开展水生态环境或生物完整性长时间跨度和趋势变化的研究越来越受到关注,其中更多是鱼类完整性的研究[28-32],也有学者利用浮游生物开展研究[6,33]. 分析趋势变化中,因各研究区的具体情况相异,研究者应用不同的方式建立参照状态,确定参数赋分和评价标准划分规则. 基础数据充足的区域,研究常利用历史调查资料和早期研究数据建立参照状态[30,32];采用非连续赋分法[34]对各参数进行赋分;以Karr[1]早期应用于鱼类完整性的等分法,将IBI最小到最大值的分布范围等分(4~6个等级)建立评价标准,分析研究区中长期尺度下完整性的变化趋势. 对于缺少历史数据的区域,学者采用现场调查法或参照位点法建立参照状态,取各年参数数值分布的最优值作为参照状态,同样采用非连续赋分法和等分法进行参数赋分和评价标准划分[6,28,33]. 基于参照状态法的完整性评价,由于赋分和评价标准的绝对性,趋势分析中可以避免偏差的产生,但目前各个国家都存在不同程度历史数据不足的情况,以历史调查资料和早期研究数据确定参照状态并依据参照状态划定参数期望值和赋分范围是完整性指数方法发展面临的现实阻碍. 对于如中国一样生物和生境等历史性数据相对匮乏的国家,难以应用以上方法建立参照状态并开展趋势分析,更多的研究还需要基于现有调查结果,基于参照位点法开展评价,这种方法也是目前应用较广泛的方法.
3.2 完整性评价在趋势分析中存在的问题
参数赋分环节,Karr[1]的非连续赋分方法仍被部分学者用于完整性分析. 非连续赋分法基于绝对概念建立,其依据位点参数相对参照状态值的偏离程度划分评价等级(3级划分,赋分1、3、5分),绝对的赋分方式不会产生类似参照位点95th(分值会随参照位点年度参数值的变化而变化)相对赋分方法所产生的偏差,但间断赋分法被很多学者认为主观性强,参数赋分受主观因素影响,同时敏感性和识别力较连续赋分方法低. 美国环境保护署(US EPA)于2016年重新修正了参数记分方法,修改原非连续赋分法,采用连续赋分法进行完整性评价[26]. 目前,非连续赋分法仍然在完整性研究中被学者持续使用,但此方法正在被连续赋分方法取代. 0~10连续赋分的参数记分方法具有更高的灵敏性和分辨力,正向参数记分基于参照组的95th分值赋分,反向参数基于受损组95th分值赋分,赋分基于的95th分值是参数记分的基础. 连续赋分法实现了对参数更精确的赋分,可更有效地区分参数值间的差异[25,35]. 但当记分建立于调查可得的位点数据,而位点的生物状况发生改善性(或衰退性)的变化,在进行时间纵向趋势分析时,以上变化将体现为年度间赋分基础值的改变,随之直接影响参数记分值在年度间的可比性及IBI值的可比性,即使评价是基于同样的位点和同样的核心参数,这种影响仍不可避免. 由于IBI整个方法都建立在相对概念的基础上,为确保评价结果的可比性,在参数赋分环节,校正偏差是非常必要的.
评价标准建立环节,除早期基于非连续赋分法的评价标准等分法,基于参照点和所有位点分值分布(25th和95th)的评价标准建立方法也被广泛应用于完整性评价[27,36]. 评价标准的划分多基于当年参照点/所有位点IBI分值分布,即参照点25th或所有位点95th确定. 这种方法相对早期非连续赋分——最高分等分的方式更科学、客观,也更依赖参照状态的确定. 因为这种评分划分基于相对分值(即以当年调查点位IBI值),当点位水生态状况随时间发生变化,进行时间纵向趋势分析时,评价标准同样会由于相对分值的变化而产生改变.IBI各项参数记分和评价标准的变化(如生态改善引起各位点参数值升高),会导致年度记分基础发生不同程度的变化,这些变化在一定程度上会掩盖流域内水生生物变化的趋势和程度. 综上所述,当利用参照位点法进行IBI评价时,需要对各年份的参数记分和评价标准划分环节产生的偏差进行修正,以初始年为基础,对各年度数据进行校准,以消除因记分基础改变对趋势分析产生的干扰.
完整性评价理想的状态是实现基于历史数据法确定参照状态进行完整性分析,在不考虑自然因素引起的变异因素干扰下,可以无需偏差校正直接开展趋势分析,但此方法的局限性也非常明显,且难以解决. 因此,历史数据不同程度的缺乏,使得参照位点法被更广泛地采用. 然而,对水环境压力及水生态都处于恢复性变化的水体基于参照位点法开展完整性评价分析,利用初期参照点的数据持续进行评价显然不合理,因为基于相对干扰较大(较后续年份)的参照位点建立的评价方法,会对结果产生过高评价. 针对年度的完整性分析,利用当年参照点的调查数据建立评价更为科学、客观. 然而,如上文所述,开展时间纵向趋势分析时就必须对参数赋分和标准划分环节产生的偏差进行校正.
3.3 完整性评价在趋势分析中偏差的校正
完整性评价的原理是以建立的参照状态(参照状态下各参数的分值)作为基准,以受损位点参数分值与参照状态的偏离度(分值差异)建立评价标准进行评价. 当建立的参照状态(基于历史数据)参数值是固定值,即不在时间尺度产生变化时,参数赋分计算和标准建立中就不存在变量,如Zhu等[30]在长江流域的鱼类IBI年际变化研究所述方法,直接以早期长江流域鱼类的历史数据建立参照状态和确定评价标准,此类研究方法中完整性评价不存在时间尺度上的偏差,无需校正. 反之,当建立的参照状态(基于参照位点法)的参数值不可避免会随时间发生变化时,此类评价中参数赋分计算就存在变量,在变量作用下,完整性评价将产生偏差,类似方法在黄河[28]、辽河[37]、松花江[38]、太湖[39]、淮河[10]中都曾被采用.
本文提出校正系数r进行以上偏差校正的方法,r可量化表征IBI各参数以当年基础值进行赋分在趋势分析中产生的偏差大小. 本文研究数据显示,2015年5项正向参数的r值均高于1.0,分布于1.1~4.9之间;反向参数r值为0.7,参数值呈现不同程度的偏差(r=1.0表示无偏差),这也表明以底栖动物各特征参数(总分类单元数,EPT丰度和EPT分类单元数,敏感种丰度和分类单元%)直接分析,生物呈现明显正向变化,但未经校正时IBI未能显示一致的结果. 校正前,2015年优良位点仅为20.0%(比初始年低10.0%),最差位点达30.0%(比初始年高13.3%),反向变化的结果与特征参数值的正向变化相悖;经r和标准校正后,2015年优良位点达33.3%(比初始年高3.3%),最差位点为0%(比初始年低16.7%),IBI结果与特征参数的变化趋势相吻合,验证了校正方法的可行性. 与本研究中参数赋分及标准建立原则一致的IBI研究方法,同样可利用校正方法进行分析,如采用最佳期望值[29]、比值法[28]、参照状态95th[38]等基于参照状态分值分布进行参数归一化赋分的IBI评价,可利用校正系数对计算中的变量(参照状态分值)进行校正;基于参照点25th[39]和所有位点95th[28,38]法建立评价标准的IBI也可以初始年数据进行评价标准的校正分析. 综上所述,基于校正的评价方法可以初步解决IBI评价(基于参照位点法)应用于趋势分析时所产生的偏差干扰问题. 但是,本文提出的校正系数法也存在一定的缺陷,未考虑和排除由于自然因素变化产生的偏差,进行长时间跨度的评价时,这些偏差将连同本研究捕捉的干扰类偏差一同被校正系数校正.
3.4 完整性评价的展望
流域生态完整性评价中参照状态的确定是必要前提[40],在完整性评价应用于跟踪长期趋势变化的研究中,除了本研究关注的偏差因素和校正,缺少对参照条件自然变化的评估和校正也一直受到学者们的关注. 在没有人为干扰变化的区域,受到自然气候变异性干扰,IBI分值仍然会产生变化,这类干扰阻碍了利用完整性研究人为干扰的发展和应用[41],特别是在长期趋势分析中,如何有效区分人为和自然因素的影响[42],去除非人为因素在IBI分值中产生的偏差,尚有待方法学的系统研究.
基于参照位点法的完整性评价,除了趋势分析中偏差的影响,还存在对结果过高评估的问题. 目前,学者多依据现有调查数据的高分值分布建立并设置参照状态参数分值,参照状态的高分表示在调查区域可以获得的最优得分[43]. 调查区域的最优状态与理论的参照状态通常是存在差距的,因此当参照状态是基于受到一定程度干扰的位点数据所建立,评价就存在过高评估[28,43]. 这是一个回避不了的问题,而且目前开展的很多相关研究都在一定程度上存在这样的情况,如长江[30]、黄河[28]、辽河[37]、松花江[38]、太湖[39]、淮河[10]等,建议在评价时对研究区开展完整性最优区域的状态和保护目标间匹配度的评估分析,基于评估分析结果对评价标准和结果进行调整,降低和修正过高评价.
完整性评价方法应用至今,由于方法本身的开放性得到广泛应用,但也因而产生难以标准化的方法缺陷,参照状态的确定、参数赋分、评价标准划分等环节的标准,在方法学上仍然有待研究. 对于如何从方法学上解决完整性分析在长期趋势分析中因各种因素产生的偏差仍是值得关注的研究方向.
4 结论
1) 利用IBI(基于参照位点法)进行趋势分析,参数赋分和评价标准划分环节会产生不可忽略的偏差,开展时间纵向趋势研究时必须对参数赋分和评价标准环节的偏差进行校正.
2) 总分类单元数、EPT分类单元数、EPT密度等5项正向参数的r主要分布在1.1~4.9之间;反向参数HBI的r分布在0.7~0.8之间,参数赋分基础值(参照组/受损组参数的95th)在年际间变化显著. 评价标准划分环节的影响表现为评估结果在最差区域和优良区域占比的差异,差异达到10.0%~13.4%,偏差对评价结果产生明显的干扰.
3) 经校正后,IBI评价达良好-优状态的区域占比较校正前高4.3%~13.3%,评价为很差状态的区域占比较校正前低6.7%~30.0%,校正后结果与生物特征参数实际恢复性变化相符,未校正的结果存在相悖情况.
4) 以校正系数和初始年数据对完整性评价关键环节中的偏差进行修正. 基于校正的评价方法可以初步解决将完整性评价应用于趋势分析所产生的偏差干扰问题.
致谢:衷心感谢黑龙江省生态环境监测中心及相关各市生态环境监测中心、吉林省生态环境监测中心及相关各市生态环境监测中心、内蒙古自治区呼伦贝尔生态环境监测站在基础调查数据方面的支持和帮助!