The effects of groundwater table and flood irrigation strategies on soil water and salt dynamics and reed water use in the Yellow River Delta, China

Vegetation management in shallow groundwater table environments requires an understanding of the interactions between the physical and biological factors that determine root-zone soil salinization and moisture. In this study, the effects of groundwater depth and flood irrigation strategies on water and salt dynamics and reed water use were analyzed in the shallow groundwater region of the Yellow River Delta in China using the HYDRUS-1D model. The results indicated that there is a conflict between water, salt stress, and reed water use with variations in groundwater depth. A water table depth of 3.5 m is the minimum limit to maintain a safe level of soil salinity, but at this depth, the environmental stress on reeds is worsened by the decrease in soil water storage. Maintaining the flood pulses on the wetland, especially during May, may be critical for restoring the reed wetland in the Yellow River Delta.     

Dynamics of soil organic matter in primary and secondary forest succession on sandy soils in The Netherlands: An application of the ROMUL model

We applied the simulation model ROMUL of soil organic matter dynamics in order to analyse and predict forest soil organic matter (SOM) changes following stand growth and also to identify gaps of data and modelling problems. SOM build-up was analysed (a) from bare sand to forest soil during a primary succession in Scots pine forest and (b) on mature forest soil under Douglas fir plantations as an example of secondary succession in The Netherlands. As some of the experimental data were unreliable we compiled a set of various scenarios with different soil moisture regime, initial SOM pools and amount and quality of above and below ground litter input. This allowed us to find the scenarios that reflect the SOM dynamics more realistically. In the Scots pine forest, total litter input was estimated as 0.50 kg m⁻² year⁻¹. Two scenarios were defined for the test runs: (a) forest floor moisture regimes—‘dry, mesic and hydric’ and (b) augmenting a root litter pool with three ratios of needles and branches to roots: 1:1, 1:1.5 and 1:2.0. The scenario finally compiled had the following characteristics: (a) climate for dry site with summer drought and high winter moisture of forest floor; (b) a litter input of 0.25 kg m⁻² year⁻¹ above ground and 0.50 kg m⁻² year⁻¹ below ground; (c) a low nitrogen and ash content in all litter fall fractions. The test runs for the estimation of the initial SOM pools and the amount and proportion of above and below ground litter fall were also performed in the Douglas fir plantation. The inputs of above ground litter tested in various combinations were 0.30 and 0.60 kg m⁻² year⁻¹, and below ground litter 0.30, 0.60 and 0.90 kg m⁻² year⁻¹. The scenario that fitted the experimental data had an SOM pool of 20–25 kg m⁻², an aboveground litter input of 0.6 kg m⁻² year⁻¹and a below ground litter input of 0.9 kg m⁻² year⁻¹. The long-term simulation corresponded well with the observed patterns of soil organic matter accumulation associated with the forest soil development in primary and secondary succession. During primary succession in Scots pine forest on dry sand there is a consistent accumulation of a raw humus forest floor. The soil dynamics in the Douglas fir plantation also coincide with the observed patterns of SOM changes during the secondary succession, with SOM decreasing significantly under young forest, and SOM being restored in the older stands.     

Integrating soil carbon cycling with that of nitrogen and phosphorus in the watershed model SWAT: Theory and model testing

In this paper we describe and test a sub-model that integrates the cycling of carbon (C), nitrogen (N) and phosphorus (P) in the Soil Water Assessment Tool (SWAT) watershed model. The core of the sub-model is a multi-layer, one-pool soil organic carbon (S_C) algorithm, in which the decomposition rate of S_C and input rate to S_C (through decomposition and humification of residues) depend on the current size of S_C. The organic N and P fluxes are coupled to that of C and depend on the available mineral N and P, and the C:N and N:P ratios of the decomposing pools. Tillage explicitly affects the soil organic matter turnover rate through tool-specific coefficients. Unlike most models, the turnover of soil organic matter does not follow first order kinetics. Each soil layer has a specific maximum capacity to accumulate C or C saturation (S_x) that depends on texture and controls the turnover rate. It is shown in an analytical solution that S_x is a parameter with major influence in the model C dynamics. Testing with a 65-yr data set from the dryland wheat growing region in Oregon shows that the model adequately simulates the S_C dynamics in the topsoil (top 0.3 m) for three different treatments. Three key model parameters, the optimal decomposition and humification rates and a factor controlling the effect of soil moisture and temperature on the decomposition rate, showed low uncertainty as determined by generalized likelihood uncertainty estimation. Nonetheless, the parameter set that provided accurate simulations in the topsoil tended to overestimate S_C in the subsoil, suggesting that a mechanism that expresses at depth might not be represented in the current sub-model structure. The explicit integration of C, N, and P fluxes allows for a more cohesive simulation of nutrient cycling in the SWAT model. The sub-model has to be tested in forestland and rangeland in addition to agricultural land, and in diverse soils with extreme properties such high or low pH, an organic horizon, or volcanic soils.