ALBERT

Description: We use traveltime data of local earthquakes and controlled sources observed by a large, temporary, amphibious seismic network to reveal the anatomy of the southcentral Chilean subduction zone (37–39°S) between the trench and the magmatic arc. At this location the giant 1960 earthquake (M = 9.5) nucleated and ruptured almost 1000 km of the subduction megathrust. For the three-dimensional tomographic inversion we used 17,148 P wave and 10,049 S wave arrival time readings from 439 local earthquakes and 94 shots. The resolution of the tomographic images was explored by analyzing the model resolution matrix and conducting extensive numerical tests. The downgoing lithosphere is delineated by high seismic P wave velocities. High v p/v s ratio in the subducting slab reflects hydrated oceanic crust and serpentinized uppermost oceanic mantle. The subducting oceanic crust can be traced down to a depth of 80 km, as indicated by a low velocity channel. The continental crust extends to approximately a 50-km depth near the intersection with the subducting plate. This suggests a wide contact zone between continental and oceanic crust of about 150 km, potentially supporting the development of large asperities. Eastward the crustal thickness decreases again to a minimum of about a 30-km depth. Relatively low v p/v s at the base of the forearc does not support a large-scale serpentinization of the mantle wedge. Offshore, low v p and high v p/v s reflect young, fluid-saturated sediments of forearc basins and the accretionary prism.

Description: Providing quantitative microzonation results that can be taken into account in urban land-use plans is a challenging task that requires collaborative efforts between the seismological and engineering communities. In this study, starting from the results obtained by extensive geophysical and seismological investigations, we propose and apply an approach to the Gubbio basin (Italy) that can be easily implemented for cases of moderate-to-low ground motion and that takes into account not only simple 1D, but also more complicated 3D effects. With this method, the sites inside the basin are classified by their fundamental resonance frequencies, estimated from the horizontal-to-vertical spectral ratio applied to noise recordings (HVNSR). The correspondence between estimates of the fundamental frequency from this method and those derived from earthquake recordings was verified at several calibration sites. The amplification factors used to correct the response spectra are computed by the ratio between the response spectra at sites within the basin and the response spectra at a hard-rock site using data from two seismic transects. Empirical amplification functions are then assigned to the fundamental frequencies after applying an interpolation technique. The suitability of the estimated site-specific correction factors for response spectra was verified by computing synthetic response spectra for stations within the basin, starting from the synthetic recording at a nearby rock station, and comparing them with observed ones.

Description: A temporal seismic network recorded local seismicity along a 130 km long segment of the transpressional dextral strike-slip Liquiñe-Ofqui fault zone (LOFZ) in southern Chile. Seventy five shallow crustal events with magnitudes up to M(tief)w 3.8 and depths shallower than 25 km were observed in an 11-month period mainly occurring in different clusters. Those clusters are spatially related to the LOFZ, to the volcanoes Chaitén, Michinmahuida and Corcovado, and to active faulting on secondary faults. Further activity along the LOFZ is indicated by individual events located in direct vicinity of the surface expression of the LOFZ. Focal mechanisms were calculated using deviatoric moment tensor inversion of body wave amplitude spectra which mostly yield strike-slip mechanisms indicating a NE–SW direction of the P-axis for the LOFZ at this latitude. The seismic activity reveals the present-day activity of the fault zone. The recent M(tief)w 6.2 event near Puerto Aysén, Southern Chile at 45.4°S on April 21, 2007 shows that the LOFZ is also capable of producing large magnitude earthquakes and therefore imposing significant seismic hazard to this region.

Description: We thank Dahlen & Nolet for the comments (DN05) on our paper (HH05). There are many points of agreement, as we think is clear from HH05, but we respectfully continue to differ in opinion on some fundamental aspects of the finite frequency sensitivity kernels known as 'banana doughnut' kernels—hereinafter BDKs, as per the original nomenclature of Dahlen et al.—and their benefit to global tomography. In contrast to DN05's summary statement, HH05's main concern about BDKs is not the effect of uncertainty in the earthquake source signature or origin time. HH05 argue that (i) the evaluation of sensitivity kernels in simple media has limitations for the interpretation of broad-band signals by means of (linearized) finite frequency tomography; (ii) finite frequency kernels are (indeed) oscillatory, but in general heterogeneity their structure will be complex and different from BD features; (iii) the resolved length scales of model variations are induced by the spectral scales present in the data, which makes the notion of 'hole' irrelevant; and (iv) with the need for 'damping' (regularization) and without a basis that matches properly the multi-scale aspects of finite frequency sensitivity, ray theory or finite frequency theory inversions are likely to yield results that are practically the same.

Description: The 26 December 2004 Sumatra-Andaman earthquake of Mw 9.3 triggered a massive tsunami in the Indian Ocean. We here report on observations of the Indian Ocean tsunami at broadband seismic stations located on islands in the area. The tsunami induces long-period (〉1000 s) signals on the horizontal components of the sensor. Frequency-time analysis shows that the long-period signals cannot be due to seismic surface waves, but that it arrives at the expected time of the tsunami. The waveforms are well correlated to tide gauge observations at a location where both observations are available. To explain the signals we favour tilt due to coastal loading but we cannot at the present stage exclude gravitational effects. The density of broadband stations is expected to increase rapidly in the effort of building an earthquake monitoring system. They may unexpectedly become useful tsunami detectors as well.

Description: Magnetotelluric (MT) data from 66 sites along a 45-km-long profile across the San Andreas Fault (SAF) were inverted to obtain the 2-D electrical resistivity structure of the crust near the San Andreas Fault Observatory at Depth (SAFOD). The most intriguing feature of the resistivity model is a steeply dipping upper crustal high-conductivity zone flanking the seismically defined SAF to the NE, that widens into the lower crust and appears to be connected to a broad conductivity anomaly in the upper mantle. Hypothesis tests of the inversion model suggest that upper and lower crustal and upper-mantle anomalies may be interconnected.We speculate that the high conductivities are caused by fluids and may represent a deep-rooted channel for crustal and/or mantle fluid ascent. Based on the chemical analysis of well waters, it was previously suggested that fluids can enter the brittle regime of the SAF system from the lower crust and mantle. At high pressures, these fluids can contribute to fault-weakening at seismogenic depths. These geochemical studies predicted the existence of a deep fluid source and a permeable pathway through the crust. Our resistivity model images a conductive pathway, which penetrates the entire crust, in agreement with the geochemical interpretation. However, the resistivity model also shows that the upper crustal branch of the high-conductivity zone is locatedNEof the seismically defined SAF, suggesting that the SAF does not itself act as a major fluid pathway. This interpretation is supported by both, the location of the upper crustal highconductivity zone and recent studies within the SAFOD main hole, which indicate that pore pressures within the core of the SAF zone are not anomalously high, that mantle-derived fluids are minor constituents to the fault-zone fluid composition and that both the volume of mantle fluids and the fluid pressure increase to the NE of the SAF.We further infer from the MT model that the resistive Salinian block basement to the SW of the SAFOD represents an isolated body, being 5–8 km wide and reaching to depths 〉7 km, in agreement with aeromagnetic data. This body is separated from a massive block of Salinian crust farther to the SW. The NE terminus of resistive Salinian crust has a spatial relationship with a near-vertical zone of increased seismic reflectivity ∼15 km SW of the SAF and likely represents a deep-reaching fault zone.

