Earthquake Catalog Research Articles

Source mechanism of the 2023 Ms 5.5 earthquake in Subei, Gansu Province revealed by relocated aftershocks and InSAR: complement to the ‘shallow slip deficit’ of the eastern boundary of the Altyn Tagh fault

The Ms 5.5 earthquake struck on 24 October 2023, in Subei County, Gansu Province, China, occurring along the eastern segment of the Altyn Tagh fault. It raises the question of whether this earthquake is linked to the ongoing shortening slip rate along this segment or triggered by other seismic events. Analyzing the fault geometry of the Subei earthquake and understanding the significance of the weakening activity rate for seismic hazards in neighboring regions is crucial. The surface deformation from small- and medium-sized earthquakes (magnitudes less than Mw5.5) is often subtle, and the coseismic deformation detected by interferometric synthetic aperture radar (InSAR) is vulnerable to atmospheric disturbances, leading to significant measurement errors. Moreover, inaccuracies in the regional crustal velocity structure can cause errors in earthquake localization based on seismic data. These challenges complicate the establishment of a rupture model for seismogenic faults and hinder the inversion of fault slip models. To overcome these limitations, we employed the time-series InSAR stacking method and aftershock relocation to determine the fault geometry of the Subei earthquake. A two-step inversion method was utilized to ascertain both the fault geometry and slip distribution. Our modeling indicates that the 2023 Subei earthquake had a thrust mechanism with a component of strike-slip. The rupture did not reach the surface, with the maximum fault slip measuring 0.45 m at a depth of 2.5–3.5 km. The fault dips westward, and the moment magnitude is calculated at 5.4. This earthquake is associated with the ongoing weakening of the left-lateral strike-slip rupture along the Altyn Tagh fault in the Subei region. Furthermore, retrograde thrust tectonics significantly contribute to the absorption of accumulated stress during this process.Our findings highlight the potential of utilizing time-series InSAR images to enhance earthquake catalogs with geodetic observations, offering valuable data for further studies of the earthquake cycle and active tectonics. This approach is also applicable in other tectonically active regions, enhancing understanding of seismic hazards and risk assessment.

Geophysical Coupling Before Three Earthquake Doublets Around the Arabian Plate

In this study, we analysed lithospheric, atmospheric, and top-side ionospheric magnetic field data six months before the three earthquake doublets occurred in the last ten years around the Arabian tectonic plate. They occurred in 2014, close to Dehloran (Iran), in 2018, offshore Kilmia (Yemen) and in 2022, close to Bandar-e Lengeh (Iran). For all the cases, we considered the equivalent event in terms of total released energy and mean epicentral coordinates. The lithosphere was investigated by calculating the cumulative Benioff strain with the USGS earthquake catalogue. Several atmospheric parameters (aerosol, SO2, CO, surface air temperature, surface latent heat flux humidity, and dimethyl sulphide) have been monitored using the homogeneous data from the MERRA-2 climatological archive. We used the three-satellite Swarm constellation for magnetic data, analysing the residuals after removing a geomagnetic model. The analysis of the three geo-layers depicted an interesting chain of lithosphere, atmosphere, and ionosphere anomalies, suggesting a geophysical coupling before the Dehloran (Iran) 2014 earthquake. In addition, we identified interesting seismic accelerations that preceded the last 20 days, the Kilmia (Yemen) 2018 and Bandar-e Lengeh (Iran) 2022 earthquake doublets. Other possible interactions between the geolayers have been observed, and this underlines the importance of a multiparametric approach to properly understand a geophysical complex topic as the preparation phase of an earthquake.

Seismotectonic Map of The Sinai Triple Junction

