Seismic Catalog Research Articles

Exploring Data Set Bias and Decision Support with Predictive Uncertainty Through Bayesian Approximations and Convolutional Neural Networks

Abstract Individual seismic catalogs can contain multiscale observations from fault level to global scales and associated waveforms from discrete events reflect crustal structure across many different scales and locations. Seismic network aperture, geographic location, and observation distance may not provide informative guidance or intuition on how different catalogs will behave across models trained under different conditions. We rely on uncertainty to provide guardrails for when to trust model decisions, but understanding when our uncertainty is trustworthy is an open challenge. Here, we explore Bayesian approximation methods for assigning predictive uncertainty in seismic event classification problems. We find that computationally expensive Bayesian approximations do not outperform simple ensemble methods. We also find that when exploiting multiple seismic event catalogs, joint training with data from all the catalogs combined with Bayesian approximations and supervised training for classification can obscure bias and result in less robust uncertainty while also not providing substantial performance benefits compared to training individual models for each catalog.

Intraplate seismicity in southwestern Norway: Enhanced catalogue highlights diffusive earthquake occurrence linked to inherited weakness zones

Summary Intraplate earthquakes in stable continental regions exhibit diverse characteristics in terms of timing, spatial distribution, and magnitude. They are often unexpected, and their underlying physical mechanisms are not well understood. This complexity is particularly apparent in Norway, where seismicity is mostly localised on the continental margin and coastal areas. Various studies have attempted to explain the causes of seismicity in Norway by invoking different sources of stress, ranging from regional stress due to ridge push to local effects such as topography or deglaciation. In this study, we revisit these questions by investigating the distribution of seismicity in southwestern Norway using an enhanced earthquake catalogue. To achieve this, we revised the Norwegian National Seismic Network seismic catalogue from 2000 to 2023 and built a new catalogue using machine-learning-based techniques on data from a temporary seismic deployment in the region. Thanks to the increased station density during this deployment, we were also able to calculate new fault plane solutions that consistently showed a WNW-ESE direction for the most compressive axis. Furthermore, we demonstrate that seismicity in southwestern Norway, while diffuse, tends to be localised around the major crustal shear zones of the region, such as the Bergen Arc Shear Zone and the Hardangerfjord Shear Zone.

Seismic and GNSS strain-based probabilistic seismic hazard evaluation for northern Chile using DAS Magnitude Scale

BackgroundProbabilistic Seismic Hazard Assessment (PSHA) is a leading methodology for determining key ground motion parameters such as Peak Ground Acceleration (PGA) and spectral acceleration (SA), essential for structural design. This approach uses extensive earthquake data, typically spanning over a century, leveraging frequency and magnitude statistics. However, long-term ground shaking probabilities may not always be accurately captured by traditional data-driven methods. To address these limitations, this study develops a PSHA map for Northern Chile using both seismic and GNSS (Global Navigation Satellite System) data. A curated homogeneous earthquake catalog, based on the advanced seismic moment magnitude scale Mwg(Das Magnitude Scale), replaces the traditional Mw scale to ensure superior accuracy, particularly for intermediate and smaller earthquakes.ResultsUsing the earthquake catalog, seismicity parameters ‘a’ and ‘b’ from the Gutenberg-Richter relationship were derived. Seismogenic modeling and Ground Motion Models (GMMs) were applied to estimate ground motion probabilities for a 475-year return period. Additionally, a PSHA map was constructed using GNSS strain rates, translating velocity-derived strain rates into seismic moment rates and ground shaking probabilities for seismic source zones. Comparative analyses revealed higher PGA values from GNSS strain data compared to seismic catalog data. GNSS strain data proved invaluable for refining seismic segmentation in Northern Chile, enhancing the precision of PSHA calculations.ConclusionsA PSHA map for Northern Chile, synthesizing seismic catalog data and GNSS strain rates using a Logic Tree-based algorithm, has been developed for a 475-year return period. This map provides a critical tool for generating seismic hazard assessments aligned with building codes and emergency planning protocols. By integrating GNSS strain rates and seismic data, this study advances the reliability and accuracy of long-term ground shaking predictions.

Does Subsurface Fault Geometry Affect Aleatory Variability in Modeled Strike-Slip Fault Behavior?

