Zircon U–Pb age, trace element and Hf isotope composition of Kongling terrane in the Yangtze Craton: refining the timing of Palaeoproterozoic high-grade metamorphism

U–Pb age, trace element and Hf isotope compositions of zircon were analysed for a metasedimentary rock and two amphibolites from the Kongling terrane in the northern part of the Yangtze Craton. The zircon shows distinct morphological and chemical characteristics. Most zircon in an amphibolite shows oscillatory zoning, high Th/U and ¹⁷⁶Lu/¹⁷⁷Hf ratios, high formation temperature, high trace element contents, clear negative Eu anomaly, as well as HREE-enriched patterns, suggesting that it is igneous. The zircon yields a weighted mean ²⁰⁷Pb/²⁰⁶Pb age of 2857 ± 8 Ma, representing the age of the magmatic protolith. The zircon in the other two samples is metamorphic. It has low Th/U ratios, low trace element concentrations, variable HREE contents (33.8 ≥ Lu_N≥2213; 14.7 ≤ Lu_N/Sm_N ≤ 354) and ¹⁷⁶Lu/¹⁷⁷Hf ratios (0.000030–0.001168). The data indicate that the zircon formed in the presence of garnet and under upper amphibolite facies conditions. The metamorphic zircon yields a weighted mean ²⁰⁷Pb/²⁰⁶Pb age of 2010 ± 13 Ma. These results combined with previously obtained Palaeoproterozoic metamorphic ages suggest a c. 2.0 Ga Palaeoproterozoic collisional event in the Yangtze Craton, which may result from the assembly of the supercontinent Columbia. The zircon in two samples yields weighted mean two-stage Hf model ( T _DM2) ages of 3217 ± 110 and 2943 ± 50 Ma, respectively, indicating that their protoliths were mainly derived from Archean crust.     

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.     

Isotope fractionation of neon during stepheating extraction?: a comment on ‘Re‐interpretation of the existence of a primitive plume under Australia based on neon isotope fractionation during step heating’ by Gautheron and Moreira (2003)

Using a Rayleigh distillation fractionation model, we calculate that the maximum isotope fractionation potentially achievable is less than 5% during the early stages of gas release from a sample. Our calculation corrects the erroneous conclusions of Gautheron and Moreira (2003), who re‐interpreted the plume‐like neon isotopic compositions found in metasomatic apatite from a south‐eastern Australian xenolith (Matsumoto et al., 1997) to be the result of Rayleigh‐type isotope fractionation of originally MORB‐type neon during stepheating gas extraction. We stress that the modelling of neon isotopic fractionation by Gautheron and Moreira (2003) is incorrect, and that the finding of a plume‐like neon isotopic composition in the apatite by Matsumoto et al. (1997) remains a quite valid and robust conclusion.     

A critical state model for sands dependent on stress and density

An elastoplastic model for sands is presented in this paper, which can describe stress–strain behaviour dependent on mean effective stress level and void ratio. The main features of the proposed model are: (a) a new state parameter, which is dependent on the initial void ratio and initial mean stress, is proposed and applied to the yield function in order to predict the plastic deformation for very loose sands; and (b) another new state parameter, which is used to determine the peak strength and describe the critical state behaviour of sands during shearing, is proposed in order to predict simply negative/positive dilatancy and the hardening/softening behaviour of medium or dense sands. In addition, the proposed model can also predict the stress–strain behaviour of sands under three-dimensional stress conditions by using a transformed stress tensor instead of ordinary stress tensor. Copyright © 2004 John Wiley & Sons, Ltd.     

The Siberian Solar Radio Telescope: the current state of the instrument,observations, and data

The Siberian Solar Radio Telescope (SSRT) is one of the world's largest solar radio heliographs. It commenced operation in 1983, and since then has undergone several upgrades. The operating frequency of the SSRT is 5.7 GHz. Since 1992 the instrument has had the capability to make one-dimensional scans with a high time resolution of 56 ms and an angular resolution of 15 arc sec. Making one of these scans now takes 14 ms. In 1996 the capability was added to make full, two-dimensional images of the solar disk. The SSRT is now capable of obtaining images with an angular resolution of 21 arc sec every 2 min. In this paper we describe the main features and operation of the instrument, particularly emphasizing issues pertaining to the imaging process and factors limiting data quality. Some of the data processing and analysis techniques are discussed. We present examples of full-disk solar images of the quiet Sun, recorded near solar activity minimum, and images of specific structures: plages, coronal bright points, filaments and prominences, and coronal holes. We also present some observations of dynamic phenomena, such as eruptive prominences and solar flares, which illustrate the high-time-resolution observations that can be done with this instrument. We compare SSRT observations at 5.7 GHz, including computed `light curves', both morphologically and quantatively, with observations made in other spectral domains, such as 17 GHz radio images, Hα filtergrams and magnetograms, extreme-ultraviolet and X-ray observations, and dynamic radio spectra.     

