Monitoring Etna volcanic plumes using a scanning LiDAR

In this paper, we use data obtained from LiDAR measurements during an ash emission event on 15 November 2010 at Mt. Etna, in Italy, in order to evaluate the spatial distribution of volcanic ash in the atmosphere. A scanning LiDAR system, located at 7?km distance from the summit craters, was directed toward the volcanic vents and moved in azimuth and elevation to analyse different volcanic plume sections. During the measurements, ash emission from the North East Crater and high degassing from the Bocca Nuova Crater were clearly visible. From our analysis we were able to: (1) evaluate the region affected by the volcanic plume presence; (2) distinguish volcanic plumes containing spherical aerosols from those having non-spherical ones; and (3) estimate the frequency of volcanic ash emissions. Moreover, the spatial distribution of ash mass concentration was evaluated with an uncertainty of about 50?%. We found that, even during ash emission episodes characterised by low intensity like the 15 November 2010 event, the region in proximity of the summit craters should be avoided by air traffic operations, the ash concentration being greater than 4?×?10^?3?g/m³. The use of a scanning permanent LiDAR station may usefully monitor the volcanic activity and help to drastically reduce the risks to aviation operations during the frequent Etna eruptions.     

Electron accepting capacity of dissolved organic matter as determined by reaction with metallic zinc

Information about the chemical electron accepting capacity (EAC) of dissolved organic matter (DOM) is scarce owing to a lack of applicable methods. We quantified the electron transfer from metallic Zn to natural DOM in batch experiments at DOC concentrations of 10–100 mg-C L^− 1 and related it to spectroscopic information obtained from UV-, synchronous fluorescence, and FTIR- spectroscopy. The electron donating capacity of DOM and pre-reduced DOM was investigated using Fe(CN)₆³⁻ as electron acceptor. Presence of DOM resulted in release of dissolved Zn, consumption of protons, and slower release of hydrogen compared to reaction of metallic Zn with water at pH 6.5. Comparison with reaction stoichiometry confirmed that DOM accepted electrons from metallic Zn. The release of dissolved Zn was dependent on pH, DOC concentration, ionic strength, and organic matter properties. The reaction appeared to be completed within about 24 h and was characterized by pseudo first order kinetics with rate constants of 0.5 to 0.8 h^− 1. EAC per mass unit of carbon ranged from 0.22 mmol g^− 1 C to 12.6 mmol g^− 1 C. Depending on the DOM, a calculated 28–127% of the electrons transferred from metallic Zn to DOM could be subsequently donated to Fe(CN)₆³⁻. EAC decreased with DOC concentration, and increased with aromaticity, carboxyl, and phenolic content of the DOM. The results indicate that an operationally defined EAC of natural DOM can be quantified by reaction with metallic Zn and that DOM properties control the electron transfer. Shortcomings of the method are the coagulation and precipitation of DOM during the experiment and the production of hydrogen and dissolved Zn by reaction of metallic Zn with water, which may influence the determined EAC.     

Conservation of Mediterranean seascapes: analyses of existing protection schemes

Marine protected areas (MPAs) are aimed at managing and protecting marine environments. Their design, however, often disregards both a thorough knowledge of the distribution of habitats and assemblages and the use of proper experimental evaluations of the efficacy of MPAs by comparing protected vs. unprotected zones. About 200 MPAs have been recently instituted in the Mediterranean area, but the evidence of their efficacy is scant. The MPA of Torre Guaceto (Southern Adriatic Sea, Italy) is one of the rare cases of effective protection enforcement. The reserve was instituted more than 10 years ago, a period currently considered as sufficient to show responses by organisms to protection. The MPA is divided into a C zone, the general reserve, where many activities are permitted, a B zone, the partial reserve where restrictions increase and two A zones, the integral reserve where access is prohibited. The goals of the paper were to map the distribution of benthic assemblages to assess if they were properly represented in the differently protected zones, and to test the efficacy of protection by quantifying possible differences between the assemblages in two control areas and in the two A zones, where human impact is completely excluded. The analysis of habitat and assemblage distribution within the MPA showed that the zones with total protection do not include most valuable environmental types. Most of the considered variables (i.e. cover of substratum, number of taxa, and average abundance of the most common taxa) were not significantly different in and out of the A zones, at each time of sampling. Results, however, suggested a possible effect of protection in modifying patterns of abundance of sponges under Cystoseira canopy (more abundant in the fully protected zone). In the subtidal habitat, differences were found in the structure of the whole assemblage and in the abundance of encrusting coralline red algae (more abundant outside the fully protected area). Notwithstanding the correct general methodology employed in the study, a lack of statistical power could have a role in preventing the detection of ecologically relevant effects of protection. In some instances, data pooling allowed a discrimination between cases where there was clearly no effect of protection and cases where there might be. On this basis, the optimization of this experimental design should be considered in further studies. In any case, if the goals of MPAs have not been clearly stated, efficacy of protection might prove very difficult to test even with the use of sound experimental designs.     

