晚中新世柴达木盆地低偏心率时期倾角驱动的干湿变化

邓欣宜，聂军胜，任雪萍

兰州大学资源环境学院，西部环境教育部重点实验室， 兰州 730000

人类目前生活的地质时期是1.17万年至今的全新世，这是一个地球轨道偏心率较低、太阳辐射振幅变化较小的间冰期[1-2]。在工业化革命后温室气体持续增排的背景下，到21世纪末大气CO2浓度可能上升到（500～600）×10-6[3]，远高于北半球冰盖形成时的阈值。未来全球气候可能会继续变暖，演变成北半球夏季无冰的“单极冰室模式”[4-5]。晚中新世是一个距今较近的温暖期，该时期仅在南极形成了大规模冰盖，北半球尚未形成大规模的稳定冰盖[6-7]。因此，分析过去暖期的干湿变化，如晚中新世的低偏心率期气候干湿变化规律及其对地球轨道参数的响应机制，对预测未来亚洲地区气候干湿变化具有重要参考价值[8-16]。

近几十年来，国内外许多学者基于中国黄土[17-25]、石笋[26-27]、湖泊[24,28-29]、海洋沉积物[30-33]等载体，对与当前天文背景（低偏心率下的全新世）相似的温暖间冰期，如深海氧同位素11和19阶段（MIS 11和MIS 19）等，开展了广泛的研究，为剖析区域和全球气候演化的驱动机制奠定了重要基础。以MIS11为例[34]，东亚高分辨率石笋氧同位素记录的轨道尺度降水变化具有明显的2万年周期，支持了北半球夏季太阳辐射对亚洲夏季风降水的控制[26-27]，而中国黄土磁化率记录的降水变化存在显著10万年周期，这与北半球冰量的周期性变化一致[14,35-37]。此外，模拟结果表明，在低偏心率时期，倾角变化是控制气候变化的主要因素。例如，Wu 等[38]利用大气、生物、海洋和海冰地球系统模型（LOVECLIM）模拟了南北两个半球海温和海冰在417～511 ka期间对岁差和倾角的不同响应，发现低偏心率时期（417～460 ka ）倾角变化发挥主导作用。同时，第四纪约1～2 Ma时期大气环流模型（GCM）的结果也表明，低偏心率时期倾角变化控制着南北半球高纬度冰量和气候变化[39]。上述研究主要集中在北半球有冰时期，北半球没有永久冰盖时低偏心率时期气候变化研究还很薄弱，气候变化如何响应轨道参数变化的机制尚不清楚，这些都限制了对亚洲内陆环境未来变化的预测[40]。

以往对中国北方干旱-半干旱地区高分辨率气候变化的研究偏重从较长尺度理解轨道尺度气候干湿变化规律，少有针对低偏心率这一特殊短时期轨道尺度干湿变化规律和驱动机制研究[41-43]，制约了对干旱区低偏心率时期轨道尺度干湿变化规律和驱动机制的认识。例如，在柴达木盆地、贵德盆地、天水盆地、兰州盆地基于磁化率或Rb/Sr重建的古气候变化记录共同表明，晚渐新世（28.1～24.1 Ma）、早中新世（21.5～17.2 Ma）和晚中新世（约 14～7 Ma）降水变化以显著的10万年周期为主，这归因于偏心率调节的降水变化[11,16,41-43]。因此，为解答亚洲干旱区低偏心率时期气候变化如何响应地球轨道参数，亟需获取对气候变化更为敏感地区的高分辨率古气候记录。

柴达木盆地是重建低偏心率时期高分辨率轨道尺度气候变化历史的理想场所[42,44]。首先，柴达木盆地处于东亚季风区、亚洲内陆干旱区和青藏高原高寒区这三大自然区的交汇位置，对气候变化尤为敏感[45-46]。其次，柴达木盆地内发育巨厚的新生代河湖相地层，沉积速率高，较完整地记录了气候变化的详细信息[47-50]；位于其东北部的大红沟剖面地层出露较为连续，晚中新世时段平均沉积速率高达约34 cm/ka，提供了高分辨率环境干湿变化研究的良好地质载体[51-52]。本研究选取柴达木盆地大红沟剖面河湖相沉积地层为研究对象，运用频率磁化率指标，重建北半球无冰的晚中新世（9～12 Ma）高分辨率干湿变化历史，对比轨道参数，探究在低偏心率时期该区干湿变化规律、主导周期及其对地球轨道参数的响应机制。

