A universal equation for predicting rainfall-erosion losses : an aid to conservation farming in humid regions.

Cover title.

A mass-conservative control volume-finite element method for solving Richards’ equation in heterogeneous porous media

We present a mass-conservative vertex-centred finite volume method for efficiently solving the mixed form of Richards’ equation in heterogeneous porous media. The spatial discretisation is particularly well-suited to heterogeneous media because it produces consistent flux approximations at quadrature points where material properties are continuous. Combined with the method of lines, the spatial discretisation gives a set of differential algebraic equations amenable to solution using higher-order implicit solvers. We investigate the solution of the mixed form using a Jacobian-free inexact Newton solver, which requires the solution of an extra variable for each node in the mesh compared to the pressure-head form. By exploiting the structure of the Jacobian for the mixed form, the size of the preconditioner is reduced to that for the pressure-head form, and there is minimal computational overhead for solving the mixed form. The proposed formulation is tested on two challenging test problems. The solutions from the new formulation offer conservation of mass at least one order of magnitude more accurate than a pressure head formulation, and the higher-order temporal integration significantly improves both the mass balance and computational efficiency of the solution.

3-D kinematical conservation laws (KCL): Evolution of a surface in R-3 - in particular propagation of a nonlinear wavefront

3-D KCL are equations of evolution of a propagating surface (or a wavefront) Omega(t), in 3-space dimensions and were first derived by Giles, Prasad and Ravindran in 1995 assuming the motion of the surface to be isotropic. Here we discuss various properties of these 3-D KCL.These are the most general equations in conservation form, governing the evolution of Omega(t) with singularities which we call kinks and which are curves across which the normal n to Omega(t) and amplitude won Omega(t) are discontinuous. From KCL we derive a system of six differential equations and show that the KCL system is equivalent to the ray equations of 2, The six independent equations and an energy transport equation (for small amplitude waves in a polytropic gas) involving an amplitude w (which is related to the normal velocity m of Omega(t)) form a completely determined system of seven equations. We have determined eigenvalues of the system by a very novel method and find that the system has two distinct nonzero eigenvalues and five zero eigenvalues and the dimension of the eigenspace associated with the multiple eigenvalue 0 is only 4. For an appropriately defined m, the two nonzero eigenvalues are real when m > 1 and pure imaginary when m < 1. Finally we give some examples of evolution of weakly nonlinear wavefronts.

An Application of 3-D Kinematical Conservation Laws: Propagation of a 3-D Wavefront

Three-dimensional (3-D) kinematical conservation laws (KCL) are equations of evolution of a propagating surface Omega(t) in three space dimensions. We start with a brief review of the 3-D KCL system and mention some of its properties relevant to this paper. The 3-D KCL, a system of six conservation laws, is an underdetermined system to which we add an energy transport equation for a small amplitude 3-D nonlinear wavefront propagating in a polytropic gas in a uniform state and at rest. We call the enlarged system of 3-D KCL with the energy transport equation equations of weakly nonlinear ray theory (WNLRT). We highlight some interesting properties of the eigenstructure of the equations of WNLRT, but the main aim of this paper is to test the numerical efficacy of this system of seven conservation laws. We take several initial shapes for a nonlinear wavefront with a suitable amplitude distribution on it and let it evolve according to the 3-D WNLRT. The 3-D WNLRT is a weakly hyperbolic 7 x 7 system that is highly nonlinear. Here we use the staggered Lax-Friedrichs and Nessyahu-Tadmor central schemes and have obtained some very interesting shapes of the wavefronts. We find the 3-D KCL to be suitable for solving many complex problems for which there presently seems to be no other method capable of giving such physically realistic features.

Eigenvalues of kinematical conservation laws (KCL) based 3-D weakly nonlinear ray theory (WNLRT)

System of kinematical conservation laws (KCL) govern evolution of a curve in a plane or a surface in space, even if the curve or the surface has singularities on it. In our recent publication K. R. Arun, P. Prasad, 3-D kinematical conservation laws (KCL): evolution of a surface in R-3-in particular propagation of a nonlinear wavefront, Wave Motion 46 (2009) 293-311] we have developed a mathematical theory to study the successive positions and geometry of a 3-D weakly nonlinear wavefront by adding an energy transport equation to KCL. The 7 x 7 system of equations of this KCL based 3-D weakly nonlinear ray theory (WNLRT) is quite complex and explicit expressions for its two nonzero eigenvalues could not be obtained before. In this short note, we use two different methods: (i) the equivalence of KCL and ray equations and (ii) the transformation of surface coordinates, to derive the same exact expressions for these eigenvalues. The explicit expressions for nonzero eigenvalues are important also for checking stability of any numerical scheme to solve 3-D WNLRT. (C) 2010 Elsevier Inc. All rights reserved.

