1985–南海诸岛珊瑚礁区热压力时空变化研究分析.docx

? ? 1985–南海诸岛珊瑚礁区热压力时空变化研究分析 ? ? 陆永强，陈正华*，余克服*，何鑫，张威，兰思香 ( 1. 广西大学 海洋学院，广西 南宁 530004；2. 广西大学 广西南海珊瑚礁研究重点实验室，广西 南宁 530004；3. 广西大学珊瑚礁研究中心，广西 南宁 530004) 1 引言 珊瑚礁是世界上生物多样性最丰富的生态系统，被称为海洋里的“热带雨林”，对人类和海洋生物有着重要的意义[1–2]。近几十年来全球气候变暖造成海表温度（SST）上升[3]，导致大规模珊瑚礁白化和死亡[4]。SST超过夏季最高温度1℃，珊瑚便会受到热压力影响排出体内共生的藻类（或共生虫黄藻失去体内色素），珊瑚露出白色的碳酸钙骨架，这种现象称为白化[5–6]。若热压力增强或持续时间过长，珊瑚就会死亡[7]。珊瑚白化对珊瑚礁的生物多样性和生态系统服务功能产生了严重的影响，预计21 世纪末全球暖化将导致珊瑚礁减少75%[8]。 我国珊瑚礁主要分布在南海诸岛，面积约为30 000 km2，占世界珊瑚礁面积的5%[9–10]。研究发现南海诸岛珊瑚礁正处于快速退化状态[11–12]，如西沙群岛2005–2006 年活珊瑚平均覆盖度为65%～70%[13]，到2015 年西沙群岛活珊瑚平均覆盖度为16.3%[11]。南海诸岛珊瑚礁是易受热压力影响的地区之一[14]。有研究表明南海诸岛珊瑚礁的热压力、白化和死亡与厄尔尼诺（El Ni?o）存在紧密的联系[15]，El Ni?o 一方面通过大气桥理论与南海进行大气环流[16]，另一方面吕宋海峡海流将El Ni?o 与南海相连[17]。珊瑚白化与热压力有直接关系[18]，在南海诸岛珊瑚礁退化的背景下，亟需对南海诸岛珊瑚礁区热压力及El Ni?o 对其影响展开研究。对于珊瑚礁热压力和白化的研究，传统研究手段大多集中在部分岛礁，南海诸岛珊瑚礁分布零散，研究人员无法定期访问所有珊瑚礁区域[19]。遥感具有覆盖范围广、时间空间分辨率高、成本低、可实现连续观测等特点，随着遥感技术的迭代升级，大范围的珊瑚礁区监测可以通过遥感卫星实现[20–21]。 本文将使用美国国家海洋和大气管理局（NOAA）珊瑚礁监测中心（CRW）的1985–2019 年每日5 km分辨率SST 数据集[22–23]，对南海诸岛珊瑚礁区的热压力进行时空变化分析，基于热压力评估珊瑚礁白化和死亡的风险。一方面可以探究珊瑚礁对大规模热压力的抵抗能力[24]，另一方面可以找出珊瑚礁热胁迫最小的地点，以便找出合适珊瑚礁区避难所[25–26]。珊瑚白化与El Ni?o 显著相关[27]，本文将对珊瑚礁区热压力和El Ni?o 进行相关性分析，量化El Ni?o 对南海诸岛珊瑚礁区热压力变化的影响。对南海诸岛珊瑚礁热压力的研究，有利于科学管理南海珊瑚礁，推进珊瑚礁生态系统的可持续发展[28]。 2 数据来源和方法 2.1 研究区域 珊瑚礁数据从全球珊瑚礁分布（Global Distribution of Coral Reefs）网站获取，网址/datasets/1。南海诸岛珊瑚礁远离陆地受到人类活动影响较小，本文将南海诸岛珊瑚礁作为研究区。为了准确评估珊瑚礁区的热压力时空变化，本文将南海珊瑚礁进行20 km 缓冲扩展来表示南海珊瑚礁研究区[29]。缓冲后的珊瑚礁区主要分布在东沙群岛、西沙群岛、中沙群岛（黄岩岛）和南沙群岛（图1a）共4 851个像元，因此本文便将南海珊瑚礁研究区分为东沙珊瑚礁区（20°20′～21°10′N，116°26′～117°10′E）、西沙珊瑚礁区（15°32′～17°18′N，110°59′～112°58′E）、中沙珊瑚礁区（14°53′～15°25′N，117°29′～118°03′E）和南沙珊瑚礁区（7°03′～12°13′N，109°23′～118°02′E）。 图1 南海诸岛珊瑚礁研究区Fig. 1 Location of coral reef regions in the South China Sea islands 2.2 数据来源 2.2.1 海表温度 海表温度数据选取自美国国家大气与海洋中心管理局珊瑚礁观测计划中心( /index.php)的两组数据集：（1）1985 年1 月1 日至2019 年12 月31 日的每日更新5 km SST 数据集，该数据集是CRW 在2014 年发布的每日5 km 分辨率的海表温度产品，时间跨度是1985 年至今[23]；（2）3.1 version 5 km 分辨率的最热月平均数据（The Maximum Monthly Mean, MMM），MMM 代表1985–2012 年每个像元最热月平均气候SST 值[30]。 2.2.2 ONI 本文选取热带太平洋（5°N～