More Like this

1. Teleseismic Converted Waves Image of the Indo-Burma Subduction System Across the Bangladesh-India-Myanmar (BIMA) Array

   Ajala, R.; Persaud, P.; Steckler, M.; Sandvol, E.; Islam, M.; Gaherty, J.; Than, O.; Oo, K.; Akhter, S.; Ni, J. (January 2020, American Geophysical Union, Fall Meeting 2020)

   null (Ed.)

   Recent GPS studies show that the Indo-Burma subduction system is locked with the implication of a potential large-magnitude earthquake. To inform better seismic hazard models in the region, we need an improved understanding of the crustal structure and the dynamics of the Indo-Burma subduction system. The Bangladesh-India-Myanmar (BIMA) tripartite project deployed 60 broadband seismometers across the subduction system and have been continuously recording data for ~2 years. In this study, we computed receiver functions from 30 high-quality earthquakes (M≥5.9) with epicentral distances between 30º and 90º recorded by the array. The algorithm utilized ensures the uniqueness of the seismic model and provides an uncertainty estimate of every converted wave amplitude. We stacked all the receiver functions produced at each station along the entire transect to generate a cross-sectional model of the average crustal structure. The level of detail in the image is improved by computing higher frequency receiver functions up to 4 Hz. The results represent some of the strongest constraints on crustal structure across the subduction system. Beneath the Neogene accretionary prism's outer belt, we observe a primary conversion associated with the Ganges Brahmaputra Delta that ranges in depth from ~10 km near the deformation front up to ~12 km at the eastern boundary. From the eastern end of the Neogene accretionary prism to the Sagaing Fault, we image the Indian subducting slab and the Central Myanmar basin. The depth-extent of seismicity associated with the Wadati-Benioff zone is consistent with the locations of primary conversions from the subducting plate. We further verify the converted phases of the slab by analyzing azimuthal moveout variations. The Central Myanmar basin is roughly bowl-shaped in cross-section with a maximum thickness of ~15 km about halfway between the Kabaw and Sagaing faults. The average crustal thickness beneath the Ganges-Brahmaputra delta is ~20 km, most likely representing a transitional crust formed from thinning of the continental crust intruded and underplated by igneous rocks. In contrast, the average thickness of the continental crust beneath the Central Myanmar basin is ~40 km. Our results provide a baseline model for future geophysical investigations of the Indo-Burma subduction zone. 

   more » « less
2. Estimating site amplification variability in Yangon, Myanmar, from a dense nodal seismic array

   https://doi.org/10.1177/87552930251352843  

   Islam, Md_Mohimanul; Persaud, Patricia; Thant, Myo; Win, Khin; Sandvol, Eric (July 2025, Earthquake Spectra)

   We explore a novel acquisition geometry that can be used to estimate the linear component of site amplification using a dense nodal seismic network installed in Yangon, Myanmar’s largest city. The city is surrounded by several seismically active faults, including the Sagaing Fault, which is capable of generating M_w > 7.0 earthquakes. As part of the Irrawaddy delta system, this densely populated city sits on young water-saturated alluvium that is likely to amplify earthquake ground motions. Assessing site response is crucial for understanding the seismic hazard potential to minimize the loss of property and lives. Using a dense seismic array comprised of 110 three-component nodes, we estimated the frequency-dependent site amplification pattern of Yangon from regional (Lg) and local (Sg) seismic phases. Since this acquisition geometry is not sensitive to Q or geometric spreading, this approach provides a fast and cost-effective way to estimate the linear component of site response as a function of frequency. Our Lg and Sg site response results identify regions with high site amplification that have significantly greater seismic hazard risks for regional and local distance earthquakes. We observed consistent site response characteristics between both Lg and Sg phases. Site amplification patterns correlate well with the surficial geology and subsurface structure beneath the city. De-amplification is observed across all frequencies at stations located above an anticlinal structure composed of older Pliocene rocks (i.e. the Irrawaddy Formation). Conversely, highly amplified areas correspond to younger Pleistocene to recent alluvial plains consisting of loose, unconsolidated alluvium. We found a dominant horizontal-to-vertical spectral ratio (HVSR) peak at ∼1.0 Hz from ambient noise, likely corresponding to the thickness of unconsolidated sediments. We suggest that the growing number of nodal networks worldwide can be used to estimate frequency-dependent site amplification, addressing key data gaps in seismic hazard assessment. 

   more » « less
