1. Introduction
Cyclones are extreme atmospheric conditions that have significantly contributed, over the course of time, to the evolution of climate change and its impact on the Earth’s system. It is well known [
1] that cyclone events are triggered by high temperature air masses following a rotational motion developed around a low pressure center. They are driven by a strong thermal source, which continuously supplies energy in the form of latent heat, so they generate spiral motions of the atmosphere, causing rotational dynamics in the geographic region concerned. These events are mainly formed in tropical latitudes; hence, the name
tropical cyclones. According to their speeds of motion, they assume particular characteristic conditions; moreover, they are cataloged in classes of strength, according to the Saffir–Simpson scale [
1,
2]. A hurricane is defined as a tropical cyclone that has reached the highest class, corresponding to category 5.
The strong pressure differences among air masses rotating around the low pressure center generates turbulence in the surrounding regions [
3,
4,
5]. In this work, we studied the characteristic properties of turbulence at different spatiotemporal scales within the regions of hurricane
Faraji, the first tropical cyclone of the 2021 season that reached category 5 on the Saffir–Simpson scale; it developed in the Indian Ocean between 4 and 13 February 2021. We would like to point out that the word “turbulence” is used hereinafter in the generic sense of “random superposition of different length scales”. The actual mechanisms that can produce this small-scale formation (turbulent cascade, local instabilities, convection, etc.) are not explored in detail in our investigation. The study uses data and imagery collected by the
SEVIRI radiometer [
6], installed on the
Meteosat Second Generation-8 (IODC) (MSG-8), an
EUMETSAT geostationary satellite. With these products, remote sensing techniques were used to obtain a series of
N images, made of matrices of
pixels, where
and
are, respectively, the resolution in longitude and latitude, in a given period of acquisition. Each pixel (representing a total area of 1 Km
for SEVIRI) is assigned a brightness temperature value, so that the image is a matrix representing the temperature field of the cyclone observed at a given time. Starting from this scalar field, we carried out the analysis of the spatiotemporal dynamics of cyclones using the
proper orthogonal decomposition (POD) technique [
7,
8,
9,
10]. Thus, the original field was decomposed on the basis of empirical eigenfunctions extracted directly from the data by diagonalizing the autocorrelation matrix of the temperature field, integrated on the whole domain (Fredholm problem of the second kind). The spatial and temporal components were factorized in the POD expansion by attributing the spatial dependency to the eigenfunctions and the temporal evolution to the coefficients of the development.
The POD technique has several advantages over other traditional techniques (such as Fourier or Wavelet analyses): (a) it does not assume any special behavior at the boundaries of the domain; (b) the empirical basis is optimal, in the sense that it captures (on average) the maximum possible energy on each mode, which means that the POD spectrum has the most rapid decrease as a function of the discrete number of the eigenfunctions; (c) the POD eigenfunctions are orthonormal and the temporal coefficients are decorrelated from each other. When studying turbulent motions, especially in the non-homogeneous case, the POD decomposition acquires a special value, as it is well known that turbulence in shear layers is characterized by the presence of “coherent structures”, namely energy-containing, velocity field structures that persist in the flow for several nonlinear turnover times. As a consequence of the property (b), the POD eigenfunctions in this framework are correlated with the turbulent coherent structures. In the case of statistically homogeneous isotropic turbulence, the coherent structures are plane waves; in fact, the POD decomposition applied to the periodic fields yields, as a result, the Fourier expansion.
However, concerning differences with expansions in which the basis is fixed a priori, as with the Fourier or Wavelet decompositions, no characteristic scale can be intrinsically associated directly with the POD eigenfunctions, in the sense that each eigenfunction does not have, in principle, a single characteristic scale, but it may be composed of many of them, all mixed together. Moreover, since the basis is extracted empirically directly from the data, the energy spectra of two different datasets cannot be easily compared directly, since the eigenfunctions themselves, and not only the coefficients, contain the information concerning the energetic budget.
Through the POD analysis, and from the analysis of the spectra, we infer that the turbulent dynamics of cyclones appear to be characterized by three different slopes of the spectra, at three different ranges of characteristic scales that we call macroscale, mesoscale, and small-scale, each one likely associated to different atmospheric properties within the cyclone and characterized by a different evolution.
The remainder of the manuscript is structured as follows. In the following section, we briefly describe the technique used for obtaining the data through a remote sensing analysis and, hence, we recall the main properties of the POD decomposition, which will be used throughout the rest of the article. Then, we present the results of the analysis carried out on the Faraji cyclone. Finally, we present our conclusions and outlooks from the present study.
3. Results
In
Figure 2, we show a contour plot of the time-averaged temperature field
, as a function of the
(latitude and longitude) coordinates, to have a physical representation of the average structure of hurricane Faraji. The hotter central zone coincides with the “eye” of the cyclone, surrounded by colder zones. One can notice that high-temperature spots are however present, not only in the middle of the domain, but also in the upper part of it, yielding a complicated picture of the structure of this hurricane.
To analyze the existence of a direct relationship between the order
j of the POD eigenfunctions and the spatial scales—not defined a-priori because the order of the eigenfunction is connected to the energy content of the mode, not to the characteristic scales—we calculated the Fourier spectrum, shell-integrated in the wave vector space
, for each snapshot of the temperature field. To this aim, we computed the two-dimensional fast Fourier transform (FFT), of the temperature field. Then we divided the
plane into two-dimensional concentric circular regions (shells). The
i-th shell had an internal radius equal to
and an external radius equal to
. The width of all the shells was therefore equal to
. Then, we attributed the
i-th shell with a spectral energy equal to the sum of the energies corresponding to values of
, such that:
, with
. Then, we computed all contributions of energy coming from different shells as a function of
. Finally, the shell-integrated spectrum was averaged at the different times, to evaluate the mode energy content in the dynamics factorizing the spatial dependence. We plotted this spectrum for the Faraji cyclone in
Figure 3.
From the plot in
Figure 3, one can observe that the spectrum exhibits three distinct slopes; thus, allowing us to make a first classification of the spatial dependency of the eigenfunctions into three different zones, which we associate to three corresponding ranges of characteristic scales: macroscales, mesoscales, and small-scales. Such ranges of scales correspond, according to the limits of the fits, to scales larger than ∼100 km for the macroscales, between ∼22.6 and ∼100 km for the mesoscales, and smaller than ∼22.6 km for small-scales. In the plot, we highlighted the spectral indices of the fitting power law
for the macroscales, mesoscales, and small-scales, which go approximately as
, respectively.
The slopes of the three different zones and their limits have been determined in the following way: suppose we have a spectrum of as a function of , which behaves as a power law in three distinct zones (let us call them Z1, Z2, and Z3). By taking the logarithms of both E and , the spectrum in the three zones behaves like a straight line , where y represents , x is , is the slope of the power law, and is a parameter depending on the energy contained at some arbitrary value of . In order to determine the spectral slopes and the values of at which the transitions among the different zones happen, we separately consider zones Z1 and Z2 first, then zones Z2 and Z3.
We fit the first two zones with a linear function
, defined as:
where
and
are the slopes of the two different power laws in Z1 and Z2, respectively,
is the value where the change in the slope between the two zones takes place, and
and
are the constant parameters of the straight line fits. By imposing that the two linear functions are to match their values at
, we can see that not all the parameters are independent. Therefore, we can express one of those parameters, for instance,
, as a function of the remaining ones (for instance,
,
, and
in this case). This can be obtained through the relation:
By substituting this relation in the expression of
, the final form of the fitting function is:
Finally, we fit the quantities:
versus
in Z1 and Z2 against the function defined in Equation (
9). This is a four-parameter fit, which yields the values of the slope
, the value
where the change in the slope happens, and the parameters
and
, as well as the value of the slope in the second zone
through the relation (
8).
After determining the slopes and the break in the spectrum in zones Z1 and Z2, we repeat the procedure for zones Z2 and Z3, to find the slope in Z3 (the one in Z2 is of course the same obtained in the first part of the procedure) and the second spectral break. This procedure gave us the values of and the point in which the separation among macro-, meso-, and small-scales takes place, described above.
The three-scale separation hypothesis was further confirmed by constructing the spectra, shell-integrated in the Fourier space, corresponding to each of the POD eigenfunctions. In
Figure 4, we show the example of three eigenfunctions, corresponding to
,
,
, which we associate to the macro-, meso- and small-scales, respectively, for the Faraji cyclone. On the left panels, we show the contour plots of the eigenfunctions as a function of longitude and latitude, while on the right panels, we show the shell-integrated spectra of the three eigenfunctions. It is clear that, with increasing
j (from top to bottom), the characteristic scales become smaller. This is also visible from the spectra shown in the right panels, in which the peaks corresponding to the typical length scales are clearly visible and they shift continuously towards higher wavenumbers, that is, smaller scales, as
j increases.
In
Figure 4, we print the maximum value of
, namely the wave vector of the peak of the spectra, which identifies the characteristic length scale
ℓ of the eigenfunction (
), normalized to the total dimension of the image
. Although the peaks are rather well visible in the spectra, they do not correspond to a real delta function, which is a single plane wave in Fourier space. This happens essentially for several effects: (1) the data (as, consequently, the eigenfunctions) are not periodic; therefore, some spurious energy is found at the smaller scales (higher wave vectors) due to this effect; (2) we compute a discrete Fourier transform, that is, when the energy is not exactly concentrated on a single wave vector
k, which corresponds to a integer multiple of the fundamental wavelength
, the energy itself is spread over the modes close to
k; (3) as it is visible in the contour plots of the eigenfunctions, the structures observed may have different characteristic scales along the longitude and latitude directions, which can produce a broadening of the peaks when the spectra are shell-integrated over the two directions. Nevertheless, we observed that the trend corresponding to a continuous shift of the peaks of the spectra towards smaller characteristic length scales is present for all values of the POD eigenfunction order
j.
Afterwards, we carried out the same kind of analysis on the temporal component of the POD spectrum, namely, we computed energy spectra for the temporal evolution of the
j-th POD coefficient
, by Fourier-transforming in time the function
for all the values of
j. As done for the spatial analysis, we plotted the time evolution of the POD coefficients
,
and
, along with the corresponding spectra in
Figure 5.
It is clearly visible from the plots in the left panels that the number of oscillations in time increases for increasing values of j, meaning that higher js correspond to higher frequencies of the POD coefficients. This is even visible in a more clear way by looking at the spectra of the coefficients, plotted on the right panels. Again, the spectra exhibit rather well defined peaks, which move towards higher frequencies with increasing values of j. We also notice that the peak does not correspond to a delta-like function. Once again, this may depend on the fact that the time behavior of the coefficients (as visible in the left part of the figure) is not really periodic, and that the true frequencies are not necessarily integer multiples of the fundamental frequency , P being the total observation period, given by: min.
As seen in
Figure 5, we found that there was a monotonic shift towards the highest frequencies in the maximum values of the temporal spectra of the POD coefficients with increasing values of their order
j. We therefore identified for each value of
j the peak in the frequency spectra and plotted its value against
j (not shown here), and we found that there is a linear relation between the two quantities, as already found for the spatial scales in the Fourier spectra in
. Therefore, we can easily associate, through a linear best-fit procedure, a value of frequency to each POD mode with index
j, in analogy with what we did for the length scales. Starting from this result, in order to better characterize how the energy of the field is distributed as a function of the frequency, we evaluated the average energy spectrum of the POD coefficients for the temperature field, inferred from the remote sensing analysis, of the Faraji cyclone, and we plotted it in
Figure 6, on a double logarithmic scale, as a function of the frequency
. In the plot, we calculated the energy as the eigenvalue
, since this quantity represents the energy on the
j-th mode averaged in time. The frequency, instead, is computed from the peaks of the Fourier spectra of each POD coefficient
inferred through the best-fit procedure explained above.
From the plot, it is clearly visible that, also in this case, we can distinguish three different zones, each one approximately described by a power law with different scalings, which correspond approximately to the same POD modes observed in the spatial spectra, confirming that a linear correspondence between spatial scales, temporal scales, and POD order of the eigenvalues/eigenfunctions exists. Therefore, we can assume that, in the frequency domain, the subdivision we made for the spatial scales in macro-, meso-, and small-scales, may hold. However, in the case of the frequency spectrum, the differences among the three zones are not so evident as in the cases of the spatial scales. The slopes of the power laws of the energy E as a function of the distribution of the frequencies are all quite close to −2, which means that the energy content decays continuously, as for all three characteristic scales. The characteristic time scale t corresponding to the macro-, meso-, and small-scales are, respectively: h, and h.
The determination of the slopes in the frequency spectrum obtained through the POD coefficients and the separation among the three zones were carried out in the same way as we already did in
Figure 3 for the averaged spectrum of temperature as a function of the spatial scales.
Finally, exploiting the linearity of the POD development, we can use Equation (
6) to reconstruct the partial sums of the POD modes at different times for values of
j in the macro-, meso-, and small-scales, to study how the turbulence evolves in the three different ranges of scales. The basic idea of this approach is that the turbulent evolution in the three zones is likely due to different mechanisms, operating at different length scales. We show examples of such reconstructions in
Figure 7, which represents the spatial structures in the three different zones corresponding to the macroscales (panel a), the mesoscales (panel b), and the small-scales (panel c), for the data at 21.00 UTC on 8 February 2021.
From the reconstructions, the scale separation makes it possible to distinguish the differences in the internal structure of the cyclone and follow the evolution in time. With the large scales (panel a), we can identify the macroscale structures of the cyclone, while in the mesoscale (panel b), some filament structures remain, all rotating around the eye of the cyclone and with a tendency to expand towards the upper-right part of the cyclone. Some traces of those filament structures persist at smaller scales (panel c). However, here, the most evident phenomenon is the presence of small-scale spots, mainly spread in the zone of the cyclone external to the eye.
4. Discussion
In this work, we showed that, using data acquired from satellite missions, with the help of specific observational, acquisition, and processing procedures, it is possible to obtain significant results to understand the turbulence in extreme cyclonic phenomena. In particular, by exploiting remote sensing techniques, we captured a certain number of snapshots that were present in a given region of the globe at specific time instants, assigning temperature values to the pixels making up the reference matrix, from which it was possible to acquire a large amount of information on the physical nature of the phenomenon.
Specifically, the use of images acquired by sensors on geostationary satellites, such as SEVIRI, made it possible to apply a cyclone analysis from a purely dynamic point of view to an observational scalar field, such as the brightness temperature. We point out that, in this kind of investigation, the information along the line-of-sight is integrated, yielding a two-dimensional dataset at different sensing times. We used the POD technique, an analysis widely employed to study turbulent regimes in fluids, to extract useful information about the most energetic structures present in the evolution of the temperature field inside the Faraji hurricane, a category 5 cyclone on the Saffir–Simpson scale.
This technique makes it possible to separate the spatial and temporal behaviors in the evolution of the temperature field of the cyclone, starting from temperature variations, and to understand how these can influence the energy distributions at different scales. By extracting the POD eigenfunctions and coefficients through a suitably written numerical code, we found a correspondence between spatial scales and temporal scales in the POD modes. In particular, we found that, when increasing the order of the POD eigenfunction, both the typical length and time scales of the structures became smaller. From the differences in the slopes of the Fourier spectra in both the wave vector and frequency domains, it was possible to infer the presence of three different zones, which we called macro-, meso-, and small-scales, characterized by power laws with different spectral slopes. However, the differences in the slopes were much more evident for the spectra in the wave vector space than in the frequency domain, where, instead, the scaling of the energy was approximately proportional to .
There was another important observation concerning the reconstruction of the temperature field at the three different length scales—the macro-, meso-, and small-scales were mainly localized at different positions in the cyclone, which shows that their evolution was likely due to different dynamical mechanisms.
Different spectral slopes are commonly observed in the atmosphere. We cite, among others, the recent articles by Larsén et al. (2016) [
21], Nosov et al. (2019) [
22], and Shikhovtsev et al. (2019, 2021) [
23,
24]. However, a direct comparison of our results with the ones obtained in those papers was not easy to carry out because of the different conditions in the methods applied. In our case, we obtained the spectra of the temperature field in a hurricane, which is a highly nonstationary phenomenon, whilst those studies were related to measures of wind speed in standard atmospheric conditions. Moreover, they performed direct measurements of such quantities, whilst we obtain the temperature fields as indirect measures from the remote sensing analysis. Finally, our spectra were obtained by using the POD analysis, which is quite different from the standard Fourier analysis from which the spectra were obtained in those articles. However, a comparison of our work with the results of those articles can be very interesting and deserves a dedicated investigation.
Although this preliminary study supplies important information about the possibility of understanding, in more detail, the origin and evolution of a cyclone, drawing definitive conclusions concerning the latter point requires more extensive studies. First, the same analysis should be repeated on several different cyclones, with the aim of understanding whether the features found in the case of the Faraji hurricane are specific to this particular cyclone, or it is possible to extend this dynamical behavior to other cyclones of the same category, or of different categories, on the Saffir–Simpson scale. Moreover, it would be necessary to follow the time evolution of the structures in the three different zones of the cyclone to understand whether it is possible to understand the origins of the different dynamical behaviors that seem to characterize them. All of those studies will be the subjects of forthcoming papers.