Near-Surface Atmospheric Turbulence in the Presence of a Squall Line above a Forested and Deforested Region in the Central Amazon

: Squall lines (SLs) are convective systems that cause heavy precipitation and consequently modify the atmospheric thermodynamic structure near the surface. SLs generated along the northern coast of Brazil and their effect upon atmospheric structure during their westward displacement into the Amazon are studied. Satellite imagery was employed to identify an SL above two experimental sites in the central Amazon and to characterize differences in the near-surface turbulent and ozone exchange during the passage of the SLs. The two sites, which are separated by about 100 km, feature contrasting vegetation. One site is tall canopy rainforest and the other is deforested. From our case study, it is noted that: equivalent potential temperature signiﬁcantly drops, principally in the forested region; the average near-surface wind speed increases 5 fold; the skewness of vertical wind velocity becomes considerably negative; signiﬁcant increases in turbulence intensity are observed. These changes suggest the presence of strong downdrafts generated by the SL. Shear production and dissipation rate of turbulent kinetic energy are considerably larger during the SL when compared to periods with absence of SL. In this study, we show that SLs are capable of modifying the vertical organization of the turbulence over forested and deforested areas, leading to changes in certain chemical processes that occur near the surface. To the best of our knowledge, this study represents a ﬁrst in demonstrating that near-surface turbulent ﬂow in the Amazon region is modiﬁed by the presence of SLs.


