Next Article in Journal
Testing Jump-Diffusion in Epileptic Brain Dynamics: Impact of Daily Rhythms
Previous Article in Journal
A Compositional Model of Consciousness Based on Consciousness-Only
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Stochastic and Self-Organisation Patterns in a 17-Year PM10 Time Series in Athens, Greece

by
Dimitrios Nikolopoulos
1,*,
Aftab Alam
2,
Ermioni Petraki
3,
Michail Papoutsidakis
1,
Panayiotis Yannakopoulos
3 and
Konstantinos P. Moustris
4
1
Department of Industrial Design and Production Engineering, University of West Attica, GR-12244 Aigaleo, Greece
2
Centre for Earthquake Studies, National Centre for Physics, Islamabad 44000, Pakistan
3
Department of Informatics and Computer Engineering, University of West Attica, GR-12233 Aigaleo, Greece
4
Department of Mechanical Engineering, University of West Attica, Petrou Ralli & Thivon 250, GR-12244 Aigaleo, Greece
*
Author to whom correspondence should be addressed.
Submission received: 5 February 2021 / Revised: 24 February 2021 / Accepted: 25 February 2021 / Published: 5 March 2021

Abstract

:
This paper utilises statistical and entropy methods for the investigation of a 17-year PM10 time series recorded from five stations in Athens, Greece, in order to delineate existing stochastic and self-organisation trends. Stochastic patterns are analysed via lumping and sliding, in windows of various lengths. Decreasing trends are found between Windows 1 and 3500–4000, for all stations. Self-organisation is studied through Boltzmann and Tsallis entropy via sliding and symbolic dynamics in selected parts. Several values are below −2 (Boltzmann entropy) and 1.18 (Tsallis entropy) over the Boltzmann constant. A published method is utilised to locate areas for which the PM10 system is out of stochastic behaviour and, simultaneously, exhibits critical self-organised tendencies. Sixty-six two-month windows are found for various dates. From these, nine are common to at least three different stations. Combining previous publications, two areas are non-stochastic and exhibit, simultaneously, fractal, long-memory and self-organisation patterns through a combination of 15 different fractal and SOC analysis techniques. In these areas, block-entropy (range 0.650–2.924) is significantly lower compared to the remaining areas of non-stochastic but self-organisation trends. It is the first time to utilise entropy analysis for PM10 series and, importantly, in combination with results from previously published fractal methods.
Data Set License: license under which the dataset is made available (CC0, CC-BY, CC-BY-SA, CC-BY-NC, etc.)

1. Introduction

Urban air pollution is a related problem for the sustainable environment, human society and economy [1,2]. Imposed by human intervention, urban air pollution is now a sustained part of the urban environment, governed by mechanisms such as emission, dilution and transportation [3]. Significant urban pollutants are the sulphur dioxides (SO2), nitrogen oxides (NOx), ozone (O3) and the respirable particulate matter [4,5,6,7]. The term particulate matter (PM) refers to a mixture of solid particles and liquid droplets hovering in the air close to the ground level (human respiratory height) or even higher. Airborne particulate matter represents a complex mixture of organic and inorganic substances. Some particles, such as dust, dirt, soot or smoke, are large or dark enough to be seen with the naked eye. Other particles are so small that they can only be detected using an electron microscope. Moreover, PM particles are not all of the same shape. Due to these differences and variations between PM particles, the term aerodynamic diameter is introduced. The aerodynamic diameter of a particle is defined as that of a sphere of density 1 g / c m 3 (density of water), which settles in still air at the same velocity as the particle in question [8]. Taking into account this definition, PM with an aerodynamic diameter less than or equal to 10 μ m is referred as PM10, PM with an aerodynamic diameter less than or equal to 5.0 μ m is referred as PM5, PM with an aerodynamic diameter less than or equal to 2.5 μ m is referred as PM2.5, etc.
PM come in many different sizes and shapes and can be made up of hundreds of different chemicals. Some of them are emitted directly from a source, such as construction sites, unpaved roads, fields, smokestacks, fires, etc. Most PM is formed in the atmosphere as a result of complex chemical reactions of substances such as sulphur dioxides (SO2) and nitrogen oxides (NOx). These pollutants are emitted from power plants, industries and vehicles. It is obvious that the size of PM (especially the one of PM10) is crucial to cause adverse health effects in people living in large cities and urban areas in general. This is due to the easiness by which PM10 enter the human body through the respiratory system and the upper respiratory tract. During the last decades, many scientists all over the world have investigated the health impact due to PM exposure [9,10,11,12,13,14,15,16,17,18,19]. Furthermore, the scientific research of the topic focuses on understanding the variations of urban air pollutants with time, by investigating the contribution of the environmental and meteorological parameters to the corresponding variations and modelling and forecasting future alterations [20,21,22,23,24,25,26].
Alternative approaches in PM10 investigation involve concepts that have their origins in the science of complex dynamical systems. Diverging techniques based on chaos have been applied to investigate the inter-state transitions of PM and several other urban air pollution systems. In particular, Kai et al. [27] reported two power-laws in three air pollution time series in Shanghai, China, indicating different self-organised critical (SOC) states. Liu et al. [28] employed detrended fluctuation analysis (DFA) and multifractals to investigate the temporal variations of the concentrations of SO2, NO2 and PM10 in the city of Shanghai, China. They reported that the corresponding air pollution was a self-organised critical (SOC) process. They observed different SOC behaviour within the SO2, NO2 and PM10 concentration time series, that have been attributed to differences in the corresponding power-law relations. Varotsos et al. [29] observed power-law correlations and persistent behaviour of deseasonalised ozone concentration time series with lags between one week and five years. Lee et al. [30] reported scaling invariance and decreasing box dimension function in hourly average of O3 time series concentrations during a whole year by analysing their data-series with monofractals. Windsor and Toumi [31] reported, through the method of R/S analysis, high persistency of hourly ozone concentration time series extending up to 400 days. Similarly, Weng et al. [32] confirmed through (R/S) analysis the nonlinearity, fractality, persistency and long-lasting characteristics of deseasonalised ground level peak ozone concentration time series in southern Taiwan. Xue et al. [33] concluded that R/S analysis, DFA, power-law spectral analysis, linear unbiased estimates, correlation analyses and power spectrum analysis are sophisticated approaches to analyse and describe the distribution, periodicity and trends of air pollutant concentrations. Similarly, Dong et al. [34] employed DFA and multifractal DFA in PM2.5 and PM10 data series and reported the statistical self-affinity, long memory and multifractality of their pollutant time series. Finally, Yuval and Broday [35] employed continuous wavelet transform fractal analysis and found that the predictability of air pollutants is consistent with that of meteorological variables in the short scale.
Up to now, the reporting team has employed different chaotic techniques to investigate and delineate the critical inter-change phases of PM10 pollutants in the Greater Athens Area (GAA) in Greece [36,37]. Detrended fluctuation analysis (DFA), rescaled-range ( R / S ) analysis and fractal-dimension (FD) analysis via three different methods (Higuchi, Katz and Sevcik) have been employed in lengthy 17-year PM10 time series in GAA. Most importantly, a new novel method has been introduced and applied in PM10 time series that is based on the combination of two or more chaotic techniques [36,37]. In [36], a combination of four techniques is reported ( R / S analysis versus FD analysis via the methods of Katz, Higuchi and Sevcik). In [37], a combination of two techniques is employed (DFA versus R / S analysis); however, most significantly, a combination of 13 different chaos analysis techniques is reported. The combination of many different chaos analysis techniques have been employed and checked in critical electromagnetic time series [38] and critical radon in groundwater series [39]. The importance of the above approach is that it employs sliding windows forwarded one sample, covering, in that manner, a fine analysis of the time series.
As follows from the available scientific data of the topic, no investigation has targeted the entropy of urban air-pollution time series. Thus, the present paper focuses on the use of entropy analysis methods for the investigation of lengthy PM10 time series in GAA. Two significant measures of entropy are employed, namely the Boltzmann entropy (BE) and the Tsallis entropy (TE). Both techniques are employed in sliding windows and via symbolic dynamics [40,41,42], in selected parts. It is significant to note here that, since entropy measures the level of organisation of a system, a low order of fractality is related to a low order of entropy and vice versa. Hence, in the consensus of the above-mentioned publications of the reporting team, the use of entropy is a significant alternative viewpoint of the internal non-linear dynamics of the PM10 time series generation system in GAA. Based on a very recent paper of the team [43], which studies this time series with pure statistical methods, the present paper employs statistical techniques to analyse the stochastic tendencies, as well. Similar to the previous studies, sliding windows are employed in the statistical analysis. This yields to analysis of monthly and annually periods of the PM10 system together with other stochastic tendencies. The underlying concept is that the PM10 of GAA is probably both stochastic and non-linear since it has pure stochastic-statistical parts and others with chaos behaviour and long memory. The study combines the stochastic and entropy techniques, in a similar approach as with the fractal and long-memory techniques. Where the PM10 system is out of stochastic behaviour, it is analysed for self-organised tendencies. Several periods are identified and discussed.
The analysis is performed on measurements recorded by the Department of Environmental Protection of the Hellenic Ministry of Environment and Energy at a sampling frequency of 1 m i n 1 . As with the previous two related papers [36,37], this paper employs the same sub-series of the whole dataset, namely the average daily concentrations during 2001–2018 from five different PM10 monitoring stations within the GAA [44]. This sub-series is considered more suitable for long-term analysis because it averages the intra-day variations and minimises, consequently, the overall bias to the estimations. It should also be mentioned that the analysis of PM10 air-pollution systems is a major concern of the scientific community with studies originating from many cities all over the world. Especially for Athens, the capital of Greece, many papers depict that PM10 are among the most harmful pollutants for its inhabitants [7,8,10,11,12].
In the next sections, the research methodology, results and outcomes of this paper are given and discussed. Section 2 presents the experimental methods, namely the area of study and the instrumentation of measurements. Section 3 provides the mathematical background that is further needed in data analysis. More analytically, Section 3.1 provides the statistical methods that are needed to study the stochastic behaviour of the series. Section 3.2 is a brief outline of the methods that are employed in the analysis of self-organisation and fractals. Section 3.3 provides general information on entropy while the specific methods that are employed are given in Section 3.4 and Section 3.5. The results and discussion are presented in Section 4. Due to the diversity and extent of the results, a step-by-step approach is utilised. In this view, the main contributions of this paper are the following: (A) The statistical trends of all the series are delineated. Towards this, the time evolution of the first four statistical moments is calculated through one- and two-month sliding windows and lumping of annual windows. Together, the box-plot-outlier statistical behaviour of the different stations are also reported. Through this procedure and for all station data:(A1) decreasing trends are reported during the last two decades; (A2) noteworthy discrepancies from the Gaussian behaviour are addressed for one- and two-month data analyses; (A3) the annually organised data follow the Gauss distribution; and (A4) deviations of the general statistical patterns are found in association with several outliers. Potential sources are outlined. Importantly. the LYK station is, unexpectedly, found to have similar patterns as ARI and MAR stations. (B) The time evolution of the self-organisation patterns of all stations via Boltzmann and Tsallis entropy is determined through sliding windows. The time-evolution plots are reported, and, from these, the segments with entropy values below critical thresholds are found. The dates of 66 two-month windows are tabulated. (C) The statistical patterns above critical values are found for all stations. Towards this, a previously published statistical methodology, based on the residual concentrations, is used. The corresponding plots are given and the corresponding ASCII data are extracted. (D) The statistical trends of all series are combined with the corresponding self-organisation patterns. Towards this, a previously published two-step meta-analysis methodology is applied to the ASCII data of Points (B) and (C). The dates of nine two-month windows are tabulated with: (D1) statistical patterns above critical values; and (D2) entropy measures below critical values, importantly, from the data of at least three stations simultaneously. (E) The results of Point (D) are compared with the results of the previous publications on the same times series. Of extreme importance is that the dates of three two-month windows are located that are common with the previous publications. The dates of these two-month windows are significant because 15 different techniques provide evidence of critical self-organisation and fractal trends. (F) For the two most significant areas of Point (E), block entropy analyses are reported versus the analyses of the remaining identified parts with SOC patterns. Significantly lower block entropy values are reported for the two areas identified with the combination of 15 techniques, compared to the remaining areas.

