承载力最早起源于古希腊时代的工程力学领域, 是指建筑物不发生破坏时的最大荷载能力^[1]。随着资源短缺、生态破坏等问题的出现, 生态学家Park、Burgess和Odum首次在生态环境领域提出了承载力(Carrying Capacity)的概念^[2-3], 研究领域从种群承载力、土地承载力逐渐拓展到区域承载力和生态承载力^[4]。生态承载力不同于单一要素的承载力, 它既考虑自然资源与环境的可持续性, 同时也关注人类社会经济与生态系统完整性、协调性与相关性, 重点关注人类经济社会发展压力与资源环境的关系, 并结合相关要素的动态变化以及管理目标具体分析^[5-6]。目前, 在国际国内, 生态承载力并没有形成科学统一的内涵^[7-8]。随着复合生态系统理论的提出, 生态承载力被认为是在资源合理开发利用和环境良性循环的条件下复合生态系统的承载能力与压力的反映^[7-9], 并逐渐成为评估区域生态环境质量的重要指标, 广泛应用到陆地生态系统功能评价中^[7]。

海岸带是海陆相互作用的过渡地带, 是人类活动密集的生态脆弱带, 也是社会经济最为活跃的区域。全球海岸带以地球表面积的20%养育着全球约50%的人口^[10], 且仍有继续上涨的趋势^[11]。随着人类活动向海洋的不断拓展, 全球范围约41%的海岸带受到人类活动的强烈影响^[12]。我国海岸带地区聚集了大约60%的人口和三分之二的大中型城市。近些年来, 海岸带地区的城市化和工业化进程加速, 然而区域发展缺乏科学规划, 导致过度及不良空间竞争, 资源环境压力大, 区域产业结构不合理, 污染排放总量不断增加, 对近海生态系统产生了日益严重的影响^[13-14]。国内外学者针对海岸带污染胁迫与响应、海域综合评价、生态脆弱性、生态敏感性、生态系统健康等方面开展了一系列研究^[15-19], 其工作重点主要是对海岸带城市环境质量进行承载力的评价, 或者是对地质环境、水产养殖、滨海旅游承载力等某一特定方面进行研究^[20-24]。

然而, 海岸带生态系统属于典型的“社会-经济-自然”复合生态系统, 它的承载力涉及到生态、经济、社会等很多方面, 单一指标、静态分析往往难以评价复合生态系统的的持续发展问题, 所以需要构建指标体系进行多目标分析以及动态综合评价^[25]。因此, 本文提出综合生态承载力(CECC)概念, 指的是在一定时空尺度下, 特定区域的自然资源、环境质量、生产消费、基础设施、科技创新等方面共同支撑区域可持续发展的能力。通过构建P(压力)-S(状态)-R(响应)概念模型, 分析社会-经济-自然不同要素与生态系统综合承载力的关系, 以连云港为例定量评价城市生态承载力, 为沿海地区科学合理地进行产业结构优化调整及海陆环境管理提供技术支撑。

1 研究区域与方法 1.1 研究区概况

连云港位于江苏省东北部, 北纬33°59′—35°07′、东经118°24′—119°48′之间。东濒黄海, 与韩国、日本隔海相望, 北与山东日照市接壤, 西与山东临沂市毗邻, 海域面积6677 km², 海岸线长度204.817 km。

连云港经济发展位居江苏省前列, 2015年连云港地区生产总值、工业增加值、工业用电量、社会消费品零售总额等增幅全省第一, 工业投资增幅全省第三^[26]。连云港作为长江三角洲的重点发展区域, 北接渤海湾, 西依大陆桥, 处于连接新亚欧大陆桥产业带、亚太经济圈、环渤海经济圈和长三角经济圈的“十字结点”位置, 同时又是国家“一带一路”建设重点发展区域, 在全省乃至全国的地理位置与经济地位十分重要^[27]。随着《江苏沿海地区发展规划(2009—2020年)》等战略规划的实施, 连云港的发展迎来了新的战略机遇, 同时也面临资源与环境的严峻挑战^[28]。

1.2 综合生态承载力评价模型 1.2.1 模型构建的假设

本文采用的P-S-R概念模型基于海岸带环境压力、海岸带城市经济社会发展现状以及人类和生态系统共同做出的响应三者间的因果关系链条, 突出展现人类活动对海岸带生态系统的影响。模型的建立主要基于以下三个假设：

