ALBERT

Notes: Dissolved organic nitrogen and carbon (DOC) are significant in the C and N cycle in terrestrial ecosystems. Little is known about their dynamics in the field and the factors regulating their concentrations and fluxes. We followed the fluxes and concentrations of the two in a Norway spruce (Picea abies (L.) Karst.) forest ecosystem in Germany from 1995 to 1997 by sampling at fortnightly intervals. Bulk precipitation, throughfall, forest floor percolates from different horizons and soil solutions from different depths were analysed for major ions, dissolved organic N and DOC. The largest fluxes and concentrations were observed in percolates of the Oi layer, which contain amino N and amino sugar N as the major components. The average ratio of dissolved organic C to N in forest floor percolates corresponded to the C/N ratio of the solid phase. Concentrations and fluxes were highly dynamic with time and decreased with depth. The largest fluxes in forest floor percolates occurred when the snow melted. The concentrations and fluxes of dissolved organic N were significantly correlated with DOC, but the correlation was weak, indicating different mechanisms of release and consumption. The dynamics of dissolved organic N and DOC in forest floor percolates were not explained by pH and ionic strength of the soil solution nor by the water flux, despite large variations in these. Furthermore, the release of these fractions from the forest floor was not related to the quality and amount of throughfall. Concentrations of dissolved organic N in forest floor percolates increased with soil temperature, while temperature effects on DOC were less pronounced, but their fluxes from the forest floor were not correlated with temperature. In the growing season concentrations of both dissolved organic N and C in forest floor percolates decreased with increasing intensity of throughfall. Thus, the average throughfall intensity was more important than the amount of percolate in regulating their concentrations in forest floor percolates. Our data emphasize the role of dissolved organic N and DOC in the N and C cycle of forest ecosystems.

Notes: The migration of colloidal particles in the unsaturated zone may enhance leaching of sorbing contaminants to surface waters and groundwaters. This paper describes a simple model of particle leaching and translocation based on the dual-porosity model MACRO. The model includes descriptions of processes such as ‘source-limited’ particle detachment due to the kinetic energy of rain, replenishment of the depleted store of particles up to a maximum value determined by the amount of dispersible clay in the soil, and convective transport of particles in macropores subject to a filter sink term varying as a function of pore water velocity. Estimates of model parameter values are obtained by comparing simulations with measurements of particle concentrations in the discharge from tile drains made in a silty clay soil in southwest Sweden. Calibrated in this way, the model is shown to reproduce satisfactorily the observations. A sensitivity analysis suggests that the critical parameters in the model are those related to particle filtering during vertical transport in soil macropores.

Notes: The stabilization of organic matter in soil by interaction with aluminium (Al) or allophane is important in maintaining soil quality, and in retarding the decomposition of soil organic matter. Complexation of Al by soil organic matter may also ameliorate Al toxicity. Here we use 13C-NMR spectroscopy to assess the interaction of soil organic matter with both Al and allophane in two poorly drained podzols containing only trace amounts of iron. The 13C-NMR spectrum of the subsoil of the allophane-rich One Tree Point podzol shows an intense peak at 179 p.p.m., assigned to carbon in carboxylic acids. This peak shifts to 177 p.p.m. after removal of allophane (11% of the soil mass) by treatment with HF. We infer that the carboxyl groups in the organic matter are bonded to structural Al on the surface of allophane spherules. In the non-allophanic Te Kopuru podzol, on the other hand, the organic matter apparently interacts with Al ions in the soil solution. This soil also has more aromatic carbon and fewer carbons in carboxyl and carbohydrate structures than the allophanic sample. There is an indication that allophane stabilizes carbohydrate groups as well as carboxyl groups.

Notes: Soil profiles under beech, spruce and a grassland have been analysed to study the evolution of natural n-alkanes in pollution-free ecosystems. The soils had all developed on granitic bedrock, at an altitude of 1300–1500 m in the region of Mont-Lozère (southern Massif-Central, France). In contrast to the grassland soil, the two forest soils both possessed a well-developed acidic moder humus-type horizon. This could be subdivided as follows: fresh litter (OL), fragmentation (OF) and humification (OH) layers; two litters, one fresh (OL1) and one old (OL2) could actually be distinguished in the beech forest soil. The n-alkane signature of the parent plants was preserved in the top litter. Immediately underneath, in the OF layer(s) the original n-alkane signatures were progressively but rapidly replaced by a common signature composed of n-C27 and n-C25 with larger proportions of the former than of the latter. These two hydrocarbons were most probably produced in situ by fungi. These results appear to illustrate the action of soil microorganisms which metabolize the inherited n-alkanes and produce new compounds of the same family. Unlike the alkanes and the low molecular weight fatty acids ≤ C20 (which increase greatly in the OL2 layer under beech as a result of intense microbial activity), the heavy fatty acids (〉 C20) show no significant change in the organic horizon.

Notes: This paper reports a procedure for determining the content of strongly fixed NH4+ in soil. The procedure consists of a Kjeldahl digestion followed by an acid attack of the residue with a 5 m HF:1 m HCl solution. Distillations after each of the two treatments recover different forms of NH4+. The procedure was tested on fine earth (〈 2 mm) and skeleton (〉 2 mm) fractions of two forest soils developed on sandstone parent material. In both soil fractions we evaluated three different forms of NH4+-N: (i) Kjeldahl, (ii) non-exchangeable and (iii) micaceous. The last is located in the interlayer of mica flakes larger than 50 μm that resist the Kjeldahl digestion and is considered strongly fixed. The total NH4+-N content of a soil is obtained by the summation of the Kjeldahl and the micaceous NH4+-N. In the soils under consideration, the micaceous form prevails in the skeleton because this fraction is richer in micas of sand size (〉 50 μm). Following the proposed procedure, we found that micas (muscovite and biotite) contain about 3000 mg kg–1 of NH4+-N in the interlayer. The presence of micaceous NH4+-N in soil is generally ignored because the skeleton is usually excluded from analyses, and the micas larger than 50 μm cannot be dissolved by the Kjeldahl treatments. The micaceous NH4+ is the least extractable form of NH4+-N, and we infer that it is the least available to plants.

Notes: Hardsetting of soil comprises two contrasting, unfavourable physical conditions at different water contents, namely extreme hardness when dry and complete slaking when wet. This behaviour was changed by adding small amounts of ferrihydrite and aluminium hydroxide to the soil. The aggregation and its stability in the wet state increased, whereas the tensile strength of the dry soils decreased. Both resulted mainly by making the fine particles less mobile by aggregation. In some cases adding large amounts of Al hydroxide caused a resuspension of soil particles by reversing the charge, resulting in instability when wet and high stability when dry. We postulate that aggregation results from the interaction between negatively charged siliceous surfaces and positively charged oxides. In addition we think that newly formed Si–O–Fe bonds play a role in binding particles together as a result of a positive relation between the amount of ferrihydrite added and the oxalate-soluble Si. We conclude that hardsetting resulting from weak structure is caused by lack of aggregating agents such as metal hydroxides.

Notes: Simple models describing nitrogen processes are required both to estimate nitrogen mineralization in field conditions and to predict nitrate leaching at large scales. We have evaluated such a model called LIXIM, which allows calculation of nitrogen mineralization and leaching from bare soils, assuming that these are the dominant processes affecting N in bare soil. LIXIM is a layered, functional model, with a 1-day time step. Input data consist of frequent measurements of water and mineral N contents in soil cores, standard meteorological data and simple soil characteristics. The nitrate transport is simulated using the ‘mixing-cells’ approach. The variations in N mineralization with temperature and moisture are accounted for, providing calculation of the ‘normalized time’. An optimization routine is used to estimate the actual evaporation and the N mineralization rates that provide the best fit between observed and simulated values of water and nitrate contents in all measured soil layers.The model was evaluated in two field experiments (on loamy and chalky soils) including treatments, lasting 9–20 months. The water and nitrate contents in soil were satisfactorily simulated in both sites, and all treatments, including a 15N tracer experiment performed in the loamy soil. In the chalky soil, the calculated water balance agreed well with drainage results obtained in lysimeters and independent estimates of evaporation. At both sites, N mineralization was reduced by the incorporation of crop residues (wheat or oilseed rape straw); the amounts of nitrogen immobilized varied between 20 and 35 kg N ha−1. In the treatments without crop residues, the mineralization rate followed first-order kinetics (against normalized time) in the loamy soil, and zero-order kinetics in the chalky soil. In the latter soil, the mineralization kinetics calculated in situ were close to the kinetics measured in laboratory conditions when both were expressed against normalized time.

Notes: Predictions of changes in soil solution chemistry and exchangeable cations which occur on ash deposition after slash burning are complex and may be facilitated by the use of chemical models. Multi-ion sorption in the topsoil of an Amazonian Acrisol was studied by sequentially adding small amounts of electrolytes to soil and mixtures of soil and ash in batch experiments. A chemical equilibrium model that included inorganic complexation, multiple cation exchange and sparingly soluble salts (aluminium hydroxide and magnesian calcites) was used to interpret the results. The model predicted well the pH and sorption values in all experiments in which there was no addition of ash. The model suggested that cation exchange was the main process determining concentrations of soil solutions in all cases where neutral salt solutions were added, and that proton buffering was achieved by the dissolution of Al(OH)3 which was followed by Al3+ adsorption. Calculation of ion activity products in solutions from various batch experiments in soil + ash mixtures suggested that magnesian calcites of differing solubility may be in equilibrium with the activities of Mg2+ and Ca2+ in solution. An incongruent dissolution of Mg resulted in less soluble magnesian calcites in the ash. The model estimated satisfactorily the pH and the sorption of ions for all experiments with differing ash additions to the soil. Most of the Ca and significant amounts of Mg added in the ash are expected to remain for a long time in the soil and may determine the Ca and Mg status of the soil solution, primarily controlled by principles of solubility products.

Notes: To understand the process and the kinetics of potassium release from the clay interlayer in natural and arable soils in more detail, I tested the hypotheses that large, monovalent cations, especially NH4+ and Cs+, can reduce the release rates of K+ which is exchanged by Ca2+, even if these monovalent cations are present in concentrations of only a few μm. Percolation experiments were carried out with different illitic soil materials, some containing vermiculite, with 5 m m CaCl2 at pH 5.8 and 20°C, in some cases for over 7000 h. NH4+ and Cs+ both caused a large decrease in the rate at which K+ was released, Cs+ especially. Suppression began at 5 μm NH4+ Blocking by 20 μm NH4+ was easily reversible: the release rates readily increased when NH4+ was omitted from the exchange solution. Blocking by 2 μm Cs+ was equal to approximately 90% of that at 10 μm Cs+. Larger concentrations of Cs+ than 10 μm did not further reduce release but rather caused a slight increase, probably because of enhanced exchange of K+ by Cs+ without exfoliation of the interlayer space. Blocking by Cs+ was not reversible within 〉 7000 h of percolation by 5 m m CaCl2. The blocking effect was reproduced in several different soil materials using 10 μm Cs+ but was most pronounced in vermiculite-rich samples. As NH4+ is present in most arable soils, at least in concentrations of a few μm, I conclude that the observed effects are of significance in the K dynamics processes in soils, for example near the roots of plants. Further, very small concentrations of Cs+ in exchange solutions containing a large background of Ca2+ appear to be useful for suppressing K+ release from the interlayer in laboratory studies, probably without significantly altering the exchange at outer mineral surfaces.

Notes: Burning vegetation produces partly charred plant material which subsequently could contribute to the highly refractory proportion of soil organic matter. The presence of charred organic carbon (COC) was investigated in 17 horizons originating from nine soils from Germany and the Netherlands using a suite of complementary methods (high-energy ultraviolet photo-oxidation, scanning electron microscopy, solid-state 13C nuclear magnetic resonance, lignin analysis by CuO oxidation). Charred organic carbon could not be detected in the A horizons of an Alisol and a Gleysol, but it contributed up to 45% of the organic carbon and up to about 8 g kg–1 of the soil in a range of grey to black soils (Cambisol, Luvisol, Phaeozem, Chernozem and Greyzem). All these soils have chernozemic soil properties (dark colour, A–C profile, high base saturation, bioturbation). A 10-km colour sequence of four chernozemic soils, which were very similar in chemical and physical properties, showed a strong relation between colour and the content of COC. This suggests that the COC affects mainly soil colour in the sequence studied. Finely divided COC seems to be a major constituent of many chernozemic soils in Germany. These results suggest that besides climate, vegetation and bioturbation, fire has played an important role in the pedogenesis of chernozemic soils.

Notes: Accessibility of adsorbates to internal active sites of soils depends mainly on the porous structure of the material. We aimed to determine the distribution of pore sizes in soils with varied physico-chemical properties, by combining adsorption of gases and mercury porosimetry. Microporosity was studied by physical adsorption of N2 at 77 K and CO2 at 273 K; mercury intrusion porosimetry allowed us to evaluate the macroporosity; and mesoporosity was determined by capillary condensation of N2 and mercury porosimetry. The soils investigated were essentially macroporous, with volumes between 0.33 and 0.73 cm3 g–1; the maxima in the differential pore-size distribution were in the range 1500–4000 nm. Volumes of meso- or micropores were always less than 10% of macropore volumes. Calculations based on the theory of Dubinin and the αs-method (for N2 at 77 K) provided, generally, coincident results. In a soil containing much organic matter, N2 adsorption was only one-ninth that of CO2 adsorption, showing that N2 adsorption into the narrow micropores of organic matter was kinetically restricted. When accessibility to micropores was not restricted, the total volume of micropores could be deduced from N2 adsorption, whereas CO2 measured exclusively the narrowest microporosity.

Notes: A study was conducted to determine the influence of switchgrass roots on the mobility, adsorption–desorption and mineralization of atrazine in Cullen clay loam and Emporia loamy sand soils. Bromide and atrazine distribution profiles in the leachates indicated greater preferential movement in columns with roots than in columns without roots. Larger concentrations of atrazine were detected at lower depths of Emporia soil with switchgrass roots than without. Adsorption of atrazine was greater in Cullen than in Emporia soil and conformed to Freundlich isotherms. In both Cullen and Emporia soils, adsorption and desorption were not different between soil with or without switchgrass roots. After 84 days of incubation, less than 6% of the applied atrazine was mineralized in the Cullen soil and 2% in Emporia soil. Mineralization was greater in the Cullen soil than in the Emporia soil at 42, 56, 70 and 84 days of sampling. The presence of switchgrass roots did not affect the mineralization of atrazine in Emporia soil. The presence of switchgrass roots caused preferential movement of atrazine, but did not affect its adsorption and mineralization in either soil type.

