Early Pleistocene records of carbonate burial and terrestrial input in the Timor Sea and their paleoclimatic implications
-
摘要:
深海碳酸盐埋藏作为地球表层碳库的重要碳汇，在地质历史时期与大气二氧化碳浓度变化息息相关。古记录重建显示，大气二氧化碳浓度并不总是与全球平均气温具有较好的对应关系，其他因素可能在过去全球变化中起到了重要的作用。本文以澳大利亚西北岸外帝汶海IODP U1482站钻孔2～1.07 Ma沉积物为研究材料，测试其碳酸盐和主微量元素含量，探讨深海碳酸盐埋藏的影响因素。结果显示，碳酸盐含量与指示河流输入的钾元素含量、指示生产力的铀元素含量和底栖有孔虫δ13C以及指示风尘输入的log(Zr/Rb)等记录的长期变化趋势均可分为2～1.63 Ma、1.63～1.31 Ma和1.31～1.07 Ma等3个阶段，可能与沃克环流和哈德莱环流调控的印尼—澳洲地区的干湿条件有关。1.31 Ma之后钾元素含量的持续降低和碳酸盐含量的增加可能揭示哈德莱环流的加强导致了澳洲西北内陆的干旱化趋势。早更新世U1482站碳酸盐含量在轨道时间尺度上主要受以河流输入为主的陆源沉积物稀释的影响；频谱分析显示其具有显著的19 ka和约29 ka变化周期，可能指示该研究时段内除了岁差周期外，以斜率周期为主导的冰期-间冰期旋回对区域降水和陆源沉积物输入的调控作用。
Abstract:Deep-sea carbonate burial is an important carbon sink of the Earth’s surface carbon reservoir. Its variations are closely related to that of atmospheric carbon dioxide. Paleo-environment reconstructions show that changes in atmospheric CO2 concentration do not always correlate well with past temperature changes, implying that other factors may contribute to the past global change. We examined deep-sea sediments spanning the interval of 2～1.07 Ma retrieved from IODP Site U1482 located offshore northwestern Australia in the Timor Sea. Carbonate and elemental contents were analyzed to investigate the factors that influenced carbonate burial. The early Pleistocene records at Site U1482 indicate that all of our records including carbonate content, potassium content (a proxy of regional precipitation), log(Zr/Rb) (a proxy of aeolian dust input), benthic foraminiferal δ13C, and uranium content (proxies of paleoproductivity) show long-term changes punctuated at 1.63 Ma and 1.31 Ma, possibly related to the precipitation pattern over Indonesia-Australia modulated by Walker and Hadley circulations. The continuous decrease in carbonate content and increase in potassium content after 1.31 Ma likely implicate that the intensified aridification in the northwestern Australian hinterland was caused by strengthening of the Hadley circulation. Dilution by terrestrial sediments, mainly of riverine origin, is the predominant factor modifying carbonate content on orbital timescales. Spectral analysis shows that the carbonate record was dominated by 19 ka and 29 ka cycles, likely suggesting the effect of the obliquity-paced glacial/interglacial cyclicity on the regional precipitation and terrestrial input in addition to the precessional control.
-
Keywords:
- carbonate content /
- regional precipitation /
- early Pleistocene /
- Timor Sea
-
南华北盆地位于华北板块东南部，是在华北古生代克拉通基础上发育起来的的中、新生代叠合盆地。盆地晚古生代沉积背景与华北古陆一致，海水自东北向西南侵进，石炭—二叠纪发育一套以海陆交互相为主的具陆源碎屑含煤建造特征的内陆盆地沉积[1]。作为我国东部重要的油气远景区之一，南华北盆地石炭—二叠系沉积地层记录完整，煤系地层发育，是该区重要的页岩气生储有利层系。但是，海陆交互相地层砂岩和泥岩频繁交互，垂向上具有较强的非均质性，泥页岩的发育受到沉积相带的控制，泥页岩中有机质的富集和保存也受到沉积环境的古气候、古水深、古盐度和氧化还原特性等多方面的影响。因此，恢复煤系地层的古沉积环境在研究本区有机质富集规律、确定页岩气有利层段等方面有重要意义。
稀土元素化学性质稳定，在风化剥蚀和沉积成岩过程中不易受变质作用影响，可以被定量地保存下来，因此已被广泛运用于确定沉积源区、研究古气候特征和沉积水体条件等诸多方面[2-5]。本研究对南华北盆地太康隆起上古生界山西组和太原组岩心进行系统取样，通过泥岩样品稀土元素在垂向上的变化特征，探讨本区石炭—二叠系含煤地层沉积过程中的古水体条件、古气候条件及其对有机质富集的影响。
1. 区域地质背景
南华北盆地西起豫西隆起，东临郯庐断裂带，北接渤海湾盆地，南至秦岭大别造山带，大地构造上分为豫西隆起区、开封坳陷、太康隆起、周口坳陷、长山隆起、信阳—合肥坳陷、淮北隆起等7 个一级构造单元[6]。太康隆起是南华北盆地的次级构造单元，东为徐蚌隆起，西接豫西隆起，南与周口坳陷相邻，北为开封坳陷，区内发育新近系、古近系、三叠系、二叠系、石炭系、奥陶系、寒武系等地层,经历了嵩阳、五台、晋宁、加里东、海西、印支、燕山、喜马拉雅运动等多次构造运动，其内部构造格局及地层保存等多受东秦岭－大别造山带及郯庐断裂系影响[7]。晚石炭世，华北板块与西伯利亚I板块发生对接碰撞，古地势演化为北高南低，海水由东南方向入侵，在华北地区形成了广阔的陆表海环境，沉积了一套准碳酸盐台地相和三角洲—潟湖潮坪相的暗色砂泥岩、灰岩和煤层，是华北地区主要的成煤期之一[8]。在此背景下，本区形成了连片发育的泥页岩和煤层，具备一定的烃源岩基础。
2. 样品采集及测试方法
本次研究所选剖面位于南华北盆地太康隆起西部开封市尉氏县，所有样品采自上古生界参数井尉参1井，具体井位见图1。本文以上古生界煤系地层—山西组和太原组为主要研究对象，另从上下整合接触地层本溪组和下石盒子组中选择性采样，作为垂向连续变化趋势参照对比，自本溪组顶部至下石盒子组共取样28件，其中本溪组3件、太原组6件、山西组13件、下石盒子组6件。样品均为泥岩或粉砂质泥岩，未受风化作用影响，具体采样位置见图2。样品分析由中石化胜利油田研究院石油地质测试中心完成，稀土元素和微量元素（Ni和Cr）采用NexION300D 等离子体质谱仪进行测定，岩石中有机碳分析采用有机碳分析仪（CS-600）依据GB/T 19145—2003测定。分析结果见表1。
![]()
图 1 研究区位置及早二叠世岩相古地理图 （据曾秋楠，2019[9]修改）a. 太康隆起区域构造；b. 太康隆起山西组沉积相；c. 尉参1井上古生界岩性及采样位置。Figure 1. Palaeogeographic map and location map of Early Permian of the study area表 1 太康隆起上古生界泥岩稀土元素质量分数表Table 1. REE content of the Upper Paleozoic samples on Taikang Uplift层位 样品
编号La Ce Pr Nd Sm Eu Gd Tb Dy Ho Er Tm Yb Lu Y ∑REE LREE HREE L/H TOC 下石盒子组 X-1 47.4 106 12.8 51.5 7.82 1.2 5.4 0.822 4.56 0.974 3.15 0.599 3.99 0.559 24.6 271.374 226.72 20.054 11.31 0.12 X-2 10.2 23.9 3.11 12.5 2.84 1.09 2.76 0.597 3.8 0.746 2.26 0.434 2.81 0.428 19.7 87.175 53.64 13.835 3.88 X-3 54.4 104 12 47.1 8.93 1.92 7.52 1.18 6.13 1.23 3.49 0.63 4.09 0.606 35.4 288.626 228.35 24.876 9.18 0.61 X-4 75.7 147 17.4 68.8 11.5 1.85 8.16 1.2 6.23 1.27 4 0.741 4.78 0.69 36.8 386.121 322.25 27.071 11.90 0.55 X-5 27.1 61.5 7.58 27.8 3.42 0.472 3.06 0.561 3.68 0.848 2.84 0.576 3.78 0.548 20.9 164.665 127.872 15.893 8.05 0.13 X-6 53.1 98.3 11.3 40.6 6.36 1.17 5.44 0.94 5.44 1.13 3.4 0.638 4.13 0.596 31 263.544 210.83 21.714 9.71 0.14 山西组 S-1 92.3 177 19.9 77.4 15.1 3.48 10.8 1.64 9.15 1.89 5.61 1.02 6.51 0.94 50.1 472.84 385.18 37.56 10.26 1.22 S-2 86.1 161 19.1 73.1 14.4 2.65 10.2 1.45 7.36 1.49 4.4 0.767 4.62 0.691 44.3 431.628 356.35 30.978 11.50 0.11 S-3 87.4 163 19.1 65 5.41 0.429 5.15 0.738 4.14 0.962 3.27 0.633 4.03 0.555 23.6 383.417 340.339 19.478 17.47 S-4 49.8 91.5 10.8 39.4 6.13 0.95 5.43 0.979 6.09 1.34 3.99 0.777 5 0.707 33.9 256.793 198.58 24.313 8.17 0.13 S-5 113 201 24.1 97.9 16.2 2.41 11.9 1.89 9.81 2.17 6.53 1.13 7.31 1.08 65 561.43 454.61 41.82 10.87 S-6 70.3 134 15 54.8 7.95 0.96 6.24 0.983 5.09 1.06 3.26 0.586 3.77 0.53 28.6 333.129 283.01 21.519 13.15 0.26 S-7 48.7 82.7 9.9 40.7 11 1.61 6.37 0.931 4.93 1.02 2.9 0.508 3.33 0.491 27.3 242.39 194.61 20.48 9.50 S-8 68.9 130 15 54.1 8 1.23 6.21 0.953 5.29 1.14 3.71 0.648 4.35 0.616 31.9 332.047 277.23 22.917 12.10 S-9 35.7 63.6 7.59 33.5 10.5 1.52 5.84 0.894 4.45 0.855 2.31 0.422 2.76 0.394 23 193.335 152.41 17.925 8.50 S-10 54.8 104 12 45 8.1 1.45 6.86 1.09 5.9 1.19 3.56 0.608 3.95 0.582 32.8 281.89 225.35 23.74 9.49 1.04 S-11 74.5 144 17 65.5 11.8 2.01 10.5 1.75 9.21 1.75 5.01 0.844 5.39 0.773 48.6 398.637 314.81 35.227 8.94 1.84 S-12 93.5 187 21.8 84.4 15.3 2.44 11.4 1.75 9.49 1.93 5.81 1.02 6.34 0.904 50.6 493.684 404.44 38.644 10.47 1.46 S-13 81.8 165 19 72.6 12.9 2.14 10.8 1.76 9.27 1.85 5.51 0.966 6.23 0.913 51.3 442.039 353.44 37.299 9.48 1.94 太原组 T-1 69.2 138 16 60.8 11.4 1.66 8.69 1.4 7.73 1.52 4.45 0.77 4.79 0.681 40.1 367.191 297.06 30.031 9.89 T-2 73.2 143 16.5 63 11.4 1.91 9.84 1.6 8.77 1.68 4.77 0.803 5.15 0.721 45.2 387.544 309.01 33.334 9.27 2.37 T-3 79.9 151 17.4 64.8 11.2 1.67 9.47 1.55 8.35 1.64 4.82 0.826 5.23 0.746 46.3 404.902 325.97 32.632 9.99 T-4 59.1 117 13.5 52 8.86 1.86 7.52 1.15 5.62 1.04 2.93 0.476 3.02 0.428 29.5 304.004 252.32 22.184 11.37 2.18 T-5 66.1 121 13.4 50.2 8.67 1.8 7.28 1.18 6.51 1.27 3.75 0.651 4.2 0.623 37.5 324.134 261.17 25.464 10.26 1.57 T-6 64.3 117 13.4 49.2 8.08 1.51 6.9 1.1 5.97 1.16 3.39 0.593 3.88 0.544 34.1 311.127 253.49 23.537 10.77 0.97 本溪组 B-1 6.8 18 1.97 8.5 2.69 0.678 3.38 0.939 6.85 1.4 4.06 0.772 5 0.72 37.3 99.059 38.638 23.121 1.67 0.13 B-2 15.3 122 8.92 39.7 7.95 0.801 4.71 0.871 5.88 1.31 4.32 0.786 5.25 0.769 31.3 249.867 194.671 23.896 8.15 0.32 B-3 40.4 96.3 8.6 36 9.17 1.97 6.38 1.08 6.24 1.25 3.65 0.704 4.93 0.726 31.8 249.2 192.44 24.96 7.71 0.52 北美页岩 32.00 73.00 7.90 33.00 5.70 1.24 5.20 0.65 5.80 1.04 3.40 0.50 3.10 0.48 27.00 173.21 152.84 20.37 7.50 球粒陨石 0.31 0.81 0.12 0.60 0.20 0.07 0.26 0.05 0.32 0.07 0.21 0.32 0.21 0.03 1.96 3.58 2.11 1.48 1.43 PAAS 38.2 79.6 8.83 33.09 5.55 1.08 4.66 0.774 4.68 0.991 2.85 0.405 2.82 0.433 注：表中 稀土元素的数据单位为10−6，TOC的单位为%，LREE表示轻稀土元素含量（La-Eu），HREE表示重稀土元素含量（Gd-Lu），∑REE=LREE+HREE+Y元素含量，L/H=LREE/HREE。 ![]()
图 2 南华北盆地太康隆起上古生界泥岩球粒陨石标准化稀土元素分配模式a. 下石盒子组，b. 山西组，c. 太原组，d. 本溪组。Figure 2. Chondrite-normalized REE patterns of mudstone samples in Upper Paleozoic, Taikang Uplift, southern North China Basin3. 稀土地球化学特征
3.1 含量特征
测试共检测出La、Ce、Pr、Eu、Y等15种稀土元素（表1）。测试结果表明，本区上古生界本溪组、太原组、山西组和下石盒子组泥岩稀土元素具有不同的地球化学特征。本溪组泥岩稀土总量（∑REE）明显相对较小，为99.06×10−6～249.87×10−6，平均值199.38×10−6；太原组∑REE为304.00×10−6～404.90×10−6，平均值349.82×10−6；山西组∑REE平均值略高于太原组，为371.02×10−6，但波动幅度明显较大，为193.34×10−6～561.43×10−6，最高值与最低值之差达368.09×10−6；下石盒子组∑REE为87.18×10−6～386.12×10−6，平均值243.58×10−6，相较山西组有所下降，但仍明显高于本溪组（表2）。
表 2 太康隆起上古生界不同层段泥岩稀土元素质量分数Table 2. Comparison of REE contents in mudstone from different formations on Taikang Uplift层位 ∑REE ∑LREE ∑HREE 范围 均值 范围 均值 范围 均值 下石盒子组 87.18～386.12 243.58 53.64～322.25 194.94 13.84～27.07 20.57 山西组 193.34～561.43 371.02 152.41～454.61 303.1 17.93～41.82 28.61 太原组 304.00～404.90 439.82 252.32～325.97 283.17 22.18～33.33 27.86 本溪组 99.06～249.87 199.38 38.64～194.67 141.92 23.12～24.96 23.99 上古生界泥岩整体表现出明显的REE富集特征，高于大陆上地壳（UCC）（稀土平均含量∑REE为146.4×10−6），也高于北美页岩（173.21×10−6）。垂向上本溪组和太原组之间存在一期REE富集程度加剧事件，太原组和山西组泥岩存在REE超常富集现象（ΣREE分别为UCC标准值的2.12倍和2.27倍），至下石盒子组该趋势波动下降。各稀土元素在垂向上的变化规律与∑REE基本一致，表明沉积过程中具有相同的物源[10]。
3.2 地球化学参数
研究沉积物稀土元素分配模式可以采用球粒陨石标准化[11]、北美页岩（NASC）标准化[12]和澳大利亚后太古宙平均页岩（PAAS）[13]标准化等多个方式。本文采用球粒陨石标准化研究样品相对地球原始物质的分异程度，以此讨论沉积物源区特征；采用北美页岩标准化研究沉积过程中的分异程度。
轻、重稀土元素比值∑LREE/∑HREE可以在一定程度上反映二者的分异程度。比值越大，轻稀土越富集，重稀土越亏损。研究区本溪组、太原组、山西组和下石盒子组4个层位的∑LREE/∑HREE平均值分别为5.84、10.26、10.76和9.00，具有明显的轻稀土元素富集、重稀土元素亏损的特征。轻稀土元素的变化特征基本决定了稀土总量的变化特征。除本溪组外，太原组、山西组和下石盒子组的∑LREE/∑HREE比值与上地壳的∑LREE/∑HREE比值（9.5）相近，与上地壳稀土元素分异情况相似，表明沉积过程中主要受到陆源碎屑物质的影响。
（La/Yb）N和（Ce/Yb）N可以反映稀土元素标准化图解中曲线的倾斜程度。经球粒陨石标准化后，研究区本溪组泥岩（La/Yb）N为0.92～5.55，平均值为2.82，（Ce/Yb）N为0.93～6.02，平均值为4.01；太原组泥岩（La/Yb）N为9.63～13.26，平均值为10.82，（Ce/Yb）N为7.20～10.04，平均值为7.91；山西组泥岩（La/Yb）N为6.75～14.69，平均值为10.29，（Ce/Yb）N为4.74～10.49，平均值为7.39；下石盒子组泥岩（La/Yb）N为2.46～10.73，平均值为7.30，（Ce/Yb）N为2.21～7.97，平均值为5.67。4个层位泥岩样品的轻、重稀土元素分异程度均较大，其中太原组和山西组分异程度相近，明显高于本溪组和下石盒子组（表3）。
表 3 太康隆起上古生界泥岩稀土元素地球化学参数Table 3. REE geochemical parameters of samples from Upper Paleozoic, Taikang Uplift层位 样品
编号La/
CeLa/
YbLa/
SmGd/
YbSm/
Nd（La/
Sm）N（La/
Yb）N（Ce/
Yb）N（La/
Lu）N（Gd/
Yb）N（La/
Yb）SδCeS δEuS δCeN δEuN δCeanom 下石盒子组 X-1 0.45 11.88 6.06 1.35 0.15 3.91 8.05 6.89 8.21 1.09 1.15 0.94 0.81 1.02 0.60 0.00 X-2 0.43 3.63 3.59 0.98 0.23 2.32 2.46 2.21 2.31 0.79 0.35 0.92 1.71 1.01 1.27 0.01 X-3 0.52 13.30 6.09 1.84 0.19 3.93 9.01 6.59 8.69 1.49 1.29 0.89 1.03 0.97 0.76 −0.05 X-4 0.51 15.84 6.58 1.71 0.17 4.25 10.73 7.97 10.62 1.38 1.53 0.88 0.84 0.96 0.62 −0.04 X-5 0.44 7.17 7.92 0.81 0.12 5.11 4.86 4.22 4.79 0.65 0.69 0.93 0.64 1.02 0.48 0.01 X-6 0.54 12.86 8.35 1.32 0.16 5.39 8.71 6.17 8.62 1.06 1.25 0.87 0.87 0.96 0.65 −0.05 山西组 S-1 0.52 14.18 6.11 1.66 0.20 3.94 9.60 7.05 9.50 1.34 1.37 0.90 1.20 0.98 0.89 −0.04 S-2 0.53 18.64 5.98 2.21 0.20 3.86 12.62 9.03 12.06 1.78 1.81 0.86 0.96 0.95 0.71 −0.06 S-3 0.54 21.69 16.16 1.28 0.08 10.42 14.69 10.49 15.24 1.03 2.10 0.87 0.36 0.95 0.26 −0.05 S-4 0.54 9.96 8.12 1.09 0.16 5.24 6.75 4.74 6.82 0.88 0.96 0.86 0.72 0.94 0.54 −0.06 S-5 0.56 15.46 6.98 1.63 0.17 4.50 10.47 7.13 10.13 1.31 1.50 0.84 0.76 0.92 0.57 −0.08 S-6 0.52 18.65 8.84 1.66 0.15 5.71 12.63 9.22 12.84 1.34 1.81 0.90 0.60 0.98 0.44 −0.04 S-7 0.59 14.62 4.43 1.91 0.27 2.86 9.91 6.44 9.60 1.55 1.42 0.82 0.84 0.90 0.63 −0.10 S-8 0.53 15.84 8.61 1.43 0.15 5.56 10.73 7.75 10.82 1.15 1.53 0.88 0.77 0.96 0.57 −0.05 S-9 0.56 12.93 3.40 2.12 0.31 2.19 8.76 5.97 8.77 1.71 1.25 0.84 0.85 0.92 0.63 −0.08 S-10 0.53 13.87 6.77 1.74 0.18 4.36 9.40 6.83 9.11 1.40 1.34 0.88 0.85 0.97 0.63 −0.05 S-11 0.52 13.82 6.31 1.95 0.18 4.07 9.36 6.93 9.33 1.57 1.34 0.88 0.79 0.96 0.59 −0.04 S-12 0.50 14.75 6.11 1.80 0.18 3.94 9.99 7.65 10.01 1.45 1.43 0.90 0.81 0.99 0.60 −0.03 S-13 0.50 13.13 6.34 1.73 0.18 4.09 8.89 6.87 8.67 1.40 1.27 0.91 0.80 1.00 0.59 −0.03 太原组 T-1 0.50 14.45 6.07 1.81 0.19 3.92 9.79 7.47 9.83 1.47 1.40 0.90 0.73 0.99 0.54 −0.03 T-2 0.51 14.21 6.42 1.91 0.18 4.14 9.63 7.20 9.83 1.54 1.38 0.90 0.79 0.98 0.59 −0.04 T-3 0.53 15.28 7.13 1.81 0.17 4.60 10.35 7.49 10.36 1.46 1.48 0.88 0.71 0.96 0.53 −0.05 T-4 0.51 19.57 6.67 2.49 0.17 4.30 13.26 10.04 13.36 2.01 1.90 0.90 1.00 0.99 0.74 −0.03 T-5 0.55 15.74 7.62 1.73 0.17 4.92 10.66 7.47 10.27 1.40 1.52 0.89 0.99 0.97 0.74 −0.06 T-6 0.55 16.57 7.96 1.78 0.16 5.13 11.23 7.82 11.44 1.44 1.61 0.87 0.89 0.95 0.66 −0.06 本溪组 B-1 0.38 1.36 2.53 0.68 0.32 1.63 0.92 0.93 0.91 0.55 0.13 1.07 0.99 1.17 0.73 0.06 B-2 0.13 2.91 1.92 0.90 0.20 1.24 1.97 6.02 1.93 0.72 0.28 2.27 0.57 2.49 0.43 0.44 B-3 0.42 8.19 4.41 1.29 0.25 2.84 5.55 5.06 5.39 1.05 0.79 1.13 1.13 1.23 0.84 0.05 注：表中数据单位为10−6；N表示采用球粒陨石标准化（据Boynton, 1984），S表示采用北美页岩标准化（据Haskin等,1968），P表示采用PAAS标准化（据McLenenan, 1989）；δCe=CeN/（LaN×PrN）1/2；δEu=EuN/（SmN×GdN）1/2；δCeanom=lg[3CeN/（2LaN+ωNdN）]。 经球粒陨石标准化后可知，除本溪组泥岩样品表现出明显的Ce正异常以外（δCeN平均值为1.63），太原组、山西组和下石盒子组泥岩样品均表现为极弱的Ce负异常（δCeN平均值分别为0.97、0.95和0.99），相对于球粒陨石分异不明显。4个层位的泥岩样品相对于球粒陨石均表现为明显的Eu负异常（δEuN平均值由下至上分别为0.67、0.63、0.59和0.73）（表3）。
北美页岩标准化的下石盒子组（La/Yb）S的平均值为1.04，与北美页岩（1.00）基本相近；太原组和山西组的（La/Yb）S的平均值分别为1.55和1.47，明显高于北美页岩，表明轻、重稀土元素分异程度较大，轻稀土元素相对富集；本溪组的（La/Yb）S的平均值为0.40，表明轻稀土元素亏损，重稀土元素相对富集。δCeS由下至上平均值分别为1.49、0.89、0.87和0.91，即太原组、山西组和下石盒子组泥岩样品相对于北美页岩整体呈较弱的Ce负异常，本溪组为较弱的正异常。δEuS由下至上平均值分别为0.90、0.85、0.79和0.98，即本溪组、太原组和山西组泥岩样品相对于北美页岩整体呈较弱的Eu负异常，下石盒子组泥岩样品相对于北美页岩分异不明显（表3）。
3.3 分布模式
由球粒陨石标准化稀土元素配分模式图（图2）可以看出，从太原组到下石盒子组的配分曲线均显示La-Sm段曲线斜率相对较大，Gd-Er段曲线平缓，Eu和Tm呈谷状负异常，表明轻稀土元素富集且分异程度较大，重稀土元素相对亏损，分异程度较低，基本为典型的沉积岩稀土配分特征。
北美页岩标准化元素配分模式图显示，本溪组样品的REE配分曲线表现为较为平缓的左倾型（图3d），表明本溪组轻稀土元素相对亏损，重稀土元素相对富集，Ce异常变化范围较大，整体表现为较弱的Ce负异常（δCeS：0.14～1.97，平均值0.90），并具有明显的Eu负异常（δEuS：0.21～0.61，平均值0.40）。太原组、山西组和下石盒子组的REE分布特征均表现为平坦的NASC标准化配分模式，其中山西组相比太原组具有较弱的Ce负异常和较强的Eu负异常（图3b、3c），下石盒子组，轻稀土元素的富集程度和轻、重稀土元素的分异程度均略低（图3a）。
![]()
图 3 南华北盆地太康隆起上古生界泥岩北美页岩标准化稀土元素分配模式a. 下石盒子组，b. 山西组，c. 太原组，d. 本溪组。Figure 3. NASC-normalized REE patterns of Upper Paleozoic mudstone samples, Taikang Uplift, southern North China Basin4. 稀土元素地质意义
4.1 成岩作用影响的讨论
前人研究表明，成岩作用对泥岩REE存在一定的影响，可能改变Ce的异常值，使δCeN与δEuN呈良好的负相关、δCeN与∑REE呈良好的正相关[14]。研究区4个层位泥岩样品的δCeN、δEuN和∑REE之间不存在明显的相关性（图4），表明成岩作用对样品REE的影响十分有限，δCeN与（La/Sm）N之间较差的相关性也表明样品的Ce异常可以代表原始沉积物的信息，利用稀土元素恢复古沉积环境具有较大的可靠度。
![]()
图 4 研究区泥岩δCeN-δEuN、δCeN-∑REE、δCeN-（La/Sm）N相关性图解Figure 4. Diagrams of δCeN-δEuN, δCeN-∑REE and δCeN-（La/Sm）N for Upper Paleozoic mudstone samples of the study area4.2 稀土元素的古环境意义
4.2.1 氧化还原性
沉积物的Ce异常可以反映沉积水体的氧化—还原条件。前文已述δCeS由下至上平均值分别为1.49、0.89、0.87和0.91，即太原组、山西组和下石盒子组泥岩样品相对于北美页岩整体呈较弱的Ce负异常，本溪组为较弱的正异常。一般认为δCeS负异常是海相沉积环境的指标，但在边缘海、浅海及封闭海域中Ce异常一般不明显，不发生亏损，在开阔海域则会发生严重亏损[15]。因此，可以认为研究区泥岩受到一定程度的海水影响，太原期和山西期影响更大，整体为弱还原—弱氧化环境。
此外，Ce异常指数还可以用δCeanom表示，依据Elderfield H等提出的定义，其计算公式为：
δCeanom=lg[3ω（Ce）N/（2ω（La）N+ω（Nd）N）]
当δCeanom＜−0.1时，表明Ce亏损，反映沉积水体为氧化环境；当δCeanom＞−0.1时，表明Ce富集，反映沉积水体为还原环境[16]。
研究区本溪组δCeanom为0.05～0.44，平均值为0.18；太原组δCeanom为−0.06～−0.03，平均值为−0.04；山西组δCeanom为−0.08～−0.03，平均值为−0.05；下石盒子组δCeanom为−0.05～0.01，平均值为−0.02。4个层位的δCeanom均大于−0.1，表明各段地层氧化还原条件差别不大，整体形成于缺氧的还原环境。其中，太原组沉积时海侵范围扩大，沉积水体变深，沉积环境整体比较稳定，更利于有机质的保存，进而形成富有机质泥岩；山西组含煤地层沉积期发生过多期区域性海侵海退事件，沉积环境相对动荡，沉积物阶段性暴露于地表，受到一定程度的氧化，导致呈δCeanom弱负异常，并且波动范围明显较大。
4.2.2 沉积速率
前人研究表明，REE的分异程度可以反映沉积速率[15]。由于各稀土元素在电价、吸附性等方面的差异，稀土元素会随沉积水体的变化发生分异，主要表现为轻、重稀土元素之间的分异，以及铈元素和铕元素与其他元素之间的分异。这是因为在海（湖）水环境中，黏土等细碎屑质和悬浮物是吸附REE的主要载体，当悬浮物在海水中停留时间较长时，REE有充分的时间被分解，进入海水的REE也有充足的时间被细粒物质或悬浮物吸附，导致REE强烈分异，沉积物的轻、重稀土元素相应地出现富集或亏损，Ce也会随沉积水体的氧化还原条件发生选择性分异。当沉降速率较快时，REE随吸附物快速沉降，与海水发生交换的机会少，其分异程度相应较弱，沉积物的REE页岩标准化分配曲线平缓，Ce分异不明显或呈弱负异常，曲线斜率（可以用（La/Yb）S的值表示）一般在1左右。
由（La/Yb）S来指示的REE分异程度可以看出，本溪组的沉积速率明显较高，（La/Yb）S平均值为0.40；太原组样品的（La/Yb）S为1.38～1.90，平均值为1.55，反映了太原组的沉积速率显著降低，为低能、缺氧环境，与太原组沉积期海侵面积扩大，沉积水体变深的构造背景相符；山西组沉积速率波动幅度大，整体与太原组相近，（La/Yb）S为0.96～2.10，平均值为1.47；下石盒子组（La/Yb）S为0.35～1.53，平均值为1.04，表明沉积速率随环境变化再次升高。各层段残余有机碳含量（TOC）的差别也证实了这一点。水体中黏土等细碎屑质和悬浮物除了易吸附REE以外，也是吸附有机质的主要载体，因此有机质具有随沉积速率变大而减小的趋势[17]。本溪组沉降速率大，泥岩TOC平均值为0.32（n=3）；太原组和山西组泥岩TOC则分别为1.77（n=4）和1.00（n=8），表明沉积速率较低时，有机质能够充分被吸附，并随细粒物质的沉降而保存下来。这一特征在垂向剖面上也很明显，太原组和山西组泥岩样品的TOC高值与（La/Yb）S低值有良好的对应关系（图5）。
![]()
图 5 太康隆起尉参1井稀土元素质量分数及地球化学参数垂向分布图Figure 5. Vertical distribution of REE contents and geochemical parameters of Well Weican-1 in Taikang Uplift4.3 沉积源区构造背景的讨论
稀土元素由于其自身稳定的化学性质，后期风化作用、沉积分选、成岩作用和蚀变作用对沉积物中稀土元素的组成影响较弱，可以在一定程度上保存源岩的稀土特征[18]。REE特征相近代表物源相同或存在一定的相关性，因此，可以利用沉积岩稀土元素特征对物源区构造环境进行讨论。
从本区样品球粒陨石标准化的稀土元素分配与典型构造背景模式的对比可知，太原组、山西组和下石盒子组典型泥岩样品的REE分布曲线特征整体相近，显示物源具有大陆边缘构造背景的亲和性；本溪组典型泥岩样品以轻微的Ce正异常和Eu异常不明显为特征区别于另外3个层位，但整体仍具有大陆边缘构造背景的亲和性（图6a）。本次测试中所有样品均出现显著的铥（Tm）负异常，与常规稀土元素分配模式不符。铥元素是稀土元素中含量最少的元素，可以认为是实验误差所致。由于其化学性质稳定，不涉及关键化学参数，不影响其他元素数据的使用。
![]()
图 6 太康隆起上古生界泥岩源岩构造背景判别a. 泥岩球粒陨石标准化稀土元素分布曲线与典型构造背景砂岩的对比，b. （La/Yb）N-ΣREE判别图 （底图据Bhatia，1985），c. 泥岩δEuP-（Gd/Yb）P散点图 （底图据McLennan等，1993[20]）；d. 泥岩Ni-Cr散点图 （底图据Taylor等，1986[21]）。Figure 6. Tectonic settings of the mudstone samples in Upper Paleozoic, Taikang Uplifta. comparison of chondrite normalized REE patterns of mudstone with various tectonic settings，b.（La/Yb）N-ΣREE diagram （base map after Bhatia, 1985），c. Scattered diagram of δEup-（Gd/Yb）p of the mudstones samples （base map after McLennan et al. 1993），d. Scattered diagram of Ni-Cr of the sampled mudstones （base map after Taylor et al. 1986）Sm/Nd比值也可以反映物质来源。Sm和Nd均为不相容元素，通常在地质地壳作用过程中不发生分异。不同性质地球物质的Sm/Nd比值变化范围较小，如地幔为0.260～0.375，大洋玄武岩为0.234～0.425，海水为0.211，页岩为0.209。本区除本溪组样品的Sm/Nd比值相对较高以外（平均值为0.26），太原组、山西组和下石盒子组泥岩样品的Sm/Nd比值相近，主要分布在0.16～0.23，平均值分别为0.17、0.18和0.17，反映物源为上地壳。
(La/Yb）N-ΣREE判别图可以进一步揭示源岩的属性[19]。太原组和山西组样品分布点高度集中，大部分落在花岗岩区和花岗岩与沉积岩重叠区，山西组另有部分样品落在沉积岩区，表明太原期和山西期具有较为统一的物源，沉积物源主要来自上地壳长英质源区，但二者之间仍存在一定的差异。本溪组样品主要落在玄武岩区，表明物源受海相环境控制；下石盒子组样品点较为分散，横跨多个区域，揭示了较为动荡的沉积环境和混杂物源，与沉积速率的变化相符。
稀土元素比值中，δEuP和（Gd/Yb）P可以用来判定源岩性质和形成地质年代。一般认为显生宙形成的花岗岩类由于富含富钾长石，通常具有Eu弱亏损、HREE强烈亏损和（Gd/Yb）P＜2的特征。δEuP-（Gd/Yb）P图解表明，研究区石炭—二叠系泥岩样品的（Gd/Yb）P为0.41～1.51，落在后太古宇区域；山西组泥岩样品中δEuP＜0.85的样品相对其他组明显较多，表明山西组有部分物源是形成于后太古宙的花岗岩类，与（La/Yb）N-ΣREE图解所示结论相符。结合微量元素Ni-Cr图解，可以进一步验证太原组、山西组和下石盒子组泥岩的母岩主要形成于后太古宙，因为后太古宙泥岩相较于太古宙具有更高的长英质组分含量，而镁铁质元素，尤其是Ni和Co的含量更低[22]。
区域资料显示，晚石炭世中期本溪组沉积时，南华北地区为南、北、西三面高，向东、东北倾斜的斜坡环境，沉积中心在研究区东北方的邯郸—肥城以北地区。本次研究采集的本溪组样品均位于本溪组顶层，对应晚石炭世末期，该时期海水退缩。早二叠世初再次发生海侵，由东向南、向西的海侵范围进一步扩大，沉积中心来到研究区所在的郑州—临沂一带，沉积环境为典型的浅水陆表海，发育海相灰岩和煤层交互出现的沉积序列，这一点在尉参1井钻遇地层岩性上有所体现。至早二叠世末，受兴蒙洋向华北板块俯冲加剧的影响，华北北缘的古阴山再度上升，沉积范围大规模向东南退缩。受海西运动影响，盆地整体上表现为由弱伸展背景的克拉通内坳陷向克拉通内陆（含煤）坳陷转变。
结合上古生界泥岩稀土元素对物源属性的指示和已有的区域构造演化认识，可以认为太康隆起地区太原组和山西组具有不同的物源亚区，太原组沉积物源可能来自秦岭造山带中的华北板块南缘带，山西组沉积物源可能还包括北秦岭构造带[23]，母岩物质均源自上地壳，为花岗岩和沉积岩混合物源。
4.4 稀土元素与有机质富集关系的讨论
稀土元素在垂向上的变化规律可以反映沉积水体的变化，而沉积水体成分的变化又受到古气候的直接影响，因此，可以用稀土元素的变化特征来讨论古气候的演变[24]。一般认为，ΣREE的高值代表陆源碎屑输入大，指示温湿气候条件，低值代表陆源输入少，指示干旱气候[10]。研究区各层位样品的ΣREE与TOC均具有一定的正相关性（图7），垂向上TOC高值段与ΣREE高值段对应良好（图5），表明稀土元素丰度高的泥岩主要对应富有机质层段。根据有机质的形成与保存条件可知，残余有机碳质量分数高的层段指示温暖湿润的气候环境，低值段指示干冷或干燥的气候环境[25]。
![]()
图 7 研究区泥岩TOC-∑REE相关性图解Figure 7. Diagrams of TOC-∑REE for mudstone samples of the study area根据ΣREE、(La/Yb）S和TOC垂向上的变化规律，结合区域沉积背景，尉参1井本溪组顶部至下石盒子组底部的连续剖面可以分为5个层段（图5）。
（1）本溪组顶部铝土泥岩段沉积时期（对应井段2 809～2 822 m），∑REE值和（La/Yb）S明显较低，表明沉积水体浅，沉降速率大，陆源物质注入不丰富，表现出相对干旱的气候特征。有机碳含量也相应较低，为全井段最低处，本溪组顶部泥岩宏观表现为浅灰色—灰色。
（2）向上至太原组上部泥岩段沉积初期（2 774～2 809 m），∑REE值和（La/Yb）S稳步增大，且二者变化趋势相同，表明沉积水体加深的同时，物源注入加大，沉降速率降低致使REE发生充分分异，表现为温暖湿润的气候条件。有机质在这一时期得到充分的富集，表现为在太原组泥岩的TOC为全井段碳质泥岩和煤以外的最大值，向上随（La/Yb）S的减小稳定降低。
（3）太原组顶部碳质泥岩到山西组中部煤层以下的泥岩连续发育层段（2 718～2 774 m），∑REE值和（La/Yb）S不再具有相同的变化趋势，太原组顶部和山西组底部灰色泥岩夹灰岩层段，∑REE值整体较高，波动较明显，但（La/Yb）S曲线平缓无波动，小于下伏层段，表明物源丰富且稳定，沉降速率升高，稀土元素分异不明显，有机质含量最高值低于下伏层段，但整体较高，基本保持在约2%，指示温暖湿润的气候环境。
（4）山西组中部灰黑色含煤泥页岩段（2 702～2 718 m）∑REE值和（La/Yb）S变化趋势相近，波动明显但范围较小，∑REE值整体较下伏地层略有降低，（La/Yb）S值整体较下伏地层略有升高，表明物源丰富，沉降速率大，水平面相对下伏地层沉积时期略有降低，泥岩TOC含量较高并发育3 m厚煤层，结合区域沉积背景，代表了温湿气候条件下浅水三角洲沉积体系中的分流间沼泽环境。
（5）山西组上部碎屑岩层段（俗称大占砂岩、香炭砂岩和小紫泥岩段）（对应井深2 670～2 702 m），∑REE值和（La/Yb）S变化趋势一致，波动范围大，整体均高于下伏泥岩段，但TOC降低，表明沉积环境较为动荡，结合区域沉积背景，可以认为该时期为经历了一次海侵—海退变化的三角洲平原环境。
5. 结论
（1）太康隆起上古生界泥岩中稀土元素总量较高，高于大陆上地壳（UCC）和北美页岩。∑LREE/∑HREE、（La/Yb）N和（Ce/Yb）N等化学参数表明本溪组稀土元素特征与太原组、山西组和下石盒子组明显不同，其∑REE相对较低，稀土元素分异程度较低，具有明显的δCeS正异常和（La/Yb）S负异常，轻稀土元素相对亏损，重稀土元素相对富集。太原组和山西组稀土元素总量高，轻、重稀土元素分异程度相近，明显高于本溪组和下石盒子组，分配模式为轻稀土元素富集且分异程度较大，重稀土元素相对亏损且分异程度较低，具有较弱的Ce负异常和较弱的Eu负异常。
（2）泥岩的Ce异常值反映本溪组至下石盒子组氧化还原条件相近，整体形成于缺氧的还原环境，受到一定程度的海水影响。太原组沉积时海侵范围扩大，沉积水体变深，沉积环境整体比较稳定；山西组含煤地层沉积期发生过多期区域性海侵海退事件，沉积环境相对动荡。
（3）稀土元素分布模式、Sm/Nd比值和（La/Yb）N-ΣREE判别图等显示太原组和山西组样品分布点高度集中，大部分落在花岗岩区和花岗岩与沉积岩重叠区，沉积物源主要来自上地壳长英质源区。结合δEuP-（Gd/Yb）P图解可以认为太原组和山西组泥岩的母岩形成于后太古宙，主要为活动性大陆边缘抬升基地物源。
（4）（La/Yb）S 可以表示REE的分异程度，进一步反映沉积速率。本溪组的（La/Yb）S平均值最小，太原组和山西组明显增大，且山西组的（La/Yb）S波动范围相对更大，下石盒子组的（La/Yb）S平均值减小，表明沉积速率先降低再增大。结合泥岩TOC在垂向上的变化特征可知，研究层位自下而上沉积时期气候条件经历了由冷干向温湿的转变，太原组沉降速率低，沉积环境稳定，气候条件温湿，是最有利于有机质富集和保存的层段。
-
![]()
图 1 研究区概况和IODP U1482站位图
a： 印澳地区现代表层海水温度[28]与洋流方向[11,14]；LC=Leeuwin Current (鲁汶流)；使用Ocean Data View制图[29]； b：2023年1月表层海水叶绿素浓度；c： 2022年7月表层海水叶绿素浓度[30]。
Figure 1. Regional oceanographic settings and the location of IODP Site U1482
a:Map of modern sea surface temperature[28] and ocean currents[11,14] in the Indo-Pacific region, created using Ocean Data View[29], with LC representing the Leeuwin Current; b: sea surface chlorophyll concentration in January, 2023; c: sea surface chlorophyll concentration in July, 2022[30].
![]()
图 2 IODP U1482站2～1.07 Ma期间碳酸盐含量与其他指标对比
a： 碳酸盐含量，b： log(Ca/Ti)， c：钾含量， d： log(Terr/Ca)， e： log(Zr/Rb)， f： 铀含量， g： 底栖有孔虫δ13C [14]， h： 底栖有孔虫δ18O [14]。
Figure 2. Comparison of carbonate content and other proxies from IODP Site U1482 during 2～1.07 Ma
a: Carbonate content, b: log(Ca/Ti), c: potassium content, d: log(Terr/Ca), e: log(Zr/Rb), f: uranium content, g: benthic foraminiferal δ13C[14], h: benthic foraminiferal δ18O[14].
![]()
图 3 IODP U1482站2～1.07 Ma期间碳酸盐含量与钾、铀元素记录的周期相关性以及相位关系
a：碳酸盐含量频谱(灰色)与钾元素含量频谱（黑色）， b： 碳酸盐含量频谱（灰色）与铀元素含量频谱（黑色），c： 碳酸盐含量与钾元素含量的频谱相关性与相位关系，d： 碳酸盐含量与铀元素含量的频谱相关性与相位关系。频谱分析采用Redfit-X软件[37]。
Figure 3. Spectral coherence and phase relationship between carbonate content and elemental records of potassium and uranium from IODP Site U1482 during 2～1.07 Ma
a:Spectra of carbonate content (gray) and potassium content (black), b: spectra of carbonate content (gray) and uranium content (black), c: coherence and phase relationship between carbonate content and potassium content,d: coherence and phase relationship between carbonate content and uranium content. Spectral analyses were performed using Redfit-X[37].
![]()
图 4 IODP U1482站有孔虫壳体和碎片扫描电镜照片[8]
a—c：Planulina Wuellerstorfi， d—f：Trilobatus Trilobus。样品取自A孔岩芯145.12～145.17 m CCSF处。
Figure 4. Scanning electron microscope images of foraminiferal tests and fragments of IODP Site U1482
Tests and fragments of Planulina Wuellerstorfi (a-c) and Trilobatus Trilobus (d-f) are retrieved from Hole A between depths of 145.12 m and 145.17 m CCSF.
![]()
图 5 IODP U1482站2～1.07 Ma期间碳酸盐、钾和铀元素含量及其堆积速率
a： 碳酸盐含量（灰色线）与堆积速率（黑色线），b： 钾元素含量（灰色线）与堆积速率（黑色线），c： 铀元素含量（灰色线）与堆积速率（黑色线），d： U1482站线性沉积速率。
Figure 5. Mass accumulation rates and contents of carbonate, potassium, and uranium of IODP Site U1482 during 2～1.07 Ma
a: Carbonate content (grey) and mass accumulation rate (black), b: potassium content (grey) and mass accumulation rate (black), c: uranium content (grey) and mass accumulation rate (black), d: linear sedimentation rate at Site U1482.
![]()
图 6 IODP U1482站2～1.07 Ma期间碳酸盐含量和陆源输入记录与岁差的相位关系
0相位设在岁差最大值。深灰色箭头长度与外圈半径之比为各记录在23 ka周期上与岁差的相关度。浅灰色扇形指代相位误差。6月21日夏季半球太阳辐射梯度为23°N与23°S之间太阳辐射量的差值。交叉频谱分析采用Redfit-X软件[37]。
Figure 6. Relationship in phase between carbonate content and terrestrial input records on the precessional band from IODP Site U1482 during 2～1.07 Ma
Phase zero was set at the precessional maxima. The ratio of the length of each vector to the radius of the outer circle corresponds to coherency. Shadings denote uncertainties of phase relationship. Summer intertropical insolation gradient on June 21 is the difference of June 21 insolation between at 23°N and 23°S. Cross-spectral analyses were performed using Redfit-X[37].
![]()
图 7 U1482站和邻近站位古海洋记录与太平洋经向和纬向温度梯度对比
a：U1482站与U1483站钾元素含量对比[48]，b： U1482站铀元素含量与U1483站总氮含量对比[48]，c： U1460站[53]与U1463表层海水温度[54]记录，d：赤道太平洋纬向温度梯度[55]，e： 南海经向温度梯度[52,56]。
Figure 7. Comparison of paleoceanographic records of Site U1482 and its neighboring sites with meridional and zonal temperature gradients of the Pacific Ocean
a: Comparison in potassium contents between Sites U1482 and U1483[48], b: comparison between U1482 uranium content and U1483 total nitrogen content[48], c: sea surface temperature records from Sites U1460[53] and U1463[54], d: zonal temperature gradient in the equatorial Pacific Ocean[55], e: meridional temperature gradient in the South China Sea[52,56].
-
[1] Zeebe R E, Westbroek P. A simple model for the CaCO3 saturation state of the ocean: the “Strangelove, ” the “Neritan, ” and the “Cretan” Ocean[J]. Geochemistry, 2003, 4(12): 1104.
[2] Ridgwell A. A Mid Mesozoic Revolution in the regulation of ocean chemistry[J]. Marine Geology, 2005, 217(3-4):339-357. doi: 10.1016/j.margeo.2004.10.036
[3] Zhang Y G, Pagani M, Liu Z H, et al. A 40-million-year history of atmospheric CO2[J]. Philosophical Transactions of the Royal Society A, 2013, 371(2001):20130096. doi: 10.1098/rsta.2013.0096
[4] Ridgwell A, Zeebe R E. The role of the global carbonate cycle in the regulation and evolution of the earth system[J]. Earth and Planetary Science Letters, 2005, 234(3-4):299-315. doi: 10.1016/j.jpgl.2005.03.006
[5] Raymo M E, Ruddiman W F, Froelich P N. Influence of late Cenozoic mountain building on ocean geochemical cycles[J]. Geology, 1988, 16(7):649-653. doi: 10.1130/0091-7613(1988)016<0649:IOLCMB>2.3.CO;2
[6] Berner R A, Lasaga A C, Garrels R M. The carbonate-silicate geochemical cycle and its effect on atmospheric carbon dioxide over the past 100 million years[J]. American Journal of Science, 1983, 283(7):641-683. doi: 10.2475/ajs.283.7.641
[7] Pagani M, Liu Z H, LaRiviere J, et al. High Earth-system climate sensitivity determined from Pliocene carbon dioxide concentrations[J]. Nature Geoscience, 2010, 3:27-30. doi: 10.1038/ngeo724
[8] Rosenthal Y, Holbourn A E, Kulhanek D K, et al. Western Pacific Warm Pool[M]. College Station: International Ocean Discovery Program, 2018.
[9] Straume E O, Nummelin A, Gaina C, et al. Climate transition at the Eocene–Oligocene influenced by bathymetric changes to the Atlantic–Arctic oceanic gateways[J]. Proceedings of the National Academy of Sciences of the United States of America, 2022, 119(17):e2115346119.
[10] Wei J L, Liu H L, Zhao Y, et al. Simulation of the climate and ocean circulations in the middle Miocene climate optimum by a coupled model FGOALS-g3[J]. Palaeogeography, Palaeoclimatology, Palaeoecology, 2023, 617:111509. doi: 10.1016/j.palaeo.2023.111509
[11] Gordon A L. Oceanography of the Indonesian seas[J]. Oceanography, 2005, 18(4):13. doi: 10.5670/oceanog.2005.18
[12] Cane M A, Molnar P. Closing of the Indonesian seaway as a precursor to east African aridification around 3-4 million years ago[J]. Nature, 2001, 411(6834):157-162. doi: 10.1038/35075500
[13] He Y X, Wang H Y, Liu Z H. Development of the Leeuwin Current on the northwest shelf of Australia through the Pliocene-Pleistocene period[J]. Earth and Planetary Science Letters, 2021, 559:116767. doi: 10.1016/j.jpgl.2021.116767
[14] Chen Y X, Xu J, Liu J, et al. Climatic and tectonic constraints on the Plio–Pleistocene evolution of the Indonesian Throughflow intermediate water recorded by benthic δ18O from IODP site U1482[J]. Quaternary Science Reviews, 2022, 295:107666. doi: 10.1016/j.quascirev.2022.107666
[15] Gordon A L, Susanto R D, Vranes K. Cool Indonesian throughflow as a consequence of restricted surface layer flow[J]. Nature, 2003, 425(6960):824-828. doi: 10.1038/nature02038
[16] Kuhnt W, Holbourn A, Hall R, et al. Neogene history of the Indonesian throughflow[M]//Clift P, Kuhnt W, Wang P, et al. Continent‐Ocean Interactions Within East Asian Marginal Seas. Washington: American Geophysical Union, 2004: 299-320.
[17] Talley L D, Sprintall J. Deep expression of the Indonesian Throughflow: Indonesian intermediate water in the South Equatorial Current[J]. Journal of Geophysical Research:Oceans, 2005, 110(C10):C10009.
[18] Sprintall J, Wijffels S E, Molcard R, et al. Direct estimates of the Indonesian Throughflow entering the Indian Ocean: 2004-2006[J]. Journal of Geophysical Research:Oceans, 2009, 114(C7):C07001.
[19] Rosenthal Y, Linsley B K, Oppo D W. Pacific Ocean heat content during the past 10, 000 years[J]. Science, 2013, 342(6158):617-621. doi: 10.1126/science.1240837
[20] 李铁刚, 熊志方, 贾奇. 晚中新世以来印度洋-太平洋暖池水体交换过程及其气候效应[J]. 海洋科学进展, 2020, 38(3):377-389 doi: 10.3969/j.issn.1671-6647.2020.03.001 LI Tiegang, XIONG Zhifang, JIA Qi. Water exchange between western Pacific warm pool and Indian warm pool and its climatic effects since the late Miocene[J]. Advances in Marine Science, 2020, 38(3):377-389.] doi: 10.3969/j.issn.1671-6647.2020.03.001
[21] Furnas M J, Carpenter E J. Primary production in the tropical continental shelf seas bordering northern Australia[J]. Continental Shelf Research, 2016, 129:33-48. doi: 10.1016/j.csr.2016.06.006
[22] Marin M, Feng M. Intra-annual variability of the North West Shelf of Australia and its impact on the Holloway Current: excitement and propagation of coastally trapped waves[J]. Continental Shelf Research, 2019, 186:88-103. doi: 10.1016/j.csr.2019.08.001
[23] Condie S A, Dunn J R. Seasonal characteristics of the surface mixed layer in the Australasian region: implications for primary production regimes and biogeography[J]. Marine and Freshwater Research, 2006, 57(6):569-590. doi: 10.1071/MF06009
[24] Pei R J, Kuhnt W, Holbourn A, et al. Monitoring Australian Monsoon variability over the past four glacial cycles[J]. Palaeogeography, Palaeoclimatology, Palaeoecology, 2021, 568:110280. doi: 10.1016/j.palaeo.2021.110280
[25] Sarim M, Xu J, Zhang P, et al. Late quaternary clay mineral and grain-size records from northwest Australia and their implications for paleoclimate, ocean currents, and paleodrainage of the Bonaparte basin[J]. Palaeogeography, Palaeoclimatology, Palaeoecology, 2023, 610:111353. doi: 10.1016/j.palaeo.2022.111353
[26] Hesse P P, McTainsh G H. Australian dust deposits: modern processes and the Quaternary record[J]. Quaternary Science Reviews, 2003, 22(18-19):2007-2035. doi: 10.1016/S0277-3791(03)00164-1
[27] Kuhnt W, Holbourn A, Schönfeld J, et al. Cruise report Sonne 257, WACHEIO - Western australian climate history from eastern Indian ocean sediment archives[R]. Darwin-Fremantle: Institut für Geowissenschaften, Christian-Albrechts-Universitat Kiel, 2017.
[28] Boyer T P, Garcia H E, Locarnini R A, et al. World Ocean Atlas 2018: sea surface temperature[DB/OL]. NOAA National Centers for Environmental Information. [2022-10-24]. https://www.ncei.noaa.gov/archive/accession/NCEI-WOA18.
[29] Schlitzer R. Ocean data view[CP/DK]. [2023-10-25]. https://odv.awi.de.
[30] NASA. SeaWiFS Mission page[DB/OL]. [2023-04-07]. https://oceancolor.gsfc.nasa.gov/.
[31] Murgese D S, de Deckker P. The Late Quaternary evolution of water masses in the eastern Indian Ocean between Australia and Indonesia, based on benthic foraminifera faunal and carbon isotopes analyses[J]. Palaeogeography, Palaeoclimatology, Palaeoecology, 2007, 247(3-4):382-401.
[32] de Deckker P, Barrows T T, Rogers J. Land–sea correlations in the Australian region: post-glacial onset of the monsoon in northwestern Western Australia[J]. Quaternary Science Reviews, 2014, 105:181-194. doi: 10.1016/j.quascirev.2014.09.030
[33] Holbourn A, Kuhnt W, Kawamura H, et al. Orbitally paced paleoproductivity variations in the Timor Sea and Indonesian Throughflow variability during the last 460 kyr[J]. Paleoceanography, 2005, 20(3):PA3002.
[34] Lo Giudice Cappelli E, Holbourn A, Kuhnt W, et al. Changes in Timor Strait hydrology and thermocline structure during the past 130 ka[J]. Palaeogeography, Palaeoclimatology, Palaeoecology, 2016, 462:112-124. doi: 10.1016/j.palaeo.2016.09.010
[35] Liu W, Baudin F, Moreno E, et al. Comparison of 240 ka long organic carbon and carbonate records along a depth transect in the Timor Sea: primary signals versus preservation changes[J]. Paleoceanography, 2014, 29(5):389-402. doi: 10.1002/2013PA002539
[36] de Vleeschouwer D, Dunlea A G, Auer G, et al. Quantifying K, U, and Th contents of marine sediments using shipboard natural gamma radiation spectra measured on DV JOIDES Resolution[J]. Geochemistry, Geophysics, Geosystems, 2017, 18(3):1053-1064. doi: 10.1002/2016GC006715
[37] Ólafsdóttir K B, Schulz M, Mudelsee M. REDFIT-X: cross-spectral analysis of unevenly spaced paleoclimate time series[J]. Computers & Geosciences, 2016, 91:11-18.
[38] 汪品先. 西太平洋边缘海的冰期碳酸盐旋回[J]. 海洋地质与第四纪地质, 1998, 18(1):1-11 WANG Pinxian. Glacial carbonate cycles in western Pacific marginal seas[J]. Marine Geology & Quaternary Geology, 1998, 18(1):1-11.]
[39] Feely R A, Sabine C L, Lee K, et al. In situ calcium carbonate dissolution in the Pacific Ocean[J]. Global Biogeochemical Cycles, 2002, 16(4):1144.
[40] Fischer G, Wefer G. Use of Proxies in Paleoceanography: Examples from the South Atlantic[M]. Berlin: Springer, 1999: 255-284.
[41] Sulpis O, Jeansson E, Dinauer A, et al. Calcium carbonate dissolution patterns in the ocean[J]. Nature Geoscience, 2021, 14(6):423-428. doi: 10.1038/s41561-021-00743-y
[42] Kuhnt W, Holbourn A, Xu J, et al. Southern Hemisphere control on Australian monsoon variability during the late deglaciation and Holocene[J]. Nature Communications, 2015, 6:5916. doi: 10.1038/ncomms6916
[43] Zhang P, Xu J, Holbourn A, et al. Obliquity induced latitudinal migration of the Intertropical Convergence Zone during the past ~410 kyr[J]. Geophysical Research Letters, 2022, 49(21):e2022GL100039. doi: 10.1029/2022GL100039
[44] Liu L W, Chen J, Chen Y, et al. Variation of Zr/Rb ratios on the Loess Plateau of Central China during the last 130000 years and its implications for winter monsoon[J]. Chinese Science Bulletin, 2002, 47(15):1298-1302. doi: 10.1360/02tb9288
[45] Schofield A. Uranium Content of Igneous Rocks of Australia 1:5000000 Maps—Explanatory Notes and Discussion[M]. Camberra: Geoscience Australia, 2009: 20.
[46] McManus J, Berelson W M, Klinkhammer G P, et al. Authigenic uranium: relationship to oxygen penetration depth and organic carbon rain[J]. Geochimica et Cosmochimica Acta, 2005, 69(1):95-108. doi: 10.1016/j.gca.2004.06.023
[47] Auer G, Hauzenberger C A, Reuter M, et al. Orbitally paced phosphogenesis in Mediterranean shallow marine carbonates during the middle Miocene Monterey event[J]. Geochemistry Geophysics Geosystems, 2016, 17(4):1492-1510. doi: 10.1002/2016GC006299
[48] Zhang Y, Andrade T, Ravelo A C, et al. Aridification of northwest Australia and nutrient decline in the Timor sea during the 40 Kyr world[J]. Paleoceanography and Paleoclimatology, 2023, 38(10):e2023PA004683. doi: 10.1029/2023PA004683
[49] Auer G, Petrick B, Yoshimura T, et al. Intensified organic carbon burial on the Australian shelf after the Middle Pleistocene transition[J]. Quaternary Science Reviews, 2021, 262:106965. doi: 10.1016/j.quascirev.2021.106965
[50] Lea D W, Pak D K, Spero H J. Climate impact of late quaternary equatorial pacific sea surface temperature variations[J]. Science, 2000, 289(5485):1719-1724. doi: 10.1126/science.289.5485.1719
[51] Etourneau J, Schneider R, Blanz T, et al. Intensification of the Walker and Hadley atmospheric circulations during the Pliocene-Pleistocene climate transition[J]. Earth and Planetary Science Letters, 2010, 297(1-2):103-110. doi: 10.1016/j.jpgl.2010.06.010
[52] Herbert T D, Peterson L C, Lawrence K T, et al. Tropical ocean temperatures over the past 3.5 million years[J]. Science, 2010, 328(5985):1530-1534. doi: 10.1126/science.1185435
[53] Petrick B, Martínez-García A, Auer G, et al. Glacial Indonesian throughflow weakening across the mid-Pleistocene climatic transition[J]. Scientific Reports, 2019, 9(1):16995. doi: 10.1038/s41598-019-53382-0
[54] Smith R A, Castañeda I S, Groeneveld J, et al. Retracted: Plio-Pleistocene Indonesian throughflow variability drove eastern Indian ocean sea surface temperatures[J]. Paleoceanography and Paleoclimatology, 2020, 35(10):e2020PA003872. doi: 10.1029/2020PA003872
[55] Wara M W, Ravelo A C, Delaney M L. Permanent El Niño-like conditions during the Pliocene warm period[J]. Science, 2005, 309(5735):758-761. doi: 10.1126/science.1112596
[56] Li L, Li Q Y, Tian J, et al. A 4-Ma record of thermal evolution in the tropical western Pacific and its implications on climate change[J]. Earth and Planetary Science Letters, 2011, 309(1-2):10-20. doi: 10.1016/j.jpgl.2011.04.016
[57] Stuut J B W, De Deckker P, Saavedra-Pellitero M, et al. A 5.3-million-year history of monsoonal precipitation in northwestern Australia[J]. Geophysical Research Letters, 2019, 46(12):6946-6954. doi: 10.1029/2019GL083035
[58] Christensen B A, Renema W, Henderiks J, et al. Indonesian Throughflow drove Australian climate from humid Pliocene to arid Pleistocene[J]. Geophysical Research Letters, 2017, 44(13):6914-6925. doi: 10.1002/2017GL072977
[59] Fedorov A V, Brierley C M, Lawrence K T, et al. Patterns and mechanisms of early Pliocene warmth[J]. Nature, 2013, 496(7443):43-49. doi: 10.1038/nature12003
[60] Chen H J, Xu Z K, Lim D, et al. Geochemical records of the provenance and silicate weathering/erosion from the eastern Arabian sea and their responses to the Indian summer monsoon since the mid-Pleistocene[J]. Paleoceanography and Paleoclimatology, 2020, 35(4):e2019PA003732. doi: 10.1029/2019PA003732
[61] Beaufort L, de Garidel-Thoron T, Mix A C, et al. ENSO-like forcing on oceanic primary production during the Late Pleistocene[J]. Science, 2001, 293(5539):2440-2444. doi: 10.1126/science.293.5539.2440
[62] Di Nezio P N, Timmermann A, Tierney J E, et al. The climate response of the Indo-Pacific warm pool to glacial sea level[J]. Paleoceanography, 2016, 31(6):866-894. doi: 10.1002/2015PA002890
[63] Windler G, Tierney J E, DiNezio P N, et al. Shelf exposure influence on Indo-Pacific Warm Pool climate for the last 450, 000 years[J]. Earth and Planetary Science Letters, 2019, 516:66-76. doi: 10.1016/j.jpgl.2019.03.038
[64] Liautaud P, Huybers P. Uniformitarian prediction of early-Pleistocene atmospheric CO2[J]. Geophysical Research Letters, 2022, 49(20):e2022GL100304. doi: 10.1029/2022GL100304
-
期刊类型引用(5)
1. 王超群，李凤杰，王佳. 川东北地区二叠系吴家坪组地球化学特征及其古环境意义. 天然气地球科学. 2025(01): 183-195 . 
百度学术
2. 蔡倩茹，王金铎，张关龙，宋智华，王圣柱，熊峥嵘，倪胜利. 准噶尔盆地东北缘上石炭统巴山组烃源岩沉积环境分析及物源示踪——来自泥岩地球化学的证据. 石油实验地质. 2024(01): 146-157 . 百度学术
3. 苑广尧，李凤杰. 柴达木盆地北缘古近系稀土元素地球化学特征及其地质意义. 天然气地球科学. 2023(08): 1374-1384 . 百度学术
4. 常海亮，张宏伟，杜春彦，祝朝辉，段超，侯雪刚. 鹿邑凹陷太原组-山西组泥岩元素地球化学特征及对沉积环境的指示. 科学技术与工程. 2022(25): 10845-10856 . 百度学术
5. 马万里，江小青，李璇，曾亮，杨平，马锦龙. 柴达木盆地西北缘上干柴沟组泥岩地球化学特征与古环境古气候意义. 矿物岩石地球化学通报. 2021(05): 1166-1180 . 百度学术
其他类型引用(2)
下载:
计量
- 文章访问数: 52
- HTML全文浏览量: 8
- PDF下载量: 32
- 被引次数: 7
