Mercury indicating inflow zones and ruptures along the Wenchuan Ms 8.0 earthquake fault

During the Wenchuan Fault Scientific Drilling Project,we determined the values of total mercury(HgT)and gaseous elemental mercury(GEM) from drilled cores and drilling mud,respectively.Geochemical analysis shows HgT values ranging from 0.24 to 6.45 ng/g for the Penguan complex and from 2.90 to 137.54 ng/g for T3 sediment.The average levels of HgT for the Penguan complex and T3 sediments are 1.81 ± 0.26 ng/g and23.96 ± 4.80 ng/g,respectively.Major anomalous peaks of HgT appear at depth of 614,731,993 and 1,107 m,which correspond to the long-term high seismic activity during crustal deformation in response to tectonic stresses.Gaseous elemental mercury dissolved in drilling mud was also analyzed.We found fluid inflow zones with high GEM at depths of 590-750 m,suggesting that fluid-filled ruptures exist in the LMS fault zone.It indicates that mercury provides geochemical evidence for inflow zones and ruptures/fault zones in the Wenchuan Ms 8.0 earthquake fault.

Surface rupture and hazard characteristics of the Wenchuan Ms 8.0 earthquake, Sichuan, China

Longmen Shan is located the special joint be-tween Tibetan Plateau inland in the west and Yangtze craton in the east. Consisting of a se-ries of parallel imbricated thrust, it develops, from the west to the east, the Maoxian- Wenchuan, Beichuan-Yingxiu and Pengxian- ﹡This research was supported by China National Natural Science Foundation grant 40841010, 40972083 and China National Science and Technology supporting Plan Foundation grant 2006BAC13B02-07, 2006BAC13B01-604. Guanxian faults. The Wenchuan earthquake is a thrust with strike-slip type, and thre surface ruptures are located on the Beichuan-Yingxiu fault zone and Pengxian-Guanxian fault zone. The surface rupture on the Beichuan-Yingxiu fault shows the thrust and dextral slip charac-teristic. The maximum vertical displacement of the surface rupture is about 10.3 m and the maxi-mum right-lateral displacement is about 5.85m. Though the vertical displacements and the hori-zontal displacements in the different segments have certain differences, as a whole, the ratio of the vertical displacement and the horizontal dis-placement is close to 1:1. The surface rupture on the Pengxian-Guanxian fault shows thrust and dextral characteristic. The rates of vertical dis-placements and the horizontal displacements ones on the most other segments are between 1:3 and 1:2. So the Beichuan- Yingxiu fault is a dextral-slip and thrust fault and the average ver-tical displacement is equal to the average hori-zontal displacement, while the Pengxian- Guan xian fault is thrust fault with a little dextral-slip component. The total intensity area above Ⅵ de-gree of the Wenchuan earthquake is about 333000 km2. The high earthquake intensity line stretches to N40-50°E along Longmen Shan tectonic belt. The rate of the long axis and the minor is betw- een 8:1 and 10:1. Three Ⅺ intensity regions are isolated in distribution. It presents a multipoint instantaneous characteristic of the rupture.

Possible precursory anomalies in ground water level associated with the Wenchuan Ms 8.0 earthquake in 2008,Sichuan,China

We have examined the water level data from 16 wells in Sichuan province,China,recorded before the Wenchuan Ms8.0 earthquake occurred in 2008.We found that the data of 5 among these wells exhibit possible precursory anomalies,which are respectively named the Chuan No.08 well in Deyang,Chuan No.ll well in Pujiang,Chuan No.13 well in Luzhou,Chuan No.22 well in Qionglai and Beichuan.In time durations,these anomalies are of long-,intermediate-and short-term signals which are primarily distributed on the northeast(NE)trending Longmenshan fault zone as well as the parallel Huayingshan fault zone.It seems that the variations of the well levels on the Huayingshan fault zone imply compression while those on the Longmenshan fault zone reflect extension.These anomalies occurred first in the areas outside the epicenter region of the Wenchuan Ms8.0 earthquake,then migrated to the seismogenic fault zone.And they migrated from southwest to northeast in agreement with the rupture process of the main shock.

Airblast evolution initiated by Wangjiayan landslides in the M_(s)8.0 Wenchuan earthquake and its destructive capacity analysis

