Leaf-Level Field Spectroscopy to Discriminate Invasive Species ( Psidium guajava L. and Hovenia dulcis Thunb.) from Native Tree Species in the Southern Brazilian Atlantic Forest

: Invasive species are known to have potential advantages over the native community and can be expressed in their leaf functional traits. Thus, leaf-level traits with spectral reﬂectance can provide valuable insights for distinguishing invasive trees from native trees in complex forest environments. We conducted ﬁeld spectroscopy measurements in a subtropical area, where we also collected trait data for 12 functional traits of invasive ( Psidium guajava and Hovenia dulcis ), and native species ( Psidium cattleianum and Luehea divaricata ). We found that photosynthetic pigments were responsible for the greatest interspeciﬁc variability, especially in the green region of the spectrum at 550 nm, therefore contributing to detection of invasive species. In addition, according to LDA and stepwise procedures, the most informative reﬂectance spectra were concentrated in the visible range that is closely related to pigment absorption features. Furthermore, we aimed to understand the leaf optical properties of the target invasive species by using a combination of narrow bands and linear regression models. P. guajava showed high correlations with speciﬁc leaf area, Car/Chl and relative water content. H. dulcis had a strong correlation with water content, speciﬁc leaf area and Chla/Chlb. Overall, this methodology proved to be appropriate for discriminating invasive trees, although parameterization by species is necessary.


Introduction 1.Background
From the last century onwards, plant invasions have been increasing worldwide [1,2] and have been considered a major threat to global biodiversity [3].Invasive alien species (IAS) are introduced by humans outside of their natural range and cause negative ecological, economic, and social impacts on the novel environment [4][5][6].Invasive trees are usually adaptive and competitively distinct in relation to co-occurring native species [7].Therefore, they can alter the richness, composition and abundance of the community and also produce changes in ecosystem functioning [8].

IAS in Protected Areas
Invasive species have been reported causing risks in protected areas all over the world [9].In Brazil, the establishment of protected areas was a main strategy for biodiversity conservation [10].However, in addition to habitat fragmentation and the presence of disturbed environments [11], these areas are also currently being threatened by biological invasions [12].Outside the protected areas, land uses also need attention, as they mostly determine a strong propagule pressure and contribute to the success of plant invasions in these areas [13][14][15][16][17].
The occurrence of invasive species has increased in forest ecosystems in recent years, and negatively affects species composition, forest microclimate, and soil chemistry [8].Additive or synergistic effects of habitat disturbance and species invasions, such as competition for light and nutrients, result in today's landscapes being dominated by exotic species, which is a direct consequence of this competitive exclusion [18].Changes caused by invasive plants that disrupt ecosystem functioning and promote the invasion and establishment of new invasive plants create a feedback loop between habitat disturbance and species invasion, i.e., "invasive collapse" of disturbed ecosystems [19,20].
IAS represent an ongoing challenge for the management of protected areas associated with high economic costs [20][21][22].In this regard, remote sensing can potentially contribute as a methodological approach for the detection of invasive species and their monitoring over large and hard-to-reach areas, including forests.

Leaf Optical Properties
Remote sensing (RS) in the optical domain, which can be applied to a range of environmental studies, focuses on the spectral properties of the leaf in relation to its biochemical content (chlorophyll, water, dry matter) and its anatomical structure [23].The diffuse reflectance of a leaf, modified by its internal properties, may contain features that are useful for mapping the functional properties of leaves [24].In a healthy green leaf, the spectral region from 430 nm to 660 nm is characterized by a strong chlorophyll a absorption peak.In contrast, the absorption range of chlorophyll b is more intense (450 nm and 650 nm).The carotenoids have an absorption peak at 450 nm.In the near infrared (NIR) region, the reflectance increases dramatically (700 nm to 1200 nm), with a reflectance peak at 900 nm.Properties extracted from leaves and their reflectance spectra are input parameters and serve to feed canopy models.Combinations of narrow bands, such as the proposition of spectral indices, consist of a classical approach in remote sensing and are widely used for extracting information related to photosynthetic pigment content, water content, and dry mass, mainly due to their simplicity, and because they usually produce good results for individual datasets [23].

