Irpinia Region Research Articles

Bayesian Focal Mechanism Estimation from P-, S-Wave Amplitudes, and Polarities for a Microearthquake Sequence in Irpinia, Italy

ABSTRACT The P-wave peak distribution in azimuth and takeoff angles, corrected for the distance through empirical attenuation laws, can reveal the radiation pattern amplitude of the source (Tarantino et al., 2019). This piece of information, jointly combined with the available polarities and S-wave peak amplitudes, can provide the full focal mechanism, constraining the solution also when only a few seismic station records are available, that is, in the case of microseismicity. We proposed a new technique, named P-, S-wave amplitude, and polarities (P-SAP), designed to compute the focal mechanism by jointly inverting the P-, S-wave amplitude ratios and P-wave polarities in a Bayesian framework and employing the octree strategy (Fang et al., 1996) to explore the space of possible solutions. The outputs are strike, dip, and rake angles of the most likely triplet (principal and auxiliary planes) with the related uncertainties, as well as other multiple solutions, if present. We tested the methodology to synthetic data, and we applied it to a microseismic sequence that occurred in Irpinia region, southern Italy. A background microseismicity occurs in a volume delimited by the faults activated during the 1980 M 6.9 Irpinia earthquake. This faults system is complex and composed of northwest–southeast-striking normal faults along the Apennines chain and an approximately east–west-oriented strike-slip fault, deep-seated in the Potenza area. A network of three-component accelerometers and velocimeters are currently deployed to monitor the area. The major topic of this work is validating the proposed P-SAP methodology for microseismicity studies. We also inferred the optimum stress tensor of the sequence, confirming that the microseismicity is controlled by the regional stress field and can reveal characteristics useful to highlight behaviors of larger-scale seismicity.

Regulation of deep carbon degassing by gas-rock-water interactions in a seismic region of Southern Italy

This study is focused on fluids characterization and circulations through the crust of the Irpinia region, an active seismic zone in Southern Italy, that has experienced several high-magnitude earthquakes, including a catastrophic one in 1980 (M = 6.9 Ms). Using isotopic geochemistry and the carbon‑helium system in free and dissolved volatiles in water, this study aims to explore the processes at depth that can alter pristine chemistry of these natural fluids. Gas-rock-water interactions and their impact on CO2 emissions and isotopic composition are evaluated using a multidisciplinary model that integrates geochemistry and regional geological data. By analyzing the He isotopic signature in the natural fluids, the release of mantle-derived He on a regional scale in Southern Italy is verified, along with significant emissions of deep-sourced CO2. The proposed model, supported by geological and geophysical constraints, is based on the interactions between gas, rock, and water within the crust and the degassing of deep-sourced CO2. Furthermore, this study reveals that the Total Dissolved Inorganic Carbon (TDIC) in cold waters results from mixing between a shallow and a deeper carbon endmember that is equilibrated with carbonate lithology. In addition, the geochemical signature of TDIC in thermal carbon-rich water is explained by supplementary secondary processes, including equilibrium fractionation between solid, gas, and aqueous phases, as well as sinks such as mineral precipitation and CO2 degassing. These findings have important implications for developing effective monitoring strategies for crustal fluids in different geological contexts and highlight the critical need to understand gas-water-rock interaction processes that control fluid chemistry at depths that can affect the assessment of the CO2 flux in atmosphere. Finally, this study highlights that the emissions of natural CO2 from the seismically active Irpinia area are up to 4.08·10+9 mol·y−1, which amounts is in the range of worldwide volcanic systems.

Microseismicity clustering and mechanic properties reveal fault segmentation in southern Italy

We analyze the clustered microseismicity recorded over more than one decade in the Irpinia region, an active seismic area in Southern Italy. We studied the spatiotemporal properties of rate, source parameters and slipping mechanisms of the clustered seismicity in relationship with local medium properties and the deformation signal. Spatial pattern of both medium properties and source parameters indicates a segmentation along strike of the fault network with evidence of a high- and a low-coupling segments separated by a region with diffused microseismicity and high b-value. The high-coupling segment shows high fraction of clustered seismicity, low b-value, high relative coseismic slip, low Vp/Vs (1.7–1.75) and near-repeating earthquakes. We find that the seismicity in this high-coupling segment is triggered by non-tectonic loading, such as the recharge of a karst aquifer, suggesting a system close to critical conditions. The high-coupling segment can be identified as a potential strong motion area for future earthquakes, while the intermediate segment can mark the edges of an extensional transfer zone, which can act as a discontinuity for slip propagation. The low-coupling segment is the same where the Ms. 6.9, 1980 Irpinia earthquake nucleated. This segment shows today both relatively high b-values and high coseismic slip, but above all it shows temporal dependencies in its properties, which leads us to pay attention to it for the future spatiotemporal evolution of microseismicity in this area.

