湖泊生态系统灾变早期预警方法技术方案

本发明专利技术公开了一种湖泊生态系统突变早期预警方法，该方法采用生态系统突变检测结果可以获取湖泊生态系统能质出现突变即灾变的时间节点，由此选取生态系统能质突变临界点，结合由生态系统总磷、总氮和能质时间序列构建的湖泊生态系统灾变预警模型，计算得到相应的总氮、总磷浓度即为太湖生态系统预警阈值，其结果更能适用于实际湖泊环境。本发明专利技术可以解决传统基于临界放缓理论的扰动模拟试验方法难以适应现实生态系统干扰因子的问题，是一种更准确、更实用的湖泊生态系统灾变预警方法。

【技术实现步骤摘要】
湖泊生态系统灾变早期预警方法
本专利技术属于水生态系统预警领域，特别涉及一种湖泊生态系统突变早期预警方法。
技术介绍
大型浅水湖泊作为一种复杂的、独特的湖泊生态系统，水体受上边界(水-气界面)和下边界(水-土界面)影响极大，泥水界面相互作用强烈，沉积较为缓慢，物质迁移受风浪影响较大，内源营养盐爆发释放风险高，水环境承载力较低，水生生态系统更脆弱，对污染响应比较敏感，是极易发生富营养化的水体之一。湖泊富营养化后会导致一系列的生态系统异常响应，包括微生物生物量与生产力增加、附着生物逐步增加，蓝藻水华频发、湖泊水质恶化等，当外界驱动条件超过某一阈值时，生态系统可能会由于其结构不稳定和功能退化而产生剧烈响应，发生生态系统灾变。若生态系统发生灾变，藻类生物量剧增，营养盐循环速率加快，水生植物的急剧衰退，水体透明度下降，生物趋于小型化和更短生命周期的生物占优势，整个湖泊生态系统也会出现生物多样性下降、生态系统趋于不稳定的现象。湖泊生态系统灾变的发生往往让人始料未及，产生严重后果，并且难以恢复，人们无法确定湖泊生态系统是否在安全范围内运行，更难以预料生态系统的未来发展方向，若能在生态系统发生灾变前，通过提取生态系统发生灾变的早期信号来预警生态系统状态是否不断接近灾变阈值，并采取相应的措施阻止对人类社会产生不利影响的灾变发生，将具有非常重要的现实意义。生态系统重大突变的早期信号提取存在很大的困难，这主要是由于生态系统发生突变前，系统状态可能只呈现微小的变化，很难察觉和监测。现有湖泊生态系统重大突变预警研究主要基于临界放缓理论，通过测量系统恢复时间，计算系统状态变量的自相关系数、方差、偏度，分析功率谱，观察频繁波动等突变检测方法在时间维度上进行生态系统灾变的早期预警。然而，基于临界放缓理论的系统恢复时间的测量多是通过扰动模拟试验实现，并结合相关模拟模型表征生态系统内部的压力-响应关系，由此得到的结论难以在实际生态系统中得到验证，这是因为现实生态系统的干扰因子远较模拟实验复杂。另一方面，在自相关系数、方差、偏度的计算和功率谱分析过程中，数据获取困难和数据处理过程带入的误差和主观性、外界干扰形态的不确定性等，都会成为造成生态系统突变早期预警方法局限性原因，使得预警信号经常出现误报和漏报。长时间序列数据能够显示生态系统发展过程，有较强的连续性，能够更详细地反映生态系统对外界压力的响应状态。因此，获取长期高分辨率的实际监测数据，将模型模拟、控制实验和现实案例分析结合起来，从多个角度来提取生态系统灾变的预警信号，可以对可能发生的灾变进行更有效地预警。
技术实现思路
本专利技术利用长时间序列湖泊实际监测数据，基于能质理论，建立湖泊生态系统灾变预警模型，获得针对特定湖泊的生态系统灾变预警阈值，解决实验室扰动模拟实验难以适用于实际湖泊环境的问题，为实现“一湖一策”的针对性管理提供依据本专利技术的目的是基于长时间序列湖泊实际监测数据，提供一种进行亚热带大型浅水湖泊生态系统在灾变预警的方法，该方法可用于湖泊流域战略环境影响评价中支持累积性风险评价工作，也可以用于环境保护部门环境风险管理工作的决策支持。本专利技术的技术方案，一种湖泊生态系统突变早期预警方法，包括以下步骤：(1)步骤一、计算生态系统能质和结构能质：利用生态系统各类有机体或有机质的干重含量、信息含量计算生态系统能质和生态系统结构能质；(2)步骤二、以总磷和总氮作为湖泊生态系统灾变的预警指标，连续t检验得到的湖泊生态系统能质突变点为灾变临界值；(3)步骤三、利用湖泊生态系统长时间序列总磷、总氮和能质数据进行多元非线性曲面拟合，构建湖泊生态系统灾变预警模型，得到湖泊生态系统灾变总磷浓度和总氮浓度的预警阈值；通过生态系统能质突变结果获取湖泊生态系统能质出现突变即灾变的时间节点，由此选取生态系统能质突变临界点，计算得到相应的总氮、总磷浓度即为湖泊生态系统预警阈值。其中，能质是指能量使生态系统从有组织的、远离平衡状态达到相对于它所处环境的热力学平衡状态时所能做到的最大有用功；它代表了该系统目前状态与热力学平衡状态之间的距离，表征湖泊生态系统所含的生物量及其所携带的信息量，是预示系统演化发展的一个目标函数。进一步地，所述各类有机体或有机质包括：浮游植物、浮游动物、底栖动物和大型水生植物。进一步地，采用公式a计算所述步骤一中生态系统能质：其中，Wi为生态系统第i种有机体或有机物信息含量，即权重转换因子，J/mg；Ci为生态系统第i种有机体或有机物的干重含量，mg/L；n为生态系统有机体种类总数；Wi值是根据生态系统中某一种生物有机体所含的无重复基因的个数来确定的，不同物种有机体所含的无重复基因数不同，甚至相差很大。对于一般的生态系统而言，系统对环境的响应有使能质变大的趋势，其值越大，表明湖泊对外做功的能力越强，湖泊越健康。进一步地，采用公式b计算所述步骤二中生态系统结构能质：其中，Ct为生态系统中有机体总的干重含量。结构能质可表征湖泊生态系统利用环境资源的能力，反映湖泊生态系统的多样性和复杂性。一般而言，其值越大，湖泊生态系统结构就越复杂。更进一步地，在所述Ex和Exst计算中，所述浮游植物、浮游动物、底栖动物和大型水生植物的权重依次为3.4、36.6、32.5和58.2。进一步地，所述步骤二中突变检测的参数：Targetp＝0.2,cutofflength＝6,Huberparameter＝1。进一步地，所述步骤三中使用湖泊生态系统灾变预警模型分析生态系统能质变化过程，当总磷、总氮小于预警阈值时，生态系统能质随营养盐浓度的增加而增加，即越多的营养盐输入，生态系统健康状态越好；当总磷、总氮等于预警阈值时，生态系统能质为最大值，达到生态系统灾变预警状态；当超过灾变预警阈值后，随着营养盐的增加，水环境生态系统能质出现递减状态，更为频繁和更高强度的蓝藻水华，导致水生态系统群落结构遭到破坏，大型沉水植被退化，生态系统自我恢复能力丧失，从而引发浅水湖泊草型清水稳态向藻型浊水稳态的转变，继而发生生态系统灾变。进一步地，所述湖泊生态系统突变早期预警方法适用于亚热带浅水湖泊。能质：指能量使生态系统从有组织的、远离平衡状态达到相对于它所处环境的热力学平衡状态时所能做到的最大有用功；它代表了该系统目前状态与热力学平衡状态之间的距离，表征湖泊生态系统所含的生物量及其所携带的信息量，是预示系统演化发展的一个目标函数结构能质：表征湖泊生态系统利用环境资源的能力，反映湖泊生态系统的多样性和复杂性。一般而言，其值越大，湖泊生态系统结构就越复杂。生态系统灾变：湖泊富营养化后会导致一系列的生态系统异常响应，包括微生物生物量与生产力增加、附着生物逐步增加，蓝藻水华频发、湖泊水质恶化等，当外界驱动条件超过某一阈值时，生态系统可能会由于其结构不稳定和功能退化而产生剧烈响应，发生生态系统灾变。若生态系统发生灾变，藻类生物量剧增，营养盐循环速率加快，水生植物的急剧衰退，水体透明度下降，生物趋于小型化和更短生命周期的生物占优势，整个湖泊生态系统也会出现生物多样性下降、生态系统趋于不稳定的现象。本专利技术将能质理论用于湖泊生态系统灾变预警，该方法的数据基础是长时间序列的实际监测数据，其结果更能适用于实际湖泊环境。本方法可以解决传统基于临界放缓理论的本文档来自技高网...

