Biostatistical Models of Antibiotic Residues in Soil Water Due to the Use of Poultry Manure

The large-scale and intensified development of livestock and poultry breeding has resulted in a large number of veterinary antibiotics being used in the aquaculture industry and increasing year by year. Accumulation of antibiotics in the organism.

To satisfy the increasing demand for vegetables, despite the poverty of coastal soils and land pressure, farmers tend to intensify production by using mineral and organic fertilizers and pesticides.

Today, poultry manure is frequently used as fertilizer in the study area. And poultry manure has been used as effective fertilizers for centuries. It is very low, after the body absorbs a small part, it undergoes metabolic reactions such as hydroxylation, cleavage, and glucuronidation to produce inactive products, and about 60% to 90% pass through feces as it is Excretion of urine and excrement of livestock and poultry have undoubtedly become one of the main sources of environmental pollution by antibiotics.

The drug design of antibiotics is mainly to kill pathogenic pathogens. Once they enter the environment, they will inevitably pose a potential threat to other organisms in the environment. Firstly, the resistant pathogenic bacteria or variant pathogens in livestock and poultry are produced and continuously released into the environment. Secondly, livestock and poultry continue to excrete these antibiotics or their metabolites into the environment, so that drug resistant pathogens and mutant pathogens in the environment are continuously generated. Both of these in turn stimulated the producers to increase the dosage and renew the drug varieties, which resulted in a vicious circle of “drug-contaminated environment → emergence of drug-resistant or mutated pathogens → increased dosages → environmental pollution”.

Fertilization by poultry manure has shown an important variance in soil water chemical characteristics. The state for regional soil moisture reserve is the strategic storage of water resources in the district.

The distribution of soil moisture directly affects the supply of groundwater resources, determines the amount of water, which is absorbed from the soil and evaporated by the earth’s surface plants, plays a decisive factor for plants’ productivity, and also is regarded as strategic factors influencing the ecological environment security, the economic development and the people’s lives in arid and semiarid areas [1].

The intensive use of poultry manure and other animal feces are a significant environmental risk to soil water. Presence of antibiotic and soil water physicochemical properties played key roles in degradation of numerous molecules and other processing. Fertilization is the commonest managing agricultural soils, and for a long time, intensive farming appealed to fertilizer to increase yields.

Then livestock and poultry excrement will constitute one of the main reasons for the non-point source pollution of antibiotics in China. Therefore, based on the research results in recent years, combining the use of antibiotics in the livestock and poultry industry and the status of residues, this paper analyzes the fate and environmental risks of livestock and poultry excrement in soil, and puts forward the corresponding mathematical model.

Darcy law is a fundamental theoretical method to describe the motion law of soil moisture, therefore a variety of Richards equations are deduced. For the nonlinear partial differential equations, the previous research method is to discuss their definite solutions and we commonly can acquire their numerical solution through the numerical method.

Whether the analytical solution of Richards equations, which describes the change of soil moisture content with the change of time and space position in Darcy’s law, has been the expectation. If we substitute some empirical representations of hydraulic conductivity and water diffusivity into Richards equations, the exact solution of Richards equation on soil moisture content, soil depth and time is of great significance. Furthermore, the greater parts of antibiotic were found in the soil water by poultry manure.

According to the literature, few studies have dealt with the dynamics of antibiotic residues in soil water due to the use of poultry manure. In the past few years, many powerful methods to construct exact solutions of nonlinear evolution equations have been established and developed such as the homogeneous balance method [2-4], the (G′/G) -expansion method [5,6], the exp-function method [7,8] and so on. One of the most effective and direct methods for constructing exact solutions of nonlinear differential equations is the (G′/G) - expansion method.

The (G′/G) - expansion method, first introduced by Wang et al. [5], has been widely used to search for various exact solutions of NLEEs [9-11]. The (G′/G) - expansion method is based on the explicit linearization of nonlinear differential equations for traveling waves with a certain substitution which leads to a second-order differential equation with constant coefficients [14-16]. Finding an exact solution for Richards equation, by using the (G′ /G) - expansion method, is the main goal of the present study.

First of all, we introduce a form of Richards equations as follows:

where D(θ) denotes water diffusivity; K(θ) denotes hydraulic conductivity; t denotes time; θ denotes soil moisture content; x, y, z denote coordinate axes.

If the soil moisture content is lower than the saturated (unsaturated) moisture content with little change, we take as D(θ) = a , where a is a constant.

Many researchers have committed themselves to estimating soil hydraulic conductivity, as a result, various empirical representations of hydraulic conductivity are proposed. We assume that [12,13] unsaturated hydraulic conductivity is calculated by using the Libardi method, that is

where β is a constant; K₀ and θ₀ are the values of K and θ during steady-state infiltration, respectively.

Next, we have intend to simplify equation (2.1), in other words, here we only consider the case that soil moisture flows in the vertical direction, and therefore we have

