谢琴, 陈家胜, 刘晓静, 刘秀铭. (2022). 甘肃张掖白垩纪早Aptian期黑色页岩揭示的气候岁差驱动及其意义* [J]. 古地理学报, 24(2): 375-388.
Xie Qin, Chen Jia-Sheng, Liu Xiao-Jing, Liu Xiu-Ming. (2022). Detection of climatic precession variability of black shale in the early Aptian of Cretaceous in Zhangye of Gansu Province and its significance[J]. Journal Of Palaeogeography, 24(2): 375-388.   
Doi: 10.7605/gdlxb.2022.02.026
Permissions
Copyright©2022, 《古地理学报》编辑部
版权所有.《古地理学报》编辑部
甘肃张掖白垩纪早Aptian期黑色页岩揭示的气候岁差驱动及其意义*
谢琴1,2, 陈家胜1,2, 刘晓静1,2, 刘秀铭1,2,3
1 福建师范大学地理研究所,福建福州 350007
2 福建师范大学湿润亚热带生态—地理过程教育部重点实验室,福建福州 350007
3 麦考瑞大学环境与地理系,澳大利亚悉尼 NSW2109

通讯作者简介 陈家胜,男,1985年生,福建师范大学地理科学学院助理研究员,硕士生导师,主要从事环境磁学、土壤与环境等研究。E-mail: chenjsh04@fjnu.edu.cn。

第一作者简介 谢琴,女,1996年生,福建师范大学硕士研究生,主要从事白垩纪古气候研究。E-mail: xieqincwnu@163.com。

收稿日期:2021-06-12 修回日期:2021-11-16
基金资助:*福建省公益类科研院所专项项目(编号: 2019R1002-5)、国家自然科学基金项目(编号: 41602190)、福建师范大学创新团队项目(编号: IRTL1705)联合资助
摘要

白垩纪中期是温室气候的典型代表,在该时期地球经历了深刻的环境变化。有研究表明白垩纪中期欧洲低纬(约 20°~30°N)地区海相地层的黑色页岩呈现斜率信号,同一时期同样纬度的中国地区以陆相沉积为主,但陆相低纬黑色页岩的轨道尺度周期研究较少。因此本研究选择甘肃张掖南台子白垩纪早 Aptian期下沟组黑色页岩(古纬度约 23°N)作为研究对象,以 CaCO3含量作为高分辨率古气候替代性指标,利用平均频谱拟合差和年代标尺优化法分析了南台子早 Aptian期下沟组黑色页岩段的轨道周期变化,并利用频谱分析方法对调谐后的时间域序列进行分析。结果显示南台子下沟组黑色页岩段以岁差周期为主导,与欧洲低纬地区相关研究显著不同。这一气候变化的探索不仅提供了中国西北地区白垩纪轨道尺度气候变化的新证据,而且也有助于更全面认识白垩纪中期低纬地区气候变化的驱动机制。

关键词: 岁差; 碳酸钙; 黑色页岩; TimeOpt; ASM
中图分类号:P532 文献标志码:A 文章编号:1671-1505(2022)02-0375-14
Detection of climatic precession variability of black shale in the early Aptian of Cretaceous in Zhangye of Gansu Province and its significance
Xie Qin1,2, Chen Jia-Sheng1,2, Liu Xiao-Jing1,2, Liu Xiu-Ming1,2,3
1 Institute of Geography,Fujian Normal University,Fuzhou 350007,China
2 Key Laboratory of Humid Sub-tropical Eco-geographical Process of Ministry of Education, Fujian Normal University,Fuzhou 350007,China
3 Department of Environment and Geography,Macquarie University,Sydney NSW2109,Australia

About the corresponding author Chen Jia-Sheng,born in 1985,is an assistant researcher and master tutor of College of Geography Science,Fujian Normal University. He is currently engaged in environmental magnetism,soil and environment research and so on. E-mail: chenjsh04@fjnu.edu.cn.

About the first author Xie Qin,born in 1996,is a master degree candidate of the Fujian Normal University. She is currently engaged in Cretaceous paleoclimate research. E-mail: xieqincwnu@163.com.

Fund:Special Project of Fujian Public Welfare Research Institute(No.2019R1002-5),the National Natural Science Foundation of China(No.41602190),and the Innovation Research Team Fund of Fujian Normal University(No. IRTL1705)
Abstract

The earth experienced profound environmental changes at the mid-Cretaceous when a typical representative of the greenhouse climate occurred. In the mid-Cretaceous,studies have shown that black shale in low-latitude(ca.20°~30°N)marine strata in Europe records presents obliquity signals,continental deposits are prevalent in inland China at the same latitude and during the same period,but there are few studies on the orbital-scale periodicity of continental low-latitude black shales. In this study,CaCO3 content is used as an alternative high-resolution indicator of paleoclimate,and average spectral misfit(ASM)and time scale optimization(TimeOpt)analyses are used to explore the orbital periodicity of the black shale in the Xiagou Formation of the early Aptian of Cretaceous in Nantaizi (palaeoatitude of ca 23°N),Zhangye of Gansu Province, China. The result shows that the climate is dominated by precession and eccentricity cycles,which is significantly different from the result in Europe. The exploration of this climate change not only provides new evidence for the cretaceous orbital-scale climate change in the Northwest China,but also contributes to a more comprehensive understanding of the driving mechanism of climate change in the mid-Cretaceous.

Key words: precession; CaCO3; black shale; TimeOpt; ASM

开放科学(资源服务)标志码(OSID)

1 概述

地球轨道参数的周期性变化影响不同纬度和时间尺度上的太阳辐射量的周期性变化, 进而导致气候的周期性变化, 因此在地球表层沉积系统中会记录着地球轨道的周期信号(米兰科维奇旋回)(Milankovitch, 1941; Hinnov, 2000)。将地质记录中的米兰科维奇旋回信号与标准天文理论曲线对比, 可以更准确地认识地质记录的旋回与天文周期之间的联系, 确定气候的主导周期以及地层的相对年代(Hinnov, 2000; Hinnov, 2013), 对认识地质演化历史具有重要意义。

根据Ruddiman(2001)所绘制的近30万年以来不同纬度和不同季节太阳辐射量变化图可知夏季高纬和低纬的太阳辐射量主要以岁差周期为主; 汪品先等(2018)所绘制的2 Ma以来的高低纬度月均太阳辐射量的频谱分析表明, 低纬(15° N)地区显示强烈的岁差周期; Short等(1991)通过数值模拟表明中低纬度地区受到~20 ka岁差周期影响大, 岁差周期受到偏心率的调控作用明显。以上研究表明低纬气候过程有强烈的岁差周期, 岁差受偏心率周期调控。

第四纪的气候驱动机制研究较为成熟。低纬地区的非洲季风是受岁差影响的典型例子, 如热带大西洋ODP659站点(18° 05'N, 21° 02'W)风尘含量交叉频谱分析的结果表明北非季风呈现明显的岁差和偏心率周期(Tiedemann et al., 1994); 被吹到热带大西洋的非洲干涸湖泊的硅藻记录着~23 ka的岁差周期(Ruddiman, 2001)。中低纬区域的亚洲季风区石笋 δ 18O 数据记录都表明主导周期为岁差周期, 如: 中国新疆科桑洞(42° 52'N, 81° 45'E)(Cheng et al., 2012)、湖北三宝洞(31° 40'N, 110° 26'E)和林竹洞(31° 31'N, 110° 19'E)(Cheng et al., 2016)、江苏葫芦洞(32° 30'N, 119° 10'E)(Cheng et al., 2009)、贵州董哥洞(25° 17'N, 108° 5'E)(Wang et al., 2005)、云南小白龙洞(24° 12'N, 103° 21'E)(Cai et al., 2015), 印度北部Bittoo洞(30° 47'N, 77° 46'E)(Kathayat et al., 2016)等。同时, 中国广东省湖光岩玛珥湖(29° 9'N, 110° 17'E)的孢粉记录(Wang et al., 2007)和赤道印度洋MD900963站点(5° N, 73° E)的海水表层生产力频谱分析都显示出明显的岁差周期(Beaufort et al., 1997)。从上述研究来看, ~20 ka岁差周期确实对低纬地区影响较大。

