Subduction and vertical coastal motions in the eastern Mediterranean

Convergence in the eastern Mediterranean of oceanic Nubia with Anatolia and the Aegean is complex and poorly understood. Large volumes of sediment obscure the shallow structure of the subduction zone, and since much of the convergence is accommodated aseismically, there are limited earthquake data to constrain its kinematics. We present new source models for recent earthquakes, combining these with field observations, published GPS velocities and reflection-seismic data to investigate faulting in three areas: the Florence Rise, SW Turkey and the Pliny and Strabo Trenches. The depths and locations of earthquakes reveal the geometry of the subducting Nubian plate NE of the Florence Rise, a bathymetric high that is probably formed by deformation of sediment at the surface projection of the Anatolia-Nubia subduction interface. In SW Turkey, the presence of a strike-slip shear zone has often been inferred despite an absence of strike-slip earthquakes. We show that the GPS-derived strain-rate field is consistent with extension on the orthogonal systems of normal faults observed in the region and that strike-slip faulting is not required to explain observed GPS velocities. Further SW, the Pliny and Strabo Trenches are also often interpreted as strike-slip shear zones, but almost all nearby earthquakes have either reverse-faulting or normal-faulting focal mechanisms. Oblique convergence across the trenches may be accommodated either by a partitioned system of strike-slip and reverse faults or by oblique slip on the Aegean-Nubia subduction interface. The observed late-Quaternary vertical motions of coastlines close to the subduction zone are influenced by the interplay between: (1) thickening of the material overriding the subduction interface associated with convergence, which promotes coastal uplift; and (2) subsidence due to extension and associated crustal thinning. Long-wavelength gravity data suggest that some of the observed topographic contrasts in the eastern Mediterranean are supported by mantle convection. However, whether the convection is time dependent and whether its pattern moves relative to Nubia are uncertain, and its contribution to present-day rates of vertical coastal motions is therefore hard to constrain. The observed extension of the overriding material in the subduction system is probably partly related to buoyancy forces arising from topographic contrasts between the Aegean, Anatolia and the Mediterranean sea floor, but the reasons for regional variations are less clear.AH is supported by a Shell Exploration studentship. This study forms part of the NERC- and 1069 ESRC-funded project "Earthquakes Without Frontiers" under grant NEJ02001X/1, and was 37 partly 1070 funded by the NERC grant "Looking inside the Continents from Space"

Microseismicity of the Mid-Atlantic Ridge at 7°S-8°15′S and at the Logatchev Massif oceanic core complex at 14°40′N-14°50′N

Lithospheric formation at slow spreading rates is heterogeneous with multiple modalities, favoring symmetric spreading where magmatism dominates or core complex and inside corner high formation where tectonics dominate. We report microseismicity from three deployments of seismic networks at the Mid-Atlantic Ridge (MAR). Two networks surveyed the MAR near 7 degrees S in the vicinity of the Ascension transform fault. Three inside corner high settings were investigated. However, they remained seismically largely inactive and major seismic activity occurred along the center of the median valley. In contrast, at the Logatchev Massif core complex at 14 degrees 45N seismicity was sparse within the center of the median valley but concentrated along the eastern rift mountains just west of the serpentine hosted Logatchev hydrothermal vent field. To the north and south of the massif, however, seismic activity occurred along the ridge axis, emphasizing the asymmetry of seismicity at the Logatchev segment. Focal mechanisms indicated a large number of reverse faulting events occurring in the vicinity of the vent field at 3-5 km depth, which we interpret to reflect volume expansion accompanying serpentinization. At shallower depth of 2-4 km, some earthquakes in the vicinity of the vent field showed normal faulting behavior, suggesting that normal faults facilitates hydrothermal circulation feeding the vent field. Further, a second set of cross-cutting faults occurred, indicating that the surface location of the field is controlled by local fault systems

Joint Inversion of Active and Passive Seismic Data in Central Java