Introduction
Approximately two-thirds of the global precipitation occurs in tropical regions, primarily caused by Mesoscale Convective Systems (MCSs) ( [1][2][3], among others). MCSs are characterized as organized propagating storms that can last for at least several hours resulting in intense rain events, strong winds and, in many cases, hail storms [4]. Squall lines (SLs) are a common type of MCSs, in which a series of storms are arranged in a linear manner, at scales greater than 100 km [4,5]. Such SLs can be identified using satellite images and their displacement can be tracked. In the Amazon, SLs often originate along the northern Atlantic coast. These storms may either propagate very large distances (greater than 1000 km) from the region of origin or they may dissipate close to the coast [6][7][8].
Several studies have shown that SLs can cause intense downdrafts in the Amazon region [9][10][11][12]. It is already known that these downdrafts produce a rapid decrease in surface values of the equivalent potential temperature (θ e ) ( [13,14], among others). Furthermore, multiple authors have documented that such downdrafts can alter Ozone (O 3 ) levels close to the surface [15][16][17][18]. Scala et al. [16] used data from the second NASA Amazon Boundary Layer Experiment [19], to estimate O 3 transport during an SL occurrence. They observed that updrafts and downdrafts altered O 3 levels in both the lower and upper troposphere. Using data from the Large-Scale Biosphere-Atmosphere experiment [10], Betts et al. [17] showed that surface measurements of θ e and O 3 were altered by convective downdrafts. This downward movement, as they describe, had the ability to transport air from higher altitudes to the surface, with a higher O 3 concentration and a lower θ e relative to values at the surface. Recent experimental studies have shown that downdraft-mediated O 3 transport from the free atmosphere to the atmospheric boundary layer (ABL) can increase the surface levels of O 3 by 20 to 30 ppbv [18,20,21].
In the literature, there are multiple studies that address additional aspects, characteristics, and environmental conditions associated with the occurrence of SLs in the Amazon ( [6,9,22,23], among others). Negron-Juarez et al. [23] showed that patches of collapsed trees had a relationship with the so-called "blowdowns", which are generated during the passage of SLs.
Various studies employing experimental observations or numerical models have shown that the presence of SLs intensify the turbulence in its vicinity [5,[24][25][26]. However, to date, there is no study detailing the role of SLs in the turbulent exchange processes and in the structure of turbulence at the forest-atmosphere interface, such as the Amazon region. Therefore, this work presents the first detailed case-study investigating how the presence of SL-generated-downdrafts alters the processes of turbulent exchange and transport of O 3 [27] over forested and deforested surfaces.

Experimental Site
The data used in this work were collected at the Cuieiras (2 • 36.11 S-60 • 12.56 W) and Manacapurú (3 • 12.80 S-60 • 35.92 W) sites, hereafter designated T0k and T3, respectively ( Figure 1). The site T3 is approximately 70 km from the city of Manaus, State of Amazonas, Brazil, in a grassland region of 2.5 by 2 km, surrounded by old-growth forest [28]. The data were obtained as part of the US Atmospheric Radiation Measurement (ARM, https://www. arm.gov/, accessed on 1 March 2021) program, during the GoAmazon (Observing and Modeling of the Green Ocean Amazon) experiment [28].
The T0k site is located approximately 60 km Northwest of the city of Manaus in an area of dense tropical rainforest with a complex topography composed by valleys and plateaus, with a maximum height difference of 60 m. The plateaus have a maximum tree height of about 30-40 m and the soil is dominated by clay. In the valleys, the vegetation grows over sandy soils and is less dense as well as having a lower canopy height [18]. The experimental data were collected between March and December 2014, during the GoAmazon Project experimental campaigns [27].

Satellite Imagery, Ozone, and Meteorological Variables
Satellite imagery data are used to identify and monitor the propagation of SLs. The GOES-13 (Geostationary Operational Environmental Satellites) channel 4 enhanced infrared images were obtained from the website of the Comprehensive Larga Array-data Stewardship System of National Oceanic and Atmospheric Administration (CLASS-NOAA, http://www.avl.class.noaa.gov, accessed on 1 March 2021). The satellite images used to identify convective clouds and SL are in the infrared band often referred to as cloud top temperature and have been employed in numerous studies of deep atmospheric convection [29][30][31].
SLs were identified and classified following [6]. This classification depends on the horizontal displacement of the SLs into the continent, being: (1) Coastal squall lines (CSL): those whose horizontal propagation to the interior of the continent but not exceeding 170 km. (2) Type 1 Squall lines (SL1) have horizontal displacement between 170 and 400 km.
(3) Type 2 Squall lines (SL2) travel more than 400 km inland. In this study,an SL2, which occurred on 20 April 2014, was investigated as a case study. The SL2 selected formed on the northern coast of Brazil and traveled towards the sites under study (≈1350 km). It should be noted that the intense rainfall associated with the SLs impairs the performance of sonic anemometors used to estimate turbulent intensity near the surface. As a result, we have selected this SL2 for our case study because during its passage all surface sensors, located in both experimental sites, performed well, thus providing good characterization with these observations. At site T0k, measurements of O 3 (model 49i, Thermo Fischer Scientific Inc., Waltham, MA, USA) were performed at a sampling rate of 1 Hz at 40 m above ground in the K34 tower (located at 2.6028 • S-60.2098 • W; see Figure 1). The O 3 analyzer was installed inside a container with a controlled environment, located at ground level, approximately 5 m from the tower. The air was pumped at a rate of 12 L/min, from the inlet to the analyzer using a 9.8 mm diameter teon tube [18].
At the T3 site, measurements of O 3 concentration were performed with an ultraviolet light absorbing gas analyzer (the same model was utilized at the T0k site). The instrument generates a new measurement every 4 s, and was installed at a height of 3.5 m above the ground. To ensure good data quality, the O 3 analyzer was calibrated prior to its installation and regularly reset. The processing of O 3 data was carried out following [21].
The micrometeorological variables measured at T0k (at the K34 tower) were the three wind speed components (u, v, and w) and the virtual temperature for the heights of 1.5, 7.0, 13.5, 18.4, 22.1, 24.5, 31.6, 34.9, 40.4, and 48.2 m, utilizing a sampling rate of 20 Hz (model CSAT3, Campbell Scientific, Inc., Logan, UT, USA). We emphasize that fast response data from a sonic anemometer were measured only at the T0k site.
At site T3, the air temperature was measured at 2.5 m above ground (microwave radiometer, MWRP), as well as profiles of variance and skewness of the vertical wind velocity using a DLPROF (Doppler LiDAR profiles). This device generates temporal and spatial estimates every 1 s, with a vertical resolution of 30 m. For more details regarding the instrumentation used at the T0k and T3 sites and, in general, during the GoAmazon project, see [28,32].
To identify downdraft events, the equivalent potential temperature (θ e ) was calculated from the data of air temperature, water vapor mixing ratio and barometric pressure [33]: with (L v ), the latent heat of vaporization, (c pd ), the heat capacity of dry air, (R d ), the gas constant of dry air, (p 0 ), a reference pressure of 1000 hPa, and (r), the water vapor mixing ratio.

Turbulent Parameters
The data from our case study (April 20) were visually inspected to eliminate faulty data before estimating turbulent parameters. In addition, spike detection and elimination was applied in the fast response data to correct common faults in the time series of turbulent variables [34]. The time interval used to estimate the turbulent statistical moments was 5 min for the night-time period [35] and 30 min for the daytime period [36].
The mean horizontal wind speed was obtained starting from the equation U = √ u 2 + v 2 where u and v are horizontal components of the wind. The standard deviation of the vertical wind, w, was defined as σ w = (w 2 ) 1/2 . We also calculated the turbulent fluxes of momentum −ρw u and sensible heat, H = ρc p T v w . Where ρ is the density of air; c p is the specific heat of the air at constant pressure; T v is the virtual temperature, in addition to higher-order statistical moments, such as the skewness of w, Sk w = w 3 /σ 3 w [37]. The turbulent kinetic energy (TKE) was also estimated using a simplified form of the balance equation, which assumes horizontal homogeneity and absence of advection [37]: where e is the TKE, p is the atmospheric pressure, θ is the mean potential temperature, θ v is the virtual potential temperature, the terms with ( ) represent the fluctuation of the turbulent quantity with respect to its mean according to the "Reynolds mean" procedure [37]. Term I represents the temporal variation of TKE, term II is the production of TKE by vertical wind shear, term III is TKE production due to buoyancy, terms IV and V are the turbulent and pressure correlation of TKE, respectively, and the term VI expresses the dissipation of TKE into heat. In this work, only the terms II, III, and VI of the TKE balance equation (Equation (2)) were used to investigate the SLs role in both TKE production via wind shear (I) and buoyancy (II), and TKE dissipation per unit mass ( ). To calculate the vertical gradient of mean wind velocity in term II, we used the wind speed values provided by the sonic anemometers located at heights between 34.9 and 48.2 m. The vertical momentum flux was calculated at the 40.4 m level. The value of the term III was calculated using data from the sonic anemometer placed at a height of 48.2 m. To obtain the dissipation rate of TKE per unit mass, , the properties of the inertial subdomain of the power spectra established by Kolmogorov [38] were used and the validity of the Taylor Hypothesis [37] was accepted, in order to derive the following equation: where f is the frequency, S u is the power spectral density of the horizontal component of wind speed in the direction of the mean wind, U is the mean wind speed, and α = 0.55 is the Kolmogorov constant [39]. To select only S u values that were within the inertial subrange of the developed turbulence, values of the frequency range 1.0 > f < 2.5 Hz were used [40].

Case Study of a Squall Line
Using 30-min GOES satellite images for 2014, it was possible to observe that between January and December 2014, 13 SL2 events occurred. The months of March and April showed the highest frequency of occurrence of this phenomenon. These months are during the rainy season of the central Amazon region, when the intertropical convergence zone is well established over the N-NE coast of Brazil [41].
After identifying all SL2 events during the selected period, one was selected for the case study (on 20 April). The main reason for using only one case, is that for a better understanding of the SL role in organizing turbulence close to the surface we needed a profile of sonic anemometers working reasonably well during the occurrence of these events. Unfortunately, this only happened in just one case. The SL's events are accompanied by heavy rains, which impairs a good measurement. This SL2 formed around 18 UTC (UTC = local time + 4 h), extending from the Brazilian state of Maranhão ≈ 42 • W to Guyana ≈ 59 • W (Figure 2a, red line). The SL2 propagated inland toward the continent during the early morning. It is also noted that at this moment its intensity was reduced (Figure 2b). Still in Figure 2b, it is possible to observe that there were some intense convective cells (red circles) southwest of SL2. At 0715 UTC (Figure 2c) it is observed that there was a integration of SL2 with these convective cells, making SL2 more intense and well defined. Between 0715 and 0815 UTC (Figure 2d), it is noted that the first branch, which we will call branch 1, arrived at the T3 and T0k sites. This SL2 continues its propagation inland (Figure 2e), with the second branch, called here branch 2, arriving at around 1640 UTC to the both sites ( Figure 2f). This SL2 continued to propagate into the continent and lost its linear definition as of 18 UTC. Based on the methodology of [6], the SL2 event selected for this study had an approximate lifetime of 22 h from its origin until its dissipation. The propagation mean velocity from its origin (in the north coast of South America) around 1800 UTC, until its arrival into the central Amazon region (≈1350 km apart), around 0800 UTC was ≈27 m s −1 , with length of 825 km and width of 200 km. Beginning here, we discuss how this SL2 event occurrence may have modulated the organization of the atmospheric turbulence, thermodynamics, and concentration of superficial levels of O 3 in both sites analyzed.

Surface Values of Ozone, Wind Speed, and Temperature during the Squall Line
During the SL2 passage above both sites, T3 and T0k, there was an increase in wind speed (U) and O 3 concentration and at the same time a decrease in θ e at around 0815 and 1640 UTC (Figure 3a). The latter suggests the occurrence of downdrafts above the sites during these events. It is believed that the first downdraft occurred when SL2-branch 1 was above T3 (at 0745 UTC) and T0k (0815 UTC). The second downdraft probably occurred when SL2 branch 2 was overhead at both sites (1640 UTC). Note that the U values ( Figure 3a) were quite similar for the two sites, oscillating between 0 and 2 m s −1 during the absence of downdrafts and above 4 m s −1 during the occurrence of downdrafts. The maximum values of U were: 5.15 and 6.46 m s −1 during branch1 for T0k and T3 sites, respectively. For the 30 min period before these values were 1.10 and 1.13 m s −1 . During branch 2, the maximum values of U were 4.60 and 6.11 m s −1 for T0k and T3, respectively, for the 30 min period prior, these values were 1.19 and 1.32 m s −1 , that is, the U values increase around five times during the presence of SL2. Figure 3b shows that the first observed downdraft produced a decrease in θ e values of approximately 6 and 10 K, for T3 and T0k, respectively. During the second downdraft the fall of θ e was considerably stronger above T0k (≈16 K) than above T3 (6 K). We believe that the proximity of the T3 site to the very wide Solimões River, played an important role. Being in the proximity of a large water body generates greater thermal inertia (in this case T3), which means that there will be smaller decreases in θ e during downdraft episodes, when compared to the T0k site, which is located in a region of primary forest relatively distant from wide rivers. Figure 3c shows that during the first downdraft, O 3 values increased by 13 and 15 ppbv, for T3 and T0k, respectively. However, it is not possible to observe significant increases in the O 3 concentration close to the surface during the occurrence of the second downdraft. We believe that on this occasion the O 3 concentration gradient between the air transported by the downdraft and that at the observation level was small, resulting in no evident changes in the O 3 concentration at both sites. Attention should be drawn to the fact that after the occurrence of the first downdraft, the T0k site showed a decrease in the O 3 concentration (≈10 ppbv), while in the T3 site the O 3 concentration practically did not change, fluctuating very close to 13 ppbv. This contrasting response can be explained by the fact that the T0k site is located in a forest area, in which the oxidation of sesquiterpenes (emitted by vegetation) represent a sink for O 3 during both daytime and nighttime. According to Freire et al. [27], this oxidation sink would account for about 40% of O 3 loss overnight at T0k.
Several studies have already shown that, during the occurrence of downdrafts, it is possible to observe an increase in the superficial levels of O 3 in the Amazonian forest [15,17,18,21,[42][43][44]. According to Betts et al. [17], this increase is probably due to the action of an air parcel, rich in O 3 , from the middle and upper troposphere, moving via the downdraft towards the ground. This would account for the O 3 concentration increase at the two sites during the occurrence of the first analyzed downdraft. In addition, these downdrafts are also capable of transporting cool dry air from the upper atmosphere to near the surface and, as a consequence, produce a decrease in surface θ e [15,17]. According to Dias-Júnior et al. [21], once this denser air reaches the atmospheric layers near the surface it diverges and propagates horizontally in the form of density current causing increases in the U values. These results are consistent with those of Melo et al. [43] through results of numerical simulations. Figure 4a,b show the temporal variation of the vertical profiles of SK w and σ w respectively, measured above the T3 site. These profiles extend to 1 km above the measurement level with a vertical resolution of 30 m. It is possible to observe that around 0740 and 1640 UTC (these times are delimited by gray dashed lines) SK w and σ w were close to [−3] and [+3 m s −1 ], respectively. These values indicate the presence of strong downdrafts and intense turbulence at the same times that the two branches of SL2 are located exactly above the T3 site.  Figure 5 shows temporal variations of the vertical profiles of σ w now for the region within and immediately above the forest canopy at T0k site. It is evident that the values of σ w increased considerably when branch 1 (around 0815 UTC, Figure 5a) and branch 2 (around 1640 UTC, Figure 5b) of SL2 were passing over the T0k site. Note that intense turbulence is observed only during the presence of the downdrafts, in the absence of an SL and in particular SL2 it is difficult to notice the presence of strong turbulence at depths greater than 0.5 h (where h is the canopy height).

Skewness and Variance of the Vertical Wind Velocity
It is noted that the presence of an SL accompanied by its downdrafts can potentially increase the turbulence intensity in almost the entire observed atmospheric boundary layer ( Figure 4). In addition, above forested areas such downdrafts generate a significant increase in turbulence intensity in almost the entire interior of the forest canopy ( Figure 5). Santana et al. [45], used an extensive dataset from different Amazonian experimental sites and obtained average profiles of various statistical moments associated with the atmospheric turbulence. In their analysis, they demonstrated that the higher density of leaf area of the forest canopy is responsible for filtering practically all the turbulence that would reach the region located below z < 0.5 h. However, here we show that such filtering is insufficient to prevent the turbulent action in the presence of an SL.

Sensible Heat Fluxes
We now pay special attention to the turbulent vertical profiles observed within and immediately above the forest canopy at site T0k. Figure 6 show sensible heat flux (H) values for the periods in which branch 1 (Figure 6a) and branch 2 (Figure 6b) of SL2 passed above the site T0k. It is possible to have an overall view of the temporal variation of the vertical profile of H during the displacement of SL2 above the experimental site. The gray dashed lines indicate the time when SL2 passed above the site.  Several studies have shown that H values above the Amazon rainforest are positive during the daytime and negative at nighttime [39,[45][46][47], corroborating the results found here for situations in the absence of SL2.
In the presence of strong winds, in this case due to the arrival of SL2, (Figure 3a) it was observed that the magnitude of H increased substantially during the nocturnal period. Fitzjarrald and Moore [48] observed H values 10 times higher in conditions of strong winds when compared to the average for all nocturnal periods above their experimental site, the Duke Reserve (02 • 55 51 S, 59 • 58 59 W) in the Amazon very close to Manaus. Dias-Júnior et al. [39] have observed that for situations of high winds and high values of σ w , H corresponds to about 72% of the total turbulent heat flux observed for the entire analyzed nighttime period at the Rebio-Jarú (Amazon) site. Santana et al. [45] observed that the highest values of H occur during the presence of strong winds for both daytime and nighttime conditions.
In the present study, we observe that at the forest-atmosphere interface and within the canopy H oscillates between positive and negative values in the presence of an SL. This leads to an important observation regarding the general formulations of the fluxgradient relationship, as studied in Monin-Obukhov's Similarity Theory [36,37]. One of its assumptions is that the sensible heat flux is approximately constant with the height in the atmospheric surface layer (ASL), a condition that is clearly not verified in the presence of SLs as shown above.

Multiresolution Analysis and Atmospheric Turbulent Regimes
In this section, we focus on the event that occurred during nighttime (branch 1). The data provided by sonic anemometer and sonic thermometer located at the highest available measuring level above the forest canopy were used. At this height, H did not change sign and became strongly negative during the occurrence of the SL2 during the nighttime period (Figure 6a). Therefore, it would be possible to observe the effect of the SL arrival on the H values using multiresolution decomposition [49]. Figure 7 shows the multiresolution w'T'-cospectra for the period of 3 h before the arrival of branch 1 (blue line) and during the passage of branch 1 (red dashed line) at site T0k. The cospectra were calculated using 2 15 data points, which is equivalent to approximately 28 min of data recorded at 20 Hz.
The cospectral peak time scale will be considered as the time scale associated with the highest sensible heat flux [50]. It is observed that the cospectral peak of H for the 3 h period before branch 1 presents a time scale of approximately 13 s. On the other hand, the temporal scale of the cospectral peak during the presence of branch 1 was of approximately 52 s. The larger time scale associated to the spectral peak during the SL2 passage is associated with the presence of larger vortices compared to those that frequently populate the Amazon nocturnal boundary layer in the absence of the SL2. In addition, these vortices present a larger temporal scale and certainly a larger spatial scale being more energetic (greater amplitude of the cospectra), and thus justifying the greater heat exchanges between the forest and the atmosphere during the presence of the SL2 above the site. Figure 8 shows the relationship between σ w and the mean horizontal wind speed (U) obtained through the fast response data measured at the highest level above the T0k site. Each point corresponds to a 5 min average of data, as suggested by [35]. The blue circles and red asterisks correspond to periods of absence and presence of SL2, respectively. The black dashed line indicates the wind speed threshold below which only points belonging to the SL2-absence period were found. Values above this threshold only occurred in the presence of SL2. The value of the wind speed threshold (v L ) splits the available data into two subsets: data with weak and strong turbulence regimes, which have been determined according to the methodology used by Dias-Júnior et al. [39]. In this way, all situations where the wind speed was greater than 2.8 m s −1 were associated with the strong turbulence regime.
In Figure 8, a clear distinction can be noted between the points corresponding to the strong and weak turbulence regimes. The slope of the line corresponding to the weak turbulence regime (blue solid line) is clearly smaller than the slope of the line corresponding to the strong turbulence regime (red line).  It is observed that during the presence of SL2 only the occurrence of the strong turbulence regime (red asterisks) was observed, whereas in the absence of SL2 only the occurrence of weak turbulence regime was verified. A physical explanation for the shift from a weak to strong turbulence regime and vice-versa was provided by Sun et al. [51]. They showed that in the weak turbulence regime, the increase of TKE near the surface is controlled by the energy used to increase the potential turbulent energy (PTE) through the buoyancy flux. Under such conditions, TKE does not increase significantly with U. Still according to Sun et al. [51], in the strong turbulence regime, the potential temperature gradient below a height z is suppressed by the presence of large turbulent vortices, and wind shear generates TKE directly without being consumed by the increase in PTE. This change in the mechanism of PTE consumption in the strong turbulence regime generates a considerable increase of TKE with respect to U. Similar results were found here, as shown on Figures 6, 7 and 9.  Figure 9 shows some TKE balance terms, such as the TKE production rates by wind shear (black line) and buoyancy (blue line), as well as the TKE dissipation rate per unit mass (red dashed line), obtained during the passages of the first (Figure 9a) and the second (Figure 9b) branch of SL2 above the site T0k, respectively. It can be observed that the values of TKE production by wind shear and TKE dissipation are significantly higher during the periods in which SL2 passes above the site, than those in the absence of SL2. This behavior is not observed for the buoyancy TKE production values, which only slightly increase during the passage of SL2. It was also noted that the magnitude of the TKE production by wind shear and TKE dissipation rate were greater during branch 1, than during branch 2 of the SL2 above the T0k site. In other words, the production and dissipation of TKE were stronger during the nighttime period than during daytime period. These results suggest that TKE production by wind shear is the dominant mechanism in the generation of turbulence during the passage of the squall line.

Conclusions
We have studied a squall line (SL2) that originated in the north coast of South America and traveled (≈1350 km) into the central Amazon region. This SL2 caused a significant impact in the vertical exchange processes between the surface and the atmosphere at both, a forested (T0k) and a pasture (T3) site. With the arrival of SL2, the presence of downdrafts were observed, and well characterized by measurements at the surface, mainly: (1) the values of σ w became considerably larger, both in the forested region and in the deforested region; (2) sudden drops in potential equivalent temperature values (θ e ); (3) significant increases in surface ozone levels (O 3 ) and horizontal wind velocity, U, (five times higher), which persisted during the passage of SL2 above the experimental sites. It was also observed that SL2 has the potential to transport O 3 from the high-medium troposphere downwards into the atmospheric boundary layer.
This study represents a first over the Amazon region, showing that the presence of an SL can produce significant changes in the organization of atmospheric turbulence near the surface. It was noted that the turbulent vortices associated with the presence of SL2 have a considerably longer time scale and are more energetic compared to those that are present in the absence of SL2. Sensible heat exchanges during the presence of SL2, both at the forest-atmosphere interface and within the forest canopy, oscillated between positive and negative values during both nighttime and daytime.
It was also observed that the presence of the SL2 above the site modifies the background turbulence regime, causing a transition from a weak to a strong turbulent state. Furthermore, we showed that the TKE production is almost exclusively of mechanical origin. A better understanding of these processes of interaction between forest-atmospheresquall line are fundamental for obtaining more realistic parameterizations to represent the surface processes in numerical models, since the presence of such lines is quite common in the Amazon.

Conflicts of Interest:
The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.