白垩纪中期(Aptian-Turonian: 125— 89 Ma)是地质历史中的极热期, 深海与极地表层的极端温度(15~20 ℃)(Huber et al., 2002; Jenkyns et al., 2004; Friedrich et al., 2012)以及热带海水表层温度(33 ℃到35~36 ℃)(Wilson et al., 2002; Forster et al., 2007; Bornemann et al., 2008)都大大超过了现代值; 大气CO2浓度是工业革命前的4~10倍(Berner and Kothavala, 2001); 赤道到两极的温度梯度比现在更低(Littler et al., 2011)。

白垩纪中期全球广泛发育了黑色页岩(Schlanger and Jenkyns, 1976)。Aptian期的黑色页岩的沉积尤为显著, 且在欧洲海相地层中对其研究程度高。在Aptian早期, 意大利中部Piobbico岩心(古纬度约20° N; Tiraboschi et al., 2009; 图 1)的黑色页岩层段和邻近地层记录的主导周期为斜率周期(Huang et al., 2010); 意大利北部的Cismon APTICORE(古纬度约20° N; 图 1)黑色页岩层的频谱分析显示主导周期为短偏心率、斜率和岁差周期(Li et al., 2008); 墨西哥东北部Santa Rosa峡谷剖面(古纬度约20° ~30° N, 图 1)的磁学参数频谱分析结果均表明短偏心率与斜率周期信号强烈, 岁差周期较弱(Li et al., 2008); 北大西洋DSDP398站点(古纬度介于20° ~30° N之间; 图 1)的磁学参数频谱分析结果表明偏心率周期较强(Li et al., 2008); 古纬度为约30° N的法国南部Vocontian盆地(图 1)岩心的自然伽马射线数据EHA分析结果表明岁差、斜率和偏心率都存在, 但偏心率和斜率周期信号较强(Beil et al., 2020)。以上结果都是基于欧洲海相地层的研究, 并且Aptian期的低纬度海相地层(20° ~30° N)的主导周期多数以斜率周期为主。

图 1 120 Ma全球板块分布以及相关研究区域的地理位置
该图源自于网站(www.odsn.de)。1. 意大利Piobbico岩心; 2. Cismon APTICORE岩心; 3. 墨西哥Santa Rosa剖面; 4. 北大西洋DSDP398站点; 5. 法国南部的Vocontian盆地。Sc和Sm分别代表白垩纪中期四川盆地(潘忠习等, 1999)和思茅盆地 (Jiang et al., 2001)的位置。五角星代表研究区位置Fig.1 Map of global plate distribution at 120 Ma and geographical location of relevant study areas

欧洲海相白垩纪地层出露广泛, 研究程度较高, 亚洲的中国区域陆相白垩系分布较广, 类型复杂多样(Xi et al., 2019)。中国西北地区陆相白垩纪地层在六盘山盆地、酒泉盆地等区域广泛分布(Xi et al., 2019)。酒泉盆地研究程度较高, 并已经将盆地内下沟组黑色页岩段的年代确定为早Aptian期(Li et al., 2013; Suarez et al., 2013)。古纬度约为23° N(Frost et al., 1995)的西北甘肃张掖南台子地区也发育很好的白垩纪陆相沉积, 地层与酒泉盆地划分一致, 存在相同时间段相似的黑色页岩(杨雨, 1997)。理论上说, 其气候过程极有可能呈现出强烈的岁差周期, 如果确实如此, 这与欧洲白垩纪中期的低纬海相地层所呈现的斜率周期不同。因此有必要深入对南台子下沟组黑色页岩段气候周期开展实验与分析, 补充认识白垩纪中期低纬地区陆相沉积的轨道周期变化规律。

2 地质背景

白垩纪甘肃位于走廊潮水盆地, 发育内陆开阔盆地湖相泥质— 泥灰质沉积, 该地区的气候以潮湿亚热带— 热带气候为主(张兰生和方修琦, 2012)。早期将巨厚黑色页岩划为赤金堡组的顶部, 作为赤金堡组与下沟组的分界线(杨雨, 1997)。Li等(2013)将黑色页岩划为下沟组的底部, 根据锆石U-Pb定年、植物化石的年龄参考以及与OAE1a事件相联系, 共同确定黑色页岩的年代为早Aptian期(124— 125 Ma)。另外Suarze等(2013)在下沟组湖相地层通过碳同位素对比也确定黑色页岩的年代为早Aptian期。因此河西走廊下沟组黑色页岩年代应属于白垩纪早 Aptian期。

研究区地处甘肃省西北部、河西走廊中段、祁连山北麓(图 2-a)。研究剖面位于临泽县倪家营乡南台子村彩色丘陵景区西侧(图 2-b)。彩色丘陵是由下沟组和中沟组所组成(张荷生和崔振卿, 2007), 黑色页岩位于下沟组底部(图 2-c)。为了便于采样, 在黑色页岩的平行剖面开展研究, 将该剖面命名为冰沟剖面(图 2-d)。

图 2 甘肃张掖南台子下沟组黑色页岩段Google Earth图和野外实景照片
a— 研究区及邻近地区的地理图, 浅蓝色虚线代表河西走廊位置, 红色五角星代表研究区位置(修改自Liu et al., 2017); b— 张掖南台子彩色丘陵与下沟组黑色页岩段的Google Earth图, 实际采样点坐标:38° 57'N, 99° 59'E; c— b图中巨厚黑色页岩地层实景; d— 冰沟剖面黑色页岩实景, Ⅱ 段为研究层位Fig.2 Google Earth map and field photos of the black shale of Xiagou Formation in Nantaizi, Zhangye of Gansu Province, China

冰沟剖面发育了4段黑色页岩, 这4段页岩之间间隔相似, 厚度相当, 与黑色页岩段相邻的黄色段主要由砾石和砂岩组成。黑色页岩Ⅱ 段(38° 57'N, 99° 59'E)出露更加完整且更具代表性(图 2-d), 因此在黑色页岩Ⅱ 段进行采样。黑色页岩Ⅱ 段总厚29.85 m, 按照岩性主要分为3段(图 3):剖面下段(0~11.75 m)主要为黑色和灰色的泥岩和页岩互层; 中段(11.75~17.35 m)主要为黄色、棕黄色、暗灰色粉砂岩, 偶夹灰色、黑色泥页岩; 上段(17.35~29.85 m)主要为黑色和灰色泥页岩, 偶夹黄色粉砂条带。本次对黑色页岩Ⅱ 段从下至上进行采样(图 3), 首先去除表层风化物, 在黑色和灰色区域以0.05 m等间距连续采样, 黄色粉砂岩处根据岩性变化将采样间距设为0.1~0.3 m, 在Ⅱ 段共采样523个。

