A numerical discrete element study on the influence of an embedded viscoelastic-plastic layer at different viscosity values

The major aim of this study was to examine the influence of an embedded viscoelastic-plastic layer at different viscosity values on accretionary wedges at subduction zones. To quantify the effects of the layer viscosity, we analysed the wedge geometry, accretion mode, thrust systems and mass transport pattern. Therefore, we developed a numerical 2D 'sandbox' model utilising the Discrete Element Method. Starting with a simple pure Mohr Coulomb sequence, we added an embedded viscoelastic-plastic layer within the brittle, undeformed 'sediment' package. This layer followed Burger's rheology, which simulates the creep behaviour of natural rocks, such as evaporites. This layer got thrusted and folded during the subduction process. The testing of different bulk viscosity values, from 1 × 10**13 to 1 × 10**14 (Pa s), revealed a certain range where an active detachment evolved within the viscoelastic-plastic layer that decoupled the over- and the underlying brittle strata. This mid-level detachment caused the evolution of a frontally accreted wedge above it and a long underthrusted and subsequently basally accreted sequence beneath it. Both sequences were characterised by specific mass transport patterns depending on the used viscosity value. With decreasing bulk viscosities, thrust systems above this weak mid-level detachment became increasingly symmetrical and the particle uplift was reduced, as would be expected for a salt controlled forearc in nature. Simultaneously, antiformal stacking was favoured over hinterland dipping in the lower brittle layer and overturning of the uplifted material increased. Hence, we validated that the viscosity of an embedded detachment strongly influences the whole wedge mechanics, both the respective lower slope and the upper slope duplex, shown by e.g. the mass transport pattern.

Physical properties, velocity and attenuation characteristics, and mineralogy of gabbros from ODP Hole 176-735B

Laboratory compressional wave (Vp) and shear wave (Vs) velocities were measured as a function of confining pressure for the gabbros from Hole 735B and compared to results from Leg 118. The upper 500 m of the hole has a Vp mean value of 6895 m/s measured at 200 MPa, and at 500 meters below seafloor (mbsf), Vp measurements show a mean value of 7036 m/s. Vs mean values in the same intervals are 3840 m/s and 3857 m/s, respectively. The mean Vp and Vs values obtained from log data in the upper 600 m are 6520 and 3518 m/s, respectively. These results show a general increase in velocity with depth and the velocity gradients estimate an upper mantle depth of 3.32 km. This value agrees with previous work based on dredged samples and inversion of rare element concentrations in basalts dredged from the conjugate site to the north of the Atlantis Bank. Laboratory measurements show Vp anisotropy ranging between 0.4% and 8.8%, with the majority of the samples having values less than 3.8%. Measurements of velocity anisotropy seem to be associated with zones of high crystal-plastic deformation with predominant preferred mineral orientations of plagioclase, amphiboles, and pyroxenes. These findings are consistent with results on gabbros from the Hess Deep area and suggest that plastic deformation may play an important role in the seismic properties of the lower oceanic crust. In contrast to ophiolite studies, many of the olivine gabbros show a small degree of anisotropy. Log derived Vs anisotropy shows an average of 5.8% for the upper 600 m of Hole 735B and tends to decrease with depth where the overburden pressure and the age of the crustal section suggests closure of cracks and infilling of fractures by alteration minerals. Overall the results indicate that the average shear wave splitting in Hole 735B might be influenced by preferred structural orientations and the average value of shear wave splitting may not be a maximum because structural dips are <90°. The maximum fast-wave orientation values could be influenced by structural features striking slightly oblique to this orientation or by near-field stress concentrations. However, flexural wave dispersion analyses have not been performed to confirm this hypothesis or to indicate to what extent the near-field stresses may be influencing shear wave propagation. Acoustic impedance contrasts calculated from laboratory and logging data were used to generate synthetic seismograms that aid in the interpretation of reflection profiles. Several prominent reflections produced by these calculations suggest that Fe-Ti oxides and shear zones may contribute to the reflective nature of the lower oceanic crust. Laboratory velocity attenuation (Q) measurements from below 500 m have a mean value of 35.1, which is consistent with previous vertical seismic profile (VSP) and laboratory measurements on the upper 500 m.