ALBERT

Description: The Central Andes of southern Peru, Bolivia, Argentina and Chile (between 12°S and 42°S) comprise the largest orogenic plateau in the world associated with abundant arc volcanism, the Central Andean Plateau, as well as multiple segments of flat-slab subduction making this part of the Earth a unique place to study various aspects of active plate tectonics. The goal of this continental-scale ambient noise tomography study is to incorporate broad-band seismic data from 20 seismic networks deployed incrementally in the Central Andes from 1994 May to 2012 August, to image the vertically polarized shear wave velocity (Vsv) structure of the South American Cordillera. Using dispersion measurements calculated from the cross-correlation of 330 broad-band seismic stations, we construct Rayleigh wave phase velocity maps in the period range of 8–40 s and invert these for the shear wave velocity (Vsv) structure of the Andean crust. We provide a dispersion misfit map as well as uncertainty envelopes for our Vsv model and observe striking first-order correlations with our shallow results (∼5 km) and the morphotectonic provinces as well as subtler geological features indicating our results are robust. Our results reveal for the first time the full extent of the mid-crustal Andean low-velocity zone that we tentatively interpret as the signature of a very large volume Neogene batholith. This study demonstrates the efficacy of integrating seismic data from numerous regional broad-band seismic networks to approximate the high-resolution coverage previously only available though larger networks such as the EarthScope USArray Transportable Array in the United States.

Description: The International Geomagnetic Reference Field (IGRF) is a very frequently used model of the Earth’s main magnetic field by both the science community and the industry. This model is updated every five years. We present here the second generation of the GRIMM magnetic field model that was derived to contribute to the IGRF-11. The model has been developed from a newly, reprocessed CHAMP satellite data set covering nearly 10 years. It has a temporal and spatial resolution significantly improved compared to previous models.

Description: The ability of any satellite gravity mission concept to monitor mass transport processes in the Earth system is typically tested well ahead of its implementation by means of various simulation studies. Those studies often extend from the simulation of realistic orbits and instrumental data all the way down to the retrieval of global gravity field solution time-series. Basic requirement for all these simulations are realistic representations of the spatio-temporal mass variability in the different sub-systems of the Earth, as a source model for the orbit computations. For such simulations, a suitable source model is required to represent (i) high-frequency (i.e., subdaily to weekly) mass variability in the atmosphere and oceans, in order to realistically include the effects of temporal aliasing due to non-tidal high-frequency mass variability into the retrieved gravity fields. In parallel, (ii) low-frequency (i.e., monthly to interannual) variability needs to be modelled with realistic amplitudes, particularly at small spatial scales, in order to assess to what extent a new mission concept might provide further insight into physical processes currently not observable. The new source model documented here attempts to fulfil both requirements: Based on ECMWF’s recent atmospheric reanalysis ERA-Interim and corresponding simulations from numerical models of the other Earth system components, it offers spherical harmonic coefficients of the time-variable global gravity field due to mass variability in atmosphere, oceans, the terrestrial hydrosphere including the ice-sheets and glaciers, as well as the solid Earth. Simulated features range from sub-daily to multiyear periods with a spatial resolution of spherical harmonics degree and order 180 over a period of 12 years. In addition to the source model, a de-aliasing model for atmospheric and oceanic high-frequency variability with augmented systematic and random noise is required for a realistic simulation of the gravity field retrieval process, whose necessary error characteristics are discussed. The documentation of the updated ESA Earth System Model (updated ESM) for gravity mission simulation studies is organized as follows: The characteristics of the updated ESM along with some basic validation is presented in Volume 1. A detailed comparison to the original ESA ESM (Gruber et al., 2011) is provided in Volume 2, while Volume 3 contains the description of a strategy to derive realistic errors for the de-aliasing model of high-frequency mass variability in atmosphere and ocean.

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: This paper presents the application of two different active contour models for the segmentation of high-resolution TerraSAR-X data. Both methods facilitate the detection of land-water-boundaries in semi-automated procedures and can be used to delineate flood extent and to map open water surfaces in general. For the extraction of smooth water bodies amplitude thresholding approaches are quite common and often applied. For rough water bodies however the application of amplitude thresholding methods is not successful. This paper demonstrates the poten-tials and limitations of active contour models for mapping both smooth and rough water bodies in high-resolution SAR data. Examples of both different segmentation methods are presented.

