Diffuse CO<Subscript>2</Subscript> degassing at Vesuvio,Italy

At Vesuvio, a significant fraction of the rising hydrothermal–volcanic fluids is subjected to a condensation and separation process producing a CO₂–rich gas phase, mainly expulsed through soil diffuse degassing from well defined areas called diffuse degassing structures (DDS), and a liquid phase that flows towards the outer part of the volcanic cone. A large amount of thermal energy is associated with the steam condensation process and subsequent cooling of the liquid phase. The total amount of volcanic–hydrothermal CO₂ discharged through diffuse degassing has been computed through a sequential Gaussian simulation (sGs) approach based on several hundred accumulation chamber measurements and, at the time of the survey, amounted to 151 t d^–1. The steam associated with the CO₂ output, computed assuming that the original H₂O/CO₂ ratio of hydrothermal fluids is preserved in fumarolic effluents, is 553 t d^–1, and the energy produced by the steam condensation and cooling of the liquid phase is 1.47×10¹² J d^–1 (17 MW). The location of the CO₂ and temperature anomalies show that most of the gas is discharged from the inner part of the crater and suggests that crater morphology and local stratigraphy exert strong control on CO₂ degassing and subsurface steam condensation. The amounts of gas and energy released by Vesuvio are comparable to those released by other volcanic degassing areas of the world and their estimates, through periodic surveys of soil CO₂ flux, can constitute a useful tool to monitor volcanic activity.Editorial responsibility: H. Shinohara     

Benthic fluxes of mercury species in a lagoon environment (Grado Lagoon,Northern Adriatic Sea,Italy)

The role of the major biogeochemical processes in Hg cycling at the sediment–water interface was investigated in the Grado Lagoon (Northern Adriatic Sea). This wetland system has been extensively contaminated from the Idrija Hg Mine (Slovenia) through the Isonzo River suspended load carried by tidal fluxes. Three approaches were used to study the sediment–water exchange of total Hg (THg), methylmercury (MeHg), reactive Hg (RHg) and dissolved gaseous Hg (DGHg): (1) estimation of diffusive fluxes from porewater and overlying water concentrations, (2) measurements of benthic fluxes using a deployed light benthic chamber in situ and (3) measurements of benthic fluxes during oxic–anoxic transition with a laboratory incubation experiment. The THg solid phase, ranging between 9.5 and 14.4 μg g⁻¹, showed slight variability with depth and time. Conversely, MeHg contents were highest (up to 21.9 ng g⁻¹) at the surface; they tended to decrease to nearly zero concentration with depth, thus suggesting that MeHg production and accumulation occur predominantly just below the sediment–water interface. Porewater MeHg concentrations (0.9–7.9 ng L⁻¹, 0.15–15% of THg) varied seasonally; higher contents were observed in the warmer period. The MeHg diffusive fluxes (up to 17 ng m⁻² day⁻¹) were similar to those in the nearby Gulf of Trieste [Covelli, S., Horvat, M., Faganeli, J., Brambati, A., 1999. Porewater distribution and benthic flux of mercury and methylmercury in the Gulf of Trieste (Northern Adriatic Sea). Estuar. Coast. Shelf Sci. 48, 415–428], although the lagoon sediments contained four-fold higher THg concentrations. Conversely, the THg diffusive fluxes in the lagoon (up to 110 ng m⁻² day⁻¹) were one- to two-fold higher than those previously estimated for the Gulf of Trieste. The diurnal MeHg benthic fluxes were highest in summer at both sites (41,000 and 33,000 ng m⁻² day⁻¹ at the fishfarm and in the open lagoon, respectively), thus indicating the influence of temperature on microbial processes. The diurnal variations of dissolved THg and especially MeHg were positively correlated with O₂ and inversely with DIC, suggesting an important influence of benthic photosynthetic activities on lagoon benthic Hg cycling, possibly through the production of organic matter promptly available for methylation. The results from the dark chamber incubated in the laboratory showed that the regeneration of dissolved THg was slightly affected by the oxic–anoxic transition. Conversely, the benthic flux of MeHg was up to 15-fold higher in sediments overlain by O₂ depleted waters. In the anoxic phase, the MeHg fluxes proceeded in parallel with Fe fluxes and the methylated form reached approximately 100% of dissolved THg. The MeHg is mostly released into overlying water (mean recycling efficiency of 89%) until the occurrence of sulphide inhibition, due to scavenging of the available Hg substrate for methylation. The results suggest that sediments in the Grado Lagoon, especially during anoxic events, should be considered as a primary source of MeHg for the water column.     

