Rift Propagation Research Articles

The influence of firn layer material properties on surface crevasse propagation in glaciers and ice shelves

Abstract. Linear elastic fracture mechanics (LEFM) models have been used to estimate crevasse depths in glaciers and to represent iceberg calving in ice sheet models. However, existing LEFM models assume glacier ice to be homogeneous and utilize the mechanical properties of fully consolidated ice. Using depth-invariant properties is not realistic as the process of compaction from unconsolidated snow to firn to glacial ice is dependent on several environmental factors, typically leading to a lower density and Young's modulus in upper surface strata. New analytical solutions for longitudinal-stress profiles are derived using depth-varying properties based on borehole data from the Ronne Ice Shelf and are used in an LEFM model to determine the maximum penetration depths of an isolated crevasse in grounded glaciers and floating ice shelves. These maximum crevasse depths are compared to those obtained for homogeneous glacial ice, showing the importance of including the effect of the upper unconsolidated firn layers on crevasse propagation. The largest reductions in the penetration depth ratio were observed for shallow grounded glaciers, with variations in Young's modulus being more influential than firn density (maximum differences in crevasse depth of 46 % and 20 %, respectively), whereas firn density changes resulted in an increase in penetration depth for thinner floating ice shelves (95 %–188 % difference in crevasse depth between constant and depth-varying properties). Thus, our study shows that the firn layer can increase the vulnerability of ice shelves to fracture and calving, highlighting the importance of considering depth-dependent firn layer material properties in LEFM models for estimating crevasse penetration depths and predicting rift propagation.

Rift propagation signals the last act of the Thwaites Eastern Ice Shelf despite low basal melt rates

Abstract Rift propagation, rather than basal melt, drives the destabilization and disintegration of the Thwaites Eastern Ice Shelf. Since 2016, rifts have episodically advanced throughout the central ice-shelf area, with rapid propagation events occurring during austral spring. The ice shelf's speed has increased by ~70% during this period, transitioning from a rate of 1.65 m d−1 in 2019 to 2.85 m d−1 by early 2023 in the central area. The increase in longitudinal strain rates near the grounding zone has led to full-thickness rifts and melange-filled gaps since 2020. A recent sea-ice break out has accelerated retreat at the western calving front, effectively separating the ice shelf from what remained of its northwestern pinning point. Meanwhile, a distributed set of phase-sensitive radar measurements indicates that the basal melting rate is generally small, likely due to a widespread robust ocean stratification beneath the ice–ocean interface that suppresses basal melt despite the presence of substantial oceanic heat at depth. These observations in combination with damage modeling show that, while ocean forcing is responsible for triggering the current West Antarctic ice retreat, the Thwaites Eastern Ice Shelf is experiencing dynamic feedbacks over decadal timescales that are driving ice-shelf disintegration, now independent of basal melt.

The orogenic bridge theory: towards a predictive tool for past and future plate tectonics

Wegener’s Continental Drift Theory has laid the foundations of modern plate tectonics. However, despite decades of work and studies around the globe, modern plate tectonics still does not explain all the datasets acquired up to now and is well overdue for a major update. We propose a new theory, the orogenic bridge theory, which partly builds on the Continental Drift Theory and modern plate tectonics and reconciles them with the idea put forward by a competing theory, the Land Bridge Theory (or Isthmian Links). The orogenic bridge theory states that the style of continental rifting is directly controlled by preexisting orogenic structures. On the one hand, preexisting orogens trending parallel to an opening rift facilitate breakup and rift propagation and control the strike and geometry of rift-related structures. This endmember has already been broadly studied worldwide. On the other hand, orogens oriented orthogonal (or highly oblique) to the opening rift will act as strong barriers forcing the rift to step, therefore delaying or impeding breakup and rift propagation and localizing the formation of major-offset transform faults. In the present contribution, we review the evidence in favor of a correlation between rift-orthogonal orogens and major transforms and discuss some of the main implications of the orogenic bridge theory.

Rift Propagation Interacting With Pre‐Existing Microcontinental Blocks