(1) 模型由压力指标和单向承压指标共同支撑, 使承载状态保持稳定；(2)模型假设在某一连续的时间段内, 环境、资源、经济要素呈线性变化, 国家和区域政策导向相对稳定；(3)模型假设生态承载力在有单向外力的支持下可以保持健康发展, 当某一方向的外力不足时, 生态承载力降低, 生态系统出现退化。

1.2.2 指标体系的构建

综合生态承载力指标体系应体现连云港作为海岸带城市的发展特点以及社会-经济-自然子系统间的耦合关系, 但任何一个承载力评价指标体系都难以完全涵盖影响该城市承载力大小的各项因素。本文遵循科学性、全面性、层次性、独立性、动态性、可操作性等基本原则, 构建P-S-R模型指标体系。

在P-S-R模型框架内, 海岸带城市生态承载力具体指标的选取主要取决于以下三个方面：

一是压力层面, 反映的是人类社会经济活动对海岸带生态系统造成的负荷, 主要包括当前海岸带城市面临的环境问题及人口压力等, 如万元GDP工业COD排放强度、人口密度等。

二是状态层面, 反映的是生态系统结构功能以及区域环境质量的现状, 主要包括当前海岸带城市的资源环境现状、经济社会发展水平、重点海洋产业的发展能力等, 如海洋渔业总产值、港口货物吞吐量等。

三是响应层面, 反映的是人类面临环境问题时, 为保持和提高生态系统稳定性、服务功能及承载能力采取的必要措施, 主要包括当前国家和地方政府制定的海岸带社会及科技支撑政策, 如污水处理率、研发经费投入占GDP比重、高新技术产业产值^[29]等。

在具体指标的筛选过程中, 为了确保指标选取的科学性、系统性, 避免各指标间的重复性, 在考虑数据可得性的同时, 参考国内外大量文献, 通过实地调研以及专家研讨, 充分考虑海岸带城市的环境压力、资源供给及经济发展特点, 将人类经济社会活动与生态环境变化有机结合在一起, 最终确定了4层23个具体指标(表 1)。

1.2.3 评估数学模型

国内外对生态承载力进行定量评估的模型有海岸带综合管理周期模型、净初级生产力模型、生态足迹模型、状态空间模型等^[30-31]。其中, 状态空间模型可以综合评价社会、经济、自然等多方面的影响因素, 适合本文基于P-S-R概念模型建立的综合评价体系, 因此选用状态空间模型进行数据的定量描述。

状态空间模型中的承载状态点(图 1)(图中A、B、C点), 表示一定时间段内生态承载力的不同状况。任何低于该曲面的点(如图中B、C点)代表某一特定资源环境组合下, 该区域处于平衡或可载状态, 而任何高于该曲面的点(如图中A点)则表明人类活动导致该区域处于超载状态, 生态承载力评估的最终目的就是要在可持续发展的原则下, 使曲面上任意一点的配置达到均衡^[32-33]。

1.2.4 可载性判断

通过比较空间原点到该状态点和状态点在X_max O Y_max上的投影之间的矢量模的大小判定空间内一点同曲面X_max O Y_max间的位置关系, 该区域的综合生态承载能力可由矢量在坐标轴上的模计算得出^[34]。由于同一指标年际变化相对稳定, 且指标之间存在一定的相互作用关系, 因此, 本文拟采用指标的平均值作为其理想值。基于状态空间模型的海岸带生态承载力的实际值与理想值的计算公式如下：

实际值：

(1)

理想值：

(2)

式中：m为指标的个数；ω_i表示第i个指标的权重；χ_i表示第i个指标的值；χ_i表示第i个指标的平均值；CECC表示综合生态承载力实际值, CECC′表示综合生态承载力理想值。

根据状态空间模型, 将生态承载力状况分为3个级别, 即可载、满载、超载。生态承载力综合判定标准如下：

当CECC＞CECC′时, 可载；

当CECC′-r≤CECC≤CECC′+r时, 满载；

当CECC＜CECC′时, 超载。

其中, 根据海岸带区域综合承载力已有的容差研究结果^[21], 结合评价区域调研及专家建议, 将容差r取0.1。

1.3 数据来源与处理 1.3.1 数据来源

本文通过选取2005—2014年10年间连云港市的社会、经济、环境等方面的数据构建综合生态承载力评价模型, 为确保评估体系的类比性及可靠性, 本文的数据主要通过以下几个渠道获得：生态环境数据主要来源于《江苏省环境统计年鉴》、《中国城市统计年鉴》及部分以往的研究数据；社会经济数据主要来源于《江苏省统计年鉴》、《连云港市统计年鉴》、《连云港市国民经济与社会发展统计公报》以及各年度的《连云港市政府工作报告》等。

