Some comments on the Oort cloud

The arguments used by Lyttleton to prove the nonexistence of the Oort cloud are reviewed, and it is shown that Oort's hypothesis remains consistent with observation. The 1950 model of the cloud cannot be correct and, by use of the results from a number of more recent papers, an improved model is described and compared with observations. It is emphasized that comparison of the predictions of theory with observations should concentrate on thea-distribution, as the 1/a-distribution masks much of the detailed structure of the theory. An order of magnitude argument is given which shows that 20% of so-called new comets have passed through the planetary system before, and the implications of this to the statistics of near-parabolic comet orbits are briefly investigated.     

Reply to the comment by Addo Van Pul and John H. Van Boxel on ‘a first order closure scheme to describe counter gradient momentum transport in plant canopies’ by Li,Miller and Lin

Correspondence

Reply to the comment by Addo Van Pul and John H. Van Boxel on a first order closure scheme to describe counter gradient momentum transport in plant canopies by Li, Miller and Lin     

Monazite U-Pb dating of staurolite grade metamorphism in pelitic schists

A study of the occurrence of and relations between rare-earth element (REE) minerals in pelitic schists indicates that monazite forms at or near the P and T of the staurolite isograd. Samples at staurolite grade from the Silurian Perry Mountain Formation in the Rumford quadrangle of Maine yield monazite in sufficient quantities to permit accurate dating of the metamorphic events forming the monazites. The bulk chemistry of the metapelites, as seen in the major element abundances and REE patterns, does not vary significantly across the study area. Thus the appearance and disappearance of REE phases is assumed to reflect changes in metamorphic grade. In a sample from the biotite zone, scanning electron microscope and microprobe studies show allanite and monazite intimately associated on a 10 m scale. The texture suggest that metastable detrital monazite breaks down, distributing its REE components to allanite. From samples below staurolite grade in which monazite is not present, our observations suggest that REEs are partitioned into allanite. At or near the staurolite isograd monazite forms as a metamorphic mineral, initiating its role as a geochronometer. Garnet-biotite geothermometry on samples at this grade from this and other studies places constraints on the minimum temperature necessary to form monazite: 525° C±25°C at 3.1±0.25 kbar. A total of 15 separates from nine schist samples ranging up to sillimanite grade have been dated. Each date is remarkably concordant, even though petrologic and textural studies by previous workers have shown that the rocks in the area have been affected by at least three metamorphic episodes. Calculations indicate insignificant Th disequilibrium in these monazites. The conditions associated with the metamorphic events suggest that monazite remains closed to lead loss provided that subsequent metamorphisms are at or below sillimanite grade. Two distinct metamorphic events are resolved, one at around 400 Ma and one at about 370 Ma. The latter was due to thermal effects of a nearby pluton that yields concordant monazite ages of 363 Ma. This work suggests that in addition to dating plutonism and high-grade metamorphism, monazite should be viewed as a reliable geochronometer for moderate metamorphism of pelitic schists.     

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.     

Gas maturity and alteration systematics across the Western Canada Sedimentary Basin from four mud gas isotope depth profiles

Carbon isotope and molecular compositions of Mississippian to Upper Cretaceous mud gases have been examined from four depth profiles across the Western Canada Sedimentary Basin (WCSB). The profiles range from the shallow oil sands in the east (R₀ = 0.25) to the very mature sediments in the overthrust zone to the west (R₀ = 2.5). In the undisturbed WCSB, δ¹³C₁–δ¹³C₂ and δ¹³C₂–δ¹³C₃ cross-plots show three maturity and alteration trends: (1) pre-Cretaceous gas sourced from type II kerogen; (2) Cretaceous Colorado Group gas; and (3) Lower Cretaceous Mannville Group biodegraded gas. A fourth set of distinctly different maturity trends is recognized for Lower Cretaceous gas sourced from type III kerogen in the disturbed belt of the WCSB. Displacement of these latter maturity trends to high δ¹³C₂ values suggests that the sampled gas was trapped after earlier formed gas escaped, probably as a result of overthrusting. Unusually ¹³C-enriched gas (δ¹³C₁ = −34‰, δ¹³C₂ = −13‰, and δ¹³C₃ = 0‰), from the Gething Formation in the disturbed belt, is the result of late stage gas cracking in a closed system. In general, gas maturity is consistent with the maturity of the host sediments in the WCSB, suggesting that migration and mixing of gases was not pervasive on a broad regional and stratigraphic scale. The ‘Deep Basin’ portion of the WCSB is an exception. Here extensive cross-formational homogenization of gases has occurred, in addition to updip migration along the most permeable stratigraphic units.     

