Cyclostratigraphy of the Ordovician Saergan Formation in the Tarim Basin

塔里木盆地中—上奥陶统萨尔干组是一套滞留盆地相沉积地层^[1⁃2],是我国著名的海相烃源岩之一^[3⁃5]。前人从岩石学、地球化学、地球物理学、层序地层学等方面对萨尔干组页岩进行广泛的研究^[1,6⁃11],并为萨尔干组黑色页岩的形成提出两种可能的沉积模式,分别为海侵模式和海退模式^[12]。海侵模式认为萨尔干组沉积于明显的海侵过程,海水逐渐淹没碳酸盐台地的同时,富含有机质的沉积物在台地上堆积形成黑色页岩；海退模式则倾向萨尔干组形成于海退过程,该过程使盆地与开阔海洋沟通变差,底层水循环减弱,导致富含有机质的沉积物形成黑色页岩^[12]。全球广泛分布的与萨尔干组同时期的富含有机质黑色页岩似乎证实了海侵模式^[13⁃14],但柯坪地区不同剖面间碳同位素的变化趋势与该模式相矛盾^[12]。由大湾沟组灰岩到萨尔干组黑色页岩的岩性突变被认为是海平面下降使得盆地与开阔海洋沟通变差,底层水缺氧导致,这是海退模式的主要证据。但是,上奥陶统最底部Nemagraptus gracilis笔石带并未在柯坪地区四石场剖面识别出,关键生物带的缺失不支持这种解释^[16]。因此,岩性、生物和化学地层学等结果的不吻合,使得准确判断萨尔干组沉积模式更为困难^[12]。高精度的地质年代可以约束萨尔干组的沉积时限,实现高精度的区域地层对比,为萨尔干组的沉积模式提供更高分辨率的年代学证据。

旋回地层学通过识别和解释地层中保存的由地球轨道力（如偏心率、斜率和岁差旋回）控制的（准）周期性沉积旋回信号,建立精度可达20~400 ka的天文年代标尺,从而精确约束各种地质事件及地质过程的绝对年代和持续时间^[17]。前人在中—晚奥陶世旋回地层学研究中取得显著成果^[18⁃19],在此期间发生的气候变冷和碳同位素漂移等重大事件均被证明可能与地球轨道旋回控制的古气候变化有关^[19⁃20]。目前,~1.2 Ma超长斜率旋回已经被证明是奥陶纪气候和冰川—海平面变化的重要驱动因素之一^[18⁃19]。

前人已经对塔里木盆地大湾沟剖面萨尔干组进行详细的旋回地层学研究,建立的高精度天文年代标尺,为中—晚奥陶世提供了高分辨率的年代地层框架^[18]。本文对塔里木盆地苏巴什沟剖面中—上奥陶统萨尔干组进行系统样品采集并进行质量磁化率测试,利用获得的数据开展详细的旋回地层学研究,识别剖面地层中记录的米兰科维奇旋回,利用地质历史时期稳定的405 ka长偏心率旋回建立高精度的天文年代标尺,旨在：（1）实现40万年偏心率时间尺度上萨尔干组的地层对比；（2）确定萨尔干组的沉积模式；（3）探讨萨尔干组沉积时期海平面变化驱动机制。

塔里木盆地位于中国的西北部,面积约占5.6×10⁵ km²,是世界上最大的内陆盆地之一。其内部主要被塔克拉玛干沙漠所覆盖,露头仅见于西北部和东北部地区（图1b）。塔里木盆地发源于前震旦纪陆壳基底^[10],古生代位于冈瓦纳大陆和劳伦大陆之间的古特提斯洋^[22]。中—晚奥陶世,塔里木盆地位于北半球冈瓦纳大陆西北缘,毗邻西伯利亚、华南、华北以及澳大利亚板块^[23⁃24]（图1a）。前人将塔里木盆地分为十二个一级构造单元,分别为三个隆起（塔北隆起、中央隆起和塔南隆起）、五个坳陷（库车坳陷、北部坳陷、西南坳陷、东南坳陷和塘古孜巴斯坳陷）和四个边缘断隆构造（柯坪断隆、库鲁克塔格断隆、铁克力克断隆和阿尔金山断隆）^[25]。研究剖面所在的柯坪地区位于中央隆起的西部,奥陶纪属于斜坡相沉积环境^[1,12]。该地区沉积层序所对应的海进—海退旋回结构在塔里木盆地范围内可以很好地对比,并且与全球海平面变化曲线类似,反映了该地区沉积地层的发育受到全球海平面变化的控制^[10]。

