Tsallis q-Stat and the Evidence of Long-Range Interactions in Soil Temperature Dynamics

The complexities in the variations of soil temperature and thermal diffusion poses a physical problem that requires more understanding. The quest for a better understanding of the complexities of soil temperature variation has prompted the study of the q-statistics in the soil temperature variation with the view of understanding the underlying dynamics of the temperature variation and thermal diffusivity of the soil. In this work, the values of Tsallis stationary state q index known as q-stat were computed from soil temperature measured at different stations in Nigeria. The intrinsic variations of the soil temperature were derived from the soil temperature time series by detrending method to extract the influences of other types of variations from the atmosphere. The detrended soil temperature data sets were further analysed to fit the q-Gaussian model. Our results show that our datasets fit into the Tsallis Gaussian distributions with lower values of q-stat during rainy season and around the wet soil regions of Nigeria and the values of q-stat obtained for monthly data sets were mostly in the range 1.2≤q≤2.9 for all stations, with very few values q closer to 1.2 for a few stations in the wet season. The distributions obtained from the detrended soil temperature data were mostly found to belong to the class of asymmetric q-Gaussians. The ability of the soil temperature data sets to fit into q-Gaussians might be due and the non-extensive statistical nature of the system and (or) consequently due to the presence of superstatistics. The possible mechanisms responsible this behaviour was further discussed.


