基于GEOROC大数据分析地壳厚度地球化学指标

摘要

前人研究认为,火山岩中部分地球化学指标与岩浆弧地壳厚度之间存在一定的相关性,并通过统计主量元素K2O、CaO和Na2O指标及微量元素Ce/Y、Sm/Yb、Dy/Yb、Sr/Y、La/Yb指标与地壳厚度之间关系,约束地质史上某些区域的地壳厚度发展和变化.本文基于GEOROC数据库,以SiO2含量57％和火山岩年龄23Ma为界,将全球火山岩数据分成年轻-壳源(＞57％,＜ 23Ma)、年轻-幔源(＜57％,＜23Ma)、古老-壳源(＞57％,＞23Ma)和古老-幔源(＜57％,＞23Ma)四个数据集,并通过核函数估计方法获得了各个地球化学指标与地壳厚度的归一化联合概率密度分布图.本文统计结果表明,年轻-幔源火山岩中的K2O含量分布与壳源火山岩呈现指数正相关关系、CaO含量分布于地壳厚度呈现线性负相关关系,年轻-壳源火山岩中Ce/Y、La/Yb和Sm/Yb与现今地壳厚度有指数正相关关系.由以上5种地化指标建立的回归方程确定系数R2均大于0.7,可以认为相关关系显著.本文认为幔源岩浆在穿透地壳到达地表过程中,地壳厚度控制了富K壳源物质进入地幔熔体和富Ca矿物结晶分异过程,导致了火山岩中K2O和CaO含量的相关变化;而下地壳部分熔融形成的壳源岩浆,不同深度压力控制了残留相矿物比例,导致Ce/Y、La/Yb和Sm/Yb体现出与地壳厚度的相关性.本文建立的回归函数是基于大量数据概率密度分布的统计分析得出的,由于离群数据普遍存在,回溯历史地壳厚度变化需要大量数据统计支撑,否则难以获得可靠的结果.