Spatial distribution of soil heavy metal pollution estimated by different interpolation methods: Accuracy and uncertainty analysis

Mapping the spatial distribution of contaminants in soils is the basis of pollution evaluation and risk control. Interpolation methods are extensively applied in the mapping processes to estimate the heavy metal concentrations at unsampled sites. The performances of interpolation methods (inverse distance weighting, local polynomial, ordinary kriging and radial basis functions) were assessed and compared using the root mean square error for cross validation. The results indicated that all interpolation methods provided a high prediction accuracy of the mean concentration of soil heavy metals. However, the classic method based on percentages of polluted samples, gave a pollution area 23.54-41.92% larger than that estimated by interpolation methods. The difference in contaminated area estimation among the four methods reached 6.14%. According to the interpolation results, the spatial uncertainty of polluted areas was mainly located in three types of region: (a) the local maxima concentration region surrounded by low concentration (clean) sites, (b) the local minima concentration region surrounded with highly polluted samples; and (c) the boundaries of the contaminated areas. (C) 2010 Elsevier Ltd. All rights reserved.

Effects of nitrogen fertilization on soil respiration in temperate grassland in Inner Mongolia, China

Nitrogen addition to soil can play a vital role in influencing the losses of soil carbon by respiration in N-deficient terrestrial ecosystems. The aim of this study was to clarify the effects of different levels of nitrogen fertilization (HN, 200 kg N ha(-1) year(-1); MN, 100 kg N ha(-1) year(-1); LN, 50 kg N ha(-1) year(-1)) on soil respiration compared with non-fertilization (CK, 0 kg N ha(-1) year(-1)), from July 2007 to September 2008, in temperate grassland in Inner Mongolia, China. Results showed that N fertilization did not change the seasonal patterns of soil respiration, which were mainly controlled by soil heat-water conditions. However, N fertilization could change the relationships between soil respiration and soil temperature, and water regimes. Soil respiration dependence on soil moisture was increased by N fertilization, and the soil temperature sensitivity was similar in the treatments of HN, LN, and CK treatments (Q (10) varied within 1.70-1.74) but was slightly reduced in MN treatment (Q (10) = 1.63). N fertilization increased soil CO2 emission in the order MN > HN > LN compared with the CK treatment. The positive effects reached a significant level for HN and MN (P < 0.05) and reached a marginally significant level for LN (P = 0.059 < 0.1) based on the cumulative soil respiration during the 2007 growing season after fertilization (July-September 2007). Furthermore, the differences between the three fertilization treatments and CK reached the very significant level of 0.01 on the basis of the data during the first entire year after fertilization (July 2007-June 2008). The annual total soil respiration was 53, 57, and 24% higher than in the CK plots (465 g m(-2) year(-1)). However, the positive effects did not reach the significant level for any treatment in the 2008 growing season after the second year fertilization (July-September 2008, P > 0.05). The pairwise differences between the three N-level treatments were not significant in either year (P > 0.05).

Drivers of soil net nitrogen mineralization in the temperate grasslands in Inner Mongolia, China

Soil net nitrogen mineralization (NNM) of four grasslands across the elevation and precipitation gradients was studied in situ in the upper 0-10 cm soil layer using the resin-core technique in Xilin River basin, Inner Mongolia, China during the growing season of 2006. The primary objectives were to examine variations of NNM among grassland types and the main influencing factors. These grasslands included Stipa baicalensis (SB), Aneulolepidum Chinense (AC), Stipa grandis (SG), and Stipa krylovii (SK) grassland. The results showed that the seasonal variation patterns of NNM were similar among the four grasslands, the rates of NNM and nitrification were highest from June to August, and lowest in September and October during the growing season. The rates of NNM and nitrification were affected significantly by the incubation time, and they were positively correlated with soil organic carbon content, total soil nitrogen (TN) content, soil temperature, and soil water content, but the rates of NNM and nitrification were negatively correlated with available N, and weakly correlated with soil pH and C:N ratio. The sequences of the daily mean rates of NNM and nitrification in the four grasslands during the growing season were AC > SG > SB > SK, and TN content maybe the main affecting factors which can be attributed to the land use type.

Relationships between precipitation, soil water and groundwater at Chongling catchment with the typical vegetation cover in the Taihang mountainous region, China

