Mitscherlich s Law: Sum of two exponential Processes; Conclusions 2009, 1 st July

Similar documents
Pearson s Chi-Square Test Modifications for Comparison of Unweighted and Weighted Histograms and Two Weighted Histograms

3.1 Random variables

Analytical Solutions for Confined Aquifers with non constant Pumping using Computer Algebra

ELASTIC ANALYSIS OF CIRCULAR SANDWICH PLATES WITH FGM FACE-SHEETS

F-IF Logistic Growth Model, Abstract Version

ON INDEPENDENT SETS IN PURELY ATOMIC PROBABILITY SPACES WITH GEOMETRIC DISTRIBUTION. 1. Introduction. 1 r r. r k for every set E A, E \ {0},

( ) [ ] [ ] [ ] δf φ = F φ+δφ F. xdx.

A NEW VARIABLE STIFFNESS SPRING USING A PRESTRESSED MECHANISM

Multiple Criteria Secretary Problem: A New Approach

A new approach in classical electrodynamics to protect principle of causality

Title. Author(s)Y. IMAI; T. TSUJII; S. MOROOKA; K. NOMURA. Issue Date Doc URL. Type. Note. File Information

Nuclear size corrections to the energy levels of single-electron atoms

Stress Intensity Factor

CALCULATING THE NUMBER OF TWIN PRIMES WITH SPECIFIED DISTANCE BETWEEN THEM BASED ON THE SIMPLEST PROBABILISTIC MODEL

Liquid gas interface under hydrostatic pressure

Absorption Rate into a Small Sphere for a Diffusing Particle Confined in a Large Sphere

Deflection of light due to rotating mass a comparison among the results of different approaches

MAGNETIC FIELD AROUND TWO SEPARATED MAGNETIZING COILS

COMPUTATIONS OF ELECTROMAGNETIC FIELDS RADIATED FROM COMPLEX LIGHTNING CHANNELS

On the integration of the equations of hydrodynamics

Euclidean Figures and Solids without Incircles or Inspheres

Application of Parseval s Theorem on Evaluating Some Definite Integrals

Gradient-based Neural Network for Online Solution of Lyapunov Matrix Equation with Li Activation Function

Solving Some Definite Integrals Using Parseval s Theorem

Application of homotopy perturbation method to the Navier-Stokes equations in cylindrical coordinates

Bounds on the performance of back-to-front airplane boarding policies

FREE TRANSVERSE VIBRATIONS OF NON-UNIFORM BEAMS

Lecture 8 - Gauss s Law

Duality between Statical and Kinematical Engineering Systems

Mathematical Model of Magnetometric Resistivity. Sounding for a Conductive Host. with a Bulge Overburden

New problems in universal algebraic geometry illustrated by boolean equations

An Application of Fuzzy Linear System of Equations in Economic Sciences

This is a very simple sampling mode, and this article propose an algorithm about how to recover x from y in this condition.

Δt The textbook chooses to say that the average velocity is

Conservative Averaging Method and its Application for One Heat Conduction Problem

APPLICATION OF MAC IN THE FREQUENCY DOMAIN

working pages for Paul Richards class notes; do not copy or circulate without permission from PGR 2004/11/3 10:50

Auchmuty High School Mathematics Department Advanced Higher Notes Teacher Version

On Polynomials Construction

EM Boundary Value Problems

A Comparison and Contrast of Some Methods for Sample Quartiles

A New Approach to General Relativity

3.6 Applied Optimization

Hammerstein Model Identification Based On Instrumental Variable and Least Square Methods

Central Coverage Bayes Prediction Intervals for the Generalized Pareto Distribution

Chapter 5 Linear Equations: Basic Theory and Practice

Section 8.2 Polar Coordinates

A dual-reciprocity boundary element method for axisymmetric thermoelastodynamic deformations in functionally graded solids

A New Method of Estimation of Size-Biased Generalized Logarithmic Series Distribution

A Backward Identification Problem for an Axis-Symmetric Fractional Diffusion Equation

On a quantity that is analogous to potential and a theorem that relates to it

Scattering in Three Dimensions

Surveillance Points in High Dimensional Spaces

CHAPTER 3. Section 1. Modeling Population Growth

Using Laplace Transform to Evaluate Improper Integrals Chii-Huei Yu

SOLVING THE VISCOUS COMPOSITE CYLINDER PROBLEM BY SOKOLOV S METHOD

