Time-dependent simulations of steady C-type shocks

Using a time-dependent multifluid, magnetohydrodynamic code, we calculated the structure of steady perpendicular and oblique C-type shocks in dusty plasmas. We included relevant processes to describe mass transfer between the different fluids, radiative cooling by emission lines and grain charging, and studied the effect of single- and multiple-sized grains on the shock structure. Our models are the first of oblique fast-mode molecular shocks in which such a rigorous treatment of the dust grain dynamics has been combined with a self-consistent calculation of the thermal and ionization structures including appropriate microphysics. At low densities, the grains do not play any significant rôle in the shock dynamics. At high densities, the ionization fraction is sufficiently low that dust grains are important charge and current carriers and, thus, determine the shock structure. We find that the magnetic field in the shock front has a significant rotation out of the initial upstream plane. This is most pronounced for single-sized grains and small angles of the shock normal with the magnetic field. Our results are similar to previous studies of steady C-type shocks showing that our method is efficient, rigorous and robust. Unlike the method employed in the previous most detailed treatment of dust in steady oblique fast-mode shocks, ours allow a reliable calculation even when chemical or other conditions deviate from local statistical equilibrium. We are also able to model transient phenomena.     

Generalized polytropes: Stellar structure with a boltzmann factor

In order to study the structure of a chemically homogeneous star in equilibrium, a density profile of the form T ^N exp(–µm(–)/kT) is suggested. As for polytropes, qualitative aspects of the resulting stellar model can be discussed analytically. In particular it is shown that one reobtains forN=3 Eddington's standard model, whereas forN<3 nearly polytropic models result. WhenN>3, the effective polytropic index does vary appreciably over the star. Numerical results indicate that the proposed density profile is quite reasonable in view of the simplicity of the model. From a comparison of the degree of precision of a polytropic approximation with that of the newly proposed model it follows that the new approximation is definitely better than the polytropic one. It is suggested that the model may be useful to study the structure of stellar clouds, clusters and (spherical) galaxies.Now at Department of Applied Mathematics, Queen Mary College, University of London, England.     

The oblateness effect on the solar radiation incident at the top of the atmospheres of the outer planets

Calculations of the daily solar radiation incident at the top of the atmospheres of Jupiter, Saturn, Uranus, and Neptune, with and without the effect of the oblateness, are presented in a series of figures illustrating the seasonal and latitudinal variation of the ratio of both insolations. It is shown that for parts of the summer, the daily insolation of an oblate planet is increased, the zone of enhanced solar radiation being strongly dependent upon the obliquity, whereas the rate of increase is fixed by both the flattening and the obliquity. In winter, the oblateness effect results in a more extensive polar region, the daily solar radiation of an oblate planet always being reduced when compared to a spherical planet. In addition, we also numerically studied the mean daily solar radiation. As previously stated by A.W. Brinkman and J. McGregor (1979, Icarus, 38, 479–482), it is found that in summer the horizon plane is tilted toward the Sun for latitudes less than the subsolar point, but is titled away from the Sun beyond this latitude. It follows that the mean summer daily insolation is increased between the equator and the subsolar point, but decreased poleward of the above-mentioned limit. In winter, however, the horizon plane is always tilted away from the Sun, causing the mean winter daily insolation to be reduced. The partial gain of the mean summertime insolation being much smaller than the loss during winter season evidently yields a mean annual daily insolation which is decreased at all latitudes.     

Corrections to the improved lunar ephemeris

Based primarily upon the formation of new conditional equations using analytical partial derivatives of the moon's mean elements, meridian circle observations of the moon from 1952–67 have been examined to determine corrections to the constants of lunar theory and to the fundamental coordinate system (FK4). With certain exceptions, the new corrections are in agreement with those published earlier by the author. Systematic corrections to FK4 are surprisingly large, although in agreement with some other recent determinations. New corrections to the lunar ephemeris, resulting from the discussion, are also presented.     

Approximate spectral analysis by least-squares fit

An approximate method of spectral analysis called successive spectral analysis based upon the mean-quadratic approximation of an empirical function by generalised trigonometric polynomial with both unknown frequencies and coefficients is developed. A few quotations describing some properties of the method as well as one of the possible methods for numerical solution are given.     

Models for the formation of binary and millisecond radio pulsars

The peculiar combination of a relatively short pulse period and a relatively weak surface dipole magnetic field strength of binary radio pulsars finds a consistent explanation in terms of (i) decay of the surface dipole component of neutron-star magnetic fields on a timescale of (2–5) × 10⁶ yr, in combination with (ii) spin-up of the rotation of the neutron star during a subsequent mass-transfer phase. The four known binary radio pulsars appear to fall into two different categories. Two of them, PSR 0655 + 64 and PSR 1913 + 16, have short orbital periods (<25 h) and high mass functions, indicating companion masses 0.7M⊙ (∼1 (± 0.3) M⊙ and 1.4 M⊙, respectively). The other two, PSR 0820 + 02 and PSR 1953 + 29, have long orbital periods (117d), nearly circular orbits, and low, almost identical mass functions of about 3×10^-3 M⊙, suggesting companion masses of about 0.3M⊙. It is pointed out that these two classes of systems are expected to be formed by the later evolution of binaries consisting of a neutron star and a normal companion star, in which the companion was (considerably) more massive than the neutron star, or less massive than the neutron star, respectively. In the first case the companion of the neutron star in the final system will be a massive white dwarf, in a circular orbit, or a neutron star in an eccentric orbit. In the second case the final companion to the neutron star will be a low-mass (∼ 0.3 M⊙) helium white dwarf in a wide and nearly circular orbit. In systems of the second type the neutron star was most probably formed by the accretion-induced collapse of a white dwarf. This explains in a natural way why PSR 1953 + 29 has a millisecond rotation period and PSR 0820 + 02 has not. Among the binary models proposed for the formation of the 1.5-millisecond pulsar, the only ones that appear to be viable are those in which the companion disappeared by coalescence with the neutron star. In such models the companion may have been a red dwarf of mass 0.03M⊙, a neutron star, or a massive (>0.7M⊙) white dwarf. Only in the last-mentioned case is a position of the pulsar close to the galactic plane a natural consequence. In the first-mentioned case the progenitor system most probably was a cataclysmic-variable binary in which the white dwarf collapsed by accretion.     