Introduction
The complexities in soil thermal diffusion and the complex thermal interaction between the soil and the atmosphere could be responsible for the sporadic changes in the soil temperature measurements. Different factors are responsible for the variations in the soil temperature and these factors can contribute greatly to the complex variations of the soil temperature [1][2][3]. Dynamical complexity, a phenomenon which has been attributed to the dynamics of the atmospheric parameters, have been studied and characterised in different ways [4][5][6]. Many researchers have attributed the complexity of the atmosphere as contributing factors to the complexity of soil temperature variations [7][8][9]. However, the underlining dynamics of both systems are different and as such, the intrinsic complex dynamics of the soil temperature propagation and the various factors affecting the thermal diffusivity of the soil are the major determinants of its internal temperature dynamics. The complexities in soil dynamical processes in aspects such as the soil hydrology have been evaluated using the entropic approach and it has been able to reveal the variations of soil dynamics resulting from varying hydrological processes [10,11], which usually affects the temperature variation [1]. The other intrinsic factors affecting the soil thermal diffusivity include soil vegetation cover [1,3,8,12], soil colour, soil texture and soil organic carbon [1,3,7,8,12].
As described earlier, most atmospheric temperature systems are easily characterized by dynamical complexity. This framework is based on the study of complexity in these systems considering the Tsallis statistical distribution. Since statistical distributions are particularly important in the understanding the complexity of a dynamical system. The dynamical complexity of atmospheric and many environmental systems are mainly characterized by superstatistics, which can be seen as a generalization of many other distributions that may also characterize the atmospheric systems. For example, the symmetric q-Gaussian can be easily derived into a superstatistical distribution [13][14][15]. Superstatistics can describe fluctuations in natural systems that follows from a χ 2 -distribution like the q-Gaussian [14,16]. This explains the reason why heavy tailed distributions such as the atmospheric temperature and rainfall has been found to assume superstatistical characteristics [17,18]. Besides the symmetrical q-Gaussian, the superstatistical distributions can also be derived from the log-normal distributions [14], which have been found to characterize the atmospheric humidity and precipitation [19]. Relating the distributions describing natural systems such as the atmosphere to the presence of superstatistics could be a good reason for the possible presence of q-Gaussians. This because, like the ability to derive the superstatistics from q-Gaussians, the q-Gaussians can also be derived from super statistics. However, the asymmetric q-Gaussian cannot be easily derived from superstatistical distributions, as the processes can be extremely complicated.
Alternatively, q-Gaussians can be derived from either Tsallis entropy or gamma distributions [20,21]. Particularly, gamma distributions allow for the formulation of asymmetric q-Gaussians. An extensive derivation of asymmetric q-Gaussian from gamma distributions can be found in Budini (2015) where he Generally classified the asymmetric q-Gaussians [21]. In this investigation, our interest in the consideration of q-Gaussians stems from the fact that gamma distributions have been observed in satellite rain rate measurements from tropical regions [22].
In this work, we have investigated the complexity of soil temperature and thermal diffusivity dynamics for a better understanding of the statistical mechanics of its variation. An attempt is made to investigate the dynamical processes in the soil temperature variation by considering the Tsallis formalism for q-Gaussian distribution, as we attempt to fit statistical distribution from the measured soil temperature data into Tsallis q-Gaussian. A good number of studies have considered the use of the Tsallis formalism in fitting experimental data into q-Gaussians [23][24][25][26][27][28][29][30]. To achieve this, we evaluate the non-extensive statistical mechanics in soil temperature time series measured from four different locations in Nigeria. The nonextensive statistical mechanics was investigated for soil temperature in this work with the following questions in consideration: i.
Does the Tsallis q-stat parameter vary seasonally for soil temperature? ii. Is there a possible influence of changes in factors such as rainfall on the variations of the q-stat? iii. Is there an influence of vegetational changes such as changes in vegetational cover soil texture and other intrinsic properties?
Our ability to answer these questions will not only shed more light on soil temperature dynamics but also provide a better understanding of the viability of the application of nonextensive statistical mechanics and Tsallis statistics to the characterization and modelling of soil temperature and soil thermal diffusivity in the future.
As mentioned before the q-Gaussian can be derived from Tsallis entropy in which the Tsallis q-index plays a vital role. The Tsallis q-index which was introduced by Tsallis in 1988 provides the generalization of the Boltzmann-Gibbs statistical mechanics. The Tsallis entropy can be derived as follows [20,28,33]: Given that the Shannon entropy from Boltzmann-Gibbs statistical Mechanics is: The non-extensive generalization of the Shannon entropy is therefore given as P i is the probability density function of a discrete random variable x i such that and the q logarithm ln q (x) is given as The parameter q is the Tsallis q-index, which varies for different systems based on the variations in their dynamical characteristics.
Putting expression (4) into (2) we arrive at Tsallis entropy in the form: where p is the probability of each micro-configurations of the series, w is the total number of micro-configurations in the system and q is the Tsallis q-stat (which must be a real number). This Tsallis entropy expression is such that; for q = 1 we obtain a special case of Tsallis entropy which results back into Boltzmann-Gibbs statistical mechanics where q = 1 results in non-extensive statistical mechanics. Considering the continuous case for the random variable x, S q maximizes with the probability as follows: and such that, This gives the Tsallis distribution in the form: Combining Equations (11) and (14) The normalization constant Z q is given by: Considering that the Lagrange multiplier β for (8) is given as: we can show that: and From established theories [14,15], the q-Gaussian for a positive number β can be generalized as follows: where the normalizing constant C q for different values of q expressed as follows: For q < 1, the support of G q (β; x) is compact since the density vanishes for |x| > 1/ (1 − q)β.
In the case of asymmetric q-Gaussians Budini (2015) introduced the asymmetric Poissonian q-Gaussian G p >1 (x q, a, β) relating the function to Tsallis index q, the shape parameters α and α and the asymmetric parameter a [21].
The introduction of the Poissonian distribution in Ref. [21] is anchored on the introduction of two independent gamma functions Y 1 and Y 2 , whose joint probability can be defined as: It is assumed the scale θ the same, while α and α are the shape parameters of the functions Y 1 and Y 2 respectively.
An asymmetric q-Gaussian was derived as a function of functions Y 2 and Y 2 in which the random variables of X as follows: The resulting distribution P(x) for X, which is a function of Y 1 and Y 2 , is given as such that, the normalization constant is C αα given as The maximum value of this P(x) above can be given as (20) In which the asymmetry condition that in the limit βx 1, the power law behaviour emerges as, such that, in the limit βx −1 we have The derived function asymmetric Poissonian function G p >1 (x q, a, β) can be written as follows [21]: The Tsallis q-index is related to the shape parameters α and α as follows: where α = α in this asymmetric case. However, if α = α and a = 0, we have the symmetric case and Equation (7) reduces to The normalization constant in the symmetric case transforms to C α = 2 (2α−1) Γ 2 (α)/Γ(2α). Considering the background on asymmetric q-Gaussian given by Budini [21] and the presence of Gamma distributions in natural systems like precipitation [22], we therefore consider that the presence of asymmetric q-Gaussian in natural system such as soil temperature, which is controlled by various forcing mechanisms is hugely possible. In previous work, we have shown that other systems such as the ionosphere have demonstrated such phenomenon [27].

