Revisiting hotspots and continental breakup—Updating the classical three-arm model

ABSTRACT Classic models proposed that continental rifting begins at hotspots—domal uplifts with associated magmatism—from which three rift arms extend. Rift arms from different hotspots link up to form new plate boundaries, along which the continent breaks up, generating a new ocean basin and leaving failed arms, termed aulacogens, within the continent. In subsequent studies, hotspots became increasingly viewed as manifestations of deeper upwellings or plumes, which were the primary cause of continental rifting. We revisited this conceptual model and found that it remains useful, though some aspects require updates based on subsequent results. First, the rift arms are often parts of boundaries of transient microplates accommodating motion between the major plates. The microplates form as continents break up, and they are ultimately incorporated into one of the major plates, leaving identifiable fossil features on land and/or offshore. Second, much of the magmatism associated with rifting is preserved either at depth, in underplated layers, or offshore. Third, many structures formed during rifting survive at the resulting passive continental margins, so study of one can yield insight into the other. Fourth, hotspots play at most a secondary role in continental breakup, because most of the associated volcanism reflects plate divergence, so three-arm junction points may not reflect localized upwelling of a deep mantle plume.

Faulting of rocks as a critical energy process

Faulting of rocks is a dominant earth process that governs small-scale fracturing, formation of tectonic plate boundaries, and earthquakes occurrence1–4. Since the 18th century, the mechanical settings for rock faulting were commonly analyzed with the Coulomb criterion5 that offers empirical, useful tools for scientific and engineering applications1,6–12. Here we revisit the processes of rock faulting by an alternative approach that incorporates elastic energy, strain-state, and three-dimensional deformation; these mechanical fundamentals are missing in Coulomb criterion. We propose that a stressed rock-body fails as two conditions are met: (1) The elastic energy generated by the loading system equals or exceeds a critical energy intensity that is required for the faulting process; (2) The internal strain of the stressed rock-body due to slip and dilation along the developing faults equals the strain-state created by the loading system to maintain physical continuity13,14. Our simulations reveal that meeting these energy and strain conditions requires an orthorhombic, polymodal fault geometry that is similar to natural and experimental fault systems15–20. The application of our formulation to hundreds of rock-mechanics experiments11,21–28 provides a new, comprehensive benchmark for rock-faulting.

When did the North Anatolian fault reach southern Marmara, Turkey?

We dated syntectonic calcites on fault planes from the southern branch of the western North Anatolian fault (NAF) in northern Turkey using U-Th geochronology. We selected strike-slip faults that are kinematically related to the current regional strain field. The isotopic ages cluster around different periods during the past ~700 k.y. The most prominent cluster peak of 510.5 ± 9.5 ka (1σ) is consistent with the maximum cumulative strike-slip offset data and tectonic plate motions measured by GPS data, highlighting the fact that the present configuration of the NAF in the southern Marmara region started at ca. 500 ka or earlier. These new isotopic ages, combined with previous considerations of regional tectonics, reveal that faulting along the western NAF initiated primarily in the southern Marmara region at least a few hundred thousand years earlier than the timing suggested for the northern branch of the western NAF. This study presents an innovative approach to constrain the timing of initiation of currently active fault segments along the NAF in southern Marmara. U-Th geochronology of fault-hosted calcite thus has a wide application in determining absolute ages of fault episodes in wider shear zones along plate boundaries.

The Borborema Strike-Slip Shear Zone System (NE Brazil): Large-Scale Intracontinental Strain Localization in a Heterogeneous Plate

