The Vrancea seismic zone and its analogue in the Banda Arc, eastern Indonesia

It is now widely, although not universally, accepted that the Carpathian orogen marks the site of an arc–continent collision that followed the subduction of a now vanished small ocean basin. Seismic tomography has defined a high-velocity anomaly in the upper mantle similar to those associated with subduction zones worldwide. There is, however, no recognisable Wadati–Benioff Zone (WBZ), and intermediate-depth seismicity is confined to a relatively small, roughly cylindrical and vertically elongated region beneath the extreme southeastern corner of the mountain chain. There is no consensus in the published studies as to the origin of this ‘Vrancea Zone’.

The Banda Sea region of eastern Indonesia has sometimes been cited as an analogue for the Pannonian/Transylvanian basin and the enclosing Carpathian orocline, but at first sight the patterns of seismicity appear very different. Intermediate depth seismic activity defines a subducted slab that dips north, south and west beneath the Banda Sea, a configuration explained as a consequence of the rapid expansion of the sea during roll-back subduction. If the similar scenario proposed for the Carpathians is correct, then it is the absence of a Carpathian WBZ that is actually anomalous. Closer examination of Banda Arc seismicity shows that it can be divided into two parts, these being a scoop-shaped WBZ and an adjacent ‘Damar Zone’ of much more intense intermediate-depth activity. At its eastern end the Damar Zone merges with the WBZ, but in the west there is evidence for separation from it. A plausible explanation of this pattern is that a lower layer of the downgoing slab is peeling away from the remainder.

The Banda/Australia collision is now almost complete and the activity in the WBZ proper can be expected to decrease. Damar Zone activity, on the other hand, may persist for a much longer period, migrating towards the foreland as the detaching layer separates from the remainder of the subducted lithosphere. In a few million years the seismicity of the Banda region could well resemble the present day seismicity of the Carpathian orogen.     

Slow subduction of old lithosphere in the lesser antilles

The Lesser Antilles subduction zone is an extreme case of the subduction of old (~ 90 m.y.) lithosphere at a slow (~ 2 cm/y) convergence rate. Focal mechanisms of the largest earthquakes in the area have been obtained using body and surface wave data. During the time period (1950–1978) studied the subduction seismicity appears to represent primarily intraplate rather than interplate deformation. All three large (magnitude seven) earthquakes were from intraplate normal faults; no large thrust faulting earthquakes and few small ones occurred. These observations suggest that the plate boundary is largely decoupled, that subduction is at least partially aseismic, and that the downgoing slab is in a state of extension.     

大洋平板俯冲的数值模拟再现:洋–陆汇聚速率影响

利用地球动力学数值模拟方法探讨了洋-陆汇聚时,大洋岩石圈的绝对俯冲速率和上覆大陆岩石圈的向洋绝对逆冲速率对俯冲模式的影响,尤其是上覆大陆的向洋绝对逆冲速率与平板俯冲之间的关系。模型结果显示,对于年龄为40 Ma的含正常洋壳厚度的大洋岩石圈,在初始俯冲角度为现今洋–陆俯冲平均倾角的极小值(19°)条件下,低速大洋俯冲(绝对俯冲速率≤3 cm/a)且上覆大陆岩石圈向洋绝对逆冲速率≥1 cm/a时,具备形成平板俯冲的条件。当中–高速大洋俯冲(绝对俯冲速度3 cm/a)时,在上覆大陆的绝对逆冲速率不小于俯冲速率时可以形成平板俯冲。当增加初始俯冲角度到平均倾角的极大值(36°)时,仅在低速大洋俯冲(绝对俯冲速率≤3 cm/a)且绝对逆冲速率达到10 cm/a时(自然界中基本不存在),才有可能出现平板俯冲,其他情况均表现为陡俯冲。我们的模拟结果表明:(1)较高的大洋岩石圈绝对俯冲速率更容易克服板间耦合作用力而有利于陡俯冲形成;(2)较高的上覆大陆绝对逆冲速率更有利于俯冲板片弯曲而趋向于平板俯冲形成;(3)上覆大陆朝向海沟的逆冲速率会在俯冲板片下方产生水平向陆的地幔流,绝对逆冲速率越大该地幔流越强烈,导致作用于板片下表面的水平剪切分量越大而有利于板片弯折和平板俯冲发生;(4)初始俯冲角度的增加对平板俯冲的形成起到强烈抑制作用。这些能被现今平板俯冲,如具相似洋–陆汇聚速率条件的南美洲西海岸平板俯冲实例所验证。     