ABSTRACT The nonplanar geometry of faults influences their seismotectonic behavior, affecting the initiation, propagation, and termination of individual earthquakes as well as the stress–slip relationship and probability of multisegment rupture. Consequently, computer simulations that aim to resolve the earthquake rupture process and make predictions about a fault’s future behavior should incorporate nonplanar fault geometries. Although surface traces of faults can be mapped with high accuracy, a key challenge is to define a fault’s detailed subsurface geometry due to a general lack of data. This raises the question of which geometry to use. Does it matter which subsurface geometry is utilized in earthquake rupture simulations, as long as at least the fault trace is considered? How different is the simulated fault behavior for faults that share the same surface trace but different subsurface geometries? Using the physics-based earthquake-cycle simulator MCQsim, we generate seismic catalogs for 100 × 20 km strike-slip faults, assuming variations in fault surface trace, subsurface geometry, and strength distribution. We investigate how the long-term fault behavior—in the form of magnitude–frequency distribution, earthquake interevent time, and maximum earthquake size—is affected by fault geometry and strength distribution. We find that the simulated behavior of strike-slip faults with identical fault traces is interchangeable—even if their subsurface fault geometries differ. Implementing the fault trace constrains possible fault geometries to a level that makes the long-term behavior indistinguishable—at least for strike-slip faults with “known” strength distribution and length-to-width aspect ratios that are equal or larger than what we used here. The fault trace can provide a satisfactory representation of subsurface geometry for assessing long-term fault behavior.

Seismic velocity changes from repetitive seismicity at Mauna Loa prior to and during its 2022 eruption

AbstracMauna Loa’s short-lived eruption from late November to early December 2022 marked the culmination of nearly a decade of elevated seismic activity and geodetic inflation. The volcano has been monitored by a network of permanent, short period and broadband seismometers. I used the continuous waveform data from that network starting in 2012 to generate a catalog of seismicity that enhances the US Geological Survey Hawaiian Volcano Observatory’s public seismic catalog with four times the number of earthquakes, which were then grouped by waveform similarity. Analysis of subtle delays in the timing of arrivals of scattered waves between pairs of earthquakes in this catalog yields a history of small changes in the shallow seismic velocity structure of the volcano. Seismic velocities have been shown at other volcanoes to change during unrest and eruption. My results show a decrease in seismic velocity centered on the summit beginning in September 2022, corresponding to the onset of a vigorous precursory swarm of seismic activity and shallow inflation. During the eruption itself, I observe large changes due likely to dike opening along the northeast rift zone and deflation of the summit reservoir. However, seismic velocity changes associated with non-volcanic sources such as ground shaking from large earthquakes and meteorological influences at seasonal and diurnal time scales are also observed, and these dominate the velocity changes prior to the eruption. Proper accounting of these effects will be a requirement for use in real-time monitoring, and this work serves as a starting point in that endeavor for Mauna Loa.

Exploring the Gravitational Impact of Spring Tides on Earthquake Occurrence in Central Asia

Abstract In this study, we analyzed variations in seismic activity in relation to the syzygy of the Sun, Moon, and Earth, corresponding to the semilunar or spring tide (ST) cycle. We focused on two regions within the Central Asian Orogenic Belt: the North Tien Shan region (NTSR) and the Xinjiang Uygur tectonic province (XUTP), using data from two regional earthquake catalogs. For the NTSR, we used data from the Seismological Experimental and Methodical Expedition in Almaty, Kazakhstan, covering 1 January 1970 to 31 July 2022. For the XUTP, we utilized data from the Data Management Centre of China National Seismic Network at the Institute of Geophysics China Earthquake Administration, spanning 1 January 1987 to 30 June 2019. We applied the nonparametric Kuiper’s test to evaluate cyclic variations in seismic activity, focusing on different magnitude thresholds in the regional catalog. For both regions, the original catalog as well as the catalogs of mainshocks and aftershocks were analyzed. The nearest-neighbor declustering method was employed to decompose the seismic catalogs. Kuiper’s test statistics revealed varying trends in the occurrence times for mainshocks and aftershocks during ST cycles. The test showed “the good-to-fit” to a uniform distribution for the timing of mainshocks within the ST cycle. Kuiper’s test revealed notably consistent trends in aftershock occurrence times during ST cycles with a significance level α of 0.001 or less in both regions analyzed. These results suggest a potential correlation between aftershock activity and the gravitational effects of the ST cycle. The most substantial decrease in aftershocks occurs during the micromoon ST cycle, when the distance between the Moon and the Earth exceeds 400,000 km.

Seismicity and Present‐Day Crustal Deformation in the Southern Puna Plateau