AbstractRift propagation is a 3D thermo‐mechanical process that often precedes continental breakup. Pre‐existing microcontinental blocks and the associated lithospheric strength heterogeneities influence the style of rift propagation. Interestingly, some rifts propagate into pre‐existing blocks and eventually cut through them (e.g., the Zhongsha Block and the Reed Bank), while others bypass these microcontinental blocks forming distinct overlapping rift branches (e.g., the East African Rift System). In this study, we use 3D numerical models to investigate the interaction between microcontinental blocks and rift propagation under different far‐field extension rates. In doing so, we assess the impact of mantle lithospheric thicknesses and lower crustal rheology on the style of rift propagation. Our models reproduce the two types of rift propagation, characterized by propagating rifts that either split or bypass the pre‐existing microcontinental blocks. We find that lithospheric thickness exerts dominant control, while lower crustal rheology of microcontinental blocks and the extension rate have less effect on rift propagation. Our model results can explain rift propagation patterns, block rotation and strong lithospheric thinning in the South China Sea, the East African Rift System, and the Woodlark Basin.

Metamorphic Inheritance, Lower‐Crustal Earthquakes, and Continental Rifting

AbstractThe Malawi Rift is localized within Precambrian amphibolite‐granulite facies metamorphic belts, bounded by up to 150 km long border faults, and generates earthquakes throughout ∼40 km thick crust. Rift‐related faults are inferred to exploit pre‐existing weaknesses that allow rifting of otherwise dry and strong crust. It is unclear what these weaknesses are, and how localization into weak zones can be reconciled with strength required for lower crustal seismicity. We present results of mineral equilibria modeling, which indicate that Proterozoic metamorphism generated dry crust dominated by a quartz‐feldspar assemblage that is metastable at current conditions. For rift propagation to be possible at current cool thermal gradients and in mechanically strong, dry quartzofeldspathic rocks, mid‐ to lower‐crustal strain must be localized into relatively weak, inherited shear zones that deform primarily by aseismic, viscous creep. These shear zones are embedded within high‐strength crust, and interaction between creeping shear zones and enveloped or surrounding rocks may locally increase stress and trigger frictional, seismic slip at mid‐ to lower‐crustal depths. Over time, this interaction may produce a fracture network that allows infiltration of fluids. We therefore suggest that during rifting of previously deformed and metamorphosed crust, major faults are most likely to grow from below, with their location and orientation prescribed by underlying inherited viscous shear zones. In this case, fluids may infiltrate and locally weaken metastable lower crust, including allowing time‐dependent fluid‐driven seismicity and local partial melting, but length‐scales of this weakening is limited by the scale of the permeability network.

Ocean Coupling Limits Rupture Velocity of Fastest Observed Ice Shelf Rift Propagation Event

AbstractThe Antarctic ice sheet is buttressed by floating ice shelves that calve icebergs along large fractures called rifts. Despite the significant influence exerted by rifting on ice shelf geometry and buttressing, the scarcity of in situ observations of rift propagation contributes considerable uncertainty to understanding rift dynamics. Here, we report the first‐ever seismic recording of a multiple‐kilometer rift propagation event. Remote sensing and seismic recordings reveal that a rift in the Pine Island Glacier Ice Shelf extended 10.53 km at a speed of 35.1 m/s, the fastest known ice fracture at this scale. We simulate ocean‐coupled rift propagation and find that the dynamics of water flow within the rift limit the propagation rate, resulting in rupture two orders of magnitude slower than typically predicted for brittle fracture. Using seismic recordings of the elastic waves generated during rift propagation, we estimate that ocean water flows into the rift at a rate of at least 2,300 m3/s during rift propagation and causes mixing in the subshelf cavity. Our observations support the hypotheses that large ice shelf rift propagation events are brittle, hydrodynamically limited, and exhibit sensitive coupling with the surrounding ocean.

A step forward to understanding the development of volcanotectonic rifts: the structure of the Fremrinamar Fissure Swarm (Iceland)

