Evaluation of the adaptive potential of three non-native multipurpose species for soil rehabilitation

Soil rehabilitation involves restoring their capacity to provide goods and services; however, there is a lack of information regarding the survival and growth responses of individual species to local environmental variations. This study evaluated the adaptation capacity of three non-native multiple-use species in agroforestry design, in a degraded semiarid temperate climate zone located in the Tula watershed in Mexico, assessing their survival, growth of morphological variables, and death risk for use as soil restoration plants. According to the slope of the terrain, the species were planted in an agroforestry arrangement of lines and ridges. Plants were selected according to their multipurpose potential, including their taxonomy, form, function, and use, as well as their soil and climatic conditions. Survival rates for all three species were below acceptable levels (<49%)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$(< 49\%)$$\end{document}. Neither Senna multiglandulosa nor Sedum dendroideum adapted well to the soil and climate of the site, since their survival rates were below acceptable (∼17<49%)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$(\sim 17<49\mathrm{\%})$$\end{document}. Aloe sp. adapted to the edaphoclimatic conditions of the site because its survival rate was excellent (∼97≥90%)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$(\sim 97 \ge 90\%)$$\end{document}. Climate conditions during the autumn and winter affected growth rates of Senna multiglandulosa and Sedum dendroideum. Aloe sp. growth rates were only affected by climatic conditions during the winter months. Senna multiglandulosa's final height reduced death risk by 1.50%\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1.50\%$$\end{document}. Sedum dendroideum's final diameter, initial canopy cover and canopy growth rate reduced the risk of death by 55.1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$55.1$$\end{document}, 5.1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$5.1$$\end{document} and 30.7%\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$30.7\%$$\end{document}, respectively. The height growth rate of aloe and initial canopy cover reduced the risk of death by 93.9 and 2.4%, respectively. Plant species evaluated showed varying levels of adaptation. Aloe sp. was the species that adapted best to the site's soil and climate.