Notes: Infiltration and sheet flow depend mainly on hydraulic conductivity under nearly saturated conditions which can vary in space and time. We have studied the variation in relation to the topography at different stages of a cultivation year to understand the origin of sheet flow over the land where maize is grown. The sampling was done to a complete factorial design with two landscape positions, upslope and downslope, two sites in relation to wheeling, no and under wheeling, and three dates in relation to soil surface changes and agricultural operations, June and October 1995 and April 1996. Four measurements were randomly distributed on each of the four sites and each of the three dates. The hydraulic conductivity was measured with a tension disc infiltrometer controlling tension at the soil surface. Steady-state infiltration rates were measured at four soil water tensions: 0.1, 0.3, 0.6 and 1.5 kPa. We also measured bulk density and soil moisture content. Hydraulic conductivity decreased during the year at all tensions by a factor of 2–3 where there was no wheeling, and a factor of 4–6 under wheeling. This decrease was greater downslope than upslope, particularly under traffic. The topographic effect was less than that of time but significant. Whereas the hydraulic conductivity was similar over the field at the beginning of the crop cycle (June), it differed according to the landscape position at the end (April), especially under traffic. These topographic variations of the hydraulic conductivity may be due to slight variations of soil components and hydric conditions in the particular case studied. They are important in relation to the characteristics of rainfall and runoff and help explain the spatial and temporal distributions of sheet flow and erosional processes that themselves may contribute to a redistribution of soil particles and induce feedback effects on sheet flow and infiltration.〈section xml:id="abs1-2"〉〈title type="main"〉Variations topographiques et saisonnières de la conductivité hydraulique en condition non saturée mesurée à l’échelle parcellaire par infiltrométrie Résumé La répartition des pluies entre infiltration et ruissellement dépend principalement de la conductivité hydraulique proche de la saturation qui varie à la fois dans l’espace et le temps. Les variations de la conductivité hydraulique selon la topographie sont analysées dans l’objectif de mieux comprendre la distribution spatiale du ruissellement dans le versant. Ces variations spatiales ont été quantifiées durant une année culturale sur une parcelle de maïs, en relation avec la topographie de la parcelle et la succession des opérations culturales. Le plan d’échantillonnage est un plan factoriel complet prenant en compte la position topographique, amont et aval, la localisation relative aux traces de roue, avec et sans trace de roue, et trois dates relatives à la dynamique des états de surface et aux opérations culturales, juin et octobre 1995 et avril 1996. Quatre mesures de la conductivité hydraulique ont été réalisées sur chaque site avec un infiltromètre à succion appliquéà la surface du sol. Les flux infiltrés ont été mesurés à l’équilibre et à quatre potentiels: –0.1, –0.3, –0.6 et –1.5 kPa. Ces mesures ont été complétées par un suivides teneurs en eau et des densités apparentes sur l’ensemble du versant et de la période d’étude. La conductivité hydraulique décroît au cours de l’année quel que soit le potentiel. Elle est divisée par un facteur deux à trois durant toute la période d’étude sans trace de roue, par un facteur quatre à six avec trace de roue. Cette diminution de la conductivité hydraulique dépend de la position topographique: elle est plus forte en aval qu’en amont, en particulier sur passage de roue. Tandis que la conductivité hydraulique apparaît relativement homogène sur le champ en début de cycle cultural, elle varie fortement de l’amont à l’aval en fin de cycle cultural, en particulier sur trace de roue. Ces variations spatiales de la conductivité hydraulique semblent liées à de faibles variations des constituants du sol qui interagissent avec des conditions hydriques très différentes. Ces variations de la conductivité hydraulique sont importantes au regard des caractéristiques des pluies et du ruissellement sur ce site. Elles expliquent pour partie la distribution spatiale et temporelle du ruissellement et des processus d’érosion qui eux-même peuvent contribuer à la redistribution des constituants dans le versant et induire des effets retour sur les processus de ruissellement et d’infiltration.

Notes: Heavy density fractions of soil contain organic matter tightly bound to the surface of soil minerals. The chemical composition and ecological meaning of non-metabolic decomposition products and microbial metabolites in organic–mineral bonds is poorly understood. Therefore, we investigated the heavy fraction (density 〉 2 g cm–3) from the topsoil of a Gleysol (Bainsville, Ottawa, Canada). It accounted for 952 g kg–1 of soil and contained 19 g kg–1 of organic C. Pyrolysis-field ionization mass spectra showed intensive signals of carbohydrates, and phenols and lignin monomers, alkylaromatics (mostly aromatic) N-containing compounds, and peptides. These classes of compound have been proposed as structural building blocks of soil organic matter. In comparison, the light fraction (density 〉 2 g cm–3) was richer in lignin dimers, lipids, sterols, suberin and fatty acids which clearly indicate residues of plants and biota. To confirm the composition and stability of mineral-bound organic matter, we also investigated the heavy fraction (density 〉 2.2 g cm–3) from clay-, silt- and sand-sized separates of the topsoil of a Chernozem (Bad Lauchstädt, Germany). These heavy size separates differed in their mass spectra but were generally characterized by volatilization maxima of alkylaromatics, lipids and sterols at about 500°C. We think that the observed high-temperature volatilization of these structural building blocks of soil organic matter is indicative of the organic–mineral bonds. Some unexpected low-temperature volatilization of carbohydrates, N-containing compounds, peptides, and phenols and lignin monomers was assigned to hot-water-extractable organic matter which accounted for 7–27% of the carbon and nitrogen in the heavy fractions. As this material is known to be mineralizable, our study indicates that these constituents of the heavy density fractions are degradable by micro-organisms and involved in the turnover of soil organic matter.

Notes: A wavelet is a compact analysing kernel that can be moved over a sequence of data to measure variation locally. There are several families of wavelet, and within any one family wavelets of different lengths and therefore smoothness and their corresponding scaling functions can be assembled into a collection of orthogonal functions. Such an assemblage can then be applied to filter spatial data into a series of independent components at varying scales in a single coherent analysis. The application requires no assumptions other than that of finite variance. The methods have been developed for processing signals and remote imagery in which data are abundant, and they need modification for data from field sampling. The paper describes the theory of wavelets. It introduces the pyramid algorithm for multiresolution analysis and shows how it can be adapted for fairly small sets of transect data such as one might obtain in soil survey. It then illustrates the application using Daubechies’s wavelets to two soil transects, one of gilgai on plain land in Australia and the other across a sedimentary sequence in England. In both examples the technique revealed strongly contrasting local features of the variation that had been lost by averaging in previous analyses and expressed them quantitatively in combinations of both scale and magnitude. Further, the results could be explained as the spatial effects of change in topography or geology underlying the variation in the soil.〈section xml:id="abs1-2"〉〈title type="main"〉Analyse et éclairissement sur la variation du sol en utilisant les ondelettes RésuméUne ondelette est un noyau compact d’analyse qu’on peut passer sur une séquence de données pour quantifier la variation localement. Plusieurs familles d’ondelettes existent. Chaque famille est caractérisée par des fonctions d’échelle de longueurs d’ondes et de degrés de lissage différents, le tout constituant un ensemble de fonctions orthogonales. L’application de ces fonctions sur des données spatiales est une méthode d’analyse unique et cohérente qui permet de filtrer des données spatiales en identifiant des composantes indépendentes à différentes échelles. L’application n’éxige que l’hypothése d’une variance bornée. Ces méthodes ont été developpées pour le traitement de signaux qui contiennent énormement de données comme des images de télédétection. Elles necessitent des modifications avant de les appliquer sur des données d’échantillonage provenant du terrain. Cet article décrit la théorie d’ondelettes. Il introduit l’algorithme pyramidal pour une analyse à quelques résolutions. Il montre ensuite comment on peut adapter l’algorithme pour un ensemble de données peu nombreuses comme celles qu’on peut obtenir lors d’une prospection du sol sur des transects.Les resultats montrent l’application des ondelettes de Daubechies à deux transects pédologiques, lepremier situé sur une plaine marquée de gilgaï en Australie, et le second à travers une séquence de sediments jurassiques en Angleterre. Dans les deux exemples, la technique a révélé des contrastes bien marqués de certain traits locaux qui ont été totalement occultés dans des analyses antérieures basées sur le calcul des moyennes. La technique des ondelettes permet aussi une quantification de ces traits locaux, en fonction de l’échelle d’observation et de la magnitude. De plus, ces résultats peuvent être interprétés par les effets d’une variation spatiale de la topographie ou de la géologie qui sont à l’origine de la variation du sol.

Notes: Data from two Podzol O and E horizons, sampled in 1-cm layers at 13 points within 2 m × 2 m plots, were used to test the hypothesis that the composition of hydrogen ions (H) and aluminium (Al) adsorbed to the solid-phase soil organic matter (SOM) determines pH and Al solubility in organic-rich acidic forest soils. Organically adsorbed Al was extracted sequentially with 0.5 m CuCl2 and organically adsorbed H was determined as the difference between total acidity titrated to pH 8.2 and Al extracted in 0.5 m CuCl2. The quotient between fractions of SOM sites binding Al and H (NAl/NH) is shown to determine the variation in pH and Al solubility. It is furthermore shown that models in which pH and Al solubility are linked via a pH-dependent solubility of an Al hydroxide and in which cation exchange between Al3+ and Ca2+, rather than cation exchange between Al3+ and H+, is the main pH-buffering process cannot be used to simulate pH or Al solubility in O and E horizons. The fraction of SOM sites adsorbing Al increased by depth in the lower O and throughout the E horizon at the same magnitude as sites adsorbing H decreased. The fraction of sites binding the cations Ca2+ + Mg2+ + K+ + Na+ remained constant. It is suggested that a net reaction between Al silicates (proton acceptors) and protonated functional groups in SOM (proton donors) is the long-term chemical process determining the composition of organically adsorbed H and Al in the lower part of the O and in the E horizon of Podzols. Thus, in the long term, pH and Al solubility are determined by the interaction between organic acidity and Al alkalinity.

Notes: A method is presented for the prediction and attenuation of long-period water-layer multiples based on an adaptive multichannel lattice algorithm. The multichannel approach can be viewed as a generalization of two-dimensional linear prediction. The multichannel linear least-squares prediction problem is reviewed briefly and the performance and limitations of the algorithm are demonstrated on two different marine data sets with different properties of the simple and peg-leg multiple system. The algorithm works well even on problematic data sets and is very easy to apply.

Notes: Gravity observations made over the Groningen gasfield from 1978 onwards are available. The existing gravity data were reanalysed and a new survey was carried out. The observed gravity changes obtained during four surveys, spanning an 18-year period, were compared with the gravity effect due to natural gas extraction computed from the reservoir model and the production data. The random error in the gravity values is small enough to detect the effect of gas extraction after a few years (up to 2 μgal/year). The trends obtained from observed and calculated gravity changes agree well within their expected error margins after statistical data snooping. Due to the inadequate measurement set-up of the previous surveys, systematic errors present in the gravity data hampered the use of this data for refinement of the reservoir and the production models. However, with improved gravity monitoring, in particular a very systematic network set-up and well-calibrated instruments, this should now be feasible. By using the network analysis from the 1996 measurements for the planning of a network in space and time, most of the current problems can be avoided.

Notes: We analyse the geophone orientation and coupling in a data set from the North Sea. Based on the polarization of the water-break on the sea-floor, we have derived processing algorithms for determining the receiver orientation for gimballed and non-gimballed geophone systems. For a gimballed system, the problem reduces to a simple horizontal rotation. However, for a non-gimballed system, where all three geophone axes may vary due to varying acquisition conditions such as dipping sea-floor, twisting of recording cable, etc., the three orientation angles cannot be found directly from the recorded displacement vectors. Using the data redundancy within a common-receiver gather, a robust two-stage method is derived for the non-gimballed system in which all three orientations can initially be unknown. Testing on the North Sea data set acquired with a gimballed system shows that the three-component geophones in the data set are orientated satisfactorily within an error of 5°. However, there are some undesirable cross-couplings between the vertical and horizontal geophones, which results in leakage of shear-wave energy from the horizontal components to the vertical components.

Notes: The great majority of data processing methods for magnetotelluric measurements are based on an average Fourier spectral analysis to estimate either the transfer function or the coherence function. This assumes that the time-series data are stationary. The aim of this paper is to present a preprocessing method able to extract the stationary and coherent part of the original signals. The practitioner can then apply the magnetotelluric method of his choice to these new data. This preliminary data sorting is done in four steps: (i) slow drifts are eliminated using a high-pass filter whose cut-off frequency is determined by an iterative procedure. Based on run tests, this procedure is also able to remove segments with non-independent samples in the time series; (ii) non-stationary segments are eliminated after band-pass filtering; (iii) non-coherent segments are eliminated before spectral analysis; (iv) the impedance tensor value is then retained, at a given frequency, only if the signals are coherent. This preprocessing method was tested on the simplest, but still used, magnetotelluric method which uses only two field components, and it was found that the average resistivity standard deviation decreased significantly from 14.6 Ωm without sorting to 8.6 Ωm after sorting.

Notes: Prestack depth migrations based on wavefield extrapolation may be computationally expensive, especially in 3D. They are also very dependent on the acquisition geometry and are not flexible regarding the geometry of the imaging zone. Moreover, they do not deal with all types of wave, considering only primary reflection events through the model. Integral approaches using precalculated Green's functions, such as Kirchhoff migration and Born-based imaging, may overcome these problems. In the present paper, both finite-difference traveltimes and wavefront construction are used to obtain asymptotic Green's functions, and a generalized diffraction tomography is applied as an example of Born-based acoustic imaging. Target-orientated imaging is easy to perform, from any type of survey and subselection of shot/receiver pairs. Multifield imaging is possible using Green's functions that take into account, for instance, reflections at model boundaries. This may help to recover parts of complex structures which would be missing using a paraxial wave equation approach. Finally, a numerical evaluation of the resolution, or point-spread, function at any point of the depth-migrated section provides valuable information, either at the survey planning stage or for the interpretation.

Notes: We propose a new method for removing sea-surface multiples from marine seismic reflection data in which, in essence, the reflection response of the earth, referred to a plane just above the sea-floor, is computed as the ratio of the plane-wave components of the upgoing wave and the downgoing wave. Using source measurements of the wavefield made during data acquisition, three problems associated with earlier work are solved: (i) the method accommodates source arrays, rather than point sources; (ii) the incident field is removed without simultaneously removing part of the scattered field; and (iii) the minimum-energy criterion to find a wavelet is eliminated.Pressure measurements are made in a horizontal plane in the water. The source can be a conventional array of airguns, but must have both in-line and cross-line symmetry, and its wavefield must be measured and be repeatable from shot to shot. The problem is formulated for multiple shots in a two-dimensional configuration for each receiver, and for multiple receivers in a two-dimensional configuration for each shot. The scattered field is obtained from the measurements by subtracting the incident field, known from measurements at the source. The scattered field response to a single incident plane wave at a single receiver is obtained by transforming the common-receiver gather to the frequency–wavenumber domain, and a single component of this response is obtained by Fourier transforming over all receiver coordinates. Each scattered field component is separated into an upgoing wave and a downgoing wave using the zero-pressure condition at the water-surface. The upgoing wave may then be expressed as a reflection coefficient multiplied by the incident downgoing wave plus a sum of scattered downgoing plane waves, each multiplied by the corresponding reflection coefficient. Keeping the upgoing scattered wave fixed, and using all possible incident plane waves for a given frequency, yields a set of linear simultaneous equations for the reflection coefficients which are solved for each plane wave and for each frequency. To create the shot records that would have been measured if the sea-surface had been absent, each reflection coefficient is multiplied by complex amplitude and phase factors, for source and receiver terms, before the five-dimensional Fourier transformation back to the space–time domain.

Notes: A new algorithm for tomographic inversion of traveltimes of reflected and refracted seismic waves is developed. The inversion gives interface configurations and velocity distributions in layers. The important features of the algorithm are: (a) the inclusion of shot time delays in the list of unknown parameters; (b) the regularization is applied in such a way that the most probable model is characterized by the similarity of neighbouring interfaces. As the problem under consideration is non-linear, several iterations are necessary in order to obtain the final model. In the case of a very inexact initial model, a ‘layer-by-layer’ inversion strategy is recommended as a first inversion step. The inversion program is supplied with a user interface, thanks to which one can: (a) pick interactively and identify seismic traveltimes; (b) build and edit depth/velocity models; and (c) display calculated traveltime curves and compare them with picked traveltimes as well as with the original seismic sections. The efficiency of the inversion software developed is illustrated by a numerical example and a field example in which shallow seismic data are considered. Application to wide-aperture reflection/refraction profiling (WARRP) data is also possible.

