Calcic amphiboles in calc-alkaline and alkaline magmas: thermobarometric and chemometric empirical equations valid up to 1,130°C and 2.2 GPa

The following article presents constraints of the stability of Mg-rich (Mg/(Mg + Fe²⁺) > 0.5) calcic amphibole in both calc-alkaline and alkaline magmas, testing of previous thermobarometers, and formulation of new empirical equations that take into consideration a large amount of literature data (e.g. more than one thousand amphibole compositions among experimental and natural crystals). Particular care has been taken in choosing a large number of natural amphiboles and selecting quality experimental data from literature. The final database of experimental data, composed of 61 amphiboles synthesized in the ranges of 800–1,130°C and 130–2,200 MPa, indicates that amphibole crystallization occurs in a horn-like P–T stability field limited by two increasing curves (i.e. the thermal stability and an upper limit), which should start to bend back to higher pressures. Among calcic amphiboles, magnesiohornblendes and tschermakitic pargasites are only found in equilibrium with calc-alkaline melts and crystallize at relatively shallow conditions (P up to ~1 GPa). Kaersutite and pargasite are species almost exclusively found in alkaline igneous products, while magnesiohastingsite is equally distributed in calc-alkaline and alkaline rocks. The reliability of previous amphibole applications was checked using the selected experimental database. The results of this testing indicate that none of the previous thermobarometers can be successfully used to estimate the P, T and fO₂ in a wide range of amphibole crystallization conditions. Multivariate least-square analyses of experimental amphibole compositions and physico-chemical parameters allowed us to achieve a new thermobarometric model that gives reasonably low uncertainties (T ± 23.5°C, P ± 11.5%, H₂O_melt ± 0.78wt%) for calc-alkaline and alkaline magmas in a wide range of P–T conditions (up to 1,130°C and 2,200 MPa) and ∆NNO values (±0.37 log units) up to 500 MPa. The ^AK-^[4]Al relation in amphibole can be readily used to distinguish crystals of calc-alkaline liquids from those of alkaline magmas. In addition, several chemometric equations allowing to estimate the anhydrous composition of the melts in equilibrium with amphiboles of calc-alkaline magmas were derived.     

Assessment of RegCM4 simulated inter-annual variability and daily-scale statistics of temperature and precipitation over Mexico

The skill of a regional climate model (RegCM4) in capturing the mean patterns, interannual variability and extreme statistics of daily-scale temperature and precipitation events over Mexico is assessed through a comparison of observations and a 27-year long simulation driven by reanalyses of observations covering the Central America CORDEX domain. The analysis also includes the simulation of tropical cyclones. It is found that RegCM4 reproduces adequately the mean spatial patterns of seasonal precipitation and temperature, along with the associated interannual variability characteristics. The main model bias is an overestimation of precipitation in mountainous regions. The 5 and 95 percentiles of daily temperature, as well as the maximum dry spell length are realistically simulated. The simulated distribution of precipitation events as well as the 95 percentile of precipitation shows a wet bias in topographically complex regions. Based on a simple detection method, the model produces realistic tropical cyclone distributions even at its relatively coarse resolution (dx = 50 km), although the number of cyclone days is underestimated over the Pacific and somewhat overestimated over the Atlantic and Caribbean basins. Overall, it is assessed that the performance of RegCM4 over Mexico is of sufficient quality to study not only mean precipitation and temperature patterns, but also higher order climate statistics.     

A focal plane detector design for a wide-band Laue-lens telescope

The energy range above 60 keV is important for the study of many open problems in high energy astrophysics such as the role of Inverse Compton with respect to synchrotron or thermal processes in GRBs, non thermal mechanisms in SNR, the study of the high energy cut-offs in AGN spectra, and the detection of nuclear and annihilation lines. Recently the development of high energy Laue lenses with broad energy bandpasses from 60 to 600keV have been proposed for a Hard X ray focusing Telescope (HAXTEL) in order to study the X-ray continuum of celestial sources. The required focal plane detector should have high detection efficiency over the entire operative range, a spatial resolution of about 1mm, an energy resolution of a few keV at 500keV and a sensitivity to linear polarization. We describe a possible configuration of the focal plane detector based on several CdTe/CZT pixelated layers stacked together to achieve the required detection efficiency at high energy. Each layer can operate both as a separate position sensitive detector and polarimeter or work with other layers to increase the overall photopeak efficiency. Each layer has a hexagonal shape in order to minimize the detector surface required to cover the lens field of view. The pixels would have the same geometry so as to provide the best coupling with the lens point spread function and to increase the symmetry for polarimetric studies.     

