Palaeoenvironmental significance of the mammalian faunas of Italy since the Pliocene

The evolution of mammalian communities is a significant tool for reconstructing past environments and climate. In this paper, a palaeoecological and palaeoenvironmental reconstruction based on the study of the Italian mammal faunas ranging from Middle Pliocene to the Holocene is presented using the cenogram method and quantification of temperatures based on arvicolid species richness. These analyses reveal open and arid conditions during glacial periods and less open and more humid conditions during interglacials, with some differences between north and central–south Italy. Northern and Southern communities reflect similar trends in the evolution of palaeoenvironments and in the variation of the temperatures, but differences in estimated temperatures, which indicate a more temperate climate for southern Italy throughout the period, are considered here. Copyright © 2002 John Wiley & Sons, Ltd.     

Dysaerobic conditions during Heinrich events 4 and 5: Evidence from phosphorus distribution in a North Atlantic deep-sea core

Reactive phosphorus undergoes diagenetic transformation once transferred into marine sediments. The degree of regeneration and redistribution of phosphorus depends on early diagenetic and environmental conditions, which may be linked to larger scale phenomena, such as bottom water circulation, water column ventilation, and organic carbon flux. Phosphorus phases of the <50-μm-sized fraction of deep-sea sediments from core SU 90-09 (North Atlantic, 43°31′N, 30°24′W, 3375 m below sea level) have been analyzed using a sequential extraction technique (SEDEX method) to reconstruct phosphorus geochemistry during Heinrich events 4 and 5. Comparison with Holocene samples from the same site indicates that postdeposition diagenetic transformation has not affected phosphorus distribution in the deep part of the sediments. Total and reactive phosphorus average 0.40 ± 0.04 mg/g and 0.30 ± 0.05 mg/g, respectively, and are comparable to values found in analog deep-sea environments in the North Atlantic. Detrital phosphorus, the phase linked to igneous- and metamorphic-derived material, sharply increases during Heinrich events and covaries with the ice-rafted debris record, whereas authigenic and Fe-bound phosphorus phases, both influenced by redox conditions, decrease or even disappear. These findings suggest that during the deposition of Heinrich layers (HLs), environmental parameters hampered the precipitation of these phases. Large freshwater discharges in relation to iceberg surges may have provoked a temporary stratification of the water column. Accordingly, dysaerobic conditions in the sediments may have fostered the loss of dissolved phosphorus from the sediments to the water column, in a direct and rapid response to the changed conditions. Decreasing trends in organic matter elemental ratios (total organic carbon/organic phosphorus) and Rock-Eval oxygen index values, along with the presence of partly authigenic dolomite and ankerite within HLs, also support this assumption.     

Non-linear crustal corrections in high-resolution regional waveform seismic tomography

We compare 3-D upper mantle anisotropic structures beneath the North American continent obtained using standard and improved crustal corrections in the framework of Non-linear Asymptotic Coupling Theory (NACT) applied to long period three component fundamental and higher mode surface waveform data. Our improved approach to correct for crustal structure in high-resolution regional waveform tomographic models goes beyond the linear perturbation approximation, and is therefore more accurate in accounting for large variations in Moho topography within short distances as observed, for instance, at ocean–continent margins. This improved methodology decomposes the shallow-layer correction into a linear and non-linear part and makes use of 1-D sensitivity kernels defined according to local tectonic structure, both for the forward computation and for the computation of sensitivity kernels for inversion. The comparison of the 3-D upper mantle anisotropic structures derived using the standard and improved crustal correction approaches shows that the model norm is not strongly affected. However, significant variations are observed in the retrieved 3-D perturbations. The largest differences in the velocity models are present below 250 km depth and not in the uppermost mantle, as would be expected. We suggest that inaccurate crustal corrections preferentially map into the least constrained part of the model and therefore accurate corrections for shallow-layer structure are essential to improve our knowledge of parts of the upper mantle where our data have the smallest sensitivity.     

Origin and emplacement of ultramafic–mafic intrusions in the Erro-Tobbio mantle peridotite (Ligurian Alps, Italy)

