The FFT Moving Average (FFT-MA) Generator: An Efficient Numerical Method for Generating and Conditioning Gaussian Simulations

A fast Fourier transform (FFT) moving average (FFT-MA) method for generating Gaussian stochastic processes is derived. Using discrete Fourier transforms makes the calculations easy and fast so that large random fields can be produced. On the other hand, the basic moving average frame allows us to uncouple the random numbers from the structural parameters (mean, variance, correlation length, ... ), but also to draw the randomness components in spatial domain. Such features impart great flexibility to the FFT-MA generator. For instance, changing only the random numbers gives distinct realizations all having the same covariance function. Similarly, several realizations can be built from the same random number set, but from different structural parameters. Integrating the FFT-MA generator into an optimization procedure provides a tool theoretically capable to determine the random numbers identifying the Gaussian field as well as the structural parameters from dynamic data. Moreover, all or only some of the random numbers can be perturbed so that realizations produced using the FFT-MA generator can be locally updated through an optimization process.     

Virtual MHD Jets on Grids

As network performance has outpaced computational power and storage capacity, a new paradigm has evolved to enable the sharing of geographically distributed resources. This paradigm is known as Grid computing and aims to offer access to distributed resource irrespective of their physical location. Many national, European and international projects have been launched during the last years trying to explore the Grid and to change the way we are doing our everyday work. In Ireland, we have started the CosmoGrid project that is a collaborative project aimed to provide high performance super-computing environments. This will help to address complex problems such as magnetohydrodynamic outflows and jets in order to model and numerically simulate them. Indeed, the numerical modeling of plasma jets requires massive computations, due to the wide range of spatial-temporal scales involved. We present here the first jet simulations and their corresponding models that could help to understand results from laboratory experiments.     

Spinel lherzolite xenoliths from the Premier kimberlite (Kaapvaal craton, South Africa): Nature and evolution of the shallow upper mantle beneath the Bushveld complex

Coarse-grained, granular spinel lherzolites xenoliths from the Premier kimberlite show evidence of melt extraction and metasomatic enrichment, documenting a complex history for the shallow mantle beneath the Bushveld complex. Compositions of orthopyroxene, clinopyroxene and spinel indicate equilibration within the spinel–peridotite facies of the upper mantle, at depths from 80 to 100 km and temperatures from 720 to 850 °C. Bulk compositions have lower Mg-number [atomic 100 Mg/(Mg + Fe*)] than previously studied spinel peridotites from Premier, and have higher Mg/Si than low-temperature coarse grained garnet lherzolites, suggesting shallower melting conditions or metasomatic enrichment. Clinopyroxene in one sample is highly LREE-depleted indicating very minor modification of a residue of 20% melt extraction, whereas the calculated REE pattern for a melt in equilibrium with a mildly LREE-depleted sample is similar to MORB or late Archean basalt, possibly related to the Bushveld Complex. Bulk and mineral compositions suggest minimal refertilization by silicate melts in four out of six samples, but REE patterns indicate introduction of a LIL-enriched component that may be related to kimberlite.     

Organic matter and metamorphic history of CO chondrites

The metamorphic grades of a series of eight CO3 chondrites (ALHA77307, Colony, Kainsaz, Felix, Lancé, Ornans, Warrenton and Isna) have been quantified. The method used was based on the structural grade of the organic matter trapped in the matrix, which is irreversibly transformed by thermal metamorphism. The maturation of the organic matter is independent with respect to the mineralogical context and aqueous alteration. This metamorphic tracer is thus valid whatever the chemical class of chondrites. Moreover, it is sensitive to the peak metamorphic temperature.The structural grade of the organic matter was used along with other metamorphic tracers such as petrography of opaque minerals, Fa and Fs silicate composition in type I chondrules, presolar grains and noble gas (P3 component) abundance. The deduced metamorphic hierarchy and the attributed petrographic types are the following: ALHA77307 (3.03) < Colony (3.1) < Kainsaz (3.6) < Felix (3.6 (1)) < Ornans (3.6 (2)) < Lancé (3.6 (3)) < Warrenton (3.7 (1)) < Isna (3.7 (2)).For most metamorphosed objects, the peak metamorphic temperature can be estimated using a geothermometer calibrated with terrestrial metasediments [Beyssac O., Goffe B., Chopin C., and Rouzaud J. N. (2002) Raman spectrum of carbonaceous material in metasediments: a new geothermometer. J. Metamorph. Geol., 20, 859-871]. A value of 330 °C was obtained for Allende (CV chondrite), Warrenton and Isna, consistent with temperatures estimated from Fe diffusion [Weinbruch S., Armstrong J., and Palme H. (1994). Constraints on the thermal history of the Allende parent body as derive from olivine-spinel thermometry and Fe/Mg interdiffusion in olivine. Geochim. Cosmochim. Acta58(2), 1019-1030.], from the Ni content in sulfide-metal assemblages [Zanda B., Bourot-Denise M., and Hewins R. (1995) Condensate sulfide and its metamorphic transformations in primitive chondrites. Meteorit. Planet. Sci.30, A605.] and from the d₀₀₂ interlayer spacing in poorly graphitized carbon [Rietmeijer, F., and MacKinnon, I. (1985) Poorly graphitized carbon as a new cosmothermometer for primitive extraterrestrial materials. Nature, 315, 733-736]. The trapped noble gas and C content appear to be sensitive but not precise metamorphic tracers, indicating that the “Ornans paradox” does not exist. Major problems with the current petrologic types derived from Induced ThermoLuminescence are pointed out.     