Notes: Non-aliased integral (Kirchhoff-type) transformations for forward and downward wavefield extrapolations in inhomogeneous media with interfaces are described. Special weights are computed to compensate for operator aliasing and finite-aperture effects, even when the data are spatially aliased and irregularly sampled. Basic components of the algorithm, such as Green's function computation, can be replaced by alternative solutions in conjunction with ray-tracing methods. Applications of this algorithm to model and real data in both two and three dimensions are discussed in terms of its impact on seismic modelling, multiple prediction and prestack imaging.

Notes: Attenuation of seismic waves, quantified by the seismic quality factor Q, holds important information for seismic interpretation, due to its sensitivity to rock and fluid properties. A recently published study of Q, based on surface seismic reflection data, used a modified spectral ratio approach (QVO), but both source and receiver responses were treated as isotropic, based on simple raypath arguments. Here, this assumption has been tested by computing apparent attenuation generated by frequency-dependent directivity of typical marine source and receiver arrays and acquisition geometries. Synthetic wavelet spectra were computed for reflected rays, summed over the first Fresnel zone, from the base of a single interval, 50–3000 m thick and velocity 2000 m/s, overlying a 2200 m/s half-space, and for offsets of 71–2071 m. The source and receiver geometry were those of an actual survey.The modelled spectra are clearly affected by directivity, most strongly because of surface ghosts. In general, the strong high-frequency component, produced by the array design, leads to apparently negative attenuation in individual reflection events, though this is dependent on offset and target depth. For shallow targets (less than 400–500 ms two-way traveltime (TWT) depth), apparent Q-values as extreme as −50 to −100 were obtained. For deeper target depths, the directivity effect is far smaller. The implications of the model study were tested on real data. QVO was applied to 20 true-spectrum-processed CMPs, in a shallow (405–730 ms TWT) and a deeper (1000–1300 ms TWT) interval, firstly using a measured far-field source signature (effectively isotropic), and secondly using computed directivity effects instead. Mean interval Q−1-values for the deeper interval, 0.029 ± 0.011 and 0.027 ± 0.018 for conventional and directional processing, respectively, suggested no directivity influence on attenuation estimation. For the shallow interval (despite poor spectral signal-to-noise ratios and hence scattered attenuation estimates), directional processing removed directivity-generated irregularities from the spectral ratios, resulting in an improvement from Q−1int = −0.036 ± 0.130 to a realistic Q−1int = 0.012 ± 0.030: different at 94% confidence level. Equivalent Q-values are: for the deeper interval, 35 and 37 for conventional and directional processing, respectively, and −28 and 86 for the shallow interval.These results support the conclusions of the model studies, i.e. that source/receiver directivity has a negligible effect except for shallow targets (e.g. TWT depth ≤ 500 ms) imaged with conventional acquisition geometry. In such cases directivity corrections to spectra are strongly recommended.

Notes: Modelling of AVO signatures for reservoir characterization requires VS estimation from other available logs when shear-wave data are not available. We tested various models for predicting VS from P-wave velocity, porosity and shale volume measured in well logs. Effective medium models which characterize the pore space in terms of ellipsoidal inclusions were compared with statistical VP–VS regressions. The inclusion models were calibrated by non-linear minimization of the difference between model-predicted velocities and actual measured velocities. The quality of the VS prediction was quantified in terms of the rms error by comparison with shear-wave data in wells where both VP and VS were measured. The linear regressions were found to be more robust and the rms error in the prediction was comparable to effective medium model-based predictions.

Notes: Multichannel filtering to obtain wavefield separation has been used in seismic processing for decades and has become an essential component in VSP and cross-well reflection imaging. The need for good multichannel wavefield separation filters is acute in borehole seismic imaging techniques such as VSP and cross-well reflection imaging, where strong interfering arrivals such as tube waves, shear conversions, multiples, direct arrivals and guided waves can overlap temporally with desired arrivals. We investigate the effects of preprocessing (alignment and equalization) on the quality of cross-well reflection imaging wavefield separation and we show that the choice of the multichannel filter and filter parameters is critical to the wavefield separation of cross-well data (median filters, f–k pie-slice filters, eigenvector filters). We show that spatial aliasing creates situations where the application of purely spatial filters (median filters) will create notches in the frequency spectrum of the desired reflection arrival. Eigenvector filters allow us to work past the limits of aliasing, but these kinds of filter are strongly dependent on the ratio of undesired to desired signal amplitude. On the basis of these observations, we developed a new type of multichannel filter that combined the best characteristics of spatial filters and eigenvector filters. We call this filter a ‘constrained eigenvector filter’. We use two real data sets of cross-well seismic experiments with small and large well spacing to evaluate the effects of these factors on the quality of cross-well wavefield separation. We apply median filters, f–k pie-slice filters and constrained eigenvector filters in multiple domains available for these data sets (common-source, common-receiver, common-offset and common-midpoint gathers). We show that the results of applying the constrained eigenvector filter to the entire cross-well data set are superior to both the spatial and standard eigenvector filter results.

Notes: For many geophysical and astrophysical applications the relationship between physical bulk properties and porosity is needed. As part of the preparation for the ‘Rosetta’ comet rendezvous mission, a simple model for the porosity dependence of the elastic properties of granular media, i.e. the elastic moduli and the propagation velocities of elastic waves, has been developed based on textural properties and the contact stiffness of the constituent particles. It is shown that the derived relationships fit very well with sandstone data. The model is also consistent with data for snow and ice and is in agreement with the transformation mechanisms from snow to ice. A short review shows the relevance to another physical bulk property, the thermal conductivity.

Notes: Seismic velocities (Vp and Vs) of compressional (P-) and shear (S-) waves are important parameters for the characterization of marine sediments with respect to their sedimentological and geotechnical properties. P- and S-wave velocity data of near-surface marine sediments (upper 9 m) of the continental slope of the Barents Sea are analysed and correlated to sedimentological and geotechnical properties. The results show that the S-wave velocity is much more sensitive to changes in lithology and mechanical properties than the P-wave velocity, which is characterized by a narrow range of values. The correlation coefficients between S-wave velocity and silt and clay content, wet bulk density, porosity, water content and shear strength are higher than 0.5 while the correlation coefficients of P-wave velocity and the same parameters are always lower than 0.4. Although the relationship between Vs and clay content has been widely described, the data show that Vs is better correlated with silt content than with clay content for the sediments of the area investigated. However, they show different trends. While Vs increases with increasing clay content, it decreases with increasing silt content.

Notes: In order to perform resistivity imaging, seismic waveform tomography or sensitivity analysis of geophysical data, the Fréchet derivatives, and even the second derivatives of the data with respect to the model parameters, may be required. We develop a practical method to compute the relevant derivatives for 2.5D resistivity and 2.5D frequency-domain acoustic velocity inversion. Both geophysical inversions entail the solution of a 2.5D Helmholtz equation. First, using differential calculus and the Green's functions of the 2.5D Helmholtz equation, we strictly formulate the explicit expressions for the Fréchet and second derivatives, then apply the finite-element method to approximate the Green's functions of an arbitrary medium. Finally, we calculate the derivatives using the expressions and the numerical solutions of the Green's functions. Two model parametrization approaches, constant-point and constant-block, are suggested and the computational efficiencies are compared. Numerical examples of the derivatives for various electrode arrays in cross-hole resistivity imaging and for cross-hole seismic surveying are demonstrated. Two synthetic experiments of resistivity and acoustic velocity imaging are used to illustrate the method.

Notes: Seismic exploration in mountainous areas imposes serious compromises on both acquisition and processing. Access restrictions usually result in profiles that are not straight and are not recorded along the true dip direction (if there is a true dip direction!). Processing constraints often result in very poor approximate corrections for elevations and for deviations from a straight line. Most fundamentally, 2D acquisition and processing assumes that the earth is 2D; this assumption is often seriously violated in mountainous areas. While we cannot efficiently correct 2D seismic data for the effects of a fully 3D subsurface, we can improve the data quality in thrust areas where the assumption of 2D subsurface variation is reasonable. We do this in a series of small steps, which improves the accuracy of several approximations made in processing 2D land data.

Notes: During processing of magnetotelluric (MT) data, acquired in a survey carried out in southern Italy, a problem was encountered, connected with the so-called ‘dead band’ of the MT signal (around 1 Hz). In the apparent resistivity curves of some MT soundings, a V-shaped minimum appeared, centred on the dead-band frequency. This phenomenon turned out to originate from a strong artificial source and was not due to a downward bias of the robust processing techniques adopted. The source distance from the MT sounding locations was such that the V-shaped minimum fell precisely in the dead-band frequency range. Theoretical considerations about fields generated by an electric dipole led us to the probable identification of the source as the d.c.-powered railway between Naples and Bari.

Notes: Analysis of prestack P-wave seismic data yields information about both the P- and S wave properties of the earth. An anticipated advantage of having two measurements (P and S) is that they can be combined into a new measurement that is less sensitive to lithology variations and more sensitive to fluid effects. The amplitude-variation-with-offset (AVO) gradient is one such measure that is often used qualitatively as a fluid indicator. The gradient always becomes softer (more negative) when hydrocarbon replaces brine in the pore spaces but the overall AVO response is dominated by the lithology. Fluid effects are expressed primarily by the normal-incidence P-wave response and only secondarily by the offset dependence. The gradient often does not function as an effective fluid indicator. This is partially due to the fact that the gradient is roughly twice as sensitive to S- than to P-wave properties. More importantly, effective random noise in the CMP gathers introduces a strong correlation between the AVO intercept and gradient and, hence, between the measured P- and S-wave properties. This correlation in the AVO attributes corresponds to a significant error in the estimation of the S-wave properties and can dominate the measurements from many of the popular AVO techniques. A simple method to minimize the effect of this noise-induced correlation is to stack the data. The stack corresponds to a coordinate rotation in elastic space with the stack amplitudes measured along one of the new axes and the other (unmeasured) axis naturally tending to line up with the noise and thus suppressing it. Fluid effects cause the data to move roughly perpendicular to this noise trend. The stack axis is then in the direction of the fluid effect. The stack thus combines both the P- and S-wave (normal and oblique incidence) information into a single measurement which can be made to optimally suppress background noise and highlight fluid effects. A major consequence of this interpretation is the simplicity of both prospect identification and quantitative amplitude analysis.

Notes: This paper describes a large-scale reservoir characterization experiment carried out in Oman in 1991 which comprised the acquisition, processing and interpretation of a 28.4 km2 3D multicomponent seismic experiment over the Natih field. The objective of the survey was to obtain information on the fracture network present in the Natih carbonates from shear-wave anisotropy. Shear-wave anisotropy in excess of 20% time splitting was encountered over a large part of the survey. The seismic results are confirmed by geological and well data but provide additional qualitative information on fracturing where this was not available before. Regions of stronger and weaker shear-wave anisotropy appear to be fault-bounded. The average well flow rates (which are fracture-dominated) within such blocks correlate with the average anisotropy of the blocks. The further observation that the anisotropy is largest in the fracture gas cap of the reservoir suggests that shear waves can provide a direct hydrocarbon indicator for fractured rock.

Notes: It is shown that the electromagnetic (EM) spectrum is characterized by strong amplitude-modulated transmitters operating in the target bandwidth of transient electromagnetic (TEM) measurements. As these transmitters cause significant noise in TEM soundings, it is mandatory to band-limit the input signals to improve the signal-to-noise ratio and thereby the depth of exploration. Band-limitation will distort the TEM responses, which leads to erroneous inversion results if the applied low-pass filters are not accounted for in the inversion scheme. We incorporate the low-pass filters in the inversion scheme and test the inversion approach on theoretical and field data. Inversion of band-limited theoretical responses results in recovery of erroneous resistivity models if the filters are not included in the inversion scheme. By contrast, inversion of band-limited theoretical and field data, for which the applied low-pass filters are included in the inversion scheme, leads to recovery of similar resistivity models, independent of the applied cut-off frequencies.

Notes: Non-linear least-squares inversion operates iteratively by updating the model parameters in each step by a correction vector which is the solution of a set of normal equations. Inversion of geoelectrical data is an ill-posed problem. This and the ensuing suboptimality restrict the initial model to being in the near vicinity of the true model. The problem may be reduced by introducing damping into the system of equations. It is shown that an appropriate choice of the damping parameter obtained adaptively and the use of a conjugate-gradient algorithm to solve the normal equations make the 1D inversion scheme efficient and robust. The scheme uses an optimal damping parameter that is dependent on the noise in the data, in each iterative step. The changes in the damping and relative residual error with iteration number are illustrated. A comparison of its efficacy over the conventional Marquardt and simulated annealing methods, tested on Inman's model, is made. Inversion of induced polarization (IP) sounding is obtained by inverting twice (true and modified) DC apparent resistivity data. The inversion of IP data presented here is generic and can be applied to any of the IP observables, such as chargeability, frequency effect, phase, etc., as long as these observables are explicitly related to the DC apparent resistivity. The scheme is used successfully in inverting noise-free and noisy synthetic data and field data taken from the published literature.

Notes: The aim of this study was to investigate whether bovine faeces could replace rumen liquor collected from fistulated sheep in the determination of organic matter digestibility (OMD) of forages in vitro. The first experiment compared two inocula that had been prepared by mixing 83 and 333 g of cattle faeces with 1 litre of artificial saliva according to the first stage of the in vitro technique of Tilley and Terry (Journal of the British Grassland Society,18, 104–111, 1963). Inocula made from 333 g of faeces gave higher values of OMD than that prepared from 83 g of faeces and was used in subsequent studies. When a 48-h acid pepsin digestion, the second stage of the Tilley and Terry (1963) technique, was included, the OMD values and the ease of filtering of undigested residues were increased. The second experiment compared digestibility determined with inocula produced from faeces with that produced from rumen liquor. The OMD of a ryegrass hay determined with an inoculum prepared from faeces was 0·468 compared with 0·528 when using an inoculum prepared from rumen liquor. The third experiment investigated the accuracy of estimating the OMD of eight contrasting forages determined using rumen liquor from three sheep (y) and faeces from two cows (x). All regressions were significant (P 〈 0·001) and had residual standard deviations (r.s.d.) of between ±0·019 and ±0·022. The final study compared faeces collected from cattle fed with either hay or concentrates. The OMD of forage was higher from the hay-fed animals, but the r.s.d. values of the regressions were similar. It was concluded that cattle faeces shows potential as an alternative to rumen liquor collected from rumen-fistulated sheep for use in the in vitro digestibility assay of forages.

Notes: The effect of botanical composition and nutrient availability on the relative allocation of biomass to stems and leaves in a permanent upland pasture in the central Pyrenees was assessed. Six short-term and medium-term fertilizer treatments (nitrogen and phosphorus) were applied to a meadow to create large differences in the proportions of the different species and a wide range of herbage nutrient status. The above-ground herbage dry-matter components (green leaves, sheaths and stems for grasses and dicotyledons, and senescent material) were measured for the first growth cycle.The leaf mass depended mainly on the herbage N status, whereas the stem mass depended both on the botanical composition and on the herbage N and P status. During spring growth, the proportion of leaves in the above-ground dry matter decreased faster in plots that had the highest nutrient status or that were composed of species characteristic of nutrient-rich (Festuca-poor) habitats. Application of fertilizer decreased the proportion of leaves both for grasses and dicotyledons, but to a larger extent for grasses. There was a single relationship between the proportion of leaves in the above-ground dry matter and the total mass of above-ground dry matter, irrespective of the sampling date, the botanical composition or the herbage nutrient status. This statistical relationship resulted from (i) a faster increase in stem mass than leaf mass for a given botanical composition when the herbage nutritional status increased, (ii) a greater stem mass when the sward was composed of species usually found in nutrient-rich habitats.

