Predicting current and future spatial community patterns of plant functional traits

Community-level patterns of functional traits relate to community assembly and ecosystem functioning. By modelling the changes of different indices describing such patterns - trait means, extremes and diversity in communities - as a function of abiotic gradients, we could understand their drivers and build projections of the impact of global change on the functional components of biodiversity. We used five plant functional traits (vegetative height, specific leaf area, leaf dry matter content, leaf nitrogen content and seed mass) and non-woody vegetation plots to model several indices depicting community-level patterns of functional traits from a set of abiotic environmental variables (topographic, climatic and edaphic) over contrasting environmental conditions in a mountainous landscape. We performed a variation partitioning analysis to assess the relative importance of these variables for predicting patterns of functional traits in communities, and projected the best models under several climate change scenarios to examine future potential changes in vegetation functional properties. Not all indices of trait patterns within communities could be modelled with the same level of accuracy: the models for mean and extreme values of functional traits provided substantially better predictive accuracy than the models calibrated for diversity indices. Topographic and climatic factors were more important predictors of functional trait patterns within communities than edaphic predictors. Overall, model projections forecast an increase in mean vegetation height and in mean specific leaf area following climate warming. This trend was important at mid elevation particularly between 1000 and 2000 m asl. With this study we showed that topographic, climatic and edaphic variables can successfully model descriptors of community-level patterns of plant functional traits such as mean and extreme trait values. However, which factors determine the diversity of functional traits in plant communities remains unclear and requires more investigations.

Using species richness and functional traits predictions to constrain assemblage predictions from stacked species distribution models

Aim: Modelling species at the assemblage level is required to make effective forecast of global change impacts on diversity and ecosystem functioning. Community predictions may be achieved using macroecological properties of communities (MEM), or by stacking of individual species distribution models (S-SDMs). To obtain more realistic predictions of species assemblages, the SESAM framework suggests applying successive filters to the initial species source pool, by combining different modelling approaches and rules. Here we provide a first test of this framework in mountain grassland communities. Location: The western Swiss Alps. Methods: Two implementations of the SESAM framework were tested: a "Probability ranking" rule based on species richness predictions and rough probabilities from SDMs, and a "Trait range" rule that uses the predicted upper and lower bound of community-level distribution of three different functional traits (vegetative height, specific leaf area and seed mass) to constraint a pool of environmentally filtered species from binary SDMs predictions. Results: We showed that all independent constraints expectedly contributed to reduce species richness overprediction. Only the "Probability ranking" rule allowed slightly but significantly improving predictions of community composition. Main conclusion: We tested various ways to implement the SESAM framework by integrating macroecological constraints into S-SDM predictions, and report one that is able to improve compositional predictions. We discuss possible improvements, such as further improving the causality and precision of environmental predictors, using other assembly rules and testing other types of ecological or functional constraints.

Variation in habitat suitability does not always relate to variation in species' plant functional traits.

Habitat suitability models, which relate species occurrences to environmental variables, are assumed to predict suitable conditions for a given species. If these models are reliable, they should relate to change in plant growth and function. In this paper, we ask the question whether habitat suitability models are able to predict variation in plant functional traits, often assumed to be a good surrogate for a species' overall health and vigour. Using a thorough sampling design, we show a tight link between variation in plant functional traits and habitat suitability for some species, but not for others. Our contrasting results pave the way towards a better understanding of how species cope with varying habitat conditions and demonstrate that habitat suitability models can provide meaningful descriptions of the functional niche in some cases, but not in others.

Plant functional and phylogenetic turnover correlate with climate and land use in the Western Swiss Alps

Understanding the relative importance of historical and environmental processes in the structure and composition of communities is one of the longest quests in ecological research. Increasingly, researchers are relying on the functional and phylogenetic β-diversity of natural communities to provide concise explanations on the mechanistic basis of community assembly and the drivers of trait variation among species. The present study investigated how plant functional and phylogenetic β-diversity change along key environmental and spatial gradients in the Western Swiss Alps. Methods Using the quadratic diversity measure based on six functional traits: specific leaf area (SLA), leaf dry matter content (LDMC), plant height (H), leaf carbon content (C), leaf nitrogen content (N), and leaf carbon to nitrogen content (C/N) alongside a species-resolved phylogenetic tree, we relate variations in climate, spatial geographic, land use and soil gradients to plant functional and phylogenetic turnover in mountain communities of the Western Swiss Alps. Important findings Our study highlights two main points. First, climate and land use factors play an important role in mountain plant community turnover. Second, the overlap between plant functional and phylogenetic turnover along these gradients correlates with the low phylogenetic signal in traits, suggesting that in mountain landscapes, trait lability is likely an important factor in driving plant community assembly. Overall, we demonstrate the importance of climate and land use factors in plant functional and phylogenetic community turnover, and provide valuable complementary insights into understanding patterns of β-diversity along several ecological gradients.