Agricultural Land Use Effects on Climate over China as Simulated by a Regional Climate Model

The Abdus Salam International Centre for Theoretical Physics (ICTP) Regional Climate Model version 3 (RegCM3) is used to investigate the climate effects of land use change related to agriculture over China. The model is driven by the European Center for Medium-range Weather Forecast 40-yr Re-Analysis (ERA40)data. Two sets of experiments for 15 yr (1987-2001) are conducted, one with the potential vegetation cover and the other the agricultural land use (AG). The results show that the AG effects on temperature are weak over northern China while in southern China a significant cooling is found in both winter (December-January-February) and summer (June-July-August). The mean cooling in the sub-regions of South China (SC) in winter and the sub-regions of Southeast (SE) China in summer are found to be the greatest,up to 0.5℃ and 0.8℃, respectively. In general, the change of AG leads to a decrease of annual mean temperature by 0.5-1℃ in southern China. Slight change of precipitation in western China and a decrease of precipitation in eastern China are simulated in winter, with the maximum reduction reaching -7.5% over SE. A general decrease of precipitation over northern China and an increase over southern China are simulated in summer,in particular over SE where the increase of precipitation can be up to 7.3%. The AG effects on temperature and precipitation show strong interannual variability. Comparison of the climate effects between AG and the present-day land use (LU) is also performed. In southern China, the ratio of temperature (precipitation)changes caused by AG and LU is greater than (closer to) the ratio of the number of grid cells with changed vegetation cover due to AG and LU variations.     

Development and validation of a regional coupled atmosphere lake model for the Caspian Sea Basin

We present a validation analysis of a regional climate model coupled to a distributed one dimensional (1D) lake model for the Caspian Sea Basin. Two model grid spacings are tested, 50 and 20 km, the simulation period is 1989–2008 and the lateral boundary conditions are from the ERA-Interim reanalysis of observations. The model is validated against atmospheric as well as lake variables. The model performance in reproducing precipitation and temperature mean seasonal climatology, seasonal cycles and interannual variability is generally good, with the model results being mostly within the observational uncertainty range. The model appears to overestimate cloudiness and underestimate surface radiation, although a large observational uncertainty is found in these variables. The 1D distributed lake model (run at each grid point of the lake area) reproduces the observed lake-average sea surface temperature (SST), although differences compared to observations are found in the spatial structure of the SST, most likely as a result of the absence of 3 dimensional lake water circulations. The evolution of lake ice cover and near surface wind over the lake area is also reproduced by the model reasonably well. Improvements resulting from the increase of resolution from 50 to 20 km are most significant in the lake model. Overall the performance of the coupled regional climate—1D lake model system appears to be of sufficient quality for application to climate change scenario simulations over the Caspian Sea Basin.     

Numerical simulation of the tsunami generated by a past catastrophic landslide on the volcanic island of Ischia, Italy

The island of Ischia, Gulf of Naples, Italy, like many other volcanic islands is affected by mass failures, that are mainly related to secondary volcanic processes such as slope steepening and seismic shaking. The block resurgence of its main relief, Mount Epomeo, has been recognised to contribute cyclically to mass instability and cause landslides, that occasionally may reach the sea and start tsunamis. In this work we explore the consequences of the Ischia Debris Avalanche (IDA), a flank collapse that occurred in historical times, and involved the whole Mount Epomeo edifice including its submarine portion, and that may have caused gigantic sea waves affecting all the coasts of Ischia and of the Gulf of Naples. The IDA and the generated tsunami have been taken as the worst-case scenario for the occurrence of a new tsunami in the area. They have been simulated through numerical codes developed and maintained by the University of Bologna. The simulation shows that the IDA-induced tsunami attacks severely all the coasts of the Gulf of Naples with the highest waves obtained for the island of Ischia, the island of Capri and the peninsula of Sorrento. The propagation pattern of the IDA tsunami can be used to get hints on the impact that such an event may have had on early populations habiting Gulf of Naples, but also to get clues on the area that could be most severely hit by a tsunami generated by a smaller-scale landslide that may occur in the same source zone.     

A peak acceleration threshold for soil liquefaction: lessons learned from the 2012 Emilia earthquake (Italy)