Description: The passive margin of the South Atlantic shows typical features of a rifted volcanic continental margin, encompassing seaward dipping reflectors, continental flood basalts and high-velocity/density lower crust at the continent–ocean transition, probably emplaced during initial seafloor spreading in the Early Cretaceous. The Springbok profile offshore western South Africa is a combined transect of reflection and refraction seismic data. This paper addresses the analysis of the seismic velocity structure in combination with gravity modelling and isostatic modelling to unravel the crustal structure of the passive continental margin from different perspectives. The velocity modelling revealed a segmentation of the margin into three distinct parts of continental, transitional and oceanic crust. As observed at many volcanic margins, the lower crust is characterised by a zone of high velocities with up to 7.4 km/s. The conjunction with gravity modelling affirms the existence of this body and at the same time substantiated its high densities, found to be 3100 kg/m3. Both approaches identified the body to have a thickness of about 10 km. Yet, the gravity modelling predicted the transition between the high-density body towards less dense material farther west than initially anticipated from velocity modelling and confirmed this density gradient to be a prerequisite to reproduce the observed gravity signal. Finally, isostatic modelling was applied to predict average crustal densities if the margin was isostatically balanced. The results imply isostatic equilibrium over large parts of the profile, smaller deviations are supposed to be compensated regionally. The calculated load distribution along the profile implies that all pressures are hydrostatic beneath a depth of 45 km.

Description: We propose a new rapid procedure for determining the energy magnitude Me for shallow events from broadband teleseismic P-wave signals within the distance range 20°–98°. To accomplish this task, we compute spectral amplitude decay functions for different periods using numerical simulations based on the reference Earth model AK135Q. By means of these functions, we correct the spectra of the teleseismic recordings for the propagation path effects, and calculate the radiated seismic energy ES, and hence Me. We use cumulative P-wave windows for simulating a real- or near real-time procedure and test it for 61 shallow earthquakes. The results show that our approach is able to provide a rapid and reliable Me determination within 7–15 minutes after the earthquake origin time, and is therefore suitable for implementation in rapid response systems.

Description: Detailed images of the lithosphere beneath the western Bohemian Massif were obtained by analysis of more than 8500 P receiver functions. At the intersection of Regensburg-Leipzig-Rostock zone and Eger Rift, crustal thickness decreases to 26 km from approx. 31 km in the surroundings. The receiver functions display a positive phase at about 6 s delay time and a strong negative phase at 7 to 8 s, which coincides with an area of Moho updoming, CO2 mantle-derived degassing and earthquake swarm activity. These phases can be modeled by a velocity increase at 50 km and a velocity decrease at 65 km depth. The velocity decrease, observed over an area of 5300 km2, gives evidence for local asthenospheric updoming and/or a confined body of partial melt, which might be the cause for high CO2 mantle fluid flow and earthquake swarm activity in this recently nonvolcanic, intracontinental rift area.

Description: Compound-specific isotope analysis has become an important tool in environmental studies and is an especially powerful way to evaluate biodegradation of hydrocarbons. Here, carbon isotope ratios of light hydrocarbons were used to characterise in-reservoir biodegradation in the Gullfaks oil field, offshore Norway. Increasing biodegradation, as characterised, for example, by increasing concentration ratios of Pr/n-C17 and Ph/n-C18, and decreasing concentrations of individual light hydrocarbons were correlated to 13C-enrichment of the light hydrocarbons. The δ13C values of C4 to C9n-alkanes increase by 7–3‰ within the six oil samples from the Brent Group of the Gullfaks oil field, slight changes (1–3‰) being observed for several branched alkanes and benzene, whereas no change (〈1‰) in δ13C occurs for cyclohexane, methylcyclohexane, and toluene. Application of the Rayleigh equation demonstrated high to fair correlation of concentration and isotope data of i- and n-pentane, n-hexane, and n-heptane, documenting that biodegradation in reservoirs can be described by the Rayleigh model. Using the appropriate isotope fractionation factor of n-hexane, derived from laboratory experiments, quantification of the loss of this petroleum constituent due to biodegradation is possible. Toluene, which is known to be highly susceptible to biodegradation, is not degraded within the Gullfaks oil field, implying that the local microbial community exhibits rather pronounced substrate specificities. The evaluation of combined molecular and isotopic data expands our understanding of the anaerobic degradation processes within this oil field and provides insight into the degradative capabilities of the microorganisms. Additionally, isotope analysis of unbiodegraded to slightly biodegraded crude oils from several oil fields surrounding Gullfaks illustrates the heterogeneity in isotopic composition of the light hydrocarbons due to source effects. This indicates that both source and also maturity effects have to be well constrained when using compound-specific isotope analysis for the assessment of biodegradation.

Description: In this study we present the new tomographic code ANITA which provides 3-D anisotropic P and isotropic S velocity distribution based on P and S traveltimes from local seismicity. For the P anisotropic model, we determine four parameters for each parameterization cell. This represents an orthorhombic anisotropy with one predefined direction oriented vertically. Three of the parameters describe slowness variations along three horizontal orientations with azimuths of 0°, 60°, and 120°, and one is a perturbation along the vertical axis. The nonlinear iterative inversion procedure is similar to that used in the LOTOS code. We have implemented this algorithm for the updated data set of central Java, part of which was previously used for the isotropic inversion. It was obtained that the crustal and uppermost mantle velocity structure beneath central Java is strongly anisotropic with 7–10% of maximal difference between slow and fast velocity in different directions. In the forearc (area between southern coast and volcanoes), the structure of both isotropic and anisotropic structure is strongly heterogeneous. Variety of anisotropy orientations and highly contrasted velocity patterns can be explained by a complex block structure of the crust. Beneath volcanoes we observe faster velocities in vertical direction, which is probably an indicator for vertically oriented structures (channels, dykes). In the crust beneath the middle part of central Java, north to Merapi and Lawu volcanoes, we observe a large and very intense anomaly with a velocity decrease of up to 30% and 35% for P and S models, respectively. Inside this anomaly E-W orientation of fast velocity takes place, probably caused by regional extension stress regime. In a vertical section we observe faster horizontal velocities inside this anomaly that might be explained by layering of sediments and/or penetration of quasi-horizontal lenses with molten magma. In the mantle, trench parallel anisotropy is observed throughout the study area. Such anisotropy in the slab entrained corner flow may be due to presence of B-type olivine having predominant axis parallel to the shear direction, which appears in conditions of high water or/and melting content.