We analysed all the Holocene structures defining the Fremrinamar Fissure Swarm (FFS), in the Northern Volcanic Zone of Iceland, through the interpretation of aerial photos, orthomosaics and Digital Surface Models (DSMs), and field surveys. We measured the strike, dip, length and kinematics of 761 normal faults and reconstructed the slip profile of 76 main faults (length >2 km), with the purpose of evaluating the overall direction of along-axis rift propagation. We also measured the strike of 146 eruptive fissures and 1,128 extension fractures. A total of 421 faults dip towards the east and 340 dip towards the west, mainly striking N0°-10°E. Maximum fault length is 14.2 km, and W-dipping faults are longer than E-dipping faults. The majority of eruptive fissures strike N10°-20°E, and are concentrated in the southern part of the FFS, around the Fremrinamar central volcano. Extension fractures mainly strike N0°-10°E, with a maximum length of 2,508 m. We evaluated the variation of strike, fracture density and spacing along the FFS, and observed a change of its trend from NNE-SSW in the central-southern part, to NNW-SSE in the northern part. We interpret this evidence as the effect of the intersection with the Grimsey Lineament. The tapering of fault slip profiles indicates a main northward propagation of the rift, and thus of the deformation, interpreted as the effect of lateral propagation of dykes from the magma chamber below the central volcano towards the north. Such interpretation is also supported by the distribution of normal faults, vertical offset and dilation values, and also by the rift width, which tend to decrease towards the north.

Ancient silicic volcanism and incipient rift propagation during Gondwana breakup; geophysical insights from the Bumbeni ridge, South Africa

A detailed understanding of the association between magmatism and rift initiation is key to defining the prevailing tectonic conditions active during the initial break-up of Gondwana. The Bumbeni Complex in northern KwaZulu-Natal, South Africa represents an early Cretaceous silicic volcanic centre exposed at the southwestern termination of a linear zone of magnetic anomalies which are buried beneath sedimentary cover. Newly acquired high resolution aeromagnetic data combined with legacy seismic and borehole data is analysed to define the structural architecture of the basement beneath the Zululand Basin in northern KwaZulu-Natal. Multiple prominent positive and negative anomalies ranging from −200 nT to +300 nT based on International Geomagnetic Reference Field (IGRF) and corrected Reduced to Pole (RTP) data are resolved and form a distinct geophysical lineament. Geophysical and borehole core data indicate that the lineament represents a chain of silicic central volcanoes developed as a ridge feature, the Bumbeni Ridge, along the axial zone of an incipient rift structure. Detailed geological mapping and high resolution geophysical data are analysed and compared with analogous contemporary continental rift environments, as well as published data within the adjacent offshore Northern Natal Valley. Results are presented which propose a continental rift evolutionary model for Gondwana breakup in the Zululand Basin in northern KwaZulu-Natal.

Origin of paired extension-compression during rotational rifting: An early Paleogene example from the northeast Atlantic region and its implications

Abstract The formation of paired extension-compression (PEC) postulated by rotational kinematics of rift propagation is demonstrated by analogue models but rarely observed in nature. In our study of the early Paleogene continental rift in the northeast Atlantic, a PEC is proposed based on the northeastward propagation and the coeval compression at the rift tip. The propagation is deduced from tectono-magmatic trends, including along-axis development of magmatism, and migration of tectonic faulting inward and toward the rift tip. Where this propagation terminated, we documented coeval extension and compression in the form of a horst-and-graben system (H&G) and V-shaped anticline (VA), respectively. Given their structural characteristics and spatiotemporal relationships with the rift, their origin is best illustrated by a three-stage model: (1) Rifting initiated at the site of mantle upwelling and propagated northeastward in the Paleocene. (2) The rift tip was stalled by an elevated mafic-ultramafic body at the Barents Shelf, which led to forward projection of the rift’s driving force to create the H&G and the VA (PEC), dissipating the along-axis force component. (3) Domination of axis-perpendicular components then promoted orthogonal extension and sheared margin development. Our study suggests that PEC has a crucial role in both termination of propagation and rift-mode conversion.

Simulating the processes controlling ice-shelf rift paths using damage mechanics

Abstract Rifts are full-thickness fractures that propagate laterally across an ice shelf. They cause ice-shelf weakening and calving of tabular icebergs, and control the initial size of calved icebergs. Here, we present a joint inverse and forward computational modeling framework to capture rifting by combining the vertically integrated momentum balance and anisotropic continuum damage mechanics formulations. We incorporate rift–flank boundary processes to investigate how the rift path is influenced by the pressure on rift–flank walls from seawater, contact between flanks, and ice mélange that may also transmit stress between flanks. To illustrate the viability of the framework, we simulate the final 2 years of rift propagation associated with the calving of tabular iceberg A68 in 2017. We find that the rift path can change with varying ice mélange conditions and the extent of contact between rift flanks. Combinations of parameters associated with slower rift widening rates yield simulated rift paths that best match observations. Our modeling framework lays the foundation for robust simulation of rifting and tabular calving processes, which can enable future studies on ice-sheet–climate interactions, and the effects of ice-shelf buttressing on land ice flow.

