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: A dairy system simulator, Dairy_sim, was designed to assess the interactions between climate and management in spring-calving milk production systems based on the grazing of grass pastures. The simulator comprises three main components: a grass herbage growth model, an intake and grazing behaviour model, and a nutrient demand model. The simulator was initially parameterized using the Irish National Dairy Blueprint. Sensitivity analysis indicated that the simulator was most sensitive to stocking rate, milk output per cow and nitrogen fertilizer inputs, but less sensitive to other variables. Field data from four grazing systems were used to test the simulator and it was concluded that Dairy_sim was suitable for evaluating the interaction of climate and management for rotational grazing dairy systems based on perennial ryegrass pastures with Friesian cows. The simulator, Dairy_sim, was then used to evaluate the effects of the regional climates of Ireland on system management. The results indicated that, between regions, herbage production at the same input of nitrogen may vary proportionally by 0·10 and that the length of the grazing season may vary by 0·25. It was concluded that the simulator could be a useful tool for developing region-specific dairy production blueprints.

Notes: Cycling of potassium in grassland systems has received relatively little attention in research and practice in recent years. Balanced nutrient systems require consideration of nutrients other than nitrogen (N). Potassium (K) is needed in large amounts and is closely related to N nutrition. In intensive dairy farming, surpluses of K arise from the input of concentrates and fertilizer and are returned to the grassland and may lead to increasing K content in the soil. Organic farming, on the other hand, is characterized by limitations in input of nutrient sources and quantities. Leaching of K from grassland is usually low, but high levels of available soil K, high K input from fertilizer or at urine patches lead to increasing losses. High K inputs have a negative influence on Mg and Ca uptake by plants and can cause accelerated leaching of these cations. High levels of K have been associated with inducing nutrition-related dairy cow health problems such as milk fever (hypocalcaemia) and grass tetany (hypomagnesaemia). This review gives an overview of the cycling of potassium and related cations in grassland systems especially with regard to leaching losses and identifies limitations to knowledge.

Notes: Advancing maize crop maturity is associated with changes in ear-to-stover ratio which may have consequences for the digestibility of the ensiled crop. The apparent digestibility and nitrogen retention of three diets (Early, Mid and Late) containing maize silages made from maize of advancing harvest date [dry matter (DM) contents of the maize silages were 273, 314 and 367 g kg−1 for the silages in the Early, Mid and Late diets respectively], together with a protein supplement offered in sufficient quantities to make the diets isonitrogenous, were measured in six Holstein–Friesian steers in an incomplete Latin square design with four periods. Dry-matter intake of maize silage tended to be least for the Early diet and greatest for the Medium diet (P = 0·182). Apparent digestibility of DM and organic matter did not differ between diets. Apparent digestibility of energy was lowest in the Late diet (P = 0·057) and the metabolizable energy concentrations of the three silages were calculated as 11·0, 11·1 and 10·6 MJ kg−1 DM for the Early, Medium and Late diets respectively (P = 0·068). No differences were detected between diets in starch digestibility but the number of undamaged grains present in the faeces of animals fed the Late diet was significantly higher than with the Early and Mid diets (P = 0·006). The apparent digestibility of neutral-detergent fibre of the diets reduced significantly as silage DM content increased (P = 0·012) with a similar trend for the apparent digestibility of acid-detergent fibre (P = 0·078). Apparent digestibility of nitrogen (N) was similar for the Early and Mid diets, both being greater than the Late diet (P = 0·035). Nitrogen retention did not differ between diets. It was concluded that delaying harvest until the DM content is above 300 g kg−1 can negatively affect the nutritive value of maize silage in the UK.

