Characterization of Firebrands Released From Different Burning Tree Species

The number, dimensions, and initial velocity of the firebrands released from burning Quercus suber, Eucalyptus globulus, Quercus robur, and Pinus pinaster trees were analyzed in laboratory experiments using a particle image velocimetry system. Additionally, the flame height, tree mass decay, vertical flow velocity, and temperature at the top of the trees were measured during the experiments. The relationship between the various parameters was analyzed and a good connection was found. The specimens burnt were mostly young trees, so large particles (e.g., pine cones, thick trunk barks, branches) were not included in this study as they were not present. Actually, the firebrands produced in the laboratory tests, mainly burning leaves, had a cross-sectional area of <1,600 mm2, having the potential to cause short distance spotting (up to tens of meters). Quercus trees are often considered to have a lower fire risk than eucalyptus or pine trees. However, in this study, Quercus suber and Quercus robur were the species that produced more firebrands, both in terms of number and total volume. The tests with Quercus suber were the only ones using specimens from an adult tree, confirming the great importance of the age of trees in the propensity to release firebrands. The results obtained with Quercus robur confirmed the high tendency of this species to originate spot fires at a short distance. Thus, these results are of great relevance to afforestation plans and to evaluating the risk of the presence of these species in wildland–urban interface areas.


INTRODUCTION
Spotting is one of the most relevant and common mechanisms of fire spread in wildfires, especially in extreme burning conditions. Its importance derives from the increased difficulty in firefighting caused by spotting; the serious consequences in terms of accidents, infrastructures, and other economic losses caused; and the great difficulty with regard to the prediction of this phenomenon. Moreover, spotting can dramatically increase the fire spread rate, being in some circumstances the main factor affecting it (Rothermel, 1983). Indeed, despite the several studies performed on this topic, the prediction of spotting is still very basic, which clearly jeopardizes the firefighting operations and threatens fire safety.
The analysis of the spot fires mechanism is frequently divided into four stages : (1) release of firebrands, (2) lofting of the firebrands in the thermal plume, (3) downwind transport of the firebrands, and (4) ignition of the fuel recipient where the firebrands land. Most of the work seeking the determination of the maximum spotting distance is carried out in the second and third phases (e.g., Tarifa et al., 1967;Albini, 1979;Ellis, 2000;Hall et al., 2015) considering that firebrands, after leaving the thermal plume, travel at a velocity defined by the terminal velocity and the wind speed. The determination of the probability of ignition of landing firebrands was also the focus of some experimental studies. Statistical analysis based on tests with different combinations of firebrand types and fuel bed recipients (e.g., Ganteaume et al., 2009;Viegas et al., 2014;Ellis, 2015) or in a shower of firebrands artificially generated that hit building façades or other relevant structures (e.g., Manzello et al., 2011;Manzello and Suzuki, 2014;Oliveira et al., 2014) was produced.
The work performed on the release of firebrands is less frequent, and consequently, there is a great lack of knowledge on this phase of the process. The release of firebrands depends on several factors, the most important ones being the fire intensity, the airflow direction, velocity and turbulence, the characteristics of the source of the firebrands (e.g., tree type, species, age, etc.), and the characteristics of the firebrands (type, shape and dimensions, terminal velocity, moisture content, rate of combustion, etc.). Research on the firebrand generation can be divided into two areas, namely (1) firebrands generated by burning structures, which was initially investigated by Waterman (1969), paving the way for more recent research such as Suzuki et al. (2012) or Foote et al. (2011); and (2) firebrands generated by vegetation, which is the main topic of this work and is an even less investigated problem.
Six relevant works addressing this topic are as follows: (1) Gould et al. (2008) performed four field experiments in plots with the size of 200 × 200 m and, among other analyses, studied the firebrand generation and spotting behavior resulting from a fire front spreading in a eucalyptus forest.
Besides the successful validation of the model developed by Ellis (2000) to predicting maximum spotting distance of fibrous firebrands, these authors confirmed the importance of the fire intensity, and consequently of the convection column, in the release of the firebrands. It was also stated that younger fuels, including young trees, produce fewer firebrands than mature specimens. (2 and 3) Manzello et al. (2007) analyzed the size and mass distribution of firebrands generated from burning Douglasfir (Pseudotsuga menziesii) trees with heights of 2.6 and 5.2 m by collecting the firebrands in pans methodically distributed on the ground. The average firebrand size reported was 3 and 4 mm in diameter and 40 and 53 mm in length for the 2.6m tall and 5.2-m tall trees, respectively. Later, Manzello et al. (2008) carried out similar tests with 3.6-m high Korean pine trees, which did not result in a significant number of firebrands produced. However, the dominant firebrands collected had a cylindrical shape with a diameter of 5 mm and a length of 40 mm, on average. (4) Almeida et al. (2014) used a methodology based on particle image velocimetry (PIV) to analyze the firebrands lofted during the burning of pieces of the bark of eucalyptus trees in several scenarios with different arrangements of the barks tested. The methodology then used was quite similar to the one described herein. (5) Tohidi et al. (2015) presented a mechanical firebrand break-off model for determining the formation of cylindrical firebrands produced during the burning of coniferous trees, stating the importance of the limb break-off process on the release of the particles. For validation, the laboratory data presented in Manzello et al. (2007) and Manzello et al. (2008), previously mentioned, were used. (6) Filkov et al. (2017) studied the firebrand production in a prescribed fire in a pine forest by following a methodology described in Houssaimi et al. (2015). Around 80% of the firebrands collected in pans distributed on the ground, mostly from bark fragments (≈70%) and twigs (≈30%), had a crosssectional area between 50 and 200 mm 2 .
The present work intends to contribute to mitigating these knowledge gaps. The main focus is studying the release of firebrands by typical tree species in the Mediterranean area, namely cork trees (Quercus suber), oak trees (Quercus robur), pine trees (Pinus pinaster), and eucalyptus trees (Eucalyptus globulus). These species were chosen among the other tree species with major representation in the Mediterranean forests (San-Miguel-Ayanz et al., 2016) for their role in the spotting mechanisms. The eucalyptus species is by far the most associated one to spotting in several fire events such as the Black Saturday Australian Fires (Cruz et al., 2012) or the fires on October 15, 2017 that occurred in Central Portugal (Viegas et al., 2019). Spot fires driven from pine species also have some references, such as the 2010 Mount Carmel Forest Fire (Israel) described by Kutiel (2012). The role of the oak and cork trees is not as recognized. However, its potential has been well-established in several fire events, for example: (a) Fire of Cicouro-Portugal  where an episode of oak leaves spotting drastically increased the rate of fire spread, surprising a group of firefighters, causing two fatalities; (b) Fire of Pedrógão Grande of June 2017, in Central Portugal , where the control of the fire in Regadas was lost due to spotting of cork oak leaves across a 10-m wide road.
The main aim of this work is to understand the relative role of the mentioned tree species in the short-distance spotting process. Spotting by large firebrands is not addressed in this work since the intense airflow convection column required to naturally produce such particles (Gould et al., 2008) was not achieved in the laboratory experiments. On the other hand, the trees tested did not have large particles such as thick barks or cones available to burn and potentially produce medium-to long-distance spotting, but, given the large number of particles produced, they have a great potential to produce a large number of spot fires. The most common firebrands observed in the laboratory tests carried out were burning leaves released from the trees. Many reports (e.g., Viegas et al., 2013Viegas et al., , 2019Ribeiro et al., 2014) on large forest fires mention spotting by leaves as an essential mechanism TABLE 1 | Main data related to the experimental conditions (m 0 -initial mass of the trees; m f -residual mass of the trees after the burning tests; H-tree height; FMC-fuel moisture content in wet basis; T-air temperature; RH-air relative humidity). of the fire spread. However, the majority of the research on spot fires addresses rigid artificial shapes (like spheres, cylinders, and disks) or natural shapes (cones, barks, and twigs) and usually neglect other potential firebrands like leaves. Clements (1977) is one of the rare exceptions where the terminal velocity of leaves from 17 different species were determined. This author found terminal velocities values of 1.54, 1.57, and 1.80 m·s −1 for the leaves of Quercus species (Q. phellos, Q. nigra, and Q. falcata, respectively) and 2.90, 4.00, and 4.14 m·s −1 for Pinus trees needles (P. echinate, P. elliottii, and P. taeda, respectively).