Subduction zone stresses: Constraints from mechanics and from topographic and geoid anomalies

There are, in principle, direct relations between several important phenomena associated with subduction zones: the depth of oceanic trenches, the magnitude of the net force from trenches acting on subducting plates, the distribution and fault plane orientations of earthquakes, the magnitude of stresses on subduction faults, the bathymetry of back-arc regions, and the magnitudes of gravity and geoid anomalies. These phenomena are related through the stresses transmitted through surface and subducted lithosphere, and are associated with the mass anomaly of the subducted lithosphere. Quantitative estimates suggest that observed trench depths imply a trench pull force on subducting plates which is comparable to the ridge push force but much less than the excess weight of the subducted lithospheric slab. It is further suggested that either the mass anomaly of subducted lithosphere is much less than would be expected on the basis of conventional thermal and compositional models or that (a) a large resistance acts on the upper part of slabs due to high-stress corner flow, and (b) the mass anomaly of the slab is 70–90% compensated either by a broad 1 km-deep back-arc depression or a low density mantle wedge above the slab or both.     

Slab dehydration and basalt petrogenesis in subduction systems involving very young oceanic lithosphere

Compositions of post-Miocene basalts erupted in the Garibaldi and Central America volcanic arcs exhibit significant correlations with the age of the subducted plate. In general, SiO₂, Al₂O₃, CaO, V, and (Sr/P)_N decrease and FeO, MgO, TiO₂ and Na₂O increase as the age of the subducted plate decreases. Variations in CaO/Al₂O₃, SiO₂, (Sr/P)_N, and Ba are compatible with lesser slab input, and hence less hydrous melting conditions in the mantle wedge in segments of the arcs overlying the youngest oceanic lithosphere. This interpretation is supported by comparison with peridotite melting experiments, which suggest higher melt pressures and temperatures in the mantle wedge above very young oceanic lithosphere. These observations point to a model in which dehydration of the downgoing slab occurs at shallow depths in subduction systems involving oceanic lithosphere younger than about 20 Ma. Because young oceanic lithosphere is relatively warm, little post-subduction heating is required to produce metamorphic reactions that release slab volatiles. Geodynamic models indicate most volatile-liberating reactions will occur within the seismogenic zone in oceanic lithosphere younger than 20 Ma, thus limiting the volatile flux beneath the arc and encouraging drier, higher temperature and higher pressure melting conditions in the mantle wedge in comparison to typical arc systems. Liberation of volatiles in the downgoing plate is strongly dependant on the shear stress on the fault, but is predicted to occur within the seismogenic zone for shear stresses greater than 33 MPa. Similarly, early loss of volatiles is predicted over a wide range of convergence rates, plate dips, and convergence angles. These results are shown to be robust for realistic ranges of slab dip, convergence angle, and shear stress, suggesting that volatile-poor melt generation is a characteristic of modern and ancient arc systems that involve subduction of young oceanic lithosphere.     

The latest geodynamics in Asia: Synthesis of data on volcanic evolution,lithosphere motion,and mantle velocities in the Baikal-Mongolian region

