A simple model for predicting soil temperature in snow-covered and seasonally frozen soil: model description and testing

Microbial processes in soil are moisture, nutrient and temperature dependent and, consequently, accurate calculation of soil temperature is important for modelling nitrogen processes. Microbial activity in soil occurs even at sub-zero temperatures so that, in northern latitudes, a method to calculate soil temperature under snow cover and in frozen soils is required. This paper describes a new and simple model to calculate daily values for soil temperature at various depths in both frozen and unfrozen soils. The model requires four parameters average soil thermal conductivity, specific beat capacity of soil, specific heat capacity due to freezing and thawing and an empirical snow parameter. Precipitation, air temperature and snow depth (measured or calculated) are needed as input variables. The proposed model was applied to five sites in different parts of Finland representing different climates and soil types. Observed soil temperatures at depths of 20 and 50 cm (September 1981-August 1990) were used for model calibration. The calibrated model was then tested using observed soil temperatures from September 1990 to August 2001. R-2-values of the calibration period varied between 0.87 and 0.96 at a depth of 20 cm and between 0.78 and 0.97 at 50 cm. R-2 -values of the testing period were between 0.87 and 0.94 at a depth of 20cm. and between 0.80 and 0.98 at 50cm. Thus, despite the simplifications made, the model was able to simulate soil temperature at these study sites. This simple model simulates soil temperature well in the uppermost soil layers where most of the nitrogen processes occur. The small number of parameters required means, that the model is suitable for addition to catchment scale models.

Incorporation of water vapour transfer in the JULES Land Surface Model: implications for key soil variables and land surface fluxes

This study focuses on the mechanisms underlying water and heat transfer in upper soil layers, and their effects on soil physical prognostic variables and the individual components of the energy balance. The skill of the JULES (Joint UK Land Environment Simulator) land surface model (LSM) to simulate key soil variables, such as soil moisture content and surface temperature, and fluxes such as evaporation, is investigated. The Richards equation for soil water transfer, as used in most LSMs, was updated by incorporating isothermal and thermal water vapour transfer. The model was tested for three sites representative of semi-arid and temperate arid climates: the Jornada site (New Mexico, USA), Griffith site (Australia) and Audubon site (Arizona, USA). Water vapour flux was found to contribute significantly to the water and heat transfer in the upper soil layers. This was mainly due to isothermal vapour diffusion; thermal vapour flux also played a role at the Jornada site just after rainfall events. Inclusion of water vapour flux had an effect on the diurnal evolution of evaporation, soil moisture content and surface temperature. The incorporation of additional processes, such as water vapour flux among others, into LSMs may improve the coupling between the upper soil layers and the atmosphere, which in turn could increase the reliability of weather and climate predictions.

An easily implemented agro-hydrological procedure with dynamic root simulation for water transfer in the crop-soil system: validation and application

Models for water transfer in the crop-soil system are key components of agro-hydrological models for irrigation, fertilizer and pesticide practices. Many of the hydrological models for water transfer in the crop-soil system are either too approximate due to oversimplified algorithms or employ complex numerical schemes. In this paper we developed a simple and sufficiently accurate algorithm which can be easily adopted in agro-hydrological models for the simulation of water dynamics. We used a dual crop coefficient approach proposed by the FAO for estimating potential evaporation and transpiration, and a dynamic model for calculating relative root length distribution on a daily basis. In a small time step of 0.001 d, we implemented algorithms separately for actual evaporation, root water uptake and soil water content redistribution by decoupling these processes. The Richards equation describing soil water movement was solved using an integration strategy over the soil layers instead of complex numerical schemes. This drastically simplified the procedures of modeling soil water and led to much shorter computer codes. The validity of the proposed model was tested against data from field experiments on two contrasting soils cropped with wheat. Good agreement was achieved between measurement and simulation of soil water content in various depths collected at intervals during crop growth. This indicates that the model is satisfactory in simulating water transfer in the crop-soil system, and therefore can reliably be adopted in agro-hydrological models. Finally we demonstrated how the developed model could be used to study the effect of changes in the environment such as lowering the groundwater table caused by the construction of a motorway on crop transpiration. (c) 2009 Elsevier B.V. All rights reserved.

