Spectral Observations of PM10 Fluctuations in the Hilbert Space

*Thomas Plocoste and Rudy Calif*

### **Abstract**

During the last 20 years, many megacities have experienced air pollution leading to negative impacts on human health. In the Caribbean region, air quality is widely affected by African dust which causes several diseases, particularly, respiratory diseases. This is why it is crucial to improve the understanding of PM10 fluctuations in order to elaborate strategies and construct tools to predict dust events. A first step consists to characterize the dynamical properties of PM10 fluctuations, for instance, to highlight possible scaling in PM10 density power spectrum. For that, the scale-invariant properties of PM10 daily time series during 6 years are investigated through the theoretical Hilbert frame. Thereafter, the Hilbert spectrum in time-frequency domain is considered. The choice of theoretical frame must be relevant. A comparative analysis is also provided between the results achieved in the Hilbert and Fourier spaces.

**Keywords:** PM10 data, empirical mode decomposition, Hilbert spectral analysis, time-frequency representation, Fourier space

#### **1. Introduction**

Generally, the concentration of air pollutants varies and is impacted by the local pollutant emission levels and meteorological and topographical conditions [1, 2]. Particulate matter (PM) is a complex mixture of elemental and organic carbon, ammonium, nitrates, sulfates, mineral dust, trace elements, and water [3]. PM with an aerodynamic diameter of <10 μm, i.e., PM10, are well known for their impact on human health [4]. Many studies have highlighted that exposure to PM increases the number of hospital admissions for cardiovascular disease, acute bronchitis, asthma attacks, respiratory disease, and congestive heart failure [5–8]. In the Caribbean area, one of the main emitters of PM10 is from large-scale sources, i.e., African dust [9]. Knowledge of the dynamics of PM10 process is crucial to elaborate strategies and construct tools to predict dust events. The time-frequency distribution of a signal provides information about how the spectral content of a signal evolves with time, thus providing an ideal tool to dissect, analyze, and interpret nonstationary signals [10]. Contrary to classical methods, the need of a timefrequency representation (TFR) is stemmed from the inadequacy of either time domain or frequency domain analysis to fully describe the nature of nonstationary signals [10]. In literature, there are numerous methods to obtain energy density as a function of time and frequency simultaneously as the short-time Fourier transform (STFT), Hilbert-Huang transform (HHT), and wavelet transform (WT) [10–12].

In this study, the scaling properties of PM10 data are firstly analyzed, and then the TFR is investigated. In order to highlight the performance of the Hilbert space, an analysis of PM10 data was also performed in the Fourier space.

computing of power spectral density (PSD). The PSD separates and measures the amount of variability occurring in different frequency bands. In this study, PSD are

In order to investigate the scaling properties of PM10 data, classically the discrete Fourier transform of the times series considered is computed. The expression of Fourier transform *X*(*f*) for a process *x*(*t*) is recalled here. An *N* point-long interval is used to construct the value at frequency domain point *f*, *Xf* [17]:

*X f* ð Þ¼ <sup>ð</sup>þ*<sup>T</sup>*

�*T*

Consequently the power spectral density *E*(*f*) is estimated by computing the

To determine the scale invariance of a given time series in a joint amplitudefrequency space, the Hilbert-Huang transform [19, 20] is performed. HHT can be summarized in two steps: (i) empirical mode decomposition (EMD) and (ii) Hilbert spectral analysis (HSA). Empirical mode decomposition is a powerful tool to separate a nonlinear and nonstationary time series into a sum of intrinsic mode functions (IMF) without a priori basis as required by traditional Fourier-based method [19–21]. An IMF must satisfy the following two conditions: (i) the difference between the number of local extrema and the number of zero-crossings must be zero or one, and (ii) the local maxima and the envelope defined by the local minima are close to zero. Therefore, the original signal *x*(*t*) is decomposed into a sum of *n*1

*x t*ðÞ¼ <sup>X</sup>*<sup>n</sup>*�<sup>1</sup>

*<sup>C</sup>*~*m*ðÞ¼ *<sup>t</sup>*

1 *π P*

*m*¼1

To obtain a physically significant IMF, this selection process must be stopped by a certain criterion. For more details, EMD decomposition is widely described in the

To characterize the time-frequency energy distribution from the original signal *x* (*t*), HSA is applied on each obtained IMF component *Cm*(*t*) to extract the instantaneous amplitude and frequency [19, 24]. The Hilbert transform is defined by:

> ðþ<sup>∞</sup> �∞

with *P* the Cauchy principal value [24, 25]. We can specify an analytical signal *z*

*Cm t* <sup>0</sup> ð Þ

*x t*ð Þ*e*

ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi <sup>R</sup>*e*<sup>2</sup>ð Þþ *X f* ð Þ <sup>I</sup>*m*<sup>2</sup>ð Þ *X f* ð Þ <sup>q</sup>