Bringing NEOs into focus

Duncan Steel and Mark Bailey report on the well-attended December 2001 G Discussion Meeting whose subject was near-Earth objects.     

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.     

Dublin days: the 2003 NAM

Mark E Bailey and Tom Ray give an overview of the National Astronomy Meeting in Dublin, where record numbers of astronomers went to talk, listen and network.     

Orbitolinidae and Alveolinidae (Foraminiferida) from the uppermost Albian-lower Cenomanian of Monti d'Ocre (Abruzzi, Italy)

A biostratigraphic study carried out in the Monti d'Ocre area, Abruzzi, Central Apennines, allowed us to recognize Orbitolina (Conicorbitolina) moulladei ‘Strata 5 (1985) 1’, Praealveolina iberica Reichel and Praealveolina simplex Reichel in uppermost Albian–lower Cenomanian shelf-edge deposits of the Fossato Machè succession. These foraminifers have now been found for the first time in the Apennines of central Italy; their finding is quite important from a palaeobiogeographic viewpoint, as it contributes to the improvement of our knowledge on facies distribution in the circum-Mediterranean regions during the Cretaceous Period. In the study area, the coeval Monte Rotondo and Monte Orsello sections also crop out; these are characterized by bauxite deposits and stratigraphic gaps reflecting episodes of emergence on the carbonate platform. The Monte Rotondo and Monte Orsello sections accumulated in a platform back-reef environment; consequently, in this sector of the Monti d'Ocre area, the depositional environment shifted from a back-reef westward and southward to a shelf-edge northward, during the latest Albian–early Cenomanian. Owing to synsedimentary tectonics, the area investigated underwent differential subsidence: westward and southward, wide areas were uplifted and subjected to emergence, karstification and bauxite accumulation, whereas sedimentation continued in a shelf-edge environment in the north-eastern area.     

An accurate and robust multigrid algorithm for 2D forward resistivity modelling

We present an adaptation of the full multigrid algorithm in DC resistivity modelling in an effort to increase its accuracy. There is a great difficulty with conventional multigrid solvers in representing the physics of an arbitrary distribution of electrical conductivity on a very coarse grid. In general, conventional rectangular finite‐difference or 5‐point approximations of Poisson's equation cannot represent, at a coarse grid level, the effective anisotropy on a coarse scale which results from fine structure in the model. An exception to this generalization occurs when the principal axes of structural anisotropy are aligned with the coordinate axis. Additional and similarly generated problems arise when a coarse cell is obliged to represent fine structure containing very high conductivity contrasts. We have developed an adaptation of the usual resistive‐network representation of the continuum, which avoids some of these problems, and have compared it with the traditional resistive network currently used. The network adaptation consists of replacing the usual 5‐point Laplacian operator stencil used on the finite‐difference grid with a 9‐point stencil, and the conductivity scalar with a 6‐parameter conductivity parametrization. This parametrization permits representation of arbitrarily orientated anisotropy as well as more complex behaviour related to high conductivity contrasts. The importance of multigrid solvers does not lie in their speed at forward modelling (which is comparable with other methods), but rather in their potential for inverse modelling. Inverse solvers which proceed by refinement of an initially very coarse solution can, in principle, take time only linearly proportional to the number of gridpoints in the final desired model.