Remote Sensing for the Detection, Mapping, and Monitoring of IAS
Remote sensing has been used in recent years to detect invasive plants [25][26][27].IAS detection and mapping is possible due to the dominance of invasive species that can form homogeneous and large patches in the area, which contrasts the seasonal phenology of IAS and the native plants, and biochemical, physiological, or structural traits that can distinguish the IAS from the native community [28].However, the platform and sensor used for detection must be compatible with the target species and/or the environment characteristics.Multiscale approaches, such as integrating data from different platforms and sensors, are potential methodologies for the discrimination and mapping of IAS [29][30][31].
Phenology-based approaches are widely used, for example, to map IAS according to flowering and leaf senescence [32][33][34].On the other hand, structural traits tend to be spatially separable within a native ecosystem using processing techniques based on image segmentation.However, their application is indicated when the IAS present aggressive and dominant characteristics [35], or they form groups that are distinct from the native community [36,37].
Recent studies, such as those of Omer [38], Tesfamichael [39], and Iqbal [40], have used methods based on field spectroscopy to discriminate invasive species.Omer [38] used continuous wavelet analysis and feature reduction techniques to discriminate five invasive plant species.Tesfamichael [39] investigated the potential of remote sensing in identifying native, non-native, and hybrid Tamarix species in South Africa.Classification of species at the leaf and canopy level was performed using field spectroradiometer data.Iqbal [40] used a hyperspectral field sensor to discriminate invasive plant species from adjacent native species in two protected areas in Pakistan.Spectral separability was calculated using the Jeffries-Matusita distance index based on selected wavebands.
In tropical and subtropical forest ecosystems, mapping target IAS seems to be complex as they are spread out in the woods and most have morphological traits that are similar to the native species.Therefore, methods that not only use structural traits, but also spectral data should be used [41,42].Furthermore, to increase the potential for discrimination between native and invasive plants, physiological characteristics linked to plant functioning and growth must be included [32].Invasive tree species can be detected directly via the forest canopy [43] or indirectly by measuring the leaves and simulating the canopy [44,45].
Functional traits at leaf-level, such as photosynthetic pigments [46], nitrogen concentration [47] and moisture content [48], are directly related to physiological functions such as photosynthesis, respiration, and transpiration.These characteristics can be reflected in the optical properties of the leaves and, when related to the spectral signatures, can contribute to improving the discrimination and mapping of IAS [49,50].
Given the importance of the topic, the present study fills a gap related to remote detection of IAS in complex forest ecosystems.Moreover, it intends to assess invasive forest species that exhibit sparse invasion characteristics and sometimes form clusters; however, these species are not characterized as monodominant.We pursued the following objectives: 1) determine leaf functional traits and spectral characterization of native and invasive trees; 2) discriminate invasive species from native species based on their reflectance values via field spectroscopy; and 3) understand the relationship between functional characteristics of invasive plants and their spectral patterns (optical types).

Target Species
We selected two invasive species, namely Hovenia dulcis Thunb.(Japanese raisin tree) and Psidium guajava L. (Guava), and two co-occurring native species, namely Luehea divaricata Mart.& Zucc.(Whips-horse) and Psidium cattleianum Sabine (Strawberry guava).Native-invasive species pairs were chosen according to morphological and taxonomic parameters (most similar to plant comparisons), and occurrence distribution (frequencies and population densities).(See Table S1 and Figure S1 in Supplementary Materials).

Target Species
We selected two invasive species, namely Hovenia dulcis Thunb.(Japanese raisin tree) and Psidium guajava L. (Guava), and two co-occurring native species, namely Luehea divaricata Mart.& Zucc.(Whips-horse) and Psidium cattleianum Sabine (Strawberry guava).Native-invasive species pairs were chosen according to morphological and taxonomic parameters (most similar to plant comparisons), and occurrence distribution (frequencies and population densities).(See Table S1 and Figure S1 in Supplementary Materials).