Airblasts,as one common phenomenon accompanied by rapid movements of landslides or rock/snow avalanches,commonly result in catastrophic damages and are attracting more and more scientific attention.To quantitatively analyze the intensity of airblast initiated by landslides,the Wangjiayan landslide,occurred in the Wenchuan earthquake,is selected here with the landslide propagation and airblast evolution being studied using FLUENT by introducing the Voellmy rheological law.The results reveal that:(1)For the Wangjiayan landslide,its whole travelling duration is only 12 s with its maximum velocity reaching 36 m/s at t=10 s;(2)corresponding to the landslide propagation,the maximum velocity,28 m/s,of the airblast initiated by the landslide also appears at t=10 s with its maximum pressure reaching594.8 Pa,which is equivalent to violent storm;(3)under the attack of airblast,the load suffered by buildings in the airblast zone increases to 1300 Pa at t=9.4 s and sharply decreased to-7000 Pa as the rapid decrease of the velocity of the sliding mass at t=10 s,which is seriously unfavorable for buildings and might be the key reason for the destructive collapse of buildings in the airblast zone of the Wangjiayan landslide.

Analysis of debris flow control effect and hazard assessment in Xinqiao Gully,Wenchuan M_(s)8.0 earthquake area based on numerical simulation

Xinqiao Gully is located in the area of the 2008 Wenchuan M_(s)8.0 earthquake in Sichuan province,China.Based on the investigation of the 2023"6-26"Xinqiao Gully debris flow event,this study assessed the effectiveness of the debris flow control project and evaluated the debris flow hazards.Through field investigation and numerical simulation methods,the indicators of flow intensity reduction rate and storage capacity fullness were proposed to quantify the effectiveness of the engineering measures in the debris flow event.The simulation results show that the debris flow control project reduced the flow intensity by41.05%to 64.61%.The storage capacity of the dam decreases gradually from upstream to the mouth of the gully,thus effectively intercepting and controlling the debris flow.By evaluating the debris flow of different recurrence intervals,further measures are recommended for managing debris flow events.

Identifying drivers of urban landuse changes in the Wenchuan earthquake- affected area by using night-time light data

To learn the process of urban land evolution before and after an earthquake is vital to formulate the urban reconstruction control policies and recovery measures in the earthquake-stricken areas.However,spatiotemporal evolution and its driving factors of urban land in earthquake-prone areas remains limited due to the scarcity of ground observation data.This research,leveraging night-time light remote sensing imagery and land cover data,conducted a comprehensive analysis of the long-term evolution characteristics of urban land in earthquake-prone areas.It introduced methodologies for assessing the socio-economic impact and the primary natural environmental factors driving urban land evolution in these regions.To validate the proposed methods,the 2008 Wenchuan earthquake-affected area in China was selected as a representative study area.The results indicated that the average Digital Number(DN)values in socio-economically impacted areas showed a trend of rising,falling,and then rising again after the earthquake.DN values in three types of damaged areas including Type Ⅱ,Type Ⅲ,and Type Ⅳ exceeded pre-earthquake levels.The analysis of determinative factors influencing urban land evolution revealed that slope and elevation were key elements in controlling urban land expansion before the earthquake,whereas factors such as slope,elevation,lithology,and faults had a stronger influence on urban land expansion after the earthquake.It can be seen that,in view of the differences in the natural conditions of regions for post-disaster reconstruction,the local government need to actively adjust and adapt to urban spatial planning,so as to leverage the scale effect of large-scale inputs of funds,facilities,human resources and other factors after the disaster,thus enhancing resilience and recovery efficiency in response to disaster impacts.

Surface Rupture and Co-seismic Displacement Produced by the Ms 8.0 Wenchuan Earthquake of May ^(12)th,2008,Sichuan,China:Eastwards Growth of the Qinghai-Tibet Plateau

An earthquake of Ms 8 struck Wenchuan County, western Sichuan, China, on May 12^th, 2008 and resulted in long surface ruptures （〉300 km）. The first-hand observations about the surface ruptures produced by the earthquake in the worst-hit areas of Yingxiu, Beichuan and Qingchuan, ascertained that the causative structure of the earthquake was in the central fault zones of the Longmenshan tectonic belt. Average co-seismic vertical displacements along the individual fault of the Yingxiu-Beiehuan rupture zone reach 2.514 m and the cumulative vertical displacements across the central and frontal Longmenshan fault belt is about 5-6 m. The surface rupture strength was reduced from north of Beichuan to Qingchuan County and shows 2-3 m dextral strike-slip component. The Wenchuan thrust-faulting earthquake is a manifestation of eastward growth of the Tibetan Plateau under the action of continuous convergence of the Indian and Eurasian continents.