Notes: The objective of this study was to investigate the relationship between level of sward height reduction (SHR) and short-term intake of herbage by lactating dairy cows offered swards differing in initial bulk density (BD). Three experiments were carried out in which cows were presented with swards representing different levels of SHR (nominally described as ungrazed, low, moderate and high). Experiments 1 and 2 differed with respect to initial sward BD [ungrazed sward 1·7 vs. 2·5 kg dry matter (DM) m−3 respectively]. Experiment 3 investigated the interaction between BD and SHR.In each experiment, sixteen Holstein/Friesian cows (fitted with excreta collection bags) grazed for a 1-h period in 200-m2 plots that had been grazed on the previous day to predetermined sward surface heights (SSH) (= levels of SHR). Herbage intake rates were assessed from changes in live weight, with a correction for insensible weight loss (IWL). Biting rates were recorded from visual observation. High levels of SHR were associated with a significant reduction in SSH, herbage mass and leaf fraction, and a significant increase in sward BD, stem and dead fraction, and DM concentration.Herbage intake, expressed either as DM intake per bite or DM intake per hour, declined as level of SHR increased from low to high. The level of SHR generally had no effect on biting rate. Intake rates varied from 1·9 to 4·4 kg DM h−1, whereas DM intake bite−1 ranged from 0·5 to 1·3 g. Pooled regression analysis identified SSH (P 〈 0·001; r2 = 0·94) as the principal determinant of DM intake bite−1. The regression equation was not significantly improved by the addition of terms for leaf fraction, BD, or herbage mass. In Experiment 3, a significant interaction between level of SHR and sward BD was observed. It is concluded that the principal factor controlling intake (g DM bite−1 or kg DM h−1), as swards are progressively grazed down, is SSH, but at a high level of SHR, sward BD also influences intake bite−1.

Notes: The severity of damage by certain invertebrate species to white clover (Trifolium repens) leaves on the main stolons of plants grown in mixed perennial ryegrass/white clover swards was examined in a field experiment in the Netherlands in which two cutting frequencies (high and low) and three white clover cultivars (Retor, Alice and Gwenda) were used. The damage to the leaves was described in terms of the numbers of damaged leaves and the extent of that damage (slight 〈20%, moderate 20–50% and heavy 〉50%). The relationships between leaf damage and sward characteristics (white clover content, above-ground biomass and sward height) were evaluated throughout the growing season.Over the whole experimental period, 23·7% and 27·4% of the total number of leaves produced per stolon were damaged by slugs and weevils in the low- and high-frequency cutting treatments respectively. High-frequency cutting increased the number of leaves in the total leaf damage and moderate leaf damage categories by 21·4% and 34·8%, respectively, compared with the low-frequency cutting. The cv. Retor (medium-leaved) experienced the most severe damage by invertebrates. It had much higher leaf damage than cvs Alice (large-leaved) and Gwenda (small-leaved) at either cutting frequency, both in the total number of damaged leaves and in the different damage categories. Differences among cultivars in the number of damaged leaves and relative leaf damage occurred primarily in spring, late summer and autumn, but did not differ during the early- and mid-summer months.This study indicates that variations in leaf damage among clover cultivars were associated with differences in measured sward characteristics. Both the number of damaged leaves and the relative leaf damage were strongly negatively correlated with white clover content and biomass in spring, late summer and autumn under each cutting treatment. White clover content and biomass explained 65%, 59% and 50% of the variation in the number of damaged leaves in spring, late summer and autumn, respectively, and 58%, 57% and 45% of the variation in relative leaf damage in these three periods. Thus, sward characteristics may play a role in regulating the severity of invertebrate damage to clover leaves in addition to the primary effects of HCN.

Notes: In a field experiment carried out over 3 years, the nitrate content of herbage from perennial ryegrass (Lolium perenne) swards increased exponentially with nitrogen application rate, but herbage nitrate content appeared to reach potentially dangerous concentrations only when nitrogen application rates were greater than those needed to stimulate dry-matter production. Thus, on average over all the harvests, maximum yield could be obtained with annual application rates of 400 kg N ha–1 (six applications of 67 kg N ha–1) for perennial ryegrass and 300 kg N ha–1 (six applications of 50 kg N ha–1) for perennial ryegrass/white clover (Trifolium repens) swards, whereas the mean nitrate concentrations were 3340 and 2929 mg NO3 kg–1 dry matter (DM) respectively. Nitrate content, however, varied considerably from harvest to harvest, reaching maxima of 9345 mg NO3 kg–1 DM at 400 kg N ha–1 for perennial ryegrass and 6255 mg NO3 kg–1 DM at 300 kg N ha–1 for perennial ryegrass/white clover. The nitrate content of herbage from perennial ryegrass/white clover swards was always greater than that of perennial ryegrass swards receiving the same rate of nitrogen application, even though in the herbage from the mixed sward the nitrate content of white clover was usually less than half that of the perennial ryegrass component. The physical environment did not have a clearly interpretable effect on nitrate content, although herbage harvested in May had a much lower nitrate content than that harvested at any other time of the season. It was not possible to find a single multiple regression equation relating herbage nitrate content to nitrogen application and to other environmental variables that explained more than 60% of the variance in herbage nitrate, but it is suggested that, by reducing the later-season nitrogen applications from 67 to 50 and finally to 33 kg N ha–1 for perennial ryegrass and from 50 to 33 kg N ha–1 for perennial ryegrass/white clover, it would be possible to achieve over 90% of the maximum yield while reducing average nitrate content to 〈40% of that at maximum yield, with no samples containing more than 2300 mg NO3 kg –1 DM.

Notes: The effect of management of maize residues on the population of insects inimical to the establishment of alfalfa, Medicago sativa L., was examined 1 month after sowing. Alfalfa was sown in early and late April, and late May for 3 consecutive years in Centre County, Pennsylvania, USA. The sowings were made by conservation tillage (direct drilled into maize residue with minimal disturbance of soil) into three different maize residues after (1) silage, (2) silage-plus-winter rye cover and (3) grain. A fourth sowing was by conventional tillage (spring ploughing and harrowing maize residue) after grain. A pesticide, carbofuran (granular formulation) was applied at sowing to half of the plots. Methiocarb bait, a molluscicide, was applied as a split application to the same plots, one-half at sowing and the remainder 2 weeks later. Alfalfa plots in the silage-plus-rye maize residues were colonized with fewer insects than the other residue treatments. Excessive growth of rye in early spring that was not successfully suppressed by herbicide treatment produced vigorous rye plants and fewer alfalfa seedlings. Consequently, silage plus rye had the lowest yield of alfalfa in early April sowings in 2 out of 3 years. Insects known to feed on alfalfa, such as tarnished plant bugs, Lygus lineolaris (Palsot de Beauvois), the plant bug, Plagiognathus politus Uhler, pea aphids, Acyrthosiphon pisum (Harris), and meadow spittlebugs, Philaenus leucophthalmus spumarius (L.), were collected by sweep net 1 month after sowing and less often in the silage-plus-rye treatment. These species were present in greater numbers in the other maize residues that had significantly more alfalfa forage.The insidious flower bug, Orius insidious (Say), and a damsel bug, Nabis americoferus Carayon, were collected in significantly greater numbers in the early April sowings, which corresponded with the peak populations of pea aphid. The potato leafhopper, Empoasca fabae (Harris), was most abundant in plots sowed in late April or late May.Pesticide treatment applied at the time of sowing had very little effect on numbers of insects collected by sweep net 1 month later in 1986 and 1988. However, pesticide treatment significantly increased yield by 280 kg ha−1 in 1987, even though sweep net collections of insects were not reduced by the pesticide treatment. Therefore, the beneficial effect of the pesticide could not be explained on the basis of the insects collected.The highest yields of alfalfa were obtained from the early April sowing into maize residues. This coincided with the time when the majority of alfalfa pests were less abundant than in later sowings; fewer pests were found on the sowings into silage-plus-rye residue. Also, when the rye forage yield was combined with the alfalfa yield, this became the most productive system.

Notes: Quantitative traits and allozymes were used in two experiments to clarify the ecotypic differentiation of natural cocksfoot (Dactylis glomerata) populations from north-west Spain. Thirty-nine populations belonging to the four main ecotypes in Galicia were cultivated in two field experiments. In the first experiment, there were significant differences between ecotypes for heading date, flag leaf length, growth habit and vigour. Coastal populations were the least infected by rust. In the second experiment, allozyme markers presumed to be diagnostic for ssp. glomerata vs. izcoi (TO1 1·00, TO1 1·03) were found in all populations studied (ten tetraploids and one diploid izcoi). This suggests that gene flow occurs between coastal and interior populations. Culm and panicle lengths of tetraploids (without application of fertilizer) were inside the izcoi range in 1997, but some plants from all populations exceeded it in 1998. It is concluded that ssp. izcoi also occurs on the coast. Populations with a high level of complementary seasonal growth were detected.

Notes: The trace gas ethylene affects plant growth and atmospheric chemistry and it interferes with soil restoration. In soil ethylene is simultaneously produced and consumed by different microorganisms. The effects of land use and soil moisture conditions on processes leading to an accumulation of ethylene are still unclear. We measured the rates at which montane and lowland soils from Austria produced and consumed ethylene over a range of water tensions and oxygen supply. Complete anaerobiosis (waterlogging, zero tension) favoured ethylene production, whereas ethylene degradation rates were greatest in soils at 30 kPa water tension. Soils from the lowland region of eastern Austria produced ethylene at rates of up to 12 pmol C2H4 g–1 h–1 under anaerobic conditions, and they consumed ethylene at rates reaching 231 pmol C2H4 g–1 h–1, after addition of 20 μl l–1 ethylene. Deciduous forest soils consumed ethylene fastest. Ethylene formed rapidly and was also consumed rapidly in soils rich in humus and total nitrogen. Soils taken from the mountains both produced and consumed ethylene more rapidly than lowland soils did. Production rates reached 146 pmol C2H4 g–1 h–1 under anaerobic conditions. Spruce forest soils produced significantly more ethylene than pastures. Ethylene formation was negatively correlated with soil pH. In montane soils ethylene production was related to the availability of simple carbon sources, expressed by the amount of extractable glucose equivalents. Maximum ethylene degradation amounted to 895 pmol g–1 h–1. Most of the soils were net sinks for ethylene at a water tension of 30 kPa and drier.

Notes: Transport of gas in soil depends on the paths available through the porous system. We determined the rate of diffusion of air through a core of silty-clay soil experimentally. After the experiment, the core was impregnated with resin and sectioned at 100 μm. From the images of the sections, the porous structure was reconstructed in three dimensions with a 100-μm resolution and gas diffusion was simulated in it. We computed the self-diffusion to describe gas diffusion. The experiment and simulation were compared using tortuosity, which is the ratio between the self-diffusion coefficient calculated in air and the coefficient calculated in the reconstructed sample. The results showed that the experimentally measured tortuosity (2.3) and the numerically calculated tortuosity (1.75) were similar. This suggests that the self-diffusion propagator is useful for simulating self-diffusion in a numerical three-dimensional reconstructed sample and that the scale chosen for the reconstruction of the sample (100 μm) was suitable for this particular soil. Résumé Les propriétés de transfert en phase gazeuse des sols ont étéétudiées sur un sol limono-argileux. A l’issue de l’expérience, l’échantillon a été imprégné avec une résine polyester puis sectionné en coupes sériées. A l’aide des images numérisées des coupes, l’échantillon a été reconstruit en trois dimensions. Nous avons ensuite simulé le processus de diffusion gazeuse dans l’échantillon reconstruit à l’aide du propagateur de diffusion. Expérience et simulation ont été comparées à l’aide de la tortuosité, rapport entre le coefficient de diffusion gazeuse dans l’air et celui dans le sol. Les résultats montrent que la tortuosité mesurée expérimentalement (2,3) et la tortuosité calculée numériquement (1,75) sont proches. Ceci prouve que l’utilisation du propagateur de diffusion est justifiée pour la simulation du processus de self-diffusion dans un échantillon tridimensionnel reconstruit et que l’échelle choisie pour la reconstruction (100 μm) s’avère pertinente pour l’étude de la diffusion gazeuse dans ce sol.

Notes: It is common practice to fit mathematical models to radionuclide activity–depth profiles in soils in order to quantify rates of vertical migration through the soil profile. We have fitted six such models to 21 different activity–depth profiles of radiocaesium (137Cs) derived from Chernobyl and determined relations between the models and the values of their parameters. The advection and dispersion parameters obtained using four solutions to the advection–dispersion equation (each based on different initial and boundary conditions or different simplifications) are in good agreement. We further develop a relation between parameter values obtained using the advection–dispersion models and those determined by a simpler exponential function of the form Aexp(–Bt) where t is the time and A and B are parameters to be estimated. One of the advection–dispersion models proved to be significantly better than the others in terms of goodness-of-fit, versatility and ease of use. A simple model, using calculations based on measured characteristics of the activity–depth profile, was shown to accord well with parameters derived from more complex models based on statistical curve fitting. We have also evaluated the ‘residence time’ or ‘compartmental’ model approach to characterizing radionuclide activity–depth profiles. We relate such models to a numerical solution of a simple advection equation, and we show that apparent dispersion in compartmental models is an artefact of numerical dispersion, which can be quantified by the Courant condition. For activity profiles that have a significant advection component, using solutions to the advection–dispersion equation, we have observed a strong positive correlation between advection and dispersion in the profile.

Notes: Immobilization of nitrogen (N) in decomposing straw varies between soils, and the objective of this study was to identify the mechanisms responsible. Internode segments of wheat straw were incubated in Denmark and in Scotland in arable soils fertilized with NH4NO3, labelled with 15N, for periods up to 1 year. Straw was recovered from the soils periodically and analysed for microbial biomass and different forms of N using chemical methods and CPMAS 15N NMR spectroscopy.The total N content of the straw increased, as long as the soil was not too wet, such that there was overall immobilization. This was accompanied by a rapid increase in the content of amino acid N and to a lesser extent of glucosamine N and a concomitant decrease in the carbohydrate content of the straw. Using direct and plate counts for bacterial and ergosterol content for fungal estimation, we found that fungal biomass was much greater than that of bacteria. This correlated with the forms of N in the straw as determined by CPMAS 15N NMR, which showed spectra that were more typical of fungi than of bacteria. It seems that immobilization of N is primarily caused by fungi as they decompose the straw.

Notes: Collecting soil data is time-consuming and costly, often exceeding practical possibilities. A methodology for the delineation of soil mapping units in an alluvial plain of Western Peloponnese, Greece, was investigated. A detailed soil survey of an area of 300 ha was used to obtain the basic soil data for evaluating the performance of the proposed methodology. The methodology consists of the following steps: (a) data collection from borings and representative soil profiles, (b) definition of the soil mapping units in the study area, (c) determination of the range of the diagnostic variables for each mapping unit from field observations and statistical analysis of the analytical data from representative soil profiles, (d) determination of the class of each diagnostic variable by observation at a network of boring points, (e) subjective assignment of numerical values to soil variables at the bore points, (f) estimation of the values of each soil variable at the points of a regular grid using the interpolation methods kriging and inverse squared distance, (g) application of the fuzzy set theory to the interpolated data and the production of thematic fuzzy maps, and (h) validation of the results through a number of independent test borings. The results obtained show that the proposed methodology can produce soil maps of recent alluvial plains with acceptable accuracy and cost.

