Cutting the Greenness Index into 12 Monthly Slices: How Intra-Annual NDVI Dynamics Help Decipher Drought Responses in Mixed Forest Tree Species

: We studied the correspondence between historical series of tree-ring width (TRW) and the normalized difference vegetation index (NDVI, i.e., greenness index) values acquired monthly over an entire year by unmanned aerial vehicles. Dendrochronological techniques revealed differentiated responses between species and seasonality. Pinus engelmannii Carri è re and Juniperus deppeana Steudel were affected by warm temperatures (TMAX) during the winter prior to growth and benefited from precipitation (PP) during the seasons prior to the spring period. The standardized precipitation–evapotranspiration index (SPEI) confirmed the high sensitivity of P. engelmannii to drought (r = 0.7 SPEI). Quercus grisea Liebm. presented a positive association with PP at the beginning and end of its growth season. Monthly NDVI data at the individual tree level in the three species (NDVI ~0.37–0.48) statistically confirmed the temporal differences. Q. grisea showed a drastic decrease during the dry season (NDVI = 0.1) that had no impact on drought sensitivity in the same period, according to the climate-TRW relationship. We conclude that a relationship is plausible between the crown greenness index and radial growth, although more extended temporal windows of the NDVI should be explored. Differences in susceptibility to drought found among the species would presumably have implications for the composition of these forests under drought scenarios.


Introduction
Climate change has caused changes in the biogeography of ecosystems [1], and rates of modification in species composition and structure are expected to increase in intensity and speed [2], with a consequent conversion of vegetation [3].With their heterogeneity of species and ages, mixed forests make an ideal laboratory for understanding these processes of change in growth rates, phenological transitions, and relationships with the environment.
In particular, regions under drought stress, in which changes in seasonal precipitation [4,5] are associated with increasingly frequent heat waves [6], provide an excellent opportunity to improve our understanding of the susceptibility of tree species to these phenomena.These climatic adversities produce alterations in plant community composition, with implications for reduced resilience and recovery rates and thus for ecosystem functioning.Drought is recognized as an environmental disaster (see [7]) comprising different hydrometeorological variables and environmental factors, not to mention the stochastic nature of water demand [8].
However, these complex ecological mechanisms represent a research challenge given the multifactorial nature of the drivers underlying the forest response to the interconnectedness of biotic and abiotic factors [9].For example, in response to drought stress, the tree limits its nutrient supply, with a reduction in radial growth [10], defoliation, and Remote Sens. 2024, 16, 389 2 of 20 crown decline [11].In other words, although the phenomenon is multifactorial, climatic adversity can be manifested in symptoms of reduced radial growth and modification of photosynthetic activity with implications for the dynamics of carbon allocation [9].
An alternative recently pursued by the scientific community is the integration of dendroecological data with remote sensing [12][13][14], assuming an interconnection between tree structures explained by the rates of increase recorded in the growth rings and vegetation indices.However, testing these connections has been challenging due to several inconsistencies that require further research [15].For example, the limited temporal and spatial resolutions of satellite images act to increase their uncertainty in spectral resolution, especially in natural stands where there is a mix of species of uneven ages and variations in density [16].As a current response, the advent of unmanned aerial vehicle (UAV or drone) technology has come to resolve these inconsistencies through the capacity of these devices to monitor near-real-time variations in vegetation indices at the individual tree level with a high level of reliability and flexibility [17].
One of the most widely used indices is the NDVI as a proxy for photosynthetically active biomass; this greenness index is derived from measurements of the optical reflectance of solar light in the red and near-infrared wavelengths, for which reason it is very sensitive to the conditions of the ecosystem.However, certain inconsistencies or divergences, such as foliage production over growth ring increase, still require clarification [18,19].Thus, the combination of dendrochronological data and drones offers a promising alternative for refining our understanding of tree growth in the face of climatic variation and is now clearly positioned as a research strategy on the research agenda [9,20].
In practice, however, it is clear that scientific understanding of the conjunction of the TRW series with the NDVI poses three problems that merit further analysis to strengthen their dependence [21]: First, the spatial coverage of the tree-ring data must correspond spatially with the spectral data; second, the NDVI data are assumed to be equally correlated with TRW throughout the year, regardless of the time of sampling.Changes in phenology and combinations of species, as well as the discontinuous seasonality of satellite NDVI and the scarcity of continuous data at the tree level, produce a differential response in the seasonality of the NDVI-TRW relationship.Third, most of the studies that use the NDVI in forestry lack periodic replications that cover at least one phenological cycle in the same area, for which reason conducting continuous periodic flights from UAVs provides a comparative and novel advantage of seasonality that can contribute to scientific knowledge.
Consequently, a strategy is required to identify the contribution of the NDVI to seasonal radial wood growth, which has implications for the productivity of forest ecosystems.That is to say, the conjunction of continuous spectral data from UAVs with in situ field information on trunk growth as a product of photosynthesis is assumed to be an optimal tool to achieve a refined understanding of species growth in the face of hydroclimatic variations.In the future, this algorithm could be useful for adapting the inconsistencies found with satellite technologies [19,21].
Moreover, to our knowledge, there is no dendroecological research featuring multiple species converged in a single stand, despite the great diversity of these ecotones.Above all, these studies have excluded leafy species and tend to focus on species that are easily dated [22].Nevertheless, from an ecological perspective, hard-to-date species such as those of Quercus, Arbutus, and Juniperus have been reported as useful for fuel, shelter, and food for fauna, among other ecological connections [23].
In this study, we evaluated the coupling of NDVI-growth relationships by comparing multispecies tree-ring data with monthly NDVI series at the individual tree level.We hypothesized that the responses would be differential as a function of species and their seasonal responses to climate.
We posed the following specific research questions: (1) What is the sensitivity to climatic variables and drought of the three species in a mixed forest?(2) What is the variation in the NDVI over the year among the species?(3) Is there a plausible correspondence between the NDVI and tree-ring width?

