Double-planed deep seismic zone and upper-mantle structure in the Northeastern Japan Arc

Summary. The ScSp wave converted from the ScS wave at the boundary between the descending lithospheric slab and the mantle above it was clearly observed from a nearby deep earthquake with magnitude 7.7 at some stations of the seismic network of Tohoku University which covers the Tohoku District, the northeastern part of Honshu, Japan. By applying the three-dimensional seismic-ray tracing method, the location of this boundary was determined from the difference in arrival time between the ScS and ScSp waves. The result shows that the upper boundary of the descending slab lies exactly on the upper plane of the double-planed deep seismic zone found in the Northeastern Japan Arc.
There is an additional evidence that the boundary is located on the upper plane of the double-planed deep seismic zone. The hypocentre distribution of intermediate-depth earthquakes located by the small-scale seismic-array observation is extremely different from that obtained by the relatively large-scale seismic network. The discrepancy in the distribution of hypocentres of the same earthquake independently located is well explained by the inclined lithospheric slab model derived from the difference in arrival time between the ScS and ScSp waves.
The earthquakes with reverse faulting or with down-dip compressional stresses occur at the upper boundary of the descending slab. Within the descending slab, the earthquakes with down-dip extensional stresses also occur in a very narrow zone from 30 to 40 km below the dipping boundary in the depth range from 50 to about 200 km, and these shocks form the lower plane of the double-planed deep seismic zone.     

Lengths of intermediate and deep seismic zones and temperatures in downgoing slabs of lithosphere

Summary. If intermediate and deep earthquakes occur in the coldest portions of the downgoing slabs of lithosphere, then different lengths of seismic zones represent different temperatures in the slabs. As the slab descends through the aesthenosphere, it warms primarily by conduction of heat through its upper and lower surfaces. Isotherms are advected downwards to distances approximately proportional both to the rate of subduction and to the square of the thickness of the lithosphere. Consequently, the lengths of seismic zones should be approximately proportional to the product of the rates and the squares of these thicknesses. As these thicknesses are approximately proportional to the square root of the age of the lithosphere, the lengths ought therefore to be approximately proportional to the product of the convergence rates times the ages. Although there is considerable scatter, observed lengths are approximately proportional to such products, and are not simply related to the rate, the age or the thickness alone. The data crudely fit the relationship: length = rate × age/10. Using this relationship, we infer that the Philippine sea and Pacific plates move slowly, if at all, with respect to one another and that the Farallon plate may have been too young to be subducted to a great distance beneath western North America in the Palaeogene. Calculations of temperatures at the depths of the deepest events suggest that these cut-off temperatures increase from about 600 ± 100°C at 200 km to 830 ± 50°C at 650 km depth, but the cut-off potential temperature is approximately constant. Assuming that the strength is a thermally activated parameter, and using the activation energy for olivine, a crude estimate of activation volume is obtained from the dependence of the cut-off temperature on depth.     

A hypothesis for the seismogenesis of a double seismic zone

The seismogenesis of a double seismic zone, in particular the lower layer of a double seismic zone, has not been adequately explained in the literature. On the basis of seismic data and geothermal structures along three well-studied cross-sections in the Kuril-Kamchatka and Japan subduction zones, we investigate the temperature/pressure conditions associated with seismogenic structures of the double seismic zones. the corresponding T/P loci seem to suggest that earthquakes observed in the lower layer and in the lower part (below approximately 130 ± 20 km) of the top layer of a double seismic zone were caused by metastable phase transition-a mechanism similar to that responsible for deep-focus earthquakes only at lower temperature/pressure conditions. Under this hypothesis, the wedge-shaped configuration of a double seismic zone is interpreted to represent the loci of the kinetic boundary of the phase transition. According to theoretical/experimental studies and the constraints imposed by our observations, a likely candidate for such a phase transition is the metastable Al-rich enstatite decomposing into the assemblage of Al-poor enstatite plus garnet. Earthquakes in the upper part of the top layer were most probably due to conventional mechanisms such as dehydration of subducted materials and/or facies change from basalt to eclogite. That the top layer involves more than one seismogenic mechanism is also implied by the distinct behaviour of seismicity in the vicinity of 130 ± 20 km. Because the presence of deviatoric stress is critical to the reaction rate of a metastable phase transition, it is inferred that single seismic zones are also caused by the same mechanisms, except that the implicit layer of a supposed double seismic zone is missing, due to the insufficient amount of appropriate metastable minerals or to the lack of appropriate deviatoric stresses in the source region.     

