Research
Research
📚RPNet: Robust P-wave first-motion polarity determination using deep learning (Han et al., 2025; SRL)
We present RPNet, a robust deep-learning model for P-wave first-motion polarity determination aimed at deriving earthquake focal mechanism solutions. RPNet integrates advanced deep learning techniques, including inception modules, attention mechanisms, and Monte Carlo dropout, to improve prediction accuracy and stability with quantified uncertainty. We conduct benchmark tests against four existing models using datasets from the western United States and Italy. The results address key issues such as sensitivity to misaligned P-wave onsets, class imbalance between up and down polarities, and handling of “unknown” labels. To enhance robustness, we apply intensive random time shifts of P-wave onsets within a ±0.5-second range during training. To address the class imbalance, RPNet is trained with an equal amount of data for both up and down polarities. The performance of RPNet is evaluated through three tests. First, we test RPNet on a dataset comprising 10% of the data excluded from training, achieving 99% recall for both up and down polarities. Next, we apply RPNet to independent Hi-net data from Japan, where it demonstrates superior generalization compared to previous deep-learning models, achieving recall rates above 97.5% even under time shifts of up to ±0.40 seconds. Finally, RPNet is tested on the 2016 Kumamoto earthquake sequence, where its automatically derived focal mechanism solutions closely match those manually derived by the Japan Meteorological Agency, outperforming the previous model in Kagan angle and polarity misfit. These results highlight RPNet’s potential as a reliable tool for automating focal mechanism derivation across diverse regions and waveform conditions without the need for additional optimization.
Reference
Han, J., S, Kim, & D.-H. Sheen (2025), RPNet: Robust P-wave first-motion polarity determination using deep learning, Seismological Research Letters.
Data and Resources
RPNet package
📚Research Earthquake Catalog of the Southern Korean Peninsula (Han et al., 2024; SRL)
A seismicity catalog spanning 2012–2021 is proposed for the inland and near‐coastal areas of the southern Korean Peninsula (SKP). Using deep learning (DL) techniques combined with conventional methods, we developed an integrated framework for compiling a comprehensive seismicity catalog. The proposed DL‐based framework allowed us to process, within a week, a large volume of data (spanning 10 yr) collected from more than 300 seismic stations. To improve the framework’s performance, a DL picker (i.e., EQTransformer) was retrained using the local datasets from the SKP combined with globally obtained data. A total of 66,858 events were detected by phase association using a machine learning algorithm, and a DL‐based event discrimination model classified 29,371 events as natural earthquakes. We estimate source information more precisely using newly updated parameters for locations (a 1D velocity model and station corrections related to the location process) and magnitudes (a local magnitude equation) based on data derived from the application of the DL picker. Compared with a previous catalog, the proposed catalog exhibited improved statistical completeness, detecting 21,475 additional earthquakes. With the newly detected and located earthquakes, we observed the relative low seismicity in the northern SKP, and the linear trends of earthquakes striking northeast–southwest (NE–SW) and northwest–southeast (NW–SE) with a near‐right angle between them. In particular, the NE–SW trend corresponds to boundaries of major tectonic regions in the SKP that potentially indicates the development of fault structures along the boundaries. The two predominant trends slightly differ to the suggested optimal fault orientations, implying more complex processes of preexisting geological structures. This study demonstrates the effectiveness of the DL‐based framework in analyzing large datasets and detecting many microearthquakes in seismically inactive regions, which will advance our understanding of seismotectonics and seismic hazards in stable continental regions.
References
Han, J., Joo Seo, K., Kim, S., Sheen, D. H., Lee, D., & Byun, A. H. (2024). Research Catalog of Inland Seismicity in the Southern Korean Peninsula from 2012 to 2021 Using Deep Learning Techniques. Seismological Research Letters, 95(2A), 952-968.
Han, J., Seo, K. J., Byun, A. H., Kim, S., Sheen, D. H., & Lee, D. (2023, December) Deep Learning-Based Framework toward Comprehensive Seismic Catalogs: Optimization on the Southern Korean Peninsula. In AGU Fall Meeting Abstracts (S31G-0418).
Data and Resources
Earthquake catalog & Deep learning model file
📚Seismicity near the eastern Denali fault (Han et al., 2024; YGS)
We studied earthquakes near Burwash Landing, Yukon. Using data from temporary and permanent seismic stations, we enhanced the understanding of both regional and local earthquakes. The study used deep learning and template matching to effectively detect earthquakes, even from noisy data. Following detection, seismic parameters, earthquake location, and magnitude were estimated and refined. The analysis revealed 103 local earthquakes, with 28 located in an area of geothermal resource potential. Notable small-magnitude earthquakes were observed near Bock’s Creek fault. No earthquakes were observed on the Denali fault during the study period. The existence of active faults strike-parallel to the Denali fault suggests that local permeable structures may exist in the area. Regional observations detected 46 432 regional earthquakes in 13 years, but none along a section of the Denali fault near Burwash Landing, Yukon, which we interpret as a seismic gap.
References
Han, J., Dettmer, J., Gosselin, J. M., Gilbert, H., Biegel, K., & Kim, S. (2024). Seismicity near the eastern Denali fault from temporary and long-term seismic recordings. In: Yukon Exploration and Geology Technical Papers 2023, L.H. Weston and Purple Rock Inc. (eds.), Yukon Geological Survey,p. 37–50.
📚Temporary seismic arrays for monitoring volcanoes in the South Korea (Han et al., 2024; Geosci. J.)
Temporary seismic networks on Mt. Halla and Ulleung Island volcanoes were deployed, which employ broadband and geophone arrays to monitor potential volcanic activities and to estimate high-resolution magmatic structures beneath these volcanoes. The purpose of this paper is to introduce these networks and present early results through basic seismic analyses, suggesting the potential for future comprehensive seismological studies. The array in Mt. Halla volcano consists of five broadband sensors (JH array), and it has been operational around the Baengnokdam summit crater since October 2020. There was an additional linear geophone array (HL array) installed in September 2021 for detailed shallow subsurface imaging. Ulleung Island volcano had been under observation for two years since June 2021 with a network of nine broadband sensors (UL array) along its coast and in the Nari crater basin, complemented by a 52-geophone array (UG array) deployed in May 2022 for high-resolution subsurface studies. Despite the noisy environments typical of temporary setups, power spectral density analyses confirmed the quality of data as comparable to established reference noise models in permanent stations. Our study aimed to initiate studies uncovering seismic activities and structures beneath Mt. Halla and Ulleung Island volcanoes, specifically regarding volcanic activity. This approach detected no clear sign of volcanic seismicity on both islands, suggesting a period of magmatic dormancy. Seismic velocity variation (dv/v) analyses further indicated that environmental factors, rather than volcanic processes, influenced the changes in the physical properties of the underground structures. Conversely, the receiver function analysis and ambient noise data processing hinted at the presence of complex subsurface structures, potentially indicative of volcanic features, such as partial melting. Despite the lack of direct evidence for active magmatic processes, the collected seismic data provides a crucial baseline for future monitoring and a deeper understanding of the magmatic and tectonic dynamics beneath these volcanoes, offering valuable insights for ongoing volcanic research.
References
Han, J., Han, J., Heo, D., Kim, S., Lee, S., Koh, M. H., ... & Ahn, U. (2024). Data and early results from temporary seismic arrays for monitoring and investigating magmatic processes beneath Mt. Halla and Ulleung Island volcanoes, South Korea. Geosciences Journal, 28(5), 761-780.
📚Seismic event and phase detection using deep learning (Han et al., 2023; Geosci. J.)
Deep learning (DL) methods have a high potential for earthquake detection applications because of their high efficiency at processing measurement data, such as picking seismic phases. However, the performance of DL methods must be evaluated to ensure that they can replace conventional methods so that full automation can be achieved. State-of-art DL methods incorporate advanced techniques and train with large global datasets to enhance their earthquake detection capabilities. In this study, we tested a representative DL model on the 2016 Gyeongju earthquake sequence in the Korean Peninsula and compared the results with a previously established catalog and with the results of the conventional Short Time Average/Long Time Average (STA/LTA) method. The DL model demonstrated reasonable improvements in efficiency and performance by detecting more and smaller earthquakes within a much shorter running time than the other methods. In addition, the DL algorithms generally provided precise pickings of P- and S-wave phases. The DL model showed good generalization because it appropriately detected earthquakes in the study area that were not included in the training dataset. However, our results did suggest possible errors that should be accounted for, such as inconsistent phase picking, missing large earthquakes, and detecting non-natural earthquake signals. From the result of tests, local optimization may be important for realizing fully automatic earthquake monitoring, such as retraining with a local dataset, fine-tuning, or transfer learning. In addition, incorporating post-processing techniques such as phase association and discrimination into the DL framework is necessary.
References
Han, J., Kim, S., Sheen, D. H., Lee, D., Lee, S. J., Yoo, S. H., & Park, D. (2023). Seismic event and phase detection using deep learning for the 2016 Gyeongju earthquake sequence. Geosciences Journal, 27(3), 285-295.
On 3 May 2020, an ML 3.1 earthquake occurred in Haenam, southwestern Korea, in an area devoid of recorded seismicity since instrumental observations began in 1978. Careful examination of the temporal occurrence of seismicity, and the magnitude distribution of the sequence before and after the ML 3.1 earthquake, indicates typical swarm-like behavior. The earthquake swarm started with an ML 0.6 event on 26 April 2020, intensified up to 3 May 2020, and abruptly terminated with an ML 1.0 event on 9 May 2020. The Pusan National University Geophysics Laboratory (PNUGL) deployed a temporary seismic array with eight three-component short-period instruments to monitor the short-lived bursts of seismicity. During the monitoring campaign, we detected > 700 microearthquakes by applying a matched-filter technique to the combined dataset produced by PNUGL, the Korea Meteorological Administration, and the Korea Institute of Ocean Science and Technology. We determined earthquake parameters for 299 earthquakes that were detected at four or more seismic stations. We also determined the focal mechanism solutions of the 10 largest earthquakes in the swarm using first-motion polarities with S/P ratios. The focal mechanism, hypocentral depth, and stress orientation of the largest earthquake in the sequence were also determined using waveform inversions. The distribution of earthquake hypocenters, together with focal mechanism solutions, indicates that the earthquake swarm activated deeply-buried faults (~20 km) oriented either NNE-SSW or WNW-ESE. We also report details of the temporary seismic monitoring network, including the instrumentation, detection of microearthquakes, and variations in event-detection threshold influenced by anthropogenic and natural noise fluctuations. We also discuss the limitations associated with lowering the detection threshold of microearthquakes by increasing the number of seismic stations or by adopting advanced event-detection techniques.
References
Han, J., Seo, W., Kim, H. J., Kim, W. Y., Won, D., Chung, J. I., & Kim, K. H. (2020). Monitoring a short-lived earthquake swarm during April–May 2020 in Haenam, Korea, and its preliminary results. Geosciences Journal, 25, 43-57.
The magnitude 5.4 earthquake occurred at the northern part of Pohang city on 15 November 2017. Pusan National University operated temporary seismic stations to monitor the aftershock of Gyeongju earthquake which occurred on 12 September 2016 and found there was microseismic activity near the Pohang city, therefore additional stations were deployed 7 days before the Pohang earthquake. During the operation, totally 3,211 earthquakes including mainshock were detected and located. In this study, using the S wave arrival time of 12 stations which can be corrected the seismometer’s orientation and POHB which was operated by KIGAM(Korean Institute of Geoscience and Mineral Resources), the shear wave anisotropy in Pohang area was investigated through shear wave splitting method. Through eigenvalue method and cluster analysis method, fast wave direction and delayed time were determined and totally 6,231 splitting results were obtained. Before analysis the spatial variation, to find temporal variation of anisotropy, the station PHA2 which was operated by KMA(Korea Meteorological Administration) from January 2009 to November 2017(before the mainshock) was used for analysis shear wave anisotropy. As a result, the anisotropy which is related to the geothermal power plant was began to be observed after the mud loss and first stimulation period. Also, after the second stimulation, the anisotropy which is related to fault structure was began to be observed. Then, to find out the spatial patterns of shear wave anisotropy, the spatial average of fast wave direction was calculated using the TESSA. As a result, the patterns of anisotropy were mainly divided into two trends of directions, N30°E and E-W. In the former case, fast wave direction matches with the trend of fault, therefore anisotropy is related to structure. And the latter case, it matches with the direction of maximum horizontal stress, therefore anisotropy is related to the stress around this area.
References
Han, J., (2020). Study on Shear Wave Anisotropy of Pohang Area (Master’s thesis). Pusan National University.