METHODOLOGY General Conditions and Characterization of Fuels
The experimental work presented herein was carried out at the Forest Fires Research Laboratory (LEIF: Laboratório de Estudos sobre Incêndios Florestais), located in Lousã (Coimbra-Portugal), of Association for the Development of Industrial Aerodynamics (ADAI: Associação para o Desenvolvimento da Aerodinâmica Industrial). The young trees were collected from forests close to the LEIF in the early hours on the day of the experiments. The samples of Quercus suber were collected 2 days before the experiments. Quercus suber is a protected species in Portugal, but a license (reference 53521/2014/DCNF-C/DLAP) for cutting a mature specimen that was about to fall on a road was obtained. Due to the impossibility of cutting young specimens, the branches of this mature tree were used as samples of Quercus suber. For simplicity, the authors will mention each Quercus suber branch as a tree and use this differentiation during the discussion of results.
The tests were performed on a single day. The 10 previous days were normally sunny with an average temperature of 19.7 • C (from 12 to 31 • C), an average wind velocity of 6.2 km·h −1 (from 0 to 22 km.h −1 ), and an accumulated precipitation of 1.0 mm (data from www.wunderground.com). So, the trees collected presented a moisture content value that allowed a sustained burning. On the day of the experiments, the fuels were stored in the laboratory and protected from the sun with the indoor environmental conditions  Table 1 shows the moisture content values of the different trees used. The moisture content was determined immediately before each burning test by using a moisture analyzer (AnD ML-50 0.1%) using samples composed of leaves and branches (<3 mm in diameter) in a percentage distribution representative of each species. A brief characterization of the structure of the trees used is presented in Table 2. This determination was carried out by detaching the leaves and branches of one representative tree (i.e., roughly medium dimensions) from each species. The main trunk of the trees was not considered in the analysis. The images of the trees used in the experiments are presented in Figure 1.

