Lunar motion: Theory and observations

Results of several fits of the lunar theory ELP 2000-82B and of Moons' theory of libration are presented. The theories are fitted both to JPL numerical integrations and to LLR observations     

Consolidation of double‐layered ground with vertical drains

Vertical drains are usually installed in subsoil consisting of several layers. Due to the complex nature of the problem, over the past decades, the consolidation properties of multi‐layered ground with vertical drains have been analysed mainly by numerical methods. An analytical solution for consolidation of double‐layered ground with vertical drains under quasi‐equal strain condition is presented in this paper. The main steps for the computation procedure are listed. The convergence of the series solution is discussed. The comparisons between the results obtained by the present analytical method and the existing numerical solutions are described by figures. The orthogonal relation for the system of double‐layered ground with vertical drains is proven. Finally, some consolidation properties of double‐layered ground with vertical drains are analysed. Copyright © 2001 John Wiley & Sons, Ltd.     

Three dimensional simulation of fluid flow in X‐ray CT images of porous media

A numerical scheme is developed in order to simulate fluid flow in three dimensional (3‐D) microstructures. The governing equations for steady incompressible flow are solved using the semi‐implicit method for pressure‐linked equations (SIMPLE) finite difference scheme within a non‐staggered grid system that represents the 3‐D microstructure. This system allows solving the governing equations using only one computational cell. The numerical scheme is verified through simulating fluid flow in idealized 3‐D microstructures with known closed form solutions for permeability. The numerical factors affecting the solution in terms of convergence and accuracy are also discussed. These factors include the resolution of the analysed microstructure and the truncation criterion. Fluid flow in 2‐D X‐ray computed tomography (CT) images of real porous media microstructure is also simulated using this numerical model. These real microstructures include field cores of asphalt mixes, laboratory linear kneading compactor (LKC) specimens, and laboratory Superpave gyratory compactor (SGC) specimens. The numerical results for the permeability of the real microstructures are compared with the results from closed form solutions. Copyright © 2004 John Wiley & Sons, Ltd.     

Relict mountain slope deposits of northern Portugal: facies,sedimentogenesis and environmental implications

This paper deals with the formation processes and the palaeoenvironmental significance of relict slope deposits located on the uppermost part of the north Portugal mountains. For this purpose, seven key sites representative of the different lithofacies have been selected and analysed in detail. The data show that three main dynamic processes are responsible for the emplacement of regional fossil slope deposits: runoff, debris flows and dry grain flows. The ubiquity of these processes and the lack of frost‐related features or landforms do not support the existence of severe Pleistocene climates in this part of the lberian Peninsula as postulated by previous work. Pedological data gathered at one of the study sites show that a subalpine environment was probably present at 700–800 m altitude between 29 and 14 kyr. Using data from the Pyrenees Mountains, a 6.5 to 12°C depression in mean annual temperature has been tentatively postulated for this Pleniglacial period. Copyright © 2003 John Wiley & Sons, Ltd.     

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.     

Effect of spatial variation characteristics on contouring of design storm depth

The ordinary kriging method, a geostatistical interpolation technique, was applied for developing contour maps of design storm depth in northern Taiwan using intensity–duration–frequency (IDF) data. Results of variogram modelling on design storm depths indicate that the design storms can be categorized into two distinct storm types: (i) storms of short duration and high spatial variation and (ii) storms of long duration and less spatial variation. For storms of the first category, the influence range of rainfall depth decreases when the recurrence interval increases, owing to the increasing degree of their spatial independence. However, for storms of the second category, the influence range of rainfall depth does not change significantly and has an average of approximately 72 km. For very extreme events, such as events of short duration and long recurrence interval, we do not recommend usage of the established design storm contours, because most of the interstation distances exceed the influence ranges. Our study concludes that the influence range of the design storm depth is dependent on the design duration and recurrence interval and is a key factor in developing design storm contours. Copyright © 2003 John Wiley & Sons, Ltd.     

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.     

The influence of a non‐associated flow rule on the calculation of the factor of safety of soil slopes

This paper presents a method that incorporates a non‐associated flow rule into the limit analysis to investigate the influence of the dilatancy angle on the factor of safety for the slope stability analysis. The proposed method retain's the advantage of the upper bound method, which is simple and has no stress involvement in the calculation of the energy dissipation and the factor of safety. Copyright © 2001 John Wiley & Sons, Ltd.     

Abundance and isotopic composition of noble gases in metal and graphite of the Bohumilitz IAB iron meteorite

Abstract— Abundances and isotopic compositions of noble gases in metal and graphite of the Bohumilitz IAB iron meteorite were measured. The abundance ratios of spallogenic components in metal reveal a ³He deficiency which is due to the diffusive loss of parent isotopes, that is, tritium (Tilles, 1963; Schultz, 1967). The diffusive loss likely has been induced by thermal heating by the Sun during cosmic‐ray exposure (～160 Ma; Lavielle et al, 1999). Thermal process such as impact‐induced partial loss may have affected the isotopic composition of spallogenic Ne. The ¹²⁹Xe/¹³¹Xe ratio of cosmogenic components in the metal indicates an enhanced production of epi‐thermal neutrons. The abundance ratios of spallogenic components in the graphite reveal that it contained small amounts of metal and silicates. The isotopic composition of heavy noble gases in graphite itself was obtained from graphite treated with HF/HCl. The isotopic composition of the etched graphite shows that it contains two types of primordial Xe (i.e., Q‐Xe and El Taco Xe). The isotopic heterogeneity preserved in the Bohumilitz graphite indicates that the Bohumilitz graphite did not experience any high‐temperature event and, consequently, must have been emplaced into the metal at subsolidus temperatures. This situation is incompatible with an igneous model as well as the impact melting models for the IAB‐IIICD iron meteorites as proposed by Choi et al. (1995) and Wasson et al (1980).     

Low‐energy helium ion irradiation‐induced amorphization and chemical changes in olivine: Insights for silicate dust evolution in the interstellar medium

Abstract— We present the results of irradiation experiments aimed at understanding the structural and chemical evolution of silicate grains in the interstellar medium. A series of He+ irradiation experiments have been performed on ultra‐thin olivine, (Mg,Fe)₂SiO₄, samples having a high surface/volume (S/V) ratio, comparable to the expected S/V ratio of interstellar dust. The energies and fluences of the helium ions used in this study have been chosen to simulate the irradiation of interstellar dust grains in supernovae shock waves. The samples were mainly studied using analytical transmission electron microscopy. Our results show that olivine is amorphized by low‐energy ion irradiation. Changes in composition are also observed. In particular, irradiation leads to a decrease of the atomic ratios O/Si and Mg/Si as determined by x‐ray photoelectron spectroscopy and by x‐ray energy dispersive spectroscopy. This chemical evolution is due to the differential sputtering of atoms near the surfaces. We also observe a reduction process resulting in the formation of metallic iron. The use of very thin samples emphasizes the role of surface/volume ratio and thus the importance of the particle size in the irradiation‐induced effects. These results allow us to account qualitatively for the observed properties of interstellar grains in different environments, that is, at different stages of their evolution: chemical and structural evolution in the interstellar medium, from olivine to pyroxene‐type and from crystalline to amorphous silicates, porosity of cometary grains as well as the formation of metallic inclusions in silicates.