ALBERT

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: 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: 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: 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: 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: It is common practice in the seismological community to use, especially for large earthquakes, the moment magnitude Mw as a unique magnitude parameter to evaluate the earthquake’s damage potential. However, as a static measure of earthquake size, Mwdoes not provide direct information about the released seismicwave energy and its high frequency content,which is the more interesting information both for engineering purposes and for a rapid assessment of the earthquake’s shaking potential. Therefore, we recommend to provide to disaster management organizations besides Mw also sufficiently accurate energy magnitude determinations as soon as possible after large earthquakes. We developed and extensively tested a rapid method for calculating the energy magnitude Me within about 10–15 min after an earthquake’s occurrence. The method is based on pre-calculated spectral amplitude decay functions obtained from numerical simulations of Green’s functions. After empirical validation, the procedure has been applied offline to a large data set of 767 shallow earthquakes that have been grouped according to their type of mechanism (strike-slip, normal faulting, thrust faulting, etc.). The suitability of the proposed approach is discussed by comparing our rapid Me estimates with Mw published by GCMT as well as with Mw and Me reported by the USGS. Mw is on average slightly larger than our Me for all types of mechanisms. No clear dependence on source mechanism is observed for our Me estimates. In contrast, Me from the USGS is generally larger than Mw for strike-slip earthquakes and generally smaller for the other source types. For ∼67 per cent of the event data set our Me differs ≤ ±0.3 magnitude units (m.u.) from the respective Me values published by the USGS. However, larger discrepancies (up to 0.8 m.u.) may occur for strike-slip events. A reason of that may be the overcorrection of the energy flux applied by the USGS for this type of earthquakes. We follow the original definition of magnitude scales, which does not apply a priori mechanism corrections to measured amplitudes, also since reliable fault-plane solutions are hardly available within 10–15 min after the earthquake origin time. Notable is that our uncorrected Me data show a better linear correlation and less scatter with respect to Mw than Me of the USGS. Finally, by analysing the recordings of representative recent pairs of strong and great earthquakes, we emphasize the importance of combining Mw and Me in the rapid characterization of the seismic source. They are related to different aspects of the source and may differ occasionally even more than 1 m.u. This highlights the usefulness and importance of providing these two magnitude estimates together for a better assessment of an earthquake’s shaking potential and/or tsunamigenic potential.

Description: The Sumatran margin suffered three great earthquakes in recent years (Aceh-Andaman 26 December 2004 Mw = 9.1, Nias 28 March 2005 Mw = 8.7, Bengkulu 12 September 2007 Mw = 8.5). Here we present local earthquake data from a dense, amphibious local seismic network covering a segment of the Sumatran margin that last ruptured in 1797. The occurrence of forearc islands along this part of the Sumatran margin allows the deployment of seismic land-stations above the shallow part of the thrust fault. In combination with ocean bottom seismometers this station geometry provides high quality hypocentre location for the updip end of the seismogenic zone in an area where geodetic data are also available. In this region, the Investigator Fracture Zone (IFZ), which consists of 4 sub-ridges, is subducted below the Sunda plate. This topography appears to influence seismicity at all depth intervals. A well-defined linear streak of seismicity extending from 80 to 200 km depth lies along the prolongation of closely spaced IFZ sub-ridges. More intermediate depth seismicity is located to the southeast of this string of seismicity and is related to subducted rough oceanic seafloor. The plate interface beneath Siberut Island which ruptured last in 1797 is characterised by almost complete absence of seismicity.

Description: We report on a rare example of aseismic response of a creeping fault to the earthquake cycle of a nearby megathrust. Interferometric synthetic aperture radar (InSAR) is used to detect and analyze shallow creep of two crustal faults at Mejilones Peninsula, Northern Chile, located in the hanging wall of the 2007 Mw 7.7 Tocopilla subduction earthquake. We generate two independent time series of surface deformation spanning ∼3.5 yr of late interseismic and ∼1.5 yr early postseismic deformation associated with this event. The analysis reveals creep on the Mejillones fault as well as on a previously unmapped fault to the west of the Mejillones fault. The InSAR deformation maps and distributed slip models obtained from the data reveal that fault creep reversed between the interseismic and postseismic periods. Given the regional stress field perturbations due to interseismic and coseismic deformation, we argue that the observed shallow creep and its slip reversal are directly linked to the megathrust seismic cycle. Moreover, from similar eastward dips but opposite slip directions of the two faults, we infer that fault strength must be very low and that the kinematics is controlled by crustal flexure associated with the seismic cycle on the underlying megathrust. ⺠New evidence of crustal fault creep in Northern Chile from InSAR. ⺠New evidence that sense of fault creep is governed by megathrust earthquake cycles. ⺠New evidence that stress induced by megathrust earthquake alters the sense of creep.

Description: The Himalaya and the Tibetan Plateau are uplifted by the ongoing northward underthrusting of the Indian continental lithosphere below Tibet resulting in lithospheric stacking. The layered structure of the Tibetan upper mantle is imaged by seismic methods, most detailed with the receiver function method. Tibet is considered as a place where the development of a future craton is currently under way. Here we study the upper mantle from Germany to northern Sweden with seismic S receiver functions and compare the structure below Scandinavia with that below Tibet. Below Proterozoic Scandinavia, we found two low-velocity zones on top of each other, separated by a high-velocity zone. The top of the upper low-velocity zone at about 100 km depth extends from Germany to Archaean northern Sweden. It agrees with the lithosphere-asthenosphere boundary (LAB) below Germany and Denmark. Below Sweden it is known as the 8°discontinuity, or as a mid-lithospheric discontinuity (MLD), similar to observations in North America. Seismic tomography places the LAB near 200 km in Scandinavia, which is close to the top of our deeper low-velocity zone. We also observed the bottom of the asthenosphere (the Lehmann discontinuity) deepening from 180 km in Germany to 260 km below Sweden. Remnants of old subduction in the upper about 100 km below Scandinavia and Finland are known from controlled source seismic experiments and local earthquake studies. Recent tomographic studies indicate delamination of the lithosphere below southern Scandinavia and northern Germany. We are suggesting that the large-scale layered structure in the Scandinavian upper mantle may be caused by processes similar to the ongoing lithospheric stacking in Tibet.