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.     

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.     

Observational evidence for solar coronal decimetric radio pulsations with very short periods

Using the decimetric (700–1500 MHz) radio spectrometer and the synchronous observational system with high temporal resolution at four frequencies (1420, 2130, 2840 and 4260 MHz) of Yunnan Observatory, two rare events were observed on 2001 June 24 and 1990 July 30. The former was a small radio burst exhibiting pulsations with short periods (about 29, 40 and 100 ms) in the impulsive phase. The latter was a large radio burst, which at 2840 MHz produced radio pulsations with period of about 30 ms. This paper focuses on pulsations with very short periods in the range of 29–40 ms. The mechanism of generation of such pulsations may be modulation of radio radiation by the periodic trains of whistler packets originating in unstable regions of the corona. Alternatively, these pulsations can be attributed to wave-wave non-linear interactions of electrostatic upper hybrid waves driven by beams of precipitating electrons in flaring loops.     

Prediction Test for the Two Extremely Strong Solar Storms in October 2003

In late October and early November 2003, a series of space weather hazard events erupted in solar-terrestrial space. Aiming at two intense storm (shock) events on 28 and 29 October, this paper presents a Two-Step method, which combines synoptic analysis of space weather–`observing’ and quantitative prediction – ‘palpating’, and uses it to test predictions. In the first step, ‘observing’, on the basis of observations of the source surface magnetic field, interplanetary scintillation (IPS) and ACE spacecraft, we find that the propagation of the shock waves is asymmetric and northward relative to the normal direction of their solar sources due to the large-scale configuration of the coronal magnetic fields, and the Earth is located near the direction of the fastest speed and greatest energy of the shocks. Being two fast ejection shock events, the fast explosion of extremely high temperature and strong magnetic field, and background solar wind velocity as high as 600 and 1000 km s⁻¹, are also helpful to their rapid propagation. According to the synoptic analysis, the shock travel times can be estimated as 21 and 20 h, which are close to the observational results of 19.97 and 19.63 h, respectively. In the second step, ‘palpating’, we adopt a new membership function of the fast shock events for the ISF method. The predicted results here show that for the onset time of the geomagnetic disturbance, the relative errors between the observational and the predicted results are 1.8 and 6.7%, which are consistent with the estimated results of the first step; and for the magnetic disturbance magnitude, the relative errors between the observational and the predicted results are 4.1 and 3.1%, respectively. Furthermore, the comparison among the predicted results of our Two-Step method with those of five other prevailing methods shows that the Two-Step method is advantageous in predicting such strong shock event. It can predict not only shock arrival time, but also the magnitude of magnetic disturbance. The results of the present paper tell us that understanding the physical features of shock propagation thoroughly is of great importance in improving the prediction efficiency.     

Automated analysis of eclipsing binary light curves – II. Statistical analysis of OGLE LMC eclipsing binaries

In the first paper of this series, we presented EBAS – Eclipsing Binary Automated Solver, a new fully automated algorithm to analyse the light curves of eclipsing binaries, based on the ebop code. Here, we apply the new algorithm to the whole sample of 2580 binaries found in the Optical Gravitational Lensing Experiment (OGLE) Large Magellanic Cloud (LMC) photometric survey and derive the orbital elements for 1931 systems. To obtain the statistical properties of the short-period binaries of the LMC, we construct a well-defined subsample of 938 eclipsing binaries with main-sequence B-type primaries. Correcting for observational selection effects, we derive the distributions of the fractional radii of the two components and their sum, the brightness ratios and the periods of the short-period binaries. Somewhat surprisingly, the results are consistent with a flat distribution in log P between 2 and 10 d. We also estimate the total number of binaries in the LMC with the same characteristics, and not only the eclipsing binaries, to be about 5000. This figure leads us to suggest that  (0.7 ± 0.4)  per cent of the main-sequence B-type stars in the LMC are found in binaries with periods shorter than 10 d. This frequency is substantially smaller than the fraction of binaries found by small Galactic radial-velocity surveys of B stars. On the other hand, the binary frequency found by Hubble Space Telescope ( HST ) photometric searches within the late main-sequence stars of 47 Tuc is only slightly higher and still consistent with the frequency we deduced for the B stars in the LMC.     

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.     

OSIRIS: AO-assisted integral-field spectroscopy at the Keck Observatory

OSIRIS (OH-Suppressing Infra-Red Integral-field Spectrograph) is a new facility instrument for the Keck Observatory. After seeing first light in February 2005, OSIRIS is currently undergoing commissioning. OSIRIS provides the capability of performing three-dimensional spectroscopy in the near-infrared z, J, H, and K bands at the resolution limit of the Keck II telescope, which is equipped with adaptive optics and a laser guide star. The science case for OSIRIS is summarized, and the instrument and associated data reduction software are described.     

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%.