Introduction
Rehabilitating a soil is the process of restoring its capacity to provide goods and services again, but it is important to take into consideration that a rehabilitated soil is not the same as it was before degradation. Intervention strategies aim to restore the productivity of the area rather than to restore its original structure (Peri et al. 2021;Chazdon 2008). Seeding or planting of woody plants, soil retention and conservation tools, or a combination of the two, are the most common practices worldwide (Mata-Balderas et al. 2014). However, each ecosystem has its own dynamics of disturbance and subsequent recovery, depending on the type, intensity, and frequency of disturbance. Thus, there are no universally effective strategies, but rather complementary techniques that must be evaluated in order to determine which are the most effective methods to recover an ecosystem based on its particular properties (Ceccon 2013; Duarte et al. 2018;López-Barrera 2015). Rehabilitation activities are designed to ensure that plants achieve high levels of survival and growth during the field establishment step (Navarro et al. 2006).
Many factors influence the survival of species at an early stage and the success of plantations, such as structural characteristics, morphology, soil properties, nutrient availability, Communicated by Agustin Merino. climate, soil preparation, planting dates, humidity, weed control, site-specific characteristics as well as many other factors (Aguilos et al. 2020;Davis and Jacobs 2005;Navarro et al. 2006;Bernaola-Paucar et al. 2015;Palacios-Romero et al. 2017). Identification of the limiting environmental factors in an area of interest prior to planting is important to determine what morphological and physiological characteristics the vegetative material should possess for fast rooting and adaptation to site conditions (Navarro et al. 2006;Landis et al. 2010). Despite the fact that diameter is recognized as an important attribute in nursery plants in order to ensure survival (Tsakaldimi et al. 2013), as well as tolerance of adverse climate conditions (Prieto et al. 2011), which is due to its relationship with root volume and taproots (Jacobs et al. 2009), morphological characteristics are not the only factors that explain the variation in responses after planting.
The Mexican ecosystem is highly heterogeneous, creating an environment that promotes great biological diversity and an important natural capital, and is among the five nations with the greatest degree of biological and cultural diversity in the world (Challenger et al. 2009;Domínguez et al. 2009;Sarukhán et al. 2015). However, natural capital has suffered alarming deterioration due to the conversion of habitats to other soil uses and their consequent fragmentation and degradation (Challenger et al. 2009;Sarukhán et al. 2015). The natural vegetation cover of the country was only 54% in 1993, and 38% in 2002, of which approximately 50% was degraded (Challenger et al. 2009). More than 45% of soil degradation in Mexico is associated to human activity (Cotler et al. 2007), with rates higher in semiarid regions with low annual precipitation, which has a ranged of 260 to 600 mm (Zamora et al. 2020). Tula watershed is located in the states of Hidalgo, Tlaxcala, and Mexico, and it has two main climate types: temperate subhumid and semiarid (Köppen classification modified by García, 1964). Some of the semiarid climate regions are experiencing population growth, an expedition of mining concessions for the exploitation of construction materials, and agricultural and urban increase, resulting in the degradation of soils because of increasing demands of natural resources (INEGI, 2017;Zamora et al. 2020).
Given these serious alterations to Mexico's ecosystems, it is important to develop conservation and rehabilitation strategies. Rehabilitating ecosystems can prevent and reverse the loss of biodiversity and restore ecosystem services (Ceccon 2013). For degraded ecosystems to be restored, it is necessary to reintroduce species that can thrive in the edaphoclimatic conditions of these regions. Incorporating shrubs and succulents as Senna multiglandulosa, Sedum dendroideum and Aloe sp. under agroforestry practices can lead to sustainable soil management, as well as the rehabilitation of degraded areas. Among the factors considered in selecting these species were their ability to grow in semiarid to temperate climates, their resistance to drought and low temperatures, their ability to form soil, their effectiveness in erosion control, and their use in Mexico as alternative crops (Gazca and Benavides 2018; Kumar et al. 2022;Gutiérrez 2009).
Survival analysis considers data on the number of deaths of each species to calculate the likelihood of survival (Crawley 2007), which indicates the species' ability to adapt to different habitat conditions (Sigala et al. 2015). Cox proportional hazards regression model (Bradburn et al. 2003;Cox 1972) predicts a species' death probability at a given time t based on certain morphological variables. Using survival analysis and Cox proportional hazards modeling, this study estimated the probability of survival for each species and determined their risk of death based on their morphological characteristics. During the early stages of any field trial, it is crucial to evaluate the progress of the trial, which is why forest plantations should receive more attention at this stage in the research process (Aguilos et al. 2020).
This study aimed to evaluate the adaptation capacity of three non-native multipurpose species (Senna multiglandulosa, Sedum dendroideun, and Aloe sp.) in a degraded area with a semiarid temperate climate, located in the Tula watershed in Mexico, for restoring degrades soils through agroforestry by assessing their survival, growth of morphological variables, and death risk.

Geographical location of study area
The study area was the Tula watershed located in the states of Mexico, Hidalgo, and Tlaxcala ( Fig. 1 Zamora et al. 2020, Buendía-Espinoza et al. 2022). There are two main climates within the watershed: semiarid temperate (28% of the total area) and temperate sub-humid (72% of the total area). In the temperate semiarid climate, annual precipitation ranges from 500 to 600 mm per year. The predominant soils are haplic phaeozem and eutrophic cambisol, both of medium texture (loamy). In the sub-humid temperate climate, annual precipitation ranges from 600 to 1000 mm per year with rainfall in summer and the percentage of winter rainfall is 5 to 10.2% of the annual total. The precipitation in the driest month is less than 40 mm and the predominant soils are medium-textured haplic phaeozem and fine-textured pelic vertisol (INEGI 2017;García-CONABIO, 1998;INIFAP and CONABIO 1995;Zamora et al. 2020, Buendía-Espinoza et al. 2022).

Experimental site
Experimental site was established in Ejido San Felipe Teotitlán, municipality of San Felipe Teotitlán, located in the semiarid temperate zone within the Tula watershed at an elevation of 2 473 m and an area of 1 680 m 2 (56 × 30 m). The average temperature is 14 °C, with thermal oscillations in the monthly averages of 6.6 °C, the coldest month (January) has an average temperature of 11.1 °C, and the warmest month (May) has an average temperature of 17.7 °C. Rainfall ranges from 500 to 600 mm annually, of which approximately 83% falls between June and October (INEGI 2009).