Numerical Integration

International Journal of Mathematical Archive-3(12), 2012, Available online through ISSN

Magnetometer Calibration Algorithm Based on Analytic Geometry Transform Yongjian Yang, Xiaolong Xiao1,Wu Liao

1 Similarity Analysis

An Exact Solution of Navier Stokes Equation

THE LAPLACE EQUATION. The Laplace (or potential) equation is the equation. u = 0. = 2 x 2. x y 2 in R 2

When two numbers are written as the product of their prime factors, they are in factored form.

Chaos and bifurcation of discontinuous dynamical systems with piecewise constant arguments

ac p Answers to questions for The New Introduction to Geographical Economics, 2 nd edition Chapter 3 The core model of geographical economics

6 PROBABILITY GENERATING FUNCTIONS

Energy Levels Of Hydrogen Atom Using Ladder Operators. Ava Khamseh Supervisor: Dr. Brian Pendleton The University of Edinburgh August 2011

On absence of solutions of a semi-linear elliptic equation with biharmonic operator in the exterior of a ball

What Form of Gravitation Ensures Weakened Kepler s Third Law?

7.2.1 Basic relations for Torsion of Circular Members

4/18/2005. Statistical Learning Theory

Study on Application of New Theorem of Kinetic Energy

Hypothesis Test and Confidence Interval for the Negative Binomial Distribution via Coincidence: A Case for Rare Events

To Feel a Force Chapter 7 Static equilibrium - torque and friction

Coupled Electromagnetic and Heat Transfer Simulations for RF Applicator Design for Efficient Heating of Materials

TESTING THE VALIDITY OF THE EXPONENTIAL MODEL BASED ON TYPE II CENSORED DATA USING TRANSFORMED SAMPLE DATA

Estimation of the Correlation Coefficient for a Bivariate Normal Distribution with Missing Data

A matrix method based on the Fibonacci polynomials to the generalized pantograph equations with functional arguments

Math 2263 Solutions for Spring 2003 Final Exam

Physics 2A Chapter 10 - Moment of Inertia Fall 2018

Web-based Supplementary Materials for. Controlling False Discoveries in Multidimensional Directional Decisions, with

Question 1: The dipole

PROBLEM SET #1 SOLUTIONS by Robert A. DiStasio Jr.

Solution of a Spherically Symmetric Static Problem of General Relativity for an Elastic Solid Sphere

Modeling of High Temperature Superconducting Tapes, Arrays and AC Cables Using COMSOL

Solving Problems of Advance of Mercury s Perihelion and Deflection of. Photon Around the Sun with New Newton s Formula of Gravity

MEASURING CHINESE RISK AVERSION

Problem 1. Part b. Part a. Wayne Witzke ProblemSet #1 PHY 361. Calculate x, the expected value of x, defined by

Research Article On Alzer and Qiu s Conjecture for Complete Elliptic Integral and Inverse Hyperbolic Tangent Function

2. Electrostatics. Dr. Rakhesh Singh Kshetrimayum 8/11/ Electromagnetic Field Theory by R. S. Kshetrimayum

Recent Advances in Chemical Engineering, Biochemistry and Computational Chemistry

Electrodynamic Forces in Steel Strip during Induction Heating

Bayesian Analysis of Topp-Leone Distribution under Different Loss Functions and Different Priors

Berkeley Math Circle AIME Preparation March 5, 2013

Numerical Inversion of the Abel Integral Equation using Homotopy Perturbation Method

Math 1105: Calculus I (Math/Sci majors) MWF 11am / 12pm, Campion 235 Written homework 3

arxiv: v1 [math.co] 4 May 2017

On the Comparison of Stability Analysis with Phase Portrait for a Discrete Prey-Predator System

Classical Worm algorithms (WA)

Why Professor Richard Feynman was upset solving the Laplace equation for spherical waves? Anzor A. Khelashvili a)

Transcription:

