Vapor plumes in a tropical wet forest Spotting the invisible evaporation

. Forest evaporation exports a vast amount of water vapor from land ecosystems into the atmosphere. Meanwhile, evaporation during rain events is neglected or considered of minor importance in dense ecosystems. Air convection moves the water vapor upwards leading to the formation of large invisible vapor plumes, while the identiﬁcation of visible vapor plumes has not yet been studied. This work describes the formation process of vapor plumes in a tropical wet forest as evidence of evaporation processes happening during rain events. In the dry season of 2018 at La Selva Biological Station (LSBS) in Costa Rica it was possible to spot visible vapor plumes within the forest canopy. The combination of time-lapse videos at the canopy top with conventional meteorological measurements along the canopy proﬁle allowed us to identify the driver conditions required for this process to happen. This phenomenon happened only during rain events. Visible vapor plumes during the daytime occurred when the following three conditions are accomplished: presence of precipitation ( P ), air convection, and a lifting condensation level value smaller than 100 m at 43 m ( z lcl . 43


Introduction
Forest cover in tropical regions is endangered by deforestation (Curtis et al., 2018;Rosa et al., 2016), compromising the evaporation flux from land. Forest evaporation is a mixture of water vapor originating from water intercepted on plant surfaces, soil water, and plant transpiration (Roberts, 1999;Savenije, 2004;Shuttleworth, 1993). Forest evaporation is considered of major importance as a regional and local cooling system (Ellison et al., 2017) as a result of its capacity to recycle the atmospheric moisture at different timescales (van der Ent and Savenije, 2011). The water vapor originating from evaporation at the surface is horizontally transported in the atmosphere by advection (Lavers et al., 2015;Strong et al., 2007), where forest presence at continental scale induces the "biotic pump mechanism" that favors the maintenance of similar precipitation amounts between inland and coastal environments (Makarieva and Gorshkov, 2007;Makarieva et al., 2013a). Meanwhile, vertical transport is linked to wind shear (Chen et al., 2015) and convection (Trzeciak et al., 2017) that, in large ecosystems, influence the formation of convective clouds at the top of the atmospheric boundary layer (Fuentes et al., 2016;Manoli et al., 2016). This process plays an important role in the formation of precipitation in tropical basins (Adams et al., 2011;van der Ent and Savenije, 2011) because of the contribution of water vapor originating from local evaporation (Brubaker et al., 1993).
Evaporation is usually neglected or considered of minor importance during rain events in dense forest ecosystems (Klaassen et al., 1998). This is because during rainfall the vapor pressure deficit is close to zero (Bosveld and Bouten, 2003;Loescher et al., 2005;Mallick et al., 2016), reducing the atmospheric water demand and stopping the transpiration process (Gotsch et al., 2014). However, the increment of evaporation with the size of rain events suggests that evaporation also occurs during the events and not only after-Published by Copernicus Publications on behalf of the European Geosciences Union.
wards (Allen et al., 2020). This has been evidenced by discrepancies found between modeled and measured evaporation rates in tropical forests (Schellekens et al., 2000). When it rains, part of the precipitation is intercepted and evaporated directly to the atmosphere (David et al., 2006) even when vapor pressure deficit and available radiation are low (Lankreijer et al., 1999). Under high-humidity conditions, a portion of the precipitation can evaporate after a raindrop splashes on the canopy or the forest floor. This process is known as "splash-droplet evaporation" (Dunin et al., 1988;Dunkerley, 2009;Murakami, 2006) and is based on the principle that raindrop size increases with rain intensity. Consequently, when larger drops hit the surface (e.g., ground, leaves, and branches), this allows the formation of smaller rain droplets that can be easily evaporated after the splash. This process has been identified as the main source of evaporation that explains the difference between intercepted water and measured evaporation in studies carried out in banana plants (Bassette and Bussière, 2008) and eucalyptus plantations (Dunin et al., 1988).
Forest evaporation produces coherent structures of water vapor called plumes, cells, or rolls (Couvreux et al., 2010). Plumes of water vapor have been identified above forest ecosystems during the daytime with a high-resolution scanning Raman lidar technique (Cooper et al., 2006;Kao et al., 2000). These plumes reached heights above the canopy up to 100 m, depicting their importance as water vapor providers at local scale. This phenomenon has been studied in astrophysics (Berg et al., 2016;Sparks et al., 2019), vulcanology (Kern et al., 2017;Sioris et al., 2016), and regional and global meteorology (Herman et al., 2017;Knoche and Kunstmann, 2013;Wang, 2003;Wright et al., 2017). However, to the best of our knowledge, little attention has been paid to the small events observed during rain events. Additionally, Couvreux et al. (2010) highlighted the lack of sampling techniques able to characterize the occurrence of these plumes close to the surface. Visible vapor plumes are classified as ascending clouds formed by clusters of tiny particles of water in liquid form (Spellman, 2012). This characteristic makes it difficult to measure them with sophisticated systems based on 3D wind components (e.g., eddy-covariance systems) that are developed to measure water in gas form (Foken et al., 2012a). These types of measurements are sensitive to rainy and high-humidity conditions (Camuffo, 2019;Foken et al., 2012b;Kelton and Bricout, 1964;Moncrieff et al., 2005;Mauder and Zeeman, 2018;Peters et al., 1998), making it difficult to use them to identify the occurrence of visible vapor plumes in forested ecosystems. This mismatch between measurement systems and target phenomena underlines the need to identify the conditions under which visible vapor plumes are formed. This type of constraint requires an innovative data analysis approach, which is the focus of this paper. This work aims (1) to test an innovative approach to link visual information and conventional meteorological data describing a local hydrological phenomenon. It also aims (2) to identify the meteorological conditions when visible vapor plumes are present in a tropical wet forest and it tries (3) to explain the processes involved in the formation of these plumes. The data analysis is based on conventional meteorological data vertically distributed along the forest canopy layer and time-lapse videos during daytime conditions.

Study site
The monitoring was carried out at La Selva Biological Station (LSBS) on the Caribbean lowlands of Costa Rica (10 • 26 00 N, 83 • 59 00 W). This station registered a mean annual precipitation of 4351 mm yr −1 , a mean annual temperature of 26.3 • C, and a mean daily temperature difference of 9.5 • C. A short dry season occurs at LSBS between February and April every year, and it is characterized by a reduction in precipitation without the vegetation experiencing a soil water deficit Lieberman and Lieberman, 1987;Loescher et al., 2005). LSBS is covered by a matrix of old-growth and secondary forests, small forest plantations, and experimental permanent plots with mixed tree species (Fig. 1). All instrumentation was placed within the area of 1.0 ha of the major research infrastructure plot (MRI plot), located within an old-growth forest on the upper terrace of the Puerto Viejo river . The MRI plot is situated in the upper section of a small hill facing southwest towards an affluent of the Puerto Viejo river. The soil is classified as Andic Humidotropept with a clay and organic matter content of 35 % and 23 %, respectively (Sollins et al., 1994). The tree density in 2017 was 371 trees ha −1 of individuals with a tree diameter bigger than 10 cm. The palm Welfia regia H. Wendl and the tree Pentaclethra macroloba (Willd.) Kuntze are the most abundant species with 56 and 43 trees ha −1 , respectively. The average leaf area index (LAI) in 2005 was 3.56 m 2 m −2 (Tang et al., 2012). The plot is located within a stable forest plot in terms of changes in canopy height and tree biomass fixation (Dubayah et al., 2010).

