ALBERT

Description: In this thesis, the systematic, situation-oriented selection of approaches to sustainability assessment and effects of selection on assessment results are investigated. The central focus lies on the practice-oriented design of a framework to support selection decisions as well as the necessary criteria and scales for the systematic, quantifiable description of assessment approaches and assessment situations within such a framework. Sustainability assessments are important instruments for the derivation of goals, strategies and measures for shaping sustainable development in all domains. They provide decision-makers in science, industry, politics and society with vital answers to sustainability-related questions that arise in the most diverse contexts. Numerous different assessment approaches are available for carrying out sustainability assessments within such assessment situations. Because of the multitude and diversity of assessment situations and approaches, not every approach is fitting for every situation. In current practice, the fit between approaches and assessment situations is not, or only insufficiently, taken into account when selecting sustainability assessment approaches. Furthermore, no systematic studies have yet been conducted on the effects of approach selection on assessment results. The central result of this work is a concept for the situation-oriented selection support of sustainability assessment approaches based on a Multi-Criteria Decision-Making framework. With the framework, "fitness scores" are calculated, which are used to quantify and operationalize the fit between assessment approaches and assessment situations. With the developed concept, different assessment approaches are selected and exemplarily applied within a use case. Hereby, the effects of approach selection on assessment results are examined. On the basis of the knowledge gained with regard to approach selection, framework development and application, the potentials and limits of assessment approaches widely fitting for diverse assessment situations are finally derived.

Description: On 12 September 2007, an Mw 8.4 earthquake occurred within the southern section of the Mentawai segment of the Sumatra subduction zone, where the subduction thrust had previously ruptured in 1833 and 1797. Traveltime data obtained from a temporary local seismic network, deployed between December 2007 and October 2008 to record the aftershocks of the 2007 event, was used to determine two-dimensional (2-D) and three-dimensional (3-D) velocity models of the Mentawai segment. The seismicity distribution reveals significant activity along the subduction interface and within two clusters in the overriding plate either side of the forearc basin. The downgoing slab is clearly distinguished by a dipping region of high Vp (8.0 km/s), which can be a traced to ∼50 km depth, with an increased Vp/Vs ratio (1.75 to 1.90) beneath the islands and the western side of the forearc basin, suggesting hydrated oceanic crust. Above the slab, a shallow continental Moho of less than 30 km depth can be inferred, suggesting that the intersection of the continental mantle with the subducting slab is much shallower than the downdip limit of the seismogenic zone despite localized serpentinization being present at the toe of the mantle wedge. The outer arc islands are characterized by low Vp (4.5–5.8 km/s) and high Vp/Vs (greater than 2.0), suggesting that they consist of fluid saturated sediments. The very low rigidity of the outer forearc contributed to the slow rupture of the Mw 7.7 Mentawai tsunami earthquake on 25 October 2010.

Description: On 2012 May 20 and 29, two damaging earthquakes with magnitudes Mw 6.1 and 5.9, respectively, struck the Emilia-Romagna region in the sedimentary Po Plain, Northern Italy, causing 26 fatalities, significant damage to historical buildings and substantial impact to the economy of the region. The earthquake sequence included four more aftershocks with Mw ≥ 5.0, all at shallow depths (about 7–9 km), with similar WNW–ESE striking reverse mechanism. The timeline of the sequence suggests significant static stress interaction between the largest events. We perform here a detailed source inversion, first adopting a point source approximation and considering pure double couple and full moment tensor source models. We compare different extended source inversion approaches for the two largest events, and find that the rupture occurred in both cases along a subhorizontal plane, dipping towards SSW. Directivity is well detected for the May 20 main shock, indicating that the rupture propagated unilaterally towards SE. Based on the focal mechanism solution, we further estimate the co-seismic static stress change induced by the May 20 event. By using the rate-and-state model and a Poissonian earthquake occurrence, we infer that the second largest event of May 29 was induced with a probability in the range 0.2–0.4. This suggests that the segment of fault was already prone to rupture. Finally, we estimate peak ground accelerations for the two main events as occurred separately or simultaneously. For the scenario involving hypothetical rupture areas of both main events, we estimate Mw = 6.3 and an increase of ground acceleration by 50 per cent. The approach we propose may help to quantify rapidly which regions are invested by a significant increase of the hazard, bearing the potential for large aftershocks or even a second main shock.

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: In the framework of the Dead Sea Integrated Research project (DESIRE), 59 seismological stations were deployed in the region of the Dead Sea Basin. Twenty of these stations recorded data of sufficiently high quality between May and September 2007 to be used for ambient seismic noise analysis. Empirical Green’s functions are extracted from cross-correlations of long term recordings. These functions are dominated by Rayleigh waves, whose group velocities can be measured in the frequency range from 0.1 to 0.5 Hz. Analysis of positive and negative correlation lags of the Green’s functions makes it possible to identify the direction of the source of the incoming energy. Signals with frequencies higher than 0.2 Hz originate from the Mediterranean Sea, while low frequencies arrive from the direction of the Red Sea. Travel times of the extracted Rayleigh waves were measured between station pairs for different frequencies, and tomographically inverted to provide independent velocity models. Four such 2D models were computed for a set of frequencies, all corresponding to different sampling depths, and thus together giving an indication of the velocity variations in 3D extending to a depth of 10 km. The results show low velocities in the Dead Sea Basin, consistent with previous studies suggesting up to 8 km of recent sedimentary infill in the Basin. The complex structure of the western margin of the Basin is also observed, with sedimentary infill present to depths not exceeding 5 km west of the southern part of the Dead Sea. The high velocities associated with the Lisan salt diapir are also observed down to a depth of ~5 km. The reliability of the results is confirmed by checkerboard recovery tests.

Description: The Sumatra-Andaman earthquake with a magnitude of 9.3, and the subsequent destructive tsunami which caused more than 225 000 fatalities in the region of the Indian Ocean, happened on 26 December 2004. Less than one month later, the United Nations (UN) World Conference on Disaster Reduction took place in Kobe, Japan to commemorate the 1995 Kobe earthquake. The importance of preparedness and awareness on regional, national and community levels with respect to natural disasters was discussed during this meeting, and resulted in the approval of the Hyogo Declaration on Disaster Reduction. Based on this declaration the UN mandated the Intergovernmental Oceanographic Commission (IOC) of UNESCO (United Nations Education, Science and Cultural Organization), taking note of its over 40 years of successful coordination of the Pacific Tsunami Warning System (PTWC), to take on the international coordination of national early-warning efforts for the Indian Ocean and to guide the process of setting up a Regional Tsunami Early Warning System for the Indian Ocean.

Description: Trees, as long living plants, are governed by environmental and/or climate changes within their habitat. Their growth rings record to a large extent the temporal dynamics of these changes either directly or through tree physiological reactions. They render the highest time resolution thus far possible for environmental or climate reconstructions of the past 10,000 years (exactly dated, annually resolved, see: [828], [994]). Trees are a substantial part of the human environment with a high socio-economic value. Their large geographical extension over various regions of the world, including those with greatest population densities but also marginal areas allows to gain unique informations about local and regional consequences of global climate change.

