Controls on low-pressure anatexis

Low-pressure anatexis, whereby rocks melt in place after passing through the andalusite stability field, develops under more restricted conditions than does low-pressure metamorphism. Our thermal modelling and review of published work indicate that the following mechanisms, operating alone, may induce anatexis in typical pelitic rocks without inducing wholesale melting in the lower crust: (i) magmatic advection by pervasive flow; (ii) crustal-scale detachment faulting; and (iii) the presence of a high heat-producing layer. Of these, only magmatic advection by pervasive flow and crustal-scale detachment faulting have been shown quantitatively to provide sufficient heat to cause widespread melting. Combinations of the above mechanisms with pluton-scale magmatic advection, shear heating, removal of the lithospheric mantle, or with each other provide additional means of developing suitable high temperatures at shallow crustal levels to generate low-pressure anatexis.     

Testing earthquake prediction methods: «The West Pacific short-term forecast of earthquakes with magnitude MwHRV ≥ 5.8»

Analyzing the tables and probability maps posted by Yan Y. Kagan and David D. Jackson in April 2002–September 2004 at http://scec.ess.ucla.edu/~ykagan/predictions_index.html and the catalog of earthquakes for the same period, the conclusion is drawn that the underlying method could be used for prediction of aftershocks, while it does not outscore random guessing when main shocks are considered.     

3-D seismic structure of the source area of the 1993 Latur, India, earthquake and its implications for rupture nucleations

The Latur earthquake (M_w 6.1) of 29 September 1993 is a rare stable continental region (SCR) earthquake that occurred on a previously unknown blind fault. In this study, we determined detailed three-dimensional (3-D) P- and S-wave velocity (V_p, V_s) and Poisson's ratio (σ) structures by inverting the first P- and S-wave high-quality arrival time data from 142 aftershocks that were recorded by a network of temporary seismic stations. The source zone of the Latur earthquake shows strong lateral heterogeneities in V_p, V_s and σ structures, extending in a volume of about 90 × 90 × 15 km³. The mainshock occurred within, but near the boundary, of a low-V_p, high-V_s and low-σ zone. This suggests that the structural asperities at the mainshock hypocenter are associated with a partially fluid-saturated fractured rock in a previously unknown source zone with intersecting fault surfaces. This might have triggered the 1993 Latur mainshock and its aftershock sequence. Our results are in good agreement with other geophysical studies that suggest high conductivity and high concentration of radiogenic helium gas beneath the source zone of the Latur earthquake. Our study provides an additional evidence for the presence of fluid related anomaly at the hidden source zone of the Latur earthquake in the SCR and helps us understand the genesis of damaging earthquakes in the SCR of the world.     

Aftershocks series monitoring of the September 18, 2004 M = 4.6 earthquake at the western Pyrenees: A case of reservoir-triggered seismicity?

On September 18, 2004, a 4.6 mbLg earthquake was widely felt in the region around Pamplona, at the western Pyrenees. Preliminary locations reported an epicenter less than 20 km ESE of Pamplona and close to the Itoiz reservoir, which started impounding in January 2004. The area apparently lacks of significant seismic activity in recent times. After the main shock, which was preceded by series of foreshocks reaching magnitudes of 3.3 mbLg, a dense temporal network of 13 seismic stations was deployed there to monitor the aftershocks series and to constrain the hypocentral pattern. Aftershock determinations obtained with a double-difference algorithm define a narrow epicentral zone of less than 10 km², ESE–WNW oriented. The events are mainly concentrated between 3 and 9 km depth. Focal solutions were computed for the main event and 12 aftershocks including the highest secondary one of 3.8 mbLg. They show mainly normal faulting with some strike-slip component and one of the nodal planes oriented NW–SE and dipping to the NE. Cross-correlation techniques applied to detect and associate events with similar waveforms, provided up to 33 families relating the 67% of the 326 relocated aftershocks. Families show event clusters grouped by periods and migrating from NW to SE. Interestingly, the narrow epicentral zone inferred here is located less than 4 km away from the 111-m high Itoiz dam. These hypocentral results, and the correlation observed between fluctuations of the reservoir water level and the seismic activity, favour the explanation of this foreshock–aftershock series as a rapid response case of reservoir-triggered seismicity, burst by the first impoundment of the Itoiz reservoir. The region is folded and affected by shallow dipping thrusts, and the Itoiz reservoir is located on the hangingwall of a low angle southward verging thrust, which might be a case sensible to water level fluctuations. However, continued seismic monitoring in the coming years is mandatory in this area to infer more reliable seismotectonic and hazard assessments.     

Two-stage metamorphic evolution of the Bemarivo Belt of northern Madagascar: constraints from reaction textures and in situ monazite dating

