Late Glacial ice advances in southeast Tibet

The sensitivity of Tibetan glacial systems to North Atlantic climate forcing is a major issue in palaeoclimatology. In this study, we present surface exposure ages of erratic boulders from a valley system in the Hengduan Mountains, southeastern Tibet, showing evidence of an ice advance during Heinrich event 1. Cosmogenic nuclide analyses (¹⁰Be and ²¹Ne) revealed consistent exposure ages, indicating no major periods of burial or pre-exposure. Erosion-corrected (3 mm/ka) ¹⁰Be exposure ages range from 13.4 to 16.3 ka. This is in agreement with recalculated exposure ages from the same valley system by [Tschudi, S., Schäfer, J.M., Zhao, Z., Wu, X., Ivy-Ochs, S., Kubik, P.W., Schlüchter, C., 2003. Glacial advances in Tibet during the Younger Dryas? Evidence from cosmogenic ¹⁰Be, ²⁶Al, and ²¹Ne. Journal of Asian Earth Sciences 22, 301–306.]. Thus this indicates that local glaciers advanced in the investigated area as a response to Heinrich event 1 cooling and that periglacial surface adjustments during the Younger Dryas overprinted the glacial morphology, leading to deceptively young exposure ages of certain erratic boulders.     

Southern Ocean warming and increased ice shelf basal melting in the twenty-first and twenty-second centuries based on coupled ice-ocean finite-element modelling

We utilise a global finite-element sea ice–ocean model (FESOM), focused on the Antarctic marginal seas, to analyse projections of ice shelf basal melting in a warmer climate. Ice shelf–ocean interaction is described using a three-equation system with a diagnostic computation of temperature and salinity at the ice–ocean interface. A tetrahedral mesh with a minimumhorizontal resolution of 4 km and hybrid vertical coordinates is used. Ice shelf draft, cavity geometry, and global ocean bathymetry have been derived from the RTopo-1 data set. The model is forced with the atmospheric output from two climate models: (1) the Hadley Centre Climate Model (HadCM3) and (2) Max Planck Institute’s ECHAM5/MPI-OM coupled climate model. Results from experiments forced with their twentieth century output are used to evaluate the modelled present-day ocean state. Sea ice coverage is largely realistic in both simulations; modelled ice shelf basal melt rates compare well with observations in both cases, but are consistently smaller for ECHAM5/MPI-OM. Projections for future ice shelf basal melting are computed using atmospheric output for the Intergovernmental Panel on Climate Change (IPCC) scenarios E1 and A1B. In simulations forced with ECHAM5 data, trends in ice shelf basal melting are small. In contrast, decreasing convection along the Antarctic coast in HadCM3 scenarios leads to a decreasing salinity on the continental shelf and to intrusions of warm deep water of open ocean origin. In the case of the Filchner–Ronne Ice Shelf (FRIS), this water reaches deep into the cavity, so that basal melting increases by a factor of 4 to 6 compared to the present value of about 90 Gt/year. By the middle of the twenty-second century, FRIS becomes the dominant contributor to total ice shelf basal mass loss in these simulations. Our results indicate that the surface freshwater fluxes on the continental shelves may be crucial for the future of especially the large cold water ice shelves in the Southern Ocean.     

The competing models for the origin and internal evolution of granitic pegmatites in the light of melt and fluid inclusion research

In this paper we discuss the main petrogenetic models for granitic pegmatites and how these models have evolved over time. We suggest that the present state of knowledge requires that some aspects of these models to be modified, or absorbed into newer ones. Pegmatite formation and internal evolution have long supposed the need for highly water- and flux-enriched magmas to explain the differences between pegmatites and other intrusives of similar major element composition. Compositions and textural characteristics of fluid and melt inclusions in pegmatite minerals provide strong evidence for such magmas. Furthermore, we show that melt inclusion research has increased the number of potential flux components, which may include H₂O, OH^?, CO₂, HCO ₃ ^? , CO ₃ ^2? , SO ₄ ^2? , PO ₄ ^3? , H₃BO₃, F , and Cl, as well as the elements Li, Na, K, Rb, Cs, and Be, herein described as melt structure modifiers. In this paper we emphasize that the combined effect which these components have on the properties of pegmatite melts is difficult to deduce from experimental studies using only a limited number of these components. The combination and the amount of the different magmatic species, together with differences in the source region, and variations in pressure and temperature cause the great diversity of the pegmatites observed. Some volatile species, such as CO ₃ ^2? and alkalis, have the capacity to increase the solubility of H₂O in silicate melt to an extraordinary degree, to the extent that melt-melt-fluid immiscibility becomes inevitable. It is our view that the formation of pegmatites is connected with the complex interplay of many factors.     