1.3.2 数据预处理

指标层中的各项指标相互之间在数据特征、数量级及量纲上存在着差异, 为了保证指标间可以进行对比, 运用极差法对指标进行无量纲化处理。根据指标内容, 将指标分为成本性指标与效益性指标两大类, 成本性指标为压力层下的指标, 与其属性数据成反比关系, 效益性指标与其属性数据成正比关系, 包括状态层和响应层下的指标。具体计算过程如下：

假设有m个指标的数据集(X_i)_m, 对于成本性指标, 有：

(3)

对于效益性指标, 有：

(4)

式中, Y_i为指标归一化化值, X_i为指标初始值, max(X_i)和min(X_i)分别为指标数据的最大值和最小值。当Y_i=0时, 指标状态最差, 当Y_i=1时, 指标状态最好。

1.3.3 指标权重的确定

某一指标的权重是指该指标在整体评估中的相对重要水平。常用的权重确定方法分为主观赋权法和客观赋权法两类。两种方法各有优劣^[34]。因此, 本文采用主客观组合赋权的方法, 通过德尔菲法与熵值法分别计算指标权重, 在吸纳专家经验和建议的同时, 结合数据本身的属性, 从而尽量避免两种方法的缺点, 充分发挥其优势, 计算其组合权重。

(1) 熵值法。熵值法赋权已相对成熟, 在此不再对其具体过程赘述。

(2) 德尔菲法。具体的实现过程如下：

① 从环境科学与生态学领域邀请10位专家进行首轮打分, 对分数进行统计分析, 采访部分专家听取其意见和建议。然后, 将首轮专家打分的统计分析结果整理出详细报告, 再次征询专家意见并依据反馈报告修改意见, 经过若干轮反馈, 专家意见逐渐趋向统一。

② 计算权重。计算专家对各指标赋分的平均值, 再对其进行归一化处理, 具体计算过程如下：

(5)

(6)

式中, d_i表示10位专家对第i个指标赋分的平均值, x_ji表示第j个专家对第i个指标所赋分值。

本研究结合上述两种方法, 通过组合权重的计算方法, 保证指标权重的主观性与客观性, 计算方法如下：

(7)

式中, 组合权重为w_j^*, 熵值法所赋权重为w_jm, 德尔菲法所赋权重为w_jn, α为熵值法权重占组合权重的比例, (1-α)为德尔菲法权重占组合权重的比例。根据迟国泰等人^[35]的研究结果, 最佳的组合权重结果是主观权重和客观权重各占50%, 本文α取0.5。

1.3.4 数据分析方法

(1) 贡献度判断

为具体分析引起连云港综合生态承载力变化的原因, 本文引入“贡献度”的概念。所谓“贡献度”, 就是指不同因素层指标的承载力占综合承载力的百分比, 可用于分析导致承载力满载或超载的具体因素。因素层的指标贡献度越大, 对生态承载力有促进作用；反之, 对生态承载力有抑制性作用。公式如下：

(8)

(9)

式中, χ_i表示第i个指标的值, j表示因素层中具体指标的个数, CECC_k表示因素层指标的承载力, R_k表示因素层指标对综合生态承载力的贡献度。

(2) 相关与因子分析

相关分析是一种研究两个或两个以上同等地位随机变量间的相关关系的统计分析方法, 可以体现总体间的因果关系, 其结果表示变量间相互关系的密切程度。本研究通过相关分析, 探究综合生态承载力与因素层指标间的依赖性, 分析影响承载力的关键因素。

因子分析是指研究从变量群中提取共性因子的统计技术, 可用于确定多变量之间的联系。因子分析的结果将相关比较密切的几个变量归位同一类因子并找到有代表性的变量。通过因子分析, 我们可以确定变量与因子之间的相关系数, 从而反映变量对于总体的重要性。

本研究中, 相关性与因子分析均通过SPSS 22.0软件进行分析计算, 计算过程不再赘述^[36]。

2 结果与分析 2.1 综合生态承载力分析

本文以连云港市为评价单元, 计算2005—2014年间的综合生态承载力, 得到连云港市10年间综合生态承载力的变化趋势, 结果如图 2所示。