The geodynamic evolution of the Sinai Triple Junction, a highly deformed and seismically active area, is controlled by the Red Sea rift, Gulf of Suez and Aqaba-Dead Sea conjunctions. However, the driving forces for the focusing deformation at crustal depths beneath this area are still ambiguous. Here, we provide an updated seismotectonic map of the area relying on updated seismological and geodetic datasets. A homogenized earthquake catalog has been compiled from well-located earthquakes (> Mw 2.0) by the Egyptian Seismic Network and International Seismological Center in the period between 1990 and 2020.We calculated the average b-value along three seismogenic zones including Gulf of Aqaba, northern Red Sea and Gulf of Suez that amount to 1.1, 0.99 and 0.97, respectively. Additionally, we complied and updated a comprehensive P-wave-based database for the fault plane solutions in the area for events with Mw> 3.5 till 2023. Furthermore, a unified velocity field for the region as well as slip-rate and locking-depth at the active fault segments were estimated from a consistent geodetic dataset from peer-reviewed GPS velocities between 1999 and 2018. Results indicate a dominant NNE left-lateral strike-slip fault with normal component along the Gulf of Aqaba. Pure NW-SE to WNW-ESE dip-slip normal faulting, associated with a strike-slip component in some cases, is dominating the northern and central parts of the Gulf of Suez, whereas pure normal dip-slip movement with an NNE–SSW extension in a horizontal direction is observed in the southern part of the gulf. The estimated slip-rate and locking-depths at the Aqaba fault segments falls between 4.8-4.9 mm/yr and 8-12 km, respectively.

A Proposal for a High-Frequency Earthquake Magnitude (m3Hz) for Seismic Hazard and Rapid Damage Assessment

Abstract Obtaining a size estimate of an earthquake that well represents its potential impact is of primary importance in seismology and earthquake engineering. The magnitude scales currently used in earthquake catalogs, mainly Ms and Mw, are not designed to represent the variability of shaking due to variations in high-frequency radiation. Conversely, ML, which best describes the seismic-wave energy released by an earthquake, is unable to provide an accurate representation of events with magnitudes exceeding approximately 6.5. Therefore, building on and extending the concept of high-frequency magnitude (Atkinson and Hanks, 1995), we propose in this article a new high-frequency magnitude scale m3Hz. This magnitude is estimated by also considering correction factors for the crustal model and site effects. We apply the new magnitude scale m3Hz to two data sets (Central Italy and Japan) and validate its ability to capture high-frequency radiation effects through comparison with available source parameters. Finally, a procedure for estimating m3Hz is also proposed for preinstrumental earthquakes. This procedure addresses a limitation of previous magnitude estimation techniques, such as Me, and is validated for four events for which both macroseismic and instrumental intensity data are available.

Relocation of the 8 September 2023 High Atlas, Morocco, Earthquake Aftershock Sequence

ABSTRACT The earthquake that occurred on 8 September 2023, with a magnitude of 6.8, was the most destructive earthquake event in Morocco in the past decade. This earthquake took place in the Al Haouz region, located in the western part of the High Atlas Mountain range. To better understand what caused and triggered this earthquake, the earthquake catalogs including P and S arrival times were collected from the Moroccan seismic network and combined with regional data from the International Seismological Centre. The mainshock and aftershocks were relocated by using iLoc, a state-of-the-art single-event location algorithm, and then by the multiple event location double-difference algorithm, hypoDD. The improved earthquake relocations using iLoc and the double-difference methods provide sharper lineation of seismicity and agree well with tomographic images of the earthquake zone. The seismicity distribution and the focal mechanism of the mainshock indicate that the earthquake sequence has occurred along the South Atlas fault system.

Classifying Seismic Events Linked to Solar Activity: A Retrospective LSTM Approach Using Proton Density

The influence of solar activity on seismic activity is a subject of debate. Previous studies have shown that there is sometimes a correlation and sometimes a contradiction between solar activity maxima and large earthquakes. Long-term memory neural network is used to study the relationship between solar activity and seismic activity. This study emphasizes retrospective classification rather than direct prediction, refining the LSTM architecture to maximize classification accuracy and processing data from the Solar and Heliospheric Observatory and the U.S. Geological Survey earthquake catalogs. A declustering technique is used to select large seismic events and weighted learning corrects for class imbalances. The LSTM model accurately classified earthquakes (84.47%) and proton density variations. The results support the theory that solar activity, in particular proton density, can anticipate earthquake events.

The temporal and spatial evolution characteristics of induced seismicity in the Changning shale gas field based on dense array

