Plate Boundary Research Articles

Rapid shear zone weakening during subduction initiation

Subduction zones play a pivotal role in the mechanics of plate tectonics by providing the driving force through slab pull and weak megathrusts that facilitate the relative motion between tectonic plates. The initiation of subduction zones is intricately linked to the accumulation of slab pull and development of weakness at plate boundaries and, by consequence, the largest changes in the energetics of mantle convection. However, the transient nature of subduction initiation accompanied by intense subsequent tectonic activity, leaves critical evidence poorly preserved and making subduction initiation difficult to constrain. We overcome these limitations through a comprehensive analysis focused on Puysegur, a well-constrained extant example of subduction initiation offshore South Island, New Zealand. Through time-dependent, three-dimensional thermo-mechanical computations and quantitative comparison to new geophysical and geological observations, including topography, stratigraphy, and seismicity, we demonstrate that subduction initiation develops with a fast strain weakening described with a small characteristic displacement ([Formula: see text] 4 to 8 km). Potential physical mechanisms contributing to the strain weakening are explored and we find that the observed fast weakening may arise through a combination of grain-size reduction within the lower lithosphere and fluid pressurization at shallower depths. With the shared commonality in the underlying physics of tectonic processes, the rapid strain weakening constrained at Puysegur offers insights into the formation of the first subduction during early Earth and the onset of plate tectonics.

Spreading ridge migration enabled by plume-ridge de-anchoring

It has long been recognised that spreading ridges are kept in place by competing subduction forces that drive plate motions. Asymmetric strain rates pull spreading ridges in the direction of the strongest slab pull force, which partially explains why spreading ridges can migrate vast distances. However, the interaction between mantle plumes and spreading ridges plays a relatively unknown role on the evolution of plate boundaries. Using a numerical model of mantle convection, we show that plumes with high buoyancy flux (>3000 kg/s) can capture spreading ridges within a 1000 km radius and anchor them in place. Exceptionally high buoyancy fluxes may fragment the overriding plate into smaller plates to accommodate more efficient plate motion. If the plume buoyancy flux wanes below 1000 kg/s the ridge may be de-anchored, leading to rapid ridge migration rates when combined with asymmetric plate boundary forces. Our results show that plume-ridge de-anchoring may have contributed to the rapid migration of the SE Indian Ridge from 43 million years ago (Ma) due to waning buoyancy flux from the Kerguelen plume, supported by magma flux estimates and radiogenic isotope geochemistry of eruption products. The plume-ridge de-anchoring mechanism we have identified has global implications for the evolution of plate boundaries near mantle plumes.

Source Modeling of Deep Plate Boundary 2021 Miyagi Earthquake (Mw7.0) Employing Modified Semi-Empirical Technique with Site Effects: A Step Forward Towards Hazard Mitigation

ABSTRACT The source modeling of 2021 Miyagi earthquake (Mw7.0) offers a significant tool to understand the initial evaluation of M9 earthquake cycle in this region. This article seeks to simulate the 2021 Miyagi earthquake (Mw7.0) using the modified semi-empirical technique (MSET) after incorporating site effects determined using the Horizontal to Vertical spectral ratio technique. We propose the best-fitting source model of this earthquake from a spectrum of rupture model parameters using MSET. We believe that this effort is the first to use MSET to model this earthquake and will provide significant contribution for seismic hazard assessment of the Miyagi region.

Opening of the North Atlantic Ocean and the rise of Scandinavian mountains

Large-scale topography is usually associated with tectonic plate boundaries, but the Scandinavian mountains (Scandes) are located far from any active tectonic setting, and their origin is unknown. We demonstrate that the Precambrian lithospheric structure of Fennoscandia controlled both the Cenozoic ocean opening and mountain rise in the North Atlantic region. Our new seismic receiver function analysis reveals a block of thick continental crust formed by Proterozoic crustal stacking. We propose that this block created a wide continental shelf bounded by two transform fault systems during continental breakup. This geometry resulted in the formation of lithospheric steps at the continental margin on both sides of the stacked crustal structure, coinciding with the highest Southern and Northern Scandes. We propose that edge-driven convection at these steps caused the mountain rise. This study presents a general model for the formation of high elevation behind passive margins.

Precision Seismicity of the San Andreas Fault System in Central California: Insights into Active Tectonics and Earthquake Physics