从图 2中可以看出, 连云港市海岸带综合生态承载力总体呈现波动上升趋势, 10年间的平均增幅达到10.4%。其中, 2007—2008年涨幅最大, 同比增长46.59%, 2006—2007年的涨幅也达到了26.64%。2008—2011年的四年时间里其承载力出现不稳定波动, 并呈现小幅下降, 2011年后承载力空间重新呈现上升趋势并趋于稳定, 但增长速度有所减慢, 这种波动式阶梯逻辑斯蒂上升模式符合多数城市的发展模式^[37]。

根据综合生态承载力判定标准, 设置生态承载力理想值的误差线为0.1, 可以得出, 2005—2007年连云港海岸带生态承载力小于理想值, 处于超载状态, 2008—2011年的生态承载力基本处于平衡状态附近, 2012年后的生态承载力大于理想值, 处于可载状态。

2.2 贡献度判断

本文通过贡献度分析, 分别探讨因素层6个指标影响综合承载能力的差异。

压力层两个影响因素属于成本性指标，其贡献度如图 3所示。10年间连云港海岸带环境压力的贡献度呈现先上升后下降的态势, 2009年环境压力的贡献度最大。这说明随着初期沿海工业的不断发展, 各类污染物的排放成为破坏生态系统平衡的重要因素。随着“十一五规划”的实施, 倡导沿海地区转变经济增长方式，实现可持续发展。2007年后连云港市大力发展节能环保产业、推广低碳经济, 开发清洁能源, 加强水源保护及海洋综合管理并疏散中心城区人口, 使污染排放对生态承载力的影响逐渐减弱。但是, 2009年后环境压力的整体贡献程度进一步下降, 说明环境压力综合指标逐渐减小, 各类环境污染物排放再次增加, 抑制了生态系统承载力的上升趋势。另一方面, 海岸带人口聚集度的贡献度在2006年小幅上升后, 总体呈下降趋势, 2012年后趋于平稳。在2008年之前, 海岸带人口聚集度大于海岸带环境压力的贡献度, 说明2008年之前发展压力主要来源于人口, 而在2009年以后连云港逐渐控制人口增长, 虽然旅游人次有所增加, 但总体趋于稳定。随着沿海工业的迅速发展, 2008年后生态系统的主要压力来源逐渐变为环境压力并且始终大于人口聚集度, 再次印证各类工业污染物的排放成为限制连云港生态环境可持续发展的重要因素。

状态层两个影响因素的贡献度如图 4所示。2005—2014年连云港市海岸带经济发展水平尤其是港口经济的贡献度呈现上升趋势, 2014年达到29.19%。而海岸可利用资源的贡献度呈波动下降趋势, 从2005年的32.08%下降到2014年的9.84%, 其中以人均可利用水资源的减少最为显著。这说明连云港作为沿海港口城市, 其海洋经济产业及海岸带资源拥有重要地位, 必须合理开发并加以保护。在2010年之前, 海岸带可利用资源的贡献度大于海岸带经济发展水平, 说明2010年之前连云港的发展以海岸带自然资源为依托, 为整体经济发展提供了重要保障, 但是随着海岸带可利用资源的减少, 其经济发展模式逐渐转变, 已经不单纯依赖于资源支撑, 逐渐转向科技支撑。

图 5显示了响应层两个影响因素的贡献度, 结果表明, 10年间连云港社会支撑条件的贡献度先上升再缓慢下降又上升, 2008年最高为32.2%, 2010年为15.39%, 2014年又上升到23.21%, 说明2008年之前连云港大力发展城市基础设施建设, 拓展城市空间, 通过积极对接国家推动江苏沿海开发的战略部署, 广泛招商引资, 城市化建设取得突破性进展, 社会支撑条件迅速提升, 但是随着基础设施的日益完善, 其对于综合承载力的贡献度逐渐下降。2012年后, 由于注重污染物综合治理与利用、外资引进等方面的工作, 社会支撑条件的贡献度小幅上涨。而科技支撑条件的贡献度呈持续上升趋势, 2014年达到13%。连云港针对环境污染和生态破坏的压力, 在国家和区域政策的支持下, 进行了经济转型和产业结构调整, 注重高新技术产业的发展, 推行科技强市的战略方针, 其创新研发投入与高新技术产业产值逐渐增加, 为生态承载力的提升打下了基础, 但是目前的整体科技支撑条件仍相对薄弱, 如研发经费投入与高新技术产业产值占比不高，有进一步提升的空间。

2.3 综合生态承载力的关键因素分析