图 3 甘肃张掖南台子下沟组黑色页岩Ⅱ 段的岩性地层
a— 黑色页岩Ⅱ 段的下段; b— 黑色页岩Ⅱ 段的中段与上段; c— 泥页岩; d— 页岩; e— 砂岩; f— 泥岩; g— 黑色页岩Ⅱ 段岩性柱Fig.3 Lithologic stratigraphy of the black shale section Ⅱ of Xiagou Formation in Nantaizi, Zhangye of Gansu Province, China

3 方法
3.1 碳酸钙含量测定

碳酸钙(CaCO3)含量采用气量法测定, 其实验步骤如下: 首先将0.5 g样品放入烘箱(< 65 ℃)烘干, 用玛瑙研钵将样品研磨至均匀; 其次使用碳酸钙简易装置测试, 加过量的浓度为10%的稀盐酸于待测样品中, 等待盐酸与样品中的CaCO3完全反应; 最后读取CaCO3与盐酸反应后产生的二氧化碳气压值, 根据该气压值进行换算就可得到样品的CaCO3含量。本研究将全部样品进行CaCO3含量测试。

3.2 总有机碳和总有机氮含量测定

总有机碳(TOC)和总有机氮(TN)含量测定的实验步骤如下(Suarez et al., 2013; Zhang et al., 2016): 首先将1~2 g样品研磨成粉状, 加入过量0.5 mol的盐酸去除无机碳。其次将样品与过量盐酸反应直至把无机碳全部去除, 再用去离子水将样品洗涤多次直至pH值达到中性。最后将样品放置在45~50 ℃烘箱干燥后, 重新研磨至均质状。本研究以0.25 m间隔选取153个样品测试TOC。土壤总有机碳和总有机氮使用碳氮元素分析仪(Elementar Vario EL Ⅲ , 德国)测量。以上实验在福建师范大学湿润亚热带山地生态国家重点实验室完成。

3.3 旋回地层学研究方法

3.3.1 频谱分析(MTM和EHA)

Multi-Taper Method(MTM)分析是将深度域信号转化为频率信号, 并识别数据中的周期信号(Thomson, 1982)。但MTM仅能反映一段深度(时间)内平均谱结构, 不能反应频率随深度(时间)的变化情况(Meyers et al., 2001)。Evolutive Harmonic Analysis(EHA)是在数据窗口上使用滑动MTM方法来分析深度(时间)域数据序列, 可得到动态空间频率的谱结构, 该方法可以识别地层数据序列中的主导频率周期的沉积速率变化情况(Meyers et al., 2001)。本研究有关MTM和EHA的分析全部根据R软件“ astrochron” 包下的mtm和eha方法计算。

3.3.2 理论天文轨道周期

本研究的天文参数模型来自La2004(Laskar et al., 2004)天文解决方案(提供过去250 Ma到未来的250 Ma的参数模型), La2004方案已经考虑地月的潮汐耗散与太阳系的混沌运动等影响因素, 可建立不确定性为0.2%的浮动天文年代标尺。作者利用La2004天文解决方案计算早 Aptian(120— 125 Ma)时期的偏心率、斜率和岁差周期变化理论值(图 4), 经过频谱分析后得到5个偏心率周期(400.08 ka、131.61 ka、123.48 ka、99.03 ka、94.36 ka)、3个斜率周期(46.74 ka、36.77 ka、35.72 ka)和3个岁差周期(22.23 ka、21.06 ka、17.99 ka)。本研究关于黑色页岩段的旋回分析将以频谱分析得到的理论天文周期为基准。整个冰沟剖面的旋回周期分析采用的理论周期为1600 ka、405 ka、100 ka(Huang et al., 2010)。

图 4 甘肃张掖南台子下沟组黑色页岩Ⅱ 段早 Aptian期(120— 125 Ma)的岁差、斜率和偏心率的频谱分析结果Fig.4 Spectrum analysis results of precession, obliquity, and eccentricity cycles of the black shale section Ⅱ of Xiagou Formation in Nantaizi, Zhangye of Gansu Province, China during early Aptian(120-125 Ma)

3.3.3 平均频谱拟合差(ASM)

Average Spectral Misfit(ASM)方法能够准确识别和校正地层中的轨道信号。它通过计算每个沉积速率下对应的地层轨道信号与理论天文信号的平均差值(ASM值), 并利用蒙特卡洛模拟测试无轨道信号的零假设可信度, 从而得到最优沉积速率(Meyers and Sageman, 2007)。预处理后的CaCO3数据进行ASM分析, 在整个ASM计算中用到的理论天文周期为123.48 ka、94.36 ka、36.77 ka、22.23 ka、17.99 ka。置信度水平设为0.95%, 设置100个等间距沉积速率, 沉积速率范围设为1~100 cm/ka, 进行10 000次蒙特卡洛模拟。本研究有关ASM的分析全部根据R软件“ astrochron” 包下的asm方法计算。

3.3.4 年代标尺优化法(TimeOpt)

Time-scale Optimization(TimeOpt)综合考虑地层数据周期间层级关系和地层数据计算的轨道信号调制属性, 识别气候主导周期, 从而根据地层厚度计算最优沉积速率(Meyers, 2015; Meyers, 2019)。TimeOpt方法的核心是对天文信号的2个属性进行评估, 即岁差和偏心率的频谱功率与岁差信号的偏心率幅度调制属性的集中度, 这2个属性都是根据线性概率统计模型来评估的。包络线回归模型评估通过带通滤波以及希尔伯格变换提取的岁差振幅包络线与偏心率模型之间的拟合程度, 利用最小二乘估计确定每个沉积速率下所对应的最佳拟合, 每个沉积速率下“ 拟合” 质量是根据偏心率模型和岁差包络线之间的Pearson相关系数来量化的, 从而得到相关性 renvolpoe2; 频谱功率回归模型是评估经过沉积速率校正后的时间校准序列的天文目标周期(岁差与偏心率)频谱功率集中度, 拟合质量是通过时间校准序列与天文模型序列之间的相关性确定, 得到 rspectral2; 最终的拟合度 ropt2结合了振幅包络线与频谱功率评估的信息, 通过蒙特卡洛模拟对无天文信号的零假设进行正式评估, 确定 ropt2的统计显著性。Meyers等(2015)对OPD 926B站点的底栖有孔虫 δ 18O 数据进行TimeOpt方法分析, 在功率谱图中表现出了强烈的斜率信号。当有强烈的斜率信号时, 在频谱功率图中仍然能识别出斜率周期。因此TimeOpt方法不仅能够识别偏心率与岁差, 也能识别出斜率周期。

对整个冰沟剖面(图 2-d)进行旋回周期分析时, 首先提取该剖面的红度值, 将红度值标准化, 并且进行倾角(60° )校正。随后将红度值进行预处理(间隔设置为0.22 m)以及TimeOpt分析, 在整个TimeOpt计算中, 线性概率回归模型使用超长偏心率(1600 ka)、长偏心率周期和短偏心率周期(405 ka 、100 ka), Taner滤波器带宽范围设置为1/250~1/650 cycles/ka(roll-off rate设置为103); 设置100个等间距沉积速率, 沉积速率范围设置为1~40 cm/ka。

对黑色页岩段的CaCO3含量和TOC含量数据进行预处理(CaCO3含量的间隔设置为0.05 m, TOC含量的间隔设置为0.1 m)。预处理后进行TimeOpt分析, 在整个TimeOpt计算中, 线性概率回归模型使用4个偏心率周期(131.61 ka、123.48 ka、99.03 ka、94.36 ka)和3个岁差周期(22.23 ka、21.06 ka、17.99 ka), 岁差包络线利用Taner过滤器提取, 其滤波器带宽范围设置为0.035~0.065 cycles/ka(roll-off rate设置为103); 设置100个等间距沉积速率, CaCO3沉积速率范围为1~100 cm/ka, TOC的沉积速率范围设置为1~20 cm/ka。本研究有关TimeOpt的分析全部根据R的“ astrochron” 包下的TimeOpt和TimeOptSim方法计算。