Nonlinear Dynamics for Local Fractional Burgers Equation Arising in Fractal Flow

The local fractional Burgers’ equation (LFBE) is investigated from the point of view of local fractional conservation laws envisaging a nonlinear local fractional transport equation with a linear non-differentiable diffusion term. The local fractional derivative transformations and the LFBE conversion to a linear local fractional diffusion equation are analyzed.

Wave-activity conservation laws and stability theorems for semi-geostrophic dynamics: part 1. pseudomomentum-based theory

There exists a well-developed body of theory based on quasi-geostrophic (QG) dynamics that is central to our present understanding of large-scale atmospheric and oceanic dynamics. An important question is the extent to which this body of theory may generalize to more accurate dynamical models. As a first step in this process, we here generalize a set of theoretical results, concerning the evolution of disturbances to prescribed basic states, to semi-geostrophic (SG) dynamics. SG dynamics, like QG dynamics, is a Hamiltonian balanced model whose evolution is described by the material conservation of potential vorticity, together with an invertibility principle relating the potential vorticity to the advecting fields. SG dynamics has features that make it a good prototype for balanced models that are more accurate than QG dynamics. In the first part of this two-part study, we derive a pseudomomentum invariant for the SG equations, and use it to obtain: (i) linear and nonlinear generalized Charney–Stern theorems for disturbances to parallel flows; (ii) a finite-amplitude local conservation law for the invariant, obeying the group-velocity property in the WKB limit; and (iii) a wave-mean-flow interaction theorem consisting of generalized Eliassen–Palm flux diagnostics, an elliptic equation for the stream-function tendency, and a non-acceleration theorem. All these results are analogous to their QG forms. The pseudomomentum invariant – a conserved second-order disturbance quantity that is associated with zonal symmetry – is constructed using a variational principle in a similar manner to the QG calculations. Such an approach is possible when the equations of motion under the geostrophic momentum approximation are transformed to isentropic and geostrophic coordinates, in which the ageostrophic advection terms are no longer explicit. Symmetry-related wave-activity invariants such as the pseudomomentum then arise naturally from the Hamiltonian structure of the SG equations. We avoid use of the so-called ‘massless layer’ approach to the modelling of isentropic gradients at the lower boundary, preferring instead to incorporate explicitly those boundary contributions into the wave-activity and stability results. This makes the analogy with QG dynamics most transparent. This paper treats the f-plane Boussinesq form of SG dynamics, and its recent extension to β-plane, compressible flow by Magnusdottir & Schubert. In the limit of small Rossby number, the results reduce to their respective QG forms. Novel features particular to SG dynamics include apparently unnoticed lateral boundary stability criteria in (i), and the necessity of including additional zonal-mean eddy correlation terms besides the zonal-mean potential vorticity fluxes in the wave-mean-flow balance in (iii). In the companion paper, wave-activity conservation laws and stability theorems based on the SG form of the pseudoenergy are presented.

Observational and energetics constraints on the non-conservation of potential/Conservative Temperature and implications for ocean modelling

This paper seeks to elucidate the fundamental differences between the nonconservation of potential temperature and that of Conservative Temperature, in order to better understand the relative merits of each quantity for use as the heat variable in numerical ocean models. The main result is that potential temperature is found to behave similarly to entropy, in the sense that its nonconservation primarily reflects production/destruction by surface heat and freshwater fluxes; in contrast, the nonconservation of Conservative Temperature is found to reflect primarily the overall compressible work of expansion/contraction. This paper then shows how this can be exploited to constrain the nonconservation of potential temperature and entropy from observed surface heat fluxes, and the nonconservation of Conservative Temperature from published estimates of the mechanical energy budgets of ocean numerical models. Finally, the paper shows how to modify the evolution equation for potential temperature so that it is exactly equivalent to using an exactly conservative evolution equation for Conservative Temperature, as was recently recommended by IOC et al. (2010). This result should in principle allow ocean modellers to test the equivalence between the two formulations, and to indirectly investigate to what extent the budget of derived nonconservative quantities such as buoyancy and entropy can be expected to be accurately represented in ocean models.

ENO adaptive method for solving one-dimensional conservation laws