Records of environmental and climatic changes during the late Holocene from Svalbard: palaeolimnology of Kongressvatnet

A multi-core, multidisciplinary palaeolimnological study of the partially varved sediment of a deep, meromictic, arctic lake, Kongressvatnet (Svalbard, Western Spitsbergen), provides a record of environmental and climatic changes during last ca. 1800 years. The chronology of sedimentation was established using several dating techniques (¹³⁷Cs, ²¹⁰Pb, varve counts, palaeomagnetic correlation). A multiproxy record of palaeolimnological variability was compiled based on sedimentation rates, magnetic properties, varve thickness, organic matter, geochemistry, pigments from algal and photosynthetic bacteria, mineralogy and biological assemblages (diatoms, Cladocera). The major features recognised in our master core K99-3 include a shift in sediment source and supply (magnetic measurements, geochemistry) probably caused by glaciological changes in the catchment around 38–32 cm core depth (AD 700–820). Additional environmental changes are inferred at 20–18, 8–4.5 and 3–2 cm (AD ca. 1160–1255; 1715–1880; 1940–1963, respectively). During the past ca. 120 years a prominent sedimentological change from brownish-grey, partly laminated silt-clay (varves) to black organic-rich deposits was observed. From AD 1350 to AD1880 the sediment is comprised of a continuous sequence of varves, whereas the earlier sediments are mostly homogeneous with only a few short intercalated laminated sections between AD 860 and 1350. Sedimentation and accumulation rates increased during the last 30 years (modern warming). Pigment concentrations are very low in the lower ca. 32 cm of the core (AD 820) probably because of the high turbidity high energy environment. The high sulphur content in the uppermost 32 cm of sediment has given rise to two horizontally stratified populations of sulphur anaerobic photosynthetic bacteria, as inferred from their specific carotenoids. These bacteria populations are much more abundant during the Little Ice Age (LIA) than during warmer periods (e.g., during the Medieval Warm Period and 20th century). Diatoms are lacking from the core base up to 18 cm (ca. AD 1255); at this level, species indicative of mesotrophic water are present, whereas from 17 cm to the top of the core, oligotrophic taxa such as Staurosira construens/S. pinnata complex dominate, indicating extended ice coverage and more oligotrophic waters during the LIA. The concentration of Cladocera subfossil remains (dominated by Chydorus) are relatively high in the deepest sections (54–32 cm), whereas the upper 32 cm are characterized by a very low concentration of remains, possibly because of the strongly anoxic conditions, and in this upper sediment section rotifer resting eggs become prevalent. We interpret these changes as responses to climate forcing through its impact on glacial melt water, lake ice cover duration and mainly redox conditions in deep water. The observed changes suggest that at least some of our recorded changes may parallel the Greenland Ice core, although our study added more details about the inferred climatic changes. Further aspects are discussed, such as catchment processes, glacial activity, duration of the Medieval Warm Period, the Little Ice Age, local human activity, and limnology.     

Mechanisms of bubble coalescence in silicic magmas

Bubble coalescence is an important process that strongly affects magmatic degassing. Without coalescence, bubbles remain isolated from one another in the melt, severely limiting gas release. Despite this fact, very little has been done to identify coalescence mechanisms from textures of magmatic rocks or to quantify the dynamics of bubble coalescence in melts. In this paper, we present a systematic study of bubble-coalescence mechanisms and dynamics in natural and experimentally produced bubbly rhyolite magma. We have used a combination of natural observations aided by high-resolution X-ray computed tomography, petrological experiments, and physical models to identify different types of bubble?Cbubble interaction that lead to coalescence on the timescales of magma ascent and eruption. Our observations and calculations suggest that bubbles most efficiently coalesce when inter-bubble melt walls thin by stretching rather than by melt drainage from between converging bubble walls. Orders of magnitude are more rapid than melt drainage, bubble wall stretching produces walls thin enough that inter-bubble pressure gradients may cause the melt wall to dimple, further enhancing coalescence. To put these results into volcanogical context, we have identified magma ascent conditions where each coalescence mechanism should act, and discuss the physical conditions for preserving coalescence structures in natural pumice. The timescales we propose could improve volcanic eruption models, which currently do not account for bubble coalescence. Although we do not address the effect of shear strain on bubble coalescence, the processes discussed here may operate in several different eruption regimes, including vesiculation of lava domes, post-fragmentation frothing of vulcanian bombs, and bubbling of pyroclasts in conduits.     

