Leaﬁng Patterns and Drivers across Seasonally Dry Tropical Communities

: Investigating the timing of key phenological events across environments with variable seasonality is crucial to understand the drivers of ecosystem dynamics. Leaf production in the tropics is mainly constrained by water and light availability. Identifying the factors regulating leaf phenology patterns allows e ﬃ ciently forecasting of climate change impacts. We conducted a novel phenological monitoring study across four Neotropical vegetation sites using leaf phenology time series obtained from digital repeated photographs (phenocameras). Seasonality di ﬀ ered among sites, from very seasonally dry climate in the caatinga dry scrubland with an eight-month long dry season to the less restrictive Cerrado vegetation with a six-month dry season. To unravel the main drivers of leaf phenology and understand how they inﬂuence seasonal dynamics (represented by the green color channel (Gcc) vegetation index), we applied Generalized Additive Mixed Models (GAMMs) to estimate the growing seasons, using water deﬁcit and day length as covariates. Our results indicated that plant-water relationships are more important in the caatinga, while light (measured as day-length) was more relevant in explaining leaﬁng patterns in Cerrado communities. Leaﬁng behaviors and predictor-response relationships (distinct smooth functions) were more variable at the less seasonal Cerrado sites, suggesting that di ﬀ erent life-forms (grasses, herbs, shrubs, and trees) are capable of overcoming drought through speciﬁc phenological strategies and associated functional traits, such as deep root systems in trees.