State of stress in the Southern Tyrrhenian subduction zone from fault-plane solutions

In this paper we present revised locations and original focal mechanisms computed for intermediate and deep earthquakes that occurred within the Southern Tyrrhenian subduction zone between 1988 and 1994, in order to improve our knowledge of the state of stress for this compressional margin. In particular, we define the stress distribution within a large portion of the descending slab, between 40 and about 450 km depth. The seismicity distribution reveals a continuous 40–50 km thick slab that abruptly increases its dip from subhorizontal in the Ionian Sea to a constant 70° dip in the Tyrrhenian. We computed focal mechanisms for events with magnitudes ranging from 2.7 and 5.7, obtaining the distribution of P - and T -axes for many events for which centroid moment tensor (CMT) solutions are not available, thus enabling the sampling of a larger depth range compared to previous studies. We define three portions of the slab characterized by different distributions of P - and T -axes. A general down-dip compression is found between 165 and 370 km depth, whereas in the upper part of the slab (40–165 km depth) the fault-plane solutions are strongly heterogeneous. Below 370 km the P -axes of the few deep events located further to the north have a shallower dip and are not aligned with the 70° dipping slab, possibly suggesting that they belong to a separated piece of subducted lithosphere. There is a good correspondence between the depth range in which the P -axes plunge closer to the slab dip (∼ 70°) and the interval characterized by the highest seismic energy release (190–370 km).     

The dynamical origin of subduction zone topography

Summary. Subduction zones are expressed topographically by long linear oceanic trenches flanked by a low outer rise on the seaward side and an island arc on the landward side. This topographic structure is reflected in free air gravity anomalies, suggesting that much of the topography originates from dynamical forces applied at the base of the crust. We have successfully reproduced the general topographic features of subduction zones by supposing that the stresses generated by the bending of the viscous lower lithosphere as it subducts are transmitted through the thin elastic upper portion of the lithosphere. The trench is due to a zone of extensional flow (associated with low pressure) in the upper part of the viscous lithosphere.
The stresses in the subducting slab are computed using a finite element technique, assuming a Maxwell viscoelastic constitutive relation. Various dips (10 to 90°) were investigated, as well as depth dependent and non-Newtonian (power law, n = 3) viscosities. Observed subduction zone dimensions are well reproduced by these models. The effective viscosity required at mid-depth in the lithosphere is about 6 × 10²² P. This low value is probably due to the stress dependence of the effective viscosity. However, these models also show that the topography of the subduction zone depends primarily upon the geometry of the subducting slab (dip, radius of curvature of the bend) rather than upon its rheology. Shear stresses beneath the trench reach maxima of approximately 50 MPa. An interesting feature of some solutions is a dynamically supported bench or platform between the trench and island arc.     

Lithoprobe seismic reflection profiling across Vancouver Island: results from reprocessing

Summary. Multichannel seismic reflection sections recorded across Vancouver Island have revealed two extensive zones of deep seismic reflections that dip gently to the northeast, and a number of moderate northeasterly dipping reflections that can be traced to the surface where major faults are exposed. Based on an integrated interpretation of these data with information from gravity, heat flow, seismicity, seismic refraction, magnetotelluric and geological studies it is concluded that the lower zone of gently dipping reflections is due to underplated oceanic sediments and igneous rocks associated with the current subduction of the Juan de Fuca plate, and that the upper zone represents a similar sequence of accreted rocks associated with an earlier episode of subduction. The high density/high velocity material between the two reflection zones is either an underplated slab of oceanic lithosphere or an imbricated package of mafic rocks. Reprocessing of data from two of the seismic lines has produced a remarkable image of the terrane bounding Leech River fault, with its dip undulating from >60° near the surface to 20° at 3 km depth and ∼38° at 6 km depth.     

