Quantitative measurement of channel-block hydraulic interactions by experimental saturation of a large, natural, fissured rock mass

The hydrodynamic behavior of fissured media relies on the relationships between a few very conductive fractures (channels) and the remaining low-conductivity fractures and matrix (blocks). We made a quantitative measurement of those relationships and their effect on water drainage and storage in a 19,000 m3 natural reservoir consisting of karstified limestones. This reservoir was artificially saturated with water by closing a water gate on the main outlet during a varying time (delta t) fixed by the operator. The water gate was completely or partly closed until the water pressure reached a particular specified value. If the water gate was left completely closed long enough, the water pressure was fixed by the elevation of temporary outlets at the site boundaries. The water elevation within the reservoir was monitored by means of pressure cells located on single fractures representative of the bedding plane and the two families of fractures of the massif network. The comparison of pressure variations with the network geometry allows us to identify a typical double permeability characterized by a few very conductive channels along 10 vertical faults. These channels limit blocks consisting of low-conductivity bedding planes and a rather impervious matrix. Depending on the closure interval, delta t, of the water gate, the total volume of water stored in the reservoir can vary from 0.8 m3 (delta t = 5 minutes) to 18.6 m3 (delta t = 2 days). Such a variance of storage versus closure time is explained by the reservoir's double permeability that is characterized by blocks that saturate much more slowly than channels. If plotted versus time, this injected volume fits a power relationship, according to the saturation state of the blocks. In less than 34 minutes, storage is close to zero in the blocks and to 1.6 to 2 m3 in the channels. For closing times shorter than 1 hour, only 1% of the volume that flows in the channels is stored into the blocks. Depending on the water pressure and for a given delta t = 3000 minutes, the volume of water stored is controlled by the geometry of the joint network and of the aquifer boundaries. Such an experiment shows that the flow is concentrated in about 10% of the fractured network (channels). The water storage that takes place in the 90% remaining fractures (blocks) depends mainly on time during which pressure remains high into the 10% channels. The accurate modeling of such typical double-permeability media needs a careful study of the geometry of the channels whose narrowings modify the flow and induce dam effects that maintain a high pressure gradient with surrounding blocks.     

Morphological analysis of deep-seated gravitational slope deformation (DSGSD) in the western part of the Argentera massif. A morpho-tectonic control?

The western part of the Argentera–Mercantour massif (French Alps) hosts very large currently active landslides responsible of many disorders and risks to the highly touristic valleys of the Mercantour National Park and skiing resorts. A regional scale mapping of gravitational deformations has been compared to the main geo-structures of the massif. A relative chronology of the events has been established and locally compared to absolute ¹⁰Be dating obtained from previous studies. Two types of large slope destabilisations were identified as follows: deep-seated landslides (DSL) that correspond to rock volumes bounded by a failure surface, and deep-seated gravitational slope deformations (DSGSD) defined as large sagging zones including gravitation landforms such as trenches and scarps or counterscarps. Gravitational landforms are mainly collinear to major N140°E and N020°E tectonic faults, and the most developed DSGSD are located in areas where the slope direction is comparable to the orientation of faults. DSL are mostly included within DSGSD zones and located at the slopes foot. Most of DSL followed a similar failure evolution process according to postglacial over steepened topographies and resulting from a progressive failure growing from the foot to the top of the DSGSD that lasts over a 10 ky time period. This massif-scale approach shows that large-scale DSGSD had a peak of activity from the end of the last deglaciation, to approximately 7000 years bp. Both morphologic and tectonic controls can be invoked to explain the gravitational behaviour of the massif slopes.     

On Geoseismic Noise and Helioseismic Oscillations

Izvestiya, Physics of the Solid Earth - Abstract—In their paper published in Izvestiya, Physics of the Solid Earth, G.A. Sobolev et al. (2020) discussed the results of their study. Firstly,...     

The aftershock dynamics of the Sumatra-Andaman earthquake

