@misc{GholamrezaieScheckWenderothBottetal.2019, author = {Gholamrezaie, Ershad and Scheck-Wenderoth, Magdalena and Bott, Judith and Heidbach, Oliver and Strecker, Manfred}, title = {3-D crustal density model of the Sea of Marmara}, series = {Postprints der Universit{\"a}t Potsdam Mathematisch-Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam Mathematisch-Naturwissenschaftliche Reihe}, number = {737}, issn = {1866-8372}, doi = {10.25932/publishup-43466}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-434661}, pages = {785 -- 807}, year = {2019}, abstract = {Abstract. The Sea of Marmara, in northwestern Turkey, is a transition zone where the dextral North Anatolian Fault zone (NAFZ) propagates westward from the Anatolian Plate to the Aegean Sea Plate. The area is of interest in the context of seismic hazard of Istanbul, a metropolitan area with about 15 million inhabitants. Geophysical observations indicate that the crust is heterogeneous beneath the Marmara basin, but a detailed characterization of the crustal heterogeneities is still missing. To assess if and how crustal heterogeneities are related to the NAFZ segmentation below the Sea of Marmara, we develop new crustal-scale 3-D density models which integrate geological and seismological data and that are additionally constrained by 3-D gravity modeling. For the latter, we use two different gravity datasets including global satellite data and local marine gravity observation. Considering the two different datasets and the general non-uniqueness in potential field modeling, we suggest three possible "end-member" solutions that are all consistent with the observed gravity field and illustrate the spectrum of possible solutions. These models indicate that the observed gravitational anomalies originate from significant density heterogeneities within the crust. Two layers of sediments, one syn-kinematic and one pre-kinematic with respect to the Sea of Marmara formation are underlain by a heterogeneous crystalline crust. A felsic upper crystalline crust (average density of 2720 kgm⁻³) and an intermediate to mafic lower crystalline crust (average density of 2890 kgm⁻³) appear to be cross-cut by two large, dome-shaped mafic highdensity bodies (density of 2890 to 3150 kgm⁻³) of considerable thickness above a rather uniform lithospheric mantle (3300 kgm⁻³). The spatial correlation between two major bends of the main Marmara fault and the location of the highdensity bodies suggests that the distribution of lithological heterogeneities within the crust controls the rheological behavior along the NAFZ and, consequently, maybe influences fault segmentation and thus the seismic hazard assessment in the region.}, language = {en} } @article{GholamrezaieScheckWenderothBottetal.2019, author = {Gholamrezaie, Ershad and Scheck-Wenderoth, Magdalena and Bott, Judith and Heidbach, Oliver and Strecker, Manfred}, title = {3-D crustal density model of the Sea of Marmara}, series = {Solid Earth}, volume = {10}, journal = {Solid Earth}, publisher = {Copernicus Publ.}, address = {G{\"o}ttingen}, issn = {1869-9510}, doi = {10.5194/se-10-785-2019}, pages = {785 -- 807}, year = {2019}, abstract = {Abstract. The Sea of Marmara, in northwestern Turkey, is a transition zone where the dextral North Anatolian Fault zone (NAFZ) propagates westward from the Anatolian Plate to the Aegean Sea Plate. The area is of interest in the context of seismic hazard of Istanbul, a metropolitan area with about 15 million inhabitants. Geophysical observations indicate that the crust is heterogeneous beneath the Marmara basin, but a detailed characterization of the crustal heterogeneities is still missing. To assess if and how crustal heterogeneities are related to the NAFZ segmentation below the Sea of Marmara, we develop new crustal-scale 3-D density models which integrate geological and seismological data and that are additionally constrained by 3-D gravity modeling. For the latter, we use two different gravity datasets including global satellite data and local marine gravity observation. Considering the two different datasets and the general non-uniqueness in potential field modeling, we suggest three possible "end-member" solutions that are all consistent with the observed gravity field and illustrate the spectrum of possible solutions. These models indicate that the observed gravitational anomalies originate from significant density heterogeneities within the crust. Two layers of sediments, one syn-kinematic and one pre-kinematic with respect to the Sea of Marmara formation are underlain by a heterogeneous crystalline crust. A felsic upper crystalline crust (average density of 2720 kgm⁻³) and an intermediate to mafic lower crystalline crust (average density of 2890 kgm⁻³) appear to be cross-cut by two large, dome-shaped mafic highdensity bodies (density of 2890 to 3150 kgm⁻³) of considerable thickness above a rather uniform lithospheric mantle (3300 kgm⁻³). The spatial correlation between two major bends of the main Marmara fault and the location of the highdensity bodies suggests that the distribution of lithological heterogeneities within the crust controls the rheological behavior along the NAFZ and, consequently, maybe influences fault segmentation and thus the seismic hazard assessment in the region.}, language = {en} } @article{GomezGarciaMeessenScheckWenderothetal.2019, author = {Gomez-Garcia, Angela Maria and Meeßen, Christian and Scheck-Wenderoth, Magdalena and Monsalve, Gaspar and Bott, Judith and Bernhardt, Anne and Bernal, Gladys}, title = {3-D Modeling of Vertical Gravity Gradients and the Delimitation of Tectonic Boundaries: The Caribbean Oceanic Domain as a Case Study}, series = {Geochemistry, geophysics, geosystems}, volume = {20}, journal = {Geochemistry, geophysics, geosystems}, number = {11}, publisher = {American Geophysical Union}, address = {Washington}, issn = {1525-2027}, doi = {10.1029/2019GC008340}, pages = {5371 -- 5393}, year = {2019}, abstract = {Geophysical data acquisition in oceanic domains is challenging, implying measurements with low and/or nonhomogeneous spatial resolution. The evolution of satellite gravimetry and altimetry techniques allows testing 3-D density models of the lithosphere, taking advantage of the high spatial resolution and homogeneous coverage of satellites. However, it is not trivial to discretise the source of the gravity field at different depths. Here, we propose a new method for inferring tectonic boundaries at the crustal level. As a novelty, instead of modeling the gravity anomalies and assuming a flat Earth approximation, we model the vertical gravity gradients (VGG) in spherical coordinates, which are especially sensitive to density contrasts in the upper layers of the Earth. To validate the methodology, the complex oceanic domain of the Caribbean region is studied, which includes different crustal domains with a tectonic history since Late Jurassic time. After defining a lithospheric starting model constrained by up-to-date geophysical data sets, we tested several a-priory density distributions and selected the model with the minimum misfits with respect to the VGG calculated from the EIGEN-6C4 data set. Additionally, the density of the crystalline crust was inferred by inverting the VGG field. Our methodology enabled us not only to refine, confirm, and/or propose tectonic boundaries in the study area but also to identify a new anomalous buoyant body, located in the South Lesser Antilles subduction zone, and high-density bodies along the Greater, Lesser, and Leeward Antilles forearcs.}, language = {en} } @article{KaiserCacaceScheckWenderoth2013, author = {Kaiser, Bj{\"o}rn Onno and Cacace, Mauro and Scheck-Wenderoth, Magdalena}, title = {3D coupled fluid and heat transport simulations of the Northeast German Basin and their sensitivity to the spatial discretization - different sensitivities for different mechanisms of heat transport}, series = {Environmental earth sciences}, volume = {70}, journal = {Environmental earth sciences}, number = {8}, publisher = {Springer}, address = {New York}, issn = {1866-6280}, doi = {10.1007/s12665-013-2249-7}, pages = {3643 -- 3659}, year = {2013}, abstract = {Based on a numerical model of the Northeast German Basin (NEGB), we investigate the sensitivity of the calculated thermal field as resulting from heat conduction, forced and free convection in response to consecutive horizontal and vertical mesh refinements. Our results suggest that computational findings are more sensitive to consecutive horizontal mesh refinements than to changes in the vertical resolution. In addition, the degree of mesh sensitivity depends strongly on the type of the process being investigated, whether heat conduction, forced convection or free thermal convection represents the active heat driver. In this regard, heat conduction exhibits to be relative robust to imposed changes in the spatial discretization. A systematic mesh sensitivity is observed in areas where forced convection promotes an effective role in shorten the background conductive thermal field. In contrast, free thermal convection is to be regarded as the most sensitive heat transport process as demonstrated by non-systematic changes in the temperature field with respect to imposed changes in the model resolution.}, language = {en} } @article{IbarraLiuMeessenetal.2019, author = {Ibarra, Federico and Liu, Sibiao and Meeßen, Christian and Prezzi, Claudia Beatriz and Bott, Judith and Scheck-Wenderoth, Magdalena and Sobolev, Stephan Vladimir and Strecker, Manfred}, title = {3D data-derived lithospheric structure of the Central Andes and its implications for deformation: Insights from gravity and geodynamic modelling}, series = {Tectonophysics : international journal of geotectonics and the geology and physics of the interior of the earth}, volume = {766}, journal = {Tectonophysics : international journal of geotectonics and the geology and physics of the interior of the earth}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0040-1951}, doi = {10.1016/j.tecto.2019.06.025}, pages = {453 -- 468}, year = {2019}, abstract = {We present a new three-dimensional density model of the Central Andes characterizing the structure and composition of the lithosphere together with a geodynamic simulation subjected to continental intraplate shortening. The principal aim of this study is to assess the link between heterogeneities in the lithosphere and different deformation patterns and styles along the orogen-foreland system of the Central Andes. First, we performed a 3D integration of new geological and geophysical data with previous models through forward modelling of Bouguer anomalies. Subsequently, a geodynamic model was set-up and parametrized from the previously obtained 3D structure and composition. We do not find a unambigous correlation between the resulting density configuration and terrane boundaries proposed by other authors. Our models reproduce the observed Bouguer anomaly and deformation patterns in the foreland. We find that thin-skinned deformation in the Subandean fold-and thrust belt is controlled by a thick sedimentary layer and coeval underthrusting of thin crust of the foreland beneath the thick crust of the Andean Plateau. In the adjacent thick-skinned deformation province of the inverted Cretaceous extensional Santa Barbara System sedimentary strata are much thinner and crustal thickness transitions from greater values in the Andean to a more reduced thickness in the foreland. Our results show that deformation processes occur where the highest gradients of lithospheric strength are present between the orogen and the foreland, thus suggesting a spatial correlation between deformation and lithospheric strength.}, language = {en} } @misc{ZieglerHeidbachReineckeretal.2016, author = {Ziegler, Moritz O. and Heidbach, Oliver and Reinecker, John and Przybycin, Anna M. and Scheck-Wenderoth, Magdalena}, title = {A multi-stage 3-D stress field modelling approach exemplified in the Bavarian Molasse Basin}, series = {Postprints der Universit{\"a}t Potsdam Mathematisch-Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam Mathematisch-Naturwissenschaftliche Reihe}, number = {556}, issn = {1866-8372}, doi = {10.25932/publishup-40980}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-409806}, pages = {18}, year = {2016}, abstract = {The knowledge of the contemporary in situ stress state is a key issue for safe and sustainable subsurface engineering. However, information on the orientation and magnitudes of the stress state is limited and often not available for the areas of interest. Therefore 3-D geomechanical-numerical modelling is used to estimate the in situ stress state and the distance of faults from failure for application in subsurface engineering. The main challenge in this approach is to bridge the gap in scale between the widely scattered data used for calibration of the model and the high resolution in the target area required for the application. We present a multi-stage 3-D geomechanical-numerical approach which provides a state-of-the-art model of the stress field for a reservoir-scale area from widely scattered data records. Therefore, we first use a large-scale regional model which is calibrated by available stress data and provides the full 3-D stress tensor at discrete points in the entire model volume. The modelled stress state is used subsequently for the calibration of a smaller-scale model located within the large-scale model in an area without any observed stress data records. We exemplify this approach with two-stages for the area around Munich in the German Molasse Basin. As an example of application, we estimate the scalar values for slip tendency and fracture potential from the model results as measures for the criticality of fault reactivation in the reservoir-scale model. The modelling results show that variations due to uncertainties in the input data are mainly introduced by the uncertain material properties and missing S-Hmax magnitude estimates needed for a more reliable model calibration. This leads to the conclusion that at this stage the model's reliability depends only on the amount and quality of available stress information rather than on the modelling technique itself or on local details of the model geometry. Any improvements in modelling and increases in model reliability can only be achieved using more high-quality data for calibration.}, language = {en} } @article{ZieglerHeidbachReineckeretal.2016, author = {Ziegler, Moritz O. and Heidbach, Oliver and Reinecker, John and Przybycin, Anna M. and Scheck-Wenderoth, Magdalena}, title = {A multi-stage 3-D stress field modelling approach exemplified in the Bavarian Molasse Basin}, series = {Solid earth}, volume = {7}, journal = {Solid earth}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {1869-9510}, doi = {10.5194/se-7-1365-2016}, pages = {1365 -- 1382}, year = {2016}, language = {en} } @article{KaiserCacaceScheckWenderothetal.2011, author = {Kaiser, Bjoern Onno and Cacace, Mauro and Scheck-Wenderoth, Magdalena and Lewerenz, Bjoern}, title = {Characterization of main heat transport processes in the Northeast German Basin constraints from 3-D numerical models}, series = {Geochemistry, geophysics, geosystems}, volume = {12}, journal = {Geochemistry, geophysics, geosystems}, number = {13}, publisher = {American Geophysical Union}, address = {Washington}, issn = {1525-2027}, doi = {10.1029/2011GC003535}, pages = {17}, year = {2011}, abstract = {To investigate and quantify main physical heat driving processes affecting the present-day subsurface thermal field, we study a complex geological setting, the Northeast German Basin (NEGB). The internal geological structure of the NEGB is characterized by the presence of a relatively thick layer of Permian Zechstein salt (up to 5000 m), which forms many salt diapirs and pillows locally reaching nearly the surface. By means of three-dimensional numerical simulations we explore the role of heat conduction, pressure, and density driven groundwater flow as well as fluid viscosity related effects. Our results suggest that the regional temperature distribution within the basin results from interactions between regional pressure forces as driven by topographic gradients and thermal diffusion locally enhanced by thermal conductivity contrasts between the different sedimentary rocks with the highly conductive salt playing a prominent role. In contrast, buoyancy forces triggered by temperature-dependent fluid density variations are demonstrated to affect only locally the internal thermal configuration. Locations, geometry, and wavelengths of convective thermal anomalies are mainly controlled by the permeability field and thickness values of the respective geological layers.}, language = {en} } @article{RodriguezPicedaScheckWenderothBottetal.2022, author = {Rodriguez Piceda, Constanza and Scheck-Wenderoth, Magdalena and Bott, Judith and Gomez Dacal, Maria Laura and Cacace, Mauro and Pons, Michael and Prezzi, Claudia and Strecker, Manfred}, title = {Controls of the Lithospheric Thermal Field of an Ocean-Continent Subduction Zone}, series = {Lithosphere / Geological Society of America}, volume = {2022}, journal = {Lithosphere / Geological Society of America}, number = {1}, publisher = {GeoScienceWorld}, address = {McLean}, issn = {1941-8264}, doi = {10.2113/2022/2237272}, pages = {26}, year = {2022}, abstract = {In an ocean-continent subduction zone, the assessment of the lithospheric thermal state is essential to determine the controls of the deformation within the upper plate and the dip angle of the subducting lithosphere. In this study, we evaluate the degree of influence of both the configuration of the upper plate (i.e., thickness and composition of the rock units) and variations of the subduction angle on the lithospheric thermal field of the southern Central Andes (29 degrees-39 degrees S). Here, the subduction angle increases from subhorizontal (5 degrees) north of 33 degrees S to steep (similar to 30 degrees) in the south. We derived the 3D temperature and heat flow distribution of the lithosphere in the southern Central Andes considering conversion of S wave tomography to temperatures together with steady-state conductive thermal modeling. We found that the orogen is overall warmer than the forearc and the foreland and that the lithosphere of the northern part of the foreland appears colder than its southern counterpart. Sedimentary blanketing and the thickness of the radiogenic crust exert the main control on the shallow thermal field (<50km depth). Specific conditions are present where the oceanic slab is relatively shallow (<85 km depth) and the radiogenic crust is thin. This configuration results in relatively colder temperatures compared to regions where the radiogenic crust is thick and the slab is steep. At depths >50km, the temperatures of the overriding plate are mainly controlled by the mantle heat input and the subduction angle. The thermal field of the upper plate likely preserves the flat subduction angle and influences the spatial distribution of shortening.}, language = {en} } @article{CherubiniCacaceScheckWenderothetal.2013, author = {Cherubini, Yvonne and Cacace, Mauro and Scheck-Wenderoth, Magdalena and Moeck, Inga and Lewerenz, Bj{\"o}rn}, title = {Controls on the deep thermal field - implications from 3-D numerical simulations for the geothermal research site Groß Sch{\"o}nebeck}, series = {Environmental earth sciences}, volume = {70}, journal = {Environmental earth sciences}, number = {8}, publisher = {Springer}, address = {New York}, issn = {1866-6280}, doi = {10.1007/s12665-013-2519-4}, pages = {3619 -- 3642}, year = {2013}, abstract = {The deep thermal field in sedimentary basins can be affected by convection, conduction or both resulting from the structural inventory, physical properties of geological layers and physical processes taking place therein. For geothermal energy extraction, the controlling factors of the deep thermal field need to be understood to delineate favorable drill sites and exploitation compartments. We use geologically based 3-D finite element simulations to figure out the geologic controls on the thermal field of the geothermal research site Gro Schonebeck located in the E part of the North German Basin. Its target reservoir consists of Permian Rotliegend clastics that compose the lower part of a succession of Late Carboniferous to Cenozoic sediments, subdivided into several aquifers and aquicludes. The sedimentary succession includes a layer of mobilized Upper Permian Zechstein salt which plays a special role for the thermal field due to its high thermal conductivity. Furthermore, the salt is impermeable and due to its rheology decouples the fault systems in the suprasalt units from subsalt layers. Conductive and coupled fluid and heat transport simulations are carried out to assess the relative impact of different heat transfer mechanisms on the temperature distribution. The measured temperatures in 7 wells are used for model validation and show a better fit with models considering fluid and heat transport than with a purely conductive model. Our results suggest that advective and convective heat transport are important heat transfer processes in the suprasalt sediments. In contrast, thermal conduction mainly controls the subsalt layers. With a third simulation, we investigate the influence of a major permeable and of three impermeable faults dissecting the subsalt target reservoir and compare the results to the coupled model where no faults are integrated. The permeable fault may have a local, strong impact on the thermal, pressure and velocity fields whereas the impermeable faults only cause deviations of the pressure field.}, language = {en} } @article{MeessenSippelScheckWenderothetal.2018, author = {Meessen, Christian and Sippel, Judith and Scheck-Wenderoth, Magdalena and Heine, C. and Strecker, Manfred}, title = {Crustal structure of the andean foreland in Northern Argentina}, series = {Journal of geophysical research : Solid earth}, volume = {123}, journal = {Journal of geophysical research : Solid earth}, number = {2}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2169-9313}, doi = {10.1002/2017JB014296}, pages = {1875 -- 1903}, year = {2018}, abstract = {Previous thermomechanical modeling studies indicated that variations in the temperature and strength of the crystalline crust might be responsible for the juxtaposition of domains with thin-skinned and thick-skinned crustal deformation along strike the foreland of the central Andes. However, there is no evidence supporting this hypothesis from data-integrative models. We aim to derive the density structure of the lithosphere by means of integrated 3-D density modeling, in order to provide a new basis for discussions of compositional variations within the crust and for future thermal and rheological modeling studies. Therefore, we utilize available geological and geophysical data to obtain a structural and density model of the uppermost 200km of the Earth. The derived model is consistent with the observed Bouguer gravity field. Our results indicate that the crystalline crust in northern Argentina can be represented by a lighter upper crust (2,800kg/m(3)) and a denser lower crust (3,100kg/m(3)). We find new evidence for high bulk crustal densities >3,000kg/m(3) in the northern Pampia terrane. These could originate from subducted Puncoviscana wackes or pelites that ponded to the base of the crystalline crust in the late Proterozoic or indicate increasing bulk content of mafic material. The precise composition of the northern foreland crust, whether mafic or felsic, has significant implications for further thermomechanical models and the rheological behavior of the lithosphere. A detailed sensitivity analysis of the input parameters indicates that the model results are robust with respect to the given uncertainties of the input data.}, language = {en} } @article{SpoonerScheckWenderothGoetzeetal.2019, author = {Spooner, Cameron and Scheck-Wenderoth, Magdalena and G{\"o}tze, Hans-J{\"u}rgen and Ebbing, J{\"o}rg and Hetenyi, Gyoergy}, title = {Density distribution across the Alpine lithosphere constrained by 3-D gravity modelling and relation to seismicity and deformation}, series = {Solid earth}, volume = {10}, journal = {Solid earth}, number = {6}, publisher = {Copernicus}, address = {G{\"o}ttingen}, organization = {AlpArray Working Grp}, issn = {1869-9510}, doi = {10.5194/se-10-2073-2019}, pages = {2073 -- 2088}, year = {2019}, abstract = {The Alpine orogen formed as a result of the collision between the Adriatic and European plates. Significant crustal heterogeneity exists within the region due to the long history of interplay between these plates, other continental and oceanic blocks in the region, and inherited crustal features from earlier orogenies. Deformation relating to the collision continues to the present day. Here, a seismically constrained, 3-D structural and density model of the lithosphere of the Alps and their respective forelands, derived from integrating numerous geoscientific datasets, was adjusted to match the observed gravity field. It is shown that the distribution of seismicity and deformation within the region correlates well to thickness and density changes within the crust, and that the present-day Adriatic crust is both thinner and denser (22.5 km, 2800 kg m(-3) ) than the European crust (27.5 km, 2750 kg m(-3)). Alpine crust derived from each respective plate is found to show the same trend, with zones of Adriatic provenance (Austro-Alpine unit and Southern Alps) found to be denser and those of European provenance (Helvetic zone and Tauern Window) to be less dense. This suggests that the respective plates and related terranes had similar crustal properties to the present-day ones prior to orogenesis. The model generated here is available for open-access use to further discussions about the crust in the region.}, language = {en} } @misc{KuehnKempkadeLuciaetal.2017, author = {K{\"u}hn, Michael and Kempka, Thomas and de Lucia, Marco and Scheck-Wenderoth, Magdalena}, title = {Dissolved CO2 storage in geological formations with low pressure, low risk and large capacities}, series = {Energy procedia}, volume = {114}, journal = {Energy procedia}, publisher = {Elsevier}, address = {Amsterdam}, issn = {1876-6102}, doi = {10.1016/j.egypro.2017.03.1607}, pages = {4722 -- 4727}, year = {2017}, abstract = {Geological CO2 storage is a mitigation technology to reduce CO2 emissions from fossil fuel combustion. However, major concerns are the pressure increase and saltwater displacement in the mainly targeted deep groundwater aquifers due to injection of supercritical CO2. The suggested solution is storage of CO2 exclusively in the dissolved state. In our exemplary regional case study of the North East German Basin based on a highly resolved temperature and pressure distribution model and a newly developed reactive transport coupling, we have quantified that 4.7 Gt of CO2 can be stored in solution compared to 1.5 Gt in the supercritical state.}, language = {en} } @unpublished{ScheckWenderothSchmeisserMuttietal.2013, author = {Scheck-Wenderoth, Magdalena and Schmeißer, Dieter and Mutti, Maria and Kolditz, Olaf and H{\"u}nges, Ernst and Schultz, Hans-Martin and Liebscher, Axel and Bock, Michaela}, title = {Geoenergy - new concepts for utilization of geo-reservoirs as potential energy sources}, series = {Environmental earth sciences}, volume = {70}, journal = {Environmental earth sciences}, number = {8}, publisher = {Springer}, address = {New York}, issn = {1866-6280}, doi = {10.1007/s12665-013-2877-y}, pages = {3427 -- 3431}, year = {2013}, language = {en} } @article{CacaceKaiserLewerenzetal.2010, author = {Cacace, Mauro and Kaiser, Bjoern Onno and Lewerenz, Bjoern and Scheck-Wenderoth, Magdalena}, title = {Geothermal energy in sedimentary basins : what we can learn from regional numerical models}, issn = {0009-2819}, doi = {10.1016/j.chemer.2010.05.017}, year = {2010}, abstract = {Understanding the interactions between the different processes that control the geothermal and fluid flow fields in sedimentary basins is crucial for exploitation of geothermal energy. Numerical models provide predictive and feasible information for a correct assessment of geothermal resources especially in areas where data acquisition is demanding. Here, we present results from numerical efforts to characterize the thermal structure and its interaction with the fluid system for the area of the North East German Basin (NEGB). The relative impact of the different (diffusive and advective) processes affecting the hydrothermal setting of the basin are investigated by means of three- dimensional numerical simulations. Lithospheric-scale numerical models are evaluated to understand the specific thermal signature of the relevant factors influencing the present-day conductive geothermal field in the NEGB. Shallow and deep structural controls on the thermal configuration of the basin are addressed and quantified. Interaction between the resulting thermal field and the active fluid system is investigated by means of three-dimensional simulations of coupled fluid flow and heat transport. Factors influencing stability and reliability of modeling predictions are discussed. The main effort is to build a physically consistent model for the basin which integrates the impacts of thermal gradients on the regional fluid regime and their coupling with the main geological units defining the basin.}, 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{DegenSpoonerScheckWenderothetal.2021, author = {Degen, Denise and Spooner, Cameron and Scheck-Wenderoth, Magdalena and Cacace, Mauro}, title = {How biased are our models?}, series = {Geoscientific model development : an interactive open access journal of the European Geosciences Union}, volume = {14}, journal = {Geoscientific model development : an interactive open access journal of the European Geosciences Union}, number = {11}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {1991-959X}, doi = {10.5194/gmd-14-7133-2021}, pages = {7133 -- 7153}, year = {2021}, abstract = {Geophysical process simulations play a crucial role in the understanding of the subsurface. This understanding is required to provide, for instance, clean energy sources such as geothermal energy. However, the calibration and validation of the physical models heavily rely on state measurements such as temperature. In this work, we demonstrate that focusing analyses purely on measurements introduces a high bias. This is illustrated through global sensitivity studies. The extensive exploration of the parameter space becomes feasible through the construction of suitable surrogate models via the reduced basis method, where the bias is found to result from very unequal data distribution. We propose schemes to compensate for parts of this bias. However, the bias cannot be entirely compensated. Therefore, we demonstrate the consequences of this bias with the example of a model calibration.}, language = {en} } @article{CherubiniCacaceBloecheretal.2013, author = {Cherubini, Yvonne and Cacace, Mauro and Bl{\"o}cher, Guido and Scheck-Wenderoth, Magdalena}, title = {Impact of single inclined faults on the fluid flow and heat transport - results from 3-D finite element simulations}, series = {Environmental earth sciences}, volume = {70}, journal = {Environmental earth sciences}, number = {8}, publisher = {Springer}, address = {New York}, issn = {1866-6280}, doi = {10.1007/s12665-012-2212-z}, pages = {3603 -- 3618}, year = {2013}, abstract = {The impact of inclined faults on the hydrothermal field is assessed by adding simplified structural settings to synthetic models. This study is innovative in carrying out numerical simulations because it integrates the real 3-D nature of flow influenced by a fault in a porous medium, thereby providing a useful tool for complex geothermal modelling. The 3-D simulations for the coupled fluid flow and heat transport processes are based on the finite element method. In the model, one geological layer is dissected by a dipping fault. Sensitivity analyses are conducted to quantify the effects of the fault's transmissivity on the fluid flow and thermal field. Different fault models are compared with a model where no fault is present to evaluate the effect of varying fault transmissivity. The results show that faults have a significant impact on the hydrothermal field. Varying either the fault zone width or the fault permeability will result in relevant differences in the pressure, velocity and temperature field. A linear relationship between fault zone width and fluid velocity is found, indicating that velocities increase with decreasing widths. The faults act as preferential pathways for advective heat transport in case of highly transmissive faults, whereas almost no fluid may be transported through poorly transmissive faults.}, language = {en} } @article{FreymarkBottCacaceetal.2019, author = {Freymark, Jessica and Bott, Judith and Cacace, Mauro and Ziegler, Moritz 0. and Scheck-Wenderoth, Magdalena}, title = {Influence of the Main Border Faults on the 3D Hydraulic Field of the Central Upper Rhine Graben}, series = {Geofluids}, journal = {Geofluids}, publisher = {Wiley-Hindawi}, address = {London}, issn = {1468-8115}, doi = {10.1155/2019/7520714}, pages = {21}, year = {2019}, abstract = {The Upper Rhine Graben (URG) is an active rift with a high geothermal potential. Despite being a well-studied area, the three-dimensional interaction of the main controlling factors of the thermal and hydraulic regime is still not fully understood. Therefore, we have used a data-based 3D structural model of the lithological configuration of the central URG for some conceptual numerical experiments of 3D coupled simulations of fluid and heat transport. To assess the influence of the main faults bordering the graben on the hydraulic and the deep thermal field, we carried out a sensitivity analysis on fault width and permeability. Depending on the assigned width and permeability of the main border faults, fluid velocity and temperatures are affected only in the direct proximity of the respective border faults. Hence, the hydraulic characteristics of these major faults do not significantly influence the graben-wide groundwater flow patterns. Instead, the different scenarios tested provide a consistent image of the main characteristics of fluid and heat transport as they have in common: (1) a topography-driven basin-wide fluid flow perpendicular to the rift axis from the graben shoulders to the rift center, (2) a N/NE-directed flow parallel to the rift axis in the center of the rift and, (3) a pronounced upflow of hot fluids along the rift central axis, where the streams from both sides of the rift merge. This upflow axis is predicted to occur predominantly in the center of the URG (northern and southern model area) and shifted towards the eastern boundary fault (central model area).}, language = {en} } @article{RodriguezPicedaScheckWenderothCacaceetal.2022, author = {Rodriguez Piceda, Constanza and Scheck-Wenderoth, Magdalena and Cacace, Mauro and Bott, Judith and Strecker, Manfred}, title = {Long-Term Lithospheric Strength and Upper-Plate Seismicity in the Southern Central Andes, 29 degrees-39 degrees S}, series = {Geochemistry, geophysics, geosystems}, volume = {23}, journal = {Geochemistry, geophysics, geosystems}, number = {3}, publisher = {American Geophysical Union}, address = {Washington}, issn = {1525-2027}, doi = {10.1029/2021GC010171}, pages = {22}, year = {2022}, abstract = {We examined the relationship between the mechanical strength of the lithosphere and the distribution of seismicity within the overriding continental plate of the southern Central Andes (SCA, 29 degrees-39 degrees S), where the oceanic Nazca Plate changes its subduction angle between 33 degrees S and 35 degrees S, from subhorizontal in the north (<5 degrees) to steep in the south (similar to 30 degrees). We computed the long-term lithospheric strength based on an existing 3D model describing variations in thickness, density, and temperature of the main geological units forming the lithosphere of the SCA and adjacent forearc and foreland regions. The comparison between our results and seismicity within the overriding plate (upper-plate seismicity) shows that most of the events occur within the modeled brittle domain of the lithosphere. The depth where the deformation mode switches from brittle frictional to thermally activated ductile creep provides a conservative lower bound to the seismogenic zone in the overriding plate of the study area. We also found that the majority of upper-plate earthquakes occurs within the realm of first-order contrasts in integrated strength (12.7-13.3 log Pam in the Andean orogen vs. 13.5-13.9 log Pam in the forearc and the foreland). Specific conditions characterize the mechanically strong northern foreland of the Andes, where seismicity is likely explained by the effects of slab steepening.}, language = {en} } @article{CacaceBloecherWatanabeetal.2013, author = {Cacace, Mauro and Bl{\"o}cher, Guido and Watanabe, Norihiro and M{\"o}ck, Inga and B{\"o}rsing, Nele and Scheck-Wenderoth, Magdalena and Kolditz, Olaf and H{\"u}nges, Ernst}, title = {Modelling of fractured carbonate reservoirs - outline of a novel technique via a case study from the Molasse Basin, southern Bavaria, Germany}, series = {Environmental earth sciences}, volume = {70}, journal = {Environmental earth sciences}, number = {8}, publisher = {Springer}, address = {New York}, issn = {1866-6280}, doi = {10.1007/s12665-013-2402-3}, pages = {3585 -- 3602}, year = {2013}, abstract = {Fluid flow in low-permeable carbonate rocks depends on the density of fractures, their interconnectivity and on the formation of fault damage zones. The present-day stress field influences the aperture hence the transmissivity of fractures whereas paleostress fields are responsible for the formation of faults and fractures. In low-permeable reservoir rocks, fault zones belong to the major targets. Before drilling, an estimate for reservoir productivity of wells drilled into the damage zone of faults is therefore required. Due to limitations in available data, a characterization of such reservoirs usually relies on the use of numerical techniques. The requirements of these mathematical models encompass a full integration of the actual fault geometry, comprising the dimension of the fault damage zone and of the fault core, and the individual population with properties of fault zones in the hanging and foot wall and the host rock. The paper presents both the technical approach to develop such a model and the property definition of heterogeneous fault zones and host rock with respect to the current stress field. The case study describes a deep geothermal reservoir in the western central Molasse Basin in southern Bavaria, Germany. Results from numerical simulations indicate that the well productivity can be enhanced along compressional fault zones if the interconnectivity of fractures is lateral caused by crossing synthetic and antithetic fractures. The model allows a deeper understanding of production tests and reservoir properties of faulted rocks.}, language = {en} } @article{ScheckWenderothCacaceMaystrenkoetal.2014, author = {Scheck-Wenderoth, Magdalena and Cacace, Mauro and Maystrenko, Yuriy Petrovich and Cherubini, Yvonne and Noack, Vera and Kaiser, Bjoern Onno and Sippel, Judith and Bjoern, Lewerenz}, title = {Models of heat transport in the Central European Basin System: Effective mechanisms at different scales}, series = {Marine and petroleum geology}, volume = {55}, journal = {Marine and petroleum geology}, publisher = {Elsevier}, address = {Oxford}, issn = {0264-8172}, doi = {10.1016/j.marpetgeo.2014.03.009}, pages = {315 -- 331}, year = {2014}, abstract = {Understanding heat transport in sedimentary basins requires an assessment of the regional 3D heat distribution and of the main physical mechanisms responsible for the transport of heat. We review results from different 3D numerical simulations of heat transport based on 3D basin models of the Central European Basin System (CEBS). Therefore we compare differently detailed 3D structural models of the area, previously published individually, to assess the influence of (1) different configurations of the deeper lithosphere, (2) the mechanism of heat transport considered and (3) large faults dissecting the sedimentary succession on the resulting thermal field and groundwater flow. Based on this comparison we propose a modelling strategy linking the regional and lithosphere-scale to the sub-basin and basin-fill scale and appropriately considering the effective heat transport processes. We find that conduction as the dominant mechanism of heat transport in sedimentary basins is controlled by the distribution of thermal conductivities, compositional and thickness variations of both the conductive and radiogenic crystalline crust as well as the insulating sediments and by variations in the depth to the thermal lithosphere-asthenosphere boundary. Variations of these factors cause thermal anomalies of specific wavelength and must be accounted for in regional thermal studies. In addition advective heat transport also exerts control on the thermal field on the regional scale. In contrast, convective heat transport and heat transport along faults is only locally important and needs to be considered for exploration on the reservoir scale. The general applicability of the proposed workflow makes it of interest for a broad range of application in geosciences including oil and gas exploration, geothermal utilization or carbon capture and sequestration issues. (C) 2014 Elsevier Ltd. All rights reserved.}, language = {en} } @article{KaiserCacaceScheckWenderoth2013, author = {Kaiser, Bj{\"o}rn Onno and Cacace, Mauro and Scheck-Wenderoth, Magdalena}, title = {Quaternary channels within the Northeast German Basin and their relevance on double diffusive convective transport processes - constraints from 3-D thermohaline numerical simulations}, series = {Geochemistry, geophysics, geosystems}, volume = {14}, journal = {Geochemistry, geophysics, geosystems}, number = {8}, publisher = {American Geophysical Union}, address = {Washington}, issn = {1525-2027}, doi = {10.1002/ggge.20192}, pages = {3156 -- 3175}, year = {2013}, abstract = {The internal geological structure of the Northeast German Basin (NEGB) is affected by intense salt diapirism and by the presence of several stratified aquifer complexes of regional relevance. The shallow Quaternary to late Tertiary freshwater aquifer is separated from the underlying Mesozoic saline aquifers by an embedded Tertiary clay enriched aquitard (Rupelian Aquitard). An important feature of this aquitard is that hydraulic connections between the upper and lower aquifers do exist in areas where the Rupelian Aquitard is missing (hydrogeological windows). Three-dimensional thermohaline numerical simulations are carried out to investigate the effects of such hydrogeological windows in the Rupelian Aquitard on the resulting groundwater, temperature, and salinity distributions. Numerical results suggest that hydrogeological windows act as preferential domains of hydraulic interconnectivity between the different aquifers at depth and enable vigorous heat and mass transport which causes a mixing of warm and saline groundwater with cold and less saline groundwater within both aquifers. In areas where the Rupelian Aquitard confines the Mesozoic aquifer, dissolved solutes from major salt structures are transported laterally giving rise to plumes of variable salinity content ranging from few hundreds of meters to several tens of kilometers. Furthermore, destabilizing thermal buoyancy forces may overwhelm counteracting stabilizing salinity induced forces offside of salt domes. This may result in buoyant upward groundwater flow transporting heat and mass to shallower levels within the same Mesozoic Aquifer.}, language = {en} } @article{NoackScheckWenderothCacace2012, author = {Noack, Vera and Scheck-Wenderoth, Magdalena and Cacace, Mauro}, title = {Sensitivity of 3D thermal models to the choice of boundary conditions and thermal properties: a case study for the area of Brandenburg (NE German Basin)}, series = {Environmental earth sciences}, volume = {67}, journal = {Environmental earth sciences}, number = {6}, publisher = {Springer}, address = {New York}, issn = {1866-6280}, doi = {10.1007/s12665-012-1614-2}, pages = {1695 -- 1711}, year = {2012}, abstract = {Based on newly available data of both, the structural setting and thermal properties, we compare 3D thermal models for the area of Brandenburg, located in the Northeast German Basin, to assess the sensitivity of our model results. The structural complexity of the basin fill is given by the configuration of the Zechstein salt with salt diapirs and salt pillows. This special configuration is very relevant for the thermal calculations because salt has a distinctly higher thermal conductivity than other sediments. We calculate the temperature using a FEMethod to solve the steady state heat conduction equation in 3D. Based on this approach, we evaluate the sensitivity of the steady-state conductive thermal field with respect to different lithospheric configurations and to the assigned thermal properties. We compare three different thermal models: (a) a crustal-scale model including a homogeneous crust, (b) a new lithosphere-scale model including a differentiated crust and (c) a crustal-scale model with a stepwise variation of measured thermal properties. The comparison with measured temperatures from different structural locations of the basin shows a good fit to the temperature predictions for the first two models, whereas the third model is distinctly colder. This indicates that effective thermal conductivities may be different from values determined by measurements on rock samples. The results suggest that conduction is the main heat transport mechanism in the Brandenburg area.}, language = {en} } @article{SpoonerScheckWenderothCacaceetal.2020, author = {Spooner, Cameron and Scheck-Wenderoth, Magdalena and Cacace, Mauro and G{\"o}tze, Hans-J{\"u}rgen and Luijendijk, Elco}, title = {The 3D thermal field across the Alpine orogen and its forelands and the relation to seismicity}, series = {Global and planetary change}, volume = {193}, journal = {Global and planetary change}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0921-8181}, doi = {10.1016/j.gloplacha.2020.103288}, pages = {14}, year = {2020}, abstract = {Temperature exerts a first order control on rock strength, principally via thermally activated creep deformation and on the distribution at depth of the brittle-ductile transition zone. The latter can be regarded as the lower bound to the seismogenic zone, thereby controlling the spatial distribution of seismicity within a lithospheric plate. As such, models of the crustal thermal field are important to understand the localisation of seismicity. Here we relate results from 3D simulations of the steady state thermal field of the Alpine orogen and its forelands to the distribution of seismicity in this seismically active area of Central Europe. The model takes into account how the crustal heterogeneity of the region effects thermal properties and is validated with a dataset of wellbore temperatures. We find that the Adriatic crust appears more mafic, through its radiogenic heat values (1.30E-06 W/m3) and maximum temperature of seismicity (600 degrees C), than the European crust (1.3-2.6E-06 W/m3 and 450 degrees C). We also show that at depths of < 10 km the thermal field is largely controlled by sedimentary blanketing or topographic effects, whilst the deeper temperature field is primarily controlled by the LAB topology and the distribution and parameterization of radiogenic heat sources within the upper crust.}, language = {en} } @article{BallingMaystrenkoScheckWenderoth2013, author = {Balling, Philipp and Maystrenko, Yuriy and Scheck-Wenderoth, Magdalena}, title = {The deep thermal field of the Gl{\"u}ckstadt Graben}, series = {Environmental earth sciences}, volume = {70}, journal = {Environmental earth sciences}, number = {8}, publisher = {Springer}, address = {New York}, issn = {1866-6280}, doi = {10.1007/s12665-013-2750-z}, pages = {3505 -- 3522}, year = {2013}, abstract = {With this paper, we assess the present-day conductive thermal field of the Glueckstadt Graben in NW Germany that is characterized by large salt walls and diapirs structuring the graben fill. We use a finite element method to calculate the 3D steady-state conductive thermal field based on a lithosphere-scale 3D structural model that resolves the first-order structural characteristics of the graben and its underlying lithosphere. Model predictions are validated against measured temperatures in six deep wells. Our investigations show that the interaction of thickness distributions and thermal rock properties of the different geological layers is of major importance for the distribution of temperatures in the deep subsurface of the Glueckstadt Graben. However, the local temperatures may result from the superposed effects of different controlling factors. Especially, the upper sedimentary part of the model exhibits huge lateral temperature variations, which correlate spatially with the shape of the thermally highly conductive Permian salt layer. Variations in thickness and geometry of the salt cause two major effects, which provoke considerable lateral temperature variations for a given depth. (1) The "chimney effect" causes more efficient heat transport within salt diapirs. As a consequence positive thermal anomalies develop in the upper part and above salt structures, where the latter are covered by much less conductive sediments. In contrast, negative thermal anomalies are noticeable underneath salt structures. (2) The "thermal blanketing effect" is caused by thermally low conductive sediments that provoke the local storage of heat where these insulating sediments are present. The latter effect leads to both local and regional thermal anomalies. Locally, this translates to higher temperatures where salt margin synclines are filled with thick insulating clastic sediments. For the regional anomalies the cumulative insulating effects of the entire sediment fill results in a long-wavelength variation of temperatures in response to heat refraction effects caused by the contrast between insulating sediments and highly conductive crystalline crust. Finally, the longest wavelength of temperature variations is caused by the depth position of the isothermal lithosphere-asthenosphere boundary defining the regional variations of the overall geothermal gradient. We find that a conductive thermal model predicts observed temperatures reasonably well for five of the six available wells, whereas the steady-state conductive approach appears not to be valid for the sixth well.}, language = {en} } @misc{SippelMeessenCacaceetal.2017, author = {Sippel, Judith and Meeßen, Christian and Cacace, Mauro and Mechie, James and Fishwick, Stewart and Heine, Christian and Scheck-Wenderoth, Magdalena and Strecker, Manfred}, title = {The Kenya rift revisited}, series = {Postprints der Universit{\"a}t Potsdam : Mathematisch Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsdam : Mathematisch Naturwissenschaftliche Reihe}, number = {644}, issn = {1866-8372}, doi = {10.25932/publishup-41822}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-418221}, pages = {45 -- 81}, year = {2017}, abstract = {We present three-dimensional (3-D) models that describe the present-day thermal and rheological state of the lithosphere of the greater Kenya rift region aiming at a better understanding of the rift evolution, with a particular focus on plume-lithosphere interactions. The key methodology applied is the 3-D integration of diverse geological and geophysical observations using gravity modelling. Accordingly, the resulting lithospheric-scale 3-D density model is consistent with (i) reviewed descriptions of lithological variations in the sedimentary and volcanic cover, (ii) known trends in crust and mantle seismic velocities as revealed by seismic and seismological data and (iii) the observed gravity field. This data-based model is the first to image a 3-D density configuration of the crystalline crust for the entire region of Kenya and northern Tanzania. An upper and a basal crustal layer are differentiated, each composed of several domains of different average densities. We interpret these domains to trace back to the Precambrian terrane amalgamation associated with the East African Orogeny and to magmatic processes during Mesozoic and Cenozoic rifting phases. In combination with seismic velocities, the densities of these crustal domains indicate compositional differences. The derived lithological trends have been used to parameterise steady-state thermal and rheological models. These models indicate that crustal and mantle temperatures decrease from the Kenya rift in the west to eastern Kenya, while the integrated strength of the lithosphere increases. Thereby, the detailed strength configuration appears strongly controlled by the complex inherited crustal structure, which may have been decisive for the onset, localisation and propagation of rifting.}, language = {en} } @article{SippelMeessenCacaceetal.2017, author = {Sippel, Judith and Meessen, Christian and Cacace, Mauro and Mechie, James and Fishwick, Stewart and Heine, Christian and Scheck-Wenderoth, Magdalena and Strecker, Manfred}, title = {The Kenya rift revisited}, series = {Solid earth}, volume = {8}, journal = {Solid earth}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {1869-9510}, doi = {10.5194/se-8-45-2017}, pages = {45 -- 81}, year = {2017}, abstract = {We present three-dimensional (3-D) models that describe the present-day thermal and rheological state of the lithosphere of the greater Kenya rift region aiming at a better understanding of the rift evolution, with a particular focus on plume-lithosphere interactions. The key methodology applied is the 3-D integration of diverse geological and geophysical observations using gravity modelling. Accordingly, the resulting lithospheric-scale 3-D density model is consistent with (i) reviewed descriptions of lithological variations in the sedimentary and volcanic cover, (ii) known trends in crust and mantle seismic velocities as revealed by seismic and seismological data and (iii) the observed gravity field. This data-based model is the first to image a 3-D density configuration of the crystalline crust for the entire region of Kenya and northern Tanzania. An upper and a basal crustal layer are differentiated, each composed of several domains of different average densities. We interpret these domains to trace back to the Precambrian terrane amalgamation associated with the East African Orogeny and to magmatic processes during Mesozoic and Cenozoic rifting phases. In combination with seismic velocities, the densities of these crustal domains indicate compositional differences. The derived lithological trends have been used to parameterise steady-state thermal and rheological models. These models indicate that crustal and mantle temperatures decrease from the Kenya rift in the west to eastern Kenya, while the integrated strength of the lithosphere increases. Thereby, the detailed strength configuration appears strongly controlled by the complex inherited crustal structure, which may have been decisive for the onset, localisation and propagation of rifting.}, language = {en} } @misc{GholamrezaieScheckWenderothSippeletal.2018, author = {Gholamrezaie, Ershad and Scheck-Wenderoth, Magdalena and Sippel, Judith and Strecker, Manfred}, title = {Variability of the geothermal gradient across two differently aged magma-rich continental rifted margins of the Atlantic Ocean}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-409493}, pages = {19}, year = {2018}, abstract = {Abstract. The aim of this study is to investigate the shallow thermal field differences for two differently aged passive continental margins by analyzing regional variations in geothermal gradient and exploring the controlling factors for these variations. Hence, we analyzed two previously published 3-D conductive and lithospheric-scale thermal models of the Southwest African and the Norwegian passive margins. These 3-D models differentiate various sedimentary, crustal, and mantle units and integrate different geophysical data such as seismic observations and the gravity field. We extracted the temperature-depth distributions in 1 km intervals down to 6 km below the upper thermal boundary condition. The geothermal gradient was then calculated for these intervals between the upper thermal boundary condition and the respective depth levels (1, 2, 3, 4, 5, and 6 km below the upper thermal boundary condition). According to our results, the geothermal gradient decreases with increasing depth and shows varying lateral trends and values for these two different margins. We compare the 3-D geological structural models and the geothermal gradient variations for both thermal models and show how radiogenic heat production, sediment insulating effect, and thermal lithosphere-asthenosphere boundary (LAB) depth influence the shallow thermal field pattern. The results indicate an ongoing process of oceanic mantle cooling at the young Norwegian margin compared with the old SW African passive margin that seems to be thermally equilibrated in the present day.}, language = {en} } @article{GholamrezaieScheckWenderothSippeletal.2018, author = {Gholamrezaie, Ershad and Scheck-Wenderoth, Magdalena and Sippel, Judith and Strecker, Manfred}, title = {Variability of the geothermal gradient across two differently aged magma-rich continental rifted margins of the Atlantic Ocean}, series = {Solid Earth}, volume = {9}, journal = {Solid Earth}, number = {1}, publisher = {Copernicus}, address = {G{\"o}ttingen}, issn = {1869-9529}, doi = {10.5194/se-9-139-2018}, pages = {139 -- 158}, year = {2018}, abstract = {Abstract. The aim of this study is to investigate the shallow thermal field differences for two differently aged passive continental margins by analyzing regional variations in geothermal gradient and exploring the controlling factors for these variations. Hence, we analyzed two previously published 3-D conductive and lithospheric-scale thermal models of the Southwest African and the Norwegian passive margins. These 3-D models differentiate various sedimentary, crustal, and mantle units and integrate different geophysical data such as seismic observations and the gravity field. We extracted the temperature-depth distributions in 1 km intervals down to 6 km below the upper thermal boundary condition. The geothermal gradient was then calculated for these intervals between the upper thermal boundary condition and the respective depth levels (1, 2, 3, 4, 5, and 6 km below the upper thermal boundary condition). According to our results, the geothermal gradient decreases with increasing depth and shows varying lateral trends and values for these two different margins. We compare the 3-D geological structural models and the geothermal gradient variations for both thermal models and show how radiogenic heat production, sediment insulating effect, and thermal lithosphere-asthenosphere boundary (LAB) depth influence the shallow thermal field pattern. The results indicate an ongoing process of oceanic mantle cooling at the young Norwegian margin compared with the old SW African passive margin that seems to be thermally equilibrated in the present day.}, language = {en} } @misc{GholamrezaieScheckWenderothSippeletal.2018, author = {Gholamrezaie, Ershad and Scheck-Wenderoth, Magdalena and Sippel, Judith and Strecker, Manfred}, title = {Variability of the geothermal gradient across two differently aged magma-rich continental rifted margins of the Atlantic Ocean}, series = {Postprints der Universit{\"a}t Potsadm : Mathematisch-Naturwissenschaftliche Reihe}, journal = {Postprints der Universit{\"a}t Potsadm : Mathematisch-Naturwissenschaftliche Reihe}, number = {621}, issn = {1866-8372}, doi = {10.25932/publishup-41821}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-418210}, pages = {20}, year = {2018}, abstract = {The aim of this study is to investigate the shal- low thermal field differences for two differently aged pas- sive continental margins by analyzing regional variations in geothermal gradient and exploring the controlling factors for these variations. Hence, we analyzed two previously pub- lished 3-D conductive and lithospheric-scale thermal models of the Southwest African and the Norwegian passive mar- gins. These 3-D models differentiate various sedimentary, crustal, and mantle units and integrate different geophysi- cal data such as seismic observations and the gravity field. We extracted the temperature-depth distributions in 1 km intervals down to 6 km below the upper thermal boundary condition. The geothermal gradient was then calculated for these intervals between the upper thermal boundary condi- tion and the respective depth levels (1, 2, 3, 4, 5, and 6 km below the upper thermal boundary condition). According to our results, the geothermal gradient decreases with increas- ing depth and shows varying lateral trends and values for these two different margins. We compare the 3-D geologi- cal structural models and the geothermal gradient variations for both thermal models and show how radiogenic heat pro- duction, sediment insulating effect, and thermal lithosphere- asthenosphere boundary (LAB) depth influence the shallow thermal field pattern. The results indicate an ongoing process of oceanic mantle cooling at the young Norwegian margin compared with the old SW African passive margin that seems to be thermally equilibrated in the present day.}, language = {en} }