Slip on a mapped normal fault for the 28th December 1908 Messina earthquake (Mw 7.1) in Italy.

The 28th December 1908 Messina earthquake (Mw 7.1), Italy, caused >80,000 deaths and transformed earthquake science by triggering the study of earthquake environmental effects worldwide, yet its source is still a matter of debate. To constrain the geometry and kinematics of the earthquake we use elastic half-space modelling on non-planar faults, constrained by the geology and geomorphology of the Messina Strait, to replicate levelling data from 1907-1909. The novelty of our approach is that we (a) recognise the similarity between the pattern of vertical motions and that of other normal faulting earthquakes, and (b) for the first time model the levelling data using the location and geometry of a well-known offshore capable fault. Our results indicate slip on the capable fault with a dip to the east of 70° and 5 m dip-slip at depth, with slip propagating to the surface on the sea bed. Our work emphasises that geological and geomorphological observations supporting maps of capable non-planar faults should not be ignored when attempting to identify the sources of major earthquakes.

Orogen-scale uplift in the central Italian Apennines drives episodic behaviour of earthquake faults.

Many areas of the Earth's crust deform by distributed extensional faulting and complex fault interactions are often observed. Geodetic data generally indicate a simpler picture of continuum deformation over decades but relating this behaviour to earthquake occurrence over centuries, given numerous potentially active faults, remains a global problem in hazard assessment. We address this challenge for an array of seismogenic faults in the central Italian Apennines, where crustal extension and devastating earthquakes occur in response to regional surface uplift. We constrain fault slip-rates since ~18 ka using variations in cosmogenic 36Cl measured on bedrock scarps, mapped using LiDAR and ground penetrating radar, and compare these rates to those inferred from geodesy. The 36Cl data reveal that individual faults typically accumulate meters of displacement relatively rapidly over several thousand years, separated by similar length time intervals when slip-rates are much lower, and activity shifts between faults across strike. Our rates agree with continuum deformation rates when averaged over long spatial or temporal scales (104 yr; 102 km) but over shorter timescales most of the deformation may be accommodated by <30% of the across-strike fault array. We attribute the shifts in activity to temporal variations in the mechanical work of faulting.