Notes: The effect of offering a total mixed ration of silage and concentrate (proportionately 0·44 silage) system [indoor feeding system (IF)] was compared with grazing at a high daily herbage allowance with a low level of concentrate supplementation [early grazing system (EG)] in early spring on the performance of spring-calving dairy cows in Ireland. Sixty-four spring-calving Holstein–Friesian dairy cows (mean calving date, 2 February) were allocated to one of two systems between 16 February and 4 April 2004. An equal number of primiparous and multiparous cows were assigned to each system. The dairy cows on the IF system were housed for a 7-week period and offered a diet of 10·9 kg DM cow−1 d−1 (s.d. 2·3) of concentrate, the remainder of the diet was 8·6 kg DM cow−1 d−1 (s.d. 1·9) of grass silage. The dairy cows on the EG system were offered a mean daily herbage allowance of 15·1 kg DM cow−1 d−1 (s.d. 3·7) and were supplemented with 3·0 kg DM cow−1 d−1 (s.d. 1·0) of concentrate. There was no difference in milk yield between the two systems but the cows in the EG system had a higher milk protein concentration (2·9 g kg−1) and a higher milk protein yield than in the IF system. Milk fat concentration was higher for cows in the IF than EG system (3·0 g kg−1). There was no difference in total daily dry-matter intake between the systems, measured in week 6 of the study. Mean live weight of the cows in the IF system was greater than in the EG system. The results of the study suggest that a slightly greater performance can be achieved by a system offering a high daily herbage allowance to spring-calving dairy cows in early lactation compared with a system offering a total mixed ration containing a high proportion of concentrate with grass silage.

Notes: The relative competitive ability of Artemisia frigida Willd. (C3, semi-shrub) and Cleistogenes squarrosa (Trin.) Keng (C4, perennial bunchgrass), grown in a greenhouse, was studied under three nitrogen fertilizer levels, three defoliation intensities and five competitive replacement series in a factorial arrangement. Relative yield total and the respective crowding coefficients for each species were calculated with reference to aboveground, below-ground and total biomass. The index of aggressivity was also determined relative to total biomass. Cleistogenes squarrosa was a consistently weak competitor at all nitrogen fertilizer and clipping levels compared with A. frigida. Artemesia frigida had a higher relative yield when grown in mixtures with C. squarrosa than in monocultures, whereas the yield of C. squarrosa was lower in mixtures with A. frigida than in monocultures. The relative competitive ability of A. frigida over C. squarrosa declined significantly at different clipping intensities with increasing nitrogen fertilizer. The interaction of nitrogen fertilizer with clipping was significant, although the clipping effect on aggressivity was not significant.

Notes: An experiment was conducted to determine the genetic variation among diploid perennial ryegrass (Lolium perenne L.) varieties for sward structural characteristics considered to be important for intake by cattle. Assessments were made between June and September in 2000 and 2001. Six varieties (Abergold, Respect, Agri, Herbie, Barezane and Barnhem) were subjected to a cutting experiment where swards were cut after 3 to 4 weeks of regrowth during the growing season. The variables, measured in three 2-week periods, were herbage mass of dry matter (DM), sward surface height (SSH), bulk density, proportion of green leaf, tiller density, tiller weight, extended tiller height, length of sheath and length of leaf blade. Significant differences among varieties were found in both years for herbage mass of DM, SSH, bulk density, proportion of green leaf, tiller density, tiller weight and length of sheath. The results show that there is significant genetic variation among diploid perennial ryegrass varieties for sward characteristics important for intake during grazing.

Notes: The effects of one severe winter-grazing of lucerne were studied over 3 years in an experiment in the Ebro Valley, Spain. In this region the crop is harvested six to seven times per season and winter grazing is a traditional practice. On average, winter-grazing reduced the yield at the first harvest in spring by 200 kg dry matter (DM) ha−1. This limited yield reduction of 0·06 was accompanied by an increase in the proportion of lucerne in the herbage DM from 0·54 to 0·62, a reduction in the proportion of weeds from 0·39 to 0·36, and a reduction in the proportion of dead material from 0·06 to 0·02. The crude protein concentration and the in vitro DM digestibility increased by 20 g kg−1 DM and 0·03, respectively. The traditional practice, i.e. of grazing lucerne with sheep once in the winter season, results in only a limited reduction in yield in the spring. In addition, the spring crop has a higher nutritive value.

