Numerical simulation of the final stages of terrestrial planet formation

Three representative numerical simulations of the growth of the terrestrial planets by accretion of large protoplanets are presented. The mass and relative-velocity distributions of the bodies in these simulations are free to evolve simultaneously in response to close gravitational encounters and occasional collisions between bodies. The collisions between bodies, therefore, arise in a natural way and the assumption of expressions for the relative velocity distribution and the gravitational collision cross section is unnecessary. These simulations indicate that the growth of bodies with final masses approaching those of Venus and the Earth is possible, at least for the case of a two-dimensional system. Simulations assuming an initial uniform distribution of orbital eccentricities on the interval from 0 to e_max are found to produce final states containing too many bodies with masses which are too small when e_max < 0.10, while simulations with e_max > 0.20 result in too many catastrophic collisions between bodies thus preventing rapid accretion of planetary-size bodies. The e_max = 0.15 simulation ends with a state surprisingly similar to that of the present terrestrial planets and, therefore, provides a rough estimate of the range of radial sampling to be expected for the terrestrial planets.     

Stable isotope (O and C) and pollen trends in eastern Lake Erie,evidence for a locally-induced climatic reversal of Younger Dryas age in the Great Lakes basin

A cool period from about 11000 to 10 500 BP (11 to 10.5 ka) is recognized in pollen records from the southern Great Lakes area by the return of Picea and Abies dominance and by the persistence of herbs. The area of cooling appears centred on the Upper Great Lakes. A high-resolution record (ca. 9 mm/y) from a borehole in eastern Lake Erie reveals, in the same time interval, this pollen anomaly, isotope evidence of meltwater presence (a — 3 per mil shift in ¹⁸O and a +1.1 per mil shift in ¹³C), increased sand, and reduced detrital calcite content, all suggesting concurrent cooling of Lake Erie. The onset of cooling is mainly attributed to the effect of enhanced meltwater inflow on the relatively large upstream Main Lake Algonquin during the first eastward discharge of glacial Lake Agassiz. Termination of the cooling coincides with drainage of Lake Algonquin, and is attributed to loss of its cooling effectiveness associated with a substantial reduction in its surface area. It is hypothesized that the cold extra inflow effectively prolonged the seasonal presence of lake ice and the period of spring overturn in Lake Algonquin. The deep mixing would have greatly increased the thermal conductive capacity of this extensive lake, causing suppression of summer surface lakewater temperatures and reduction of onshore growing-degree days. Alternatively, a rapid flow of meltwater, buoyed on sediment-charged (denser) lakewater, may have kept the lake surface cold in summer. Other factors such as wind-shifted pollen deposition and possible effects from the Younger Dryas North Atlantic cooling could have contributed to the Great Lakes climatic reversal, but further studies are needed to resolve their relative significance.Contribution to Climo Locarno — Past and Present Climate Dynamics; Conference September 1990, Swiss Academy of Sciences — National Climate ProgramGeological Survey of Canada Contribution 58 890     

Noble gases in the Allende and Abee meteorites and a gas-rich mineral fraction: investigation by stepwise heating

Noble gases in three meteoritic samples were examined by stepwise heating, in an attempt to relate peaks in the outgassing curves to specific minerals: NeKrXe in Allende (C3V) and an Allende residue insoluble in HF-HCl, and Xe in Abee (E4). In Allende, chromite and carbon contain most of the trapped Ne (²⁰Ne/²²Ne ≈ 8.7) and anomalous Xe enriched in light and heavy isotopes, and release it at ~850°C (bulk meteorite) or 1000°C (residue). Mineral Q, containing most of the trapped Ar, Kr, Xe as well as some Ne (²⁰Ne/²²Ne ≈ 10.4), releases its gases mainly between 1200 and 1600°C, well above the release temperatures of organic polymers (300–500°) or amorphous carbon (800–1000°). The high noble-gas release temperature, ready solubility in oxidizing acids, and correlation with acid-soluble Fe and Cr all point to an inorganic rather than carbonaceous nature of Q.All the radiogenic ¹²⁹Xe is contained in HCl, HF-soluble minerals, and is distributed as follows over the peaks in the release curve: Attend 1000° (75%), 1300° (25%); Abee (data of Hohenberg and Reynolds, 1969) ~850° (15%), 1100° (60%), 1300° (25%). No conclusive identifications of host phases can yet be given; possible candidates are troilite and silicates for Allende, and djerfisherite, troilite and silicates for Abee.Mineral Q strongly absorbs air xenon, and releases some of it only at 800–1000°C. Dilution by air Xe from Q and other minerals may explain why temperature fractions from bulk meteorites often contain less ^124–130Xe for a given enrichment in heavy isotopes than does xenon from etched chromitecarbon samples, although chromite-carbon is the source of the anomalous xenon in either case. Air xenon contamination thus is an important source of error in the derivation of fission xenon spectra.     

Evolution of the Bonneville shoreline scarp in west-central Utah: Comparison of scarp-analysis methods and implications for the diffusion model of hillslope evolution

