ALBERT

Notes: Soil minerals are known to influence the biological stability of soil organic matter (SOM). Our study aimed to relate properties of the mineral matrix to its ability to protect organic C against decomposition in acid soils. We used the amount of hydroxyl ions released after exposure to NaF solution to establish a reactivity gradient spanning 12 subsoil horizons collected from 10 different locations. The subsoil horizons represent six soil orders and diverse geological parent materials. Phyllosilicates were characterized by X-ray diffraction and pedogenic oxides by selective dissolution procedures. The organic carbon (C) remaining after chemical removal of an oxidizable fraction of SOM with NaOCl solution was taken to represent a stable organic carbon pool. Stable organic carbon was confirmed as older than bulk organic carbon by a smaller radiocarbon (14C) content after oxidation in all 12 soils. The amount of stable organic C did not depend on clay content or the content of dithionite–citrate-extractable Fe. The combination of oxalate-extractable Fe and Al explained the greatest amount of variation in stable organic C (R2 = 0.78). Our results suggest that in acid soils, organic matter is preferentially protected by interaction with poorly crystalline minerals represented by the oxalate-soluble Fe and Al fraction. This evidence suggests that ligand exchange between mineral surface hydroxyl groups and negatively charged organic functional groups is a quantitatively important mechanism in the stabilization of SOM in acid soils. The results imply a finite stabilization capacity of soil minerals for organic matter, limited by the area density of reactive surface sites.

Notes: Organo-mineral associations stabilize soil organic matter, though the mechanisms by which they do so are unclear. We used particle-size fractions 〈 6.3 μm of two soils to examine the importance of Fe oxides, short-range order Al silicates and the surface areas of minerals and micropores on the formation of organo-mineral associations. In the subsoil Fe oxides were most strongly statistically correlated with the mineral-bound organic carbon. We therefore assume that they are the most important substrates for the formation of organo-mineral associations. There is no indication that this is caused by physical protection of organic matter in their micropores (〈 2 nm). In the Haplic Podzol, dithionite–citrate–bicarbonate-soluble short-range order Al silicates may also play a role. Fe oxide particles were calculated to offer specific surface areas of ∼ 200 m2 g−1 (goethite) and ∼ 800 m2 g−1 (ferrihydrite), corresponding to crystal diameters of only a few nm. We assume that the resulting large amount of oxide-specific reactive surface sites (conditionally charged hydroxyl groups) is responsible for their dominant role as sorbents. With maximum C loadings of 1.3 mg C per m2 Fe oxide for the Dystric Cambisol and 1.1 mg C per m2 Fe oxide + short-range order Al silicates for the Haplic Podzol the subsoils of both soils seem to have reached saturation with respect to organic matter sorption. In contrast to subsoil horizons, organo-mineral associations from topsoils contain much larger amounts of organic matter. Here a larger C loading on Fe oxides or a greater importance of other sorbents in addition to the oxides must be assumed.

Notes: If biospheric sinks, such as soil organic carbon, are to be used to meet obligations for greenhouse gas emission reduction under the Kyoto Protocol, the permanence of these sinks needs to be considered. Further, since only direct human-induced carbon sinks can be included, and sinks resulting from indirect and natural effects cannot be used, there is a pressing need to separate direct human-induced effects from indirect and natural effects. Since these effects also influence the permanence of soil organic stocks, this paper attempts to synthesize existing knowledge in soil science, and use models to examine the likely influence of direct, indirect and natural effects on the permanence of soil organic carbon stocks.

Notes: The formation of basic aluminium sulphate (BAS) minerals [(K,Na)nAlx(OH)y(SO4)z] has often been invoked to explain sulphate retention in soils. These minerals have not yet, however, been directly observed in the soil. We extracted the clay fractions of Andosols intercepting large inputs of volcanogenic sulphur dioxide and acids (HCl, HF), simply by dispersing clays with Na+-resins in deionized water without any other chemical treatment. Clay fractions concentrate 39–63% of total sulphur content of soil. Transmission electron microscopy coupled with energy-dispersive analysis revealed the presence of BAS particles, appearing as nodules and spheres. These particles have an equivalent diameter smaller than 0.2 µm. They have an Al:S ratio close to 2.2 and 3.8 and are possibly amorphous aluminite or basaluminite, respectively. They seem to have been formed in microenvironments enriched in sulphate, but also in fluoride anions. Their formation seems to have been enhanced by the combination of large inputs of acids and SO2 and an effective Al supply from weathering of volcanic glass.