Small-scale mantle flows and induced lithospheric stress near island arcs

Summary. This paper explores the middle ground between complex thermally-coupled viscous flow models and simple corner flow models of island arc environments. The calculation retains the density-driven nature of convection and relaxes the geometrical constraints of corner flow, yet still provides semianalytical solutions for velocity and stress. A novel aspect of the procedure is its allowance for a coupled elastic lithosphere on top of a Newtonian viscous mantle. Initially, simple box-like density drivers illustrate how vertical and horizontal forces are transmitted through the mantle and how the lithosphere responds by trench formation. The flexural strength of the lithosphere spatially broadens the surface topography and gravity anomalies relative to the functional form of the vertical flow stresses applied to the plate base. I find that drivers in the form of inclined subducting slabs cannot induce self-driven parallel flow; however, the necessary flow can be provided by supplying a basal drag of 1–5 MPa to the mantle from the oceanic lithosphere. These basal drag forces create regional lithospheric stress and they should be quantifiable through seismic observations of the neutral surface. The existence of a shallow elevated phase transition is suggested in two slab models of 300 km length where a maximum excess density of 0.2 g cm⁻³ was needed to generate an acceptable mantle flow. A North New Hebrides subduction model which satisfies flow requirements and reproduces general features of topography and gravity contains a high shear stress zone (75 MPa) around the upper slab surface to a depth of 150 km and a deviatoric tensional stress in the back arc to a depth of 70 km. The lithospheric stress state of this model suggests that slab detachment is possible through whole plate fracture.     

The role of soil processes in determining mechanisms of slope failure and hillslope development in a humid-tropical forest eastern Puerto Rico

Translational failures, with associated downslope earthflow components and shallow slides, appear to be the primary mechanism of hillslope denudation in the humid tropical forests of the mountains of eastern Puerto Rico. In-situ weathering of quartz diorite and marine-deposited volcaniclastics produces residual soil (saprolite; up to 21 m deep) / weathered rock profiles. Discontinuous zones of contrasting density and permeability particularly in quartz-diorite slopes at 0.5 m, and between 3 and 7 m, create both pathways and impedances for water that can result in excess pore pressures and, ultimately, aid in determining the location of failure planes and magnitudes of slope failures. In combination with relict fractures which create planes of weakness within the saprolite, and the potential significance of tensile stresses in the upper zone of saprolite (hypothesized to be caused by subsurface soil creep), shear failure can then occur during or after periods of heavy rainfall.Results of in-situ shear-strength testing show negative y-intercepts on the derived Mohr-Coulomb failure envelopes (approximately 50% of all tests) that are interpreted as apparent tensile stresses. Observation of tension cracks 1–2 m deep support the test data. Subsurface soil creep can cause extension of the soil and the development of tensile stresses along upper-slope segments. Shear-strength data support this hypothesis for both geologic types. Apparent values of maximum and mean tensile stress are greatest along upper slopes (16.5 and 6.29 kPa). Previously documented maximum rates of downslope movement coincided with local minima of shear strength, and the shear-strength minimum for all tests was located near 0.5 m below land surface, the shallow zone of contrasting permeabilities. These results indicate that subsurface soil creep, a slow semi-continuous process, may exert a profound influence on rapid, shallow slope failures in saprolitic soils.Data indicate that cove slopes in quartz diorite tend to be the most unstable when saturation levels reach 75%. Deep failures (7 m deep) appear the most critical but not the most frequent because pore pressure build-up will occur more rapidly in the upper perched zone of translocated clays before reaching the lower zone between 3 and 7 m. Frequent shallow failures could reduce the probability of deeper failures by removing overburden and reducing shear stress at depth. Deep failures are more likely to result from storm events of great duration and intensity.Sixty-six ‘naturally occurring’ and more than 100 ‘road-related’ landslides were mapped. Forest elevations exceed 1000 m, but the majority of these failures were found between 600 and 800 m in elevation. This appears to be the area where there is sufficient concentration of subsurface water to result in excess pore pressures. The high percentage of slope failures in the 600–800-m range, relative to the percentage at higher elevations, suggests that differences in soil-water processes are responsible for the form of these mountain slopes. Steep linear segments are maintained at higher elevations. Slope angles are reduced in the 600–800-m range by frequent shallow slides, creating a largely concave surface. In combination, slope segments above 800 m, and those between 600 and 800 m, produce the characteristic form of the mountains of eastern Puerto Rico.     