Experimental design
The monitoring was carried out on the MRI plot in the highest tower (43 m), which is located within a depression of the forest canopy (Fig. 1). Along the vertical axis of the tower the air temperature ( • C) and relative humidity (%) were measured with HOBO ® smart sensors (part code: S-THB-M008). The sensors were located at 2, 8, and 43 m height, placed at a distance of 1.5 m from the tower, and protected with a radiation shield (HOBO ® part code: RS-3) of 10 cm diameter. The use of radiation shields together with conventional air temperature sensors allows for keeping a mean absolute error during the daytime in warm tropical ecosystems below 0.3 • C (da Cunha, 2015; Terando et al., 2017). Also, the shelter provided by the forest canopy for the measurements carried out at 2 m helps one record similar temperatures to the surrounding near-surface environment (Lundquist and Huggett, 2008). The measurement of minimum air temperatures or nighttime temperatures does not require the cover of the radiation shield to keep low biases (< 0.5 • C) on the mean air temperature due to the reduced or total absence of solar radiation (Terando et al., 2017). At the highest point of the tower, the precipitation (mm min −1 ) was recorded with a Davis ® rain gauge. Soil temperature ( • C) was measured in two different locations at 5 and 15 cm depth with a soil temperature sensor (HOBO ® part code: TMC20-HD). Soil moisture ( , m 3 m −3 ) was measured at the same locations as soil temperature at 5 cm depth with an ECH 2 O ® EC sensor. Soil temperature was recorded with a four-channel data logger (HOBO ® part code: U12-008) and the other sensors with a USB Micro Station (HOBO ® part code: H21-USB). Meteorological data collected along the tower, as well as soil temperature data were recorded with averages of 1 and 5 min, respectively. All data were summarized in 5 min time intervals for the analysis. A Bushnell ® Natureview ® Essential HD camera (12 megapixels) was installed at the top of the tower facing northwest.