Description: As part of the South American Geodynamic Activities project we observed the present day deformation field in the territories of Chile and Argentina using the Global Positioning System. The results clearly show that the earthquake cycle dominates the contemporary surface deformation of the central and southern Andes. Compared to geological timescales, the transient elastic deformation related to subduction earthquakes presents a short-term signal which can be explained by interseismic, coseismic, and postseismic phases of interplate thrust earthquakes. We constructed the Andean Elastic Dislocation Model (AEDM) in order to subtract the interseismic loading from the observed velocities. The estimated parameters of the AEDM, and the amount and depth of coupling between the subducting Nazca and overriding South American Plates, represent long-term features and show that the seismogenic interface between both plates is fully locked and that the depth of coupling increases from north to south. The prominent signals in the residual velocity field (i.e. observed velocities minus AEDM) are obviously due to postseismic relaxation processes; they are visible in the area of the 1995 Mw 8.0 Antofagasta earthquake and in the area of the 1960 Mw 9.5 Valdivia earthquake. Although postseismic deformations, compared to geologic timescales, are short-term signals, those signals are valuable constraints on important long-term features of Andean evolution, i.e., the viscosity of the upper mantle and lower crust. The observed surface data are best fitted with a three-dimensional finite element model in which we incorporate a mantle viscosity of 4 × 1019 Pa s. The most obvious long-term deformation signal is manifested in the back-arc of the subduction zone where the Brazilian Shield thrusts beneath the Subandean zone. The style and amount of backarc shortening changes along strike of the orogen, increasing from zero in the south (latitude 〈 −38° S) to values in the order of 10 mm yr−1 close to the Bolivian Orocline. In the fore-arc, whilst we see indications for long-term E-W extension, we did not find any apparent slip partitioning. In addition to this long-term signal, we suggest that the asymmetry of interseismic and coseismic deformation may lead to tectonic structures in the fore-arc. If the coseismic deformation does not release all of the accumulated deformation, then, over many earthquake cycles, part of the interseismic deformation may be transformed into permanent long-term plastic deformation.

Description: The lithospheric structure of the Aegean region is investigated by analysis of Rayleigh-wave fundamental mode dispersion measurements. Isotropic 1-D models for almost 100 two-station ray paths across the region display distinct variations in the Moho depth and crustal S-wave velocities. The descending slab of the subducting African plate can be resolved down to 120 km depth beneath the volcanic arc. Three different regions are distinguished in terms of Moho depth: (1) The forearc, with large crustal thicknesses between 38 and 48 km and an average of 43 km, (2) the northern Aegean, with an average Moho depth of 28 km and (3) the southern Aegean (central volcanic arc, i.e. Cyclades, and Sea of Crete) with an even thinner crust of around 25 km. Lateral variations in structure between 25 and 55 km depth indicate a marked difference between the western and eastern forearc, collocated with pronounced changes in trench and slab geometry as well as published deformation rates. S velocities between 25 and 55 km depth are low everywhere beneath the forearc but increase from the Peleponnesus to Crete. An abrupt change occurs between western and central Crete in terms of the visibility of the Aegean Moho and the seismic structure of the lithospheric mantle wedge: An Aegean mantle wedge with S velocities above 4.4 km s−1 is only observed to the east of central Crete, whereas to the west velocities of less than 4.0 km s−1 occur down to the plate contact. These low velocities above the slab may indicate the presence of a melange of metamorphic rocks at the depths. A low-velocity asthenospheric layer is observed beneath the Sea of Crete and the Cyclades below 40 km depth, between the thinned lithosphere above and the slab below. The observed radial anisotropy in the northern part of the Aegean is likely to be due to preferred orientation of anisotropic minerals within the lower crust, possibly caused by lateral ductile flow associated with recent lithospheric extension.

Description: The North Anatolian Fault Zone (NAFZ) below the Sea of Marmara forms a “seismic gap” where a major earthquake is expected to occur in the near future. This segment of the fault lies between the 1912 Ganos and 1999 İzmit ruptures and is the only NAFZ segment that has not ruptured since 1766. To monitor the microseismic activity at the main fault branch offshore of Istanbul below the Çınarcık Basin, a permanent seismic array (PIRES) was installed on the two outermost Prince Islands, Yassiada and Sivriada, at a few kilometers distance to the fault. In addition, a temporary network of ocean bottom seismometers was deployed throughout the Çınarcık Basin. Slowness vectors are determined combining waveform cross correlation and P wave polarization. We jointly invert azimuth and traveltime observations for hypocenter determination and apply a bootstrap resampling technique to quantify the location precision. We observe seismicity rates of 20 events per month for M 〈 2.5 along the basin. The spatial distribution of hypocenters suggests that the two major fault branches bounding the depocenter below the Çınarcık Basin merge to one single master fault below ∼17 km depth. On the basis of a cross-correlation technique we group closely spaced earthquakes and determine composite focal mechanisms implementing recordings of surrounding permanent land stations. Fault plane solutions have a predominant right-lateral strike-slip mechanism, indicating that normal faulting along this part of the NAFZ plays a minor role. Toward the west we observe increasing components of thrust faulting. This supports the model of NW trending, dextral strike-slip motion along the northern and main branch of the NAFZ below the eastern Sea of Marmara.

Description: Here we present the results of local source tomographic inversion beneath central Java. The data set was collected by a temporary seismic network. More than 100 stations were operated for almost half a year. About 13,000 P and S arrival times from 292 events were used to obtain three-dimensional (3-D) Vp, Vs, and Vp/Vs models of the crust and the mantle wedge beneath central Java. Source location and determination of the 3-D velocity models were performed simultaneously based on a new iterative tomographic algorithm, LOTOS-06. Final event locations clearly image the shape of the subduction zone beneath central Java. The dipping angle of the slab increases gradually from almost horizontal to about 70°. A double seismic zone is observed in the slab between 80 and 150 km depth. The most striking feature of the resulting P and S models is a pronounced low-velocity anomaly in the crust, just north of the volcanic arc (Merapi-Lawu anomaly (MLA)). An algorithm for estimation of the amplitude value, which is presented in the paper, shows that the difference between the fore arc and MLA velocities at a depth of 10 km reaches 30% and 36% in P and S models, respectively. The value of the Vp/Vs ratio inside the MLA is more than 1.9. This shows a probable high content of fluids and partial melts within the crust. In the upper mantle we observe an inclined low-velocity anomaly which links the cluster of seismicity at 100 km depth with MLA. This anomaly might reflect ascending paths of fluids released from the slab. The reliability of all these patterns was tested thoroughly.

Description: We employ P to S converted waveforms to investigate effects of the hot mantle plume on seismic discontinuities of the crust and upper mantle. We observe the Moho at depths between 13 and 17 km, regionally covered by a strong shallow intracrustal converted phase. Coherent phases on the transverse component indicate either dipping interfaces, 3-D heterogeneities or lower crustal anisotropy. We find anomalies related to discontinuities in the upper mantle down to the transition zone evidently related to the hot mantle plume. Lithospheric thinning is confirmed in greater detail than previously reported by Li et al., and we determine the dimensions of the low-velocity zone within the asthenosphere with greater accuracy. Our study mainly focuses on the temperature-pressure dependent discontinuities of the upper mantle transition zone. Effects of the hot diapir on the depths of mineral phase transitions are verified at both major interfaces at 410 and 660 km. We determine a plume radius of about 200 km at the 660 km discontinuity with a core zone of about 120 km radius. The plume conduit is located southwest of Big Island. A conduit tilted in the northeast direction is required in the upper mantle to explain the observations. The determined positions of deflections of the discontinuities support the hypothesis of decoupled upper and lower mantle convection.

Description: Using three different short-period electromagnetic sensors with resonance frequencies of 1 Hz (Mark L4C-3D), 2 Hz (Mark L-22D), and 4.5 Hz (I/O SM-6), coupled with three digital acquisition system, the portable data acquisition system (PDAS) Teledyne Geotech, the refraction technology (REFTEK) 72A, and the Earth Data Logger PR6-24 (EDL), the effect of the seismic instruments on the horizontal-to-vertical spectral ratio (H/V) using seismic noise for frequencies less than 1 Hz has been evaluated. For all possible sensors–acquisition system pairs, the background seismic signal and instrumental self-noise power spectral densities have been calculated and compared. The results obtained when coupling the short-period sensors with different acquisition systems show that the performance of the considered instruments at frequencies 〈1 Hz strongly depends upon the sensor–acquisition system combination and the gain used, with the best performance obtained for sensors with the lowest resonance frequency. For all acquisition systems, it was possible to retrieve correctly the H/V peak down to 0.1–0.2 Hz by using a high gain and a 1-Hz sensor. In contrast, biased H/V spectral ratios were retrieved when low-gain values were considered. Particular care is required when using 4.5-Hz sensors, because they may not even allow the fundamental resonance frequency peak to be reproduced.