Post‐Jurassic brittle deformation of Lusitanian Basin in the context of W Iberia evolution

AbstractA detailed structural analysis of the fracture network exposed in the Jurassic strata is used to reconstruct the Lusitanian Basin's brittle tectonic history related to the Meso‐Cenozoic paleostress trajectories of the Iberian plate. Structural analysis is made by high‐resolution virtual outcrop models and orthophoto mosaics, along with information obtained in the field. The paleostress regime is determined based on the fault‐slip inversion method. Structural features are predominantly NNE‐SSW, NE‐SW and NW‐SE‐trending extensional fractures, including joints, veins, normal faults, and ~E‐W‐oriented strike‐slip faults. These structures remained active in the early basin evolution and were repeatedly reactivated by shearing and contraction. The chronological succession and paleostress reconstruction revealed three tectonic regimes (i) NE‐oriented extension, (ii) NE‐oriented strike‐slip and (iii) NW‐shortening. The first stress regime was driven by the North Atlantic rift propagation in the Iberia's west and northwest margins in the Late Jurassic–Early Cretaceous. The younger stress states involve reactivation and inversion of pre‐existing fractures by Africa–Europe convergence since the Late Cretaceous. The findings are consistent with the regional stress field which the Iberian plate has experienced since the Meso‐Cenozoic.

Strain and stress gradients through the backarc regions of Miocene western Japan: A new type of arc-parallel extension

The overriding lithosphere is occasionally stretched parallel to the consuming plate boundary. Such arc-parallel extension manifests itself mostly in forearc regions. Here, we report unique arc-parallel extension, which affected only the backarc regions with arc-parallel strain and stress gradients. The extension took place in the backarc regions of the northern Ryukyu and western SW Japan arcs in the Middle Miocene. The extension was evidenced by NW-SE trending normal faults in the extensional zone and paleostresses that were determined from dike orientations. To determine the age of and spatial extent of this zone, we analyzed fault-slip data from the Amakusa region, northern Ryukyu arc, where early middle Miocene intrusions allowed us to judge the magmatism to be older than the faults that affected Cretaceous and Eocene formations. The slip senses of the faults have been controversial in Amakusa. As a result, the region was found to have belonged to the extensional zone. The extensional zone is 300–400 km long and ∼100 km wide. We also found that the arc-parallel extension resulted in crustal strain with a southwestward increasing trend along the zone. The strain magnitudes showed a negative correlation with distance from the Okinawa Trough. This correlation is explained by the increasing stress magnitude toward the trough. We suggest that the stalled rift propagation of the northern Okinawa Trough caused the unique arc-parallel extension in the Middle Miocene.

Unravelling rift propagation from Cenozoic subsidence patterns: Examples from the Dangerous Grounds, South China Sea

AbstractHow seafloor spreading in marginal basins controls subsidence patterns in continental margins under rift propagation settings remains poorly understood. Herein, based on detailed geological interpretation and backstripping of 37 high‐resolution multichannel seismic reflection sections, we calculate the Cenozoic tectonic subsidence across the Dangerous Grounds (DG), South China Sea, in order to study interactions between seafloor spreading and the propagation of rift basins. Our results show that rapid subsidence occurred in the eastern Dangerous Grounds (EDG) from the Palaeocene to early Oligocene, then migrated westwards to the middle Dangerous Grounds (MDG) during the early Oligocene to early Miocene, and finally reached the western Dangerous Grounds (WDG) during the early to middle Miocene. We interpret this western migration of intensive subsidence as the response of the inherited pre‐Cenozoic heterogeneous lithospheric structure of South China. Thermal re‐equilibration of the lithosphere then dominated the rapid post‐spreading subsidence throughout the whole DG from the middle Miocene. An anomalous subsidence deficit occurred in the EDG during the early Oligocene to early Miocene, and in the MDG during the early to middle Miocene; this phenomenon might be related to small‐scale mantle convection associated with the propagation of seafloor spreading from east to southwest.

Glaciological history and structural evolution of the Shackleton Ice Shelf system, East Antarctica, over the past 60 years

