Contaminants in sediments: Using robust regression for grain-size normalization

Sediment contaminant concentrations usually show an inverse correlation with grain size. This can cause difficulties in distinguishing real differences in contamination from artifacts caused by variations in sediment texture. To overcome this, regression analysis is frequently used to remove the dependency of concentrations on grain size. However, least squares regression lines can be affected markedly by the presence of a small number of unusual samples in the dataset. These outliers may represent samples which are more severely contaminated or which were derived from areas with different underlying geology. They can be removed semi-manually, but robust regression methods such as least absolute values provide a convenient and objective alternative. The methods are illustrated using an example dataset of metal contaminants in sediments from the Humber Estuary, United Kingdom. Least squares regression on the complete dataset yields a rather poor grain size normalization for several elements. By contrast, least absolute values regression produces results very similar to those obtained by least squares regression after careful manual removal of outliers, but it avoids the need for subjective judgments of which data points to omit from the analysis. The intercepts of several of the fitted regression lines were non-zero, indicating that regression-based normalization is preferable to methods based on ratios.     

Display of acoustic reflection profiles using a dot matrix plotter

Acoustic reflection profiling data display is traditionally done with the aid of a facsimile type of recorder. It is not uncommon to record the unprocessed acoustic data on a tape recorder for subsequent playback through a laboratory computer. This still involves the use of some sort of facsimile recorder for the ultimate display of profiles. This paper presents the results of a study to adapt a high-speed digital dot matrix plotter for the ultimate display in place of the conventional facsimile recorder. Because a minicomputer drives the display directly, a host of signal conditioning procedures are permitted, with the final display being generated in real time. Algorithms are developed to control the marking density, allow adaptive threshold control, bottom tracking, automatic gain control, and de-emphasis of water column boundary reverberation. These techniques are just a few of the many that can be employed since the computer can readily be carried on a large ship in deep water, or a small vessel in a harbour. Shallow water is the difficult case for high energy acoustic sources because the water column boundaries behave much like an excited acoustic cavity. For this reason, a section of seismic profile is shown which was obtained with a 7·5 kHz pinger in only 8 m of water in Narragansett Bay. This research was partiallysupported by the Division of Computer Research of the National Science Foundation.     

The anomalous tides near Broad Sound

Observations of tidal current and height, in conjunction with theoretical mathematical models are used to investigate the propagation of the tide near Broad Sound, a narrowing estuary situated on a wide section of continental shelf toward the southern end of the Great Barrier Reef. The observations indicate that the dense offshore reefs severely inhibit tidal flow, with the result that tides flood toward Broad Sound from the north and from the south, along the main lagoon. There is a local magnification of the semi-diurnal tides within Broad Sound itself.Models of flow across reefs confirm the effectiveness of dense, shallow, and broad reefs in acting as a barrier to the tide. The diffraction of tides through large gaps in the reef is modelled using conformal mapping techniques and with the inclusion of energy leakage, the diffraction model predicts magnification of the semi-diurnal tidal heights by a factor of about 4 and a phase lag of 3 h on the shelf near Broad Sound, these values being consistent with observation. The observed convergence of the tide close to, and within Broad Sound itself is consistent with the proximity of the semi-diurnal tidal period to the natural period for flow in Broad Sound, considered as a narrowing estuary. This results in further amplification, by an additional factor of about 1.5, so that the tides in Broad Sound are increased by a factor of between 5 and 6, altogether, compared with those elsewhere on the east Australian coast.     

The formation and chronology of the PAT 91501 impact-melt L chondrite with vesicle-metal-sulfide assemblages

The L chondrite Patuxent Range (PAT) 91501 is an 8.5-kg unshocked, homogeneous, igneous-textured impact melt that cooled slowly compared to other meteoritic impact melts in a crater floor melt sheet or sub-crater dike [Mittlefehldt D. W. and Lindstrom M. M. (2001) Petrology and geochemistry of Patuxent Range 91501 and Lewis Cliff 88663. Meteoritics Planet. Sci. 36, 439-457]. We conducted mineralogical and tomographic studies of previously unstudied mm- to cm-sized metal-sulfide-vesicle assemblages and chronologic studies of the silicate host. Metal-sulfide clasts constitute about 1 vol.%, comprise zoned taenite, troilite, and pentlandite, and exhibit a consistent orientation between metal and sulfide and of metal-sulfide contacts. Vesicles make up ∼2 vol.% and exhibit a similar orientation of long axes. ³⁹Ar-⁴⁰Ar measurements probably date the time of impact at 4.461 ± 0.008 Gyr B.P. Cosmogenic noble gases and ¹⁰Be and ²⁶Al activities suggest a pre-atmospheric radius of 40-60 cm and a cosmic ray exposure age of 25-29 Myr, similar to ages of a cluster of L chondrites. PAT 91501 dates the oldest known impact on the L chondrite parent body. The dominant vesicle-forming gas was S₂ (∼15-20 ppm), which formed in equilibrium with impact-melted sulfides. The meteorite formed in an impact melt dike beneath a crater, as did other impact melted L chondrites, such as Chico. Cooling and solidification occurred over ∼2 h. During this time, ∼90% of metal and sulfide segregated from the local melt. Remaining metal and sulfide grains oriented themselves in the local gravitational field, a feature nearly unique among meteorites. Many of these metal-sulfide grains adhered to vesicles to form aggregates that may have been close to neutrally buoyant. These aggregates would have been carried upward with the residual melt, inhibiting further buoyancy-driven segregation. Although similar processes operated individually in other chondritic impact melts, their interaction produced the unique assemblage observed in PAT 91501.     

