Multivariate Intrinsic Random Functions for Cokriging

In multivariate geostatistics, suppose that we relax the usual second-order-stationarity assumptions and assume that the component processes are intrinsic random functions of general orders. In this article, we introduce a generalized cross-covariance function to describe the spatial cross-dependencies in multivariate intrinsic random functions. A nonparametric method is then proposed for its estimation. Based on this class of generalized cross-covariance functions, we give cokriging equations for multivariate intrinsic random functions in the presence of measurement error. A simulation is presented that demonstrates the accuracy of the proposed nonparametric estimation method. Finally, an application is given to a dataset of plutonium and americium concentrations collected from a region of the Nevada Test Site used for atomic-bomb testing.     

Thermophysical properties of the α–β–γ polymorphs of Mg<Subscript>2</Subscript>SiO<Subscript>4</Subscript>: a computational study

Thermophysical properties of the various polymorphs (i.e. α-, β- and γ) of Mg₂SiO₄ were computed with the CRYSTAL06 code within the framework of CO-LCAO-GTF approach by using the hybrid B3LYP density functional method. Potential wells were calculated through a symmetry preserving, variable cell-shape structure relaxation procedure. Vibrational frequencies were computed at the long-wavelength limit corresponding to the center of the Brillouin zone (k → 0). Thermodynamic properties were estimated through a semiclassical approach that combines B3LYP vibrational frequencies for optic modes and the Kieffer’s model for the dispersion relation of acoustic modes. All computed values except volume (i.e. electronic energy, zero point energy, optical vibrational modes, thermal corrections to internal energy, standard state enthalpy and Gibbs free energy of reaction, bulk modulus and its P and T derivatives, entropy, C _V, C _P) are consistent with available experimental data and/or reasonable estimates. Volumes are slightly overestimated relative to those determined directly by X-ray diffraction. A set of optimized volumetric properties that are consistent with the other semiclassical properties of the phases α, β and γ have been derived by optimization procedure such that the calculated boundaries for the α/β and β/γ equilibria have the best overall agreement with the experimental data for these transitions. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.

Steerable wavelet analysis of CMB structures alignment

This paper reviews the application of a novel methodology for analysing the isotropy of the universe by probing the alignment of local structures in the CMB. The strength of the proposed methodology relies on the steerable wavelet filtering of the CMB signal. One the one hand, the filter steerability renders the computation of the local orientation of the CMB features affordable in terms of computation time. On the other hand, the scale-space nature of the wavelet filtering allows to explore the alignment of the local structures at different scales, probing possible different phenomena. We present the WMAP first-year data analysis recently performed by the same authors (Wiaux et al.), where an extremely significant anisotropy was found. In particular, a preferred plane was detected, having a normal direction with a northern end position at (θ, ) = (34°, 331°), close to the northern end of the CMB dipole axis. In addition, a most preferred direction was found in that plane, with a northern end direction at (θ, ) = (71°, 91°), very close to the north ecliptic pole. This result synthesised for the first time previously reported anomalies identified in the direction of the dipole and the ecliptic poles axes. In a forthcoming paper (Vielva et al.), we have extended our analysis to the study of individual frequency maps finding first indications for discarding foregrounds as the origin of the anomaly. We have also tested that the preferred orientations are defined by structures homogeneously distributed in the sky, rather than from localised regions. We have also analysed the WMAP 3-year data, finding the same anomaly pattern, although at a slightly lower significance level.     

Water temperature modelling in a controlled environment: comparative study of heat budget equations

One of the challenges when modelling a complex variable such as water temperature in rivers is that it can be difficult to determine the sources of error and to ensure that the simulations are truly representative of the reality. Therefore, a heat budget study was completed in a controlled environment, which excluded advection and bottom fluxes but enabled observation of all the other fluxes. A 21.42 m³ pool was installed and insulated to limit heat exchange through the sides and bottom. All the major energy fluxes were monitored for a 50‐day period. Different equations for individual heat budget terms were tested to determine their ability to reproduce the observations. This experiment also permitted to assess the relative importance of each component of the heat budget. Performance of each semi‐empirical equation was determined by comparing predictions and measured values. It was thus possible to choose the formulae that best represented the measured heat exchange processes, while understanding the limits of some of the semi‐empirical representations of heat exchange processes. The results highlight the importance of radiative terms into the heat budget because they controlled the major sources and sinks. The study also showed the importance of the wind function determination into the calculation of latent heat flux. The resulting water temperature model returned simulated hourly water temperature with an overall root mean square error of 0.71 °C/h and a modified Nash–Sutcliffe coefficient of 0.97. Copyright © 2012 John Wiley & Sons, Ltd.     

Removing attitude-related variations in the line-of-sight for kinematic GPS positioning

We are presenting a method for removing attitude-related variations (ATTRVs) in dynamic 1 Hz GPS positioning. The ATTRVs are separated into vertical and horizontal components. These result from the translational and rotational motions of the vehicle, which is a marine research vessel in our case. We have developed new observation equations that use corrected pseudoranges and carrier phases to account for ATTRVs. In the present contribution, we are only focusing on the vertical signals. These modeled ATTRVs are included as corrections in the line-of-sight (LoS) to each GPS satellite in the RINEX data sets using a new software called RNXATTCOR. Precise IGS sp3-orbits are used as inputs together with precise lever arm coordinates of the onboard GPS antennas, observations from the marine inertial navigation system and a priori 3D position of the vessel. The corrected RINEX data sets are then processed using kinematic processing or sequential processing in Precise Point Positioning (PPP) mode. The reduction of the standard deviation from a running mean in the final processed GPS time series is about 95%. The current method is being proposed for marine geodesy science applications such as along-track calibration/validation of altimetry products, and also for applications related to offshore kinematic high precision GPS positioning such as drilling, offshore platforms stability, marine pipeline positioning, earthquake monitoring and tsunami detection.     