4 结果
4.1 古气候替代性指标

研究沉积地层气候数据的轨道周期需要选择合理古气候替代性指标。CaCO3含量和TOC含量都是常用的古气候替代性指标。根据前人研究表明TOC含量是指示古湖泊的初级生产力的良好指标(陈发虎等, 2001), 在同样条件下, 湖泊的生产力越高, 则沉积物中的有机质含量高, TOC含量也就越高。但TOC的来源复杂多样, 来源不同的TOC所携带的气候信息不同。TOC/TN值大小能明显区分湖泊中的TOC是湖泊水生植物还是受陆源植物碎屑影响而形成的(钱君龙等, 1997; Meyers and Ishiwatari, 1993; Meyers et al., 1998), 陆源有机质的TOC/TN值为20~30, 湖泊自身的浮游藻类和大型水生植物的TOC/TN值分别为4~10和10~20之间(Meyers and Lallier-vergé s, 1999)。从图 5可知, 黑色页岩Ⅱ 段的TOC/TN值都小于20, 因此黑色页岩Ⅱ 段沉积物的TOC主要来源于湖泊本身。影响湖泊水生植物生长的主要因素是温度与营养盐, 温度升高, 植物光合作用加强, 水生植物大量繁殖, 因此湖泊生产力提高(杨伦庆等, 2009)。在生物生产力高的湖泊中, 大量水生植物生长所需的光合作用会消耗大量CO2, 进而导致CaCO3发生沉淀(宋磊等, 2012)。因此湖泊中生物生产力的提高与CaCO3含量的增加具有良好的相关性(Hodell et al., 1999; 陈敬安等, 2002)。从图 5中可见TOC与CaCO3含量曲线变化相近, 并且在随深度变化趋势以及幅度上具有较好的可对比性。因此TOC与CaCO3的同时富集主要是因为湖泊水生植物的大量繁殖所导致的。与TOC 相比, CaCO3含量的分辨率较高, 对后期的周期分析更有益, 且CaCO3含量测试更便捷、快速和有效。因此综合考虑, 最终将CaCO3含量作为古气候替代性指标。

图 5 甘肃张掖南台子下沟组黑色页岩Ⅱ 段CaCO3含量、TOC含量和TOC/CN值变化
a— 岩性随深度变化曲线; b— CaCO3含量随深度变化曲线; c— TOC含量随深度变化曲线; d— TOC/TN随深度变化曲线Fig.5 Variation of content of CaCO3 and TOC, TOC/CN values of the black shale section Ⅱ of Xiagou Formation in Nantaizi, Zhangye of Gansu Province, China

4.2 黑色页岩研究区浮动年代标尺

地质样品的红色与氧化铁的含量有关, Fe3+代表相对干旱的氧化环境; 灰黑色与有机质和碳酸盐的含量有关, 代表相对还原的湿润环境(刘秀铭等, 2014)。张掖南台子下沟组冰沟剖面的影像(图 6-a)进行红度(a* )提取, 将倾角(60° )进行校正, 可初步建立a* 变化曲线(图 6-a)。将红度值代入TimeOpt分析, 根据包络回归模型以及频谱功率回归模型确定红色的包络线拟合曲线和黑色的频谱功率拟合曲线(图 6-b), 包络线拟合曲线显示了在沉积速率为29.7 cm/ka下最大的 renvolope2值为0.58, 频谱功率拟合曲线显示了在沉积速率为12. 6cm/ka下最大的 rspectral2的值0.16。 renvolope2与 rspectral2的乘积确定了在沉积速率为28. 6cm/ka下最大的 ropt2的值0.114(图 6-c)。使用与数据相同的lag-1自相关系数的AR1替代值(ρ ARI=0.80)进行2000次蒙特卡罗模拟, 得出沉积速率28.6 cm/ka下观测的 ropt2=0.114的ρ -value值为0.005, 这表明可以在99.5%的置信水平上拒绝非天文驱动的零假设(图 6-d)。图 6-e显示在最佳沉积速率28.6 cm/ka时天文目标周期(红色虚线)与地层周期(黑色曲线)的拟合情况, 目标周期与地层周期在超长偏心率、长偏心率处拟合良好, 短偏心率拟合较弱, 表示在地层中超长偏心率与长偏心率周期信号较强, 短偏心率周期信号较弱。TimeOpt重建的超长偏心率模型(黄线)与超长偏心率振幅包络线(蓝线)的拟合结果如图 6-a所示, 2条曲线波动大致相同, 表明二者的拟合情况较好, 并且从图 6-a中能明显看出红度值(红线)的变化波动与长偏心率振幅大小(黑线)一致性良好, 由此看出冰沟剖面的40万年长偏心率周期显著, 并且在地质历史当中40万年碳周期较为常见(汪品先等, 2018)。根据TimeOpt方法计算的浮动年代标尺如图 6-a所示, 整个冰沟剖面的持续沉积时间为2.884 Ma。因此, 将a* 代入TimeOpt方法得到南台子冰沟剖面的沉积速率为28.6 cm/ka, 沉积持续时间为2884 ka, 且具有显著的40万年长偏心率周期。

图 6 甘肃张掖南台子下沟组冰沟剖面红度(a* )的TimeOpt分析
a— 红度(a* )变化曲线(红色), 带通滤波提取的长偏心率信号(黑色), 通过希尔伯格变换得到的超长偏心率幅包络线(蓝色)和TimeOpt重构的偏心率模型(黄色)。b— 相关性r2envolope (红点)和r2spectral(黑色曲线)随沉积速率的变化。c综合振幅包络拟合与频谱功率拟合的r2opt随沉积速率的变化曲线。d— 2000次蒙特卡洛模拟结果, 用以评估 r2opt=0.114的置信水平。e— a* 序列周期图, TimeOpt给定的沉积速率为28.6 cm/ka(黑色曲线为线性频谱, 灰色曲线为对数频谱)。黄色阴影区域代表计算长偏心率振幅包络线的带通滤波器(1/250~1/650 cycles/ka), 垂直的红色虚线表示超长偏心率(E1)、长偏心率(E2)和短偏心率(e)目标周期Fig.6 TimeOpt analysis of a* of the Binggou section of Xiagou Formation in Nantaizi, Zhangye of Gansu Province, China

4.3 下沟组黑色页岩Ⅱ 段主导周期

4.3.1 ASM计算结果

甘肃张掖南台子下沟组黑色页岩Ⅱ 段的CaCO3含量序列ASM分析结果如图 7所示。MTM频谱分析揭示了CaCO3含量序列中有2个超过95%置信水平的天文信号(图 7-a), 将这2个显著天文信号与理论天文信号拟合匹配, 利用蒙特卡洛模拟计算地层轨道信号与天文目标信号的平均差值即ASM值, 当二者一致性最好时即ASM的值最小, 对应的最优沉积速率为17.074 cm/ka(图 7-b)。17.074 cm/ka的沉积速率对应的零假设水平(Ho-SL)为 0.066%, 小于1%的临界值, 通过了置信度检验, 表明 CaCO3 含量序列中有99.934%的可能性记录了米兰科维奇信号(图 7-b)。图 7-c显示了经最优沉积速率17.074 cm/ka校正后的CaCO3含量序列的周期信号(黑色虚线)与理论周期信号(红色实线)拟合情况, 结果显示CaCO3含量序列中的周期信号与理论偏心率和岁差拟合情况良好, 说明地层中的偏心率和岁差信号显著。同时, 根据最优沉积速率17.074 cm/ka计算29.85 m长的黑色页岩段的沉积持续时间为175 ka。因此, 根据ASM分析得知黑色页岩Ⅱ 段的主导周期为偏心率与岁差周期。

