引言

国土空间是人类活动参与下空间域沿时间轴发展形成的动态、多维、复杂的人地关系时空系统，土地是其载体和空间依托^[1]。当前国土空间格局中存在空间开发失衡、资源约束趋紧、区域差距扩大、产业布局欠妥、“三生”空间冲突加剧、陆海开发缺乏统筹等问题^[2-4]，阻碍人与自然和谐发展。

国土空间开发适宜性由一定地域空间的资源环境承载力、经济发展基础与潜力决定，指的是国土空间承载城镇化和工业化发展的适宜程度，包含人类社会安全与公正、经济效率和生态安全三方面要义^[5]，可以明确识别未来土地保护的制约因素和发展机遇，从而有效促进国土空间开发格局优化，构建合理高效的开发利用格局，为解决上述问题提供科学依据。目前，国土空间开发适宜性是编制国土空间规划的重要支撑依据之一^[6]。国土空间开发适宜性评价视角主要分为农业、城镇建设两个方面，其中建设用地适宜性评价内容逐渐由微观的施工难易程度发展为宏观的土地利用需求和建设用地潜在数量供需分析^{[7, 8]}。学者从自然条件、社会经济发展水平、生态环境安全、区位条件等多角度选取评价指标^[9]，利用统计数据、地图信息数据、观测数据、社会调查等通过木桶原理、“千层饼”模型^[10]、层次分析法^{[11, 12]}、专家打分法^[13]、矩阵匹配法^[14]、多要素加权叠加分析^{[15, 16]}、矩阵列联分析^[17]等方法对评价单元进行综合分析。

沿海城市是地球表面陆域与海域区位优势的集合体，也是目前人口以及社会经济活动较为集中的区域，全球沿海地区以地球表面20% 的面积承载着全球约50% 的人口^[18]。沿海城市土地利用既是可持续社会经济发展的基础和保障，又影响着陆域、河口、河岸生态环境^[19]。江苏沿海地区处于我国沿海、长江和陇海兰新线三大生产力布局主轴线的交汇区域，具有优越的地理区位条件和国家战略支持。虽然与全国平均水平相比，江苏沿海地区经济社会发展基础较好，但与江苏省发展水平相比仍存在一定差距。2018年南通、盐城、连云港沿海三市土地面积占全省30.95%，人口占比27.09%，但是对江苏省地区生产总值的贡献率仅为18.02%，根据第三次土地调查数据，2018年沿海三市国土空间开发强度为19.40%，单位建设用地面积经济产出仅为全省水平的59.66%，国土空间开发利用与其经济发展之间存在明显的不匹配，国土空间开发利用需要进一步优化。

目前，国外海岸带研究侧重于技术突破、综合管理及公众参与，研究重点主要有海岸带综合管理、海岸带海陆交互作用、滩涂生物技术、浅海滩涂自然保护及保护区建设、滩涂生态及生态开发模式、滩涂围垦水利工程设施^[20]等。国内学者主要着眼于两方面：一是海岸带土地利用结构特征、动态监测、时空动态特征及其驱动因素，并探究海岸带土地利用变化带来的生态环境影响^[21]；二是脆弱性评估、景观生态评估、生态系统服务价值评估等生态环境现状及价值视角的评价研究^[22-24]。

以往有关国土空间开发适宜性评价以及沿海城市土地利用的相关研究中，许多学者采用不同的方式进行陆海统筹评价，为研究提供了理论及实践基础。冯佳凝等^[25]通过统筹陆域与海域评价分类方案，综合考虑土地资源、水资源、生态环境、区位优势以及海洋资源各方面因素，对南通市国土空间开发适宜性进行评价。刘柄麟等^[26]通过对海岸带陆域及海域分别进行国土空间开发适宜性评价，并建立判别矩阵识别功能空间，进而为海岸带空间分区管控提供参考。然而，沿海城市的国土空间开发适宜性评价仍然缺乏对人海相互作用的关注，陆域评价体系中缺少海域国土空间指标因子。考虑到人类活动通过改变土地利用方式对海洋生态及环境造成影响，海洋灾害等自然因素对国土空间开发构成限制作用，海洋资源及产业对沿海区域经济社会发展发挥支撑作用，本文以海岸带类型较为独特、海涂比例占比高的江苏省沿海为例开展研究，从自然环境适宜性、资源环境承载力两个维度构建国土空间开发适宜性评价指标体系，融合地图信息数据和统计数据，在陆域国土空间开发适宜性体系中增加对海域因素的考虑，进行沿海城市国土空间开发适宜性评价。