Notes: Soil organic carbon (SOC) stocks in different terrestrial ecosystems have already been calculated. However, SOC stocks of grassy verges and grass-covered urban areas have never been reported. In this study the total grass-covered area of verges and urban areas in different agro-pedological regions of Flanders, Belgium was calculated and SOC stocks estimated. The total grassy area of verges along roads, waterways and railways was calculated to be 18 027 ha, 9530 ha for gardens and parks and 2360 ha for recreations areas. Total SOC stock for grass-covered verges, extrapolated to a depth of 1 m, was estimated to be 3520 kt SOC. Total SOC stock for grassy vegetation in urban areas was 1738 kt SOC. These total SOC stocks equate, respectively, to 0·10 and 0·05 of the total SOC stock in Flemish grassland.

Notes: The survival of pathogenic and spoilage micro-organisms in soil and on grass fertilized with spiked anaerobic digestion residue (ADR) was investigated in a climate chamber during periods of up to 56 d. In addition, the survival of these organisms over time was investigated during ensiling of grass at 390 g dry matter (DM) kg−1 or 610 g DM kg−1. Micro-organisms included in these studies were: Clostridium tyrobutyricum, Salmonella serotype Typhimurium, Listeria monocytogenes, Campylobacter jejuni, Escherichia coli, Talaromyces emersonii, Byssochlamys nivea, Porcine parvovirus and Swine vesicular disease virus. Soil and grass still contained high numbers of E. coli, Cl. tyrobutyricum and T. emersonii (1·9–5·4 log10) 49 and 56 d after fertilization with spiked ADR. Listeria monocytogenes and S. Typhimurium were generally found in the samples. This indicates that, within this time span, there is a risk of silage contamination by bacteria, moulds and viruses present in ADR spread on grassland. An increase in DM content of the crop decreased its ensilability as measured by pH and short-chain fatty acid content. However, no clear differences were found in survival of pathogenic and spoilage micro-organisms between the two silages with different DM contents, regardless of storage time. The lack of moisture and oxygen was probably sufficient to cause the reduction in vegetative bacteria in the 610 g DM kg−1 silage. However, bacterial and fungal spores and the viruses studied were not significantly inactivated by ensiling at these high DM contents and could, therefore, pose a health risk to farm animals fed on the silage if present in ADR applied to crops prior to ensilage.

Notes: Cattle and sheep can create and maintain a mixture of relatively tall and short patches in grass swards through selective grazing. In swards that are grazed by cattle this heterogeneous structure can result in the frequencies of height measurements having a skewed distribution that has variously been better described by the double-normal distribution the gamma distribution and the Weibull distribution than by the more common normal distribution. The fit of these statistical distributions, and the adequacy of the potentially useful log-normal distribution, to sward height frequencies were tested in sown temperate swards grazed by sheep and compared within a single sward. It was concluded that the single-normal and Weibull distributions were inadequate and that overall the log-normal and gamma distributions had the best fit to the measurements.

Notes: Descriptions of anthropocentric phenologies of flower colour were obtained for all thirteen British mesotrophic grassland communities, defined by the National Vegetation Classification. They were obtained by combining information on species abundance with data on flowering time and flower colour. Although the results are crude estimates, they clearly demonstrate that grassland communities differ in their floral phenologies with yellow-flowered species being significantly most abundant, and blue-flowered species being least abundant than the colours of other flowering species in British grasslands. The peak months for flowering were June and July. Significant interactions were observed between community type and flower colour; community type and month; and flower colour and month. These interactions may be explained by ecological differences between the communities, for example hay meadows flowered earlier than other grassland types. It is argued that other differences between communities in their floral phenologies, for example the absence of blue flowers from the agriculturally improved MG7 and MG11 swards, may support the hypothesis that the limited availability of pollination niches may affect plant community structure.