Species distribution models reveal apparent competitive and facilitative effects of a dominant species on the distribution of tundra plants

Abiotic factors are considered strong drivers of species distribution and assemblages. Yet these spatial patterns are also influenced by biotic interactions. Accounting for competitors or facilitators may improve both the fit and the predictive power of species distribution models (SDMs). We investigated the influence of a dominant species, Empetrum nigrum ssp. hermaphroditum, on the distribution of 34 subordinate species in the tundra of northern Norway. We related SDM parameters of those subordinate species to their functional traits and their co-occurrence patterns with E. hermaphroditum across three spatial scales. By combining both approaches, we sought to understand whether these species may be limited by competitive interactions and/or benefit from habitat conditions created by the dominant species. The model fit and predictive power increased for most species when the frequency of occurrence of E. hermaphroditum was included in the SDMs as a predictor. The largest increase was found for species that 1) co-occur most of the time with E. hermaphroditum, both at large (i.e. 750 m) and small spatial scale (i.e. 2 m) or co-occur with E. hermaphroditum at large scale but not at small scale and 2) have particularly low or high leaf dry matter content (LDMC). Species that do not co-occur with E. hermaphroditum at the smallest scale are generally palatable herbaceous species with low LDMC, thus showing a weak ability to tolerate resource depletion that is directly or indirectly induced by E. hermaphroditum. Species with high LDMC, showing a better aptitude to face resource depletion and grazing, are often found in the proximity of E. hermaphroditum. Our results are consistent with previous findings that both competition and facilitation structure plant distribution and assemblages in the Arctic tundra. The functional and co-occurrence approaches used were complementary and provided a deeper understanding of the observed patterns by refinement of the pool of potential direct and indirect ecological effects of E. hermaphroditum on the distribution of subordinate species. Our correlative study would benefit being complemented by experimental approaches.

Predicting spatial patterns of plant biodiversity : from species to communities

