Mainshock Mw Research Articles

Kinematic reconstruction of active tectonic and halokinetic structures in the 2021 NW Palagruža earthquake area (Central Adriatic)

In March 2021, a compressional earthquake sequence (mainshock Mw 5.2) occurred in the Central Adriatic Sea, offshore Croatia. The struck area is characterized by a complex tectonic settings due to due interaction between tectonic and halokinetic structures. Former studies in this region, mostly based on geophysical and seismological data, do not provide a comprehensive description of the geological complexities of the area, caused by the presence of different types of active structures. Following the interpretation and depth conversion of a set of seismic reflection profiles, we present a kinematic restoration to obtain the shortening rates for the last 6 Myr. We highlight the presence of three main types of active structures: i) shallow thrusts and related folds deforming the seafloor, ii) deep thrusts promoting large-scale deformation, iii) halokinetic structures deforming, at least, the Messinian. We highlight how the structural setting of the area is more complex than previously interpreted, with a possible decoupling of shallow and deep structures and an interaction between tectonic and halokinetic structures. This work opens new perspectives about the role of halokinetic processes in active seismogenic areas.

The 2013–2018 Matese and Beneventano Seismic Sequences (Central–Southern Apennines): New Constraints on the Hypocentral Depth Determination

The Matese and Beneventano areas coincide with the transition from the central to the southern Apennines and are characterized by both SW- and NE-dipping normal faulting seismogenic structures, responsible for the large historical earthquakes. We studied the Matese and Beneventano seismicity by means of high-precision locations of earthquakes spanning from 29 December 2013 to 4 September 2018. Events were located by using all of the available data from temporary and permanent stations in the area and a 1D computed velocity model, inverting the dataset with the Velest code. For events M > 2.8 we used P- and S-waves arrival times of the strong motion stations located in the study area. A constant value of 1.83 for Vp/Vs was computed with a modified Wadati method. The dataset consists of 2378 earthquakes, 18,715 P- and 12,295 S-wave arrival times. We computed 55 new fault plane solutions. The mechanisms show predominantly normal fault movements, with T-axis trends oriented NE–SW. Only relatively small E–W trending clusters in the eastern peripheral zones of the Apenninic belt show right-lateral strike-slip kinematics similar to that observed in the Potenza (1990–1991) and Molise (2002 and 2018) sequences. These belong to transfer zones associated with differential slab retreat of the Adriatic plate subduction beneath the Apennines. The Matese sequence (December 2013–February 2014; main shock Mw 5.0) is the most relevant part of our dataset. Hypocentral depths along the axis of the Apenninic belt are in agreement with previous seismological studies that place most of the earthquakes in the brittle upper crust. We confirm a general deepening of seismicity moving from west to the east along the Apennines. Seismicity depth is controlled by heat-flow, which is lower in the eastern side, thus causing a deeper brittle–ductile transition.

Aseismic transient during the 2010\u20132014 seismic swarm: evidence for longer recurrence of M\u2009\u2265\u20096.5 earthquakes in the Pollino gap (Southern Italy)?

In actively deforming regions, crustal deformation is accommodated by earthquakes and through a variety of transient aseismic phenomena. Here, we study the 2010–2014 Pollino (Southern Italy) swarm sequence (main shock MW 5.1) located within the Pollino seismic gap, by analysing the surface deformation derived from Global Positioning System and Synthetic Aperture Radar data. Inversions of geodetic time series show that a transient slip, with the same mechanism of the main shock, started about 3–4 months before the main shock and lasted almost one year, evolving through time with acceleration phases that correlate with the rate of seismicity. The moment released by the transient slip is equivalent to MW 5.5, significantly larger than the seismic moment release revealing therefore that a significant fraction of the overall deformation is released aseismically. Our findings suggest that crustal deformation in the Pollino gap is accommodated by infrequent “large” earthquakes (MW ≥ 6.5) and by aseismic episodes releasing a significant fraction of the accrued strain. Lower strain rates, relative to the adjacent Southern Apennines, and a mixed seismic/aseismic strain release are in favour of a longer recurrence for large magnitude earthquakes in the Pollino gap.

Structure of a normal seismogenic fault zone in carbonates: The Vado di Corno Fault, Campo Imperatore, Central Apennines (Italy)

