Wetlands and low-gradient topography are associated with longer hydrologic transit times in Precambrian Shield headwater catchments

The estimation of hydrologic transit times in a catchment provides insights into the integrated effects of water storage, mixing dynamics, and runoff generation processes. There has been limited effort to estimate transit times in southern boreal Precambrian Shield landscapes, which are characteristically heterogeneous with surface cover including till, thin soils, bedrock outcrops, and depressional wetland features that play contrasting hydrologic roles. This study presents approximately 3.5 years of precipitation and streamflow water isotope data and estimates mean transit times (MTTs) and the young water fraction (p_y) across six small catchments in the Muskoka-Haliburton region of south-central Ontario. The main objectives were to define a typical range of MTTs for headwater catchments in this region and to identify landscape variables that best explain differences in MTTs/p_y using airborne light detection and ranging and digital terrain analysis. Of the transit time distributions, the two parallel linear reservoir and gamma distributions best describe the hydrology of these catchments, particularly because of their ability to capture more extreme changes related to events such as snowmelt. The estimated MTTs, regardless of the modelling approach or distribution used, are positively associated with the percent wetland area and negatively with mean slope in the catchments. In this landscape, low-gradient features such as wetlands increase catchment scale water storage when antecedent conditions are dryer and decrease transit times when there is a moisture surplus, which plausibly explains the increases in MTTs and mean annual runoff from catchments with significant coverage of these landscape features.     

Simulations of multimode bolometric interferometers

We describe a procedure for the numerical modelling of astronomical interferometers, with particular relevance to far-infrared and submillimetre wavelengths. The scheme is based on identifying a set of modes that carry power from the sky to the detector. The procedure is extremely general, and can be used to model scalar or vector fields, in any state of coherence and polarization, the only limitation being that the propagation of a coherent field through the system be described by an integral transform, a constraint that is in practise always met.
We present simulations of ideal, multimode two-dimensional interferometers, and show that the modal theory reproduces the correct behaviour of both Michelson and Fizeau interferometers. We calculate simulated visibility data for a multimode bolometric Michelson interferometer, with a synthesized source, and produce a dirty map, recovering the original source with the usual artefacts associated with interferometers.     

Comet nuclei: Morphology and implied processes of surface modification

Surface morphology and related issues for nuclei of three comets: Halley, Borrelly and Wild 2, are considered in the paper. Joint consideration of publications and results of our analysis of the comets’ images led to conclusions, partly new, partly repeating conclusions published by other researchers. It was found that typical for all three nuclei is the presence of rather flat areas: floors of craters and other depressions, mesas and terraces. This implies that flattening surfaces or planation is a process typical for the comet nuclei. Planation seems to work through the sublimation-driven slope collapse and retreat. This requires effective sublimation so this process should work only when a comet is close to the Sun and if on the nucleus there are starting slopes, steep and high enough to support the “long-distance” avalanching of the collapsing material. If the surface had no starting slopes, then instead of planation, the formation of pitted-and-hilly surfaces should occur. An example of this could be the mottled terrain of the Borelly nucleus. Both ways of the sublimational evolution on the nucleus surface should lead to accumulation of cometary regolith. The thickness of the degassed regolith is not known, but it is obvious that in surface depressions, including the flat-floor ones, it should be larger compared with nondepression areas. This may have implications for the in situ study of comets by the Deep Impact and Rosetta missions.Our morphological analysis puts constraints on the applicability of the popular “rubble-pile comet nucleus” hypothesis (Weissman, 1986. Are cometery nuclei primordial rubble piles? Nature 320, 242-244.). We believe that the rubble pile hypothesis can be applicable to the blocky Halley nucleus. The Borelly and Wild 2 nuclei also could be rubble piles. But in these cases the “rubbles” have to be either smaller than 30-50 m (a requirement to keep lineament geometry close to ideal), or larger than 1-2 km (a requirement to form the rather extended smooth, flat surfaces of mesa tops and crater floors). Another option is that the Borelly and Wild 2 nuclei are not rubble piles.In relation to surface morphology we suggest that three end-member types of the comet nuclei may exist: (1) impact cratered “pristine” bodies, (2) non-cratered fragments of catastrophic disruption, and (3) highly Sun-ablated bodies. In this threefold classification, the Wild 2 nucleus is partially ablated primarily cratered body. Borrelly is significantly ablated and could be either primarily cratered or not-cratered fragment. Halley is certainly partially ablated but with the available images it is difficult to say if remnants of impact craters do exist on it.Recently published observations and early results of analysis of the Tempel 1 nucleus images taken by Deep Impact mission are in agreement with our conclusions on the processes responsible for the Halley, Borrelly and Wild 2 nuclei morphologies. In particular, we have now more grounds to suggest that decrease in crater numbers and increase of the role of smooth flat surfaces in the sequence Wild 2?Tempel 1?Borelli reflects a progress in the sublimational degradation of the nucleus surface during comet passages close to the Sun.     

Boundary curves for families of planar orbits

For monoparametric familiesf(x,y)=c of planar orbits, created by a planar potentialV(x,y), we introduce the notion of the family boundary curves (FBC). All members of the familyf(x,y)=c are traced in an allowable region of thexy plane, defined by the corresponding FBC, with total energyE=E(c) varying along the family. Family boundary curves are also found for two-parametric familiesf(x,y,b)=c. The relation of equilibrium points and asymptotic orbits, possibly possessed by the potentialV(x,y), to be FBC is studied.     

