Analysis and Estimation of Natural Processes with Nonhomogeneous Spatial Variation Using Secondary Information

Natural processes encountered in mining, hydrogeologic, environmental, etc. applications usually are poorly known because of scarcity of data over the area of interest. Therefore, stochastic estimation techniques are the tool of choice for a careful accounting of the heterogeneity and uncertainty involved. Within such a framework, a better utilization of all available data concerning the process of interest and all other natural processes related to it, is of primary importance. Because many natural processes show complicated spatial trends, the hypothesis of spatial homogeneity cannot be invoked always, and the more general theory of intrinsic spatial random fields should be employed. Efficient use of secondary information in terms of the intrinsic model requires that suitable permissibility criteria for the generalized covariances and cross-covariances are satisfied. A set of permissibility criteria are presented for the situation of two intrinsic random fields. These criteria are more general and comprehensive than the ones currently available in the geostatistical literature. A constrained least-square technique is implemented for the inference of the generalized covariance and cross-covariance parameters, and a synthetic example is used to illustrate the methodology. The numerical results show that the use of secondary information can lead to significant reductions in the estimation errors.     

The role of preconditioning in the solution to FE coupled consolidation equations by Krylov subspace methods

The repeated solution in time of the linear system arising from the finite element integration of coupled consolidation equations is a major computational effort. This system can be written in either a symmetric or an unsymmetric form, thus calling for the implementation of different preconditioners and Krylov subspace solvers. The present paper aims at investigating when either a symmetric or an unsymmetric approach should be better used. The results from a number of representative numerical experiments indicate that a major role in selecting either form is played by the preconditioner rather than by the Krylov subspace method itself. Two other important issues addressed are the size of the time integration step and the possible lumping of the flow capacity matrix. It appears that ad hoc block constrained preconditioners provide the most robust algorithm independently of the time step size, lumping, and symmetry. Copyright © 2008 John Wiley & Sons, Ltd.     

Recent Trends Towards Oligotrophication of the Northern Adriatic: Evidence from Chlorophyll <Emphasis Type="Italic">a</Emphasis> Time Series

The results of the updated and quality-checked data base of field observations on chlorophyll a (Chl a) collected in the period 1970–2007 in the Northern Adriatic Sea are presented. From the last decade, SeaWiFS satellite information was also considered. Results demonstrate a global tendency towards Chl a reduction in the period of investigation, which is more marked in the eutrophic area under the influence of the Po River. In the rest of the basin, which presents meso- or oligotrophic characteristics, long-term changes are more difficult to detect. The long-term field dataset can be divided into two periods: the last decade characterized by the strong decrease observed in the whole northern Adriatic and the earlier period with no or slight increase. The recent substantial reduction of Chl a concentrations is confirmed all over the basin (−0.11 mg m⁻³ year⁻¹) from satellite-derived information. Results are consistent with recently evidenced decrease in concentrations of phosphate and ammonia and point to the existence of oligotrophication in the Northern Adriatic. Results indicate forcefully that the still common perception of the Adriatic Sea as a very eutrophic basin is no longer appropriate, at least for its northern part and in recent years.     

COCHISE: the first light of the Italian millimetre telescope at Concordia (Dome C, Antarctica)

COCHISE (Cosmological Observations at Concordia with High-sensitivity Instrument for Source Extraction) is a 2.6 m telescope located on the high Antarctic Plateau near the Italian–French Concordia Base. The telescope is mainly devoted to Cosmological observations, able to operate between 200 μm and 3 mm of wavelength. In this paper we describe the main characteristics of the instrument. We also report on the first light, obtained during summer 2010–2011: this result marks the beginning of millimetre astrophysical observations at Concordia. Responsivity, noise equivalent temperature and field of view of the instrument are reported. At present COCHISE is the largest telescope located at Concordia. Beside the scientific expectations, the use of this kind of instrument in the Antarctic environment poses technological aspects of relevant interest: thus COCHISE can be considered as a pathfinder for future Antarctic telescopes.     

On the use of infinite elements in the analysis of two-dimensional layered elastic systems via discretized integral equations