Bond rolling resistance and its effect on yielding of bonded granulates by DEM analyses

A discrete element modelling of bonded granulates and investigation on the bond effect on their behaviour are very important to geomechanics. This paper presents a two‐dimensional (2‐D) discrete element theory for bonded granulates with bond rolling resistance and provides a numerical investigation into the effect of bond rolling resistance on the yielding of bonded granulates. The model consists of mechanical contact models and equations governing the motion of bonded particles. The key point of the theory is that the assumption in the original bond contact model previously proposed by the authors (55th CSCE‐ASCE Conference, Hamilton, Ont., Canada, 2002; 313–320; J. Eng. Mech. (ASCE) 2005; 131 (11):1209–1213) that bonded particles are in contact at discrete points, is here replaced by a more reliable assumption that bonded particles are in contact over a width. By making the idealization that the bond contact width is continuously distributed with the normal/tangential basic elements (BE) (each BE is composed of spring, dashpot, bond, slider or divider), we establish a bond rolling contact model together with bond normal/tangential contact models, and also relate the governing equations to local equilibrium. Only one physical parameter β needs to be introduced in the theory in comparison to the original bond discrete element model. The model has been implemented into a 2‐D distinct element method code, NS2D. Using the NS2D, a total of 86 1‐D, constant stress ratio, and biaxial compressions tests have been carried out on the bonded granular samples of different densities, bonding strengths and rolling resistances. The numerical results show that: (i) the new theory predicts a larger internal friction angle, a larger yielding stress, more brittle behaviour and larger final broken contact ratio than the original bond model; (ii) the yielding stress increases nonlinearly with the increasing value of β, and (iii) the first‐yield curve (initiation of bond breakage), which define a zone of none bond breakage and which shape and size are affected by the material density, is amplified by the bond rolling resistance in analogous to that predicted by the original bond model. Copyright © 2006 John Wiley & Sons, Ltd.     

A multi‐level parallelized substructuring‐frontal solution for coupled thermo/hydro/mechanical problems in unsaturated soil

This paper proposes a multi‐level parallelized substructuring–frontal combined algorithm for the analysis of the problem of thermo/hydraulic/mechanical behaviour of unsaturated soil. Temperature, displacement, pore water pressure and pore air pressure are treated as the primary variables in a non‐linear analysis. Details are given firstly of the substructuring–frontal combined approach. The incorporation of the algorithm in a multi‐level parallel strategy is then discussed. The parallel processing can thus be carried out at different substructural levels. The method thus developed impacts, in a positive way, on both computer storage requirement and execution time. Copyright © 2003 John Wiley & Sons, Ltd.     

On the most general accurate solutions for Buchdahl's fluid spheres

The most general accurate solutions for the Buchdhal fluid sphere were obtained and matched with the Schwarzchild's exterior solution at the pressure free interface. Various parameters of the solutions were so adjusted that the energy density, pressure and temperature were positive and decreasing away from the centre, and the velocity of sound was less than unity throughout the spheroid model. Using this procedure, the maximum mass of the fluid sphere with a surface density of2 × 10¹⁴gm cm^-3 was determined to be 3.82MΘ and 4.57 MΘ for strong and weak energy conditions respectively. PACS number: 0402, 0402J, 0440D, 95301 This revised version was published online in July 2006 with corrections to the Cover Date.     

Vertical stress distribution due to an arbitrarily shaped foundation using triangulation and an analytical expression for a triangular loaded region

Using the basic Boussinesq's equation, the expression for the vertical stress distribution (σ_z) underneath any point on the ground surface due to a general triangular loaded region in a preferred orientation with a linearly varied loading has been successfully derived. When the triangle is not in a preferred orientation, a simple axis transformation is required and the expression will be equally applicable. Based on this expression, σ_z due to an arbitrarily shaped loaded foundation can simply be determined by first triangulating the loaded area and summing up the contributions from each generated triangular region. The procedures for triangulating and calculating the stress distribution can be simply automated through computer programs.