Abstract The North Anatolian Fault Zone (NAFZ) is a prominent tectonic structure with a significant impact on the observed active deformation in Türkiye. Detailed knowledge of the seismic anisotropy in the crust and mantle along this nascent shear deformation zone provides insights into the kinematics associated with past and present tectonic events. We employed teleseismic earthquakes observed by the Dense Array North Anatolia seismic network to map 3‐ D variations in crustal and mantle anisotropy in/around the NW segment of the NAFZ. To achieve this, we first performed a harmonic decomposition analysis of P‐receiver functions. The results were then used as a priori information to conduct an anisotropic receiver function inversion with the Neighborhood Algorithm that enabled imaging of the actual orientation and geometry of anisotropic structures. SKS splitting measurements are further used to make a comparison between the anisotropic behavior of crustal and mantle structures. Crustal anisotropy parameters estimated in our analyses/models well identify the signature of deformation caused by accumulated strain in the earthquake cycle through the strike of shallow cutting faults in the brittle crust beneath the NAFZ. Diffuse intense anisotropic energy at lower crustal depths was attributed to lattice preferred orientation of crystals or partially molten lenses elongated along the shear direction. Strong harmonic energy variations beneath the northern part of the Istanbul Zone likely reflect imprints of LPO‐originated frozen fabric at shallow depths (0–20 km) associated with the palaeotectonic Odessa Shelf, Intra‐Pontide Suture Zones or remnants of the Tethys Ocean.
SUMMARY Information on fault zone structure is essential for our understanding of earthquake mechanics, continental deformation and seismic hazard. We use the scattered seismic wavefield to study the subsurface structure of the North Anatolian Fault Zone (NAFZ) in the region of the 1999 İzmit and Düzce ruptures using data from an 18-month dense deployment of seismometers with a nominal station spacing of 7 km. Using the forward- and back-scattered energy that follows the direct P-wave arrival from teleseismic earthquakes, we apply a scattered wave inversion approach and are able to resolve changes in lithospheric structure on a scale of 10 km or less in an area of about 130 km by 100 km across the NAFZ. We find several crustal interfaces that are laterally incoherent beneath the surface strands of the NAFZ and evidence for contrasting crustal structures either side of the NAFZ, consistent with the presence of juxtaposed crustal blocks and ancient suture zones. Although the two strands of the NAFZ in the study region strike roughly east–west, we detect strong variations in structure both north–south, across boundaries of the major blocks, and east–west, parallel to the strike of the NAFZ. The surface expression of the two strands of the NAFZ is coincident with changes on main interfaces and interface terminations throughout the crust and into the upper mantle in the tomographic sections. We show that a dense passive network of seismometers is able to capture information from the scattered seismic wavefield and, using a tomographic approach, to resolve the fine scale structure of crust and lithospheric mantle even in geologically complex regions. Our results show that major shear zones exist beneath the NAFZ throughout the crust and into the lithospheric mantle, suggesting a strong coupling of strain at these depths.
Abstract Eighty‐two broadband seismic stations of the Superior Province Rifting Earthscope Experiment (SPREE) collected 2.5 years of continuous seismic data in the area of the high gravity anomaly associated with the Midcontinent Rift (MCR). Over 100 high‐quality teleseismic earthquakes were used for crustal P wave receiver function analysis. Our analysis reveals that the base of the sedimentary layer is shallow outside the MCR, thickens near the flanks where gravity anomalies are low, and shallows again in the MCR's center where the gravity anomalies peak. This pattern is similar to that found from local geophysical studies and is consistent with reverse faulting having accompanied the cessation of rifting at 1.1 Ga. Intermittent intracrustal boundaries imaged by our analysis might represent the bottom of the MCR's mostly buried dense volcanic layers. Outside the MCR, the Moho is strong, sharp, and relatively flat, both beneath the Archean Superior Province and the Proterozoic terranes to its south. Inside the MCR, two weaker candidate Mohos are found at depths up to 25 km apart in the rift's center. The intermediate layer between these discontinuities tapers toward the edges of the MCR. The presence of this transitional layer is remarkably consistent along the strike of the MCR, including beneath its jog in southern Minnesota, near the Belle Plaine Fault. We interpret these results as evidence for extensive underplating as a defining characteristic of the rift, which remains continuous along the Minnesota jog, where due to its orientation, it is minimally affected by the reverse faulting that characterizes the NNE striking parts of the rift.
We present shear-wave splitting analyses of SKS and SKKS waves recorded at sixteen Superior Province Rifting Earthscope Experiment (SPREE) seismic stations on the north shore of Lake Superior, as well as fifteen selected Earthscope Transportable Array instruments south of the lake. These instruments bracket the Mid-Continent Rift (MCR) and sample the Superior, Penokean, Yavapai and Mazatzal tectonic provinces. The data set can be explained by a single layer of anisotropic fabric, which we interpret to be dominated by a lithospheric contribution. The fast S polarization directions are consistently ENE-WSW, but the split time varies greatly across the study area, showing strong anisotropy (up to 1.48 s) in the western Superior, moderate anisotropy in the eastern Superior, and moderate to low anisotropy in the terranes south of Lake Superior. We locate two localized zones of very low split time (less than 0.6 s) adjacent to the MCR: one in the Nipigon Embayment, an MCR-related magmatic feature immediately north of Lake Superior, and the other adjacent to the eastern end of the lake, at the southern end of the Kapuskasing Structural Zone (KSZ). Both low-splitting zones are adjacent to sharp bends in the MCR axis. We interpret these two zones, along with a low-velocity linear feature imaged by a previous tomographic study beneath Minnesota and the Dakotas, as failed lithospheric branches of the MCR. Given that all three of these branches failed to propagate into the Superior Province lithosphere, we propose that the sharp bend of the MCR through Lake Superior is a consequence of the high mechanical strength of the Superior lithosphere ca. 1.1 Ga.
The Midcontinent Rift has characteristics of a large igneous province, causing geologists to rethink some long-standing assumptions about how this giant feature formed.
Key Points A series of triaxial compression tests were performed on methane hydrate-bearing sediments and hydrate-free sediments with different amounts of fines content and at three levels of density The rise in fines content and decrease in void ratio enhanced the peak shear strength and promoted dilation behavior of methane hydrate-bearing sediments The presence of methane hydrate increased the stress ratios at the critical state of methane hydrate-bearing sediments with various amounts of fines content