Peat bank growth,Holocene palaeoecology and climate history of South Georgia (sub-Antarctica), based on a botanical macrofossil record

Botanical macrofossil analysis of a more than 9000 years old, radiocarbon dated peat sequence of a moss peat bank from South Georgia, shows a clear evolution in the vegetation. Seven ecological phases could be distinguished and they can be interpreted in terms of climate development during the Holocene. Until 2200 years ago, Warnstorfia fontinaliopsis was the dominant moss species pointing to a wet environment. Lower numbers of this species in association with the presence of drier species are assumed to indicate drier periods, such as occurring between ca 6000–5200 and 4400–3400 cal yr BP. The most prominent and definitive vegetation change took place around 2200 cal yr BP. A Polytrichum–Chorisodontium moss peat bank was formed, which is still growing there today. The forcing mechanism for this vegetation change is thought to be a temperature decrease, rather than a precipitation decrease. This conclusion is mainly based on the fact that, today, moss peat banks have their optimal occurrence range in the maritime Antarctic, a region were the mean annual temperature is ca 4 °C lower than on South Georgia. The remarkable change in the moss bank vegetation at 2200 cal yr BP raises the question whether this moment was only a short climatic deterioration, or a definitive change to a cooler and wetter climate after a Holocene climatic optimum period.     

A feather hydrogen isoscape for Mexico

Developing useful biological isoscapes for areas of the world is a priority. This is the case for Mexico that hosts a large percentage of North America's Neotropical migrant birds. Here we investigated the use of House Sparrow (Passer domesticus) feathers to create a spatially explicit feather deuterium isoscape for that country using samples (n = 461) that were collected across Mexico. Considerable and useful spatial hydrogen isotopic structure was observed, suggesting that isotopes may be a potential forensic tool for evaluating origins of Mexican derived fauna and flora. The most positive feather δD values occurred in the northeast and most negative in the south-central part of the country, roughly matching δD patterns observed in groundwater. A weak negative isotopic relationship was found with altitude in both the Pacific and Atlantic drainage systems. The most parsimonious model describing isotopic spatial variation in feathers between 300 and 3000 m a.s.l. included groundwater δD (δD_gw; precipitation proxy), sex, amount of precipitation, and the coefficient of variation in amount of precipitation. Overall, δD_gw was a poor predictor of sparrow δD_f values for all of Mexico. However, this relationship was considerably strengthened when we considered sex separately, removed the Baja peninsula from our sample, and considered the Atlantic and Pacific drainage basins separately. The strongest relationship between δD_gw and δD_f was found for female sparrows in the Atlantic drainage basin (r² = 0.464). We recommend that researchers interested in inferring origins of migratory birds and other animals in Mexico create species specific isotopic basemaps that may be guided by the isotopic patterns we have observed for House Sparrows and groundwater.     

Surface complexation effects on phosphate adsorption to ferric iron oxyhydroxides along pH and salinity gradients in estuaries and coastal aquifers

Non-conservative behavior of dissolved inorganic phosphate (DIP) in estuaries is generally ascribed to desorption from iron and aluminum (hydr)oxides with increasing salinity. Here, we assess this hypothesis by simulating the reversible adsorption of phosphate onto a model oxide (goethite) along physico-chemical gradients representative of surface and subsurface estuaries. The simulations are carried out using a surface complexation model (SCM), which represents the main aqueous speciation and adsorption reactions of DIP, plus the ionic strength-dependent coulombic interactions in solution and at the mineral-solution interface. According to the model calculations, variations in pH and salinity alone are unlikely to explain the often reported production of DIP in surface estuaries. In particular, increased aqueous complexation of phosphate by Mg²⁺ and Ca²⁺ ions with increasing salinity is offset by the formation of ternary Mg-phosphate surface complexes and the drop in electrical potential at the mineral-water interface. However, when taking into account the downstream decrease in the abundance of sorption sites, the model correctly simulates the observed release of DIP in the Scheldt estuary. The sharp increase in pH accompanying the admixing of seawater to fresh groundwater should also cause desorption of phosphate from iron oxyhydroxides during seawater intrusion in coastal aquifers. As for surface estuaries, the model calculations indicate that significant DIP release additionally requires a reduction in the phosphate sorption site density. In anoxic aquifers, this can result from the supply of seawater sulfate and the subsequent reductive dissolution of iron oxyhydroxides coupled to microbial sulfate reduction.