From a synthesis of data on volcanic evolution,movement of the lithosphere,and mantle velocities in the Baikal-Mongolian region,we propose a comprehensive model for deep dynamics of Asia that assumes an important role of the Gobi,Baikal,and North Transbaikal transition-layer melting anomalies.This layer was distorted by lower-mantle fluxes at the beginning of the latest geodynamic stage(i.e.in the early late Cretaceous) due to avalanches of slab material that were stagnated beneath the closed fragments of the Solonker,Ural-Mongolian paleoceans and Mongol-Okhotsk Gulf of Paleo-Pacific.At the latest geodynamic stage,Asia was involved in east-southeast movement,and the Pacific plate moved in the opposite direction with subduction under Asia.The weakened upper mantle region of the Gobi melting anomaly provided a counterflow connected with rollback in the Japan Sea area.These dynamics resulted in the formation of the Honshu-Korea flexure of the Pacific slab.A similar weakened upper mantle region of the North Transbaikal melting anomaly was associated with the formation of the Hokkaido-Amur flexure of the Pacific slab,formed due to progressive pull-down of the slab material into the transition layer in the direction of the Pacific plate and Asia convergence.The early—middle Miocene structural reorganization of the mantle processes in Asia resulted in the development of upper mantle low-velocity domains associated with the development of rifts and orogens.We propose that extension at the Baikal Rift was caused by deviator flowing mantle material,initiated under the moving lithosphere in the Baikal melting anomaly.Contraction at the Hangay orogen was created by facilitation of the tectonic stress transfer from the Indo-Asian interaction zone due to the low-viscosity mantle in the Gobi melting anomaly.     

Deep slab subduction and dehydration and their geodynamic consequences: Evidence from seismology and mineral physics

We present new pieces of evidence from seismology and mineral physics for the existence of low-velocity zones in the deep part of the upper mantle wedge and the mantle transition zone that are caused by fluids from the deep subduction and deep dehydration of the Pacific and Philippine Sea slabs under western Pacific and East Asia. The Pacific slab is subducting beneath the Japan Islands and Japan Sea with intermediate-depth and deep earthquakes down to 600 km depth under the East Asia margin, and the slab becomes stagnant in the mantle transition zone under East China. The western edge of the stagnant Pacific slab is roughly coincident with the NE–SW Daxing'Anling-Taihangshan gravity lineament located west of Beijing, approximately 2000 km away from the Japan Trench. The upper mantle above the stagnant slab under East Asia forms a big mantle wedge (BMW). Corner flow in the BMW and deep slab dehydration may have caused asthenospheric upwelling, lithospheric thinning, continental rift systems, and intraplate volcanism in Northeast Asia. The Philippine Sea slab has subducted down to the mantle transition zone depth under Western Japan and Ryukyu back-arc, though the seismicity within the slab occurs only down to 200–300 km depths. Combining with the corner flow in the mantle wedge, deep dehydration of the subducting Pacific slab has affected the morphology of the subducting Philippine Sea slab and its seismicity under Southwest Japan. Slow anomalies are also found in the mantle under the subducting Pacific slab, which may represent small mantle plumes, or hot upwelling associated with the deep slab subduction. Slab dehydration may also take place after a continental plate subducts into the mantle.     

On the thermal regime of some tectonic units in a continental collision environment in Romania

An analysis of the evolution of the main tectonic units in Romania shows that the thermal regime of the lithosphere should be derived according, on the one hand, to the particular tectonic interactions the various tectonic units have been involved in and, on the other hand, to the investigated depth interval. A steady-state conduction model of the crustal temperature field based on the heat flow distribution and information on the structure is presented for the Romanian territory. It shows large lateral thermal variations between tectonic units, as a result of different geological and thermal histories. Within the frame of a complex modelling of the thermal evolution of the lithosphere in the East Carpathians, as proposed in this study, the zero-order thermal effects of the pre-Miocene oceanic subduction of the Eurasian plate are evaluated. The deep thermal structure of the subducted slab is derived and shown to be compatible with the velocity structure of the lithosphere and the intermediate-depth seismicity of the Vrancea area.     

The Deep Seismicity of the Tyrrhenian Sea*