1 研究区概况

柴达木盆地是中国西北内陆地区、青藏高原东北部一个封闭的山间断陷盆地，位于35～39°N、90～99°E（图1）。盆地被周围高大山系所环绕，西北、东北和南部分别为阿尔金山、祁连山-南山和昆仑山系，周围山脉海拔范围高达4 000～5 000 m，盆地内部海拔为2 700～3 000 m。气候方面，盆地当前大气环流模式主要受西风环流控制[53]，而盆内以高寒大陆性气候为主，终年干旱少雨，多大风天气；降水主要集中在夏季，年均降水量小于150 mm；蒸发量相对较强（年均2 000～3 000 mm）。构造上，盆地被东北部祁连山逆冲断裂、西北部阿尔金走滑断裂、南部东昆仑断裂和东部鄂拉山走滑断裂4个典型断裂带包围，在盆地的形成和演化中发挥着重要作用[54]。盆内部新生代河湖相地层分布广泛，大量沉积物质主要由周边山地剥蚀和风化作用提供，富含新生代环境演化过程的关键信息[47]。

大红沟剖面（37°31′N、95°09′E）位于柴达木盆地东北部，出露较为连续，全长6 200 m（图1）。剖面地层由老到新依次划分为路乐河、下干柴沟、上干柴沟、下油砂山、上油砂山和狮子沟组，主要为河湖相沉积。其中，本文研究的9～12 Ma段地层主要包括下干柴沟组（3 250～3 500 m）上部和上干柴沟组大部（3 500～4 400 m）。下干柴沟组以曲流河相沉积为主，岩性包括砂岩、泥质粉砂岩夹杂砾岩等，上干柴沟组则是以曲流河-滨湖相沉积为主，岩性以棕红色泥岩和黄绿色砂岩为主[55]。

图1 研究区及采样剖面位置图Fig.1 The study area showing sampling section location

目前对于大红沟剖面地层古地磁年代还存在三种不同的年代模式，第一种观点结合孢粉、介形虫、叶化石和介形虫等将剖面年代界定为古新世——晚中新世（约52～7 Ma）[48,50]；第二种观点结合新发现的红沟动物化石群将该剖面年代界定为晚渐新世——上新世（约25～5 Ma），与地层中发现的新近纪孢粉时代较为吻合[51]，第三种观点则将地层年代界定为早中新世——上新世（约21～5 Ma）[52]。其中，后面两种年代模式都是结合晚中新世哺乳动物化石群得到，年代主要差别在于路乐河组，而其他时段古地磁极性柱对比较好，说明剖面地层古地磁结果较为可靠[51-52]。此外，在第三种年代模式下的古气候研究发现，大红沟剖面降水在中中新世暖期增强，与区域乃至全球气候变化较为一致，也说明该年代模式是相对可靠的。因此本研究基于第三种年代模式进行高分辨率古气候重建。根据该年代模型已有的古地磁控制点[52]建立了9～12 Ma时期对应的时间标尺（4个古地磁年代控制点对应年代分别是 12.049、11.056、9.786和 9.105 Ma，结合控制点线性内插可得到详细的年代序列），为高分辨率轨道尺度干湿变化历史的重建提供了良好的年代基础[52]。

2 材料及方法