Experimental Apparatus and Procedures
The experimental apparatus can be seen in Figure 2. The tests were carried out in a burning trees platform with an area of 1.5 × 1.5 m and a height of 0.15 m, supported by three load cells with 0.1 g of precision each, connected to a computer for automatic mass recording.
The burning of the trees was only sustainable through the addition of a complementary fuel under the canopy. Since the use of bushes or other wild fuels as an ignition fuel would produce firebrands that would disturb the final analysis, denatured alcohol (DA) was used as a complementary fuel.  Thus, in each burning test, a container with a diameter of 0.75 cm filled with 0.35 L of DA (ethanol 90% v/v) was placed under the tree. In order to check the eventual effect of the DA on the final results and to correct the mass loss, the authors carried out some experiments with a similar methodology but with no trees, thus exclusively burning 0.35 L of DA. As expected, these tests did not produce any firebrand and allowed for the determination of the constant mass-loss rate of the burning alcohol, which was found to be equal to 2.8 g/s.
Using a supporting structure, an S-Pitot tube (6 mm of internal diameter) and a thermocouple (1-mm thick) were placed over the platform, 1.8 m above the ground. This equipment was connected to a computer in order to record the vertical velocity and temperature of the convective airflow. Consistently, the mass, the temperature, and the pressure were registered at a frequency of 2 Hz during the experiments.
For each test, the tree specimen was placed vertically on the platform after measuring the fuel moisture content (FMC) of the samples of leaves and branches < 3 mm in diameter. Ignition was achieved by dropping a match on the DA, which immediately started to burn, spreading the flames to the tree canopy. The DA continued to burn with a constant intensity/flame height until after the test was over.