Effects of carbon dioxide on the searching behaviour of the root-feeding clover weevil Sitona lepidus (Coleoptera : Curculionidae)

The respiratory emission of CO2 from roots is frequently proposed as an attractant that allows soil-dwelling insects to locate host plant roots, but this role has recently become less certain. CO2 is emitted from many sources other than roots, so does not necessarily indicate the presence of host plants, and because of the high density of roots in the upper soil layers, spatial gradients may not always be perceptible by soil-dwelling insects. The role of CO2 in host location was investigated using the clover root weevil Sitona lepidus Gyllenhall and its host plant white clover (Trifolium repens L.) as a model system. Rhizochamber experiments showed that CO2 concentrations were approximately 1000 ppm around the roots of white clover, but significantly decreased with increasing distance from roots. In behavioural experiments, no evidence was found for any attraction by S. lepidus larvae to point emissions of CO2, regardless of emission rates. Fewer than 15% of larvae were attracted to point emissions of CO2, compared with a control response of 17%. However, fractal analysis of movement paths in constant CO2 concentrations demonstrated that searching by S. lepidus larvae significantly intensified when they experienced CO2 concentrations similar to those found around the roots of white clover (i.e. 1000 ppm). It is suggested that respiratory emissions of CO2 may act as a 'search trigger' for S. lepidus, whereby it induces larvae to search a smaller area more intensively, in order to detect location cues that are more specific to their host plant.

Factors affecting distribution of vegetation types on abandoned cropland in the hilly-gullied Loess Plateau region of China

A study was conducted in the forest-steppe region of the Loess Plateau to provide insight into the factors affecting the process of vegetation establishment, and to provide recommendations for the selection of indigenous species in order to speed up the succession process and to allow the establishment of vegetation more resistant to soil erosion. Four distinctive vegetation types were identified, and their distribution was affected not only by the time since abandonment but also by other environmental factors, mainly soil water and total P in the upper soil layers. One of the vegetation types, dominated by Artemisia scoparia, formed the early successional stage after abandonment while the other three types formed later successional stages with their distribution determined by the soil water content and total P. It can be concluded that the selection of appropriate species for introduction to accelerate succession should be determined by the local conditions and especially the total P concentration and soil water content.

The effects of scene heterogeneity on soil moisture retrieval from passive microwave data

The s–x model of microwave emission from soil and vegetation layers is widely used to estimate soil moisture content from passive microwave observations. Its application to prospective satellite-based observations aggregating several thousand square kilometres requires understanding of the effects of scene heterogeneity. The effects of heterogeneity in soil surface roughness, soil moisture, water area and vegetation density on the retrieval of soil moisture from simulated single- and multi-angle observing systems were tested. Uncertainty in water area proved the most serious problem for both systems, causing errors of a few percent in soil moisture retrieval. Single-angle retrieval was largely unaffected by the other factors studied here. Multiple-angle retrievals errors around one percent arose from heterogeneity in either soil roughness or soil moisture. Errors of a few percent were caused by vegetation heterogeneity. A simple extension of the model vegetation representation was shown to reduce this error substantially for scenes containing a range of vegetation types.

Decomposition in soil and chemical characteristics of pollen