A boundary integral equation approach is presented for the two-dimensional plane strain analysis of horizontally layered elastic systems resting on a rough and rigid base. The validity of Somigliana identity for unbounded layers is discussed and an Infinite Boundary Element is proposed in order to take into account the lateral unboundedness of the physical problem in the discretized equations. This element is implemented in a computer code based on a successive stiffness solution procedure. A series of numerical tests concerning single- and multi-layered problems is illustrated, and the performance of the proposed solution method is compared with that of more traditional boundary element and finite element techniques.     

Direct,partitioned and projected solution to finite element consolidation models

Direct, partitioned, and projected (conjugate gradient‐like) solution approaches are compared on unsymmetric indefinite systems arising from the finite element integration of coupled consolidation equations. The direct method is used in its most recent and computationally efficient implementations of the Harwell Software Library. The partitioned approach designed for coupled problems is especially attractive as it addresses two separate positive definite problems of a smaller size that can be solved by symmetric conjugate gradients. However, it may stagnate and when converging it does not prove competitive with a global projection method such as Bi‐CGSTAB, which may take full advantage of its flexibility in working on scaled and reordered equations, and thus may greatly improve its computational performance in terms of both robustness and convergence rate. The Bi‐CGSTAB superiority to the other approaches is discussed and demonstrated with a few representative examples in two‐dimensional (2‐D) and three‐dimensional (3‐D) coupled consolidation problems. Copyright © 2002 John Wiley & Sons, Ltd.     

Mantle-derived carbon in Hercynian granites. Stable isotopes signatures and C/He associations in the thermomineral waters,N-Portugal

Na–HCO₃–CO₂-rich thermomineral waters issue in the N of Portugal, within the Galicia-Trás-os-Montes region, linked to a major NNE-trending fault, the so-called Penacova-Régua-Verin megalineament. Along this tectonic structure different occurrences of CO₂-rich thermomineral waters are found: Chaves hot waters (67 °C) and also several cold (16.1 °C) CO₂-rich waters. The δ²H and δ¹⁸O values of the thermomineral waters are similar to those of the local meteoric waters. The chemical composition of both hot and cold mineral waters suggests that water–rock reactions are mainly controlled by the amount of dissolved CO₂ (g) rather than by the water temperature. Stable carbon isotope data indicate an external CO₂ inorganic origin for the gas. δ¹³C_CO2 values ranging between ? 7.2‰ and ? 5.1‰ are consistent with a two-component mixture between crustal and mantle-derived CO₂. Such an assumption is supported by the ³He/⁴He ratios measured in the gas phase, are between 0.89 and 2.68 times the atmospheric ratio (Ra). These ratios which are higher than that those expected for a pure crustal origin (≈ 0.02 Ra), indicating that 10 to 30% of the He has originated from the upper mantle. Release of deep-seated fluids having a mantle-derived component in a region without recent volcanic activity indicates that extensive neo-tectonic structures originating during the Alpine Orogeny are still active (i.e., the Chaves Depression).     

Mantle rheology: radial and lateral viscosity variations inferred from microphysical creep laws

Inferences on the rheology of the mantle based on theoretical and experimental rate equations for steady state creep are discussed and compared with results from geophysical models. The radial increase of viscosity by one to three orders of magnitude across the mantle, required by inversion of postglacial rebound and geodynamic data, is confirmed by microphysical models based on the estimation of continuous and discontinuous changes of creep parameters with depth. The upper mantle (viscosity 10²⁰–10²¹ Pa s) is likely to show non-Newtonian rheology (power-law creep) for average grain sizes larger than 0.1 mm as an order of magnitude. Given the variability of both grain size and stress conditions, local regions of linear rheology can be present. The rheology of transition zone and lower mantle (viscosity 10²²–10²⁴ Pa s) cannot be definitely resolved at present. Estimation of creep parameters leads to possible nonlinear or mixed rheology, if grain sizes are not lower than 0.1 mm and flow conditions can be approximated by a constant strain rate of about 10⁻¹⁵ s⁻¹. This conclusion can be modified by different flow conditions (e.g. a decrease in strain rate or constant viscous dissipation). Furthermore, experiments on fine-grained garnetites and perovskite analogues have shown that diffusion creep is predominant at laboratory conditions. However, the pressure dependence of creep in these phases is unknown, and therefore direct extrapolation to lower mantle conditions is necessarily speculative. Lateral variations of viscosity, largest in the upper and lowermost mantle (up to 2–4 orders of magnitude) are predicted by models based on lateral temperature anomalies derived from seismic tomographic models.