Microclimatic Warming Leads to a Decrease in Species and Growth Form Diversity: Insights From a Tropical Alpine Grassland

Due to warming, changes in microclimatic temperatures have shifted plant community structure and dynamics in tundra and alpine regions. The directionality and magnitude of these changes are less known for tropical alpine ecosystems. To understand the likely trajectory of these shifts in the Andes, we conducted a warming experiment in the northern Andes—using open-top chambers (OTC). In this study, we ask (1) how do OTCs affect air and soil temperatures in microclimates of tropical alpine regions, year-round and during the dry season? (2) What are the effects of 7 years of warming on (a) the aboveground biomass (AGB) and (b) the plant taxonomic and growth form diversity? We installed five monitoring blocks in 2012 at ca. 4,200 m asl with 20 OTCs and 50 control plots randomly distributed within each block. We measured AGB, plant community diversity, and growth form diversity between 2014 and 2019. After 7 years of warming, we found significant increases in mean monthly (+0.24°C), daily (+0.16°C), and night air temperatures (+0.33°C) inside the OTCs, and the OTCs intensified microclimatic conditions during the dry season. Additionally, OTCs attenuated extreme temperatures—particularly in the soil—and the number of freezing events. AGB significantly increased in OTCs, and by 2019, it was 27% higher in OTCs than in control. These changes were driven mainly by a progressive increment of tussock grasses such as Calamagrostis intermedia, typical of lower elevations. The increase of tussocks led to a significant decrease in species diversity and evenness inside OTCs, but not in species richness after accounting by sampling time. Furthermore, cushions and herbs decreased inside OTCs. Our results show that experimental warming using OTCs in equatorial regions leads to decreased daily thermal amplitude and night temperatures rather than the level of increase in mean temperatures observed in temperate regions. The increase of tussocks and decrease in diversity of species and growth forms due to prolonged modifications in microclimatic temperature might be a step toward shrub-dominated ecosystems. Further research on this topic would help understand shifts in growth form dominance and the direction and rate of change of the system.