Description: Results obtained from S and P receiver functions produced a clear image of the top and bottom of the subducting Nazca lithosphere beneath northern Chile. Using data from the teleseismic events recorded at 15 permanent Integrated Plate Boundary Observatory Chile (IPOC) stations, we obtained new constraints on the geometry and thickness of the descending Nazca lithosphere. We observed the subducted crust of the Nazca plate at depths ranging from 50 km beneath the Coastal Cordillera down to 110 km beneath the Western Cordillera. We found significant along‐strike variations in the geometry of the Nazca plate beneath northern Chile. On closer inspection, it appears that the oceanic Nazca plate is divided into two distinct segments as it descends beneath the continental South American plate. The transition from the relatively steeper (∼23°) and deeper slab to the north of 21°S to the flatter southern segment (∼19°) is shown reasonably clearly by our data. This feature could well be associated with variations in the curvature of the plate margin and the geometry of the Chile trench, which is mainly curved to the north of 21°S. We have also mapped the continental Moho of the South American plate at depths ranging between 60 and 70 km to the east of the Longitudinal Valley. Beneath the Coastal Cordillera, this boundary becomes invisible, probably due to the serpentinization of the forearc mantle wedge that reduces the velocity in the uppermost mantle. The base of the subducted Nazca plate was clearly identified as a sharp boundary in the results obtained from the P and S receiver functions. The thickness of the subducted oceanic Nazca plate, which has an age of ∼50 My, is estimated to be ∼50 km. Although this thickness is consistent with that predicted by thermal gradients, the explanation of the sharpness of the lithosphere‐asthenosphere boundary may require another mechanism such as hydration or melting.

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: An early detection of the presence of rupture directivity plays a major role in the correct estimation of ground motions and risks associated to the earthquake occurrence. We present here a simple method for a fast detection of rupture directivity, which may be additionally used to discriminate fault and auxiliary planes and have first estimations of important kinematic source parameters, such as rupture length and rupture time. Our method is based on the inversion of amplitude spectra from P-wave seismograms to derive the apparent duration at each station and on the successive modelling of its azimuthal behaviour. Synthetic waveforms are built assuming a spatial point source approximation, and the finite apparent duration of the spatial point source is interpreted in terms of rupture directivity. Since synthetic seismograms for a point source are calculated very quickly, the presence of directivity may be detected within few seconds, once a focal mechanism has been derived. The method is here first tested using synthetic datasets, both for linear and planar sources, and then successfully applied to recent Mw 6.2–6.8 shallow earthquakes in Peloponnese, Greece. The method is suitable for automated application and may be used to improve kinematic waveform modelling approaches.

Description: Multibeam bathymetry revealed the occurrence of numerous craterlike depressions, so-called pockmarks, on the sea floor of the Hammerfest Basin and the Loppa High, south-western Barents Sea. To investigate whether these pockmarks are related to ongoing gas seepage, microbial processes associated with methane metabolism were analyzed using geochemical, biogeochemical and microbiological techniques. Gravity cores were collected along transects crossing individual pockmarks, allowing a direct comparison between different locations inside (assumed activity center), on the rim, and outside of a pockmark (reference sites). Concentrations of hydrocarbons in the sediment, particularly methane, were measured as headspace (free) gas, and in the occluded and adsorbed gas fraction. Down to a depth of 2.6 m below sea floor (mbsf) sulfate reduction rates were quantified by radiotracer incubations. Concentrations of dissolved sulfate in the porewater were determined as well. Neither the sulfate profiles nor the gas measurements show any evidence of microbial activity or active fluid venting. Methane concentrations and sulfate reduction rates were extremely low or even below the detection limit. The results show that the observed sediment structures are most likely paleo-pockmarks, their formation probably occurred during the last deglaciation.

Description: On 2001 May 7, following unintentional water injection, a moderate size induced earthquake struck the Ekofisk oil field, North Sea. Despite of its relatively moderate magnitude, clear low-frequency waveforms could be recorded up to more than 2000 km epicentral distance, suggesting a slow rupture at very shallow depth and wave propagation through low-velocity shallow structures. The event poses a rare opportunity to constrain rupture velocity, duration and rise time of a superficial M 〉 4 event occurring on a horizontal plane in soft, water-saturated sediments. Two previous studies discussed the earthquake point source finding vertical dip-slip focal mechanisms with opposite senses of P and T axes. A further investigation was thus required to provide a basis for a deeper discussion of the failure dynamics. We significantly improve the used data set, test different earth models and derive a point source as well as a kinematic rupture model. We carefully discuss parameter uncertainties and effects related to shallow sources and wave propagation through different crustal structures to resolve the previous controversy. We additionally provide a kinematic rupture model, based on apparent source times derived from Rayleigh and Love waves. The waveforms resolve a predominant unilateral rupture along a horizontal plane at about 2 km depth. We derive an unusually slow rupture, consequence of a slow rupture velocity of about 500 m s –1 and a long rise time of about 7 s. An independent modelling of GPS- based static displacements allows to confirm the focal mechanism polarity and to locate the centroid at the eastern side of the field, resulting in a much larger seismic moment in comparison with dynamic seismic moment. The rupture directivity is confirmed by the relative location of the centroid with respect to the epicentre, which is set at the site of water injection.

Description: The Dead Sea Transform (DST) is a major left-lateral strike-slip fault that accommodates the relative motion between the African and Arabian plates, connecting a region of extension in the Red Sea to the Taurus collision zone in Turkey over a length of about 1100 km. The Dead Sea Basin (DSB) is one of the largest basins along the DST. The DSB is a morphotectonic depression along the DST, divided into a northern and a southern sub-basin, separated by the Lisan salt diapir. We report on a receiver function study of the crust within the multidisciplinary geophysical project, DEad Sea Integrated REsearch (DESIRE), to study the crustal structure of the DSB. A temporary seismic network was operated on both sides of the DSB between 2006 October and 2008 April. The aperture of the network is approximately 60 km in the E—W direction crossing the DSB on the Lisan peninsula and about 100 km in the N—S direction. Analysis of receiver functions from the DESIRE temporary network indicates that Moho depths vary between 30 and 38 km beneath the area. These Moho depth estimates are consistent with results of near-vertical incidence and wide-angle controlled-source techniques. Receiver functions reveal an additional discontinuity in the lower crust, but only in the DSB and west of it. This leads to the conclusion that the internal crustal structure east and west of the DSB is different at the present-day. However, if the 107 km left-lateral movement along the DST is taken into account, then the region beneath the DESIRE array where no lower crustal discontinuity is observed would have lain about 18 Ma ago immediately adjacent to the region under the previous DESERT array west of the DST where no lower crustal discontinuity is recognized.

Description: The Van (Eastern Anatolia, Turkey) earthquake occurred on Sunday, October 23, 2011 with a moment magnitude of 7.2. The tectonics of this region is characterized by strike–slip faulting on the Bitlis Suture Zone, and thrusting in the Zagros fold and thrust belt. Using high-rate (1 second) GPS data from permanent GNSS stations from the CORS-TR network, co-seismic displacements of eleven stations were determined using precise point positioning during this earthquake. We used the time series of coordinate changes for fourteen CORS-TR stations, and calculated the crust movements before and after the earthquake. According to the PPP solutions computed using high frequency GPS data to determine the co-seismic motions of stations, we conclude for the Van earthquake an occurrence time of 10:41:22 (UTC). No pre-seismic horizontal movement of stations at the level more than 5 mm before the earthquake could be observed. That means that no kinematic warning or prediction before the earthquake exists. Along an east–west horizontal line north of the Van Sea with a length of about 100 km, the northern part of this line experienced extension of 0.2–1 ppm in a NW–SE direction. The southern part experienced N–S shortening of 0.5–1.5 ppm. The N–S shortening we estimated geodetically matches well with the N–S shortening and thrust focal mechanism derived independently using seismic data by the USGS. Co-seismic surface displacements derived from the GPS data are consistent with the teleseismic source model given by the USGS. The geodetic source model derived from the GPS data reproduces the same moment magnitude and centroid as the teleseismic model, but shows a higher spatial resolution of the slip distribution. We also analyzed the post-seismic surface displacements derived from the GPS data within the first two weeks after the mainshock. No reasonable slip distribution on the co-seismic fault plane could be found, indicating that the sources for the early post-seismic deformation might come from the widely scattered aftershocks.

Description: Rapid improvements in telemetry technology and the general decrease in communication costs have raised a growing interest in low-cost wireless sensing units. This is especially the case for structural monitoring purposes, where they are becoming a more valuable alternative to conventional wired monitoring system. The main advantages associated with the use of wireless sensing unit include a considerable decrease in installation costs, decentralization of data analysis, and the possibility of broadening the functional capabilities by exploiting the use, at the same time and place, of different sensors. In this work, the design of a low-cost wireless sensing unit able both to collect, analyze, store, and communicate data and estimated parameters is presented. The suitability of a network of these low-cost wireless instruments for monitoring the vibration characteristics and dynamic properties of strategic civil infrastructures is validated during a ambient vibration recording field test on the Fatih Sultan Mehmet Bridge in Istanbul, Turkey.

