Earthquake Damage Research Articles

Post-earthquake fire performance of partially encased composite steel and concrete columns

Considering the potential severity of losses caused by post-earthquake fires compared to earthquakes alone, this study employed a three-dimensional Finite Element (FE) model established in ABAQUS to investigate the response of partially encased composite (PEC) steel and concrete columns under post-earthquake fire, serving as preliminary analysis for subsequent experimental studies. The validity of the FE model was verified using results from existing seismic and fire loading tests. In simulating the hysteresis behavior of PEC columns, this study evaluated the impacts of concrete discrete cracking and mechanical models on the simulation outcomes, thereby refining the modeling technique. A typical PEC column was selected as the focus of this investigation. Sequential earthquake-fire simulations were conducted using a data transfer technique to assess the fire resistance performance of different seismic damaged columns with different axial compression load ratios. The results indicate that the behavior of damaged PEC columns evolves based on initial damage, and the failure exhibits overall bending accompanied by local buckling of the steel flange. As the axial load ratio increases, the fire resistance time of seismic-damaged PEC columns significantly decreases. When the axial compression ratio is increased from 0.2 to 0.5, the fire resistance time of columns with damage factor D within 0.6 decreases by 59.07% on average. This sensitivity to axial pressure is particularly pronounced when earthquake damage is severe. Moreover, it finds that the lateral drift ratio under seismic action should be controlled within 2% to ensure that the post-earthquake fire resistance will not be significantly reduced.

Post-earthquake damage probability assessment of high-speed railway simply supported bridge based on traffic capacity

Railway transportation is an important channel for emergency rescue and disaster relief. Traditionally, the method of evaluating the damage degree of bridges under different seismic intensities from the perspective of structural damage cannot be directly used to judge the traffic capacity of trains on bridges after earthquakes. It is impossible to guarantee the running safety and affect the transportation of materials after the earthquake. Because the damage state of track structure directly affects the running safety of trains, this study combines the relationship between the damage of track structure and the deformation of bridge structure after earthquake, and puts forward the transverse relative displacement of girder end (TRDGE) as the damage index of train running function of railway bridge after earthquake. Then, the OpenSees-TRBF (Train-Railway-Bridge-Foundation Coupled System Dynamic Analysis Platform) co-simulation method is used to carry out the train-track-bridge coupling dynamic calculation under different seismic damage states, and the running safety analysis database under different TRDGEs is established. Then, the Particle Swarm Optimization (PSO) algorithm is used to classify the bridge damage index of traffic capacity through different target safety speeds. Based on the probabilistic seismic demand model, the exceedance probability curve of traffic capacity damage under different seismic intensities is established. Finally, based on the curve, the running safety risk domain of trains on the bridge under different seismic intensities is obtained with a confidence level of 95 %. This study proposes a post-earthquake damage assessment method for railway bridges based on the assurance of running safety probabilities, evaluating the damage degree of the bridge with the residual traffic capacity of the bridge after the earthquake. This study can be applied to assess earthquake damage on existing railway bridges and provide references for the safe speed of trains on the bridge after an earthquake.

Comprehensive review of AI and ML tools for earthquake damage assessment and retrofitting strategies

Comprehensive review of AI and ML tools for earthquake damage assessment and retrofitting strategies

A novel tool for seismic response analysis of tunnel in multilayered media based on kinematic earthquake source

Earthquake damage of underground structure indicates that the seismic response of tunnels is affected by the earthquake source, propagation medium, and local tunnel-surrounding rock system. Current research on the seismic response of tunnels primarily focuses on local tunnel systems, posing a challenge in comprehensively assessing the impact of earthquakes on tunnel group. This study introduces a theoretical seismological approach for the seismic response analysis of tunnels and presents a novel methodology for assessing the response. The proposed approach uses a kinematic seismic source as the earthquake source, assuming the geological strata as a layered medium and a tunnel as a point within a semi-infinite space. The wavenumber-integral method was used to solve for the seismic wave field and strain. To reflect the followability of underground structure to the surrounding rock deformation, we improved it and derived an analytical formula for the strain. The 2022 Menyuan earthquake was considered as a case study; the seismic damage to tunnels along the Lan-Xin Railway was analyzed using the proposed approach. The results confirm that there is a strong correlation between tunnel damage and free-field strain induced by earthquakes. The theoretical influence range of earthquake on tunnel was investigated using the diameter deformation ratio of the tunnel as an indicator. For tunnels embedded in a grade III rock mass, the theoretical range was approximately 7.0 km on both sides of the fault. For tunnels embedded in a grade V rock mass, the theoretical range can reach up to 25.9 km. The results also revealed that tunnels on the hanging-wall side of the fault had a more significant impact than those on the foot-wall side, although their affected areas were smaller.