AbstractThe Southern Puna plateau in the central Andes has a complicated tectonic history that includes episodes of distributed shortening and extension, lithospheric delamination, uplift and Quaternary backarc volcanism. In this study, the upper crustal structure and present‐day deformation in this area is investigated using a new regional earthquake catalog derived with a deep‐learning‐based phase picker. Results show abundant strike‐slip seismicity at shallow depths in the eastern Southern Puna plateau that reveals active fault systems in the area and indicates N‐S extension/E‐W compression that changes orientation and relative magnitude from north to south. A broad zone of seismic quiescence in the western plateau may indicate a zone of upper crustal decoupling from large‐scale deformation. The region separating the western and eastern plateau exhibits a complex stress field that can be related to the boundary of east/west oriented middle‐to‐lower crustal flow in the main volcanic arc. Southeast of the plateau in the Sierras Pampeanas, crustal seismicity deepens and is dominated by conjugate reverse faulting structures associated with the direction of plate convergence. Vp and Vs seismic velocity models of the upper crust obtained through local earthquake tomography with the improved seismic catalog show low‐velocity anomalies near intermontane basins, except in the Antofagasta basin where a high‐velocity anomaly possibly represents shallow intrusive component of Quaternary basaltic volcanism. Below the Cerro Galan caldera, an upper crustal 10‐day long earthquake swarm is observed which may indicate local stress perturbations from fluids at the top of the crustal magmatic system that feeds this volcano.

Probabilistic Seismic Hazard Assessment for the North China Plain Earthquake Belt: Sensitivity of Seismic Source Models and Ground Motion Prediction Equations

In this study, a multi-source data fusion method was proposed for the development of a Hybrid seismic hazard model (HSHM) in China by using publicly available data of the 5th Seismic Ground Motion Parameter Zoning Map (NSGM) and historical seismic catalogues and integrating with modern ground motion prediction equations (GMPEs). This model incorporates the characteristics of smoothed seismicity and areal sources for regional seismic hazard assessment. The probabilistic seismic hazard for the North China Plain earthquake belt was investigated through sensitivity analysis related to the seismicity model and GMPEs. The analysis results indicate that the Hybrid model can produce a consistent result with the NSGM model in many cases. However, the NSGM model tends to overestimate hazard values in locations where no major events have occurred and underestimate hazard values in locations where major events have occurred. The Hybrid model can mitigate the degree of such biases. Compared to the modern GMPEs, the GMPE with epicentral distance measures significantly underestimate the seismic hazard under near-field and large-magnitude scenarios. In addition, a comparison of the uniform hazard spectra (UHS) obtained by the models, with China’s design spectrum, shows that the current design spectrum is more conservative than the calculated UHS.

The catalogue of 1987–2023 earthquakes in the western (French) Alps north of 43.5°N

The western Alpine belt is the focal point of moderate but constant seismic activity. Numerous geodynamic and seismological studies underline the many scientific questions linked to Alpine dynamic processes and the associated rate of deformation resulting from late continental collision phases, and the assessment of seismic hazard and its associated risks. Seismic monitoring in the northern French Alps was launched in 1987 with the installation of the first seismological network, hosted by the Earth Science Observatory at Grenoble Alpes University. Since then, the network has been updated over the years, keeping pace with technical developments in seismic instrumentation, monitoring and data management. The homogenization of information described in this scientific data paper is motivated by the need for a single instrumental seismic catalogue, consisting in arrival times picks and hypocenters using the same velocity model. The final catalogue SISMalp-1987–2023 (CC-BY-4.0) contains 50,822 weak-to-moderate earthquakes with local magnitude between −1.8 and 5.0; it can serve as a basis for all general seismotectonic and risk assessment studies in the region.

Seismic Response of Hectometer‐Scale Fracture Systems to Hydraulic Stimulation in the Bedretto Underground Laboratory, Switzerland