Understanding the distribution and composition of species assemblages and being able to predict them in space and time are highly important tasks io investigate the fate of biodiversity in the current global changes context. Species distribution models are tools that have proven useful to predict the potential distribution of species by relating their occurrences to environmental variables. Species assemblages can then be predicted by combining the prediction of individual species models. In the first part of my thesis, I tested the importance of new environmental predictors to improve species distribution prediction. I showed that edaphic variables, above all soil pH and nitrogen content could be important in species distribution models. In a second chapter, I tested the influence of different resolution of predictors on the predictive ability of species distribution models. I showed that fine resolution predictors could ameliorate the models for some species by giving a better estimation of the micro-topographic condition that species tolerate, but that fine resolution predictors for climatic factors still need to be ameliorated. The second goal of my thesis was to test the ability of empirical models to predict species assemblages' characteristics such as species richness or functional attributes. I showed that species richness could be modelled efficiently and that the resulting prediction gave a more realistic estimate of the number of species than when obtaining it by stacking outputs of single species distribution models. Regarding the prediction of functional characteristics (plant height, leaf surface, seed mass) of plant assemblages, mean and extreme values of functional traits were better predictable than indices reflecting the diversity of traits in the community. This approach proved interesting to understand which environmental conditions influence particular aspects of the vegetation functioning. It could also be useful to predict climate change impacts on the vegetation. In the last part of my thesis, I studied the capacity of stacked species distribution models to predict the plant assemblages. I showed that this method tended to over-predict the number of species and that the composition of the community was not predicted exactly either. Finally, I combined the results of macro- ecological models obtained in the preceding chapters with stacked species distribution models and showed that this approach reduced significantly the number of species predicted and that the prediction of the composition is also ameliorated in some cases. These results showed that this method is promising. It needs now to be tested on further data sets. - Comprendre la manière dont les plantes se répartissent dans l'environnement et s'organisent en communauté est une question primordiale dans le contexte actuel de changements globaux. Cette connaissance peut nous aider à sauvegarder la diversité des espèces et les écosystèmes. Des méthodes statistiques nous permettent de prédire la distribution des espèces de plantes dans l'espace géographique et dans le temps. Ces modèles de distribution d'espèces, relient les occurrences d'une espèce avec des variables environnementales pour décrire sa distribution potentielle. Cette méthode a fait ses preuves pour ce qui est de la prédiction d'espèces individuelles. Plus récemment plusieurs tentatives de cumul de modèles d'espèces individuelles ont été réalisées afin de prédire la composition des communautés végétales. Le premier objectif de mon travail est d'améliorer les modèles de distribution en testant l'importance de nouvelles variables prédictives. Parmi différentes variables édaphiques, le pH et la teneur en azote du sol se sont avérés des facteurs non négligeables pour prédire la distribution des plantes. Je démontre aussi dans un second chapitre que les prédicteurs environnementaux à fine résolution permettent de refléter les conditions micro-topographiques subies par les plantes mais qu'ils doivent encore être améliorés avant de pouvoir être employés de manière efficace dans les modèles. Le deuxième objectif de ce travail consistait à étudier le développement de modèles prédictifs pour des attributs des communautés végétales tels que, par exemple, la richesse en espèces rencontrée à chaque point. Je démontre qu'il est possible de prédire par ce biais des valeurs de richesse spécifiques plus réalistes qu'en sommant les prédictions obtenues précédemment pour des espèces individuelles. J'ai également prédit dans l'espace et dans le temps des caractéristiques de la végétation telles que sa hauteur moyenne, minimale et maximale. Cette approche peut être utile pour comprendre quels facteurs environnementaux promeuvent différents types de végétation ainsi que pour évaluer les changements à attendre au niveau de la végétation dans le futur sous différents régimes de changements climatiques. Dans une troisième partie de ma thèse, j'ai exploré la possibilité de prédire les assemblages de plantes premièrement en cumulant les prédictions obtenues à partir de modèles individuels pour chaque espèce. Cette méthode a le défaut de prédire trop d'espèces par rapport à ce qui est observé en réalité. J'ai finalement employé le modèle de richesse en espèce développé précédemment pour contraindre les résultats du modèle d'assemblage de plantes. Cela a permis l'amélioration des modèles en réduisant la sur-prédiction et en améliorant la prédiction de la composition en espèces. Cette méthode semble prometteuse mais de nouveaux tests sont nécessaires pour bien évaluer ses capacités.

Turnover of plant lineages shapes herbivore phylogenetic beta diversity along ecological gradients.

Understanding drivers of biodiversity patterns is of prime importance in this era of severe environmental crisis. More diverse plant communities have been postulated to represent a larger functional trait-space, more likely to sustain a diverse assembly of herbivore species. Here, we expand this hypothesis to integrate environmental, functional and phylogenetic variation of plant communities as factors explaining the diversity of lepidopteran assemblages along elevation gradients in the Swiss Western Alps. According to expectations, we found that the association between butterflies and their host plants is highly phylogenetically structured. Multiple regression analyses showed the combined effect of climate, functional traits and phylogenetic diversity in structuring butterfly communities. Furthermore, we provide the first evidence that plant phylogenetic beta diversity is the major driver explaining butterfly phylogenetic beta diversity. Along ecological gradients, the bottom up control of herbivore diversity is thus driven by phylogenetically structured turnover of plant traits as well as environmental variables.

Improving plant functional groups for dynamic models of biodiversity: at the crossroads between functional and community ecology

The pace of on-going climate change calls for reliable plant biodiversity scenarios. Traditional dynamic vegetation models use plant functional types that are summarized to such an extent that they become meaningless for biodiversity scenarios. Hybrid dynamic vegetation models of intermediate complexity (hybrid-DVMs) have recently been developed to address this issue. These models, at the crossroads between phenomenological and process-based models, are able to involve an intermediate number of well-chosen plant functional groups (PFGs). The challenge is to build meaningful PFGs that are representative of plant biodiversity, and consistent with the parameters and processes of hybrid-DVMs. Here, we propose and test a framework based on few selected traits to define a limited number of PFGs, which are both representative of the diversity (functional and taxonomic) of the flora in the Ecrins National Park, and adapted to hybrid-DVMs. This new classification scheme, together with recent advances in vegetation modeling, constitutes a step forward for mechanistic biodiversity modeling.