The deep seismicity of the Tyrrhenian Sea is analysed using data from a new instrumental catalogue of the seismicity of the Italian area. We use algorithms for the determination of absolute and relative hypocentral locations and for the evaluation of the geometry and coherence of the state of stress within the subducting slab. The structure of the Benioff zone, although simpler than previously indicated, reveals anomalous traits both in the seismicity distribution and in the stress geometry, confirming that standard subduction models cannot be applied in the Tyrrhenian region. The velocity anomaly and the location of few isolated events indicate that the subducted slabs extend to the north along the Apenninic chain approximately to the latitude of the Irpinia region, in Central Italy (? 42°N).     

Deep slab subduction and dehydration and their geodynamic consequences: Evidence from seismology and mineral physics

We present new pieces of evidence from seismology and mineral physics for the existence of low-velocity zones in the deep part of the upper mantle wedge and the mantle transition zone that are caused by fluids from the deep subduction and deep dehydration of the Pacific and Philippine Sea slabs under western Pacific and East Asia. The Pacific slab is subducting beneath the Japan Islands and Japan Sea with intermediate-depth and deep earthquakes down to 600 km depth under the East Asia margin, and the slab becomes stagnant in the mantle transition zone under East China. The western edge of the stagnant Pacific slab is roughly coincident with the NE–SW Daxing'Anling-Taihangshan gravity lineament located west of Beijing, approximately 2000 km away from the Japan Trench. The upper mantle above the stagnant slab under East Asia forms a big mantle wedge (BMW). Corner flow in the BMW and deep slab dehydration may have caused asthenospheric upwelling, lithospheric thinning, continental rift systems, and intraplate volcanism in Northeast Asia. The Philippine Sea slab has subducted down to the mantle transition zone depth under Western Japan and Ryukyu back-arc, though the seismicity within the slab occurs only down to 200–300 km depths. Combining with the corner flow in the mantle wedge, deep dehydration of the subducting Pacific slab has affected the morphology of the subducting Philippine Sea slab and its seismicity under Southwest Japan. Slow anomalies are also found in the mantle under the subducting Pacific slab, which may represent small mantle plumes, or hot upwelling associated with the deep slab subduction. Slab dehydration may also take place after a continental plate subducts into the mantle.     

Model calculations of the thermal fields of subducting lithospheric slabs and partial melting

A new numerical model that simulates a downgoing slab is used to study the conditions required to produce melting on its upper surface. Models with dip angles of 26.6° and 45°, rates of subduction of 0.7 and 5.6 cm y⁻¹, varying heat sources and rising material from the top of the slab are included. The results indicate that melting will not be greatly affected by dip angle, though the rate of subduction and the amount of shearstrain heating are important. When melting occurs, material rising from the top of the slab may produce high heat flow values at the surface of the earth on the continental side of the ocean trench, if the process continues sufficiently long. The sinking slab produces a positive gravity anomaly on the continental side of subduction, which is reduced in amplitude when rising material is present.     

Why Mt Etna?

The Etna volcano is located in an apparently anomalous position on the hinge zone of the Apennines subduction and its Na-alkaline geochemistry does not favour a magma source from the deep slab as indicated for the Aeolian K-alkaline magmatism. The steeper dip of the regional foreland monocline at the front of the Apennines in the Ionian Sea than in Sicily, implies a larger rollback of the subduction hinge in the Ionian Sea. Moreover, the lengthening of the Apennines arc needs extension parallel to the arc. Therefore, the larger southeastward subduction rollback of the Ionian lithosphere with respect to the Hyblean plateau in Sicily, should kinematically produce right-lateral transtension and a sort of vertical 'slab window' which might explain (i) the Plio-Pleistocene alkaline magmatism of eastern Sicily (e.g. the Etna volcano) and (ii) the late Pliocene to present right lateral transtensional tectonics and seismicity of eastern Sicily. The area of transfer of different dip and rollback occurs along the inherited Mesozoic passive continental margin between Sicily and the oceanic Ionian Sea, i.e. the Malta escarpment.     

