南极沿海地区不同断面的表层雪化学和稳定同位素组成的空间变化

《Water》：Spatial Variations in Surface Snow Chemistry and Stable Isotope Composition Along Different Transects of the Antarctic Coastal Zone

　　

本研究阐明了南极沿海地区表层雪化学和稳定同位素组成的空间分布特征及其影响因素。研究人员分析了八个表层雪化学离子断面和五个表层雪稳定同位素断面的数据，以研究海盐离子（Na⁺、Cl^?、Mg²⁺、K⁺）、非海盐离子（NO₃^?、甲磺酸（MSA）、非海盐硫酸盐（nssSO₄^2?）、非海盐钙（nssCa²⁺））以及稳定同位素δ¹⁸O的分布。化学离子浓度表现出明显的源特征。主成分分析表明，南极沿海地区离子源的排序为：海盐气溶胶输送，其次是海洋生物活动和地壳粉尘，然后是大气化学反应。所有断面的数据显示，δ¹⁸O随距海岸距离增加而降低的总体模式，但δ¹⁸O–距离斜率（?6.84至?1.98‰/100 km）、δ¹⁸O–海拔斜率（?0.97至?0.20‰/100 m）和δ¹⁸O–温度斜率（0.61至1.14‰/°C）存在显著区域差异。水汽源区与输送路径、局地地形抬升以及局地温度场共同影响沿海地区δ¹⁸O的分异。研究结果表明，南极沿海地区的雪化学和同位素组成受海洋源输送、生物地球化学循环以及极地地形和气候过程的控制。

**研究背景与问题**
南极沿海地区位于冰、海、气相互作用的前沿，受南极绕极流、局地海洋气旋及冰架融化等因素影响，成为全球气候变化响应最为敏感的区域之一。极地雪冰中的化学离子和稳定同位素是重建现代气候过程与古环境变化的重要指标。海盐离子（如Na⁺、Cl^?、Mg²⁺、K⁺）主要源自海水或海冰，受海冰范围和气象条件强烈影响；非海盐离子（如甲磺酸（MSA）、非海盐硫酸盐（nssSO₄^2?）、非海盐钙（nssCa²⁺））则主要来自海洋生物活动和地壳粉尘；NO₃^?主要源自平流层。稳定同位素δ¹⁸O和δD可用于指示气候的时间和空间变化，而氘盈余（d-excess）则有助于识别水汽源区。然而，现有研究大多集中于单一断面或局地区域，缺乏对不同沿海区域化学离子和稳定同位素变化的系统比较分析。特别是局地过程（如吹雪、海冰融化）与大尺度环流的协同效应尚未得到充分阐明。为此，研究人员利用南极沿海八个表层雪化学离子断面和五个稳定同位素断面的数据，系统探讨了海洋气溶胶输入、大气环流输送、局地地形及生物活动等因素对雪化学特征的综合影响。该研究为增强南极冰盖沿海地区气候与环境变化的观测网络提供了关键数据，并为深入理解极地海-陆相互作用及全球气候系统耦合过程奠定了科学基础。论文发表在《Water》。

**关键技术方法**
研究人员主要采用以下三种关键技术方法：（1）**富集因子（EF）**：以Na⁺为海盐参照，计算Cl^?、K⁺、Mg²⁺、Ca²⁺和SO₄^2?相对于海水的富集程度，以区分海盐与非海盐来源；（2）**主成分分析（PCA）**：对六个化学离子断面（排除DD1和PV断面）进行降维处理，提取控制化学组成变化的主要因子；（3）**δ¹⁸O–距离、δ¹⁸O–海拔、δ¹⁸O–温度斜率**：通过线性回归计算五个断面的斜率，评估不同沿海区域稳定同位素分馏的差异。数据来源包括：CHINARE第18次南极考察（2001–2002）期间自中山站至LGB69断面采集的17个表层雪样品（化学离子和稳定同位素分析），以及ITASE数据库（http://www2.umaine.edu/itase）和Pangaea数据库（https://doi.pangaea.de）中的七个化学离子断面和四个同位素断面数据。

**研究结果**
**3.1 化学离子的空间分布与来源分析**
**3.1.1 化学离子的空间分布**
海盐离子（Na⁺、Cl^?、Mg²⁺、K⁺）浓度随海拔变化呈现显著区域差异。在东南极PEL区域ZL断面的低海拔段（0–600 m），海盐离子浓度达到最高（Na⁺: 970.78–1458.92 ppb；Cl^?: 952.48–1659.08 ppb）；在东南极VL区域DD1断面，海盐离子浓度总体最低（Na⁺: 4.6–437 ppb）。东南极各断面（JJ、ZL、TT）的海盐离子浓度均随海拔呈指数衰减，符合海盐离子的典型变化特征。西南极沿海区域海盐离子浓度随海拔显著降低，PV断面浓度远高于GD断面；西南极沿海断面海盐离子浓度总体高于东南极。南极半岛（AP）DD2断面海盐离子浓度持续较高，无明显趋势。非海盐离子（NO₃^?、MSA、SO₄^2?、Ca²⁺）的空间分布也差异显著：NO₃^?浓度范围为8–200 ppb，最小值在ZL断面（8.55 ppb），最大值在FF断面（207.4 ppb），且不随海拔呈系统变化；SO₄^2?浓度范围为10–450 ppb，JJ断面浓度较低，但高海拔区nssSO₄^2?贡献较高；MSA浓度范围为2–150 ppb，ZL断面值最高（2.77–145.82 ppb），可能与强生物活动有关；Ca²⁺和nssCa²⁺浓度普遍较低且稳定，表明长距离大陆粉尘输入。