1 沿海城市区域人地关系作用机理及适宜性评价方法 1.1 沿海城市区域人地关系作用机理

资源环境承载力具有基础性、有限性、空间可调控性、动态变化性、区域关联性五维特征^[27]，是合理安排人类社会经济活动的重要参考依据。资源环境承载力是区域资源环境与人类社会两大系统之间相互耦合关系的体现，除受资源环境要素的保有量影响外，人类社会对资源的利用方式及效率也对承载力具有显著影响。自然环境适宜性则更强调自然资源本底条件对土地开发利用的限制及约束。本研究将沿海城市人地关系之间的作用划分为三大作用力，主要包括自然环境中地质灾害及海洋灾害等“损害力”，基本农田、生态保护红线、规划用途管制区等“约束力”，资源环境系统对开发建设的“支撑力”，三方面因素共同影响区域国土空间开发的适宜性，其作用机理如图 1所示。支撑力表示人类根据生产生活活动对自然资源的需求进行国土空间开发，但为保证自然资源系统各要素供给及生态保护等重要功能，又要求开发建设活动合理适度，同时提升面临环境问题的响应能力；损害力表示沿海城市自然环境中的各类灾害发生对系统正常功能产生的负面影响；通过划定基本农田、生态保护红线等刚性约束手段，对国土空间开发产生约束力，有利于减少自然环境中各类灾害发生的概率，从而减少其对自然资源系统可持续性及人类经济社会活动产生的负面影响。

1.2 国土空间开发适宜性评价方法

沿海城市的土地开发利用强度、区域结构和开发方式等被自然及生态条件约束。沿海城市经济发展迅速、人口集聚度高，该区域的国土空间开发相比于陆域其他地域单元，除了受到如地形、灾害等陆域自然环境因素影响以外，还受到海洋资源环境（如海平面上升导致的岸线侵蚀、岸线资源保护、海洋生态系统等）的限制，与此同时人类对海岸带的开发利用活动（如化学需氧量、氨氮、工业废水等污染物的排放）也给沿海城市环境施加了一定的压力。基于此，本研究从自然环境适宜性、资源环境承载力两方面构建评价指标体系，而后通过叠加得到国土空间开发适宜性评价结果。

自然环境适宜性评价指标主要对国土空间开发适宜性起限制性作用，综合考虑土地资源、海洋资源和生态环境等多方面因素构建指标体系^[28]，自然环境适宜性评价指标体系如表 1所示，对以下评价指标使用木桶原理，各栅格得分计算公式为：

式中，f(x) 为自然环境适宜性得分；x_i为栅格图斑各评价因子得分。

资源环境承载力对国土空间开发适宜性起到上限约束作用，代表区域所能承载的最大经济社会活动。江苏沿海地区土地资源相对丰富，地处“一带一路”倡议建设、长江经济带发展、长三角一体化发展三大国家战略交汇区域，近年来在加快建设海洋强省等政策背景下具有极大的开发潜力。同时，沿海地区由于位于陆海交互作用地带，生态价值显著，区域内分布有丹顶鹤、麋鹿等珍稀动物保护区，海州湾国家级海洋公园，蛎蚜山国家级海洋公园，小洋口国家级海洋公园等自然保护地，因此江苏沿海地区人类社会经济活动同时应考虑对自然生态功能的影响。

