1. |
- Castaño, Nicolas, et al.
(author)
-
Contrasting drivers of aboveground woody biomass and aboveground woody productivity in lowland forests of Colombia
- 2024
-
In: Ecography. - 0906-7590 .- 1600-0587. ; 2024
-
Journal article (peer-reviewed)abstract
- The relative importance of abiotic and biotic factors in shaping forest biomass stocks and fluxes remains a controversial issue. Here, using data gathered from 39 1 ha plots located in flooded and terra firme mature tropical lowland forests of the Amazon and Orinoquia regions of Colombia, we evaluated the importance of climate, soil fertility, and flooding, as well as tree taxonomic/phylogenetic diversity and forest structural properties, in determining the aboveground biomass stocks (AGB; Mg ha−1) and aboveground woody productivity (AWP; Mg ha−1 year−1). Using information-theoretic multimodel inference and variance partitioning we found that forest structural features such as the number of trees with diameter at breast height ≥ 70 cm, and wood density, are the main drivers of variation in AGB. However, taxonomic diversity also contributes to AGB because it is associated with more large trees in these forests. In contrast, the key drivers of AWP in these forests were soil P and Mg concentrations, with no significant effects of diversity indices. These findings emphasize the need to include major soil cations other than N and P (e.g. Mg) in experimental studies to improve our understanding about the extent to which soil fertility can modulate increases in forest AWP due to climate change. Terra firme forests had higher AGB stocks than flooded forests, but both had similar AWP; and we found similar results for the drivers of AGB and AWP between flooded and terra firme forests. Our results provide limited evidence for strong effects of plant diversity on AGB or AWP. Therefore, we call for caution on generalizations of nature-based initiatives aiming to preserve diversity based on maximizing carbon stocks and productivity, due to the complex nature of the processes controlling carbon accumulation and carbon fluxes in tropical forests.
|
|
2. |
- Davies, Stuart J., et al.
(author)
-
ForestGEO: Understanding forest diversity and dynamics through a global observatory network
- 2021
-
In: Biological Conservation. - : Elsevier BV. - 0006-3207. ; 253
-
Journal article (peer-reviewed)abstract
- ForestGEO is a network of scientists and long-term forest dynamics plots (FDPs) spanning the Earth's major forest types. ForestGEO's mission is to advance understanding of the diversity and dynamics of forests and to strengthen global capacity for forest science research. ForestGEO is unique among forest plot networks in its large-scale plot dimensions, censusing of all stems ≥1 cm in diameter, inclusion of tropical, temperate and boreal forests, and investigation of additional biotic (e.g., arthropods) and abiotic (e.g., soils) drivers, which together provide a holistic view of forest functioning. The 71 FDPs in 27 countries include approximately 7.33 million living trees and about 12,000 species, representing 20% of the world's known tree diversity. With >1300 published papers, ForestGEO researchers have made significant contributions in two fundamental areas: species coexistence and diversity, and ecosystem functioning. Specifically, defining the major biotic and abiotic controls on the distribution and coexistence of species and functional types and on variation in species' demography has led to improved understanding of how the multiple dimensions of forest diversity are structured across space and time and how this diversity relates to the processes controlling the role of forests in the Earth system. Nevertheless, knowledge gaps remain that impede our ability to predict how forest diversity and function will respond to climate change and other stressors. Meeting these global research challenges requires major advances in standardizing taxonomy of tropical species, resolving the main drivers of forest dynamics, and integrating plot-based ground and remote sensing observations to scale up estimates of forest diversity and function, coupled with improved predictive models. However, they cannot be met without greater financial commitment to sustain the long-term research of ForestGEO and other forest plot networks, greatly expanded scientific capacity across the world's forested nations, and increased collaboration and integration among research networks and disciplines addressing forest science.
|
|
3. |
- Gómez-Correa, Luisa F., et al.
(author)
-
Canopy tree mortality depends on the proportion of crown exposed to sunlight, but this effect varies with species' wood density
- 2023
-
In: Biotropica. - 0006-3606 .- 1744-7429. ; 55, s. 1136-1147
-
Journal article (peer-reviewed)abstract
- Understanding what drives changes in tree mortality as well as the covariates influencing trees' response is a research priority to predict forest responses to global change. Here, we combined drone photogrammetry and ground-based data to assess the influence of crown exposure to light (relative to total crown area), growth deviations (relative to conspecifics), tree size, and species' wood density (as a surrogate for light-demanding and shade-tolerant life-history strategies) on the mortality of 984 canopy trees in an Amazon terra firme forest. Trees with lower wood density were less prone to die when their proportion of crown was more exposed to sunlight, but this relationship with relative crown exposure weakened and slightly reversed as wood density increased. Trees growing less than their species average had higher mortality, especially when the species' wood density decreased. The role of wood density in determining the survival of canopy trees under varying light conditions indicates differential responses of light-demanding versus shade-tolerant species. Our results highlight the importance of accounting for life-history strategies, via plant functional types, in vegetation dynamic models aiming to predict forest demography under a rapidly changing climate. Abstract in Spanish is available with online material.
|
|
4. |
- Hudson, Lawrence N, et al.
(author)
-
The database of the PREDICTS (Projecting Responses of Ecological Diversity In Changing Terrestrial Systems) project
- 2017
-
In: Ecology and Evolution. - : John Wiley & Sons. - 2045-7758. ; 7:1, s. 145-188
-
Journal article (peer-reviewed)abstract
- The PREDICTS project-Projecting Responses of Ecological Diversity In Changing Terrestrial Systems (www.predicts.org.uk)-has collated from published studies a large, reasonably representative database of comparable samples of biodiversity from multiple sites that differ in the nature or intensity of human impacts relating to land use. We have used this evidence base to develop global and regional statistical models of how local biodiversity responds to these measures. We describe and make freely available this 2016 release of the database, containing more than 3.2 million records sampled at over 26,000 locations and representing over 47,000 species. We outline how the database can help in answering a range of questions in ecology and conservation biology. To our knowledge, this is the largest and most geographically and taxonomically representative database of spatial comparisons of biodiversity that has been collated to date; it will be useful to researchers and international efforts wishing to model and understand the global status of biodiversity.
|
|
5. |
- Hülsmann, Lisa, et al.
(author)
-
Latitudinal patterns in stabilizing density dependence of forest communities
- 2024
-
In: Nature. - 0028-0836 .- 1476-4687. ; 627, s. 564-571
-
Journal article (peer-reviewed)abstract
- Numerous studies have shown reduced performance in plants that are surrounded by neighbours of the same species1,2, a phenomenon known as conspecific negative density dependence (CNDD)3. A long-held ecological hypothesis posits that CNDD is more pronounced in tropical than in temperate forests4,5, which increases community stabilization, species coexistence and the diversity of local tree species6,7. Previous analyses supporting such a latitudinal gradient in CNDD8,9 have suffered from methodological limitations related to the use of static data10–12. Here we present a comprehensive assessment of latitudinal CNDD patterns using dynamic mortality data to estimate species-site-specific CNDD across 23 sites. Averaged across species, we found that stabilizing CNDD was present at all except one site, but that average stabilizingCNDD was not stronger toward the tropics. However, in tropical tree communities, rare and intermediate abundant species experienced stronger stabilizing CNDD than did common species. This pattern was absent in temperate forests, which suggests that CNDD influences species abundances more strongly in tropical forests than it does in temperate ones13. We also found that interspecific variation in CNDD, which might attenuate its stabilizing effect on species diversity14,15, was high but not significantly different across latitudes. Although the consequences of these patterns for latitudinal diversity gradients are difficult to evaluate, we speculate that a more effective regulation of population abundances could translate into greater stabilization of tropical tree communities and thus contribute to the high local diversity of tropical forests.
|
|
6. |
- Leite, Melina de Souza, et al.
(author)
-
Major axes of variation in tree demography across global forests
- 2024
-
In: Ecography. - 0906-7590 .- 1600-0587.
-
Journal article (peer-reviewed)abstract
- The future trajectory of global forests is closely intertwined with tree demography, and a major fundamental goal in ecology is to understand the key mechanisms governing spatio-temporal patterns in tree population dynamics. While previous research has made substantial progress in identifying the mechanisms individually, their relative importance among forests remains unclear mainly due to practical limitations. One approach to overcome these limitations is to group mechanisms according to their shared effects on the variability of tree vital rates and quantify patterns therein. We developed a conceptual and statistical framework (variance partitioning of Bayesian multilevel models) that attributes the variability in tree growth, mortality, and recruitment to variation in species, space, and time, and their interactions – categories we refer to as organising principles (OPs). We applied the framework to data from 21 forest plots covering more than 2.9 million trees of approximately 6500 species. We found that differences among species, the species OP, proved a major source of variability in tree vital rates, explaining 28–33% of demographic variance alone, and 14–17% in interaction with space, totalling 40–43%. Our results support the hypothesis that the range of vital rates is similar across global forests. However, the average variability among species declined with species richness, indicating that diverse forests featured smaller interspecific differences in vital rates. Moreover, decomposing the variance in vital rates into the proposed OPs showed the importance of unexplained variability, which includes individual variation, in tree demography. A focus on how demographic variance is organized in forests can facilitate the construction of more targeted models with clearer expectations of which covariates might drive a vital rate. This study therefore highlights the most promising avenues for future research, both in terms of understanding the relative contributions of groups of mechanisms to forest demography and diversity, and for improving projections of forest ecosystems.
|
|
7. |
- Medina-Vega, José A., et al.
(author)
-
Tropical tree ectomycorrhiza are distributed independently of soil nutrients
- 2024
-
In: Nature Ecology and Evolution. - 2397-334X. ; 8, s. 400-410
-
Journal article (peer-reviewed)abstract
- Mycorrhizae, a form of plant–fungal symbioses, mediate vegetation impacts on ecosystem functioning. Climatic effects on decomposition and soil quality are suggested to drive mycorrhizal distributions, with arbuscular mycorrhizal plants prevailing in low-latitude/high-soil-quality areas and ectomycorrhizal (EcM) plants in high-latitude/low-soil-quality areas. However, these generalizations, based on coarse-resolution data, obscure finer-scale variations and result in high uncertainties in the predicted distributions of mycorrhizal types and their drivers. Using data from 31 lowland tropical forests, both at a coarse scale (mean-plot-level data) and fine scale (20 × 20 metres from a subset of 16 sites), we demonstrate that the distribution and abundance of EcM-associated trees are independent of soil quality. Resource exchange differences among mycorrhizal partners, stemming from diverse evolutionary origins of mycorrhizal fungi, may decouple soil fertility from the advantage provided by mycorrhizal associations. Additionally, distinct historical biogeographies and diversification patterns have led to differences in forest composition and nutrient-acquisition strategies across three major tropical regions. Notably, Africa and Asia’s lowland tropical forests have abundant EcM trees, whereas they are relatively scarce in lowland neotropical forests. A greater understanding of the functional biology of mycorrhizal symbiosis is required, especially in the lowland tropics, to overcome biases from assuming similarity to temperate and boreal regions.
|
|
8. |
- Needham, Jessica F., et al.
(author)
-
Demographic composition, not demographic diversity, predicts biomass and turnover across temperate and tropical forests
- 2022
-
In: Global Change Biology. - : Wiley. - 1354-1013 .- 1365-2486. ; 28, s. 2895-2909
-
Journal article (peer-reviewed)abstract
- The growth and survival of individual trees determine the physical structure of a forest with important consequences for forest function. However, given the diversity of tree species and forest biomes, quantifying the multitude of demographic strategies within and across forests and the way that they translate into forest structure and function remains a significant challenge. Here, we quantify the demographic rates of 1961 tree species from temperate and tropical forests and evaluate how demographic diversity (DD) and demographic composition (DC) differ across forests, and how these differences in demography relate to species richness, aboveground biomass (AGB), and carbon residence time. We find wide variation in DD and DC across forest plots, patterns that are not explained by species richness or climate variables alone. There is no evidence that DD has an effect on either AGB or carbon residence time. Rather, the DC of forests, specifically the relative abundance of large statured species, predicted both biomass and carbon residence time. Our results demonstrate the distinct DCs of globally distributed forests, reflecting biogeography, recent history, and current plot conditions. Linking the DC of forests to resilience or vulnerability to climate change, will improve the precision and accuracy of predictions of future forest composition, structure, and function.
|
|
9. |
- Perng, Bo Hao, et al.
(author)
-
Mapping distribution of woody plant species richness from field rapid assessment and machine learning
- 2024
-
In: Taiwania. - 0372-333X. ; 69, s. 1-15
-
Journal article (peer-reviewed)abstract
- Sustainable forest management needs information on spatial distribution of species richness. The objectives of this study were to understand whether knowledge, method, and effort of a rapid assessment affected accuracy and consistency in mapping species richness. A simulation study was carried out with nine 25–50 ha census plots located in tropical, subtropical, and temperate zones. Each forest site was first tessellated into non-overlapping cells. Rapid assessment was conducted in all cells to generate a complete coverage of proxies of the underlying species richness. Cells were subsampled for census, where all plant individuals were identified to species in these census cells. An artificial neural network model was built using the census cells that contain rapid assessment and census information. The model then predicted species richness of cells that were not censused. Results showed that knowledge level did not improve the accuracy and consistency in mapping species richness. Rapid assessment effort and method significantly affected the accuracy and consistency. Increasing rapid assessment effort from 10 to 40 plant individuals could improve the accuracy and consistency up to 2.2% and 2.8%, respectively. Transect reduced accuracy and consistency by up to 0.5% and 0.8%, respectively. This study suggests that knowing at least half of the species in a forest is sufficient for a rapid assessment. At least 20 plant individuals per cell is recommended for rapid assessment. Lastly, a rapid assessment could be carried out by local communities that are familiar with their forests; thus, further supporting sustainable forest management.
|
|
10. |
- Piponiot, Camille, et al.
(author)
-
Distribution of biomass dynamics in relation to tree size in forests across the world
- 2022
-
In: New Phytologist. - : Wiley. - 0028-646X .- 1469-8137. ; 234, s. 1664-1677
-
Journal article (peer-reviewed)abstract
- Tree size shapes forest carbon dynamics and determines how trees interact with their environment, including a changing climate. Here, we conduct the first global analysis of among-site differences in how aboveground biomass stocks and fluxes are distributed with tree size. We analyzed repeat tree censuses from 25 large-scale (4–52 ha) forest plots spanning a broad climatic range over five continents to characterize how aboveground biomass, woody productivity, and woody mortality vary with tree diameter. We examined how the median, dispersion, and skewness of these size-related distributions vary with mean annual temperature and precipitation. In warmer forests, aboveground biomass, woody productivity, and woody mortality were more broadly distributed with respect to tree size. In warmer and wetter forests, aboveground biomass and woody productivity were more right skewed, with a long tail towards large trees. Small trees (1–10 cm diameter) contributed more to productivity and mortality than to biomass, highlighting the importance of including these trees in analyses of forest dynamics. Our findings provide an improved characterization of climate-driven forest differences in the size structure of aboveground biomass and dynamics of that biomass, as well as refined benchmarks for capturing climate influences in vegetation demographic models.
|
|
11. |
- Zuleta, Daniel, 1990, et al.
(author)
-
Damage to living trees contributes to almost half of the biomass losses in tropical forests
- 2023
-
In: Global Change Biology. - : Wiley. - 1354-1013 .- 1365-2486. ; 29, s. 3409-3420
-
Journal article (peer-reviewed)abstract
- Accurate estimates of forest biomass stocks and fluxes are needed to quantify global carbon budgets and assess the response of forests to climate change. However, most forest inventories consider tree mortality as the only aboveground biomass (AGB) loss without accounting for losses via damage to living trees: branchfall, trunk breakage, and wood decay. Here, we use ~151,000 annual records of tree survival and structural completeness to compare AGB loss via damage to living trees to total AGB loss (mortality + damage) in seven tropical forests widely distributed across environmental conditions. We find that 42% (3.62 Mg ha−1 year−1; 95% confidence interval [CI] 2.36–5.25) of total AGB loss (8.72 Mg ha−1 year−1; CI 5.57–12.86) is due to damage to living trees. Total AGB loss was highly variable among forests, but these differences were mainly caused by site variability in damage-related AGB losses rather than by mortality-related AGB losses. We show that conventional forest inventories overestimate stand-level AGB stocks by 4% (1%–17% range across forests) because assume structurally complete trees, underestimate total AGB loss by 29% (6%–57% range across forests) due to overlooked damage-related AGB losses, and overestimate AGB loss via mortality by 22% (7%–80% range across forests) because of the assumption that trees are undamaged before dying. Our results indicate that forest carbon fluxes are higher than previously thought. Damage on living trees is an underappreciated component of the forest carbon cycle that is likely to become even more important as the frequency and severity of forest disturbances increase.
|
|
12. |
- Zuleta, Daniel, 1990, et al.
(author)
-
Drought-induced mortality patterns and rapid biomass recovery in a terra firme forest in the Colombian Amazon
- 2017
-
In: Ecology. - 0012-9658 .- 1939-9170. ; 98, s. 2538-2546
-
Journal article (peer-reviewed)abstract
- Extreme climatic events affecting the Amazon region are expected to become more frequent under ongoing climate change. In this study, we assessed the responses to the 2010 drought of over 14,000 trees ≥10 cm dbh in a 25 ha lowland forest plot in the Colombian Amazon and how these responses varied among topographically defined habitats, with tree size, and with species wood density. Tree mortality was significantly higher during the 2010–2013 period immediately after the drought than in 2007–2010. The post-drought increase in mortality was stronger for trees located in valleys (+243%) than for those located on slopes (+67%) and ridges (+57%). Tree-based generalized linear mixed models showed a significant negative effect of species wood density on mortality and no effect of tree size. Despite the elevated post-drought mortality, aboveground biomass increased from 2007 to 2013 by 1.62 Mg ha1 yr1 (95% CI 0.80–2.43 Mg ha1 yr1). Biomass change varied among habitats, with no significant increase on the slopes (1.05, 95% CI 0.76 to 2.85 Mg ha1 yr1), a significant increase in the valleys (1.33, 95% CI 0.37–2.34 Mg ha1 yr1), and a strong increase on the ridges (2.79, 95% CI 1.20–4.21 Mg ha1 yr1). These results indicate a high carbon resilience of this forest to the 2010 drought due to habitat-associated and interspecific heterogeneity in responses including directional changes in functional composition driven by enhanced perfor-mance of drought-tolerant species that inhabit the drier ridges.
|
|
13. |
- Zuleta, Daniel, 1990, et al.
(author)
-
Importance of topography for tree species habitat distributions in a terra firme forest in the Colombian Amazon
- 2020
-
In: Plant and Soil. - : Springer Science and Business Media LLC. - 0032-079X .- 1573-5036. ; 450, s. 133-149
-
Journal article (peer-reviewed)abstract
- Aims: To test the relative importance of topography versus soil chemistry in defining tree species-habitat associations in a terra firme Amazonian forest. Method: We evaluated habitat associations for 612 woody species using alternative habitat maps generated from topography and soil chemistry in the 25-ha Amacayacu Forest Dynamics Plot, Colombian Amazon. We assessed the ability of each habitat map to explain the community-level patterns of species-habitat associations using two methods of habitat randomization and different sample size thresholds (i.e., species’ abundance). Results: The greatest proportion of species-habitat associations arose from topographically-defined habitats (55% to 63%) compared to soil chemistry-defined (19% to 40%) or topography plus soil chemistry-defined habitats (18% to 42%). Results were robust to the method of habitat randomization and to sample size threshold. Conclusions: Our results demonstrate that certain environmental factors may be more influential than others in defining forest-level patterns of community assembly and that comparison of the ability of different environmental variables to explain habitat associations is a crucial step in testing hypotheses about the mechanisms underlying assembly. Our results point to topography-driven hydrological variation as a key factor structuring tree species distributions in what are commonly considered homogeneous Amazonian terra firme forests.
|
|
14. |
- Zuleta, Daniel, 1990, et al.
(author)
-
Individual tree damage dominates mortality risk factors across six tropical forests
- 2022
-
In: New Phytologist. - : Wiley. - 0028-646X .- 1469-8137. ; 233, s. 705-721
-
Journal article (peer-reviewed)abstract
- The relative importance of tree mortality risk factors remains unknown, especially in diverse tropical forests where species may vary widely in their responses to particular conditions. We present a new framework for quantifying the importance of mortality risk factors and apply it to compare 19 risks on 31 203 trees (1977 species) in 14 one-year periods in six tropical forests. We defined a condition as a risk factor for a species if it was associated with at least a doubling of mortality rate in univariate analyses. For each risk, we estimated prevalence (frequency), lethality (difference in mortality between trees with and without the risk) and impact (‘excess mortality’ associated with the risk, relative to stand-level mortality). The most impactful risk factors were light limitation and crown/trunk loss; the most prevalent were light limitation and small size; the most lethal were leaf damage and wounds. Modes of death (standing, broken and uprooted) had limited links with previous conditions and mortality risk factors. We provide the first ranking of importance of tree-level mortality risk factors in tropical forests. Future research should focus on the links between these risks, their climatic drivers and the physiological processes to enable mechanistic predictions of future tree mortality.
|
|
15. |
- Zuleta, Daniel, 1990, et al.
(author)
-
Interspecific and intraspecific variation of tree branch, leaf and stomatal traits in relation to topography in an aseasonal Amazon forest
- 2022
-
In: Functional Ecology. - : Wiley. - 0269-8463 .- 1365-2435. ; 36, s. 2955-2968
-
Journal article (peer-reviewed)abstract
- Tropical forest responses to variation in water availability, which are critical for understanding and predicting the effects of climate change, depend on trait variation among trees. We quantified interspecific and intraspecific variation in 18 branch, leaf and stomatal traits for 19–72 dominant tree species along a local topographic gradient in an aseasonal Amazon terra firme forest, and tested trait relationships with tree size, elevation, and species' topographic associations. We further tested whether correlation and coordination of traits vary among trees, among species and/or among trees within species. Intraspecific trait variation was substantial and exceeded interspecific variation in 10 of 18 traits. For leaf acquisition traits, intraspecific variation was mainly related to tree topographic elevation, while most of the variation in branch, leaf and stomatal traits was related to tree size. Interspecific variation showed no clear relationships with species' habitat association. Although trait correlations and coordinations were generally maintained among trees and among species, bivariate relationships varied among trees within species, across habitat association classes and across tree size classes. Our results demonstrate the magnitude and importance of intraspecific trait variation in tropical trees, especially as related to tree size. Furthermore, these results suggest that previous findings relating interspecific variation with topographic association in seasonal forests do not necessarily generalize to aseasonal forests. Read the free Plain Language Summary for this article on the Journal blog.
|
|
16. |
- Zuleta, Daniel, 1990, et al.
(author)
-
Local and regional determinants of vascular epiphyte mortality in the Andean mountains of Colombia
- 2016
-
In: Journal of Ecology. - : Wiley. - 0022-0477 .- 1365-2745. ; 104, s. 841-849
-
Journal article (peer-reviewed)abstract
- We present the first large-scale assessment of vascular epiphyte mortality in the neotropics. Our goals were to explore the primary types of vascular epiphyte death and to identify local and regional determinants of epiphyte mortality in natural forests located 60 to 2900m a.s.l. in the Colombian Andes. Based on two consecutive annual surveys, we followed the fate of 4247 epiphytes to estimate the epiphyte mortality rate on 116 host trees at nine sites. A logistic regression analysis for proportional data with a binomial distribution of the error was applied to determine the probability of epiphyte death in relation to local and regional explanatory variables. The overall epiphyte mortality rate was 7.5±1.1% year-1 (mean ± standard error). Non-mechanical factors, such as desiccation, accounted for a mortality rate of 1.9±0.3% year-1. Mechanical factors, such as falling branches, accounted for a mortality rate of 5.6±1.1% year-1. According to generalized linear modelling analyses, both local and regional factors played key roles in determining epiphyte mortality. The actual evapotranspiration (regional factor) and the mean epiphyte attachment height (local factor) were both consistently positively associated with the probability of epiphyte death. Additional variables identified as possible determinants of the epiphyte mortality were the temperature seasonality, annual temperature range, the height and number of branches of the tree and the abundance of large trees (DBH ≥10cm). Synthesis. The recorded high mortality rate indicates that natural epiphyte assemblages must be highly dynamic to avoid local extinction of species. Our study identifies actual evapotranspiration as an important driver of epiphyte mortality, and we highlight its importance in determining the fate tropical epiphyte communities may experience if evapotranspiration increases due to climate change. We hope our study addresses the paucity of research on non-tree growth forms, typically ignored in vegetation dynamics, and encourages their inclusion in future studies that investigate the function of tropical ecosystems. We present the first large-scale assessment of vascular epiphyte mortality in the neotropics. A recorded high mortality rate indicates that natural epiphyte assemblages must be highly dynamic to avoid local extinction of species. We hope this study will encourage the inclusion of non-tree growth forms in future studies that investigate the function of tropical ecosystems. Photography by Samuel Monsalve.
|
|