2. Experimental Methods

The urban complex of Athens and Piraeus (the port of Athens city) is a large urban area in the Attica Peninsula which includes Athens city, Piraeus city and their suburbs. It extends mainly to the Attica basin and within the administrative boundaries of the five regional units of the Attica Region: Central Sector of Athens, Southern Sector of Athens, Northern Sector of Athens, Western Sector of Athens and Piraeus. In a similar meaning, the term greater Athens area (GAA) is often used. According to the Greek National Census 2011, the GAA has a population of about 2.6 million inhabitants, which corresponds to 69% of the total population of the Attica Region [45]. The GAA basin is surrounded by high mountains. More specifically, on the west side of GAA, the mountain Egaleo (468 m a.s.l.) is located. To the north-northwest of GAA, there is Parnitha Mountain (1,413 m a.s.l.); to the north-north-east, there is Penteli Mountain (1,109 m a.s.l.); and on the east side, there is Ymittos Mountain (1,026 m a.s.l.). The only opening of GAA is located in the south region, the seaside of Saronic Gulf, which is a part of the Aegean Sea. The climate of the GAA is temperate and characterised as a typical Mediterranean climate. In general, sunny days are very common, even in winter. Rains occur mainly from October to April but over the whole year the rainfall amount is very low and does not exceed 400-450 m m . The mean minimum monthly temperature is about 7.0 °C and the mean maximum monthly temperature is about 31.8 °C [46]. Finally, the prevailing winds are south direction winds resulting in the coexistence of the surrounding mountains a major air pollution problem due to the bad air circulation and air pollutants dispersion [46].
The PM10 mean daily concentrations used in this work came from five air pollution monitoring stations (Figure 1 ) located in the GAA. These stations belong to the air pollution monitoring network which operates under the auspices of the Hellenic Ministry of Environment and Energy (HMEE). Pollutants are measured on a continuous basis 24 h a day. The response time of the automatic analysers is of the order of 1 min, i.e. each analyser gives a value approximately every minute. With a microprocessor, located in each automatic station and connected to the automatic analysers, the mean hourly values are calculated every hour. These hourly concentration values are transferred to the server of the service, through a telephone line, and in this way it is possible to continuously monitor the levels of air pollution in the area. Furthermore, the measuring technology of PM10 is based on the method known as absorption of β -radiation. The calibration of PM10 analysers is based on the absorption of β -radiation by a standard sample of known concentration [47]. Figure 1 depicts the map of the GAA and the location of the five air pollution monitoring stations. Table 1 shows the characteristics of the five examined air pollution monitoring stations.

3. Mathematical Methods

3.1. Statistical Methods

At first, statistics are used to explore the general trends of the mean daily PM10 concentrations recorded by the five examined air pollution monitoring stations within the GAA. To investigate the statistical trends of the data, sliding windows of different lengths are employed. In particular, windows of length 32 (approximately one-month data) are moved forwards with a step of 32 samples (one-month analysis lumped one month forwards) and with a step of 1 (one-month analysis slid-glided one day forwards). The latter analysis provides very fine tuning of the PM10 concentrations series, since every 32-length window contains 31 common samples with the previous and the following windows and only one different sample due to the sliding process [37]. The 32-length window with a one-step movement has been applied with success in a recent publication of the team for PM10 air pollution system [37]. As a step further, windows of 64 samples (approximately two-month duration) are utilised as well and moved one sample forwards, in the consensus of another recent publication of the team, also for PM10 time series [37]. In each window, the statistical parameters checked are the average value, standard deviation, kurtosis and skewness of the series, namely the first four statistical moments of the data under investigation. Moreover, windows of 365 samples (annual data) are also employed that are moved forwards 365 samples (lumping process). For this analysis, additionally, the p-value corresponding to the Shapiro–Wilk test is also derived to explore if the corresponding annual data follow the Gaussian distribution. The statistical tests are derived with R package.

3.2. Chaos and Self-Organisation Methods

Several physical systems exhibit complex non-linear behaviour described by chaos. While chaos leads to states that are very sensitive to the initial conditions of the system, there can be critical states (phase space solutions) that are combined in a manner that the system will pass through a series of such states via a complex solution curve that will emerge inevitably. Such states are found in fractal and self-organised systems.
Fractal systems reflect their chaotic behaviour when dilated, translated or rotated in space. Depending on their mathematical properties, these systems are characterised as self-affine or self-similar. Such systems are fractals because the self-affinity and the self-similarity characterise all parts of the system, in the manner that every part is a large- or small-scale imitation or representation of the whole system. Due to this property, the fractal systems can be investigated by analysing their parts. It is notable that scaling and fractal behaviour are associated with the system’s long memory [48,49] and complexity [48], where the complexity of a system indicates the existence of linear mechanisms and orderliness [48]. It is this complex orderliness that emerges as self-organisation of a system. When these self-organised states are connected with the corresponding chaotic critical fractal states, they represent the so-called self-organised critical (SOC) states. SOC and critical fractal states are closely related and also associated with the long memory of the system. For this reason, SOC analysis can reveal fractality and long memory and vice versa [50]. Systems in SOC states can show strong long-lasting associations connecting the past, presence and future of the system. In such critical states, the system will collapse through a series of states where its current state (presence) is affected not only by a previous state (deterministic approach) but also by a series of quite past states (long-lasting interactions of the past), and it also affects (inevitably) its future states in a long-lasting manner.

3.3. Entropy Analysis

Complex nonlinear dynamical systems can be viewed as information capacitors containing discrete series of messages [40,41]. Messaging is achieved by partitioning a continuous time series into a finite number of cells [42]. This procedure is often referred as “symbolic dynamics” [40,42]. Under symbolic dynamics, the actual continuous complex dynamics are viewed under finite precision and can be studied with the concept of entropy [40,50]. The related analysis can be implemented in sequential windows of a finite number of cells [50] (entropy analysis versus time) or in separated series blocks of symbols (block entropy analysis). There are several metrics of entropy that can be effectively utilised, to analyse the complexity of a system [51]. This is because large entropy values imply the existence of many different kinds of patterns and, therefore, high complexity and low organisation. On the other hand, decreased values of entropy (local minima), indicate increased organisation of patterns (local maxima) and, therefore, low order of complexity of a system.
In this paper, two types of entropy measures are employed, namely the Boltzmann entropy and the Tsallis entropy. These two measures are utilised as entropy analysis versus time and as block entropy analysis. They are described in detail below.

3.4. Entropy Analysis Versus Time

