Garnet and spinel in fertile and depleted mantle: insights from thermodynamic modelling

We performed thermodynamic calculations based on model and natural peridotitic compositions at pressure and temperature conditions relevant to the Earth’s upper mantle, using well-established free energy minimization techniques. The model is consistent with the available experimental data in Cr-bearing peridotitic systems and can therefore be used to predict phase relations and mineral compositions in a wide range of realistic mantle compositions. The generated phase diagrams for six different bulk compositions, representative of fertile, depleted and ultra-depleted peridotitic mantle, shown that the garnet + spinel stability field is always broad at low temperatures and progressively narrows with increasing temperatures. In lithospheric sections with hot geotherms (ca. 60 mW/m²), garnet coexists with spinel across an interval of 10–15 km, at ca. 50–70 km depths. In colder, cratonic, lithospheric sections (e.g. along a 40 mW/m² geotherm), the width of the garnet–spinel transition strongly depends on bulk composition: In fertile mantle, spinel can coexist with garnet to about 120 km depth, while in an ultra-depleted harzburgitic mantle, it can be stable to over 180 km depth. The formation of chromian spinel inclusions in diamonds is restricted to pressures between 4.0 and 6.0 GPa. The modes of spinel decrease rapidly to less than 1 vol % when it coexists with garnet; hence, spinel grains can be easily overlooked during the petrographical characterization of small mantle xenoliths. The very Cr-rich nature of many spinels from xenoliths and diamonds from cratonic settings may be simply a consequence of their low modes in high-pressure assemblages; thus, their composition does not necessarily imply an extremely refractory composition of the source rock. The model also shows that large Ca and Cr variations in lherzolitic garnets in equilibrium with spinel can be explained by variations of pressure and temperature along a continental geotherm and do not necessarily imply variations of bulk composition. The slope of the Cr# [i.e. Cr/(Cr + Al)mol] isopleths in garnet in equilibrium with spinel changes significantly at high temperatures, posing serious limitations to the applicability of empirical geobarometric methods calibrated on cratonic mantle xenoliths in hotter, off-craton, lithospheric mantle sections.     

Seismic study of the Mesozoic carbonate basement around Mt. Somma–Vesuvius, Italy