Hovenia dulcis Thunberg
Hovenia dulcis, commonly known as Japanese raisin tree, is native to East Asia and grows in China, Japan, North Korea, South Korea, Thailand, and Vietnam [58].It is described as a deciduous fast-growing tree that can reach 25 m in height, is tolerant of shade though prefers habitats with high incidences of light [59,60].Leaves are alternate, simple, 10-15 cm long and 7-12 cm in width.It reproduces sexually by seed, and fruits are The Target IAS of the Study Hovenia dulcis Thunberg Hovenia dulcis, commonly known as Japanese raisin tree, is native to East Asia and grows in China, Japan, North Korea, South Korea, Thailand, and Vietnam [58].It is described as a deciduous fast-growing tree that can reach 25 m in height, is tolerant of shade though prefers habitats with high incidences of light [59,60].Leaves are alternate, simple, 10-15 cm long and 7-12 cm in width.It reproduces sexually by seed, and fruits are produced in large quantities [61].It is considered to be an invasive species in forest ecosystems in South America, especially in the Atlantic Forest domain [62][63][64].The species was brought to southern Brazil and introduced to rural properties motivated by economic reasons, shade and wood uses [61].As a consequence, it has spread beyond cultivation areas and has become a growing problem in subtropical deciduous seasonal forests, often being found in the canopies of secondary forest fragments [64][65][66].H. dulcis can outcompete other species for light and nutrients [67] and promotes changes in structure and in species composition in the plant community [68].

Psidium guajava L.
Psidium guajava (commonly known as guava) is a tropical tree from southern Mexico and northern South America that is under cultivation around the world and has become invasive in southern Brazil [69].Adult guava trees grow to 3-8 m in height.Leaves are opposite, simple, 5-15 cm long and 3-7 cm wide [70].It reproduces sexually and fruits are fleshy, edible, and produce large quantities of seeds [71][72][73].The species is well adapted to a wide range of environmental conditions, and it can form dense monocultures and produce allelopathic effects that displace native plants [74].

Spectral Measurements in the Field
Field sampling was conducted in QCSP in the spring of 2020 when plants were actively growing and reaching greater photosynthetic capacity.The average field temperature and relative humidity were 32 • C and 49%, respectively.We collected data from six mature trees of each species.All individuals were sampled under similar solar radiation/light conditions, i.e., measurements occurred at the same time (they started at 10:00 am and ended at 1:30 pm) on cloudless and windless days in a homogeneous area with 20 years of secondary forest succession (Figure 1b,d).
Spectral data collection was conducted in three campaigns in the spring, resulting in the collection of 24 subsamples.Care was taken in the design of the sampling to represent the different light incidence directions when collecting the samples.For this purpose, six individuals (samples) per species were selected and only leaves from the upper third of the canopy were collected.The subsample was distributed among the quadrants (north, south, east, and west), with a replicate of four leaves for all variables analyzed, thus making a total of 384 measured leaves.Sampling for leaf functional traits was conducted simultaneously in the same field.

Leaf Sampling and Laboratory Processing
We analyzed 12 leaf functional traits for both native and invasive species, which were related to photosynthetic pigments, water content and vegetation structure (Table 1).These traits at leaf-level are relevant as they represent ecological strategies that can be correlated with measurable leaf spectral properties [50].For the study, the same six individuals of each tree species previously selected were used for the measurements.In all cases, leaf sampling was conducted on the middle third of the canopy, and thirty-two healthy leaves were collected per individual.Four leaves were used in each analysis.To avoid the degradation of photosynthetic pigments, leaf samples were immediately wrapped in aluminum foil, frozen in liquid nitrogen, and stored on thermal box in the field until they could be transferred to a −80 • C freezer in the Plant Biotechnology Lab at Santa Maria Federal University.Fresh leaf tissues (50 mg) were then homogenized in liquid nitrogen, incubated at 65 • C with dimethyl sulfoxide (DMSO) until the pigments were completely extracted, as per Hiscox and Israelstam [75], and estimated with Lichtenthaler's formula [76].Concentrations of chlorophyll a, chlorophyll b and carotenoids were quantified using a spectrophotometer (Spectrophotometer VM5, Celm E-205D, Bel Engineering, Monza, Italy) based on their absorbance at 663, 645 and 470 nm, respectively.
Leaf samples for water content and for structural measurements were stored in humidified bags in the field and kept in black plastic bags on thermal box to prevent wilting during transport.Leaf samples were first measured with a scanner/portable leaf area meter (Portable Leaf Area Meter -AM300, ADC BioScientific Ltd, Hoddesdon, UK) to determine the leaf area.Leaves were fresh weighed and then oven-dried at ±65 • C until constant weight before recording dry mass.We calculated leaf mass per area (LMA; g cm −2 ) as leaf dry mass divided by leaf area.Specific leaf area (SLA; cm 2 g −1 ) was calculated as area per unit mass.