AbstractWe performed a series of hydraulic stimulations at 1.1 km depth in the Bedretto underground laboratory, Switzerland, as part of an overall research strategy attempting to understand induced seismicity on different scales. Using an ultra‐high frequency seismic network we detect seismic events as small as Mw < −4, revealing intricate details of a complex fracture network extending over 100 m from the injection sites. Here, we outline the experimental approach and present seismic catalogs as well as a comparative analysis of event number per injection, magnitudes, b‐values, seismogenic index and reactivation pressures. In our first‐order seismicity analysis, we could make the following observations: The rock volume impacted by the stimulations in different intervals differs significantly with a lateral extent from a few meters to more than 150 m. In most intervals multiple fractures were reactivated. The seismicity typically propagates upwards toward shallower depth on parallel oriented planes that are consistent with the stress field and seem to a large extent associated with preexisting open fractures. This experiment confirms the diversity in seismic behavior independent from the injection protocol. The overall seismicity patterns demonstrate that multi‐stage stimulations using zonal isolation allow developing an extended fracture network in a 3D rock volume, which is necessary for enhanced geothermal systems. Our stimulations covering two orders of magnitude in terms of injected volume will give insights into upscaling of induced seismicity from underground laboratory scale to field scale.

Improvement of the 2007–2015 Earthquake Catalog Along the 300 km Long Postglacial Merasjärvi–Stuoragurra Fault Complex in Northern Fennoscandia Using Automatic Event Detection

We present an updated and validated seismic catalog for the northern Fennoscandian region, focusing on postglacial faults from the Merasjärvi fault system in the southwest to the Iešjávri fault system in the northeast. This work involved a comprehensive review of continuous waveforms derived from open datasets from 2007 to 2015 and processed using the Regressive ESTimator algorithm. The primary objective was to refine the delineation of seismicity along the above-mentioned postglacial faults and highlight their seismic potential. Our analysis revealed distinct waveform patterns originating primarily from two main sources: approximately 15% were associated with areas mapped as postglacial faults, and the remainder of the events outside these areas, 89%, were concentrated in areas with active mines. Compared to previously reported events in the Fennoscandian Earthquake Catalogue (FENCAT), we observed a 22% increase in seismic activity within postglacial fault zones. These results demonstrate that the Regressive ESTimator algorithm not only improves the detection of tectonic seismicity but also effectively identifies seismic signals resulting from mining activities in the study area. The Merasjärvi, Lainio–Suijavaara, Palojärvi, and Maze and Iešjávri fault systems appear to form a continuous deformation complex of approximately 300 km long, which we propose naming the Merasjärvi–Stuoragurra fault complex.

Monthly Maximum Magnitude Prediction in the North–South Seismic Belt of China Based on Deep Learning

The North–South Seismic Belt is one of the major regions in China where strong earthquakes frequently occur. Predicting the monthly maximum magnitude is of significant importance for proactive seismic hazard defense. This paper uses seismic catalog data from the North–South Seismic Belt since 1970 to calculate and extract multiple seismic parameters. The monthly maximum magnitude is processed using Variational Mode Decomposition (VMD) with sample segmentation to avoid information leakage. The decomposed multiple modal data and seismic parameters together form a new dataset. Based on these datasets, this paper employs four deep learning models and four time windows to predict the monthly maximum magnitude, using prediction accuracy (PA), False Alarm Rate (FAR), and Missed Alarm Rate (MR) as evaluation metrics. It is found that a time window of 12 generally yields better prediction results, with the PA for Ms 5.0–6.0 earthquakes reaching 77.27% and for earthquakes above Ms 6.0 reaching 12.5%. Compared to data not decomposed using VMD, traditional error metrics show only a slight improvement, but the model can better predict short-term trends in magnitude changes.

Characteristics of the seismogenic zone in an arc-continent collision belt: insights from seismic b values in Eastern Taiwan

Eastern Taiwan overlies a suture zone between the Eurasian Plate and the Philippine Sea Plate and is characterized by frequent earthquakes, often resulting in significant disasters. Notably, the region exhibits characteristics such as a high frequency of earthquakes and a short recurrence period for intense seismic events. While prior research has explored seismic b values across various periods in Taiwan, detailed investigations of the b value in the eastern region are lacking. This study employs the earthquake catalog compiled by the Taiwan Central Weather Administration to analyze spatial–temporal variations in b values in eastern Taiwan. The analysis encompasses seismic events occurring between January 1996 and June 2019. The seismic catalog is divided into three distinct time periods related to large seismic events: period I, 1996–2003 (the Chengkung earthquake); period II, 2003–2013 (the Ruisui earthquake); and period III, 2013–2019 (the Hualien earthquake). Our results indicate that most seismic events with a magnitude greater than 6 are associated with low b values. The overall b value increases during period II and then decreases substantially during period III. Although the estimated b values changed slightly, but the uncertainty in b values remained stable in this study. The epicenters of large earthquakes often overlap with areas with lower b values, especially in plate suture zones, which means that areas with lower b values usually have a higher probability of larger earthquakes. Given the extremely high potential for a catastrophic earthquake, mitigating measures should be adopted at all times.Graphical