Functional homogenization of bumblebee communities in alpine landscapes under projected climate change

Background: Bumblebees represent an active pollinator group in mountain regions and assure the pollination of many different plant species from low to high elevations. Plant-pollinator interactions are mediated by functional traits. Shift in bumblebee functional structure under climate change may impact plant-pollinator interactions in mountains. Here, we estimated bumblebee upward shift in elevation, community turnover, and change in functional structure under climate change. Method: We sampled bumblebee species at 149 sites along the elevation gradient. We used stacked species distribution models (S-SDMs) forecasted under three climate change scenarios (A2, A1B, RCP3PD) to model the potential distribution of the Bombus species. Furthermore, we used species proboscis length measurements to assess the functional change in bumblebee assemblages along the elevation gradient. Results: We found species-specific response of bumblebee species to climate change. Species differed in their predicted rate of range contraction and expansion. Losers were mainly species currently restricted to high elevation. Under the most severe climate change scenarios (A2), we found a homogenization of proboscis length structure in bumblebee communities along the elevation gradient through the upward colonization of high elevation by species with longer proboscides. Conclusions: Here, we show that in addition to causing the shift in the distribution of bumblebee species, climate change may impact the functional structure of communities. The colonization of high elevation areas by bumblebee species with long proboscides may modify the structure of plant-pollination interaction networks by increasing the diversity of pollination services at high elevation.

Heritability of left ventricular structure and function in Caucasian families.

AIMS: The aim of this study was to investigate the heritability as well as genetic and environmental correlations of left ventricular (LV) structural and functional traits in complex pedigrees of a Caucasian population. METHODS AND RESULTS: We randomly recruited 459 white European subjects from 52 families (50% women; mean age 45 years). LV structure was measured by M-mode and 2D echocardiography and LV function was measured by conventional Doppler and tissue Doppler imaging (TDI). Other measurements included blood pressure, anthropometric, and biochemical measurements. We estimated the heritability of LV traits while adjusting for covariables, including sex, age, body height and weight, systolic and diastolic blood pressures, and heart rate. With full adjustment, heritability of LV mass was 0.23 (P= 0.025). The TDI-derived mitral annular velocities Ea and Aa showed moderate heritability (h(2)= 0.36 and 0.53, respectively), whereas the mitral inflow A peak had weak heritability (h(2) = 0.25) and the E peak was not heritable (h(2) = 0.11). We partitioned the total phenotypic correlation when it reached significance, into a genetic and an environmental component. The genetic correlations were 0.61 between the E and Ea peaks and 0.90 between the A and Aa peaks. CONCLUSION: Our study demonstrated moderate heritability for LV mass as well as the mitral annular Ea and Aa peaks. We also found significant genetic correlations between the E and Ea peaks and between the A and Aa peaks. Our current findings support the ongoing research to map and detect genetic variants that contribute to the variation in LV mass and other LV structural and functional phenotypes.

Functionally and phylogenetically diverse plant communities key to soil biota

Recent studies assessing the role of biological diversity for ecosystem functioning indicate that the diversity of functional traits and the evolutionary history of species in a community, not the number of taxonomic units, ultimately drives the biodiversity-ecosystem-function relationship. Here, we simultaneously assessed the importance of plant functional trait and phylogenetic diversity as predictors of major trophic groups of soil biota (abundance and diversity), six years from the onset of a grassland biodiversity experiment. Plant functional and phylogenetic diversity were generally better predictors of soil biota than the traditionally used species or functional group richness. Functional diversity was a reliable predictor for most biota, with the exception of soil microorganisms, which were better predicted by phylogenetic diversity. These results provide empirical support for the idea that the diversity of plant functional traits and the diversity of evolutionary lineages in a community are important for maintaining higher abundances and diversity of soil communities.

On the fate of sexual traits under asexuality.

