Origin of short period comets

Reasons for interest in the origin of short-period comets and the difficulties of computing their long-term dynamcal evolution are reviewed. The relative advantages of a source region in an extended inner core of the Oort cloud or a compact comet belt just beyond the planetary system are finely balanced, and it is premature to consider the problem solved. A complication is that some comets belonging to the Jupiter family may be part of a time-dependent system, possibly the remains of a giant comet such as Chiron which could have been part of the system 104 yr ago. The origin of short-period comets plays a pivotal role in many areas of solar system science: planet formation, the source of water (possibly life) on the terrestrial planets, the cratering record on the terrestrial planets and satellites of the outer planets, and the environmental impact posed by massive bodies and their decay products in the Earth's near-space environment.     

Thermal balance of the mid-latitude F2-region at night

The thermal balance of the plasma in the night-time mid-latitude F2-region is examined using solutions of the steady-state O⁺ and electron heat balance equations. The required concentrations and field-aligned velocities are obtained from a simultaneous solution of the time-dependent O⁺ continuity and momentum equations.The results demonstrate the systematic trend for the O⁺ temperature to be 10–20 K greater than the electron temperature during the night at around 300 km, as observed at St. Santin by Bauer and Mazaudier. It is shown that frictional heating between the O⁺ and neutral gases is the cause of the O⁺ temperature being greater than the electron temperature; the greater the importance of frictional heating in the thermal balance the greater is the difference in the O⁺ and electron temperatures. A study is made of the roles played in the thermal balance of the plasma by the thermal conductivity of the O⁺ and electron gases; collisional heat transfer between O⁺ electrons and neutrals; frictional heating between the O⁺ and neutral gases; and advection and convection due to field-aligned O⁺ and electron motions. The results of the study show that, at around 300 km, electron cooling by excitation of the fine structure of the ground state of atomic oxygen plays a major role in the thermal balance of the electrons and, since the temperature of the ions is little affected by this electron cooling process, in determining the difference between the ion and electron temperatures.     

Interhemispheric flow of thermal plasma in a closed magnetic flux tube at mid-latitudes under sunspot minimum conditions

The coupled H⁺ and O⁺ time-dependent continuity and momentum equations are solved within a region of the L = 3 magnetic flux tube lying between (and including) the F2-layers of conjugate hemispheres. The method of solution is an extended and modified version of the Murphy et al. (1976) method. The model is used to study the coupling between the F2-layers of conjugate hemispheres during magnetically quiet periods.The results of the calculations strongly indicate that the protonosphere acts as a reservoir, with variable H⁺ content, which prevents direct coupling between the F2-layers of conjugate hemispheres. However there is generally a significant interhemispheric flow of plasma. This flow is caused by conditions in the summer and winter topside ionospheres and it maintains continuity in the plasma concentration within the protonosphere. There are times when the direction of flow is from the winter hemisphere to the summer hemisphere. It is suggested that maintenance of the winter F2-layer at night is not assisted directly by the F2-layer of the conjugate summer hemisphere.It is shown that during the first few days of protonosphere replenishment after a magnetic storm there is an upflow of H⁺ in the topside ionosphere at all times in the summer hemisphere. There is also an upflow of H⁺ during the daytime in both hemispheres. A comparison with the results obtained when the interhemispheric H⁺ flux is held permanently at zero shows that both F2-layers are little affected by the interhemispheric H⁺ flux. Nevertheless both F2-layers are affected by the H⁺ tube content of the protonosphere. When the H⁺ flux at 1000 km in one hemisphere is much greater than the H⁺ flux at 1000 km in the conjugate hemisphere, there is a corresponding signature in the interhemispheric H⁺ flux.The results suggest that there is insufficient time between magnetic storms for complete replenishment of the protonosphere to occur.     

Comet Ikeya-Zhang rises for the NAM

A possible naked-eye comet that may have been important in early cometary theory is announced by D J Asher , M E Bailey , A Christou , J McFarland , M F Muir and P P Rafferty .
Early indications sugest that Comet C/2002 (Ikeya-Zhang), discovered on 1 February 2002, may brighten to naked-eye visibility in late March 2002. It has also been suggested that it may be identical to one of the brighter comets of the 16th or 17th centuries, namely C/1532 R1 or C/1661 C1. The first of these, observed for more than 100 days towards the end of 1532, played an important role in the development of cometary theory. The second, although identified by Halley as having an orbit similar to that of the comet C/1532 R1, was not seen on its predicted return in 1788/1789 and so was presumably unrelated. Here we present long-term orbital integrations of C/2002 C1 which suggest that it orginated from the Oort cloud, and will be ejected again, within ˜0.3 Myr. There is a chance of 10–20% that it will end its life by falling into the Sun during a Halley-type phase of cometary evolution. The discovery of Ikeya-Zhang so closeto perigee by two amateur astonomers highlights the need for surveys covering both hemispheres to discover long-period and intermediate-period comets on Earth-crossing orbits.     

Formation and pathways of North Atlantic Deep Water in a coupled ice–ocean model of the Arctic–North Atlantic Oceans