Notes: The contention that agricultural grassland within the UK is spatially and structurally uniform across wide differences in livestock farming intensity is examined. Total nitrogen (N) input was used as a surrogate for management intensity. Farms with average annual N inputs of 〉200 kg ha−1 were categorized as highly intensive, those with N inputs of 50–200 kg N ha−1 as moderately intensive and those with inputs of 〈50 kg N ha−1 were categorized as extensive. Four farms within each management category were selected in two discrete regions: one in south-west England, typifying a landscape dominated by agricultural grassland, the other in south-east England with a more mixed-farming landscape. Specific N input and management data, and sward botanical composition and structure, were obtained on four fields per farm. The south-west region had a higher proportion of farms with dairy cattle than in the south-east England region. The average N input and stocking rates were higher within the intensive-management category farms in south-west England than in south-east England. Old permanent grasslands contained variable amounts of sown species, such as Lolium perenne, and were dominated by a few (generally 〈5) non-sown grass species. Grassland that received 〉75 kg N ha −1 year−1 contained 〈7 plant species and 〈3 forb species m−2. Grassland communities with 〉12 plant species and 〉5 forb species m−2 were only found in study fields that received 〈50 kg N ha−1 year−1. Grasslands with 〉10 forb species m−2 were found only in sites receiving 〈15 kg N ha−1. Of importance at a landscape scale was ubiquity of species-poor and similar grassland plant communities across the management intensity range within and between livestock farms. Grassland in the extensive management category was similar in structure in terms of sward height and/or sward density to more intensively managed grasslands. This study revealed, albeit in two lowland sample regions in the UK, the ubiquity of grass-dominated, species-poor and structurally uniform grasslands, irrespective of apparently broad differences in farming intensity. The consequences of this spatial uniformity on grassland biodiversity are likely to be profound.

Notes: The effects of offering a range of grass silages and mixtures of grass and maize silages on the intake of beef cattle were studied. Four grass silages (GS1, GS2, GS3 and GS4) were used. Grass silage 1 was ensiled from a second regrowth in mid-late September and treated with an inoculant additive. Grass silages 2, 3 and 4 were ensiled, without additive, from a primary regrowth harvested in early July, late May and mid-June respectively. Wilting periods were 8, 30, 36 and 36 h for GS1, GS2, GS3 and GS4 respectively. Grass silages 1, 2 and 3 were precision chopped and ensiled in bunker silos, while GS4 was ensiled in round bales. The DM content (g kg−1) and starch concentration (g kg−1 DM) of the three maize silages (MS1, MS2 and MS3) used were 256 and 128, 256 and 184, and 402 and 328 for MS1, MS2 and MS3 respectively. Seventy-two Charolais and Limousin cross-bred steers were used in a changeover design with two 4-week periods. The study consisted of sixteen treatments incorporating the four grass silages fed alone and with the three maize silages arranged as a 4 × 4 factorial design. The grass silage and maize silage mixtures were offered in a ratio of 0·60:0·40 (DM basis) once daily using individual Calan gates. All silages were offered ad libitum with 3 kg per head per day of a concentrate supplement. Dry matter and metabolizable energy (ME) intakes were highest with diets based on grass silage GS4 compared with diets containing the other grass silages. Metabolizable energy intakes of diets containing no maize silage, and those based on MS1 and MS2, were similar (P 〉 0·05) but lower than that of diets containing MS3. Only limited increases were found in DM and ME intakes with the inclusion of maize silage in grass silage-based diets while offering high-quality grass silage (assessed in terms of DM content, and fibre and N concentrations) promoted high voluntary intakes.

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: The fermentation quality of small-bale silage and haylage for feeding to horses in Sweden, and using a conventional high-density hay baler, was investigated in two experiments. Treatments studied were use of additives (inoculants containing lactic acid bacteria and a chemical additive consisting of hexamethylenetetramine, sodium nitrate, sodium benzoate and sodium propionate), the influence of dry-matter (DM) content of wilted herbage and the effect of number of stretch film layers on fermentation pattern and aerobic stability. All silages and haylages were made from predominantly Timothy swards and were well fermented as indicated by low levels of ammonia and butyric acid. Values of pH were higher and concentrations of organic acids were lower in haylages than in the silages. This was not considered to be indicative of a poor fermentation in the haylage but of a restricted fermentation due to the high DM content of the herbage. The additives enhanced aerobic storage stability because of inhibition of mould growth. The only statistically significant effect of varying the number of stretch film layers was a higher content of CO2 inside the bales when ten layers of stretch film were applied compared with six layers.