塔里木盆地萨尔干组出露厚度为6.0~13.0 m,岩性为夹少量灰岩薄层或透镜体的黑色碳质泥页岩,局部含有少量泥质条带^[15]。在柯坪地区沿北东—南西向出露于苏巴什沟以东、四石厂以西、沙井子—三岔口以北以及乌什以南等地^[26]。萨尔干组产出四个笔石带,从下到上依次为Pterograptus elegans带、Didymograptus murchisoni带、Dicellograptus vagus带和Nemagraptus gracilis带^[16],其时代为奥陶纪达瑞威尔晚期到桑比早期（图1c）。本文研究的苏巴什沟地区（地理坐标为40°35′27″ N,78°57′20″ E）是所有产出萨尔干组沉积序列中最西端的剖面^[1]。苏巴什沟剖面从底部到顶部依次出露大湾沟组、萨尔干组和坎岭组。其中萨尔干组厚度为9.0 m（图2a）,岩性主要为黑色页岩夹薄层灰岩透镜体,与下伏大湾沟组薄层灰岩和上覆坎岭组灰岩界线清晰且均为整合接触,夹于萨尔干组上部和下部的薄层灰岩差异显著,下部可见泥粒灰岩或颗粒灰岩,上部主要为含生屑泥质灰岩^[1]。N. gracilis带首次出现层位（first appearance datum,FAD）为桑比阶的底界^[12],苏巴什沟剖面N. gracilis带FAD大致在萨尔干组6.1 m处^[16]（图2a）。

Figure 2.  Lithostratigraphy, biostratigraphy, cyclostratigraphy, and magnetic susceptibility data series

大湾沟剖面（40°43′16″ N,79°31′51″ E）在2000年被推荐为上奥陶统底界全球辅助层型剖面^[14],是目前萨尔干组实测厚度最大、研究程度最高的剖面^[1]（图1b）。该剖面位于苏巴什沟剖面东北约50 km处的柯坪县印干村西北侧,该剖面的萨尔干组实测厚度为12.0 m,由黑色页岩夹灰黑色薄层或透镜状泥屑灰岩组成。沿大湾沟剖面,海相地层从大坪阶到凯迪阶中部连续出露,依次为大湾沟组、萨尔干组、坎岭组和印干组（图1c）。在大湾沟剖面的萨尔干组黑色页岩中还可见草莓状黄铁矿^[1]。大湾沟剖面萨尔干组的旋回地层学结果已经由前人发表^[18],作为本研究对比依据。

磁化率（magnetic susceptibility,MS）是衡量一个物质被磁化难易程度的无量纲岩石磁学参数。通常情况下,经测试的样品磁化率为岩石中各种磁性矿物（顺磁性、反磁性和铁磁性矿物）的磁化率值总和^[27]。磁化率值的变化可以作为陆源输入的指标^[28⁃29]。在海相地层中,当日照量处于较高（低）阶段时,气温升高（降低）,降水增加（减少）,地表径流增强（减弱）,因此岩石磁化率值较高（低）。磁化率的变化还能够记录冰期—间冰期旋回以及更高频的气候振荡等^[30⁃31]。