Collection of Leaf-Level Spectral Data
For leaf spectral characterization of the target species, field spectroscopy was performed using a handheld portable spectroradiometer (ASD FieldSpec ® , Malvern Panalytical, Malvern, UK) within the 300-1200 nm range, with spectral resolution of 3 nm.However, the spectral range of 400-900 nm was delimited for this work due to noise observed in the equipment.One exposed branch from the higher third of the canopy for each tree was removed and four leaves were then collected.Immediately (less than 10 min), leaf samples were overlapped and placed on a black background and leaf-level spectral reflectance signatures were recorded for each sample.
The spectroradiometer was placed vertically over the target (leaf) at a distance of 10 cm, and data was collected by 25 degrees of field of view (FOV for field measurement).The positioning of the spectroradiometer was performed by moving in constant azimuthal movement (90 • ).The distance between the target and the object was kept constant to ensure that no reflectance came from the surroundings of the leaf.The spectroradiometer was calibrated with a Lambertian reference plate before each sample measurement.For each target, 10 readings were collected from the reference plate and leaf sample to obtain the average spectrum.All measurements were collected with RS3 software (ASD), and ViewSpecPro software (ASD) was used for post processing of spectra files.During the field work, meteorological information, such as temperature, humidity, wind speed and solar radiation inclination associated with each individual, was also collected.A thermo hygro-anemometer and a clinometer were used.

Statistical Analyses
Our study was analyzed via the following three steps: (1) analysis of the functional traits of the native and invasive trees using leaf traits measured in the field; (2) analysis of spectral reflectance; and (3) integrated analysis of leaf functional traits and spectral reflectance (Figure 2).All the analyses were performed using Microsoft Excel (Microsoft, Redemond, Whashington, USA) and R (R Core Team 2020).A descriptive statistical analysis of all 12 traits of the target species with a 95% confidence interval (CI) is also provided.Additionally, the Kruskal-Wallis test was performed to compare the species, followed by Dunn's post hoc test.
spectral band ratios were calculated for each variable.This procedure was performed for each of the invasive species and permitted the identification of the wavelength with the best correlation for each trait.Finally, simple linear regression (SLR) was used to determine the significant functional optical properties based on the most significant functional characteristics of the leaves according to the best simple ratio correlation (r).The performance of the models was validated using the coefficient of determination (r²), root mean square error (RMSE), and residual sum of squares (RSS).To evaluate the spectral signatures of the species, the analytical technique was used.Linear discriminant analysis (LDA) was performed to distinguish the spectra with two subsets; one containing 70% of the spectral data to build the model and one containing the remaining 30% of the spectral data used for validation.Then, stepwise procedures were used to select important wavelengths that best explain the differences among species.This process selects the variables one by one in accordance with the p-value < 0.05 until no variable can be entered.The leave-one-out cross-validation (LOOCV) was used to estimate the model generalization, which is a disposition of k-fold cross-validation, where k is the number of examples in the dataset (n).
To test the correlation between leaf characteristics and their respective spectral reflectance for invasive species, a simple relationship was investigated to estimate the biophysical and biochemical parameters using an interactive correlation environment (ICE) according to Ogashawara [77].The method was adapted for Microsoft Excel and consists of a correlation matrix to select the single band ratio.First, a correlation matrix was created for each independent variable in terms of reflectance values with a resolution of 1 nm.All spectral band ratios were calculated for each variable.This procedure was performed for each of the invasive species and permitted the identification of the wavelength with the best correlation for each trait.Finally, simple linear regression (SLR) was used to determine the significant functional optical properties based on the most significant functional characteristics of the leaves according to the best simple ratio correlation (r).The performance of the models was validated using the coefficient of determination (r 2 ), root mean square error (RMSE), and residual sum of squares (RSS).

Leaf Functional Traits for Native and Invasive Species
We measured twelve leaf traits that were closely related to plant photosynthesis in terms of water content, biochemical content, and leaf structure and physiology based on dry mass.The statistical summary of the analyses can be found in the Supplementary Materials (see Table S2 in Supplementary Materials).