Experimental design
An agroforestry design based on contour planting was applied to the field experiment as a means of increasing organic matter production, reducing runoff, and contributing to the creation of a microclimate (Ceccon 2013). For this purpose, six 30 m-long trenches, 30 cm deep and wide, were dug perpendicular to the slope of the land in order to reduce runoff and sediment dragging, where 156 plants of Aloe sp. were planted each meter. Each trench was planted with one Senna multiglandulosa plant every ten m to provide nitrogen to the soil and to protect it. A total of 24 plants were planted in the six trenches. An internal ridge was created between the trenches, where 288 plants of Sedum dendroide species were transplanted in a real frame arrangement of 2 m in order to establish a microclimate and protect the soil. During the period July 30 to August 11, 2020, excavation and transplanting activities were conducted (Fig. 2).

Target species
Based on edaphoclimatic conditions at the experimental site as well as their multipurpose potential, taking into account taxonomy, form, function, and uses, multipurpose plants were selected.
Senna multiglandulosa, known as retama, grows at altitudes between 2100 and 4000 m from Central Mexico to Guatemala (INECOL 1997). This species is able to thrive in semiarid to temperate climates, is drought and cold temperature tolerant, prefers loamy to sandy soils, even if they contain high levels of stoniness, tolerates poor soils and tepetose soils, requires medium to high moisture levels, and is useful for erosion control (Gazca and Benavides 2018).
Aloe sp., also known as aloe, is a semiperennial succulent plant of the Liliaceae family and has succulent leaves arranged in the shape of roses. It is native to Africa, but has spread throughout Mexico, thriving in a wide variety of habitats, including forests and deserts. Aloe sp. grows in areas with rainfall between 200 and 600 mm, slowing down in areas with rainfall below this amount, but is resistant to drought and thrives at altitudes ranging from 0 to 2500 m. This species grows best in loamy soil, but can also thrive in sandy soils or soils with very little organic matter and a pH of alkaline-neutral. The thermal regime of this plant is from 18 to 27 °C, and at temperatures below 5 °C, it shows damage (Pedroza and Gómez 2006).
Sedum dendroideum, also known as Siempreviva, is a shrubby plant in the family Crassulaceae, which is native to Mexico, Guatemala, and Central America (Pérez-Calix 2008). The species is found in semi-warm, semiarid and cold climates, and in a wide range of habitats, including rocky areas, crags, slopes, forests, and scrub with altitudes between 1350 and 2750 m above sea level (Gutiérrez 2009). Plants can grow without difficulty at maximum temperatures, but during the winter they must be kept dry in order to prevent damage since they store water in their succulent leaves (López 2019). When temperatures are low during frosty conditions, this species is capable of surviving temperatures as low as 0 °C (Toogood 2007).

Vegetation sampling and measurement of variables
During a period of 12 months beginning in September 2020 and ending in August 2021, the survival and morphological variables of the species were evaluated monthly. In each evaluation, the survival of each seedling was assessed, assigning a 1 or 0 to dead or live plants, respectively, and measuring the height, diameter, and canopy cover of the seedlings. The height of each plant (cm) was measured from the base of the plant to the apex of its dominant vertical shoot using a flexometer. The stem diameter (cm) was measured at ground level with a vernier. The coverage (cm 2 ) of the seedlings was calculated using the square method, which involves placing a 45.7 cm × 61 cm acrylic grid divided into 1 cm × 1 cm quadrants over the seedling and counting the number of squares occupied by the seedling. Squares with more than 50% covered were considered units.

Evaluation of species survival, vegetative growth, and mortality risk in the field
Based on the data collected on morphological variables and plant survival, the following indicators were calculated:

Survival estimation
The survival rate was calculated based on the Kaplan-Meier estimator, Ŝ KM , also known as the limit product estimator, a nonparametric statistical method that takes censoring into account when studying the period during which an event must occur (Kaplan and Meier 1958). This estimator can be used when the individual times of the censored and uncensored individuals are known, so that survival rates can be calculated when an individual dies or reaches the end of their follow-up period. The Kaplan-Meier estimator is defined as: where r t i indicates the live plants, d t i indicates the dead plants, and t i < t represents the time at which the measurement was made (Kaplan and Meier 1958;Rivas-Ruiz et al. 2014;Gallardo et al. 2016).
According to Kaplan and Meier (1958), survival function is defined as S(t) = P(t ≥ t), t ≥ 0 , where S(t) represents the probability of death occurring at time T equal to or greater than time t . The status of each plant (alive or dead) at the end of the evaluation period and its duration (months) were considered.
is the distribution function of T . S(t) decreases at a given rate at each instant, which represents the probability (risk) of an individual's death at that point (Kaplan and Meier 1958). SPSS v25 software was used for the analysis. Rodríguez-Echeverry and Leiton (2020) indicate that a survival rate exceeding 90% is excellent, between 70 and 89% is acceptable, between 50 and 69% is marginal, and less than 49% is inacceptable.

Estimation of death risks
Cox proportional hazards regression (Cox 1972) was used to examine the impact of planting site and production system on seedling morphology variables. The proportional hazards model used was as follows: where h i (t) is the risk of death for individual i at time t , derived by multiplying the unspecified baseline hazard function, h o (t) , by an exponential function of k covariates (Allison 1995). The model estimates a coefficient for each factor or covariate and tests the null hypothesis that H 0 ∶ = 0 using the Chi-square test. When such a coefficient is negative, it indicates that the risk of death is reduced with increasing covariates, whereas a positive value indicates the opposite (Williams 2008). SPSS v25 software was used for the analysis.

Estimation of vegetative growth
The height growth rate (HGR, cm month −1 ), diameter growth rate (DGR, mm month −1 ), and canopy growth rate (CGR, cm 2 month −1 ) were calculated using the following equations (Griscom et al. 2005): where H, D, andC represent the diameter, height, or canopy, respectively, at t 1 at the beginning of the experiment and at t 2 at the end of the experiment on a monthly basis. Species response at the planting site was assessed over four periods of analysis (Spring, summer, autumn, and winter).

Survival estimation
All species established in the experimental area survived 44.1% after a period of 12 months (Table 1). Considering the censored data, which consisted of plants that did not succumb to the death event during the analysis, the Aloe sp. had the highest survival estimate of 97.4%, followed by Senna multiglandulosa and Sedum dendroideum with survival rates of 17% each (Table 1). Both Sedum dendroideum and Senna multiglandulosa species achieved median survival rates of 50% at 7 and 8 months following establishment, respectively. The survival rate of Aloe sp. plants during the evaluation period did not exceed the established median survival limits, since 97.4% of the plants survived (Fig. 3).
As shown in Fig. 4, each species has a different risk of death. Sedum dendroideum recorded a death risk of 10% from October to December 2020; however, this increased to 50% from March to August 2021. The species Aloe had an estimated death rate of approximately 2% in April (month 8), which remained unchanged throughout the assessment. Senna multiglandulosa had an increased mortality risk in April 2021.
Cox regression results indicate that the predictor variable final height (finhei) has a positive effect on survival in the species Senna multiglandulosa; that is, a 1 cm increase in height reduces the risk of death by 1.5% (1-0.985), provided all other variables remain constant (Table 2).

