Technological and infrastructure collaborative seismic research in Western Mexico

In February and March 2014, Spanish, Mexican and British scientists and technicians explored the western margin of Mexico, a region with a high occurrence of large earthquakes (> Mw = 7.5) and tsunami generation, on board the British Royal Research Ship James Cook. This successful joint cruise, named TSUJAL, was made possible thanks to a cooperative agreement between NERC and CSIC as part of the Ocean Facilities Exchange Group (OFEG), a major forum of European oceanographic institutions for the exchange of ship time, equipment and personnel. A dense geophysical data set was acquired using for the first time 6 km length seismic streamer facilities from Spain’s Consejo Superior de Investigaciones Cientificas (CSIC), usually operating in the Spanish RV Sarmiento de Gamboa, onboard the British RRS James Cook by solving all mechanical, electrical and electronic problems. The RRS James Cook in turn provides the seismic source and the acoustic, hullmounted echosounder operated by the British Natural Environment Research Council (NERC). Multiscale seismic and echosounder images unravel the subduction geometry, nature of the crust, and evidence faults and mass wasting processes. The data are crucial to estimating fault seismic parameters, and these parameters are critical to carrying out seismic hazard in Mexico, especially when considering largemagnitude earthquakes (Mw 8.0), and to constrain tsunami models.Peer Reviewe

Fast Bayesian Optimal Experimental Design for Seismic Source Inversion

We develop a fast method for optimally designing experiments in the context of statistical seismic source inversion. In particular, we efficiently compute the optimal number and locations of the receivers or seismographs. The seismic source is modeled by a point moment tensor multiplied by a time-dependent function. The parameters include the source location, moment tensor components, and start time and frequency in the time function. The forward problem is modeled by elastodynamic wave equations. We show that the Hessian of the cost functional, which is usually defined as the square of the weighted L2 norm of the difference between the experimental data and the simulated data, is proportional to the measurement time and the number of receivers. Consequently, the posterior distribution of the parameters, in a Bayesian setting, concentrates around the "true" parameters, and we can employ Laplace approximation and speed up the estimation of the expected Kullback-Leibler divergence (expected information gain), the optimality criterion in the experimental design procedure. Since the source parameters span several magnitudes, we use a scaling matrix for efficient control of the condition number of the original Hessian matrix. We use a second-order accurate finite difference method to compute the Hessian matrix and either sparse quadrature or Monte Carlo sampling to carry out numerical integration. We demonstrate the efficiency, accuracy, and applicability of our method on a two-dimensional seismic source inversion problem

Compilation of parameterized seismogenic sources in Iberia for the SHARE European-scale seismic source model.

Abstract: SHARE (Seismic Hazard Harmonization in Europe) is an EC-funded project (FP7) that aims to evaluate European seismic hazards using an integrated, standardized approach. In the context of SHARE, we are compiling a fully-parameterized active fault database for Iberia and the nearby offshore region. The principal goal of this initiative is for fault sources in the Iberian region to be represented in SHARE and incorporated into the source model that will be used to produce seismic hazard maps at the European scale. The SHARE project relies heavily on input from many regional experts throughout the Euro-Mediterranean region. At the SHARE regional meeting for Iberia, the 2010 Working Group on Iberian Seismogenic Sources (WGISS) was established; these researchers are contributing to this large effort by providing their data to the Iberian regional integrators in a standardized format. The development of the SHARE Iberian active fault database is occurring in parallel with IBERFAULT, another ongoing effort to compile a database of active faults in the Iberian region. The SHARE Iberian active fault database synthesizes a wide range of geological and geophysical observations on active seismogenic sources, and incorporates existing compilations (e.g., Cabral, 1995; Silva et al., 2008), original data contributed directly from researchers, data compiled from the literature, parameters estimated using empirical and analytical relationships, and, where necessary, parameters derived using expert judgment. The Iberian seismogenic source model derived for SHARE will be the first regional-scale source model for Iberia that includes fault data and follows an internationally standardized approach (Basili et al., 2008; 2009). This model can be used in both seismic hazard and risk analyses and will be appropriate for use in Iberian- and European-scale assessments

On the simulation of the seismic energy transmission mechanisms