Since the late 1960s, seismicity of the San Andreas Fault system in central California has been under intense observation by U. S. Geological Survey’s Northern California Seismic Network (NCSN). Designed to provide accurate focal depth control and a low magnitude of catalog completeness, the network maps fault, both large and small across the full breadth of this 400-kilometer-long portion of the plate boundary, including many that are known only from their seismicity. Beginning in 1984, with the advent of routine digital recording, new opportunities were created to significantly enhance the spatial resolution of seismicity. This led to the co-evolution of method that take full advantage of digital waveforms and new insights into the nature of seismicity and faults. Collectively referred to as precision seismicity, discoveries included repeating earthquakes and highly organized and stable distributions of seismicity on major faults. Long term observation of seismicity before and after the six largest earthquakes over the last half-century underscore the outstanding challenges for earthquake predictability and theories of earthquake nucleation.

Paleotectonic setting of basaltoid volcano-plutonic belts with porphyry copper deposits

The article considers geotectonic position and geological features of basaltoid volcano-plutonic belts (VPBs) associated with porphyry copper and gold-copper deposits. Three types of the belts, the oceanic, perioceanic, and riftogenic ones, are distinguished among these VPBs. The oceanic and peri-oceanic belts were incorporated into the composition of ensimatic island arcs and formed at a late stage of their development, after inversion of the basaltoid troughs. The riftogenic belts appeared in island-arc oceanic settings at a late stage of filling the spreading zones. VPBs of the first two types with porphyry gold-copper deposits are located over the convergent boundaries of tectonic plates in an above-subduction position, while belts of the third type form extended linear zones within individual rift structures and contain smallscale porphyry copper objects with poor gold and molybdenum contents. As examples of the structures, the Sunda, Philippine archipelago, New British, Solomon, Alaska, and Antilles ensimatic island arcs and the Irendyk and Novoalekseevsky belts of rift zones of the Southern Urals are described. Volcanic-plutonic associations participating in the structure of these VPBs (including those productive for porphyry copper mineralization), as well as structural-lithological complexes of their basement, are characterized. It is noted that repeated magmatism manifestations during the long-term development of the ensimatic island arcs lead to appearance of several basaltoid VPBs with porphyry gold-copper deposits, while basaltoid belts with porphyry copper objects in the rift settings are formed only once.

Crustal Structure of Southern California from Velocity and Attenuation Tomography

Recent studies of Southern California have employed velocity and attenuation tomography to elucidate two aspects of crustal structure. Low-frequency velocity tomography reveals large-scale heterogeneity that can be approximated by a discrete set (𝐾 ≤ 10) of tectonic regions, each characterized by an isostatically balanced lithospheric column reflecting the composition and tectonic history of the region. The boundaries between tectonic regions are typically localized structures expressed at the surface by major faults, topographic fronts, and geochemical transitions. The efficacy with which tomographic regionalization (TR) represents the subsurface geography of major tectonic units in Southern California indicates that the TR idealization works surprisingly well in this part of the Pacific‐North America plate boundary. High‐frequency attenuation tomography in Southern California supports the hypothesis that the decay of P and S waves propagating through the upper and middle crust is dominated by elastic scattering from small-scale heterogeneities of high fractal dimension, rather than by anelastic dissipation. The amplitude of the scattering varies systematically among the tectonic regions and correlates with the degree of recent faulting and deformation within a region. These results motivate speculation that small-scale crustal heterogeneities responsible for high-frequency attenuation are generated within the seismogenic zone through a dynamic interplay between distributed deformation and localized fracturing.

Precise Relocation and Source Characterization of Two Earthquake Sequences in Yellow Sea Using Regional Lg-Wave Observations

