Investigating the scale of structural controls on chlorinated hydrocarbon distributions in the fractured-porous unsaturated zone of a sandstone aquifer in the UK

Contaminant migration behaviour in the unsaturated zone of a fractured porous aquifer is discussed in the context of a study site in Cheshire, UK. The site is situated on gently dipping sandstones, adjacent to a linear lagoon historically used to dispose of industrial wastes containing chlorinated solvents. Two cores of more than 100 m length were recovered and measurements of chlorinated hydrocarbons (CHCs), inorganic chemistry, lithology, fracturing and aquifer properties were made. The results show that selecting an appropriate vertical sampling density is crucial both to providing an understanding of contaminant pathways and distinguishing whether CHCs are present in the aqueous or non-aqueous phase. The spacing of such sampling should be on a similar scale to the heterogeneity that controls water and contaminant movement. For some sections of the Permo-Triassic aquifer, significant changes in lithology and permeability occur over vertical distances of less than 1 m and samples need to be collected at this interval, otherwise considerable resolution is lost, potentially leading to erroneous interpretation of data. At this site, although CHC concentrations were high, the consistent ratio of the two main components of the plume (tetrachloroethene and trichloroethene) provided evidence of movement in the aqueous phase rather than in dense non-aqueous phase liquid (DNAPL).     

Comment on: 'The 3D shear experiment over the Natih field in Oman: the effects of fracture-filling fluids on shear propagation' by C.M. van der Kolk, W.S. Guest and J.H.H.M. Potters

We find errors in theory and application in the paper by van der Kolk et al. (Geophysical Prospecting 49 , 179–197 (2001)) that invalidate their arguments for the cause of the reduction in the velocity of the vertically propagating slower split shear-wave in the gas cap over the Natih field in Oman. We suggest that existing theories are adequate and can explain the anomaly.     

Discrete-element modelling of detachment folding

A two‐dimensional, discrete‐element modelling technique is used to investigate the initiation and growth of detachment folds in sedimentary rocks above a weak décollement level. The model depicts the sedimentary rocks as an assemblage of spheres that obey Newton's equations of motion and that interact with elastic forces under the influence of gravity. Faulting or fracturing between neighbouring elements is represented by a transition from repulsive–attractive forces to solely repulsive forces. The sedimentary sequence is mechanically heterogeneous, consisting of intercalated layers of markedly different strengths and thicknesses. The interlayering of weak and strong layers within the sedimentary rocks promotes the localization of flexural flow deformation within the weak layers. Even with simple displacement boundary conditions, and straightforward interlayering of weak and strong layers, the structural geometries that develop are complex, with a combination of box, lift‐off and disharmonic detachment fold styles forming above the décollement. In detail, it is found that the modelled folds grow by both limb rotation and limb lengthening. The combination of these two mechanisms results in uplift patterns above the folds that are difficult, or misleading, to interpret in terms of simple kinematic models. Comparison of modelling results with natural examples and with kinematic models highlights the complexities of structural interpretation in such settings.     

A numerical model of the partitioning of trace chemical solutes during drop freezing

Partitioning of volatile chemicals among the gas, liquid, and solid phases during freezing of liquid water in clouds can impact trace chemical distributions in the troposphere and in precipitation. We describe here a numerical model of this partitioning during the freezing of a supercooled liquid drop. Our model includes the time-dependent calculation of the coupled processes of crystallization kinetics, heat transport, and solute mass transport, for a freezing hydrometeor particle. We demonstrate the model for tracer partitioning during the freezing of a 1000 μm radius drop on a 100 μm ice substrate, under a few ambient condition scenarios. The model effectively simulates particle freezing and solute transport, yielding results that are qualitatively and quantitatively consistent with previous experimental and theoretical work. Results suggest that the ice shell formation time is governed by heat loss to air and not by dendrite propagation, and that the location of ice nucleation is not important to freezing times or the effective partitioning of chemical solutes. Even for the case of nucleation at the center of the drop, we found that dendrites propagated rapidly to form surface ice. Freezing then proceeded from the outside in. Results also indicate that the solid-liquid interfacial surface area is not important to freezing times or the effective partitioning of chemical solutes, and that the rate aspects of trapping are more important than equilibrium solid-liquid partitioning to the effective partitioning resulting from freezing.     

Electrodynamic Disaggregation: Does it Affect Apatite Fission-Track and (U-Th)/He Analyses?