Abstract. The discovery of Antarctica's deepest subglacial trough beneath the Denman Glacier, combined with high rates of basal melt at the grounding line, has caused significant concern over its vulnerability to retreat. Recent attention has therefore been focusing on understanding the controls driving Denman Glacier's dynamic evolution. Here we consider the Shackleton system, comprised of the Shackleton Ice Shelf, Denman Glacier, and the adjacent Scott, Northcliff, Roscoe and Apfel glaciers, about which almost nothing is known. We widen the context of previously observed dynamic changes in the Denman Glacier to the wider region of the Shackleton system, with a multi-decadal time frame and an improved biannual temporal frequency of observations in the last 7 years (2015–2022). We integrate new satellite observations of ice structure and airborne radar data with changes in ice front position and ice flow velocities to investigate changes in the system. Over the 60-year period of observation we find significant rift propagation on the Shackleton Ice Shelf and Scott Glacier and notable structural changes in the floating shear margins between the ice shelf and the outlet glaciers, as well as features indicative of ice with elevated salt concentration and brine infiltration in regions of the system. Over the period 2017–2022 we observe a significant increase in ice flow speed (up to 50 %) on the floating part of Scott Glacier, coincident with small-scale calving and rift propagation close to the ice front. We do not observe any seasonal variation or significant change in ice flow speed across the rest of the Shackleton system. Given the potential vulnerability of the system to accelerating retreat into the overdeepened, potentially sediment-filled bedrock trough, an improved understanding of the glaciological, oceanographic and geological conditions in the Shackleton system are required to improve the certainty of numerical model predictions, and we identify a number of priorities for future research. With access to these remote coastal regions a major challenge, coordinated internationally collaborative efforts are required to quantify how much the Shackleton region is likely to contribute to sea level rise in the coming centuries.

Hydrothermal Reservoir and Electrical Anisotropy Investigated by Magnetotelluric Data, Case Study of Asal Rift, Republic of Djibouti

At the center of the Republic of Djibouti, an eroded rift called Asal is located where tectonic and magmatic activities can be observed at the surface. Multiple studies were carried out with different exploration methods, such as structural, geophysical and hydrogeological, to understand rifting processes and characterize the subsurface of this rift. Among these subsurface exploration methods, the deep geoelectrical structures need to be better defined with the magnetotelluric (MT) method to better delineate the deep resistivity structures. With the objective of improving our understanding of the deep rift structure, magnetotelluric (MT) data acquired in the Asal rift were analyzed and inverted to build a 2D electrical conductivity model of the hydrothermal system. To achieve this, a dimensionality analysis of the MT data along a 2D profile perpendicular to the rift axis was carried out. Results of this analysis justify the approximation of 2D conductivity structure. Then, 2D inversion models were achieved to build models of the conductive structures. Dimensionality analysis results revealed the existence of electrical anisotropy. Consistent correlation between geoelectric strike and electrical anisotropy direction was suggested. Electrical anisotropy direction determined from the ellipticity of the phase tensor for the short periods was interpreted as the consequence of tectonic activity and horizontal deformation of the rift. Moreover, electrical anisotropy direction for the long periods was assumed to be related to the effects of combined magmatic-tectonic activities with predominant magma/dyke intrusion, which implies the vertical deformation and the subsidence of the rift and may imply the alignment of Olivine. Moreover, the variation and rotation of paleo and recent stress fields direction of plate motion in Asal rift located at the junction of three diverging plates—Arabia, Nubia and Somalia—over geological time can generate both magmatic and tectonic activities which in turn can induce a preferred direction of electrical anisotropy which is the direction of the highest conductivity. While the north-south electrical anisotropy direction is parallel to the direction of Red Sea Rift propagation, the north-east electrical anisotropy direction is aligned with the extension direction between Arabia and Somalia plates. Results of the 2D inversion models presented for the Asal rift allowed to identify two superimposed conductive units close to the surface and are interpreted as a shallow aquifer and a wide potential hydrothermal system. These conductive mediums are overlying a relatively resistive medium. The latter is associated with a magmatic system likely containing hot and/or partly molten rocks. The 2D conductivity model developed in this study could be considered as conceptual model of Asal rift prior to modeling multiphase fluid flow and heat transfer and/or could be used to identify the hydrothermal system for future drilling target depth of geothermal exploration.

