Options
Tinti, Elisa
Loading...
Preferred name
Tinti, Elisa
Email
elisa.tinti@ingv.it
Staff
staff
ORCID
Scopus Author ID
6507918462
64 results
Now showing 1 - 10 of 64
- PublicationOpen AccessVariability in synthetic earthquake ground motions caused by source variability and errors in wave propagation models(2019-06-24)
; ; ; ; ;; ; Numerical 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.283 28 - PublicationOpen AccessRock and fault rheology explain differences between on fault and distributed seismicityAnalysis of seismicity can illuminate active fault zone structures but also deformation within large volumes of the seismogenic zone. For the Mw 6.5 2016-2017 Central Italy seismic sequence, seismicity not only localizes along the major structures hosting the mainshocks (on-fault seismicity), but also occurs within volumes of Triassic Evaporites, TE, composed of alternated anhydrites and dolostones. These volumes of distributed microseismicity show a different frequency-magnitude distribution than on-fault seismicity. We interpret that, during the sequence, shear strain-rate increase, and fluid overpressure promoted widespread ductile deformation within TE that light-up with distributed 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 hydro-fracturing and faulting of dolostones. Our results indicate that ductile crustal deformation can cause distributed microseismicity, which obeys to different scaling laws than on-fault seismicity occurring on structures characterized by elasto-frictional stick-slip behaviour.
24 6 - PublicationRestrictedRupture process of the 2007 Niigata-ken Chuetsu-oki earthquake by non-linear joint inversion of strong motion and GPS data(2008-08-23)
; ; ; ; ;Cirella, A.; Istituto Nazionale di Geofisica e Vulcanologia, Sezione Roma1, Roma, Italia ;Piatanesi, A.; Istituto Nazionale di Geofisica e Vulcanologia, Sezione Roma1, Roma, Italia ;Tinti, E.; Istituto Nazionale di Geofisica e Vulcanologia, Sezione Roma1, Roma, Italia ;Cocco, M.; Istituto Nazionale di Geofisica e Vulcanologia, Sezione Roma1, Roma, Italia; ; ; We image the rupture history of the 2007 Niigata-ken Chuestu-oki (Japan) earthquake by a nonlinear joint inversion of strong motion and GPS data, retrieving peak slip velocity, rupture time, rise time and slip direction. The inferred rupture model contains two asperities; a small patch near the nucleation and a larger one located 10÷15 km to the south-west. The maximum slip ranges between 2.0 and 2.5 m and the total seismic moment is 1.6×1019 Nm. The inferred rupture history is characterized by rupture acceleration and directivity effects, which are stable features of the inverted models. These features as well as the source-to-receiver geometry are discussed to interpret the high peak ground motions observed (PGA is 1200 gals) at the Kashiwazaki-Kariwa nuclear power plant (KKNPP), situated on the hanging-wall of the causative fault. Despite the evident source effects, predicted PGV underestimates the observed values at KKNPP by nearly a factor of 10.277 352 - PublicationRestrictedThe role of shale content and pore-water saturation on frictional properties of simulated carbonate faults(2021)
; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; 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.38 2 - PublicationRestrictedEarthquake fracture energy inferred from kinematic rupture models on extended faults(2005)
; ; ; ;Tinti, E.; Istituto Nazionale di Geofisica e Vulcanologia, Sezione Roma1, Roma, Italia ;Spudich, P.; U.S. Geological Survey ;Cocco, M.; Istituto Nazionale di Geofisica e Vulcanologia, Sezione Roma1, Roma, Italia; ; We estimate fracture energy on extended faults for several recent earthquakes by retrieving dynamic traction evolution at each point on the fault plane from slip history imaged by inverting ground motion waveforms. We define the breakdown work (Wb) as the excess of work over some minimum traction level achieved during slip. Wb is equivalent to "seismological" fracture energy (G) in previous investigations. Our numerical approach uses slip velocity as a boundary condition on the fault. We employ a three-dimensional finite difference algorithm to compute the dynamic traction evolution in the time domain during the earthquake rupture. We estimate Wb by calculating the scalar product between dynamic traction and slip velocity vectors. This approach does not require specifying a constitutive law and assuming dynamic traction to be collinear with slip velocity. If these vectors are not collinear, the inferred breakdown work depends on the initial traction level. We show that breakdown work depends on the square of slip. The spatial distribution of breakdown work in a single earthquake is strongly correlated with the slip distribution. Breakdown work density and its integral over the fault, breakdown energy, scale with seismic moment according to a power law (with exponent 0.59 and 1.18, respectively). Our estimates of breakdown work range between 4e+5 and 2e+7 J/m2 for earthquakes having moment magnitudes between 5.6 and 7.2. We also compare our inferred values with geologic surface energies. This comparison might suggest that breakdown work for large earthquakes goes primarily into heat production.191 28 - PublicationOpen AccessThe Role of Fault Rock Fabric in the Dynamics of Laboratory Faults(2022)
; ; ; ; ; ; ;; ; Fault stability is inherently linked to the frictional and healing properties of fault rocks and associated fabrics. Their complex interaction controls how the stored elastic energy is dissipated, that is, through creep or seismic motion. In this work, we focus on the relevance of fault fabrics in controlling the reactivation and slip behavior of dolomite-anhydrite analog faults. We designed a set of laboratory experiments where we first develop fault rocks characterized by different grain size reduction and localization at normal stresses of σN = 15, 35, 60, and 100 MPa and second, we reload and reactivate these fault rocks at the frictional stability transition, achieved at σN = 35 MPa by reducing the machine stiffness. If normal stress is lowered this way, reactivation occurs with relatively large stress drops and large peak-slip velocities. Subsequent unstable behavior produces slow stick-slip events with low stress drop and with either asymmetric or Gaussian slip velocity function depending on the inherited fault fabric. If normal stress is raised, deformation is accommodated within angular cataclasites promoting stable slip. The integration of microstructural data (showing brittle reworking of preexisting textures) with mechanical data (documenting restrengthening and dilation upon reactivation) suggests that frictional and chemically assisted healing, which is common in natural faults during the interseismic phase, can be a relevant process in developing large instabilities. We also conclude that fault rock heterogeneity (fault fabric) modulates the slip velocity function and thus the dynamics of repeating stick-slip cycles.30 60 - PublicationOpen Access
62 66 - PublicationOpen AccessOn the evolution of elastic properties during laboratory stick-slip experiments spanning the transition from slow slip to dynamic rupture(2016-12)
; ; ; ; ; ; ; ; ; ; ; The physical mechanisms governing slow earthquakes remain unknown, as does the relationship between slow and regular earthquakes. To investigate the mechanism(s) of slow earthquakes and related quasi‐dynamic modes of fault slip we performed laboratory experiments on simulated fault gouge in the double direct shear configuration. We reproduced the full spectrum of slip behavior, from slow to fast stick slip, by altering the elastic stiffness of the loading apparatus (k) to match the critical rheologic stiffness of fault gouge (kc). Our experiments show an evolution from stable sliding, when k > kc, to quasi‐dynamic transients when k ~ kc, to dynamic instabilities when k < kc. To evaluate the microphysical processes of fault weakening we monitored variations of elastic properties. We find systematic changes in P wave velocity (Vp) for laboratory seismic cycles. During the coseismic stress drop, seismic velocity drops abruptly, consistent with observations on natural faults. In the preparatory phase preceding failure, we find that accelerated fault creep causes a Vp reduction for the complete spectrum of slip behaviors. Our results suggest that the mechanics of slow and fast ruptures share key features and that they can occur on same faults, depending on frictional properties. In agreement with seismic surveys on tectonic faults our data show that their state of stress can be monitored by Vp changes during the seismic cycle. The observed reduction in Vp during the earthquake preparatory phase suggests that if similar mechanisms are confirmed in nature high‐resolution monitoring of fault zone properties may be a promising avenue for reliable detection of earthquake precursors.400 132 - PublicationOpen AccessAsperity size and neighboring segments can change the frictional response and fault slip behavior: insights from laboratory experiments and numerical simulations(2024)
; ; ; ; ; ; ; ; ; ; ; ; ; Accurate assessment of the rate and state friction parameters of rocks is essential for producing realistic earthquake rupture scenarios and, in turn, for seismic hazard analysis. Those parameters can be directly measured on samples, or indirectly based on inversion of coseismic or postseismic slip evolution. However, both direct and indirect approaches require assumptions that might bias the results. Aiming to reduce the potential sources of bias, we take advantage of a downscaled analog model reproducing megathrust earthquakes. We couple the simulated annealing algorithm with quasi-dynamic numerical models to retrieve rate and state parameters reproducing the recurrence time, rupture duration and slip of the analog model, in the ensemble. Then, we focus on how the asperity size and the neighboring segments’ properties control the seismic cycle characteristics and the corresponding variability of rate and state parameters. We identify a tradeoff between (a-b) of the asperity and (a-b) of neighboring creeping segments, with multiple parameter combinations that allow mimicking the analog model behavior. Tuning of rate and state parameters is required to fit laboratory experiments with different asperity lengths. Poorly constrained frictional properties of neighboring segments are responsible for uncertainties of (a-b) of the asperity in the order of per mille. Roughly one order of magnitude larger uncertainties derive from asperity size. Those results provide a glimpse of the variability that rate and state friction estimates might have when used as a constraint to model fault slip behavior in nature.61 22 - PublicationOpen AccessOn the 2016-2017 Central Italy Seismic Sequence: Uncertainty for Source Parameters of Mainshocks and Revised Catalogue of Moment Tensors(2019-12)
; ; ; ; ; ; ; ; ; The Central Italy seismic sequence began on August 24th, 2016, and was marked by three mainshocks in two months culminated with the Mw 6.5, October 30th, 2016, event. Location, depth and prevalent normal faulting mechanisms indicate that the sequence originated in the shallow crust of the Apennine chain where the current extensional regime overprints contractional structures. Structural complexity plays a major role in fault segmentation and interaction in this region, with important consequences on seismic behavior and mechanics of earthquake faulting. This complexity is evidenced by the co-existence of fault planes with different focal mechanisms in the same area. Here we analyze the robustness of moment tensor solutions for the three mainshocks of the 2016-2017 Central Italy sequence. In particular, we study the effect of number and distribution of the inverted stations and employed wave speed model (1D and 3D) with the goal of providing more reliable estimates of the source parameters (strike, dip, rake and Mw) and corresponding uncertainties. The latter are estimated by performing a bootstrap analysis on hundreds of solutions computed by varying the distribution of stations for 1D and 3D velocity models. Moreover, we report on reviewed source geometries of the Central Italy sequence as retrieved by moment tensor analysis by integrating the actual TDMT revised catalogue (http://terremoti.ingv.it/tdmt) for M4+, with new updated solutions based on a new Italian 3D wave speed model. The realization of a complete moment tensor catalogue, in addition to the estimate of uncertainties associated to the computed focal planes and Mw for the three mainshocks, can contribute to explain the complexity of the seismogenic processes active in the Central Apennines and help in understanding the main features of this seismic sequence.21 3