Description: Combined P and S receiver functions from seismograms of teleseismic events recorded at 65 temporary and permanent stations in the Aegean region are used to map the geometry of the subducted African and the overriding Aegean plates. We image the Moho of the subducting African plate at depths ranging from 40 km beneath southern Crete and the western Peloponnesus to 160 km beneath the volcanic arc and 220 km beneath northern Greece. However, the dip of the Moho of the subducting African plate is shallower beneath the Peloponnesus than beneath Crete and Rhodes and flattens out beneath the northern Aegean. Observed P-to-S conversions at stations located in the forearc indicate a reversed velocity contrast at the Moho boundary of the Aegean plate, whereas this boundary is observed as a normal velocity contrast by the S-to-P conversions. Our modeling suggests that the presence of a large amount of serpentinite (more than 30%) in the forearc mantle wedge, which generally occurs in the subduction zones, may be the reason for the reverse sign of the P-to-S conversion coefficient. Moho depths for the Aegean plate show that the southern part of the Aegean (crustal thickness of 20–22 km) has been strongly influenced by extension, while the northern Aegean Sea, which at present undergoes the highest crustal deformation, shows a relatively thicker crust (25–28 km). This may imply a recent initiation of the present kinematics in the Aegean. Western Greece (crustal thickness of 32–40 km) is unaffected by the recent extension but underwent crustal thickening during the Hellenides Mountains building event. The depths of the Aegean Moho beneath the margin of the Peloponnesus and Crete (25–28 and 25–33 km, respectively) show that these areas are also likely to be affected by the Aegean extension, even though the Cyclades (crustal thickness of 26–30 km) were not significantly involved in this episode. The Aegean lithosphere-asthenosphere boundary (LAB) mapped with S receiver functions is about 150 km deep beneath mainland Greece, whereas the LAB of the subducted African plate dips from 100 km beneath Crete and the southern Aegean Sea to about 225 km under the volcanic arc. This implies a thickness of 60–65 km for the subducted African lithosphere, suggesting that the Aegean lithosphere was not significantly affected by the extensional process associated with the exhumation of metamorphic core complexes in the Cyclades.

Description: It has recently been shown that correlations of seismic noise can contain significant information about the Green's function along the station profile. Using an array of 38 temporary broad-band stations located in Finland between 1998 September and 1999 March, we study the resulting 703 noise correlations to understand how they are influenced by the directivity of the noise field. The latter information is obtained through f-k analysis of data from two permanent seismic arrays in Germany and Norway and from a subset of stations of the array in Finland. Both types of analysis confirm that the characteristic of the seismic noise is strongly frequency-dependent. At low frequencies (0.02-0.04 Hz), we observe diffuse noise and/or randomly distributed sources. In contrast, the noise is strongly direction-dependent and not fully diffuse in the intermediate period ranges (0.04-0.25 Hz) which correspond to the first and second micro seismic peak, created at the Irish and Scottish coast and the western coast of Norway In this frequency interval the noise is sufficiently close to a plane wave to introduce systematic errors in group velocities for station pairs which are not parallel to the direction of the dominant incident noise. Phase velocities calculated by slant stack over many traces are however correct, independently of profile direction. In the high-frequency band (0.25-1.0 Hz), the situation is a mix between the low-frequency and the intermediate frequency cases. Average phase velocities and individual group velocities from well-oriented profiles are in excellent agreement with results from Rayleigh wave studies of the area.

Description: This study focuses on the present-day deformation mechanisms of the south central Chile margin, at the transition zone between two megathrust earthquake segments defined from historical data: the Valdivia and Concepción sectors. New GPS data and finite-element models with complex geometries constrained by geophysical data are presented to gain insight into forearc kinematics and to address the role of upper plate faults on contemporary deformation. GPS vectors are heterogeneously distributed in two domains that follow these two earthquake segments. We find that models which simulate only interseismic locking on the plate interface fail to reproduce surface deformation in the entire study area. In the Concepción domain, models that include a crustal-scale fault in the upper plate better reproduce the GPS observations. In the Valdivia domain, GPS data show regional-scale vertical axis rotations, which could reflect postseismic deformation processes at the edge of the Mw 9.5 earthquake that ruptured in 1960 and/or activity of another crustal fault related to motion of a forearc sliver. Our study suggests that upper plate faults in addition to earthquake cycle transients may exert an important control on the surface velocity of subduction zone forearcs.

Description: The volcanic arc of the Hellenic subduction zone with its four volcanic centers is of major relevance when evaluating the seismovolcanic hazard for the Aegean region. We present results from a 22-station temporary seismic network (CYCNET) in the central Hellenic Volcanic Arc (HVA). CYCNET recordings allow to analyze the level and spatio-temporal evolution of microseismic activity in this region for the first time. A total of 2175 events recorded between September 2002 and July 2004 are analyzed using statistical methods, cluster analysis and relative relocation techniques. We identify distinct regions with significantly varying spatio-temporal behavior of microseismicity. A large portion of the seismic activity within the upper crust is associated with the presence of islands representing horst structures that were generated during the major Oligocene extensional phase. In contrast, the central part of the Cyclades metamorphic core complex remains aseismic considering our magnitude threshold of 1.8 except one spot where events occur swarm-like and with highly similar waveforms. The highest activity in the study area was identified along the SW–NE striking Santorini–Amorgos zone. Within this zone the submarine Columbo volcano exhibits strong temporal variations of seismic activity on a high background level. This activity is interpreted to be directly linked to the magma reservoir and therein the migration of magma and fluids towards the surface. NE of Columbo where no volcanic activity has yet been reported we observe a similar seismicity pattern with small-scaled activity spots that might represent local pathways of upward migrating fluids or even developing volcanic activity within this zone of crustal weakness. In contrast, the Santorini and Milos volcanic complexes do not show significant temporal variations and low to moderate background activity, respectively. Relating our results to the distribution of historical earthquakes and the GPS-derived horizontal velocity field we conclude that the Santorini–Amorgos zone is presently in the state of right-lateral transtension reflecting a major structural boundary of the volcanic arc subdividing it into a seismically and volcanically quiet western and an active eastern part.

