Changes in siliceous paleoproductivity during the early Middle Pleistocene in the Northwest Bay of Bengal
-
摘要: 西北孟加拉湾位于印度夏季风的核心降水区,是研究印度季风活动及其对环境影响的理想地区。通过对西北孟加拉湾马哈纳迪盆地IODP353航次U1446站位岩芯生物硅含量和沉积通量的分析,重建了研究区早更新世末—中更新世早期(920~370 ka)的硅质生物生产力演变历史,探讨了印度夏季风降水和径流活动对研究区古生产力变化的影响机制。结果显示,西北孟加拉湾在920~370 ka期间,硅质生产力总体上表现为中更新世过渡期(Mid-Pleistocene Transition,MPT)较低,MPT之后呈上升的趋势;具有明显的冰期-间冰期旋回性,冰期时下降,间冰期上升,且在MPT之后波动幅度增大。对比该站位已有的δ18Osw和Rb/Ca比值等夏季风降水和径流指标发现,印度夏季风增强时硅质生产力上升,减弱时硅质生产力下降,推测印度夏季风的增强导致孟加拉湾降水以及周边河流径流量剧增,河流提供的营养物质使得西北孟加拉湾表层硅质生物勃发,硅质生产力迅速升高。而印度夏季风减弱时,营养盐的减少抑制了表层浮游硅质生物的繁殖,从而使得硅质生产力下降。U1446站位生物硅通量的频谱和小波分析也显示,研究区的硅质生产力受轨道周期变化的调控。Abstract: The Northwest Bay of Bengal (BoB) is located in the core precipitation zone of the Indian Summer Monsoon (ISM), and an ideal area to study Indian monsoon activity and its impact on the environment. We reconstructed the history of siliceous productivity evolution in the study area during late Early Pleistocene to early Middle Pleistocene (920~370 ka) by analyzing biogenic silica content and sediment fluxes of cores from Site U1446 of IODP Expedition 353 in the Mahanadi Basin, Northwest BoB, and explored the mechanism of ISM precipitation and runoff activities on the paleoproductivity changes in the study area. Results indicate that siliceous productivity reflects a trend from low in the Mid-Pleistocene Transition (MPT) to high after MPT in the northwest BoB during 920~370 ka in an obvious glacial-interglacial cycle. It decreased during the glacial period and increased during the interglacial period, and the fluctuation increased after the MPT. A comparison of existing summer monsoon precipitation and runoff indicators such as δ18Osw and Rb/Ca ratio at this site reveals that siliceous productivity increased (decreased) when the ISM strengthened (weakened). It is hypothesized that the enhanced ISM led to a dramatic increase in runoff from the Indian peninsula and precipitation in the BoB, and the nutrients from rivers led to a proliferation of surface siliceous organisms and a rapid upsurge in siliceous productivity in the northwest BoB. However, the weakened ISM and the reduction of nutrients inhibited the flourishing of planktonic siliceous organisms, which then decreased the surface productivity. Spectral and wavelet analyses of biogenic silica fluxes at Site U1446 also indicate that siliceous productivity was modulated by orbital cycle variations in this study area.
-
近年来,随着微区分析技术的快速发展,应用激光剥蚀、离子探针等对生物壳体的原位与高分辨分析受到广泛重视。原位微区分析是利用电子束轰击或激光束剥蚀进行采样,实现原位元素—同位素组成分析。自1991年Perkins等[1]使用1 064 nm激光分析贝壳的Mg、Mn、Sr等元素以来,应用激光剥蚀等离子体质谱仪(LA-ICP-MS)分析双壳类贝壳的研究越来越多[2-6]。与传统溶液法相比,激光剥蚀的原位微区分析更高效,可以更准确、连续反映生物体不同生长时期的元素和同位素组成,对建立高分辨率的环境记录具有重要的研究价值。
牡蛎是广温、广盐性的双壳纲软体动物,其钙质外骨骼具有明显的生长纹层,记录了生物体由老到新的生长时间序列,纹层的微量元素和同位素组成与周围水体存在物理化学上的动态平衡[7-9],可以作为反映气候和环境的良好替代性指标[3],如贝壳的Mg/Ca比等常用于指示水体温度[9-15];Sr/Ca比可以指示水体光照程度[16-17]或水体盐度[18-22];Na/Ca比常指示水体盐度[21-24]。传统研究手段主要通过微钻取样、磷酸法提取CO2及分析牡蛎壳的地球化学组成[20-22, 25-28],但受取样分辨率的制约,生物壳体的环境示踪研究还有很大发展空间。
本文选取南黄海海域现代长牡蛎壳体为研究对象,尝试建立生物壳体原位微区高分辨率的化学记录。长牡蛎(Crassostrea gigas)又名太平洋牡蛎,属于珍珠贝目牡蛎科牡蛎属,生长于半咸水环境,广泛分布于中国东南沿海5~40 m水深的潮下带和浅水区[29]。个体寿命可以超过15年,最适宜的盐度范围是20‰~26‰,生长水温是15~25 ℃,高于28 ℃或低于6 ℃时生长缓慢或停止[28-31]。其壳体几乎全由低镁方解石组成[31-32],在韧带部位具有明显的周期性生长层[25],牡蛎壳的年增长层通常由夏季形成的白垩质方解石(Ch,chalky calcite)较厚亮层和冬季形成的叶片状方解石(Fo,foliated calcite)较薄暗层构成[25, 27]。其壳层发育主要受控于环境温度、水深、食物来源等外部条件,是研究中国东部沿海地区河口浅海环境特征和海陆相互作用的理想材料[28, 30-31]。白垩质方解石质地较软,呈亮白色,含有宽孔,通常与春夏季节的快速生长有关[28],易受到成岩、溶解和重结晶作用的影响,通常不适合古气候重建;叶片状方解石由密集的方解石板条组成,质地较硬,呈暗灰色,更能抵抗后期环境改造的影响,其同位素或元素组成可以反映原始的环境条件[11-12]。本文应用LA-ICP-MS分析长牡蛎左壳韧带区域不同生长层的元素组成,探讨其蕴含的环境信息。
1. 材料和方法
1.1 样品来源
南黄海海域牡蛎滩出露于粉砂淤泥质潮滩之上,面积约3.4 km2[30, 33]。据吕四海洋站1968—2003年的实测资料,该水域冬季表层平均水温4.9 ℃,夏季表层平均水温为27.3 ℃,盐度一般为27‰~30‰,海水悬沙浓度一般为0.2~0.7 g/L[34],属中潮海岸,波浪小,年均波高低于0.3 m,大潮平均潮差大于4 m,潮流属往复流,平均流速大于1 m/s[35]。
研究样品为一个已死亡长牡蛎的左壳,于2008年11月采自江苏省海门市近岸的现代南黄海海域牡蛎滩(32°09.265′N、121°32.162′E,图1)。壳体保存完好,全长约20 cm,通体黄白色,鳞片坚厚,背缘有虫蚀痕,韧带槽宽大,韧带部位、闭壳肌痕和生长缘呈同一直线,可直接通过生长层计数确定其生长年限[25, 28, 36]。
1.2 实验方法
通过将牡蛎左壳(图2A)依次以清水冲洗、干燥、沿最大生长线切割、抛光,获得长45 mm、宽10~20 mm、厚5 mm的牡蛎壳韧带部位切片(图2B)。通过肉眼和显微镜下观察,基本以一个白垩质方解石(Ch)层和一个叶片状方解石(Fo)层计作一组生长层[25, 27-28, 36],结合后续的地球化学方法进行调整,共划分14组生长层,单层宽1~4 mm,由老到新依次标记为1—14(图2C)。
原位微区测试分析应用同济大学海洋地质国家重点实验室的LA-MC-ICP-MS完成。激光剥蚀系统为Australian Scientific Instruments公司的Resolution M50L,激光光源为193 nm ArF准分子气体激光器,串接的ICP-MS为Agilent 7900型ICP-MS。外标标样为美国国家标准技术研究院NIST SRM 610人工合成玻璃,内插质控标样为NIST SRM 612玻璃,用以校正ICP-MS灵敏度漂移[37]。
在获取壳体生长年轮信息的基础上,根据壳体剖面内特征生长层的位置,选用线扫描的剥蚀方式,控制线扫描速度为15 μm/s、激光束斑为40 μm、激光能量密度为2.5 J/cm2、激光频率为6 Hz,测线沿壳体生长方向垂直于壳体剖面,自内向外、由新到老,顺叶面方解石层进行展布,进行6条连续短线剥蚀,测线总长34 725.7 μm,轨迹见图2C。线扫描方向与壳生长方向相反。总线扫时间为2 700 s,选取12~2 320 s为有效测试数据范围,并采用测试前20 s和结束后的载气平均响应值作为空白背景。
1.3 数据校正
由ICP-MS直接获取的数据为各质量数的测试响应值(单位为counts per second,即每秒的计数值,记作cps),在扣除背景空白、获取有效的计数率信号cps后,研究采用Longerich的定量方法[38]对cps进行校正,得到元素比值,具体公式如下:
$$ C_{\rm{sam }}^i/C_{\rm{sam }}^{i{\rm s}} = \frac{{{\rm{cps}}_{\rm{sam} }^i}}{{{\rm{cps}}_{\rm{sam} }^{i{\rm s}}}}\Bigg/\left( {\frac{{{\rm{cps}}_{\rm{rm}}^i}}{{{\rm{cps}}_{\rm{rm}}^{i{\rm s}}}}} \right) \times \frac{{C_{\rm{rm}}^i}}{{C_{\rm{rm}}^{i{\rm s}}}} $$ (1) 公式(1)中,上标i表示待测元素,is代表内标元素,下标sam为待测样品,rm为标准样品,C表示元素含量。本文中测试样品sam为牡蛎壳,标准样品rm为NIST 610,内标元素is选用43Ca(牡蛎壳主要成分为CaCO3)。标样推荐值由地质环境标准物质数据库获取(GeoReM,http://georem.mpch-mainz.gwdg.de/),元素比值分析结果以摩尔比(mmol/mol, μmol/mol)形式表示。
2. 结果与讨论
2.1 壳体微量元素组成
LA-ICP-MS线剥蚀模式测试长牡蛎壳体的实测元素18个,将测试元素稳定变化阶段的平均响应值记作cpssam,各元素测试时的空白背景平均信号响应值记作cpsblk,以各元素cpsblk的3σ(3倍标准差)作为样品检出限[37-38]。根据前期研究[39],当cpssam<3σ时,可以认为元素低于LA-ICP-MS检测限;3σblk≤cpssam<10σblk,元素可被检出,但仅能进行定性或半定量分析;cpssam≥10σblk,元素可进行数据定量分析。
本研究中,长牡蛎壳体的LA-ICP-MS线剥蚀信号(图3)存在3种变化模式:
(1)模式1:信号稳定且持续波动,涉及7个元素。Na、Mg、Ca、Sr的信号强度高,cpssam远超10σ,可以进行定量计算。B、P、Mn的信号同样存在持续变化,cpssam为3~10σ,可以定性分析变化趋势。
(2)模式2:剥蚀信号平均响应值cpssam相对较高,曲线上表现为个别层位出现集中单个小峰,其余90%的cpssam略高于空白,涉及6种元素:Al、K、Zn、Cu、Ba、Pb。它们的cpssam为3~10σ,仅能定性分析。
(3)模式3:V、Cr、Rb、Cs、Hg这5种元素含量较低,cpssam低于3σ,无法被LA-ICP-MS检出。
2.2 壳体微量元素与钙比值的定量分析
考虑到测试标样基体不匹配[37, 40]、碳酸盐样品的多离子团原子团干扰[4]的影响,研究中采用微量元素与Ca比值代替其绝对含量进行研究。从比值看,长牡蛎壳体中元素含量分布趋势相近。可以定性分析的元素中,K、P、Al、B、Zn、Mn 6种元素与Ca的比值为0.01~1 mmol/mol,依次为K/Ca>P/Ca>Al/Ca>B/Ca>Zn/Ca>Mn/Ca;Cu、Ba、Pb 3种元素与Ca的比值均低于10 μmol/mol,表现为Cu/Ca>Ba/Ca>Pb/Ca。
对于3种可定量的微量元素,以43Ca为内标元素,对Na/Ca、Mg/Ca、Sr/Ca通过局部加权回归散点平滑(LOESS,locally weighted scatterplot smoothing)进行信号过滤处理,提取低频率的比值信息,其中LOESS拟合比例为0.05。在整个生长周期内,Sr/Ca和Na/ Ca沿生长方向递减,而Mg/Ca反向递增;在生长后期的7~8个生长层中,Na、Mg、Sr的含量均出现了较为规律的周期性波动(图4)。
图 4 牡蛎样品不同生长层中Na/Ca、Mg/Ca、Sr/Ca比值变化细线表示原始元素比值,粗线表示LOESS平滑结果(黄线为Sr/Ca,蓝线为Na/Ca,红线为Mg/Ca);蓝直线表示线性拟合趋势,红色短竖线表示峰谷横坐标位置;灰色条带表示Sr/Ca与Mg/Ca低值与Na/Ca高值同期的层位,紫色条带表示Sr/Ca与Mg/Ca低值与Na/Ca低值同期的层位。Figure 4. Variations of Sr/Ca, Na/Ca and Mg/Ca ratios in the growth bands of shell carbonateThin colored lines: Original element ratio results; Thick colored lines: LOESS smoothing results(yellow: Sr/Ca, blue: Na/Ca, red: Mg/Ca); Thin blue straight lines: Linear fitting results; Red short vertical lines: the horizontal axis of the trough and peak; Grey bars shows lower values of Sr/Ca and Mg/Ca and peak values of Na/Ca, while purple bars shows the lower values of Sr/Ca, Mg/Ca and Na/Ca.2.2.1 Na/Ca比值变化趋势
研究样品的Na/Ca比值为10.23~16.19 mmol/mol,中位数及均值分别为12.5和13.6 mmol/mol。在整个生长周期中,Na/Ca比值为7.1~24.4 mmol/mol,呈现明显递减趋势(yNa/Ca=0.38x+7.13, R2=0.724 0);每期生长层的叶片方解石与白垩方解石中的Na/Ca比值波动变化,高值多落在灰色致密的叶片方解石层。
已发表[8, 22-23]的文献中,牡蛎等(Crassostrea gigas, Crassostrea virginica, Aricaica islandica)壳体Na/ Ca比值为5~30 mmol/mol,且叶片方解石层具有较高的Na/ Ca比值。本研究中的C. gigas壳体Na/ Ca平均为13.6±4.7 mmol/mol(n=8 460),变化范围为7.1~24.4 mmol/mol,与文献报道基本一致。
由于碳酸盐矿物晶格中Ca2+间的空隙决定了离子替代和共沉淀的元素种类,进而决定了不同元素对方解石的亲和性差异。对于无法进行正常的结晶沉淀的Na+来说,更可能在晶体缺陷中存在[41]。Toshihiro等[42]研究发现,各种生物骨骼微观结构中的大部分Na浓度均匀,表明环境温度和生物骨骼微观结构和钙化率对生物碳酸盐骨骼中的Na浓度影响很小;但Na/Ca比值会随埋藏年龄增加而降低,这表明在埋藏成岩作用期间生物成因碳酸钙中的结构取代的Na容易被浸出,并被释放回环境水体。Rucker等[22]发现褶牡蛎C. virginica壳体中的Na、Mg、Mn、Sr与水体盐度有密切联系(相关因子为±0.3)。
因此,本研究中Na/ Ca比值的变化原因与生物本身的渗透调节机制、盐度反应[7, 21]以及长牡蛎的方解石壳的结构变化(叶片方解石和白垩方解石交替)有密切联系。生物通过调控生物碳酸盐中的Na摄入,影响了壳体中Na+的潜在储存浓度;而牡蛎壳的两种叶片方解石和白垩方解石的两种结构交替变化,影响了壳体中Na+的实际储存含量,从而造成了壳体内不同方解石生长层之间Na/ Ca比值的差异。
2.2.2 Sr/Ca与Mg/Ca比值变化特征
本研究中的C. gigas壳体Sr/Ca比为0.76±0.12 mmol/mol(n=8 460),与已发表的牡蛎C. gigas(0.83 mmol/mol,n=268[43];0.82 mmol/mol,n=271[19]),C. virginica(0.86 mmol/mol,n=10[20]),Malleus albus(0.81 mmol/mol,n=1[44])和Saxostrea gigas(1.06 mmol/mol,n=1[44])有较好的一致性。在整个生长周期中,Sr/Ca比波动递减,每期生长层中的谷值多落在灰色致密的叶片方解石层中,对应Na/ Ca的相对高值;第6~14期Sr/Ca变化相对稳定,显著低于早期生长层;第14期(末次)生长层晚期的Sr/Ca比异常升高,可能受到后期成岩作用影响。
Mg/Ca比一般为10.2~16.2 mmol/mol,中位数及均值分别为8.3和7. 9 mmol/mol。在整个生长周期中,Mg/Ca比呈上升趋势,各期叶片方解石层均具有较低的Mg/Ca比值。已发表[13-14, 19-20, 43-45]的成体和幼体牡蛎(C. gigas, C. virginica, Malleus albus, Saxostrea gigas)壳体Mg/Ca比为5.5~27 mmol/mol,同样叶片方解石层具有较低的Mg/Ca比值。
研究认为,Sr/Ca与Mg/Ca比值在壳体中差异的决定性因素为矿物晶体结构。C. gigas壳体为低镁方解石,Ca含量为37%~40%[46],低镁方解石属三方晶系结构,Ca是六次配位,六次配位的Sr2+、Ca2+和Mg2+的半径分别为1.13×10−10、1.00×10−10和0.72×10−10 m[47],因此,Sr2+和Mg2+可以类质同象的形式替换Ca2+。一般较小的阳离子有利于六次配位[47],Ca2+会优先和Mg2+替换,因此,Sr在壳体中含量显著小于Mg。
Sr/Ca与Mg/Ca比值在各生长层中差异主要反映周围水体环境特征,即牡蛎通过方解石壳的结构变化(叶片方解石和白垩方解石交替)响应环境变化,反映在壳体不同生长层的微量元素组成上。光合生物Tridacna gigas壳中,高度分辨的Sr/Ca时间序列追踪了日照周期[16-17]。C. gigas壳中的 Mg、Sr和Na浓度之间的正相关性会随着环境温度和盐度的变化而变化[19-22]。Higuera-Ruiz等[23]观察到C. gigas壳体生长层中存在每两周一次的Mg/Ca、Sr/Ca和Na/Ca周期变化,说明壳体生长受到潮汐控制,在水温较高且盐度较高的月份具有广泛的Mg/Ca、Sr/Ca高值与Na/Ca低值,并提出壳体中Mg/Ca的最小值与大潮期间冷水流入具有直接性关联。Ullmann等[19]研究表明C. gigas壳体白垩层中的Mg/Ca比值显著高于叶片层,且与氧同位素负相关,说明Mg/Ca比值与温度的高相关性;而Sr/Ca的季节变化由代谢过程控制。
本研究区水域属中潮海岸[35],波浪小,盐度为27‰~30‰,受潮汐影响显著。因此,本研究牡蛎壳体中Sr/Ca、Mg/Ca比值的变化主要受控于壳体矿物成分和潮汐控制的环境水体温盐状况。牡蛎壳的叶片方解石和白垩方解石两种结构交替变化,壳体内不同方解石生长层Sr/Ca、Mg/Ca比值差异可以指示牡蛎生长水体环境变化。
2.3 壳体Mg/Ca比值对环境海水温度的指示
方解石在固、液相中的Mg/Ca比值与温度密切相关[48]。在10‰~35‰的盐度范围内,淡水与咸水的混合对方解石Mg/Ca比值的影响微乎其微[14],因此,通常假设Mg和Ca浓度在盐度>10‰的环境中是保守的[9, 18],故双壳类钙质壳体的Mg/Ca比主要反映了环境水体温度,如贻贝[9]、扇贝[10]、牡蛎[11-14]等。Tynan等[15]总结前人研究的9种不同水体环境下双壳类壳体Mg/Ca比值与温度的关系,发现尽管不同属种对环境水体的响应方式、敏感程度不同,相同物种在不同环境下的响应敏感程度也有差异,但Mg/Ca比值与海水温度均为正相关变化,且Mg/Ca的高值对应于较温暖的时期。研究区海门近岸水文数据[33, 35]表明:研究区冬季的平均水温为4.9 ℃,夏季为27.4 ℃,季节间温差ΔT≈22.5 ℃(图5黄色阴影区所示)。归纳总结4种基于牡蛎壳体Mg/Ca比值对生长环境温度的模型并通过理论计算,探讨不同温度校准公式在本研究的适用性情况(图5)。
图 5 依据牡蛎Mg/Ca比值的经验校正公式获得的海水表面温度比较红线: Mouchi公式[13]计算结果,橙线: Surge公式[14],蓝线: Tynan开放海洋公式[15],绿色: Tynan河口公式[15];实线代表各校正公式建立时适用的生长层范围,虚线代表适用区外延的结果;黄色阴影区表示本研究区平均海水温度器测值的变化范围[33, 35]。Figure 5. Comparison of reconstructed sea surface temperatures based on the Mg/Ca ratios in oyster shell carbonatesMouchi Equation[13] is in red, Surge Equation[14] in orange, Tynan-ocean Equation[15] in blue, Tynan-estuary Equation[15] in green; each solid line represents the applicable zone corresponding to different growing layers, respectively, while the dash line denotes the layers lack of constraints; yellow shaded area denotes mean sea surface temperature range gauged by instruments[33, 35].Mouchi等[13]研究了海洋环境下Crassostrea gigas长牡蛎壳体Mg/Ca比值与温度的相关性,并给出关系式:T=3.77 Mg/Ca+1.88(其中,温度单位为摄氏度℃,Mg/Ca比值单位为mmol/mol),温度校正公式适用于幼年生长层(≤2岁),温度为4.8~23.0 ℃。采用Mouchi(2013)海洋环境下的Mg/Ca-T校正公式(图5红线)时,若将所有生长层的Mg/Ca比值数据(n=8 270)代入,获得平均温度为31.6 ℃,温度变化为10.7~69.7 ℃,远超实际水温,说明该温度校正公式不适用于本研究中牡蛎壳体第1—14期所有生长层的数据(图5红虚线)。根据Mouchi所推荐的幼年(≤2岁)生长层范围,选择第1~2期生长层的Mg/Ca数据(n=717,图5红实线),获得季节性振幅之间的范围为17.6±1.7 ℃(夏季)和11.7±1.2 ℃(冬季),第1—2期生长层平均温度14.4±3.3 ℃,温度为10.7~22.7 ℃,季节间温差12.0 ℃。尽管第1—2期生长层的TMg/Ca符合Mouchi公式的温度范围,但由此获得的夏季平均温度较真实器测值(27.4 ℃)低9.8 ℃,冬季平均温度较真实器测值(4.9 ℃)高6.8 ℃,季节间温差较实际低10.5 ℃。
Surge等[14]研究了河口环境下Crassostrea virginica褶牡蛎壳体Mg/Ca比值与温度的相关性(T=1.39×Mg/Ca-0.35),该温度校正公式适用于成年生长层(中晚期),温度变化为19.1~31.4 ℃。采用Surge(2008)河口环境的Mg/Ca-T校正公式(图5橙线)时,选择第6~14期成年生长层的Mg/Ca数据(n=4 450),获得季节性温度范围为14.9±2.2 ℃(夏季)和9.6±1.6 ℃(冬季),第6—14期生长层(图5橙实线)平均温度11.9±3.2 ℃,温度变化为5.8~24.7 ℃,季节间温差18.9 ℃,公式不适用的第1—5期生长层Mg/Ca-SST在图5中用橙色虚线表示。由Surge公式获得的TMg/Ca温度变化范围及季节性温差均与真实器测值(4.9~27.4 ℃,ΔT≈ 22.5 ℃)非常接近。
Tynan等[15]提出海洋环境下Saccostrea glomerata悉尼岩牡蛎壳体Mg/Ca比值与温度的关系式为:T=1.41×Mg/Ca-3.25,温度为11~30 ℃;河口环境的T=1.23×Mg/Ca +2.90,温度为10~26 ℃,适用于幼年生长层(≤2年)。采用Tynan(2016)海洋环境的Mg/Ca-T校正公式(图5蓝线)时,所有生长层获得的季节性温度为11.1±2.5 ℃(夏季)和4.9±2.2 ℃(冬季);总平均温度7.9±3.8 ℃,温度变化范围为0.05~22.2 ℃,季节间温差ΔT≈21.1 ℃。由该公式获得的夏季平均温度较实际(27.4 ℃)低16.3 ℃,冬季平均温度与实际水温(4.9 ℃)吻合,季节间温差仅比实际温差略低1.4 ℃。采用Tynan(2016)河口环境的Mg/Ca-T校正公式(图5绿线)时,所有生长层获得的季节性温度为15.4±2.2 ℃(夏季)和10.0±1.9 ℃(冬季);总平均温度2.6±3.3 ℃,变化范围为5.8~25.1 ℃,季节间温差ΔT≈19.3 ℃。由该公式获得的夏季平均温度较器测值(27.4 ℃)低12.0 ℃,冬季平均温度较真实器测值(4.9 ℃)稍高0.9 ℃,季节间温差较实际低3.2 ℃。
综合分析以上4种模型,本文认为:
Mouchi公式中,仅第1—2期生长层的Mg/Ca-SST(指Mg/Ca比值计算获得的海水温度)处于器测温度范围内(图5红实线),但壳体早期生长层厚度变化大、Mg/Ca比值变化较为不稳定;而第3—14期生长层的Mg/Ca-SST(图5红虚线)绝大部分均超过30 ℃,显著偏高,与实际情况不符。研究区水体环境受到潮汐及近岸淡水输入影响,水体盐度(27‰~30‰)较开阔海洋环境(海水平均盐度约为34.7‰)低,故基于开阔海环境下Crassostrea gigas的Mouchi公式在此并不适用。
两种Tynan校正公式计算所得的Mg/Ca-SST大部分都落在研究区的海水温度变化范围内,但其温度校正公式仅适用于幼年生长层,实际有效部分仅存在于第1—2期生长层(图5绿实线、蓝实线),其对应的Mg/Ca-SST均低于10 ℃,显著偏低,与实际情况不符。基于海洋和河口环境下Saccostrea glomerata的两种Tynan公式均不适用本研究的原因,可归咎于Saccostrea属牡蛎与Crassostrea属牡蛎响应环境水体变化的敏感程度差异。
Surge公式能够覆盖Mg/Ca比值规律变化的第6—14期成年生长层(图5橙实线),计算得到的海水温度在冬夏季的平均温度、变化范围、季节性温差方面与研究区器测数据较为一致,符合研究区的水体温盐环境。我们推测,产生这一情况的原因在于受潮汐及河流淡水输入影响的研究区水体环境具有较高的水温和较大的季节性温盐波动[15],与Surge公式的盐度河口环境更为相似。因此,虽然Surge公式是基于与本文不同的牡蛎种(Crassostrea virginica)所获得的温度校正公式,但它却获得了与研究区近岸海水温度器测值更相符的结果。
3. 结论
(1)牡蛎是研究中国东部沿海地区河口浅海环境特征和海陆相互作用的良好载体,应用LA-ICP-MS分析技术可以揭示出长牡蛎壳体中高分辨率且连续的微量元素信息。
(2)环境水体物理化学性质的变化对壳体生长速率影响较大,牡蛎通过方解石壳的结构变化(叶片方解石和白垩方解石交替)响应环境变化,可反映在壳体不同生长层的微量元素组成上,本研究中发现对牡蛎季节性生长过程指示密切的元素有Na、Sr、Mg。
(3)牡蛎壳体中Mg/Ca可适用于该生物生活水体的温度重建。通过比较文献中现有的4种牡蛎壳体Mg/Ca比值古水温重建模型,发现由Surge等(2008)公式计算所得的近岸海水温度在冬夏季平均温度、变化范围以及季节性温差方面,更为符合研究区的海水温盐情况,而其他3种计算模型的结果与实际情况出入较大。产生上述结果的原因可能是研究区受到潮汐及近岸淡水输入影响,具有较高的水温和温盐波动。
(4)本文对于运用LA-ICP-MS分析技术开展河口海岸—边缘海地区高分辨率的壳体化石古环境重建具有借鉴意义。存在一些问题在未来的工作中值得进一步完善,包括牡蛎壳体年生长层的精细划分和对比、与现代观测数据的结合以及LA-ICP-MS分析技术基体匹配和激光线扫条件进一步优化等。
致谢:感谢博士后连尔刚、博士生杨承帆给予的指导和帮助。
-
图 1 孟加拉湾的表层环流及周边主要河流
橙色箭头代表西南夏季风环流,蓝色箭头代表东北冬季风环流。图形使用Ocean Data View绘制。
Figure 1. Surface circulation in the Bay of Bengal and major rivers nearby
Orange arrows represent the southwest summer monsoon currents and blue arrows represent the northeast winter monsoon currents. The figure was created by Ocean Data View.
图 2 孟加拉湾1995—2012年冬季、夏季表层海水温盐图以及硅酸盐浓度对比图
图形使用Ocean Data View绘制,数据来源于World Ocean Atlas (ODV: World Ocean Atlas 2013 (awi.de))。
Figure 2. Salinity, temperature, and comparison of silicate concentration of surface sea water in the Bay of Bengal during winter and summer of 1955—2012
Figures are created by Ocean Data View and data are from World Ocean Atlas (ODV: World Ocean Atlas 2013 (awi.de)).
图 5 季风指标与生物硅通量的比较
a. U1446站位δ18Osw值,b. U1446站位Rb/Ca比值,c. Tex86H 换算得到U1446站位表层古水温变化,d. U1446站位生物硅通量,e. U1445站位生物硅通量。a-c数据来自于Clemens等[27]以及NOAA官网(https://www.ncei.noaa.gov/access/paleo-search/study/33213),e数据来源于Lee等[11]。
Figure 5. Comparison of monsoon indicators to biogenic silica fluxes
a: δ18Osw at Site U1446, b: Rb/Ca ratio at Site U1446, c: Variations in surface paleotemperature at Site U1446 from Tex86H converted data, d: biogenic silica fluxes (Data in a–c are from Clemens et al [27] and the NOAA official website (https://www. ncei.noaa.gov/access/paleo-search/study/33213)) , e: data are from Lee et al[11].
图 8 生物硅含量及通量和地球轨道参数的比较
a. 岁差,b. 北纬15°夏季(7月)月均太阳辐射,c. 偏心率,d. 沉积速率,e. 生物硅通量。a-c数据来源于Berger[47],d-e数据均来自U1446孔。
Figure 8. Comparison of biogenic silica content and fluxes and the Earth orbital parameters
a: precession, b: average monthly solar radiation in summer (July) at 15°N, c: eccentricity, d: sedimentation rate, e: biogenic silica fluxes. Data in a–c are derived from Berger [47], d–e data are from Site U1446.
-
[1] 汪品先, 翦知湣, 赵泉鸿, 等. 南海演变与季风历史的深海证据[J]. 科学通报, 2003, 48(23):2549-2561 doi: 10.1360/03wd0156 WANG Pinxian, JIAN Zhimin, ZHAO Quanhong, et al. Evolution of the South China Sea and monsoon history revealed in deep-sea records [J]. Chinese Science Bulletin, 2003, 48(23): 2549-2561. doi: 10.1360/03wd0156
[2] Felix M. A comparison of equations commonly used to calculate organic carbon content and marine palaeoproductivity from sediment data [J]. Marine Geology, 2014, 347: 1-11. doi: 10.1016/j.margeo.2013.10.006
[3] Xu Y H, Wang L, Lai Z K, et al. The biogenic silica variation and paleoproductivity evolution in the eastern Indian Ocean during the past 20 0000 a [J]. Acta Oceanologica Sinica, 2019, 38(1): 78-84. doi: 10.1007/s13131-019-1372-z
[4] Nelson D M, Tréguer P, Brzezinski M A, et al. Production and dissolution of biogenic silica in the ocean: revised global estimates, comparison with regional data and relationship to biogenic sedimentation [J]. Global Biogeochemical Cycles, 1995, 9(3): 359-372. doi: 10.1029/95GB01070
[5] Madhupratap M, Gauns M, Ramaiah N, et al. Biogeochemistry of the Bay of Bengal: physical, chemical and primary productivity characteristics of the central and western Bay of Bengal during summer monsoon 2001 [J]. Deep Sea Research Part II:Topical Studies in Oceanography, 2003, 50(5): 881-896. doi: 10.1016/S0967-0645(02)00611-2
[6] Cortese G, Gersonde R, Hillenbrand C D, et al. Opal sedimentation shifts in the World Ocean over the last 15 Myr [J]. Earth and Planetary Science Letters, 2004, 224(3-4): 509-527. doi: 10.1016/j.jpgl.2004.05.035
[7] Volk T, Hoffert M I. Ocean carbon pumps: analysis of relative strengths and efficiencies in ocean-driven atmospheric CO2 changes[M]//Sundquist E T, Broecker W S. The Carbon Cycle and Atmospheric CO2: Natural Variations Archean to Present. Washington: American Geophysical Union, 1985: 99-110.
[8] Prasanna Kumar S, Nuncio M, Ramaiah N, et al. Eddy-mediated biological productivity in the Bay of Bengal during fall and spring intermonsoons [J]. Deep Sea Research Part I:Oceanographic Research Papers, 2007, 54(9): 1619-1640. doi: 10.1016/j.dsr.2007.06.002
[9] Phillips S C, Johnson J E, Giosan L, et al. Monsoon-influenced variation in productivity and lithogenic sediment flux since 110 ka in the offshore Mahanadi Basin, northern Bay of Bengal [J]. Marine and Petroleum Geology, 2014, 58: 502-525. doi: 10.1016/j.marpetgeo.2014.05.007
[10] Prasanna Kumar S, Muraleedharan P M, Prasad T G, et al. Why is the Bay of Bengal less productive during summer monsoon compared to the Arabian Sea? [J]. Geophysical Research Letters, 2002, 29(24): 2235.
[11] Lee J, Kim S, Khim B K. A paleoproductivity shift in the northwestern Bay of Bengal (IODP Site U1445) across the Mid-Pleistocene transition in response to weakening of the Indian summer monsoon [J]. Palaeogeography, Palaeoclimatology, Palaeoecology, 2020, 560: 110018. doi: 10.1016/j.palaeo.2020.110018
[12] Ota Y, Kuroda J, Yamaguchi A, et al. Monsoon-influenced variations in plankton community structure and upper-water column stratification in the western Bay of Bengal during the past 80 ky [J]. Palaeogeography, Palaeoclimatology, Palaeoecology, 2019, 521: 138-150. doi: 10.1016/j.palaeo.2019.02.020
[13] Clark P U, Archer D, Pollard D, et al. The middle Pleistocene transition: characteristics, mechanisms, and implications for long-term changes in atmospheric pCO2 [J]. Quaternary Science Reviews, 2006, 25(23-24): 3150-3184. doi: 10.1016/j.quascirev.2006.07.008
[14] LaFond E C. Bay of Bengal[M]//Fairbridge R W. The Encyclopedia of Oceanography. New York: Van Nostrand Reinhold Co. , 1966: 110-118.
[15] Curray J R, Emmel F J, Moore D G, et al. Structure, tectonics, and geological history of the northeastern Indian ocean[M]//Nairn A E M, Stehli F G. The Ocean Basins and Margins: The Indian Ocean. New York: Springer, 1982: 399-450.
[16] Subrahmanyam V, Subrahmanyam A S, Murty G P S, et al. Morphology and tectonics of Mahanadi Basin, northeastern continental margin of India from geophysical studies [J]. Marine Geology, 2008, 253(1-2): 63-72. doi: 10.1016/j.margeo.2008.04.007
[17] Chatterjee A, Shankar D, Shenoi S S C, et al. A new atlas of temperature and salinity for the North Indian Ocean [J]. Journal of Earth System Science, 2012, 121(3): 559-593. doi: 10.1007/s12040-012-0191-9
[18] Boyer T P, Antonov J I, Baranova O K, et al. World ocean database 2013[R]. Silver Spring: National Oceanographic Data Center Ocean Climate Laboratory, 2013.
[19] Da Silva R, Mazumdar A, Mapder T, et al. Salinity stratification controlled productivity variation over 300 ky in the Bay of Bengal [J]. Scientific Reports, 2017, 7(1): 14439. doi: 10.1038/s41598-017-14781-3
[20] Varkey M J, Murty V S N, Suryanarayana A. Physical oceanography of the bay of Bengal and Andaman sea [J]. Oceanography and Marine Biology:An Annual Review, 1996, 34: 1-70.
[21] Miranda J, Baliarsingh S K, Lotliker A A, et al. Long-term trend and environmental determinants of phytoplankton biomass in coastal waters of northwestern Bay of Bengal [J]. Environmental Monitoring and Assessment, 2020, 192(1): 55. doi: 10.1007/s10661-019-8033-8
[22] Shankar D, McCreary J P, Han W, et al. Dynamics of the East India Coastal Current: 1. Analytic solutions forced by interior Ekman pumping and local alongshore winds [J]. Journal of Geophysical Research:Oceans, 1996, 101(C6): 13975-13991. doi: 10.1029/96JC00559
[23] Shetye S R, Shenoi S S C, Gouveia A D, et al. Wind-driven coastal upwelling along the western boundary of the Bay of Bengal during the southwest monsoon [J]. Continental Shelf Research, 1991, 11(11): 1397-1408. doi: 10.1016/0278-4343(91)90042-5
[24] Clemens S C, Kuhnt W, LeVay L J, et al. Indian monsoon rainfall[C]//Proceedings of the International Ocean Discovery Program. College Station, TX: International Ocean Discovery Program, 2016: 1-15.
[25] Mortlock R A, Froelich P N. A simple method for the rapid determination of biogenic opal in pelagic marine sediments [J]. Deep Sea Research Part A. Oceanographic Research Papers, 1989, 36(9): 1415-1426. doi: 10.1016/0198-0149(89)90092-7
[26] Wang D, Ding X, Bassinot F. Observations of contrasted glacial-interglacial dissolution of foraminifera above the lysocline in the Bay of Bengal, northeastern Indian Ocean [J]. Acta Oceanologica Sinica, 2021, 40(1): 155-161. doi: 10.1007/s13131-021-1821-3
[27] Clemens S C, Yamamoto M, Thirumalai K, et al. Remote and local drivers of Pleistocene South Asian summer monsoon precipitation: a test for future predictions [J]. Science Advances, 2021, 7(23): eabg3848. doi: 10.1126/sciadv.abg3848
[28] Lisiecki L E, Raymo M E. A pliocene-pleistocene stack of 57 globally distributed benthic δ18O records [J]. Paleoceanography, 2005, 20(1): PA1003.
[29] 王汝建, 肖文申, 向霏, 等. 北冰洋西部表层沉积物中生源组分及其古海洋学意义[J]. 海洋地质与第四纪地质, 2007, 27(6):61-69 WANG Rujian, XIAO Wenshen, XIANG Fei, et al. Distribution pattern of biogenic components in surface sediments of the western arctic ocean and their paleoceanographic implications [J]. Marine Geology & Quaternary Geology, 2007, 27(6): 61-69.
[30] Berger W H. Biogenous deep sea sediments: production, preservation and interpretation[M]//Riley J P, Chester R. Treatise on Chemical Oceanography 5. San Diego, California: Academic, 1976: 265-388.
[31] Ragueneau O, Tréguer P, Leynaert A, et al. A review of the Si cycle in the modern ocean: recent progress and missing gaps in the application of biogenic opal as a paleoproductivity proxy [J]. Global and Planetary Change, 2000, 26(4): 317-365. doi: 10.1016/S0921-8181(00)00052-7
[32] 同济大学海洋地质系. 古海洋学概论[M]. 上海: 同济大学出版社, 1989: 1-316 Department of Marine Geology, Tongji University. Introduction to Paleoceanography[M]. Shanghai: Tongji University Press, 1989: 1-316.
[33] Liu S F, Zhang H, Cao P, et al. Paleoproductivity evolution in the northeastern Indian Ocean since the last glacial maximum: evidence from biogenic silica variations [J]. Deep Sea Research Part I:Oceanographic Research Papers, 2021, 175: 103591. doi: 10.1016/j.dsr.2021.103591
[34] Achyuthan H, Deshpande R D, Rao M S, et al. Stable isotopes and salinity in the surface waters of the Bay of Bengal: implications for water dynamics and palaeoclimate [J]. Marine Chemistry, 2013, 149: 51-62. doi: 10.1016/j.marchem.2012.12.006
[35] Paul J T, Ramaiah N, Sardessai S. Nutrient regimes and their effect on distribution of phytoplankton in the Bay of Bengal [J]. Marine Environmental Research, 2008, 66(3): 337-344. doi: 10.1016/j.marenvres.2008.05.007
[36] Chakrapani G J, Subramanian V. Factors controlling sediment discharge in the Mahanadi River basin, India [J]. Journal of Hydrology, 1990, 117(1-4): 169-185. doi: 10.1016/0022-1694(90)90091-B
[37] Georg R B, West A J, Basu A R, et al. Silicon fluxes and isotope composition of direct groundwater discharge into the Bay of Bengal and the effect on the global ocean silicon isotope budget [J]. Earth and Planetary Science Letters, 2009, 283(1-4): 67-74. doi: 10.1016/j.jpgl.2009.03.041
[38] Laruelle G G, Roubeix V, Sferratore A, et al. Anthropogenic perturbations of the silicon cycle at the global scale: key role of the land–ocean transition [J]. Global Biogeochemical Cycles, 2009, 23(4): GB4031.
[39] Datta D K, Gupta L P, Subramanian V. Distribution of C, N and P in the sediments of the Ganges-Brahmaputra-Meghna river system in the Bengal basin [J]. Organic Geochemistry, 1999, 30(1): 75-82. doi: 10.1016/S0146-6380(98)00203-4
[40] Sarma V V S S, Rao G D, Viswanadham R, et al. Effects of freshwater stratification on nutrients, dissolved oxygen, and phytoplankton in the Bay of Bengal [J]. Oceanography, 2016, 29(2): 222-231. doi: 10.5670/oceanog.2016.54
[41] Short D A, Mengel J G, Crowley T J, et al. Filtering of milankovitch cycles by earth's geography [J]. Quaternary Research, 1991, 35(2): 157-173. doi: 10.1016/0033-5894(91)90064-C
[42] Wang P X, Tian J, Lourens L J. Obscuring of long eccentricity cyclicity in Pleistocene oceanic carbon isotope records [J]. Earth and Planetary Science Letters, 2010, 290(3-4): 319-330. doi: 10.1016/j.jpgl.2009.12.028
[43] Araya-Melo P A, Crucifix M, Bounceur N. Global sensitivity analysis of the Indian monsoon during the Pleistocene [J]. Climate of the Past, 2015, 11(1): 45-61. doi: 10.5194/cp-11-45-2015
[44] Erb M P, Jackson C S, Broccoli A J. Using single-forcing GCM simulations to reconstruct and interpret Quaternary climate change [J]. Journal of Climate, 2015, 28(24): 9746-9767. doi: 10.1175/JCLI-D-15-0329.1
[45] Kim S, Takahashi K, Khim B K, et al. Biogenic opal production changes during the Mid-Pleistocene Transition in the Bering Sea (IODP Expedition 323 Site U1343) [J]. Quaternary Research, 2014, 81(1): 151-157. doi: 10.1016/j.yqres.2013.10.001
[46] Iwasaki S, Takahashi K, Kanematsu Y, et al. Paleoproductivity and paleoceanography of the last 4.3 Myrs at IODP Expedition 323 Site U1341 in the Bering Sea based on biogenic opal content [J]. Deep Sea Research Part II:Topical Studies in Oceanography, 2016, 125-126: 145-154. doi: 10.1016/j.dsr2.2015.04.005
[47] Berger A. Orbital variations and insolation database[R]. Boulder CO, USA: National Oceanic and Atmospheric Administration, 1992.
-
期刊类型引用(2)
1. 常健玮,刘国群,颜雯婷,刘洋,徐铖,陈敏,秦华. 牡蛎壳粉替代石灰对橘园土壤酸化改良及果实品质的影响. 中国土壤与肥料. 2024(10): 136-141 . 百度学术
2. 林梵宇,尹希杰,黄威,黄杰超,梁毓娜. 利用微区XRF技术的大洋固体矿产成分快速无损检测. 海洋地质与第四纪地质. 2021(01): 223-232 . 本站查看
其他类型引用(1)