Fault plane orientations of microearthquakes at Mt. Etna from the inversion of P-wave rise times

A crucial point in the analysis of tectonic earthquakes occurring in a volcanic area is the inference of the orientation of the structures along which the ruptures occur. These structures represent zones of weakness which could favor the migration of melt toward the surface and the assessment of their geometry is a fundamental step toward efficient evaluation of volcanic risk. We analyzed a high-quality dataset of 171 low-magnitude, tectonic earthquakes that occurred at Mt. Etna during the 2002–2003 eruption. We applied a recently developed technique aimed at inferring the source parameters (source size, dip and strike fault) and the intrinsic quality factor Q_p of P waves from the inversion of rise times. The technique is based on numerically calibrated relationships among the rise time of first P waves and the source parameters for a circular crack rupturing at a constant velocity. For the most of the events the directivity source effect did not allow us to constrain the fault plane orientation. For a subset of 45 events with well constrained focal mechanisms we were able to constrain the “true” fault plane orientation. The level of resolution of the fault planes was assessed through a non linear analysis based on the random deviates technique. The significance of the retrieved fault plane solutions and the fit of the assumed source model to data were assessed through a χ-square test. Most of the retrieved fault plane solutions agree with the geometrical trend of known surface faults. The inferred source parameters and Q_p are in agreement with the results of previous studies.     

Entropy production and coarse graining of the climate fields in a general circulation model

We extend the analysis of the thermodynamics of the climate system by investigating the role played by processes taking place at various spatial and temporal scales through a procedure of coarse graining. We show that the coarser is the graining of the climatic fields, the lower is the resulting estimate of the material entropy production. In other terms, all the spatial and temporal scales of variability of the thermodynamic fields provide a positive contribution to the material entropy production. This may be interpreted also as that, at all scales, the temperature fields and the heating fields resulting from the convergence of turbulent fluxes have a negative correlation, while the opposite holds between the temperature fields and the radiative heating fields. Moreover, we obtain that the latter correlations are stronger, which confirms that radiation acts as primary driver for the climatic processes, while the material fluxes dampen the resulting fluctuations through dissipative processes. We also show, using specific coarse-graining procedures, how one can separate the various contributions to the material entropy production coming from the dissipation of kinetic energy, the vertical sensible and latent heat fluxes, and the large scale horizontal fluxes, without resorting to the full three-dimensional time dependent fields. We find that most of the entropy production is associated to irreversible exchanges occurring along the vertical direction, and that neglecting the horizontal and time variability of the fields has a relatively small impact on the estimate of the material entropy production. The approach presented here seems promising for testing climate models, for assessing the impact of changing their parametrizations and their resolution, as well as for investigating the atmosphere of exoplanets, because it allows for evaluating the error in the estimate of their thermodynamical properties due to the lack of high-resolution data. The findings on the impact of coarse graining on the thermodynamic fields on the estimate of the material entropy production deserve to be explored in a more general context, because they provide a way for understanding the relationship between forced fluctuations and dissipative processes in continuum systems.     

Climate entropy budget of the HadCM3 atmosphere–ocean general circulation model and of FAMOUS, its low-resolution version

The entropy budget is calculated of the coupled atmosphere–ocean general circulation model HadCM3. Estimates of the different entropy sources and sinks of the climate system are obtained directly from the diabatic heating terms, and an approximate estimate of the planetary entropy production is also provided. The rate of material entropy production of the climate system is found to be ～50 mW m^?2 K^?1, a value intermediate in the range 30–70 mW m^?2 K^?1 previously reported from different models. The largest part of this is due to sensible and latent heat transport (～38 mW m^?2 K^?1). Another 13 mW m^?2 K^?1 is due to dissipation of kinetic energy in the atmosphere by friction and Reynolds stresses. Numerical entropy production in the atmosphere dynamical core is found to be about 0.7 mW m^?2 K^?1. The material entropy production within the ocean due to turbulent mixing is ～1 mW m^?2 K^?1, a very small contribution to the material entropy production of the climate system. The rate of change of entropy of the model climate system is about 1 mW m^?2 K^?1 or less, which is comparable with the typical size of the fluctuations of the entropy sources due to interannual variability, and a more accurate closure of the budget than achieved by previous analyses. Results are similar for FAMOUS, which has a lower spatial resolution but similar formulation to HadCM3, while more substantial differences are found with respect to other models, suggesting that the formulation of the model has an important influence on the climate entropy budget. Since this is the first diagnosis of the entropy budget in a climate model of the type and complexity used for projection of twenty-first century climate change, it would be valuable if similar analyses were carried out for other such models.     

