@phdthesis{RodriguezPiceda2022,
  author    = {Rodriguez Piceda, Constanza},
  title     = {Thermomechanical state of the southern Central Andes},
  doi       = {10.25932/publishup-54927},
  url       = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-549275},
  school      = {Universit{\"a}t Potsdam},
  pages     = {xx, 228},
  year      = {2022},
  abstract  = {The Andes are a ~7000 km long N-S trending mountain range developed along the South American western continental margin. Driven by the subduction of the oceanic Nazca plate beneath the continental South American plate, the formation of the northern and central parts of the orogen is a type case for a non-collisional orogeny. In the southern Central Andes (SCA, 29°S-39°S), the oceanic plate changes the subduction angle between 33°S and 35°S from almost horizontal (< 5° dip) in the north to a steeper angle (~30° dip) in the south. This sector of the Andes also displays remarkable along- and across- strike variations of the tectonic deformation patterns. These include a systematic decrease of topographic elevation, of crustal shortening and foreland and orogenic width, as well as an alternation of the foreland deformation style between thick-skinned and thin-skinned recorded along- and across the strike of the subduction zone. Moreover, the SCA are a very seismically active region. The continental plate is characterized by a relatively shallow seismicity (< 30 km depth) which is mainly focussed at the transition from the orogen to the lowland areas of the foreland and the forearc; in contrast, deeper seismicity occurs below the interiors of the northern foreland. Additionally, frequent seismicity is also recorded in the shallow parts of the oceanic plate and in a sector of the flat slab segment between 31°S and 33°S. The observed spatial heterogeneity in tectonic and seismic deformation in the SCA has been attributed to multiple causes, including variations in sediment thickness, the presence of inherited structures and changes in the subduction angle of the oceanic slab. However, there is no study that inquired the relationship between the long-term rheological configuration of the SCA and the spatial deformation patterns. Moreover, the effects of the density and thickness configuration of the continental plate and of variations in the slab dip angle in the rheological state of the lithosphere have been not thoroughly investigated yet. Since rheology depends on composition, pressure and temperature, a detailed characterization of the compositional, structural and thermal fields of the lithosphere is needed. Therefore, by using multiple geophysical approaches and data sources, I constructed the following 3D models of the SCA lithosphere: (i) a seismically-constrained structural and density model that was tested against the gravity field; (ii) a thermal model integrating the conversion of mantle shear-wave velocities to temperature with steady-state conductive calculations in the uppermost lithosphere (< 50 km depth), validated by temperature and heat-flow measurements; and (iii) a rheological model of the long-term lithospheric strength using as input the previously-generated models. The results of this dissertation indicate that the present-day thermal and rheological fields of the SCA are controlled by different mechanisms at different depths. At shallow depths (< 50 km), the thermomechanical field is modulated by the heterogeneous composition of the continental lithosphere. The overprint of the oceanic slab is detectable where the oceanic plate is shallow (< 85 km depth) and the radiogenic crust is thin, resulting in overall lower temperatures and higher strength compared to regions where the slab is steep and the radiogenic crust is thick. At depths > 50 km, largest temperatures variations occur where the descending slab is detected, which implies that the deep thermal field is mainly affected by the slab dip geometry. The outcomes of this thesis suggests that long-term thermomechanical state of the lithosphere influences the spatial distribution of seismic deformation. Most of the seismicity within the continental plate occurs above the modelled transition from brittle to ductile conditions. Additionally, there is a spatial correlation between the location of these events and the transition from the mechanically strong domains of the forearc and foreland to the weak domain of the orogen. In contrast, seismicity within the oceanic plate is also detected where long-term ductile conditions are expected. I therefore analysed the possible influence of additional mechanisms triggering these earthquakes, including the compaction of sediments in the subduction interface and dehydration reactions in the slab. To that aim, I carried out a qualitative analysis of the state of hydration in the mantle using the ratio between compressional- and shear-wave velocity (vp/vs ratio) from a previous seismic tomography. The results from this analysis indicate that the majority of the seismicity spatially correlates with hydrated areas of the slab and overlying continental mantle, with the exception of the cluster within the flat slab segment. In this region, earthquakes are likely triggered by flexural processes where the slab changes from a flat to a steep subduction angle. First-order variations in the observed tectonic patterns also seem to be influenced by the thermomechanical configuration of the lithosphere. The mechanically strong domains of the forearc and foreland, due to their resistance to deformation, display smaller amounts of shortening than the relatively weak orogenic domain. In addition, the structural and thermomechanical characteristics modelled in this dissertation confirm previous analyses from geodynamic models pointing to the control of the observed heterogeneities in the orogen and foreland deformation style. These characteristics include the lithospheric and crustal thickness, the presence of weak sediments and the variations in gravitational potential energy. Specific conditions occur in the cold and strong northern foreland, which is characterized by active seismicity and thick-skinned structures, although the modelled crustal strength exceeds the typical values of externally-applied tectonic stresses. The additional mechanisms that could explain the strain localization in a region that should resist deformation are: (i) increased tectonic forces coming from the steepening of the slab and (ii) enhanced weakening along inherited structures from pre-Andean deformation events. Finally, the thermomechanical conditions of this sector of the foreland could be a key factor influencing the preservation of the flat subduction angle at these latitudes of the SCA.},
  language  = {en}
}
@article{SpoonerScheckWenderothCacaceetal.2022,
  author    = {Spooner, Cameron and Scheck-Wenderoth, Magdalena and Cacace, Mauro and Anikiev, Denis},
  title     = {How Alpine seismicity relates to lithospheric strength},
  series = {International journal of earth sciences},
  volume    = {111},
  journal   = {International journal of earth sciences},
  number    = {4},
  publisher = {Springer},
  address   = {Berlin ; Heidelberg},
  issn      = {1437-3254},
  doi       = {10.1007/s00531-022-02174-5},
  pages     = {1201 -- 1221},
  year      = {2022},
  abstract  = {Despite the amount of research focussed on the Alpine orogen, different hypotheses still exist regarding varying spatial seismicity distribution patterns throughout the region. Previous measurement-constrained regional 3D models of lithospheric density distribution and thermal field facilitate the generation of a data-based rheological model of the region. In this study, we compute the long-term lithospheric strength and compare its spatial variation to observed seismicity patterns. We demonstrate how strength maxima within the crust (similar to 1 GPa) and upper mantle (> 2 GPa) occur at temperatures characteristic of the onset of crystal plasticity in those rocks (crust: 200-400 degrees C; mantle: similar to 600 degrees C), with almost all seismicity occurring in these regions. Correlation in the northern and southern forelands between crustal and lithospheric strengths and seismicity show different patterns of event distribution, reflecting their different tectonic settings. Seismicity in the plate boundary setting of the southern foreland corresponds to the integrated lithospheric strength, occurring mainly in the weaker domains surrounding the strong Adriatic plate. In the intraplate setting of the northern foreland, seismicity correlates to modelled crustal strength, and it mainly occurs in the weaker and warmer crust beneath the Upper Rhine Graben. We, therefore, suggest that seismicity in the upper crust is linked to weak crustal domains, which are more prone to localise deformation promoting failure and, depending on the local properties of the fault, earthquakes at relatively lower levels of accumulated stress than their neighbouring stronger counterparts. Upper mantle seismicity at depths greater than modelled brittle conditions, can be either explained by embrittlement of the mantle due to grain-size sensitive deformation within domains of active or recent slab cooling, or by dissipative weakening mechanisms, such as thermal runaway from shear heating and/or dehydration reactions within an overly ductile mantle. Results generated in this study are available for open access use to further discussions on the region.},
  language  = {en}
}
@article{BarbotWeiss2021,
  author    = {Barbot, Sylvain and Weiss, Jonathan R.},
  title     = {Connecting subduction, extension and shear localization across the Aegean Sea and Anatolia},
  series = {Geophysical journal international},
  volume    = {226},
  journal   = {Geophysical journal international},
  number    = {1},
  publisher = {Blackwell},
  address   = {Oxford [u.a.]},
  issn      = {0956-540X},
  doi       = {10.1093/gji/ggab078},
  pages     = {422 -- 445},
  year      = {2021},
  abstract  = {The Eastern Mediterranean is the most seismically active region in Europe due to the complex interactions of the Arabian, African, and Eurasian tectonic plates. Deformation is achieved by faulting in the brittle crust, distributed flow in the viscoelastic lower-crust and mantle, and Hellenic subduction, but the long-term partitioning of these mechanisms is still unknown. We exploit an extensive suite of geodetic observations to build a kinematic model connecting strike-slip deformation, extension, subduction, and shear localization across Anatolia and the Aegean Sea by mapping the distribution of slip and strain accumulation on major active geological structures. We find that tectonic escape is facilitated by a plate-boundary-like, translithospheric shear zone extending from the Gulf of Evia to the Turkish-Iranian Plateau that underlies the surface trace of the North Anatolian Fault. Additional deformation in Anatolia is taken up by a series of smaller-scale conjugate shear zones that reach the upper mantle, the largest of which is located beneath the East Anatolian Fault. Rapid north-south extension in the western part of the system, driven primarily by Hellenic Trench retreat, is accommodated by rotation and broadening of the North Anatolian mantle shear zone from the Sea of Marmara across the north Aegean Sea, and by a system of distributed transform faults and rifts including the rapidly extending Gulf of Corinth in central Greece and the active grabens of western Turkey. Africa-Eurasia convergence along the Hellenic Arc occurs at a median rate of 49.8mm yr(-1) in a largely trench-normal direction except near eastern Crete where variably oriented slip on the megathrust coincides with mixed-mode and strike-slip deformation in the overlying accretionary wedge near the Ptolemy-Pliny-Strabo trenches. Our kinematic model illustrates the competing roles the North Anatolian mantle shear zone, Hellenic Trench, overlying mantle wedge, and active crustal faults play in accommodating tectonic indentation, slab rollback and associated Aegean extension. Viscoelastic flow in the lower crust and upper mantle dominate the surface velocity field across much of Anatolia and a clear transition to megathrust-related slab pull occurs in western Turkey, the Aegean Sea and Greece. Crustal scale faults and the Hellenic wedge contribute only a minor amount to the large-scale, regional pattern of Eastern Mediterranean interseismic surface deformation.},
  language  = {en}
}
@misc{RubeyBruneHeineetal.2017,
  author    = {Rubey, Michael and Brune, Sascha and Heine, Christian and Davies, D. Rhodri and Williams, Simon E. and M{\"u}ller, R. Dietmar},
  title     = {Global patterns in Earth's dynamic topography since the Jurassic},
  series = {Postprints der Universit{\"a}t Potsdam : Mathematisch Naturwissenschaftliche Reihe},
  journal   = {Postprints der Universit{\"a}t Potsdam : Mathematisch Naturwissenschaftliche Reihe},
  number    = {623},
  issn      = {1866-8372},
  doi       = {10.25932/publishup-41824},
  url       = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-418241},
  pages     = {899 -- 919},
  year      = {2017},
  abstract  = {We evaluate the spatial and temporal evolution of Earth's long-wavelength surface dynamic topography since the Jurassic using a series of high-resolution global mantle convection models. These models are Earth-like in terms of convective vigour, thermal structure, surface heat-flux and the geographic distribution of heterogeneity. The models generate a degree-2-dominated spectrum of dynamic topography with negative amplitudes above subducted slabs (i.e. circum-Pacific regions and southern Eurasia) and positive amplitudes elsewhere (i.e. Africa, north-western Eurasia and the central Pacific). Model predictions are compared with published observations and subsidence patterns from well data, both globally and for the Australian and southern African regions. We find that our models reproduce the long-wavelength component of these observations, although observed smaller-scale variations are not reproduced. We subsequently define "geodynamic rules" for how different surface tectonic settings are affected by mantle processes: (i) locations in the vicinity of a subduction zone show large negative dynamic topography amplitudes; (ii) regions far away from convergent margins feature long-term positive dynamic topography; and (iii) rapid variations in dynamic support occur along the margins of overriding plates (e.g. the western US) and at points located on a plate that rapidly approaches a subduction zone (e.g. India and the Arabia Peninsula). Our models provide a predictive quantitative framework linking mantle convection with plate tectonics and sedimentary basin evolution, thus improving our understanding of how subduction and mantle convection affect the spatio-temporal evolution of basin architecture.},
  language  = {en}
}
@phdthesis{Schuetz2013,
  author    = {Sch{\"u}tz, Felina},
  title     = {Surface heat flow and lithospheric thermal structure of the northwestern Arabian Plate},
  url       = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus-69622},
  school      = {Universit{\"a}t Potsdam},
  year      = {2013},
  abstract  = {The surface heat flow (qs) is paramount for modeling the thermal structure of the lithosphere. Changes in the qs over a distinct lithospheric unit are normally directly reflecting changes in the crustal composition and therewith the radiogenic heat budget (e.g., Rudnick et al., 1998; F{\"o}rster and F{\"o}rster, 2000; Mareschal and Jaupart, 2004; Perry et al., 2006; Hasterok and Chapman, 2011, and references therein) or, less usual, changes in the mantle heat flow (e.g., Pollack and Chapman, 1977). Knowledge of this physical property is therefore of great interest for both academic research and the energy industry. The present study focuses on the qs of central and southern Israel as part of the Sinai Microplate (SM). Having formed during Oligocene to Miocene rifting and break-up of the African and Arabian plates, the SM is characterized by a young and complex tectonic history. Resulting from the time thermal diffusion needs to pass through the lithosphere, on the order of several tens-of-millions of years (e.g., Fowler, 1990); qs-values of the area reflect conditions of pre-Oligocene times. The thermal structure of the lithosphere beneath the SM in general, and south-central Israel in particular, has remained poorly understood. To address this problem, the two parameters needed for the qs determination were investigated. Temperature measurements were made at ten pre-existing oil and water exploration wells, and the thermal conductivity of 240 drill core and outcrop samples was measured in the lab. The thermal conductivity is the sensitive parameter in this determination. Lab measurements were performed on both, dry and water-saturated samples, which is labor- and time-consuming. Another possibility is the measurement of thermal conductivity in dry state and the conversion to a saturated value by using mean model approaches. The availability of a voluminous and diverse dataset of thermal conductivity values in this study allowed (1) in connection with the temperature gradient to calculate new reliable qs values and to use them to model the thermal pattern of the crust in south-central Israel, prior to young tectonic events, and (2) in connection with comparable datasets, controlling the quality of different mean model approaches for indirect determination of bulk thermal conductivity (BTC) of rocks. The reliability of numerically derived BTC values appears to vary between different mean models, and is also strongly dependent upon sample lithology. Yet, correction algorithms may significantly reduce the mismatch between measured and calculated conductivity values based on the different mean models. Furthermore, the dataset allowed the derivation of lithotype-specific conversion equations to calculate the water-saturated BTC directly from data of dry-measured BTC and porosity (e.g., well log derived porosity) with no use of any mean model and thus provide a suitable tool for fast analysis of large datasets. The results of the study indicate that the qs in the study area is significantly higher than previously assumed. The new presented qs values range between 50 and 62 mW m⁻². A weak trend of decreasing heat flow can be identified from the east to the west (55-50 mW m⁻²), and an increase from the Dead Sea Basin to the south (55-62 mW m⁻²). The observed range can be explained by variation in the composition (heat production) of the upper crust, accompanied by more systematic spatial changes in its thickness. The new qs data then can be used, in conjunction with petrophysical data and information on the structure and composition of the lithosphere, to adjust a model of the pre-Oligocene thermal state of the crust in south-central Israel. The 2-D steady-state temperature model was calculated along an E-W traverse based on the DESIRE seismic profile (Mechie et al., 2009). The model comprises the entire lithosphere down to the lithosphere-asthenosphere boundary (LAB) involving the most recent knowledge of the lithosphere in pre-Oligocene time, i.e., prior to the onset of rifting and plume-related lithospheric thermal perturbations. The adjustment of modeled and measured qs allows conclusions about the pre-Oligocene LAB-depth. After the best fitting the most likely depth is 150 km which is consistent with estimations made in comparable regions of the Arabian Shield. It therefore comprises the first ever modelled pre-Oligocene LAB depth, and provides important clues on the thermal state of lithosphere before rifting. This, in turn, is vital for a better understanding of the (thermo)-dynamic processes associated with lithosphere extension and continental break-up.},
  language  = {en}
}
@article{KoulakovSobolevWeberetal.2006,
  author    = {Koulakov, Ivan and Sobolev, Stephan Vladimir and Weber, Bernd and Oreshin, Sergey and Wylegalla, Kurt and Hofstetter, Rami},
  title     = {Teleseismic tomography reveals no signature of the Dead Sea Transform in the upper mantle structure},
  series = {Earth and planetary science letters},
  volume    = {252},
  journal   = {Earth and planetary science letters},
  number    = {1-2},
  publisher = {Elsevier},
  address   = {Amsterdam},
  issn      = {0012-821X},
  doi       = {10.1016/j.epsl.2006.09.039},
  pages     = {189 -- 200},
  year      = {2006},
  abstract  = {We present results of a tomographic inversion of teleseismic data recorded at 48 stations of a temporary network which was installed in the area of the Dead Sea Transform (DST) and operated for 1 yr in the framework of the multidisciplinary DESERT Project. The 3366 teleseismic P and PKP phases from 135 events were hand picked and corrected for surface topography and crustal thickness. The inversion shows pronounced low-velocity anomalies in the crust, beneath the DST, which are consistent with recent results from local-source tomography. These anomalies are likely related to the young sediments and fractured rocks in the fault zone. The deeper the retrieved anomalies are quite weak. Most prominent is the high-velocity strip-like anomaly striking SE-NW. We attribute this anomaly to the inherited heterogeneity of lithospheric structure, with a possible contribution by the shallow Precambrian basement east of the DST and to lower crustal heterogeneity reported in this region by other seismic studies. We do not observe reliable signature of the DST in the upper mantle structure. Some weak indications of low-velocity anomalies in the upper mantle beneath the DST may well result from the down-smearing of the strong upper crustal anomalies. We also see very little topography of the lithosphere-asthenosphere boundary beneath the DST, which would generate significant horizontal velocity variations. These results are consistent with predictions from a recent thereto-mechanical model of the DST. Our tomographic model provides some indication of hot mantle flow from the deeper upper mantle rooted in the region of the Red Sea. However, resolution tests show that this anomaly may well be beyond resolution of the model. (c) 2006 Elsevier B.V. All rights reserved.},
  language  = {en}
}