skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Crustal structure variations in south-central Mexico from receiver functions
Summary Mexico has a complex geological history that is typified by the distinctive terranes that are found in the south-central region. Crustal thickness variations often correlate with geological terranes that have been altered by several processes in the past, for example aerial or subduction erosion, underplating volcanic material or rifting but few geophysical studies have locally imaged the entire continental crust in Mexico. In this paper, the thickness of three layers of the crust in south-central Mexico is determined. To do this, we use P- and S-wave receiver functions (RF) from 159 seismological broad-band stations. Thanks to its adaptive nature, we use an empirical mode decomposition (EMD) algorithm to reconstruct the RFs into intrinsic mode functions (IMF) in order to enhance the pulses related to internal discontinuities within the crust. To inspect possible lateral variations, the RFs are grouped into quadrants of 90°, and their amplitudes are mapped into the thickness assuming a three-layer model. Using this approach, we identify a shallow sedimentary layer with a thickness in the range of 1–4 km. The upper-crust was estimated to be of a few kilometers (<10 km) thick near the Pacific coast, and thicker, approximately 15 km in central Oaxaca and under the Trans-Mexican Volcanic Belt (TMVB). Close to the Pacific coast, we infer a thin crust of approximately 16 ± 0.9 km, while in central Oaxaca and beneath the TMVB, we observe a thicker crust ranging between 30 and 50 km ± 2.0 km. We observe a crustal thinning, of approximately 6 km, from central Oaxaca (37 ± 1.9 km) towards the Gulf of Mexico, under the Veracruz Basin, where we estimate a crustal thickness of 31.6 ± 1.9 km. The boundary between the upper and lower crust in comparison with the surface of the Moho do not show significant variations other than the depth difference. We observe small crustal variations across the different terranes on the study area, with the thinnest crust located at the Pacific coast and Gulf of Mexico coast. The thickest crust is estimated to be in central Oaxaca and beneath the TMVB.  more » « less
Award ID(s):
1645063
PAR ID:
10121955
Author(s) / Creator(s):
 ;  ;  ;  ;  
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Geophysical Journal International
Volume:
219
Issue:
3
ISSN:
0956-540X
Page Range / eLocation ID:
p. 2174-2186
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract The crustal structure in south‐central Alaska has been influenced by terrane accretion, flat slab subduction, and a modern strike‐slip fault system. Within the active subduction system, the presence of the Denali Volcanic Gap (DVG), a ∼400 km region separating the active volcanism of the Aleutian Arc to the west and the Wrangell volcanoes to the east, remains enigmatic. To better understand the regional tectonics and the nature of the volcanic gap, we deployed a month‐long north‐south linear geophone array of 306 stations with an interstation distance of 1 km across the Alaska Range. By calculating multi‐component noise cross‐correlation and jointly inverting Rayleigh wave phase velocity and ellipticity across the array, we construct a 2‐D shear wave velocity model along the transect down to ∼16 km depth. In the shallow crust, we observe low‐velocity structures associated with sedimentary basins and image the Denali fault as a narrow localized low‐velocity anomaly extending to at least 12 km depth. About 12 km, below the fold and thrust fault system in the northern flank of the Alaska Range, we observe a prominent low‐velocity zone with more than 15% velocity reduction. Our velocity model is consistent with known geological features and reveals a previously unknown low‐velocity zone that we interpret as a magmatic feature. Based on this feature's spatial relationship to the Buzzard Creek and Jumbo Dome volcanoes and the location above the subducting Pacific Plate, we interpret the low‐velocity zone as a previously unknown subduction‐related crustal magma reservoir located beneath the DVG. 
    more » « less
  2. SUMMARY The Ecuadorian forearc is a complex region of accreted terranes with a history of large megathrust earthquakes. Most recently, a Mw 7.8 megathrust earthquake ruptured the plate boundary offshore of Pedernales, Ecuador on 16 April 2016. Following this event, an international collaboration arranged by the Instituto Geofisico at the Escuela Politécnica Nacional mobilized a rapid deployment of 65 seismic instruments along the Ecuadorian forearc. We combine this new seismic data set with 14 permanent stations from the Ecuadorian national network to better understand how variations in crustal structure relate to regional seismic hazards along the margin. Here, we present receiver function adaptive common conversion point stacks and a shear velocity model derived from the joint inversion of receiver functions and surface wave dispersion data obtained through ambient noise cross-correlations for the upper 50 km of the forearc. Beneath the forearc crust, we observe an eastward dipping slow velocity anomaly we interpret as subducting oceanic crust, which shallows near the projected centre of the subducting Carnegie Ridge. We also observe a strong shallow positive conversion in the Ecuadorian forearc near the Borbon Basin indicating a major discontinuity at a depth of ∼7 km. This conversion is not ubiquitous and may be the top of the accreted terranes. We also observe significant north–south changes in shear wave velocity. The velocity changes indicate variations in the accreted terranes and may indicate an increased amount of hydration beneath the Manabí Basin. This change in structure also correlates geographically with the southern rupture limit of multiple high magnitude megathrust earthquakes. The earthquake record along the Ecuadorian trench shows that no event with a Mw >7.4 has ruptured south of ∼0.5°S in southern Ecuador or northern Peru. Our observations, along with previous studies, suggest that variations in the forearc crustal structure and subducting oceanic crust may influance the occurrence and spatial distribution of high magnitude seismicity in the region. 
    more » « less
  3. The New England Appalachians provide a fascinating window into a host of fundamental geological problems. These include the modification of crustal and mantle lithospheric structure via orogenesis, terrane accretion, and continental rifting, the evolution of individual terranes through processes such as channel flow and ductile extrusion, and the causes and consequences of the Northern Appalachian Anomaly (NAA), a prominent geophysical anomaly in the upper mantle. Recent and ongoing deployments of dense seismic arrays in New England are providing images of the crust and upper mantle in unprecedented detail, allowing us to address both new and longstanding science questions. These deployments include the Seismic Experiment for Imaging Structure beneath Connecticut (SEISConn, 2015-2019), the New England Seismic Transects (NEST, 2018-present), and the GEology of New England via Seismic Imaging Studies (GENESIS, 2022-present) arrays. Here we present results from these experiments that are shedding new light on the tectonic evolution of New England and the ways in which structures and processes in the upper mantle can affect the structure of the overlying lithosphere. These include detailed new images of crustal architecture beneath central and southern New England, including a sharp transition from thick (~48 km) crust Laurentia terranes to thin (~32 km) crust beneath Appalachian terranes. The character of this offset beneath the SEISConn and NEST arrays suggests an overlap of two Moho boundaries, forming an overthrust-type structure that may have resulted from reactivation of faults during the compression and shortening associated with the formation of the hypothesized Acadian Altiplano. Beneath SEISConn, there is evidence for multiple relict structures preserved in the lithosphere from past episodes of terrane accretion and suturing, as well as anisotropic layering that constrains the kinematics of past lithospheric deformation events. Beneath the NEST line in central New England, we infer a relatively shallow (~80 km) lithosphere-asthenosphere boundary above the NAA upper mantle geophysical anomaly, providing evidence for lithospheric thinning above a presumed asthenospheric upwelling. Finally, preliminary results suggest layered crustal anisotropy beneath the GENESIS array, perhaps corresponding to a past episode of channel flow in the mid-crust. 
    more » « less
  4. The New England Appalachians provide a fascinating window into a host of fundamental geological problems. These include the modification of crustal and mantle lithospheric structure via orogenesis, terrane accretion, and continental rifting, the evolution of individual terranes through processes such as channel flow and ductile extrusion, and the causes and consequences of the Northern Appalachian Anomaly (NAA), a prominent geophysical anomaly in the upper mantle. Recent and ongoing deployments of dense seismic arrays in New England are providing images of the crust and upper mantle in unprecedented detail, allowing us to address both new and longstanding science questions. These deployments include the Seismic Experiment for Imaging Structure beneath Connecticut (SEISConn, 2015-2019), the New England Seismic Transects (NEST, 2018-present), and the GEology of New England via Seismic Imaging Studies (GENESIS, 2022-present) arrays. Here we present results from these experiments that are shedding new light on the tectonic evolution of New England and the ways in which structures and processes in the upper mantle can affect the structure of the overlying lithosphere. These include detailed new images of crustal architecture beneath central and southern New England, including a sharp transition from thick (~48 km) crust Laurentia terranes to thin (~32 km) crust beneath Appalachian terranes. The character of this offset beneath the SEISConn and NEST arrays suggests an overlap of two Moho boundaries, forming an overthrust-type structure that may have resulted from reactivation of faults during the compression and shortening associated with the formation of the hypothesized Acadian Altiplano. Beneath SEISConn, there is evidence for multiple relict structures preserved in the lithosphere from past episodes of terrane accretion and suturing, as well as anisotropic layering that constrains the kinematics of past lithospheric deformation events. Beneath the NEST line in central New England, we infer a relatively shallow (~80 km) lithosphere-asthenosphere boundary above the NAA upper mantle geophysical anomaly, providing evidence for lithospheric thinning above a presumed asthenospheric upwelling. Finally, preliminary results suggest layered crustal anisotropy beneath the GENESIS array, perhaps corresponding to a past episode of channel flow in the mid-crust. 
    more » « less
  5. 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