In this paper, we performed microseismicity detection and location using the deep learning method and obtained a high-precision earthquake catalog in the Changning gas field, which is one of the largest shale gas production demonstration areas in China. We found that the spatial and temporal characteristics of seismicity in the region are indicative of its correlation with industrial operations. The distribution of earthquakes at depth reflected variations in reservoir depth and provided valuable constraints on it. The horizontal layered distribution of earthquakes at depth is due to the formation of fracture surfaces from interconnected fractures near the reservoir during hydraulic fracturing (HF) operations, which clearly demonstrates how HF operations impact seismicity. We suggested that the horizontally layered distribution was driven by two fundamental mechanisms: reactivation of pre-existing faults during HF and injection of high-pressure fluids into the reservoir, leading to fracture creation. Several ML ≥4 earthquakes, which did not occur on well-defined seismogenic faults, may have been triggered by pore elastic coupling resulting from regional stress accumulation and fluid injection. Significantly, the ML 4.9 seismic sequence occurring at the basement indicates that fracking has reactivated and promoted pre-existing faults, highlighting the need for further investigation into potential seismic hazards in the region.

Forecasting future earthquakes with deep neural networks: application to California

SUMMARY We use the spatial map of the logarithm of past estimated released earthquake energies as input of fully convolutional networks (FCN) to forecast future earthquakes. This model is applied to California and compared with an elaborated version of the epidemic type aftershock sequence (ETAS) model. Our long-term earthquake forecast simulations show that the FCN model is close to the ETAS model in forecasting earthquakes with $M \ge 3.0,\,\,4.0,\,\,{\rm{and\,\,}}5.0$ according to the Molchan diagram. Moreover, training and implementing the FCN model is 2000–4000 times faster than calibrating the ETAS model and generating its probabilistic forecasts. The FCN model is straightforward in terms of its neural network structure and feature engineering. It does not require extensive knowledge of statistical seismology or the analysis of earthquake catalogue completeness. Using the earthquake catalogue with $M \ge 0$ as FCN input can enhance the model's performance in some time–magnitude forecasting windows.

Assessing the Adequacy of Earthquake Catalog Sampling for Long-Term Seismicity in Low-to-Moderate Seismic Regions: A Geodetic Perspective

Abstract Seismic hazard assessment in low-to-moderate seismicity regions can benefit from the knowledge of surface deformation rates to better constrain earthquake recurrence models. This, however, amounts to assuming that the known seismicity rate, generally observed over historical times (i.e., up to a few centuries in Europe), provides a representative sample of the underlying long-term activity. We here investigate how this limited sampling can affect the estimated seismic hazard and whether it can explain the disagreement between the seismic moment loading rate as seen by nowadays Global Navigation Satellite Systems (GNSS) measurements and the seismic moment release rate by past earthquakes, as is sometimes observed in regions with limited activity. We approach this issue by running simulations of earthquake time series over very long timescales that account for temporal clustering and the known magnitude–frequency distribution in such regions, and that those are constrained to a seismic moment rate balance between geodetic and seismicity estimates at very long timescales. We show that, in the example of southeastern Switzerland, taken here as a case study, this sampling issue can indeed explain this disagreement, although it is likely that other phenomena, including aseismic deformation and changes in strain rate due to erosional and/or glacial rebound, may also play a significant role in this mismatch.

Improvements and Heterogeneities of the Global Centroid Moment Tensor Catalog

Abstract Earthquake catalogs are heterogeneous, especially those developed over long time spans. Changes in seismological monitoring, which provides the records on which these catalogs are based, are common. Typically, instruments and networks become more sensitive over time, allowing for the detection and characterization of smaller earthquakes. In pursuit of improvement, new methods for routine data analysis are occasionally introduced, modifying the procedures for catalog compilation. The resulting heterogeneities may not be evident to users, but they should be unveiled and considered in any application of the catalog, especially in statistical seismology, which analyzes large earthquake data sets. The Global Centroid Moment Tensor catalog is considered the most homogeneous database of global seismicity. However, a detailed analysis of its heterogeneities has been lacking. This work reviews changes in the catalog’s development from 1976 to 2023 and reveals how these have caused improvements and heterogeneities in the resulting data. Several periods are distinguished, separated by milestones in the methods employed for moment tensor inversion and catalog compilation, as well as by the advent of global broadband monitoring in 2004. These changes are shown to have caused variations in the catalog’s completeness and in the determinations of centroid depths, scalar seismic moments, and moment tensors. The magnitude of completeness is measured here in detail, both temporally and spatially. It has decreased over the years and shows spatial variations within each period, correlated to regional differences in network monitoring and compilation biases. Moment tensor determinations have been significantly different since 2004, resulting in a different frequency distribution of rake angles and a different dependence of the double-couple component as a function of rake. This work is expected to benefit all future uses of the catalog, enabling better characterization of seismicity properties and improved building and testing of models for earthquake occurrence.

