Stable isotopes of pedogenic carbonates from the Somma–Vesuvius area,southern Italy,over the past 18 kyr: palaeoclimatic implications

Stable isotopes were measured in the carbonate and organic matter of palaeosols in the Somma–Vesuvius area, southern Italy in order to test whether they are suitable proxy records for climatic and ecological changes in this area during the past 18000 yr. The ages of the soils span from ca. 18 to ca. 3 kyr BP. Surprisingly, the Last Glacial to Holocene climate transition was not accompanied by significant change in δ¹⁸O of pedogenic carbonate. This could be explained by changes in evaporation rate and in isotope fractionation between water and precipitated carbonate with temperature, which counterbalanced the expected change in isotope composition of meteoric water. Because of the rise in temperature and humidity and the progressive increase in tree cover during the Holocene, the Holocene soil carbonates closely reflect the isotopic composition of meteoric water. A cooling of about 2°C after the Avellino eruption (3.8 ka) accounts for a sudden decrease of about 1‰ in δ¹⁸O of pedogenic carbonate recorded after this eruption. The δ¹³C values of organic matter and pedogenic carbonate covary, indicating an effective isotope equilibrium between the organic matter, as the source of CO₂, and the pedogenic carbonate. Carbon isotopes suggest prevailing C₃ vegetation and negligible mixing with volcanogenic or atmospheric CO₂. Copyright © 2000 John Wiley & Sons, Ltd.     

Evaluating overland flow sediment transport capacity

An equation for evaluating the sediment transport capacity of overland flow is a necessary part of a physically based soil erosion model describing sediment detachment and transport as distributed processes. At first, for the hydraulic conditions of small-scale and large-scale roughness, the sediment transport capacity relationship used in the WEPP model is calibrated by Yalin and Govers' equation. The analysis shows that the transport coefficient K_t depends on the Shields parameter, Y, according to a semi-logarithmic (Yalin) or a linear (Govers) equation. The reliability of the semi-logarithmic equation is verified by Smart's, and Aziz and Scott's experimental data. Then the Low's formula, whose applicability is also proved by Smart's, and Aziz and Scott's data, is transformed as a stream power equation in which a stream power coefficient, K_SP, depending on Shields parameter, slope, sediment and water-specific weight, appears. A relationship between transport capacity and effective stream power is also proposed. Finally, the influence of rainfall on sediment transport capacity and the prediction of critical shear stress corresponding to overland flow are examined. © 1998 John Wiley & Sons, Ltd.     

Predicting the equilibrium bed slope in natural streams using a stochastic model for incipient sediment motion

In this paper, a method to predict the equilibrium bed slope in natural streams based on the incipient motion criterion is proposed. The method is based on the criterion suggested by Gessler to simulate the grain size distribution of the armour coat using the concept of critical shear stress of a sediment mixture. In particular, a different expression of the probability for a single particle size to be part of the armour coat is firstly proposed; then, a simple two‐steps criterion is suggested to estimate the safety factor required by the proposed approach. The method is applied in three different Italian regions (Calabria, Basilicata, and Tuscany) and required several field campaigns involving the survey of 251 stream reaches upstream of grade‐control structures. The area including five Calabrian streams is firstly used to test the suitability of the two‐steps approach. In this area, considering the detailed information about the particle size of bed material, the effect of using a simplified grain size distribution for each torrent is also checked. The method is then applied within the two additional areas (Basilicata and Tuscany) in order to check the geographic influence on the prediction of the stable longitudinal profile. A comparison between measured and estimated values of the equilibrium bed slope showed a good agreement for the Calabrian streams, where the number of the investigated stream reaches is greater, and satisfactory results for the torrents located in Basilicata and Tuscany where the field data set is more limited. The overall results encourage the extension of the proposed method to additional field data and suggest the use of this approach by hydraulic designers in order to stabilize the longitudinal profile of natural streams. Copyright © 2011 John Wiley & Sons, Ltd.     

Discriminating the long distance dispersal of fine ash from sustained columns or near ground ash clouds: The example of the Pomici di Avellino eruption (Somma-Vesuvius,Italy)

Ash samples from tephra layers correlated with the Pomici di Avellino (Avellino Pumice) eruption of Somma-Vesuvius were collected in distal archives and their composition and particle morphology investigated in order to infer their behaviour of transportation and deposition. Differences in composition and particle morphologies were recognised for ash particles belonging to the magmatic Plinian and final phreatomagmatic phases of the eruption. The ash particles were dispersed in opposite directions during the two different phases of the eruption, and these directions are also different from that of coarse-grained fallout deposits. In particular, ash generated during magmatic phase and injected in the atmosphere to form a sustained column shows a prevailing SE dispersion, while ash particles generated during the final phreatomagmatic phase and carried by pyroclastic density currents show a general NW dispersion. These opposite dispersions indicate an ash dispersal influenced by both high and low atmosphere dynamics. In particular, the magmatic ash dispersal was first driven by stratospheric wind towards NE and then the falling particles encountered a variable wind field during their settling, which produced the observed preferential SE dispersal. The wind field encountered by the rising ash clouds that accompanied the pyroclastic density currents of the final phreatomagmatic phase was different with respect to that encountered by the magmatic ash, and produced a NW dispersal. These data demonstrate how ash transportation and deposition are greatly influenced by both high and low atmosphere dynamics. In particular, fine-grained particles transported in ash clouds of small-scale pyroclastic density currents may be dispersed over distances and cover areas comparable with those injected into the stratosphere by Plinian, sustained columns. This is a point not completely addressed by present day mitigation plans in case of renewal of activity at Somma-Vesuvius, and can yield important information also for other volcanoes potentially characterised by explosive activity.     

