Origins and energetics of maar volcanoes: examples from the ultrapotassic Sabatini Volcanic District (Roman Province,Central Italy)

Maar volcanoes represent a common volcano type which is produced by the explosive interaction of magma with external water. Here, we provide information on a number of maars in the ultrapotassic Sabatini Volcanic District (SVD, Roman Province) as young as ∼90 ka. The SVD maars are characterised in terms of crater and ejecta ring morphologies, eruptive successions and magma compositions, in light of the local substrate settings, with the aim of assessing magma–water interaction conditions, eruption energetics and genetic mechanisms. Feeder magmas spanned the whole SVD differentiation trend from trachybasalts–shoshonites to phonolites. From the ejected lithic fragments from aquifer rocks, the range of depth of magma–water explosive interaction is estimated to have been mostly at ∼400–600 m below ground level, with a single occurrence of surficial interaction in palustrine–lacustrine environment. In particular, the interaction with external water may have triggered the explosive behaviour of poorly differentiated magmas, whereas it may have acted only as a late controlling factor of the degree of fragmentation and eruption style for the most differentiated magma batches during low-flux ascent in an incipiently fragmented state. Crater sizes, ejecta volumes and ballistic data allow a reconstruction of the energy budget of SVD maar-forming eruptions. Erupted tephra volumes from either monogenetic or polygenetic maars ranged 0.004–0.07 km³ during individual maar-forming eruptions, with corresponding total magma thermal energies of 8 × 10¹⁵–4 × 10¹⁷ J. Based on energy partitioning and volume balance of erupted magmas and lithic fractions vs. crater holes, we consider the different contributions of explosive excavation of the substrate vs. subsidence in forming the SVD maar craters. Following available models based on crater sizes, highly variable fractions (5–50%) of the magma thermal energies would have been required for crater excavation. It appears that subsidence may have played a major role in some SVD maars characterised by low lithic contents, whilst substrate excavation became increasingly significant with increasing degrees of aquifer fragmentation.     

Non‐impact origin of the crater‐like structures in the Gilf Kebir area (Egypt): Implications for the geology of eastern Sahara

Abstract— Several small crater‐like structures occur in Gilf Kebir region (SW Egypt). It has been previously suggested that they could be the result of meteoritic impacts. Here we outline the results of our geological and geophysical survey in the area. The proposed impact origin for these structures is not supported by our observations and analyses, and we suggest an alternative interpretation. The crater‐like structures in Gilf Kebir area are likely related to endogenic processes typical of hydrothermal vent complexes in volcanic areas which may reflect the emplacement of subvolcanic intrusives.     

Offshore Occurrences of Heavy-Mineral Placers,Northwest Galicia,Spain

ABSTRACT

Results of offshore water-jet sediment sampling demonstrate the potential for interstadial placer deposits of industrial minerals in the shore-face to inner-shelf of Cape Ortegal, NW Spain. Placer deposition is interpreted to have occurred during glacial interstadials since the Last Glacial Maximum, and deposits sampled at 24 to 42 m. BSL were formed during the Younger Dryas, and subsequently preserved as a result of very rapid rising sea-levels. The sediments are almost entirely terrestrially derived with a particularly high heavy mineral content (Sp. Gr. >2.9), of more than 75% in all samples. Reworked garnet from earlier-sedimentary deposits constitutes an important fraction of the placer deposits identified but is generally absent in other sediments. Indications of the presence of +5 m-thick paleo-placer deposits with sections include:

• Placer garnet + spinel grades increase with depth in those sediments at 24 to 42 m BSL to the ocean-floor.

• The d50 of sediments with background grades decrease with depth while those of placer deposits do not display any significant trend.

• Garnets in non-mineralized sediments are almost entirely angular, whereas in placer enriched sediments they vary in form from well-rounded to angular. Some grains in the latter display surface morphologies typical of reworked sediment while none of the former do.