�2*πift*dt (1)

*E f* ð Þ¼ j j *X f* ð Þ <sup>2</sup> (3)

*Cm*ð Þþ*t rn*ð Þ*t* (4)

*<sup>t</sup>* � *<sup>t</sup>*<sup>0</sup> dt<sup>0</sup> (5)

(2)

estimated through the Fourier and Hilbert spaces.

*DOI: http://dx.doi.org/10.5772/intechopen.88279*

*Spectral Observations of PM10 Fluctuations in the Hilbert Space*

Thus, the analytical expression of *X*(*f*) is [18]

∣*X f* ð Þ∣ ¼

*3.1.1 Fourier analysis*

following expression:

*3.1.2 Hilbert analysis*

literature [19–23].

**93**

IMF modes with the residual *rn*(*t*):

for each IMF mode *Cm*(*t*) with

This chapter is organized as follows. Section 2 presents PM10 data analyzed in this study. Section 3 describes the methods applied in order to investigate PM10 dynamics. Section 4 comments on the results obtained and then discusses them.

### **2. Experimental data**

Guadeloupe archipelago is a French West Indies island located in the middle of the Caribbean basin, i.e., 16.25°N latitude and 61.58°W longitude, which experiences a tropical and humid climate [13, 14]. The time series analyzed here belong to Guadeloupe air quality network which is managed by the Gwad'Air agency (http:// www.gwadair.fr/). PM10 concentrations are measured at Pointe-à-Pitre (16.2422°N 61.5414°W) using the Thermo Scientific tapered element oscillating microbalance (TEOM) models 1400ab and 1400-FDMS. Hourly PM10 concentrations were sampled during the period from 1 January 2005 to 31 December 2010. We processed these data into daily average concentrations. In total, there are 2150 daily averaged data points available continuously for 6 years. **Figure 1** displays PM10 daily signal illustrating huge fluctuations and thus indicating a strong variability. These strong oscillations observed in the middle of each year are attributed to PM10 related to dust outbreaks coming from the African coast from May to September [9]. For the rest of the year, PM10 is mainly generated by anthropogenic pollution [15].

**Figure 1.**

*Illustration of PM10 daily average concentrations between 2005 and 2010, highlighting intermittent burst events with huge fluctuations.*

#### **3. Methods**

#### **3.1 Scaling analysis (1D representation)**

The description of natural phenomena by the study of statistical scale laws is not recent [16]. Self-similarity of complex systems has been widely observed in nature and is the simplest form of scale invariance. A scale invariance can be detected by

computing of power spectral density (PSD). The PSD separates and measures the amount of variability occurring in different frequency bands. In this study, PSD are estimated through the Fourier and Hilbert spaces.

## *3.1.1 Fourier analysis*

In this study, the scaling properties of PM10 data are firstly analyzed, and then the TFR is investigated. In order to highlight the performance of the Hilbert space,

This chapter is organized as follows. Section 2 presents PM10 data analyzed in this study. Section 3 describes the methods applied in order to investigate PM10 dynamics. Section 4 comments on the results obtained and then discusses them.

Guadeloupe archipelago is a French West Indies island located in the middle of the Caribbean basin, i.e., 16.25°N latitude and 61.58°W longitude, which experiences a tropical and humid climate [13, 14]. The time series analyzed here belong to Guadeloupe air quality network which is managed by the Gwad'Air agency (http:// www.gwadair.fr/). PM10 concentrations are measured at Pointe-à-Pitre (16.2422°N 61.5414°W) using the Thermo Scientific tapered element oscillating microbalance (TEOM) models 1400ab and 1400-FDMS. Hourly PM10 concentrations were sampled during the period from 1 January 2005 to 31 December 2010. We processed these data into daily average concentrations. In total, there are 2150 daily averaged data points available continuously for 6 years. **Figure 1** displays PM10 daily signal illustrating huge fluctuations and thus indicating a strong variability. These strong oscillations observed in the middle of each year are attributed to PM10 related to dust outbreaks coming from the African coast from May to September [9]. For the rest of the year, PM10 is mainly generated by anthropogenic pollution [15].

The description of natural phenomena by the study of statistical scale laws is not recent [16]. Self-similarity of complex systems has been widely observed in nature and is the simplest form of scale invariance. A scale invariance can be detected by

**2008**

*Illustration of PM10 daily average concentrations between 2005 and 2010, highlighting intermittent burst*

**2009**

**2010**

**2011**

an analysis of PM10 data was also performed in the Fourier space.

**2. Experimental data**

*Functional Calculus*

**3. Methods**

*events with huge fluctuations.*

**Figure 1.**

**92**

**3.1 Scaling analysis (1D representation)**

**2005**

**PM10 (µg/m3**

**)**

**2006**

**2007**