苏巴什沟剖面萨尔干组厚9.0 m,野外共计采集701个样品。马璐等^[1]对大湾沟剖面和苏巴什沟剖面萨尔干组灰岩夹层进行沉积微相分析,结果表明这些灰岩形成于斜坡上部的浅水灰岩沉积区,经重力流作用后被运移到此处,属于异源钙质浊积岩,因此本文将苏巴什沟剖面萨尔干组0.99~1.02 m、1.89~1.92 m、2.91~2.94 m、7.80~7.82 m、8.40~8.42 m采集的灰岩数据移除,从而排除异源钙质浊积岩的影响（图2）。合理的采样间距可以确保采集的样品数据记录所需天文周期的最高频率信号^[32]。根据Nemagraptus gracilis等笔石生物年代框架的约束^[15]和对大湾沟剖面萨尔干组页岩0.19 cm/ka沉积速率的估计^[18],对苏巴什沟剖面萨尔干组以1~2 cm间距采样。将浊积岩移除后的沉积序列共8.9 m,对野外采集样品进行碎样装塑料方盒,使用中国地质大学（北京）古地磁实验室MFK1-FA卡帕桥磁化率仪测得数据。大湾沟剖面的磁化率数据源于Fang et al.^[18]。

剖面移除灰岩后的数据序列为8.9 m,共测得691组数据,对其进行间隔为0.013 m的线性插值。使用Matlab软件中的smooth.m函数,采用“moving”方法计算并去除大于1.3 m的长周期和趋势,保留高频旋回信号^[32]。为了识别米兰科维奇信号,对插值、平滑处理后的等间距深度域序列进行多窗口频谱（Multi-Taper Method,MTM）分析^[33],同时对预处理后的序列进行滑动频谱分析来确定不同深度的周期分布特征^[27]。利用稳定的405 ka周期建立年龄模型,在Matlab中通过depthtotime.m函数^[27]将磁化率深度域序列转换为时间序列。对获得的时间域序列进行频谱分析和带通滤波分析的天文旋回检验。除滤波分析在R^[34]的“Astrochron”软件包中的“Bandpass”函数^[35]下执行外,其余分析均在Matlab中完成。

奥陶纪虽然没有可靠的天文解决方案,但利用地质历史时期稳定的405 ka长偏心率周期^[36⁃37],可以实现天文校准。Fang et al.^[38]在华北板块上奥陶统平凉组中,识别出的斜率周期为30.6 ka,岁差周期为19.6 ka和16.3 ka。塔里木盆地大湾沟剖面的中—晚奥陶世萨尔干组记录的短偏心率、斜率和岁差周期分别为100 ka、31 ka和~20 ka^[18]。前人计算出,458.4 Ma时地球斜率周期为33.2±3.8 ka,岁差周期为20.8±1.5 ka和17.0±1.1 ka^[39]。Svensen et al.^[40]从挪威奥斯陆地区获得斜率周期为30.3 ka,岁差周期为19.2 ka和16.3 ka,并得到了高精度锆石U-Pb年代学的支持。前人基于405 ka长偏心率校准获得的奥陶纪轨道周期与估算的天文轨道周期基本一致（表1）,因此本文采用405 ka长偏心率周期进行校准。

苏巴什沟剖面萨尔干组磁化率值整体在（1~8）×10^-8 m³/kg之间变化,平均值为4.36×10^-8 m³/kg。在0~3.5 m厚度范围内磁化率值先迅速降低,此后磁化率值缓慢下降,波动范围较小,整体处于低值状态；在3.5~5.8 m厚度范围内,磁化率值逐步升高后保持稳定,整体为趋势平稳的高值；5.8~9.0 m处磁化率值较前一段有所降低,磁化率值的变化幅度也较明显（图2a）。大湾沟剖面萨尔干组磁化率值的变化范围为（1~12）×10^-6 SI,平均值为5.5×10^-6 SI（图2b）。大湾沟剖面磁化率的变化趋势整体与苏巴什沟剖面的磁化率序列相似。苏巴什沟剖面在0~3.5 m处与大湾沟剖面2.4~6.0 m处磁化率均呈降低趋势,苏巴什沟剖面在3.5~5.8 m处与大湾沟剖面6.0~8.8 m处磁化率值均迅速达到高值并呈波动变化,在N.gracilis笔石的首现层位处,磁化率均为低值,以及苏巴石沟剖面6.0~9.0 m与大湾沟剖面9.0~10.5 m均具有两处明显峰值（图2）。苏巴什沟剖面与大湾沟剖面^[18]具有相似的变化趋势,磁化率的高值对应了黑色页岩,磁化率的低值对应了夹杂薄层灰岩的页岩。

