红树林湿地是全球四大高生产力海洋生态系统之一，平均生产力(以碳计)为2 500 mg/(m²·d)，其凋落物约占进入海洋陆源有机质的11%^{[1, 2]}。因此说，红树林湿地中的红树林自身在维护“陆—海—气”系统中的碳循环及近岸生态系统平衡和生物多样性方面发挥着重要作用^{[3, 4]}。但是，在全球气候变暖情景下，海平面上升势必对长于潮间带的红树林生存空间造成威胁，为与海平面变化保持同步，会迫使红树林发生向陆迁移和群落兴衰等问题^[5]。红树林迁移效应的驱动力源于相对海平面上升速率和红树林潮滩地表高程抬升速率的大小关系^[6]。当前，国外潮滩地表高程变化测量的一种可信、准确的方法是利用地表高程测量仪^[7]，但在国内尚未普遍推广；鉴于我国目前广泛使用²¹⁰Pb沉积速率计算的背景下，进一步深化沉积物的压实作用对沉积速率和地表高程抬升速率的影响不失为一种行而有效的方法。以往鲜有研究把红树林潮滩地表高程抬升速率与沉积速率区别开来，更有甚者可能将其混为一谈。因此，红树林潮滩地表高程抬升速率与沉积速率之间的关系研究已成为解决问题的关键。

²¹⁰Pb已被国内外广泛用于湖泊、河流和海洋现代沉积速率的研究中^{[6, 8, 9, 10, 11, 12, 13, 14]}。由于上覆沉积物的压实作用致使各层深的孔隙度和密度存在差异，不同程度的改变了²¹⁰Pb的垂向分布特征^[15]；尤其是对淤泥/淤泥质粉砂等压缩率较大的沉积物，更要注意沉积物的压实作用^[16]。累加质量深度是消除孔隙度和密度差异的常用压实校正方法^{[6, 8, 17, 18]}，采用相同的压实校正标准，便于横向比较不同海区的沉积速率。有鉴于此，本文基于广西典型港湾红树林区8根短柱状沉积物，在考虑压实作用的基础上计算红树林潮滩沉积速率，通过与地表高程抬升速率进行对比研究，进而计算当前潮滩地表高程抬升速率；并结合有关北部湾相对海平面变化资料，揭示当前海平面上升对广西红树林向陆/向海迁移的驱动机制。

1.1 样品采集

2011年5月，趁低潮时将直径10 cm的有机玻璃管缓缓打入广西典型红树林区潮滩沉积物中，共采集了8根短柱状样(图 1)。其中，SJC、HXL和JXW分别位于茅尾海的老鼠簕、秋茄和桐花树群落内(图 1a)；DDH位于丹兜海的白骨壤群落内；YLW01、YLW02、YLW03和YLW04分别位于英罗湾的红海榄、秋茄、木榄和桐花树群落内(图 1b、表 1)。在样品获取过程中尽量保持沉积物的原始状态，并利用美国Trimble 公司的DSM 212H 双信标DGPS 定位设备进行现场定位，平面定位精度优于1 m。将柱状样进行现场密封后带回实验室，按2 cm间隔进行分样，将子样品密封冷冻直至实验室分析。样品在采集、运输、保存和测定的过程中，均严格按照《海洋调查规范》等相关要求和规定进行。

图 1 研究区域和采样站位分布图 Fig. 1 Maps showing the study area and coring location

图选项

取约20 g湿样品于称量盒内，称重后一并放入烘箱，在105℃~110℃恒温下烘10 h以上。将烘干后的称量盒和样品取出，放入干燥器内冷却后再称重，将两次平行试验的误差控制在0.5%内，以此计算沉积物的含水率。

对入管蜡封并放置3星期的约5 g干样品进行²¹⁰Pb和²²⁶Ra放射性比活度测试，采用EG&G Ortec公司生产的高纯锗低本底γ能谱仪完成。以46.5 keV(²¹⁰Pb)处的能量峰来计算总²¹⁰Pb比活度(²¹⁰Pb_tot)，以351.92 keV(²¹⁴Pb，²²⁶Ra的子体)处的能量峰来计算本底²¹⁰Pb比活度(²¹⁰Pb_sup)，其差值即为过剩²¹⁰Pb(²¹⁰Pb_ex)的比活度^[19]。