Description of the Study Area
We selected a strategic area where species of three genera cohabit in one hectare located in the Sierra Tarahumara (Figure 1).This region is considered an important source of environmental goods and services for the communities in the area [24].The study site is characterized by young trees in a mixed forest.These forests are physiognomically dominated by the genera Pinus and Quercus, with varying proportions and species compositions, in association with species of the genera Arbutus and Juniperus [25].
We selected a strategic area where species of three genera cohabit in one hectare located in the Sierra Tarahumara (Figure 1).This region is considered an important source of environmental goods and services for the communities in the area [24].The study site is characterized by young trees in a mixed forest.These forests are physiognomically dominated by the genera Pinus and Quercus, with varying proportions and species compositions, in association with species of the genera Arbutus and Juniperus [25].
The studied species belong to the genera Pinus, Quercus, and Juniperus.Pinus engelmannii Carrière is a frost-resistant species found in a wide distribution range from 1200 to 3000 m asl in elevation.Its wood is used for commercial purposes [23].Quercus grisea Liebm. is distributed from the southern United States of America to central Mexico and can tolerate both cold damage and drought stress [26].Its wood can be used for producing cellulose.Juniperus deppeana Steudel is distributed from the southern United States of America to the states of Michoacán and Veracruz in Mexico, at elevations of 1200 to 2900 m asl.It is a drought-and fireresistant species that can tolerate compacted, nutrient-poor, and alkaline soils [23].It is not a commercial species, but it is used as fuel in mountainous areas and its wood can be utilized for the production of fences or small items of furniture of high durability.The studied species belong to the genera Pinus, Quercus, and Juniperus.Pinus engelmannii Carrière is a frost-resistant species found in a wide distribution range from 1200 to 3000 m asl in elevation.Its wood is used for commercial purposes [23].Quercus grisea Liebm. is distributed from the southern United States of America to central Mexico and can tolerate both cold damage and drought stress [26].Its wood can be used for producing cellulose.Juniperus deppeana Steudel is distributed from the southern United States of America to the states of Michoacán and Veracruz in Mexico, at elevations of 1200 to 2900 m asl.It is a drought-and fire-resistant species that can tolerate compacted, nutrient-poor, and alkaline soils [23].It is not a commercial species, but it is used as fuel in mountainous areas and its wood can be utilized for the production of fences or small items of furniture of high durability.

