Detection of Extreme Phenomena in the Stable Boundary Layer over the Amazonian Forest

: We apply different methods for detection of extreme phenomena (EP) in air-turbulent time series measured in the nocturnal boundary layer above the Amazon forest. The methods used were: (a) a Morlet complex wavelet transform, which is often used in analysis of non-linear application processes. Through the use of the wavelet, it is possible to observe a phase singularity that involves a strong interaction between an extensive range of scales; (b) recurrence plot tests, which were used to identify a sudden change between different stable atmospheric states. (c) statistical analysis of early-warning signals, which verify simultaneous increases in the autocorrelation function and in the variance in the state variable; and (d) analysis of wind speed versus turbulent kinetic energy to identify different turbulent regimes in the stable boundary layer. We found it is adequate to use a threshold to classify the cases of strong turbulence regime, as a result of the occurrence of EP in the tropical atmosphere. All methods used corroborate and indicate synergy between events that culminate in what we classify as EP of the stable boundary layer above the tropical forest.


Introduction
Many of the techniques used for the analysis of micrometeorological data do not take into account situations where extreme phenomena (EP) are manifested. Here, EP are defined as sudden variations that occur simultaneously in several time series of thermodynamic and chemical species variables, measured at the forest-atmosphere interface. It is known that the manifestation of these phenomena in the Amazonian nocturnal boundary layer (NBL) is generally associated with strong turbulent exchange between vegetation and the atmosphere [1][2][3][4][5]. The detection and prediction of EP in geophysics, ecology, economy and other knowledge fields have already received attention [6][7][8][9][10], although these EP have different natures. The occurrence of an EP might be characterized as an emergent phenomenon in non-linear processes as discussed in a general way by [11][12][13] and in specific applications by [14][15][16], among others. Such processes are often associated with the emergence of critical slowing down, where a system is close to a critical tipping point and recovery rates of equilibrium states decrease, which have been observed in distinct scientific domains. Weng and Lau, and Lau and Weng [17,18] studied the occurrence of EP in the tropical atmosphere of Indonesia. They observed that during the outbreak of the EP, phase singularities (that involves a strong interaction between an extensive range of scales) appear in diagrams produced with Morlet's wavelet transform.
The recurrence plot (RP) technique is also useful to identify the occurrence of EP [19,20]. According to [20] the occurrence of these phenomena is always associated with sudden changes in the dynamics of the system, producing two qualitatively distinct states. They noticed changes on the texture patterns of the recurrence plot during the appearance of an EP. Trulla et al. [19] also verified the appearance of recurrence points (marked with darker colors in the RP) with the proximity of an EP. In addition to the aspects mentioned above, some studies reveal that the occurrence of such phenomena in the stable boundary layer are likely related to a transition that marks a change in turbulent regimes [21], or beyond transitions between regimes as verified by [22][23][24]. These authors mention the transition from weak to strong turbulence thresholds (referred here as λ 1 ). However, some recent results have revealed the need to define a second mode of transition threshold between nocturnal turbulence regimes (λ 2 ), associated with a threshold of transition between a regime of strong turbulence (regime 2 according to the classification of [25]) to another regime called here "very strong turbulence". Values greater than λ 2 can be observed in the SBL during the presence of convective clouds and strong downdrafts above the Amazon rainforest [3,26,27]. The complexity of the biosphere-atmosphere interaction processes are amplified by the presence of external forcing, such as convective clouds [26]. It is known that during the top-down movement of air parcels arising from higher levels of the atmosphere, which are colder and drier, going into the forest canopy, hotter air (and other particles inside the forest canopy) are being ejected from the forest canopy into the lower layers of the atmosphere.
The possibility of the existence of a second transition threshold was also reported by [28] in their analysis of the variability of carbon dioxide concentration statistics in relation to the average wind speed for the tropical forest region. Acevedo et al. [29], using data from experimental sites located in mid-latitudes, show that for situations of very strong turbulence regime there is an atmospheric condition in which they call "fully coupled state". Farias [21] reported that a very strong turbulence regime was associated with the occurrence of an EP. The association between EP warning conditions and turbulence regime changes might be used to improve turbulence mixing parameterizations in numerical models to better quantify the fluxes as well as the chemistry of reactive gases in the NBL, among other environmental indicators [4,30,31].
This work intends to show that during the occurrence of an EP in the stable boundary layer over the Amazon forest, it is possible to use different techniques both for the detection and for a better characterization of these phenomena. It should be noted that the occurrence of these phenomena seems to converge to the same condition of variability (to a singular point), such as: (i) the emergence of phase singularities in wavelets analyses, (ii) sudden variations in sensible heat and momentum fluxes near the surface and (iii) changes of the turbulent regimes (transition between the thresholds λ 1 and λ 2 ).