CBM and CO2-ECBM related sorption processes in coal: A review

This article reviews the state of research on sorption of gases (CO₂, CH₄) and water on coal for primary recovery of coalbed methane (CBM), secondary recovery by an enhancement with carbon dioxide injection (CO₂-ECBM), and for permanent storage of CO₂ in coal seams.Especially in the last decade a large amount of data has been published characterizing coals from various coal basins world-wide for their gas sorption capacity. This research was either related to commercial CBM production or to the usage of coal seams as a permanent sink for anthropogenic CO₂ emissions. Presently, producing methane from coal beds is an attractive option and operations are under way or planned in many coal basins around the globe. Gas-in-place determinations using canister desorption tests and CH₄ isotherms are performed routinely and have provided large datasets for correlating gas transport and sorption properties with coal characteristic parameters.Publicly funded research projects have produced large datasets on the interaction of CO₂ with coals. The determination of sorption isotherms, sorption capacities and rates has meanwhile become a standard approach.In this study we discuss and compare the manometric, volumetric and gravimetric methods for recording sorption isotherms and provide an uncertainty analysis. Using published datasets and theoretical considerations, water sorption is discussed in detail as an important mechanisms controlling gas sorption on coal. Most sorption isotherms are still recorded for dry coals, which usually do not represent in-seam conditions, and water present in the coal has a significant control on CBM gas contents and CO₂ storage potential. This section is followed by considerations of the interdependence of sorption capacity and coal properties like coal rank, maceral composition or ash content. For assessment of the most suitable coal rank for CO₂ storage data on the CO₂/CH₄ sorption ratio data have been collected and compared with coal rank.Finally, we discuss sorption rates and gas diffusion in the coal matrix as well as the different unipore or bidisperse models used for describing these processes.This review does not include information on low-pressure sorption measurements (BET approach) to characterize pore sizes or pore volume since this would be a review of its own. We also do not consider sorption of gas mixtures since the data base is still limited and measurement techniques are associated with large uncertainties.     

Solubility of CO2 in a Ca-rich leucitite: effects of pressure,temperature, and oxygen fugacity

The solubility of carbon dioxide in a Ca-rich leucitite has been investigated as a function of pressure (0.1–2.0 GPa), temperature (1200–1600°C), and oxygen fugacity. The experiments were done in a rapid-quench internally-heated pressure vessel (0.1 GPa) and a piston cylinder (0.5–2.0 GPa). The leucitite glass, previously equilibrated at NNO, and silver oxalate were loaded in Fe-doped Pt capsules (oxidized conditions) and graphitelined Pt capsules (reduced conditions). Secondary Ion Mass Spectrometry and bulk carbon analyses were used to determine the amount of dissolved carbon. Speciation of carbon was characterized by Fourier transform microinfrared spectroscopy. At oxidized conditions, only CO₃ ^2- is observed as a dissolved species. The solubility is high with CO₂ contents in the melt attaining 6.2 wt% at 2.0 GPa and 1350°C. The solubility increases with pressure and shows a significant negative temperature dependence. An excellent correlation is obtained when the data are fit to a model, based on the simplified solubility reaction CO₂ ^(vapor)+O^2-(melt)CO₃ ^2-(melt), which describes the solubility of CO₂ as a function of pressure, temperature and fCO₂. At reduced conditions, the amount of carbon dissolved is significantly lower, and CO₃ ^2- is still the only species present in the melt. If the solubility model established at oxidized conditions is applied, the carbon dissolved appears to be essentially a function of fCO₂ alone although divergence increases in a consistent manner with pressure and temperature. This could suggest a low but significant solubility of CO with a positive temperature dependence or a departure of the calculated fluid compositions determined by the equation of state from the actual ones. The strong preferential solubility of carbon in its oxidized C⁴⁺ form, even at reduced conditions, implies that ascending melts with high CO₂ solubility can experience significant oxidation through degassing. This could reconcile the oxidized nature of some Ca-rich alkaline magmas with more reduced mantle source regions.     

Distribution of222Rn over the north Pacific: Implications for continental influences

The atmospheric distribution of²²²Rn over the north Pacific is simulated with a three-dimensional chemical tracer model using meteorological input from the NASA-GISS general circulation model (4°×5° resolution). Radon-222 (half-life 3.8 days) is a tracer of continental air. Model results are in good agreement with measurements from ships and aircraft. Strong Asian influence is found throughout the tropospheric column over the north Pacific in spring, reflecting a combination of frequent convection over the continent, strong westerly winds at altitude, and subsidence over the ocean. In summer, the upper troposphere over the north Pacific is heavily affected by deep convection over China; however, Asian influences at the surface are then at their yearly minimum. In winter, strong Asian influence is found near the surface but not at high altitudes. Transport of American air over the Pacific is important only at low latitudes. American sources account for 11% of total²²²Rn in the model at Midway, 30% at Mauna Loa and 59% at Oahu. Results for Hawaii indicate two seasonal peaks of American influence, one in summer and one in winter. The tropical western Pacific is particularly remote from continental influences year round.