Application of micropolar plasticity to post failure analysis in geomechanics

A micropolar elastoplastic model for soils is formulated and a series of finite element analyses are employed to demonstrate the use of a micropolar continuum in overcoming the numerical difficulties encountered in application of finite element method in standard Cauchy–Boltzmann continuum. Three examples of failure analysis involving a deep excavation, shallow foundation, and a retaining wall are presented. In all these cases, it is observed that the length scale introduced in the polar continuum regularizes the incremental boundary value problem and allows the numerical simulation to be continued until a clear collapse mechanism is achieved. The issue of grain size effect is also discussed. Copyright © 2004 John Wiley & Sons, Ltd.     

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.     

Comet nuclei: Morphology and implied processes of surface modification

Surface morphology and related issues for nuclei of three comets: Halley, Borrelly and Wild 2, are considered in the paper. Joint consideration of publications and results of our analysis of the comets’ images led to conclusions, partly new, partly repeating conclusions published by other researchers. It was found that typical for all three nuclei is the presence of rather flat areas: floors of craters and other depressions, mesas and terraces. This implies that flattening surfaces or planation is a process typical for the comet nuclei. Planation seems to work through the sublimation-driven slope collapse and retreat. This requires effective sublimation so this process should work only when a comet is close to the Sun and if on the nucleus there are starting slopes, steep and high enough to support the “long-distance” avalanching of the collapsing material. If the surface had no starting slopes, then instead of planation, the formation of pitted-and-hilly surfaces should occur. An example of this could be the mottled terrain of the Borelly nucleus. Both ways of the sublimational evolution on the nucleus surface should lead to accumulation of cometary regolith. The thickness of the degassed regolith is not known, but it is obvious that in surface depressions, including the flat-floor ones, it should be larger compared with nondepression areas. This may have implications for the in situ study of comets by the Deep Impact and Rosetta missions.Our morphological analysis puts constraints on the applicability of the popular “rubble-pile comet nucleus” hypothesis (Weissman, 1986. Are cometery nuclei primordial rubble piles? Nature 320, 242-244.). We believe that the rubble pile hypothesis can be applicable to the blocky Halley nucleus. The Borelly and Wild 2 nuclei also could be rubble piles. But in these cases the “rubbles” have to be either smaller than 30-50 m (a requirement to keep lineament geometry close to ideal), or larger than 1-2 km (a requirement to form the rather extended smooth, flat surfaces of mesa tops and crater floors). Another option is that the Borelly and Wild 2 nuclei are not rubble piles.In relation to surface morphology we suggest that three end-member types of the comet nuclei may exist: (1) impact cratered “pristine” bodies, (2) non-cratered fragments of catastrophic disruption, and (3) highly Sun-ablated bodies. In this threefold classification, the Wild 2 nucleus is partially ablated primarily cratered body. Borrelly is significantly ablated and could be either primarily cratered or not-cratered fragment. Halley is certainly partially ablated but with the available images it is difficult to say if remnants of impact craters do exist on it.Recently published observations and early results of analysis of the Tempel 1 nucleus images taken by Deep Impact mission are in agreement with our conclusions on the processes responsible for the Halley, Borrelly and Wild 2 nuclei morphologies. In particular, we have now more grounds to suggest that decrease in crater numbers and increase of the role of smooth flat surfaces in the sequence Wild 2?Tempel 1?Borelli reflects a progress in the sublimational degradation of the nucleus surface during comet passages close to the Sun.     

Isotope hydrological study of mean transit time in the granitic Strengbach catchment (Vosges massif,France): application of the FlowPC model with modified input function

Measurements of ¹⁸O concentrations in precipitation, soil solution, spring and runoff are used to determine water transit time in the small granitic Strengbach catchment (0·8 km²; 883–1146 m above sea level) located in the Vosges Mountains of northeastern France. Water transit times were calculated by applying the exponential, exponential piston and dispersion models of the FlowPC program to isotopic input (rainfall) and output (spring and stream water) data sets during the period 1989–95. The input function of the model was modified compared with the former version of the model and estimated by a deterministic approach based on a simplified hydrological balance. The fit between observed and calculated output data showed marked improvements compared with results obtained using the initial version of the model. An exponential piston version of the model applied to spring water indicates a 38·5 month mean transit time, which suggests that the volume in the aquifer, expressed in water depth, is 2·4 m. A considerable thickness (>45 m) of fractured bedrock may be involved for such a volume of water to be stored in the aquifer. Copyright © 2005 John Wiley & Sons, Ltd.     

Green Lake Landslide and other giant and very large postglacial landslides in Fiordland,New Zealand

Green Lake Landslide is an ancient giant rock slide in gneiss and granodiorite located in the deeply glaciated Fiordland region of New Zealand. The landslide covers an area of 45 km² and has a volume of about 27 km³. It is believed to be New Zealand's largest landslide, and possibly the largest landslide of its type on Earth. It is one of 39 known very large (10⁶–10⁷ m³) and giant (≥10⁸ m³) postglacial landslides in Fiordland discussed in the paper. Green Lake Landslide resulted in the collapse of a 9 km segment of the southern Hunter Mountains. Slide debris moved up to 2.5 km laterally and 700 m vertically, and formed a landslide dam about 800 m high, impounding a lake about 11 km long that was eventually infilled with sediments. Geomorphic evidence supported by radiocarbon dating indicates that Green Lake Landslide probably occurred 12 000–13 000 years ago, near the end of the last (Otira) glaciation. The landslide is described, and its geomorphic significance, age, failure mechanism, cause, and relevance in the region are discussed, in relation to other large landslides and recent earthquake-induced landslides in Fiordland. The slope failure occurred on a low-angle fault zone undercut by glacial erosion, and was probably triggered by strong shaking (MM IX–X) associated with a large (≥ M 7.5–8) earthquake, on the Alpine Fault c. 80 km to the northwest. Geology was a major factor that controlled the style and size of Green Lake landslide, and in that respect it is significantly different from most other gigantic landslides. Future large earthquakes on the Alpine Fault in Fiordland are likely to trigger more very large and giant landslides across the region, causing ground damage and devastation on a scale that has not occurred during the last 160 years, with potentially disastrous effects on towns, tourist centres, roads, and infrastructure. The probability of such an event occurring within the next 50 years may be as high as 45%.     

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.     

Porphyroblast inclusion-trail orientation data: eppure non son girate*!

Extensive examination of large numbers of spatially orientated thin sections of orientated samples from orogens of all ages around the world has demonstrated that porphyroblasts do not rotate relative to geographical coordinates during highly non-coaxial ductile deformation of the matrix subsequent to their growth. This has been demonstrated for all tectonic environments so far investigated. The work also has provided new insights and data on metamorphic, structural and tectonic processes including: (1) the intimate control of deformation partitioning on metamorphic reactions; (2) solutions to the lack of correlation between lineations that indicate the direction of movement within thrusts and shear zones, and relative plate motion; and (3) a possible technique for determining the direction of relative plate motion that caused orogenesis in ancient orogens.