对苏巴什沟剖面萨尔干组深度域磁化率序列进行频谱分析和滑动频谱分析,识别到多个高置信优势峰谱,分别在0.741~1.283 m、0.204~0.267 m、0.057~0.071 m和0.035~0.039 m具有优势峰值,其置信度均大于95%,在滑动窗口频谱分析中也可以找到对应的能量信号（图3a,b）。以上显著谱峰与Fang et al.^[18]在大湾沟剖面识别的0.7~1.0 m、0.20~0.26 m、0.042~0.055 m和0.030~0.034 m的峰值接近。其比值约为21.14∶7.71∶5.71∶1.86∶1.11∶1,与前人预测的奥陶纪地球轨道参数周期吻合（表1）,这些峰值可能分别代表长偏心率（E）、短偏心率（e）、斜率（O）和岁差（P）周期（图3a,b）。这表明苏巴什沟剖面萨尔干组的沉积可能受到了天文轨道力的控制。

Figure 3.  Spectrum analysis and evolutive Fast Fourier Transform

根据MTM频谱分析的结果,对苏巴什沟剖面深度域磁化率数据序列进行带通滤波,以提取地层中记录的长偏心率信号。苏巴什沟剖面滤波带宽为1.0±0.2旋回/m,大湾沟剖面滤波带宽为1.2±0.3旋回/m。滤波结果共识别出10个长偏心率旋回（图2a）。利用405 ka年龄模型进行天文调谐,将磁化率的深度域序列转化为时间域序列。为了检验旋回周期识别的合理性以及天文校准的可靠性,将校准后得到的磁化率时间序列进行频谱分析。苏巴什沟剖面时间域的频谱分析结果显示在405 ka、100 ka、31.2 ka、30.6 ka、21.8 ka、20.6 ka和20.0 ka处有显著波峰（图3c,d）。因此,根据405 ka年龄模型的天文校准,将剖面深度域序列转化为时间域序列,建立了苏巴什沟剖面萨尔干组~3.9 Ma的“浮动”天文年代标尺。

柯坪地区和北部坳陷区阿瓦提断陷的萨尔干组是在快速海进背景下,发育于滞留盆地相的一套欠补偿的凝缩段沉积^[1⁃2]。Fang et al.^[18]依据N.gracilis笔石带的FAD对应于桑比阶底界（458.4 Ma）,将大湾沟剖面的浮动天文年代标尺校准到国际地质年代标尺上。因此,将N.gracilis笔石带的首现层位作为基准,将两个剖面的天文年代标尺共同校准到国际地质年代表上进行剖面间对比（图4）。对比结果显示,不论是405 ka长偏心率周期的良好对应关系还是磁化率数据序列的变化趋势,均证实了两个剖面间萨尔干组高精度旋回地层对比的可靠性。

Figure 4.  The constructed astronomical time scale and 405 ka filtering outputs

根据校准结果显示,大湾沟剖面在462.5~461.5 Ma的磁化率值为高值,直到~460.1 Ma两个剖面的磁化率值一直呈缓慢下降趋势,振幅较小。此后两个剖面的磁化率值均呈上升状态,并且呈波动状变化（图4）。磁化率的降低已被证明是由于日照量增大,海平面升高所导致^[43]。值得注意的是,萨尔干组上、下部夹层灰岩透镜体不同的微相也指示了海平面的上升。由西南至东北的方向,苏巴什沟和大湾沟剖面的萨尔干组下部的灰岩薄层微相中可见粗砂屑颗粒堆积层,指示海侵前相邻的浅海碳酸盐岩台地被机械搬运异地快速沉积成岩,上部的含生屑泥状灰岩可能代表海平面上升后在滞留盆地内的正常沉积^[1]。萨尔干组沿柯坪—阿克苏一线经历了从无到有,从薄到厚的变化,同时下部灰岩透镜体夹层数量呈现大湾沟剖面多,苏巴什沟少的特点。夹层数量呈现大湾沟剖面多,苏巴什沟少的特点。前人在大湾沟的黑色页岩中观察到了草莓状黄铁矿^[18],而在苏巴什沟剖面目前尚未发现,说明大湾沟剖面萨尔干组沉积于还原环境,指示大湾沟剖面相对苏巴什沟剖面萨尔干组更靠近下部,水体更深。