Wave-cut pluvial shoreline scarps are ideal natural experiments in hillslope evolution because the ages of these scarps are often precisely known and because they form with a range of heights, alluvial textures, and microclimates (i.e., orientation). Previous work using midpoint-slope methods on pluvial scarps in the Basin and Range concluded that scarp evolution is nonlinear and microclimatically controlled. The purpose of this study was to further examine the influence of scarp height, texture and microclimate in an attempt to calibrate a nonlinear model of scarp evolution. To do this, over 150 profiles of the Bonneville shoreline in the adjacent Snake and Tule Valleys, west-central Utah were collected and analyzed by fitting the entire scarp profile to diffusion-equation solutions, taking into account uncertainty in the initial scarp angle. In contrast to previous studies, this analysis revealed no evidence for nonlinearity or microclimatic control. To understand the reason for this discrepancy, we undertook a systematic study of the accuracy of each scarp-analysis method. The midpoint-slope-inverse method was found to yield biased results, with systematically higher diffusion ages for young, tall scarps. The slope-offset method is unbiased but has limited resolution because it requires many scarp profiles to yield a single diffusion age. A method that incorporates the full scarp profile and uncertainty in the initial scarp angle was found to be the most accurate technique. The application of the full-scarp method to the Bonneville shoreline supports the use of a linear diffusion model for scarps up to 20 m in height. Scarp orientation had no discernable effect on diffusivity values. Soil texture was found to have a weak but significant inverse relationship with diffusivity values.     

Virtual vegetables and adopted sheep: ethical relation, authenticity and Internet-mediated food production technologies

Associating ideas related to cultural geographies of 'quality' food production and consumption with recent discussion of Internet technologies, virtuality and simulation, this paper examines two newly constituted Internet enterprises offering customers 'virtual' experiences of food production, and 'quality' food products in which they have invested. The paper shows how these enterprises are associated with particular types of geographical and ethical relation, and a search for 'authentic', connected, relationships with food, farming, locality and 'nature'. These relationships are constituted in an assemblage of things held together by flows of food, products, money and electronic communication. The paper ends by suggesting a number of theoretical ideas that might usefully be explored in empirical research into Internet-mediated food production.     

Mineralogy and origin of stichtite in chromite-bearing serpentinites

Stichtite, a rare (14 known localities worldwide) hydrated carbonate-hydroxide of Mg and Cr with ideal formula Mg₆Cr₂ (OH)₁₆ CO₃ · 4H₂O, occurs exclusively in Cr-rich serpentinites of ophiolites or greenstone belts. Physical properties (hardness = 1.5–2, specific gravity = 2.16–2.2, perfect basal [0001] cleavage, grain size commonly < 100 μm) resemble talc, but the mineral has an attractive purple to lilac color; chemical analyses demonstrate it to be a non-silicate. Stichtite generally occurs as irregular to rounded masses (< 1 cm – 30 cm across) and as veinlets (< 1 mm – > 2 cm wide) within serpentinite. Macroscopic and microscopic textures, such as crosscutting veinlets and stringers, demonstrate that stichtite formation invariably post-dated serpentinization. In some specimens stichtite surrounds relict grains of Cr-rich spinel; in others stichtite has completely replaced euhedral or subhedral chromites. Chemical analyses of stichtites reveal substantial substitution of Al and Fe³⁺ for Cr in specimens from many localities, reflecting a possible compositional continuum between stichtite and rhombohedral polymorphs hydrotalcite (Mg₆Al₂ (OH)₁₆ CO₃ · 4H₂O) and pyroaurite (Mg₆Fe₂ (OH)₁₆ CO₃ · 4H₂O). We report the first electron microprobe analyses of stichtites from seven localities, and summarize all available published chemical data. Stichtites very likely inherited part of their trivalent cation chemistry from precursor Cr-rich spinels, but stichtite growth apparently post-dated characteristic “ferritchromit” alteration, as demonstrated by the depletion of Al and enrichment in Fe³⁺ in stichtite relative to primary chromite core compositions. Stichtite appears to form by reaction between serpentine and altered chromite, during addition of substantial fluid, either as separate H₂O and CO₂ phases, or as a mixed volatile phase. Such reactions must involve removal of substantial SiO₂, possibly by transport and remote deposition of silica by throughgoing aqueous and carbonic fluid. Received: 4 April 1996 / Accepted: 16 September 1996     

Time-dependent, coupled, Ekman boundary layer solutions incorporating Stokes drift

Recent technological advances in current measuring devices has resulted in a large observational database related to wind-driven motions in the upper ocean mixed layer. This has served to highlight the fact that transient motions make up a substantial contribution of the resulting Ekman currents. At the same time, certain discrepancies have emerged between the observed angular deflections of the steady-state currents from the surface wind stress, both at the surface and at sub-surface depths, which cannot be reconciled using the classical Ekman model. This paper seeks to tackle these two issues.First a general analytical method is presented for solving the time dependent horizontal momentum Ekman equations. Analysis of the unsteady terms that arise from simple special cases shows how the evolution proceeds through three stages. At early times, the Coriolis acceleration is insignificant, and the current is unidirectional and deepens through downward diffusion of momentum. Later Coriolis acceleration deflects the current vectors in the upper layers, whilst downward diffusion of momentum continues to deepen the layer. Finally, once diffusion has penetrated down to the depth of the steady-state current, then the transients decay on the inertial or diffusive timescale, depending upon the boundary conditions of the particular problem.In the second half of the paper, a new steady-state model is developed that includes the effects of wind-generated waves, through the action of their Stokes drift on the planetary vorticity. Comparisons between observations and the theoretical predictions, demonstrate that inclusion of the Stokes drift is the key to reconciling the discrepancies in the angular deflections of the steady-state currents. This leads to the conclusion that Ekman layer currents are significantly influenced by the surface waves.