Deposits of clastic carbonate-dominated (calciclastic) sedimentary slope systems in the rock record have been identified mostly as linearly-consistent carbonate apron deposits, even though most ancient clastic carbonate slope deposits fit the submarine fan systems better. Calciclastic submarine fans are consequently rarely described and are poorly understood. Subsequently, very little is known especially in mud-dominated calciclastic submarine fan systems. Presented in this study are a sedimentological core and petrographic characterisation of samples from eleven boreholes from the Lower Carboniferous of Bowland Basin (Northwest England) that reveals a >250 m thick calciturbidite complex deposited in a calciclastic submarine fan setting. Seven facies are recognised from core and thin section characterisation and are grouped into three carbonate turbidite sequences. They include: 1) Calciturbidites, comprising mostly of highto low-density, wavy-laminated bioclast-rich facies; 2) low-density densite mudstones which are characterised by planar laminated and unlaminated muddominated facies; and 3) Calcidebrites which are muddy or hyper-concentrated debrisflow deposits occurring as poorly-sorted, chaotic, mud-supported floatstones. These

Phd by thesis

Source parameters for historical earthquakes worldwide are compiled to develop a series of empirical relationships among moment magnitude ( M ), surface rupture length, subsurface rupture length, downdip rupture width, rupture area, and maximum and average displacement per event. The resulting data base is a significant update of previous compilations and includes the additional source parameters of seismic moment, moment magnitude, subsurface rupture length, downdip rupture width, and average surface displacement. Each source parameter is classified as reliable or unreliable, based on our evaluation of the accuracy of individual values. Only the reliable source parameters are used in the final analyses. In comparing source parameters, we note the following trends: (1) Generally, the length of rupture at the surface is equal to 75% of the subsurface rupture length; however, the ratio of surface rupture length to subsurface rupture length increases with magnitude; (2) the average surface displacement per event is about one-half the maximum surface displacement per event; and (3) the average subsurface displacement on the fault plane is less than the maximum surface displacement but more than the average surface displacement. Thus, for most earthquakes in this data base, slip on the fault plane at seismogenic depths is manifested by similar displacements at the surface. Log-linear regressions between earthquake magnitude and surface rupture length, subsurface rupture length, and rupture area are especially well correlated, showing standard deviations of 0.25 to 0.35 magnitude units. Most relationships are not statistically different (at a 95% significance level) as a function of the style of faulting: thus, we consider the regressions for all slip types to be appropriate for most applications. Regressions between magnitude and displacement, magnitude and rupture width, and between displacement and rupture length are less well correlated and have larger standard deviation than regressions between magnitude and length or area. The large number of data points in most of these regressions and their statistical stability suggest that they are unlikely to change significantly in response to additional data. Separating the data according to extensional and compressional tectonic environments neither provides statistically different results nor improves the statistical significance of the regressions. Regressions for cases in which earthquake magnitude is either the independent or the dependent parameter can be used to estimate maximum earthquake magnitudes both for surface faults and for subsurface seismic sources such as blind faults, and to estimate the expected surface displacement along a fault for a given size earthquake.

New Empirical Relationships among Magnitude, Rupture Length, Rupture Width, Rupture Area, and Surface Displacement

This essential reference for graduate students and researchers provides a unified treatment of earthquakes and faulting as two aspects of brittle tectonics at different timescales. The intimate connection between the two is manifested in their scaling laws and populations, which evolve from fracture growth and interactions between fractures. The connection between faults and the seismicity generated is governed by the rate and state dependent friction laws - producing distinctive seismic styles of faulting and a gamut of earthquake phenomena including aftershocks, afterslip, earthquake triggering, and slow slip events. The third edition of this classic treatise presents a wealth of new topics and new observations. These include slow earthquake phenomena; friction of phyllosilicates, and at high sliding velocities; fault structures; relative roles of strong and seismogenic versus weak and creeping faults; dynamic triggering of earthquakes; oceanic earthquakes; megathrust earthquakes in subduction zones; deep earthquakes; and new observations of earthquake precursory phenomena.

The Mechanics of Earthquakes and Faulting

Full-waveform inversion FWI is a challenging data-fitting procedure based on full-wavefield modeling to extract quantitative information from seismograms. High-resolution imaging at half the propagated wavelength is expected. Recent advances in high-performance computing and multifold/multicomponent wide-aperture and wide-azimuth acquisitions make 3D acoustic FWI feasible today. Key ingredients of FWI are an efficient forward-modeling engine and a local differential approach, in which the gradient and the Hessian operators are efficiently estimated. Local optimization does not, however, prevent convergence of the misfit function toward local minima because of the limited accuracy of the starting model, the lack of low frequencies, the presence of noise, and the approximate modeling of the wave-physics complexity. Different hierarchical multiscale strategiesaredesignedtomitigatethenonlinearityandill-posedness of FWI by incorporating progressively shorter wavelengths in the parameter space. Synthetic and real-data case studies address reconstructing various parameters, from VP and VS velocities to density, anisotropy, and attenuation. This review attempts to illuminate the state of the art of FWI. Crucial jumps, however, remain necessary to make it as popular as migration techniques. The challenges can be categorized as 1 building accurate starting models with automatic procedures and/or recording low frequencies, 2 defining new minimization criteria to mitigate the sensitivity of FWI to amplitude errors and increasing the robustness of FWI when multiple parameter classes are estimated, and 3 improving computational efficiency by data-compression techniquestomake3DelasticFWIfeasible.