**3.1.2 化学离子的来源分析**
七个断面（除DD1）的富集因子（EF）分析显示，Cl^?和Mg²⁺的EF值接近1，符合海水来源，主要来自海盐气溶胶；SO₄^2?和Ca²⁺的EF值显著大于1，表明非海水富集；K⁺在JJ、TT和FF断面EF值远大于1，暗示额外来源。主成分分析（PCA）对六个断面（排除DD1和PV）的结果显示，PC1方差贡献率最大（53.00%–73.54%），主要受海盐离子控制；PC2贡献率次之（14.42%–31.10%），主要由海洋生物活动（MSA、nssSO₄^2?）和地壳粉尘（nssCa²⁺）指示物主导；PC3贡献率最小（11.40%–13.62%），与大气化学反应产物NO₃^?相关。不同断面PCA模式存在差异：DD2断面PC1中生物/陆地指示物（nssCa²⁺、MSA、SO₄^2?、Ca²⁺）的载荷超过海盐离子，可能与南极半岛靠近巴塔哥尼亚粉尘源区和威德尔海生物生产力高有关。

**3.2 δ¹⁸O的空间变化**
五个沿海断面的δ¹⁸O均随距海岸距离增加而降低，但变化速率差异显著（?6.84至?1.98‰/100 km）。东南极北缘的ZL'和VV断面δ¹⁸O下降最陡（?6.84至?5.06‰/100 km），其中ZL'断面在距海岸7.4–20.6 km范围内δ¹⁸O升高，可能与陡坡强抬升有关。VL区域的II断面下降较缓（?1.98‰/100 km，R²=0.28），反映区域过程的复杂性。δ¹⁸O–海拔斜率范围为?0.97至?0.20‰/100 m，ZL'和VV断面斜率最陡（?0.75和?0.97‰/100 m），SS断面最缓（?0.20‰/100 m，R²=0.19）。δ¹⁸O–温度斜率范围为0.61至1.14‰/°C，其中ZL'断面为0.80‰/°C，VV断面为0.89‰/°C，II断面为0.96‰/°C，GD'断面为1.14‰/°C，SS断面为0.61‰/°C，总体与经典Dansgaard关系一致，表明温度是主要控制因素。

**讨论**
海盐离子的空间分布反映了海洋气溶胶随距离的衰减，东南极海盐离子浓度高于西南极，与海冰动态和海洋气溶胶输送的区域差异相关。非海盐离子NO₃^?无沿海趋势，支持平流层来源假设；ZL断面MSA浓度极高（达146 ppb），与2001/2002年夏季异常大气环流导致的海冰融化及浮游植物大量繁殖有关。PCA结果证实海盐气溶胶沉积是沿海区域的主导过程，但断面间差异受生物热点距离、粉尘输送路径和局地微气候影响。δ¹⁸O的空间变化受水汽源区与输送路径、地形抬升和局地温度场共同控制：ZL'、VV和GD'断面δ¹⁸O–距离斜率最大，受海洋气团影响强烈；II断面变化缓慢，与较稳定的局地环流有关；SS断面温度敏感性低（0.61‰/°C），可能反映局地二次蒸发或气团混合。

**结论**
本研究分析了南极沿海八个表层雪化学断面和五个稳定同位素断面的空间分布与来源特征，主要结论如下：
- 东南极海盐离子（Na⁺、Cl^?、Mg²⁺、K⁺）浓度总体高于西南极，与海冰动态和海洋气溶胶输送的区域差异相关。非海盐离子（NO₃^?、MSA、SO₄^2?、Ca²⁺）呈现不同的源特征。
- 富集因子（EF）表明Cl^?和Mg²⁺主要来自海洋气溶胶，nssSO₄^2?和nssCa²⁺来源相反，K⁺兼具海盐和非海盐来源。主成分分析（PCA）进一步表明，南极沿海离子源排序为：海盐气溶胶输送，其次为海洋生物活动和长距离输送矿物粉尘，然后为大气化学反应。
- 所有断面均显示δ¹⁸O随距海岸距离增加而降低的总体模式，但变化速率差异显著（?6.84至?1.98‰/100 km）。II、SS和GD'断面δ¹⁸O随距离波动剧烈，表明局地微气候和地形对稳定同位素分布有强烈影响。南极沿海断面（ZL'、VV、GD'）表现出最强的同位素海拔效应，而南极半岛SS断面响应较弱。
- II和GD'断面较高的δ¹⁸O–温度斜率（0.96至1.14‰/°C）表明这些区域δ¹⁸O对温度变化异常敏感，可能与强海洋性气候下稳定降水保存了原始温度信号有关。SS断面温度敏感性较低（0.61‰/°C），表明局地气团混合或蒸发效应增强。
本研究存在以下局限性：表层雪数据集跨越较长时期，难以消除短期气象事件（如暴风雪）的影响，并可能引入年际气候变率（如ENSO）；部分断面（如DD1和FF）的采样日期不明确，可能产生季节性偏差；δ¹⁸O变率归因于局地过程的定性分析，未进行轨迹或模型分析。
本研究刻画了南极沿海环境中物质的空间分异格局及其来源，为研究极地生物地球化学循环提供了新数据和理论框架。所提出的“局地过程主导区域差异”范式可用于预测全球气候变化影响下的极地放大效应。