Environmental shifts and life-history changes may result in formerly adaptive traits becoming non-functional or maladaptive. In the absence of pleiotropy and other constraints, such traits may decay as a consequence of neutral mutation accumulation or selective processes, highlighting the importance of natural selection for adaptations. A suite of traits are expected to lose their adaptive function in asexual organisms derived from sexual ancestors, and the many independent transitions to asexuality allow for comparative studies of parallel trait maintenance versus decay. In addition, because certain traits, notably male-specific traits, are usually not exposed to selection under asexuality, their decay would have to occur as a consequence of drift. Selective processes could drive the decay of traits associated with costs, which may be the case for the majority of sexual traits expressed in females. We review the fate of male and female sexual traits in 93 animal lineages characterized by asexual reproduction, covering a broad taxon range including molluscs, arachnids, diplopods, crustaceans and eleven different hexapod orders. Many asexual lineages are still able occasionally to produce males. These asexually produced males are often largely or even fully functional, revealing that major developmental pathways can remain quiescent and functional over extended time periods. By contrast, for asexual females, there is a parallel and rapid decay of sexual traits, especially of traits related to mate attraction and location, as expected given the considerable costs often associated with the expression of these traits. The level of decay of female sexual traits, in addition to asexual females being unable to fertilize their eggs, would severely impede reversals to sexual reproduction, even in recently derived asexual lineages. More generally, the parallel maintenance versus decay of different trait types across diverse asexual lineages suggests that neutral traits display little or no decay even after extended periods under relaxed selection, while extensive decay for selected traits occurs extremely quickly. These patterns also highlight that adaptations can fix rapidly in natural populations of asexual organisms, in spite of their mode of reproduction.

Genetics of decayed sexual traits in a parasitoid wasp with endosymbiont-induced asexuality.

Trait decay may occur when selective pressures shift, owing to changes in environment or life style, rendering formerly adaptive traits non-functional or even maladaptive. It remains largely unknown if such decay would stem from multiple mutations with small effects or rather involve few loci with major phenotypic effects. Here, we investigate the decay of female sexual traits, and the genetic causes thereof, in a transition from haplodiploid sexual reproduction to endosymbiont-induced asexual reproduction in the parasitoid wasp Asobara japonica. We take advantage of the fact that asexual females cured of their endosymbionts produce sons instead of daughters, and that these sons can be crossed with sexual females. By combining behavioral experiments with crosses designed to introgress alleles from the asexual into the sexual genome, we found that sexual attractiveness, mating, egg fertilization and plastic adjustment of offspring sex ratio (in response to variation in local mate competition) are decayed in asexual A. japonica females. Furthermore, introgression experiments revealed that the propensity for cured asexual females to produce only sons (because of decayed sexual attractiveness, mating behavior and/or egg fertilization) is likely caused by recessive genetic effects at a single locus. Recessive effects were also found to cause decay of plastic sex-ratio adjustment under variable levels of local mate competition. Our results suggest that few recessive mutations drive decay of female sexual traits, at least in asexual species deriving from haplodiploid sexual ancestors.

Functional and molecular evolutionary analysis of the "Brevix Radix" gene family in mono-and dicotyledonous plants

