Magma Dynamics, Crystallisation and Chemical Differentiation of the 1959 Kilauea Iki Lava Lake, Hawaii, revisited

Abstract

Using constraints from an extensive database of geological and geochemical observations along with results from fluid mechanical studies of convection in magma chambers, we identify the main physical processes at work during the solidification of the 1959 Kilauea Iki lava lake. In turn, we investigate their quantitative influence on the crystallization and chemical differentiation of the magma, and on the development of the internal structure of the lava lake. In contrast to previous studies, vigorous stirring in the magma, driven predominately by the descent of dense crystal-laden thermal plumes from the roof solidification front and the ascent of buoyant compositional plumes due to the in situ growth of olivine crystals at the floor, is predicted to have been an inevitable consequence of very strong cooling at the roof and floor. The flow is expected to have caused extensive but imperfect mixing over most of the cooling history of the magma, producing minor compositional stratification at the roof and thermal stratification at the floor. The efficient stirring of the large roof cooling is expected to have resulted in significant internal nucleation of olivine crystals, which ultimately settled to the floor. Additional forcing due to either crystal sedimentation or the ascent of gas bubbles is not expected to have increased significantly the amount of mixing. In addition to convection in the magma, circulation driven by the convection of buoyant interstitial melt in highly permeable crystal-melt mushes forming the roof and the floor of the lava lake is envisaged to have produced a net upward flow of evolved magma from the floor during solidification. In the floor zone, mush convection may have caused the formation of axisymmetric chimneys through which evolved magma drained from deep within the floor into the overlying magma and potentially the roof. We hypothesize that the highly evolved, pipe-like 'vertical olivine-rich bodies' (VORBs) [Helz, 1980 Bull. Volcanol. 43 (1980) 675] observed in the floor zone of the lake are fossil chimneys. In the roof zone, buoyant residual liquid both produced at the roof solidification front and gained from the floor as a result of incomplete convective mixing is envisaged to have percolated or 'leaked' into the overlying highly-permeable cumulate, displacing less buoyant interstitial melt downward. The results from Rayleigh fractionation-type models formulated using boundary conditions based on a quantitative understanding of the convection in the magma indicate that most of the incompatible element variation over the height of the lake can be explained as a consequence of a combination of crystal settling and the extensive but imperfect convective mixing of buoyant residual liquid released from the floor solidification front. The remaining chemical variation is understood in terms of the additional influences of mush convection in the roof and floor on the vertical distribution of incompatible elements. Although cooling was concentrated at the roof of the lake, the floor zone is found to be thicker than the roof zone, implying that it grew more quickly. The large growth rate of the floor is explained as a consequence of a combination of the substantial sedimentation of olivine crystals and more rapid in situ crystallization due to both a higher liquidus temperature and enhanced cooling resulting from imperfect thermal and chemical mixing.