²¹⁰Pb定年: 选用常量初始浓度(CIC)模式^[14]。

沉积物干密度：前人研究表明，高度压实的页岩干密度介于2.55~2.65 g/cm³；因此，沉积物干密度可通过成熟的经验公式进行估算^[10]：

BD = (PD·D) / (D+PD·(1-D))

式中：BD 是沉积物干密度(g/cm³)； PD 是页岩干密度(孔隙度近似为0)，本文取2.60 g/cm³； D 是沉积物含水率(%)。

压实作用: 是指沉积物沉积后，由于上覆沉积物不断加厚，在重荷压力下所发生的脱水、孔隙度降低、体积缩小、密度增大的地质过程。由于含水率变化引起沉积物密度差异，压实校正被广泛应用于²¹⁰Pb剖面^[6]。本文，各层位子样品利用柱样的平均干密度进行校正，即公式(1)：

式中：CI_x 是层位x的压缩间隔长度 (cm)；BD_x 是层位x的干密度 (g/cm³)； BD_avg 是柱样的平均干密度 (g/cm³)；I 是原始层位间隔，2 cm。

把各层位进行压实校正后的CI_x 累加在一起，得到层位x的校正深度： CI₁+CI₂+CI₃+…+CI_x。

广西典型红树林区8根短柱状样中²¹⁰Pb_tot的比活度并非呈现出理想状态下随深度呈指数衰减的趋势，局部层位偶现异常的增减波动，但总体向下减少的趋势还是比较明显(图 2)。而SJC、DDH和YLW03柱中表层10 cm内的²¹⁰Pb_tot 却呈现相对稳定/衰减的异常现象，这应与沉积物表层混合作用有关^[20]。鉴于广西海岸带非封闭体系，受浪、潮、流以及人类活动的影响显著，因此选用CIC模式定年更合理一些。它适用于沉积物主要来源于表层侵蚀产物，²¹⁰Pb 含量明显受物源影响，即沉积物增加同时导致相应²¹⁰Pb 增加的沉积系统^[21]。

图 2 210Pbtot比活度的垂向分布特征 Fig. 2 Vertical distributions of 210Pbtot activities in the sediment cores

图选项

未考虑压实作用下，选用深层沉积物中的²¹⁰Pb_tot作为²¹⁰Pb_sup^[19]，其²¹⁰Pb_ex比活度的垂向分布特征如图 3所示。将²¹⁰Pb_ex与深度进行指数拟合，拟合系数R²介于0.58~0.89之间，拟合程度良好。计算可得，SJC、HXL、JXW、DDH、YLW01、YLW02、YLW03和YLW04柱的平均沉积速率分别为1.04、0.52、0.52、0.19、0.33、0.31、0.26和0.22 cm/a。其中，位于茅尾海钦江口附近SJC、HXL、JXW柱的沉积速率明显高于周边无大型河流入海的丹兜海DDH和英罗湾YLW01、YLW02、YLW03和YLW04柱，可见河流输沙入海是影响潮滩沉积速率的主要原因。

图 3 210Pbex随深度的垂向分布及其指数拟合曲线 Fig. 3 Vertical distributions of 210Pbex activities and its exponential curves fitting with depth

图选项

沉积物中含水率自上而下递减趋势明显(图 4)，表层沉积物中的含水率高达40%~50%，而深层沉积物中的含水率低至20%~30%，沉积物压实作用显著；且深层沉积物中的含水率趋于相对稳定。这源于上覆沉积物不断加厚，在重荷压力下所发生的脱水、孔隙度降低等一系列的地质过程。基于线性回归模型对含水率随深度的变化进行拟合，拟合系数R²介于0.38~0.83，表明含水率随深度多呈线性衰减趋势。

图 4 柱状沉积物中含水率的垂向分布及其线性回归 Fig. 4 Vertical distributions of moisture contents and its simple linear regressions with depth