The equations of relative motion in the orbital reference frame

The analysis of relative motion of two spacecraft in Earth-bound orbits is usually carried out on the basis of simplifying assumptions. In particular, the reference spacecraft is assumed to follow a circular orbit, in which case the equations of relative motion are governed by the well-known Hill–Clohessy–Wiltshire equations. Circular motion is not, however, a solution when the Earth’s flattening is accounted for, except for equatorial orbits, where in any case the acceleration term is not Newtonian. Several attempts have been made to account for the \(J_2\) effects, either by ingeniously taking advantage of their differential effects, or by cleverly introducing ad-hoc terms in the equations of motion on the basis of geometrical analysis of the \(J_2\) perturbing effects. Analysis of relative motion about an unperturbed elliptical orbit is the next step in complexity. Relative motion about a \(J_2\)-perturbed elliptic reference trajectory is clearly a challenging problem, which has received little attention. All these problems are based on either the Hill–Clohessy–Wiltshire equations for circular reference motion, or the de Vries/Tschauner–Hempel equations for elliptical reference motion, which are both approximate versions of the exact equations of relative motion. The main difference between the exact and approximate forms of these equations consists in the expression for the angular velocity and the angular acceleration of the rotating reference frame with respect to an inertial reference frame. The rotating reference frame is invariably taken as the local orbital frame, i.e., the RTN frame generated by the radial, the transverse, and the normal directions along the primary spacecraft orbit. Some authors have tried to account for the non-constant nature of the angular velocity vector, but have limited their correction to a mean motion value consistent with the \(J_2\) perturbation terms. However, the angular velocity vector is also affected in direction, which causes precession of the node and the argument of perigee, i.e., of the entire orbital plane. Here we provide a derivation of the exact equations of relative motion by expressing the angular velocity of the RTN frame in terms of the state vector of the reference spacecraft. As such, these equations are completely general, in the sense that the orbit of the reference spacecraft need only be known through its ephemeris, and therefore subject to any force field whatever. It is also shown that these equations reduce to either the Hill–Clohessy–Wiltshire, or the Tschauner–Hempel equations, depending on the level of approximation. The explicit form of the equations of relative motion with respect to a \(J_2\)-perturbed reference orbit is also introduced.     

The gravitational perturbation spectrum in linear satellite theory

A new method for calculating the perturbation spectrum in the framework of Kaula's linear satellite theory (LST) is introduced. The novelty of this approach consists in using recent results on the spectral decomposition of the perturbation frequencies in LST to provide a closed formulation for the amplitude and the phase of each line in the perturbation spectrum. The theory presented here can be applied to perturbations in the elements or in the radial and transverse directions due to the geopotential or to the tides. Separate algorithms are developed for application to orbits with circulating or frozen perigee.     

Kinematic evolution of thrust-related structures in the Umbro-Romagnan parautochthon (northern Apennines, Italy)

In the internal part of the Umbro-Marchean-Romagnan Apennines, the foredeep clastic wedge constituting the Neogene part of the sedimentary cover is completely detached from the underlying Mesozoic–Palaeogene succession. The resulting (Umbro-Romagnan) parautochthon consists of tectonostratigraphic units with a general geometry of broad synclinal blocks separated by narrow faulted anticlines.
Thrust-related structures observed in the field require thrust ramp propagation to have occurred within already folded rocks; therefore, they cannot be restored using simple fault-bend fold or fault-propagation folding models. Evidence for a passive fold origin in the studied rocks suggests that an early detachment folding episode preceded ramp propagation. The latter was facilitated by the enhanced thickness of incompetent material in the cores of detachment anticlines, which became the preferential sites where thrust ramps cut up-section. Depending on the trajectory of such thrust ramps, different types of fault-related structures could develop. Hanging-wall anticlines which give way to monoclinal structures higher up in the section are associated with listric thrust ramps, whereas hanging wall monoclines approximately parallel to the underlying fault surface are associated with straight-trajectory ramps.
This kinematic evolution, which occurred partly during syn-depositional compression, also accounts for the observed lithofacies distribution. The latter reflects an early differentiation of the foredeep trough into sub-basins that are progressively younger towards the foreland. The detachment anticlines that originally bounded such sub-basins were the site of later thrust propagation, leading to a tectonic juxtaposition of different tectonostratigraphic units consisting of broad NW-SE elongate synclinal blocks.     

Structural trends and palaeogeography of the central and western Sicily belt: new insights

