Influence of the zone of weakness on dip angle and shear heating of subducted slabs

We study the importance of the zones of weakness and the pattern of downgoing flow in steady-state models of subducting lithosphere, which interacts mechanically and thermally with the ambient mantle. The non-linear system of governing equations consists of (i) the momentum equation in stream function formulation and (ii) the steady-state heat transfer equation including conduction and advection of heat and dissipation. A finite element method has been applied to this system. We consider the viscosity to be a non-linear function of both the temperature and the stream function. In steady-state two-dimensional (2D) flow, the stream function isolines follow material trajectories. They are used to follow the top of the subducting slab, which because of its possible increase in water content, is assumed to have a lower viscosity. The zone of weakness has been thus obtained in the self-consistent fashion since the stream function as well as the temperature are the output from our modeling and no a priori assumptions about the shape of the bending lithosphere are taken into account. It was shown that several orders decrease of viscosity in the zone of weakness is required to obtain the dip angle of about 45°. If the decrease of viscosity is not sufficient enough, the subducted slab either sinks almost vertically or does not exhibit a plate-like behavior. We have also demonstrated that shear heating can unrealistically increase at the zone of weakness for fast subductions if decrease of viscosity is underestimated.     

Controls on modern alluvial fan processes in the central Alps,northern Italy

Alluvial fan development in Alpine areas is often affected by catastrophic sedimentary processes associated with extreme ?oods events, causing serious risks for people living on the fans. Hazard assessment in these areas depends on proper identi?cation of the dominant sedimentary processes on the fans. Data from a set of 209 alluvial fans from the central Alps of Italy are presented in this paper and analysed with the help of various statistical techniques (linear regression, principal components analysis, cluster analysis, discriminant analysis and logistic regression). First, we used modern sedimentary facies and historical records (?ood events since 15th century), to distinguish between the two dominant sedimentary processes on alluvial fans: debris ?ows and stream?ows. Then, in order to analyse the main controls on past and present fan processes, 36 morphological, geological and land‐use variables were analysed. As with observations for arid‐environment fans, catchment morphology is the most in?uential factor in the study area, whereas geology and land use are minor controls. The role of climatic change and landsliding within the catchments also seems to be very important and is discussed. Statistical techniques also help in differentiating groups of alluvial fans by sets of controlling factors, including stage and type of evolution. Finally, by using discriminant analysis and logistic regression, we classi?ed alluvial fans according to the dominant sedimentary process, with a success rate ranging between 75 and 92 per cent. Copyright © 2004 John Wiley & Sons, Ltd.     

The effect of dissolved barium on biogeochemical processes at cold seeps

A numerical model was applied to investigate and quantify the biogeochemical processes fueled by the expulsion of barium and methane-rich fluids in the sediments of a giant cold-seep area in the Derugin Basin (Sea of Okhotsk). Geochemical profiles of dissolved Ba²⁺, Sr²⁺, Ca²⁺, SO₄²⁻, HS⁻, DIC, I⁻ and of calcium carbonate (CaCO₃) were fitted numerically to constrain the transport processes and the kinetics of biogeochemical reactions. The model results indicate that the anaerobic oxidation of methane (AOM) is the major process proceeding at a depth-integrated rate of 4.9 μmol cm⁻² a⁻¹, followed by calcium carbonate and strontian barite precipitation/dissolution processes having a total depth-integrated rate of 2.1 μmol cm⁻² a⁻¹. At the low seepage rate prevailing at our study site (0.14 cm a⁻¹) all of the rising barium is consumed by precipitation of barite in the sedimentary column and no benthic barium flux is produced. Numerical experiments were run to investigate the response of this diagenetic environment to variations of hydrological and biogeochemical conditions. Our results show that relatively low rates of fluid flow (<∼5 cm a⁻¹) promote the dispersed precipitation of up to 26 wt% of barite and calcium carbonate throughout the uppermost few meters of the sedimentary column. Distinct and persistent events (several hundreds of years long) of more vigorous fluid flow (from 20-110 cm a⁻¹), instead, result in the formation of barite-carbonate crusts near the sediment surface. Competition between barium and methane for sulfate controls the mineralogy of these sediment precipitates such that at low dissolved methane/barium ratios (<4-11) barite precipitation dominates, while at higher methane/barium ratios sulfate availability is limited by AOM and calcium carbonate prevails. When seepage rates exceed 110 cm a⁻¹, barite precipitation occurs at the seafloor and is so rapid that barite chimneys form in the water column. In the Derugin Basin, spectacular barite constructions up to 20 m high, which cover an area of roughly 22 km² and contain in excess of 5 million tons of barite, are built through this process. In these conditions, our model calculates a flux of barium to the water column of at least 20 μmol cm⁻² a⁻¹. We estimate that a minimum of 0.44 × 10⁶ mol a⁻¹ are added to the bottom waters of the Derugin Basin by cold seep processes, likely affecting the barium cycle in the Sea of Okhotsk.     