Mitschelich s Law: Sum of two exponential Pocesses; Conclusions 29, st July Hans Schneebege Institute of Statistics, Univesity of Elangen-Nünbeg, Gemany Summay It will be shown, that Mitschelich s fomula, the cop yield as function of the fetilize, can be patitioned into two pats. One comes fom the fetilize in the soil, the othe fom the fetilize, supplied fom outside. As a consequence, the quantity of the fetilize, which is extacted fom the soil by the cop, can be computed. Coespondingly new fetilizing can be planned. Intoduction One poblem in fetilizing, on which many authos have woked, ecently was epoted on by Giebel et al. (26): Until now one may be able to deive a pognosis upon the spatial aveage of soil mineal-n content of a field. Howeve, no appoaches ae available that would allow the estimation of the spatial distibution of plant-available N ove an entie field within the same o a following yea. Nevetheless, adequate epesentation of spatial N distibution within a field is impotant to site-specific and envionmentally sustainable N fetilization. This esult of spatial aveages of soil mineal-n content of a field was obtained by Giebel et al. with samples of gound-dillings and a following analysis. Anothe poblem is the question, what pat of the cop-yield comes fom the plant-available fetilize in the soil, and what pat is caused by the given fetilize. This poblem is solved hee in quite a diffeent way without technical devices, only with analytical methods and the fact, that the Mitschelich- pocess can be patitioned into two exponential pocesses. One of these pocesses is the so-called standadized Mitschelich pocess, the othe is an exponential dying pocess. Fom the quantity of lost fetilize in the soil consequences fo futue fetilizing can be dawn. The Law of Mitschelich (99) is an exponential gowing pocess with hoizontal asymptote. This sot of pocesses is widely known in natual sciences: In biology we find it as Mendel s law of genetics, in mechanics it is the genealization of Hooke s law (Schneebege 25): the fist (linea) tem of the Taylo-seies of this genealization is Hooke s law; in electical engineeing chaging of a battey follows this law, etc. Mitschelich fist fomulated the law in fom of a diffeential equation: d / dx = ' = b(a ) () whee x is the quantity of fetilize, y the cop-yield and (as usual in statistics) the hypothetical value of y; a is the maximum value of the yield (we exclude ove-fetilization!) and b> is the most impotant exponential facto. Fomula () means, that the gowing of the cop is popotional to (a- ). With bounday condition = ) = c we have as solution of equation () bx = c + (a c)( e ) (2)

The bounday condition = ) = c means, that c is the cop-yield with only the plantavailable fetilize in the soil, with no supplementay fetilize, i.e. x=. See the cuve in figue. MATERIALS AND METHODS Fo ou analysis we use the esults of an expeiment with winte-wheat at the Univesity of Technology München/Weihenstephan in 22, given in table : Fo evey of the =5 values x of fetilize, s=4 values of cop-yield y s ae given. Fo claity in the figues not the expeimental points ( x, y s ) with the geat vaiance of the y-values, but ( x, y ) ae given as stas; see also the values y in table. Table Yield y (winte-wheat in kg/ha) fo given fetilize x (N in kg/ha) 2 3 4 5 6 7 8 s x 8 2 4 6 6 8 y s 4.8 73.8 79.6 84.5 97.7 82.9 96.2 2.6 2 58.6 89.9 94.7 78. 5.5 89.4 92.7 2. 3 62.7 76.4 84.4 78.2 8.4 9. 97.6 3.4 4 5.8 79.3 88.5 94. 93.4 96.4 93.5 92.5 y 53.475 79.85 86.8 83.725 94.25 94.425 95..5 continuation 9 2 3 4 5 s x 8 2 2 22 22 24 26 y s 2.8 96.4 3.5 98.6 2.7 2.2 8. 2 92. 2. 94.8 3.5 5. 4.7 6.5 3 4.7 2. 4.8 93.7 4. 97.7. 4 96. 8.4 3. 7.4 4.6 78. 3.3 y 98.925 2.225.55.8 4.75 95.65 6.975 The thee paametes a,b,c of the Mitschelich-cuve (2) ae estimated with the method of Least Squaes of Gauss, i.e. a, b and c ae solutions of the condition 2 2 ( y s (c + (a c)( e ))) = (y s ) s s Min (3) This poblem is solved iteatively with the non-linea Simplex-Method of Nelde and Mead (965). The esult is a=3., b=.7583, c=53.2 (a c=59.9) (4).7583 x = 53.2 + 59.9( e ) (5) It can easily be shown, that computation with the values of y (instead of y) yields the same esult. The cuve, plotted in figue, is extapolated to negative values of x (dotted line). We get x = d = 83.8 as solution of ) =. d(>) is the plant-available fetilize in the soil, which

