Search

Paleomagnetic data constrain paleogeographic motion of rocks relative to the Earth's spin axis, which is a sum of plate motion relative to the mantle and true polar wander. Discerning between these effects is challenging for studies aiming to reconstruct paleo-plate motions from deformed orogenic terranes. Here, we study the paleolatitudinal drift history of the Lhasa terrane of southern Tibet that migrated from the northern Gondwana to the southern Eurasian margin between late Triassic and early Cretaceous time. Previous work identified a 180 Ma near-equatorial Lhasa latitude and assumed near-constant paleolatitudinal drift. Large-scale true polar wander at this time, however, which has been argued for in previous work, requires highly variable Lhasa plate motion rates relative to Gondwana. Here, we test whether the alternative interpretation of constant plate motion rates provides a better prediction of paleomagnetic data. To this end, we present a new paleomagnetic pole from ∼155 Ma volcanics (here dated by U/Pb zircon) of the Lhasa terrane. Our pole comprises site mean directions from 46 lavas, passes a fold test, and is supported by an extensive rock magnetic and microscopic analysis that reveals no evidence of remagnetization. Our results give an average direction of D±ΔDx=337.1°±2.6°, I±ΔIx=−13.8°±4.9°, and a corresponding paleopole position at λp=45.3°N, φp=295.3°E with K=69.8 and A95=2.5°, predicting a similar near-equatorial, paleolatitude as at 180 Ma. This paleolatitudinal standstill is consistent with predicted paleolatitudes from global APWPs that account for Jurassic TPW and assume a constant Lhasa-Gondwana (India) oceanic spreading rate of ∼8 cm/a during the 215-130 Ma opening of the Neotethys and closure of the Mesotethys oceans. Contemporaneous arc magmatism on the Lhasa terrane was previously interpreted to indicate that the Lhasa terrane was located above a subduction zone during its northward journey to Tibet: our results show that if our reconstruction is accurate, then trench migration rates also must have been ∼8 cm/a. Our results are consistent with the range of rates of global long-lived trench migration that have been reconstructed in other studies, but they are much higher than measured in most modern systems. Our results more completely document the Lhasa terrane's separation from Gondwana and provide a straightforward methodology of using paleomagnetic data to test plate kinematic scenarios, overcoming the problem of the unknown effect of true polar wander on paleomagnetic data from deformed terranes.

The high-pressure metamorphic Nevado-Filábride Complex (NFC) in the Betics mountain range of southeastern Spain exhibits continental and ocean-derived tectonic units, which are key for understanding the geodynamic evolution of the Western Mediterranean. We address the current debate in the definition of tectonic units, the emplacement of (ultra)mafic rocks, and the timing of burial metamorphism by conducting a structural study combined with single grain fusion 40Ar/39Ar dating of white micas in structurally critical outcrops of the eastern Sierra de Los Filábres. One older 40Ar/39Ar age population (38–27 Ma) is found at distance from the main shear zones in the relics of an early foliation, while a younger 40Ar/39Ar population (22–12 Ma) is dominant in the vicinity of these shear zones, where the early foliation is obliterated. Both age groups are interpreted as the record of deformation or fluid-induced recrystallization during distinct fabric-forming events, while alternative scenarios are discussed. A key observation is the presence of an ophiolitic mélange, which—together with new and published geochronological data—allows for a new tectonic hypothesis. This considers Paleogene subduction beneath a Jurassic oceanic lithosphere, followed by the continued subduction of NFC and overlying ophiolites below the Alpujárride Complex. Exhumation during westward slab roll-back led to the formation of an extensional detachment system that obliquely cut nappe contacts. Although the timing constraints for high pressure-low temperature (HP-LT) metamorphism in the NFC remain inconclusive, the new tectonic hypothesis provides a solution that can account for both Paleogene and Miocene ages of HP-LT metamorphism.