Characterization of the Particles Released Using the PIV System
For determining the particles released by the trees, a system of PIV was used. The software used was DynamicStudio, version 15, from Dantec. A Speedsense 1,040 camera and a Frame Grabber 1,031 from National Instruments (National Instruments Corp., Switzerland) 1 were used.
The PIV device was located 5.8 m away from the tree and the platform and 3.3 m above the ground, so as to minimize parallax errors, and horizontally oriented to the smoke plume. The image selected for analysis had an area of 0.7 × 0.7 m, thus capturing images in the height range of 2.95-3.65 m (Figure 3). The height of the PIV capturing area was chosen considering an expected maximum flame height of about 2.5 m, according to previous tests performed by Almeida et al. (2014), thus avoiding the interference of the flame in the analysis of the PIV images. Since there were no predominant horizontal airflows, the flames observed tended to be vertical. Therefore, the flame height coincided with the flame length and consisted of the vertical distance between the DA level and the top of the flame. Eventual detached flames were not considered.
As can be seen in Figure 3, the PIV image-capturing system targeted the central axis of the smoke column at a height of 3.3 m. In the middle of the smoke column, the firebrands have a strong vertical upward component (Figure 4). After being lofted, the firebrands, commonly extinct or less glowing, start a descending movement out of the smoke column and thus were not able to be captured by the PIV system. Thus, almost all the particles captured in the PIV images are glowing and being lofted in the convective column of smoke.
The authors assumed a radial symmetry of the tree and of the characteristics of the smoke plume. The PIV images were acquired at a frequency of 143 Hz, in double frame, using the Shadow Processing tool to identify the firebrands (Figure 4). The domain of interest in the PIV images was the analysis area of size 0.7 × 0.7 m (square 2 of Figure 3). This methodology was previously tested with good results, allowing for an accurate particle velocity determination through the comparison of two consecutive frames, thus reducing the likelihood of firebrands being incorrectly identified. Since the brightness of the incandescent particles highlights them in the images taken, the process of tracking is facilitated without the need for complementary illumination systems (laser or other). Thus, not only the number and area of the firebrands but also their velocity can be accurately determined using the double frame in the Cross-Correlation tool.
Due to the limitation of the computer storage capacity, and because of the large number of images being captured, the PIV images were acquired during 44 s in each test. Thus, the capture of images by the PIV system was initiated only when the first release of firebrands was observed. This process was started typically by an initial increase in the release rate of the firebrands and, after reaching a maximum value at time t M , the number of firebrands produced started to decrease. In order to increase the accuracy of the method, the authors used the data of the time interval t M ± 5 s, which was designated as the "10 s-peak period." Thus, the authors ensured that the comparison between the various burning tests was made at the same stage of the process around the maximum release period. Since the period of 10 s provided in this method is limited in duration and consequently in sample size, when appropriate, the entire PIV image acquisition period of 44 s was used in the analysis of the results.
In Figure 5, on the right of the burning tree, it is possible to see a vertical pole with 25-cm markers that were used as a scale to estimate the flame height, using frames captured from the video recordings. The same scale was used to calibrate the PIV images.
The firebrands released were divided into five classes based on average cross-sectional area A av : 10< A av ≤50 mm 2 ; 50< A av ≤100 mm 2 ; 100< A av ≤500 mm 2 ; 500< A av ≤1,000 mm 2 ; A av >1,000 mm 2 . Firebrands with A av ≤10 mm 2 were not considered in this analysis since their potential to cause new ignitions is very low.
It was observed that the majority of firebrands released were the leaves of the burning trees. Thus, since the authors did not find in the literature typical values of the terminal velocity for the leaves of the species used in these tests, a set of tests with that purpose was performed using a methodology very similar to that used by Clements (1977). Ten leaves of each species, randomly collected, were dropped from a height of 9 m in the absence of wind. The time t fall that each leaf took to fall in the last 5 m, when the falling velocity is practically constant, was found using video footage. The terminal velocity V t was determined by the quotient between the 5 m distance and the time t fall elapsed on this final 5 m. The resulting V t values ( Table 3) were not significantly different from those found by Clements (1977) for the species of the same taxonomic genre.