Description: In 2004 and 2005 a passive seismic experiment was carried out in the northern and northeastern part of the Bohemian Massif (Sudetes) to study the lithospheric structure. We present results from Ps and Sp receiver function analyses. With one exception, Moho depth at stations in the northwestern part of the study area varies between 28 and 32 km. Thicker crust up to 35 km was mapped toward the south (Moldanubian unit) and toward the east (Moravo–Silesian and Brunovistulian units) confirming results from previous active seismic measurements. There exists a relatively sharp step in Moho depth between units of the central Sudetes (~ 30 km) and the Moravo–Silesian unit (~ 35 km). The vp/vs ratios inverted from primary and multiple Moho Ps conversions hint for different crustal compositions of the units. Toward the Carpathian thrust we have no clear indications for any crustal root or slab beneath the western Carpathians. However, our data suggests a deepening of the Moho or at least a complicated crust–mantle transition in this area. Additional Ps phases were observed between 6 and 10 s delay time in the Sudetes. These phases cannot be explained by Moho reverberations, but are most probably caused by low velocity zones in the middle crust or lithospheric mantle as shown by modeling of theoretical receiver functions. The stations showing these abnormal phases are located in the area of Permo-Carboniferous basins on probably Teplá–Barrandian crust. Therefore we assume that the phases hint at a mid-crustal low velocity zone between 16 and 20 km depth, which is interpreted as a felsic solidified magma reservoir of the Permo-Carboniferous volcanism beneath the Sudetic Basins. Sp receiver functions show phases with negative polarity at 9 to 12 s lead time on average, which we interpret as lithosphere–asthenosphere boundary at about 80 to 110 km depth.

Description: We use local earthquake data observed by the amphibious, temporary seismic MERAMEX array to derive spatial variations of seismic attenuation (Qp) in the crust and upper mantle beneath Central Java. The path-averaged attenuation values (t∗) of a high quality subset of 84 local earthquakes were calculated by a spectral inversion technique. These 1929 t∗-values inverted by a least-squares tomographic inversion yield the 3D distribution of the specific attenuation (Qp). Analysis of the model resolution matrix and synthetic recovery tests were used to investigate the confidence of the Qp-model. We notice a prominent zone of increased attenuation beneath and north of the modern volcanic arc at depths down to 15 km. Most of this anomaly seems to be related to the Eocene–Miocene Kendeng Basin (mainly in the eastern part of the study area). Enhanced attenuation is also found in the upper crust in the direct vicinity of recent volcanoes pointing towards zones of partial melts, presence of fluids and increased temperatures in the middle to upper crust. The middle and lower crust seems not to be associated with strong heating and the presence of melts throughout the arc. Enhanced attenuation above the subducting slab beneath the marine forearc seems to be due to the presence of fluids.

Description: For a long time the root mean square (RMS) error has been used in the EM community: - to characterize data fit for a particular model; - as a criterion to compare several models obtained from inversion. The RMS error appears to be a natural choice since we usually tackle inverse problems in a least-squares sense. Over the years, RMS became a customary criterion and gained ultimate significance. However, on the hunt for low RMS values, one often needs to introduce subjectivity by arbitrarily adjusting error floors or masking “bad” data without referring to the assumptions behind RMS. In this contribution, we revisit basic assumptions behind RMS, demonstrate its deficiency and propose alternative ways, which may provide more insight into our data and allow a more comprehensive assessment of the quality of the modelling result/resistivity model.

Description: We analyse data from seismic stations surrounding the Alboran Sea between Spain and North Africa to constrain variations of the lithosphere–asthenosphere boundary (LAB) in the region. The technique used is the receiver function technique, which uses S-to-P converted teleseismic waves at the LAB below the seismic stations. We confirm previous data suggesting a shallow (60–90 km) LAB beneath the Iberian Peninsula and we observe a similarly shallow LAB beneath the Alboran Sea where the lithosphere becomes progressively thinner towards the east. A deeper LAB (90–100 km) is observed beneath the Betics, the south of Portugal and Morocco. The structure of the LAB in the entire region does not seem to show any indication of subduction related features. We also observe good P receiver function signals from the seismic discontinuities at 410 and 660 km depth which do not indicate any upper-mantle anomaly beneath the entire region. This is in agreement with the sparse seismic activity in the mantle transition zone suggesting the presence of only weak and regionally confined anomalies.

Description: We determined a high-resolution 3-D S-wave velocity model for a 26 km × 12 km area in the northern part of the basin of Santiago de Chile. To reach this goal, we used microtremor recordings at 125 sites for deriving the horizontal-to-vertical (H/V) spectral ratios that we inverted to retrieve local S-wave velocity profiles. In the inversion procedure, we used additional geological and geophysical constraints and values of the thickness of the sedimentary cover already determined by gravimetric measurements, which were found to vary substantially over short distances in the investigated area. The resulting model was derived by interpolation with a kriging technique between the single S-wave velocity profiles and shows locally good agreement with the few existing velocity profile data, but allows the entire area, as well as deeper parts of the basin, to be represented in greater detail. The wealth of available data allowed us to check if any correlation between the S-wave velocity in the uppermost 30 m (v30S) and the slope of topography, a new technique recently proposed by Wald and Allen, exists on a local scale. We observed that while one lithology might provide a greater scatter in the velocity values for the investigated area, almost no correlation between topographic gradient and calculated v30S exists, whereas a better link is found between v30S and the local geology. Finally, we compared the v30S distribution with the MSK intensities for the 1985 Valparaiso event, pointing out that high intensities are found where the expected v30S values are low and over a thick sedimentary cover. Although this evidence cannot be generalized for all possible earthquakes, it indicates the influence of site effects modifying the ground motion when earthquakes occur well outside of the Santiago basin.

Description: We implement the effects of gravitational self-attraction and loading (SAL) into a global baroclinic ocean circulation model and investigate effects on sea level patterns, ocean circulation, and density distributions. We compute SAL modifications as an additional force on the water masses at every time step by decomposing the field of ocean bottom pressure anomalies into spherical harmonic functions and then applying Love numbers to account for the elastic properties of the solid Earth. Considering SAL in the postprocessing turns out to be insufficient, especially in coastal waters and on subweekly time scales, where SAL modifies local sea level by around 0.6–0.8 cm on average; in the open ocean, changes mostly remain around 0.3 cm. Modifications of water velocities as well as of heat and salt distributions are modeled, yet they are small. Simple parameterizations of SAL effects currently used in a number of ocean circulation models suffer from the process's inhomogeneity in space and time. These parameterizations improve the modeled sea level patterns but fail to reproduce SAL impacts on circulation and density distributions. We therefore suggest to explicitly consider the full SAL effect in ocean circulation models, especially when investigating sea level variations faster than around 4 days.

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: Frequent landsliding is one of the greatest natural hazards facing the inhabitants of Central Asia's Fergana Basin and the surrounding mountain ranges. Active tectonics in the region is rapidly building the Tien Shan, one of the highest mountain ranges on Earth, and the extreme topographic relief promotes frequent landslide activity, which causes major losses of life and property. In southwestern Kyrgyzstan alone, on average 10 people die and seven houses are destroyed each year in these sudden and rapidly moving landslides.

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: A series of linked marine and land studies have recently targeted the Sumatra subduction zone, focusing on the 2004 and 2005 plate boundary earthquake ruptures in Indonesia. A collaborative research effort by scientists from the United Kingdom (UK Sumatra Consortium), Indonesia, United States, France, and Germany is focusing on imaging the crustal structure of the margin to examine controls on along-strike and updip earthquake rupture propagation. The fundamental science objective is to examine how margin architecture and properties control earthquake rupture location and propagation.

Description: Late Miocene to Quaternary volcanic rocks from the frontal arc to the back-arc region of the Central Volcanic Zone in the Andes show a wide range of δ11B values (+4 to −7 ‰) and boron concentrations (6 to 60 ppm). Positive δ11B values of samples from the volcanic front indicate involvement of a 11B-enriched slab component, most likely derived from altered oceanic crust, despite the thick Andean continental lithosphere, and rule out a pure crust-mantle origin for these lavas. The δ11B values and boron concentrations in the lavas decrease with increasing depth of the Wadati-Benioff Zone. This across-arc variation in δ11B values and decreasing B/Nb ratios from the arc to the back-arc samples are attributed to the combined effects of boron-isotope fractionation during progressive dehydration in the slab and a steady decrease in slab-fluid flux toward the back arc, coupled with a relatively constant degree of crustal contamination as indicated by similar Sr, Nd and Pb isotope ratios in all samples. Three-component mixing calculations for slab-derived fluid, the mantle wedge and the continental crust based on B, Sr and Nd isotope data indicate that the slab-fluid component dominates the boron composition of the fertile mantle and that the primary arc magmas were contaminated by an average addition of 15 to 30% crustal material. Modeling of fluid-mineral boron-isotope fractionation as a function of temperature shows that dehydration reactions liberate continuously changing fluid compositions from the slab during progressive subduction. A combination of a boron-isotope fractionation model and a temperature model for the Central Andean subduction zone fits the across-arc variation in δ11B and we conclude that the boron-isotope composition of arc volcanic rocks, especially in island arcs, is dominated by changing δ11B-composition of boron-rich slab-fluids during progressive dehydration. Owing to the decrease in slab-derived fluid flux crustal contamination becomes more important toward the back-arc. Because of the boron-isotope fractionation effect, across-arc variations in δ11B need not necessarily reflect different mixing proportions between boron derived from the slab-fluid and the mantle wedge.