图 7 甘肃张掖南台子下沟组黑色页岩Ⅱ 段CaCO3含量序列ASM分析
a— MTM频谱分析图。黑色曲线代表频谱功率随频率变化, 蓝色虚线表示95%的置信度水平。b— 零假设置信水平。红色虚线指示最优沉积速率17.074 cm/ka, 黑色虚线指示的为设置信度1%的临界值。c— 最优沉积速率校正后的显著地层频率(黑色虚线)与文目标频率(红色实线)的拟合Fig.7 ASM analysis of CaCO3 centent sequence in the black shale section Ⅱ of Xiagou Formation in Nantaizi, Zhangye of Gansu Province, China

4.3.2 TimeOpt计算结果

TimeOpt方法对CaCO3含量数据进行分析的结果如图 8。包络线拟合曲线显示了在沉积速率为54.62 cm/ka下最大的 renvolope2值为0.99, 频谱功率拟合曲线显示了在沉积速率为18.74 cm/ka下最大的 rspectral2的值0.24(图 8-a)。 renvolope2与 rspectral2的乘积确定了最大的 ropt2的值0.228, 对应的沉积速率为18.74 cm/ka(图 8-b)。使用与数据相同的lag-1自相关系数的AR1替代值(ρ ARI=0.65)进行2000次蒙特卡罗模拟, 得出在沉积速率18.74 cm/ka下观测的 ropt2值为 0.228时的ρ -value值为0.005, 表明有99.5%的可能性拒绝非天文驱动的零假设(图 8-c)。图 8-d显示在最佳沉积速率18.74 cm/ka时天文目标周期(红色虚线)与地层周期(黑色曲线)的对比情况, 目标周期与地层周期在岁差与偏心率处一致性良好, 而斜率周期的频谱功率较弱, 表示在地层中岁差与短偏心率周期信号较强, 斜率周期信号弱。岁差振幅包络线(图 8-e的红色曲线)与TimeOpt重建的偏心率模型(图 8-e的黑色曲线)的对比结果(图 8-e)以及二者的交叉图中的2条曲线(图 8-f)近乎重合, 表明二者的一致性良好。时间域CaCO3含量曲线(蓝色)与所提取的岁差滤波(黑色)清晰地反映CaCO3与岁差的对应关系较好(图 8-g), 岁差的低值对应着CaCO3含量的低值, 由此看出CaCO3含量的岁差周期明显。不仅如此, 这也反映出了偏心率对岁差的调制作用, 当偏心率(红色曲线)幅度小时, 岁差的变化幅度也小。因此, 将CaCO3含量数据代入TimeOpt方法得到南台子下沟组黑色页岩的主导周期都为岁差与短偏心率周期。

图 8 甘肃张掖南台子下沟组黑色页岩Ⅱ 段CaCO3含量数据的TimeOpt分析
a— 相关性r2envolope (红点)和r2spectral (黑色曲线)随沉积速率的变化; b— 综合振幅包络拟合与频谱功率拟合的r2opt随沉积速率的变化曲线; c— 2000次蒙特卡洛模拟结果, 用以评估r2opt=0.228的置信水平; d— CaCO3含量序列周期图, TimeOpt给定的沉积速率为18.74cm/ka(频谱); 黄色阴影区域代表计算岁差振幅包络线的带通滤波器(0.035~0.065 cycles/ka); e— CaCO3含量序列岁差振幅包络线(红色曲线)与在偏心率模型(黑色曲线)的对比图; f— CaCO3含量序列岁差振幅包络线与TimeOpt重构的偏心率模型的交叉图; g— 时间域 CaCO3含量序列(蓝色曲线), 带通滤波提取的岁差信号(黑色)和通过希尔伯格变换得到的数据振幅包络线(红色); h— TOC含量序列周期图, TimeOpt给定的沉积速率为11.60 cm/ka(黑色曲线为线性频谱, 灰色曲线为对数频谱); 黄色阴影区域代表计算岁差振幅包络线的带通滤波器(0.035~0.065 cycles/ka), 垂直的红色虚线表示短偏心率(e)和岁差(p)目标周期Fig.8 TimeOpt analysis of CaCO3 content of the black shale section Ⅱ of Xiagou Formation in Nantaizi, Zhangye of Gansu Province, China

冰沟剖面黑色页岩持续沉积时间为1.575 Ma(图 6-a的灰色阴影), 每个黑色页岩段与其邻近的黄色砂岩段的持续沉积时间为400 ka, 其中旋回Ⅱ (图 6-a)总长为100 m, 平均沉积速率25 cm/ka。旋回Ⅱ 黑色页岩部分总厚约30 m, 沉积持续时间为160 ka(图 8-e), 沉积速率为18.74 cm/ka, 沉积速率比整个冰沟剖面的平均沉积速率(28.0 cm/ka)低。砂岩段总厚70 m, 沉积持续时间为240 ka, 沉积速率为29.2 cm/ka, 砂岩段沉积速率比整个冰沟剖面的平均沉积速率高。

另外, 将TOC数据用TimeOpt方法得到的周期频谱图如图 8-h所示, 地层数据的偏心率周期与天文目标周期一致性良好, 岁差拟合度较差以及斜率的频谱功率较弱, 表明主导周期为偏心率周期, 而斜率与岁差周期信号微弱。这也能反映CaCO3含量指标更具优势, CaCO3含量周期图(图 8-d)显示有显著的偏心率与岁差周期, 而 TOC含量周期图(图 8-h)只显示有明显的偏心率周期, TOC含量没有显著的岁差周期可能与TOC含量分辨率不高、反应不灵敏有关, 因而只显示出较大的偏心率周期。

4.3.3 频谱分析

南台子下沟组黑色页岩Ⅱ 段的天文浮动年代标尺是根据TimeOpt经沉积速率调谐后所建立的, 将深度域CaCO3含量序列转化为时间域序列, 黑色页岩Ⅱ 段的沉积持续时间约160 ka(图 9-a)。CaCO3时间域序列的MTM频谱分析(图 9-b)显示在99.77 ka、25.75 ka和17.74 ka处有显著谱峰, 置信度水平分别在90%以上。时间域序列的EHA分析(图 9-b)结果发现, 99.77 ka、25.75 ka和17.74 ka信号在整个EHA中呈现连续的功率强的红色信号, 表明短偏心率和岁差周期显著。综上, 张掖南台子下沟组黑色页岩剖面具有显著的短偏心率周期与岁差周期。

图 9 甘肃张掖南台子下沟组黑色页岩Ⅱ 段调谐后年代数据的MTM和EHA结果
a— CaCO3含量随时间的变化曲线。b— MTM和EHA分析结果图:MTM图中红色虚线分别解释为90%的置信度水平, 黑色虚线上标注的数字分别代表 99.77 ka、25.75 ka和17.74 ka; EHA图中带有“ e” 和“ p” 标签的白色虚线分别代表解释的短偏心率和岁差周期Fig.9 Result of MTM and EHA with tuned age data of the black shale section Ⅱ of Xiagou Formation in Nantaizi, Zhangye of Gansu Province, China