Estimation of vegetative growth
Tables 5, 6, 7, and 8 present the average growth rates of plant height, stem diameter, and canopy cover for each species during the evaluation period. During the autumn cycle, Sedum dendroideum and Aloe sp. both had negative growth rates in height and canopy cover, with − 0.29 cm month −1 and − 7.37 cm 2 month −1 and − 0.53 cm month −1 and − 4.26 cm 2 month −1 , respectively. Senna multiglandulosa demonstrated a positive growth rate in height (0.87847 cm month −1 ) and a negative growth rate in canopy cover, the latter of which was the highest at − 22.80 cm 2 month −1 . Senna multiglandulosa and Sedum dendroideum showed positive growth rates at the stem diameter level (Table 5).
During the second period, which corresponds to the winter cycle, all three species showed negative growth rates for height, stem diameter, and negative canopy cover, with Senna multiglandulosa and Sedum dendroideum showing the greatest decreases in stem diameter and canopy cover, with − 0.005067 cm month −1 and − 10.44 cm 2 month −1 , respectively. Species Aloe sp. showed the greatest reduction in height, down − 1.15235 cm month −1 (Table 6).
In the springer cycle, which corresponds to the third period, all three species showed positive growth rates for height, stem diameter, and canopy cover, with the exception of Sedum dendroideum, which again displayed a decrease in height (− 0.44872 cm month −1 ). Senna multiglandulosa had the highest canopy cover growth rate with 59.5392 cm 2 month −1 in this period (Table 7).
In the fourth period, which corresponds to the summer cycle, all three species had positive growth rates in height, stem diameter, and canopy cover, with the exception of Sedum dendroideum, which showed a decrease of − 0.02442 cm month −1 in stem growth rate. Aloe  sp. and Senna multiglandulosa achieved the greatest growth rates in terms of height and canopy cover, measuring 0.94737 cm month −1 and 46.8333 cm 2 month −1 , respectively.

Survival estimation
Plant survival is related to rainfall regime (del Campo 2002;Alloza 2003), so its scarcity during the early period following planting is the risk factor that generates the greatest risk for survival (Alloza and Vallejo 1999), particularly in those environments where vegetative activity can begin as early as the winter months. The average survival rate of the three species at the experimental site was below acceptable levels ( < 49% , Rodríguez-Echeverry and Leiton 2020). According to Omary (2011) and Vallejo et al. (2012), plant survival can be affected by rainfall, slope, soil properties such as humidity, temperature, pH, electrical conductivity, nutrients, as well as planting or transplanting techniques (Ortega et al. 2006), water availability optimization processes used during plant establishment, and functional traits of plants developed in water-limited environments (Hernández et al. 2010). In this study, all three species were transplanted on the same soil substrate, so any differences in their survival rates may have been the result of the incidence of various factors which contribute to the establishment stage of plants, such as environmental conditions, plant management, morphology, physiology, and genetics (South 2000;Chen and Klinka 1998). The best example of adaptive power here is aloe, which is a Crassulacean Acid Metabolism (CAM) plant that maintains a series of regulatory mechanisms in its metabolism to adapt to seasonal changes (Honda et al 2000). Under stress conditions, aloe maintains this photosynthetic pathway of CAM plants and switches to the C3 metabolic pathway once conditions improve. Also, aloe has a thick cuticle that makes it more resistant to cold temperatures (Pedroza and Gómez 2006). Sedum dendroideum is also a CAM plant, but it has a thin cuticle which does not allow it to withstand cool season changes, whereas Senna multiglandulosa is a C3 plant, whose photosynthetic efficiency varies depending on the species, the phenological stage, and the aridity of the environment (Rundel et al. 1999). Senna multiglandulosa and Sedum dendroideum showed survival rates below acceptable levels ( ∼ 17% < 49% , Rodríguez-Echeverry and Leiton 2020), while Aloe sp. displayed excellent survival rates ( ∼ 97 ≥ 90%, Rodríguez-Echeverry and Leiton 2020), indicating adaptability to the site's edaphoclimatic conditions. Senna multiglandulosa and Sedum dendroideum showed the highest mortality during the months of March, April and May, likely due to excessive evapotranspiration and rainfall distribution in the study area, which averages 500 to 600 mm of rainfall per year, 83% of which falls between June and October (INEGI 2009). Margolis and Brand (1990) and Navarro et al. (2006) suggest that reforestation may be successful if plants are able to survive the establishment phase, which usually lasts less than two years, however the first year is crucial, as the plants begin to absorb water and nutrients from the soil once they come into contact with it. Therefore, factors that affect the plant's water status at the time of planting are critical to its survival (Burdett 1990;Heiskanen and Rikala 2000). Sigala et al. (2015) report that Cox proportional hazards regression can be used to estimate the proportion of species mortality risk based on morphological variables in conjunction with the site-specific edaphoclimatic conditions. Senna multiglandulosa's final height (finhei) is the most important predictor of survival, that is, it negatively affects mortality, which is consistent with Silva et al. (2009) and Stovall et al. (2019), who reported that juvenile trees have smaller canopy covers to allocate more resources to growth in height but are more likely to survive. Sedum dendroideum's final diameter (findia) was the most significant factor in reducing its risk of death, followed by its canopy cover growth rate (CGR) and initial canopy cover (inidos). The final diameter (findia) reduced mortality by approximately 55%, which is consistent with Fontes (1999), who suggests that some species allocate more resources to developing their diameter as they develop, whereas early in their development, their diameter is greater than that of other species. The growth rate of canopy cover reduced (CGR) the risk of death by approximately 31%, which is in agreement with Moser et al. (2007), who assert that canopy leaf tissue, which catches sunlight and transpires water, plays an important role in establishing and growing new plants, as well as determining their geographic and temporal distribution. For aloe species, height growth rate (HGR) was the most important predictor variable in reducing the risk of death by 93.9%, followed by initial canopy cover (inican) by 2.4%. These findings are consistent with the findings of Valdecantos (2001), who states that height increases are typically positive under medium-high survival rates indicating good climatic conditions. Navarro et al. (2006) point out that a plantation's success is determined by the plant's ability to grow and develop according to seasonal conditions and the species' capabilities, which are traditionally measured in terms of survival and growth. During the autumn cycle, when the rainy season ended and the cold season began, Senna multiglandulosa achieved a high survival rate and a positive growth rate in height (HGR: 0.87847 cm month −1 ) and diameter (DGR: 0.046528 cm month −1 ), however, its canopy cover area (CGR: -22.8056 cm 2 month −1 ) declined. Seedlings of tree and shrub species generally develop narrow canopy covers to maximize their height and, to a lesser extent, diameter growth, and in response to light demands (Martínez-Sánchez 2008). Under subhumid conditions with summer rainfall in Chapultepec forest, Mexico City, Mexico, Gazca and Benavides (2018) also reported positive growth rates in height and diameter for Senna multiglandulosa. Moreover, Rico and Bachman (2006), as well as Terrones et al. (2004), recognize this species as fast-growing.