图选项

由于含水率变化引起的沉积物密度差异，压实作用被用于校正²¹⁰Pb剖面^[6]。本文，将各层位沉积物利用柱样的平均干密度进行校正，再把各层位进行压实校正后的CI_x 累加在一起，得到层位x的校正深度。将²¹⁰Pb_ex 的比活度与校正深度进行指数拟合(图 5)，拟合系数R²介于0.63~0.94之间，拟合程度较未考虑压实作用有大幅提高。计算得，SJC、HXL、JXW、DDH、YLW01、YLW02、YLW03和YLW04柱的压实沉积速率分别为0.78、0.44、0.52、0.16、0.31、0.26、0.24和0.22 cm/a。

图 5 210Pbex随校正深度的垂向分布及其指数拟合曲线 Fig. 5 Vertical distributions of 210Pbex activities and its exponential curves fitting with depth correction

图选项

为了比较压实校正前/后沉积速率的差异，将8根短柱中的沉积速率绘于图 6(a)中。显而易见，未压实校正沉积速率均大于压实校正沉积速率，约是压实校正沉积速率的1.00~1.34倍(平均1.12倍)。将由压实校正前/后沉积速率换算出的沉积物底部年龄绘于图 6(b)中，由于压实校正年龄应更接近于实际沉积年龄，故SJC、HXL、JXW、DDH、YLW01、YLW02和YLW03柱压实校正的沉积物年龄相比未压实校正的年龄分别增加了31、27、0、39、10、52和14年，而YLW04柱的年龄却减少了4年。

图 6 压实校正前/后沉积速率和沉积物年龄的关系 Fig. 6 Calculated sedimentation rates and sediment ages before/after compaction correction

图选项

考虑压实作用下，由压实沉积速率和CI_x压缩间隔反推各年份不同含水率下的理论沉积速率(非真实沉积速率)，以及无压实作用下的表层沉积速率(图 7)。由于源自含水率的推算，除表层沉积速率为当年实际沉积速率外，其余各年份的沉积速率仅反映当年实际沉积速率被不同程度压实后的现状，属于理论沉积速率，它与含水率的时间演变完全一致，用于研究压实作用下的浅层沉积物下沉速率。

图 7 压实作用下各年份的理论沉积速率 Fig. 7 Temporal changes of theoretical sedimentation rates under compaction

图选项

为了查清潮滩地表高程抬升速率与沉积速率的关系，将未考虑/考虑压实作用两种情况下的红树林潮滩沉积速率和地表高程抬升速率的示意图绘于图 8中(每层代表一个年度内的沉积物累积厚度)。

图 8 压实校正前/后红树林潮滩沉积速率和地表高程抬升速率示意图 Fig. 8 Diagrammatic sketches of sedimentation rates and rising rates of seawater-sediment interface before/after compaction correction

图选项

(1) 未考虑压实作用下：地表高程抬升速率=未压实沉积速率(图 8a阴影部分)；

(2) 考虑压实作用下：由于新一年度的沉积物对其下方的沉积物起到再压实作用，地表高程抬升速率 (图 8c阴影部分)要小于红树林潮滩表层沉积速率(图 8b阴影部分)；地表高程抬升速率等于“表层沉积速率—表层沉降速率”。通过对比相邻两个年度的沉积物压实模型(图 8右)发现，地表高程抬升速 率(图 8c阴影部分)等于底层压实沉积速率(图 8d阴影部分)。

3.1 相对海平面上升速率

全球海平面上升速率：据IPCC(2007)^[22]，1961~2003年间全球平均海平面以每年~1.8 mm(1.3~2.3)的平均速率上升；其中，1993~2003年间，全球平均海平面以每年约3.1 mm(2.4~3.8)的平均速率上升。在未考虑气候—碳循环反馈的不确定性(也未包括冰盖流量变化的整体效应)前提下，预测2090~2099年较1980~1999年的海平面上升18~59 cm。