Large-scale strike-slip faults are fundamental tectonic elements of the continental lithosphere. They constitute plate boundaries (continental transforms), separate terranes with contrasting geological histories within accretionary orogens, or accommodate heterogeneous deformation in intracontinental settings. In ancient orogens, where deeper levels of the crust are exposed, these faults are expressed as shear zones materialized by up to tens of km-wide mylonitic belts. The Borborema shear zone system in northeastern Brazil is one of the largest and best-exposed intracontinental strike-slip shear zone systems in the world, cropping out over 250,000 km2. Here, we review its main geophysical, structural, petrologic, and geochronologic characteristics and discuss the factors controlling its development. This complex continental scale shear zone system is composed of a set of NE- to NNE-trending dextral shear zones from which there are two major E-trending dextral shear zones with horse-tail terminations into the transpressional belt branch, as well as several smaller E-trending dextral and NE-trending dextral and sinistral shear zones. The major shear zones are marked by extensive linear or curvilinear magnetic gradients, implying their continuation at depth. The major shear zones are materialized by migmatite to amphibolite-facies mylonites, but the entire system shows evidence of late deformation at lower temperatures. The system developed during the late stages of the Neoproterozoic Brasiliano (Pan-African) orogeny (mainly from 590 to 560 Ma), postdating by more than 20 Ma the main stage of contractional deformation. Localization of strike-slip shearing in this intraplate setting was controlled by rheological contrasts between blocks with distinct Paleoproterozoic histories, the presence of preorogenic Neoproterozoic rifts, the craton geometry, and zones of enhanced magmatic activity, highlighting the importance of rheological heterogeneity in controlling shear zone nucleation and evolution.

The Horizontal Kinematics of the North Island of New Zealand

<p>The advantages and disadvantages of the 'displacement' approach and the 'strain' approach to the analysis of repeated geodetic surveys for crustal deformation are discussed and two methods of geodetic strain analysis are described in detail. Repeated geodetic surveys in the central North Island show i) secular widening of the Taupo Volcanic Zone (TVZ) at 7 mm y-1 without significant transcurrent motion ii) north-south dextral motion at 14 mm y-1 and east-west narrowing at 4 mm y-1 across the northern end of the North Island Shear Belt iii) 3.1 m extension at 135' across a 15 km-wide region north of Lake Taupo, and adjacent zones of compressive rebound all associated with the 1922 Taupo Earthquakes. From the epicentral distribution and horizontal strain pattern a 15 km-square fault dipping 40' and striking parallel to the TVZ is inferred for the 1922 earthquakes. The seismic moment, 1.3 x 10 26 dyne cm, and the stress drop, 134 bars, are abnormally high for the TVZ. Widening of the TVZ is considered to be back-arc spreading. The spreading axis is postulated to extend northeast into the Havre Trough via a north-south dextral transform; and southwest into the Waverley Fault Zone and Waimea Depression via the sinistral reverse Raetihi Transform. Deformation of the North Island is not homogeneous. Fault zones are idealized as line plate boundaries and four plates -Indian, Central, Kermadec and Pacific - are postulated to account for the deformation. The Indian-Pacific macroplate pole is adopted and non-unique positions and rotation rates for the remaining poles are determined from geodetic strain data and the geometry of plate interactions. The Central Plate is moving away from the Indian Plate at the back-arc spreading axis; the Kermadec Plate is moving dextrally with respect to the Central Plate at the North Island Shear Belt which accommodates most of the transcurrent component of motion between the Indian and Pacific plates in the North Island and gives almost pure subduction of the Pacific Plate under the Kermadec Plate at the Hikurangi Margin.</p>

The Horizontal Kinematics of the North Island of New Zealand

<p>The advantages and disadvantages of the 'displacement' approach and the 'strain' approach to the analysis of repeated geodetic surveys for crustal deformation are discussed and two methods of geodetic strain analysis are described in detail. Repeated geodetic surveys in the central North Island show i) secular widening of the Taupo Volcanic Zone (TVZ) at 7 mm y-1 without significant transcurrent motion ii) north-south dextral motion at 14 mm y-1 and east-west narrowing at 4 mm y-1 across the northern end of the North Island Shear Belt iii) 3.1 m extension at 135' across a 15 km-wide region north of Lake Taupo, and adjacent zones of compressive rebound all associated with the 1922 Taupo Earthquakes. From the epicentral distribution and horizontal strain pattern a 15 km-square fault dipping 40' and striking parallel to the TVZ is inferred for the 1922 earthquakes. The seismic moment, 1.3 x 10 26 dyne cm, and the stress drop, 134 bars, are abnormally high for the TVZ. Widening of the TVZ is considered to be back-arc spreading. The spreading axis is postulated to extend northeast into the Havre Trough via a north-south dextral transform; and southwest into the Waverley Fault Zone and Waimea Depression via the sinistral reverse Raetihi Transform. Deformation of the North Island is not homogeneous. Fault zones are idealized as line plate boundaries and four plates -Indian, Central, Kermadec and Pacific - are postulated to account for the deformation. The Indian-Pacific macroplate pole is adopted and non-unique positions and rotation rates for the remaining poles are determined from geodetic strain data and the geometry of plate interactions. The Central Plate is moving away from the Indian Plate at the back-arc spreading axis; the Kermadec Plate is moving dextrally with respect to the Central Plate at the North Island Shear Belt which accommodates most of the transcurrent component of motion between the Indian and Pacific plates in the North Island and gives almost pure subduction of the Pacific Plate under the Kermadec Plate at the Hikurangi Margin.</p>