ABSTRACT Two significant earthquakes of magnitude ML 4.6 and 4.5 occurred on 18 January and 3 December 2021 in the central region of the Yellow Sea, respectively. The earthquakes occurred beneath the Gunsan sedimentary basin at about 10 km depth with a strike-slip faulting mechanism on nodal planes striking northwest–southeast (NW–SE) and north-northeast–south-southwest (NNE–SSW). Despite a lack of close-by seismographic stations, we successfully utilized regional Lg-wave observations on both coasts of the sea—the Korean peninsula on the east and eastern China on the west. For nine earthquakes in two event sequences, the Lg-wave differential travel times of the nearby event pairs at the common station are carefully measured using the waveform cross-correlation technique. The double-difference earthquake relocation method is employed to obtain precise relative epicentral locations using the Lg correlation measurements. Relocated epicenters align along the NW–SE direction, indicating that the nodal plane striking the same direction is the likely fault plane on which both sequences occurred. This is the first case reported in the literature in which the causative fault plane has been identified for earthquakes in the central Yellow Sea region. It has an important implication for current regional tectonics; it favors neither old tectonic features trending NE–SW (Qianliyan uplift) nor the north–south alignment of significant earthquakes in the region along the Amur plate boundary. The Lg waves from the earthquake sequences are dominant seismic signals on all three-component records at stations in 160–550 km and allowed us to analyze source properties of the two largest earthquakes using the empirical Green’s function approach. Azimuthal variations of the source corner frequencies suggest that earthquake rupture likely propagated toward southeast (125°) along the fault plane, supporting the aftershock relocation results.

Rupture direction of paleoearthquakes on the Alpine Fault, New Zealand, as recorded by curved slickenlines

We observed and further exhumed curved slickenlines on fault planes associated with paleo-surface rupture of the Alpine Fault, New Zealand’s ∼30 mm/yr continental transform plate boundary. Dynamic rupture modeling indicates that the geometry of such curvature provides a record of past earthquake rupture directions. We focused our efforts on three sites that span a region known to variably halt or allow passage of past earthquakes (an “earthquake gate”) to contribute rupture direction constraints to the fault’s spatiotemporally rich paleoseismic record. At Hokuri Creek and Martyr River, we observed both convex-up and convex-down curved slickenlines on and adjacent to principal slip surfaces, indicating past ruptures from both the northeast and southwest of these locations. At Martyr River, relationships suggest that the most recent event (inferred to correlate to 1717 CE) ruptured from the southwest. Our results demonstrate the utility of curved slickenlines as a valuable new paleoseismological tool for determining past rupture directions, applicable to surface-rupturing faults globally.

Breakdown of laminar regime in flat plate boundary layer seen in numerical solutions of the Navier–Stokes equations obtained with 16th-order scheme

Numerical solutions of the Navier–Stokes equations describing the onset and the development of the boundary layer instability of a subsonic flow about a finite width flat plate are described. They were obtained via exciters-free calculations with the 16th-order multioperators-based scheme optimized to resolve very small solutions scales. Some relevant details of the scheme are briefly outlined. The occurrence of the Tollmien–Schlichting waves near the leading edge, their downstream propagation and breakdowns some distance away are displayed. The role of the oblique waves seen in the solutions is emphasized. The nonlinear stage in the form of structured vortices followed by the solutions randomization are described.

Tectonic Inversion of Intracontinental Basins and the Conundrum of Plate‐Scale Stress Propagation: The Case of the Jatobá Basin, NE Brazil

AbstractWe chose the Jatobá Basin (JB) in the Borborema Province in northeastern Brazil to address two outstanding problems, basin tectonic inversion and stress propagation, because (a) evidence of tectonic inversion is conspicuous in the JB, and (b) the basin sits on intracontinental northeastern Brazil, far from plate boundaries and the most obvious stress sources. We found robust evidence of oblique reverse faulting inside the basin, putting sediments of the JB's depocenter hundreds of meters above the hard granitic host basement. We also found sound evidence of hydrographic flow reversal, that is, earlier basement‐to‐basin runoff, followed by basin‐to‐basement runoff. This change means that formerly the basement stood higher than the basin, and later the basin was uplifted (inverted) to currently stand hundreds of meters higher than the basement, especially where the gravity anomaly map shows the deepest depocenter. Coarse conglomerates, sourced in the uplifted basin and blanketing the adjacent basement for kilometers, are only cut by the current hydrographic network. From these main lines of evidence, we infer recent basin tectonic inversion (Quaternary, i.e., <2.6 Ma?). The WSW‐ENE orientation of the maximum compressive stress needed to invert the basin likely results from Mid‐Atlantic and Andean pushes. We conclude that lithospheric plates can deform in their cores because of stresses propagating from plate boundaries to distant plate interiors, where the lithosphere has been locally weakened and strain has concentrated.

The Curie Surface and Lithospheric Thermal Structure in Mongolia‐Baikal Region