Description: Seismic anisotropy was investigated by measuring shear-wave splitting at 19 broadband stations in Greenland. We examined mostly SKS and SKKS phases, but also some PKS and depth phases of SKS (e.g. pSKS, sSKS) for deep events. Splitting parameters (fast polarization and time delay) were determined for these phases. The fast polarizations at nine sites in southern Greenland are quite uniformly oriented about N–NE. Two sites in central northern Greenland show a similar geometry to southern Greenland. Similar fast polarizations in southern and central northern Greenland suggest continuity of structural fabric beneath large parts of Greenland. This coherent pattern extends across a number of geological provinces of varying age and suggests a common cause of anisotropy not related to the bitwise formation of the Greenland continental block. Four sites in an east–west oriented belt crossing central Greenland show varying fast polarizations and suggest a separate process causing the anisotropy there, which may indicate that these processes are not currently active. The overall pattern of anisotropy in our results, with the exception of variations across central Greenland, is similar to results obtained from Rayleigh waves. The irregular geometry of splitting across central Greenland may be related to the impact of the Iceland plume at ∼ 60 Ma. Reported splitting time delays range from 0.4 to 1.4 s with an average of 0.8 s, which can generally not be explained by crustal anisotropy alone. If confined to a lithosphere of thickness on the order of 100 km, time delays of up to 1.4 s indicate anisotropy of up to about 6%, assuming that the a crystallographic axis of olivine is preferentially contained in the horizontal plane. We suggest that the anisotropy beneath Greenland is located mainly in the upper mantle but some contributions from the crust and lower mantle may be present.

Description: We report on a receiver function study of the crust and upper mantle within DESERT, a multidisciplinary geophysical project to study the lithosphere across the Dead Sea Transform (DST). A temporary seismic network was operated on both sides of the DST between 2000 April and 2001 June. The depth of the Moho increases smoothly from about 30 to 34-38 km towards the east across the DST, with significant north-south variations east of the DST. These Moho depth estimates from receiver functions are consistent with results from steep- and wide-angle controlled-source techniques. Steep-angle reflections and receiver functions reveal an additional discontinuity in the lower crust, but only east of the DST. This leads to the conclusion that the internal crustal structure east and west of the DST is different. The P to S converted phases from both discontinuities at 410 and 660 km are delayed by 2 s with respect to the IASP91 global reference model. This would indicate that the transition zone is consistent with the global average, but the upper mantle above 410 km is 3-4 per cent slower than the standard earth model.

Description: TOPO-EUROPE addresses the 4-D topographic evolution of the orogens and intra-plate regions of Europe through a multidisciplinary approach linking geology, geophysics, geodesy and geotechnology. TOPO-EUROPE integrates monitoring, imaging, reconstruction and modelling of the interplay between processes controlling continental topography and related natural hazards. Until now, research on neotectonics and related topography development of orogens and intra-plate regions has received little attention. TOPO-EUROPE initiates a number of novel studies on the quantification of rates of vertical motions, related tectonically controlled river evolution and land subsidence in carefully selected natural laboratories in Europe. From orogen through platform to continental margin, these natural laboratories include the Alps/Carpathians–Pannonian Basin System, the West and Central European Platform, the Apennines–Aegean–Anatolian region, the Iberian Peninsula, the Scandinavian Continental Margin, the East-European Platform, and the Caucasus–Levant area. TOPO-EUROPE integrates European research facilities and know-how essential to advance the understanding of the role of topography in Environmental Earth System Dynamics. The principal objective of the network is twofold. Namely, to integrate national research programs into a common European network and, furthermore, to integrate activities among TOPO-EUROPE institutes and participants. Key objectives are to provide an interdisciplinary forum to share knowledge and information in the field of the neotectonic and topographic evolution of Europe, to promote and encourage multidisciplinary research on a truly European scale, to increase mobility of scientists and to train young scientists. This paper provides an overview of the state-of-the-art of continental topography research, and of the challenges to TOPO-EUROPE researchers in the targeted natural laboratories.

Description: We examine a 24-hour period of active San Andreas Fault (SAF) tremor and show that this tremor is largely composed of repeated similar events. Utilizing this similarity, we locate the subset of the tremor with waveforms similar to an identified low frequency earthquake (LFE) “master template,” located using P and S wave arrivals to be ∼26 km deep. To compensate for low signal-to-noise, we estimate event-pair differential times at “clusters” of nearby stations rather than at single stations. We find that the locations form a near-linear structure in map view, striking parallel to the SAF and near the surface trace. Therefore, we suggest that at least a portion of the tremor occurs on the deep extension of the fault, likely reflecting shear slip, similar to subduction zone tremor. If so, the SAF may extend to the base of the crust, ∼10 km below the deepest regular earthquakes on the fault.

Description: The Seychelles plateau is a prime example of a microcontinent, yet mechanisms for its creation and evolution are poorly understood. Recently acquired teleseismic data from a deployment of 26 stations on 18 islands in the Seychelles are analyzed to study upper mantle seismic anisotropy using SKS splitting results. Strong microseismic noise is attenuated using a polarization filter. Results show significant variation in time delays (dt =0.4–2.4 s) and smooth variations in orientation (f = 1569 where f is the polarization of the fast shear wave). The splitting results cannot be explained by simple asthenospheric flow associated with absolute plate motions. Recent work has suggested that anisotropy measurements for oceanic islands surrounding Africa can be explained by mantle flow due to plate motion in combination with density-driven flow associated with the African superswell. Such a mechanism explains our results only if there are lateral variations in the viscosity of the mantle. It has been suggested that the Seychelles are underlain by a mantle plume. Predictions of flow-induced anisotropy from plumelithosphere interaction can explain our results with a plume possibly impinging beneath the plateau. Finally, we consider lithospheric anisotropy associated with rifting processes that formed the Seychelles. The large variation in the magnitude of shear wave splitting over short distances suggests a shallow source of anisotropy. Fast directions align parallel to an area of transform faulting in the Amirantes. Farther from this area the orientation of anisotropy aligns in similar directions as plate motions. This supports suggestions of transpressive deformation during the opening of the Mascarene basin.

Description: A detailed knowledge of the thickness of the lithosphere in the north Atlantic is an important parameter for understanding plate tectonics in that region. We achieve this goal with as yet unprecedented detail using the seismic technique of S-receiver functions. Clear positive signals from the crust-mantle boundary and negative signals from a mantle discontinuity beneath Greenland, Iceland and Jan Mayen are observed. According to seismological practice, we call the negative phase the lithosphere-ashenosphere boundary (LAB). The seismic lithosphere under most of the Iceland and large parts of central Greenland is about 80 km thick. This depth in Iceland is in disagreement with estimates of the thickness of the elastic lithosphere (10-20 km) found from postglacial rebound data. In the region of flood basalts in eastern Greenland, which overlies the proposed Iceland plume track, the lithosphere is only 70 km thick, about 10 km less than in Iceland which is located directly above the proposed plume. At the western Greenland coast, the lithosphere.

Description: We describe results of an active-source seismology experiment across the Chilean subduction zone at 38.2◦S. The seismic sections clearly show the subducted Nazca plate with varying reflectivity. Below the coast the plate interface occurs at 25 km depth as the sharp lower boundary of a 2–5 km thick, highly reflective region, which we interpret as the subduction channel, that is, a zone of subducted material with a velocity gradient with respect to the upper and lower plate. Further downdip along the seismogenic coupling zone the reflectivity decreases in the area of the presumed 1960 Valdivia hypocentre. The plate interface itself can be traced further down to depths of 50–60 km below the Central Valley. We observe strong reflectivity at the plate interface as well as in the continental mantle wedge. The sections also show a segmented forearc crust in the overriding South American plate. Major features in the accretionary wedge, such as the Lanalhue fault zone, can be identified. At the eastern end of the profile a bright west-dipping reflector lies perpendicular to the plate interface and may be linked to the volcanic arc.