Data
The data used in this work were collected during the rainy seasons of the years 2013 (from January to May) and 2014 (from January to May and December). Data from the dry seasons of the two years were also analyzed, but events of EP were not detected. The site is located 60 km north of the city of Manaus at Cuieiras reserve (2 • 36.11 S-60 • 12.56 W, Figure 1), in a typical Amazonian terra firme forest [32]. In this work, five cases of EP occurrence were detected and analyzed, in 2 May 2013;  24 March, 13 April, 16 May and 8 December 2014. However, only the case of 13 April 2014, from 01 to 05 UTC (UTC = local time + 4 h) is presented in the manuscript. For other cases, we performed EWS tests and recurrence plot analysis for wind speed and temperature data. These cases are presented in the supplementary material. Time series of horizontal and vertical components of wind velocity (u, v, w), and the sonic virtual air temperature (T V ), measured at 48.15 m above ground with 3-D sonic anemometers (model CSAT3, Campbell Scientific Inc., Logan, UT, USA) at a rate of 20 Hz were used in the analyses. Also, ozone mixing ratio (O 3 ), humidity and pressure data measured at a height of 35 m were used, with a sampling rate of 1 Hz. To match this, the 3-D sonic anemometer data were later resampled to the sampling frequency of 1 Hz. For cloud cover analyses, available radar images, obtained every 12 min, were used. For more details see [26].

Methods
The methods used in this work are briefly presented in the following four topics. The first method uses time-frequency diagram of the phases of Morlet's wavelet coefficients, as used by [17]. The second deals with the analysis and classification of turbulence regimes in the tropical atmosphere. The third uses recurrence analysis applied to turbulent signals and was based on the work of [20]. The fourth deals with the early-warning signals and their possible applicability to EP that occur in the tropical atmosphere. Other detailed description of the methods used herein, are provided by [17], explaining the emergence of a phase singularity in a range of various scales during convective processes in the Indonesian atmosphere; by [25] for analyses of turbulence regimes; by [33] explaining the applications and limitations of early-warning signals in dynamic systems; and by [20] for a better understanding of the methods of recurrence analysis applied for complex systems.

