Topography and Traits Modulate Tree Performance and Drought Response in a Tropical Forest

Predicting drought responses of individual trees in tropical forests remains challenging, in part because trees experience drought differently depending on their position in spatially heterogeneous environments. Specifically, topography and the competitive environment can influence the severity of water stress experienced by individual trees, leading to individual-level variation in drought impacts. A drought in 2015 in Puerto Rico provided the opportunity to assess how drought response varies with topography and neighborhood crowding in a tropical forest. In this study, we integrated 3 years of annual census data from the El Yunque Chronosequence plots with measurements of functional traits and LiDAR-derived metrics of microsite topography. We fit hierarchical Bayesian models to examine how drought, microtopography, and neighborhood crowding influence individual tree growth and survival, and the role functional traits play in mediating species’ responses to these drivers. We found that while growth was lower during the drought year, drought had no effect on survival, suggesting that these forests are fairly resilient to a single-year drought. However, growth response to drought, as well as average growth and survival, varied with topography: tree growth in valley-like microsites was more negatively affected by drought, and survival was lower on steeper slopes while growth was higher in valleys. Neighborhood crowding reduced growth and increased survival, but these effects did not vary between drought/non-drought years. Functional traits provided some insight into mechanisms by which drought and topography affected growth and survival. For example, trees with high specific leaf area grew more slowly on steeper slopes, and high wood density trees were less sensitive to drought. However, the relationships between functional traits and response to drought and topography were weak overall. Species sorting across microtopography may drive observed relationships between average performance, drought response, and topography. Our results suggest that understanding species’ responses to drought requires consideration of the microenvironments in which they grow. Complex interactions between regional climate, topography, and traits underlie individual and species variation in drought response.


