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. Abstract The along‐strike variations of the velocity, thickness, and dip of subducting slabs and the volcano distribution have been observed globally. It is, however, unclear what controls the distribution of volcanoes and the associated magma generation. With the presence of nonuniform volcanism, the Aleutian‐Alaska subduction zone (AASZ) is an ideal place to investigate subduction segmentation and its relationship with volcanism. Using full‐wave ambient noise tomography, we present a high‐resolution 3‐D shear wave velocity model of the AASZ for the depths of 15–110 km. The velocity model reveals the distinct high‐velocity Pacific slab, the thicker, flatter, and more heterogeneous Yakutat slab, and the northeasterly dipping Wrangell slab. We observe low velocities within the uppermost mantle (at depth <60 km) below the Aleutian arc volcanoes, representing partial melt accumulation. The large crustal low‐velocity anomaly beneath the Wrangell volcanic field suggests a large magma reservoir, likely responsible for the clustering of volcanoes. The Denali volcanic gap is above an average‐velocity crust but an extremely fast mantle wedge, suggesting the lack of subsurface melt. This is in contrast with the lower‐velocity back‐arc mantle beneath the adjacent Buzzard Creek‐Jumbo Dome volcanoes to the east. The back‐arc low velocities associated with the Pacific, the eastern Yakutat, and the Wrangell slabs may reflect subduction‐driven mantle upwelling. The structural variation of the downgoing slabs and the overriding plate explains the change of volcanic activity along the AASZ. Our findings demonstrate the combined role of the subducting slab and the overriding plate in controlling the characteristics of arc magmatism. 
    more » « less
  4. Abstract The impact of past tectonic events on the formation and modification of continental lithosphere remains as an open question of fundamental importance. Eastern North America provides a complete record of supercontinent assembly and breakup over the past 1.3 Ga, serving as a natural laboratory for our understanding of continental crust and mantle lithosphere and for integrating geologic and geophysical observations. In this study, we used teleseismic Ps receiver functions to image the detailed distribution of crustal thickness beneath eastern North America. The radial‐component receiver functions were calculated from seismic waveforms recorded by a total of 659 broadband stations during 2010–2017, yielding a high‐resolution image of Moho depth distribution. The depths of the Moho and intracrustal layers vary within and across the major tectonic units. Specifically, there are distinct differences in crustal thickness between the northern and southern Grenville Province. A dipping intracrustal feature can be seen within the central Grenville Province, with the depth increasing eastward from 5 to 27 km. The Moho depth decreases southeastward across the Grenville‐Appalachian boundary, with a sharp Moho offset of up to 12–15 km in the central segment and a more gradual variation to the north and south. The thickness difference between the southern and northern Grenville‐aged crusts suggests different tectonic and/or exhumation histories during and after the Grenville Orogeny. The low‐angle eastward dipping crustal feature is interpreted to be a Grenville‐aged collisional structure. Differences in the steepness of the Moho offset along the strike of Appalachians probably reflect variation of the steepness of the subsurface boundary between Laurentia and accreted terranes with different intensities of postorogenic modification. The observed spatial relation between the geologically defined tectonic boundaries and crustal thickness variations provides new constraints on the depth extent of the tectonic units within the crust. 
    more » « less
  5. 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