Geophysical Studies in the Southwest Pacific : Primarily Studies of Crustal Structure between New Zealand and Antarctica

<p>Geophysical data - primarily magnetic field measurements, bathymetry, and seismicity data - are presented for the area between New Zealand and Antarctica from approximately 145[degrees]W to 155[degrees]E. The data are used to determine the structure of the Pacific-Antarctic boundary, the oceanic part of the Pacific plate and the area of intersection of the Indian, Pacific and Antarctic plates. In the southwest Pacific basin the magnetic anomalies are very clear and an extensive pattern of anomaly lineations with some offsets is mapped. The magnetic anomalies show that the uniform Pacific basin area formed between about 83 and 63 mybp. The Pacific-Antarctic boundary is shown to differ either side of about 175[degrees]W. To the east it is a relatively uniform aseismic spreading ridge, offset by some transform faults. West of 175[degrees]W, to 161[degrees]E, the boundary consists of a seismically active zone of disturbed bathymetry and magnetic anomalies striking about N.70[degrees]W. The zone, the Pacific-Antarctic fracture zone, probably consists of several fractures striking about N45[degrees]W. The area between the Pacific-Antarctic boundary and the southwest Pacific basin represents the interval 10 to -55 mybp, and only in the east are anomaly lineations clear. The Indian-Antarctic Pacific triple junction is near 61.5[degrees]S, 161[degrees]E and is a stable ridge-fault-fault junction; the Indian-Antarctic boundary being the ridge. Plate tectonics is applied to the area and the structure is shown to fit, and be explained by a different rotation pole for each of the major intervals indicated by the structure, i.e. 0-10 mybp, 10-63 mybp and 63-80 mybp. The poles, with rotation rates deduced from the magnetic anomalies, are used to reconstruct the position of New Zealand relative to Antarctica at 80 mybp. The two continents probably started to separate at close to 83 mybp. The times of the major changes of structure and plate movement in the area are shown to coincide with major plate movement changes in the southwest Pacific area and in the rest of the world. A new method for determining poles of rotation, based only on epicentre locations is presented, The method is applied to independently determine the Indian-Pacific, Pacific-Antarctic and Indian-Antarctic poles. The poles should form a consistent. set and they do. The method yields effectively instantaneous poles, is quantitative, and is applicable to most plate boundaries. Earthquake magnitude-frequency relationship b values for the plate boundaries in the area are determined. Comparisons with results from elsewhere indicate an association of high b with high temperature and conversely. Several factors which have previously been suggested as determining b value are shown to not be determinants. A revised and extended magnetic reversal time scale based on model studies of the southwest Pacific basin anomalies is presented. Other model studies indicate that a magnetized layer thickness of at least 2 km is probable. Variations of anomaly amplitudes are studied. A detailed study of the application of numerical correlation techniques to magnetic anomalies is presented. It is concluded that horizontal scale variations and discontinuities in profiles can be critical. Methods for over-coming some of the problems, and for determining quantitative error estimates, are. given. The methods, and conclusions, are applicable to any correlation problem.</p>

Geophysical Studies in the Southwest Pacific : Primarily Studies of Crustal Structure between New Zealand and Antarctica