Tomographic images of the upper mantle below central Europe and the Mediterranean

Results from delay time tomography of the European-Mediterranean upper mantle are discussed and where possible interpreted in terms of geodynamic processes. Slab-like positive velocity anomalies of which the locations correlate well with deeper seismicity are found beneath Spain, the Tyrrhenian basin, and the Aegean. These structures are interpreted as images of subducted slabs. Large aseismic regions with positive velocity anomalies are found beneath the Western Mediterranean, Italy, the Alps, Dinarides, the Pannonian basin, northern Greece, and the Aegean. These anomalies can also be linked to subducted lithosphere. From the anomaly patterns it is deduced that subduction occurred below the Western Mediterranean and along both sides of the Adriatic micro-plate. Beneath the Dinarides and northern Greece the velocity structures suggest detachment of the slab from the surface.     

Structure and dynamics of the Somma-Vesuvius volcanic complex

Summary We review recently obtained results about the velocity structure of the Somma-Vesuvius (Southern Italy) volcanic complex and present an interpretation of structural features, both at local and regional scale, and of the local seismicity. The local structure of Somma-Vesuvius is reviewed, referring to three depth ranges; i.e. shallow (0–5 km), intermediate (5–15 km) and deep (from 15 km to the upper mantle). The shallow velocity structure is inferred by the joint inversion of shot and local earthquake arrival time data. The main feature pointed out by this inversion is a high-velocity anomaly at the crater axis extending down to a depth of about 5 km. This anomaly can be explained with the presence of residual magma crystallised in the shallow conduits, which accumulated during the last eruptive cycles. The local seismicity is strongly clustered around this anomaly, due to the focusing effect of the rigidity contrast. The space-time seismicity pattern at Somma-Vesuvius is the result of the superposition of background seismicity, mainly due to gravitational instability of the volcanic edifice and to small external stress perturbations, with intense episodic earthquake swarms possibly due to magmatic or hydrothermal activity into the shallow system. The velocity structure in the 10–15 km depth range is characterized by the presence of a low-velocity layer, which has been independently confirmed by multi-channel seismic reflection data and P-Sv conversions from teleseismic waveforms. The study of the deep structure was performed by regional tomography with teleseisms; it confirmed the presence of a low-velocity anomaly underneath the volcano, which appears to have roots at greater depths. The regional structure between the Thyrrenian and the Adriatic sea has been inferred by tomographic inversion of teleseismic arrival times. The main result from this study which is very important for geodynamic interpretations is the first evidence for a continuous subducting slab under the Apennines, in an area where previous models hypothesized a slab window. Received March 3, 2000 revised version accepted July 4, 2001     

Late Cretaceous-Paleogene deepwater basin of North Afghanistan and the Central Pamirs: Issue of Hindu Kush earthquakes

The within-Iranian backarc basins, including the largest Sebzawar Basin, opened in the Mid-Cretaceous. Spreading in this basin was completed by the end of the Cretaceous. The basin closed in the Eocene with the formation of subduction zones and volcanic-plutonic belts. Data on North Afghanistan and the Central Pamirs have allowed us to reconstruct the eastern continuation of the Sebzawar Basin up to the west of the Central Pamirs. No fragments of oceanic crust are retained in Afghanistan and the Pamirs, but by analogy with the Sebzawar Basin, thick Paleogene flysch sequences and volcanic-plutonic complexes indicate setting of the active margin and subduction. It is suggested that the belt of mantle seismicity that extends for 550 km to the south of the Central Pamirs is related to the plunging and deformation of the lithosphere once underlying the Cretaceous-Paleogene basin. The extremely vigorous seismicity of the Hindu Kush megasource at the western termination of the seismic belt is caused by a number of specific tectonic features that predetermined the early onset of plunging of the subducted sheet (slab). In the megasource, the slab sank to a depth of 300 km and became vertical; its active deformation has proceeded up to the present. In the eastern part of the seismic belt, the slab started to plunge much later and therefore has retained a gentle slope, so that the depth of the hypocenters is shallower (down to 200 km), and earthquakes are less strong.     

