Age and early metamorphic history of the Sanbagawa belt: Lu–Hf and P–T constraints from the Western Iratsu eclogite

Two distinct age estimates for eclogite-facies metamorphism in the Sanbagawa belt have been proposed: (i) c.  120–110 Ma based on a zircon SHRIMP age for the Western Iratsu unit and (ii) c.  88–89 Ma based on a garnet–omphacite Lu–Hf isochron age from the Seba and Kotsu eclogite units. Despite the contrasting estimates of formation ages, petrological studies suggest the formation conditions of the Western Iratsu unit are indistinguishable from those of the other two units—all ∼20 kbar and 600–650 °C. Studies of the associated geological structures suggest the Seba and Western Iratsu units are parts of a larger semi-continuous eclogite unit. A combination of geochronological and petrological studies for the Western Iratsu eclogite offers a resolution to this discrepancy in age estimates. New Lu–Hf dating for the Western Iratsu eclogite yields an age of 115.9 ± 0.5 Ma that is compatible with the zircon SHRIMP age. However, petrological studies show that there was significant garnet growth in the Western Iratsu eclogite before eclogite facies metamorphism, and the early core growth is associated with a strong concentration of Lu. Pre-eclogite facies garnet (Grt1) includes epidote–amphibolite facies parageneses equilibrated at 550–650 °C and ∼10 kbar, and this is overgrown by prograde eclogite facies garnet (Grt2). The Lu–Hf age of c.  116 Ma is strongly skewed to the isotopic composition of Grt1 and is interpreted to reflect the age of the pre-eclogite phase. The considerable time gap ( c.  27 Myr) between the two Lu–Hf ages suggests they may be related to separate tectonic events or distinct phases in the evolution of the Sanbagawa subduction zone.     

Pressure-induced incipient amorphization of α-quartz and transition to coesite in an eclogite from Antarctica: a first record and some consequences

Monocrystalline quartz inclusions in garnet and omphacite from various eclogite samples from the Lanterman Range (Northern Victoria Land, Antarctica) have been investigated by cathodoluminescence (CL), Raman spectroscopy and imaging, and in situ X‐ray (XR) microdiffraction using the synchrotron. A few inclusions, with a clear‐to‐opalescent lustre, show ‘anomalous’ Raman spectra characterized by weak α‐quartz modes, the broadening of the main α‐quartz peak at 465 cm^?1, and additional vibrations at 480–485, 520–523 and 608 cm^?1. CL and Raman imaging indicate that this ‘anomalous’α‐quartz occurs as relicts within ordinary α‐quartz, and that it was preserved in the internal parts of small quartz inclusions. XR diffraction circular patterns display irregular and broad α‐quartz spots, some of which show an anomalous d‐spacing tightening of ～2%. They also show some very weak, hazy clouds that have d‐spacing compatible with coesite but not with α‐quartz. Raman spectrometry and XR microdiffraction characterize the anomalies with respect to α‐quartz as (i) a pressure‐induced disordering and incipient amorphization, mainly revealed by the 480–485 and 608‐cm^?1 Raman bands, together with (ii) a lattice densification, evidenced by d‐spacing tightening; (iii) the cryptic development of coesite, 520–523 cm^?1 being the main Raman peak of coesite and (iv) Brazil micro‐twinning. This ‘anomalous’α‐quartz represents the first example of pressure‐induced incipient amorphization of a metastable phase in a crustal rock. This issue is really surprising because pressure‐induced amorphization of metastable α‐quartz, observed in impactites and known to occur between 15 and 32 GPa during ultrahigh‐pressure (UHP) experiments at room temperature, is in principle irrelevant under normal geological P–T conditions. A shock (due to a seism?) or a local overpressure at the inclusion scale (due to expansion mismatch between quartz and its host mineral) seem the only geological mechanisms that can produce such incipient amorphization in crustal rocks. This discovery throws new light on the modality of the quartz‐coesite transition and on the pressure regimes (non‐lithostatic v. lithostatic) during high‐pressure/UHP metamorphism. In particular, incipient amorphization of quartz could favour the quartz‐coesite transition, or allow the growth of metastable coesite, as already experimentally observed.     

Interplay between erosion and tectonics in the Western Alps

Bedrock fission‐track analysis, high‐resolution petrography and heavy mineral analyses of sediments are used to investigate the relationships between erosion and tectonics in the Western Alps. Along the Aosta Valley cross‐section, exhumation rates based on fission‐track data are higher in the fault‐bounded western block than in the eastern block (0.4–1.5 vs. 0.1–0.3 mm yr⁻¹). Erosion rates based on the analysis of bed‐load in the Dora Baltea drainage display the same pattern and have similar magnitudes in the relative sub‐basins (0.4–0.7 vs. 0.04–0.08 mm yr⁻¹). Results highlight that climate, relief and lithology are not the controlling factors of erosion in the Western Alps. The main driving force behind erosion is instead tectonics that causes the differential upward motion of crustal blocks.     

Self-consistent response of a galactic disc to vertical perturbations