Description: Structural features of volcanic and hydrothermal systems can be used to infer the location of magma chambers or productive geothermal areas. The Hengill volcanic triple-junction complex has a well-developed geothermal system, which is being exploited to extract hot fluids that are used for electrical power and heat production. In the framework of the I-GET project, a 4-month temporary seismological network including seven high-dynamic broadband instruments was deployed and 1D transient electromagnetic soundings (TEM) and 3D magnetotelluric (MT) surveys were performed to improve the understanding of the relationships between structural features, seismic activity and fluid production at the Hengill geothermal system. The MT and TEM data set are analysed elsewhere. The analysis of the seismological data set allowed the detection and classification of more than 600 earthquakes, among which long-period (LP) earthquakes were observed for the first time in this area. This work focuses first on a joint inversion for the 3D velocity structure and determination of the locations of the hypocentres from about 250 local volcano-tectonic earthquakes with clear P- and S-wave arrival times. The results confirm those from earlier tomography studies in this area. Integrating the seismic velocity and resistivity models in a semi-quantitative approach by cross-plotting the resistivity model with the velocity ratio VP/VS delineates a structural body with a high seismic velocity ratio and low resistivity that is interpreted as the main heat source of the geothermal system.

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: Ambient noise tomography is applied to the significant data resources now available across Tibet and surrounding regions to produce Rayleigh wave phase speed maps at periods between 6 and 50 s. Data resources include the permanent Federation of Digital Seismographic Networks, five temporary U.S. Program for Array Seismic Studies of the Continental Lithosphere (PASSCAL) experiments in and around Tibet, and Chinese provincial networks surrounding Tibet from 2003 to 2009, totaling ∼600 stations and ∼150,000 interstation paths. With such a heterogeneous data set, data quality control is of utmost importance. We apply conservative data quality control criteria to accept between ∼5000 and ∼45,000 measurements as a function of period, which produce a lateral resolution between 100 and 200 km across most of the Tibetan Plateau and adjacent regions to the east. Misfits to the accepted measurements among PASSCAL stations and among Chinese stations are similar, with a standard deviation of ∼1.7 s, which indicates that the final dispersion measurements from Chinese and PASSCAL stations are of similar quality. Phase velocities across the Tibetan Plateau are lower, on average, than those in the surrounding nonbasin regions. Phase velocities in northern Tibet are lower than those in southern Tibet, perhaps implying different spatial and temporal variations in the way the high elevations of the plateau are created and maintained. At short periods (〈20 s), very low phase velocities are imaged in the major basins, including the Tarim, Qaidam, Junggar, and Sichuan basins, and in the Ordos Block. At intermediate and long periods (〉20 s), very high velocities are imaged in the Tarim Basin, the Ordos Block, and the Sichuan Basin. These phase velocity dispersion maps provide information needed to construct a 3-D shear velocity model of the crust across the Tibetan Plateau and surrounding regions.

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: Amplitude ratio of 30 short-period conspicuous P5KP and PKPab phases from five intermediate depth or deep events in Fiji-Tonga recorded at European stations around 150° distance shows a mean value two to three times the ratio of the synthetic amplitudes obtained by the normal-mode theory (and ak135 model) or by full-wave theory (and PREM). There is a large variance in the results, also observed in five amplitude ratios from one event in Argentina observed at temporary stations in China around 156°. Global recordings of three major deep earthquakes in Fiji, Bonin, and Western Brazil observed at ASAR, WRA, and ZRNK arrays, at 59 North America stations and at six South Pole stations displayed conspicuous P4KP and PcP (or ScP) phases. The amplitude ratio values of P4KP vs P(S)cP are sometimes almost one order of magnitude larger than the corresponding values of the synthetics. In both cases, arrival times and slowness values (corrected for ellipticity and station elevation) at the distances up to 23° beyond the A cutoff point predicted by ray theory match both the synthetics, suggesting the observations are the AB branch of PmKP (m = 4, 5) around 1 Hz. In disagreement to ray theory, no reliable BC branch is observed neither on the recordings nor on the normal-mode synthetics. The high amplitude ratio values cannot be explained by realistic perturbations of the velocity or attenuation values of the global models in the proximity of the core-to-mantle boundary (CMB). We speculate that the focusing effects and/or strong scattering most likely associated to some anomalous velocity areas of the lowermost mantle are responsible for that. The results suggest limitations of the previous evaluations of the short-period attenuation in the outer core from PmKP amplitudes (m ≥ 3), irrespective of the fact that they are obtained by using ray theory, normal-mode or full-wave synthetics. Attempts to use PmKP arrival times in order to refine velocity structure in the proximity of CMB should be also regarded with care if the propagation times have been computed with ray theory.

Description: The territory of Lithuania and adjacent areas of the East European Craton have always been considered a region of low seismicity. Two recent earthquakes with magnitudes of more than 5 in the Kaliningrad District (Russian Federation) on 21 September 2004 motivated re-evaluation of the seismic hazard in Lithuania and adjacent territories. A new opportunity to study seismicity in the region is provided by the PASSEQ (Pasive Seismic Experiment) project that aimed to study the lithosphere–asthenosphere structure around the Trans-European Suture Zone. Twenty-six seismic stations of the PASSEQ temporary seismic array were installed in the territory of Lithuania. The stations recorded a number of local and regional seismic events originating from Lithuania and adjacent areas. This data can be used to answer the question of whether there exist seismically active tectonic zones in Lithuania that could be potentially hazardous for critical industrial facilities. Therefore, the aim of this paper is to find any natural tectonic seismic events in Lithuania and to obtain more general view of seismicity in the region. In order to do this, we make a manual review of the continuous data recorded by the PASSEQ seismic stations in Lithuania. From the good quality data, we select and relocate 45 local seismic events using the well-known LocSAT and VELEST location algortithms. In order to discriminate between possible natural events, underwater explosions and on-shore blasts, we analyse spatial distribution of epicenters and temporal distribution of origin times and perform both visual analysis of waveforms and spectral analysis of recordings. We show that the relocated seismic events can be grouped into five clusters (groups) according to their epicenter coordinates and origin and that several seismic events might be of tectonic origin. We also show that several events from the off-shore region in the Baltic Sea (at the coasts of the Kaliningrad District of the Russian Federation) are non-volcanic tremors, although the origin of these tremor-type events is not clear.

Description: Seismological models of upper mantle structure provide important constraints on the Earth"s convection system. Resolving the details of the upper mantle discontinuities is important for modelling the composition of the mantle and for understanding the effect that the discontinuities may have on mantle convection. Recently, numerous permanent and temporary seismic stations and networks have been set up around the world. It is possible to get the seismic records for the research needs from data management systems like IRIS, GEOFON, GEOSCOPE, FREESIA, etc. The use of seismograms collected from a large number of stations and earthquakes around the world enable us to study the global and the regional structure of the Earth. In this work, the receiver function technique (e.g. Owens et. Al., 1995) is applied to study the upper mantle structure in the northwest Pacific subduction zone and in the Hawaiian hotspot area. In the northwest Pacific, the Pacific plate is subducted into the upper mantle to more than 600 km depth, indicated by seismicity. In Hawaii, the volcanic edifice of the Hawaiian Islands and seamounts are believed to result from the passage of the oceanic lithosphere over a stationary mantle hotspot (Wilson, 1963; Morgan, 1971; Morgan et. al., 1995). In both regions the upper mantle structure is affected by the cold and warm materials, respectively. To study the extension of the temperature anomaly is important for understanding the Earth"s convection system. The olivine component of the mantle material is intensively studied in laboratories (e.g. Ito and Takahashi, 1989; Irifune, 1987). With increasing temperature and pressure, the olivine crystal undergoes a series of phase transformations which will result in a variation of the seismic structure. The effect of the temperature anomaly on the main upper mantle discontinuities will be discussed in chapter 2. Recently, the receiver function technique is increasingly applied to investigate the upper mantle discontinuities. To isolate the upper mantle conversion phases, newly developed moveout correction and migration methods are applied to separately distributed seismic stations as well as station arrays. The receiver function method used in this study will be introduced in chapter 3. In chapter 4 and 5, receiver function studies in the area of the northwest Pacific subduction zone and the Hawaiian mantle plume are presented. Regional tectonic background and the previous seismological works in these two areas will be first introduced in each chapter, and followed by description of data, processing steps, results and interpretations. In chapter 6, I will summarize the observations of the 410 and 660 topography in the northwest Pacific subduction zone and in the area around the Hawaiian mantle plume.