【技术保护点】
一种湖泊生态系统灾变早期预警方法，其特征在于，包括以下步骤：(1)步骤一、计算生态系统能质和结构能质：利用生态系统各类有机体或有机质的干重含量、信息含量计算生态系统能质和生态系统结构能质；(2)步骤二、以总磷和总氮作为湖泊生态系统灾变的预警指标，连续t检验得到的湖泊生态系统能质突变点为灾变临界值；(3)步骤三、利用湖泊生态系统长时间序列总磷、总氮和能质数据进行多元非线性曲面拟合，构建湖泊生态系统灾变预警模型，得到湖泊生态系统灾变总磷浓度和总氮浓度的预警阈值；通过生态系统能质突变结果获取湖泊生态系统能质出现突变即灾变的时间节点，由此选取生态系统能质突变临界点，计算得到相应的总氮、总磷浓度即为湖泊生态系统预警阈值。

【技术特征摘要】
1.一种湖泊生态系统灾变早期预警方法，其特征在于，包括以下步骤：(1)步骤一、计算生态系统能质和结构能质：利用生态系统各类有机体或有机质的干重含量、信息含量计算生态系统能质和生态系统结构能质；(2)步骤二、以总磷和总氮作为湖泊生态系统灾变的预警指标，连续t检验得到的湖泊生态系统能质突变点为灾变临界值；(3)步骤三、利用湖泊生态系统长时间序列总磷、总氮和能质数据进行多元非线性曲面拟合，构建湖泊生态系统灾变预警模型，得到湖泊生态系统灾变总磷浓度和总氮浓度的预警阈值；通过生态系统能质突变结果获取湖泊生态系统能质出现突变即灾变的时间节点，由此选取生态系统能质突变临界点，计算得到相应的总氮、总磷浓度即为湖泊生态系统预警阈值。2.根据权利要求1所述的方法，其特征在于，所述各类有机体或有机质包括：浮游植物、浮游动物、底栖动物和大型水生植物。3.根据权利要求2所述的方法，其特征在于，采用公式a计算所述步骤一中生态系统能质：其中，Wi为生态系统第i种有机体或有机物信息含量，即权重转换因子，J/mg；Ci为生态系统第i种有机体或有机物的干重含量，mg/L；n为生态系统有机体种类总数；Wi值是根据生态系统中某一种生物有机体所含的无重复基因的个数来确定的，不同物种有机体所含的无重复基因数不同。4.根据权...

【专利技术属性】
技术研发人员：郭飞，胡艳，孙福红，郄玉，廖海清，冯承莲，王凡凡，吴丰昌，
申请(专利权)人：中国环境科学研究院，
类型：发明
国别省市：北京,11