Dextral-Slip Thrust Faulting and Seismic Events of the Ms 8.0 Wenchuan Earthquake,Longmenshan Mountains,Eastern Margin of the Tibetan Plateau

Dextral-slip thrust movement of the Songpan-Garze terrain over the Sichuan block caused the Ms 8.0 Wenchuan earthquake of May 12, 2008 and offset the Central Longmenshan Fault （CLF） along a distance of -250 km. Displacement along the CLF changes from Yingxiu to Qingchuan. The total oblique slip of up to 7.6 m in Yingxiu near the epicenter of the earthquake, decreases northeastward to 5.3 m, 6.6 m, 4.4 m, 2.5 m and 1.1 m in Hongkou, Beichuan, Pingtong, Nanba and Qingchuan, respectively. This offset apparently occurred during a sequence of four reported seismic events, EQ1-EQ4, which were identified by seismic inversion of the source mechanism. These events occurred in rapid succession as the fault break propagated northeastward during the earthquake. Variations in the plunge of slickensides along the CLF appear to match these events. The Mw 7.5 EQ1 event occurred during the first 0-10 s along the Yingxiu-Hongkou section of the CLF and is characterized by 1.7 m vertical slip and vertical slickensides. The Mw 8.0 EQ2 event, which occurred during the next 10-42 s along the Yingxiu-Yanziyan section of the CLF, is marked by major dextralslip with minor thrust and slickensides plunging 25°-35° southwestward. The Mw 7.5 EQ3 event occurred during the following 42-60 s and resulted in dextral-slip and slickensides plunging 10° southwestward in Beichuan and plunging 73° southwestward in Hongkou. The Mw 7.7 EQ4 event, which occurred during the final 60-95 s along the Beichuan-Qingchuan section of the CLF, is characterized by nearly equal values of dextral and vertical slips with slickensides plunging 45°-50° southwestward. These seismic events match and evidently controlled the concentrations of landslide dams caused by the Wenchuan earthquake in Longmenshan Mountains.

Analysis of the Wenchuan Ms8.0 Earthquake's Co-seismic Stress and Displacement Change by Using the Finite Element Method

The variation of in situ stress before and after earthquakes is an issue studied by geologists. In this paper, on the basis of the fault slip dislocation model of Wenchuan Ms8.0 earthquake, the changes of co-seismic displacement and the distribution functions of stress tensor around the Longmen Shan fault zone are calculated. The results show that the co-seismic maximum surface displacement is 4.9 m in the horizontal direction and 6.5 m in the vertical direction, which is almost consistent with the on-site survey and GPS observations. The co-seismic maximum horizontal stress in the hanging wall and footwall decreased sharply as the distance from the Longmen Shan fault zone increased. However, the vertical stress and minimum horizontal stress increased in the footwall and in some areas of the hanging wall. The study of the co-seismic displacement and stress was mainly focused on the long and narrow region along the Longmen Shan fault zone, which coincides with the distribution of the earthquake aftershocks. Therefore, the co-seismic stress only affects the aftershocks, and does not affect distant faults and seismic activities. The results are almost consistent with in situ stress measurements at the two sites before and after Wenchuan Ms8.0 earthquake. Along the fault plane, the co-seismic shear stress in the dip direction is larger than that in the strike direction, which indicates that the faulting mechanism of the Longmen Shan fault zone is a dominant thrust with minor strike-slipping. The results can be used as a reference value for future studies of earthquake mechanisms.

Ground Surface Ruptures and Near-Fault,Large-Scale Displacements Caused by the Wenchuan Ms8.0 Earthquake Derived from Pixel Offset Tracking on Synthetic Aperture Radar Images