Dendrochronological Data and Processing
At least 14 dominant trees per species were randomly sampled for the dendrochronological analysis.From each tree, two increment cores were taken using a Pressler borer (Ø = 5.1 mm) as a strategy to reduce bias due to the eccentricity of the trunk pith.Diameter at breast height (DBH) and total height (TH) were also recorded (Table 1).The samples were mounted on wooden rails and dried at around 18 • C. The extracted cores were polished with different grades of sandpaper (from 80 to 1000 in grain size) to highlight their features.The growth rings were counted and dated using a standard dendrochronological technique [27], which consists of measuring the total ring width (TRW) to an accuracy of 0.01 mm using the Velmex measuring system (Velmex Inc., Bloomfield, NY, USA).Dating was checked using the COFECHA program, which compares the series of each tree with a master chronology [28].Standardization was then performed to eliminate biological and geometric growth trends not associated with climatic variables, using a negative exponential model to obtain the ring width indices (RWIs) for each species.Standardization and chronologies were developed using the dplR library [29][30][31].Finally, statistics were calculated, obtaining the mean (Mean), first-order autocorrelation (AC), mean correlation among trees (Rbar), and expressed population signal (EPS).The quality of the chronologies was evaluated through the EPS value, a statistic considered to reflect well-replicated chronologies (>0.85) [32].
Following the dendroecological convention, as part of the statistical analysis [27,28], Pearson correlation indices (r) were calculated to assess the influence of local climate on radial growth.For this, correlations were performed between the climatic variables of maximum temperature (TMAX), minimum temperature (TMIN), and precipitation (PP), and the residual annual growth series of each species.Monthly climatic data obtained from the KNMI-Climate Explorer website (http://climexp.knmi.nl/,accessed on 29 August 2023) of the CRU TS 4.07 base were used for the 0.5 • grid centered on the site at 27.0-27.5• N and 107.0-107.5 • W. The correlation analysis window was from August of the previous year to July of the growth year.
The SPEI [33] was used to determine the influence of drought severity and duration on radial growth.This is a robust strategy in dendroecology because it allows estimation of the influence of cumulative drought, given the ability of trees to respond differentially according to the seasonality of the water stress [34].The SPEI was obtained from the website http://sac.csic.es/spei/index.html(accessed on 30 August 2023) for the same grid as mentioned above.Values were obtained on time scales from 1 to 24 months.Finally, correlations were calculated between these values and the residual growth series for each species.

Remote Sensing Data and Procedures
To record the phenology of an entire year starting from the growth season, images were acquired monthly from April 2021 to April 2022, preferably on sunny days, around the 15th day of each month.One UAV flight per session was conducted in an east-to-west direction at around midday as a strategy to avoid changes attributable to the diurnal variation in reflectance since this time of day is reported to be ideal for recording spectral data [15].To eliminate radiometric interference that could reduce the potential interpretation of the images and distort the results of the remote sensing analysis, we decided to perform a periodic calibration of the spectral sensor, following the procedure referred to by Daniels et al. [35].This was fundamental to preserving the spectral quality of the images at different weather stations, and thus sustaining the capability for NDVI monitoring [36].
The equipment used was a DJI Phantom 4 multispectral quadcopter (Dajiang Innovation Technology Co., Shenzhen, China) equipped with an RGB digital camera and a set of five multispectral cameras, covering blue (450 nm ± 16 nm), green (560 nm ± 16 nm), red (650 nm ± 16 nm), red edge (730 nm ± 16 nm), and near-infrared (840 nm ± 26 nm) spectra, all with a global obturator of 2 MP [37].During each flight mission, at least 846 photographs were taken in JPG format at an altitude of 60 m, with the camera fixed at an angle of 90 • (e.g., the nadir angle) and overlaps of 80 and 70% in the horizontal and vertical orientations, respectively.Figure 2 illustrates the workflow followed to acquire and process the images for subsequent linkage with dendroecological data.The UAV employs an integrated georeferencing system during the flights and, therefore, we do not consider it necessary to use the real-time-kinematic (RTK) system, since this procedure allows vertical and horizontal accuracies of ±0.1 and ±0.3 m, respectively (see [37]).
interpretation of the images and distort the results of the remote sensing analysis, we decided to perform a periodic calibration of the spectral sensor, following the procedure referred to by Daniels et al. [35].This was fundamental to preserving the spectral quality of the images at different weather stations, and thus sustaining the capability for NDVI monitoring [36].
The equipment used was a DJI Phantom 4 multispectral quadcopter (Dajiang Innovation Technology Co., Shenzhen, China) equipped with an RGB digital camera and a set of five multispectral cameras, covering blue (450 nm ± 16 nm), green (560 nm ± 16 nm), red (650 nm ± 16 nm), red edge (730 nm ± 16 nm), and near-infrared (840 nm ± 26 nm) spectra, all with a global obturator of 2 MP [37].During each flight mission, at least 846 photographs were taken in JPG format at an altitude of 60 m, with the camera fixed at an angle of 90° (e.g., the nadir angle) and overlaps of 80 and 70% in the horizontal and vertical orientations, respectively.Figure 2 illustrates the workflow followed to acquire and process the images for subsequent linkage with dendroecological data.The UAV employs an integrated georeferencing system during the flights and, therefore, we do not consider it necessary to use the real-time-kinematic (RTK) system, since this procedure allows vertical and horizontal accuracies of ±0.1 and ±0.3 m, respectively (see [37]).As part of the processing, individual trees were manually geolocated on the orthomosaic, and delimitation of the crown area was complemented with the ForestTools algorithm of the statistical software R Studio, version 4.3.1 (The R Project for Statistical Computing, Vienna, Austria).This tool also helps to delimit the tree crowns through the variable window filter algorithm.Using the zone statistics tool of QGIS version 3.16.15,descriptive NDVI statistics for each tree of the genera P. engelmannii, Q. grisea, and J. deppeana were extracted from the center of each canopy, with the mean selected as the variable of interest.
The reflectance level was calculated according to the wavelength of each band through the multispectral orthophotos at the tree level.To calculate the multispectral indices, the semiautomatic classification add-on was run in the free software QGIS.This tool allows the calculation of reflectance in a given area according to the values of each band.The photosynthetic capacity of the vegetation was obtained through the NDVI.Although a variety of indices exist, this indicator stands out as a measure of chlorophyll abundance and absorption of photosynthetically active radiation and is thus an indicator of the potential biomass captured by the tree [38].The calculation was performed with the raster calculator of the QGIS software using the following equation: where NDVI = normalized difference vegetation index, NIR = near-infrared band (840 ± 26 nm), and RED = red band (=650 ± 16 nm).
To rule out autocorrelation between the monthly NDVI readings and the species studied, we subjected the data to an analysis of variance (ANOVA) statistical test with software R Studio, version 4.3.1 (R Core Team, Vienna, Austria.Available on: https://www.R-project.org/,accessed on 3 October 2022).Using the significant intra-annual and interspecies differences, it was possible to attribute variations in the NDVI among dates and species to changes inherent to the tree phenology as a biological explanation that strengthens the link between photosynthetic activity in the canopy and radial growth.Finally, to evaluate the differences between the different factors (species and dates), a Bonferroni post hoc test was performed as a strategy that allowed control of the simultaneous confidence level for the full set of each interval (see Table A1).

Results
The species present at the site had a well-replicated common period in their chronologies between 1975 and 2019.Pinus engelmannii was the longest-lived species, followed by J. deppeana and Q. grisea (Table 1).The mean total tree-ring width ranged from 1.20 mm for Q. grisea to 2.13 mm for P. engelmannii.Juniperus deppeana recorded a mean growth of 1.40 mm per year.Juniperus deppeana and P. engelmannii had higher proportions of earlywood and latewood (93% and 7%; 83% and 17%, respectively) compared with Q. grisea (16% and 84%).Figure A1 presents the interannual variability in growth, showing similar behaviors in some specific years.All three species presented a decrease in growth in the years 1996, 2002, and 2011, while an increase in species growth occurred in the years 1997 and 2010.These increases and decreases in growth were represented by narrow and wide growth rings, respectively.
Regarding the dendrochronological statistics, the highest values of the first-order autocorrelation (AC) were presented by J. deppeana (0.40).On the other hand, the highest mean inter-tree correlation (Rbar) and the best expressed population signal (EPS) were observed in P. engelmannii, at 0.66 and 0.99, respectively (Table 2).All species showed an EPS value suitable for dendroclimatic analysis (>0.85) [39].

Variable Climate-Growth Relationships
The behavior of the species in relation to the climatic variables is shown in Figure 3a.Pinus engelmannii presented the highest correlation with the climatic variables.This species was found to be susceptible to TMAX in the period from September of the previous year to April of the growth year, except for December and February, with correlation values of r = −0.22 to r = −0.39.This same species showed a positive correlation with TMIN in December of the previous year (r = 0.24).Finally, rainfall during the previous growth season (from August of the previous year to May of the current year) had a positive influence on the radial growth of P. engelmannii, with correlation values ranging from r = 0.24 to r = 0.45.−0.22 to r = −0.39.This same species showed a positive correlation with TMIN in December of the previous year (r = 0.24).Finally, rainfall during the previous growth season (from August of the previous year to May of the current year) had a positive influence on the radial growth of P. engelmannii, with correlation values ranging from r = 0.24 to r = 0.45.
Juniperus deppeana showed a negative correlation between radial growth and TMAX prior to the growth season (January r = −0.28 and February r = −0.29 of the current year), and a positive correlation with PP for some months of the winter period prior to the growth season (October of the previous year r = 0.26 and January of the current year r = 0.37).
Finally, Q. grisea showed a negative correlation between radial growth and TMAX in the initial months of growth (March and April of the current year, r = −0.22 to −0.25) and a positive correlation with PP in November of the previous year (r = 0.29), and at the beginning of the growth season in March, April, and May of the current year (r = 0.25, r = 0.34, and r = 0.22, respectively).The growth of the three species showed a differentiated response to the drought index among species (Figure A2).Pinus engelmannii was shown to be the species most sensitive to the drought index, with high correlation values of up to 0.7 for the period October Juniperus deppeana showed a negative correlation between radial growth and TMAX prior to the growth season (January r = −0.28 and February r = −0.29 of the current year), and a positive correlation with PP for some months of the winter period prior to the growth season (October of the previous year r = 0.26 and January of the current year r = 0.37).
Finally, Q. grisea showed a negative correlation between radial growth and TMAX in the initial months of growth (March and April of the current year, r = −0.22 to −0.25) and a positive correlation with PP in November of the previous year (r = 0.29), and at the beginning of the growth season in March, April, and May of the current year (r = 0.25, r = 0.34, and r = 0.22, respectively).
The growth of the three species showed a differentiated response to the drought index among species (Figure A2).Pinus engelmannii was shown to be the species most sensitive to the drought index, with high correlation values of up to 0.7 for the period October to December on a scale of 14 to 16 months.Juniperus deppeana and Q. grisea both showed a moderate response to drought.On the one hand, J. deppeana presented a higher association on the scale of 8 to 10 months duration for June and July, while moderate associations were identified in Q. grisea on the scale of 2 to 6 months for May and June.

NDVI and Statistical Analysis
A total of 13 monthly flights were conducted in the period April 2021 to April 2022 as a strategy to capture phenological variations over a year.Although trees tend to respond seasonally to drought, here, we refined the flight period to every thirty days, assuming the ability to record obvious variations in the NDVI during this period.To ensure consistent imagery with the UAV, monthly calibration of the sensors was conducted prior to drone flight.The product-specific application "DJI Assistant 2" was used in this process (https://www.dji.com/mx/downloads/softwares/assistant-dji-2,accessed on 7 April 2021).The NDVI was obtained at the individual tree level for 270 trees, with a higher proportion of P. engelmannii, J. deppeana, and Q. grisea found, as is typical in the forests of the Sierra Madre Occidental.The analysis was classified by genus and monthly trend plots were generated (Table 3 and Figure 3b).The highest values were obtained for P. engelmannii, with a mean of 0.48 and minimum and maximum values of 0.18 and 0.68, respectively.The mean NDVI values found for Q. grisea and J. deppeana were 0.45 and 0.37, with minimum values of 0.30 and 0.14 and maximum values of 0.58 and 0.61, respectively (Table 3).The average DBH of the 270 trees per species was greater in Q. grisea than in P. engelmannii and J. deppeana (15.86 ± 0.90 cm, 11.37 ± 0.13 cm, and 7.94 ± 0.19 cm, respectively) and the average total height was greater in Q. grisea than in P. engelmannii and J. deppeana (6.55 ± 0.16 m, 6.43 ± 0.06 m, and 3.63 ± 0.05 m, respectively; Table 4).For those NDVI values that may reach saturation and are not attributable to phenological stages, it is recommended to apply the solution algorithm proposed by Liu et al. [40] (but see [41]).On the other hand, the highest mean NDVI values were found for the species P. engelmannii, Q. grisea, and J. deppeana (0.48, 0.45 ± 0.02, and 0.37 ± 0.01, respectively; Table 4).Regarding the NDVI for P. engelmannii and J. deppeana, smoothed inter-monthly oscillations were observed, being of lesser magnitude for J. deppeana (Figure 3b).For Q. grisea, the NDVI decreased drastically from April to June, which can be explained by the gradual loss of leaves.
The analysis of variance revealed statistically significant differences among species, between sampling dates, and in the species × sampling date interaction, proving that the monthly NDVI variations differ significantly for the two sampling dates, and among the species and species dynamics throughout the year (Table 5, Figure A3).The NDVI statistical values differed among the three genera in May, June, July, August, October, and December 2021, as well as in March 2022 (Figure A4).These differences occurred among one of the three species for April, September, and November 2021, and in January, February, and April 2022 (Figure A4).In April 2021 and 2022, Q. grisea and J. deppeana showed similar behavior, while statistical differences were found for the species P. engelmannii (2021: p = 0.0002 and 0.0020, and 2022: p = 0.0360 and <0.0001, respectively).At the end of summer (September), and in the winter period (November 2021 and January and February 2022), the species Q. grisea and P. engelmannii behaved similarly and no significant differences were found between them, while J. deppeana behaved significantly differently from the other species (Figures A4 and A5).
Significant statistical differences were found in the monthly NDVI values per species throughout the year.Figure 3b shows the variations found for the species studied.For Q. grisea in the period May to July, there were very marked differences corresponding to the highest SPEI correlations occurring in the same period of limited soil moisture.

Discussion
Tree vigor is usually reported using qualitative visual criteria in the upper canopy (e.g., current percentage of crown transparency or defoliation; see [9]), conferring an inherent bias to the subjectivity of the assessors, as well as their partial perspective of the canopy.However, it is important to note that remote sensing through NDVI per se does not account for tree vigor as such but instead presents reflectance rates that should not be considered as true photosynthetic rates.
Our study makes an advantageous contribution with the numerical intra-annual assessment provided by the NDVI, continuously over the same tree and from an aerial perspective of the crown.This is most useful in mixed forests, where species coexistence provides multiple ecological connections that are still worthy of further study.For example, species composition is a useful management strategy that decision-makers should have in their portfolio of tools to address climate change [42].
The relationship between trunk growth and photosynthetic activity is the main mechanism that scientists seek to understand; that is, the interdependence of trunk TRW data with those of canopy NDVI helps to discern responses to predicted drought scenarios and their implications for vegetation productivity rates, as seen in [43].The results found here show that NDVI data taken with drones at the individual tree level are sensitive to the biomass growth recorded by the trunk, so the spatial and spectral resolution seems to contribute significantly to knowledge in this area.The NDVI has shown strong connections with TRW, so the literature gives it the ability to spatially predict the trend of the RWI (see [21]).
Furthermore, there is a clear seasonal association between TRW and the NDVI, attributed to phenological responses and interspecies drought stress.Thus, it is demonstrated that NDVI values obtained continuously with drones are suitable for monitoring vegetation dynamics in ecosystems such as the one studied here.Notably, it has been demonstrated that the NDVI has a seasonal relationship with tree radial growth; therefore, the optimal time windows for the use of monthly NDVI data in association with the TRW of the species studied have been identified.In this case, the radial growth trends in Figure A1 show evident episodes of recovery and synchronized disturbances.The studied area is droughtlimited, as shown by the correlations of TMAX and SPEI in Figures 3a and A2, and has been previously documented as such [44].Thus, the results of the NDVI dynamics in Figure 3b support that the drought factor differentially affects the species.
Juniperus deppeana showed the greatest variation in NDVI values throughout the year (Figure 3b); in contrast, P. engelmannii had the most consistent mean NDVI values, which is attributed to the fact that it retains its leaves throughout the year.Quercus grisea showed the greatest variation in NDVI, with the lowest values in the dry months (Figure 3b).
The drastic reduction in the NDVI in Q. grisea during the dry season is explained by the loss of foliage, and, although a reduction in radial growth during the dry season is plausible, the species did not show sensitivity to drought (Figure A2).One possible interpretation is that the oaks can redistribute their carbon stocks as a drought adaptation strategy [45].This is consistent with the fact that growth begins prior to the formation of new foliage [46].Furthermore, given the loss of foliage, evapotranspiration rates are reduced with a consequent and progressive saving of carbon stocks [47].However, this growth in earlywood formation could potentially lead to carbohydrate depletion, making this species susceptible to dieback in the event of severe drought.The minimum carbohydrate thresholds that species must meet to survive should therefore be the subject of future study [48].
Pinus engelmannii had the highest NDVI values, followed by J. deppeana and finally Q. grisea, evidencing differentiated rates in the canopy as a photosynthetic driver [49,50].In the cases of P. engelmannii and J. deppeana, high NDVI values (Figure 3b) in spring suggest correspondence with earlywood formation, which is consistent with Michelot et al. [45].With the available data, we therefore speculate that there is a direct transfer of newly synthesized assimilates for radial growth.Moreover, the photosynthetic capacity of these species during the pre-growth seasons seems to reinforce carbon reserves with which to form the new tree ring in spring [51].Likewise, the increase in the NDVI during autumn (Figure 3b) is a strategy that strengthens cell lignification before winter.However, these results are inconclusive, and other variables such as age, a factor that affects carbohydrate distribution and efficiency, still require consideration.
The distribution of carbohydrates in tree structures can make it difficult to match NDVI and TRW values [52], for which reason it is recommended to extend observations of phenology and xylogenesis, as well as to extend the spatial and temporal resolution of the remote sensor.In other words, the delay of carbon allocation in radial stands may take longer than the period of analysis utilized here, depending on the physiological conditions of the tree, including its distribution priorities throughout its structures.It has been documented, for example, that trees allocate reserves in times of stress [53].Thus, the UAV data were operationally valuable and allowed us to elucidate how reflectance and its association with TRW diverge among species (Figure 3b, Table 3).Our NDVI series should therefore be temporally robust, and it is recommended to progressively continue the flight missions beyond one year, given the occurrence of interannual ecophysiological processes, as well as to subdivide the growth rings sub-annually.
Phenology changes according to differences in carbon distribution, which is inherent to each particular species [45], especially in our case, where edaphic and climatic conditions did not vary.From the available data, P. engelmannii and J. deppeana showed a synchronized trend of NDVI values, with pronounced reductions in winter and summer.These findings could be associated with sensitivity to TMAX during the same periods (Figure 3a); i.e., the warm conditions during winter and summer trigger reactivation of the cambium [54], decreasing the availability of carbohydrates in the aerial part of the tree (Figure 3a).Under this premise, intra-annual growth does not have a time lag with the phenology of these species, although more temporal NDVI data are required.Quercus grisea, a contrasting species in terms of phenology and anatomy, seems to show a different carbon assimilation strategy that is reflected in seasonal delays in phenology and radial growth [55].
This approach remains to be replicated in other ecotones, since it has been shown in agricultural sciences [56,57] and viniculture [58,59] that further replication efforts are still required to optimize the potential uses of this type of sensor.For example, our methodology could be used in agricultural crops where drought is a limiting factor [60], albeit with the pertinent adjustments made to the flight configurations.
Variations in climatic conditions throughout the year can affect the quality of NDVI data acquired by drones.This is a limitation that must be taken into account.Maintaining the integrity and reliability of the data is critical, so the best conditions for image acquisition, including complementary robust sensors, should be sought.
Photosynthesis is a function of the ecophysiological processes occurring in the tree, which in turn correspond to its seasonality [61].In other words, carbohydrate translocation may be temporally out of phase between the TRW and photosynthesis, the understanding of which becomes more challenging during periods of drought stress.Therefore, groundtruth photosynthetic measurements should be estimated at the canopy, and eventually at the leaf level, and should also be complemented with subsequent investigations.For example, vessel lumen size is directly related to the capacity for hydraulic conductivity [62], highlighting the need for studies of wood anatomy and formation (xylogenesis) to refine metrics and associate these with changes in NDVI rates across seasons (Figure 3b).
The overlapping of canopies in very dense forests is another limitation of our methodology, which could be overcome with radar sensors, for example, although the cost remains to be assessed.

Conclusions
The results of this study helped to address the three research questions posed and were of ecological significance.Inclusion of the NDVI as an explanatory factor contributed to deciphering responses to the climate experienced by mixed forests.Of particular note is the novel inclusion of periodic flights on a monthly scale over a year covering the different phenological stages of the same trees, which refined certain ecological mechanisms related to drought that have not been addressed in previous studies.
Vulnerability to drought was differentiated into two divergent groups: on the one hand, evergreen conifers suggested the correspondence of crown phenology with radial growth, presumably with no pauses in carbohydrate allocation, while on the other, the oaks are less susceptible to drought and have no correspondence between the greenness index and treering growth, which is attributed to delays in C distribution; i.e., P. engelmannii and J. deppeana seem to depend on crown phenology, whereas the radial growth of Q. grisea is governed by its C reserves.We could speculate that radial growth was mainly photosynthesis-dependent, although with different strategies presented between coniferous and porous trees.We therefore recommend complementing these analyses with additional data pertaining to anatomy, xylogenesis, and age.
The NDVI at the individual tree level is an ally of dendroecology in terms of explaining drought.Thus, multiple proxies seem to amplify the outstanding responses of trees to adverse climatic conditions.According to these results, changes in species composition may occur in mixed forests if drought duration and severity are magnified.In addition, changes in carbon fluxes and productivity rates could be modified.However, it is advisable

Figure 1 .
Figure 1.Map of study site location, UAV flight paths, and sampled tree crowns per species.Figure 1. Map of study site location, UAV flight paths, and sampled tree crowns per species.

Figure 1 .
Figure 1.Map of study site location, UAV flight paths, and sampled tree crowns per species.Figure 1. Map of study site location, UAV flight paths, and sampled tree crowns per species.

Figure 2 .
Figure 2. Workflow of the procedure followed to generate and process the NDVI data at the individual tree level for twelve months and to link the data to the tree rings to decipher the relationship between crown activity and radial activity.The images were processed and analyzed using digital photogrammetric and computer vision techniques with the open-source software OpenDroneMap (ODM version: 2.8.4;Cleveland Metroparks, Cleveland, OH, USA, https://github.com/Open-DroneMap/ODM,accessed on 18 January 2022) version 2.8.4.LTS.This software is open access and its robust operation utilizes the structure from motion and multi-view stereo (SfM and MVS) algorithms, obtaining 3D point clouds of 1000-20,000 points m −2 , resulting in a multispectral orthomosaic of good accuracy.As part of the processing, individual trees were manually geolocated on the orthomosaic, and delimitation of the crown area was complemented with the ForestTools algorithm of the statistical software R Studio, version 4.3.1 (The R Project for Statistical

Figure 2 .
Figure 2. Workflow of the procedure followed to generate and process the NDVI data at the individual tree level for twelve months and to link the data to the tree rings to decipher the relationship between crown activity and radial activity.The images were processed and analyzed using digital photogrammetric and computer vision techniques with the open-source software OpenDroneMap (ODM version: 2.8.4;Cleveland Metroparks, Cleveland, OH, USA, https://github.com/OpenDroneMap/ODM,accessed on 18 January 2022) version 2.8.4.LTS.This software is open access and its robust operation utilizes the structure from motion and multi-view stereo (SfM and MVS) algorithms, obtaining 3D point clouds of 1000-20,000 points m −2 , resulting in a multispectral orthomosaic of good accuracy.As part of the processing, individual trees were manually geolocated on the orthomosaic, and delimitation of the crown area was complemented with the ForestTools algorithm of the statistical software R Studio, version 4.3.1 (The R Project for Statistical Computing, Vienna, Austria).This tool also helps to delimit the tree crowns through the variable window filter algorithm.Using the zone statistics tool of QGIS version 3.16.15,descriptive NDVI statistics for each tree of the genera P. engelmannii, Q. grisea, and J. deppeana were extracted from the center of each canopy, with the mean selected as the variable of interest.

Figure 3 .
Figure 3. (a) Climate-radial growth correlations in the studied species.The months of the previous and current year are represented by lower-and upper-case letters, respectively.Horizontal lines denote the p < 0.05 significance level.(b) Monthly boxplot of NDVI at the species level for the annual period of the flight.The values correspond to the NDVI calculated from flights conducted from August to December 2021 and from January to July 2022.

Figure 3 .
Figure 3. (a) Climate-radial growth correlations in the studied species.The months of the previous and current year are represented by lower-and upper-case letters, respectively.Horizontal lines denote the p < 0.05 significance level.(b) Monthly boxplot of NDVI at the species level for the annual period of the flight.The values correspond to the NDVI calculated from flights conducted from August to December 2021 and from January to July 2022.

Figure A1 .
Figure A1.Residual chronologies of the total ring width indices of the studied species.

Figure A2 .
Figure A2.Pearson correlations (r) between the standardized precipitation and evapotranspiration index (SPEI) and residual tree-ring width chronologies for each species: (a) J. deppeana, (b) P. engelmannii, and (c) Q. grisea.SPEI was obtained on a scale from 1 to 24 months.

Figure A1 . 20 Figure A1 .
Figure A1.Residual chronologies of the total ring width indices of the studied species.

Figure A2 .
Figure A2.Pearson correlations (r) between the standardized precipitation and evapotranspiration index (SPEI) and residual tree-ring width chronologies for each species: (a) J. deppeana, (b) P. engelmannii, and (c) Q. grisea.SPEI was obtained on a scale from 1 to 24 months.

Figure A2 .
Figure A2.Pearson correlations (r) between the standardized precipitation and evapotranspiration index (SPEI) and residual tree-ring width chronologies for each species: (a) J. deppeana, (b) P. engelmannii, and (c) Q. grisea.SPEI was obtained on a scale from 1 to 24 months.

Figure A3 .
Figure A3.Graphs of the goodness of fit of the model.(a) Linearity, (b)Homogeneity of variance and (c) Normality of residuals.The green line denotes a normal distribution.

Figure A3 .
Figure A3.Graphs of the goodness of fit of the model.(a) Linearity, (b)Homogeneity of variance and (c) Normality of residuals.The green line denotes a normal distribution.

Figure A4 .
Figure A4.Post hoc test plot per species.Values correspond to NDVI calculated from flights conducted in the years 2021-2022.Green, pink, and blue colors correspond to Q. grisea, P. engelmannii, and J. deppeana, respectively.Rectangles indicate those months in which there were no significant differences between all the genera studied.

Figure A4 .
Figure A4.Post hoc test plot per species.Values correspond to NDVI calculated from flights conducted in the years 2021-2022.Green, pink, and blue colors correspond to Q. grisea, P. engelmannii, and J. deppeana, respectively.Rectangles indicate those months in which there were no significant differences between all the genera studied.

Figure A5 .
Figure A5.Diagram of NDVI variations per genus throughout the sampling year.Green, pink, and blue colors correspond to the genera Q. grisea, P. engelmannii, and J. deppeana, respectively.

Figure A5 .
Figure A5.Diagram of NDVI variations per genus throughout the sampling year.Green, pink, and blue colors correspond to the genera Q. grisea, P. engelmannii, and J. deppeana, respectively.

Table 1 .
Dasometric statistics of the study species for the dendrochronological sampling.

Table 2 .
Dendrochronological statistics of the studied trees.
Values are means ± SE.Abbreviations: AC = first-order autocorrelation; Rbar = mean correlation among trees; and EPS = expressed population signal.

Table 3 .
Descriptive statistics of NDVI per species on the different sampling dates in the study area.
N = number of trees, Mean = mean NDVI value, SD = standard deviation, Min = minimum, Max = maximum.

Table 4 .
Descriptive statistics of dasometric variables per species in the study site.

Table 5 .
ANOVA of the monthly NDVI and per species in the study site.