Instrumentation Development for In Situ 40Ar/39Ar Planetary Geochronology

The chronology of the Solar System, particularly the timing of formation of extra‐terrestrial bodies and their features, is an outstanding problem in planetary science. Although various chronological methods for in situ geochronology have been proposed (e.g., Rb‐Sr, K‐Ar), and even applied (K‐Ar), the reliability, accuracy, and applicability of the ⁴⁰Ar/³⁹Ar method makes it by far the most desirable chronometer for dating extra‐terrestrial bodies. The method however relies on the neutron irradiation of samples, and thus a neutron source. Herein, we discuss the challenges and feasibility of deploying a passive neutron source to planetary surfaces for the in situ application of the ⁴⁰Ar/³⁹Ar chronometer. Requirements in generating and shielding neutrons, as well as analysing samples are described, along with an exploration of limitations such as mass, power and cost. Two potential solutions for the in situ extra‐terrestrial deployment of the ⁴⁰Ar/³⁹Ar method are presented. Although this represents a challenging task, developing the technology to apply the ⁴⁰Ar/³⁹Ar method on planetary surfaces would represent a major advance towards constraining the timescale of solar system formation and evolution.     

Statistical analysis of galaxy surveys – III. The non-linear clustering of red and blue galaxies in the 2dFGRS

We present measurements of the higher order clustering of red and blue galaxies as a function of scale and luminosity made from the two-degree field galaxy redshift survey (2dFGRS). We use a counts-in-cells analysis to estimate the volume-averaged correlation functions,     , as a function of scale up to the order of   p = 5  , and also the reduced void probability function. Hierarchical amplitudes are constructed using the estimates of the correlation functions:     . We find that (i) red galaxies display stronger clustering than blue galaxies at all orders measured; (ii) red galaxies show values of   S _p   that are strongly dependent on luminosity whereas blue galaxies show no segregation in   S _p   within the errors; this is remarkable given the segregation in the variance; (iii) the linear relative bias shows the opposite trend to the hierarchical amplitudes, with little segregation for the red sequence and some segregation for the blue; (iv) faint red galaxies deviate significantly from the 'universal' negative binomial reduced void probabilities followed by all other galaxy populations. Our results show that the characteristic colour of a galaxy population reveals a unique signature in its spatial distribution. Such signatures will hopefully further elucidate the physics responsible for shaping the cosmological evolution of galaxies.     

Scales and causes of heterogeneity in bars in a large multi‐channel river: Río Paraná, Argentina

To date, published studies of alluvial bar architecture in large rivers have been restricted mostly to case studies of individual bars and single locations. Relatively little is known about how the depositional processes and sedimentary architecture of kilometre‐scale bars vary within a multi‐kilometre reach or over several hundreds of kilometres downstream. This study presents Ground Penetrating Radar and core data from 11, kilometre‐scale bars from the Río Paraná, Argentina. The investigated bars are located between 30 km upstream and 540 km downstream of the Río Paraná – Río Paraguay confluence, where a significant volume of fine‐grained suspended sediment is introduced into the network. Bar‐scale cross‐stratified sets, with lengths and widths up to 600 m and thicknesses up to 12 m, enable the distinction of large river deposits from stacked deposits of smaller rivers, but are only present in half the surface area of the bars. Up to 90% of bar‐scale sets are found on top of finer‐grained ripple‐laminated bar‐trough deposits. Bar‐scale sets make up as much as 58% of the volume of the deposits in small, incipient mid‐channel bars, but this proportion decreases significantly with increasing age and size of the bars. Contrary to what might be expected, a significant proportion of the sedimentary structures found in the Río Paraná is similar in scale to those found in much smaller rivers. In other words, large river deposits are not always characterized by big structures that allow a simple interpretation of river scale. However, the large scale of the depositional units in big rivers causes small‐scale structures, such as ripple sets, to be grouped into thicker cosets, which indicate river scale even when no obvious large‐scale sets are present. The results also show that the composition of bars differs between the studied reaches upstream and downstream of the confluence with the Río Paraguay. Relative to other controls on downstream fining, the tributary input of fine‐grained suspended material from the Río Paraguay causes a marked change in the composition of the bar deposits. Compared to the upstream reaches, the sedimentary architecture of the downstream reaches in the top ca 5 m of mid‐channel bars shows: (i) an increase in the abundance and thickness (up to metre‐scale) of laterally extensive (hundreds of metres) fine‐grained layers; (ii) an increase in the percentage of deposits comprised of ripple sets (to >40% in the upper bar deposits); and (iii) an increase in bar‐trough deposits and a corresponding decrease in bar‐scale cross‐strata (<10%). The thalweg deposits of the Río Paraná are composed of dune sets, even directly downstream from the Río Paraguay where the upper channel deposits are dominantly fine‐grained. Thus, the change in sedimentary facies due to a tributary point‐source of fine‐grained sediment is primarily expressed in the composition of the upper bar deposits.     