Rift propagation in south Tibet controlled by under-thrusting of India: a case study of the Tangra Yumco graben (south Tibet)

Active graben systems in south Tibet and the Himalaya are the surface expression of ongoing east–west extension, although the cause and spatiotemporal evolution of normal faulting is a still a matter of debate. We reconstructed the exhumation history driven by normal faulting in the southern Tangra Yumco graben using new thermochronological data. The Miocene cooling history of the footwall of the main graben-bounding fault is constrained by zircon (U–Th)/He ages (16.7 ± 1.0 to 13.3 ± 0.6 Ma), apatite fission track ages (15.9 ± 2.1 to 13.0 ± 2.1 Ma) and apatite (U–Th)/He ages (7.9 ± 0.4 to 5.3 ± 0.3 Ma). Thermo-kinematic modelling of the data indicates that normal faulting began 19.0 ± 1.1 Ma at a rate ofc.0.2 km myr−1and accelerated toc.0.4 km myr−1atc.5 Ma. In the northern Tangra Yumco rift, remodelling of published data shows that faulting startedc.5 myr later at 13.9 ± 0.8 Ma. The age difference and the distance of 130 km between these two sites indicates that rifting and normal faulting propagated northward at an average rate ofc.25 km myr−1. As this rate is similar to the Miocene convergence rate between India and south Tibet, we argue that the under-thrusting of India beneath Tibet exerted an important control on the propagation of rifts in south Tibet.Supplementary material: Figures S1, S2, and Table S1 with details on apatite fission track analysis are available athttps://doi.org/10.6084/m9.figshare.c.6198584

Rotational Extension Promotes Coeval Upper Crustal Brittle Faulting and Deep‐Seated Rift‐Axis Parallel Flow: Dynamic Coupling Processes Inferred From Analog Model Experiments

AbstractThe lower parts of warm, thick continental crust can flow in a ductile fashion to accommodate thinning of the upper brittle crust during extension. Naturally occurring continental rifts with a rift‐axis parallel deformation gradient imply an underlying rotational component. In such settings, rift‐parallel crustal flow transports material perpendicular to the direction of rifting. We use analog experiments to investigate rotational rifting and coeval crustal flow. To test the effect of rift‐axis parallel flow on rift evolution, we use different gravitational loads resulting in a range of horizontal pressure gradient magnitudes which drive horizontal lower‐crustal flow. The use of (three dimensional) 3D Digital Volume Correlation techniques on X‐ray CT data combined with 3D Digital Image Correlation techniques applied to topographic stereo images provides detailed insights on the contemporaneous evolution of ductile flow patterns and brittle rift structures, respectively. Our results depict a complex flow field in the ductile lower crust during rotational rifting with: (a) extension‐parallel horizontal inward flow and vertical upward flow that compensates thinning of the brittle upper crustal layer; (b) rift‐axis parallel lateral flow, that compensates greater amounts of thinning further away from the rotation axis; and (c) different degrees of mechanical coupling between the brittle and viscous layers that change during rift propagation. Our analog experiments provide insights into ductile lower crustal flow patterns during rift evolution. The results emphasize the three dimensionality of rifting, which is an important effect that should be considered when estimating the amount of crustal extension from two dimensional (2D) cross sections.

Seismotectonics and active faulting of Usangu basin, East African rift system, with implications for the rift propagation

Seismicity and focal mechanisms within the NE-SW extending Usangu basin of the Western branch, north to northeast of the Rungwe Volcanic Province (RVP) southern Tanzania were determined from data recorded on temporary SEGMeNT and permanent AfricaArray seismic stations. The aim of the study has been to investigate stress regime and evaluates the role played by the Usangu basin in accommodating the relative movement between the two sub-basins of the Rukwa rift and the Malawi rift. Results show that most of the seismicity (~70%) concentrates along the Usangu border fault. A lesser amount of seismicity (~30%) is found within the basin and in the southwestern end of the basin towards the RVP. These findings suggest that the Usangu border fault, intrabasinal and transfer faults within the basin are associated with the seismic activities. This implies ongoing active deformation within the basin, which is accommodated by slip on the Usangu border fault system and seismicity associated with intrabasinal faults could indicate the temporal strain migration from the border fault. Although it has previously been suggested that the Eastern branch propagates southeasterly towards the coast, however, given the dense and long deployment of seismic station within the study area the possibility of second propagation of the Eastern branch towards the RVP in the Western branch cannot be ruled out. The seismicity therefore indicates that the Western branch could likely link with the Eastern branch via a wide zone of deformation combining the Mamizi-Kilombero-Ruhuhu fault zone and the zone of Mtera-Ruaha-Usangu rift graben, Focal mechanisms along the Usangu border fault and transfer faults within the basin show a combination of normal faulting and strike slip motion indicating extension with sinistral and dextral components of motion, consistent with palaeostress data and the strike slip deformation within the basin through accommodating the relative movement between the southern-most sub-basin of the Rukwa rift and northern-most sub-basin of the Malawi rift.

