Finite-fault slip models for the 15 April 1979 (Mw 7.1) Montenegro earthquake and its strongest aftershock of 24 May 1979 (Mw 6.2)

We revisit the April 1979 Montenegro earthquake sequence to invert for finite-fault slip models for the mainshock of 15 April 1979 (Mw 7.1) and of the strongest aftershock of 24 May 1979 (Mw 6.2) using P, SH and SV waveforms, retrieved from IRIS data center. We also used body waveform modelling inversion to confirm the focal mechanism of the mainshock as a pure thrust mechanism and rule out the existence of considerable strike slip component in the motion. The mainshock occurred along a shallow (depth 7 km), low angle (14°) thrust fault, parallel to the coastline and dipping to the NE. Our preferred slip distribution model for the mainshock indicates that rupture initiated from SE and propagated towards NW, with a speed of 2.0 km/s. Moment was released in a main slip patch, confined in an area of L ∼ 50 km × W ∼ 23 km. The maximum slip (∼ 2.7 m) occurred ∼ 30 km to the NW of the hypocenter (location of rupture initiation). The average slip is 49 cm and the total moment release over the fault is 4.38e19 Nm. The slip model adequately fits the distribution of the Mw ≥ 4.3 aftershocks, as most of them are located in the regions of the fault plane that did not slip during the mainshock. The 24 May 1979 (Mw 6.2) strongest aftershock occurred ∼ 40 km NW of the mainshock. Our preferred slip model for this event showed a characteristic two-lobe pattern, where each lobe is ∼ 7.5 × 7.5 km2. Rupture initiated in the NW lobe, where the slip obtained its maximum value of 45 cm, very close to the hypocenter, and propagated towards the south-eastern lobe where it reached another maximum value — for this lobe — of 30 cm, approximately 10 km away from the hypocenter. To indirectly validate our slip models we produced synthetic PGV maps (Shake maps) and we compared our predictions with observations of ground shaking from strong motion records. All comparisons were made for rock soil conditions and in general our slip models adequately fit the observations especially at the closest stations where the shaking was considerably stronger. Through the search of the parameter space for our inversions we obtained an optimum location for the mainshock at 42.04°N and 19.21° E and we also observed that better fit to the observations was obtained when the fault was modeled as a blind thrust fault.