Style of Plate Spreading Derived from the 2008-2014 Velocity Field Across the Northern Volcanic Zone of Iceland

Abstract:

The Northern Volcanic Zone (NVZ) of Iceland is a subaerial part of the divergent boundary between the North-American and Eurasian Plates. At this latitude, the full spreading between the plates is accommodated by the NVZ. We derived the plate boundary velocity field from GPS campaign and continuous measurements between 2008 and 2014, a time period free of any magma intrusion. Average velocities were estimated in the ITRF08 reference frame. The overall extension is consistent with 18 mm/yr in the 104°N direction spreading, in accordance with the MORVEL2010 plate motion model. We find that a 40km-wide band along the plate boundary accommodates about 75% of the full plate velocities. Within this zone, the average strain rate is approximately 0.35 μstrain/yr. The deformation field and the strain rate are, however, much affected by other sources of deformations in the NVZ. These include magmatic sources at the most active volcanic centers, glacial rebound near the ice-caps and geothermal power-plant water extraction. Magmatic sources include a shallow magma chamber deflation under Askja caldera, as well as under Þeistareykir and eventual deep magma inflation north of Krafla volcano. Vatnajökull ice cap melting causes large uplift and outward displacements in the southern part of the NVZ. The two geothermal power-plants near Krafla are inducing local deflations. Our GPS velocities show a 35° change in the direction of the plate boundary axis north of Askja volcano that we infer to be linked to the geometric arrangement of volcanic systems within the NVZ.

We use a simple arctangent model to describe the plate spreading to provide constraints on the location and the locking depth of the spreading axis. For that purpose we divided the area in short overlapping segments having the same amount of GPS points along the plate spreading direction and inverted for the location of the center of the spreading axis and locking depth. With this simple model we can account for most of the plate spreading related deformation in the NVZ. It appears that the locking depth is not uniform all along the length of the plate boundary, with a deeper locking depth in the low activity volcanic systems and a shallower locking depth in the more active volcanic systems of Krafla and Askja.