Water Content
In general, the water content differed between species according to the Kruskal-Wallis test (Table 2).However, after post hoc testing, the significant differences were found for the EWT trait (Table 3).The native species L. divaricata had the highest LWC for the total dataset and the highest intraspecific variability.Nonetheless, when comparing the means, the invasive H. dulcis presented the highest mean value among the species (Figure 3a).The native species showed greater FMC as well as greater intraspecific variability.In contrast, invasive species had lower values for data amplitude and lower intraspecific variability (Figure 3b).The results also showed that, in general, both invasive species presented lower EWT values when compared to the corresponding native species, and less intraspecific variability (Figure 3c).The traits related to vegetation structure showed significant difference between species (Table 4).After the post hoc test (Table 5), the greatest differences were observed among the groups: invasive species H. dulcis and P. guajava; invasive species H. dulcis and native P. cattleianum; followed by the native species L. divaricata and P. cattleinum.

chi-
For leaf chlorophyll a content (23.75 ± 8.98 µg cm −2 ), high values were indicated for L. divaricata (Chla = 30.38µg cm −2 ) and P. guajava (Chla = 25.96µg cm −2 ), then for H. dulcis (21.45 µg cm −2 ) and P. cattleianum (15.55 µg cm −2 ; Figure 5a).For chlo-rophyll b, high contents were also observed for L. divaricata (Chlb = 9.46 µg cm −2 ) and P. guajava (Chlb = 8.88 µg cm −2 ), followed by P. cattleianum (Chlb = 6.42 µg cm −2 ) and H. dulcis (Chlb = 5.9 µg cm −2 ; Figure 5b).Again, L. divaricata (Chl = 115.3µg cm −2 ; Chl a + b = 63.8 µg cm −2 ) and P. guajava (Chl = 103.8µg cm −2 ; Chl a + b = 57.4µg cm −2 ) showed greater Chl and Chl a + b contents (Figures 5c and 6c).In general, IAS presented less intraspecific variability, on average, than native species.The lowest Chla/Chlb rates were associated with the Psidium congeneric species, such as P. cattleianum (Chla/Chlb = 2.58) and P. guajava (Chla/Chlb = 2.85; Figure 6a).Our results showed that although these species are described as initial secondary, in general, they could present shade-tolerant behavior.This suggests that P. guajava is also capable of extending the range of invasion from degraded areas and forest edges, where it usually occurs, to the understory of subtropical forests and, consequently, this should alert land managers to the risk.In contrast, H. dulcis showed the highest Chla/Chlb rates (Chla/Chlb = 2.58).This was expected as this IAS is associated with light demand characteristics, and, therefore, its invasion process is favored by the degradation and disturbance of forest ecosystems The lowest Chla/Chlb rates were associated with the Psidium congeneric species, such as P. cattleianum (Chla/Chlb = 2.58) and P. guajava (Chla/Chlb = 2.85; Figure 6a).Our results showed that although these species are described as initial secondary, in general, they could present shade-tolerant behavior.This suggests that P. guajava is also capable of extending the range of invasion from degraded areas and forest edges, where it usually occurs, to the understory of subtropical forests and, consequently, this should alert land managers to the risk.In contrast, H. dulcis showed the highest Chla/Chlb rates (Chla/Chlb = 2.58).This was expected as this IAS is associated with light demand characteristics, and, therefore, its invasion process is favored by the degradation and disturbance of forest ecosystems Generally, leaf carotenoid content is found in the same proportion as chlorophyll a in the PSI and PSII photosystems.However, we found carotenoids to be around 2.5 times more abundant than Chl a.The highest contents were described for L. divaricata (Car = 11.49µg cm −2 ) and P. guajava, which also had the highest Car/Chl ratio (Car = 11.59 µg cm −2 ; Car/Chl = 0.18; Figures 5d and 6b).

Leaf Spectral Reflectance
Differences in the spectral feature of species occur mainly in magnitude, as observed for the species shown here in Figure 7a.Two spectral regions showed a greater difference in features between species, the region of the spectral domain of green and NIR.In the green range, a gradual increase was observed for all species (Figure 7b), which was even higher for the invasive species.Furthermore, there was an approximation between the native species (P.cattleianum and L. divaricata) and, consequently, a greater distance in terms of magnitude, compared to the invasive species (P.guajava and H. dulcis).Results also showed that, in general, the NIR region of the spectrum was able to distinguish the invasive species with deciduous characteristics.The highest reflectance values were for the invasive H. dulcis (55%) when compared to the native L. divaricata (45%; Figure 7c).