Grain boundary and volume diffusion experiments in yttrium aluminium garnet bicrystals at 1,723 K: a miniaturized study

Yb-Y inter-diffusion along a single grain boundary of a synthetic yttrium aluminium garnet (YAG) bicrystal has been studied using analytical transmission electron microscopy (ATEM). To investigate the diffusion, a thin-film containing Yb as the diffusant was deposited perpendicular to the bicrystal grain boundary by pulsed laser deposition (PLD). Structural properties and their change with time in both the diffusant source and the grain boundary are reported. The diffusion profiles are incorporated in a numerical diffusion model, which is applied to determine the grain boundary diffusion coefficient, D _gb , at 1.723 K it is equal to 3 × 10⁻¹⁵ m²/s. We find that grain boundary diffusion is 4.85 orders of magnitude faster than volume diffusion, which was determined from the same diffusion experiment. This result is discussed in the context of special versus general grain boundaries. Finally, we successfully tested the capability of synchrotron-based nano-X-ray fluorescence analysis to map micro-chemical patterns in two dimensions with sub-micrometre resolution.     

Growth of multilayered polycrystalline reaction rims in the MgO–SiO<Subscript>2</Subscript> system,part II: modelling

Part I of this contribution (Gardés et al. in Contrib Mineral Petrol, 2010) reported time- and temperature-dependent experimental growth of polycrystalline forsterite-enstatite double layers between single crystals of periclase and quartz, and enstatite single layers between forsterite and quartz. Both double and single layers displayed growth rates decreasing with time and pronounced grain coarsening. Here, a model is presented for the growth of the layers that couples grain boundary diffusion and grain coarsening to interpret the drop of the growth rates. It results that the growth of the layers is such that (Δx)² ∝ t ^1−1/n , where Δx is the layer thickness and n the grain coarsening exponent, as experimentally observed. It is shown that component transport occurs mainly by grain boundary diffusion and that the contribution of volume diffusion is negligible. Assuming a value of 1 nm for the effective grain boundary width, the following Arrhenius laws for MgO grain boundary diffusion are derived: log D _gb,0^Fo (m²/s) = −2.71 ± 1.03 and E _gb^Fo = 329 ± 30 kJ/mol in forsterite and log D _gb,0^En (m²/s) = 0.13 ± 1.31 and E _gb^En = 417 ± 38 kJ/mol in enstatite. The different activation energies are responsible for the changes in the enstatite/forsterite thickness ratio with varying temperature. We show that significant biases are introduced if grain boundary diffusion-controlled rim growth is modelled assuming constant bulk diffusivities so that differences in activation energies of more than 100 kJ/mol may arise. It is thus important to consider grain coarsening when modelling layered reaction zones because they are usually polycrystalline and controlled by grain boundary transport.     

Experimental growth of åkermanite reaction rims between wollastonite and monticellite: evidence for volume diffusion control

Growth rates of monomineralic, polycrystalline åkermanite (Ca₂MgSi₂O₇) rims produced by solid-state reactions between monticellite (CaMgSiO₄) and wollastonite (CaSiO₃) single crystals were determined at 0.5 GPa dry argon pressure, 1,000–1,200°C and 5 min to 60 h, using an internally heated pressure vessel. Inert Pt-markers, initially placed at the monticellite–wollastonite interface, indicate symmetrical growth into both directions. This and mass balance considerations demonstrate that rim growth is controlled by transport of MgO. At 1,200°C and run durations between 5 min and 60 h, rim growth follows a parabolic rate law with rim widths ranging from 0.4 to 16.3 μm indicating diffusion-controlled rim growth. The effective bulk diffusion coefficient \( D_{\text{eff,MgO}}^{\text{Ak}} \) is calculated to 10^?15.8±0.1 m² s^?1. Between 1,000°C and 1,200°C, the effective bulk diffusion coefficient follows an Arrhenius law with E _a = 204 ± 18 kJ/mol and D ₀ = 10^?8.6±1.6 m² s^?1. Åkermanite grains display a palisade texture with elongation perpendicular to the reaction interface. At 1,200°C, average grain widths measured normal to elongation, increase with the square root of time and range from 0.4 to 5.4 μm leading to a successive decrease in the grain boundary area fraction, which, however, does not affect \( D_{\text{eff,MgO}}^{\text{Ak}} \) to a detectible extent. This implies that grain boundary diffusion only accounts for a minor fraction of the overall chemical mass transfer, and rim growth is essentially controlled by volume diffusion. This is corroborated by the agreement between our estimates of the effective MgO bulk diffusion coefficient and experimentally determined volume diffusion data for Mg and O in åkermanite from the literature. There is sharp contrast to the MgO–SiO₂ binary system, where grain boundary diffusion controls rim growth.     

