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.     

Sampling the vertical particle flux in the upper water column using a large diameter free-drifting NetTrap adapted to an Indented Rotating Sphere sediment trap

Further development of the large, surface-tethered sediment trap (NetTrap) employed as part of the MedFlux program is described whereby the large collection capacity of the NetTrap is combined with an Indented Rotating Sphere/Sample Carousel (IRSC) sediment trap (IRSC–NT). This trap is capable of collecting particle flux either in a time series or settling velocity mode; settling velocity mode allows the collection of particles that fall within discrete settling velocity intervals. During short field deployments in the Mediterranean Sea the IRSC–NT configured in the settling velocity mode successfully collected unpoisoned samples for chemical and microbiological experiments. In addition to the development of the IRSC–NT, particle-settling behavior above and below the swimmer-excluding IRS valve was tested during on-deck experiments using a specially constructed water-tight trap. Chemical analyses of settling materials (published elsewhere) suggested that separation of particles by settling velocity was achieved. However, due to the motion of the ship, it was not possible to directly measure particle-settling velocities within the trap. Particle release from the IRS did not bias the apparent settling velocity spectrum. Rotation of the IRS did not engender turbulence at the surface of the sphere or within the skewed funnel below. Tests of different ball designs over the course of the MedFlux program showed that a “ridge and saddle” pattern was optimal for efficiently transferring particles under the IRS seal while still reducing swimmer entrance to the collection funnel. The large size of the IRSC–NT did not prevent it from drifting effectively with the current. Several modifications of the present design are proposed that should improve the accuracy of the settling velocity measurements.     

Evidence for silicate dissolution on Mars from the Nakhla meteorite

Veins containing carbonates, hydrous silicates, and sulfates that occur within and between grains of augite and olivine in the Nakhla meteorite are good evidence for the former presence of liquid water in the Martian crust. Aqueous solutions gained access to grain interiors via narrow fractures, and those fractures within olivine whose walls were oriented close to (001) were preferentially widened by etching along [001]. This orientation selective dissolution may have been due to the presence within olivine of shock‐formed [001](100) and [001]{110} screw dislocations. The duration of etching is likely to have been brief, possibly less than a year, and the solutions responsible were sufficiently cool and reducing that laihunite did not form and Fe liberated from the olivine was not immediately oxidized. The pores within olivine were mineralized in sequence by siderite, nanocrystalline smectite, a Fe‐Mg phyllosilicate, and then gypsum, whereas only the smectite occurs within augite. The nanocrystalline smectite was deposited as submicrometer thick layers on etched vein walls, and solution compositions varied substantially between and sometimes during precipitation of each layer. Together with microcrystalline gypsum the Fe‐Mg phyllosilicate crystallized as water briefly returned to some of the veins following desiccation fracturing of the smectite. These results show that etching of olivine enhanced the porosity and permeability of the nakhlite parent rock and that dissolution and secondary mineralization took place within the same near‐static aqueous system.