珊瑚礁是由珊瑚生物群落(珊瑚虫)聚集的石灰石保护结构形成的, 这些生活在热带水域的群体作为一个统一的有机体, 通常接近地表并生长了成百上千年。珊瑚礁生态系统是生物多样性和生态价值最高的海洋生态系统之一, 它古老、珍贵又多姿多彩, 常被誉为“海洋中的热带雨林”或“海上长城”^[1-2]。在近几十年里, 由于全球气候变化, 海洋酸化和人类的不合理活动^[3], 世界上许多地方珊瑚礁的退化速度远超其生长速度, 并呈加速退化的趋势发展^[4]。如我国南海的珊瑚礁面积约有37200 km², 占世界珊瑚礁面积的5%, 居东南亚地区第二位^[5-9], 但南海珊瑚礁生物资源呈现出十分严重的衰退趋势, 珊瑚礁损害率甚至高达90%以上^[10]。南海诸岛珊瑚礁不仅具有渔业、旅游业和海岸保护等价值和服务, 还因其得天独厚的地理位置而对领土、领海和海洋专属经济区的划分起着无可替代的作用^[11], 于是保护和修复珊瑚礁生态系统无疑成为当今国际生态保护的研究热点。

为探求珊瑚礁生态系统各部分的运行机制和相互作用, 借助模型手段进行研究的方法颇为盛行。珊瑚礁生态系统模型研究是根据生态系统内生物及其周围环境的特性, 通过计算机等手段构建出模型, 对珊瑚礁、珊瑚礁生物种群群落及生存环境进行追踪和检测, 一定程度上有助于了解和掌握系统动态, 可为保护和修复珊瑚礁生态系统提供可行的理论与实践指导。国际上对珊瑚礁生态系统建模研究起始于20世纪70年代, 随着珊瑚礁退化造成的负面影响日益显著, 系统建模研究发展迅速。Macintyre等^[12]为激发对珊瑚礁地质调查中进行模拟的有效性讨论, 于1974年首次采用图形模型描绘了珊瑚礁生态系统中碳流动模式的位置和特征, 为未来珊瑚礁生态系统的模拟分析研究奠定了基础。Lidz等^[13]构建了迎风边缘模型, 为海平面波动期间碳酸盐岩台地的进积提供了一个重要而独特的运行机制, 为背礁露台盆地的填充做出重要贡献。van der Laan和Bradbury^[14]使用一个简单的人工渔礁系统, 通过元胞自动机实现了模拟荆棘冠海星对珊瑚产生破坏的现象, 成功研究了二者捕食与幼虫扩散和生存的影响关系。在珊瑚礁礁体建模研究的基础上, 近些年来珊瑚礁生物如珊瑚共生体、珊瑚虫、珊瑚鱼等珊瑚物种的模型日益增多, Mcclanahan^[15]、Jokiel, Martinelli^[16]、Andres和Rodenhouse^[17]等国际学者分别对珊瑚生物的物种分布、种间关系、种群大小及生长速度等进行监测和评估, 旨在保护和提高珊瑚礁生态系统的生物多样性。国内珊瑚礁生态系统建模研究起步较晚, 珊瑚礁礁体与波浪变形传播^[18-21]等方面深受学者关注, 投入使用的数理模型相对较多。

在珊瑚礁礁体模型、种群模型和生态系统模型研究中, 各种建模方法与技术复杂多样, 模型模拟结果差异较大, 为了能够准确把握珊瑚礁生态系统模型研究的趋势, 更好地解决珊瑚礁生态系统退化诊断的建模工作, 有必要对当前有关珊瑚礁模型进行归纳分析。文献计量是一种基于数理统计的定量分析方法, 是以科学文献的外在特征和数量为研究对象, 采用数学与统计学方法来研究文献的分布结构、数量关系和变化规律等, 用以描述、评价和预测科学技术的现状与发展趋势^[22-23]。为了解各学科的研究状况、热点及趋势, 文献计量方法已被多学科广泛应用^[24]。因此本文借助文献计量方法为辅助手段, 着重对当前国内外用以探究珊瑚礁生态系统的多种模型研究进行定量的统计分析, 旨在了解掌握目前珊瑚礁生态系统模型研究的热点及发展趋势, 为国内珊瑚礁生态修复的后续发展提供文献参考。

1 数据来源与研究方法 1.1 数据来源