Morlet Wavelet Transform
The wavelet transform decomposes a given signal s(t) into elementary functions ψ β,α (t) derived from a function called the mother wavelet. This decomposition is carried out by translation and dilatation of the "mother wavelet" and displaced along the time to generate a family of "daughter wavelets", which are expressed as function of two parameters: one for the scales (α), and other for the time position (β) [34,35]: where α, β ∈ R, with α > 0. A normalization factor (α) −1/2 is introduced in Equation (1) to maintain the energy of the mother wavelet across the entire wavelet family used. The wavelet transform of a real signal s(t) is defined in relation to the analysis wavelet ψ(t) as [18]: where ψ * is the complex conjugate of ψ defined in the "time-scale" real semi-plane (β, α). According to [35] it is possible to obtain the phase φ(β, α) of the signal, in an instant β and on a scale α according to the expression: In this work, Morlet complex wavelet was used, given by: Among the advantages of using Morlet wavelet, the following stand out: (a) its complex character that facilitates the analysis by scale of the signal phase [17]; (b) the possibility of detecting singularities observing the convergence of continuous phase lines in time-frequency phase spaces (which we will call here phase singularities, as defined by [17]; (c) allows the detection of frequency structure of non-linear dynamic processes containing quasi-periodicity and bifurcations [18].
Similar to the methods used by [17] in their analyzed cases of convective processes in the Indonesian atmosphere, we use Morlet wavelet transform to obtain the time signal phase, to construct an abstract singularity detection space (a phase space), which can be identified by the confluence of lines of the same phase in a time-scale diagram. Furthermore, in that work, several different examples of applications in geophysics are brought up, of what has been called "subharmonic cascade or period doubling route to chaos" in the specialized literature on non-linear dynamics [12] (chapters 13 and 14). By applying this methodology to our experimental data, we were able to characterize the region where there is an occurrence of EP (region where there is a phase singularity).

Turbulence Regimes
Historically, a phase space diagram has been used to classify different turbulence regimes that occur in atmospheric night environments. Such phase space consisted of a diagram proposed by [25], in which a velocity scale associated with turbulent kinetic energy ( ) was plotted on the ordinate axis as a function of the average wind speed (V M = (u 2 + v 2 ) 1/2 ) measured in the NBL, plotted on the abscissa axis. In this work this diagram was used to assess the threshold values of horizontal wind speed, λ 1 and λ 2 , which are associated with the speeds observed during the transitions between the weak to strong turbulence regimes and between the strong to very strong turbulence regimes, respectively. The values of V TKE and V M in this analysis correspond to averages of 5 min of data [24,36].
To determine λ 1 , mean values and standard deviations of V TKE were plotted against each interval of 0.5 m s −1 of V M over the entire available dataset. It was observed that in the region of V M = 2.2 m s −1 , an abrupt change in the mean values of V TKE occurred. In this way all situations in which the wind speed was above λ 1 = 2.2 m s −1 but less than λ 2 were associated with strong turbulent regime.
The second transition threshold, λ 2 , was determined from a plot of V M versus kurtosis of the concentration of carbon dioxide (K C ) (not shown here). We tried to find the limit value in which a change in trend in K C values occurs (growth starts to decrease), in relation to the average wind speed, considering that λ 2 represents the breaking point from which the basic state of the atmospheric system becomes disturbed, generating high variability in all the quantities measured at the surface. This new atmospheric state is characterized by average wind speeds higher than previously works conducted above Amazon forest [24,28]. According to [37], a significant increase in kurtosis of a turbulent quantity (from 2 to 3 times) provides a significant measure of the turbulence's intermittency. In addition, it can be particularly useful for "successfully distinguishing different types of turbulence" [37].
In our dataset, we found a significant change in K C values for horizontal wind speed above 4.6 m s −1 , and from this point on, λ 2 = 4.6 m s −1 is used to identify the beginning of the "very strong" turbulence regime.
For the analysis of the turbulent fluxes of sensible heat (H = ρc p w T ) and momentum (τ = ρw u ) associated with each regime, the covariance technique was used for each 5-min of data.
The equivalent potential temperature (θ E ) was calculated according to the thermodynamic equations used by [27] for the tropical atmospheric boundary layer.
where θ is the potential temperature, q V is the mixing ratio of water vapor and T LCL is the temperature at the lifted condensation level [27,38,39]. The value of θ E was calculated by aggregating contributions of quantities such as pressure and humidity [38] (Chapter 11) [40] (Chapter 4). It provides important information about the convective adjustment of the atmosphere in these high humidity environments [41].

Recurrence Plots for the Analysis of Turbulent Time Series
Recurrence plots are three-dimensional graphs, in which the time series to be analyzed is plotted on the x and y axes and the third dimension are the colors that vary between darker/colder tones (which represents high recurrence) and lighter/hotter tones (low recurrence). If recurring points intensify near an effective transition between two qualitatively distinct states (such as different turbulent regimes), changes in texture patterns are expected to identify this imminent transition. According to [20], the recurrence is a particularly useful method to analyze intervals of time in which the system "repeats certain phases of its trajectory".
For the construction of recurrence graphs, it is necessary to define some important parameters, such as the limit distance ( ) between two points in different trajectories, so that they are considered to be recurring, and the time delay (τ) used in the construction of the phase diagram, in its turn providing valuable information about the trajectory of the dynamic system over time [42]. The representation of each point x i (at time i) in relation to all other points x j (at time j) of the phase space is given mathematically by: is Heaviside function and the is a threshold distance, . is the norm and M is embedding dimension, determined by the false nearest neighbor algorithm [43]. The Heaviside step function delimits values between zero and one. The value zero indicates non-recurring and is represented by lighter/hotter tones and the value one refers to the maximum possible recurrence and is marked with darker/colder tones. The values associated with intermediate colors are incorporated as a point approaches or distances from that point in relation to which it is desired to estimate the recurrence. According to [42] the Heaviside function Θ is defined by: .
The value zero indicates non-recurring and is represented by lighter/hotter tones and the value one refers to the maximum possible recurrence and is marked with darker/colder tones. The values associated with intermediate colors are incorporated as a point approaches or distances from that point in relation to which it is desired to estimate the recurrence.
Here, were applied the methods presented in detail in [20], using the software provided by [43,44]. For the construction of the recurrence diagrams, the time lag τ = 8 s and the threshold = 0.1 for wind speed and temperature and = 0.6 for the O 3 mixing ratio were used. This " " threshold was carefully chosen to be not too small to exclude recurring points and not too large to include points of doubtful recurrence [20].

Early-Warning Signals
Many dynamical systems have transition points between different equilibrium states. Locating and understanding these transition points is particularly important to assist in the interpretation of phenomena that occur in nature, such as those that occur in the tropical atmosphere. One of the possibilities to analyze these transition points is using tools such as early-warning signals (EWS) [6,33,45,46]. These tools can be used for emphasizing changes in the properties of statistical categories such as an increase in the values of the autocorrelation function and variance. In this work, we used as a starting point the "early-warning signals toolbox" for detecting critical transitions in time series, provided by [33], whose main characteristics are: 1. One state variable must be selected. To verify which variable best suits as indicator, several tests were performed with the available variables (horizontal and vertical wind velocity, atmospheric pressure, O 3 mixing ratio and θ E , among others). 2. The autocorrelation at first lag (ACF1) was determined. According to [47] (p. 331), the autocorrelation function, ρ, is defined by: where z(t) is a stationary function at time t, which depends only on the difference (t 2 − t 1 ) between the two instants t 2 and t 1 in which z(t 1 ) and z(t 2 ) values have been measured. 3. The variance of signal was calculated. Mathematically, the unbiased estimator of variance can be given by: where µ is the mean value, z t is the value of the t elements and n is the data number of elements. The autocorrelation function and variance were calculated on a moving window with the length of 10% of the data series. The increment between successive moving windows was one point (1 s).
The objective is to verify whether before the occurrence of an EP a gradual increase of both the autocorrelation function and the variance of the turbulent quantities occurs. Gradual and simultaneous increases of the variance and the autocorrelation function of a state variable are precursors of a transition emergence between two physically distinct states [33,45,46].
To determine the robustness of the estimated trends (increase or decrease) associated with the autocorrelation function and the variance, the coefficient of Kendall [tau] was used. The Kendall [tau] indicator, as proposed by [48], considers a numerical sequence of "n" points which, in our case, is the number of points used in the autocorrelation function or in variance data series. Each numerical value in this series should be compared to all other subsequent values in the series, in order to form matching and discordant pairs.
The coefficient of Kendall [tau] varies between the limits of −1 to +1, being negative for a sequence of discordant pairs (decreasing trend) and positive for a sequence of matching pairs in which both members of an observation are larger than the pairs of the preceding observation (positive value). The value of Kendall tau is given by the ratio of the resulting score value among the n points of the data series in the form: where c is the number of concordant pairs and d is the number of discordant pairs [48]. We emphasize that the use of wavelet analysis and EWS had the same objective, i.e., to identify the presence of EP. Both methods were used so that there was no doubt regarding the occurrence of these phenomena in our time series. The use of recurrence plot allowed us to precisely locate the emergence of EP. Finally, the phase space of the turbulent regimes were used to show that during the occurrence of an EP there is a transition from a strong turbulent regime to another very strong turbulent regime.
In summary, the following steps were followed in this study, in order to confirm situations of EP occurrence in the tropical NBL: 1. All data were submitted to the quality tests proposed by [49] to verify spurious data in the available time series. 2. Nights with strong variations and sudden changes in turbulent signals of air temperature (>2 K), wind speed (>4.6 m s −1 ) and increases in O 3 -mixing ratio (>10 ppbv) were identified. 3. Selection of nights where there was a transition between the strong to very strong turbulence regimes (λ 2 ). 4. Selection of an appropriate state variable: a procedure consisting of tests to verify which of the available variables present the best results when submitted to the EWS tests (variables showing positive values for Kendall tau). For our study, three variables were used: wind speed, equivalent potential temperature and O 3 -mixing ratio. 5. Definition of the breaking point: the place where the EP starts, using recurrence diagrams (changes in texture patterns are expected to identify this imminent transition), and phase diagrams using Morlet wavelet coefficients (emergence of a phase singularity along the scales from 8192 s to 256 s ). Several studies have shown similar simultaneous variations in the presence of deep convection [3,4,26,27]. These strong variations in the time series can be detected by a series of other techniques that will be shown below. It is observed that the trends associated with the autocorrelation function and the variance are in accordance with what is expected for detection of early warnings before abrupt changes (such as the occurrence of EP) in several kinds of temporal processes [33,50], since they tend to grow as an EP approaches. The positive values of Kendall [tau] coefficients qualify the robustness of the indicators along the time series of the state variable [33,50].

Strong Variations in the Turbulent Time Series and Early-Warning Signals
The strong variations observed in the thermodynamics time series and O 3 -mixing ratio near the surface (Figure 2) and the increase in the autocorrelation function and variance values (Figure 2a,c,e) were associated with occurrence of an EP and whose origin is possibly related to the action of deep convective clouds upon the area of the experimental site. Figure 3 shows the cloud cover in the UTC time corresponding (approximately) to the occurrence of an EP for two (13 April and 16 May 2014) of the five days in which the phenomenon was observed. The radar images for other three days (not shown here) also show an extensive cloud cover above the Cuieiras reserve, during the EP occurrence.
The phenomenon studied here involves the action of convective clouds [27,51]. According to [27], "the fall of θ E and rise of O 3 -mixing ratio are tightly coupled but exactly out of phase, as bursts of downdraft air reach the surface, indicating that the same vertical transport process is responsible for both changes". Also, we observed the occurrence of sudden variations in the pressure signals that became more evident at 35 m height (Figure 2d) (a trend already described by [29] for another set of meteorological data), a few moments before the outbreak of an EP. Such variations in pressure signals could be associated with downdrafts and updrafts acting over the experimental site, as mentioned by [52] in their study on aerosols evolution cycles above the Amazon forest.

Wavelet and Recurrence Plot Analyses
Morlet wavelet analysis was used to construct a phase space similar to that performed by [17].  (Figure 4d). Figure 4c shows the existence of a main phase singularity along the D 3 axis, which encompasses a range of time scales, from the highest scale (8192 s ≈ 2.27 h) to the lowest verified scale (256 s ≈ 4.27 min). Also it is possible to verify that there is a drop in temperature that occurs simultaneously with the maximum wind speed (Figure 4b, close to 0330 UTC, indicated by the arrow D 3 ) and a strong interaction among several physical scales (Figure 4d), with a brief detectable bifurcation below the 512 s scale (≈8.53 min). It is possible to notice the emergence of a phase singularity in a range of various scales, what was discussed by [17] in their analyzed cases of convective processes in the Indonesian atmosphere. They are expressions of non-linear dynamics that can be found in tropical atmospheric dynamics. Figure 5 presents three recurrence graphs corresponding to the three turbulent variables (V, θ E and O 3 mixing ratio, respectively) selected for this study. The time interval used for the recurrence analyzes was from 01 to 05 UTC. A change of dynamics in the data series can be observed around 0330 UTC for all three turbulent variables used. A sudden change in the dynamics of the system can be observed around 0330 UTC, in which there is an intensification of hotter bands in the recurrence graph for the three turbulent variables used. After this time (around 0330 UTC) the structure resembling a chess board that individualize recurrence processes are no longer verified. These changes in texture patterns provide indications of the occurrence of a transition between two qualitatively distinct states.  It is important to note the structures in colder colors (preceding the EP) that can identify the recurrence process [19,20,42] are suppressed after 0330 UTC. Also, at the same time (0330 UTC) there is the existence of a phase singularity (Figure 4c,d). Therefore, it can be assumed that a tipping point (for this particular case) occurred around 0330 UTC.

Turbulence Regimes
To better explain some interesting features of distinct turbulent regimes observed in the NBL, Figure 6a shows the intensity of the turbulence (V TKE ) as a function of V M . For the construction of Figure 6, data from all 5 investigated nights were used, each point corresponding to averages of 5 min. It is possible to notice the presence of several turbulence regimes, which are clustered in 4 regions: around 3 distinct fitted segments (0A, AB and BC) and within the area delimited by the ellipse D, which are present in the Figure 6a. The black circles represent occurrences of the regime of weak turbulence (regime 1, as defined by [25]); the green diamonds represent occurrences of a top-down type of turbulent regime (regime 3, as defined by [25]); the blue triangles and the magenta squares represent occurrences of two kinds of strong turbulence regime (regime 2) that occur in the absence or presence of extreme phenomena (EP's), respectively. Here, the strong turbulence regime (regime 2, as defined by [25]) was split into two sub-sets with different slopes, the first one representing a type of strong turbulence regime (wind speed greater than λ 1 = 2.2 m s −1 and less than λ 2 = 4.6 m s −1 ), not associated with EP, and the magenta squares marking the points of strong turbulence regime in which the value of V M was greater than λ 2 = 4.6 m s −1 and corresponding to situations in which the wind speed was maximum for the analyzed nights. These points ("magenta" squares) were associated with the occurrence of an EP and whose origin is possibly related to the action of deep convective clouds near the experimental site. The occurrence of a very strong turbulence regime, by itself, seems not enough for the generation of an EP. However, it seems that all EP's such as those studied here could be classified as belonging to very strong turbulent regime. Figure 6b shows the behavior of the flux H over each one of the turbulence regimes presented in the Figure 6a. It is possible to verify that the more intense turbulence regimes are also responsible for the higher intensities of the sensible heat fluxes, which reach −80 W m −2 in the very strong turbulence regime. In Figure 6c, the intensity of the vertical flux of u-momentum follows a growing trend for the three turbulence regimes identified in Figure 6a. The red lines shown in Figure 6b,c correspond to the average values of H-turbulent fluxes and momentum flux obtained at each 0.5 m s −1 interval for the wind speed and their respective standard deviations. However, for the very strong turbulence regime, an interval of 0.8 m s −1 was used due to less data. It is noted that this increasing tendency in momentum flux is not very evident during the weak turbulence regime (black circles), but it becomes clear with the occurrence of the strong turbulence regime (blue triangles), following the same trend during the occurrence of the "very strong" turbulence regime, up to the 7.0 m s −1 observed upper limit of V M velocity.
These results corroborate those presented by [24] for data from Rebio Jarú reserve in South-Western Amazon in which the authors found that there are approximately 40-fold increases in the turbulent sensible heat fluxes during the occurrence of strong turbulence regimes, as compared to weak turbulence regimes. Here, however, it is possible to verify that this increase in H values reaches its maximum value after the second transition (λ 2 ), i.e., between the strong turbulence regime and a "very strong" turbulence regime.
After the transition from the strong to the very strong turbulence regime, there is a reduction in the intensity of H values with higher mean horizontal wind speeds. This tendency of reduction in the fluxes was also observed by [28]. These authors were able to obtain subsidies that could have led them to the proposition of the existence of a second type of regime of strong turbulence above the Amazon forest. They found a maximum horizontal wind speed of 3.5 m s −1 for nighttime data measured at 39.4 and 46.0 m (in the Uatumã reserve, northeast of Manaus, where the Amazon Tall Tower Observatory is located-ATTO). They also used V TKE proposed by [25] to study the variability of CO 2 mixing ratio and suggested the existence of two different patterns of CO 2 there. But they did not associate their results with the occurrence of EP's above the experimental site. It should be noted that for the case of momentum flux shown in Figure 6c, its intensity does not decrease beyond the λ 2 values, showing only a small reduction in its growth rate immediately after V M ≈ 4.6 m s −1 . Possibly, the higher increase trend is because vertical and horizontal velocity components continue to grow until the climax of the EP occurrence, thus raising the value of momentum flux. In the case of H, the gradual increase of the wind velocity, associated with the variations in the other turbulent quantities, acts at first to increase the intensity of the fluxes; however, after a certain threshold (λ 2 ), the wind speed starts to act to reduce the average value of the fluxes instead of increasing it, as was also verified by [28].
Therefore, it is observed that extreme phenomena are atmospheric processes that drive transitional changes in the turbulent regimes and where several phenomena are synchronized, such as: (i) a clear separation between two regions, one of them with smaller air density (before the occurrence of EP) and the other with greater air density (during the occurrence of EP) (Figure 2f) a region where there is a phase lines convergence ("singularities", in the terminology proposed by [17] (Figure 4); (ii) In this region the state variables (in this case, the equivalent potential temperature) validate the criteria formulated in the methodology of [50] (Figure 2c). These transitional changes observed in measurements made at 48 m above the surface, are also noted at other heights, located both above and inside the forest canopy (not shown here), with dramatic shifts in their values during the moment of occurrence of the EP.

Conclusions
In this work, experimental data measured in the nocturnal boundary layer, above the Amazon rainforest, were used to analyze strong variations that occur simultaneously in time series (wind speed, temperature, ozone, air humidity, among others). These variations were called extreme phenomena as they significantly affect the turbulent motion in the local boundary layer. The outbreak of these phenomena occurred during the presence of strongly convective clouds, with their respective downdrafts and, consequently, producing strong and well-located increases in pressure variance and air density. Also, it has been shown that these extreme phenomena can be detected by different techniques. One of them analyzed early-warning signals, where there was an increase in the autocorrelation function and in the variance of the state variables immediately before the onset of an extreme phenomena. Morlet wavelet transform is also an efficient tool to identify the beginning of extreme phenomena, since there occurs an intense phase singularity at different frequencies. Also, the recurrence plot technique shows a strong change in atmospheric dynamics, i.e., a transition between qualitatively different atmospheric states. Such a transition is generally expected during the outbreak of an extreme phenomena. Finally, during the occurrence of the extreme phenomena a new category of turbulence regime emerges, classified here as "very strong", where the turbulent fluxes of heat and moment undergo changes completely different from those observed in the absence of the extreme phenomena. All these results demonstrate that during the outbreak of the extreme phenomena, there is a synchronization of events that converge to the same condition of variability.