Part 2: Regional Crustal Structure

Fig 1: Map of iceland showing tectonic features, stations (triangles) of the B array and earthquake epicenters (squares). Note that the epicenter of the earthquake on the Kolbeinsey Ridge (arrow) is off the top of the map. (PostScript version of map).


Fig 2: Low-pass filtered, vertical component record section, reduced to 6.5 km/s, of a Vatnajokumm earthquake observed on the B96 array. Note clear P, PmP, PP, and PPP phases. (PostScript version of plot).


Fig 3: Unfiltered, transverse-horizontal component record section, reduced to 3.8 km/s, of a Vatnajokumm earthquake observed on the B96 array. Note clear S and SmS phases. (PostScript version of plot).


Fig 4: 2-D model of compressional velocity along a NS profile extending from Vatnajokull though Bardardallur to Eyjafjordur, based on traveltimes of P, PmP, PP and PPP phases. The Moho declines from a depth of about 25 km in the north to 31 km in the south. (Color PostScript, B&W PostScript versions of plot, and RTMOD mesh file).


Summary

Design. The B96 array recorded several microearthquakes from source regions north and south of the array (Fig 1). The hypocentral parameters for most of these events have been determined by using arrival times from the SIL and U. Iceland arrays. Record sections for these events often show very clear arrivals of both compressional and shear wave phases. We use these data to determine the regional crustal and upper mantle structure.

The earthquakes fall into four different source regions:

  1. The Kolbeinsey Ridge, 290-360 km north of the B96 array;
  2. Eyjafjordur and Skjalfandi, 70-130 km north of the array;
  3. Vatnajokull, 70-150 km south of the array; and
  4. Myrdalsjokull, 210-280 km south of the array.
Events from the Kolbeinsey Ridge and Myrdalsjokull have clear mantle arrivals (Pn and Sn) which constrain the sub-Moho compressional and shear wave velocities. Events from Eyjafjordur, Skjalfandi and Vatnajokull have clear crustal reverberations (PP, PmP, SS, SmS), which constrain the crustal velocities and Moho deopth.

Analysis Arrival times of the major seismic phases were picked "by eye", using a variety of filter settings (Fig 2 and 3). In general, the P and phases are clearly present on the unfiltered seismograms and were usually picked from them directly. The exception is the S phase from Kolbeinsey Ridge and Myrdalsjokull, which were picked from the low-pass filtered seismograms. The PP, PPP, PmP, SS and SmS phases were typically picked from 0.5-2 Hz low-pass filtered seismograms, although traveltimes were compared to the unfiltered records and adjusted for time-shifts introduced by the filtering.

A 2D compressional velocity model extending from Vatnajokull, north along the B96 array to Eyjafjordur was created using traveltime data from the Vatnajokull and Eyjafjordur earthquakes (Fig 4). The model of shallow structure from the previous section was used to constrain the structure beneath the array, The P traveltime constrain the mid crustal velocities, PmP-P traveltime differences constrain crustal thickness, and PP-P and PPP-P traveltime differences constrain the velocities in the slow, near-surface material. The resulting model fits the data well, except that the P traveltimes from two two neighboring Vatnajokull earthquakes cannot be exactly reconciled (possibly due to errors in hypocentral location). This discrepancy results in a small (0.2s) static anomaly for event (1260156).

Apparent velocities were computed by the least-quares fitting of straignt lines to the linear portions of the traveltime curves (i.e. the short range part of PmP and SmS was omitted). For events in the 70-130 km range, the Pg and Sg rays (with apparent velocities of 6.58-6.86 km/s and 3.78-3.89 km/s, respectively) trurn in the mid-crust. A lower-crustal Pg (with an apparent velocity of 7.34 km/s) is observed from a more distant source (Myrdalsjokull). Clear Pn and Sn waves are also observed from the Myrdalsjokull and Kolbeinsey Ridge regions, with apparent velocities in the 7.66-8.00 km/s and 3.96-4.37 km/s, respectively.

Velocity ratios were then computed by simple division. This technique for determining the compressional to shear velocity ratio is insensitive to uncertainties in source depth and origin time. Ratios were in the range:

region sampledP/S velocity ratiophases
upper crust1.74-1.83Pg and Sg
lower crust1.68-1.79PmP and SmS
uppermost mantle1.85-1.93Pn and Sn

The velocity ratio for the mantle is generally higher than the crust.

Results

  1. Crustal thickness varies from 31 km at the southern end of the model to 25 km at the northen end, with most of the variation occuring at the southern end of the B96 array, at a latitude range of 65:20 - 65:30. Crustal thickness at the northern end is in agreement with estimates previously derived from the EW FIRE profile (see map).
  2. Compressional to shear wave velocity ratios for the lower crust, as deduced from PmP and SmS phase velocities, are in the 1.68-1.79 range typical of normal rock (as contrasted to partial melt).
  3. The near-surface velocities at the southern end of the profile, which are constrained by PP-P differential traveltimes, are significantly lower than those determined from the shot data for central part of the model. This difference is probably due to the removal of unconsolidated surface material in the central part (i.e. Bardadalur valley) by glacial scouring.
  4. Mantle-refracted phases (Pn and Sn) are detected from Myrdalsjokull and the Kolbeinsey Ridge. Both compressional and shear velocities are higher (by about 4%) for Myrdalsjokull (which cross central Iceland) than for Kolbeinsey Ridge (which propagate parallel to that plate boundary). Furthermore, the compressional to shear wave velocity ratio is lower from Kolbeinsey than Myrdalsjokull (1.85 compared to 1.93). Thus the uppermost mantle sampled by the waves from Myrdalsjokull appears relatively cool. Melting in central Iceland must therefore begin several tens of kilometers below Moho.

Supporting data