CO<Subscript>2</Subscript> and He degassing at El Chichón volcano,Chiapas, Mexico: gas flux,origin and relationship with local and regional tectonics

During 2007–2008, three CO₂ flux surveys were performed on El Chichón volcanic lake, Chiapas, Mexico, with an additional survey in April 2008 covering the entire crater floor (including the lake). The mean CO₂ flux calculated by sequential Gaussian simulation from the lake was 1,190 (March 2007), 730 (December 2007) and 1,134 g m⁻² day⁻¹ (April 2008) with total emission rates of 164 ± 9.5 (March 2007), 59 ± 2.5 (December 2007) and 109 ± 6.6 t day⁻¹ (April 2008). The mean CO₂ flux estimated from the entire crater floor area was 1,102 g m⁻² day⁻¹ for April 2008 with a total emission rate of 144 ± 5.9 t day⁻¹. Significant change in CO₂ flux was not detected during the period of survey, and the mapping of the CO₂ flux highlighted lineaments reflecting the main local and regional tectonic patterns. The ³He/⁴He ratio (as high as 8.1 R _A) for gases in the El Chichón crater is generally higher than those observed at the neighbouring Transmexican Volcanic Belt and the Central American Volcanic Arc. The CO₂/³He ratios for the high ³He/⁴He gases tend to have the MORB-like values (1.41 × 10⁹), and the CO₂/³He ratios for the lower ³He/⁴He gases fall within the range for the arc-type gases. The high ³He/⁴He ratios, the MORB-like CO₂/³He ratios for the high ³He/⁴He gases and high proportion of MORB-CO₂ (M = 25 ±15%) at El Chichón indicate a greater depth for the generation of magma when compared to typical arc volcanoes.     

Geochemical evidence of the renewal of volcanic activity inferred from CO<Subscript>2</Subscript> soil and SO<Subscript>2</Subscript> plume fluxes: the 2007 Stromboli eruption (Italy)

On 27 February 2007, a new eruption occurred on Stromboli which lasted until 2 April. It was characterized by effusive activity on the Sciara del Fuoco and by a paroxysmal event (15 March). This crisis represented an opportunity for us to refine the model that had been developed previously (2002–2003 eruption) and to improve our understanding of the relationship between the magmatic dynamics of the volcano and the geochemical variations in the fluids. In particular, the evaluation of the dynamic equilibrium between the volatiles (CO₂ and SO₂) released from the magma and the corresponding fluids discharged from the summit area allowed us to evaluate the level of criticality of the volcanic activity. One of the major accomplishments of this study is a 4-year database of summit soil CO₂ flux on the basis of which we define the thresholds (low–medium–high) for this parameter that are empirically based on the natural volcanological evolution of Stromboli. The SO₂ fluxes of the degassing plume and the CO₂ fluxes emitted from the soil at Pizzo Sopra la Fossa are also presented. It is noteworthy that geochemical signals of volcanic unrest have been clearly identified before, during and after the effusive activity. These signals were found almost simultaneously in the degassing plume (SO₂ flux) and in soil degassing (CO₂ flux) at the summit, although the two degassing processes are shown to be clearly different. The interpretation of the results will be useful for future volcanic surveillance at Stromboli.     

Analysis of the liquefaction phenomena in the village of Vittorito (L��Aquila)

