The double solid reactant method: II. An application to the shallow groundwaters of the Porto Plain,Vulcano Island (Italy)

This paper presents an example of application of the double solid reactant method (DSRM) of Accornero and Marini (Environmental Geology, 2007a), an effective way for modeling the fate of several dissolved trace elements during water–rock interaction. The EQ3/6 software package was used for simulating the irreversible water–rock mass transfer accompanying the generation of the groundwaters of the Porto Plain shallow aquifer, starting from a degassed diluted crateric steam condensate. Reaction path modeling was performed in reaction progress mode and under closed-system conditions. The simulations assumed: (1) bulk dissolution (i.e., without any constraint on the kinetics of dissolution/precipitation reactions) of a single solid phase, a leucite-latitic glass, and (2) precipitation of amorphous silica, barite, alunite, jarosite, anhydrite, kaolinite, a solid mixture of smectites, fluorite, a solid mixture of hydroxides, illite-K, a solid mixture of saponites, a solid mixture of trigonal carbonates and a solid mixture of orthorhombic carbonates. Analytical concentrations of major chemical elements and several trace elements (Cr, Mn, Fe, Ni, Cu, Zn, As, Sr and Ba) in groundwaters were satisfactorily reproduced. In addition to these simulations, similar runs for a rhyolite, a latite and a trachyte permitted to calculate major oxide contents for the authigenic paragenesis which are comparable, to a first approximation, with the corresponding data measured for local altered rocks belonging to the silicic, advanced argillic and intermediate argillic alteration facies. The important role played by both the solid mixture of trigonal carbonates as sequestrator of Mn, Zn, Cu and Ni and the solid mixture of orthorhombic carbonates as scavenger of Sr and Ba is emphasized.

A novel methodological approach for land subsidence prediction through data assimilation techniques

Computational Geosciences - Anthropogenic land subsidence can be evaluated and predicted by numerical models, which are often built over deterministic analyses. However, uncertainties and...     

Chronology of the Mediterranean sea-level highstand during the Last Interglacial: a critical review of the U/Th-dated deposits

Relative sea-level (RSL) evolution during Marine Isotopic Stage (MIS) 5 in the Mediterranean basin is still not fully understood despite a plethora of morphological, stratigraphic and geochronological studies carried out on highstand deposits of this area. In this review we assembled a database of 323 U/Th-dated samples (e.g. corals, molluscs, speleothems) which were used to chronologically constrain RSL evolution within MIS 5. The application of strict geochemical criteria to the U/Th samples indicates that only ~33% of data available for the Mediterranean Sea can be considered ‘reliable’. Most of these data (~65%) refer to the MIS 5e highstand, while only ~17% could be related to the MIS 5a. No attribution to MIS 5c can be unequivocally supported. Nevertheless, the resulting framework does not allow us to define a satisfactory RSL trend during the MIS 5e highstand and subsequent MIS 5 substages. Overall, the proposed selection of reliable/unreliable data would be useful for detecting areas where MIS 5 substage attributions are not supported by confident U/Th chronological data and thus the related reconstructions need to be revised. In this regard, the resulting framework calls for a reappraisal and re-examination of the Mediterranean records with advanced geochronological methodologies.     

Models of interacting pairs of thin,quasi-geostrophic vortices: steady-state solutions and nonlinear stability

We study pairwise interactions of elliptical quasi-geostrophic (QG) vortices as the limiting case of vanishingly thin uniform potential vorticity ellipsoids. In this limit, the product of the vertical extent of the ellipsoid and the potential vorticity within it is held fixed to a finite non-zero constant. Such elliptical “lenses” inherit the property that, in isolation, they steadily rotate without changing shape. Here, we use this property to extend both standard moment models and Hamiltonian ellipsoidal models to approximate the dynamical interaction of such elliptical lenses. By neglecting non-elliptical deformations, the simplified models reduce the dynamics to just four degrees of freedom per vortex. For simplicity, we focus on pairwise interactions between identical elliptical vortices initially separated in both the horizontal and vertical directions. The dynamics of the simplified models are compared with the full QG dynamics of the system, and show good agreement as expected for sufficiently distant lenses. The results reveal the existence of families of steadily rotating equilibria in the initial horizontal and vertical separation parameter space. For sufficiently large vertical separations, equilibria with varying shape exist for all horizontal separations. Below a critical vertical separation (stretched by the constant ratio of buoyancy to Coriolis frequencies N / f), comparable to the mean radius of either vortex, a gap opens in horizontal separation where no equilibria are possible. Solutions near the edge of this gap are unstable. In the full QG system, equilibria at the edge of the gap exhibit corners (infinite curvature) along their boundaries. Comparisons of the model results with the full nonlinear QG evolution show that the early stages of the instability are captured by the Hamiltonian elliptical model but not by the moment model that inaccurately estimates shorter-range interactions.     

