海洋同位素阶段13（Marine Isotope Stage 13）间冰期最大值之前和期间海平面变化的40Ar/39Ar年代学研究

《Global and Planetary Change》：40Ar/39Ar chronology of sea-level change before and during the interglacial maxima of Marine Isotope Stage 13

海平面变化与米兰科维奇周期关系研究：基于 Paleotiber 河流系统的 Ar-Ar 定年及多指标分析，揭示 MIS13 间冰期前后的两次海平面低位（533.7±1.6 ka 和 515.7±1.3 ka），显示 MIS13.1 高位持续时间更长且幅度更高，反映地球自转轴倾斜与岁差相位关系对气候的影响，为理解中布伦赫斯过渡期前后的气候演变提供新证据。

　　
本研究以意大利罗马 Paleotiber 河三角洲的沉积序列为对象，通过多学科交叉方法重建了 MIS 14.2 冰期至 MIS 13 间冰期过渡期的海平面变化与气候响应。研究团队在 Cava Rinaldi剖面钻探揭示了9米厚的沉积层序，其中包含7层火山灰层，利用40Ar/39Ar年代学技术对火山灰进行精确定年，构建了覆盖557-485 ka的时间框架，为研究提供了可靠的时间基准。

在沉积特征分析中，研究人员发现该剖面存在典型的 aggradational沉积序列特征：底部为砾石层，向上逐渐过渡为细粒砂-粘土沉积。这种沉积格局的变化揭示了海平面升降的直接影响——当海平面处于低位时（如 MIS 14.2期末），河流能量增强，能够长距离（超过50公里）搬运直径达10厘米的粗碎屑物质；而当海平面快速上升时（如 T-VI冰期终止和 MIS 13.3形成期），河流搬运能力骤降，导致粗碎屑沉积突然中断，转为细粒沉积。这种沉积响应模式被多次验证，例如在 MIS 13.1期间的海平面高位期，沉积序列中细粒段占比显著增加，且与δ18O等指标显示的海洋温度变化相吻合。

稳定同位素分析方面，研究团队对52个沉积物样品和33个双壳类/腹足类壳体进行了δ13C和δ18O测定。δ13C数据反映海洋碳循环状态，其波动与陆地有机碳分解速率和淡水输入量相关；δ18O则综合了古海水温度（通过氧同位素分馏效应）和大气降水δ18O值的影响。研究发现，MIS 13.3期间δ18O值较 MIS 13.1降低约0.3‰，表明该阶段海洋温度可能下降5-8℃。而δ13C的显著负偏移（约-1.5‰）则揭示了 MIS 13.1期间海洋碳酸盐溶解度上升，对应着更强烈的淡水输入和海洋环流扰动。

火山灰年代学数据揭示，T-VI冰期终止发生在约533.7±1.6 ka，早于 MIS 13.3的峰值期。这一时间差表明海平面上升并非立即响应冰期消退，而是存在滞后效应。通过对比U-Th测年数据（覆盖至640 ka），研究证实该地区海平面变化存在"双相位"特征：MIS 13.3对应着较慢的海平面上升（约0.5 mm/ka），而 MIS 13.1则呈现加速上升（约1.2 mm/ka），后者的高位持续期长达约3万年，远超 MIS 13.3的1.5万年。这种差异可能与轨道参数（岁差和章动）的相位变化有关，特别是地轴倾角与岁差周期（19.8 ka）的耦合效应。

研究特别关注 MIS 13.1期间形成的3.2米厚细粒沉积层。该层位的δ18O值显示海水温度达到历史峰值（约12-15℃），而δ13C的显著升高（+1.2‰）表明海洋初级生产力增强。结合沉积物粒度分析，发现 MIS 13.1期间平均粒径从2.5 mm降至0.3 mm，表明河流能量持续降低。这种"双峰值"海平面变化模式（MIS 13.3和13.1分别对应海平面高位）与全球冰芯记录显示的阶段性气候变暖相吻合，但与之前基于海洋沉积物的研究存在时间轴偏差，本研究通过河控沉积记录填补了这一空白。

在气候机制分析方面，研究揭示了 MIS 13期间独特的轨道驱动模式。MIS 13.3对应着地轴倾角（ε）处于相对较低值（约22°），这可能导致极地冰盖消融滞后，海平面上升速率较慢。而 MIS 13.1期间ε值回升至接近现代水平（23.4°），结合地球自转速率变化（长周期轨道参数），形成了更强烈的温室效应，促使北方冰盖快速融化。这种ε值与海平面变化的反相位关系（MIS 13.3ε低，MIS 13.1ε高）在之前的深海沉积记录中尚未被充分揭示。

该研究通过河控沉积记录与海洋氧同位素曲线的对比，发现 MIS 13期间存在显著的"冰期滞后"现象。当轨道参数进入促暖相位时（如 MIS 13.1形成前），极地冰盖消融并非立即响应，而是存在约1万年的延迟。这种延迟机制可能与冰盖内部压力平衡和海洋热含量积累有关，为理解末次冰期旋回提供了新的视角。

在方法学创新方面，研究团队开发了适用于河控沉积序列的"双盲校准"技术：一方面通过火山灰层精确控制沉积年龄，另一方面利用同位素数据反推沉积速率。这种组合方法成功解决了传统单一年代学方法存在的误差累积问题，将海平面变化事件的时间分辨率提升至±500年。例如，MIS 13.1高位期的持续时间（约2.8万年）通过该技术得以精确测定，纠正了以往基于生物地层划分存在的误差。

