A new multiscale scheme for computing statistical moments in single phase flow in heterogeneous porous media

In this work we develop a new multiscale procedure to compute numerically the statistical moments of the stochastic variables which govern single phase flow in heterogeneous porous media. The technique explores the properties of the log-normally distributed hydraulic conductivity, characterized by power-law or exponential covariances, which shows invariance in its statistical structure upon a simultaneous change of the scale of observation and strength of heterogeneity. We construct a family of equivalent stochastic hydrodynamic variables satisfying the same flow equations at different scales and strengths of heterogeneity or correlation lengths. Within the new procedure the governing equations are solved in a scaled geology and the numerical results are mapped onto the original medium at coarser scales by a straightforward rescaling. The new procedure is implemented numerically within the Monte Carlo algorithm and also in conjunction with the discretization of the low-order effective equations derived from perturbation analysis. Numerical results obtained by the finite element method show the accuracy of the new procedure to approximated the two first moments of the pressure and velocity along with its potential in reducing drastically the computational cost involved in the numerical modeling of both power-law and exponential covariance functions.     

The story of galaxy evolution in full colour

A Holy Grail of modern astronomy is understanding the origin of Edwin Hubble's morphological sequence of galaxy types. What made some collapsing gas clouds turn into elegant spiral systems like our own Milky Way, whereas others became smooth, featureless ellipticals? More fundamentally, does the taxonomic scheme introduced by Hubble in the 1920s have any physical relevance? The Space Telescope that bears Hubble's name is providing answers to these and other questions in the context of modern theories of structure formation.
Richard S Ellis, Roberto G Abraham, Jarle Brinchmann and Felipe Menanteau show how modern telescopes producing colourful images are pointing the way to understanding galaxy evolution.     

Reconstructing fluvial bar surfaces from compound cross‐strata and the interpretation of bar accretion direction in large river deposits

The interpretation of fluvial styles from the rock record is based for a significant part on the identification of different types of fluvial bars, characterized by the geometric relationship between structures indicative of palaeocurrent and surfaces interpreted as indicative of bar form and bar accretion direction. These surfaces of bar accretion are the boundaries of flood‐related bar increment elements, which are typically less abundant in outcrops than what would be desirable, particularly in large river deposits in which each flood mobilizes large volumes of sediment, causing flood‐increment boundary surfaces to be widely spaced. Cross‐strata set boundaries, on the other hand, are abundant and indirectly reflect the process of unit bar accretion, inclined due to the combined effect of the unit bar surface inclination and the individual bedform climbing angle, in turn controlled by changes in flow structure caused by local bar‐scale morphology. This work presents a new method to deduce the geometry of unit bar surfaces from measured pairs of cross‐strata and cross‐strata set boundaries. The method can be used in the absence of abundant flood‐increment bounding surfaces; the study of real cases shows that, for both downstream and laterally accreting bars, the reconstructed planes are very similar to measured bar increment surfaces.     

Quantum Instability of Magnetized Stellar Objects

The equations of state for degenerate electron and neutron gases are studied in the presence of magnetic fields. After including quantum effects in the investigation of the structural properties of these systems, it is found that some hypermagnetized stars can be unstable according to the criterion of stability of pressures. Highly magnetized white dwarfs should collapse producing a supernova type Ia, while superstrong magnetized neutron stars cannot stand their own magnetic field and must implode, too. A comparison of our results with a set of the available observational data of some compact stars is also presented, and the agreement between this theory and observations is verified.     

Degradation of rope properties under increasing monotonic load

This paper describes the development of a computational model to predict the response of synthetic-fiber ropes under both monotonic and cyclic loads. The model addresses the issue of damage to a rope cross-section and its effect on rope response. Isotropic damage is assumed and is therefore characterized by a scalar quantity known as the damage index. The damage index is used to quantify the amount of deterioration that takes place in a damaged rope throughout its loading history. For the current study, data obtained from static tension tests are used to determine the evolution relationship for the damage index. This relationship is based on the stiffness ratio between experimental data and simulation of rope response without any source of damage and is assumed to be described by an exponential equation. Numerical simulations of damaged rope behavior compare well with experimental results.     

