Multiple suppression by 2D filtering in the parabolic τ–p domain: a wave-equation-based method1

The filter for wave-equation-based water-layer multiple suppression, developed by the authors in the x-t, the linear τ-p, and the f-k domains, is extended to the parabolic τ-2 domain. The multiple reject areas are determined automatically by comparing the energy on traces of the multiple model (which are generated by a wave-extrapolation method from the original data) and the original input data (multiples + primaries) in τ-p space. The advantage of applying the data-adaptive 2D demultiple filter in the parabolic τ-p domain is that the waves are well separated in this domain. The numerical examples demonstrate the effectiveness of such a dereverberation procedure. Filtering of multiples in the parabolic τ-p domain works on both the far-offset and the near-offset traces, while the filtering of multiples in the f-k domain is effective only for the far-offset traces. Tests on a synthetic common-shot-point (CSP) gather show that the demultiple filter is relatively immune to slight errors in the water velocity and water depth which cause arrival time errors of the multiples in the multiple model traces of less than the time dimension (about one quarter of the wavelet length) of the energy summation window of the filter. The multiples in the predicted multiple model traces do not have to be exact replicas of the multiples in the input data, in both a wavelet-shape and traveltime sense. The demultiple filter also works reasonably well for input data contaminated by up to 25% of random noise. A shallow water CSP seismic gather, acquired on the North West Shelf of Australia, demonstrates the effectiveness of the technique on real data.     

Merging eruption datasets: building an integrated Holocene eruptive record for Mt Taranaki,New Zealand

Acquiring detailed eruption frequency datasets for a volcano system is essential for realistic eruption forecasts. However, accurate datasets are inherently difficult to compile, even if one or more well-dated eruption records are available. A single record typically under-represents the eruption frequency, while combining two or more records may result in an overrepresentation. Although glass compositions have proven to be successful in tephrochronological studies of dominantly rhyolitic tephras; microlitic growth and thin glass shards inhibit their application to andesitic tephras. A method consisting of a combination of two techniques for correlating syn-eruptive deposits is demonstrated on data from the typical andesitic stratovolcano of Mt. Taranaki, New Zealand. Firstly, tentative matches are identified using the radiocarbon age and associated error of each event. Secondly, the compositions of titanomagnetite micro-phenocrysts are used as an independent check, and shown to be a useful correlation tool where age data is available. Using two lake-core records containing tephra layers in an overlapping time-frame, the radiocarbon age-correlation procedure suggested 31 tephra matches. Geochemistry data were available for 15 of these pairs. In three of these cases, the titanomagnetite compositions did not match. Hence, these “paired” tephras were from compositionally distinct magmas and therefore likely represent separate events. An additional three matches were reassigned within the temporal uncertainty limits of the dating procedure, based on better geochemical pairing. The final combined dataset suggests that there have been at least 138 separate ash fall-producing eruptions between 96 and 10 150 years B.P. from Taranaki. Using the combined dataset the mixture of Weibulls renewal model forecasts a probability of 0.52 for an eruption occurring in the next 50 years at this volcano. The present annual eruption probability is estimated at 1.6%. This likelihood is almost double that obtained when relying on a single stratigraphic record.     

First cosmogenic 10Be age constraint on the timing of Younger Dryas glaciation and ice cap thickness,western Scottish Highlands

We use cosmogenic ¹⁰Be surface exposure age techniques at a locality close to Rannoch Moor, western Scottish Highlands, in order to establish the age and chronology of its most recent glaciation. Glacial erratics and an in situ bedrock quartz vein sampled from this site—the summit of Beinn Inverveigh—have yielded zero‐erosion exposure ages of 12.9 ± 1.5 ka to 11.6 ± 1.0 ka, implying complete ice cover of the mountain during the Younger Dryas, or Loch Lomond Stadial. These results fit closely with published ¹⁴C dates that bracket the maximum (lateral) extent of ice cap outlet glaciers, and are the first internally consistent ages to specifically address this period of glaciation in Scotland. Furthermore, the dates imply that previous palaeoglaciological reconstructions for this area may have underestimated both the thickness of the former ice cap and, by implication, its volume. © British Geological Survey/Natural Environment Research Council copyright 2007. Reproduced with the permission of BGS/NERC. Published by John Wiley & Sons, Ltd.     

A gravitational kinetic theory for planetesimals

An analytical theory is developed for the velocity evolution of nonaccreting planetesimal populations, based on the Boltzmann and Fokker-Planck equations. Adapting Shkarofsky's calculation of plasma viscosities, the rate of increase in random velocities due to gravitational encounters between planetesimals of equal mass is found to be one-third to one-half Safronov's result. Comparison with Wetherill's numerical experiments suggests that the Fokker-Planck equation underestimates the effectiveness of encounters and that Safronov's value is approximately correct. For populations of nonuniform sizes, the Fokker-Planck equation indicates an efficient redistribution of energy from the largest bodies to the smaller ones. By conserving angular momentum, the rate of radial spreading of orbits is also derived.     

A neural network approach to classification of sidescan sonarimagery from a midocean ridge area

A neural-network approach to classification of sidescan-sonar imagery is tested on data from three distinct geoacoustic provinces of a midocean-ridge spreading center: axial valley, ridge flank, and sediment pond. The extraction of representative features from the sidescan imagery is analyzed, and the performance of several commonly used texture measures are compared in terms of classification accuracy using a backpropagation neural network. A suite of experiments compares the effectiveness of different feature vectors, the selection of training patterns, the configuration of the neural network, and two widely used statistical methods: Fisher-pairwise classifier and nearest-mean algorithm with Mahalanobis distance measure. The feature vectors compared here comprise spectral estimates, gray-level run length, spatial gray-level dependence matrix, and gray-level differences. The overall accurate classification rates using the best feature set for the three seafloor types are: sediment ponds, 85.9%; ridge flanks, 91.2%; and valleys, 80.1%. While most current approaches are statistical, the significant finding in this study is that high performance for seafloor classification in terms of accuracy and computation can be achieved using a neural network with the proper combination of texture features. These are preliminary results of our program toward the automated segmentation and classification of undersea terrain     

Noise source level density due to surf. I. Monterey Bay, CA

Ambient noise measurements made in Monterey Bay, CA, in 1981 were reduced by estimations of wave-breaking noise and the residual noise was combined with modeled transmission loss (TL) to estimate the spectral source level of surf-generated noise. A Hamilton geoacoustic model of the coastal environment was derived and used in a finite-element parabolic equation propagation-loss model to obtain TL values. Estimates of both the continuous, or local, and discrete components of wave-breaking noise intensity were subtracted from the total measured noise field to determine the contribution due to surf only. Surf breaking on a uniform 12.5-km linear section of beach near Ft. Ord was found to be the dominant source of surf-generated noise. Estimated noise source level densities for heavy surf at Ft. Ord beach varied from 138 dB ref. 1 μPa Hz^-1/2 m at 1 m from the source at 50 Hz to 107 dB at 1 kHz, with a slope of about -5 dB per octave. Although these results must be considered as preliminary, since they are based on a small number of measurements, they may he useful for prediction of ambient noise in other littoral regions