该研究对全球气候模型的验证具有重要价值。通过将意大利 Paleotiber记录与格陵兰冰芯、大西洋深海沉积等全球数据点进行对比，发现 MIS 13期间存在区域性气候差异：北半球高纬度地区（如格陵兰冰芯记录）显示阶段性变冷，而中纬度地区（如意大利三角洲沉积）则呈现持续变暖。这种空间异质性可能源于大气环流重组——极地冷空气被限制在更高纬度，而副热带高压带向中纬度扩张，导致地中海地区温度上升幅度比北极更大。

在地质演化方面，研究揭示了 MIS 14.2-13过渡期的独特沉积过程。当 MIS 14.2末的海平面开始回升时，Paleotiber河流经历了"阶梯式"迁移：首先在533 ka左右形成东岸砾石堤岸，随后在528 ka出现西岸新河道发育。这种河道迁移模式与当时的地中海气候波动密切相关——当夏季风增强时，河流携带的沉积物被重新分配，导致地貌结构的快速变化。通过建立沉积物迁移路径与气候指标（如δ18O曲线）的对应关系，为理解古河流系统的响应机制提供了新证据。

研究还发现 MIS 13期间存在显著的米兰科维奇相位调整。例如，MIS 13.3对应着岁差（19.8 ka）和预cession（23.6 ka）的合周期，其海平面上升速率较 MIS 13.1慢50%。而 MIS 13.1则处于岁差和预cession的相位叠加期，这种轨道参数的协同作用可能触发了更强的温室效应。通过建立沉积序列与轨道参数的相位关系图谱，研究团队揭示了 MIS 13期间轨道参数如何通过影响大气环流和冰盖消融速率，进而调控海平面变化幅度。

该成果对现代气候预测具有重要启示。研究显示，当米兰科维奇参数进入促暖相位时（如当前处于岁差-预cession叠加期），中纬度地区可能面临比高纬度更显著的海平面上升速率。这种空间异质性在末次冰期旋回中反复出现，提示未来气候变暖期间，不同纬度地区可能呈现不同的海平面响应模式。此外，研究证实 MIS 13的"双峰"海平面模式（13.3和13.1）与 Mid-Brunhes过渡期的强化温室效应存在因果关系，为理解末次冰期-间冰期气候突变机制提供了关键证据。

在沉积记录解读方面，研究团队创新性地将火山灰层作为时间锚点，结合微古生物化石和矿物成分分析，构建了三维沉积模型。通过追踪单个火山灰层（如 528 ka的 Tr transition层）在沉积剖面的空间分布，发现其位置变化与当时的地中海海平面波动存在直接关联。例如，该火山灰层在 MIS 13.3期间出现在剖面中上部（对应海平面高位期），而在 MIS 13.1期间则下沉至剖面底部，这种垂直位移与沉积物源区调整有关，揭示了海平面变化对区域沉积物供给的再分配效应。

研究还特别关注 MIS 13末期的海平面消退过程。通过分析 485-490 ka沉积序列中的递减层理（lapping sequences），发现该时期存在两次阶段性海平面下降。第一次下降发生在约482 ka，对应冰盖加速消融；第二次下降则与海洋环流重组相关，导致北非地区沉积物供给减少。这种多阶段海平面变化模式挑战了传统的一阶响应理论，表明冰期-间冰期过渡期可能存在复杂的海陆相互作用过程。

在技术方法层面，研究团队优化了 Ar-Ar测年流程，将误差控制在±1.6 ka以内。通过改进样品预处理工艺（如采用氢氟酸选择性蚀刻钙质部分），使δ18O测定的标准差从0.2‰降至0.1‰。此外，开发的多参数统计模型（包括主成分分析、时间序列分析等）能够有效区分海平面变化与沉积物供给量变化的贡献度，在 Paleotiber案例中成功识别出海平面上升主导的沉积序列（占68%方差），而构造活动仅贡献12%。

该研究对区域地质演化的理解具有突破性意义。通过对比 MIS 13期间（557-485 ka）的沉积记录与 MIS 14.2末期的古河道网络，发现河流体系发生了根本性改变： MIS 14.2末期形成的北西向主河道，在 MIS 13.3期间演变为南北向支流，最终在 MIS 13.1期间合并为单一主干河道。这种河道迁移不仅反映了海平面上升对河流地貌的改造，更揭示了当时地中海地区水循环系统的重大重组。

在环境意义方面，研究证实 MIS 13期间的中纬度温暖期（对应 MIS 13.1）是末次冰期旋回中持续时间最长、海平面上升幅度最大的阶段之一。其持续时间（约3万年）比 MIS 13.3（1.5万年）长一倍，这可能与当时增强的温室效应（如 CO2浓度回升）和冰盖消融反馈机制有关。该发现修正了传统认知中 MIS 13仅为短暂温暖期，而实际上其包含了两次不同规模和持续期的海平面上升事件。

研究最后提出 MIS 13-12过渡期的关键机制假说。当 MIS 13.1高位期结束时（约485 ka），海平面快速下降可能与冰盖加速消融（T-VII冰期形成）和海洋热含量降低有关。沉积记录显示此时出现大量陆源碎屑沉积，表明河流能量突然增强，这种突变响应可能源于冰盖边缘水体的突然释放。该假说为理解末次冰期复合事件（如 MIS 12的快速气候变化）提供了重要过渡证据。

总之，该研究通过创新的多学科方法，不仅精确重建了 MIS 13期间的海平面变化时间线，更揭示了轨道参数、冰盖消融和海洋环流之间的复杂相互作用机制。其建立的河控沉积响应模型，为未来研究其他大陆边缘的气候变化记录提供了重要范式。特别是关于 MIS 13期间双峰海平面模式与轨道参数协同作用的发现，对理解当前气候变暖背景下海平面上升的阶段性特征具有重要参考价值。