Factors controlling the vegetation distribution and coal‐forming environments in a strike‐slip basin. The Pennsylvanian Peñarroya‐Belmez‐Espiel Basin,southern Spain

Coal‐forming environments require humid to perhumid conditions. Tectonics governs the size, location and availability of coal seams developed in such environments. While large Pennsylvanian paralic basins generated thick and continuous coal seams, many other small coeval basins, which were tectonically active, developed a puzzling succession, with carbonaceous deposits that varied in size, thickness and the nature of the coal‐forming flora. This study, conducted in the Peñarroya‐Belmez‐Espiel coalfield, a Variscan strike‐slip basin in the south of Spain, provides insights into this subject. The coal seams analysed, generated in different depositional environments, have quantitatively different palynological assemblages. Lacustrine coals are dominated by lycopsids; distal alluvial plain/marginal lacustrine coals are dominated by sphenophytes and tree ferns, and middle alluvial fan coals are dominated by sphenophytes, tree ferns and lycopsids. This means that when conditions were favourable for peat accumulation, peat accumulated regardless of the nature of the available flora.     

Geology and paleomagnetism of El Potrero pluton, Baja California: Understanding criteria for timing of deformation and evidence of pluton tilt during batholith growth

The 102 Ma El Potrero pluton, in the western foothills of Sierra San Pedro Mártir, in north-central Baja California, was emplaced during a long period of contractional deformation bracketed between 132 and 85 Ma that affected this segment of the Peninsular Ranges Batholith. The pluton records regional and emplacement related deformation manifested by: (1) a solid-state fabric developed on its eastern contact, which is produced by eastward lateral pluton expansion; (2) cleavage triple point zones in the host-rock NW and SE of the pluton; (3) subhorizontal ductile shear zones indicative of top-to-the-east transport; (4) magmatic and tectonic foliations parallel to regional structural trends and regional shear zones; (5) variable axial ratios of microgranitoid enclaves close to pluton–wall rock contacts; (6) evidence of brittle-emplacement mechanisms in the western border of the pluton, which contrast with features indicating mainly ductile mechanisms toward the east; and, (7) markedly discordant paleomagnetic directions that suggest emplacement in an active tectonic setting. The overall mean for 9 accepted paleomagnetic sites is Dec = 34.6°, I = 25.7° (k = 88.3, α₉₅ = 5.5°), and is deviated  35° with respect to the reference cratonic direction. This magnetization is interpreted to indicate a combination of tilt due to initial drag during vertical diapiric ascent (or westward lateral-oblique expansion) of the adjacent San Pedro Mártir pluton and later rotation ( 15°) by Rosarito Fault activity in the southwest; this rotation may have occurred as eastward contraction acted to fill the space emptied by the ascending San Pedro Mártir pluton. The Rosarito fault may have tilted several plutons in the area (Sierra San Pedro Mártir, El Potrero, San José, and Encinosa). Magnetic susceptibility fabrics for 13 sites reflect mostly emplacement-related stress and regional stress. Paleomagnetic data and structural observations lead us to interpret the El Potrero pluton as a syntectonic pluton, emplaced within a regional shear zone delimited by the Main Mártir Thrust and the younger Rosarito Fault.     

Ray-traced slant factors for mitigating the tropospheric delay at the observation level

