On the statistical mechanics of distributed seismicity

In order to understand the underlying physics of distributed seismicity better we have considered a 2-D array of slider blocks connected by springs and interacting via static friction with a surface. There is no driving plate in this model. The time evolution of the system is found from numerical simulations in a cellular automata formulation. Energy is conserved and is the single control parameter. The distribution of energies in the springs is found to obey a modified Maxwell-Boltzmann statistics. It is found that the number-size statistics of clusters of unstable sliding blocks is identical to those in percolation clusters in the site-to-site percolation model. There is a well-defined critical point when unstable blocks become connected across the array. It has been previously suggested that distributed seismicity in a seismic zone is the percolation backbone of a 3-D percolation cluster. The fact that low-level seismicity satisfies the Gutenberg-Richter frequency-magnitude relation and is nearly constant in time also suggests that this background seismicity is similar to thermally induced noise.     

The Hatton Bank continental margin—I. Shallow structure from two-ship expanding spread seismic profiles

Summary. The Hatton Bank passive continental margin exhibits thick seaward dipping reflector sequences which consist of basalts extruded during rifting between Greenland and Rockall Plateau. Multichannel seismic reflection profiling across the margin reveals three reflector wedges with a maximum thickness near 7 km, extending from beneath the upper continental slope to the deep ocean basin. We present results of the velocity structure within the dipping reflector sequences at eight locations across the margin, interpreted by synthetic seismogram modelling a set of multichannel expanding spread profiles parallel to the margin. At the top of some reflector sequences, we observe a series of 100 m thick high- and low-velocity zones, which are interpreted as basalt flows alternating with sediments or weathered and rubble layers. At the profile locations, the base of the dipping reflectors correlates with P -wave velocities near 6.5 km s⁻¹. However, elsewhere the reflectors appear to extend significantly deeper than the inferred 6.5 km s⁻¹ velocity contour, indicating that the velocity structure may not be controlled solely by lithological boundaries but also by metamorphic effects. Shear-waves were observed on two lines, permitting the calculation of Poisson's ratio. The decrease in Poisson's ratio from 0.28 to near 0.25 in the upper 5 km of crust may also indicate the effect of metamorphism on seismic properties, or alternatively may be explained by crack closure under load.     

Crustal structure across the Grand Banks–Newfoundland Basin Continental Margin – II. Results from a seismic reflection profile

New multichannel seismic reflection data were collected over a 565 km transect covering the non-volcanic rifted margin of the central eastern Grand Banks and the Newfoundland Basin in the northwestern Atlantic. Three major crustal zones are interpreted from west to east over the seaward 350 km of the profile: (1) continental crust; (2) transitional basement and (3) oceanic crust. Continental crust thins over a wide zone (∼160 km) by forming a large rift basin (Carson Basin) and seaward fault block, together with a series of smaller fault blocks eastwards beneath the Salar and Newfoundland basins. Analysis of selected previous reflection profiles (Lithoprobe 85-4, 85-2 and Conrad NB-1) indicates that prominent landward-dipping reflections observed under the continental slope are a regional phenomenon. They define the landward edge of a deep serpentinized mantle layer, which underlies both extended continental crust and transitional basement. The 80-km-wide transitional basement is defined landwards by a basement high that may consist of serpentinized peridotite and seawards by a pair of basement highs of unknown crustal origin. Flat and unreflective transitional basement most likely is exhumed, serpentinized mantle, although our results do not exclude the possibility of anomalously thinned oceanic crust. A Moho reflection below interpreted oceanic crust is first observed landwards of magnetic anomaly M4, 230 km from the shelf break. Extrapolation of ages from chron M0 to the edge of interpreted oceanic crust suggests that the onset of seafloor spreading was ∼138 Ma (Valanginian) in the south (southern Newfoundland Basin) to ∼125 Ma (Barremian–Aptian boundary) in the north (Flemish Cap), comparable to those proposed for the conjugate margins.     

The relationship between seismic velocity and density in the continental crust — a useful constraint?

Summary . Plots of seismic velocity and density of rock samples show that a range of densities is possible for rocks of each seismic velocity and vice versa. although a single linear relationship is often assumed in crustal gravity calculations. Because of the scatter, whenever rocks of known seismic velocity are converted to density using this relationship, a reduction is made to the resolving power of the resulting gravity calculation. If these rocks reach thicknesses of more than a few kilometres, then the uncertainties become significant when compared with the size of commonly observed gravity anomalies. Examples are considered from the North Sea, Mississippi and Carolina Trough. It is concluded that the use of a seismic velocity measurement as the only indication of rock density does not provide a useful constraint when attempting to reproduce observed gravity variations. An appropriate model for isostatic compensation is probably the most important factor for successful predictions of crustal structure on the basis of gravity data.     

Seismic imaging of the laterally varying D" region beneath the Cocos Plate

We use an axisymmetric, spherical Earth finite difference algorithm to model SH -wave propagation through cross-sections of laterally varying lower mantle models beneath the Cocos Plate derived from recent data analyses. Synthetic seismograms with dominant periods as short as 4 s are computed for several models: (1) a D" reflector 264 km above the core–mantle boundary with laterally varying S -wave velocity increases of 0.9–2.6 per cent, based on localized structures from a 1-D double-array stacking method; (2) an undulating D" reflector with large topography and uniform velocity increase obtained using a 3-D migration method and (3) cross-sections through the 3-D mantle S -wave velocity tomography model TXBW. We apply double-array stacking to assess model predictions of data. Of the models explored, the S -wave tomography model TXBW displays the best overall agreement with data. The undulating reflector produces a double Scd arrival that may be useful in future studies for distinguishing between D" volumetric heterogeneity and D" discontinuity topography. Synthetics for the laterally varying models show waveform variability not observed in 1-D model predictions. It is challenging to predict 3-D structure based on localized 1-D models when lateral structural variations are on the order of a few wavelengths of the energy used, particularly for the grazing geometry of our data. Iterative approaches of computing synthetic seismograms and adjusting model characteristics by considering path integral effects are necessary to accurately model fine-scale D" structure.     

The seismic structure of thinned continental crust in the northern Bay of Biscay

Summary. The stretching and thinning of the continental crust, which occurs during the formation of passive continental margins, may cause important changes in the velocity structure of such crust. Further, crust attenuated to a few kilometres' thickness, can be found underlying 'oceanic' water depths. This paper poses the question of whether thinned continental crust can be distinguished seismically from normal oceanic crust of about the same thickness. A single seismic refraction line shot over thinned continental crust as part of the North Biscay margin transect in 1979 was studied in detail. Tau— p inversion suggested that there are differences between oceanic and continental crust in the lower crustal structure. This was confirmed when synthetic seismograms were calculated. The thinned continental crust (β± 7.0) exhibits a two-gradient structure in the non-sedimentary crust with velocities between 5.9 and 7.4 km s⁻¹; an upper 0.8 s⁻¹ layer overlies a 0.4 s⁻¹ layer. No layer comparable to oceanic layer 3 was detected. The uppermost mantle also contains a low-velocity zone.