Large-Depth Ground-Penetrating Radar for Investigating Active Faults: The Case of the 2017 Casamicciola Fault System, Ischia Island (Italy)

We conducted large-depth Ground-Penetrating Radar investigations of the seismogenic Casamicciola fault system at the volcanic island of Ischia, with the aim of constraining the source characteristics of this active and capable fault system. On 21 August 2017, a shallow (hypocentral depth of 1.2 km), moderate (Md = 4.0) earthquake hit the island, causing severe damage and two fatalities. This was the first damaging earthquake recorded on the volcanic island of Ischia from the beginning of the instrumental era. Our survey was performed using the Loza low-frequency (15–25 MHz) GPR system calibrated by TDEM results. The data highlighted variations in the electromagnetic signal due to the presence of contacts, i.e., faults down to a depth larger than 100 m below the surface. These signal variations match with the position of the synthetic and antithetic active fault system bordering the Casamicciola Holocene graben. Our study highlights the importance of employing large-depth Ground-Penetrating Radar geophysical techniques for investigating active fault systems not only in their shallower parts, but also down to a few hundred meters’ depth, providing a contribution to the knowledge of seismic hazard studies on the island of Ischia and elsewhere.

Probabilistic damage evaluation of isolated steel box-girder bridge excited by near-field earthquakes

Modern civilization regards bridges as lifelines, and their failure due to significant earthquakes in the past has raised concerns about their seismic vulnerability. The effect of near-fault earthquakes on the bridge vulnerability is a significant problem that has not been well explored. The present work utilizes fragility analysis to assess the probabilistic seismic risk evaluation of a continuous isolated steel box girder bridge exposed to near-field ground motions (NFGMs). Probabilistic seismic demand models (PSDMs) of the bridge are explored to reliably associate key demand and intensity measures to capture the seismic performance of the system. This study examines a steel box girder bridge in Guwahati, India, isolated with lead rubber (LRB) bearings and a friction pendulum system (FPS). To achieve a comparable seismic response, each bearing system is designed for a similar isolation period. The PGV/PGA ratio is a useful metric for classifying ground motion suites and estimating their frequency content by classifying the input excitations into low- and high-frequency motions. To perform probabilistic seismic risk assessments, two ensembles of forward directivity records with PGV: PGA >160 cm/s/g (high frequency) and PGV: PGA <160 cm/s/g (low frequency) are considered. Moreover, one ensembles of near-field incorporating the fling effect and far-field are also utilized. The failure probability under the seismic hazards for isolated bridges is assessed utilizing different damage measures (DM) such as maximum displacement ductility of pier (DDP), displacement of deck (MDD), drift ratio of pier (PDR), and rotational ductility (RDP). The near-field earthquake damage probability is significantly impacted by the PGV/PGA ratio, as a high PGV/PGA ratio results in a greater probability of damage than a low PGV/PGA ratio. In conclusion, guidelines are provided that will help the design engineers for designing bridges with LRB and FPS isolators for seismic isolation in the event of varying ground motion characteristics.

Seismic Fragility Analysis of a Multi-Tower Super High-Rise Building Under Near-Fault Ground Motions

ABSTRACT Earthquake damage to engineering structures often occur in near-fault areas. Near-fault ground motion is characterized by high-energy velocity pulses, long action period, and strong destructiveness. This paper forces on the seismic analysis of a multi-tower super high-rise building that subjected to near-fault ground motions. The interactions between each tower makes the internal forces and dynamic responses much complex and totally different from single-tower super high-rise buildings. However, currently, multi-tower buildings in China are usually regarded as multiple single-tower buildings, and then separately designed based on the relevant provisions in seismic code for single-tower buildings. The mutual influence between the towers is almost not taken into account. Seismic response history analysis was carried on a three-tower super high-rise building using finite element method. Various near-fault and far-fault ground motions and high-energy velocity pulses were used as inputs. Through calculation of inter-story drift ratio and maximum displacement, the quantitative seismic fragility analysis was conducted for the main tower and sub-tower structures. The results show that the responses from near-fault ground motions to the structure are much greater than those from far-fault ground motions. The damage probabilities of various damage levels for near-fault ground motions are 12.5% higher than far-fault ground motions averagely. A further analysis reveals that the velocity pulses are the main factor causing structural damage in near-fault ground motions.

