Clinopyroxene/liquid trace element partitioning in natural trachyte–trachyphonolite systems: insights from Campi Flegrei (southern Italy)

Trace element partition coefficients between clinopyroxenes and associated glassy matrix (^Cpx/L D) have been determined for 13 REE, HFSE^4+,5+, U, Th, Sr, Pb, Sc and V from combined LA-ICP-MS/EMP analyses in selected trachytes and trachyphonolites from Campi Flegrei. Composition of clinopyroxene and glass is pretty homogeneous in the trachyphonolites, pointing to an overall attainment of the equilibrium conditions. In trachytes, conversely, phases show some compositional heterogeneity (due to the presence of clinopyroxene xenocrysts) that requested a more careful petrographic and geochemical inspection of the samples to assess the equilibrium clinopyroxene composition. In the trachyte clinopyroxenes, REE are compatible from Nd to Lu (^Cpx/L D up to 2.9), like Y, Ti, Sc and V. The ^Cpx/L D for Eu is lower than those of the adjacent REE, highlighting Eu²⁺ contribution. High D values are also shown by U, Th, Pb, Zr, Hf, Nb and Ta relatively to basaltic and andesitic systems, whereas the D _Sr is roughly similar to that found for less evolved magmas. Trachyphonolites are characterized by an overall decrease of the ^Cpx/L D for highly-charged cations (with the exception of V), and by a slight increase of D _Sr. REE are still compatible from Nd to Lu (^Cpx/L D up to 2.1), like Ti, Y, Sc and V. This variation is also predicted for REE and Y by models based on the elastic strain theory, being consistent with the slightly lower polymerization degree estimated for the trachyphonolites. However, the observed ^Cpx/L D _(REE,Y) are matched by the modelled ones only considering very low T (≤825°C), which are believed unlikely. This mismatch cannot be attributed to effects induced by the water-rich composition of the trachyte–trachyphonolite suite, since they would lower the observed ^Cpx/L D _(REE,Y). Moreover, the anomalous inflections of measured ^Cpx/L D for HREE suggests some crystal-chemical control, such as the entrance of these elements in a site distinct from M2. It is concluded that the large ^Cpx/L D determined for trachytes and trachyphonolites are likely induced by hitherto unconstrained changes of the Z³⁺ activities related to the composition of melt and/or solid. All these considerations strongly highlight the importance of a direct characterization of trace element partitioning in natural samples from magmatic systems poorly characterized by experimental studies. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Analysis of Titan CH4 3.3 μm upper atmospheric emission as measured by Cassini/VIMS

After molecular nitrogen, methane is the most abundant species in Titan’s atmosphere and plays a major role in its energy budget and its chemistry. Methane has strong bands at 3.3 μm emitting mainly at daytime after absorption of solar radiation. This emission is strongly affected by non-local thermodynamic equilibrium (non-LTE) in Titan’s upper atmosphere and, hence, an accurate modeling of the non-LTE populations of the emitting vibrational levels is necessary for its analysis. We present a sophisticated and extensive non-LTE model which considers 22 CH₄ levels and takes into account all known excitation mechanisms in which they take part. Solar absorption is the major excitation process controlling the population of the v₃-quanta levels above 1000 km whereas the distribution of the vibrational energy within levels of similar energy through collisions with N₂ is also of importance below that altitude. CH₄-CH₄ vibrational exchange of v₄-quanta affects their population below 500 km. We found that the ν₃ → ground band dominates Titan’s 3.3 μm daytime limb radiance above 750 km whereas the ν₃ + ν₄ → ν₄ band does below that altitude and down to 300 km. The ν₃ + ν₂ → ν₂, the 2ν₃ → ν₃, and the ¹³CH₄ν₃ → ground bands each contribute from 5% to 8% at regions below 800 km. The ν₃ + 2ν₄ → 2ν₄and ν₂ + ν₃ + ν₄ → ν₂ + ν₄ bands each contribute from 2% to 5% below 650 km. Contributions from other CH₄ bands are negligible. We have used the non-LTE model to retrieve the CH₄ abundance from 500 to 1100 km in the southern hemisphere from Cassini-VIMS daytime measurements near 3.3 μm. Our retrievals show good agreement with previous measurements and model results, supporting a weak deviation from well mixed values from the lower atmosphere up to 1000 km.     