An Early Investigation Of The Striated Tail Of Comet Hale-Bopp (C/1995 O1)

The Sekanina-Farrell particle fragmentation model for the striated tails of dust comets is successfully applied to two images of comet Hale-Bopp to study the motions of 12 striae in a time span of March 12–15, 1997. There is evidence for recurring outbursts with a periodicity of 11^h21^m, consistent with results based on analysis of dust jets. The ejecta in all the striae appear to have been released from one source on the nucleus between the end of January and the second half of February 1997, some 60 to 40 days before perihelion. The parent particles were subjected to a radiation pressure acceleration of β_p ≃ 0.55 and their fragmentation lifetimes in 11 of the 12 striae were practically constant and equal to 13–15 days, when normalized to 1 AU from the Sun. Brief analysis of Watanabe et al.'s measurements of striae on their images from March 5–9, 1997 shows even shorter fragmentation lifetimes for the parent particles, mostly about 7–11 days at1 AU. This revised version was published online in July 2006 with corrections to the Cover Date.     

HST/NICMOS Observations of M82

The irregular galaxy M82 is known as the archetypal starburst galaxy. Its proximity (3.5 Mpc) makes this galaxy an ideal laboratory for studying the properties of its starburst. The detailed morphology of the [FeII] 1.644 μm and emission Pa_α (at 1.87 μm) is revealed by the NICMOS images. The peak of the 2.2 μm continuum brightness(evolved population) lies very close to the dynamical centre. Most of the Pa_α emission (which traces the young population) is distributed in a ring of star formation (with a `hole' lacking line emission at the centre of the galaxy). These observations support the scenario in which the starburst in M82 is propagating outwards. It has long been suggested that the [FeII] emission in starburst galaxies can be used as a measure of supernova (SN) activity. M82 shows a large number of radio supernova remnants (SNRs), approximately 50, lying in the plane of the galaxy. The comparison of the positions of the bright compact [FeII] emitting regions with the location of the radio SNRs shows that there is no one-to-one spatial correspondence between the two emissions, suggesting that the radio and [FeII] emissions trace two populations of SNRs with different ages. Young (a few hundred years) SNRs are best traced by their radio emission, whereas the [FeII] stage lasts for at least a few 10⁴ yr. The compact [FeII] sources contribute only some 20 % of the total [FeII] emission observed in M82. However, much of the remaining unresolved [FeII] emission in the plane of the galaxy may arise from SNRs that expanded and merged into a general interstellar medium within a few 10⁴ yr. Presumably, as much as 70% of the total extinction-corrected [FeII]1.644 μm in M82 is associated with SNRs. The extended and diffuse [FeII] component in M82 seems to be related with the superwind above and below the disc of the galaxy. This revised version was published online in July 2006 with corrections to the Cover Date.     

Regolith Surface Reflectance: A New Attempt to Model

The reflectance coefficient of the regolith layer of celestial bodies has been studied in relation to the physical properties of regolith particles (size, refractive index, and packing density) on the basis of an accurate numerical radiative-transfer algorithm for a semi-infinite flat layer. Using the geometric-optics approximation, we have found that a shape mixture of randomly oriented spheroids can successfully model the single-scattering phase function of independent soil grains. In order to take into account the effect of packing density in a regolith layer, the concept of the so-called static structure factor was used. The main effect of increasing packing density is to suppress the forward-scattering peak of the phase function and to increase the albedo of the reflecting surface. We also investigated the influence of fine dust on the reflected light. An addition of small particles not only increases the surface albedo, but also changes the brightness profile and enhances the backscattering. Although the problem of unique solution, which is inherent in the retrieval of the properties of a medium from the measurements of the intensity of light scattered by this media, cannot be removed in the proposed model, the procedure used here, in contrast to widely used approximations, allows us to fit observational data with a set of real characteristics of the regolith. Semiempirical approaches are able to fit the measurements well with a small number of free parameters, but they do not explicitly contain crucial physical characteristics of the regolith such as grain sizes or the refractive index. We compared the numerical solution of the radiative-transfer equation with the Hapke approximation, which is most often used by investigators. The errors introduced by the Hapke model are small only for near-isotropic scattering by isolated particles. However, independent regolith grains are known to scatter light mainly in the forward direction.     

Integral field spectroscopy with the GEMINI multiobject spectrographs

We describe the integral field unit (IFU) which converts the Gemini Multiobject Spectrograph (GMOS) installed on the Gemini-North telescope to an integral field spectrograph,which produces spectra over a contiguous field of view of 7 × 5 arcsec with spatial sampling of 0.2 arcsecover the wavelength range 0.4-1.0 μm.GMOS is converted to this mode by the remote insertion of the IFU into thebeam in place of the masks used for the multiobject mode. A separate fieldof half the area of the main field, but otherwise identical, is alsoprovided to improve background subtraction. The IFU contains 1500lenslet-coupled fibres and was the first facility of any type for integralfield spectroscopy employed on an 8/10 m telescope.We describe the design, construction and testing of the GMOS IFU and present measurements of the throughput both in the laboratory and at the telescope. We compare these with a theoretical prediction made before construction started. All are in good agreement with each other, with the on-telescope throughput exceeding 60% (averaged over wavelength). Finallywe show an example of data obtained during commissioning to illustrate the power of the device.