本文环境承载力评价采用“压力—状态—响应”（P—S—R）模型构建指标体系。关于人海关系，已有学者进行相关研究，刘天宝等^[29]将人海关系的理解分为三种类型，分别为资源与生态领域的海洋（海洋的开发与保护等）、空间与权力的海洋（地缘政治思想等）、文化观念的海洋（海洋文明与海洋观）。本文从人海辩证关系思考海洋，赋予“社会性”的特征，海洋可以为人类所利用，提供资源、空间等发展要素，同时构成人类活动的约束。

国内学者在资源环境承载力评价指标方面考虑的因素主要包括水土资源、经济、社会、生态方面，选取的指标包括经济[人均国内生产总值（GDP）、第二、三产业比重、城镇化水平、人均收入等]、资源（水资源量、人口密度、单位产值用水/电量等）、环境（污染物排放量、污染治理情况、森林/绿化覆盖率）^{[30, 31]}。本文结合已有研究，选择如表 2所示的指标进行资源环境承载力评价。压力层面用以反映经济社会活动对资源环境及生态系统造成的压力，选择万元GDP化学需氧量排放量、万元GDP工业废水排放量、人口密度等作为评价指标。状态层面用以反映研究区经济发展结构水平及资源环境现状，该层指标主要包括经济社会发展情况、经济结构等。响应层面用以反映面对环境问题时，区域社会经济及生态环境系统为减少环境问题所带来的负面效应所能给予的支撑能力，该层指标主要包括工业废水处理率、废气治理设施处理能力、一般工业固体废物综合利用率、林地覆盖率、城市绿化覆盖率及实际使用外资额。

资源环境承载力评价中压力层数据为负向因子，状态层、响应层为正向因子，因此对负向因子和正向因子采取不同的标准化方法。

负向指标标准化公式如下：

正向因子标准化公式如下：

式中，y_max为该项因子中的最大值；y_min为该项因子中的最小值；y_i为该项第i个因子值；y_i为y_i标准化后的值。

对标准化后的各层因子进行主成分分析，以各因子旋转后的方差贡献率为权重分别计算压力层、状态层、响应层的综合得分，然后使用聚类分析将各层结果划分为4类，由高到低分别赋值5、4、3、2。由于资源环境承载能力受到压力、状态、响应三层的共同影响，因此对各层结果进行线性加权求和，叠加得到最终的资源环境承载能力得分。

最终，为了表征自然环境适宜性与资源环境承载能力的交互作用，对二者进行乘积运算，得到最终的国土空间开发适宜性得分。

1.3 研究数据

本文从自然环境适宜性、资源环境承载力两个方面构建评价指标体系，研究数据包括第二次土地调查数据、DEM数据、江苏省土地限制性评价数据、经济社会统计数据、岸线底质类型数据。

江苏省第二次土地调查数据包含江苏省各市各区县土地利用现状图斑，及地类名称、用地分类、规划用地等属性数据。

江苏省DEM栅格数据，该数据栅格大小为25m× 25m，通过ArcGIS10.3软件处理获得江苏省坡度栅格数据。

江苏省岸线底质类型数据，包含基岩岸线、人工岸线、砂质岸线、淤泥质岸线、河口岸线、具有自然海岸形态和生态功能岸线长度、分布等属性。岸线侵蚀数据来源于《2018年江苏省海洋灾害公报》，岸线侵蚀区根据年侵蚀速率建立1000m缓冲区。根据《海洋生态保护红线划定技术指南》，重要砂质岸线保护范围为以砂质岸滩高潮线至向陆一侧500m或第一个永久性构筑物或防护林，考虑到河口岸线、具有自然海岸形态和生态功能岸线的生态保护价值，将砂质岸线、河口岸线和具有自然海岸形态和生态功能岸线建立500m缓冲区作为岸线保护区。

地质灾害数据，包含滑坡泥石流灾害、地面沉降、活动断层、岩溶塌陷、地质遗迹、崩塌灾害等属性数据，来源于2017年江苏省土地限制性评价。

