As one of the most serious geological disasters, earthquake is of sudden and destructive characteristic. Therefore, it is of great significance to earthquake monitoring and early warning. The phenomenon of surface thermal infrared radiation enhancement is a common precursor of moderate and strong earthquakes and has been used as an important reference information for early warning and short term prediction. A variety of explanations have been given to understand internal mechanism of the above phenomenon, in which the stress-induced heating hypothesis is widely accepted and has been confirmed in the laboratory rock mechanical loading experiments, that is, under ideal conditions in the laboratory, the rock heats up when it is pressed and cools down when it is stretched. Under field conditions in practice, however, weak seismic precursors of thermal anomalies are often interfered by various environmental factors(solar radiation, atmospheric movement and human activities, etc.), and it has not been investigated whether the corresponding relationship between the above crustal compression-extension motions and thermal radiation anomalies can be observed under field conditions. The earth's gravity field, as one of the basic physical fields of the earth, contains the density distribution of crustal structure, which can be served to study the migration of the earth's material, the deformation of the crust and the change of the stress field. In this paper, we use GRACE gravity and MODIS thermal infrared remote sensing data to verify the stress-induced heat hypothesis in the field with Wenchuan earthquake as the time node. Firstly, the crustal mass density obtained by GRACE satellite was compared with thermal infrared radiation. Then, the gravity anomalies extraction method based on maximum shear strain and in-situ temperature method were used to obtain the gravity anomalies and thermal anomalies respectively. Furthermore, the correlation between the two anomalies before the earthquake was detected from the time scale and space scale respectively, and the consistency analysis between the above anomalies and the spatial distribution of the tectonic fault zone was carried out. For this purpose, two important indicators i.e., anomaly intensity and anomaly distribution, were established in time domain and space domain, respectively. The following conclusions could be drawn: 1)The stress-induced heating hypothesis can be verified by remote sensing in field conditions. The warming zone of the crust(positive thermal offset index)corresponds to the compression zone, and the cooling zone(negative thermal offset index)corresponds to the stretching zone. The consistency of positive and negative variation between the crustal mass density and thermal offset index is 88.9%, which provides field observation evidence for the stress-induced heating hypothesis. 2)The spatio-temporal variation of gravity anomalies and thermal anomalies before earthquake has strong correlation. In the time domain, there is a strong correlation between the gravity anomalies and the thermal anomalies, which shows that the intensity of the two anomalies suddenly increases synchronously and reaches the maximum simultaneously three months before the earthquake. In the spatial domain, gravity anomalies mostly occur at the junction of positive and negative values of thermal offset index, which indicates that the spatial distribution of gravity anomalies and thermal anomalies also has a certain correlation. In addition, the two anomalies appear to be distributed along the fault zone for many times, which shows that they are closely related to tectonic activities.
The terrain in southeastern Tibet is steep and the valleys are crisscrossed. Since the Quaternary, glacial ice and debris have blocked the course of the Yarlung Tsangpo River and its tributary river valleys to form giant dammed lakes, and the huge flood deposits formed by the dammed lake outburst floods are often associated with moraines, ice water deposits, lacustrine deposits, aeolian sand or other running water sediments to form complex river valley accumulation landforms. Different types of sediments in alpine and canyon areas are similar in morphology, structure and fabric, and are difficult to distinguish. Grain size and morphological characteristics are the most important structural characteristics of sediment, and the distribution rules are controlled by many factors such as sedimentary environment, physical properties of detrital material, transporting medium and transporting mode, etc., which is an important proxy index for restoring paleoclimate and inverting paleoenvironment. However, the relevant research on identifying sediment types in alpine valley area of southeast Tibet by grain size and morphology index is still in the exploratory stage. In order to understand the particle size characteristics and spatial differentiation laws of outburst flood sediments and the micromorphological characteristics of particle surfaces, we collected 33 samples of Holocene flood retention sediments preserved along the river within about 350km from the outlet of the Jiacha Gorge in the middle reaches of the Yarlung Tsangpo River to Pai Town, and measured them with Malvern 3000 laser diffraction particle size meter and Zeiss Signma scanning electron microscope, combined with digital geomorphology(DEM)data extracted river channel width and steepness coefficient. The features of spatial distribution law of particle size are analyzed, and the following understanding is obtained. The particle size of outburst flood retention deposits is characterized on the whole by fine-silty sand(2.57~5.18Φ)with poor sorting, positive skew and narrow peak state. Two end element models are obtained: The main peak of EM1 terminal element is 3.16Φ, with an average percentage content of 42.7%, which may represent the alluvial characteristics of higher energy of outburst floods in alpine valley areas, and the main peak of EM2 terminal elements is 2.06Φ with an average percentage content of 55.6%, which can be used to indicate the accumulation process of the outburst flood lag deposits. Affected by the width of the river, the EM1 content has a tendency to increase downstream, while EM2 has the opposite trend. The surface microstructure of quartz particles in the outburst flood lag deposits is mainly characterized by mechanical scratches, shell-like fractures, upturn cleavage and cleavage steps, with low structural maturity, mostly angular shape, and rare denudation pores of chemical origin. As a typical representative of climbing sand dunes in the valley area of the semi-humid monsoon area, the genesis of the dunes is of great guiding significance for revealing the source of sand dunes in the valley area of the alpine valley area, identifying paleoflood deposit and aeolian deposit, distinguishing aeolian deposit and paleoflood slackwater deposits on both sides of the riverbank, and windbreak and sand fixation engineering in the Yarlung Tsangpo River. By comparing the particle size and surface micromorphology characteristics of the known outburst flood deposits of the Yarlung Tsangpo River, we believe that the sand source of the Fozhang dunes is mainly from the outburst flood deposits and was transformed later by wind forces.
The occurrence of earthquakes is closely related to the crustal tectonic movement and the migration of earth mass, which consequently cause the changes of the earth‘s gravitational field. Global time-varying gravity field data obtained by GRACE gravity satellite can be used to detect pre-seismic gravity anomalies. For example, gravity signals caused by several large earthquakes, such as the 2005 MW8.6 Indonesia earthquake, the 2010 MW8.8 Chile earthquake and the 2011 MW9.0 Japan earthquake, have been successfully extracted using GRACE data. However, previous studies on GRACE satellite-based seismic gravity changes focused more on the dynamics of the co-seismic gravity field than on the pre-seismic gravity anomalies which are of great significance for the early warning of earthquakes. Moreover, the commonly adopted difference disposal of the gravity field with the gravity field of adjacent months or the average gravity field of many years when obtaining gravity anomalies cannot effectively remove the inherent north-south stripe noise in GRACE data. On the contrary, it is more likely to cause the annihilation of the medium-high order information in GRACE gravity field model, which results in the loss of some gravity information related to tectonic activities. To explore the pre-seismic gravity anomalies in a more refined way, this study proposes a method of characterizing gravity variation based on the maximum shear strain of gravity, inspired by the concept of crustal strain. In other words, the gravity strain tensor is obtained by further calculating the second-order gradient of the increment of disturbance potential after the removal of hydrological disturbance, and then the maximum shear strain of gravity is ultimately generated to characterize the pre-earthquake tectonic activities. Then, to better understand the seismogenic process of the fault zone by further extracting the pre-earthquake anomalous changes, the data of the maximum shear strain time series are analyzed in this study by means of the offset index K to describe the gravity anomaly. Because the maximum shear strain is calculated by the second-order gradient of GRACE gravity field, this method can suppress the stripe noise better than the difference disposal, thus effectively improving the sensitivity of gravity anomaly detection. The exploratory experiments are carried out in the Tibetan plateau and its surrounding area, which locates among the Pacific Ocean, the Indian Ocean and Eurasia, with the highest altitude, most complex topography and frequent strong earthquakes. Ultimately, the Wenchuan earthquake and Nepal earthquake were used as an example to complete the extraction of pre-earthquake gravity anomaly information by the above method, and the pre-earthquake tectonic activity of the fault zones was analyzed. The results show that a large area of gravity anomalies consistent with the spatial distribution of the fault zone appeared on the Longmenshan fault zone during the half a year before the earthquake, and the maximum anomalous value appeared within 50km from the epicenter, while no anomalies appeared during the non-earthquake period. In addition, compared with the traditional methods, the proposed method has a better ability to extract anomaly information of gravity field, which provides a new idea for understanding the dynamic mechanism of large earthquakes using GRACE data.
Two earthquakes with magnitude larger than 7.0 occurred in 2008 and 2014 on the southwestern end of the Altyn Tagh Fault, which is located in the northwestern borderland of Tibetan plateau. Occurrences of these two earthquakes provide important insights into regional geodynamics and potential seismic risk. Layered viscoelastic model is employed in the paper to study the interaction between these two events. We find that most of aftershocks were triggered by coseismic stress produced by the 2008 Yutian earthquake, and the effect of this earthquake is insignificant on the occurrence of the 2014 Yutian earthquake. However, stress transfer by viscoelastic relaxation of postseismic deformation is in favor of occurrence of the 2014 Yutian earthquake. The coseismic and postseismic stress transfer produced by the 2014 Yutian earthquake leads to stress increasing on the western segment of the Altyn Tagh Fault. Since the occurrence time of the last major earthquake on the western segment of the Altyn Tagh Fault is tens of years ago, it should have accumulated large moment deficit on the fault segment. The Altyn Tagh Fault should be considered as a fault with high potential seismic risk.
Using double-difference precise location results of seismic S wave data from 2000 to 2012, distribution and variation of dynamic moving average Qs value are imaged in Hainan Island and adjacent regions through seismic attenuation tomography. We explore regional crustal S-wave attenuation characteristics of temporal and spatial variation combined with seismic activity, volcanoes and heat distribution, and GPS baseline changes. The results show that:(1)Attenuation imaging shows that there is a significant lateral heterogeneity in Qs distribution in the study area and the Qs values are high in the central region and low in the north and south regions of Hainan Island. Low-Qs areas are baseically converged to the north of Wangwu-Wenjiao Fault, to the south of Jianfeng-Diaoluo Fault, and to the east of Puqian-Qinglan Fault. Earthquakes are basically converged to the edge of the transitional regions of high and low Qs values. Heat flow sites and volcanoes zones are located in low-Qs area in the region. (2)There is a strong correlation between dynamic moving average Qs value, seismic activity and Luzhou-Qiongzhong GPS baseline. From 2000 to 2008, the average Qs value of the study area is relatively high, the seismic activity is strong, and Luzhou-Qiongzhong GPS baseline is decreasing. From 2008 to 2012, the average Qs value of the study area shows a downward trend, the seismic activity is weak, and the Luzhou-Qiongzhong GPS baseline displays an elongation trend.
In the paper,we use 296334 Pn rays from 27777 earthquakes recorded by 1354 seismic stations in the eastern Asia region,an area spanning from 15癗-60癗 and 60癊-145癊,to invert the Pn velocity structure of uppermost mantle and its lateral variations. We also discuss the tectonic implication of the lateral variation based on the Pn velocity structure. Inversion results show that the mean velocity of Pn is 8.03km/s in the eastern Asia region. The lateral variation spans from -0.42km/s to 0.41km/s. Two different tectonic regions can be identified according to distributions of Pn velocity and its lateral variation. In the region from the boundary of the Pacific plate to the east of 108癊,i.e.the eastern China,Pn velocity is relatively low. In the region located at west of 108癊,Pn velocity is high. The different distribution of Pn velocity is consistent to the regional geodynamic background. The dominant geodynamic effect is subduction of the Pacific plate in the eastern China and collision of the Indian Ocean plate in the western China. At the same time,Pn velocity is relatively low in active tectonic regions and high in the stable tectonic regions. The Pn velocity also shows negative correlation to the regional seismicity. The higher of the seismicity,the lower of the Pn velocity,and vice versa. The negative correlation indicates the relationship of lateral variation of Pn velocity and tectonic deformation with seismicity level. The region with low Pn velocity usually is a region with high tectonic deformation and intensive seismicity.