Description: The volcanism responsible for creating the chain of the Hawaiian islands and seamounts is believed to mark the passage of the oceanic lithosphere over a mantle plume1,2. In this picture hot material rises from great depth within a fixed narrow conduit to the surface, penetrating the moving lithosphere3. Although a number of models describe possible plume–lithosphere interactions4, seismic imaging techniques have not had sufficient resolution to distinguish between them. Here we apply the S-wave ‘receiver function’ technique to data of three permanent seismic broadband stations on the Hawaiian islands, to map the thickness of the underlying lithosphere. We find that under Big Island the lithosphere is 100–110 km thick, as expected for an oceanic plate 90–100 million years old that is not modified by a plume. But the lithosphere thins gradually along the island chain to about 50–60 km below Kauai. The width of the thinning is about 300 km. In this zone, well within the larger-scale topographic swell, we infer that the rejuvenation model5 (where the plume thins the lithosphere) is operative; however, the largerscale topographic swell is probably supported dynamically.

Description: We examine shear-wave splitting of SKS waveforms collected by a temporary array of 68 stations in the region of the Dead Sea basin. The observed splitting parameters exhibit systematic variations along a dense, EW-trending 60 km profile across the basin. The delay times vary significantly between 1.0 and 2.8 seconds with smaller values in the very center of the profile. The fast polarizations are oriented more-or-less parallel to the strike of the Dead Sea transform fault and vary between −10 and 20 degrees with respect to North. Finite-frequency waveform modeling reveals that the source-region of the small-scale lateral variations is likely located within the crust. The modeling further shows that purely isotropic velocity variations affect shear-wave splitting: To a large degree, the observed variations of splitting parameters can be explained by the sedimentary fill of the basin and its low isotropic seismic velocities, whereas the mantle is uniformly anisotropic. Our study indicates that precaution must be taken when interpreting short-scale lateral variations of shear wave splitting in terms of anisotropic structures in the crust or upper mantle.

Description: In sequence of the GFZ93 high resolution gravity models (Gruber Th., et al, 1993) a new model, named GFZ95A, which is complete to degree and order 360 of a spherical harmonic series was computed. The model is based on new data sets, which were collected during the last months. This new data promises a major step towards a more precise high resolution gravity model. Especially from new available data over CIS (Community of Independent States) major progress can be expected with respect to the former models, which were based on predicted data in this area.

Description: We performed a teleseismic P wave tomography study using seismic events at both teleseismic and regional distances, recorded by a temporary seismic array in the Argentine Puna Plateau and adjacent regions. The tomographic images show the presence of a number of positive and negative anomalies in a depth range of 20–300 km beneath the array. The most prominent of these anomalies corresponds to a low-velocity body, located in the crust, most clearly seen in the center of the array (27°S, 67°W) between the Cerro Peinado volcano, the Cerro Blanco caldera and the Farallon Negro in the east. This anomaly (southern Puna Magmatic Body) extends from the northern most part of the array and follows the line with the highest density of stations towards the south where it becomes smaller. It is flanked by high velocities on the west and the east respectively. On the west, the high velocities might be related to the subducted Nazca plate. On the northeast the high velocity block coincides with the position of the Hombre Muerto basin in the crust and could be indicating an area of lithospheric delamination where we detected a high velocity block at 100 km depth on the eastern border of the Puna plateau, north of Galan. This block might be related to a delamination event in an area with a thick crust of Paleozoic metamorphic rocks at the border between Puna and Eastern Cordillera. In the center of the array the Southern Puna magmatic body is also flanked by high velocities but the most prominent region is located on the east and is interpreted as part of the Sierras Pampeanas lithosphere with high velocities. The position of the Sierras Pampeanas geological province is key in this area as it appears to limit the extension of the plateau towards the south.

Description: The transition between the lithosphere and the asthenosphere is subject to numerous contemporary studies as its nature is still poorly understood. The thickest lithosphere is associated with old cratons and platforms and it has been shown that seismic investigations may fail to image the lithosphere-asthenosphere boundary in these areas. Instead, several recent studies have proposed a mid-lithospheric discontinuity of unknown origin existing under several cratons. In this study we investigate the Tanzania craton in East Africa which is enclosed by the eastern and western branches of the East African Rift System. We present evidence from S receiver functions for two consecutive discontinuities at depths of 50–100 km and 140–200 km, which correspond to significant S wave velocity reductions under the Tanzania craton and the Albert and Edward rift segments. By comparison with synthetic waveforms we show that the lower discontinuity coincides with the LAB exhibiting velocity reductions of 6–9%. The shallower interface reveals a velocity drop that varies from 12% beneath the craton to 24% below the Albert-Edward rift. It is interpreted as an infiltration front marking the upper boundary of altered lithosphere due to ascending asthenospheric melts. This is corroborated by computing S velocity variations based on xenolith samples which exhibit a dense system of crystallized veins acting as pathways of the infiltrating melt. Mineral assemblages in these veins are rich in phlogopite and pyroxenite which can explain the reduced shear wave velocities. Melt infiltration represents a suitable mechanism to form a mid-lithospheric discontinuity within cratonic lithosphere that is underlain by anomalously hot mantle.

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: Clear S-to-P converted waves from the crust–mantle boundary (Moho) and lithosphere–asthenosphere boundary (LAB) have been observed on the eastern part of the Dead Sea Basin (DSB), and are used for the determination of the depth of the Moho and the LAB. A temporary network consisting of 18 seismic broad-band stations was operated in the DSB region as part of the DEad Sea Integrated REsearch project for 1.5 years beginning in September 2006. The obtained Moho depth (∼35 km) from S-to-P receiver functions agrees well with the results from P-to-S receiver functions and other geophysical data. The thickness of the lithosphere on the eastern part of the DSB is about 75 km. The results obtained here support and confirm previous studies, based on xenolith data, geodynamic modeling, heat flow observations, and S-to-P receiver functions. Therefore, the lithosphere on the eastern part of the DSB and along Wadi Araba has been thinned in the Late Cenozoic, following rifting and spreading of the Red Sea. The thinning of the lithosphere occurred without a concomitant change in the crustal thickness and thus an upwelling of the asthenosphere in the study area is invoked as the cause of the lithosphere thinning.

Description: In the past decade the analysis of seismic noise has become an efficient tool to recover the Green's function between pairs of receivers by cross-correlation of seismic traces. Most studies focus on the investigation of the surface wave component of the ambient noise. Several attempts to recover the body wave part of the Green's function have been documented. In this paper I present the results of cross-correlation of seismic noise and the retrieval of refracted and reflected P-waves along a seismic line in the Karoo region (Republic of South Africa). Body wave refractions (direct phases) and reflections have been observed in the Green's functions derived from ambient noise records of up to 60 hours. The results are compared with shot gathers from a controlled source experiment (borehole explosions), carried out along the same line. The significant potential of ambient noise analysis, especially with respect to P-wave reflections will be shown and discussed.

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 magnitude-8.8 Maule (Chile) earthquake of 27 February 2010 ruptured a segment of the Andean subduction zone megathrust that has been suspected to be of high seismic potential1,2,3,4,5,6. It is the largest earthquake to rupture a mature seismic gap in a subduction zone that has been monitored with a dense space-geodetic network before the event. This provides an image of the pre-seismically locked state of the plate interface of unprecedentedly high resolution, allowing for an assessment of the spatial correlation of interseismic locking with coseismic slip. Pre-seismic locking might be used to anticipate future ruptures in many seismic gaps6,7,8,9,10,11,12, given the fundamental assumption that locking and slip are similar. This hypothesis, however, could not be tested without the occurrence of the first gap-filling earthquake. Here we show evidence that the 2010 Maule earthquake slip distribution correlates closely with the patchwork of interseismic locking distribution as derived by inversion of global positioning system (GPS) observations during the previous decade. The earthquake nucleated in a region of high locking gradient and released most of the stresses accumulated in the area since the last major event in 1835. Two regions of high seismic slip (asperities) appeared to be nearly fully locked before the earthquake. Between these asperities, the rupture bridged a zone that was creeping interseismically with consistently low coseismic slip. The rupture stopped in areas that were highly locked before the earthquake but where pre-stress had been significantly reduced by overlapping twentieth-century earthquakes. Our work suggests that coseismic slip heterogeneity at the scale of single asperities should indicate the seismic potential of future great earthquakes, which thus might be anticipated by geodetic observations.

