Volcanological evolution of the Rivi–Capo Volcanic Complex at Salina,Aeolian Islands: magma storage processes and ascent dynamics

Critical to understanding explosive eruptions is establishing how accurately representative pyroclasts are of processes during magma vesiculation and fragmentation. Here, we present data on densities, and vesicle size and number characteristics, for representative pyroclasts from six silicic eruptions of contrasting size and style from Raoul volcano (Kermadec arc). We use these data to evaluate histories of bubble nucleation, coalescence, and growth in explosive eruptions and to provide comparisons with pumiceous dome carapace material. Density/vesicularity distributions show a scarcity of pyroclasts with ～65–75 % vesicularity; however, pyroclasts closest to this vesicularity range have the highest bubble number density (BND) values regardless of eruptive intensity or style. Clasts with vesicularities greater than this 65–75 % “pivotal” vesicularity range have decreasing BNDs with increasing vesicularities, interpreted to reflect continuing bubble growth and coalescence. Clasts with vesicularities less than the pivotal range have BNDs that decrease with decreasing vesicularity and preserve textures indicative of processes such as stalling and open system degassing prior to vesiculation in a microlite-rich magma, or vesiculation during slow ascent of degassing magma. Bubble size distributions (BSDs) and BNDs show variations consistent with 65–75 % representing the vesicularity at which vesiculating magma is most likely to undergo fragmentation, consistent with the closest packing of spheres. We consider that the observed vesicularity range may reflect the development of permeability in the magma through shearing as it flows through the conduit. These processes can act in concert with multiple nucleation events, generating a situation of heterogeneous bubble populations that permit some regions of the magma to expand and bubbles to coalesce with other regions in which permeable networks are formed. Fragmentation preserves the range in vesicularity seen as well as any post-fragmentation/pre-quenching expansion which may have occurred. We demonstrate that differing density pyroclasts from a single eruption interval can have widely varying BND values corresponding to the degree of bubble maturation that has occurred. The modal density clasts (the usual targets for vesicularity studies) have likely undergone some degree of bubble maturation and are therefore may not be representative of the magma at the onset of fragmentation.     

Bismuth‐melt trails trapped in cassiterite–quartz veins

Native bismuth in the form of metallic melt has been considered instrumental to the formation of some metallic ore deposits via a mechanism dubbed the “Liquid Bismuth Collector Model.” Here, we provide petrographical documentation of trail‐forming, μm‐sized blebs of native bismuth in cassiterite–quartz veins from the Santa Bárbara greisen Sn deposit in the Rondônia tin province of northern Brazil. These inclusions suggest the trapping of a Bi melt that took place during vein formation, in a mechanism similar to the entrapment of fluid inclusions.     

Kinematics and significance of a poly-deformed crustal-scale shear zone in central to south-eastern Madagascar: the Itremo–Ikalamavony thrust

International Journal of Earth Sciences - Across the crystalline basement of Madagascar, late Archaean rocks of the Antananarivo Block are tectonically overlain by Proterozoic, predominantly...     

Towards a predictive model to assess the natural position of the Posidonia oceanica seagrass meadows upper limit

The upper portion of the meadows of the protected Mediterranean seagrass Posidonia oceanica occurs in the region of the seafloor mostly affected by surf-related effects. Evaluation of its status is part of monitoring programs, but proper conclusions are difficult to draw due to the lack of definite reference conditions. Comparing the position of the meadow upper limit with the beach morphodynamics (i.e. the distinctive type of beach produced by topography and wave climate) provided evidence that the natural landwards extension of meadows can be predicted. An innovative model was therefore developed in order to locate the region of the seafloor where the meadow upper limit should lie in natural conditions (i.e. those governed only by hydrodynamics, in absence of significant anthropogenic impact). This predictive model was validated in additional sites, which showed perfect agreement between predictions and observations. This makes the model a valuable tool for coastal management.     

Roles of methane and carbon dioxide in geological processes on Mars

We discuss in this paper possible roles of methane and carbon dioxide in geological processes on Mars. These volatiles in the martian crust may migrate upward from their sources either directly or via various traps (structural, sedimentary, ground ice, gas hydrates). They are then likely emitted to the atmosphere by seepage or through diverse vent structures. Though gas hydrates have never been directly detected on Mars, theoretical studies favor their presence in the crust and polar caps; they could have played an important role as significant gas reservoirs in the subsurface. The martian gas hydrates would possibly be a binary system of methane and carbon dioxide occupying clathrate cavities. Landforms such as mud volcanoes with well-known linkage to gas venting are extensively distributed on Earth, and methane is the primary gas involved. Thus, identification of these landforms on Mars could suggest that methane and possibly carbon dioxide have contributed to geological processes of the planet. For example, we present a newly identified field in Chryse Planitia where features closely resembling terrestrial mud volcanoes occur widely, though with no observable activity. We also present results of a preliminary search for possible recent or present-day, methane-emission zones in the regions over which enrichments of atmospheric methane have been reported.     

