Fast and flexible interpolation via PUM with applications in population dynamics

Alessandra De Rossi,Roberto Cavoretto,Emma Perracchione
DOI: https://doi.org/10.48550/arXiv.1512.03934
2015-12-13
Abstract:In this paper the Partition of Unity Method (PUM) is efficiently performed using Radial Basis Functions (RBFs) as local approximants. In particular, we present a new space-partitioning data structure extremely useful in applications because of its independence from the problem geometry. Moreover, we study, in the context of wild herbivores in forests, an application of such algorithm. This investigation shows that the ecosystem of the considered natural park is in a very delicate situation, for which the animal population could become extinguished. The determination of the so-called sensitivity surfaces, obtained with the new fast and flexible interpolation tool, indicates some possible preventive measures to the park administrators.
Numerical Analysis
What problem does this paper attempt to address?
The main problem that this paper attempts to solve is how to perform interpolation calculations efficiently and flexibly and apply them to population dynamics, especially to study the population dynamics of wild herbivores in forests. Specifically: 1. **Efficient and Flexible Interpolation Method**: The paper proposes a Partition of Unity Method (PUM) based on Radial Basis Functions (RBFs) and introduces a new spatial partitioning data structure - a block - based partitioning structure. This method can perform efficient interpolation calculations in regions of any geometric shape without depending on the specific geometric form of the problem. 2. **Population Dynamics Application**: The paper applies the above - mentioned interpolation method to study the wild herbivore population in a natural park in northern Italy. By constructing and analyzing the so - called "sensitivity surfaces", researchers can determine which small changes in the parameters of the ecosystem may increase the extinction risk of the herbivore population. ### Specific Problem Description - **Improvement of Interpolation Method**: Traditional methods such as kd - trees are less efficient when dealing with complex geometric shapes, while the block - based partitioning structure proposed in the paper can significantly improve computational efficiency. For example, when interpolating 107 sensitivity points, it only takes 0.31 seconds using the new method, while it takes 5.02 seconds using kd - trees. - **Ecosystem Stability Analysis**: Through numerical simulation, researchers have found that the current number of herbivore populations is very close to the separatrix surface between stability and extinction. This means that even a small change in environmental conditions may cause the herbivore population to enter the extinction area. ### Solution Suggestions To prevent the extinction of herbivores, the paper proposes the following measures: - **Reduce Mortality (\(\mu\))**: By establishing winter safe shelters, the mortality of herbivores can be reduced. - **Improve Nutrient Absorption (\(\alpha^{-1}\))**: Plant more nutritious grass plants to increase the nutrients that herbivores obtain from grass. These measures aim to keep the current parameter values away from the sensitivity surfaces, thereby maintaining the stability of the ecosystem. ### Summary The core contribution of the paper lies in developing an efficient and flexible interpolation tool and successfully applying it to the study of actual ecosystems, providing a scientific basis and feasible management suggestions for protecting wild herbivores.