Vegetation cover plays an important role in the process of evaporation and infiltration. To explore the relationships between precipitation, soil water and groundwater in Taihang mountainous region, China, precipitation, soil water and water table were observed from 2004 to 2006, and precipitation, soil water and groundwater were sampled in 2004 and 2005 for oxygen-18 and deuterium analysis at Chongling catchment. The soil water was sampled at three sites covered by grass (Carex humilis and Carex lanceolata), acacia and arborvitae respectively. Precipitation is mainly concentrated in rainy seasons and has no significant spatial variance in study area. The stable isotopic compositions are enriched in precipitation and soil water due to the evaporation. The analysis of soil water potential and isotopic profiles shows that evaporation of soil water under arborvitae cover is weaker than under grass and acacia, while soil water evaporation under grass and acacia showed no significant difference. Both delta O-18 profiles and soil water potential dynamics reveal that the soil under acacia allows the most rapid infiltration rate, which may be related to preferential flow. In the process of infiltration after a rainstorm, antecedent water still takes up over 30% of water in the topsoil. The soil water between depths of 0-115 cm under grass has a residence time of about 20 days in the rainy season. Groundwater recharge from precipitation mainly occurs in the rainy season, especially when rainstorms or successive heavy rain events happen.

Probabilistic modelling of soil moisture dynamics of irrigated cropland in the North China Plain

A probabilistic soil moisture dynamic model is used to estimate the soil moisture probability distribution and plant water stress of irrigated cropland in the North China Plain. Soil moisture and meteorological data during the period of 1998 to 2003 were obtained from an irrigated cropland ecosystem with winter wheat and maize in the North China Plain to test the probabilistic soil moisture dynamic model. Results showed that the model was able to capture the soil moisture dynamics and estimate long-term water balance reasonably well when little soil water deficit existed. The prediction of mean plant water stress during winter wheat and maize growing season quantified the suitability of the wheat-maize rotation to the soil and climate environmental conditions in North China Plain under the impact of irrigation. Under the impact of precipitation fluctuations, there is no significant bimodality of the average soil moisture probability density function.

Effects of inoculated Microcoleus vaginatus on the structure and function of biological soil crusts of desert

Microcoleus vaginatus Gom., the dominant species in biological soil crusts (BSCs) in desert regions, plays a significant role in maintaining the BSC structure and function. The BSC quality is commonly assessed by the chlorophyll a content, thickness, and compressive strength. Here, we have studied the effect of different proportions of M. vaginatus, collected from the Gurbantunggut Desert in northwestern China, on the BSC structure and function under laboratory conditions. We found that when M. vaginatus was absent in the BSC, the BSC coverage, quantified by the percentage of BSC area to total land surface area, was low with a chlorophyll a content of 4.77 x 10(-2) mg g(-1) dry soil, a thickness of 0.86 mm, and a compressive strength of 12.21 Pa. By increasing the percentage of M. vaginatus in the BSC, the BSC coverage, chlorophyll a content, crust thickness, and compressive strength all significantly increased (P < 0.01). The maximum chlorophyll a content (13.12 mg g(-1)dry soil), the highest crust thickness, and the compressive strength (1.48 mm and 36.60 Pa, respectively) occurred when the percentage of inoculated M. vaginatus reached 80% with a complex network of filaments under scanning electron microscope. The BSC quality indicated by the above variables, however, declined when the BSC was composed of pure M. vaginatus (monoculture). In addition, we found that secretion of filaments and polymer, which stick sands together in the BSC, increased remarkably with the increase of the dominant species until the percentage of M. vaginatus reached 80%. Our results suggest that not only the dominant species but also the accompanying taxa are critical for maintaining the structure and functions of the BSC and thus the stability of the BSC ecosystems.

Stable Carbon Isotopic Composition in Soil Organic Carbon and C3/C4 Source Variations at the Haibei Alpine Meadow

Stable carbon isotopes of organic matter originated from different soil layers (0～5 cm, 5～15 cm, 15～25 cm, 25～35 cm, 35～50 cm, 50～65 cm) were investigated in the Haibei Alpine Meadow Ecosystem Research Station of the Chinese Academy of Sciences. The preliminary results indicated that δ13C values of soil organic matter increased with increased soil depth. δ13C of soil organic carbon in 0～5 cm layer showed the lowest value, -25.09‰; while 50～65 cm soil layer possessed the lowerδ13C value, -13.87‰. Based on mass balance model of stable isotopes, it was proposed that the percentage of C4 carbon source tend to increase with increased soil depth. The preliminary study indicated that alpine meadow might have undergone a successive process from C4-dominated community to C3-dominated one. However, changing δ13C values in atmospheric CO2 overtime and different processes of soil organic carbon formation (or eluviation) might somewhat contribute to increasing δ13C values. In this case, mass balance model would underestimate C3 community and overestimate C4 community.

