Numerical manifold method for dynamic nonlinear analysis of saturated porous media

It is well known that the Babuska–Brezzi stability criterion or the Zienkiewicz–Taylor patch test precludes the use of the finite elements with the same low order of interpolation for displacement and pore pressure in the nearly incompressible and undrained cases, unless some stabilization techniques are introduced for dynamic analysis of saturated porous medium where coupling occurs between the displacement of solid skeleton and pore pressure. The numerical manifold method (NMM), where the interpolation of displacement and pressure can be determined independently in an element for the solution of u–p formulation, is derived based on triangular mesh for the requirement of high accurate calculations from practical applications in the dynamic analysis of saturated porous materials. The matrices of equilibrium equations for the second‐order displacement and the first‐order pressure manifold method are given in detail for program coding. By close comparison with widely used finite element method, the NMM presents good stability for the coupling problems, particularly in the nearly incompressible and undrained cases. Numerical examples are given to illustrate the validity and stability of the manifold element developed. Copyright © 2006 John Wiley & Sons, Ltd.     

A major sulphur isotope event at c. 510 Ma: a possible anoxia–extinction–volcanism connection during the Early–Middle Cambrian transition?

A new approach to constraining seawater δ³⁴S and sulphate concentration using francolite‐bound sulphate reveals an abrupt increase in δ³⁴S to +50‰ around the Early–Middle Cambrian boundary. Such high δ³⁴S values are best explained by increased rates of pyrite burial due to ocean anoxia coupled with an increased sensitivity of the ocean sulphate reservoir to perturbations due to low sulphate concentrations of 500–700 μgL^?1. We argue that the spread of anoxic waters at this time was partly the result of greenhouse warming related to the eruption of the Kalkarindji Large Igneous Province of northern Australia and that it triggered the collapse of early metazoan reef ecosystems during the latest Early Cambrian. Mass extinctions of the last 260 Myr have all coincided with enhanced volcanic activity, while several are also associated with positive shifts in seawater δ³⁴S. Extending this correlation back in time further implicates volcanically induced climate change as a major determining factor in biosphere evolution. Terra Nova, 18, 257–263, 2006     

Mezen Bay—a key area for understanding Weichselian glaciations in northern Russia

Sediment successions in coastal cliffs around Mezen Bay, southeastern White Sea, record an unusually detailed history of former glaciations, interstadial marine and fluvial events from the Weichselian. A regional glaciation model for the Weichselian is based on new data from the Mezen Bay area and previously published data from adjacent areas. Following the Mikulinian (Eemian) interglacial a shelf‐centred glaciation in the Kara Sea is reflected in proglacial conditions at 100–90 ka. A local ice‐cap over the Timan ridge existed between 75 and 65 ka. Renewed glaciation in the Kara Sea spread southwestwards around 60 ka only, interrupted by a marine inundation, before it advanced to its maximum position at about 55–50 ka. After a prolonged ice‐free period, the Scandinavian ice‐sheet invaded the area from the west and terminated east of Mezen Bay about 17 ka. The previously published evidence of a large ice‐dammed lake in the central Arkhangelsk region, Lake Komi, finds no support in this study. Copyright © 2003 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.     

Multistep matched-field inversion for broad-band data from ASIAEX2001

This paper discusses the results of geoacoustic inversion carried out using explosive charge data from the Asian Seas International Acoustic Experiment (ASIAEX) East China Sea (ECS) Experiment. A multifrequency incoherent matched-field inversion processor and a genetic algorithm (GA) are used for the inversion. A multistep matched field inversion approach is presented, which makes use of the varying sensitivities of wave fields at various frequencies to reduce the inversion problem into a sequence of smaller inversions with fewer unknowns to estimate at each stage. Different parameters are estimated using data at different frequencies according to their sensitivities. Inversion results for different areas in the ECS region are summarized and compared with core data.     

A one-dimensional model for the interaction between continental-scale ice sheets and atmospheric stationary waves

 The great continental ice sheets of the Pleistocene represented a significant topographic obstacle to the westerly winds in northern midlatitudes. This work explores how consequent changes in the atmospheric stationary wave pattern might have affected the shape and growth of the ice sheets themselves. A one dimensional (1-D) model is developed which permits an examination of the types and magnitudes of the feedbacks that might be expected. When plausible temperature perturbations are introduced at the ice-sheet margin which are proportional to the stationary wave amplitude, the equilibrium shape of the ice sheet is significantly altered, and depends on the sign of the perturbation. The proposed feedback also affects the response of the ice sheet to time-varying climate forcing. The results suggest that the evolution of a continental-scale ice sheet with a land-based margin may be significantly determined by the changes it induces in the atmospheric circulation. Received: 1 October 1999 / Accepted: 17 July 2000     