To implement entropy analysis versus time, the following steps are followed:
  • The time series is divided into windows of equal size, n.
  • The data of each window are sub-grouped in N equal-spaced bins.
  • In each window, the number of series samples, m j , j [ 1 , N ] , is counted that have values between zones, j 1 and j. The process is iterated for all possible j, i.e., it is implemented for every zone. If N > n , there will be zones that will not contain series values, namely m i = 0 .
  • The probability, P j , in zone j is calculated as P j = m j n . If N > , there will be zones with zero probability, since m j is zero in such cases.
  • In every time series window i , i [ 1 , n ] , the Boltzmann entropy is calculated as
    H ( i ) = k j = 1 . . N P j · l n ( P j )
    where k is the Boltzmann constant, i.e., k = 1.38064852 · 10 23 m 2 k g s 2 K 1 . Boltzmann entropy H ( i ) is a measure of uncertainty and gives the average amount of information needed to predict the actual distribution of measurements in window i. Tsallis entropy is calculated as
    S q ( i ) = k 1 q 1 ( 1 j = 1 . . N p j q )
    where k is the Boltzmann constant and q is a real number measuring the non-extensivity of the system [51]. k is taken as 1.8 according to previous publications [40,41].
The procedure described in this section is a sliding-window technique. Thus, the present version of entropy analysis (versus time) refers to sliding windows. The reader should note that other researchers refer to this method as the gliding window technique [42,52].

3.5. Block Entropy Analysis

In block entropy analysis, the symbolic sub-sequences created by the symbolic organisation of a part of a time series are called words. The words are generated by selecting λ different letters from an alphabet and adopting certain block (word)-size, n, i.e., numbers of sequential letters that will be treated as a whole. Depending on λ and n, the maximum number, N, of different words in the selected alphabet can be determined. For example, for a λ = 2 lettering, a threshold C is adopted. Each value above this threshold is then symbolised by a letter, e.g., 1, while each below by another, e.g., 0 [40,52]. Under this procedure, the initial time series of, for instance, length L = 20 , is transformed into 11001010111000101010 in this λ = 2 lettering. The λ = 2 -letter symbolic time series can be read in discrete sets (lumping technique), of n = 2 blocks, as (11|00|10|10|11|10|00|10|10|10|), where each block is one of the N = 2 2 = 4 different words, i.e., (00, 01, 10,11), in this lettering [40]. Additionally, the above λ = 2 letter sequence may be organised in blocks of n = 3 letters with maximum of N = N 3 , 2 = n λ = 3 2 = 9 different words, i.e., (000, 001, 010, 100, 110, 011, 010, 001, 111). Other sequences of words may be generated as well.
In general [40,52], a L-length time series is transformed, through symbolic dynamics, into a symbolic time series sequence, [ A 1 , A 2 . . . A n . . . A L ] , composed by λ different letters, [ A 1 , A 2 . . . A λ ] , from a λ -length alphabet. Symbolic time series sequences are re-organised in n-sized word blocks composed by letters of the alphabet [ A 1 , A 2 . . . A λ ] . In linguistics, the word size is unconfined and, hence, linguistic words contain some or, potentially, all letters. On the other hand, in symbolic dynamics, the words are of fixed length n, n λ , and are chosen from N max = N n , λ = n λ different fixed-sized words in the [ A 1 , A 2 . . . A λ ] alphabet. In this way, the symbolic time series are re-organised as
. . . A 1 . . . A n B 1 A n + 1 . . . A 2 n B 2 . . . A i n + 1 . . . A ( i + 1 ) n B i + 1 . . .
blocks, where i is the consecutive number of the block, i.e., i = 1 . . . Total number of blocks . In lumping, the n word blocks are sequentially independent sequentially dependent in sliding process. The probability of occurrence of a block, [ A 1 , A 2 . . . A n ] , of size n is calculated as
p ( n ) ( A 1 , A 2 . . . A n ) = Number of occurences of block [ A 1 , A 2 . . . A n ] Total number of blocks
In view of block entropy, Boltzmann block entropy is calculated as [40,52]
H ( n ) = k ( A 1 , A 2 , . . . A n ) p ( n ) ( A 1 , A 2 , . . . A n ) lnp ( n ) ( A 1 , A 2 , . . . A n )
where p ( n ) ( A 1 , A 2 . . . A n ) is the probability of occurrence of a sub-sequence of length n and k is the Boltzmann constant. The entropy H ( n ) is also a measure of uncertainty and gives the average amount of information necessary to predict a certain sub-sequence of length n [40].
Tsallis block entropy is calculated as [40,52]
S q ( n ) = k 1 q 1 ( 1 ( A 1 , A 2 , . . . A n ) [ p ( n ) ( A 1 , A 2 . . . A n ) ] q )
where p ( n ) ( A 1 , A 2 . . . A n ) is the probability of occurrence of block [ A 1 , A 2 . . . A n ] and, as above, k is the Boltzmann constant and q = 1.8 . As already mentioned, a high level of organisation is indicated when low values of Tsallis entropy are produced.
All types of block entropy analysis in this paper are calculated with λ = 2, n = 2–7, namely with two-letter symbolisation and 2–7 fixed-sized words.

4. Results and Discussion