The aftershocks of the catastrophic Sumatra-Andaman earthquake of December 26, 2004 (M = 9.0) are analyzed in the general context of the theory of critical phenomena. The analysis relies on the idea that, according to this theory, critical transitions have two key properties. The first is that the intensity of the fluctuations in a dynamical system monotonically increases with the approach of the bifurcation point, so that at a certain time instant, a sufficiently strong internal pulse initiates the catastrophe. This transition can be treated as spontaneous. The second property is that the reactance of the dynamical system drastically increases on the approach of the bifurcation. Even a weak external perturbation in the near-threshold interval can result in a catastrophe. In this case, it is reasonable to refer to the critical transition as an induced transition. The aftershocks of the Sumatra-Andaman earthquake are likely to demonstrate the typical features of induced seismicity. First, the strongest aftershock (M = 7.2) occurred 3 h 20 min after the main shock. It could have probably been induced by the round-trip seismic echo. Second, it was found that the spectral density of the aftershock sequence significantly increases at about ～0.3 mHz, which is close to the frequency of the spheroidal mode ₀S₂. This suggests that the spheroidal oscillations of the Earth, which are excited by the main seismic shock, modulate the aftershock activity. Both hypotheses are supported by the analysis of the aftershocks of the Tohoku earthquake of March 11, 2011 (M = 9.0).     

On magnetic precursors of earthquakes

A simple classification is proposed for magnetic precursors of earthquakes. The high-latitude observations of the so-called Big Magnetic Pulses (BMP)—isolated high-amplitude magnetic pulses that sporadically arise against the undisturbed magnetic background—are used for searching for a probable correlation between BMPs and earthquakes. It is shown by the superposed epoch analysis that seismic activity slightly increases after BMPs. A wide maximum in the number of the earthquakes is observed within the first hour after the BMP. Thus, probably, a new type of magnetic precursors has been found. However, the mechanism of the relationship between BMPs and earthquakes remains open to discussion.     

On the Leontovich boundary condition in geoelectromagnetism

The paper is dedicated to the 70th anniversary of the formulation of the Leontovich impedance boundary condition. This boundary condition is applied in geophysics for sounding of the lithosphere and diagnostics of the magnetosphere by electromagnetic waves. The use of the Leontovich condition is shown to open new possibilities for investigating the anharmonicity of MHD oscillations of the magnetosphere. The differential impedance equation, which is derived with the aid of the Leontovich boundary condition, can prove useful in the investigation of the electroconductivity of the lithosphere by the method of inductive sounding. It is shown within the framework of a simple model that the method of the Leontovich parabolic equation makes it possible to find a correction to the Leontovich boundary condition caused by the fact that, in general, the relation between the horizontal components of an electromagnetic field is not local. It is shown that the anharmonicity of MHD oscillations of the magnetosphere coupled with the nonlocality of the boundary condition on the Earth’s surface can lead to an apparent nonlinearity of the surface impedance calculated in accordance with the classical technique of magnetotelluric sounding.     

Correlation between Pc1 electromagnetic activity and earthquakes

This work addresses the experimental study of an important aspect in the general problem of lithosphere-magnetosphere coupling, namely the correlation between the Pc1 electromagnetic waves (0.2–5 Hz) and the earthquakes. Using long time series of measurements from the catalogues of Pc1 and earthquakes, we revealed a new effect: the diurnal Pc1 activity in the middle latitudes is statistically higher the lower the diurnal global seismic activity. We assume that the influence of earthquake-related acoustic waves on the upper atmosphere suppresses the activity of the Pc1 waves.     

Dependence of the aftershock flow on the main shock magnitude

Previously, we predicted and then observed in practice the property of aftershocks which consists in the statistically regular clustering of events in time during the first hours after the main shock. The characteristic quasi-period of clustering is three hours. This property is associated with the cumulative action of the surface waves converging to the epicenter, whereas the quasi-period is mainly determined by the time delay of the round-the-world seismic echo. The quasi-period varies from case to case. In the attempt to find the cause of this variability, we have statistically explored the probable dependence of quasi-period on the magnitude of the main shock. In this paper, we present the corresponding result of analyzing global seismicity from the USGS/NEIC earthquake catalog. We succeeded in finding a significant reduction in the quasiperiod of the strong earthquakes clustering with growth in the magnitude of the main shock. We suggest the interpretation of this regularity from the standpoint of the phenomenological theory of explosive instability. It is noted that the phenomenon of explosive instability is fairly common in the geophysical media. The examples of explosive instability in the radiation belt and magnetospheric tail are presented. The search for the parallels in the evolution of explosive instability in the lithosphere and magnetosphere of the Earth will enrich both the physics of the earthquakes and physics of the magnetospheric pulsations.