Magma generation and crustal accretion as evidenced by supra-subduction ophiolites of the Albanide–Hellenide Subpelagonian zone

Abstract Ophiolites of the Mirdita–Subpelagonian zone form a nearly continuous belt in the Albanide–Hellenide orogen, including mid‐ocean ridge basalt (MORB) associations in the western Mirdita sector and supra‐subduction zone (SSZ) complexes, with prevalent island arc tholeiitic (IAT) and minor boninitic affinities in the eastern part of the belt (i.e. eastern Mirdita, Pindos, Vourinos). In addition, basalts with geochemical features intermediate between MORB and IAT (MORB/IAT) are found in the central Mirdita and in the Aspropotamos sequence (Pindos). These basalts alternate with pure MORB and are cut by boninitic dykes. The distinctive compositional characteristics of the mafic magmas parental to the different ophiolitic suites can be accounted for by partial melting of mantle sources progressively depleted by melt extractions. Partial melting processes (10–20%) of lherzolitic sources generated pure MORB, leaving clinopyroxene‐poor lherzolite as a residuum. Approximately 10% water‐assisted partial melting of this latter source, in an SSZ setting, may in turn generate basalts with MORB/IAT intermediate characteristics, whereas IAT basalts and boninites may have been derived from 10–20% and 30% partial melting, respectively, of the same source variably enriched by subduction‐derived fluids. In addition, boninites may also have been derived by comparatively lower degrees of hydrated partial melting of more refractory harzburgitic sources. A generalized petrologic model based on mass balance calculations between bulk rock and mineral compositions, indicate that most of the intrusives (from ultramafic cumulates to gabbronorites and plagiogranites), as well as sheeted dykes and volcanics (from basalts to rhyodacites) forming the bulk crustal section of the SSZ ophiolites, may be accounted for by shallow fractional crystallization from low‐Ti picritic parental magmas very similar in composition to IAT picrites from Pacific intraoceanic arcs. The most appropriate tectono‐magmatic model for the generation of the SSZ Tethyan ophiolites implies low velocity plate‐convergence of the intraoceanic subduction and generation of a nascent arc with IAT affinity and progressive slab roll‐back, mantle diapirism and extension from the arc axis to the forearc region, with generation of MORB/IAT intermediate basalts and boninitic magmas.     

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.     

Estimation and spatial interpolation of rainfall intensity distribution from the effective rate of precipitation

Great emphasis is being placed on the use of rainfall intensity data at short time intervals to accurately model the dynamics of modern cropping systems, runoff, erosion and pollutant transport. However, rainfall data are often readily available at more aggregated level of time scale and measurements of rainfall intensity at higher resolution are available only at limited stations. A distribution approach is a good compromise between fine-scale (e.g. sub-daily) models and coarse-scale (e.g. daily) rainfall data, because the use of rainfall intensity distribution could substantially improve hydrological models. In the distribution approach, the cumulative distribution function of rainfall intensity is employed to represent the effect of the within-day temporal variability of rainfall and a disaggregation model (i.e. a model disaggregates time series into sets of higher solution) is used to estimate distribution parameters from the daily average effective precipitation. Scaling problems in hydrologic applications often occur at both space and time dimensions and temporal scaling effects on hydrologic responses may exhibit great spatial variability. Transferring disaggregation model parameter values from one station to an arbitrary position is prone to error, thus a satisfactory alternative is to employ spatial interpolation between stations. This study investigates the spatial interpolation of the probability-based disaggregation model. Rainfall intensity observations are represented as a two-parameter lognormal distribution and methods are developed to estimate distribution parameters from either high-resolution rainfall data or coarse-scale precipitation information such as effective intensity rates. Model parameters are spatially interpolated by kriging to obtain the rainfall intensity distribution when only daily totals are available. The method was applied to 56 pluviometer stations in Western Australia. Two goodness-of-fit statistics were used to evaluate the skill—daily and quantile coefficient of efficiency between simulations and observations. Simulations based on cross-validation show that kriging performed better than other two spatial interpolation approaches (B-splines and thin-plate splines).     

AVES: A high performance computer cluster array for the INTEGRAL satellite scientific data analysis

