Source author record

Guy Ouillon

Guy Ouillon appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

physics.geo-ph cond-mat.mtrl-sci cond-mat.stat-mech nlin.AO physics.data-an

Catalog footprint

What is connected

7works

5topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2019arXiv

Forecasting the rates of future aftershocks of all generations is essential to develop better earthquake forecast models

Currently, one of the best performing and most popular earthquake forecasting models rely on the working hypothesis that: "locations of past background earthquakes reveal the probable location of future seismicity". As an alternative, we present a class of smoothed seismicity models (SSMs) based on the principles of the Epidemic Type Aftershock Sequence (ETAS) model, which forecast the location, time and magnitude of all future earthquakes using the estimates of the background seismicity rate and the rates of future aftershocks of all generations. Using the Californian earthquake catalog, we formulate six controlled pseudo-prospective experiments with different combination of three target magnitude thresholds: 2.95, 3.95 or 4.95 and two forecasting time horizons: 1 or 5 year. In these experiments, we compare the performance of:(1) ETAS model with spatially homogenous parameters or GETAS (2) ETAS model with spatially variable parameters or SVETAS (3) three declustering based SSMs (4) a simple SSM based on undeclustered data and (5) a model based on strain rate data, in forecasting the location and magnitude of all (undeclustered) target earthquakes during many testing periods. In all conducted experiments, the SVETAS model comes out with consistent superiority compared to all the competing models. Consistently better performance of SVETAS model with respect to declustering based SSMs highlights the importance of forecasting the future aftershocks of all generations for developing better earthquake forecasting models. Among the two ETAS models themselves, accounting for the optimal spatial variation of the parameters leads to strong and statistically significant improvements in forecasting performance.

preprint2013arXiv

Automatic reconstruction of fault networks from seismicity catalogs including location uncertainty

