High resolution climatic and environmental changes in Shandong Peninsula during the last glacial maximum
-
摘要:
末次冰盛期是末次冰期以来地球表面最寒冷的时期。通过对山东半岛地下画廊溶洞石笋DXHL3进行高精度230Th定年和高分辨率氧-碳稳定同位素分析，重建了该地区23.2～18.8 ka时期高分辨率夏季风气候环境变化历史。结果发现，山东半岛末次冰盛期夏季风气候存在显著的千年尺度波动；REDFIT分析表明，氧、碳同位素值均存在显著的2.2 ka周期。末次冰盛期山东半岛总体处于冷干状态，具体又可细分为4个阶段，即冷干（23.2～22 ka）-相对暖湿（22～21 ka）-冷干（21～19 ka）-相对暖湿（19～18.8 ka）。石笋DXHL3的δ18O变化趋势总体上与亚洲季风区内其他石笋记录相似，也与西太平洋暖池婆罗洲的石笋记录类似。这反映了末次冰盛期西太平洋暖池气候环境变化可能对山东半岛的夏季风气候变化产生了重要影响，推测西太平洋暖池表层海水温度变化引起的大气环流和黑潮的变化，是将气候环境变化信号从低纬地区传递到中高纬地区的纽带。
Abstract:The last glacial maximum is the most cold-dry period on the Earth's surface since the last glaciation. We reconstructed high-resolution summer monsoon changes over the Shandong Peninsula during the period of 23.2～18.8 ka using precisely-dated high-resolution oxygen and carbon isotopes (δ18O and δ13C). Results show that the summer monsoon climate in the last glacial maximum in Shandong Peninsula had significant millennial scale fluctuations, and REDFIT analysis shows a significant 2.2 ka cycle in both oxygen and carbon isotopes, which was in a cold-dry state in general, and can be subdivided into four stages, namely, cold and dry (23.2～22 ka), relatively warm and wet (22～21 k), cold and dry (21～19 ka), and relatively warm and wet (21～19 ka). The newly obtained δ18O record is generally similar to other records of stalagmite δ18O from monsoonal China, as well as the stalagmite δ18O records from northern Borneo in the Western Pacific Warm Pool (WPWP). The climate and environment changes in the WPWP might have exerted a substantial impact on summer monsoon climate over the Shandong Peninsula during the last glacial maximum. It is speculated that temperature change at sea surface of WPWP affected the atmospheric circulation and Kuroshio, which should be the main mechanism sending climate and environment signals from low latitude to middle-high latitude zones in the northern hemisphere.
-
热带海气相互作用在不同时间尺度上影响了全球气候变化。海气相互耦合作用的变化，首先会影响海洋生物圈，尤其是大洋上层的浮游植物。在低纬海区，浮游植物的活跃程度可以通过初级生产力来进行量化，并与海洋表层的风应力动态相关[1]。自Beaufort等成功建立热带大洋初级生产力转换方程，并开创性地应用到赤道太平洋5个柱状样岩芯中以来[2, 3]，利用颗石藻特征属种Florisphaera profunda相对百分含量恢复古生产力，在热带太平洋及其边缘海得到了广泛的应用[4-18]。到目前为止，热带太平洋的古生产力研究大多显示了冰期高、间冰期低的变化特征[4, 6, 8, 10]，众多证据表明该变动受到北半球高纬地区冰量变动的远程调控[4, 5, 8]。
然而，由于海气相互作用的复杂性，对于热带西太平洋海区古生产力在冰期—间冰期时间尺度上变化的驱动机制，到目前为止依然存在争议。通过对热带大洋多个岩芯古生产力的对比以及频谱分析，Beaufort等提出晚更新世古生产力的变化受到冰期旋回以及长期ENSO (El Niño-Southern Oscillation)过程的影响[3]；这一观点在随后得到海水表层温度和其他古生产力等古海洋学记录证据的支持[8, 9, 13, 19]。然而，也有学者持不同的看法，例如，Zhang等提出西赤道太平洋生产力变动的主控因素是来自东亚大陆中部的风尘输入[20]。随着研究的深入开展，发现热带西太平洋初级生产力的受控因素较为复杂，可能是长期ENSO过程以及东亚季风共同调控的结果。但是由于研究指标的局限性，往往侧重某一方面进行讨论。Li等认为尽管东亚冬季风携带来的营养物质能够在一定程度上影响生产力，但冰期—间冰期水体结构的变动对生产力的作用不容忽视，并提出冰期热带太平洋类似于现代El Niño状态，热带西太平洋温跃层上升，因此加强了下部营养物质的向上输送，有利于初级生产力的增加[9]。Xu等则利用热带西太平洋70万年以来东亚冬季风风尘记录及溶解铁与生产力指标的对比，进一步支持了东亚风尘带来的营养物质对冰期生产力的刺激作用[21]。
热带西菲律宾海位于西太平洋暖池北部边缘，由于同时受东亚季风和热带ENSO过程的影响，因此对气候变化的响应非常敏感。东亚冬季风在1—3月爆发，可以引起颗石藻的繁盛[22]；同时，ENSO过程也可以造成热带西太平洋水体结构的变动，从而影响生产力的变化[3, 23]。近年来热带西太平洋古海洋学的发展，提供了众多指示东亚季风强弱以及长期ENSO过程的记录，为进一步开展西菲律宾海冰期—间冰期时间尺度的古生产力研究提供了条件。在西菲律宾海选取合适的岩芯，通过与前人已有的研究工作的对比，有利于进一步查明热带西太平洋古生产力变化的驱动机制。
1. 材料和方法
1.1 研究材料
本研究使用的MD06-3047柱状样岩芯(17°00.44′N、124°47.93′E，水深2510m，柱长8.9m)，于2006年由法国极地研究所的Marion-Dufresne考察船取自吕宋岛以东海域本哈姆海台(图 1)。本哈姆海台是沿着中央海盆洋脊方向的海底高原，约在45~50Ma期间形成[24]。钻孔位置位于吕宋岛以东约240km。由于吕宋岛东部大陆架很窄(平均少于10km)，且岩芯位置距离陆架相对较远，因此海平面变化导致的来自吕宋岛的陆源输入对于本站位的影响可以忽略[25]。吕宋岛以及邻近海域主要受亚热带东亚季风气候的控制。
![]()
图 1 MD06-3047岩芯站位图(根据ODV软件制图)Figure 1. Map of the sampling location in the west Philippine Sea(from Schilitzer, R., Ocean Data View, odv.awi.de, 2017; available at http://odv.awi.de)岩芯以黄褐色粉砂质黏土为主，未见浊流层以及滑坡沉积。MD06-3047岩芯所在海区现代溶跃面深度约3400m [26, 27]，沉积物样品中出现翼足类壳体，表明钙质微体化石保存情况良好[28]。MD06-3047岩芯年代模式根据底栖有孔虫C. wuellerstorfi的δ18O记录[19]与LR04氧同位素标准曲线[29]进行对比建立。岩芯上部330cm用于本研究工作，按照2cm间隔取样，用于颗石藻特征种属鉴定，时间跨度为260ka至晚全新世，分辨率约为1.5ka。
1.2 Florisphaera profunda——有效的古生产力指标
颗石藻是一种广泛生存在海洋透光带的浮游植物，对光照和营养物质变动响应敏感。相对于透光带上部，下透光带相对光线不足但营养物质较高。在热带大洋，下透光带属种主要以F. profunda为主，而大部分颗石藻属种生活在上部透光带[30]。F. profunda占颗石藻群落的相对百分含量与营养跃层的变动密切相关[31]。营养跃层是指海洋上层水体营养盐浓度(如NO3-, PO43-)发生显著变动的水层。Li等根据ODV数据库的现代调查数据，绘制了西菲律宾海区Ph05-5站位附近温度、盐度、磷酸盐、硝酸盐等参数的深度剖面曲线，认为研究区的年均营养跃层深度约在200m水深[7]。当营养跃层变浅，颗石藻种群以上部透光带属种为主，而当营养跃层加深的时候，颗石藻种群中下透光带属种F. profunda的比例则大大增加。在低纬度开放大洋，营养跃层的深度主要受到风力强度的影响。风力强时，上部水层混合加强，下部营养物质上涌至上部透光带；反之，当风力减弱的时候，上部水体混合作用减弱，透光带营养物质输入减少[5]。营养跃层深度变化与颗石藻种群组成的这种关系，已成为追踪热带大洋营养跃层深度变化的重要指标[6, 10, 31-33]。同时，通过校正F. profunda相对百分含量(Fp%)与初级生产力(PP)的关系，建立了有效的古生产力替代性指标[2]。F. profunda相对丰度与初级生产力的函数关系如下：
$$ {\rm{PP}} = 617 - [279{\rm{log}}({\rm{Fp}} + 3)] $$ (1) 其中，PP是年均生产力(gC ·m-2·a-1)，Fp是F. profunda的相对百分含量×100。该方程基于印度洋低纬海区大量表层沉积物样品建立。研究表明，利用Fp得到的古生产力变化与其他古生产力替代性指标相一致[2]。
1.3 实验方法
取沉积物样品1~5mg放至载玻片上，用洗瓶滴加一滴蒸馏水使其分散，用圆滑的牙签侧面来回涂抹至均匀，自然晾干后，在中间滴适量中性树胶，用盖玻片进行固定。将制备的玻片在烘箱中50℃烘干。颗石鉴定工作在同济大学海洋地质国家重点实验室，利用颗石自动鉴定系统对玻片进行常见种属统计。统计方法如下：在Leica偏光显微镜(LEICA DM6000 B)下放大500倍，随机选取40个颗石分布均匀清楚的视域，利用数字摄像头(Spot Inlight FireWire，200万像素)进行拍照。图片进行适当处理后，运行Syraco4.0软件分析统计颗石数目和种属鉴定。该自动鉴定的方法已经成功地应用于不同海区多个岩芯中[3, 10, 11, 32, 34]。通过计算获得颗石藻下透光带属种F. profunda占全部颗石个数的百分含量，并按照方程(1)计算获得260ka以来初级生产力变化。
2. 结果
2.1 F. profunda相对丰度变化
260ka以来，F. profunda百分含量为45%~93%，平均值为67%, 其在MIS 8，MIS 6，MIS 4以及MIS 2的平均值分别为58.1%、66.6%、62.4%、66.6%，而在MIS 7，MIS 5，MIS 3以及全新世的平均值分别为63%、72.5%、68.9%、69.6%。尽管F. profunda百分含量在间冰期平均值相对冰期平均值略高，但通过与代表冰期旋回变化的底栖有孔虫δ18O曲线进行对比(图 2)，发现两者之间相关性并不好，因此认为F.profunda%在过去26万年来冰期-间冰期旋回特征不明显。相对于冰期旋回时间尺度来说，Fp%在千年尺度上变化特征更为明显。比如，Fp%在MIS 5a, MIS 5b, MIS 5c出现低值，分别为59.4%、57.6%、51.9%，而在MIS 5d, MIS 5e则呈现高值，分别为68.4%、80.5%，变化幅度高达28%，指示了在MIS 5期内，营养跃层发生了重大的变动，其在MIS 5d之前相对较深，且波动较少，而在MIS 5c—MIS 5a期间营养跃层变化幅度很大，并出现营养跃层异常浅的现象。在冰期，MIS 4—MIS 2期以及MIS 6期，Fp%的变化幅度甚至更大，最高可达40%，指示了研究区营养跃层在千年尺度和亚轨道尺度上的振荡。此外，F.profunda相对百分含量还具有长期变化的趋势，整体可以划分为两个阶段，在150~260ka期间，即MIS 8后期至MIS 6中期，整体相对较低，而在150ka至晚全新世，整体有所升高，表明西菲律宾海营养跃层自MIS 8期以来呈现了由浅变深的整体变化趋势。
![]()
图 2 260 ka以来MD06-3047岩芯F. profunda相对百分含量以及初级生产力变化曲线MD06-3047岩芯的年代模式由底栖有孔虫C. wuellerstorfi的氧同位素曲线对比LR04标准曲线而建立[19]Figure 2. Records of variations in the abundance of F. profunda and primary productivity for the Core MD06-3047 in the last 260 kaAge model of MD06-3047 was established based on benthic foraminiferal C. wuellerstorfi δ18O stratigraphy by correlating to LR04 δ18O stack [19]利用PAST软件对26万年以来F.profunda%进行单频谱分析，结果显示了较强的轨道周期(55, 28, 19, 14ka)和千年周期(7.5, 4, 3ka)。其中55ka可能是10ka偏心率周期与41ka斜率周期叠加造成，28ka则可能是由斜率周期(41ka)与岁差周期(19ka)合成。
2.2 初级生产力变化
MD06-3047岩芯在研究时间段260ka以来，生产力的变动范围为64~148gC·m-2·a-1，平均值为104gC·m-2·a-1; 初级生产力的变动趋势与Fp%的变化趋势正好相反(图 2)。其在冰期，即MIS 8，MIS 6，MIS 4以及MIS 2的平均值相对较高，分别为120、104.7、111.4、104.4gC·m-2·a-1, 而在间冰期，即MIS 7，MIS 5, MIS 3以及全新世相对较低，平均值分别为110.2、94.1、100、98.9gC·m-2·a-1；此外，初级生产力在MIS 5a, MIS 5b以及MIS 5c出现峰值，分别为116.1、119.7、131.8gC·m-2·a-1，而在MIS 5d和MIS 5e则出现低值，分别为99.9、80.9gC·m-2·a-1。此外，初级生产力表现出长期的变化趋势，即在MIS 8期至155ka期间，初级生产力整体较高，平均值为96.8gC·m-2·a-1，自155ka开始，初级生产力出现快速降低，随后短暂升高，并再次降低，自MIS 5d之后，生产力有所回升，并出现多次峰值，但是155ka至全新世晚期，生产力整体相对较低，平均值为116.4gC·m-2·a-1。
3. 讨论
3.1 相邻海区不同岩芯初级生产力对比
将MD06-3047岩芯获取的初级生产力结果与邻近的站位进行对比，发现不同站位的古生产力记录虽然在整体上存在冰期高—间冰期低的变化趋势，但是各站位的生产力变化范围以及变化幅度之间差别较大。整体看来位于苏禄海区的MD97-2141站位的初级生产力相对较高，变化范围为80~224gC·m-2·a-1, 在过去200ka以来古生产力平均值为136.8gC·m-2·a-1[4]。MD06-3050岩芯次之，在研究时段260ka以来生产力变化范围为69.6~225.5gC·m-2·a-1，平均值120gC·m-2·a-1[32]。MD06-3047岩芯初级生产力变化幅度以及平均值略低于MD06-3050岩芯。Ph05-5岩芯的初级生产力最低，变化幅度为63~112.5gC·m-2·a-1，平均值为83gC·m-2·a-1[35]。说明尽管研究区域相近，不同站位的初级生产力受控因素仍存在较大的差别。由于颗石藻是钙质浮游生物，其在海底的保存受到溶解作用的影响，颗石藻的种属会随着深部溶解作用的加强而有所变化，诸如Gephyrocapsa spp, Emiliania huxleyi等易溶种属相对含量会有所减少，而F.profunda, Calcidiscus leptoporus等抗溶种相对含量则有所增加[36]。Ph05-5岩芯所在站位的现代水深为3382m，接近西太平洋海区现代溶跃面深度(约3400m)[26, 27]，因此，该站位受到溶解作用较强，导致颗石藻易溶种属减少，F.profunda相对百分含量较高，从而导致初级生产力的评估可能受到影响，相对其他站位显示了低值。
此外造成不同站位岩芯初级生产力差别的另外一个重要原因是鉴定方法的差异。Ph05-5站位是人工镜下鉴定，而MD06-3047岩芯和MD06-3050岩芯利用颗石自动鉴定系统进行种属统计，由于自动鉴定系统无法像人工鉴定那样进行视域旋转，可能会漏检到一部分F.profunda，导致人工鉴定的结果得到的古生产力数值往往低于自动鉴定结果[10, 32]，因此，在对不同站位数据结果进行对比时，要注意不同方法所产生的影响。
值得一提的是，尽管MD97-2141岩芯所在水深为3633m，比位于本哈姆海台的MD06-3047，MD06-3050以及Ph05-5岩芯均要深，受溶解作用应该最强，且采用人工鉴定的方法，但是其初级生产力平均值以及变化幅度最高，因此推测在过去的200ka，苏禄海区表层初级生产力整体高于西菲律宾海，这可能与苏禄海区受较强东亚季风有关[5]。
3.2 西菲律宾海区初级生产力变化的受控因素
已有的研究表明热带西太平洋气候条件受东亚季风和ENSO过程的共同影响[7, 8]。更新世以来东亚季风强度以及ENSO过程的发生频率也有了较为深入的研究[21, 25, 37, 38]，然而，到目前为止，在不同时期东亚季风与ENSO过程对热带西太平洋生产力的影响，孰重孰轻？如何评价？仍旧是热带西太平洋海区古气候研究的难点。借助于最近发表的热带东西太平洋温度梯度的记录[28]以及东亚冬季风强度变动指标[21]，可以尝试探讨热带西菲律宾海古生产力变化的影响因素。
Beaufort等通过对赤道太平洋和印度洋9个高质量岩芯的颗石藻初级生产力记录进行频谱分析，发现了显著的30ka周期和岁差周期，提出类ENSO过程是驱动晚更新世低纬大洋初级生产力变化的重要机制[3]。对MD06-3047岩芯F.profunda%进行频谱分析，同样发现了比较显著的30ka周期和岁差周期(图 3)，表明该岩芯受到热带ENSO过程的影响。在现代热带太平洋，由于信风的作用，形成了自东向西的北赤道流，大量暖水遇到菲律宾岛屿等地形阻隔后积聚在热带西太平洋，因此，热带西太平洋具有相对较厚的温跃层/营养跃层，相对较高的海水表层温度，并与东赤道太平洋温跃层以及表层海水温度形成明显的梯度。当发生El Niño时，赤道信风显著减弱，北赤道流减弱，西太平洋暖池的部分温暖海水向东输送，因此，西太平洋温跃层变浅，东太平洋温跃层加深，且表层海水温度有所升高，ΔSST则相对减少；反之，在La Niña年间，东西太平洋ΔSST有所增加。古海洋学研究表明，热带太平洋存在长期的类ENSO过程[28, 39-43]，当东西太平洋ΔSST较大时，认为热带太平洋处于长期类La Niña状态，反之则为类El Niño状态。
![]()
图 3 MD06-3047站F. profunda百分含量与初级生产力频谱分析(图中标出80%和90%置信曲线)Figure 3. The spectral analysis for F. profunda percentageX-axis indicates frequency, Y-Axis stands for power从图 4中可以看出在MIS 8后期至MIS 5a, 即260~80ka，热带东西太平洋温度差值(ΔSST)与Fp%和PP表现出一致性，当ΔSST较小时，F. profunda百分含量相对较低，指示营养跃层较浅，对应于较高的初级生产力，热带西太平洋处于类El Niño的状态；当ΔSST较大时，F. profunda百分含量相对较高，指示营养跃层加深，恰好对应于较低的初级生产力，此时热带西太平洋对应于类La Niña状态。而在该时段，相同岩芯中能够指示东亚冬季风强弱的伊利石/蒙脱石指标，与PP以及Fp%曲线相似性较低，因此，认为在MIS 5a之前，长期的ENSO过程相对于东亚冬季风来说，对热带西太平洋古生产力变化的影响可能更为显著。
![]()
图 4 西太平洋边缘海区260ka以来不同站位F.profunda百分含量对比MD97-2141(8°47′N、121°17′E, 水深3633m，苏禄海区)资料据文献[4]；Ph05-5(16°3′N、124°21′E, 水深3382m, 本哈姆高原)资料据文献[35]；MD06-3050(15°57.0943′N, 124°46.7747′E; 水深2967m，本哈姆高原)资料据文献[5, 33, 35]Figure 4. Comparison of F.profunda% record at different sites in Western Pacific marginal seas in the last 260ka然而，自MIS 5末期，东西太平洋ΔSST的变化与Fp%呈现相反的趋势，即ΔSST相对较低的时候，如MIS 3期，Fp%值相对较高，指示了营养跃层较深，初级生产力相对较低，这与类El Niño状态时较浅的温跃层/营养跃层相矛盾，同理，当ΔSST较大时，指示了热带西太平洋为长期的类La Niña状态，如MIS 2期，此时，热带西太平洋应呈现较深的温跃层/营养跃层，然而实际上Fp%显示了低值，指示了较浅的营养跃层，对应着较高的初级生产力。基于以上讨论，推测在MIS 5a至晚全新世期间，长期ENSO过程并非影响生产力变化的主要因素。实际上，基于单个有孔虫个体的Mg/Ca比的海水表层温度指示了在末次冰期冰盛期(LGM)热带太平洋ENSO活动有所减弱，进一步支持了我们的推论[41]。
相反，在MIS 5a至末次冰消期，伊利石/蒙脱石比值与初级生产力之间表现了较高的相似性变化(图 5)。伊利石/蒙脱石的高值，指示了东亚冬季风的加强。中尺度铁肥实验表明，对包括热带太平洋在内的很多大洋海区，海洋上层常因缺铁等营养因素而使浮游植物的生长受到限制[44]。有证据表明，在赤道太平洋中部，冰期大陆矿物悬浮颗粒(风尘)的通量最大可以超过间冰期该通量的2倍[45]。冰期，风尘的增加，会带来更多的铁，很可能会刺激热带赤道太平洋的生物生产力。Zhang等的研究支持了这一观点，认为过去550ka以来，西赤道太平洋生物生产力增加的主要驱动者很可能是来自东亚中部大陆的风尘通量[20]。此外，Kim等对西太平洋沉积物捕获器样品的分析表明，在La Niña强盛的时期，颗粒物通量的增加与季风加强造成的上部水体的混合加强有关[23]。东亚冬季风的加强一方面使得研究区上部水体混合增加，营养跃层变浅，向上输送的营养物质增加；另一方面加强了风尘的输入，刺激了上部水体浮游植物的生长。尽管不能定量地评估风尘通量对生产力的影响，但可以肯定的是，东亚冬季风加强以及大陆的干旱条件，在冰期向海洋输送了大量的大陆风尘，进而影响热带西太平洋生产力的增加[21]。苏禄海MD97-2141岩芯也显示了生产力在冰期增加，认为与东亚冬季风的加强有关，并且提出在过去的8万年间，苏禄海区存在8个显著的初级生产力增加事件，对应于该时期中国黄土记录的东亚季风的加强[4]。而当伊利石/蒙脱石出现低值时，即MIS3期，东亚冬季风相对较弱，上部水体混合减弱，风尘输入也有所减少，对应于上部水体相对较低的初级生产力。
![]()
图 5 过去260ka以来西菲律宾海MD06-3047岩芯初级生产力及F. profunda相对百分含量变化与其他古环境指标对比图自上而下分别为东亚冬季风替代指标——黏土矿物伊利石/蒙脱石变化曲线[21]；基于经验公式计算得到的初级生产力变化曲线(本研究)；指示热带太平洋长期ENSO过程的热带东西太平洋海水表层温度梯度异常变化曲线(相对于全新世SST数值)[28]；颗石藻下透光带属种F. profunda相对百分含量变化曲线(本研究)Figure 5. Correlation between primary productivity and abundance of F. profunda of Core MD06-3047 for the past 260ka, compared with long term ENSO proxy and East Asian Winter Monsoon proxyThe above curve is illite/smectite ratio in clay mineral fraction of core MD06-3047 that could indicate the intensity of the East Asian Winter Monsoon [21]; Zonal SST gradient anomaly relative to the Holocene values, which represents the SST difference of the western and eastern Pacific based on the Mg/Ca-based SSTs [28]东亚季风与ENSO过程之间存在着复杂的耦合过程[46]，一定程度上解释了MD06-3047岩芯呈现较不明显的冰期-间冰期旋回变化以及较强的千年尺度变化。通过多指标对比研究表明，在260~80ka，西菲律宾海区初级生产力受到长期ENSO过程下营养跃层/温跃层变动的影响相对较强；而在80 ka以来，东亚冬季风相对于长期ENSO过程来说，对该区初级生产力的影响可能更为显著。
4. 结论
对位于西菲律宾海的MD06-3047岩芯，进行颗石藻化石鉴定统计，获取F. profunda百分含量，并基于此恢复了该海区26万年以来的初级生产力和营养跃层的变化历史，发现初级生产力冰期旋回的变动特征并不明显。频谱分析结果显示西菲律宾海上层水体结构具有明显的轨道周期和千年尺度变化周期，说明该区的初级生产力的受控因素相对较为复杂。
(1) 26万年到8万年之间，东西太平洋温度梯度(ΔSST)与初级生产力变动具有较强的相似性，指示了长期ENSO过程对西菲律宾海初级生产力的影响相对于东亚季风对初级生产力的影响更为显著，认为在MIS7初期，MIS6期前半段以及MIS5中后期，热带太平洋呈现类El Niño状态，对应研究区较浅的营养跃层/温跃层，初级生产力相对较高；其余时段热带太平洋则处于长期的类La Niña状态，对应较深的营养跃层/温跃层，初级生产力相对较低。
(2) 8万年以来，西菲律宾海初级生产力受东亚冬季风的影响加强或者ENSO状态的减弱，可能掩盖了长期ENSO过程对生产力影响的信号。冰期时东亚冬季风加强，一方面可以加强上部水体的混合作用，另一方面增加了来自东亚中部的风尘物质的输入，两者均有利于初级生产力的增加；而在间冰期，东亚冬季风减弱，上部水体混合作用减少，同时风尘物质输入量也显著降低，因此初级生产力相对降低。
-
![]()
图 2 石笋DXHL3与年龄模式
a: 石笋剖面图，其中黑色短线代表测年点位置； b: 石笋的年龄-深度模式图；c为a中进行显微镜观察部分（黑框内）。
Figure 2. Staglagmite DXHL3 and its age model
a: Profile of stalagmite, where short black line indicates the location of the dating point; b: Age-depth plot of stalagmite; c: the microscopic part of (a) (the black box).
![]()
图 3 石笋DXHL3的δ18O和δ13C记录及与极地冰芯、婆罗洲北部石笋及中国季风区内其他石笋δ18O 记录对比
a：格陵兰冰芯NGRIP[17]； b-c：地下画廊溶洞石笋DXHL3； d：婆罗洲北部石笋SSC01[26]；e：豪猪洞石笋HZZ11[18]；f：葫芦洞石笋MSD[19]； g：南极冰芯EDML[17]。各石笋记录均显示了年龄及误差；图中黑色圆点和水平误差棒代表相应记录的定年点及误差；黑色箭头代表变化趋势；所有数据均进行了平滑滤波处理，灰色虚线代表的是原数据，黑色实线为三点滑动平均后数据。
Figure 3. Comparison of the DXHL3 δ18O and δ13C records with δ18O records of ice cores from polar areas, speleothems from northern Borneo, and other stalagmites in monsoonal China
a: Greenland ice core NGRIP[17]; b-c: underground gallery cave stalagmite DXHL3; d: North Borneo stalagmite SSC01[26]; e: Haozhu cave stalagmite HZZ11[18]; f: Hulu cave stalagmite MSD[19]; g: Antarctic ice core EDML[17]. All stalagmite records show age and error. Black dots and horizontal error bars represent the corresponding recorded dating points and errors. Black arrows indicate the changing trend. All data are smoothed and filtered; the gray dashed line are the original data and the solid black line are the three-point sliding average.
![]()
图 4 石笋DXHL3的氧（a）和碳（b）同位素REDFITF分析结果
MC表示Monte Carlo置信水平线。
Figure 4. REDFIT analysis of the oxygen (a) and carbon (b) isotopes of stalagmite DXHL3
MC: Monte Carlo confidence level.
表 1 石笋DXHL3的U-230Th年代数据
Table 1 U-230Th dating results of stalagmite DXHL3
样品号 深度/mm 238U/10−9 232Th/10−12 实测 δ234U [230Th/238U]活度比 [230Th/232Th]/10−6 未校正年龄/aBP 校正年龄/aBP δ234Ui DXHL-1 3.5 24467±21 2220±8 3015.5±4.5 0.67890±0.00072 123563±448 19774±33 19774±33 3188.8±4.7 DXHL-3 7.0 22657±28 635±22 3147.5±5.5 0.73498±0.00097 432661±14833 20794±42 20794±42 3338.0±5.8 DXHL-4 10.0 19816±26 56±7 3189.6±6.4 0.77094±0.00106 4531392±550191 21654±48 21654±48 3390.9±6.8 DXHL-5 26.5 18180±25 165±16 3120.4±6.8 0.77806±0.00114 1413001±138238 22270±53 22270±53 3323.2±7.2 DXHL-6 31.0 29792±43 128±21 3118.2±8.0 0.78114±0.00124 3004717± 499744 22379±61 22379±61 3321.8±8.6 DXHL-7 59.0 27930±33 323±15 3115.6±6.5 0.78600±0.00099 1120882±50279 22545±49 22545±49 3320.6±6.9 DXHL-8 64.2 24744±32 621±23 3051.2±6.2 0.77506±0.00111 509515±18740 22591±51 22591±51 3252.4±6.6 DXHL-9 90.0 23807±23 66±20 3054.3±4.6 0.77846±0.00116 4609934±1361855 22680±46 22680±46 3256.5±5.0 DXHL-10 111.0 35390±39 293±18 3186.8±5.5 0.80779±0.00136 1609388±97936 22793±53 22793±53 3398.8±5.9 DXHL-2 132.8 31713±32 50±9 3113.7±5.1 0.80362±0.00091 8416357±1550677 23106±42 23106±42 3323.8±5.4 注: λ230=9.1599×10−6 a−1 ；λ234=2.8263×10−6 a−1 ；λ238=1.55125×10−10 a−1 。δ234U=([234U/238U]活度比−1)×1000；234U初始值是根据230Th年龄获得的，即δ234U初始值=δ234U测量值×eλ234×t。校正230Th年龄假设初始的230Th/232Th摩尔比为(4±2)×10−6。 
下载: 导出CSV
-
[1] Clark P U, Dyke A S, Shakun J D, et al. The last glacial Maximum [J]. Science, 2009, 325(5941): 710-714. doi: 10.1126/science.1172873
[2] Yokoyama Y, Lambeck K, De Deckker P, et al. Timing of the Last Glacial Maximum from observed sea-level minima [J]. Nature, 2000, 406(6797): 713-716. doi: 10.1038/35021035
[3] 杨石岭, 董欣欣, 肖举乐. 末次冰盛期以来东亚季风变化历史: 中国北方的地质记录[J]. 中国科学:地球科学, 2019, 62(8):1181-1192 doi: 10.1007/s11430-018-9254-8 YANG Shiling, DONG Xinxin, XIAO Jule. The East Asian Monsoon since the Last Glacial Maximum: Evidence from geological records in northern China [J]. Science China:Earth Sciences, 2019, 62(8): 1181-1192. doi: 10.1007/s11430-018-9254-8
[4] 侯光良, 许长军, 兰措卓玛, 等. 末次冰盛期中国人类活动的响应与适应[J]. 热带地理, 2018, 38(6):819-827 HOU Guangliang, XU Changjun, LANCUO Zhuoma, et al. The response and adaptation of Chinese human activities to the Last Glacial Maximum [J]. Tropical Geography, 2018, 38(6): 819-827.
[5] 王照波, 张剑, 王江月, 等. 山东蒙山第四纪冰川遗迹光释光测年研究及冰期划分与对比[J]. 地质论评, 2017, 63(3):694-702 WANG Zhaobo, ZHANG Jian, WANG Jiangyue, et al. A study on quaternary glaciation moraines in mount Mengshan, Shandong Province with optically stimulated luminescence (OSL) and comparison of glacial period division [J]. Geological Review, 2017, 63(3): 694-702.
[6] Yang S L, Ding Z L, Li Y Y, et al. Warming-induced northwestward migration of the East Asian monsoon rain belt from the Last Glacial Maximum to the mid-Holocene [J]. Proceedings of the National Academy of Sciences of the United States of America, 2015, 112(43): 13178-13183. doi: 10.1073/pnas.1504688112
[7] 李曼玥, 张生瑞, 许清海, 等. 华北平原末次冰盛期以来典型时段古环境格局[J]. 中国科学:地球科学, 2019, 62(8):1279-1287 doi: 10.1007/s11430-018-9264-2 LI Manyue, ZHANG Shengrui, XU Qinghai, et al. Spatial patterns of vegetation and climate in the North China Plain during the Last Glacial Maximum and Holocene climatic optimum [J]. Science China:Earth Sciences, 2019, 62(8): 1279-1287. doi: 10.1007/s11430-018-9264-2
[8] 杨利荣, 岳乐平, 弓虎军. 呼伦贝尔沙地末次冰盛期晚期至全新世风成沙表面矿物特征及环境意义[J]. 地理研究, 2015, 34(6):1066-1076 YANG Lirong, YUE Leping, GONG Hujun. The environmental implication from microscopic texture of eolian sand of Hulun Buir duneland centred on late last glacial maximum and Holocene [J]. Geographical Research, 2015, 34(6): 1066-1076.
[9] Chen Q M, Cheng X, Cai Y J, et al. Asian summer monsoon changes inferred from a stalagmite δ18O record in central China during the last glacial period [J]. Frontiers in Earth Science, 2022, 10: 863829. doi: 10.3389/feart.2022.863829
[10] Cheng H, Edwards R L, Sinha A, et al. The Asian monsoon over the past 640, 000 years and ice age terminations [J]. Nature, 2016, 534(7609): 640-646. doi: 10.1038/nature18591
[11] 程海, 张海伟, 赵景耀, 等. 中国石笋古气候研究的回顾与展望[J]. 中国科学:地球科学, 2019, 62(10):1489-1513 doi: 10.1007/s11430-019-9478-3 CHENG Hai, ZHANG Haiwei, ZHAO Jingyao, et al. Chinese stalagmite paleoclimate researches: a review and perspective [J]. Science China:Earth Sciences, 2019, 62(10): 1489-1513. doi: 10.1007/s11430-019-9478-3
[12] 耿庆明, 丁长慎, 耿升明. 山东沂水县水文地质条件与分区探讨[J]. 山东师大学报:自然科学版, 1990, 5(2):66-70 GENG Qingming, DING Changshen, GENG Shengming. A study of hydrogeologic conditions and zones in Yishui county of Shandong Province [J]. Journal of Shandong Normal University:Natural Science Edition, 1990, 5(2): 66-70.
[13] 迟宏. 近千年以来山东半岛高分辨率气候环境变化的石笋记录[D]. 鲁东大学硕士学位论文, 2016. CHI Hong. Climate change in the past 1000 years revealed by stalagmite from Kaiyuan Cave in Shandong Peninsula, China[D]. Master Dissertation of Ludong University, 2016.
[14] Shen C C, Wu C C, Cheng H, et al. High-precision and high-resolution carbonate 230Th dating by MC-ICP-MS with SEM protocols [J]. Geochimica et Cosmochimica Acta, 2012, 99: 71-86. doi: 10.1016/j.gca.2012.09.018
[15] Shen C C, Cheng H, Edwards R L, et al. Measurement of attogram quantities of 231Pa in dissolved and particulate fractions of seawater by isotope dilution thermal ionization mass spectroscopy [J]. Analytical Chemistry, 2003, 75(5): 1075-1079. doi: 10.1021/ac026247r
[16] Schulz M, Mudelsee M. REDFIT: estimating red-noise spectra directly from unevenly spaced paleoclimatic time series [J]. Computers & Geosciences, 2002, 28(3): 421-426.
[17] Lemieux-Dudon B, Blayo E, Petit J R, et al. Consistent dating for Antarctic and Greenland ice cores [J]. Quaternary Science Reviews, 2010, 29(1-2): 8-20. doi: 10.1016/j.quascirev.2009.11.010
[18] Zhang H B, Griffiths M L, Huang J H, et al. Antarctic link with East Asian summer monsoon variability during the Heinrich Stadial–Bølling interstadial transition [J]. Earth and Planetary Science Letters, 2016, 453: 243-251. doi: 10.1016/j.jpgl.2016.08.008
[19] Wang Y J, Cheng H, Edwards R L, et al. A high-resolution absolute-dated late Pleistocene monsoon record from Hulu Cave, China [J]. Science, 2001, 294(5550): 2345-2348. doi: 10.1126/science.1064618
[20] 王庆, 周厚云, 迟宏, 等. 最近千年来山东半岛西部气候环境变化的石笋δ18O、δ13C记录(I)[J]. 海洋地质与第四纪地质, 2015, 35(5):135-142 WANG Qing, ZHOU Houyun, CHI Hong, et al. The stalagmite records of climate and environment change on the western Shandong Peninsula during the past 1000 years: 18O and 13C values (I) [J]. Marine Geology & Quaternary Geology, 2015, 35(5): 135-142.
[21] Tan L C, Liu W, Wang T L, et al. A multiple-proxy stalagmite record reveals historical deforestation in central Shandong, northern China [J]. Science China Earth Sciences, 2020, 63(10): 1622-1632. doi: 10.1007/s11430-019-9649-1
[22] Xue G, Cai Y J, Ma L, et al. A new speleothem record of the penultimate deglacial: Insights into spatial variability and centennial-scale instabilities of East Asian monsoon [J]. Quaternary Science Reviews, 2019, 210: 113-124. doi: 10.1016/j.quascirev.2019.02.023
[23] 李倩, 李广雪, 张强, 等. 末次冰消期临沂石笋微量元素记录及其气候意义[J]. 中国海洋大学学报, 2018, 48(11):100-107 LI Qian, LI Guangxue, ZHANG Qiang, et al. Stalagmite trace element records and implications for paleoclimate change in Linyi during the Last Deglaciation [J]. Periodical of Ocean University of China, 2018, 48(11): 100-107.
[24] Li Q, Li G X, Chen M T, et al. East Asian summer monsoon variations during the last deglaciation, recorded from a stalagmite at Linyi, northern China [J]. Quaternary International, 2017, 464: 327-335.
[25] Liu Z Y, Wen X Y, Brady E C, et al. Chinese cave records and the East Asia summer monsoon [J]. Quaternary Science Reviews, 2014, 83: 115-128. doi: 10.1016/j.quascirev.2013.10.021
[26] Partin J W, Cobb K M, Adkins J F, et al. Millennial-scale trends in west Pacific warm pool hydrology since the Last Glacial Maximum [J]. Nature, 2007, 449(7161): 452-455. doi: 10.1038/nature06164
[27] 黄荣辉, 李维京. 夏季热带西太平洋上空的热源异常对东亚上空副热带高压的影响及其物理机制[J]. 大气科学, 1988, 12(S1):107-116 doi: 10.3878/j.issn.1006-9895.1988.t1.08 HUANG Ronghui, LI Weijing. Influence of heat source anomaly over the western tropical Pacific on the subtropical high over East Asia and its physical mechanism [J]. Chinese Journal of Atmospheric Sciences, 1988, 12(S1): 107-116. doi: 10.3878/j.issn.1006-9895.1988.t1.08
[28] Huang R H, Sun F Y. Impacts of the tropical western Pacific on the East Asian Summer Monsoon [J]. Journal of the Meteorological Society of Japan. Ser. II, 1992, 70(1B): 243-256. doi: 10.2151/jmsj1965.70.1B_243
[29] 黄荣辉, 陈际龙, 黄刚, 等. 中国东部夏季降水的准两年周期振荡及其成因[J]. 大气科学, 2006, 30(4):545-560 doi: 10.3878/j.issn.1006-9895.2006.04.01 HUANG Ronghui, CHEN Jilong, HUANG Gang, et al. The quasi-biennial oscillation of summer monsoon Rainfall in China and its cause [J]. Chinese Journal of Atmospheric Sciences, 2006, 30(4): 545-560. doi: 10.3878/j.issn.1006-9895.2006.04.01
[30] Li D, Tan L C, Cai Y J, et al. Is Chinese stalagmite δ18O solely controlled by the Indian summer monsoon? [J]. Climate Dynamics, 2019, 53(5-6): 2969-2983. doi: 10.1007/s00382-019-04671-x
[31] Wu J Y, Wang Y J, Cheng H, et al. An exceptionally strengthened East Asian summer monsoon event between 19.9 and 17.1 ka BP recorded in a Hulu stalagmite [J]. Science in China Series D:Earth Sciences, 2009, 52(3): 360-368. doi: 10.1007/s11430-009-0031-1
[32] Caley T, Roche D M, Renssen H. Orbital Asian summer monsoon dynamics revealed using an isotope-enabled global climate model [J]. Nature Communications, 2014, 5: 5371. doi: 10.1038/ncomms6371
[33] Cai Y J, Tan L C, Cheng H, et al. The variation of summer monsoon precipitation in central China since the last deglaciation [J]. Earth and Planetary Science Letters, 2010, 291(1-4): 21-31. doi: 10.1016/j.jpgl.2009.12.039
[34] Jiang X Y, He Y Q, Shen C C, et al. Stalagmite-inferred Holocene precipitation in northern Guizhou Province, China, and asynchronous termination of the climatic Optimum in the Asian monsoon territory [J]. Chinese Science Bulletin, 2012, 57(7): 795-801. doi: 10.1007/s11434-011-4848-6
[35] Selvaraj K, Chen C T A, Lou J Y, et al. Holocene weak summer East Asian monsoon intervals in Taiwan and plausible mechanisms [J]. Quaternary International, 2011, 229(1-2): 57-66. doi: 10.1016/j.quaint.2010.01.015
[36] Zhao K, Wang Y J, Edwards R L, et al. Contribution of ENSO variability to the East Asian summer monsoon in the late Holocene [J]. Palaeogeography, Palaeoclimatology, Palaeoecology, 2016, 449: 510-519. doi: 10.1016/j.palaeo.2016.02.044
[37] Zhong W, Cao J Y, Xue J B, et al. A 15, 400-year record of climate variation from a subalpine lacustrine sedimentary sequence in the western Nanling Mountains in South China [J]. Quaternary Research, 2015, 84(2): 246-254. doi: 10.1016/j.yqres.2015.06.002
[38] Fairbanks R G. A 17, 000-year glacio-eustatic sea level record: influence of glacial melting rates on the Younger Dryas event and deep-ocean circulation [J]. Nature, 1989, 342(6250): 637-642. doi: 10.1038/342637a0
[39] 业治铮, 汪品先. 南海晚第四纪古海洋学研究[M]. 青岛: 青岛海洋大学出版社, 1992. YE Zhizheng, WANG Pinxian. Research on Late Quaternary Paleoceanography of the South China Sea[M]. Qingdao: Qingdao Ocean University Press, 1992.
[40] Sawada K, Handa N. Variability of the path of the Kuroshio ocean current over the past 25, 000 years [J]. Nature, 1998, 392(6676): 592-595. doi: 10.1038/33391
[41] Jian Z M, Wang P X, Saito Y, et al. Holocene variability of the Kuroshio Current in the Okinawa Trough, northwestern Pacific Ocean [J]. Earth and Planetary Science Letters, 2000, 184(1): 305-319. doi: 10.1016/S0012-821X(00)00321-6
[42] Liu Y G, Fu Y X, Du D W, et al. Paleoceanographic records in the sedimentary cores from the middle Okinawa Trough [J]. Chinese Science Bulletin, 2003, 48(1): 74-81. doi: 10.1360/csb2003-48-1-74
[43] 何炽鹏, 李冬玲, 李正, 等. 末次冰盛期以来冲绳海槽古海洋环境研究进展[J]. 第四纪研究, 2020, 40(6):1577-1587 doi: 10.11928/j.issn.1001-7410.2020.06.17 HE Chipeng, LI Dongling, LI Zheng, et al. Research progress in the paleoenvironment around the Okinawa Trough since the Last Glacial Maximum [J]. Quaternary Sciences, 2020, 40(6): 1577-1587. doi: 10.11928/j.issn.1001-7410.2020.06.17
[44] Xu D K, Lu H Y, Wu N Q, et al. Asynchronous marine-terrestrial signals of the last deglacial warming in East Asia associated with low- and high-latitude climate changes [J]. Proceedings of the National Academy of Sciences of the United States of America, 2013, 110(24): 9657-9662. doi: 10.1073/pnas.1300025110
[45] Li D W, Yu M, Jia Y H, et al. Gradually cooling of the yellow sea warm current driven by tropical pacific subsurface water temperature changes over the past 5 kyr [J]. Geophysical Research Letters, 2021, 48(10): e2021GL093534.
-
期刊类型引用(0)
其他类型引用(1)
计量
- 文章访问数: 206
- HTML全文浏览量: 35
- PDF下载量: 26
- 被引次数: 1