Notes: Fundamental knowledge about the complex processes during the decomposition, mineralization and transfer of residue organic matter in soils is essential to assess risks of changes in agricultural practices. In a double tracer (13C, 15N) experiment the effect of maize straw on the mineralization dynamics and on the distribution of maize-derived organic matter within particle size fractions was investigated. Maize straw (a C4 plant) labelled with 15N was added to soils (13.2 g dry matter kg–1 soil) which previously had grown only C3 plants, establishing two treatments: (i) soil mixed with maize straw (mixed), and (ii) soil with maize straw applied on the surface (surface). Samples were incubated in the laboratory at 14°C for 365 days. The size fractions (〉 200 μm, 200–63 μm, 63–2 μm, 2–0.1 μm and 〈 0.1 μm), obtained after low-energy sonication (0.2 kJ g–1), were separated by a combination of wet-sieving and centrifuging. The mineralization of maize C was similar in the two treatments after one year. However, decomposition of maize particulate organic matter (predominantly in the fraction 〉 200 μm) was significantly greater in the mixed treatment, and more C derived from the maize was associated with silt- and clay-sized particles. A two-component model fitted to the data yielded a rapidly mineralizable C pool (about 20% of total C) and a slowly mineralizable pool (about 80%). Generally, the size of the rapidly mineralizable C pool was rather small because inorganic N was rapidly immobilized after the addition of maize. However, the different mean half-lives of the C pools (rapidly decomposable mixed 0.035 years, and surface-applied 0.085 years; slowly decomposable mixed 0.96 years, and surface-applied 1.7 years) showed that mineralization was delayed when the straw was left on the surface. This seems to be because there is little contact between the soil microflora and plant residues. Evidently, the organic matter is more decomposed and protected within soil inorganic compounds when mixed into the soil than when applied on the soil surface, despite similar rates of mineralization.

Notes: Mercury porosimetry data can be interpreted in terms of soil structure using ideas drawn from (i) network modelling and percolation theory and (ii) fractal geometry. We linked mercury intrusion to soil structure quantified by image analysis within a relevant common pore radius scale. We compared (i) three independent methods for computing fractal dimensions of the matrix and of the solid–pore interface, namely fitted square boxes method and pore chord distribution on scanning electron microscope images of soil thin sections, and mercury porosimetry, and (ii) two independent methods for characterizing pore connectivity (image analysis) and percolation process (pressure threshold from mercury porosimetry). The results from analyses of the pore size distribution by mercury porosimetry differed from those from the image analysis. Mercury intrusion is controlled by both the connectivity of the pore space network and locally by pore throats leading to larger pore bodies. By contrast, image analysis is unaffected by pore connectivity and measures pore bodies. On the other hand, the chord length method might not adequately capture the scaling properties of the solid–pore interface, whereas the mercury porosimetry data were also difficult to interpret in terms of fractal geometry because of the effects of pore connectivity. However, fractal dimension values of both the solid phase and the solid–pore interface increased as a function of clay content, whereas both percolation probability values and throat radius values at the mercury percolation threshold decreased. The results show the merit of applying both fractals and percolation theory for determining structural parameters relevant to mercury and water transport in soil.

Notes: Soil microbes produce exudates which upon drying become water-repellent, thus altering hydraulic properties. The influence of microbial activity caused by adding plant nutrients on the hydraulic characteristics of soil aggregates is reported. Soil aggregates were collected from a field that had been fertilized with different amounts of nitrogen. Aggregates were also incubated with different nutrient treatments in the laboratory. Their sorptivity, hydraulic conductivity and water repellency were measured with a new device. Adding nitrogen was found to decrease sorptivity and hydraulic conductivity because of increased water repellency in the field. In the laboratory studies, the addition of nutrients caused severe water repellency in the soil aggregates. Respiration studies identified a large increase in biological activity following nutrient amendment which produces water-repellent materials.

Notes: Solution cadmium (Cd) concentrations and sorption and desorption of native and added Cd were studied in a range of New Zealand soils. The concentration of Cd in solution and the concentrations and patterns of native soil Cd desorbed and added Cd sorbed and desorbed varied greatly between the 29 soils studied. Correlation analysis revealed that pH was the most dominant soil variable affecting solution Cd concentration and sorption and desorption of native and added Cd in these soils. However, organic matter, cation exchange capacity (CEC) and total soil Cd were also found to be important. Multiple regression analysis showed that the log concentration of Cd in solution was strongly related to soil pH, organic matter and total Cd, which in combination explained 76% of the variation between soils. When data from the present study were combined into a single multiple regression with soil data from a previously published study, the equation generated could explain 81% of the variation in log Cd solution concentration. This reinforces the importance of pH, organic matter and total Cd in controlling solution Cd concentrations. Simple linear regression analysis could at best explain 53% of the total variation in Cd sorption or desorption for the soils studied. Multiple regression analysis showed that native Cd desorption was related to pH, organic matter and total Cd, which in combination explained 85% of the variation between soils. For sorption of Cd (from 2 μg Cd g–1 soil added), pH and organic matter in combination explained 75% of the variation between soils. However, for added Cd desorption (%), pH and CEC explained 77%. It is clear that the combined effects of a range of soil properties control the concentration of Cd in solution, and of sorption and desorption of Cd in soils. The fraction of potentially desorbable added Cd in soils could also be predicted from a soil’s Kd value. This could have value for assessing both the mobility of Cd in soil and its likely availability to plants.

Notes: The magnitude of radiocaesium fixation by micaceous clay minerals is affected by their transformation, which depends on weathering in soil. The net retention of radiocaesium traces was quantified by sorption–desorption experiments in the various horizons of four sandy soils forming an acid brown earth–podzol weathering sequence derived from sandy sediments and characterized by marked changes in mineral composition. The features of the 2:1 minerals of the four soils, resulting from an aluminization process in depth and a desaluminization process towards the surface, had a strong influence on Cs+ fixation. Beneath the desaluminization front, which deepens from the acid brown earth to the podzol, hydroxy interlayered vermiculite was dominant and the 137Cs+ fixation was the weakest. At the desaluminization front depth, vermiculite was responsible for the strongest 137Cs+ fixation. In the upper layers, smectite appeared in the podzolized soils and the 137Cs+ fixation decreased. The magnitude in Cs+ fixation therefore appeared as a tracer of the transformation process affecting the 2:1 clay minerals in the acid brown earth–podzol weathering sequence. This magnitude was positively correlated with the vermiculite content of the studied soil materials estimated by the rubidium saturation method.

Notes: Rhizosphere processes strongly influence the availability of phosphorus (P) to plants. Organic ligands that are exuded from the root surface mobilize phosphorus by dissolution of P minerals or by desorption of adsorbed phosphate. We developed a mechanistic model to study the mobilization of phosphate sorbed on goethite by the exudation of citrate and consequent uptake of phosphate by the root. The use of a model allows the effects of the organic anion and pH on P desorption to be separated. The model is also used to predict concentration profiles developing around the root for phosphate, citrate (with or without accounting for degradation) and pH, providing insight into the processes that occur in the rhizosphere. Results of model calculations show that with larger rates of citrate exudation, greater P availability is predicted. Exudation at a rate of 0.5 μmol citrate m–1 root day–1, which is in the range found for P-deficient plants, increased P availability almost 2-fold at fairly large phosphate loading of goethite (1.9 μmol m–2) and almost 30-fold at small phosphate loading (1.3 μmol m–2). Competitive adsorption causes a much greater relative increase in the phosphate concentration in solution at small than at large phosphate loading, which explains this result. Simultaneous acidification of the rhizosphere results in a smaller P mobilization than at a fixed pH of 5, as a result of the pH dependence of phosphate adsorption in the presence of citrate. Sorption of citrate increases its persistence against microbial decay, and hence has a positive effect on the mobilization of adsorbed phosphate.

Notes: Standard procedure for dispersing natural soils for particle size fractionation can be too aggressive for soil containing coal and other organic particles from coal industries. We have investigated ultrasonic dispersion for the latter in four soils differing in pedogenesis (Phaeozem, Podzol, reclaimed mine soils), carbon content (27.5–138.6 g kg–1), clay content (80–153 g kg–1) and sources of particles (airborne coal dust, combustion residues, lignite particles). As we found previously for natural soils, the ultrasonic energy needed for complete dispersion varies between 450 and 500 J ml–1, but the resulting particle size distributions differ from those obtained by standard textural analysis. This is probably related to the different properties of native soil organic matter and coal and combusted particles. Coal and soot particles may partly resist oxidation with hydrogen peroxide, depending on material and particle size. The diameter of lignite particles, remaining after oxidation, is overestimated in sedimentation analysis by a factor of 1.66. Sand-sized lignite particles can be disrupted by ultrasonication and redistributed to finer particle size fractions. The ultrasonic dispersion and particle size fractionation procedure can be applied to soils containing coal and combusted particles, but caution is needed in interpreting the results if they contain large proportions of coal particles.

Notes: The resistance of soil humic material in soil seems to depend on its molecular characteristics. We have studied 12 chemically modified humic samples in mineral soil to identify molecular characteristics connected with their resistance to biodegradation. The treatments, used to introduce changes in the composition, reactivity and N content of the original humic acid (HA), consisted of acetylation, amidation, methylation, nitration, oximation, sulfonation, acid and alkaline hydrolysis, hydrogen peroxide oxidation and fixation of ammonia and of stearic acid. The relations between respirometric data (total mineralization coefficient and mineralization rates at different stages of the 85-day mineralization curve) and a series of HA characteristics (wet chemical analyses, spectroscopic data and relative yields of different pyrolysis products) were examined by correlation and by multiple regression models.The structural characters of HA most directly related to the susceptibility to biodegradation were the concentration of O–alkyl structures and oxygen content mainly in phenolic structures. The aliphatic/aromatic ratio showed no significant correlation with the resistance of organic matter. The fact that the yields of lignin-derived pyrolysis compounds were directly related to carbon mineralization suggested that the structural contribution of domains consisting of preserved biomacromolecules is more important than the neat aromaticity of the HA.The multiple regression models suggested that molecular size of the HA had the greatest influence in the early mineralization stages, whereas the characteristics of the C backbone (concentration of O–alkyl and alkyl carbons) had a greater influence in the advanced transformation stages in the soil.The lack of significant correlations with N concentration indicated that native N forms as well as N-containing groups introduced in peat HA did not have a measurable trophic effect on the biodegradation. The disordered macromolecular structure of the HAs seems to play a greater role in their resistance to biodegradation than the relative amounts of their major structural constituents.

Notes: A simple equation to describe sorption of anions by goethite would be useful as a means of characterizing batches of goethite and in studies of plant uptake of anions from the sorbed form. A suitable relationship between solution concentration (c) of phosphate or citrate and their sorption (S) by goethite at a constant pH or at different pH values is 〈displayedItem type="mathematics" xml:id="di-fml-1" numbered="no"〉〈mediaResource alt="image" href="urn:x-wiley:13510754:EJSS203:EJSS_203_m1"/〉 where b is a parameter, SMax is the maximum sorption, and a is a parameter at constant pH. In the middle range of sorption (from about 30% to about 70% of maximum sorption) this equation approximated to a Tempkin equation, but the full equation is more useful as it applied over the whole sorption range. The values of a varied with pH. This variation could be explained by changes in the electric potential of the adsorbing surfaces and in the degree of dissociation of the anions. The parameter a could therefore be replaced by a function of pH. The effects were consistent with formation of bidentate phosphate complexes and tridentate citrate complexes with the goethite surface.

Notes: The composition of soil organic matter (SOM) is influenced by land use and fertilization. We studied changes in the SOM in a long-term field experiment on a sandy Podzoluvisol. The control plot and four combinations of manurial treatments of the experiment were selected: one with mineral fertilizer only and three combinations of organic manure with mineral fertilizer: cattle manure + NPK, cattle manure + PK and straw + NPK. The SOM was extracted by sodium pyrophosphate solution (pH = 10) and hot water (100°C). The extracts were analysed by Fourier-Transform Infrared (FT-IR) spectroscopy and gel permeation chromatography (GPC). The FT-IR spectra from sodium pyrophosphate extracts indicate that composition of SOM is indeed influenced by different fertilization. The C=O band at 1710 cm–1 in the samples of the plots fertilized with cattle manure has the highest absorption intensity, whereas the material from the plot fertilized with straw + NPK has the least intense. The GPC analyses of the extracts showed that adding cattle manure + NPK increased the molecular size of SOM in comparison with the control plot. The analysis of hot-water extracts with FT-IR showed no significant differences in functional groups, but GPC chromatograms distinguished features in molecular size distribution. Fertilization with cattle manure increased the molecular size of the SOM in comparison with the control, but the differences in content of carboxylic groups and molecular weight were detected in sodium pyrophosphate extracts only.

Notes: The dynamics of water flow in soils influences the transport behaviour of solutes. Transport of bromide and herbicides through undisturbed soil columns was investigated under conditions of unsaturated steady-state and transient water flow. Effective transport parameters were obtained from fitting the convection–dispersion equation to curves of concentration against cumulative drainage, and these enabled us to interpret the observed behaviour. Under both steady-state and transient flow bromide and herbicides were transported through similar parts of the pore volume of a homogeneous single grain soil (Bv horizon). However, in aggregated Ah and Ap horizons preferential transport occurred during transient flow but not during steady-state flow. For preferential flow the mean transport volume seemed to depend on the prevailing pore system and the fraction of preferentially flowing water. Solute leaching was more efficient under steady-state than under transient flow for bromide in all soils and for herbicides in the Bv horizon. However, when transient flow caused preferential transport, herbicide loss was greater under transient flow than during steady-state flow. Under preferential flow conditions a three-step herbicide concentration development recurred in successive drainage events. This behaviour was not observed for the non-reactive tracer. It seemed to be caused by sorption. A steady-state model with cumulative drainage as independent variable instead of time can predict the transport of non-reactive and adsorbed solutes in homogeneous soils without features of preferential flow. Otherwise constant effective input parameters cannot be assessed a priori.

Notes: The surface morphology of quartz grains can indicate the degree of weathering of soil material. We have compared two methods of assessing the relative weathering of soils on the basis of differences in the surface morphology of quartz sand grains in a catena of soils in Rwanda. One method is based on the presence or absence of surface features indicative of weathering or freshness, while the other uses the frequency and size of dissolution etch pits. A ranking of relative weathering could be obtained using the first technique for the slightly and somewhat weathered soils but not for the weathered soils. On the other hand, weathering trends and differences between the horizons studied were detected in weathered soils using the second method. The introduction of more specific definitions of the weathering classes used in the latter method leads to clear improvement of the inter-observer reproducibility of the weathering classification.The surface features on the quartz grains suggest that the soil at the summit is less weathered than the other soils of the sequence. Quartz grains from the well-drained soils on the slopes, which are subjected to more intense leaching and thus to stronger chemical weathering, have more triangular etch pits and chatter marks. In the imperfectly drained soils in the valley bottom quartz grains are less etched because dissolution is inhibited by the oversaturation in silica of the drainage waters.

