ABSTRACT
In Amazonian terra-firme non inundated forests, local floristic composition and species occurrence are explained by water availability as determined by topographic conditions. Topographic complexity can render these conditions quite variable across the landscape and the effects on plant ecological responses are difficult to document. We used a set of topographically defined hydrological metrics to evaluate community composition and single-species responses of four plant groups [pteridophytes (ferns and lycophytes), Melastomataceae, palms (Arecaceae) and Zingiberales] to topographic conditions in the middle Juruá River region, in western Brazilian Amazonia. The area spans two geological formations (Içá and Solimões) with contrasting topography. River terraces are also found along the main rivers in the area. Local topographic conditions were approximated by height above the nearest drainage (HAND), slope, and Strahler´s drainage order, all obtained from a SRTM digital elevation model (DEM). Data were analyzed using linear and generalized linear mixed models and regression trees. HAND was most successful in explaining floristic composition for all plant groups, except for Melastomataceae, and was more important in the hilly Içá formation than in the Solimões. Individual occurrences of 57% species were predicted by at least one of the topographic variables, suggesting a marked habitat specialization along topographic gradients. For these species, response models using SRTM-DEM-derived variables gave similar results than models using field-measured topography only. Our results suggest that topographical variables estimated from remote sensing can be used to predict local variation in the structure of plant communities in tropical forests.
KEYWORDS:
DEM; geological formations; HAND; vegetation mapping
RESUMO
Nas florestas de terra firme não inundáveis da Amazônia, a composição florística e a ocorrência de espécies podem ser explicadas pela disponibilidade hídrica relacionada com a topografia. Dada a complexidade topográfica, a disponibilidade de água pode ser bastante variável e seus efeitos na resposta das plantas, difícil de documentar. Neste estudo avaliamos as respostas individuais de espécie de quatro grupos de plantas [pteridófitas (samambaias e licófitas), Melastomataceae, palmeiras (Arecaceae) e Zingiberales] às condições topográficas na região do médio Rio Juruá, no oeste da Amazônia brasileira. A área abrange duas formações geológicas (Içá e Solimões) com topografias contrastantes. Terraços fluviais também são encontrados ao longo dos rios principais. As condições topográficas foram medidas usando a altura acima da drenagem mais próxima (HAND), declividade e ordem de drenagem de Strahler, todas obtidas a partir de um modelo digital de elevação SRTM-DEM. Os dados foram analisados usando modelos lineares generalizados mistos e árvores de regressão. HAND foi a principal variável explicativa da composição florística para todos os grupos de plantas, exceto Melastomataceae, tendo maior efeito na formação Içá do que na Solimões. Ocorrências individuais de 57% das espécies foram explicadas por pelo menos uma das variáveis, sugerindo uma especialização marcada de habitat ao longo de gradientes topográficos. Para essas espécies, modelos usando variáveis derivadas do SRTM-DEM deram resultados semelhantes aos modelos usando apenas a topografia medida em campo, o que indicam que variáveis topográficas derivadas do SRTM-DEM podem ser usadas para prever variações locais na estrutura de comunidades de plantas em florestas tropicais.
PALAVRAS-CHAVE:
DEM; formações geológicas; HAND; mapeamento de vegetação
INTRODUCTION
Answering general questions on plant species distribution in Amazonia has been hampered by the lack of detailed information on how the species niche is constrained by environmental conditions (Tuomisto et al. 2019Tuomisto, H.; Van doninck, J.; Ruokolainen, K.; Moulatlet, G.M.; Figueiredo, F.O.G.; Sirén, A.; et al. 2019. Discovering floristic and geoecological gradients across Amazonia. Journal of Biogeography, 46: 1734-1748.). In particular, soil nutrient concentration and water availability are important determinants of plant species distributions and floristic composition. Several studies show that occurrence and abundance of Amazonian plant species vary along gradients of soil nutrient concentration (Ca, Mg, K, Na) and available phosphorus (Tuomisto et al. 1998Tuomisto, H.; Poulsen, A.D.; Moran, R.C. 1998. Edaphic distribution of some species of the fern genus Adiantum in western Amazonia. Biotropica, 30: 392-399.; 2016; Costa et al. 2005Costa, F.R.C.; Magnusson, W.E.; Luizao, R.C. 2005. Mesoscale distribution patterns of Amazonian understorey herbs in relation to topography, soil and watersheds. Journal of Ecology, 93: 863-878.; Higgins et al. 2011Higgins, M.A.; Ruokolainen, K.; Tuomisto, H.; Llerena, N.; Cardenas, G.; Phillips, O.L.; et al. 2011. Geological control of floristic composition in Amazonian forests. Journal of Biogeography, 38: 2136-2149.; Zuquim et al. 2012Zuquim, G.; Tuomisto, H.; Costa, F.R.C.; Prado, J.; Magnusson, W.E.; Pimentel, T.; et al. 2012. Broad scale distribution of ferns and lycophytes along environmental gradients in central and northern Amazonia, Brazil. Biotropica, 44: 752-762.; Baldeck et al. 2016Baldeck, C.A.; Tupayachi, R.; Sinca, F.; Jaramillo, N.; Asner, G.P. 2016. Environmental drivers of tree community turnover in western Amazonian forests. Ecography, 39: 1089-1099.; Cámara-Leret et al. 2017Cámara-Leret, R.; Tuomisto, H.; Ruokolainen, K.; Balslev, H.; Munch Kristiansen, S. 2017. Modelling responses of western Amazonian palms to soil nutrients. Journal of Ecology, 105: 367-381.). Floristic composition also varies along local topographical gradients from hilltops to valley bottoms (Tuomisto et al. 1995; Tuomisto and Poulsen 2000Tuomisto, H.; Poulsen, A.D. 2000. Pteridophyte diversity and species composition in four Amazonian rain forests. Journal of Vegetation Science, 11: 383-396.; Poulsen et al. 2006Poulsen, A.D.; Tuomisto, H.; Balslev, H. 2006. Edaphic and floristic variation within a 1-ha plot of lowland Amazonian rain forest. Biotropica, 38: 468-478.) because topographic units (i.e., bottomlands, slopes and hilltops or plateaus) differ substantially in water availability (Nobre et al. 2011Nobre, A.D.; Cuartas, L.A.; Hodnett, M.; Rennó, C.D.; Rodrigues, G.; Silveira, A.; et al. 2011. Height Above the Nearest Drainage - a hydrologically relevant new terrain model. Journal of Hydrology, 404: 13-29.). As topography and soil nutrients can be correlated due to topographically controlled soil formation, plant species may respond to both, as they are exposed to geochemical sedimentary layers with different nutrient concentration at different topographical positions (Tuomisto and Ruokolainen 1994, Vormisto et al. 2004Vormisto, J.; Tuomisto, H.; Oksanen, J. 2004. Palm distribution patterns in Amazonian rainforests: What is the role of topographic variation? Journal of Vegetation Science, 15: 485-494.; Baldeck et al. 2013; Chauvel et al. 1987Chauvel, A.; Lucas, Y.; Boulet, R. 1987. On the genesis of the soil mantle of the region of Manaus, Central Amazonia, Brazil. Experientia, 43: 234-241.; Osher and Buol 1998Osher, L.J.; Buol, S.W. 1998. Relationship of soil properties to parent material and landscape position in eastern Madre de Dios, Peru. Geoderma, 83: 143-166.; Chadwick and Asner 2016Chadwick, K.D.; Asner, G.P. 2016. Tropical soil nutrient distributions determined by biotic and hillslope processes. Biogeochemistry, 127: 273-289.). Yet floristic changes in topographic gradients are less obvious than soil gradients due to the absence of a metric that allows cross-scale comparisons of topographic conditions.
In Amazonia large areas are seasonally inundated on the river floodplains, which makes them ecologically and floristically distinct from areas of terra firme (Campbell et al. 1986Campbell, D.G.; Daly, D.C.; Prance, G.T.; Maciel, U.N. 1986. Quantitative ecological inventory of terra firme and várzea tropical forest on the Rio Xingu, Brazilian Amazon. Brittonia, 38: 369-393.; Balslev et al. 1987Balslev, H.; Luteyn, J.; Øllgaard, B.; Holm-Nielsen, L.B. 1987. Composition and structure of adjacent unflooded and floodplain forest in Amazonian Ecuador. Opera Botanica, 92: 37-57.; Gentry 1988Gentry, A.H. 1988. Changes in plant community diversity and floristic composition on environmental and geographical gradients. Annals of the Missouri Botanical Garden, 75: 1-34.; Wittmann et al. 2013Wittmann, F.; Householder, E.; Piedade, M.T.F.; de Assis, R.L.; Schöngart, J.; Parolin, P.; et al. 2013. Habitat specifity, endemism and the neotropical distribution of Amazonian white-water floodplain trees. Ecography, 36: 690-707.). In terra firme forests, environmental gradients are often defined by topographic variation. Topography, in turn, is a proxy for several environmental conditions that are important determinants of species distribution (Vormisto et al. 2004Vormisto, J.; Tuomisto, H.; Oksanen, J. 2004. Palm distribution patterns in Amazonian rainforests: What is the role of topographic variation? Journal of Vegetation Science, 15: 485-494.; Moeslund et al. 2013Moeslund, J.E.; Arge, L.; Bøcher, P.K.; Dalgaard, T.; Odgaard, M.V.; Nygaard, B.; et al. 2013. Topographically controlled soil moisture is the primary driver of local vegetation patterns across a lowland region. Ecosphere, 4: 1-26.). Especially drainage has a clear relation to local topography, determining water availability as lower areas have higher water availability due to proximity to the water table and due to runoff of rainwater from higher areas (Fan et al. 2019Fan, Y.; Clark, M.; Lawrence, D.M.; Swenson, S.; Band, L.E.; Brantley, S.L.; et al. 2019. Hillslope hydrology in global change research and earth system modeling. Water Resources Research, 55: 1737-1772.). These characteristics make it possible to use topography as a practical surrogate of water availability on the terrain (Haitjema and Mitchell-Bruker 2005Haitjema, H.M.; Mitchell-Bruker, S. 2005. Are water tables a subdued replica of the topography? Ground Water, 43: 781-786.; Rennó et al. 2008Rennó, C.D.; Nobre, A.D.; Cuartas, L.A.; Soares, J.V.; Hodnett, M.G.; Tomasella, J.; et al. 2008. HAND, a new terrain descriptor using SRTM-DEM: Mapping terra-firme rainforest environments in Amazonia. Remote Sensing of Environment, 112: 3469-3481.; Costa et al. 2022).
Direct measurements of soil water availability or water table depth are time-consuming and very much affected by weather conditions during and just before measurement. In fact, few studies on plant species distribution in Amazonian terra firme have correlated such measurements with plant species distribution (e.g., Jirka et al. 2007Jirka, S.; McDonald, A.J.; Johnson, M.S.; Feldpausch, T.R.; Couto, E.G.; Riha, S.J. 2007. Relationships between soil hydrology and forest structure and composition in the southern Brazilian Amazon. Journal of Vegetation Science, 18: 183-194.; Guimaraes et al. 2021Guimaraes, A.F.; de Souza, C.R.; Rosa, C.; Paulo dos Santos, J.; Teixeira, L.A.F.; Zanzini, L.P.; et al. 2021. Small-scale environmental variations drive vegetation structure and diversity in Amazon riverine forests. Flora, 283: 151916.). Instead, topographic variation has been used as a surrogate for local hydrological conditions when addressing changes in local floristic composition, either by measuring topographic variation in the field (Tuomisto et al. 1995Tuomisto, H.; Ruokolainen, K.; Kalliola, R.; Linna, A.; Danjoy, W.; Rodriguez, Z. 1995. Dissecting Amazonian biodiversity. Science, 269: 63-66.; Tuomisto and Poulsen 2000Tuomisto, H.; Poulsen, A.D. 2000. Pteridophyte diversity and species composition in four Amazonian rain forests. Journal of Vegetation Science, 11: 383-396.; Drucker et al. 2008Drucker, D.P.; Costa, F.R.C.; Magnusson, W.E. 2008. How wide is the riparian zone of small streams in tropical forests? A test with terrestrial herbs. Journal of Tropical Ecology, 24: 65-74.; Costa et al. 2009Costa, F.R.C.; Guillaumet, J.-L.; Lima, A.P.; Pereira, O.S. 2009. Gradients within gradients: The mesoscale distribution patterns of palms in a central Amazonian forest. Journal of Vegetation Science, 20: 69-78.), or by making a simple (and subjective) classification of topographic units (e.g. plateaus, slopes and valleys) into hydromorphic environments (Kahn and de Castro 1985Kahn, F.; de Castro, A. 1985. The palm community in a forest of central Amazonia, Brazil. Biotropica, 17: 210-216.; Svenning 1999Svenning, J.-C. 1999. Recruitment of tall arborescent palms in the Yasuní National Park, Amazonian Ecuador: are large treefall gaps important? Journal of Tropical Ecology, 15: 355-366.; Drucker et al. 2008). In the absence of field measured topography, topographic metrics derived from digital elevation models (DEMs) have been applied as determinants of floristic composition. Out of these metrics, the algorithm HAND (height above the nearest drainage, Rennó et al. 2008Rennó, C.D.; Nobre, A.D.; Cuartas, L.A.; Soares, J.V.; Hodnett, M.G.; Tomasella, J.; et al. 2008. HAND, a new terrain descriptor using SRTM-DEM: Mapping terra-firme rainforest environments in Amazonia. Remote Sensing of Environment, 112: 3469-3481.) has stood out as a proxy for topographically defined hydrological conditions in terrains with relatively homogeneous soil conditions (Moulatlet et al. 2014Moulatlet, G.M.; Costa, F.R.C.; Rennó, C.D.; Emilio, T.; Schietti, J. 2014. Local hydrological conditions explain floristic composition in lowland Amazonian forests. Biotropica, 46: 395-403.; Schietti et al. 2014Schietti, J.; Emilio, T.; Rennó, C.D.; Drucker, D.P.; Costa, F.R.C.; Nogueira, A.; et al. 2014. Vertical distance from drainage drives floristic composition changes in an Amazonian rainforest. Plant Ecology & Diversity, 7: 241-253.), but in terrains with highly heterogeneous soil conditions, its applicability has not been evaluated yet.
Remote sensing data may help to fill the gap in environmental information needed to map species distribution at different scales (Rocchini et al. 2016Rocchini, D.; Boyd, D.S.; Féret, J.-B.; Foody, G.M.; He, K.S.; Lausch, A.; et al. 2016. Satellite remote sensing to monitor species diversity: potential and pitfalls. Remote Sensing in Ecology and Conservation, 2: 25-36.). However, the association of environmental processes with species distribution from local to regional scales needs remote sensing data coupled with data on species occurrences that covers a large part of the environmental heterogeneity (Tuomisto et al. 2019Tuomisto, H.; Van doninck, J.; Ruokolainen, K.; Moulatlet, G.M.; Figueiredo, F.O.G.; Sirén, A.; et al. 2019. Discovering floristic and geoecological gradients across Amazonia. Journal of Biogeography, 46: 1734-1748.). We used an extensive dataset of field measured topography in the western Brazilian Amazon to assess the ability of HAND and two other DEM-derived topographic metrics (i.e., the estimated topographic conditions) as proxies of local hydrological conditions to explain floristic composition and plant species distribution across two geological formations. We focused on four plant groups (pteridophytes, palms, Melastomataceae and Zingiberales) and addressed community and individual species responses to both field and DEM-derived topography. We specifically asked: i) whether species distributions of the different plant groups can be explained by the estimated topographic conditions; ii) what is the response magnitude and direction (i.e., positive, or negative) of individual species to the estimated topographic conditions; and iii) to which degree the species responses can be compared to local topographic measurements, i.e., if DEM-derived variables can replace topographic field measurements as useful metrics to be associated with floristics.
MATERIAL AND METHODS
Study Area
The study was conducted in western Amazonia in Brazil along the rivers Juruá and Tarauacá, with field sampling extending across 500 km (Figure 1). The study area is covered by continuous rainforest at an elevation range of 65-200 m.a.s.l. Mean annual rainfall is 2200 mm and annual mean temperature is 27 °C, with temperatures as low as 15 ºC occurring during annual cold spells (Karger et al. 2017Karger, D.N.; Conrad, O.; Böhner, J.; Kawohl, T.; Kreft, H.; Soria-Auza, R.W.; et al. 2017. Climatologies at high resolution for the earth’s land surface areas. Scientific Data, 4: 170122.). The area is in a complex geological setting. The Içá and Solimões formations underlie most of the terrain, and river terraces are present along both rivers (Table 1). The Içá Formation consists of relatively nutrient poor and loamy to sandy sediments deposited during the Pliocene to Pleistocene, and the terrain is typically steep and hilly (Hoorn and Wesselingh 2011Hoorn, C.; Wesselingh, F. 2011. Amazonia, Landscape and Species Evolution: A Look into the Past. John Wiley & Sons, Chichester, 883p.). The Solimões Formation consists of relatively nutrient rich clay sediments deposited under semi-marine or lacustrine conditions during the Miocene (Hoorn and Wesselingh 2011), and the topography is generally flat or undulating. Alluvial terraces extend along the main rivers that consist of sediments of varying textures in abandoned floodplains due to lateral migration of the river channel. The differences in soil nutrient concentration of these geological formations have already been shown to explain the main gradients in plant species composition (Higgins et al. 2011Higgins, M.A.; Ruokolainen, K.; Tuomisto, H.; Llerena, N.; Cardenas, G.; Phillips, O.L.; et al. 2011. Geological control of floristic composition in Amazonian forests. Journal of Biogeography, 38: 2136-2149.; Tuomisto et al. 2016Tuomisto, H.; Moulatlet, G.M.; Balslev, H.; Emilio, T.; Figueiredo, F.O.G.; Pedersen, D.; et al. 2016. A compositional turnover zone of biogeographical magnitude within lowland Amazonia. Journal of Biogeography, 43: 2400-2411.).
Location of the study area in western Amazonia (Brazil). Sampling transects (white dots) were located along the rivers Juruá and Tarauacá. The dark line shows the limit between the Içá and Solimões formations, as determined by Higgins et al. (2011Higgins, M.A.; Ruokolainen, K.; Tuomisto, H.; Llerena, N.; Cardenas, G.; Phillips, O.L.; et al. 2011. Geological control of floristic composition in Amazonian forests. Journal of Biogeography, 38: 2136-2149.). Background shading represents elevation differences as indicated by the SRTM DEM. Lighter areas are of higher elevation than darker areas.
Environmental variables measured along transects set in each of three geological settings along the middle Juruá River (western Brazilian Amazon). Mean (minimum - maximum) values are shown.
Sampling design
A total of 71 line transects were established between March and June of 2012. The transects were placed in sites where Landsat image interpretation suggested that primary terra-firme forest was sufficiently close to the river to be reached by foot and inventoried within one day of trekking. The exact location and compass bearing of transects were chosen to include as much as possible of the local topographical heterogeneity. Each transect was 500 m long and consisted of 20 subunits of 25 m x 5 m. The topography of each transect was measured using a clinometer at 20-25-m intervals and whenever there was a significant change in slope. When transects were so flat that the measurement error was estimated to be larger than the actual slope of the terrain, topography was not measured, and the slope was assumed to equal zero. Transect coordinates were obtained in the field with a hand-held GPS.
Floristic data
We use the same floristic data as Tuomisto et al. (2016Tuomisto, H.; Moulatlet, G.M.; Balslev, H.; Emilio, T.; Figueiredo, F.O.G.; Pedersen, D.; et al. 2016. A compositional turnover zone of biogeographical magnitude within lowland Amazonia. Journal of Biogeography, 43: 2400-2411.). Floristic inventories were made at the 25-m subunit resolution within each transect. For pteridophytes, we registered all individuals with at least one green leaf (leafy stem in the case of lycophytes) longer than 10 cm, including epiphytes, hemiepiphytes and climbers if they had such leaves less than two meters above ground. All Melastomataceae individuals with post-cotyledon leaves were included. All Zingiberales individuals with a minimum height of 5 cm were recorded and, in the case of clonal species, bunches of leaves separated by at least 20 cm were considered as separate individuals. All palm individuals higher than 5 cm were included, but palm seedlings that could not be identified to species level were excluded. Each ramet in a clonal or colonial species was counted as an individual. Representative voucher specimens were deposited in the herbaria of the University of Turku, Finland (pteridophytes and Melastomataceae), University of Aarhus, Denmark (palms), Instituto Nacional de Pesquisas da Amazônia, Brazil (palms and Zingiberales) and Instituto de Botanica - Universidade de São Paulo, Brazil (pteridophytes and Melastomataceae). Pteridophytes and Melastomataceae were inventoried in all 71 transects, palms in 40 and Zingiberales in 39 (Supplementary Material, Table S1).
Topographic data
To place all transects on the same elevational and topographical framework, we used the digital elevation model (DEM) of the Shuttle Radar Topographic Mission (SRTM) with resolution of 1 arc sec (~30 m at the equator) to obtain elevation data for each 25-m subunit. Before extracting the SRTM elevation values for each subunit, we compared the SRTM-derived topographic profile of each transect with the corresponding field-measured true profile. This was done to account for possible errors associated with the measurements of transect coordinates. If there was a mismatch, we adjusted transect georeferencing to obtain the best possible match of the SRTM profile with the field-measured topography using an algorithm that optimized the search for initial points with similar topographic profiles. The initial SRTM profiles were based on the field-measured coordinates at the start and end points of each transect, as transects were assumed to be straight lines without major deviations. The coordinate correction procedure iteratively recalculated the start and end coordinates of each transect, such that the sum of the squared differences between field-measured and SRTM-derived topographic profiles were minimized. This was done with the constraint that the new coordinates could not be displaced by more than two SRTM pixels (approximately 60 m) from the original coordinates. In the next step, the new coordinates were used to assign coordinates for each 25-m subunit and to extract SRTM elevation values by bilinear interpolation. With this procedure, we had SRTM-derived topographic profiles that matched the topographies measured in the field. For five of the 71 transects (transects 780, 794, 806, 807 and 808) our algorithm failed to match with the topographic profiles within the range of 60 m and, therefore, these transects were not included in the final analysis. The R code for the algorithm can be found at (https://github.com/gamamo/Jurua-Hydro).
Based on the adjusted transect georeferencing, we calculated a single HAND value for each transect subunit. HAND is a topographic measure that equals the vertical distance (in meters) between a point of interest and the nearest watercourse (or point where the water table is assumed to be near the surface). Low HAND values correspond to moister sites close to the water table, and high HAND values to drier sites further above the water table. Since HAND is independent of absolute elevation, its values are comparable over large areas. We calculated HAND using a remote sensing-based algorithm developed by Rennó et al. (2008Rennó, C.D.; Nobre, A.D.; Cuartas, L.A.; Soares, J.V.; Hodnett, M.G.; Tomasella, J.; et al. 2008. HAND, a new terrain descriptor using SRTM-DEM: Mapping terra-firme rainforest environments in Amazonia. Remote Sensing of Environment, 112: 3469-3481.) which combines the STRM-DEM with a drainage network. We set as drainage network for each geological formation the one produced by Banon (2013Banon, L. 2013. Árvores de decisão aplicadas à extração automática de redes de dranagem. Master´s dissertation, Instituto Nacional de Pesquisas Espaciais, Brazil. 87p. (http://urlib.net/8JMKD3MGP7W/3FB9RFE).
http://urlib.net/8JMKD3MGP7W/3FB9RFE...
), which uses a decision tree based on several topographic attributes extracted from the SRTM-DEM, including accumulated area, slope and curvatures, to determine the best drainage network for a given area.
To complement the HAND estimates, we derived two other topographic variables related to the local hydrological conditions, namely slope and stream order of the nearest drainage (according to Strahler 1957Strahler, A.N. 1957. Quantitative analysis of watershed geomorphology. Eos, Transactions American Geophysical Union, 38: 913-920.) from the SRTM-DEM. Slope is related to the overland and subsurface flow of water, and therefore it affects potential soil moisture, which is generally higher on gentler slopes. Strahler’s stream order is an estimate of the size of a water channel based on the same drainage network that was generated to calculate the HAND. In Amazonian forests, water table depth varies not only along topographic profiles, but also according to the distance to rivers of different drainage orders. Low values of Strahler`s drainage order are found often near the headwaters and indicate intermittent waterflow dominated by ground water streamflow. Strahler`s drainage order increases downstream and larger values indicate more continuous waterflow, the magnitude of which depends on the interplay between infiltration from large rivers and ground water depth (Miguez-Macho and Fan 2012Miguez-Macho, G.; Fan, Y. 2012. The role of groundwater in the Amazon water cycle: 1. Influence on seasonal streamflow, flooding and wetlands. Journal of Geophysical Research. Atmospheres, 117: 1-30.).
Data Analysis
The floristic responses (community-level analysis) to the topographic predictors were assessed using regression models developed for each of the four plant groups separately. To use community data as the response variable, we reduced the dimensionality of the data through NMDS (non-metric multidimensional scaling) to two dimensions. A floristic dissimilarity matrix for each NMDS was calculated using the extended Sørensen dissimilarity. This is based on the classical Sørensen similarity index, but it uses intermediate plots as steppingstones to calculate dissimilarities between plots that do not share species, which facilitates relating long compositional gradients to environmental variables (De’ath 1999De’ath, G. 1999. Extended dissimilarity: a method of robust estimation of ecological distances from high beta diversity data. Plant Ecology, 144: 191-199.).
We modeled the relationship between the NMDS axes and the set of topographic variables using linear mixed models (LMM) with gaussian distribution. Each of the two NMDS axes for each of the four plant groups was used separately as a response variable. We built models that included as explanatory variables HAND, drainage order and slope, and models that included only field measured topography as explanatory variable. To assess how the set of topographic variables explained the distribution of each individual species, we used generalized linear mixed models (GLMM) with binomial distribution (presence-absence data). Species occurrences were the response variables, and the quadratic forms of HAND, slope and Strahler’s drainage order were the independent terms. We also built GLMMs that included only field measured topography as explanatory variable. For each individual species modelling, only transects within the same geological formation where a species had been registered were included in the species model to avoid overpredictions. We restricted our analysis to those species considered as frequent in the dataset, with at least 20 occurrences. Transect identity and transect affiliation to geological formation were set as random terms in all LMMs. The first aimed to account for the autocorrelation of sub-samples within transects; the second, to account for the marked differences in the environmental conditions of each geological formation. In GLMMs, only transect identity was set as a random term because the modelled species were rare in the three geological formations, so geological differences were not a relevant factor in the models.
In both community (LMM) and species-level models (GLMM), all variables were standardized by subtracting the mean and dividing by the standard deviation prior to the analyses to compare the effect sizes based on model coefficients. For the models that contained HAND, slope and Strahler’s order as variables, we applied a model selection procedure that started with a full model and then, models with different combination of the variables were compared using the dredge function of the R package MuMIn (Barton 2018Barton, K. 2018. Mu-MIn: Multi-model inference. R package. (http://r-forge. r-project. org/projects/mumin/).
http://r-forge. r-project. org/projects/...
). The best model was selected based on the Akaike Information Criteria (AIC). When more than one model had AIC < 2, the best model was determined as an average of their estimated standardized coefficients (Anderson and Burnham 2004Anderson, D.; Burnham, K. 2004. Model Selection and Multi-Model Inference . v. 63. 2nd ed. Springer-Verlag, New York, 488p.), using the model.avg function of the R package MuMIn (Barton 2018). Model evaluation was done by comparing the model’s AIC with a null model assuming no explanatory variables. Then, we evaluated the effect of the different variables in explaining either community composition (NMDS axes) or individual species occurrence in each model by comparing the magnitude and significance of the standardized beta coefficients of the final models. Coefficient values closer to 1 or -1 indicate stronger effect, whereas coefficient values closer to 0 indicate non-significant effects. The likelihood-ratio based pseudo- R2 was reported individually for each model, with a maximum value of 1. LMM and GLMM using field measured topography as explanatory variable were also evaluated by comparing the model’s AIC with a null model assuming no explanatory variable. All analyses were performed in the R statistical software (R Core Team 2020R Core Team. 2020. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. (http://www.R-project.org/).
http://www.R-project.org/...
).
To evaluate the hierarchical importance of DEM-derived topographic variables in structuring plant communities in each geological formation, we used distance-based multivariate regression trees (MRT) (De’Ath 2002). We built MRT separately for each species group in each of the three geological surfaces (Solimões Formation, Içá Formation, and fluvial terrace). The two NMDS axes were set as dependent variables, and HAND, slope, and Strahler´s order as independent variables. MRTs were obtained by repeatedly splitting the observation units into two clusters that were determined by a break point in an environmental variable. Each split aimed to minimize the floristic differences between subunits in the same cluster, and at each level in the MRT, the variable that gives the minimum within-group sum of squared distance to the group mean was selected. This process was repeated with each sub-group of the previous step until each observation unit formed its own cluster.
We calculated dissimilarities using the Bray-Curtis index for occurrence data including all the species of each taxonomic group using the R package vegan (Oksanen et al. 2018). We used cross-validation to select the MRT with the smallest error (De’Ath 2002). The explained variation of the MRT is given by the residual error of the overall tree. MRTs were calculated using the functions of the R package mvpart (Therneau and Atkinson 2013).
RESULTS
Floristic composition
HAND was relevant to explain community composition of all plant groups except Melastomataceae (Figure 2a; Table 2). The effect was stronger for palms (b = -0.12) than for the other groups. Slope was only significant for pteridophytes (b = 0.04), and Strahler’s drainage order was only significant for palms (b = 0.05). When NMDS axis 2 was the response variable (Figure 2b), HAND was still a significant variable for the same plant groups and had a strong effect for palms (b = -0.12). Strahler’s drainage order was significant for palms (b = 0.10) and Zingiberales (b = -0.03). Slope had a significant effect for palms (b = -0.06) and Zingiberales (b = -0.03).
Output of the LMM for plant community data from the middle Juruá River (western Brazilian Amazon). Standardized slope coefficients and the standard error (horizontal bars) of the LMMs are shown for all variables used to model species occurrence with the NMDS axis 1 (A) and 2 (B). Coefficients different from zero are associated with significant effects. SDO = Strahler’s drainage order.
Summary of linear mixed models for plant communities along the middle Juruá River (western Brazilian Amazon) showing the results of model selection of the response of community composition (MDS1 and MDS2) to DEM-derived topographic metrics [HAND, slope and Strahler’s drainage order (SDO)] and to field measured topography (Topo). The fit of the best models, alongside the null model assuming no explanatory variables (intercept-only), are shown for every response variable. Null models are shown separately for models with DEM-derived variables (Null) and for models with field measured topography (Null Topo). Response variables with null values are those for which the null model AIC was lower than for any model that included the variable.
Individual species responses
Individual species responses to DEM-derived topographic variables were variable (Figure 3; Supplementary Material, Table S2). Pteridophytes (81%) and palms (80%) had the highest proportion of species that responded to the variation in HAND. Pteridophyte and Melastomataceae species responded mostly negatively, i.e., the probability of species occurrence decreased with increasing HAND values (Figure 3). For palm species, 55% had a positive response, i.e., species occurrence probability increased with increasing HAND values. Slope and Strahler’s drainage order were significant in several models for all plant groups (Figure 3; Supplementary Material, Table S2). Model fit, as measured by the pseudo-R2, varied from 0.05 - 0.45 in pteridophytes, 0.05 - 0.63 in Zingiberales, 0.05 - 0.45 in palms and 0.01 - 0.39 in Melastomataceae (Supplementary Material, Table S2). Out of the species restricted to a single geological formation, only one (Calathea neblinense - Zingiberales) had a positive relationship with HAND.
Output of GLMM for plant species occurrence along the middle Juruá River (western Brazilian Amazon). Dot size indicates the slope coefficient of the GLMM for all variables used to model species occurrence. Coefficients with higher variation are associated with stronger effects, either negative (decreasing species occurrence probability; red circles) or positive (increasing species occurrence probability; blue circles) of each variable. The tile plot to the right of each plot indicates in which geological formation each species was recorded (IF = Içá Formation; RT = river terraces; SF = Solimões Formation). SDO = Strahler’s drainage order. The numbers 1 and 2 after each variable’s name refer to the linear and squared terms, respectively. Species full names can be found in Table S2. This figure is in color in the electronic version.
Congruence with field measured topography
Overall, LMM and GLMM with field-measured topography as the sole explanatory variable had similar AIC values than models with DEM-derived topographic variables (Table 3, Supplementary Material, Figures S1, S2; Table S2). At the species level, the difference between the AICs, as assessed with an ANOVA, was not significant for any plant group (pteridophytes: p = 0.985, Zingiberales: p = 0.958, palms: p = 0.964, Melastomataceae: p = 0.927) (Supplementary Material, Figure S3).
Regression trees
Regression trees (Figure 4) showed that, in the river terraces, the main floristic differences for each species group were explained by drainage (pteridophytes, palms and Zingiberales, MRT residual errors = 0.87, 0.78 and 0.62, respectively) and slope (Melastomataceae, MRT residual error = 0.85). HAND was the variable that caused the main division of the MRT for the Içá Formation for all plant groups (MRT residual error: pteridophytes = 0.84, Zingiberales = 0.93, palms = 0.66, Melastomataceae = 0.94). In the Solimões Formation, the main division for all plant groups was associated with HAND (pteridophytes and Melastomataceae, MRT residual error = 0.88 and 0.96, respectively) and drainage (Zingiberales and palms, MRT residual error = 0.82 and 0.51, respectively).
Distance-based multivariate regression trees showing clusters of transects established along 500 km of the middle Juruá River (western Brazilian Amazon). At each split of the tree, the dataset was divided into two according to the variables with higher importance value (upper nodes) in explaining differences in species composition. The value of each variable used to separate communities is shown. Branches could be subdivided several times. The number of subunits (n) that share similar species composition based on the Bray-Curtis similarity index is shown for each branch. Drainage = Strahler’s drainage order.
DISCUSSION
HAND explained species composition of all plant groups, suggesting that this variable detected the influence of hydrological conditions as determined by topographic gradients. The underlying mechanisms to explain species distribution along topographic gradients are possibly related to species tolerance to desiccation or waterlogging. Desiccation is often associated with the distance to the water table on hilltops, while waterlogging is more common in areas close to the water table in local topographic bottomlands (Lopez and Kursar 1999Lopez, O.R.; Kursar, T.A. 1999. Flood tolerance of four tropical tree species. Tree Physiology, 19: 925-932.; Parent et al. 2008Parent, C.; Capelli, N.; Berger, A.; Crèvecoeur, M.; Dat, J.F. 2008. An overview of plant responses to soil waterlogging. Plant Stress, 2: 20-27.; Oliveira et al. 2018Oliveira, R.S.; Costa, F.R.C.; van Baalen, E.; de Jonge, A.; Bittencourt, P.R.; Almanza, Y.; et al. 2018. Embolism resistance drives the distribution of Amazonian rainforest tree species along hydro-topographic gradients. New Phytologist, 221:1457-1465.; Fontes et al. 2020Fontes, C.G.; Fine, P.V.A.; Wittmann, F.; Bittencourt, P.R.L.; Piedade, M.T.F.; Higuchi, N.; et al. 2020. Convergent evolution of tree hydraulic traits in Amazonian habitats: implications for community assemblage and vulnerability to drought. New Phytologist, 228: 106-120.; Garcia et al. 2022Garcia, M.N.; Hu, J.; Domingues, T.F.; Groenendijk, P.; Oliveira, R.S.; Costa, F.R.C. 2022. Local hydrological gradients structure high intraspecific variability in plant hydraulic traits in two dominant central Amazonian tree species. Journal of Experimental Botany, 73: 939-952.).
Slope was also a significant variable to explain the floristic composition of pteridophytes (NMDS 1) and palms and Zingiberales (NMDS 2). The response of these plant groups to slope might reflect the strong slope gradient of the study area. Along topographic profiles, the slope angle indicates the capacity of the soil to retain moisture that infiltrates from precipitation (Rennó et al. 2008Rennó, C.D.; Nobre, A.D.; Cuartas, L.A.; Soares, J.V.; Hodnett, M.G.; Tomasella, J.; et al. 2008. HAND, a new terrain descriptor using SRTM-DEM: Mapping terra-firme rainforest environments in Amazonia. Remote Sensing of Environment, 112: 3469-3481.). On steep slopes, high lateral and sub superficial runoff diminish water infiltration, which is a limiting factor to plants that do not have deeper roots to anchor and access moisture from deeper soil layers (Fan et al. 2019Fan, Y.; Clark, M.; Lawrence, D.M.; Swenson, S.; Band, L.E.; Brantley, S.L.; et al. 2019. Hillslope hydrology in global change research and earth system modeling. Water Resources Research, 55: 1737-1772.). On the other hand, steep slopes do not accumulate as much litter as flat areas, so the lack of litter facilitates the establishment of seedlings of plants that do not have large nutrient reserves, such as the pteridophytes (Rodrigues and Costa 2012Rodrigues, F.R. de O.; Costa, F.R.C. 2012. Litter as a filter of emergence for herbaceous seedlings and sporophytes in central Amazonia. Journal of Tropical Ecology, 28: 445-452.).
We found a significant effect of Strahler´s order in explaining floristic composition for palms (NMDS 1 and NMDS 2) and for Zingiberales (NMDS 2) in models where HAND also had a significant effect, suggesting that species responses to the modelled hydrological conditions were also related to the complementary effect of HAND in the gradient of drainage order formed along a watershed. Close to the watershed exit, bottomlands get inundated by rainwater because water infiltration is low when the soils are saturated due to the inundations of the large rivers (Miguez-Macho and Fan 2012Miguez-Macho, G.; Fan, Y. 2012. The role of groundwater in the Amazon water cycle: 1. Influence on seasonal streamflow, flooding and wetlands. Journal of Geophysical Research. Atmospheres, 117: 1-30.). At the headwaters, the influence of low Strahler’s order rivers on the soil saturation is dependent on terrain slope. While flat terrains can get frequently moist, in bottomlands of the areas with more pronounced slopes, the soil does not stay saturated after a rainfall because the small rivers have relatively low influence on soil saturation (Fan et al. 2013Fan, Y.; Li, H.; Miguez-Macho, G. 2013. Global patterns of groundwater table depth. Science, 339: 940-943.).
The occurrence of 57 % of the species in all plant groups was predicted using the DEM-derived variables. The models met some of our expectations. Species such as Mickelia nicotianifolia, Campyloneurum fuscosquamatum (pteridophytes), Astrocaryum ulei (palm) and Clidemia epiphytica (Melastomataceae) are commonly found in bottomlands and were correctly modelled to be negatively related to HAND, while Oenocarpus bataua (palm), commonly found on hilltops in western Amazonia, were correctly modelled to be positively related to HAND. Also, Calathea altisssima (ginger) has been described as common species on steep slopes and our models predicted its occurrence correctly. Habitat specialization related to topographic position was reported before for many palms (Svenning 1999Svenning, J.-C. 1999. Recruitment of tall arborescent palms in the Yasuní National Park, Amazonian Ecuador: are large treefall gaps important? Journal of Tropical Ecology, 15: 355-366.; Vormisto et al. 2004Vormisto, J.; Tuomisto, H.; Oksanen, J. 2004. Palm distribution patterns in Amazonian rainforests: What is the role of topographic variation? Journal of Vegetation Science, 15: 485-494.), pteridophytes (Tuomisto and Ruokolainen 1994Tuomisto, H.; Ruokolainen, K. 1994. Distribution of Pteridophyta and Melastomataceae along an edaphic gradient in an Amazonian rain forest. Journal of Vegetation Science, 5: 25-34.; Tuomisto et al. 1995; 1998; Zuquim et al. 2012Zuquim, G.; Tuomisto, H.; Costa, F.R.C.; Prado, J.; Magnusson, W.E.; Pimentel, T.; et al. 2012. Broad scale distribution of ferns and lycophytes along environmental gradients in central and northern Amazonia, Brazil. Biotropica, 44: 752-762.), Melastomataceae (Tuomisto and Ruokolainen 1994) and Zingiberales (Costa et al. 2005Costa, F.R.C.; Magnusson, W.E.; Luizao, R.C. 2005. Mesoscale distribution patterns of Amazonian understorey herbs in relation to topography, soil and watersheds. Journal of Ecology, 93: 863-878.). These previous studies have used field-measured topographical variables, which accurately reflect slope and relative topographic position, but do not contain information about distance to water table and cannot be correlated across study areas. Our models using DEM-derived topography metrics were similar to models using field measured topography as explanatory variables at both community and species level, which is consistent with the idea that specialization to local topography might be related to species tolerance to water availability. However, as soil conditions were not evaluated in our study, the usefulness of DEM-variables as surrogates of field measured topography still requires further investigation.
Floristic patterns differed between geological formations. In the Içá Formation, HAND was the most important predictor of floristic composition for all the plant groups. This result may have been an effect of the Içá Formation being hillier, with a wider range of environmental conditions from lower to higher topographic positions relative to the Solimões Formation or the river terraces. Water availability likely changes for plants growing in different parts of this topographic gradient. Moreover, other soil conditions related to water availability may explain species distribution. Soils of the Içá Formation and river terraces in high topographic positions have coarser texture (Tuomisto et al. 2016Tuomisto, H.; Moulatlet, G.M.; Balslev, H.; Emilio, T.; Figueiredo, F.O.G.; Pedersen, D.; et al. 2016. A compositional turnover zone of biogeographical magnitude within lowland Amazonia. Journal of Biogeography, 43: 2400-2411.) and may dry out faster than those in lower topographic levels. In the Solimões Formation and river terraces, finer sediments prevail, and the water table is superficial, so that water availability is not a limiting factor for plants in dry periods. This would explain why drainage order was so important to explain floristic composition in the river terraces.
Applicability of DEM-derived variables
Our DEM-derived hydrological variables cannot be causally related to the physiological requirements of plants. Actual measurements of water table depth and soil moisture during sufficiently long periods of time would allow to characterize the role of hydrological conditions more accurately across topographic gradients in the different geological formations. However, long term field measurements across large areas in remote parts of Amazonia are logistically impracticable and of limited use for environmental monitoring (Frappart et al. 2019Frappart, F.; Papa, F.; Güntner, A.; Tomasella, J.; Pfeffer, J.; Ramillien, G.; et al. 2019. The spatio-temporal variability of groundwater storage in the Amazon River Basin. Advances in Water Resources, 124: 41-52.). Algorithms like HAND make it possible to remove the effect of regional topography and to retain only the local topography, so local hydrological conditions as estimated by HAND can be an efficient surrogate measure to describe hillslope hydrology (Fan et al. 2019Fan, Y.; Clark, M.; Lawrence, D.M.; Swenson, S.; Band, L.E.; Brantley, S.L.; et al. 2019. Hillslope hydrology in global change research and earth system modeling. Water Resources Research, 55: 1737-1772.) that can then be associated with species distribution. As the HAND eliminates the effect of absolute elevation values, the obtained results can be assumed to be quantitatively comparable to other parts of Amazonia.
The HAND algorithm has been used for ecological applications previously (Banon et al. 2019Banon, G.P.R.; Banon, G.J.F.; Villamarín, F.; Arraut, E.M.; Moulatlet, G.M.; Rennó, C.D.; et al. 2019. Predicting suitable nesting sites for the Black caiman (Melanosuchus niger Spix 1825) in the Central Amazon basin. Neotropical Biodiversity, 5: 47-59.; Zuquim et al. 2021Zuquim, G.; Tuomisto, H.; Chaves, P.P.; Emilio, T.; Moulatlet, G.M.; Ruokolainen, K.; et al. 2021. Revealing floristic variation and map uncertainties for different plant groups in western Amazonia. Journal of Vegetation Science, 32: e13081.). Its applicability, however, is still limited by the need of validation points. HAND needs validation because the SRTM does not always accurately reflect topography due to complex interaction with the forest canopy, in a way that highly dissected valleys may remain unmapped (Valeriano et al. 2006Valeriano, M.M.; Kuplich, T.M.; Storino, M.; Amaral, B.D.; Mendes Jr., J.N.; Lima, D.J. 2006. Modeling small watersheds in Brazilian Amazonia with shuttle radar topographic mission-90 m data. Computers & Geosciences, 32: 1169-1181.). Our results show the usefulness of HAND and other DEM-derived metrics as predictors of floristic composition, but edaphic conditions need to be properly controlled to allow further interpretations. Thus, an assessment of the correlation between HAND and other environmental variables would help to disentangle its effect as a proxy for local hydrological conditions in Amazonian landscapes.
CONCLUSIONS
Our study showed that DEM-derived topographic variables can be useful to detect species preferences for local hydrological conditions defined by topographic gradients. The selected DEM-derived topographic variables were able to detect changes in floristic composition and individual species occurrences. The integration of topographic modelling of hydrological conditions in species modelling frameworks represents a significant advance in the mapping of Amazonian biodiversity, as knowledge on species tolerance to water- and soil-related environmental gradients is also critical for predicting how the viability and distribution of species might be affected by environmental changes.
ACKNOWLEDGMENTS
We are thankful to the people who made our expedition along the Juruá River possible and fruitful: the boat crew (José Rodrigues da Silva Filho, Radisvom Teixeira de Lima and Perla Maciel da Silva) and the local field assistants (especially José Francisco Gomes da Cruz). Research permits were granted by the Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) (process # 001693/2011-5), Centro Estadual de Unidades de Conservação (CEUC)/Secretaria do Meio Ambiente e Desenvolvimento Sustentável do Estado do Amazonas (SDS) (authorization # 047/2011) and Sistema de Autorização e Informação em Biodiversidade (SISBIO) (license # 27998-3). This work was funded by the University of Turku Graduate School (grant to GMM), the Academy of Finland (grants # 139959 and 273737 to HT), Danish National Research Council (grant # 9040-00136B to HB), European Community (FP7-ENB-2007-1; EU Grant Agreement 212631 to HB), CNPq (fellowship to TE), Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP) (project # 2016/13462-0, grant # 2016/18925-9 to LB), and Coordenação de Aperfeiçoamento de Pessoal de Nivel Superior (CAPES) (finance code 001). GMM acknowledges the postdoctoral grant from the Project SEP-CONACYT CB-2017-2018 (# A1-S-34563).
REFERENCES
- Anderson, D.; Burnham, K. 2004. Model Selection and Multi-Model Inference . v. 63 2nd ed. Springer-Verlag, New York, 488p.
- Baldeck, C.A.; Tupayachi, R.; Sinca, F.; Jaramillo, N.; Asner, G.P. 2016. Environmental drivers of tree community turnover in western Amazonian forests. Ecography, 39: 1089-1099.
- Balslev, H.; Luteyn, J.; Øllgaard, B.; Holm-Nielsen, L.B. 1987. Composition and structure of adjacent unflooded and floodplain forest in Amazonian Ecuador. Opera Botanica, 92: 37-57.
- Banon, G.P.R.; Banon, G.J.F.; Villamarín, F.; Arraut, E.M.; Moulatlet, G.M.; Rennó, C.D.; et al 2019. Predicting suitable nesting sites for the Black caiman (Melanosuchus niger Spix 1825) in the Central Amazon basin. Neotropical Biodiversity, 5: 47-59.
- Banon, L. 2013. Árvores de decisão aplicadas à extração automática de redes de dranagem Master´s dissertation, Instituto Nacional de Pesquisas Espaciais, Brazil. 87p. (http://urlib.net/8JMKD3MGP7W/3FB9RFE).
» http://urlib.net/8JMKD3MGP7W/3FB9RFE - Barton, K. 2018. Mu-MIn: Multi-model inference R package. (http://r-forge. r-project. org/projects/mumin/).
» http://r-forge. r-project. org/projects/mumin/ - Cámara-Leret, R.; Tuomisto, H.; Ruokolainen, K.; Balslev, H.; Munch Kristiansen, S. 2017. Modelling responses of western Amazonian palms to soil nutrients. Journal of Ecology, 105: 367-381.
- Campbell, D.G.; Daly, D.C.; Prance, G.T.; Maciel, U.N. 1986. Quantitative ecological inventory of terra firme and várzea tropical forest on the Rio Xingu, Brazilian Amazon. Brittonia, 38: 369-393.
- Chadwick, K.D.; Asner, G.P. 2016. Tropical soil nutrient distributions determined by biotic and hillslope processes. Biogeochemistry, 127: 273-289.
- Chauvel, A.; Lucas, Y.; Boulet, R. 1987. On the genesis of the soil mantle of the region of Manaus, Central Amazonia, Brazil. Experientia, 43: 234-241.
- Costa, F.R.C.; Magnusson, W.E.; Luizao, R.C. 2005. Mesoscale distribution patterns of Amazonian understorey herbs in relation to topography, soil and watersheds. Journal of Ecology, 93: 863-878.
- Costa, F.R.C.; Guillaumet, J.-L.; Lima, A.P.; Pereira, O.S. 2009. Gradients within gradients: The mesoscale distribution patterns of palms in a central Amazonian forest. Journal of Vegetation Science, 20: 69-78.
- De’ath, G. 1999. Extended dissimilarity: a method of robust estimation of ecological distances from high beta diversity data. Plant Ecology, 144: 191-199.
- Drucker, D.P.; Costa, F.R.C.; Magnusson, W.E. 2008. How wide is the riparian zone of small streams in tropical forests? A test with terrestrial herbs. Journal of Tropical Ecology, 24: 65-74.
- Fan, Y.; Li, H.; Miguez-Macho, G. 2013. Global patterns of groundwater table depth. Science, 339: 940-943.
- Fan, Y.; Clark, M.; Lawrence, D.M.; Swenson, S.; Band, L.E.; Brantley, S.L.; et al 2019. Hillslope hydrology in global change research and earth system modeling. Water Resources Research, 55: 1737-1772.
- Fontes, C.G.; Fine, P.V.A.; Wittmann, F.; Bittencourt, P.R.L.; Piedade, M.T.F.; Higuchi, N.; et al 2020. Convergent evolution of tree hydraulic traits in Amazonian habitats: implications for community assemblage and vulnerability to drought. New Phytologist, 228: 106-120.
- Frappart, F.; Papa, F.; Güntner, A.; Tomasella, J.; Pfeffer, J.; Ramillien, G.; et al 2019. The spatio-temporal variability of groundwater storage in the Amazon River Basin. Advances in Water Resources, 124: 41-52.
- Garcia, M.N.; Hu, J.; Domingues, T.F.; Groenendijk, P.; Oliveira, R.S.; Costa, F.R.C. 2022. Local hydrological gradients structure high intraspecific variability in plant hydraulic traits in two dominant central Amazonian tree species. Journal of Experimental Botany, 73: 939-952.
- Gentry, A.H. 1988. Changes in plant community diversity and floristic composition on environmental and geographical gradients. Annals of the Missouri Botanical Garden, 75: 1-34.
- Guimaraes, A.F.; de Souza, C.R.; Rosa, C.; Paulo dos Santos, J.; Teixeira, L.A.F.; Zanzini, L.P.; et al 2021. Small-scale environmental variations drive vegetation structure and diversity in Amazon riverine forests. Flora, 283: 151916.
- Haitjema, H.M.; Mitchell-Bruker, S. 2005. Are water tables a subdued replica of the topography? Ground Water, 43: 781-786.
- Higgins, M.A.; Ruokolainen, K.; Tuomisto, H.; Llerena, N.; Cardenas, G.; Phillips, O.L.; et al 2011. Geological control of floristic composition in Amazonian forests. Journal of Biogeography, 38: 2136-2149.
- Hoorn, C.; Wesselingh, F. 2011. Amazonia, Landscape and Species Evolution: A Look into the Past John Wiley & Sons, Chichester, 883p.
- Jirka, S.; McDonald, A.J.; Johnson, M.S.; Feldpausch, T.R.; Couto, E.G.; Riha, S.J. 2007. Relationships between soil hydrology and forest structure and composition in the southern Brazilian Amazon. Journal of Vegetation Science, 18: 183-194.
- Kahn, F.; de Castro, A. 1985. The palm community in a forest of central Amazonia, Brazil. Biotropica, 17: 210-216.
- Karger, D.N.; Conrad, O.; Böhner, J.; Kawohl, T.; Kreft, H.; Soria-Auza, R.W.; et al 2017. Climatologies at high resolution for the earth’s land surface areas. Scientific Data, 4: 170122.
- Lopez, O.R.; Kursar, T.A. 1999. Flood tolerance of four tropical tree species. Tree Physiology, 19: 925-932.
- Miguez-Macho, G.; Fan, Y. 2012. The role of groundwater in the Amazon water cycle: 1. Influence on seasonal streamflow, flooding and wetlands. Journal of Geophysical Research. Atmospheres, 117: 1-30.
- Moeslund, J.E.; Arge, L.; Bøcher, P.K.; Dalgaard, T.; Odgaard, M.V.; Nygaard, B.; et al 2013. Topographically controlled soil moisture is the primary driver of local vegetation patterns across a lowland region. Ecosphere, 4: 1-26.
- Moulatlet, G.M.; Costa, F.R.C.; Rennó, C.D.; Emilio, T.; Schietti, J. 2014. Local hydrological conditions explain floristic composition in lowland Amazonian forests. Biotropica, 46: 395-403.
- Nobre, A.D.; Cuartas, L.A.; Hodnett, M.; Rennó, C.D.; Rodrigues, G.; Silveira, A.; et al 2011. Height Above the Nearest Drainage - a hydrologically relevant new terrain model. Journal of Hydrology, 404: 13-29.
- Oliveira, R.S.; Costa, F.R.C.; van Baalen, E.; de Jonge, A.; Bittencourt, P.R.; Almanza, Y.; et al 2018. Embolism resistance drives the distribution of Amazonian rainforest tree species along hydro-topographic gradients. New Phytologist, 221:1457-1465.
- Osher, L.J.; Buol, S.W. 1998. Relationship of soil properties to parent material and landscape position in eastern Madre de Dios, Peru. Geoderma, 83: 143-166.
- Parent, C.; Capelli, N.; Berger, A.; Crèvecoeur, M.; Dat, J.F. 2008. An overview of plant responses to soil waterlogging. Plant Stress, 2: 20-27.
- Poulsen, A.D.; Tuomisto, H.; Balslev, H. 2006. Edaphic and floristic variation within a 1-ha plot of lowland Amazonian rain forest. Biotropica, 38: 468-478.
- R Core Team. 2020. R: A language and environment for statistical computing R Foundation for Statistical Computing, Vienna, Austria. (http://www.R-project.org/).
» http://www.R-project.org/ - Rennó, C.D.; Nobre, A.D.; Cuartas, L.A.; Soares, J.V.; Hodnett, M.G.; Tomasella, J.; et al 2008. HAND, a new terrain descriptor using SRTM-DEM: Mapping terra-firme rainforest environments in Amazonia. Remote Sensing of Environment, 112: 3469-3481.
- Rocchini, D.; Boyd, D.S.; Féret, J.-B.; Foody, G.M.; He, K.S.; Lausch, A.; et al 2016. Satellite remote sensing to monitor species diversity: potential and pitfalls. Remote Sensing in Ecology and Conservation, 2: 25-36.
- Rodrigues, F.R. de O.; Costa, F.R.C. 2012. Litter as a filter of emergence for herbaceous seedlings and sporophytes in central Amazonia. Journal of Tropical Ecology, 28: 445-452.
- Schietti, J.; Emilio, T.; Rennó, C.D.; Drucker, D.P.; Costa, F.R.C.; Nogueira, A.; et al 2014. Vertical distance from drainage drives floristic composition changes in an Amazonian rainforest. Plant Ecology & Diversity, 7: 241-253.
- Strahler, A.N. 1957. Quantitative analysis of watershed geomorphology. Eos, Transactions American Geophysical Union, 38: 913-920.
- Svenning, J.-C. 1999. Recruitment of tall arborescent palms in the Yasuní National Park, Amazonian Ecuador: are large treefall gaps important? Journal of Tropical Ecology, 15: 355-366.
- Tuomisto, H.; Ruokolainen, K. 1994. Distribution of Pteridophyta and Melastomataceae along an edaphic gradient in an Amazonian rain forest. Journal of Vegetation Science, 5: 25-34.
- Tuomisto, H.; Poulsen, A.D. 2000. Pteridophyte diversity and species composition in four Amazonian rain forests. Journal of Vegetation Science, 11: 383-396.
- Tuomisto, H.; Poulsen, A.D.; Moran, R.C. 1998. Edaphic distribution of some species of the fern genus Adiantum in western Amazonia. Biotropica, 30: 392-399.
- Tuomisto, H.; Ruokolainen, K.; Kalliola, R.; Linna, A.; Danjoy, W.; Rodriguez, Z. 1995. Dissecting Amazonian biodiversity. Science, 269: 63-66.
- Tuomisto, H.; Moulatlet, G.M.; Balslev, H.; Emilio, T.; Figueiredo, F.O.G.; Pedersen, D.; et al 2016. A compositional turnover zone of biogeographical magnitude within lowland Amazonia. Journal of Biogeography, 43: 2400-2411.
- Tuomisto, H.; Van doninck, J.; Ruokolainen, K.; Moulatlet, G.M.; Figueiredo, F.O.G.; Sirén, A.; et al 2019. Discovering floristic and geoecological gradients across Amazonia. Journal of Biogeography, 46: 1734-1748.
- Valeriano, M.M.; Kuplich, T.M.; Storino, M.; Amaral, B.D.; Mendes Jr., J.N.; Lima, D.J. 2006. Modeling small watersheds in Brazilian Amazonia with shuttle radar topographic mission-90 m data. Computers & Geosciences, 32: 1169-1181.
- Vormisto, J.; Tuomisto, H.; Oksanen, J. 2004. Palm distribution patterns in Amazonian rainforests: What is the role of topographic variation? Journal of Vegetation Science, 15: 485-494.
- Wittmann, F.; Householder, E.; Piedade, M.T.F.; de Assis, R.L.; Schöngart, J.; Parolin, P.; et al 2013. Habitat specifity, endemism and the neotropical distribution of Amazonian white-water floodplain trees. Ecography, 36: 690-707.
- Zuquim, G.; Tuomisto, H.; Costa, F.R.C.; Prado, J.; Magnusson, W.E.; Pimentel, T.; et al 2012. Broad scale distribution of ferns and lycophytes along environmental gradients in central and northern Amazonia, Brazil. Biotropica, 44: 752-762.
- Zuquim, G.; Tuomisto, H.; Chaves, P.P.; Emilio, T.; Moulatlet, G.M.; Ruokolainen, K.; et al 2021. Revealing floristic variation and map uncertainties for different plant groups in western Amazonia. Journal of Vegetation Science, 32: e13081.
-
CITE AS:
Moulatlet, G.M.; Rennó, C. D.; Figueiredo, F.O.G.; Ruokolainen, K.; Banon, L.; Emilio, T.; Balslev, H.; Tuomisto, H. 2022. The role of topographic-derived hydrological variables in explaining plant species distributions in Amazonia. Acta Amazonica 52: 218-228.
SUPPLEMENTARY MATERIAL
(only available in the electronic version)
Moulatlet et al. The role of topographic-derived hydrological variables in explaining plant species distributions in Amazonia
Comparison between generalized linear mixed models built at species-level with DEM-derived topographic metrics (Hydrology) and with field measured topography (Topography) from transects along the middle Juruá River (western Brazilian Amazon). A - Distribution of AIC values for each plant group; B - Comparison between AICs for each modelled species. Pearson’s R = 0.99.
Comparison between generalized linear mixed models built at species-level with DEM-derived topographic metrics (Hydrology) and with field measured topography (Topography) for each species of each plant group sampled in transects along the middle Juruá River (western Brazilian Amazon). The results of Hydrological vs. Topographic models are connected for each individual species.
Presence of one selected species (vertical bars) of each plant group sampled in the set of transects along the middle Juruá River (western Brazilian Amazon). The transect (250 m total length) is divided in 20 subunits of 25 m each. Solid lines correspond to the field measured topographic profiles. Dashed lines correspond to the calculated HAND values. Species occurrence probabilities according to the field measured topography and to HAND are shown in the right column. Both variables are standardized by subtracting the mean and dividing by the standard deviation.
Number of plant species used in each model and the corresponding number of transects and subunits within transects along the middle Juruá River (western Brazilian Amazon). Values within parenthesis indicate the number of species significantly associated to at least one hydrological variable in GLMM.
Summary of generalized linear mixed models (GLMM) for individual plant species sampled in transects along the middle Juruá River (western Brazilian Amazon). Coefficients of each model (at logit scale) are shown with the respective standard errors (SE), p-values and pseudo-R2. The numbers 1 and 2 after each variable’s name refer to the linear and squared terms, respectively. AIC values are shown for the models with DEM-derived topographic variables (AIC DEM) and for the models with field measured topography (AIC Topo).
Edited by
ASSOCIATE EDITOR:
Publication Dates
-
Publication in this collection
12 Sept 2022 -
Date of issue
Jul-Sep 2022
History
-
Received
05 Jan 2022 -
Accepted
16 June 2022