Geochemistry of coral from Papua New Guinea as a proxy for ENSO ocean–atmosphere interactions in the Pacific Warm Pool

A Porites sp. coral growing offshore from the Sepik and Ramu Rivers in equatorial northern Papua New Guinea has yielded an accurate 20-year history (1977–1996) of sea surface temperature (SST), river discharge, and wind-induced mixing of the upper water column. Depressions in average SSTs of about 0.5–1.0 °C (indicated by coral Sr/Ca) and markedly diminished freshwater runoff to the coastal ocean (indicated by coral δ¹⁸O, δ¹³C and UV fluorescence) are evident during the El Niño – Southern Oscillation (ENSO) events of 1982–1983, 1987 and 1991-1993. The perturbations recorded by the coral are in good agreement with changes in instrumental SST and river discharge/precipitation records, which are known to be diagnostic of the response of the Pacific Warm Pool ocean–atmosphere system to El Niño. Consideration of coastal ocean dynamics indicates that the establishment of northwest monsoon winds promotes mixing of near-surface waters to greater depths in the first quarter of most years, making the coral record sensitive to changes in the Asian–Australian monsoon cycle. Sudden cooling of SSTs by 1°C following westerly wind episodes, as indicated by the coral Sr/Ca, is consistent with greater mixing in the upper water column at these times. Furthermore, the coral UV fluorescence and oxygen isotope data indicate minimal contribution of river runoff to surface ocean waters at the beginning of most years, during the time of maximum discharge. This abrupt shift in flood-plume behaviour appears to reflect the duration and magnitude of northwest monsoon winds, which tend to disperse flood plume waters to a greater extent in the water column when wind-mixing is enhanced. Our results suggest that a multi-proxy geochemical approach to the production of long coral records should provide comprehensive reconstructions of tropical paleoclimate processes operating on interannual timescales.     

Thallium isotopes as a potential tracer for the origin of cratonic eclogites

Cratonic eclogites are inferred to originate either from subducted ocean crust or mantle melts accreted onto the roots of continents. These models have different implications for the growth of continents, but it is currently difficult to determine the origin of individual eclogite suites.Upper ocean crust altered at low temperatures and marine sediments both display high thallium (Tl) concentrations and strongly fractionated Tl isotope signatures relative to the ambient upper mantle. In this study we carry out the first examination of the suitability of Tl isotopes as a tracer for an ocean-crust origin of cratonic eclogites. We have analysed the Tl isotope composition of clinopyroxene and garnet in six eclogites from the Kaalvallei and Bellsbank kimberlite pipes in South Africa. Minerals were pre-cleaned with an HCl leaching technique and the leachates display variably light Tl isotope ratios. These most likely reflect low-temperature hydrothermal alteration occurring after eruption of the kimberlite that carried the eclogites to the surface.The leached mineral pairs all display identical Tl isotope ratios, strongly suggesting that the source of the analysed Tl is identical for each mineral pair. It is, however, not possible to exclude the possibility that the analysed Tl originates from kimberlitic material that was not removed by the cleaning procedure.Only one of the six samples exhibits a Tl isotope composition different from ambient mantle. Assuming that the Tl isotope signatures indeed represent the eclogite minerals and not any form of contamination, the Tl isotope composition in this sample is consistent with containing a minor component (<3%) of ocean crust altered at low temperatures.Thallium isotopes may become one of the most sensitive indicators for the presence of low-T altered ocean crust because of the stark contrast in Tl concentration and isotopic composition between the mantle and altered ocean crust. In fact, no other chemical or isotopic tracer could have provided an indication that any of the samples studied here had a subduction origin. However, much work is still required before it becomes clear if Tl isotope measurements are a viable means to establish the origin of cratonic eclogites.     

On the extraction of extended structure from Herschel-SPIRE scanning observations in the presence of 1/f noise

We present results from a study of the impact of uncorrelated 1/ f noise on the extraction of spatial structure, on a range of scales, from sky mapping observations made using the Herschel-SPIRE (the spectral and photometric imaging receiver) photometer in the scan-map mode. These studies were carried out using a detailed instrument simulator, and the output reduced using the map-making algorithm to be implemented in the SPIRE data pipeline. The influence of source size scale, telescope-scanning rate and 1/ f noise knee frequency is investigated, and operational bounds to the expected losses are presented, using the case of zero 1/ f (white) noise as a benchmark. Both cross-linked and non-cross-linked observing options are studied. The results presented here represent the best current estimate of the sensitivity of the SPIRE photometer to emission on arbitrary scales. The data presented are general and scalable to any SPIRE observation made using the scanning mode.     

Peak Frequency Dynamics in Solar Microwave Bursts