We investigate the formation process and pathways of deep water masses in a coupled ice–ocean model of the Arctic and North Atlantic Oceans. The intent is to determine the relative roles of these water masses from the different source regions (Arctic Ocean, Nordic Seas, and Subpolar Atlantic) in the meridional overturning circulation. The model exhibits significant decadal variability in the deep western boundary current and the overturning circulation. We use detailed diagnostics to understand the process of water mass formation in the model and the resulting effects on the North Atlantic overturning circulation. Particular emphasis is given to the multiple sources of North Atlantic Deep Water, the dominant deep water masses of the world ocean. The correct balance of Labrador Sea, Greenland Sea and Norwegian Sea sources is difficult to achieve in climate models, owing to small-scale sinking and convection processes. The global overturning circulation is described as a function of potential temperature and salinity, which more clearly signifies dynamical processes and clarifies resolution problems inherent to the high latitude oceans. We find that fluxes of deep water masses through various passages in the model are higher than observed estimates. Despite the excessive volume flux, the Nordic Seas overflow waters are diluted by strong mixing and enter the Labrador Sea at a lighter density. Through strong subpolar convection, these waters along with other North Atlantic water masses are converted into the densest waters [similar density to Antarctic Bottom Water (AABW)] in the North Atlantic. We describe the diminished role of salinity in the Labrador Sea, where a shortage of buoyant surface water (or excess of high salinity water) leads to overly strong convection. The result is that the Atlantic overturning circulation in the model is very sensitive to the surface heat flux in the Labrador Sea and hence is correlated with the North Atlantic Oscillation. As strong subpolar convection is found in other models, we discuss broader implications.     

Numerical Considerations for Lagrangian Stochastic Dispersion Models: Eliminating Rogue Trajectories,and the Importance of Numerical Accuracy

When Lagrangian stochastic models for turbulent dispersion are applied to complex atmospheric flows, some type of ad hoc intervention is almost always necessary to eliminate unphysical behaviour in the numerical solution. Here we discuss numerical strategies for solving the non-linear Langevin-based particle velocity evolution equation that eliminate such unphysical behaviour in both Reynolds-averaged and large-eddy simulation applications. Extremely large or ‘rogue’ particle velocities are caused when the numerical integration scheme becomes unstable. Such instabilities can be eliminated by using a sufficiently small integration timestep, or in cases where the required timestep is unrealistically small, an unconditionally stable implicit integration scheme can be used. When the generalized anisotropic turbulence model is used, it is critical that the input velocity covariance tensor be realizable, otherwise unphysical behaviour can become problematic regardless of the integration scheme or size of the timestep. A method is presented to ensure realizability, and thus eliminate such behaviour. It was also found that the numerical accuracy of the integration scheme determined the degree to which the second law of thermodynamics or ‘well-mixed condition’ was satisfied. Perhaps more importantly, it also determined the degree to which modelled Eulerian particle velocity statistics matched the specified Eulerian distributions (which is the ultimate goal of the numerical solution). It is recommended that future models be verified by not only checking the well-mixed condition, but perhaps more importantly by checking that computed Eulerian statistics match the Eulerian statistics specified as inputs.     

Alternating field characteristics of pseudo-single-domain (2–14 μm) and multidomain magnetite

We report normalized AF demagnetization curves of anhysteretic remanences (ARM's) produced by 1-, 10- and 40-Oe steady fields and of saturation isothermal remanence (IRM_s) in a suite of dispersed, unannealed magnetite powders with median sizes of 2, 4, 6, 10 and 14 μm (pseudo-single-domain or PSD size range) and 100 μm (multidomain or MD size). Interpreted in the light of the domain structure test first proposed by Lowrie and Fuller [12], the relative stability trend of curves for the 2 μm sample is of single-domain (SD) type, the 1-Oe ARM being most resistant to demagnetization followed by the 10-Oe and 40-Oe ARM's and IRM_s. For the 100-μm sample, the trend is exactly reversed and is of MD-type. In the 4–14 μm samples, hitherto undescribed transitional trends between SD-type and MD-type occur. At 6 μm, 1-Oe, 10-Oe and 40-Oe ARM's preserve an SD-type trend but for all AF's > 75 Oe, IRM_s is more resistant than any of these remanences. At 10 μm, this trend is unmistakable, and only at 14 μm do the 1-Oe, 10-Oe and 40-Oe ARM curves merge. We conclude (1) that the Lowrie-Fuller test distinguishes between small MD grains enhanced by PSD remanence and large MD grains lacking PSD remanence, rather than between SD and MD structures per se, and (2) that in the PSD transition region from 6 to 14 μm in magnetite, IRM_s changes over to MD-type relative stability around 6 μm, whereas 10-Oe and 40-Oe ARM's achieve an MD-type trend around 14 μm, in accord with the predicted field dependence of the PSD threshold size.Our theoretical interpretation assumes that the intrinsic (internal field) coercive force spectra of weak-field and strong-field remanences are identical but that the observed (external field) spectrum is shifted to lower fields as a result of the internal demagnetizing field — NJ_r of the remanence J_r. The effect is slight for weak-field J_r's but substantial for IRM_s. Since all coercivities, high as well as low, are shifted, the result of the Lowrie-Fuller test is determined simply by the shape of the intrinsic coercivity spectrum or the corresponding AF demagnetization curve. Depending on the model of self-demagnetization used, either subexponential or sublinear AF decay curves of weak-field remanence will automatically lead to an MD-type trend, whereas by either model the decay curves that characterize SD and PSD remanences (decaying slowly initially and then more rapidly) will always produce and SD-type trend.