Since the start of production in 1968 in the Groningen gas field (Netherlands) considerable land subsidence (>30 cm) has occurred above the field. Variability in reservoir compaction has led to earthquakes on reactivated Mesozoic age reservoir faults. Even though the impact of this seismicity (MW ≤3.6) on society has been large, due to substantial structural damage to buildings, surface deformation induced by the co-seismic slip has been too small to detect using geodetic data. It is possible that differential compaction across faults is not only accommodated by seismic slip, but also by aseismic slip (e.g., creep). Aseismically slipping reservoir faults would relax the stresses in the reservoir and, thus, reduce the severity of the seismicity. In this study we explore the potential occurrence of aseismic slip on the reservoir faults. We perform a sensitivity analysis to investigate whether aseismic slip on the different reservoir faults has the capacity to produce detectable surface signals. We use the analytical Okada (1992) model of slip on a discrete dislocation in a uniform elastic half-space to simulate the deformations originating from slip on a wide range of fault geometries, representing the variability in the field. Unsurprisingly, laterally extensive faults with strong compaction contrasts across them (large differential slip magnitudes) produce the largest surface signals. To determine which potentially aseismically slipping faults produce surface signals that could be detectable in persistent scatterer InSAR time series, we analyze the surface patterns for large differential displacements across large length scales, since InSAR observations are most sensitive to spatially extensive patterns with high spatial gradients. We use the results of the sensitivity analysis to guide our search for patterns originating from aseismically slipping reservoir faults in PS-InSAR time series data of the Groningen area. First results show that these specific patterns are rare, indicating that the amount of aseismic slip is limited. For faults lacking surface signals related aseismic slip, the results of sensitivity analysis are used to determine upper limits for the aseismic differential slip magnitudes.

Plate reconstruction studies show that the Neotethys Ocean was closing due to convergence of Africa and Eurasia towards the end of the Cretaceous. The period around 75 Ma reflects the onset of continental collision between the two plates, although convergence was still mainly accommodated by subduction, with the Neotethys slab subducting beneath Eurasia. Africa was separated from the rapidly north moving Indian plate by the Owen oceanic transform in the northeast. The rest of the plate was surrounded by mid-ocean ridges. Geologic observations in large basins show that Africa was experiencing continent-wide rifting related to northeast-southwest extension. We aim to quantify the forces and related paleostresses associated with this tectonic setting. To constrain these forces, we use the latest plate kinematic reconstructions, while balancing horizontal gravitational stresses, plate boundary forces and the plate's interaction with the underlying mantle. The contribution of dynamic topography to horizontal gravitational stresses is based on recent mantle convection studies. We model intraplate stresses and compare them with the strain observations. We find that slab pull, horizontal gravitational stresses and transform shear tractions in general acted with the same orientation as the absolute motion of the African plate 75 Ma. Both the balance between these three and the other, resistive, forces, and the fit to strain observations require the net slab pull, as experienced by the plate, to be low, pointing to the absence of a mature continuous Neotethys subduction zone at the time. This corresponds well to reconstructions of micro-continents interfering with the Neotethyan subduction.