Constraints in 3D gravity inversion

A three-dimensional (3D) inversion program is developed to interpret gravity data using a selection of constraints. This selection includes minimum distance, flatness, smoothness and compactness constraints, which can be combined using a Lagrangian formulation. A multigrid technique is also implemented to resolve separately large and short gravity wavelengths. The subsurface in the survey area is divided into rectangular prismatic blocks and the problem is solved by calculating the model parameters, i.e. the densities of each block. Weights are given to each block depending on depth, a priori information on density and the density range allowed for the region under investigation. The present computer code is tested on modelled data for a dipping dike and multiple bodies. Results combining different constraints and a weight depending on depth are shown for the dipping dike. The advantages and behaviour of each method are compared in the 3D reconstruction. Recovery of geometry (depth, size) and density distribution of the original model is dependent on the set of constraints used. From experimentation, the best combination of constraints for multiple bodies seems to be flatness and a minimum volume for the multiple bodies. The inversion method is tested on real gravity data from the Rouyn-Noranda (Quebec) mining camp. The 3D inversion model for the first 10 km is in agreement with the known major lithological contacts at the surface; it enables the determination of the geometry of plutons and intrusive rocks at depth.     

Grain-size-sensitive flow and shear-stress enhancement at the brittle–ductile transition of the continental crust

Localized shear zones along low-angle normal faults have been identified in regions of extension at the brittle-ductile transition of the continental crust. The possibility of the strain localizing at a depth of 10 km is interpreted here as a consequence of an increase in the equivalent shear stress applied to the flow of the lower crust. This enhancement of the flow stress is seen as a prerequisite for the triggering of brittle deformation mechanisms leading to strain localization. The lower crust rheology used to examine this stress increase is strain-rate, temperature and grain-size dependent, due to the coupling of dislocation and diffusion creep. The model structure proposed consists of a top layer, the upper crust, gliding rigidly above a bottom layer, the lower crust, which deforms in simple shear. During a short time interval (1400 years), the equivalent shear stress is found to increase by a factor of up to 3 (67 MPa for anorthite and 17 MPa for quartz). For anorthite, this stress could explain the activation of a Mohr-Coulomb failure with a friction coefficient of 0.2, which is reasonable at the depth of 10 km. Dislocation creep is activated during a rapid change in the prescribed velocity, whereas diffusion creep dominates if the velocity is held constant, highlighting the importance of grain-size sensitivity for lower crustal rheology.     

Early Pleistocene climate changes in the central Mediterranean region as inferred from integrated pollen and planktonic foraminiferal stable isotope analyses

Vegetation inherited from a Pliocene subtropical climate evolved through obliquity oscillations and global cooling leading to modern conditions. An integrated, highly time-resolved record of pollen and stable isotopes (δ¹⁸O and δ¹³C of Globigerina bulloides) was obtained to understand vegetation responses to Early Pleistocene climate changes. Continental and marine responses are compared in the Central Mediterranean region with a particular consideration of environmental changes during anoxic events.Pollen data illustrate vegetation dynamics as follows: [1] development of mesothermic elements (warm and humid conditions); [2] expansion of mid- and high-altitude elements (cooler but still humid conditions); and [3] strengthening of steppe and herb elements (cooler and dry conditions). These successions correlate with precession. δ¹⁸O variations recorded by Globigerina bulloides define two cycles (MIS 43-40) related to obliquity. At northern low- to mid-latitudes, the pollen signal records temperature and wetness changes related to precession even during global climate changes induced by obliquity. This may result in unexpected increasing wetness during glacial periods, which has to be considered specific to the Central and Eastern Mediterranean region. Lastly, an analysis of anoxic events reveals that enhanced runoff is indicated by increasing frequency of the riparian trees Liquidambar and Zelkova.     

Precipitation water storage capacity in a temperate mixed oak–beech canopy

The storage capacity of a temperate mixed oak–beech stand was investigated as a function of stand density and species composition. Measurements were performed in selected zones delimited by three neighbouring trees. Three independent approaches were compared: (i) a spraying laboratory experiment to estimate the water storage on foliage before and after dripping; (ii) a mechanistic model describing rainfall partitioning within the forest canopy and providing estimates of foliage storage capacities; and (iii) linear regression analyses to evaluate the canopy (foliage + branches) storage capacity using the relationship between throughfall and rainfall. Good agreement was generally observed between the laboratory experiment and the mechanistic model estimates, while estimations from the regression method tended to exceed those from the other approaches. Storage capacity estimates ranged from 0·22 mm to 0·80 mm for pure oak zones, from 0·24 mm to 1·12 mm for mixed zones and from 0·53 mm to 1·17 mm for pure beech zones. The increase of storage capacity with increasing proportion of beech in the canopy resulted from higher beech LAI compared with oak. Similarly, for mixed and pure beech canopies, storage capacity was higher for high density zones than for low density zones as a result of the increase in LAI with increasing local basal area; in contrast, for pure oak, the storage capacity was not related to basal area because of the lower shade‐tolerance of this species compared with beech. Copyright © 2008 John Wiley & Sons, Ltd.     

The Monitor project: JW 380 – a 0.26-, 0.15-M⊙, pre-main-sequence eclipsing binary in the Orion nebula cluster

We present a general recipe for constructing N -body realizations of galaxies comprising near spherical and disc components. First, an exact spherical distribution function for the spheroids (halo and bulge) is determined, such that it is in equilibrium with the gravitational monopole of the disc components. Second, an N -body realization of this model is adapted to the full disc potential by growing the latter adiabatically from its monopole. Finally, the disc is sampled with particles drawn from an appropriate distribution function, avoiding local-Maxwellian approximations. We performed test simulations and find that the halo and bulge radial density profile very closely match their target model, while they become slightly oblate due to the added disc gravity. Our findings suggest that vertical thickening of the initially thin disc is caused predominantly by spiral and bar instabilities, which also result in a radial re-distribution of matter, rather than scattering off interloping massive halo particles.