AbstractThe Mongolia‐Baikal region, located between the Siberia, North China, and Tarim cratons, exhibits the most pronounced lithospheric accretion and extension since the Phanerozoic. Despite its distance from plate boundaries, the region's notable tectonic activity underscores its importance for studying intraplate lithospheric deformation. However, its lithospheric thermal structure, particularly at high resolution, remains inadequately characterized. By integrating magnetic anomaly constraints with heat flow measurements, we have obtained a high‐resolution lithospheric thermal structure for the Mongolia‐Baikal region, revealing lithospheric thicknesses ranging from 60 to 180 km. The Curie depth and mantle heat flow also show diverse variations across different regions. Our results suggest that elevated mantle heat flow in the Baikal Rift Zone, coupled with a uniform Curie depth and minimal lithospheric thinning, likely reflects lithospheric weakening from mantle upwelling, facilitating rifting in the early stages of the rift's formation. Additionally, some localized lithospheric thermal variations identified in our results may be closely linked to the formation of the Hangai Plateau and the low magnetic anomalies observed in the Mongolia‐Okhotsk suture.

The Law of the Total Pressure Loss in Supersonic Streamwise Corner Flow

Abstract The total pressure loss is a critical issue in the design of internal flow channels for hypersonic vehicles. In the flow through the streamwise corner region, the adjacent walls induce secondary flows that exacerbate the viscous effects, thereby intensifying the total pressure loss. This study investigates the impact of the corner’s dihedral angle and the compressibility of the flow on the total pressure loss in the streamwise corner region. The results indicated that as the dihedral angle increases, the distribution of the total pressure contour in the corner region became significantly distorted. The proportion of the boundary layer area within the corner region concurrently increased as the dihedral angle deviates from 90 degrees. Furthermore, we deduce that the total pressure loss is the largest in the supersonic flat plate boundary layer. This paper normalized the total pressure loss coefficient for different dihedral angles in the corner region and established an empirical relationship for the distribution of normalized total pressure loss coefficients along the flow path. It is discovered that under supersonic inflow conditions, the distribution law of the total pressure loss in the corner region exhibits higher predictive accuracy.

When did the Dead Sea fault become a transform?

This study re-evaluates the ~20 Myr development of the Dead Sea Fault System (DSFS) and its tectonic definition as a transform plate boundary. The DSFS conveys sinistral displacement between the Arabian-Sinai plates: ~105 km along its ~400 km-long southern segment (Gulf of Aqaba-Eilat to the Hula basin); ~90 km and 4–16 km along the central and northern segments (~190 km long each, across Lebanon, western Syria, and southern Turkey). A review of previous studies, combined with new seismological data analysis, associates the northward displacement decline with obstacles along the DSFS propagation path. During the Miocene, DSFS propagated up to the NW-trending Irbid rift (1st obstacle) and splayed NW towards the Mediterranean and NE along the Late Cretaceous Palmyra fold-thrust belt (2nd obstacle). Its reactivation uplifted the Hermon and the Anti-Lebanon mountain ranges. Northward DSFS propagation into the cold and rigid Aleppo plateau lithosphere (3rd obstacle) was stalled until the early Pliocene (~5 Ma), when volcanism and ongoing regional tectonic forcing enabled the DSFS to shift to the Yammouneh fault and rupture through the Missyaf-Ghab branch farther north (central and northern segments, respectively). During the Pleistocene-recent, connection of the DSFS with the ophiolite belt and East Anatolian Fault System (EAFS) along the Bitlis suture zone (4th obstacle) has not yet been established. Seismological data show a clear separation between the EAFS and the DSFS, while seismicity is scattered across the Aleppo plateau and the central and northern DSFS segments. In contrast, seismicity is localized along the southern DSFS segment. Our findings suggest that, at present, the DSFS has still not made a structural, seismologic, and tectonic connection with the EAFS. Hence, we redefine the DSFS as a pre-transform and suggest its interaction with the EAFS is a world-class example of a fault-fault-fault triple junction in the making.

Limit Cycle Oscillation of Elastic Plate in Supersonic Flow