ABSTRACTIn contrast to animals, plants cannot move from their place of birth and, therefore, need to adapt to their particular habitat in order to survive. Thus, plant development is remarkably plastic, making plants an ideal system for the isolation of genes that account for intraspecific natural variation and possibly environmental adaptation. However, to date, this approach mostly identified null alleles and missed mutations with subtle effects. For instance, BREVIS RADIX (BRX) has been isolated as a key regulator of root growth through a naturally occurring loss-of-function allele in the Arabidopsis thaliana accession Uk-1 and is the founding member of a highly-conserved plant-specific gene family.In this work, we show that a strong selective pressure is acting on the BRX gene family and dates back before the monocot-dicot divergence. However, functional diversification is observed mainly in dicotyledon BRX family genes and is correlated with acceleration in the evolutionary rates in the N-terminal regions. Population genetic data revealed that BRX is highly conserved across Arabidopsis accessions and presents signatures of adaptation. Interestingly, a seven amino acid deletion polymorphism in BRX sequence was found in a few accessions, which seems to be responsible for their enhanced primary root growth. Nevertheless, BRX might not only be active in the root, as suggested by its expression in the shoot. Indeed, leaves and cotyledons of brx mutants are significantly smaller than wild- type. This phenotype is a direct consequence of the absence of BRX function in the shoot rather than an indirect effect of an altered root system growth. Interestingly, cotyledons of brx plants reflect the same physiological defects as the root. Moreover, phenotypes in BRX gain-of-function plants, such as epinastic leaves and increased epidermal cell size, could be associated with an increase in leaf brassinosteroid content.Collectively, these results indicate that BRX contributes to local adaptation by ubiquitously regulating plant growth, probably through the modulation of brassinosteroid biosynthesis.RÉSUMÉContrairement à la plupart des animaux, les plantes ne peuvent se mouvoir et doivent ainsi s'adapter à leur environnement pour survivre. Pour cette raison, elles représentent un système idéal pour l'identification de gènes contribuant à la variation naturelle intra- spécifique, ainsi qu'à l'adaptation. Cependant, cette approche a, jusqu'à présent, surtout permis d'isoler des allèles nuls et non des mutations conférant des effets plus subtiles. C'est le cas du gène Β REVIS RADIX (BRX), un régulateur clé de la croissance racinaire, qui a été identifié grâce à un allèle non-fonctionnel présent dans l'accession naturelle d'Arabidopsis thaliana Uk-1. BRX et ses homologues des plantes mono- et dicotylédones forment une famille très conservée et spécifique aux plantes.Dans ce travail, nous démontrons que la famille de gènes BRX est soumise à une forte pression de sélection qui remonte avant la divergence entre mono- et dicotylédones. Cependant, une diversification fonctionnelle a été observée chez les gènes des dicotylédones et corrèle avec une accélération de la vitesse d'évolution dans leur région N- terminale. Une analyse génétique de différentes accessions naturelles d'Arabidopsis a révélé que BRX est hautement conservé et présente des signatures d'adaptation. Remarquablement, un polymorphisme de délétion de sept acides aminés a été détecté dans quelques accessions et a pour conséquence une plus forte croissance de la racine primaire. Néanmoins, il semble que le rôle de BRX ne se limite pas qu'à la racine, comme indiqué par son expression dans les parties aériennes de la plante. En effet, les mutants brx présentent des cotylédons et des feuilles significativement plus petits que le type sauvage, une conséquence directe de l'absence d'activité de BRX dans ces organes. Nous avons aussi noté que les cotylédons des mutants brx, à l'instar des racines, ont une perception altérée de l'auxine et peuvent être complémentés par l'application exogène de brassinostéroïdes. De plus, dans des plantes présentant un gain de fonction BRX, les feuilles sont épinastiques et les cellules de leur épiderme plus grandes. Ces phénotypes sont accompagnés d'une augmentation de la concentration de brassinostéroïdes dans les feuilles. Conjointement, ces résultats démontrent que BRX contribue à une adaptation locale de la plante par la régulation générale de sa croissance, probablement en modulant la biosynthèse des brassinostéroïdes.

Plant traits co-vary with altitude in grasslands and forests in the European Alps

Biological traits that are advantageous under specific ecological conditions should be present in a large proportion of the species within an ecosystem, where those specific conditions prevail. As climatic conditions change, the frequency of certain traits in plant communities is expected to change with increasing altitude. We examined patterns of change for 13 traits in 120 exhaustive inventories of plants along five altitudinal transects (520-3100 m a.s.l.) in grasslands and in forests in western Switzerland. The traits selected for study represented the occupation of space, photosynthesis, reproduction and dispersal. For each plot, the mean trait values or the proportions of the trait states were weighted by species cover and examined in relation to the first axis of a PCA based on local climatic conditions. With increasing altitude in grasslands, we observed a decrease in anemophily and an increase in entomophily complemented by possible selfing; a decrease in diaspores with appendages adapted to ectozoochory, linked to a decrease in achenes and an increase in capsules. In lowlands, pollination and dispersal are ensured by wind and animals. However, with increasing altitude, insects are mostly responsible for pollination, and wind becomes the main natural dispersal vector. Some traits showed a particularly marked change in the alpine belt (e.g., the increase of capsules and the decrease of achenes), confirming that this belt concentrates particularly stressful conditions to plant growth and reproduction (e.g. cold, short growing season) that constrain plants to a limited number of strategies. One adaptation to this stress is to limit investment in dispersal by producing capsules with numerous, tiny seeds that have appendages limited to narrow wings. Forests displayed many of the trends observed in grasslands but with a reduced variability that is likely due to a shorter altitudinal gradient.