In recent years, considerable attention has been paid to research and development methods able to assess the seismic energy propagation on the territory. The seismic energy propagation is strongly related to the complexity of the source and it is affected by the attenuation and the scattering effects along the path. Thus, the effect of the earthquake is the result of a complex interaction between the signal emitted by the source and the propagation effects. The purpose of this work is to develop a methodology able to reproduce the propagation law of seismic energy, hypothesizing the "transmission" mechanisms that preside over the distribution of seismic effects on the territory, by means of a structural optimization process with a predetermined energy distribution. Briefly, the approach, based on a deterministic physical model, determines an objective correction of the detected distributions of seismic intensity on the soil, forcing the compatibility of the observed data with the physical-mechanical model. It is based on two hypotheses: (1) the earthquake at the epicentre is simulated by means of a system of distortions split into three parameters; (2) the intensity is considered coincident to the density of elastic energy. The optimal distribution of the beams stiffness is achieved, by reducing the difference between the values of intensity distribution computed on the mesh and those observed during four regional events historically reported concerning the Campania region (Italy)

Broadband Records of Earthquakes in Deep Gold Mines and a Comparison with Results from SAFOD, California

For one week during September 2007, we deployed a temporary network of field recorders and accelerometers at four sites within two deep, seismically active mines. The ground-motion data, recorded at 200 samples/sec, are well suited to determining source and ground-motion parameters for the mining-induced earthquakes within and adjacent to our network. Four earthquakes with magnitudes close to 2 were recorded with high signal/noise at all four sites. Analysis of seismic moments and peak velocities, in conjunction with the results of laboratory stick-slip friction experiments, were used to estimate source processes that are key to understanding source physics and to assessing underground seismic hazard. The maximum displacements on the rupture surfaces can be estimated from the parameter Rv, where v is the peak ground velocity at a given recording site, and R is the hypocentral distance. For each earthquake, the maximum slip and seismic moment can be combined with results from laboratory friction experiments to estimate the maximum slip rate within the rupture zone. Analysis of the four M 2 earthquakes recorded during our deployment and one of special interest recorded by the in-mine seismic network in 2004 revealed maximum slips ranging from 4 to 27 mm and maximum slip rates from 1.1 to 6:3 m=sec. Applying the same analyses to an M 2.1 earthquake within a cluster of repeating earthquakes near the San Andreas Fault Observatory at Depth site, California, yielded similar results for maximum slip and slip rate, 14 mm and 4:0 m=sec

On the determination of the earthquake slip distribution via linear programming techniques

The description that one can have of the seismic source is the mani- festation of an imagined model, obviously outlined from Physic Theories and supported by mathematical methods. In that context, the modelling of earthquake rupture consists in finding values of the parameters of the selected physics-mathematical model, through which it becomes possible to reproduce numerically the records of earthquake effects on the Earths surface. We present and test a Linear Programming (LP) inversion in dual form, for reconstructing the kinematics of the rupture of large earthquakes through space-time seismic slip distribution on finite faults planes

Robust Hydraulic Fracture Monitoring (HFM) of Multiple Time Overlapping Events Using a Generalized Discrete Radon Transform

In this work we propose a novel algorithm for multiple-event localization for Hydraulic Fracture Monitoring (HFM) through the exploitation of the sparsity of the observed seismic signal when represented in a basis consisting of space time propagators. We provide explicit construction of these propagators using a forward model for wave propagation which depends non-linearly on the problem parameters - the unknown source location and mechanism of fracture, time and extent of event, and the locations of the receivers. Under fairly general assumptions and an appropriate discretization of these parameters we first build an over-complete dictionary of generalized Radon propagators and assume that the data is well represented as a linear superposition of these propagators. Exploiting this structure we propose sparsity penalized algorithms and workflow for super-resolution extraction of time overlapping multiple seismic events from single well data

Locating earthquakes with surface waves and centroid moment tensor estimation

Traditionally, P wave arrival times have been used to locate regional earthquakes. In contrast, the travel times of surface waves dependent on source excitation and the source parameters and depth must be determined independently. Thus surface wave path delays need to be known before such data can be used for location. These delays can be estimated from previous earthquakes using the cut-and-paste technique, Ambient Seismic Noise tomography, and from 3D models. Taking the Chino Hills event as an example, we show consistency of path corrections for (>10 s) Love and Rayleigh waves to within about 1 s obtained from these methods. We then use these empirically derived delay maps to determine centroid locations of 138 Southern California moderate-sized (3.5 > M_w > 5.7) earthquakes using surface waves alone. It appears that these methods are capable of locating the main zone of rupture within a few (~3) km accuracy relative to Southern California Seismic Network locations with 5 stations that are well distributed in azimuth. We also address the timing accuracy required to resolve non-double-couple source parameters which trades-off with location with less than a km error required for a 10% Compensated Linear Vector Dipole resolution