Search

The earthquake Magnitude-Frequency-Distribution, FMD, is usually modelled with the Gutenberg-Richter relation law, where the b-value controls the relative rate of small and large earthquakes. b-value has been documented to show an inverse dependence on differential stress, it increases with the fault roughness or during fluid-induced earthquakes. For some seismic sequences a near real-time characterization of the b-value has been used to discriminate between foreshocks and aftershocks. Here we examine the influence on b-values of different lithologies hosting earthquakes.In general, seismicity not only localizes along the major structures where mainshocks nucleate, but it can be also distributed within volumes of the seismogenic layer characterized by different lithologies. For the Mw 6.5 2016–2017 Central Italy seismic sequence, the lithology can be properly defined by seismic reflection profiles. Here the fractured carbonate of the Apennines, located at almost 1-2 km and 4-6 km of depth, are characterized by b-values ranging between 1.3 and 1.4 that can be diagnostic of brittle dominated deformation. At 2-4 km and 6-10 km of depth, the Triassic Evaporites showing a bimodal brittle-ductile deformation and compartmentalized fluid overpressure (documented in deep boreholes) are linked to high b-values, in the range of 1.5-1.65 reaching 1.80 for clustered swarms. Between 10-12 km of depth the phyllosilicate rich basement, with its predominant velocity strengthening behaviour, is hosting small magnitude earthquakes with b-values around 1.4. Our results indicate that away from the large earthquake faults, characterized by a stress dependent elasto-frictional rheology, FMD are strongly controlled by rock lithology and style of deformation.

In the last 15 years, activities for geo-energy production are associated to subsurface fluid injection in enhanced geothermal systems, for enhanced oil recovery, for the disposal of wastewater or for carbon dioxide capture and storage. In several regions, M>3 earthquakes occurred following fluid injection, and some of these earthquakes have caused extensive damage, putting geo-energy production projects at risk of being discontinued. Evaluating the conditions under which fluid injection can induce earthquakes is therefore important to safeguard local infrastructures and to ensure continuity of geo-energy projects. To shed light on the effect of fluid injection on a fault located in the proximity of a reservoir, we implemented into the Q-DYN seismic cycle simulator the fluid diffusion equation (one-way coupling). We ran models of seismic cycles on a rate-and-state-dependent fault under a quasi-dynamic approximation, and we developed a systematic study to assess how fault frictional heterogeneities, the stress state of the fault upon injection, the timing of injection relative to the phase of the seismic cycle and factors controlling fluid flow, i.e. permeability, porosity, flow-rate, influence fault slip behavior and earthquake magnitude. Our results show that localized pore-pressure perturbations allow us to gain deeper physical insight into the propagation and arrest of earthquake ruptures and that changes in the fault physical properties can promote a spectrum of fault slip behavior and recorded magnitudes.

Serpentinites are poly-mineralic rocks distributed almost ubiquitously in active tectonic regions worldwide. They are composed of rheologically weak (lizardite and crysotile) and strong (e.g., magnetite and pyroxene) phases. In particular, lizardite typically shows low friction coefficients and is supposed to localise deformation along weak shear zones characterized by aseismic behaviour. Major faults hosting serpentinite lithologies are characterised by seismic activity, tremors, and other slip modes. We advance the hypothesis that low strain domains, which are enriched in rheologically strong phases, can act as potential site of nucleation of unstable slip as the result of the velocity-dependent rheology of magnetite-rich serpentinites. Through an experimental and microstructural approach, we explore the different mechanisms whose interplay controls the complex behaviour of these lithologies.For this study we collected natural samples of lizardite-magnetite rich serpentinites within the low strain domains of the Elba Island ophiolites (Italy). Rocks were characterised, powdered, and deformed in a set of shear experiments at four different normal stresses (25, 50, 75 and 100 MPa) in the biaxial apparatus BRAVA. The experiments consist of an initial phase of sliding at 10 μm/s, a slide-hold-slide test, and two series of velocity stepping (sliding velocity from 0.1 to 300 μm/s). Fundamental parameters to quantify the frictional properties of serpentinites are individuated in the (a-b) value, the critical slip distance Dc, and the critical stiffness kc, which is derived by their combination.The material shows friction values of ~0.4 with velocity weakening behaviour and negative frictional healing. The module of the negative (a-b) parameter increases neatly with decreasing sliding velocity while Dc decreases, causing kc to rise. At low velocities (< 3 μm/s) sliding is unstable and the fault undergoes stick-slip behaviour. This is explained by the increase of the critical stiffness to values higher than the loading system stiffness. This systematic change of mechanical properties and fault slip behaviours with sliding velocity is interpreted to be the result of the time-dependent arrangement of grains in a heterogeneous experimental fault architecture.Back-scattered SEM images of the principal slip zones of recovered samples support this hypothesis. Elasto-frictional behaviour is controlled by the build-up of a partial (granular) load-bearing framework of strong magnetite grains, while visco-frictional rheology is controlled by the (phyllosilicatic) anastomosing and foliated lizardite matrix. At low sliding velocities, the granular phase interacts creating force chains thus promoting frictional instabilities. At higher velocities, dilation promotes the activity of throughgoing weaker phyllosilicate planes thus favouring stable slip.Our experiments shed light on the role of fault rock heterogeneity in nucleating dynamic slip in nature as well as in controlling the slip mode during earthquakes or slow-slip events in serpentinite terrains.

After many large earthquakes aftershocks activity can reach depths greater than the base of the seismogenic zone that is defined by background seismic activity. This observation is generally explained by strain rate induced embrittlement associated with the increase of post-mainshock strain rate, which favors a transition from ductile to brittle behavior. However, the underlying physical processes are not well understood. Here we integrate geological and geophysical data for the 2016–2017 Central Italy seismic sequence with laboratory experiments to provide a geological and physical interpretation for the post-mainshock transient deepening of the base of the seismogenic zone.The base of the seismogenic zone in the central-northern Apennines is set typically at 9-10 kilometers and corresponds to the top of the phyllitic basement. Structural studies on exhumed basement rocks highlight a heterogeneous basement fabric consisting of competent, 10 to 200 m wide, quartz-rich lenses surrounded by an interconnected and frictionally weak phyllosilicate-rich matrix. The matrix controls the overall rheology of the basement due to its interconnectivity, and promotes aseismic deformation because its rate-strengthening behavior.Following each major (Mw > 5.5) event of the 2016–2017 sequence, a dramatic and abrupt increase in seismic rate is observed below 10 km, hence within the basement. Here we document the presence of seismicity clusters made of more than 4 earthquakes and characterized by small magnitude (Mw < 2.5), small dimensions (< 500 meters), small temporal duration (< 14 days) and a swarm-like behavior. Furthermore, these clusters are often represented by multiple or repeating earthquakes with a cross correlation coefficient higher than 0.7 for all the three components. These observations suggest that the increase of shear stressing rate within the basement is responsible for deepening of seismicity. To further explore this idea, we performed laboratory experiments on rocks from exhumed outcrops of basement rocks. We found that shear stressing rate promotes accelerated creep on the phyllosilicate-rich matrix and dynamic instabilities on the quartz-rich gouge belonging to the lenses. Our integrated analysis suggests that the mainshocks of the 2016-2017 seismic sequence promoted an increase of shear stressing rate within the basement allowing the phyllosilicate-rich matrix to creep faster hence favoring the loading and the repeated failures of locked seismogenic patches represented by the quartz-rich lenses.