Notes: Natural wetlands are a significant source of atmospheric methane, an important greenhouse gas. Compared with numerous papers on measurements of methane emission from natural wetland surfaces, there are few reports on methane configuration and distribution within wetland soil profiles. By using a newly designed gas sampler, we succeeded in collecting free-phase gas from beneath the water table down to 120 cm in a peat. The volumetric percentage of methane in the gas phase increased with depth and was generally more than 50% beneath the zone within which the water table fluctuates. The volume of the gas phase in the peat beneath the water table was estimated to be from 0 to 19% with significant variation with depth, suggesting uneven distribution of gas bubbles. Using the volume ratio of the gas and liquid phases and methane concentration data in the gas phase, as well as assuming that methane was in equilibrium (based on Henry's Law between the two phases), we calculated that ∼60% of the methane accumulates in the form of bubbles. These results suggest the importance of ebullition in methane emission, which might be a major cause for the reportedly large variation of methane emission in both space and time. Most importantly, our results show the need to consider gaseous-phase methane for understanding the production, transport and emission mechanisms of methane in wetlands, which has been overlooked to date.

Notes: Soil microbial biomass and microbial quotient (the ratio of soil microbial biomass to soil organic carbon) are considered to be useful as rapidly responding indicators of perturbations of soil properties. In this paper we will use a well-tested model (the continuous-quality theory) to analyse these variables in a Swedish 35-year-old field experiment with a black fallow, crop with no N addition, crop with calcium nitrate addition, and six treatments with organic amendments: straw, green manure, peat, farmyard manure, sawdust and sewage sludge.The model predicts correctly that the amount of microbial biomass increases for all the treatments with organic amendments compared with the black fallow treatment. The microbial biomass quotient increases also for all the amended treatments, except peat and sewage sludge, and decreases for the other treatments. The microbial biomass and microbial quotient increase with both the amounts of organic matter added (crop residues and amendments) and the quality of the added matter. However, to fully explain the observations it is also necessary to have an increasing microbial mortality with substrate quality. Moreover, short-term observations can be misleading with respect to both the magnitude and direction of long-term changes in biomass and related variables. Special attention must be paid to such amendments as sewage sludge, where contaminants such as heavy metals may determine process rates. We find no relation between microbial biomass or microbial quotient and yields.

Notes: The mobility of hydrophobic organic compounds (HOCs) in soils can be influenced by the presence of dissolved organic matter (DOM). While numerous studies have determined interactions of HOCs with humic and fulvic acids, only few data exist on the partitioning of HOCs to natural, non-fractionated DOM as it occurs in soil solutions. In this study, DOM was extracted from 17 soil samples with a broad range of chemical and physical properties, originating from different land uses. The partition coefficients of pyrene to DOM were determined in all soil extracts and for two commercial humic acids using the fluorescence quenching method. For the soil extracts, log KDOC values ranged from 3.2 to 4.5 litres kg−1. For the Aldrich and Fluka humic acids, log KDOC was 4.98 and 4.96 litres kg−1, respectively, thus indicating that they are not representative for soil DOM. After excluding these two values, the statistical analysis of the data showed a significant negative correlation between log KDOC and pH. This was also shown for one sample where the pH was adjusted to values ranging from 3 to 9. A multiple regression analysis suggested that ultraviolet absorbance at 280 nm (an indicator for aromaticity) and the E4:E6 ratio (an indicator for molecular weight) had additional effects on log KDOC. The results indicate that the partitioning of pyrene to DOM is reduced at alkaline pH, probably due to the increased polarity of the organic macromolecules resulting from the deprotonation of functional groups. Only within a narrow pH range was the KDOC of pyrene mainly related to the aromaticity of DOM.

Notes: This paper describes the characteristic shrinkage curve of a soil as it dries. It also introduces a new method to measure specific soil moisture states such as permanent wilting point and field capacity using the shrinkage curve. The model considers the soil fabric as a non-rigid, aggregated and unsaturated soil water medium. Its functional parameters (defined as pedohydral parameters) are the independent parameters of the shrinkage curve, which defines two major aspects of the structure of the fabric: (i) the volumetric functional elements of the soil fabric, and (ii) the arrangement of solid, water and air as functions of water content.Two African tropical soils were characterized and analysed according to a structural model with parameters of their shrinkage curve. The parameters of four soil horizons for the two soil types were determined continuously along their shrinkage curves. They were then used as descriptive variables representing the soil's hydrostructural behaviour in a canonical variate analysis, the results of which showed the A horizons to be distinct from the lower horizons, which appear as continua down the profiles. Analysis of the results reveals the importance of the clay and of the iron/clay ratio in the hydrostructural properties of the two soils. Transitional points of the shrinkage curve matched well some particular moisture states of the water potential curve, such as wilting points and field capacities. A more precise method of calculating water-holding characteristics and other structural properties such as air capacity and swelling index can be obtained from the characteristic shrinkage curve.

Notes: The organic matter (OM) of soils with andic properties has long been considered highly stable because of the presence of Al–humus complexes and sorption of organic ligands onto amorphous compounds. In this study, we characterized soils under different land use regimes located within an amphibolitic massif close to Santiago de Compostela (Spain), where soils with andic properties are present. Slash and burn agriculture was a common practice in the area until the second half of the 20th century. Thereafter, modern agriculture was progressively introduced into the area (AGR soils), and the rest of the land was either reforested or abandoned (FOR soils). We found that the mean organic C content of AGR soils (48.7 g kg−1) was ∼ 50% that of FOR soils (94.2 g kg−1). Mean soil pH was significantly greater (P 〈 0.05) in the AGR than in the FOR soils (4.95 compared with 4.63), which is attributed to liming and Ca-phosphate fertilization of the former. Mean concentrations of the Al forms studied (extractable with CuCl2, sodium pyrophosphate, ammonium oxalate, or NaOH) were significantly smaller (P 〈 0.01) in AGR (1.4, 4.9, 9.3, 11.0 g kg−1, respectively) than in FOR soils (3.9, 10.2, 16.5, 17.9 g kg−1, respectively). The results show the vulnerability of the OM and Al–humus complexes in these soils to modern agricultural practices, which has led to the attenuation – and in some cases even the disappearance – of andic soil properties in a relatively short time (〈 30 years) following changes in land use/management. We propose the inclusion of the formative element ‘andic’ in the criteria for the definition of Umbrisol subunits; this would avoid the abrupt discontinuity observed in the current World Reference Base classification.

Notes: Based on an acoustic assumption (that the shear-wave velocity is zero) and a dispersion relationship, we derive an acoustic wave equation for P-waves in tilted transversely isotropic (TTI) media (transversely isotropic media with a tilted symmetry axis). This equation has fewer parameters than an elastic wave equation in TTI media and yields an accurate description of P-wave traveltimes and spreading-related attenuation. Our TTI acoustic wave equation is a fourth-order equation in time and space. We demonstrate that the acoustic approximation allows the presence of shear waves in the solution. The substantial differences in traveltime and amplitude between data created using vertical transversely isotropic (VTI) and TTI assumptions is illustrated in examples.

Notes: A high-resolution P-wave seismic reflection survey was conducted in the area of Lambro Park within the city of Milan (northern Italy). Several high permeability channel sequences were identified and the stratigraphy of the infill was accurately mapped. The seismic signature imaged a previously undefined unit and also established the lateral correlation of some depositional units within the three major aquifer groups along the survey line.The complexity of the acoustic framework and the reduced depth of the source location limited effective elastic-wave generation. The geological setting also caused propagation of strong coherent noise patterns. Further interference, observed in the recorded data, was due to the traffic noise from the nearby highway. The attenuation of the undesired events required the design of specific filters and their multistep implementation.The results of forward modelling based on borehole information and of noise tests were crucial factors in the design of the processing parameters and in the stratigraphic interpretation of the final stacked section.

Notes: The subsurface geometry of the Sebastián Vizcaíno Basin is obtained from the 2D inversion of gravity data, constrained by a density-versus-depth relationship derived from an oil exploration deep hole. The basin accumulated a thick pile of marine sediments that evolved in the fore-arc region of the compressive margin prevalent along western North America during Mesozoic and Tertiary times. Our interpretation indicates that the sedimentary infill in the Sebastián Vizcaíno Basin reaches a maximum thickness of about 4 km at the centre of a relatively symmetric basin. At the location of the Suaro-1 hole, the depth to the basement derived from this work agrees with the drilled interface between calcareous and volcaniclastic members of the Alisitos Formation. A sensitivity analysis strongly suggests that the assumed density function leads to a nearly unique solution of the inverse problem.

Notes: We describe a method to invert a walkaway vertical seismic profile (VSP) and predict elastic properties (P-wave velocity, S-wave velocity and density) in a layered model looking ahead of the deepest receiver. Starting from Bayes's rule, we define a posterior distribution of layered models that combines prior information (on the overall variability of and correlations among the elastic properties observed in well logs) with information provided by the VSP data. This posterior distribution of layered models is sampled by a Monte-Carlo method. The sampled layered models agree with prior information and fit the VSP data, and their overall variability defines the uncertainty in the predicted elastic properties. We apply this technique first to a zero-offset VSP data set, and show that uncertainty in the long-wavelength P-wave velocity structure results in a sizable uncertainty in the predicted elastic properties. We then use walkaway VSP data, which contain information on the long-wavelength P-wave velocity (in the reflection moveout) and on S-wave velocity and density contrasts (in the change of reflectivity with offset). The uncertainty of the look-ahead prediction is considerably decreased compared with the zero-offset VSP, and the predicted elastic properties are in good agreement with well-log measurements.

Notes: The existence of non-zero reflectivity for zero-offset PS waves in horizontally layered media is discussed. Field measurements have suggested the occurrence of this phenomenon. We show that, in some cases, if anisotropy is considered in the problem, this reflectivity can be predicted. By using an approximated formulation to determine the coefficient of reflection in arbitrarily anisotropic media, it is clear that, for some elastic arrangements of the subsurface, reflected energy is associated with the converted wave resulting from normal P-wave incidence.

Notes: A simple scaling relationship is shown to facilitate comparison, correlation and integration of data recorded using the common experimental configurations in electrical and electromagnetic depth sounding. Applications of the scheme to field data from typical geological and landfill environments show that it is robust and, where transient electromagnetic (TEM) data are available, enables easy identification and quantification of electrical static shift (galvanic distortion) in magnetotelluric and direct current (DC) sounding curves. TEM-based procedures are suggested for both the direct removal of static shift in DC sounding curves and effective joint data inversion with the most-squares criterion in the presence of static shift. A case study of aquifer characterization using sounding data from borehole sites in the Vale of York in England shows that static shift is a common problem in this glacial-covered terrain and demonstrates the effectiveness of the proposed joint DC-TEM inversion strategy in handling distorted soundings.

Notes: A three-dimensional (3D) electrical resistivity modelling code is developed to interpret surface and subsurface data. Based on the integral equation, it calculates the charge density caused by conductivity gradients at each interface of the mesh, allowing the estimation of the potential everywhere without the need to interpolate between nodes. Modelling generates a huge matrix, made up of Green's functions, which is stored by using the method of pyramidal compression. The potential is compared with the analytical and the numerical solutions obtained by finite-difference codes for two models: the two-layer case and the vertical contact case. The integral method is more accurate around the source point and at the limits of the domain for the potential calculation using a pole-pole array. A technique is proposed to calculate the sensitivity (Jacobian) and Hessian matrices in 3D. The sensitivity is based on the derivative with respect to the block conductivity of the potential computed using the integral equation; it is only necessary to compute the electrical field at the source location. A direct extension of this technique allows the determination of the second derivatives. The technique is compared with the analytical solutions and with the calculation of the sensitivity according to the method using the inner product of the current densities calculated at the source and receiver points. Results are very accurate when the Green's function that includes the source image is used. The calculation of the three components of the electric field on the interfaces of the mesh is carried out simultaneously and quickly, using matrix compression.

Notes: In a sedimentary environment, layered models are often capable of representing the actual geology more accurately than smooth minimum structure models. Furthermore, interval thicknesses and resistivities are often the parameters to which non-geophysicist experts can relate and base decisions on when using them in waste site remediation, groundwater modelling and physical planning.We present a laterally constrained inversion scheme for continuous resistivity data based on a layered earth model (1D). All 1D data sets and models are inverted as one system, producing layered sections with lateral smooth transitions. The models are regularized through laterally equal constraints that tie interface depths and resistivities of adjacent layers. Prior information, e.g. originating from electric logs, migrates through the lateral constraints to the adjacent models, making resolution of equivalences possible to some extent. Information from areas with well-resolved parameters will migrate through the constraints in a similar way to help resolve the poorly constrained parameters. The estimated model is complemented by a full sensitivity analysis of the model parameters, supporting quantitative evaluation of the inversion result.Examples from synthetic 2D models show that the model recognition of a sublayered 2D wedge model is improved using the laterally constrained inversion approach when compared with a section of combined 1D models and when compared with a 2D minimum structure inversion. Case histories with data from two different continuous DC systems support the conclusions drawn from the synthetic example.

Notes: We analyse the problem of radiation of seismic waves by a vibroseis source when the baseplate is subject to flexure. A theoretical model is proposed to account for baseplate flexure, generalizing the well-known model of the vibroseis source of Sallas and Weber, which was developed for a rigid plate. Using the model proposed, we analyse the effect of flexure on the properties of seismic waves. We show that the flexure does not contribute to the far-field and mainly affects the readings of the reference accelerometer that is used to measure the force applied to the ground; these readings generally become dependent on the location of the sensor on the plate. For muddy and sandy soils, the effect of flexure on baseplate-acceleration measurements is nonetheless pronounced at the high end of the vibroseis frequency band only (∼100 Hz), and is negligible at all frequencies for stiffer soils. The corresponding phase lags introduced by the flexural vibrations at high frequencies lead to errors in the traveltime measurements (through the cross-correlation function) of up to 0.6 ms for muddy soils and less for denser soils. We show the existence of an optimal position of the reference sensor on the baseplate and also propose a general method of eliminating the phase lag due to the baseplate flexure in acceleration measurements.

Notes: Scalar radiomagnetotelluric measurements were carried out on a contaminated test area close to the Brazi Refinery in Romania in order to detect and to monitor a 1 m thick oil layer expected at 5 m depth. Radio transmitters broadcasting in a frequency range from 10 kHz to 300 kHz were selected to observe the apparent resistivity and the phase data associated with the E- and B-polarizations. They were located parallel and perpendicular to the assumed strike direction of the contamination plume. The data were interpreted by a 2D inversion technique from which the conductivity structure of the area was derived. The 2D inversion models of all profiles on the contaminated area show a poor-conductivity zone above the groundwater table which could be associated with the oil contamination.A first attempt was also made to monitor the contaminated layer: the radiomagnetotelluric measurements were repeated on the same profiles a year later, but this time in a dry period, not in a rainy one. The 2D inversion results of the measurements in the dry period indicate that the high-resistivity layer moved closer to the surface. Additional reference measurements were then carried out on a non-contaminated area situated at a distance from the refinery, in the opposite direction to the flow of the groundwater. These reference measurements were used for the derivation of the unperturbed geology and they were also compared with the measurements of the contaminated test area. There is a significant difference in the frequency dependences of the apparent resistivities of the reference and contaminated areas, which could indicate a contamination at shallow depth. The 2D inversion results show the increase of resistivity at a depth of about 5 m beneath the contaminated area where the oil contamination is expected according to the information from the boreholes.

Notes: Computer simulations are used to calculate the elastic properties of model cemented sandstones composed of two or more mineral phases. Two idealized models are considered – a grain-overlap clay/quartz mix and a pore-lining clay/quartz mix. Unlike experimental data, the numerical data exhibit little noise yet cover a wide range of quartz/cement ratios and porosities. The results of the computations are in good agreement with experimental data for clay-bearing consolidated sandstones.The effective modulus of solid mineral mixtures is found to be relatively insensitive to microstructural detail. It is shown that the Hashin–Shtrikman average is a good estimate for the modulus of the solid mineral mixtures. The distribution of the cement phase is found to have little effect on the computed modulus–porosity relationships. Numerical data for dry and saturated states confirm that Gassmann's equations remain valid for porous materials composed of multiple solid constituents. As noted previously, the Krief relationship successfully describes the porosity dependence of the dry shear modulus, and a recent empirical relationship provides a good estimate for the dry-rock Poisson's ratio.From the numerical computations, a new empirical model, which requires only a knowledge of system mineralogy, is proposed for the modulus–porosity relationship of isotropic dry or fluid-saturated porous materials composed of multiple solid constituents. Comparisons with experimental data for clean and shaly sandstones and computations for more complex, three-mineral (quartz/dolomite/clay) systems show good agreement with the proposed model over a very wide range of porosities.

Notes: With the increasing use of permanently installed seismic installations, many of the issues in time-lapse seismic caused by the lack of repeatability can be reduced. However, a number of parameters still influence the degree of reliability of 4D seismic data. In this paper, the specific impact of seawater velocity variations on time-lapse repeatability is investigated in a synthetic study. A zero-lag time-lapse seabed experiment with no change in the subsurface but with velocity changes in the water column is simulated. The velocity model in the water column is constant for the baseline survey while the model for the repeat survey is heterogeneous, designed from sea salinity and temperature measurements in the West of Shetlands. The difference section shows up to 80% of residual amplitude, which highlights the poor repeatability. A new dynamic correction which removes the effect of seawater velocity variations specifically for permanent installations is developed. When applied to the synthetic data, it reduces the difference residual amplitude to about 3%. This technique shows substantial improvement in repeatability beyond conventional time-lapse cross-equalization.

Notes: A new method for time-lapse signal separation and enhancement using singular-value decomposition is presented. Singular-value decomposition is used to separate a 4D signal into its constituent parts: common geology, time-lapse response and noise. Synthetic tests which demonstrate the advantages of the singular-value decomposition technique over traditional differencing methods are also presented. This signal separation and enhancement technique is used to map out both the original and moved oil–water contacts across the Nelson Field. The singular-value decomposition technique allows the oil–water contact to be mapped across regions which would have been missed using traditional differencing methods. In particular, areas toward the edges of the field are highlighted by the technique. The oil–water contact is observed to move upwards across the field, with the largest movements being associated, as anticipated, with natural production. The results obtained are broadly consistent with those predicted by the reservoir simulator model. Singular-value decomposition is demonstrated to be a useful tool for enhancing the time-lapse signal and for gaining confidence in areas where traditional differencing fails.

Notes: This work looks at the application of neural networks in geophysical well-logging problems and specifically their utilization for inversion of nuclear downhole data. Simulated neutron and γ-ray fluxes at a given detector location within a neutron logging tool were inverted to obtain formation properties such as porosity, salinity and oil/water saturation. To achieve this, the forward particle-radiation transport problem was first solved for different energy groups (47 neutron groups and 20 γ-ray groups) using the multigroup code EVENT. A neural network for each of the neutron and γ-ray energy groups was trained to re-produce the detector fluxes using the forward modelling results from 504 scenarios. The networks were subsequently tested on unseen data sets and the unseen input parameters (formation properties) were then predicted using a global search procedure. The results obtained are very encouraging with formation properties being predicted to within 10% average relative error. The examples presented show that neural networks can be applied successfully to nuclear well-logging problems. This enables the implementation of a fast inversion procedure, yielding quick and reliable values for unknown subsurface properties such as porosity, salinity and oil saturation.

Notes: Marine gravimeters mounted on stabilized platforms are commonly used in aircraft to perform airborne gravity measurements. The role of the stabilized platform is to level the sensor mechanically, whatever the aircraft attitude. However, this compensation is generally insufficient due to the sensitivity of modern gravity sensors. Correcting the offlevel error requires that an offlevel correction calculated from positioning data be added to gravimeter measurements, which complicates not only the processing, but also the assessment of precision and resolution.This paper is a feasibility study describing the levelling of a completely strapped-down LaCoste and Romberg gravimeter for airborne gravimetry operation, by means of GPS positioning data. It focuses on the calculation of the sensor offlevel correction needed for the complete gravity data processing. The precision of the offlevel correction that can be achieved, in terms of GPS data precision and gravity wavelengths, is theoretically studied and estimated using the gravity and GPS data acquired during the Alpine Swiss French airborne gravity survey carried out in 1998 over the French Western Alps. While a 1 cm precision of GPS-determined baseline coordinates is sufficient to achieve a 5 mGal precision of the offlevel correction, we maintain that this precision has to reach 1 mm to ensure a 1 mGal precision of the offlevel correction at any wavelength.Without a stabilized platform, the onboard instrumentation becomes significantly lighter. Furthermore, the correction for the offlevel error is straightforward and calculated only from GPS data. Thus, the precision and the resolution of airborne gravity surveys should be estimated with a better accuracy.

Notes: Dip-moveout (DMO) correction is often applied to common-offset sections of seismic data using a homogeneous isotropic medium assumption, which results in a fast execution. Velocity-residual DMO is developed to correct for the medium-treatment limitation of the fast DMO. For reasonable-sized velocity perturbations, the residual DMO operator is small, and thus is an efficient means of applying a conventional Kirchhoff approach. However, the shape of the residual DMO operator is complicated and may form caustics. We use the Fourier domain for the operator development part of the residual DMO, while performing the convolution with common-offset data in the space–time domain. Since the application is based on an integral (Kirchhoff) method, this residual DMO preserves all the flexibility features of an integral DMO. An application to synthetic and real data demonstrates effectiveness of the velocity-residual DMO in data processing and velocity analysis.

Notes: The finite-difference method applied to the full 3D wave equation is a rather time-consuming process. However, in the 2.5D case, we can take advantage of the medium symmetry. By taking the Fourier transform with respect to the out-of-plane direction (the symmetry axis), the 3D problem can be reduced to a repeated 2D problem. The third dimension is taken into account by a sum over the corresponding wave-vector component. A criterion for where to end this theoretically infinite sum derives from the stability conditions of the finite-difference schemes employed. In this way, the computation time of the finite-difference calculations can be considerably reduced. The quality of the modelling results obtained with this 2.5D finite-difference scheme is comparable to that obtained using a standard 3D finite-difference scheme.

Notes: For non-linear kinematic inversion of elastic anisotropy parameters and related investigations of the sensitivity of seismic data, the derivatives of the wavespeed (phase velocity and group velocity) with respect to the individual elastic moduli are required. This paper presents two analytic methods, called the eigenvalue and eigenvector methods, to compute the derivatives of the wavespeeds for wave propagation in a general anisotropic medium, which may be defined by up to 21 density-normalized elastic moduli. The first method employs a simple and compact form of the eigenvalue (phase velocity) and a general form of the group velocity, and directly yields general expressions of the derivatives for the three wave modes (qP, qS1, qS2). The second method applies simple eigenvector solutions of the three wave modes and leads to other general forms of the derivatives. These analytic formulae show that the derivatives are, in general, functions of the 21 elastic moduli as well as the wave propagation direction, and they reflect the sensitivity of the wavespeeds to the individual elastic moduli. Meanwhile, we give results of numerical investigations with some examples for particular simplified forms of anisotropy. They show that the eigenvalue method is suitable for the qP-, qS1- and qS2-wave computations and mitigates the singularity problem for the two quasi-shear waves. The eigenvector method is preferable to the eigenvalue method for the group velocity and the derivative of the phase velocity because it involves simpler expressions and independent computations, but for the derivative of the group velocity the derivative of the eigenvector is required. Both methods tackle the singularity problem and are applicable to any degree of seismic anisotropy for all three wave modes.

Notes: The phase error between the real phase shift and the Gazdag background phase shift, due to lateral velocity variations about a reference velocity, can be decomposed into axial and paraxial phase errors. The axial phase error depends only on velocity perturbations and hence can be completely removed by the split-step Fourier method. The paraxial phase error is a cross function of velocity perturbations and propagation angles. The cross function can be approximated with various differential operators by allowing the coefficients to vary with velocity perturbations and propagation angles. These variable-coefficient operators require finite-difference numerical implementation. Broadband constant-coefficient operators may provide an efficient alternative that approximates the cross function within the split-step framework and allows implementation using Fourier transforms alone. The resulting migration accuracy depends on the localization of the constant-coefficient operators. A simple broadband constant-coefficient operator has been designed and is tested with the SEG/EAEG salt model. Compared with the split-step Fourier method that applies to either weak-contrast media or at small propagation angles, this operator improves wavefield extrapolation for large to strong lateral heterogeneities, except within the weak-contrast region. Incorporating the split-step Fourier operator into a hybrid implementation can eliminate the poor performance of the broadband constant-coefficient operator in the weak-contrast region. This study may indicate a direction of improving the split-step Fourier method, with little loss of efficiency, while allowing it to remain faster than more precise methods such as the Fourier finite-difference method.

Notes: Multi-refractor imaging is a technique for constructing a single two-dimensional image of a number of refractors by stacking multiple convolved and cross-correlated reversed shot records. The method is most effective with high-fold data that have been obtained with roll-along acquisition programs because the stacking process significantly improves the signal-to-noise ratios. The major advantage of the multi-refractor imaging method is that all the data can be stacked to maximize the signal-to-noise ratios before the measurement of any traveltimes. However, the signal-to-noise ratios can be further increased if only those traces that have arrivals from the same refractor are used, and if the correct reciprocal times or traces are employed.A field case study shows that multi-refractor imaging can produce a cross-section similar to the familiar reflection cross-section with substantially higher signal-to-noise ratios for the equivalent interfaces.

Notes: Combined time-lapse reservoir simulation and seismic modelling has been performed on both 1D and 3D models of a channelized turbidite reservoir. The models have been built using core, log, laboratory and seismic data from the Nelson Field (central North Sea) as a template. Oil and water movement in the main channels, channel margins and interchannel regions is investigated, with a particular focus being the effect of poor net-to-gross. The analysis confirms that saturation effects dominate the response whilst stress-sensitivity effects play a minor role. The trough–peak signature in the seismic difference volumes formed by the sweep of the water can be continued and mapped slightly further than the channel margins. This characteristic 4D signature remains roughly intact, despite the complicated depositional architecture, and accurately delineates the area of moved fluid, but it needs additional calibration to combat the detrimental influence of the low net-to-gross. Signal strength is largely dependent on reservoir quality, but is also compromised by the net-to-gross, fluid distribution and, more critically, by the exact timing of the seismic survey. For example, a region of bypassed oil zone remains undetected as it forms early during the production. This work demonstrates that to understand fully the 4D signature at a quantitative level requires adequate knowledge of the fluid properties, but also, more critically, the geology.

Notes: In the central North Sea ‘Gannet-A’ field, a 50 ft oil rim is overlain by a gas cap of variable thickness. Oil is produced from horizontal wells which initially produced dry oil, but as the field became more mature, a significant water cut was seen in several wells. A dedicated 4D seismic monitor survey was acquired in order to assess the remaining distribution of oil reserves. By forward modelling the synthetic seismic response to parameters such as contact movement and residual saturations (using 2D and 3D wedge models), and comparing the results with real seismic data, we are able to decipher the contact movements across the field. It is shown that, in one part of the field, the increased water cut is caused primarily by the vertical displacement of the entire oil rim into the initial gas cap. This oil-rim displacement produces a very different 4D seismic response from the case of a static gas–oil contact and rising oil–water contact (normal production). As a result of these observations, we are able to optimize field production by both re-perforation of existing wells and by drilling sidetracks into the displaced rim: a brown-field development opportunity that might otherwise be missed.

Notes: The cross-calibration of different vintage data is an important prerequisite in attempting to determine the time-lapse seismic effects induced by hydrocarbon production in a reservoir. This paper reports the preprocessing and cross-calibration procedures adopted to modify the data of four seismic vintages (1982, 1989, 1992 and 1999) from the Oseberg field in the North Sea, for optimal conditions for a time-lapse seismic amplitude analysis. The final results, in terms of time-lapse variations, of acoustic impedance and of amplitude-versus-offset, are illustrated for selected data sets. The application of preprocessing to each individual vintage data set reduces the effects of the different acquisition and noise conditions, and leads to consistency in the amplitude response of the four vintages. This consistency facilitates the final amplitude cross-calibration that is carried out using, as reference, the Cretaceous horizon reflections above the Brent reservoir. Such cross-calibration can be considered as vintage-consistent residual amplitude correction.Acoustic impedance sections, intercept and gradient amplitude-versus-offset attributes and coherent amplitude-versus-offset estimates are computed on the final cross-calibrated data. The results, shown for three spatially coincident 2D lines selected from the 1982, 1989 and 1999 data sets, clearly indicate gas-cap expansion resulting from oil production. Such expansion is manifested as a decrease in acoustic impedance and a modification of the amplitude-versus-offset trends in the apical part of the reservoir.

Notes: This paper presents a new algorithm for estimating non-minimum-phase seismic wavelets by using the second- and higher-order statistics (HOS) of the wavelets. In contrast to many, if not most, of the HOS-based methods, the proposed method does not need to assume that subsurface seismic reflectivity is a non-Gaussian, statistically independent and identically distributed random process. The amplitude and phase spectra of the wavelets are estimated, respectively, using the second-order statistics (SOS) and third-order moment (TOM) of the wavelets, which will, in turn, be derived from the HOS of the seismic traces. In our approach, the wavelets can be ‘calculated’ from seismic traces efficiently; no optimization or inversion is necessarily required. Very good results have been obtained by applying this method to both synthetic and real-field data sets.

Notes: We study the propagation of elastic waves that are generated in a fluid-filled borehole surrounded by a cracked transversely isotropic medium. In the model studied the anisotropy and borehole axes coincide. To obtain the effective elastic moduli of a cracked medium we have applied Hudson's theory that enables the determination of the overall properties as a function of the crack orientation in relation to the symmetry axis of the anisotropic medium. This theory takes into account the hydrodynamic mechanism of the elastic-wave attenuation caused by fluid filtration from the cracks into a porous matrix.We have simulated the full waveforms generated by an impulse source of finite length placed on the borehole axis. The kinematic and dynamic parameters of the compressional, shear and Stoneley waves as functions of the matrix permeability, crack orientation and porosity were studied. The modelling results demonstrated the influence of the crack-system parameters (orientation and porosity) on the velocities and amplitudes of all wave types. The horizontally orientated cracks result in maximal decrease of the elastic-wave parameters (velocities and amplitudes).Based on the fact that the shear- and Stoneley-wave velocities in a transversely isotropic medium are determined by different shear moduli, we demonstrate the feasibility of the acoustic log to identify formations with close to horizontal crack orientations.

Notes: The development of an electromagnetic numerical modelling scheme for a magnetic dipole in an arbitrary casing segment in an inhomogeneous conductivity background has been difficult, due to the very high electrical conductivity and magnetic permeability contrasts between the steel casing and the background medium. To investigate the effect of steel casing efficiently, we have developed an accurate but simple finite-element modelling scheme to simulate electromagnetic fields in a medium of cylindrically symmetric conductivity structures. In order to preserve the cylindrical symmetry in the resulting electromagnetic fields, a horizontal loop current source is used throughout. One of the main advantages of the approach is that the problem is scalar when formulated using the azimuthal electric field, even if the casing is both electrically conductive and magnetically permeable. Field calculations have been made inside the cased borehole as well as in another borehole which is not cased. Careful analyses of the numerical modelling results indicate that the anomaly observed in a cross-borehole configuration is sensitive enough to be used for tomographic imaging.

Notes: We know much about the influence of management on stocks of organic matter in subtropical soils, yet little about the influence on the chemical composition. We therefore studied by CPMAS 13C NMR spectroscopy the composition of the above-ground plant tissue, of the organic matter of the whole soil and of silt- and clay-size fractions of the topsoil and subsoil of a subtropical Acrisol under grass and arable crops. Soil samples were collected from three no-till cropping systems (bare soil; oats−maize; pigeon pea + maize), each receiving 0 and 180 kg N ha−1 year−1, in a long-term field experiment. Soil under the original native grass was also sampled. The kind of arable crops and grass affected the composition of the particulate organic matter. There were no differences in the composition of the organic matter in silt- and clay-size fractions, or of the whole soil, among the arable systems. Changes were observed between land use: the soil of the grassland had larger alkyl and smaller aromatic C contents than did the arable soil. The small size fractions contain microbial products, and we think that the compositional difference in silt- and clay-size fractions between grassland and the arable land was induced by changes in the soil's microbial community and therefore in the quality of its biochemical products. The application of N did not affect the composition of the above-ground plant tissue nor of the particulate organic matter and silt-size fractions, but it did increase the alkyl C content in the clay-size fraction. In the subsoil, the silt-size fraction of all treatments contained large contents of aromatic C. Microscopic investigation confirmed that this derived from particles of charred material. The composition of organic matter in this soil is affected by land use, but not by variations in the arable crops grown.

Notes: Oxidative treatment can isolate a stable organic matter pool in soils for process studies of organic matter stabilization. Wet oxidation methods using hydrogen peroxide are widely used for that purpose, but are said to modify poorly crystalline soil constituents. We investigated the effect of a modified NaOCl oxidation (pH 8) on the mineral composition of 12 subsoils (4.9–38.2 g organic C kg−1) containing varying amounts of poorly crystalline mineral phases, i.e. 1.1–20.5 g oxalate-extractable Fe kg−1, and of different phyllosilicate mineralogy. Post-oxidative changes in mineral composition were estimated by (i) the determination of elements released into the NaOCl solution, (ii) the difference in dithionite- and oxalate-extractable Si, Al and Fe, and (iii) the specific surface areas (SSAs) of the soils. The NaOCl procedure reduced the organic C concentrations by 12–72%. The amounts of elements released into the NaOCl extracts were small (≤ 0.14 g kg−1 for Si, ≤ 0.13 g kg−1 for Al, and ≤ 0.03 g kg−1 for Fe). The SSA data and the amounts of dithionite- and oxalate-extractable elements suggest that the NaOCl oxidation at pH 8 does not attack pedogenic oxides and hydroxides and only slightly dissolves Al from the poorly crystalline minerals. Therefore, we recommend NaOCl oxidation at pH 8 for the purpose of isolating a stable organic matter pool in soils for process studies of organic matter stabilization.

Notes: A well-developed podzol hydrosequence that has been partially covered with drift sand, and partially subjected to improved drainage, provides new insights into the causes of variation in soil organic matter chemistry in such soils. While E horizons invariably move towards a dominance of aliphatic components reflecting residual accumulation, the chemistry of organic matter in well-drained B horizons is determined mainly by decaying roots, which are transformed by microorganisms to humus aggregates. In poorly drained, stratified B horizons, humus coatings dominate and the chemistry is very close to that of dissolved organic carbon. When a sand cover inhibits the supply of fresh litter, microbial decomposition in the A horizon causes a shift in chemistry towards that of the E horizon. Similarly, upon improved drainage and removal of complexed metals from the top of the B horizon, microbial decomposition of all palatable organic matter in the top of the B horizon causes a shift towards E-horizon chemistry. This is probably the mechanism by which most E horizons in podzols are formed, and not by re-solution. Marked chemical changes upon improved drainage may take only decades. During microbial decay, small polysaccharide-derived pyrolysis products (mainly furans, furaldehydes and acetic acid) remain abundant due to the contribution of microbial sugars. Both micromorphology and factor analysis on quantified results of pyrolysis-gas chromatography/mass spectrometry contribute significantly to the interpretation of the humus chemistry of these profiles and thus to our understanding of soil genesis. Organic chemistry of the investigated podzols can be understood only in the context of their genesis.

Notes: Water repellency can be a significant factor in soil physical behaviour, but little is known about the depth dependence of the contact angle of field soils. We investigated contact angles and wetting properties as a function of depth for a wide range of agricultural and forest soils in Germany. The agricultural soils ranged from silty to sandy texture (six profiles), and the forest soils ranged from sandy to loamy texture (eight profiles). Contact angles (CA) were measured with the Wilhelmy plate method (WPM). In most of the soils, advancing WPM contact angles were considerably greater than 0°  and they varied irregularly with depth. In general, sandy soils had larger WPM contact angles than silty soils. From the relation of the contact angle with texture and pH the quality of soil organic matter (SOM) was considered as more important for the wetting properties than the total amount of soil organic carbon (SOC). Finally, it was found that for soils with intermediate sand contents either under agricultural or forest use, the kind of land use seemed not to influence CA. Coarse-textured sandy soils that were used only as forest sites were more hydrophobic than silty soils which were exclusively used as agricultural soils. We conclude that a coarse texture favours, in combination with other factors (mainly pH), hydrophobic SOM.

Notes: Aggregate breakdown due to rainfall action causes crusting and interrill erosion. Erodibility is seemingly determined by the capacity of surface aggregates to resist the effects of rainfall. In this paper, we evaluated the relevance of an aggregate stability measurement, which comprises three treatments, in order to characterize aggregate breakdown dynamics. Two cultivated soils were studied: a clay loam slightly sensitive to erosion and a more susceptible silt loam. We compared the size distributions of microaggregates produced by the three aggregate stability treatments with the results from a rainfall simulation. The behaviour of four initial aggregate size classes (〈 3 mm, 3–5 mm, 5–10 mm and 10–20 mm) was also compared to study the influence of the initial aggregate size on the nature of resulting aggregates. The mean weight diameter was from 200 to 1400 µm for the silt loam and from 600 to 7000 µm for the clay loam. The two experiments – aggregate stability measurements and aggregate breakdown dynamics under rainfall – yielded similar results. Qualitative analysis showed that for both soils the sizes of fragments produced by breakdown with the aggregate stability tests and under rainfall were similar and seemed to be qualitatively independent of the size of initial aggregates. We first schematized the structural organization of aggregates in cultivated horizons with a simple hierarchical model at two levels: (i) 〈 250 µm microaggregates and (ii) 〉 250 µm macroaggregates made by the binding together of microaggregates. We then developed a model of aggregate breakdown dynamics under rainfall which gives, for various rainfall durations, the size distributions of resulting fragments on the basis of aggregate stability measurements. We obtained a correlation coefficient, r, of 0.87 for the silt loam and of 0.91 for the clay loam, showing that the experimental and predicted mass percentages were linearly related for each size fraction.

Notes: The use of ultrasonic energy for the dispersion of aggregates in studies of soil organic matter (SOM) fractionation entails a risk of redistribution of particulate organic matter (POM) to smaller particle-size fractions. As the mechanical strength of straw also decreases with increasing state of decomposition, it can be expected that not all POM will be redistributed to the same extent during such dispersion. Therefore, we studied the redistribution of POM during ultrasonic dispersion and fractionation as a function of (i) dispersion energy applied and (ii) its state of decomposition. Three soils were dispersed at different ultrasonic energies (750, 1500 and 2250 J g−1 soil) or with sodium carbonate and were fractionated by particle size. Fraction yields were compared with those obtained with a standard particle-size analysis. Undecomposed or incubated (for 2, 4 or 6 months) 13C-enriched wheat straw was added to the POM fraction (0.25–2 mm) of one of the soils before dispersion and fractionation. Dispersion with sodium carbonate resulted in the weakest dispersion and affected the chemical properties of the fractions obtained through its high pH and the introduction of carbonate. The mildest ultrasonic dispersion treatment (750 J g−1) did not result in adequate soil dispersion as too much clay was still recovered in the larger fractions. Ultrasonic dispersion at 1500 J g−1 soil obtained a nearly complete dispersion down to the clay level (0.002 mm), and it did not have a significant effect on the total amount of carbon and nitrogen in the POM fractions. The 2250 J g−1 treatment was too destructive for the POM fractions since it redistributed up to 31 and 37%, respectively, of the total amount of carbon and nitrogen in these POM fractions to smaller particle-size fractions. The amount of 13C-enriched wheat straw that was redistributed to smaller particle-size fractions during ultrasonic dispersion at 1500 J g−1 increased with increasing incubation time of this straw. Straw particles incubated for 6 months were completely transferred to smaller particle-size fractions. Therefore, ultrasonic dispersion resulted in fractionation of POM, leaving only the less decomposed particles in this fraction. The amounts of carbon and nitrogen transferred to the silt and clay fractions were, however, negligible compared with the total amounts of carbon and nitrogen in these fractions. It is concluded that ultrasonic dispersion seriously affects the amount and properties of POM fractions. However, it is still considered as an acceptable and appropriate method for the isolation and study of SOM associated with silt and clay fractions.

Notes: Soils are increasingly viewed as a potential sink for atmospheric carbon. However, their use to meet CO2 emission reductions is problematic for there are knowledge gaps regarding the mechanisms involved in the sequestration of organic carbon (OC). There is evidence which suggests that OC concentrations are controlled by the mineralogy and related specific surface area (SSA) of a given soil. The goal of this study was to examine the importance of sorptive mechanisms on OC stabilization. The objectives were (i) to determine the SSA, clay mineralogy and dithionite- and oxalate-extractable Fe and Al concentrations of several soils, and (ii) to analyse how these variables are related to OC concentrations. Five soils were sampled and analysed: two Umbrisols, a Stagnic Acrisol, an Anthrosol/Vertisol/Gleysol-Chernozem and a Gleysol (FAO terminology), all located in Hesse, Germany. Oxalate-extractable Fe and Al were found to be the best predictors of OC concentrations in the soils examined. Specific surface area correlated significantly with the OC content of the B and C horizons of one Umbrisol and the entire profile of the Anthrosol/Vertisol/Gleysol-Chernozem. The relationship between SSA and OC concentrations is likely to be restricted to certain soils and might be a product of the sorptive capacity of Fe and Al oxides. We can assume that the available mineral surface area on oxides is a limiting factor in terms of a soil's capacity to sequester organic carbon. As such, attention should be paid to soil mineralogy and how this might limit the use of soils as a sink for atmospheric CO2.

Notes: In the Ecuadorian Cordillera, the hardened volcanic ashes (cangahuas) account for 15% of the cultivated area. The soil resulting from the fragmentation of these materials, generally by heavy machinery, shows an apparent stable millimetric structure. However, this new structure is highly susceptible to disintegration under rain, because it contains no organic matter and little clay, and the material itself is readily eroded in consequence.We studied the evolution of soil aggregate stability in two factorial experiments during five cultivation cycles with two kinds of soil preparation and five fertilization treatments. The aggregate stability was not influenced by either kind of soil preparation, nor by large additions of cattle manure (80 t ha−1) or green manure (10 t ha−1), nor by growing a perennial grass. The variation in the aggregate stability seemed to depend on the components inherited from the original volcanic material: in the plots with larger clay content, and with swelling clay minerals, the aggregates were less stable than those composed of isometric fine silt particles.

Notes: The water retention characteristic provides the traditional data set for the derivation of a soil's pore-size distribution. However, the technique employed to achieve this requires that assumptions be made about the way pores interconnect. We explore an alternative approach based on stray field nuclear magnetic resonance (STRAFI-NMR) to probe the water-filled pores of both saturated and unsaturated soils, which does not require information relating to pore connectivity. We report the relative size distributions of water-occupied pores in saturated and unsaturated samples of two sets of glass beads of known particle size, two sands, and three soils (a silty loam, a sandy loam and a loamy sand), using measurements of the NMR T1 proton relaxation time of water. The T1 values are linearly related to pore size and consequently measured T1 distributions provide a measure of the pore-size distribution. For both the sands and the glass beads at saturation the T1 distributions are unimodal, and the samples with small particle sizes show a shift to small T1 values indicating smaller voids relative to the samples with larger particles. Different matric potentials were used to reveal how the water-occupied pore-size distribution changes during drainage. These changes are inconsistent with, and demonstrate the inadequacies of, the commonly employed parallel-capillary tube model of a soil pore space. We find that not all pores of the same size drain at the same matric potential. Further, we observe that the T1 distribution is shifted to smaller values beyond the distribution at saturation. This shift is explained by a change in the weighted average of the relaxation rates as the proportion of water in the centre of water-filled pores decreases. This is evidence for the presence of pendular structures resulting from incomplete drainage of pores. For the soils the results are similar except that at saturation the T1 distributions are bimodal or asymmetrical, indicative of inter-aggregate and intra-aggregate pore spaces. We conclude that the NMR method provides a characterization of the water-filled pore space which complements that derived from the water retention characteristic and which can provide insight into the way pore connectivity impacts on drainage.

Notes: Soil response to acid and sulphur inputs is influenced largely by the soil's physico-chemical properties. We studied the effects of such depositions in two types of Andosols exposed to volcanogenic emission (Masaya, Nicaragua), namely Eutric Andosols rich in allophanic constituents, and Vitric Andosols rich in volcanic glass. Small mineral reserves and large contents of secondary short-range ordered minerals indicate a more advanced weathering of the Eutric than the Vitric Andosols. Strong correlations between soil specific surface and oxalate-extractable Al, Si and Fe contents highlight the predominant contribution of short-range ordered minerals to surface area.Both types of Andosols showed a decrease in pH upon acid input. Sulphur deposition increased the soil's S content to 5470 mg S kg−1. However, the acid neutralizing capacity of the soil solid phase (ANCs) was not significantly affected by the acid and S inputs. Non-exchangeable (mineral reserve) and exchangeable cations and total contents of sulphur and phosphorus dictate most of the ANCs variation. In the Vitric Andosols, mineral reserves contributed up to 97% to these four additive pools, whereas the exchangeable cations accounted for 1–4%. In the Eutric Andosols, the contribution of mineral reserves was less (71–92%), but the exchangeable cation content was greater (1–20%), whereas the contribution of sulphur and phosphorus was significant at 1–15% and 2–7%, respectively. The main process involved in H+ consumption is mineral weathering in Vitric Andosols and ion exchange in Eutric Andosols.

Notes: Sulphate sorption on to the surface of short-range ordered minerals and precipitation of Al-hydroxy sulphate contribute to the acid neutralizing capacity of soils. The correct measurement of total inorganic sulphate is thus essential in soils that are accumulating SO42– anions. We extracted SO42– by various solutions, namely 0.005 m Ca(NO3)2, 0.016 m KH2PO4, 0.5 m NH4F and 0.2 m acidic NH4-oxalate (pH 3), from Vitric and Eutric Andosols exposed to prolonged deposition of acid and SO2 from an active volcano (Masaya, Nicaragua). We attributed sulphate extractable by KH2PO4 (20–3030 mg kg−1) to anion-exchangeable SO42–, which was much smaller than NH4F- and oxalate-extractable SO42– (400–9680 and 410–10 480 mg kg−1, respectively). Our results suggest the occurrence of a sparingly soluble Al-hydroxy-mineral phase extractable by both NH4F and oxalate. The formation of Al-hydroxy minerals would result from the combination of enhanced weathering caused by strong acid loading and simultaneous occurrence of large SO42– concentrations in soil solution. Oxalate extracted slightly more inorganic SO42– than did NH4F, this additional amount of SO42– correlating strongly with oxalate-extractable Si and Fe contents. Preferential occlusion of SO42– by short-range ordered minerals, especially ferrihydrite, explains this behaviour. If we exclude the contribution of occluded sulphate then oxalate and NH4F mobilize similar amounts of SO42– and are believed to mobilize all of the inorganic SO42– pool.

Notes: Topsoil samples from cultivated and adjacent non-cultivated fields on three major agricultural soils in North Cameroon were fractionated into particle-size fractions that were analysed subsequently for their C and 13C contents. The aim was to obtain further insight into the dynamics of soil organic matter (SOM) in relation to land use in Cameroon. Since organic carbon contents of the fractions were often very small, samples and analyses were extensively replicated to obtain robust statistical estimates of observed differences. For each soil type, differences in δ13C values between fields could be related to changes in the input and decomposition of organic matter arising from soil type, land management and, for example, the nature and abundance of weeds. Turnover of organic matter appeared to be fastest in the sand fraction, which is in line with results from earlier studies. In the finer fractions, clear differences in reaction to changes in input and decomposition were observed, that seem to be linked to differences in clay mineralogy. The results illustrate that SOM in the various fractions is much less stable and more strongly affected by changes in land use than might be assumed on the basis of changes in total SOM contents alone. At the same time, they demonstrate the relevance of 13C isotope analyses of SOM for studies on the impact of land use on these savannah soils with little SOM that are highly susceptible to degradation.

Notes: The effects of past applications of farmyard manure (FYM, applied from 1942 to 1967), metal-contaminated sewage sludge (applied from 1942 to 1961) and mineral fertilizer (NPK, applied from 1942 until now) on the microbial biomass and community structure in a sandy loam, arable soil from the Woburn Market Garden Experiment, UK, were investigated in 1998. Concentrations of Cu, Ni and Zn in soils which previously received sewage sludge were less than current European Union (EU) limits, but the soil Cd concentration was more than twice the permitted limit. Organic-C concentration in the FYM-treated soil and contaminated soils was about twice that of NPK-treated soil. The initial microbial biomass-C and estimates of total bacterial numbers by acridine orange direct count were significantly (P 〈 0.05) greater in the FYM-treated soil compared with the NPK-treated and the most contaminated soils. Total phospholipid fatty acid (PLFA) concentration (another measure of biomass) was significantly greater in the FYM-treated soil compared with either the low or high metal-contaminated soils, both of which contained similar PLFA concentrations. In the metal-contaminated soils, in contrast, fluorescent Pseudomonas counts, as a percentage of total plate counts, were at least 1.5 times greater than in the uncontaminated soils. The concentrations of these microbial parameters were significantly (P 〈 0.05) less in the NPK soil than in all the other treatments. Biomass-C as a percentage of organic-C was also significantly (P 〈 0.05) greater in the uncontaminated soils compared with the metal-contaminated soils. Biomass specific respiration rates in the metal-contaminated soils were c. 1.5 times those in the FYM-treated soil. In the metal-contaminated soils, the concentration of mono-unsaturated and hydroxy-fatty acids (derived from phospholipids), and lipopolysaccharide hydroxy-fatty acids (all indicative of Gram-negative bacteria) were significantly (P 〈 0.05) greater than branched fatty acids (indicative of Gram-positive bacteria). Furthermore, Gram-negative counts were 62–68% greater than Gram-positive counts in the metal-contaminated soils. Branched fatty acid concentration was significantly (P 〈 0.05) greater in the FYM-treated soil than in the metal-contaminated soils. Gram-positive counts were also 63% greater than Gram-negative counts in the FYM-treated soil. We found that effects of the relatively small heavy metal concentration caused measurable decreases in soil microbial biomass-C concentrations, acridine orange direct counts and Gram-positive counts. There were also increases in biomass specific respiration rates, and the microbial community had changed substantially, nearly 40 years after the metal inputs ceased. We conclude that, at the very least, the current EU permitted limits for heavy metals in agricultural soils should not be relaxed.

Notes: We used the wavelet transform to quantify the performance of models that predict the rate of emission of nitrous oxide (N2O) from soil. Emissions of N2O and other soil variables that influence emissions were measured on soil cores collected at 256 locations across arable land in Bedfordshire, England. Rate-limiting models of N2O emissions were constructed and fitted to the data by functional analysis. These models were then evaluated by wavelet variance and wavelet correlations, estimated from coefficients of the adapted maximal overlap discrete wavelet transform (AMODWT), of the fitted and measured emission rates.We estimated wavelet variances to assess whether the partition of the variance of modelled rates of N2O emission between scales reflected that of the data. Where the relative distribution of variance in the model is more skewed to coarser scales than is the case for the observation, for example, this indicates that the model predictions are too smooth spatially, and fail adequately to represent some of the variation at finer scales. Scale-dependent wavelet correlations between model and data were used to quantify the model performance at each scale, and in several cases to determine the scale at which the model description of the data broke down. We detected significant changes in correlation between modelled and predicted emissions at each spatial scale, showing that, at some scales, model performance was not uniform in space. This suggested that the influence of a soil variable on N2O emissions, important in one region but not in another, had been omitted from the model or modelled poorly. Change points usually occurred at field boundaries or where soil textural class changed.We show that wavelet analysis can be used to quantify aspects of model performance that other methods cannot. By evaluating model behaviour at several scales and positions wavelet analysis helps us to determine whether a model is suitable for a particular purpose.

Notes: Crack development is predominant in soil structure formation. A number of fracture mechanics models have been applied to soil to describe cracking, but most are not applicable for soil in a wet, plastic state. We address this weakness by applying a new elastic–plastic fracture mechanics approach to describe crack formation in plastic soil. Samples are fractured using a deep-notch (modified four-point) bend test, with data on load transmission, sample bending, crack growth, and crack-mouth opening collected to assess the crack-tip opening angle (CTOA). CTOA provides a powerful parameter for describing soil cracking since it can be induced by soil shrinkage (an easily measured parameter) and can be used to describe elastic–plastic fracture in numerical approximations, such as finite element modelling. The test variables we studied were the direction of the applied consolidation stress, clay content, and pore water salinity. All samples were formed by consolidating soil slurry one-dimensionally with a 120-kPa vertical effective stress. Tests on pure kaolinite showed that the direction of the consolidation stress did not affect CTOA, which was 0.23 ± 0.02 m m−1 for specimens cut both in a horizontal and in a vertical direction to the applied stress. Soil clay content had a marked influence, however, with silica sand:kaolinite mixtures by weight of 20:80 and 40:60 reducing CTOA to 0.14 ± 0.02 m m−1 and 0.12 ± 0.01 m m−1, respectively. These smaller values of CTOA indicate that less strain is required to induce fracture when the amount of clay is less. Salinity (0.5 m NaCl) caused a reduction in the CTOA of pure kaolinite from 0.23 ± 0.02 m m−1 to 0.17 ± 0.03 m m−1.

Notes: In soil carbon dynamics, the role of physicochemical interactions between organic matter and minerals is not well understood nor quantified. This paper examines the interactions between soil organic matter and poorly crystalline aluminosilicates in a volcanic ash soil on La Réunion in the southern tropics. The soil examined is a profile composed of a surface soil (L-Ao-E-Bh) overlying four buried horizons (horizons 2Bw, 3Bw, 4Bw, 5Bw) that have all developed from successive tephra deposits. Non-destructive spectroscopy (XRD, FTIR and NMR of Si and Al) showed that the mineralogical composition varies from one buried horizon to another. Further, we show that buried horizons characterized by large amounts of crystalline minerals (feldspars, gibbsite) have the least capacity to store organic matter and the fastest carbon turnover. In contrast, buried horizons containing much poorly crystalline material (proto-imogolite and proto-imogolite allophane, denoted LP-ITM) store large amounts of organic matter which turns over very slowly. To understand the mechanism of interactions between LP-ITM and organic matter better, we focused on a horizon formed exclusively of LP-ITM. We demonstrate, using Δ14C and δ13C values, that even though LP-ITM is extraordinarily effective at stabilizing organic matter, C linked to LP-ITM is still in dynamic equilibrium with its environment and cycles slowly. Based on Δ14C values, we estimated the residence time of organic C as ∼ 163 000 years for the most stabilized subhorizon, a value that is comparable to that for organic carbon stabilized in Hawaiian volcanic soils. However, this calculation is likely to be biased by the presence of microcharcoal. We characterized the organo-mineral binding between organic matter and LP-ITM by 27Al NMR, and found that the organic matter is not only chelated to LP-ITM, but it may also limit the polymerization of mineral phases to a stage between proto-imogolite and proto-imogolite allophane. Our results demonstrate the important role of poorly crystalline minerals in the storage of organic C, and show that mineral and organic compounds have to be studied simultaneously to understand the dynamics of organic C in the soil.

Notes: The most southerly ombrotrophic peat bogs in Europe are in Galicia (northwest Spain). The humic matter in them originates from chemical processes in anaerobic conditions. We studied the acid properties of fulvic acids and humic acids isolated from two peat horizons of an ombrotrophic peat bog by potentiometric titration. Solutions containing 25, 50 and 100 mg l−1 of each humic substance were titrated at ionic strengths 0.005 m, 0.01 m and 0.1 m (with KNO3 as the inert electrolyte). Charge curves were analysed with a Donnan model to determine the intrinsic proton binding parameters. The concentration of the humic substance affected the charge curves more significantly at pH exceeding 6, and tended to disappear at greater concentrations. The proton binding conditional constants decreased with increasing ionic strength, this effect being more significant in the carboxylic groups with less affinity for protons. The proton binding constant of the carboxyl groups in a fulvic acid was one order of magnitude less than the value for the corresponding humic acid, whereas for the phenolic groups the values for both fractions were similar. The total content of acid groups was approximately 2 mol kg−1 greater in the fulvic fraction than in the humic fraction. Both humic fractions from the lower horizon contained more acid groups than those from the upper horizon, mainly because the content of carboxyl groups increases with soil depth. Therefore, the humic substances in the lower horizon of the peat will be more negatively charged, which will affect their solubility and the binding of metal ions.

Notes: Risk assessment of groundwater pollution requires quantitative information on the release kinetics of pollutants and organic matter from contaminated soil. We applied a new experimental design for column outflow experiments to investigate the release of polycyclic aromatic hydrocarbons (PAHs) and dissolved organic matter under water-saturated flow conditions. We used materials originating from a soil contaminated with non-aqueous phase liquids. To distinguish between release at equilibrium and release limited by mass transfer, we used two flow velocities and multiple interruptions to the flow. We quantified release and transport parameters by inverse numerical simulation of the individual breakthrough curves, applying a model based on the advection–dispersion equation including non-equilibrium and non-linear sorption. Release of the dissolved organic C takes place in two steps. Initially, a large amount of readily available organic matter is released. This first flush is followed by an outflow with typical characteristics of rate-limited release: larger concentrations in slower flow and increased concentrations after interruptions. The breakthrough of PAHs responds neither to the different flow velocities nor to the interruptions. We hypothesize that release of PAHs from the contaminated material is governed by dissolution at equilibrium according to Raoult's law. The boundary conditions of the experimental design, i.e. the flow velocities and multiple interruptions, enable us to distinguish between release at equilibrium and that which is rate-limited. Also, the response of the breakthrough behaviour to the boundary conditions can be used to estimate inversely effective release parameters.

Notes: Partitioning carbon (C) pools presents theoretical and practical challenges in soil modelling. In most current models, organic matter (OM) is commonly divided into a number of pools using data-fitting techniques. It is often not possible to measure the properties of each pool of soil C simulated in this way. Furthermore, there is no effective way in current soil models to properly simulate inputs of OM with variable qualities in terms of decomposability. Here, we describe a simple model, taking the OM as a whole and assuming a rate constant function changing exponentially with time, to simulate C decomposition in a multi-component OM pool. The model requires fewer parameters to be estimated than current multi-component models, and the simulated properties of OM are measurable. When changing the quantity and the quality of OM input, the model produces results similar to multi-component models, but avoids the difficulty and uncertainty of OM partitioning.

Notes: The estimation of soil carbon content is of pressing concern for soil protection and in mitigation strategies for global warming. This paper describes the methodology developed and the results obtained in a study aimed at estimating organic carbon contents (%) in topsoils across Europe. The information presented in map form provides policy-makers with estimates of current topsoil organic carbon contents for developing strategies for soil protection at regional level. Such baseline data are also of importance in global change modelling and may be used to estimate regional differences in soil organic carbon (SOC) stocks and projected changes therein, as required for example under the Kyoto Protocol to the United Nations Framework Convention on Climate Change, after having taken into account regional differences in bulk density.The study uses a novel approach combining a rule-based system with detailed thematic spatial data layers to arrive at a much-improved result over either method, using advanced methods for spatial data processing. The rule-based system is provided by the pedo-transfer rules, which were developed for use with the European Soil Database. The strong effects of vegetation and land use on SOC have been taken into account in the calculations, and the influence of temperature on organic carbon contents has been considered in the form of a heuristic function. Processing of all thematic data was performed on harmonized spatial data layers in raster format with a 1 km × 1 km grid spacing. This resolution is regarded as appropriate for planning effective soil protection measures at the European level. The approach is thought to be transferable to other regions of the world that are facing similar questions, provided adequate data are available for these regions. However, there will always be an element of uncertainty in estimating or determining the spatial distribution of organic carbon contents of soils.

Notes: Stones on the surface of the soil enhance infiltration and protect the soil against erosion. They are often removed in modern mechanized agriculture, with unfortunate side-effects. We evaluated experimentally the influence of surface stones on infiltration, runoff and erosion under field conditions using a portable rainfall simulator on bare natural soil in semi-arid tropical India, because modernization and mechanization often lead to removal of these stones in this region. Four fields with varied cover of stones from 3 to 65% were exposed to three rainfall intensities (48.5, 89.2 and 136.8 mm hour−1). Surface stones retarded surface runoff, increased final infiltration rates, and diminished sediment concentration and soil loss. The final infiltration ranged from 26 to 83% of rainfall when the rainfall intensity was 136.8 mm hour−1. The reduction in runoff and soil erosion and increase in infiltration were more pronounced where stones rested on the soil surface than where they were buried in the surface layer. The sediment yield increased from 2 g l−1 for 64.7% stone cover with rainfall of 48.5 mm hour−1 to 70 g l−1 for 3.5% stone cover with rain falling at 136.8 mm hour−1. The soil loss rate was less than 2 t ha−1 hour−1 for the field with stone cover of 64.7% even when the rainfall intensity was increased to 136.8 mm hour−1. The effects of stones on soil loss under the varied rainfall intensities were expressed mathematically. The particles in the sediment that ran off were mostly of silt size.

Notes: Phosphorus loss from land can be a major factor affecting surface water quality. We studied P-release mechanisms in wetland soils that had been drained and cultivated for four decades and then re-flooded. We measured redox, pH and solution composition in two sites in the field and in four peat and calcareous soils incubated in biogeochemical microcosms. The redox and pH measurements during the 120 days of incubation and the resulting soil solution composition indicated that the main process leading to P release is reductive dissolution of ferric hydroxides on which P was adsorbed and in which P was occluded. The molar Fe:P ratio increased with period of reduction from below 1 in the first week of re-flooding to 15–60 after 120 days. This suggests an increased P-retention capacity upon reoxidation of the soil solution, whether within the soil profile or in the drainage canals. Prolonged flooding of the calcite-poor, gypsum-rich peat soils increased the oversaturation of soil solutions with respect to hydroxyapatite and occasionally β-Ca3(PO4)2(c), indicating that in spite of the large Ca concentration, the rate of Ca-P precipitation was insufficient to maintain the saturation status of the Ca-P system. In the calcareous soils the Ca-P system effectively controlled the P activity in soil solution throughout the incubation period. In both cases the precipitation of Ca-P minerals could be an important P-retention mechanism.

Notes: Stable microaggregates can physically protect occluded soil organic matter (SOM) against decomposition. We studied the effects of agricultural management on the amount and characteristics of microaggregates and on SOM distribution in a marine loam soil in the Netherlands. Three long-term farming systems were compared: a permanent pasture, a conventional-arable system and an organic-arable system. Whole soil samples were separated into microaggregates (53–250 µm), 20–53 µm and 〈 20 µm organo-mineral fractions, sand and particulate organic matter, after complete disruption of macroaggregates. Equal amounts of microaggregates were isolated, irrespective of management. However, microaggregates from the pasture contained a larger fraction of total soil organic C and were more stable than microaggregates from the two arable fields, suggesting greater SOM stabilization in microaggregates under pasture. Moreover, differences in the relative contribution of coarse silt (〉 20 µm) versus fine mineral particles in the microaggregates of the different management systems demonstrate that different types of microaggregates were isolated. These results, in combination with micromorphological study of thin sections, indicate that the great earthworm activity under permanent pasture is an important factor explaining the presence of very stable microaggregates that are relatively enriched in organic C and fine mineral particles. Despite a distinctly greater total SOM content and earthworm activity in the organic- versus the conventional-arable system, differences in microaggregate characteristics between both arable systems were small. The formation of stable and strongly organic C-enriched microaggregates seems much less effective under arable conditions than under pasture. This might be related to differences in earthworm species' composition, SOM characteristics and/or mechanical disturbance between pasture and arable land.

Notes: Design and analysis of land-use management scenarios requires detailed soil data. When such data are needed on a large scale, pedotransfer functions (PTFs) could be used to estimate different soil properties. Because existing regression-based PTFs for estimating cation exchange capacity (CEC) do not, in general, apply well to arid areas, this study was conducted (i) to evaluate the existing models and (ii) to develop neural network-based PTFs for predicting CEC in Aridisols of Isfahan in central Iran. As most researches have found a significant correlation between CEC and soil organic matter content (OM) and clay content, we also used these two variables for modelling of CEC. We tested several published PTFs and developed two neural network algorithms using multilayer perceptron and general regression neural networks based on a set of 170 soil samples. The data set was divided into two subsets for calibration and testing of the models. In general, the neural network-based models provided more reliable predictions than the regression-based PTFs.

Notes: The distribution of iron-cyanide complexes between ferrocyanide, [FeII(CN)6]4–, and ferricyanide, [FeIII(CN)6]3–, in soils on contaminated sites depends on the redox potential, EH. We carried out microcosm experiments in which ferrocyanide (20 mg l−1) was added to an uncontaminated moderately acidic subsoil (pH 5.2), and varied the EH of the soil suspension between 200 and 700 mV over up to 109 days. Ferrocyanide and ferricyanide were analysed by capillary isotachophoresis. At redox potentials ranging from 400 to 700 mV, small amounts of iron-cyanide complexes were adsorbed, and ferrocyanide was almost completely oxidized to ferricyanide. Decreasing EH to 200 mV led to nearly complete removal of iron-cyanide complexes from solution, and the complexes were not mobilized after subsequent aeration (EH 〉 350 mV). Under weakly to moderately reducing conditions (EH ≈ 200 mV), iron-cyanide complexes were removed from solution by precipitation, which occurred, presumably in the form of e.g. Fe2[FeII(CN)6], Fe4[FeII(CN)6]3 or Mn2[FeII(CN)6], after reductive dissolution of Mn and Fe oxides. Four different sets of geochemical model calculations were carried out. The species distribution between ferrocyanide and ferricyanide in solution was predicted reliably under varying pH and redox conditions when iron-cyanide complex concentrations and Fe concentrations, excluding Fe bound in iron-cyanide complexes, were used in model calculations. In model calculations on the fate of iron-cyanide complexes in soil, adsorption reactions must be considered, especially under oxidizing conditions. Otherwise, the calculated iron-cyanide complex concentrations are larger than those actually measured.

Notes: Water movement in soil can be described accurately at the local scale, provided that soil hydraulic properties can be determined with precision. Traditional methods for characterizing soil are often time consuming, and large areas cannot be sampled easily. We present a simple method for overcoming these difficulties. It is easy to implement and cheap. It is known as the Beerkan method, and it relies on particle-size analysis, dry bulk density and simple infiltration tests in cylinders. We describe the experimental protocol and the method of data analysis, leading to the estimation of parameters describing hydraulic properties. Shape parameters depend on soil texture and are derived from particle-size data. Normalization parameters depend on soil structure. They are derived by inverse modelling and optimization from the infiltration tests. The theoretical background relies on the sorptivity concept and scaled forms of the infiltration equation. The formalism for one- and three-dimensional analysis is described. We assess the accuracy of the method using published data and simulated values, showing the soundness of the approach. For the purpose of illustration, we implemented a simple optimization technique on two bounding cases.

Notes: Fires in boreal forests frequently convert organic matter in the organic layer to black carbon, but we know little of how changing fire frequency alters the amount, composition and distribution of black carbon and organic matter within soils, or affects podzolization. We compared black carbon and organic matter (organic carbon and nitrogen) in soils of three Siberian Scots pine forests with frequent, moderately frequent and infrequent fires.Black carbon did not significantly contribute to the storage of organic matter, most likely because it is consumed by intense fires. We found 99% of black carbon in the organic layer; maximum stocks were 72 g m−2. Less intense fires consumed only parts of the organic layer and converted some organic matter to black carbon (〉 5 g m−2), whereas more intense fires consumed almost the entire organic layer. In the upper 0.25 m of the mineral soil, black carbon stocks were 0.1 g m−2 in the infrequent fire regime.After fire, organic carbon and nitrogen in the organic layer accumulated with an estimated rate of 14.4 g C m−2 year−1 or 0.241 g N m−2 year−1. Maximum stocks 140 years after fire were 2190 g organic C m−2 and 40 g N m−2, with no differences among fire regimes. With increasing fire frequency, stocks of organic carbon increased from 600 to 1100 g m−2 (0–0.25 m). Stocks of nitrogen in the mineral soil were similar among the regimes (0.04 g m−2). We found that greater intensities of fire reduce amounts of organic matter in the organic layer but that the greater frequencies may slightly increase amounts in the mineral soil.

Notes: The arrangement of a plant's roots in the soil determines the ability of the plant to resist uprooting. We have investigated the influence of root morphology on anchorage using simple patterns of root systems and numerical simulation. The form and mechanical properties of roots were derived from results found in the literature. Major parameters determining soil characteristics, root patterns and strength were varied so that their influence could be evaluated. The design of the experimental method we used generated an optimal number of configurations of different root systems, the tensile resistances of which were calculated by two-dimensional finite element analysis. We could quantify the influence of specific parameters, e.g. branching angle, number of lateral roots and soil cohesion, as well as global parameters such as total contact area, basal diameter and volume of the whole root system. We found that the number of roots and the diameter of roots were major components affecting the resistance to uprooting. The combination of topology and biomass explained 70% of the variation of tensile resistance.

Notes: The inability of physical and chemical techniques to separate soil organic matter into fractions that have distinct turnover rates has hampered our understanding of carbon (C) and nutrient dynamics in soil. A series of soil organic matter fractionation techniques (chemical and physical) were evaluated for their ability to distinguish a potentially labile C pool, that is ‘recent’ root and root-derived soil C. ‘Recent’ root and root-derived C was operationally defined as root and soil C labelled by 14CO2 pulse labelling of rye grass–clover pasture growing on undisturbed cores of soil. Most (50–94%) of total soil + root 14C activity was recovered in roots.Sequential extraction of the soil + roots with resin, 0.1 m NaOH and 1 m NaOH allocated ‘recent’ soil + root 14C to all fractions including the alkali-insoluble residual fraction. Approximately 50% was measured in the alkali-insoluble residue but specific activity was greater in the resin and 1 m NaOH fractions. Hot 0.5 m H2SO4 hydrolysed 80% of the 14C in the alkali-insoluble residue of soil + roots but this diminished specific activity by recovering much non-14C organic matter. Pre-alkali extraction treatment with 30% H2O2 and post-alkali treatment extractions with hot 1 m HNO3 removed organic matter with a large 14C specific activity from the alkali-insoluble residue.Density separation failed to isolate a significant pool of ‘recent’ root-derived 14C. The density separation of 14C-labelled roots, and roots remixed with non-radioactive soil, showed that the adhesion of soil particles to young 14C-labelled roots was the likely cause of the greater proportion of 14C in the heavy fraction.Simple chemical or density fractionations of C appear unsuitable for characterizing ‘recent’ root-derived C into fractions that can be designated labile C (short turnover time).

Notes: The physical characteristics of the soil at the root–soil interface are crucial because they determine both physical aspects of root function such as water and nutrient uptake and the microbial activity that is most relevant to root growth. Because of this we have studied how root activity modifies the structure and water retention characteristic of soil adjacent to the root for maize, wheat and barley. These plants were grown in pots for a 6-week growth period, then the soil adjacent to the root (rhizosphere soil) and bulk soil aggregates were harvested. These soil aggregates were then saturated and equilibrated at matric potentials between −600 kPa and saturation, and the water retention characteristics were measured. From subsamples of these aggregates, thin sections were made and the porosity and pore-size distributions were studied with image analysis. Both image analysis and estimates of aggregated density showed that the rhizosphere soil and bulk soil had similar porosities. Growing different plants had a small but significant effect on the porosity of the soil aggregates. Image analysis showed that for all the plant species the structure of the rhizosphere soil was different to that of the bulk soil. The rhizosphere soil contained more larger pores. For maize and barley, water retention characteristics indicated that the rhizosphere soil tended to be drier at a given matric potential than bulk soil. This effect was particularly marked at greater matric potentials. The difference between the water retention characteristics of the bulk and rhizosphere soil for wheat was small. We compare the water retention characteristics with the data on pore-size distribution from image analysis. We suggest that differences in wetting angle and pore connectivity might partly explain the differences in water retention characteristic that we observed. The impact of differences between the water retention properties of the rhizosphere and bulk soil is discussed in terms of the likely impact on root growth.

Notes: In highly weathered tropical conditions, soil organic matter is important for soil quality and productivity. We evaluated the effects of deforestation and subsequent arable cropping on the qualitative and quantitative transformation of the humic pool of the soil at three locations in Nigeria. Cultivation reduced the humic pool in the order: acetone-soluble hydrophobic fraction (HE) 〉 humic acid (HA) 〉 humin (HU) 〉 fulvic acid (FA), but not to the same degree at all three sites. The C and N contents, as well as the C/N ratios of humic extracts, were large and not substantially influenced by land use. The δ13C values of the humic extracts were invariably more negative in forested soils thereby showing a dilution of δ13C signature with cultivation from C3 to C4 plants. The δ13C values of apolar HE fractions were generally more negative, indicating a reduced sensitivity compared with other humic fractions to turnover of crop residues. The contents of hydrophobic constituents (alkyl and aromatic C), as revealed by cross-polarization magic angle spinning (CPMAS) 13C-NMR spectroscopy, in HA, FA and HU were generally 〈 50%, with the exception of larger hydrophobicity in HU in the forested soil at Nsukka and HA in that at Umudike. The HE fraction contained significantly more apolar constituents, and consequently had a larger intrinsic hydrophobicity than the other humic fractions. The larger reduction of apolar humic constituents than of the less hydrophobic humic fractions, when these soils were deforested for cultivation, indicates that at those sites the stability of accumulated organic matter is to be ascribed mainly to the selective preservation of hydrophobic compounds.

Notes: The value of nested sampling for exploring the spatial structure of univariate variation of the soil has been demonstrated in several studies and applied to practical problems. This paper shows how the method can be extended to the multivariate case. While the extension is simple in theory, in practice the direct estimation of covariance components by equating mean-square matrices with their expectation will often lead to estimates that are not positive semidefinite. This paper discusses solutions to this problem for balanced and unbalanced sample designs. In the balanced case there is a residual maximum likelihood (REML) estimator that will find estimates of covariance components that maximize an overall likelihood on the condition that all components are positive semidefinite (p.s.d.). This is possible because the condition is met if the differences of successive mean-square matrices are positive semidefinite, and this constraint can be incorporated into an algorithm. This does not hold for unbalanced designs. In this paper the problem was solved for unbalanced designs by scaling covariance components that were not p.s.d. to the nearest p.s.d. matrix according to a Euclidean distance.These methods were applied to data from three surveys, two with balanced and one with unbalanced sampling. Different patterns of scale-dependence of the correlation of soil properties were found. For example, at Ginninderra Experimental Station in Australia the soil water content and bulk density were correlated significantly, with the correlation increasing with distance to 56 m, but at longer distances the properties were not significantly correlated. By contrast, the pH of the soil and the available P content showed correlation that increased with distance. The implications of these results for planning more detailed sampling, both for prediction and for investigation of processes, are discussed.

Notes: Soils and sediments contain only small amounts of organic matter, and large concentrations of paramagnetic metals can give poor solid-state nuclear magnetic resonance (NMR) spectra of organic matter. Pretreatment of samples with hydrofluoric acid (HF) dissolves significant proportions of the mineral matrix and extracts paramagnetic elements. We investigated the effects of 10% HF treatment on the stable isotope content of carbon (C) and nitrogen (N) of organic matter from soils, composts and shales. Additionally we inferred molecular and isotopic characteristics of lost materials from calculations of isotope mass balances. Treatment with HF enriched C and N in mineral samples substantially (factors 2.5–42.4), except for Podzol B horizons (1.1–1.7) and organic material (1.0–1.3). After treatment most of the C (59.7–91.7%) and N (53.7–86.6%) was recovered, although changing C/N ratios often indicated a preferential loss of N-rich material. Isotope ratios of C and N in the remaining material became more negative when net alterations exceeded 0.3‰. The isotope ratios of the lost material contained more 13C (1–2‰) and 15N (1–4‰) than the initial organic matter. Acid hydrolysis typically removes proteins, amino acids and polysaccharides, all of which are enriched in 13C, and in the case of proteins and amino acids, enriched in 15N as well. We conclude that HF treatment released fresh, soluble, probably microbial, biomass in addition to carbohydrates. Net changes of the bulk chemical composition of organic matter were small for most soils, size fractions and plant material, but not for samples containing little organic matter, or those rich in easily soluble organic matter associated with iron oxides, such as Podzol B horizons.

Notes: Aromatic acids can reach the soil from direct anthropogenic activities or, indirectly, from the degradation of many aromatic compounds, such as pesticides or polycyclic aromatic hydrocarbons. Because of the anionic character of aromatic acids at the pH of most soil and sediment environments, they are expected to move rapidly through the soil profile and to pose a great risk of ground water contamination. We designed batch and column leaching tests to characterize the behaviour of three aromatic acids differing in their chemical structures, picloram (4-amino-3,5,6-trichloropicolinic acid), phthalic acid (2,2-benzenedicarboxylic acid), and salicylic acid (2-hydroxybenzoic acid), in four European soils with different physicochemical characteristics. Batch experiments revealed that the persistence of the three acids in soil:water suspensions decreased in the order: picloram ≫ phthalic acid 〉 salicylic acid, and their dissipation curves were relatively independent of soil type. Sorption by the soils, their clay-size fractions and model sorbents indicated much greater affinity of soil constituents for salicylic acid than for picloram or phthalic acid, most likely due to the ability of salicylic acid to form bidentate complexes with positively charged soil components. The extent of leaching of the aromatic acids in hand-packed soil columns decreased in the order: picloram (90–96%) 〉 phthalic acid (25–90%) 〉 salicylic acid (0–37%), which was consistent with the sorption and persistence results of the batch tests. The organic C content, the amount of small-size pores, and the initial concentration of aromatic acid in soil appeared to be important factors influencing the leaching patterns of phthalic acid and salicylic acid in the soils studied, but did not greatly influence the leaching pattern of picloram. Sorption and leaching of polar aromatic acids in soil can therefore vary considerably depending on the structural characteristics of the aromatic acid or soil type.

Notes: Preferential flow, as it bypasses the soil matrix, can greatly enhance the leaching of chemicals. When a soil is drained there is the risk that such short-circuiting results in more or less direct passage of polluting chemicals from the soil to the groundwater. If the groundwater table is shallow the chemicals could be transferred back into the surface soil by hydraulic lift through roots and subsequent release by exudation or from decaying plant residues and again become exposed to leaching by preferential flow, thus strongly enhancing the chance of export via the drains.We investigated the leaching of bromide in a tile-drained arable field over 2 years of crop rotation. The site was a former wetland, artificially drained a century ago for agriculture. Bromide was applied over 1.6 ha at a dosage of 10 g Br per m2 in August 1995 after the harvest of wheat. During the 2 years 18% of the applied bromide was exported via the drainage system, most of it in preferential flow events and more than half of it in a single winter storm 5 months after the application. Within 7 months 56% of the applied tracer was leached out of the main root zone into the groundwater. Subsequently the tracer re-emerged in water taken up by sugar beet in the following season. The beet accumulated 50% of the initially applied bromide in their leaves and released it again after harvest when the leaves were left as green manure on the field. Our results show that this recycling of solutes to the topsoil can have an important influence on their leaching as the solutes are thus again exposed to preferential transport into drains in the course of preferential flow events.

Notes: We tested the Rothamsted Carbon Model (RothC) against three long-term (27–28 years) experimental sites on Thai upland soils in order to see how this widely used ‘temperate’ soil carbon turnover model performed in a typical farming region in the tropics. We were able to verify – over a much longer period than had been examined in previous studies – that RothC performs well in a tropical region in plots used for continuous cropping experiments of maize and cassava without organic matter application. However, the model overestimated soil organic carbon (SOC) in some plots to which large amounts of organic matter (rice straw or cassava stalks) were applied. This overestimate could not be attributed to errors in estimating either the amount of C input to the soil or the ratio of decomposable plant materials to resistant plant materials entering the soil. Among many factors affecting SOC dynamics (e.g. weather conditions, soil characteristics, etc.), which are different in tropical regions from temperate regions, we conclude that the activity of soil fauna might be a major factor which makes the performance of RothC worse where much organic matter was applied. We suggest that care should be taken when applying RothC to tropical soils with large amounts of added organic matter.

Notes: Intact lysimeters (50 cm diameter, 70 cm deep) of silt loam soil under permanent grassland were used to investigate preferential transport of phosphorus (P) by leaching immediately after application of dairy effluent. Four treatments that received mineral P fertilizer alone (superphosphate at 45 kg P ha−1 year−1) or in combination with effluent (at ∼ 40–80 kg P ha−1 year−1) over 2 years were monitored. Losses of total P from the combined P fertilizer and effluent treatments were 1.6–2.3 kg ha−1 (60% of overall loss) during eight drainage events following effluent application. The rest of the P lost (40% of overall loss) occurred during 43 drainage events following a significant rainfall or irrigation compared with 0.30 kg ha−1 from mineral P fertilizer alone. Reactive forms of P (mainly dissolved reactive P: 38–76%) were the dominant fractions in effluent compared with unreactive P forms (mainly particulate unreactive P: 15–56%). In contrast, in leachate following effluent application, particulate unreactive P was the major fraction (71–79%) compared with dissolved reactive P (1–7%). The results were corroborated by 31P nuclear magnetic resonance analysis, which showed that inorganic orthophosphate was the predominant P fraction present in the effluent (86%), while orthophosphate monoesters and diesters together comprised up to 88% of P in leachate. This shows that unreactive P forms were selectively transported through soil because of their greater mobility as monoesters (labile monoester P and inositol hexakisphosphate) and diesters. The short-term strategies for reducing loss of P after application of dairy effluent application should involve increasing the residence time of applied effluent in the soil profile. This can be achieved by applying effluent frequently in small amounts.

Notes: Seismic wavefield reconstruction is posed as an inversion problem where, from inadequate and incomplete data, we attempt to recover the data we would have acquired with a denser distribution of sources and receivers. A minimum weighted norm interpolation method is proposed to interpolate prestack volumes before wave-equation amplitude versus angle imaging.Synthetic and real data were used to investigate the effectiveness of our wavefield reconstruction scheme when preconditioning seismic data for wave-equation amplitude versus angle imaging.

Notes: We estimate the concentration of gas hydrate and free gas at an area located to the north of the Knipovich Ridge (western Svalbard margin). The method is based on P-wave velocities computed by reflection tomography applied to multicomponent ocean-bottom seismometer data. The tomographic velocity field is fitted to theoretical velocities obtained from a poro-elastic model based on a Biot-type approach (the interaction between the rock frame, gas hydrate and fluid is modelled from first physical principles). We obtain average hydrate concentrations of 7% and maximum free-gas saturations of 0.4% and 9%, depending on the saturation model.

Notes: Fluid-flow simulators used in the oil industry model the movement of fluids through a porous reservoir rock. These simulators either ignore coupling between the flow and concurring deformation of the solid rock frame or take it into account approximately, in the so-called loose or staggered-in-time mode. In contrast to existing simulators, the one we describe here fully couples two-phase (oil and water) flow to subsurface deformation and simultaneously accounts for all relevant physical phenomena. As such, our flow simulator inherently links time-dependent fluid pressures, saturations, permeabilities and flow velocities to stresses in the whole subsurface. These stresses relate to strains through the non-linear theory of elasticity, allowing us to model time-lapse changes in seismic velocities and anisotropy. The velocity variations manifest themselves in time shifts and reflection amplitudes that are conventionally measured from 4D seismic data. Changes in anisotropy produce time-dependent shear-wave splitting that can be used for monitoring the horizontal stresses.

Notes: Many joint inversion schemes use 1D forward modelling in the integrated interpretation of various geophysical data. In extending the joint inversion approach to the investigation of 2D structures, the discretization of the model parameters and the appropriate choice of the forward-modelling procedure play a very important role. In this paper, a hybrid seismic–geoelectric joint inversion method is proposed for the investigation of 2D near-surface geological structures. The electric and seismic models are coupled together through the use of common boundaries between the adjacent layers. Assuming a 2D model composed of homogeneous layers with curved boundaries, a fast ray-tracing algorithm is used for the calculation of refraction seismic traveltime data. In the geoelectric forward modelling, a locally 1D approximation is used. The boundary surfaces are written in the form of series expansion; the inversion algorithms are formulated for the expansion coefficients and the petrophysical parameters as unknowns. Two versions of the inversion method are proposed: in versions A and B, interval-wise constant functions and Chebyshev polynomials are, respectively, used as basis functions of the series expansion. The versions are tested by means of synthetic and in situ measured data. The tests show that both methods are stable and accurate.

Notes: Ghawar, the largest oilfield in the world, produces oil from the Upper Jurassic Arab-D carbonate reservoir. The high rigidity of the limestone–dolomite reservoir rock matrix and the small contrast between the elastic properties of the pore fluids, i.e. oil and water, are responsible for the weak 4D seismic effect due to oil production. A feasibility study was recently completed to quantify the 4D seismic response of reservoir saturation changes as brine replaced oil. The study consisted of analysing reservoir rock physics, petro-acoustic data and seismic modelling. A seismic model of flow simulation using fluid substitution concluded that time-lapse surface seismic or conventional 4D seismic is unlikely to detect the floodfront within the repeatability of surface seismic measurements. Thus, an alternative approach to 4D seismic for reservoir fluid monitoring is proposed. Permanent seismic sensors could be installed in a borehole and on the surface for passive monitoring of microseismic activity from reservoir pore-pressure perturbations. Reservoir production and injection operations create these pressure or stress perturbations. Reservoir heterogeneities affecting the fluid flow could be mapped by recording the distribution of epicentre locations of these microseisms or small earthquakes. The permanent borehole sensors could also record repeated offset vertical seismic profiling surveys using a surface source at a fixed location to ensure repeatability. The repeated vertical seismic profiling could image the change in reservoir properties with production.

Notes: For multivalued traveltime computation on dense grids, we propose a wavefront-orientated ray-tracing (WRT) technique. At the source, we start with a few rays which are propagated stepwise through a smooth two-dimensional (2D) velocity model. The ray field is examined at wavefronts and a new ray might be inserted between two adjacent rays if one of the following criteria is satisfied: (1) the distance between the two rays is larger than a predefined threshold; (2) the difference in wavefront curvature between the rays is larger than a predefined threshold; (3) the adjacent rays intersect. The last two criteria may lead to oversampling by rays in caustic regions. To avoid this oversampling, we do not insert a ray if the distance between adjacent rays is smaller than a predefined threshold. We insert the new ray by tracing it from the source. This approach leads to an improved accuracy compared with the insertion of a new ray by interpolation, which is the method usually applied in wavefront construction. The traveltimes computed along the rays are used for the estimation of traveltimes on a rectangular grid. This estimation is carried out within a region bounded by adjacent wavefronts and rays. As for the insertion criterion, we consider the wavefront curvature and extrapolate the traveltimes, up to the second order, from the intersection points between rays and wavefronts to a gridpoint. The extrapolated values are weighted with respect to the distances to wavefronts and rays. Because dynamic ray tracing is not applied, we approximate the wavefront curvature at a given point using the slowness vector at this point and an adjacent point on the same wavefront. The efficiency of the WRT technique is strongly dependent on the input parameters which control the wavefront and ray densities. On the basis of traveltimes computed in a smoothed Marmousi model, we analyse these dependences and suggest some rules for a correct choice of input parameters. With suitable input parameters, the WRT technique allows an accurate traveltime computation using a small number of rays and wavefronts.

Notes: We present a new least-squares migration method called least-squares wave-path migration. The proposed method combines an iterative conjugate-gradient solver with a stationary-phase wave-path migration operator. Numerical tests demonstrate that (i) least-squares wave-path migration is computationally more efficient than and almost as accurate as Kirchhoff least-squares migration, and (ii) many of the artefacts seen in wave-path migration images are suppressed after several conjugate-gradient iterations. Previous results have shown that 3D wave-path migration is up to 100 times faster than a standard 3D Kirchhoff migration, but sometimes at the cost of reduced quality. With the proposed least-squares wave-path migration method, the image quality in wave-path migration can be improved at an acceptable increase in computational cost.

Notes: Imaging pre-salt reflections for data acquired from the coastal region of the Red Sea is a task that requires prestack migration velocity analysis. Conventional post-stack time processing lacks the lateral inhomogeneity capability, necessary for such a problem. Prestack migration velocity analysis in the vertical time domain reduces the velocity–depth ambiguity that usually hampers the performance of prestack depth-migration velocity analysis. In prestack τ-migration velocity analysis, the interval velocity model and the output images are defined in τ (i.e. vertical time). As a result, we avoid placing reflectors at erroneous depths during the velocity analysis process and thus avoid inaccurately altering the shape of the velocity model, which in turn speeds up the convergence to the true model. Using a 1D velocity update scheme, the prestack τ-migration velocity analysis produces good images of data from the Midyan region of the Red Sea. For the first seismic line from this region, only three prestack τ-migration velocity analysis iterations were required to focus pre-salt reflections in τ. However, the second line, which crosses the first line, is slightly more complicated and thus required five iterations to reach the final, reasonably focused, τ-image. After mapping the images for the two crossing lines to depth, using the final velocity models, the placements of reflectors in the two 2D lines were consistent at their crossing point. Some errors occurred due to the influence of out-of-plane reflections on 2D imaging. However, such errors are identifiable and are generally small.

Notes: We have developed a method for imaging magnetic data collected for mineral exploration to yield the following structural information: depth, model type (structural index) and susceptibility. The active nature of mineral exploration data requires we derive the structural information from a robust quantity: we propose that the first- or second-order analytic-signal amplitude is suitably stable. The procedure is to normalize the analytic-signal amplitude by the peak value and then use non-linear inversion to estimate the depth and the structural index for each anomaly. In our field example, different results are obtained depending on whether we inverted for the first- or second-order analytic-signal amplitude. This is probably because the two-dimensional contact, thin sheet or horizontal cylinder models we have assumed are not appropriate. In cases such as these, when our model assumptions are not correct, the results should not be interpreted quantitatively, but they might be useful for giving a qualitative indication of how the structure might vary.With a priori information, it is possible to assume a model type (i.e. set the structural index) and generate estimates of the depth and susceptibility. These data can then be gridded and imaged. If a contact is assumed, the susceptibility contrast is estimated; for the dike model, the susceptibility-thickness is estimated; for the horizontal cylinder, the susceptibility-area is estimated. To emphasize that the results are dependent on our assumed model, we advocate prefixing any derived quantity by the term ‘apparent’.

Notes: An approach to determining the effective elastic moduli of rocks with double porosity is presented. The double-porosity medium is considered to be a heterogeneous material composed of a homogeneous matrix with primary pores and inclusions that represent secondary pores. Fluid flows in the primary-pore system and between primary and secondary pores are neglected because of the low permeability of the primary porosity. The prediction of the effective elastic moduli consists of two steps. Firstly, we calculate the effective elastic properties of the matrix with the primary small-scale pores (matrix homogenization). The porous matrix is then treated as a homogeneous isotropic host in which the large-scale secondary pores are embedded. To calculate the effective elastic moduli at each step, we use the differential effective medium (DEM) approach. The constituents of this composite medium – primary pores and secondary pores – are approximated by ellipsoidal or spheroidal inclusions with corresponding aspect ratios.We have applied this technique in order to compute the effective elastic properties for a model with randomly orientated inclusions (an isotropic medium) and aligned inclusions (a transversely isotropic medium). Using the special tensor basis, the solution of the one-particle problem with transversely isotropic host was obtained in explicit form.The direct application of the DEM method for fluid-saturated pores does not account for fluid displacement in pore systems, and corresponds to a model with isolated pores or the high-frequency range of acoustic waves. For the interconnected secondary pores, we have calculated the elastic moduli for the dry inclusions and then applied Gassmann's tensor relationships. The simulation of the effective elastic characteristic demonstrated that the fluid flow between the connected secondary pores has a significant influence only in porous rocks containing cracks (flattened ellipsoids). For pore shapes that are close to spherical, the relative difference between the elastic velocities determined by the DEM method and by the DEM method with Gassmann's corrections does not exceed 2%. Examples of the calculation of elastic moduli for water-saturated dolomite with both isolated and interconnected secondary pores are presented. The simulations were verified by comparison with published experimental data.

Notes: Shear-wave statics in marine seismic exploration data are routinely too large to be estimated using conventional techniques. Near-surface unconsolidated sediments are often characterized by low values of Vs and steep velocity gradients. Minor variations in sediment properties at these depths correspond to variations in the shear-wave velocity and will produce significant static shifts. It is suggested that a significant proportion of the shear-wave statics solution can be estimated by performing a separate high-resolution survey to target near-surface unconsolidated sediments. Love-wave, shear-wave refraction and geotechnical measurements were individually used to form high-resolution near-surface shear-wave velocity models to estimate the shear-wave statics for a designated survey line. Comparisons with predicted statics revealed that shear-wave statics could not be estimated using a velocity model predicted by substituting geotechnical measurements into empirical relationships. Empirical relationships represent a vast simplification of the factors that control Vs and are therefore not sufficiently sensitive to estimate shear-wave statics. Refraction measurements are potentially sensitive to short-wavelength variations in sediment properties when combined with accurate navigational data. Statics estimated from Love-wave data are less sensitive, and sometimes smoothed in appearance, since interpreted velocity values represent an average both laterally and vertically over the receiver array and the frequency–depth sensitivity range, respectively.For the survey site, statics estimated from near-surface irregularities using shear-wave refraction measurements represent almost half the total statics solution. More often, this proportion will be greater when bedrock relief is less.

Notes: Parsimonious post-stack migration is extended to three dimensions. By tracing single rays back along each incident wave direction (as determined by a local slant stack at the receivers), the ray tracing can be embedded in the migration. This approach significantly reduces the computer time and disk space needed because it is not necessary to build and save image time maps; 3D migration can be performed on a workstation or personal computer rather than using a supercomputer or cluster.The location of a reflector in the output image is defined by tracing a zero-offset ray to the one-way traveltime (the image condition); the orientation of the reflector is defined as a surface perpendicular to the raypath. The migration impulse response operator is confined to the first Fresnel zone around the estimated reflection point, which is much smaller than the large isochronic surface in traditional Kirchhoff depth migration. Additional efficiency is obtained by applying an amplitude threshold to reduce the amount of data to be migrated. Tests on synthetic data show that the proposed implementation of parsimonious 3D post-stack Kirchhoff depth migration is at least two orders of magnitude faster than traditional Kirchhoff migration, at the expense of slightly degraded migration image coherence. The proposed migration is expected to be a useful complement to conventional time migrations for fast initial imaging of subsurface structures and for real-time imaging of near-offset sections during data acquisition for quality control.

Notes: The semi-automated detection of objects has been quite successful in detecting various types of seismic object, such as chimneys. The same technique can be applied successfully to detect faults in 3D seismic data. We show that several different attributes – among others, similarity, frequency and curvature, all of which potentially enhance the visibility of faults – can be combined successfully by an artificial neural network. This results in a fault ‘probability’ cube in which faults are more continuous and noise is suppressed compared with single-attribute cubes. It is believed that the fault-cube can be improved further by applying image-processing techniques to enhance the fault prediction.

Notes: Filters for migrated offset substacks are designed by partial coherence analysis to predict ‘normal’ amplitude variation with offset (AVO) in an anomaly free area. The same prediction filters generate localized prediction errors when applied in an AVO-anomalous interval. These prediction errors are quantitatively related to the AVO gradient anomalies in a background that is related to the minimum AVO anomaly detectable from the data. The prediction-error section is thus used to define a reliability threshold for the identification of AVO anomalies. Coherence analysis also enables quality control of AVO analysis and inversion. For example, predictions that are non-localized and/or do not show structural conformity may indicate spatial variations in amplitude–offset scaling, seismic wavelet or signal-to-noise (S/N) ratio content. Scaling and waveform variations can be identified from inspection of the prediction filters and their frequency responses. S/N ratios can be estimated via multiple coherence analysis.AVO inversion of seismic data is unstable if not constrained. However, the use of a constraint on the estimated parameters has the undesirable effect of introducing biases into the inverted results: an additional bias-correction step is then needed to retrieve unbiased results. An alternative form of AVO inversion that avoids additional corrections is proposed. This inversion is also fast as it inverts only AVO anomalies. A spectral coherence matching technique is employed to transform a zero-offset extrapolation or near-offset substack into P-wave impedance. The same technique is applied to the prediction-error section obtained by means of partial coherence, in order to estimate S-wave velocity to P-wave velocity (VS/VP) ratios. Both techniques assume that accurate well ties, reliable density measurements and P-wave and S-wave velocity logs are available, and that impedance contrasts are not too strong. A full Zoeppritz inversion is required when impedance contrasts that are too high are encountered. An added assumption is made for the inversion to the VS/VP ratio, i.e. the Gassmann fluid-substitution theory is valid within the reservoir area. One synthetic example and one real North Sea in-line survey illustrate the application of the two coherence methods.

Notes: A local plane-wave approach of generalized diffraction tomography in heterogeneous backgrounds, equivalent to Kirchhoff summation techniques when applied in seismic reflection, is re-programmed to act as repeated synthetic aperture radar (SAR) imaging for seismic prestack depth migration. Spotlight-mode SAR imaging quickly provides good images of the electromagnetic reflectivity of the ground via fast Fourier transform (FFT)-based signal processing. By calculating only the Green's functions connecting the aircraft to the centre of the illuminated patch, scattering structures around that centre are also recovered. SAR technology requires us to examine seismic imaging from the local point of view, where the quantity and quality of the available information at each image point are what are important, regardless of the survey geometry. When adapted to seismics, a local image of arbitrary size and sampling is obtained by FFT of seismic energy maps in the scattering wavenumber domain around each node of a pre-calculated grid of Green's functions. These local images can be used to generate a classic prestack depth-migrated section by collecting only their centres. However, the local images also provide valuable information around the centre, as in SAR. They can therefore help to pre-analyse prestack depth migration efficiently, and to perform velocity analysis at a very low cost. The FFT-based signal-processing approach allows local, efficient and automatic control of anti-aliasing, noise and resolution, including optimized Jacobian weights. Repeated local imaging could also be used to speed up migration, with interpolation between local images associated with a coarse grid of Green's functions, as an alternative to interpolation of Green's functions. The local images may, however, show distortions due to the local plane-wave approximation, and the velocity variations across their frame. Such effects, which are not necessarily a problem in SAR, should be controlled and corrected to further enhance seismic imaging. Applications to realistic models and to real data show that, despite the distortion effects, the local images can yield similar information to prestack depth migration, including common-image-point gathers for velocity analyses and AVO/AVA effects, at a much lower cost when a small target is considered.

Notes: Vibroseis is a method that imparts coded seismic energy into the ground. The energy is recorded with geophones and then processed using the known (coded) input signal. The resulting time-domain representation of vibroseis data is an impulsive wavetrain with wavelet properties consistent with the coded input signal convolved with the earth's reflectivity series. Historically, vibratory seismic surveys collect data from one source location at a time, summing one or more sources at each location. We present a method of designing orthogonal sweeps using the concept of combisweeps. The orthogonal sweeps allow simultaneous recording and later separation of two or more unique source locations. Orthogonality of sweeps permits separation of the data into unique source-location field records by a conventional correlation procedure. The separation power of the orthogonal sweeps is demonstrated by a comparison between separated data and data acquired with one vibrator. Separation noise was at a negligible level for our demonstration data sets when two vibrators were located 50 m to 200 m apart. Coincident generation and recording of two vibroseis sweeps at different locations would allow almost double the amount of data to be recorded for a given occupation time and requires only half the storage medium.

Notes: Shales are a major component of sedimentary basins, and they play a decisive role in fluid flow and seismic-wave propagation because of their low permeability and anisotropic microstructure. Shale anisotropy needs to be quantified to obtain reliable information on reservoir fluid, lithology and pore pressure from seismic data, and to understand time-to-depth conversion errors and non-hyperbolic moveout. A single anisotropy parameter, Thomsen's δ parameter, is sufficient to explain the difference between the small-offset normal-moveout velocity and vertical velocity, and to interpret the small-offset AVO response. The sign of this parameter is poorly understood, with both positive and negative values having been reported in the literature. δ is sensitive to the compliance of the contact regions between clay particles and to the degree of disorder in the orientation of clay particles. If the ratio of the normal to shear compliance of the contact regions exceeds a critical value, the presence of these regions acts to increase δ, and a change in the sign of δ, from the negative values characteristic of clay minerals to the positive values commonly reported for shales, may occur. Misalignment of the clay particles can also lead to a positive value of δ. For transverse isotropy, the elastic anisotropy parameters can be written in terms of the coefficients W200 and W400 in an expansion of the clay-particle orientation distribution function in generalized Legendre functions. For a given value of W200, decreasing W400 leads to an increase in δ, while for fixed W400, δ increases with increasing W200. Perfect alignment of clay particles with normals along the symmetry axis corresponds to the maximum values of W200 and W400, given by 〈inlineGraphic alt="inline image" href="urn:x-wiley:00168025:GPR495:GPR_495_mu1" location="equation/GPR_495_mu1.gif"/〉 and 〈inlineGraphic alt="inline image" href="urn:x-wiley:00168025:GPR495:GPR_495_mu2" location="equation/GPR_495_mu2.gif"/〉. A comparison of the predictions of the theory with laboratory measurements shows that most shales lie in a region of the (W200, W400)-plane defined by W400/W200≤Wmax400/Wmax200.