New results on the pressure–temperature–time evolution, deduced from conventional geothermobarometry and in situ U‐Th‐total Pb dating of monazite, are presented for the Bemarivo Belt in northern Madagascar. The belt is subdivided into a northern part consisting of low‐grade metamorphic epicontinental series and a southern part made up of granulite facies metapelites. The prograde metamorphic stage of the latter unit is preserved by kyanite inclusions in garnet, which is in agreement with results of the garnet (core)‐alumosilicate‐quartz‐plagioclase (inclusions in garnet; GASP) equilibrium. The peak metamorphic stage is characterized by ultrahigh temperatures of ～900–950 °C and pressures of ～9 kbar, deduced from GASP equilibria and feldspar thermometry. In proximity to charnockite bodies, garnet‐sillimanite‐bearing metapelites contain aluminous orthopyroxene (max. 8.0 wt% Al₂O₃) pointing to even higher temperatures of ～970 °C. Peak metamorphism is followed by near‐isothermal decompression to pressures of 5–7 kbar and subsequent near‐isobaric cooling, which is demonstrated by the extensive late‐stage formation of cordierite around garnet. Internal textures and differences in chemistry of metapelitic monazite point to a polyphasic growth history. Monazite with magmatically zoned cores is rarely preserved, and gives an age of c. 737 ± 19 Ma, interpreted as the maximum age of sedimentation. Two metamorphic stages are dated: M1 monazite cores range from 563 ± 28 Ma to 532 ± 23 Ma, representing the collisional event, and M2 monazite rims (521 ± 25 Ma to 513 ± 14 Ma), interpreted as grown during peak metamorphic temperatures. These are among the youngest ages reported for high‐grade metamorphism in Madagascar, and are supposed to reflect the Pan‐African attachment of the Bemarivo Belt to the Gondwana supercontinent during its final amalgamation stage. In the course of this, the southern Bemarivo Belt was buried to a depth of >25 km. Approximately 25–30 Myr later, the rocks underwent heating, interpreted to be due to magmatic underplating, and uplift. Presumably, the northern part of the belt was also affected by this tectonism, but buried to a lower depth, and therefore metamorphosed to lower grades.     

The structure and dynamics of 2-dimensional fluids in swelling clays

The interlayer pores of swelling 2:1 clays provide an ideal 2-dimensional environment in which to study confined fluids. In this paper we discuss our understanding of the structure and dynamics of interlayer fluid species in expanded clays, based primarily on the outcome of recent molecular modelling and neutron scattering studies. Counterion solvation is compared with that measured in bulk solutions, and at a local level the cation-oxygen coordination is found to be remarkably similar in these two environments. However, for the monovalent ions the contribution to the first coordination shell from the clay surfaces increases with counterion radius. This gives rise to inner-sphere (surface) complexes in the case of potassium and caesium. In this context, the location of the negative clay surface charge (i.e. arising from octahedral or tetrahedral substitution) is also found to be of major importance. Divalent cations, such as calcium, eagerly solvate to form outer-sphere complexes. These complexes are able to pin adjacent clay layers together, and thereby prevent colloidal swelling. Confined water molecules form hydrogen bonds to each other and to the clays' surfaces. In this way their local environment relaxes to close to the bulk water structure within two molecular layers of the clay surface. Finally, we discuss the way in which the simple organic molecules methane, methanol and ethylene glycol behave in the interlayer region of hydrated clays. Quasi-elastic neutron scattering of isotopically labelled interlayer CH₃OD and (CH₂OD)₂ in deuterated clay allows us to measure the diffusion of the CH₃- and CH₂-groups in both clay and liquid environments. We find that in both the one-layer methanol solvates and the two-layer glycol solvates the diffusion of the most mobile organic molecules is close to that in the bulk solution.     

2D Self-Potential tomographies for studying groundwater flows in the Varco d''Izzo landslide (Basilicata, southern Italy)

A geoelectrical monitoring activity has been carried out to improve the geological and hydrogeological knowledge about the Varco d'Izzo landslide (Potenza, Basilicata, Southern Apennine, Italy), an active rotational–translational slide evolved in earth-flow. In this work we have focused on the Self-Potential (SP) method by applying three different SP measuring techniques and combining modern technologies for data acquisition and new methods for tomographic inversion. A SP map and three static SP tomographies have been carried out to better analyse the groundwater circulation system and to better reconstruct the geometry of the landslide body. In the accumulation zone, which is the area more exposed to the geomorphological activity, a new SP measuring strategy has been applied. This strategy, based on time-continuous 2D SP tomographies, helps identify water flow changes in subsurface by studying the time series of SP tomographic images. The analysis of time-dependent changes of water infiltration in near surface is the key to better understand the hydrogeological processes underlying the ground instability phenomena. The time-lapse analysis of tomographic images has allowed us not only to investigate the correlation between the temporal changes of SP signals and rainfall events, but also to quantify the range of these changes. The modification of the distribution of the SP source accumulation zones is associated with the dynamics of the groundwater flows. These preliminary results allow us to consider the SP tomographic method as a tool for geophysical monitoring of landslide areas and encourage to develop new measuring systems for near-real time applications.