Discrimination of Reflectance Spectra
For the discrimination between native and invasive species, the model created using LDA achieved an overall accuracy of 97% via the 70/30 method, which indicates accurate discrimination.While LD1 clearly separated the native L. divaricata, LD2 was effective for the invasive H. dulcis.For the discrimination of the invasive P. guajava and the native P. cattleianum, the combination of the first and second discriminants was necessary (Figure 8).

Discrimination of Reflectance Spectra
For the discrimination between native and invasive species, the model created using LDA achieved an overall accuracy of 97% via the 70/30 method, which indicates accurate discrimination.While LD1 clearly separated the native L. divaricata, LD2 was effective for the invasive H. dulcis.For the discrimination of the invasive P. guajava and the native P. cattleianum, the combination of the first and second discriminants was necessary (Figure 8).
For the discrimination between native and invasive species, the model created using LDA achieved an overall accuracy of 97% via the 70/30 method, which indicates accurate discrimination.While LD1 clearly separated the native L. divaricata, LD2 was effective for the invasive H. dulcis.For the discrimination of the invasive P. guajava and the native P. cattleianum, the combination of the first and second discriminants was necessary (Figure 8).The confusion matrix of the best-performing LDA model (Table 8) demonstrates the ability of generalization for the total set of samples, which achieved 95% classification accuracy.The stepwise method was applied to identify the most informative variables that contributed to distinguishing IAS from native species.Thus, 32 spectral bands from the visible range were selected as the most significant for discriminating the species in this study (Table 9).Moreover, these bands are where pigment absorption features can be detected.

Leaf Optical Properties of the Invasive Target Species
The most significant simple ratio (SR) was individually selected for each leaf trait of each IAS.The most-significant optical properties for each group of leaf traits (water content, vegetation structure and photosynthetic pigments) are shown in in Table 10 and Figure 9.The total set of leaf traits analyzed, and their respective SR can be found in the Supplementary Materials (see Tables S2 and S3, and Figures S2 and S3).9), P. guajava expressed the strongest correlations in NIR.The NIR reflectance ratio between 817 and 802 nm was correlated with LWC (r = 0.91), thus it was one of the best bands to predict the water content in P. guajava leaves (Figure 9a).However, we know that this result may not be exclusively related to water content.For SLA, a narrow waveband of NIR (R822/R801) showed good accuracy (r² = 0.95) (Figure 9c).For Car/Chl, the spectral behavior of the species suggests an increase in reflectance at 839 nm and 785 nm) (r = 0.75), and an increase in carotenoids in leaves at the same time (Figure 9e).Table 10.Summary of the most-significant leaf optical properties and leaf characteristics of two invasive tree species in QCSP.LWC (%) = liquid water content; SLA (cm 2 .g−1 ) = specific leaf area; Chla/Chlb = chlorophyll a:b ratio; Car/Chl = carotenoid:total chlorophyll ratio; SR= simple ratio, n = number of samples, r = simple ratio correlation; r 2 = coefficient of determination; RMSE = root mean square error; CV = coefficient of variation; RSS = residual sum of squares.Significant results at α = 0.01, and at a 99% confidence level (p-value < 0.01).According to the simple linear regression (Figure 9), P. guajava expressed the strongest correlations in NIR.The NIR reflectance ratio between 817 and 802 nm was correlated with LWC (r = 0.91), thus it was one of the best bands to predict the water content in P. guajava leaves (Figure 9a).However, we know that this result may not be exclusively related to water content.For SLA, a narrow waveband of NIR (R822/R801) showed good accuracy (r 2 = 0.95) (Figure 9c).For Car/Chl, the spectral behavior of the species suggests an increase in reflectance at 839 nm and 785 nm) (r = 0.75), and an increase in carotenoids in leaves at the same time (Figure 9e).