INTRODUCTION
In tropical alpine ecosystems, macroclimatic conditions and large-scale climatological phenomena play an important but secondary role in limiting plant life compared to microclimatic conditions (Sklenář et al., 2015). Microclimates are subtly affected by the topography, substrate, and aspect of the landscape and are also modified by plants (Anthelme et al., 2014;Opedal et al., 2015). While regional climate warming predictions for tropical alpine regions forecast an increase in temperature between 1 and 5 • C with mean values of 2 and 3.5 • C (Vuille et al., 2018), at a small scale, these projected increases in temperatures may differ . In tropical alpine ecosystems, the interaction between microclimatic change and its effect on vegetation is not well understood. One potential strategy is that plants might alter important functional characteristics to cope with this change (Körner, 2003) or become less dominant or even locally extinct. Therefore, it is vital to examine how changes in microclimatic temperatures affect plant diversity and productivity to further assess the consequences of climate warming on tropical alpine ecosystems.
Temperature, humidity and seasonality gradients due to elevation and latitude (Weberbauer, 1945;Troll, 1959;Sarmiento, 1986;Cuesta et al., 2017Cuesta et al., , 2020 define different ecological niches in tropical alpine systems. To fill this mosaic of niches, tropical alpine plants developed various morphological and physiological adaptations (e.g., Cleef, 1978;Smith and Young, 1987;Ramsay, 2001;Ramsay and Oxley, 2001). For example, plant community composition and growth forms change gradually as elevation increases as a coping response to harsh environmental conditions particularly, cold temperatures, where facilitation mechanisms prevail over competition (Sklenář and Ramsay, 2001;Anthelme et al., 2012).
To empirically examine the effects of climate change on mountain ecosystems, experimental studies shorten the space-time scale (i.e., warming experiments and transplant experiments) and link global warming with the dynamic changes of plant communities and their effect on ecosystem functions (Tito et al., 2020). Among the benefits that space for time experimentation provides is the potential to compare with studies that look at changes over a prolonged timeframe. Passive warming experiments use various methodologies and structures to simulate the temperatures and conditions predicted for the future. While similar studies have used open-side chambers (Pancotto et al., 2020), International Tundra Experiment (ITEX) corners, and circular chambers (Rustad et al., 2001), the method that has been successful in comparable systems is the open-top chamber (OTC) designed by ITEX (Marion et al., 1997). This OTC has a hexagonal shape made of polycarbonate sides hinged together to have a smaller area at the top of the structure than at the bottom. It allows airflow above, below and between the corners of the chamber yet retains warmth when shortwave radiation enters and delays the escape of longwave radiation, as would occur in a greenhouse, passively and minimally increasing the area's temperature.
Currently, passive warming studies have been reported mainly for the tundra and Northern hemisphere alpine ecosystems and two in southern South America (Rustad et al., 2001;Sierra-Almeida and Cavieres, 2010;Pancotto et al., 2020;Fazlioglu and Wan, 2021). These studies have focused on understanding vegetation community responses to climate warming, including changes in biomass, plant phenology, freezing resistance mechanisms, species distribution shifts, increased dominance of certain species, and changes in the alpha and beta diversity of the ecosystem (Klein et al., 2004;Elmendorf et al., 2012;Liu et al., 2018;Yang et al., 2018). Studies that have used OTCs have reported temperature increases in mean air temperatures of 1.7 • C for tropical alpine ecosystems (Lasso et al., 2021) and increases in daily soil temperatures of up to 0.4 • C (Klein et al., 2004(Klein et al., , 2007Pancotto et al., 2020), as well as changes in the number of freezing events and their intensities (Sierra-Almeida and Cavieres, 2010). Presently, only one study has reported the effect of experimental warming on plant communities of tropical alpine ecosystems (Lasso et al., 2021), but no study has reported the long-term (over three years) effects of experimental warming in these systems.
The observed impacts of passive warming experiments on tundra/alpine vegetation are shifts in the abundance of dominant species, dominant growth forms, and changes in aboveground biomass (AGB). Studies that successfully undertook experimental warming (Klein et al., 2004;Yang et al., 2018) attributed the change in species composition to an increase in competition between current dominant species and species that migrated from lower elevations, which resulted in a decrease in species richness. For example, Klein et al. (2004) reported that warming decreased species richness by five to 14 species, depending on habitat and grazing history, in one year and by nine to 15 species in another year. Regarding growth forms, Lasso et al. (2021) reported shifts in growth form dominance, attributed to an increase in graminoids and shrubs, observed after three years of warming in the tropical alpine region of Colombia. Therefore, although it seemed that the upward migration of new species would increase species richness (Steinbauer et al., 2018(Steinbauer et al., , 2020, these experimental studies suggest that increase may be temporary. Moreover, the effect of climate warming could be the loss of species and growth forms characteristic of tundra/alpine ecosystems. Regarding AGB, studies have reported decreases (by a maximum of 23%; Fu et al., 2013), increases (by 72 g/m 2 in 2000 and 102 g/m 2 in 2001 in one site; Klein et al., 2007), or no effect on AGB or aboveground productivity due to experimental warming (Liu et al., 2018). Given these results, it is important to assess the direction and magnitude of the effect of climate warming on tropical alpine vegetation.
In this study, we analyze for the first time the effects of a prolonged passive warming experiment (7 years) on plant communities in the páramo habitats of Ecuador. By implementing a warming experiment using OTCs installed in the Yanacocha Reserve in 2012, we seek to answer the following questions: (1) How do OTCs affect air and soil temperatures in the microclimate of tropical alpine regions year-round and during the dry season? Furthermore, how are these temperatures affected in the dry season? (2) What are the effects of 7 years of continuous passive experimental warming on (a) the AGB and (b) species diversity, and (c) growth form diversity of tropical alpine ecosystems?
In terms of the effects of the OTCs on microclimate temperatures, given the results reported from multiple studies (Klein et al., 2004(Klein et al., , 2007Sierra-Almeida and Cavieres, 2010), we expected daily mean air temperature to increase a minimum of 0.4 • C and a smaller increase in soil temperature to a maximum of 0.4 • C. We expected these differences to be particularly noticeable in the dry season. Regarding biodiversity, we expected to find a decrease in species diversity and evenness due to a dominance of certain species characteristic of lower elevations. Given the gradual change in species and growth form diversity as elevation increases, we expected changes in growth form diversity. For example, we expected an increase in graminoid and shrub species because their growth form might be more resilient to changes in microclimatic conditions or, in the case of shrubs, higher temperatures result in favorable conditions for them (Harte and Shaw, 1995;Cruz-Maldonado et al., 2021;Lasso et al., 2021). While there has not been a clear trend for AGB under warming conditions, we expected an increase in AGB due to the previously mentioned increase in shrubs and graminoids that may contribute more biomass per square meter than other growth forms.

Study Area
The Yanacocha Reserve, located about 15 km northwest of Quito, the capital of Ecuador (Figure 1A), was created by the Jocotoco Foundation, a World Land Trust Partner spanning 1,200 ha (960 ha of montane evergreen forest and 240 ha of páramo). The study area, embedded in the eastern flank of the Pichincha Volcanic Complex, is located in the western mountain range of the Andes in Ecuador. The western mountain range is found in Ecuador's north and central regions and is geologically more recent than the eastern mountain range (Hughes and Pilatasig, 2002). It is influenced climatically by air masses that originate in the Pacific Ocean (Vuille et al., 2000). Precipitation varies seasonally-there is more precipitation between January and May and October to December, whereas June through September is the dry season (Navarro-Serrano et al., 2020). The average annual precipitation is 1,500 mm year −1 (Rugel, 2019) and less than 250 mm for the dry season (Sklenář et al., 2008). Therefore, a decrease in monthly rainfall during the dry season, coupled with increased UV radiation due to cloud-free skies leads to augmented thermal convection, resulting in higher maximum and minimum daily temperatures. Soils within the Yanacocha Reserve study area are Andosols, typical of the high Andes (WRB, 2015) with well-defined layers of volcanic ash intercalated with the organic layers. The topsoil layer in both forest and páramo varies from 15 to 18 cm in depth containing more than 7% of soil organic matter (Calderón-Loor et al., 2020).

Experimental Design
In 2012, the Yanacocha Reserve became the location of a longterm ecological research site to monitor vegetation changes at different elevations and under different scenarios (Teunissen van Manen et al., 2019Calderón-Loor et al., 2020;Llerena-Zambrano et al., 2021). For example, the Yanacocha Reserve houses two of the GLORIA monitoring summits for the Pichincha Volcanic Complex Tovar et al., 2020). In July of 2012, a passive warming experiment was established. It started only with control plots using a hierarchical block design and later, after the installation of the OTCs in 2012 (Figure 1B), the monitoring of the experimental plots started in 2013.
Five monitoring blocks-installed at 4,200 m asl (−0.13518 Lat, −78.57458 Long)-contained ten units of 6 m × 6 m  Figure 1C). Furthermore, OTCs were placed randomly in four units per block, covering the biodiversity and AGB experimental plots. Therefore, four units per block had an OTC sub-unit, making a total of 20 OTC sub-units and 50 control sub-units among all five blocks . The blocks were placed at a minimum distance of 5 m from each other to reduce spatial autocorrelation. Through incidence-based species accumulation and rarefaction curves-using the functions specaccum and rarecurve both from the Vegan package in R (Oksanen et al., 2020)-the saturation of the curves showed an appropriate sampling effort (Supplementary Figure 1).
To confirm that temperature caused the effects observed in the OTCs and that any change in biomass or diversity would not be attributed to a possible reduction in herbivory-caused by the physical blocking of the OTC walls (Marion et al., 1997)six herbivore exclusion units were installed among the five blocks in 2012. The exclusions consisted of 1-meter-tall fences buried 30 cm below ground and covered the area of a whole unit within a block. Before installing the exclusion, we cleared all small mammals: genera Thomasomys, Cryptotis, Sylvilagus, and Microryzomys. To prepare the area, we trapped and then released them outside of the exclusions.

Field Methods
Control and OTC biodiversity was monitored once a year between 2014 and 2019 (Table 1) using permanent the 1 m × 1 m plots, sub-divided in 10 cm × 10 cm squares to determine species cover (percent) and composition. In addition, we collected control AGB samples from one of the two designated AGB plots located within one of the sub-units, following Calderón-Loor et al. (2020), once per year between 2012 and 2019 (Table 1 and Figure 1C)-through two strips of 10 cm × 1 m at the top and the bottom of the plot. In the OTC sub-units, we collected AGB from the outer limit of the biodiversity plot, where 1 m × 10 cm strips were harvested outward instead of inward, as was the case for the control plots ( Figure 1C). AGB samples were tagged, bagged, and placed inside plastic coolers to preserve the integrity of the collected samples before laboratory processing. In sum, the monitoring included assessing the effect of temperature on the following variables over time: AGB, which refers explicitly to live plus dead plant material collected above ground level, and vegetation composition.
Thermochron ibuttons R (DS 1922L-F5#) were installed in the soil (at 10 cm in depth) and at 10 cm above the ground to measure soil and air temperature, respectively, in 18 plots (9 for both control and OTCs), making a total of 36 loggers, within blocks 1 through 3. The loggers recorded the temperature every hour. We also placed Hygrochron R temperature and humidity data loggers (DS 1923-F5#) when installing the blocks, but they failed after the second year of sampling. Therefore, for this study, we focused on temperature data only.

Laboratory and Herbarium Methods
We first weighed the AGB samples in the field with a manual scale to obtain fresh weight. Then we transferred the samples to a laboratory where we sorted them between above ground necromass and live mass. We distinguished between necromass and live mass based on the color and rigidity of the plantsyellow/brown leaves for necromass and rigid green leaves for live mass (Calderón et al., 2013). Samples were then dried at 60 • C for 72 h in an oven until they reached a stable weight. Once the samples were dried, we weighed them again separately. This study's statistical analysis used the sum of the dry aboveground necromass and live mass. At the QCA Herbarium, we processed and deposited vouchers of plants identified in the field and confirmed the plant identifications (Supplementary Table 1). Identification of plant vouchers were done in the QCA Herbarium where samples were also deposited. We identified 51 out of 70 morphospecies to species level. The rest of the 19 morphospecies were left at their most substantiated level and grouped into family groups. For example, the "Apiaceae group" housed individuals that look like species within the Oreomyrrhis and Niphogeton genera but were not identified up to species. The "Asteraceae group" included individuals that look like species from the Aphanactis, Gamochaeta, and Achyrocline genera. The "Poaceae group" included most species from the Poaceae family except C. intermedia and Calamagrostis fibrovaginata Laegaard, grouped in their own Calamagrostis group (Supplementary Table 1).

Trait Data
We obtained the growth form data for each species and morphospecies from the GLORIA-Andes database (Muriel et al., unpublished data), which uses a revised version of the Ramsay and Oxley classification (Ramsay and Oxley, 1997). We used the following growth form categories: cushion, herb, rosette, shrub, and tussocks.

Temperature Analysis
Temperature data was downloaded from all loggers twice a year from June 2012 to 2019 to obtain a complete yearly log of hourly temperatures. However, due to a malfunctioning in some of the loggers, gaps in the annual sequence occurred in 2015, 2016, and 2018. Therefore, we bootstrapped the dataset, with the boot function in the boot package in R, to resample the years 2012, 2013, 2014, 2017, and 2018-using 1,000 bootstrap replicates-to minimize the impact of gaps in the dataset (Canty and Ripley, 2021). With this resampled data, we calculated summary statistics-including freezing events and intensities and the mean, minimum and maximum temperatures for air and soil, for every day of the year, of the dry season and night temperatures for OTC and control loggers-as well as the 95% confidence intervals (Crowley, 1992;Dixon, 1993).
We estimated statistics for the dry season using temperature data from June to September. Night temperatures were estimated using the values between 6 PM and 6 AM. Finally, we calculated freezing times by counting the number of times in 24 h the temperature was below zero and the intensity of the freezing period by counting the consecutive hours the temperature was below zero in 24 h (Sierra-Almeida and Cavieres, 2010).
We assessed the difference between control and OTC plots for every month of the year with T-tests on the temperature data (not bootstrapped) and Wilcoxon two-sample test for freezing lengths (Sokal and Rohlf, 2012)-using the stat v3.3 package in base R (R Core Team, 2021). Afterward, we used a linear mixed-effects model to assess the effect of the OTC on temperature, where the treatment interacts with the position of the logger (either air or soil), and the random effects are the hour of the day and the logger number. We used the lme4 package in R (Bates et al., 2014). We ran this model for the total temperature dataset, dry season temperatures, and night temperatures separately, using the bootstrapped data.

Assessing the Effect of Warming
Aboveground biomass was estimated as the sum of dry aboveground necromass and live mass in kilograms per square meter and then tested for normality using Q-Q plots in R. Since we did not obtain normality from these data, the dataset was log-transformed for our analysis. No transformations were done for any of the other variables since they all met the normality assumption.
Like with temperature, we compared all diversity metrics and AGB between control and OTC plots through two methods. First, to show differences when we do not account for time and other factors, we used the Wilcoxon two-sample test to compare warmed and control plots. Then, we used generalized linear mixed-effects models and general linear mixed-effects models to analyze the effects of the warming treatment and time (year) on these response variables using the glmer and lmer functions from the lme4 package in R (Bates et al., 2014). For AGB we tested the following three models: (a) biomass from control plots against OTC plots, (b) biomass from exclusion plots against control plots, and (c) biomass from exclusion plots against OTC plots. Finally, we used a Poisson distribution for the species richness model.

Microclimate
Compared to controls, daily mean air temperature increased in OTCs by 0.16 • C ( Table 2). During the dry season, the OTCs had the most increase in air temperature for the daily average maximum by 2.33 • C ( Table 2). At night, the mean temperature increased for OTCs compared to controls by 0.33 • C, and the highest increase was in the average maximum temperature by 1.31 • C ( Table 2). On the contrary, the OTCs significantly decreased the daily mean soil temperature by −0.08 • C ( Table 2). The highest increase in soil temperature was for the daily average minimum by 0.32 • C ( Table 2).
The OTCs significantly increased the mean monthly air temperature by an average of 0.24 • C except in September and October, at the end of the dry season (Figure 2A). In contrast, the monthly mean and maximum soil temperature in OTCs were lower than controls from July to October and April to December, respectively ( Figure 2B). The OTCs attenuated soil monthly minimums and maximums in most of the year except for the height of the wet season in February and March ( Figure 2B).
There were significantly fewer freezing events and of less intensity in OTCs than in control in soil and air during the dry season, but OTCs and control had a similar number of freezing events in the air during the whole year (Supplementary  Tables 3, 4). Temperature loggers in the soil recorded zero freezing events during the dry season (daily and nightly; Supplementary Table 3).
The biggest difference in AGB within the same year between OTC and control was in 2019 with 27.9% (p = 0.007). However, AGB decreased over time in both OTCs and control plots as shown by the negative values of the mixed-effects model coefficients ( Figure 3B and Supplementary Table 5). Exclusion plots (mean = 1.24 kg/m 2 ; sd = ± 0.16) did not show significant differences with OTC plots (mean = 1.11 kg/m 2 ; sd = ± 0.16) or control plots (mean = 0.93 kg/m 2 ; sd = ± 0.16; Figure 3B and Supplementary Table 6).
In the mixed-effects model the years 2016, 2017, and 2019 had lower diversity values than 2015 ( Figure 4B and Supplementary  Tables 7, 8).
Evenness in the OTCs (mean = 0.46; sd = ± 0.17) was significantly lower than in control plots (mean = 0.54; sd = ± 0.16) over the entire study period (W = 3,434; p = 0.0088; Figure 4A). This pattern was further confirmed when time was held as a fixed factor in the mixed model, and evenness for 2016, 2017, and 2019 was lower than in 2015 (Supplementary Table 8 and Figure 4B).
There was no significant difference in species richness between control plots and OTCs ( Figure 4A). Only in 2019, species richness was significantly lower in OTCs (mean = 9.6; sd = ± 2.4) than in control (mean = 12.19; sd = ± 2.7) as was suggested by the Wilcoxon test (p = 0.015) and visualized by the mixed-effects model (Supplementary Table 7 and Figure 4B). Nevertheless, by 2019, five out of the 57 recorded species [e.g., Astragalus geminiflorus Bonpl., Azorella pedunculata (Spreng.) Mathias & Constance] were no longer present inside the OTCs, and six more demonstrated a progressive reduction in the warmed plots (Supplementary Table 1). On the other hand, only two species appeared in the OTC plots as the warming experiment progressed without being present in the control plots, Bartsia laticrenata Benth. and G. acostae Cuatrec.

Growth Forms
Growth form diversity was significantly lower (p = 0.0026) in the OTC plots (mean = 2.36; sd = ± 0.7) than in control plots (mean = 2.73; sd = ± 0.71) but decreased in both treatments over time (Figure 5A). This difference holds after accounting by sampling time in the mixed-effects model (Figure 5B and Supplementary Table 9). Growth form evenness was also significantly lower (p = 0.003) in OTC plots (mean = 0.56; sd = ± 0.2) than in control plots (mean = 0.66; sd = ± 0.17) over the entire study period and decreased over time in both ( Figure 5B). However, after accounting by time in the mixedeffects model, evenness in OTCs was not significantly different from control, and all times had lower evenness compared to 2014 ( Figure 5B and Supplementary Table 9).
Percent cover of cushion growth form was significantly lower in OTCs than in control plots, while tussock cover was significantly higher over the entire study period ( Figure 6A and Supplementary Table 10). Over time, cushion species cover decreased in the OTCs but increased in control plots ( Figure 6B and Supplementary Table 10). Meanwhile, percent cover of tussock species (e.g., Calamagrostis group) increased over time in OTCs by 16.5% but kept similar cover values in control plots ( Figure 6B and Supplementary Table 10).

The Role of Temperature in Microclimatic Conditions
In our study, the OTCs changed the microthermal conditions, potentially leading to shifts in AGB, plant community composition, and growth form diversity. Generally, air temperature and soil temperatures had different responses to the OTCs, where air temperature on average increased and soil temperature on average decreased. Specifically, the thermal conditions changed by first, a general increase in daily and monthly air temperature with the highest increase in means occurring for night temperatures. Second, the dry season accentuated the differences in microclimatic conditions between control and OTC plots, for both air and soil temperatures, where soil temperature in OTCs was mainly lower than control, and air temperature in OTCs had its highest increase in temperature compared to controls. Third, extreme temperatures were attenuated, including reducing freezing events. Daily (+0.16 • C) and monthly (+0.24 • C) mean air temperature increased but in a lower amount than what we hypothesized and less than observed in tundra and subtropical grassland environments. For subtropical grasslands, Buhrmann et al. (2016) reported an increase of 1.6 to 2.4 • C in air temperature during the day in OTCs compared to controls and 0.0 to 0.6 • C during the night across all four seasons. In our study, night temperatures in OTCs were significantly different from control and showed higher mean air temperatures by an average of +0.33 • C (p < 0.001), which falls within the range of change in temperate regions. Our results show that in equatorial regions, larger changes occur for night temperatures than from daily temperatures, different from temperate and subtropical alpine areas. Additionally, lower values of temperature change than Lasso et al. (2021) in the Colombian páramos may result from differences in OTC design, such as the type of fiberglass used in their warming experiment.
Our study found smaller differences in the daily mean soil temperature between OTCs and control than in the mean air temperature. This result was similar to the results reported by Pancotto et al. (2020) in southern Patagonia. However, a metaanalysis of 32 experimental sites across four temperate biomes: tundra, low tundra, grassland, and forest found a broader range, of 0.3 to 6.0 • C, for mean soil warming temperature increases (Rustad et al., 2001). The studies analyzed in that meta-analysis used different warming methods, not only OTCs, which, along with latitudinal differences, may contribute to the differences from our results. For example, in subtropical grasslands, the absolute maximum temperature in OTCs increased by 0.9 to 1.8 • C while the absolute minimum was increased by 0.3 to 0.8 • C (Buhrmann et al., 2016). Meanwhile, we obtained a mean maximum temperature increase of 0.55 • C and an increase in mean minimum temperature of 0.4 • C.
The OTCs stabilized soil temperatures by keeping the monthly minimums and maximums consistent throughout the year with maximum temperatures lower than those in control plots. This reduction in the soil thermal range may reduce plants' ability to tolerate extreme events under normal conditions (Sierra-Almeida and Cavieres, 2010). Further, soil temperature attenuation also reduces the frequency and intensity of freezing events (Supplementary Table 3). In this study, the mean monthly soil minimum temperature for OTCs was above freezing. This attenuation could be caused by greater evapotranspiration during the day, which could generate a buffer that reduces heat transfer through convection during the night (Sierra-Almeida and Cavieres, 2010). Consequently, the thermal amplitude experienced by plants is reduced, which would alter their physiological response to extreme conditions (Sklenář et al., 2016).
The dry season analyses revealed an important pattern in the temperature trends observed in the OTCs, including the intensification of temperature and shifts in trends around the start and end of the dry season. In our results, at the beginning of the dry season (June), air temperatures in OTCs intensified, obtaining their maximum difference from controls in July by 0.7 • C. During June and July, the OTC trapped higher amounts of UV radiation, causing this peak in temperature (Marion and Pidgeon, 1992). Additionally, monthly maximum air temperature in OTCs only fell below the control plots in January and February, during the wet season, where days have higher humidity, cloud cover, and lower radiation (Vuille et al., 2003;Letts and Mulligan, 2005). Studies in tundra found that passive warming treatments increased mean growing season air temperature by 1 to 3 • C, a warming range predicted and observed for tundra regions (Walker et al., 2006). Therefore, variations in monthly precipitation and UV radiation over a year-long period are important in describing OTC temperature trends.

Relationship Between Biomass, Biodiversity, and Growth Forms
Over the study period, OTCs showed higher values of AGB but lower values of species and growth form diversity and evenness than control plots. This pattern seems to be linked to an increase in the vegetation cover of tussocks and a decrease of cushion and herb species in the OTCs, leading to a less diverse community. While AGB reduced in both control and OTCs over the study period, AGB was highest in exclusion plots followed by OTCs, and the lowest means of AGB were found in controls.
By the end of our 7-year experiment, OTCs had 27% more AGB than control plots, as we hypothesized, and which is consistent with other experiments in the alpine meadows of the Hedguan mountains (Yang et al., 2018) and Patagonia (Pancotto et al., 2020). However, other studies have reported no warming impact or decreases in AGB. For example, Fu et al. (2013) found that AGB significantly decreased with experimental warming at 4,313 m asl. This study used different methods than ours and had a much shorter experimental time of 1 year, potentially explaining the difference from our results. Although OTCs generally had more AGB than controls, we observed an overall decrease of AGB for both controls and OTCs. This contrasts with results reported by Calderón-Loor et al. (2020) in the same páramos of the Yanacocha Reserve for the period 2012-2014 in which AGB increased, despite showing similar values of AGB to ours for 2014. Finally, the higher amount of AGB in exclusion plots than OTCs shows that the OTCs did not entirely exclude herbivory by small mammals, as was also observed by Barrio et al. (2021). However, herbivory intensity is dependent on site and OTC design. In our study, given the visual evidence of herbivory in OTCs, more AGB in OTCs may be more likely related to the changes in thermal conditions rather than the potential exclusion effect of herbivores by the OTCs. However, further research on herbivory intensity related to OTCs in tropical alpine regions would help clarify this topic.
The higher AGB in OTCs (by 9%) in our study area is mainly related to increased vegetation cover of tussocks (e.g., C. intermedia), which increased by 16.5% by the end of the experiment and it is a pattern that has been observed in other studies as well. For example, in subtropical alpine grasslands, annual graminoid biomass production significantly increased by nearly 20% in OTCs (Buhrmann et al., 2016). Similarly, Na et al. (2011) found an increase in AGB inside OTCs after 2 years of warming in the Tibetan Plateau, where sedges were responsible for 90% of the total biomass gain. Lastly, in the páramos of Colombia, graminoid and shrub coverage increased inside OTCs in one out of the two study sites (Lasso et al., 2021).
In our study, an increase in tussocks in OTCs could be due to three reasons: (1) tussocks have specific adaptations that allow them to become dominant, (2) the drier environments created by OTCs favor tussocks, and (3) at this point, tussocks are favored but may be replaced later by shrubs and woody plants. First, tussock grasses have broad thermal niches  and are highly tolerant to freezing (Márquez et al., 2006;Rada et al., 2019). Second, they compete by being taller, shading out shade-intolerant species, which leads to local/plot extinctions (Wong, 1991), dispersing by wind , and by being photosynthetically efficient (i.e., reduced leaf area index and higher leaf dry matter content (Leegood, 2002;Cuesta, unpublished data). Next, a secondary effect of the OTCs is that they create drier environments (e.g., Lasso et al., 2021), and tussocks are resistant to cold and resistant to drought (Monteiro et al., 2011;Rada et al., 2019). Therefore, tussock grasses seem suited to take advantage of conditions created by increased temperatures coupled with increased water limitation.
Lastly, we hypothesized that both graminoid and shrub coverage would increase but only the increase in graminoid cover was significant. In their monitoring studies, Porro et al. (2019) and Hamid et al. (2020) found that the vegetation cover of shrubs and graminoids increased with warmer temperatures. However, in our study, we found that shrubs as a group did not increase, but that individual species, such as C. jussieui, showed an increase in cover. The magnitude of this trend was small and may depend on more prolonged warming exposure periods and higher warming temperatures. Previous studies on the tropical upper forest line suggest that woody plants cannot grow adequately at soil temperatures below 5 or 6 • C (Hoch and Körner, 2003;Körner and Paulsen, 2004;Hoch and Korner, 2005;Wiley and Helliker, 2012). Therefore, in the long term, as the world continues to get warmer in the future, shrubs [e.g., Hypericum laricifolium Juss., Hesperomeles ferruginea (Pers.) Benth.] and trees (e.g., Polylepis pauta Hieron., Gynoxys acostae) might outcompete tussock grasses, and there is evidence for this in the paleo records (Flantua et al., 2019). Here our hypothesis was partially met, so given our results, tussock dominance could be a temporary state on the way to a forested landscape.
The increase in percent coverage of tussock species and some shrubs (e.g., C. intermedia, C. jussieui) in OTCs is related to the downward trend observed for diversity and evenness of species and growth form, as we hypothesized. This downward trend is consistent with observations made in a 14-year monitoring study in the Apennines (Porro et al., 2019). Furthermore, Walker et al. (2006) found that warming increased the cover of deciduous shrubs and graminoids and decreased species diversity and evenness. While our results for species evenness and diversity and growth form diversity were significant, species richness was not, yet the downward trend for species richness was similar to that found in previous studies in alpine/tundra regions (Klein et al., 2004;Yang et al., 2018). For example, after the fourth sampling time, five out of the 57 recorded species (e.g., A. geminiflorus, A. pedunculata) were no longer present inside the OTCs, and six more evidenced a progressive reduction in the warmed plots (mainly herbs and cushions). This result was contrary to that of Steinbauer et al. (2018), who found increased species richness in nine regions across Europe over 145 years, mostly related to temperature change. However, after 20 years of monitoring, Steinbauer et al. (2020) showed that previously recorded species richness increases in the Alps were slowing down due to the disappearance of species, specifically coldadapted species. In our study, decreasing species richness was more in line with this trend. Yet, the slower rate of decrease compared to other OTC studies may be due to the OTCs raising the daily mean temperature less than the increased mean temperature in those studies.
Besides the increase in tussock percent cover, we also observed decreased cushion (e.g., Plantago sericea Ruiz & Pav.), herb (e.g., Geranium humboldtii Willd. ex Spreng.), and rosette cover in the OTCs. Cushion plants are stress-tolerant, which allows them to be a dominant growth form at the highest elevations, along with rosettes, which decreased in both control and OTC plots in our study (Sklenář et al., 2010;Anthelme et al., 2017;Rada et al., 2019). However, when subjected to increased temperatures, Pancotto et al. (2020) found that their study subject Astelia pumila (G. Forst.) R. Br., a cushion plant, decreased their photosynthetic efficiency. Thus, the increase in temperature could be a potential explanation for the decrease in cushion plant coverage in our site. The reduction in cushion plant cover, together with the overall decrease of rosettes, which in tropical alpine ecosystems are characteristic of high elevations  (Sklenář and Balslev, 2007;Rada et al., 2019), in OTC and control plots over time, may suggest that a macro-scale process may be already occurring in the Yanacocha Reserve. One potential explanation for this pattern is that rosettes are already migrating upward. On the Antisana volcano, about 63 km southeast of our study site, Moret  Time is an essential variable in warming experiments. OTCs accelerate the rate to visualize change compared to prolonged comparative studies (Porro et al., 2019;Steinbauer et al., 2020;Anthelme et al., 2021), and our study found significant changes after 7 years of sampling. However, studies in non-tropical ecosystems found differences in a much shorter timeframe; for example, studies in tundra found changes in plant communities after two growing seasons (Walker et al., 2006). The distinct seasonality found in tundra and sub-tropical ecosystems, contrasting from the tropics, may make the tropics more resilient to temperature changes and may reveal the effects of climate warming over more time. Furthermore, those studies reported much higher temperature changes than what we obtained in our OTCs, which could also be contributing to the longer response times.

CONCLUSION AND FUTURE RESEARCH
Given the rapid rate of global temperature rise, it is crucial to understand its effect on ecosystems that support millions of people and are home to a diverse and unique group of life forms. In this study, the first long-term OTC study in a tropical alpine ecosystem, we experimentally increased the temperature of vegetation plots to assess the effects of warming on AGB and taxonomic and growth form diversity. While we hypothesized that air temperature would increase like other studies in similar regions, we obtained a smaller increase than expected. However, OTCs did modify the microclimatic conditions that plants were experiencing. Two limitations of this study in characterizing the microclimatic conditions and their relationship to tropical alpine vegetation are the lack of humidity data together with photosynthetically active radiation and the use of one study site. Having multiple experimental sites would allow for more robustness and the feasibility of upscaling our results into a regional context. Shifts in thermal conditions caused by the OTCs altered AGB, taxonomic, and growth form diversity, confirming the importance of small-scale temperature on shaping plant community composition and vegetation physiognomy in alpine ecosystems. Within OTCs, we found a trend toward the dominance of one or a few species and growth forms, specifically tussock grasses, leading to a decrease in species and growth form diversity and evenness, specifically tussock grasses. The dominance of tussocks may be a transitional phase; in the long term, as warming continues to increase, shrubs may outcompete the tussocks.
We obtained two unexpected results from this study: (1) a decrease in total biomass in both control and OTCs, and (2) a decrease in rosette coverage in both OTCs and control plots. Further research is needed to elucidate the overall decrease in AGB in the Yanacocha Reserve under control and experimental conditions. Measuring necromass decomposition rates together with soil properties (moisture, respiration, and N mineralization rates), belowground biomass, and soil organic carbon will help improve our understanding of the mechanisms driving AGB changes over time. Additionally, AGB contributions by growth forms are needed to understand shifts in growth form dominance and their links to resource acquisition traits. Unlike in shrub páramo or Andean forests, demographic patterns in the grass páramos-such as species mortality and recruitment rates-have not been studied yet. Furthermore, because many species have clonal reproduction, it is challenging to determine individuals and perform population dynamics studies in tropical alpine systems. In our study, analyses on population dynamics would further the results reported here by providing information on immigration from lower elevations and possible extinction processes.

DATA AVAILABILITY STATEMENT
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

AUTHOR CONTRIBUTIONS
FC, PM, and SD conceived the experiment. SD, FC, PM, and CT built the general structure of the datasets. SD, FC, and CT analyzed the data and wrote the manuscript. SD, FC, PM, RJ, ES, and EP collected field data. All co-authors commented on and/or approved the manuscript. All authors contributed to the article and approved the submitted version.