We analyze the dynamics of the broadband frequency spectrum of 338 microwave bursts observed in the years 2001?–?2002 with the Owens Valley Solar Array. A subset of 38 strong microwave bursts that show a single spectral maximum are studied in detail. Our main goal is to study changes in spectral peak frequency ν _pk with time. We show that, for a majority of these simple bursts, the peak frequency shows a high positive correlation with flux density – it increases on the rise phase in ≈83% of 24 bursts where it could be cleanly measured, and decreases immediately after the peak time in ≈62% of 34 bursts. This behavior is in qualitative agreement with theoretical expectations based on gyrosynchrotron self-absorption. However, for a significant number of events (≈30?–?36%) the peak frequency variation is much smaller than expected from self-absorption, or may be entirely absent. The observed temporal behavior of ν _pk is compared with a simple model of gyrosynchrotron radio emission. We show that the anomalous behavior is well accounted for by the effects of Razin suppression, and further show how an analysis of the temporal evolution of ν _pk can be used to uniquely determine the relative importance of self-absorption and Razin suppression in a given burst. The analysis technique provides a new, quantitative diagnostic for the gyrosynchrotron component of solar microwave bursts. Applying this analysis technique to our sample of bursts, we find that in most of the bursts (60%) the spectral dynamics of ν _pk around the time of peak flux density is caused by self-absorption. On the other hand, for a significant number of events (≈70%), the Razin effect may play the dominant role in defining the spectral peak and dynamics of ν _pk, especially on the early rise phase and late decay phase of the bursts.     

Reconstructing deep‐time histories from integrated thermochronology: An example from southern Baffin Island,Canada

We present a multi‐chronometric approach for reconstructing deep‐time thermal histories using southern Baffin Island as a case study. This continuous thermal history begins with the Palaeoproterozoic Trans‐Hudson Orogeny and is derived from inverse and forward models that integrate thermochronometers spanning some 500°C: new apatite U–Pb ages and K‐feldspar ⁴⁰Ar/³⁹Ar multi‐diffusion domain data, published (U–Th)/He zircon ages and new multi‐kinetic fission‐track results. Integration of data from a wider temperature range reduces ambiguities in thermal‐history modelling and permits us to constrain the timing of geological processes including, extended post‐orogenic cooling, enhanced later Proterozoic cooling, and then episodic burial and exhumation in the Palaeozoic–Mesozoic.     

Reflection moveout inversion in azimuthally anisotropic media: accuracy, limitations and acquisition

Parameter estimation from the elliptical variations in the normal-moveout (NMO) velocity in azimuthally anisotropic media is sensitive to the angular separation between the survey lines in 2D, or equivalently, the source-to-receiver azimuth in 3D, and to the set of azimuths used in the inversion procedure. The accuracy in estimating the orientation of an NMO ellipse, in particular the parameter α, is also sensitive to the magnitude of anisotropy. On the other hand, the accuracy in estimating the semi-axes of the NMO-velocity ellipse is about the same for any magnitude of anisotropy.   To invert for the NMO ellipse parameters at least three NMO-velocity measurements along distinct azimuth directions are needed. In order to maximize the accuracy and stability in parameter estimation, it is best to have the azimuths for the three source-to-receiver directions 60° apart. Having more than three distinct source-to-receiver azimuths (e.g. full azimuthal coverage) provides a useful data redundancy that enhances the quality of the estimates.   In order to maximize quality in the inversion process, it is recommended to design the seismic data acquisition such that it contains small sectors (≤10°) with adequate fold and offset distribution.   Using three NMO-velocity measurements, 60° apart, an azimuthally anisotropic layer overlain by an azimuthally isotropic overburden (as might occur for fractured reservoirs) should have a relative thickness (in time) with respect to the total thickness at least equal to the ratio of the error in the NMO (stacking) velocity to the interval anisotropy of the fractured layer. Coverage along more than three azimuths, however, improves this limitation, which is imposed by Dix differentiation, by at most 50%, depending on the number of observations (NMO velocities) that enter the inversion procedure.     

Harzburgite to lherzolite and back again: metasomatic processes in ultramafic xenoliths from the Wesselton kimberlite, Kimberley, South Africa