Traits
Additionally, the other invasive species, H. dulcis, showed a strong correlation with LWC (r = 0.92) and SLA (r = 0.88) at the same reflectance ratio (R706/R532) in the visible region of the spectrum (Figure 9b,d).Lastly, Chla/Chlb was best correlated with R818 (NIR) and R769 (red) (Figure 9f).Analyzing the increase in reflectance in NIR, we can relate it to a decrease in chlorophyll a content and, consequently, an increase in chlorophyll b.

Discussion
The detection and mapping of IAS still represents a major challenge in remote sensing, especially in complex forest ecosystems [42,78].Tropical and subtropical forests, mainly in the Atlantic Forest domain, still have not been studied in depth [79].In this paper, the use of field spectrometry for native and invasive tree species discrimination was developed based on leaf functional traits and proved to be appropriate.

Characterization of Leaf Functional Traits and Spectral Behavior
The optical properties of leaves are related to their biochemical composition and structure depending on the species and the phenological age of the leaves [28].Considering that the spectral curves of green leaves are generally similar in shape, the difference is mainly in magnitude, as was observed for the species in our study.The optical properties of water are well known [80], and EWT and FMC are two different ways to define leaf water content.As described by [81], these two attributes are perfectly correlated when LMA is constant.However, this was not observed since LMA was not constant for our species.In general, IAS presented lower FWC, FMC and EWT when compared to the native species.Mainly, low FMC values may suggest that invasive species are more susceptible to fire at leaf and canopy levels, thus increasing the risk of wildfires in invaded areas [80][81][82][83].
Leaf reflectance proved successful for generating LMA and SLA estimates [84].However, there appears to be little agreement between physical and empirical bases of the methods, beyond which spectral wavelengths fit best for estimation [85].Quantitative information on LMA provides a better understanding of the taxonomy of functional groups, regulation of physiological mechanisms and ecosystem functioning [86][87][88].On the other hand, SLA is considered one of the main functional traits that drive plant differentiation since it is directly related to the efficiency of water use and, therefore, consists of a variable potential for the spectral characterization of vegetation [89][90][91][92][93]. Overall, deciduous species showed higher SLA than evergreen species; however, when the species were classified based on their origins, the IAS had higher SLA values, thus summarizing ecological strategies such as those related to the acquisitive end of the leaf economic spectrum and faster growing species [94].The values obtained for LMA in this study confirm what was described by [28], i.e., that deciduous species are associated with lower LMA (H.dulcis and L. divaricata) when compared to evergreen species (P.guajava and P. cattleianum).
Photosynthetic pigments indicate that there is potential to differentiate species [80].A marked difference was noticed in the green region at 550 nm, as mentioned by [95].Additionally, the near-infrared reflectance plateau at 850 nm indicates the wavelength region with the greatest reflectance as well as the greatest differences between species [28].The red edge, around 750 nm, demonstrates a distinct reflectance peak between the co-competing invasive and native species (H.dulcis × L. divaricata and P. guajava × P. cattleianum).

Spectral Discrimination of IAS from Native Species
Our results show that the VIS wavelengths are an important spectral region for discriminating the species, especially for pigment absorption features [28,96,97].Moreover, differences in pigment content in the visible range are detected when comparing invasive and native species at the same site.[40] discriminated more plant species pairs (invasive species from the native species) within VIS regions (84%) than in NIR regions (60%).Our study also confirms the limitations of NIR as mentioned by [98] for 26 tree species in a tropical dry forest in Costa Rica.The reflectance in NIR has shown to be more promising for differentiating tree functional types based on leaf phenology (evergreen and deciduous species) than for differentiating IAS from the native community.

Leaf Optical Properties of IAS
Understanding the relationship between leaf functional traits and optical types of plants is still a knowledge gap [99].Thus, we sought to fill this gap mainly on IAS in this study.Psidium guajava expressed the strongest correlations for LWC and SLA in NIR.Previous studies also proposed models based on wavelengths located in the NIR region for SLA [87,100] and for LWC [101].Carotenoid and Car/Chl ratio contents represent indicators of photosynthetic activity and photoprotective mechanisms in plants [97].Recently, [102] proposed two models exclusively in the VIS region (500, 660 and 700 nm) and accurately described Car/Chl changes in the range from 0.15 (r 2 = 0.87) to 0.6 (r 2 = 0.82), requiring no species-specific parameterization.Although the response obtained in our study may suggest the activation of photoprotection mechanisms by the species [97].Hovenia dulcis, however, was optically significant for LWC and SLA in a narrow waveband of the visible region (R706/R530).Additionally, we found that the species seems to be water-efficient, even though water absorption is low in this region [103], and high correlations may be associated with leaf chlorophyll content [104] and structural aspects [28].Even so, we must consider that the species did not have excessive losses in terms of water content, and the increase in visible reflectance reflects the increase in liquid water content.