The Erro-Tobbio peridotites (Voltri Massif, Ligurian Alps) represent subcontinental lithospheric mantle tectonically exhumed during Permo–Mesozoic extension of the Europe–Adria lithosphere. Previous studies have shown that exhumation started during Permian times, and occurred along kilometer-scale lithospheric shear zones which enhanced progressive deformation and recrystallization from spinel- to plagioclase-facies conditions. Ongoing field and petrologic investigations have revealed that the peridotites experienced, during uplift, a composite history of diffuse melt migration and multiple episodes of ultramafic–mafic intrusions. In this paper we present the results of field, structural and petrologic–geochemical investigations into a sector of the Erro-Tobbio peridotite unit that preserves well this multiple intrusion history. Melt impregnation in the peridotites is evidenced by significant plagioclase enrichment and crystallization of unstrained orthopyroxene replacing kinked mantle olivine and clinopyroxene; impregnating melts were thus opx-saturated. Melt–rock interaction caused chemical changes in mantle minerals (e.g. Al decrease and REE increase in cpx; Ti and Cr# enrichment in spinel). Nevertheless, clinopyroxenes still exhibit LREE depletion (Ce_N/Sm_N = 0.006–0.011), indicating a depleted signature for the percolating melts. Melt impregnation was thus related to diffuse porous flow migration of depleted MORB-type melt fractions that modified their compositions towards opx saturation by mantle–melt interaction during ascent. The impregnated peridotites are intruded by a hectometer-scale stratified cumulate body, mostly consisting of troctolites and plagioclase wehrlites, showing gradational, interfingered contacts with the host mantle rocks. Subsequent intrusion events are revealed by the occurrence of olivine gabbros as decameter-wide lenses, variably thick (centimeter- to meter-scale) dykes and thin dykelets, which crosscut both the peridotite foliation and the magmatic layering in the cumulates. Overall, major and trace element compositions of minerals in the intrusives indicate that they represent variably differentiated cumulus products crystallized from rather primitive N-MORB-type aggregated melts. Slightly more evolved compositions are shown by olivine gabbros, relative to the troctolites and plagioclase wehrlites of the cumulate body. Mineral chemistry features (e.g. the Fo–An correlation and high Na, Ti, Mg# in cpx) indicate that the studied intrusive rocks crystallized at moderate pressure conditions (3–5 kbar, i.e. 9–15 km depth). Our study thus points to a progressive transition from porous flow melt migration to emplacement of magmas in fractures, presumably related to progressive change of lithospheric mantle rheology during extension-related uplift and cooling.     

Seismic response of multi‐supported structures by proper orthogonal decomposition

The seismic analysis of structures is usually carried out considering the ground motion as fully‐correlated in space and determining the structural response by pseudo‐deterministic methods such as the response spectrum technique. Actually, the partial correlation of the seismic acceleration may influence heavily the behaviour of spatially extended structures, such as bridges, viaducts or pipelines. In order to take its partial correlation into account, the seismic ground motion is schematized as a stochastic process dependent on time and on space; the hypotheses of stationarity and homogeneity are used to obtain simple and general results. The influence of the partial correlation of the seismic ground motion on the structural response is investigated by introducing suitable Equivalent Spectra. The acceleration of the support‐points of the structure is represented by the Proper Orthogonal Decomposition (POD), defining the modes of the earthquake. The method is formulated for any kind of multi‐degree‐of‐freedom system and is applied, as a case study, to an ideal single‐storey multi‐supported frame with an axially rigid beam. In the case of two supports, the POD decouples the pseudo‐static and the dynamic contributions to the structural response. This property is preserved for structural systems with many supports, where only the lower modes of the earthquake, usually the first two POD modes, are responsible for the structural response. Copyright © 2003 John Wiley & Sons, Ltd.     

Monte Carlo study of the arrival time distribution of particles in extensive air showers in the energy range 1–100 TeV

A detailed simulation of vertical showers in atmosphere produced by primary gammas and protons, in the energy range 1–100 TeV, has been performed by means of the FLUKA Monte Carlo code, with the aim of studying the time structure of the shower front at different detector heights. It turns out that the time delay distribution can be fitted using few parameters coincident with the distribution central moments. Such parameters exhibit a smooth behaviour as a function of energy. These results can be used both for detector design and for the interpretation of the existing measurements. Differences in the time structure between gamma and proton induced showers are found and explained in terms of the nonrelativistic comonent of extensive air showers.     

A harmonic analysis of lunar topography

A global lunar topographic map has been derived from existing Earth-based and orbital observations supplemented in areas without data by a linear autocovariance predictor. Of 2592 bins, each 5° square, 1380 (64.7% by area) contain at least one measurement. A spherical harmonic analysis to degree 12 yields a mean radius of 1737.53 ± 0.03 km (formal standard error) and an offset of the center of figure of 1.98±0.06 km toward (19±2)°S, (194±1)°E. A Bouguer gravity map, derived from a 12-degree free-air gravity model and the present topography data, is presented for an elevation of 100 km above the mean surface. It is confirmed that the low-degree gravity harmonics are determined primarily by surface height variations and only secondarily by lateral density variations.     

Density distribution functions and correlations of earthquake parameters

The statistical analysis of the source parameters of 9 earthquake sets of different types (aftershocks, scattered events, swarms) and of different seismic regions shows that the density distribution function (ddf) of the linear dimension l of a fault is represented by a negative power law, as well as the ddf of the static stress drop σ and of the scalar seismic moment M_o. It is then suggested, and tentatively verified, that also the ddf of the root mean square ground acceleration, defined as a function of l and σ, may be represented by a negative power law and that, at least in the cases examined, it scales like the ddf of σ. It is seen that the variability of the static stress drop is significant from one region to another, as is well known, but it seems remarkable also in the same seismic region (in particular in California, σ varies by several orders of magnitude) and in the different sets of events of a given region (as observed again for California). It is hypothesized that a correlation, although weak, between the stress drop and the linear dimension of a fault exists and the analyses seem not to contradict that σ may be a decreasing function of l. Finally, it is suggested that the seismicity of a region may be represented two-dimensionally as a function of the ddf of the stress drop and of the linear dimension of a fault instead of the classic b and b_o values.