Predicted soil organic carbon stocks and changes in Kenya between 1990 and 2030

Under the United Nations Framework Convention on Climate Change (UNFCCC), Non-Annex 1 countries such as Kenya are obliged to report green house gas (GHG) emissions from all sources where possible, including those from soils as a result of changes in land use or land management. At present, the convention encourages countries to estimate emissions using the most advanced methods possible, given the country circumstances and resources. Estimates of soil organic carbon (SOC) stocks and changes were made for Kenya using the Global Environment Facility Soil Organic Carbon (GEFSOC) Modelling System. The tool conducts analysis using three methods: (1) the Century general ecosystem model; (2) the RothC soil C decomposition model; and (3) the Intergovernmental Panel on Climate Change (IPCC) method for assessing soil C at regional scales. The required datasets included: land use history, monthly mean precipitation, monthly mean minimum and maximum temperatures for all the agro-climatic zones of Kenya and historical vegetation cover. Soil C stocks of 1.4-2.0 Pg (0-20 cm), compared well with a Soil and Terrain (SOTER) based approach that estimated similar to .8-2.0 Pg (0-30 cm). In 1990 48% of the country had SOC stocks of < 18 t C ha(-1) and 20% of the country had SOC stocks of 18-30 t C ha(-1), whereas in 2000 56% of the country had SOC stocks of < 18 t C ha(-1) and 31% of the country had SOC stocks of 18-30 t C ha(-1). Conversion of natural vegetation to annual crops led to the greatest soil C losses. Simulations suggest that soil C losses remain substantial throughout the modelling period of 1990-2030. All three methods involved in the GEFSOC System estimated that there would be a net loss of soil C between 2000 and 2030 in Kenya. The decline was more marked with RothC than with Century or the IPCC method. In non-hydric soils the SOC change rates were more pronounced in high sandy soils compared to high clay soils in most land use systems. (C) 2007 Elsevier B.V. All rights reserved.

Evaluation of two soil carbon models using two Kenyan long term experimental datasets

RothC and Century are two of the most widely used soil organic matter (SOM) models. However there are few examples of specific parameterisation of these models for environmental conditions in East Africa. The aim of this study was therefore, to evaluate the ability of RothC and the Century to estimate changes in soil organic carbon (SOC) resulting from varying land use/management practices for the climate and soil conditions found in Kenya. The study used climate, soils and crop data from a long term experiment (1976-2001) carried out at The Kabete site at The Kenya National Agricultural Research Laboratories (NARL, located in a semi-humid region) and data from a 13 year experiment carried out in Machang'a (Embu District, located in a semi-arid region). The NARL experiment included various fertiliser (0, 60 and 120 kg of N and P2O5 ha(-1)), farmyard manure (FYM - 5 and 10 t ha(-1)) and plant residue treatments, in a variety of combinations. The Machang'a experiment involved a fertiliser (51 kg N ha(-1)) and a FYM (0, 5 and 10 t ha(-1)) treatment with both monocropping and intercropping. At Kabete both models showed a fair to good fit to measured data, although Century simulations for treatments with high levels of FYM were better than those without. At the Machang'a site with monocrops, both models showed a fair to good fit to measured data for all treatments. However, the fit of both models (especially RothC) to measured data for intercropping treatments at Machang'a was much poorer. Further model development for intercrop systems is recommended. Both models can be useful tools in soil C Predictions, provided time series of measured soil C and crop production data are available for validating model performance against local or regional agricultural crops. (C) 2007 Elsevier B.V. All rights reserved.

Indicators of nitrate in wetland surface and soil-waters: interactions of vegetation and environmental factors

This paper describes a new bio-indicator method for assessing wetland ecosystem health: as such, the study is particularly relevant to current legislation such as the EU Water Framework Directive, which provides a baseline of the current status Of Surface waters. Seven wetland sites were monitored across northern Britain, with model construction data for predicting, eco-hydroloplical relationships collected from five sites during 1999, Two new sites and one repeat site were monitored during 2000 to provide model test data. The main growing season for the vegetation, and hence the sampling period, was May-August during both years. Seasonal mean concentrations of nitrate (NO3-) in surface and soil water samples during 1999 ranged from 0.01 to 14.07 mg N 1(-1), with a mean value of 1.01 mg N 1(-1). During 2000, concentrations ranged from trace level (<0.01 m- N 1(-1)) to 9.43 mg N 1(-1), with a mean of 2.73 mg N 1(.)(-1) Surface and soil-water nitrate concentrations did not influence plant species composition significantly across representative tall herb fen and mire communities. Predictive relationships were found between nitrate concentrations and structural characteristics of the wetland vegetation, and a model was developed which predicted nitrate concentrations from measures of plant diversity, canopy structure and density of reproductive structures. Two further models, which predicted stem density and density of reproductive structures respectively, utilised nitrate concentration as one of the independent predictor variables. Where appropriate, the models were tested using data collected during 2000. This approach is complementary to species-based monitoring, representing a useful and simple too] to assess ecological status in target wetland systems and has potential for bio-indication purposes.