ESenTRy: an on-site earthquake early warning system based on the instrumental modified Mercalli intensity

Earthquake early warning systems (EEWSs) are real-time seismology-based applications that purpose to reduce earthquake damage by taking some precautions for environments that are likely to be damaged. These systems are divided into two categories as on-site and regional. In this study, an on-site EEWS called ESenTRy, consisting of an accelerometer, GPS and communication components, has been designed. In this scope, a Windows operating system-based software has been developed using the .NET Framework and C# language. The developed software has the ability to read the Güralp Compressed Format strong motion data in real-time, compute the strong motion indices such as destructive intensity, real-time intensity and instrumental modified Mercalli intensity, issue an alarm in the target area and alert the defined Modbus clients. To test the reliability of ESenTRy, two simulations have been performed using the strong motion data belonging to the magnitude 7.7 Kahramanmaraş-Pazarcık earthquake that occurred in Türkiye on February 6, 2023, at 01:17:32 UTC and magnitude 5.9 Düzce-Gölyaka earthquake that occurred in Türkiye on December 23, 2022, at 01:08:15 UTC. In the test applications, depending on the distance between the stations and earthquake epicenters, the time difference between the Level-2 alarm trigger times and S-wave arrival times varied between 1.5 and 16 s, while the time difference between the Level-2 alarm trigger times and P wave arrival times varied between 0.5 and 29.5 s. The flowchart diagrams of the developed algorithm and results obtained from the simulations are presented in detail.

Reappraising the seismogenic potential of a low-strain rate region: Active faulting in the eastern Siena Basin (southern Tuscany, Italy)

We investigated the active tectonics and earthquake potential of the eastern Siena Basin, a slowly deforming portion of southern Tuscany in the inner Northern Apennines. This region hosts several historical settlements and valuable cultural heritage, but also frequent background seismicity and rare damaging earthquakes in the Mw range 5.0–6.2. We describe in detail an active, capable, and seismogenic fault system that we identified in the eastern Siena Basin, a few kilometers south-east of the city of Siena, thanks to the presence of an active quarry (Cava Capanni) that exploits travertines of Middle Pleistocene-Holocene age. Travertines are unique rock masses that may preserve living evidence of active and seismogenic faulting, thus providing remarkable seismotectonic insight. The active fault system consists of at least two segments rupturing travertines younger than 45 ka, with a cumulative vertical displacement of 111 cm, and an estimated minimum slip rate of 0.02–0.03 mm/y. We maintain that this displacement is the result of at least three coseismic movements accompanied by clastic dykes injected within the fault damage zone due to liquefaction phenomena. The fault system is seen to extend east of the quarry, affecting Pliocene and Mesozoic deposits.The Cava Capanni fault system is evidence of a poorly understood but potentially seismogenic tectonic mechanism of regional extent. Its orientation and kinematics are compatible with the activity of faults that are oriented obliquely or orthogonally to the main chain axis, in contrast with the setting of the axial and outer zone of the Northern Apennines, where extension and compression are accommodated by Apennines-parallel faults.

Multi-Source Feature-Fusion Method for the Seismic Data of Cultural Relics Based on Deep Learning.

The museum system is exposed to a high risk of seismic hazards. However, it is difficult to carry out seismic hazard prevention to protect cultural relics in collections due to the lack of real data and diverse types of seismic hazards. To address this problem, we developed a deep-learning-based multi-source feature-fusion method to assess the data on seismic damage caused by collected cultural relics. Firstly, a multi-source data-processing strategy was developed according to the needs of seismic impact analysis of the cultural relics in the collection, and a seismic event-ontology model of cultural relics was constructed. Additionally, a seismic damage data-classification acquisition method and empirical calculation model were designed. Secondly, we proposed a deep learning-based multi-source feature-fusion matching method for cultural relics. By constructing a damage state assessment model of cultural relics using superpixel map convolutional fusion and an automatic data-matching model, the quality and processing efficiency of seismic damage data of the cultural relics in the collection were improved. Finally, we formed a dataset oriented to the seismic damage risk analysis of the cultural relics in the collection. The experimental results show that the accuracy of this method reaches 93.6%, and the accuracy of cultural relics label matching is as high as 82.6% compared with many kinds of earthquake damage state assessment models. This method can provide more accurate and efficient data support, along with a scientific basis for subsequent research on the impact analysis of seismic damage to cultural relics in collections.