We present a model for the lattice vibrational density of states of MgTiO3 (geikielite) and FeTiO3 (ilmenite) that predicts thermodynamic properties, in agreement with observational data. The model is based on Kieffer’s method combined with spectroscopic data. For both substances experimental data sets are influenced by non-stoichiometry. For geikielite that affects the volume, whereas for ilmenite volume and bulk modulus are affected. We show that Kieffer’s method enables predicting bulk moduli in pressure–temperature space. We demonstrate that intrinsic anharmonicity or electronic effects significantly affect the heat capacity of ilmenite, whereas that is not the case for geikielite. We use Kieffer’s method to derive multiple-Einstein models, from which we demonstrate that thermodynamic properties are insignificantly influenced by dispersion in Grüneisen, mode-q and anharmonicity parameters for both substances. We show that our results enable predicting thermodynamic properties and shear modulus of the solid solution formed from geikielite and ilmenite. Geikielite and ilmenite are added to our thermodynamic database for the system MgO–SiO2–FeO, to enable modeling phase stability and physical properties of titanium-rich reservoirs in the Earth’s Moon.

A hitherto unknown Neoproterozoic orogenic system, the Saharides, is described in North Africa. It formed during the 900–500-Ma interval. The Saharides involved large subduction accretion complexes occupying almost the entire Arabian Shield and much of Egypt and parts of the small Precambrian inliers in the Sahara including the Ahaggar mountains. These complexes consist of, at least by half, juvenile material forming some 5 million km(2) new continental crust. Contrary to conventional wisdom in the areas they occupy, evolution of the Saharides involved no continental collisions until the end of their development. They formed by subduction and strike-slip stacking of arc material mostly by precollisional coastwise transport of arc fragments rifted from the Congo/Tanzania cratonic nucleus in a manner very similar to the development of the Nipponides in east Asia, parts of the North American Cordillera and the Altaids. The Sahara appears to be underlain by a double orocline similar to the Hercynian double orocline in western Europe and northwestern Africa and not by an hypothetical “Saharan Metacraton.” The method we develop here may be useful to reconstruct the structure of some of the Precambrian orogenic belts before biostratigraphy became possible.

Subduction zones are pivotal for the recycling of Earth’s outer layer into its interior. However, the conditions under which new subduction zones initiate are enigmatic. Here, we constructed a transdisciplinary database featuring detailed analysis of more than a dozen documented subduction zone initiation events from the last hundred million years. Our initial findings reveal that horizontally forced subduction zone initiation is dominant over the last 100 Ma, and that most initiation events are proximal to pre-existing subduction zones. The SZI Database is expandable to facilitate access to the most current understanding of subduction zone initiation as research progresses, providing a community platform that establishes a common language to sharpen discussion across the Earth Science community., Despite numerous advances in our understanding of subduction since the theory of plate tectonics was established, the mechanisms of subduction zone initiation remain highly controversial. Here, the authors present a transdisciplinary and expandable community database of subduction zone initiation events in the last 100 Ma, which establishes a clear direction for future research.

To determine a paleopole, the paleomagnetic community commonly applies a loosely defined set of quantitative data filters that were established for studies of geomagnetic field behavior. These filters require costly and time-consuming sampling procedures, but whether they improve the precision and influence the position of paleopoles has not yet been systematically analyzed. In this study, we performed a series of experiments on four datasets which consist of 73–125 lava sites with 6–7 samples per lava. The datasets are from different regions and ages, and are large enough to represent paleosecular variation, yet include demonstrably unreliable paleomagnetic directions. We show that the systematic application of data filters based on within-site scatter (a maximum angular deviation filter on individual directions, a k-cutoff, a minimum number of samples per site, and eliminating the farthest outliers per site) cannot identify unreliable directions. We find instead that excluding unreliable directions relies on the subjective interpretation of the expert, highlighting the importance of making all data available following the FAIR principles. In addition, data filters that decrease the number of sites even have an adverse effect; they decrease the precision of the paleopole. Between-site scatter often outweighs within-site scatter, and when collecting paleomagnetic poles, the extra efforts put into collecting multiple samples per site are more effectively spent on collecting more single-sample sites.