One-Year Seismic Survey of the Tectonic CO2-Rich Site of Mefite d'Ansanto (Southern Italy): Preliminary Insights in the Seismic Noise Wavefield.

A passive seismic experiment is carried out at the non-volcanic highly degassing site of Mefite d'Ansanto located at the northern tip of the Irpinia region (southern Italy), where the 1980 MS 6.9 destructive earthquake occurred. Between 2020 and 2021, background seismic noise was recorded by deploying a broadband seismic station and a seismic array composed of seven 1 Hz three-component sensors. Using two different array configurations, we were allowed to explore in detail the 1-20 Hz frequency band of the seismic noise wavefield as well as Rayleigh wave phase velocities in the 400-800 m/s range. Spectral analyses and array techniques were applied to one year of data showing that the frequency content of the signal is very stable in time. High frequency peaks are likely linked to the emission source, whereas at low frequencies seismic noise is clearly correlated to meteorological parameters. The results of this study show that small aperture seismic arrays probe the subsurface of tectonic CO2-rich emission areas and contribute to the understanding of the link between fluid circulation and seismogenesis in seismically active regions.

EASOt-AP: An open-source MATLAB package to estimate the seismic moment, rupture radius, and stress-drop of earthquakes from time-dependent P-wave displacements

We developed a MATLAB package, to rapidly obtain the apparent EArthquake SOurce time function using the Average of near-source P-wave displacement records (EASOt-AP) and to estimate the source parameters of the seismic moment, the rupture radius, and the average static stress drop. The algorithm implemented in this package is based on a rapid and straightforward methodology that is recently developed by Zollo et al. (2021). To this purpose, EASOt-AP models a LPDT curve, which is the average P-wave displacement versus time in the logarithm scale. In this paper, the performance of the EASOt-AP has been shown by evaluating the strong motion data recorded in the near-source range of small magnitude events (2 = Mw ≤ 3) that occurred in the Irpinia region, southern Italy. The source parameters follow a self-similar, constant-stress-drop scaling with a relatively low average stress drop of about 01-0.2 MPa. This Tool is designed for an easy use, and all steps of data processing are performed completely within the numerical environment of MATLAB.

The MS 6.9, 1980 Irpinia Earthquake from the Basement to the Surface: A Review of Tectonic Geomorphology and Geophysical Constraints, and New Data on Postseismic Deformation

The MS 6.9, 1980 Irpinia earthquake occurred in the southern Apennines, a fold and thrust belt that has been undergoing post-orogenic extension since ca. 400 kyr. The strongly anisotropic structure of fold and thrust belts like the Apennines, including late-orogenic low-angle normal faults and inherited Mesozoic extensional features besides gently dipping thrusts, result in a complex, overall layered architecture of the orogenic edifice. Effective decoupling between deep and shallow structural levels of this mountain belt is related to the strong rheological contrast produced by a fluid-saturated, shale-dominated mélange zone interposed between buried autochthonous carbonates—continuous with those exposed in the foreland to the east—and the allochthonous units. The presence of fluid reservoirs below the mélange zone is shown by a high VP/VS ratio—which is a proxy for densely fractured fluid-saturated crustal volumes—recorded by seismic tomography within the buried autochthonous carbonates and the top part of the underlying basement. These crustal volumes, in which background seismicity is remarkably concentrated, are fed by fluids migrating along the major active faults. High pore fluid pressures, decreasing the yield stress, are recorded by low stress-drop values associated with the earthquakes. On the other hand, the mountain belt is characterized by substantial gas flow to the surface, recorded as both distributed soil gas emissions and vigorous gas vents. The accumulation of CO2-brine within a reservoir located at hypocentral depths beneath the Irpinia region is not only interpreted to control a multiyear cyclic behavior of microseismicity, but could also play a role in ground motions detected by space-based geodetic measurements in the postseismic period. The analysis carried out in this study of persistent scatterer interferometry synthetic aperture radar (PS-InSAR) data, covering a timespan ranging from 12 to 30 years after the 1980 mainshock, points out that ground deformation has affected the Irpinia earthquake epicentral area in the last decades. These ground motions could be a result of postseismic afterslip, which is well known to occur over years or even decades after a large mainshock such as the 23 November 1980, MS 6.9 earthquake due to cycles of CO2-brine accumulation at depth and its subsequent release by Mw ≥ 3.5 earthquakes, or most likely by a combination of both. Postseismic afterslip controls geomorphology, topography, and surface deformation in seismically active areas such as that of the present study, characterized by ~M 7 earthquakes. Yet, this process has been largely overlooked in the case of the 1980 Irpinia earthquake, and one of the main aims of this study is to fill such the substantial gap of knowledge for the epicentral area of some of the most destructive earthquakes that have ever occurred in Italy.

