ALBERT

Description: The TAMMAR segment of the Mid-Atlantic Ridge forms a classic propagating system centred about two degrees south of the Kane Fracture Zone. The segment is propagating to the south at a rate of 14 mm yr –1 , 15 per cent faster than the half-spreading rate. Here, we use seismic refraction data across the propagating rift, sheared zone and failed rift to investigate the crustal structure of the system. Inversion of the seismic data agrees remarkably well with crustal thicknesses determined from gravity modelling. We show that the crust is thickened beneath the highly magmatic propagating rift, reaching a maximum thickness of almost 8 km along the seismic line and an inferred (from gravity) thickness of about 9 km at its centre. In contrast, the crust in the sheared zone is mostly 4.5–6.5 km thick, averaging over 1 km thinner than normal oceanic crust, and reaching a minimum thickness of only 3.5 km in its NW corner. Along the seismic line, it reaches a minimum thickness of under 5 km. The PmP reflection beneath the sheared zone and failed rift is very weak or absent, suggesting serpentinisation beneath the Moho, and thus effective transport of water through the sheared zone crust. We ascribe this increased porosity in the sheared zone to extensive fracturing and faulting during deformation. We show that a bookshelf-faulting kinematic model predicts significantly more crustal thinning than is observed, suggesting that an additional mechanism of deformation is required. We therefore propose that deformation is partitioned between bookshelf faulting and simple shear, with no more than 60 per cent taken up by bookshelf faulting.

Description: Estimating small-scale V P / V S variations at depth can be a powerful tool to infer lithology and hydration of a rock, with possible implications for frictional behavior. In principle, from the differential arrival times of P and S phases from a set of spatially clustered earthquakes, an estimate of the local V P / V S can be extracted, because the V P / V S is the scaling factor between the P and S differential times for each pair of earthquakes. We critically review the technique proposed by Lin and Shearer (2007) , in which the mean value over all stations is subtracted from the differential arrival times of each pair of events in order to make the method independent of a priori information on origin times. The demeaned differential P and S arrival times are plotted on a plane, and the V P / V S ratio is estimated by fitting the points on this plane. We tested the method by both theoretical analysis and numerical simulations of P and S travel times in several velocity models. We found that the method returns exact values of V P / V S only in the case of a medium with homogeneous V P / V S , whereas, when a V P / V S gradient is present, the estimates are biased as an effect of systematic differences between P and S takeoff angles. We demonstrated that this bias arises from the demeaning of the arrival times over the stations. In layered models with V P / V S decreasing with depth, we found that V P / V S is overestimated or underestimated, respectively, for takeoff angles larger or smaller than 90°. Moreover, we calculated analytically the dependence of this bias on the takeoff angles. Our simulations also showed that the difference between the calculated and the expected V P / V S is reduced for simple horizontally layered velocity structures (〈0.06), whereas it is 0.27 in a more realistic velocity model mimicking a subduction zone.

Description: The very shallow part of subduction megathrusts occasionally hosts tsunami earthquakes, with unusually slow rupture propagation. The aftershock sequence of the 2010 M w  8.8 Maule earthquake, offshore Chile, provides us with the opportunity to study systematic changes in source properties for smaller earthquakes within a single segment of a subduction zone. We invert amplitude spectra for double-couple moment tensors and centroid depths of 71 aftershocks of the Maule earthquake down to magnitudes M w  4.0. In addition, we also derive average source durations. We find that shallower earthquakes tend to have longer normalized source durations on average, similar to the pattern observed previously for larger magnitude events. This depth dependence is observable for thrust and normal earthquakes. The normalized source durations of normal-faulting earthquakes are at the lower end of those for thrust earthquakes, probably because of the higher stress drops of intraplate earthquakes compared to interplate earthquakes. We suggest from the similarity of the depth dependence of normal and thrust events and between smaller and larger magnitude earthquakes that the depth-dependent variation of rigidity, rather than frictional conditional stability at the plate interface, is primarily responsible for the observed pattern. Tsunami earthquakes probably require both low rigidity and conditionally stable frictional conditions; the presence of long-duration moderate-magnitude events is therefore a helpful but not sufficient indicator for areas at risk of tsunami earthquakes. Online Material: Table of earthquake inversion results for the 2010 M w  8.8 Maule, Chile, aftershock sequence.