Deformed archaeological remains at Lilybaeum in Western Sicily (southern Italy): possible ground signatures of a missed large earthquake

AbstractArchaeoseismic analysis performed in Western Sicily points to deformed archeological remains at Lilybaeum, a Punic coastal city founded in 397 B.C. at the Island’s westernmost edge. Starting from the direct observation of deformed ruins, an interdisciplinary work strategy, which included field-structural analysis, drone-shot high-resolution aerial photogrammetry, and geophysical prospecting, was employed to investigate whether the identified deformations may represent the ground effects of a previously unknown large earthquake in the area. Among the unearthed remains, some mosaics and a stone-paved monumental avenue show evidence of tectonic deformation, being fractured, folded, and uplifted. The trend of folding and fracturing is consistent with the NNW–SSE oriented tectonic max stress axis to which Western Sicily is currently subjected. Displacement along a fracture deforming the Decumanus Maximus together with the finding of a domino-type directional collapse, enable us to interpret the observed deformation as the ground signature of a coseismic slip. The seismic rupture occurred along a previously unmapped deformation front that fits well within the seismotectonic context of Western Sicily. Measured offset, geophysical prospecting, and age-constraints all suggest the possibility that a highly-energetic earthquake nucleated in the area following a coseismic rupture along a NE–SW trending back-verging reverse fault towards the end of the fourth century CE. Since seismic catalogs do not provide evidence of such a large earthquake, this event might represent a missed entry in the historical seismic record. This finding provides constraints to redefine the seismic hazard of Western Sicily, a region where recurrence-time intervals for large earthquakes are still unknown.

Machine-learning based location of the 2021 MW 7.4 Maduo, Qinghai, China earthquake sequence: Insight into intraplate seismogenesis

On 22 May 2021, an MW 7.4 earthquake occurred in Maduo County, Qinghai Province, China, which is located on the Kunlun Mountain Pass-Jiangcuo fault inside the Bayan Har block, providing a good opportunity to investigate seismogenesis of large intraplate earthquakes. We analyze two years of continuous seismic data from June 2021 to June 2023, which were recorded at 34 portable seismic stations of the MaduoArray deployed in the source zone by our group. The LOC-FLOW workflow of automatic detection and location is applied to construct a complete and high-precision seismic catalog for the region, which includes machine-learning phase picking (PhaseNet), earthquake phase association (REAL), velocity model updating and station correction (VELEST), absolute earthquake location (HypoInverse), and relative location (HypoDD). As a result, 78,832 Maduo aftershocks and other local earthquakes are detected and relocated precisely. Our results show that the length of the Maduo aftershock zone is ∼170 km, which is mainly distributed along the NWW-SEE oriented Kunlun Mountain Pass-Jiangcuo fault, and there is a horsetail bifurcation feature at the eastern end of the aftershock sequence. The seismogenic fault is nearly vertical, and local seismicity occurs on both sides of the fault. Our results also show that there is no seismic gap or aftershock sparse area in the region. Previous studies have revealed a low-velocity and high-conductivity anomaly below the source zone, reflecting fluids ascending from the lower crustal flow. These results provide new insights into the cause of the 2021 Maduo earthquake.

A Unified Seismicity Catalog Development for Saudi Arabia: Multi-Network Fusion and Machine Learning-Based Anomaly Detection

A Unified Seismicity Catalog Development for Saudi Arabia: Multi-Network Fusion and Machine Learning-Based Anomaly Detection

Interplay of tectonic and dynamic processes shaping multilayer extensional system in southern-central Apennines

High-hazard seismic zones can remain silent over centuries with meager seismicity rates challenging our understanding of seismic processes. We focus on the comprehensive analysis of cascading episodes of swarms and seismic sequences following the 2009 L’Aquila mainshock (MW 6.3) in the southern-central Apennine that previously experienced ~ M7 earthquakes. We enhance the seismic catalog, unmasking low-magnitude seismicity down to completeness magnitude ML ~ 0, and we unveil that the microseismicity might be secondarily triggered by the L’Aquila mainshock, influencing the frictional properties in the nearby fault zones or opening fault valves generating the intense seismic activity detected from 2009 to 2013. The diffusivity, observed in the most seismic episodes, and the high Vp/Vs values (> 1.88) indicate fluid circulation promoting multilayered extensional seismicity within 11–15 km and 16–23 km depth ranges. Mapping the 3D distribution of seismicity alongside geological data reveals an evident tectonic influence, unveiling unknown geometric aspects and providing the first evidence of a NNE-dipping deformation zone bounding at depths of 11–15 km the overlying fault system. Deeper seismicity suggests a mantellic CO2 ascending shape. These findings enrich the literature on tectonic seismic swarms in extensional domains, providing essential constraints on fluid involvement in the seismic processes and contributing to forthcoming discussions on the seismotectonic setting in high-seismic-risk areas of the Apennines of Italy.