Earthquake early warning feasibility in the Campania region (southern Italy) and demonstration system for public school buildings

We present the results of a feasibility study of an earthquake early warning system (EEWS) for the Campania region (southern Italy) using schools as specific targets. The study considered the seismogenic zones as sources of potential earthquakes for the area, the Italian accelerometric network as the recording network for seismic event occurrence, and the performances of the software platform PRESTo Plus for data analysis and processing. We analyze the distribution of lead-times for all possible threatening seismic sources for each municipality in the region under study by extracting the lead-time value corresponding to the 5th, 10th and 25th percentiles of the distributions. We discuss the results for the 5th percentile in order to analyze the worst-case scenario: in the case of a single site, the lead-time is expected to be larger than this value in the 95 % of the cases. Since the population distribution in Campania is uneven and most of the people live nearby the coast, whilst the most destructive earthquakes occur along the Apennine chain, we can conclude that an efficient EEWS can allow most of the schools in the area to undertake some mitigating actions. The testing of the EEWS was carried out in the high school ITIS ‘E. Majorana’, located at Somma Vesuviana, about 80 km from the seismogenic Irpinia region. For this purpose, the Sentinel, an actuator made up of low-cost hardware (i.e., Arduino®), was developed in close cooperation with students and teachers of the school to receive alert messages from the PRESTo Plus platform and warn the school users in case of a seismic event. The EEWS and the Sentinel were successfully tested during some blind drills performed during normal school activities.

Moment Tensor Resolution: Case Study of the Irpinia Seismic Network, Southern Italy

Abstract We present an analysis of the reliability of focal mechanisms obtained through moment tensor inversion at the Irpinia Seismic Network, southern Italy. Our analysis is based on the methods proposed by Zahradnik and Custodio (2012) and Sokos and Zahradnik (2013). We present two different studies: (1) we compute maps of theoretical focal mechanism resolution for the Irpinia region and (2) we study the reliability of the solutions obtained from waveform inversion of five regional earthquakes. Theoretically, we find that when data error is the dominant source of error, focal mechanism resolution is better close to the spots of higher station density rather than at the center of the network. Using real data, we were able to successfully study four of the five regional events, in spite of the large source–station distances (up to ∼280 km) and significant azimuthal gaps (>319°). We used variance reduction, double‐couple percentage, signal‐to‐noise ratio, condition number, and focal mechanism variability index to assess the quality of the solutions. Our quality assessment is validated by comparison with independent focal mechanism solutions.

Comparison between Differential SAR interferometry and ground measurements data in the displacement monitoring of the earth-dam of Conza della Campania (Italy)

This paper presents the application of Differential Synthetic Aperture Radar (SAR) interferometry (DInSAR) algorithms for the precise monitoring of earth dams ground deformation. The test site selected for this study is the Conza Dam, located in the southern Apennines (Italy), very close to the epicenter of the big earthquake (Mw=6.9) which took place on November 23th, 1980, striking the Irpinia region. Among the multiple advanced DInSAR techniques developed by the SAR community during the last decade, this work uses the so-called SUBSOFT software, developed by the Remote Sensing Laboratory (RSLab) group from the Universitat Politècnica de Catalunya (UPC), which is based on the use of Coherent Pixels Technique (CPT) algorithm. The analysis is carried out using 51 ENVIronmental SATellite-Advanced Synthetic Aperture Radar (ENVISAT-ASAR) images, corresponding to the period from 29th of November 2002 to the 30th of July 2010. In this framework, ground displacements recorded by a network of conventional ground-based sensors are also available for the same temporal span. Indeed, the embankment dam is well instrumented to measure internal settlements by means of cross-arms placed in six different cross-sections, and superficial displacements by means of targets for precise leveling. A statistical analysis has been performed to carry out a better comparison between the measurements obtained with the conventional field sensors and the interferometric data. The high agreement between final DInSAR displacements and in-situ instrumental data, demonstrates the reliability of such technique for the precise monitoring of civil infrastructures, and concretely, in dams with a high exposure factor and its consequent risk.

