Environmental changes in biodiversity hotspot ecosystems of
South Ecuador: RESPonse and feedback effECTs (FOR2730)

Abstract:

An understanding of the genetic diversity and structure of plant species is essential in order to comprehend the degree of biodiversity loss and to develop successful restoration programs. Handroanthus is an important genus that presents one of the most valuable timbers of South America. Handroanthus chrysanthus is an important species distributed in Central and South America. Microsatellite markers are not previously developed for this species. Ten microsatellites for Handroanthus chrysanthus developed using high-throughput sequencing are presented here. The usefulness of these microsatellite loci for the genetic analysis of subspecies H. chrysanthus subsp. chrysanthus (distributed in coastal dry forests) and subspecies H. chrysanthus subsp. meridionalis (distributed in premontane moist forests) is analyzed. At least eight polymorphic microsatellites are useful for each subspecies, seven of which can be used in both subspecies.

Abstract:

In the past two decades, the Amazon-exposed, tropical montane rain forests in south Ecuador experienced increasing deposition of reactive N mainly from Amazonian forest fires, episodic Ca and Mg inputs from Saharan dust, and a low but constant P deposition from unknown sources. To explore the response of this tropical, perhumid ecosystem to nutrient inputs, we established in 2007 a Nutrient Manipulation Experiment (NUMEX). Since 2008, we have applied 50 kg ha−1 year−1 of N as urea, 10 kg ha−1 year−1 of P as NaH2PO4·H2O, 50 kg ha−1 year−1 of N + 10 kg ha−1 year−1 of P and 10 kg ha−1 year−1 of Ca as CaCl2·H2O in a randomized block design at 2000 m a.s.l. in a natural forest of the south Ecuadorian Andes. Previous studies have shown that alkali and alkaline earth metals had beneficial effects on the functioning of N and P co-limited tropical forests occurring on acidic soils. Therefore, we determined the response of all major aqueous ecosystem fluxes of K, Ca, Mg and Na to nutrient amendments, to understand how increasing atmospheric deposition would affect their cycling in the future. Additions of N and P decreased K leaching from the organic layer and in the mineral soil, thus tightening K cycling. This suggests that increasing future N and P availability may result in K limitation in the long term. The leaching of Ca and Mg from the canopy increased in response to amendments of N and P and we observed an enhanced uptake of these nutrients also if Ca was amended alone. Although N was applied as urea, acidity of soil solutions and leaching of K, Ca, Mg and Na did not increase following separate N amendments. In spite of the acid soils and of its low cation-exchange competitivity, Na included in the P fertilizer was only partly leached from the organic layer. We suggest that it was probably required to cover an unmet Na demand of the soil fauna. Our results demonstrate the major role in the functioning of the tropical montane forests played by K, Ca and Mg as potential future growth-limiting elements and increasingly required nutrients in response to rising N and P availability, while they also support the importance of Na as a functional element in these ecosystems.

Abstract:

Summary Community trait assembly in highly diverse tropical rainforests is still poorly understood. Based on more than a decade of field measurements in a biodiversity hotspot of southern Ecuador, we implemented plant trait variation and improved soil organic matter dynamics in a widely used dynamic vegetation model (the Lund-Potsdam-Jena General Ecosystem Simulator, LPJ-GUESS) to explore the main drivers of community assembly along an elevational gradient. In the model used here (LPJ-GUESS-NTD, where NTD stands for nutrient-trait dynamics), each plant individual can possess different trait combinations, and the community trait composition emerges via ecological sorting. Further model developments include plant growth limitation by phosphorous (P) and mycorrhizal nutrient uptake. The new model version reproduced the main observed community trait shift and related vegetation processes along the elevational gradient, but only if nutrient limitations to plant growth were activated. In turn, when traits were fixed, low productivity communities emerged due to reduced nutrient-use efficiency. Mycorrhizal nutrient uptake, when deactivated, reduced net primary production (NPP) by 61–72% along the gradient. Our results strongly suggest that the elevational temperature gradient drives community assembly and ecosystem functioning indirectly through its effect on soil nutrient dynamics and vegetation traits. This illustrates the importance of considering these processes to yield realistic model predictions.

Abstract:

90% of the 1.2 billion people living with less than 1 USD/day, depend on the utilization of forests. The dry forests of the Tumbes-Piura-region at the coast of Southwest Ecuador are one of the 25 global hotspots of biodiversity with an immense number of endemic species. Nonetheless, they are endangered by anthropogenic impact. In this study we surveyed 89 households of five different districts to analyze the extent and importance of the different forest uses in the region. The results showed that forest utilization is practised by 97% of the households. With a share of 31% it contributes the 2nd highest amount to the mean total household income, only surpassed by livestock farming (36%). The most important use of the forests is silvopasture, which represents not only an important source of cash income, but also a substantial contribution to subsistence.

Abstract:

The amount of carbon stored in deadwood is equivalent to about 8per cent of the global forest carbon stocks1. The decomposition of deadwood is largely governed by climate2–5 with decomposer groups—such as microorganisms and insects—contributing to variations in the decomposition rates2,6,7. At the global scale, the contribution of insects to the decomposition of deadwood and carbon release remains poorly understood7. Here we present a field experiment of wood decomposition across 55 forest sites and 6 continents. We find that the deadwood decomposition rates increase with temperature, and the strongest temperature effect is found at high precipitation levels. Precipitation affects the decomposition rates negatively at low temperatures and positively at high temperatures. As a net effect—including the direct consumption by insects and indirect effects through interactions with microorganisms—insects accelerate the decomposition in tropical forests (3.9% median mass loss per year). In temperate and boreal forests, we find weak positive and negative effects with a median mass loss of 0.9per cent and −0.1per cent per year, respectively. Furthermore, we apply the experimentally derived decomposition function to a global map of deadwood carbon synthesized from empirical and remote-sensing data, obtaining an estimate of 10.93.2petagram of carbon per year released from deadwood globally, with 93per cent originating from tropical forests. Globally, the net effect of insects may account for 29per cent of the carbon flux from deadwood, which suggests a functional importance of insects in the decomposition of deadwood and the carbon cycle.

Abstract:

Tropical dry forests are an intricate ecosystem with special adaptations to periods of drought. Arbuscular mycorrhizal fungi (AMF) are essential for plant survival in all terrestrial ecosystems but might be of even greater importance in dry forests as plant growth is limited due to nutrient and water deficiency during the dry season. Tropical dry forests in Ecuador are highly endangered, but studies about AMF communities are scarce. We investigated the AMF community of a premontane semi-deciduous dry forest in South Ecuador during the dry season. We estimated AMF diversity, distribution, and composition of the study site based on operational taxonomic units (OTUs) and compared the results to those from the tropical montane rainforest and páramo in South Ecuador. OTU delimitation was based on part of the small ribosomal subunit obtained by cloning and Sanger sequencing. Nearly all OTUs were Glomeraceae. The four frequent OTUs were Glomus, and comparison with the MaarjAM database revealed these to be globally distributed with a wide range of ecological adaptations. Several OTUs are shared with virtual taxa from dry forests in Africa. Ordination analysis of AMF communities from the tropical dry and montane rainforests in South Ecuador revealed a unique AMF community in the dry forest with only few overlapping OTUs. Most OTUs that were found in both dry and rainforests and on the two continents were globally distributed Glomus.

Abstract:

Screening species-rich communities for the variation in functional raits along environmental gradients may help understanding the abiotic drivers of plant performance in a mechanistic way. We investigated tree leaf trait variation along an elevation gradient (1000–3000 m) in highly diverse neotropical montane forests to test the hypothesis that elevational trait change reflects a trend toward more conservative resource use strategies at higher elevations, with interspecific trait variation decreasing and trait integration increasing due to environmental filtering. Analysis of trait variance partitioning across the 52 tree species revealed for most traits a dominant influence of phylogeny, except for SLA, leaf thickness and foliar Ca, where elevation was most influential. The community-level means of SLA, foliar N and Ca, and foliar N/P ratio decreased with elevation, while leaf thickness and toughness increased. The contribution of intraspecific variation was substantial at the community level in most traits, yet smaller than the interspecific component. Both within-species and between-species trait variation did not change systematically with elevation. High phylogenetic diversity, together with small-scale edaphic heterogeneity, cause large interspecific leaf trait variation in these hyper-diverse Andean forests. Trait network analysis revealed increasing leaf trait integration with elevation, suggesting stronger environmental filtering at colder and nutrient-poorer sites.

Abstract:

Theory predicts positive effects of species richness on the productivity of plant communities through complementary resource use and facilitative interactions between species. Results from manipulative experiments with tropical tree species indicate a positive diversity–productivity relationship (DPR), but the existing evidence from natural forests is scarce and contradictory. We studied forest aboveground productivity in more than 80 humid tropical montane oldgrowth forests in two highly diverse Andean regions with large geological and topographic heterogeneity and related productivity to tree diversity and climatic, edaphic and stand structural factors with a likely influence on productivity. Main determinants of wood production in the perhumid study regions were elevation (as a proxy for temperature), soil nutrient (N, P and base cation) availability and forest structural parameters (wood specific gravity, aboveground biomass). Tree diversity had only a small positive influence on productivity, even though tree species numbers varied largely (6–27 species per 0.04 ha). We conclude that the productivity of highly diverse Neotropical montane forests is primarily controlled by thermal and edaphic factors and stand structural properties, while tree diversity is of minor importance.

Abstract:

A new satellite-based algorithm for rainfall retrieval in high spatio-temporal resolution fo Ecuador is presented. The algorithm relies on the precipitation information from the Integrated Multi-SatEllite Retrieval for the Global Precipitation Measurement (GPM) (IMERG) and infrared (IR) data from the Geostationary Operational Environmental Satellite-16 (GOES-16). It wa developed to (i) classify the rainfall area (ii) assign the rainfall rate. In each step, we selected the most important predictors and hyperparameter tuning parameters monthly. Between 19 April 2017 and 30 November 2017, brightness temperature derived from the GOES-16 IR channels and ancillary geo-information were trained with microwave-only IMERG-V06 using random forest (RF). Validation was done against independent microwave-only IMERG-V06 information not used for training. The validation results showed the new rainfall retrieva technique (multispectral) outperforms the IR-only IMERG rainfall product. This offers using the multispectral IR data can improve the retrieval performance compared to single-spectrum IR approaches. The standard verification scored a median Heidke skill score of ~0.6 for the rain area delineation and R between ~0.5 and ~0.62 for the rainfall rate assignment, indicating uncertainties for Andes’s high elevation. Comparison of RF rainfall rates in 2 km2 resolution with daily rain gauge measurements reveals the correlation of R = ~0.33.

Abstract:

Background: Since the 1990’s, afforestation programs in the páramo have been implemented to offset carbon emissions through carbon sequestration, mainly using pine plantations. However, several studies have indicated that after the establishment of pine plantations in grasslands, there is an alteration of carbon pools including a decrease of the soil organic carbon (SOC) pool. The aim of this study is to investigate the impact of the establishment of pine plantations on the carbon stocks in different altitudes of the páramo ecosystem of South Ecuador. Results: At seven locations within an elevational gradient from 2780 to 3760 m a.s.l., we measured and compared carbon stocks of three types of land use: natural grassland, grazed páramo, and Pinus patula Schlltdl. & Cham. plantation sites. For a more accurate estimation of pine tree carbon, we developed our own allometric equations. There were significant (p < 0.05) differences between the amounts of carbon stored in the carbon pools aboveground and belowground for the three types of land use. In most of the locations, pine plantations revealed the highest amounts of aboveground and belowground carbon (55.4 and 6.9 tC/ha) followed by natural grassland (23.1 and 2.7 tC/ha) and grazed páramo sites (9.1 and 1.5 tC/ha). Concerning the SOC pools, most of the locations revealed significant lower values of plantations’ SOC in comparison to natural grassland and grazed páramo sites. Higher elevation was associated with lower amounts of pines’ biomass. Conclusions: Even though plantations store high amounts of carbon, natural páramo grassland can also store substantial amounts above and belowground, without negatively affecting the soils and putting other páramo ecosystem services at risk. Consequently, plans for afforestation in the páramo should be assessed case by case, considering not only the limiting factor of elevation, but also the site quality especially affected by the type of previous land use.

Abstract:

Morphological traits provide the interface between species and their environment. For example, body size affects the fitness of individuals in various ways. Yet especially for ectotherms, the applicability of general rules of interspecific clines of body size and even more so of other morphological traits is still under debate. Here we tested relationships between elevation (as a proxy for temperature) and productivity with four ecologically relevant morphological traits of orthopteran assemblages that are related to fecundity (body size), dispersal (wing length), jumping ability (hind femur length), and predator detection (eye size). We measured traits of 160 orthopteran species that were sampled along an extensive environmental gradient at Mt. Kilimanjaro (Tanzania), spanning elevations from 790 to 4,410 m above sea level (a.s.l.) with different levels of plant productivity. For traits other than body size, we calculated the residuals from a regression on body length to estimate the variation of traits irrespective of body size. Bayesian analyses revealed that mean body size of assemblages, as well as the means of relative wing length, hind femur length, and eye size, decreased with increasing elevation. Body size and relative eye size also decreased with increasing productivity. Both phylogenetic relationships, as well as species‐specific adaptations, contributed to these patterns. Our results suggest that orthopteran assemblages had higher fecundity and better dispersal and escape abilities, as well as better predator detection at higher temperatures (low elevations) than at low temperatures (high elevations). Large body sizes might be advantageous in habitats with low productivity because of a reduced risk of starvation. Likewise, large eye size might be advantageous because of the ability to detect predators in habitats with low vegetation cover, where hiding possibilities are scarce. Our study highlights that changes in temperature and productivity not only lead to interspecific changes in body size but are also related to independent changes of other morphological traits that influence the ecological fit of organisms in their environment.

Abstract:

Habitat loss and landscape fragmentation are important drivers of changes in biodiversity. In fragmented landscapes, bird species are able to use multiple forest patches that may therefore share an important portion of the regional biodiversity. In turn, these patches are linked through their shared bird diversity (i.e. species-habitat networks). Identifying the importance of nodes (e.g. forest patches) in species-habitat networks is increasingly important to improve conservation planning. Within this network approach, the relative importance of patches for birds can be identified via centrality indices-measures of the magnitude of shared bird diversity of patches across the entire network (i.e. patch centrality). Here, we tested for changes in patch centrality within bird species-patch networks of two habitat guilds, forest specialists and generalists, in relation to patch area, patch shape irregularity and within-patch habitat characteristics across 15 Polylepis woodland patches in a páramo landscape of southern Ecuador. Patch centrality for forest specialists decreased with greater influence of surrounding páramo plants, i.e. an increasing proportion of bunch-grasses and small shrubs, in the within-patch habitat and was unaffected by either patch area or patch shape irregularity. On the other hand, patch centrality for generalists was positively influenced by patch shape irregularity but was unaffected by patch area or the influence of surrounding páramo plants in the within-patch habitat. Patch centrality reveals that the relative importance of Polylepis woodlands lies in their habitat quality. Forest specialists are dependent on mature Polylepis woodland patches, while generalists benefit from the natural irregular shape of the woodlands. Finally, a species-habitat network approach facilitates the recognition of important Polylepis patches and their characteristics for conservation of the Andean bird community at a landscape scale.

Abstract:

The Random Forest (RF) algorithm, a decision-tree-based technique, has become a promising approach for applications addressing runoff forecasting in remote areas. This machine learning approach can overcome the limitations of scarce spatio-temporal data and physical parameters needed for process-based hydrological models. However, the influence of RF hyperparameters is still uncertain and needs to be explored. Therefore, the aim of this study is to analyze the sensitivity of RF runoff forecasting models of varying lead time to the hyperparameters of the algorithm. For this, models were trained by using (a) default and (b) extensive hyperparameter combinations through a grid-search approach that allow reaching the optimal set. Model performances were assessed based on the R2, %Bias, and RMSE metrics. We found that: (i) The most influencing hyperparameter is the number of trees in the forest, however the combination of the depth of the tree and the number of features hyperparameters produced the highest variability-instability on the models. (ii) Hyperparameter optimization significantly improved model performance for higher lead times (12- and 24-h). For instance, the performance of the 12-h forecasting model under default RF hyperparameters improved to R2 = 0.41 after optimization (gain of 0.17). However, for short lead times (4-h) there was no significant model improvement (0.69 &lt; R2 &lt; 0.70). (iii) There is a range of values for each hyperparameter in which the performance of the model is not significantly affected but remains close to the optimal. Thus, a compromise between hyperparameter interactions (i.e., their values) can produce similar high model performances. Model improvements after optimization can be explained from a hydrological point of view, the generalization ability for lead times larger than the concentration time of the catchment tend to rely more on hyperparameterization than in what they can learn from the input data. This insight can help in the development of operational early warning systems.

Abstract:

Tropical mountain ecosystems are threatened by climate and land-use changes. Their diversity and complexity make projec- tions how they respond to environmental changes challenging. A suitable way are trait-based approaches, by distinguishing between response traits that determine the resistance of species to environmental changes and efect traits that are relevant for species’ interactions, biotic processes, and ecosystem functions. The combination of those approaches with land surface models (LSM) linking the functional community composition to ecosystem functions provides new ways to project the response of ecosystems to environmental changes. With the interdisciplinary project RESPECT, we propose a research framework that uses a trait-based response-efect-framework (REF) to quantify relationships between abiotic conditions, the diversity of functional traits in communities, and associated biotic processes, informing a biodiversity-LSM. We apply the framework to a megadiverse tropical mountain forest. We use a plot design along an elevation and a land-use gradient to collect data on abiotic drivers, functional traits, and biotic processes. We integrate these data to build the biodiversity- LSM and illustrate how to test the model. REF results show that aboveground biomass production is not directly related to changing climatic conditions, but indirectly through associated changes in functional traits. Herbivory is directly related to changing abiotic conditions. The biodiversity-LSM informed by local functional trait and soil data improved the simulation of biomass production substantially. We conclude that local data, also derived from previous projects (platform Ecuador), are key elements of the research framework. We specify essential datasets to apply this framework to other mountain ecosystems.

Abstract:

Tropical Mountain Forest (TMF) provides important ecological functions like evapotranspiration (ET) that supplies moisture and energy to the atmosphere. ET observations are scarce and difficult to accomplish particularly in areas of high heterogeneity where TMF are. Remote sensing (RS) allows to quantify and to determine ET spatial variation at the landscape level. Detail imaginary improves high spatial variability retrieval. Thought the greater detail introduces cast shadows by trees which hamper image interpretation. The objective of this study is to characterize ET estimation for the TMF of the southern Ecuadorian Andes by combining meteorological data with high-resolution satellite images. Shadows from high resolution images were masked out by applying focal statistics. The analysis included two meteorological periods typical of the area; a wet period when rain prevails and a dry period when precipitation is more sporadic. The reference evapotranspiration (ET0) was calculated using the FAO-Penman Montheid method by applying data obtained from an automatic weather station. The enhanced vegetation index (EVI) was derived from 2 m resolution WorldView2 satellite images. Results showed a lower ET mean value during the wet period: 1.54 mm day−1 compared to 2.37 mm day−1. Two forest types, differentiated from its structural composition and topographical position (ravine and ridge), marked ET spatial variation. Ravine forest that has a more dense and closed canopy showed higher ET values for both meteorological conditions. A comparison between ET estimations and ET field measurements from a scintillometer device showed a good agreement (coefficient of correlation r = 0.89) that proves the validity of the method. This study demonstrates that the application of high spatial resolution improves ET estimation in TMF especially when shadows are removed. Also, emphasizes the importance of analysing spatial heterogeneity to properly assess ecosystem water flux terms.

Abstract:

Forest plantations with exotic species in southern Ecuador have mostly been characterized as having negative externalities in both ecological and economic terms. After ten years of research in mountain forests in southern Ecuador on restoration and reforestation, several native tree species with good growth responses (e.g. Handroanthus chrysanthus, Cedrela montana and Juglans neotropica) have been identified in open field conditions. Some species (e.g. Podocarpus oleifolius and P. sprucei) were able to grow under the shelter provided by Pinus and Alnus trees.

Abstract:

Ratios of stable oxygen isotopes in tree rings (δ18O) are a valuable proxy for reconstructing past climates. Such reconstructions allow us to gain better knowledge of climate dynamics under different (eg warmer) environmental conditions, which also forms the basis for effective risk management. The latter aspect is particularly relevant for our study site on the western flanks of the Andes in Southern Ecuador, since the region is frequently affected by droughts and heavy precipitation events during the rainy season (January to April), leading to enormous social and economic losses. In particular, we focus on precipitation amounts and moisture source regions as they are known to influence the δ18O signature of tree rings. Moisture source regions are based on 240 h backward trajectories that were calculated with the trajectory model LAGRANTO for the rainy seasons 2008 to 2017. A moisture source diagnostic was applied to the air parcel pathways. The resulting moisture source regions were analyzed by calculating composites based on precipitation amounts, season, and calendar year. The precipitation amounts were derived from data of a local Automatic Weather Station (AWS). The analysis confirms that our study site receives its moisture both, from the Atlantic and the Pacific Oceans. Heavy precipitation events are linked to higher moisture contributions from the Pacific, and local SST anomalies along the coast of Ecuador are of higher importance than those off the coast toward the central Pacific. Moreover, we identified increasing moisture contributions from the Pacific over the course of the rainy season. This change and also rain amount effects are detectable in preliminary data of δ18O variations in tree rings of Bursera graveolens. These signatures can be a starting point for investigating atmospheric and hydroclimatic processes, which trigger δ18O variations in tree rings, more extensively in future studies.

Abstract:

(1) We investigated the variation in tree fine root traits and their functional diversity along a local topographic gradient in a Neotropical montane forest to test if fine root trait variation along the gradient is consistent with the predictions of the root economics spectrum on a shift from acquisitive to conservative traits with decreasing resource supply. (2) We measured five fine root functional traits in 179 randomly selected tree individuals of 100 species and analyzed the variation of single traits (using Bayesian phylogenetic multi‐level models) and of functional trait diversity with small scale topography. (3) Fine roots exhibited more conservative traits (thicker diameters, lower specific root length and nitrogen concentration) at upper slope compared to lower slope positions, but the largest proportion of variation (40‐80%) was explained by species identity and phylogeny. Fine root functional diversity decreased towards the upper slopes. (4) Our results suggest that local topography and the related soil fertility and moisture gradients cause considerable small‐scale variation in fine root traits and functional diversity along tropical mountain slopes, with conservative root traits and greater trait convergence being associated with less favorable soil conditions due to environmental filtering. We provide evidence of a high degree of phylogenetic conservation in fine root traits.

Abstract:

The response of organic carbon (C) concentrations in ecosystem solutions to environmental change affects the release of dissolved organic matter (DOM) from forests to surface and groundwaters. We determined the total organic C (TOC) concentrations (filtered <1–7 μm) and the ratios of TOC/dissolved organic nitrogen (DON) concentrations, electrical conductivity (EC), and pH in all major ecosystem solutions of a tropical montane forest from 1998 to 2013. The forest was located on the rim of the Amazon basin in Ecuador and experienced increasing numbers of days with >25°C, decreasing soil moisture, and rising nitrogen (N) deposition from the atmosphere during the study period. In rainfall, throughfall, mineral soil solutions (at the 0.15- and 0.30-m depths), and streamflow, TOC concentrations and fluxes decreased significantly from 1998 to 2013, while they increased in stemflow. TOC/DON ratios decreased significantly in rainfall, throughfall, soil solution at the 0.15-m depth, and streamflow. Based on Δ14C values, the TOC in rainfall and mineral soil solutions was 1 year old and that of litter leachate was 10 years old. The pH in litter leachate decreased with time, that in mineral soil solutions increased, while those in the other ecosystem solutions did not change. Thus, reduced TOC solubility because of lower pH values cannot explain the negative trends in TOC concentrations in most ecosystem solutions. The increasing TOC concentrations and EC in stemflow pointed at an increased leaching of TOC and other ions from the bark. Our results suggest an accelerated degradation of DOM, particularly of young DOM, associated with the production of N-rich compounds simultaneously to changing climatic conditions and increasing N availability. Thus, environmental change increased the CO2 release to the atmosphere but reduced DOM export to surface and groundwater.

Abstract:

Litterfall is the most easily and most frequently measured part of net primary productivity (NPP) of forests. It has been shown that litterfall accounts for about one third of total NPP and thus serves as a proxy for the total productivity of forests. Moreover, litterfall carries nutrients from the forest canopy to the soil and therefore is also a major vector of nutrient cycling. I reviewed the published literature about litterfall rates and chemical properties in Andean forests and found reports from 44 forest sites, which I evaluated together with unpublished data from 12 sites in a lower montane forest in Ecuador. I found many more reports from tropical (52 sites) than from temperate Andean forests (4 sites). In the humid tropical north Andes, litterfall showed a hump-shaped elevational distribution. It increased from premontane to lower montane forests and decreased to upper montane forests. The tropical lower montane forest had a similar productivity than tropical lowland forests. The litterfall of the temperate southern beech forests was similar to that of the tropical upper montane forests. The C/N and C/P ratios of litterfall decreased with increasing elevation, while the N/P ratios were not correlated with elevation. This illustrates that the forests become increasingly nutrient efficient with increasing elevation, while there is no indication of a general change in the kind of nutrient limitation. There was a negative correlation between litterfall and soil organic layer thickness (r = −0.61, p < 0.001) illustrating that the organic matter input via litterfall is a less important driver of organic matter accumulation on top of the mineral soil than other, mainly abiotic properties including temperature and soil waterlogging. My evaluation suggests that there are systematic relationships between abiotic conditions and litterfall, which could be used to predict litterfall in the Andes. However, the data coverage particularly of the southern Andes (Bolivia, Chile, Argentina), the Andean dry forests, and the widespread tree plantations is poor. The observed elevational influence on litterfall in the humid tropical Andes suggests that the forest productivity will likely respond to climate change driving the vegetation belts to higher elevation with an unknown overall effect on C sequestration of these forests.

Abstract:

Interest in atmosphere-surface flux modeling over the mountainous regions of the globe has increased recently, with a major focus on the prediction of water, carbon and other functional indicators in natural and disturbed conditions. However, less research has been centered on exploring energy fluxes (net radiation; sensible, latent and soil heat) and actual evapotranspiration (ETa) over the Neotropical Andean biome of the páramo. The present study assesses the implementation and parameterization of a state-of-art Land-Surface Model (LSM) for simulation of these fluxes over two representative páramo catchments of southern Ecuador. We evaluated the outputs of the LSM Community Land Model (CLM ver. 4.0) with (i) ground-level flux observations from the first (and highest) Eddy Covariance (EC) tower of the Northern Andean páramos; (ii) spatial ETa estimates from the energy balance-based model METRIC (based on Landsat imagery); and (iii) derived ETa from the closure of the water balance (WB). CLM’s energy predictions revealed a significant underestimation on net radiation, which impacts the sensible and soil heat fluxes (underestimation), and delivers a slight overestimation on latent heat flux. Modeled CLM ETa showed acceptable goodness-of-fit (Pearson R = 0.82) comparable to ETa from METRIC (R = 0.83). Contrarily, a poor performance of ETa WB was observed (R = 0.46). These findings provide solid evidence on the CLM’s accuracy for the ETa modeling, and give insights in the selection of other ETa methods. The study contributes to a better understanding of ecosystem functioning in terms of water loss through evaporative processes, and might help in the development of future LSMs’ implementations focused on climate / land use change scenarios for the páramo.