Convective instabilities in a variable viscosity fluid cooled from above

Whether in the mantle or in magma chambers, convective flows are characterized by large variations of viscosity. We study the influence of the viscosity structure on the development of convective instabilities in a viscous fluid which is cooled from above. The upper and lower boundaries of the fluid are stress-free. A viscosity dependence with depth of the form ν₀ + ν₁ exp(?γ.z) is assumed. After the temperature of the top boundary is lowered, velocity and temperature perturbations are followed numerically until convective breakdown occurs. Viscosity contrasts of up to 10⁷ and Rayleigh numbers of up to 10⁸ are studied.For intermediate viscosity contrasts (around 10³), convective breakdown is characterized by the almost simultaneous appearance of two modes of instability. One involves the whole fluid layer, has a large horizontal wavelength (several times the layer depth) and exhibits plate-like behaviour. The other mode has a much smaller wavelength and develops below a rigid lid. The “whole layer” mode dominates for small viscosity contrasts but is suppressed by viscous dissipation at large viscosity contrasts.For the “rigid lid” mode, we emphasize that it is the form of the viscosity variation which determines the instability. For steep viscosity profiles, convective flow does not penetrate deeply in the viscous region and only weak convection develops. We propose a simple method to define the rigid lid thickness. We are thus able to compute the true depth extent and the effective driving temperature difference of convective flow. Because viscosity contrasts in the convecting region do not exceed 100, simple scaling arguments are sufficient to describe the instability. The critical wavelength is proportional to the thickness of the thermal boundary layer below the rigid lid. Convection occurs when a Rayleigh number defined locally exceeds a critical value of 160–200. Finally, we show that a local Rayleigh number can be computed at any depth in the fluid and that convection develops below depth z_r (the rigid lid thickness) such that this number is maximum.The simple similarity laws are applied to the upper mantle beneath oceans and yield estimates of 5 × 10¹⁵?5 × 10¹⁶ m² s^?1 for viscosity in the thermal boundary layer below the plate.     

Mutual replacement reactions in alkali feldspars I: microtextures and mechanisms