Monitoring period
All environmental variables were monitored between 24 January 2018 and 26 March 2018. The camera was installed to collect photographs above the canopy between 21 and 25 March 2018. The photographs were set to be collected continuously from 05:00 to 18:30 local time (UTC−6). However, the light conditions determined which images were selected for analysis (see Appendix A). These pictures were used to determine when the vapor plumes were visible at the MRI plot. The photographs were classified into three conditions ( Fig. 2): -Clear view includes all the pictures with clear and cloudy sky where the canopy is clearly visible and there is neither mist nor plumes present ( Fig. 2a and b).
-Mist and fog include the presence of a homogeneous blurry view of the canopy. The blurriness of each picture varies depending on the humidity conditions. Special care was taken to prevent the erroneous classification of photographs affected by a fogged-up lens. This category is subsequently called "mist" (Fig. 2c).

Data analysis
Data processing and analysis was performed with the opensource software R (R Core Team, 2017). All temperatures were converted from K to • C. Superficial soil temperature (T s.0 , • C) was estimated with Eq. (1) (Holmes et al., 2008). This equation describes the diurnal variations of soil temperature as sine waves depending on the 24 h moving averages of soil temperature at 5 cm depth (T s.5 , • C). The daily amplitude of air temperature (T A , • C) is defined as the difference between T s.5 and the air temperature at 2 m (T 2 m ). The oscillations are determined by the damping depth (ν, m), which is calculated with Eq.
(2). The depth difference between the T s.0 and T s.5 is defined as z b (m). The sine pattern depends on the angular frequency (ω, s −1 ), time (t) in s, and φ (-) as a constant for phase change. Equation (3) is used to determine ω with τ (s) as the wave period. Equation (2) calculates ν with the soil thermal diffusivity (η, m 2 s −1 ) and ω. Equation (4) (Nakshabandi and Kohnke, 1965) is used to determine η, where ρ s is the soil bulk density of 0.76 g cm −3 (Sollins et al., 1994) for the experimental plot, c s is the specific heat for clay soils (837.36 W kg −1 • C −1 ), and k is the soil thermal conductivity of 1.58 W m −1 • C −1 (Pielke, 2013). The last two parameters were chosen according to the soil water conditions during the monitoring period, which was close to soil field capacity (see Appendix B).
The virtual potential temperature (θ v , • C) of the air was calculated to take into account the variation in the adiabatic lapse rate due to changes in pressure (Barr et al., 1994;Stull, 1988Stull, , 2017. For saturated (cloudy) air conditions Eq. (5) calculates the θ v based on the water vapor mixing ratio (ψ s ) of the saturated air, the liquid water mixing ratio (ψ L ), and the virtual temperature (θ ). The parameters ψ s and ψ L were determined with Eqs. (7) and (8), respectively. These equations require one to know the mass of the liquid water in the air (m liq.air ), the mass of the water vapor in the air (m vap.air ), and the mass of the dry air (m dry.air ). Due to the lack of instrumentation to estimate the mass of liquid water in the air, we used a fixed value of 0.05 g m −3 . This value corresponds to the liquid water content (LWC) in the air reported by Thompson (2007) for continental fog events. The selection of this value was based on (1) the similarity between the vapor plumes and fog and (2) the fact that both types of events occur close to the ground surface. The variables m vap.air and m dry.air were determined using the saturation and actual vapor pressures of the air (Stull, 2017). The virtual temperature was estimated with Eq. (6), where d is the dry adiabatic lapse rate near the surface (0.0098 • C m −1 ), z is the height above the ground in meters, and T z is the air temperature at the same heights.
Convection can be identified by evaluating the temperature gradient ( θ v z ) due to the absence of wind profile measurements to determine the atmospheric stability parameter along the tower. Values of θ v z > 0 are linked to stable stratification; meanwhile, θ v z < 0 shows an unstable stratification (Stull, 2017), which will drive convection.
The condensation of vapor close to the forest canopy can be identified by calculating the lifting condensation level (z lcl ) in meters with Eq. (9). This equation determines the elevation at which a parcel of air condensates allowing the formation of clouds. This equation uses the difference between air temperature (T z ) and dew point temperature (T dew.z ) at one specific height (z), divided by the difference between d and the dew point temperature lapse rate ( dew ) (Stull, 2017).
An estimation of the evaporation during the monitored period was retrieved from Jiménez- . This data set is used only as a reference for the evaporation process during the monitoring period on the same site. This is because this quantification has limitations accomplishing the Monin-Obukhov similarity theory (MOST) for complex terrains (Breedt et al., 2018). So, it is based only on the vertical transport of water vapor, neglecting the advected energy of the forest canopy.

Results and discussion
The monitoring period experienced a diurnal variation in air temperature along the vertical profile of the canopy, with a temperature difference of more than 10 • C at 43 m and less than 7 • C at 2 m height (Fig. 3). The highest temperatures were registered at 43 m height reaching more than 30 • C, decreasing in magnitude towards the forest floor. These peak temperatures were recorded around noon with differences up to 5 • C between the air temperatures at 43 and 2 m height. The T s.0 oscillates between 20.7 and 25.4 • C. The amplitude of the oscillation increased with the sunniest days, but the daily difference did not exceed 4 • C. The maximum value was 0.47 m 3 m −3 during the heavy rains, almost reaching the saturation point for clay soils of 0.50 m 3 m −3 (Saxton and Rawls, 2006). The minimum was recorded after the driest period just before the rains on 24 March 2018 (0.42 m 3 m −3 ), getting close to soil field capacity for clay soils (Saxton and Rawls, 2006). Evaporation always occurred during daytime on all sampling days (Fig. 3). During the four sunny days, the evaporation was larger than 5 mm d −1 , with a contribution of more than 1.0 mm d −1 from 8 m height and no more than 0.7 mm d −1 from 2 m height . In contrast, on 24 March 2018 the continuous rains summed up to 58.7 mm d −1 and the evaporation was estimated as 1.8 mm d −1 at 43 m and only 0.2 mm d −1 at 2 m height (Table 1).
During the visual monitoring with the field camera, clear view conditions were predominant along four days (Fig. 3). These days were characterized by sunny conditions with temperatures above the 25 • C, no large rain events, and a decreasing trend in soil moisture. These days were characterized by cumulus clouds crossing the sky above the forest canopy in the daytime. Any water vapor ascending from the forest canopy needed to reach a height of more than 100 m to form visible vapor plumes (Fig. 3). Also, on 24 March 2018 it was possible to identify three short periods with   (Table A1). These mist events were linked with superficial soil temperatures higher than 2 • C with respect to air temperature. Finally, the vapor plumes were visible only during rainy conditions on 24 March 2018 (videos available at Jiménez-Rodríguez et al., 2019b). Soil temperatures during this day were warmer than the air column along the forest canopy (Fig. 3).
Evaporation during sunny days provided the conditions to form vapor plumes as those described by Cooper et al. (2006) and Kao et al. (2000). The evaporation peaks during these days occurred around noon, registering a z lcl higher than 500 m (Fig. 3), which is the height required to form visible clouds. This is the reason why is not possible to see the vapor rising from the surface. The vapor plumes were visible on the day with continuous precipitation (24 March 2018). On this day, the z lcl dropped beneath 100 m because during rain events the θ v of all the air columns dropped quickly. This drop kept the θ v beneath the superficial soil temperature, allowing a localized convection event. This convection process forced the evaporated water to move upwards, forming buoyant clouds close to the forest surface. The evaporation during rain events is the result of the splash-droplet evaporation process (Murakami, 2006;Dunkerley, 2009), which can provide water vapor as a consequence of the fragmentation of raindrops when hitting the surface.
Energy convection plays an important role in forest ecosystems during the nighttime (Bosveld et al., 1999). This is a consequence of the mass transport capacity of the intermittent nocturnal convective fluxes (Cooper et al., 2006). The convection process is forced by the ground heat flux (Jacobs et al., 1994), which is enhanced by the larger soil moisture in clay soil, which increases the soil heat capacity (Abu-Hamdeh, 2003). A coupled canopy system enables sensible heat and water vapor transport from the soil to the atmosphere just above the canopy layer (Göckede et al., 2007). This facilitates the generation of the convection process, allowing the ascending warm air to cool down at the canopy top and condensate, forming the visible water vapor plumes. The condensation releases heat (Goosse, 2015), driving the convection. Vapor plumes are always present as a consequence of the moisture exchange between the surface and the atmosphere (Lawford, 1996), where evaporation from land cover with enough water supply provides the required air moisture (Kao et al., 2000). However, the conditions needed to form a visible buoyant cloud close to the surface require a big difference in air temperature over height. The temperature gradient at 43, 8, and 2 m is negative during plumes and mist conditions; meanwhile clear view conditions have a larger range with more positive values (see Appendix C).
The visible vapor plumes can be spotted on the canopy depressions surrounding the tower (Fig. 1). These depressions are characterized by a low leaf area index and shorter canopy height, which translates into areas with low potential to produce transpiration during rain events. This implies that the main source of water vapor is linked to water evaporated from wet surfaces and soil evaporation, while transpiration may contribute to a lesser extent. Visible vapor plumes are the result of the condensation of water vapor rising from a warmer surface. When a column of warm humid air reaches the dew point temperature, the water vapor condensates around aerosols in the air allowing the formation of clouds (Stull, 2017). In this regard, there are different sources of aerosols at LSBS. One source is linked to wind carrying aerosols from nearby agricultural land uses (Loescher et al., 2004). A second source is linked to convective rains that characterize the dry season at LSBS. Additionally, the presence of bioparticles cannot be discarded (e.g., airborne bacteria, fungi, pollen, and plant fragments) as a source of aerosols from the forests (Huffman et al., 2013;Valsan et al., 2015). The high intensity rains may induce a bioparticles burst from the forest canopy. These bioparticles have been reported in Australia (Bigg et al., 2015), India (Valsan et al., 2015), Mexico (Rodriguez-Gomez et al., 2020), and the Amazon (Pöschl et al., 2010). Also, convective rains transport a portion of the required aerosols for the condensation process from the free troposphere into the boundary layer and later they form clouds . Meanwhile, the splashdroplet evaporation process (Murakami, 2006) provides the main source of water vapor after raindrops hit the canopy and soil surfaces. As plumes are not stagnant and continue moving upwards thanks to air convection, the water vapor is removed from the understory towards higher altitudes. The water condensation at the canopy level drastically reduced the volume of water vapor due to the phase change (Makarieva et al., 2013b). This allowed the ambient air to remain unsaturated allowing the continuous provision of water vapor by the splash-droplet evaporation process continuously more water vapor.
Cloud formation usually happens high above the surface boundary layer where the forest canopy is located, but the available information of cloud formation close to the forest canopy is scarce. The temperature gradient ( θ v z ) at 43, 8, and 2 m is negative during plumes and mist conditions, while clear view conditions have a larger range with more positive gradients. Lifting condensation level is a key element that allowed us to differentiate between plumes and mist conditions (see Appendix C). The combination of variables such as z lcl , θ v z , and P allows one to identify the formation of vapor plumes in tropical wet forests (Fig. 4). The z lcl is the height in the atmosphere at which a parcel of moist air becomes saturated if it experiences a forced ascent (Stull, 2017). It provides an estimate of the height at which the clouds can be formed. The temperature gradient is an indicator of how easily a parcel of air can be lifted (Spellman, 2012) and can be used as a proxy of the atmospheric stability. During unstable atmospheric conditions ( θ v z < 0) it is easier for the air parcels to move upwards than under stable conditions ( θ v z > 0). Finally, precipitation saturates the air column and provides the water vapor after the splash-droplet evaporation process on the canopy and forest floor surfaces.
During the full monitoring period at LSBS, only 1.4 % of our study period accomplished the conditions required for the formation of visible vapor plumes (precipitation, z lcl < 100 m, and 0 > θ v z > −1). These conditions differ from those needed to form mist. In a tropical wet forest in Costa Rica, fog and mist formation happens before sunrise (Allen et al., 1972). However, fog does not involve the upward con-  vective flux needed for vapor plumes, while mist is affected by this upward convective flux but without rain (Stull, 2017). Vapor plumes are buoyant cloud formations with an identifiable shape (Spellman, 2012), which is the main characteristic that allows differentiation between fog and mist events. Mist and fog are formed by microscopic water droplets floating in the air, which can reduce the visibility to less than 1 km in the case of fog or a lesser extent with mist (Spellman, 2012). This paper described the formation of visible vapor plumes based on photographs as a visual indication of a process that is usually invisible to the human eye. The occurrence of this phenomenon under rainy conditions makes it difficult to quantify its contribution to the forest evaporation with current measuring techniques. Vapor plumes occurring during rainy days compromise the performance of more sophisticated instruments that are highly sensitive to rain or mist conditions (Centre for Atmospheric Science, 2020; Mauder and Zeeman, 2018). Instruments such as sonic anemometers (e.g., CSAT3, CSAT3B) and open-path CO 2 /H 2 O analyzers (e.g., LI-7500) are strongly affected by high humidity and rainfall (Campbell Scientific Inc., 2017Foken et al., 2012a;LI-COR, 2016;Moncrieff et al., 2005). The presence of rain causes departures from the measurements, increasing the sonic speed (Camuffo, 2019;Kelton and Bricout, 1964;Peters et al., 1998) or blocking the face of the transducers (Campbell Scientific Inc., 2017), causing a frequency loss during rain events (Zhang et al., 2016). The eddy-covariance technique is considered the standard measurement for determining atmospheric fluxes; however, it is dependent on fully turbulent transport over a homogeneous surface (Foken et al., 2012a). This means that the localized nature of the visible vapor plumes makes measuring them very susceptible to sensor placement, complicating its monitoring using eddycovariance systems located high above the canopy. Additionally, measuring devices based on 3D wind components (e.g., eddy-covariance systems) are developed to measure water in gas form (Foken et al., 2012a) and are not intended to measure visible vapor plumes that are ascending clusters of tiny water particles (Spellman, 2012).
The description of the formation process of visible vapor plumes provides a first step towards understanding this phenomenon within forest hydrology. This description helps identify when this phenomenon occurs, allowing one to screen existing data sets in other tropical research sites to analyze its frequency of occurrence. However, it is important to test if the conditions required to form visible vapor plumes are the same in other latitudes and ecosystems. Also, new developments in air temperature monitoring techniques, such as distributed temperature sensing (Euser et al., 2014;Heusinkveld et al., 2020;Izett et al., 2019;Schilperoort et al., 2018) or thermal infrared imagery (Costa et al., 2019;Egea et al., 2017;Lapidot et al., 2019;Nieto et al., 2019), may contribute to accurately quantify the contribution of visible vapor plumes as local recyclers of forest evaporation. These methods are suitable alternatives to eddy-covariance systems that are sensitive to rainy conditions when visible vapor plumes occur.
Understanding the formation process is a prerequisite before the quantification of such a complex process. Further studies aiming to analyze the occurrence of visible vapor plumes in more detail will need to consider the conditions that give origin to this phenomenon: air convection, precipitation presence, and lifting condensation level at the top of the canopy. While the quantification of its contribution to the hydrological cycle have to overcome the limitations of current measuring techniques. The identification of air convection should be based on direct measurements of ground heat flux (e.g., soil heat flux plates), a more detailed air temperature profile along the forest canopy (e.g., using distributed temperature sensing), and multiple wind speed measurements along the canopy profile (e.g., at least one per canopy layer plus one above the canopy). This set of measurements will help identify air convection and advection within the canopy structure. The wind measurements should be carried pairing sonic and cup anemometers at the same locations, allowing one to overcome the limitations of sonic anemometers during rain events (Mauder and Zeeman, 2018), when liquid water covers the ultrasonic transducers.
The evaporation contribution to the local hydrological cycle by visible vapor plumes requires a detailed quantification of the latent heat flux (ρλE) above and below the canopy. The use of net radiometers at different heights (same locations as wind speed measurements) will complement the detailed air temperature and wind speed measurements. It is important to underline that some experimental sites worldwide accomplish the equipment requirements mentioned above (FLUXNET, 2020), opening the opportunity to reanalyze their data sets towards the identification of the conditions needed for the formation of visible vapor plumes. Also, these sites provide an opportunity to quantify the bias that eddycovariance systems have due to the existence of this phenomenon. Direct measurements of atmospheric water (gas and liquid phase) can be achieved with closed-path gas analyzers (e.g., LI-7000DS-LI-COR, EC155-Campbell Sci., FMA-Los Gatos Research), allowing one to determine the total water content in the air. These measurements will benefit from combining high-resolution infrared images from above and below the canopy, allowing one to study the spatial distribution of the phenomenon. These images will provide information under day and night conditions, helping one to identify the splash-droplet evaporation process at canopy and ground level when the view field is focused towards specific locations of the forest canopy. Finally, further research can search for the detailed source of vapor with the implementation of the direct measurement of water stable isotopes using mass spectrometers or cavity output spectroscopy. This type of research can provide more insights into the effect of vapor plumes on the microclimate of forest ecosystems. Moreover, the occurrence of this phenomenon in other vegetation types may be addressed to understand the main drivers and the role played in local hydrological systems.

Conclusions
Visual monitoring captured the formation of visible vapor plumes close to the surface boundary layer of a tropical wet forest during rainy conditions. These visible plumes are the visual evidence of evaporation processes happening during rain events, where the splash-droplet evaporation process provides the required water vapor to form visible vapor plumes. This water vapor is part of the intercepted water evaporated from the forest floor and plant surfaces since transpiration is likely reduced by the low vapor pressure deficit but not stopped. It is raised by air convection, driven by warm soil temperatures, finally condensing close to the forest canopy due to the drop in the virtual potential air temperature along the forest air column. Consequently, this phenomenon can be identified in tropical wet forest when precipitation occurs, the lifting condensation level at 43 m height (z lcl ) is lower than 100 m, and the temperature gradient ( θ v z ) at 2 m height is between 0 and −1 • C m −1 . Contrary to the vapor plumes, mist appears when no precipitation occurs (P = 0), z lcl at 43 m is less than 10 m, and θ v z is less than −1 • C m −1 . This work also brings attention to the forest evaporation role during rain events, where little information is available. The exploratory nature of this work opened new research opportunities aiming to improve the setup to monitor this phenomenon and provide a further accurate quantification of the contribution within the local hydrology.
Appendix A: Time-lapse videos detailed information Author contributions. The project conceptualization and funding acquisition was carried out by CDJR and MCG. The investigation and data curation was carried out by AdPGA and CDJR. Data analysis was performed by BS and CDJR. Finally, the project administration, writing of the original manuscript, and data visualization was carried out by CDJR with input from all coauthors.
Competing interests. The authors declare that they have no conflict of interest.
Acknowledgements. This work was carried out with a fellowship from the Organization for Tropical Studies (Glaxo Centroamerica Fellowship-Fund 502), with the aid of a scholarship from PINN-MICITT Costa Rica (contract: PED-032-2015-1), and the aid of grant 863.15.022 from The Netherlands Organization for Scientific Research (NWO). Also, NASA's funding NNX12AN43H and 80NSSC18K0708 provided the Leaf Area Index data sets. Special thanks to Bernal Matarrita, Orlando Vargas, Wagner López, Danilo Brenes, Diego Dierick, Enrique Castro, and Marisol Luna for their help and advice at the research station, to all the staff of the OTS for their willingness to support our project, and to Shigeki Murakami and one anonymous reviewer who helped to improve the paper. Review statement. This paper was edited by Laurent Pfister and reviewed by Shigeki Murakami and one anonymous referee.