Description: Local earthquake data from a dense temporary seismological network in the southern Dead Sea area have been analysed within the project DESIRE (Dead Sea Integrated Research Project). Local earthquakes are used for the first precise image of the distribution of the P-wave velocity and the vP/vS ratios. 65 stations registered 655 local events within 18 months of observation time. A subset of 530 well-locatable events with 26 730 P- and S-arrival times was used to calculate a tomographic model for the vP and vP/vS distribution. Since the study area is at first-order 2-D, a gradual approach was chosen, which compromised a 2-D inversion followed by a 3-D inversion. The sedimentary basin fill is clearly imaged through high vP/vS ratios and low vP. The basin fill shows an asymmetric structure with average depth of 7 km at the western boundary and depth between 10 and 14 km at the eastern boundary. This asymmetry is reflected by the vertical strike-slip eastern border fault, and the normal faulting at the western boundary, caused by the transtensional deformation within the last 5 Myr. Within the basin fill the Lisan salt diapir is imaged through low vP/vS ratios, reflecting its low fluid content. The extensions were determined to 12 km in E–W and 17 km in N–S direction while its depth is 5–6 km. The thickness of the pre-basin sediments below the basin fill cannot be derived from the tomography data—it is estimated to less than 3 km from former investigations. Below the basin, down to 18 km depth very low P-wave velocities and low vP/vS ratios are observed—most likely caused by fluids from the surrounding crust or the upper mantle.

Description: For the first time, ambient noise tomography is used to clearly image the magma chamber beneath Lake Toba caldera, one of the largest Quaternary calderas on Earth. Using data from 40 seismic stations deployed between May and October 2008 around Lake Toba, empirical Green's functions are extracted from long term cross-correlations of continuous records. These functions are dominated by Rayleigh waves, whose group velocities can be measured in the period range from 2.5 to 12 seconds. Arrival times of these waves are picked for a given period and inverted using 2-D tomography to calculate lateral variations in velocity for the given period. This was done for six different periods, which all correspond to different sampling depths. Thus the six 2-D models presented together provide information on velocity variations with depth. The results show a low-velocity body coincident with the Lake Toba caldera, representing the magma chamber under the volcano. The chamber is observed to have a complex 3-D geometry, with at least two separate sub-chambers underlying the caldera. Other results include a deep low velocity body, possibly another magma chamber, south west of the lake with an upper limit of ∼7 km depth. The maximum depth to which this body reaches could not be resolved. The Sumatra Fault marks a velocity contrast, but only down to depths not greater than 5 km. The reliability of the results was further confirmed by checkerboard recovery tests.

Description: During the TOR-1 passive seismic experiment in 1996/97, a maximum of 139 temporary seismograph stations were operating over the Sorgenfrei-Tornquist Zone (STZ) in an area extending from northern Germany through Denmark to central Sweden. One of the objectives was to study horizontal anisotropy directions in the subcrustal lithosphere and asthenosphere across the Trans-European Suture Zone. To achieve this goal,broad-band and intermediate-period (5 s) data of the TOR-1 stations and additional stations of permanent networks (GRSN, GEOFON) were analysed for splitting of SKS and SKKS phases. As a result of the relatively dense station spacing, the method offers good lateral resolution of anisotropy.Preliminary results suggest that the directions of the fast horizontal S wave velocity are affected by the STZ. In central Europe and southern Sweden, far away from the STZ, fast S wave directions are approximately E-W while they turn more northerly closer to the STZ where they are approximately parallel to the trend of the STZ. No significant shear wave splitting was observed north of 57 degr. N and east of 14 degr. E. Small delay times between 0.2 and 0.5 s observed at the northernmost TOR-1 station T40S and T60S may be controlled by anisotropy in a thickened crust. The mantle contribution of horizontal anisotropy within the STZ is probably constrained to an approximately 60-km-thick zone in the depth range between 70 and 300 km. The observations are consistent with a model where azimuthally anisotropy is not governed by present-day mantle flow in the asthenosphere, but rather is frozen into the subcrustal lithosphere during the last episode of tectonic activity.

Description: The plate tectonic concept has nearly revolutionized the geo-sciences in the last decades. Many geodynamical, geophysical and geological phenomena could be understood for the first time. The investigation of the present-day global crustal kinematics is a crucial contribution which can be delivered by space geodesy only. Since the beginning of the MERIT-Campaign continuous good quality laser ranging data to LAGEOS are available. Using these data a precise estimation of present-day global-scale tectonic motion is possible. During the WEGENER/MEDLAS project the worldwide tracking network has been densified by a number of fixed and mobile stations in the southern Europe and eastern Mediterranean. Having now for a number of sites three to four occupations within four years significant site motions can be extracted. The derived motions are essential input information for deformation analyses and the interpretation of plate tectonic parameters.

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: We investigate microseismic activity at the convergent plate boundary of the Hellenic subduction zone on- and offshore south-eastern Crete with unprecedented precision using recordings from an amphibian seismic network. The network configuration consisted of up to eight ocean bottom seismometers as well as five temporary short-period and six permanent broadband stations on Crete and surrounding islands. More than 2,500 local and regional events with magnitudes up to M L = 4.5 were recorded during the time period July 2003–June 2004. The magnitude of completeness varies between 1.5 on Crete and adjacent areas and increases to 2.5 in the vicinity of the Strabo trench 100 km south of Crete. Tests with different localization schemes and velocity models showed that the best results were obtained from a probabilistic earthquake localization using a 1-D velocity model and corresponding station corrections obtained by simultaneous inversion. Most of the seismic activity is located offshore of central and eastern Crete and interpreted to be associated with the intracrustal graben system (Ptolemy and Pliny trenches). Furthermore, a significant portion of events represents interplate seismicity along the NNE-ward dipping plate interface. The concentration of seismicity along the Ptolemy and Pliny trenches extends from shallow depths down to the plate interface and indicates active movement. We propose that both trenches form transtensional structures within the Aegean plate. The Aegean continental crust between these two trenches is interpreted as a forearc sliver as it exhibits only low microseismic activity during the observation period and little or no internal deformation. Interplate seismicity between the Aegean and African plates forms a 100-km wide zone along dip from the Strabo trench in the south to the southern shore-line of Crete in the north. The seismicity at the plate contact is randomly distributed and no indications for locked zones were observed. The plate contact below and north of Crete shows no microseismic activity and seems to be decoupled. The crustal seismicity of the Aegean plate in this area is generally confined to the upper 20 km in agreement with the idea of a ductile deformation of the lower crust caused by a rapid return flow of metamorphic rocks that spread out below the forearc. In the region of the Messara half-graben at the south coast of central Crete, a southward dipping seismogenic structure is found that coalesces with the seismicity of the Ptolemy trench at a depth of about 20 km. The accretionary prism south of Crete indicated by the Mediterranean Ridge showed no seismic activity during the observation period and seems to be deforming aseismically.

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 present a series of 1-D shear velocity models for the Sumatran Forearc and Arc derived from Rayleigh wave group dispersion in noise correlation functions from vertical and pressure records from an onshore–offshore seismic deployment. The 1-D models represent the crustal structure of the downgoing Indian Plate, the accretionary prism and the arc. There is a progression in shear velocity across the forearc to the arc associated with thickening of the accretionary prism and the development of an arc crust. The velocity structure inferred for the upper 20 km based on path averages between stations on the accretionary prism has velocities consistent with a thick sediment package in agreement with estimates of depth to the plate boundary determined from active source experiments. We also find low Indian Plate shear velocities, 〈4 km s−1 to 25 km depth beneath our station locations on the downgoing plate. These low seismic velocities are consistent with at least 14–24 per cent serpentinization of the oceanic crust and upper mantle of the downgoing plate. This high degree of serpentinization, may weaken the plate interface and explain the segmentation observed in the great Sumatran thrust earthquakes if the serpentinization is localized. The success of this study suggests that future onshore–offshore seismic deployments will be able to utilize this method.

