厦门湾红树林生境适宜性评估及修复潜力分析
2021-03-02王玉玉俞炜炜马志远晁碧霄陈光程胡文佳
张 典,王玉玉,俞炜炜,马志远,晁碧霄,陈光程,陈 彬,胡文佳*
(1.福建省海洋生态保护与修复重点实验室,福建 厦门 361005; 2.自然资源部第三海洋研究所,福建 厦门 361005;3.北京林业大学生态与自然保护学院,北京 100083)
红树林是生长在热带、亚热带潮间带的木本植物群落,为人类社会、沿岸和海洋生态系统提供了多样的生态系统服务[1]。我国红树林主要分布于南部和东南部沿海滩涂[2],然而由于种种因素的影响,过去50 a间红树林面积大幅减少[3-4]。我国从90年代开始致力于红树林的生态保护和修复,现阶段红树林的保护与修复更是受到空前的关注和重视,彻底扭转了红树林面积持续下降的趋势[5-7]。
在红树林修复的过程中,适宜的环境条件和恰当的宜林地选址是保证修复效率的关键[8]。国内外众多学者已经对红树林宜林环境指标开展了大量研究[9-13],但是利用模型手段定量探索多种环境因子交互作用下对红树林分布影响的研究还较为少见。此外,我国造林设计的基础科学依据较为薄弱,对宜林滩涂的潜力分析不充分是导致部分区域红树林修复成效不理想的主要原因之一[14]。目前传统上红树林修复区及其面积的确定多基于单项或多项宜林条件的经验估算[11,15-16],难以在空间上实现精细化的定量预测。因此,深入分析区域尺度上红树林与复杂环境因子间的定量关系,从区域角度科学评估和确定优先修复地点,准确分析红树林生境适宜性和修复潜力,对红树林生态修复的规划与选址具有重要的科学意义。
早期研究物种分布与环境间的定量关系往往依据经验进行判断,存在一定的主观性。随着研究手段的多样化和科技水平的提升,利用物种分布模型(Species Distribution Models, SDMs)对物种的适宜分布空间进行评价成为了近年来的研究热点[17]。SDMs通过总结归纳物种出现、潜在生境要求与环境变量之间的关系,以此预测物种在分布区域出现或聚集的概率,能够利用不完整的信息对物种的潜在栖息地和分布范围进行较为全面的预测[18]。其中最大熵模型(Maximum Entropy, MaxEnt)由于预测效果优异、数据要求低等特点,被广泛应用于物种保护区规划[19-20]、物种空间分布对气候变化的响应[21-22]、珍稀动植物和入侵物种的潜在分布预测[23-25]等方面。该模型基于最大熵理论,在已知样本点和对应环境变量的基础上,通过拟合具有熵值最大的概率分布对物种的潜在生境做出预测[26]。近年来,MaxEnt模型的研究对象已扩展到湿地或潮间带植被,例如入侵物种互花米草(Spartinaalterniflora)[27]以及其他滨海湿地植物[28],但目前有关红树林的研究案例还十分罕见[12]。
基于此,本研究以福建省厦门湾作为研究区,利用MaxEnt模型开展区域尺度上的红树林生境适宜性评估,并叠加LULC数据中不同地类的可修复性,从而估算区域红树林修复潜力。厦门湾作为福建省重要海湾,接近中国大陆红树林分布的最北界,是福建省红树林自然分布的典型地区,也是红树林生态修复的热点区域[16,32]。通过本研究可判断厦门湾红树林生境适宜性的分布格局,识别潜在修复区的地理位置,定量估算潜在修复区的面积。研究结果可为厦门湾红树林的保护和修复行动提供科学基础,并为其他区域或类似尺度上的红树林修复空间规划提供可借鉴的方法。
1 数据与方法
1.1 研究区域概况
本研究以福建省厦门湾作为研究区,包括晋江市围头角至龙海市镇海角连线以西、九龙江口紫泥镇以东海域及周边陆地,总面积约为1 082 km2。厦门湾红树林主要分布于九龙江口浮宫片区(海门岛+浮宫)、海城片区(东园+海城+玉枕洲+大涂州)、甘文片区(紫泥+甘文农场)以及厦门市思明、湖里、集美、海沧和翔安5个区,总面积约为373 hm2[33-34]。在20世纪70—80年代由于厦门湾内的伐木种草、围塘养殖和围海造田等活动曾导致红树林面积急剧减少,随着1988年九龙江口红树林省级自然保护区的成立以及1997年后红树林修复工程的兴起,红树林面积呈现出持续增长的趋势[32,35-36]。考虑到红树林的生态特点,本研究沿海岸线分别向陆方向和向海方向设置各2 km距离的缓冲区,作为模型的建模评估范围,大致包括向陆高程7 m以内、向海水深6 m以内所有区域,完全覆盖厦门湾现有红树林分布范围及其毗邻区域(图1)。
图1 厦门湾研究区Fig. 1 Study area in Xiamen Bay图中灰色(图5、7蓝色)部分为厦门湾海域范围,下同。
1.2 数据来源及处理
1.2.1 红树林分布数据的获取及处理 基于2018年福建省海洋生态红线数据和2017年厦门湾红树林的调查站点数据确定了湾内红树林主要分布区[33-34],利用2018年ESRI World Imagery遥感影像对主要分布区内的红树林小斑进行目视解译得到现存红树林小斑共131个,面积总计328.67 hm2,并在2018—2019年通过现场勘查进行验证,表明目视解译提取的红树林小斑精度超过85%。在ArcGIS10.5软件内采用渔网格工具对识别出的红树林小斑进行重采样,采样精度为50 m×50 m,最终获得1 348个厦门湾红树林“出现点”数据,将所有出现点的经纬度坐标在Excel里导出为csv格式作为MaxEnt模型的生物分布输入数据(图2)。
图2 2018年厦门湾红树林分布概况Fig. 2 Distribution of mangroves in Xiamen Bay in 2018由于尺度原因,厦门岛内五缘湾及筼筜湖附近的水域未在图上显示。
1.2.2 环境数据的获取及处理 温度、盐度、地形、潮汐、水质、潮滩沉积物以及风浪作用等是影响红树林生长分布的主要环境因素[37-38]。本研究根据红树林适生环境条件特点,兼顾厦门湾的评价尺度、环境特征和数据的可获得性,筛选出生物气候变量、地形、水质、沉积物、水文等5组共36个环境变量用于模型建立(表1)。
表1 厦门湾红树林MaxEnt模型输入的环境变量列表
利用ArcGIS10.5软件的空间分析模块计算以上环境变量数据,同时采用软件中的“扩展”分析工具将各类环境变量填充至研究区建模范围,所有环境变量数据均重采样至250 m×250 m的分辨率下,并在ArcGIS10.5软件内转换为整形的ASCII格式,作为输入MaxEnt模型的环境变量数据集。
在红树林修复潜力计算中需采用LULC数据,本研究的LULC数据来源于国家基础地理信息中心全国地理信息资源目录服务系统(http://www.webmap.cn/),数据精度为30 m,在ArcGIS10.5软件中裁切至研究区范围以便开展后续运算。
1.3 MaxEnt模型建立
1.3.1 模型原理及参数设置 MaxEnt模型基于最大熵原理,即在充分考虑已知信息而不对未知情况做任何主观假设的前提下预测一个随机事件的概率分布,此时信息熵最大的分布状态为最接近真实的状态[26]。最大熵问题是一个带约束的最优化问题,当已知x输出为y时,对于给定的训练数据集以及特征函数fi(x,y),i=1, 2, ...,n,最大熵模型求解方程为:
对各种地理水纹记号进行汇总、划分、归类后,引入视觉传达艺术的设计方法,并融入我国传统水文化中的创意元素,进行图形设计,填补地理水纹记号的空缺;然后依据我国水利信息化的要求,进行数字化处理;最后再对整理后的结果进行处理,向水利信息化靠拢,最终形成一套完整的系统。
∑yP(y|x)=1
(1)
式(1)中:H(P)为条件熵,P(y|x)为条件概率分布假设,log取自然对数,使用拉格朗日乘子法,将带约束的最优化之原始问题转换为无约束的最优化之对偶问题进行求解[40]。
本研究使用MaxEnt 3.4.1 版本软件开展评估研究,MaxEnt模型参数设置选择75%的厦门湾红树林分布数据作为训练数据用于模型建立,剩下的25%作为测试数据用于模型检验[25]。底质类型和流速数据选择分类型(Categorical),其余环境变量皆为连续型(Continuous)。在模型建立时选择默认特征组合,取样方法为随机抽取,模型重复计算次数设置为10次,其它设置选项保持默认。模型输出结果中,每个栅格的像元值代表该栅格内红树林的分布概率,其像元值的取值范围在0~1区间,像元值越大表示该栅格内红树林的适生程度越高。栅格数据结果在ArcGIS 10.5内进行后续分析,采用自然断点法对栅格预测结果分类,获得生境适宜性等级划分的界值。
1.3.2 模型结果检验 模型预测一般会产生过高估计(假阳性)和过低估计(假阴性)两类错误[17],本研究同时选择AUC值、Kappa系数和真实技巧统计值(True Skill Statistic, TSS)3种统计方法,能够较为全面地评估模型的准确性[22]。
AUC值通过模型中自带的ROC (Receiver Operating Characteristic)曲线分析方法得到,其曲线以假阳性率为横坐标,真阳性率为纵坐标,曲线下的面积即为AUC值[41]。取值范围在0.50到1.00之间,当AUC值大于0.75时,模型被认为可用[42]。Kappa和TSS值通过3.5.3版本的R语言进行计算,其中Kappa值表示的是预测正确的样本个数与预测错误的样本个数之间的比值,TSS值则为实测样本的净预测成功率, Kappa和TSS值的取值范围都在-1.0~+1.0之间,其值越接近+1.0表明模型预测效果越好,其值小于等于0时则表明预测结果差于随机预测结果[43-44]。
1.4 红树林修复潜力计算方法
在红树林的潜在适宜分布区中,一些用地类型事实上无法用于开展生态修复(如人造地表),而另外一些用地类型具有很高的修复潜力(如湿地),此外部分人工用地类型如养殖虾塘等目前被认为可以通过退塘还林修复为红树林生境[45]。因此,在红树林潜在适生区范围内我们可以根据用地类型判断各栅格内实际开展红树林修复的可能性。本研究将生境适宜性评价结果与用地类型的修复可能性结合,建立了红树林修复潜力计算方法。首先通过德尔菲法依据专家经验判断不同用地类型可修复/转换为红树林林地的可能性大小,对每种用地类型设置转换系数(表2),在ArcGIS10.5软件中利用“重分类”工具将转换系数值赋予各用地类型。再将生境适宜性评价结果与转换系数相乘,从而得到每个栅格内的红树林修复潜力估算值,估算值越大代表修复潜力越高。最后将潜在修复区内的现存红树林面积扣除,即可定量估算区域内理论上可用于红树林修复的用地面积。通过叠加厦门湾红树林相关保护区和生态修复工程的位置[34],可进一步分析现有保护与修复的空间格局并提出对策和建议。
表2 厦门湾LULC-红树林转换系数赋值
2 结果与讨论
2.1 厦门湾红树林分布与环境因子的关系
2.1.1 环境变量贡献分析 环境变量贡献度计算结果表明,对红树林分布影响最重要的环境因子类型包括温度(24.6%)、地形(27.2%)和水质(27.1%)。刀切图的结果(图3)显示,在单项环境变量中,离岸线欧式距离、盐度、无机氮浓度、活性磷酸盐浓度、镉浓度、最湿月降水量、冬季海表温度等环境变量作为唯一变量参与模型时,对模型提供的信息最多。其中离岸线的欧式距离包含了更多与其他环境变量无关的独立信息,如果剔除该变量,则会对模型结果产生较大的影响。
图3 基于刀切法的模型正则化训练增益Fig. 3 Training gain of model regularization based on jackknife
2.1.2 红树林适生环境因子响应曲线分析 响应曲线的分析结果表明,温度、降水、水质等环境变量与红树林分布概率多呈单峰模式,将存在概率大于等于0.5以上的环境变量范围作为厦门湾红树林分布的最适环境变量区间。选取各变量组内的代表性环境因子(最冷季平均温度、冬季海表温度、最湿月降水量、盐度、无机氮浓度、活性磷酸盐浓度、镉浓度、水深),绘制厦门湾红树林分布的环境变量响应关系曲线。结果如图4所示,影响厦门湾红树林分布的代表性环境变量适宜区间为:最冷季度平均温度12.95~13.62 ℃,冬季海表温度17.05~17.44 ℃,最湿月降水量166.57~172.21 mm,盐度6.02~17.62,无机氮浓度1.176~1.541 mg/L,活性磷酸盐浓度0.103~0.154 mg/L,镉浓度0.071~0.185 μg/L,水深1.99~6.67 m。
图4 厦门湾红树林分布对代表性环境因子的响应曲线Fig. 4 Response curve of mangrove distribution to representative environmental factors in Xiamen Bay
2.2 厦门湾红树林生境适宜性评估结果
模型检验结果表明AUC平均值为0.98,最大Kappa值为0.53,TSS平均值为0.63,提示模型预测结果可信。采用自然断点法将红树林按照生境适宜性评估值划分为高适宜区(>0.58)、中适宜区(>0.39~0.58)、低适宜区(>0.24~0.39)和不适宜区(≤0.24)4类。预测得到厦门湾高适宜区面积约为1 525.0 hm2,中适宜区和低适宜区面积分别为2 318.8 hm2和3 412.5 hm2,高适宜区主要位于九龙江河口、同安湾下潭尾、海沧东屿、厦门岛西南部、同安丙洲岛、同安官浔溪和漳州开发区沿岸滩涂(图5)。各行政区之间比较,红树林潜在适宜生境分布面积最大的分别为漳州龙海市、厦门翔安区、厦门海沧区和厦门同安区(图6)。
图5 厦门湾红树林生境适宜性分布图Fig. 5 Map for habitat suitability of mangroves in Xiamen Bay1.海沧沿岸,2.下潭尾,3.大涂洲。
图6 厦门湾各行政市(区)红树林潜在适宜生境面积比较Fig. 6 Comparison of potential suitable habitats for mangroves in sub-region of Xiamen Bay
2.3 厦门湾红树林修复潜力预测结果
根据修复潜力计算结果, 将预测得到的红树林潜在修复进一步区分为高修复潜力区(>0.50)、中修复潜力区(>0.32~0.50)、低修复潜力区(>0.19~0.32)和不可修复区(≤0.19)4类,其中高修复潜力区面积达571.0 hm2,中修复潜力区面积达1 053.3 hm2,低修复潜力区面积2 470.5 hm2。理论潜在修复可能性较高的区域主要位于九龙江河口、同安湾下潭尾、海沧南部沿岸、海沧东屿沿岸滩涂、同安丙洲岛和官浔溪(图7)。统计高修复潜力区中的用地类型占比,可见水体和湿地这两种用地类型的比例最高,分别达到46.1%和33.2%。
图7 厦门湾红树林潜在修复区分布Fig. 7 Potential areas for mangrove restoration in Xiamen Bay1.海沧沿岸,2.下潭尾,3.大涂洲。
从各行政区的潜在修复区分布来看,潜在修复区集中在龙海市、翔安区、海沧区和同安区,其中高修复潜力区主要分布于漳州龙海市(95.86%)、厦门翔安区(2.10%)、厦门同安区(0.71%)和厦门海沧区(0.58%)。对比潜在适宜生境,厦门湾潜在修复区的总面积减少了43.2%,意味着有接近一半的潜在适生区因为人类活动开发占据而造成永久或半永久损失,无法再用于红树林的修复。厦门湾内各行政区的潜在适宜生境损失率各不相同,厦门市的潜在适宜生境损失程度明显高于其他两个地级市,达到58.10%,其中最高损失率出现在厦门市同安区(95.68%)和海沧区(92.39%)。扣除本研究解译的现存红树林分布,在厦门湾内还有约406.4 hm2的高修复潜力区和1 001.2 hm2的中修复潜力区,可以考虑用于未来红树林修复的选址,主要分布于九龙江口和下潭尾现存红树林附近、海沧区南部沿岸、火烧屿南部、安海湾南部沿岸、丙洲岛—官浔溪、厦门岛西南部等地(表3)。
表3 厦门湾现存红树林和理论可修复面积对比
2.4 厦门湾红树林适生条件分析
红树林的生长分布与环境关系密切,其中海滩的地理条件和气候因素是决定红树林引种是否成功的最重要因素[46]。本研究MaxEnt模型输出结果显示温度和地形因子对于红树林分布的影响最大,除此之外,我们还发现在小区域尺度上水质因子也是影响红树林分布的一个重要因素,这一点与Abd-El Monsef 等(2017)的研究[11]一致。
温度是控制红树林分布的主要因素之一。厦门湾现有真红树共7种,其中分布最多的为本土种的秋茄(Kandeliacandel)和白骨壤(Avicenniamarina)[47]。按照红树植物对气温适应范围的分类,秋茄和白骨壤都属抗低温广布种[9]。相关实验研究结果表明抗低温广布种的叶片半致死温度低于-6 ℃,生物学临界温度在10 ℃左右[46,48]。MaxEnt模型输出结果表明厦门湾红树林对最冷季平均温度的适生区间为12.95~13.62 ℃,当气温低于或超过阈值范围时,红树林存在概率会迅速下降。其结果可以与传统生理学研究的结果互相参照。
根据模型输出结果,厦门红树林水深的适生阈值在1.99~6.67 m,最适峰值为4.62 m,该适生阈值可作为造林选址的参考。在以往的研究中,学者们往往以野外试种或实验室模拟的方式来探索红树林的滩面高程和宜林水位线[49-52]。然而受到各地不同潮汐特点的制约,国内对红树林宜林水位线的划分和定义标准不尽相同[49,52]。张乔民等(2001)建议利用潮汐浸淹频率代替特征潮位更能直观反应潮汐变化对红树林的影响[9],然而其浸淹频率难以在全湾的研究区域上以较为精确的数值形式表示。本研究选用海图水深数据,以理论最低潮面为基准面,可以修正测得水深中的潮高,更便于不同研究区域之间的比较。
红树林对于水体中营养盐和重金属的净化容量是有一定限度的,在一般的研究中并未将营养盐或重金属等环境变量作为影响红树林的分布的限制因子[37]。然而本研究发现在单项环境变量中,营养盐和重金属浓度均会对红树林分布产生一定的影响。本研究模型输出结果表明红树林最适生长的活性磷酸盐浓度范围为0.103~0.154 mg/L,无机氮浓度为1.176~1.541 mg/L,镉浓度为0.071~0.185 μg/L。在各种重金属离子中,铅、镉等重金属离子容易富集在红树植物的根部相对难以迁移,其中秋茄幼苗对土壤中镉的耐受临界值为30 μg/g[53]。红树林对氮、磷等营养盐的富集效率受到树种的影响,同一条件下各红树树种对营养盐的去除效率由高到低为海桑(Sonneratiacaseolaris)>桐花树(Aegicerascorniculatum)>木榄(Bruguieragymnorrhiza)>秋茄[54-55]。根据本研究结果,红树林适生的营养盐浓度较高,这可能与红树林适生区域集中在河口等营养盐浓度较高的区域有关,也从侧面反映了红树林对营养盐的耐受程度较高。因此,我们建议在实际的保护和修复过程中要对红树林分布区的水质环境进行有效监控,对于天然红树林的污水净化开发利用持保守态度,对人工林则要根据其树种类型和规模充分考虑其净化容量,提高净化效率。
2.5 研究方法的适用性分析
本研究应用MaxEnt模型探索区域尺度上的红树林潜在适宜生境和修复区的分布,取得了良好的模拟效果。厦门湾内红树林的主要分布区或修复区如九龙江河口、下潭尾湿地公园、海沧沿岸等区域均在模型输出结果的高适宜区域内。除此之外,在本研究数据集中未记录红树林分布点的一些区域亦预测了较高的适生概率,例如翔安区九溪河口、安海湾沿岸海域、漳州开发区沿岸等区域,经文献资料的查阅和对比,这些区域均现存有零星红树林或是曾有红树林分布[47,56]。
有学者曾采用传统的经验评估方法预测厦门的潜在宜林地,林鹏等(2005)曾在红树林宜林地现状调查的基础上,预估厦门市红树林宜林滩涂面积在108~198.4 hm2[16];胡杰(2016)通过文献查阅和实地考察,结合GIS技术评估出厦门市秋茄和白骨壤宜林地的面积分别为216.12 hm2和259.34 hm2[15]。根据本研究的结果,如将预测得到的高修复潜力区和中修复潜力区纳入红树林潜在修复区,则可统计得到厦门市红树林潜在修复区面积约为217.17 hm2,扣除现有红树林面积21.00 hm2,尚有约196.17 hm2的理论修复面积。传统方法评估的结果与本研究利用生物分布模型分析得到的结果较为近似,说明模型结果具有一定的参考价值。
本研究因受限于研究尺度和数据分辨率,影响了厦门岛沿岸五缘湾、筼筜湖等部分区域的模拟精度,在这些区域内未能清晰识别红树林适生地点。在今后的研究中可考虑提升水文、地形、沉积物等数据的质量和精度,从而进一步提升模型的整体性能。另外就数据获取时间而言,为了保证数据的一致性,水质等数据只收集到大潮时期的数据,且多集中于春夏季节,如果能有更多时期的数据来反映其季节性的变化,对提升模型准确度可能会更有帮助。此外,本研究中的厦门湾红树林修复潜力主要考虑是自然适生条件与土地利用空间格局之间的关系,未能包括物种入侵、航运、污染等所有可能影响红树林修复的干扰因子,在未来的研究中可进一步补充。
2.6 厦门湾红树林保护与修复潜力分析与对策
厦门湾内共有龙海市九龙江口红树林省级自然保护区和福建厦门珍稀海洋物种国家级自然保护区两处红树林相关的保护区,面积总计336.79 km2,已基本形成了有效的红树林保护网络,覆盖了相当一部分潜在适生区。厦门湾内已开展的生态修复工程包括筼筜湖生态修复工程、下潭尾滨海湿地公园、五缘湾湿地公园、海沧湾岸线整治工程、环东海域红树林绿化项目、安海湾环境综合整治工程、角屿红树林生态园等六处,大部分落在本研究预测得到的理论潜在修复区内。在现有的保护区和修复工程之外仍存在一些理论潜在修复区,建议未来厦门湾内的红树林修复选址可考虑泉州市安海湾沿岸、龙海市九龙江河口、漳州港沿岸和附近岛屿。在开展修复前需对拟选址区域的环境条件进行充分评估,以提高红树林造林的成功率。
有研究表明1980年以来97.6%的中国红树林损失可归因于历史形成的围塘养殖[57],相当一部分养殖鱼塘具备恢复红树林的基本生态条件,因此有学者提出退塘还林应成为中国红树林修复的主要来源[5,7]。本研究表明厦门湾红树林理论可修复区域的用地类型主要是水体和湿地,与遥感影像的比较发现水体类型内包括了大量的养殖坑塘,从侧面印证了养殖虾塘等存在较高的修复可能性。此外,我们发现红树林潜在适宜生境的损失比例与区域经济发展状况密切相关。各市损失比例为厦门(58.1%)>泉州(38.2%)>漳州(34.0%),厦门市的比例远远高于漳州市和泉州市。厦门市是三市中城镇化率和人均地区生产总值最高的地区[58-60],根据近60 a来厦门湾围填海区域的变化趋势,厦门市的围填海强度远远高于漳州市和泉州市[61],我们推测这可能是人类活动强度影响红树林修复潜力的内在因素。根据本研究的结果,厦门市各行政区内有大量红树林潜在适宜生境已被其他用途占据,相较而言漳州市和泉州市可能具备更好的红树林修复前景。
3 结论
(1)在厦门湾研究区域内,红树林的分布主要受到温度、地形和水质的影响,其中单项环境变量的最优适生条件为:最冷季度平均温度12.95~13.62 ℃,冬季海表温度17.05~17.44 ℃,最湿月降水量166.57 ~172.21 mm,盐度6.02~17.62,无机氮浓度1.176~1.541 mg/L,活性磷酸盐浓度0.103~0.154 mg/L,镉浓度0.071~0.185 μg/L,水深1.99~6.67 m。
(2)生境适宜性较高的区域主要分布于九龙江河口、翔安区下潭尾、海沧区东屿、厦门岛西南部、同安区丙洲岛和官浔溪、漳州开发区沿岸滩涂,高适宜区面积达1 525.0 hm2。高修复潜力区主要分布于九龙江河口、翔安区的下潭尾、海沧区东屿和南部沿岸、同安区丙洲岛和鳄鱼屿、翔安区九溪河口、安海湾、漳州开发区沿岸滩涂,高修复潜力区面积达571.0 hm2。与现存红树林分布相对照,理论上还有约406.4 hm2的高修复潜力区和1 001.2 hm2的中修复潜力区可能用于红树林生态修复。
(3)在厦门湾内约有43.2%的潜在红树林适生区因为人类活动开发占据而造成永久或半永久的损失,无法再用于红树林的修复,其损失率与社会经济发展程度密切相关。厦门湾内现有的保护区和修复工程基本覆盖了预测得到的潜在修复区,未来厦门湾内的红树林修复选址可考虑泉州市安海湾沿岸、龙海市九龙江河口、漳州港沿岸和附近岛屿。