Conclusions
It was shown that field spectroscopy is a potential method for correct discrimination of IAS from their co-occurring native species in the subtropical Brazilian Atlantic Forest.Although this method requires target species parameterization, these results are consistent, and the technique may be extended to other IAS and sites within the same subtropical forest ecosystem.This study is also important for the development of Dossel models with the goal of increasing the space-time resolution, and can be used with airborne or space-based remote sensing.We emphasize the importance of adding species leaf traits to the spectral database, as this may advance the expansion of knowledge in the field.Finally, the outcomes observed in this study may contribute to monitoring and improving IAS control practices in protected areas.
From this study we can conclude that: 1. Linear discriminant analysis (LDA) can be used to accurately discriminate Psidium guajava and Hovenia dulcis from Psidium cattleianum and Luehea divaricata; 2.
The greatest discrimination for IAS is located in the VIS region, specifically in the red (705 nm) and green regions (553 nm), which are known for being highly sensitive to pigment content variation.This suggests improved separability in chlorophyll absorption pits, as it also may suggest that the difference in anthocyanin content may enhance discrimination between species; 3.
For both IAS, the LWC and SLA showed similar behavior; 4.
P. guajava correlated to Car/Chl (R = 0.75) in NIR (R839/R785), which may suggest photoprotection activation.This provides an efficient process of acclimatization to environments normally unfavorable for other species; 5.
H. dulcis was best correlated to Chla/Chlb at R818/R769 (R = 0.7), which may suggest that photosynthetic activity is maintained even under conditions of high luminosity and temperature; 6.
IAS showed good correlation with Car/Clh and Chla/Chlb.Therefore, we may extend the space-time resolution based on data from orbital and suborbital platforms.

Figure 1 .
Figure 1.Images show: (a) Quarta Colonia State Park (QCSP); (b) closer view of the study site; (c) aerial view image of woody cover and QCSP head office; and (d) secondary forest canopy.

Figure 1 .
Figure 1.Images show: (a) Quarta Colonia State Park (QCSP); (b) closer view of the study site; (c) aerial view image of woody cover and QCSP head office; and (d) secondary forest canopy.

Figure 2 .
Figure 2. Schematic diagram illustrating flow of methodology.Detailed data processing and analyses are described in Section 2.6.Overall, steps represent data acquisition, data organization, processing and analysis, and results.

Figure 7 .
Figure 7. Leaf-level spectral reflectance curves of the four target species: Psidium cattleianum (blue), Psidium guajava (green), Luehea divaricata (red) and Hovenia dulcis (yellow) with spectral ranges at 400-900 nm (a); zoom in green spectral range (b); and zoom in NIR spectral range (c).Native species are represented by the dashed lines, and invasive species by the solid lines.

Figure 7 .
Figure 7. Leaf-level spectral reflectance curves of the four target species: Psidium cattleianum (blue), Psidium guajava (green), Luehea divaricata (red) and Hovenia dulcis (yellow) with spectral ranges at 400-900 nm (a); zoom in green spectral range (b); and zoom in NIR spectral range (c).Native species are represented by the dashed lines, and invasive species by the solid lines.

Table 1 .
Synthesis of traits and functions recorded for plant species at Quarta Colonia State Park during the 2020 field season.
The evergreen species (P.cattleianum and P. guajava) in this study showed greater leaf mass per area (LMA; Figure

Table 8 .
Confusion matrix obtained from the LDA-LOOCV species classification model showing invasive and native species.Numbers in the diagonal of the matrix are the number of samples predicted correctly.Were: P. cattleianum-Psidium cattleianum; P.guajava-Psidium guajava; H. dulcis-Hovenia dulcis; L. divaricata-Luehea divaricata.

Table 9 .
Most-significant spectral bands selected by stepwise procedures.Blue, green and red are the wavelength regions.
According to the simple linear regression (Figure