@misc{KloseChaparroSchillingetal.2020, author = {Klose, Tim and Chaparro, M. Carme and Schilling, Frank and Butscher, Christoph and Klumbach, Steffen and Blum, Philipp}, title = {Fluid flow simulations of a large-scale borehole leakage experiment}, series = {Zweitver{\"o}ffentlichungen der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, journal = {Zweitver{\"o}ffentlichungen der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, number = {1}, issn = {1866-8372}, doi = {10.25932/publishup-57353}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-573539}, pages = {23}, year = {2020}, abstract = {Borehole leakage is a common and complex issue. Understanding the fluid flow characteristics of a cemented area inside a borehole is crucial to monitor and quantify the wellbore integrity as well as to find solutions to minimise existing leakages. In order to improve our understanding of the flow behaviour of cemented boreholes, we investigated experimental data of a large-scale borehole leakage tests by means of numerical modelling using three different conceptual models. The experiment was performed with an autoclave system consisting of two vessels bridged by a cement-filled casing. After a partial bleed-off at the well-head, a sustained casing pressure was observed due to fluid flow through the cementsteel composite. The aim of our simulations is to investigate and quantify the permeability of the cement-steel composite. From our model results, we conclude that the flow occurred along a preferential flow path at the cement-steel interface. Thus, the inner part of the cement core was impermeable during the duration of the experiment. The preferential flow path can be described as a highly permeable and highly porous area with an aperture of about 5 mu m and a permeability of 3 . 10(-12) m(2) (3 Darcy). It follows that the fluid flow characteristics of a cemented area inside a borehole cannot be described using one permeability value for the entire cement-steel composite. Furthermore, it can be concluded that the quality of the cement and the filling process regarding the cement-steel interface is crucial to minimize possible well leakages.}, language = {en} } @misc{RodriguezPicedaScheckWenderothGomezDacaletal.2020, author = {Rodriguez Piceda, Constanza and Scheck Wenderoth, Magdalena and Gomez Dacal, Maria Laura and Bott, Judith and Prezzi, Claudia Beatriz and Strecker, Manfred}, title = {Lithospheric density structure of the southern Central Andes constrained by 3D data-integrative gravity modelling}, series = {Zweitver{\"o}ffentlichungen der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, journal = {Zweitver{\"o}ffentlichungen der Universit{\"a}t Potsdam : Mathematisch-Naturwissenschaftliche Reihe}, number = {7}, issn = {1866-8372}, doi = {10.25932/publishup-56262}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-562628}, pages = {29}, year = {2020}, abstract = {The southern Central Andes (SCA) (between 27 degrees S and 40 degrees S) is bordered to the west by the convergent margin between the continental South American Plate and the oceanic Nazca Plate. The subduction angle along this margin is variable, as is the deformation of the upper plate. Between 33 degrees S and 35 degrees S, the subduction angle of the Nazca plate increases from sub-horizontal (< 5 degrees) in the north to relatively steep (similar to 30 degrees) in the south. The SCA contain inherited lithological and structural heterogeneities within the crust that have been reactivated and overprinted since the onset of subduction and associated Cenozoic deformation within the Andean orogen. The distribution of the deformation within the SCA has often been attributed to the variations in the subduction angle and the reactivation of these inherited heterogeneities. However, the possible influence that the thickness and composition of the continental crust have had on both short-term and long-term deformation of the SCA is yet to be thoroughly investigated. For our investigations, we have derived density distributions and thicknesses for various layers that make up the lithosphere and evaluated their relationships with tectonic events that occurred over the history of the Andean orogeny and, in particular, investigated the short- and long-term nature of the present-day deformation processes. We established a 3D model of lithosphere beneath the orogen and its foreland (29 degrees S-39 degrees S) that is consistent with currently available geological and geophysical data, including the gravity data. The modelled crustal configuration and density distribution reveal spatial relationships with different tectonic domains: the crystalline crust in the orogen (the magmatic arc and the main orogenic wedge) is thicker (similar to 55 km) and less dense (similar to 2900 kg/m(3)) than in the forearc (similar to 35 km, similar to 2975 kg/m(3)) and foreland (similar to 30 km, similar to 3000 kg/m(3)). Crustal thickening in the orogen probably occurred as a result of stacking of low-density domains, while density and thickness variations beneath the forearc and foreland most likely reflect differences in the tectonic evolution of each area following crustal accretion. No clear spatial relationship exists between the density distribution within the lithosphere and previously proposed boundaries of crustal terranes accreted during the early Paleozoic. Areas with ongoing deformation show a spatial correlation with those areas that have the highest topographic gradients and where there are abrupt changes in the average crustal-density contrast. This suggests that the short-term deformation within the interior of the Andean orogen and its foreland is fundamentally influenced by the crustal composition and the relative thickness of different crustal layers. A thicker, denser, and potentially stronger lithosphere beneath the northern part of the SCA foreland is interpreted to have favoured a strong coupling between the Nazca and South American plates, facilitating the development of a sub-horizontal slab.}, language = {en} } @article{RolphOverduinRavensetal.2022, author = {Rolph, Rebecca and Overduin, Pier Paul and Ravens, Thomas and Lantuit, Hugues and Langer, Moritz}, title = {ArcticBeach v1.0}, series = {Frontiers in Earth Science}, volume = {10}, journal = {Frontiers in Earth Science}, publisher = {Frontiers Media}, address = {Lausanne}, issn = {2296-6463}, doi = {10.3389/feart.2022.962208}, pages = {19}, year = {2022}, abstract = {In the Arctic, air temperatures are increasing and sea ice is declining, resulting in larger waves and a longer open water season, all of which intensify the thaw and erosion of ice-rich coasts. Climate change has been shown to increase the rate of Arctic coastal erosion, causing problems for Arctic cultural heritage, existing industrial, military, and civil infrastructure, as well as changes in nearshore biogeochemistry. Numerical models that reproduce historical and project future Arctic erosion rates are necessary to understand how further climate change will affect these problems, and no such model yet exists to simulate the physics of erosion on a pan-Arctic scale. We have coupled a bathystrophic storm surge model to a simplified physical erosion model of a permafrost coastline. This Arctic erosion model, called ArcticBeach v1.0, is a first step toward a physical parameterization of Arctic shoreline erosion for larger-scale models. It is forced by wind speed and direction, wave period and height, sea surface temperature, all of which are masked during times of sea ice cover near the coastline. Model tuning requires observed historical retreat rates (at least one value), as well as rough nearshore bathymetry. These parameters are already available on a pan-Arctic scale. The model is validated at three study sites at 1) Drew Point (DP), Alaska, 2) Mamontovy Khayata (MK), Siberia, and 3) Veslebogen Cliffs, Svalbard. Simulated cumulative retreat rates for DP and MK respectively (169 and 170 m) over the time periods studied at each site (2007-2016, and 1995-2018) are found to the same order of magnitude as observed cumulative retreat (172 and 120 m). The rocky Veslebogen cliffs have small observed cumulative retreat rates (0.05 m over 2014-2016), and our model was also able to reproduce this same order of magnitude of retreat (0.08 m). Given the large differences in geomorphology between the study sites, this study provides a proof-of-concept that ArcticBeach v1.0 can be applied on very different permafrost coastlines. ArcticBeach v1.0 provides a promising starting point to project retreat of Arctic shorelines, or to evaluate historical retreat in places that have had few observations.}, language = {en} } @article{KloseChaparroSchillingetal.2020, author = {Klose, Tim and Chaparro, M. Carme and Schilling, Frank and Butscher, Christoph and Klumbach, Steffen and Blum, Philipp}, title = {Fluid flow simulations of a large-scale borehole leakage experiment}, series = {Transport in Porous Media}, volume = {136}, journal = {Transport in Porous Media}, number = {1}, publisher = {Springer}, address = {New York}, issn = {0169-3913}, doi = {10.1007/s11242-020-01504-y}, pages = {125 -- 145}, year = {2020}, abstract = {Borehole leakage is a common and complex issue. Understanding the fluid flow characteristics of a cemented area inside a borehole is crucial to monitor and quantify the wellbore integrity as well as to find solutions to minimise existing leakages. In order to improve our understanding of the flow behaviour of cemented boreholes, we investigated experimental data of a large-scale borehole leakage tests by means of numerical modelling using three different conceptual models. The experiment was performed with an autoclave system consisting of two vessels bridged by a cement-filled casing. After a partial bleed-off at the well-head, a sustained casing pressure was observed due to fluid flow through the cementsteel composite. The aim of our simulations is to investigate and quantify the permeability of the cement-steel composite. From our model results, we conclude that the flow occurred along a preferential flow path at the cement-steel interface. Thus, the inner part of the cement core was impermeable during the duration of the experiment. The preferential flow path can be described as a highly permeable and highly porous area with an aperture of about 5 mu m and a permeability of 3 . 10(-12) m(2) (3 Darcy). It follows that the fluid flow characteristics of a cemented area inside a borehole cannot be described using one permeability value for the entire cement-steel composite. Furthermore, it can be concluded that the quality of the cement and the filling process regarding the cement-steel interface is crucial to minimize possible well leakages.}, language = {en} } @article{RodriguezPicedaScheckWenderothGomezDacaletal.2020, author = {Rodriguez Piceda, Constanza and Scheck Wenderoth, Magdalena and Gomez Dacal, Maria Laura and Bott, Judith and Prezzi, Claudia Beatriz and Strecker, Manfred}, title = {Lithospheric density structure of the southern Central Andes constrained by 3D data-integrative gravity modelling}, series = {International journal of earth sciences}, volume = {110}, journal = {International journal of earth sciences}, number = {7}, publisher = {Springer}, address = {New York}, issn = {1437-3254}, doi = {10.1007/s00531-020-01962-1}, pages = {2333 -- 2359}, year = {2020}, abstract = {The southern Central Andes (SCA) (between 27 degrees S and 40 degrees S) is bordered to the west by the convergent margin between the continental South American Plate and the oceanic Nazca Plate. The subduction angle along this margin is variable, as is the deformation of the upper plate. Between 33 degrees S and 35 degrees S, the subduction angle of the Nazca plate increases from sub-horizontal (< 5 degrees) in the north to relatively steep (similar to 30 degrees) in the south. The SCA contain inherited lithological and structural heterogeneities within the crust that have been reactivated and overprinted since the onset of subduction and associated Cenozoic deformation within the Andean orogen. The distribution of the deformation within the SCA has often been attributed to the variations in the subduction angle and the reactivation of these inherited heterogeneities. However, the possible influence that the thickness and composition of the continental crust have had on both short-term and long-term deformation of the SCA is yet to be thoroughly investigated. For our investigations, we have derived density distributions and thicknesses for various layers that make up the lithosphere and evaluated their relationships with tectonic events that occurred over the history of the Andean orogeny and, in particular, investigated the short- and long-term nature of the present-day deformation processes. We established a 3D model of lithosphere beneath the orogen and its foreland (29 degrees S-39 degrees S) that is consistent with currently available geological and geophysical data, including the gravity data. The modelled crustal configuration and density distribution reveal spatial relationships with different tectonic domains: the crystalline crust in the orogen (the magmatic arc and the main orogenic wedge) is thicker (similar to 55 km) and less dense (similar to 2900 kg/m(3)) than in the forearc (similar to 35 km, similar to 2975 kg/m(3)) and foreland (similar to 30 km, similar to 3000 kg/m(3)). Crustal thickening in the orogen probably occurred as a result of stacking of low-density domains, while density and thickness variations beneath the forearc and foreland most likely reflect differences in the tectonic evolution of each area following crustal accretion. No clear spatial relationship exists between the density distribution within the lithosphere and previously proposed boundaries of crustal terranes accreted during the early Paleozoic. Areas with ongoing deformation show a spatial correlation with those areas that have the highest topographic gradients and where there are abrupt changes in the average crustal-density contrast. This suggests that the short-term deformation within the interior of the Andean orogen and its foreland is fundamentally influenced by the crustal composition and the relative thickness of different crustal layers. A thicker, denser, and potentially stronger lithosphere beneath the northern part of the SCA foreland is interpreted to have favoured a strong coupling between the Nazca and South American plates, facilitating the development of a sub-horizontal slab.}, language = {en} } @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{MarkowskaFohlmeisterTrebleetal.2019, author = {Markowska, Monika and Fohlmeister, Jens Bernd and Treble, Pauline C. and Baker, Andy and Andersen, Martin S. and Hua, Quan}, title = {Modelling the C-14 bomb-pulse in young speleothems using a soil carbon continuum model}, series = {Geochimica et cosmochimica acta : journal of the Geochemical Society and the Meteoritical Society}, volume = {261}, journal = {Geochimica et cosmochimica acta : journal of the Geochemical Society and the Meteoritical Society}, publisher = {Elsevier}, address = {Oxford}, issn = {0016-7037}, doi = {10.1016/j.gca.2019.04.029}, pages = {342 -- 367}, year = {2019}, abstract = {The 'bomb-pulse' method is a chronological approach to further constrain the age of speleothems that grew between 1950 CE - present. Establishing dependable chronological constraints is crucial for modern calibration studies of speleothems to instrumental climate records, which provides the basis for paleoclimate interpretations. However, a large unknown is how 14C is transferred from the atmosphere to any individual speleothem owing to the site-specific residence times of organic matter above cave systems. Here, we employ the bomb-pulse method to build chronologies from 14C measurements in combination with a new unsaturated zone C model which considers C decomposition as a continuum, to better understand unsaturated zone 14C dynamics. The bomb-pulse curves of eight speleothems from southern Australia in three contrasting climatic regions; the semi-arid Wellington Caves site, the mediterranean Golgotha Cave site and the montane Yarrangobilly Caves site, are investigated. Overall, the modelled 14C bomb-pulse curves produce excellent fits with measured 14C speleothem data (r2 = 0.82-0.99). The C modelling reveals that unsaturated zone C is predominately young at the semi-arid site, with a weighted-mean residence time of 32 years and that tree root respiration is likely an important source of vadose CO2. At the montane site, ∼39\% of C is young (<1 years), but the weighted-mean C ages are older (145-220 years). The mediterranean site has very little contribution from young C (<12\%: 0-1 years), with weighted-mean ages between 157 and 245 years, likely due to greater adsorption of organic matter in the upper vadose zone during matrix flow, and remobilisation of C from young syngenetic karst. New end members for low speleothem Dead Carbon Proportion (DCP) are identified (2.19\% and 1.65\%, respectively) for Australian montane and semi-arid zone speleothems, where oversupply of modern CO2 in the vadose zone leads to lower DCP. It is also demonstrated that DCP can be quite variable over small time scales, that processes may be difficult to untangle and a constant DCP assumption is likely invalid. DCP variability over time is mainly controlled by the changes vadose zone CO2, where vegetation regeneration, wild-fires and karst hydrology play an important role.}, language = {en} } @article{AichnerMakhmudovRajabovetal.2019, author = {Aichner, Bernhard and Makhmudov, Zafar and Rajabov, Iljomjon and Zhang, Qiong and Pausata, Francesco Salvatore R. and Werner, Martin and Heinecke, Liv and Kuessner, Marie L. and Feakins, Sarah J. and Sachse, Dirk and Mischke, Steffen}, title = {Hydroclimate in the Pamirs Was Driven by Changes in Precipitation-Evaporation Seasonality Since theLast Glacial Period}, series = {Geophysical research letters}, volume = {46}, journal = {Geophysical research letters}, number = {23}, publisher = {American Geophysical Union}, address = {Washington}, issn = {0094-8276}, doi = {10.1029/2019GL085202}, pages = {13972 -- 13983}, year = {2019}, abstract = {The Central Asian Pamir Mountains (Pamirs) are a high-altitude region sensitive to climatic change, with only few paleoclimatic records available. To examine the glacial-interglacial hydrological changes in the region, we analyzed the geochemical parameters of a 31-kyr record from Lake Karakul and performed a set of experiments with climate models to interpret the results. delta D values of terrestrial biomarkers showed insolation-driven trends reflecting major shifts of water vapor sources. For aquatic biomarkers, positive delta D shifts driven by changes in precipitation seasonality were observed at ca. 31-30, 28-26, and 17-14 kyr BP. Multiproxy paleoecological data and modelling results suggest that increased water availability, induced by decreased summer evaporation, triggered higher lake levels during those episodes, possibly synchronous to northern hemispheric rapid climate events. We conclude that seasonal changes in precipitation-evaporation balance significantly influenced the hydrological state of a large waterbody such as Lake Karakul, while annual precipitation amount and inflows remained fairly constant.}, language = {en} } @phdthesis{Senftleben2020, author = {Senftleben, Robin}, title = {Earth's magnetic field over the last 1000 years}, doi = {10.25932/publishup-47315}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-473150}, school = {Universit{\"a}t Potsdam}, pages = {xii, 104}, year = {2020}, abstract = {To investigate the reliability and stability of spherical harmonic models based on archeo/-paleomagnetic data, 2000 Geomagnetic models were calculated. All models are based on the same data set but with randomized uncertainties. Comparison of these models to the geomagnetic field model gufm1 showed that large scale magnetic field structures up to spherical harmonic degree 4 are stable throughout all models. Through a ranking of all models by comparing the dipole coefficients to gufm1 more realistic uncertainty estimates were derived than the authors of the data provide. The derived uncertainty estimates were used in further modelling, which combines archeo/-paleomagnetic and historical data. The huge difference in data count, accuracy and coverage of these two very different data sources made it necessary to introduce a time dependent spatial damping, which was constructed to constrain the spatial complexity of the model. Finally 501 models were calculated by considering that each data point is a Gaussian random variable, whose mean is the original value and whose standard deviation is its uncertainty. The final model arhimag1k is calculated by taking the mean of the 501 sets of Gauss coefficients. arhimag1k fits different dependent and independent data sets well. It shows an early reverse flux patch at the core-mantle boundary between 1000 AD and 1200 AD at the location of the South Atlantic Anomaly today. Another interesting feature is a high latitude flux patch over Greenland between 1200 and 1400 AD. The dipole moment shows a constant behaviour between 1600 and 1840 AD. In the second part of the thesis 4 new paleointensities from 4 different flows of the island Fogo, which is part of Cape Verde, are presented. The data is fitted well by arhimag1k with the exception of the value at 1663 of 28.3 microtesla, which is approximately 10 microtesla lower than the model suggest.}, language = {en} } @phdthesis{Brune2018, author = {Brune, Sascha}, title = {Modelling continental rift dynamics}, doi = {10.25932/publishup-43236}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-432364}, school = {Universit{\"a}t Potsdam}, pages = {192}, year = {2018}, abstract = {Continental rift systems open up unique possibilities to study the geodynamic system of our planet: geodynamic localization processes are imprinted in the morphology of the rift by governing the time-dependent activity of faults, the topographic evolution of the rift or by controlling whether a rift is symmetric or asymmetric. Since lithospheric necking localizes strain towards the rift centre, deformation structures of previous rift phases are often well preserved and passive margins, the end product of continental rifting, retain key information about the tectonic history from rift inception to continental rupture. Current understanding of continental rift evolution is based on combining observations from active rifts with data collected at rifted margins. Connecting these isolated data sets is often accomplished in a conceptual way and leaves room for subjective interpretation. Geodynamic forward models, however, have the potential to link individual data sets in a quantitative manner, using additional constraints from rock mechanics and rheology, which allows to transcend previous conceptual models of rift evolution. By quantifying geodynamic processes within continental rifts, numerical modelling allows key insight to tectonic processes that operate also in other plate boundary settings, such as mid ocean ridges, collisional mountain chains or subduction zones. In this thesis, I combine numerical, plate-tectonic, analytical, and analogue modelling approaches, whereas numerical thermomechanical modelling constitutes the primary tool. This method advanced rapidly during the last two decades owing to dedicated software development and the availability of massively parallel computer facilities. Nevertheless, only recently the geodynamical modelling community was able to capture 3D lithospheric-scale rift dynamics from onset of extension to final continental rupture. The first chapter of this thesis provides a broad introduction to continental rifting, a summary of the applied rift modelling methods and a short overview of previews studies. The following chapters, which constitute the main part of this thesis feature studies on plate boundary dynamics in two and three dimension followed by global scale analyses (Fig. 1). Chapter II focuses on 2D geodynamic modelling of rifted margin formation. It highlights the formation of wide areas of hyperextended crustal slivers via rift migration as a key process that affected many rifted margins worldwide. This chapter also contains a study of rift velocity evolution, showing that rift strength loss and extension velocity are linked through a dynamic feed-back. This process results in abrupt accelerations of the involved plates during rifting illustrating for the first time that rift dynamics plays a role in changing global-scale plate motions. Since rift velocity affects key processes like faulting, melting and lower crustal flow, this study also implies that the slow-fast velocity evolution should be imprinted in rifted margin structures. Chapter III relies on 3D Cartesian rift models in order to investigate various aspects of rift obliquity. Oblique rifting occurs if the extension direction is not orthogonal to the rift trend. Using 3D lithospheric-scale models from rift initialisation to breakup I could isolate a characteristic evolution of dominant fault orientations. Further work in Chapter III addresses the impact of rift obliquity on the strength of the rift system. We illustrate that oblique rifting is mechanically preferred over orthogonal rifting, because the brittle yielding requires a lower tectonic force. This mechanism elucidates rift competition during South Atlantic rifting, where the more oblique Equatorial Atlantic Rift proceeded to breakup while the simultaneously active but less oblique West African rift system became a failed rift. Finally this Chapter also investigates the impact of a previous rift phase on current tectonic activity in the linkage area of the Kenyan with Ethiopian rift. We show that the along strike changes in rift style are not caused by changes in crustal rheology. Instead the rift linkage pattern in this area can be explained when accounting for the thinned crust and lithosphere of a Mesozoic rift event. Chapter IV investigates rifting from the global perspective. A first study extends the oblique rift topic of the previous chapter to global scale by investigating the frequency of oblique rifting during the last 230 million years. We find that approximately 70\% of all ocean-forming rift segments involved an oblique component of extension where obliquities exceed 20°. This highlights the relevance of 3D approaches in modelling, surveying, and interpretation of many rifted margins. In a final study, we propose a link between continental rift activity, diffuse CO2 degassing and Mesozoic/Cenozoic climate changes. We used recent CO2 flux measurements in continental rifts to estimate worldwide rift-related CO2 release, which we based on the global extent of rifts through time. The first-order correlation to paleo-atmospheric CO2 proxy data suggests that rifts constitute a major element of the global carbon cycle.}, language = {en} } @article{vanSchaikBronstertdeJongetal.2014, author = {van Schaik, N. Loes M. B. and Bronstert, Axel and de Jong, S. M. and Jetten, V. G. and van Dam, J. C. and Ritsema, C. J. and Schnabel, Susanne}, title = {Process-based modelling of a headwater catchment in a semi-arid area: the influence of macropore flow}, series = {Hydrological processes}, volume = {28}, journal = {Hydrological processes}, number = {24}, publisher = {Wiley-Blackwell}, address = {Hoboken}, issn = {0885-6087}, doi = {10.1002/hyp.10086}, pages = {5805 -- 5816}, year = {2014}, abstract = {Subsurface stormflow is thought to occur mainly in humid environments with steep terrains. However, in semi-arid areas, preferential flow through macropores can also result in a significant contribution of subsurface stormflow to catchment runoff for varying catchment conditions. Most hydrological models neglect this important subsurface preferential flow. Here, we use the process-oriented hydrological model Hillflow-3D, which includes a macropore flow approach, to simulate rainfall-runoff in the semi-arid Parapunos catchment in Spain, where macropore flow was observed in previous research. The model was extended for this study to account for sorptivity under very dry soil conditions. The results of the model simulations with and without macropore flow are compared. Both model versions give reasonable results for average rainfall situations, although the approach with the macropore concept provides slightly better results. The model results for scenarios of extreme rainfall events (>13.3mm30min(-1)) however show large differences between the versions with and without macropores. These model results compared with measured rainfall-runoff data show that the model with the macropore concept is better. Our conclusion is that preferential flow is important in controlling surface runoff in case of specific, high intensity rainfall events. Therefore, preferential flow processes must be included in hydrological models where we know that preferential flow occurs. Hydrological process models with a less detailed process description may fit observed average events reasonably well but can result in erroneous predictions for more extreme events. Copyright (c) 2013 John Wiley \& Sons, Ltd.}, language = {en} } @article{MoradiCarminatiVetterleinetal.2011, author = {Moradi, Ahmad B. and Carminati, Andrea and Vetterlein, Doris and Vontobel, Peter and Lehmann, Eberhard and Weller, Ulrich and Hopmans, Jan W. and Vogel, Hans-J{\"o}rg and Oswald, Sascha Eric}, title = {Three-dimensional visualization and quantification of water content in the rhizosphere}, series = {New phytologist : international journal of plant science}, volume = {192}, journal = {New phytologist : international journal of plant science}, number = {3}, publisher = {Wiley-Blackwell}, address = {Hoboken}, issn = {0028-646X}, doi = {10.1111/j.1469-8137.2011.03826.x}, pages = {653 -- 663}, year = {2011}, abstract = {Despite the importance of rhizosphere properties for water flow from soil to roots, there is limited quantitative information on the distribution of water in the rhizosphere of plants. Here, we used neutron tomography to quantify and visualize the water content in the rhizosphere of the plant species chickpea (Cicer arietinum), white lupin (Lupinus albus), and maize (Zea mays) 12 d after planting. We clearly observed increasing soil water contents (h) towards the root surface for all three plant species, as opposed to the usual assumption of decreasing water content. This was true for tap roots and lateral roots of both upper and lower parts of the root system. Furthermore, water gradients around the lower part of the roots were smaller and extended further into bulk soil compared with the upper part, where the gradients in water content were steeper. Incorporating the hydraulic conductivity and water retention parameters of the rhizosphere into our model, we could simulate the gradual changes of h towards the root surface, in agreement with the observations. The modelling result suggests that roots in their rhizosphere may modify the hydraulic properties of soil in a way that improves uptake under dry conditions.}, language = {en} } @phdthesis{Schroeder2015, author = {Schr{\"o}der, Sarah}, title = {Modelling surface evolution coupled with tectonics}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-90385}, school = {Universit{\"a}t Potsdam}, pages = {viii, 129}, year = {2015}, abstract = {This study presents the development of 1D and 2D Surface Evolution Codes (SECs) and their coupling to any lithospheric-scale (thermo-)mechanical code with a quadrilateral structured surface mesh. Both SECs involve diffusion as approach for hillslope processes and the stream power law to reflect riverbed incision. The 1D SEC settles sediment that was produced by fluvial incision in the appropriate minimum, while the supply-limited 2D SEC DANSER uses a fast filling algorithm to model sedimantation. It is based on a cellular automaton. A slope-dependent factor in the sediment flux extends the diffusion equation to nonlinear diffusion. The discharge accumulation is achieved with the D8-algorithm and an improved drainage accumulation routine. Lateral incision enhances the incision's modelling. Following empirical laws, it incises channels of several cells width. The coupling method enables different temporal and spatial resolutions of the SEC and the thermo-mechanical code. It transfers vertical as well as horizontal displacements to the surface model. A weighted smoothing of the 3D surface displacements is implemented. The smoothed displacement vectors transmit the deformation by bilinear interpolation to the surface model. These interpolation methods ensure mass conservation in both directions and prevent the two surfaces from drifting apart. The presented applications refer to the evolution of the Pamir orogen. A calibration of DANSER's parameters with geomorphological data and a DEM as initial topography highlights the advantage of lateral incision. Preserving the channel width and reflecting incision peaks in narrow channels, this closes the huge gap between current orogen-scale incision models and observed topographies. River capturing models in a system of fault-bounded block rotations reaffirm the importance of the lateral incision routine for capturing events with channel initiation. The models show a low probability of river capturings with large deflection angles. While the probability of river capturing is directly depending on the uplift rate, the erodibility inside of a dip-slip fault speeds up headward erosion along the fault: The model's capturing speed increases within a fault. Coupling DANSER with the thermo-mechanical code SLIM 3D emphasizes the versatility of the SEC. While DANSER has minor influence on the lithospheric evolution of an indenter model, the brittle surface deformation is strongly affected by its sedimentation, widening a basin in between two forming orogens and also the southern part of the southern orogen to south, east and west.}, language = {en} }