1980-2020年滇池生态脆弱性评价及主要驱动因子.pdf

（湖泊科学），（）：年滇池生态脆弱性评价及主要驱动因子徐灵芝，潘继征，李 勇，华跃洲，李清濯，阳 振，何尚卫，杜成栋，（：苏州科技大学环境科学与工程学院，苏州）（：江苏省环境科学与工程重点实验室，苏州）（：中国科学院南京地理与湖泊研究所，湖泊与环境国家重点实验室，南京）（：滇池湖泊生态系统云南省野外科学观测研究站，昆明）（：广东省佛山市南海区苏州科技大学环境研究院，佛山）（：中国科学院大学，北京）摘 要：人类活动的加剧和经济社会的发展导致滇池开发强度持续增加，滇池生态系统结构与功能受到严重影响，评估滇池的生态脆弱性程度与主要驱动因子是明晰滇池生态系统现状和问题、实现其精准治理和可持续发展的首要任务。基于“暴露程度 敏感程度 适应程度”模型（模型），选取 个准则 个要素 个指标构建滇池生态系统脆弱性评估指标体系，利用逼近理想解排序法（）进行权重方案比选，并通过鲁棒性检验分析，确定计算权重的最优方案。通过分析、年这 个时间段滇池的生态脆弱性，识别出影响滇池生态系统的主要驱动因子，以期为滇池未来生态保护与修复方向的确定提供参考。结果显示，年滇池生态脆弱性呈现先增加后降低的趋势，生态脆弱度最高的是 年（），属于中度脆弱。影响滇池生态系统的主要因素为敏感程度指标，其次为暴露程度指标。在暴露程度方面，影响生态系统的主要驱动因子逐渐从单一的工业污染向工农业的复合污染转变，年工业废水排放量为主要驱动因子，年建设用地面积是主要胁迫因素，年化肥施用量（折纯）为影响滇池生态系统的主要因素；在敏感程度方面，影响生态系统的主要驱动因子逐渐从水质指标向水生态指标转变，年水体恢复能力是主要因素，年与 年土著鱼类种数为主要驱动因子。基于上述分析，为提升滇池生态系统健康水平，应采取改善湖区生境条件、修复水生生物种群、建设农业面源污染控制体系等对策。关键词：滇池；生态脆弱性评价；模型；逼近理想解排序法；鲁棒性检验；驱动因子分析 ，（：，）（：，）（：，）（：，）（：，）（：，）：收稿；收修改稿。江苏省高校水处理技术与材料协同创新中心项目和滇池湖泊生态系统云南省野外科学研究观测站项目联合资助。通信作者；：。徐灵芝等：年滇池生态脆弱性评价及主要驱动因子 ，（），：（），（），（）；，（），：；滇池是我国西南地区第一大湖，是人类活动与发展的热点区域和生态敏感区。以来，随着人类活动的加剧和经济社会的发展，滇池开发强度持续增加，导致湖体水质恶化、动植物栖息地退化、水生生物多样性降低等生态环境问题。自“九五”以来，滇池在污染控制与生态修复等方面开展了一系列水环境治理行动，但目前滇池生态系统功能仍未显著改善，给流域生态安全带来了巨大影响。生态脆弱性研究始于 世纪初期美国学者 将生态过渡带（）的概念引入生态学研究领域。一般而言，生态脆弱性是指生态系统由于暴露在自然和社会环境压力下容易受到干扰的状态，可以认为是暴露程度、敏感程度和适应程度的函数。通过分析生态脆弱性评估结果，可以识别关键的驱动因子，对开展针对性的区域保护和发展至关重要。目前生态脆弱性评估已逐渐从单一生态要素向自然、经济、社会复合的综合化方向发展。同时，生态脆弱性评估也已构建一系列定量评价模型，包括二元统计建模算法、压力 状态 响应（）模型、驱动力 压力 状态 影响 响应（）模型、生态敏感性 生态压力 恢复能力模型、社会生态脆弱性模型等。其中，暴露程度敏感程度适应程度模型（，模型）的评价指标建立过程逐级递进、流程规范清晰，可以揭示自然经济社会的复合影响，近年脆弱性评价也逐渐趋于以 为主要思路的评价体系。基于此，本文以滇池生态系统作为研究对象，建立暴露程度敏感程度适应程度评价指标体系，通过逼近理想解排序法（法）与鲁棒性检验选取最优权重，构建滇池生态脆弱性 评估模型，探究 年滇池脆弱性的变化趋势，从暴露程度、敏感程度和适应程度 个角度分析影响滇池生态系统的驱动因子，为滇池及流域后续生态修复方向及其他湖泊脆弱性评价提供参考依据。研究区概况与数据来源 研究区概况滇池为长江流域金沙江水系的高原淡水湖泊，位于云南省昆明市主城区的南部（，）。滇池南北长约 ，东西宽约 ，湖面面积共 ，库容 亿，蒸发量 亿 。滇池流域为典型的高原盆地地貌，滇池地处流域最低点，是流域唯一汇水和出水通道，入湖河流源近流短（图）。（湖泊科学），（）图 滇池概况 数据来源 以来为保护沿湖农田及村舍，昆明市水利局于滇池开展了环湖防浪堤修建项目，切断了水体与陆地生态连续性，致使沿湖湖滨带生态环境急剧恶化，根据现有的水体数据显示，滇池 年水体营养盐含量呈逐步攀升的趋势；随着 年综合治理滇池的“八五”计划和十年规划的出台，昆明市政府开始了漫长的滇池污染治理工程，但 年滇池水体、及仍呈现波动增加的趋势；直至 年对滇池 条入湖河流污染的治理管控、入草海污染物总量的大幅度削减、滇池水生态修复措施的逐步推进，年水质趋稳向好，水体营养盐含量呈下降趋势。基于以上的人类活动、滇池治理程度以及滇池水环境变化，故将评估阶段分为 个阶段：年、年、年和 年。由于 年时间久远，此时还未开展详细的水质与生态调查，同时该时间段社会经济数据缺失量较多，因此仅对 年、年、年这 个阶段进行脆弱性的研究与外界驱动因子的识别及分析。其中 年作为反映滇池受人类活动破坏剧烈的阶段，年作为滇池污染初步治理的阶段，年作为滇池综合治理的阶段。对滇池生态脆弱性进行评价及分析，其中所用到的数据来源情况如表 所示。其中，环境数据与生态数据采用的是全湖平均值。环保法规条例的颁布情况为定性指标，为便于后续权重的计算，参照云南省河湖库渠健康评价指南（试行）中渠系管理水平的赋分进行定量，无环保法规条例、有法规条例但未完善、有法规条例且已完善分别赋分、和 分。年通过并出台的滇池保护条例为第一部真正意义上针对滇池管理的条例，但由于存在水环境质量标准不明确、管理机构不合理、执法力度不够等问题，其体系尚未完全；年中华人民共和国环境影响评价法与滇池保护条例修订版出台，意味着滇池生态保护体系的法制程度已基本完善。因此以 和 年为分界点：年前，环保法规条例的颁布情况赋分 分；年赋分 分；年赋分 分。表 数据来源情况 数据类型数据名称数据来源自然数据年降雨量、年气温、年日照时间云南省统计年鉴昆明市统计年鉴经济数据人均密度、建设用地面积、耕地面积、化肥施用量（折纯）、工业废水排放量、工业废气排放量、工业固废排放量、人均、第三产业比重、环保投资云南省统计年鉴、昆明市统计年鉴、五年规划环境数据水体恢复能力（）、水质综合污染指数（、）、沉积物重金属潜在危害指数（、）昆明市环境状况公报、滇池管理局生态数据水生植被覆盖度、浮游植物密度、浮游植物多样性、浮游动物多样性、底栖动物密度、底栖动物多样性、土著鱼类种数昆明市环境状况公报、滇池管理局社会数据环保法规条例颁布情况云南省河湖库渠健康评价指南（试行）徐灵芝等：年滇池生态脆弱性评价及主要驱动因子 研究方法 指标体系构建本文基于评估指标的完整性、代表性、可比性和可获取性原则，根据滇池的实际情况，从暴露程度敏感程度 适应程度 方面选择适用性强、代表性高的指标，最终构建了 个要素层 个指标层的滇池生态脆弱性评价体系（表）。表 滇池生态脆弱性评价体系 准则层要素层指标层指标内容与要素层关系暴露程度自然胁迫因子年降雨量 正向年气温 正向年日照时间 正向人为胁迫因素人口密度（人）正向建设用地面积 正向耕地面积 正向化肥施用量（折纯）正向工业废水排放量 万 正向工业废气排放量 万 正向工业固废排放量 万 正向敏感程度水质水体恢复能力（）（）负向水质综合污染指数（、）（）正向沉积物沉积物重金属潜在危害指数（、）正向水生态水生植被覆盖度 负向浮游植物密度（）正向浮游植物多样性负向浮游动物多样性负向底栖动物密度（以寡毛类为主）（）正向底栖动物多样性负向土著鱼类种数 种负向适应程度社会经济发展人均 元正向第三产业比重 正向社会响应环保法规条例颁布情况正向环保投资 亿元正向 暴露程度指生态系统受到外界压力的胁迫程度，主要考虑自然胁迫和人文胁迫因素。滇池作为旱雨季分明的高原湖泊，降雨量、气温、日照时间会通过影响滇池流域水文条件从而影响滇池水环境水生态，故选取降雨量、气温、日照时间作为自然胁迫因子；同时，由于滇池流域环湖开发、贴线建设缺乏有效约束，形成了“环滇发展、环湖开发”的空间形态，与经济发展相关的人口密度、建设用地面积、耕地面积、化肥施用量（折纯）、工业废水排放量、工业废气排放量、工业固废排放量是影响滇池生态脆弱性的重要因素，故将这些因素作为人文胁迫因子。敏感程度反映生态系统内部因素对外界干扰的响应程度，主要包括水质、沉积物与水生态。水质指标是影响湖泊生态脆弱性的基本因素之一，而水体恢复能力与水质综合污染指数是文献中普遍用来评估湖泊水质的指标；滇池作为浅水湖泊，沉积物中高毒性的重金属会影响大型底栖动物的生存与水生植物的根系，改变水体大型底栖动物与水生植物的优势群落，故采用沉积物重金属潜在危害指数评估沉积物；水生态指标是决定生态脆弱性的重要因素，湖泊中水生植被、浮游动植物、底栖动物以及鱼类往往是水生态系统 （湖泊科学），（）的重要组成部分，这些生物丰度以及生物多样性是指示湖泊生态环境的重要因子，且滇池土著鱼类发生显著变化，故以水生植被覆盖度、浮游植物密度、浮游植物多样性、浮游动物多样性、底栖动物密度（以寡毛类为主）、底栖动物多样性、土著鱼类种数作为水生态指标。适应程度指生态系统受到外界胁迫后再恢复的能力，在 模型中主要体现人为的干预（社会经济发展和社会响应）。社会经济发展是间接性表征人为干预、提升滇池生态系统功能的指标，人均 是社会经济发展程度的重要指标，第三产业比重意味着依靠自然资源的社会经济能够的可持续发展程度，第三产业的发展可以降低污染物排放，故将人均 和第三产业比重作为评估社会经济发展的指标。滇池自“九五”以来一直是昆明乃至全国保护的重点湖泊，连续出台各类滇池生态环保保护条例，且环保投资力度较大，故将环保法规条例颁布情况和环保投资作为评估社会响应的指标。指标权重确认为保证各时段滇池生态脆弱性结果的准确性，削弱主观选择权重方法的不确定性，本文通过 法确定权重的最优方案，再利用鲁棒性检验分析该方案的权重，提高结果的可信度。权重方法比选 法能够准确有效地识别理想解，使用广泛，且计算简单方便，因此运用 法选择最优权重计算方案。本文在通过极值法对各指标进行标准化处理的基础上，运用层次分析法（法）、模糊层次分析法（法）、熵权法、灰色关联分析法（法）以及模糊灰色关联分析法（法，包括 欧几里得距离法、熵权法、法）计算不同指标的权重，最后利用 法计算各方法到正负理想解的贴近程度（附表），欧几里得距离法为所有权重评价方法中排序最高的计算方式。权重检验性分析 综合指数的鲁棒性分析是指评估模型参数改变的情况下，检验指标是否能够保持相同排名的方法。虽然本文选择了 法以削弱主观选择权重方法的不确定性，但数据标准化的选择上仍存在一定的主观性，为降低这种不确定性因素带来的影响，利用鲁棒性方法检验权重计算的最优方案，以提高结果的精确度。本文共采用 种不同的标准化方法（极值法、均值化法、四分位距法和 法）、种不同的权重方法（法、法、熵权法、法和 法），共形成 种滇池生态系统各指标权重的排名。根据各指标在 种权重中出现的频数，建立排名频率变化范围矩阵，纵坐标为 欧几里得距离法的各指标排序，频数表示经不同标准化方法及权重方法计算后该排名出现的次数（附表），对角线位置相对其他区域频数较高，鲁棒性检验较强，表明 欧几里得距离法具有较高的可信度与稳定性。因此，最后使用 欧几里得距离法确定权重，各指标最终权重值如表 所示。生态脆弱性评价模型目前 模型常用的计算方法为：（）（）（）（）（）式中，为生态脆弱度，范围为，越趋近于，生态系统越脆弱；、和 分别为暴露程度综合指数、敏感程度综合指数和适应程度综合指数；、和 分别为暴露程度、敏感程度和适应程度中各指标的标准化值；为权重值；、和 分别为暴露程度、敏感程度和适应程度所选的指标数。为明确 年各时间段滇池生态系统脆弱性的变化，需对评价结果进行分级。目前脆弱性评价的分级有等间距法、数轴法、总分频率法、自然断裂法。数轴法与总分频率法呈现正态分布，其中心位置需存在大量的落点，自然断点法因多与 联用需存在较多空间样点。由于本文的滇池生态脆弱性评价仅以 个时间段进行研究且空间样点较少，数轴法、总分频率法、自然断裂法不适用本文的脆弱性分级，因此使用等间距法进行分级评价。具体分级依据：，），微度脆弱；，），轻度脆弱；，），中度脆弱；，），重度脆弱；，极度脆弱。徐灵芝等：年滇池生态脆弱性评价及主要驱动因子表 滇池生态脆弱性评估各指标权重 准则层要素层指标层名称权重名称权重名称权重暴露程度自然胁迫因子年降雨量年气温年日照时间人为胁迫因子人口密度建设用地面积耕地面积化肥施用量（折纯）工业废水排放量工业废气排放量工业固废排放量敏感程度水质水体恢复能力水质综合污染指数沉积物沉积物重金属潜在危害指数水生态水生植被覆盖度浮游植物密度浮游植物多样性浮游动物多样性底栖动物密度底栖动物多样性土著鱼类种数适应程度社会经济发展人均 第三产业比重社会响应环保法规条例颁布情况环保投资 驱动因子识别单参数分析法能够确定每个指标对生态脆弱性的贡献程度，基于此可以提出相应的生态治理和保护建议。计算方法如下：图 滇池生态脆弱性评价结果（年适应程度综合指数为，是相对数）（）式中，为指标 对生态脆弱性的贡献程度；为指标 的标准值；为指标 的权重。结果与讨论 脆弱性评价结果滇池 年各时间段生态脆弱度如图 所示。年滇池生态系统脆弱程度呈先增加后减少的趋势，其中生态脆弱度最高的是 年（），属于中度脆弱；其次为 年（），属于轻度脆弱；最后是 年（），属于轻度脆弱。其中，暴露程度综合指数总体略有上升，敏感程度综合指数先显著增加后又 （湖泊科学），（）下降，适应程度综合指数呈线性上升趋势。对比 与 年的暴露程度、敏感程度与适应程度综合指数，年的生态系统高敏感程度综合指数是造成该时间段内生态系统为中度脆弱的主要原因。相比 年，年敏感程度综合指数由 上升至，增加了约 倍，后因 年水环境污染治理与生态修复等活动的开展，改善了水生态环境，敏感性综合指数由 降至。驱动因子识别不同时间段，影响滇池生态脆弱度的主要驱动因子也存在差异（图）。年，滇池生态脆弱度受暴露程度指标的显著影响，其次为敏感程度指标。暴露程度方面，主要受人为胁迫因子的影响，其中工业废水排放量和工业固废排放量是主要驱动因子；敏感程度方面，主要受水质的影响，其中水体恢复能力是主要驱动因子。年，滇池生态脆弱度受敏感程度指标的显著影响，其次为暴露程度指标。敏感程度方面，主要受水生态的影响，其中土著鱼类种数、水生植被覆盖度和底栖动物密度是主要驱动因子；暴露程度方面，主要受人为胁迫因子的影响，其中建设用地面积和耕地面积是主要驱动因子。年，滇池生态脆弱度受敏感程度指标的显著影响，其次为暴露程度指标。敏感程度方面，主要受水生态的影响，其中土著鱼类种数、浮游植物密度和浮游动物多样性是主要驱动因子；暴露程度方面，主要受人为胁迫因子的影响，其中化肥施用量（折纯）和人口密度是主要驱动因子。图 年准则层（）、要素层（）、暴露程度指标层（）、敏感程度指标层（）和适应程度指标层（）对滇池生态脆弱度的贡献程度 （），（），（），（），（）总体而言，年，滇池生态脆弱度的主要影响因素由暴露程度指标转变为敏感程度指标，由人为胁迫因子转变为水生态因子。暴露程度方面，主要驱动因子由工业废水排放量、建设用地面积转变为化肥施用量（折纯）和人口密度；敏感程度方面，主要驱动因子由水体恢复能力转变为土著鱼类种数、浮游植物密度和浮游动物多样性；适应程度方面，虽然其不是主要影响因子，但是随着滇池流域经济社会的发展以及对生态环境建设的重视，人均、第三产业比重、环保法规条例颁布情况和环保投资这些指标的贡献程度也呈现增加趋势。驱动因子分析通过对 年人为胁迫因子与水质指标的整体趋势进行分析，明确指标的整体变化规律。徐灵芝等：年滇池生态脆弱性评价及主要驱动因子 年滇池人为胁迫因素的变化趋势中（图），建设用地面积、人口密度、化肥施用量（折纯）呈现逐年递增的趋势；工业废气排放量 年前逐年递增，之后呈上下波动的状态；工业废水排放量 年呈攀升趋势，年逐年递减；耕地面积曲线为库兹涅茨曲线（倒 曲线），年达到最大值；工业固废排放量曲线呈现上下波动的状态，最大值位于 年。图 年滇池建设用地面积与人口密度（）、耕地面积与化肥施用量（折纯）（）、工业废水排放量与工业废气排放量（）和工业固废排放量（）的变化特征 （），（），（），（）年滇池全湖 年均值先上升再下降，年呈现波动攀升的趋势，年达到最大值，为 ，年呈现下降的趋势（图）。年滇池 年均值与 表现出相似的变化曲线，年呈现波动攀升的趋势，年达到最大值，为 ；年呈现下降的趋势（图）。年全湖 年平均值波动式下降，年达到最小值，为 （图）。年全湖 年均值波动式下降，于 年达到最小值，为 ，但近 年滇池全湖年平均 浓度为 ，全湖仍存在藻华的风险（图）。年主要驱动因子分析 暴露程度方面，该阶段工业废水排放量与工业固废排放量是主要驱动因子。滇池流域此时处于城镇化发展初期，工业发展较快，同时由于社会对生态环境的忽视导致环保体系的缺失，工业固废、废水等污染物乱排乱放现象频发。敏感程度方面，该阶段水体恢复能力（）是主要驱动因子。滇池呈现富营养化，至 滇池为重度富营养。年滇池草海部分区域水体溶解氧已降至 以下，最低不到 ，草海水体溶解氧仍呈逐级递减的趋势，造成了水体缺氧的状况。但是此时滇池由于早期较好的水质，水体营养盐浓度相对较低，水质综合污染指数在水质因素中占比较低，故溶解氧成为 年影响滇池生态系统的主要水质因素。（湖泊科学），（）图 年（）、（）、（）、（）的年际变化 （），（），（）（）年主要驱动因子分析 暴露程度方面，该阶段建设用地面积、耕地面积与工业固废排放量是主要驱动因子。在此期间，全国正处于城镇化快速发展阶段，滇池流域城镇化水平持续上升，建设用地面积持续至增加。同时，根据图，耕地面积与工业固废排放量的最大值均位于 年。年，国务院将曲靖市的寻甸回族彝族自治县划分为昆明市管辖，这直接导致了昆明市耕地面积与工业固废贮存量的骤然增加，工业固废排放量也因贮存量的增加而升高。根据 年中国人口和就业统计年鉴，年末贵州省总人口数为 万人，与 年同比减少了 万人。同时昆明市五华区启动外来务工人员教育服务管理体系，解决了外来务工人员的教育及就业问题，因此 年末人口数的急剧增加可能与因该教育服务管理体系的启动而迁移至昆明市的贵州省人口有关。敏感程度方面，该阶段土著鱼类种数、水生植被覆盖度和底栖动物密度是主要驱动因子。滇池进行大范围的围湖造堤，导致水生生物栖息地遭到破坏，但水生态指标对环境的响应具有滞后性，故 年水生动植物逐渐成为影响生态脆弱性的显著因素。该阶段滇池水环境处于高营养盐浓度、高藻密度水平，不断从清水草型向浊水藻型退化，生物逐渐由敏感种向耐污种转变，受水质恶化及微囊藻的毒害作用影响，水生植被覆盖度和土著鱼类种数明显减少，高藻密度也为耐污种底栖动物提供了丰富的食物来源。年主要驱动因子分析 暴露程度方面，该阶段化肥施用量（折纯）、人口密度是主要驱动因子。耕地面积自 年后逐年递减，与环滇池生态保护规划 年以及滇池湖滨“四退三还一护”生态建设工作指导意见的通知等生态保护政策的颁布有关；年后滇池流域工业废水及固废的排放量因环保意识的提高、生态保护政策的颁布，自 年后逐年递减。与之相对的是由于化肥施用等导致的农业面源污染存在难监测、难控制、难管理的问题，正在成为滇池流域的主要污染源之一；同时随着社会经济的持续发展，人口不断增加，人口密度逐渐成为影响滇池生态脆弱度的主要因素。敏感程度方面，该阶段土著鱼类种数、浮游植物密度和浮游动物多样性是主要驱动因子。虽然 年间已展开诸多土著鱼类增殖放流的活动，但由于 年刚建立滇池土著鱼类保护研究科技创新团队，鱼类驯化技术、增殖放流及种子工程还待完善，且鱼类的生态恢复需要较长的时间尺度，故土著鱼类种徐灵芝等：年滇池生态脆弱性评价及主要驱动因子数仍维持在极低水平。年滇池全年平均藻密度约为 ，达到了重度水华水平，与 年相比，年浮游植物密度较高。年随着以鱼控藻项目的开展，以鲢、鲫为主的经济鱼类喜食大型浮游动物，同时滇池水体的富营养化程度较高，均导致浮游动物物种趋于单一化、浮游动物多样性下降。建议根据滇池生态脆弱性评价结果，近年来土著鱼类种数、浮游植物密度和化肥施用量（折纯）是滇池生态系统功能紊乱的主要影响因素。因此，未来对滇池的保护建议从以下几个方面出发：）改善湖区生境条件，修复水生生物种群。实施底质改良与地形重塑等措施，利用自然恢复方式恢复滇池湖体水生植被，构建湿地、水域、河漫滩、跌水等丰富多样的鱼类产卵场和栖息繁衍地生境，科学制定增殖放流种类与数量，恢复土著水生生物种群。）建设农业面源污染控制体系。推广测土配方施肥技术，调整滇池流域种植业生产结构，扩大绿色种植、循环农业规模，在有效削减化肥施用量的同时因地制宜地构建源头过程末端全过程的农业面源污染控制体系，促进农业绿色可持续发展。结论）从时间分布来看，年滇池生态系统脆弱程度呈先升高后降低的趋势，其中生态脆弱度最高的是 年（），属于中度脆弱；其次为 年（），属于轻度脆弱；最后是 年（），属于轻度脆弱。其中，暴露程度综合指数总体略有上升，脆弱程度综合指数先显著增加后又下降，适应程度综合指数呈线性上升趋势。）滇池生态脆弱性受敏感程度指标的显著影响，其次为暴露程度指标。暴露程度方面，影响生态系统的主要驱动因子逐渐从单一的工业污染向工农业的复合污染转变，年工业废水排放量为主要驱动因子，年建设用地面积是主要胁迫因素，年化肥施用量（折纯）为影响滇池生态系统的主要外界因素。敏感程度方面，影响生态系统的主要驱动因子逐渐从水质指标向水生态指标转变，年水体恢复能力是主要内部因素，年与 年土著鱼类种数为主要驱动因子。）根据滇池生态系统存在的问题与滇池已开展的生态修复，结合 年滇池生态脆弱性评价的主要驱动力，改善湖区生境条件、修复水生生物种群、建设农业面源污染控制体系是提升滇池生态系统健康水平的主要对策。致谢：本文的数据获取得到了滇池管理局的大力协助，在数据处理、分析及评估过程得到了杨齐老师、周晓辉、徐双、华映肖等人的大力帮助，在此表示诚挚的谢意。附录附表和见电子版（：）。参考文献 ，（）：朱广伟，许海，朱梦圆等 三十年来长江中下游湖泊富营养化状况变迁及其影响因素 湖泊科学，（）：，（）：郑丙辉，曹晶，王坤等 水质较好湖泊环境保护的理论基础及中国实践 湖泊科学，（）：，（）：杨枫，许秋瑾，宋永会等 滇池流域水生态环境演变趋势、治理历程及成效 环境工程技术学报，（）：，（）：何佳，徐晓梅，杨艳等 滇池水环境综合治理成效与存在问题 湖泊科学，（湖泊科学），（）（）：，（）：，（）：，：，（）：，（）：，“”，：，：，：，：李丹 赣南地区生态脆弱性评价及其驱动力研究学位论文 南昌：江西师范大学，：，（），：，：，：，：，：，：鲍姗姗 罗源湾海湾生态系统脆弱性评价学位论文 厦门：厦门大学，徐君，李贵芳，王育红 生态脆弱性国内外研究综述与展望 华东经济管理，（）：，：，（）：杨澄宇，代超，伊璇等 基于正交设计及 模型的湖泊流域总量控制 以滇池流域为例 中国环境科学，（）：李波，曹正浩 滇中引水工程水资源配置方案研究 水利水电快报，（）：徐波，孟端星，沈卫等 近代以来云南生态环境变迁及其社会动因研究 昆明学院学报，（）：何杰，李徽 滇池流域海绵城市建设研究 昆明学院学报，（）：吕忠梅，刘超 水污染治理的环境法律观念更新与机制创新 从滇池污染治理个案出发 时代法学，（）：沈寿文“分工型”立法体制下省市两级地方性法规的分工 基于云南省若干地方性法规的实证研究 云南大学学报：社会科学版，（）：，：，（）：李平星，陈诚 基于 模型的经济发达地区生态脆弱性评价 以太湖流域为例 生态环境学报，（）：，（），（）：徐晓梅，吴雪，何佳等 滇池流域水污染特征（年）及防治对策 湖泊科学，（）：，（）：赵世民，王道玮，李晓铭等 滇池及其河口沉积物中重金属污染评价 环境化学，（）：，：，（）：吴小节，彭韵妍，汪秀琼 中国生态文明发展状况的时空徐灵芝等：年滇池生态脆弱性评价及主要驱动因子演变与驱动因素 干旱区资源与环境，（）：辛岭，安晓宁 我国农业高质量发展评价体系构建与测度分析 经济纵横，（）：，（）：邓雪，李家铭，曾浩健等 层次分析法权重计算方法分析及其应用研究 数学的实践与认识，（）：陈曦，曾亚武，刘伟 基于模糊层次分析法的农村水库大坝安全二级模糊综合评价 水利水电技术，（）：赵航 基于 模型的生态脆弱性评价研究学位论文 北京：中国地质大学（北京），：赵旭 基于模糊灰色关联法的湖南城市土地可持续利用时空特征分析 土壤通报，（）：郭靖，郭生练，胡安焱 基于 法的水文模型多目标参数自动优选方法研究 水电能源科学，（）：，（）：石育中，杨新军，王婷 陕南秦巴山区可持续生计安全评价及其鲁棒性分析 地理研究，（）：，（）：，？，：，：李恒 滇池植被变迁和生态条件的关系 云南大学学报：自然科学版，（）：朱英豪 昆明市土地利用的生态系统价值评价 云南财经大学学报，（）：丁宏翔 年昆明市工业污染物排放状况分析与研究 环境科学导刊，（）：熹微 多措并举构建维护农民工薪酬权益 创造，（）：