A continuum‐discrete hydromechanical analysis of granular deposit liquefaction

A coupled continuum‐discrete hydromechanical model was employed to analyse the liquefaction of a saturated loose deposit of cohesionless particles when subjected to a dynamic base excitation. The pore fluid flow was idealized using averaged Navier–Stokes equations and the discrete element method was employed to model the solid phase particles. A well established semi‐empirical relationship was utilized to quantify the fluid–particle interactions. The conducted simulations revealed a number of salient micro‐mechanical mechanisms and response patterns associated with the deposit liquefaction. Space and time variation of porosity was a major factor which affected the coupled response of the solid and fluid phases. Pore fluid flow was within Darcy's regime. The predicted response exhibited macroscopic patterns consistent with experimental results and case histories of the liquefaction of granular soil deposits. Copyright © 2004 John Wiley & Sons, Ltd.     

Metamorphic phase relations in orthopyroxene-bearing granitoids: implication for high-pressure metamorphism and prograde melting in the continental crust

In this work, the factors controlling the formation and preservation of high-pressure mineral assemblages in the metamorphosed orthopyroxene-bearing metagranitoids of the Sandmata Complex, Aravalli-Delhi Mobile Belt (ADMB), northwestern India have been modelled. The rocks range in composition from farsundite through quartz mangerite to opdalite, and with varying K₂O, Ca/(Ca + Na)_rock and FeO_tot + MgO contents. A two stage metamorphic evolution has been recorded in these rocks.
An early hydration event stabilized biotite with or without epidote at the expense of magmatic orthopyroxene and plagioclase. Subsequent high-pressure granulite facies metamorphism (∼15 kbar, ∼800 °C) of these hydrated rocks produced two rock types with contrasting mineralogy and textures. In the non-migmatitic metagranitoids, spectacular garnet ± K-feldspar ± quartz corona was formed around reacting biotite, plagioclase, quartz and/or pyroxene. In contrast, biotite ± epidote melting produced migmatites, containing porphyroblastic garnet incongruent solids and leucosomes.
Applying NCKFMASHTO T–M (H₂O) and P–T pseudosection modelling techniques, it is demonstrated that the differential response of these magmatic rocks to high-pressure metamorphism is primarily controlled by the scale of initial hydration. Rocks, which were pervasively hydrated, produced garnetiferous migmatites, while for limited hydration, the same metamorphism formed sub-solidus garnet-bearing coronae. Based on the sequence of mineral assemblage evolution and the mineral compositional zoning features in the two metagranitoids, a clockwise metamorphic P–T path is constrained for the high-pressure metamorphic event. The finding has major implications in formulating geodynamic model of crustal amalgamation in the ADMB.     

Long-Period Variations of the Eccentricity Vector Valid also for Near Circular Orbits around a Non-Spherical Body

This paper studies the long period variations of the eccentricity vector of the orbit of an artificial satellite, under the influence of the gravity field of a central body. We use modified orbital elements which are non-singular at zero eccentricity. We expand the long periodic part of the corresponding Lagrange equations as power series of the eccentricity. The coefficients characterizing the differential system depend on the zonal coefficients of the geopotential, and on initial semi-major axis, inclination, and eccentricity. The differential equations for the components of the eccentricity vector are then integrated analytically, with a definition of the period of the perigee based on the notion of “free eccentricity”, and which is also valid for circular orbits. The analytical solution is compared to a numerical integration. This study is a generalization of (Cook, Planet. Space Sci., 14, 1966): first, the coefficients involved in the differential equations depend on all zonal coefficients (and not only on the very first ones); second, our method applies to nearly circular orbits as well as to not too eccentric orbits. Except for the critical inclination, our solution is valid for all kinds of long period motions of the perigee, i.e., circulations or librations around an equilibrium point.     

The tidal stripping of satellites

We present an improved analytic calculation for the tidal radius of satellites and test our results against N -body simulations.
The tidal radius in general depends upon four factors: the potential of the host galaxy, the potential of the satellite, the orbit of the satellite and the orbit of the star within the satellite . We demonstrate that this last point is critical and suggest using three tidal radii to cover the range of orbits of stars within the satellite. In this way we show explicitly that prograde star orbits will be more easily stripped than radial orbits; while radial orbits are more easily stripped than retrograde ones. This result has previously been established by several authors numerically, but can now be understood analytically. For point mass, power-law (which includes the isothermal sphere), and a restricted class of split power-law potentials our solution is fully analytic. For more general potentials, we provide an equation which may be rapidly solved numerically.
Over short times (≲1–2 Gyr ∼1 satellite orbit), we find excellent agreement between our analytic and numerical models. Over longer times, star orbits within the satellite are transformed by the tidal field of the host galaxy. In a Hubble time, this causes a convergence of the three limiting tidal radii towards the prograde stripping radius. Beyond the prograde stripping radius, the velocity dispersion will be tangentially anisotropic.     

Elastic solutions for stresses in a transversely isotropic half‐space subjected to three‐dimensional buried parabolic rectangular loads

