Application of Multiple Point Geostatistics to Non-stationary Images

Simulation of flow and solute transport through aquifers or oil reservoirs requires a precise representation of subsurface heterogeneity that can be achieved by stochastic simulation approaches. Traditional geostatistical methods based on variograms, such as truncated Gaussian simulation or sequential indicator simulation, may fail to generate the complex, curvilinear, continuous and interconnected facies distributions that are often encountered in real geological media, due to their reliance on two-point statistics. Multiple Point Geostatistics (MPG) overcomes this constraint by using more complex point configurations whose statistics are retrieved from training images. Obtaining representative statistics requires stationary training images, but geological understanding often suggests a priori facies variability patterns. This research aims at extending MPG to non-stationary facies distributions. The proposed method subdivides the training images into different areas. The statistics for each area are stored in separate frequency search trees. Several training images are used to ensure that the obtained statistics are representative. The facies probability distribution for each cell during simulation is calculated by weighting the probabilities from the frequency trees. The method is tested on two different object-based training image sets. Results show that non-stationary training images can be used to generate suitable non-stationary facies distributions.     

High-resolution mm interferometry and the search for massive protostellar disks: the case of Cep-A HW2

The direct detection of accretion onto massive protostars through rotating disks constitutes an important tile in the massive-star-formation-theory mosaic. This task is however observationally very challenging. A very interesting example is Cepheus A HW2. The properties of the molecular emission around this YSO seems to suggest the presence of a massive rotating disk (cf. Patel et al. in Nature 437:109, 2005). We have carried out sub-arcsec-resolution PdBI observations of high-density and shock tracers such as SO₂, SiO, CH₃CN, and CH₃OH towards the center of the outflow. A detailed analysis of the spatial distribution and of the velocity field traced by all observed species leads us to conclude that, on a ∼700 AU scale, the Cep-A “disk” is actually the result of the superposition of multiple hot-core-type objects, at least one of them ejecting an outflow at a small angle with respect to the line of sight. Together with the well-known large-scale outflow ejected by HW2, this setup makes for a very complex spatial and kinematic picture. Based on observations carried out with the IRAM Plateau de Bure Interferometer. IRAM is supported by INSU/CNRS (France), MPG (Germany) and IGN (Spain).     

The 2D Steady Hydraulic Head Field Surrounding a Pumping Well in a Finite Heterogeneous Confined Aquifer

We present a second-order analytic solution [in terms of a heterogeneous log-transmissivity Y(r) = ln T(r)] for the hydraulic head field in a finite 2D confined heterogeneous aquifer under steady radial flow conditions assuming fixed head boundary conditions at the well and at a circular exterior boundary. The solution may be used to obtain the gradient used in calculation of solute transport to a well in a heterogeneous transmissivity field. The solution, obtained using perturbation methods coupled with Green's function techniques, leads us to postulate a more general form of the head for arbitrarily large-variance fields and may be used to obtain moment relations between the log-transmissivity and head under convergent flow conditions when Y(r) is expressed as a random space function. We present expressions for the mean head field when the log-transmissivity is Gaussian and conditioned on the transmissivity value at the well for an arbitrary ln T covariance. Finally, we look at the effect of parameter variations on the mean head behavior and present numerical simulations verifying the second-order mean head expressions.     

An ocean observation system for monitoring the affects of climate change on the ecology and sustainability of pelagic fisheries in the Pacific Ocean

Climate change presents an emerging challenge to the sustainable management of tuna fisheries, and robust information is essential to ensure future sustainability. Climate and harvest affect tuna stocks, populations of non-target, dependent species and the ecosystem. To provide relevant advice we need an improved understanding of oceanic ecosystems and better data to parameterise the models that forecast the impacts of climate change. Currently ocean-wide data collection in the Pacific Ocean is primarily restricted to oceanographic data. However, the fisheries observer programs that operate in the region offer an opportunity to collect the additional information on the mid and upper trophic levels of the ecosystem that is necessary to complement this physical data, including time-series of distribution, abundance, size, composition and biological information on target and non-target species and mid trophic level organisms. These observer programs are in their infancy, with limited temporal and spatial distribution but recent international and national policy decisions have been made to expand their coverage. We identify a number of actions to initiate this monitoring including: consolidating collaborations to ensure the use of best quality data; developing consistency between sub-regional observer programmes to ensure that they meet the objectives of ecosystem monitoring; interrogating of existing time series to determine the most appropriate spatial template for monitoring; and exploring existing ecosystem models to identify suitable indicators of ecosystem status and change. The information obtained should improve capacity to develop fisheries management policies that are resilient and can be adapted to climate change.     

Lower Carboniferous post-orogenic granites in central-eastern Sierra de Velasco,Sierras Pampeanas,Argentina: U–Pb monazite geochronology,geochemistry and Sr–Nd isotopes