Seismic-induced rockfalls and landslide dam following the October 30, 2016 earthquake in Central Italy

On October 30, 2016, a seismic event and its aftershocks produced diffuse landslides along the SP 209 road in the Nera River Gorge (Central Italy). Due to the steep slopes and the outcropping of highly fractured and bedded limestone, several rockfalls were triggered, of which the main event occurred on the slope of Mount Sasso Pizzuto. The seismic shock acted on a rock wedge that, after an initial slide, developed into a rockfall. The debris accumulation blocked the SP 209 road and dammed the Nera River, forming a small lake. The river discharge was around 3.6 m³/s; the water overtopped the dam and flooded the road. By a preliminary topographic survey, we estimated that the debris accumulation covers an area of about 16,500 m², while the volume is around 70,000 m³. The maximum volume occupied by the pre-existing talus mobilized by the rockfall is about 20% of the total volume. Besides blocking the road, the rockfall damaged a bridge severely, while, downstream of the dam, the water flow caused erosion of a road embankment. A rockfall protection gallery, a few hundred meters downstream of the dam, was damaged during the event. Other elastic nets and rigid barriers were not sufficient to protect the road from single-block rockfalls, with volumes around 1–2 m³. Considering the geological and geomorphological conditions, as well as the high seismicity and the socioeconomic importance of the area, a review of the entire rockfall protection systems is required to ensure protection of critical infrastructure and local communities.     

The influence of slope gradient and gully channel on the run-out behavior of rockslide-debris flow: an analysis on the Verghereto landslide in Italy

Landslides - Rockslide-debris flow may cause catastrophic damages because of its high speed and long run-out distance. The influences of slope gradient and gully channel on the entrainment,...     

Mapping landslide susceptibility with frequency ratio,statistical index,and weights of evidence models: a case study in northern Iran

The main objective of this study is to investigate potential application of frequency ratio (FR), weights of evidence (WoE), and statistical index (SI) models for landslide susceptibility mapping in a part of Mazandaran Province, Iran. First, a landslide inventory map was constructed from various sources. The landslide inventory map was then randomly divided in a ratio of 70/30 for training and validation of the models, respectively. Second, 13 landslide conditioning factors including slope degree, slope aspect, altitude, plan curvature, stream power index, topographic wetness index, sediment transport index, topographic roughness index, lithology, distance from streams, faults, roads, and land use type were prepared, and the relationships between these factors and the landslide inventory map were extracted by using the mentioned models. Subsequently, the multi-class weighted factors were used to generate landslide susceptibility maps. Finally, the susceptibility maps were verified and compared using several methods including receiver operating characteristic curve with the areas under the curve (AUC), landslide density, and spatially agreed area analyses. The success rate curve showed that the AUC for FR, WoE, and SI models was 81.51, 79.43, and 81.27, respectively. The prediction rate curve demonstrated that the AUC achieved by the three models was 80.44, 77.94, and 79.55, respectively. Although the sensitivity analysis using the FR model revealed that the modeling process was sensitive to input factors, the accuracy results suggest that the three models used in this study can be effective approaches for landslide susceptibility mapping in Mazandaran Province, and the resultant susceptibility maps are trustworthy for hazard mitigation strategies.     

Crash-testing the cauldron code for joint lensing and dynamics analysis of early-type galaxies

We apply the joint lensing and dynamics code for the analysis of early-type galaxies, 'Combined Algorithm for Unified Lensing and Dynamics ReconstructiON ( cauldron )', to a rotating N -body stellar system with dark matter halo which significantly violates the two major assumptions of the method, i.e. axial symmetry supported by a two-integral distribution function. The goal is to study how cauldron performs in an extreme case, and to determine which galaxy properties can still be robustly recovered. Three data sets, corresponding to orthogonal lines of sight, are generated from the N -body system and analysed with the identical procedure followed in the study of real lens galaxies, adopting an axisymmetric power-law total density distribution. We find that several global properties of the N -body system are recovered with remarkable accuracy, despite the fact that the adopted power-law model is too simple to account for the lack of symmetry of the true density distribution. In particular, the logarithmic slope of the total density distribution is robustly recovered to within less than 10 per cent (with the exception of the ill-constrained very inner regions), the inferred angle-averaged radial profile of the total mass closely follows the true distribution, and the dark matter fraction of the system (inside the effective radius) is correctly determined within ∼10 per cent of the total mass. Unless the line-of-sight direction is almost parallel to the total angular momentum vector of the system, reliably recovered quantities also include the angular momentum, the   V /σ  ratio and the anisotropy parameter δ. We conclude that the cauldron code can be safely and effectively applied to real early-type lens galaxies, also providing reliable information for the systems that depart significantly from the method's assumptions.