This paper describes a case-history of liquefaction occurred near the village of Vittorito after the April 6, 2009 L’Aquila earthquake (moment magnitude M_w = 6.3), approximately 45 km far from the epicentre. In the document, first, an estimation of the seismic motion in the area has been made. Thereafter, the performed geotechnical investigation is described, followed by the application of some fast assessment criteria for the occurrence of liquefaction, recently proposed by the new Italian Building Code. A careful assessment of all the parameters involved in conventional Seed and Idriss (1971) liquefaction analyses is considered. The cyclic resistance ratio CRR is evaluated by cone penetration tests CPT and by in situ seismic dilatometer tests SMDT; in the latter case CRR is evaluated by different empirical correlations with shear wave velocity V_s and horizontal stress index K_D. Analytical data confirmed the observed occurrence of the liquefaction in Vittorito, even if the acceleration field in the area, produced by the L’Aquila earthquake, was very low.     

Yields of suspended sediment and dissolved solids from the Andean basins of Ecuador

Abstract

Water discharge and suspended and dissolved sediment data from three rivers (Napo, Pastaza and Santiago) in the Ecuadorian Amazon basin and a river in the Pacific basin (Esmeraldas) over a 9-year period, are presented. This data set allows us to present: (a) the chemical weathering rates; (b) the erosion rates, calculated from the suspended sediment from the Andean basin; (c) the spatio-temporal variability of the two regions; and (d) the relationship between this variability and the precipitation, topography, lithology and seismic activity of the area. The dissolved solids load from the Esmeraldas basin was 2 × 10⁶ t year^-1, whereas for the Napo, Pastaza and Santiago basins, it was 4, 2 and 3 × 10⁶ t year^-1, respectively. For stations in the Andean piedmont of Ecuador, the relationship between surface sediment and the total sediment concentration was found to be close to one. This is due to minimal stratification of the suspended sediment in the vertical profile, which is attributed to turbulence and high vertical water speeds. However, during the dry season, when the water speed decreases, sediment stratification appears, but this effect can be neglected in the sediment flux calculations due to low concentration rates. The suspended sediment load in the Pacific basin was 6 × 10⁶ t year^-1, and the total for the three Amazon basins was 47 × 10⁶ t year^-1. The difference between these contributions of the suspended sediment load is likely due to the tectonic uplift and the seismic and volcanic dynamics that occur on the Amazon side.

Editor Z.W. Kundzewicz

Citation Armijos, E., Laraque, A., Barba, S., Bourrel, L., Ceron, C., Lagane, C., Magat, P., Moquet, J.-S., Pombosa, R., Sondag, F., Vauchel, P., Vera, A., and Guyot, J.L., 2013. Yields of suspended sediment and dissolved solids from the Andean basins of Ecuador. Hydrological Sciences Journal, 58 (7), 1478–1494.     

The Ricker Hills Tillite provides evidence of Oligocene warm-based glaciation in Victoria Land, Antarctica

The relationship between the Ricker Hills Tillite (RHT), which represents the northernmost outcrop of lithified continental glacial deposits in Victoria Land, is discussed with respect to the glacial landscape assemblage of the Ricker Hills, a nunatak at the internal border of the Transantarctic Mountains. A warm-based ice sheet deposited the tillite and induced syn- to post-depositional glacial deformation under wet conditions both of the tillite and of the bedrock. The thickness of the ice sheet on the nunatak is estimated to have been 600 m, at most. The area had been deeply eroded before deposition of the RHT as documented by the low elevation of tillite outcrops located in overdeepened depressions of the nunatak. Micropaleontological analysis evidences only the presence of Permian to Jurassic palynomorphs. X-ray diffraction and SEM–EDS analyses of clay minerals in the RHT indicate continental chemical weathering under wet conditions after the RHT deposition. As documented by clay mineral assemblage variation in CRP drillholes, the progressive cooling of the Antarctic continent allowed chemical weathering in “warm” conditions until the Late Oligocene period in southern Victoria Land, leading to a chronological constrain for RHT deposition. Conservatively estimating the sea level to have been between the tillite outcrops and the erosional trimline limiting horns in the Ricker Hills, at the time of RHT deposition, we suggest that the maximum uplift of the area would not have exceeded 900–1500 m since at least Late Oligocene.