Optimising optimal image subtraction

Difference imaging is a technique for obtaining precise relative photometry of variable sources in crowded stellar fields and, as such, constitutes a crucial part of the data reduction pipeline in surveys for microlensing events or transiting extrasolar planets. The Optimal Image Subtraction (OIS) algorithm of Alard & Lupton (1998) permits the accurate differencing of images by determining convolution kernels which, when applied to reference images with particularly good seeing and signal‐to‐noise (S/N), provide excellent matches to the point‐spread functions (PSF) in other images of the time series to be analysed. The convolution kernels are built as linear combinations of a set of basis functions, conventionally bivariate Gaussians modulated by polynomials. The kernel parameters, mainly the widths and maximal degrees of the basis function model, must be supplied by the user. Ideally, the parameters should be matched to the PSF, pixel‐sampling, and S/N of the data set or individual images to be analysed. We have studied the dependence of the reduction outcome as a function of the kernel parameters using our new implementation of OIS within the IDL‐based TRIPP package. From the analysis of noise‐free PSF simulations of both single objects and crowded fields, as well as the test images in the ISIS OIS software package, we derive qualitative and quantitative relations between the kernel parameters and the success of the subtraction as a function of the PSF widths and sampling in reference and data images and compare the results to those of other implementations found in the literature. On the basis of these simulations, we provide recommended parameters for data sets with different S/N and sampling. (© 2007 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Martian basalt (shergottite) Queen Alexandra Range 94201 and lunar basalt 15555: A tale of two pyroxenes

Abstract— We present here ion microprobe analyses of rare earth and other selected trace and minor elements in pyroxenes of shergottite Queen Alexandra Range 94201 and lunar basalt 15555. Pyroxene zonation patterns record the crystallization histories of these two basaltic samples from Mars and the Moon, respectively, and allow a comparison of mafic melt evolution on these two planetary bodies. Elemental abundances and trends in pyroxenes of these two rocks indicate that their minerals formed by continuous, closed system fractional crystallization of their respective parent melts. This further supports the idea that QUE 94201 closely represents the composition of a true Martian basaltic melt (McSween et al., 1996). The main differences in pyroxene elemental zonation patterns in these two objects are attributed to earlier crystallization of whitlockite in QUE 94201 (i.e., before the Fe-rich pyroxenes) than in 15555 (after the Fe-rich pyroxenes). The size of Eu anomalies in pyroxenes of QUE 94201 is intermediate between that in pyroxenes of 15555 and the other shergottites and may imply that fO₂ conditions during crystallization of this Martian basalt were significantly more reducing than for other shergottites, although not quite as reducing as for lunar basalts. Cerium anomalies appear to be less prevalent in pyroxenes of QUE 94201 than other Antarctic shergottites and could be indicative of lesser degree of weathering in the Antarctic.     

Developing a modern pollen–climate calibration data set for Norway

Bjune, A. E., Birks, H. J. B., Peglar, S. M. & Odland, A. 2010: Developing a modern pollen–climate calibration data set for Norway. Boreas, Vol. 39, pp. 674–688. 10.1111/j.1502‐3885.2010.00158.x. ISSN 0300‐9483. Modern pollen–climate data sets consisting of modern pollen assemblages and modern climate data (mean July temperature and mean annual precipitation) have been developed for Norway based on 191 lakes and 321 lakes. The original 191‐lake data set was designed to optimize the distribution of the lakes sampled along the mean July temperature gradient, thereby fulfilling one of the most critical assumptions of weighted‐averaging regression and calibration and its relative, weighted‐averaging partial least‐squares regression. A further 130 surface samples of comparable taphonomy, taxonomic detail and analyst became available as a result of other projects. These 130 samples, all from new lakes, were added to the 191‐lake data set to create the 321‐lake data set. The collection and construction of these data sets are outlined. Numerical analyses involving generalized linear modelling, constrained ordination techniques, weighted‐averaging partial least‐squares regression, and two different cross‐validation procedures are used to asses the effects of increasing the size of the calibration data set from 191 to 321 lakes. The two data sets are used to reconstruct mean July temperature and mean annual precipitation for a Holocene site in northwest Norway and a Lateglacial site in west‐central Norway. Overall, little is to be gained by increasing the modern data set beyond about 200 lakes in terms of modern model performance statistics, but the down‐core reconstructions show less between‐sample variability and are thus potentially more plausible and realistic when based on the 321‐lake data set.     

Sedimentary characteristics of paleo-aeolian dune sands of Salawusu Formation in the Salawusu River Valley

The Salawusu Formation of Milanggouwan section in Salawusu River Valley in-cludes 7 layers of paleo-mobile dune sands, and 4 layers of paleo-fixed and semi-fixed dune sands. Their structures have been observed and their grain size, surface textural features and several main chemical elements have been analyzed. The results showed that: 1) Some of the aeolian structural characteristics of these dune sands are similar to that of the recent dune sands. 2) They are also similar to the recent dune sands in grain size components, and parameters of Mz, σ, Sk and Kg, as well as in several main chemical components. 3) The scattergrams of Mz-σ and SiO2–Al2O3+TOFE and the probability curves of grain size showed that these paleo-dune sands are different from paleosols and fluvio-lacustrine facies, but are consistency with recent dune sands. 4) Quartz sands have well roundness and surface textural features such as dish-shaped pits, crescent-shaped pits, pockmarked pits, upturned cleavage plates, siliceous precipitates and siliceous crevasses, indicating that they had been carried for a long time by the wind. As the 11 layers of paleo-dune sands possess the aeolian characteristics in structure, grain size, surface textural features and chemical elements, the origin of their formation should be attributed to eolation.     

