Tectonics and seismicity in the Northern Apennines driven by slab retreat and lithospheric delamination

Understanding how long-term subduction dynamics relates to the short-term seismicity and crustal tec tonics is a challenging but crucial topic in seismotectonics. We attempt to address this issue by linking long-term geodynamic evolution with short-term seismogenic deformation in the Northern Apennines. This retreating subduction orogen displays tectonic and seismogenic behaviors on various spatiotemporal scales that also characterize other subduction zones in the Mediterranean area. We use visco-elasto-plastic seismo-thermo-mechanical (STM) modeling with a realistic 2D setup based on available geological and geophysical data. The subduction dynamics and seismicity are coupled in the numerical modeling, and driven only by buoyancy forces, i.e., slab pull. Our results suggest that lower crustal rheology and lithospheric mantle temperature modulate the crustal tectonics of the Northern Apennines, as inferred by previous studies. The observed spatial distribution of upper crustal tectonic regimes and surface displacements requires buoyant, highly ductile material in the subduction channel beneath the internal part of the orogen. This allows protrusion of the asthenosphere in the lower crust and lithospheric delamination associated with slab retreat. The resulting surface velocities and principal stress axes generally agree with present-day observations, suggesting that slab delamination and retreat can explain the dynamics of the orogen. Our simulations successfully reproduce the type and overall distribution of seismicity with thrust faulting events in the external part of the orogen and normal faulting in its internal part. Slab temperatures and lithospheric mantle stiffness affect the cumulative seismic moment release and spatial distribution of upper crustal earthquakes. The properties of deep, sub-crustal material are thus shown to influence upper crustal seismicity in an orogen driven by slab retreat, even though the upper crust is largely decoupled from the lithospheric mantle. Our simulations therefore highlight the effect of deep lower crustal rheologies, self-driven subduction dynamics and mantle properties in controlling shallow deformation and seismicity