Notes: An experiment was conducted to examine the response of herbage grazed by dairy cows to sodium fertilizer applied with or without sulphur fertilizer. The residual effects of applying sodium fertilizer in the previous year were also evaluated. The application of sodium or sulphur fertilizer did not affect herbage growth or height, but the application of sulphur fertilizer increased crude protein content of herbage. The increase in sodium content with application of sodium fertilizer was small, but was greater when sodium fertilizer had also been applied in the previous year. Herbage potassium was increased when sodium fertilizer was applied in the year of the experiment, but only if sodium had not been applied in the previous year. Application of sulphur fertilizer increased herbage sulphur content and reduced the contents of boron, chromium, molybdenum and nickel. Cows grazing pasture that had received sodium fertilizer had increased milk yields and the content of lactose in milk, whereas those grazing pasture that had received sulphur fertilizer application had reduced milk yields and the content of milk fat.

Notes: The content of ergot alkaloids (ergovaline and chanoclavine), and their production in October 1996 and during the period May–September 1997, were investigated in seventeen ecotypes of perennial ryegrass (Lolium perenne L.) and in nineteen ecotypes of meadow fescue (Festuca pratensis Huds.), naturally infected with Neotyphodium spp. The ecotypes were collected in the north-eastern part of the Czech Republic. In 1996 the content of ergovaline in the ecotypes of perennial ryegrass ranged from 0·00 to 2·73 μg g–1 dry matter (DM) (one cut), and in 1997 from 0·00 to 4·65 μg g–1 DM (five cuts). In meadow fescue the content of ergovaline varied from 0·00 to 0·61 μg g–1 DM (one cut) in 1996, and in 1997 from 0·00 to 2·31 μg g–1 DM (five cuts). The content of chanoclavine (investigated in 1997 in four cuts only) in perennial ryegrass ranged between 0·00 and 3·39 μg g–1 DM, and in meadow fescue between 0·00 and 2·26 μg g–1 DM. Most ecotypes of L. perenne reacted to the high temperature and heavy rainfall in June and July of 1997 with an enhanced production of ergovaline, whereas the content of chanoclavine was not changed. Such reaction to stress conditions was not observed in the ecotypes of F. pratensis. Large differences in the production of both ergot alkaloids between different ecotypes of both plant species were observed.

Notes: Seeds of sainfoin (Onobrychis viciifolia Scop.), alfalfa (Medicago sativa L.) and Astragalus adsurgens Pall. were put into and recovered from earth orbit from China in 1994. The isoenzymes in leaves and flowers and amino acids in leaves of the first post-flight generation of plants were analysed. Germination of their seeds under NaCl and polyethylene glycol (PEG) stress were tested. The results showed that space flight caused some changes in the pattern of peroxidase and esterase in sainfoin and A. adsurgens, and variation in amylase patterns in leaves of alfalfa. The total amino acid content increased in sainfoin and alfalfa leaves. Under salt and water stress, the proportion of progeny seeds of alfalfa and A. adsurgens, respectively, germinating after 7 and 18 days’ imbibition was not markedly different between treatment and control. However, the progeny seeds of sainfoin which had been in orbit exhibited higher tolerance to NaCl and PEG stress during germination.

Notes: A simple model of depletion by a grazing ruminant was tested at high depletion levels and for different sizes of feeding station. The model divides an initially uniform sward profile into grazing horizons, dependent on bite depth, and assumes a constant within-horizon bite area.Two grazing experiments were conducted using cattle. Uniform areas of oats and alfalfa herbage were grazed individually to a high level of depletion (≈150 bites m–2). Feeding station areas were 0·16, 0·30 and 0·53 m2. Initial sward heights were 10 and 20 cm in oats and 20 cm in alfalfa. Size of area did not significantly affect the observed number of bites removed per square metre, the mean residual herbage height or mass, or the proportion of each grazing horizon depleted, derived from the frequency distribution of residual heights. The mean residual height for all treatment combinations could be explained by assuming a ‘take half’ rule for mean bite depth, and allowing for the derived proportion of the area of each grazing horizon depleted.A simple bite placement simulator was used to generate, for a single grazing horizon and for a given maximum potential area of a bite, the expected relationship among mean effective area of a bite, the proportion of the area of the horizon grazed and the number of bites removed per unit area. The simulator mimics a loosely systematic grazing style. The observed bite numbers and the derived proportions of grazing horizon depleted can be reconciled if the within-horizon mean effective bite area is not constant but declines as predicted by the bite placement simulator. The implications for the shape of the gain function within a feeding station are discussed.

Notes: Variability in the response to salinity within Chloris gayana (Rhodes grass) germplasm was evaluated under field conditions, and random amplified polymorphic DNA (RAPD) markers were used to assess genetic relatedness among cultivars/accessions. RAPD analysis showed a clustering of cultivars of known relatedness: cv. Pioneer and accessions Local and Trancas (derived from an old Pioneer pasture established in saline soil) belonged to the same cluster, Katambora to another and tetraploid Boma and Callide could be further separated, Boma belonging to a fourth, distant cluster. Field experiments were laid out in two types of plots: control [with electrical conductivity of the saturation extract (EC) = 3·64 dS m−1] and saline (EC=13·10 dS m−1) and two experiments were carried out: one to evaluate the effects of salinity on emergence and establishment, and the other, with a uniform number of plants per plot, to evaluate yield under saline conditions. All cultivars/accessions had salinity-associated decreases in dry-matter (DM) production during the establishment phase. After this stage, 1-year DM yield was similar in all cultivars within each salinity level and production in the saline plots was significantly lower than in controls only in cv. Callide and accession Trancas. Second-year production in the non-saline plots increased by 30% on average over the previous year, whereas an average 40% reduction was observed in the saline plots. Thus, salinity had a negative effect on Rhodes grass establishment and persistence. The cultivars could not be ranked unequivocally by production under saline conditions, but tetraploids Boma and Callide may be said to be less tolerant than the rest on the basis of an increased proportion of dead leaves and decreased number of stolons observed in the saline plot.

Notes: Three experiments were carried out to study the effects of feeding lucerne silage (wilted to give different dry-matter (DM) contents) and ventilated hay to dairy cows on milk production, milk quality, milk-renneting properties, clostridial spore content and the quality of cheese prepared from the milk. The lucerne, cut at vegetative or early-bud stages of maturity, was harvested from alternate windrows and conserved as silage or artificially dried hay. The lucerne was wilted until it reached different DM contents of 550, 360 and 432 g kg–1 in the three experiments, harvested, chopped with a self-loading forage wagon and ensiled in low and narrow clamps made up of transferable prefabricated walls. The organic acid content, pH, yeast and mould counts of the lucerne silage suggested that there was no aerobic deterioration.In each experiment, fifty Italian Friesian lactating cows were divided into two groups and fed two maize silage-based rations for 6 weeks, which only differed in the lucerne forage [silage (S) vs. ventilated hay (H)], in a cross-over experimental design. The lucerne in the rations represented 35%, 23% and 24% of the DM of the rations for the three experiments. The microbiological profiles of the ration were influenced more by the maize silage than by the lucerne silage.Individual daily DM intakes were similar for the two treatments in Experiments 1 and 3 (on average 18·7 kg in Experiment 1 and 20·3 kg in Experiment 3) and slightly lower for S cows in comparison to H cows in Experiment 2 (18·0 vs. 19·0 kg). Milk yields of S and H cows were 21·0 and 20·8, 20·0 and 20·6 (P 〈 0·01), and 28·4 and 27·9 kg d–1 in Experiments 1, 2 and 3 respectively. Milk composition was similar for all the experiments for the two treatments, except that the protein content was lower and the fat content was higher in the silage treatment than in the hay. The renneting properties and microbiology of the milk were not influenced by the introduction of lucerne silage into the rations, although the season in which it was consumed had a greater effect on the microbiological content, in terms of standard bacterial counts, proteolytic, coli and lactic acid bacteria, and clostridia spores. The clostridial spore counts were always very low (〈 400 per litre), thus fulfilling the requirements for top-quality milk for Grana cheese production. In the third experiment, the quality of Grana Padano cheese produced was examined, and no differences between treatments were observed after 12 months of maturation.These results show that lucerne silage can be included in the ration of dairy cows instead of ventilated lucerne hay, which is considered to be the top-quality hay available for the production of milk destined for Grana cheese, without any negative effects on milk and cheese quality.

Notes: This study assesses near-infrared reflectance spectroscopy (NIRS) as a rapid and reliable method for estimating the clover content of clover/grass mixtures. NIRS calibrations were derived both for red clover (Trifolium pratense) and white clover (T. repens) mixtures. To maximize variability, the set of samples for analysis was selected from field plots receiving a wide range of nitrogen fertilizer application and harvested over 2 years at weekly intervals throughout the growing season. The samples were scanned using a NIR-Systems model 5000 monochromator. A total of 183 white clover/grass and 282 red clover/grass samples were used to calibrate and cross-validate the equations derived. The NIRS calibrations obtained from this study have the potential to be used in research on clover/grass mixtures as well as for advisory work.

Notes: On two occasions in 1995 the effects of severity of conditioning and pressing on the drying rate of Italian (Lolium multiflorum Lam.; IR) and perennial (Lolium perenne L.; PR) ryegrass were assessed. The cut herbage was either left unconditioned (0C) or conditioned by passing through a laboratory scale macerator three (3C) or six (6C) times after which the herbage was placed into wire mesh trays and either left unpressed (0P) or pressed (P). The twelve treatment combinations (two species × three conditioning levels × two pressings) were replicated three times on each of the two occasions. On each occasion the trays plus herbage were weighed at hourly intervals over an ≈6-h period and dry matter (DM) of the herbage was estimated from the weight change. Pressing treatment gave, on average, a negative effect on the drying rate of the herbage. Conditioning significantly (P 〈 0·001) increased the drying rate of the herbage. Compared with 0C herbage, the greatest increase in drying rate was obtained with the 3C treatment; the additional conditioning imposed by the 6C treatment produced only a small further increase in drying rate. Drying constants calculated for each treatment combination showed that IR dried 1·45 times faster than PR and that, compared with 0C, herbage drying rate was 2·18 and 2·38 times greater as a consequence of the 3C and 6C treatments respectively. Pressed (P) herbage dried at a rate of 0·90 of that of the unpressed herbage (0P).

Notes: Five methods for increasing the botanical diversity of permanent grassland, either by sowing site-specific species-rich grass/forb seed mixtures (strip-seeding; or over-sowing after sward disturbance by light harrowing, partial rotary cultivation or turf removal), or by introducing transplanted plug plants, were compared with a control treatment in replicated field experiments on six farm sites in Environmentally Sensitive Areas (ESAs) in England and Wales. Effects on herbage production under hay cutting in July and on botanical composition were recorded in the two subsequent years. Turf removal before sowing was the only treatment that significantly reduced herbage production; this treatment also had the greatest effect on increasing botanical diversity (to a mean of twenty-eight plant species per site compared with fifteen species for the control 2 years after sowing). The least successful establishment of sown species resulted from light harrowing before sowing; the rotary-cultivated and strip-seeded treatments increased species diversity, although by less than turf removal. Successful establishment of introduced species was greatest on sites having a low soil nutrient status. Species that established successfully from seed on most sites and treatments included the grasses Alopecurus pratensis, Cynosurus cristatus, Festuca rubra and Phleum pratense, and the forbs Achillea millefolium, Leucanthemum vulgare, Plantago lanceolata and Prunella vulgaris; in addition, Centaurea nigra, Hypochoeris radicata and Lotus corniculatus were also established by one or more methods on most sites. Lychnis flos-cuculi established successfully on mesotrophic sites, and Medicago lupulina on calcareous sites. Several species failed to establish at all or most sites where they were sown, e.g. Helianthemum nummularium, Pimpinella saxifraga and Rhinanthus minor. Most transplanted plug-plant species established successfully in the short term, but many failed to persist or their frequency in the sward remained low; exceptions included A. millefolium and P. lanceolata. The results are discussed in relation to the requirements for management to further the objectives of ESAs and agri-environmental schemes.

Notes: The expected reduction in the use of fertilizer nitrogen (N) on grassland in the Netherlands has led to renewed interest in white clover. Therefore, the performance of a newly sown perennial ryegrass/white clover sward on clay soil was assessed during 4 consecutive years. The experiment consisted of all combinations of two defoliation systems, i.e. one or two silage cuts per year (S1, S2), spring N application rate, i.e. 0 or 50 kg ha−1 year−1 (N0, N50), and the management system, i.e. rotational grazing and cutting, or cutting only (RGC, CO). The overall mean white clover cover was 30%. All treatments affected white clover cover, which was 8% higher with S2 than with S1, 6% higher with N0 than with N50 and 12% higher with CO than with RGC. The overall mean annual dry-matter (DM) yield (13·1 t ha−1 year−1) was significantly affected only by the management system: in two relatively wetter years, the annual DM yield was 1·19 t ha−1 higher with RGC than with CO, whereas there was no difference in two relatively drier years. Nitrogen application increased the DM yield in the first cut by 7·0 kg kg−1 N applied, but had no significant effect on the annual DM yield. Herbage quality was not affected by the experimental treatments. The average in vitro organic matter digestibility was 0.801, and the average crude protein content was 193 g kg−1 DM. With the expected reduction in the use of fertilizer N, perennial ryegrass/white clover swards should be seriously considered as an alternative option to perennial ryegrass swards on these clay soils.

Notes: The effects of NH4+−N, NO3−-N or urea-N addition on N uptake by perennial ryegrass (Lolium perenne), 15N recovery and pH of four limed and unlimed soils were determined in a pot experiment over 10 weeks. The best form of N in terms of herbage N uptake and fertilizer recovery differed between the soils. Recovery of applied 15N in herbage was greatest for the soil with the lowest pH and highest organic matter content, and overall recovery in the soil–plant system was influenced more by soil type than by N form. There was always an apparent mineralization of soil N when perennial ryegrass was present. Soil pH changes ranged from +0·11 to −0·58 units for the unlimed soils and from +0·03 to −1·06 units for the limed soils. The use of NO3−-N rather than NH4+-N or urea-N avoided further acidification. Lime increased herbage N only from the soils with the lowest pH values. Although the acidifying effect of NH4+-fertilizer was alleviated, liming may increase nitrification and possibly N loss via denitrification and/or leaching in the field.The major influence of soil type on herbage N uptake, pH response to the treatments and fertilizer recoveries implies that liming and fertilizer N management decisions should consider soil characteristics, such as organic matter, clay contents and pH.

Notes: The extent of semi-natural grassland has diminished considerably across lowland landscapes of England and Wales during the second half of the twentieth century. Locating, describing and evaluating the dwindling cover has been a major challenge for conservationists. A concentrated vegetation survey effort at grassland sites has been mounted within different parts of Britain since the late 1970s. Plant community recognition has benefited considerably from the development of the contemporary National Vegetation Classification, and its widespread adoption permits national inventory of comparable vegetation data.Findings of a range of surveys (ninety-eight in total), undertaken between 1978 and 1996 in England and Wales covering different forms of unimproved lowland grassland, are collated and reviewed. Vegetation data were abstracted from internally published survey reports. Calcicolous and neutral grasslands have been covered more thoroughly than acidic and wet or marshy grasslands. Cover data are summarized at community level. Overall estimates from survey results indicate that there are some 27 500–40 000 ha of calcicolous grassland, 7500–15 000 ha of unimproved neutral pasture and hay meadow, 8000–15 000 ha of acidic grassland and 9000–17 500 ha of wet grassland in lowland England and Wales; these represent only 1–2% of the cover of permanent lowland grassland. Some communities have additional representation in heathlands, mires and upland environments.Although they require further refinement, the cover data for individual communities provide a context for assessing priorities in site-based and agri-environment conservation programmes. It is concluded that, as well as arresting further depletion, it will be necessary to restore and expand lowland grassland habitats to counteract the negative impacts of fragmentation and isolation of various community types, such as the Centaureo–Cynosuretum, which is widely but thinly distributed. Habitat rehabilitation schemes also need to assimilate local patterns of community diversity characteristic of both wet and dry grasslands. It is suggested that reversal of the recent successional trends that followed relaxation of grazing at certain sites might produce a more appropriate balance in the relative cover of coarse tall grasslands and fine short turf. Vegetation surveys provide a source of spatial data for identifying local aggregations of semi-natural grassland remnants.