Accurate assessment of rate and state friction parameters is essential for producing realistic rupture scenarios and, in turn, for seismic hazard analysis. Those parameters can be directly measured in the laboratory, with experimental apparati that reproduce fault conditions in nature. Alternatively, indirect estimates (i.e., inversion) of rate and state parameters are based on postseismic slip evolution studies and numerical modeling. Both direct and indirect approaches require a series of assumptions that might bias the results.Here we take advantage of a downscaled analog model reproducing experimentally megathrust earthquakes. The analog model shares many characteristics of real subduction zones, although being intentionally oversimplified with respect to nature. This allows reducing the number of potential sources of bias (e.g., fault geometry and asperity size). We perform five analog models with a single, rectangular asperity of different lengths embedded in a nearly velocity neutral matrix. We focus on two different physical conditions, namely the along-strike asperity length and the asperity to neighboring segments length ratio, and study systematically how they tune the model seismic behavior. Then, by coupling quasi-dynamic numerical models with the simulated annealing algorithm, we retrieve rate and state parameters that allow reproducing both the recurrence time, rupture duration and slip amplitude of the analog model, in ensemble. We identify a tradeoff between (a-b) of the asperity and (a-b) of neighboring creeping segments, with multiple combinations that allow mimicking the analog model behavior and variability. We also identify a negative correlation between (a-b) of the asperity and asperity size, with Dc remaining relatively constant within the investigated asperity size range. When estimating (a-b), poorly constrained properties of neighboring segments are responsible for uncertainties in the order of per mille. Roughly one order of magnitude larger uncertainties derive from asperity size. Those results provide a first order assessment of the variability that rate and state friction estimates retrieved for nature conditions might have when used as constraint to model fault slip behavior.

