CrystalMom: a new model for the evolution of crystal size distributions in magmas with the quadrature-based method of moments

Nucleation and growth of crystals, and the resulting crystal size distribution, play a fundamental role in controlling the physical properties of magmas and consequently the dynamics of the eruptions. In the past decades, laboratory experiments demonstrated that size and shape of crystals strongly control the physical properties of magma and lava. Additionally, natural and experimental samples are usually characterized in terms of their crystal size distribution to link it with physical processes that are not directly observable, such as cooling or decompression mechanisms. In this paper, we present CrystalMoM, a new predictive model, based on the quadrature-based method of moments, developed for studying the kinetic of crystallization in volcanic systems. The quadrature-based method of moments, well established in the field of chemical engineering, represents a mesoscale modelling approach that rigorously simulates the space–time evolution of a distribution of particles, by considering its moments. The method is applied here, for the first time, for studying the equilibrium/disequilibrium crystallization in magma, modelling the temporal evolution of the moments of a crystal size distribution. The model, verified against numerical and experimental data, represents a valuable tool to infer the cooling and decompression rates from the crystal size distribution observed in natural samples.     

Paleomagnetic dating of the most recent silicic eruptive activity at Pantelleria (Strait of Sicily)

We report on the paleomagnetism of ten sites in the products of the most recent silicic eruptive cycle of Pantelleria, Strait of Sicily. Previously radiometrically dated at 5–10 ka, our comparison with proxies for geomagnetic field directions allows us to narrow considerably the time window during which these eruptions occurred. The strongly peralkaline composition causes the magmas to have low viscosities, locally resulting in strong agglutination of proximal fall deposits. This allows successful extraction of paleomagnetic directions from the explosive phases of eruptions. One of our sites was located in the Serra della Fastuca fall deposit, produced by the first explosive event of the eruptive cycle. The other nine sites were located in the most recent explosive (pumice fall and agglutinate from Cuddia del Gallo and Cuddia Randazzo) and effusive (Khaggiar lava) products. The (very similar) paleomagnetic directions gathered from eight internally consistent sites were compared to reference geomagnetic field directions of the last 5–10 ka. Directions from Cuddia del Gallo agglutinate and Khaggiar flows translate into 5.9- to 6.2-ka ages, whereas the Fastuca pumices yield a slightly older age of 6.2–6.8 ka. Hence, the most recent silicic eruptive cycle lasted at most a millennium and as little as a few centuries around 6.0 ka. Paleomagnetically inferred ages are in good agreement with published (and calibrated by us) ¹⁴C dates from paleosols/charcoals sampled below the studied volcanic units, whereas K/Ar data are more scattered and yield ∼30% older ages. Our data show that the time elapsed since the most recent silicic eruptions at Pantelleria is comparable to the quiescence period separating the two latest volcanic cycles.     

Assessment of water vapor retrievals from a GPS receiver network

We present an assessment of a GPS receiver operational network to produce accurate integrated precipitable water vapour (IPWV) during a two-week field experiment carried out in Central Italy around the city of Rome, where different instruments were operative. This experimental activity provided an excellent opportunity to compare the GPS products with independent measurements provided by ground-based and space-based sensors and to evaluate their quality in terms of absolute accuracy of IPWV, analyzing also the spatial scale of GPS estimates. For instance, the assimilation into Numerical Weather Prediction models of IPWV provided by a GPS network or its exploitation in space geodesy applications to correct tropospheric effects requires an accuracy in the order of 0.1 cm to be ascribed to IPWV observations. In this work, we assessed that the accuracy for GPS IPWV estimates is 0.07 cm. Moreover, this experiment has pointed out strengths and limitations of an operational network for the water vapor estimation, such as a proper receiver distribution to achieve the desired spatial resolution and a coverage of GPS stations in both flat and mountains regions.     

Macroseismic effects highlight site response in Rome and its geological signature