Description: This study presents P - and S -wave velocity variations for the upper mantle in southern Scandinavia and northern Germany based on teleseismic traveltime tomography. Tectonically, this region includes the entire northern part of the prominent Tornquist Zone which follows along the transition from old Precambrian shield units to the east to younger Phanerozoic deep sedimentary basins to the southwest. We combine data from several separate temporary arrays/profiles (276 stations) deployed over a period of about 15 yr and permanent networks (31 stations) covering the areas of Denmark, northern Germany, southern Sweden and southern Norway. By performing an integrated P - and S -traveltime analysis, we obtain the first high-resolution combined 3-D V P and V S models, including variations in the V P / V S ratio, for the whole of this region of study. Relative station mean traveltime residuals vary within ±1 s for P wave and ±2 s for S wave, with early arrivals in shield areas of southern Sweden and later arrivals in the Danish and North German Basins, as well as in most of southern Norway. In good accordance with previous, mainly P -velocity models, a marked upper-mantle velocity boundary (UMVB) is accurately delineated between shield areas (with high seismic mantle velocity) and basins (with lower velocity). It continues northwards into southern Norway near the Oslo Graben area and further north across the Southern Scandes Mountains. This main boundary, extending to a depth of at least 300 km, is even more pronounced in our new S -velocity model, with velocity contrasts of up to ±2–3 per cent. It is also clearly reflected in the V P / V S ratio. Differences in this ratio of up to about ±2 per cent are observed across the boundary, with generally low values in shield areas to the east and relatively higher values in basin areas to the southwest and in most of southern Norway. Differences in the V P / V S ratio are believed to be a rather robust indicator of upper-mantle compositional differences. For the depth interval of about 100–300 km, thick, depleted, relatively cold shield lithosphere is indicated in southern Sweden, contrasting with more fertile, warm mantle asthenosphere beneath most of the basins in Denmark and northern Germany. Both compositional and temperature differences seem to play a significant role in explaining the UMVB between southern Norway and southern Sweden. In addition to the main regional upper-mantle velocity contrasts, a number of more local anomaly features are also outlined and discussed.

Description: To investigate temporal seismic velocity changes due to earthquake related processes and environmental forcing in Northern Chile, we analyse 8 yr of ambient seismic noise recorded by the Integrated Plate Boundary Observatory Chile (IPOC). By autocorrelating the ambient seismic noise field measured on the vertical components, approximations of the Green's functions are retrieved and velocity changes are measured with Coda Wave Interferometry. At station PATCX, we observe seasonal changes in seismic velocity caused by thermal stress as well as transient velocity reductions in the frequency range of 4–6 Hz. Sudden velocity drops occur at the time of mostly earthquake-induced ground shaking and recover over a variable period of time. We present an empirical model that describes the seismic velocity variations based on continuous observations of the local ground acceleration. The model assumes that not only the shaking of large earthquakes causes velocity drops, but any small vibrations continuously induce minor velocity variations that are immediately compensated by healing in the steady state. We show that the shaking effect is accumulated over time and best described by the integrated envelope of the ground acceleration over the discretization interval of the velocity measurements, which is one day. In our model, the amplitude of the velocity reduction as well as the recovery time are proportional to the size of the excitation. This model with two free scaling parameters fits the data of the shaking induced velocity variation in remarkable detail. Additionally, a linear trend is observed that might be related to a recovery process from one or more earthquakes before our measurement period. A clear relationship between ground shaking and induced velocity reductions is not visible at other stations. We attribute the outstanding sensitivity of PATCX to ground shaking and thermal stress to the special geological setting of the station, where the subsurface material consists of relatively loose conglomerate with high pore volume leading to a stronger nonlinearity compared to the other IPOC stations.