In this paper we describe a new computing system array, designed, built and now used at the Space Astrophysics and Planetary Institute (IAPS) in Rome, Italy, for the INTEGRAL Space Observatory scientific data analysis. This new system has become necessary in order to reduce the processing time of the INTEGRAL data accumulated during the more than 9 years of in-orbit operation. In order to fulfill the scientific data analysis requirements with a moderately limited investment the starting approach has been to use a ??cluster?? array of commercial quad-CPU computers, featuring the extremely large scientific and calibration data archive on line.     

Mechanisms linking active rock glaciers and impounded surface water formation in high‐mountain areas

Rock glaciers are slowly flowing mixtures of debris and ice occurring in mountains. They can represent a reservoir of water, and melting ice inside them can affect surface water hydrochemistry. Investigating the interactions between rock glaciers and water bodies is therefore necessary to better understand these mechanisms. With this goal, we elucidate the hydrology and structural setting of a rock glacier–marginal pond system, providing new insights into the mechanisms linking active rock glaciers and impounded surface waters. This was achieved through the integration of waterborne geophysical techniques (ground penetrating radar, electrical resistivity tomography and self‐potentials) and heat tracing. Results of these surveys showed that rock glacier advance has progressively filled the valley depression where the pond is located, creating a dam that could have modified the level of impounded water. A sub‐surface hydrological window connecting the rock glacier to the pond was also detected, where an inflow of cold and mineralised underground waters from the rock glacier was observed. Here, greater water contribution from the rock glacier occurred following intense precipitation events during the ice‐free season, with concomitant increasing electrical conductivity values. The outflowing dynamic of the pond is dominated by a sub‐surface seepage where a minor fault zone in bedrock was found, characterised by altered and highly‐fractured rocks. The applied approach is evaluated here as a suitable technique for investigating logistically‐complex hydrological settings which could be possibly transferred to wider scales of investigation. Copyright © 2017 John Wiley & Sons, Ltd.     

Precursory slope distress prior to the 2010 Mount Meager landslide,British Columbia

In 2010, the south flank of Mount Meager failed catastrophically, generating the largest (53 ± 3.8 × 10⁶ m³) landslide in Canadian history. We document the slow deformation of the edifice prior to failure using archival historic aerial photographs spanning the period 1948–2006. All photos were processed using Structure from Motion (SfM) photogrammetry. We used the SfM products to produce pre-and post-failure geomorphic maps that document changes in the volcanic edifice and Capricorn Glacier at its base. The photographic dataset shows that the Capricorn Glacier re-advanced from a retracted position in the 1980s then rapidly retreated in the lead-up to the 2010 failure. The dataset also documents 60 years of progressive development of faults, toe bulging, and precursory failures in 1998 and 2009. The 2010 collapse was conditioned by glacial retreat and triggered by hot summer weather that caused ice and snow to melt. Meltwater increased pore water pressures in colluvium and fractured rocks at the base of the slope, causing those materials to mobilize, which in turn triggered several secondary failures structurally controlled by lithology and faults. The landslide retrogressed from the base of the slope to near the peak of Mount Meager involving basement rock and the overlying volcanic sequence. Elsewhere on the flanks of Mount Meager, large fractures have developed in recently deglaciated areas, conditioning these slopes for future collapse. Potential failures in these areas have larger volumes than the 2010 landslide. Anticipated atmospheric warming over the next several decades will cause further loss of snow and glacier ice, likely producing additional slope instability. Satellite- and ground-based monitoring of these slopes can provide advanced warning of future landslides to help reduce risk in populated regions downstream.     

Vent temperature trends at the Vulcano Fossa fumarole field: the role of permeability

