Refine
Year of publication
Language
- English (21)
Is part of the Bibliography
- yes (21)
Keywords
- Northeast German Basin (3)
- Argentine margine (2)
- European basin system (2)
- Geomechanical Model (2)
- Long-Term (2)
- Middle Strand (2)
- NW Turkey (2)
- North Anatolian Fault (2)
- Northwestern Anatolia (2)
- Pull-Apart (2)
Institute
- Institut für Geowissenschaften (21) (remove)
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.
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.
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.
Dissolved CO2 storage in geological formations with low pressure, low risk and large capacities
(2017)
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.
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.
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.
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.
The Kenya rift revisited
(2017)
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.
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.