TY - JOUR A1 - Kaiser, Björn Onno A1 - Cacace, Mauro A1 - Scheck-Wenderoth, Magdalena T1 - 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 JF - Environmental earth sciences N2 - 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. KW - Mesh convergence KW - Conduction KW - Advection KW - Convection KW - Thermal field KW - Northeast German Basin Y1 - 2013 U6 - https://doi.org/10.1007/s12665-013-2249-7 SN - 1866-6280 SN - 1866-6299 VL - 70 IS - 8 SP - 3643 EP - 3659 PB - Springer CY - New York ER - TY - JOUR A1 - Kaiser, Bjoern Onno A1 - Cacace, Mauro A1 - Scheck-Wenderoth, Magdalena A1 - Lewerenz, Bjoern T1 - Characterization of main heat transport processes in the Northeast German Basin constraints from 3-D numerical models JF - Geochemistry, geophysics, geosystems N2 - 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. KW - advection KW - convection KW - coupled fluid and heat transport KW - numerical simulations KW - Northeast German Basin KW - salt structures Y1 - 2011 U6 - https://doi.org/10.1029/2011GC003535 SN - 1525-2027 VL - 12 IS - 13 PB - American Geophysical Union CY - Washington ER - TY - JOUR A1 - Rodriguez Piceda, Constanza A1 - Scheck-Wenderoth, Magdalena A1 - Bott, Judith A1 - Gomez Dacal, Maria Laura A1 - Cacace, Mauro A1 - Pons, Michael A1 - Prezzi, Claudia A1 - Strecker, Manfred T1 - Controls of the Lithospheric Thermal Field of an Ocean-Continent Subduction Zone BT - the Southern Central Andes JF - Lithosphere / Geological Society of America N2 - 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. Y1 - 2022 U6 - https://doi.org/10.2113/2022/2237272 SN - 1941-8264 SN - 1947-4253 VL - 2022 IS - 1 PB - GeoScienceWorld CY - McLean ER - TY - JOUR A1 - Cherubini, Yvonne A1 - Cacace, Mauro A1 - Scheck-Wenderoth, Magdalena A1 - Moeck, Inga A1 - Lewerenz, Björn T1 - Controls on the deep thermal field - implications from 3-D numerical simulations for the geothermal research site Groß Schönebeck JF - Environmental earth sciences N2 - 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. KW - Thermal field KW - Coupled fluid and heat transport KW - Faults KW - Groß beta Schönebeck Y1 - 2013 U6 - https://doi.org/10.1007/s12665-013-2519-4 SN - 1866-6280 SN - 1866-6299 VL - 70 IS - 8 SP - 3619 EP - 3642 PB - Springer CY - New York ER - TY - JOUR A1 - Cacace, Mauro A1 - Kaiser, Bjoern Onno A1 - Lewerenz, Bjoern A1 - Scheck-Wenderoth, Magdalena T1 - Geothermal energy in sedimentary basins : what we can learn from regional numerical models N2 - 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. Y1 - 2010 UR - http://www.sciencedirect.com/science/journal/00092819 U6 - https://doi.org/10.1016/j.chemer.2010.05.017 SN - 0009-2819 ER - TY - JOUR A1 - Spooner, Cameron A1 - Scheck-Wenderoth, Magdalena A1 - Cacace, Mauro A1 - Anikiev, Denis T1 - How Alpine seismicity relates to lithospheric strength JF - International journal of earth sciences N2 - 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. KW - lithosphere KW - strength KW - rheology KW - 3D-Model KW - Alps KW - seismicity Y1 - 2022 U6 - https://doi.org/10.1007/s00531-022-02174-5 SN - 1437-3254 SN - 1437-3262 VL - 111 IS - 4 SP - 1201 EP - 1221 PB - Springer CY - Berlin ; Heidelberg ER - TY - JOUR A1 - Degen, Denise A1 - Spooner, Cameron A1 - Scheck-Wenderoth, Magdalena A1 - Cacace, Mauro T1 - How biased are our models? BT - a case study of the alpine region JF - Geoscientific model development : an interactive open access journal of the European Geosciences Union N2 - 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. Y1 - 2021 U6 - https://doi.org/10.5194/gmd-14-7133-2021 SN - 1991-959X SN - 1991-9603 VL - 14 IS - 11 SP - 7133 EP - 7153 PB - Copernicus CY - Göttingen ER - TY - JOUR A1 - Cherubini, Yvonne A1 - Cacace, Mauro A1 - Blöcher, Guido A1 - Scheck-Wenderoth, Magdalena T1 - Impact of single inclined faults on the fluid flow and heat transport - results from 3-D finite element simulations JF - Environmental earth sciences N2 - 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. KW - Hydrothermal field KW - 3-D numerical simulations KW - Inclined faults KW - Fault zone KW - Coupled fluid flow and heat transport KW - Finite elements Y1 - 2013 U6 - https://doi.org/10.1007/s12665-012-2212-z SN - 1866-6280 SN - 1866-6299 VL - 70 IS - 8 SP - 3603 EP - 3618 PB - Springer CY - New York ER - TY - JOUR A1 - Freymark, Jessica A1 - Bott, Judith A1 - Cacace, Mauro A1 - Ziegler, Moritz 0. A1 - Scheck-Wenderoth, Magdalena T1 - Influence of the Main Border Faults on the 3D Hydraulic Field of the Central Upper Rhine Graben JF - Geofluids N2 - 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). Y1 - 2019 U6 - https://doi.org/10.1155/2019/7520714 SN - 1468-8115 SN - 1468-8123 PB - Wiley-Hindawi CY - London ER - TY - JOUR A1 - Rodriguez Piceda, Constanza A1 - Scheck-Wenderoth, Magdalena A1 - Cacace, Mauro A1 - Bott, Judith A1 - Strecker, Manfred T1 - Long-Term Lithospheric Strength and Upper-Plate Seismicity in the Southern Central Andes, 29 degrees-39 degrees S JF - Geochemistry, geophysics, geosystems N2 - 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. KW - subduction zone KW - Andes KW - rheology KW - seismicity KW - flat-slab Y1 - 2022 U6 - https://doi.org/10.1029/2021GC010171 SN - 1525-2027 VL - 23 IS - 3 PB - American Geophysical Union CY - Washington ER -