Fracture and surface crust development in a Holocene pahoehoe lava flow on the Island of Tenerife, Canaries

Abstract An almost horizontal pahoehoe surface in a Holocene plagioclase basanite lava on Tenerife displays three scales of fracture within the surface crust. An early-formed set of large-scale fractures divides up the surface into an orthogonal set of rectangular slabs with dimensions of several metres and depths of 10–12 cm. The shortest slab dimension is parallel to the flow emplacement direction, inferred from a strong surface lineation. The slabs are domed with the centre an average of 9.6 cm (with range 4–19.6 cm) above the edges of the slabs. Profiles of the slabs normal and perpendicular to the margins and through the crest indicate that they can be described by a power law in which the deflection of the slab, h , is related to the distance from the crest, x , with an exponent between 2 and 3. Analysis of joints within the slabs indicates two smaller scale networks. An intermediate scale joint network bounds blocks with rectilinear to polygonal shapes in plan-view and has a characteristic mean spacing of 24.2 cm (range 10.5–48 cm). The major fractures in this set are normal and parallel to the slab margins. A smaller-scale joint network bounds polygonal equant blocks in plan-view and has characteristic spacing of 6.4 cm (range 3.7–10.5 cm). A model of cooling from the pahoehoe surface is used to constrain the growth of the crust and timing of fracture development. The large-scale slabs are attributed to localised accumulation of gas beneath the growing crust causing buoyant forces. The tensile stresses caused by uplift are sufficient to form the large-scale fractures after 2 or 3 h of cooling. The intermediate scale fracture network is attributed to the flexure of the slab crust. The smaller scale polygonal joint network is related to the build up of isotropic tensile stresses in the cooling slab crust due to thermal contraction with fracture development being promoted by the flexure of the slabs. An analysis of the slab deformation indicates that lava crust is weak. The weakness is explained by division of the crust into three zones: an outer zone with small scale joints that cause negligible strength, a middle zone of elastic behaviour in which stress can build up, and a lower zone of plastic deformation. The crustal slabs display profiles similar to that expected in a bending elastic plate. The deformation of the 10-cm-thick crust can be explained if the elastic zone was about 2-cm thick. This result agrees with an independent calculation of elastic zone thickness based on the position of the brittle–ductile transition being located at the 600°C isotherm at a depth of about 2 cm when the crustal slabs were rifted apart.

[1]  Laszlo P. Keszthelyi,et al.  Calculated effect of vesicles on the thermal properties of cooling basaltic lava flows , 1994 .

[2]  P. Roeder,et al.  Olivine-liquid equilibrium , 1970 .

[3]  G. Walker,et al.  Structure, and origin by injection of lava under surface crust, of tumuli, “lava rises”, “lava-rise pits”, and “lava-inflation clefts” in Hawaii , 1991 .

[4]  C. W. Gear,et al.  Numerical initial value problem~ in ordinary differential eqttations , 1971 .

[5]  P. C. Grew,et al.  Stress corrosion theory of crack propagation with applications to geophysics , 1977 .

[6]  R. Sparks,et al.  The 1975 sub-terminal lavas, mount etna: a case history of the formation of a compound lava field , 1976 .

[7]  D. Neuville,et al.  Thermodynamic and rheological properties of rhyolite and andesite melts , 1993 .

[8]  H. Huppert,et al.  Convection and crystallization in magma cooled from above , 1990 .

[9]  D. Dingwell Recent experimental progress in the physical description of silicic magma relevant to explosive volcanism , 1998, Geological Society, London, Special Publications.

[10]  D. Dingwell,et al.  Structural relaxation in silicate melts and non-Newtonian melt rheology in geologic processes , 1989 .

[11]  Richard J Goldstein,et al.  Turbulent convection in a horizontal layer of water , 1973, Journal of Fluid Mechanics.

[12]  M. P. Ryan,et al.  The viscosity of synthetic and natural silicate melts and glasses at high temperatures and 1 bar (105 pascals) pressure and at higher pressures , 1986 .

[13]  A. Provost,et al.  Rates of solidification of Apollo 11 basalt and Hawaiian tholeiite , 1972 .

[14]  Paul Budkewitsch,et al.  Modelling the evolution of columnar joints , 1994 .

[15]  P. L. Hancock,et al.  Brittle microtectonics: principles and practice , 1985 .

[16]  D. L. Peck,et al.  The Formation of Columnar Joints in the Upper Part of Kilauean Lava Lakes, Hawaii , 1968 .

[17]  J. P. Kauahikaua,et al.  Emplacement and inflation of pahoehoe sheet flows: observations and measurements of active lava flows on Kilauea volcano, Hawaii , 1994 .

[18]  D. L. Peck Cooling and vesiculation of Alae lava lake, Hawaii , 1978 .

[19]  A. McBirney,et al.  Properties of some common igneous rocks and their melts at high temperatures , 1973 .

[20]  J. Martí,et al.  Basanite-phonolite lineages of the Teide-Pico Viejo Volcanic Complex, Tenerife, Canary Islands , 1998 .

[21]  A. Aydin,et al.  Effect of thermal regime on growth increment and spacing of contraction joints in basaltic lava , 1993 .