The Afragola settlement near Vesuvius,Italy: The destruction and abandonment of a Bronze Age village revealed by archaeology,volcanology and rock-magnetism

Public works in progress in the Campanian plain north of Somma-Vesuvius recently encountered the remains of a prehistoric settlement close to the town of Afragola. Rescue excavations brought to light a Bronze Age village partially destroyed and buried by pyroclastic density currents (PDCs) of the Vesuvian Pomici di Avellino eruption (3.8 ¹⁴C ka BP) and subsequently sealed by alluvial deposits. Volcanological and rock-magnetic investigations supplemented the excavations.Careful comparison between volcanological and archaeological stratigraphies led to an understanding of the timing of the damage the buildings suffered when they were struck by a series of PDCs. The first engulfed the village, located some 14 km to the north of the inferred vent, and penetrated into the dwellings without causing major damage. The buildings were able to withstand the weak dynamic pressure of the currents and deviate their path, as shown by the magnetic fabric analyses. Some later collapsed under the load of the deposits piled up by successive currents. Stepwise demagnetization of the thermal remanent magnetization (TRM) carried by potsherds embedded in the deposits yields deposition temperatures in the order of 260–320 °C, fully consistent with those derived from pottery and lithic fragments from other distal and proximal sites. The fairly uniform temperature of the deposits is here ascribed to the lack of pervasive air entrainment into the currents. This, in turn, resulted from the lack of major topographical obstacles along the flat plain.The coupling of structural damage and sedimentological analyses indicates that the currents were not destructive in the Afragola area, but TRM data indicate they were still hot enough to cause death or severe injury to humans and animals. The successful escape of the entire population is apparent from the lack of human remains and from thousands of human footprints on the surface of the deposits left by the first PDCs. People were thus able to walk barefoot across the already emplaced deposits and escape the subsequent PDCs. The rapid cooling of the deposits was probably due to both their thinness and heat dissipation due to condensation of water vapour released in the mixture by magma–water interaction.     

A full-scale study of Darcy-Weisbach friction factor for channels vegetated by riparian species

In this article, an open channel flow resistance equation, deduced applying dimensional analysis and incomplete self-similarity condition for the flow velocity distribution, was tested using measurements carried out in a full-scale channel equipped with three types of riparian plants (Salix alba L., Salix caprea L. and Alnus glutinosa L.). In the experimental channel, having banks lined with boulders, the vegetation branches were anchored in a concrete bottom. For each species, the measurements were carried out with plants having different amounts of leaves, different plant density and plant area index. The relationship between the scale factor Γ of the velocity profile and the Froude number was separately calibrated by measurements carried out without and with vegetation. The component of Darcy-Weisbach friction factor corresponding to the riparian vegetation f_v was calculated as the difference between the measured friction factor value (channel grain roughness + vegetation) and that calculated for the channel without vegetation in the same hydraulic conditions. Using these f_v values, the relationship between the scale factor Γ and the Froude number was calibrated. In this last relationship, a scaling coefficient a varying with the investigated vegetation type was introduced. This coefficient, as expected, gives the highest friction factor values for vegetation having branches with leaves. The theoretical flow resistance law, coupled with the relationship for estimating the Γ function having a scaling coefficient different for each investigated vegetation type, allowed an accurate estimate of the Darcy-Weisbach friction factor (errors less than or equal to 20% for 82.6% of the investigated cases). Finally, for the investigated vegetation species that are characterized by a condition with few leaves or leafless, the scaling coefficient a resulted strongly related to the bending stiffness. This analysis demonstrated that the highest Darcy-Weisbach friction factors correspond to vegetation species characterized by the highest values of bending stiffness. The friction factor values calculated for this last condition are characterized by errors that were less than or equal to ±20% for 90.6% of cases.     

Laboratory determination of mechanical and hydraulic properties of chemically grouted sands

The purpose of this paper is to investigate mechanical and hydraulic properties of sands treated with mineral-based grouts through the results of a laboratory test programme consisting of unconfined compression tests (UCS), triaxial bender element tests (BeT) and constant flow permeability tests in triaxial apparatus. An improved apparatus was set up for obtaining high quality, multiple grouted specimens from a single column. Two selected natural sands having different grain sizes were grouted with two mineral-based silica grouts, resulting in different levels of improvement. The behaviour of the sands treated by mineral grouts, in terms of strength, initial stiffness and permeability, was compared with that exhibited by more traditional silicate grouts. The results of this study indicate that sands treated with mineral grouts result in higher strengths, higher initial shear modulus and lower permeability values than the sands treated with the silicate solution. The effect of grout type, effective confining pressure, and sand particle-size on small-strain shear modulus of grouted sand specimens was evaluated. Based on test results, the small strain shear modulus increment from treated to untreated specimens has been correlated with the unconfined compressive strength, obtaining a unique relationship regardless of grout type and grain-size of tested sands.     

