A global search inversion for earthquake kinematic rupture history: Application to the 2000 western Tottori, Japan earthquake

We present a two-stage nonlinear technique to invert strong motions records and geodetic data to retrieve the rupture history of an earthquake on a finite fault. To account for the actual rupture complexity, the fault parameters are spatially variable peak slip velocity, slip direction, rupture time and risetime. The unknown parameters are given at the nodes of the subfaults, whereas the parameters within a subfault are allowed to vary through a bilinear interpolation of the nodal values. The forward modeling is performed with a discrete wave number technique, whose Green’s functions include the complete response of the vertically varying Earth structure. During the first stage, an algorithm based on the heat-bath simulated annealing generates an ensemble of models that efficiently sample the good data-fitting regions of parameter space. In the second stage (appraisal), the algorithm performs a statistical analysis of the model ensemble and computes a weighted mean model and its standard deviation. This technique, rather than simply looking at the best model, extracts the most stable features of the earthquake rupture that are consistent with the data and gives an estimate of the variability of each model parameter. We present some synthetic tests to show the effectiveness of the method and its robustness to uncertainty of the adopted crustal model. Finally, we apply this inverse technique to the well recorded 2000 western Tottori, Japan, earthquake (Mw 6.6); we confirm that the rupture process is characterized by large slip (3-4 m) at very shallow depths but, differently from previous studies, we imaged a new slip patch (2-2.5 m) located deeper, between 14 and 18 km depth

Rupture process of the 2007 Niigata-ken Chuetsu-oki earthquake by non-linear joint inversion of strong motion and GPS data

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

Using geophysical data inversion to constrain earthquake dynamics: a study on dynamically consistent source time functions.

Earthquake kinematic models are often used to retrieve the main parameters of the causative dynamic rupture process. These models are usually obtained through the inversion of seismograms and geodetic data and they can be used as boundary conditions in dynamic modeling to calculate the traction evolution on the fault. Once traction and slip time histories are inferred at each point on the fault plane, it is feasible to estimate the dynamic and breakdown stress drop, the strength excess and the slip weakening distance (Dc). However the measure of these quantities can be biased by the adopted parametrization of kinematic source models. In this work we focus our attention on the importance of adopting source time functions (STFs) compatible with earthquake dynamics to image the kinematic rupture history on a finite fault. First, we compute synthetic waveforms, through a forward modeling, to evaluate the effects of STFs on the ground motion and on the radiated energy. Therefore, adopting different STFs, we perform kinematic inversion of strong motion and GPS data, using a new non linear two-stages search algorithm (Piatanesi et al., 2007) . We have quantitatively verified that the chioce of STFs affects ground motion time histories within the frequency band commonly used in kinematic inversion and that the inferred peak slip velocity and rise time strongly change among the inverted models. These differences has a dramatic impact when kinematic models are used to infer dynamic traction evolution. The shape of the slip weakening curve, the ratio between Dc and the final slip and the dynamic stress drop distribution are remarkably affected by the assumed STFs. We recommend the adoption in kinematic inversions of source time functions that are compatible with earthquake dynamics

Joint Inversion of GPS and Strong Motion Data for Earthquake Rupture Models

Ricostruzione del processo di rottura cosismico su faglia finita attraverso l’inversione congiunta di dati sismologici e geodetici. Implementazione e validazione di una nuova tecnica di inversione non lineare, di tipo global search, per l’inversione congiunta di dati GPS e dati strong motion. Analisi statistica dell’ensemble dei modelli di rottura esplorati dall’algoritmo di inversione. Analisi sulla consistenza dinamica dei modelli cinematici di rottura. Applicazioni: (1) Test Sintetici atti a validare la capacità di risoluzione e robustezza della tecnica sviluppata; (2) Analisi del terremoto di Tottori (2000); (3) Analisi del terremoto di Niigata (2007); (4) Determinazione di scenari di scuotimento in aree di interesse prioritario e strategico (terremoto dell’Irpinia, 1980)

On the use of ground-motion simulations within ShakeMap methodology: application to the 2008 Iwate-Miyagi Nairiku (Japan) and 1980 Irpinia (Italy) earthquakes