3. New Insights into the Crustal Structure of the San Fernando Valley, California, from a Dense Nodal Seismic Array

   https://doi.org/10.1785/0220240473  

   Juárez-Zúñiga, Alan; Persaud, Patricia (May 2025, Seismological Research Letters)

   Abstract The San Fernando Valley (SFV), a densely populated region in Southern California, has high earthquake hazard due to a complex network of active faults and the amplifying effects of the sedimentary basin. Since the devastating 1994 Mw 6.7 Northridge earthquake, numerous studies have examined its structure using various geological and geophysical datasets. However, current seismic velocity models still lack the resolution to accurately image the near-surface velocity structure and concealed or blind faults, which are critical for high-frequency wavefield simulations and earthquake hazard modeling. To address these challenges, we develop a 3D high-resolution shear-wave velocity model for the SFV using ambient noise data from a dense array of 140 seismic nodes and 10 Southern California Seismic Network stations. We also invert gravity data to map the basin geometry and integrate horizontal-to-vertical spectral ratios and aeromagnetic data to constrain interfaces and map major geological structures. With a lateral resolution of 250 m near the basin center, our model reveals previously unresolved geological features, including the detailed geometry of the basin and previously unmapped structure of faults at depth. The basin deepens from the Santa Monica Mountains in the south to approximately 4 km near its center and 7 km in the Sylmar sub-basin at the basin’s northern margin. Strong velocity contrasts are observed across major faults, at the basin edges, and in the basin’s upper 500 m, for which we measure velocities as low as 200 m/s. Our high-resolution model will enhance ground-motion simulations and earthquake hazard assessments for the SFV and has implications for other urban areas with high seismic risk. 

   more » « less
4. Dense Seismic Recordings of the 2023 Kahramanmaraş Earthquake Sequence in Southeastern Türkiye

   https://doi.org/10.1785/0220240152  

   Zor, Ekrem; Peng, Zhigang; Ergin, Mehmet; Sandvol, Eric; Sevim, Fatih; Tapırdamaz, Mustafa Cengiz; Mach, Phuc; Yalvaç, Oğuzhan; Tarancıoğlu, Adil; Koşma, Mustafa Kemal; et al (February 2025, Seismological Research Letters)

   Abstract The devastating 6 February 2023 Kahramanmaraş earthquake sequence in southeastern Türkiye started with a moment magnitude (Mw) 7.8 earthquake, for which the initial rupture broke the Sakçagöz segment near Nurdağı and then jumped into a bilateral rupture along multiple segments of the Eastern Anatolian fault zone (EAFZ). This complicated rupture was followed nine hours later by an Mw 7.6 event near Ekinözü. To better understand the spatiotemporal evolution of aftershocks, site amplification, and the structural and tectonic framework of the EAFZ in this diffuse triple junction, we deployed a dense seismometer array covering both aftershock zones for nearly four months. The main Eastern Anatolian Seismic Temporary (EAST) array includes 125 nodal, 10 broadband, and 6 strong-motion seismic stations distributed around the rupture zone. An additional linear array of 73 nodal stations was also installed across the Pazarcık segment of the EAFZ and the Sakçagöz segment near the Mw 7.8 epicenter to record fault-zone waves for ∼30 days. This article shows example recordings and the EAST array geometry, preliminary research results, and the metadata related to all of the stations in this array. A deep-learning-based phase picking for one month of continuous recording yielded millions of seismic phase readings and tens of thousands of aftershock locations after phase associations. We also give examples of both local and teleseismic waveforms recorded by the nodal arrays, which can be used for subsequent high-resolution earthquake relocation, imaging of crustal structures, and fault-zone imaging. 

   more » « less
5. Seismic Hazard Analyses From Geologic and Geomorphic Data: Current and Future Challenges

   https://doi.org/10.1029/2018TC005365  

   Morell, Kristin_D; Styron, Richard; Stirling, Mark; Griffin, Jonathan; Archuleta, Ralph; Onur, Tuna (October 2020, Tectonics)

   Abstract The loss of life and economic consequences caused by several recent earthquakes demonstrate the importance of developing seismically safe building codes. The quantification of seismic hazard, which describes the likelihood of earthquake‐induced ground shaking at a site for a specific time period, is a key component of a building code, as it helps ensure that structures are designed to withstand the ground shaking caused by a potential earthquake. Geologic or geomorphic data represent important inputs to the most common seismic hazard model (probabilistic seismic hazard analyses, or PSHAs), as they can characterize the magnitudes, locations, and types of earthquakes that occur over long intervals (thousands of years). However, several recent earthquakes and a growing body of work challenge many of our previous assumptions about the characteristics of active faults and their rupture behavior, and these complexities can be challenging to accurately represent in PSHA. Here, we discuss several of the outstanding challenges surrounding geologic and geomorphic data sets frequently used in PSHA. The topics we discuss include how to utilize paleoseismic records in fault slip rate estimates, understanding and modeling earthquake recurrence and fault complexity, the development and use of fault‐scaling relationships, and characterizing enigmatic faults using topography. Making headway in these areas will likely require advancements in our understanding of the fundamental science behind processes such as fault triggering, complex rupture, earthquake clustering, and fault scaling. Progress in these topics will be important if we wish to accurately capture earthquake behavior in a variety of settings using PSHA in the future. 

   more » « less