Estimation of vegetative growth
During this same period, Sedum dendroideum species showed a high survival rate and a positive diameter growth rate (DGR: 0.01888825 cm month −1 ), but negative height (HGR: − 0.292636 cm month −1 ) and canopy cover growth rates (CGR: − 7.377292 cm 2 month −1 ). Negative canopy growth rates, which contribute to initial survival, but subsequent growth is positively related to lower canopy cover levels (Yang et al. 2011), may have affected plant diameter and height by decreasing photosynthesis (Moser et al. 2007;Olivas et al. 2013). During this time period, Aloe sp. had a high survival rate, but a negative growth rate in height (HGR: − 0.53082 cm month −1 ) and canopy area (CGR: − 4.26496 cm 2 month −1 ). Although the planting season for Aloe sp. extends into August, its negative growth rates in height and canopy cover could be a result of early frosts, which are a major problem during this period due to the plant's slow rooting habit, low growth tendencies and susceptible leaves (Pedroza and Gomez 2006).
During the winter period, Senna multiglandulosa again maintained an excellent survival rate, while Sedum dendroideum obtained an acceptable survival rate, but their growth rates were negative for height (− 0.80694 and − 0.832245 cm month −1 , respectively), diameter (− 0.00042 and − 0.005067 cm month −1 ), and canopy cover (− 10.4444 and − 4.357831 cm 2 month −1 ). Gazca and Benavides (2018) also reported negative growth rates for Senna multiglandulosa during the autumn cycle in Chapultepec forest, Mexico City, Mexico, under subhumid conditions and summer rainfall. Negative growth rates in canopy cover could have led to negative growth rates in heights and diameters since plants' establishment and growth depend on leaf tissue that performs photosynthesis (Moser et al. 2007;Olivas et al. 2013). In plants, leaf area is affected by hydrological, biogeochemical, and biophysical processes (Peduzzi et al. 2012), as well as by solar radiation interception, photoassimilate conversion efficiency, use of water and nitrogen, and temperature regulation (de Freitas et al. 2007; Barrios and Cobo 2004). The canopy's negative growth rate may also be attributed to the high levels of stress caused by the distribution of rainfall at the study site, which resulted in, at best, regrowth, while in severe cases, a significant reduction in height and diameter (Burdett 1990;Heiskanen and Rikala 2000).
During the same period, Aloe sp. plants also survived well, but their growth rates were negative in height (− 1.15235 cm month −1 ) and canopy cover (− 4.88889 cm 2 month −1 ). These results are in agreement with Pedroza and Gómez (2006) who explain that Aloe plants are restricted in their growth and development by frost, mainly at temperatures below 7 °C, causing severe damage at temperatures below 0 °C.
During the seasons of spring and summer, when the rainy season begins, Senna multiglandulosa y Sedum dendroideum obtained survival rates below acceptable levels ( ∼ 17% ) at the end of the spring season. Navarro et al. (2006) and Vallejo et al. (2012) note that high mortality of established plants in semiarid areas and degraded dry sites is primarily determined by rainfall amount and spatial distribution, evapotranspiration rates, the length of the growing season, and soil characteristics, including depth and moisture retention capacity. Nevertheless, Senna multiglandulosa showed positive growth rates in terms of height (HGR: 8.94444 cm month −1 -spring, and HGR: 2.75 cm month −1 -summer), diameter (DGR: 0.143611 cm month −1 -spring, and DGR: 0.103333 cm month −1 -summer), and canopy cover (CGR: 59.5392 cm 2 month −1 -spring, and CGR: 46.8333 cm 2 month −1 -summer). Growth rates for height, stem diameter, and canopy cover might have been positive because the rainy seasons have been linked to evapotranspiration (Navarro et al. 2006) and plant survival (Padilla and Pugnaire 2007), which facilitate a deeper root system (Grantz et al. 1998;Padilla and Pugnaire 2007).
During the spring period, Sedum dendroideum showed positive growth rates for diameter (DGR: 0.053577 cm month −1 ) and canopy cover (CGR: 0.76403 cm 2 month −1 ), but negative growth rates for height (HGR: − 0.44872 cm month −1 ). However, during the summer, Sedum dendroideum exhibited positive growth rates in terms of height (HGR: 0.50239 cm month −1 ) and canopy cover (CGR: 2.82863 cm 2 month −1 ) and slightly negative growth rates in terms of diameter (DGR: − 0.02442 cm month −1 ). Growth rates in diameter and height may have been positive due to canopy growth, as some species invest more resources in diameter growth during development compared to their beginnings, and as species progress through more advanced stages of development, they become taller (Fontes 1999). During these two seasons, Aloe sp. plants showed positive growth rates in height (HGR: 1.81499 cm month −1 -spring, and HGR: 0.94737 cm month −1 -summer), as well as in canopy cover (CGR: 3.705972 cm 2 month −1 -spring, and CGR: 6.730263 cm 2 month −1 -summer), corroborating what Pedroza and Gómez (2006) said that Aloe sp. plants resume their growth gradually after leaving the frost season with rainfall under 400 mm.

Conclusions
Based on this study, it's possible to determine whether a species can adapt to edaphoclimatic conditions where it's established, taking into account survival, vegetative growth, and death risk, so that agroforestry can be recommended as a land management method on degraded soils that benefits society as a whole as well as rural development. Among the three plant species evaluated, Aloe sp showed the greatest ability to adapt to the edaphoclimatic conditions of the study area. The seasonal conditions of the site affect the growth rates of morphological variables of the species evaluated in varying degrees. Climate conditions during autumn and winter affected the height, stem diameter, and canopy diameter of Senna multiglandulosa and Sedum dendroideum species. Plant height and canopy cover of Aloe sp. species were affected by the winter weather conditions. Each species' morphological variables determined its death risk. For Senna multiglandulosa species, the final height variable reduced the death risk while for Sedum dendroideum species, the final diameter, initial canopy cover, and growth rate in canopy cover reduced the death risk. Variable height growth rate and initial canopy cover in Aloe species reduced death risk.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.