National and international seismic codes and recommendations provide criteria for liquefaction exclusion based on a peak ground acceleration (PGA) threshold value. In this paper, after a brief review of the procedures and the values suggested in those documents, a database of liquefaction case histories was created, exploiting the background data used in the most relevant verification charts, currently employed in research and professional practice. This dataset was used to identify, on the basis of simple statistical analyses, a PGA threshold on the free ground surface below which liquefaction is unlikely to occur, regardless of the geological site conditions. The calculated value, which is on the order of 0.07–0.1 g, based on the model employed to fit the data, was analyzed in light of information collected during the 2012 Emilia seismic sequence in Italy, which produced many liquefaction events triggered by low acceleration values. The case history of the Emilia earthquake advises setting a PGA threshold for code and recommendations at the lower probability level of occurrence, in the order of 1 %.     

Stability and chemical equilibrium of amphibole in calc-alkaline magmas: an overview, new thermobarometric formulations and application to subduction-related volcanoes

This work focuses on a rigorous analysis of the physical–chemical, compositional and textural relationships of amphibole stability and the development of new thermobarometric formulations for amphibole-bearing calc-alkaline products of subduction-related systems. Literature experimental results (550–1,120°C, <1,200 MPa, −1 ≤ ΔNNO ≤ +5), H₂O–CO₂ solubility models, a multitude of amphibole-bearing calc-alkaline products (whole-rocks and glasses, representing 38 volcanoes worldwide), crustal and high-P (1–3 GPa) mantle amphibole compositions have been used. Calcic amphiboles of basalt-rhyolite volcanic products display tschermakitic pargasite (37%), magnesiohastingsite (32%) and magnesiohornblende (31%) compositions with aluminium number (i.e. Al# = ^[6]Al/Al_T) ≤ 0.21. A few volcanic amphiboles (~1%) show high Al# (>0.21) and are inferred to represent xenocrysts of crustal or mantle materials. Most experimental results on calc-alkaline suites have been found to be unsuitable for using in thermobarometric calibrations due to the high Al# (>0.21) of amphiboles and high Al₂O₃/SiO₂ ratios of the coexisting melts. The pre-eruptive crystallization of consistent amphiboles is confined to relatively narrow physical–chemical ranges, next to their dehydration curves. The widespread occurrence of amphiboles with dehydration (breakdown) rims made of anhydrous phases and/or glass, related to sub-volcanic processes such as magma mixing and/or slow ascent during extrusion, confirms that crystal destabilization occurs with relatively low T–P shifts. At the stability curves, the variance of the system decreases so that amphibole composition and physical–chemical conditions are strictly linked to each other. This allowed us to retrieve some empirical thermobarometric formulations which work independently with different compositional components (i.e. Si^*, Al_T, Mg^*, ^[6]Al^*) of a single phase (amphibole), and are therefore easily applicable to all types of calc-alkaline volcanic products (including hybrid andesites). The Si^*-sensitive thermometer and the fO₂–Mg^* equation account for accuracies of ±22°C (σ_est) and 0.4 log units (maximum error), respectively. The uncertainties of the Al_T-sensitive barometer increase with pressure and decrease with temperature. Near the P–T stability curve, the error is <11% whereas for crystal-rich (porphyritic index i.e. PI > 35%) and lower-T magmas, the uncertainty increases up to 24%, consistent with depth uncertainties of 0.4 km, at 90 MPa (~3.4 km), and 7.9 km, at 800 MPa (~30 km), respectively. For magnesiohornblendes, the ^[6]Al^*-sensitive hygrometer has an accuracy of 0.4 wt% (σ_est) whereas for magnesiohastingsite and tschermakitic pargasite species, H₂O_melt uncertainties can be as high as 15% relative. The thermobarometric results obtained with the application of these equations to calc-alkaline amphibole-bearing products were finally, and successfully, crosschecked on several subduction-related volcanoes, through complementary methodologies such as pre-eruptive seismicity (volcano-tectonic earthquake locations and frequency), seismic tomography, Fe–Ti oxides, amphibole–plagioclase, plagioclase–liquid equilibria thermobarometry and melt inclusion studies. A user-friendly spreadsheet (i.e. AMP-TB.xls) to calculate the physical–chemical conditions of amphibole crystallization is also provided.     

Lessons learned from two case histories of seismic microzonation in Italy

The prediction of the variability of the seismic ground motion in a given built-up area is considered an effective tool to plan appropriate urban development, to undertake actions on seismic risk mitigation and to understand the damage pattern caused by a strong-motion event. The procedures for studying the seismic response and the seismic microzonation of an urban area are well established; nevertheless, some controversial points still exists and are discussed here. In this paper, the selection of a reference input motion, the construction of a subsoil model and the seismic response analysis procedures are discussed in detail, based on the authors’ experience in two Italian case histories: the seismic microzonation of the city of Benevento, which was a predictive study, and the simulation of seismic response and damage distribution in the village of San Giuliano di Puglia, which was a retrospective analysis.