Amount of Asian lithospheric mantle subducted during the India/Asia collision

Body wave seismic tomography is a successful technique for mapping lithospheric material sinking into the mantle. Focusing on the India/Asia collision zone, we postulate the existence of several Asian continental slabs, based on seismic global tomography. We observe a lower mantle positive anomaly between 1100 and 900 km depths, that we interpret as the signature of a past subduction process of Asian lithosphere, based on the anomaly position relative to positive anomalies related to Indian continental slab. We propose that this anomaly provides evidence for south dipping subduction of North Tibet lithospheric mantle, occurring along 3000 km parallel to the Southern Asian margin, and beginning soon after the 45 Ma break-off that detached the Tethys oceanic slab from the Indian continent. We estimate the maximum length of the slab related to the anomaly to be 400 km. Adding 200 km of presently Asian subducting slab beneath Central Tibet, the amount of Asian lithospheric mantle absorbed by continental subduction during the collision is at most 600 km. Using global seismic tomography to resolve the geometry of Asian continent at the onset of collision, we estimate that the convergence absorbed by Asia during the indentation process is ~ 1300 km. We conclude that Asian continental subduction could accommodate at most 45% of the Asian convergence. The rest of the convergence could have been accommodated by a combination of extrusion and shallow subduction/underthrusting processes. Continental subduction is therefore a major lithospheric process involved in intraplate tectonics of a supercontinent like Eurasia.     

Two-dimensional simulations of surface deformation caused by slab detachment

Detachment of the deeper part of subducted lithosphere causes changes in a subduction zone system which may be observed on the Earth's surface. Constraints on the expected magnitudes of these surface effects can aid in the interpretation of geological observations near convergent plate margins where detachment is expected. In this study, we quantify surface deformation caused by detachment of subducted lithosphere. We determine the range of displacement magnitudes which can be associated with slab detachment using numerical models. The lithospheric plates in our models have an effective elastic thickness, which provides an upper bound for rapid processes, like slab detachment, to the surface deformation of lithosphere with a more realistic rheology. The surface topography which develops during subduction is compared with the topography shortly after detachment is imposed. Subduction with a non-migrating trench system followed by detachment leads to a maximum surface uplift of 2–6 km, while this may be higher for the case of roll-back preceding detachment. In the latter situation, the back-arc basin may experience a phase of compression after detachment. Within the context of our elastic model, the surface uplift resulting from slab detachment is sensitive to the depth of detachment, a change in friction on the subduction fault during detachment and viscous stresses generated by sinking of the detached part of the slab. Overall, surface uplift of these magnitudes is not diagnostic of slab detachment since variations during ongoing subduction may result in similar vertical surface displacements.     

西藏及其周围地区地壳、地幔地震层析成像——印度板块大规模俯冲于西藏高原之下的证据

文中介绍有关西藏—喜马拉雅碰撞带的一项地震层析成像研究。根据一个用天然地震数据产生的全球波速模型 ,印度板块有可能以近水平状俯冲于整个西藏高原之下至 16 5～ 2 6 0km深度。西藏岩石圈具有低波速地壳和高波速下岩石圈 (75～ 12 0km深 )。在 12 0～ 16 5km深度范围 ,西藏岩石圈与俯冲的印度板块之间有一层低速软流圈物质。高原中部从地表到 310km深处有一低速体 ,说明地幔物质有可能穿过俯冲板块的脆弱部位上隆。这些结果以及野外实测的地壳缩短值说明高原的抬升得助于印度板块的近水平俯冲。我们推论俯冲印度板块的升温上浮以及上覆软流层的存在是造成西藏高原高海拔抬升以及内部地表仍相对平坦的主要原因。2 0 0 1年 1月 2 6日在印度西部发生的毁灭性大地震有可能是俯冲应力在印度板块后缘薄弱处引发的岩石圈大断裂。     