Between 1994 and 2010, we completed 16 thermal surveys of Vulcano’s Fossa fumarole field (Aeolian Islands, Italy). In each survey, between 400 and 1,200 vent temperatures were collected using a thermal infrared thermometer from distances of ～1?m. The results show a general decrease in average vent temperature during 1994–2003, with the average for the entire field falling from ～220°C in 1994 to ～150°C by 2003. However, between 2004 and 2010, we witnessed heating, with the average increasing to ～190°C by 2010. Alongside these annual-scale field-wide trends, we record a spatial re-organisation of the fumarole field, characterised by shut down of vent zones towards the crater floor, matched by rejuvenation of zones located towards the crater rim. Heating may be expected to be associated with deflation because increased amounts of vaporisation will remove volume from the hydrothermal system Gambino and Guglielmino (J Geophys Res 113:B07402, 2008). However, over the 2004–2010 heating period, no ground deformation was observed. Instead, the number of seismic events increased from a typical rate of 37 events per month during 1994–2000 to 195 events per month during 2004–2010. As part of this increase, we noticed a much greater number of high-frequency events associated with rock fracturing. We thus suggest that the heating event of 2004–2010 was the result of changed permeability conditions, rather than change in the heat supply from the deeper magmatic source. Within this scenario, cooling causes shut down of lower sectors and re-establishment of pathways located towards the crater rim, causing fracturing, increased seismicity and heat flow in these regions. This is consistent with the zone of rejuvenation (which lies towards and at the rim) being the most favourable location for fracturing given the stress field of the Fossa cone Sch?pa et al. (J Volcanol Geotherm Res 203:133–145, 2011); it is also the most established zone, having been active at least since the early twentieth century. Our data show the value of deploying multi-disciplinary geophysical campaigns at degassing (fumarolic) hydrothermal systems. This allows more complete and constrained understanding of the true heat loss dynamics of the system. In the case study presented here, it allows us to distinguish true heating from apparent heating phases. While the former are triggered from the bottom-up, i.e. they are driven by increases in heat supply from the magmatic source, the latter are triggered from the top-down, i.e. by changing permeability conditions in the uppermost portion of the system to allow more efficient heat flow over zones predisposed to fracturing.     

Latest Miocene to Quaternary horizontal and vertical displacement rates during simultaneous contraction and extension in the Southern Apennines orogen, Italy

Foreland contraction and hinterland extension in the Southern Apennines orogen of Italy produced a complex spatial and temporal pattern of vertical and horizontal displacement. Remarkably, Late Miocene to mid-Pleistocene foreland migration of the contractional front at ∼16 mm yr⁻¹ was not accompanied by uplift and the frontal thrust belt remained at or below sea level. Only following a mid-Pleistocene reduction in horizontal displacement did the frontal thrust belt and foreland begin uplift at ∼0.5 mm yr⁻¹, a rate that increased to ∼1 mm yr⁻¹ after 125 ka. Although the extensional hinterland experienced net subsidence during formation of the Tyrrhenian basin, an extensional transition zone adjacent to the frontal thrust belt records sustained uplift at ∼0.3 mm yr⁻¹. The interaction of preexisting crustal structure and deep tectonic processes resulted in time-integrated displacement rates suggesting steady-state deformation for periods of 10⁶ years. Displacement rate changes were abrupt and occurred over intervals of 10⁵ years or less.     

Active thrust tectonics in western Sicily (southern Italy): the 1968 Belice earthquake sequence

The 1968 Belice earthquake sequence, characterized by six main shocks with 5 < M < 5.4, represents the strongest seismic event recorded in western Sicily in historical times. The epicentral area is located in the Belice Valley, a region lacking any topographic lineament likely to result from a fault with significant offsets of any kind. Instrumental data show that hypocentres of the major shocks are distributed along a roughly N-dipping plane extending from about 36 km to 1 km depth. Fault plane solutions show pure thrusting mechanisms on N-dipping, ENE-trending planes, or oblique slip with a right-lateral component of motion along steep WSW-dipping planes, both as a result of approximate N-S shortening. The observed destruction indicates that isoseismal areas are elongated in an ENE direction. Similarly, the epiccntral distribution of events with M ≥ 4 outlines a roughly elliptical ENE-elongated area located 20 km NW of the Sciacca-Rocca Ficuzza thrust front. This ENE-striking structure, representing the regional morphotectonic feature closest to the epicentral area, consists of two main imbricate fan systems. In the southernmost system, Quaternary deposits (tentatively dated as 1.0–0.7 Myr old) are involved in a large ramp anticline uplifting them to a maximum altitude of 346 m. The occurrence of Holocene lacustrine piggy-back basins on the rear of this structure also indicates late Quaternary activity of the underlying thrust. Seismological, structural and morphotectonic observations suggest that multiple ruptures might have occurred during the 1968 sequence on a blind crustal thrust ramp located beneath the epicentral area. Slip propagated southwards along the shallow ramp-flat system characterizing the thin-skinned foreland fold and thrust belt of southwestern Sicily, being dispersed in flexural folding processes and diffuse strain along this path.