Figure 2 and Figure 3 present the trends of PM10 concentration series according to the sliding window technique (please see Section 3.1). Each sub-figure corresponds to a station of Table 1 characterised by the HMEE as urban-traffic (UT) or suburban-background (SB), accordingly. It is important to mention that UT stations are characterised by a continuously built-up urban area with medium and high buildings in the manner that their pollution levels are mainly determined by the emissions from the nearby sources of traffic (roads, motorways and highways). On the other hand, SB stations are affected indirectly from all kind of sources of pollutants. Thus, they are not influenced by a certain source of pollution, but rather, collectively, by all nearby sources of pollution (industry, traffic, combustion, etc.). Referring to the data in Figure 2 and Figure 3, it should be noted that each data point corresponds to a segmented part of the series with n = 32 samples. Since the measurements are conducted on a daily basis, the n = 32 segment is the nearest power of two window to a monthly analysis. In this sense, Figure 2 and Figure 3 show the time evolution of, approximately, monthly data. The step in these figures is 32; therefore, there is a full coverage of the whole time series via lumping.
All sub-figures of Figure 2 and Figure 3 show a clear decreasing trend during the last two decades (approximately Windows 40–130 depending on the completeness of data). This is more clearly observed in the example case of the THR station shown in Figure 4. The analysis of this figure is in segments of n = 365 samples with a step of 365, namely, the analysis is annual via lumping. The first order approximation of this trend is P M c = 0.76 · y + 1573 ( P M c is the PM10 concentration and y is the year) and quantifies this decreasing trend. The low value of the square of the Spearman’s correlation coefficient of this first order approximation ( r 2 = 0.5) is indicative of the multifaceted process of the PM10 generating system. This is also observed in Figure 5 and Figure 6 where the time evolution of the first four statistical moments of all investigated series is presented. When the data are analysed annually through lumping (Figure 4), the corresponding data follow the normal distribution since the p-value of the Shapiro–Wilk test is very low in all cases (even for the first peak which emerges mainly due to the start of the windowing process). On the other hand, when the data are analysed in segments of n = 64 (approximately two-month segmented data) glided with a step of 1, the observations provide a different viewpoint. Indeed, in Figure 5 and Figure 6, several skewness values are above 1, indicating that the corresponding parts of approximately two-month data (n = 64) do not follow the normal distribution. This is reinforced by the fact that there are several kurtosis values well above the normality kurtosis threshold value of 3. These findings imply, as mentioned above, that the process is multifaceted, and, hence, other mechanisms are present. These mechanisms, on the one hand, set the PM10 series away from normality. On the other hand, as indicated in two recent publications of the team [36,37], chaotic and SOC characteristics may also determine specific parts of the time series. In these time series parts the system is well out of the stochastic behaviour and, as a result, away of normality.
In the above consensus, it should be stressed that the decreasing tendency in the average PM10 values is also observed in Figure 5 and Figure 6 where the sliding window approach of step one is followed. The finer texture is due to the dense sliding window analysis. Despite that, the lumping process also delineates the stochastic trends. As already mentioned in Section 3.3, there are cases where the lumping process is advantageous. As can be observed from these figures, there is significant decreasing trend between Windows 1 and 4000 for the AGP (Figure 5a), THR (Figure 5c), ARI (Figure 6a) and MAR (Figure 6b) series and between Windows 1 and 3500 for the LYK (Figure 5b) series. However, when a whole year of data are included in the analysis, the partial trends are suppressed leading the several variance-inducing trends to diminish. All the above trends are presented integrally in Figure 7. This combination box and whiskers plot shows clearly that the PM10 data have many outliers and wide interquartile ranges (IQR). This is an alternative indication that the PM10 generation system is not only stochastic in nature. The same figure also leads to another observation noted in a very recent publication of the team [43]: ARI, LYK and MAR stations have similar behaviour. This behaviour can also be outlined by observing Figure 5b in comparison to Figure 6 and Figure 2b versus Figure 3. According to these observations, there is indication that the LYK should be characterised as UT and not as a SB.
From the above data, it is evident that the statistical analysis explains partly the variations of the PM10 concentrations in GAA and, hence, may not provide a full outline of the involved physical mechanisms. This is because the related physical phenomena are not stochastic only. As analysed in recent publications of the team [36,37], chaos and long memory are present in various parts of the PM10 series. As described in Section 1 for urban pollution time series and in previous publications of the team for other environmental time series, the presence of chaos and long memory in time series may be associated with self-organised critical (SOC) behaviour and vice versa. In this consensus and following the discussion of Figure 2Figure 6, Figure 8, Figure 9 and Figure 10 present the results of the entropy analysis versus time for windows of length n = 64 and a step of 1. The reader may recall that the use of the sliding window technique of these two figures is identical to the one of Figure 5 and Figure 6. As implemented in the analysis of the recent related publications of the team [36,37], the existing measurements are treated jointly, namely the missing values are neglected, because, as explained in these publications, this data treatment is advantageous for SOC analysis. For other types of analysis, e.g., for smart systems, other approximations are also valid [53].
The majority of the n = 64 sample windows of all examined time series exhibit medium to low variations of the Boltzmann entropy (BE) and Tsallis entropy (TE). Generally, this is an alternative indication that the PM10 system has stochastic behaviour, and, hence, it is in non-critical phases. For this reason, the statistical analysis of the data of these periods is adequate and provides useful outcomes. However, there are periods where BE, TE or both are considerably lower compared to the values of the remaining windows. Similar findings have been reported for the same time series in [36,37] through fractal and long-memory analysis of the data series. It is of importance that the periods of low BE and TE correspond to SOC phases of the system. These SOC phases have been associated with periods of high fractality and long memory [36,37] and are of significance for the environmental PM10 generating system. Indeed, BE is the fundamental measure of randomness [40]. Therefore, if BE is low, then the average amount of information necessary to predict a forthcoming window of length n = 64 is small. On the contrary, in the windows where the system is governed by stochastic processes, BE is elevated, and, hence, there is need for a large amount of information to predict future states. Moreover, there is a great uncertainty in predicting steps in the future provided the history of the present state of the system. There is another significant property of systems in SOC phases that differentiates their behaviour when compared to the stochastic ones. Importantly, the systems in SOC phases or, generally, in self-organisation define their future states not only from their present states but also by their long-term history [36,37]. These systems are non-Markovian in nature, contrary to the stochastic systems where their future is a probabilistic projection of their present states only. The above findings are also verified by the results of the TE. Tsallis entropy drops to lower values in certain windows. In these cases, there is a drop in the number of different kinds of patterns of the system [40]. This implies self-organisation. This finding in conjunction with the results of BE means that the states where BE or TE drop refer to SOC states, characteristic of a chaotic behaviour with high level of prediction [36,37]. It is important to also note that, especially Tsallis entropy, is very sensitive in recognising long-range interactions, long memory and multi-fractals [40]. The current outcomes, hence, imply the presence of fractality and long memory in the environmental PM10 generating system. This verifies the results of the previous studies of the team [36,37,40]. All the above show how important are the dynamics of entropy for the investigation of environmental time series by providing noteworthy measures of their complexity.
As already mentioned, the team has utilised, in four publications [36,37,38,39], combinations of 2–13 different chaotic techniques to analyse time series for hidden fractal and long-memory trends. Two of these publications referred to PM10 time series. As reported, it is not only the identification of critical fractal and long-memory epochs that provide valuable results, but rather the identification of common epochs by different techniques. The issue is also discussed in other publications of the team [50]. In the above consensus, a further step is implemented, namely the combination of the two entropic measures reported in Figure 8, Figure 9 and Figure 10. According to the data presented thus far, it is evident that the PM10 air-pollution system has both stochastic and long-lasting trends. Since the PM10 data exhibit both type of trends, it is proper that the two entropy analysis techniques are combined further with a stochastic type technique.
In the above sense, in [39], chaotic techniques are combined with the so-called relative concentration data, which refer to the residuals of the concentration data. According to the authors, the residual of each PM10 concentration is the term r e s i d u a l ( t ) = c ( t ) c m a ( t ) , where c ( t ) is the daily PM10 concentration and c m a ( t ) is the seven-day moving average. The daily concentrations ( c ( t ) ) of each station are indicated as PM10 in the first sub-plot of each one of three-column multi-plot of Figure 8, Figure 9 and Figure 10. As discussed by Alam et al. [39] (and the references therein), this technique has a long history in successfully delineating anomalous periods of stochastic data (e.g., radon in soil, groundwater and atmosphere). The anomalous periods are usually identified as those with r e s i d u a l ( t ) exceeding the limits of ± σ , ± 2 σ or even ± 3 σ , where σ is the standard deviation of the whole data. Following this viewpoint, the facts presented in the penultimate paragraph of the Introduction and the results presented thus far, it is evident that the PM10 air-pollution environmental system has both stochastic and long-lasting parts, the latter where entropy drops (see also discussion in previous publications of the team [40,41] and the references therein). In the consensus described above, the parameter r e s i d u a l ( t ) is, indeed, a measure of the stochastic tendencies of the PM10 time series. In this viewpoint, the critical limit for r e s i d u a l ( t ) within a certain window to characterise this as a stochastic anomaly could be selected as ± σ , ± 2 σ or even ± 3 σ . In the opinion of the authors, ± σ is adequate for the analysis. This is because the ± σ selection is a fair compromise between the several hidden trends of the data in comparison to the normal distribution trend, i.e. the trend needed so that σ is a meaningful measure of randomness-stochastic behaviour. Following this logic, if a daily PM10 relative concentration ( r e s i d u a l ( t ) ) exceeds ± σ , it is considered as an anomaly from a stochastic tendency ( σ is the standard deviation of the whole data). Accepting the limit for the stochastic trends of the data, if cut-off limits for entropy (BE and TE) are set, as described in the previous publications of the team [36,37,38,39], the combination of the three techniques (one statistical and two entropic) can be implemented.
Figure 11 presents the r e s i d u a l ( t ) as a function of time for all stations in Table 1. It is emphasised again that r e s i d u a l ( t ) is a linear function of the daily PM10 concentrations and that it is a measure of stochastic behaviour of the data. For integrity, both figures also present the ± 1.96 · σ , which is the value corresponding to the 95% confidence interval. The non-stochastic phases of the PM10 air-pollution system are those corresponding to out-of-thresholds values. Clearly, if the PM10 system might exhibit low-entropy trends, it has to be only in the non-stochastic parts, hence where it is out of the stochastic thresholds. Of course, this is only a necessary condition. The sufficient condition for a low entropy-high complexity state (also a low fractal state) is the one where the system is both out of stochastic behaviour and in a low entropy state. Towards this, the daily PM10 concentrations for which r e s i d u a l ( t ) is out of thresholds are extracted in ASCII format as a function of the time moments for which the exceedance is found. A similar procedure is followed for the exceedances of BE and TE of Figure 8, Figure 9 and Figure 10. The corresponding cut-off limits are selected as −2.0 for BE and 1.18 for TE. These limits are arbitrarily set as in the previous related publications [36,37,38,39] in accordance to the ASCII data of the temporal variations of BE and TE. This part of the analysis is described as Step 2 meta-analysis of the data. Thereafter, the out-of-thresholds time series are searched for common dates (Step 3 meta-analysis). Summarising the steps for the analysis are as follows: (a) Calculate TE, BE and r e s i d u a l time series. (b) Set cut-off limits and calculate the out-of-thresholds time series (Step 2 meta-analysis). (c) Locate the common dates for which BE, TE and r e s i d u a l out-of-thresholds time series are commonly addressed for every possible combination (Step 3 meta-analysis). For this study, the possible combinations are BE versus TE and RC, BE versus TE and TE versus RC.
Table 2 presents the common dates for every station of this study in accordance to the cut-off limits mentioned above and the meta-analysis of the ASCII data of the calculations of Figure 8, Figure 9, Figure 10 and Figure 11. As in the previous publications of the team, the presented dates are those that correspond to the start time of the window of analysis, considering that this describes the whole window. It is interesting that, under a similar calculation methodology (Step 2 and 3 meta-analyses) as in the previous related papers of the team [36,37], several common areas are identified. This is because, in this paper, only three techniques are joined together, whereas, in [37], 13 different techniques are joined. It should be emphasised that, as mentioned above, the identification of common dates via different techniques provides far more evidence regarding the justification that, during these dates, the PM10 environmental system is, indeed, in critical states (here SOC states). As analysed above and in the related papers of the reporting team, during these critical states, the PM10 system is in a chaotic and SOC regime that corresponds to an extreme solution path of the space phase volume, where the system behaves in a predictable non-Markovian manner, such that the present state of the system is not only determined by its previous states but also by its long history, and, most importantly, it determines also its future states, in a fractal and SOC behaviour, i.e., its presence interferes with its long-term history and future.
It can be observed that certain dates are commonly found by more than one station. For this reason, Table 3 presents only the dates that are commonly found by at least three stations. Six dates are commonly identified by four stations. For these dates, there is increased evidence that the underlying environmental system that generates the PM10 variations is in SOC states. This is a very significant finding because the dates of Table 3 show clearly that the SOC regime is common to at least three stations. That clearly indicates that are other mechanisms that lead the system to critical predictable states. As reported by Nikolopoulos et al. [37], such mechanisms are:
  • combined, continuously altering, power-law fractal or SOC state interactions of the PM10 air pollution system with several diverse sources of pollutants such as industrial processes, vehicular emissions and energy production from power stations, coupled with complicated physical and chemical processes [27,34];
  • different long-range temporal correlations for the small and large fluctuations of the PM10 time series [34,54], for example, different types of persistency or anti-persistency at the small and large time scales [28];
  • scale-free and self-affine fractal PM10 systems with the power-law scaling and long-term memory [54] recognised as a footprint of SOC behaviour [28];
  • non-linear dynamics governing the PM10 system defined by low-dimensional chaos [55]; and
  • exogenous forces (e.g., meteorological processes and photochemical reactions) can amplify or reduce the linear and non-linear mechanisms of the interactions among the ambient pollutants, which can further complicate the temporal structure of the dynamics of ambient pollutant concentrations, leading to complex behaviour of the air pollution system [56].