Strain localization and migration during the pulsed lateral propagation of the Shire Rift Zone, East Africa

We investigate the spatiotemporal patterns of strain accommodation during multiphase rift evolution in the Shire Rift Zone (SRZ), East Africa. The NW-trending SRZ records a transition from magma-rich rifting phases (Permian-Early Jurassic: Rift-Phase 1 (RP1), and Late Jurassic-Cretaceous: Rift-Phase 2 (RP2)) to a magma-poor phase in the Cenozoic (ongoing: Rift-Phase 3 (RP3)). Our observations show that although the rift border faults largely mimic the pre-rift basement metamorphic fabrics, the rift termination zones occur near crustal-scale rift-orthogonal basement shear zones (Sanangoe (SSZ) and the Lurio shear zones) during RP1-RP2 period. In RP3, the RP1-RP2 sub-basins were largely abandoned, and the rift axes migrated northeastward (rift-orthogonally) into the RP1-RP2 basin margin, and northwestward (strike-parallel) ahead of the RP2 rift-tip. The northwestern RP3 rift-axis side-steps across the SSZ with a rotation of border faults across the shear zone, and terminates farther northwest at another regional-scale shear zone. We suggest that over the multiple pulses of tectonic extension and strain migration in the SRZ, pre-rift basement fabrics acted as: 1) favorably-oriented zones of mechanical strength contrast that localized the large rift faults, and 2) mechanical ‘barriers’ that refracted and possibly, temporarily halted the lateral propagation of the rift zone. Further, the cooled RP1-RP2 mafic dikes localized later-phase deformation in the form of border fault hard-linking transverse faults that exploited strength contrasts within the dike clusters and served as mechanically-strong zones that arrested some of the RP3 fault-tips. Overall, we argue that during pulsed rift propagation, inherited crustal strength anisotropies may serve as both strain-localizing, refracting, and ‘strain barrier’ tectonic structures.

Early Evolution of the Adelaide Superbasin

Continental rifts have a significant role in supercontinent breakup and the development of sedimentary basins. The Australian Adelaide Superbasin is one of the largest and best-preserved rift systems that initiated during the breakup of Rodinia, yet substantial challenges still hinder our understanding of its early evolution and place within the Rodinian supercontinent. In the past decade, our understanding of rift and passive margin development, mantle plumes and their role in tectonics, geodynamics of supercontinent breakup, and sequence stratigraphy in tectonic settings has advanced significantly. However, literature on the early evolution of the Adelaide Superbasin has not been updated to reflect these advancements. Using new detrital zircon age data for provenance, combined with existing literature, we examine the earliest tectonic evolution of the Adelaide Superbasin in the context of our modern understanding of rift system development. A new maximum depositional age of 893 ± 9 Ma from the lowermost stratigraphic unit provides a revised limit on the initiation of sedimentation and rifting within the basin. Our model suggests that the basin evolved through an initial pulse of extension exploiting pre-existing crustal weakness to form half-grabens. Tectonic quiescence and stable subsidence followed, with deposition of a sourceward-shifting facies tract. Emplacement and extrusion of the Willouran Large Igneous Province occurred at c. 830 Ma, initiating a new phase of rifting. This rift renewal led to widespread extension and subsidence with the deposition of the Curdimurka Subgroup, which constitutes the main cyclic rift sequence in the Adelaide Superbasin. Our model suggests that the Adelaide Superbasin formed through rift propagation to an apparent triple junction, rather than apical extension outward from this point. In addition, we provide evidence suggesting a late Mesoproterozoic zircon source to the east of the basin, and show that the lowermost stratigraphy of the Centralian Superbasin, which is thought to be deposited coevally, had different primary detrital sources.