地磁场是地球的基本物理场之一，起源于地球液态外核，能够反映地球内部结构和运动状态，同时环绕在地球周围，保护着地球宜居环境。研究地磁场变化规律在认识地球深部动力机制、地表过程和环境、现代地磁场特征和变化趋势等方面具有重要意义。大量不同时空分布的高质量观测数据是全面认识地磁场区域和全球变化特征的必要条件。近现代的卫星观测数据满足这一要求，但覆盖时间短。考古磁学可以把地磁场观测记录回溯至近几千年甚至全新世，在建立地磁场区域和全球模型方面具有重要贡献。然而，目前全球考古磁学数据分布严重不均，大部分数据集中在欧洲和中东地区，北半球低纬度地区和南半球数据稀少，制约了对地磁场形态特征和变化规律的准确认识。 

　　东南亚地区正是上述全球数据匮乏区之一。中科院地质与地球物理研究所蔡书慧副研究员联合美国加州大学圣迭戈分校Lisa Tauxe教授、美国芝加哥伊利诺大学考古学家Mitch Hendrickson副教授、澳大利亚¹⁴C定年专家Quan Hua博士和柬埔寨考古项目组专家等，对位于东南亚柬埔寨北部的一处古代冶铁遗址（Tonle Bak）开展了考古磁学研究，获得了该地区首批地磁场全矢量数据。研究结果不仅大大改善了全球考古磁学数据结构，还记录到地磁场负异常，为认识低纬度地区地磁场特征及核幔边界结构提供了新思路，并为今后相关区域考古磁学工作提出了新方向。 

　　Tonle Bak遗址位于著名的吴哥窟古城以东～135 km（图1a），考古调查结果显示它所在的Phnom Dek地区很可能为吴哥王朝扩张时期很重要的铁金属供应地，在该地区发现了十几处冶铁遗址，本文采样区Tonle Bak遗址是其中一处（图1b）。本研究对出土的3个冶炼窑址（Furnace1，Furnace2，Furnace3）和3个发掘沟堆积剖面（Trench1，Trench3，Trench4）系统采集了定向样品（窑底或窑壁固结烧土）、非定向样品（烧土、吹风管、炉渣）和碳屑（图2）。碳屑的¹⁴C定年结果将遗址年龄范围限定为公元1034-1391年。窑址定向样品记录的地磁场方向结果结合邻近地区已发表数据显示地磁场方向在公元1200-1300年间发生突变（图3a），即使变化速率下限也高于该地区现代地磁场方向平均变化速率。同时样品强度结果显示公元1100-1300年地磁场强度存在不对称“V”形变化，最低值出现在公元1200-1300年（图3b），与方向结果对应。 

　　 

　　　　图1 （a）研究区位置图。红色五角星为采样位置，绿色圆形为吴哥窟位置。其余符号为本文用到的已发表数据采点位置，详见图注。（b）采样遗址高程图。图中标明了本文采样的3个窑址和3个发掘沟位置 

 

图2 （a）代表性窑址（F3）定向样品采样照片。图中字母代表定向样品采集位置；（b）代表性发掘沟（T1）非定向样品和碳屑采样照片。蓝色数字为非定向样品采集位置，旁边括号中为样品记录的古强度值，红色数字代表碳屑采集位置 

　 

图3 （a）地磁场方向结果。红色为本文新增结果，黑色为泰国已发表结果，蓝色虚线代表地磁场方向移动轨迹；（b）地磁场强度结果。红色五角星为本文新增结果，其余符号为已发表数据，详见图注。桔色、灰色和黄绿色实线分别为中国考古强度参考曲线和全球模型CALS10k.2和SHAWQ2k对采样位置预测结果 

　　Tarduno et al.（2015）在非洲南部的考古磁学研究发现类似的地磁异常，公元1300年左右地磁场方向突变伴随地磁场强度降低，研究区位于非洲下部低剪切速度地幔异常体（LLSVP）陡峭边界的正上方，因此作者将这一磁异常归因于核幔边界处地形异常引起的磁通量溢出（magnetic flux expulsion）产生的反向涡流（reversal flux patch）。本研究获得的柬埔寨考古磁学结果与非洲南部结果非常相似，且两个研究区都位于全球模型预测的赤道附近地磁场负异常区（EGA，图4a），指示二者可能有相似的成因。对比EGA覆盖范围内其它地区地磁场强度结果，发现各个地区地磁场强度在公元1100-1400年都有降低的趋势，只是降低的幅度和最低值出现的时间存在差异（图4b-e），指示低纬度地区核幔边界处可能普遍存在磁通量溢出现象。 

　　  

　　　图4 （a）全球模型SHAWQ2k预测的公元1250年地表地磁场强度分布；（b-e）非洲、东亚-东南亚、南美和西南太平洋地区公元1000-1500年地磁场强度数据，数据来源详见图注  

　　由于各个地区已有数据质量和数量不一致，目前对各地区磁异常形态特征还没有足够准确的认识。它们是同时出现还是相互之间有传递关系？是否反映该区域核幔边界结构比较复杂？它们在时间上是否具有周期性？与现代地磁场的显著特征——南大西洋异常（SAA）——有何关联？对认识现代地磁场变化趋势有何启示？回答上述问题，需要更多低纬度地区不同时间的高质量数据，同时也需要结合核幔边界结构和流体运动状态方面的相关模拟计算工作来解释其形成机制和变化规律。 

　　研究成果发表于国际顶级学术期刊PNAS (蔡书慧, Doctor R, Tauxe L, Hendrickson M, Hua Q, Leroy S, Phon K. Archaeomagnetic results from Cambodia in Southeast Asia: Evidence for possible low-latitude flux expulsion[J]. Proceedings of the National Academy of Sciences, 2021, 18(11), DOI: 10.1073/pnas.2022490118)。本研究受国家自然科学基金委、美国国家科学基金会、中国科学院人才项目和研究所重点部署项目等资助。