Changes in Allometric Attributes and Biomass of Forests and Woodlands across an Altitudinal and Rainfall Gradient : What Are the Implications of Increasing Seasonality due to Anthropogenic Climate Change ?

Canonical correspondence analysis and linear regressions were used to relate height, diameter, and dispersion measurements of 36,380 stems from 197 species recorded in 2,341 plots against both climatic and landscape variables. Above groundbiomass increased in wetter and cooler locations that ameliorate the seasonal rainfall deficits. Taller and greater diameter trees with lower wood densities occur at higher altitudes. Differences between locations are based on a change in the composition of species rather than a change in the allometric properties within a species. The results support the hydraulic limitation and species packing hypotheses. These interrelationships may be affected by the interactions of fire frequency and drought which are a common feature of much of the study area. Under current climate change scenarios it is likely that therewill be a reduction in above ground biomass, the number of stems per hectare, average height, average diameter, and basal area due to increasing seasonality of rainfall, temperatures, and the intensity and frequency of fires. The largest of trees are likely to be removed early due to their inability to cope with increased drought stress. The results suggest a marked reduction in carbon storage will occur across the study region in eastern Australia.


Introduction
Central to functional ecology is the ability to link changes in species traits to ecosystem processes [1].The allometric traits of trees determine the capture of resources, affect the transport of those resources throughout the plant, and also affect resistance to mechanical damage.These traits are indicative of survival versus growth life history trade-offs [2][3][4].Differences in allometry potentially structure coexistence between species [5].However, the way in which these allometric traits change across environmental gradients has been poorly investigated in any detail across diverse habitats [4,6].Cumulatively, allometric traits also structure the above ground biomass obtainable by a species and collectively the community in which it occurs.Therefore, knowledge of the spatial patterns of allometric traits allows understanding of changes in ecosystem function [2,7,8].
Altitude is an easily measurable surrogate for a number of environmental factors and thus is often used to measure change within assemblages including structure, composition, and species traits.Altitude has been shown to correlate with variation in tree allometric properties due to the concomitant changes in rainfall, air, and soil temperatures and differences in nutrient availability and soil chemistry [9].Increasing altitude not only affects climate but also is associated with a decrease in atmospheric pressure which directly affects photosynthesis [10].Biomass has been found to peak in cool wet environments associated with higher altitudes [9].In contrast, Larjavaara and Muller-Landau [8] found above ground biomass declined with increasing altitude.Culmsee et al. [9] noted a number of changes in forest structure including a decrease in tree height with increasing altitude.Lines et al. [4] while investigating intra-and interspecific changes in allometric traits with altitude discovered that taller species were associated with cooler higher altitude regions.Conversely, within a species, individuals were taller in warmer low altitudes [4].Slik et al. [11] attributed stunting in height with altitude to lower growth and productivity at colder temperatures.
Slik et al. [11] found that both basal area and above ground biomass were positively correlated with annual rainfall but negatively with increasing rainfall seasonality.Within moist and wet tropical forests of the Amazon Basin above ground biomass was positively correlated with increasing precipitation of the driest quarter [12].Keith et al. [13] found a peak in biomass in wet and cool regions.In an intercontinental study Banin et al. [3] correlated higher above ground biomass with seasonal climates where rainfall in the wettest period was highest.They suggested that higher rainfall within the wettest period allowed the soil profile to remain moist through the drier season and that the seasonal rainfall regime reduced hydraulic constraints on trees allowing them to achieve greater biomass [3].Culmsee et al. [9] while investigating a number of forest types showed that neither temperature nor precipitation of the driest quarter was consistently correlated with above ground biomass suggesting that only in relatively dry systems is precipitation a limiting factor.Support for the limiting role of precipitation on biomass in arid and semiarid systems was also found.
Above ground biomass accumulation and basal area are influenced by the presence of large diameter trees within a site [14,15].This led Midgley et al. [14] to propose the packing hypothesis which states that plot basal area and biomass will be strongly related to size of the largest individual because larger individuals pack space more efficiently.Stegen et al. [6] argued that if the largest tree was strongly correlated with biomass it is fundamentally through the influence of climate on the maximum individual tree mass.In open forests the abundance of trees may be below their theoretical maximum density and their canopies more open due to resource limitation allowing small individuals more access to resources, thus affecting the size distribution of stems and their density across the community [16].
Most investigations into above ground biomass have been carried out within tropical closed forests (e.g., [3,7,9,11,[17][18][19][20]) with investigations into these attributes within temperate woodlands being infrequent [5,13,21].Tropical closed forest studies have largely focused on competition for light being a major limiting factor structuring allometric attributes of species and associated above ground biomass.Woodlands and forests away from tropical regions will have reduced light due to the more oblique nature of solar radiation; however, they are still likely to be decoupled from the constraints of limited light.Thus other environmentally limiting factors are likely to be of importance [3,4,16].
Aiba and Nakashizuka [5] argued that relationships between tree properties and ecology are sensitive to which species are studied and argued that the number of species needs to be large enough to allow for a general understanding of patterns at the community level.Most above ground biomass studies have been based on correlations restricted to single forest types, limited in geographic extent and with few samples [4,6].Lines et al. [4] also argued that many studies have underused computation methods that effectively separate covarying correlates.
Using an extensive dataset of many species collected across a large regional area I attempt to assess the major factors influencing tree AGB and associated allometric properties.What factors are associated with AGB within temperate forests and woodlands within a semiarid to temperate summer dominant rainfall zone is investigated.It is assessed whether altitude, temperature, and rainfall are the major influences as has been shown within studies of tropical forests worldwide.Furthermore it is investigated whether there is evidence for the packing hypothesis of Midgley et al. [14] and whether the relationships between tree height and altitude exist within eastern Australian woodlands and forests.The implications of the results are discussed in terms of predicted anthropogenic climate change within the region.

Study Area.
The study was undertaken primarily within the New England Tablelands, Nandewar, and Brigalow Belt South Bioregions [22] within New South Wales (Figure 1; approximately 8 million hectares).In the region of investigation 2,341 sites were placed almost exclusively within the national reserve system (national parks, nature reserves, state conservation areas, community conservation areas, and aboriginal areas) and on privately managed reserves (Nature Conservation Trust Agreements and Voluntary Conservation Agreements).Such locations where chosen in order to minimise structural disturbance factors and different management objectives.Seventy-seven private and public conservation areas were surveyed and sites placed randomly within prechosen strata.These strata were generally based on rock type, altitude, aspect, and landform position allowing for replication within each conservation area.No locations were under grazing or any other primary production and the focus of land management was conservation (maintain or improve current natural condition).Sample sites were chosen for minimal postcolonisation anthropogenic disturbance.Any locations with obvious overstorey removal were excluded from the primary dataset.
Altitudes ranged from 122 m above sea level at Narran Lakes (30 ∘ 31  , 151 ∘ 30  ) to 1,371 m above sea level at Butterleaf northeast of Glen Innes (29 ∘ 30  , 152 ∘ 00  ).The study area is dominated by the Dividing Range which runs southwest along the eastern edge of the region with an associated high altitude tableland area (over 800 m a.s.l).Altitudes fall to the west levelling to open plains.Some isolated high altitude areas occur in western locations away from the main range (Mt Kaputar 1,300 m; Warrumbungle Range 1,100 m).The most northerly sites were within the Yelarbon area (28 ∘ 42  , 140 ∘ 47  ) and the most southerly nearby Biddon (31 ∘ 38  , 148 ∘ 48  ).Annual mean temperatures range from 11.7 ∘ C to 20.5 ∘ C and the annual mean rainfall from 422 mm to 1,140 mm.Around 60% of annual rainfall occurs from November to March; thus rainfall is seasonal with a marked summer peak and dry winters.Seasonality increases from south to north.Rainfall is proportionately higher along the far east of the Dividing Range as it receives additional orographic precipitation while the western areas are increasingly within rain shadow.Generally taller forests and closed forests occur along the eastern escarpment with open forests and woodlands dominating in intermediate altitudes and open woodlands and shrublands common within the extreme western and lower altitude locations.

Site Survey.
A random tree with a minimum diameter at breast height (130 cm) (DBH) of 10 cm was chosen as a central tree at each location.From this tree 15 of the closest trees of over 10 cm DBH (130 cm) were recorded (16 trees inclusive of the central stem).The size of these tree sites was unbounded and dependent on the distance to the furthest tree recorded.For each tree the identity of the species and distance to the first tree (meters) were recorded.At each location information was gathered on location and altitude with reference to Geographical Positioning System (GPS; Garman Map 62s; Garmin International, Olathe, KA, USA) along with physiography, soil drainage, and depth.Sites were surveyed over a period of six years (2006 to 2012).Nomenclature follows that given by PlantNET (http://plantnet.rbgsyd.nsw.gov.au/accessed December 2012).An average of all stem distances from the central stem was used for calculations of above ground biomass in addition to that provided by the furthest stem as a form of standardisation.To assess the effect of clumping, the evenness of the distribution of distances within each stem plot was calculated using PAST [23].Evenness essentially measures the dispersion of stems.If all stems have the same distance from the central stem then a perfectly even score would result (+1); however if stems are dispersed unevenly (most toward the centre a few further away or vice versa) within the stem plot the evenness score will be much reduced (towards 0).The most uneven plots were excluded from the analysis.Evenness was an attribute that was tested in its own right and in combination with the AGB scores.

Explanatory Variables.
The above ground biomass (AGB) of all species was calculated based on the dry forest equation provided by Chave et al. [18], where  is the oven-dry wood specific gravity (grams per cubic meter),  is the diameter (in centimeters) of the tree at breast height (130 cm), and  is tree height (in meters).
The resulting estimate of AGB is in kilograms [18].Wood density was derived from the average of available sources for each species ( [24]; http://www.timberanswers.com/accessed February 2013).Where only moist wood density was available the figure was converted to the oven-dry equivalent using the conversion supplied by Chave et al. [19] where  is the radius, Wood density is considered to be a taxonomically conservative trait [2,3,6,17,25]; thus for species where this trait was unavailable the average from the species closest relatives was assigned.AGB was calculated for each individual tree.The above ground biomass scores were then extrapolated to kilograms per hectare based on the size of the original plot and the combined scores of all individual trees.A figure for the combined stem m 2 per hectare was calculated based on the area of the stem plot raised to the size of a hectare.
International Journal of Ecology depression; scored 1 to 6), soil depth (deep to skeletal; scored 1 to 3), drainage (waterlogged to well drained; scored 1 to 4) original plot size (all of which were measured in the field), and stem plot evenness.Thus a total of 43 variables were calculated or scored for each site.

Statistical Analyses.
Species data was presented in the form of a matrix of species by sites.The first matrix was populated by the AGB of each individual.Three additional matrices were created where the AGB was based on the size of the original plot, the original plot scaled to one hectare based on the furthest stem from the centre, and the original plot scaled to one hectare using the average distance to all stems from the central tree."Environmental" data was presented as a matrix of factors by sites.The CANOCO5 package [27] was used to test the effects of physical, locational, and climatic variables on the unscaled AGB data and the two datasets scaled to one hectare.This was done in order to test the robustness of results and to quantify the effect of scaling to one hectare.As the response data is compositional and the gradient lengths were 9.4 SD units long, unimodal canonical correspondence analysis (CCA) was considered the most appropriate analysis technique by the advisor module within CANOCO5.Significant environmental variables were chosen by the manual forward selection procedure in CANOCO5, followed by Monte Carlo permutation tests (1000 iterations).Forward selection allows for the removal of collinear/redundant variables from the regression model and subsequent reanalysis of the contribution of the remaining unselected variables as each new variable is added to the model.Holm's correction was used in order to adjust the significance values for potential inflated family-wise type I errors, thus reducing the chances of colinearity and false significance [27].Only variables which achieved a  value of less than 0.05 after Holm's correction were included within the model.The effect of significant variables on AGB was tested using regression modelling.Generalised additive modelling (GAM) was used to model the single attribute of AGB against the constrained ordination space.GAM models of differing model complexity starting with the null model through differing degrees of freedom were generated.The GAM model with highest parsimony was tested via Akaike Information Criterion.Simple linear regression was also performed on AGB and specific gravity per plot to test conceptual models of the contribution of larger diameter trees and changes in wood density across the environmental gradient sampled.The evenness of plots was also regressed against AGB and specific gravity to test if there may have been significant effect of differences in dispersion on the results obtained.

Results
A total of 197 species were identified over 10 cm DBH within sites and 36,334 individuals were measured for diameter, height, and dispersion.The most abundant genera were Eucalyptus and Acacia which accounted for 67 and 27 species, respectively.1. Axes are the first two axes of the CCA analysis.
There was no significant effect of evenness on AGB or specific gravity, at each site across the dataset.There was a high degree of congruence between the CCA results for the scaled and unscaled AGB datasets.Within the scaled dataset original plot size was added as a variable and found to be nonsignificant.Evenness (test of clumping) also did not appear as a significant variable in any of the analyses.The nonsignificance of plot size and evenness are likely due to the size and extent of the dataset which significantly downweights the influence of differences in clumping of trees across sites.Thus it was considered that the scaled to one hectare analysis could be used without undue influence of scaling error.Thirty-six of the 43 environmental variables were found to be of significance after Holms correction ( ≤ 0.05) in relation to variation in AGB (Figure 2).These 36 significant variables accounted for 11.7% of the total inertia of the scaled dataset.The first two axes of the CCA model explained 22.7% of the variance within the dataset.MaTWaP and PWP (seasonal summer attributes) were of greatest explanatory power for AGB within the dataset (Table 1).The explained variation in AGB was directly explored using GAM (Figure 2).Overall the 12 variables with greatest explanatory power (as assessed by the arbitrary cut-off point of 0.3% explained variance; Table 1) indicate that an increase in moisture within the soil (MMILQ, MMIWaQ), greater precipitation (PWP, AMR), higher altitudes but lower radiation (AMR, HPR), and lower temperatures (AMT, MaTWaP, and MeTWeQ) are associated with an increase in AGB (Table 1; Figure 2).
Despite the incorporation of a significant number of climatic variables the spatial variables of latitude, longitude, and altitude were found to be of significance.These spatial variables are accounting for unmeasured environmental attributes.Two significant variables of minor explanatory power were site physical characteristics (physiography and drainage).There is a slight positive relationship between physiography and drainage variables indicating greater AGB at lower topographic positions but on well drained sites.Simple linear regression showed a relationship (14.5% of variance) between the AGB of the largest stem within a site and the AGB of the entire site (Figure 3).Allometric traits are known to be phylogenetically conservative [2]; therefore some additional graphical and regression techniques have been applied.A negative relationship (29% of variance) between the average SpecGrav (wood density) for all Eucalyptus and altitude exists (Figure 4).Linear regression was performed on average Eucalyptus height records for each site excluding all other genera to control for phylogeny (Figure 5).In addition the most abundantly recorded Myrtaceae taxa that occurred across the greatest altitudinal range were regressed against altitude (Figure 5).Overall there is an increase in the height of Eucalyptus species with increasing altitude.In contrast, there is a slight negative relationship between the heights of widespread individual species across the altitudinal gradient (Figure 5).
The proportional distribution of AGB across major taxonomic groups and altitude is shown in Figure 6.Partitioning the contribution of major taxonomic groups shows clearly that there is a change in the dominance of broad taxonomic groups over the altitudinal gradient.Softwood species (Cupressaceae: Callitris endlicheri and C. glaucophylla) provide approximately 30% of the total AGB between 400 and 800 m altitude.Fabaceae (primarily Acacia) is most important in terms of AGB at lowest altitudes (<800 m).
Eucalyptus, though prevalent across all altitudinal zones, accounts for 80% of all the AGB with Callitris accounting almost no AGB at the highest altitudes [28].The importance of Callitris also decreases at the lowest altitudes with an increase in both Eucalyptus and Casuarinaceae (Allocasuarina and Casuarina).The increase in the prominence of Casuarinaceae is due to the prominence of Allocasuarina luehmannii and Casuarina cristata at lower altitudes.Rutaceae are of minor prominence within the region at lower altitudes and are replaced by Proteaceae at higher altitudes.The results potentially indicate greater taxonomic diversity within the overstorey within the mid altitudes and a weakening of the biomass dominance of Eucalyptus within this zone.

Discussion
Latitude, longitude, and altitude were all found to be of significance in the CCA model explaining variance in select allometric traits (Figure 2).Despite 15 climatic and two physical variables being of significance in the CCA model it is obvious that there were still unmeasured traits that covary with these spatial variables.It is difficult to predict what these might be; in relation to altitude they could include a reduction in photosynthesis due to the decrease in atmospheric pressure [10].As the dividing range runs southwest the combination of latitude and longitude would be acting along this geographic feature.It is unknown what latitude and longitude may account for, though the angle at which solar radiation occurs would vary along a latitudinal gradient.Latitude has featured as a significant variable in angiosperm wood traits in a number of investigations in a wide variety of communities and locations [2,19,29].AGB was primarily correlated with an increase in moisture availability, particularly during the driest period and with an increase in the moisture index (Figures 2 and 3).Higher mean moisture index values indicate a surplus of water (precipitation minus evaporation) available for retention in the landscape.However annual mean rainfall was not found to be a significant factor in the CCA model potentially indicating that the amount of rainfall is less informative than when that rainfall occurs and how much is retained after evaporation.Although AGB was negatively correlated to the annual mean radiation (AMRad), the highest AMRad is associated with the most northerly and western locations which are the driest.AGB was largely unresponsive to changes in temperature and other measures of radiation.Although higher temperatures were uninformative, lower temperatures are correlated with some increase in AGB (Figure 3).These results are in agreement with a number of investigations that have shown moisture availability to be a better predictor than temperature and solar radiation in drier climates [3,11,[30][31][32].The results also concur with Keith et al. [13] and Culmsee et al. [9] who found biomass to peak in cool wet environments but do not match those of Banin et al. [3].The system investigated here is a temperate one rather than tropical and moisture may be retained during dry seasons due to an increase in the covariables associated with increasing moisture index.Overall AGB was highest where moisture was available all year round (+ve to PDQ, PDP, PWP, HPMI, and MMILQ) either through direct rainfall and/or water retention in the landscape.These results are more akin to the results of Slik et al. [11] who suggested that stable wet climates were important for accumulation of AGB.Current climate change scenarios for the region predict that, while mean annual rainfall may not diminish significantly seasonality of rainfall will increase along with a reduction in the mean moisture index [33].This will lead to a reduction in AGB and therefore reduce the ability of these forests and woodlands to store carbon.
It is logical to assume a relationship between wood density and the number of stems as both scores are used to calculate AGB (see Section 2), though Stegen et al. [7] have shown that this is not a consistent rule across spatial scales  or different forest types.Through similar reasoning it can be assumed that height and diameter should correlate strongly with AGB due to their importance in the calculations; this was not the case here.Therefore, at least in this region the number of stems and the density of wood (specific gravity) are of primary importance to overall AGB.Further analysis by linear regression of the biomass of the largest tree against the total site biomass showed a small but significant positive relationship (explaining 14.5% of the variance).As discussed by Midgley et al. [14] it is expected that a negative relationship would result if competition was limiting AGB and that a positive relationship signifies that competition is weak.Though only 14.5% of the variance in overall AGB is explained by the AGB of the largest stem there is some support for the species packing hypothesis and weak competition.
Stegen et al. [7] suggested that disturbance and the temporal trends associated with recovery can decouple the importance of tree size and tree density relationships with AGB.Niklas et al. [34] also argued that disturbance rates were important in determining maximum plant size and that disturbance, not competition, had a central role in determining community basal area and how biomass accumulates.It is possible that the lack of disturbances or the type of those disturbances and the ability of the species to cope with these disturbances have increased the link between AGB and the density of stems and wood density.The most prevalent landscape level disturbances within the study region are fire and drought.Fire and drought are disturbance events that cause disproportional mortality by removing both larger stems and the youngest individuals; both events are independent of resource availability and competition.Thus, independently fire and drought are likely to promote intermediate sized stems and, depending on frequency and severity, to decouple competitive interactions between size classes [35,36].Hence, fire and drought may diminish the importance height and International Journal of Ecology diameter in calculations of AGB and reduce the correlation between the largest stem and site AGB.Evidence suggests that interactions of fire and drought may cause synergistically greater disturbance within the study region exacerbating effects on species [37,38].It is predicted that within the region there will be an increase in the regularity and severity of drought and fire [33].Though not part of this study it would be worth investigating how wood density and AGB vary across communities that experience different fire frequencies.
Many allometric properties are phylogenetically conservative [2,3,6,17,25].Thus, it is necessary to view the taxonomic distribution of the values obtained in order to control for this confounding factor.Partitioning the contribution of major taxonomic groups shows clearly that there is a change in the dominance of broad taxonomic groups over the altitudinal gradient (Figure 6).Softwood species provide approximately a significant component of the total AGB between 400 and 800 m altitude largely at the expense of Eucalyptus ([28], Figure 6).The increasing prevalence of Eucalyptus and the decrease in the prominence of Callitris could account for much of the increase in AGB with altitude.Similarly Eucalyptus generally have greater stem widths and height than Callitris and Allocasuarina; this disparity may also account for the decline in height and stem widths with decreasing altitude.To remove this potential confounding factor Eucalyptus was isolated for analysis but still showed the same overall patterns with increasing altitude.The average wood density (specific gravity) of Eucalyptus stems decreases with increasing altitude (Figure 4) indicating that these correlations are not simply an artefact of increased prevalence of softwood species at lower altitudes and therefore not due to phylogeny.This result suggests that lower altitudes sites (lower rainfall and moisture index and increased temperatures and radiation) favour Eucalyptus species with higher wood density (specific gravity).Higher altitudes sites (increased rainfall and mean moisture index, lower temperatures, and radiation) favour species with lower wood densities (specific gravity).
There is an increase in the height of all Eucalyptus stems (interspecies comparison) with increasing altitude (Figure 5).This was not found during intraspecies comparisons of the most widespread taxa (Figure 5).Within individual species that occurred across a wide range of altitudes there was a very slight decrease in height with increasing altitude across all examined in contrast to the overall pattern of increased average height within stands.Within Eucalyptus, species with lower wood density and taller growth habits are favoured at higher altitudes (wetter and cooler) while shorter species with higher wood density are favoured at lower altitudes (drier and hotter).This supports the findings of Slik et al. [11] who showed that drought disproportionately increased the mortality of large trees with low wood densities and also supports the hydraulic limitation hypothesis [3,4].The results are also in agreement with those presented by Lines et al. [4] who showed that landscape level trends in observed tree heights and other allometric properties were primarily a result of changing species composition along gradients rather than changes within species.Species show comparatively little adaptability in their allometric characteristics with changes in climate conditions [4,21].Species have a poor capacity to adapt to predicted climate change and a change in composition is the most likely outcome in the short term.

Conclusion
AGB increases in wetter and cooler locations.This increase is related to the amelioration of seasonal rainfall deficits that occur in this summer dominant rainfall zone.Taller and greater diameter trees with lower wood densities occur at higher altitudes.These changes in AGB, though not linked to phylogeny, are based on a change in the composition of species rather than a change in the allometric properties within a species.The data supports the hydraulic limitation hypothesis and provides some evidence for species packing.AGB at a site is associated closely with an increase in the density of stems at a site and the average wood density of trees and less so by the influence of the largest tree and average stem height and diameter.Within this region these interrelationships are likely affected by interactions with fire frequency.Fire is a common feature of much of the study area; fire frequency and related variables require further investigation in this context.Both drought and fire are likely to become more frequent and severe with climate change.Drought will increasingly become the most limiting factor within the region and thus, singularly and in combination with fire, will drive changes in allometric properties and subsequently species interactions.Based on the evidence provided in this investigation significant changes in species composition and the distribution of allometric traits will occur across the landscape with changing climate.These changes will include a reduction in above ground biomass, the number of stems per hectare, the average height, average diameter, and basal area.The largest of trees are likely to be removed early due to their inability to cope with increased drought stress.A flow on of these affects will be a marked reduction in carbon storage across north eastern Australia.

Figure 2 :
Figure 2: Generalised additive model based contour plot for AGB per plot against the 12 most significant environmental variables (minimum of 0.3% explained variance) within the CCA model.Codes given in Table1.Axes are the first two axes of the CCA analysis.

6 InternationalFigure 3 :Figure 4 :
Figure 3: Contribution of the above ground biomass (AGB) of the largest tree per plot and the total AGB of the plot.

Figure 5 :
Figure5: Change in average tree height of all Eucalyptus species within site (interspecific variation) and the average height of five selected widespread species (intraspecific variation) against the altitude gradient.Note there is a small negative change in the average tree heights of individual species with increasing altitude.Overall there is an increase in the height of the summed Eucalyptus scores for each site.

Figure 6 :
Figure 6: Variation in the contribution of major taxonomic groups to above ground biomass (AGB) across an altitudinal gradient in eastern Australia.

Table 1 :
Codes and descriptions of significant variables ( ≤ 5%) listed in order of decreasing importance (contribution) from forward selection in canonical correspondence analysis.