Most important, however, are the dates that are presented in Table 4. These are common dates during which the following constraints are simultaneously met for the PM10 system: (a) is out of stochastic behaviour; (b) shows low entropy values both for BE and TE; and (c) presents critical behaviour commonly found with different techniques (Points (a) and (b) from this study). This is an uncommon finding and should be emphasised. Of extreme importance are the first two dates. During these dates, 13 + 2 = 15 different techniques coincide (two from this study, BE and TE). For these two dates, there is very strong evidence that the PM10 system exhibits critical fractal and SOC behaviour. Date 3 of Table 4 is of lesser significance compared to the other two, because there is coincidence only from four techniques (including the two of this study, BE and TE). The reader should emphasise the fact that it is not trivial to find SOC and fractal traces simultaneously. In the related literature, the analysis is usually confined to visual observations and some statistical analysis of the data. Even in the literature that focuses on fractal and SOC trends of environmental data (please see the Introduction and the references cited therein), to the knowledge of the authors, the majority of the related papers report findings with some techniques, but not in combination. The fact that the findings of this study, have coincidences with the outcomes of previous papers of the team justifies the validity of the whole approach, and, most significantly, it is in agreement to the related theory that predicts coincidence between critical fractal and SOC trends. It should be noted that this interesting finding is only addressed in the three dates of Table 4 and not in the others of Table 3. This can be attributed to the multi-parameter nature of the related processes (combined effect of, traffic, industry, wind, temperature gradients, rainfall alterations, particulate dust particles, heating of houses, etc.). The non-linear inter-relations of these related parameters lead the PM10 system out of stochastic behaviour force it to enter into non-linearity and, in some cases, deterministic chaos and SOC states. The reader should also recall that the dates of Table 4 correspond to the beginning dates of segmented-windowed data, therefore, the coincidence is wider in time.
It is evident from the above that the first two dates of Table 4 are the most significant. For this reason, block entropy analysis is further applied in the two-month windows that correspond to these dates. As implied in Section 3.5, it is significant to emphasise in relation that block entropy provides the local SOC trends of the time series. This is because the actual local variations are re-organised as ones and zeros and, thereafter, as words of different letters. As discussed in the references cited in Section 3.5, the technique is very advantageous when comparing time series segments of different hidden-internal patterns. To allow comparisons of the aforementioned two-month segment areas with other parts of the PM10 data, it is crucial to locate proper PM10 parts. The window parts that correspond to the dates of Table 3, but excluding the ones of Table 4, can be located for which the PM10 system: (a) is out of stochastic behaviour; and (b) exhibits some kind of SOC trends. To generate a wide window area, all the windows of Table 3, except those of Table 4, are joined together in an area, hereafter called the "net" area. To discriminate this with the previous areas, the windows-areas of the first two dates of Table 4 are hereafter called "Area 1" and "Area 2". Table 5 presents the results of block entropy analysis utilising the same metrics as in the entropy analysis versus time, i.e., Boltzmann block entropy (BE) and Tsallis block entropy (TE).
All values of Table 5 have a ±0.007 calculation uncertainty. This was estimated by biasing the start of the window of calculation by ±10 samples, where the value of 10 was arbitrarily selected. Observing Table 5, the following information can be supported: (a) BE and TE block entropy values saturate above n = 5 letters. This saturation has been reported by others [51]. (b) The five different stations of the investigation exhibit differentiations in the BE and TE block entropy values; however, the different values are of similar ranges. (c) TE has significantly lower values compared to BE. The most significant observation is the following: (d) BE and TE entropy values of the net-area are significantly higher than the ones of Areas 1 and 2. This implies a fact of extreme importance. The self-organisation of the PM10 system in Areas 1 and 2 is significantly higher than the one of the net area. The reader should recall that the net area also has SOC patterns. Therefore, this identification of higher organisation in Areas 1 and 2 provides very well based evidence that the system in Areas 1 and 2 is, actually, in the most significant SOC states. Since these areas are also identified in the previous publications of the team, there are 15 different methodologies that provide this evidence. It is very important that in Areas 1 and 2 the system is not only in chaos and SOC states but in highly predictable states. Therefore, it is in solutions where its presence is defined not only by its past but also by its long-term history, and, significantly for non-Markovian systems, it defines its long-term future. This important finding should be emphasised and is a significant outcome of this research.
Despite the increased evidence found for the two dates of Areas 1 and 2, all the dates of Table 3 are important, as well as the results of other dates reported in the previous publications of the team [36,37]. The above facts are indicative of the advancement of the new methodology described in the previous publications [36,37,38,39]. It is significant that the technique has been checked in other environmental systems (electromagnetic radiation disturbances and radon in soil variations). Three related publications of the team clearly show that the PM10 environmental system can be described by similar techniques. It is also important to emphasise that this is not the first research to address fractal and SOC trends in PM10 air-pollution variations. As presented in the several papers cited in the Section 1, many researchers have reported such trends. However, this is the first time to combine so many different techniques and locate the most important time series parts found from very lengthy time series. Of course the fact that only two areas are commonly found does not lower the validity of the reported data, because there are many more areas with SOC trends given in Table 3. Moreover, the paper analyses the stochastic trends of the data in a systematic manner. Importantly, many areas of the whole time series exhibit stochastic behaviour, and this is the trivial hypothesis when working with environmental air-pollution time series. However, it is also studied very systematically that there are also non-stochastic fractal and SOC phases which cannot be neglected and should be taken into account in the related studies. Finalising, it is an issue of further research to investigate more data and other air-pollution time series.

5. Conclusions

In the following, the most significant findings of the present paper are summarised:
  • Statistical and entropy analysis methods are employed for the study of a 17-year PM10 time series recorded from five stations in Athens, Greece. The purpose is to investigate further the stochastic trends of the series and explore if non-stochastic periods of low entropy values exist.
  • The stochastic trends are analysed in monthly, two-month and annual windows via lumping and sliding windows. Decreasing time-evolution trends were found: (a) between Windows 40 and 130 for n = 32 and a step of 32 (lumping process); and (b) between Windows 1 and 4000 for the AGP, THR, ARI and MAR series and between Windows 1 and 3500 for the LYK series, all for n = 64 and a step of 1 (sliding window process). For Case (b), periods with high variance, skewness above 1 and kurtosis well above 3 are found. Deviations from the Gaussian distribution are addressed in various parts and, therefore, non-statistical behaviour. However, when the data are analysed for n = 365 and a step of 365 (annual analysis through lumping), they follow the normal distribution. Several outliers are also found. ARI, LYK and MAR stations have similar IQR behaviour, indicating that LYK should be characterised as UT. The discrepancies from the stochastic behaviour are attributed to the multiple facets involved in physical procedures. The related mechanisms are discussed.
  • Self organisation is investigated via Boltzmann and Tsallis entropy. Sliding windows of n = 64 and a step of 1 are employed and symbolic dynamics in selected parts. The majority of the n = 64 windows of all examined time series exhibit medium to low variations of both entropy values. Several low entropy parts are identified with Boltzmann entropy over Boltzmann constant below -2.0 entropy and Tsallis entropy over the Boltzmann constant below 1.18. The implications of the findings are discussed.
  • A previously published combination two-step method is utilised to locate areas for which the PM10 system is out of stochastic behaviour and, simultaneously, exhibits critical self-organised patterns.According to previous publications, the non-stochastic periods are taken those with r e s i d u a l ( t ) out of ± σ . These periods are located and extracted to ASCII data. The two-step methodology is applied to those periods and the critical entropy periods of Conclusion (b) above. Sixty-six different periods of two-month duration are found in the time series for various dates. From these, nine periods are common to at least three different stations. This is very significant because it is associated with internal dynamics of the GAA basin.
  • Searching for enhanced evidence of SOC and fractal trends, as well as long memory, the findings of this study are further compared to previous publications for the same series. Two areas are found for which the series is non-stochastic and exhibit, simultaneously, fractal, long-memory and self-organisation patterns through a combination of 15 different fractal and SOC analysis techniques.
  • In the two most significant areas, block entropy analysis is further applied. For comparison, block entropy is also utilised in the remaining identified parts with SOC patterns. Block entropy values are in the range 0.650–2.924 for words of 2–7 letters. Block entropy saturates above five letters. Significantly lower block entropy values are found for the two areas identified with the combination of 15 techniques, compared to the remaining areas.
  • This is the first time to utilise entropy analysis for PM10 series and, importantly, in combination with results from previously published fractal methods. This combination approach is expected to be applied in the future to other environmental series, e.g., in environmental ozone and PM2.5 series, as well as in more seismic related series, such as radon in groundwater and soil and M H z k H z environmental electromagnetic disturbances. Despite the limited number of critical windows found, the enhanced evidence, on the one hand, and the increased outcomes, on the other hand, which are expected to be received in more dense time series, make the present approach a significant study tool. Future expectations are the extent of the methodology with multifractal techniques to account for the multifractality found in nature, as well as combining the present techniques with more advanced statistical approaches.

Author Contributions