Algorithmic Identification of the Precursory Scale Increase Phenomenon in Earthquake Catalogs

Abstract The precursory scale increase (Ψ) phenomenon describes the sudden increase in rate and magnitude in a precursory area AP, at precursor time TP, and with precursor magnitude MP prior to the upcoming large earthquake with magnitude Mm. Scaling relations between the Ψ variables form the basis of the “Every Earthquake a Precursor According to Scale” (EEPAS) earthquake forecasting model. EEPAS is a well-established space–time point process model that forecasts large earthquakes in the medium term, that is, the coming months to decades, depending on Mm. In Aotearoa New Zealand, EEPAS contributes to hybrid models for public earthquake forecasting and to the source model of time-varying seismic hazard models, including the latest revision of the National Seismic Hazard Model. The Ψ phenomenon was recently shown not to be unique for a given earthquake, with smaller precursory areas AP associated with larger precursor times TP and vice versa. This trade-off between AP and TP has also been found for the spatial and temporal distributions of the EEPAS models. Detailed analysis of the Ψ phenomenon has so far been limited by the manual and labor-intensive procedure of identifying Ψ in earthquake catalogs. Here, we introduce two algorithms to automatically detect Ψ and apply them to real and simulated earthquake catalog data. By randomizing the catalog and removing aftershocks, we confirm that the Ψ phenomenon is a feature of space–time earthquake clustering prior to major earthquakes. Multiple Ψ identifications confirm the trade-off between AP and TP, and the scaling relations for both real and simulated catalogs are consistent with the original scaling relations on which EEPAS is based. We identify opportunities for future work to refine the algorithms and apply them to physics-based simulated catalogs to enhance the understanding of Ψ. A better understanding of Ψ has the potential to improve forecasting of large upcoming earthquakes.

Microseismicity at the Time of a Large Creep Event on the Calaveras Fault is Unresponsive to Stress Changes

The potential relationship between surface creep and deeper geological processes is unclear, even on one of the world’s best-studied faults. From June to August 2021, a large creep event with surface slip of more than 16 mm was recorded on the Calaveras fault in California, part of the San Andreas fault system. This event initially appeared to be accompanied by along-fault migration of seismicity, suggesting it penetrated to depth. Other studies have suggested that surface creep events are likely a shallow feature, unrelated to deep seismicity. To provide more detail on the relationship between earthquakes, surface creep, and potential aseismic slip at seismogenic depth, we tripled the number of earthquakes in the Northern California Earthquake Catalog in the region of the creep event for all of 2021. This was accomplished by implementing earthquake detection techniques based on both template matching (EQCorrscan) and AI-based automatic earthquake phase picking (PhaseNet). After manual inspection, the detected earthquakes were first located using Hypoinverse and subsequently relocated via GrowClust. Our enhanced catalog indicates that the spatiotemporal pattern of earthquakes here is not strongly influenced by the creep event and is better explained by structural heterogeneity than transient stress changes, indicating a decoupling of seismicity rate and surficial creep on this major fault.

Onset of Aftershocks: Constraints on the Rate-and-State Model

Abstract Aftershock rates typically decay with time t after the mainshock according to the Omori–Utsu law, R(t)=K(c+t)−p, with parameters K, c, and p. The rate-and-state (RS) model, which is currently the most popular physics-based seismicity model, also predicts an Omori–Utsu decay with p = 1 and a c-value that depends on the size of the coseismic stress change. Because the mainshock-induced stresses strongly vary in space, the c-value should vary accordingly. Short-time aftershock incompleteness (STAI) in earthquake catalogs has prevented a detailed test of this prediction so far, but the newly developed a-positive method for reconstructing the true earthquake rate now allows its testing. Using previously published slip models, we calculate the coseismic stress changes for the six largest mainshocks in Southern California in recent decades and estimate the maximum shear as a scalar proxy of the coseismic stress tensor. Aftershock rates reconstructed for events in different stress ranges show that the rates follow a power law with p = 1 independent of stress with no clear sign of a c-value. The onset of the power-law decay is abrupt and more delayed in areas with smaller stress changes. The observations do not necessarily contradict the RS model, as STAI limits the resolution for early aftershocks, and the RS model can reproduce the observations for specific Aσ values. However, the observations lead to strong constraints, namely Aσ<10 kPa and a power-law decay of the background rate with distance to the fault, with exponent 2.7.

