Estimation of the viscosities of silicate liquids at high pressure from measurements of oxygen diffusivity

The dynamics and evolution of a magma ocean depend significantly on the structure and physical properties of silicate liquids at high pressure. Viscosity is particularly important because of its effect on the dynamics of convection and crystal setting. Traditionally, the viscosities of silicate liquids have been measured at high pressure by falling sphere viscometry but, due to some limitations of this technique, few measurements have been made at pressures exceeding 2.5 GPa. An alternative approach, which has previously been proven up to 2 GPa, is to use high pressure measurements of oxygen self-diffusivity (D) to estimate viscosity (eta) from the Eyring relationship eta = kT/D lambda (where k is the Boltzmann constant, T is absolute temperature, and lambda is the diffusive jump distance). We have performed preliminary experiments on Na2Si4O9 liquid in order to explore the feasibility of this method up to pressures in excess of 10 GPa. Diffusion couples were prepared from glass starting materials with one half of the sample enriched in O-18 and the other half containing the natural abundance (0.2 wt percent). Diffusion experiments were performed at 1600 to 1800 C at pressures in the range of 2.5 to 10 GPa for times up to 6 minutes using a 1200 ton multianvil apparatus. Oxygen diffusivities were derived from the resulting O-18 concentration profiles, which were analyzed using an ion microprobe. At 1800 C, oxygen diffusivities increase continuously from 1 x 10(exp -10) m(exp -2) s(exp -1) at 2.5 Gpa to 5 x 10(exp -10) m(exp -2) s(exp -1) at 10 GPa. From the Eyring relationship, the viscosity of Na2SiO9 liquid is predicted to decrease by approx. 0.7 log units over this pressure range at 1800 C. These trends, which can be related to changes in the structure of the liquid at high pressure, are in agreement with results of molecular dynamics calculations. These results suggest that it should be possible to estimate the viscosities of silicate liquids up to at least 15 GPa and 2200 C by this method.