Sediments deposited during glacial interstadials at depths of up to 130 m BSL may represent important, undeveloped exploration targets for offshore placer deposits of industrial minerals. Industrial mineral exploration should extend the depths of investigations to those corresponding to major interstadial periods.     

Seismic response of masonry-infilled steel frames via multi-scale finite-element analyses

Effects of masonry infills on the seismic vulnerability of steel frames is studied through multi-scale numerical modelling. First, a micro-modelling approach is utilized to define a homogenized masonry material, calibrated on experimental tests, which is used for modelling the nonlinear response of a one-story, single span, masonry-infilled portal under horizontal loads. Based on results of the micro-model, the constitutive behavior of a diagonal strut macro-element equivalent to the infill panel is calibrated. Then, the diagonal strut is used to model infill panels in the macro-scale analysis of a multi-span multi-story infilled moment-resisting (MR) steel frame. The seismic vulnerability of the MR frame is evaluated through a nonlinear static procedure. Numerical analyses highlight that infills may radically modify the seismic response and the failure mechanism of the frame, hence the importance of the infill correct modelling.     

Displacement capacity of masonry piers: parametric numerical analyses versus international building codes

The nonlinear behaviour of masonry piers loaded in their plane is investigated by parametric numerical simulations. Each pier has a cantilever scheme, is loaded by a constant axial load and is subjected to an increasing horizontal displacement at the top. The macro-modelling approach is used to perform numerical analyses, adopting two different constitutive laws: a total strain crack model and a plastic model. The numerical model is calibrated on a block-masonry type for which experimental tests are available in literature. Parametric numerical simulations are performed by varying the aspect-ratio and the compression level, in order to assess the influence of such parameters on both shear strength and displacement capacity. By comparing numerical results with formulas of international codes, a good agreement for the shear strength is obtained, while significant differences are observed for the displacement capacity, which is influenced by both parameters. The authors propose a simple empirical formula for the displacement capacity, obtained by fitting the numerical results. The expression can be useful in the practical design for considering the influence of aspect-ratio and compression level, currently neglected by building codes.     

Ground motion observations and simulation for local earthquakes in the Campi Flegrei volcanic area

Ground motion produced by low magnitude earthquakes can be used to predict peak values in high seismic risk areas where large earthquakes data are not available. In the present work 20 local earthquakes (M_D∈[?0.3, 2.2]) occurred in the Campi Flegrei caldera during the last decade were analyzed. We followed this strategy: empirical relations were used to calibrate synthetic modeling, accounting for the source features and wave propagation effects. Once the source and path parameters of ground motion simulation were obtained from the reference data set, we extrapolated scenarios for stronger earthquakes for which real data are not available. The procedure is structured in two steps: (1) evaluation of ground motion prediction equation for Campi Flegrei area and assessment of input parameters for the source, path and site effects in order to use the finite fault stochastic approach (EXSIM code); (2) simulation of two moderate-to-large earthquake scenarios for which only historical data or partial information are available (M_w4.2 and M_w5.4). The results show that the investigated area is characterized by high attenuation of peak amplitude and not negligible site effects. The stochastic approach has revealed a good tool to calibrate source, path and site parameters on small earthquakes and to generate large earthquake scenario. The investigated magnitude range represents a lower limit to apply the stochastic method as a calibration tool, due to the small size of involved faults (fault length around 200/300 m).     

Strong motion recorded during the Emilia 2012 thrust earthquakes (Northern Italy): a comprehensive analysis