Provenance of Pleistocene Rhine River Middle Terrace sands between the Swiss–German border and Cologne based on U–Pb detrital zircon ages

Detrital zircon U–Pb age distributions derived from samples representing ancient or relatively young large-scale continental drainage networks are commonly taken to reflect the geochronological evolution of the tapped continental area. Here, we present detrital zircon U–Pb ages and associated heavy mineral data from Pleistocene Rhine River Middle Terrace sands and equivalents between the Swiss–German border and Cologne in order to test the commonly assumed Alpine provenance of the material. Samples from eight localities were analyzed for their heavy mineral assemblages. Detrital zircon U–Pb ages were determined by laser ablation inductively coupled mass spectrometry on selected samples from five locations along the Rhine River. The zircon age populations of all samples show a similar distribution, their main peaks being between 300 and 500 Ma. Minor age populations are recognized at 570 and 1,070 Ma. The 300–400 Ma maximum reflects the Variscan basement drained by or recycled into the Rhine River and its tributaries. The 400–500 Ma peak with predominantly Early Silurian ages points to Baltica or to the mid-German crystalline rise as original sources. One distinct peak at c. 570 Ma probably represents input from Cadomian terranes. The Precambrian U–Pb ages are compatible with derivation from sources in Baltica and in northern Gondwana. The heavy mineral populations of Middle Terrace sands and equivalents are characterized to a variable extend by garnet, epidote, and green hornblende. This association is often referred to as the Alpine spectrum and is considered to be indicative of an Alpine provenance. However, hornblende, epidote, and garnet are dominant heavy minerals of collisional orogens in general and may also be derived from Variscan and Caledonian units or from intermittent storage units. A remarkable feature of the detrital zircon age distribution in the Rhine River sediments from the Swiss–German border to Cologne is the absence of ages younger than 200 Ma and in particular of any ages reflecting the Alpine orogeny between c. 100 and 35 Ma. Sediments from rivers draining the equally collisional Himalaya orogen contain detrital zircons as young as 20 Ma. Our results question the assumption that Pleistocene Rhine River sediments were directly derived from the Alps. The lag time between the formation and deposition age of the youngest zircon in the studied Pleistocene Rhine River deposits is 200 Ma. Together with the absence of Alpine zircon ages, this stresses that detrital zircon age data from ancient sedimentary units found in poorly understood tectonic or paleogeographic settings need to be interpreted with great care, one could miss an entire orogenic cycle.     

Analytical Solution for Testing Debris Avalanche Numerical Models

—We present here the analytical solution of a one-dimensional dam-break problem over inclined planes. This solution is used to test a numerical model developed for debris avalanches. We consider a dam with infinite length in one direction where material is released from rest at the initial instant. We solve analytically and numerically the depth-averaged long-wave equations derived in a topography-linked coordinate system. The numerical and analytical solutions provide for a Coulomb-type friction law at the base of the flow. The analytical solution is obtained by using the method of characteristics and describes the flow over a constant slope, provided that the angle is higher than the friction angle. The numerical model utilizes a finite-difference method based on a Godunov-type scheme. Comparison between analytical and numerical results illustrates the remarkable stability and precision of the numerical method as well as its ability to deal with strong discontinuities.     

Seismic wave velocity and anisotropy of serpentinized peridotite in the Oman ophiolite

Shallow seismic measurements in harzburgite from the Oman ophiolite performed in a zone where the maximum horizontal anisotropy is expected (vertical foliation and horizontal lineation) point to a dominant dependence of seismic properties on fracturing.

Optical microscopy studies show that microcracks are guided by the serpentine (lizardite) penetrative network oriented subparallel to the harzburgite foliation and subperpendicular to the mineral lineation, and that serpentine (lizardite) vein filling has a maximum concentration of (001) planes parallel to the veins walls. The calculated elastic properties of the oriented alteration veins filled with serpentine in an anisotropic matrix formed by oriented crystals of olivine and orthopyroxene are compared with seismic velocities measured on hand specimens.

Laboratory ultrasonic data indicate that open microcracks are closed at 100 MPa pressure, e.g. (J. Geophys. Res. 65, (1960) 1083) and (Proc. ODP Sci. Results Leg 118, (1990) 227). Above this pressure, laboratory measurements and modeling show that P-compressional and S-shear wave velocities are mainly controlled by the mineral preferred orientation. Veins sealed with serpentine are effective in slightly lowering P and S velocities and increasing anisotropy. The penetrative lizardite network does not affect directly the geometry of seismic anisotropy, but contributes indirectly in the fact that this network controls the microcrack orientations.

Comparison between seismic measurements of peridotite and gabbro in the same conditions suggest that P- and S-waves anisotropies are a possible discriminating factor between the two lithologies in the suboceanic lithosphere.