Notes: Burns’s equation for describing solute movement through soil is attractive because it is simple and predicts adequately in many instances. However, the assumptions implicit in it are not inconsistent with preferential solute flow. We have explored the consequences of this by leaching initially resident chloride and surface-applied tritium and nitrate through 250-mm-long intact cores of a silt loam soil. The applied flow rates of 3 and 5 mm h−1 (realistic rainfall intensities) produced unsaturated soil conditions, except near the base where free water dripped out. Burns’s equation described the movement of the three solutes fairly successfully with the water content parameter having values between 0.29 and 0.48, similar to the actual volumetric water content of 0.47.The leaching of resident chloride to 450-mm-deep mole drains in the field was also successfully simulated using Burns’s equation. However, simulation of the leaching of bromide applied to the soil surface as a solid salt was problematic. This resulted from uncertainty as to whether to treat the application as a pulse input to the flux or resident concentration. The observed behaviour fell about midway between the simulations for these contrasting initial conditions.

Notes: Pollution of ground water caused by excessive and uncontrolled use of nitrogen fertilizer is worrying. A recent example of such pollution has been observed in an agricultural basin in the province of Nevsehir, Turkey, where up to 900 kg ha−1 nitrogen fertilizer is used for growing potatoes in sandy soils under irrigation. Using nitrogen fertilizer in amounts that guarantee large yields without polluting ground water is essential. We present results of field experiments and numerical simulations involving 15N-labelled nitrogen fertilizer leaching. In the field, we monitored the movement of water and the distributions of nitrogen species within the soil–water–plant continuum. The detailed dynamics of the nitrogen cycle within the system were simulated. Simulations included calibration and validation of the nitrogen version of the LEACHM model (LEACHN, version 3) and long-term applications of the model. The model’s predictions of nitrogen fluxes under long-term use of fertilizer and irrigation were analysed. Nearly half of the applied ammonium-N was converted to nitrate-N during the growing season. With increasing additions of N the rate of plant uptake declined, while leaching increased significantly, and the fraction of nitrogen remaining in the soil profile increased only moderately. In long-term applications, a significant fraction of the applied fertilizer tended to accumulate after the first year in soil as the residual nitrogen not taken up by the crop. Accumulated residual nitrogen is converted to nitrate-N and leached rapidly from the soil profile during the wet season following the harvest. To reduce leaching of the residual nitrate, the rates, frequencies and timings of fertilizer application and irrigation must be scheduled in accordance with the plant growth periods and the hydraulic regime of the soil.

Notes: Two sequential extractions with unbuffered 0.1 m BaCl2 were done to study the release of salt-exchangeable H+ and Al from mineral horizons of five Podzols and a Cambisol. Released Al was found to have a charge close to 3+ in all horizons and in both extractions. This finding was supported by the near-equality of the titrated exchangeable acidity (EAT) and the sum of exchangeable acids (EA = He + 3Ale, calculated from the pH and Al concentration of the extract). The ratio between EA of the second and the first extraction was over 0.50 in the Bs2 and C horizons and smaller in the other horizons. H+ was assumed to be in equilibrium with weak acid groups, and the modified Henderson–Hasselbach equation, pKHH = pH − n log (α/(1 − α)), was used to explain pH of the extract. The degree of dissociation (α) was calculated as the ratio between effective and potential cation exchange capacity. Value of the empirical constant n was found to be near unity in most horizons. When the monoprotic acid dissociation was assumed in all horizons, pKHH had the same value in both extractions. For Al3+, two equilibrium models were evaluated, describing (i) complexation reactions of Al3+ with soil organic matter, and (ii) equilibrium with Al(OH)3. Apparent equilibrium constants were written as (i) pKo = xpH − pAl3+, and (ii) log Qgibbs= log Al3+ − 3log H+. The two extractions gave an average reaction stoichiometry x close to 2 in all horizons. Results suggest that an equilibrium with organic Al complexes can be used to express dissolved Al3+, aluminium being apparently bound to bidentate sites. The value of log Qgibbs was below the solubility of gibbsite (log Kgibbs = 8.04) in many horizons. In addition, log Qgibbs of the second extraction was greater than that of the first extraction in all horizons except the C horizon. This indicates that equilibrium with Al(OH)3 cannot explain dissolved Al3+ in the soils. We propose that the models of pKHH and pKo can be used to simulate exchangeable H+ and Al3+ in soil acidification models.

Notes: The conformational structure of dissolved humic substances is an important property that controls the reactivity of humus in the soil solution. High performance size-exclusion chromatography was used here to study the changes in molecular size of different humic substances brought about by addition of mineral (HCl) and monocarboxylic (formic, acetic, propionic, and butyric) acids. The CPMAS-NMR spectra showed that humic substances had varying chemical composition and that the ratio of hydrophilic to hydrophobic carbon (HI/HB) was greater for a humic acid from soil than for ones from oxidized coal and lignite. All humic substances showed a decrease in UV absorbance of chromatographic peaks when treated with either HCl or monocarboxylic acids. This was due to the hypochromic effect by which the absorptivity of associated molecules is decreased when they are separated. We attributed the molecular separation upon acid treatment to the formation of intermolecular hydrogen bonding that alters the original conformation stabilized mainly by weaker hydrophobic interactions. Addition of organic acids not only further decreased peak absorbances of humic acids but also caused their shift to larger elution volumes, indicating a larger conformational disruption than with HCl. The extent of the molecular size changes showed a relation to the number of carbons of monocarboxylic acids and to the HI/HB ratios of humic materials. The larger the carbon content of organic acids and the smaller the HI/HB ratio of humic materials, the larger was the decrease of the average molecular size of humic acids. These results suggest that dissolved humic substances associate predominantly by hydrophobic forces and that the apolar components of humic substances largely control their aggregation and reactivity in the environment.

Notes: Soil water retention curves are needed to describe the availability of soil water to plants and to model movement of water through unsaturated soils. Measuring these characteristics is time-consuming, labour-intensive and therefore expensive. This study was conducted to develop and evaluate functions based on neural networks to predict soil water retention characteristics. Dutch and Scottish data sets were available; they contained data on 178 and 165 soil horizons, respectively. A series of three neural networks (A, B and C) was developed. Neural network A had the following input variables: topsoil, bulk density, organic matter, clay, silt and sand content. In addition neural network B had matric potential as input, and network C included soil structural data expressed as the upper and lower boundary of the ped-size class. Neural network A had three output variables: the volumetric water content at matric potentials of 0, –100 and –15 000 hPa. Both models B and C had volumetric water content, at the matric potential given as input, as output variable. The networks were tested against independent data that were extracted from the original sets of soil profiles. Accuracy of the predictions was quantified by the root of the mean squared difference (RMSE) between the measured and the predicted water contents, and the coefficient of determination (R2). For network A the RMSE varied for the three estimated water contents from 0.0264 to 0.0476 cm3 cm–3, and R2 varied from 0.80 to 0.93 for the individual model outputs. Networks B and C had an RMSE of 0.0435 and 0.0426 cm3 cm–3, respectively. For both networks, R2 was 0.89. The neural networks performed somewhat better than previous regression functions, but the improvements were not significant.

Notes: Andisols can absorb large amounts of phosphorus rapidly, and then release it slowly, yet the mechanisms by which they retain P and release it for plant growth are poorly understood. Ligand exchange of organic compounds from Al–humic complexes by P and/or Si release – due to breakdown of allophanic microstructure to provide sorption sites – might account for the retention of P, but its extent is not known. We applied a soil column flow-through technique to quantify the release of anions and organic carbon (C) associated with P sorption by two andic soils, and we related the anion release to possible mechanisms for the retention of P. Phosphate (H2PO4–, HPO42–) sorption and concurrent anion desorption were obtained by passing a 1-g P 1–1 (32 mmol KH2PO4 in 1 mm CaCl2) solution through the soil columns (25 cm3). Total dissolved P, Fe, Al, S, Ca, Mg, K, Mn, organic C and pH were determined in the eluent. Changes in eluent pH and the patterns of the retention of P and corresponding concentrations of Al, Si and organic C in the eluent were similar for the two Andisols. The general pattern and changes in pH of the eluent coincided with changes in the patterns of release of organic C and Si and the rate of P retention. Release of silica accounted for 〈 6% of the P sorbed and had only a minor role in P retention in these two Andisols. Release of organic C, however, accounted on a molar basis for 40% and 83%, respectively, of the P sorbed. Direct measurements of the pH of the eluent and release of anions and organic C concurrent to P retention contribute to rapid assessment of the controlling mechanisms of P retention. The results indirectly confirm the hypothesis of ligand exchange of solution P with organic complexes held on allophanic surfaces. The organic C release, however, is not specifically related to either the fast or the slow P retention phase. The shift in the controlling P retention reaction associated with a change from the fast to the slow P retention phase is clearly indicated by an abrupt change of the pH of the eluent. This shift, in previous studies identified graphically by a change in slope of the P sorption isotherm, can be identified directly by measuring the pH of the matrix.

Notes: Radiocaesium fixation in soils is reported to occur on frayed edge sites of micaceous minerals. The weathering of mica in acid soils may therefore influence the Cs+ fixation process and thereby the mobility of the radiopollutant. We produced a laboratory weathering model biotite → trioctahedral vermiculite → oxidized vermiculite → hydroxy interlayered vermiculite (HIV) and quantified the Cs+ fixation of each mineral both in a fixed K+–Ca2+ background and in acid conditions. The transformation process was achieved through K depletion by Na-tetraphenylboron, oxidation with Br2 and Al-intercalation using NaOH and AlCl3. In a constant K+–Ca2+ background, vermiculite fixed 92–95% of the initial 137Cs+ contamination while biotite and HIV fixed only 18–33%. In acid conditions, the interlayer occupancy by either potassium (biotite) or hydroxy-Al groups (HIV) strongly limited Cs+ fixation to 1–4% of the initial 137Cs+ contamination. Cs+ fixation occurred on vermiculitic sites associated with micaceous wedge zones. Though both oxidized and trioctahedral vermiculites fixed similar Cs+ amounts in a constant K+–Ca2+ background (92–95%), the oxidized vermiculite retained much more radiocaesium in acid conditions (78–84% against 54–59%), because of its dioctahedral character.

Notes: Soil organic matter can be intimately associated with mineral particles of various sizes. For structural studies, soil organic matter can be isolated in particle size fractions after complete dispersion of the aggregates by ultrasonication. The ultrasonic dispersion energy necessary for complete dispersion was investigated in three A and two B horizons originating from four soils differing in pedogenesis (Gleysol, Phaeozem, Podzol, Alisol), organic C (4.2–34.5 g kg–1) and clay content (24–294 g kg–1). Calorimetric calibration of five probe-type ultrasonifiers revealed that the actual energy output from an instrument can depart widely from its nominal output, and that this discrepancy varies from instrument to instrument. Calorimetric calibration is therefore essential for consistency and comparisons between laboratories. Between 450 and 500 J ml–1 of ultrasonic dispersion energy was enough to disperse completely all samples investigated. The particle size distributions obtained were close to those from standard analysis, except for smaller yields (–20 to –80 g kg–1) of sand size fractions, which suggests that dispersion by ultrasound is more effective. Based on total C, C:N ratio and distribution of dissolved C, no detachment of soil organic matter from primary organomineral complexes and no redistribution between particle size fractions could be detected in the range 30–590 J ml–1 of dispersion energy.

Notes: Understanding the natural variation of carbon within the soil, and between soil types, is crucial to improve predictive models of carbon cycling in high and mid-latitude ecosystems in response to global warming. We measured the carbon isotope distributions (12C, 13C and 14C) in soil organic matter (SOM) from Podzols, Brown Podzolic soils and Stagnohumic Gleysols from the British uplands, which were then compared with the total amounts and turnover of carbon in these soils. We did so by sampling at 2-cm intervals down six profiles of each soil type. The average amount of carbon stored in the top 28 cm of the Stagnohumic Gleysols is twice that of the other two soils. The 13C content and 14C age show a general increase with depth in all soils, and there is also a significant correlation between isotopic variation and the main pedogenic features. The latter suggests that soil-forming processes are significant in determining the carbon isotope signatures retained in SOM. Organic matter formed since 1960 is not found below 5 cm in any of the soils. Evidently organic detritus in the surface layers (LF and Oh) is rapidly mineralized. This accords with our modelled net annual C fluxes which show that more than 80% of the CO2 emanating from these soils is derived from the top 5 cm of each profile. Although these soils contain much carbon, they do not appear to assimilate and retain SOM rapidly. The mean residence time of most of their carbon is in the 2–50 years range, so the soils are fairly ineffective sinks for excess CO2 in the atmosphere. Under the predicted future ‘greenhouse’ climate, likely to favour more rapid microbial decomposition of organic materials, these soils are a potential source of CO2 and are therefore likely to accelerate global warming.

Notes: We improve two aspects of the modelling scheme for the simulation of electromagnetic radio waves, based on the Fourier pseudospectral method.When there are large contrasts in the material properties, use of the standard algorithm (regular grid) causes a series of artefacts, as, for instance, ringing and acausal events. These problems, due to the non-locality of the differential operator, are solved by using the staggered Fourier method (staggered grid).Realistic radiation patterns can be obtained from simple combinations of magnetic and electric sources. If the directivity pattern of the antenna is known, from either a finite-difference simulation or an analytic evaluation or an experimental characterization, it can then be simulated by a composite-source concept. This effective source is implemented in the modelling algorithm by means of a perturbation technique, which first computes the intensity and directional spectra of the single electromagnetic sources. Their location is optimized to obtain the best fit with a minimum number of sources. The approach is, in principle, valid for the far-field radiation pattern of the antenna.

Notes: It is well known that the migrated image given by a Kirchhoff-type (diffraction-stack) migration with limited aperture is always accompanied by some events which depend on the migration aperture. Although these events may severely affect the quality of migration, they have been studied only in 2D cases. Here, the events due to the migration aperture in 3D situations are investigated using a new method of analysing the reconstructed wavefield. It is found that a finite migration aperture results in a reconstructed wavefield with two components. One comes from the tangent points and curves between the traveltime surfaces of reflected and point-diffracted rays and is independent of the migration aperture, and the other is from the boundary of the migration aperture and depends strongly on the location and size as well as on the shape of the migration aperture. It is this last component that describes the aperture effect in migration. If the migration aperture is not sufficiently large, and if the input for migration is not zero on the boundary of the migration aperture, the boundary component may partially or totally cancel the migration signal. Furthermore, for synthetic data, the aperture effect cannot be eliminated by enlarging the migration aperture because, except for the common-shotpoint data, the aperture effect always exists however large the migration aperture becomes. This leads to the conclusion that the published Kirchhoff-type operators are not the exact inverse operators of the Fresnel–Kirchhoff integral if the input data are synthetic.