Seismic Retrofitting Techniques for Existing Masonry Buildings

In many seismically active regions of the world there are large numbers of masonry buildings. Most of these buildings have not been designed for seismic loads. Recent earthquakes have shown that many of these buildings are seismically vulnerable and should be considered for retrofitting. Different conventional and unconventional retrofitting techniques are available to increase the strength and/or ductility of unreinforced masonry (URM) walls. This paper reviews and discusses seismic retrofitting of masonry walls with emphasis on the conventional techniques. Retrofitting procedures are discussed with regard to a case study: a stone masonry building in Irpinia region, damaged by the 1980 earthquake. The interventions are evaluated by means of finite elements with a macroelement obtained with a homogenization technique. Linear and nonlinear procedures are compared, and peculiarities of each procedure are shown.

Effects of Surface Geology on Seismic Ground Motion Deduced from Ambient-Noise Measurements in the Town of Avellino, Irpinia Region (Italy)

The effects of surface geology on ground motion provide an important tool in seismic hazard studies. It is well known that the presence of soft sediments can cause amplification of the ground motion at the surface, particularly when there is a sharp impedance contrast at shallow depth. The town of Avellino is located in an area characterised by high seismicity in Italy, about 30 km from the epicentre of the 23 November 1980, Irpinia earthquake (M = 6.9). No earthquake recordings are available in the area. The local geology is characterised by strong heterogeneity, with impedance contrasts at depth. We present the results from seismic noise measurements carried out in the urban area of Avellino to evaluate the effects of local geology on the seismic ground motion. We computed the horizontal-to-vertical (H/V) noise spectral ratios at 16 selected sites in this urban area for which drilling data are available within the first 40 m of depth. A Rayleigh wave inversion technique using the peak frequencies of the noise H/V spectral ratios is then presented for estimating Vs models, assuming that the thicknesses of the shallow soil layers are known. The results show a good correspondence between experimental and theoretical peak frequencies, which are interpreted in terms of sediment resonance. For one site, which is characterised by a broad peak in the horizontal-to-vertical spectral-ratio curve, simple one-dimensional modelling is not representative of the resonance effects. Consistent variations in peak amplitudes are seen among the sites. A site classification based on shear-wave velocity characteristics, in terms of Vs30, cannot explain these data. The differences observed are better correlated to the impedance contrast between the sediments and basement. A more detailed investigation of the physical parameters of the subsoil structure, together with earthquake data, are desirable for future research, to confirm these data in terms of site response.

A study to constrain the geometry of an active fault in southern Italy through borehole breakouts and downhole logs

Identification of an active fault and the local versus regional present-day stress field in the Irpinia region (southern Apennines) have been performed along a 5900 m deep well (San Gregorio Magno 1) by a detailed breakout and geophysical log analysis. The selected area is characterized by diffuse low magnitude seismicity, although in historical times moderate to large earthquakes have repeatedly struck it. On 23rd November 1980 a strong earthquake ( M = 6.9) nucleated on a 38-km long normal fault, named Irpinia fault, producing the first unequivocal historical surface faulting ever documented in Italy. The analysis of stress-induced wellbore breakouts shows a direction of minimum horizontal stress N18°±24°, fairly consistent with the regional stress trend (N44°±20°). The small discrepancy between our result and the regional stress orientation might be related to the influence of local stress sources such as variations of the Irpinia fault plane orientation and the presence of differently oriented active shear zones. This paper shows for the first time a detailed analysis on the present-day stress along a well to identify the Irpinia fault at depth and constrain its geometry.

Inferences on the source mechanisms of the 1930 Irpinia (Southern Italy) earthquake from simulations of the kinematic rupture process