We introduce the Anisotropic Clustering of Location Uncertainty Distributions (ACLUD) method to reconstruct active fault networks on the basis of both earthquake locations and their estimated individual uncertainties. After a massive search through the large solution space of possible reconstructed fault networks, we apply six different validation procedures in order to select the corresponding best fault network. Two of the validation steps (cross-validation and Bayesian Information Criterion (BIC) process the fit residuals, while the four others look for solutions that provide the best agreement with independently observed focal mechanisms. Tests on synthetic catalogs allow us to qualify the performance of the fitting method and of the various validation procedures. The ACLUD method is able to provide solutions that are close to the expected ones, especially for the BIC and focal mechanismbased techniques. The clustering method complemented by the validation step based on focal mechanisms provides good solutions even in the presence of a significant spatial background seismicity rate. Our new fault reconstruction method is then applied to the Landers area in Southern California and compared with previous clustering methods. The results stress the importance of taking into account undersampled sub-fault structures as well as of the spatially inhomogeneous location uncertainties.

preprint2013arXiv

The Barycentric Fixed Mass Method for Multifractal Analysis

We present a novel method to estimate the multifractal spectrum of point distributions. The method incorporates two motivated criteria (barycentric pivot point selection and non-overlapping coverage) in order to reduce edge effects, improve precision and reduce computation time. Implementation of the method on synthetic benchmarks demonstrates the superior performance of the proposed method compared with existing alternatives routinely used in the literature. Finally, we use the method to estimate the multifractal properties of the widely studied growth process of Diffusion Limited Aggregation and compare our results with recent and earlier studies. Our tests support the conclusion of a genuine but weak multifractality of the central core of DLA clusters, with Dq decreasing from 1.75+/-0.01 for q=-10 to 1.65+/-0.01 for q=+10.

preprint2012arXiv

New empirical tests of the multifractal Omori law for Taiwan

We report new tests on the Taiwan earthquake catalog of the prediction by the Multifractal Stress Activation (MSA) model that the p-value of the Omori law for the rate of aftershocks following a mainshock is an increasing function of its magnitude Mm. This effort is motivated by the quest to validate this crucial prediction of the MSA model and to investigate its possible dependence on local tectonic conditions. With careful attention to the long-term as well as short-term time-dependent magnitude completeness of the Taiwan catalog, and with the use of three different declustering techniques, we confirm the universality of the prediction p(Mm) = (0.09 \pm 0.03) \times Mm + (0.47 \pm 0.10), valid for the SCEC Southern California catalog, the Harvard-CMT worldwide catalog, the JMA Japan catalog and the Taiwan catalog. The observed deviations of the two coefficients of the p(Mm) linear dependence on Mm from catalog to catalog are not significant enough to correlate meaningfully with any tectonic features.

preprint2010arXiv

Segmentation of Fault Networks Determined from Spatial Clustering of Earthquakes

We present a new method of data clustering applied to earthquake catalogs, with the goal of reconstructing the seismically active part of fault networks. We first use an original method to separate clustered events from uncorrelated seismicity using the distribution of volumes of tetrahedra defined by closest neighbor events in the original and randomized seismic catalogs. The spatial disorder of the complex geometry of fault networks is then taken into account by defining faults as probabilistic anisotropic kernels, whose structures are motivated by properties of discontinuous tectonic deformation and previous empirical observations of the geometry of faults and of earthquake clusters at many spatial and temporal scales. Combining this a priori knowledge with information theoretical arguments, we propose the Gaussian mixture approach implemented in an Expectation-Maximization (EM) procedure. A cross-validation scheme is then used and allows the determination of the number of kernels that should be used to provide an optimal data clustering of the catalog. This three-steps approach is applied to a high quality relocated catalog of the seismicity following the 1986 Mount Lewis ($M_l=5.7$) event in California and reveals that events cluster along planar patches of about 2 km$^2$, i.e. comparable to the size of the main event. The finite thickness of those clusters (about 290 m) suggests that events do not occur on well-defined euclidean fault core surfaces, but rather that the damage zone surrounding faults may be seismically active at depth. Finally, we propose a connection between our methodology and multi-scale spatial analysis, based on the derivation of spatial fractal dimension of about 1.8 for the set of hypocenters in the Mnt Lewis area, consistent with recent observations on relocated catalogs.

preprint2007arXiv

Automatic Reconstruction of Fault Networks from Seismicity Catalogs: 3D Optimal Anisotropic Dynamic Clustering

We propose a new pattern recognition method that is able to reconstruct the 3D structure of the active part of a fault network using the spatial location of earthquakes. The method is a generalization of the so-called dynamic clustering method, that originally partitions a set of datapoints into clusters, using a global minimization criterion over the spatial inertia of those clusters. The new method improves on it by taking into account the full spatial inertia tensor of each cluster, in order to partition the dataset into fault-like, anisotropic clusters. Given a catalog of seismic events, the output is the optimal set of plane segments that fits the spatial structure of the data. Each plane segment is fully characterized by its location, size and orientation. The main tunable parameter is the accuracy of the earthquake localizations, which fixes the resolution, i.e. the residual variance of the fit. The resolution determines the number of fault segments needed to describe the earthquake catalog, the better the resolution, the finer the structure of the reconstructed fault segments. The algorithm reconstructs successfully the fault segments of synthetic earthquake catalogs. Applied to the real catalog constituted of a subset of the aftershocks sequence of the 28th June 1992 Landers earthquake in Southern California, the reconstructed plane segments fully agree with faults already known on geological maps, or with blind faults that appear quite obvious on longer-term catalogs. Future improvements of the method are discussed, as well as its potential use in the multi-scale study of the inner structure of fault zones.

preprint2003arXiv

Are Aftershocks of Large Californian Earthquakes Diffusing?

We analyze 21 aftershock sequences of California to test for evidence of space-time diffusion. Aftershock diffusion may result from stress diffusion and is also predicted by any mechanism of stress weakening. Here, we test an alternative mechanism to explain aftershock diffusion, based on multiple cascades of triggering. In order to characterize aftershock diffusion, we develop two methods, one based on a suitable time and space windowing, the other using a wavelet transform adapted to the removal of background seismicity. Both methods confirm that diffusion of seismic activity is very weak, much weaker than reported in previous studies. A possible mechanism explaining the weakness of observed diffusion is the effect of geometry, including the localization of aftershocks on a fractal fault network and the impact of extended rupture lengths which control the typical distances of interaction between earthquakes.