本文以Web of Science数据库中的Web of Science ^TM核心合集为国外文献数据源, 通过数据库检索相关文献。检索方式为“主题检索”, 检索式为“coral reef”and “model”, 文献类型为“article”, 搜索年份截止在2018年10月(包括10月), 共得到4057条符合条件的结果。将结果按照相关性进行排序, 选取前2000个结果作为数据进行热点统计分析。

相比起国外较为丰富的研究现状, 国内对于珊瑚礁生态系统模型构建的研究尚在萌芽之中。在中国期刊全文数据库(CNKI)中, 选择主题为“珊瑚礁”并且“模型”的期刊文献, 搜索年限为至今, 仅得到34篇相关文献。

1.2 研究方法

本文借助文献计量为数理统计手段, 对国内外珊瑚礁生态系统模型研究的现状以及热点态势进行分析, 并进行可视化处理。利用Bibexcel软件不仅可以完成多种计量分析工作且便于与其他分析软件进行数据交换的特点, 将关键词的词频进行统计并生成相关矩阵。具体方法为：将国外文献按照相关性排序后选取前2000篇以及国内34篇题录分别导入Bibexcel软件中, 导出文献的高频关键词。国外研究去除最高频次并筛选出词频大于13的关键词共计57个, 合并其中的同、近义词, 如:“biodiversity”、“diversity”都是多样性的意思, 故统一用“biodiversity”来表示, 合并后共计43个高频关键词。国内研究则选取了9个频次大于等于2的关键词, 利用Bibexcel软件将关键词进行两两配对。统计2000篇国际文献里高频关键词的共现频次形成43×43的共词矩阵, 随后将矩阵导入到Ucinet中生成##h文件, 再用Netdraw做出国际珊瑚礁生态系统模型研究高频关键词共现网络可视图。按照关键词的基本属性对关键词进行聚类分析, 得出系统聚类分析表。

2 结果分析 2.1 发文量及各阶段内容分析

以Web of Science数据库检索的期刊文献数量作为国外数据来源, 日期截止到2018年10月, 共得到4057条结果；在CNKI中检索到34篇国内与珊瑚礁生态系统模型研究有关的文献, 综合得到国际年度发文量如图 1所示。

2.1.1 发文量分析

由图 1可以看出, 就国外发文量而言, 尽管个别年份发文数量有些许下降, 但文献数量总体上呈大幅度增长的趋势, 研究具体可分为五个阶段。第一阶段(1974—1990)萌芽阶段：1974年起, 美国学者发布了首篇利用概念模型模拟珊瑚礁地质环境的文章, 开辟了借助模型探究珊瑚礁生态系统的先河, 随后陆续有学者加入模型研究行列, 由于模型研究正处于探索阶段, 发文量稀少。第二阶段(1991—1997)缓慢发展阶段：随着珊瑚礁模型研究的不断深入探索, 相关研究的发文量较前一阶段有明显增多, 但增长速率较低, 且呈现出停滞状态。发文量最多的是美国, 发文量64篇, 其次是澳大利亚、加拿大, 发文量分别为39篇和15篇。第三阶段(1998—2005)稳定发展阶段：随着国际社会对珊瑚礁生态系统的聚焦, 有关珊瑚礁生态系统珊瑚共生体及生境模拟模型研究明显增多, 其中发文量较多的国家有美国、澳大利亚和法国。第四阶段(2006—2016)波动增长阶段：由于珊瑚白化退化现象日益严重, 珊瑚礁生态系统的退化引起了各国学者的强烈关注, 大量相关模型研究涌现, 其中美国的发文量仍处于领先地位, 十年间发文量已高达1211篇占国际总发文量的44%, 其次是英国和澳大利亚。第五阶段(2016年至今)加速增长阶段：2016年至今虽发文量有下降形势, 但珊瑚礁生态系统的保护和修复已成为全球瞩目的研究领域, 备受各国重视, 致使各国发文量普遍大幅度提高, 该领域的研究进入了加速发展阶段。发文量超过100的国家有四个, 依次是美国、澳大利亚、英国和法国, 美国在加勒比海地区珊瑚礁的研究成果始终遥遥领先。

经过图 1对比可以明显看出, 我国在这方面起步较晚, 对于珊瑚礁生态系统模型的研究相对匮乏。1997年才发表第一篇文章, 同时相关研究的发展不稳定。虽然我国致力于珊瑚礁生态系统的修复工作已持续十余年, 但针对构建模型进行生态系统研究的内容仍然匮乏, 直到2016年以后, 随着科技的发展时代的进步, 研究手段和技术都有了很大提高, 发文量才有所上升, 但仍远远少于国际发文量。2018年我国的发文量达到了最高值, 共有9篇文献, 占国内相关研究总发文量的26.4%。