Description: Following the ongoing debate about suitability of short-period sensor for seismic noise measurements at frequencies lower than 1 Hz, in this study we compare recordings from two different seismometers (Güralp CMG-3ESPC and Mark L4C-3D) installed side by side in the GeoForschungsZentrum laboratory. The comparison carried out in terms of Power Spectral Density and coherency analysis shows an excellent agreement between the short-period and the broad-band recordings in the frequency band 0.2–20 Hz. Therefore, this result highlights that with a calibrated short-period sensor one can obtain the same results that would be obtained by using a broad-band seismometer in the band of engineering interest.

Description: The interdisciplinary Dead Sea Rift Transect (DESERT) project that was conducted in Israel, the Palestine Territories and Jordan has provided a rich palette of data sets to examine the crust and uppermost mantle beneath one of Earth’s most prominent fault systems, the Dead Sea Transform (DST). As part of the passive seismic component, thirty broad-band sensors were deployed in 2000 across the DST for roughly one year. During this deployment, we recorded 115 teleseismic earthquakes that are suitable for a fundamental mode Rayleigh wave analysis at intermediate periods (35–150 s). Our initial analysis reveals overall shear velocities that are reduced by up to 4 per cent with respect to reference Earth model PREM. To the west of the DST, we find a seismically relatively fast but thin lid that is about 80 km thick. Towards the east, shallow seismic velocities are low while a deeper low velocity zone is not detected. This contradicts the currently favoured thermomechanical model for the DST that predicts lithospheric thinning through mechanical erosion by an intruding plume from the Red Sea. On the other hand, our current results are somewhat inconclusive regarding asthenosphere velocities east of the DST due to the band limitation of the recording equipment in Jordan.

Description: A set of seismological stations was deployed in the Central Andes region along a ~600 km long profile at 21°S between Chile and Bolivia and operated for a period of almost two years, from March 2002 to January 2004. Here we present the results of the tomographic inversion for P-wave velocity anomalies, based on teleseismic data recorded at the stations. The reliability of the results has been checked by a series of synthetic tests. The tomographic images show high-velocities on the west of the profile that are indicative of cold material from the fore-arc. A low-velocity anomaly is detected at the border between the fore- and the volcanic are where the Quebrada Blanca seismic anomaly was previously described. This anomaly might be related to the presence of fluids that originate at the cluster of earthquakes at a depth of ~100 km in the subducted plate. A strong low-velocity anomaly is detected beneath the entire Altiplano plateau and part of the Eastern Cordillera, in agreement with previous receiver function results. The Brazilian Shield is thought to be responsible for the strong high-velocity anomaly underneath the Interandean and Subandean regions.

Description: On 8 January 2006, an intermediate-depth earthquake occurred at the western part of the Hellenic trench close to the island of Kythera (southern Greece). This is the first intermediate-depth earthquake in the broader Aegean area that has produced such an extensive set of useful recordings, as it was recorded by the main permanent seismological networks and numerous acceleration sensors operating in Greece, as well as by EGELADOS, a large-scale temporary amphibian broadband seismological network deployed in the southern Aegean area. An effort to combine all the available data (broadband velocity and acceleration sensor) was made to study the properties of ground-motion attenuation of this earthquake. The combination of both types of data revealed interesting properties of the earthquake wave field, which would remain hidden if only one type of data was used. Moreover, the data have been used for a validation of existing peak ground-motion empirical prediction relations and the preliminary study of the very inhomogeneous attenuation pattern of the southern Aegean intermediate-depth events at both near- and far-source distances

Description: We used data from both permanent and temporary seismic networks on Iceland and Greenland to investigate the crustal thickness by partly reinterpreting earlier data (P receiver functions) and adding S receiver functions to better constrain the results. We obtained good conversions from the Moho and also crustal multiples in both Iceland and Greenland. The central ice covered part of Greenland has an average crustal thickness of 40 km, typical for a craton. At the edges of Greenland the crustal thickness decreases to 30–40 km. On the east coast of Greenland, where the track of the Iceland plume is thought to have affected the lithosphere, the crustal thickness is only 24–26 km. In contrast to previous studies, we find that the crustal thickness in the east and the northwest coastal regions of Iceland is more than 40 km, similar to beneath the active volcanic region. In the southwest region of Iceland and along the mid-ocean ridge, the crustal thickness is only 25 km or less. Also in contrast to earlier receiver function interpretations, which deduced a broad crust-mantle transition zone for Iceland, we find indications for a normal, sharper Moho beneath a number of sites

Description: Data from 90 permanent broad-band stations spread over central and eastern Europe were analysed using Ps receiver functions to study the crustal and upper-mantle structure down to the mantle transition zone. Receiver functions provide valuable information on structural features, which are important for the resolution of European lithospheric dynamics. Moho depths vary from less than 25 km in extensional areas in central Europe to more than 50 km at stations in eastern Europe (Craton) and beneath the Alpine–Carpathian belt. A very shallow Moho depth can be observed at stations in the Upper Rhine Graben area (ca. 25 km), whereas, for example, stations in the SW Bohemian Massif show a significantly deeper Moho interface at a depth of 38 km. Vp/Vs ratios vary between 1.60 and 1.96, and show no clear correlation to the major tectonic units, thus probably representing local variations in crustal composition. Delayed arrivals of converted phases from the mantle transition zone are observed at many stations in central Europe, whereas stations in the cratonic area show earlier arrivals compared with those calculated from the IASP91 Earth reference model. Differential delay times between the P410s and P660s phases indicate a thickened mantle transition zone beneath the eastern Alps, the Carpathians and the northern Balkan peninsula, whereas the transition zone thickness in eastern and central Europe agrees with the IASP91 value. The thickening of the mantle transition zone beneath the eastern Alps and the Carpathians could be caused by cold, deeply subducted oceanic slabs.

Description: The Andes were formed by Cenozoic tectonic shortening of the South American plate margin overriding the subducting Nazca Plate. Using coupled, thermo-mechanical, numerical modeling of the dynamic interaction between subducting and overriding plates, we searched for factors controlling the intensity of the tectonic shortening. From our modeling, constrained by geological and geophysical observations, we infer that the most important factor was fast and accelerating (from 2 to 3 cm yr(hoch)-1) westward drift of the South American Plate, whereas possible changes in the convergence rate were not as important. Other important factors are the crustal structure of the overriding plate and the shear coupling at the plate interface. The model in which the South American Plate has a thick (40-45 km at 35 Ma) crust and relatively high friction coefficient (0.05) at the Nazca-South American plate interface generates more than 300 km of tectonic shortening over the past 35 million years and replicates well the crustal structure and evolution of the high Central Andes, However, modeling does not confirm that possible climate-controlled changes to the sedimentary trench-fill during the last 30 million years might have significantly influenced the upperplate shortening rate. The model with initially thinner (less than 40 km) continental crust and a lower friction coefficient (less than 0.015) results in less than 40 ikm of shortening in the South American Plate, replicating the situation in the Southern Andes. During upper-plate deformation, the processes that cause a reduction in lithospheric strength and an increase in interplate coupling are particularly important. The most significant of these processes appears to be: (1) delamination of the lower crust and mantle lithosphere, driven by gabbro-eclogite transformation in the thickening lower crust, and (2) mechanical failure of the foreland sediments. The modeling demonstrates that delaminating lithosphere interacts with subduction-zone corner flow, influencing both the rate of tectonic shortening and magmatic-arc productivity, and suggests an anti-correlation between these two parameters. Our model also predicts that the down-dip limit of the frictional coupling domain between the Nazca and South American Plates should be ~15-20 km deeper in the Southern Andes (south of 28° S) compared to the high Central Andes, which is consistent with GPS and seismological observations.

