南黄海牡蛎壳元素组成的原位微区分析及环境指示

宋佳泽，黄湘通，杨守业，亓利剑

同济大学海洋地质国家重点实验室，上海 200092

近年来，随着微区分析技术的快速发展，应用激光剥蚀、离子探针等对生物壳体的原位与高分辨分析受到广泛重视。原位微区分析是利用电子束轰击或激光束剥蚀进行采样，实现原位元素—同位素组成分析。自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 2008 年11 月江苏省海门市近岸牡蛎滩采样位置(A)、研究区2008 年(B)与2003 年(C)的卫星影像Fig.1 Sampling location （A） and satellite image （B/C） of oyster bank in the offshore of Haimen City, Jiangsu Province.Google Earth satellite images show the environment change from 2003 to 2008

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 和结束后的载气平均响应值作为空白背景。

图2 长牡蛎样品左壳(A)、壳体韧带槽切片(B)与韧带槽切片剖面上的生长层与线扫描轨迹(C)Fig.2 Photos showing the Pacific oyster left shell （A） used in the analysis, cross section of the shell ligamental area （B） and shell growth layers and trajectory of line-scanning （C）

1.3 数据校正

由ICP-MS 直接获取的数据为各质量数的测试响应值（单位为counts per second，即每秒的计数值，记作cps），在扣除背景空白、获取有效的计数率信号cps 后，研究采用Longerich 的定量方法[38]对cps进行校正，得到元素比值，具体公式如下：

公式（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σ，可以定性分析变化趋势。

图3 LA-ICP-MS 线剥蚀测试长牡蛎壳体样品的3 种变化模式图Fig.3 Three typical variation models of elemental abundances in the Pacific oyster shell under LA-ICP-MS line-scan mode

（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）。

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-10m[47]，因此，Sr2+和Mg2+可以类质同象的形式替换Ca2+。一般较小的阳离子有利于六次配位[47]，Ca2+会优先和Mg2+替换，因此，Sr 在壳体中含量显著小于Mg。

图4 牡蛎样品不同生长层中Na/Ca、Mg/Ca、Sr/Ca 比值变化Fig.4 Variations of Sr/Ca, Na/Ca and Mg/Ca ratios in the growth bands of shell carbonate

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）。

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 ℃。

图5 依据牡蛎Mg/Ca 比值的经验校正公式获得的海水表面温度比较Fig.5 Comparison of reconstructed sea surface temperatures based on the Mg/Ca ratios in oyster shell carbonates

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-ICPMS 分析技术可以揭示出长牡蛎壳体中高分辨率且连续的微量元素信息。

（2）环境水体物理化学性质的变化对壳体生长速率影响较大，牡蛎通过方解石壳的结构变化（叶片方解石和白垩方解石交替）响应环境变化，可反映在壳体不同生长层的微量元素组成上，本研究中发现对牡蛎季节性生长过程指示密切的元素有Na、Sr、Mg。

（3）牡蛎壳体中Mg/Ca 可适用于该生物生活水体的温度重建。通过比较文献中现有的4 种牡蛎壳体Mg/Ca 比值古水温重建模型，发现由Surge 等（2008）公式计算所得的近岸海水温度在冬夏季平均温度、变化范围以及季节性温差方面，更为符合研究区的海水温盐情况，而其他3 种计算模型的结果与实际情况出入较大。产生上述结果的原因可能是研究区受到潮汐及近岸淡水输入影响，具有较高的水温和温盐波动。

（4）本文对于运用LA-ICP-MS 分析技术开展河口海岸—边缘海地区高分辨率的壳体化石古环境重建具有借鉴意义。存在一些问题在未来的工作中值得进一步完善，包括牡蛎壳体年生长层的精细划分和对比、与现代观测数据的结合以及LA-ICPMS 分析技术基体匹配和激光线扫条件进一步优化等。

致谢：感谢博士后连尔刚、博士生杨承帆给予的指导和帮助。