2.1.2 研究内容分析

总体上看, 国外用于研究珊瑚礁生态系统的模型范围广、用途多, 该领域研究内容丰富, 囊括了珊瑚-珊瑚共生体-生物种群-群落等模型的构建, 用以对珊瑚礁生态系统组成部分进行模拟、监测、追踪、评估和预测, 各阶段致力的研究内容有所不同。萌芽阶段学者主要致力于对珊瑚礁成礁环境的探索, 通过建模对珊瑚礁地质地貌性质、成礁环境进行模拟和研究。缓慢增长阶段在前人研究的基础上, 加大了珊瑚礁成礁环境的探究力度, 模型构建手段日益成熟, 同时珊瑚礁生物群落的组成、生态分布、生物生产力和营养结构等研究增多, 运用数学模拟等手段对群落的形成、演化和动力过程、变化趋势进行研究成为该领域的主要研究方向。稳定发展阶段该领域有了新的关注点, 此时珊瑚退化现象逐渐涌现, 各国纷纷投入到珊瑚退化的研究中。失去共生藻的珊瑚因缺失营养物质和新陈代谢的能力大范围白化、退化, 为达到通过保护和繁殖共生藻缓解珊瑚白化现象的目的, 模拟藻类生存环境的仿真模型骤增, 用于追踪和监测珊瑚共生藻的模型层出不穷。珊瑚礁生态系统的剧烈恶化使得保护和修复珊瑚礁生态系统成为学者热议话题。由于珊瑚礁生态系统是高度复杂的系统, 其特征在于生物和非生物成分之间的非线性关系, 研究人员将研究重点转移到珊瑚礁生物和非生物成分的保护上。珊瑚、珊瑚共生体、珊瑚礁鱼类、底栖环境等模型模拟研究成果丰硕, 珊瑚礁的生物群落物种复杂度及生物多样性的研究成为该领域的新焦点。由于世界各地的珊瑚礁都承受着由异常的温水事件带来巨大压力, 导致大规模珊瑚褪色和大范围的珊瑚礁破坏, 所以近阶段以来, 全球气候变化、海水温度变化对珊瑚礁生态系统的影响机制研究明显增多, 用于预测和判断未来气候变化会给珊瑚礁生态系统带来的胁迫和风险的模型逐渐增多, 珊瑚礁生态系统风险、威胁、胁迫和修复的研究成为该领域的主要发展趋势。

2.2 高频关键词分析

关键词是对文章的高度概括和精炼, 将筛选出的国际2000篇和国内34篇相关文献利用Bibexcel软件导出高频关键词, 由于关键词较多, 为保证研究的共性, 国外分析选择出词频大于13的关键词共计57个, 合并其中的同、近义词, 得到共计43个高频关键词, 国内分析则选取频次大于2的共计9个。这些高频关键词在一定程度上分别反映了国内外珊瑚礁生态系统仿真模拟研究的核心和热点。利用Bibexcel将从Web of Science数据库中提取高频关键词结果如表 1所示。

根据表 1可以看出：国外关于“珊瑚礁生态模型”研究的文献中, 出现最多的关键词为Coral reef(珊瑚礁), 此外climate change (气候变化)、Coral bleaching(珊瑚白化)、Coral reef fish(珊瑚礁鱼类)、Great Barrier Reef (大堡礁)和biodiversity(生物多样性)五个关键词出现频率较大, 频次均超过50占总统计高频词的25.2%, 说明国外对于珊瑚礁生态系统建立的模型主要集中用于珊瑚白化以及生态系统的物种多样性上。气候变化是珊瑚白化的原因之一, 海水升温会使珊瑚虫释放其体内的虫黄藻。虫黄藻既为珊瑚提供营养, 又给珊瑚带来色彩, 因此虫黄藻被释放后珊瑚就会出现不同程度的“白化”^[25-26]。丰富的珊瑚礁资源为珊瑚礁鱼类提供了栖息及繁育场所, 但珊瑚白化对珊瑚礁鱼类的生境造成改变, 为珊瑚礁生态系统的生物多样性带来不容小觑的影响。因出现频次介于50和20(包括20)之间的关键词共19个, 占总统计高频词的36.3%。根据Connectivity(连通性)、ocean acidification(海洋酸化)、resilience(弹性)、calcification(钙化)等关键词可以看出研究人员对于珊瑚礁礁体的结构和组成投入了大量研究, 甚至构建人工礁, 以期通过建立一个稳定、抵抗风浪、提供庇护所的岩体结构, 促使珊瑚快速安全的附着、生长和繁殖^[27-28]。频次低于20但大于13的关键词共有268个, 占全部的16.1%。随着珊瑚礁生态系统修复研究的不断深入, Hydrodynamics(流体力学, 水动力学)、Water quality(水质)、Fisheries(捕鱼业)、Coral Disease(珊瑚疾病)的研究不断增多, 说明用以探究珊瑚礁生物生境特性如营养物质、生存动力的模型愈加丰富, 珊瑚礁生态系统修复研究范围日益广泛。