Description: The Mw = 9.3 Sumatra earthquake of 26 December 2004 generated a tsunami that affected the entire Indian Ocean region and caused approximately 230 000 fatalities. In the response to this tragedy the German government funded the German Indonesian Tsunami Early Warning System (GITEWS) Project. The task of the GEOFON group of GFZ Potsdam was to develop and implement the seismological component. In this paper we describe the concept of the GITEWS earthquake monitoring system and report on its present status. The major challenge for earthquake monitoring within a tsunami warning system is to deliver rapid information about location, depth, size and possibly other source parameters. This is particularly true for coast lines adjacent to the potential source areas such as the Sunda trench where these parameters are required within a few minutes after the event in order to be able to warn the population before the potential tsunami hits the neighbouring coastal areas. Therefore, the key for a seismic monitoring system with short warning times adequate for Indonesia is a dense real-time seismic network across Indonesia with densifications close to the Sunda trench. A substantial number of supplementary stations in other Indian Ocean rim countries are added to strengthen the teleseismic monitoring capabilities. The installation of the new GITEWS seismic network – consisting of 31 combined broadband and strong motion stations – out of these 21 stations in Indonesia – is almost completed. The real-time data collection is using a private VSAT communication system with hubs in Jakarta and Vienna. In addition, all available seismic real-time data from the other seismic networks in Indonesia and other Indian Ocean rim countries are acquired also directly by VSAT or by Internet at the Indonesian Tsunami Warning Centre in Jakarta and the resulting "virtual" network of more than 230 stations can jointly be used for seismic data processing. The seismological processing software as part of the GITEWS tsunami control centre is an enhanced version of the widely used SeisComP software and the well established GEOFON earthquake information system operated at GFZ in Potsdam (http://geofon.gfz-potsdam.de/db/eqinfo.php). This recently developed software package (SeisComP3) is reliable, fast and can provide fully automatic earthquake location and magnitude estimates. It uses innovative visualization tools, offers the possibility for manual correction and re-calculation, flexible configuration, support for distributed processing and data and parameter exchange with external monitoring systems. SeisComP3 is not only used for tsunami warning in Indonesia but also in most other Tsunami Warning Centres in the Indian Ocean and Euro-Med regions and in many seismic services worldwide.

Description: We analyse broad-band SS waveformdata recorded by several networks in Europe with sources mainly in the west Pacific to study the underside reflections of teleseismic SS waves in the lithosphere and the upper mantle beneath eastern Asia and the NW Pacific ocean. SS bounce points sample a corridor from the Aleutian, Kamchatka and Japan subduction zones through the North China Craton and Central Asian Orogenic Belt to the Tibetan plateau. The corridor passes through different tectonic units such as subduction zones, an old continental shield, a fold belt and a high plateau. We investigate the seismic structure of the lithosphere and the mantle transition zone beneath the different geotectonic units along the profile and infer the correlation of geodynamic processes at different depths.We explore the short period frequency content in the SS waveform data and use moveout correction and common midpoint stack to acquire profiles with high lateral and depth resolution from the crust to the mantle transition zone. Clear SS precursors of the 410 and 660 km discontinuities show the effects of the interaction between the subducted oceanic lithosphere and the mantle transition zone beneath the NW Pacific subduction zones. A low-velocity layer has also been detected beneath the 410 km discontinuity and can be traced along the entire profile. Due to the improved resolution acquired by the method presented here we have been able to study the shallower structures such as the Moho and the lithosphere–asthenosphere boundary by SS precursors. The continental Moho can be clearly seen along this corridor. The depth variation agrees well with earlier receiver function results. We also see negative reflectors along the profile at varying depths, which can be interpreted as the lithosphere–asthenosphere boundary.

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: Noise measurements were carried out at 381 sites in the Cologne area (Germany) using both short period and broad band sensors. The large number of data allowed both assesment of the influence of different sensors in the site response estimation and to compare the widely used H/V technique with the recently proposed Fourier Phase Spectral Method (FPSM). The results show that short period sensors are able to reliably retrieve site effects at frequencies well below their corner frequencies. Moreover, the H/V method should be preferred to the FPSM in determining the fundamental resonance frequency of soils. Finally, a map showing the resonance frequency distribution in the studied area was drawn using the results obtained applying the H/V technique.

Description: Earthquake early warning systems should provide reliable warnings as quickly as possible with a minimum number of false and missed alarms. Wireless meshed networks, coupled with low-cost seismometers for monitoring, evaluation, and information about seismic vibrations in space and time are introducing a new generation of warning infrastructures for mega-cities. The use of a cooperative method for signal analysis makes it possible to distinguish earthquakes (with a certain minimal magnitude) from other ground shaking in a city. The paper gives a short overview of our approach for developing decentralized early warning systems and an evaluation based on experiences gained from model investigations, testbeds in Berlin, and prototype installations in Istanbul.

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: Two passive seismic experiments have been carried out across the Trans European Suture Zone (TESZ) from northern Germany to southern Sweden (TOR) and across the Proterozoic-Archaean suture in Finland (SVEKALAPKO) to improve our understanding of the processes involved in the creation of the European continent. Teleseismic earthquakes recorded by the studies of the crust-mantle, and upper mantle seismic discontinuities with the receiver function method. Along the TOR network the depth to the Moho increases from 30 km at the southern edge of the profile to 40 km at the Elbe Line. Between the Elbe Line and TESZ the Moho branches off and whereas the deeper branch continues at 40 km depth to the TESZ a second branch appears at 30.35 km depth. The upper branch descends north of the TESZ to below 55 km under the northern end of the TOR profile. The crustal thickening north of the TESZ is accompanied by an increase in average Vp/Vs values, appearance of intracrustal conversion zones and north dipping features which we interpret as remnants of the subduction and subsequent collision between Avalonia and Baltica. In southern Finland beneath the SVEKALAPKO network the Moho starts in the south at the depth of 40-45 km, plunges to about 65 km depth south of the Archaean-Proterozoic suture. This deepening of the Moho is coincident with a north dipping intracrustal structure apparently related to the subduction and collision and of the Proterozoic and Archaean provinces in Proterozoic. North of the line of the suture the Moho rises smoothly to 45-50 km depth in the Archaean province. Along the TOR profile, 410 and 660 discontinuities were hard to detect. However, manyfold stacking of receiver functions revealed that the conversions from the two discontinuities arrive more or less accordingto IASP91 predicted time. Across the SVEKALAPKO network 410 and 660 discontinuities arrive markedly earlier than IASP91 theoretical arrival times. In particular north of the Archaean-Proterozoic suture in Finland the 410 and 660 km conversions arrive about 2s earlier, indicating about 5 per cent higher average upper mantle velocities and lower temperatures than what IASP91 global model predicts. Test

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: S receiver functions obtained from seismograms of teleseismic events recorded at 78 European permanent broad-band stations are used to estimate the thickness of the European lithosphere. Our results provide new, independent information about the lithospheric thickness beneath the Precambrian platform of Eastern Europe and the Phanerozoic platform of central Europe. Detailed high-resolution images of the lithosphere–asthenosphere boundary (LAB) reveal indications for a typical continental lithosphere of about 100 km thickness beneath a majority of stations within Central Europe, whereas in the vicinity of the Trans-European Suture Zone (TESZ), the lithosphere thickens to about 130 km. A relatively thin lithosphere of 80 km was found beneath the Upper Rhine Graben region suggesting that the Cenozoic extension affects the whole lithosphere. No clear signal from the LAB was detected beneath the Alps and Carpathians. The LAB Sp phase might be disturbed by complicated structure due to ongoing collision/subduction in these regions, or the data are not yet sufficiently dense. A relatively thicker lithosphere of about 120 km was found beneath the SW part of the Bohemian Massif that was formed during the Variscan orogeny.We found an LAB depth of about 190 km near a single station located in the Vrancea area/Eastern Carpathians, which is characterized by the occurrence of intermediate deep earthquakes. Beneath the stations located in the Precambrian platform of Eastern Europe, the LAB deepens to approximately 〉200 km, even though the converted phase from the LAB is not as sharp as found beneath other stations located in Central Europe or even is missing.

Description: The GLATIS project (Greenland Lithosphere Analysed Teleseismically on the Ice Sheet) with collaborators has operated a total of 16 temporary broadband seismographs for periods from 3 months to 2 years distributed over much of Greenland from late 1999 to the present. The very first results are presented in this paper, where receiver-function analysis has been used to map the depth to Moho in a large region where crustal thicknesses were previously completely unknown. The results suggest that the Proterozoic part of central Greenland consists of two distinct blocks with different depths to Moho. North of the Archean core in southern Greenland is a zone of very thick Proterozoic crust with an average depth to Moho close to 48 km. Further to the north the Proterozoic crust thins to 37–42 km. We suggest that the boundary between thick and thin crust forms the boundary between the geologically defined Nagssugtoqidian and Rinkian mobile belts, which thus can be viewed as two blocks, based on the large difference in depth to Moho (over 6 km). Depth to Moho on the Archean crust is around 40 km. Four of the stations are placed in the interior of Greenland on the ice sheet, where we find the data quality excellent, but receiver-function analyses are complicated by strong converted phases generated at the base of the ice sheet, which in some places is more than 3 km thick.

Description: We use Global Positioning System (GPS) velocities and kinematic Finite Element models (FE-models) to infer the state of locking between the converging Nazca and South America plates in South-Central Chile (36[degree sign]S - 46[degree sign]S) and to evaluate its spatial and temporal variability. GPS velocities provide information on earthquake-cycle deformation over the last decade in areas affected by the megathrust events of 1960 (Mw= 9.5) and 2010 (Mw= 8.8). Our data confirm that a change in surface velocity patterns of these two seismotectonic segments can be related to their different stages in the seismic cycle: Accordingly, the northern (2010) segment was in a final stage of interseismic loading whereas the southern (1960) segment is still in a postseismic stage and undergoes a prolonged viscoelastic mantle relaxation. After correcting the signals for mantle relaxation, the residual GPS velocity pattern suggests that the plate interface accumulates slip deficit in a spatially and presumably temporally variable way towards the next great event. Though some similarity exist between locking and 1960 coseismic slip, extrapolating the current, decadal scale slip deficit accumulation towards the ~ 300-yr recurrence times of giant events here does neither yield the slip distribution nor the moment magnitude of the 1960 earthquake. This suggests that either the locking pattern is evolving in time (to reconcile a slip deficit distribution similar to the 1960 earthquake) or that some asperities are not persistent over multiple events. The accumulated moment deficit since 1960 suggests that highly locked patches in the 1960 segment are already capable of producing a M ~ 8 event if triggered to fail by stress transfer from the 2010 event.

Description: Determining the relative orientation of the horizontal components of seismic sensors is a common problem that limits data analysis and interpretation for several acquisition setups, including linear arrays of geophones deployed in borehole installations or ocean bottom seismometers deployed at the seafloor. To solve this problem we propose a new inversion method based on a complex linear algebra approach. Relative orientation angles are retrieved by minimizing, in a least-squares sense, the l2-norm between the complex traces (hodograms) of adjacent pairs of sensors. This methodology can be applied without restrictions only if the wavefield recorded by each pair of sensors is very similar. In most cases, it is possible to satisfy this condition by low-pass filtering the recorded waveforms. The main advantage of our methodology is that, in the complex domain, the relative orientations of seismic sensors can be viewed as a linear inverse problem, which ensures that the preferred solution corresponds to the global minimum of a misfit function. It is also possible to use simultaneously more than one independent data set (other seismic events) to better constrain the solution of the inverse problem. Furthermore, by a computational point of view, our method results faster than the relative orientation methods based on waveform cross-correlation. After several tests on synthetic data sets we applied successfully our methodology to different types of real data. These applications include the alignment of borehole sensors relative to a Vertical Seismic Profiling (VSP) acquisition and the orientation of Ocean Bottom Seismometers (OBS) relative to a neighbouring land station of known orientation. Using land stations, the absolute orientation of OBS can be retrieved. Finally, as a last application, we checked the correct orientation for land stations of a seismological array in Germany.

Description: An automatic procedure is presented to retrieve rupture parameters for large earthquakes along the Sunda arc subduction zone. The method is based on standard array analysis and broadband seismograms registered within 30°–100° epicentral distance. No assumptions on source mechanism are required. By means of semblance the coherency of P waveforms is analysed at separate large-aperture arrays. Waveforms are migrated to a 10°×10° wide source region to study the spatio-temporal evolution of earthquakes at each array. The multiplication of the semblance source maps resulting at each array increases resolution. Start, duration, extent, direction, and propagation velocity are obtained and published within 25 min after the onset of the event. First preliminary results can be obtained even within 16 min. Their rapid determination may improve the mitigation of the earthquake and tsunami hazard. Real-time application will provide rupture parameters to the GITEWS project (German Indonesian Tsunami Early Warning System). The method is applied to the two M8.0 Sumatra earthquakes on 12 September 2007, to the M7.4 Java earthquake on 2 September 2009, and to major subduction earthquakes that have occurred along Sumatra and Java since 2000. Obtained rupture parameters are most robust for the largest earthquakes with magnitudes M≥8. The results indicate that almost the entire seismogenic part of the subduction zone off the coast of Sumatra has been ruptured. Only the great Sumatra event in 2004 and the M7.7 Java event on 17 July 2006 could reach to or close to the surface at the trench. Otherwise, the rupturing was apparently confined to depths below 25 km. Major seismic gaps seem to remain off the coast of Padang and the southern tip of Sumatra.

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: The main structures of the mid-European lithosphere, crossing Europe from the Iberian peninsula to the Bohemian Massif, are predominantly formed by the Variscan orogeny. To investigate the anisotropy of the mantle at the transition zone between the two Variscan units, the Saxothuringicum and the Moldanubicum, we carried out a field experiment in SE Germany in 1995-1996: 23 mobile broad-band stations were installed for 6 months in the Vogtland-Oberpfalz-Bavarian Forest area. The station profile crossed the suture zone of the Saxothuringicum and Moldanubicum near the KTB borehole (German Continental Deep Drilling Program). With a mean station spacing of about 10 km? we intended to obtain a high lateral resolution of the anisotropy parameters and to resolve possible changes when passing the suture zone. The analysis of the observed birefringence of SKS phases shows E-W directions for the fast polarization. Therefore, the directions deviate only slightly from the strike of the Hercynian mountain belt and from the direction of the absolute plate motion in that region. Indications for the transition zone come from a rotation of the fast polarization direction from 86 degrees +/- 13 degrees in the northern part of the profile (Saxothuringicum) to 110 degrees +/- 15 degrees in the southern part (Moldanubicum) as well as from strong variations of the splitting parameters with respect to the azimuths of the incoming waves in the middle of the profile. We interpret these variations as an expression of a complex mantle structure formed either by several anisotropic layers with inclined symmetry axes in at least one layer or by a model consisting of inhomogeneous anisotropic layers. A comparison of the azimuthal variations of the splitting parameters in the middle of the profile with those observed at the Grafenberg station GRAl-situated in the central part of the transition zone approximately 100 km to the west-shows remarkable differences, which may reflect lateral variations in the direction and inclination of the symmetry axes in the transition zone even on a small scale. Both observations-the change in the fast polarization direction from the northern to the southern part as well as the variations with respect to different azimuths in the middle of the profile-suggest that the transition zone between the Saxothuringicum and the Moldanubicum continues down into the upper mantle.

Description: In this essay we provide a short introduction to and overview of the basics of stable isotope geochemistry and its common application in petroleum geochemistry. We identify the processes that are responsible for the carbon and hydrogen isotopic compositions of biological and geological organic matter and indicate the utility of stable isotopes in oil-source rock correlations. Stable isotope analyses are also exploited in the investigation of different alteration processes within oils and petroleum reservoirs. State of the art work is presented, and future research needs are identified.

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.