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: The Maule earthquake (2010 February 27, M w 8.8, Chile) broke the subduction megathrust along a previously locked segment. Based on an international aftershock deployment, catalogues of precisely located aftershocks have become available. Using 23 well-located aftershocks, we calibrate the classic teleseismic backprojection procedure to map the high-frequency seismic radiation emitted during the earthquake. The calibration corrects traveltimes in a standard earth model both with a static term specific to each station, and a ‘dynamic’ term specific to each combination of grid point and station. The second term has been interpolated over the whole slipping area by kriging, and is about an order of magnitude smaller than the static term. This procedure ensures that the teleseismic images of rupture development are properly located with respect to aftershocks recorded with local networks and does not depend on accurate hypocentre location of the main shock. We track a bilateral rupture propagation lasting ~160 s, with its dominant branch rupturing northeastwards at about 3 km s –1 . The area of maximum energy emission is offset from the maximum coseismic slip but matches the zone where most plate interface aftershocks occur. Along dip, energy is preferentially released from two disconnected interface belts, and a distinct jump from the shallower belt to the deeper one is visible after about 20 s from the onset. However, both belts keep on being active until the end of the rupture. These belts approximately match the position of the interface aftershocks, which are split into two clusters of events at different depths, thus suggesting the existence of a repeated transition from stick-slip to creeping frictional regime.

Description: We focus on the relation between seismic and total postseismic afterslip following the Maule M w 8.8 earthquake on 2010 February 27 in central Chile. First, we calculate the cumulative slip released by aftershock seismicity. We do this by summing up the aftershock regions and slip estimated from scaling relations. Comparing the cumulative seismic slip with afterslip models we show that seismic slip of individual aftershocks exceeds locally the inverted afterslip model from geodetic constraints. As the afterslip model implicitly contains the displacements from the aftershocks, this reflects the tendency of afterslip models to smear out the actual slip pattern. However, it also suggests that locally slip for a number of the larger aftershocks exceeds the aseismic slip in spite of the fact that the total equivalent moment of the afterslip exceeds the cumulative moment of aftershocks by a large factor. This effect, seen weakly for the Maule 2010 and also for the Tohoku 2011 earthquake, can be explained by taking into account the uncertainties of the seismicity and afterslip models. In spite of uncertainties, the hypocentral region of the Nias 2005 earthquake is suggested to release a large fraction of moment almost purely seismically. Therefore, these aftershocks are not driven solely by the afterslip but instead their slip areas have probably been stressed by interseismic loading and the mainshock rupture. In a second step, we divide the megathrust of the Maule 2010 rupture into discrete cells and count the number of aftershocks that occur within 50 km of the centre of each cell as a function of time. We then compare this number to a time-dependent afterslip model by defining the ‘afterslip to aftershock ratio’ (ASAR) for each cell as the slope of the best fitting line when the afterslip at time t is plotted against aftershock count. Although we find a linear relation between afterslip and aftershocks for most cells, there is significant variability in ASAR in both the downdip and along-strike directions of the megathrust. We compare the spatial distribution of ASAR with the spatial distribution of seismic coupling, coseismic slip and Bouguer gravity anomaly, and in each case we find no significant correlation.

Description: Earthquakes are commonly located by linearized inversion of discrete arrival time picks made from signals recorded at a network of seismic stations. If mis-picks are made, these will contribute to the location, therefore causing potential bias. For data recorded by a dense seismic array, direct imaging methods can be applied instead. We describe the ‘coalescence microseismic mapping’ method, which is a bridge between the two approaches and will operate with seismic data recorded continuously on a sparse array. By continuously mapping scalar signals derived from the envelope of seismic arrivals we derive robust estimates of the spatio-temporal coordinates of the origins of seismic events. Noisy data are migrated away from the correct origin, so do not contribute to errors in location. The method is rooted in a Bayesian formulation of event location traveltime inversion, allows imaging of source locations and has the capacity to handle errors in modelled traveltimes. It has the advantage of working with any 3-D velocity model, which therefore may include anisotropy. It also automatically incorporates both P - and S -wave data. A multiresolution grid search leads to an efficient implementation, with a search over a larger domain including joint inversion for location and velocity structure possible where warranted by the data quality. We discuss the theory and implementation of this method and illustrate it with real data from microseismic events in Iceland caused by melt intrusion in the crust.