Notes: Borehole-wall imaging is currently the most reliable means of mapping discontinuities within boreholes. As these imaging techniques are expensive and thus not always included in a logging run, a method of predicting fracture frequency directly from traditional logging tool responses would be very useful and cost effective. Artificial neural networks (ANNs) show great potential in this area. ANNs are computational systems that attempt to mimic natural biological neural networks. They have the ability to recognize patterns and develop their own generalizations about a given data set. Neural networks are trained on data sets for which the solution is known and tested on data not previously seen in order to validate the network result. We show that artificial neural networks, due to their pattern recognition capabilities, are able to assess the signal strength of fracture-related heterogeneity in a borehole log and thus fracture frequency within a borehole. A combination of wireline logs (neutron porosity, bulk density, P-sonic, S-sonic, deep resistivity and shallow resistivity) were used as input parameters to the ANN. Fracture frequency calculated from borehole televiewer data was used as the single output parameter. The ANN was trained using a back-propagation algorithm with a momentum learning function. In addition to fracture frequency within a single borehole, an ANN trained on a subset of boreholes in an area could be used for prediction over the entire set of boreholes, thus allowing the lateral correlation of fracture zones.

Notes: Long-offset transient electromagnetics (LOTEM) is now regarded as a suitable electrical method for deep exploration along with magnetotellurics (MT). In this method, the vertical magnetic-field impulse response and, occasionally, the horizontal electric-field step response of a grounded-wire source on the surface of the earth are measured. Here, these two responses are computed for 3D models of three deep resistivity structures of interest in hydrocarbon exploration: (i) a faulted graben in a resistive basement rock at a depth of 4 km beneath a conductive overburden; (ii) a facies change in a resistive layer buried at a depth of 2 km in the conductive overburden above a resistive basement; and (iii) an anticlinal uplift of a resistive layer at a depth of 1 km in the conductive overburden above a resistive basement. The results show that the sensitivity of the electric-field response to model perturbation is generally greater than that of the magnetic-voltage response. Further, the electric-field sensitivity is confined to early and intermediate times while that of the magnetic-voltage response is confined to intermediate and late times. Hence it is recommended that both electric and magnetic recordings are made in a LOTEM survey so that the final results can be presented as apparent-resistivity curves derived from the two responses jointly as well as separately.

Notes: The surface nuclear magnetic resonance (SNMR) method has been tested at a site in Haldensleben, northern Germany, to assess the suitability of this new method for groundwater exploration and environmental investigations. More information is obtained by SNMR, particularly with respect to aquifer parameters, than with other geophysical techniques. SNMR measurements were carried out at three borehole locations, together with 2D and 1D direct current geoelectrics, as well as ground-penetrating radar, and well logging (induction log, gamma-ray log and pulsed neutron-gamma log). Permeabilities were calculated from the grain-size distributions of core material determined in the laboratory. It is demonstrated that the SNMR method is able to detect groundwater and the results are in good agreement with other geophysical and hydrogeological data. Using the SNMR method, the water content of the unsaturated and saturated zones (i.e. porosity of an aquifer) can be reliably determined. This information and resistivity data permit in situ determination of other aquifer parameters. Comparison of the SNMR results with borehole data clearly shows that the water content determined by SNMR is the free or mobile water in the pores. The permeabilities estimated from the SNMR decay times are similar to those derived from sieve analysis of core material. Thus, the combination of SNMR with geoelectric methods promises to be a powerful tool for studying aquifer properties.

Notes: Since the work of Postma and Backus, much has been learned about elastic constants in vertical transversely isotropic (VTI) media when the anisotropy is due to fine layering of isotropic elastic materials. Nevertheless, there has continued to be some uncertainty about the possible range of Thomsen's anisotropy parameters ε and δ for such media. We use both Monte Carlo studies and detailed analysis of Backus' equations for both two- and three-component layered media to establish the results presented. We show that ε lies in the range −3/8 〈inlineGraphic alt="leqslant R: less-than-or-eq, slant" extraInfo="nonStandardEntity" href="urn:x-wiley:00168025:GPR163:les" location="les.gif"/〉 ε 〈inlineGraphic alt="leqslant R: less-than-or-eq, slant" extraInfo="nonStandardEntity" href="urn:x-wiley:00168025:GPR163:les" location="les.gif"/〉 ½[〈v2p〉〈v−2p〉−1], for finely layered media; smaller positive and all negative values of ε occur for media with large fluctuations in the Lamé parameter λ in the component layers. We show that δ can also be either positive or negative, and that for constant density media, sign (δ) = sign (〈v−2p〉 − 〈v−2s〉〈v2s/v2p〉). Monte Carlo simulations show that among all theoretically possible random media, positive and negative δ are equally likely in finely layered media. (Of course, the δs associated with real earth materials may span some smaller subset of those that are theoretically possible, but answering this important question is beyond our present scope.) Layered media having large fluctuations in λ are those most likely to have positive δ. This is somewhat surprising since ε is often negative or a small positive number for such media, and we have the general constraint that ε − δ 〉 0 for layered VTI media. Since Gassmann's results for fluid-saturated porous media show that the mechanical effects of fluids influence only the Lamé parameter λ, not the shear modulus μ, these results suggest that small positive δ occurring together with small positive ε (but somewhat larger than δ) may be indicative of changing fluid content in a layered earth.

Notes: The multifold acquisition principle was applied to a borehole radar survey, performed in a granitic site (Grimsel Test Site, Switzerland). Two multifold coverage acquisitions (40-fold and 20-fold) were carried out in a subhorizontal borehole. Instrumental drifts (transmission time and sampling frequency fluctuations) were corrected in order to remove shifts observed on CMP gathers and to optimize velocity analysis and trace stacking. Computation of velocity spectra was adapted in order to take into account the features of the medium investigated (homogeneous velocity, various reflector orientations). The NMO velocities were then interpreted as angles between reflectors and the survey line. The processing, based on the computation of several constant velocity stacked sections performed with different NMO velocities, leads to better results than the standard DMO + NMO processing. The signal-to-noise ratio of the stacked profile is improved in comparison with the single-fold section, which results from a standard acquisition. From a practical point of view, the implementation of a multifold radar survey within a borehole is difficult but a greater investigation range is obtained, more reflectors are detected and the mapping of geological discontinuities is improved.

Notes: Electromagnetic (EM) fields radiated from the transmission lines of the New Zealand electricity grid have been digitally recorded at test sites near the Tokaanu geothermal field. Amplitudes and phases of the 50 Hz signals (and the odd harmonics up to 450 Hz) were determined using a software implementation of a phase-locked filter. These data were then analysed to determine the components of the magnetotelluric impedance tensor and the corresponding apparent resistivities and phases. At most sites, there was sufficient variation in the elliptical polarization of the EM fields to enable the impedance tensors to be determined in full. Sites where the EM data had been affected by near-source effects were identified by having large vertical magnetic field components and by being closer to a power line source than about 3–5 skin depths. With the test measurements, the north-eastern part of the Tokaanu geothermal field was successfully delineated giving low resistivities (〈 5 Ωm) on the inside and higher resistivities on the outside, in agreement with the Schlumberger array DC apparent resistivities. The small size of the 50 Hz magnetotelluric equipment and its portable nature make this method of resistivity measurement suitable for reconnaissance resistivity mapping in places with difficult access.

Notes: Parameter estimation from the elliptical variations in the normal-moveout (NMO) velocity in azimuthally anisotropic media is sensitive to the angular separation between the survey lines in 2D, or equivalently, the source-to-receiver azimuth in 3D, and to the set of azimuths used in the inversion procedure. The accuracy in estimating the orientation of an NMO ellipse, in particular the parameter α, is also sensitive to the magnitude of anisotropy. On the other hand, the accuracy in estimating the semi-axes of the NMO-velocity ellipse is about the same for any magnitude of anisotropy.  To invert for the NMO ellipse parameters at least three NMO-velocity measurements along distinct azimuth directions are needed. In order to maximize the accuracy and stability in parameter estimation, it is best to have the azimuths for the three source-to-receiver directions 60° apart. Having more than three distinct source-to-receiver azimuths (e.g. full azimuthal coverage) provides a useful data redundancy that enhances the quality of the estimates.  In order to maximize quality in the inversion process, it is recommended to design the seismic data acquisition such that it contains small sectors (≤10°) with adequate fold and offset distribution.  Using three NMO-velocity measurements, 60° apart, an azimuthally anisotropic layer overlain by an azimuthally isotropic overburden (as might occur for fractured reservoirs) should have a relative thickness (in time) with respect to the total thickness at least equal to the ratio of the error in the NMO (stacking) velocity to the interval anisotropy of the fractured layer. Coverage along more than three azimuths, however, improves this limitation, which is imposed by Dix differentiation, by at most 50%, depending on the number of observations (NMO velocities) that enter the inversion procedure.

Notes: Obtaining high-resolution images of the geology and hydrogeology of the subsurface in the depth range from ground level to 50 m is one of the major challenges of modern geophysics. The methods which are commonly used (such as compressional-wave surveys and ground-penetrating radar) often suffer from adverse effects caused by the near-surface conditions, changes in water saturation and various sources of noise. This paper demonstrates some of the advantages offered by the use of shear-wave seismology and by the combination of shear- and compressional-wave seismic methods in shallow subsurface investigations.Multicomponent shallow seismic tests were carried out at four different sites to examine the effectiveness of different acquisition geometries under a variety of near-surface geological conditions. Near-surface conditions encountered at the sites included thick clays, clay/sand sequences overlying Chalk, mudstone overlying granodiorite bedrock and landfill material.Under all conditions, shear-wave data acquisition was found to have advantages over compressional-wave acquisition for the investigation of the shallow subsurface. Shear head waves, being unaffected by water saturation, achieved penetration to greater depths at a site in Crewkerne, Dorset where compressional head-wave penetration was limited to the near-surface layers. Better vertical resolution was achieved at shallow depths using shear-wave reflection energy at a landfill site. Shear-wave reflections from shallow interfaces were in some cases less affected by noise compared with the equivalent compressional-wave reflections. Combinations of shear- and compressional-wave data recording allowed the measurement of a Poisson's ratio log and gave indications of seismic anisotropy at two sites where dipping clay layers were present.

Notes: A laboratory study was carried out to investigate the influence of confining stress on compressional- and shear-wave velocities for a set of rock samples from gas-producing sandstone reservoirs in the Cooper Basin, South Australia. The suite of samples consists of 22 consolidated sublitharenites with helium porosity ranging from 2.6% to 16.6%. We used a pulse-echo technique to measure compressional- and shear-wave velocities on dry samples (cylindrical 4.6 × 2 cm) at room temperature and at elevated confining stress (≤ 60 MPa). Compressional- and shear-wave velocities in samples increase non-linearly with confining stress. A regression equation of the form V = A − Be−DP gives a good fit to the measured velocities with improved prediction of velocities at high confining stresses compared with equations suggested by other studies. The predicted microcrack-closure stresses of the samples show values ranging from 70 MPa to 95 MPa and insignificant correlation with porosity, permeability or clay content. There is a positive correlation between change in velocity with core porosity and permeability, but this association is weak and diminishes with increasing confining stress. Experimental results show that pore geometry, grain-contact type, and distribution and location of clay particles may be more significant than total porosity and clay content in describing the stress sensitivity of sandstones at in situ reservoir effective stress. The stress dependence of Cooper Basin sandstones is very large compared with data from other studies. The implication of our study for hydrocarbon exploration is that where the in situ reservoir effective stress is much less than the microcrack-closure stress of the reservoir rocks, the variation of reservoir effective stress could cause significant changes in velocity of the reservoir rocks. The velocity changes induced by effective stress in highly stress-sensitive rocks can be detected at sonic-log and probably surface-seismic frequencies.

Notes: The applicability of seismic refraction profiling for the detection of velocity inversion, which is also known as a low-velocity layer (LVL), is investigated with the aid of synthetic seismogram computations for a range of models. Our computational models focus on the inherent ambiguities in the interpretation of first-arrival time delays or ‘skips’ in terms of LVL model parameters. The present modelling results reveal that neither the measure nor even the existence of a shadow zone and/or a time shift (skip) in first arrivals is necessarily indicative of an LVL. Besides attenuation effects, the cap-layer velocity gradient is a critical parameter, determining the termination point of the cap-layer diving wave and thus the time skip.We suggest that shallow LVLs can be delineated more reliably by traveltime and amplitude modelling of coherent phases reflected from their top and bottom boundaries, often clearly observed in the pre- and near-critical ranges in seismogram sections of refraction profiling experiments with a close receiver spacing. We demonstrate the applicability of this approach for a field data set of a refraction profile in the West Bengal Basin, India. The inferred LVL corresponds to the Gondwana sediments underlying the higher-velocity layer of the Rajmahal Traps. This interpretation is consistent with the data from a nearby well in the region.

Notes: The determination of small-scale velocity anomalies (from tens to a few hundreds of metres) is a major problem in seismic exploration. The impact of such anomalies on a structural interpretation can be dramatic and conventional techniques such as tomographic inversion or migration velocity analysis are powerless to resolve the ambiguity between structural and velocity origins of anomalies. We propose an alternative approach based on stochastic modelling of numerous anomalies until a set of models is found which can explain the real data. This technique attempts to include as much a priori geological information as possible. It aims at providing the interpreter with a set of velocity anomalies which could possibly be responsible for the structural response. The interpreter can then choose one or several preferred models and pursue a more sophisticated analysis. The class of retained models are all equivalent in terms of data and therefore represent the uncertainty in the model space.The procedure emulates the real processing sequence using a simplified scheme. Essentially, the technique consists of five steps: 1 Interpretation of a structural anomaly in terms of a velocity anomaly with its possible variations in terms of position, size and amplitude. 2 Drawing a model by choosing the parameters of the anomaly within the acceptable range. 3 Modelling the traveltimes in this model and producing the imaging of the reflected interface. 4 Comparing the synthetic data with the real data and keeping the model if it lies within the data uncertainty range. 5 Iterate from step 2.In order to avoid the high computational cost inherent in using statistical determinations, simplistic assumptions have been made:• The anomaly is embedded in a homogeneous medium: we assume that the refraction and the time shift due to the anomaly have a first-order effect compared with ray bending in the intermediate layers.• We model only the zero-offset rays and therefore we restrict ourselves to structural problems.• We simulate time migration and so address only models of limited structural complexity.These approximations are justified in a synthetic model which includes strong lateral velocity variations, by comparing the result of a full processing sequence (prestack modelling, stack and depth migration) with the simplified processing. This model is then used in a blind test on the inversion scheme.

Notes: The relative merits of a 48-gun, 9324 cu. in. (153 litre) airgun array and a 200 kg explosive source are considered for the purposes of long-range (0–400 km) refraction seismic work, with particular reference to traveltime modelling. Theoretical source calculations indicate that in the frequency range 2.5–12.0 Hz, the airgun source will produce an RMS pressure ∼ 8% of that produced by the explosive source and an initial burst pressure ∼17% of that produced by the explosive source. Observed data support these calculations at short ranges and illustrate the greater attenuation of the airgun signal with range due to its lack of very low frequency (〈 5 Hz) content. At short offsets, the airgun array provides a preferable seismic source to the explosives, due to densely spaced shots and a consistent waveform resulting in excellent trace-to-trace coherence. With increasing offsets, it may be necessary to stack the airgun data to enhance its signal-to-noise ratio: here we use a 4-fold stack. Large explosive shots, although more powerful, produce a less consistent waveform and are more widely spaced due to operational constraints. The offset at which airguns provide a preferable source is dependent on the ambient noise. This practical comparison of real sources demonstrates that, even without advanced processing, a well-tuned airgun array may provide a preferable source to explosives at offsets up to 160 km, under favourable experimental conditions.