Fluid–structure interaction of an elastic plate with piezoelectric elements, turbulent freestream flow at Mach 2.5, and a pressurized cavity is investigated computationally and correlated with a recent experiment. The pressure field on the plate is measured using pressure-sensitive paint, and the structural response is observed using the measured voltage of a piezoelectric patch. The measurements show a dominant frequency of oscillation which indicates the likely onset of flutter and a postflutter limit cycle oscillation (LCO). A computational investigation is conducted to study the effects of static pressure differential, temperature differential, cavity pressure coupling, and plate boundary conditions on the linear flutter onset condition and the nonlinear postflutter LCO characteristics. Rivets that connect the plate to the supporting structure are modeled as local constraint in the in-plane direction, and their effect is investigated. Computations show that the coupling between the cavity acoustic and plate structural modes is necessary for flutter onset in the wind-tunnel conditions. Correlation between computed and measured aerodynamic pressure shows reasonable agreement in amplitude and frequency. Computational results are obtained using Piston Theory and also potential flow aerodynamics. Computed and measured pressure LCO mode shapes are extracted and correlated.

Experimental investigation of seeping gas film cooling effect at downstream wall in hypersonic flow

Generating a seeping gas film through porous material to cover the downstream wall can achieve a large-area thermal protection from aerodynamic heating for the windward surfaces of hypersonic vehicles. This study explores the cooling efficiency of the seeping gas film downstream of the porous wall through experimental research. It introduces a continuous high-temperature wind tunnel heat flux calculation method and investigates the influence of cooling gas injection rate and blowing length on downstream wall cooling efficiency. Additionally, an engineering evaluation of the thermal protection effect of the seeping gas film on the downstream wall of the flat plate boundary layer is conducted. Experimental findings demonstrate that cooling efficiency decreases rapidly along the flow direction, and an increase in the cooling gas injection rate enhances the cooling efficiency downstream of the porous wall, while the attenuation rate of cooling efficiency in the downstream region accelerates with higher injection rates. Moreover, a rise in injection rate results in a reduction of the effectiveness-cost ratio in the upstream region, with the seeping gas film’s effectiveness-cost ratio being relatively high at an injection rate of 0.2 %. Longer blowing lengths enhance the cooling efficiency at the end of the porous wall but accelerate the attenuation rate of cooling efficiency along the flow direction. Conversely, under constant blowing length, a higher injection rate contributes to a faster attenuation rate of cooling efficiency downstream. Assessment of the thermal protection effect for different combinations of blowing lengths (20 ∼ 60 mm) and injection rates (0.1 %∼0.5 %) reveals that layouts featuring shorter blowing lengths and higher injection rates demonstrate superior cooling effects downstream of the porous wall when supplied with the same mass flow rate.

Compact Eddy Structure of Turbulence Aft of an Inclined Slender Body

The compact eddy structure and the correlation structure of three distinct turbulent boundary layers were investigated, with a view to understanding the turbulent inflow seen by a stern-mounted propulsion tail rotor. The first is of a zero-pressure-gradient, equilibrium, flat plate boundary layer, and the second and third are on the lee and port sides of a body of revolution (BOR) inclined at a 5° angle of attack. The BOR is an axially symmetric body with an ellipsoidal nose, straight midbody, and 1.17-diam-long stern ramp with a stern-mounted tail rotor. Measurements on the BOR were taken at a body-diameter-based Reynolds number of 600,000 and a rotor advance ratio of J=1.27. Particle image velocimetry was taken at two principal locations at the stern of the BOR. Distortion of the correlation function was evident in the compact eddy structure of the BOR boundary layer on both the lee and port sides. On the port side, the turbulent modes were dominated by the embedded shear layer, whereas the lee side showed more significant pressure gradient effects in the eddy structure.

Optimizing steel plate shear wall performance: influence of infill plate thickness and connection length to vertical boundary elements