Comparing sampling needs for variograms of soil properties computed by the method of moments and residual maximum likelihood

It has been generally accepted that the method of moments (MoM) variogram, which has been widely applied in soil science, requires about 100 sites at an appropriate interval apart to describe the variation adequately. This sample size is often larger than can be afforded for soil surveys of agricultural fields or contaminated sites. Furthermore, it might be a much larger sample size than is needed where the scale of variation is large. A possible alternative in such situations is the residual maximum likelihood (REML) variogram because fewer data appear to be required. The REML method is parametric and is considered reliable where there is trend in the data because it is based on generalized increments that filter trend out and only the covariance parameters are estimated. Previous research has suggested that fewer data are needed to compute a reliable variogram using a maximum likelihood approach such as REML, however, the results can vary according to the nature of the spatial variation. There remain issues to examine: how many fewer data can be used, how should the sampling sites be distributed over the site of interest, and how do different degrees of spatial variation affect the data requirements? The soil of four field sites of different size, physiography, parent material and soil type was sampled intensively, and MoM and REML variograms were calculated for clay content. The data were then sub-sampled to give different sample sizes and distributions of sites and the variograms were computed again. The model parameters for the sets of variograms for each site were used for cross-validation. Predictions based on REML variograms were generally more accurate than those from MoM variograms with fewer than 100 sampling sites. A sample size of around 50 sites at an appropriate distance apart, possibly determined from variograms of ancillary data, appears adequate to compute REML variograms for kriging soil properties for precision agriculture and contaminated sites. (C) 2007 Elsevier B.V. All rights reserved.

The analysis of ranked observations of soil structure using indicator geostatistics

Structure is an important physical feature of the soil that is associated with water movement, the soil atmosphere, microorganism activity and nutrient uptake. A soil without any obvious organisation of its components is known as apedal and this state can have marked effects on several soil processes. Accurate maps of topsoil and subsoil structure are desirable for a wide range of models that aim to predict erosion, solute transport, or flow of water through the soil. Also such maps would be useful to precision farmers when deciding how to apply nutrients and pesticides in a site-specific way, and to target subsoiling and soil structure stabilization procedures. Typically, soil structure is inferred from bulk density or penetrometer resistance measurements and more recently from soil resistivity and conductivity surveys. To measure the former is both time-consuming and costly, whereas observations made by the latter methods can be made automatically and swiftly using a vehicle-mounted penetrometer or resistivity and conductivity sensors. The results of each of these methods, however, are affected by other soil properties, in particular moisture content at the time of sampling, texture, and the presence of stones. Traditional methods of observing soil structure identify the type of ped and its degree of development. Methods of ranking such observations from good to poor for different soil textures have been developed. Indicator variograms can be computed for each category or rank of structure and these can be summed to give the sum of indicator variograms (SIV). Observations of the topsoil and subsoil structure were made at four field sites where the soil had developed on different parent materials. The observations were ranked by four methods and indicator and the sum of indicator variograms were computed and modelled for each method of ranking. The individual indicators were then kriged with the parameters of the appropriate indicator variogram model to map the probability of encountering soil with the structure represented by that indicator. The model parameters of the SIVs for each ranking system were used with the data to krige the soil structure classes, and the results are compared with those for the individual indicators. The relations between maps of soil structure and selected wavebands from aerial photographs are examined as basis for planning surveys of soil structure. (C) 2007 Elsevier B.V. All rights reserved.

Effect of manure application on crop yield and soil chemical properties in a long-term field trial of semi-arid Kenya