The Vado di Corno Fault Zone (VCFZ) is an active extensional fault cutting through carbonates in the Italian Central Apennines. The fault zone was exhumed from ∼2 km depth and accommodated a normal throw of ∼2 km since Early-Pleistocene. In the studied area, the master fault of the VCFZ dips N210/54° and juxtaposes Quaternary colluvial deposits in the hangingwall with cataclastic dolostones in the footwall. Detailed mapping of the fault zone rocks within the ∼300 m thick footwall-block evidenced the presence of five main structural units (Low Strain Damage Zone, High Strain Damage Zone, Breccia Unit, Cataclastic Unit 1 and Cataclastic Unit 2). The Breccia Unit results from the Pleistocene extensional reactivation of a pre-existing Pliocene thrust. The Cataclastic Unit 1 forms a ∼40 m thick band lining the master fault and recording in-situ shattering due to the propagation of multiple seismic ruptures. Seismic faulting is suggested also by the occurrence of mirror-like slip surfaces, highly localized sheared calcite-bearing veins and fluidized cataclasites. The VCFZ architecture compares well with seismological studies of the L'Aquila 2009 seismic sequence (mainshock MW 6.1), which imaged the reactivation of shallow-seated low-angle normal faults (Breccia Unit) cut by major high-angle normal faults (Cataclastic Units).

Wairau Basin and fault connections across Cook Strait, New Zealand: seismic and geological evidence

Cook Strait separates North Island and South Island of New Zealand. It contains the Wairau and other sedimentary basins, which have sedimentary thicknesses of up to 4 km. The Strait overlies part of the Australia–Pacific plate boundary where subduction below the North Island changes to strike-slip motion through the South Island. Strike-slip faults are well documented in the South Island and less so in the North Island. Cook Strait has been surveyed by a variety of shallow penetration seismic in the past, and this has resulted in various models and uncertainty about how major faults join across the strait. The shallow-penetration data image numerous fault traces close to the sea bed, most of which can be traced for only a short distance. Gaps between individual fault traces are common. In recent years, a growing body of research, government and oil company seismic data has imaged structures much more deeply and is the subject of this paper. The Wairau Basin structural style is one dominated by normal faulting, with expansion of basin area and progressive eastwards shift in depocentres with time. A steep structural gradient or faulted margin to basement subcrop and outcrop on the North Island margin indicates a significant north-northwest to northwest component of fault control in this area of Neogene development. The seismic interpretation of this deeper data suggests there are continuous seismogenic zones underlying the discontinuous young seafloor fault scarps. Such discontinuities have led to a perception that Cook Strait faults are not connecting through, that fault rupturing across Cook Strait is improbable. Fault connections across Cook Strait along seismogenic zones have recently been interpreted from historical records. These indicate that the October 1848 earthquakes involved mainshock (Mw 7.4–7.7) displacement on the Awatere Fault in the South Island and probable aftershock displacement on the Ohariu and Wellington faults in the North Island. The January–February 1855 earthquakes (mainshock Mw ∼8.2–8.4) occurred on the Wairarapa Fault in southern North Island, Wharekauhau, Nicholson Bank faults in Cook Strait, Vernon and Awatere faults in the South Island, and Needles Fault off the NE Marlborough coast. These interpretations are supported by recent seismic data that indicate underlying basin-defining seismogenic zones connect across Cook Strait, with movements that date back at least to the top Paleogene.

Reassessment of the source of the 1976 Friuli, NE Italy, earthquake sequence from the joint inversion of high-precision levelling and triangulation data

SUMMARY In this study, we revisit the mechanism of the 1976 Friuli (NE Italy) earthquake sequence (main shocks Mw 6.4, 5.9 and 6.0). We present a new source model that simultaneously fits all the available geodetic measurements of the observed deformation. We integrate triangulation measurements, which have never been previously used in the source modelling of this sequence, with high-precision levelling that covers the epicentral area. We adopt a mixed linear/non-linear optimization scheme, in which we iteratively search for the best-fitting solution by performing several linear slip inversions while varying fault location using a grid search method. Our preferred solution consists of a shallow north-dipping fault plane with assumed azimuth of 282° and accommodating a reverse dextral slip of about 1 m. The estimated geodetic moment is 6.6 × 1018 Nm (Mw 6.5), in agreement with seismological estimates. Yet, our preferred model shows that the geodetic solution is consistent with the activation of a single fault system during the entire sequence, the surface expression of which could be associated with the Buia blind thrust, supporting the hypothesis that the main activity of the Eastern Alps occurs close to the relief margin, as observed in other mountain belts. The retrieved slip pattern consists of a main coseismic patch located 3–5 km depth, in good agreement with the distribution of the main shocks. Additional slip is required in the shallower portions of the fault to reproduce the local uplift observed in the region characterized by Quaternary active folding. We tentatively interpret this patch as postseismic deformation (afterslip) occurring at the edge of the main coseismic patch. Finally, our rupture plane spatially correlates with the area of the locked fault determined from interseismic measurements, supporting the hypothesis that interseismic slip on the creeping dislocation causes strain to accumulate on the shallow (above ∼10 km depth) locked section. Assuming that all the long-term accommodation between Adria and Eurasia is seismically released, a time span of 500–700 years of strain-accumulating plate motion would result in a 1976-like earthquake.