Description: Due to the recent eruptive and highly disruptive volcanic events in 2010 in Iceland, scientific and societal interest is overwhelming in gaining as much information as possible about the volcanic structures and processes to enhance the understanding of the partially glacier‐covered Eyjafjallajökull and Katla volcanic systems. Numerous petrological, geochemical and geophysical investigations of these systems have already been published. However, to date no electrical or electromagnetic data have been acquired on these two volcanoes to attempt to image the resistivity structure beneath and around them, although electromagnetic methods are far more sensitive to fluid distribution (in this case partial melt) than any other geophysical method. In July 2011, a pilot study took place to collect broadband magnetotelluric (MT) data around the Eyjafjallajökull. The MT data are supplemented with transient electromagnetic (TEM) measurements. This data set is the first one collected at these volcanic systems and will complement the existing geophysical data. Very fresh data and preliminary results from the experiment will be shown.

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: Throughout the Phanerozoic, the primary inorganic marine calcium carbonate mineralogy oscillated between calcite and aragonite, reflecting changes in seawater chemistry. These variations in seawater composition also appear to have influenced the evolution of calcifying organisms. However, the processes controlling these mineralogical and biological changes are poorly constrained. Previous work has focused mainly on the Mg/Ca ratio in seawater as the primary driver. Here, we examine the role of dissolved SO4 in these processes by performing controlled laboratory precipitation experiments and geochemical modeling of evaporite formation. We show that an increase in dissolved SO4 decreases the Mg/Ca ratio at which calcite is destabilized and aragonite becomes the dominant CaCO3 polymorph. Our data suggest that the Mg/Ca and SO4 thresholds for formation of calcite seas are significantly lower than previous estimates and are mutually dependent. This shows that our understanding of Phanerozoic changes in seawater chemistry and the models relating this to primary CaCO3 mineralogy need reevaluation.