Garnets from phlogopite harzburgite xenoliths from the Wesselton kimberlite show zoning from low-Ca harzburgitic cores to rims with lherzolitic Ca-Cr relations. Garnet cores are depleted in Y and HREE, but have sinuous REE patterns enriched in the MREE. Rimwards increase in Ca and decrease in Cr and Mg is accompanied by increases in Zr, Y, Ti and HREE. Secondary replacement rims on some garnets consist of garnet with low Ca and Cr, but high Mg, Ti and HREE. The zoning, and the secondary replacement rims, are attributed to different stages of a metasomatic process that has converted harzburgite to lherzolite, at temperatures near 1000 °C. Modelling of zoning profiles suggests that the process can be divided into three parts: (a) Inwards diffusion of Ca, Zr and Y over periods of 10,000–30,000 years, from a fluid depleted in Ti, Ga and Y; (b) formation of overgrowths high in Ca, Zr, Y and Ti, followed by annealing over periods of several thousand years; (c) formation of secondary reaction rims of low-Ca garnet, on very short timescales prior to eruption. The sinuous REE patterns of the garnet cores are regarded as “primary” features reflecting an ancient metasomatic event superimposed on a depleted protolith. The high Zr/Y, Zr/Ti and Zr/Ca of the fluids corresponding to stage (a) are ascribed to the presence of phlogopite and garnet in the matrix near the fluid source (presumed to be a melt, possibly a kimberlite precursor), leading to the development of concentration fronts in the percolating fluid. The overgrowths of stage (b) appear to coincide with the precipitation of phlogopite in the rock. The low Ca of the fluid responsible for the secondary replacement rims of stage (c) may reflect the late precipitation of clinopyroxene or Ca-carbonate as part of the metasomatic assemblage. These processes have significantly modified the modal, major- and trace-element composition of the mantle volume sampled by the Wesselton kimberlite, within <1 Ma of eruption. Recognition of such effects and their distribution in time and space is essential to understanding of the evolution of the subcontinental lithospheric mantle. Received: 11 February 1998 / Accepted: 24 June 1998     

MoisturEC: A New R Program for Moisture Content Estimation from Electrical Conductivity Data

Noninvasive geophysical estimation of soil moisture has potential to improve understanding of flow in the unsaturated zone for problems involving agricultural management, aquifer recharge, and optimization of landfill design and operations. In principle, several geophysical techniques (e.g., electrical resistivity, electromagnetic induction, and nuclear magnetic resonance) offer insight into soil moisture, but data‐analysis tools are needed to “translate” geophysical results into estimates of soil moisture, consistent with (1) the uncertainty of this translation and (2) direct measurements of moisture. Although geostatistical frameworks exist for this purpose, straightforward and user‐friendly tools are required to fully capitalize on the potential of geophysical information for soil‐moisture estimation. Here, we present MoisturEC, a simple R program with a graphical user interface to convert measurements or images of electrical conductivity (EC) to soil moisture. Input includes EC values, point moisture estimates, and definition of either Archie parameters (based on experimental or literature values) or empirical data of moisture vs. EC. The program produces two‐ and three‐dimensional images of moisture based on available EC and direct measurements of moisture, interpolating between measurement locations using a Tikhonov regularization approach.     

The Dual-Domain Porosity Apparatus: Characterizing Dual Porosity at the Sediment/Water Interface

The characterization of pore-space connectivity in porous media at the sediment/water interface is critical in understanding contaminant transport and reactive biogeochemical processes in zones of groundwater and surface-water exchange. Previous in situ studies of dual-domain (i.e., mobile/less-mobile porosity) systems have been limited to solute tracer injections at scales of meters to hundreds of meters and subsequent numerical model parameterization using fluid concentration histories. Pairing fine-scale (e.g., sub-meter) geoelectrical measurements with fluid tracer data over time alleviates dependence on flowpath-scale experiments, enabling spatially targeted characterization of shallow sediment/water interface media where biogeochemical reactivity is often high. The Dual-Domain Porosity Apparatus is a field-tested device capable of variable rate-controlled downward flow experiments. The Dual-Domain Porosity Apparatus facilitates inference of dual-domain parameters, i.e., mobile/less-mobile exchange rate coefficient and the ratio of less mobile to mobile porosity. The Dual-Domain Porosity Apparatus experimental procedure uses water electrical conductivity as a conservative tracer of differential loading and flushing of pore spaces within the region of measurement. Variable injection rates permit the direct quantification of the flow-dependence of dual-domain parameters, which has been theorized for decades but remains challenging to assess using existing experimental methodologies.     

Statistical techniques for the classification of chromites in diamond exploration samples

A database of 1267 quality-screened major- and trace-element analyses of chromites (s.l.) from kimberlites, lamproites, ultramafic lamprophyres (UMLs) and crustal sources (‘greenstones’, including ophiolites, gabbros, basalts and komatiites) has been subjected to statistical analysis, in order to derive discriminants for use in diamond exploration. The techniques used included nearest-neighbour analysis, CART (classification and regression trees) and MARS (multivariate adaptive regression splines). The results show that both CART and MARS approaches can correctly discriminate kimberlite/lamproite chromites from UML/‘greenstone’ chromites at levels near 90%. Discrimination into the four classes separately is achieved at levels of ca. 80% by MARS and > 70% by CART; much of the misclassification is between the kimberlite and lamproite classes. These results probably represent the maximum probable level of discrimination on chemical criteria, given that ascending magmas may sample both mantle and crustal rocks. The CART approach produces a classification tree that requires no further computation to classify a given grain; the MARS approach requires the use of a simple software package. Tests on known samples illustrate the high level of accuracy of the methods in an exploration context, as well as the useful petrogenetic conclusions that can be drawn from some ‘misclassifications’.