We study the self-consistent, linear response of a galactic disc to vertical perturbations, as induced, say, by a tidal interaction. We calculate the self-gravitational potential corresponding to a non-axisymmetric, self-consistent density response of the disc using the Green's function method. The response potential is shown to oppose the perturbation potential because the self-gravity of the disc resists the imposed potential, and this resistance is stronger in the inner parts of a galactic disc. For the   m = 1  azimuthal wavenumber, the disc response opposes the imposed perturbation up to a radius that spans a range of 4–6 disc scalelengths, so that the disc shows a net warp only beyond this region. This physically explains the well known but so far unexplained observation that warps typically set in beyond this range of radii. We show that the inclusion of a dark matter halo in the calculation only marginally changes (by ∼10 per cent) the radius for the onset of warps. For perturbations with higher azimuthal wavenumbers, the net signature of the vertical perturbations can only be seen at larger radii – for example, beyond 7 exponential disc scalelengths for   m = 10  . Also, for the high- m cases, the magnitude of the negative disc response due to the disc self-gravity is much smaller. This is shown to result in corrugations of the mid-plane density, which explains the puzzling scalloping with   m = 10  detected in H  i in the outermost regions ∼30 kpc in the Galaxy.     

Asymmetry of 56Ni ejecta and polarization in the type-IIP supernova 2004dj

I study the question of whether the asymmetry of ⁵⁶Ni ejecta that results in the asymmetry of the Hα emission line at the nebular epoch of the type-IIP supernova SN 2004dj can account for the recently detected polarization of the supernova radiation. I have developed a model of the Hα profile and luminosity with nonthermal ionization and excitation in a spherically symmetric envelope for an asymmetric bipolar ⁵⁶Ni distribution. I have calculated the polarized radiation transfer against the background of the recovered electron density distribution. The observed polarization is shown to be reproduced at the nebular epoch around day 140 for the same parameters of the envelope, and the ⁵⁶Ni distribution for which the evolution of the Hα luminosity and profile is explained. Yet the model polarization decreases with time more slowly than is observed. The origin of the additional component responsible for the early polarization on day 107 is discussed.     

P‐wave transmission across fractures with nonlinear deformational behaviour

Stress wave attenuation across fractured rock masses is a great concern of underground structure safety. When the wave amplitude is large, fractures experience nonlinear deformation during the wave propagation. This paper presents a study on normal transmission of P‐wave across parallel fractures with nonlinear deformational behaviour (static Barton–Bandis model). The results show that the magnitude of transmission coefficient is a function of incident wave amplitude, nondimensional fracture spacing and number of fractures. Two important indices of nondimensional fracture spacing are identified, and they divide the area of nondimensional fracture spacing into three parts (individual fracture area, transition area and small spacing area). In the different areas, the magnitude of transmission coefficient has different trends with nondimensional fracture spacing and number of fractures. In addition, the study reveals that under some circumstances, the magnitude of transmission coefficient increases with increasing number of fractures, and is larger than 1. Copyright © 2006 John Wiley & Sons, Ltd.     

Thick and thin models of the evolution of carbon dioxide on Mars

We present results from a new simulation code that accounts for the evolution of the reservoirs of carbon dioxide on Mars, from its early years to the present. We establish a baseline model parameter set that produces results compatible with the present (i.e., Patm?6.5 mbar with permanent CO₂ ice cap) for a wide range of initial inventories. We find that the initial inventory of CO₂ broadly determines the evolutionary course of the reservoirs of CO₂. The reservoirs include the atmosphere, ice cap, adsorbed CO₂ in the regolith, and carbonate rocks. We track the evolution of the free inventory: the atmosphere, ice cap and regolith. Simulations begin at 4.53 Gyr before present with a rapid loss of free inventory to space in the early Noachian. Models that assume a relatively small initial inventory (?5 bar) have pronounced minima in the free inventory of CO₂ toward the end of the Noachian. Under baseline parameters, initial inventories below ∼4.5 bar result in a catastrophic loss of the free inventory to space. The current free inventory would be then determined by the balance between outgassing, sputtering losses and chemical weathering following the end of the late bombardment. We call these “thin” models. They generically predict small current free inventories in line with expectations of a small present CO₂ ice cap. For “thick” models, with initial inventories ?5 bar, a surplus of 300-700 mbar of free CO₂ remains during the late-Noachian. The histories of free inventory in time for thick models tend to converge within the last 3.5 Gyr toward a present with an ice cap plus atmospheric inventory of about 100 mbar. For thick models, the convergence is largely due to the effects of chemical weathering, which draws down higher free inventories more rapidly than the low. Thus, thick models have ?450 mbar carbonate reservoirs, while thin models have ?200 mbar. Though both thick and thin scenarios can reproduce the current atmospheric pressure, the thick models imply a relatively large current CO₂ ice cap and thin models, little or none. While the sublimation of a massive cap at a high obliquity would create a climate swing of greenhouse warming for thick models, under the thin model, mean temperatures and pressures would be essentially unaffected by increases in obliquity.