A complex seismic sequence characterised by two thrust earthquakes of magnitudes M \(_\mathrm{L}\) 5.9 and M \(_\mathrm{L}\) 5.8 occurred on May 20 and 29, 2012, respectively, and activated the central portion of the Ferrara Arc structure beneath the Po Plain in northern Italy. The sequence, referred to as Emilia 2012, was recorded by the Italian Strong Motion Network, the Irpinia Network, the Friuli Venezia Giulia Network and 15 temporary stations installed by the Civil Protection Department. In this study, we compile and analyse a large dataset that contains 3,273 waveforms from 37 \(M_\mathrm{L} \ge 4.0\) seismic events. The main aim of this paper is to characterise the ground motion induced by the Emilia 2012 seismic sequence and compare it with other data in the Italian strong motion database and to the recent Ground Motion Prediction Equations (GMPEs) developed for northern Italy, all of Italy and Europe. This is achieved by (1) the computation and analysis of the strong motion parameters of the entire Emilia Strong Motion Dataset (ESMD) and (2) a comprehensive investigation of the May 29 event recordings in terms of time–frequency analysis, the ground motion parameters and the response spectra. This detailed analysis was made possible by the temporary Civil Protection Department stations that were installed soon after the May 20 event at several municipalities in the epicentral area. Most of the recordings are characterised by low-frequency content and long durations, which is a result of the thick sedimentary cover that is typical of the Po Plain. The distributions of the observed horizontal peak ground accelerations and velocities (PGAs and PGVs) with distance are generally consistent with the GMPEs. This is particularly true for the data from M \(_\mathrm{L} \ge \) 5.0 (M \(_\mathrm{W}\ge \) 5.0) events, though the data are scattered at distances beyond approximately 60–70 km and show faster attenuation than the European GMPEs. The horizontal components for the May 29 event at two near-fault stations (Mirandola and San Felice sul Panaro) are overestimated by all of the analysed GMPEs. In contrast, the vertical components, which played an important role in the shaking near the source, are underestimated. The May 29 event produced intense velocity pulses on the horizontal components and the highest peak ground acceleration ever recorded in Italy on the vertical component of the Mirandola near-fault station. The ground motion recordings contained in the ESMD significantly enrich the Italian strong motion database. They contribute new information about (1) the possibility of exceeding the largest recorded PGA in Italy, (2) the development of a spectral design that takes into account the role of the vertical component and the extreme variability of the near-fault ground shaking, and (3) the characterisation of the ground motions in deep sedimentary basins.     

A Middle‐Late Triassic 40Ar/39Ar age for the Paasselkä impact structure (SE Finland)

Abstract– ⁴⁰Ar/³⁹Ar dating of recrystallized feldspar glass particles separated from clast‐rich impact melt rocks from the approximately 10 km Paasselkä impact structure (SE Finland) yielded a Middle to Late Triassic (Ladinian‐Karnian) pseudo‐plateau age of 228.7 ± 3.0 (3.4) Ma (2σ). This new age makes Paasselkä the first known Triassic impact structure dated by isotopic methods on the Baltic Shield. The new Paasselkä impact age is, within uncertainty, coeval with isotopic ages recently obtained for the Lake Saint Martin impact structure in Canada, indicating a new Middle to Late Triassic impact crater population on Earth. The comparatively small crater size, however, suggests no relationship between the Paasselkä impact and a postulated extinction event at the Middle/Late Triassic boundary.     

Managing environmental problems in Cuban karstic aquifers

The study area is located along the Dead Sea Rift, the climate is considered arid in its southern margin near the Dead Sea, which is the lowest water reservoir found on the globe (412 m BSL), to semiarid in its northern part. During the last few decades, the water resources became depleted limiting the natural development of the agricultural settlements, which are the most common type of communities in the region. Previous studies suggested that a large amount of freshwater is lost as the result of salinization processes, which occur when fresh groundwater from the mountain aquifers, flow into the saline clastic Neogene aquifer complex. In order to comprehend this complex system, a detailed outlining of the regional hydrogeological system is essential. Since there are no boreholes, which penetrate the aquiferous rock sequences within the Rift, it was necessary to interpolate a large variety of data from several fields of geosciences. The methods applied included geological mapping, geophysical modeling based on interpretation of seismic profiles and geochemical modeling based on chemical and isotopic analysis of runoff, sediments and groundwater. The combined modeling based on results from the different types of analyses implied to several conclusions relevant to the regional water management policy: (1) groundwater becomes saline as it flows from the margins of the Rift to its center. Therefore, it is recommended to exploit it along the foothills of the rift escarpment. (2) Geophysical modeling indicated that the foothills and the Karstic mountain aquifer extend into the subsurface of the valley and can be farther exploited (up to 15 mm³ per annum) by relatively shallow wells. (3) Several mechanisms of groundwater salinization were deciphered: (a) the dense vertical faulting systems act as potential conduits for saline water, which flow-up from deep-seated sources and penetrate into the fresh aquifers. (b) Fresh groundwater in the clastic aquifer complex is rare, furthermore, two evaporates bodies were encountered (Auja and Zaharat el Qurein), also acting as sources for fresh water salinization. (c) Although the quantity of runoff recharge to the Jordan Valley aquifer complex is negligible, the increase in its salt-content (TDS) turns this negligible freshwater recharge to a significant contributor of salts.     