Data and Methods
In this work, the measured soil temperature data were obtained from the National Space Research and Development Agency (NASRDA) and Centre for Atmospheric Research (CAR) data depositories in Nigeria. The data sets were measured using the sensors deployed in different weather stations at different regions of Nigeria.
The measured temperature data characterizing the dynamics of soil temperature variation tend to possess sporadic variations in time. The time series measurement for soil temperature variation can be seen in Figure 1. The main diagram in Figure 1 is characterised by daily peaks due to diurnal variation, which dominates the intrinsic dynamics. Considering the purpose of our study, the internal variation of the soil temperature dynamics is required. Therefore, to obtain the intrinsic dynamics of the soil temperature we performed a detrending analysis on the soil temperature as to remove the diurnal effects follows [34]: We let the measured time series be x (t i ) with a length v and the fixed-point average temperature values of x f ix(t j ) , where i = 1, 2, 3 . . . . . . v represents the observed time series, and j = 1, 2, 3 . . . ..v, such that the detrended or diurnal variation reduced time is given as The detrending filter produces soil temperature time series x (t i ) as seen in the inner diagram in Figure 1. Further details on this process can be found in Ref. [34]. Similar analyses were carried on other time series to obtain the detrended time series before further analysis.
The datasets were further subjected to stationarity test, as this will ensure the sustainability of the datasets in modelling. The Dicky Fuller methods combined with the Kwiatkowski Philip Schmidt and Shinn (KPSS)-tests were conducted on the detrended soil temperature data sets [35][36][37].
The next step is to analyse our detrended time series data to obtain the probability distribution function (PDF) from which we obtain the values of the stationary state Tsallis q-index (q-stat).
Considering the theory of q-statistics [23][24][25][26], the probability distribution for the energy states {E i } is given by the expression and for the continuous state X in {X}, the probability density function becomes The expression P(x) represents the probability of the dynamics of the microscopic elements x j in states {X}.
The stationary state q-index q stat can be derived based on the above theory of the from the probability distribution functions (29) where PDF (Z) ≡ is the probability distribution with Z as the probability element, coefficients A q , β q denote the normalization constants and q = q stat is the stationary state in this case for which q ≤ 3. Based on the q-Gaussian equation, we can estimate the value of q from time series z i with length i, that is, Z : {z 1 z 2 . . . . . . , z i } estimate the optimum value of q stat from the series z i using the probability density function (PDF) above.
To fit the data to the probability density function, the values X i was divided into windows with width of small value Z with a [Z min , Z max ] for which the frequency of Z; that fall within the Z intervals can be obtained. The normalized histogram for PDF(Z) were obtained, such that, the probability density function PDF{P i (Z i )} represents stationary PDF function.
The resulting fitting would give the corresponding value of Tsallis q-index. To obtain the power relation of this probability density function, we go further obtain the linear relationship between log q [P(Z)] and Z 2 . This procedure was carried out for monthly soil temperature datasets measured from the following locations: Port-Harcourt To evaluate for the degree of asymmetry in the q-Gaussian based on Ref. [21], we obtain the shape parameter ∝ by evaluating the gradient S 1 of linear fit of the normalized probability function H (Hnorm) plotted against the corresponding absolute value of |x| based on (21). Similarly, we obtain the shape parameter ∝ by evaluating the gradient S 2 of linear fit of the normalized probability function H (Hnorm) plotted against the corresponding absolute value of |−x| based in (22). The coefficient of asymmetry a, which can be obtained from the shape parameters ∝ and ∝ are related to the corresponding slopes as follows: Recall that, Therefore, the asymmetric parameter a is related to the gradients as follows: The values of the asymmetric parameters were evaluated from the detrended soil temperature measured for different months of the year, based on this procedure.

