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.

40Ar-39Ar dating of rocks drilled at sites 60-458 and 60-459

40Ar-39Ar dating of a high-MgO bronzite andesite from near the top of basement drilled at Site 458 shows the characteristic symptoms of artificially disturbed samples - i.e., an inverse staircase-type age spectrum, approximate linearity on an isochron plot, and concordance between total fusion age and isochron age. From conclusions based on other artificially disturbed samples (Ozima et al., 1979), we suggest that the reference isochron age (33.6 Ma) approximates the age of the sample. A basalt from deeper in Hole 458 gives an isochron age of 19.1 ± 0.2 Ma, which is slightly younger than the plateau age of 21.4 ± 1.0 Ma. Both ages are, however, considerably younger than the age of fossils in the overlying sediments (30 - 34 Ma). The age discrepancy may be explained if the 40Ar-39Ar age represents the age of secondary minerals, which formed later. No useful age data were obtained from a basalt sample recovered from Hole 459B.