Notes: Factors reducing yield of organic white clover seed production were investigated in Denmark. In thirty-six fields from different regions, density of flower heads, and weevils of the genera Apion and Hypera, hatching from flower heads, were assessed. In fifteen of these fields, additional measurements were made in order to calculate potential yield and yield-reducing factors. Flower heads had a mean of ninety-two florets per head, of which proportionately 0·59 were intact, 0·08 were not pollinated and 0·34 damaged by weevils; 0·16 of the florets contained larvae of Apion species. In intact florets the mean number of seeds was 3·6. Density and size of flower heads were correlated with low density of flower heads leading to significantly more florets per head. The number of seeds per intact floret was correlated with the proportion of unpollinated florets. The number of insect-damaged florets was influenced by both larvae of Apion in and the numbers of Hypera nigrirostris (Fabr.) hatched from flower heads. One larva of H. nigrirostris caused approximately ten times the damage of a larva of Apion. Densities of Apion and Hypera were influenced by location, with fields with adjacent weevil sources being most at risk. An estimate of potential seed yield showed that, even under good harvest conditions, only a small part of potential seed yield is realized and that under suboptimal harvest conditions seed loss is high. The main factors responsible for low yields in Danish organic white clover seed production were identified as unfavourable harvest conditions and weevil damage.

Notes: Demographic and biomass allocation patterns of Agropyron cristatum were measured on steppe and dune sites in the Hunshandake desert of North China in 2001 and 2002. Total plant population, reproductive shoot densities and its differentiation rates were significantly higher in the steppe sites in both years. Plant heights for both vegetative and reproductive shoots were greater in the year with the higher rainfall. The dune sites had a higher biomass allocation to vegetative shoots and roots, while the steppe sites had a higher biomass allocation to reproductive shoots and seed production. It is suggested that the population demography and biomass allocation of the species responded to the differences in the soil variables in the steppe and dune sites.

Notes: The objective was to determine how the frequency and intensity of defoliation of tall oat grass [Arrhenatherum elatius (L.) Presl.], affected (i) annual dry-matter (DM) accumulation, (ii) seasonal distribution of the accumulated forage, and (iii) morphological traits related to persistence and vigour. Data were collected twice a year throughout three growing years, beginning the second year after establishment. Two frequencies (A: each time modal height of plants reached 20 cm; B: each time modal height of plants reached 40 cm) and two defoliation intensities (I: 5 cm; II: 10 cm stubble) were arranged in a completely randomized block design with a fully factorial combination and four replicates. Measurements were made of annual DM accumulation, seasonal distribution of accumulated DM and morphological traits related to persistence and vigour, i.e. number of tillers m−2, number of tussocks m−2, crown diameter, and crown area of tussocks. The infrequent defoliation (B) increased annual DM accumulation besides maintaining levels of persistence and vigour. The most productive treatment (BI) had a significantly higher annual DM accumulation than treatments AI and AII. Frequent and severe defoliations (treatment AI) led to plant depletion, which reduced crown diameter, number of tillers m−2 and crown area at the end of the third year compared with infrequent but severe defoliation (treatment BI). Annual DM accumulation was intermediate for less intense defoliations (AII), which produced forage during winter, but resulted in less vigorous plants by the end of the experiment. Crown area was a good indicator of tiller number, particularly under frequent defoliation.

Notes: Field survey and satellite image processing methods were used to estimate the total available forage over an area of 95 034 ha in north-eastern Syria, and to assess grazing impact on the area. The above-ground plant biomass was measured by a quadrat method at three sites in each of eight vegetation classes. Available forage was measured by excluding woody parts of shrubs from the whole aerial plant parts. The total above-ground plant biomass and available forage were estimated by extrapolating the measured point data to the whole target area using classified vegetation data by satellite image processing. Grazing impact was assessed by calculating the differences between the total available forage at the end of growing season and the end of dry season. The values for the estimated total available forage (s.e. of mean) in the area were 55 628 000 (12 920 000) kg DM and 30 007 000 (2 437 000) kg DM at the end of growing season and dry season respectively. Although the area of the cereal fields covered only 0·315 of the area, about 0·69 and 0·82 of the available forage existed in the harvested cereal fields at the ends of growing season and dry season respectively. The integration of cereal fields and rangeland is a normal land use system for livestock management in the area. The higher cover of herbaceous vegetation types showed higher grazing impacts which reduced the total available forage at the end of the growing season by 0·817 (0·199) at the end of the dry season. Although these dense herbaceous vegetation types could possibly produce more available forage, they would incur more intensive grazing impact. On the contrary, lighter grazing impact would occur with a higher cover of shrub vegetation types. The importance of maintaining plant cover over the rangeland area to protect the land against soil erosion is stressed.