From oceanic plateaus to allochthonous terranes: Numerical modelling

Large segments of the continental crust are known to have formed through the amalgamation of oceanic plateaus and continental fragments. However, mechanisms responsible for terrane accretion remain poorly understood. We have therefore analysed the interactions of oceanic plateaus with the leading edge of the continental margin using a thermomechanical–petrological model of an oceanic-continental subduction zone with spontaneously moving plates. This model includes partial melting of crustal and mantle lithologies and accounts for complex rheological behaviour including viscous creep and plastic yielding. Our results indicate that oceanic plateaus may either be lost by subduction or accreted onto continental margins. Complete subduction of oceanic plateaus is common in models with old (> 40 Ma) oceanic lithosphere whereas models with younger lithosphere often result in terrane accretion. Three distinct modes of terrane accretion were identified depending on the rheological structure of the lower crust and oceanic cooling age: frontal plateau accretion, basal plateau accretion and underplating plateaus.Complete plateau subduction is associated with a sharp uplift of the forearc region and the formation of a basin further landward, followed by topographic relaxation. All crustal material is lost by subduction and crustal growth is solely attributed to partial melting of the mantle.Frontal plateau accretion leads to crustal thickening and the formation of thrust and fold belts, since oceanic plateaus are docked onto the continental margin. Strong deformation leads to slab break off, which eventually terminates subduction, shortly after the collisional stage has been reached. Crustal parts that have been sheared off during detachment melt at depth and modify the composition of the overlying continental crust.Basal plateau accretion scrapes oceanic plateaus off the downgoing slab, enabling the outward migration of the subduction zone. New incoming oceanic crust underthrusts the fractured terrane and forms a new subduction zone behind the accreted terrane. Subsequently, hot asthenosphere rises into the newly formed subduction zone and allows for extensive partial melting of crustal rocks, located at the slab interface, and only minor parts of the former oceanic plateau remain unmodified.Oceanic plateaus may also underplate the continental crust after being subducted to mantle depth. (U)HP terranes are formed with peak metamorphic temperatures of 400–700 °C prior to slab break off and subsequent exhumation. Rapid and coherent exhumation through the mantle along the former subduction zone at rates comparable to plate tectonic velocities is followed by somewhat slower rates at crustal levels, accompanied by crustal flow, structural reworking and syndeformational partial melting. Exhumation of these large crustal volumes leads to a sharp surface uplift.     

The Gibraltar Arc seismogenic zone (part 1): Constraints on a shallow east dipping fault plane source for the 1755 Lisbon earthquake provided by seismic data, gravity and thermal modeling

The Great Lisbon earthquake of 1755 with an estimated magnitude of 8.5–9.0 is the most destructive earthquake in European history, yet the source region remains enigmatic. Recent geophysical data provide compelling evidence for an active east dipping subduction zone beneath the nearby Gibraltar Arc. Marine seismic data in the Gulf of Cadiz image active thrust faults in an accretionary wedge, above an east dipping decollement and an eastward dipping basement. Tomographic and other data support subduction and rollback of a narrow slab of oceanic lithosphere beneath the westward advancing Gibraltar block.Although, no instrumentally recorded seismicity has been documented for the subduction interface, we propose the hypothesis that this shallow east dipping fault plane is locked and capable of generating great earthquakes (like the Nankai or Cascadia seismogenic zones). We further propose this east dipping fault plane to be a candidate source for the Great Lisbon earthquake of 1755. In this paper we use all available geophysical data on the deep structure of the Gulf of Cadiz–Gibraltar region for the purpose of constraining the 3-D geometry of this potentially seismogenic fault plane. To this end, we use new depth processed seismic data, have interpreted all available published and unpublished time sections, examine the distribution of hypocenters and perform 2-D gravity modeling. Finally, a finite-element model of the forearc thermal structure is constructed to determine the temperature distribution along the fault interface and thus the thermally predicted updip and downdip limits of the seismogenic zone.