5 讨论:低纬岁差驱动

白垩纪中期欧洲海相地层(古纬度约20° ~30° N)以斜率周期为主导。但是相似的纬度, 本研究区(古纬度约23° N)以岁差和偏心率为主导周期。主导周期与欧洲海相地层不同。

第四纪以来低纬地区的洞穴石笋氧同位素数据、孢粉记录以及赤道大西洋、南海等记录的气候变化主要以~20 ka岁差周期为主, 低纬区域显示出了显著的~20 ka岁差周期。理论上, 岁差主要控制着中低纬度的太阳辐射量的变化, 偏心率影响较为微弱(Cheng et al., 2021), 例如东亚季风区的石笋记录都呈现出岁差信号显著, 而偏心率信号相对微弱(Wang et al., 2016; 程海等, 2020)。但在地质记录中发现岁差与偏心率周期同样显著的案例, 并将偏心率显著的原因解释为偏心率对岁差的调谐作用, 如地中海所形成的浅色的石灰岩序列是因为受到偏心率调谐后的弱季风降水无法阻断地中海循环, 从而形成浅色石灰岩序列; 深色腐泥层序列的形成是在偏心率调谐下的强季风降水导致更多的淡水资源流入地中海, 从而阻断地中海循环, 在缺氧条件下形成腐泥层序列, 该腐泥层序列的岁差与偏心率周期都非常显著(Ruddiman, 2001)。研究区不仅显示出显著的岁差信号, 同时也显示了强烈的偏心率信号(图 6-d), 并将这一现象解释为偏心率对岁差的调谐作用。因此南台子下沟组黑色页岩Ⅱ 段显示的岁差和偏心率周期与低纬太阳辐射相关。

青藏高原约55~50 Ma开始隆升(Searle et al., 1987; Jiang et al., 2001), Jiang等(2001)根据沙漠的纬向分布特点以及沙漠盛行风的相关记录认为在青藏高原隆起前东亚是由行星风系所控制。白垩纪亚欧非大陆的位置与形态(图 1)接近于第四纪以来的亚欧非大陆, 相当于现代欧洲— 亚洲— 非洲大陆, 西侧与美洲大陆相邻, 因此这样的地形条件下中低纬度区域可能存在大规模的季风以及由季风带来的大范围降水(程海等, 2020)。四川盆地与思茅盆地(图 1)的纬度与研究区相近, 属于干旱的沙漠气候(江新胜和李玉文, 1996; Hasegawa et al., 2012)。本研究区黑色页岩年代属于青藏高原未隆起的Aptian早期(124 Ma), 位于23° N, 以河湖相沉积为主, 存在大量黑色页岩, 黑色页岩的存在代表着湿润环境(戴霜等, 2012)。研究区紧靠副热带高压控制下的沙漠干旱区南缘(图 1), 当时的张掖地区应该与今天的非洲撒哈拉沙漠南缘较为接近, 其气候型为热带干湿季气候, 与非洲季风的气候过程相当, 可能受到热带辐合带南北移动和海陆热力差异的影响。但这还需要古气候模拟等更多证据支撑。

6 结论

本研究对白垩纪中期低纬中国甘肃张掖南台子下沟组黑色页岩段进行分析, 探讨了白垩纪低纬度与轨道周期间的关系, 得出主要结论如下:

1)以CaCO3含量作为高分辨率气候替代性指标, 用ASM、TimeOpt和频谱分析方法计算得到南台子黑色页岩的主导周期都为岁差和偏心率。

2)Aptian期欧洲低纬海相黑色页岩揭示的气候主导周期的多为斜率周期, 而处于相同纬度的中国陆相南台子黑色页岩的主导周期为岁差周期。

(责任编辑 李新坡; 英文审校 陈吉涛)