Chemical and isotopic equilibrium between CO2 and CH4 in fumarolic gas discharges: Generation of CH4 in arc magmatic-hydrothermal systems

The chemical and isotopic composition of fumarolic gases emitted from Nisyros Volcano, Greece, and of a single gas sample from Vesuvio, Italy, was investigated in order to determine the origin of methane (CH₄) within two subduction-related magmatic-hydrothermal environments.Apparent temperatures derived from carbon isotope partitioning between CH₄ and CO₂ of around 340°C for Nisyros and 470°C for Vesuvio correlate well with aquifer temperatures as measured directly and/or inferred from compositional data using the H₂O-H₂-CO₂-CO-CH₄ geothermometer. Thermodynamic modeling reveals chemical equilibrium between CH₄, CO₂ and H₂O implying that carbon isotope partitioning between CO₂ and CH₄ in both systems is controlled by aquifer temperature.N₂/³He and CH₄/³He ratios of Nisyros fumarolic gases are unusually low for subduction zone gases and correspond to those of midoceanic ridge environments. Accordingly, CH₄ may have been primarily generated through the reduction of CO₂ by H₂ in the absence of any organic matter following a Fischer-Tropsch-type reaction. However, primary occurrence of minor amounts of thermogenic CH₄ and subsequent re-equilibration with co-existing CO₂ cannot be ruled out entirely. CO₂/³He ratios and δ¹³C_CO2 values imply that the evolved CO₂ either derives from a metasomatized mantle or is a mixture between two components, one outgassing from an unaltered mantle and the other released by thermal breakdown of marine carbonates. The latter may contain traces of organic matter possibly decomposing to CH₄ during thermometamorphism.     

Ar / Ar dating of high-pressure rocks from the Ligurian Alps: Evidence for a continuous subduction–exhumation cycle

We present ³⁹Ar–⁴⁰Ar dating of phengite, muscovite and paragonite from a set of mafic and metasedimentary rocks sampled from the high-pressure (HP) metaophiolites of the Voltri Group (Western Alps) and from clasts in the basal layer conglomerates from the Tertiary molasse which overlie the high-pressure basement. The white mica-bearing rocks display peak eclogitic and blueschist-facies parageneses, locally showing complex greenschist-facies replacement textures. The internal discordance of age spectra is proportional to the chemical complexity of the micas. High-Si phengites from eclogite clasts record a ³⁹Ar–⁴⁰Ar age of ca. 49 Ma for the eclogite stage and ca. 43 Ma for the blueschist retrogression; phengites from a blueschist basement sample yield an age of ca. 40 Ma; low-Si muscovite from a metasediment dates the formation of the greenschist paragenesis at ca. 33 Ma. Our data indicate that the analyzed samples reached high-pressure conditions at different times over a time-span of c.a. 10 Ma. Subduction was continuing during exhumation and blueschist retrograde re-equilibration of higher-pressure, eclogite-facies rocks. This process kept the isotherms depressed, allowing the older HP-rocks to escape thermal re-equilibration. Our results, added to literature data, fit a tectonic model of a subduction–exhumation cycle, with different tectonic slices subducted at different times from Early Eocene until the Eocene–Oligocene boundary.     

Origin and evolution of Pliocene–Pleistocene granites from the Larderello geothermal field (Tuscan Magmatic Province, Italy)

