Intermittent Plurisink Model and the Emergence of Complex Heterogeneity Patterns: A Simple Paradigm for Explaining Complexity in Soil Chemical Distributions

The spatial complexity of the distribution of organic matter, chemicals, nutrients, and pollutants has been demonstrated to have multifractal nature. This fact supports the possibility of existence of some emergent heterogeneity structure built under the evolution of the system. The aim of this paper is providing a consistent explanation of the mentioned results via an extremely simple model.


Introduction: Searching Explanations for Soil Heterogeneity
Heterogeneity and complexity are ubiquitous at all scales in soil and hydrologic systems.Nowadays, new technologies are of an invaluable help for providing a great number of highly calibrated field measurements.One can get a huge amount of data from computer tomography of soil samples at microscopic scales, digital terrain catchments of landscapes, and river basins among many other technological tools.Then mathematical tools are needed to analyse and interpret those data as well to construct models to predict.However, along the way needed to get such a final purpose, scientists also need to understand why the heterogeneity is produced and what the organizing principles that might underlie the heterogeneity and complexity are (McDonnell et al. [1]).Also it is encouraged to explore the scaling behaviour of heterogeneity and the emergent properties in soil and hydrologic systems.
In this paper we are mainly interested in some aspects concerning the heterogeneity in the soil scenario.
Certainly we believe that the above understanding, besides providing coherence to science, also may be useful to get the practical purpose itself.In the case of soil, an illustrating example supporting this and the issues mentioned above is the study of soil texture heterogeneity.On one hand, Multifractal Analysis of fine granulometry soil data obtained by laser diffraction techniques provides information about the scaling behaviour of particle size distribution (PSD) heterogeneity (Montero [2]).In a second step, models able to replicate the heterogeneity formerly shown may be useful for prediction purposes (Martín and García-Gutiérrez [3]).The answer to why such heterogeneity exists, however, is not an easy issue since different sources of heterogeneity should be expected.In Frisch and Sornette [4] and Sornette [5], it is suggested that the fractal behaviour might be the result of a natural mixing of simple multiplicative process that takes place along the fragmentation of different particles, also pointing out that there is no accepted theoretical explanation.Recently fragmentation algorithms were proposed to replicate the multifractal nature of soil PSD (Martín et al. [6]).In 2 Journal of Chemistry this respect it is needed to say that any partial, but coherent, explanation should help to understand the (possibly several ones) organizing principles involved.
On the other hand, the spatial complexity of the distribution of organic matter, chemicals, nutrients, and pollutants has been studied by different authors (Kravchenko et al. [7], Lehmann et al. [8]).Multifractal Analysis has been successfully used to study the spatial variability of chemicals and organic matter contents, which is characterized by the generalized fractal dimensions (Kravchenko et al. [7]).Searching why such structured heterogeneity exists, a reducionist approach based in the description of transport equations in soil, seems an unlikely choice to describe the emergence of such complex pattern across the spatial scales.On the contrary, an explanation based on the fact that many complex systems in nature evolve in an intermittent burst-like way rather than in a smooth gradual manner (Rodríguez-Iturbe and Rinaldo [9]) would be more adequate.Further, such a kind of structured heterogeneity is commonly interpreted as the result of chaos or self-organization which leads to the emergent structure built under the evolution of the system (Sornette [5]).The aim of this paper is to provide a small contribution via an extremely simple model, which gives a consistent explanation to the mentioned results on spatial variability of chemicals or pollutants in soil.
The paper is organized as follows.In Section 2.1 the model is presented and in Section 2.2 the entropy scaling analysis method is described.Section 3 is devoted to analysing the results obtained in different simulations and their discussion.

The Model.
Let us suppose that  is a soil area squareshaped.Suppose further that at any of the four corners there is a sink  ( = 1, 2, 3, 4) randomly acting in an intermittent manner.Suppose each sink  acting with relative frecuency   .A pollutant deposit ("pollutant seed") is supposedly located in an arbitrary point of the square.When a given sink  acts, its suction action is able to attract the pollutant matter to another point reducing the distance to the sink in a factor   < 1, where the pollutant rests until a new (or the same) sink acts.This factor reflects the mean value of the suction power of the respective sinks.However, the "flying" pollutant matter leaves a unit of pollutant at any point where the pollutant "rests" along its travelling.
Although a much more sophisticated model might be constructed for a more realistic performance under the same essential idea, we rather prefer to emphasize how complexity may appear under quite simple and natural actions evolving in time.

