A large non-double-couple component of a tectonic earthquake indicates that its rupture likely was complex and likely involved multiple faults. Detailed source models of such earthquakes can add to our understanding of earthquake source complexity. The 2007 Martinique earthquake in the Caribbean Sea is one of the largest recent earthquakes with a known large non-double-couple component. It was an intermediate depth intraslab earthquake within the South America plate where it is subducting beneath the Caribbean plate. We applied potency density tensor inversion (PDTI) to teleseismic P waves generated by the 2007 Martinique earthquake to model its source processes and focal mechanism distribution. We identified two focal mechanisms: a strike-slip mechanism with a north–south T axis, and a down-dip extension (DDE) mechanism with an east–west T axis. Rupture by the DDE mechanism was predominant in the northern part of the source region and strike-slip rupture in the southern part. These two focal mechanisms had approximately parallel P axes and approximately orthogonal T axes. The seismic moments released by both types of rupture were almost equal. These results indicate that the 2007 Martinique earthquake had a large non-double-couple component. We identified five sub-events with two predominant directions of rupture propagation: two strike-slip sub-events propagated to the southeast and three DDE sub-events propagated to the east. Although the directions of propagation were consistent for each focal mechanism, each sub-event appears to have occurred in isolation. For example, the rupture of one DDE sub-event propagated from the edge of the source region back towards the hypocentre. Complex ruptures that include multiple sub-events may be influenced by high pore fluid pressure associated with slab dehydration. Our results show that PDTI can produce stable estimates of complex seismic source processes and provide useful information about the sources of complex intermediate depth intraslab earthquakes for which fault geometry assumptions are difficult.

The 2016 Kaikoura earthquake, New Zealand, ruptured more than a dozen faults, making it difficult to prescribe a model fault for analysing the event by inversion. To model this earthquake from teleseismic records, we used a potency density tensor inversion, which projects multiple fault slips onto a single model fault plane, which reduced the non-uniqueness due to the uncertainty in selecting the faults’ orientations. The resulting distribution of potency-rate density tensors is consistent with observed surface ruptures. In its initial stage, the rupture propagated northeastward primarily at shallow depths, and the rupture propagated northeastward at deep depths beneath a gap in reported surface ruptures. The main rupture phase started in the northeastern part of the Kekerengu fault after 50 s and propagated bilaterally to the northeast and southwest. The non-double-couple component grew to a large fraction of the source elements as the rupture went through the junction of the Jordan Thrust and the Papatea fault, which suggests that the rupture branched into both faults as it back-propagated toward the southwest. The potency density tensor inversion sheds new light on the irregular evolution of this earthquake, which produced a fault rupture pattern of unprecedented complexity. Our source model should provide new insights into source process of the 2016 Kaikoura earthquake (e.g., back-rupture propagation), which should prompt research to determine a more realistic model with segmented faults using near-field data.

The 2016 Kaikoura earthquake, New Zealand, ruptured more than a dozen faults, making it difficult to prescribe a model fault for analyzing the event by inversion. To model this earthquake from teleseismic records, we used a potency density tensor inversion, which projects multiple fault slips onto a single model fault plane, reducing the non-uniqueness due to the uncertainty in selecting the faults’ orientations. The resulting distribution of potency-rate density tensors is consistent with observed surface ruptures. In its initial stage, the rupture propagated northeastward primarily at shallow depths. Later, the rupture propagated northeastward at greater depths beneath a gap in reported surface ruptures. The main rupture phase started in the northeastern part of the Kekerengu fault after 50 s and propagated bilaterally to the northeast and southwest. The non-double-couple component grew to a large fraction of the source elements as the rupture went through the junction of the Jordan Thrust and the Papatea fault, which suggests that the rupture branched into both faults as it back-propagated toward the southwest. The potency density tensor inversion sheds new light on the irregular evolution of this earthquake, which produced a fault rupture pattern of unprecedented complexity. Our source model of the 2016 Kaikoura earthquake (e.g., back-rupture propagation) could prompt research to determine a more realistic model with segmented faults using near-field data.