Knowledge of the thermal history of the Linzizong Group (69–47 Ma) within the Gangdese arc is critical for interpreting the geologic evolution and isotope-based paleoaltimetric results of the southern Lhasa terrane of the Tibetan Plateau. Here, we combine results from geochronologic and thermochronologic studies of this group (divided into Dianzhong, Nianbo, and Pana formations upsection) and the structurally overlying Qianggeren granite (∼52 Ma) in the Linzhou basin. Whole rock 40Ar/39Ar ages of volcanic rocks from the stratigraphically lower Dianzhong and Nianbo formations are ∼10 Myr younger than their corresponding zircon U-Pb ages, suggesting a thermal disturbance of the argon system. Zircon (U-Th)/He ages (ZHe, 63 dates) range from 54 to 24 Ma, and apatite (U-Th-Sm)/He ages (AHe, 43 dates) range from 27 to 4 Ma. Inverse modeling of the thermochronologic data from the Qianggeren granite indicates rapid cooling between 42 and 26 Ma, possibly induced by movement of the Gulu-Hamu thrust. Positive correlations between ZHe ages and effective uranium and other geologic observations provide evidence that the Dianzhong and Nianbo formations were heated 300 °C at 54–50 Ma, and that the entire Linzizong Group was variably heated to 130–170 °C at 42–27 Ma. These findings, together with published geochronologic and thermochronologic data, suggest that abundant 50–45 Ma zircon fission track and ZHe ages from the Gangdese arc likely reflect conductive cooling of the Gangdese arc after a ∼52 Ma magmatic flare-up episode rather than rapid regional exhumation. Areas far from Cenozoic faults and deeply incised river valleys in southern Tibet have experienced only ∼3 km of exhumation since 45 Ma, consistent with the establishment of a low-relief, plateau-like physiography by Eocene time. The (hydro)thermal events are also manifested by widespread calcite recrystallization and δ18O and Δ47 alteration in most carbonates from the Linzizong Group, which must be taken into account in past and future carbonate-based paleoaltimetric studies.

Global reorganizations of tectonic plates may be caused by a trigger such as a continental collision or a rising mantle plume. However, whether and how such a trigger propagates through a plate circuit remains unclear. Here we use a plate kinematic model to quantify relative motions between the African and Eurasian plates following a plume-induced plate motion change that triggered formation of a new subduction zone within the Neotethys Ocean at 105 Ma. We constrain the plate kinematic model by geomagnetic intensity variations recorded in Atlantic quiet zone crust that formed during the Cretaceous Normal Superchron (126–83 Ma), during which magnetic reversals were absent. We find that convergence rate changes between Africa and Eurasia are well explained by the initiation and arrest of the plume-induced subduction zone. Our plate kinematic model also reveals that the plate acceleration that followed upon subduction initiation changed the Africa–Eurasia convergence direction, which in turn was accommodated by subsequent subduction initiation about 85 Ma in the Alpine region that then triggered a cascade of regional tectonic events. This plate tectonic chain reaction illustrates how changes in plate motion, underpinned by mantle dynamics, may self-perpetuate through a plate circuit.

The India–Asia collision zone is the archetype to calibrate geological responses to continent–continent collision, but hosts a paradox: there is no orogen-wide geological record of oceanic subduction after initial collision around 60–55 Ma, yet thousands of kilometers of post-collisional subduction occurred before the arrival of unsubductable continental lithosphere that currently horizontally underlies Tibet. Kinematically restoring incipient horizontal underthrusting accurately predicts geologically estimated diachronous slab break-off, unlocking the Miocene of Himalaya–Tibet as a natural laboratory for unsubductable lithosphere convergence. Additionally, three endmember paleogeographic scenarios exist with different predictions for the nature of post-collisional subducted lithosphere but each is defended and challenged based on similar data types. This paper attempts at breaking through this impasse by identifying how the three paleogeographic scenarios each challenge paradigms in geodynamics, orogenesis, magmatism or paleogeographic reconstruction and identify opportunities for methodological advances in paleomagnetism, sediment provenance analysis, and seismology to conclusively constrain Greater Indian paleogeography.