In the last decades, it has been observed that faults can slip both by slow aseismic creep and seismic events (i.e., earthquakes). Between these two slip modes, a wide variety of fault slip behavior can be observed, including low-frequency earthquakes, slow slip events and tremors. This wide variety of slip modes can radiate seismic energy at different frequencies whose content may be linked to the physical mechanisms at play. In the laboratory, it is possible to reproduce the entire spectrum of fault slip modes by modulating the loading stiffness of the apparatus depending on the critical fault rheologic stiffness (i.e. k/kc). This technique allows us to study, under controlled laboratory conditions, the acoustic signature of different fault slip modes to infer the physical mechanisms at their origin. To shed light on the nucleation mechanisms and seek for reliable precursors to failure of different slip modes, we performed friction experiments on powders that differ for granulometry and grain shape (i.e., glass beads with a grainsize < 150 µm; and quartz powders MinUSil with an average grain size of 10.5 µm), to simulate fault gouge. The experiments were conducted in a double direct shear configuration, instrumented with an array of piezoelectric sensors to record continuously Acoustic Emissions (AEs) at high recording rate (~10MHz). The experiments are performed at a constant displacement rate of 10 µm/s and using a spring to reduce the apparatus stiffness k, to match the critical fault rheological stiffness (kc). Following this procedure we obtain slow slip events (i.e., k = kc) and fast events (i.e. k

Seismic waves contain information about the earthquake (EQ) source and many forms of noise deriving from the seismometer, anthropogenic effects, background noise associated with ocean waves, and microseismic noise. Separating the noise from the EQ signal is a critical first step in EQ physics and seismic waveform analysis. However, this is difficult because optimal parameters for filtering noise typically vary with time and may strongly alter the shape of the waveform. A few recent works have employed Deep Learning (DL) model for seismic denoising, among which we have taken as a benchmark Deep Denoiser and SEDENOSS. These models turn the noisy trace into a 2D signal (spectrograms) within the model to denoise the traces, making the process pretty heavy. We propose a novel DL-powered seismic denoising algorithm based on Diffusion Models (DMs), keeping the signal in 1D. DMs are the latest trend in Machine Learning (ML), having revolutionized the application fields of audio and image processing for denoising (DiffWave), synthesis (Stable Diffusion), and sequence modeling (STARS). The training of DMs proceeds by polluting a signal with noise until the signal has completely vanished into noise, then reversing the process by iterative denoising, conditioned on the latent signal representation. This makes DMs the ideal tool for seismic traces cleaning, as the model naturally learns from seismic sequences by denoising, which aligns the ML training procedure and the final task objective. In a preliminary evaluation, we used the Stanford Earthquake Dataset (STEAD); our proposed Diffusion-based Seismic Denoiser (DiffSD) outperforms the state-of-the-art DL methods on the Signal Noise Ratio (SNR), Scale-Invariant Source to Distortion Ratio (SI-SDR), and Source to Distortion Ratio (SDR) metrics. DiffSD also yields qualitatively pleasing EQ traces out of visual inspection in time and frequency. Finally, DiffSD proceeds from regenerating clean EQ signals from noise, which opens the way to data-driven EQ sequence generations, potentially instrumental to further study and dataset augmentations.

Geological observations reveal that earthquakes nucleate, propagate, and arrest in complex fault zones whose structural heterogeneity depends on the tectonic loading, geometry, lithology, rheology, presence of fluids, and strain localization processes. These fault zones can host a wide range of fault slip behaviors (e.g., creep, aseismic- and slow-slip events, afterslip, and earthquakes). This implies that the environment in which earthquakes occur is diverse, and that different physical and chemical processes can be involved during the coseismic dynamic rupture.Earthquakes are generated by rupture propagation and slip within fault cores and dissipate the stored elastic and gravitational strain energy in fracture and frictional processes in the fault zone (from microscale - less than a millimeter - to macroscale - centimeters to kilometers) and in radiated seismic waves. Understanding this energy partitioning is fundamental in earthquake mechanics to describe dynamic fault weakening and causative rupture processes operating over different spatial and temporal scales.The energy dissipated in earthquake rupture propagation is commonly called fracture energy (G) or breakdown work (Wb). Here we discuss these two parameters, and we review fracture energy estimates from seismological, modeling, geological, and experimental studies and show that fracture energy scales with fault slip and earthquake size. Our analysis confirms that seismological estimates of fracture energy and breakdown work are comparable and scale with seismic slip. The inferred scaling laws show modest deviations explained in terms of epistemic uncertainties. The original collection of fracture energy estimates from laboratory experiments confirms the scaling with slip over a slip range of more than 10 decades. Fracture energy associated with breaking of intact rocks is larger than for precut specimens and might suggest differences between the role of fracture and friction, or a different size of the rupture front zone. It is important to recall that fault products after deformation in the laboratory correspond to fault products observed in nature, and acoustic emissions recorded in the laboratory can be processed as seismic waves on a natural fault. We conclude that although material-dependent constant fracture energies are important at the microscale for fracturing grains of the fault zone, they are negligible with respect to the macroscale processes governing rupture propagation on natural faults.In this study we discuss the scaling of fracture energy and breakdown work with slip, and we propose different interpretations relying on different processes characterizing complex fault zones. Our results suggest that, for earthquake ruptures in natural faults, the estimates of G and Wb are consistent with a macroscale description of the causative processes.Reconciling observations and results from laboratory experiments and numerical modeling with geological observations can be done, provided that we accept the evidence that earthquakes can occur in a variety of geological settings and fault zone structures governed by different physical and chemical processes.

Understanding the dynamics of microearthquakes is a timely challenge to solve current paradoxes in earthquake mechanics, such as the stress drop and fracture energy scaling with seismic moment. Dynamic modeling of microearthquakes induced by fluid injection is also relevant for studying rupture propagation following a stimulated nucleation. We study the main features of unstable dynamic ruptures caused by fluid injection on a target preexisting fault (50m x 50m) generating a Mw=1 event. The selected fault is located in the Bedretto Underground Laboratory (Swiss Alps) at ≈1000m depth. We perform fully dynamic rupture simulations and model seismic wave propagation in 3D by adopting a linear slip-weakening law. We use the distributed multi-GPU implementation of SeisSol on the supercomputer Marconi100. Stress field and fault geometry are well constrained by in-situ observations, allowing us to minimize the a priori imposed parameters. We investigate the scaling relations of stress drop, slip-weakening distance (Dc) and fracture energy (Gc) focusing on their role in governing dynamics of rupture propagation and arrest for a target Mw=1 induced earthquake. We explore different homogenous conditions of frictional parameters, and we show that the spontaneous arrest of the rupture is possible in the modeled stress regime, by assuming a high ratio between stress excess and dynamic stress drop (the fault strength parameter S), characterizing the fault before the fluid pressure change. The rupture arrest of modeled induced earthquakes depends on the heterogeneity of dynamic parameters due to the spatially variable effective normal stress. Moreover, for a fault with high S values (not ready to slip), small variations of Dc (≈0.5÷1.2 mm) can drive the rupture from self-arrested to run-away. Studying dynamic interactions (stress transfer) among slipping points on the rupturing fault provides insight on the breakdown process zone and shear stress evolution at the crack tip leading to failure. The inferred spatial dimension of the cohesive zone in our models is nearly ~0.3-0.4m, with a maximum slip of ~0.6 cm. Finally, we compare stress drop and fracture energy estimated from synthetic waveforms with assumed dynamic parameters. Our results suggest that meso-scale processes near the crack-tip affect rupture dynamics of micro-earthquakes.

The Maltese archipelago lies in the centre of the Sicily Channel (Central Mediterranean), a crustal domain subject to extension since the Miocene and characterised by the presence of several active graben systems in the Plio-Quaternary. This region is affected by diffuse seismicity, in the form of isolated swarms interspersed with long pauses of quiescence and events of modest magnitude, in general not exceeding 4.0. The key features of seismicity are poorly constrained by unfavourable geographical conditions and the scarcity of seismic stations in the area between southern Sicily and Tunisia. Until a few years ago, the only Maltese station operating (WDD of the MedNet network) was the only source of information to characterise local Maltese seismicity. On the other hand, the localisations of major events (M > 4) made by the other monitoring agencies (e.g. Istituto nazionale di Geofisica e Vulcanologia - INGV) are based on traveltime readings from seismograms recorded at great distances from Malta (> 100 km) and are therefore affected by considerable formal errors. Moreover, the only published catalogue of offshore seismicity around the Maltese islands (Seismic Monitoring and Research Group, University of Malta) is based on single station location at WDD up to 2014. This is characterised by considerable uncertainties in hypocentral parameters, in particular a lack of depth information. In consequence, therefore, it is very difficult to obtain information on the geometry and kinematics of active tectonic structures, both offshore and possibly on.In recent years, there has been a considerable increase in the number of seismic stations on the Maltese archipelago: there are currently eight active stations, that make up the Malta Seismic Network (MSN). The study we present is an example of the application of 3-D localisation techniques of Maltese seismicity, which benefits from the recent implementation of the local seismic network. In particular, we focus on the swarm occurring offshore in the period September-October 2020, and characterised by an unusual number of events (> 100), including a main event (M > 4.0) that was strongly felt over the archipelago. We have handpicked P- and S-wave traveltimes for all the events using recordings from the MSN and a selected number of Italian stations according to quality criteria, and inverted them using the Hypoellipse code. The precise earthquake localisations allow us to obtain details of the structures activated during the swarm, together with a more detailed insight into the time evolution of the sequence. Accompanying these analyses, we calculated the moment tensor solutions for some of the largest events (M > 3) and performed a spectral analysis to distinguish the waveform characteristics of events occurring at shallow depths (< 10 km) from those nucleating in the mid-crust (> 15 km depth), as a function of different wavepaths through the crust. The prospect is that in the future it will be possible to better constrain the seismotectonics of the Maltese archipelago and have a more accurate picture of the seismogenic potential of active faults in this sector of the Mediterranean.

Fault zone properties can change significantly during the seismic cycle in response to stress changes, microcracking and wall rock damage. Lab experiments show consistent changes in elastic properties prior to and after lab earthquakes (EQ) and previous works show that machine learning/deep learning (ML/DL) techniques are successful for capturing such changes. Here, we apply DL techniques to assess whether similar changes occur during the seismic cycle of tectonic EQ. The main motivation is to generalize lab-based findings to tectonic faulting, to predict failure and identify precursors. The novelty is that we use EQ traces as probing signals to estimate the fault state.We train DL model to distinguish foreshocks, aftershocks and time to failure of the Mw 6.5 2016 Norcia EQ in central Italy, October 30th 2016. We analyze a 25-second window of 3-component data around the P- and S-wave arrivals for events near the Norcia fault with M>0.5 and ±2 months before/after the Norcia mainshock. Normalized waveforms are used to train a Convolutional Neural Network (CNN). As a first task we divide events into two classes (foreshocks/aftershocks), and then refine the classification as a function of time-to-failure (TTF) for the mainshock. Our DL model perform very well for TTF classification into 2, 4, 8, or 9-classes for the 2 months before/after the mainshock. We explore a range of seismic ray paths near, through, and away from the Norcia mainshock fault zone. Model performance exceeds 90% for most stations. Waveform investigations show that wave amplitude is not the key factor; other waveform properties dictate model performance. Models derived from seismic spectra, rather than time-domain data, are equally good. We challenged the model in several ways to confirm the results. We found reduced performance in training the model with the wrong mainshock time and by omitting data immediately before/after the mainshock. Foreshock/aftershock identification is significantly degraded also by removing high frequencies (filtering seismic data above 25 Hz). We tested data from different years to understand seasonality at individual stations for the time period September to December and removed these effects. Comparing these seasonality effects defined from noise with our EQ results shows that foreshocks/aftershocks for the 2016 Norcia mainshock are well resolved. Training with data containing EQ offers a huge increase in classification performance over noise only, proving that EQ signals are the sole that enable assessing timing as a function of the fault status. To confirm our results and understand which stations are able to detect changes of fault properties we perform a further test cleaning the signals from the seasonality by confounding the DL with a shuffled noise (adversarial training).We conclude that DL is able to recognize variations in the stress state and fracture during the seismic cycle. The model uses EQ-induced changes in seismic attenuation to distinguish foreshocks from aftershocks and time to failure. This is an important step in ongoing efforts to improve EQ prediction and precursor identification through the use of ML and DL.

Tectonic faults are often assumed to slip either slow due to stable, velocity-strengthening frictional behavior, or fast as a result of velocity-weakening friction leading to dynamic (seismic) rupture. As a consequence, velocity-strengthening faults may be regarded as intrinsically stable as they do not spontaneously nucleate seismic events. However, recent laboratory and in-situ experiments of fluid injection challenged such assumptions. Here we present a fully coupled hydro-mechanical fault model in which stable, rate-strengthening frictional behavior is combined with dynamic weakening due to rapid poroelastic effects, allowing unstable (seismic) slip to occur on nominally stable faults. In our numerical experiments, fluid injection reduces the effective normal stress and frictional resistance, thus bringing the fault to failure. The onset of fault failure is controlled by competing mechanisms of shear-induced dilation and shear-enhanced compaction, which cause fault weakening and the propagation of a slow-slip transient from the fluid injection point. When a critical size of the slow slip patch is reached, dynamic rupture eventually nucleates at the slow-slip event front and propagates beyond the fluid pressure perturbed region. Further numerical experiments indicate that, when the fault is critically stressed, the growth of the aseismic patch – prior to dynamic rupture – occurs in a few seconds, whereas at lower stress levels, the aseismic slip phase propagates slowly over hundreds of seconds. These results predict that poroelastic compaction and fluid pressurization can cause the transition from aseismic slow-slip to fast seismic slip and the propagation of dynamic rupture on velocity-strengthening faults. In particular, they demonstrate that compaction-induced fluid pressurization can overcome the initial phase of shear-induced dilatancy, thus allowing the propagation of dynamic rupture in the form of pulse-like pore-pressure waves. The implication that earthquake rupture may nucleate on rate-strengthening faults, presently considered to be nominally stable, requires a re-evaluation of seismic hazard in many areas, particularly in the case of fluid injection in enhanced geothermal systems and CO2 storage.

Analysis of seismicity can illuminate active fault zone structures but also pervasive deformation occurring within large volumes of the seismogenic zone. Here we show that, for the Mw 6.5 2016–2017 Central Italy seismic sequence, seismicity not only localizes along SW-dipping normal faults, but also occurs within larger volumes of Triassic Evaporites, TE, composed of alternated anhydrites and dolostones. This off-fault and diffuse microseismicity shows a different frequency-magnitude distribution than on-fault seismicity along the structures hosting the largest events. We interpret that, during the sequence, shear strain-rate increase promoted widespread ductile deformation within TE that light-up with diffuse microseismicity. This interpretation is supported by field and laboratory observations showing that TE background ductile deformation is complex and dominated by distributed failure and folding of the anhydrites associated with boudinage fracturing and faulting of dolostones. Our results indicate that ductile crustal deformation can cause off-fault and diffuse microseismicity, which obeys to different scaling laws than on-fault seismicity on structures characterized by elasto-frictional stick-sip behaviour.

The 2016 Central Italy earthquake sequence is characterized by remarkable rupture complexity, including highly heterogeneous slip across multiple faults in an extensional tectonic regime. The dense coverage and high quality of geodetic and seismic data allow to image intriguing details of the rupture kinematics of the largest earthquake of the sequence, the Mw 6.5 October 30th, 2016 Norcia earthquake, such as an energetically weak nucleation phase. Several kinematic models suggest multiple fault planes rupturing simultaneously, however, the mechanical viability of such models is not guaranteed.Using 3D dynamic rupture and seismic wave propagation simulations accounting for two fault planes, we constrain 'families' of spontaneous dynamic models informed by a high-resolution kinematic rupture model of the earthquake. These families differ in their parameterization of initial heterogeneous shear stress and strength in the framework of linear slip weakening friction.First, we dynamically validate the kinematically inferred two-fault geometry and rake inferences with models based on only depth-dependent stress and constant friction coefficients. Then, more complex models with spatially heterogeneous dynamic parameters allow us to retrieve slip distributions similar to the target kinematic model and yield good agreement with seismic and geodetic observations. We discuss the consistency of the assumed constant or heterogeneous static and dynamic friction coefficients with mechanical properties of rocks at 3-10 km depth characterizing the Italian Central Apennines and their local geological and lithological implications. We suggest that suites of well-fitting dynamic rupture models belonging to the same family generally exist and can be derived by exploiting the trade-offs between dynamic parameters.Our approach will be applicable to validate the viability of kinematic models and classify spontaneous dynamic rupture scenarios that match seismic and geodetic observations at the same time as geological constraints.

Moment tensor inversions of broadband velocity data are usually managed by adopting Green’s functions for 1D layered seismic wavespeed models. This assumption can impact on source parameter estimat...

SUMMARYOn 24 August 2016 at 01:36 UTC a ML6.0 earthquake struck several villages in central Italy, among which Accumoli, Amatrice and Arquata del Tronto. The earthquake was recorded by about 350 seismic stations, causing 299 fatalities and damage with macroseismic intensities up to 11. The maximum acceleration was observed at Amatrice station (AMT) reaching 916 cm s–2 on E–W component, with epicentral distance of 15 km and Joyner and Boore distance to the fault surface (RJB) of less than a kilometre. Motivated by the high levels of observed ground motion and damage, we generate broad-band seismograms for engineering purposes by adopting a hybrid method. To infer the low frequency seismograms, we considered the kinematic slip model by Tinti et al . The high frequency seismograms were produced using a stochastic finite-fault model approach based on dynamic corner-frequency. Broad-band synthetic time-series were therefore obtained by merging the low and high frequency seismograms. Simulated hybrid ground motions were compared both with the observed ground motions and the ground-motion prediction equations (GMPEs), to explore their performance and to retrieve the region-specific parameters endorsed for the simulations. In the near-fault area we observed that hybrid simulations have a higher capability to detect near source effects and to reproduce the source complexity than the use of GMPEs. Indeed, the general good consistency found between synthetic and observed ground motion (both in the time and frequency domain), suggests that the use of regional-specific source scaling and attenuation parameters together with the source complexity in hybrid simulations improves ground motion estimations. To include the site effect in stochastic simulations at selected stations, we tested the use of amplification curves derived from HVRSs (horizontal-to-vertical response spectra) and from HVSRs (horizontal-to-vertical spectral ratios) rather than the use of generic curves according to NTC18 Italian seismic design code. We generally found a further reduction of residuals between observed and simulated both in terms of time histories and spectra.

During the 2016–2017 Central Italy earthquake sequence, a series of moderate to large earthquakes M > 5 occurred near the Amatrice and Norcia towns. These events are recorded on a dense seismic network, providing relevant observational evidence of complex earthquakes in time and space. In this work, we used this substantial data set to study the ground-motion characteristics of the Norcia earthquake M6.5 on October 30, 2016, through a broadband ground-motion simulation. Three-component broadband seismograms are generated to cover the entire frequency band of engineering interest. Low and high frequencies are computed considering the heterogeneous slip rupture model of Scognamiglio et al. (2018) [1]. High frequencies are calculated using a stochastic approach including P, SV, and SH waves, while low frequencies are obtained through a forward simulation of the kinematic model at the various stations. To predict earthquake-induced ground-motions in the area, we adopted region-specific attenuation and source scaling parameters derived by Malagnini et al. (2011) [2]. Ground-motion parameters, including peak ground acceleration (PGA), peak ground velocity (PGV) and spectral amplitudes, are calculated at the selected sites adopting physics-based parameters to understand better the earthquake fault rupture, the wave propagation, and their impacts on the seismic hazard assessment in the region. We showed that combining the fault rupture history over the entire frequency spectrum of engineering interest, the attenuation characteristics of the seismic wave propagation and the properly defined site responses can improve the prediction of ground-motions and time histories, especially in near seismic sources.

The 2016 Central Italy sequence showed a remarkable complexity involving multiple faults. Highly heterogeneous slip distributions were inferred from kinematic finite source inversions. The coverage and quality of the geodetic and seismic data allow resolving high-resolution details of rupture kinematics of the largest event of the sequence, the Mw 6.5 30 October 2016 Norcia earthquake. Composite fault rupture models suggest that two fault planes may have slipped simultaneously. Nevertheless, kinematic modeling cannot assess the mechanic viability of such multiple fault plane models.Using SeisSol, a software package for simulating wave propagation and dynamic rupture based on the arbitrary high-order accurate derivative discontinuous Galerkin method, we therefore try to generate spontaneous dynamic ruptures models compatible with the two fault planes constrained by kinematic inversions. To this end, we adopt a simple slip-weakening friction law with spatially variable dynamic friction and initial strength parameters along multiple faults, compatible with the slip distributions found in the literature. Although we do not to aim to explore the full parameter space, our approach allows testing the feasibility of kinematic models in conjunction with successfully generating spontaneous dynamic rupture scenarios matching seismic and geodetic observations with geological constraints. Such linking enhances and validates the physical implications of kinematic earthquake source inversion.

The 2019 Ridgecrest earthquake sequence culminated in the largest seismic event in California since the 1999 Mw 7.1 Hector Mine earthquake. Here, we combine geodetic and seismic data to study the rupture process of both the 4 July Mw 6.4 foreshock and the 6 July Mw 7.1 mainshock. The results show that the Mw 6.4 foreshock rupture started on a northwest-striking right-lateral fault, and then continued on a southwest-striking fault with mainly left-lateral slip. Although most moment release during the Mw 6.4 foreshock was along the southwest-striking fault, slip on the northwest-striking fault seems to have played a more important role in triggering the Mw 7.1 mainshock that happened ∼34 hr later. Rupture of the Mw 7.1 mainshock was characterized by dominantly right-lateral slip on a series of overall northwest-striking fault strands, including the one that had already been activated during the nucleation of the Mw 6.4 foreshock. The maximum slip of the 2019 Ridgecrest earthquake was ∼5 m, located at a depth range of 3–8 km near the Mw 7.1 epicenter, corresponding to a shallow slip deficit of ∼20%–30%. Both the foreshock and mainshock had a relatively low-rupture velocity of ∼2 km/s, which is possibly related to the geometric complexity and immaturity of the eastern California shear zone faults. The 2019 Ridgecrest earthquake produced significant stress perturbations on nearby fault networks, especially along the Garlock fault segment immediately southwest of the 2019 Ridgecrest rupture, in which the coulomb stress increase was up to ∼0.5 MPa. Despite the good coverage of both geodetic and seismic observations, published coseismic slip models of the 2019 Ridgecrest earthquake sequence show large variations, which highlight the uncertainty of routinely performed earthquake rupture inversions and their interpretation for underlying rupture processes.

We study the October 30th 2016 Norcia earthquake (MW 6.5) to retrieve the rupture history by jointly inverting seismograms and coseismic GPS displacements obtained by dense local networks. The adopted fault geometry consists of a main normal fault striking N155°and dipping 47° belonging to the Mt. Vettore-Mt. Bove fault system (VBFS) and a secondary fault plane striking N210° and dipping 36° to the NW. The coseismic rupture initiated on the VBFS and propagated with similar rupture velocity on both fault planes. Up-dip from the nucleation point, two main slip patches have been imaged on these fault segments, both characterized by similar peak-slip values (~3 m) and rupture times (~3 s). After the breakage of the two main slip patches, coseismic rupture further propagated southeastward along the VBFS, rupturing again the same fault portion that slipped during the August 24th earthquake. The retrieved coseismic slip distribution is consistent with the observed surface breakages and the deformation pattern inferred from InSAR measurements. Our results show that three different fault systems were activated during the October 30th earthquake. The composite rupture model inferred in this study provides evidences that also a deep portion of the NNE-trending section of the Olevano-Antrodoco-Sibillini (OAS) thrust broke co-seismically, implying the kinematic inversion of a thrust ramp. The obtained rupture history indicates that, in this sector of the Apennines, compressional structures inherited from past tectonics can alternatively segment boundaries of NW-trending active normal faults or break co-seismically during moderate-to-large magnitude earthquakes.

The presence of weak phyllosilicates in mature carbonate fault zones has been invoked to explain weak faults. However, the relation between frictional strength, fault stability, mineralogical composition, and fabric of fault gouge, composed of strong and weak minerals, is poorly constrained. We used a biaxial apparatus to systematically shear different mixtures of shale (68% clay, 23% quartz and 4% plagioclase) and calcite, as powdered gouge, at room temperature, under constant normal stresses of 30, 50, 100 MPa and under room-dry and pore fluid-saturated conditions, i.e. CaCO3-equilibrated water. We performed 30 friction experiments during which velocity-stepping and slide-hold-slide tests were employed to assess frictional stability and to measure frictional healing, respectively. Our frictional data indicate that the mineralogical composition of fault gouges significantly affects frictional strength, stability, and healing as well as the presence of CaCO3-equilibrated water. Under room-dry condition, the increasing shale content determines a reduction in frictional strength, from μ = 0.71 to μ = 0.43, a lowering of the healing rates and a transition from velocity-weakening to velocity-strengthening behavior. Under wet condition, with increasing shale content we observe a more significant reduction in frictional strength (μ = 0.65–0.37), a near-zero healing and a velocity strengthening behavior. Microstructural investigations evidence a transition from localized deformation promoted by grain size reduction, in calcite-rich samples, to a more distributed deformation with frictional sliding along clay-enriched shear planes in samples with shale content greater than 50%. For faults cutting across sedimentary sequences composed of carbonates and clay-rich sediments, our results suggest that clay concentration and its ability to form foliated and interconnected networks promotes important heterogeneities in fault strength and slip behavior.

The 2016 central Italy seismic sequence consists so far of a series of moderate‐to‐large earthquakes activating within a few months along a 60‐km‐long and Apenninic‐trending normal‐fault system. Regrettably, the high vulnerability of the local infrastructure and the shallowness of the largest events (depth around 8 km) resulted in 299 casualties and more than 20,000 homeless, with great difficulties in the disaster management. The sequence evolved around its largest event ( M w 6.5, 30 October) that occurred right in the middle of a fault system already activated two months before with a first M w 6.0 mainshock (on 24 August) located to the south near the town of Amatrice. Then, another M w 5.9 mainshock occurred just four days before the largest mainshock (26 October) at the northernmost extent of the sequence, near the town of Visso. We analyze the space–time evolution of the first four months of seismic activity through the relocation of ∼26,000 earthquakes and the kinematic source models of the three mainshocks. All the main events nucleated at the base of a southwest‐dipping normal‐fault system segmented by the presence of crosscutting compressional structures. The presence of these inherited faults, separating diverse geological domains, appears to modulate evolution of the sequence interfering with coseismic slip distribution and fault segments interaction. Several secondary antithetic and synthetic faults are located at a shallow depth (

SUMMARYNumerical simulations of earthquake ground motions are used both to anticipate the effects of hypothetical earthquakes by forward simulation and to infer the behaviour of the real earthquake source ruptures by the inversion of recorded ground motions. In either application it is necessary to assume some Earth structure that is necessarily inaccurate and to use a computational method that is also inaccurate for simulating the wavefield Green's functions. We refer to these two sources of error as ‘propagation inaccuracies’, which might be considered to be epistemic. We show that the variance of the Fourier spectrum of the synthetic earthquake seismograms caused by propagation inaccuracies is related to the spatial covariance on the rupture surface of errors in the computed Green's functions, which we estimate for the case of the 2009 L'Aquila, Italy, earthquake by comparing erroneous computed Green's functions with observed L'Aquila aftershock seismograms (empirical Green's functions). We further show that the variance of the synthetic seismograms caused by the rupture variability (aleatory uncertainty) is related to the spatial covariance on the rupture surface of aleatory variations in the rupture model, and we investigate the effect of correlated variations in Green's function errors and variations in rupture models. Thus, we completely characterize the variability of synthetic earthquake seismograms induced by errors in propagation and variability in the rupture behaviour. We calculate the spectra of the variance of the ground motions of the L'Aquila main shock caused by propagation inaccuracies for two specific broad-band stations, the AQU and the FIAM stations. These variances are distressingly large, being comparable or in some cases exceeding the data amplitudes, suggesting that the best-fitting L'Aquila rupture model significantly overfits the data and might be seriously in error. If these computed variances are typical, the accuracy of many other rupture models for past earthquakes may need to be reconsidered. The results of this work might be useful in seismic hazard estimation because the variability of the computed ground motion, caused both by propagation inaccuracies and variations in the rupture model, can be computed directly, not requiring laborious consideration of multiple Earth structures.

We discuss the debated issue of scale dependence in earthquake source mechanics with the goal of providing supporting evidence to foster the adoption of a coherent interpretative framework. We examine the heterogeneous distribution of source and constitutive parameters during individual ruptures and their scaling with earthquake size. We discuss evidence that slip, slip-weakening distance and breakdown work scale with seismic moment and are interpreted as scale dependent parameters. We integrate our estimates of earthquake stress drop, computed through a pseudo-dynamic approach, with many others available in the literature for both point sources and finite fault models. We obtain a picture of the earthquake stress drop scaling with seismic moment over an exceptional broad range of earthquake sizes (−8

Laboratory and theoretical studies provide insight into the mechanisms that control earthquake nucleation, when fault slip velocity is slow (v∼0.001 cm/s) to dynamic stick slip (v> 1 cm/s). The full range of fault slip modes was achieved by controlling the ratioK=k/kc, wherekis the elastic loading stiffness andkcis the fault zone critical rheologic stiffness. We show thatKcontrols the transition from slow-and-silent slip (K> 0.9) to fast-and-audible (K< 0.7,v= 3 cm/s, slip duration 0.003 s) slip events. Microstructural observations show that with accumulated strain, deformation concentrates in shear zones containing sharp shear planes made of nanoscale grains, which favor the development of frictional instabilities. Once this fabric is established, fault fabric does not change significantly with slip velocity, and fault slip behavior is mainly controlled by the interplay between the rheological properties of the slipping planes and fault zone stiffness.

We present the revised Time Domain Moment Tensor (TDMT) catalogue for earthquakes with M_L larger than 3.6 of the first month of the ongoing Amatrice seismic sequence (August 24th - September 25th). Most of the retrieved focal mechanisms show NNW–SSE striking normal faults in agreement with the main NE-SW extensional deformation of Central Apennines. We also report a preliminary finite fault model analysis performed on the larger aftershock of this period of the sequence (M_w 5.4) and discuss the obtained results in the framework of aftershocks distribution.

Temporal changes in seismic velocity during the earthquake cycle have the potential to illuminate physical processes associated with fault weakening and connections between the range of fault slip behaviours including slow earthquakes, tremor and low-frequency earthquakes1. Laboratory and theoretical studies predict changes in seismic velocity before earthquake failure2; however, tectonic faults fail in a spectrum of modes and little is known about precursors for those modes3. Here we show that precursory changes of wave speed occur in laboratory faults for the complete spectrum of failure modes observed for tectonic faults. We systematically altered the stiffness of the loading system to reproduce the transition from slow to fast stick–slip and monitored ultrasonic wave speed during frictional sliding. We find systematic variations of elastic properties during the seismic cycle for both slow and fast earthquakes indicating similar physical mechanisms during rupture nucleation. Our data show that accelerated fault creep causes reduction of seismic velocity and elastic moduli during the preparatory phase preceding failure, which suggests that real-time monitoring of active faults may be a means to detect earthquake precursors.

SUMMARY In this study, we investigate the rupture history of the 2009 April 6 (Mw 6.1) L’Aquila normal faulting earthquake by using a non-linear inversion of strong motion, GPS and DInSAR data. Both the separate and joint inversion solutions reveal a complex rupture process and a heterogeneous slip distribution. Slip is concentrated in two main asperities: a smaller shallow patch of slip located updip from the hypocentre and a second deeper and larger asperity located southeastwards along-strike direction. The key feature of the source process emerging from our inverted models concerns the rupture history, which is characterized by two distinct stages. The first stage begins with rupture nucleation and with updip propagation at relatively high (∼4.0 km s−1), but still subshear, rupture velocity. The second stage starts nearly 2.0–2.5 s after nucleation and it is characterized by the along-strike rupture propagation. The largest and deeper asperity fails during this stage of the rupture process. The rupture velocity is larger in the updip than in the along-strike direction. The updip and along-strike rupture propagation are separated in time and associated with a Mode II and a Mode III crack, respectively. The comparison between the source models inferred in this study with the Poisson ratio anomalies in the crustal volume containing the fault plane allows the interpretation of the delay in along-strike rupture propagation in terms of a structural control of the rupture history. Our results show that the L’Aquila earthquake featured a very complex rupture, with strong spatial and temporal heterogeneities suggesting a strong frictional and/or structural control of the rupture process.

On 24 August 2016 a magnitude ML 6.0 occurred in the Central Apennines (Italy) between Amatrice and Norcia causing nearly 300 fatalities. The main shock ruptured a NNW-SSE striking, WSW dipping normal fault. We invert waveforms from 26 three-component strong motion accelerometers, filtered between 0.02 and 0.5 Hz, within 45 km from the fault. The inferred slip distribution is heterogeneous and characterized by two shallow slip patches updip and NW from the hypocenter, respectively. The rupture history shows bilateral propagation and a relatively high rupture velocity (3.1 km/s). The imaged rupture history produced evident directivity effects both N-NW and SE of the hypocenter, explaining near-source peak ground motions. Fault dimensions and peak slip values are large for a moderate-magnitude earthquake. The retrieved rupture model fits the recorded ground velocities up to 1 Hz, corroborating the effects of rupture directivity and slip heterogeneity on ground shaking and damage pattern.

On 6 April 2009, a magnitude Mw = 6.1 earthquake struck the Abruzzi region in central Italy. Despite its moderate size, the earthquake caused more than 300 fatalities and partially destroyed the city of L'Aquila and many surrounding villages. The mainshock was preceded by an earthquake swarm that started at the end of 2008. The largest earthquakes of the swarm included an Mw = 4.0, which occurred on 2009/03/30 at 13:38:26 (UTC), and Mw = 3.9 and Mw = 3.5 events that occurred on 2009/04/05 at 20:48 and 22:39 (UTC), respectively. By the end of November 2009, more than 16,000 aftershocks with ML ≥ 0.5 had been recorded by the INGV seismic network (Figure 1). Current advances in data transmission and communication yield high-quality broadband velocity and strong-motion waveforms in near real time. These data are all crucial for rapid determination of earthquake source parameters ( e.g. , fault geometry, focal depth, and seismic moment). For the L'Aquila mainshock, the velocimeter data of the Italian National Seismic Network (INSN, code IV), MedNet (code MN, station PDG), the North-East Italy Broadband Network (code NI, stations ACOM and PALA), and the SudTirol Province (code SI, station KOSI) were available in real time. In the following days, the strong-motion data of the RAN network (Rete Accelerometrica Nazionale) and, in addition, the displacement data recorded by the Istituto Nazionale di Geofisica e Vulcanologia (INGV) GPS network (Anzidei et al. 2009) also become available. In this study we present a retrospective analysis of the rapid source parameters determination procedure developed at INGV (Scognamiglio et al. 2009) as applied to the L'Aquila seismic sequence. Our approach consists of two stages: the near real-time determination of the seismic moment tensor, which is already routinely performed for all ML ≥ 3.5 earthquakes …

SUMMARY In this study, we aim to understand the dependence of the critical slip weakening distance (D c) on the final slip (D tot) during the propagation of a dynamic rupture and the consistency of their inferred correlation. To achieve this goal we have performed a series of numerical tests suitably designed to validate the adopted numerical procedure and to verify the actual capability in measuring D c. We have retrieved two kinematic rupture histories from spontaneous dynamic rupture models governed by a slip weakening law in which a constant D c distribution on the fault plane as well as a constant D c/D tot ratio are assumed, respectively. The slip velocity and the shear traction time histories represent the synthetic ‘real’ target data which we aim to reproduce. We use a 3-D traction-at-split nodes numerical procedure to image the dynamic traction evolution by assuming our modelled slip velocity as a boundary condition on the fault plane. We assume a regularized Yoffe function as source time function in our modelling attempts and we measure the critical slip weakening distance from the inferred traction versus slip curves at each point on the fault. We compare the inferred values with those of the target dynamic models. Our numerical tests show that fitting the slip velocity functions of the target models at each point on the fault plane is not enough to retrieve good traction evolution curves and to obtain reliable measures of D c. We find that the estimation of D c is very sensitive to any small variation of the slip velocity function. An artificial correlation between D c/D tot is obtained when a fixed shape of slip velocity is assumed on the fault (i.e. constant rise time and constant time for positive acceleration) which differs from that of the target model. We point out that the estimation of fracture energy (breakdown work) on the fault is not affected by biases in measuring D c.

We propose a new source-time function, to be used in kinematic modeling of ground-motion time histories, which is consistent with dynamic propagation of earthquake ruptures and makes feasible the dynamic interpretation of kinematic slip models. This function is derived from a source-time function first proposed by Yoffe (1951), which yields a traction evolution showing a slip-weakening behavior. In order to remove its singularity, we apply a convolution with a triangular function and obtain a regularized source-time function called the regularized Yoffe function. We propose a parameterization of this slip-velocity time function through the final slip, its duration, and the duration of the positive slip acceleration ( Tacc ). Using this analytical function, we examined the relation between kinematic parameters, such as peak slip velocity and slip duration, and dynamic parameters, such as slip-weakening distance and breakdown-stress drop. The obtained scaling relations are consistent with those proposed by Ohnaka and Yamashita (1989) from laboratory experiments. This shows that the proposed source-time function suitably represents dynamic rupture propagation with finite slip-weakening distances.

We study the dynamic traction and the slip velocity evolution within the cohesive zone during the propagation of a dynamic rupture using rate- and state-dependent constitutive laws. We solve the elastodynamic equation for a 2-D in-plane crack using a finite difference algorithm. We show that rate and state constitutive laws allow a quantitative description of the dynamic rupture growth. We confirm the findings of previous studies that slip weakening (SW) is a characteristic behavior of rate and state friction. Our simulations show that the state variable evolution controls the slip acceleration and the slip-weakening behavior. These modeling results help in understanding the physical interpretation of the breakdown process and the weakening mechanisms. We compare the time histories of slip velocity, state variable and total dynamic traction to investigate the temporal evolution of slip acceleration and stress drop during the breakdown time. Because the adopted analytical expression for the state variable evolution controls the slip velocity time histories, we test different evolution laws to investigate slip duration and the healing mechanisms. We show that the classic slowness or slip laws do not yield fast restrengthening or self-healing, although they appropriately describe rupture initiation, propagation and the long-term restrengthening during the interseismic period. Self-healing rupture mode, yielding to short slip durations, has been obtained for homogeneous faults by modifying the evolution law introducing a fast restrengthening of dynamic traction immediately after the weakening phase. In this study, we discuss how the direct effect of friction and the friction behavior at high slip rates affect the weakening and healing mechanisms.

On May 20, 2012 (02:03:53 UTC), an Mw 5.86 (Ml 5.9) earthquake struck the Pianura Padana Emiliana region (northern Italy), causing five deaths and damage to several villages and to the towns of Ferrara and Modena. The mainshock was preceded, three hours earlier, by a Mw 3.98 (Ml 4.1) foreshock, which almost co-located with the main event. After the main event, the seismic sequence included six earthquakes with magnitudes >5.0. The biggest aftershock was located about 12 km west of the first mainshock, and was a Mw 5.66 (Ml 5.8) earthquake that occurred on May 29, 2012 (07:00:03 UTC); this can be considered as a second mainshock. After this event, the official death toll of the seismic sequence was 17 people. Moreover, there had been severe damage to the economy of the region and there were 13,000 homeless. [...]

Near-fault seismic recordings for recent earthquakes (Chi Chi earthquake, 1999, and Parkfield earthquake, 2004) show the high spatial heterogeneity of ground motion. This variability is controlled by fault geometry, rupture complexity, and also by wave propagation and site effects. Nowadays, the number of available records in the near-source region is still not enough to infer a robust parameterization of the ground motion and to retrieve multiparametric predictive equations valid at close distances from the fault. The use of a synthetic approach may help to overcome this limitation and to study the strong ground motion variability. In this article we focus on ground-motion dependence on different earthquakes breaking the same fault, as it has been rarely recorded by instruments. We model seismic scenarios from different rupture models of a fault similar to the 1980 Irpinia, Italy, earthquake source ( M w 6.9). A discrete wavenumber/finite element technique is used to compute full-wave displacement and velocity time series in the low-frequency band (up to 2 Hz). We investigate the variability of the ground motion as a function of different source parameters (rupture velocity, slip distribution, nucleation point, and source time function), whose values depend on the state of knowledge of the physical model driving the process. The probability density functions of the simulated ground-motion parameters, such as displacement response spectrum and peak ground velocity, are used to identify particular scenarios that match specific engineering requests.

[1] We image the rupture history of the 2009 L'Aquila (central Italy) earthquake using a nonlinear joint inversion of strong motion and GPS data. This earthquake ruptured a normal fault striking along the Apennines axis and dipping to the SW. The inferred slip distribution is heterogeneous and characterized by a small, shallow slip patch located up-dip from the hypocenter (9.5 km depth) and a large, deeper patch located southeastward. The rupture velocity is larger in the up-dip than in the along-strike direction. This difference can be partially accounted by the crustal structure, which is characterized by a high velocity layer above the hypocenter and a lower velocity below. The latter velocity seems to have affected the along strike propagation since the largest slip patch is located at depths between 9 and 14 km. The imaged slip distribution correlates well with the on-fault aftershock pattern as well as with mapped surface breakages.

This work has been funded by the 2005–2007 DPC-S4 and 2008–2010 DPC-S3 Italian Civil Protection contracts.

We present a conceptual model for the effective critical friction distance for fault zones of finite width. A numerical model with 1D elasticity is used to investigate implications of the model for shear traction evolution during dynamic and quasi-static slip. The model includes elastofrictional interaction of multiple, parallel slip surfaces, which obey rate and state friction laws with either Ruina (slip) or Dieterich (time) state evolution. A range of slip acceleration histories is investigated by imposing perturbations in slip velocity at the fault zone boundary and using radiation damping to solve the equations of motion. The model extends concepts developed for friction of bare surfaces, including the critical friction distance L , to fault zones of finite width containing wear and gouge materials. We distinguish between parameters that apply to a single frictional surface, including L and the dynamic slip weakening distance d o , and those that represent slip for the entire fault zone, which include the effective critical friction distance, D cb , and the effective dynamic slip weakening distance D o . A scaling law for D cb is proposed in terms of L and the fault zone width. Earthquake source parameters depend on net slip across a fault zone and thus scale with D cb , D o , and the slip at yield strength D a . We find that D a decreases with increasing velocity jump size for friction evolution via the Ruina law, whereas it is independent of slip acceleration rate for the Dieterich law. For both laws, D a scales with fault zone width and shear traction exhibits prolonged hardening before reaching a yield strength. The parameters D cb and D o increase roughly linearly with fault zone thickness. This chapter and Chapter 7 in the volume discuss the problem of reconciling laboratory measurements of the critical friction distance with theoretical and field-based estimates of the effective dynamic slip weakening distance.