热力学计算模拟对初始月幔结构的约束

月球岩浆洋结晶形成的初始月球内部结构是其后续演化过程的开端,其结晶过程受月球岩浆洋的初始深度和物质组成这两个参数的制约。由于缺少直接来自月球深部的岩石样品,目前关于月球岩浆洋演化过程的探讨主要依赖实验和计算模拟手段。岩浆洋模型中形成的月壳厚度是否与探测结果一致是月球岩浆洋演化模型合理性的重要约束。最新的GRAIL(Gravity Recovery and Interior Laboratory)探测数据推算月壳厚度为34~43km,低于阿波罗时期认为的约70km,这对已有的月球岩浆洋演化模型提出了挑战。本文采用并修正FXMOTR程序包,针对月球岩浆洋在不同的初始深度和物质组成情况下的结晶过程,进行了一系列热力学计算模拟。通过量化月球岩浆洋的初始深度和物质组成对月壳厚度的影响,结合关于月球内部微量元素分配的研究结果,对比了月球岩浆洋结晶后期的残余熔体与原始克里普组分(urKREEP)的成分。本文的模拟结果显示,一个全月幔熔融且初始成分为月球初始上月幔组成(LPUM)的岩浆洋将在其深部结晶2.5%石榴子石,形成的月壳厚度符合GRAIL的约束,并且结晶出了合适的urKREEP成分。在此模型的基础上获取了月球初始的内部成分和密度结构,并对后期月幔翻转(Overturn)的程度进行了探讨。