综上所述,本文建立了中—晚奥陶世大湾沟剖面和苏巴什沟剖面的沉积古地理相对位置,并根据海平面升降关系,确定大湾沟剖面在~462.5 Ma时先开始沉积萨尔干组（图5a）,随着海侵的进行,海平面不断上升,苏巴什沟剖面直到~461.5 Ma时萨尔干组才开始沉积（图5b）,晚于大湾沟剖面~1.03 Ma。此后海平面短时间内虽略有下降,但整体仍为上升趋势,两个剖面萨尔干组均连续沉积至~457.7 Ma（图5c）。上述分析表明大湾沟和苏巴什沟剖面的萨尔干组在明显的海侵中先后接受沉积（图5）,支持了前人所提出的海侵模式^[12]。

Figure 5.  The sedimentary model in the Kalpin area during the Middle⁃Upper Ordovician

生物地层学、层序地层学以及气候模拟等方面的大量研究均证明了中奥陶世冰川的存在^[44⁃49]。塔里木盆地苏巴什沟剖面萨尔干组黑色页岩是一套连续的中奥陶世海相地层沉积,为研究斜率超长周期对该时期冰川—海平面变化的驱动机制解释提供了理想的材料。大湾沟和苏巴什沟两个剖面斜率周期的振幅调制记录（图6d,e）显示~30.6 ka斜率旋回的振幅变化具有~1.2 Ma的主要周期。~462.5 Ma大湾沟剖面萨尔干组开始沉积时期,此阶段磁化率值整体偏高,对应了~1.2 Ma斜率周期的振幅增大阶段,恰好与萨尔干组沉积处于海平面上升时期相吻合（图6a,b）；随后直到~460.3 Ma时,磁化率值一直呈现降低的趋势,对应了~1.2 Ma斜率周期的振幅减小阶段。~459.9 Ma萨尔干组的磁化率值升高,斜率振幅增大,此时对应海平面上升时期。~1.2 Ma斜率长周期与全球海平面变化具有共同的相位关系,~1.2 Ma斜率长周期的最低值对应了海平面的低位。在~1.2 Ma斜率周期的低值阶段,地轴斜率长期处于较弱的状态,振幅波动较小,季节性反差不大,夏季较为凉爽,冬季较为温和,利于高纬地区冰川逐年积累,最终导致海平面的降低^[19,50]；相反,当~1.2 Ma斜率周期处于高值阶段时,振幅波动较大,季节性增强,炎热的夏季和寒冷的冬季不利于高纬地区冰川的积累,从而使得海平面升高。综上所述,萨尔干组沉积时期的海平面变化受到了~1.2 Ma超长斜率旋回的控制。

Figure 6.  Sea⁃level change and obliquity amplitude modulation

（1） 塔里木盆地苏巴什沟剖面地层记录了完整的米兰科维奇信号,其中包括405 ka长偏心率、100 ka短偏心率、30.6~31.2 ka斜率以及20.0~21.8 ka岁差旋回。

（2） 利用405 ka长偏心率周期进行天文年代校准,建立了苏巴什沟剖面萨尔干组的天文年代标尺,以Nemagraptus gracilis笔石带的首现层位作为等时界面,将大湾沟剖面和苏巴什沟剖面萨尔干组在405 ka长偏心率时间尺度上进行了高精度对比,证明了中—晚奥陶世时期萨尔干组的沉积模式为海侵模式。

（3） 中奥陶世海平面变化或受控于~1.2 Ma超长斜率旋回。