Environmental controls on clay mineralogy of an Early Jurassic mudrock (Blue Lias Formation,southern England)

We present a dataset including clay mineralogy, gamma-ray spectrometry, organic matter content and magnetic susceptibility of the Hettangian to lowest Sinemurian successions of Dorset and Somerset, southern UK (Blue Lias Formation, Bristol Channel and Wessex basins). In both areas, the clay assemblages comprise predominantly detrital illite, kaolinite and illite/smectite mixed layers. Clays probably originated from the erosion of the Hercynian massifs, the relative proportions of kaolinite and illite being modulated by arid-humid climatic fluctuations. The organic matter (OM) content (types II to IV) ranges up to 12% in both areas. A clear stratigraphical trend in clay mineral assemblages is apparent in Somerset, whereas in Dorset sharp contrasts between adjacent horizons and a greater dilution by carbonate mask the long-term evolution. Correlations between both areas based on similar vertical trends in clay mineral abundance support the suggestion of a hiatus within the angulata Zone of the Dorset succession. As expected, the kaolinite/illite ratio correlates with the Th/K ratio deduced from gamma-ray spectrometry. However, significant departures from the correlation occur in OM-rich intervals, suggesting that Th may be partly adsorbed on to OM particles. Surprisingly, high magnetic susceptibility correlates with abundant kaolinite, not with Fe-rich clays, indicating either that kaolinite is accompanied by a soil-inherited magnetisable phase (possibly iron oxide) or that illite-rich rocks are more strongly diluted by carbonate than are kaolinite-rich strata.     

An automated salt‐tracing gauge for flow‐velocity measurement

This article introduces the SVG (salt‐velocity gauge), a novel automated technique for measuring flow velocity by means of salt tracing. SVG allows a high measuring rate (up to one every 2 seconds), short control section length (down to 10 cm), high accuracy (+[sol ]?1·5 cm s^?1), and unbiased calculation of the mean velocity in experimental conditions with turbulent, supercritical flow. A few cubic centimetres of saturated salt solution (NaCl) are injected into the flow at regular time intervals using a programmable solenoid valve. The tracer successively passes two conductivity probes placed a short distance downstream. The transformation of the signal between the two probes is modelled as a one‐dimensional diffusion wave equation. Model calibration gives an estimation of the mean velocity and the diffusion for each salt plume. Two implementations of the SVG technique are described. The first was an outdoors simulated rainfall experiment in Senegal (conductivity probes at 40 cm apart, 8 Hz measurement rate, salt injections at 10 second intervals). Mean velocity was estimated to range between 0·1 and 0·3 m s^?1. The second was a laboratory‐based flume experiment (conductivity probes at 10 cm apart, 32 Hz, salt injections at 2 second intervals). Another SVG with probes at 34 cm apart was used for comparison. An acoustic Doppler velocimeter (ADV) was also used to give an independent assessment of velocity. Using the 10 cm salt gauge, estimated mean velocity ranged from 0·6 to 0·9 m s^?1 with a standard deviation of 1·5 cm s^?1. Comparisons between ADV, 10 cm SVG and 34 cm SVG were consistent and demonstrated that the salt‐tracing results were unbiased and independent of distance between probes. Most peaks were modelled with r² > 90 per cent. The SVG technology offers an alternative to the dye‐tracing technique, which has been severely criticized in the literature because of the wide interval of recommended values for the correction factor α to be applied to the timings. This article demonstrates that a fixed value of α is inappropriate, since the correction factor varies with velocity, diffusion and the length of the control section. Copyright © 2005 John Wiley & Sons, Ltd.     

Using radon-222 to study coastal groundwater/surface-water interaction in the Crau coastal aquifer (southeastern France)

Radon has been used to determine groundwater velocity and groundwater discharge into wetlands at the southern downstream boundary of the Crau aquifer, southeastern France. This aquifer constitutes an important high-quality freshwater resource exploited for agriculture, industry and human consumption. An increase in salinity occurs close to the sea, highlighting the need to investigate the water balance and groundwater behavior. Darcy velocity was estimated using radon activities in well waters according to the Hamada “single-well method” (involving comparison with radon in groundwater in the aquifer itself). Measurements done at three depths (7, 15 and 21 m) provided velocity ranging from a few mm/day to more than 20 cm/day, with highest velocities observed at the 15-m depth. Resulting hydraulic conductivities agree with the known geology. Waters showing high radon activity and high salinity were found near the presumed shoreline at 3,000 years BP, highlighting the presence of ancient saltwater. Radon activity has also been measured in canals, rivers and ponds, to trace groundwater discharges and evaluate water balance. A model of the radon spatial evolution explains the observed radon activities. Groundwater discharge to surface water is low in pond waters (4 % of total inputs) but significant in canals (55 l/m²/day).     

Coseismic and post-seismic signatures of the Sumatra 2004 December and 2005 March earthquakes in GRACE satellite gravity

The GRACE satellite mission has been measuring the Earth's gravity field and its temporal variations since 2002 April. Although these variations are mainly due to mass transfer within the geofluid envelops, they also result from mass displacements associated with phenomena including glacial isostatic adjustment and earthquakes. However, these last contributions are difficult to isolate because of the presence of noise and of geofluid signals, and because of GRACE's coarse spatial resolution (>400 km half-wavelength). In this paper, we show that a wavelet analysis on the sphere helps to retrieve earthquake signatures from GRACE geoid products. Using a wavelet analysis of GRACE geoids products, we show that the geoid variations caused by the 2004 December ( M _w= 9.2) and 2005 March ( M _w= 8.7) Sumatra earthquakes can be detected. At GRACE resolution, the 2004 December earthquake produced a strong coseismic decrease of the gravity field in the Andaman Sea, followed by relaxation in the area affected by both the Andaman 2004 and the Nias 2005 earthquakes. We find two characteristic timescales for the relaxation, with a fast variation occurring in the vicinity of the Central Andaman ridge. We discuss our coseismic observations in terms of density changes of crustal and upper-mantle rocks, and of the vertical displacements in the Andaman Sea. We interpret the post-seismic signal in terms of the viscoelastic response of the Earth's mantle. The transient component of the relaxation may indicate the presence of hot, viscous material beneath the active Central Andaman Basin.     