Late Cenozoic tectono‐stratigraphic sequences of the Crotone Basin: insights on the geodynamic history of the Calabrian arc and Tyrrhenian Sea

The Crotone Basin was generated in the late Cenozoic as a forearc basin of the Ionian arc‐trench system. New data are gained through detailed field mapping, high‐resolution stratigraphic analysis of a key area and examination of offshore well data and seismic reflection profiles. Major unconformities divide the basin fill into major sequences, which reveal a three‐stage internal organization thought to reflect geodynamic events of the Calabrian arc and backarc area closely. The first stage is characterized by extensional block faulting and uplift followed by rapid drowning during high subsidence and transtension in the basin along a major NNW‐ to NW‐striking fault system. This stage is interpreted to reflect resumption of rollback after an episode of slab tearing triggered by transitory docking of continental lithosphere in the trench. The initial uplift is inferred to reflect decoupling and rebound after the transitory coupling phase. The second stage is characterized by increased subsidence and continued extension/transtension. This trend presumably reflects a decreasing rate of rollback resulting from a tendency towards viscous coupling after acceleration of slab downwelling. The third stage is characterized by short‐lived transpression along major shear zones and local inversion of former basins. This is inferred to reflect entrance into the trench of buoyant continental lithosphere, resulting in significant deceleration of slab rollback and consequently a break in, or slowing of, backarc extension, and predominance of the effects of compression related to Africa–Europe convergence. Overall, the above evolution resulted in the formation of a progressively narrower and rapidly retreating slab, inducing extreme rates of backarc extension, and may have played a critical role in determining the intermittent nature of the backarc rifting.     

The deep seismicity of the Tyrrhenian Sea

The study reappraises the deep seismicity of the Tyrrhenian Sea. Careful examination of the quality of reported hypocentres shows that the earthquakes define a zone dipping NW, about 200 km along strike, 50 km thick, and reaching a depth of about 500 km. The zone is slightly concave to the NW at a depth of 300 km, but, contrary to many previous reports, is not tightly concave, nor are there significant spatial gaps in the seismicity, which is effectively continuous with depth. Seismicity is, however, concentrated in the depth interval 250–300 km, where the dip of the seismic zone changes from 70° (above 250 km) to a more gentle dip of 45° at greater depths. Seven fault-plane solutions are available for the largest earthquakes in this depth interval, all of them consistent with a P -axis down the dip of the seismic zone, and all of them requiring movement on faults out of the plane of the subducting slab.
Two deep earthquakes near Naples lie well outside the main zone of activity; for one of which a fault-plane solution is available that has a P -axis not aligned with the dip of the seismic zone. The tightly concave slab-geometry favoured by other reports is supported mainly by the location of these events near Naples, which we think may represent deformation in a separate, probably shallower dipping, piece of subducted lithosphere.
The lack of shallow seismicity, and particularly of thrust faulting earthquakes, at the surface projection of the Benioff zone suggests that active subduction has ceased. Estimates of the convergence rate responsible for subduction in the last 10 Myr far exceed the present convergence rate of Africa and Eurasia, suggesting that the subduction was related instead to the stretching and thinning of the crust in the Tyrrhenian Sea.     