The 12 May 2008 Wenchuan Ms8.0 earthquake produced surface displacements along the causative fault, the Yingxiu-Beichuan Fault, which are up to several meters near the fault. Because of the large gradient, satellite synthetic aperture radar （SAR） interferometric data are strongly incoherent; the usual SAR interferometry method does not allow such displacements to be measured. In the present study, we employed another approach, the technique based on pixel offset tracking, to solve this problem. The used image data of six tracks are from the Advanced Land Observing Satellite, Phased Array type L-band Synthetic Aperture Radar （ALOS/PALSAR） dataset of Japan. The results show that the entire surface rupture belt is 238 km long, extending almost linearly in a direction of 42°north-east. It is offset left laterally by a north-west-striking fault at Xiaoyudong, and turns at Gaochuan, where the rupture belt shifts toward the south by 5 km, largely keeping the original trend. In terms of the features of the rupture traces, the rupture belt can be divided into five sections and three types. Among them, the Beichuan-Chaping and Hongkou-Yingxiu sections are relatively complex, with large widths and variable traces along the trend. The Pingtong-Nanba and Qingping-Jingtang sections appear uniform, characterized by straight traces and small widths. West of Yingxiu, the rupture traces are not clear. North of the rupture belt, surface displacements are 2.95 m on average, mostly 2-3.5 m, with 7-9 m the maximum near Beichuan. South of the rupture belt, the average displacement is 1.75 m, dominated by 1-2 m, with 3-4 m at a few sites. In the north, the displacements in the radar line of sight are of subsidence, and in the south, they are uplifted, in accordance with a right-slip motion that moves the northern wall of the fault to the east, and the southern wall to the west, respectively. Along the Guanxian-Jiangyou Fault, there is a uplift zone in the radar line of sight, which is 66 km long, 1.5-6 km wide, and has vertical displacements of approximately 2 m, but no observable rupture traces.

Characteristics and factors that influenced damage to dams in the M_s 8.0 Wenchuan earthquake

Based on raw data from dams damaged in the Wenchuan earthquake, including many that were severely damaged, characteristics and factors that influenced the damage are discussed in this paper. Findings from this study include： severely damaged dams were densely distributed along the seismologic fault; small dams, especially small earth-rock dams, had the most serious damage that was caused by a variety of factors; the most serious damage was caused by seismic waves; damage was aggregated by aftershocks; and the extent of the damage patterns increased with the seismic intensity. Damage patterns varied in different intensity zones and cracking was the most common type of damage. Most of the dams had a good base with relatively high bearing capacity, and the walls of the earth-rock dams were mostly of clay soil. This type of base and body material mitigated some of the damage to dams. Reservoir maintenance and other factors also have a significant impact on the seismic safety of the dam. Finally, some recommendations to reduce seismic damage to dams are proposed.

Post-earthquake economic resilience and recovery efficiency in the border areas of the Tibetan Plateau: A case study of areas affected by the Wenchuan Ms 8.0 Earthquake in Sichuan, China in 2008

The border areas of the Tibetan Plateau and the neighboring mountainous areas have a high incidence of earthquakes with a magnitude greater than Ms 5.0, as well as having a dense distribution of geological disasters such as collapses, landslides, and debris flows. Revealing the post-disaster economic development and recovery process is very important for enhancing disaster prevention and response capacity in order to formulate control policies and recovery methods for post-disaster economic reconstruction based on economic resilience. Using long-term socioeconomic data and the autoregressive integrated moving average (ARIMA) model, this paper calculated the economic resilience index of the areas most severely affected by the Wenchuan Earthquake of 2008 and adopted the improved variable returns to scale (VRS) date envelopment analysis (DEA) model and the Malmquist productivity index to analyze the efficiency and effect of annual post-disaster recovery. The results show that: (1) the economic resilience index of the areas most severely affected by the Wenchuan Earthquake was 0.877. The earthquake resulted in a short-term economic recession in the affected areas, but the economy returned to pre-quake levels within two years. In addition, the industrial economy was less resilient than agriculture and the service industry. (2) The comprehensive economic recovery efficiency of the disaster-stricken area in the year following the disaster was 0.603. The comprehensive efficiency, the pure technical efficiency, and the scale efficiency of the plain and hilly areas were significantly greater than those of the plateau and mountain areas. (3) The annual fluctuation in total factor productivity (TFP) following the disaster was considerable, and the economic recovery efficiency decreased significantly, resulting in a short-term economic recession. The TFP index returned to steady state following decreases of 33.7% and 15.2%, respectively, in the two years following the disaster. (4) The significant decrease in the post-disaster recovery efficiency was caused mainly by technological changes, and the renewal of the production system was the leading factor in determining the economic resilience following the disaster. With the decline in the scale of economic recovery following the earthquake, long-term economic recovery in the disaster-stricken areas depended mainly on pure technical efficiency, and the improvement in the latter was the driving force for maintaining the long-term growth of the post-disaster economy. Therefore, according to the local characteristics of natural environment and economic system, the disaster-stricken areas need to actively change and readjust their economic structures. At the same time, attention should be paid to updating the production system to enhance the level of technological progress and give full play to the scale effects of large-scale capital, new facilities, human resources, and other investment factors following the disaster so as to enhance the impact of economic resilience and recovery efficiency in response to the disaster.