经济社会发展数据主要包含国土开发强度（建设用地面积占比）、年末常住人口、地区生产总值、人口密度、城镇化率（非农人口占比）、第二产业比重、第三产业比重、进出口总额、城镇居民人均可支配收入、人均可利用水资源量、实际使用外资额、城市绿化覆盖率、森林覆盖率等。数据均来自2019年南通、连云港、盐城统计年鉴和2018年南通、连云港、盐城水资源公报。

环境数据主要包括环境压力和环境治理两个方面，环境压力数据为2018年化学需氧量、氨氮、总氮、工业废水、二氧化硫、氮氧化物、烟粉尘排放量，一般工业固体废物产生量；环境治理数据包括工业废水处理量、废气治理设施处理能力、一般工业固体废物综合利用率。数据来源于2018年江苏省环境排放统计数据。

2 研究结果分析 2.1 自然环境适宜性整体较好

江苏沿海地区自然环境适宜性适宜、较适宜、一般适宜、较不适宜、不适宜区域面积分别为6734.45 km²、7286.70 km²、2241.42 km²、278.90 km²、14 914.48 km²，分别占江苏沿海地区土地总面积的21.41%、23.16%、7.13%、0.89%、47.41%（表 3）。整体而言，江苏沿海地区自然环境条件较好，适宜及较适宜区域占比44.57%，不适宜区域主要受到基本农田和生态红线的限制。适宜区域主要集中连片分布在连云港沿海、盐城南部和南通中部地区。

2.2 资源环境承载力呈现“南高北低”特征

该部分对压力层、状态层、响应层分别进行主成分分析，以各因子旋转后方差贡献率为权重计算得分，然后利用聚类分析将结果划分成4类，由高到低分别赋值5、4、3、2。

通过主成分分析法提取出三个主成分，对方差的累计贡献率为85.195%。由压力层旋转成分矩阵可以看出，万元GDP的工业废水排放量、二氧化硫排放量以及一般工业固体废物产生量对评价结果产生的影响较大（表 4）。

通过主成分分析法对状态层各因子进行因子分析，三个主成分对方差的累计贡献率84.411%。状态层各项指标中人均GDP、进出口总额、人均消费支出、城镇居民人均可支配收入对状态层得分结果起着较大的影响作用，而二、三产业比重对结果的影响较为轻微，该结果的可能原因是江苏省沿海区域的经济结构差异不大，而经济总量、开放程度和人均收入及消费水平存在一定程度的差异，对结果产生较为显著的影响（表 5）。

通过主成分分析法对响应层各因子进行因子分析，三个主成分对方差的累计贡献率为73.736%。根据响应层旋转成分矩阵可以看出，江苏沿海地区工业废水处理率、森林覆盖率对得分结果影响较显著，这两项指标分别表征了人类社会和自然环境对污染物的转化吸收能力（表 6）。

根据以上分析，江苏沿海地区各区县资源环境承载力评价结果如表 7所示。

2.3 国土空间开发适宜性评价

江苏沿海地区自然环境条件较好，国土空间开发受到地质灾害、海岸侵蚀等自然环境因素影响较小，不适宜区域主要受到基本农田和生态红线的限制。适宜区域主要集中连片分布在连云港沿海、盐城南部和南通中部地区。

资源环境承载力方面，江苏沿海城市整体呈现出南高北低的格局。南通市资源环境压力较江苏沿海其他地区而言较大，盐城市整体压力情况较好，但是滨海县受各类污染物排放量影响资源环境压力在沿海各区县中最大。资源环境承载力状态层评价结果表现出明显的南高北低、核心高外围低的空间分布特征，与经济发展现状呈现良好的对应关系。响应层方面江苏沿海整体呈现南高北低的空间特征，海州区由于工业废水处理率、固体废物综合利用率以及森林覆盖率较高，其应对环境问题能力较强，响应层得分较高。

