A transmission electron microscope investigation of shock metamorphism in olivine of the Ilafegh 013 chondrite

Abstract— Ilafegh 013 is a slightly weathered, type-3 ordinary chondrite that is unusual in that it is highly shocked. Olivine grains in chondrules were studied in detail by transmission electron microscopy (TEM) to characterize the dislocations, fractures and recrystallization textures. Only c dislocations, in the (010) glide plane, seem to be activated by the shock wave. Fractures are more numerous in the TEM than at the scale of the optical microscope. They can be divided into two groups: (a) Open fractures with injections of molten metal, sulphide, and in some cases, phases containing Ca, Al, Na and K. (b) Very straight and narrow cracks parallel to {130} planes without injected material and consisting of imperfectly relaxed subgrain boundaries composed of two families of dislocations (Burgers vectors a and c ). Both types of dislocations are roughly in screw orientation (i.e., the subgrain boundaries have a pronounced twist character). It is suggested that they formed during the shock in order to relax stresses around shear fractures (fractures with displacements of their walls parallel to the fracture plane). Recrystallization occurs in areas where the dislocation density is high. Since diaplectic glass was not detected, it is suggested that recrystallization resulted from the large elastic energy of the lattice defects stored in the deformed crystals.     

A model for conductive faults with non-matching grids

In this paper, we are interested in modeling single-phase flow in a porous medium with known faults seen as interfaces. We mainly focus on how to handle non-matching grids problems arising from rock displacement along the fault. We describe a model that can be extended to multi-phase flow where faults are treated as interfaces. The model is validated in an academic framework and is then extended to 3D non K-orthogonal grids, and a realistic case is presented.     

Rockfall rebound: comparison of detailed field experiments and alternative modelling approaches

The accuracy of rockfall trajectory simulations mainly rests on the calculation of the rebound of fragments following their impact on the slope. This paper is dedicated to the comparative analysis of two rebound modelling approaches currently used in rockfall simulation using field experiments of single rebounds. The two approaches consist in either modelling the rock as a single material point (lumped mass approach) or in explicitly accounting for the fragment shape (rigid body approach). A lumped mass model accounting for the coupling between translational and rotational velocities and introducing a slope perturbation angle was used. A rigid body approach modelling the rocks as rigid locally deformable (in the vicinity of the contact surface) assemblies of spheres was chosen. The comparative analysis of the rebound models shows that both of them are efficient with only a few parameters. The main limitation of each approach are the calibration of the value of the slope perturbation (‘roughness’) angle, for the lumped mass approach, and the estimation of the rock length and height from field geological and historical analyses, for the rigid body approach. Finally, both rebound models require being improved in a pragmatic manner to better predict the rotational velocities distribution. Copyright © 2012 John Wiley & Sons, Ltd.     

On regularized time varying gravity field models based on grace data and their comparison with hydrological models

Determination of spherical harmonic coefficients of the Earth’s gravity field is often an ill-posed problem and leads to solving an ill-conditioned system of equations. Inversion of such a system is critical, as small errors of data will yield large variations in the result. Regularization is a method to solve such an unstable system of equations. In this study, direct methods of Tikhonov, truncated and damped singular value decomposition and iterative methods of ν, algebraic reconstruction technique, range restricted generalized minimum residual and conjugate gradient are used to solve the normal equations constructed based on range rate data of the gravity field and climate experiment (GRACE) for specific periods. Numerical studies show that the Tikhonov regularization and damped singular value decomposition methods for which the regularization parameter is estimated using quasioptimal criterion deliver the smoothest solutions. Each regularized solution is compared to the global land data assimilation system (GLDAS) hydrological model. The Tikhonov regularization with L-curve delivers a solution with high correlation with this model and a relatively small standard deviation over oceans. Among iterative methods, conjugate gradient is the most suited one for the same reasons and it has the shortest computation time.     

Possible long-term decline in impact rates: 1. Martian geological data