Three-dimensional ray tracing through a numerical weather model has been applied to a global precise point positioning (PPP) campaign for modeling both the elevation angle- and azimuth-dependence of the tropospheric delay. Rather than applying the ray-traced slant delays directly, the delay has been parameterized in terms of slant factors, which are applied in a similar manner to traditional mapping functions, but which can account for the azimuthal asymmetry of the delay. Five strategies are considered: (1) Vienna Mapping Functions 1 (VMF1) and estimation of a residual zenith delay parameter; (2) VMF1, estimation of a residual zenith delay and estimation of two tropospheric gradient parameters; (3) three-dimensional ray-traced slant factors and estimation of a residual zenith delay; (4) using only ray-traced slant factors and no estimation of any tropospheric parameters and; (5) using both ray-traced slant factors and estimating a residual zenith delay and two tropospheric gradient parameters. The use of the ray-traced slant factors (solution 3) showed a 3.8% improvement in the repeatability of the up component when compared to the assumption of a symmetric atmosphere (solution 1), while the estimation of two tropospheric gradient parameters gave the best results showing an 7.6% improvement over solution 1 in the up component. Solution 4 performed well in the horizontal domain, allowing for sub-centimeter repeatability but the up component was degraded due to deficiencies in the modeling of the zenith delay, particularly for stations located at equatorial latitudes. The magnitude of the differences in the mean coordinates between solution 2 and solution 3, and the strong correlation with the differences between the north component and the ray-traced gradients (coefficient of correlation of 0.83), as well as the impact of observation geometry on the gradient solution indicate that the use of the ray-traced slant factors could have an implication on the realization of reference frames. The estimated tropospheric products from the PPP solutions were compared to those derived from ray tracing. For the zenith delay, a root mean square (RMS) of 5.4 mm was found, while for the gradient terms, a correlation coefficient of 0.46 for the N–S and 0.42 for the E–W was found for the north–south and east–west components, suggesting that there are still important differences in the gradient parameters which could be due to either errors in the NWM or to non-tropospheric error sources leaking into the PPP-estimated gradients.     

Identifying Controls on Water Chemistry of Tropical Cloud Forest Catchments: Combining Descriptive Approaches and Multivariate Analysis

We investigated controls on the water chemistry of a South Ecuadorian cloud forest catchment which is partly pristine, and partly converted to extensive pasture. From April 2007 to May 2008 water samples were taken weekly to biweekly at nine different subcatchments, and were screened for differences in electric conductivity, pH, anion, as well as element composition. A principal component analysis was conducted to reduce dimensionality of the data set and define major factors explaining variation in the data. Three main factors were isolated by a subset of 10 elements (Ca²⁺, Ce, Gd, K⁺, Mg²⁺, Na⁺, Nd, Rb, Sr, Y), explaining around 90% of the data variation. Land-use was the major factor controlling and changing water chemistry of the subcatchments. A second factor was associated with the concentration of rare earth elements in water, presumably highlighting other anthropogenic influences such as gravel excavation or road construction. Around 12% of the variation was explained by the third component, which was defined by the occurrence of Rb and K and represents the influence of vegetation dynamics on element accumulation and wash-out. Comparison of base- and fast flow concentrations led to the assumption that a significant portion of soil water from around 30 cm depth contributes to storm flow, as revealed by increased rare earth element concentrations in fast flow samples. Our findings demonstrate the utility of multi-tracer principal component analysis to study tropical headwater streams, and emphasize the need for effective land management in cloud forest catchments.     

A new computational strategy for solving two‐phase flow in strongly heterogeneous poroelastic media of evolving scales

We develop a new computational methodology for solving two‐phase flow in highly heterogeneous porous media incorporating geomechanical coupling subject to uncertainty in the poromechanical parameters. Within the framework of a staggered‐in‐time coupling algorithm, the numerical method proposed herein relies on a Petrov–Galerkin postprocessing approach projected on the Raviart–Thomas space to compute the Darcy velocity of the mixture in conjunction with a locally conservative higher order finite volume discretization of the nonlinear transport equation for the saturation and an operator splitting procedure based on the difference in the time‐scales of transport and geomechanics to compute the effects of transient porosity upon saturation. Notable features of the numerical modeling proposed herein are the local conservation properties inherited by the discrete fluxes that are crucial to correctly capture the fingering patterns arising from the interaction between heterogeneity and nonlinear viscous coupling. Water flooding in a poroelastic formation subject to an overburden is simulated with the geology characterized by multiscale self‐similar permeability and Young modulus random fields with power‐law covariance structure. Statistical moments of the poromechanical unknowns are computed within the framework of a high‐resolution Monte Carlo method. Numerical results illustrate the necessity of adopting locally conservative schemes to obtain reliable predictions of secondary recovery and finger growth in strongly heterogeneous deformable reservoirs. Copyright © 2011 John Wiley & Sons, Ltd.