A detailed analysis of the earthquake effects on the urban area of Rome has been conducted for the L??Aquila sequence, which occurred in April 2009, by using an online macroseismic questionnaire. Intensity residuals calculated using the mainshock and four aftershocks are analyzed in light of a very accurate and original geological reconstruction of the subsoil of Rome based on a large amount of wells. The aim of this work is to highlight ground motion amplification areas and to find a correlation with the geological settings at a subregional scale, putting in evidence the extreme complexity of the phenomenon and the difficulty of making a simplified model. Correlations between amplification areas and both near-surface and deep geology were found. Moreover, the detailed scale of investigation has permitted us to find a correlation between seismic amplification in recent alluvial settings and subsiding zones, and between heard seismic sound and Tiber alluvial sediments.     

Ostreopsis cf. ovata bloom in the northern Adriatic Sea during summer 2009: ecology, molecular characterization and toxin profile

Intense blooms of the benthic dinoflagellate Ostreopsis cf. ovata have occurred in the northern Adriatic Sea since 2006. These blooms are associated with noxious effects on human health and with the mortality of benthic organisms because of the production of palytoxin-like compounds. The O. cf. ovata bloom and its relationships with nutrient concentrations at two stations on the Conero Riviera (northern Adriatic Sea) were investigated in the summer of 2009. O. cf. ovata developed from August to November, with the highest abundances in September (1.3 × 10⁶ cells g⁻¹ fw corresponding to 63.8 × 10³ cells cm⁻²). The presence of the single O. cf. ovata genotype was confirmed by a PCR assay. Bloom developed when the seawater temperature was decreasing. Nutrient concentrations did not seem to affect bloom dynamics. Toxin analysis performed by high resolution liquid chromatography-mass spectrometry revealed a high total toxin content (up to 75 pg cell⁻¹), including putative palytoxin and all the ovatoxins known so far.     

Sub-seismic scale fracture pattern and in situ permeability data in the chalk atop of the Krempe salt ridge at Lägerdorf, NW Germany: Inferences on synfolding stress field evolution and its impact on fracture connectivity

Chalk is exposed in the Heidestrasse quarry at Lägerdorf, at the top of the NE-SW trending Krempe salt ridge. Structural data indicate the presence of two joint sets, striking almost parallel and perpendicular to the salt ridge, respectively, and of a set of conjugate extensional faults and fault zones striking NW-SE, i.e. almost perpendicular to the salt ridge. Within the overall NW-SE trend of joints and faults, strike variations occur from the massive chalk exposed in the lower half of the quarry, to the overlying layered chalk. A large variability characterizes the normalized spacing of both joint sets, which does not show any clear trend with layer dip. In situ measurements indicate that the cross-sectional permeability of tight joints increases 1-2 orders of magnitude with respect to the undeformed chalk. We propose that joint and fault azimuthal variability resulted from changes through time of the stress ellipsoid at the top of the salt ridge, while joint spacing variability is associated with the weak mechanical influence of bedding in chalk. Azimuthal variability improves fracture connectivity and, hence, permeability and fluid flow.     

Microstructural Investigation of Naturally Deformed Leucogneiss from an Alpine Shear Zone (Southern Calabria – Italy)

In order to reveal any correlation between the preferred orientation of quartz and progressive intensity of strain, a suite of deformed leucogneiss collected within a ductile shear zone outcropping in Calabria (southern Italy) was investigated. Based on the microstructural approach and matrix-clast relationships, the method applied here may be useful for appraising the connection between deformation mechanisms of constituent minerals and bulk textural properties of naturally sheared rocks. Accordingly, quartz c-axis orientation patterns were determined by image-assisted analysis. Results revealed a correlation between finite strain and textural features: As strain increases, the matrix develops at the expense of the clast counterpart, which is instead progressively reduced in both size and amount.     

Transient processes in Stromboli’s shallow basaltic system inferred from dolerite and magmatic breccia blocks erupted during the 5 April 2003 paroxysm