Conceptualization, D.N.; Data curation, D.N. and E.P.; Formal analysis, D.N., M.P. and P.Y.; Investigation, D.N., A.A., E.P. and K.P.M.; Methodology, D.N. and E.P.; Project administration, K.P.M.; Resources, D.N., A.A. and K.P.M.; Software, D.N.; Supervision, D.N.; Visualization, D.N. and A.A.; Writing-original draft, D.N., A.A., E.P., M.P., P.Y. and K.P.M.; Writing-review & editing, D.N., E.P., M.P., P.Y. and K.P.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Chan, C.K.; Yao, X. Air pollution in mega cities in China. Atmos. Environ. 2008, 42, 1–42. [Google Scholar] [CrossRef]
  2. Fang, M.; Chan, C.K.; Yao, X.H. Managing air quality in a rapidly developing nation: China. Atmos. Environ. 2009, 43, 79–86. [Google Scholar] [CrossRef]
  3. Liu, Y.; Wang, T. Worsening urban ozone pollution in China from 2013 to 2017 – Part 1: The complex and varying roles of meteorology. Atmos. Chem. Phys. 2020, 20, 6305–6321. [Google Scholar] [CrossRef]
  4. Hand, J.L.; Prenni, A.; Copeland, S.; Schichtel, B.A.; Malm, W.C. Thirty years of the Clean Air Act Amendments: Impacts on haze in remote regions of the United States (1990–2018). Environment. Atmos. Environ. 2020, 243, 117865. [Google Scholar] [CrossRef]
  5. Thomas, D. Geochemical precursors to seismic activity. Pure Appl. Geophys. 1988, 126, 241–252. [Google Scholar] [CrossRef]
  6. Cíardenas Rodríguez, M.; Dupont-Courtade, L.; Oueslati, W. Air pollution and urban structure linkages: Evidence from European cities. Renew. Sust. Energ. Rev. 2016, 53, 1–9. [Google Scholar] [CrossRef]
  7. Lee, C. Impacts of urban form on air quality: Emissions on the road and concentrations in the US metropolitan areas. J. Environ. Manage. 2019, 246, 196–202. [Google Scholar] [CrossRef]
  8. Hinds, W.C. Aerosol Technology: Properties, Behavior, and Measurement of Airborne Particles, 2nd ed.; Wiley-VCH: Hoboken, NJ, USA, 1999; p. 53. [Google Scholar]
  9. Lee, Y.; Lee, J.; Kim, H.; Ha, E. Effects of PM10 on mortality in pure COPD and asthma-COPD overlap: Difference in exposure duration, gender, and smoking status. Sci. Rep. 2020, 10, 2402. [Google Scholar] [CrossRef] [Green Version]
  10. Hwang, J.; Bae, H.; Choi, S.; Yi, H.; Ko, B.; Kim, N. Impact of air pollution on breast cancer incidence and mortality: A nationwide analysis in South Korea. Sci. Rep. 2020, 5392, 770–808. [Google Scholar] [CrossRef]
  11. Zanobetti, A.; Schwartz, J. The effect of fine and coarse particulate air pollution on mortality: A national analysis. Health Perspect. 2009, 6, 898–903. [Google Scholar] [CrossRef] [Green Version]
  12. World Health Organisation. Air Quality Guidelines, 2nd ed.; WHO Regional Office for Europe: Copenhagen, Denmark, 2000; pp. 1–40. [Google Scholar]
  13. Lee, S.; Jang, M.; Oh, S.; Lee, J.; Suh, M.; Park, M. Associations between Particulate Matter and Otitis Media in Children: A Meta-Analysis. Int. J. Environ. Res. Public Health 2020, 17, 4604. [Google Scholar] [CrossRef] [PubMed]
  14. Zhang, F.; Luo, L.; Wang, Z.; Zhang, W.; Li, C.; Qiu, Z.; Huang, D. Estimation of the Effects of Air Pollution on Hospitalization Expenditures for Asthma. Int. J. Health Serv. 2020, 50, 100–109. [Google Scholar] [CrossRef]
  15. Tsoli, S.; Ploubidis, G.B.; Kalantzi, O.I. Particulate air pollution and birth weight: A systematic literature review. Atmos, Pollut. Res. 2019, 10, 1084–1122. [Google Scholar] [CrossRef]
  16. Braithwaite, I.; Zhang, S.; Kirkbride, J.; Osborn, D.; Hayes, J. Air Pollution (Particulate Matter) Exposure and Associations with Depression, Anxiety, Bipolar, Psychosis and Suicide Risk: A Systematic Review and Meta-Analysis. Environ. Health Perspect. 2019, 127, 12602. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Thompson, J.E. Airborne Particulate Matter. J. Occup. Environ. Med. 2018, 60, 392–423. [Google Scholar] [CrossRef] [PubMed]
  18. Nasser, Z.; Salameh, P.; Nasser, W.; Abou Abbas, L.; Elias, E.; Leveque, A. Outdoor particulate matter (PM) and associated cardiovascular diseases in the Middle East. Int. J. Occup. Med. Environ. Health 2015, 28, 641–661. [Google Scholar] [CrossRef]
  19. Goldberg, M. Particulate Air Pollution and Daily Mortality: Who Is at Risk? J. Aerosol. Med. 1996, 9, 43–53. [Google Scholar] [CrossRef]
  20. Ibarra-Berastegi, G.; Elias, A.; Barona, A.; Saenz, J.; Ezcurra, A.; de Argando na, J. From diagnosis to prognosis for forecasting air pollution using neural networks: Air pollution monitoring in Bilbao. Environ. Model. Softw. 2008, 23, 622–637. [Google Scholar] [CrossRef]
  21. Kourentzes, N.; Barrow, D.; Crone, S. Neural network ensemble operators for time series forecasting. Expert Syst. Appl. 2014, 41, 4235–4244. [Google Scholar] [CrossRef] [Green Version]
  22. Pérez, P.; Trier, A.; Reyes, J. Prediction of PM2.5 concentrations several hours in advance using neural networks in Santiago, Chile. Atmos. Environ. 2000, 34, 1189–1196. [Google Scholar] [CrossRef]
  23. Moustris, K.P.; Ziomas, I.C.; Paliatsos, A.G. 3-Day-Ahead Forecasting of Regional Pollution Index for the pollutants NO2, CO, SO2 and O3 using Artificial Neural Networks in Athens, Greece. Water Air Soil. Pollut. 2010, 224, 29–43. [Google Scholar] [CrossRef]
  24. Moustris, K.P.; Larisssi, I.K.; Nastos, P.T.; Koukouletsos, K.V.; Paliatsos, A.G. Development and Application of Artificial Neural Network Modeling in Forecasting PM10 Levels in a Mediterranean City. Water Air Soil Pollut. 2013, 224, 1–11. [Google Scholar] [CrossRef]
  25. Cacciola, M.; Pellicanó, D.; Megali, G.; Lay-Ekuakille, A.; Versaci, M.; Morabito, F.C. Aspects about Air Pollution prediction on urban environment. In Proceedings of the 4th Imeko TC19 Symposium on Environmental Instrumentation and Measurements Protecting Environment, Climate Changes and Pollution Control, Lecce, Italy, 3–4 June 2013; pp. 15–20. [Google Scholar]
  26. Catalano, M.; Galatioto, F.; Bell, M.; Namdeo, A.; Bergantino, A. Improving the prediction of air pollution peak episodes generated by urban transport networks. Environ. Sci. Policy 2016, 60, 69–83. [Google Scholar] [CrossRef] [Green Version]
  27. Kai, S.; Chun-qiong, L.; Nan-shan, A.; Xiao-hong, Z. Using three methods to investigate time-scaling properties in air pollution indexes time series. Nonlinear Anal. Real World Appl. 2008, 9, 693–707. [Google Scholar] [CrossRef]
  28. Liu, Z.; Wang, L.; Zhu, H. A time-scaling property of air pollution indices: A case study of Shanghai. China Atmos. Pollut. Res. 2015, 6, 886–892. [Google Scholar] [CrossRef]
  29. Varotsos, C.; Ondov, J.; Efstathiou, M. Scaling properties of air pollution in Athens, Greece and Baltimore. Md. Atmos. Environ. 2005, 39, 4041–4047. [Google Scholar] [CrossRef]
  30. Lee, C.K.; Juang, L.C.; Wang, C.C.; Liao, Y.Y.; Yu, C.C.; Liu, Y.C.; Ho, D.S. Scaling characteristics in ozone concentration time series (OCTS). Chemosphere 2006, 62, 934–946. [Google Scholar] [CrossRef] [PubMed]
  31. Windsor, H.L.; Toumi, R. Scaling and persistence of UK pollution. Atmos. Environ. 2001, 35, 4545–4556. [Google Scholar] [CrossRef]
  32. Weng, Y.C.; Chang, N.B.; Lee, T.Y. Nonlinear time series analysis of ground-level ozone dynamics in Southern Taiwan. J. Environ. Manag. 2008, 87, 405–414. [Google Scholar] [CrossRef]
  33. Xue, Y.; Pan, W.; Lu, W.z.; He, H.D. Multifractal nature of particulate matters (PMs) in Hong Kong urban air. Sci. Total Environ. 2015, 532, 744–751. [Google Scholar] [CrossRef]
  34. Dong, Q.; Wang, Y.; Li, P. Multifractal behavior of an air pollutant time series and the relevance to the predictability. Environ. Pollut 2017, 222, 444–457. [Google Scholar] [CrossRef]
  35. Broday, Y.D.M. Studying the time scale dependence of environmental variables predictability using fractal analysis. Environ. Sci. Technol. 2010, 44, 4629–4634. [Google Scholar] [CrossRef]
  36. Nikolopoulos, D.; Moustris, K.; Petraki, E.; Koulougliotis, D.; Cantzos, D. Fractal and Long-Memory Traces in PM10 Time Series in Athens, Greece. Environments 2019, 6, 29. [Google Scholar] [CrossRef] [Green Version]
  37. Nikolopoulos, D.; Moustris, K.; Petraki, E.; Cantzos, D. Long-memory traces in PM10 time series in Athens, Greece: Investigation through DFA and R/S analysis. Meteorol. Atmospheric Phys. 2020. [Google Scholar] [CrossRef]
  38. Nikolopoulos, D.; Petraki, E.; Yannakopoulos, P.; Priniotakis, G.; Voyiatzis, I.; Cantzos, D. Long-Lasting Patterns in 3 kHz Electromagnetic Time Series after the ML = 6.6 Earthquake of 2018-10-25 near Zakynthos, Greece. Geosciences 2020, 10, 235. [Google Scholar] [CrossRef]
  39. Alam, A.; Wang, N.; Zhao, G.; Mehmood, T.; Nikolopoulos, D. Long-lasting patterns of radon in groundwater at Panzhihua, China: Results from DFA, fractal dimensions and residual radon concentration. Geochem. J. 2019, 53, 341–358. [Google Scholar] [CrossRef]
  40. Petraki, E.; Nikolopoulos, D.; Fotopoulos, A.; Panagiotaras, D.; Nomicos, C.; Yannakopoulos, P.; Kottou, S.; Zisos, A.; Louizi, A.; Stonham, J. Long-range memory patterns in variations of environmental radon in soil. Anal. Methods 2013, 5, 4010–4020. [Google Scholar] [CrossRef]
  41. Nikolopoulos, D.; Petraki, E.; Vogiannis, E.; Chaldeos, Y.; Giannakopoulos, P.; Kottou, S.; Nomicos, C.; Stonham, J. Traces of self-organisation and long-range memory in variations of environmental radon in soil: Comparative results from monitoring in Lesvos Island and Ileia (Greece). J. Radioana.l Nucl. Chem. 2014, 299, 203–219. [Google Scholar] [CrossRef]
  42. Karamanos, K.; Nicolis, G. Symbolic Dynamics and Entropy Analysis of Feigenbaum Limit Sets. Chaos Soliton. Frac. 1999, 10, 1135–1150. [Google Scholar] [CrossRef]
  43. Moustris, K.; Petraki, E.; Ntourou, K.; Priniotakis, G.; Nikolopoulos, D. Spatiotemporal Evaluation of PM10 Concentrations within the Greater Athens Area, Greece. Trends, Variability and Analysis of a 19 Years Data Series. Environments 2020, 53, 85. [Google Scholar] [CrossRef]
  44. PM10 Full Dataset between 2001–2018 in Greece. Available online: https://ypen.gov.gr/perivallon/poiotita-tis-atmosfairas/dedomena-metriseon-atmosfairikis-rypansis/ (accessed on 1 February 2021).
  45. Authority, H.S. Population-Housing Census; Hellenic Statistical Authority: Piraeus, Greece, 2011. [Google Scholar]
  46. Climatic Data for Selected Stations in Greece. Available online: http://www.hnms.gr/emy/en/climatology/climatology_city?perifereia=Attiki&poli=Athens_Hellinikon/ (accessed on 1 February 2021).
  47. Amaral, S.; Carvalho, J.; Costa, M.; Pinheiro, C. An Overview of Particulate Matter Measurement Instruments. Atmosphere 2015, 6, 1327–1345. [Google Scholar] [CrossRef] [Green Version]
  48. Morales, I.O.; Landa, O.; Fossion, R.; Frank, A. Scale invariance, self-similarity and critical behaviour in classical and quantum system. J. Phys. Conf. Ser. 2012, 380, 12020. [Google Scholar] [CrossRef] [Green Version]
  49. Musa, M.; Ibrahim, K. Existence of long memory in ozone time series. Sains Malays. 2012, 41, 1367–1376. [Google Scholar]
  50. Petraki, E. Electromagnetic Radiation and Radon-222 Gas Emissions as Precursors of Seismic Activity. Ph.D. Thesis, Department of Electronic and Computer Engineering, Brunel University, London, UK, 2016. [Google Scholar]
  51. Kalimeri, M.; Papadimitriou, C.; Balasis, G.; Eftaxias, K. Dynamical complexity detection in pre-seismic emissions using non-additive Tsallis entropy. Physica A 2008, 387, 1161–1172. [Google Scholar] [CrossRef]
  52. Karamanos, K. Entropy analysis of substitutive sequences revisited. J Phys. A Math. Gen. 2001, 34, 9231–9241. [Google Scholar] [CrossRef]
  53. Izonin, I.; Tkachenko, R.; Verhun, V.; Zub, K. An approach towards missing data management using improved GRNN-SGTM ensemble method. Eng. Sci. Technol. Int. J. 2020. [Google Scholar] [CrossRef]
  54. Lu, W.; Xue, Y.; di He, H. Detrended fluctuation analysis of particle number concentrations on roadsides in Hong Kong. Build. Environ. 2014, 82, 580–S587. [Google Scholar] [CrossRef]
  55. Khokhlov, V.N.; Glushkov, A.V.; Loboda, N.S.; Bunyakova, Y.Y. Short-range forecast of atmospheric pollutants using non-linear prediction method. Atmos. Environ. 2008, 42, 7284–7292. [Google Scholar] [CrossRef]
  56. Yu, H.L.; Lin, Y.C.; Sivakumar, B.; Kuo, Y. A study of the temporal dynamics of ambient particulate matter using stochastic and chaotic techniques. Atmos. Environ. 2015, 69, 37–45. [Google Scholar] [CrossRef]