MODFLOW-Style parameters in underdetermined parameter estimation

In this article, we discuss the use of MODFLOW-Style parameters in the numerical codes MODFLOW_2005 and MODFLOW_2005-Adjoint for the definition of variables in the Layer Property Flow package. Parameters are a useful tool to represent aquifer properties in both codes and are the only option available in the adjoint version. Moreover, for overdetermined parameter estimation problems, the parameter approach for model input can make data input easier. We found that if each estimable parameter is defined by one parameter, the codes require a large computational effort and substantial gains in efficiency are achieved by removing logical comparison of character strings that represent the names and types of the parameters. An alternative formulation already available in the current implementation of the code can also alleviate the efficiency degradation due to character comparisons in the special case of distributed parameters defined through multiplication matrices. The authors also hope that lessons learned in analyzing the performance of the MODFLOW family codes will be enlightening to developers of other Fortran implementations of numerical codes.     

An alternative approach to transverse and profile terrain curvature

Terrain curvature is one of the most important parameters of land surface topography. Well-established methods used in its measurement compute an index of plan or profile curvature for every single cell of a digital elevation model (DEM). The interpretation of these outputs may be delicate, especially when selected locations have to be analyzed. Furthermore, they involve a high level of simplification, contrasting with the complex and multiscalar nature of the surface curvature itself. In this paper, we present a new method to assess vertical transverse and profile curvature combining real-scale visualization and the possibility to measure these two terrain derivatives over a large range of scales. To this purpose, we implemented a GIS tool that extracts longitudinal and transverse elevation profiles from a high-resolution DEM. The performance of our approach was compared with some of the most commonly used methods (ArcMap, Redlands, CA, USA; ArcSIE, Landserf) by analyzing the terrain curvature around charcoal production sites in southern Switzerland. The different methods produced comparable results. While conventional methods quickly summarize terrain curvature in the form of a matrix of values, they involve a loss of information. The advantage of the new method lies in the possibility to measure and visualize the shape and size of the curvature, and to obtain a realistic representation of the average curvature for different subsets of spatial points. Moreover, the new method makes it possible to control the conditions in which the index of curvature is calculated.     

Hysteretic cyclic response of concrete columns reinforced with smooth bars

The application of smooth (plain) bars in reinforced concrete (RC) construction has been abandoned since the 1970s; however, there are many old reinforced concrete buildings in the world whose construction is based on this old style that are now in need of structural seismic rehabilitation according to the requirements of present day seismic rehabilitation codes. The focus of this study concerns the investigation of the hysteretic cyclic response of RC columns with smooth bars. The results of six column specimens having a variety of details for overlapping splices of longitudinal bars while experiencing two different levels of axial loads under cyclic loading reversals are presented. Through analysis of test observations and the obtained experimental results, it is attempted to clarify major aspects of hysteretic response for RC columns with smooth bars, from a seismic assessment point of view. The hysteretic force–drift responses of columns are deeply investigated and a new concept explaining the flag shape form of the hysteretic response is presented. Furthermore, the rocking response of columns is predicted with a new formulation that assumes an internal compression strut inside the column body as a consequence of rocking that originated from high base rotations. Finally, a simple hysteresis rule is proposed which is the result of considering the combination of two springs in parallel to provide the total hysteretic response as the summation of rocking hysteretic and bottom anchor (smooth bar) hysteretic responses.     

Motif Discovery on Seismic Amplitude Time Series: The Case Study of Mt Etna 2011 Eruptive Activity

Algorithms searching for similar patterns are widely used in seismology both when the waveforms of the events of interest are known and when there is no a priori-knowledge. Such methods usually make use of the cross-correlation coefficient as a measure of similarity; if there is no a-priori knowledge, they behave as brute-force searching algorithms. The disadvantage of these methods, preventing or limiting their application to very large datasets, is computational complexity. The Mueen–Keogh (MK) algorithm overcomes this limitation by means of two optimization techniques—the early abandoning concept and space indexing. Here, we apply the MK algorithm to amplitude time series retrieved from seismic signals recorded during episodic eruptive activity of Mt Etna in 2011. By adequately tuning the input to the MK algorithm we found eight motif groups characterized by distinct seismic amplitude trends, each related to a different phenomenon. In particular, we observed that earthquakes are accompanied by sharp increases and decreases in seismic amplitude whereas lava fountains are accompanied by slower changes. These results demonstrate that the MK algorithm, because of its particular features, may have wide applicability in seismology.