Intracrystal microtextures formed by a process of mutual replacement in alkali feldspars record fluid–rock reactions that have affected large volumes of the Earth’s crust. Regular, ≤1 μm-scale ‘strain-controlled’ perthitic microtextures coarsen, by up to 10³, by a dissolution–reprecipitation process, producing microporous patch or vein perthites on scales >100 μm. We have developed earlier studies of such reactions in alkali feldspar cm-scale primocrysts in layered syenites from the Klokken intrusion, South Greenland. We present new hyperspectral CL, SEM images, and laser ICPMS analytical data, and discuss the mechanism of such replacement reactions. The feldspars grew as homogeneous sodic sanidines which unmixed and ordered by volume diffusion during cooling into the microcline field at ~450°C, giving regular, fully coherent ‘braid’ cryptoperthite. At ≤450°C the crystals reacted with a circulating post-magmatic aqueous fluid. The braid perthite behaved as a single reactant ‘phase’ which was replaced by two product phases, incoherent subgrains of low albite and microcline, with micropores at their boundaries. The driving force for the reactions was coherency strain energy, which was greater than the surface energy in the subgrain mosaic. The external euhedral crystal shapes and bulk major element composition of the primocrysts were unchanged but they became largely pseudomorphs composed of subgrains usually with the ‘pericline’ and ‘adularia’ habits (dominant {110} and subordinate {010} morphology) characteristic of low T growth. The subgrains have an epitactic relationship with parent braid perthite. Individual subgrains show oscillatory zoning in CL intensity, mainly at blue wavelengths, which correlates with tetrahedral Ti. Regular zoning is sometimes truncated by irregular, discordant surfaces suggesting dissolution, followed by resumption of growth giving regular zoning. Zones can be traced through touching subgrains, of both albite and microcline, for distances up to ~500 μm. At ≤340°C, the microcline subgrains underwent a third stage of unmixing to give straight lamellar film perthites with periodicities of ~1 μm, which with further cooling became semicoherent by the development of spaced misfit dislocations. Sub-grain growth occurred in fluid films that advanced through the elastically strained braid perthite crystals, which dissolved irreversibly. Braid perthite was more soluble than the strain-free subgrain mosaics which precipitated from the supersaturated solution. Some volumes of braid texture have sharp surfaces that suggest rapid dissolution along planes with low surface energies. Others have complex, diffuse boundaries that indicate a phase of coherent lamellar straightening by volume diffusion in response to strain relief close to a slowly advancing interface. Nucleation of strain-free subgrains was the overall rate-limiting step. To minimise surface energy subgrains grew with low energy morphologies and coarsened by grain growth, in fluid films whose trace element load (reflected in the oscillatory zoning) was dictated by the competitive advance of subgrains over a range of a few tens of mm. The cross-cutting dissolution surfaces suggest influxes of fresh fluid. Removal of feldspar to give 2 vol% porosity would require a feldspar:fluid ratio of ~1:26 (by wt). The late reversion to strain-controlled exsolution in microcline subgrains is consistent with loss of fluid above 340°C following depressurization of the intrusion. A second paper (Part II) describes trace element partitioning between the albite and microcline subgrains, and discusses the potential of trace elements as a low-T geothermometer. This paper and the Part II are dedicated in memory of J.V. Smith and W.L. Brown, both of whom died in 2007, in acknowledgement of their unrivalled contributions to the study of the feldspar minerals over more than half a century.     

The breakdown of potassium feldspar at high water pressures

The equilibrium position of the reaction between sanidine and water to form “sanidine hydrate” has been determined by reversal experiments on well characterised synthetic starting materials in a piston cylinder apparatus. The reaction was found to lie between four reversed brackets of 2.35 and 2.50 GPa at 450 °C, 2.40 and 2.59 GPa at 550 °C, 2.67 and 2.74 GPa at 650 °C, and 2.70 and 2.72 GPa at 680 °C. Infrared spectroscopy showed that the dominant water species in sanidine hydrate was structural H₂O. The minimum quantity of this structural H₂O, measured by thermogravimetric analysis, varied between 4.42 and 5.85 wt% over the pressure range of 2.7 to 3.2 GPa and the temperature range of 450 to 680 °C. Systematic variation in water content with pressure and temperature was not clearly established. The maximum value was below 6.07 wt%, the equivalent of 1 molecule of H₂O per formula unit. The water could be removed entirely by heating at atmospheric pressure to produce a metastable, anhydrous, hexagonal KAlSi₃O₈ phase (“hexasanidine”) implying that the structural H₂O content of sanidine hydrate can vary. The unit cell parameters for sanidine hydrate, measured by powder X-ray diffraction, were a = 0.53366 (±0.00022) nm and c = 0.77141 (±0.00052) nm, and those for hexasanidine were a = 0.52893 (±0.00016) nm and c = 0.78185 (±0.00036) nm. The behaviour and properties of sanidine hydrate appear to be analogous to those of the hydrate phase cymrite in the equivalent barium system. The occurrence of sanidine hydrate in the Earth would be limited to high pressure but very low temperature conditions and hence it could be a potential reservoir for water in cold subduction zones. However, sanidine hydrate would probably be constrained to granitic rock compositions at these pressures and temperatures. Received: 6 May 1997 / Accepted: 2 October 1997