本剖面采样平均间隔1 m，共采集古环境样品999个，样品平均时间分辨率为3 ka。我们首先将散装样品在恒温烘箱（约40℃）中烘干，用研钵将样品研磨粉碎至无明显颗粒状（不破坏磁性矿物)，同时将粉末样品装入一个边长为2 cm的立方体（无磁性）塑料盒中压实固定并称重。然后，在远离电磁场干扰的环境中测试每个样品的磁化率值。测试仪器采用Barington MS2 磁化率仪，分别获得低频（470 Hz）和高频（4 700 Hz）磁化率（χlf和χhf），并重复测量2次，取其平均值。由平均高低频磁化率之间的差值计算（χfd=χlf-χhf）可得频率磁化率（χfd）。频率磁化率通常反映超顺磁（SP）和单畴（SSD）临界点附近颗粒含量的变化[56]，这些颗粒的形成主要受气候变化相关的成壤作用控制[57]。同时，以往基于黄土高原现代表层土壤中χfd与气候变化（降水和温度）关系的研究指出，χfd代表成壤或风化过程产生的超顺磁颗粒中亚铁磁性矿物含量的变化，与降水变化的相关性较高，χfd越高表明降水相对越多[58-60]。

对于河湖相沉积物，频率磁化率指标蕴含的信息更加复杂，可能受到其他因素的影响，比如沉积后的还原溶解作用及搬运过程中基岩碎屑物质加入的影响等。但是，大红沟剖面晚中新世以河流相和边缘湖相沉积为主，非封闭性湖泊，还原作用相对弱，磁性矿物发生溶解的可能性较小。另外，该剖面热退磁方法揭示出该序列沉积物中的主要载磁性矿物包括磁铁矿和赤铁矿[52]。而9～12 Ma期间大红沟剖面物质来源相对稳定[52]，说明碎屑物质来源变化不大，加入到沉积物中的基岩碎屑物质没有明显变化，不影响轨道周期的研究。因此我们认为在轨道尺度上，其他因素引起的误差有限，不足以影响对轨道周期变化的讨论。因此，本研究采用大红沟剖面沉积序列中的χfd记录来指代柴达木盆地东北部9～12 Ma的降水变化。

3 结果与讨论

将大红沟剖面频率磁化率log10χfd曲线与同时期地球轨道参数对比（图2），结果表明，在低偏心率和低进动（岁差）振幅时期（9.2～9.4 、 9.6～9.8 、11.2～11.4 Ma），大红沟剖面频率磁化率变化以明显的4万年周期为主，与倾角周期一致。同时，从相位关系上看，柴达木盆地低偏心率期（9.2～9.4 Ma和9.6～9.8 Ma），气候干湿变化与倾角变化具有相同的相位变化关系，倾角较大时，气候较湿润（图2）；然而，在另外一段低偏心率期（11.2～11.4 Ma），气候变化与倾角变化大致处于反相位关系（图2），造成差异的原因可能是干湿变化序列没有经过年代调谐，轨道尺度上年代控制存在一定误差。上述结果表明，低偏心率下柴达木盆地干湿变化受倾角变化控制。

图2 大红沟剖面9~12 Ma频率磁化率记录和地球轨道参数对比图a——c.地球轨道参数偏心率（红线）、岁差（绿线）和倾角（蓝线）曲线[2]， d.频率磁化率 (χfd) 4 万年周期滤波曲线， e.Log10χfd曲线，f.低偏心率时期频率磁化率小波分析图，白线及箭头指示频率0.025。Fig.2 Comparison in frequency magnetic susceptibility (χfd) recorded from the Dahonggou section with astronomical orbital parameter during 9 to 12 Maa-c.The earth’s orbital parameters of eccentricity (red line), precession (green line), and obliquity [2]; d.40 ka Gaussian bandpass filtered output of χfd;e.log10χfd curve; f.the wavelet transform of χfd during the period of low eccentricity, white line and arrows point at the frequency of 0.025.

前人研究发现，在晚中新世的低偏心率和低进动振幅时段，全球冰量变化和其他西风区气候变化也响应4万年周期性波动，与我们的结果一致，支持倾角驱动。例如，根据南海 ODP1146钻孔中晚中新世底栖有孔虫δ18O的记录，冰量变化在此时期（约9.2～9.7 Ma）以约4万年周期主导[6,61-62]。此外，有学者基于磁性地层和旋回地层学在西班牙东北部的普拉多剖面晚中新世（9.1～10.3 Ma）边缘湖相地层的研究发现，在偏心率最低值时（如9.6～9.8 Ma，约9.2～9.4 Ma）的地层旋回厚度明显高于全剖面平均岩层厚度，推测在进动振幅较小的时期，倾角变化可能对气候变化起主导作用[10]。