Results and Discussion
In this work, the soil temperature time series was analysed to understand the statistical mechanics of the thermal diffusion processes and the energy transfer dynamics in the soil. The test for stationarity shows that the data sets are stationary. The Augmented Dickey Fuller (ADF) test for the data sets gave p values exceedingly less than 0.01, which gives us the confidence to reject the null hypothesis and classify the data sets as stationary and the KPSS test showed unity for all the values. Both cases of test results from the ADF test and the KPSS test show that we can reject the null hypothesis that the soil temperature data is stationary [35][36][37][38]. The stationarity of the data implies that our time series data can preserve the q-Gaussian model.
The monthly Tsallis q-stat was computed from monthly temperature data from different stations in Nigeria. This was carried out by dividing the annual soil temperature data into monthly segments and each of the segments was fitted into q-Gaussian. The values of q-stat obtained showed that the detrended soil temperature data generally satisfy the non-extensive statistical mechanics, as the q-index values were greater than 1 in all cases. The values of q-stat obtained from the q-Gaussian were in the range 1.2 ≤ q ≤ 2.9 (See Tables 1-4). The observed ability of the atmospheric temperature data to fit the q-Gaussian might also be a consequence of the presence of superstatistics. As mentioned earlier, physical systems that exhibit non-extensive statistical mechanics have been linked with superstatistics [13,14,33].   These observed results imply that the intrinsic properties of soil in different locations such as the thermal diffusivity, moisture variation, changes in vegetational cover might be responsible for the changes in the soil temperature dynamics and the variation in q-stat, which are found to be in a particular range for different locations. (See Figures 2-9). We found that the measured soil temperature data sets fit into q-Gaussians with qstat values ranging from 1.2 ≤ q ≤ 2.8. However, it is possible that due to the system reorganization between the rainy season and the dry season, the soil temperature in the savannah region or in the dry season (for southern region) were mostly found to fit into the q-Gaussian model when the data sets were binned into distributions with larger width or class intervals. The distributions were found to assume q-Gaussian more accurately with less noise with larger class intervals during dry season. For instance, we found that for Yola and Makurdi stations, when the soil temperature time series distributions are subdivided into smaller class intervals, the resulting distribution could only fit into the q-Gaussian model from the month of May and June, respectively, and cannot fit into the q-Gaussian model after the rainy season from November (see Tables 3 and 4). The observed possibility of smaller range classes, which allows for a distribution with larger number of bins in the wet season might be because of the presence of rain and high humidity in the raining season, modifying the micro changes in temperature variation. The raining season in Nigeria varies for different locations with longer rainy seasons in the south and shorter rainy seasons in the northern part of Nigeria. A good number of investigations have established that the soil water content is a strong determinant factor in soil temperature regulation [1,10].  The variations in soil texture for different locations might be another factor affecting the nature of the distributions. The soil texture is one of the characteristics of the soil that strongly affects the soil thermal diffusivity [2,7]. The soil structure in the southern part of Guinee savannah where Yola and Makurdi measuring stations are located which is also known as derived Savannah region [39] is mostly characterized by combination of loosely packed aeolian sandy soil and rocky (quartz) sediments, with high silica content [40]. The thermal conduction of soils with high silica can be highly consistent with minimal changes in temperature [7]. The thermal diffusivity of loosely packed soil structure in some areas of the Savannah may also affect the thermal conduction, as conduction is lower in loosely packed soil. As a result, the micro-changes soil temperature occurs in lower magnitude, due to high degree of aeration and conduction gaps resulting from the characteristically large pores found in the soil structure in this region. With less rain in the Savannah region the soil temperature tends to be more consistent. Also, the temperature could be more consistent due to the soil texture. A combination of these characteristics might be responsible for the near zero micro-changes in the soil temperature, hence, the small distributions in class intervals and larger value of q-Gaussian. This characterizes the dynamics of soil temperature in arid/Savannah region and possibly other regions with less rainfall and soil moisture differently, compared to other regions with a higher soil humidity and soil moisture content differently from seasons with rainfall in the same region. This is because the presence of moisture could greatly increase the soil thermal conductivity.  The values of q were suspected to change with the weather conditions, this might be responsible for the variations in q values especially for lower values of q below 1.3 and sometimes close to unity close to wet season in some locations. The lower values of q-stat being recorded in the wet season might be due the situation where the conductivity of the soil varies due to the moisture content. The atmosphere and other environmental systems have been characterized with a mixture of both chaotic and stochastic variations and sometimes with high degrees of freedom. These variations usually result in radical changes in the dynamical component of these environmental systems [41].   Another interesting aspect of our results show that the q-Gaussians fitted from the soil temperature data classified under the extended group of asymmetric q-Gaussians (see , as derived in earlier (see also Ref. [21]). These asymmetric properties could be associated with the seasonal variations in energy transfer between the air and soil interface. The symmetrical changes in the q-Gaussian distributions obtained from the soil temperature data sets can be associated with the inconsistencies in the atmosphere-soil energy transfer interface. However, this asymmetric property might not be totally linked to changes in location as it could be found in the southern region for both Nsukka and Port-Harcourt and toward the northern region at Makurdi (see Figures 2-7). However, the asymmetry in the q-Gaussian found at Yola towards the northern/middle belt region of Nigeria appeared to be lower compared to the q-Gaussians fitted for the other three locations (see Figures 8 and 9). These inconsistencies may range from the inconsistencies due to vegetational cover, soil carbon and other parameters that could affect the thermal diffusivity of the soil from one location to the other. The effect of these varying factors on the soil-atmospheric energy transfer could be responsible for the asymmetric properties of soil temperature q-Gaussian.     Figures 10-13 right panels). The values of asymmetry a tends to vary with season, as the values were lower in month of February, which falls within the dry season in this region, when the q-Gaussians were closer to symmetry. However, the values of asymmetry a were higher in the months of April at Yola and July at Makurdi which fall within early and deep rainy season in this region, respectively.
The asymmetry holds, since all values of a are greater than zero, that is, where a > 0, then all the probability distributions presented here are asymmetrical. However, the degree of asymmetry becomes higher with higher values of |a|. In the cases where the values of |a| → 0 for a given probability distribution, then such a distribution tends towards symmetry.
The asymmetric parameters computed from the time series showed that the distributions obtained from detrended soil temperature were found to be more symmetrical in the dry season month, especially in the months of January and February with higher degrees of asymmetry during the rainy season from April, particularly for Savannah regions, where there is a clear distinction between the rainy season and dry season (see Figures 10-13). These observations demonstrate the possible effect of humidity and rainfall as one of the possible driving forces responsible for changes in symmetrical positions of the q-Gaussian distributions. The shift in the symmetry shows that for different seasons the soil thermal system undergoes self-reorganisation. This changes in the symmetry during the rainy season can be attributed to the direct influence of the rainfall, which is a major atmospheric driver influencing the soil moisture content and hence the temperature regulation [10,41].