利用Bibexcel将从CNKI中提取出高频关键词可知, 除最高频率的珊瑚礁以外, 其他关键词出现了2次分别为：数学模型、南海、承载力、钙质土、波浪破碎、波浪增水、波浪传播和波浪传播变形。从高频词统计词可以看出国内珊瑚礁模型研究与珊瑚礁成礁环境和水文动力学相关。

2.3 高频关键词的可视化分析

共词分析法是通过对某一学科领域的研究文献中的关键词共同出现在一篇文献中的现象进行分析, 从而判断该领域主题之间相互关系的研究方法^[29]。利用Bibexcel软件将国际上的43个关键词及词频进行两两配对, 形成43×43的共词矩阵, 随后将矩阵导入到Ucinet中生成##h文件, 再用Netdraw做出国际珊瑚礁生态系统模型研究的高频关键词共现网络可视图, 如图 2所示完成共词分析。网络图中, 节点大小表示关键词的中心度, 节点越大中心度越高, 说明该关键词与其他词共现的次数越多；节点间的连线表示两组关键词的关系, 线条越粗表示关系越密切, 两组共现次数越多^[30-31]。

如图 2所示, 与coral reef(珊瑚礁)共线次数较多的有climate change(气候变化)、coral bleaching(珊瑚白化)、biodiversity(生物多样性)等词, 说明这些研究方向是珊瑚礁生态系统模型研究的热门方向。气候变化对珊瑚礁生态系统的环境质量和内部生物生境产生直接影响, 造成珊瑚退化, 从而严重威胁珊瑚礁生态系统生物种群的生存和繁殖, 大幅度降低了生态系统内物种丰富度, 因此生态系统内生物多样性的追踪调查研究是当前研究领域的热点。其他关键词也由连线相互牵引形成了联系紧密、纵横交错的网状关系图。

国内研究因数量较少, 研究方向相对较分散, 并未出现共词现象。

2.4 高频关键词聚类分析

聚类分析法是一种以“物以类聚”为基本原理的多元统计方法。此方法能将大量数据按各自的特性来进行合理的分类, 将不同特性的数据分类到不同的类或者簇, 所以聚类分析后的同一个簇中的对象有很大的相似性, 而不同簇间的对象有很大的相异性^[32]。对国外文献中提取的43个高频关键词进行聚类分析, 按照关键词的不同属性共得到四个群集如表 2所示。

从表 2可以看出, 聚类分析后可以将高频关键词分为4个群集。群集1主要囊括了珊瑚礁生态系统中生活的各种生物因子和系统内各物种相互作用的种间关系, 如Acropora(鹿角珊瑚)、macroalgae(大型藻类)、Symbiodinium(共生藻)、Coral reef fish(珊瑚礁鱼)、competition(竞争)、Symbiosis(互利共生)等。珊瑚礁生物群落的种类组成丰富多彩, 是构成珊瑚礁生态系统的主体, 生物群落多样性很大程度上反映了生态系统的多样性。珊瑚礁上的生物分布具有明显的地方性, 生物间的生存竞争可以有效限制种群的扩散。群集2包括了珊瑚礁生态系统生物的生存环境因子以及环境变化引发的珊瑚疾病等研究, 如water quality(水质)、Habitat (生境)、resilience (水质)、Connectivity(连通性)和Coral bleaching(珊瑚白化)等。随珊瑚退化的日益加重, 为关注和保护珊瑚礁生态系统, 对珊瑚及生物的生活条件监测和模拟的模型研究日益受学者关注。群集3主要是针对特定区域如各地的海洋保护区的研究, 像Hawaii (夏威夷)、Red Sea (红海)、Great Barrier Reef(大堡礁)等。群集4是针对于气候变化对珊瑚礁影响的研究, 如global warming(全球变暖)、Sea surface temperature(海平面温度)等, 近些年来全球气候变化及一系列极端天气引发的珊瑚礁生态系统的破坏和损害引发学者高度重视, 海水的升温和降温, 海平面的升降变化等都会对生存于浅海的珊瑚礁产生巨大的负面影响。利用模型评估和预测气候变化对珊瑚礁生态系统的影响可以有效应对未来可能发生的损坏, 以便提出有效的防范和修复措施。

国内研究较少, 不便于聚类分析。但从承载力、钙质土表示国内对珊瑚礁成礁环境方面有所涉及；从波浪破碎、波浪增水、波浪传播和波浪传播变形说明国内对于珊瑚礁水文动力学研究较为关注。

3 结论与讨论

通过对珊瑚礁生态系统建模研究相关文献的统计与分析, 并对结果进行总结归纳, 得到如下结论：

按照时间顺序对发文量进行分析, 可将国外针对该学科的探究分为五个阶段：萌芽阶段, 缓慢发展阶段, 稳定发展阶段、波动增长阶段和加速增长阶段。各阶段涉及到的研究重点和热点有所不同, 但珊瑚礁成礁环境生态系统的生物多样性建模研究持续是热点。美国在该领域的研究始终处于遥遥领先的地位, 其次是澳大利亚、加拿大等拥有富饶珊瑚礁资源的国家。国内研究起步较晚, 发文量也相对匮乏, 构建的模型主要涉及珊瑚礁的礁体地形和波浪变形传播方面。国外在1974—2018共45年间发表了4057篇文献, 国内则在1997—2018共22年间发表了34篇文献, 可见国外发文量的增长幅度约是国内发文量的90倍, 因此国内应大力加强构建珊瑚礁生态系统模型的研究力度。

从高频关键词统计分析来看, 除Coral reef(珊瑚礁)外, climate change (气候变化)、Coral bleaching(珊瑚白化)、Coral reef fish(珊瑚礁鱼类)、Great Barrier Reef (大堡礁)和biodiversity(生物多样性)出现次数较多。国内文献中除珊瑚礁外, 波浪破碎、波浪增水、波浪传播和波浪传播变形等频率相对较高。透过这些高频关键词, 可以看出国内外用于研究珊瑚礁生态系统的模型各有侧重：国外主要致力于调查珊瑚礁退化原因及监测生态系统生物多样性的研究上, 而国内则将重点放在珊瑚礁成礁环境与波浪传播方面。

从高频关键词共现网络可视图可以看出：国际研究中climate change(气候变化)、biodiversity(生物多样性)、coral bleaching(珊瑚白化)等是珊瑚礁生态系统模型研究的核心与热点。而国内由于起步晚研究少, 高频关键词较为分散, 无共词现象。

从高频关键词聚类分析结果来看：国外的模型研究主要分为四个方面, 一是针对生态系统内生物种群及多样性的模拟研究；二是对于珊瑚礁生物生境质量和珊瑚疾病的调查研究；三是对于海洋保护区等特定区域的监测研究；四是针对气候变化给珊瑚礁生态系统带来影响的评估和预测研究。

综上所述, 当前全世界珊瑚退化、珊瑚礁生态环境恶化形势逼人, 珊瑚礁生态系统亟待修复。采用建模手段研究珊瑚礁生态系统可有效为保护和修复生态提供科学的理论和技术指导。珊瑚礁生态系统建模研究成果丰富, 当前国内外对于此研究的进展差别较大。国外通过计算机建模^[33]、动力学模拟^[34]、高光谱遥感和辐射传输建模^[35]、数学模型^[36-37]、基于模糊逻辑^[33]的珊瑚礁建模等技术手段, 对珊瑚礁生态系统各组成部分进行模拟、评估和修复等工作。国外的模型研究内容较全面, 如模拟珊瑚礁礁盘结构、监测珊瑚生物种群、修复珊瑚礁底栖环境等, 珊瑚礁生态系统的生物多样性研究持续是该领域的核心和热点, 同时用于预测珊瑚礁生态系统的风险、胁迫和修复的模型研究也在逐渐增多。而国内用建模方法研究珊瑚礁生态系统的研究进展较为落后, 常用的建模方法是构建数理模型, 主要用于对珊瑚礁礁体与波浪传播变形方面的数值模拟研究。相比之下, 我国在该领域的研究尚显薄弱, 还处于初步探索阶段, 所以我国应不仅局限于当前研究方向, 还应将重点尽快转到对珊瑚礁生态系统风险、威胁、胁迫和修复的评估和预测上, 相信今后我国珊瑚礁的生态修复工作会取得更加长远的成绩。