Sediments and pollution in the Northern Adriatic Sea: a statistical analysis

Data concerning grain-size distribution, mineralogy, geochemistry, and chlorinated hydrocarbon concentrations, obtained from 246 surface sediment samples collected in the Northern Adriatic sea (Italian area), provide the basis for a multivariate statistical analysis which evaluates the influence of riverine inputs, the differences among various depositional environments and the presence of polluted zones.Cluster analysis applied to grain size has identified nine groups of lithotypes associated with the morphology and the hydrology of the basin.In the Gulf of Venice, principal component analysis on clay minerals and carbonates has permitted recognition of two main areas located north and south of the Brenta river. More subareas correspond to other lithological types.Principal component analysis of nine metal (Hg, Cd, Cu, Pb, Cr, Zn, Fe, Ni, and Co) concentrations reveals that the studied sea area is more complex from the geochemical point of view. Three geochemical zones related to riverine inputs, different lithotypes and specific forms of contamination are identified. High pollution of both halocarbons and heavy metals is associated with fine grain size.     

Multiphase model for inclusion‐reinforced geostructures: Application to rock‐bolted tunnels and piled raft foundations

A multiphase model is proposed to describe the mechanical behaviour of geomaterials reinforced by linear inclusions. This macroscopic approach considers the reinforced soil or rock mass as the superposition of continuous media. Equations of motion and constitutive laws of the model are first derived. Its implementation in a finite element computer code is then detailed. A modified implicit algorithm for elastoplastic problems is proposed. The model and its implementation are fully validated for rock‐bolted tunnels (comparison with scale model experiments) and piled raft foundations (comparison with the classical ‘hybrid method’). The Messeturm case history is finally presented to assess the handiness of the approach for real structures. Copyright © 2001 John Wiley & Sons, Ltd.     

High-pressure,low-temperature metamorphism in Alpujarride Units of southeastern Betics (Spain)Métamorphisme de haute pression et basse température dans les unités Alpujarrides du Sud-Est de la cordillère Bétique (Espagne)

We present the first occurrences of high-pressure, low-temperature ferro-magnesiocarpholite-bearing mineral assemblages associated to quartz segregations in the Alpujarride units of southeastern Betics (Sierra de Almagro, Sierra de los Pinos and Sierra Cabrera). Thermobarometric results show that the carpholite-bearing rocks underwent the same P–T conditions in the three outcrops, i.e. 8–10 kbar, 350–400 °C. Metamorphic and structural data allow us to conclude that these rocks belong to the same Alpujarride unit. In the Sierra de Almagro, tectonic units with carpholite-bearing rocks overlie low-pressure, low-temperature Alpujarride units, then forming a stack with an inverted tectono-metamorphic sequence, as observed in the central and western part of the Alpujarride complex. The preservation of carpholite-bearing assemblages in these rocks implies that no significant temperature increase occurred during the exhumation history. To cite this article: G. Booth-Rea et al., C. R. Geoscience 334 (2002) 857–865.     

Image analysis and circular statistics for shape-fabric analysis: applications to lithified ignimbrites

 Computer-assisted image analysis can be successfully used to derive quantitative textural data on pyroclastic rock samples. This method provides a large number of different measurements such as grain size, particle shape and 2D orientation of particle main axes (directional- or shape-fabric) automatically and in a relatively short time. Orientation data reduction requires specific statistical tests, mainly devoted to defining the kind of particle distribution pattern, the possible occurrence of preferred particle orientation, the confidence interval of the mean direction and the degree of randomness with respect to pre-assigned theoretical frequency distributions. Data obtained from image analysis of seven lithified ignimbrite samples from the Vulsini Volcanic District (Central Italy) are used to test different statistics and to provide insight about directional fabrics. First, the possible occurrence of a significant deviation from a theoretical circular uniform distribution was evaluated by using the Rayleigh and Tukey χ ² tests. Then, the Kuiper test was performed to evaluate whether or not the observation fits with a unimodal, Von Mises-like theoretical frequency distribution. Finally, the confidence interval of mean direction was calculated. With the exception of one sample (FPD10), which showed a well-developed bimodality, all the analysed samples display significant anisotropic and unimodal distributions. The minimum number of measurements necessary to obtain reasonable variabilities of the calculated statistics and mean directions was evaluated by repeating random collections of the measured particles at increments of 100 particles for each sample. Although the observed variabilities depend largely on the pattern of distribution and an absolute minimum number cannot be stated, approximately 1500–2000 measurements are required in order to get meaningful mean directions for the analysed samples. Received: 9 April 1996 / Accepted: 26 December 1996     

Foldbelts with early salt withdrawal and diapirism: Physical model and examples from the northern Gulf of Mexico and the Flinders Ranges, Australia