Soil organic carbon storage and soil CO2 flux in the alpine meadow ecosystem

High-resolution sampling, measurements of organic carbon contents and C-14 signatures of selected four soil profiles in the Haibei Station situated on the northeast Tibetan Plateau, and application of C-14 tracing technology were conducted in an attempt to investigate the turnover times of soil organic carbon and the soil-CO2 flux in the alpine meadow ecosystem. The results show that the organic carbon stored in the soils varies from 22.12x10(4) kg C hm(-2) to 30.75x10(4) kg C hm(-2) in the alpine meadow ecosystems, with an average of 26.86x10(4) kg C hm(-2). Turnover times of organic carbon pools increase with depth from 45 a to 73 a in the surface soil horizon to hundreds of years or millennia or even longer at the deep soil horizons in the alpine meadow ecosystems. The soil-CO2 flux ranges from 103.24 g C m(-2) a(-1) to 254.93 gC m(-2) a(-1), with an average of 191.23 g C m(-2) a(-1). The CO2 efflux produced from microbial decomposition of organic matter varies from 73.3 g C m(-2) a(-1) to 181 g C m(-2) a(-1). More than 30% of total soil organic carbon resides in the active carbon pool and 72.8%. 81.23% of total CO2 emitted from organic matter decomposition results from the topsoil horizon (from 0 cm to 10 cm) for the Kobresia meadow. Responding to global warming, the storage, volume of flow and fate of the soil organic carbon in the alpine meadow ecosystem of the Tibetan Plateau will be changed, which needs further research.

The thermo-visco-plastic instability analysis of saturated soil

A theoretical analysis of instability of saturated soil is presented considering the simple shearing of a heat conducting thermo-visco-plastic material. It is shown that the instability is mainly the consequence of thermal softening which overcomes the strain hardening and the other type of instability is controlled by strain softening. The effects of other factors such as permeability to the instability are discussed in this paper.

On the shear instability of saturated soil

In this paper. the dynamic instability of simple shear of saturated soil is discussed. The governing equations are obtained based on mixture theory in which the inertia effect and the compressibility of grains are considered. Perturbation method is used to analyze and it is shown that two types of instability may exist. One of them is dominated by pore-pressure-softening, while the other by strain-softening.

An Analytical Solution for Three-Dimensional Diffraction of Plane P-Waves by a Hemispherical Alluvial Valley with Saturated Soil Deposits

An analytical solution for the three-dimensional scattering and diffraction of plane P-waves by a hemispherical alluvial valley with saturated soil deposits is developed by employing Fourier-Bessel series expansion technique. Unlike previous studies, in which the saturated soil deposits were simulated with the single-phase elastic theory, in this paper, they are simulated with Biot's dynamic theory for saturated porous media, and the half space is assumed as a single-phase elastic medium. The effects of the dimensionless frequency, the incidence angle of P-wave and the porosity of soil deposits on the surface displacement magnifications of the hemispherical alluvial valley are investigated. Numerical results show that the existence of a saturated hemispherical alluvial valley has much influence on the surface displacement magnifications. It is more reasonable to simulate soil deposits with Biot's dynamic theory when evaluating the displacement responses of a hemispherical alluvial valley with an incidence of P-waves.

Comparison of the Dynamic Response of a Saturated Soil Between Two Models

By comparing the dynamic responses of saturated soil to Biot's and Yamamoto's models, the properties of the two models have be pointed out. First of all, an analysis has been made for energy loss of each model from the basic equations. Then the damping of elastic waves in coarse sand and fine sand with loading frequency and soil's parameters have been calculated and the representation of viscous friction and Coulomb friction in the two models has been concluded. Finally, the variations of loading wave damping and stress phase angles with water depth and soil's parameters have been obtained as loading waves range in ocean waves.

Numerical Study on Rainfall Infiltration in Rock-Soil Blocks

A mathematical model for the rain infiltration in the rock-soil slop has been established and solved by using the finite element method. The unsteady water infiltrating process has been simulated to get water content both in the homogeneous and heterogeneous media. The simulated results show that the rock blocks in the rock-soil slop can cause the wetting front moving fast. If the rain intensity is increased, the saturated region will be formed quickly while other conditions are the same. If the rain intensity keeps a constant, it is possible to accelerate the generation of the saturated region by properly increasing the vertical filtration rate of the rock-soil slop. However, if the vertical filtration rate is so far greater than the rain intensity, it will be difficult to form the saturated region in the rock-soil slop. The numerical method was verified by comparing the calculation results with the field test data.