The Steel Plate Shear Wall (SPSW) is widely used as a lateral force resisting system in areas prone to high seismic activity. This is because of its exceptional ductility and ability to dissipate energy. Previous studies have indicated that disconnecting the infill plate from the Vertical Boundary Element (VBE) can decrease the demands on the column. The observed outcomes lead to a decline in both the energy dissipation capability and the lateral load-bearing capacity. The primary goal of this investigation is to mitigate stress on the column while concurrently reducing both the energy dissipation and lateral force capacities. To accomplish this, a partial connection will be established between the infill plate and the Vertical Boundary Element. Thus, the overarching objective of this research is to identify the optimal connection length that minimizes the demands on the column while simultaneously curtailing the reduction in energy dissipation and lateral force capacities. Using Abaqus software, twenty models with different infill plate thicknesses and lengths of connection between the infill plate and Vertical Boundary Element were analysed. To determine the optimal length of connection between the infill plate and Vertical Boundary Element, the energy dissipation capacity, lateral load capacity, and column stresses were evaluated and compared. The maximum dissipation of energy occurs at a connection length of 75% between the infill plate and Vertical Boundary Element and an infill plate thickness of 12 mm. In contrast, the greatest lateral load capacity is achieved at connection lengths of 50% between a 25 mm infill plate thickness and a Vertical Boundary Element weight of 25%. The column stress reaches its minimum value when the infill plate measures 20 and 25 mm, and it reaches its maximum value when the infill plate measures 6 mm and 12 mm. By increasing the thickness of the infill plate and minimizing its connection to the Vertical Boundary Elements, the column stress in Vertical Boundary Elements can be decreased, while simultaneously enhancing lateral load capacity and maintaining an adequate level of energy dissipation.

Not trench-parallel but trench-normal source fault of the 1994 Hokkaido Toho-oki earthquake as revealed by the aftershock relocation using HypoDD

In the southwestern area of the Kurile Islands subduction zone, three large earthquakes including the Mw8.2 1969 and the Mw8.3 1994 Hokkaido Toho-oki earthquakes and the Mw7.8 1975 tsunami earthquake occurred in almost the same location. The three main shocks and their aftershocks were re-determined simultaneously using a double-difference earthquake location method, HypoDD. As a result, the 1969 and the 1975 main shocks and their aftershocks were concentrated near the upper surface of the subducting Pacific plate, indicating that both events were plate boundary megathrust earthquakes as shown by previous studies. On the other hand, the 1994 main shock and its aftershocks were distributed within the Pacific plate indicating that the 1994 event was an intraslab event as reported by previous studies. The result that interplate earthquakes and intraslab earthquakes could be appropriately distinguished is strong evidence that the accuracy of hypocenter determination is sufficiently high. The centroid moment tensor solution of the 1994 main shock has two nodal planes which are parallel or normal to the Kurile Trench. According to the high-resolution hypocenter distribution obtained in the present study, the aftershocks appear to locate on the nodal plane normal to the trench, suggesting that the source fault ruptured by the 1994 main shock was not parallel but normal to the Kurile Trench. Moreover, we found that the Coulomb failure stress change (ΔCFS) produced by the trench-normal fault plane is more likely to trigger the largest aftershock than ΔCFS produced by the trench-parallel fault plane.Graphical abstract

Miocene Alkaline Basaltic Magmatism in Northeastern Tibetan Plateau: Implications for Mantle Evolution and Plateau Outward Growth

AbstractThe widespread Cenozoic alkaline magmatism within and around the Tibetan Plateau offers a prime opportunity to probe the nature of the mantle at the depths where basalt magmas originate. The close temporal and spatial relationship between volcanism and regional strike‐slip fault systems also helps better understand the geodynamics of outward growth of the plateau in response to the continued India‐Asia convergence. We present a comprehensive study of the deeply sourced alkaline basalts formed along the Kunlun strike‐slip fault with the aim of understanding their petrogenesis and the composition of mantle sources beneath the northeastern Tibetan Plateau. High Nb/U and Ce/Pb ratios and relatively depleted bulk‐rock Sr‐Nd‐Pb isotope compositions corroborate the mantle origin of these alkaline basalts. Homogeneous and low 87Sr/86Sr of clinopyroxene indicates negligible crustal contamination during magmatic evolution. Low δ26Mg in the alkaline basalts and positive correlations with Hf/Sm and Ti/Ti* indicate that the basalts were derived from mantle that was metasomatized by melts derived from sedimentary carbonates during the Paleo‐Tethyan seafloor subduction. Based on 40Ar/39Ar dating results, it appears that the alkaline basaltic magmatism in the northeastern Tibetan Plateau occurred simultaneously with Kunlun strike‐slip faulting. These observations suggest that the India‐Asia convergence must have reactivated ancient subduction plate boundaries and resulted in strike‐slip faulting along these suture zones within and around the Tibetan Plateau. The eruption of low‐volume and deeply rooted alkaline basalts may have been controlled by fractures associated with the strike‐slip fault systems.