Late Miocene evaporite geochemistry of Lorca and Fortuna basins (Eastern Betics,SE Spain): Evidence of restriction and continentalization

The Lorca and Fortuna basins are two intramontane Neogene basins located in the eastern Betic Cordillera (SE Spain). During the Late Tortonian—Early Messinian, marine and continental evaporites precipitated in these basins as a consequence of increased marine restriction and isolation. Here we show a stratigraphic correlation between the evaporite records of these basins based on geochemical indicators. We use SO₄ isotope compositions and Sr isotopic ratios in gypsum, and halite Br contents to characterize these units and to identify the marine or continental source of the waters feeding the evaporite basins. In addition, we review the available chronological information used to date these evaporites in Lorca (La Serrata Fm), including a thick saline deposit, that we correlate with the First Evaporitic Group in Fortuna (Los Baños Fm). This correlation is also supported by micropalaeontological data, giving a Late Tortonian age for this sequence. The Second Evaporitic Group, (Chicamo Fm), and the Third Evaporitic Group (Rambla Salada Fm) developed only in Fortuna during the Messinian. According to the palaeogeographical scheme presented here, the evaporites of the Lorca and Fortuna basins were formed during the Late Tortonian—Early Messinian, close to the Betic Seaway closure. Sulphate isotope compositions and Sr isotopic ratios of the Ribera Gypsum Mb, at the base of the Rambla Salada Fm (Fortuna basin), match those of the Late Messinian selenite gypsum beds in San Miguel de Salinas, in the near Bajo Segura basin (40 km to the East), and other Messinian Salinity Crisis gypsum deposits in the Mediterranean. According to these geochemical indicators and the uncertainty of the chronology of this unit, the assignment of the Rambla Salada Fm to the MSC cannot be ruled out.     

A high speed near infrared photometer for lunar occultation studies

A recently developed near infrared high speed photometer intended for lunar oc-cultation studies is described. The primary scientific objective is to reach milli arc second levels of angular resolution so that circumstellar structure of the occulted sources can be resolved. Near infrared sky brightness close to the lunar limb is also studied. Angular diameter derived from the observed occultation of IRC +20169 is presented and system performance discussed.     

Normalization and the detection of integrability: The generalized Van Der Waals potential

Deprit and Miller have conjectured that normalization of integrable Hamiltonians may produce normal forms exhibiting degenerate equilibria to very high order. Several examples in the class of coupled elliptic oscillators are known. In order to test the utility of normalization as a detector of integrability we normalize, to high order, a perturbed Keplerian system known to have several integrable limits; the generalized van der Waals Hamiltonian for a hydrogen atom. While the separable limits give rise to high order degeneracy we find a non-separable, integrable limit for which the normal form does not exhibit degeneracy. We conclude that normalization may, in certain cases, indicate integrability but is not guaranteed to uncover all integrable limits.     

X-ray transient sources

The galactic distribution and physical nature of X-ray transient sources is investigated. Two types of transients are considered. The observational data on 41 X-ray transient sources are given, and the average parameters of hard and soft X-ray transients are estimated.     

Pencil-beam surveys for trans-neptunian objects: Limits on distant populations

Two populations of minor bodies in the outer Solar System remain particularly elusive: Scattered Disk Objects and Sedna-like objects. These populations are important dynamical tracers, and understanding the details of their spatial- and size-distributions will enhance our understanding of the formation and on-going evolution of the Solar System. By using newly-derived limits on the maximum heliocentric distances that recent pencil-beam surveys for trans-neptunian objects were sensitive to, we determine new upper limits on the total numbers of distant SDOs and Sedna-like objects. While generally consistent with populations estimated from wide-area surveys, we show that for magnitude-distribution slopes of α ? 0.7-1.0, these pencil-beam surveys provide stronger upper limits than current estimates in literature.