Introduction
Temporal patterns of leaf replacement, or leaf exchange phenology, are of major importance in understanding ecosystem processes such as carbon, water, and energy exchanges, controlling seasonal cycles of photosynthetic activity [1][2][3]. A high heterogeneity of leafing patterns has been observed for the tropics, as climate seasonality varies widely in the intensity and length of the dry season [1,4]. Identifying the main factors that regulate leaf phenology in the tropics is thus crucial to understand ecosystem dynamics and to efficiently forecast climate change impacts on these increasingly threatened environments.
Water and light availability are considered the main abiotic cues regulating the timing and periodicity of leaf production in tropical vegetations [4][5][6][7][8][9]. For instance, water availability regulates the length of the growing season and the phenological synchrony among cerrado species [4]. Seasonal communities with longer dry seasons exhibit limited photosynthetic activity and a marked annual periodicity of leaf flushing and senescence patterns, favoring a larger proportion of deciduous species [1,[10][11][12]. In fact, within seasonally dry tropical forests, nearly all plant individuals lose leaves in a synchronic deciduous behavior during the dry season [1,[13][14][15]. In areas with less pronounced dry seasons and elevated total annual rainfall, species and individuals may display different degrees of deciduousness, establishing communities with a wide range of leafing behaviors that may change according to soil moisture conditions, topography [8,16], and the intensity of the dry season ( [4] and references therein).
In the tropics, variation in day-length, and consequently in available light, has major impacts towards higher latitudes, triggering the input of new leaves [6][7][8]. In tropical seasonal environments such as semi-deciduous forests and savannas day-length triggers early leaf flushing, anticipating the onset of the rainy season [4,8,[17][18][19]. Plant lifeforms also have direct implications in determining species adaptations to water and light availability, and thereby on leafing phenology [19][20][21]. Trees can produce new leaves during the dry season due to their capacity to reach accumulated groundwater sources (deeper rooting systems), allowing plant growth when photosynthetically active radiation is at its maximum [4,20,22,23]. Contrarily, given their shallow root system and architectural constraints, grasses are highly dependent on precipitation and rapidly respond to rainfall events [20,24].
Large scale studies (e.g., [12,25]) have contributed to the understanding of ecosystem productivity constraints worldwide and land surface phenology has been an important tool to determine seasonality of leaf reflectance. Nonetheless, these studies consider general vegetation cover types at broader spatial scales and use summarized data sets (monthly temporal scales) and often overlook the need for ground validation and the pursue of processes and mechanisms underlying phenology at finer spatial and temporal scales [26,27].
Long-term phenological time series in the tropics are scarce, stimulating the search for reliable plant monitoring methods [28,29]. Imagery based on digital cameras, or "phenocams" has been considered an alternative method to successfully and continuously monitor plant greening across landscapes [30,31]. Phenocams can track temporal changes in leafing patterns across tropical seasonal environments [32,33], while reducing human labor, increasing accuracy by eliminating observer subjectivity, improving temporal resolution to an hourly/daily basis, and improved spatial resolution, which allow simultaneous multiple-sites monitoring [30,[33][34][35].
Here, we conducted a novel phenological multi-site near-surface monitoring study using digital cameras to describe community leafing patterns in four Neotropical seasonally dry environments that differ in their degree of seasonality and vegetation structure: a seasonally dry tropical forest (SDTF), locally known as caatinga, and three different Cerrado (Neotropical savanna) vegetation types-two woody savannas (cerrado and dense cerrado) and a grassy shrubland savanna (cerrado shrubland). We combined time series of leaf phenology obtained by phenocams and local environmental variables observed at a fine temporal scales (daily) to describe the patterns of leaf flushing and senescence and to evaluate the constraints imposed by water, light, and temperature in all four vegetation types. Specifically, we investigated: (i) whether the timing and length of the growing season differs across the four seasonally dry tropical communities and what are the drivers of phenological transitions in each environment; and (ii) whether the growing season varies according to the vegetation structure (woody versus grassy) and degree of seasonality, determined by the water deficit observed during the dry season. We expected the growing season to differ between Cerrado and caatinga sites, with the leafing season responding to rainfall events in the caatinga, while Cerrado communities should be influenced by daylength (light availability) due to the milder and shorter dry season. We also expected leafing patterns of Cerrado communities to differ between woody and grassy vegetation types according to [36], while the leafing pattern for the caatinga community should be restricted to the onset of the rainy season.

Study Sites
Sites were geographically distributed across two main vegetation domains [37] or biomes [38]: the Caatinga, or xeric shrubland, and the Cerrado, which includes grasslands, savannas, and shrublands ( Figure 1). We studied four different sites: one from the caatinga and three different Cerrado vegetation types. Vegetation communities differ in their seasonality: the caatinga site has a very constraint climate with an 8-month dry season, while the Cerrado sites show less restricted water conditions, with a 6-month dry season. Sites also differed in regard to vegetation structure, from a dominant herbaceous component with scattered shrubs (caatinga and cerrado shrubland) to woody formations of cerrado sensu stricto (woody cerrado and dense cerrado). To easily identify the four study sites and vegetation types, we will hereafter refer to them as follows: caatinga, cerrado shrubland, woody cerrado, and dense cerrado. Our study sites can be sorted, according to their degree of seasonality and woody cover, from caatinga to cerrado shrubland to woody cerrado to dense cerrado (Table 1).

Figure 1.
Geographic distribution of study sites, with biome distribution (background color) following [38], and the vegetation types monitored in this study: caatinga, cerrado shrubland, woody cerrado and dense cerrado. Geographic distribution of study sites, with biome distribution (background color) following [38], and the vegetation types monitored in this study: caatinga, cerrado shrubland, woody cerrado and dense cerrado. Our first site represents a xeric, sclerophyllous, seasonally dry vegetation located in the semi-arid region in northeast Brazil ( Figure 1). It may be denominated as a "xeric shrubland biome" [38] or considered as part of the largest continuous seasonally dry tropical forest (STDF) region in the Americas but is locally known as caatinga [39]. The study area comprises approximately 600 ha, 390 m a.s.l, and belong to the Reserva Legal da Embrapa Semiárido, northeast Brazil [40]. The climate is classified as semiarid [41], with mean annual temperature of 26.2 • C (according to the experimental station of Bebedouro, 10 km from the site). Based on climatic data from 1970 to 2014, total mean annual precipitation is 510 mm, distributed mainly from January to April. However, during our three-year study, mean annual precipitation was circa 260 mm and mean annual temperature was 27 • C (Figure 2b,c) characterizing a particularly dry period. Local vegetation was composed of xerophilous trees and shrubs and a herbaceous layer of species adapted to xeric conditions. The most common families were Fabaceae, Euphorbiaceae, Poaceae, and Cactaceae, and the discontinuous canopy reaches a height of circa 5 m [40].

Cerrado
The other three sites are considered part of the Cerrado biome/domain, or Cerrado sensu lato [42], and also included in the Grassland, Savanna, and Shrubland biome as proposed by [38] (Figure 1). All three sites are located in São Paulo, southeastern Brazil: two in the municipality of Itirapina and one in Santa Rita do Passa Quatro ( Figure 1). The three sites are under the same climatic regime, classified as humid subtropical, with a cool dry winter and a hot wet summer [41], with a 6 month dry season from April to September (mean monthly rainfall below 60 mm [5,7]) and wet season from October to March (Table 2).

•
Cerrado shrubland-the cerrado shrubland (campo sujo) is characterized by a dominant grassy herbaceous layer with scattered shrubs and small trees [24,42]. This site is part of the Itirapina Ecological Station, with an area of 2300 ha and an altitude of 700 m. Previous vegetation surveys found that 79% of the species from the herbaceous-shrubland layer belong to Asteraceae, Fabaceae, Poaceae, and Cyperaceae, and 21% of the small tree species were mainly from Fabaceae, Myrtaceae, and Melastomataceae [43]. As reported by [44], local historic climatic data show a total mean annual precipitation of 1524 mm and mean temperature of 20.7 • C. During the three years we monitored the area, the mean annual precipitation was 1272 mm and a mean annual temperature 23.5 • C (Figure 2e,f). Historic climatic information was made available by the Centro de Recursos Hídricos e Estudos Ambientais (CRHEA-EESC/USP), located about 15 km from the study sites.
• Woody cerrado-this site is located in a private area of 260 hectares and 700 m a.s.l. adjacent to the cerrado shrubland site. This woody cerrado is a remnant of the original Cerrado that once covered the entire region, and therefore, faces the same climatic regime as the cerrado shrubland site. The vegetation is a cerrado sensu stricto, which is characterized by a dominant woody layer (trees and shrubs ranging from 3 to 12 m in height) with discontinuous crown cover and a scattered herbaceous layer [44]. According to local meteorological data (2011-2015), mean annual total precipitation was 1478 mm and mean annual temperature 22.9 • C (Figure 2h,i). The most common plant families are Myrtaceae, Fabaceae, and Malpighiaceae [44], and the vegetation has been classified as semi-deciduous based on long-term leaf exchange strategies [4].

•
Dense cerrado-this site is a dense woody cerrado that belongs to the Reserva Ecológica Pé-de-Gigante (PEG), located within the Parque Estadual do Vassununga, in Santa Rita do Passa Quatro, SP. The PEG reserve comprises a continuous area of 1060 ha and 649 m a.s.l. covered by a heterogeneous landscape of Cerrado vegetation types, from open grasslands to dense woody cerrado. As reported by Batalha et al. [24], local climatic data (1986-1995) presented a total mean annual rainfall of 1499 mm and mean annual temperature of 21.5 • C. Our local climatic data (2013 to 2015) showed a total mean annual precipitation of 1150 mm and a mean annual temperature of 22.5 • C (Figure 2k,l). The study site is a transition from woody cerrado to cerradão [42], which we classified as dense cerrado; it is characterized by a less discontinuous top canopy, nearly no herbaceous layer, and a high density of shrubs and trees [45]. Dominant woody layer reaches 10 to 15 m height is composed mainly of Ptedoron pubecens, Copaifera langsdofii, and Anadenanthera peregrina var. falcata (all Fabaceae) [46].

Near-Surface Phenology Monitoring
For each study site, a digital hemispherical Mobotix Q 24 lens camera (Mobotix AG-Germany) was placed in a tower attached to an extension arm facing northeast (east for the caatinga site) at a mean vertical distance of 10 m (6 m for the caatinga site) from the top of the canopy. At the caatinga and dense cerrado study sites, the local flux towers measure 16 and 26 m, respectively; the phenological tower at woody cerrado measures 18 meters. At the cerrado shrubland, the camera system was placed at 4 m from the vegetation in a landscape standing point [see 33]. The cameras were supplied by a 12 V battery that was charged by solar panels. Cameras were set to automatically take a sequence of five JPEG images (at 1280 × 960 pixels of resolution) in the first 10 minutes of each hour, from 6:00 to 18:00 h (UC−3; Universal Time Coordinated), on a daily basis, as described in Alberton et al. [32]. The camera systems were installed at different times during the monitored years (see Table 1), so we analyzed the period ranging from 2013 to 2015 for all sites except the woody cerrado, where the time series started in 2011 (Table 1). Due to several environment-related energy supply issues throughout the monitoring years for the woody cerrado site, we developed an algorithm to fill in gaps of more than seven days with no recorded images, using the R programming language [47] (see Supplementary material for more details).
Raw images were initially screened visually to remove photographs where canopy view was obstructed. The remaining images were analyzed as described by Richardson et al. [48], Ahrends et al. [49], and Alberton et al. [32]. We carried out field expeditions to identify plant species and compile a list of the main species captured by the field of view of each camera (see Supplementary material). We then matched the identified individuals with the crowns observed in the images. This information, associated with visual image inspections, was essential to properly classify the individual crowns into plant functional groups based on leaf exchange strategies and to confirm the phenological transition dates.
Regions of interest (ROIs) were defined in the images of each site to capture the community as well as separate deciduous species (for cerrado areas only; based on [4]), and grass versus woody cover (for cerrado shrubland only). The community ROIs refer to the entire community present in the image (ROI = entire image excluding the tower and bare soil patches); deciduous ROIs refer to all deciduous crowns observed in the image (cerrado shrubland = 16.16%; woody cerrado = 5.44%; dense cerrado = 21.5%); grass ROIs refer to the herbaceous cover observed in the cerrado shrubland (23.70%). The proportion of deciduous species and grass cover was determined by the percentage of pixels in relation to the total pixels from the community ROI.
We analyzed all ROIs in terms of the contribution of the relative brightness of red, green, and blue color channels [RGB chromatic coordinates; 50]. The normalized RGB chromatic coordinate (RGBcc) index is the most suitable index to detect leaf color changes, such as leaf development and

Near-Surface Phenology Monitoring
For each study site, a digital hemispherical Mobotix Q 24 lens camera (Mobotix AG-Germany) was placed in a tower attached to an extension arm facing northeast (east for the caatinga site) at a mean vertical distance of 10 m (6 m for the caatinga site) from the top of the canopy. At the caatinga and dense cerrado study sites, the local flux towers measure 16 and 26 m, respectively; the phenological tower at woody cerrado measures 18 m. At the cerrado shrubland, the camera system was placed at 4 m from the vegetation in a landscape standing point (see [33]). The cameras were supplied by a 12 V battery that was charged by solar panels. Cameras were set to automatically take a sequence of five JPEG images (at 1280 × 960 pixels of resolution) in the first 10 min of each h, from 6:00 to 18:00 h (UC−3; Universal Time Coordinated), on a daily basis, as described in Alberton et al. [32]. The camera systems were installed at different times during the monitored years (see Table 1), so we analyzed the period ranging from 2013 to 2015 for all sites except the woody cerrado, where the time series started in 2011 (Table 1). Due to several environment-related energy supply issues throughout the monitoring years for the woody cerrado site, we developed an algorithm to fill in gaps of more than seven days with no recorded images, using the R programming language [47] (see Supplementary material for more details).
Raw images were initially screened visually to remove photographs where canopy view was obstructed. The remaining images were analyzed as described by Richardson et al. [48], Ahrends et al. [49], and Alberton et al. [32]. We carried out field expeditions to identify plant species and compile a list of the main species captured by the field of view of each camera (see Supplementary material). We then matched the identified individuals with the crowns observed in the images. This information, associated with visual image inspections, was essential to properly classify the individual crowns into plant functional groups based on leaf exchange strategies and to confirm the phenological transition dates.
Regions of interest (ROIs) were defined in the images of each site to capture the community as well as separate deciduous species (for cerrado areas only; based on [4]), and grass versus woody cover (for cerrado shrubland only). The community ROIs refer to the entire community present in the image (ROI = entire image excluding the tower and bare soil patches); deciduous ROIs refer to all deciduous crowns observed in the image (cerrado shrubland = 16.16%; woody cerrado = 5.44%; dense cerrado = 21.5%); grass ROIs refer to the herbaceous cover observed in the cerrado shrubland (23.70%). The proportion of deciduous species and grass cover was determined by the percentage of pixels in relation to the total pixels from the community ROI.
We analyzed all ROIs in terms of the contribution of the relative brightness of red, green, and blue color channels (RGB chromatic coordinates; [50]). The normalized RGB chromatic coordinate (RGB cc ) index is the most suitable index to detect leaf color changes, such as leaf development and photosynthetic dynamics, and the most efficient index to suppress light variations [32,[48][49][50][51][52][53][54][55]. We thus calculated the normalized index of the green color channel (Gcc) according to the formulas: The Gcc index was calculated separately for each of the five hourly images taken each day. One single measurement was then extracted by taking the 90th percentile of all values observed every three days, to minimize interferences related to lighting changes caused by weather, season, and time of day [55].

Defining the Growing Season
We extracted phenological metrics to define the growing season and to detect phenological transitions [56]. Time series of tropical communities are often more difficult to characterize temporally, as they are mostly rain-green systems [56]. To overcome problems related to the detection of green-up responses coupled with precipitation events, which result in noisier observations, we propose, for the first time, the use confidence intervals of curve derivatives to identify changes in phenology. To do so, we first fitted Generalized Additive Mixed Models (GAMM) to the Gcc community time series (see [4]). Time (sequence of days) was used as an independent smoother variable. Assuming that errors are not independent, a condition inherent to the structure of time series, we assumed an auto-regressive moving average (ARMA) correlation to our model residuals. We then generated 10,000 independent simulations of fitted curve and calculated the derivative of each simulated curve at each time step, thus building confidence intervals for the rate of changes (derivatives significantly different from zero) along the curve. A significance level of 10% was chosen due to the high variability associated with daily phenocam observation. From the derivatives, we were able to extract phenological metrics, or phenological transitions, that corresponded to significant periods of increases in greenness values, followed by a short period of no change, and then a significant period of decrease in greenness. Based on previous classifications of land surface phenology time series [57,58] and according to the traditional classification of directly observed phenological variables [7], we defined the phenological metrics as: Start of growing season (SOS)-represents the beginning of the growing season and is measured as the first day detected by a significant derivative from the total seasonal amplitude on the left side of the curve.
End of growing season (EOS)-represents the end of the growing season and is measured as the last day detected by a significant derivative from the total seasonal amplitude on the right side of the curve.
Length of growing season (LOS)-represents the duration of the growing season and is calculated as the difference between SOS and EOS.
Peak of growing season (POS)-represents the highest greenness percentage, generally associated to the highest percentage of species flushing new leaves. It is measured as the highest value of the seasonal curve.

Environmental Cues of Leafing Phenology
Photoperiod, temperature, water availability, and water-energy interactions [59] were chosen as the potential environmental factors to characterize leafing phenological patterns across sites. The set of explanatory variables was selected based on a priori assumptions and on information extracted from the literature.
At the caatinga and dense cerrado study sites, environmental data was measured with sensors installed in the local flux towers where phenocams were set up. A CR1000 datalogger (Campbell Scientific Inc.) was used to calculate averages of 30-s measurements at each 30-min interval between May 2013 and December 2015 for the caatinga site and between August 2013 and December 2015 in the dense cerrado site. For the other locations, a U30-Hobo meteorological station, installed in the phenological tower located at each site, collected all environmental variables.
Maximum and minimum daily temperatures were used to calculate the temperature amplitude along the year: Rainfall effects were tested using the following variables: cumulative precipitation, from the last 30 days (Rainfall_cum_30), and lagged rainfall events from the last seven (Rainfall_Lag7) and 30 (Rainfall_Lag30) days.
Water and energy interactions were calculated to evaluate the effects of the dry season intensity using the cumulative water deficit (cwd) according to James et al. [60]: where cwdn is the computed cumulative water deficit for a given time (n), Pn is the precipitation of the day (mm), and ETn is the evapotranspiration of the day (mm day −1 ). When a positive value of (P − ET) is identified, zero is used to set a starting point. Every time cwdn ≥ 0, then cdwn is set to zero again (more details in [59][60][61]). Interactions between evapotranspiration demand, which is induced by the amount of energy in the system (able to heat and/or stress the plant), and water availability in the soil are explained by the climatic water balance (see [59]). When evaporative demand is higher than the amount of water available, a deficit occurs, which explains the cwd (a measure of absolute drought) [59][60][61]. Evapotranspiration was calculated by the eddy covariance method in the caatinga and dense cerrado sites. For cerrado shrubland and woody cerrado, we used the evapotranspiration value of the CRHEA-EESC/USP meteorological station (located 10 km away from the sites). Day-length (the number of sunlight hours) was chosen to represent photoperiod seasonality and was calculated using the latitude of each location (geosphere package in R). The environmental variables were aggregated into three-day windows using the same percentile approach applied to the greenness series [55].

Data Analysis
Phenological metrics were extracted from the community Gcc time series and compared among vegetation types. To investigate the relationship between leafing phenology and environmental cues, we fitted separate Generalized Additive Mixed Models (GAMMs) using the community Gcc time series, leaf exchange strategy, and dominant life-forms as dependent variables (Y) and the abiotic environmental data as explanatory variables (X). We used GAMM as they are able to capture nonlinear relationships across multiple covariates, reliably handle structural time series in ecological studies and vegetation modeling [62][63][64].
The general formula of a GAM is: Where g is a specified link function, (µ) the expected mean response, and f(X) are smooth functions of the covariates (X) [65]. Basically, GAMs aim to maximize the quality of the prediction (µ) by estimating non-parametric functions of the predictor variables (X) through connections with the dependent variable that are established by a link function (g) [63,66]. The Generalized Additive Mixed Model (GAMM) approach adds a random factor into the general GAM formula. Since we are dealing with time series (Gcc) and assuming that errors are not independent, an auto-regressive moving average (ARMA) correlation structure was nested into the models.
Initially, scatterplots among covariates were verified to identify linear and non-linear relationships. After running multicollinearity inspection, we used a full model for each dependent variable where we entered all chosen explanatory variables. Variable selection, which means a smooth component selection when working with GAMMs, was carried using a shrinkage approach, denominated double penalty, which relies on the fact that the space of the smoother (a spline basis) can be decomposed into components: range space and null space [66]. The method penalizes both components shrinking them to zero. By introducing a penalty in the null space, the GAMM smoothing estimator can select the most significant terms from the model (for more details about the method, see [66]). By analyzing the approximate smooth estimators from the model, we can identify terms that are not statistically significant (significance with p values < 0.001) and remove them altogether. Finally, to understand how each term is interacting with the response variable, we interpreted plots of the partial fits of each smoother parameter and evaluated F-tests of model outcomes.

Detecting Growing Seasons Across Sites
The Gcc index calculated for the community time series tracked the seasonality of leafing patterns in all sites ( Figure 3). Using the derivatives, we were able to find the timing of leafing phenological transitions for all annual cycles ( Table 2). In general, growing seasons differed among vegetation types mainly in relation to the start and end dates, season length, and shape of the seasonal curves, which are related to the characteristic leafing response of each community. smooth component selection when working with GAMMs, was carried using a shrinkage approach, denominated double penalty, which relies on the fact that the space of the smoother (a spline basis) can be decomposed into components: range space and null space [66]. The method penalizes both components shrinking them to zero. By introducing a penalty in the null space, the GAMM smoothing estimator can select the most significant terms from the model [for more details about the method, see 66]. By analyzing the approximate smooth estimators from the model, we can identify terms that are not statistically significant (significance with p values < 0.001) and remove them altogether. Finally, to understand how each term is interacting with the response variable, we interpreted plots of the partial fits of each smoother parameter and evaluated F-tests of model outcomes.

Detecting Growing Seasons Across Sites
The Gcc index calculated for the community time series tracked the seasonality of leafing patterns in all sites ( Figure 3). Using the derivatives, we were able to find the timing of leafing phenological transitions for all annual cycles ( Table 2). In general, growing seasons differed among vegetation types mainly in relation to the start and end dates, season length, and shape of the seasonal curves, which are related to the characteristic leafing response of each community.  The start of the growing season (SOS) occurred earlier for cerrado vegetation types in comparison to the caatinga, and ranged from day of the year (DOY) 147 to DOY 239 ( Table 2). The woody cerrado had the earliest SOS (mean DOY = 163), occurring mid-June (mid-dry season). SOS in the dense cerrado and cerrado shrubland were similar and occurred in August (mean DOY 221 and 233, respectively). SOS in the caatinga occurred by the end of October (mean DOY 296), after the first rainfall events. All communities reached the peak of their growing season (POS) about 78 days after their SOS. Caatinga and dense cerrado were the first to reach their POS (69 and 67 days, respectively), while the cerrado shrubland and the woody cerrado reached their POS later (80 to 96 days, respectively). For the two wooded cerrado vegetation types, POS occurred during the dry-to-rainy season transition, contrary to the caatinga and cerrado shrubland, where POS dates occurred within the rainy season.
The end of the growing season (EOS) occurred faster and earlier in the caatinga vegetation, and the difference between the first (2013/2014 = 183) and the second (2014/2015 = 214) cycles was of about a month (31 days; Table 2, Figure 3). The length of the growing season (LOS) was also the shortest in the caatinga, with a mean LOS of 274 days. Additionally, green-down was more abrupt and EOS occurred immediately after the end of the rainy season (Figure 3a). Conversely, cerrado vegetation types showed a slower transition and extended green-down curve, reaching the lowest values from the mid to the end of the dry season (Figure 3b-d). The woody cerrado showed different EOS across years, ranging from mid-May (DOY 134 in 2012/21013) to mid-June (DOY 170 in 2011/2012), and LOS also varied from 314 to 342 days. Only one LOS was recorded for the cerrado shrubland and the dense cerrado sites (a single complete cycle), and both sites had the longest LOS, with 345 days.

Model Predictions
Predictor variables were selected as smooth terms in most of the models. Two of the explanatory variables had linear relationships when modeled against the ROI data: Rainfall_lag7 in the caatinga and Tamp in the woody cerrado. In general, the GAMMs produced models with medium to high explanation power (R 2 from 0.29 to 0.88). Day-length and cwd were the most recurrent predictors, with the highest statistical significance values across all vegetation types (Table 3). However, other variables were also significant in explaining leaf phenology according to the GAMMs (Figures 3 and 4). The most important relationships are described in more detail below. Table 3. Approximate significance of the smooth terms used in the GAMMs. Effective degrees of freedom (edf) and F-test values are given for each of the following variables: day-length, cumulative water deficit (cwd), precipitation lagged seven days (Raifall_lag7), cumulative precipitation from the last 30 days (Rainfall_cum_30), and thermal amplitude (Tamp). The coefficient of determination (R 2 ) of each model is also shown. All variables were significant, with p values < 0.001. comu = community; dec = deciduous; grass = grassy layer. ROI: region of interest. For caatinga, the variables that best explained community greening (Gcc) were cwd and rainfall (lag7 and cum30) (Figure 4a,b). The partial fit between predictors and Gcc demonstrated that vegetation was sensitive to low values of precipitation and to cwd values between −10 and −5, which caused an increase in Gcc values (Figure 4b). Day-length, also included in the model as a significant variable (Table 3), started to affect the community Gcc after 12.5 h of sunlight. For the three cerrado vegetation sites, day-length was the variable that best explained community Gcc patterns ( Table 3). The same values, ranging from 11.5 to 12.0 hours, were observed as a sensitive threshold for phenological transitions across all vegetation types (Figure 4e, 4h, and 4K). Cwd resulted in different smooth functions among cerrado sites, where cerrado shrubland and dense cerrado were strongly sensitive to cwd with green-down values between -150 and -100 ( Figure  4f and 4l), while woody cerrado was less sensitive (Figure 4i).

Site
Cumulative rainfall (Rainfall_cum_30) was significant for both cerrado shrubland and dense cerrado ( Table 3), indicating that changes in Gcc can occur in both vegetation types when values of accumulated precipitaton rates are between 25 to 30 mm (Figure 4g and 4m). The Gcc time series from the woody cerrado community was the only one that interacted with thermal amplitude (Tamp-3.558, 17.05, P < 0.001, Table 3). There was a positive linear relationship between the covariates, where Gcc values increased as Tamp increased ( Figure 4J).
The ROIs separating deciduous leaf exchange strategies and life forms (grassy layer) interacted with different variables across vegetation types. We compared the deciduous ROIs among sites, but since the overall caatinga community displays a deciduous behavior, we used the community ROI for the caatinga (and not the deciduous ROI) to carry out such comparisons. The relationship between the community Gcc patterns and precipitation events were higher in the caatinga, followed by the cerrado shrubland, woody cerrado, and, finally, dense cerrado ( Figure 5). Particularly, the caatinga vegetation showed a pattern of fast response to water availability, with pulses of green-up following rainfall events (Figure 5a). Among the cerrado vegetation types, the Gcc curve of the cerrado shrubland was better linked to precipitation in comparison to curves from the woody and dense cerrados, mainly when considering the Gcc curve oscillations (Figure 5b to 5d).
The deciduous ROI curves from the cerrado sites showed variable patterns. The deciduous shrub layer from the cerrado shrubland followed a pattern similar to the community pattern, with an overlap of green-up and green-down periods (Figure 5b). Conversely, the deciduous ROI curve from the woody cerrado showed green-up within the rainy season, with high values of Gcc still during the For the three cerrado vegetation sites, day-length was the variable that best explained community Gcc patterns ( Table 3). The same values, ranging from 11.5 to 12.0 h, were observed as a sensitive threshold for phenological transitions across all vegetation types (Figure 4e,h,k). Cwd resulted in different smooth functions among cerrado sites, where cerrado shrubland and dense cerrado were strongly sensitive to cwd with green-down values between −150 and −100 (Figure 4f,l), while woody cerrado was less sensitive (Figure 4i).
Cumulative rainfall (Rainfall_cum_30) was significant for both cerrado shrubland and dense cerrado ( Table 3), indicating that changes in Gcc can occur in both vegetation types when values of accumulated precipitaton rates are between 25 to 30 mm (Figure 4g,m). The Gcc time series from the woody cerrado community was the only one that interacted with thermal amplitude (Tamp-3.558, 17.05, P < 0.001, Table 3). There was a positive linear relationship between the covariates, where Gcc values increased as Tamp increased (Figure 4j).
The ROIs separating deciduous leaf exchange strategies and life forms (grassy layer) interacted with different variables across vegetation types. We compared the deciduous ROIs among sites, but since the overall caatinga community displays a deciduous behavior, we used the community ROI for the caatinga (and not the deciduous ROI) to carry out such comparisons. The relationship between the community Gcc patterns and precipitation events were higher in the caatinga, followed by the cerrado shrubland, woody cerrado, and, finally, dense cerrado ( Figure 5). Particularly, the caatinga vegetation showed a pattern of fast response to water availability, with pulses of green-up following rainfall events (Figure 5a). Among the cerrado vegetation types, the Gcc curve of the cerrado shrubland was better linked to precipitation in comparison to curves from the woody and dense cerrados, mainly when considering the Gcc curve oscillations (Figure 5b-d).
The deciduous ROI curves from the cerrado sites showed variable patterns. The deciduous shrub layer from the cerrado shrubland followed a pattern similar to the community pattern, with an overlap of green-up and green-down periods (Figure 5b). Conversely, the deciduous ROI curve from the woody cerrado showed green-up within the rainy season, with high values of Gcc still during the dry season (Figure 5c). In the dense cerrado, green-up occurred during the dry season before the first rainfall events, reaching its peak in the transition between the dry and rainy seasons (Figure 5d). The grass ROIs had green-up curves similar to the community and to the deciduous species but peaked during the rainy season. The grassy component also showed marked fluctuations within the rainy season, apparently following individual rainfall events, even though we observed the highest peak during the dry season (Figure 5b).  (Figure 5c). In the dense cerrado, green-up occurred during the dry season before the first rainfall events, reaching its peak in the transition between the dry and rainy seasons (Figure 5d). The grass ROIs had green-up curves similar to the community and to the deciduous species but peaked during the rainy season. The grassy component also showed marked fluctuations within the rainy season, apparently following individual rainfall events, even though we observed the highest peak during the dry season (Figure 5b). Among the cerrado vegetation types, we observed that leafing strategy interacted with different environmental variables. In the cerrado shrubland, cwd and cumulative rainfall (Rainfall_Cum_30) interacted exclusively with the green-up pattern of the deciduous ROI (Figure 6c,d). Cwd values between −100 and −150 mm led to green-down events, and green-up events only began after cwd values over −100 mm. The variables day-length and Tamp also significantly affected the green-up and green-down patterns. Tamp between 5 and 10 °C increased Gcc, which decreased after a Tamp over 10 °C. (Figure 6a,b). In the woody cerrado, the variables cwd and Tamp significantly affected the deciduous ROI model, with green-up being positively affected by the decrease in water defict (below −100 mm) and negatively affected by Tamp over 10 °C (Figure 6e,f). In the dense cerrado, the curve inflection of deciduous ROI was sensitive to day-lengths between 12.0 to 12.5 hours (Figure 6g). Cumulative rainfall (Rainfall_Cum_30) and cwd showed peaks of increasing green-up when values ranged from 15 to 30 mm and −75 to −50 mm, respectively (Figure 6h,j).
The grass layer from the cerrado shrubland was significantly related to day-length, where a positive linear relationship was established with day-lenth over 11.5 hours, while with Tamp the curve showed an inverse sigmoid relationship (Figure 6a,b). Among the cerrado vegetation types, we observed that leafing strategy interacted with different environmental variables. In the cerrado shrubland, cwd and cumulative rainfall (Rainfall_Cum_30) interacted exclusively with the green-up pattern of the deciduous ROI (Figure 6c,d). Cwd values between −100 and −150 mm led to green-down events, and green-up events only began after cwd values over −100 mm. The variables day-length and Tamp also significantly affected the green-up and green-down patterns. Tamp between 5 and 10 • C increased Gcc, which decreased after a Tamp over 10 • C. (Figure 6a,b). In the woody cerrado, the variables cwd and Tamp significantly affected the deciduous ROI model, with green-up being positively affected by the decrease in water defict (below −100 mm) and negatively affected by Tamp over 10 • C (Figure 6e,f). In the dense cerrado, the curve inflection of deciduous ROI was sensitive to day-lengths between 12.0 to 12.5 h (Figure 6g). Cumulative rainfall (Rainfall_Cum_30) and cwd showed peaks of increasing green-up when values ranged from 15 to 30 mm and −75 to −50 mm, respectively (Figure 6h,j).
The grass layer from the cerrado shrubland was significantly related to day-length, where a positive linear relationship was established with day-lenth over 11.5 h, while with Tamp the curve showed an inverse sigmoid relationship (Figure 6a,b). Remote Sens. 2019, 11, x FOR PEER REVIEW 13 of 22

Communities Leafing Patterns
The Gcc index extracted from the digital images allowed us to track leafing seasonality in the studied tropical seasonal communities and determine the growing season using phenological metrics detected from the derivatives of fitted generalized additive mixed models. As expected, Gcc curves coincided with the rainy season, leading to differences in growing season and phenological transitions among vegetations, especially between the caatinga and the cerrado vegetation types. For the woody and dense cerrado, the beginning of the growing season (SOS) and their peak (POS) occurred during the dry season, while for the caatinga and the cerrado shrubland, leaf flushing began with the first rainfall events and reached its peak within the rainy season. According to our results, environmental responses of community leafing patterns change gradually across these vegetation types, with the caatinga as the most water-responsive vegetation, followed by the cerrado shrubland, which shows important day-length and water-related responses, and finally both wooded cerrado communities, where day-length was the main variable affecting the onset of the growing season.
The shape of growing season curves differed among vegetation types and were clearly influenced by water availability throughout the year, leading to a variation in the length of the growing season (LOS): much shorter in the caatinga in comparison to the cerrado sites. Green-up and green-down curves were more abrupt for the caatinga, which may be explained by the predominance and synchronicity of seasonally deciduous species that occur in these environments [10,65]. Conversely, the mild dry season observed for the Cerrado favors a different proportion of deciduous, semi-deciduous, and evergreen species. Green cover in Cerrado is maintained even during the dry season and Gcc curves change throughout the rainy season. Our findings support the direct observations of leaf exchange patterns that have been previously determined at the woody cerrado site and the classification of the Cerrado vegetation as a semi-deciduous vegetation [4,67]. Community green-down patterns for the shrubland and woody cerrado sites started at the end of the rainy season (represented by the last Gcc oscillation before the descending pattern of the terminal derivative), while in the dense cerrado green-down occurred from the second Gcc oscillation after the peak. Such decreasing values of Gcc (from the peak date until the beginning of senescence) correspond to the period of leaf maturation. Thus, Gcc index is sensitive to leaf color changes [68] and there is a gradual color change with the ageing of community leaves (see Methods).

Communities Leafing Patterns
The Gcc index extracted from the digital images allowed us to track leafing seasonality in the studied tropical seasonal communities and determine the growing season using phenological metrics detected from the derivatives of fitted generalized additive mixed models. As expected, Gcc curves coincided with the rainy season, leading to differences in growing season and phenological transitions among vegetations, especially between the caatinga and the cerrado vegetation types. For the woody and dense cerrado, the beginning of the growing season (SOS) and their peak (POS) occurred during the dry season, while for the caatinga and the cerrado shrubland, leaf flushing began with the first rainfall events and reached its peak within the rainy season. According to our results, environmental responses of community leafing patterns change gradually across these vegetation types, with the caatinga as the most water-responsive vegetation, followed by the cerrado shrubland, which shows important day-length and water-related responses, and finally both wooded cerrado communities, where day-length was the main variable affecting the onset of the growing season.
The shape of growing season curves differed among vegetation types and were clearly influenced by water availability throughout the year, leading to a variation in the length of the growing season (LOS): much shorter in the caatinga in comparison to the cerrado sites. Green-up and green-down curves were more abrupt for the caatinga, which may be explained by the predominance and synchronicity of seasonally deciduous species that occur in these environments [10,65]. Conversely, the mild dry season observed for the Cerrado favors a different proportion of deciduous, semi-deciduous, and evergreen species. Green cover in Cerrado is maintained even during the dry season and Gcc curves change throughout the rainy season. Our findings support the direct observations of leaf exchange patterns that have been previously determined at the woody cerrado site and the classification of the Cerrado vegetation as a semi-deciduous vegetation [4,67]. Community green-down patterns for the shrubland and woody cerrado sites started at the end of the rainy season (represented by the last Gcc oscillation before the descending pattern of the terminal derivative), while in the dense cerrado green-down occurred from the second Gcc oscillation after the peak. Such decreasing values of Gcc (from the peak date until the beginning of senescence) correspond to the period of leaf maturation. Thus, Gcc index is sensitive to leaf color changes [68] and there is a gradual color change with the ageing of community leaves (see Methods).
The interannual variation in SOS and EOS dates was higher in the caatinga (SOS differed 28 days between years and EOS differed 31 days) and in the woody cerrado (SOS differed 38 days between years and EOS differed 36 days). For the caatinga, the Atlantic Dipole phenomenon (a temperature variation of the Atlantic Ocean) causes high interannual variability in precipitation rates [69], which would lead to shifts in seasonality dynamics of the growing season among years. According to our models and given that plant species from the caatinga are influenced mainly by water availability, there is an immediate response in both reproductive and vegetative phenology after rainfall events, given these events determine the onset of the rainy season [70,71]. Regarding the woody cerrado, studies usually report that leaf flushing (SOS) begins at the end of the dry season or during the dry-to-rainy season transition [4,11,[72][73][74]. Among woody cerrado species, our results indicate two main leaf flushing strategies: species either respond to the first rains that occur during the dry-to-rainy season transition or leaf flushing begins still during the dry season (see [4] for further discussions). Similarly, leaf flushing of shrubland cerrado species began during the dry season; however, their POS dates differed from the woody cerrado community as it occurred during the rainy season, probably due to the dominant grass layer.
Our model predicted non-linear relationships between phenology and the recurrent significant variables: day-length and cwd. The skewness pattern of greening in the caatinga indicates there is a link between time and duration of the irregular leafing cycles, which are restricted to the rainy season and follow the precipitation pulses that occur in the system. Thresholds in plant-water relationships indicate that Gcc changes are stimulated by low rates of accumulated precipitation, which explains the fast response of leaf flushing observed for caatinga species.
Day-length was the most common variable explaining Gcc patterns for cerrado communities. Increasing day-length (longer days) triggers leaf flushing of most cerrado species in the communities studied, starting early in the dry season. Photoperiod has been reported to trigger leaf flushing in seasonal dry tropical vegetation [8,19,[75][76][77], but studies in the Brazilian caatinga are still lacking [78].
The community leaf flushing peak occurring before the onset of the rainy season supports the hypothesis that cerrado woody species have access to underground water, due to deep root systems, or are able to rehydrate themselves by using water from leaves of the previous cycle [79,80]. Consistent photoperiod signals have been suggested to trigger and synchronize leaf flushing in the tropics, even though species show different leaf exchange strategies [4,8,76,81] even in non-seasonal systems [7]. Community leaf flushing before the rainy season would also avoid the higher herbivory rates that occur during the rainy season and the leaching-induced nutrient loss during leaf development, given nutrient uptake would take place before the rains [4,82,83].
The cwd relationships also explained the seasonality of leaf production in all four seasonally dry tropical communities investigated in this study. Its relationship with the Gcc time series in the caatinga was subtle when compared to the relationships observed in the cerrado sites, which corroborates that caatinga species are sensitive to the slightest increase in water availability. Cwd showed similar patterns across the three cerrado sites, with green-down occurring when water deficits were high (−150 mm to −100 mm) and green-up with cwd around −50 mm. Cumulative water deficits are essential in understanding plant physiological responses given that it integrates water and energy interactions [59]. Additionally, cwd is related to drought stress experienced during the dry period [60,61]. Global scale studies have proposed a threshold of 2000 mm of annual rainfall in tropical forests in order to sustain water demands and keep the evergreen state of species throughout the year [12,84]. In Cerrado species, the dry season has been stated as a key factor triggering vegetative phenological transitions, affecting leaf abscission and leaf emergence [85].
The 30-day cumulative rainfall (Rainfall_cum_30), which was present in all models, explained the community leafing patterns in all sites except for the woody cerrado. In comparison with the caatinga, the cerrado sites needed twice as much cumulative rainfall to stimulate increases in Gcc. These results indicate the role water availability plays in leafing patterns of cerrado species, where day-length acts as a first trigger, stimulating leaf unfolding, and rainfall intensity acts as a driver of leaf expansion, thus avoiding drought limitations during drier periods [86,87]. Contrarily, in the woody cerrado, thermal amplitude (Tamp) was the variable that best explained the community Gcc. Tamp was also positively correlated with rates of green-up and senescence in plant communities across a snow-free mountain, which included grasslands and cerrado areas [88]. At the woody cerrado site, Tamp reached 20 • C during the dry-to-wet season transition and was also related to earlier patterns of leaf flushing.
Anticipating leaf flushing may provide major carbon gains, since species with early leaf emergence increase carbon assimilation rates once the rainy season begins [85,89,90], resulting in a longer and more productive growing season.

Deciduous Leaf Exchange Strategies and Life Forms
The deciduous component showed different patterns across the four communities: in the caatinga, all species were deciduous, while the cerrado shrubland was composed of a deciduous woody component and a grassy layer and, in the two wooded cerrado sites, the deciduous component decreased-deciduous ROI cover was around 6% at the woody cerrado and 20% at the dense cerrado. The proportion of deciduousness is reflected in the degree of seasonality of the community Gcc time series across sites, from the marked seasonal caatinga, followed by the cerrado shrubland with a deciduous woody component to the woody cerrado. The factors affecting the proportion of the deciduous component coincide with the results found for the entire community, where woody cerrado species are less water-dependent in comparison with caatinga and cerrado shrubland species.
The caatinga was the most synchronous plant community, with all deciduous species responding simultaneously to cumulative rainfall. Such high synchronicity suggests that a common climatic factor triggers leaf phenological transitions, which could be explained by the strong water-response relationship. Opportunistic drought deciduous species (ODDS) are mainly driven by water availability with leaf flushing following the first rainfall events [91]. Seasonally harsh climatic conditions, characteristic of the caatinga, suggests a dominance of ODDS, with synchronic leaf flushing occurring after the first rainfall events at the end of the dry season [10,11,16,91]. The interannual variation of rainfall amount and distribution likely determines community leaf flushing and its intensity [92].
However, plant species from relatively more moderate environments, such as the cerrado, may be able to cope with water limitations and produce new leaves still during the dry season [93]. For instance, the deciduous leafing behavior from the woody cerrado was more affected by day-length than water availability. Cerrado trees use deep roots to exploit deeper soil water sources, which enables them to flush new leaves before the rainy season [23,76,94]. Schedule drought deciduous species (SDDS) are those able to endure drought conditions by maximizing their carbon gain without experiencing stress and are mainly driven by photoperiod, with high synchronicity and low interannual variability [4,22,91,95]. The persistent greening during the dry season indicates a longer resistance to drought limitations and a dominance of SDDS, even though the main drivers for deciduousness in the woody cerrado were cwd and Tamp.
In the cerrado shrubland, leaf phenology of the grassy layer was not predicted by water relationships, contradicting our expectations regarding the temporal niche separation hypothesis [36] and the water constraints imposed on phenological responses among grass species [43,96]. The grass ROI model was partially (50%) explained by day-length and Tamp. The reponse suggests the influence of invasive African grasses, particularly Melinis minutiflora and Urochloa brizantha, which are extremely abundant in the area [43,97]. The isolated greening event that occurred during the dry season, observed in both yearly cycles, could be an indication of the effect invasive species had on our results. The anticipated green-up event may have caused the low correlation with rainfall covariates from the model, as leafing of invasive species began before the rainy season and, thus, before native graminoids. The anticipation suggests the priority effects hypothesis, which states that an earlier establishment in the growing season benefits invasive species as they will use resources first [98], and explains the phenological behavior of such species [98]. Additionally, the rapid evolution of phenology may be a widespread and essential process during the expansion of many invasive species [99]. However, studies comparing the phenological behavior of invasive and native grasses are still lacking, hindering adequate management actions.
Our results demonstrated that phenological leaf exchange dynamics in seasonally constrained environments tend to be more synchronic, since caatinga and grass species are adapted to flush during a limited but favorable period for plant development (rainy season). As the environment becomes less constrained by climate (mild dry season), leaf exchange patterns become less restricted to the rainy season [100], allowing leaf flushing to begin in the dry season.

Conclusions
We found that water and light are the most important predictors for leaf phenological patterns across seasonally dry environments (caatinga and cerrado). For caatinga, short cumulative seven-day rainfall events were able to initiate the growing season and leaf flushing. Day-length and cwd were the most influential factors across cerrado sites, regardless of leaf exchange strategies. Degree of seasonality, which is related to the length and intensity of the dry season, reflected on community leafing patterns and on what drives the deciduous component of these communities. While deciduous species from the xeric caatinga were responsive to cumulative water availability, the woody deciduous species of cerrado were responsive to day-length. The absence of rainfall covariates in the grass model suggests the interference of invasive grasses. Understanding leaf phenology mechanisms are essential to effectively forecast future Earth systems. The results reported here contribute to detect leaf flushing and senescence transitions, their main drivers, and the thresholds affecting them. Thus, our findings contribute to understand the dynamics of tropical seasonal communities across areas with varying water dependence. The effects of anthropogenic climate changes in the tropics are predicted to result in longer and more intense dry seasons [89,90]. The increase of aridity during the dry season has been reported to cause greater impacts over productivity than the excess of water, seasonally or annually, could cause [61,101]. The impact of more intense dry seasons would be bigger in deciduous species as they are more sensitive to seasonal climatic variations. A major impact is expected for the deciduous caatinga species given the strong water-dependence. Cerrado species, triggered by day-length, would continue to flush new leaves in the dry season; however, dryness could affect the full development of leaves and shorten the growing season, which would anticipate senescence. Non-deciduous species, a significant part of cerrado communities [4], would also be affected since the increase of aridity reduces C gain, affecting productivity [61,83,101].

Conflicts of Interest:
The authors declare no conflict of interest.