Description: Global earthquake locations are often associated with very large systematic travel-time residuals even for clear arrivals, especially for regional and near-regional stations in subduction zones because of their strongly heterogeneous velocity structure. Travel-time corrections can drastically reduce travel-time residuals at regional stations and, in consequence, improve the relative location accuracy. We have extended the shrinking-box source-specific station terms technique to regional and teleseismic distances and adopted the algorithm for probabilistic, nonlinear, global-search location. We evaluated the potential of the method to compute precise relative hypocentre locations on a global scale. The method has been applied to two specific test regions using existing P - and pP -phase picks. The first data set consists of 3103 events along the Chilean margin and the second one comprises 1680 earthquakes in the Tonga-Fiji subduction zone. Pick data were obtained from the GEOFON earthquake bulletin, produced using data from all available, global station networks. A set of timing corrections varying as a function of source position was calculated for each seismic station. In this way, we could correct the systematic errors introduced into the locations by the inaccuracies in the assumed velocity structure without explicitly solving for a velocity model. Residual statistics show that the median absolute deviation of the travel-time residuals is reduced by 40–60 per cent at regional distances, where the velocity anomalies are strong. Moreover, the spread of the travel-time residuals decreased by ~20 per cent at teleseismic distances (〉28°). Furthermore, strong variations in initial residuals as a function of recording distance are smoothed out in the final residuals. The relocated catalogues exhibit less scattered locations in depth and sharper images of the seismicity associated with the subducting slabs. Comparison with a high-resolution local catalogue reveals that our relocation process significantly improves the hypocentre locations compared to standard locations.

Description: We infer seismic azimuthal anisotropy from ambient-noise-derived Rayleigh waves in the wider Vienna Basin region. Cross-correlations of the ambient seismic field are computed for 1953 station pairs and periods from 5 to 25? s to measure the directional dependence of interstation Rayleigh-wave group velocities. We perform the analysis for each period on the whole data set, as well as in overlapping 2°-cells to regionalize the measurements, to study expected effects from isotropic structure, and isotropic–anisotropic trade-offs. To extract azimuthal anisotropy that relates to the anisotropic structure of the Earth, we analyse the group velocity residuals after isotropic inversion. The periods discussed in this study (5–20? s) are sensitive to crustal structure, and they allow us to gain insight into two distinct mechanisms that result in fast orientations. At shallow crustal depths, fast orientations in the Eastern Alps are S/N to SSW/NNE, roughly normal to the Alps. This effect is most likely due to the formation of cracks aligned with the present-day stress-field. At greater depths, fast orientations rotate towards NE, almost parallel to the major fault systems that accommodated the lateral extrusion of blocks in the Miocene. This is coherent with the alignment of crystal grains during crustal deformation occurring along the fault systems and the lateral extrusion of the central part of the Eastern Alps.

Description: We analyze S-receiver functions to investigate the variations of lithospheric thickness below the entire region of Turkey and surroundings. The teleseismic data used here have been compiled combining all permanent seismic stations which are open to public access. We obtained almost 12 000 S-receiver function traces characterizing the seismic discontinuities between the Moho and the discontinuity at 410 km depth. Common-conversion-points stacks yield well-constrained images of the Moho and of the lithosphere–asthenosphere boundary (LAB). Results from previous studies suggesting shallow LAB depths between 80 and 100 km are confirmed in the entire region outside the subduction zones. We did not observe changes of LAB depths across the North and East Anatolian Faults. To the east of Cyprus, we see indications of the Arabian LAB. The African plate is observed down to about 150 km depth subducting to the north and east between the Aegean and Cyprus with a tear at Cyprus. We also observed the discontinuity at 410 km depth and a negative discontinuity above the 410, which might indicate a zone of partial melt above this discontinuity.

Description: We analyze S-receiver functions to investigate variations of lithospheric thickness below the entire region of Turkey and surrounding areas. The teleseismic data used here have been compiled combining all permanent seismic stations which are open to public access. We obtained almost 12 000 S-receiver function traces characterizing the seismic discontinuities between the Moho and the discontinuity at 410 km depth. Common-conversion-point stacks yield well-constrained images of the Moho and of the lithosphere–asthenosphere boundary (LAB). Results from previous studies suggesting shallow LAB depths between 80 and 100 km are confirmed in the entire region outside the subduction zones. We did not observe changes in LAB depths across the North and East Anatolian faults. To the east of Cyprus, we see indications of the Arabian LAB. The African plate is observed down to about 150 km depth subducting to the north and east between the Aegean and Cyprus with a tear at Cyprus. We also observed the discontinuity at 410 km depth and a negative discontinuity above the 410, which might indicate a zone of partial melt above this discontinuity.