Focusing patterns of seismicity with relocation and collapsing

Seismicity is generally concentrated on faults or in fault zones of varying, sometimes complex geometry. An earthquake catalog, compiled over time, contains useful information about this geometry, which can help understanding the tectonics of a region. Interpreting the geometrical distribution of events in a catalog is often complicated by the diffuseness of the earthquake locations. Here, we explore a number of strategies to reduce this diffuseness and hence simplify the seismicity pattern of an earthquake catalog. These strategies utilize information about event locations contained in their overall catalog distribution. They apply this distribution as an a priori constraint on relocations of the events, or as an attractor for each individual event in a collapsing scheme, and thereby focus the locations. The latter strategy is not a relocation strategy in a strict sense, although event foci are moved, because the movements are not driven by data misfit. Both strategies simplify the seismicity pattern of the catalog and may help to interpret it. A synthetic example and a real-data example from an aftershock sequence in south west Iceland are presented to demonstrate application of the strategies. Entropy is used to quantify their effect.     

Ellipticity corrections for seismic phases

The advent of broad-band seismology has meant that use is being made of a wide range of seismic phases, for many of which ellipticity corrections have not been readily available. In particular, when many seismic phases are used in location schemes, it is important that the systematic effects of ellipticity are included for each phase.
An efficient and effective procedure for constructing ellipticity corrections is to make use of the ray-based approach of Dziewonksi & Gilbert (1976), as reformulated by Doornbos (1988), in conjunction with the rapid evaluation of traveltimes and slownesses for a given range using the tauspline procedure of Buland & Chapman (1983).
Ellipticity coefficients have been tabulated for a wide range of seismic phases and are available in electronic form. The ellipticity correction procedures have been extended to include an allowance for diffraction phenomena, for example P _diff, S _diff diffracted along the core-mantle boundary. Corrections for additional phases can be generated by building the ellipticity coefficients from suitable combinations of the coefficients for different phase segments.     

Emplacement and arrest of sheets and dykes in central volcanoes

Sheet intrusions are of two main types: local inclined (cone) sheets and regional dykes. In Iceland, the inclined sheets form dense swarms of (mostly) basaltic, 0.5–1 m thick sheets, dipping either at 20–50° or at 75–90° towards the central volcano to which they belong. The regional dykes are (mostly) basaltic, 4–6 m thick, subvertical, subparallel and form swarms, less dense than those of the sheets but tens of kilometres long, in the parts of the volcanic systems that are outside the central volcanoes. In both types of swarms, the intrusion intensity decreases with altitude in the lava pile. Theoretical models generally indicate very high crack-tip stresses for propagating dykes and sheets. Nevertheless, most of these intrusions become arrested at various crustal depths and never reach the surface to supply magma to volcanic eruptions. Two principal mechanisms are proposed to explain arrest of dykes and sheets. One is the generation of stress barriers, that is, layers with local stresses unfavourable for the intrusion propagation. The other is mechanical anisotropy whereby sheet intrusions become arrested at discontinuities. Stress barriers may develop in several ways. First, analytical solutions for a homogeneous and isotropic crust show that the intensity of the tensile stress associated with a pressured magma chamber falls off rapidly with distance from the chamber. Thus, while dyke and sheet injection in the vicinity of a chamber may be favoured, dyke and sheet arrest is encouraged in layers (stress barriers) at a certain distance from the chamber. Second, boundary-element models for magma chambers in a mechanically layered crust indicate abrupt changes in tensile stresses between layers of contrasting Young’s moduli (stiffnesses). Thus, where soft pyroclastic layers alternate with stiff lava flows, as in many volcanoes, sheet and dyke arrest is encouraged. Abrupt changes in stiffness between layers are commonly associated with weak and partly open contacts and other discontinuities. It follows that stress barriers and discontinuities commonly operate together as mechanisms of dyke and sheet arrest in central volcanoes.     

Geochemical stratigraphy of the lavas of the Solund-Stavfjord Ophiolite Complex,W. Norway,and magma-chamber dynamics

Bulletin of Volcanology - The stratigraphy and geochemistry of the uppermost 200–300&;nbsp;m of the metabasalt sequence of the Solund-Stavfjord Ophiolite Complex of western Norway has been...     

Rayleigh wave tomography in the North Atlantic: high resolution images of the Iceland, Azores and Eifel mantle plumes