Several alternative gravity forward modelling methodologies and associated numerical codes with their own advantages and limitations are available for the Solid Earth community. With the upcoming state-of-the-art lithosphere density models and accurate global gravity field data sets it is vital to understand the opportunities and limitations of the various approaches. In this paper, we discuss the four widely used techniques: global spherical harmonics (GSH), tesseroid integration (TESS), triangle integration (TRI), and hexahedral integration (HEX). A constant density shell benchmark shows that all four codes can produce similar precise gravitational potential fields. Two additional shell tests were conducted with more complicated density structures: lateral varying density structures and a Moho density interface between crust and mantle. The differences between the four codes were all below 1.5 percent of the modeled gravity signal suitable for reproducing satellite-acquired gravity data. TESS and GSH produced the most similar potential fields (< 0.3 percent). To examine the usability of the forward modelling codes for realistic geological structures, we use the global lithosphere model WINTERC-G, that was constrained, among other data, by satellite gravity field data computed using a spectral forward modeling approach. This spectral code was benchmarked against the GSH and it was confirmed that both approaches produce similar gravity solution with negligible differences between them. In the comparison of the different WINTERC-G-based gravity solutions, again GSH and TESS performed best. Only short-wavelength noise is present between the spectral and tesseroid forward modelling approaches, likely related to the different way in which the spherical harmonic analysis of the varying boundaries of the mass layer is performed. The Spherical harmonic basis functions produces small differences compared to the tesseroid elements especially at sharp interfaces, which introduces mostly short-wavelength differences. Nevertheless, both approaches (GSH and TESS) result in accurate solutions of the potential field with reasonable computational resources. Differences below 0.5 percent are obtained, resulting in residuals of 0.076 mGal standard deviation at 250 km height. The biggest issue for TRI is the characteristic pattern in the residuals that is related to the grid layout. Increasing the resolution and filtering allows for the removal of most of this erroneous pattern, but at the expense of higher computational loads with respect to the other codes. The other spatial forward modelling scheme HEX has more difficulty in reproducing similar gravity field solutions compared to GSH and TESS. These particular approaches need to go to higher resolutions, resulting in enormous computation efforts. The hexahedron-based code performs less than optimal in the forward modelling of the gravity signature, especially of a lateral varying density interface. Care must be taken with any forward modelling software as the approximation of the geometry of the WINTERC-G model may deteriorate the gravity field solution.

Following hundreds of millions of years of subduction in all circum-Pacific margins, the Pacific Plate started to share a mid-ocean ridge connection with continental Antarctica during a Late Cretaceous south Pacific plate reorganization. This reorganization was associated with the cessation of subduction of the remnants of the Phoenix Plate along the Zealandia margin of East Gondwana, but estimates for the age of this cessation from global plate reconstructions (∼86 Ma) are significantly younger than those based on overriding plate geological records (105–100 Ma). To find where this discrepancy comes from, we first evaluate whether incorporating the latest available marine magnetic anomaly interpretations change the plate kinematic estimate for the end of convergence. We then identify ways to reconcile the outcome of the reconstruction with geological records of subduction along the Gondwana margin of New Zealand and New Caledonia. We focus on the plate kinematic evolution of the Phoenix Plate from 150 Ma onward, from its original spreading relative to the Pacific Plate, through its break-up during emplacement of the Ontong Java Nui Large Igneous Province into four plates (Manihiki, Hikurangi, Chasca, and Aluk), through to the end of their subduction below East Gondwana, to today. Our updated reconstruction is in line with previous compilations in demonstrating that as much as 800–1100 km of convergence occurred between the Pacific Plate and Zealandia after 100 Ma, which was accommodated until 90–85 Ma. Even more convergence occurred at the New Zealand sector owing to spreading of the Hikurangi Plate relative to the Pacific Plate at the Osbourn Trough, with the most recent age constraints suggesting that spreading may have continued until 79 Ma. The end of subduction below most of East Gondwana coincides with a change in relative plate motion between the Pacific Plate and East Gondwana from westerly to northerly, of which the cause remains unknown. In addition, the arrival of the Hikurangi Plateau in the subduction zone occurred independent from, and did not likely cause, the change in Pacific Plate motion. Finally, our plate reconstruction suggests that the previously identified geochemical change in the New Zealand arc around 105–100 Ma that was considered evidence of subduction cessation, may have been caused by Aluk-Hikurangi ridge subduction instead. The final stages of convergence before subduction cessation must have been accommodated by subduction without or with less accretion. This is common in oceanic subduction zones but makes dating the cessation of subduction from geological records alone challenging.