Description: The addition of clay to growing media aims at constant supply of potassium, phosphorus and micro nutrients, pH-buffering, improvement of rewettability, and cohesion of growing media. The identification of suitable clays and their classification is a prerequisite for product and cultivation safety. A range of different clays highly variable in their mineral parameters were selected for experiments on nutrient buffering (P) and Mn toxicity, their ability to improve the rewettability and binding capacity of growing media. From batch experiments and growth trials it was derived that a threshold value for the sum of exchangeable and easy reducible Mn in clays for growing media is not justified, as even very high Mn contents in clay were not phytotoxic. The P binding capacity of clays was strongly correlated with the oxalate extractable Fe and Al content. A newly developed capillary rise method (WOK) was used to characterize the rewettability of growing media. The speed of rewetting mainly depends on the fineness of the amended clay. Surface free energy (SFE) data of the growing media indicate that those with a good rewettability show also high values for SFE. Compared to the kaolinitic and illitic clay amendments, bentonites show no significant increase in the SFE. Surfaces of clay minerals exhibiting originally polar and hydrophilic surfaces, can render hydrophobic when coated with weakly or non-polar organic matter moieties. Dissolved organic carbon (DOC) sorption was found to be positively correlated with the specific surface area (SSA), cation exchange capacity (CEC) and amount of dithionite extractable Al and Fe. Clays containing expandable clay minerals with high CEC and SSA (e.g. smectites) and those rich in Al- and Fe-oxides seem to be less effective for improving rewettability, whereas an addition of non-expandable clays with lower SSA, CEC (e. g. kaolinitic and illitic clays) and amorphous oxide content appears more promising. New insights on the adsorption of DOM on clay mineral`s surfaces will be obtained in the third year by chemical analysis of the surface with X- ray photoelectron spectroscopy (XPS), surface topography analysis and contact angle measurements. A new method for the determination of the binding capacity of clays in blocking media is currently being tested and will be studied together with determinations on K-dynamics.

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: A unique and very interesting earthquake of magnitude Mw 7.2 occurred in the Van region of Turkey on October 23, 2011 that caused a heavy loss of human lives and properties. The earthquake occurred on a blind oblique thrust fault oriented towards the NE–SW direction and dipping towards NW as evidenced by focal mechanism solution and aftershock distribution. In this study, we analyzed the seismogenesis and earthquake triggering during this sequence with the help of estimated seismological parameters (b-value of frequency–magnitude relation, p-value of aftershocks temporal decay and D-value of fractal dimension), 2D mapping of b- and p-values, 3D mapping of b-value and coseismic Coulomb stress modeling. The estimated seismic b-value equal to 0.89 reveals that the mainshock occurred in a highly stressed region and sequence comprised larger magnitude aftershocks due to the presence of large size asperities within the rupture zone. The normal estimate of p-value (0.98) suggests a tectonic genesis of the aftershocks sequence. The estimated D-value equal to 1.80 reveals that rupture propagated in a two-dimensional plane filled up by fractures. The spatial 2D and 3D mapping of seismic b-value suggests that the Van earthquake originated in a highly heterogeneous fractured rock matrix with fluid intrusions into it at deeper depth beneath the mainshock hypocenter region. The estimated coseismic Coulomb stress using the variable slip model for depth range 0–30 km exhibits a ‘butterfly’ pattern and most of the aftershocks fall (90%) in the region of enhanced Coulomb stress. This suggests that most of the aftershock activities have been triggered by transfer of positive Coulomb stress due to coseismic slip of the mainshock. The results estimated in the present study have potential useful implications in future seismic hazard assessment and risk mitigation in Van and the surrounding regions.

Description: We investigate the source of the Mw 7.6 Padang earthquake by inverting three-component global positioning system (GPS) data and broad-band regional seismic-displacement waveforms. The earthquake involved oblique-reverse slip either on an E–W, south-dipping plane, or on a N–S, west-dipping plane. Finite-fault inversions indicate that the rupture primarily propagated downdip and southwest from the hypocentre, with a scalar seismic moment between 3.4 and 3.7 × 1020 N m (Mw 7.62–7.65). Analysis of the seismic and geodetic data do not allow for unique identification of the causative focal plane. Aftershock patterns strongly suggest the E–W plane was the causative focal plane, but aligned geological structures in the downgoing plate favour the N–S plane. In either case, this unusually large intraslab earthquake has moved closer to failure the deepest portion of the overlying megathrust, which last ruptured during a great earthquake more than two centuries ago and is late in its seismic cycle.

Description: We present an attenuation model for midperiod Rayleigh waves in Central Asia and surrounding regions. This model is defined by maps of attenuation coefficient across the region of study in the period band 14–24 s. The model is constructed to characterize the regional variations in attenuation of seismic waves in the crust, which are related to the tectonic history of the studied territory, to calibrate the regional surface-wave magnitude scale, and to extend the teleseismic 'surface-wave magnitude – body wave magnitude' (Ms-mb) discriminant to regional distances. The construction of the model proceeds in three stages. The first stage in model construction is the measurement of Rayleigh wave spectral amplitudes. We collected and processed waveform data for 200 earthquakes occurring from 2003 to 2006 inside and around Eurasia, and used records of about 135 broadband permanent and temporary stations. This data set provided a sufficient number of spectral amplitude measurements between 14 and 24 s periods for the construction of two–dimensional tomographic maps of attenuation coefficients. At the second stage of the work, the integral of attenuation coefficients along given paths is estimated using both inter-station measurements and single-station measurements corrected for source and receiver terms. The third stage includes the refining of source parameters, recalculation of attenuation coefficient integrals after this refinement, grooming of resulting coefficients, and multistage tomographic inversion of the data. Tomographic maps for the set of periods from 14 to 24 s, which exhibit clear correlation with geology and tectonics of the territory under study, were obtained. Validation of these maps using the inter-station measurements confirms their accuracy in predicting the observations.

Description: [...] In this thesis it is shown that the epoch-differenced ionospheric delay correction is sufficient for estimating the tropospheric delay, e.g., the Zenith Total Delay (ZTD), from SF GPS data. Based on this result, the Satellite-specific Epoch-differenced Ionospheric Delay model (SEID) was developed. In the SEID model the ionospheric corrections for SF data are generated from the observations of surrounding reference stations equipped with DF receivers. With the derived ionospheric corrections and the SF data, pseudo L2 data are generated, which can be processed using existing GPS processing software packages without any changes. [...]

Description: Three-component wide-angle seismic data acquired in southern Tibet during Project INDEPTH show strong P-to-S converted reflections from reflectors that are aligned at a depth of ∼15 kilometers beneath the northern Yadong-Gulu rift. These converted reflections are locally higher in amplitude than the corresponding P-wave reflections. Modeling of reflection mode conversion as a function of incidence angle indicates that this condition obtains for a reflector that is a solid over fluid interface; it is not typical of a solid-solid interface. The likely candidates for a fluid trapped within the crystalline crust of southern Tibet are granitic magma and water (brine).

Description: An ascorbic acid extraction at pH 7.5 has been examined to assess the influence of reaction conditions (pH, ascorbic acid concentration) on the dissolution of Fe from synthetic 2-line ferrihydrite and from other Fe-bearing minerals. The method was highly selective for Fe in ferrihydrite with only small amounts of Fe extracted from other (oxyhydr)oxides or clays. The labile Fe extracted from the synthetic 2-line ferrihydrite stored as a slurry decreased with time, and high resolution microscopy showed that the older materials formed networked aggregates that slow down the dissolution. The apparent rate constant for the dissolution of fresh 2-line ferrihydrite (similar to 10(-3) s(-1)) was an order of magnitude larger than that for aged suspensions (similar to 10(-4) s(-1)). Fresh 2-line ferrihydrite that was filtered, freeze-dried, frozen and freeze-dried, and stored in the dry was even less readily dissolved (apparent rate constants similar to 10(-6) s(-1)). These ferrihydrites also contained networked aggregates and, additionally, appear to occur as granular aggregates (visible to the naked eye) during the early stages of dissolution. Storage of filtered, freeze-dried, and frozen and freeze-dried ferrihydrites, whether in water or air, produced similar dissolution behaviour because aggregation caused by de-watering decreases the labile Fe content and is not reversed by re-hydration. The determination of labile Fe in ferrihydrite requires that natural samples should be collected, stored and extracted wet. The most aged samples dissolved by parabolic dissolution kinetics indicate that the rates of dissolution were controlled by the diffusion of reactant into the internal porosity of aggregates. (C) 2010 Elsevier B.V. All rights reserved.

Description: The high elevation of the southern Puna plateau, the widespread melting of its crust, the gap in intermediate depth seismicity and the recent eruptions of ignimbrite complexes can be explained by delamination of the lithospheric mantle beneath it. To test this hypothesis, an array consisting of 73 broad band and short period seismic stations was deployed in the region for a period of 2 years starting in 2007. We inverted the data using the two plane wave approach and obtained 1-D and 3-D Rayleigh wave phase velocities. Our dispersion curve shows that at short periods (〈70 s) the phase velocities are slightly higher than those of the Tibetan plateau and lower than those of the Anatolian plateau. At periods of 100–140 s we observe a low velocity zone that might be remnant hot asthenosphere below a flat slab (7–10 Ma). We estimate the average continental lithosphere thickness for the region to be between 100 and 130 km. Our three-dimensional Rayleigh wave phase velocities show a high velocity anomaly at low frequencies (0.007, 0.008, and 0.009 Hz) slightly to the north of Cerro Galan. This would be consistent with the hypothesis of delamination in which a piece of lithosphere has detached and caused upwelling of hot asthenosphere, which in turn caused widespread alkaline-collision related volcanism. This interpretation is also corroborated by our shear wave velocity model, where a high velocity anomaly beneath the northern edge of Cerro Galan at 130 km depth is interpreted as the delaminated block on top of the subducting Nazca slab.

Description: The PISCO’94 (Proyecto de Investigatión Sismólogica de la Cordillera Occidental, 1994) seismological network of 31 digital broad band and short-period three-component seismometers was deployed in northern Chile between the Coastal Cordillera and the Western Cordillera. More than 5300 local seismic events were observed in a 100 day period. A subset of high-quality P and S arrival time data was used to invert simultaneously for hypocenters and velocity structure. Additional data from two other networks in the region could be included. The velocity models show a number of prominent anomalies, outlining an extremely thickened crust (about 70 km) beneath the forearc region, an anomalous crustal structure beneath the recent magmatic arc (Western Cordillera) characterized by very low velocities, and a high-velocity slab. A region of an increased Vp /Vs ratio has been found directly above the Wadati-Benioff zone, which might be caused by hydration processes. A zone of lower than average velocities and a high Vp /Vs ratio might correspond to the asthenospheric wedge. The upper edge of the Wadati-Benioff zone is sharply defined by intermediate depth hypocenters, while evidence for a double seismic zone can hardly be seen. Crustal events between the Precordillera and the Western Cordillera have been observed for the first time and are mainly located in the vicinity of the Salar de Atacama down to depths of about 40 km.

Description: The seismicity data file used for this study is represented by the earthquake catalogue CENEC for Europe north of 44°N (Grünthal et al. 2009a). This paper describes in detail how this homogeneous data file in terms of moment magnitudes Mw (with Mw greater than 3.5) has been derived. The degree of harmonization achieved in CENEC is quantitatively analysed in Grünthal et al. (2009b).

Description: A mobile volcano fast response system (VFRS) that can be used for volcano monitoring in case of volcanic unrest anywhere in the world is currently under development in Germany. The main goal of the project called Exupéry is to provide the communication technology for stations in the field and an expert system that collects all data from various sources, assembles them in a database, and allows users to assess the data through one common web GIS interface. The system also includes an integrated automatic alert level including the alert level estimation in order to characterize the activity state of the volcano. The web GIS interface serves as a decision support system to assist scientists and local authorities in deciding how to react in the case of volcanic unrest.

Description: The Laramide orogeny (80-50 Ma) resulted in thick-skinned deformation of the western United States, more than 700 km inboard of the plate boundary where the Farallon Plate was subducting below North America. Most studies conclude that this event was the result of low-angle or flat subduction of the Farallon plate, whereby horizontal compressive stress from the shallow slab produced inboard crustal compression and shortening. However, it is still not clear what factors caused the Farallon plate to shallow prior to Laramide time or how stress was transferred from the flat slab to the continental interior. Three hypotheses have been proposed for triggering flat subduction: (1) an increase in the westward velocity of the North American plate; (2) subduction of an buoyant oceanic plateau with abnormally thick crust and possibly a low density harzburgite mantle lithosphere layer; and (3) slab suction produced by subduction-induced mantle wedge flow and enhanced by the presence of thick Colorado Plateau lithosphere in the backarc. In this study, we use numerical models to study the development of low-angle subduction below a continental plate with a structure similar to that of the western US. The two-dimensional, plain strain models use the SOPALE code, in which Arbitrary Lagrangian-Eulerian finite element techniques are used to compute the coupled thermal-mechanical evolution of the lithosphere-upper mantle system. We first assess what factors are needed to dynamically develop low-angle subduction. We find that the main control is the continental velocity, with enhanced slab shallowing as the continental velocity increases. In order to create a section of horizontal (i.e., flat) subduction, a further requirement is the presence of an oceanic plateau with a low-density harzburgite layer. The slab suction force seems to be less capable of creating a flat slab than the other two factors. This may be due to the low viscosity in the mantle wedge in our models (10^19~10^20Pa s), which cannot produce a sufficiently large hydrodynamic force to flatten the oceanic plate. Future work will examine variations in the strength of both the continental plate and the interface between the continent and low-angle oceanic plate, in order to explore the relationship between flat subduction and Rocky Mountain foreland deformation during Laramide time.

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: Frequency-domain helicopter-borne electromagnetic (HEM) data sets are commonly interpreted using one-dimensional (1-D) modelling and inversion tools. In many cases this approach is valid (e.g., horizontally layered targets and groundwater applications) but there are areas of higher dimension that are not recovered correctly applying 1-D methods. In recent years multi-dimensional forward modelling codes became available. As there is no analytic solution to verify (or falsify) the obtained numerical solutions, comparison with 1-D values as well as amongst various two-dimensional (2-D) and three-dimensional (3-D) codes is essential. Values obtained over the background structure and at the centre of a large (a few hundred metres edge length) structure should match 1-D values. At the edges of the structure – i.e., at higher dimensional parts – only a comparison of different codes is possible. The more codes – especially if based on different methods and/or written by different programmers – agree the more reliable is the obtained synthetic data set.

Description: The seismic signature of the Moho from which geologic and tectonic evolution hypotheses are derived is to a large degree a result of the seismic methodology which has been used to obtain the image. Seismic data of different types, passive source (earthquake) broad-band recordings, and controlled source seismic refraction, densely recorded wide-angle deep seismic reflection, and normal incidence reflection (using VibroseisTM, explosives, or airguns), have contributed to the description of the Moho as a relatively complex transition zone. Of critical importance for the quality and resolution of the seismic image are the acquisition parameters, used in the imaging experiments. A variety of signatures have been obtained for the Moho at different scales generally dependent upon bandwidth of the seismic source. This variety prevents the development of a single universally applicable interpretation. In this way source frequency content, and source and sensor spacing determine the vertical and lateral resolution of the images, respectively. In most cases the different seismic probes provide complementary data that gives a fuller picture of the physical structure of the Moho, and its relationship to a petrologic crust–mantle transition. In regional seismic studies carried out using passive source recordings the Moho is a relatively well defined structure with marked lateral continuity. The characteristics of this boundary change depending on the geology and tectonic evolution of the targeted area. Refraction and wide-angle studies suggest the Moho to be often a relatively sharp velocity contrast, whereas the Moho in coincident high quality seismic reflection images is often seen as the abrupt downward decrease in seismic reflectivity. The origin of the Moho and its relation to the crust–mantle boundary is probably better constrained by careful analysis of its internal details, which can be complex and geographically varied. Unlike the oceanic Moho which is formed in a relatively simple, well understood process, the continental Moho can be subject to an extensive variety of tectonic processes, making overarching conclusions about the continental Moho difficult. Speaking very broadly: 1) In orogenic belts still undergoing compression and active continental volcanic arcs, the Moho evolves with the mountain belt, 2) In collapsed Phanerozoic orogenic belts the Moho under the collapse structure was formed during the collapse, often by a combination of processes. 3) In regions having experienced widespread basaltic volcanism, the Moho can result from underplated basalt and basaltic residuum. In Precambrian terranes the Moho may be as ancient as the formation of the crust, in others Precambrian tectonic and magmatic processes have reset it. We note that seismic reflection data in Phanerosoic orogens as well as from Precambrian cratonic terranes often show thrust type structures extending as deep as the Moho, and suggest that even where crust and mantle xenoliths provide similar age of formation dates, the crust may be semi-allochothonous

Description: The integrated plate boundary in Chile (IPOC) combines 15 broadband stations with strong-motion sensors, GPS, strain sensors and magneto-telluric stations. The Chilean subduction zone setting provides a high background rate of seismicity (crustal, intermediate depth, and plate interface) in a region with exceptionally low ambient noise, particularly at higher frequencies. We have deployed seismic mini-arrays in the vicinity of IPOC stations PB02 and PB07, and installed a third array to the east of these stations near the village of Quillagua, such that all three arrays form a triangle. Each array has 10 elements and an aperture in the km range. The study area lies just to the north of the northern boundary of the rupture area of the Tocopilla earthquake of 2007 Mw=7.7) and just above or slightly to the east of the downdip limit of plate interface seismicity. Installing the mini-arrays in the area of the existing IPOC has the following advantages: • Independent knowledge of background structure and seismicity from existing and ongoing studies. • Should any transients or other unusual signals be found in the array data, we can look for anomalous signals in geodetic and MT recordings, which will help to narrow down possible underlying mechanisms.

Description: A local, broadband, seismic network of four observatory-quality stations (KTB-NET) was operated during the drilling of the KTB hole, within the framework of the interdisciplinary German Continental Deep Drilling Program (KTB). The aim was to investigate the seismic activity with regard to the tectonic stress field and to compare it with data from in situ measurements in the 9.1-km deep borehole sections (bottom temperature of 260°C and heat flow of 82–85 mW/m2). From October 1990 to November 1995, over 80 local microearthquakes with magnitudes from 0.2 to 2.8 ML were recorded: eight small events by the KTB-NET only and four earthquake swarms with 73 events by the KTB-NET and stations of the Vogtland/Western Bohemia networks. Six of the small events are located within or close to the KTB-NET. The swarm events occurred at the southwestern extension of the Ohre rift, in an area 20 km north of the drill site, which is revealed to be part of the Vogtland/Western Bohemia seismotectonic unit, characterized by swarm activity. The hypocenters are limited to the upper 13 km of the crust, with a distinct concentration between 10 and 12 km. All types of fault plane solutions are found, but at depths greater than 8 km, reversed faulting mechanisms predominate. P axes are very uniformly oriented in a NNW-SSE direction, corresponding to the well-known regional stress orientation in central Europe and in agreement with the special in situ stress measurements of the KTB program. The focal mechanism of a ML=1.2 event induced by a fluid injection experiment fits into the results obtained from the natural events. Possible indications for the brittle-ductile transition are discussed in view of the observed earthquake depth and focal mechanism distributions.

Description: Analysis of seismic anisotropy in the crust and the uppermost mantle gives lots of information about the ambient mantle flow, stress state, and the dynamic processes inside the Earth. Thus, seismic anisotropy and its main distinctive features beneath the southeastern Mediterranean region are studied through the analysis of teleseismic shear-wave splitting observed at six broadband seismic stations belonging to the GEOFON and the MedNet. Although the number of the recording stations is small; a total of 495 splitting parameters are obtained, which revealed significant variations in the observed fast polarization directions beneath the study area. The stations in northern Egypt and Cyprus show fast velocity directions oriented roughly N–S to NNE–SSW, coincident with many previous results. A slightly different splitting pattern comprising NE–SW fast polarization directions is observed in the stations located along the Dead Sea fault in the southeastern Mediterranean; which are consistent with the current strike-slip motion between Africa and Arabia. In addition, NW–SE fast polarization directions are recognized in the latter group. The observed delay times vary greatly but their averages lie between 0.35 and 1 s. Although large-scale mechanisms, such as the absolute plate motion of Africa and Arabia towards Eurasia and the differential motion between Arabia and Africa can be invoked to predominantly explain the origin of anisotropic features, we suggest that density-driven flow in the asthenosphere is a possible additional cause of the wide range of the splitting pattern observed beneath some stations.

Description: Geophysical experiments next year in Romania may provide insight into a common but short-lived seismic process that can be observed and understood at only one spot on Earth at present. About 150 stations will be set up in the Vrancea area in the southeast Carpathian Mountains to, in effect, record the terminal phase of the detachment of a subducting slab of oceanic lithosphere. This is a major regional tomographic study using a large number of broadband seismometers, which will operate for 6 months. Images will be used for hazard assessment as well as for a delineation of detachment history.Active subduction of oceanic lithosphere at convergent plate boundaries involves earthquakes, magmatism, metamorphism, and deformation—some of the most vivid manifestations of any plate tectonic process. The initiation and termination of subduction, however, remains relatively poorly understood. When convergence of lithospheric plates ceases and the suction force of the subducting plate becomes negligible, the subducting slab moves into an almost vertical position. If subduction occurs in an arcuate geometry, the slab is likely to be segmented.

Description: For almost 10 years the KTB superdeep drilling project has offered an excellent field laboratory for adapting seismic techniques to crystalline environments and for testing new ideas for interpreting seismic reflections in terms of lithological or textural properties of metamorphic rock units. The seismic investigations culminated in a three-dimensional (3-D) reflection survey on a 19 × 19 km area with the drill site at its center. Interpretation of these data resulted in a detailed, structural model of the German Continental Deep Drilling Program (KTB) location with dominant, steep faults in the upper crust. The 3-D reflection survey was part of a suite of seismic experiments, ranging from wide-angle reflection and refraction profiles to standard vertical seismic profiles (VSP) and more sophisticated surface-to-borehole observations. It was predicted that the drill bit would meet the most prominent, steeply dipping, crustal reflector at a depth of about 6500–7000 m, and indeed, the borehole penetrated a major fault zone in the depth interval between 6850 and 7300 m. This reflector offered the rare opportunity to relate logging results, reflective properties, and geology to observed and modeled data. Post-Variscan thrusting caused cataclastic deformation, with partial, strong alterations within a steeply dipping reverse fault zone. This process generated impedance contrasts within the fault zone on a lateral scale large enough to cause seismic reflections. This was confirmed by borehole measurements along the whole 9.1 km deep KTB profile. The strongest, reflected signals originated from fluid-filled fractures and cataclastic fracture zones rather than from lithological boundaries (i.e., first-order discontinuities between different rock types) or from texture- and/or foliation-induced anisotropy. During the interpretation of seismic data at KTB several lessons were learned: Conventional processing of two-dimensional (2-D) reflection data from a presite survey showed predominantly subhorizontal layering in the upper crust with reflectivity striking in the Variscan direction. Drilling, however, revealed that all rock units are steeply dipping. This confirms that surface common depth point (CDP) seismics strongly enhances subhorizontal reflectivity and may thus produce a very misleading crustal image. Although this was shown for synthetic examples earlier, the KTB provides the experimental proof of how crucial this insight can be.

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: Soil moisture is a key state variable that controls runoff formation, infiltration and partitioning of radiation into latent and sensible heat. However, the experimental characterisation of near surface soil moisture patterns and their controls on runoff formation remains a challenge. This subject was one aspect of the BMBF-funded OPAQUE project (operational discharge and flooding predictions in head catchments). As part of that project the focus of this dissertation is on: (1) testing the methodology and feasibility of the Spatial TDR technology in producing soil moisture profiles along TDR probes, including an inversion technique of the recorded signal in heterogeneous field soils, (2) the analysis of spatial variability and temporal dynamics of soil moisture at the field scale including field experiments and hydrological modelling, (3) the application of models of different complexity for understanding soil moisture dynamics and its importance for runoff generation as well as for improving the prediction of runoff volumes. To fulfil objective 1, several laboratory experiments were conducted to understand the influence of probe rod geometry and heterogeneities in the sampling volume under different wetness conditions. This includes a detailed analysis on how these error sources affect retrieval of soil moisture profiles in soils. Concerning objective 2 a sampling strategy of two TDR clusters installed in the head water of the Wilde Weißeritz catchment (Eastern Ore Mountains, Germany) was used to investigate how well “the catchment state” can be characterised by means of distributed soil moisture data observed at the field scale. A grassland site and a forested site both located on gentle slopes were instrumented with two Spatial TDR clusters that consist of up to 39 TDR probes. Process understanding was gained by modelling the interaction of evapotranspiration and soil moisture with the hydrological process model CATFLOW. A field scale irrigation experiment was carried out to investigate near subsurface processes at the hillslope scale. The interactions of soil moisture and runoff formation were analysed using discharge data from three nested catchments: the Becherbach with a size of 2 km², the Rehefeld catchment (17 km²) and the superordinate Ammelsdorf catchment (49 km²). Statistical analyses including observations of pre-event runoff, soil moisture and different rainfall characteristics were employed to predict stream flow volume. On the different scales a strong correlation between the average soil moisture and the runoff coefficients of rainfall-runoff events could be found, which almost explains equivalent variability as the pre-event runoff. Furthermore, there was a strong correlation between surface soil moisture and subsurface wetness with a hysteretic behaviour between runoff soil moisture. To fulfil objective 3 these findings were used in a generalised linear model (GLM) analysis which combines state variables describing the catchments antecedent wetness and variables describing the meteorological forcing in order to predict event runoff coefficients. GLM results were compared to simulations with the catchment model WaSiM ETH. Hereby were the model results of the GLMs always better than the simulations with WaSiM ETH. The GLM analysis indicated that the proposed sampling strategy of clustering TDR probes in typical functional units is a promising technique to explore soil moisture controls on runoff generation and can be an important link between the scales. Long term monitoring of such sites could yield valuable information for flood warning and forecasting by identifying critical soil moisture conditions for the former and providing a better representation of the initial moisture conditions for the latter.zeige weniger

Description: Abflussentwicklung, Infiltration und die Umverteilung von Strahlung in latenten und sensiblen Wärmestrom werden maßgeblich durch die Bodenfeuchte der vadosen Zone gesteuert. Trotz allem, gibt s wenig Arbeiten die sich mit der experimentellen Charakterisierung der Bodenfeuchteverteilung und ihre Auswirkung auf die Abflussbildung beschäftigen. Der Fokus dieser Dissertation wurde darauf ausgerichtet: (1) die Methode des Spatial TDR und deren Anwendbarkeit einschließlich der Inversion des TDR Signals in heterogenen Böden zu prüfen, (2) die Analyse der räumlichen und zeitlichen Dynamik der Bodenfeuchte auf der Feldskala einschließlich Feldexperimenten und hydrologischer Modellierung, (3) der Aufbau verschiedener Modellanwendungen unterschiedlicher Komplexität um die Bodenfeuchtedynamiken und die Abflussentwicklung zu verstehen und die Vorhersage des Abflussvolumens zu verbessern. Um die Zielsetzung 1 zu erreichen, wurden verschiedene Laborversuche durchgeführt. Hierbei wurde der Einfluss der Sondenstabgeometrie und verschiedener Heterogenitäten im Messvolumen bei verschiedenen Feuchtegehalten untersucht. Dies beinhaltete eine detaillierte Analyse wie diese Fehlerquellen die Inversion des Bodenfeuchteprofils beeinflussen. Betreffend der Zielsetzung 2, wurden 2 TDR-Cluster in den Quellgebieten der Wilden Weißeritz installiert (Osterzgebirge) und untersucht, wie gut der Gebietszustand mit räumlich hochaufgelösten Bodenfeuchtedaten der Feldskala charakterisiert werden kann. Um die Interaktion zwischen Evapotranspiration und Bodenfeuchte zu untersuchen wurde das hydrologische Prozessmodell CATFLOW angewendet. Ein Beregnungsversuch wurde durchgeführt um die Zwischenabflussprozesse auf der Hangskala zu verstehen. Die Interaktion zwischen Bodenfeuchte und Abflussentwicklung wurde anhand von drei einander zugeordneten Einzugsgebieten analysiert. Statistische Analysen unter Berücksichtigung von Basisabfluss, Bodenvorfeuchte und verschiedenen Niederschlagscharakteristika wurden verwendet, um auf das Abflussvolumen zu schließen. Auf den verschiedenen Skalen konnte eine hohe Korrelation zwischen der mittleren Bodenfeuchte und dem Abflussbeiwert der Einzelereignisse festgestellt werden. Hierbei konnte die Bodenfeuchte genauso viel Variabilität erklären wie der Basisabfluss. Im Hinblick auf Zielsetzung 3 wurden “Generalised liner models” (GLM) genutzt. Dabei wurden Prädiktorvariablen die den Gebietszustand beschreiben und solche die die Meteorologische Randbedingungen beschreiben genutzt um den Abflussbeiwert zu schätzen. Die Ergebnisse der GLMs wurden mit Simulationsergebnissen des hydrologischen Gebietsmodells WaSiM ETH verglichen. Hierbei haben die GLMs eindeutig bessere Ergebnisse geliefert gegenüber den WaSiM Simulationen. Die GLM Analysen haben aufgezeigt, dass die verwendete Messstrategie mehrerer TDR-Cluster in typischen funktionalen Einheiten eine viel versprechende Methode ist, um den Einfluss der Bodenfeuchte auf die Abflussentwicklung zu verstehen und ein Bindeglied zwischen den Skalen darstellen zu können. Langzeitbeobachtungen solcher Standorte sind in der Lage wichtige Zusatzinformationen bei der Hochwasserwarnung und -vorhersage zu liefern durch die Identifizierung kritischer Gebietszustände für erstere und eine bessere Repräsentation der Vorfeuchte für letztere.zeige weniger

Description: The guidebook contains a thematic description an background information for each field day. The overall setting and development of the Saxo-Thuringian zone within the Variscan orogen is presented in Kroner and Romer (2010). This material is not included in this guide book (because of copyright reasons) and will be handed out to field trip participants as separate handout.

Description: In eastern Turkey, the ongoing convergence of the Arabian and African plates with Eurasia has resulted in the westward extrusion of the Anatolian Plate. To better understand the current state and the tectonic history of this region, we image crust and uppermost mantle structure with ambient noise tomography. Our study area extends from longitudes of 32° to 44°E. We use continuous data from two temporary seismic deployments, our 2006–2008 North Anatolian Fault Passive Seismic Experiment and the 1999–2001 Eastern Turkey Seismic Experiment, as well as from additional seismographs in the region. We compute daily cross-correlations of noise records between all station pairs and stack them over the entire time period for which they are available, as well as in seasonal subsets, to obtain interstation empirical Green's functions. After selecting interstation cross-correlations with high signal-to-noise ratios and measuring interstation phase velocities, we compute phase velocity maps at periods ranging from 8 to 40 s. At all periods, the phase velocity maps are similar for winter and summer subsets of the data, indicating that seasonal variations in noise sources do not bias our results. Across the study area, we invert the phase velocity estimates for shear velocity as a function of depth. The shear velocity model, which extends to 50 km depth, highlights tectonic features apparent at the surface: the Eastern Anatolian Plateau is a prominent low-velocity anomaly whereas the Kırşehir Massif has relatively fast velocities. There is a large velocity jump across the Inner Tauride Suture/Central Anataolian Fault Zone throughout the crust whereas the North Anatolian Fault does not have a consistent signature. In addition, in the southeastern part of our study area, we image a high velocity region below 20 km depth which may be the northern tip of the underthrusting Arabian Plate.

Description: The deep structure of the Upper Rhine Graben (URG), a continental rift in SW Germany and E France, is still poorly known. This deficit impedes a full understanding of the geodynamic evolution of this prominent rift. We study the lithosphere–asthenosphere structure using teleseismic waveforms obtained from the passive broad-band TIMO project across the central URG. The recovered, crust-corrected traveltime residuals relative to the iasp91 earth model are tiny (mostly less than 0.2–0.3 s). The average measured slowness (〈1 s deg−1) and backazimuth (〈5°) deviations are also very small and do not show any systematic wave front anomalies. These observed perturbation values are smaller than expected ones from synthetic 3-D ray tracing modelling with anomalies exceeding 2–3 per cent seismic velocity in the mantle. Thus there is no significant hint for any deep-seated anomaly such as a mantle cushion, etc. This result means that the rifting process did not leave behind a lower lithospheric signature, which could be clearly verified with high-resolution teleseismic experiments. The only significant traveltime perturbation at the central URG is located at its western side in the upper crust around a known geothermal anomaly. The upper crustal seismic anomaly with traveltime delays of 0.2–0.3 s cannot be explained with increased temperature alone. It is possibly related to a zone of highly altered granite. In the west of our network a traveltime anomaly (0.6–0.7 s delay) related with the Eifel plume is confirmed by the TIMO data set.

Description: The kinetics and mechanisms of nanoparticulate amorphous calcium carbonate (ACC) crystallization to calcite, via vaterite, were studied at a range of environmentally relevant temperatures (7.5-25 degrees C) using synchrotron-based in situ time-resolved Energy Dispersive X-ray Diffraction (ED-XRD) in conjunction with high-resolution electron microscopy, ex situ X-ray diffraction and infrared spectroscopy. The crystallization process occurs in two stages; firstly, the particles of ACC rapidly dehydrate and crystallize to form individual particles of vaterite; secondly, the vaterite transforms to calcite via a dissolution and reprecipitation mechanism with the reaction rate controlled by the surface area of calcite. The second stage of the reaction is approximately 10 times slower than the first. Activation energies of calcite nucleation and crystallization are 73 +/- 10 and 66 +/- 2 kJ mol(-1), respectively. A model to calculate the degree of calcite crystallization from ACC at environmentally relevant temperatures (7.5-40 degrees C) is also presented.

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: Fifteen earthquakes (Mw 4.1–6.4) occurring at ten major segments of the Sumatran Fault Zone (SFZ) were analyzed to identify their respective fault planes. The events were relocated in order to assess hypocenter uncertainty. Earthquake source parameters were determined from three-component local waveforms recorded by IRIS-DMC and GEOFON broadband lA networks. Epicentral distances of all stations were less than 10°. Moment tensor solutions of the events were calculated, along with simultaneous determination of centroid position. Joint analysis of hypocenter position, centroid position, and nodal planes produced clear outlines of the Sumatran fault planes. The preferable seismotectonic interpretation is that the events activated the SFZ at a depth of approximately 14–210 km, corresponding to the interplate Sumatran fault boundary. The identification of this seismic fault zone is significant to the investigation of seismic hazards in the region.