The coseismic ground deformations of the 1997 Umbria-Marche earthquakes: a lesson for the development of new GPS networks

After the occurrence of the two main shocks Mw=5.7 (00.33 GMT) and Mw=6.0 (09:40 GMT) on September 26, 1997, which caused severe damages and ground cracks in a wide area of the Umbria Marche region, the Istituto Nazionale di Geofisica in cooperation with the Istituto Geografico Militare Italiano set out to detect the coseismic ground deformation and reoccupied the available geodetic monuments placed across the epicentral area, belonging to the first order Italian GPS network IGM95 and to the Tyrgeonet network. The comparison between the pre and post-earthquakes coordinate set, the latter obtained from the surveys performed in the early days of October 1997 in the Umbria Marche earthquake area, showed maximum displacements values at the closest stations to the epicentres, up to 14.0±1.8 and 24.0±3.0 cm in the horizontal and vertical components, respectively. The availability of the IGM95 stations allowed geodetic data to be translated into relevant geophysical results. For the first time in Italy, the evaluation of post-earthquake coordinates at 13 vertices provided the estimation of a significant deformation field associated with a seismic sequence. Unfortunately, the same actions could not be applied to the October 14, 1997, Mw=5.6 Sellano earthquake, whose epicentre was located a few tens of km south of the previous ones, due to a lack of available geodetic vertices of Tyrgeonet and IGM95 networks in the surroundings of the epicentral zone. This fact, which prevented the estimation of coseismic deformation and seismic source modelling for this earthquake, clarified the need to set up tailor made GPS networks devoted to geophysical applications, able to capture a possible coseismic signal, but also interseismic and post-seismic signals, at the surface of the Earth’s crust at the scale of the expected magnitudes and fault length. Here we show and discuss the development of the Discrete GPS and Continuous GPS (CGPS) networks in the Italian region started since the early 1990s, which greatly increased after the 1997 Umbria Marche earthquakes, and the insights gained from this action which can be also integrated as Global Observing Strategy to monitor our Environment from Earth and Space.

Co-seismic displacements associated to the Molise (Southern Italy) earthquake sequence of October–November 2002 inferred from GPS measurements

The 2002 earthquake sequence of October 31 and November 1 (main shocks Mw = 5.7) struck an area of the Molise region in Southern Italy. In this paper we analyzed the co-seismic deformation related to the Molise seismic sequence, inferred from GPS data collected before and after the earthquake, that ruptured a rather deep portion of crust releasing a moderate amount of seismic energy with no surface rupture. The GPS data have been reduced using two different processing strategies and softwares (Bernese and GIPSY) to have an increased control over the result accuracy, since the expected surface displacements induced by the Molise earthquake are in the order of the GPS reliability. The surface deformations obtained from the two approaches are statistically equivalent and show a displacement field consistent with the expected deformation mechanism and with no rupture at the surface. In order to relate this observation with the seismic source, an elastic modeling of fault dislocation rupture has been performed using seismological parameters as constraints to the model input and comparing calculated surface displacements with the observed ones. The sum of the seismic moments (8.9 × 10 17 Nm) of the two main events have been used as a constraint for the size and amount of slip on the model fault while its geometry has been constrained using the focal mechanisms and aftershocks locations. Since the two main shocks exhibit the same fault parameters (strike of the plane, dip and co-seismic slip), we modelled a single square fault, size of 15 km × 15 km, assumed to accommodate the whole rupture of both events of the seismic sequence. A vertical E–W trending fault (strike = 266°) has been modeled, with a horizontal slip of 120 mm. Sensitivity tests have been performed to infer the slip distribution at depth. The comparison between GPS observations and displacement vectors predicted by the dislocation model is consistent with a source fault placed between 5 and 20 km of depth with a constant pure right-lateral strike-slip in agreement with fault slip distribution analyses using seismological information. The GPS strain field obtained doesn't require a geodetic moment release larger than the one inferred from the seismological information ruling out significant post-seismic deformation or geodetic deformation released at frequencies not detectable by seismic instruments. The Molise sequence has a critical seismotectonic significance because it occurred in an area where no historical seismicity or seismogenic faults are reported. The focal location of the sequence and the strike-slip kinematics of main shocks allow to distinguish it from the shallower and extensional seismicity of the southern Apennines being more likely related to the decoupling of the southern Adriatic block from the northern one.