碳酸盐 δ¹⁸O 古温度计:从氧同位素反演古温度怎么做
碳酸盐(有孔虫、贝壳、石笋、碳酸盐岩等)的 δ¹⁸O 随结晶温度变化,是重建古温度与古水体最经典的地球化学工具之一。 本文梳理其原理、常用换算方程、计算步骤与必须交代的前提假设与局限,供古气候/古海洋研究设计与数据解读参考。
一、原理:温度控制的碳酸盐–水氧同位素分馏
碳酸盐结晶时,其 δ¹⁸O 相对成矿水的分馏随温度升高而减小(温度越高、碳酸盐 δ¹⁸O 越偏负)。 这一"碳酸盐–水氧同位素温度计"由 Urey(1947)提出。关键前提:碳酸盐 δ¹⁸O 常以 VPDB 报告、水 δ¹⁸O 以 VSMOW 报告, 套用方程前必须正确换算两套标尺(见 δ¹⁸O 标度换算)。
二、常用方程(以 Kim & O'Neil 1997 为例)
文献有多套实验/理论/经验标定方程。低温无机方解石最常被引用的是 Kim & O'Neil(1997):
1000 · ln α(方解石–水) = 18.03 · (10³/T) − 32.42 (T 为开尔文温度 K)
- α 为方解石与水之间的 ¹⁸O/¹⁶O 分馏因子;该式基于 10、25、40°C 合成实验标定,适用约 10–40°C。
- 文石需校正:25 ℃ 下文石相对方解石富集 ¹⁸O 约 0.6–0.8‰——Tarutani 等(1969)实测 0.6‰,Kim 等(2007)以新的文石–水标定配 Kim 与 O'Neil(1997)方解石–水标定得 0.78‰。用方解石方程处理文石数据须相应扣除;更稳妥的做法是直接用文石–水方程 1000·lnα = 17.88·(10³/T) − 31.14(Kim 等 2007,0–40 ℃)。
- 其他常见标定:O'Neil 等 1969、Coplen 2007(Devils Hole 天然标定)、Watkins 等 2014(pH/生长速率动力学)、Daëron 等 2019 等(Friedman 与 O'Neil 1977 为分馏因子汇编,非原始标定)。不同方程给出的温度差异并不小:Coplen(2007)指出沿用此前发表的分馏因子会使 Devils Hole 方解石的反演温度比实测地下水温低 8 ℃;Daëron 等(2019)测得的平衡 ¹⁸α(方解石/水) 比传统标定高约 1.5‰(25 ℃ 附近约合 7 ℃)。因此报告必须写明所用方程及其 AFF 前提,跨研究比较前先统一方程。
三、计算步骤
- 测定碳酸盐 δ¹⁸O(VPDB),并换算到与水一致的 VSMOW 标尺以计算 α。
- 确定或假设成矿水 δ¹⁸O(VSMOW)——这是最大的不确定来源(见下)。
- 由碳酸盐与水的 δ¹⁸O 求 1000·ln α,代入方程反解温度 T。
- 据矿物类型(文石/方解石)与生物类群做必要的分馏/"生命效应"校正。
四、必须交代的前提与局限
| 问题 | 影响 | 应对 |
|---|---|---|
| 成矿水 δ¹⁸O 未知 | 温度反演强依赖水 δ¹⁸O 假设,海水 δ¹⁸O 随地质时代与冰量变化 | 独立约束水 δ¹⁸O;或用团簇同位素 Δ₄₇ 无需假设水值 |
| 盐度/局地水 δ¹⁸O | δ¹⁸O_sw 与盐度的关系是局地性的、随水团与淡水收支而异,用全球平均值处理近岸、河口或表层样品会产生系统偏差 | 优先用就地实测或区域 δ¹⁸O_sw–盐度关系(如 LeGrande 与 Schmidt 2006 的全球网格数据集),并注明所用关系式 |
| 冰量效应 | 大陆冰盖优先封存 ¹⁶O、抬高残余海水 δ¹⁸O(无冰海洋常取 δ¹⁸O_sw ≈ −1.2‰ SMOW(Zachos 等 2001 的温标即按此计算),亦有文献取 −1.0‰,取值不同会直接平移反演温度) | 冰期区间需做冰量校正,无法约束时结果不可靠 |
| 成岩蚀变 | 海底重结晶会重置原生信号,且可隐性发生——Sexton 等(2006)指出反射光下呈霜状(frosty)的浮游有孔虫壳体虽形态完整,δ¹⁸O 已系统偏正、反演温度偏低,这正是白垩纪—古近纪"凉热带悖论"的成因 | 优选反射光下呈玻璃状(glassy)的壳体,配合 Δ₄₇ 与微区筛查 |
| 动力学效应/生命效应 | 不止生物碳酸盐:Watkins 等(2014)在无机方解石生长实验中即观察到 >2‰ 的非平衡氧同位素效应(25 ℃ 附近约合 10 ℃),且随 pH 与晶体生长速率系统变化;生物类群还叠加共生藻光合、个体发育等效应 | 选合适类群与对应校正、注明所用校正与出处;对高生长速率或高 pH 体系应谨慎解读,或改用 Δ₄₇ |
五、方法学提示
- 酸解分馏因子(AFF)一致性:碳酸盐经磷酸消解后实测的是 CO₂,换算回碳酸盐 δ¹⁸O 需用酸解分馏因子 α(CO₂(酸)–方解石)。IUPAC/IAEA 现行共识取 25 ℃ 方解石 AFF = 1.01025("公认值、假定精确");而 Kim 与 O'Neil(1997)当年用的是 1.01050,两者算出的碳酸盐 δ¹⁸O 相差约 0.24‰(25 ℃ 附近约合 1 ℃)。因此把按现行标准归一化的数据直接代入 1997 年原式会引入系统偏差,跨方程比较前须按同一 AFF 重算。另须注意 1.01025 与 1.010254 来源不同、不可混称——后者是 Kim、Coplen 与 Horita(2015)温度式在 25 ℃ 的计算值。
- 团簇同位素 Δ₄₇:可独立给出结晶温度、无需假设水 δ¹⁸O,并能反算水值、辅助成岩筛查,是近年常用的互补/进阶手段(属专门方法)。
- 交叉验证:磷酸盐 δ¹⁸O 等独立代用指标可与碳酸盐 δ¹⁸O 相互印证。
六、检测平台
碳酸盐 δ¹³C·δ¹⁸O 采用磷酸法在线气体制备—同位素比值质谱(GasBench II + MAT 253 Plus)测定,以 VPDB 报告。 核素科技提供碳酸盐(有孔虫、贝壳、石笋、碳酸盐岩等)δ¹³C·δ¹⁸O 检测,可用于古温度/古环境重建的数据获取; 古温度反演所需的水 δ¹⁸O 假设、方程选择与校正由研究者按体系确定,具体送样量与前处理请与技术顾问确认。
参考文献
- Kim S.-T., O'Neil J.R. (1997) Equilibrium and nonequilibrium oxygen isotope effects in synthetic carbonates. Geochimica et Cosmochimica Acta 61(16):3461-3475.(本页方解石–水方程 1000·lnα=18.03·10³/T−32.42 的原始出处;10/25/40 ℃ 合成实验标定,实测值;该文明载所用酸解分馏因子为 1.01050) https://doi.org/10.1016/S0016-7037(97)00169-5
- Kim S.-T., O'Neil J.R., Hillaire-Marcel C., Mucci A. (2007) Oxygen isotope fractionation between synthetic aragonite and water: Influence of temperature and Mg2+ concentration. Geochimica et Cosmochimica Acta 71(19):4704-4715.(文石–水方程 1000·lnα=17.88(10³/T)−31.14,0–40 ℃,实测值;并给出 25 ℃ 文石–方解石分馏 0.78‰) https://doi.org/10.1016/j.gca.2007.04.019
- Tarutani T., Clayton R.N., Mayeda T.K. (1969) The effect of polymorphism and magnesium substitution on oxygen isotope fractionation between calcium carbonate and water. Geochimica et Cosmochimica Acta 33(8):987-996.(25 ℃ 文石较方解石富集 ¹⁸O 约 0.6‰ 的实测原始出处;本站未取得原文全文,数值经 Kim 与 O'Neil 1997、Kim 等 2007 原文逐字转引核实) https://doi.org/10.1016/0016-7037(69)90108-2
- Coplen T.B. (2007) Calibration of the calcite-water oxygen-isotope geothermometer at Devils Hole, Nevada, a natural laboratory. Geochimica et Cosmochimica Acta 71(16):3948-3957.(天然标定实测值;指出沿用旧分馏因子会使反演温度比实测地下水温低 8 ℃——本页"不同方程差异并不小"的依据。本站打开的是 USGS 官方著录与摘要,未取得期刊原文全文) https://pubs.usgs.gov/publication/70030756
- Daëron M., Drysdale R.N., Peral M., Huyghe D., Blamart D., Coplen T.B., Lartaud F., Zanchetta G. (2019) Most Earth-surface calcites precipitate out of isotopic equilibrium. Nature Communications 10:429.(平衡 ¹⁸α(方解石/水) 较传统标定高约 1.5‰;并证实生物成因方解石 Δ₄₇ 亦存在细微非平衡偏差。实测值,开放获取) https://pmc.ncbi.nlm.nih.gov/articles/PMC6347637/
- Watkins J.M., Hunt J.D., Ryerson F.J., DePaolo D.J. (2014) The influence of temperature, pH, and growth rate on the δ18O composition of inorganically precipitated calcite. Earth and Planetary Science Letters 404:332-343.(无机方解石生长实验实测 >2‰ 的非平衡氧同位素效应,随 pH 与生长速率系统变化。开放获取版为 OSTI 托管的作者稿) https://www.osti.gov/servlets/purl/1466154
- Camin F. 等 (2025) Stable Isotope Reference Materials and Scale Definitions — Outcomes of the 2024 IAEA Experts Meeting. Rapid Communications in Mass Spectrometry 39(14):e10018.(逐字重申 25 ℃ 方解石酸解分馏因子 1.01025 为"公认值、假定精确"并为 IUPAC 推荐。约定共识值,非实测值。开放获取) https://eprintspublications.npl.co.uk/10390/1/eid10390.pdf
- LeGrande A.N., Schmidt G.A. (2006) Global gridded data set of the oxygen isotopic composition in seawater. Geophysical Research Letters 33(12):L12604.(1°×1° 全球海水 δ¹⁸O 网格数据集,由观测结合区域 δ¹⁸O–盐度关系构建,文献汇编值;据以说明 δ¹⁸O_sw–盐度关系的局地性。本站打开的是著录与摘要,未取得原文全文) https://doi.org/10.1029/2006GL026011
- Zachos J., Pagani M., Sloan L., Thomas E., Billups K. (2001) Trends, Rhythms, and Aberrations in Global Climate 65 Ma to Present. Science 292(5517):686-693.(其 Fig.2 温标按无冰海洋 δ¹⁸O_sw ≈ −1.2‰ SMOW 计算——本页冰量效应取值的依据。假设约定值,非实测值) https://doi.org/10.1126/science.1059412
- Sexton P.F., Wilson P.A., Pearson P.N. (2006) Microstructural and geochemical perspectives on planktic foraminiferal preservation: "Glassy" versus "Frosty". Geochemistry, Geophysics, Geosystems 7(12):Q12P19.(海底成岩改造使热带浮游有孔虫 δ¹⁸O 系统偏正、古温度偏低。实测值。本站打开的是著录与摘要,未取得原文全文) https://doi.org/10.1029/2006GC001291
- Urey H.C. (1947) The thermodynamic properties of isotopic substances. Journal of the Chemical Society 562-581.(碳酸盐氧同位素受结晶温度控制这一预言的原始出处,理论推导值;本站未取得原文全文,卷页与归属经 Daëron 等 2019 开放全文的参考文献表逐字核实) https://doi.org/10.1039/jr9470000562
常见问题
碳酸盐 δ¹⁸O 为什么能当古温度计?
碳酸盐结晶时其 δ¹⁸O 相对成矿水的氧同位素分馏随温度变化——温度越高、碳酸盐相对水越贫 ¹⁸O。测得碳酸盐 δ¹⁸O 并已知或假设成矿水 δ¹⁸O,代入标定方程(如 Kim & O'Neil 1997:1000·lnα=18.03·10³/T−32.42,T 为开尔文)即可反解温度。注意碳酸盐 δ¹⁸O 常以 VPDB 报告、水以 VSMOW 报告,须先正确换算标尺。
用碳酸盐 δ¹⁸O 反演温度最大的坑是什么?
最大不确定来自成矿水 δ¹⁸O:它不被独立记录,且海水 δ¹⁸O 随地质时代与冰量变化(大陆冰盖封存 ¹⁶O 会抬高海水 δ¹⁸O)。其次是成岩重结晶重置信号、部分生物碳酸盐的"生命效应"偏离无机平衡、以及矿物学差异(文石较方解石约 +0.6~0.8‰)。因此报告须写明所用水 δ¹⁸O 假设、标定方程与校正。
团簇同位素 Δ₄₇ 和传统 δ¹⁸O 古温度计什么关系?
Δ₄₇(团簇同位素)测定碳酸盐中 ¹³C–¹⁸O 成键的富集程度,主要受结晶温度控制、不需假设成矿水 δ¹⁸O(但 Daëron 等 2019 指出生物成因方解石的 Δ₄₇ 也存在细微非平衡偏差)、不需假设成矿水 δ¹⁸O,因而能独立给出温度并反算水 δ¹⁸O,还可辅助识别成岩蚀变,是传统 δ¹⁸O 温度计的互补/进阶手段(属专门方法)。
相关内容
首次发布:2026-07-17