We present electrical resistivity models of the crust and upper mantle from two‐dimensional (2‐D) inversion of magnetotelluric (MT) data collected in the Rio Grande rift, New Mexico, USA. Previous geophysical studies of the lithosphere beneath the rift identified a low‐velocity zone several hundred kilometers wide, suggesting that the upper mantle is characterized by a very broad zone of modified lithosphere. In contrast, the surface expression of the rift (e.g., high‐angle normal faults and synrift sedimentary units) is confined to a narrow region a few tens of kilometers wide about the rift axis. MT data are uniquely suited to probing the depths of the lithosphere that fill the gap between surface geology and body wave seismic tomography, namely the middle to lower crust and uppermost mantle. We model the electrical resistivity structure of the lithosphere along two east‐west trending profiles straddling the rift axis at the latitudes of 36.2 and 32.0°N. We present results from both isotropic and anisotropic 2‐D inversions of MT data along these profiles, with a strong preference for the latter in our interpretation. A key feature of the anisotropic resistivity modeling is a broad (~200‐km wide) zone of enhanced conductivity (<20 Ωm) in the middle to lower crust imaged beneath both profiles. We attribute this lower crustal conductor to the accumulation of free saline fluids and partial melt, a direct result of magmatic activity along the rift. High‐conductivity anomalies in the midcrust and upper mantle are interpreted as fault zone alteration and partial melt, respectively.
We develop and apply a technique for receiver functions that is analogous to “velocity analysis” in seismic reflection processing, in which a velocity model is found directly from the data. In the case of receiver functions, which represent waves converted from
- NSF-PAR ID:
- 10370924
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Geochemistry, Geophysics, Geosystems
- Volume:
- 20
- Issue:
- 2
- ISSN:
- 1525-2027
- Page Range / eLocation ID:
- p. 974-991
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
more » « less
Abstract We conduct a joint inversion of teleseismic receiver functions and Rayleigh wave phase velocity dispersion from both ambient noise and earthquakes using data from 79 seismic stations in southern Africa, which is home to some of the world's oldest cratons and orogenic belts. The area has experienced two of the largest igneous activities in the world (the Okavango dyke swarm and Bushveld mafic intrusion) and thus is an ideal locale for investigating continental formation and evolution. The resulting 3‐D shear wave velocities for the depth range of 0–100 km and crustal thickness measurements show a clear spatial correspondence with known geological features observed on the surface. Higher than normal mantle velocities found beneath the southern part of the Kaapvaal craton are consistent with the basalt removal model for the formation of cratonic lithosphere. In contrast, the Bushveld complex situated within the northern part of the craton is characterized by a thicker crust and higher crustal
V p /V s but lower mantle velocities, which are indicative of crustal underplating of mafic materials and lithospheric refertilization by the world's largest layered mafic igneous intrusion. The thickened crust and relatively low elevation observed in the Limpopo belt, which is a late Archean collisional zone between the Kaapvaal and Zimbabwe cratons, can be explained by eclogitization of the basaltic lower crust. The study also finds evidence for the presence of a stalled segment of oceanic lithosphere beneath the southern margin of the Proterozoic Namaqua‐Natal mobile belt. -
Subduction of the very young (<15 Myr old) oceanic lithosphere of the Nazca plate in central to southern Colombia is observationally related to an unusually high and unusually variable amount of intermediate (>50 km) depth seismicity. From 2010 through 2019 89% of central and southern Colombia’s 11,466 intermediate depth events occurred between 3.5°N and 5.5°N, highlighting these unusual characteristics of the young slab. In addition, morphologic complexity and possible tears characterize the Nazca slab in Colombia and complicate mantle flow in the region. Prior SKS-phase shear-wave splitting results indicate sub-slab anisotropy is dominated by plate motion parallel-to-subparallel orientations in the region, suggesting the young slab has entrained a relatively thick portion of the sub-slab mantle. These observations suggest the subduction of young lithosphere has significant effects on both the overlying and underlying asthenosphere in the Colombia subduction zone. Here we use more than 10 years of data to calculate receiver functions for the Red Sismológica Nacional de Colombia’s network of broadband seismometers. These receiver functions allow us to tie these prior observations of the Colombia subduction zone to distinct, structural features of the slab. We find that the region of high seismicity corresponds to a low seismic velocity feature along the top of the subducting plate between 3.5°N and 5.5°N that is not present to the south. Moderately elevated P-wave velocity to S-wave velocity ratios are also observed within the slab in the north. This feature likely represents hydrated slab mantle and/or uneclogitized oceanic crust extending to a deeper depth in the north of the region which may provide fluids to drive slab seismicity. We further find evidence for a thick layer of material along the slab’s lithosphere-asthenosphere boundary characterized by spatially variable anisotropy. This feature likely represents entrained asthenosphere at the base of the plate sheared by both the overlying plate and complex flow related to proposed slab tears just north and south of the study region. These observations highlight how structural observations provide key contextual constraints on short-term (seismogenic) and long-term (anisotropic fabric) dynamic processes in the Colombia subduction zone. Plain-language Summary The Nazca oceanic plate is very young (<15 million years old) where it is pulled or subducted beneath the South America plate in central and southern Colombia. Earthquakes occurring in the subducted Nazca plate at depths greater than 50 km are nearly 9x more common in central Colombia than in southern Colombia. The subducted Nazca plate also has a complex shape in this region and may have been torn both in northern Colombia and to the south near the Colombia-Ecuador border. The slow flow of mantle rock beneath the subducted plate is believed to be affected by this and earlier studies have inferred this flow is mostly in the same direction as the subducting plate's motion. We have used 10+ years of data to calculate receiver functions, which can detect changes in the velocity of seismic waves at the top and bottom of the subducted plate to investigate these features. We found that the Nazca plate is either hydrated or has rocks with lower seismic velocities at its top in the central part of Colombia where earthquakes are common. We also find that a thick layer of mantle rock at the base of the subducted plate has been sheared.more » « less
-
more » « less
SUMMARY To evaluate the plate flexure model for the formation of the Appalachian Basin, we investigate the extent to which crustal structure beneath and surrounding the basin was modified by the Palaeozoic orogenic events that created the basin. We jointly invert receiver functions and surface wave dispersion measurements to obtain 1-D crustal Vs profiles for 261 seismic stations located within and around the basin. The average crustal thickness for the region is 44 km, and the crust gradually thins to the east, consistent with previous studies. Four areas of anomalous crust are identified with respect to the eastward thinning of the crust. An area of thick crust is found along the Grenville Front on the western side of the Appalachian Basin where the crust thickens by ∼5–10 km. Moho depths of up to 54 km in this region likely result from suture-thickening. The crust is thinner beneath the Neoproterozoic Scranton rift by ∼5–7 km, coincident with a ∼40 mGal Bouguer gravity high. Across the Neoproterozoic Rome Trough, the crust thins by ∼4–5 km, coincident with a ∼10 mGal Bouguer gravity high. Density models for these rifts show that the rift-related crustal thinning is sufficient to explain the gravity anomalies. The Vs models obtained for stations in the rifts indicate little, if any, mafic layering in the mid-crust and only a modest amount of mafic layering in the lower crust. In the northwestern portion of the Appalachian Basin in northeastern Ohio and northwestern Pennsylvania within the Elzevir block, another area of anomalously thick crust (50–52 km) is found. This region is not associated with any known tectonic structures or boundaries or a gravity anomaly. The lower ∼5–10 km of the crust in this region is characterized by high (>3.9 km s−1) shear wave velocities and thus appears to be mafic. The origin of anomalous crustal structure in all four areas is best attributed to Precambrian tectonic events that predate the formation of the Appalachian Basin, indicating that the crystalline crust beneath and surrounding the basin was not significantly affected by the Palaeozoic basin-forming orogenic events, a finding which supports the use of plate flexure models for understanding basin formation.
-
more » « less
Abstract The East African Rift System provides a rare location in which to observe a wide scope of rifting states. Well‐defined active narrow rifting in the Main Ethiopian Rift (MER) transitions to incipient extension and eventually pre‐rifted lithosphere through the northwestern flank of the Ethiopian Plateau (EP). Although the MER is well studied, the off‐axis region has received less attention. We develop Rayleigh wave phase velocity maps,
Ps receiver functions, and H‐κ stack surfaces, and jointly invert these data using a trans‐dimensional, hierarchical Bayesian inversion algorithm to create shear velocity profiles across the MER and EP. All shear velocities observed are slower than the PREM global average, a reflection of the elevated temperatures that persist from plume impingement. In the EP, we find a shallow mantle slow shear velocity lineament parallel to the MER axis, amidst otherwise faster shear velocities. The crust is shallow in the MER, and also in the northwestern‐most EP flank. Thicker crust found elsewhere throughout the plateau is caused by crustal underplating and flood basalt emplacement. Shear velocities more reduced than the already low regional average, in concert with surficial volcanic features, geodetic observations, and slow P‐ and S‐wave anomalies, support off‐axis extension in the Ethiopian plateau, requiring reevaluation of the localization of continental breakup in the narrow MER.