The geology of the Sicilian mainland is summarized by N–S geological sections. A continuous late Cenozoic orogenic belt through central and western Sicily resulted from a complex deformative history, recorded by several tectonic events. The deformation mainly involved the sedimentary cover of the old African continental margin, formed in a large basinal area, bordered at its southern margin by a shallow-water carbonate environment attached to Gondwana. The orogenic belt involves a complex architecture of thrust systems, of different size, geometry and palaeogeographical origin. Deformation, which mainly developed in the earlier stages of thrusting in the basinal rock assemblages, mainly gave rise to a stack of three different duplex structures, respectively, composed of Palaeozoic, Mesozoic–Palaeogene and Neogene strata. Large-scale clockwise rotation of the thrusts predated transpressional movements in the hinterland during the latest Miocene to Pliocene. High- angle reverse faults, with lateral components, modified earlier tectonic contacts within the allochthons. Contemporaneous southwards- directed imbrications affected the external southern areas, progressively incorporating foreland and piggyback basirts. The stratigraphic relationships of basin-fills to the tectonic structures reveals that reactivation processes have been active during the last Plio-Pleistocene.     

A viscoelastic shear zone model of compressional and extensional plate boundaries

A model is proposed describing the mechanical evolution of a shear zone along compressional and extensional plate boundaries, subject to constant strain rate. The shear zones are assumed as viscoelastic with Maxwell rheology and with elastic and rheological parameters depending on temperature and petrology. Stress and strain are computed as functions of time and depth. For both kinds of boundaries the model reproduces the existence of a shallow seismogenic zone, characterized by a stress concentration. The thickness of the seismogenic layer is evaluated considering the variations of shear stress and frictional strength on faults embedded in the shear zone. Assuming that a fault dislocation takes place, the brittle-ductile transition is assumed to occur at the depth at which the time derivative of total shear stress changes from positive to negative values. The effects of different strain rates and geothermal gradients on the depth of the brittle-ductile transition are studied. The model predictions are consistent with values inferred from seismicity data of different boundary zones.     

Polyphase thrusting and dyke emplacement in the central Southern Alps (Northern Italy)

The Triassic succession of the central Southern Alps (Italy) is stacked into several units bounded by south-verging low-angle thrust faults, which are related to two successive steps of crustal shortening. The thrust surfaces are cut by high-angle extensional and strike-slip faults, which controlled the emplacement of hypabissal magmatic intrusions that post-date thrusts motions. Intrusion ages based on SHRIMP U–Pb zircon dating span between 42 ± 1 and 39 ± 1 Ma, suggesting close time relationships with the earliest Adamello intrusion stages and, more in general, with the widespread calc-alkaline magmatism described in the Southern Alps. Fission-track ages of magmatic apatites are indistinguishable from U–Pb crystallization ages of zircons, suggesting that the intrusion occurred in country rocks already exhumed above the partial annealing zone of apatite (depth < 2–4 km). These data indicate that the central Southern Alps were already structured and largely exhumed in the Middle Eocene. Although we describe minor faults affecting magmatic bodies and local reactivations of older structures, no major internal deformations have occurred in the area after the Bartonian. Neogene deformations were instead concentrated farther south, along the frontal part of the belt.     

Numerical SPH analysis of debris flow run-out and related river damming scenarios for a local case study in SW China

A smoothed particle hydrodynamics (SPH) numerical modeling method implemented for the forward simulation of propagation and deposition of flow-type landslides was combined with different empirical geomorphological index approaches for the assessment of the formation of landslide dams and their possible evolution for a local case study in southwestern China. The SPH model was calibrated with a previously occurred landslide that formed a stable dam impounding the main river, and it enabled the simulation of final landslide volumes, and the spatial distribution of the resulting landslide deposits. At four different sites on the endangered slope, landslides of three different volumes were simulated, respectively. All landslides deposited in the main river, bearing the potential for either stable impoundment of the river and upstream flooding scenarios, or sudden breach of incompletely formed or unstable landslide dams and possible outburst floods downstream. With the empirical indices, none of the cases could be identified as stable formed landslide dam when considering thresholds reported in the literature, showing up the limitations of these indices for particular case studies of small or intermediate landslide volumes and the necessity to adapt thresholds accordingly for particular regions or sites. Using the occurred benchmark landslide as a reference, two cases could be identified where a complete blockage occurs that is more stable than the reference case. The other cases where a complete blockage was simulated can be considered as potential dam-breach scenarios.