In order to investigate the scaling properties of PM10 data, classically the discrete Fourier transform of the times series considered is computed. The expression of Fourier transform *X*(*f*) for a process *x*(*t*) is recalled here. An *N* point-long interval is used to construct the value at frequency domain point *f*, *Xf* [17]:

$$X(f) = \int\_{-T}^{+T} x(t)e^{-2\pi f t}dt\tag{1}$$

Thus, the analytical expression of *X*(*f*) is [18]

$$|X(f)| = \sqrt{\mathcal{Re}^2(X(f)) + \mathcal{Im}^2(X(f))}\tag{2}$$

Consequently the power spectral density *E*(*f*) is estimated by computing the following expression:

$$E(f) = |\mathcal{X}(f)|^2 \tag{3}$$

#### *3.1.2 Hilbert analysis*

To determine the scale invariance of a given time series in a joint amplitudefrequency space, the Hilbert-Huang transform [19, 20] is performed. HHT can be summarized in two steps: (i) empirical mode decomposition (EMD) and (ii) Hilbert spectral analysis (HSA). Empirical mode decomposition is a powerful tool to separate a nonlinear and nonstationary time series into a sum of intrinsic mode functions (IMF) without a priori basis as required by traditional Fourier-based method [19–21]. An IMF must satisfy the following two conditions: (i) the difference between the number of local extrema and the number of zero-crossings must be zero or one, and (ii) the local maxima and the envelope defined by the local minima are close to zero. Therefore, the original signal *x*(*t*) is decomposed into a sum of *n*1 IMF modes with the residual *rn*(*t*):

$$\varkappa(t) = \sum\_{m=1}^{n-1} \mathcal{C}\_m(t) + r\_n(t) \tag{4}$$

To obtain a physically significant IMF, this selection process must be stopped by a certain criterion. For more details, EMD decomposition is widely described in the literature [19–23].

To characterize the time-frequency energy distribution from the original signal *x* (*t*), HSA is applied on each obtained IMF component *Cm*(*t*) to extract the instantaneous amplitude and frequency [19, 24]. The Hilbert transform is defined by:

$$
\tilde{C}\_m(t) = \frac{1}{\pi} P \int\_{-\infty}^{+\infty} \frac{C\_m(t')}{t - t'} \,\mathrm{d}t' \tag{5}
$$

with *P* the Cauchy principal value [24, 25]. We can specify an analytical signal *z* for each IMF mode *Cm*(*t*) with

$$Z\_m(t) = \mathbf{C}\_m(t) + j\tilde{\mathbf{C}}\_m(t) = A\_m(t)e^{j\phi\_m(t)}\tag{6}$$

where *Am*ðÞ¼ *t* ∣*Zm*ð Þ*t* ∣ ¼ ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi *Cm*ð Þ*<sup>t</sup>* <sup>2</sup> <sup>þ</sup> *<sup>C</sup>*~*m*ð Þ*<sup>t</sup>* <sup>2</sup> q describes an amplitude and *<sup>φ</sup>m*ðÞ¼ *<sup>t</sup>* argð Þ¼ *<sup>z</sup>* arctan *<sup>C</sup>*~*m*ð Þ*<sup>t</sup> Cm*ð Þ*<sup>t</sup>* � h represents the phase function of IMF modes. Consequently, the instantaneous frequency *ωm*ð Þ*t* is defined from the phase *φm*ð Þ*t* by

$$
\rho\_m(t) = \frac{1}{2\pi} \frac{d\rho\_m(t)}{dt} \tag{7}
$$

SPEC*x*ð Þ¼ *<sup>t</sup>; <sup>f</sup>* j j *Sx*ð Þ *<sup>t</sup>; <sup>f</sup>* <sup>2</sup> (12)

j j *x t*ð Þ <sup>2</sup>

ðþ<sup>∞</sup> �∞

*dτ* is the STFT of *x*(*t*), *w*(*τ*) is a window

dt (13)

<sup>⩽</sup>*f*⩽4*:*<sup>57</sup> � <sup>10</sup>�<sup>5</sup> Hz which

where *Sx*ð Þ¼ *<sup>t</sup>; <sup>f</sup>* <sup>Ð</sup> <sup>þ</sup><sup>∞</sup>

at point (*t*,*f*) [12]:

resolution.

**4. Results**

*E f* ð Þ� *f*

**Figure 2.**

**95**

*significant only in the Hilbert space.*

**4.1 Scaling properties**

*3.2.2 Hilbert spectrum*

�<sup>∞</sup> *<sup>x</sup>*ð Þ*<sup>τ</sup> <sup>w</sup>*ð Þ *<sup>τ</sup>* � *<sup>t</sup> <sup>e</sup>*�*jf <sup>τ</sup>*

ðþ<sup>∞</sup> �∞ ðþ<sup>∞</sup> �∞

*Spectral Observations of PM10 Fluctuations in the Hilbert Space*

*DOI: http://dx.doi.org/10.5772/intechopen.88279*

(e.g., Hanning, rectangular, Hamming), *t* is time, and *f* is frequency.

As depicted in Eq. (13), SPEC roughly describes the energy density of the signal

SPEC*x*ð Þ *t; f* dtdf ¼

The SPEC has been applied successfully in various research fields [12, 35–37]. The main advantages of SPEC are an easily understanding interpretation, and it allows a fast computation. However, the main drawback of SPEC is the same as that of the STFT [12]. Indeed, there is a trade-off between time and frequency

The Hilbert spectrum (HS) is a joint time-frequency representation introduced by [19]. It is important to notice that the two important tools (i.e., EMD and HS) for exploratory analysis of the data are provided by HSA method. This approach was applied successfully in various research fields as fault diagnosis for rolling bearing [11], turbulence [38], environment [34, 39], and geophysics [40], to cite a few.

In order to identify the presence of scaling in PM10 time series, the PSD is estimated in the Hilbert and Fourier spaces. **Figure 2** depicts the power spectral density provided by the Hilbert transform and the Fourier transform. On this figure, we try to detect a power law behavior of the form *<sup>h</sup>*ð Þ� *<sup>ω</sup> <sup>ω</sup>*�*β<sup>h</sup>* and

corresponds to time scales 6*:*1 hours⩽*T*⩽55*:*4 days, a power law behavior is clearly

**PM10**

**10-9 10-8 10-7 10-6 10-5 f, (Hz)**

*The spectrum of PM10 time series in the Hilbert space and the Fourier space. A power law behavior is*

and Fourier spaces. On the frequency range 2*:*<sup>09</sup> � <sup>10</sup>�<sup>7</sup>

**10<sup>2</sup>**

**10<sup>4</sup>**

**S(f), H( )** **10<sup>6</sup>**

�*β<sup>f</sup>* where *β<sup>h</sup>* and *β<sup>f</sup>* are, respectively, the spectral exponents in the Hilbert

**H( ) Hilbert spectrum slope H( ) -1.02 0.10 S(f) Fourier spectrum**

Thus, the original signal *x*(*t*) can be expressed as

$$\mathbf{x}(t) = \mathcal{R}\mathbf{e} \sum\_{m=1}^{N} \mathcal{A}\_{m}(t) e^{j\varphi\_{m}(t)} = \mathcal{R}\mathbf{e} \sum\_{m=1}^{N} \mathcal{A}\_{m}(t) e^{j\int\_{-\infty}^{t} a\_{m}(t)dt} \tag{8}$$

where R*e* is a part real [19, 20, 26].

Due to the simultaneous representation of frequency modulation and amplitude modulation, the HHT can be considered as a generalization of the Fourier transform [19, 20]. The energy in a time-frequency space is designated as the Hilbert spectrum with *<sup>H</sup>*ð Þ¼ *<sup>ω</sup>; <sup>t</sup>* <sup>A</sup><sup>2</sup> ð Þ *ω; t* . The Hilbert marginal spectrum *h*(*ω*) is defined by

$$h(\boldsymbol{\alpha}) = \frac{1}{T} \int\_0^T H(\boldsymbol{\alpha}, t) \mathbf{d}t \tag{9}$$

where *T* is the total data length. The Hilbert spectrum *H*(*ω*,*t*) gives a measure of amplitude from each frequency and time, while the marginal spectrum *h*(*ω*) gives a measure of the total amplitude from each frequency [27]. As a result, the marginal spectrum can be compared to the Fourier spectrum [19, 20].

In conclusion, for a scale-invariant process, the Fourier *E*(*f*) and the Hilbert *h*(*ω*) spectral densities obtained follow a power law over a range of frequencies:

$$E(f) \sim f^{-\beta\_f} \tag{10}$$

$$h(w) \sim w^{-\beta\_h} \tag{11}$$

where *f* and *ω* are the frequencies and *β<sup>f</sup>* and *β<sup>h</sup>* are the spectral exponents, respectively, in the Fourier and Hilbert spaces. It reveals the scale-free memory effect as a power law dependence of the frequency distribution. Consequently, *β<sup>f</sup>* and *β<sup>h</sup>* contain information about the degree of stationarity of the studied parameter [16, 28, 29]:


#### **3.2 Time-frequency representation (2D representation)**

#### *3.2.1 Spectrogram*

The spectrogram (SPEC) of a signal *x*(*t*) is defined as the squared magnitudes of the STFT as shown in Eq. (12) [12]:

*Spectral Observations of PM10 Fluctuations in the Hilbert Space DOI: http://dx.doi.org/10.5772/intechopen.88279*

$$\text{SPEC}\_{\mathbf{x}}(t, f) = |\mathcal{S}\_{\mathbf{x}}(t, f)|^2 \tag{12}$$

where *Sx*ð Þ¼ *<sup>t</sup>; <sup>f</sup>* <sup>Ð</sup> <sup>þ</sup><sup>∞</sup> �<sup>∞</sup> *<sup>x</sup>*ð Þ*<sup>τ</sup> <sup>w</sup>*ð Þ *<sup>τ</sup>* � *<sup>t</sup> <sup>e</sup>*�*jf <sup>τ</sup> dτ* is the STFT of *x*(*t*), *w*(*τ*) is a window (e.g., Hanning, rectangular, Hamming), *t* is time, and *f* is frequency.

As depicted in Eq. (13), SPEC roughly describes the energy density of the signal at point (*t*,*f*) [12]:

$$\int\_{-\infty}^{+\infty} \int\_{-\infty}^{+\infty} \mathbf{SPEC\_x}(t, f) \mathbf{dt} \mathbf{df} = \int\_{-\infty}^{+\infty} |x(t)|^2 \mathbf{dt} \tag{13}$$

The SPEC has been applied successfully in various research fields [12, 35–37]. The main advantages of SPEC are an easily understanding interpretation, and it allows a fast computation. However, the main drawback of SPEC is the same as that of the STFT [12]. Indeed, there is a trade-off between time and frequency resolution.

#### *3.2.2 Hilbert spectrum*

*Zm*ðÞ¼ *<sup>t</sup> Cm*ðÞþ*<sup>t</sup> jC*~*m*ðÞ¼ *<sup>t</sup> Am*ð Þ*<sup>t</sup> <sup>e</sup>*

ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi *Cm*ð Þ*<sup>t</sup>* <sup>2</sup> <sup>þ</sup> *<sup>C</sup>*~*m*ð Þ*<sup>t</sup>* <sup>2</sup>

quently, the instantaneous frequency *ωm*ð Þ*t* is defined from the phase *φm*ð Þ*t* by

1 2*π*

*<sup>j</sup>φm*ð Þ*<sup>t</sup>* <sup>¼</sup> <sup>R</sup><sup>e</sup>

Due to the simultaneous representation of frequency modulation and amplitude modulation, the HHT can be considered as a generalization of the Fourier transform [19, 20]. The energy in a time-frequency space is designated as the Hilbert spectrum

where *T* is the total data length. The Hilbert spectrum *H*(*ω*,*t*) gives a measure of amplitude from each frequency and time, while the marginal spectrum *h*(*ω*) gives a measure of the total amplitude from each frequency [27]. As a result, the marginal

In conclusion, for a scale-invariant process, the Fourier *E*(*f*) and the Hilbert *h*(*ω*)

ð Þ *ω; t* . The Hilbert marginal spectrum *h*(*ω*) is defined by

*dφm*ð Þ*t*

X *N*

A*m*ð Þ*t e j* Ð*t*

*m*¼1

*ωm*ðÞ¼ *t*

A*m*ð Þ*t e*

*<sup>h</sup>*ð Þ¼ *<sup>ω</sup>* <sup>1</sup> *T* ð*T* 0

spectral densities obtained follow a power law over a range of frequencies:

*E f* ð Þ� *f*

where *f* and *ω* are the frequencies and *β<sup>f</sup>* and *β<sup>h</sup>* are the spectral exponents, respectively, in the Fourier and Hilbert spaces. It reveals the scale-free memory effect as a power law dependence of the frequency distribution. Consequently, *β<sup>f</sup>* and *β<sup>h</sup>* contain information about the degree of stationarity of the studied parameter

• If 1< *β<sup>f</sup>* or *β<sup>h</sup>* <3, the process is nonstationary with increments stationary.

• Spectral analysis has been widely applied in various research fields [30–34].

The spectrogram (SPEC) of a signal *x*(*t*) is defined as the squared magnitudes of

spectrum can be compared to the Fourier spectrum [19, 20].

• If *β<sup>f</sup>* or *β<sup>h</sup>* <1, the process is stationary.

• If *β<sup>f</sup>* or *β<sup>h</sup>* >1, the process is nonstationary.

**3.2 Time-frequency representation (2D representation)**

q

*Cm*ð Þ*<sup>t</sup>* � h

Thus, the original signal *x*(*t*) can be expressed as

X *N*

*m*¼1

*x t*ðÞ¼ Re

where R*e* is a part real [19, 20, 26].

where *Am*ðÞ¼ *t* ∣*Zm*ð Þ*t* ∣ ¼

*<sup>φ</sup>m*ðÞ¼ *<sup>t</sup>* argð Þ¼ *<sup>z</sup>* arctan *<sup>C</sup>*~*m*ð Þ*<sup>t</sup>*

*Functional Calculus*

with *<sup>H</sup>*ð Þ¼ *<sup>ω</sup>; <sup>t</sup>* <sup>A</sup><sup>2</sup>

[16, 28, 29]:

*3.2.1 Spectrogram*

**94**

the STFT as shown in Eq. (12) [12]:

*<sup>j</sup>φm*ð Þ*<sup>t</sup>* (6)

�<sup>∞</sup> *<sup>ω</sup>m*ð Þ*<sup>t</sup>* dt (8)

describes an amplitude and

dt (7)

*H*ð Þ *ω; t* dt (9)

�*β<sup>f</sup>* (10)

*<sup>h</sup>*ð Þ� *<sup>ω</sup> <sup>ω</sup>*�*β<sup>h</sup>* (11)

represents the phase function of IMF modes. Conse-

The Hilbert spectrum (HS) is a joint time-frequency representation introduced by [19]. It is important to notice that the two important tools (i.e., EMD and HS) for exploratory analysis of the data are provided by HSA method. This approach was applied successfully in various research fields as fault diagnosis for rolling bearing [11], turbulence [38], environment [34, 39], and geophysics [40], to cite a few.

#### **4. Results**

#### **4.1 Scaling properties**

In order to identify the presence of scaling in PM10 time series, the PSD is estimated in the Hilbert and Fourier spaces. **Figure 2** depicts the power spectral density provided by the Hilbert transform and the Fourier transform. On this figure, we try to detect a power law behavior of the form *<sup>h</sup>*ð Þ� *<sup>ω</sup> <sup>ω</sup>*�*β<sup>h</sup>* and *E f* ð Þ� *f* �*β<sup>f</sup>* where *β<sup>h</sup>* and *β<sup>f</sup>* are, respectively, the spectral exponents in the Hilbert and Fourier spaces. On the frequency range 2*:*<sup>09</sup> � <sup>10</sup>�<sup>7</sup> <sup>⩽</sup>*f*⩽4*:*<sup>57</sup> � <sup>10</sup>�<sup>5</sup> Hz which corresponds to time scales 6*:*1 hours⩽*T*⩽55*:*4 days, a power law behavior is clearly

#### **Figure 2.**

*The spectrum of PM10 time series in the Hilbert space and the Fourier space. A power law behavior is significant only in the Hilbert space.*

weak energy (in blue). This highlights the energy activity related to PM10 concentrations during the study period. Such an approach gives the possibility of tracking the evolution of PM10 data spectral content in time, which is typically represented by variations of the amplitudes and frequencies of the components from which the

*Spectral Observations of PM10 Fluctuations in the Hilbert Space*

*DOI: http://dx.doi.org/10.5772/intechopen.88279*

On **Figure 4**, strong energies are observed throughout the years with slight fluctuations on the frequency range 0⩽*f*⩽<sup>1</sup> <sup>10</sup><sup>6</sup> Hz. For *<sup>f</sup>*>4 <sup>10</sup><sup>6</sup> Hz, strong energies are also noticed in the middle of each year and at the beginning of 2010 with more fluctuations. In **Figure 5**, energy distributions are more localized. On the frequency range 0⩽*f*⩽<sup>1</sup> <sup>10</sup><sup>6</sup> Hz, we can observe the influence of small-scale event on energy behavior. As noticed, this energy may be weak or null. As an example, the impact of a general strike in early 2009 that paralyzed Guadeloupean archipelago at least 2 months is highlighted by zero energy due to the lack of PM10 sources, i.e., industrial activity and road traffic. For *<sup>f</sup>*>1 <sup>10</sup><sup>6</sup> Hz, one can see more precisely energy variation related to dust events from mesoscale to large scale. Contrary to SPEC, HS clearly illustrates localized energy fluctuations due to small-scale event. In fact, the STFT makes an assumption that any signal as piecewise stationary and uses suitable window function to produce the short-time spectral characteristics of the signal. However, in **Spectrogram**

signal is composed [46].

**2005**

**2005**

**0**

**1**

**2**

**3**

**Frequency (Hz)**

**4**

**5**

**10-6**

**2006**

**2007**

**2008**

*Hilbert spectrum of PM10 time series with a color gradient from strong energy (in red) to weak energy (in blue).*

**2009**

**2010**

**2011**

**Figure 4.**

**Figure 5.**

**97**

**0**

**1**

**2**

**3**

**4**

**5**

**2006**

**2007**

**2008**

*Spectrogram of PM10 time series with a color gradient from strong energy (in red) to weak energy (in blue).*

**Hilbert Spectrum**

**2009**

**2010**

**2011**

**Figure 3.**

*The spectrum of wind speed in the Hilbert space and the Fourier space. A power law behavior is observed in both spaces.*

noticed in the Hilbert space with an estimated spectral exponent *β<sup>h</sup>* = 1.02 0.10. *β<sup>h</sup>* is equal to 1 power law scaling observed in the mesoscale range [41]. In the Fourier space, this power law is not significant. This is due to the existence of intermittent dust events with huge fluctuations in PM10 data (see **Figure 1**). Indeed, the Fourier transform is a linear asymptotic approach which requires high-order harmonic components to mimic nonlinear and nonstationary process [42]. Thus, the highorder harmonics may lead an artificial energy transfer flux from a large scale (low frequency) to a small scale (high frequency) in the Fourier space. Consequently, the Fourier-based spectrum may be contaminated by this artificial energy flux [42]. The artificial energy transfer may give a less steep power spectrum as we observed in **Figure 2**. By contrast, combined with the EMD method, HSA has very local abilities both in physical and spectral domains and does not require any higherorder harmonic components to simulate the nonlinear and nonstationary events. As a consequence, HSA method may provide a more accurate scaling exponent and singularity spectrum [42].

According to [43], wind speed dominates the amount of pollutant dispersion in the atmospheric boundary layer. In addition, this meteorological parameter could also transport PM10 from large-scale sources, i.e., African dust [9]. To complete our results, we used hourly wind speed measurements provided by the French weather office (Météo France Guadeloupe) located at Abymes (16.2630°N 61.5147°W). PM10 and wind speed measurements are very close, i.e., ≈8.1 km of distance, and performed at the center of the island under the same atmospheric conditions [2]. **Figure 3** illustrates the PSD provided by the Hilbert transform and the Fourier transform for wind speed data. This time, a power law behavior is observed in both spaces on the same frequency range 3*:*<sup>54</sup> <sup>10</sup><sup>7</sup> <sup>⩽</sup>*f*⩽1*:*<sup>36</sup> <sup>10</sup><sup>4</sup> Hz which corresponds to time scales 2*:*1 hours⩽*T*⩽32*:*7 days. Contrary to PM10 which is a passive scalar, wind speed is a vector quantity. The estimated spectral exponents are identical with, respectively, 0.89 0.06 and 0.90 0.12 in the Hilbert and Fourier spaces. As PM10, spectral exponent values are also close to 1. For wind speed, at low frequencies, a spectrum close to the 1 power law is likely occurs close to a rough surface, due to a strong interaction between the mean flow vorticity and the fluctuating vorticity [44, 45].

#### **4.2 Time-frequency domain**

The TFR in the Fourier and Hilbert spaces are, respectively, illustrated in **Figures 4** and **5**. Both figures show a color gradient from strong energy (in red) to *Spectral Observations of PM10 Fluctuations in the Hilbert Space DOI: http://dx.doi.org/10.5772/intechopen.88279*

weak energy (in blue). This highlights the energy activity related to PM10 concentrations during the study period. Such an approach gives the possibility of tracking the evolution of PM10 data spectral content in time, which is typically represented by variations of the amplitudes and frequencies of the components from which the signal is composed [46].

On **Figure 4**, strong energies are observed throughout the years with slight fluctuations on the frequency range 0⩽*f*⩽<sup>1</sup> <sup>10</sup><sup>6</sup> Hz. For *<sup>f</sup>*>4 <sup>10</sup><sup>6</sup> Hz, strong energies are also noticed in the middle of each year and at the beginning of 2010 with more fluctuations. In **Figure 5**, energy distributions are more localized. On the frequency range 0⩽*f*⩽<sup>1</sup> <sup>10</sup><sup>6</sup> Hz, we can observe the influence of small-scale event on energy behavior. As noticed, this energy may be weak or null. As an example, the impact of a general strike in early 2009 that paralyzed Guadeloupean archipelago at least 2 months is highlighted by zero energy due to the lack of PM10 sources, i.e., industrial activity and road traffic. For *<sup>f</sup>*>1 <sup>10</sup><sup>6</sup> Hz, one can see more precisely energy variation related to dust events from mesoscale to large scale. Contrary to SPEC, HS clearly illustrates localized energy fluctuations due to small-scale event. In fact, the STFT makes an assumption that any signal as piecewise stationary and uses suitable window function to produce the short-time spectral characteristics of the signal. However, in

#### **Figure 4.**

noticed in the Hilbert space with an estimated spectral exponent *β<sup>h</sup>* = 1.02 0.10. *β<sup>h</sup>* is equal to 1 power law scaling observed in the mesoscale range [41]. In the Fourier space, this power law is not significant. This is due to the existence of intermittent dust events with huge fluctuations in PM10 data (see **Figure 1**). Indeed, the Fourier transform is a linear asymptotic approach which requires high-order harmonic components to mimic nonlinear and nonstationary process [42]. Thus, the highorder harmonics may lead an artificial energy transfer flux from a large scale (low frequency) to a small scale (high frequency) in the Fourier space. Consequently, the Fourier-based spectrum may be contaminated by this artificial energy flux [42]. The artificial energy transfer may give a less steep power spectrum as we observed in **Figure 2**. By contrast, combined with the EMD method, HSA has very local abilities both in physical and spectral domains and does not require any higherorder harmonic components to simulate the nonlinear and nonstationary events. As a consequence, HSA method may provide a more accurate scaling exponent and

*The spectrum of wind speed in the Hilbert space and the Fourier space. A power law behavior is observed in*

**10-8 10-6 10-4 f, (Hz)**

**Wind speed**

**H( ) Hilbert spectrum slope H( ) -0.89 0.06 S(f) Fourier spectrum slope S(f) -0.90 0.12**

According to [43], wind speed dominates the amount of pollutant dispersion in the atmospheric boundary layer. In addition, this meteorological parameter could also transport PM10 from large-scale sources, i.e., African dust [9]. To complete our results, we used hourly wind speed measurements provided by the French weather office (Météo France Guadeloupe) located at Abymes (16.2630°N 61.5147°W). PM10 and wind speed measurements are very close, i.e., ≈8.1 km of distance, and performed at the center of the island under the same atmospheric conditions [2]. **Figure 3** illustrates the PSD provided by the Hilbert transform and the Fourier transform for wind speed data. This time, a power law behavior is observed in both

sponds to time scales 2*:*1 hours⩽*T*⩽32*:*7 days. Contrary to PM10 which is a passive scalar, wind speed is a vector quantity. The estimated spectral exponents are identical with, respectively, 0.89 0.06 and 0.90 0.12 in the Hilbert and Fourier spaces. As PM10, spectral exponent values are also close to 1. For wind speed, at low frequencies, a spectrum close to the 1 power law is likely occurs close to a rough surface, due to a strong interaction between the mean flow vorticity and the

The TFR in the Fourier and Hilbert spaces are, respectively, illustrated in **Figures 4** and **5**. Both figures show a color gradient from strong energy (in red) to

<sup>⩽</sup>*f*⩽1*:*<sup>36</sup> <sup>10</sup><sup>4</sup> Hz which corre-

singularity spectrum [42].

**Figure 3.**

*Functional Calculus*

*both spaces.*

fluctuating vorticity [44, 45].

**4.2 Time-frequency domain**

**96**

spaces on the same frequency range 3*:*<sup>54</sup> <sup>10</sup><sup>7</sup>

**10-2**

**10<sup>0</sup>**

**10<sup>2</sup>**

**S(f), H( )** **10<sup>4</sup>**

**10<sup>6</sup>**

*Spectrogram of PM10 time series with a color gradient from strong energy (in red) to weak energy (in blue).*

**Figure 5.** *Hilbert spectrum of PM10 time series with a color gradient from strong energy (in red) to weak energy (in blue).*

reality, most of air pollution signals are usually nonstationary [9, 14, 47]. The Fourier transform-based technique treats the signal as a sum of predefined basis functions. If the analyzing signal is well matched with the bases, it performs better; otherwise the performance is degraded [10]. Here, the SPEC highlight energy fluctuations linked to PM10 coming from African dust between May and September (large-scale sources) [9] and from the eruption of Soufrière on Montserrat in February 2010 (mesoscale sources) [48]. However, the SPEC does not detect energy fluctuation related to anthropogenic pollution, i.e., local sources. This shows HS is a robust method in timefrequency domain. Indeed, based on the EMD method, this TFR is fully data adaptive, and the signal decomposition is performed without any predefined basis functions. These results confirm the superiority of HS over STFT in TFR.

EMD empirical mode decomposition IMF intrinsic mode function HSA Hilbert spectral analysis HHT Hilbert-Huang transform TFR time-frequency representation

*DOI: http://dx.doi.org/10.5772/intechopen.88279*

*Spectral Observations of PM10 Fluctuations in the Hilbert Space*

STFT short-time Fourier transform

HS Hilbert spectrum

WT wavelet transform

*E*(*f*) Fourier spectral density

*C*(*t*) intrinsic mode function component

*r*(*t*) residual of the intrinsic mode function *φ* phase function of the intrinsic mode function

\*† and Rudy Calif <sup>2</sup>†

\*Address all correspondence to: karusphere@gmail.com

† These authors contributed equally.

provided the original work is properly cited.

)

1 Department of Research in Geoscience, KaruSphère SASU, Abymes, Guadeloupe,

© 2019 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/ by/3.0), which permits unrestricted use, distribution, and reproduction in any medium,

2 EA 4539—LaRGE (Laboratoire de Recherche en Géosciences et Energies), Département de Physique, Université des Antilles, Pointe-à-Pitre, Guadeloupe,

*ω* instantaneous frequency (Hz)

*N* total length of a sequence *x*(*t*) particulate matter signal (μg/m<sup>3</sup>

*f* frequency (Hz) *β* spectral exponent A instantaneous amplitude

*h*(*ω*) Hilbert spectrum

*j* scale index

**Author details**

Thomas Plocoste<sup>1</sup>

France

France

**99**

**Nomenclature**