A site amplification factor map of the Iberian Peninsula and the Balearic Islands

Considering the lack of site effect cartography in the Iberian Peninsula region to be related to the strong motion data, we present a site amplification factor map covering the Iberian Peninsula and the Balearic Islands. Steps leading to produce the map have included the classification of geological units from the 1:1,000,000 scale Geological Map of Spain into six site classes characterized on a seismic response basis, and the calculation, for every site class, of both short-period and mid-period amplification factors. In order to test the validity and applicability of the map, we have calculated the synthetic seismic intensity in the particular localities where several past earthquakes were felt, supposing a point source approximation. The synthetic intensities have been obtained with and without site amplification factors, resulting in two types of synthetic isoseismals maps that have been compared with the corresponding observed isoseismals. As the amplification factors obtained are only applicable to the linear domain, the synthetic intensities greater than VII are only illustrative. A main conclusion has been drawn about the fact that synthetic isoseismals with site amplification factor approach to the related observed isoseismals in a higher degree than synthetic isoseismals without site amplification factor. In addition, the resemblances between the synthetic isoseismals produced with site amplification factors and the corresponding observed isoseismals have been found to be more evident at shorter epicentral distances, provided that intensities are lower than VIII, and therefore, we remain in the linear domain.     

Asteroid 2002NY40 as a source of meteorite-dropping bolides

The existence of asteroidal meteoroid streams capable of producing meteorite-dropping bolides has long being invoked, but evidence is scarce. Recent modelling of previously reported associations suggests that the time-scales to keep the orbital coherence of these streams producing meteorites are too short. We present an unequivocal association between near earth object (NEO) 2002NY40 and at least one bright fireball detected over Finland in 2006 August. Another two additional fireballs recorded from Spain and Finland seem to be related, together producing a fireball-producing stream (β Aquarids). On the basis of historical data, the 2006 finding suggests the existence of a meteoroid complex capable of producing meteorites. Taking into account present time-scales for orbital decoherence, if 2002NY40 has large meteoroids associated with it, such behaviour would be the consequence of a relatively recent asteroidal fragmentation. Supporting our claim, the heliocentric orbits of two recently discovered NEOs, 2004NL8 and 2002NY40, were found to exhibit a good similarity to each other and also to the orbits of the three bolides. The fireball spectra of the two Finish bolides showed that the chemical abundances of these objects are consistent with the main elements found in chondrites. This result is consistent with the probable Low iron, Low metal (LL) chondritic mineralogy of asteroid 2002NY40. Consequently, this asteroid may be delivering LL chondrites to the Earth. Additional fireball reports found in the literature suggest that the associated β Aquarid complex may have been delivering meteorites to the Earth during, at least, the last millennium.     

Recent evolution of surge-related events and assessment of coastal flooding risk on the eastern coasts of the English Channel

This paper is based on statistical analysis of hourly tide measurements for some 285 equivalent full years from the stations of Weymouth, Bournemouth, Portsmouth, Newhaven, Dover and Sheerness in the UK, and of Cherbourg, Le Havre, Dieppe, Boulogne, Calais and Dunkirk in France. For each tidal value, surge heights have been determined and correlated with hourly or three-hourly wind and air pressure data from nearby meteorological stations. Major surges in the area are generally produced by storms associated with wind from north-west or south-west that tend to push oceanic water into the Channel. Recent medium-term climate evolution does not seem to increase the flooding risk at French stations, where surge-related winds tend to decrease in frequency and speed (Cherbourg, Dieppe and Boulogne) or show little change (Le Havre). However, the long-term risk of flooding will increase through the loss in land elevation due to a continuation of the local relative sea-level rise, especially if this effect will be enhanced by an acceleration in the global sea-level rise predicted by climatic models. The northern side of the Channel (Weymouth, Bournemouth and Portsmouth) is mainly exposed to southerly winds that show variable trends. It is also apparently affected by strong subsidence trends during the last two decades. If lasting, such trends can only increase long-term flooding risk. The flooding risk has not increased near the eastern end of the Channel. The duration of significant cyclonic events tends to decrease near Cherbourg but tends to increase near Weymouth, with no conclusive trends in other stations (Portsmouth, Calais and Dunkirk), where extreme surges may occur also in relatively high-air-pressure situations. In conclusion, medium-term coastal flooding risk seems to increase especially at Weymouth, Bournemouth and Portsmouth, and also, but less so, at Le Havre and Sheerness. In addition, few extreme surges occurred during the last decades at the time of spring high tide, which would seem to be a fortunate coincidence or, in some cases, an effect of tide–surge interaction. The risk of occurrence of less favourable random events in the near future is therefore of concern, and flood potential would greatly increase if the global sea-level rise expected in the near future is also considered.     