Number and Cross-Sectional Area of the Firebrands Released
The total number of firebrands released during the capture period of the PIV images is presented in Figure 6. In order to allow for the extension of the use of these data in future works, the number of particles released is denoted as those released per unit of area (N·m −2 ). Once again, due to data storage limitations of the system, the capture of PIV images was started after fire ignition, when the operators considered the burning process to already be in sustained development. Therefore, the first point    of the graph in Figure 6 does not initiate at t = 0 s, when the ignition occurred.
The distribution of the total number of firebrands released during the capture period of the PIV images, grouped in classes based on their cross-sectional area A av , is presented in Figure 7.
As previously mentioned, the analysis was limited to the 10 speak period, allowing for a fair comparison among the different species used.
Each pair of tests on the same species showed a very similar tendency in terms of the number of firebrands released, except  in the oak trees tests, where, for test OAK01, the number of firebrands counted was much higher than that for OAK02. By observing the OAK02 test video, it was possible to notice that during the capture period of the PIV images, this tree burned mostly on one side of the canopy unlike the other tests, where the burning was more comprehensive. Therefore, the release of firebrands in the OAK02 test may have been affected. Regarding the distribution of firebrands by the classes based on the crosssectional area A av , the results show a satisfactory reproducibility in the tendency obtained for each pair of tests, even for both the oak trees tests.
It is possible to state that the maximum number of firebrands was released by the cork trees and by the oak trees (excluding OAK02), followed by the pine trees and finally the eucalyptus trees. On the other side, the most represented area classes A av in terms of the number of firebrands released are [10, 50] mm 2 and [100, 500] mm 2 . This figure also shows that Cork tree tests have a higher tendency to produce larger firebrands, while tests with the other species drove mainly toward the production of smaller particles within the range of 10<A av ≤50 mm 2 . It is worth noting that the cork trees were collected 2 days before the tests while the other trees were cut on the day of the experiments. Although the storage of the cork trees occurred under controlled conditions and these trees did not show visible degradation, the decay of FMC (

Temperature and Vertical Airflow Velocity
In the graphs in Figure 8, the values of temperature T and vertical airflow velocity V f , measured at a height of 1.8 m and the values of the velocity of the particles V p determined using the PIV images at the heights of 2.95-3.65 m are represented as a function of the time elapsed since the ignition up to 75 saver. The number of firebrands released is represented in the same figure to facilitate the establishment of relationships. The velocity and temperature values of the CORK02 test were not obtained as the Pitot tube did not work properly during this experiment. The sampling of the values of the velocity of the firebrands is limited to the 44s acquisition time period for the PIV images for all tests. This period is marked in the figure by the two vertical solid lines. The dotted vertical lines in the graphs limit the 10 s-peak period of the particles released.
As can be verified, the values of the velocity of the firebrands do not differ considerably from the values of airflow velocity registered at the same time, following an approximate trend as the number of firebrands released. It would be expected that the V p values were consistently smaller than the V f values because of gravity and drag effects. It is worth noting that the values of V f at a height of 1.8 m were not measured exactly at the same position as that of the firebrands (2.95-3.65 m high). Besides, the values of V f correspond only to the vertical component of the flow velocity, while V p is the modulus of the two-dimensional velocity of the particle and is limited to the PIV image-capturing area. This could be the reason why occasionally the values of V p are larger than V f . The average velocity of the particles in the smoke plume should be given by the difference between the airflow velocity and the terminal velocity of the particle V t (Equation 1).
During the burning tests, the authors observed that most of the firebrands released were leaves or parts of leaves. Some twigs were also released but in much smaller quantities. Therefore, the authors used the terminal velocity values measured for the leaves of each species (Table 3) to make a comparison with the results obtained with the PIV system. The results given in Table 4 show that the values of V t obtained in the PIV tests are of the same order of magnitude as those measured in the free fall tests. The small discrepancies can be justified by the differences in the relative measurement positions of V f and V f , as explained above.
In the case of OAK2 test, a negative value of V t was obtained possibly because the tree did not burn with radial symmetry. Thus, the central axis of the smoke plume did not coincide with the central vertical axis of the tree where the Pitot tube was located. Therefore, the V f value determined was not the highest flow value in the smoke plume, leading to apparent negative V p values.

Mass Decay
In order to estimate the mass loss of each tree during the burning test, the authors considered the mass loss of the DA, which was also placed in a pan on the weighing platform. A reference test was made using only the alcohol container under the same conditions as those of the burning tests. Assuming that the burning conditions of the alcohol in the pan were not affected by the presence of the tree, in all tests, the authors estimated the mass m of the tree by subtracting to the total mass m T the mass m DA of the alcohol, which decreased linearly with a rate of 2.8g·s −1 (r 2 = 0.99).
In order to make a fair comparison of the mass loss among all tests, the mass values were normalized using the initial value of the tree mass m 0 and its value after the burning test m f , as given in Equation (2) for the determination of the relative mass α. The relative mass α as a function of the elapsed time for all tests performed is presented in Figure 9.
During the 10 s-peak period presented in Figure 9, the relative mass decay curves can be fitted to an exponential function as indicated by Equation (3) and as shown in Figure 10.
In Equation (4), the relative mass decay coefficient "a" allows for the comparison among the several tests. As can be seen in Table 5, the OAK02 test, which was not a conventional one, and the eucalyptus trees tests were those in which the relative mass consumption was slower. The relative mass loss rateα, which is used to estimate the relative mass decay per unit of time, is defined by Equation (4). The up-flow velocity V f and the velocity of the firebrands Vp are average values determined for the 10 s-peak period with greater release of firebrands. The specific tests were previously described in Table 3. TABLE 5 | Exponential relative mass decay coefficient "a" (s −1 ) and correlation coefficient achieved in the fitting of the relative mass α of the trees by the elapsed time t during the 10 s-peak period (Figure 10), according to Equation (4 The variation in relative mass loss rateα during the 44 s of the capture period of the PIV images is shown in Figure 11. This representation of the relative mass loss rateα indicates a faster mass consumption for cork tests, OAK01, and pine tests.

Correlation Between Parameters
The mass and the volume of the firebrands released were not determined since only one PIV camera was used and consequently only a two-dimensional analysis was carried out. However, a rough association between the cross-sectional area can be made with the volume and mass of the particles. The accumulated values of the total cross-sectional area A ac (Equation 5) of the firebrands released since the beginning of image acquisition up to a certain moment t were related with the relative mass α of the trees, as presented in Figure 12.
As expected, the accumulated total cross-sectional area A ac of the firebrands released increases with the decrease in tree mass (Figure 12). The tree mass decay is dependent on three parameters: (1) the mass transferred to the firebrands released, (2) the mass lost in the water volatilization process, and (3) the mass lost to the formation of volatile gases. The relationship given in Figure 12 deals with that between the mass transferred and the firebrands released. Making the same relationship for the 10 s-peak period (Figure 13), the slope m of the tendency lines formed by each data set (Table 6) can be analyzed. Larger values of the slope m mean less mass reduction dependency on the mass associated with the A ac of the firebrands released.
Using the same relationship in Figure 13, the exponential coefficient k obtained from Equation (6) was determined as presented in Table 6.
The results show that the mass loss of burning eucalyptus trees is less driven by the production of firebrands than for other trees. This is consistent with Figure 7, where it is shown that the small firebrands were the most represented in the EUC experiments. On the opposite side, the variation of α vs. Ac for the cork trees was very strong as denoted by the lowest k value, which is consistent with the larger number of firebrands released (Figure 6, above) and with the lowest values of moisture content ( Table 1).
When equations 3 and 6 are combined using the common parameter α, Equation (7) is obtained.
In Figure 14, the evolution of A ac as a function of t is represented using the data from all tests. The dotted line presented, characterized by the linear function shown in the figure, was achieved by the overall tendency of the values. The slope of this function is close to 1, showing the close similarity between both members of Equation (7). Therefore, the relationship between the parameters "a" and k, derived from the relative mass α, and the accumulated total area of the firebrands A ac is well-established.
The flame height H f variation was also analyzed and compared with the other parameters described above. As can be seen in Figure 15, the flame height follows a trend comparable to the relative mass loss decayα.
The variations in the flame height and in the total crosssectional area of the firebrands released are presented in Figure 16. The total cross-sectional area A t represents the sum of the areas of all the firebrands captured at each instant. It is possible to see that both parameters follow the same tendency; however, the A t peaks show a delay of some seconds when compared with the H f peaks. The time difference between A t and H f peaks is hereby denoted by "delay time." In Table 7, the delay times determined for each test are presented. One of the factors affecting the delay time is the period required for the particles to take off from the tree by detachment as a consequence of fire; this time is hereby denoted by "reaction FIGURE 12 | Variation of the relative mass of the trees α as a function of the accumulated total cross-section area A ac of the released firebrands captured in the PIV images.
FIGURE 13 | Variation of the relative mass of the trees α as a function of the accumulated total cross-section area A ac of the released firebrands captured in the PIV images during the 10 s-peak period.
TABLE 6 | Linear coefficient m (mm −2 ) and exponential coefficient k (mm −2 ) and respective correlation coefficients achieved in the fitting of the mass of the trees by the estimated values of firebrands released ( Figure 11) (Table 4) and the vertical distance between the canopy and the PIV image center. The reaction time was determined by calculating the difference between the delay time and the uplift time ( Table 7). The reaction times found for the several species were very similar except for the eucalyptus trees, which require more time to release firebrands due to a change in the fire intensity. On the other hand, the reaction time of the oak trees and the pine trees was about four times lower than that of the eucalyptus trees, and about two times smaller than that of the cork trees (except for the PIN02). The authors believe that these differences are much related to the thickness and robustness of the petiole that links the leaves to the branches. However, specific tests are required to better explore these observations.

FINAL DISCUSSION OF RESULTS
The results of several tests aiming at the analysis of the firebrands released during the burning of different trees common in the Mediterranean forests were presented and preliminarily discussed. A compilation of the main results (Table 8) and an integrated discussion are carried out in this section. Contrary to what is usually assumed, the eucalyptus trees were found to be the species releasing fewer firebrands, while cork and oak trees produced more firebrands. In many situations of fuel management in the wildland-urban interface, the oak and cork trees are used to reduce the fire risk propensity, but the results of these tests and other observations show that their predisposition to produce firebrands able to cause spot fires at a shorter distance (tens of meters) cannot be neglected. These results are corroborated by the analysis that was carried out in order to understand the dependence of the mass loss on the release of firebrands. Once again, the loss of mass recorded during the burning of the cork trees was the one that proved to be most influenced by the release of particles. Eucalyptus was the species in which this dependence was the lowest.
On the other hand, the terminal velocity of the firebrands released by the cork trees was shown to be higher than that of the trees in the OAK tests, and so a larger spotting distance is expected for the firebrands produced by oak trees.
It is worth highlighting that in these experiments, with the exception of the cork trees, only young trees were tested and these results are not extended to mature trees. For example, the thicker barks of eucalyptus trees or pine cones were not available in these tests, and their relevant role in the development of spot fires in real fire events was not addressed. On the other hand, these results must be interpreted considering the limited fire intensity achieved in these experiments. In real fire events, not only would several episodes with higher intensity and thus more convective airflow velocity be expected but also the existence of meteorological wind, which would lead to the production of more and larger firebrands, would be anticipated. Still, this study allows a comparative analysis of the species with higher potential to produce firebrands that are able to originate spot fires at a distance up to some tens of meters.
It was observed that most of the firebrands released were the leaves of the trees. Since the eucalyptus trees used were the species with a larger mass percentage of leaves (Table 2) and a lower mass percentage of thick branches, a greater release of firebrands would be expected in these trees. However, it was also noted that the eucalyptus leaves require a high energy value to release firebrands, which is shown by the larger reaction time (Table 8) and by a higher value (52 and 56%, respectively) of final mass m f of fuel to burn after the tests ( Table 1).

CONCLUSIONS
In this work, the number and characteristics of the firebrands produced during the burning of four different trees were analyzed. The authors observed that a young tree belonging to Eucalyptus globulus had a lower potential to produce spot fires, followed by Pinus pinaster, Quercus robur, and Quercus suber. These trees released mostly leaves and twigs capable of producing spot fires at a short distance.
The results presented are relevant when used in fire behavior prediction models that consider the mechanism of spotting. However, the results obtained for cork, oak, and pine trees are related to young specimens and cannot be directly transposed to adult trees.
The methodology based on PIV images proved to be useful in this analysis. However, since the analysis was made on a twodimension basis, it is expected that the number and mass of firebrands determined are underestimated since firebrands can be hidden behind each other. Moreover, some firebrands may have had an early abandonment of the smoke plume and so might not have been captured by the PIV images and considered in the analysis. Thus, the previously mentioned sequence of the potential for those trees to cause spot fires is valid.
The authors of this study intend to perform additional tests with other species commonly found in the Mediterranean and non-Mediterranean forests of not only trees but also shrubs and herbaceous plants that have the potential to cause spot fires. Moreover, the authors aim to analyze the effect of the age of the tree and the effect of wind on the firebrands released during the burning of several forest species. The test results like the ones presented herein will be used as a reference to harmonize the results of tests performed in different conditions.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author/s.

AUTHOR CONTRIBUTIONS
MA: coordinated the study, designed, supervised and participated in laboratory tests, analyzed data, led the writing of the manuscript. LP: participated in the laboratory tests, performed PIV data processing, revised the manuscript. DV: supervised the study, discussed the results, revised and contributed to the writing of the manuscript. All authors contributed to the article and approved the submitted version.