In this work an efficient third order non-linear finite difference scheme for solving adaptively hyperbolic systems of one-dimensional conservation laws is developed. The method is based oil applying to the solution of the differential equation an interpolating wavelet transform at each time step, generating a multilevel representation for the solution, which is thresholded and a sparse point representation is generated. The numerical fluxes obtained by a Lax-Friedrichs flux splitting are evaluated oil the sparse grid by an essentially non-oscillatory (ENO) approximation, which chooses the locally smoothest stencil among all the possibilities for each point of the sparse grid. The time evolution of the differential operator is done on this sparse representation by a total variation diminishing (TVD) Runge-Kutta method. Four classical examples of initial value problems for the Euler equations of gas dynamics are accurately solved and their sparse solutions are analyzed with respect to the threshold parameters, confirming the efficiency of the wavelet transform as an adaptive grid generation technique. (C) 2008 IMACS. Published by Elsevier B.V. All rights reserved.

Bicomplexes and conservation laws in non-Abelian Toda models

A bicomplex structure is associated with the Leznov-Saveliev equation of integrable models. The linear problem associated with the zero-curvature condition is derived in terms of the bicomplex linear equation. The explicit example of a non-Abelian conformal affine Toda model is discussed in detail and its conservation laws are derived from the zero-curvature representation of its equation of motion.

Conservation laws for non-potential operators

A study was developed in order to build a function M invariant in time, by means of Hamiltonian's formulation, taking into account the equation associated to the problem, showing that starting from this function the equation of motion of the system with the contour conditions for non-conservative considered problems can be obtained. The Hamiltonian method is extended for these kind of systems in order to validate for non-potential operators through variational approach.

The first Brazilian municipal initiative of payments for environmental services and its potential for soil conservation

It has become evident that policies aimed at mitigating the growing water resources and water use conflicts in Brazil are crucial. The municipality of Extrema in Minas Gerais state in Brazil pioneered the first Brazilian municipal PES initiative (Conservador das Aguas program), based on the relationship between forests and the benefits they provide. This study aimed to assess soil loss in the Posses sub-basin, where the Conservador das Aguas program began. Additionally, we aimed to determine the potential that this PES initiative has for soil conservation, as well as to minimize the soil losses as a function of forest area size and location in order to propose a technical approach for implementing PES. In this sense, considering the prescribed conservation practices, land use situation, and soil cover in the Posses sub-basin, we analyzed the effectiveness of the Conservador das Aguas program before and after implementation in relation to reduced soil loss under 36 different land use and soil cover scenarios. We used a geographic information system (GIS) for spatializing and producing different information plans and the Revised Universal Soil Loss Equation (RUSLE) for estimating soil loss. As a result, we found that minimization of soil loss may be obtained by adopting pasture conservation practices. Additionally the expected average soil loss in the Posses sub-basin under conditions of land use and soil cover, before and after implementing the water conservation program was 30.63 and 7.06 Mg ha(-1) year(-1), respectively. (C) 2014 Elsevier B.V. All rights reserved.

Wide localized solutions of the parity-time-symmetric nonautonomous nonlinear Schrodinger equation

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Geographic information system-based decision support for soil conservation planning in Tajikistan

Soil erosion on sloping agricultural land poses a serious problem for the environment, as well as for production. In areas with highly erodible soils, such as those in loess zones, application of soil and water conservation measures is crucial to sustain agricultural yields and to prevent or reduce land degradation. The present study, carried out in Faizabad, Tajikistan, was designed to evaluate the potential of local conservation measures on cropland using a spatial modelling approach to provide decision-making support for the planning of spatially explicit sustainable land use. A sampling design to support comparative analysis between well-conserved units and other field units was established in order to estimate factors that determine water erosion, according to the Revised Universal Soil Loss Equation (RUSLE). Such factor-based approaches allow ready application using a geographic information system (GIS) and facilitate straightforward scenario modelling in areas with limited data resources. The study showed first that assessment of erosion and conservation in an area with inhomogeneous vegetation cover requires the integration of plot-based cover. Plot-based vegetation cover can be effectively derived from high-resolution satellite imagery, providing a useful basis for plot-wise conservation planning. Furthermore, thorough field assessments showed that 25.7% of current total cropland is covered by conservation measures (terracing, agroforestry and perennial herbaceous fodder). Assessment of the effectiveness of these local measures, combined with the RUSLE calculations, revealed that current average soil loss could be reduced through low-cost measures such as contouring (by 11%), fodder plants (by 16%), and drainage ditches (by 53%). More expensive measures such as terracing and agroforestry can reduce erosion by as much as 63% (for agroforestry) and 93% (for agroforestry combined with terracing). Indeed, scenario runs for different levels of tolerable erosion rates showed that more cost-intensive and technologically advanced measures would lead to greater reduction of soil loss. However, given economic conditions in Tajikistan, it seems advisable to support the spread of low-cost and labourextensive measures.