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.

Occurrence and type of dolomite in ODP Site 124-768 samples (Table 1)

Occurrence of deep-sea dolomites has been reported from numerous settings (for discussion see Lumsden, 1988). Different authors agree that dolomite formation in the pelagic realm is a relatively early diagenetic process (e.g., Jorgensen, 1983; Shimmield and Price, 1984; Kablanow et al., 1984; Kulm et al., 1984). Baker and Burns (1985) suggest that most of the pelagic dolomites formed within a few tens of meters below the seafloor within the zone of microbial sulfate reduction. According to Fuechtbauer and Richter (1988), dolomite can form in the deep-sea at a minimum temperature of 10°C. Other deep-sea dolomites are products of fluids derived from underlying evaporites or submarine weathering of basalts (Garrison, 1981). In some cases (Mullins et al., 1985; Dix and Mullins, 1988; Mullins et al., 1988), the existence of dolomite is linked to disconformities and its formation may have resulted from circulation of seawater through the sediment during prolonged exposure (Dix and Mullins, 1988, p. 287). At Site 768 (Fig. 1), lithified carbonate layers, some containing variable amounts of dolomite, occur below 201 mbsf (Miocene). These beds alternate with unconsolidated or semi-lithified marl layers interbedded in clays and siliciclastic turbidites. The irregular depth distribution of the limestone beds and the variation in preservation and recrystallization of the calcareous microfaunas suggest that lithification of carbonates at Site 768 not only reflects burial diagenesis as described by Garrison (1981) and others, but in part may be a selective, early diagenetic process. The different types and distribution of the dolomite additionally seem to support this assumption. The purpose of this report is to document the occurrence and textural nature of the dolomite at Site 768. Methods used were analyses of stained thin sections (Alizarin S and Ferrocyanide) and studies with the scanning electron microscope. No geochemical analyses (e.g., stable isotopes) were carried out; they will be the subject of further investigations.

Composition of Messinian sediments from the Tyrrhenian Sea

The data base for this study is represented by essentially nonevaporitic Messinian sediments recovered at ODP Sites 654, 653, 652, and 656 along the eastern Sardinian margin, and of the overlying early Pliocene oozes. Grain-size distribution, carbonate content, and microscopic observation of the sand size fractions were investigated. Messinian paleoenvironments, documented in the western Tyrrhenian Sea (ODP Sites 654 and 653), provide additional evidence supporting the deep basin desiccation model. A sharp lithologic contrast between early Pliocene pelagic oozes and latest Messinian conformable gypsiferous silts supports this model. The "lago-mare" biofacies was only occasionally observed in the shallowest site and is limited to the topmost part of the Messinian. Sites 652 and 656, lying in the deeper part of the Tyrrhenian and located on the downthrown side of an important eastward dipping fault system known as "Faglia centrale" are characterized by terrigenous sedimentation, with partly recycled minor evaporites. Of special interest is Site 652, where the thickness of the (probable) Messinian is 530 m. Sedimentary characters indicate a permanently subaqueous but nonmarine environment, with turbidites accumulating in a rapidly subsiding basin. According to the model proposed, this basin was fed by continental waters during times of maximum evaporitic draw-down, with temporary marine incursions from the west or southwest when the water level was higher. A basement ridge separated the evaporating pond from this endoreic lake located on the opposite (eastern) margin of the Tyrrhenian Basin, which was then limited to its western part. Post-Messinian reactivation of the "Faglia centrale" is necessary to account for the inversion of the relief.

Oxygen and hydrogen isotope ratios of pore waters from ODP sites in the Mediterranean Sea

Pore fluids from two ODP sites at Eastern Mediterranean mud volcanoes have been analyzed for their Cl concentration and their delta18O and deltaD isotopic composition. The Cl data span a wide range of concentrations, from extremely depleted with respect to seawater (as low as 60 mM) at the crest of Milano dome (site 970) to strongly enriched (up to 5.4 M) at Napoli dome (site 971). Chloride enrichment is known to be due to dissolving Messinian evaporites, whereas the source of the low-Cl fluid is deduced from stable isotope data presented here. The isotopic composition of the endmember fluid is found to be +10? for delta18O and -32? for deltaD for low- as well as for high-Cl waters. From this signature it can be concluded that neither gas hydrates nor meteoric water play a significant role in the freshening of the pore water. Several other processes altering the delta18O/deltaD composition of pore waters are discussed and considered to be of only negligible influence. The process characterizing the isotopic composition of the fluid is found to be clay mineral dehydration (mainly smectite-illite transformation), corresponding to a depth range of 3.5-7 km and an elevated temperature of about 120-165°C. A quantitative estimate shows that this reaction is capable of producing the observed extreme Cl depletion.

Interstitial water chemistry at DSDP Leg 79 Holes

Interstitial water analyses of samples collected at Sites 544-547 of DSDP Leg 79 are presented. In Site 547 chloride concentrations increase to almost 80% of the halite saturation values. Gypsum occurrences in the sediments immediately overlying the halite deposit can be explained in terms of migration of Ca**2+ and SO2**2- from the underlying evaporites. At shallower depths sulfate concentrations decrease rapidly as a result of sulfate reduction processes. The same processes lead to the removal of calcium in the form of calcium carbonate. At Site 547, the chloride concentration depth profile suggests a maximum of dissolved chloride which may be the result of advective flow from nearby (abput 6 km) evaporite salt diapirs.

(Table 3) Sulfur and oxygen isotopic composition of cold water soluble sulfate, and strontium content and isotopic composition of upper Miocene evaporite samples recovered from DSDP Site 23-227

The origin of three Red Sea submarine brine pools was investigated by analysis of the S and O isotope ratios of dissolved sulfate and Sr isotope ratios of dissolved Sr in the brines. Sulfur and O isotope ratios of sulfate and Sr isotope ratios of evaporitic source rocks for the brines were measured for comparison. The S, O and Sr isotope ratios of evaporites recovered from DSDP site 227 are consistent with an upper Miocene evaporites age. The Valdivia Deep brine formed by karstic dissolution of Miocene evaporites by overlying seawater and shows no signs of hydrothermal input. The Suakin Deep brines are derived from, or have isotopically exchanged with Miocene or older evaporites. There has been only minor dilution of the brine by overlying seawater. Strontium isotope ratios of Suakin brine may indicate addition of a minor (15%) amount of volcanic Sr to the brine, but there is no evidence of high temperature brine-rock interaction. The sulfate in the Atlantis II brine was apparently derived from seawater. The O isotope ratio of sulfate in the present Atlantis II brine could reflect isotopic exchange between seawater sulfate and the brine at approximately 255°C. Approximately 30% of the Sr in the Atlantis II brine is derived from the underlying basalt, probably by hydrothermal leaching. Atlantis II brine is the only known example from the Red Sea which has a significant high-temperature hydrothermal history.