Description: Receiver functions (RF) are used to investigate the upper mantle structure beneath the Eifel, the youngest volcanic area of Central Europe. Data from 96 teleseismic events recorded by 242 seismological stations from permanent and a temporary network has been analysed. The temporary network operated from 1997 November to 1998 June and covered an area of approximately 400 × 250 km2 centred on the Eifel volcanic fields. The average Moho depth in the Eifel is approximately 30 km, thinning to ca. 28 km under the Eifel volcanic fields. RF images suggest the existence of a low velocity zone at about 60–90 km depth under the West Eifel. This observation is supported by P- and S-wave tomographic results and absorption (but the array aperture limits the resolution of the tomographic methods to the upper 400 km). There are also indications for a zone of elevated velocities at around 200 km depth, again in agreement with S-wave and absorption tomographic results. This anomaly is not visible in P-wave tomography and could be due to S-wave anisotropy. The RF anomalies at the Moho, at 60–90 km, and near 200 km depth have a lateral extent of about 100 km. The 410 km discontinuity under the Eifel is depressed by 15–25 km, which could be explained by a maximum temperature increase of +200°C to +300°C. In the 3-D RF image of the Eifel Plume we also notice two additional currently unexplained conversions between 410 and 550 km depth. They could represent remnants of previous subduction or anomalies due to delayed phase changes. The lateral extent of these conversions and the depression of the 410 km discontinuity is about 200 km. The 660 km discontinuity does not show any depth deviation from its expected value. Our observations are consistent with interpretation in terms of an upper mantle plume but they do not rule out connections to processes at larger depth.

Description: Teleseismic data recorded during one and a half years are investigated with the receiver function technique to determine the crustal and upper-mantle structures underneath the highly elevated Altiplano and Puna plateaus in the central Andes. A series of converting interfaces are determined along two profiles at 21°S and 25.5°S, respectively, with a station spacing of approximately 10 km. The data provide the highest resolution gained from a passive project in this area, so far. The oceanic Nazca plate is detected down to 120 km beneath the Altiplano whereas beneath the Puna, the slab can unexpectedly be traced down to 200 km depth at longer periods. A shallow crustal low-velocity zone is determined beneath both plateaus exhibiting segmentation. In the case of the Altiplano, the segments present vertical offsets and are separated by inclined interfaces, which coincide with major fault systems at the surface. An average depth to Moho of about 70 km is determined for the Altiplano plateau. A strong negative velocity anomaly located directly below the Moho along with local crustal thinning is interpreted beneath the volcanic arc of the Altiplano plateau between 67°W and 68.5°W. A deep section of the Puna profile reveals thinning of the mantle transition zone. Although poorly resolved, the detected anomaly may suggest the presence of a mantle plume, which may constitute the origin of the anomalous temperatures at the depth of the upper-mantle discontinuities.

Description: To study the applicability of the passive seismic interferometry technique to near-surface geological studies, seismic noise recordings from a small scale 2-D array of seismic stations were performed in the test site of Nauen (Germany). Rayleigh wave Green's functions were estimated for different frequencies. A tomographic inversion of the traveltimes estimated for each frequency from the Green's functions is then performed, allowing the laterally varying 3-D surface wave velocity structure below the array to be retrieved at engineering–geotechnical scales. Furthermore, a 2-D S-wave velocity cross-section is obtained by combining 1-D velocity structures derived from the inversion of the dispersion curves extracted at several points along a profile where other geophysical analyses were performed. It is shown that the cross-section from passive seismic interferometry provides a clear image of the local structural heterogeneities that are in excellent agreement with georadar and geoelectrical results. Such findings indicate that the interferometry analysis of seismic noise is potentially of great interest for deriving the shallow 3-D velocity structure in urban areas.

Description: As part of the DEad Sea Integrated REsearch project (DESIRE) a 235 km long seismic wide-angle reflection/refraction (WRR) profile was completed in spring 2006 across the Dead Sea Transform (DST) in the region of the southern Dead Sea basin (DSB). The DST with a total of about 107 km multi-stage left-lateral shear since about 18 Ma ago, accommodates the movement between the Arabian and African plates. It connects the spreading centre in the Red Sea with the Taurus collision zone in Turkey over a length of about 1 100 km. With a sedimentary infill of about 10 km in places, the southern DSB is the largest pull-apart basin along the DST and one of the largest pull-apart basins on Earth. The WRR measurements comprised 11 shots recorded by 200 three-component and 400 one-component instruments spaced 300 m to 1.2 km apart along the whole length of the E–W trending profile. Models of the P-wave velocity structure derived from the WRR data show that the sedimentary infill associated with the formation of the southern DSB is about 8.5 km thick beneath the profile. With around an additional 2 km of older sediments, the depth to the seismic basement beneath the southern DSB is about 11 km below sea level beneath the profile. Seismic refraction data from an earlier experiment suggest that the seismic basement continues to deepen to a maximum depth of about 14 km, about 10 km south of the DESIRE profile. In contrast, the interfaces below about 20 km depth, including the top of the lower crust and the Moho, probably show less than 3 km variation in depth beneath the profile as it crosses the southern DSB. Thus the Dead Sea pull-apart basin may be essentially an upper crustal feature with upper crustal extension associated with the left-lateral motion along the DST. The boundary between the upper and lower crust at about 20 km depth might act as a decoupling zone. Below this boundary the two plates move past each other in what is essentially a shearing motion. Thermo-mechanical modelling of the DSB supports such a scenario. As the DESIRE seismic profile crosses the DST about 100 km north of where the DESERT seismic profile crosses the DST, it has been possible to construct a crustal cross-section of the region before the 107 km left-lateral shear on the DST occurred.

Description: Mt. Merapi is one of the most dangerous volcanoes in Indonesia, located within the tectonically active region of south-central Java. This study investigates how Mt. Merapi affected - and was affected by - nearby tectonic earthquakes. In 2001, a Mw6.3 earthquake occurred in conjunction with an increase in fumarole temperature at Mt. Merapi. In 2006, another Mw6.3 earthquake took place, concomitant with an increase of magma extrusion and pyroclastic flows. Here, we develop theoretical models to study the amount of stress transfer between the earthquakes and the volcano, showing that dynamic, rather than static, stress changes are likely responsible for the temporal and spatial proximity of these events. Our examination of the 2001 and 2006 events implies that volcanic activity at Mt. Merapi is influenced by stress changes related to remote tectonic earthquakes, a finding that is important for volcano hazard assessment in this densely inhabited area.