Seismic and volcanic activities in Central Java, Indonesia, the area of interest of this study, are directly or indirectly related to the subduction of the Indo-Australian plate. In the framework of the MERapi AMphibious EXperiments (MERAMEX), a network consisting of about 130 seismographic stations was installed onshore and offshore in Central Java and operated for more than 150 days. In addition, 3-D active seismic experiments were carried out offshore. In this paper, we present the results of processing combined active and passive seismic data, which contain traveltimes from 292 local earthquakes and additional airgun shots along three offshore profiles. The inversion was performed using the updated LOTOS-06 code that allows processing for active and passive source data. The joint inversion of the active and passive data set considerably improves the resolution of the upper crust, especially in the offshore area in comparison to only passive data. The inversion results are verified using a series of synthetic tests. The resulting images showan exceptionally strong low-velocity anomaly (−30 per cent) in the backarc crust northward of the active volcanoes. In the upper mantle beneath the volcanoes, we observe a low-velocity anomaly inclined towards the slab, which probably reflects the paths of fluids and partially melted materials in the mantle wedge. The crust in the forearc appears to be strongly heterogeneous. The onshore part consists of two high-velocity blocks separated by a narrow low-velocity anomaly, which can be interpreted as a weakened contact zone between two rigid crustal bodies. The recent Java Mw = 6.3 earthquake (2006/05/26-UTC) occurred at the lower edge of this zone. Its focal strike slip mechanism is consistent with the orientation of this contact

Distribution of Temperature and Strength in the Central Andean Lithosphere and Its Relationship to Seismicity and Active Deformation

We present three-dimensional (3D) models of the present-day steady-state conductive thermal field and strength distribution in the lithosphere beneath the Central Andes. Our primary objective was to investigate the influence that the structure of the Central Andean lithosphere has on its thermal and rheological state, and the relationship between the latter and the active deformation in the region. We used our previous data-driven and gravity-constrained 3D density model as starting point for the calculations. We first assigned lithology-derived thermal and rheological properties to the different divisions of the density model and defined temperature boundary conditions. We then calculated the 3D steady-state conductive thermal field and the maximum differential stresses for both brittle and ductile behaviors. We find that the thickness and composition of the crust are the main factors affecting the modeled thermal field, and consequently also the strength distribution. The orogen is characterized by a thick felsic crust with elevated temperatures and a low integrated strength, whereas the foreland and forearc are underlain by a more mafic and thinner crust with lower temperatures and a higher integrated strength. We find that most of the intraplate deformation coincides spatially with the steepest strength gradients and suggest that the high potential energy of the orogen together with the presence of rheological lateral heterogeneities produce high compressional stresses and strong strain localization along the margins of the orogen. We interpret earthquakes within the modeled ductile field to be related to the weakening effect of long-lived faults and/or the presence of seismic asperities.Fil: Ibarra, Federico. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Geociencias Básicas, Aplicadas y Ambientales de Buenos Aires. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Geociencias Básicas, Aplicadas y Ambientales de Buenos Aires; ArgentinaFil: Prezzi, Claudia Beatriz. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Instituto de Geociencias Básicas, Aplicadas y Ambientales de Buenos Aires. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Instituto de Geociencias Básicas, Aplicadas y Ambientales de Buenos Aires; ArgentinaFil: Bott, Judith. German Research Centre for Geosciences; AlemaniaFil: Scheck Wenderoth, Magdalena. German Research Centre for Geosciences; AlemaniaFil: Strecker, Manfred. Universitat Potsdam; Alemani

P and S velocity sturcture of the crust and the upper mantle beneath Central Java from local tomography inversion

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

Constraints from GPS measurements on the dynamics of deformation in Anatolia and the Aegean

We estimate the strength of the lithosphere in Anatolia and the Aegean, and the boundary forces acting upon it, using a dynamical model that treats the lithosphere as a thin fluid sheet deforming in response to variations in gravitational potential energy. This model has one free material parameter, the power law exponent, n, of the vertically averaged rheology of the lithosphere, and two parameters that specify the forces per unit length applied to its edges. Solutions to this model that best fit the velocities of 346 reliable GPS sites require an effective viscosity of the lithosphere of 1022 to 1021 Pa s at strain rates of 10 to 100 nanostrain per year. The best-fitting force at the Arabia-Anatolia boundary is consistent with the lithostatic pressure due to the high topography there, and the force at the Nubia-Aegean boundary is consistent with the contrast in lithostatic pressure across that boundary. No additional force, from “slab rollback” or basal tractions due to convection in the mantle, is required to explain the observations. These results are supported by scaling relations derived from approximate analytical solutions. The inverse relationship between the viscosity of the lithosphere and deviatoric stress produces strong slowly deforming regions in the Southern Aegean and Central Anatolia whose motions resemble those of microplates. The distribution of geodetic strain rates within the region, and the partitioning between normal and strike-slip faulting, are explained by the interplay between boundary conditions, internal variations in gravitational potential energy, and the power law rheology of the lithosphere