New Features in the pyCSEP Toolkit for Earthquake Forecast Development and Evaluation

Abstract The Collaboratory for the Study of Earthquake Predictability (CSEP) is a global community dedicated to advancing earthquake predictability research by rigorously testing probabilistic earthquake forecast models and prediction algorithms. At the heart of this mission is the recent introduction of pyCSEP, an open-source software tool designed to evaluate earthquake forecasts. pyCSEP integrates modules to access earthquake catalogs, visualize forecast models, and perform statistical tests. Contributions from the CSEP community have reinforced the role of pyCSEP in offering a comprehensive suite of tools to test earthquake forecast models. This article builds on Savran, Bayona, et al. (2022), in which pyCSEP was originally introduced, by describing new tests and recent updates that have significantly enhanced the functionality and user experience of pyCSEP. It showcases the integration of new features, including access to authoritative earthquake catalogs from Italy (Bolletino Seismico Italiano), New Zealand (GeoNet), and the world (Global Centroid Moment Tensor), the creation of multiresolution spatial forecast grids, the adoption of non-Poissonian testing methods, applying a global seismicity model to specific regions for benchmarking regional models and evaluating alarm-based models. We highlight the application of these recent advances in regional studies, specifically through the New Zealand case study, which showcases the ability of pyCSEP to evaluate detailed, region-specific seismic forecasts using statistical functions. The enhancements in pyCSEP also facilitate the standardization of how the CSEP forecast experiments are conducted, improving the reliability, and comparability of the earthquake forecasting models. As such, pyCSEP exemplifies collaborative research and innovation in earthquake predictability, supporting transparent scientific practices, and community-driven development approaches.

Active Fault Interaction in the Eastern Betic Cordillera: A Model of Coseismic and Postseismic Stress Transfer Following Historical Earthquakes in SE Spain

