It will be shown, that Mitscherlich's formula, the crop yield as function of the fertilizer, can be partitioned into two parts. One comes from the fertilizer in the soil, the other from the fertilizer, supplied from outside. As a consequence, the quantity of the fertilizer, which is extracted from the soil by the crop, can be computed. Correspondingly new fertilizing can be planned.
One problem in fertilizing, on which many authors have worked, recently was reported on by Giebel et al. (2006): "Until now one may be able to derive a prognosis upon the spatial average of soil mineral-N content of a field. However, no approaches are available that would allow the estimation of the spatial distribution of plant-available N over an entire field within the same or a following year. Nevertheless, adequate representation of spatial N distribution within a field is important to site-specific and environmentally sustainable N fertilization". This result of spatial averages of soil mineral-N content of a field was obtained by Giebel et al. with samples of ground-drillings and a following analysis.
Another problem is the question, what part of the crop-yield comes from the plant-available fertilizer in the soil, and what part is caused by the given fertilizer. This problem is solved here in quite a different way without technical devices, only with analytical methods and the fact, that the Mitscherlich- process can be partitioned into two exponential processes. One of these processes is the so-called standardized Mitscherlich process, the other is an exponential dying process. From the quantity of lost fertilizer in the soil consequences for future fertilizing can be drawn.
The Law of Mitscherlich (1909) is an exponential growing process with horizontal asymptote. This sort of processes is widely known in natural sciences: In biology we find it as Mendel's law of genetics, in mechanics it is the generalization of Hooke's law (Schneeberger 2005): the first (linear) term of the Taylor-series of this generalization is Hooke's law; in electrical engineering charging of a battery follows this law, etc.
Mitscherlich first formulated the law in form of a differential equation:
where x is the quantity of fertilizer, y the crop-yield and (as usual in statistics) the hypothetical value of y; a is the maximum value of the yield (we exclude over-fertilization!) and b>0 is the most important exponential factor. Formula (1) means, that the growing of the crop is proportional to (a-). With boundary condition (x = 0) = c we have as solution of equation (1)
The boundary condition (x = 0) = c means, that c is the crop-yield with only the plant-available fertilizer in the soil, with no supplementary fertilizer, i.e. x = 0. See the curve in figure 1.
For our analysis we use the results of an experiment with winter-wheat at the University of Technology München/Weihenstephan in 2002, given in table 1: For every of the r=15 values xr of fertilizer, s=4 values of crop-yield yrs are given. For clarity in the figures not the experimental points (xr, yrs) with the great variance of the y-values, but (xr, ) are given as stars; see also the means in table 1.
The three parameters a,b,c of the Mitscherlich-curve (2) are estimated with the method of Least Squares of Gauss, i.e. a, b and c are solutions of the condition
This problem is solved iteratively with the non-linear Simplex-Method of Nelder and Mead (1965). The result is
It can easily be shown, that computation with the values of (instead of y) yields the same result.
The curve , plotted in figure 1, is extrapolated to negative values of x (dotted line). We get x = d = -83.8 as solution of (x) = 0. -d(>0) is the plant-available fertilizer in the soil, which gives the corresponding crop-yield =c=53.2, i.e. without fertilizing from outside. With this soil-immanent fertilizer -d we can write formula (2) also in the form of the Mitscherlich-Baule function (Baule 1918). See also Frank, Beattie and Embleton (1990):
Figure 1. Partition of the Mitscherlich-curve into the two parts 01 and 02
The following theorem immediately results from formula (6):
Theorem: The general Mitscherlich-curve (x) of formula (6) with -d>0, or formula (2) with c>0, originates from the "standardized" Mitscherlich curve 01(x)=a(1-e-bx) with d=0 in formula (6) or c=0 in formula (2) by simply shifting the curve 01(x) to the left for -d >0 units (of variable x). See figure 1.
For the strategy of fertilizing in the future, to avoid too much or too less fertilizing with the known consequences, it is important to state, that according to Liebig's Law we have for the total quantity q of plant-available fertilizer in the soil: q≥-d
We can write formula (2) in the form
The parts 01 and 02 are of specific importance:
is the standardized Mitscherlich process, beginning at the origin (x,)=(0,0), i.e. c=0. So this process is generated alone by the fertilizer, supplied from outside.
is a diminishing exponential process, starting from 02(x=0)=c. So 02(x) is the crop-yield, caused by the plant-available fertilizer in the soil alone. Note the remark to formula (2)!
Remark: Processes of the form (b) are the most important "dying processes" in natural sciences (e.g. radioactive decay), biology , medicine, economics, etc. Here one better would speak of an exhaustion process (of the fertilizer in the soil).
In figure 1 the curves 01 and 02 are plotted as dashed lines.
If the quantity of fertilizer from outside is zero (i.e. x=0), we have 01=0 and the crop is = 02=c=53,2 (100 kg/ha of winter-wheat), which comes alone from the fertilizer in the soil -d=83.8 (kg/ha of N).
Now we assume, that for example a quantity of fertilizer x = x0=100 kg/ha of N is given from outside. Then we get from formula (9) the quantity of crop 01(x0), which comes from the given x0 and with formula (10) we get the quantity of crop 02(x0), which comes from the fertilizer in the soil. The total crop is (x0)=01(x0)+02(x0) We get
01(x0)=60,1 | 02(x0)=24,9 | (x0)=85,0
The problem now is: What quantity of fertilizer in the soil gives the crop 02(x0)? This is found with the auxiliary Mitscherlich process
which yields x=d1=-32,8 as solution of 1(x)=0. -d1=32.8 (kg/ha) is the soil-immanent plant-available part of fertilizer, which, together with the fertilizer x0=100 (kg/ha) from outside, yields the total crop (x0)=85,0 (100 kg/ha). See figure 2.
Note: Equation (11) can also be written in the equivalent form:
Figure 2: Example of fertilizing with x0=100 kg/ha of N. Calculation of the quantity of fertilizer -d1=32,8 kg/ha of N, supplied by the soil
A necessary supplement follows with Paper 2 (Paper 2: Mitcherlich's Law: A Supplement)
I have to thank Dr. Maidl of the University of Technology München (Weihenstephan)/Germany for providing me with the data, used in this article.
Baule, B. (1918). Zu Mitscherlichs Gesetz der physiologischen Beziehungen, Landwirtschaftliche Jahrbücher 51, 363-385.
Frank M.D.,Beattie B.R. and Embleton M.E. (1990). A Comparison of Alternative Crop Response Models, American Journal of Agricultural Economics, 72, 597-603.
Giebel, A. et al.(2006). How representatively can we sample soil mineral nitrogen? Journal of Plant Nutrition and Soil Sciences, 169, 52-59.
Mitscherlich, E.A.(1909). Das Gesetz des Minimums und das Gesetz des abnehmenden Bodenertrags, Landwirtschaftliche Jahrbücher, 38, 537-552.
Nelder, J.R. and Mead, R. (1965). A Simplex Method for function minimization. The Computer Journal, 7, 303-313.
Schneeberger, H. (2005). Theory of the stress-strain relationship of concrete and steel. In R.K.Dhir, M.D.Newlands and L.J.Csetenyi (eds.): Applications of Nanotechnology in Concrete Design. Proceedings of the International Conference at the University of Dundee, Scotland , UK, 105-112.
Download this Paper in PDF format: Mitscherlich's Law: Sum of two exponential Processes - Conclusions (PDF, 0,5 MB)