Majoritic garnet grains within shock-induced melt veins in amphibolites from the Ries impact crater suggest ultrahigh crystallization pressures between 18 and 9 GPa

Shock-induced melt veins in amphibolites from the Nördlinger Ries often have chemical compositions that are similar to bulk rock (i.e., basaltic), but there are other veins that are confined to chlorite-rich cracks that formed before the impact and these are poor in Ca and Na. Majoritic garnets within the shock veins show a broad chemical variation between three endmembers: (1) \({}^{\text{VIII}}{{\text{M}^{2+}}_3} {}^{\text{VI}}{\text{Al}}_{2} ({}^{\text{IV}}{\text{SiO}}_{4} )_{3}\) (normal garnet, Grt), (2) \({}^{\text{VIII}}{{\text{M}^{2+}}_3} {}^{\text{VI}}[{\text{M}}^{2 + } ({\text{Si,Ti}})]({}^{\text{IV}}{\text{SiO}}_{4} )_{3}\)  (majorite, Maj), and (3) \({}^{\text{VIII}}({{\text {Na} {\text M}^{2+}}_2}) {}^{\text{VI}}[ ({\text{Si,Ti}}){\text {Al}}]({}^{\text{IV}}{\text{SiO}}_{4} )_{3}\) (Na-majorite₅₀Grt₅₀), whereby M²⁺ = Mg²⁺, Fe²⁺, Mn²⁺, Ca²⁺. In particular, we observed a broad variation in ^VI(Si,Ti) which ranges from 0.12 to 0.58 cations per formula unit (cpfu). All these majoritic garnets crystallized during shock pressure release at different ultrahigh pressures. Those with high ^VI(Si,Ti) (0.36–0.58 cpfu) formed at high pressures and temperatures from amphibole-rich melts, while majoritic garnets with lower ^VI(Si,Ti) of 0.12–0.27 cpfu formed at lower pressures and temperatures from chlorite-rich melts. Furthermore, majoritic garnets with intermediate values of ^VI(Si,Ti) (0.24–0.39) crystallized from melts with intermediate contents of Ca and Na. To the best of our knowledge the ‘MORB-type’ Ca–Na-rich majoritic garnets with maximum contents of 2.99 wt% Na₂O and calculated crystallisation pressures of 16–18 GPa are the most extreme representatives ever found in terrestrial shocked materials. At the Ries, the duration of the initial contact and compression stage at the central location of impact is estimated to only ~ 0.1 s. We used a ~ 200-µm-thick shock-induced vein in a moderately shocked amphibolite to model its pressure–temperature–time (P–T–t) path. The graphic model manifests a peak temperature of ~ 2600 °C for the vein, continuum pressure lasting for ~ 0.02 s, a quench duration of ~ 0.02 s and a shock pulse of ~ 0.038 s. The small difference between the continuum pressure and the pressure of majoritic garnet crystallization underlines the usefulness of applying crystallisation pressures of majoritic garnets from metabasites for calculation of dynamic shock pressures of host rocks. Majoritic garnets of chlorite provenance, however, are not suitable for the determination of continuum pressure since they crystallized relatively late during shock release. An extraordinary glass- and majorite-bearing amphibole fragment in a shock-vein of one amphibolite documents the whole unloading path.