/pdf/an-overview-of-full-waveform-inversion-in-exploration-59sei7fzvh.pdf

An overview of full-waveform inversion in exploration geophysics

I present a finite-difference method for modeling P-SV wave propagation in heterogeneous media This is an extension of the method I previously proposed for modeling SH-wave propagation by using velocity and stress in a discrete grid The two components of the velocity cannot be defined at the same node for a complete staggered grid: the stability condition and the P-wave phase velocity dispersion curve do not depend on the Poisson's ratio, while the S-wave phase velocity dispersion curve behavior is rather insensitive to the Poisson's ratio Therefore, the same code used for elastic media can be used for liquid media, where S-wave velocity goes to zero, and no special treatment is needed for a liquid-solid interface Typical physical phenomena arising with P-SV modeling, such as surface waves, are in agreement with analytical results The weathered-layer and corner-edge models show in seismograms the same converted phases obtained by previous authors This method gives stable results for step discontinuities, as shown for a liquid layer above an elastic half-space The head wave preserves the correct amplitude Finally, the corner-edge model illustrates a more complex geometry for the liquid-solid interface As the Poisson's ratio v increases from 025 to 05, the shear converted phases are removed from seismograms and from the time section of the wave field

/pdf/p-sv-wave-propagation-in-heterogeneous-media-velocity-stress-1srlln7ms7.pdf

P-SV wave propagation in heterogeneous media: Velocity‐stress finite‐difference method

A new finite-difference (FD) method is presented for modeling SH-wave propagation in a generally heterogeneous medium. This method uses both velocity and stress in a discrete grid. Density and shear modulus are similarly discretized, avoiding any spatial smoothing. Therefore, boundaries will be correctly modeled under an implicit formulation. Standard problems (quarter-plane propagation, sedimentary basin propagation) are studied to compare this method with other methods. Finally a more complex example (a salt dome inside a two-layered medium) shows the effect of lateral propagation on seismograms recorded at the surface. A corner wave, always in-phase with the incident wave, and a head wave will appear, which will pose severe problems of interpretation with the usual vertical migration methods.

SH-wave propagation in heterogeneous media; velocity-stress finite-difference method

Probabilistic earthquake location with non-linear, global search methods allows the use of 3D models and produces comprehensive uncertainty and resolution information represented by a probability density function over the unknown hypocentral parameters We describe a probabilistic earthquake location methodology and introduce an efficient Metropolis-Gibbs, non-linear, global sampling algorithm to obtain such locations Using synthetic travel times generated in a 3D model, we examine the locations and uncertainties given by an exhaustive grid-search and the Metropolis-Gibbs sampler using 3D and layered velocity models, and by a iterative, linear method in the layered model We also investigate the relation of average station residuals to known static delays in the travel times, and the quality of the recovery of known focal mechanisms With the 3D model and exact data, the location probability density functions obtained with the Metropolis-Gibbs method are nearly identical to those of the slower but exhaustive grid-search The location PDFs can be large and irregular outside of a station network even for the case of exact data With location in the 3D model and static shifts added to the data, there are systematic biases in the event locations Locations using the layered model show that both linear and global methods give systematic biases in the event locations and that the error volumes do not include the “true” location — absolute event locations and errors are not recovered The iterative, linear location method can fail for locations near sharp contrasts in velocity and outside of a network Metropolis-Gibbs is a practical method to obtain complete, probabilistic locations for large numbers of events and for location in 3D models It is only about 10 times slower than linearized methods but is stable for cases where linearized methods fail The exhaustive grid-search method is about 1000 times slower than linearized methods but is useful for location of smaller number of events and to obtain accurate images of location probability density functions that may be highly-irregular

Probabilistic Earthquake Location in 3D and Layered Models

The nonlinear problem of inversion of seismic waveforms can be set up using least‐squares methods. The inverse problem is then reduced to the problem of minimizing a lp;nonquadratic) function in a space of many (104to106) variables. Using gradient methods leads to iterative algorithms, each iteration implying a forward propagation generated by the actual sources, a backward propagation generated by the data residuals (acting as if they were sources), and a correlation at each point of the space of the two fields thus obtained, which gives the updated model. The quality of the results of any inverse method depends heavily on the realism of the forward modeling. Finite‐difference schemes are a good choice relative to realism because, although they are time‐consuming, they give excellent results. Numerical tests performed with multioffset synthetic data from a two‐dimensional model prove the feasibility of the approach. If only surface‐recorded reflections are used, the high spatial frequency content of the ...

/pdf/two-dimensional-nonlinear-inversion-of-seismic-waveforms-35y9f8due1.pdf

Jean Virieux

Papers

An overview of full-waveform inversion in exploration geophysics

P-SV wave propagation in heterogeneous media: Velocity‐stress finite‐difference method

SH-wave propagation in heterogeneous media; velocity-stress finite-difference method

Probabilistic Earthquake Location in 3D and Layered Models

Two‐dimensional nonlinear inversion of seismic waveforms: Numerical results