New insight into the crust and upper mantle structure under Alaska

To better understand the seismic structure of the subducting Pacific plate under Alaska, we determined the three-dimensional P-wave velocity structure to a depth of approximately 200 km beneath Alaska using 438,146 P-wave arrival times from 10,900 earthquakes. In this study an irregular grid parameterization was adopted to express the velocity structure under Alaska. The number of grid nodes increases from north to south in the study area so that the spacing between grid nodes is approximately the same in the longitude direction. Our results suggest that the subducting Pacific slab under Alaska can be divided into three different parts based on its geometry and velocity structure. The western part has features similar to those in other subduction zones. In the central part a thick low-velocity zone is imaged at the top of the subducting Pacific slab beneath north of the Kenai Peninsula, which is believed to be most likely the oceanic crust plus an overlying serpentinized zone and the coupled Yakutat terrane subducted with the Pacific slab. In the eastern part, significant high-velocity anomalies are visible to 60–90 km depth, suggesting that the Pacific slab has only subducted down to that depth.     

Implications of stress concentration on a strike-slip fault in an elastic plate subject to basal shear stress

Summary We consider a long strike-slip fault in a lithosphere modelled as an elastic slab. To the base of the slab a shear stress distribution is applied which simulates the viscous drag exerted by the asthenosphere. The resulant stress on the fault plane may directly fracture the lithosphere in its brittle upper portion; alternatively it may give rise at first to a stable aseismic sliding in the lower portion. In the latter case, stress concentration due to the deep aseismic slip is the relevant feature of the pre-seismic stress acting on the upper section of the lithosphere. The two cases are examined by use of dislocation theory and their observable effects compared. Different depths of the aseismic slip zone and the presence or absence of a uniform friction on the seismic fault are allowed for. If the model is applied to the San Andreas fault region, where a steady sliding condition actually seems to be present at shallow depth, it turns out that the slip amplitudes commonly associated with large earthquakes are consistent with average basal stress values which can be substantially lower than a few bars, a value often quoted as the steady state basal stress due to a velocity gradient in the upper asthenosphere.     

Lithospheric structure of the Arabian and Eurasian collision zone in eastern Turkey from S-wave receiver functions

Crustal and upper-mantle seismic discontinuities beneath eastern Turkey are imaged using teleseismic S -to- P converted phases. Three crustal phases are observed: the Moho with depth ranging between 30 and 55 km, indicating variable tectonic regimes within this continental collision zone; an upper-crustal discontinuity at approximately 10 km depth; and various crustal low-velocity zones, possibly associated with recent Quaternary volcanism. Imaging of the upper mantle is complicated by the 3-D geometry of the region, in particular due to the Bitlis–Zagros suture zone. However, several upper-mantle S -to- P converted phase are identified as being the signature of the lithosphere–asthenosphere boundary (LAB). The inferred LAB for the Eastern Anatolian Accretionary Complex indicates that eastern Turkey has an anomalously thin (between ∼60 and 80 km) lithosphere which is consistent with an oceanic slab detachment model. The observed LAB phases for the Arabian shield and Iranian plateau indicate that lithospheric thickness for these stable regions is on the order of 100 to 125 km thick, which is typical of continental margins.     

Estimate of the stress field in Kilauea's South Flank, Hawaii