By substituting D(θ) = a and equation (2.2) into equation (2.3), hence the following semi-empirical Richards equation is obtained:

In this section, by make use of the (G′/G) -expansion method, we obtain an exact solution for the equation (2.4), however, we omit the description of the (G′/G) -expansion method. If you are interested in this method, you can refer to the reference [4].

Let then the equation (2.4) can be equivalently changed into

Using the travelling wave variable θ(z,t) =θ(ξ) and ξ =z − ωt carries out the equation (2.5) into an ordinary differential equation for θ =θ (ξ)

In order to apply the (G′/G) -expansion method, we use the Painlevé transformation v = e^βθ, or equivalently hence the equation (2.6) can be written as

Suppose that the solution of ordinary differential equation (2.7) can be expressed by a polynomial in (G′/G) as follows:

where G = G(ξ) satisfies the second order LODE in the form

where are real constants to be determined later.

According to the (G′/G) - expansion method, considering the homogeneous balance between vv" and v'v² in the equation (2.7), we get 3n +1 = 2n + 2⇒n =1, hence we can write (2.8) as

Substituting (2.10) along with (2.9) into (2.7) and collecting all terms with the same order of together, the left-hand side of (2.7) are converted into polynomial in . Setting each coefficient of each polynomial to zero, we derive a set of algebraic equations for as follows:

Solving the algebraic equations above yields

where a₀ is arbitrary constant.

By using (2.11), (2.10) can be written as

where ξ = z.

Substituting the general solutions of equation (2.9) into (2.12), we have an exact solution of the equation (2.7) as follows:

Where and a₀ are arbitrary constants.

Therefore by we have an exact solution of the equation (2.4) as follows:

According to what has been discussed above, we assume that soil moisture content, soil depth and time satisfy

where are arbitrary constants.

Data used here were collected during an internal drainage experiment, carried out on a sandy–loam Red Yellow Latosol (Typic Hapludox) of the county of Piracicaba, SP, Brazil, with a fairly homogeneous profile down to the depth of 2m [17]. They consist of soil water content values calculated from tensiometer readings, through the use of laboratory established soil water retention curves and of soil water potential heads, measured with the same mercury manometer tensiometers. The soil water retention curves were established by drying wet samples, to minimise hysteresis effects. Data are averages of 3 replicates collected during the internal drainage of three 3m by 3m plots, 10m apart, previously submitted to infiltration by ponding water on the soil surface until quasi-steady flow.

To prevent evaporation losses, the soil surface was covered with a plastic sheet during the drainage process, so that all conditions stated in Libardi et al. [18] were met. Soil-water retention curve at each soil depth was determined using undisturbed samples (45mm diameter by 45mm height) on porous plate funnels (matric potential head h from 0.26m to 1.87m) and porous plate pressure cells (matric potential head h from 3.0m to 8.0m). Soil bulk densities were also determined using volumetric cylinders of 45 mm diameter and 45mm height. Table 2 also presents the obtained values of parameters a, m, n, r, and the equation of Van Genuchten [19] as well as their r2 values; soil water content values during the redistribution time were obtained by means of these equations.

The saturated soil water contents 0 were also taken from the Van Genuchenten equation for matric potential h = 0. The determination of 0 in the field at steady infiltration is practically impossible, except with neutron probes, which have to be well calibrated and present difficulties for measurements close to the soil surface. It is important to note that it is difficult to fully saturate a soil profile in the field and that, therefore, field values of 0 are, in general, slightly lower than those measured in the laboratory.

(Tables 2 & 3) present the average raw data of soil water content and matric potential head, which are presented here for any other development of soil water movement theories. Soil water potential head H (cm of water) is considered as the sum of the gravitational head z and the matric potential head h [17].

In this model, the data were collected during an internal drainage experiment, which was carried out on a sandy–loam Red Yellow Latosol (Typic Hapludox) of the county of Piracicaba [17]. We use the data to determine the parameters in the identity (2.15), and therefore the predicted values of soil moisture content are calculated from (2.15) (Table 1 & 2).

From data in (Table 2), it revealed that volumetric soil water contents θ (m³/m³), for different redistribution times (t) and different depths (z).And the parameters of equation (2.15) for each time and soil depth were determined by using data fitting method. All of these can be seen in (Table 1). These parameters are shown in (Table 3). From these parameters, we obtain a function express in which soil water content is independent variable whereas time and soil depth are two dependent variables. And thus the predicted values of soil moisture content are also presented in (Table 4) according to equation (2.15). Compared with the actual values, the mean squared error (MSE) is 0.00077423, which suggests that equation (2.15) describes the changes of soil water content with the time and soil depth very well. According to (Table 2), the following (Figure 1) can be presented.

Table 1: The Soil Profile Characterisation.

Table 2: Volumetric Soil Water Contents θ (M³/M³), For Different Redistribution Times (T) And Different Depths (Z).