Notes: Introgression of reproductive traits from the annual, profuse flowering, ball clover (Trifolium nigrescens Viv.) into white clover (Trifolium repens L.) is one breeding strategy to improve seed yields of T. repens that must be achieved without sacrificing agronomic performance and persistency under grazing. The yield and persistency of hybrids between white clover and T. nigrescens were compared under rotational sheep grazing over three harvest years. The hybrids included the backcross (BC) 2 and 3 generations produced using white clover as the recurrent parent. The large-leaved T. repens variety Olwen, medium-leaved varieties AberDai and Menna and the small-leaved variety S184 were sown as controls. Hybrids and control varieties were sown with a perennial ryegrass companion; between April and the end of October in each harvest year the plots were rotationally grazed with sheep with clover and perennial ryegrass (DM) yield and the proportion of clover present measured over the growing season. The clover and total DM yields of the BC2 and BC3 were generally comparable with the small- and medium-leaved varieties within the experiment and significantly greater than the yields of the large-leaved variety Olwen. Throughout the 3 years of the experiment the BC2 maintained a clover content above 0·30 and comparable with the small-leaved varieties, while the clover content of the BC3 was comparable with the small- and medium-leaved varieties in the first and third harvest years. No significant difference in perennial ryegrass production was observed when grown with the backcrosses or the control varieties. Differences in stolon and growing-point density were observed at the end of the experiment with the density of the BC2 and BC3 less than the small-leaved variety S184 but, in common with the medium-leaved varieties, greater than the large-leaved variety Olwen. The implication of these results for the use of this material in future experiments and in the white clover breeding programme is discussed.

Notes: Three experiments were conducted to determine the association between leaf number per tiller at defoliation, water-soluble carbohydrate (WSC) concentration and herbage mass of juvenile ryegrass plants when grown in a Mediterranean environment. Seedlings of ryegrass were grown in nursery pots arranged side-by-side and located outside in the open-air to simulate a mini-sward in Experiments 1 and 2, and a mixture of annual ryegrass and subterranean clover (Trifolium subterraneum L.) was grown in a small plot field study in Experiment 3. Swards were defoliated mechanically with the onset of defoliation commencing within 28 d of germination. Frequency of defoliation ranged from one to nine leaves per tiller, whilst defoliation height ranged from 30 mm of pseudostem height that removed all leaf laminae in Experiment 1, to 50 mm of pseudostem height with some leaf laminae remaining post-defoliation in Experiments 2 and 3.A positive relationship between herbage mass of ryegrass, WSC concentration and leaf number per tiller at defoliation was demonstrated in all experiments. In Experiment 1, the herbage mass of leaf, pseudostem and roots of tillers defoliated at one leaf per tiller was reduced to 0·10, 0·09 and 0·06 of those tillers defoliated less frequently at six leaves per tiller. However, the reduction in herbage mass from frequent defoliation was less severe in Experiment 2 and coincided with a 0·20 reduction in WSC concentration of pseudostem compared with 0·80 measured during Experiment 1. In Experiment 3, the highest harvested herbage mass of ryegrass occurred when defoliation was nine leaves per tiller. Although the harvested herbage from this sward contained senescent herbage, the in vitro dry-matter digestibility of the harvested herbage did not differ significantly compared with the remaining treatments that had been defoliated more frequently.Leaf numbers of newly germinated ryegrass tillers in a Mediterranean environment were positively associated with WSC concentration of pseudostem and herbage mass. A minimum period of two to three leaf appearances was required to restore WSC concentrations to levels measured prior to defoliation thereby avoiding a significant reduction in herbage mass. However, maximum herbage mass of a mixed sward containing ryegrass and subterranean clover was achieved when defoliation was delayed to nine leaves per tiller.

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).