Shallow velocity model of the Northern Flank of Stromboli Volcano, deduced by high frequency surface wave dispersion

The wavefield produced by the Stromboli volcano explosion quakes shows a significant amount of surface waves. Rayleigh waves recorded by a linear array have been investigated to infer the shear-wave velocity model of the Stromboli northern flank. The group velocity dispersion curve was obtained using the multiple filter technique, while the phase velocity dispersion curve was calculated both by phase-matched filtering and performing a p– stack on the observed waveforms. Through the inversion of these curves we were able to recover the shear-wave structure to a depth of about 190 m.     

Naturally occurring radionuclides and their geochemical interactions at a geothermal site in the North German Basin

The activities of the most common, naturally occurring radionuclides ²³⁸U, ²²⁶Ra, ²¹⁰Pb, ²²⁸Ra, ²²⁸Th, and ⁴⁰K were measured by gamma-ray spectrometry in samples from reservoir rocks, geothermal fluids, and mineral precipitates at the geothermal research site Groß Schönebeck (North German Basin). Results demonstrated that the specific activity of the reservoir rock is within the range of the mean concentration in the upper earth crust of <800 Bq/kg for ⁴⁰K and <60 Bq/kg for radionuclides of the ²³⁸U and ²³²Th series, respectively. The geothermal fluid showed elevated activity concentrations (up to 100 Bq/l) for ²²⁶Ra, ²¹⁰Pb, and ²²⁸Ra, as compared to concentrations found in natural groundwater. Their concentration in filter residues even increased up to 100 Bq/g. These residues contain predominantly two different mineral phases: a Sr-rich barite (Sr, BaSO₄) and laurionite (PbOHCl), which both precipitate upon cooling from the geothermal fluid. Thereby they presumably enrich the radionuclides of Ra (by substitution of Ba) and Pb. Analysis of these precipitates further showed an increased ²²⁶Ra/²²⁸Ra ratio from around 1–1.7 during the initial months of fluid production indicating a change in fluid composition over time which can be explained by different contributions of stimulated reservoir rock areas to the overall produced fluid.     

Adaptive POD model reduction for solute transport in heterogeneous porous media

We study the applicability of a model order reduction technique to the solution of transport of passive scalars in homogeneous and heterogeneous porous media. Transport dynamics are modeled through the advection-dispersion equation (ADE) and we employ Proper Orthogonal Decomposition (POD) as a strategy to reduce the computational burden associated with the numerical solution of the ADE. Our application of POD relies on solving the governing ADE for selected times, termed snapshots. The latter are then employed to achieve the desired model order reduction. We introduce a new technique, termed Snapshot Splitting Technique (SST), which allows enriching the dimension of the POD subspace and damping the temporal increase of the modeling error. Coupling SST with a modeling strategy based on alternating over diverse time scales the solution of the full numerical transport model to its reduced counterpart allows extending the benefit of POD over a prolonged temporal window so that the salient features of the process can be captured at a reduced computational cost. The selection of the time scales across which the solution of the full and reduced model are alternated is linked to the Péclet number (P e), representing the interplay between advective and dispersive processes taking place in the system. Thus, the method is adaptive in space and time across the heterogenous structure of the domain through the combined use of POD and SST and by way of alternating the solution of the full and reduced models. We find that the width of the time scale within which the POD-based reduced model solution provides accurate results tends to increase with decreasing P e. This suggests that the effects of local-scale dispersive processes facilitate the POD method to capture the salient features of the system dynamics embedded in the selected snapshots. Since the dimension of the reduced model is much lower than that of the full numerical model, the methodology we propose enables one to accurately simulate transport at a markedly reduced computational cost.