Lunar floor-fractured craters: Modes of dike and sill emplacement and implications of gas production and intrusion cooling on surface morphology and structure

We present a temporal sequence of 1) dike emplacement, 2) sill formation and surface deformation, 3) bubble rise, foam layer formation and collapse, 4) intrusion cooling, and a synthesis of predicted deformation sequence and eruption styles. Initial lateral injection of the sill at a depth well below the upper dike tip initiates upbowing of the overburden, leveraging deformation of the crater floor melt sheet above. This is followed by lateral spreading of the sill toward the edges of the crater floor, where crater wall and rim deposit overburden inhibit further lateral growth, and the sill grows vertically into a laccolith or bysmalith, uplifting the entire floor above the intrusion. Subsidiary dikes can be emplaced in the fractures at the uplift margins and will rise to the isostatic level of the initial dike tip; if these contain sufficient volatiles to decrease magma density, eruptions can also occur. This initial phase of intrusion, sill lateral spreading and floor uplift occurs within a few hours after initial dike emplacement. During the subsequent cooling of the sill, bubbles can rise hundreds of meters to the top of the intrusion to create a foam layer; when drainage of gas bubble wall magma occurs in the foam layer, a continuous gas layer forms above the foam. Gas formation and upward migration produces an increase in sill thickness, while subsequent cooling and solidification cause a thickness decreases and subsidence. The total topographic evolution history, following an initial 2 km thick sill intrusion and floor uplift (hours), includes further floor uplift by gas formation and migration (decades; ∼30 m), followed by cooling, solidification and subsidence (∼a century; ∼350 m). An initial 2 km thick sill is predicted to have a final thickness of ∼1.7 km. This predicted sequence of events can be compared with the sequence of floor deformation and volcanism in FFCs in order to test and refine this model.