This article is composed of three independent commentaries about the state of Integrated, Coordinated, Open, Networked (ICON) principles (Goldman et al., 2021, https://doi.org/10.1002/essoar.10508554.1) in the Geomagnetism, Paleomagnetism, and Electromagnetism (GPE) section and discussion on the opportunities and challenges of adopting them. Each commentary focuses on a different topic: Global collaboration, reproducibility, data sharing and infrastructure; Inclusive equitable, and accessible science: Involvement, challenges, and support of early career, BIPOC, women, LGBTQIA+, and/or disabled researchers; Community engagement, citizen science, education, and stakeholder involvement. Data sharing practices and open repository use still varies strongly between GPE communities. Some have a long tradition of data sharing; others are only starting it. Globally, GPE leadership is strongly dominated by white males and diversity may increase through the creation of Science Equality Commissions. Improved global stakeholder involvement can increase research impacts and help fight inequalities. In all investigated topics we see promising beginnings but also recognize obstacles that include a lack of funding, a lack of understanding of diversity, and prioritizing short-term gain over long-term benefit. Nonetheless, we are hopeful that our community will embrace ICON science.

We provide an independent test for the proposed pre-breakup fits between the South American and African continental margins based on new palaeomagnetic data. In the course of modelling the opening of the South Atlantic Ocean, several reconstructions have been proposed, which correspond to different estimates of the amount of pre-drift extension in the South American and African continental margins, and amounts of displacement along intracontinental accommodation zones. The Parana–Etendeka Large Igneous Province erupted during a brief time interval between ~ 136 and ~ 131 Ma, just prior to opening of the South Atlantic, and has magmatic products on both the African (Etendeka lavas) and South American (Parana lavas) plates. Here we provide new palaeomagnetic constraints on the Etendeka volcanic rocks. These yield a declination of D ± ∆ Dx = 314.4 ± 7.0° and inclination of I ± ∆ Ix = − 42.7 ± 8.3° (Plat = 47.5, PLon = 88.9, K = 26.1, A95 = 6.3, n = 21). When our results are combined with another recently published dataset, we obtain an Etendeka palaeopole at 49.1° S, 87.6° E (D ± ∆ Dx = 315.9 ± 3.8°, I ± ∆ Ix = − 43.9 ± 4.4°, K = 23.4, A95 = 3.5, n = 75). We compare three recent models for South Atlantic Ocean opening, discuss which of the models best fits the palaeomagnetic data, and present a revised reconstruction for West Gondwana during the Early Cretaceous.

We applied a lattice vibrational technique, based on representing the vibrational density of states with multiple-Einstein frequencies, to determine consistency of data on thermophysical properties and phase diagrams in the system MgO–FeO–SiO 2 . We present analyses of these data in the temperature range between 0 and 2000 K and pressure range between 0 and 20 GPa. The result is a database containing phases relevant to the Earth upper mantle and transition zone. We show that consistency of different datasets associated with the dissociation of the ringwoodite form of Fe 2 SiO 4 depends on the crucible material that has been used to perform partitioning experiments between ringwoodite and ferropericlase, and that this results in different phase diagrams for FeSiO 3 and the post-spinel part of Mg 2 SiO 4 –Fe 2 SiO 4 . We show that the existence of a phase field coesite + ringwoodite in the phase diagram of FeSiO 3 is possible and that it might be used to fine-tune pressure scales. We demonstrate that the phase boundary between coesite and quartz is very sensitive to the low-temperature heat capacity of coesite and that heat capacity data of β-quartz are too large to be reconciled with the phase boundary between β-quartz and coesite. We compare our results with seismic data associated with the 410 km seismic discontinuity.

Paleomagnetic poles used to compute apparent polar wander paths (APWPs) are strongly dispersed, which was recently shown to cause a large fraction (>50%) of these poles to be statistically distinct from the APWP to which they contributed, suggesting that current statistical approaches overestimate paleomagnetic resolution. Here, we analyze why coeval paleopoles are so dispersed, using the paleopoles behind the most recent global APWP and a compilation of paleomagnetic data obtained from

Geodynamic modelling provides a powerful tool to investigate processes in the Earth's crust, mantle, and core that are not directly observable. However, numerical models are inherently subject to the assumptions and simplifications on which they are based. In order to use and review numerical modelling studies appropriately, one needs to be aware of the limitations of geodynamic modelling as well as its advantages. Here, we present a comprehensive yet concise overview of the geodynamic modelling process applied to the solid Earth from the choice of governing equations to numerical methods, model setup, model interpretation, and the eventual communication of the model results. We highlight best practices and discuss their implementations including code verification, model validation, internal consistency checks, and software and data management. Thus, with this perspective, we encourage high-quality modelling studies, fair external interpretation, and sensible use of published work. We provide ample examples, from lithosphere and mantle dynamics specifically, and point out synergies with related fields such as seismology, tectonophysics, geology, mineral physics, planetary science, and geodesy. We clarify and consolidate terminology across geodynamics and numerical modelling to set a standard for clear communication of modelling studies. All in all, this paper presents the basics of geodynamic modelling for first-time and experienced modellers, collaborators, and reviewers from diverse backgrounds to (re)gain a solid understanding of geodynamic modelling as a whole.

Orogens that form at convergent plate boundaries typically consist of accreted rock units that form an incomplete archive of subducted oceanic and continental lithosphere, as well as of deformed lithosphere of the former upper plate. Reading the construction of orogenic architecture forms the key to decipher the pre-orogenic paleogeographic distribution of oceans and continents, as well as bathymetric and topographic features that existed thereon such as igneous plateaus, seamounts, microcontinents, or magmatic arcs. Current classification schemes of orogens divide between settings associated with termination of subduction [continent-continent collision, continent-ocean collision (obduction)] and with ongoing subduction (accretionary orogenesis), alongside intraplate orogens. Perceived diagnostic features for such classifications, particularly of collisional orogenesis, hinge on dynamic interpretations linking downgoing plate paleogeography to upper plate deformation, plate motion changes, or magmatism. Here, we show, however, that Mesozoic-Cenozoic orogens that undergo collision almost all defy these proposed diagnostic features and behave as accretionary orogens instead. To reconstruct paleogeography of subducted and upper plates, we therefore propose an alternative approach to navigating through orogenic architecture: subducted plate units comprise nappes (or mélanges) with Ocean Plate Stratigraphy (OPS) and Continental Plate Stratigraphy (CPS) stripped from their now-subducted or otherwise underthrust lower crustal and mantle lithospheric underpinnings. Upper plate deformation and paleogeography respond to the competition between absolute motions of the upper plate and the subducting slab. Our navigation approach through orogenic architecture aims to avoid a priori dynamic interpretations that link downgoing plate paleogeography to deformation or magmatic responses in the upper plate, to provide an independent basis for geodynamic analysis. From our analysis we identify ‘rules of orogenesis' that link the rules of rigid plate tectonics with the reality of plate deformation. We use these rules for a thought experiment, in which we predict orogenic architecture that will result from subducting the present-day Indian Ocean and colliding the Somali, Madagascar, and Indian margins using a published continental drift scenario for a future supercontinent as basis. We illustrate that our inferred rules (of thumb) generate orogenic architecture that is analogous to elements of modern orogens, unlocking the well-known modern geography as inspiration for developing testable hypotheses that aid interpreting paleogeography from orogens that formed since the birth of plate tectonics.

We provide new kinematic data from the Potwar Plateau (Pakistan) to evaluate the tectonic evolution of the region during the Neogene. The plateau is bound by two major strike-slip faults in the west and the east, accommodating its southwards translation. We have recognized two Neogene deformation phases in the plateau, based on paleostress inversion and Anisotropy of Magnetic Susceptibility (AMS) tensors. The first phase lasted until the early Pliocene and was characterized by vertical minor stress and N-S compression, implying thrust tectonics. The second deformation phase is characterized by a near-vertical intermediate principal stress and near-horizontal major and minor stresses, interpreted to be associated with strike-slip tectonics since the late Pliocene. K-int vectors from 21 sites are relatively compatible with the major principal stress orientations (sigma(1)) and indicate two distinct domains. This is possibly because K-min orientations are related to compaction, whereas K-int orientations were always parallel to tectonic shortening and hence compression direction during both strike-slip (post-late Pliocene) and thrusting (pre-late Pliocene) phases. These phases are characterized by swapping of (sigma(2)) and (sigma(3)) orientations while (sigma(1)) maintained its orientation. The most prominent change occurs at the western part of the Potwar Plateau, where major principal stress directions (sigma(1)) and K-int axes fan out south-westwards. The eastern domain is dominated by NE-SW trending folds and thrust faults, which are absent in the western domain. These structural features are interpreted to be the result of the distribution of deposits of the Neoproterozoic Salt Range Formation as a substratum below the Potwar Plateau. The Salt Range Formation is very thick and widespread in the west area and almost absent in the east. This factor led to unconstrained southwards gliding of the Potwar Plateau over the salt deposits in the west as opposed to frictional sliding and substantial internal deformation in the east.