Quantin et al. [Quantin, C., Allemand, P., Mangold, N., Delacourt, C., 2004a. Icarus 172, 555-572] tabulated crater count data for 56 landslides along the walls of Valles Marineris. Under the assumption of a constant cratering rate after about 3 Gyr ago, as used in the 1999-2005 iterations of the crater chronology isochron system of Hartmann, and in the Hartmann and Neukum system, these data indicate a regularly increasing rate of landslides, which would be difficult to explain. We suggest that these data may support a decline in inner Solar System cratering rates by about a factor of 3 since 3 Gyr ago, not unlike predictions based on asteroid belt collision models. Such a decline is also supported by our review of data on lunar impact melts and glass spherules in a companion paper [Hartmann, W.K., Quantin, C., Mangold, N., 2007. Icarus 186, 11-23]. Such models produce not only a more uniform rate of landslides over the last 3 Gyr, but also a more uniform rate of resurfacing processes which also had an apparent increase under the assumption of a constant cratering rate.     

Power spectrum for the small-scale Universe

The first objects to arise in a cold dark matter (CDM) universe present a daunting challenge for models of structure formation. In the ultra small-scale limit, CDM structures form nearly simultaneously across a wide range of scales. Hierarchical clustering no longer provides a guiding principle for theoretical analyses and the computation time required to carry out credible simulations becomes prohibitively high. To gain insight into this problem, we perform high-resolution  ( N = 720³–1584³)  simulations of an Einstein–de Sitter cosmology where the initial power spectrum is   P ( k ) ∝ k ⁿ ,  with  −2.5 ≤ n ≤− 1  . Self-similar scaling is established for   n =−1  and −2 more convincingly than in previous, lower resolution simulations and for the first time, self-similar scaling is established for an   n =−2.25  simulation. However, finite box-size effects induce departures from self-similar scaling in our   n =−2.5  simulation. We compare our results with the predictions for the power spectrum from (one-loop) perturbation theory and demonstrate that the renormalization group approach suggested by McDonald improves perturbation theory's ability to predict the power spectrum in the quasi-linear regime. In the non-linear regime, our power spectra differ significantly from the widely used fitting formulae of Peacock & Dodds and Smith et al. and a new fitting formula is presented. Implications of our results for the stable clustering hypothesis versus halo model debate are discussed. Our power spectra are inconsistent with predictions of the stable clustering hypothesis in the high- k limit and lend credence to the halo model. Nevertheless, the fitting formula advocated in this paper is purely empirical and not derived from a specific formulation of the halo model.     

Adaptive mapping functions to the azimuthal anisotropy of the neutral atmosphere

The anisotropy of propagation of radio waves used by global navigation satellite systems is investigated using high-resolution observational data assimilations produced by the European Centre for Medium-range Weather Forecast. The geometry and the refractivity of the neutral atmosphere are built introducing accurate geodetic heights and continuous formulations of the refractivity and its gradient. Hence the realistic ellipsoidal shape of the refractivity field above the topography is properly represented. Atmospheric delays are obtained by ray-tracing through the refractivity field, integrating the eikonal differential system. Ray-traced delays reveal the anisotropy of the atmosphere. With the aim to preserve the classical mapping function strategy, mapping functions can evolve to adapt to high-frequency atmospheric fluctuations and to account for the anisotropy of propagation by fitting at each site and time the zenith delays and the mapping functions coefficients. Adaptive mapping functions (AMF) are designed with coefficients of the continued fraction form which depend on azimuth. The basic idea is to expand the azimuthal dependency of the coefficients in Fourier series introducing a multi-scale azimuthal decomposition which slightly changes the elevation functions with the azimuth. AMF are used to approximate thousands of atmospheric ray-traced delays using a few tens of coefficients. Generic recursive definitions of the AMF and their partial derivatives lead to observe that the truncation of the continued fraction form at the third term and the truncation of the azimuthal Fourier series at the fourth term are sufficient in usual meteorological conditions. Delays’ and elevations’ mapping functions allow to store and to retrieve the ray-tracing results to solve the parallax problem at the observation level. AMF are suitable to fit the time-variable isotropic and anisotropic parts of the ray-traced delays at each site at each time step and to provide GPS range corrections at the measurement level with millimeter accuracy at low elevation. AMF to the azimuthal anisotropy of the neutral atmosphere are designed to adapt to complex weather conditions by adaptively changing their truncations.