国土空间开发适宜性评价对自然环境适宜性和资源环境承载力两个维度的评价结果进行了综合集成，能够体现研究区域在各方面因素限制和支撑作用下适宜开发的国土空间。对自然环境适宜性、资源环境承载力得分进行乘积处理，扣除水域部分，利用自然断点法将得分结果分为4类，由高到低分别为适宜、基本适宜、基本不适宜、不适宜。如表 8所示，江苏沿海地区不适宜、基本不适宜、基本适宜、适宜区域面积分别为14 710.37 km²、3376.18 km²、6776.30 km²、5695.64 km²，分别占土地总面积的46.76%、10.73%、21.54%、18.10%。整体而言，不适宜区域与基本适宜、适宜区域占比均较大，即保护空间与发展空间均较为突出。在长三角区域各城市中，南通市适宜开发区域面积最大且占比最多，在沿海地区城市中南通市国土空间开发建设适宜性最好，盐城市其次，连云港市开发适宜性较差，其主要原因一方面是基本农田保护需要限制了国土空间开发建设活动，另一方面则是江苏沿海区域资源环境承载能力呈现为南高北低的空间格局，对国土空间开发适宜性评价结果造成影响。

3 结论与讨论 3.1 结论

（1）江苏沿海地区自然环境条件较好。江苏沿海地区自然环境适宜性适宜、较适宜、一般适宜、较不适宜、不适宜区域面积分别为6734.45 km²、7286.70 km²、2241.42 km²、278.90 km²、14 914.48 km²，分别占江苏沿海地区土地总面积的21.41%、23.16%、7.13%、0.89%、47.41%。不适宜区域主要受到基本农田和生态红线的限制。适宜区域主要集中连片分布在连云港沿海、盐城南部和南通中部地区。

（2）限制江苏沿海城市国土空间开发适宜性的因素主要为资源环境承载能力，资源环境承载力整体呈现南高北低的空间分布特征。沿海地区25个区县中资源环境承载力好、较好、较差、差等级区县个数分别为5、12、7、1个，整体承载能力较好。

（3）国土空间开发适宜性空间分异明显。江苏沿海地区不适宜、基本不适宜、基本适宜、适宜区域面积分别为14 710.37 km²、3376.18 km²、6776.30 km²、5695.64 km²，分别占土地总面积的46.76%、10.73%、21.54%、18.10%。土地空间开发适宜区域面积及比重由高到低分别为南通市、盐城市、连云港市，呈现出由南向北逐渐降低的特征。

3.2 政策建议

（1）为提高江苏沿海城市国土空间开发适宜性，需要分区域有针对性地制定方案，提升资源环境承载力。南通市崇川区、通州市可通过对万元GDP污染物排放量提出限制指标，来进一步提高资源环境承载能力，而南通市全域需进一步提高工业废水处理率和森林覆盖率进而提高资源环境承载能力。对于盐城市而言，经济发展水平、开放程度、人均收入消费水平、工业废水处理能力、森林覆盖率的提高对提高其资源环境承载力有较大意义；尤其对于滨海县而言，减少万元GDP的各类污染物排放量是提高其资源环境承载能力的当务之急。连云港市县域提高资源环境承载力的主要聚焦点应放在提高经济发展水平、开放程度、人均收入消费水平、工业废水处理能力及森林覆盖率。

（2）国土空间开发利用方面，江苏沿海地区需要以提高土地利用效率、盘活存量资源为主要途径，国土空间开发应合理布局，充分考虑适宜性。由于江苏沿海地区单位建设用地经济产出较低，而国土空间开发强度较高，因此土地利用效率需要进一步提升以缓解土地利用与经济发展之间的不匹配问题。

3.3 讨论

本文结合空间数据、统计数据，结合自然环境、社会经济等多方面因素，对江苏沿海地区进行开发适宜性评价。但是仍有三点问题需进一步深入讨论研究：第一，多要素叠加分析应用广泛，但是因子的选择、权重的确定、等级的划分对评价结果具有十分重要的影响，难免存在一定的主观意志。第二，海域数据难以统计、获取，导致研究中对海域因素考虑可能不够全面，在后续研究中可增加海域数据的收集并纳入评价体系。第三，人类社会与海洋的交互作用在沿海城市十分重要，如何通过评价方法的选择、指标体系的构建，更加科学合理地体现该交互作用还有待进一步思考完善。