Dynamic changes of gravity fields before and after the 2008 Wenchuan earthquake(Ms8.0)

The pattern evolution and dynamic mechanism of the dynamic changes of regional gravity fields occurring before and after the Wenchuan Ms8.0 earthquake are analyzed, based on five epochs of 1998 -2007 mobile gravity data from the middle-south section of the north-south seismic belt, and two epochs of field research data collected after the 2008 Wenchuan earthquake in combination with GPS data, leveling observations, and geotectonic environment data. The regional dynamic gravity changes demonstrate the effects of the eastward flow of solid matter in the Qinghai-Tibetan plateau and the preparation of the 2008 Wenchuan earthquake （2- 10 yr）. The two most meaningful gravity indicators of the Wcnchuan earthquake preparation are the positive （increasing） gravity changes occurring over many years in the southwest epicenter and the largescale gradient zone of gravity variation, with the cumulative difference between the two sides of the gradient zone of gravity exceeding 200 μGal. The positive gravity changes may facilitate a constant energy accumulation and the gradient belt may support seismic shear breakage. Overall, the gravity changes associated with the earthquake preparation indicate a pattern of accelerating increase-decelerating increase-earthquake occurrence. The Songpan-Ganzi block generally displays a negative gravity change, providing evidence for a local upwarp- ing of the deep crust-mantle and an interior expansion of the deep crust attributable to high temperatures. The viewpoint is consistent with the dilatant mechanism for earthquake preparation.

Improvement of instantaneous positioning by spatial stacking and modified sidereal filtering methods with application to the 2008 Ms8.0 Wenchuan earthquake

We obtained several displacement time series from the Sichuan permanent GPS net and processed the 1-Hz data observed during a few days before the 2008 Ms8.0 Wenchuan earthquake by double-difference instantaneous positioning technique. We filtered the data by the spatial stacking and the modified sidereal filte- ring methods to reduce correlation bias in space and time. The results indicate that these methods can improve the precision significantly.

Characteristics of the May 12,2008 Wenchuan M_S8.0 Earthquake Sequence and Discussion on Relevant Problems

The aftershock activity of the May 12, 2008 Wenchuan Ms8.0 Earthquake Sequence shows an obvious segmented feature. Most of the large aftershocks were distributed in the north and south parts of the aftershock zone. Thrusting was dominant with a small amount of strike-slip component in the south part. The aftershock activity decayed gradually, presenting the sequence features of a mainshock-aftershock pattern. The north part was the ending area of the malnshock fracture where strike-slipping was dominant, showing an obvious swarm feature. Therefore it became the major area for large aftershocks. The modulation of the earth tide on aftershock activity is remarkable; most large aftershocks occur during the period of flood and neap tide. The time period around 16：00 was the dominant occurring time for large aftershocks. The p-value, a parameter of modified Omori formula, increases gradually with time, and reaches about 1 at the end. Based on previous study, the sequence patterns, magnitude of maximum aftershock, as well as the duration of aftershock activity has been discussed. The primary results also show that the magnitude difference between the maiushock and the maximum aftershock is proportional to the rupture size of the maiushock for huge earthquakes of about Ms8.0. This means that when the magnitudes of the earthquakes are nearly the same, large rupture size corresponds to sufficient energy release.

Chongqing Beiwenquan Hot Spring Variation Associated with the May 12 2008 Ms 8.0 Earthquake in Wenchuan,Sichuan

Chongqing Beiwenquan hot springs were continuously monitored from May 1 to July 4,2008,during the May 12,2008 Ms 8.0 earthquake.On the spot in these hot springs were measured during this time. Chongqing Beiwenquan hot spring about 400 km away from Wenchuan Sichuan,had a huge response to the earthquake.After the earthquake,three higher elevation springs are cutoff,and the other four springs’ discharge has been increased nearly 50%. Temperature of these springs all reduced 1℃. Concentration of K+,Na+,Ca2+ in these hot springs has been reduced,and has an excellent linear correlation.Concentration of Mg2+ has been increased,Concentration of SO42-,F- has