Fifteen seismic reflection lines from AGIP surveys, in and around the Campanian Plain and Mt. Somma–Vesuvius (south Italy) have been interpreted. The attention has been focused to the horizon pertinent to the top of the Mesozoic carbonate sequence and the Quaternary faults dissecting it. As a matter of fact, both are very important elements for understanding the origin of the volcanic activity in the area, that often in the past, has been the topic of debates not supported by reliable data. In the study area, referring to the depth of the carbonate basement, comparison between the result achieved by the seismic prospecting and previous gravity studies has been made. It shows coherence in some areas but large discrepancy within others. Near the town of S. Anastasia, the gravity and seismic depth estimates differ as much as 1000 m or more. Furthermore, the seismic data show that the source of the greatest volcanic eruption in the area (the so-called ‘Campanian Ignimbrite') is probably not located in the Acerra depression, as suggested by other authors. A main NE–SW fault directed toward Vesuvius, considered as playing a primary role on volcanogenetic processes and previously recognised only offshore by marine seismic survey, has been now identified also inland using this new seismic information. The results presented here strengthen the hypothesis that Mt. Vesuvius is located at the crossing point of two regional Quaternary sets of fault heading NW–SE and NE–SW.     

Effect of the dilution process on the attenuation of contaminants in aquifers

This paper discusses the effects of the dilution process on levels of aquifer contamination attenuation, developing analytical equations for application in theoretical and practical cases. The key finding is that, in aquifers, the dilution process causes a reduction of the contaminant concentration. Analytical equations for the assessment of the contaminant attenuation capacity, by means of the dilution process in an aquifer, are put forward. In particular, the attenuation of the contaminant concentration in an aquifer, and thus the dilution, is described by the volumetric flow rate per unit width perpendicular to the flow direction (q _u). Moreover, a classification of q _u was developed, considering six classes of attenuation capacity due to dilution. For the equations validation, nitrate contamination attenuation due to dilution was analysed in two sample areas in Piedmont (North-West Italy). The results confirmed the hypothesis, evidencing high nitrate concentrations in areas with low q _u values, connected to low dilution rate of the contaminant in the groundwater; vice versa, low nitrate concentrations were evidenced in areas with high q _u values, due to a high nitrate dilution rate.     

Do static stress changes of a moderate‐magnitude earthquake significantly modify the regional seismic hazard? Hints from the L'Aquila 2009 normal‐faulting earthquake (Mw 6.3, central Italy)

We investigated the Coulomb stress changes in the active faults surrounding a moderate‐magnitude normal‐faulting earthquake (2009 L'Aquila, Mw 6.3) and the associated variations in the expected ground motion on regional probabilistic seismic hazard maps. We show that the static stress variations can locally increase the seismic hazard by modifying the expected mean recurrence time on neighbouring faults by up to ~290 years, with associated variations in the probability of occurrence of the maximum expected earthquake of up to ~2%. Our findings suggest that the increase in seismic hazard on neighbouring faults following moderate‐magnitude earthquakes is probably not sufficient to necessitate systematic upgrades of regional probabilistic seismic hazard maps, but must be considered to better address and schedule strategies for local‐scale mitigation of seismic risk.     

Dynamical evolution of NEAs: Close encounters, secular perturbations and resonances

We discuss the main mechanisms affecting the dynamical evolution of Near-Earth Asteroids (NEAs) by analyzing the results of three numerical integrations over 1 Myr of the NEA (4179) Toutatis. In the first integration the only perturbing planet is the Earth. So the evolution is dominated by close encounters and looks like a random walk in semimajor axis and a correlated random walk in eccentricity, keeping almost constant the perihelion distance and the Tisserand invariant. In the second integration Jupiter and Saturn are present instead of the Earth, and the 3/1 (mean motion) and v ₆ (secular) resonances substantially change the eccentricity but not the semimajor axis. The third, most realistic, integration including all the three planets together shows a complex interplay of effects, with close encounters switching the orbit between different resonant states and no approximate conservation of the Tisserand invariant. This shows that simplified 3-body or 4-body models cannot be used to predict the typical evolution patterns and time scales of NEAs, and in particular that resonances provide some fast-track dynamical routes from low-eccentricity to very eccentric, planet-crossing orbits.On leave from the Department of Mathematics, University of Pisa, Via Buonarroti 2, 56127 Pisa, Italy, thanks to the G. Colombo fellowships of the European Space Agency.     

New model to determine forces at on-bottom slender pipelines

The present paper proposes a numerical model to determine horizontal and vertical components of the hydrodynamic forces on a slender submarine pipeline lying at the sea bed and exposed to non-linear waves plus a current. The new model is an extension of the Wake II type model, originally proposed for sinusoidal waves (Soedigdo et al., 1999) and for combined sinusoidal waves and currents (Sabag et al., 2000), to the case of periodic or random waves, even with a superimposed current. The Wake II type model takes into account the wake effects on the kinematic field and the time variation of drag and lift hydrodynamic coefficients. The proposed extension is based on an evolutional analysis carried out for each half period of the free stream horizontal velocity at the pipeline. An analytical expression of the wake velocity is developed starting from the Navier–Stokes and the boundary layer equations. The time variation of the drag and lift hydrodynamic coefficients is obtained using a Gaussian integration of the start-up function. A reduced scale laboratory investigation in a large wave flume has been conducted in order to calibrate the empirical parameters involved in the proposed model. Different wave and current conditions have been considered and measurements of free stream horizontal velocities and dynamic pressures on a bottom-mounted pipeline have been conducted. The comparison between experimental and numerical hydrodynamic forces shows the accuracy of the new model in evaluating the time variation of peaks and phase shifts of the horizontal and vertical wave and current induced forces.