
MBW:Observability of Trophic Food ChainsFrom MathBioContentsExecutive SummaryIn the paper "Observer Design for Open and Closed Trophic Chains" by Z. Varga, M. Gamez, and I. Lopez the authors explore mathematical techniques for observing systems from partial information. A simple model of a tropic food chain is created, containing three species. Under the assumption that only one of these species can be observed, the authors demonstrate the circumstances under which it is possible reconstruct the state of all species in the model. OverviewMathematics UsedThe authors for this paper used dynamic modeling and the mathematics behind control theory. Along with this, Jacobian's were used to find steady states of the tropic states. Types of ModelThis paper uses a nonLotka Volterra type trophic system which means that it does not simply follow a predatorprey relationship, but rather is multi leveled and contains open and closed chains. For a discussion on Lotka Volterra type population system see: to The Dynamics of Arthropod PredatorPrey Systems Biological System StudiedThe biological system studied was a Trophic chain that involved three levels. The 0th level in this chain was the "resource", the 1st tropic level was the "producer" and the 3rd trophic level was the "primary consumer". In the closed form, dead individuals from the 1st and 2nd tropic level could be recycled back into the 0th tropic level where as in the open form there is no recycling. ContextThe motivating example is ecology, where it may be practical to monitor only a small subset of the ecosystem. Their ecological model was a trophic chain. A trophic chain is a food chain where there is a linear food chain, where each species feeds on the previous. All species die and replenish some portion of a resource, which sustains the creature at the bottom. In the author's model, there is a resource, producer and primary consumer. For example, nutrients, plants, and herbivores would fit this model. If an ecologist can construct a mathematical model of how the system behaves, physical observations may give insight into other parts of the system. The mathematics behind this procedure are borrowed from the field of control theory. An extension is if a system can be observed, it can also be controlled, although modifying the system based on observation goes beyond the scope of this paper. For another study on a trophic system with two competing predators see: Three Species Food Chain Modeling HistoryThere have been a number of publications on the application of observer design and control theory to ecology. An example is a previous paper, by Lopez, Gamez and Molnar called "Observability and observers in a food web", where they analyzed a similar problem. They looked at a simple food web. There have been application beyond ecology also by Varga, Gamez and Lopez called "Observer design for phenotypic observations of genetic processes". ModelThe model is based on three elements of a trophic chain, . These represent the nutrient, producer, and consumer.
 The malthus parameters. These force exponential decay in the absence of food.  rate at which species 1 consumes 0, rate species 2 consumers 1.  rate at which species 1, 2 increase biomass from their consumption  rate at which species 1, species 2 die and recycle to nutrient
AnalysisThe authors define two matrices:. Where A is a system of n equations. C is the observation, which is found from the definition of where i is a species in the chain. For their model, they chose which observes the resource. The asterisk denotes a nontrivial steady state. This can be found in a tedious but straightforward manner. They also define a matrix H which maps the observation to a matrix of the same dimensions as A. This also assumes that . A system can be reconstructed from observations of a single system if assuming the observation is near equilibrium. This is tedious result from control theory, that the authors use without proving. AHC must also have negative and real eigenvalues. The authors use the criteria that the characteristic polynomial should satisfy the Hurwitz criteria. The matrix A, after linearization of the system is: The authors then find the equilibrium state and linearize the system There are solutions if: This simply says that if there are too few resources, the food chain may not have a stable steady state. They then state a result from control theory. They can construct a local exponential observer: where y=h(x). A way to think about this is that an observation near equilibrium (mapped to the system by H) will exponentially decay to the true state. The authors proceed to graph the trajectories, however they graph them as a phase plot. I have graphed the solutions of the first trajectory, as a function of time. It is clear that even though the trajectories have slightly different conditions, all populations converge to the same steady state. The authors then use the same procedure to demonstrate trajectories based on observation of all three species are possible. Useful MathThere is a useful mathematical technique to place the eigenvalues of a matrix of the for ABK where A and B are known, and K is to be determined to give the desired eigenvalues. It is possible to write out all the eigenvalue equations of a system by writing: Instead of individual eigenvalues is a vector of all eigenvectors , and is a matrix with corresponding eigenvalues on the diagonals. This can be rearranged as Which can be rewritten as a block matrix equation: Rewrite again by replacing the vector by z: The vector z is in the nullspace of the matrix. It can be written as . This allows a solution for K, as . By multiplying the vector out to obtain the equation it is possible to solve for K, and select whatever eigenvalues you want: Connection to Control TheoryThe authors do not explicitly explain the connection to control theory or the motivation for equation (2). The first step is to write out the linearized system: Then write out another system, that is an estimation of the system: In system the H matrix is multiplied by the difference between the true and estimated observation of some species. Subtracting system from results in: The difference in x can be renamed as z, and the result is the local observer. Substituting in the linearized definition of y, the system becomes:
It would be great if asymptotically decayed to the true system. It is for this reason it is desired that AHC has negative real eigenvalues. Pole placement can be used to achieve this if desired. I will briefly comment on what this has to do with control theory, as Robustness Analysis of an ObserverBased Controller in a Food Web gives a good description of control theory. A system you want to control can be described as where A is the linearized model of the system. B is known, and K is chosen to make the system behave as desired. The criteria (1) will determine whether the system is controllable. This is the same criteria that determines if a system is observable. ConclusionThe authors have demonstrated that given complete information about one species, as well as a mathematical model of the interactions between the species, it is possible to determine the state of all the species in the system, under certain criteria. Using this in practice may be more difficult. Even if the model is known, determining the steady state may be challenging. The author also does not comment on over how quickly the system asymptotically decays to the true solution. Discussion of Recent CitationThe paper entitled Monitoring in a predatorprey systems via a class of high order observer design cites this paper and is by Juan Luis MataMachuca, Rafael MartinezGuerra, and Ricardo Aguilar cites this paper in particular for it's use of high order observation in order to construct a model of a predatorprey relationship. The paper explains that observer design allows for a system to be modelled in various different fields and gives examples such as population genetics, evolutionary dynamics and neuralnetworks. This paper uses a high order polynomial to understand a predatorprey relationship. The paper used various resources such as theoretical model as well as numerical results. Then this methodology was used on other various Lotka Volterra models with up to three species with success. External linksThe paper: Observer design for open and closed trophic chains 