The input to soils made by pollen and its subsequent mineralization has rarely been investigated from a soil microbiological point of view even though the small but significant quantities of C and N in pollen may make an important contribution to nutrient cycling. The relative resistance to decomposition of pollen exines (outer layers) has led to much of the focus of pollen in soil being on its preservation for archaeological and palaeo-ecological purposes. We have examined aspects of the chemical composition and decomposition of pollen from birch (Betula alba) and maize (Zea mays) in soil. The relatively large N contents, small C-to-N ratios and large water-soluble contents of pollen from both species indicated that they would be readily mineralized in soil. When added to soil and incubated at 16 degrees C an amount of C equivalent to 22-26% of the added pollen C was lost as CO2 within 22 days, with the Z. mays pollen decomposing faster. For B. alba pollen, the water-soluble fraction decomposed faster than the whole pollen and the insoluble fraction decomposed more slowly over 22 days. By contrast, there were no significant differences in the decomposition rates of the different fractions from Z. mays pollen. Solid-state C-13 nuclear magnetic resonance (NMR) revealed no gross chemical differences between the pollen of these two species, with strong resonances in the alkyl- and methyl-C region (0-45 p.p.m.) indicative of aliphatic compounds, the O-alkyl-C (60-90 p.p.m.) and the acetal- and ketal-C region (90-110 p.p.m.) indicative of polysaccharides, and the carbonyl-C region indicative of peptides and carboxylic acids. In addition, both pollens gave a small but distinct resonance at 55 p.p.m. attributed to N-alkyl-C. The resonances attributed to polysaccharides were lost completely or substantially reduced after decomposition.

Assessment of soil particle redistribution on two contrasting cultivated hillslopes

The particle size distributions of surface soils from two cultivated silty fields (Moorfield and Railway South) in Herefordshire, UK, were assessed by sampling on 20-m grids across the fields. Moorfield (8 ha) had a uniform landscape sloping mainly in a North-South direction while Railway South (12 ha) had complex undulating landscape characteristics. Samples from 3 surficial layers were also taken at 3 landscape positions at Moorfield to investigate recent (within-season) soil particle redistribution. Size fractions were determined using chemical dispersion, wet sieving (to separate the sand fractions) and laser gramilometry (for the finer fractions). The distribution of various fractions and the relationships between elevation and the various fractions suggest preferential detachment and movement of coarse to very coarse silt fractions (16-63 mu m), which were found mostly at downslope or depositional areas. Upper slope samples had higher clay to fine silt (< 16 mu m) contents than bottom slope samples. The upslope-downslope patterns of size fractions, particularly on uniformly sloping areas, of the 2 fields were similar and their deposited sediments were dominated by coarse silt fractions. Samples from 3 landscape positions at Moorfield became coarser from the less eroded summit, through the eroding side-slope to the bottom-slope depositional area. Within each of these landscape positions the top 0-2.5 cm layers were more enriched in coarse silt fractions than the bottom layers. The spatial patterns of soil particle size distributions in the 2 fields may be a result of sediment detachment and deposition caused by water erosion and tillage operations. (c) 2005 Elsevier B.V. All rights reserved.

Preparation of consistent soil data sets for modelling purposes: Secondary SOTER data for four case study areas

The common GIS-based approach to regional analyses of soil organic carbon (SOC) stocks and changes is to define geographic layers for which unique sets of driving variables are derived, which include land use, climate, and soils. These GIS layers, with their associated attribute data, can then be fed into a range of empirical and dynamic models. Common methodologies for collating and formatting regional data sets on land use, climate, and soils were adopted for the project Assessment of Soil Organic Carbon Stocks and Changes at National Scale (GEFSOC). This permitted the development of a uniform protocol for handling the various input for the dynamic GEFSOC Modelling System. Consistent soil data sets for Amazon-Brazil, the Indo-Gangetic Plains (IGP) of India, Jordan and Kenya, the case study areas considered in the GEFSOC project, were prepared using methodologies developed for the World Soils and Terrain Database (SOTER). The approach involved three main stages: (1) compiling new soil geographic and attribute data in SOTER format; (2) using expert estimates and common sense to fill selected gaps in the measured or primary data; (3) using a scheme of taxonomy-based pedotransfer rules and expert-rules to derive soil parameter estimates for similar soil units with missing soil analytical data. The most appropriate approach varied from country to country, depending largely on the overall accessibility and quality of the primary soil data available in the case study areas. The secondary SOTER data sets discussed here are appropriate for a wide range of environmental applications at national scale. These include agro-ecological zoning, land evaluation, modelling of soil C stocks and changes, and studies of soil vulnerability to pollution. Estimates of national-scale stocks of SOC, calculated using SOTER methods, are presented as a first example of database application. Independent estimates of SOC stocks are needed to evaluate the outcome of the GEFSOC Modelling System for current conditions of land use and climate. (C) 2007 Elsevier B.V. All rights reserved.