The sustainability of cereal/legume intercropping was assessed by monitoring trends in grain yield, soil organic C (SOC) and soil extractable P (Olsen method) measured over 13 years at a long-term field trial on a P-deficient soil in semi-arid Kenya. Goat manure was applied annually for 13 years at 0, 5 and 10 t ha(-1) and trends in grain yield were not identifiable because of season-to-season variations. SOC and Olsen P increased for the first seven years of manure application and then remained constant. The residual effect of manure applied for four years only lasted another seven to eight years when assessed by yield, SOC and Olsen P. Mineral fertilizers provided the same annual rates of N and P as in 5 t ha(-1) manure and initially ,gave the same yield as manure, declining after nine years to about 80%. Therefore, manure applications could be made intermittently and nutrient requirements topped-up with fertilizers. Grain yields for sorghum with continuous manure were described well by correlations with rainfall and manure input only, if data were excluded for seasons with over 500 mm rainfall. A comprehensive simulation model should correctly describe crop losses caused by excess water.

Sensitivity analysis of a catchment-scale nitrogen model

There are now considerable expectations that semi-distributed models are useful tools for supporting catchment water quality management. However, insufficient attention has been given to evaluating the uncertainties inherent to this type of model, especially those associated with the spatial disaggregation of the catchment. The Integrated Nitrogen in Catchments model (INCA) is subjected to an extensive regionalised sensitivity analysis in application to the River Kennet, part of the groundwater-dominated upper Thames catchment, UK The main results are: (1) model output was generally insensitive to land-phase parameters, very sensitive to groundwater parameters, including initial conditions, and significantly sensitive to in-river parameters; (2) INCA was able to produce good fits simultaneously to the available flow, nitrate and ammonium in-river data sets; (3) representing parameters as heterogeneous over the catchment (206 calibrated parameters) rather than homogeneous (24 calibrated parameters) produced a significant improvement in fit to nitrate but no significant improvement to flow and caused a deterioration in ammonium performance; (4) the analysis indicated that calibrating the flow-related parameters first, then calibrating the remaining parameters (as opposed to calibrating all parameters together) was not a sensible strategy in this case; (5) even the parameters to which the model output was most sensitive suffered from high uncertainty due to spatial inconsistencies in the estimated optimum values, parameter equifinality and the sampling error associated with the calibration method; (6) soil and groundwater nutrient and flow data are needed to reduce. uncertainty in initial conditions, residence times and nitrogen transformation parameters, and long-term historic data are needed so that key responses to changes in land-use management can be assimilated. The results indicate the general, difficulty of reconciling the questions which catchment nutrient models are expected to answer with typically limited data sets and limited knowledge about suitable model structures. The results demonstrate the importance of analysing semi-distributed model uncertainties prior to model application, and illustrate the value and limitations of using Monte Carlo-based methods for doing so. (c) 2005 Elsevier B.V. All rights reserved.

Moving towards a more mechanistic approach in the determination of soil heat flux from remote measurements - I. A universal approach to calculate thermal inertia

In this paper we pledge that physically based equations should be combined with remote sensing techniques to enable a more theoretically rigorous estimation of area-average soil heat flux, G. A standard physical equation (i.e. the analytical or exact method) for the estimation of G, in combination with a simple, but theoretically derived, equation for soil thermal inertia (F), provides the basis for a more transparent and readily interpretable method for the estimation of G; without the requirement for in situ instrumentation. Moreover, such an approach ensures a more universally applicable method than those derived from purely empirical studies (employing vegetation indices and albedo, for example). Hence, a new equation for the estimation of Gamma(for homogeneous soils) is discussed in this paper which only requires knowledge of soil type, which is readily obtainable from extant soil databases and surveys, in combination with a coarse estimate of moisture status. This approach can be used to obtain area-averaged estimates of Gamma(and thus G, as explained in paper II) which is important for large-scale energy balance studies that employ aircraft or satellite data. Furthermore, this method also relaxes the instrumental demand for studies at the plot and field scale (no requirement for in situ soil temperature sensors, soil heat flux plates and/or thermal conductivity sensors). In addition, this equation can be incorporated in soil-vegetation-atmosphere-transfer models that use the force restore method to update surface temperatures (such as the well-known ISBA model), to replace the thermal inertia coefficient.

Adaptation of the Integrated Nitrogen Model for Catchments (INCA) to seasonally snow-covered catchments