AbstractThe Spanish historical earthquake catalog reveals a significant number of destructive earthquakes (Mw 6.0) that have occurred along the Eastern Betic Shear Zone (EBSZ), a major active fault system at the Eastern Betic Cordillera. However, during the last century, seismic activity in this region has been characterized by the absence of large‐magnitude events, complicating yet underlining the importance of understanding active fault interaction in the EBSZ. In this work, we focus on the northeastern half of the EBSZ and apply a similar method to Yazdi et al. (2023, https://doi.org/10.1029/2023tc007917). We model a series of consecutive earthquakes ( VII; Mw 5.0) between the and centuries occurring along the Alhama de Murcia, Carrascoy and Bajo Segura faults, and assess their coseismic impact as well as explore the evolution of Coulomb stress changes due to the postseismic viscoelastic relaxation on the surrounding faults. Our results support a stress‐triggering connection between subsequent events and shed light on identifying the responsible faults for some. Our findings indicate that 1673 Mw 6.0 Orihuela and 1674 Mw 6.2 Lorca earthquakes, seemingly unrelated, could have both contributed to the forthcoming Mw 5.0 earthquake clustering along the Carrascoy fault. We show that attributing the 1829 Mw 6.6 Torrevieja earthquake to the Bajo Segura fault better describes the occurrence of earthquake sequences between the mid‐ and early centuries. Finally, we propose the Alhama de Murcia and the Carrascoy faults as responsible initiators for the early ‐century earthquake sequence in the Segura Medio valley. These findings provide valuable insights for seismic hazard modeling in the region.

Microseismicity and fault structure in the Daliangshan block within the Southeastern Tibetan Plateau

The Daliangshan block is located between the Tibetan Plateau and the South China block and has accommodated several M > 6.5 damaging earthquakes in the past ∼600 years, as well as intense tectonic deformation and complex fault structures. In this study, we analyze more than two years of continuous seismic data recorded by a recently deployed dense seismic array. We used a recently developed machine learning-based earthquake location workflow (ESPRH) to construct a high-precision earthquake catalog for the region and obtained 3539 earthquakes, which is approximately three times as many as the National Earthquake Data Center (NEDC) catalog contains. The seismicity distribution not only confirms the nature of the faults marked on the map but also delineates the detailed geometry of the unmapped faults, including the en échelon faults at the northern end of the Zemuhe Fault and the “V”-shaped conjugate fault within the Mabian Fault Zone. The Zemuhe Fault and Lianfeng Fault are prone to hosting large earthquakes according to the derived low b-value. The western side of the Daliangshan block is dominated by strike-slip faults. Combining the fault geometry presented in this paper, we observed that the fault properties on the eastern side are complex. This tectonic phenomenon is attributed to the fact that during the lateral extrusion of the southeastern edge of the Chuandian fragments, the northeastern part of the Daliangshan block was squeezed by the South China block more strongly than its southwestern part. We provide the first precise earthquake catalog for Daliangshan block, which can be used as important seismological data for regional hazard assessment and research on the southeastern (SE) margin of the Tibetan Plateau.

Large‐Scale Extensional Strain in Southern Tibet From Sentinel‐1 InSAR and GNSS Data

AbstractIn this study, we utilize C‐band Sentinel‐1 radar images from 2015 to 2022, combined with interseismic horizontal GNSS velocities, to construct large‐scale, high‐resolution, 3‐D velocity and strain rate maps over a vast region of southern Tibet. We show the distribution of prevailing dilatational strain accumulation along the seven major rift zones. Using 2‐D elastic dislocations invoking a two‐fault model in a Bayesian framework, we quantified the decadal extension rates across the seven rift zones, and we suggest a total extension rate of 18.4 ± 1.7 mm/yr, consistent with geological and geodetic estimates. The resulting strain rate maps, combined with the earthquake catalog, help us identify areas with high earthquake potential. Our study enhances our understanding of the present‐day tectonics and kinematics in southern Tibet and provides important constraints for seismic hazard assessment in this region.

On modeling seismicity in seismic hazard assessment problems

Seismicity modeling is an important part of creating General Seismic Zoning maps within the framework of a probabilistic approach. We consider the main disadvantages of individual elements of the recent seismicity models. A variant of the methodology is proposed, which, due to the improvements of those elements, should provide more accurate estimates of the future seismicity. For the first time, a stochastic seismicity model has been proposed in the form of a synthetic earthquake catalog, generated for an arbitrary conditional period and reproducing the properties of the catalog of actual earthquakes, including spatiotemporal grouping. A methodology for verifying seismicity models is proposed to check the compliance of the models with the initial data, to assess the predictive efficiency of the models, and to compare efficiency of different models.

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

Using a 1D Radially Symmetric Coda Envelope Model for Robust Moment Magnitude (Mw) Estimation in Iraq’s Tectonically Diverse Zones

ABSTRACT Robust estimation of moment magnitude (Mw) can be challenging for Iraq due to the strong lateral heterogeneity across diverse tectonic zones. We aim to improve moment magnitude estimation by investigating the reliability of using a 1D coda envelope model in diverse tectonic zones of different lateral effects and offer a way forward for reliable estimates of Mw for small events that are difficult to waveform model. Iraq comprises two main tectonic zones: (1) the Outer platform, consisting of the northwestern Zagros fold-thrust belt and the Mesopotamian foredeep, and (2) the Inner Arabian platform which is overlain by the Iraqi desert. A simple 1D coda envelope model was used because coda waves have a low sensitivity to the source and path heterogeneity. Three separate coda calibrations were conducted to investigate the robustness of a single 1D calibration to fit the country: Whole-region calibration, Zagros calibration, and Mesopotamia calibration. In the whole-region calibration, we used stations from both the Zagros and Mesopotamia zones. In the two other calibration models, we used only stations that were in those particular zones. Ground-truth reference spectra derived from the coda spectral ratio method were used to constrain high-frequency site terms. There was no drastic difference when comparing the moment magnitudes calculated from the waveform modeling and the three calibration models. The results show that the 1D coda envelope model is a reliable method even for a region with diverse tectonic zones. Hence, we recommend using the whole-region calibration model for moment magnitude estimation that provides more complete path coverage and avoids biases introduced by path correction failures. The proposed calibration is a fundamental step in updating the comprehensive earthquake catalog and probabilistic seismic hazard assessments for Iraq.