沉积物粒度是地球科学研究中应用最为广泛的指标之一,端元组分的应用则有助于理解记录在粒度信息中的不同沉积物源、沉积动力过程的匹配关系。在沉积学方面,试图进行粒度数据分解已有很长的历史,端元值法尽管存在分解潜能,或出于有效分析软件的限制,这一极具潜力的粒度分析方法并未得到广泛应用。
针对这一现状,地质地球所地球深部结构与过程研究室Greig A. Paterson副研究员与其在澳大利亚国立大学的合作者David Heslop博士,设计开发了一种基于总体特征的粒度端元分析的非参数估计方法。对比单个样品的数学分解过程,此前国际上的算法表现为求解域的非唯一性。为克服这一缺陷,新型算法首先将整个数据集作为整体进行计算,通过获得数个整体的端元函数,提高对单一样本中的具体端元组分的分析辨识度。相较于目前国际上的一些已有沉积数据分离算法,该方法在数据拟合方面具有更高的准确率。
此外,为帮助研究人员应用该算法,他们基于MATLAB的GUI研发了一个新软件(AnalySize),同时集成了国际上已有的一些粒度分析技术,包括三端元相图,CM图和一些描述性统计量。该软件具有交互式操作界面,并支持多种格式的数据输入和图件输出。
图1 200组样本的对数正态分布的分析结果及其拟合优度特征。a-c,d-f和g-i分别为2、3、4端元组分的拟合结果,HALS-NMF是新算法
图2 基于不同非参数算法的4组分拟合结果对比。a-d,拟合与真实组分的对比,e-h,各组分相对含量的对比
以上研究工作已于近期发表在国际著名地学杂志Geochem, Geophys, Geosyst(Greig A et al. New methods for unmixing sediment grain size data, Geochem, Geophys, Geosyst, 2015, 16: 4494–4506)。
Copyright © 古地磁与年代学实验室 京ICP备05029136号 地址: 北京市北土城西路19号古地磁实验室 邮编:100029 电话:010-82998366 |