以往的研究表明，轨道参数倾角的周期性变化会影响南北半球太阳辐射的时空分布，倾角可以通过调节太阳辐射的温度梯度来影响极地冰量、西风急流和海陆热差，进一步控制亚洲气候变化[15,63-65]。

在前人研究基础上，结合柴达木盆地气候变化与倾角关系分析，推测倾角可能通过控制太阳辐射影响低偏心率时期（9.2～9.8 Ma和 9.6～9.8 Ma）大红沟剖面的干湿变化。主要原因如下：首先，倾角变化影响中纬度西风环流强度和西风位置，从而调节亚洲气候变化。倾角越大，夏季北半球径向太阳辐射和温度梯度越小，西风环流越弱，位置越偏北，降水越少[15]。本研究中柴达木盆地气候较为湿润，推测该区域受西风影响较小，且西风减弱和北撤也可能增强亚洲季风环流，使柴达木盆地气候比较湿润。其次，倾角通过改变南半球经向温度梯度来影响南极海洋冰盖动力学，进一步增强了全球气候对倾角驱动的敏感性[65]。南极冰盖的扩张可以通过增强跨赤道压力梯度和潜热释放来增强东亚夏季风[66]，而柴达木盆地的高降水在9.2～9.4 Ma和9.6～9.8 Ma对应于高倾角，不利于南极冰盖增长[67]，因此推测南极冰盖对该地区降水变化的影响较小。第三，模拟和地质记录表明，倾角也可能通过直接控制北方夏季太阳辐射梯度来改变海陆热力差，再通过调节高纬大陆低压系统和低纬海洋高压系统之间的热力差异，对亚洲季风气候产生影响。因此，高倾角可对应东亚夏季风降水增强，并对周边内陆地区产生持续影响[63,68]。近期的研究也已经指出，柴达木盆地在晚中新世构造和轨道尺度上可能受到东亚夏季风的影响[12,69]。因此，基于以上讨论，我们推测在晚中新世低偏心率时期，倾角通过调节太阳辐射梯度，影响亚洲夏季风，进而控制柴达木盆地轨道尺度的干湿变化。

需要指出的是，本文对于干湿变化的推断主要是基于频率磁化率单一指标记录，该指标可以指示沉积物中20～30 nm级磁性矿物的含量，这些细颗粒磁性矿物的含量主要是通过风化成壤作用产生的。因此，该指标在河湖相沉积物中比磁化率指标具有更加明确的环境干湿变化指示意义。然而在更大的构造尺度上该指标的变化还可能受到岩性变化和构造活动等其他因素的影响，因此本文主要讨论轨道尺度的环境变化。尽管如此，建议未来还应继续开展基于其他指标和不同尺度的干湿变化研究，进一步检验本文初步结论的正确性。

4 结论

本研究以青藏高原东北部柴达木盆地大红沟剖面晚中新世（12～9 Ma）河湖相沉积地层为研究对象，在已有的古地磁年代标尺基础上，运用频率磁化率指标分析低偏心率时期柴达木盆地高分辨率干湿变化规律和主导周期，并通过对比地球轨道参数记录，探讨了该时期干湿变化与地球轨道参数变化之间的响应关系。

研究表明，晚中新世低偏心率时期（9.2～9.4 Ma，9.6～9.8 Ma和 11.2～11.4 Ma）岁差 (进动) 振幅较小，柴达木盆地干湿变化具有明显的4万年周期，说明低偏心率时期该区域干湿变化可能主要受倾角控制。其中，9.4～9.2 Ma和9.8～9.6 Ma时段降水变化与倾角变化具有同相位关系，而11.2～11.4 Ma时段降水变化和倾角变化具有反相位关系，这可能受到年代误差影响。结合前人研究，我们推测低偏心率时期的倾角变化可能通过调节太阳辐射梯度来影响柴达木盆地干湿变化。