Can beach dune ridges of the Texas Gulf Coast preserve climate signals?

A study of the evolution of North Padre Island (southern Texas Gulf Coast) dunes was carried out using LIDAR topographic data, dune vibracores through the center of the dunes, and grab samples of shoreface sand at four locations along a cross-shore profile. Grain-size analyses of the vibracores show vertical variations in shoreface sand deposition over decimeter depth intervals. A dune ridge growth model is introduced that describes the dune vertical accretion rate as a function of island progradation and freshwater lens expansion. This model allows indirect dating of the dune core samples based on a known island progradation rate (1?m/year), and height and spacing of the dunes calculated from the topographic data. A sand provenance model is also proposed that links the sand deposition in the dunes with sand sourced from various depths along the shoreface profile, depending on storm activity. We present evidence linking the changes in storm-sand deposition in the dune cores with yearly climatic fluctuations in the Gulf of Mexico associated with landfalling tropical storm activity in the period from 1942?C1965. This record of storm-induced sand variability is negatively correlated with El Ni?o-Southern Oscillation (Pacific) sea surface temperature variability, and positively correlated with North Atlantic decadal sea surface temperature variability.     

A physical analogue model to analyse interactions between tensile stresses and dissolution in carbonate slopes

A physical analogue model was developed to analyse the relationship between tensional states in rock masses, seepage and karst processes. Use was made of an experimental apparatus consisting of two hydraulic circuits realized by drilling two holes into each of two blocks of sampled limestone from the central Apennines (Italy). A static load (208.85 kg) was applied to one of the blocks in order to elicit tensile stresses within it. Physical and chemical monitoring data showed that the main process involved was temperature-dependent CaCO₃ dissolution. This process was more marked in the loaded block circuit, as Ca⁺⁺ concentration in circulated water reached 54 mg/L, whereas only 28 mg/L was reached in the unloaded one. The interaction between load and dissolution caused the observed opening of microcracks, as confirmed by further increase of water loss and by dilution in the loaded block circuit, resulting in a decrease of Ca⁺⁺ concentration. These findings were confirmed by recording additional water losses after increasing the load up to 445.05 kg. A finite difference numerical model showed that tensile stresses (max 20 kPa) within the loaded block were clustered at the intersection of the main joints with the flowpaths, thus representing points of preferential and accelerated dissolution.This study is part of a research project on analogical and numerical modelling for geological risk mitigation being conducted at CERI, Research Centre for Geological Risks of the University of Rome “La Sapienza” (Valmontone, Roma, Italy; http://www.ceri@uniroma1.it).     

Analytical expressions for three-phase generalized relative permeabilities in water- and oil-wet capillary tubes

We analyze three-phase flow of immiscible fluids taking place within an elementary capillary tube with circular cross-section under water- and oil-wet conditions. We account explicitly for momentum transfer between the moving phases, which leads to the phenomenon of viscous coupling, by imposing continuity of velocity and shear stress at fluid-fluid interfaces. The macroscopic flow model which describes the system at the Darcy scale includes three-phase effective relative permeabilities, K _{i j,r} , accounting for the flux of the ith phase due to the presence of the jth phase. These effective parameters strongly depend on phase saturations, fluid viscosities, and wettability of the solid matrix. In the considered flow setting, K _{i j,r} reduce to a set of nine scalar quantities, K _{i j,r} . Our results show that K _{i j,r} of the wetting phase is a function only of the fluid phase own saturation. Otherwise, K _{i j,r} of the non-wetting phase depends on the saturation of all fluids in the system and on oil and water viscosities. Viscous coupling effects (encapsulated in K _{i j,r} with i ≠ j) can be significantly relevant in both water- and oil-wet systems. Wettability conditions influence oil flow at a rate that increases linearly with viscosity ratio between oil and water phases.