Apatite fission-track and (U-Th)/He analyses require the liberation of intact idiomorphic apatite grains from rock samples. While routinely being carried out by mechanical methods, electrodynamic disaggregation (ED) offers an alternative approach. The high-voltage discharges produced during the ED process create localised temperature peaks (10000 K) along a narrow plasma channel. In apatite, such high temperatures could potentially reduce the length of fission tracks, which start to anneal at temperatures > 60 °C, and could also enhance He diffusion, which becomes significant at 30–40 °C over geological time scales. A comparison of fission-track analyses and (U-Th)/He ages of apatites prepared both by mechanical (jaw crusher, disk mill) and ED processing provides a way of determining whether heating during the latter method has any significant effect. Apatites from three samples of different geological settings (an orthogneiss from Madagascar, the Fish Canyon Tuff, and a muscovite-gneiss from Greece) yielded statistically identical track length distributions compared to samples prepared mechanically. Additionally, (U-Th)/He ages of apatites from a leucogranite from Morocco prepared by both methods were indistinguishable. These first results indicated that during electrodynamic disaggregation apatite crystals were not heated enough to partially anneal the fission tracks or induce significant diffusive loss of He.     

Carbonate cements in Miller field of the UK North Sea: a natural analog for mineral trapping in CO<Subscript>2</Subscript> geological storage

Miller field of the North Sea has had high concentrations of natural CO₂ for ~70 Ma. It is an ideal analog for the long-term fate of CO₂ during engineered storage, particularly for formation of carbonate minerals that permanently lock up CO₂ in solid form. The Brae Formation reservoir sandstone contains an unusually high quantity of calcite concretions; however, C and O stable isotopic signatures suggest that these are not related to the present-day CO₂ charge. Margins of the concretions are corroded, probably because of reduced pH due to CO₂ influx. Dispersed calcite cements are also present, some of which postdate the CO₂ charge and, therefore, are the products of mineral trapping. It is calculated that only a minority of the reservoired CO₂ in Miller (6–24%) has been sequestrated in carbonates, even after 70 Ma of CO₂ emplacement. Most of the CO₂ accumulation is dissolved in pore fluids. Therefore, in a reservoir similar to the Brae Formation, engineered CO₂ storage must rely on physical retention mechanisms because mineral trapping is both incomplete and slow.     

The 2dF Galaxy Redshift Survey: luminosity dependence of galaxy clustering

We investigate the dependence of the strength of galaxy clustering on intrinsic luminosity using the Anglo-Australian two degree field galaxy redshift survey (2dFGRS). The 2dFGRS is over an order of magnitude larger than previous redshift surveys used to address this issue. We measure the projected two-point correlation function of galaxies in a series of volume-limited samples. The projected correlation function is free from any distortion of the clustering pattern induced by peculiar motions and is well described by a power law in pair separation over the range     . The clustering of     galaxies in real space is well-fitted by a correlation length     and power-law slope     . The clustering amplitude increases slowly with absolute magnitude for galaxies fainter than M *, but rises more strongly at higher luminosities. At low luminosities, our results agree with measurements from the Southern Sky Redshift Survey 2 by Benoist et al. However, we find a weaker dependence of clustering strength on luminosity at the highest luminosities. The correlation function amplitude increases by a factor of 4.0 between     and −22.5, and the most luminous galaxies are 3.0 times more strongly clustered than L * galaxies. The power-law slope of the correlation function shows remarkably little variation for samples spanning a factor of 20 in luminosity. Our measurements are in very good agreement with the predictions of the hierarchical galaxy formation models of Benson et al.     

The Berkeley parasitic SETI Program

Parasitic programs for the Search for Extraterrestrial Intelligence (SETI), carried out concurrently with conventional radio astronomical observing programs, can be an attractive and cost-effective means of exploring the large multidimensional search space intrinsic to this effort. We describe a microprocessor-based automated SETI acquisition system which searches for and records spectra of narrowband signals in the IF band of an observatory receiver. Data taken with this system over 35 days at the Hat Creek Radio Observatory at 1612 MHz are discussed. Out of approximately 10⁵ spectra processed during this period, 4000 were identified by the system as containing narrowband signals and were recorded. Subsequent analysis indicates that over 3900 of these are due to local RF contamination. The remainder are undergoing further investigation.     

Revised recommended methods for analyzing crater size-frequency distributions

Impact crater populations help us to understand solar system dynamics, planetary surface histories, and surface modification processes. A single previous effort to standardize how crater data are displayed in graphs, tables, and archives was in a 1978 NASA report by the Crater Analysis Techniques Working Group, published in 1979 in Icarus. The report had a significant lasting effect, but later decades brought major advances in statistical and computer sciences while the crater field has remained fairly stagnant. In this new work, we revisit the fundamental techniques for displaying and analyzing crater population data and demonstrate better statistical methods that can be used. Specifically, we address (1) how crater size-frequency distributions (SFDs) are constructed, (2) how error bars are assigned to SFDs, and (3) how SFDs are fit to power-laws and other models. We show how the new methods yield results similar to those of previous techniques in that the SFDs have familiar shapes but better account for multiple sources of uncertainty. We also recommend graphic, display, and archiving methods that reflect computers’ capabilities and fulfill NASA's current requirements for Data Management Plans.