Figure 1. The map of the GAA and the location of the five air pollution monitoring stations.
Figure 1. The map of the GAA and the location of the five air pollution monitoring stations.
Entropy 23 00307 g001
Figure 2. Average PM10 concentrations (μg/m3) for windows of length n = 32 and a step of 32 (one-month data lumped one month forward). GAA, period 2001–2010: (a) AGP series (b) LYK; and (c) THR series. PM10 concentration measurements in μg/m3. The diagrams are produced with R package.
Figure 2. Average PM10 concentrations (μg/m3) for windows of length n = 32 and a step of 32 (one-month data lumped one month forward). GAA, period 2001–2010: (a) AGP series (b) LYK; and (c) THR series. PM10 concentration measurements in μg/m3. The diagrams are produced with R package.
Entropy 23 00307 g002
Figure 3. Average PM10 concentrations ( μ g / m 3 ) for windows of length n = 32 and a step of 32 (one month data lumped one month forward). GAA, period 2001–2010: (a) ARI series; and (b) MAR series. PM10 concentration measurements in μ g / m 3 . The diagrams are produced with R package.
Figure 3. Average PM10 concentrations ( μ g / m 3 ) for windows of length n = 32 and a step of 32 (one month data lumped one month forward). GAA, period 2001–2010: (a) ARI series; and (b) MAR series. PM10 concentration measurements in μ g / m 3 . The diagrams are produced with R package.
Entropy 23 00307 g003
Figure 4. Average PM10 concentrations ( μ g / m 3 ) and p-value of the Shapiro–Wilk test for windows of length n = 365 and a step of 365 (one year data separated over the years). GAA, period 2001–2010, THR series. The diagrams are produced with R package.
Figure 4. Average PM10 concentrations ( μ g / m 3 ) and p-value of the Shapiro–Wilk test for windows of length n = 365 and a step of 365 (one year data separated over the years). GAA, period 2001–2010, THR series. The diagrams are produced with R package.
Entropy 23 00307 g004
Figure 5. Time evolution of the first four statistical moments of segmented data of the time series for windows of length n = 64 with a step of 1. GAA, period 2001–2010: (a) AGP series; (b) LYK; and (c) THR series. Measurements in μ g / m 3 . The diagrams are produced with R package.
Figure 5. Time evolution of the first four statistical moments of segmented data of the time series for windows of length n = 64 with a step of 1. GAA, period 2001–2010: (a) AGP series; (b) LYK; and (c) THR series. Measurements in μ g / m 3 . The diagrams are produced with R package.
Entropy 23 00307 g005
Figure 6. Time evolution of the first four statistical moments of segmented data of the time series for windows of length n = 64 with a step of 1. GAA, period 2001–2010: (a) ARI series; and (b) MAR series. Measurements in μ g / m 3 . The diagrams are produced with R package
Figure 6. Time evolution of the first four statistical moments of segmented data of the time series for windows of length n = 64 with a step of 1. GAA, period 2001–2010: (a) ARI series; and (b) MAR series. Measurements in μ g / m 3 . The diagrams are produced with R package
Entropy 23 00307 g006
Figure 7. Box and whiskers plot of the full dataset PM10 concentrations in GAA in the period 2001–2010.
Figure 7. Box and whiskers plot of the full dataset PM10 concentrations in GAA in the period 2001–2010.
Entropy 23 00307 g007
Figure 8. Time evolution of Boltzmann Entropy (BE) (light blue, top) and Tsallis Entropy (TE) (green, middle) in parallel to the original time series (dark blue, bottom): (a) AGP series; and (b) THR series. Analysis in segmented data of the time series for windows of length n = 64 with a step of 1. GAA, period 2001–2018. k is the Boltzmann’s constant. Average PM10 concentrations are in μ g / m 3 . The diagrams are produced with Octave Forge package.
Figure 8. Time evolution of Boltzmann Entropy (BE) (light blue, top) and Tsallis Entropy (TE) (green, middle) in parallel to the original time series (dark blue, bottom): (a) AGP series; and (b) THR series. Analysis in segmented data of the time series for windows of length n = 64 with a step of 1. GAA, period 2001–2018. k is the Boltzmann’s constant. Average PM10 concentrations are in μ g / m 3 . The diagrams are produced with Octave Forge package.
Entropy 23 00307 g008
Figure 9. Time evolution of Boltzmann Entropy (BE) (light blue, top) and Tsallis Entropy (TE) (green, middle) in parallel to the original time series (dark blue, bottom) for the LYK series. Analysis in segmented data of the time series for windows of length n = 64 with a step of 1. GAA, period 2001–2018. k is the Boltzmann’s constant. Average PM10 concentrations are in μ g / m 3 . The diagrams are produced with Octave Forge package.
Figure 9. Time evolution of Boltzmann Entropy (BE) (light blue, top) and Tsallis Entropy (TE) (green, middle) in parallel to the original time series (dark blue, bottom) for the LYK series. Analysis in segmented data of the time series for windows of length n = 64 with a step of 1. GAA, period 2001–2018. k is the Boltzmann’s constant. Average PM10 concentrations are in μ g / m 3 . The diagrams are produced with Octave Forge package.
Entropy 23 00307 g009
Figure 10. Time evolution of Boltzmann Entropy (BE) (light blue, top) and Tsallis Entropy (TE) (green, middle) in parallel to the original time series (dark blue, bottom): (a) ARI series; and (b) MAR series. Analysis in segmented data of the time series for windows of length n = 64 with a step of 1. GAA, period 2001–2018. k is the Boltzmann’s constant.Average PM10 concentrations are in μ g / m 3 . The diagrams are produced with Octave Forge package.
Figure 10. Time evolution of Boltzmann Entropy (BE) (light blue, top) and Tsallis Entropy (TE) (green, middle) in parallel to the original time series (dark blue, bottom): (a) ARI series; and (b) MAR series. Analysis in segmented data of the time series for windows of length n = 64 with a step of 1. GAA, period 2001–2018. k is the Boltzmann’s constant.Average PM10 concentrations are in μ g / m 3 . The diagrams are produced with Octave Forge package.
Entropy 23 00307 g010
Figure 11. Time evolution of r e s i d u a l ( t ) of the PM10 concentration data as a function of time. GAA, period 2001–2018. k is the Boltzmann’s constant. Average PM10 concentrations are in μ g / m 3 . ± · σ and ± 1.96 · σ cut-off limits are presented as well.
Figure 11. Time evolution of r e s i d u a l ( t ) of the PM10 concentration data as a function of time. GAA, period 2001–2018. k is the Boltzmann’s constant. Average PM10 concentrations are in μ g / m 3 . ± · σ and ± 1.96 · σ cut-off limits are presented as well.
Entropy 23 00307 g011
Table 1. The characteristics of the examined air pollution monitoring stations within the GAA. Abbr., abbreviation; Alt., altitude above sea level (a.s.I); D.C., data completeness.
Table 1. The characteristics of the examined air pollution monitoring stations within the GAA. Abbr., abbreviation; Alt., altitude above sea level (a.s.I); D.C., data completeness.
Monitoring StationAbbr.LongitudeLatitudeAlt. (m)CharacterisationD.C.
AristotelousARI23°43’39"37°59’16"75Urban-Traffic85.8%
LykovrissiLYK23°47’19"38°04’04"234Suburban-Background89.2%
MaroussiMAR23°47’14”38°01’51”170Urban-Traffic82.5%
Agia ParaskeviAGP23°49’09”37°59’42”290Suburban-Background88.7%
ThrakomakedonesTHR23°45’29’38°08’36”550Suburban-Background77.2%
Table 2. Results of the common areas of all investigated stations exhibiting simultaneously: (a) r e s i d u a l ( t ) out of ± σ ; (b) BE value below 2.0; and (c) TE below 1.18. i/i represents the index value of the common date of each station.
Table 2. Results of the common areas of all investigated stations exhibiting simultaneously: (a) r e s i d u a l ( t ) out of ± σ ; (b) BE value below 2.0; and (c) TE below 1.18. i/i represents the index value of the common date of each station.
Monitoring Stationi/iDate
AGP12005/12/3
22005/12/5
32007/3/24
42007/3/25
52007/7/28
52007/11/6
62009/4/4
72009/4/6
72010/6/7
82010/12/28
92011/6/9
102013/5/31
112013/6/1
122013/6/2
132013/6/3
142013/6/4
152014/6/26
162014/6/27
172014/10/19
172015/1/8
182015/2/6
192015/2/2
202016/5/13
202016/7/17
ARI12007/3/25
22009/4/4
32009/4/6
42010/3/5
52010/3/6
62010/3/7
72010/3/8
82014/6/26
92014/6/27
102015/2/6
112015/2/2
112016/5/13
122016/7/17
LYK52007/7/28
12009/4/4
22010/2/22
32010/2/23
42010/2/25
52010/6/7
62014/6/26
72014/6/27
82015/1/8
92016/5/23
102016/5/24
112016/5/25
122016/5/27
MAR12007/3/24
22007/3/25
32007/6/14
52007/7/28
32009/4/4
42009/4/6
52010/2/22
62010/2/23
72010/6/7
82014/6/26
92014/6/27
102015/1/8
112015/1/22
122015/2/6
132015/2/2
142015/2/7
152016/7/22
THR12009/4/4
22009/4/6
32010/12/28
42011/6/9
52015/2/6
62015/2/2
72016/7/17
Table 3. Results of the areas of Table 2 that are commonly found by at least three stations. All investigated dates exhibit simultaneously: (a) r e s i d u a l ( t ) out of ± σ ; (b) BE value below 2.0; and (c) TE below 1.18. i/i represents the index value of the common dates.
Table 3. Results of the areas of Table 2 that are commonly found by at least three stations. All investigated dates exhibit simultaneously: (a) r e s i d u a l ( t ) out of ± σ ; (b) BE value below 2.0; and (c) TE below 1.18. i/i represents the index value of the common dates.
i/iDateMonitoring Stations
1.2007/3/25AGP, ARI, MAR
2007/7/28AGP, LYK, MAR
2.2009/4/4AGP, ARI, LYK, THR
3.2009/4/6AGP, ARI, MAR, THR
4.2010/6/7AGP, LYK, MAR
5.2014/6/26AGP, ARI, LYK, MAR
6.2014/6/27AGP, ARI, LYK, MAR
2015/1/8AGP, LYK, MAR
7.2015/2/6AGP, ARI, MAR, THR
8.2015/2/2AGP, MAR, THR
9.2016/7/7AGP, ARI, THR
Table 4. Results of the areas of Table 3 that are commonly addressed with the previous related papers of the team with other techniques.
Table 4. Results of the areas of Table 3 that are commonly addressed with the previous related papers of the team with other techniques.
i/iDateMonitoring StationsTechniquesPublication
1.2007/7/28AGP, LYK, MAR13 different fractal techniquesNikolopoulos et al. [36]
2.2010/6/7AGP, LYK, MAR13 different fractal techniquesNikolopoulos et al. [36]
3.2015/1/8AGP, LYK, MARDFA & RS analysisNikolopoulos et al. [37]
Table 5. Results of the block entropy analysis of selected areas of the PM10 time series. Clarifications are written in text.
Table 5. Results of the block entropy analysis of selected areas of the PM10 time series. Clarifications are written in text.
Area 1 Area 2 Area Net
StationLettersBETEBETEBETE
AGP21.0950.6821.1080.6501.2310.743
31.6790.8961.1430.6991.8470.921
41.8590.9271.3980.7012.3641.005
52.0950.9971.8340.9392.8471.064
62.0950.9971.8350.9382.8431.069
72.0970.9961.8320.9432.8461.062
ARI21.1960.7311.2060.7351.3280.792
31.6650.8911.4130.8011.8280.910
42.1191.0091.6780.9002.4421.039
52.1090.9971.9450.9862.8951.090
62.1110.9981.9460.9872.8971.091
72.1120.9981.9470.9872.8981.091
LYK21.0450.6511.3050.7891.2670.757
31.7080.8921.7510.9071.8870.936
41.8640.9572.0260.9752.3961.017
51.9810.9762.1381.0152.9241.094
61.9820.9772.1391.0152.2951.096
71.9810.9782.1411.0162.2971.097
MAR21.2270.7461.3230.8021.2310.740
31.8210.9431.7750.9241.7810.892
42.2111.0242.2111.0242.2810.986
52.3691.0572.2220.9932.5341.035
62.3711.0592.2210.9952.5321.037
72.3691.0622.2230.9972.5361.037
THR21.2080.7031.1570.6951.2430.754
31.8780.9411.4780.8091.8290.919
41.8410.9461.9250.9742.3971.031
51.9810.9762.0250.9932.8491.068
61.9830.9772.0270.9912.8561.069
71.9840.9752.0290.9942.8421.067
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Nikolopoulos, D.; Alam, A.; Petraki, E.; Papoutsidakis, M.; Yannakopoulos, P.; Moustris, K.P. Stochastic and Self-Organisation Patterns in a 17-Year PM10 Time Series in Athens, Greece. Entropy 2021, 23, 307. https://0-doi-org.brum.beds.ac.uk/10.3390/e23030307

AMA Style

Nikolopoulos D, Alam A, Petraki E, Papoutsidakis M, Yannakopoulos P, Moustris KP. Stochastic and Self-Organisation Patterns in a 17-Year PM10 Time Series in Athens, Greece. Entropy. 2021; 23(3):307. https://0-doi-org.brum.beds.ac.uk/10.3390/e23030307

Chicago/Turabian Style

Nikolopoulos, Dimitrios, Aftab Alam, Ermioni Petraki, Michail Papoutsidakis, Panayiotis Yannakopoulos, and Konstantinos P. Moustris. 2021. "Stochastic and Self-Organisation Patterns in a 17-Year PM10 Time Series in Athens, Greece" Entropy 23, no. 3: 307. https://0-doi-org.brum.beds.ac.uk/10.3390/e23030307

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop