Forest dynamics of a sub-xerophilous... FOREST DYNAMICS OF A SUB-XEROPHILOUS VEGETATION FORMATION IN CENTRAL PERU - CHANCHAMAYO VALLEY, PERU

1 Received on 31.01.2018 accepted for publication on 11.10.2018. 2 Universidade Nacional Agrária La Molina, Faculdade de Ciências Florestais, Lima Perú. E-mail: <soniapalacios@lamolina.edu.pe>, <rmontenegro@lamolina.edu.pe> and <reynel@lamolina.edu.pe>. 3 Center for Conservation and Sustainability, Smithsonian Conservation Biology Institute, Lima Perú. E-mail: <linaresr@si.edu>. *Corresponding author.


1.INTRODUCTION
The study of forest dynamics in the tropics has mainly focused on lowland rain forests in the Amazon revealing that their structure and dynamics are changing (Nebel et al., 2000;Lewis et al., 2004;Feeley et al., 2011).The species richness in these formations is mainly explained by their forest dynamics rather than by any other factor, showing that disturbance at smaller scales is functionally related with the species richness in a community (Phillips et al., 1994).In terms of tree growth, previous research has suggested that the availability of resources stimulates the net primary production, accelerating the dynamics of basal area and increasing the above-ground biomass (Lewis et al., 2004).However, the specific causes and consequences of such increases are still being discussed.
The study of ecological processes in forests, such as the dynamics of mortality and recruitment rates, and the changes in the floristic composition, can give us insight to and help understand the changing processes of the woody component.Studying the changes in composition improves not only our understanding about the ecology of tropical forests and their response to local and regional disturbances, but also our capacity to predict how future global changes can influence some vital services provided by these ecosystems (Feeley et al., 2011).
Examination of the changing patterns of species composition and their relation to carbon fertilization predictions, of past perturbations, of recovery from extreme drought caused by El Niño and of long-term climatic changes, showed that those changes increased the relative proportions of drought-tolerant species.These increases were determined by the occurrence of species along a soil humidity gradient and through a wider rainfall gradient (Feeley et al., 2011).Furthermore, additional research suggested that the compositional transition of tropical forests to more seasonal forest types is a consequence of increasingly dryer conditions which limit the distribution of species (Esquivel-Muelbert et al., 2017).
An optimal scenario to study the contrast between the humid flora and the flora from seasonally dry formations occurs in the central Peruvian rain forest in the Chanchamayo valley.This valley, part of the Perené River basin and located between 600 to 1200 m above sea level, harbors a sub-xerophilous floristic formation composed of small fragments embedded in a more extensive humid forest matrix.Since much of this valley is exposed to frequent slash and burn and land-use change events, it is important to study and monitor the remaining dry forest fragments over time to understand their dynamics and their contribution to diversity and conservation in these changing premontane landscapes in Central Peru.Especially considering the increase in intensity and length of dry seasons in the Amazon basin (Marengo et al., 2011).
Field observations in the study area suggested that the density of woody individuals was increasing.This could be the result of an increased recruitment rate, but not of mortality.This study sets out to prove this hypothesis.In addition, we also expected changes in the forest structure.These changes were also analyzed from a floristic composition perspective to explore species turnover.

Study area
The study area is in the Perené River Valley, San Ramon District, Chanchamayo Province, Junin Region, in Central Peru.Average annual rainfall is between 1500 and 3000 mm, whilst average annual temperature is between 15 and 19 ºC.Soils are alluvial and support mainly rainforests that have 90 to 147 tree species/ha (Marcelo-Peña and Reynel, 2014).The valley also presents scattered seasonally dry vegetation formations with species that are characteristic of Neotropical dry forests and/or savannas, such as Astronium fraxinifolium (Anacardiaceae), Curatella americana (Dilleniaceae), Machaerium pilosum (Fabaceae), Roupala montana (Proteaceae), Pseudobombax marginatum (Malvaceae), Luehea paniculata (Malvaceae), Dilodendron bipinnatum (Sapindaceae), among others.Previous studies show that these species are absent or rare in the rainforest areas of the valley (Palacios-Ramos, 2017).Many of these areas are burned regularly for crop production.
One of us (SP) stablished a permanent plot for this drier type of vegetation within in a three-hectare area in 2009 (UTM WGS84: 18L 462779 E, 8772597 N).The last recorded fires of anthropic origin in the adjacent areas were in 2005 and 2009.

Methodology
This study compares tree census data from 2009 when the permanent plot (100x60m; composed of 15 sub-plots of 20x20m) was stablished (Palacios-Ramos Forest dynamics of a sub-xerophilous... and Reynel, 2011) with the second census in 2015, after six years with three months (6.25 years).
Both censuses included the assessment of all trees with DBH (diameter at breast height, 1.30m above ground) equal or above 10 cm, following the protocol of Phillips et al. (2009).The second census included: (i) recruits (individuals that reache or passed 10cm of DBH in the 2015 census) and (ii) dead individuals.Trees that could not be referred to a previously sampled specimen were collected to be identified later in the herbarium.The samples were processed and deposited in the Universidad Nacional Agraria La Molina Herbarium (MOL: Woody plants collection).

Data Analysis
We constructed a species accumulation curve using each 20x20m sub-plot within the permanent plot (100 permutation for each year we assessed) to show species richness values of the plot according to the size of the sample unit.To assess species diversity, we calculated Fisher's Alpha index, as it allows to compare data from sample units with different abundances.In both cases we used the R programming environment (R Development Core Team, 2017) and the ´vegan´ package (Oksanen et al., 2017).
We further calculated woody plant cover values by adding the relative densities (percentage of number of trees) and the relative dominances (percentage of basal area) to analyze the horizontal structure and representativeness of the species.
We analyzed forest dynamics by using mortality and recruitment rates.The mortality rate per annum was calculated according to the exponential growth model in continuous time (Nebel et al., 2000).
where: Tm = annual mortality rate expressed in percentage; No = number of individuals recorded in the first census; Nm = number of individuals that died within the t time interval; t = time interval t expressed in years; ln = natural logarithm.
We calculated the recruitment rate per species by following an exponential function of population increase (Phillips et al., 1994;Nebel et al., 2000).
where: Tr = recruitment rate or resettlement in percentage; No = number of individuals recorded during plots setup; Nr = number of individuals resettled within time interval t; Nm = number of individuals that died during time interval t; t = time interval t in years.ln = natural logarithm.

Diversity, floristic composition and structure
The species accumulation curve for the 2015 recensus shows a stabilization trend (Figure 1).This trend is similar to the estimated values calculated from the 2009 data.Fisher's Alpha diversity was 7.39 in 2009 and 7.5 in 2015.The number of woody individuals per hectare was 642, with an increase of 21% compared to the value estimated in 2009 (532 trees per hectare).In 2009 we recorded 22 families, and added three new families in 2015, each one with one species and only one individual (Annonaceae: Guatteria sp., Vochysiaceae: Vochysia mapirensis and Lacistemataceae: Lacistema aggregatum).Six of the seven species that were represented by only one tree in 2009 kept only one individual in 2015.
Abundance distribution by diametric categories in 2015 showed more individuals in the lower diametric categories, just as in 2009, resulting in EQ-1 EQ-2  Forest dynamics of a sub-xerophilous...
an inverted-j trend.Two of the species with higher cover values kept the same distribution in both censuses (Figure 2).

Mortality and recruitment
Eleven species presented a total of 14 dead individuals in the second census.We also recorded 82 new individuals belonging to 18 species, three of them were species not recorded in the 2009 census.The relatively even distribution of dead individuals and recruits is slightly skewed towards the northwest of the plot (Figure 3).The recruitment of Physocalymma scaberrimum (17), Heteropterys intermedia (15) and Sapium glandulosum (14) represented 56% of all recruits.
There is a notorious difference between the high recruitment rates (8.82%) and the low mortality rates (0.72%).18 of the 27 species registered in 2009 showed changes in their dynamics, primarily influenced by the high rate of recruitment (Table 2).(Weberbauer, 1945;Garcia-Villacorta, 2009;

PALACIOS-RAMOS et al.
Linares -Palomino et al., 2011).All these belong to the ten species with the highest plant cover values in the plot, having maintained their representativeness since 2009 (Table 1).In our study area, they have only been registered in the sub-xerophilous formations (Palacios-Ramos and Reynel, 2011;Palacios-Ramos, 2017), except for Sapium glandulosum which is considered an ecological generalist species (Linares-Palomino et al., 2011) and also found in the adjacent more humid areas.
Studies have suggested that rainfall seasonality regulate the alpha diversity of Amazonian forests.However, the variation of diversity at local scales is the result of other factors (Ter Steege et al., 2003).For instance, even under high local humidity regimes, local edaphic conditions such as porous soils could increase hydrological stress and decrease diversity (Esquivel-Muelbert et al., 2017).A recent study in our area of interest compared soil diversity and characteristics of the different forest types and showed that the sub-xerophilous forest grew on soils different from those with humid forest growing on it.These different soils allow species to endure conditions with low organic material content, high levels of pH and interchangeable acidity (Romero, 2017).These conditions could be a limiting factor for the colonization of the adjacent humid forest species.
Another local factor that could influence the dynamics of species is fire.Fire exclusion trials in the Brazilian Cerrado showed that rainforest species are capable of colonizing savannas and compete with the local species assemblage adapted to fire (Geiger et al., 2011), suggesting that a closed canopy forest could develop on the low nutrient soil of the Brazilian Cerrado (Simon and Pennington, 2012).Time is of utmost relevance in this forest -savanna alternation.The transition from saplings to adult trees or from the savanna to forest requires long but unusual intervals without fire.The opposite transitions are favored by prolonged drought and when non-altered forests become flammable (Hoffmann et al., 2012).These statements are relevant for our study area considering that since 2005 and 2009 there have been no more fires.During those events, fires reached the plot border impacting the understory vegetation and the foliage of adjacent trees.However, many of the species in our plot, such as Roupala montana, Byrsonima spicata, Luehea paniculata, Maprounea guianense and Pseudobombax marginatum present adaptations to fire (Palacios-Ramos, 2017).And species such as Curatella americana require fire for reproduction (Penninton and Ratter, 2006).The absence of fire during these six years has not induced major changes in the horizontal structure of the forest (Table 1).The present status could be part of a cyclic succession between more open formations and closed canopy forests, although further censuses would be necessary to confirm it.
Despite the high overall number of recruits and low number of dead individuals, there is an increase of individuals in diametric classes higher than 10cm (Figure 1).However, unlike other forests in early successional states (Torres et al., 2012;Muñoz et al., 2014) this increase stabilizes above the 25cm diameter class.The smaller diametric classes kept the high number of trees and showed an inverted-j trend.This structure is characteristic of a forest with continuous recruitment and is particularly prevalent in some of the most important species of this forest (Figure 2), providing them whith good possibilities to remain in the plot over a longer term.
Curatella americana and Pseudobombax marginatum are species characteristic of the semideciduous Chiquitano forest and of the Brazilian Cerrados (Furley et al., 1988).They are also present in the subxerophilous formations in Central Peru, usually in areas with an open canopy and with predominance of grasses (Palacios-Ramos and Reynel, 2011).Long-term term Forest dynamics of a sub-xerophilous... permanence in these sub-xerophilous formations could be compromised since our plot data showed no recruitment.In addition, the increase in the density of trees per hectare of more than 20% between censuses provides an increasingly closed canopy.
Although the values we report for the mortality and recruitment rates could be strongly influenced by past locally provoked fires, they also show a floristic community with members adapted to those fires.In the dry forest, the floristic composition can be more sensitive to climatic changes associated to rainfall rather than past disturbances during relatively short timescales (Enquist and Enquist, 2011).The intensity of mortality (mortality rate per annum is lower than 5%) and the distribution of dead individuals per species demonstrate that it is a background mortality.This is part of the processes that occur in the forest community, and could be influenced by auto-shading dependence and radiation capture, the position of trees in the canopy, the successional stage of the forest (Lugo and Scatena, 1996;Melo et al., 2016), the age or size of trees, or the mortality of each species (Lugo and Scatena, 1996).The reported mortality rate has a spatial distribution with no relationship to the limits where past fires occurred (NE vertex of the plot) (Figure 3).This means that the influence of these events does not appear to be evident in the trees we evaluated (with DBH bigger or equal than 10cm).At the same time, the slight decrease in the spatial distribution of recruitment towards the NW of the plot does not seem to be related to past fire events (NE vertex of the plot) (Figure 3).
The species with the highest recruitment rate constitute only 11% of the trees in the plot, while their recruitment constitutes less than 25% of the total recruits registered in 2015.In this context, the influence of these species in this forest could be questioned and we require a longer assessment period to value their relevance.Furthermore, although the most important species showed low recruitment rates, they dominated the forest during the period assessed.

Implication for conservation
Monitoring and assessment to detect patterns of change in the floristic composition of the populations that harbor these sub-xerophylous formations is out reach for the present study.However, the existence of these formations cannot be rejected during conservation planning of the pre-montane forests of the Chanchamayo Valley and the Perené Watershed.Especially when we consider the economic value of many of the species in these forests and their adaptation to dry conditions.
The constraint of most tree taxa in the western Neotropics to the more humid areas shows low generalized tolerance to dry conditions and low capability to acclimatize to them (Esquivel-Muelbert et al., 2017).In the present context of increase in the density and length of dry season in the Amazon Watershed (Marengo et al., 2011) and the increase of temperature (Jiménez- Muñoz et al., 2013), the results we are showing could be relevant.Particularly considering that the species that occupy a dominant position in the study area regarding dominance and relative density, are species from seasonally dry areas that have morphological adaptations to survive fire (Palacios-Ramos and Reynel, 2011;Palacios-Ramos, 2017), and are already growing in soils with low content of organic material, high levels of pH and interchangeable acidity (Romero, 2017).These species could be useful for restoration efforts of similar forest systems degraded by anthropic intervention and that are dispersed throughout the valley.

5.CONCLUSIONS
The population dynamics of the sub-xerophilous formation is mainly influenced by a high recruitment rate (3.82%) and a low mortality rate (0.72%).These values are different from other studies in Neotropical seasonally dry formations.These dynamics have not produced noticeable changes in the horizontal structure nor the floristic composition.The maintenance of these conditions suggest they are favorable for the existence of this forest.
The conservation of this sub-xerophilous forest will contribute to their recovery.And if the impact of the pressure over the forest fragments decreases, it can become a source of information for research and monitoring of the seasonally dry forests in Central Peru.

Figure 2 -
Figure 2 -Abundance distribution by diameter class of selected and pooled species during the 2009 and 2015 censuses.Figura 2 -Distribuição por classes diamétricas de 2009 e 2015.
, floristic composition and forest structure Four of the five most important species in the plot (Heteropterys intermedia, Physocalymma scaberrimum, Tachigalii peruviana and Machaerium hirtum) are characteristic of seasonally dry forests and/or wooded savannas (Palacios-Ramos and Reynel, 2011; Palacios-Ramos, 2017).Luehea paniculata, Roupala montana and Byrsonima spicata are common in seasonally dry formations

Figure 3 -
Figure 3 -Spatial distribution of dead and recruited individuals in the permanent plot.Figura 3 -Distribuição espacial de mortos e recrutas na parcela.