Small angle neutron scattering by opals

Micro- and non-crystalline opals, chalcedony and flint show diffuse small angle neutron scattering (SANS). Precious opals give rise to two additional intensity maxima at very small scattering angles which are due to Bragg reflections from the closest packed non-crystalline silica spheres. A small angle texture diagram reveals that the closest packing is faulty. Synthetic non-crystalline opals yield much less intense small angle scattering due to lower contrast between silica spheres and interstitial cement or particles; in this case intensity maxima were not observed. The outer part of the scattering curves of opal-CT and microcrystalline quartz deviates from Porod's law. The specific surface of natural non-crystalline opals ranges from 0.006 to 0.018 nm^–1. In microcrystalline opals, the specific surface is about 10 times larger than in non-crystalline opals.     

Rift-related Jurassic basaltic phreatomagmatic volcanism in the central Transantarctic Mountains: precursory stage to flood-basalt effusion

 The Middle Jurassic Kirkpatrick flood basalts and comagmatic Ferrar intrusions in the Transantarctic Mountains represent a major pulse of tholeiitic magmatism related to early stages in the breakup of Gondwana. A record of the volcano-tectonic events leading to formation of this continental flood-basalt province is provided by strata underlying and only slightly predating the Kirkpatrick lavas. In the central Transantarctic Mountains, the lavas rest on widespread (≥7500 km²) tholeiitic pyroclastic deposits of the Prebble Formation. The Prebble Formation is dominated by lahar deposits and is an unusual example of a regionally developed basaltic lahar field. Related, partly fault-controlled pyroclastic intrusions cut underlying strata, and vents are represented by the preserved flanks of two small tephra cones associated with a volcanic neck. Lahar and air-fall deposits typically contain 50–60% accidental lithic fragments and sand grains derived from underlying Triassic – Lower Jurassic strata in the upper part of the Beacon Supergroup. Juvenile basaltic ash and fine lapilli consist of nonvesicular to scoriaceous tachylite, sideromelane, and palagonite, and have characteristics indicating derivation from hydrovolcanic eruptions. The abundance of accidental debris from underlying Beacon strata points to explosive phreatomagmatic interaction of basaltic magma with wet sediment and groundwater, which appears to have occurred in particular where rising magma intersected upper Beacon sand aquifers. Composite clasts in the lahar deposits exhibit complex peperitic textures formed during fine-scale intermixing of basaltic magma with wet sand and record steps in subsurface fuel-coolant interactions leading to explosive eruption. The widespread, sustained phreatomagmatic activity is inferred to have occurred in a groundwater-rich topographic basin linked to an evolving Jurassic rift zone in the Transantarctic Mountains. Coeval basaltic phreatomagmatic deposits of the Mawson and Exposure Hill Formations, which underlie exposures of the Kirkpatrick Basalt up to 1500 km to the north along strike in Victoria Land, appear to represent other parts of a regional, extension-related Middle Jurassic phreatomagmatic province which developed immediately prior to rapid outpouring of the flood basalts. This is consistent with models which assign an important role to lithospheric stretching in the generation of flood-basalt provinces. Received: 28 August 1995 / Accepted: 18 April 1996     

An experimental study in forming a long synthetic aperture at sea

The testing of a synthetic aperture technique, the ETAM algorithm, is extended and its performance for CW pseudorandom signals and broadband ship noise is examined. The results show the limitations of the technique and are of special interest for operational systems development. In the CW experiments, the transmitted signal was generated with high temporal coherence, and loss of the spatial and temporal coherence of the received signal was introduced only by the medium and the stability of the towed array. In the experiments that included the pseudorandom signal and the ship noise, the temporal coherence of the transmitted signals was deliberately chosen to be poor in order to study the effects and the performance of this algorithm with broadband signals. The related experimental results show that for received signals, which have their segments over the synthesizing period highly cross-correlated, a synthetic aperture array gain was achieved that corresponds to the length of an equivalent fully populated array