<p>Geophysical data - primarily magnetic field measurements, bathymetry, and seismicity data - are presented for the area between New Zealand and Antarctica from approximately 145[degrees]W to 155[degrees]E. The data are used to determine the structure of the Pacific-Antarctic boundary, the oceanic part of the Pacific plate and the area of intersection of the Indian, Pacific and Antarctic plates. In the southwest Pacific basin the magnetic anomalies are very clear and an extensive pattern of anomaly lineations with some offsets is mapped. The magnetic anomalies show that the uniform Pacific basin area formed between about 83 and 63 mybp. The Pacific-Antarctic boundary is shown to differ either side of about 175[degrees]W. To the east it is a relatively uniform aseismic spreading ridge, offset by some transform faults. West of 175[degrees]W, to 161[degrees]E, the boundary consists of a seismically active zone of disturbed bathymetry and magnetic anomalies striking about N.70[degrees]W. The zone, the Pacific-Antarctic fracture zone, probably consists of several fractures striking about N45[degrees]W. The area between the Pacific-Antarctic boundary and the southwest Pacific basin represents the interval 10 to -55 mybp, and only in the east are anomaly lineations clear. The Indian-Antarctic Pacific triple junction is near 61.5[degrees]S, 161[degrees]E and is a stable ridge-fault-fault junction; the Indian-Antarctic boundary being the ridge. Plate tectonics is applied to the area and the structure is shown to fit, and be explained by a different rotation pole for each of the major intervals indicated by the structure, i.e. 0-10 mybp, 10-63 mybp and 63-80 mybp. The poles, with rotation rates deduced from the magnetic anomalies, are used to reconstruct the position of New Zealand relative to Antarctica at 80 mybp. The two continents probably started to separate at close to 83 mybp. The times of the major changes of structure and plate movement in the area are shown to coincide with major plate movement changes in the southwest Pacific area and in the rest of the world. A new method for determining poles of rotation, based only on epicentre locations is presented, The method is applied to independently determine the Indian-Pacific, Pacific-Antarctic and Indian-Antarctic poles. The poles should form a consistent. set and they do. The method yields effectively instantaneous poles, is quantitative, and is applicable to most plate boundaries. Earthquake magnitude-frequency relationship b values for the plate boundaries in the area are determined. Comparisons with results from elsewhere indicate an association of high b with high temperature and conversely. Several factors which have previously been suggested as determining b value are shown to not be determinants. A revised and extended magnetic reversal time scale based on model studies of the southwest Pacific basin anomalies is presented. Other model studies indicate that a magnetized layer thickness of at least 2 km is probable. Variations of anomaly amplitudes are studied. A detailed study of the application of numerical correlation techniques to magnetic anomalies is presented. It is concluded that horizontal scale variations and discontinuities in profiles can be critical. Methods for over-coming some of the problems, and for determining quantitative error estimates, are. given. The methods, and conclusions, are applicable to any correlation problem.</p>

Catching the Main Ethiopian Rift evolving towards plate divergence

Magmatism accompanies rifting along divergent plate boundaries, although its role before continental breakup remains poorly understood. For example, the magma-assisted Northern Main Ethiopian Rift (NMER) lacks current volcanism and clear tectono-magmatic relationships with its contiguous rift portions. Here we define its magmatic behaviour, identifying the most recent eruptive fissures (EF) whose aphyric basalts have a higher Ti content than those of older monogenetic scoria cones (MSC), which are porphyritic and plagioclase-dominated. Despite these differences, calculations highlight a similar parental melt for EF and MSC products, suggesting only a different evolutionary history after melt generation. While MSC magmas underwent a further step of storage at intermediate crustal levels, EF magmas rose directly from the base of the crust without contamination, even below older polygenetic volcanoes, suggesting rapid propagation of transcrustal dikes across solidified magma chambers. Whether this recent condition in the NMER is stable or transient, it indicates a transition from central polygenetic to linear fissure volcanism, indicative of increased tensile conditions and volcanism directly fed from the base of the crust, suggesting transition towards mature rifting.