Comprehensive Risk Assessment of Farmer Households in Flood and Earthquake Prone Areas Based on Multidimensional Information Diffusion Model

In recent years, the global practice of regional natural disaster reduction and comprehensive risk prevention has promoted the integration of multiple disciplines, and the theory and technology of comprehensive risk assessment for multiple natural disaster chains urgently need innovation. This article is based on multiple natural disaster chain scenarios and sets secondary indicators for multiple influencing factors based on the first level indicators composed of basic risk factors such as hazard factors, vulnerability of disaster bearing bodies, and disaster prevention and reduction capabilities. By using multidimensional information diffusion technology, multiple hazard factor hazard combination factors, vulnerability combination factors of disaster bearing bodies, and comprehensive disaster prevention and reduction capability combination factors reflecting various natural disaster scenarios are constructed, A household comprehensive risk assessment index system for multi-dimensional information diffusion assessment of natural disaster comprehensive risk was designed by quantifying risk factors such as the intensity of comprehensive disaster causing factors, vulnerability of disaster bearing bodies, and comprehensive disaster prevention and reduction capabilities in disaster scenarios. This paper expands the traditional information diffusion comprehensive risk assessment model, makes up for the lack of information on geographical units, establishes a multi impact factor multi dimension natural disaster risk comprehensive assessment model for multi risk factors of complex disaster scenarios with small samples, obtains the theoretical value of multi-dimensional information Diffusion model comprehensive risk assessment for complex disaster scenarios, breaks through the sample size limit of natural disaster risks, and makes up for the lack of information on small samples of small probability events, Improved the reliability and accuracy of comprehensive risk assessment for multiple natural disaster chains. Then, take the disaster scenario formed by the comprehensive risk of earthquake and flood in a place in the upper reaches of the Fujiang River basin as the analysis sample, calculate the comprehensive risk assessment value of flood and earthquake damage in the study area's micro samples (rural households in villages and towns). By comparison, it is found that the comprehensive risk value assessed by the multi-dimensional information Diffusion model has good reliability and high accuracy, which is significantly better than the traditional Kriging interpolation method and geographical weighted regression method. This article provides a comprehensive risk assessment of natural disasters at the micro household level, which can better reflect the diversity and regionality of risks, help clarify the hidden risks of rural natural disasters, and provide scientific basis for government emergency management decision-making.

The impact of Kahramanmaraş (2023) earthquakes: A comparative case study for Adıyaman and Malatya

This study examines the effects of two major earthquakes of magnitude 7.7 and 7.6 that struck Kahramanmaraş on February 6th, 2023, followed by a magnitude 6.4 quake in Hatay on February 20th, which caused major damage in 11 Turkish provinces. The study focuses on Adıyaman and Malatya and uses an interdisciplinary approach to analyze the economic and environmental impacts. Primary data sources, including field visits and interviews, reveal clear labor-related challenges in both provinces, characterized by a government-induced labor shortage. In both provinces, physical capital has been severely damaged, particularly affecting small businesses, historic bazaars, and old industrial areas. The impact on businesses varies by size and location, with Adıyaman suffering more severe setbacks than other cities. The shortage of skilled labor related to the earthquake damage affects the quality of production, which can have a serious economic impact. Transportation disruptions continue to hamper supply chains and affect companies' ability to meet their export commitments. The environmental consequences, particularly the large amount of debris, pose a major challenge. The lack of a comprehensive disaster waste plan at the central government level leads to inadequate waste management. The study recommends sorting the debris at temporary sites to obtain reusable items while paying attention to the sustainability and transparency of debris management processes. In summary, this comparative case study highlights the need for tailored approaches to address the different impacts in the 11 provinces. A one-size-fits-all solution is insufficient and an individual needs assessment is needed for each province in order to implement targeted economic and environmental recovery measures.