参考文献
[1] 陈发虎, 朱艳, 李吉均, 施祺, 靳立亚, B. Wunemann. 2001. 民勤盆地湖泊沉积记录的全新世千百年尺度夏季风快速变化. 科学通报, 46(17): 1414-1419.
[Chen F H, Zhu Y, Li J J, Shi Q, Jin L Y, B. Wunemann. 2001. Rapid changes in Cenozoic centennial scale summer monsoon recorded by sediments in Minqin Basin. Chinese Science Bulletin, 46(17): 1414-1419] [文内引用:1]
[2] 陈敬安, 万国江, 汪福顺, 黄荣贵, 张峰, D. D. Zhang, R. Schmidt. 2002. 湖泊现代沉积物碳环境记录研究. 中国科学(D辑: 地球科学), 32(1): 73-80.
[Chen J A, Wan G J, Wang F S, Huang R G, Zhang F, D. D. Zhang, R. Schmidt. . 2002. Researches onenvironmental records of carbon from lake sedilments. Science in China(Series D: Earth Sciences), 32(1): 73-80] [文内引用:1]
[3] 程海, 李瀚瑛, 张旭, 张海伟, 易亮, 蔡演军, 胡永云, 石正国, 彭友兵, 赵景耀, Gayatri Kathayat, Ashish Sinha. 2020. 欧—亚—非大陆季风: 超级大陆与超级季风的雏形. 第四纪研究, 40(6): 1381-1396.
[Cheng H, Li H Y, Zhang X, Zhang H W, Yi L, Cai Y J, Hu Y Y, Shi Z G, Peng Y B, Zhao J Y, G Kathayat, A Sinha. 2020. European-Asian-African continent: an early form of supercontinent and supermonsoon. Quaternary Sciences, 40(6): 1381-1396] [文内引用:2]
[4] 戴霜, 刘学, 赵杰, 张明震, 刘俊伟, 孔立, 朱强, 黄永波. 2012. 陆地沉积物对大洋缺氧事件的响应: 六盘山群黑色页岩地球化学特征及其意义. 地学前缘, 19(4): 255-259.
[Dai S, Liu X, Zhao M Z, Liu J W, Kong L, Zhu Q, Huang Y B. 2012. The OAEs record in the terrestrial sediments: the geochemistry of blackshales in the Liupanshan Group and its paleoclimatic implications. Earth Science Frontiers, 19(4): 255-259] [文内引用:1]
[5] 江新胜, 李玉文. 1996. 中国中东部白垩纪沙漠的时空分布及其气候意义. 岩相古地理, 16(2): 42-51.
[Jiang X S, Li Y W. 1996. Spatio-temporal distribution of the cretaceous deserts in central and eastern China and its climatic significance. Sedimentary Facies and Palaeogeography, 16(2): 42-51] [文内引用:1]
[6] 刘秀铭, 吕镔, 毛学刚, 温昌辉, 俞鸣同, 郭雪莲, 陈家胜, 王涛. 2014. 风积地层中铁矿物随环境变化及其启示. 第四纪研究, 34(3): 443-457.
[Liu X M, Lü B, Mao X G, Wen C H, Yu M T, Guo X L, Chen J S, Wan T. 2014. Iron minerals of aeoliand eposits vary with environment and its significances. Quaternary Sciences, 34(3): 443-457] [文内引用:1]
[7] 潘忠习, 江新胜, 傅清平. 1999. 四川盆地白垩纪沙漠沉积磁组构特征及其古风向意义. 岩相古地理, 19(1): 14-21.
[Pan Z X, Jiang X S, Fu Q P. 1999. The variations of palaeowind direction of the cretaceous desert in the Sichuan Basin and their significance. Sedimentary Facies and Palaeogeography, 19(1): 14-21] [文内引用:1]
[8] 钱君龙, 王苏民, 薛滨, 陈如松, 柯善哲. 1997. 湖泊沉积研究中一种定量估算陆源有机碳的方法. 科学通报, 47(15): 1655-1658.
[Qian J L, Wang S M, Xue B, Chen R S, Ke S Z. 1997. A method of estimationof terrestrial organic carbon in the study of lake. Chinese Science Bulletin, 47(15): 1655-1658] [文内引用:1]
[9] 宋磊, 强明瑞, 郎丽丽, 刘星星, 王亲, 李明治. 2012. 16 ka BP共和盆地更尕海湖泊生产力演化历史. 科学通报, 57(19): 1763-1774.
[Song L, Qiang M R, Lang L L, Liu X X, Wang Q, Li M Z. 2012. Changes inpalaeoproductivity of Genggahai Lake over the past 16 ka in the Gonghe basin, northeastern Qinghai-Tibetan Plateau. Chinese Science Bulletin, 57(19): 1763-1774] [文内引用:1]
[10] 汪品先, 田军, 黄恩清, 马文涛. 2018. 地球系统与演变. 北京: 科学出版社, 490-491.
[Wang P X, Tian J, Huang E Q, Ma W T. 2018. Earth Systems and Evolution. Beijing: Science Press, 490-491] [文内引用:2]
[11] 杨伦庆, 张虎才, 类延斌, 雷国良, 张文翔, 常凤琴, 庞有智. 2009. 内蒙古西部额济纳古湖小狐山剖面有机质与碳酸盐组成及其古环境意义. 第四纪研究, 29(2): 256-267.
[Yang L Q, Zhang H C, Lei Y B, Lei G L, Zhang W X, Chang F Q, Pang Y Z. 2009. Compositions of the organic matter and carbonate in Xiaohushan section from Ejina Basin, western Inner Mongolia, China and their paleoenvironmental significance. Quaternary Sciences, 29(2): 256-267] [文内引用:1]
[12] 杨雨. 1997. 甘肃省岩石地层. 武汉: 中国地质大学出版社, 1-220.
[Yang Y. 1997. Lithostratigraphic of Gansu Province. Wuhan: China University of Geosciences Press, 1-220] [文内引用:2]
[13] 张荷生, 崔振卿. 2007. 甘肃省张掖丹霞与彩色丘陵地貌的形成与景观特征. 中国沙漠, 27(6): 942-945.
[Zhang H S, Cui Z Q. 2007. Land scape character and forming of Danxia Land form and Color Hill in Zhangye of Gansu Province. Journal of Desert Research, 27(6): 942-945] [文内引用:1]
[14] 张兰生, 方修琦. 2012. 中国古地理: 中国自然环境的形成. 北京: 科学出版社, 1-531.
[Zhang L S, Fang X Q. 2012. Palaeogeography of China: The Formation of Natural Environment in China. Beijing: Science Press, 1-531] [文内引用:1]
[15] Beaufort L, Lancelot Y, Camberlin P, Cayre O, Vincent E, Bassinot F, Labeyrie L. 1997. Insolation cycles as a major control of equatorial Indian Ocean primary production. Science, 278(5342): 1451-1454. [文内引用:1]
[16] Beil S, Kuhnt W, Holbourn A, Scholz F, Oxmann J, Wallmann K, Lorenzen J, Aquit M, Chellai E H. 2020. Cretaceous oceanic anoxic events prolonged by phosphorus cycle feedbacks. Climate of the Past, 16(2): 757-782. [文内引用:1]
[17] Berner R A, Kothavala Z. 2001. GEOCARB Ⅲ: a revised model of atmospheric CO2 over Phanerozoic time. American Journal of Science, 301(2): 182-204. [文内引用:1]
[18] Bornemann A, Norris R D, Friedrich O, Beckmann B, Schouten S, Sinninghe Damsté J S, Vogel J, Hofmann P, Wagner T. 2008. Isotopic evidence for glaciation during the cretaceous supergreenhouse. Science, 319(5860): 189-192. [文内引用:1]
[19] Cai Y, Fung I Y, Edwards R L, An Z, Cheng H, Lee J E, Tan L, Shen C C, Wang X, Day J A, Zhou W, Kelly M J, Chiang J C H. 2015. Variability of stalagmite-inferred Indian monsoon precipitation over the past 252 000 y. Proceedings of the National Academy of Sciences, 10(112): 2954. [文内引用:1]
[20] Cheng H, Edwards R L, Broecker W S, Denton G H, Kong X, Wang Y, Zhang R, Wang X. 2009. Ice age terminations. Science, 326(5950): 248-252. [文内引用:1]
[21] Cheng H, Edwards R L, Sinha A, Spötl C, Yi L, Chen S, Kelly M, Kathayat G, Wang X, Li X, Kong X, Wang Y, Ning Y, Zhang H. 2016. The Asian monsoon over the past 640 000 years and ice age terminations. Nature, 534(7609): 640-646. [文内引用:1]
[22] Cheng H, Zhang H, Cai Y, Shi Z, Yi L, Deng C, Hao Q, Peng Y, Sinha A, Li H. 2021. Orbital-scale Asian summer monsoon variations: paradox and exploration. Science China Earth Sciences, 64(4): 529-544. [文内引用:1]
[23] Cheng H, Zhang P, Spötl C, Edwards R, Cai Y, Zhang D, Sang W, Tan M, An Z. 2012. The climatic cyclicity in semiarid-arid central Asia over the past 500 000 years. Geophysical Research Letters, 39(1): L01705. [文内引用:1]
[24] Forster A, Schouten S, Moriya K, Wilson P A, Sinninghe Damsté J S. 2007. Tropical warming and intermittent cooling during the Cenomanian/Turonian oceanic anoxic event 2: sea surface temperature records from the equatorial Atlantic. Paleoceanography, 22(1): PA1219. [文内引用:1]
[25] Friedrich O, Norris R D, Erbacher J. 2012. Evolution of middle to late cretaceous oceans: a 55 My record of earth’s temperature and carbon cycle. Geology, 40(2): 107-110. [文内引用:1]
[26] Frost G M, Coe R S, Meng Z, Peng Z, Chen Y, Courtillot V, Peltzer G, Tapponnier P, Avouac J-P. 1995. Preliminary early cretaceous paleomagnetic results from the Gansu Corridor, China. Earth and Planetary Science Letters, 129(1-4): 217-232. [文内引用:1]
[27] Hasegawa H, Tada R, Jiang X, Suganuma Y, Imsamut S, Charusiri P, Ichinnorov N, Khand Y. 2012. Drastic shrinking of the Hadley circulation during the mid-cretaceous supergreenhouse. Climate of the Past, 8(4): 1323-1337. [文内引用:1]
[28] Hinnov L A. 2013. Cyclostratigraphy and its revolutionizing applications in the earth and planetary sciences. GSA Bulletin, 125(11-12): 1703-1734. [文内引用:1]
[29] Hinnov L A. 2000. New perspectives on orbitally forced stratigraphy. Annual Review of Earth and Planetary Sciences, 28(1): 419-475. [文内引用:2]
[30] Hodell D A, Brenner M, Kanfoush S L, Curtis J H, Stoner J S, Xueliang S, Yuan W, Whitmore T J. 1999. Paleoclimate of southwestern China for the past 50 000 yr inferred from lake sediment records. Quaternary Research, 52(3): 369-380. [文内引用:1]
[31] Huang C, Hinnov L, Fischer A G, Grippo A, Herbert T. 2010. Astronomical tuning of the aptian stage from Italian reference sections. Geology, 238(10): 899-903. [文内引用:2]
[32] Huber B T, Norris R D, MacLeod K G. 2002. Deep-sea paleotemperature record of extreme warmth during the cretaceous. Geology, 30(2): 123-126. [文内引用:1]
[33] Jenkyns H C, Forster A, Schouten S, Damsté J. 2004. High temperatures in the late cretaceous Arctic Ocean. Nature, 432(7019): 888-892. [文内引用:1]
[34] Jiang X S, Pan Z, Fu Q. 2001. Primary study on pattern of general circulation of atmosphere before uplift of the Tibetan Plateau in eastern Asia. Science in China Series D: Earth Sciences, 44(8): 680-688. [文内引用:2]
[35] Kathayat G, Cheng H, Sinha A, Spötl C, Edwards R L, Zhang H, Li X, Yi L, Ning Y, Cai Y, Lui W L, Breitenbach S F M. 2016. Indian monsoon variability on millennial-orbital timescales. Scientific Reports, 6(1): 1-7. [文内引用:1]
[36] Laskar J, Robutel P, Joutel F, Gastineau M, Correia A C M, Levrard B. 2004. A long-term numerical solution for the insolation quantities of the Earth. Astron Astrophys, 428(1): 261-285. [文内引用:1]
[37] Li X, Xu W, Liu W, Zhou Y, Wang Y, Sun Y, Liu L. 2013. Climatic and environmental indications of carbon and oxygen isotopes from the Lower Cretaceous calcrete and lacustrine carbonates in Southeast and Northwest China. Palaeogeography, Palaeoclimatology, Palaeoecology, 385(1): 171-189. [文内引用:2]
[38] Li Y X, Bralower T J, Montañez I P, Osleger D A, Arthur M A, Bice D M, Herbert T D, Erba E, Silva I P. 2008. Toward an orbital chronology for the early aptian oceanic anoxic event(OAE1a, ~120 Ma). Earth and Planetary Science Letters, 271(1-4): 88-100. [文内引用:3]
[39] Littler K, Robinson S A, Bown P R, Nederbragt A J, Pancost R D. 2011. High sea-surface temperatures during the early cretaceous epoch. Nature Geoscience, 4(3): 169-172. [文内引用:1]
[40] Liu Z, Liu X, Huang S. 2017. Cyclostratigraphic analysis of magnetic records for orbital chronology of the lower cretaceous Xiagou Formation in Linze, northwestern China. Palaeogeography, Palaeoclimatology, Palaeoecology, 481: 44-56. [文内引用:1]
[41] Meyers P A, Ishiwatari R. 1993. Lacustrine organic geochemistry-an overview of indicators of organic matter sources and diagenesis in lake sediments. Organic Geochemistry, 20(7): 867-900. [文内引用:1]
[42] Meyers P A, Lallier-vergés E. 1999. Lacustrine sedimentary organic matter records of late quaternary paleoclimates. Journal of Paleolimnology, 21(3): 345-372. [文内引用:1]
[43] Meyers P A, Tenzer G E, Lebo M E, Reuter J E. 1998. Sedimentary record of sources and accumulation of organic matter in Pyramid Lake, Nevada, over the past 1000 years. Limnology and Oceanography, 43(1): 160-169. [文内引用:1]
[44] Meyers S R. 2019. Cyclostratigraphy and the problem of astrochronologic testing. Earth-Science Reviews, 190: 190-223. [文内引用:1]
[45] Meyers S R. 2015. The evaluation of eccentricity-related amplitude modulation and bundling in paleoclimate data: an inverse approach for astrochronologic testing and time scale optimization. Paleoceanography, 30: 1625-1640. [文内引用:2]
[46] Meyers S R, Sageman B B. 2007. Quantification of deep-time orbital forcing by average spectral misfit. American Journal of Science, 307(5): 773-792. [文内引用:1]
[47] Meyers S R, Sageman B B, Hinnov L A. 2001. Integrated quantitative stratigraphy of the Cenomanian-Turonian bridge creek limestone member using evolutive harmonic analysis and stratigraphic modeling. Journal of Sedimentary Research, 71(4): 628-644. [文内引用:2]
[48] Milankovitch M. 1941. Kanon der Erdbestrahlung und seine Anwendung auf das Eiszeitproblem. Belgrade: Royal Serbian Academy Special Publications, 1-633. [文内引用:1]
[49] Ruddiman W F. 2001. Earth’s Climate: Past and Future. New York: W. H. Freeman and Company, 465. [文内引用:3]
[50] Schlanger S O, Jenkyns H. 1976. Cretaceous oceanic anoxic events: causes and consequences. Geologie en Mijnbouw, 55(3): 179-184. [文内引用:1]
[51] Searle M, Windley B, Coward M, Cooper D, Rex A, Rex D, Tingdong L, Xuchang X, Jan M, Thakur V. 1987. The closing of Tethys and the tectonics of the Himalaya. Geological Society of America Bulletin, 98(6): 678-701. [文内引用:1]
[52] Short D A, Mengel J G, Crowley T J, Hyde W T, North G R. 1991. Filtering of milankovitch cycles by earth’s geography. Quaternary Research, 35(2): 157-173. [文内引用:1]
[53] Suarez M B, Ludvigson G A, González L A, Al-Suwaidi A H, You H L. 2013. Stable isotope chemostratigraphy in lacustrine strata of the Xiagou Formation, Gansu Province, NW China. Geological Society London Special Publications, 382(1): 143-155. [文内引用:3]
[54] Thomson D J. 1982. Spectrum estimation and harmonic analysis. Proceedings of the IEEE, 70(9): 1055-1096. [文内引用:1]
[55] Tiedemann R, Sarnthein M, Shackleton N J. 1994. Astronomic timescale for the Pliocene Atlantic δ18O and dust flux records of Ocean Drilling Program Site 659. Paleoceanography, 9(4): 619-638. [文内引用:1]
[56] Tiraboschi D, Erba E, Jenkyns H C. 2009. Origin of rhythmic Albian black shales(Piobbico core, central Italy): calcareous nannofossil quantitative and statistical analyses and paleo-ceanographic reconstructions. Paleoceanography, 24(2): PA2222. [文内引用:1]
[57] Wang P, Li Q, Tian J, He J, Jian Z, Ma W, Dang H. 2016. Monsoon influence on planktic δ18O records from the South China Sea. Quaternary Science Reviews, 142: 26-39. [文内引用:1]
[58] Wang S Y, Houyuan L, Liu J Q, Negendank J F W. 2007. The early Holocene optimum inferred from a high-resolution pollen record of Huguangyan Maar Lake in southern China. Chinese Science Bulletin, 52(20): 2829-2836. [文内引用:1]
[59] Wang Y, Cheng H, Edwards R L, He Y, Kong X, An Z, Wu J, Kelly M J, Dykoski C A, Li X. 2005. The holocene asian monsoon: links to solar changes and north atlantic climate. Science, 308(5723): 854-857. [文内引用:1]
[60] Wilson P A, Norris R D, Cooper M J. 2002. Testing the cretaceous greenhouse hypothesis using glassy foraminiferal calcite from the core of the turonian tropics on demerara rise. Geology, 30(7): 607-610. [文内引用:1]
[61] Xi D, Wan X, Li G, Gang L. 2019. Cretaceous integrative stratigraphy and timescale of China. Science China Earth Science, 62(1): 256-286. [文内引用:2]
[62] Zhang X, Chen K, Hu D, Sha J. 2016. Mid-Cretaceous carbon cycle perturbations and oceanic anoxic events recorded in southern Tibet. Scientific Reports, 6(1): 39643. [文内引用:1]