We estimated stress and seismic strain tensors for the Kilauea volcano's south flank. the stress orientation inversion and the seismic strain calculation were performed using fault-plane solutions. the principal stress and seismic strain directions are approximately uniformly distributed in space and time during the interval covered by the data. However, the σ1, σ3 plane is approximately orthogonal to the 1, 3 plane. Therefore, a weak layer may exist beneath the south flank. σ1 has a plunge of 59° and an azimuth of 152°, with a 10° 95 per cent confidence range. We also developed a stress magnitude inversion to estimate magnitudes of boundary and interior stresses. In this inversion, the principal stress directions were taken as constraints in the seismic volume, and surface geodetic observations were used as data. the maximum magmatic pressure in Kilauea's rift zone is about 160 MPa. the direction of σ1 can be interpreted as the superposition of hydrostatic stress ( pgh ) and magmatic pressure. Without the constraint imposed by the direction of σ1, the estimated pressure is only 60MPa, the distribution of magmatic pressure may be similar to that of pgh . In contrast, the upper rift zone may be in tension. the shear stress in the rift zone is about one order of magnitude smaller than the maximum compressive stress, supporting the interpretation of magmatic flow as fluid in dikes or channels. the combination of stress orientation inversion, seismic strain calculation, and stress magnitude inversion performed in this study provides a means by which to estimate the stress state in seismic areas.     

Slab low-velocity layer in the eastern Aleutian subduction zone

Local earthquakes in the vicinity of the Alaskan Peninsula's Shumagin Islands often produce arrivals between the main P and S arrivals not predicted by standard traveltime tables. Based on traveltime and polarization, these anomalous arrivals appear to be from P -to- S conversions at the surface of the subducted Pacific Plate beneath the recording stations. The P -to- S conversion occurs at the top of a low-velocity layer which extends to at least 150 km depth and is 8 ˜ 2 per cent slower than the overlying mantle. The slab is ˜ 7 per cent faster than the mantle. The low-velocity layer contains the foci of the earthquakes in the upper plane of the double seismic zone and confines PS ray paths to lie within it. These observations indicate that layered structures persist to positions well past the surface location of the volcanic front. Reactions forming high-pressure minerals do not yield slab-like velocities until beyond the point that subduction zone magma genesis occurs. If the subducted oceanic crust forms the layer, it is subducted essentially intact.     

On the resolving power of tomographic images in the Aegean area

b
The imaging of upper mantle heterogeneity by seismic tomography is strongly limited by the uneven global distribution both of seismic recording stations and earthquake sources. This can result in a loss of resolution and significance in the final image, particularly when a sparse data set contains few ray paths which intersect at sufficiently high angles in the volume of interest. In order to investigate the theoretical resolving power of a previously published tomographic image of the Aegean area, synthetic tests of the inversion procedure using a ray-path matrix obtained in this previous study for local and teleseismic P -waves were carried out. The aim was to examine the extent to which the shape of a synthetic lithospheric slab penetrating to different depths is inherently distorted by the tomographic imaging procedure, and to compare the synthetic tomographic images with the results from the actual inversion. The distortion is found to take the form of an artificial stretching of the lithospheric slab. The maximum 'stretching factor', as indicated by the downdip displacement of the peak amplitude of the synthetic high-velocity anomaly, is found to be a factor of 2 or so, though the distortion is usually less than this. The peak amplitude of the tomographic image of a lithospheric slab is found from the inversion of traveltime data to be at depths at or below 400 km. This indicates that the high-velocity lithospheric slab in the Aegean penetrates deeper than the Benioff zone seismicity of about 200 km. However, no constraints of the maximum depth of penetration could be established with the data set used in the present work.     

Geometry and brittle deformation of the subducting Nazca Plate, Central Chile and Argentina

We use data from the Chile Argentina Geophysical Experiment (CHARGE) broad-band seismic deployment to refine past observations of the geometry and deformation within the subducting slab in the South American subduction zone between 30°S and 36°S. This region contains a zone of flat slab subduction where the subducting Nazca Plate flattens at a depth of ∼100 km and extends ∼300 km eastward before continuing its descent into the mantle. We use a grid-search multiple-event earthquake relocation technique to relocate 1098 events within the subducting slab and generate contours of the Wadati-Benioff zone. These contours reflect slab geometries from previous studies of intermediate-depth seismicity in this region with some small but important deviations. Our hypocentres indicate that the shallowest portion of the flat slab is associated with the inferred location of the subducting Juan Fernández Ridge at 31°S and that the slab deepens both to the south and the north of this region. We have also determined first motion focal mechanisms for ∼180 of the slab earthquakes. The subhorizontal T -axis solutions for these events are almost entirely consistent with a slab pull interpretation, especially when compared to our newly inferred slab geometry. Deviations of T -axes from the direction of slab dip may be explained with a gap within the subducting slab below 150 km in the vicinity of the transition from flat to normal subducting geometry around 33°S.     