Multi-Objective Optimization in Support of Life-Cycle Cost-Performance-Based Design of Reinforced Concrete Structures

Surveys on the optimum seismic design of structures reveal that many investigations focus on minimizing initial costs while satisfying performance constraints. Although reducing initial costs while complying with earthquake design codes significantly ensures occupant safety, it may still cause considerable economic losses and fatalities. Therefore, calculating potential earthquake damages over the structure’s lifetime is essential from an optimal Life-Cycle Cost (LCC) design perspective. LCC analysis evaluates economic feasibility, including construction, operation, occupancy, maintenance, and end-of-life costs. The population-based, meta-heuristic Ideal Gas Molecular Movement (IGMM) algorithm has proven effective in solving highly nonlinear mono- and multi-objective engineering problems. This paper investigates the LCC-based mono- and multi-objective optimum design of a 3D four-story concrete building structure using the Endurance Time (ET) method, which is employed for its efficiency in estimating structural responses under varying seismic hazard levels. The novelty of this work lies in integrating the ET method with the IGMM algorithm to comprehensively address both economic and performance criteria in seismic design. The results indicate that the proposed technique significantly reduces minor injury costs, rental costs, and income costs by 22%, 16%, and 16%, respectively, achieving a total reduction of 10% in all structural Life-Cycle Costs, which is considered significant.

Multi-Scale Earthquake Damaged Building Feature Set

Earthquake disasters are marked by their unpredictability and potential for extreme destructiveness. Accurate information on building damage, captured in post-earthquake remote sensing images, is critical for an effective post-disaster emergency response. The foundational features within these images are essential for the accurate extraction of building damage data following seismic events. Presently, the availability of publicly accessible datasets tailored specifically to earthquake-damaged buildings is limited, and existing collections of post-earthquake building damage characteristics are insufficient. To address this gap and foster research advancement in this domain, this paper introduces a new, large-scale, publicly available dataset named the Major Earthquake Damage Building Feature Set (MEDBFS). This dataset comprises image data sourced from five significant global earthquakes and captured by various optical remote sensing satellites, featuring diverse scale characteristics and multiple spatial resolutions. It includes over 7000 images of buildings pre- and post-disaster, each subjected to stringent quality control and expert validation. The images are categorized into three primary groups: intact/slightly damaged, severely damaged, and completely collapsed. This paper develops a comprehensive feature set encompassing five dimensions: spectral, texture, edge detection, building index, and temporal sequencing, resulting in 16 distinct classes of feature images. This dataset is poised to significantly enhance the capabilities for data-driven identification and analysis of earthquake-induced building damage, thereby supporting the advancement of scientific and technological efforts for emergency earthquake response.

Agent-based post-earthquake evacuation simulation to enhance early-stage architectural layout and non-structural design

In the indoor design process, architects make crucial decisions regarding architectural layout and the selection of non-structural elements. However, there is a lack of comprehensive consideration for human evacuation behavior, specifically in the event of earthquake evacuation, during the design process. This paper bridges this gap by presenting the application of Agent-Based Building Earthquake Evacuation Simulation (AB2E2S). The paper assesses a post-earthquake evacuation simulation prototype, which integrates an agent-based simulation technique with probabilistic earthquake damage assessment. The model is applied to the case of the engineering building at The University of Auckland, to evaluate the impact of earthquake intensity, design, and behavioral variables on Safe Evacuation Time and number of casualties. Overall, this paper demonstrates the potential of the AB2E2S prototype to inform architects and designers about the effective selection of non-structural elements and architectural layout scenarios for post-earthquake evacuation during the schematic design process.

Building information extraction and earthquake damage prediction in an old urban area based on UAV oblique photogrammetry

Building information extraction and earthquake damage prediction in an old urban area based on UAV oblique photogrammetry

Seismic Features Predict Ground Motions During Repeating Caldera Collapse Sequence

AbstractApplying machine learning to continuous acoustic emissions, signals previously deemed noise, from laboratory faults and slowly slipping subduction‐zone faults, demonstrates hidden signatures are emitted that describe physical details, including fault displacement and friction. However, no evidence currently exists to demonstrate that similar hidden signals occur during seismogenic stick‐slip on earthquake faults—the damaging earthquakes of most societal interest. We show that continuous seismic emissions emitted during the 2018 multi‐month caldera collapse sequence at the Kı̄lauea volcano in Hawai'i contain hidden signatures characterizing the earthquake cycle. Multi‐spectral data features extracted from 30 s intervals of the continuous seismic emission are used to train a gradient boosted tree regression model to predict the GNSS‐derived contemporaneous surface displacement and time‐to‐failure of the upcoming collapse event. This striking result suggests that at least some faults emit such signals and provide a potential path to characterizing the instantaneous and future behavior of earthquake faults.