Erosion of the eastern United States observed with 10Be

The problem of identifying areas of accelerated erosion in a dynamic landscape is complicated. The limited history of sediment yield measurements makes this task difficult even if geomorphic evidence is available. Beryllium-10, a cosmogenic isotope produced by cosmic rays interacting with the earth's atmosphere and surface, has chemical and physical properties that make it useful as a tracer for erosion and sediment transport processes. The rarity of the stable isotope, ⁹Be, allows ¹⁰Be to be detected with accelerator mass spectrometry in natural materials at extremely low levels. Backgrounds for rocks and sediments below 10⁵ atom per g are now attainable, a value to be compared with an average deposition rate of 1.3 × 10⁶ atom cm^?2 yr^?1. The affinity of Be for the components of soil and sediment is sufficiently high that it is effectively immobilized on contact, thereby allowing ¹⁰Be to function as a tracer of sediment transport. To a good approximation all the ¹⁰Be transport out of a drainage basin is on the sediment leaving it. The number of ¹⁰Be atoms passing the gauging station can be determined by measuring the concentration of the isotope in the sediment, if the annual sediment load is known. The ratio of the ¹⁰Be carried from the basin by the sediment to that incident upon it, called the erosion index, has been determined for 48 drainage basins within the same physiographic province, which allows them to be reasonably compared, all of which have sediment yield data. Basins located in the Atlantic coastal plain have an average index of 0.3 with the maximum observed being 0.9. Basins located between the fall line and the mountains, a region called the Piedmont, have an average value of 2.2 with individual values ranging from 0.6 to 6.7; this marked difference is thought to result from two centuries of farming on land of moderate gradient. Basins in the highland regions reflect local conditions with low indices for those in grass and timber and high indices associated with destructive land use. The data allow an estimate of the erosion index for the pre-colonial Piedmont, which then allows the pre-colonial sediment yield to be calculated. A number of basins have also been examined world wide with similar conclusions derived. An important deviation from the rule is noted for rivers that erode large regions of loess, such as the Mississippi, Hwang Ho, and Yangtze. Large aeolian deposits were laid down during the ice age in these basins, deposits that brought inherited ¹⁰Be with them and that are easily eroded.     

Exposure history of glass and breccia phases of lunar meteorite EET87521

Abstract— Glass-rich separates were prepared from a sample of the basaltic lunar meteorite EET87521 rich in dark glass. Noble gas isotopic abundances and ²⁶Al and ¹⁰Be activities were measured to find out whether shock effects associated with lunar launch helped to assemble these phases. Similar ¹⁰Be and ²⁶Al activities indicate that all materials in EET87521 had a common exposure history in the last few million years before launch. However, the glass contains much higher concentrations of trapped gases and records a much longer cosmic-ray exposure, 100 Ma–150 Ma, in the lunar regolith than does the bulk sample. The different histories show that the glass existed long before the ejection of EET87521. The trapped ⁴⁰Ar/³⁶Ar ratio of 1.6 ± 0.1 implies that the lunar exposure that produced most of the stable cosmogenic noble gases began 500 Ma ago. Cosmogenic and trapped noble gas components correlate strongly in various temperature-release fractions and phases of EET87521, which is probably because the glass contains most of the gas. The trapped solar ratios, ²⁰Ne/²²Ne = 12.68 ± 0.20 and ³⁶Ar/³⁸Ar = 5.24 ± 0.05 can be understood as resulting from a mixture consisting of ～60% solar wind and 40% solar energetic particles (SEP). All EET87521 phases show a ⁴⁰K-⁴⁰Ar gas retention age of ～3300 Ma, which is in the range of typical lunar mare basalts.     

Estimation and detection issues in matched-field processing

A conceptual framework in which the model-based, space-time acoustic signal processing procedure known as matched field processing (MFP) can be handled in a consistent manner is established. A framework for strong-signal MFP based on standard statistical estimation theory, in which MFP is regarded as essentially an estimation problem in the strong-signal regime, is developed. In the weak-signal case, the necessary requirement of detection dictates that MFP then be considered a joint detection-estimation task. It is demonstrated that, generally, MFP is essentially a space-time processing problem rather than simply an array processing (spatial processing only) procedure. An overview of the processing schemes used to date in MFP is given, showing how these methods relate to the optimal space-time structure. Weak-signal detection and estimation scenarios relevant to MFP are then noted. Present methods for dealing with the inherent instability of MFP algorithms (mismatch) are discussed. The current status of MFP is summarized, and recommendations for future research are made