Presented in this paper is a high resolution S_v-wave velocity and azimuthal anisotropy model for the upper mantle beneath the North Atlantic and surrounding region derived from the analysis of about 9000 fundamental and higher-mode Rayleigh waveforms. Much of the dataset comes from global and national digital seismic networks, but to improve the path coverage a number of instruments at coastal sites in northwest Europe, Iceland and eastern Greenland was deployed by us and a number of collaborators. The dense path coverage, the wide azimuthal distribution and the substantial higher-mode content of the dataset, as well as the relatively short path-lengths in the dataset have enabled us to build an upper mantle model with a horizontal resolution of a few hundred kilometers extending to 400 km depth. Low upper mantle velocities exist beneath three major hotspots: Iceland, the Azores and Eifel. The best depth resolution in the model occurs in NW Europe and in this area low S_v-velocities in the vicinity of the Eifel hotspot extend to about 400 km depth. Major negative velocity anomalies exist in the North Atlantic upper mantle beneath both Iceland and the Azores hotspots. Both anomalies are, above 200 km depth, 4–7% slow with respect to PREM and elongated along the mid-Atlantic Ridge. Low velocities extend to the south of Iceland beneath the Reykjanes Ridge where other geophysical and geochemical observations indicate the presence of hot plume material. The low velocities also extend beneath the Kolbeinsey Ridge north of Iceland, where there is also supporting geochemical evidence for the presence of hot plume material. The low-velocity upper mantle beneath the Kolbeinsey Ridge may also be associated with a plume beneath Jan Mayen. The anomaly associated with the Azores extends from about 25°N to 45°N along the ridge axis, which is in agreement with the area influenced by the Azores Plume, predicted from geophysical and geochemical observations. Compared to the anomaly associated with Iceland, the Azores anomaly is elongated further along the ridge, is shallower and decays more rapidly with depth. The fast propagation direction of horizontally propagating S_v-waves in the Atlantic south of Iceland correlates well with the east–west ridge-spreading direction at all depths and changes to a direction close to NS in the vicinity of Iceland.     

Structural analysis of a transform fault-rift zone junction in North Iceland

On the north coast of Iceland, the rift zone in North Iceland is shifted about 120 km to the west where it meets with, and joins, the mid-ocean Kolbeinsey ridge. This shift occurs along the Tjörnes fracture zone, an 80-km-wide zone of high seismicity, which is an oblique (non-perpendicular) transform fault. There are two main seismic lineaments within the Tjörnes fracture zone, one of which continues on land as a 25-km-long WNW-trending strike-slip fault. This fault, referred to as the Husavik fault, meets with, and joins, north-trending normal faults of the Theistareykir fissure swarm in the axial rift zone. The most clear-cut of these junctions occurs in a basaltic pahoehoe lava flow, of Holocene age, where the Husavik fault joins a large normal fault called Gudfinnugja. At this junction, the Husavik fault strikes N55°W, whereas Gudfinnugja strikes N5°E, so that they meet at an angle of 60°. The direction of the spreading vector in North Iceland is about N73°W, which is neither parallel with the strike of the Husavik fault nor perpendicular to the strike of the Gudfinnugja fault. During rifting episodes there is thus a slight opening on the Husavik fault as well as a considerable dextral strike-slip movement along the Gudfinnugja fault. Consequently, in the Holocene lava flow, there are tension fractures, collapse structures and pressure ridges along the Husavik fault, and pressure ridges and dextral pull-apart structures subparallel with the Gudfinnugja fault. The 60° angle between the Husavik strike-slip fault and the Gudfinnugja normal fault is the same as the angle between the Tjörnes fracture zone transform fault and the adjacent axial rift zones of North Iceland and the Kolbeinsey ridge. The junction between the faults of Husavik and Gudfinnugja may thus be viewed as a smaller-scale analogy to the junction between this transform fault and the nearby ridge segments. Using the results of photoelastic and finite-element studies, a model is provided for the tectonic development of these junctions. The model is based on an analogy between two offset cuts (mode I fractures) loaded in tension and segments of the axial rift zones (or parts thereof in the case of the Husavik fault). The results indicate that the Tjörnes fracture zone in general and the Husavik fault in particular, developed along zones of maximum shear stress. Furthermore, the model suggests that, as the ridge-segments propagate towards a zero-underlapping configuration, the angle between them and the associated major strike-slip faults gradually increases. This conclusion is supported by the trends of the main seismic lineaments of the Tjörnes fracture zone.     

Stress estimates from the length/width ratios of fractures

The hypothesis is advanced that, provided Young's modulus and Poisson's ratio of the rock are known, the length/width ratios of tension fractures can be used to estimate the tensile stress (assumed constant along the length of each fracture) at the time of fracture formation. The hypothesis is tested on a fissure swarm in a 10,000 year-old basaltic lava in Iceland. The length/width ratios of the fissures give the average tensile stress as of the order of a few MPa.