Extensive, mainly acidic peraluminous magmatism affected the Tuscan Archipelago and the Tuscan mainland since late Miocene, building up the Tuscan Magmatic Province (TMP) as the Northern Apennine fold belt was progressively thinned, heated and intruded by mafic magmas. Between 3.8 and 1.3 Ma an intrusive complex was built on Larderello area (Tuscan mainland) by emplacement of multiple intrusions of isotopically and geochemically distinct granite magmas. Geochemical and isotopic investigations were carried out on granites cored during drilling exploration activity on the Larderello geothermal field. With respect to the other TMP granites the Larderello intrusives can be classified as two-mica granites due to the ubiquitous presence of small to moderate amounts of F-rich magmatic muscovite. They closely resemble the almost pure crustal TMP acidic rocks and do not show any of the typical petrographic features commonly observed in the TMP hybrid granites (enclaves, patchy zoning of plagioclase, amphibole clots). On the basis of major and trace elements, as well as REE patterns, two groups of granites were proposed: LAR-1 granites (3.8–2.3 Ma) originated by biotite-muscovite breakdown, and LAR-2 granites (2.3–1.3 Ma) generated by muscovite breakdown. At least three main crustal sources (at 14–23 km depth), characterized by distinct ε_Nd(t) and ⁸⁷Sr/⁸⁶Sr values, were involved at different times, and the magmas produced were randomly emplaced at shallow levels (3–6 km depth) throughout the entire field. The partial melting of a biotite-muscovite-rich source with low ε_Nd(t) value (about −10.5) produced the oldest intrusions (about 3.8–2.5 Ma). Afterwards (2.5–2.3 Ma), new magmas were generated by another biotite-rich source having a distinctly higher ε_Nd(t) value (−7.9). Finally, a muscovite-rich source with high ε_Nd(t) (about −8.9) gave origin to the younger group of granites (2.3–1.0 Ma). The significant Sr isotope disequilibrium recorded by granites belonging to the same intrusion is interpreted, as due to the short residence time of magmas in the source region followed by their rapid transfer to the emplacement level. Partial melting was probably triggered by multiple, small-sized mafic intrusions, distributed over the last 3.8 Ma that allowed temporary overstepping of biotite- and muscovite-dehydration melting reactions into an already pre-heated crust. Dilution in time of the magmatic activity probably prevented melt mingling and homogenization at depth, as well as the formation of a single, homogeneous, hybrid pluton at the emplacement level. Moreover the high concentrations of fluxing elements (B, F, Li) estimated for the LAR granites modified melt properties by reducing solidus temperatures, decreasing viscosity and increasing H₂O solubility in granite melts. The consequences were a more efficient, fast, magma extraction and transfer from the source, and a prolonged time of crystallization at the emplacement level. These key factors explain the long-lived hydrothermal activity recorded in this area by both fossil (Plio-Quaternary ore deposits) and active (Larderello geothermal field) systems.     

Tetrahedral Plots of the Phase Relations for Basalts

The phase relations of quaternary systems are generally represented by projections onto ternary compositional planes. Such projections often obscure relationships that would only be evident in a three-dimensional tetrahedral plot. The tetrahedral plot requires that compositions of the minerals and melts be transformed into Cartesian coordinates. It is shown here how this transformation is carried out. The application is demonstrated by tetrahedral plots of experimental melt compositions of partially molten lherzolite. Furthermore, the plot can be used to evaluate whether or not a particular basaltic composition represents a primary melt. The methods are applicable to any four-component system.     

High-P (P = 1.5–1.8 GPa) blueschist from Elba: Implications for underthrusting and exhumation of continental units in the Northern Apennines

The Acquadolce Subunit on the Island of Elba, Italy, records blueschist facies metamorphism related to the Oligocene–early Miocene stages of continental collision in the Northern Apennines. The blueschist facies metamorphism is represented by glaucophane- and lawsonite-bearing metabasite associated with marble and calcschist. These rock types occur as lenses in a schistose complex representing foredeep deposits of early Oligocene age. Detailed petrological analyses on metabasic and metapelitic protoliths, involving mineral and bulk-rock chemistry coupled with P–T and P–T–X(Fe₂O₃) pseudosection modelling using PERPLE_X, show that the Acquadolce Subunit recorded nearly isothermal exhumation from peak pressure–temperature conditions of 1.5–1.8 GPa and 320–370°C. During exhumation, peak lawsonite- and possibly carpholite- or stilpnomelane-bearing assemblages were overprinted and partially obliterated by epidote-blueschist and, subsequently, albite-greenschist facies metamorphic assemblages. This study sheds new light on the tectonic evolution of Adria-derived metamorphic units in the Northern Apennines, by showing (a) the deep underthrusting of continental crust during continental collision and (b) rapid exhumation along ‘cold’ and nearly isothermal paths, compatible with syn-orogenic extrusion.