Comparison of different methods describing the peak runoff contributing areas during floods

In the last few years, the scientific community has developed several hydrological models aimed at the simulation of hydrological processes acting at the basin scale. In this context, the portion of peak runoff contributing areas represents a critical variable for a correct estimate of surface runoff. Such areas are strongly influenced by the saturated portion of a river basin (influenced by antecedent conditions) but may also evolve during a specific rainfall event. In the recent years, we have developed 2 theoretically derived probability distributions that attempt to interpret these 2 processes adopting daily runoff and flood‐peak time series. The probability density functions (PDFs) obtained by these 2 schematisations were compared for humid river basins in southern Italy. Results highlighted that the PDFs of the peak runoff contributing areas can be interpreted by a gamma distribution and that the PDF of the relative saturated area provides a good interpretation of such process that can be used for flood prediction.     

Volcanic hazard assessment at the restless Campi Flegrei caldera

Eruption forecasting and hazard assessments at the restless Campi Flegrei caldera, within the Neapolitan volcanic area, have been performed using stratigraphical, volcanological, structural and petrological data.On the basis of the reconstructed variation of eruption magnitude through time, we hypothesize that the most probable maximum expected event is a medium-magnitude explosive eruption, fed by trachytic magma. Such an eruption could likely occur in the north-eastern sector of the caldera floor that is under a tensile stress regime, when the ongoing deformation will generate mechanical failure of the rocks. A vent could open also in the western sector, at the intersection of two fault systems contemporaneously activated, as happened in the last eruption at Monte Nuovo. The eruption could likely be preceded by precursors apparent to the population, such as ground deformation, seismicity and increase in gas emissions. It will probably alternate between magmatic and phreatomagmatic phases with the generation of tephra fallout, and dilute and turbulent pyroclastic currents. During and/or after the eruption, the re-mobilization of ash by likely heavy rains, could probably generate mud flows.In order to perform a zoning of the territory in relation to the expected volcanic hazards, we have constructed a comprehensive hazard map. On this map are delimited (I) areas of variable probability of opening of a new vent, (II) areas which could be affected by variable load of fallout deposits, and (III) areas over which pyroclastic currents could flow. The areas in which a vent could likely open have been defined on the basis of the dynamics of the ongoing deformation of the caldera floor. To construct the fallout hazard map we have used the frequency of deposition of fallout beds thicker than 10 cm, the frequency of load on the ground by tephra fallout and the direction of dispersal axes of the deposits of the last 5 ka, and the limit load of collapse for the variable types of roof construction. The pyroclastic-current hazard map is based on the areal distribution and frequency of pyroclastic-current deposits of the last 5 ka.Editorial Responsibility: T. Druitt     

Timing of magma extraction during the Campanian Ignimbrite eruption (Campi Flegrei Caldera)

A core drilled within the northern part of the city of Napoli has offered the unique opportunity to observe in one single sequence the superposition of the four pyroclastic flow units emplaced during the Campanian Ignimbrite (CI) eruption. Such a stratigraphic succession has never been encountered before in natural or in man made exposures. Therefore the CI sequence was reconstructed only on the basis of stratigraphic correlations and compositional data (in literature). The occurrence of four superposed CI flows, together with all the data available (in literature) allowed us to better constrain the chemical stratigraphy of the deposit and the compositional structure of the CI magma chamber. The CI magma chamber includes two cogenetic magma layers, separated by a compositional gap. The upper magma layer was contaminated by interaction with radiogenic fluids. The two magma layers were extruded either individually or simultaneously during the course of the eruption. In the latter case they produced a hybrid magma. But no evidence of input of new geochemically and isotopically distinct magma batches just prior or during the eruption has been found. Comparison with the exposed CI deposits has permitted reconstruction of variable eruption phases and related magma withdrawal and caldera collapse episodes. The eruption was likely to have began with phreatomagmatic explosions followed by the formation of a sustained plinian eruption column fed by the simultaneous extraction from both magma layers. Towards the end of this phase the upward migration of the fragmentation surface and the decrease in magma eruption rate and/or activation of fractures formed an unstable pulsating column that was fed only by the most-evolved magma layer. This plinian phase was followed by the collapse of the eruption column and the beginning of caldera formation. At this stage expanded pyroclastic flows fed by the upper magma layer in the chamber generated. During the following major caldera collapse episode, the maximum mass discharge rate was reached and both magma layers were tapped, generating expanded pyroclastic flows. Towards the end of the eruption, only the deeper and less differentiated magma layer was tapped producing more concentrated pyroclastic flows that traveled short distances.