相对海平面上升速率：据中国海洋公报(2010)^[23]，2001年以来广西沿海的海平面总体处于历史高位，2001~2010年的平均海平面比1991~2000年的平均海平面高约22 mm，比1981~1990年的平均海平面高约48 mm。即广西沿海相对海平面平均上升速率为2.2~2.4 mm/a，可作为广西沿海相对海平面上升速率；它略低于全球平均海平面上升速率，源于区域板块构造的抬升^[24]。

在临界速率下，红树林随相对海平面上升发生同步迁移；反之，红树林迁移跟不上相对海平面上升步伐，发生衰退。

(1) 红树林后缘地貌和地层条件适合红树林生长和迁移的区域

不同学者研究获得的临界速率大小不一：① Ellison & Stoddart^[25]认为临界速率为0.8~0.9 mm/a，超过1.2 mm/a将会灭亡；② Snedaker et al.^[26]认为临界速率为1.6~1.8 mm/a；③ Ellison et al.^[27]认为临界速率为2~10 mm/a。本文选用1.2 mm/a作为红树林适宜迁移的临界速率。

(2) 红树林后缘为海堤的情况

后缘海堤限制了红树林的向陆迁移，临界速率为零。

据广西908专项实测数据：广西海岸线总长1 628.59 km，其中人工岸线长1 280.21 km (占78.6%)。受近年围塘养殖、盐场建设、港口围填以及人工岸堤修建等的影响，人工岸线长度仍呈逐年增加的趋势；即相对海平面上升速率大于地表高程抬升速率将会对红树林造成毁灭性的影响。基于未考虑/考虑压实作用两种模式，计算广西典型红树林区的相对海平面上升速率(V_{海平面上升})与地表高程抬升速率(V_地表高程)的差值，详见图 9。

图 9 压实校正前/后红树林区V海平面上升与V地表高程的差值 Fig. 9 The difference between relative sea-level rise rate and ground elevation uplift rate before/after compaction correction

图选项

(1) 未考虑压实作用：除丹兜海的DDH和英罗湾的YLW04柱外，其余6柱的相对海平面上升速率均小于地表高程抬升速率，表明滩面多处于淤积状态。这与张乔民研究获得的结论相一致^{[27, 28]}：“我国大部分红树林潮滩淤积速率接近或大于2030年前的海平面上升速率，红树林潮滩可以通过滩面淤积跟上甚至超越海平面上升，红树林面积基本上能保持稳定”^{[28, 29]}。

(2) 考虑压实作用：丹兜海的DDH和英罗湾的YLW01、YLW02、YLW03、YLW04五根短柱的相对海平面上升速率均大于潮滩地表高程抬升速率，表明该区红树林已经受到海平面上升的轻度威胁。

茅尾海的SJC、HXL和JXW，无论是否考虑压实作用，相对海平面上升速率均小于地表高程抬升速率，表明该区红树林暂未受到海平面上升的影响。可见，压实作用校正与否对地表高程抬升速率与相对海平面上升速率相当的区域尤为重要。

用²¹⁰Pb法测定红树林潮滩沉积速率/地表高程抬升速率时均应考虑压实作用的影响，尤其在含水率变化显著的潮滩。本文基于考虑/未考虑压实作用两种模式，分别计算了沉积速率与地表高程抬升速率的关系。研究发现：①无论是否考虑压实作用，计算出的茅尾海红树林区的地表高程抬升速率均大于相对海平面上升速率，这与近期茅尾海的淤积现象相符；②考虑压实作用下，英罗湾和丹兜海红树林区的地表高程抬升速率小于相对海平面上升速率；发现地表高程抬升速率等于压实沉积速率，明显低于未经压实的表层沉积速率。由于广西红树林海岸大都建有防波堤，限制了红树林向陆方向的迁移；因此，英罗湾和丹兜海的红树林正面临海平面上升的威胁。这与前人部分研究结果有所出入，主要源于未考虑沉积物压实作用，过高的估计了红树林区的沉积速率/地表高程抬升速率。可见，压实作用校正与否对地表高程抬升速率与相对海平面上升速率相当的区域尤为重要。