Description: From November 1999 through July 2000, a broadband seismic experiment was carried out at Popocatépetl Volcano to record seismic activity over a wide period range (0.04–100 s). We present an overview of the seismicity recorded during this experiment and discuss results of analyses of long-period (LP) and very-long-period (VLP) seismic signals recorded at stations nearest to the crater over a four-month interval December 1999–March 2000. Three families of LP signals (Types-I, II, and III) are identified based on distinctive waveform features observed periods shorter than 1 s, periods longer than 15 s, and within the period range 0.5–2.5 s. Type-I LP events have impulsive first arrivals and exhibit a characteristic harmonic wave train with dominant periods in the 1.4–1.9 s range during the first 10 s of signal. These events are also associated with a remarkable VLP wavelet with period near 30 s. Type-II LP events represent pairs of events occurring in rapid succession and whose signatures are superimposed. These are typically marked by slowly emergent first arrivals and by a characteristic VLP wave train with dominant period near 30 s, made of two successive wavelets whose shapes are quasi-identical to those of the VLP wavelets associated with Type-I events. Type-III LP events represent the most energetic signals observed during our experiment. These have an emergent first arrival and display a harmonic signature with dominant period near 1.1 s. They are dominated by periods in the 0.25–0.35 s band and contain no significant energy at periods longer than 15 s. Hypocentral locations of the three types of LP events obtained from phase picks point to shallow seismic sources clustered at depths shallower than 2 km below the crater floor. Observed variations in volcanic eruptive activity correlate with defined LP families. Most of the observed seismicity consists of Type-I events that occur in association with 1–3-min-long degassing bursts (“exhalations”). Eruptive activity increased in intensity in February, coinciding with an increasing occurrence of Type-II LP events. Type-III events were first observed at the end of February and during March, in coincidence with the formation of a new lava dome. Vulcanian eruptions occurred in April and May. These events typically exhibit broadband signatures extending over the full period range of the sensors and lasting 30–80 min.

Description: The Bohemian Massif is part of the Variscan belt of central Europe. We carried out a high resolution mapping of lithospheric thickness beneath central Europe by investigating 264 teleseismic events recorded at 80 broad band stations in the western Bohemian Massif with the method of S receiver function analysis. A negative phase beneath the Saxothuringian and north-eastern Teplá-Barrandian units at about 9-10 s before the S onset is interpreted as caused by the lithosphere-asthenosphere boundary (LAB) at 80 - 90 km depth. In the Moldanubian unit, the negative phase occurs at 13-15 s before the S onset, corresponding to lithospheric thickness of 120 - 130 km. The boundary between the domains is oriented E-W and probably marks the northern extension of Moldanubian Saxothuringian to the Moldanubian unit. The observed crustal/lithospheric domains could represent two distinct microplates with a relatively sharp boundary cutting through the whole lithosphere.

Description: On 12 May, a great earthquake (Ms=8.0) on the Longmenshan thrust fault rumbled through Chinas Sichuan province, killing more than 69,000 people and injuring 374,000. The Longmenshan thrust is part of the eastern border of the Tibetan Plateau, but it is not the plateaus only restless margin. An even larger earthquake (Ms=8.1) on the Kunlun fault shook northeastern Tibet in 2001, fortunately in a sparsely populated area. These massive quakes underscore the importance of understanding the tectonic response of Asia to collision by India. The International Deep Profiling of Tibet and the Himalaya (INDEPTH) program explores the dynamics of the India-Asia collision. Though many past geophysical studies have focused on the Himalayas and the southern Tibetan Plateau, the INDEPTH IV project examines the deep structure of the northeastern margin of the Tibetan Plateau.

Description: The lithosphere-asthenosphere boundary corresponds to the base of the “rigid” plates – the depth at which heat transport changes from advection in the convecting deeper upper mantle to conduction in the shallow upper mantle. Although this boundary is a fundamental feature of the Earth, mapping it has been difficult because it does not correspond to a sharp change in temperature or composition. Various definitions of the lithosphere and asthenosphere are based on the analysis of different types of geophysical and geological observations. The depth to the lithosphere-asthenosphere boundary determined from these different observations often shows little agreement when they are applied to the same region because the geophysical and geological observations (i.e., seismic velocity, strain rate, electrical resistivity, chemical depletion, etc.) are proxies for the change in rheological properties rather than a direct measure of the rheological properties. In this paper, we focus on the seismic mapping of the upper mantle high velocity lid and low velocity zone and its relationship to the lithosphere and asthenosphere. We have two goals: (a) to examine the differences in how teleseismic body-wave travel-time tomography and surface-wave tomography image upper mantle seismic structure; and (b) to summarise how upper mantle seismic velocity structure can be related to the structure of the lithosphere and asthenosphere. Surface-wave tomography provides reasonably good depth resolution, especially when higher modes are included in the analysis, but lateral resolution is limited by the horizontal wavelength of the long-period surface waves used to constrain upper mantle velocity structure. Teleseismic body-wave tomography has poor depth resolution in the upper mantle, particularly when no strong lateral contrasts are present. If station terms are used, features with large lateral extent and gradual boundaries are attenuated in the tomographic image. Body-wave models are not useful in mapping the thickness of the high velocity upper mantle lid because this type of analysis often determines wave speed perturbations from an unknown horizontal average and not absolute velocities. Thus, any feature which extends laterally across the whole region beneath a seismic network becomes invisible in the Teleseismic body-wave tomographic image. We compare surface-wave and body-wave tomographic results using southern Africa as an example. Surface-wave tomographic images for southern Africa show a strong, high velocity upper mantle lid confined to depths shallower than ~200 km, whereas body-wave tomographic images show weak high velocity in the upper mantle extending to depths of ~300 km or more. However, synthetic tests show that these results are not contradictory. The absolute seismic velocity structure of the upper mantle provided by surface wave analysis can be used to map the thermal lithosphere. Priestley and McKenzie (Priestley, K., McKenzie, D., 2006. The thermal structure of the lithosphere from shear wave velocities. Earth and Planetary Science Letters 244, 285–301.) derive an empirical relationship between shear wave velocity and temperature. This relationship is used to obtain temperature profiles from the surfacewave tomographic models of the continental mantle. The base of the lithosphere is shown by a change in the gradient of the temperature profiles indicative of the depth where the mode of heat transport changes from conduction to advection. Comparisons of the geotherms determined from the conversion of surface-wave wave speeds to temperatures with upper mantle nodule-derived geotherms demonstrate that estimates of lithospheric thickness from Vs and from the nodule mineralogy agree to within about 25 km. The Lithospheric thickness map for Africa derived from the surface-wave tomographic results shows that thick lithosphere underlies most of the Archean crust in Africa. The distribution of diamondiferous kimberlites provides an independent estimate of where thick lithosphere exists. Diamondiferous kimberlites generally occur where the lower part of the thermal lithosphere as indicated by seismology is in the diamond stability field.

Description: We use the S receiver function method to study the lithosphere at the Dead Sea Transform (DST). A temporary network of 22 seismic broad-band stations was operated on both sides of the DST from 2000 to 2001 as part of the DESERT project. We also used data from six additional permanent broad-band seismic stations at the DST and in the surrounding area, that is, in Turkey, Saudi Arabia, Egypt and Cyprus. Clear S-to-P converted phases from the crust-mantle boundary (Moho) and a deeper discontinuity, which we interpret as lithosphere- asthenosphere boundary (LAB) have been observed. The Moho depth (30-38 km) obtained from S receiver functions agrees well with the results from P receiver functions and other geophysical data. We observe thinning of the lithosphere on the eastern side of the DST from 80 km in the north of the Dead Sea to about 65 km at the Gulf of Aqaba. On the western side of the DST, the few data indicate a thin LAB of about 65 km. For comparison, we found a 90-km-thick lithosphere in eastern Turkey and a 160-km-thick lithosphere under the Arabian shield, respectively. These observations support previous suggestions, based on xenolith data, heat flow observations, regional uplift history and geodynamic modelling, that the lithosphere around DST has been significantly thinned in the Late Cenozoic, likely following rifting and spreading of the Red Sea.