Measuring Heterogeneity.
When the model is implemented a first goal is applying mathematical tools in order to parameterize heterogeneity in a reliable manner.
For simplicity let us assume that the unit square  in the plane is the support of a distribution  with highly heterogeneous features.In order to scrutinize its heterogeneity, let us consider a collection (mesh) of 2  × 2  -boxes,   = {  :  = 1, 2, . . ., 2 2 }, of side length  = 2 − , representing a partition of  for each value ,  = 1, 2, 3, . . .(see Figure 1).
When the mass (  ) inside any box   is known, the Shannon entropy (Shannon [10]) of  with respect to a fix partition   is given by provided (  ) log (  ) = 0 if (  ) = 0.
The number   (  ) is expressed in information units (bits) and its extreme values are log 2 2 , which corresponds to the most even (homogeneous) case-where all the squares have the same cumulative mass-and 0, which corresponds to the most uneven (heterogeneous) case-where the whole mass is concentrated in a single square.The Shannon entropy   (  ) is a widely accepted measure of evenness or heterogeneity in the mass distribution  at the scale level given by each partition   .In fact, it can be shown that any measure of heterogeneity with the natural properties for such goal must be a multiple of   (  ) (Khinchin [11]).
Using increasing values of  (decreasing values of ) one can obtain an increasing amount of information about the distribution as   (  ) grows to infinity.If such an increase is not erratic but rather conforms to a scaling or asymptotic behaviour of   (  ) when  ↓ 0, then the entropy or information dimension of  is defined (Rényi [12]) by means of the equation where "≈" means that −  (  ) will linearly fit log .

Results and Discussion
In order to implement the model, different set of values of   and   were selected ( = 1, 2, 3, 4).First close values of   were used under the assumption of similar intermittent frequencies, while the   values used try to investigate the effect of relative different suction powers.For any simulation the centre of the square has been chosen as initial position for the "pollutant seed." Then for any simulation a scaling entropy analysis has been made following Section 2.2.Figures 2(a) and 2(b) show two different simulations of 500 points for the same   and   values ( 1 = 0.29,  2 = 0.21,  3 = 0.29,  4 = 0.21,  1 = 0.7,  2 = 0.5,  3 = 0.7, and  4 = 0.5).The scaling analysis was made by using values  = 2 − from  = 1 to  = 6.The mass (  ) is given by the proportion of points inside any box   .The value of   (  ) is plotted against − log  and a linear fitting is implemented.The slope of the regression line gives an estimation of the entropy dimension with  2 value as coefficient of determination.It can be noticed that the physical appearance of both simulations is quite different, thus illustrating the high influence of the random effect in this case.Also the scaling analysis reveals different results ( and  2 values) for both simulated distributions.For the same   and   values, simulation of 20000 points leads to the results in Figures 3(a    Table 1 shows the results of this analysis.It is observed that the influence of the random component diminishes for increasing number of points used in the simulation.Also the     2 0,5-0,5-0,5-0,5 1,991 0,9999 0,25-0,25-0,25-0,25 0,7-0,7-0,5-0,5 1,951 0,9999 0,8-0,6-0,7-0, 6 1,924 0,9997  2 values become closer to 1. Figure 4 shows the value of the estimated entropy dimension for increasing number of points.Table 2 shows data involved in that figure. Results clearly show the emergence of a mass distribution with a well-defined structured heterogeneity that the scaling analysis reveals.In fact the robustness of the results is based on a theorem of ergodic type (Elton [13]).
Finally Table 3 shows the value of the estimated entropy dimension.The  2 values obtained reflect the scale invariance of the resulting distributions.
Smaller   values representing greater suction powers have obvious influence on the heterogeneity of the final distribution which remains parameterized by the entropy dimension.In an intuitive sense, the entropy dimension value may be interpreted as uncertainty degree.In fact it can be used together with other parameters in interpolation procedures in soil spatial variability studies (Kravchenko et al. [7]).

Conclusions
Heterogeneity is ubiquitous in many soil scenarios.In particular the spatial complexity of the distribution of organic matter, chemicals, nutrients, and pollutants is a frequent ingredient, which is in the focus of soil studies.The understanding of why the heterogeneity is produced, and what the nature of such heterogeneity is, is a need under the scientific and practical points of view.Any coherent explanation on the origin of heterogeneity should help to understand it and to choose the adequate mathematical techniques for handling it with prediction purposes.
In this paper an extremely simple model is presented, which gives a consistent explanation of the complexity of spatial variability of chemicals or pollutants in soil shown in former studies.
The results shown here strongly suggest the use of scaling methods coming from fractal geometry for the study of this kind of distributions.

Figure 1 :
Figure 1: Partition of the support  by squares of side length .