ShakeMap package uses empirical Ground Motion Prediction Equations (GMPEs) to estimate the ground motion where recorded data are not available. Recorded and estimated values are then interpolated in order to produce a shaking map associated with the seismic event of interest. The ShakeMap approach better works in regions with dense stations coverage, where the observed ground motions adequately constrain the interpolation. In poorly instrumented regions, the ground motion estimate mainly relies on the GMPE, that account only for average characteristics of source and wave propagation processes. In this study we investigated the improvement of ShakeMap in the near fault area when including synthetic estimates. We focus on the 2008, Mw 7.0, Iwate-Miyagi Nairiku (Japan) earthquake as a case study because recorded by a huge number of stations. As first we calculated the shakemaps to be used as reference maps and then removed several subsets of stations from the original data-set, replacing them with: (i) the estimations of the ground motion obtained by using a specific GMPE valid for that area, using simple source information such as the earthquake magnitude and fault geometry; (ii) the peak values from synthetic time-histories computed with a hybrid deterministic-stochastic method for extended fault, using the rupture fault model obtained from the kinematic source inversion of strong-motion records. We evaluate the deviations from the reference map and the sensitivity to the number of sites where recordings are not available. Our results show that shakemaps are more and more reliable as the coverage of stations is dense and uniformly distributed in the near-source area. Moreover, the synthetics account for propagation and source properties in a more correct way than GMPE, and largely improve the results. The hybrid maps reach good fitting levels especially when synthetics are used to integrate real data and for particular strong-motion parameters and stations’ distribution

Variability of kinematic source parameters and its implication on the choice of the design scenario

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 (Mw 6.9). A discrete wavenumber/finite element technique is used to compute fullwave 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

Variability of kinematic source parameters and its implication on the choice of the design scenario

Near-fault seismic recordings for recent earthquakes (Chi Chi earthquake, 1999; 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 near-source region is still not enough to infer a robust parameterization of the ground motion and to retrieve multi-parametric 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 paper 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 (Mw 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, 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 (SD) and peak ground velocity (PGV), have been used to identify particular scenarios that match specific engineering requests

The Manipulation of Images in Ancient Books: The Case of Sacrobosco's De Sphaera Mundi

Since the time of the invention of printing, ancient texts were graced by rich iconographic collections. These illustrations have a twofold role, decorative and explicative. However, the symbolic meaning of these images could not escape the watchful eye of the Counter-Reformation. Indeed, all science books that were bringing new truths were often censored. However, it should be noted that not only the texts were censored. Sometimes the images were also subject of such repressive actions and, in some cases, they were modified or disfigured. In this paper, we present the case of De sphaera mundi by Johannes de Sacrobosco, conserved in the library of the Astronomical Observatory of Capodimonte. The xylography that accompanies the text shows the mutilated breasts of Urania

Interleukin-12 message in a bottle

IL12 is a very potent cancer immunotherapy agent, but is difficult to harness safely if given systemically. Local gene transfer aims to confine the effects of IL12 to malignant tissues, thus avoiding toxicity. Lipid-nanoparticle mRNA achieves IL12 expression and efficacy in mouse models, opening the way to an ongoing trial

10 Hz GPS seismology for moderate magnitude earthquakes: the case of the Mw 6.3 L’Aquila (Central Italy) event

The 2009 April 6th Mw 6.3 L'Aquila destructive earthquake was successfully recorded by closely spaced 10-Hz and 1-Hz recording GPS receivers and strong motion accelerometers located above or close to the 50° dipping activated fault. We retrieved both static and dynamic displacements from Very High-Rate GPS (VHRGPS) recordings by using Precise Point Positioning kinematic analysis. We compared the GPS positions time series with the closest displacement time series obtained by doubly-integrating strong motion data, first, to assess the GPS capability to detect the first seismic arrivals (P waves) and, secondly, to evaluate the accelerometers capability to detect co-seismic offsets up to ~45 s after the earthquake occurrence. By comparing seismic and VHRGPS frequency contents, we inferred that GPS sampling rates greater than 2.5 Hz (i.e. 5 or 10 Hz) are required in the near-field of moderate magnitude events to provide “alias-free” solutions of coseismic dynamic displacements. Finally, we assessed the consistency of the dynamic VHRGPS results as a constraint on the kinematic rupture history of the mainshock. These results suggested that the high-rate sampling GPS sites in the near field can be as useful as strong motion station for earthquake source studies