地球矿物表面的轨道成像光谱学：EMIT前两年的结果与稳定性

《Remote Sensing of Environment》：Orbital imaging spectroscopy of Earth’s mineral surface: Results and stability from EMIT’s first two years

　　

研究人员评估了NASA地球矿物粉尘源调查(EMIT)轨道成像光谱仪前两年的历史矿物填图数据，开发了隐马尔可夫模型(HMM)框架以区分瞬态噪声与真实矿物变化。组成分类具有80–98%的可重复性一致性。研究人员对时间序列进行了表面矿物学变化检测测试，并引入统计方法对时间序列数据去噪，生成统一的高可靠性矿物图。可见短波红外(VSWIR)成像光谱利用太阳反射光谱探测赤铁矿、针铁矿、粘土和碳酸盐等矿物的特征电磁吸收。EMIT在国际空间站上以公顷级空间采样首次实现大尺度地表矿物学填图，绘制十类关键矿物分类及其相对丰度。由于仪器噪声、大气伪影、定位误差、观测几何、土壤水分、植被覆盖与校准漂移等影响，确定性算法在重复观测中仍会出现表观变化。该工作利用统计模型消解上述效应，回答分类精度、能否探测矿物学变化、以及如何将含噪时间序列融合为统一矿物图三个问题。研究人员在已知稳定点与预期变化点评估EMIT序列，发现矿物分类一致性达80–98%，即每5–10次过空中约1次为假变化；大量重复观测可进一步降低方差，且该方法能量化表观变化与测量噪声的一致性，为EMIT产品可靠性解释提供依据。

研究背景方面，可见短波红外(VSWIR)成像光谱技术长期用于矿物填图，从早期站点演示发展到现代广域调查，利用太阳反射率光谱探测铁氧化物、粘土和碳酸盐等矿物的特征吸收。NASA地球矿物粉尘源调查(EMIT)作为国际空间站上的成像光谱仪，以60米像素、75公里幅宽在380–2500纳米范围以7.4纳米采样，首次实现公顷级大陆尺度地表矿物学填图，绘制十类关键矿物分类及相对丰度，已被用于改进矿物粉尘源理解。尽管矿物填图算法是确定性的，但仪器噪声、大气伪影、定位误差、观测几何导致表土层不同深度矿物表观差异、太阳角改变投影阴影、土壤水分修饰液态水特征、植被覆盖掩蔽、仪器校准漂移，以及季节性降水、风成粉尘沉积、滑坡、风化等真实过程，皆可使重复观测出现分类波动；此外多吸收特征重叠与混合像元丰度反演不确定也造成解释歧义。因此需统计框架区分瞬态噪声与真实变化，评估精度、探测能力及多时相融合策略。

关键技术方法方面，研究人员选取全球39个矿物学稳定站点（每边超100公里矩形，含2022–2025年2900景EMIT场景，覆盖干旱区、稳定高原与定标场），云掩膜后对齐至统一网格构建经纬度-时间三维分类数组；采用Tetracorder算法5.27c1对大气校正后方向反射率作连续统去除吸收特征识别，聚为10类矿物（1微米组赤铁矿与针铁矿，2微米组方解石、绿泥石、白云石、石膏、伊利石+白云母、高岭石、蒙脱石、蛭石）；基于稳定点假定真值不变，以众数分类代真值，用加一平滑经验计数构建全局观测（混淆）矩阵；设极低状态转移概率（0.001%）表征地表转换千年尺度罕见性，对每像素两条独立HMM（1微米组2态、2微米组8态）用Viterbi解码求最大后验历史、前向-后向算法求边缘概率，在埃塞俄比亚Erta Ale火山站做变化检测验证。

研究结果部分，3.1检测噪声中，研究人员由稳定点经验条件概率得两组混淆矩阵：1微米组赤铁矿最稳，针铁矿因与赤铁矿可见段形状相似及混合像元易振荡；2微米组中蛭石精度98%最高，高岭石92%、蒙脱石93%、白云石近93%，石膏约84%最低（常误为蒙脱石，因同含水OH吸收及蒸发盆地共生），方解石主要误为伊利石+白云母（次生特征与碳酸盐键重叠及沉积共生），绿泥石误为伊利石+白云母与蒙脱石（混合与2.3微米大气校正伪影）。3.2隐马尔可夫模型结果中，在Erta Ale典型像素，原始1微米分类在赤/针铁矿间振荡，HMM据噪声矩阵与极低转移先验判定全时相赤铁矿；2微米组某像素在伊利石+白云母、绿泥石、方解石间振荡，HMM给伊利石+白云母66%后验且全序恒定，示范去噪融合图较单期最新像更准，弱吸收与混合边界处改善明显。3.3测量变化中，新鲜火山物质首末过空无变化；孤立单像素跳变归为混合振荡假阳性；但检出十馀邻接像素由赤铁矿转针铁矿的聚类变化，连续统去除后铁氧化物吸收变宽、反射率变暗，位于冲积扇与河床反照率变动区，机制可能为降雨沉积/侵蚀，HMM后验远出噪声界；近火口无变化点光谱浅吸收（<1%）与残差水汽/散射一致，支持无变化判据。

讨论与结论部分，研究人员指出EMIT分类一致性80–98%，用户应取全任务时序而非单景；HMM提供观测矩阵与序贯推断工具，无真实变化时退化为众数，有变化时可显式测试。方差主源为仪器噪声、定位、大气残余，混合像元使边界振荡。稳定干旱点隔离了系统噪声，但土壤水分、亚像元植被等环境噪声随地点变，通用混淆矩阵仅基线。该方法捕方差非偏倚，偏倚需地面真值评辐射/反射率/矿物图系统性误差。HMM可扩至特征深度回归（如Kalman滤波）与更多矿物类或生态时序；每点独立推断高效精确，超越一步Markov或可达更长记忆；均匀低转移概率仅为统计上界，其他领域须按历时动态化。结论重申：EMIT前两年矿物填图重复一致性80–98%，HMM框架能分瞬态噪声与真变化，多时相融合生成高可靠统一矿物图，为解释L2B产品与粉尘辐射强迫研究提供方法学基础。