In many areas of engineering practice, applied loads are not uniformly distributed but often concentrated towards the centre of a foundation. Thus, loads are more realistically depicted as distributed as linearly varying or as parabola of revolution. Solutions for stresses in a transversely isotropic half‐space caused by concave and convex parabolic loads that act on a rectangle have not been derived. This work proposes analytical solutions for stresses in a transversely isotropic half‐space, induced by three‐dimensional, buried, linearly varying/uniform/parabolic rectangular loads. Load types include an upwardly and a downwardly linearly varying load, a uniform load, a concave and a convex parabolic load, all distributed over a rectangular area. These solutions are obtained by integrating the point load solutions in a Cartesian co‐ordinate system for a transversely isotropic half‐space. The buried depth, the dimensions of the loaded area, the type and degree of material anisotropy and the loading type for transversely isotropic half‐spaces influence the proposed solutions. An illustrative example is presented to elucidate the effect of the dimensions of the loaded area, the type and degree of rock anisotropy, and the type of loading on the vertical stress in the isotropic/transversely isotropic rocks subjected to a linearly varying/uniform/parabolic rectangular load. Copyright © 2002 John Wiley & Sons, Ltd.     

Progress in Earthquake Science and Technology in China: Review and Prospects (Ⅰ)

In the article the author looks back the hard development course and great progress in earth quake science and technology in China during the last near a half of century and expounds the following 3 aspects: (1) The strong desire of the whole society to mitigate seismic disasters and reduce the effect of earthquakes on social-economic live is a great driving force to push forward the development of earthquake science and technology in China; (2) To better ensure people‘ s life and property, sustainable economic development, and social stability is an essential purpose to drive the development of earthquake science and technology in China; and (3) To insist on the dialectical connection of setup of technical system for seismic monitoring with the scientific research of earthquakes and to better handle the relation between crucial task, current scientif ic level, and the feasibility are the important principles to advance the earthquake science and technology in China. Some success and many setbacks in earthquake disaster mitigation consistently enrich our knowledge regarding the complexity of the conditions for earthquake occurrence and the process of earthquake preparation, promote the reconstruction and modernization of technical system for earthquake monitoring, and deepen the scientific research of earthquakes. During the last 5 years, the improvement and modernization of technical system for earthquake monitoring have clearly provided the technical support to study and practice of earthquake prediction and pre caution, give prominence to key problems and broaden the field of scientific research of earth quakes. These have enabled us to get some new recognition of the conditions for earthquake oc currence and process of earthquake preparation, characteristics of seismic disaster, and mecha nism for earthquake generation in China‘s continent. The progress we have made not only en courages us to enhance the effectiveness of earthquake disaster mitigation, but also provides a basis for accelerating further development of earthquake science and technology in China in the new century, especially in the 10th five-year plan. Based on the history reviewed, the author sets forth a general requirement for develop ment of earthquake science and technology in China and brings out 10 aspects to be stressed and strengthened at present and in the future. These are: upgrade and setup of the network of digitized seismic observation; upgrade and setup of the network for observation of seismic pre cursors; setup of the network for observation of strong motion; setup of the laboratories for ex periment on seismic regime; establishment of technical system for seismic information, emer gency command and urgent rescue; research on short-term and imminent earthquake predic tion; research on intermediate- and long-term earthquake prediction; research on attenuation of seismic ground motion, mechanism for seismic disaster, and control on seismic disaster; ba sic research fields related to seismology and geoscience. We expect that these efforts will signifi cantly elevate the level of earthquake science and technology in China to the advanced interna tional level, improve theories, techniques, and methods for earthquake precaution and predic tion, and enhance the effectiveness of earthquake disaster mitigation.     

Weathering process at the natural fission reactor of Bangombé

 The uranium deposits in the basin of Franceville (Gabon) host the only natural fission reactors known in the world. Unique geological conditions favoured a natural fission reaction 2 Ga ago. This was detected by anomalous isotopic compositions of uranium and rare earth elements (REE), which are produced by the fission reaction. In total, 16 reactor zones were found. Most of them are mined out. The reactor zone of Bangombé, is only 10–11 m below the surface. This site has been influenced by surface weathering processes. Six drill cores have been sampled at the site of the reactor zone of Bangombé during the course of the study and only one drill core (BAX 08) hit the core of the reactor. From these data and previous drilling campaigns, the reactor size is estimated to be 10 cm thick, 2–3 m wide and 4–6 m long. The migration of fission products can be traced by the anomalous isotope ratios of REE because of the fission process. The ¹⁴⁹Sm/¹⁴⁷Sm ratio close to the reactor zone is only 0.28 (normal: 0.92) because of the intense neutron capture of ¹⁴⁹Sm and subsequent transmutation, whereas ¹⁴⁷Sm is enriched by the fission reaction. Similar changes in isotopic patterns are detectable on other REE. The isotope ratios of Sm and Nd of whole rock and fracture samples surrounding the reactor indicate that fission-genic REE migrated only a few decimetres above and mainly below the reactor zone. Organic matter (bitumen) seems to act as a trap for fission-genic REE. Additional REE-patterns show less intense weathering with increasing depth in the log profile and support a simple weathering model. Received: 26 November 1999 · Accepted: 2 May 2000