gives the coesponding cop-yield =c=53.2, i.e. without fetilizing fom outside. With this soil-immanent fetilize d we can wite fomula (2) also in the fom of the Mitschelich- Baule function (Baule 98). See also Fank, Beattie and Embleton (99): with d = b d) = a( e ) (6) ln b a a c (<) (7) Figue. Patition of the Mitschelich-cuve into the two pats and 2 The following theoem immediately esults fom fomula (6): Theoem: The geneal Mitschelich-cuve ) of fomula (6) with d >, o fomula (2) with c>, oiginates fom the standadized Mitschelich cuve ) = a( e ) with d= in fomula (6) o c= in fomula (2) by simply shifting the cuve ) to the left fo d > units (of vaiable x). See figue. Fo the stategy of fetilizing in the futue, to avoid too much o too less fetilizing with the known consequences, it is impotant to state, that accoding to Liebig s Law we have fo the total quantity q of plant-available fetilize in the soil: q d. Patition of Mitschelich s Fomula We can wite fomula (2) in the fom = a( e ) + ce = + 2 (8)

The pats and 2 ae of specific impotance: a) =a(-e ) (9) is the standadized Mitschelich pocess, beginning at the oigin, )=(,), i.e. c=. So this pocess is geneated alone by the fetilize, supplied fom outside. b) 2 ce = () is a diminishing exponential pocess, stating fom 2 =)=c. So 2 ) is the cop-yield, caused by the plant-available fetilize in the soil alone. Note the emak to fomula (2)! Remak: Pocesses of the fom (b) ae the most impotant dying pocesses in natual sciences (e.g. adioactive decay), biology, medicine, economics, etc. Hee one bette would speak of an exhaustion pocess (of the fetilize in the soil). In figue the cuves and 2 ae plotted as dashed lines. RESULTS AND DISCUSSION If the quantity of fetilize fom outside is zeo (i.e. x=), we have = and the cop is = 2 = c = 53.2 ( kg/ha of winte-wheat), which comes alone fom the fetilize in the soil d=83.8 (kg/ha of N). Now we assume, that fo example a quantity of fetilize x = x = kg/ha of N is given fom outside. Then we get fom fomula (9) the quantity of cop o ), which comes fom the given x and with fomula () we get the quantity of cop 2 ), which comes fom the fetilize in the soil. The total cop is ) = ) ). We get + 2 ) = 6.; 2 ) = 24.9, ) = 85. The poblem now is: What quantity of fetilize in the soil gives the cop 2 )? This is found with the auxiliay Mitschelich pocess = c + (a c )( e ); = ) 24. 9 () c 2 = which yields x= d = 32. 8 as solution of ) =. d=32.8 (kg/ha) is the soil-immanent plant-available pat of fetilize, which, togethe with the fetilize x = (kg/ha) fom outside, yields the total cop ) = 85. ( kg/ha). See figue 2. Note: Equation () can also be witten in the equivalent fom: with b d ) = a( e ) (a) a d = ln < b a c

Figue 2: Example of fetilizing with x = kg/ha of N. Calculation of the quantity of fetilize d = 32. 8 kg/ha of N, supplied by the soil A necessay supplement follows with Pape 2: Mitchelich's Law: A Supplement, http://www.soil-statistic.de/mitschelich-supplement/mitschelich-supplement.html ACKNOWLEDGEMENT I have to thank D. Maidl of the Univesity of Technology München (Weihenstephan)/Gemany fo poviding me with the data, used in this aticle. REFERENCES Baule, B. (98). Zu Mitschelichs Gesetz de physiologischen Beziehungen, Landwitschaftliche Jahbüche 5, 363-385. Fank M.D.,Beattie B.R. and Embleton M.E. (99). A Compaison of Altenative Cop Response Models, Ameican Jounal of Agicultual Economics, 72, 597-63. Giebel, A. et al.(26). How epesentatively can we sample soil mineal nitogen? Jounal of Plant Nutition and Soil Sciences, 69, 52-59. Mitschelich, E.A.(99). Das Gesetz des Minimums und das Gesetz des abnehmenden Bodenetags, Landwitschaftliche Jahbüche, 38, 537-552. Nelde, J.R. and Mead, R. (965). A Simplex Method fo function minimization The Compute Jounal, 7, 33-33. Schneebege, H. (25). Theoy of the stess-stain elationship of concete and steel. In: R.K.Dhi, M.D.Newlands and L.J.Csetenyi (eds.): Applications of Nanotechnology in Concete Design. Poceedings of the Intenational Confeence at the Univesity of Dundee, Scotland, UK, 5-2.