Characterization and Evolution of Seismic Sequences in the Normal Fault Environment of the Southern Apennines

AbstractThe use of seismic catalogs enhanced through advanced detection techniques improves the understanding of earthquake processes by illuminating the geometry and mechanics of fault systems. In this study, we performed accurate hypocentral locations, source parameters estimation and stress release modeling from catalogs of microseismic sequences nucleating in the complex normal fault system of the Southern Apennines (Italy). The application of advanced location techniques resulted in the relocation of ∼30% of the earthquakes in the enhanced catalogs, with hypocenters clearly identifying local patches on kilometer‐scale structures that feature consistent orientation with the main faults of the area. When mapping the stress change on the fault plane, the inter‐event distance compared to the size of the events suggests that the dominant triggering mechanism within the sequences is static stress transfer. The distribution of events is not isotropic but dominantly aligned along the dip direction. These slip‐dominated lineations could be associated with striations related to fault roughness and could map the boundary between locked and creeping domains in Apulian platform and basement.

(Re)Discovering the Seismicity of Antarctica: A New Seismic Catalog for the Southernmost Continent

Abstract We apply a machine learning (ML) earthquake detection technique on over 21 yr of seismic data from on-continent temporary and long-term networks to obtain the most complete catalog of seismicity in Antarctica to date. The new catalog contains 60,006 seismic events within the Antarctic continent for 1 January 2000–1 January 2021, with estimated moment magnitudes (Mw) between −1.0 and 4.5. Most detected seismicity occurs near Ross Island, large ice shelves, ice streams, ice-covered volcanoes, or in distinct and isolated areas within the continental interior. The event locations and waveform characteristics indicate volcanic, tectonic, and cryospheric sources. The catalog shows that Antarctica is more seismically active than prior catalogs would indicate, examples include new tectonic events in East Antarctica, seismic events near and around the vicinity of David Glacier, and many thousands of events in the Mount Erebus region. This catalog provides a resource for more specific studies using other detection and analysis methods such as template matching or transfer learning to further discriminate source types and investigate diverse seismogenic processes across the continent.

Characterization of Micro-seismic Activity in Northern Cyprus Using Complexity and Corner Frequency Methods

Cyprus is an island country located in the eastern Mediterranean, to the south of Türkiye and the western of Syria and Lebanon, and is a popular tourist destination. Due to being surrounded by seas on all four sides, meticulous planning of rescue, assistance, and evacuation plans is necessary in the face of disasters such as earthquakes and tsunamis. Tectonically, the southern part of the island is controlled by the Cyprus Arc, while the northern part is dominated by the Kyrenia Range. The demand for raw materials for construction and industry is met through controlled quarry blasting operations carried out by open-pit quarry companies in the districts of Kyrenia and Nicosia. As a result, both natural and artificial seismic events occur in the region, and these quakes are documented in seismic catalogs by seismology centers. However, due to the low energy content of micro-seismic events and the inadequacy of seismic stations on the island, the source types of these seismic events can be misidentified in the catalogs. In this context, the study focuses on 122 seismic events with magnitudes between 0.9≤Ml≤2.7 that occurred in Northern Cyprus during the January 2018 - December 2021 period (4 years). The seismic events recorded by the station LFK, operated by Boğaziçi University Kandilli Observatory and Earthquake Research Institute Regional Earthquake-Tsunami Monitoring Center (KOERI-RETMC), were classified using Linear and Quadratic Discriminant Functions based on complexity and corner frequency methods. According to the results obtained, 10 of the 122 seismic events were identified as natural, and 96 were determined to be artificial, resulting in a general success rate of 86.89%. However, classification results for 16 seismic events were inconclusive with the methods used. As a result, more detailed secondary analyses should be conducted to accurately determine the source types of micro-seismic events, and the seismic catalogs should be updated accordingly.