In this work, we focus on improving the understanding of subsurface processes with respect to interactions with climate dynamics. We present advanced, open-source mathematical methods that enable us to investigate the influence of various model properties on the final outcomes. By relying on our approach, we have been able to showcase their importance in improving our understanding of the subsurface and highlighting the current shortcomings of currently adopted models.

The Earth’s crust beneath the Caribbean Sea formed at about 90 Ma due to large magmatic activity of a mantle plume, which brought molten material up from the deep Earth. By integrating diverse geophysical datasets, we image for the first time two fossil magmatic conduits beneath the Caribbean. The location of these conduits at 90 Ma does not correspond with the present-day Galápagos plume. Either this mantle plume migrated in time or these conduits were formed above another unknown plume.

By comparing long term lithospheric strength to seismicity patterns across the Alpine region, we show that most seismicity occurs where strengths are highest within the crust. The lower crust appears largely aseismic due to energy being dissipated by ongoing creep from low viscosities. Lithospheric structure appears to exert a primary control on seismicity distribution, with both forelands display a different distribution patterns, likely reflecting their different tectonic settings.

We have developed a first Germany-wide 3D data-based density and temperature model integrating geoscientific observations and physical processes. The model can serve as a reference for local detailed studies dealing with temperature, pressure, stress, subsidence and sedimentation. Our results help to improve subsurface utilization concepts, reveal current geomechanical conditions crucial for hazard assessment and gather information on viable resources such groundwater and deep geothermal energy.

By utilising both the observed gravity field of the Alps and their forelands and indications from deep seismic surveys, we were able to produce a 3-D structural model of the region that indicates the distribution of densities within the lithosphere. We found that the present-day Adriatic crust is both thinner and denser than the European crust and that the properties of Alpine crust are strongly linked to their provenance.

In this study we investigate groundwater flow in the deep subsurface of the Upper Rhine Graben. We make use of a 3-D numerical model covering the entire Upper Rhine Graben. The deep hydrodynamics are characterized by fluid flow from the graben flanks towards its center and in the southern half of the graben from south to north. Moreover, local heterogeneities in the shallow flow field arise from the interaction between regional groundwater flow and the heterogeneous sedimentary configuration.

The focus of the paper is to evaluate the permeability change of a matrix-fracture systems under mechanical loading und to understand the processes behind. This evaluation is based on data from laboratory experiments in comparison to 3-D numerical modelling results.

The study presented in this paper aims at reproducing findings from chemical and isotopic groundwater sample analysis along with quantifying the influence of regional (cross-boundary) flow for the area of Berlin, Germany. For this purpose we built 3-D models of the subsurface, populating them with material parameters (e.g. porosity, permeability) and solving them for coupled fluid and heat transport. Special focus was given to the setup of boundary conditions, i.e. fixed pressure at the sides.

Based on geophysical data integration and 3-D gravity modeling, we show that significant density heterogeneities are expressed as two large high-density bodies in the crust below the Sea of Marmara. The location of these bodies correlates spatially with the bends of the main Marmara fault, indicating that rheological contrasts in the crust may influence the fault kinematics. Our findings may have implications for seismic hazard and risk assessments in the Marmara region.

Accuracy of well inflow modelling in different numerical simulation approaches was compared for a multi-lateral well with laterals of varying diameter. For homogeneous cases, all simulators generally were reasonably close in terms of the total well flow (deviations smaller than 4 %). The distribution of the flow over the different laterals in a well can vary significantly between simulators (> 20 %). In a heterogeneous case with a fault the deviations between the approaches were much larger.

We examined the thermal gradient as an index of the thermal field in the Atlantic. While the thermal anomaly in the South Atlantic should be equilibrated, the thermal disturbance in the North Atlantic causes thermal effects in the present day. Characteristics of the lithosphere ultimately determine the thermal field. The thermal gradient nonlinearly decreases with depth and varies significantly both laterally and with time, which has implications for methods of thermal history reconstruction.

The paper describes theory and numerical implementation for coupled thermo–hydraulic–mechanical processes focusing on reservoir (mainly related to geothermal energy) applications.

The Kenya Rift is a zone along which the African continental plate is stretched as evidenced by strong earthquake and volcanic activity. We want to understand the controlling factors of past and future tectonic deformation; hence, we assess the structural and strength configuration of the rift system at the present-day. Data-driven 3-D numerical models show how the inherited composition of the crust and a thermal anomaly in the deep mantle interact to form localised zones of tectonic weakness.

Subsurface engineering relies on sparsely distributed data points of the stress state of the earth's crust. 3D geomechanical--numerical modelling is applied to estimate the stress state in the entire volume of a large area. We present a multi-stage approach of differently sized models which provide the stress state in an area of interest derived from few and widely scattered data records. Furthermore we demonstrate the changes in reliability of the model depending on different input parameters.

We introduce a regional 3-D structural model of the Barents Sea and Kara Sea region which is the first to combine information on five sedimentary units and the crystalline crust as well as the configuration of the lithospheric mantle. By relating the shallow and deep structures for certain tectonic subdomains, we shed new light on possible causative basin-forming mechanisms that we discuss.