Characteristics of coseismic water level changes at Tangshan well for the Wenchuan MS_8.0 earthquake and its larger aftershocks

Coseismic water level changes which may have been induced by the Wenchuan Ms8.0 earthquake and its 15 larger aftershocks （Ms〉5.4） have been observed at Tangshan well. We analyze the correlation between coseismic parameters （maximum amplitude, duration, coseismic step and the time when the coseismic reach its maximum amplitude） and earthquake parameters （magnitude, well-epicenter distance and depth）, and then compare the time when the coseismic oscillation reaches its maximum amplitude with the seismogram from Douhe seismic station which is about 16.3 km away from Tangshan well. The analysis indicates that magnitude is the main factor influencing the induced coseismic water level changes, and that the well-epicenter distance and depth have less influence. Ms magnitude has the strongest correlation with the coseismic water level changes comparing to Mw and ML magnitudes. There exists strong correlation between the maximum amplitude, step size and the oscillation duration. The water level oscillation and step are both caused by dynamic strain sourcing from seismic waves. Most of the times when the oscillations reach their maximum amplitudes are between S and Rayleigh waves. The coseismic water level changes are due to the co-effect of seismic waves and hydro-geological environments.

Temporal and Spatial Evolution of Precursory Anomalies of the Wenchuan MS8.0 Earthquake and Their Mechanical Analysis

In order to study the spatiotemporal evolution of the precursory anomalies 10 years before the Wenchuan M_S8. 0 earthquake in 2008, the epicentral distance of the precursory anomalies is calculated by using the geometric center of the rupture region and the elliptical centerline of the aftershock region. The result shows, precursor anomalies gradually increased about 2 years before the Wenchuan earthquake. The ratio of abnormal items is greater than 25% in the near source area (about twice the source scale) and 17%-24% in the remote area (about 3-5 times the source scale). There are three different stages of spatiotemporal evolution of precursory anomalies. During the α stage (including α_1 and α_2,between 700 to 3000 days before the main earthquake),the anomalies are mainly distributed in the southwest and northwest area of the Wenchuan aftershocks area. It is shown that the precursors of the far source region and the near source area have the characteristics of outward expansion. During the β stage (between 300 to 700 days before the main earthquake), the anomalies are distributed in the southwest and northern region of the aftershock region, showing a large range of anomalies. During the γ stage (including γ_1 and γ_2, 300 days before the main earthquake),the range of anomaly distribution is wide,and the anomalies are distributed in the southwest and northeast of the aftershock area. The anomalies converged to epicenter (γ_1) in the far source region and expand outwards (γ_2) in the near source region. Results of the experimental study and mechanical analysis of earthquake preparation process indicate that the three-stage characteristics of precursory anomalies in the process of earthquake preparation may be controlled by the seismogenic body,which is a form of expression in the process of earthquake preparation and a universal featureduring the earthquake preparation process,which has a certain guiding role in earthquake prediction.

Post-seismic relaxation process and vertical deformation following the 2008 Ms8.0 Wenchuan earthquake, China

The post-seismic horizontal and vertical deformations following the 2008 Ms8.0 Wenchuan earth- quake are inferred from GPS and precise leveling data. The post-seismic relaxation process is measured using GPS data from campaign stations located around the Longmenshan fault, and the derived decay time constant is 12 days. The evolution of the post-seismic vertical deformation is obtained from precise leveling data measured near the surface rupture. The results demonstrate that the hanging wall is uplifting and the foot wall is subsi- ding. The amplitude of the post-seismic deformation is lower than that of the co-seismic deformation. The re- gion with the largest post-seismic displacement is located on the leveling route between Maoxian and Beichuan on the hanging wall.

The Energy Transformation and Efficiency of the May 12,2008 Ms 8.0 Wenchuan Earthquake

The energy transformation and efficiency is now a hot topic among researches of scientific drilling into fault zones （Tanaka et al., 2006; Ma et al., 2006）. This study conducted temperature measurements and fault gouge particle analysis of borehole WFSD-1 from the Wenchuan Earthquake Fault Science Drilling Project （WFSD）, and discussed the earthquake energy budget. The research progress is illuminated as follows.