Description: The northeastern margin of the Tibetan Plateau, which includes the Qiangtang and Songpan-Ganzi terranes as well as the Kunlun Shan and the Qaidam Basin, continues to deform in response to the ongoing India–Eurasia collision. To test competing hypotheses concerning the mechanisms for this deformation, we assembled a high-quality data set of approximately 14 000 P - and 4000 S -wave arrival times from earthquakes at teleseismic distances from the International Deep Profiling of Tibet and the Himalaya, Phase IV broad-band seismometer deployments. We analyse these arrival times to determine tomographic images of P - and S -wave velocities in the upper mantle beneath this part of the plateau. To account for the effects of major heterogeneity in crustal and uppermost mantle wave velocities in Tibet, we use recent surface wave models to construct a starting model for our teleseismic body wave inversion. We compare the results from our model with those from simpler starting models, and find that while the reduction in residuals and results for deep structure are similar between models, the results for shallow structure are different. Checkerboard tests indicate that features of ~125 km length scale are reliably imaged throughout the study region. Using synthetic tests, we show that the best recovery is below ~300 km, and that broad variations in shallow structure can also be recovered. We also find that significant smearing can occur, especially at the edges of the model. We observe a shallow dipping seismically fast structure at depths of ~140–240 km, which dies out gradually between 33°N and 35°N. Based on the lateral continuity of this structure (from the surface waves) we interpret it as Indian lithosphere. Alternatively, the entire area could be thickened by pure shear, or the northern part could be an underthrust Lhasa Terrane lithospheric slab with only the southern part from India. We see a deep fast wave velocity anomaly (below 300 km), that is consistent with receiver function observations of a thickened transition zone and could be a fragment of oceanic lithosphere. In NE Tibet, it appears to be disconnected from faster wave velocities above (i.e. it is not downwelling or subducting here). Our models corroborate results of previous work which imaged a relatively slow wave velocity region below the Kunlun Shan and northern Songpan-Ganzi Terrane, which is difficult to reconcile with the hypothesis of southward-directed continental subduction at the northern margin. Wave velocities in the shallow mantle beneath the Qaidam Basin are faster than normal, and more so in the east than the west.

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: In probabilistic Bayesian inversions, data uncertainty is a crucial parameter for quantifying the uncertainties and correlations of the resulting model parameters or, in transdimensional approaches, even the complexity of the model. However, in many geophysical inference problems it is poorly known. Therefore, it is common practice to allow the data uncertainty itself to be a parameter to be determined. Although in principle any arbitrary uncertainty distribution can be assumed, Gaussian distributions whose standard deviation is then the unknown parameter to be estimated are the usual choice. In this special case, the paper demonstrates that a simple analytical integration is sufficient to marginalise out this uncertainty parameter, reducing the complexity of the model space without compromising the accuracy of the posterior model probability distribution. However, it is well known that the distribution of geophysical measurement errors, although superficially similar to a Gaussian distribution, typically contains more frequent samples along the tail of the distribution, so-called outliers. In linearized inversions these are often removed in subsequent iterations based on some threshold criterion, but in Markov chain Monte Carlo (McMC) inversions this approach is not possible as they rely on the likelihood ratios, which cannot be formed if the number of data points varies between the steps of the Markov chain. The flexibility to define the data error probability distribution in McMC can be exploited in order to account for this pattern of uncertainties in a natural way, without having to make arbitrary choices regarding residual thresholds. In particular, we can regard the data uncertainty distribution as a mixture between a Gaussian distribution, which represent valid measurements with some measurement error, and a uniform distribution, which represents invalid measurements. The relative balance between them is an unknown parameter to be estimated alongside the standard deviation of the Gauss distribution. For each data point, the algorithm can then assign a probability to be an outlier, and the influence of each data point will be effectively downgraded according to its probability to be an outlier. Furthermore, this assignment can change as the McMC search is exploring different parts of the model space. The approach is demonstrated with both synthetic and real tomography examples. In a synthetic test, the proposed mixed measurement error distribution allows recovery of the underlying model even in the presence of 6 per cent outliers, which completely destroy the ability of a regular McMC or linear search to provide a meaningful image. Applied to an actual ambient noise tomography study based on automatically picked dispersion curves, the resulting model is shown to be much more consistent for different data sets, which differ in the applied quality criteria, while retaining the ability to recover strong anomalies in selected parts of the model.

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.