We describe the mineralogy, geochemistry, and mesomicrostructure of fresh subvolcanic blocks erupted during the 5 April 2003 paroxysm of Stromboli (Aeolian Islands, Italy). These blocks represent ∼50 vol.% of the total erupted ejecta and consist of fine- to medium-grained basaltic lithotypes ranging from relatively homogeneous dolerites to strongly or poorly welded magmatic breccias. The breccia components are represented by angular fragments of dolerites entrapped in a matrix of vesiculated (lava-like to scoriae) crystal-rich (CR) basalt. All of the studied blocks are cognates with the CR basalt of the normal Strombolian activity or lavas and they are often coated by a few-centimeter thick layer of crystal-poor (CP) basaltic pumice erupted during the paroxysm. We suggest that they result from the rapid increase of pressure and related subvolcanic rock failure that occurred shortly before the 5 April 2003 explosion, when the uppermost portion of the edifice inflated and suffered brecciation as the result of the sudden rise of the gas-rich CP basalt that triggered the eruption. Dolerites and magmatic matrix of the breccias show major and trace element compositions that match those of the CR basalts erupted during normal Strombolian activity and effusive events at Stromboli volcano. Dolerites consist of (a) phenocrysts normally found in the CR basalts and (b) late-stage magmatic minerals such as sanidine, An_60-28 plagioclase, Fe–Mn-rich olivines (Fo_68-48), phlogopite, apatite, and opaque mineral pairs (magnetite and ilmenite), most of which are never found both in lava flows and scoriae erupted during the persistent explosive activity that characterizes typical Strombolian behavior. Subvolcanic crystallization of the Stromboli CR magma, leading to slowly cooled equivalents of basalts, could result from transient drainage of the magma from the summit craters to lower levels. Fingering and engulfing of the material that collapsed from the summit crater floor into the shallow basaltic system during the late evening of 28 December 2002 coupled with the short break in the summit persistent explosions between December 2002 and March 2003 permitted the CR magma pockets to solidify as dolerites, which were confined to the uppermost portion of the system and thus not involved in the ongoing flank effusive activity. Crystal size distribution of the basaltic blocks and crystallization of the finer-grained (<0.1 mm) mafic minerals of the dolerites over a time interval of ∼100 days closely agrees with the above interpretation. Vesicle filling (miarolitic cavities) locally found in some dolerites, with minerals deposited as vapor-phase crystallization is a result of continuous gas percolation through the rocks of the uppermost portion of the volcanic system. Poorly welded magmatic breccias formed during syn-eruptive processes of 5 April 2003, when the paroxysm strongly shattered the shallow subvolcanic system and many dolerite fragments were entrapped in the CR magma. In contrast, the high degree of welding between the dolerite clasts and the CR basaltic matrix in the strongly welded magmatic breccias provides a snapshot of subvolcanic intrusions of the CR basalt into the dolerite when, after a 2-month break in activity, CR magmas started to rise again to the summit craters. Blocks similar to these subvolcanic ejecta of 5 April 2003 were also erupted during previous paroxysms (e.g., 1930) suggesting that changes in the usual Strombolian activity (e.g., short breaks in the persistent mild explosions and/or flank effusive activity) lead to transient crystallization of dolerites in the shallow plumbing system.     

Dissolution/crystallization kinetics recorded in the 2002–2003 lavas of Stromboli (Italy)

On 28 December 2002, new vents opened on the flanks of Stromboli, just below the summit craters, interrupting the persistent activity of the volcano with a 7-month-long effusive eruption. We here report on the plagioclase size distribution (PlgSD) in lava samples collected following the chronology of the 2002–2003 eruption. Data reveal a linear PlgSD similar to that found in samples of normal Stromboli activity, indicating that the switch from Strombolian explosive to effusive activity is not associated with changes in texture. Nevertheless, the crystal size distribution slopes and intercepts exhibit slight sinusoidal temporal variations that are here ascribed to a magma supply mechanism able to induce “resonance” in the crystal size distribution, with an amplitude that depends on the supply rate.