A physical experiment shows that shortening applied to existing diapirs and minibasins produces anomalous structural styles that are unlike those of more typical foldbelts. Strong minibasins remain largely undeformed while weak diapirs localize contractional strain. Short diapirs form the cores to folds and thrusted folds, whereas tall diapirs are squeezed and often welded, commonly leading to the extrusion of allochthonous material. Key features of the model are observed in real examples. In the northern Gulf of Mexico passive margin, minibasins were originally separated by a polygonal pattern of deep salt ridges, with diapirs located at ridge intersections. Gravity spreading resulted in squeezed diapirs (and associated allochthonous salt) connected by variably oriented contractional, extensional, and strike-slip structures. In the Flinders Ranges convergent-margin foldbelt of South Australia, preexisting diapirs were squeezed, welded, and thrusted, with anticlines plunging away in multiple directions, so that minibasins are surrounded by highly variable structures. A different geometry is observed in La Popa Basin, Mexico, where squeezing of a linear salt wall produced a vertical weld with diapirs at the terminations, rather than the culmination. In all areas, foldbelt geometries are strongly influenced by the preestablished salt-minibasin architecture.     

The Limit Problems for the Equation of Oscillations of a Satellite

We consider the ordinary differential equation of the second order, which describes oscillations of a satellite with respect to its mass center moving along an elliptic orbit with eccentricity e. The equation has two parameters: e and μ. It is regular for 0 ≤ e < 1 and singular when e = 1. For 1 we obtain three limit problems. Their bounded solution to the first limit problem form a two-dimensional (2D) continuous invariant set with a periodic structure. Solutions to the second limit problem form 2D and 3D manifolds. The μ-depending families of odd bounded solutions are singled out. One of the families is twisted into a self-similar spiral. To obtain the limit families of the periodic solutions to the original problem match together the odd bounded solutions to the first and the second limit problem. The point of conjunction is described by the third (the basic) limit problem. The limit families are very close to prelimit ones computed in earlier studies. This revised version was published online in July 2006 with corrections to the Cover Date.     

A review of current issues in air pollution modeling and simulation

Air pollution modeling is now a mature field, and comprehensive numerical models (the chemistry-transport models) are used in many applications. This article aims at reviewing the main issues from the point of view of applied mathematics and computational physics (as viewed by the author). We address topics such as subgrid parameterization, numerical algorithms with a focus on aerosol simulation, data assimilation and inverse modeling, reduction of high-dimensional models and propagation of uncertainties. Even if this article is strictly related to air pollution modeling, many issues and methods can be extended to dispersion of tracers in other media (for instance, water). This article is based on a plenary talk given at SIAM Geosciences, Avignon, June 2005.     

Comparison of several BHT correction methods: a case study on an Australian data set

Bottom-hole temperatures (BHT) from oil exploration provide useful constraints on the subsurface thermal regime, but they need to be corrected to obtain the equilibrium temperature. In this work, we introduce several BHT correction methods and compare them using a large Australian data set of more than 650 groups of multiple BHT measurements in about 300 oil exploration boreholes. Existing and suggested corrections are classified within a coherent framework, in which methods are divided into: line/cylinder source; instantaneous/continuous heat extraction; one/two component(s). Comparisons with reservoir test temperatures show that most of the corrections lead to reliable estimates of the formation equilibrium temperature within ±10 °C, but too few data exist to perform an intercomparison of the models based on this criterion. As expected, the Horner method diverges from its parent models for small elapsed times (or equivalently large radii). The mathematical expression of line source models suffers from an unphysical delay time that also restrains their domain of applicability. The model that takes into account the difference of thermal properties between circulating mud and surrounding rocks—that is the two-component model—is delicate to use because of its high complexity. For these reasons, our preferred correction methods are the cylindrical source models. We show that mud circulation time below 10 hr has a negligible effect. The cylindrical source models rely on one parameter depending on the thermal diffusivity and the borehole radius, which are poorly constrained, but the induced uncertainty on the extrapolations remains reasonably low.     

Simulating the mobile barrier (MOSE) operation in the Venice Lagoon, Italy: global sea level rise and its implication for navigation

Three years of mobile barrier operations have been simulated with a hydrodynamic model to check the efficiency of the barriers in defending the city of Venice from flooding. The simulations have been carried out in the actual situation and with a sea-level rise of 30 and 50 cm. Moreover, the interference of the barrier operations with the ship traffic has been studied. It is found that without a security increment for the forecasted water levels, the mobile barriers cannot defend completely Venice from flooding due to the uncertainty in the forecast. With a security increment of 10 cm, the barriers work well in actual conditions but still cannot avoid flooding with a global sea-level rise. The interference with the ship traffic is acceptable under actual conditions but becomes prohibitive with a sea-level rise of 50 cm, when nearly two-thirds of the ship passages are blocked or delayed.