Decomposition in soil and chemical characteristics of pollen

The input to soils made by pollen and its subsequent mineralization has rarely been investigated from a soil microbiological point of view even though the small but significant quantities of C and N in pollen may make an important contribution to nutrient cycling. The relative resistance to decomposition of pollen exines (outer layers) has led to much of the focus of pollen in soil being on its preservation for archaeological and palaeo-ecological purposes. We have examined aspects of the chemical composition and decomposition of pollen from birch (Betula alba) and maize (Zea mays) in soil. The relatively large N contents, small C-to-N ratios and large water-soluble contents of pollen from both species indicated that they would be readily mineralized in soil. When added to soil and incubated at 16 degrees C an amount of C equivalent to 22-26% of the added pollen C was lost as CO2 within 22 days, with the Z. mays pollen decomposing faster. For B. alba pollen, the water-soluble fraction decomposed faster than the whole pollen and the insoluble fraction decomposed more slowly over 22 days. By contrast, there were no significant differences in the decomposition rates of the different fractions from Z. mays pollen. Solid-state C-13 nuclear magnetic resonance (NMR) revealed no gross chemical differences between the pollen of these two species, with strong resonances in the alkyl- and methyl-C region (0-45 p.p.m.) indicative of aliphatic compounds, the O-alkyl-C (60-90 p.p.m.) and the acetal- and ketal-C region (90-110 p.p.m.) indicative of polysaccharides, and the carbonyl-C region indicative of peptides and carboxylic acids. In addition, both pollens gave a small but distinct resonance at 55 p.p.m. attributed to N-alkyl-C. The resonances attributed to polysaccharides were lost completely or substantially reduced after decomposition.

The soil-dwelling earthworm Allolobophora chlorotica modifies its burrowing behaviour in response to carbendazim applications

Carbendazim-amended soil was placed above or below unamended soil. Control tests comprised two layers of unamended soil. Allolobophora chlorotica earthworms were added to either the upper or the unamended soil. After 72 h vertical distributions of earthworms were compared between control and carbendazim-amended experiments. Earthworm distributions in the carbendazim-amended test containers differed significantly from the ‘normal’ distribution observed in the control tests. In the majority of the experiments, earthworms significantly altered their burrowing behaviour to avoid carbendazim. However, when earthworms were added to an upper layer of carbendazim-amended soil they remained in this layer. This non-avoidance is attributed to (1) the earthworms’ inability to sense the lower layer of unamended soil and (2) the toxic effect of carbendazim inhibiting burrowing. Earthworms modified their burrowing behaviour in response to carbendazim in the soil. This may explain anomalous results observed in pesticide field trials when carbendazim is used as a control substance.

A universal agro-hydrological model for water and nitrogen cycles in the soil-crop system SMCR_N: critical update and further validation

Agro-hydrological models have widely been used for optimizing resources use and minimizing environmental consequences in agriculture. SMCRN is a recently developed sophisticated model which simulates crop response to nitrogen fertilizer for a wide range of crops, and the associated leaching of nitrate from arable soils. In this paper, we describe the improvements of this model by replacing the existing approximate hydrological cascade algorithm with a new simple and explicit algorithm for the basic soil water flow equation, which not only enhanced the model performance in hydrological simulation, but also was essential to extend the model application to the situations where the capillary flow is important. As a result, the updated SMCRN model could be used for more accurate study of water dynamics in the soil-crop system. The success of the model update was demonstrated by the simulated results that the updated model consistently out-performed the original model in drainage simulations and in predicting time course soil water content in different layers in the soil-wheat system. Tests of the updated SMCRN model against data from 4 field crop experiments showed that crop nitrogen offtakes and soil mineral nitrogen in the top 90 cm were in a good agreement with the measured values, indicating that the model could make more reliable predictions of nitrogen fate in the crop-soil system, and thus provides a useful platform to assess the impacts of nitrogen fertilizer on crop yield and nitrogen leaching from different production systems. (C) 2010 Elsevier B.V. All rights reserved.