Trace element geochemistry of CR chondrite metal

We report trace element analyses by laser ablation inductively coupled plasma mass spectrometry (LA‐ICP‐MS) of metal grains from nine different CR chondrites, distinguishing grains from chondrule interior (“interior grains”), chondrule surficial shells (“margin grains”), and the matrix (“isolated grains”). Save for a few anomalous grains, Ni‐normalized trace element patterns are similar for all three petrographic settings, with largely unfractionated refractory siderophile elements and depleted volatile Au, Cu, Ag, S. All three types of grains are interpreted to derive from a common precursor approximated by the least‐melted, fine‐grained objects in CR chondrites. This also excludes recondensation of metal vapor as the origin of the bulk of margin grains. The metal precursors were presumably formed by incomplete condensation, with evidence for high‐temperature isolation of refractory platinum‐group‐element (PGE)‐rich condensates before mixing with lower temperature PGE‐depleted condensates. The rounded shape of the Ni‐rich, interior grains shows that they were molten and that they equilibrated with silicates upon slow cooling (1–100 K h^?1), largely by oxidation/evaporation of Fe, hence their high Pd content, for example. We propose that Ni‐poorer, amoeboid margin grains, often included in the pyroxene‐rich periphery common to type I chondrules, result from less intense processing of a rim accreted onto the chondrule subsequent to the melting event recorded by the interior grains. This means either that there were two separate heating events, which formed olivine/interior grains and pyroxene/margin grains, respectively, between which dust was accreted around the chondrule, or that there was a single high‐temperature event, of which the chondrule margin records a late “quenching phase,” in which case dust accreted onto chondrules while they were molten. In the latter case, high dust concentrations in the chondrule‐forming region (at least three orders of magnitude above minimum mass solar nebula models) are indicated.     

Constraints in 3D gravity inversion

A three-dimensional (3D) inversion program is developed to interpret gravity data using a selection of constraints. This selection includes minimum distance, flatness, smoothness and compactness constraints, which can be combined using a Lagrangian formulation. A multigrid technique is also implemented to resolve separately large and short gravity wavelengths. The subsurface in the survey area is divided into rectangular prismatic blocks and the problem is solved by calculating the model parameters, i.e. the densities of each block. Weights are given to each block depending on depth, a priori information on density and the density range allowed for the region under investigation. The present computer code is tested on modelled data for a dipping dike and multiple bodies. Results combining different constraints and a weight depending on depth are shown for the dipping dike. The advantages and behaviour of each method are compared in the 3D reconstruction. Recovery of geometry (depth, size) and density distribution of the original model is dependent on the set of constraints used. From experimentation, the best combination of constraints for multiple bodies seems to be flatness and a minimum volume for the multiple bodies. The inversion method is tested on real gravity data from the Rouyn-Noranda (Quebec) mining camp. The 3D inversion model for the first 10 km is in agreement with the known major lithological contacts at the surface; it enables the determination of the geometry of plutons and intrusive rocks at depth.     

Volcanic edifice stability during cryptodome intrusion

Limit equilibrium analyses were applied to the 1980 Mount St. Helens and 1956 Bezymianny failures in order to examine the influence on stability of structural deformation produced by cryptodome emplacement. Weakening structures associated with the cryptodome include outward-dipping normal faults bounding a summit graben and a flat shear zone at the base of the bulged flank generated by lateral push of the magma. Together with the head of the magmatic body itself, these structures serve directly to localize failure along a critical surface with low stability deep within the interior of the edifice. This critical surface, with the safety coefficient reduced by 25-30%, is then very sensitive to stability condition variation, in particular to the pore-pressure ratio (r_u) and seismicity coefficient (n). For r_u=0.3, or n=0.2, the deep surface suffers catastrophic failure, removing a large volume of the edifice flank. In the case of Mount St. Helens, failure occurred within a material with angle of friction ~40°, cohesion in the range 10⁵-10⁶ Pa, and probably significant water pore pressure. On 18 May 1980, detachment of slide block I occurred along a newly formed rupture surface passing through the crest of the bulge. Although sliding of block I may have been helped by the basal shear zone, significant pore pressure and a triggering earthquake were required (r_u=0.3 and n=0.2). Detachment of the second block was guided by the summit normal fault, the front of the cryptodome, and the basal shear zone. This occurred along a deep critical surface, which was on the verge of failure even before the 18 May 1980 earthquake. The stability of equivalent surfaces at Bezymianny Volcano appears significantly higher. Thus, although magma had already reached the surface, weaker materials, or higher pore pressure and/or seismic conditions were probably required to reach the rupture threshold. From our analysis, we find that deep-seated sector collapses formed by removing the edifice summit cannot generally result from a single slide. Cryptodome-induced deformation does, however, provide a deep potential slip surface. As previously thought, it may assist deep-seated sector collapse because it favors multiple retrogressive slides. This leads to explosive depressurization of the magmatic and hydrothermal systems, which undermines the edifice summit and produces secondary collapses and explosive blasts.