Two-scale computational modelling of water flow in unsaturated soils containing irregular-shaped inclusions

The focus of this paper is two-dimensional computational modelling of water flow in unsaturated soils consisting of weakly conductive disconnected inclusions embedded in a highly conductive connected matrix. When the inclusions are small, a two-scale Richards’ equation-based model has been proposed in the literature taking the form of an equation with effective parameters governing the macroscopic flow coupled with a microscopic equation, defined at each point in the macroscopic domain, governing the flow in the inclusions. This paper is devoted to a number of advances in the numerical implementation of this model. Namely, by treating the micro-scale as a two-dimensional problem, our solution approach based on a control volume finite element method can be applied to irregular inclusion geometries, and, if necessary, modified to account for additional phenomena (e.g. imposing the macroscopic gradient on the micro-scale via a linear approximation of the macroscopic variable along the microscopic boundary). This is achieved with the help of an exponential integrator for advancing the solution in time. This time integration method completely avoids generation of the Jacobian matrix of the system and hence eases the computation when solving the two-scale model in a completely coupled manner. Numerical simulations are presented for a two-dimensional infiltration problem.

Direct shear testing on unsaturated silty soils to investigate the effects of drying and wetting on shear strength parameters at low suction

A modified conventional direct shear device was used to measure unsaturated shear strength of two silty soils at low suction values (0 ~ 50 kPa) that were achieved by following drying and wetting paths of soil water characteristic curves (SWCCs). The results revealed that the internal friction angle of the soils was not significantly affected by either the suction or the drying wetting SWCCs. The apparent cohesion of soil increased with a decreasing rate as suction increased. Shear stress-shear displacement curves obtained from soil specimens subjected to the same net normal stress and different suction values showed a higher initial stiffness and a greater peak stress as suction increased. A soil in wetting exhibited slightly higher peak shear stress and more contractive volume change behavior than that of soil in drying at the same net normal stress and suction.

Rice husk biochar and crop residue amendment in subtropical cropping soils: Effect on biomass production, nitrogen use efficiency and greenhouse gas emissions

We investigated the effect of maize residues and rice husk biochar on biomass production, fertiliser nitrogen recovery (FNR) and nitrous oxide (N2O) emissions for three different subtropical cropping soils. Maize residues at two rates (0 and 10 t ha−1) combined with three rates (0, 15 and 30 t ha-1) of rice husk biochar were added to three soil types in a pot trial with maize plants. Soil N2O emissions were monitored with static chambers for 91 days. Isotopic 15N-labelled urea was applied to the treatments without added crop residues to measure the FNR. Crop residue incorporation significantly reduced N uptake in all treatments but did not affect overall FNR. Rice husk biochar amendment had no effect on plant growth and N uptake but significantly reduced N2O and carbon dioxide (CO2) emissions in two of the three soils. The incorporation of crop residues had a contrasting effect on soil N2O emissions depending on the mineral N status of the soil. The study shows that effects of crop residues depend on soil properties at the time of application. Adding crop residues with a high C/N ratio to soil can immobilise N in the soil profile and hence reduce N uptake and/or total biomass production. Crop residue incorporation can either stimulate or reduce N2O emissions depending on the mineral N content of the soil. Crop residues pyrolysed to biochar can potentially stabilise native soil C (negative priming) and reduce N2O emissions from cropping soils thus providing climate change mitigation potential beyond the biochar C storage in soils. Incorporation of crop residues as an approach to recycle organic materials and reduce synthetic N fertiliser use in agricultural production requires a thorough evaluation, both in terms of biomass production and greenhouse gas emissions.

Efficient simulation of unsaturated flow in heterogeneous soils using an exponential integrator

- Provided a practical variable-stepsize implementation of the exponential Euler method (EEM). - Introduced a new second-order variant of the scheme that enables the local error to be estimated at the cost of a single additional function evaluation. - New EEM implementation outperformed sophisticated implementations of the backward differentiation formulae (BDF) of order 2 and was competitive with BDF of order 5 for moderate to high tolerances.