Deep reflection seismics in the Precambrian of Sweden

Summary. A total of 161 km of deep seismic profiles have been shot in the region. One profile crosses the Protogine zone in SW Sweden. Over most of the profile short, weak reflectors are seen The only area with a concentration of reflectors is in the upper two seconds between the two tectonic zones. A nearly transparent area east of the Protogine zone is interperted as a deep granite intrustion. In the Siljan impact structure where four profiles were shot, the NE part of the structure is dominated by upper crustal high amplitude reflectors. Possible causes are discussed.     

Gravitational deformation and inherited structural control on slope morphology in the subduction zone of north‐central Chile (ca. 29–33°S)

Subduction zones provide direct insight into plate boundary deformation and by studying these areas we better understand tectonic processes and variability over time. We studied the structure of the offshore subduction zone system of the Pampean flat‐slab segment (ca. 29–33°S) of the Chilean margin using seismic and bathymetric constraints. Here, we related and analysed the structural styles of the offshore and onshore western fore‐arc. Overlying the acoustic top of the continental basement, two syn‐extensional seismic sequences were recognised and correlated with onshore geological units and the Valparaíso Forearc Basin seismic sequences: (SII) Pliocene‐Pleistocene and (SI) Miocene‐Pliocene (Late Cretaceous (?) to Miocene‐Pliocene) syn‐extensional sequences. These sequences are separated by an unconformity (i.e. Valparaíso Unconformity). Seismic reflection data reveal that the eastward dipping extensional system (EI) recognised at the upper slope can be extended to the middle slope and controlled the accumulation of the older seismic package (SI). The westward dipping extensional system (EII) is essentially restricted to the middle slope. Here, EII cuts the eastward dipping extensional system (EI), preferentially parallel to the inclination of the older sequences (SI), and controlled a series of middle slope basins which are filled by the Pliocene‐Pleistocene seismic sequence (SII). At the upper slope and in the western Coastal Cordillera, the SII sequence is controlled by eastward dipping faults (EII) which are the local reactivation of older extensional faults (EI). The tectonic boundary between the middle (eastern outermost forearc block) and upper continental slope (western coastal block) is a prominent system of trenchward dipping normal fault scarps (ca. 1 km offset) that resemble a major trenchward dipping extensional fault system. This prominent structural feature can be readily detected along the Chilean erosive margin as well as the two extensional sets (EI and EII). Evidence of slumping, thrusting, reactivated faults and mass transport deposits, were recognised in the slope domain and locally restricted to some eastern dipping faults. These features could be related to gravitational effects or slope deformation due to coseismic deformation. The regional inclination of the pre‐Pliocene sequences favoured the gravitational collapse of the outermost forearc block. We propose that the structural configuration of the study area is dominantly controlled by tectonic erosion as well as the uplift of the Coastal Cordillera, which is partially controlled by pre‐Pliocene architecture.     

Stresses due to phase changes in subduction zones and an empirical equation of state for the mantle

Summary. An empirical equation of state, assuming mineralogical equilibrium, is developed for the top 700 km of the mantle. Assuming a uniform viscosity, this equation of state is used to show that the stresses due to the changes in phase induced in a descending lithospheric plate in a subduction zone are an order of magnitude larger than those due to the negative buoyancy of the slab in the asthenosphere. The stresses predicted are well within the power law creep region for likely mantle materials and so the effective viscosity will vary within the slab. Consequently the stresses will be smaller than those of 7.0 × 10⁸N/m² obtained here using uniform viscosities. These stresses are relatively compressional near the sides of the slab and tensional in the centre.