Testing of the Integrated Nitrogen model for Catchments (INCA) in a wide range of ecosystem types across Europe has shown that the model underestimates N transformation processes to a large extent in northern catchments of Finland and Norway in winter and spring. It is found, and generally assumed, that microbial activity in soils proceeds at low rates at northern latitudes during winter, even at sub-zero temperatures. The INCA model was modified to improve the simulation of N transformation rates in northern catchments, characterised by cold climates and extensive snow accumulation and insulation in winter, by introducing an empirical function to simulate soil temperatures below the seasonal snow pack, and a degree-day model to calculate the depth of the snow pack. The proposed snow-correction factor improved the simulation of soil temperatures at Finnish and Norwegian field sites in winter, although soil temperature was still underestimated during periods with a thin snow cover. Finally, a comparison between the modified INCA version (v. 1.7) and the former version (v. 1.6) was made at the Simojoki river basin in northern Finland and at Dalelva Brook in northern Norway. The new modules did not imply any significant changes in simulated NO3- concentration levels in the streams but improved the timing of simulated higher concentrations. The inclusion of a modified temperature response function and an empirical snow-correction factor improved the flexibility and applicability of the model for climate effect studies.

Simulating current global river runoff with a global hydrological model: model revisions, validation, and sensitivity analysis

Global hydrological models (GHMs) model the land surface hydrologic dynamics of continental-scale river basins. Here we describe one such GHM, the Macro-scale - Probability-Distributed Moisture model.09 (Mac-PDM.09). The model has undergone a number of revisions since it was last applied in the hydrological literature. This paper serves to provide a detailed description of the latest version of the model. The main revisions include the following: (1) the ability for the model to be run for n repetitions, which provides more robust estimates of extreme hydrological behaviour, (2) the ability of the model to use a gridded field of coefficient of variation (CV) of daily rainfall for the stochastic disaggregation of monthly precipitation to daily precipitation, and (3) the model can now be forced with daily input climate data as well as monthly input climate data. We demonstrate the effects that each of these three revisions has on simulated runoff relative to before the revisions were applied. Importantly, we show that when Mac-PDM.09 is forced with monthly input data, it results in a negative runoff bias relative to when daily forcings are applied, for regions of the globe where the day-to-day variability in relative humidity is high. The runoff bias can be up to - 80% for a small selection of catchments but the absolute magnitude of the bias may be small. As such, we recommend future applications of Mac-PDM.09 that use monthly climate forcings acknowledge the bias as a limitation of the model. The performance of Mac-PDM.09 is evaluated by validating simulated runoff against observed runoff for 50 catchments. We also present a sensitivity analysis that demonstrates that simulated runoff is considerably more sensitive to method of PE calculation than to perturbations in soil moisture and field capacity parameters.

Folsomia candida (Collembola): A "standard" soil arthropod

Folsomia candida Willem 1902, a member of the order Collembola (colloquially called springtails), is a common and widespread arthropod that occurs in soils throughout the world. The species is parthenogenetic and is easy to maintain in the laboratory on a diet of granulated dry yeast. F. candida has been used as a "standard" test organism for more than 40 years for estimating the effects of pesticides and environmental pollutants on nontarget soil arthropods. However. it has also been employed as a model for the investigation of numerous other phenomena such as cold tolerance, quality as a prey item, and effects of microarthropod grazing on pathogenic fungi and mycorrhizae of plant roots. In this comprehensive review. aspects of the life history, ecology, and ecotoxicology of F candida are covered. We focus on the recent literature, especially studies that have examined the effects of soil pollutants on reproduction in F candida using the protocol published by the International Standards Organization in 1999.

Development and assessment of a coupled crop-climate model

It is well established that crop production is inherently vulnerable to variations in the weather and climate. More recently the influence of vegetation on the state of the atmosphere has been recognized. The seasonal growth of crops can influence the atmosphere and have local impacts on the weather, which in turn affects the rate of seasonal crop growth and development. Considering the coupled nature of the crop-climate system, and the fact that a significant proportion of land is devoted to the cultivation of crops, important interactions may be missed when studying crops and the climate system in isolation, particularly in the context of land use and climate change. To represent the two-way interactions between seasonal crop growth and atmospheric variability, we integrate a crop model developed specifically to operate at large spatial scales (General Large Area Model for annual crops) into the land surface component of a global climate model (GCM; HadAM3). In the new coupled crop-climate model, the simulated environment (atmosphere and soil states) influences growth and development of the crop, while simultaneously the temporal variations in crop leaf area and height across its growing season alter the characteristics of the land surface that are important determinants of surface fluxes of heat and moisture, as well as other aspects of the land-surface hydrological cycle. The coupled model realistically simulates the seasonal growth of a summer annual crop in response to the GCM's simulated weather and climate. The model also reproduces the observed relationship between seasonal rainfall and crop yield. The integration of a large-scale single crop model into a GCM, as described here, represents a first step towards the development of fully coupled crop and climate models. Future development priorities and challenges related to coupling crop and climate models are discussed.

