Analysis of the Influence of Deforestation on the Microphysical Parameters of Clouds in the Amazon

: Studies have shown that deforestation can cause changes in energy, moisture, and precipitation flows, with implications for local and regional climate. These studies generally focus on understanding how the hydrological cycle is impacted by deforestation, but few studies have investigated these impacts on cloud microphysics in tropical forest regions. The objective of this study was to quantitatively evaluate the impacts of deforestation on the microphysical parameters of clouds, based on data extracted from active and passive orbital sensors from the TRMM satellite. The study area comprised the state of Rondônia, Brazil. The analyses of the microphysical parameters extracted from the Microwave Imager (TMI) and Precipitation Radar (PR) sensors of the 2A-CLIM and 2A25 products were performed considering a period of 14 years. The parameters analyzed were Rain Water Path (RWP), Ice Water Path (IWP), Surface Precipitation (SP), Freezing Level Height (FH), and Rainfall Type (RT). Land cover type data were extracted from the Project to Monitor Deforestation in the Legal Amazon (PMDA). Our results showed that local deforestation significantly altered the microphysical parameters of the study region. In general, the values of the microphysical parameters of the clouds in the transition areas (locations where forest pixels are neighbors to deforested pixels) were about 5–25% higher compared to forested and deforested areas associated with a higher frequency of episodes of convective rainfall possibly driven by mesoscale circulations. Correspondingly, forested areas had higher rainfall rates compared to deforested areas. Meanwhile, deforested areas had higher amounts for IWP, of around 1–16%, and FH, of around 2–8%, in relation to forested areas. Conversely, the RWP showed a decrease of around 2–20%. These results suggest that the microphysical structure of clouds has different characteristics when related to forested and deforested areas in the Amazon. This is useful for evaluation of simulations of cloud microphysical parameters in numerical models of weather and climate. evaluate the effects of using data Our results indicate that the deforestation produced statistically significant changes in the microphysical parameters of IWP, RWP, and FLH, and in SP rates. Despite this, the transition areas presented the highest values in relation to adjacent areas, which may be associated with mesoscale circulations. Among the transition zones (10 _20 km ZT), these parameters, in general, showed similar behavior, but with different magnitudes, with a tendency to decrease as they advanced towards the interior of adjacent areas (deforested and forested). We also observed that convective rains tended to increase as the deforested area increased. suggest in and from previous