Optimization with the Gradual Deformation Method

Building reservoir models consistent with production data and prior geological knowledge is usually carried out through the minimization of an objective function. Such optimization problems are nonlinear and may be difficult to solve because they tend to be ill-posed and to involve many parameters. The gradual deformation technique was introduced recently to simplify these problems. Its main feature is the preservation of the spatial structure: perturbed realizations exhibit the same spatial variability as the starting ones. It is shown that optimizations based on gradual deformation converge exponentially to the global minimum, at least for linear problems. In addition, it appears that combining the gradual deformation parameterization with optimizations may remove step by step the structure preservation capability of the gradual deformation method. This bias is negligible when deformation is restricted to a few realization chains, but grows increasingly when the chain number tends to infinity. As in practice, optimization of reservoir models is limited to a small number of iterations with respect to the number of gridblocks, the spatial variability is preserved. Last, the optimization processes are implemented on the basis of the Levenberg–Marquardt method. Although the objective functions, written in terms of Gaussian white noises, are reduced to the data mismatch term, the conditional realization space can be properly sampled.     

Estimating variation in stomatal frequency at intra-individual,intra-site,and inter-taxonomic levels in populations of the Leonardoxa africana (Fabaceae) complex over environmental gradients in Cameroon

Variation of stomatal frequency (stomatal density and stomatal index) includes genetically-based, potentially-adaptive variation, and variation due to phenotypic plasticity, the degree of which may be fundamental to the ability to maintain high water-use efficiency and thus to deal with environmental change. We analysed stomatal frequency and morphology (pore length, pore width) in leaves from several individuals from nine populations of four sub-species of the Leonardoxa africana complex. The dataset represents a hierarchical sampling wherein factors are nested within each level (leaves in individuals, individuals in sites, etc.), allowing estimation of the contribution of different levels to overall variation, using variance-component analysis. SI showed significant variation among sites (“site” is largely confounded with “sub-species”), being highest in the sub-species localized in the highest-elevation site. However, most of the observed variance was accounted for at intra-site and intra-individual levels. This variance could reflect great phenotypic plasticity, presumably in response to highly local variation in micro-environmental conditions.     

Changes in climate extremes, variability and signature on sub-Antarctic Marion Island

The ecological consequences of climate change are determined by many climate parameters, not just by the commonly investigated changes in mean temperature and rainfall. More comprehensive studies, including analyses of climate variability, extremes and aggregate changes in the climate system, can improve the understanding of the nature, and therefore possible consequences, of recent changes in climate. Here climate trends on the sub-Antarctic Marion Island are documented (between 1949 and 2003) in more detail than previously. Significant trends in biologically-relevant, and previously unexplored, parameters were observed, and the potential ecological consequences of these changes discussed. For example, the decline in precipitation experienced on the island comprises a trend for longer dry spells punctuated by fewer and smaller precipitation events. This more detailed understanding of the island’s drying trends enables more accurate predictions about its impacts, including, for example, particularly severe effects on plant species growing in soils with poor water-holding capacity. Therefore, in addition to changes in average conditions, more inclusive climate analyses should also examine trends in climatic variability and extremes, for individual climate parameters as well as for the climate system as a whole.     

Directional wave spectra at the regional scale with the KuROS airborne radar: comparisons with models

Ocean Dynamics - In situ observations, satellite observations, and regional observations from airborne remote sensing are very useful to characterize sea state evolution and related physical...     

Gravity and S-wave modelling across the Jan Mayen Ridge,North Atlantic; implications for crustal lithology

The horizontal components from fourteen Ocean Bottom Seismometers deployed along four profiles focused along the western margin of the Jan Mayen microcontinent, North Atlantic, have been modelled with regard to S-waves, based on P-wave models obtained earlier. The seismic models have furthermore been constrained by 2D gravity modelling. High V _p/V _s-ratios (2.3–7.9) within the Cenozoic sedimentary section are attributed to significant porosities, whereas V _p/V _s-ratios in the order of 1.9–2.2 for the Mesozoic and Paleozoic sedimentary rocks indicate shale-dominated lithology throughout the area. The eastern side of the Jan Mayen Ridge is interpreted as a passive, volcanic margin, based on relatively high crustal V _p/V _s-ratios (1.9), whereas lower V _p/V _s-ratios (1.75–1.8) suggest the presence of intermediate composition crust and non-volcanic margin on the western side of the ridge. In the westernmost part of the Jan Mayen Basin, slightly increased upper mantle V _p/V _s-ratios may indicate some degree of serpentization of upper mantle peridotites.