We examine here a number of parameters that define the source of the earthquake that occurred on 23rd July 1930 in Southern Italy (in the Irpinia region). Starting from the source models proposed in different studies, we have simulated the acceleration field for each hypothesized model, and compared it with the macroseismic data. We then used the hybrid stochastic-deterministic technique proposed by Zollo et al. (1997) for the simulation of the ground motion associated with the rupture of an extended fault. The accelerations simulated for several sites were associated with the intensities using the empirical relationship proposed by Trifunac and Brady (1975), before being compared with the available data from the macroseismic catalogue. A good reproduction of the macroseismic field is provided by a normal fault striking in Apenninic direction (approximately NW-SE) and dipping 55° toward the SW.

Earthquake focal mechanisms and stress inversion in the Irpinia Region (southern Italy)

The goal of this study was to estimate the stress field acting in the Irpinia Region, an area of southern Italy that has been struck in the past by destructive earthquakes and that is now characterized by low to moderate seismicity. The dataset are records of 2,352 aftershocks following the last strong event: the 23 November 1980 earthquake (M 6.9). The earthquakes were recorded at seven seismic stations, on average, and have been located using a three-dimensional (3D) P-wave velocity model and a probabilistic, non-linear, global search technique. The use of a 3D velocity model yielded a more stable estimation of take-off angles, a crucial parameter for focal mechanism computation. The earthquake focal mechanisms were computed from the P-wave first-motion polarity data using the FPFIT algorithm. Fault plane solutions show mostly normal component faulting (pure normal fault and normal fault with a strike-slip component). Only some fault plane solutions show strike-slip and reverse faulting. The stress field is estimated using the method proposed by Michael (J Geophys Res 92:357–368, 1987a) by inverting selected focal mechanisms, and the results show that the Irpinia Region is subjected to a NE–SW extension with horizontal σ3 (plunge 0°, trend 230°) and subvertical σ1 (plunge 80°, trend 320°), in agreement with the results derived from other stress indicators.

An integrated geophysical investigation of the upper crust in the epicentral area of the 1980, Ms=6.9, Irpinia earthquake (Southern Italy)

In this paper, we investigate the upper crustal structure of the Irpinia region, Southern Apennines thrust belt, Italy, through analysis and joint interpretation of gravity data, seismic reflection lines and subsurface information from many deep wells. The investigated region includes the epicentral area of the 1980 ( M s=6.9) Irpinia earthquake and is one of the Italian regions with the highest seismic hazard. The upper crustal structure is imaged by modeling a series of 15 SW-trending gravity profiles, spaced about 5 km apart, plentifully constrained by seismic reflection lines and wells, thus reducing the inherent ambiguity of the gravity modeling. Despite of the complexity of the modeled Bouguer anomalies, the application of a calibrating procedure to constrain the range of variability of the density values, as well as the use of geometric constraints, results in a good level of stability in the final density cross-sections, which in fact appear coherent both in the density values and in the geometrical features. The inferred model shows important lateral density variations that can be mostly related to NW-trending geologic structures. High-density bodies delineate carbonate platform thrust sheets and broad antiforms involving Mesozoic basinal rocks, while low-density shallow bodies are associated with Pliocene basins. In addition, important density (i.e. lithological) variations are evident along the strike of the range, the most relevant being an abrupt deepening of the Apulia Carbonate Platform in the southeastern part of the investigated region. In the epicentral region of the 1980 event, we find that the geometry of the high-density, high-velocity carbonates of the Apulia Platform appears correlated with the distribution of the aftershocks and with the P-wave velocity anomaly pattern as inferred from a previous local earthquake tomography. The structural highs of the Apulia Platform correspond to high-velocity regions, where aftershocks and coseismic slip of the mainshock are concentrated. This correlation suggests that the Apulia Carbonate Platform geometry played an important role in the rupture propagation and in the aftershock distribution.

The Deep Seismicity of the Tyrrhenian Sea*

AbstractThe deep seismicity of the Tyrrhenian Sea is analysed using data from a new instrumental catalogue of the seismicity of the Italian area. We use algorithms for the determination of absolute and relative hypocentral locations and for the evaluation of the geometry and coherence of the state of stress within the subducting slab. The structure of the Benioff zone, although simpler than previously indicated, reveals anomalous traits both in the seismicity distribution and in the stress geometry, confirming that standard subduction models cannot be applied in the Tyrrhenian region. The velocity anomaly and the location of few isolated events indicate that the subducted slabs extend to the north along the Apenninic chain approximately to the latitude of the Irpinia region, in Central Italy (˜ 42°N).