Table 3: Parameters of equation (2.15).

Table 4: The Predicted Values.

The (Figure 1) indicates that the oil water contents θ is decreasing with the increase of times, while is almost the same contents at the range of 0-14m.

(Figure 2) shows that equation (2.15) has the same tendency with the (Figure 1). Although the contents are not equal at the range of 0-14m, it can be recorded as the same contents duo to little change in this range. Since the hypothesis that the water diffusivity is constant, figure 2 obviously conforms to the actual situation.

We propose a new model equation (2.15), which describes the changes of soil water content with the time and soil depth, by finding out an exact solution of equation (2.4) through currently prevalent (G′/G) -expansion method. Even though equation (2.15) is established with the assumptions that the hydraulic conductivity satisfies exponential function and water diffusivity is a constant, to some extent, the very low mean squared error (MSE) indicates that equation (2.15) is a great reflection in change of soil water content with the time and soil depth. However, the simple assumption of water diffusivity might incur larger error in term of actual values of soil water content, thus a precise empirical representation of water diffusivity is expected to be proposed. With the rapid development of science and technology, the solutions to nonlinear differential equations would be enriched, and more excellent results on Darcy’s law might be acquired. All these results also suggest that soil water play a considerable role in the fate of antibiotic in the environment. Further research is needed to give the detection of antibiotic residues from poultry manure in soil water by the regression equations constructed in the present study.

We have no conflict of interests to disclose and the manuscript has been read and approved by all named authors.

We would like to express my gratitude to all those who helped us during the writing of this article.

1. Arikan OA, Mulbry W, Rice C (2016) The effect of composting on the persistence of four ionophores in dairy manure and poultry litter. Waste Management 54: 110–117.
2. Wang ML (1996) Exact solution for a compound KdV-Burgers equations. Phys Lett A 213: 279–287.
3. Wang ML (1995) Solitary wave solutions for variant Boussinesq equations. Phys Lett A 199: 169–172.
4. Wang ML, Zhou YB, Li ZB (1996) Application of a homogeneous balance method to exact solutions of nonlinear evolution equations in mathematical physics. Phys Lett A 216: 67–75.
5. Wang, ML, Li X, Zhang J (2008) The expansion method and traveling wave solutions of nonlinear evolution equations in mathematical physics. Phys Lett A 372: 417–423.
6. Borhanifar A (2011) Application of the - expansion method for the Zhiber-shabat equation and other related equations. Math Comput Model 54: 2109-2116.
7. Borhanifar A, Kabir MM, Maryam Vahdat L (2009) New periodic and soliton wave solutions for the generalized Zakharov system and (2+1)-dimensional Nizhnik–Novikov–Veselov system. Chaos Solitons Fractals 42: 1646–1654.
8. Borhanifar A, Kabir MM (2009) New periodic and soliton solutions by application of exp-function method for nonlinear evolution equations. Comput Appl Math 229: 158–167.
9. Aslan İ (2009) Exact and explicit solutions to some nonlinear evolution equations by utilizing the-expansion method. Appl Math Comput 215: 857–863.
10. Bekir A (2008) Application of the (G′ /G) -expansion method for nonlinear evolution equations. Phys Lett A 372: 3400–3406.
11. Zhang J, Wei X, Lu Y (2008) A generalized expansion method and its applications. Phys Lett A 372: 3653–3658.
12. Alrtimi A, Rouainia M, Haigh S (2016) Thermal conductivity of a sandy soil. Appl Therm Eng 106: 551–560.
13. Prijono S, Laksmana MTS, Suprayogo D (2016) Effects of hedgerow systems on soil moisture and unsaturated hydraulics conductivity measured by the Libardi method. Journal of degraded and mining lands management 3: 491–498.
14. Zayed EME (2009) New traveling wave solutions for higher dimensional nonlinear evolution equations using a generalized (G′ /G) -expansion method. Journal of Physics A Mathematical and Theoretical 42: 5202– 5210.
15. Li LX, Li EQ, Wang ML (2010) The (G′ /G) -expansion method and its application to traveling wave solutions of the Zakharov equation. Applied Mathematics A Journal of Chinese Universities 25: 454–462.
16. Zhang J, Jiang FL, Zhao XY (2010) An improved (G′ /G) -expansion method for solving nonlinear evolution equations. International Journal of Computer Mathematics 87: 1716–1725.
17. Libardi PL, Reichardt K (2001) Libardi’s method refinement for soil hydraulic conductivity measurement. Aust J Soil Res 39: 851–860.
18. Libardi PL, Reichardt K, Nielsen DR, Biggar JW (1980) Simple field methods for estimating soil hydraulic conductivity. Soil Science Society of America Journal 44: 3–7.
19. Vasssilis ZA, Wyseure GCL (1998) Modeling of water and nitrogen dynamics on an undisturbed soil and a restored soil after open-cast mining. Agric Water Manag 37: 21–40.