为进一步分析影响综合生态承载力的关键因素, 对因素层的6个二级指标进行Spearman秩相关分析, 结果如表 2所示, 在相关分析的过程中, 我们发现状态层的海岸带经济发展水平与响应层的科技支撑条件与综合生态承载力呈极显著正相关, 而海岸带经济发展水平、社会支撑条件和科技支撑条件三个变量相互之间可能由于存在一定的自相关性导致其存在极显著正相关关系, 在此通过因子分析进一步判断。

在进行相关性检验的基础上, 对因素层的6个指标进行因子分析, 进一步验证结果的准确性, 得出KMO=0.609, P < 0.05, 数据满足因子分析条件, 结果如表 3所示, 选取特征根λ>1, 可以得到两个主成分, 其累积贡献度为85.74%, 分析因子载荷矩阵可得海岸带经济发展水平、海岸带可利用资源、社会支撑条件和科技支撑条件具有一定的自相关性, 同为第一因子, 根据因子载荷矩阵(表 4)的计算结果, 科技支撑条件和经济发展水平的因子载荷量最大，分别为0.980和0.951, 因此在正向指标层面上, 科技支撑条件和海岸带经济发展水平对于综合生态承载力的影响居于主导地位, 而负向指标层面上, 人口聚集度对于综合生态承载力的影响并不显著。因此, 结合贡献度分析及相关分析的计算结果, 可以明确得出海岸带环境压力、科技支撑条件及经济发展水平在未来将成为制约连云港市综合生态承载力的三个关键因素, 需要在未来的政策调整中给予重点关注。

3 结论与建议

本文充分考虑人类社会经济活动对海岸带生态系统造成的各种影响, 建立了海岸带城市综合生态承载力评价的具体指标, 将P(压力)-S(状态)-R(响应)综合评估模型与状态空间模型相结合, 以连云港市为基本评价单元, 通过贡献度分析、相关分析及因子分析, 提出海岸带城市综合生态承载力及其关键影响因素评价方法, 为海岸带城市综合承载力研究提供了新的研究思路。

首先, 本文通过构建多重指标体系, 结合熵值法与德尔菲法, 对承载力的指标权重进行组合修正, 保证评价结果的准确性与科学性。研究结果表明：(1)连云港海岸带综合生态承载力呈现阶梯逻辑斯蒂上升趋势。(2)2005—2007年生态承载力处于超载状态, 2008—2011年基本处于平衡状态附近, 2012年后处于可载状态。(3)虽然2011年以来生态承载力呈现明显上升趋势, 但是其增长速率在逐渐减慢。

其次, 本文对因素层指标进行贡献度分析, 探讨不同因素对综合承载力的影响强度, 可以得出：连云港市海岸带环境压力和人口聚集度先降后升, 海岸带可利用资源波动下降, 海岸带经济发展水平、社会支撑条件以及科技支撑条件整体呈上升趋势, 是综合生态承载力上升的关键因素。结合相关分析和因子分析, 得出海岸带环境压力、科技支撑条件及经济发展水平是制约综合承载力的关键因素。

从评价结果上看, 连云港综合生态承载力变化趋势呈现“U型”, 即“先破坏-后治理”的模式。作为我国重要的沿海城市, 连云港经济社会发展模式正是我国众多海岸带城市的典型代表, 为了避免生态承载力再度下降, 同时为其他海岸带城市提供借鉴经验, 提出以下政策建议。

(1) 为减小海岸带环境压力, 应当统筹推进陆海环境综合治理, 强化主要污染物减排, 建立完善污染物排放总量控制制度, 推进沿海地区重点耗能企业的改造, 淘汰落后产能, 大力发展循环经济和节能环保产业, 以持续改善环境质量, 使生态承载力进一步提升。

(2) 针对海岸带可利用资源减少以及科技支撑能力不强的现状, 应充分运用沿海地区良好的地理位置和资源条件, 大力发展蓝色海洋经济, 优化海洋经济布局, 以连云港为核心的江苏沿海港口群建设推动临港工业发展, 合理开发促进资源集约高效利用；深化落实创新驱动发展战略, 继续保持科技支撑条件稳定提升, 增加研发经费占GDP的比重, 大力发展科技创新企业, 使沿海地区形成以创新为导向的经济发展方式。

(3) 为保障社会支撑条件继续改善, 应当进一步完善沿海交通、绿化、水利等基础设施建设, 构建规划合理、性能完善、保障有力的绿色基础设施支撑体系, 提高生态承载力。