Introduction
The Amazon Rainforest has highly heterogeneous hydroclimatology and is strongly associated with its surrounding atmospheric and hydrological cycle [1]. Covering more than 6 million km 2 , it represents approximately 50% of the world's tropical forest area [2,3], with the greatest biodiversity per area of all tropical forests in the world [4]. The average annual rainfall in the Amazon Basin is around 2300 mm year −1 [5,6], with marked spatial variability. Thus, it is possible to identify up to six regions with well-defined rainfall regimes that are associated with the local conditions of the region and the performance of the main meteorological systems that generate rain. Further details on the six rainfall regimes found in the Amazon are described by [7]. In addition, the Amazon River has an average discharge of 194,100 m 3 s −1 year −1 into the tropical Atlantic Ocean, which is equivalent to one-fifth of the world's freshwater reserves [8].
Changes in the vegetation cover of tropical forests, which occur mainly through the conversion of forest to pasture or cropland, can lead to changes in the thermodynamic characteristics of the lower atmosphere, influencing the regional and global climate [5,[9][10][11][12][13]. Studies have been carried out to assess the impacts of deforestation in the Amazon on the climate from the perspective of regional and global climate models [9,[14][15][16][17][18]. At the same time, these interactions have been investigated through comprehensive observational studies, such as "Anglo-Brazilian Climate Observations in the Amazon" (ABRACOS) [19], "Large scale Biosphere-Atmosphere experiment in Amazonia" (LBA) [20], "Amazonian Aerosol Characterization Experiment" (AMAZE-08) [21], "South American Biomass Burning Analysis" (Sambba) [22], and "Observations and Modeling of the Green Ocean Amazon Experiment" (GoAmazon 2014/5) [23]. These studies, in general, have shown that anthropogenic processes can produce visible impacts on the energy balance of the Earth's surface, including the hydrological and biogeochemical cycle, modifying the physical properties of the Earth's surface.
Regarding the analysis of the direct influence of deforestation on the intensity of precipitation in the Amazon, the study by [24] presented a binomial probability distribution with data from orbital radar and ground surface coverage from the MODIS sensor to investigate the way that rainfall frequency changes as a function of distance to the forest edge. They showed that rain is more frequent in forested areas than in deforested areas. In addition, the authors observed that the transition area between forest and deforestation presented higher occurrence of precipitation, indicating a mechanism of mesoscale circulation between forested and deforested patches. However, these analyses were restricted to elucidating the precipitation component, which has been the subject of several previous reviews. For example, ref. [25], showed that southeastern Rondônia, Brazil, is 25% drier in the less rainy season. Refs. [18,26] show that reductions in rainfall are directly proportional to the extent of deforestation. In contrast, other studies show increased rainfall in deforested areas adjacent to the forest due to mesoscale circulations [27][28][29]. Despite these studies addressing the effects of deforestation on various aspects of rainfall, there are no observational studies that demonstrate the influence of deforestation on the microphysical properties of clouds through comparative analysis.
Once an area of native forest is removed, there is an imbalance in the supply of energy, mass, and moisture flows to the atmosphere [30]. These impacts can be manifested through sudden changes in albedo, which affect evaporation and transpiration, as well as the partitioning of sensible, latent, and soil heat fluxes [30,31]. Analyses with satellite data, field measurements, and numerical models have shown that extensive deforestation in the last three decades in the state of Rondônia in southern Amazonia has altered the mechanism of precipitation, which has shifted from thermal conduction to being dynamically driven [25], in addition to a significant reduction in dry-season evapotranspiration [32].
The Amazon Rainforest, particularly during the dry season, has one of the highest concentrations of aerosol particles in the tropics [33]. These particles emitted by vegetation, and mainly by seasonal biomass burning, are linked to cloud microphysics and rain production processes, which can modify the interactions at the biosphere-atmosphere interface and produce significant effects on the net exchange of CO2, the climate, and human health [23]. Changes in concentrations and type of aerosols imply the formation of cloud drops with implications for the type of precipitation, affecting the life cycle of clouds and the rainfall regime of a region [34][35][36][37][38][39]. In general, a more polluted atmosphere tends to produce more intense precipitation events, with stronger wind gusts and more frequent electrical discharges [27,30,[40][41][42][43][44][45].
These more intense rainfall events observed in deforested regions are generally associated with differences in heat fluxes over forested and deforested areas. As turbulence predominates over deforested areas associated with a higher flow of sensible heat, the formation of mesoscale circulations, also called deforestation breezes [46], occurs, as indicated by a series of studies of different types of vegetation [24,[47][48][49]. These studies suggest that the length scale required for the development of these circulations is approximately between 10 and 100 km, although at low latitudes this spatial scale may be smaller [50].
The study of cloud microphysics is based on the analysis of the particles (liquid or solid-ice crystals) that make up the cloud, that is, physical processes within the cloud [51]. Remote sensing products, via data collected by sensors onboard satellites, provide good spatial coverage and are able to describe finer-scale physical and microphysical processes [1]. With the increasing advances in the estimation of environmental variables through remote sensing, the detection of this association is more easily achieved in comparison with macro-and mesoscale numerical models. While in situ measurements can provide direct and accurate high-resolution local data, satellites provide broader spatial and temporal coverage, allowing refined statistical analysis [52].
In this context, in 1997, the Tropical Rainfall Measurement Mission (TRMM) was launched with the objective of estimating spatialized rainfall data for the tropical region, as well as expanding and improving existing limitations in rainfall estimation techniques, such as those using the lengths of visible and infrared waves [53][54][55]. The TRMM includes, among other sensing instruments, a passive and an active microwave radiometer, which makes it possible to generate integrated information of water and ice hydrometeors, as well as the three-dimensional distribution of the vertical profile of rainfall rates, eliminating, thus, the additional restrictions hitherto existing in the recovery of precipitation [53,56].
The validations (quantitative and physical) of the TRMM precipitation recovery algorithms were carried out through field campaigns in tropical and subtropical regions. [57]. In Amazonia, these data were evaluated through the Brazilian Large-Scale Biosphere-Atmosphere-TRMM-LBA experiment [20,57]. In general, the results showed good agreement with the observed precipitation data, while for the microphysical data the TRMM tends to be underestimated by a factor of 3 when compared to the observed surface radar data being successfully applied in studies in the region Amazon [29,58,59]. However, other studies have observed precipitation underestimations due to the horizontal spacing of the product, which is not capable of representing processes with a spatial scale smaller than 0.25°, and overestimations associated with large-scale systems due to the strong emissivity of clouds. [60][61][62][63].
Some studies have been carried out to evaluate microphysical properties of clouds using TRMM data; for example, in Brazil, [64] used data from the Microwave Imager (TMI) sensor to estimate ice particles in Mesoscale Convective Systems, while [65] used data from Lightning Imaging Sensor (LIS), TMI, and Precipitation Radar (PR) to analyze microphysical characteristics of clouds in precipitating systems with and without lightning. Ref. [66] used rainfall measurement radar (PR), TMI, and visible and infra-red sensor (VIRS) data to verify the impacts of smoke from biomass burning on rainfall processes for Indonesia. Previously, [67] used LIS and PR data to analyze the relationship between ice mass and lightning density for the warm seasons of the Northern and Southern Hemispheres. Other studies with orbital data have also investigated the role of aerosols in modifying the microphysical properties of clouds [52,[67][68][69][70].
As verified by the above analyses, these studies have mainly focused on investigating the effects of land cover change on biogeochemical (aerosols, CO2) and biophysical (energy, moisture flows) flows from precipitating systems, improving our understanding of the effects of deforestation on Amazonian hydrology. To date, no studies have been conducted that investigate the impacts of deforestation on the microphysical properties of clouds through comparative statistical analysis. For example, the difference between total water and ice hydrometeors in precipitating clouds of continental areas, compared to deforested and non-deforested areas of Amazonia, is not known quantitatively. The effects of deforestation on rainfall, energy flows, and aerosols have been extensively studied, as summarized by [71]. Thus, improving this understanding depends on a deep understanding of the changes in microphysical properties on contrasting surfaces in tropical regions such as the Amazon, and this is an open scientific question.
The objective of this study was to investigate for the first time the hypothesis of variance between the content of microphysical parameters on contrasting surfaces of southwestern Amazonia; that is, we aimed to analyze through comparative statistics the manner in which changes in land cover impact the microphysical properties of clouds on a regional scale using observations from the active and passive sensors of the TRMM satellite.

Study Area
The study area corresponds to the state of Rondônia, Brazil (southwest Amazon), situated between latitudes 7°58′S and 13°4″S, and longitudes 59°5″W and 66°4″W ( Figure  1). The state covers an area of approximately 240,398 km 2 , with a population of around 1,600,000 inhabitants (IBGE, 2010). According to the Kôppen classification, the region has an Am climate (Rainy Tropical Climate), with a climatological average air temperature above 18 °C-megathermal [72]. The average annual precipitation in the study region is approximately 2000 mm (32,47). The distribution of precipitation shows strong seasonality: during the months of November to March (rainy season), the precipitation is above 200 mm/month, associated with precipitating systems such as the Intertropical Convergence Zone (ITCZ) and the South Atlantic Convergence Zone (SACZ), among others, while the months from May to September (dry season) are less rainy, with averages below 100 mm, associated with the weakening of these systems and local convection [46,73,74]. Occasionally, the state suffers the action of the phenomenon known locally as "friagem", which is related to the incursions of polar air masses (dry and cold) from southern Brazil [5,75]. Information on land cover characteristics in the study region is provided in the next section.

PMDA Product
The brazilian project to Monitoramento do Desmatamento da Floresta Amazônica Brasileira por Satélite -PRODES, henceforth Monitor Deforestation in the Legal Amazon (PMDA) is a program that monitors clear-cut deforestation in the "Legal Amazon" for the establishment of public policies. In the last three decades, it has produced estimates of the annual rates of deforestation in the region [76,77]. PMDA uses images from LANDSAT satellites, whose products have a spatial resolution of 30 m and a revisit rate of 16 days. Currently, the project uses images from the LANDSAT-8 and Operational Land Imager (OLI) sensor, as well as products from other satellites such as China-Brazil Earth Resources (CBERS 4) and Indian Remote Sensing (IRS-2) [78,79]. With this structure, the minimum area mapped by PMDA is 6.25 hectares, which minimizes the problem of cloud coverage and guarantees interoperability and reliability criteria [80,81]. Made available through the website: http://www.obt.inpe.br/OBT/assuntos/programas/amazonia/prodes [accessed on 10 April 2021].
The land cover classes adopted by PMDA are composed of Water Bodies, Native Forest, Deforested (of the "shallow cut" type, total removal of forest cover), and Non-Forest, comprised of natural savannah, regrowth, urban centers, and wetlands, among others. In 2020, according to PMDA, the forested surface coverage in the state of Rondônia corresponded to 117,856 km 2 (49%), while the deforested area covered 96,093 km 2 (39.9%). The remainder, approximately 26,447 km 2 (11%), was non-forest. Since 2000, the deforested area has increased by around 35,398 km 2 . The evolution of deforestation in the study domain can be seen in Figure 2.

Products of the 2A-CLIM and 2A25 Algorithms
The data referring to the estimates of the microphysical parameters and rainfall were extracted from the products 2A-CLIM (from the Microwave Imager-TMI sensor) and 2A25 (from the Precipitation Radar-PR) from the TRMM satellite, which relate to surface precipitation and vertical profiles of hydrometeors made available by the STORM platform: https://storm.pps.eosdis.nasa.gov/storm/, accessed on 20 March 2019.
The TMI sensor is a five-frequency passive microwave multichannel radiometer; four frequencies have dual polarization (10.65, 19.35, 37.00, and 85.50 GHz) and one has single polarization (21.30 GHz) [54,82]. The PR radar was the first rain radar to be placed in orbit. According to [54], its main objectives are to identify the three-dimensional structure of the vertical distribution of precipitation simultaneously with cloud microphysics, as well as to improve the overall accuracy of the TRMM precipitation measurement by the combined use of data from active (PR) and passive sensors (TMI and VIRS). The PR is a 13.8 GHz active radar, which records energies reflected from atmospheric and surface targets with electronic scanning, with a bandwidth of 215 km pre-boost and 247 km post-boost (Ref. [57]. The publications by [55] and [56] provide more detailed information about the types of microphysical measurements of clouds estimated by the algorithms, as well as technical analysis of the strengths and weaknesses of measurements from the TRMM satellite. From the TMI, the parameters analyzed were Rain Water Path (RWP), Ice Water Path (IWP), and Surface Precipitation (SP), while for the parameters from the PR they were Freezing Height (FLH) and Rain Type (RT). Appendix A shows the spatial variability of the parameters, while Table 1 summarizes some characteristics associated with these products and parameters; details can be found in [83][84][85]. The RWP and IWP parameters consist of the total amount of water and ice, respectively, integrated in the vertical atmospheric column, given in kg/m 2 , while the SP is the instantaneous precipitation rate on the surface, in mm/h [85].
The Rain Type (RT) parameter has three categories: Stratiform, Convective, and Others, which are classified according to two methods called Vertical Profile Method-V_Method, and Horizontal Standard Method-H_Method [86]. According to these authors, the types of rainfall by the two methods are unified and expressed by a three-digit number of the form Type (V_Method, H_Method), whose first digit indicates the main category of the type: 1: Stratiform, 2: Convective, and 3: Others, while the second and third digits indicate the type of rain according to the method used. In that study, the results were presented in percentage normalized by the area of each type of surface considering only the type of stratiform and convective rain.
The Freezing Level Height (FLH) consists of the values of the isotherm heights of 0 °C above mean sea level, given in meters [87]. These data were obtained from the 2A25 product, originally estimated from the 2A23 algorithm with PR data. The FLH is estimated from the knowledge of climatological surface temperature data with a predefined constant lapse rate provided by the National Aeronautics and Space Administration (NASA) [86,88].
FLH is calculated according to Equation (1): where Ts [K] is the interpolated climatological temperature (space/time) at mean sea level. The atmospheric temperature lapse rate is assumed to be 6.0 [K/km]. According to [86,87], sea-level temperature data were provided by NASA shortly before the launch of the TRMM satellite and consisted of 12 months of statistical data that have not been updated since.

Methodology
The data were selected considering the driest quarter in the region: June, July, and August-JJA [46,74], considering the initial interval from 2000 to 2013, which is the common period for all products. The use of data for this interval was necessary because the TRMM satellite started its downward trajectory in 2014, generating uncertainties in its measurements [88,89]. The choice of the dry period defined in this work was justified because it is the ideal period of the year to analyze the effects of deforestation on the microphysical parameters of clouds in the Amazon, since precipitation is generated mainly by local and topographic effects [42] without the influence of large-scale precipitation systems, low-frequency frontal systems, and relatively weak atmospheric instability [27,90]. Other works have analyzed the impacts of Amazon deforestation considering the dry season in the region, which is the period in which the impacts of deforestation are intensified [24,27,29,32,46,91] As mentioned, the period initially considered was from 2000 to 2013. However, we performed an additional analysis using standardized precipitation anomalies and observed that there were years with anomalies (positive and negative) of precipitation with great spatial variability, as can be seen in Figure 3. Standardization was achieved by calculating the difference between the observed value and the mean value for the period (quarter mean-JJA) and dividing the difference by the standard deviation. In this way, and in order to minimize the generated precipitation from the large-scale effects as well as the spatial variability of rainfall, we chose to work with spatially homogeneous years with negative anomalies that might be associated with the occurrence of large-scale natural phenomena such as the El Niño Southern Oscillation (ENSO) in the Tropical Pacific Ocean, and the anomalous warming of the sea surface temperature (SST) in the Tropical Atlantic Ocean [92][93][94][95][96].  To verify the relationship between deforestation and the microphysical characteristics of the clouds, it was necessary to resample the PMDA data (30 m pixel), in order to match the spatial resolution of the microphysical and rainfall data products (5 km pixel). So that all data were analyzed on a spatial resolution grid of 5 km, the nearest-neighbor method was used, due to the processing time being faster when compared to other techniques [97,98]; furthermore, because the extreme subtleties are not lost, the method is ideal for categorical, nominal, and ordinal data [99,100]. Only pixels identified as Forest and Deforested, according to PMDA, were considered in the analyses, while pixels identified as Non-Forest and Water Body, as well as pixels that did not show rain according to TRMM, were disregarded.
Studies have shown that, in the transition between forested and deforested areas, convergence zones and mesoscale circulations can be produced [24,48,50,101]. Although circulations can develop at different spatial scales, according to [102,103], they are more efficient in Transition Zones (TZs) with lengths of 10 and 20 km. TZs, according to [24], are defined as locations where forested pixels are neighbors to deforested pixels; that is, both forested and deforested pixels are considered part of this region since there are one on each side. Therefore, there is a span of at least 10 km where a single forest pixel is next to a single deforested pixel (2 × 5km pixels = 10km). Comparatively, a 20 km ZT extends one more pixel on each side (4 × 5km pixels = 10km). A schematic diagram demonstrating the concept of ZTs is shown in Figure 4c. In this way, for this research, two ZTs were defined, one with 10 km and the other with 20 km in order to verify if there are differences in the influence of deforestation in these transitory strips of surface. This was accomplished by creating a modified land use cover map from the original binary land use cover map, where Deforested pixels bordering Forest pixels were identified with a third coding. Finally, a dynamic land use cover map was produced for all the years of the study. An example of this land use map with the three surfaces, namely: (i) Forest (F), (ii) Transition (T), and (iii) Deforested (D), can be seen in Figure 4. This map of Forest, Deforested, and Transition cells was used as the basis for comparing microphysical cloud and precipitation data in the following analyses. Finally, operations were performed to delineate the study area, as well as to determine comparative statistics, such as mean and standard deviation, and to formulate graphs and tables. Subsequently, the nonparametric Wilcoxon-Mann-Whitney (WMW) test for means was applied, adopting a significance level of 5%, in which the null hypothesis was that the means between different groups were the same. This test was chosen because it does not consider the structure of the data (normality) and is an alternative to the Student t-test when the samples do not follow a normal distribution [104]. Details on the test can be found in [105,106]. For the application of the referred test, the averages of two types of surfaces were taken at a time. For example, first the differences in the averages between the Transition areas were compared with the averages of the Forest areas, then the averages of the Transition areas with those Deforested, and finally the averages of the Deforested areas with those of Forest. Finally, three groups of means were compared by the referred test.
To quantify the Relative Effect (RE) of deforestation on the microphysical parameters (Table 1), the absolute differences between the following combinations were calculated: Transition (T)-Forest (F); Transition (T)-Deforested (D) and Deforested (D)-Forest (F), where the results were divided by the surfaces F, D, and F, respectively, according to Equation (2) (2) where RE corresponds to the average relative effect given in percentage, P1 denotes the average values of the parameters of Deforested (D) and/or Transition (T) regions, and P2 denotes the average values of the parameters of Forest or Deforested regions. It is worth mentioning that the calculations were performed considering the average of all points for each type of surface.
Additionally, the degree of influence of deforestation on microphysical and rainfall parameters was verified, considering three precipitation classes. For this purpose, rainfall data were classified according to intensity using the percentile technique [107]. First, 50%, 5%, and 1% of the data were selected considering all points for each surface type (deforested, forest, and transition), so that the necessary comparisons could be made. Intense Precipitation Event (IPE) rates were defined as values above or equal to the 95th percentile (hereinafter, P95) and above the 99th percentile (hereinafter, P99), while the values considering the intervals between 45 and 55% were called P50, representing the central or average value of the data series.
In the dry season in the Amazon region, rainfall occurs mainly over elevated areas [42]. Thus, relief data were also obtained from the Shuttle Radar Topographic Mission (SRTM) project (http://www.dsr.inpe.br/topodata/, accessed on 1 August 2022), accessed on 8 January 2020, to complement the analyses, in order to filter out the influence of topography on the microphysical parameters. According to the analyses performed (not shown), surfaces with topographies above 721 m can influence the amount of precipitation and cloud microphysics in the study area. Thus, all regions (pixels) with topography above this threshold were disregarded.
In spatial data, in general, the attributes or phenomena analyzed present spatial autocorrelation or spatial correlation. Thus, it is important to verify and quantify the degree of spatial autocorrelation of the microphysical parameters used in this research (Table 1), since most statistical methods are based on the assumption that the samples are independent [105]. For this purpose, the Moran Global Index (MGI) was applied [108]. The MGI has been successfully applied in several areas of knowledge [109][110][111][112]. The MGI consists of a general measure that tests the degree of spatial autocorrelation existing in a dataset, in which the null hypothesis states that the analyzed attribute is randomly distributed among the characteristics of the study area [113]. It can range from −1 to +1, where positive values indicate direct/positive autocorrelation and negative values the opposite; a null value indicates that the data are randomly distributed [113,114]. The results showed that the MGI ranged between 0.017 and 0.038 (p-value < 0.05), which implies a weak positive spatial correlation; in other words, the spatial distribution of the microphysical processes that occur over a pixel, for example deforested, tends to be similar to the neighboring deforested pixel and different from distant pixels such as forest or transition, which is what we expected. Therefore, proceeding with a standard approach of comparative statistical analysis is acceptable and the results are valid.

Influence of Amazon Deforestation on Cloud Microphysics
The results of the influence of deforestation on precipitation and cloud microphysi-  Tables 2-5. A three-way comparison was performed, where for each event the parameter values of each of the three surface coverages were compared and statistics were recorded.

Surface Precipitation (SP)
For the first parameter, surface precipitation (SP), it can be seen from Figure 5 that the rainfall rates are different between the analyzed areas. On average, SP is higher in the transition areas compared to the other two surfaces, regardless of their intensity or extent of the TZ (Figure 5a,b). For P50 events, the SP rates of deforested areas tend to be lower compared to forest, around 22% for the 10 km TZ and ~15% for the 20 km TZ (Table 2). However, for P95 and P99 class events, the opposite was observed, with increases between ~2 and ~10% of SP rates in deforested areas compared to forested areas, with the exception of P99 of the 10 km TZ, which showed a ~5% decrease (p-value > 5%). In summary, rainfall rates were higher in the transition area compared to the other two areas, regardless of rainfall class and TZ extension. In deforested regions in relation to forested ones, rates decreased or increased depending on the observed rainfall class (intensity).

Ice Water Path (IWP)
For the microphysical parameter Ice Water Path (IWP), we observed that, similarly to the behavior of SP, the highest IWP values occurred in transition areas compared to forest and deforested areas ( Figure 6). The exception is for the P50 (ZT of 20 km), where the transition and deforested areas presented the same values, and showed a decrease of ~0.5% (p-value > 5%) in relation to deforested areas, but an increase of 15% (p-value < 5%) in relation to forest areas (Figure 6b, Table 3). However, the value of IWP in deforested areas tended to be higher in relation to forest by around ~1 to 35%, depending on the rainfall range and the TZ considered (Table 3). These results suggest that, as one moves away from forested areas towards deforested areas, IWP values tend to increase in relation to forested areas. However, the transition areas continued to have higher values in relation to the other two (Figure 6a,b). These results reinforce the presence of the mesoscale circulation discussed earlier.

Rain Water Path (RWP)
The Rain Water Path parameter (RWP) was similar to SP and to IWP; that is, the transition areas presented the highest values of RWP in relation to the values of forest and deforested areas, regardless of the rain category and the extension of the TZ (Figure  7a,b). However, when observing the deforested areas in relation to the forest, both reductions and increases were observed. For example, reductions were observed for P50 and P99 of around 1 to 20%, but the reductions observed for P99 of deforested areas in relation to forested areas were not statistically significant (Table 4). Increases were observed only for P95 class events in both TZs.

Freezing Level Height (FLH)
For the Freezing Height (FLH) parameter, when the transition and deforested areas were compared with the forest, increases between ~2 and ~9% (p-value < 5%), respectively, were observed, as shown by Figure 8 and Table 5. While the transition areas in relation to deforested areas oscillated between increases and decreases depending on the category of rainfall event and the extent of the TZ considered, these differences were not statistically significant (p-value > 5%) ( Table 5).

Rain Type (RT)
Regarding the Rain Type (RT) variable, we analyzed the frequency of stratiform and convective rainfall according to the type of surface. According to Figure 9, stratiform rainfall is more frequent in forested areas than in transition and deforested areas for both TZs (Figure 9a,b). However, convective rains are more frequent in transition regions followed by deforested ones for 10 km TZ (Figure 9a), while for 20 km TZ the deforested areas presented higher frequencies of convective rains in relation to the other two ( Figure  9b).

Discussion
The observed reductions for P50 of the deforested areas in relation to the forests can be explained by changes in the surface properties of the vegetation. Deforested areas are characterized by having higher albedo, lower roughness length, and lower rooting depth, resulting in lower evapotranspiration and latent heat flux and higher sensible heat flux. [31,33,71]. These results corroborate the pioneering studies by [14,115], who showed that changes in soil surface characteristics cause a reduction in evapotranspiration rates, and therefore induce a decrease in rainfall in Amazonia. However, as the events become intense and rare (P95 and P99), relative increases were observed in deforested areas in relation to forest, which may be associated with increased cloud cover and rain over deforested areas due to mesoscale circulations [91,103]. These results are consistent with previous studies that showed that rainfall in deforested areas is higher than in forested areas due to increased surface heating in degraded areas, which induces upward air movements, increasing moisture generating rain in deforested areas [12,27,28].
Regarding the higher SP rates observed in transition areas compared to adjacent areas, these results can be explained by the vegetation heterogeneity that produces prominent gradients in surface flows due to land use change, which in turn can trigger thermally induced mesoscale circulations [101,116,117]. This heterogeneous environment can create zones of convergence. Physically, this corresponds to deforestation-induced mesoscale circulation, similar to sea breezes carrying moist air, with convective potential energy available from the forested region to the deforested area, where it rises and precipitates [24]. Consequently, rainfall may increase, associated with mesoscale circulation as a result of initial deforestation.
This observation corroborates what has already been shown by other studies: there is some type of mesoscale circulation between the limits of forest areas and deforested areas [24,103,118]. According to [71], these circulations are responsible for bringing the forest air to the deepest boundary layer on the warmer surface, favoring the beginning of convection on the deforested side. Moreover, according to these authors, the air over the forest tends to have a higher equivalent potential temperature, associated with a lower albedo and a shallower boundary layer. By comparison, the air over the deforested region tends to have less convective inhibition, due to warmer temperatures and a deeper boundary layer. In addition, the mesoscale convergence contributes to increasing the vertical velocity which, together with the humidity and convective energy from the forest region, contributes to the formation of convective precipitation clouds.
A series of studies including meso-and large-scale models in Amazonia [102,119], Africa [103,116], and Australia [120], among others [121,122], confirms the presence of breezes induced by the heterogeneity of vegetation. Our results suggest that these mesoscale circulations are more effective when rainfall events tend to be stronger, since increased SP values were observed in deforested areas 20 km away from the edge of the transition area.
As observed, the SP of deforested areas in relation to forest areas may decrease or increase depending on the intensity of the event. Alternatively, only increases in the IWP of the deforested areas were observed in relation to the forest, regardless of the intensity of the rainfall event and the TZ (Table 3). Therefore, a plausible explanation for the successive increase in the amount of IWP in transition and deforested areas in relation to forest is that changes caused by deforestation in energy flows between the surface and the atmosphere impact cloud patterns, the Atmospheric Boundary Layer (ABL) thickness, and the intensity and duration of precipitation [71]. That is, these changes observed in the IWP are due to the action of convective clouds that are frequent in deforested areas [27]. This type of cloud is characterized by producing severe storms with great vertical development, with a greater amount of ice and a thicker ABL, influencing the duration and intensity of precipitation. Convective rains, for example, associated with a more unstable atmosphere with high concentrations of anthropogenic aerosols, are characterized by being more intense, with strong winds and more abundant electrical discharges. These favor the transport of droplets to higher levels of the atmosphere, where convection can be renewed and the amount of ice increased [35,42,123].
In the Amazonian dry season, aerosol concentrations are high due to biomass burning, influencing the formation, duration, and amount of rain [34,39,43]. For example [124] reported that the increase in biomass in the atmosphere decreases the size of cloud droplets, in turn slowing their growth and shifting the initiation of rain to higher altitudes. This configuration increases the ice production associated with supercooled drops in the mixed-phase region. In addition, there will be increased electrification and lightning due to stronger updrafts and increased latent heat at higher levels [123]. This may also explain the higher values of the amount of ice and water in the atmospheric column observed in the deforested and transition regions ( Figure 6 and Table 3).
The observed increases and decreases in FLH may be associated with the different behavior of the thermodynamic atmosphere structure caused by changes in energy flows due to deforestation. The surfaces of deforested regions, for example, are warmer and drier [35]. Accordingly, it can be considered by the boundary layer theory that the lapse rates over the forest and over the deforested regions are different; that is, the FLH value in the deforested areas will be higher than in the forest area [27,71]. These changes in thermodynamic behavior in deforested areas associated with changes in energy flows may explain the changes observed in these areas. However, the influence of other factors contributing to these changes in FLH cannot be ruled out. Ref. [87], studying the functional relationships of latent heat, latitude, and surface temperature between FLHs, found that latitude was the main controlling factor of FLH rather than temperature and latent heat. In addition, they found that the FLH depended on the time of year and responded differently between inland and coastal sites, whereas on an interannual basis, ENSO variations can play a critical role in modulating FLH [125,126] During the Amazon dry season, convective rainfall over deforested areas is characterized by convective air currents controlled mainly by diurnal convective heating [124]. This thermodynamics raises the cloud base, favoring supersaturation and droplet concentration associated with improved conditional instability [127,128]. Our results agree with the findings of [41] and [27], who observed that deforested regions have higher frequency of storm convection.
Precipitation events and deforestation in the study region of southeastern Brazilian Amazon are part of a very dynamic system. Nonetheless, this study was subject to limitations, such as the data being representative only of the afternoon and dry-season rain events. In addition, events that were not stimulated locally (meso-and large-scale precipitating systems, occurrences of cold weather, etc.) were not fully filtered in this study. However, we empirically show, despite these sources of noise, that cloud microphysical parameters are, in fact, more abundant in transition areas compared to distant forested and deforested areas due to the higher frequency of convective rainfall episodes, possibly driven by mesoscale circulations associated with surface heterogeneity. This is consistent with the conceptual model shown in Figure 10. Figure 10. Conceptual model of the effects of Amazon deforestation on cloud microphysical parameters. Surface heterogeneity due to deforestation creates horizontal temperature gradients producing precipitation events in a convective manner, possibly originating from mesoscale circulations leading to an increase/decrease in cloud microphysical parameters.

Conclusions
The main objective of this study was to evaluate the effects of the current level of deforestation in the state of Rondônia, southwestern Brazilian Amazon, on rainfall rate and cloud microphysics parameters, using data extracted from orbital sensors, for the driest quarter in the region (JJA). The results described here point to possible changes in these parameters associated with regional deforestation.
Our results indicate that the deforestation produced statistically significant changes in the microphysical parameters of IWP, RWP, and FLH, and in SP rates. Despite this, the transition areas presented the highest values in relation to adjacent areas, which may be associated with mesoscale circulations. Among the transition zones (10 _20 km ZT), these parameters, in general, showed similar behavior, but with different magnitudes, with a tendency to decrease as they advanced towards the interior of adjacent areas (deforested and forested). We also observed that convective rains tended to increase as the deforested area increased.
These results suggest that, in addition to rainfall rates, the microphysical structure of clouds has different characteristics when related to forested and deforested areas, in agreement with the changes observed in cloud patterns and precipitation amounts from previous research.
The potential impact of mesoscale circulation associated with surface heterogeneity on rainfall and cloud microphysics is critical for a proper understanding of the long-term impacts of deforestation, and for improving parameterizations in regional and global numerical models.
Deforestation in the Amazon region is a reality. If it continues, the changes in rainfall patterns will continue to evolve, and may even further harm the region's hydrological and energy regime. The understanding of specific components of Amazonian hydrology (precipitation, evapotranspiration, aerosols, etc.) is hampered by geographic limitations and relies on numerical models rather than empirical data. Thus, further studies in this area of knowledge are recommended, such as investigating the impacts of Amazon deforestation on cloud microphysics for other regions of the Amazon, performing comparative analyses with data from mesoscale numerical models, and evaluating performance with orbital data.
This study is a first examination of the regional impacts of Amazonian deforestation on cloud microphysics. The approach carried out in this work provides unprecedented and general results with analyses that had not previously been carried out. However, our research was restricted to afternoon cycle rainfall events and to dry years (negatively homogeneous rainfall anomalies). A posteriori, more in-depth analysis should be carried out to investigate these differences in the daily cycle for rainfall events in rainy years (positively homogeneous rainfall anomalies) and for normal years considering a longer period.