@article{HeckenbachBruneGlerumetal.2021,
  author    = {Heckenbach, Esther Lina and Brune, Sascha and Glerum, Anne C. and Bott, Judith},
  title     = {Is there a speed limit for the thermal steady-state assumption in continental rifts?},
  series = {Geochemistry, geophysics, geosystems : G 3 ; an electronic journal of the earth sciences},
  volume    = {22},
  journal   = {Geochemistry, geophysics, geosystems : G 3 ; an electronic journal of the earth sciences},
  number    = {3},
  publisher = {Wiley},
  address   = {Hoboken, NJ},
  issn      = {1525-2027},
  doi       = {10.1029/2020GC009577},
  pages     = {18},
  year      = {2021},
  abstract  = {The lithosphere is often assumed to reside in a thermal steady-state when quantitatively describing the temperature distribution in continental interiors and sedimentary basins, but also at active plate boundaries. Here, we investigate the applicability limit of this assumption at slowly deforming continental rifts. To this aim, we assess the tectonic thermal imprint in numerical experiments that cover a range of realistic rift configurations. For each model scenario, the deviation from thermal equilibrium is evaluated. This is done by comparing the transient temperature field of every model to a corresponding steady-state model with an identical structural configuration. We find that the validity of the thermal steady-state assumption strongly depends on rift type, divergence velocity, sampling location, and depth within the rift. Maximum differences between transient and steady-state models occur in narrow rifts, at the rift sides, and if the extension rate exceeds 0.5-2 mm/a. Wide rifts, however, reside close to thermal steady-state even for high extension velocities. The transient imprint of rifting appears to be overall negligible for shallow isotherms with a temperature less than 100 degrees C. Contrarily, a steady-state treatment of deep crustal isotherms leads to an underestimation of crustal temperatures, especially for narrow rift settings. Thus, not only relatively fast rifts like the Gulf of Corinth, Red Sea, and Main Ethiopian Rift, but even slow rifts like the Kenya Rift, Rhine Graben, and Rio Grande Rift must be expected to feature a pronounced transient component in the temperature field and to therefore violate the thermal steady-state assumption for deeper crustal isotherms.},
  language  = {en}
}
@article{NeuharthBruneWronaetal.2022,
  author    = {Neuharth, Derek and Brune, Sascha and Wrona, Thilo and Glerum, Anne and Braun, Jean and Yuan, Xiaoping},
  title     = {Evolution of rift systems and their fault networks in response to surface processes},
  series = {Tectonics},
  volume    = {41},
  journal   = {Tectonics},
  number    = {3},
  publisher = {American Geophysical Union},
  address   = {Washington},
  issn      = {0278-7407},
  doi       = {10.1029/2021TC007166},
  pages     = {22},
  year      = {2022},
  abstract  = {Continental rifting is responsible for the generation of major sedimentary basins, both during rift inception and during the formation of rifted continental margins. Geophysical and field studies revealed that rifts feature complex networks of normal faults but the factors controlling fault network properties and their evolution are still matter of debate. Here, we employ high-resolution 2D geodynamic models (ASPECT) including two-way coupling to a surface processes (SP) code (FastScape) to conduct 12 models of major rift types that are exposed to various degrees of erosion and sedimentation. We further present a novel quantitative fault analysis toolbox (Fatbox), which allows us to isolate fault growth patterns, the number of faults, and their length and displacement throughout rift history. Our analysis reveals that rift fault networks may evolve through five major phases: (a) distributed deformation and coalescence, (b) fault system growth, (c) fault system decline and basinward localization, (d) rift migration, and (e) breakup. These phases can be correlated to distinct rifted margin domains. Models of asymmetric rifting suggest rift migration is facilitated through both ductile and brittle deformation within a weak exhumation channel that rotates subhorizontally and remains active at low angles. In sedimentation-starved settings, this channel satisfies the conditions for serpentinization. We find that SP are not only able to enhance strain localization and to increase fault longevity but that they also reduce the total length of the fault system, prolong rift phases and delay continental breakup.},
  language  = {en}
}
@article{PonsSobolevLiuetal.2022,
  author    = {Pons, Micha{\"e}l and Sobolev, Stephan and Liu, Sibiao and Neuharth, Derek},
  title     = {Hindered trench migration due to slab steepening controls the formation of the Central Andes},
  series = {JGR : Solid earth},
  volume    = {127},
  journal   = {JGR : Solid earth},
  number    = {12},
  publisher = {Wiley},
  address   = {Hoboken, NJ},
  issn      = {2169-9313},
  doi       = {10.1029/2022JB025229},
  pages     = {21},
  year      = {2022},
  abstract  = {The formation of the Central Andes dates back to similar to 50 Ma, but its most pronounced episode, including the growth of the Altiplano-Puna Plateau and pulsatile tectonic shortening phases, occurred within the last 25 Ma. The reason for this evolution remains unexplained. Using geodynamic numerical modeling we infer that the primary cause of the pulses of tectonic shortening and growth of the Central Andes is the changing geometry of the subducted Nazca plate, and particularly the steepening of the mid-mantle slab segment which results in a slowing down of the trench retreat and subsequent increase in shortening of the advancing South America plate. This steepening first happens after the end of the flat slab episode at similar to 25 Ma, and later during the buckling and stagnation of the slab in the mantle transition zone. Processes that mechanically weaken the lithosphere of the South America plate, as suggested in previous studies, enhance the intensity of the shortening events. These processes include delamination of the mantle lithosphere and weakening of foreland sediments. Our new modeling results are consistent with the timing and amplitude of the deformation from geological data in the Central Andes at the Altiplano latitude. Plain Language Summary The Central Andes is a subduction-type orogeny that formed as a result of the interaction between the Nazca oceanic plate and the South American continental plate over the last 50 million years. Growth of the Andes is primarily the result of crustal shortening. Nevertheless, "geological" data compiled from previous studies have shown that phases of drastic pulsatile shortening occur at 15 and 5 Ma. In this study, we used high-resolution 2D numerical geodynamic simulations to investigate the link between oceanic and continental plate dynamics and their interaction. We find that when the oceanic plate steepens in the mantle transition zone, the trench retreat is hindered. Coupled with the weakening of the continental plate through the slab flattening and subsequent delamination of the lithospheric mantle, this leads to pulsatile shortening phases of a magnitude equivalent to that suggested by the data.},
  language  = {en}
}