Effects of passive confinement on out-of-plane robustness of unreinforced masonry infill walls

Unreinforced masonry (URM) walls are vulnerable to out-of-plane (OOP) collapse under lateral loads; however, URM structures are still popular in developing countries. Even in many countries classified as developed by the United Nations, URM walls are still commonly applied as exterior and interior infills in engineered buildings, resulting in severe OOP failure during strong earthquakes. The objectives of this study were: 1) to experimentally examine the OOP behavior and performance of URM infills with construction deficiencies commonly found in Asian countries, 2) to investigate the impact of construction quality on resistance to OOP collapse and 3) to propose a promising approach using passive confinement to improve the lateral resistance of URM infills and a method to evaluate their seismic performance. To attain these objectives, the study used two full-scale, single wythe, brick masonry specimens with specific boundary conditions at the top, a control specimen representing URM infills with construction flaws frequently found in Asian countries and a well-constructed specimen subjected to passive confinement under lateral excitation. A series of incremental dynamic shaking table tests was performed on both specimens under OOP unidirectional one-way bending. The test results revealed the positive effects of construction with passive confinement on OOP performance. The results indicated that the OOP robustness of infills can be improved without specific retrofit materials, leading to advantages such as simplicity, cost-effectiveness, and swift execution. Moreover, a theoretical method to evaluate the lateral resistance of infills was presented and verified for practical applications. This method can contribute to reducing the potential for earthquake damage in engineered buildings.

Seismic vulnerability distribution of the earthquake prone area in Central Aceh

Abstract Central Aceh is known as an earthquake-prone area due to the presence of the Sumatra Fault and secondary fault crossing the area resulting the M 6.1 earthquake in 2013. The objective of the research is to investigate seismic vulnerability of Central Aceh. Microtremor observations were conducted at 25 locations in Central Aceh analysed using the Horizontal to Vertical Spectral Ratio (HVSR). Seismic amplification is represented by the maximum values of H/V. The dominant frequency f0 representing the sediment thickness was obtained when the H/V curve is maximum. The frequency values ranged from 14.9 to 25.9 Hz, and the seismic vulnerability values varied from 0.02 to 3.1. The top soil consists of thin layers of sedimentary material dominated by tertiary rocks. The rocks are classified as Type IV and Class I based on the dominant frequency values. This result is consistent with the geology, consisting of sedimentary and volcanic rocks. Areas with a seismic vulnerability value &gt; 1.0 is considered to have a high risk of earthquake damage because they are likely to experience earthquake amplification. Areas with seismic vulnerability value &gt; 1.0 are found in West Aceh and the Geureudong volcano area.

Analysis of the Stress Field Around Concealed Active Fault From Minor Faults‐Slip Data Collected by Geological Survey: An Example in the 1984 Western Nagano Earthquake Region

AbstractEarthquakes with magnitudes of 6–7 have been reported even in various active tectonic settings where fault deformation topography have not been detected. Therefore, delineating concealed active faults generating such earthquakes is necessary to reduce earthquake damage; however, few studies exist to provide its clues regarding such faults. The 1984 Western Nagano Earthquake in Japan was a main shock with a magnitude of Mj 6.8 and a depth of 2 km at the source. Solid bedrocks are well‐exposed in the earthquake source region; however, no surface rupture have been identified, and the active fault is known to be concealed. In this study, we collected data on striations observed in fractures by geological survey around the source area of the 1984 Western Nagano Earthquake. Using the collected data, the multiple inverse method was used to estimate the stresses that affected the striation formation. Consequently, stresses similar to acting faults in this area were detected in minor faults around the known concealed active fault. This suggests that the minor faults might be part of the damage zone that has been developed around the concealed active fault. Some minor faults were recognized in Quaternary volcanic rocks, confirming that they experienced displacements recently. This study indicates the possibility of detecting concealed active faults in the bedrock by geological survey.