Influence of vegetation on the local climate and hydrology in the tropics: sensitivity to soil parameters

Land use change with accompanying major modifications to the vegetation cover is widespread in the tropics, due to increasing demands for agricultural land, and may have significant impacts on the climate. This study investigates (1) the influence of vegetation on the local climate in the tropics; (2) how that influence varies from region to region; and (3) how the sensitivity of the local climate to vegetation, and hence land use change, depends on the hydraulic characteristics of the soil. A series of idealised experiments with the Hadley Centre atmospheric model, HadAM3, are described in which the influence of vegetation in the tropics is assessed by comparing the results of integrations with and without tropical vegetation. The sensitivity of the results to the soil characteristics is then explored by repeating the experiments with a differing, but equally valid, description of soil hydraulic parameters. The results have shown that vegetation has a significant moderating effect on the climate throughout the tropics by cooling the surface through enhanced latent heat fluxes. The influence of vegetation is, however, seasonally dependent, with much greater impacts during the dry season when the availability of surface moisture is limited. Furthermore, there are significant regional variations both in terms of the magnitude of the cooling and in the response of the precipitation. Not all regions show a feedback of vegetation on the local precipitation; this result has been related both to vegetation type and to the prevailing meteorological conditions. An important finding has been the sensitivity of the results to the specification of the soil hydraulic parameters. The introduction of more freely draining soils has changed the soil-moisture contents of the control, vegetated system and has reduced, significantly, the climate sensitivity to vegetation and by implication, land use change. Changes to the soil parameters have also had an impact on the soil hydrology and its interaction with vegetation, by altering the partitioning between fast and slow runoff processes. These results raise important questions about the representation of highly heterogeneous soil characteristics in climate models, as well as the potential influence of land use change on the soil characteristics themselves.

Vegetation dynamics in abandoned crop fields on a Mediterranean island: Development of succession model and estimation of disturbance thresholds

This case study on the Sifnos island, Greece, assesses the main factors controlling vegetation succession following crop abandonment and describes the vegetation dynamics of maquis and phrygana formations in relation to alternative theories of secondary succession. Field survey data were collected and analysed at community as well as species level. The results show that vegetation succession on abandoned crop fields is determined by the combined effects of grazing intensity, soil and geological characteristics and time. The analysis determines the quantitative grazing thresholds that modify the successional pathway. Light grazing leads to dominance by maquis vegetation while overgrazing leads to phryganic vegetation. The proposed model shows that vegetation succession following crop abandonment is a complex multi-factor process where the final or the stable stage of the process is not predefined but depends on the factors affecting succession. An example of the use of succession models and disturbance thresholds as a policy assessment tool is presented by evaluating the likely vegetation impacts of the recent reform of the Common Agricultural Policy on Sifnos island over a 20-30-year time horizon. (c) 2006 Elsevier B.V. All rights reserved.

Folsomia candida (Collembola): A "standard" soil arthropod

Folsomia candida Willem 1902, a member of the order Collembola (colloquially called springtails), is a common and widespread arthropod that occurs in soils throughout the world. The species is parthenogenetic and is easy to maintain in the laboratory on a diet of granulated dry yeast. F. candida has been used as a "standard" test organism for more than 40 years for estimating the effects of pesticides and environmental pollutants on nontarget soil arthropods. However. it has also been employed as a model for the investigation of numerous other phenomena such as cold tolerance, quality as a prey item, and effects of microarthropod grazing on pathogenic fungi and mycorrhizae of plant roots. In this comprehensive review. aspects of the life history, ecology, and ecotoxicology of F candida are covered. We focus on the recent literature, especially studies that have examined the effects of soil pollutants on reproduction in F candida using the protocol published by the International Standards Organization in 1999.