The central-eastern part of the Sierra de Velasco (Sierras Pampeanas, NW Argentina) is formed by the large Huaco (40 × 30 km) and Sanagasta (25 × 15 km) granite massifs and the small La Chinchilla stock (2 × 2 km). The larger granites intrude into Ordovician metagranitoids and crosscut Devonian (?) mylonitic shear zones, whereas the small stock sharply intrudes into the Huaco granite. The two voluminous granites are biotitic-muscovitic and biotitic porphyritic syeno- to monzogranites. They contain small and rounded tonalitic and quartz-dioritic mafic microgranular enclaves. The small stock is an equigranular, zinnwaldite- and fluorite-bearing monzogranite. The studied granites are silica-rich (SiO₂ >70%), potassium-rich (K₂O >4%), ferroan, alkali-calcic to slightly calk-alkalic, and moderately to weakly peraluminous (A/CNK: 1.06–1.18 Huaco granite, 1.01–1.09 Sanagasta granite, 1.05–1.06 La Chinchilla stock). They have moderate to strong enrichments in several LIL (Li, Rb, Cs) and HFS (Nb, Ta, Y, Th, U) elements, and low Sr, Ba and Eu contents. U–Pb monazite age determinations indicate Lower Carboniferous crystallization ages: 350–358 Ma for the Huaco granite, 352.7 ± 1.4 Ma for the Sanagasta granite and 344.5 ± 1.4 Ma for the La Chinchilla stock. The larger granites have similar ?Nd values between ?2.1 and ?4.3, whereas the younger stock has higher ?Nd of ?0.6 to ?1.4, roughly comparable to the values obtained for the Carboniferous San Blas granite (?1.4 to ?1.7), located in the north of the sierra. The Huaco and Sanagasta granites have a mainly crustal source, but with some participation of a more primitive, possibly mantle-derived, component. The main crustal component can be attributed to Ordovician peraluminous metagranitoids. The La Chinchilla stock derives from a more primitive source, suggesting an increase with time in the participation of the primitive component during magma genesis. The studied granites were generated during a post-orogenic period in a within-plate setting, possibly as a response to the collapse of the previous Famatinian orogen, extension of the crust and mantle upwelling. They are part of the group of Middle Devonian–Lower Carboniferous granites of the Sierras Pampeanas. The distribution and U–Pb ages of these granites suggests a northward arc-parallel migration of this mainly post-orogenic magmatism with time.     

Multiple Iron-Rich Spinel Phases and Hematite in a Magnetic Soil Developing on Tuffite

A dusky red Latosol developing on tuffite was examined in detail in order to characterize its iron oxide mineralogy. Magnetic separation and selective chemical treatment were employed. The mineralogical analyses were carried out by powder X-ray diffraction and Mössbauer spectroscopy, with the electron microprobe technique being used to assess the chemical composition of selected grains. Results reveal that iron-rich spinel-structure phases from this pedosystem are rather variable in chemical composition and lattice parameters, and can be grouped in at least three major populations: the first group consisting of the iron-richest phase, with very low Ti and Cr and virtually no Mg content and high Al; the second group presenting Ti>Mg>>Al>Mn contents but no Cr and the third group related to the Mg-richest phase with Mg>Ti>>Al>Mn. Two major populations of hematite were also identified: one with ～15 mol% Al and other practically Al-free. The first group is thought to be a product of neoformation in an Al-rich pedogenic medium whereas the second is likely to be transformed directly from maghemite. Though chemical treatment along with physical techniques provide a powerful way to characterize iron oxides of tropical soils, some methodological difficulties still remain, particularly those concerned with (i) the characterization and quantification of maghemite in presence of relatively large amount of hematite in the soil mixture, particularly by Mössbauer spectroscopy, and (ii) the allocation of the isomorphous substituents for iron in the cation sites of the spinel structure. Possible mechanisms of alteration of iron-rich phases in this pedosystem are discussed.     

Carbon fluxes resulting from land-use changes in the Tamaulipan thornscrub of northeastern Mexico

Information on carbon stock and flux resulting from land-use changes in subtropical, semi-arid ecosystems are important to understand global carbon flux, yet little data is available. In the Tamaulipan thornscrub forests of northeastern Mexico, biomass components of standing vegetation were estimated from 56 quadrats (200 m² each). Regional land-use changes and present forest cover, as well as estimates of soil organic carbon from chronosequences, were used to predict carbon stocks and fluxes in this ecosystem.     

Semi‐Micro Reflux Procedure for Minimization of Chloride Interference by COD Determination

In this study, a semi‐micro procedure for the COD (Chemical Oxygen Demand) was undertaken to allow the measurement of effluent samples with a high chloride content. It was found that the COD values for standard and effluent samples, corresponded closely to those obtained using the standard method of Burns and Marshall. When the method was applied to samples of industrial wastewater, with a very high chloride content, the recommended procedure proved to be accurate and precise, with RSDs of 6% and 4.3% for the standard and effluent solutions, respectively. The required amount of hazardous reagents was similar to that found in the standard closed reflux method. It was also possible to determine the total amount of oxidized halides. The use of this low‐cost, laboratory‐made system made it unnecessary to acquire an expensive, imported analytical system, while maintaining a high level of accuracy in the determinations.     

Asymptotic solution for the two-body problem with constant tangential thrust acceleration

An analytical solution of the two body problem perturbed by a constant tangential acceleration is derived with the aid of perturbation theory. The solution, which is valid for circular and elliptic orbits with generic eccentricity, describes the instantaneous time variation of all orbital elements. A comparison with high-accuracy numerical results shows that the analytical method can be effectively applied to multiple-revolution low-thrust orbit transfer around planets and in interplanetary space with negligible error.     

Stream-stage response tests and their joint interpretation with pumping tests

Hydraulic head response to stream-stage variations can be used to explore the hydraulic properties of stream-aquifer systems at a relatively large scale. These stream-stage response tests, also called flooding tests, are typically interpreted using one- or two-dimensional models that assume flow perpendicular to the river. Therefore, they cannot be applied to systems that are both horizontally and vertically heterogeneous. In this work, we use the geostatistical inverse problem to jointly interpret data from stream-stage response and pumping tests. The latter tests provide flow data (which are needed to resolve aquifer diffusivity into transmissivity and storage coefficient) and may supply supplementary small-scale information. Here, we summarize the methodology for the design, execution, and joint numerical interpretation of these tests. Application to the Aznalcóllar case study allows us to demonstrate that the proposed methodology may help in responding to questions such as the continuity of aquitards, the role and continuity of highly permeable paleochannels, or the time evolution of stream-aquifer interaction. These results expand the applicability and scope of stream-stage response tests.