Conclusions
In this work, the Tsallis non-extensive statistical mechanics was investigated in soil temperature dynamics to obtain the stationary q-index for soil temperature measured from different regions of Nigeria. This was carried out by fitting the detrended soil temperature data into q-Gaussian. The analysis of soil temperature to obtain the Tsallis q-stat showed that virtually all distributions fit into the q-Gaussian model exhibiting the possible presence of non-extensive statistical mechanics. The ability of the detrended soil temperature time series to fit into Gaussian distribution demonstrates the possibility of representing the soil temperature with non-extensive statistical mechanics. It could also be because of the superstatistical nature of atmospheric systems, which are usually controlled by precipitation. As mentioned earlier, most of the physical systems that are characterized by non-extensive statistical mechanics usually fit into superstatistics. It was observed that the values of Tsallis q-stat varies with changes in the geographical belts of Nigeria as the values of q-stat varies increases towards the southern part of Nigeria. The statistical variations of soil temperature show that the monthly soil temperature distributions are dependent on seasonal variation, particularly considering the varying symmetry of the distributions. The asymmetry in the q-Gaussian, which varies seasonally has been attributed to the seasonal energy variations. The obtained result from this work shows that the Tsallis statistical mechanics can be used to determine and characterise the properties of the soil temperature and thermal diffusivity dynamics. Based on our observations, we infer that there is possible self-reorganization and dynamical forcing due to seasonal dynamical variations in the soil thermal diffusion. However, there is need for further investigations on the non-extensive dynamical properties of soil temperature for clearer understanding the drivers responsible for these variations in the future.