INTRODUCTION
Tropical rainfall regimes are predicted to change in future climate scenarios, with many parts of the tropics becoming drier (Chadwick et al., 2015;Duffy et al., 2015). Dry conditions will likely have large impacts on tropical forests: drought influences forest ecosystem structure, composition, and function (Bonal et al., 2016;Uriarte et al., 2016b), and importantly, could decrease the size of the tropical forest carbon sink (Phillips et al., 2009;Pan et al., 2011;Gatti et al., 2014). However, large uncertainties about the impacts of drought on tropical forests remain, in part due to the difficulties of experimentally manipulating moisture conditions in tropical forests (but see Nepstad et al., 2007;da Costa et al., 2010). Observational studies of forest dynamics during natural droughts provide an opportunity to learn how drought affects tropical forests, especially where long-term data has been collected over multiple years (e.g., Condit et al., 1995;Clark et al., 2010;Uriarte et al., 2016a). Understanding the impacts of recent droughts will help anticipate future changes in tropical forests caused by shifting frequency and intensity of precipitation.
Most studies of drought in tropical forests have aimed to quantify effects on carbon uptake and storage and to understand how sensitivity to drought varies across size classes and tree species. Drought increases tree mortality and reduces growth which can result in large losses of stored carbon from tropical forests (Chazdon et al., 2005;Nepstad et al., 2007;da Costa et al., 2010;Phillips et al., 2010;Gatti et al., 2014). Growth and mortality of larger trees tend to respond more strongly to drought, though this relationship is not consistent across sites (Nepstad et al., 2007;Phillips et al., 2010;Bennett et al., 2015;Powers et al., 2020). Sensitivity to drought varies drastically across species. These differences can often be explained by physiological or functional traits, such as wood density, specific leaf area, stem hydraulic conductivity, turgor loss point, rooting depth, or specific leaf area, though the strength and direction of trait relationships with drought are highly variable (Skelton et al., 2015;Uriarte et al., 2016a;Greenwood et al., 2017;O'Brien et al., 2017;Powers et al., 2020).
However, even within species and size classes, substantial unexplained variation in individual drought response exists (Condit et al., 1995;Clark et al., 2010;Phillips et al., 2010). Spatial heterogeneity in the severity of water stress experienced by individual trees linked to topography, soils, or competitive environment could explain these differences. For example, competition for water via neighborhood crowding can suppress tree growth and drive variation in mortality rates across rainfall gradients (Ruiz-Benito et al., 2013). A more intense competitive environment can amplify water shortages during drought, increasing the probability of dying (Gleason et al., 2017).
Similarly, moisture conditions are well understood to vary across topographic positions. Variation in drainage and runoff means that slopes and ridges are drier than valleys (Burt and Butcher, 1985;Western et al., 1999;Daws et al., 2002). In the northern hemisphere, southwest facing slopes receive more solar radiation and have higher rates of evapotranspiration, and so water stress is typically higher (Stephenson, 1990). Accordingly, drought-induced mortality is often higher in drier landscape positions (Fekedulegn et al., 2003;Guarín and Taylor, 2005). In wetter landscape positions, moist microsites might provide refugia for less drought-tolerant species, mitigating drought impacts (Dobrowski, 2011). In future drier climates, these sites could rescue species that might face extinction if subject to regional climates (McLaughlin et al., 2017).
Despite its likely importance, few studies have considered topographic variation in drought effects in tropical forests (but see Nakagawa et al., 2000;Silva et al., 2013;Schwartz et al., 2019). Previous research has shown that topography influences species distributions (Ashton et al., 2006;Engelbrecht et al., 2007;Bartlett et al., 2016), canopy dynamics and disturbance rates (Comita and Engelbrecht, 2009;Silva et al., 2013). Drought could amplify or diminish these differences in performance, due to moisture stress being more severe at drier topographic positions. Furthermore, variation in performance across topography should depend on species and their functional traits, as some species are more sensitive to moisture, light, or nutrient availability than others.
Studies linking tree performance, moisture, and topography have typically focused on variation across sites or plots (Fekedulegn et al., 2003;Guarín and Taylor, 2005;Ashton et al., 2006;Engelbrecht et al., 2007;Comita and Engelbrecht, 2009). However, soil moisture and tree performance both vary with microtopography, i.e., topographic relief at very small (1 m) scales (Famiglietti et al., 1998;Daws et al., 2002;Tenenbaum et al., 2006). Ecological studies of the effects of microsite topography on tree performance have typically used categorical designations of topographic position (e.g., hummock vs. hollow, Nishimua et al., 2007;dry plateau vs. wet slopes, Engelbrecht et al., 2007), but mapping and identifying these categories across a continuous landscape is difficult. New remote sensing techniques can generate digital elevation models (DEMs) at very fine scales (<1 m), which can be used to quantify microtopographic relief and can be linked to variation in soil moisture (Tenenbaum et al., 2006;Buchanan et al., 2014). To our knowledge, no study has linked these quantitative characterizations of microsite topography with variation in tree performance during drought in the tropics.
In this study, we leverage a recent drought in Puerto Rico, and link annual tree census data from second-growth forests with LiDAR-derived measures of microtopographic relief to assess how drought affects tree growth and survival, and how traits, topography, and crowding mediate drought response. Specifically, we ask: (1) What was the effect of a drought on tree growth and survival? We expected growth and survival to be reduced during drought, and for these effects to be more severe for larger trees.
(2) How does microtopographic relief affect tree growth and survival, and do its effects differ depending on rainfall? We expected that growth and survival would be lower on steeper slopes and ridges (i.e., areas with more convex curvature). Because these topographic positions tend to be drier, we also expected the effects of drought on tree growth and mortality would be amplified on steeper slopes and more convex surfaces.
(3) How does crowding affect tree growth and survival, and do its effects vary depending on rainfall? We expected that crowding would reduce growth and survival, and further predicted that the effect of crowding would be amplified during dry years, due to more intense competition for water between neighbors.
(4) How does interspecific variation in functional traits mediate the effects of rainfall, microtopography, and crowding on tree demographics? We expected that trees with traits representing more resource acquisitive strategies would have higher growth rates and lower rates of survival than trees with more conservative traits. Furthermore, we predicted that trees with acquisitive traits would be more sensitive to stressful conditions, i.e., drought, dry topographic position, and crowding.

MATERIALS AND METHODS
To address our questions, we used hierarchical Bayesian models that allowed us to examine the importance of drought, topography, and crowding for tree performance. We focused on two topographic variables that are important for moisture conditions and flow of water across surfaces: slope and curvature (Burt and Butcher, 1985). We used this analytical approach to assess if and how functional traits drive inter-specific variation in response to environmental conditions. We considered two functional traits that are associated with carbon metabolism and plant hydraulics: specific leaf area (SLA) and wood density (WD). SLA represents the investment in photosynthetic machinery (leaf surface area) relative to total investment in leaf biomass. Leaves with high SLA tend to have higher photosynthetic rates and nutrient concentrations, and high SLA species typically have higher growth rates. However, high SLA leaves tend to be shorterlived, potentially leading to shorter full-plant life spans (Reich et al., 1998;Santiago and Wright, 2007). Wood density represents the investment in wood biomass per volume of wood. Trees with higher WD are typically more shade tolerant (Valladares and Niinemets, 2008), and tend to have lower growth rates (Poorter et al., 2008). However, trees with denser wood are more resistant to cavitation (Hacke et al., 2001;Markesteijn et al., 2011) and structural damage (Everham and Brokaw, 1996;Curran et al., 2008), so species with higher wood density tend to have longer life spans, higher survival rates, and lower sensitivity to drought (Poorter et al., 2008;Phillips et al., 2010;Greenwood et al., 2017, but see Hoffmann et al., 2011.

Study Area, Tree Census, and Rainfall Data
This study was conducted with data from four forest plots, comprising the El Yunque Chronosequence Plots ( Table 1). The land-use histories and ages of these forests were determined from aerial photographs taken between 1936 and 1977 (Lugo et al., 2004): three of the plots were previously cleared for agriculture and represent a range of forest ages from 35 to 76 years since agricultural activities ceased. The fourth plot was, to our knowledge, never cleared for agriculture, but has been subject to hurricane disturbance. The plots range in elevation from 100 to 500 m above sea level, and vary in size from ∼0.5 to 1 hectare ( Table 1). The dominant soil type is Zarzal, a deep and well-drained oxisol (Mount and Lynn, 2004). All stems >1 cm diameter at breast height (dbh) have been mapped, and identified to species, and diameters and tree status (live/dead) have been recorded in September of every year since 2012. We used this data to calculate absolute diameter growth and survival for each individual tree for each census interval from 2013 to 2016.
In 2015, Puerto Rico experienced a severe meteorological drought: rainfall in El Yunque was only 2,035 mm, the second lowest recorded since 1975 (Figure 1). Annual rainfall in the region averages 3,500 mm, and rainfall was close to average in 2014 and 2016 (3,193 and 3,506 mm, respectively). For the purposes of assessing drought effects, we used growth and survival data from the 2014 (2013-2014), 2015 (2014-2015), and 2016 (2015-2016) censuses in our analyses. We used daily rainfall data from the El Verde field station to calculate the total rainfall between census dates. Rainfall since the last census was used as a predictor in all models to assess the effect of drought on tree growth and survival.

Functional Trait Data
Wood density and SLA measurements were collected on 5-10 individuals per species using standard protocols (Cornelissen et al., 2003;Swenson et al., 2012), with stem wood density for shrubs estimated from branch material following the methods outlined in Swenson and Enquist (2008). For all analyses, we used the mean trait value for each species, though we acknowledge that intraspecific trait variation may influence performance along environmental gradients (Bolnick et al., 2011). The two traits were weakly correlated (Pearson's r = −0.35, p = 0.01).

Topography Data
Topography data was derived from airborne LiDAR data, collected by the National Center for Airborne Laser Mapping in May 2011 1 . We followed standard procedures to generate a digital elevation model (DEM) at 1 m 2 resolution from LiDAR returns, using the minimum z-values of the last-return ground classified points to construct the DEM. Further details about the LiDAR data and DEM construction are in Wolf et al. (2016).
We used the 1 m LiDAR DEM to derive topographic slope and hilltop curvature, following methods in Hurst et al. (2012). Using a very fine-scale DEM to calculate topographic indices can result in local extreme values, and many studies coarsen the resolution of these DEMs to minimize the impact of these extremes (e.g., Jucker et al., 2018). Instead of coarsening the DEM's resolution, the Hurst et al. (2012) method fits a surface model to each focal grid cell based on the surrounding grid cells across a moving window of a specified length. This smooths the surface, reducing local extremes but preserving the 1-m resolution of the DEM. Slope and curvature are then calculated from the fit coefficients of the model. We used a 99 × 99 m moving window to fit this regression, such that the surface is fitted for each 1 m 2 grid cell of the DEM, but taking into account a 99 m neighborhood. This scale best fits soil moisture data collected at a nearby site (Supplementary Figure S1). Slope (m/m) and curvature (unitless) are calculated from the fitted coefficients following the methods in Hurst et al. (2012). Hilltop curvature was calculated such that positive curvature indicates valleylike microtopography and negative curvature indicates ridge-like microtopography. We used the georeferenced stem locations from the plot data to extract slope and curvature at the stem location for each tree. Slope and curvature in this dataset were not correlated (r = −0.004). Though the mean slope and curvature differ across plots, plots encompass large, overlapping ranges of values for slope and curvature (Supplementary Figure S2).

Modeling Approach
We fit hierarchical Bayesian models of annual diameter growth and survival as a function of rainfall, topography, and crowding, with a second level in the model describing the relationship between functional traits and the species-specific responses to these factors. Growth was normally distributed, as negative growth is common due to shrinkage and/or measurement error. Our model of the expected value of growth took the form: where g tsi is absolute diameter growth of individual i of species s at time t. Covariates include a binary indicator for drought/nondrought year (drought t = 0 in 2014 and 2016, drought t = 1 in 2015), slope and curvature at the precise location of each stem (slope si and curvature si ), and NCI tsi , the neighborhood crowding index, a measure of crowding. γ 1i is an individual random effect for each stem. Stem diameter (DBH tsi ) was included as a control, to account for well-known variation in growth and survival with tree size. To account for differences across plots, we included a random effect for plot (γ 2p ). Both random effects were drawn from a normal distribution with mean 0 and standard deviations σ i and σ p , respectively. Error (ε) was normally distributed. DBH, NCI, and slope were highly left-skewed and therefore log-transformed to facilitate analysis. DBH, NCI, slope, and curvature were not strongly correlated (all r < 0.16, Supplementary Table S2). We modeled the logittransformed probability of survival (p a ) with the same covariates and random effects as the growth model, as well as a predictor term for each stem's previous year's growth (g t −1,si ). Survival observations were Bernoulli distributed with probability p a . NCI is a dimensionless quantity calculated for each stem, taking into account the diameter and distance of all stems within a 10 m radius around the focal tree. Specifically, it is calculated as: Frontiers in Forests and Global Change | www.frontiersin.org where stem i has j neighbors within 10 m and d ij is the distance from stem i to its neighbor j. We used a 10 m radius as prior studies have indicated that this radius is sufficient to capture effects of crowding (Uriarte et al., 2004a). Excluding trees less than 10 m from the edge of the plot would have resulted in exclusion of a large number of individuals. For those edge trees, we scaled their NCI by the ratio of a full-size neighborhood (i.e., a 10 m radius circle) to the size of the edge tree's partial neighborhood.
We incorporated functional traits into the second level of our models to assess how interspecific trait variation mediates the effects of drought, topography, and crowding on tree demography. If functional traits capture variation in plant strategies, then species-level responses to stressful or highresources conditions should vary predictably with their trait values. We expected that functional traits might influence average growth and survival rates (β 1s ) along with species' sensitivities to drought (β 3s ), crowding (β 4s ), topography (β 5s , β 6s ), and their interactions (β 8s , β 9s , β 10s ). We did not model species-specific parameters for tree diameter (β 2 ) or the interaction between drought and diameter (β 7 ). For each covariate (k) we modeled the species-specific β ks as a normally distributed process deriving from a linear function of that species' traits: σ k is the variance in species' responses to covariate k unexplained by functional traits. Because it was strongly left-skewed, we logtransformed SLA. We centered and scaled the traits so that b k represents the mean species response to covariate k, and b k1 and b k2 represent the departure from the mean with an increase of one standard deviation of log(SLA) and WD, respectively.
In the models of both growth and survival, we included data for all species with more than 20 individuals and with available trait data, which accounted for 49% of all species, and 94% of all basal area for a total of 46 species (Supplementary Table S1). All stems were included in calculations of neighborhood crowding. The tree species Syzygium jambos (Myrtaceae) was excluded from analyses, as it is an exotic species that in recent years has been experiencing extremely high mortality due to infection by the rust Austropuccinia psidii (Burman et al., 2017). In exploratory analyses it appeared to be a major outlier and had a disproportionate influence on the magnitude and direction of relationships between the predictors and response variables, so we removed it from analyses presented here. We fit all models including and excluding the palm species Prestoea montana, as this species is often an outlier (Muscarella et al., 2013). Exclusion of P. montana did not qualitatively change the results, and here we present results including the palm species. To facilitate convergence and remove outliers, we excluded growth and NCI observations greater than two standard deviations from the mean, due to these two variables having several extreme, unrealistic values in the dataset. This resulted in 24,218 growth observations and 25,125 survival observations across all years.
We standardized all predictors by subtracting the mean and dividing by the standard deviation to facilitate model convergence and ease interpretation (Gelman and Hill, 2006).
We standardized DBH and NCI on a species-by-species basis, to avoid confounding their effects with species-specific differences in size and crowding. Other predictors were standardized across the whole dataset. We specified uninformative priors for all parameters, and estimated posterior distributions for all parameters using Markov chain Monte Carlo (MCMC) sampling implemented in JAGS (Plummer, 2003). We verified convergence visually and by ensuring the potential scale reduction statistic ( R) was equal to 1 (Gelman and Rubin, 1992). Models generally converged after 30,000 iterations. We performed posterior predictive checks by simulating predicted growth and survival for all observations (Supplementary Figures S3, S4) and calculating the R 2 and root mean squared error (rmse) between predicted and observed values. We considered β parameters to be statistically significant if their 95% credible interval did not overlap 0, and as marginally significant if their 90% credible interval did not overlap 0. All statistical analyses were conducted in R (R Core Team, 2016) with the packages rjags and R2jags (Plummer, 2016;Su and Yajima, 2018).

Growth Model
Across the study period, the average annual growth was 0.06 cm.  Figure S4). Diameter was the strongest predictor of growth, with larger trees growing more (Figure 2A and Supplementary Table S3). After diameter, drought was the most important predictor of growth. The interaction between drought and diameter was significant, such that growth of larger trees was more negatively affected by drought.
Topography and crowding variables were also significant predictors of tree growth, though less important for drought response. More crowded individuals (trees with higher NCI) had reduced growth (Figures 2A, 3C). There was no significant relationship between slope and growth. Curvature had a marginally significant positive association with growth, showing that during normal, non-drought years, trees in local valleys grow more. Furthermore, there was a significant interaction between curvature and drought such that that trees located in areas with higher curvature (i.e., valley-like microsites) grew more slowly during drought (Figures 2A, 3B). There was no significant interaction between crowding and drought, or between slope and drought, indicating that the effects of crowding and slope on growth were similar in drought vs. non-drought years (Figures 2A, 3A,C). To summarize, though crowding and slope both affected growth, they were not important for drought response, and only curvature influenced drought response.
Functional traits were correlated with species' growth responses to drought and topography, though there were  (Eq. 3) i.e., average species effects. Parameters with credible intervals that did not overlap 0 were considered to be significant. no significant trait associations with sensitivity of growth to crowding ( Table 2). Trees with dense wood had significantly less negative growth responses to drought ( Figure 4A). SLA was negatively associated with slope response, such that trees with low SLA experience a positive effect of slope on growth, whereas trees with high SLA experience a negative effect of slope ( Figure 4B). SLA also had a negative association with the curvature-drought interaction term, such that growth of species with high SLA, when located in sites with valley-like microtopography, were more negatively affected by drought (Table 2 and Figure 5).

Survival Model
The mean survival rate across the whole dataset was 93.1%, with significant inter-annual variation (91.7% in 2014, SD = 0.28; 92.5% in 2015, SD = 0.26; 95.1% in 2016, SD = 0.22; p < 0.001). The model of survival captured some observed variation (R 2 = 0.32, Supplementary Figure S5). Diameter and antecedent growth were the most important predictors of survival, with larger trees and trees that had experienced higher growth in the preceding year having a higher probability of survival. Drought did not affect survival, nor was there a significant interaction between drought and diameter (Figures 2B, 3). More crowded trees had a higher probability of survival, though this effect was not significant (Figures 2B, 3F). Survival was negatively associated with slope, meaning that trees on steeper slopes had a lower probability of survival (Figures 2B, 3D). Curvature did not have a significant effect on survival (Figures 2B, 3E). In the survival model none of the interaction terms were significant, indicating that the effects of antecedent growth, microtopography, and crowding on survival were similar, on average, across all 3 years of the study (Figures 3D-F).
There were two significant trait associations with first-level parameters in the survival model (Table 2 and Figures 4C,D). Trees with high SLA had lower average survival (i.e., the intercept term, β 1s ), and responded more positively to antecedent growth ( Table 2). Wood density was not significantly associated with any parameters in the survival model, and we found no significant trait associations with the effects of drought, microtopography, crowding, or their interactions on survival.

DISCUSSION
Drought can have large effects on tropical forests, but impacts of drought vary across species and space (Bonal et al., 2016). Spatial variability in drought impacts is often driven by differences in the local moisture conditions that individual trees experience, because of the way regional climate is filtered through topography, vegetation, and other environmental factors that vary on small scales (McLaughlin et al., 2017), and by species differences in physiology and drought response. The 2015 drought in Puerto Rico provided a unique opportunity to leverage repeat census data collected in a topographically complex tropical forest to examine the way that species differences and local filtering of climate affect drought response. Our results show that (1) some tropical forests may be resilient to a single year drought, as drought negatively affected growth but not survival, (2) microtopography and crowding influence tree performance, but their interactions with drought may not conform with predictions, and (3) functional traits can provide limited insight into the mechanisms by which drought, microtopography, and crowding affect growth and survival.

Drought Effects on Tree Performance
As expected, drought affected tree performance both directly and via interactions with microtopography. Though growth was lower in the drought year, drought did not affect survival directly or via interactions with microtopography or crowding. This lack of immediate drought-induced mortality is at odds with observations of high mortality during single-year droughts in some tropical forests (Phillips et al., 2009;Lewis et al., 2011;Doughty et al., 2015), but consistent with results from some experimental studies: in two drought experiments in the Amazon, 50-60% of rainfall was excluded, but mortality rates were low during the experiments' first years and major die-offs were only observed after about 3 years of sustained drought (Nepstad et al., 2007;da Costa et al., 2010). In 2015, El Yunque rainfall was only 56% of the annual mean, but the drought lasted for less than 1 year (Mote et al., 2017), perhaps not long enough to cause a severe die-off. However, antecedent growth was a strong predictor of survival, and drought did have a strong effect on growth, implying that survival could be reduced after a lag. Lagged effects of drought on growth and mortality lasting for multiple years have been observed in many ecosystems (Anderegg et al., 2015) including tropical forests (Lingenfelder and Newbery, 2009;Phillips et al., 2010). In our study a lag effect was not observable immediately after the drought: survival was higher in the year following the drought. Still, lag effects can last multiple years, and more years of study could reveal further effects of drought.
Tropical tree growth seems to respond to drought on more rapid time scales than survival: in the same throughfall experiments described above, growth impacts were observed in the experiments' first years (Brando et al., 2008;da Costa et al., 2010). Tree-ring and long-term monitoring studies have found strong correlations between diameter growth and rainfall in tropical trees (Brienen et al., 2010;Clark et al., 2010), and have found that growth of large trees is more severely affected by drought (Bennett et al., 2015). Accordingly, drought was the strongest predictor of growth in our model, and large trees suffered larger growth reductions during drought. Reductions in radial growth may be driven by overall declines in productivity, and/or shifts in allocation from stem growth to leaves, branches, roots, or non-structural carbohydrates (Malhi et al., 2015). Trees can use non-structural carbohydrates to maintain core metabolic  These parameters represent the trait relationships with species-specific responses to each covariate. Values provided are median estimates with 95% credible intervals in parentheses. * indicates significance (i.e., 95% credible interval did not overlap with 0), whileˆindicates marginal significance (90% credible interval does not overlap 0). NA indicates parameter was not included in model.
functions when photosynthetic rates are reduced during drought, and shift allocation during and after drought (Metcalfe et al., 2010;Doughty et al., 2014). The drivers of tree allocation of carbon are poorly understood, though they are essential for understanding the mechanisms by which drought affects trees (Malhi et al., 2015).

Topography and Crowding Effects on Tree Performance
Effects of microtopography on tree performance were evident, but did not always conform to expectations. For example, we expected that tree growth and survival would both be lower on steeper slopes. We found that survival declined with increasing slope, but that slope had no effect on growth. Contrary to our expectation that drought would further amplify variation in performance with slope, we found that the effects of slope did not vary across drought and non-drought years. This result could be explained by several possible mechanisms. First, moisture differences across slopes may be consistent regardless of rainfall. However, studies have often found that the strength and direction of the relationship between moisture and topography differ depending on rainfall (Famiglietti et al., 1998;Tenenbaum et al., 2006;Tromp-van Meerveld and McDonnell, 2006). Another possible explanation for why drought effects did not vary across slopes is that trees growing on steeper slopes could be more drought tolerant. Even if drought is more severe on steeper slopes, differences in tree performance might not be amplified if trees growing on steeper slopes are more drought tolerant. In such a case, despite experiencing more severe moisture deficits, species in dry topographic positions would not suffer more severe impacts on performance. Studies in other tropical forests have shown that local variation in water availability affects species distributions, with droughttolerant species showing a stronger affinity for dry micro-sites (Ashton et al., 2006;Engelbrecht et al., 2007;Comita and Engelbrecht, 2009), though these studies did not link this species sorting to variation in demographic rates. Our results provide some evidence for species sorting across microtopography via weak, but significant correlations between functional traits and topographic indices (Supplementary Table S2), though most species span a wide range of microtopography (Supplementary Figure S6). Trees growing at drier topographic positions may also be better acclimated to dry conditions due to individuallevel environmental plasticity, which could be reflected through intra-specific trait variation or variation in allocation to above vs. belowground structures (Nicotra et al., 2010).
Finally, topographic variation in soil depth, nutrient availability, or susceptibility to disturbance, and not in soil moisture, could drive differences in performance across slope. Soil nutrient availability is known to vary with microtopography and fine-scale variation in nutrient variability affects tree species distributions, diversity, and performance (Baraloto and Couteron, 2010;Homeier et al., 2010;Baldeck et al., 2013). If variation in nutrient availability were the key driver of differences in performance across slopes, we would not expect drought to influence variation in performance across microtopography. The elevated mortality rates on slopes may also reflect higher rates of wind disturbance, because trees on steeper slopes are more susceptible to windthrow from extreme winds (Everham and Brokaw, 1996). In general, adult trees' deep and extensive roots might make them relatively insensitive to microtopographic differences in soil moisture, which would explain why there were not many significant interactions between drought and microtopography.
However, the relationship between curvature and growth suggests that moisture plays at least some role in driving differences in performance across microtopography. We found that growth increased with increasing curvature, indicating an advantage for trees in valleys (high curvature areas) during normal years. During drought years, however, we found curvature negatively affected growth (Figure 3), despite our assumption that trees growing in areas with more positive curvature (i.e., valleys) would experience less severe moisture deficits, and thus show less sensitivity to drought. This result could be driven by species sorting along microtopography or individual acclimation as discussed above: if trees growing in valleys are more frequently of drought-sensitive species or are less acclimated to dry conditions, they might show a strong response to drought despite experiencing less severe moisture stress. Further research on how environmental conditions, including soil moisture, soil depth, and soil nutrient status, and intraspecific trait values, vary with microtopography could help tease apart these competing explanations for the relationship between microtopography, drought, and performance.
We found that crowding reduced growth but increased survival (though this effect was only marginally significant). While many studies have shown that competitive effects from crowding have a negative impact on tree performance (Uriarte et al., 2012), others have found no effect on growth and/or survival (Lasky et al., 2014;Uriarte et al., 2016a) or even positive effects (Hurst et al., 2011). Furthermore, sensitivity to crowding varies across species (Uriarte et al., 2004a,b;Canham et al., 2009). Variation in the effects of crowding across studies and our finding that crowding increased survival could have to do with the difficulty of disentangling competitive impacts from the environmental drivers of crowding. Variation in number and size of stems within and across forest stands can be driven by variation in site favorability (Clark and Clark, 2000;Alves et al., 2010;Hernández-Stefanoni et al., 2011). If this is the case, our finding of somewhat enhanced survival in more crowded neighborhoods could reflect site quality. High mortality rates in gaps relative to more crowded, closed canopy areas could also drive the association between survival and crowding. Trees in gaps tend to experience higher rates of herbivory (Richards and Coley, 2007) and may be more subject to desiccation (Denslow, 1980). These hypotheses are also consistent with our finding that crowding effects were not amplified during drought years. However, we found very weak correlations between topographic variables and crowding (Supplementary Table S2). If crowding FIGURE 5 | Relationship between SLA and species-specific estimates of curvature effect on growth. The effect of SLA on the curvature-drought interaction term means that the extent to which the effect of curvature differs in drought years depends on species' SLA. Green points show curvature effect on performance for average rainfall conditions during the non-drought years (2013,2014,2016), beige shows effects during the drought year.
is associated with site favorability, then soil nutrients or canopy gaps, but not topography, underlie variation in site favorability and crowding effects.

Interactions Between Functional Traits, Drought, Topography, and Tree Performance
Plant functional traits reflect tradeoffs in life-history strategies, which are relevant to understanding responses to drought and other disturbances and influence how trees are spatially structured within and across communities (Wright et al., 2004;Reich, 2014). SLA and wood density have frequently been used as proxies for a conservative-acquisitive tradeoff, with high wood density and low SLA indicating slower, more conservative strategies, and low wood density and high SLA representing faster, acquisitive strategies (Chave et al., 2009;Poorter et al., 2010). Our results only partially reflect this tradeoff: wood density was negatively associated with average diameter growth but had no significant relationship with average survival, while SLA was negatively associated with average survival, but had no significant relationship with growth. Because both traits integrate multiple aspects of plant physiology and function, the strength and direction of their associations with growth and survival rates have differed across studies (Yang et al., 2018). For example, while some studies have found the expected growth-survival tradeoff with SLA (Lasky et al., 2015), others have found the opposite, or no relationship between traits and performance (e.g., Poorter et al., 2008). Our results suggest that in this forest, variation in SLA and wood density do encompass at least some of the variation in average performance across species, and can be useful for understanding tradeoffs between "slow" and "fast" strategies.
Differences in life-history strategies are also reflected in inter-specific variation in performance across topography. For example, reduced survival of high SLA species on slopes and enhanced growth at high curvature suggest that high SLA is associated with low stress tolerance, but fast growth in highresource areas. Wood density was not related to variation in performance across topography for either topographic variable in either model. Many studies have linked variation in wood density to shade tolerance (Lawton, 1984;Poorter et al., 2010), and so it is possible that light, which is potentially decoupled from microtopography, is a stronger filter on wood density. Yet overall, the functional traits we considered explained little of the interspecific variation in growth and survival responses to topography. This suggests a need for more mechanistic functional traits, or a need to consider multiple traits and their relationships as part of an integrated whole-plant strategy (Yang et al., 2018).
By incorporating these easy-to-measure functional traits, we can begin to identify mechanisms by which drought affects species differentially. For example, growth of trees with high wood density was less sensitive to drought. During normal conditions, trees with high wood density typically have slower growth due to lower transpiration rates and the higher carbon cost of building denser wood (Chave et al., 2009). However, these trees can also maintain hydraulic conductivity under drier conditions (Chave et al., 2009;Reich, 2014) and so may be able to maintain photosynthetic and growth rates closer to normal during drought. Trees with high SLA tend to have higher photosynthetic rates and accordingly, higher transpiration rates, which predisposes them to more severe water loss during drought (depending on stomatal regulation, Poorter et al., 2009), so we expected that species with high SLA would suffer more severe drought effects. However, we found no significant relationships between SLA and drought response in either model. Though many other studies, including a recent global meta-analysis (Greenwood et al., 2017) have found that wood density and SLA are associated with sensitivity to drought-induced mortality (Metcalfe et al., 2010;Phillips et al., 2010;Uriarte et al., 2016a), the strength and direction of these relationships have differed across studies (Hoffmann et al., 2011;O'Brien et al., 2017).
These discrepancies likely occur because associations between these traits and drought vulnerability are indirect, and driven by to their links with other physiological traits. For example, SLA is often correlated with hydraulic traits including turgor loss point and water potential at 50% loss of conductivity (P 50 ; Mitchell et al., 2008;Xu et al., 2016; but not always; see Markesteijn et al., 2011). The relationships between WD and vessel, xylem, or other wood anatomy traits likely explain resistance to cavitation in trees with higher wood density, though these relationships are poorly understood (Phillips et al., 2010;O'Brien et al., 2017). However, linking these widely and easily measured traits to observed variation in tree species responses to drought and other environmental conditions is an important step toward building a general predictive framework for species' responses to future environmental change. Measuring hydraulic traits such as stomatal conductance, stomatal response to leaf water potential, turgor loss point, water potential at 50% loss of conductivity, or stem water potential would provide more power for predicting species' drought responses (Brodribb et al., 2003;Skelton et al., 2015;Powers et al., 2020). Though some of these measurements can be expensive and time consuming (O'Brien et al., 2017), new methods are making them easier (Maréchaux et al., 2015), providing a promising way forward to improving predictions of drought response. The modeling framework we demonstrate here, which incorporates functional traits, drought response, and environmental variability could easily be applied with hydraulic trait data, and could greatly improve our ability to predict drought responses.
This study provides insight into interactions between climate, topography, crowding, and traits that underlie individual-and species-level variation in drought response, but also highlights the complexity of predicting drought effects and some major research gaps that remain. For example, while drought influenced tree performance via effects on growth, we did not detect a direct effect of drought on survival. However, the strong effect of antecedent growth on survival suggests possible future lagged effects of drought, and the potential for repeated or prolonged droughts to have a major impact. Future studies could manipulate drought conditions (as in the throughfall exclusion experiments in the Amazon) to quantify lag effects and explore effects of repeated drought on tree performance in Caribbean forests. The functional traits we considered helped explain some of the variation in species' responses to environmental conditions, including drought, but their explanatory power was limited, suggesting a role for hydraulic and other more explicitly physiologically based functional traits. Finally, how the effects of microtopography on individual tree performance within plots scale up and relate to the effects of coarse-scale topographic relief across forest stands and landscapes is not well understood, and is an important avenue for future research.

DATA AVAILABILITY STATEMENT
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://luq.lter. network/data/luqmetadata190.

AUTHOR CONTRIBUTIONS
NBS and MU conceived and designed the study, and wrote the first draft of the manuscript. XF analyzed the rainfall data. MU and JZ collected the forest census data. RM, MNU, and NGS collected the functional trait data. NBS performed the statistical analysis. All authors contributed to the manuscript revision, read, and approved the manuscript.